<div>Hello, </div><div>I just want to update people who may have the same question in the future. After further research, my problems could be partially solved with the following changes to my input script: </div><div><br /></div><div>...</div><div>CHARGE +1<br />    LSD<br />    MULTIPLICITY 1</div><div>...</div><div>POISSON_SOLVER</div><div>PERIODIC NONE</div><div>...</div><div>&CENTER_COORDINATES<br />        &END<br />...</div><div><br /></div><div>All of this, and a sufficiently large cutoff (my research shows 5+ A) worked for me for case a) where I have a multiplicity of 1. Centering is there to prevent problems at the edge of the box with the solver.</div><div>Sadly for system b) where the multiplicity is higher, I always get the usual wavelet solver issue, even though I </div><div><br /></div><div>a) Center my atoms</div><div>b) Have a large box (more than 10 A from edges even)</div><div>c) have a large cutoff</div><div><br /></div><div>I keep getting problems with the Wavelet solver at the edge: Usually, the error message</div><div><br /></div><div> *** WARNING in ps_wavelet_methods.F:238 :: Density non-zero on the edges ***<br /> *** of the unit cell: wrong results in WAVELET solver                    ***</div><div><br /></div><div>can be treated with a larger box, it seems to me from looking around the Cp2kgroup that higher multiplicity might be difficult, but even when increasing the box to have <b>20+ A vacuum</b> in each direction I still get a problem at the edge. I attach my input script for this issue, maybe somebody can spot some obvious mistake quickly! Is this maybe a chemistry issue (getting multiplicity wrong? I just cut out a molecule off a larger structure, so at each branch of the Triazine, a Hydrogen is missing)</div><div><br /></div><div>Thank you for your patience</div><div class="gmail_quote"><div dir="auto" class="gmail_attr">Andreas Döll schrieb am Mittwoch, 4. Februar 2026 um 11:46:43 UTC+1:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Dear cp2k community,</div><div><br>I have two charged systems I want to perform a single energy calculation.</div><div><br></div><div>System a) Neutral molecule with a Lithium at some distance</div><div>System b) Charged molecule (cut out of a larger molecular structure, so 3 Hydrogens are missing) </div><div><br></div><div>For both calculations, my energy either diverges or makes jumps. <br>Usually I have performed calculations with neutral molecules, not charged systems, so I am not familiar with the proper procedure. I followed the example give here <a href="https://www.cp2k.org/exercises:2014_ethz_mmm:hfx_h2ion" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=de&q=https://www.cp2k.org/exercises:2014_ethz_mmm:hfx_h2ion&source=gmail&ust=1770392965900000&usg=AOvVaw0A0LMJivqLWiEY3Eh6xwxr">https://www.cp2k.org/exercises:2014_ethz_mmm:hfx_h2ion</a> </div><div>I defined my Charge in the DFT section, as well as using LSD option and switched to a Nonperiodic system - as far as I understand this is necessary? Do I need to define the multiplicity of my system as well ?<br><br>Cheers, </div><div>Andreas</div><div><br></div><div><br></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 visit <a href="https://groups.google.com/d/msgid/cp2k/fc486a4b-bc1d-46be-a420-43295454d310n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/fc486a4b-bc1d-46be-a420-43295454d310n%40googlegroups.com</a>.<br />