[gmx-users] Scale down/ up velocities when coupling temperature and pressure?

Yun Shi yunshi09 at gmail.com
Sat Mar 31 18:42:59 CEST 2012

Hi all,

I am just wondering about how GROMACS works when coupling temperature and

Assuming the simplest coupling method, like Berendsen, does GROMACS just
scale down/ up kinetic energy (KE) for each particle in the system, such as
KE1/KE2 = T1/T2 ? Similarly for pressure, KE1/KE2 = P1/P2? So GROMACS
scales the square root of velocities for every particle after each time

BTW, the number of degrees of freedom = 3N - Nc - Ncom according to the
manual. For conventional MD, is Nc equal to the total number of constrained
bonds (like all bonds in gromos ff, and all bonds involving hydrogens in
amber99sb ff)?

Thanks for any explanation.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20120331/9931aa0c/attachment.html>

More information about the gromacs.org_gmx-users mailing list