Thanks so much Jia<br>But Sir the atomic positions in the input is for 3 atoms, should l leave the atomic position and change the number of atom or my atomic position is wrong??<br>thanks<br>Musari A. A<br>UNAAB Nigeria<br>
<br><div class="gmail_quote">On Mon, Aug 6, 2012 at 3:34 PM, jia chen <span dir="ltr"><<a href="mailto:jiachen@princeton.edu" target="_blank">jiachen@princeton.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Dear Abolore Musari,<div><br></div><div>I guess the number of atoms should be 14, instead of 3.</div><div><br></div><div>Regards</div><div>Jia<br><br><div class="gmail_quote"><div><div class="h5">On Mon, Aug 6, 2012 at 10:24 AM, Abolore Musari <span dir="ltr"><<a href="mailto:abmus007@gmail.com" target="_blank">abmus007@gmail.com</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div class="h5">Dear QE users<br>Sir, I am try to find the lattice constant of magnetite (Fe3O4) which is a spinel ferrimagnetic material. The expt lattice constant was 15.87 and I have been trying to find the equilibrium lattice constant (celldm(1)) between 15.80 - 16.00 and all i get is a straight line graph instead of the usual exponential curve graph. I have tried all l could most esp with starting_magnetization but the results are the same. My input is pasted below pls help me to get the appropriate graph.<br>
The kind of the graph I got is attached.<br>Thanks<br><br>&control<br> calculation = 'scf',<br> restart_mode ='from_scratch',<br> prefix = 'Fe3O4',<br> pseudo_dir = '/home/abolore/Programs/Pseudos/',<br>
outdir = './tmp/'<br> /<br> &system<br> ibrav = 2, <br> celldm(1)= $15.80-16.00,<br> nat = 3, <br> ntyp = 3,<br> ecutwfc = 40,<br> ecutrho = 400,<br> starting_magnetization(1) = 0.7,<br>
starting_magnetization(2) = -0.5,<br> starting_magnetization(3) = 0.0,<br> nspin = 2,<br> occupations='smearing',<br> smearing='gaussian',<br> degauss = 0.05,<br> /<br> &electrons<br>
diagonalization = 'david',<br> mixing_mode = 'plain',<br> mixing_beta = 0.7<br> /<br>ATOMIC_SPECIES<br> Fe1 55.845 Fe.pz-nd-rrkjus.UPF<br> Fe2 55.845 Fe.pz-nd-rrkjus.UPF<br> O 16.000 O.pz-rrkjus.UPF<br>
<br>ATOMIC_POSITIONS<br>Fe1 0.125 0.125 0.125 <br>Fe2 0.500 0.500 0.500<br>O 0.2548 0.2548 0.2548<br> <br>K_POINTS (automatic)<br> 8 8 8 1 1 1<br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br><br>
<br></div></div>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br>Jia Chen<br><br><br>
</font></span></div>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br>