(unknown charset) Re: [AMBER] TI calculation using softcore potentials. simultaneous density and temperature equilibration
As the manual script shows, constant V equ is essential to place the adding
solvate to reasonable positions and regulate the temperature at the same
time. I have make a trail to directly heat the system after minimization
with constant pressure, but this work failed for some atoms with unlimited
rate (*****). I assume the heat precedure makes the system unstable when it
go along with constant pressure.
On Thu, Jul 8, 2010 at 2:07 AM, g t <sketchfoot.gmail.com> wrote:
> I'm new to MD and was hoping someone could help me with a query I had
> regarding the amber tutorial:
> TUTORIAL A9: Thermodynamic Integration using soft core
> I have followed it through and got the results I was after, but was
> about why the equilibration phase does not operate under constant volume as
> is suggested in the amber 10 manual.
> "If NTB .NE. 0, there must be a periodic boundary in the topology file.
> pressure is not used in minimization (IMIN=1, above).
> For a periodic system, constant pressure is the only way to equilibrate
> density if the
> starting state is not correct. For example, the solvent packing scheme used
> in LEaP
> can result in a net void when solvent molecules are subtracted which can
> into "vacuum bubbles" in a constant volume run. Another potential problem
> small gaps at the edges of the box. The upshot is that almost every system
> to be equilibrated at constant pressure (ntb=2, ntp>0) to get to a proper
> But be sure to equilibrate first (at constant volume) to something close to
> the final
> temperature, before turning on constant pressure."
> I notice that the tutorial states that temperature & density are being
> equilibrated together. Was this discovered through trial & error?
> will it have to be re-evaluated for every different system? Are there
> general rules when you know that you will invariably have to do two
> equilibration steps for density & temperature?
> Thank you for an assistance you may be able to provide.
> best regards,
> AMBER mailing list
CPU(CHINA PHARMACEUTICAL UNIVERSITY)
Molecular-Design & Drug-Discovery Lab
Laboratory Building A212
639 Longmian Avenue
AMBER mailing list
Received on Fri Jul 09 2010 - 20:00:03 PDT