Hi all,<div><br /></div><div>I'm trying to perform QM/MM umbrella sampling on a small molecule system in water using semiempirical QM methods in cp2k and gromacs for MD. I keep running into SCF convergence issues due to bonds to hydrogen atoms breaking in my QM region and likely leaving the QM cell since semiempirical methods don't support a periodic QM region.<br /><br />I've reduced my time step from 1fs to 0.5fs and have tried varying the size of my QM cell (even up to the full box size) and reducing the KAPPA value for my two RC constraints.<br /><br />I've done simulations on the same system using DFT and was not seeing these issues, so it's possible I'm missing a keyword that's relevant when using semiempirical methods. My cp2k input file, gromacs mdp file, and plumed file are attached. I would really appreciate any suggestions!<br /><br />Thank you,</div><div>Emma</div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/aae6a2fe-4516-4556-b2d6-c274255c6182n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/aae6a2fe-4516-4556-b2d6-c274255c6182n%40googlegroups.com</a>.<br />