I forgot to add the affiliation in my previous post. Correcting myself, and regret any inconvenience.<br><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Rajan Pandey</b> <span dir="ltr"><<a href="mailto:rajanpandey@gmail.com">rajanpandey@gmail.com</a>></span><br>
Date: Sat, Sep 3, 2011 at 4:13 PM<br>Subject: Method of computing static dielectric constant of materials<br>To: PWSCF Forum <<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>><br><br><br>Dear Quantum Espresso community,<br>
<br>I am trying to compute the static dielectric constant of materials using the methodology <br>as implemented in CP code (part of the Quantum Espresso). The method is explained<br>in the example30 of Quantum Espresso distribution. The system discussed is MgO,<br>

Ref: P.Umari and A.Pasquarello, Physical Review Letters, 89, p.157602 (2002).<br><br>The example30 mentioned above has following note:<br><br>"NOTE: the electronic dipole is defined modulo a factor (2*L=31.824i a.u.,<br>

during the MD simulation the term "ln det S" changes the Riemann<br>plane, this must be taken into account when addressing the<br>electronic dipole."<br><br>The above statement is given in the context of a cubic supercell with length L containing 64 atoms of MgO.<br>

<br>I am studying a well known system, SiO2 (alpha quartz) for sanity check. I am simulating SiO2 in a 1x1x2 <br>trigonal supercell (18 atoms) as well as in an orthorhombic (in alpha quartz case it will be tetragonal because a = b ) <br>

supercell. The ambiguity is that when I use the formula along with the "NOTE" of example30 in Quantum Espresso,<br>mentioned above, I get wrong results for static dielectric constant. However, when I do not use the <br>

"NOTE" mentioned above, and by using the formula described in example30, I get the simulated value close <br>to the experimental static dielectric constant of SiO2. This means that the method works for non-cubic unit cells too.<br>

I am not able to understand the role of the point mentioned in the NOTE (above) in general, when the formula described<br>in example30 (reiterated below), seems to work. <br>  <br>The difference d_Eps  between static and high-frequency dielectric constant is given by:<br>

<br>d_Eps=4*pi*(D2_el + D2_ion - D1_el - d1_ion)/(0.001 a.u. * Omega)<br><br>D1(2)_el = electronic contribution of the dipole at the beginning (end) of the relaxation.<br>D1(2)_ion = ionic contribution of the dipole at the beginning (end) of the relaxation.<br>

Omega = supercell volume<br><br>I shall appreciate any comments from community members.<br><br>Thanks, and regards,<br><font color="#888888"><br>Rajan<br><br>Rajan K. Pandey, Ph.D.<br><br>Advisory Research Engineer,<br>Semiconductor Research & Development Center <br>
India Systems & Technology Engineering Lab<br>IBM India Pvt. Ltd.<br>MD3 1F B354<br>Manyata Embassy  Business Park<br>Nagawara, Outer Ring Road<br>Bangalore - 560045, India<br>Phone: +91-80-28061262 <br>Mobile: +91-9901850981<br>
Email: <a href="mailto:rajapand@in.ibm.com">rajapand@in.ibm.com</a><br>
</font></div><br>