Hello CP2K Community,<div><br /></div><div>I am performing DFT calculations for obtaining the total dipole moment of my non-periodic system. I want to compare the results of different approaches including Wannier Localization, Berry phase and Voronoi Integration. I find that the first two approaches give consistent moments which are different from the Voronoi moment. The results are as follows:</div><div><br /></div><div>1) WC Localization: [-6.860, -4.215, -2.087], magnitude= 8.318</div><div>2) Berry phase: [-6.863, -4.206, -2.117], magnitude= 8.324<br /></div><div>3) Voronoi integration: [-15.754, 35.152, -37.715], magnitude=53.910</div><div><br /></div><div>My approach for the moment calculation is as follows. The Berry phase moment is directly accessible from CP2K log/output file. For the Wannier method, CP2K outputs the charge positions in a .xyz file which I use along with the corresponding charge values. Finally, in the output "*.voronoi" file, there are multiple columns three of which are "Dipole(x) Dipole(y) Dipole(z)". I perform a summation of these three columns/components over the different atoms (rows) to get the total dipole moment. Moreover, the related input section is given below:</div><div><br /></div><div>&DFT</div><div>    &LOCALIZE #WC Localization<br />      &PRINT<br />        &WANNIER_CENTERS<br />           IONS+CENTERS .TRUE.<br />        &END<br />      &END<br />    &END<br />    &PRINT<br />      &MOMENTS # Berry phase<br />         PERIODIC F<br />      &END<br />      &VORONOI # Voronoi integration<br />         OUTPUT_TEXT T<br />         SANITY_CHECK T<br />         VERBOSE T<br />         VORONOI_RADII VDW<br />      &END<br />    &END<br /></div><div>&END</div><div><br /></div><div>I wonder how I can get a Voronoi Moment consistent with the other two approaches. I searched the CP2K forum but there was no mention of a similar experience. Thanks in advance for your help and comment!</div>

<p></p>

-- <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+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/717996fe-988b-4aca-8535-c3aa240da63bn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/717996fe-988b-4aca-8535-c3aa240da63bn%40googlegroups.com</a>.<br />