<div dir="ltr">Hi, I have another question on the same topic, <div><br></div><div>in the case of using REFERENCE COAC (center of atomic charges) for calculating the total dipole of a protonated water system (bulk water + hydronium),</div>
<div>if I understood it correctly the reference would be at the sum of the nuclear charges,</div><div>I was wondering if the reference would then jump together with the proton as it hops from one water molecule to the next.</div>
<div><br></div><div>Cheers,</div><div>Wagner.</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Oct 25, 2013 at 2:56 PM,  <span dir="ltr"><<a href="mailto:hut...@pci.uzh.ch" target="_blank">hut...@pci.uzh.ch</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi<br>
<br>
Berry phase total dipoles are calculated from the determinant<br>
of the exp(-ikr) operator. For the Wannier total dipole first<br>
localized orbitals are calculated. Then the exp(-ikr) operator<br>
matrices are transformed into this new basis. Finally, the<br>
determinant is approximated by the product of the diagonal elements.<br>
This allows to sparate the dipole into molecular contributions.<br>
<br>
The two ways of calculation should give very similar total dipoles.<br>
Are you sure that you compare values with the same units?<br>
There is also the problem of 'branch', meaning that the calculated dipole<br>
is given +/- 2Pi/L.<br>
<br>
regards<br>
<br>
Juerg Hutter<br>
<br>
<br>
<br>
--------------------------------------------------------------<br>
Juerg Hutter                         Phone : <a href="tel:%2B%2B41%2044%20635%204491" value="+41446354491">++41 44 635 4491</a><br>
Physical Chemistry Institute   FAX   : <a href="tel:%2B%2B41%2044%20635%206838" value="+41446356838">++41 44 635 6838</a><br>
University of Zurich               E-mail:  <a href="mailto:hut...@pci.uzh.ch">hut...@pci.uzh.ch</a><br>
Winterthurerstrasse 190<br>
CH-8057 Zurich, Switzerland<br>
---------------------------------------------------------------<br>
<br>
-----<a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a> wrote: -----<br>
To: <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a><br>
From: labicia<br>
Sent by: <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a><br>
Date: 10/24/2013 10:54AM<br>
Subject: [CP2K:4726] difference between wannier total dipole and berry total dipole<br>
<div class="HOEnZb"><div class="h5"><br>
Good morning!<br>
<br>
We are kind of new on this subject: may somebody tell us in a few words the difference between the wannier total dipole and berry total dipole.<br>
Or at least give us some reference to document ourselves?<br>
<br>
We are simulating a periodic bulk of methane.<br>
Are both methods suitable for this kind of system?<br>
<br>
The trend of the wannier total dipole and berry total dipole are almost the same, apart for a rescaling factor (total_wannier = 4.5 * total_berry + K), with K a constant.<br>
<br>
Thank you very much,<br>
and be patient if the question is a bit trivial,<br>
Marco and Andrea<br>
  --<br>
 You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
 To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
 To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
 Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br>
 For more options, visit <a href="https://groups.google.com/groups/opt_out" target="_blank">https://groups.google.com/groups/opt_out</a>.<br>
<br>
<br>
--<br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/groups/opt_out" target="_blank">https://groups.google.com/groups/opt_out</a>.<br>
</div></div></blockquote></div><br></div>