<div dir="ltr"><div class="IVILX2C-ed-a"><div tabindex="0" class="IVILX2C-nb-P"><div style="overflow: auto;"><div style="max-height: 10000px;"><div dir="ltr">Hi everyone,<div><br><div>I'm running a simulation of the Chorismate Mutase complex in water with CP2K. The enzyme and waters are treated at the MM level, while one copy of the chorismate is treated at the QM level (PM6). When I run the system with ECOUPL NONE, everything works fine. However, when I try to switch to Coulomb coupling, the system crashes at around 950 steps. The QM temperate skyrockets and I get the error:</div><div><br></div><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;">Particles:    5766   4633 at distance [au]:     0.79662080 less than:      0.95400000; increase EMAX_SPLINE. <br>GEOMETRY wrong or EMAX_SPLINE too small </blockquote><div><br></div><div> I suspect that the error is due to the interaction of the tip3p waters with the negatively charged QM region. I tried to fix this by giving them a radius through MM_KIND RADIUS, but that unfortunately did not work. I have also tried adding LJ-interactions between the negative carboxylate and the water hydrogens (currently commented out in my input file) to no success. </div><div><br></div><div>I have uploaded input and output files to dropbox  (<a href="https://www.dropbox.com/sh/du6f32omx32d7a5/AAAiYYttGmwrDj5V_v3y3QZKa?dl=0" target="_blank" rel="nofollow" style="cursor: pointer;">https://www.dropbox.com/sh/<wbr>du6f32omx32d7a5/<wbr>AAAiYYttGmwrDj5V_v3y3QZKa?dl=0</a><wbr>) as they are too big for the forum. I would greatly appreciate some help on what causes this error and how to fix it. Please let me know if you require any additional information to help troubleshoot this issue. </div></div><div><br></div><div>Kind regards</div><div><br></div><div>Dries</div></div></div></div></div><div></div><div></div></div><div><div class="IVILX2C-ed-a"></div></div><div class="IVILX2C-nb-b"><div class="IVILX2C-nb-a IVILX2C-nb-cb"><div style="display: inline-block;"></div><div class="IVILX2C-md-a"></div></div></div></div>