Hello CP2K Users,<div><br></div><div>I got the answer to this question courtesy of one of the authors of the above-mentioned article (Dr. Kristof Bal). I am sharing it here</div><div><br></div><div><p>The charge on the counterion can
be controlled by not giving it any basis set. So, make it a "new"
atom type that has no basis assigned to it. Here's an example: </p>

<p>   
&KIND D</p>

<p>   
  ELEMENT H</p>

<p>   
  BASIS_SET NONE</p>

<p>   
  POTENTIAL GTH-PBE-q1</p>

<p>   
&END KIND</p>

<p> This gives you a
+1 ion, provided such a particle is placed in a vacuum. CP2K is an AO code, so it
can only place electrons around cores that have basis functions around them.
This is very different from PW codes like VASP. </p>

<p>If you request a neutral system
from CP2K, and include one such "D" atom in the system, you will have
a -1 slab. After all, the one "hydrogen" electron cannot be placed
around D, and must end up in the slab. </p>

<p>You can get the opposite charge
distribution by using the CORE_CORRECTION feature in CP2K to change the charge
of the core of "D". But be aware that charge bookkeeping becomes
tricky. For a +1 slab with -1 counterion, you must set CORE_CORRECTION to -2
(lower from +1 to -1) but CHARGE to +2. Why? Well, you want "one
less" electron in the slab, but also not add the extra electron you would
introduce by bringing in an extra "H". So, "not add" two
electron => CHARGE +2.</p><p>Regards,<br></p><p>Aditya</p><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Thursday, March 3, 2022 at 9:24:56 AM UTC-5 adl233 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>I apologize in advance if this is a trivial question. I am trying to model a system where a slab surface is charged in a non-periodic cell and a counterion is added to make the entire system charge neutral. I am trying to follow the procedure from the following paper: <a href="https://iopscience.iop.org/article/10.1088/1361-6595/aaa868/meta" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://iopscience.iop.org/article/10.1088/1361-6595/aaa868/meta&source=gmail&ust=1646421258528000&usg=AFQjCNE9l7xrgAyyZld1lyvMVrCmec-8jw">https://iopscience.iop.org/article/10.1088/1361-6595/aaa868/meta</a></div><div><br></div><div>However, I do not understand the exact method by which this was implemented. I read about core-correction option as well as CDFT and confused about the right way to proceed. I am copying parts of the paper relevant to this question below.</div><div><br></div><div>"The practical realization of this approach in a standard DFT code (CP2K) is as follows. The negative surface charges in this work require a positive countercharge which, in the simplest case, can be a proton. It is, however, not always straightforwardly possible to introduce gas phase ions of a specified charge into the simulation box. This can be compensated by using a DFT code that expands the Kohn–Sham orbitals in an atom-centered (localized) basis: if no basis functions are added on the counterion, no electronic density can spill over, its charge can be precisely controlled and the desired surface charge can be enforced. The method is in principle readily usable in any DFT code that uses localized basis sets (such as the here used CP2K) but has, to the best of our knowledge, not yet been described in the literature."</div><div><br></div><div>and </div><div><br></div><div>"In the setup adopted in this work, a single additional electron is considered and the countercharge (a proton) is placed at a Z position of 40 Å in a box of dimensions 16.1606 × 16.8106 × 100 Å3"</div><div><br></div><div>Regards,</div><div>Aditya</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/d568db0f-cba1-435e-a443-a0836a72ced1n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/d568db0f-cba1-435e-a443-a0836a72ced1n%40googlegroups.com</a>.<br />