<div dir="ltr">Dear ELIiot,<br>Here you are using the pz lda method which always underestimate the lattice parameter. there is no error in your calculation.<br><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Tue, Dec 3, 2013 at 5:49 PM, Elliot Menkah <span dir="ltr"><<a href="mailto:elliotsmenkah@yahoo.com" target="_blank">elliotsmenkah@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi all,<br>
<br>
I'm determining the lattice parameter of Nickel and my convergence test<br>
is giving me 3.42 Angstroms against and experimentally determined 3.52.<br>
Previous calculations by other people are showing 3.49 and 3.52<br>
Angstroms as again.<br>
<br>
What could be wrong with my parameters or calculations.<br>
<br>
Below is the input file please.<br>
<br>
<br>
#!/bin/sh<br>
# check whether echo needs the -e option<br>
#<br>
executablepath=/home/molecular/espresso-5.0.2/bin/pw.x<br>
for alat in  6.273 6.369 6.463 6.557 6.652 6.746 6.841 6.935<br>
do<br>
#input file<br>
cat > Ni-$<a href="http://alat.in" target="_blank">alat.in</a> << EOF<br>
&CONTROL<br>
        title = 'Ni-ferro-mag' ,<br>
        calculation = 'relax' ,<br>
        restart_mode = 'from_scratch' ,<br>
        outdir = './temp' ,<br>
        pseudo_dir = '/home/molecular/Elliot/pseudo' ,<br>
        prefix = 'Ni-ferromag',<br>
        tstress = .true. ,<br>
        tprnfor = .true. ,<br>
/<br>
  &SYSTEM<br>
        ibrav = 2,<br>
        celldm(1) = $alat,<br>
        nat = 1,<br>
        ntyp = 1,<br>
        ecutwfc = 60,<br>
        ecutrho = 480,<br>
        occupations='smearing',<br>
        smearing='fermi-dirac',<br>
        degauss=0.003,<br>
        nbnd=10<br>
        nspin=2,<br>
        starting_magnetization(1)=0.5,<br>
/<br>
  &ELECTRONS<br>
         conv_thr = 1.0D-8<br>
         mixing_beta = 0.2,<br>
/<br>
  &IONS<br>
/<br>
ATOMIC_SPECIES<br>
         Ni1 58.6934   Ni.pz-nd-rrkjus.UPF<br>
ATOMIC_POSITIONS (alat)<br>
        Ni1       0.000000000   0.000000000   0.000000000<br>
K_POINTS automatic<br>
   7 7 7   0 0 0<br>
<br>
EOF<br>
$executablepath < Ni-$<a href="http://alat.in" target="_blank">alat.in</a> > Ni-$alat.out<br>
grep "Final energy" Ni-$alat.out >> Ell0.txt<br>
echo "$alat" >> Lat.txt<br>
done<br>
awk '{print $4}' Ell0.txt > Nel.txt<br>
paste Lat.txt Nel.txt > Final_alat.txt<br>
<br>
rm Lat.txt Ell0.txt<br>
<br>
<br>
Thank uou.<br>
<br>
<br>
<br>
Kind regards,<br>
<br>
ELliot.<br>
<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">Thanks and Regards<br>Bramha Prasad Pandey<br>GLA University<br>Mathura (U.P)<br>INDIA.<br></div>
</div>