What do you mean by "centered" on the atoms? If you mean the electrostatic potential on a particular site, then you can use pymatgen's cube parser in the pymatgen.io.cube module. That will let you define a sphere around different atoms and get the total or average electrostatic potential in those sites. For something a little more robust, you can use the bader charge partitioning code developed by Henkelmann's group. With that you can get the bader volumes around each atom from the electron density file as the reference, and then apply it to the v_hartree file to get the the electrostatic potential in each bader volume.<br><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, June 17, 2021 at 8:39:28 AM UTC-7 mat...@gmail.com wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Unfortnunately, it is not available as an option at them moment. This has come up a few times on the forums and there have been suggestions / links to scripts that will interpolate on the cube files. Accuracy is not good enough for all applications though.<div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, 16 June 2021 at 14:08:39 UTC+1 Emily E wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear All,<div><br></div><div>The command <a href="https://manual.cp2k.org/trunk/index.html#CP2K_INPUT/FORCE_EVAL/DFT/PRINT/V_HARTREE_CUBE.html" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://manual.cp2k.org/trunk/index.html%23CP2K_INPUT/FORCE_EVAL/DFT/PRINT/V_HARTREE_CUBE.html&source=gmail&ust=1624040112960000&usg=AFQjCNExEjMHKZh9pTk6UwSY9EKSaObr0A">V_HARTREE_CUBE</a> prints out a cube file with electrostatic potentials. I am interested in generating electrostatic potentials that are centered on the atoms in the system. Is there a particular command to do this? On a related note, how does the stride keyword work or can this be used to print atom-centered electrostatic potentials?</div><div><br></div><div>Thanks!</div></blockquote></div></blockquote></div>