[gmx-users] Box exploding with freeze groups

[Chris Neale]

You could try constant volume with freeze groups, then restraints with pressure coupling.

http://www.gromacs.org/pipermail/gmx-users/2002-March/001003.html

Hello everyone,
I'm using Gromacs 3.3.3 and I'm trying to equilibrate the solvent
around a protein. At first I am freezing the C-alphas and restraining
the rest of the protein. The waters are free to move. However, when I
use these freeze groups, the simulation results in the box/solvent
exploding. The simulation doesn't crash, but the solvent takes on a
very distorted form that is not cubical (the geometry of the box) in
shape. I've seen this happen in other systems whenever I use freeze
groups. Does anyone know how to rectify this?
Thanks,
Bob Johnson