[gmx-users] bond deviation

[sadhna]

hi,
    I am trying to simulate protein in water for 1ns. I want to
equilibrate solvent before simulation because the protein is not
maintaining the desired structure due to large solvent forces.
I tried doing it by giving following in mdp file for em in water:
freezegrps           = protein SOL
freezedim            = Y Y Y N N N
 
but after starting mdrun for 1ns i am geeting fatal error;Bond deviates
more than half its own length

please tell me workaround for this problem

thanks
sadhna