<div dir="ltr">I have contacted one of the developers of DFTB3. He has also much experience with DFTB3 QM/MM water simulations, and he said that even though DFTB3 has a few problems with water, these should not cause them to fall apart. <div><br></div><div>Therefore there seem to be two possibilities:</div><div><ol><li>Something is wrong with the setup. <br></li><li>There is a bug in CP2K. </li></ol><div>I have tried my best to check if something is wrong with the setup, but I couldn't find anything. Also, I have run a lot of simulations with different settings. </div><div>If it is a bug in CP2K I cannot tell, this might be a question for the developers. </div><div><br></div><div>Best wishes,</div><div>Jadzia</div><div><br><br>On Wednesday, September 13, 2017 at 10:38:16 AM 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 everyone,<div><br></div><div>I'm trying to use DFTB3 as the QM method of a QM/MM simulation of water (which is only a test system).</div><div><br></div><div>However, my QM water molecules are falling apart. At first some H20 molecules form H3O+ and OH-, and then it goes on to form H4O++ and O--. All this happens between QM atoms only. </div><div><br></div><div>I have adjusted and switched many QM and QMMM parameters (see also my last topic <a href="https://groups.google.com/forum/#!topic/cp2k/EmRoMadjbRw" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/forum/#!topic/cp2k/EmRoMadjbRw';return true;" onclick="this.href='https://groups.google.com/forum/#!topic/cp2k/EmRoMadjbRw';return true;">here</a> which was intended to get the same simulations running with fully periodic QM/MM PBC). As I have found out, it seems that the Hubbard derivative (which is to my knowledge the DFTB3_PARAM in CP2K) is related to the problem. If I do not specify this parameter for the H and O atoms (it is set to 0 by default by CP2K) then my waters are not falling apart. The Hubbard derivatives I have used are the ones provided by <a href="http://dftb.org" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fdftb.org\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFXugspV84tTv_WctI3MTXftCJGyw';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fdftb.org\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFXugspV84tTv_WctI3MTXftCJGyw';return true;">dftb.org</a>. </div><div><br></div><div>Do you think it is fine to run the simulations without specifying the Hubbard derivatives? (i.e. using the default value of 0) </div><div>And do you have any other suggestions or comments on why this problem could arise or solved?</div><div><br></div><div>If I run QM-only simulations with the same DFTB3 settings and with the official Hubbard derivative parameters from <a href="http://dftb.org" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fdftb.org\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFXugspV84tTv_WctI3MTXftCJGyw';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fdftb.org\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFXugspV84tTv_WctI3MTXftCJGyw';return true;">dftb.org</a>, then the simulations run just fine. Only when enabling QM/MM my QM water falls apart.</div><div><br></div><div>However, I now found out that even without Hubbard derivatives I manage to let the water fall apart, namely when I change the Poisson settings of the QM region. With the following QM Poisson settings the water does not fall apart:</div><div><div><br></div><div><font size="1">    &POISSON</font></div><div><font size="1">      &EWALD</font></div><div><font size="1">       EWALD_TYPE SPME</font></div><div><font size="1">       GMAX 25</font></div><div><font size="1">       O_SPLINE 5</font></div><div><font size="1">      &END EWALD</font></div><div><font size="1">    &END POISSON</font></div></div><div><br></div><div>The entire input file is attached as well. When changing the Poisson settings to values which seem to be more optimal:</div><div><div><br></div><div><font size="1">    &POISSON</font></div><div><font size="1">      &EWALD</font></div><div><font size="1">       ALPHA .44</font></div><div><font size="1">       EWALD_TYPE SPME</font></div><div><font size="1">       GMAX 42 42 43</font></div><div><font size="1">       O_SPLINE 6</font></div><div><font size="1">      &END EWALD</font></div><div><font size="1">    &END POISSON</font></div></div><div><br></div><div>then the QM water falls apart. Which I find strange since these Poisson settings are more "sound" because they are more close to what is recommended. E.g. the box size is 42 42 53 Angstrom, and it is recommended to use 1 per Angstrom for GMAX. Also, O_SPLINE was set before to 5, which can lead to an interpolation problem to my knowledge. </div><div><br></div><div>Looking forward to hearing from you and your recommendations. </div><div><br></div><div>Best wishes,</div><div>Jadzia</div><div><br></div><div><br></div><div><br></div></div></blockquote></div></div></div>