<div dir="ltr">Hi everyone. <div><br></div><div>I have a system (protein in a water box) and I encoutered problems in MM only (via FIST) in CP2K (periodic box)<div><br></div><div>I first used xleap to prep my files, run a minimization, a heating, a equilibration a production run (around 150ns). Everything was fine with amber12 (pmemd)</div><div><br></div><div>After I switched to CP2K (2.5.1) with FIST using amber restart, topology parameters, I run into a problem around 600fs. </div><div><br></div><div>One H of the solvent (TIP3P water) slowly goes near and collides/merges with a H from the backbone of my protein and so go to geom wrong, emax_spline too small.</div><div><br></div><div> I checked FF_INFO and FF_PARAMETERS everything seems in order. When I do the DO_NONBOND keyword to true, the merging didn't append (but other problems).</div></div><div><br></div><div>I wonder how CP2K handles VdW and Electrostatic. The Electrostatic seems to win.</div><div><br></div><div>Anybody have ideas ? Do you want more infos ?</div><div><br></div><div>Thanks for every hints you can give me.</div><div><br>Rolf</div></div>