Dear Giovanni and Paolo,<br><br>Thank you very much for kind replies. I have reproduced Giovanni's result with tefield option and with the water molecule at origin. When the molecule is away from origin, weird results are produced. I look through the code in add_efield.f90 and compute_dip.f90. I think I found some possible bugs in them. The first one is not vital. In add_efield.f90, the variable, first, is set to false at the very begin of the subroutine. But this variable is used afterwards with an obvious intention that it should be true for the first time when the subroutine is called. With this small bug, the output of the sawlike potential setup is missing if the varialble, dipfield, is also true. The second possible bug is in compute_dip.f90, where r and r in rho(r) seem dismatch when calculating the electronic dipole. After matching r with rho(r) in the code, the reasonable identical results are obtains, i.e. the calculated dipole and total energy are independent of where the molecule is put in the supercell.<br>
<br>I have only limited understanding of MTP. If I understand correctly, it is emphasized that the convectional definition of polarization (P=\int r rho(r) d^3r) is ill-behaved for periodic system and only the change of P is meaningful. So is it a proper way to calculate static dipols with lberry and lelfield? BTW, what is the difference and relationship between lberry and lelfield, both dealing with MTP? Since the tefield approach is rather straighforward to me and produce reasonable result, I would stick to it for the time being.<br>
<br>Best regards<br><br>Aihua<br><br><div class="gmail_quote">On Wed, Mar 11, 2009 at 3:28 AM, <span dir="ltr"><<a href="mailto:umari@democritos.it">umari@democritos.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear Zhang,<br>
<div class="im"><br>
<br>
> I tried a test calculation of the dipole moment of water. The water molecule<br>
> is put into a 20 a.u. ^3 box. The relevant setting are<br>
><br>
> lelfield = .true.<br>
> efield = 0.01<br>
><br>
> In the output file, it reports<br>
><br>
> Electronic Dipole per cell (a.u.) -1.921145549061207<br>
> Ionic Dipole per cell (a.u.) 116.1997104738884<br>
<br>
> How is Electronic Dipole defined here? By looking into the source code,<br>
> Ionic Dipole is defined as in textbook. I also notice that total energy<br>
> difference with respect to zero-electric-field system is approximately equal<br>
> to (Electronic Dipole + Ionic Dipole) * Efield. But Electronic Dipole hardly<br>
> resembles the experimental dipole of water, which is about 1.85 D. Could<br>
> anyone explain a bit about these quantities? Thank you!!<br>
><br>
> Zhang<br>
<br>
<br>
</div>The option lelfield turns on the calculation of the electronic and<br>
ionic dipoles under the presence of an electric field.<br>
While the electronic dipole is calculated through the modern theory of<br>
the polarization formalism (MTP) which is valid also for<br>
periodic systems, the reported ionic dipole is given using the simple<br>
textbook formula. This is correct if you want to calculate dipole<br>
differences (e.g. after a structural relaxation or after the application<br>
of an electric field) but does not give directly the dipole of an<br>
isolated molecule. For calculating this, the ionic dipole should be<br>
calculated consistently through a MTP-like<br>
formalism:<br>
e.g.<br>
<br>
zeta_tot=0.d0<br>
DO na=1,nat<br>
zeta_tot=zeta_tot+(0.d0,1.d0)*tau(pdir,na)*tpi*zv(ityp(na))<br>
END DO<br>
pola_ion=dimag(log(exp(zeta_tot)))<br>
pola_ion=pola_ion*dsqrt(2.d0)/tpiba*dkfact<br>
write(stdout,*) " Ionic Dipole per cell (a.u.)",pola_ion<br>
<br>
Note that Ry atomic units are used and to compare with Hartree atomic units<br>
you should divide by a factor sqrt(2)<br>
<br>
Best regards,<br>
<br>
Paolo Umari CNR-INFM DEMOCRITOS<br>
<br>
----------------------------------------------------------------<br>
This message was sent using IMP, the Internet Messaging Program.<br>
<div><div></div><div class="h5"><br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>------------------------------------------------------------------<br>Research Fellow,<br>Dept. of Physics, National Univ. of Singapore<br>10 Kent Ridge Crescent, Singapore, 119260<br>
Tel: +65 6516 7844<br>------------------------------------------------------------------<br><br>