"Swarm relaxation": Equilibrating a large ensemble of computer simulations

Abstract

It is common practice in molecular dynamics and Monte Carlo computer simulations to run multiple, separately-initialized simulations in order to improve the sampling of independent microstates. Here we examine the utility of an extreme case of this strategy, in which we run a large ensemble of $M$ independent simulations (a "swarm"), each of which is relaxed to equilibrium. We show that if $M$ is of order $10^3$, we can monitor the swarm's relaxation to equilibrium, and confirm its attainment, within $\sim 10\bar{\tau }$, where $\bar{\tau }$ is the equilibrium relaxation time. As soon as a swarm of this size attains equilibrium, the ensemble of $M$ final microstates from each run is sufficient for the evaluation of most equilibrium properties without further sampling. This approach dramatically reduces the wall-clock time required, compared to a single long simulation, by a factor of several hundred, at the cost of an increase in the total computational effort by a small factor. It is also well-suited to modern computing systems having thousands of processors, and is a viable strategy for simulation studies that need to produce high-precision results in a minimum of wall-clock time. We present results obtained by applying this approach to several test cases.

arXiv comment12 pages. To appear in Eur. Phy. J. E, 2017

References in corpus (9)

• Theoretical perspective on the glass transition and amorphous materials
• Aging Effects in a Lennard-Jones Glass
• Free energy surface of ST2 water near the liquid-liquid phase transition
• Two-state thermodynamics of the ST2 model for supercooled water
• On the Study of Collective Dynamics in Supercooled Liquids through the Statistics of the Iso-Configurational Ensemble
• Kovacs effects in an aging molecular liquid
• Dynamical Behavior Near a Liquid-Liquid Phase Transition in Simulations of Supercooled Water
• Narayanaswamy's 1971 aging theory and material time
• Investigation of the relation between local diffusivity and local inherent structures in the Kob-Andersen Lennard-Jones model

Cited by in corpus (2)

• RUMD: A general purpose molecular dynamics package optimized to utilize GPU hardware down to a few thousand particles
• Surface tension of supercooled water nanodroplets from computer simulations