<div dir="ltr">Hello All<br>&nbsp;&nbsp; How should i calculate the real space interatomic force constants with the contribution due to the macroscopic electric field being excluded. According to Giannozzi et al (Ab initio calculation of phonons in semiconductors, Phys Rev B, v. 43, p.7231), the contribution to the dynamical matrix in polar materials can be written as a combination of an analytical term and a non-analytical term (at low k-vectors). The non-analytical term involves the Born effective charges and the dielectric constant. so, should i subtract this term from the dynamical matrix at all k-vectors (or only at certain k-vectors?) and then use q2r.x to obtain the real space force constants? Also, should i <b>not</b> specify .epsil. = true?&nbsp; I am new to pwscf so please bear with me...thanks! <br>
<br>Ravi Maranganti,<br>PhD Candidate, Mechanical Engg Dept, <br>Univ of Houston, TX <br></div>