The cat and the pendulum: area preservation of the flow of Hamiltonian systems

Symplectic transformation

Symplectic mapping

Poincaré theorem: The Hamiltonian evolution is a symplectic mapping (generalization of Liouville theorem). The theorem is shown in figure for the simple example of the mathematical pendulum

Symplectic integrators

Part II: The ergodic hypothesis

Part III: The hard-sphere model

Short history:

The study of hard-sphere systems goes back a long time. Many people see a precursor in the Roman poet and philosopher Lucretius. Another mile-stone was Daniel Bernoulli's (1700 – 1782) discussion of the pressure dependence of a hard-sphere gas (1738). Boltzmann worked on hard spheres, so did Maxwell, and many researchers since then. The phase transition in hard disks, discovered by Alder and Wainwright ^{[1]} in 1962 is an important discovery made by numerical simulations using the event-driven algorithm.
Event driven dynamics was invented by Alder and Wainwright, in 1957. Here is a cartoon of the time-evolution of four hard disks in a square box.

Event-driven Molecular Dynamics simulation for 4 disks in a box

At time t=0, each disk starts at a given position and with a given velocity. The entire time-evolution is simply the solution of Newton's equations: Each disk moves freely until it collides either with another disk or with a wall. In this simulation, time is continuous, and the time for the next "event" can be computed exactly: it is given by the minimum of N(N-1)/2 pair collision times for isolated pairs of particles, without walls, and N wall-collision times (for isolated disks). The figure on the top is SMAC fig. 2.1.

How to test ergodicity? Explain the movie

Molecular Dynamics evolution for four hard disks in a box with walls (simulation by Maxim Berman).

Chaos

The same MD simulation executed with different precisions of arithmetic

The event-driven molecular dynamics algorithm has no time-step error, and the only source of error comes from the finite precision of the arithmetic. These errors are magnified from iteration to iteration. The cartoons on the left illustrates the influence of tiny rounding errors on the dynamics. This is SMAC fig. 2.5. This extreme influence on the initial conditions is called ''chaos''. In this simulation, chaos has two consequences:

Molecular dynamics cannot really give rigorous results for the dynamics: How can we compute the exact positions of particles at t = 100000 if for t = 33, we already need double precision?

The hard-sphere system is chaotic for all radii. Chaos implies that statistical mechanics applies to a finite system of hard spheres and hard disks. This was proven in fundamental theorems by Sinai^{[2]} and Simanyi^{[3]} . The mathematical problems in these works are formidable, however, and the situation does not seem to be fully under control.

Consider a system described by a Hamiltonian dynamics

p˙n=−∂∂qn q˙n=∂∂pn
where qn(t) represents a positions and pn(t) a momenta and qn(0), pn(0) are the initial conditions. Here the total energy E=(qn(0),pn(0))
is a constant of the motion.

The cat and the pendulum: area preservation of the flow of Hamiltonian systems

Symplectic transformation

Symplectic mapping

Poincaré theorem: The Hamiltonian evolution is a symplectic mapping (generalization of Liouville theorem). The theorem is shown in figure for the simple example of the mathematical pendulum

The study of hard-sphere systems goes back a long time. Many people see a precursor in the Roman poet and philosopher Lucretius. Another mile-stone was Daniel Bernoulli's (1700 – 1782) discussion of the pressure dependence of a hard-sphere gas (1738). Boltzmann worked on hard spheres, so did Maxwell, and many researchers since then. The phase transition in hard disks, discovered by Alder and Wainwright in 1962 [1] is an important discovery made by numerical simulations using the event-driven algorithm.

Event driven dynamics was invented by Alder and Wainwright, in 1957. Here is a cartoon of the time-evolution of four hard disks in a square box.

Event_movie.jpg

Event-driven Molecular Dynamics simulation for 4 disks in a box

At time t=0, each disk starts at a given position and with a given velocity. The entire time-evolution is simply the solution of Newton's equations: Each disk moves freely until it collides either with another disk or with a wall. In this simulation, time is continuous, and the time for the next "event" can be computed exactly: it is given by the minimum of N(N-1)/2 pair collision times for isolated pairs of particles, without walls, and N wall-collision times (for isolated disks). The figure on the top is SMAC fig. 2.1.

How to test ergodicity? Explain the movie

Event_chain_box.gif

Molecular Dynamics evolution for four hard disks in a box with walls (simulation by Maxim Berman).

Chaos

^ B. J. Alder and T. E. Wainwright, Phase Transition in Elastic Disks Phys. Rev. 127, 359 (1962)

## Molecular dynamics - symplectic integrators and event driven dynamics

## Part I: Hamilton Dynamics

Consider a system described by a Hamiltonian dynamicswhere

qn(t) represents a positions andpn(t) a momenta andqn(0),pn(0) are the initial conditions. Here the total energyis a constant of the motion.

## Euler Discretization

## Symplectic and Poincaré theorem (1899)

## Symplectic integrators

## Part II: The ergodic hypothesis

## Part III: The hard-sphere model

## Short history:

The study of hard-sphere systems goes back a long time. Many people see a precursor in the Roman poet and philosopher Lucretius. Another mile-stone was Daniel Bernoulli's (1700 – 1782) discussion of the pressure dependence of a hard-sphere gas (1738). Boltzmann worked on hard spheres, so did Maxwell, and many researchers since then. The phase transition in hard disks, discovered by Alder and Wainwright^{[1]}in 1962 is an important discovery made by numerical simulations using the event-driven algorithm.Event driven dynamics was invented by Alder and Wainwright, in 1957. Here is a cartoon of the time-evolution of four hard disks in a square box.

At time t=0, each disk starts at a given position and with a given velocity. The entire time-evolution is simply the solution of Newton's equations: Each disk moves freely until it collides either with another disk or with a wall. In this simulation, time is continuous, and the time for the next "event" can be computed exactly: it is given by the minimum of

N(N-1)/2 pair collision times for isolated pairs of particles, without walls, and N wall-collision times (for isolated disks). The figure on the top is SMAC fig. 2.1.## How to test ergodicity? Explain the movie

## Chaos

The event-driven molecular dynamics algorithm has no time-step error, and the only source of error comes from the finite precision of the arithmetic. These errors are magnified from iteration to iteration. The cartoons on the left illustrates the influence of tiny rounding errors on the dynamics. This is SMAC fig. 2.5. This extreme influence on the initial conditions is called ''chaos''. In this simulation, chaos has two consequences:

t= 100000 if fort= 33, we already need double precision?^{[2]}and Simanyi^{[3]}. The mathematical problems in these works are formidable, however, and the situation does not seem to be fully under control.- Wiki Home
- Projects
- Recent Changes
- Pages and Files
- Members
- Settings

Here's how to change theNavigation page

Python

edit navigation

## Molecular dynamics - symplectic integrators and event driven dynamics

## Part I: Hamilton Dynamics

Consider a system described by a Hamiltonian dynamicsp˙n=−∂∂qnq˙n=∂∂pnwhere

qn(t) represents a positions andpn(t) a momenta andqn(0),pn(0) are the initial conditions. Here the total energyE=(qn(0),pn(0))is a constant of the motion.

## Euler Discretization

pt+Δt=pt−Δt(∂∂q)q=qt,p=ptqt+Δt=qt+Δt(∂∂p)q=qt,p=pt## Symplectic and Poincaré theorem (1899)

Symplectic transformationSymplectic mappingPoincaré theorem:The Hamiltonian evolution is a symplectic mapping (generalization of Liouville theorem). The theorem is shown in figure for the simple example of the mathematical pendulumpendulum=p22−cosq## Symplectic integrators

pt+Δt=pt−Δt(∂∂q)q=qt,p=pt+Δtqt+Δt=qt+Δt(∂∂p)q=qt,p=pt+Δt## Part II: The ergodic hypothesis

## Part III: The hard-sphere model

## Short history:

The study of hard-sphere systems goes back a long time. Many people see a precursor in the Roman poet and philosopher Lucretius. Another mile-stone was Daniel Bernoulli's (1700 – 1782) discussion of the pressure dependence of a hard-sphere gas (1738). Boltzmann worked on hard spheres, so did Maxwell, and many researchers since then. The phase transition in hard disks, discovered by Alder and Wainwright in 1962 [1] is an important discovery made by numerical simulations using the event-driven algorithm.

Event driven dynamics was invented by Alder and Wainwright, in 1957. Here is a cartoon of the time-evolution of four hard disks in a square box.

At time t=0, each disk starts at a given position and with a given velocity. The entire time-evolution is simply the solution of Newton's equations: Each disk moves freely until it collides either with another disk or with a wall. In this simulation, time is continuous, and the time for the next "event" can be computed exactly: it is given by the minimum of

N(N-1)/2 pair collision times for isolated pairs of particles, without walls, and N wall-collision times (for isolated disks). The figure on the top is SMAC fig. 2.1.## How to test ergodicity? Explain the movie

## Chaos

Phase Transition in Elastic DisksPhys. Rev. 127, 359 (1962)Dynamical systems with elastic reflections, [http://iopscience.iop.org/0036-0279/25/2/R05 Russian Mathematical Surveys 25, 137-189 (1970)]Proof of the Boltzmann-Sinai ergodic hypothesis for typical hard disk systems[http://www.springerlink.com/content/yf9uw55m3vm6utqx/ Inventiones Mathematicae 154, 123-178 (2003)]