<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">please have a look at the section : <a href="https://manual.cp2k.org/cp2k-3_0-branch/CP2K_INPUT/FORCE_EVAL/EXTERNAL_POTENTIAL.html" class="">https://manual.cp2k.org/cp2k-3_0-branch/CP2K_INPUT/FORCE_EVAL/EXTERNAL_POTENTIAL.html</a><div class=""><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On 8 Nov 2017, at 08:51, Junbo Lu <<a href="mailto:luju...@gmail.com" class="">luju...@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class="">Dear Teodoro,</div><div class=""> Thank you for your good suggestion. In option 3, you mean I can apply a physical restraint. Then how can I apply this spherical wall? By which keyword in input file?</div><div class=""> Best regards,</div><div class=""> Junbo </div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">2017-11-08 15:22 GMT+08:00 Teodoro Laino <span dir="ltr" class=""><<a href="mailto:teodor...@gmail.com" target="_blank" class="">teodor...@gmail.com</a>></span>:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word" class="">Dear Junbo,<div class=""><br class=""></div><div class="">first of all many thanks for having posted your request here: more people will be able to benefit of the on-going discussion.</div><div class=""><br class=""></div><div class="">Let me explain my point of view regarding your problem, with the understanding that they are very much in-line with Juerg’s comment.</div><div class=""><br class=""></div><div class="">First of all as you already realised, the ADBF QM/MM needs to have a constant charge and the charge is not allowed to change during the run. There is a <u class="">very good reason</u> for that: if you change the charge of the QM region you are changing the zero reference for your calculations during MD and this makes the analysis of your simulations less trivial.</div><div class=""><br class=""></div><div class="">Therefore, this is not about an "unimplemented feature” (*ignore the recommandation received in a separate list, to implement this functionality*), rather a protection to prevent users from running non-sense simulations.</div><div class=""><br class=""></div><div class="">Once we agree that we need to use the ADBF system with this “limitation", let’s move on to your questions.</div><div class=""><br class=""></div><div class="">(1) QM/MM runs do not have to be necessarily neutral. You can run QM/MM of QM charged systems (both fully periodic and non-periodic) being extremely careful to how you are decoupling the periodic images (in case of fully non-periodic calculations). You may have the corresponding counterion in the classical (MM) portion and the interaction would still be computed correctly (again be careful to decouple the QM part properly), in case QM+MM is neutral. You could also have the QM+MM charged and this would still work ok.</div><div class=""><br class=""></div><div class=""><b class="">Option 1</b>: I<u class="">f you believe that the Cl- is not important for your simulation you can freely decide to remove it from the entire system and run instead a simulation that is overall (QM+MM) charged (2+).</u></div><div class=""><br class=""></div><div class=""><b class="">Option 2</b>: <u class="">If you decide to treat the Cl- classically, be sure that the Cl- is always far from your Zn2+ center enough (at least 2 solvation shells) to avoid any interaction between the two ions that could not be treated classically (electrostatic+non-bonded).</u> In this case your overall system (QM+MM) is neutral, with your QM region being 2+ and the MM one 2-. Please, use a decoupling to be sure that the strong electrostatic coupling between the QM periodic replicas is dumped. </div><div class=""><br class=""></div><div class="">(2) Using ADBF QM/MM for both option 1 and 2 is trivial as your QM charge is not changing (and again, for fundamental reasons it should not!).</div><div class=""><b class="">Option 3</b>: If you start with Option 2, you may notice instead that Cl- may play a fundamental role in your simulations or you may already know before starting your simulations that you are interested in mapping the interaction between the Cl- and Zn2+.</div><div class="">In this case, there is only one thing that you can do.<u class=""> If you believe that the Cl- is an important part of the QM center, then you’ll have to include the Cl- into the QM region since the very beginning of your simulations.</u> This will have the major drawback that Cl- may leave the Zn2+ center with the drawback of being forced to use very large QM boxes. </div><div class="">Still, there is nothing lost or wasted in life (and in simulations too): you can apply a physical restraint (i.e. spherical wall) to prevent the Cl- from flying far and subsequently reweight the calculated free energy (assuming you are interested in free energies) for the additional constraint. This will help you in containing the size of your QM box even for extended simulation times (several hundreds of ps).</div><div class=""><br class=""></div><div class="">Personally, I like option 3.. but depending on your understanding of the system you may decide (equally good) to go for option 2 or option 1. … Possibly even for another option not necessarily considered here….</div><div class=""><br class=""></div><div class="">Hope this helps and good luck!</div><div class="">Best,</div><div class="">Teo</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""><div class=""><blockquote type="cite" class=""><span class=""><div class="">On 4 Nov 2017, at 17:07, <a href="mailto:luju...@gmail.com" target="_blank" class="">luju...@gmail.com</a> wrote:</div><br class="m_8090505173029907274Apple-interchange-newline"></span><div class=""><span class=""><div dir="ltr" class=""><div style="max-height:10000px" class=""><div dir="ltr" class="">Hi, I want to do a simulation of hydration of Zn2+ by QM/MM MD of CP2K. But there are several problems for my simulation:<br class="">(1). I do not know whether we should neutralize the system by minus ion in CP2K, such as Cl-. We can also not neutralize the system, is it right?<br class="">(2). If we add some minus ions into system, how can we define the charge of dynamic QM region in adbf QM/MM MD calculation? When a charged molecule enter the QM region, such as Cl- ion, the total charge of dynamic QM region will change. I have read the CP2K manual, I find it can not change the charge of QM region when doing MD. So this is a problem.<br class="">Can you help me for these two problems? I will be very appreciated!<br class="">With my best regards,<br class="">Junbo<br class=""></div></div><a class="m_8090505173029907274gwt-Anchor" style="display:none">显示删减掉的内容</a> <div class="m_8090505173029907274F0XO1GC-nb-P" style="display:none"><div class=""></div></div><div class=""></div><div class=""></div><div style="display:none" class=""></div><div style="display:none" class=""></div><div class=""><div class="m_8090505173029907274F0XO1GC-ed-a"></div></div><div class="m_8090505173029907274F0XO1GC-nb-b"><div class="m_8090505173029907274F0XO1GC-nb-a m_8090505173029907274F0XO1GC-nb-cb"><div class=""><div style="display:inline-block" class=""><div style="display:none" class=""></div></div> </div></div></div></div><div class=""><br class="m_8090505173029907274webkit-block-placeholder"></div>
-- <br class=""></span>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br class="">
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" target="_blank" class="">cp2k+unsubscribe@googlegroups.<wbr class="">com</a>.<span class=""><br class="">
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank" class="">cp...@googlegroups.com</a>.<br class="">
Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" class="">https://groups.google.com/<wbr class="">group/cp2k</a>.<br class="">
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" class="">https://groups.google.com/d/<wbr class="">optout</a>.<br class="">
</span></div></blockquote></div><br class=""></div></div><div class="HOEnZb"><div class="h5"><div class=""><br class="webkit-block-placeholder"></div>
-- <br class="">
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.<br class="">
To unsubscribe from this topic, visit <a href="https://groups.google.com/d/topic/cp2k/0YL4lM9Rqj8/unsubscribe" target="_blank" class="">https://groups.google.com/d/<wbr class="">topic/cp2k/0YL4lM9Rqj8/<wbr class="">unsubscribe</a>.<br class="">
To unsubscribe from this group and all its topics, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" target="_blank" class="">cp2k+unsubscribe@googlegroups.<wbr class="">com</a>.<br class="">
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank" class="">cp...@googlegroups.com</a>.<br class="">
Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" class="">https://groups.google.com/<wbr class="">group/cp2k</a>.<br class="">
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" class="">https://groups.google.com/d/<wbr class="">optout</a>.<br class="">
</div></div></blockquote></div><br class=""></div><div class=""><br class="webkit-block-placeholder"></div>
-- <br class="">
You received this message because you are subscribed to the Google Groups "cp2k" group.<br class="">
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" class="">cp2k+uns...@googlegroups.com</a>.<br class="">
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" class="">cp...@googlegroups.com</a>.<br class="">
Visit this group at <a href="https://groups.google.com/group/cp2k" class="">https://groups.google.com/group/cp2k</a>.<br class="">
For more options, visit <a href="https://groups.google.com/d/optout" class="">https://groups.google.com/d/optout</a>.<br class="">
</div></blockquote></div><br class=""></div></div></body></html>