There will be the core charge of the original atom as well - so you need neutralise the atom core (-1 in your case) then add your charge back.<div>Yes - the charge keyword will alter the number of electrons in the system so if you want these to be virtual protons not hydrogen atoms then you do need to remove the electron as well with the charge keyword.</div><div>If your overall system is charged then normal issues of charged periodic systems apply.. </div><div>Matt<br><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, 28 February 2022 at 12:37:58 UTC dr.bra...@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;">Additionally, should I in such cases change total charge od the system?<br><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, February 23, 2022 at 2:09:59 PM UTC+1 Бранислав Миловановић 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">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 <a rel="nofollow">mattwa...@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">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" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://manual.cp2k.org/cp2k-9_1-branch/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND.html%23list_CORE_CORRECTION&source=gmail&ust=1646323265347000&usg=AFQjCNHgJWz8MahvOqEZT54-o6dCE48caQ">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 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></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/21c1483b-e7ae-4491-bd49-c98de6790796n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/21c1483b-e7ae-4491-bd49-c98de6790796n%40googlegroups.com</a>.<br />