| Literature DB >> 29119281 |
Shahrazad M A Malek1, Richard K Bowles2, Ivan Saika-Voivod1, Francesco Sciortino3, Peter H Poole4.
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 [Formula: see text], we can monitor the swarm's relaxation to equilibrium, and confirm its attainment, within [Formula: see text], where [Formula: see text] 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.Keywords: Topical issue: Advances in Computational Methods for Soft Matter Systems
Mesh:
Substances:
Year: 2017 PMID: 29119281 DOI: 10.1140/epje/i2017-11588-2
Source DB: PubMed Journal: Eur Phys J E Soft Matter ISSN: 1292-8941 Impact factor: 1.890