<html><head><meta charset="ISO-8859-1"><style>body { line-height: 1.5; }blockquote { margin-top: 0px; margin-bottom: 0px; margin-left: 0.5em; }body { font-size: 10.5pt; font-family: 'Segoe UI'; color: rgb(0, 0, 0); line-height: 1.5; }</style></head><body>
<div><span></span><div style="font-family: 'Segoe UI', Tahoma; line-height: normal;">Dear Prof. Gianozzi,</div><div style="font-family: 'Segoe UI', Tahoma; line-height: normal;"><br></div><div style="font-family: 'Segoe UI', Tahoma; line-height: normal;">Thanks for your reply. While you say the occupation numbers are not actually used, do you mean that they have no effect to the integrated electron phonon coupling coefficient? Or you mean the contribution from the k points that show "NAN" occupation numbers are not counted for the integration over reciprocal space? </div><div style="font-family: 'Segoe UI', Tahoma; line-height: normal;"><br></div><div style="font-family: 'Segoe UI', Tahoma; line-height: normal;">Our Tc for MgB2 is much lower than experiment and previous QE results. We wonder whether the NAN error for occupation number is part of the reason. </div><div style="font-family: 'Segoe UI', Tahoma; line-height: normal;"><br></div><div style="font-family: 'Segoe UI', Tahoma; line-height: normal;">Thank you again for the help. </div><div style="font-family: 'Segoe UI', Tahoma; line-height: normal;"><br></div><div style="font-family: 'Segoe UI', Tahoma; line-height: normal;">Wei Liu</div></div>
<div><br></div><hr style="width: 210px; height: 1px;" color="#b5c4df" size="1" align="left">
<div><span><div style="MARGIN: 10px; FONT-FAMILY: verdana; FONT-SIZE: 10pt"><div>flying_lw@yeah.net</div></div></span></div>
<blockquote style="margin-top: 0px; margin-bottom: 0px; margin-left: 0.5em;"><div> </div><div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"><div style="PADDING-RIGHT: 8px; PADDING-LEFT: 8px; FONT-SIZE: 12px;FONT-FAMILY:tahoma;COLOR:#000000; BACKGROUND: #efefef; PADDING-BOTTOM: 8px; PADDING-TOP: 8px"><div><b>From:</b> <a href="mailto:paolo.giannozzi@uniud.it">Paolo Giannozzi</a></div><div><b>Date:</b> 2014-06-25 03:20</div><div><b>To:</b> <a href="mailto:pw_forum@pwscf.org">PWSCF Forum</a></div><div><b>Subject:</b> Re: [Pw_forum] Problem about occupation number in elph calculation</div></div></div><div><div>I don't think that those occupation numbers are actually used.</div>
<div>If they are, you get NaN everywhere</div>
<div> </div>
<div>P.</div>
<div> </div>
<div>On Tue, 2014-06-24 at 10:10 -0800, flying_lw@yeah.net wrote:</div>
<div>> Dear Quantum Espresso users,</div>
<div>> </div>
<div>> </div>
<div>> I am using QE v 5.0.2 to calculate the Tc of MgB2. The job is done</div>
<div>> normally, but I got one problem when I cheked the output</div>
<div>> file mgb2.elph.out. The occupation numbers of some k points are 'NaN'.</div>
<div>> I did the same job using QE v 5.1, all the 'NaN' become '0.0000'. I'm</div>
<div>> not sure what is happening. Will this cause something wrong in the</div>
<div>> final results?</div>
<div>> </div>
<div>> </div>
<div>> Part of the output file is as below:</div>
<div>> </div>
<div>> </div>
<div>> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div>
<div>> End of band structure calculation </div>
<div>> </div>
<div>> k = 0.0000 0.0000 0.0000 band energies (ev): </div>
<div>> </div>
<div>> -33.3325 -33.3325 -33.3154 -3.7339 5.7105 9.0432 9.0432 10.2941 </div>
<div>> 15.0631 15.0631 16.8819 </div>
<div>> </div>
<div>> occupation numbers </div>
<div>> 1.0000 1.0000 1.0000 1.0000 1.0000 0.0368 0.0368 -0.0001 </div>
<div>> -0.0000 -0.0000 -0.0000 </div>
<div>> </div>
<div>> k = 0.0000 0.0000 0.1089 band energies (ev): </div>
<div>> </div>
<div>> -33.3324 -33.3324 -33.3165 -3.5997 4.7550 9.0927 9.0927 11.5641 </div>
<div>> 15.1063 15.1063 16.8566 </div>
<div>> </div>
<div>> occupation numbers </div>
<div>> NaN NaN NaN NaN NaN NaN NaN NaN </div>
<div>> NaN NaN NaN </div>
<div>> </div>
<div>> k = 0.0000 0.0000 0.0363 band energies (ev): </div>
<div>> </div>
<div>> -33.3325 -33.3325 -33.3155 -3.7188 5.5783 9.0489 9.0489 10.4614 </div>
<div>> 15.0681 15.0681 16.8790 </div>
<div>> </div>
<div>> occupation numbers </div>
<div>> 1.0000 1.0000 1.0000 1.0000 1.0000 0.0333 0.0333 -0.0000 </div>
<div>> -0.0000 -0.0000 -0.0000 </div>
<div>> </div>
<div>> k = 0.0000 0.0000 0.1451 band energies (ev): </div>
<div>> </div>
<div>> -33.3324 -33.3324 -33.3172 -3.4981 4.2250 9.1283 9.1283 12.3377 </div>
<div>> 15.1375 15.1375 16.8387 </div>
<div>> </div>
<div>> occupation numbers </div>
<div>> NaN NaN NaN NaN NaN NaN NaN NaN </div>
<div>> NaN NaN NaN </div>
<div>> </div>
<div>> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div>
<div>> </div>
<div>> </div>
<div>> </div>
<div>> </div>
<div>> I have attached the input and output files. Thank you for your help.</div>
<div>> </div>
<div>> </div>
<div>> Best,</div>
<div>> Wei</div>
<div>> </div>
<div>> </div>
<div>> </div>
<div>> </div>
<div>> </div>
<div>> _______________________________________________</div>
<div>> Pw_forum mailing list</div>
<div>> Pw_forum@pwscf.org</div>
<div>> http://pwscf.org/mailman/listinfo/pw_forum</div>
<div> </div>
<div>-- </div>
<div> Paolo Giannozzi, Dept. Chemistry&Physics&Environment, </div>
<div> Univ. Udine, via delle Scienze 208, 33100 Udine, Italy</div>
<div> Phone +39-0432-558216, fax +39-0432-558222 </div>
<div> </div>
<div>_______________________________________________</div>
<div>Pw_forum mailing list</div>
<div>Pw_forum@pwscf.org</div>
<div>http://pwscf.org/mailman/listinfo/pw_forum</div>
</div></blockquote>
</body></html>