<div dir="ltr"><div>Hi all,</div><div><br></div><div><div>I'm wondering if it's currently possible with CP2K to apply a correction for</div><div>non-neutral periodic system as in Eq. 14 of G. Makov & M. C. Payne, Phys Rev B</div><div>51 4014 (1995), which involves determining the center of charge as well as an</div><div>integration: \int_cell dr^3 \rho(r) r^2. I assume the first part is just related</div><div>to the total system dipole calculated with respect to the origin of the</div><div>coordinate system (CP2K_INPUT / FORCE_EVAL / DFT / PRINT / MOMENTS / REFERENCE =</div><div>0), but how should I go about computing the integral? If it's not directly</div><div>supported, can someone please give some quick pointers on where and how to hook</div><div>into the each MD step to implement this? Any suggestions would be really</div><div>helpful. Thank you very much!</div></div><div><br></div><div>Best regards,</div><div><br></div><div>Peng</div></div>