<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: 0 0 0 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" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=http://input.xyz&source=gmail&ust=1680856345197000&usg=AOvVaw2PkTEOWRHNoVnY4DFZIW3B">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" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://www.cp2k.org/howto:resp&source=gmail&ust=1680856345197000&usg=AOvVaw1aNME-Vec1LVUQy1wXnVjw">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 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/b159b64d-d9ce-434c-b445-10fa4d536d81n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/b159b64d-d9ce-434c-b445-10fa4d536d81n%40googlegroups.com</a>.<br />