<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Dec 8, 2015 at 10:22 AM, Eleftheria Gkogkosi <span dir="ltr"><<a href="mailto:elefthe@mail.ntua.gr" target="_blank">elefthe@mail.ntua.gr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"> <br> For test purposes, I run F_QHA.x for the cases of alpha- and beta- Tin.<br>
 I can see that "F_vibration vs T" graphs cross above 300 K (in<br>
 consistency with the phase transition).<br>
 This does not hold for the sum E_internal+F_vibration (which I assumed<br>
 to be the Heelmholtz free energy in the first place).<br>
 So, is it safe to conclude that F_vibration is actually the Helmholtz<br>
 free energy?</blockquote><div><br></div><div>from the source code, it looks like E_internal  is the vibrational contribution to the energy, while F_vibration is the vibrational contribution to the free energy. In order to locate the phase transition, you should compare E_tot+F_vibration(T) for alpha and beta tin, where E_tot is what you get from the scf calculation.<br><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"> What does it mean when one gets negative frequencies in tjis file?<br>
 Is it a matter of poor k-point sampling or an indication of unstable<br>
 structure at a given Temperature?<br></blockquote><div><br></div><div>short answer: it depends. Sometimes the latter, sometimes it is simply an unavoidable problem of violation of the Acoustic Sum Rule in approximate calculations.<br> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
 Is it fair to erase the negative parts and proceed withy QHA for the<br>
 thermodynamical properties?<br></blockquote><div><br></div><div>if the only negative values are those corresponding to zero-frequency acoustic modes at q=0, yes, it is <br><br></div><div>Paolo<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
 Thank you very much,<br>
<span class=""><font color="#888888">--<br>
 Eleftheria Gkogkosi<br>
<br>
 Applied Mathematics & Physics Dept.<br>
 National Technical University of Athens<br>
 <a href="http://users.ntua.gr/elefthe" rel="noreferrer" target="_blank">http://users.ntua.gr/elefthe</a><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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</font></span></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div></div></div>
</div></div>