<DIV>Dear pwscf users,<BR>&nbsp;&nbsp; I have a question, quite naive. I tried to search the pdos of oxygen molecular, but I have to rely on your advice. as well as I know,the calculation procedes as follows:<BR>1) make a self-consistent calculation (input=O2.scf.in)<BR>2) run a non-self-consistent calculation (input=O2.nscf.in)<BR>3) run projwfc.x program (input=O2.pdos.in, output=O2.pdos.out)in order to produce a postscript file: O2.pdos_tot.<BR>&nbsp;&nbsp;&nbsp; After scf and nscf calculation, the program pdos was used as follows:<BR>&nbsp;&amp;INPUTPP<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = 'O2' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = '/home/olh/tmp/' ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; filpdos = 'O2.pdos',<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ngauss = 1,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; degauss = 0.02,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; DeltaE = 0.01 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Emin = -35 ,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Emax = -4 ,<BR>&nbsp;/<BR>&nbsp;&nbsp;&nbsp; However, when I calculated pdos, in order to produce a postscript file: O2.pdos_tot, I encounter a problem, output is displayed as follows:<BR>.....<BR>k =&nbsp; -0.3750000000 -0.2165063509&nbsp; 0.1678415576<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp; -30.88771 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp; -18.42806 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp; -12.12292 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp; -10.97042 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp; -10.96565 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp;&nbsp; -4.27073 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN</DIV>
<DIV>&nbsp;k =&nbsp;&nbsp; 0.0000000000 -0.4330127019 -0.1678415576<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp; -30.88770 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp; -18.42864 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp; -12.12059 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp; -10.96946 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp; -10.96552 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; e =&nbsp;&nbsp;&nbsp; -4.27538 eV<BR>&nbsp;&nbsp;&nbsp;&nbsp; psi = <BR>&nbsp;&nbsp;&nbsp; |psi|^2 =&nbsp;&nbsp; NaN</DIV>
<DIV>Lowdin Charges: </DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; Atom #&nbsp;&nbsp; 1: total charge =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN, s, p, d, f =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN&nbsp;&nbsp;&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin up&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN, s, p, d, f =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN&nbsp;&nbsp;&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin down&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN, s, p, d, f =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN&nbsp;&nbsp;&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polarization =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN, s, p, d, f =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN&nbsp;&nbsp;&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; Atom #&nbsp;&nbsp; 2: total charge =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN, s, p, d, f =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN&nbsp;&nbsp;&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin up&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN, s, p, d, f =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN&nbsp;&nbsp;&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin down&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN, s, p, d, f =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN&nbsp;&nbsp;&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polarization =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN, s, p, d, f =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN&nbsp;&nbsp;&nbsp;&nbsp; NaN<BR>&nbsp;&nbsp;&nbsp;&nbsp; Spilling Parameter:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NaN<BR>It displayed "NaN", not numerical, I am confused.<BR>Any suggestions is appreciated.</DIV>
<DIV>&nbsp;</DIV>
<DIV>Thank you in advance<BR>Lihui Ou<BR></DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;</DIV>
<DIV><BR>--<BR>
<DIV>
<DIV><FONT face="Arial" size="3">======================================</FONT></DIV>
<DIV><FONT face="Arial" size="3">Lihui&nbsp;Ou</FONT></DIV>
<DIV><FONT face="Arial" size="3">PH.D Candidate&nbsp;in Electrochemistry&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </FONT></DIV>
<DIV><FONT face="Arial" size="3">College of Chemistry and Molecular Science&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </FONT></DIV>
<DIV><FONT face="Arial" size="3">Wuhan University,<EM>430072</EM>,Hubei Province,China&nbsp;</FONT></DIV>
<DIV><EM><FONT face="Arial" size="3">E-mail:</FONT><A href="mailto:oulihui666@126.com" target="_blank"><FONT face="Arial" size="3">oulihui666@126.com</FONT></A></EM></DIV>
<DIV><FONT face="Arial" size="3">======================================</FONT></DIV></DIV></DIV>