In this paper we discuss a variation of the multiple-proposal approach to Markov chain-based Monte Carlo methods. We show that by means of intermediate states, used as pivots between the departure configuration and a number of possible arrival configurations, detailed balance is maintained and the acceptance of trial moves is significantly improved with respect to the standard Metropolis Monte Carlo approach. When applied to the case of a particle in a one-dimensional rough potential, this method results in an improved accuracy in the sampling of the distribution function. For all densities of a three-dimensional bulk system of Lennard-Jones particles, it requires less Monte Carlo steps to reach equilibrium and it allows for using wider displacement steps than Metropolis Monte Carlo. In particular, for low-intermediate densities the maximum displacement step can be set to exceptionally large values. We show how biases can be introduced to further reduce the number of operations necessary to construct the trial states while still fulfilling the detailed balance requirement. Although this method is more computationally demanding than (standard) single-proposal Monte Carlo, it samples the configuration space faster. Moreover, its multiple-proposal nature makes it amenable to parallel implementation, while requiring about half the number of operations required by multiple-try Monte Carlo.

Improving the acceptance in Monte Carlo simulations: Sampling through intermediate states / Pazzona, Federico Giovanni; Demontis, Pierfranco; Suffritti, Giuseppe Baldovino. - In: JOURNAL OF COMPUTATIONAL PHYSICS. - ISSN 0021-9991. - 295:(2015), pp. 556-568. [10.1016/j.jcp.2015.04.014]

Improving the acceptance in Monte Carlo simulations: Sampling through intermediate states

PAZZONA, Federico Giovanni;DEMONTIS, Pierfranco;SUFFRITTI, Giuseppe Baldovino
2015-01-01

Abstract

In this paper we discuss a variation of the multiple-proposal approach to Markov chain-based Monte Carlo methods. We show that by means of intermediate states, used as pivots between the departure configuration and a number of possible arrival configurations, detailed balance is maintained and the acceptance of trial moves is significantly improved with respect to the standard Metropolis Monte Carlo approach. When applied to the case of a particle in a one-dimensional rough potential, this method results in an improved accuracy in the sampling of the distribution function. For all densities of a three-dimensional bulk system of Lennard-Jones particles, it requires less Monte Carlo steps to reach equilibrium and it allows for using wider displacement steps than Metropolis Monte Carlo. In particular, for low-intermediate densities the maximum displacement step can be set to exceptionally large values. We show how biases can be introduced to further reduce the number of operations necessary to construct the trial states while still fulfilling the detailed balance requirement. Although this method is more computationally demanding than (standard) single-proposal Monte Carlo, it samples the configuration space faster. Moreover, its multiple-proposal nature makes it amenable to parallel implementation, while requiring about half the number of operations required by multiple-try Monte Carlo.
2015
Improving the acceptance in Monte Carlo simulations: Sampling through intermediate states / Pazzona, Federico Giovanni; Demontis, Pierfranco; Suffritti, Giuseppe Baldovino. - In: JOURNAL OF COMPUTATIONAL PHYSICS. - ISSN 0021-9991. - 295:(2015), pp. 556-568. [10.1016/j.jcp.2015.04.014]
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11388/160637
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 2
  • ???jsp.display-item.citation.isi??? ND
social impact