[Pw_forum] dipole moment calculation

[Aihua Zhang]

Dear all,

I tried a test calculation of the dipole moment of water. The water molecule
is put into a 20 a.u. ^3 box. The relevant setting are

lelfield = .true.
efield = 0.01

In the output file, it reports

     Electronic Dipole per cell (a.u.)   -1.921145549061207
     Ionic Dipole per cell (a.u.)    116.1997104738884

How is Electronic Dipole defined here? By looking into the source code,
Ionic Dipole is defined as in textbook. I also notice that total energy
difference with respect to zero-electric-field system is approximately equal
to (Electronic Dipole + Ionic Dipole) * Efield. But Electronic Dipole hardly
resembles the experimental dipole of water, which is about 1.85 D. Could
anyone explain a bit about these quantities? Thank you!!

Zhang

-- 
------------------------------------------------------------------
Research Fellow,
Dept. of Physics, National Univ. of Singapore
10 Kent Ridge Crescent, Singapore, 119260
Tel: +65 6516 7844
------------------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.democritos.it/pipermail/pw_forum/attachments/20090309/e0743fab/attachment.htm