<div dir="ltr"><div><div><div><div><div><div>Hi,<br><br></div>I am calculating the optimum lattice constant nd cohesive energy of Zr hcp crystal. The total energy of a sinle Zr atom I got from input file Zr.in is -5.53133015 Ry, the total energy of Zr hcp crystal I got from input file <a href="http://Zr.hcp.in">Zr.hcp.in</a> is -12.13562915 Ry, so the cohesive energy is -0.536 Ry, converted to -7.30 eV. From a reference paper, the cohesive energy should be -6.2eV. Attached below are Zr.in and <a href="http://Zr.hcp.in">Zr.hcp.in</a>. Could anyone tell me if I had made any mistake in those input files? Thanks.<br>

<br><font size="4">Zr.in:</font><br><br>&control<br>    prefix='Zr',<br>    pseudo_dir = '/home/zhang/Software/espresso_PP/',<br>    outdir='/share/scratch/tmpZr/eCut120SM0.01/',<br> /<br> &system    <br>

    ibrav=  1, celldm(1) =50, nat=  1, ntyp= 1,<br>    occupations='smearing', smearing='marzari-vanderbilt', degauss=0.01,<br>    ecutwfc =120, <br><br> /<br> &electrons<br> /<br>ATOMIC_SPECIES<br> Zr  91.224 Zr.pbe-mt_fhi.UPF<br>

ATOMIC_POSITIONS  <br>  Zr  0.000000000000 0.000000000000 0.000000000000<br>K_POINTS gamma<br><br><br><font size="4"><a href="http://Zr.hcp.in">Zr.hcp.in</a></font><br><br>&control<br>    prefix='Zr',<br>    pseudo_dir = '/home/zhang/Software/espresso_PP/',<br>

    outdir='/share/scratch/tmpZr/eCut120K18hcpSM0.01/',<br>    calculation='vc-relax'<br> /<br> &system    <br>    ibrav=  4, celldm(1) =6, celldm(3) =1.633, nat=  2, ntyp= 1,<br>    occupations='smearing', smearing='marzari-vanderbilt', degauss=0.01,<br>

    ecutwfc =120, <br> /<br> &electrons<br> /<br> &ions<br> /<br> &cell<br> /<br>ATOMIC_SPECIES<br> Zr  91.224 Zr.pbe-mt_fhi.UPF<br>ATOMIC_POSITIONS crystal <br>  Zr  0.000000000000 0.000000000000 0.000000000000<br>

  Zr  0.666666666666 0.333333333333 0.50000000000<br>K_POINTS automatic<br>18 18 36 1 1 1<br><br></div>PS The optimum lattice constant got from <a href="http://Zr.hcp.in">Zr.hcp.in</a> is very close to the data in literature. <br>

<br><br></div>I have tried different value of the smearing parameter ranged from 0.001 to 0.05, none of their cohesive energy is good. <br><br></div>If the smearing value is 0.1, simulation by <a href="http://zr.in">zr.in</a> (single atom) will stop at iteration # 10 without any error message. Why does it happen?<br>

<br></div>Thank you for your help,<br><br></div>Dongsheng<br></div>