[CP2K-user] [CP2K:18448] Inconsistent dipole moment from PRINT and total density cube file

[Chiahsin Chu]

Thanks! The printed value indeed varies with the reference point.

Best,
Jiaxin

> On Feb 8, 2023, at 11:11, Jürg Hutter <hutter at chem.uzh.ch> wrote:
> 
> Hi
> maybe you want to experiment with setting a reference point for the dipole calculation.
> CP2K_INPUT / FORCE_EVAL / DFT / PRINT / MOMENTS / REFERENCE
> 
> regards
> JH
> 
> ________________________________________
> From: cp2k at googlegroups.com <mailto:cp2k at googlegroups.com> <cp2k at googlegroups.com <mailto:cp2k at googlegroups.com>> on behalf of Chiahsin Chu <xmuchiahsin at gmail.com <mailto:xmuchiahsin at gmail.com>>
> Sent: Tuesday, February 7, 2023 11:44 PM
> To: cp2k
> Subject: [CP2K:18419] Inconsistent dipole moment from PRINT and total density cube file
> 
> Hello,
> 
> Question:
> 
> I tried to print the dipole moment of my system via the section CP2K_INPUT/FORCE_EVAL/DFT/PRINT/MOMENTS (https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/MOMENTS.html), and got the result in the standard output file. However, this value does not match the number I got from either 1) total charge density cube file (multiply the grid and the charge density and then do the summation) or 2) the jump in the Hartree potential profile (the jump here should be proportional to the cell polarisation). The values from 1) and 2) are identical, ignoring a small numerical error. I have made the unit consistent. The z-component of dipole moment in the output file is ~10e3 Debye (i.e., ~10e2 eA), and the corresponding value from total charge density cube file is ~1 eA.
> 
> My system:
> 
> There is a Pt slab covered by a thin water layer. Vacuum layers are added on both sides of the system. The surface is parallel to XY plane, and the dipole correction is added on the z direction.
> 
> My setup:
> 
> I turn the PERIODIC in the MOMENTS section off, since Berry phase moments for non uniform MOs' occupation numbers not implemented.
> 
> 
> The full input file, atomic coordinates and output file I used are attached, for the sake of completeness.
> 
> Thanks,
> 
> Jiaxin
> 
> --
> You received this message because you are subscribed to the Google Groups "cp2k" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com <mailto:cp2k+unsubscribe at googlegroups.com><mailto:cp2k+unsubscribe at googlegroups.com>.
> To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/137556ca-43bb-4997-a30b-829f5c5d5034n%40googlegroups.com<https://groups.google.com/d/msgid/cp2k/137556ca-43bb-4997-a30b-829f5c5d5034n%40googlegroups.com?utm_medium=email&utm_source=footer> <https://groups.google.com/d/msgid/cp2k/137556ca-43bb-4997-a30b-829f5c5d5034n%40googlegroups.com%3Chttps://groups.google.com/d/msgid/cp2k/137556ca-43bb-4997-a30b-829f5c5d5034n%40googlegroups.com?utm_medium=email&utm_source=footer%3E>.
> 
> -- 
> You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.
> To unsubscribe from this topic, visit https://groups.google.com/d/topic/cp2k/vKoX0UWhEVs/unsubscribe.
> To unsubscribe from this group and all its topics, send an email to cp2k+unsubscribe at googlegroups.com <mailto:cp2k+unsubscribe at googlegroups.com>.
> To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/ZR0P278MB07594239309B371D1227CC5F9FD89%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM.

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/2EFA8AA5-1530-4B7A-84D5-13A7193E2CFB%40gmail.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230215/5ef209f4/attachment-0001.htm>