AMBER Archive (2009)

Subject: RE: [AMBER] assigning velocities to each atom without heating the system

From: Ross Walker (ross_at_rosswalker.co.uk)
Date: Mon Mar 02 2009 - 14:23:22 CST


Hi Parimal,

> I am a novice to amber. I have created inpcrd and prmtop files loading
> the native structure into xleap. Now I want to perform MD. But I don't
> want to heat the system. How can I assign velocities to each atom without
> heating the system? Can we really do it in AMBER?

To follow up on what Gustavo said. If you mean by 'I don't want to heat the system' you mean you want it to REMAIN at 0 K then what Gustavo said applies. I.e. by having kinetic energy you by definition have a temperature. However, if you mean you would like to skip the 'heating' step and instead say just start your system at room temperature. What would be called a 'hot start' then you just need to set tempi=300.0. In this situation if irest=0, i.e. there are no velocities read from a restart file, then sander will generate a Boltzmann distribution of velocities that give 300.0K's worth of kinetic energy and assign these to the atoms.

Alternatively if you want to 'manually' assign velocities to each atom then you would need to modify the inpcrd file to include velocity entries and manually edit each value to give the velocity you want. This is of course tedious but I doubt this is what you meant.

Good luck,
Ross

/\
\/
|\oss Walker

| Assistant Research Professor |
| San Diego Supercomputer Center |
| Tel: +1 858 822 0854 | EMail:- ross_at_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.

_______________________________________________
AMBER mailing list
AMBER_at_ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber