<div dir="ltr">Dear CP2K experts:<div><br></div><div>      Recently I am trying to do QMMM simulation for water system. First I fixed the QM region and try to relax the MM H2O first. But after certain number (~3ps) of MD runs, it shows error info: Geometry is wrong or EMAS_SPLINE is too small. From the trajectory file, one of MM H2O moved too close to QM NH3 (ON distance < 1 angstrom). </div><div><br></div><div>      In my simulation, H2O-NH3 and H2O-N2O4 force fields are both not provided. But I define them as QM region including 4 H2O molecules between them. I used GAUSS coupling method to compute QM-MM interaction, so I think QM-MM electrostatic interaction suppose to prevent H2O and NH3 from moving too close (both O and N are negatively charged). Input files are attached.</div><div><br></div><div>      Any idea will be welcome.</div><div>      Thank you in advance!</div></div>