Hello pwscf users<br>&nbsp;I want to find the real space interatomic force constants in an insulator excluding the contribution from the macroscopic electric field. So i am calculating the phonon dispersions by specifying .epsil.=&#39;false&#39;. And then i am using q2r.x to find out the real space force constants. Then when i use matdyn.x with zasr=&#39;simple&#39; to calculate the phonon dispersions , i can see LO-TO splitting at the gamma point. What am i missing? Thanks ! <br>