<div dir="ltr">Dear CP2K user,<div><br></div><div>Recently I am using QM/MM to do some MD calculations, the box include a negative QM box and MM water solvents, therefore the net charge is -1. For example, only CO2- in the QM region surrounded by TIP3P water. As far as I understand from the Laino's work in 2005, this unphysical problem has been solved by using a finite width density of charge on MM atoms.</div><div><br></div><div>However, the MM proton still trapped into the quantum oxygen atoms, and MM oxygen, also trapped into quantum carbon atom. The input for my QMMM part is attaed below, any suggestions would be appreciated. </div><div><br></div><div>PS. When I test in pure water bulk (QM water + MM water) system, the result doesn't have this problem. So I guess if I could increase the width of MM atoms? Or apply an artifical replusive potential to avoid the trap? </div><div><br></div><div>Many thanks for your help!</div><div><br></div><div>Xiaohui</div><div><br></div><div><br></div><div><br></div><div><br></div><div><div>  &QMMM</div><div>    &CELL</div><div>      PERIODIC NONE</div><div>      ABC [angstrom] 10 10 10</div><div>    &END CELL</div><div>    &WALLS</div><div>      TYPE NONE</div><div>    &END WALLS</div><div>    &PERIODIC</div><div>      GMAX 0.5</div><div>      &POISSON</div><div>        PERIODIC NONE</div><div>        POISSON_SOLVER MULTIPOLE</div><div>      &END</div><div>      &MULTIPOLE ON</div><div>        EWALD_PRECISION 1.0E-6</div><div>        RCUT 26</div><div>      &END</div><div>    &END PERIODIC</div><div>    &QM_KIND C</div><div>      MM_INDEX 1</div><div>    &END QM_KIND</div><div>    &QM_KIND Oq</div><div>      MM_INDEX 2 3</div><div>    &END QM_KIND</div><div>    &MM_KIND H</div><div>      RADIUS 0.44</div><div>    &END MM_KIND</div><div>    &MM_KIND O</div><div>      RADIUS 0.78</div><div>    &END MM_KIND</div><div>    E_COUPL GAUSS</div><div>    USE_GEEP_LIB 10</div><div>    NOCOMPATIBILITY</div><div>  &END QMMM</div></div></div>