<div dir="ltr">You have to do the spin-polarized calculation for isolated atom. Zr has two unpaired electron in outer d-shell.<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Jun 4, 2013 at 1:02 PM, Dongsheng Zhang <span dir="ltr"><<a href="mailto:zhdsheng21@gmail.com" target="_blank">zhdsheng21@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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" target="_blank">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" target="_blank">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" target="_blank">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" target="_blank">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" target="_blank">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>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr">Regards,<br>MOHNISH,<br>-----------------------------------------------------------------<br>
Mohnish Pandey,<br>PhD Student,<br>Center for Atomic-scale Materials Design,<br>Department of Physics,<br>Technical University of Denmark<br><div>-----------------------------------------------------------------</div></div>

</div>