Molecular Dynamics
Molecular dynamics simulates matter atom by atom, integrating Newton's equations of motion under interatomic forces to watch how liquids, solids and biomolecules evolve and to compute their thermodynamic and transport properties.
Definition
Molecular dynamics is a simulation method that computes the trajectory of a system of interacting particles by numerically integrating their classical equations of motion, from which equilibrium and dynamical properties are obtained as time averages.
Scope
This area covers classical molecular dynamics: time integration of the equations of motion with symplectic integrators, the interatomic potentials and force fields that supply the forces, thermostats and barostats that realize statistical ensembles, and the closely related Monte Carlo approach to molecular simulation. It centers on method rather than on any one application domain.
Sub-topics
Core questions
- How are Newton's equations integrated stably for many interacting atoms over long times?
- How are interatomic forces modeled, from simple pair potentials to detailed force fields?
- How are temperature and pressure controlled to simulate a chosen statistical ensemble?
- How are thermodynamic and transport properties extracted from a simulated trajectory?
Key theories
- Trajectory integration
- Molecular dynamics advances positions and velocities with time-reversible symplectic integrators such as the velocity Verlet algorithm, which conserve a shadow energy and keep long simulations stable.
- Force fields and potentials
- Forces derive from interatomic potential energy functions, ranging from the Lennard-Jones pair potential for simple fluids to many-term force fields for molecules, whose accuracy sets the realism of the simulation.
- Ensembles via thermostats and barostats
- Coupling the system to thermostats and barostats modifies the dynamics so that time averages sample the canonical or isothermal-isobaric ensemble rather than the microcanonical one of bare Newtonian dynamics.
Clinical relevance
Molecular dynamics computes diffusion coefficients, viscosities, phase behavior and free energies of fluids and solids, and is a central tool in materials science, soft-matter physics and biomolecular modeling of proteins and membranes.
History
Molecular dynamics began with Alder and Wainwright's hard-sphere simulations in the late 1950s and Rahman's 1964 simulation of liquid argon with a continuous potential; faster computers and better force fields extended it from a few hundred atoms to millions and from simple liquids to biomolecules.
Key figures
- Aneesur Rahman
- Berni Alder
- Daan Frenkel
- Michael P. Allen
Related topics
Seminal works
- rahman1964
- frenkel2002
Frequently asked questions
- How does molecular dynamics differ from Monte Carlo simulation?
- Molecular dynamics follows the real-time trajectory of the system by integrating equations of motion, so it gives dynamical properties like diffusion. Monte Carlo instead samples configurations stochastically and gives equilibrium averages but no genuine time evolution.
- Why are simulation time scales so short?
- Stable integration requires time steps of about a femtosecond to resolve fast atomic vibrations, so even millions of steps cover only nanoseconds to microseconds, which is why bridging to longer biological or material processes is an ongoing challenge.