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

Lecture

Numerical Integration of EOMS

In course

CH-351: Molecular dynamics and Monte-Carlo simulation

Introduction to molecular dynamics and Monte-Carlo simulation methods. ...

Description

This lecture covers the numerical integration of equations of motion (EOMs) using methods like Position Verlet and Leap Frog. It explains the Taylor expansion for position propagation and the importance of conserving total energy in MD simulations.

Instructor

Ursula Röthlisberger

U. Röthlisberger was born in Solothurn (Switzerland). In 1988 she made her diploma in Physical Chemistry in the group of Prof. Ernst Schumacher at the University of Berne (Switzerland). Her Ph.D. thesis was done in collaboration with Dr. Wanda Andreoni at the IBM Zurich Research Laboratory in Rüschlikon. After finishing her Ph.D in 1991 she spent some time as a postdoctoral research assistant at the IBM Research Lab. From 1992-1995 she was a postdoctoral research assistant in the group of Prof. Michael L. Klein at the University of Pennsylvania in Philadelphia (USA). In 1994 she was awarded an advanced researcher fellowship (Profil 2) from the Swiss National Science Foundation. Before starting her Profile 2-fellowship she spent another year as postdoctoral research assistant in the group of Prof. Michele Parrinello at the Max-Planck-Institute for Solid State Physics in Stuttgart, Germany. In 1996 she moved as Profile 2-fellow to the ETH in Zurich, hosted by the group of Prof. Wilfred F. van Gunsteren. In 1997 she became Assistant Professor of Computer-Aided Inorganic Chemistry at the ETH Zurich.

Official source

Related concepts (15)

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.

Numerical integration

In analysis, numerical integration comprises a broad family of algorithms for calculating the numerical value of a definite integral, and by extension, the term is also sometimes used to describe the numerical solution of differential equations. This article focuses on calculation of definite integrals. The term numerical quadrature (often abbreviated to quadrature) is more or less a synonym for numerical integration, especially as applied to one-dimensional integrals.

Numerical analysis

Numerical analysis is the study of algorithms that use numerical approximation (as opposed to symbolic manipulations) for the problems of mathematical analysis (as distinguished from discrete mathematics). It is the study of numerical methods that attempt at finding approximate solutions of problems rather than the exact ones. Numerical analysis finds application in all fields of engineering and the physical sciences, and in the 21st century also the life and social sciences, medicine, business and even the arts.

Numerical methods for ordinary differential equations

Numerical methods for ordinary differential equations are methods used to find numerical approximations to the solutions of ordinary differential equations (ODEs). Their use is also known as "numerical integration", although this term can also refer to the computation of integrals. Many differential equations cannot be solved exactly. For practical purposes, however – such as in engineering – a numeric approximation to the solution is often sufficient. The algorithms studied here can be used to compute such an approximation.

Related lectures (35)

Molecular dynamics under constraints

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

Molecular Dynamics Simulations: Energy Conservation

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

Numerical Analysis: Introduction to Computational Methods

Covers the basics of numerical analysis and computational methods using Python, focusing on algorithms and practical applications in mathematics.

Adiabatic Reactor Design

Covers the design of adiabatic reactors with a focus on isomerization reactions and energy balance calculations.

Explicit Stabilised Methods: Stochastic Differential Equations