[gmx-developers] How to treat energy jumps at the cut-off?
hess at kth.se
Tue Feb 28 10:54:06 CET 2012
We are soon releasing a 4.6 pre-release with a new cut-off scheme.
This uses exact, buffered, cut-off's:
In our current git code there is a jump in the energy at the cut-off,
both for LJ and PME. For checking of energy conservation it would be better
to (also) calculate energies with the LJ and PME real-space potential
to zero. But we could also choose to only report the shifted energies.
The issue is that shifting the potential adds a (small) cost and complicates
the kernels a little. If we optionally want to shift it, we need to come
names for the non-shifted and shifted potentials (and they shouldn't be
confused with the "old" shifted interactions, which shift the force).
Do you have an opinion/suggestion?
More information about the gromacs.org_gmx-developers