Re: AMBER: restraints and constant pressure

From: Myunggi Yi <>
Date: Sun, 17 Sep 2006 10:23:33 -0400

Dr. Cheatham said,

To simplify the restraint code, rather than inducing a force (coordinate
deviation) every time the simulation coordinates are scaled (relative to
the restraint coordinates), it was decided to scale the restraint
coordinates as well. This effectively changes the restraint coordinates
at each step (every time the coordinates are scaled). Two problems may

(1) when you restart, and reference previous restraint coordinates (that
were not shifted), you will get an initially large restraint energy due to
the "shift", and

(2) if you are restraining multiple molecules, the relative positions of
the molecules will change as the simulation proceeds. Imagine the default
case where the density of the water is too low initially, the box has to
contract. If you are restraining two molecules, for example the two
strands of DNA, as the run continues, the restraint coordinates will shift
the molecules closer together which will force your MD coordinates to
cause the two strands to get closer. This is undesirable.

Best wishes,
KLB 419
Institute of Molecular Biophysics
Florida State University
Tallahassee, FL 32306
Office: (850) 645-1334
The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Wed Sep 20 2006 - 06:07:04 PDT
Custom Search