<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Dear all,<div><br></div><div>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.</div><div><br></div><div>I then get the phonon dispersion by matdyn.x and the input set as following,</div><div><div>&input</div><div>        flfrc='TE.frc',</div><div>        flfrq='TE.frq'</div><div>        asr='simple',</div><div>        q_in_band_form=.true.</div><div>        q_in_cryst_coord=.true.</div><div>/</div></div><div>6</div><div><div>6</div><div>0.5  0.5 -0.5  100   !Z</div><div>0.0  0.0  0.0  100   !G</div><div>0.0  0.0  0.5  100   !X</div><div>0.25 0.25 0.25 100   !P</div><div>0.0  0.5  0.0  100   !N</div><div>0.0  0.0  0.0  1     !G</div></div><div><br></div><div>I found that optical branches are not continuous at gamma point. I guess this is becasue the wrong LO-TO spiliting relation I got. How can I check it?</div><div><br></div><div>More interestingly, I recently tested the phonon dispersion of PbTe by supercell small displacement method and also determine the dielectric tensor and Born effective charges by DFPT at gamma point only in ph.x. I found that there was no LO-TO splitting at gamma point from the output of dynmat.x. However, if I put the dielectric tensor and Born effective charges by DFPT into Phonopy and combine the force constants obtained from small displacement method, I get the correct LO-TO splitting from Phonopy. So, where is wrong?</div><div><br></div><div><br></div><div>Thanks,</div><div>Changpeng</div><div><br></div><div><br></div><div><br></div></div></div></div></div></div></div></div></div></div>