Optimal sampling efficiency in Monte Carlo sampling with an approximate potential
Journal Article
·
· Journal of Chemical Physics
- Los Alamos National Laboratory
- U MISSOURI
Building on the work of Iftimie et al., Boltzmann sampling of an approximate potential (the 'reference' system) is used to build a Markov chain in the isothermal-isobaric ensemble. At the endpoints of the chain, the energy is evaluated at a higher level of approximation (the 'full' system) and a composite move encompassing all of the intervening steps is accepted on the basis of a modified Metropolis criterion. For reference system chains of sufficient length, consecutive full energies are statistically decorrelated and thus far fewer are required to build ensemble averages with a given variance. Without modifying the original algorithm, however, the maximum reference chain length is too short to decorrelate full configurations without dramatically lowering the acceptance probability of the composite move. This difficulty stems from the fact that the reference and full potentials sample different statistical distributions. By manipulating the thermodynamic variables characterizing the reference system (pressure and temperature, in this case), we maximize the average acceptance probability of composite moves, lengthening significantly the random walk between consecutive full energy evaluations. In this manner, the number of full energy evaluations needed to precisely characterize equilibrium properties is dramatically reduced. The method is applied to a model fluid, but implications for sampling high-dimensional systems with ab initio or density functional theory (DFT) potentials are discussed.
- Research Organization:
- Los Alamos National Laboratory (LANL)
- Sponsoring Organization:
- DOE
- DOE Contract Number:
- AC52-06NA25396
- OSTI ID:
- 956349
- Report Number(s):
- LA-UR-09-00498; LA-UR-09-498
- Journal Information:
- Journal of Chemical Physics, Journal Name: Journal of Chemical Physics; ISSN JCPSA6; ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Optimized nested Markov chain Monte Carlo sampling: theory
An efficient approach to ab initio Monte Carlo simulation
Optimized nested Markov chain Monte Carlo sampling: application to the liquid nitrogen Hugoniot using density functional theory
Conference
·
Wed Dec 31 23:00:00 EST 2008
·
OSTI ID:989792
An efficient approach to ab initio Monte Carlo simulation
Journal Article
·
Mon Jan 20 23:00:00 EST 2014
· Journal of Chemical Physics
·
OSTI ID:22255277
Optimized nested Markov chain Monte Carlo sampling: application to the liquid nitrogen Hugoniot using density functional theory
Journal Article
·
Wed Dec 31 23:00:00 EST 2008
·
OSTI ID:991274