[QE-users] Optical branches are not continuous at gamma point in phonon dispersion

[Changpeng Lin]

Dear all,

I calculated the phonon dispersion of polar materials by DFPT in ph.x and
set  epsil = .true. and trans = .true. to determine the dielectric tensor
and Born effective charges.

I then get the phonon dispersion by matdyn.x and the input set as following,
&input
        flfrc='TE.frc',
        flfrq='TE.frq'
        asr='simple',
        q_in_band_form=.true.
        q_in_cryst_coord=.true.
/
6
0.5  0.5 -0.5  100   !Z
0.0  0.0  0.0  100   !G
0.0  0.0  0.5  100   !X
0.25 0.25 0.25 100   !P
0.0  0.5  0.0  100   !N
0.0  0.0  0.0  1     !G

I found that optical branches are not continuous at gamma point (See figure
below). I guess this is becasue the wrong LO-TO spiliting relation I got.
How can I check it?
[image: image.png]

Thanks,
Changpeng
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20181114/09a61ddf/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 39545 bytes
Desc: not available
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20181114/09a61ddf/attachment.png>