Lab Home | Phone | Search
Center for Nonlinear Studies  Center for Nonlinear Studies
 Home 
 People 
 Current 
 Affiliates 
 Visitors 
 Students 
 Research 
 ICAM-LANL 
 Publications 
 Conferences 
 Workshops 
 Sponsorship 
 Talks 
 Colloquia 
 Colloquia Archive 
 Seminars 
 Postdoc Seminars Archive 
 Quantum Lunch 
 Quantum Lunch Archive 
 CMS Colloquia 
 Q-Mat Seminars 
 Q-Mat Seminars Archive 
 P/T Colloquia 
 Archive 
 Kac Lectures 
 Kac Fellows 
 Dist. Quant. Lecture 
 Ulam Scholar 
 Colloquia 
 
 Jobs 
 Postdocs 
 CNLS Fellowship Application 
 Students 
 Student Program 
 Visitors 
 Description 
 Past Visitors 
 Services 
 General 
 
 History of CNLS 
 
 Maps, Directions 
 CNLS Office 
 T-Division 
 LANL 
 
Monday, May 20, 2013
3:00 PM - 4:00 PM
CNLS Conference Room (TA-3, Bldg 1690)

Colloquium

A Verlet-type Algorithm for Improved Simulations of Langevin Equations

Niels Gronbech-Jensen
University of California, Davis

The Stormer-Verlet type numerical integrators for simulating equations of motion are well known and widely used in many contemporary scientific contexts. We review why the simple "Verlet" algorithm, which is a direct second order finite difference approximation to a second order differential, is so desirable for initial value problems with conservation properties. For Langevin dynamics, where coupling to a heat-bath is included through the fluctuation-dissipation, the nature of the conservation requirements change, and the premise for the algorithm is put into question. We present a simple re-derivation of the Stormer-Verlet algorithm, including linear friction with associated stochastic noise. We analytically demonstrate that the new algorithm correctly reproduces diffusive behavior of a particle in a flat potential, and, for a harmonic oscillator, our algorithm provides the exact Boltzmann distribution for any value of damping, frequency, and time step within the usual stability limit. The method, which is as simple as conventional Verlet schemes, is numerically tested on both low-dimensional nonlinear systems as well as more complex systems with many degrees of freedom. Finally, we discuss the opportunities and benefits of proper thermodynamic properties of the numerical integrator for simultaneous accuracy and efficiency in, e.g., Molecular Dynamics simulations.

Host: Turab Lookman