[gmx-users] Re: temperature scaling and frozen atoms
gmx3 at hotmail.com
Thu Jul 31 16:21:01 CEST 2003
Date: Thu, 31 Jul 2003 16:20:05 +0200
Content-Type: text/plain; format=flowed
I the past I have fixed several details concerning freeze groups.
According to me everything should work fine now.
I just checked a system of 4 atoms with 2 atoms frozen and 2 atoms free.
Everything works fine, even with seperate comm removal groups and seperate
temperature coupling groups, using both Berendsen and Nose-Hoover coupling.
When constraints are present in the frozen group things get a bit more
since if the starting structure does not obey the constraints perfectly one
have the whole group frozen and the correct constraint lengths at the same
The current implementation also corrects constraints in frozen groups.
Still I don't see why you can not use comm removal or seperate temperature
Note that comm removal is often not what you want when a large part of the
system is frozen,
since you might not want to restrict motion with respect to the frozen part.
A problem in one of the mdp files you showed is that the tau_t of the
thermostat is 0.02 while the timestep is 0.002. The tau_t sets the period of
the bath oscillations
and 10 timesteps per oscillation is not enough for accurate intergration.
Also you want to have rlist equal to rcoulomb and not to rcoulomb_switch.
Protect your PC - get McAfee.com VirusScan Online
More information about the gromacs.org_gmx-users