<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">here</a> which was intended to get the same simulations running with fully periodic QM/MM PBC). As I have found it, 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">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">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 the 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>