Dear Dr. Stefano and Dr.  Nicola Marzari thanks for your quick reply.<br>Actualy i have calculated the Helmholtz free energy (F) which is given as<br><br> F(V, T) = E0(V) + Fvib (V, T).<br><br>The Fvib(V, T) is calculated by PH/fqha.x utility and E0(V) is the equilibrium ground state energy which is calculated by pw.x.<br>
it is stated in a research paper>>>>For a given T and V, the equilibrium state of the crystal<br>is determined by minimising F with respect to all possible<br>degrees of freedom. The equation of state (pressure (P) versus<br>
V) of the system is obtained by equating P to minus the<br>derivative of F with respect to V at constant T, or P = − (∂F/∂V)<font size="1">T</font> <<<<<<<br><br>Sir i was tried for QHA package in QE as given paper <a href="http://arxiv.org/abs/1112.4977v1" target="_blank">http://arxiv.org/abs/1112.4977v1</a> , it is using the QHA pakage which is not well understood by me and also this was giving wrong result but  it was ok with PH/fqha.x. <br>
Please show me the ray of light regarding Helmholtz free energy (F) which is the base root of all thermodynamics property describe by one paper<span style="background-color:rgb(255,102,102)"> <a href="http://iopscience.iop.org/0953-8984/24/3/035401">http://iopscience.iop.org/0953-8984/24/3/035401</a> </span>for your kind reference.<br>
<br>Any comment in any form is appreciable.<br><br><br><br><br><br><div class="gmail_quote">2012/3/24 Nicola Marzari <span dir="ltr"><<a href="mailto:nicola.marzari@epfl.ch">nicola.marzari@epfl.ch</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
Here it is:<br>
<br>
<a href="http://arxiv.org/abs/1112.4977v1" target="_blank">http://arxiv.org/abs/1112.<u></u>4977v1</a><br>
Thermal properties of materials from ab-initio quasi-harmonic phonons<br>
Stefano Baroni, Paolo Giannozzi, Eyvaz Isaev<br>
<br>
                        nicola<div><div class="h5"><br>
<br>
<br>
On 24/03/2012 05:35, bramha pandey wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
Dear All,<br>
I want to incorporate the temperature effect in my calculation. As i<br>
know for this purpose firstly i have to compute the free energy and then<br>
this free energy is added to the ground state total energy(T=0,P=0).<br>
then by the equation of state we have to find out the lattice parameter.<br>
But i am wonder if i am trying to fit eos, it will give the equilibrium<br>
lattice parameter not the temperature dependent lattice parameter.<br>
Please let me know whether i am proceeding wrong way if so then how can<br>
i obtained my lattice parameter which is dependent upon temperature?<br>
Any type of help is highly appreciable .<br>
<br>
--<br>
Thanks and Regards<br>
Bramha Prasad Pandey<br>
GLA University, Mathura.<br>
INDIA.<br>
<br>
<br>
<br></div></div><div class="im">
______________________________<u></u>_________________<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/<u></u>mailman/listinfo/pw_forum</a><br>
</div></blockquote><span class="HOEnZb"><font color="#888888">
<br>
<br>
-- <br>
<br>
------------------------------<u></u>------------------------------<u></u>----------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
</font></span></blockquote></div><br><br clear="all"><br>-- <br>Thanks and Regards<br>Bramha Prasad Pandey<br>GLA University, Mathura.<br>INDIA.<br><br>