<div dir="ltr"><div>I think I found the CP2K setting which is causing the problem: When I remove the DFTB3_PARAM settings, then my water molecules are not falling apart anymore. </div><div><br></div><div>I will open a new topic because this issue not anymore about periodic QM/MM. </div><div><br></div><br><br>On Monday, September 11, 2017 at 5:01:15 PM UTC+2, Jadzia wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi Yingchun, <div><br></div><div>thank you so much for your answer, I'm really glad about any comments.</div><div><br></div><div>Before I switched to QM/MM I have run a lot MM simulations of the same system, and they run all fine.</div><div><br></div><div>The first E_COUPL which I wanted to try was GAUSS as you had suggest, however CP2K rejected my request telling me E_COUPL=GAUSS is not implemented for DFTB.</div><div><br></div><div>As you seemed to indicate, I also guess that the electron spillout problem might be the cause. Closely related I could also imagine that the problem might be caused by some other electrostatic related issue. </div><div><br></div><div>What would you try next to find the cause or a solution?</div><div><br></div><div>Many thanks,</div><div>Jadzia</div><div><br></div><div><br></div><div><br><br>On Monday, September 11, 2017 at 4:20:10 PM UTC+2, Yingchun Zhang wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi, Jadzia<div>  just some suggestion:</div><div>  </div><div>  1, If you want to avoid the spill out problem as you indicated in you original post, you might want to use E_COUPL = GAUSS;</div><div>  </div><div>  2, check if your pure MM run of the whole systems works.</div><div><br></div><div>Yingchun</div></div></blockquote></div></div></blockquote></div>