<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#CCCCFF" text="#330000">
    Well, I don't know about "common tangent" concept.<br>
    You can study structure stability looking at phonon dispersion for
    instance.<br>
    An instability could mean that the structure undergoes a phase
    transition (or that it's not fully relaxed)<br>
    You can also compute the free energy of your system where phonons
    will give you the harmonic part (see <a
      href="http://link.aps.org/doi/10.1103/PhysRevB.64.045123">PRL 64
      045123 (2001)</a>).<br>
    So you can look at several thing when changing the pressure, among
    them structure change, phonon, energy, free energy (not that easy).<br>
    Why not electronic density of sates?<br>
    Good luck,<br>
    <br>
    &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &Eacute;ric. <br>
    <br>
    On 07/18/2011 09:46 PM, Amin Torabi wrote:
    <blockquote
cite="mid:CAEXRF5ycO_UN1nMKvbEdBqL7dNv+rwDK0+tSG9LcFHkhpvLnfA@mail.gmail.com"
      type="cite">
      <div dir="ltr">Thanks Eric and Hongsheng for the discussion!<br>
        <br>
        Can I also have your comments on this "common tangent" 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've heard about some "mode softening", 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 moz-do-not-send="true"
              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'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're right. &nbsp;I'm so sorry for my previous misleading
            posts on<br>
            this thread. &nbsp;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'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's very expensive.<br>
              &gt; In that case you'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 moz-do-not-send="true"
              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 moz-do-not-send="true"
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 moz-do-not-send="true"
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'm also very interested in such calculations
              but don'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 moz-do-not-send="true"
                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 class="h5">Pw_forum mailing list<br>
                <a moz-do-not-send="true"
                  href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
                <a moz-do-not-send="true"
                  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,&nbsp;N6A 5B7<br>
          <span style="font-family: 'Droid Sans',arial,sans-serif;
            font-size: 13px; border-collapse: collapse; color: rgb(32,
            32, 32);">Phone:&nbsp;<a moz-do-not-send="true"
              value="+15196612111" style="color: rgb(103, 117, 58);">519-661-2111</a>&nbsp;Ext:
            87871</span>
          <div>
            <span style="font-family: 'Droid Sans',arial,sans-serif;
              font-size: 13px; border-collapse: collapse; color: rgb(32,
              32, 32);">****************************************</span></div>
        </div>
        <br>
      </div>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <br>
      <pre wrap="">_______________________________________________
Pw_forum mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>
<a class="moz-txt-link-freetext" href="http://www.democritos.it/mailman/listinfo/pw_forum">http://www.democritos.it/mailman/listinfo/pw_forum</a>
</pre>
    </blockquote>
    <br>
    <div class="moz-signature">-- <br>
      <center>
        <font color="#946e8c">
          <em>
            Be the change you wish to see in the world<br>
          </em>
          <small>
            <font color="#946e8c">
              &#8212; Mahatma Gandhi &#8212;<br>
            </font></small>
          <br>
          <a href="mailto:germaneau@gucas.ac.cn">Dr. &Eacute;ric Germaneau</a><br>
          <br>
          College of Physical Sciences<br>
          Graduate University of Chinese Academy of Sciences<br>
          Yuquan Road 19A<br>
          Beijing 100049<br>
          China<br>
          <br>
          <small>
            <em>
              <font style="color: #498848" face="Times">
                Please, if possible, don't send me MS Word or PowerPoint
                attachments<br>
                Why? See: <a
                  href="http://www.gnu.org/philosophy/no-word-attachments.html"
                  style="text-decoration:none"><font color="#498848">http://www.gnu.org/philosophy/no-word-attachments.html</font></a>
              </font></em>
          </small>
        </font></center>
    </div>
  </body>
</html>