<div dir="ltr"><br><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Kondaiah Samudrala</b> <span dir="ltr"><<a href="mailto:konda.physics@gmail.com">konda.physics@gmail.com</a>></span><br>
Date: Tue, Dec 17, 2013 at 6:49 PM<br>Subject: Phonon dispersion with LO-TO splitting<br>To: <a href="mailto:pw_forum-owner@pwscf.org">pw_forum-owner@pwscf.org</a>, <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>, <a href="mailto:pw_forum-request@pwscf.org">pw_forum-request@pwscf.org</a><br>
<br><br><div dir="ltr"><div class="gmail_quote"><br><div dir="ltr">Dear all, <div><br></div><div>I am phonon dispersion for teragonal system to check insulator to metal transition. At gamma point frequencies,  I did not get any imaginary frequencies in asr.out file, where as in complete phonon dispersion I got imaginary frequencies and sudden shifts in high symmetry directions. I am bit confused and croos check my i/p. Below I attached the i/p and o/p of calculations</div>


<div><br></div><div>AB2.asr.out :</div><div><div>     Reading Dynamical Matrix from file AB2.dynG</div><div>     ...Force constants read</div><div>     ...epsilon and Z* read</div><div>     A direction for q was not specified:TO-LO splitting will be absent</div>


<div><br></div><div>     Polarizability (A^3 units)</div><div>     multiply by 0.273227 for Clausius-Mossotti correction</div><div>        21.410069    0.000000    0.000000</div><div>         0.000000   21.410069    0.000000</div>


<div>         0.000000    0.000000   18.627445</div><div><br></div><div>     IR activities are in (D/A)^2/amu units</div><div><br></div><div># mode   [cm-1]    [THz]      IR</div><div>    1      0.00    0.0000    0.0000</div>


<div>    2      0.00    0.0000    0.0000</div><div>    3      0.00    0.0000    0.0000</div><div>    4    257.23    7.7116    2.1831</div><div>    5    257.23    7.7116    2.1831</div><div>    6    412.71   12.3727    0.0000</div>


<div>    7    412.71   12.3727    0.0000</div><div>    8    539.77   16.1819   24.4707</div><div>    9    620.07   18.5891    0.0000</div><div> </div><div>     DYNMAT       :     0.00s CPU         0.03s WALL</div><div><br>


</div><div> </div><div>   This run was terminated on:  22:47:55  16Dec2013            </div><div><br></div><div>=------------------------------------------------------------------------------=</div><div>   JOB DONE.</div>


<div>=------------------------------------------------------------------------------=</div><div><br></div><div><a href="http://AB2.ph.in" target="_blank">AB2.ph.in</a></div><div><br></div><div><div>phonons at Full</div>
<div> &inputph</div><div>   prefix='AB2',</div><div>   fildyn='AB2.dyn',</div><div>   tr2_ph=1.0d-16,</div><div>   outdir='./tmp'</div><div>   ldisp=.true.,</div><div>   nq1=4, nq2=4, nq3=4</div>


<div> /</div><div><br></div><div><br></div><div>My lattice parameters are a= 3.4 and c=6.1 \AA . Pls correct me If I did any wrong...</div><div><br></div><div>One more doubt is.. How to include LO-TO splitting in Complete phonon dispersion???? I am doing Insulator to metal transition, and metal did not contains LO-TO splitting. So, I  check LO-TO changes at transition pressure and the same I want to compare with ambient pressure . </div>


<div><br></div><div>Looking for favorable reply and Thanks in advance...</div><div><br></div><div><br></div><div>with regards</div><div>S. Appalakondaiah</div><div>Research scholar<br>
</div><div>University of Hyderabad</div><div>Hyderabad</div><div>India</div><div><br></div></div></div></div>
</div><br></div>
</div><br></div>