<div dir="ltr">Hello Marcella,<div><br></div><div>Thank you for the helpful information! On that note, I would like to ask 2 final questions:</div><div>-is there a way to perform GCE-DFT (Grand Canonical Ensemble) in CP2K? I am asking because my goal here is to include the potential of the electrode (like in experiments). I read the following work, <a href="https://www.sciencedirect.com/science/article/pii/S2451910321000636">https://www.sciencedirect.com/science/article/pii/S2451910321000636</a>, and I wanted to know if something similar is implemented in CP2K.</div><div>- I am also trying to include the solvent effects of the electrolyte with the implicit SCCS model (Andreussi) that is available on CP2K. However, it does not work with my stress tensor tag, which I had to remove. I also saw that in previous threads here, SCCS does not work with kpoints? Any insight on this matter would be immensely helpful.</div><div><br></div><div>Thank you so much for your time.</div><div><br></div><div>Hasan.</div><div><br></div><div><br></div><div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Hasan Al-Mahayni</div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Apr 6, 2023 at 4:44 AM Marcella Iannuzzi <<a href="mailto:marci.akira@gmail.com">marci.akira@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><br></div><div><br></div>Dear Hasan <div><br></div><div>The RESP feature is meant to assign charges that reproduce the DFT electrostatic potential as post processing.</div><div>The electronic density is not modified and the SCF calculation is not affected.</div><div><br></div><div>DFT in CP2K works only with a total integer charge. The constant Fermi energy method cannot be applied.</div><div><br></div><div>Regards</div><div>Marcella</div><div><br></div><div><br></div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, April 5, 2023 at 9:12:36 PM UTC+2 Hasan Al-Mahayni wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello everyone,<div><br></div><div>I have a simple (100) copper slab (100 atoms, 4 layer==> see <a href="http://input.xyz" rel="nofollow" target="_blank">input.xyz</a>); and I am trying to make the overall charge of my system equal to +1.39. Since the CHARGE tag only takes integers, I am trying to give partial charges to the atoms in the top layer to go from charge +1 to charge +1.39. I used a cp2k tutorial as reference: <a href="https://www.cp2k.org/howto:resp" rel="nofollow" target="_blank">https://www.cp2k.org/howto:resp</a> </div><div>I attach my input file ("cp2k.inp") and my output file "RESULTS". The simulation initially runs fine, the first SCF loop converges, then the simulation crashes and I get this error (See "error.txt"):</div><div> </div><div><i><b>srun: error: cdr1443: tasks 3,5,7-8,14-15,20,23-25,27,30: Segmentation fault (core dumped)</b></i><br></div><div><i><b><br></b></i></div><div>This is the tags I added:</div><div><b><i>&PROPERTIES<br> &RESP<br> &RESTRAINT<br> ATOM_LIST 76..100<br> TARGET 0.015477832<br> STRENGTH 0.0001<br> &END RESTRAINT<br> &END RESP<br>&END PROPERTIES</i></b><br></div><div><b><i><br></i></b></div><div>Does anyone know why I am getting this error? Or how can I include partial charges correctly in my system? The goal here is to have a certain Fermi Energy for my system, so if there is a better way to do this, I would love to know that.</div><div><br></div><div>Thanks for your time!</div><div><br></div><div>Hasan.</div></blockquote></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/08xyElGa3Jg/unsubscribe" target="_blank">https://groups.google.com/d/topic/cp2k/08xyElGa3Jg/unsubscribe</a>.<br>
To unsubscribe from this group and all its topics, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/b159b64d-d9ce-434c-b445-10fa4d536d81n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/b159b64d-d9ce-434c-b445-10fa4d536d81n%40googlegroups.com</a>.<br>
</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/CAECRhQzuyH6VgegfMFJ0u7-gYzF-Ak4%3DykzAcU56LZZzNKCb_w%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAECRhQzuyH6VgegfMFJ0u7-gYzF-Ak4%3DykzAcU56LZZzNKCb_w%40mail.gmail.com</a>.<br />