Remove velocity from partially frozen atoms in md-vv
authorPascal Merz <pascal.merz@me.com>
Wed, 3 Mar 2021 10:48:27 +0000 (10:48 +0000)
committerMark Abraham <mark.j.abraham@gmail.com>
Wed, 3 Mar 2021 10:48:27 +0000 (10:48 +0000)
commite3381bcf19b2ec165f571ac384b3ad0e596a7eb5
tree7923b2b8dfbdd14ef410091e200337215ec39325
parent62e69e247929952691e522eb31a7fc813e3d3f36
Remove velocity from partially frozen atoms in md-vv

md-vv would add some velocity to the frozen dimensions of partially
frozen atoms during constraining. This did not lead to wrong
trajectories, as the frozen dimensions of the positions are kept fixed
during propagation. The non-zero velocities were, however, reported in
trajectories and final configurations. They might also have lead to
slightly wrong kinetic energies, since the reported kinetic energy is
calculated after the velocities are constrained. All effects are
expected to be relatively small, since they did not accumulate, as the
velocities were regularly reset to zero once per step.

Refs #3849
docs/release-notes/2021/2021.1.rst
src/gromacs/mdlib/constr.cpp