<div dir="ltr">The code is<div><br></div><div><pre id="hterm:copy-to-clipboard-source" style="color: rgb(0, 0, 0); user-select: text;"> IF (dft_control%qs_control%gapw) THEN
CALL pw_transfer(rho0_s_gs%pw, wf_r%pw)
IF (dft_control%qs_control%gapw_control%nopaw_as_gpw) THEN
CALL pw_axpy(rho_core%pw, wf_r%pw)
END IF
ELSE
CALL pw_transfer(rho_core%pw, wf_r%pw)
END IF</pre><div><br></div><div>So for GPW it will add in a Gaussian charge with value +Zeffective centred on the atoms the same as is done when calculating electrostatics (see original paper).</div><div><br></div><div>It is quite possible to post process the results of a trajectory to add in an atom centred density.</div><div><br></div><div>Matt</div><br>On Tuesday, May 1, 2018 at 8:54:09 PM UTC+1, Ada Sedova wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Right- I know you said that. The response from David didn't seem to relate to the original question. You also said you hadn't tried it. The link to the cp2k forum discussion I added above shows that it does not in fact add it, if you are not using GAPW.<div><br></div><div>So, the question is, does it add it or not? It seems from the above, it does not. Someone familiar with the code would have to determine if this is true or not. For me, I have tens of picoseconds of simulation for systems with hundred of atoms, and re-running this to check would be a pain. I would do it if I knew for sure it would work. I'm trying to monitor the changing densities.</div><div><br></div><div> </div></div><div><br><div class="gmail_quote">On Tue, May 1, 2018 at 1:21 PM, Matt W <span dir="ltr"><<a href="javascript:" target="_blank" gdf-obfuscated-mailto="feg8P9L8AQAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">matt...@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><span><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>According to Matt, just running with the TOTAL...etc command will add the core electrons. Are you saying this works? Because a previous thread, linked above, describes this as not working unless you are using GAPW, not just GPW.</div><div><br></div></div></blockquote><div><br></div></span><div>Hi Ada,</div><div><br></div><div>no, I said it would a add in a core charge at the nuclei equal to the Z effective of the pseudo at the atom centre. </div><div><br></div><div>What is required for getting a correct topological description to get AIM descriptors? I haven't read up on it.</div><div><br></div><div>Adding in an atomic density that does not respond to the valence density could be made into a fairly simple post-processing step - is that what is needed, the [core] atom density (from an isolated atom) added to the density cube file?</div><div><br></div><div>Matt</div><div> </div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div><div>-Ada</div><div><br></div><div><br></div></div><div><div><div><br><div class="gmail_quote">On Tue, May 1, 2018 at 9:33 AM, David T <span dir="ltr"><<a rel="nofollow">amazi...@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">You are right Matt.<div>I did it for the 1st time in this article: <a href="https://pubs.acs.org/doi/abs/10.1021/jp901753p" rel="nofollow" target="_blank" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fpubs.acs.org%2Fdoi%2Fabs%2F10.1021%2Fjp901753p\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHwM8aEifhuI-KdfkRVmp3LMvfbjA';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fpubs.acs.org%2Fdoi%2Fabs%2F10.1021%2Fjp901753p\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHwM8aEifhuI-KdfkRVmp3LMvfbjA';return true;">https://pubs.acs.org/doi/abs/<wbr>10.1021/jp901753p</a></div><div><br></div><div>Sadly Henkelman and VASP persons have started using my Algorithm without citing it... but that is another story....</div><div><div><div><br><br>On Tuesday, 17 April 2018 19:19:13 UTC+1, 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">I _think_ it should add the approximate core charge for GPW (i.e. a gaussian with charge +Z centred on the atom), which should be enough for the correct topology for AIM. <br><br>I'm not speaking from experience though. May well be wrong.<br><br>On Tuesday, April 17, 2018 at 7:08:09 PM UTC+1, Ada Sedova wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">I think that only works if you're using GAPW, though, am I right? Because if you are not calculating the core electrons they won't print? At least that's according to another thread I read on this forum...</div><div><br><div class="gmail_quote">On Tue, Apr 17, 2018 at 12:49 PM, Matt W <span dir="ltr"><<a rel="nofollow">matt...@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><a href="https://manual.cp2k.org/cp2k-5_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/TOT_DENSITY_CUBE.html" rel="nofollow" target="_blank" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmanual.cp2k.org%2Fcp2k-5_1-branch%2FCP2K_INPUT%2FFORCE_EVAL%2FDFT%2FPRINT%2FTOT_DENSITY_CUBE.html\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHHuaJv-ldLNwzfRe2MZ8pEP0ed-g';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmanual.cp2k.org%2Fcp2k-5_1-branch%2FCP2K_INPUT%2FFORCE_EVAL%2FDFT%2FPRINT%2FTOT_DENSITY_CUBE.html\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHHuaJv-ldLNwzfRe2MZ8pEP0ed-g';return true;">https://manual.cp2k.org/cp2k-<wbr>5_1-branch/CP2K_INPUT/FORCE_<wbr>EVAL/DFT/PRINT/TOT_DENSITY_<wbr>CUBE.html</a><div><br></div><div>I think.<br><br>On Tuesday, April 17, 2018 at 5:14:12 PM UTC+1, Ada Sedova wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi all-<div>There was a thread from 2016 about finding a way to add core electrons to a cp2k cube file so that Bader analysis can be performed with Henkelman's code, which takes cube files, but of the Gaussian format, which includes all electrons. There was no definitive solution posted, but a suggestion was to hack the Manz DDEC code to add core electrons and print out a new cube. The Manz code has gone through some changes and a brief look at it did not reveal a rewrite option. </div><div><br></div><div>I was wondering if there have been any updates on the ability to use Henkelman's with a cp2k cube file. The other option would be to create a VASP CHGCAR file from a cp2k cube file, but the core electrons still need to be added somehow.</div><div><br></div><div>Thanks!</div><div>AS</div></div></blockquote></div></div><span><font color="#888888">
<p></p>
-- <br>
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br>
To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/IVj1WAxCAYs/unsubscribe" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/d/topic/cp2k/IVj1WAxCAYs/unsubscribe';return true;" onclick="this.href='https://groups.google.com/d/topic/cp2k/IVj1WAxCAYs/unsubscribe';return true;">https://groups.google.com/d/<wbr>topic/cp2k/IVj1WAxCAYs/<wbr>unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, send an email to <a rel="nofollow">cp2k+...@googlegroups.com</a>.<br>
To post to this group, send email to <a rel="nofollow">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.<br>
</font></span></blockquote></div><br></div>
</blockquote></div></blockquote></div></div></div></div><div><div>
<p></p>
-- <br>
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br>
To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/IVj1WAxCAYs/unsubscribe" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/d/topic/cp2k/IVj1WAxCAYs/unsubscribe';return true;" onclick="this.href='https://groups.google.com/d/topic/cp2k/IVj1WAxCAYs/unsubscribe';return true;">https://groups.google.com/d/<wbr>topic/cp2k/IVj1WAxCAYs/<wbr>unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, send an email to <a rel="nofollow">cp2k+...@googlegroups.com</a>.<br>
To post to this group, send email to <a rel="nofollow">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.<br>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><div><div>
<p></p>
-- <br>
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br>
To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/IVj1WAxCAYs/unsubscribe" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/topic/cp2k/IVj1WAxCAYs/unsubscribe';return true;" onclick="this.href='https://groups.google.com/d/topic/cp2k/IVj1WAxCAYs/unsubscribe';return true;">https://groups.google.com/d/<wbr>topic/cp2k/IVj1WAxCAYs/<wbr>unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, send an email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="feg8P9L8AQAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.<wbr>com</a>.<br>
To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="feg8P9L8AQAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.<br>
</div></div></blockquote></div><br></div>
</blockquote></div></div>