<div dir="ltr"><div><div>Hi Matt,</div><div><br></div><div>Thanks for the suggestion. I'm hoping to calculate the correction periodically</div><div>during the simulation. For my system, the excess charge will be fairly</div><div>localized, but I'm wondering why the correction will fail if this is not the</div><div>case. Thanks so much!</div><div><br></div><div>Best regards,</div><div><br></div><div>Peng</div><br>On Thursday, September 21, 2017 at 3:06:55 AM UTC-5, Matt W wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">You might want to output the electron density on a 3D grid<div><br></div><div><a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/E_DENSITY_CUBE.html" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmanual.cp2k.org%2Ftrunk%2FCP2K_INPUT%2FFORCE_EVAL%2FDFT%2FPRINT%2FE_DENSITY_CUBE.html\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGblPgGOkY7pcJMDJZsT8fMkhm20Q';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmanual.cp2k.org%2Ftrunk%2FCP2K_INPUT%2FFORCE_EVAL%2FDFT%2FPRINT%2FE_DENSITY_CUBE.html\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGblPgGOkY7pcJMDJZsT8fMkhm20Q';return true;">https://manual.cp2k.org/trunk/<wbr>CP2K_INPUT/FORCE_EVAL/DFT/<wbr>PRINT/E_DENSITY_CUBE.html</a></div><div><br></div><div>Probably you want to do this for the ideal and defect system, and subtract the two densities (<a href="https://www.cp2k.org/tools:cubecruncher" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Ftools%3Acubecruncher\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFASVQzWwPNH3VNfrO99E7VKQRsNA';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Ftools%3Acubecruncher\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFASVQzWwPNH3VNfrO99E7VKQRsNA';return true;">https://www.cp2k.org/tools:<wbr>cubecruncher</a> might be useful).</div><div><br></div><div>Hopefully this would give you a fairly localised density that you can then calculate the dipole / quadrupole moments of. If it isn't localised, then probably the correction schemes break down anyway ...</div><div><br></div><div>Matt <br><br>On Saturday, September 9, 2017 at 1:25:56 PM UTC+8, Peng wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><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></blockquote></div></div></blockquote></div></div>