Dear all,<div>I ran an ab initio molecular dynamics simulation of a molecule in solution in a cubic box, having 15.86 angstrom long side. </div><div><br /></div><div>I computed the three components of the total dipole of the box (Mu_x, Mu_y, Mu_z) in periodic conditions, thus using the Berry Phase approach, as reported in the manual at </div><div><br /></div><div><a href="https://manual.cp2k.org/trunk/CP2K_INPUT.html">CP2K_INPUT</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL.html">FORCE_EVAL</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT.html">DFT</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/LOCALIZE.html">LOCALIZE</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/LOCALIZE/PRINT.html">PRINT</a> / <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/LOCALIZE/PRINT/TOTAL_DIPOLE.html">TOTAL_DIPOLE</a></div><div>with PERIODIC TRUE</div><div><br /></div><div>The plot of the Mu_x, Mu_y and Mu_z in Debye as a function of the time (please find the plots attached) show jumps, which depend on the PBC. </div><div><br /></div><div>I read in the literature that one should "remain on the same branch of the polarization lattice" (

<i>A Beginner’s Guide to the Modern Theory of
Polarization, 

Nicola A. Spaldin</i> , https://arxiv.org/pdf/1202.1831.pdf ) and I think I should multiply my values by a factor, which depends on the box side (15.86 angstrom). </div><div>However, I cannot figure out what is this correction factor? How it depends on the box side?</div><div><br /></div><div>I would be very grateful if any of you could explain how to proceed.</div><div>Thank you very much in advance for your availability and kind support.</div><div><br /></div><div>Best regards,</div><div>Emma Rossi</div><div><br /></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/c3563bc4-535b-493f-921d-708b2ca9444bn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/c3563bc4-535b-493f-921d-708b2ca9444bn%40googlegroups.com</a>.<br />