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><br>Rajan<br><br>