<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font color="#000000">Dear Q-E developers.</font><div><font color="#000000">I am a PhD student in University of Pennsylvania. I am using QE to compute the polarizability. In the code polarz.f90, the code compute the dielectric constant first and then the polarizability. I think this is using the <span style="background-color:transparent"><font face="Arial"><span style="white-space:pre-wrap">Clausius-Mossotti relation. In this relation, there is a very general assumption, which is, the</span></font></span></font><span style="font-family:sans-serif"> material is </span><font color="#000000"><span style="background-color:transparent"><font face="Arial"><span style="white-space:pre-wrap"> </span></font></span></font><font face="sans-serif">homogeneous (you can refer to Ashcroft book or I attached a paper about this). I happened to compute a material that is not homogenous and got some strange results. I just want to know is my understanding(you are using </font><span style="color:rgb(0,0,0);font-family:Arial;white-space:pre-wrap">Clausius-Mossotti relation</span><span style="font-family:sans-serif"> ) right? or you have other theorems but the same equations.</span></div><div><span style="font-family:sans-serif">Best wishes</span></div><div><span style="font-family:sans-serif">Jiahao</span></div><div><div><img src="cid:ii_jr9bvmb00" alt="image.png" width="452" height="130"><br></div><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr">Jiahao Zhang<div>PhD student</div><div>University of Pennsylvania</div></div></div></div></div></div></div></div></div></div></div>