<html><head><meta http-equiv="Content-Type" content="text/html charset=iso-8859-1"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Any Ewald scheme (for quantum calculation this is also known as Poisson Solver) for periodical charged systems, includes automatically the compensating background charge. </div><div>The energy of a periodical charged system, without background charge, is undefined (the series 1/r is not converging for systems with a net charge different from zero - actually, in general, this series is only conditionally convergent but not absolutely convergent).</div><div><br></div><div>I don't really understand your claim (and didn't check the paper you mentioned) - but there is nothing special in cp2k compared to other codes doing Ewald (cp2k, cpmd, quantum espresso, etc..).. So the answer is:</div><div><br></div><div> there is no section to activate to have the compensating background charge to have an overall neutral simulation box : simply run a standard setup for a 3D periodic system (without decoupling) and you will get automatically this term in the energy.</div><div><br></div><div><br></div><div><br></div><div><div><div><div>On Aug 8, 2013, at 10:29 AM, Gianni Carifaccia <<a href="mailto:gian...@gmail.com">gian...@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr"><div>Dear All,</div><div><br></div><div>I need to simulate a charged QM system (I need to see what happens after protonation of my molecule) and I cannot make the whole box neutral by adding a counterion. Looking into literature, I've seen that it is possible to solve the charged box issue in cp2k through Ewald, by smearing a compensating background charge to have an overall neutral simulation box (one for all, see Sprik et al JCTC 2011, 7, 1951. <a href="http://pubs.acs.org/doi/pdf/10.1021/ct100715">http://pubs.acs.org/doi/pdf/10.1021/ct100715</a>).</div>
<div><br></div><div>My question is: which exact input section should I look for ? I looked around the Poisson solvers, but I don't feel confident I found the right place. Other Ewald is in the MM section, but I don't think I've to run a QM/MM if a I have a pure QM box. In case it is in the Poisson section, which keywords ? Further, I've not seen how to input the amount of the background charge. Is this "understood" automatically by cp2k ?</div>
<div><br></div><div>Of course, any help on how to handle this issue, and how to check things are done properly, is wellcome.</div><div><br></div><div>Thank you in advance,</div><div>Giancar</div><div><br></div></div><div><br class="webkit-block-placeholder"></div>
-- <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+uns...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/groups/opt_out">https://groups.google.com/groups/opt_out</a>.<br>
<br>
<br>
</blockquote></div><br></div></div></body></html>