<div dir="ltr"><div dir="ltr">Hi,<div><br></div><div>Thanks for your help. According to the suggestions, I have tested the convergence of dielectric constants and effective charges with respect to k mesh and they in fact converged with the size of k mesh.</div><div><br></div><div>As Prof. Marzari indicated, the phonon dispersion including LO-TO splitting might have the discontinuity becase it depends how q vector appoaching gamma point. It should be corrent. I also found some example of this kind unusal phonon dispersion in literature.</div><div><br></div><div>Thanks,</div><div>Changpeng</div></div></div><br><div class="gmail_quote"><div dir="ltr">Marzari Nicola <<a href="mailto:nicola.marzari@epfl.ch">nicola.marzari@epfl.ch</a>> 于2018年11月14日周三 下午4:32写道：<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
Note you can converge to different limits depending on direction, so a discontinuity in a plot where the direction changes angle is possible. Nicola<br>
<br>
Sent from a tiny keyboard... Contact info:<br>
<a href="http://theossrv1.epfl.ch/Main/Contact" rel="noreferrer" target="_blank">http://theossrv1.epfl.ch/Main/Contact</a><br>
<br>
On 14 Nov 2018, at 09:29, Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com" target="_blank">paulatz@gmail.com</a>> wrote:<br>
<br>
>> 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?<br>
> <br>
> It is definitely because of LO-TO splitting, if it is wrong or not, well you can check how it converges with k-points. Dielectric constants and effective charges can be VERY slow to converge with k-points. Just repeat the phonon calculation at Gamma, no need to redo all the grid each time.<br>
> <br>
> Or if you think the should be no LO-TO splitting at all (typically because of doping), just edit the first dyn matrix file (the one at gamma, i.e. dyn1) and remove all the section about epsilon and zeu, from<br>
>   "Dielectric Tensor:"<br>
> included, down to<br>
>   "Diagonalizing the dynamical matrix"<br>
> excluded.<br>
> <br>
> hth<br>
> <br>
> <br>
> -- <br>
> Lorenzo Paulatto - Paris<br>
> _______________________________________________<br>
> users mailing list<br>
> <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
_______________________________________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
</blockquote></div>