<span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 12.5px; border-collapse: collapse; "><div>Dear Huang:</div><div><span class="Apple-style-span" style="font-size: 12.5px; ">Thanks very much for your valuable discussion.</span></div>

<div>I are reading your mentioned paper:)</div></span><br><div class="gmail_quote">On Fri, Apr 29, 2011 at 5:20 PM, lfhuang <span dir="ltr">&lt;<a href="mailto:lfhuang@theory.issp.ac.cn">lfhuang@theory.issp.ac.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>
<div>






<div>Dear Wang, H:
<br><div class="im">
<font color="#444444">&gt;  Suppose a uniaxial crystal(no cubic): c&gt;a=b. Suppose there is a
<br>
</font><font color="#444444">&gt; double-degenerate polar E_u mode, in which all atoms move in the a-b plane.
<br>
</font><font color="#444444">&gt; If the wavevector around Gamma point is in the a-b plane, the degenerate E
<br>
</font><font color="#444444">&gt; mode will split into LO and TO mode.
<br>
</font><font color="#444444">&gt; Am I right?
</font></div></div><div><font color="#444444"><br>
It is probable for some ionic crystals, but still depends.
<br>
Maybe the Equ. 98 (nonanalytic part of the dynamical matrix) in Rev. Mod. Phys. 73, 515 (2001) by prof. S. Baroni etc. and the related parts in the book &quot;Dynamical Theory of Crystal Lattices&quot; by M. Born &amp; K. Huang can give you much information if these things are important/interesting to you.
<br>

<br>
In addition, is the method &quot;finite difference&quot; you mentioned exactly the &quot;frozen phonon&quot; or &quot;small displacement&quot; method? If so, the treatment like that in PWSCF should be good to obtain a satisfied LO/TO splitting, although honestly I haven&#39;t done this in the &quot;frozen phonon&quot; method. Some people propose other methods, like using a prolonged supercell, which I think, could not fully describe the long-range dipolar interaction, because of its long-range character (as prof. Stefano de Gironcoli said ~1/R^3).  
<br><div class="im">

<br>
Best Wishes!
<br>
Yours Sincerely
<br>
L. F. Huang
<br>
------
<br>
======================================================================
<br>
L.F.Huang(黄良锋) DFT and phonon physics
<br>
======================================================================
<br>
Add: Research Laboratory for Computational Materials Sciences,
<br>
Instutue of Solid State Physics,the Chinese Academy of Sciences,
<br>
P.O.Box 1129, Hefei 230031, P.R.China
<br>
Tel: 86-551-5591464-326(office)
<br>
Fax: 86-551-5591434
<br>
Our group: <a href="http://theory.issp.ac.cn/" target="_blank">http://theory.issp.ac.cn</a>
<br>
======================================================================
<br>
</div></font></div>


</div>


</div>

<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>
<br></blockquote></div><br><br clear="all"><br>-- <br>____________________________________<br>Hui Wang<br>School of physics, Fudan University, Shanghai, China<br>