RE: AMBER: AMBER : About the kinetic energy and temperature (GBSA model)

From: Ross Walker <ross.rosswalker.co.uk>
Date: Tue, 28 Mar 2006 20:35:03 -0800

Dear James,

> simulated my system again .The temperature is still very high
> ~ 5123 k . It's
> not reasonable .Why "Sander" didn't correct the
> temperature?Could you give me
> more information about the temperature regulation ,such as
> an example ?

It is very reasonable... The temperature regulation is simply a weak
thermostat, not a rigid constant. If your system's forces are too large it
can easily 'explode' and no amount of temperature regulation will save you.
I can set a bomb off in a freezer and it will still go bang...

A couple of things to try. If you haven't already done so I suggest you work
through the following tutorial.
http://www.rosswalker.co.uk/tutorials/amber_workshop/Tutorial_one/index.htm
this will give you an understanding of how to minimise a system and then
heat it up in controlled fashion. It will also give you an understanding of
shake and timesteps.

Once you have done that go back to your simulation and set ntpr=1 and ntwx=1
and run it for 1000 steps or so nstlim=1000. Then open the resulting mdcrd
file in VMD and it should be immediately obvious what is wrong with your
simulation.

All the best
Ross

/\
\/
|\oss Walker

| HPC Consultant and Staff Scientist |
| San Diego Supercomputer Center |
| Tel: +1 858 822 0854 | EMail:- ross.rosswalker.co.uk |
| http://www.rosswalker.co.uk | PGP Key available on request |

Note: Electronic Mail is not secure, has no guarantee of delivery, may not
be read every day, and should not be used for urgent or sensitive issues.


-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Apr 02 2006 - 06:10:10 PDT
Custom Search