<div dir="ltr"><div>It seems to me that the printout is misleading: the first quantity should be "polarization", not "dipole", and its units are not Debye</div><div><br></div><div>Paolo<br></div><div><br><div class="gmail_quote"><div dir="ltr">---------- Forwarded message ---------<br>From: <b class="gmail_sendername" dir="auto">Ari P Seitsonen</b> <span dir="ltr"><<a href="mailto:Ari.P.Seitsonen@iki.fi">Ari.P.Seitsonen@iki.fi</a>></span><br>Date: Mon, Nov 19, 2018 at 5:32 PM<br>Subject: Re: [QE-users] dipole moment<br>To: Quantum Espresso users Forum <<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>><br></div><br><br><br>
Dear Sabike, [your affiliation please?]<br>
<br>
   If I understand correctly your question, looking in the code in <br>
'add_efield.f90':<br>
<br>
...<br>
      tot_dipole  = -el_dipole + ion_dipole<br>
...<br>
<br>
   IF ( iverbosity > 0 ) THEN<br>
     WRITE( stdout, '(8X,"Elec. dipole ",1F15.4," Ry au, ", 1F15.4," Debye")' ) &<br>
                                             el_dipole, (el_dipole*au_debye)<br>
     WRITE( stdout, '(8X,"Ion. dipole  ",1F15.4," Ry au, ", 1F15.4," Debye")' ) &<br>
                                           ion_dipole, (ion_dipole*au_debye)<br>
   END IF<br>
<br>
   WRITE( stdout, '(8X,"Dipole       ",1F15.4," Ry au, ", 1F15.4," Debye")' ) &<br>
                   (tot_dipole* (omega/fpi)), ((tot_dipole* (omega/fpi))*au_debye)<br>
<br>
   WRITE( stdout, '(8x,"Dipole field ", 1F15.4," Ry au, ")') tot_dipole<br>
<br>
<br>
Do those lines help? So something like<br>
"Dipole = (-el_dipole + ion_dipole) * volume / (4 pi)". Was this what you <br>
were looking for?<br>
<br>
     Greetings from Lviv,<br>
<br>
        apsi<br>
<br>
-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-<br>
   Ari Paavo Seitsonen / <a href="mailto:Ari.P.Seitsonen@iki.fi" target="_blank">Ari.P.Seitsonen@iki.fi</a> / <a href="http://www.iki.fi/~apsi/" rel="noreferrer" target="_blank">http://www.iki.fi/~apsi/</a><br>
     Ecole Normale Supérieure (ENS), Département de Chimie, Paris<br>
     Mobile (F) : +33 789 37 24 25    (CH) : +41 79 71 90 935<br>
<br>
<br>
On Sat, 17 Nov 2018, Sabike Ghasemi wrote:<br>
<br>
> hi all,<br>
> i am calculating dipole moment of water by pw.x.<br>
> i do it correctly.but i don't understand "Elec. dipole" and " ion. dipole" in out put and relation them wit<br>
> dipole. can i help me ?<br>
> input :<br>
> &CONTROL<br>
>     calculation  = "scf"<br>
>     dipfield     = .TRUE.<br>
>     disk_io      = "default"<br>
>     iprint       = 1<br>
>     max_seconds  =  8.64000e+04<br>
>     nstep        = 100<br>
>     outdir       = "./tmp1/"<br>
>     prefix       = "H2Otest"<br>
>     pseudo_dir   = "./"<br>
>     restart_mode = "from_scratch"<br>
>     tefield      = .TRUE.<br>
>     title        = "H2Otest"<br>
>     tprnfor      = .TRUE.<br>
>     verbosity    = "high"<br>
> /<br>
> <br>
> &SYSTEM<br>
>     a           =  10<br>
>     degauss     =  1.00000e-02<br>
>     ecutrho     =  520<br>
>     ecutwfc     = 130<br>
>     ibrav       = 1<br>
>     nat         = 3<br>
>     ntyp        = 2<br>
>     occupations = "smearing"<br>
>     smearing    = "gaussian"<br>
>   edir        = 1<br>
>     emaxpos     =  0.9<br>
>     eopreg      =  0.2<br>
>  eamp        =  0.00000e+00<br>
> input_dft = "PBE"<br>
> vdw_corr = "Grimme-D2"<br>
> /<br>
> <br>
> &ELECTRONS<br>
>     conv_thr         =  1.00000e-06<br>
>     electron_maxstep = 200<br>
>     mixing_beta      =  7.00000e-01<br>
>     startingpot      = "atomic"<br>
>     startingwfc      = "atomic+random"<br>
> /<br>
> <br>
> K_POINTS {automatic}<br>
>  3  3  3  0 0 0<br>
> <br>
> ATOMIC_SPECIES<br>
> O      15.99940  O.pbe-rrkjus.UPF<br>
> H       1.00794  H.pbe-rrkjus.UPF<br>
> <br>
> ATOMIC_POSITIONS {angstrom}<br>
> O       5   5   5<br>
> H       5   4.209   5.563<br>
> H       5   5.791   5.563<br>
> <br>
> <br>
> <br>
> output:<br>
>  <br>
> <br>
>      iteration #  6     ecut=   130.00 Ry     beta=0.70<br>
>      Davidson diagonalization with overlap<br>
>      c_bands:  1 eigenvalues not converged<br>
>      ethr =  7.81E-08,  avg # of iterations =  6.6<br>
> <br>
>      negative rho (up, down):  0.438E-04 0.000E+00<br>
> <br>
>      Adding external electric field<br>
> <br>
>      Computed dipole along edir(3) :<br>
>         Elec. dipole          0.0027 Ry au,           0.0068 Debye<br>
>         Ion. dipole           0.0040 Ry au,         0.0101 Debye<br>
>         Dipole                0.6902 Ry au,          1.7543 Debye<br>
>         Dipole field          0.0013 Ry au<br>
> <br>
>         Potential amp.       -0.0389 Ry<br>
>         Total length         15.1178 bohr<br>
> <br>
> <br>
>      total cpu time spent up to now is      152.0 secs<br>
> <br>
>      End of self-consistent calculation<br>
> <br>
> <br>
> <br>
> best regards,<br>
> sabike ghasemi<br>
> phd student<br>
> <br>
>_______________________________________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div></div></div>