<div dir="ltr">Dear Hongsheng,<div><br></div><div>I don&#39;t understand why I should use relax instead of vc-relax. As far as I understand, when a system is under pressure, both the atomic positions inside a unit cell, as well as the dimensions of the unit cell changes. So, why should I keep the unit cell parameters fixed if I am interested in finding phase transitions?</div>

<div><br></div><div>Thanks<br><br><div class="gmail_quote">On Mon, Jul 18, 2011 at 4:04 AM, 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">

<div class="im">On 07/18/2011 09:40 AM, Amin Torabi wrote:<br>
&gt; Dear experts!<br>
&gt;<br>
&gt; Forgive me for my naive question... but I am interested in<br>
&gt; pressure-induced phase transitions of diborane. Could you please have a<br>
&gt; look at my input file, attached to this message.<br>
&gt;<br>
&gt; The following is what I think I should do: (correct me if I&#39;m wrong!)<br>
&gt;<br>
&gt; 1. doing separate vc-relax calculation at each pressure<br>
<br>
</div>By using this method, there are some limitations.  You should have known<br>
in advance the different phases between them the phase transition maybe<br>
occur when we change the pressure.  Then you can use this method to<br>
determine the phase transition path.  But I&#39;m not so sure this method<br>
can predict all of the possible phases for a system when we change the<br>
external pressure.  To do a more reasonable and complete phase<br>
transition prediction and even searching for unknown phases, you should<br>
use some third-party code like uspex for this type of job.<br>
<br>
In addition, I think you should do relax, not vc-relax calculation at<br>
each pressure.  If you do vc-relax, then the pressure imposed to the<br>
system will be released.  So you won&#39;t get the *Pressure-induced* phase<br>
transition.<br>
<div class="im"><br>
&gt; 2. comparing the optimized cell_parameters at each pressure to see when<br>
&gt; the structure changes<br>
&gt;<br>
&gt; I am wondering how I can use the enthalpy or unit cell volume at each<br>
&gt; pressure to make sure that a phase transition has occurred<br>
<br>
</div>For the other questions, see this thread on this list:  About the method<br>
to obtain the common tangent for determing the phase transition<br>
pressure.  HTH.<br>
<br>
Regards<br>
<div class="im">&gt;<br>
&gt; Thanks in advance<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; Amin Torabi<br>
&gt; Ph.D. Student<br>
&gt; Chemistry Department<br>
&gt; /The/ University /of/ Western Ontario<br>
&gt; London, On Canada, N6A 5B7<br>
&gt; Phone: <a href="tel:519-661-2111%20Ext%3A%2087871" value="+15196612111">519-661-2111 Ext: 87871</a><br>
&gt; ****************************************<br>
&gt;<br>
&gt;<br>
&gt;<br>
</div>&gt; _______________________________________________<br>
&gt; Pw_forum mailing list<br>
&gt; <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
&gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<font color="#888888"><br>
<br>
--<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>
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>
</font></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></div>