<div dir="ltr"><div dir="ltr"><div dir="ltr">On Wed, Feb 20, 2019 at 4:35 PM Juan Manuel Radear <<a href="mailto:radear82@gmail.com">radear82@gmail.com</a>> wrote:</div><div dir="ltr"><br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"># E (eV) dos(E) pdos(E)<br><div>******* 0.473E+01 0.447E+01</div></div></div></div></div></div></div></div></div></div></blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"> [...] As you can see, the problem is that the first column is not shown<div>[...] The Fermi Energy that I obtain in the scf calculation is always around -183 eV<br></div></div></blockquote><div><br></div><div>the format used for the energy is f7.3, so energies around -183 eV will simply not fit into it (7 figures - 3 decimals - 1 dot - minus sign = 2 figures left). Just increase it to f8.4:</div><div><br></div><div>$ grep f7.3 PP/src/partialdos.f90<br> WRITE (4,'(f7.3,2e11.3,14e11.3)') etev*rytoev, &<br> WRITE (4,'(f7.3,4e11.3)') etev*rytoev, (dostot(ie,is,ik), is=1,nspin), &<br> WRITE (4,'(f7.3,2e11.3,14e11.3)') etev*rytoev, ldos(ie,1,ik), &<br> WRITE (4,'(f7.3,2e11.3,14e11.3)') etev*rytoev, ldos(ie,1,ik), &<br> WRITE (4,'(f7.3,2e11.3,14e11.3)') etev*rytoev, &<br> WRITE (4,'(f7.3,4e11.3)') etev*rytoev, dostot(ie,ik), &<br><br></div><div>Paolo<br></div><div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div></div><div>If anyone could give some insight on what could be done to correct this error, I would appreciate it immensely. Here is the scf file of my most recent try:</div><div>--------------------------------------------------<br></div><div>&CONTROL<br> calculation='scf',<br> outdir='./tmp',<br> prefix='CsGeI3-v7',<br> pseudo_dir='./',<br> verbosity='high',<br> nstep = 2800,<br>/<br><br>&SYSTEM<br> ibrav=4,<br> celldm(1)=15.7947296813d0, celldm(3)=1.2693882698d0,<br> nat=15,<br> ntyp=3,<br> ecutwfc=30,<br> occupations='smearing',<br> smearing='mv',<br> degauss=0.005d0,<br>/<br><br>&ELECTRONS<br> diagonalization='david',<br> mixing_mode = 'local-TF',<br> mixing_beta = 0.1, <br> conv_thr = 1.0d-5,<br> electron_maxstep= 2000,<br>/<br><br>ATOMIC_SPECIES<br> Cs 132.905000d0 Cs.pbe-spnl-kjpaw_psl.1.0.0.UPF<br> Ge 72.640000d0 Ge.pbe-dn-kjpaw_psl.0.2.2.UPF<br> I 126.904000d0 I.pbe-n-kjpaw_psl.0.2.UPF<br> <br>ATOMIC_POSITIONS {crystal}<br> I 0.5027350000d0 0.4972750000d0 0.0042000000d0<br> I 0.5027250000d0 0.0054600000d0 0.0042000000d0<br> I 0.9945400000d0 0.4972650000d0 0.0042000000d0<br> I 0.1694016700d0 0.8306083300d0 0.3375333300d0<br> I 0.1693916700d0 0.3387933300d0 0.3375333300d0<br> I 0.6612066700d0 0.8305983300d0 0.3375333300d0<br> I 0.8360683300d0 0.1639416700d0 0.6708666700d0<br> I 0.8360583300d0 0.6721266700d0 0.6708666700d0<br> I 0.3278733300d0 0.1639316700d0 0.6708666700d0<br> Cs 0.0000000000d0 0.0000000000d0 0.0100010000d0<br> Cs 0.6666666700d0 0.3333333300d0 0.3433343300d0<br> Cs 0.3333333300d0 0.6666666700d0 0.6766676700d0<br> Ge 0.0000000000d0 0.0000000000d0 0.5394317000d0<br> Ge 0.6666666700d0 0.3333333300d0 0.8727650300d0<br> Ge 0.3333333300d0 0.6666666700d0 0.2060983700d0<br><br>K_POINTS {automatic}<br> 4 4 4 0 0 0</div><div>---------------------------------------------------------<br></div><div>And here is the PDOS input file:</div><div>--------------------------------------------------------<br></div><div>&PROJWFC<br> outdir='./tmp',<br> prefix='CsGeI3-v7',<br> Emin = -198,<br> Emax = -168,<br> degauss = 0.01,<br> DeltaE=0.01,<br> filpdos = 'CsGeI3-v7.pdos',<br>/</div><div>-----------------------------------------<br></div><div>As I mentioned before, the Fermi E is around -183, and I chose the Emin and Emax around that value, as I have seen in different tutorials. I have also tried Delta E= 0.1 to try and reduce the number of significant figures shown in the output file but with no success. <br></div><div>Thanks again in advance for your attentions and sorry for the inconveniences.<br></div><div><br></div><div><div dir="ltr" class="gmail-m_5549791023923931193gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Dr. Juan Manuel Ramirez de Arellano</div><div>Profesor de Cátedra, Departamento de Ciencias, <br></div><div>Tecnologico de Monterrey, Calle del Puente 222, 14380, Tlalpan,<br></div><div>Mexico</div><div>Email: <a href="mailto:jramirezdearellano@tec.mx" target="_blank">jramirezdearellano@tec.mx</a><br></div><div><font color="#0000ee"><u><br></u></font></div></div></div></div></div></div></div></div></div></div></div></blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
_______________________________________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div></div></div>