Got the following note when used all-bonds constrained with opls forcefield:
“You are using constraints on all bonds, whereas the forcefield has been parametrized only with constraints involving hydrogen atoms. We suggest using constraints = h-bonds instead, this will also improve performance.”
Does it make sense to have only h-bonds constrained? Why not then all-bonds constrained?