<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<title></title>
</head>
<style type="text/css">
<!--
        #send_mail { 
        font-family: "Courier New";
        background:#FFFFFF;
        color:#000000;
        }

#send_mail a {
        color:#999999
        }

#send_mail a:hover {
        color:#0000FF
        }

#send_mail div { 
        color:#000000;
        font-family: "Courier New";
        font-size:14px;
        line-height:150%;
        }

#send_mail .usersign {
        line-height:100%;
        }
-->
</style>
<body id="send_mail">
<div>


<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<title></title>
<style type="text/css">
<!--
        #send_mail { 
        font-family: "Courier New";
        background:#FFFFFF;
        color:#000000;
        }

#send_mail a {
        color:#999999
        }

#send_mail a:hover {
        color:#0000FF
        }

#send_mail div { 
        color:#000000;
        font-family: "Courier New";
        font-size:14px;
        line-height:150%;
        }

#send_mail .usersign {
        line-height:100%;
        }
-->
</style>

<div>Dear Wang, H:
<br>
<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><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 "Dynamical Theory of Crystal Lattices" 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 "finite difference" you mentioned exactly the "frozen phonon" or "small displacement" method? If so, the treatment like that in PWSCF should be good to obtain a satisfied LO/TO splitting, although honestly I haven't done this in the "frozen phonon" 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>

<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>
</font></div>


</div>


</body>
</html>