University of Cambridge > Talks.cam > DAMTP Statistical Physics and Soft Matter Seminar > Proper Statistical Sampling in Isothermal-Isobaric Discrete-Time Molecular Dynamics

Proper Statistical Sampling in Isothermal-Isobaric Discrete-Time Molecular Dynamics

Add to your list(s) Download to your calendar using vCal

  • UserProf. Oded Farago (Ben Gurion and Cambridge Chemistry)
  • ClockTuesday 22 January 2019, 13:00-14:00
  • HouseMR11, CMS.

If you have a question about this talk, please contact Professor Mike Cates.

Abstract: Molecular Dynamics simulations always involve a discretization of time, but the discrete-time behavior is increasingly different from that of the continuous-time physical equations as the time step is increased. This fact creates a dilemma for any simulation of a dynamical system: Use a small time step, resulting in dynamics that resemble continuous-time behavior at the expense of efficiency; or use a large time step that makes the simulation finish sooner at the expense of meaningful evolution. It is, therefore, essential to understand the features of different algorithms, such that optimal properties can be chosen for a given set of problems and objectives.

Our aim is to investigate and improve Molecular Dynamics simulation techniques for systems in thermal equilibrium. I will present a simple derivation of a stochastic Stormer-Verlet algorithm for the evolution of Langevin equations. The method, which is as simple as conventional Verlet schemes, has been numerically tested on both low-dimensional nonlinear systems as well as more complex molecular ensembles with many degrees of freedom. In light of the fundamental artifacts introduced by discrete time to dynamical simulations, I will provide a simple intuitive picture of the unique benefits of our algorithm that, unlike other algorithms, preserves proper configurational sampling (diffusion and Boltzmann distribution) in discrete time. I will also introduce a companion algorithm for controlling pressure in molecular ensembles, i.e., a barostat for NPT simulations.

This talk is part of the DAMTP Statistical Physics and Soft Matter Seminar series.

Tell a friend about this talk:

This talk is included in these lists:

Note that ex-directory lists are not shown.

 

© 2006-2024 Talks.cam, University of Cambridge. Contact Us | Help and Documentation | Privacy and Publicity