Université de Reims
INRA UMR FARE 614
Fractionnement des Agro Ressources
CPCB Moulin de la Housse, BP 1039,
F 51687 REIMS cedex 2
at atomistic scale (II/II)
Simulation of a
energy with a
at a given fixed
Buried side chain rotation
Hinge bending of chain
Buried side chain vibration
Water diffusive motion
Surface side chain rotation
Water hindered rotation
Values from McCammon & Harvey (1987) and Eisenberg
They make simulation system seem larger
than it is. On the 2D representation, the chain
in blue is the parent chain, all the other
chains are the images.
Molecular dynamics requires to integrate explicitly
positions of the particles in space and time to
produce a "true" dynamical trajectory (random
walk). We advance the system by some small time
t to recalculate subsequently forces and
velocities, and to fulfill stability requirement.
The most popular integrator is the Verlet algorithm.
Thermostats are needed in situations where
there might be a flux of energy in/out of the
system, e.g. to maintain a constant
temperature. A thermostat is time reversible,
deterministic and goes to the canonical
Interactions with the thermostat are non
they propagate instantaneously.
Some statistical ensembles are listed below,
they are identified by their conserved
· Microcanonical (N,V,E)
· Canonical (N,V,T)
· Grand canonical (
· Isobaric isothermal (N,P,T)
Time step size is limited by the fasted motion you
want to integrate. Time step should be
period of the movement of high frequency for
flexible molecules (usually bond stretching: CH
around 10 fs so use 1 fs). When you simulate an
atomic fluid, the time step should be comparable to
the mean time between collisions.
Cut off Distance
For a system of 1500 atoms, 4500 bonded
terms and 1.1x10
non bonded terms must be
Long range potentials (electrostatic) and also
vdW interactions are often truncated at a
finite cut off distance.
It is necessary to give each atom a velocity. If no
force, new position of atom (at t+
t) would be
determined only by velocity. Forces change the
velocity, complicating things immensely.
distribution using a
oscillators as a bath
(with a bit of help
from the large
United atom approximation is used for simulating
large melts of long chain hydrocarbons, this
model replaces each CH
group with a
united atom that has modified vdW parameters
to take account of the missing hydrogens.
How to minimize the computational effort required?
Remove the highest frequency components from the
dynamics. Bond stretch potential use stiff spring with
replace with rigid rod constraint.
Non Equilibrium Molecular Dynamics
Simplified Force fields
Energy and Entropy
Many errors in modeling can be due
to the users. Be careful in choosing a force
field adapted to its problem and in taking
into account all the parameters of
(non constant temperature dynamics, ...)
Due to the inertia of the atoms, Molecular
Dynamics enables the system to surmount
At each point (with coordinates ) on the
potential energy surface there is a well
defined "energy" . Entropy and free
energy are only defined for distinctly different
"states" e.g. A ("unfolded") and B "(folded").
State B has a lower U
and its minimum is
more probable than
However, state A has
a broader minimum
that can be occupied
in more ways.