Hi all,<br><br>I am just wondering about how GROMACS works when coupling temperature and pressure.<br><br>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 step?<br>

<br><br>-----------------------------------------------------------------<br>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)?<br>

<br>Thanks for any explanation.<br><br>Yun<br>