Molecular Dynamics Simulation

Molecular dynamics integrates Newton's equations of motion for a system of interacting atoms, generating trajectories from which structural, dynamic, and thermodynamic properties are computed.

Definition

A simulation technique that propagates atomic positions and velocities in time under the forces from a potential-energy function, producing a deterministic trajectory through phase space.

Scope

Covers the numerical integration of the classical equations of motion, time-reversible integrators such as the Verlet family, control of temperature and pressure through thermostats and barostats, treatment of long-range electrostatics and periodic boundaries, and the extraction of equilibrium and transport properties from trajectories.

Core questions

• How is Newton's equation of motion integrated stably over many time steps?
• How are temperature and pressure controlled to sample a chosen thermodynamic ensemble?
• How are long-range electrostatic interactions handled efficiently?
• How are observable properties recovered from a finite trajectory?

Key theories

Verlet integration

A simple, time-reversible, and symplectic scheme for integrating the equations of motion that conserves energy well over long simulations and underlies most molecular dynamics codes.

Ergodic sampling

Under the ergodic hypothesis, time averages along a sufficiently long trajectory equal ensemble averages, linking simulation to statistical-mechanical observables.

Mechanisms

Each step computes forces from the potential gradient, advances positions and velocities by the integrator, and applies thermostat or barostat corrections; repeating this builds a trajectory whose averages give thermodynamic and dynamic properties.

Clinical relevance

Molecular dynamics reveals conformational changes, binding events, diffusion, and reaction-adjacent processes in proteins, membranes, and materials, providing mechanistic insight that complements experiment in biophysics and drug discovery.

History

Following Alder and Wainwright's hard-sphere simulations and Rahman's continuous-potential study of liquid argon, Verlet's 1967 algorithms and the later development of thermostats, barostats, and Ewald-based electrostatics established molecular dynamics as a mature simulation discipline.

Key figures

• Loup Verlet
• Berni Alder
• Aneesur Rahman
• Herman Berendsen

Seminal works

• verlet1967
• frenkel2002

Frequently asked questions

What limits the time step in molecular dynamics?

The fastest motions, typically bond vibrations involving hydrogen, set the limit; constraining those bonds allows a time step of a few femtoseconds, while accessible total times reach microseconds or beyond.

Is molecular dynamics deterministic?

The integration is deterministic given initial conditions, but trajectories are chaotic, so tiny differences grow rapidly; meaningful results come from statistical averages rather than individual paths.