Dear Biswaranjan!<div>                             First try to get the experimental (or theoretical value done by somebody else)  value of the lattice parameter of your system. If you are using ultrasoft pseudopotential then ecutrho between 8-10 times of ecut will suffice. For fixing 'ecutwfc' you have to do "vc-relax"  calculation (see doc-input-pw for more detail) at different 'ecutwfc' then see the convergence with respect to ecut. When the difference between energies for two consecutive ecutwfc is below the minimum limit of accuracy you want then choose lower value(difference of 5 Ry. between two consecutive ecut will suffice). The same procedure you have to repeat for convergence with respect to k-points. You can write a script for doing all the calculations in one shot.<br>
<br><div class="gmail_quote">On Wed, Sep 8, 2010 at 1:24 PM,  <span dir="ltr"><<a href="mailto:bisu123@platinum.materials.iisc.ernet.in">bisu123@platinum.materials.iisc.ernet.in</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Dear all,<br>
             Can someone tell me the correct methodology to start with<br>
lattice  parameters estimation for HCP metals. I shall tell<br>
what I did for my earlier simulations involving Magnesium.<br>
<br>
First, I fixed the ecut value and then used the same ecut value to fix the<br>
ecutrho value. Then using the fixed ecut and ecutrho, I varied the<br>
k-points to obtain the K-points. in my case I obtained [12 12  9]. Then<br>
for lattice  parameter calculation, [ having fixed ecut, ecutrho, k-points<br>
of course]  I fixed  'c/a' ratio initially[ c/a= sqrt(8/3)] and varied the<br>
'a' across the ideal value of 'a'[ a=3.21 Angstroms =6.06602 at. units]<br>
and obtained the value for 'a' corresponding for the minimum energy<br>
configuration. Then using the obtained 'a' in the previous step, I varied<br>
the 'c/a' ratio across the ideal value to obtain c/a=1.58, corresponding<br>
to the minimum energy configuration. However this estimate involves<br>
significant error since the ideal value for c/a( for Mg=sqrt(8/3)=1.624).<br>
Then I used the value for 'c/a' in the previous step[ i.e 1.58] and varied<br>
'a' to find the value of 'a'.<br>
<br>
Finally, I obtained a=6.1 a.u=3.23 angstroms, c/a=1.58. The<br>
psuedopotential I used for Mg is  Mg.pw91-np-van.UPF .<br>
<br>
Thanks all,<br>
<br>
regards,<br>
Biswaranjan Dash<br>
Research Scholar,<br>
Deformation Mechanisms Modeling Group,<br>
C/0- Asst Prof S.Karthikeyan<br>
Indian Institute of Science,Bangalore<br>
India<br>
<br>
<br>
<br>
<br>
<br>
--<br>
This message has been scanned for viruses and<br>
dangerous content by MailScanner, and is<br>
believed to be clean.<br>
<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>
</blockquote></div><br><br clear="all"><br>-- <br>Regards,<br>MOHNISH,<br>-----------------------------------------------------------------<br>Mohnish Pandey<br>Y6927262,5th Year dual degree student,<br>Department of Chemical Engineering,<br>
IIT KANPUR, UP, INDIA<br>+919235721300<br>-----------------------------------------------------------------<br>
</div>