Mediaspace scheduled maintenance: Aug 25, 2026 07:00 - 12:00 AM. During this time, videos will be temporarily unavailable. Check status updates.

Lecture

Constraints in Molecular Dynamics: General Method and SHAKE Algorithm

In course

PHYS-403: Computer simulation of physical systems I

The two main topics covered by this course are classical molecular dynamics and the Monte Carlo method. ...

Description

This lecture covers the general method and SHAKE algorithm for constraints in Molecular Dynamics. It explains the evolution by position Verlet, Lagrange parameters, and linearization. The SHAKE algorithm is detailed, including partial updates and considerations for convergence.

Instructor

Alfredo Pasquarello

Alfredo Pasquarello studied physics at the Scuola Normale Superiore of Pisa and at the University of Pisa, obtaining their respective degrees in 1986. He obtained a doctoral degree at the EPFL in 1991 with a thesis on Multiphoton Transitions in Solids . Then, he moved to Bell Laboratories at Murray Hill (New Jersey), where he carried out postdoctoral research on the magnetic properties of carbon fullerenes. In 1993, he joined the Institute for Numerical Research in the Physics of Materials (IRRMA), where his activity involved first-principles simulation methods. In 1998, he was awarded the EPFL Latsis Prize for his research work on disordered silica materials. Succeeding in grant programs of the Swiss National Science Foundation, he then set up his own research group at IRRMA. In July 2003, he is appointed Professor in Theoretical Condensed Matter Physics at EPFL. Currently, he leads the Chair of Atomic Scale Simulation.

Official source

Related concepts (32)

Molecular dynamics

Molecular dynamics (MD) is a computer simulation method for analyzing the physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a fixed period of time, giving a view of the dynamic "evolution" of the system. In the most common version, the trajectories of atoms and molecules are determined by numerically solving Newton's equations of motion for a system of interacting particles, where forces between the particles and their potential energies are often calculated using interatomic potentials or molecular mechanical force fields.

Verlet integration

Verlet integration (vɛʁˈlɛ) is a numerical method used to integrate Newton's equations of motion. It is frequently used to calculate trajectories of particles in molecular dynamics simulations and computer graphics. The algorithm was first used in 1791 by Jean Baptiste Delambre and has been rediscovered many times since then, most recently by Loup Verlet in the 1960s for use in molecular dynamics. It was also used by P. H. Cowell and A. C. C.

Symplectic integrator

In mathematics, a symplectic integrator (SI) is a numerical integration scheme for Hamiltonian systems. Symplectic integrators form the subclass of geometric integrators which, by definition, are canonical transformations. They are widely used in nonlinear dynamics, molecular dynamics, discrete element methods, accelerator physics, plasma physics, quantum physics, and celestial mechanics. Symplectic integrators are designed for the numerical solution of Hamilton's equations, which read where denotes the position coordinates, the momentum coordinates, and is the Hamiltonian.

Euler method

In mathematics and computational science, the Euler method (also called the forward Euler method) is a first-order numerical procedure for solving ordinary differential equations (ODEs) with a given initial value. It is the most basic explicit method for numerical integration of ordinary differential equations and is the simplest Runge–Kutta method. The Euler method is named after Leonhard Euler, who first proposed it in his book Institutionum calculi integralis (published 1768–1870).

Molecular modelling

Molecular modelling encompasses all methods, theoretical and computational, used to model or mimic the behaviour of molecules. The methods are used in the fields of computational chemistry, drug design, computational biology and materials science to study molecular systems ranging from small chemical systems to large biological molecules and material assemblies. The simplest calculations can be performed by hand, but inevitably computers are required to perform molecular modelling of any reasonably sized system.

Related lectures (48)

Molecular dynamics under constraints

Explores molecular dynamics simulations under holonomic constraints, focusing on numerical integration and algorithm formulation.

Constraints in Molecular Dynamics: Theory and SHAKE Algorithm

Explores the theory of constraints in molecular dynamics simulations, including the SHAKE algorithm.

Molecular Dynamics Simulations: Energy Conservation

Explores energy conservation in Hamiltonian systems, numerical integration, time step choices, and constraint algorithms in molecular dynamics simulations.

Molecular dynamics and integrators

Explores molecular dynamics, integrators, trajectory generation, and error monitoring in atomistic modeling.

QM/MM Simulation: Blue Moon Ensemble Theory