<div dir="ltr">Thanks Eric and Hongsheng for the discussion!<br><br>Can I also have your comments on this &quot;common tangent&quot; concept in first order phase transitions? Is it applicable for my case?<br><br><br>And how can I use phonon calculations to help me find the phase transition? I&#39;ve heard about some &quot;mode softening&quot;, but I am not sure how it can be used in this context? Could you refer me to some article explaining the theory behind this<br>

<br>Thanks again <br>Amin<br><br><br><br><br><br><div class="gmail_quote">On Mon, Jul 18, 2011 at 9:28 PM, Hongsheng Zhao <span dir="ltr">&lt;<a href="mailto:zhaohscas@yahoo.com.cn">zhaohscas@yahoo.com.cn</a>&gt;</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;"><div class="im">On 07/19/2011 06:49 PM, Eric Germaneau wrote:<br>
&gt;<br>
&gt;<br>
&gt; On 07/18/2011 10:45 AM, Hongsheng Zhao wrote:<br>
&gt;&gt; On 07/19/2011 09:23 AM, Eric Germaneau wrote:<br>
&gt;&gt;&gt; Hey Amin,<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; The approach you proposed is ok but you should relax the cell as well.<br>
&gt;&gt;&gt; By keeping the cell fix the system does not feel the pressure you wish<br>
&gt;&gt;&gt; to apply.<br>
&gt;&gt; Not so clear about the meaning of this sentence, could you please give<br>
&gt;&gt; some more explanations? ;-(<br>
&gt; Under pressure keeping the unit cell fix is useless since the system<br>
&gt; won&#39;t feel it.<br>
&gt; It just makes sense.<br>
&gt; So, if you wish to see the effect of the pressure you have to perform<br>
&gt; vc-relax calculation.<br>
&gt; Take it that way, in thermodynamic you have 2 important ensembles,<br>
&gt; isochoric and isobaric.<br>
&gt; So either the volume is fix and you get out the corresponding pressure,<br>
&gt; or the pressure is fix and you get out the corresponding volume.<br>
&gt; Note that the volume of gas for instance decreases inversely<br>
&gt; proportionally to pressure.<br>
&gt; The pressure is the intensive variable related to the volume (extensive).<br>
&gt; You simply can not have both of them constant at the same time.<br>
<br>
</div>Yes, you&#39;re right.  I&#39;m so sorry for my previous misleading posts on<br>
this thread.  Thanks again.<br>
<div class="im"><br>
&gt; Hope it helps.<br>
&gt;&gt;<br>
&gt;&gt;&gt; You can also change the volume, relax the positions only and compute the<br>
&gt;&gt;&gt; pressure.<br>
&gt;&gt;&gt; The best is to compute the free energy as function of pressure ans so<br>
&gt;&gt;&gt; the QHA method is the best way to go.<br>
&gt;&gt; I think this method cann&#39;t do the job of structure-searching-and-finding<br>
&gt;&gt; for unknown phases when we changing the external pressure.<br>
&gt; Yes, it can but in the case of structure-searching it&#39;s very expensive.<br>
&gt; In that case you&#39;d better generate random unit cell parameters, random<br>
&gt; positions and run vc-relax.<br>
&gt; Even this can be expensive since you have to do it a lot of times ....<br>
&gt; Just take a look to this papers, PRL 97 1 (2006)<br>
</div>&gt; &lt;<a href="http://link.aps.org/doi/10.1103/PhysRevLett.97.045504" target="_blank">http://link.aps.org/doi/10.1103/PhysRevLett.97.045504</a>&gt; and JPCM 23<br>
&gt; 053201 (2011)<br>
&gt; &lt;<a href="http://stacks.iop.org/0953-8984/23/i=5/a=053201?key=crossref.44087776832bc9a0edd0dc492a67d376" target="_blank">http://stacks.iop.org/0953-8984/23/i=5/a=053201?key=crossref.44087776832bc9a0edd0dc492a67d376</a>&gt;.<br>


<div class="im">&gt; It is the main work of Professor Pickard<br>
</div>&gt; &lt;<a href="http://www.st-andrews.ac.uk/physics/condmat/pickard/Crystals%20ex%20nihilo.html" target="_blank">http://www.st-andrews.ac.uk/physics/condmat/pickard/Crystals%20ex%20nihilo.html</a>&gt;.<br>
&gt; Best.<br>
<br>
Thanks a lot for this information ;-)<br>
<br>
Regards<br>
<div class="im">&gt;<br>
&gt;&gt;<br>
&gt;&gt; Regards<br>
&gt;&gt;<br>
&gt;&gt;&gt; Note that at first you can ignore the thermal effect and compute the<br>
&gt;&gt;&gt; free energy from a regular phonon calculations.<br>
&gt;&gt;&gt; I&#39;m also very interested in such calculations but don&#39;t have experiences<br>
&gt;&gt;&gt; doing this by mean of first principal calculations.<br>
<br>
<br>
<br>
</div><div class="im">--<br>
Hongsheng Zhao &lt;<a href="mailto:zhaohscas@yahoo.com.cn">zhaohscas@yahoo.com.cn</a>&gt;<br>
School of Physics and Electrical Information Science,<br>
Ningxia University, Yinchuan 750021, China<br>
_______________________________________________<br>
</div><div><div></div><div class="h5">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>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr">Amin Torabi<br>Ph.D. Student<br>Chemistry Department<br><i>The</i> University <i>of</i> Western Ontario<br>London, On Canada, N6A 5B7<br><span style="font-family: &#39;Droid Sans&#39;,arial,sans-serif; font-size: 13px; border-collapse: collapse; color: rgb(32, 32, 32);">Phone: <a value="+15196612111" style="color: rgb(103, 117, 58);">519-661-2111</a> Ext: 87871</span><div>

<span style="font-family: &#39;Droid Sans&#39;,arial,sans-serif; font-size: 13px; border-collapse: collapse; color: rgb(32, 32, 32);">****************************************</span></div></div><br>
</div>