[gmx-developers] zero temperature II

Michael Brunsteiner mbx0009 at yahoo.com
Sun Apr 4 15:31:45 CEST 2004

dear all,
In a previous mail I mentioned a problem
with zero temperature output. I seem to have 
found what the problem is and this looks like
a bug, or at least an undocumented limitation.
When I do md with a system containing only frozen atoms 
and (not-frozen) water both ener.edr and md.log say 
that the kinetic energy and the temperature are
zero. The thermostate seems to work though, because
my system does not "explode".
If I "un-freeze" the frozen atoms the
output is back to normal. The fact that frozen atoms
are formally provided with an infinite mass
might somehow cause this problem.
Any ideas for a work-around ?
Where in the code is E_kin calculated ?

Do you Yahoo!?
Yahoo! Small Business $15K Web Design Giveaway 

More information about the gromacs.org_gmx-developers mailing list