Thank all of you very much for your help.<div><br /></div><div>Concerning what prof. Sitangshu suggestion: the raising of the cutoff over plane waves works, now I can recover the right </div><div>(semi)quadratic dispersion for the acoustic out-of-plane phonons. I used smeared occupation because if not</div><div>effective charges are computed by ph.x, giving rise to the LO-TO splitting at long wavelength (it is just a workaround). </div><div><br /></div><div>I will try, as a double check, to use the lattice parameter of the bulk system, as prof. Paulatto suggested. </div><div>Regarding the residual strain, of course it vanishes when I increase the cut off performing vc-relax from scratch. </div><div>I guess this is why my results was not as I expected.</div><div><br /></div><div>Finally my question about the possibility of obtain the correct LO behaviour at low q was related to the article</div><div>that prof. Marzari reported. I'm glad to know that the corrections will be implemented in the new version of the code.</div><div><br /></div><div>Thank you very much again!</div><div><br /></div><div><br /></div>