<div dir="ltr"><div><div>Hi Stefano,<br><br>Thanks for your reply. I am particularly looking to address an anomaly in my simulation. I've attached my scf input, ph input and q2r input files along with the generated force constant file.<br><br></div>The force constants are not symmetric (in fact anti-symmetric, for certain direction pairs) with respect to the displacement directions. For eg., if you look at line number 280 and line 800, the fcs are supposed to be equal as they are Phi_{XY} and Phi_{YX} for the same pair of atoms. I have tried this with several combinations of k-sampling and Q-sampling and I get the same anti-symetry (except when no. of Q's is very small, like Q=2). I also used a large ecutwfc so that there is no convergence issue. This problem was not there with Al lattice. The force constants were symmetric in this sense. Is there bug in the code or am I making a mistake in one of the input files? These are not negligible IFCs. Please have a look.<br><br></div><div>Thanks in advance!<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 10, 2014 at 2:44 AM, stefano de gironcoli <span dir="ltr"><<a href="mailto:degironc@sissa.it" target="_blank">degironc@sissa.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div bgcolor="#FFFFFF" text="#000000">
    <div>the IFC should be symmetric w.r.t.
      exchange of the two atoms and displacement directions as a
      consequence of being second order derivatives. other symmetries
      may apply in specific cases.<br>
      difficult to comment more without more info on what you are
      looking for.<br>
      stefano<div><div class="h5"><br>
      <br>
      <br>
      On 11/10/2014 09:16 AM, Navaneetha Krishnan wrote:<br>
    </div></div></div>
    <blockquote type="cite"><div><div class="h5">
      <div dir="ltr">
        <div>
          <div>Hello QE users,<br>
            <br>
          </div>
          I was investigating the IFCs generated by q2r.x after scf and
          ph calculations. It looks like the force constant matrix for
          certain pairs of atoms is not symmetric (with respect to the
          cartesian coordinates). For silicon, the IFC matrices have to
          be symmetric with respect to the cartesian coordinates -
          right? For aluminum, the output is symmetric, as expected. Any
          thoughts?<br>
        </div>
        <div>
          <div>
            <div><br>
            </div>
            <div>Thanks a lot!<br>
            </div>
            <div>--Navaneeth
              <div>
                <div dir="ltr">
                  <div>
                    <div>Caltech.</div>
                  </div>
                </div>
              </div>
            </div>
          </div>
        </div>
      </div>
      <br>
      <fieldset></fieldset>
      <br>
      </div></div><pre>_______________________________________________
Pw_forum mailing list
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a></pre>
    </blockquote>
    <br>
  </div>

<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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div>Navaneetha Krishnan Ravichandran,<div>Graduate Student,</div><div>Mechanical Engineering,</div><div>Caltech.</div></div></div></div>
</div>