<div>Hi Prof Lu,</div><div><br /></div><div>It is great news that it is possible to extract the reference ion densities used by MultiWFN for Hirshfield-I calculation. I am an appreciative MultiWFN user—in fact I use it in my own Hirshfeld-I implementation to create reference atom cube files on the appropriate grid—and with the reference atom densities in hand, it seems that I could use MultiWFN for the whole charge calculation, since I will be able to use those .rad files to reconstruct the weight function after the charge calculation.</div><div><br /></div><div>However, it seems that the MultiWFN reference atoms have the issue that I am trying to avoid: unattached electrons in the anions. I plotted the radial densities of the nitrogen references found in "examples/atmrad":</div><div><br /></div><div><img alt="nitrogen_rad_white.png" width="384px" height="384px" src="cid:06411367-4f14-406d-99a4-42797f5c0519" /></div><div><br /></div><div>As expected and desired, the cation densities get wider as electrons are added. However, the anion densities stay exactly the same all the way up to past an angstrom from the nucleus. The additional electron density all ends up out past the covalent radius of nitrogen.</div><div><br /></div><div>The reason for this is that N- and N2- are artificial. The electron affinity of nitrogen is near zero. N- is not really stable, or perhaps barely stable, and N2- is certainly not stable. As a result, the electrons are unattached in the reference densities.</div><div><br /></div><div>This issue was discussed extensively by Vanpoucke, Bultinck, and Van Driessche in their paper "<a href="https://doi.org/10.1002/jcc.23088">Extending Hirshfeld-I to Bulk and Periodic Materials</a>", and generated considerable discussion: <a href="https://doi.org/10.1002/jcc.23191">a comment by Thomas Manz</a>, followed by <a href="https://doi.org/10.1002/jcc.23193">a reply to that comment by the original authors</a>. It is not a simple issue with a neat solution, and the main reason I am looking around at Hirshfeld-I implementations other than my own is that I'm hoping that someone else has worked out a satisfactory alternative.</div><div><br /></div><div>In any case thanks for the information, and I love MultiWFN. Wishing you the best,</div><div>Alex<br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, January 29, 2025 at 6:42:56 AM UTC-5 Tian Lu 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;"><div>Hello, <br></div><div><br></div><div>I would like to mention that in the recent updates, Multiwfn (<a href="http://sobereva.com/multiwfn" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=http://sobereva.com/multiwfn&source=gmail&ust=1738270636915000&usg=AOvVaw1D5mlz-ZOsj-d1AztdIm8c">http://sobereva.com/multiwfn</a>) has supported calculating Hirshfeld-I charge and MBIS charge (another method based on iterative atomic spaces) based on the .molden file or electron density .cube file exported by CP2K. Both isolated and periodic systems are supported. For the Hirshfeld-I calculation, radial electron density files for almost every element at every charged states have been provided in Multiwfn package ("examples\atomrad\" subfolder) and can be directly used. I have a blog article to describe and illustrate how to perform the calculations using Multiwfn in combination with CP2K: <a href="http://sobereva.com/712" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=http://sobereva.com/712&source=gmail&ust=1738270636915000&usg=AOvVaw1F72pQs6o0ZIj2AWvsXK-R">http://sobereva.com/712</a> (currently written in Chinese, please use Google translator. English version may be available later). Please feel free to contact me if you have any problem.</div><div><br></div><div>Best regards,</div><div><br></div><div>Tian Lu<br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, January 29, 2025 at 8:05:13 AM UTC+2 Alexander Davis 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"><div>Hi,</div><div><br></div><div>I have been calculating self-consistent Hirshfeld charges using my own post-processing script, but I am thinking of switching to the CP2K implementation. However, there are two questions I would need answered.</div><div><br></div><div>One question is how it handles a well-known issue with this method, which is that it obtains the shape function from the electron densities of free ions, but sometimes these ions don't really exist. This is why my implementation is unsatisfactory: I have references for N- and N2-, but the electrons are not really attached, so the reference densities are too wide. Every self-consistent Hirshfeld implementation as to handle this somehow, I would like to know how CP2K handled it. The documentation says "This scales only the full shape function, not the added charge as in the original scheme", which is suggestive of a solution, but where can I read more about it? If the details are in one of the publications, I missed it.</div><div><br></div><div>My second question is whether I save the reference atoms as cube files, or even better, the actual weight function used for integration (which is a(r)/(a(r)+b(r)), where a and b are the reference atom densities). It's a quirk of my application, I need precisely the weight function that was used to calculate the charges.</div><div><br></div><div>Thanks,</div><div>Alex<br></div></blockquote></div></blockquote></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 visit <a href="https://groups.google.com/d/msgid/cp2k/1fc65746-417c-4ff6-8073-766b86561a9en%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/1fc65746-417c-4ff6-8073-766b86561a9en%40googlegroups.com</a>.<br />