<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"><div>Dear all:</div><div>I'm working on a perovskite system to practice doing PDOS and band structure calculations. In the past, I've been able to calculate PDOS for 2D systems with success. But in this case I have not being able to calculate correctly the PDOS. After doing a scf calculation and then the projwfc.x calculation, I obtain the following for the PDOS output file (only the first lines copied):</div><div># E (eV)  dos(E)    pdos(E)<br>*******  0.473E+01  0.447E+01<br>*******  0.373E+01  0.352E+01<br>*******  0.290E+01  0.274E+01<br>*******  0.224E+01  0.212E+01<br>*******  0.171E+01  0.162E+01<br>...</div><div>As you can see, the problem is that the first column is not shown. I have experimented with different variations of the input file. This is what I've tried:</div><div><ol><li>I have tried a scf calculation and then the projwfc calculation</li><li>I've tried scf, nscf and then the projwfc one.</li><li>I've tried relaxing the structure (originally obtained from the COD database <a href="http://www.crystallography.net/cod/">http://www.crystallography.net/cod/</a> and working with the relaxed structure.</li><li>I've tried using the original structure as obtained from the COD database.</li><li>I've tried different values of ecut and the k-point grid.</li><li>I've tried changing the DeltaE parameter because I thought it was a matter of the number if significant figures allowed in that column.<br></li></ol></div><div>The Fermi Energy that I obtain in the scf calculation is always around -183 eV<br></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_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">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></div></div></div></div></div></div></div></div>