<div dir="ltr">Thank You all for your kind help and suggestions. <div><br></div><div>I would further like to clarify a few more things from the suggestions received. </div><div><br></div><div><br></div><div>Dear Dr. Michael Mehl,</div><div><br></div><div>Quoting you from your reply, "<span style="font-size:12.8px">You must have minimum energy structure for your crystal structure at that volume.  This means minimizing the energy with respect to all atomic positions and all lattice parameters: b/a, c/a, etc., while keeping the volume fixed." </span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">I have to keep the volume of my structure fixed during the structural optimization? Is it the equilibrium volume obtained from the equation of fitting, that I have to keep fixed, while varying the other structural parameters? As you said the equilibrium lattice parameter comes very close to the experimental value given in the tutorial, but the bulk modulus value deviates beyond the 10% range in case of all the equations of state, specially for the Keane equation it deviates hugely while for the 2nd order Birch equation the deviation is the least. From the output data obtained after the fittings, the 'E_diff' values (that is the difference between the calculated total energy data and the fitted data at a given lattice parameter) are of the order of 10^-4 Ry and in case of the 2nd order Birch equation all the E_diff values are nil. </span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Dear Dr. Alberto,</span></div><div><span style="font-size:12.8px">Thank you for the link to the article. I will follow your suggestions and try again. </span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Dear Dr. Lorenzo,</span></div><div><span style="font-size:12.8px">Yes sir, the structure obtained is satisfactory and the equilibrium lattice parameter is also very close to the experimental value. </span></div><div><span style="font-size:12.8px">I have already tested the convergence of the energy with respect to the ecutwfc value, as was instructed in the tutorial. So I should keep on increasing the parameter and check every time if the experimental value is getting reproduced and stick to that ecutwfc value? Please correct me if I am wrong. I will also try the 2nd method suggested by you. </span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Thank you again to all of you. </span></div><div><span style="font-size:12.8px">Regards, </span></div><div><span style="font-size:12.8px">Sunetra Das. </span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 16 August 2017 at 23:07, Lorenzo Paulatto <span dir="ltr"><<a href="mailto:lorenzo.paulatto@impmc.upmc.fr" target="_blank">lorenzo.paulatto@impmc.upmc.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On 16/08/17 18:36, Sunetra Das wrote:<br>
>  those equation fitting match the given experimental value, not even<br>
> close.<br>
<br>
</span>Dear Sunetra,<br>
Is your structure satisfactory? Can you reproduce the experimental<br>
geometry and volume?<br>
<br>
If not, you won't be ablt to get a good bulk modulus.<br>
<br>
If yes, try increasing the ecutwfc quite a bit, or run your calculations<br>
at constant number of plane waves, which if I remember correctly means<br>
changing ecutwfc ~ alat^-2<br>
<br>
hth<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Dr. Lorenzo Paulatto<br>
IdR @ IMPMC -- CNRS & Université Paris 6<br>
phone: +33 (0)1 442 79822 / skype: paulatz<br>
www:   <a href="http://www-int.impmc.upmc.fr/~paulatto/" rel="noreferrer" target="_blank">http://www-int.impmc.upmc.fr/~<wbr>paulatto/</a><br>
mail:  23-24/423 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a></div></div></blockquote></div><br></div>