Thanks, just want to check if this is right way to do this:<div>I defined atoms labeled <b>X</b> in the coord file and then I added this section into the SUBSYS:</div><div><br>    &KIND X<br>      ELEMENT H<br>      GHOST .TRUE.<br>      BASIS_SET NONE<br>      POTENTIAL GTH-PBE-q1<br>      CORE_CORRECTION 0.12   # 3/25<br>    &END KIND</div><div><br></div><div>Is this the right way to divide for example 3<i>e</i> charge over 25 fake hydrogen atoms?</div><div><br></div><div>Branislav<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, February 21, 2022 at 4:31:49 PM UTC+1 mattwa...@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;">Multiple options - most directly:<div><a href="https://manual.cp2k.org/cp2k-9_1-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND.html#list_CORE_CORRECTION" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://manual.cp2k.org/cp2k-9_1-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND.html%23list_CORE_CORRECTION&source=gmail&ust=1645704178459000&usg=AFQjCNEKATW85_d9aDf9i6dpGdvQCJ0_Lw">https://manual.cp2k.org/cp2k-9_1-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND.html#list_CORE_CORRECTION</a></div><div>lets you change the core charge on a pseudopotential. If you make a fake hydrogen and adjust charge appropriately you'd get close to what you were after (won't be a point charge but gaussian smeared charge).<br><br></div><div>There are other routines that could load an external potential on a grid or make use of the QM/MM code in some form.</div><div>Matt</div><div><br></div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, 21 February 2022 at 15:24:34 UTC <a href data-email-masked rel="nofollow">dr.bra...@gmail.com</a> 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 cp2k users,<div><br></div><div>Is it possible to introduce ghost atoms with no basis set and with finite non-integer point charge in AIMD calculations? Something like MM point charges.</div><div><br></div><div>I want to try to simulate behaviour of some small molecules in the electric field which originates from a particular sequence of the point charges.</div><div><br></div><div>Thanks</div><div>Branislav</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 on the web visit <a href="https://groups.google.com/d/msgid/cp2k/4cd597ae-c3b1-4dd2-b061-2f26075bf057n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/4cd597ae-c3b1-4dd2-b061-2f26075bf057n%40googlegroups.com</a>.<br />