Abstract:We present a method for Monte Carlo sampling on systems with discrete variables (focusing in the Ising case), introducing a prior on the candidate moves in a Metropolis-Hastings scheme which can significantly reduce the rejection rate, called the reduced-rejection-rate (RRR) method. The method employs same probability distribution for the choice of the moves as rejection-free schemes such as the method proposed by Bortz, Kalos and Lebowitz (BKL) [Bortz et al. <a class="link-external link-http" href="http://J.Comput.Phys" rel="external noopener nofollow">this http URL</a>. 1975]; however, it uses it as a prior in an otherwise standard Metropolis scheme: it is thus not fully rejection-free, but in a wide range of scenarios it is nearly so. This allows to extend the method to cases for which rejection-free schemes become inefficient, in particular when the graph connectivity is not sparse, but the energy can nevertheless be expressed as a sum of two components, one of which is computed on a sparse graph and dominates the measure. As examples of such instances, we demonstrate that the method yields excellent results when performing Monte Carlo simulations of quantum spin models in presence of a transverse field in the Suzuki-Trotter formalism, and when exploring the so-called robust ensemble which was recently introduced in Baldassi et al. [PNAS 2016]. Our code for the Ising case is publicly available [<a class="link-external link-https" href="https://github.com/carlobaldassi/RRRMC.jl" rel="external noopener nofollow">this https URL</a>], and extensible to user-defined models: it provides efficient implementations of standard Metropolis, the RRR method, the BKL method (extended to the case of continuous energy specra), and the waiting time method [Dall and Sibani <a class="link-external link-http" href="http://Comput.Phys.Commun" rel="external noopener nofollow">this http URL</a>. 2001].

Optimized Monte Carlo Methods

Optimized population Monte Carlo

Monte Carlo methods in statistical physics: Mathematical foundations and strategies

Efficient Irreversible Monte Carlo Samplers

Extended Ensemble Monte Carlo

Unbiased Markov chain Monte Carlo methods with couplings

Multicanonical Algorithm, Simulated Tempering, Replica-Exchange Method, and All That

Combinatorial Dehn surgery on cubed and Haken 3-manifolds

A method to reduce the rejection rate in Monte Carlo Markov Chains

Geometric allocation approaches in Markov chain Monte Carlo

Non-Reversible Parallel Tempering: a Scalable Highly Parallel MCMC Scheme

An Efficient Method for Treating Conditional Monte Carlo Simulation

A new algorithm for Monte Carlo simulation of Ising spin systems

An irreversible Markov-chain Monte Carlo method with skew detailed balance conditions

A Variational Approach to Monte Carlo Renormalization Group

From quantum enhanced to quantum inspired Monte Carlo

Wide Sampling and Efficient Updating Monte Carlo Algorithms for Dimer Models

Moment-Preserving and Mesh-Adaptive Reweighting Method for Rare-Event Sampling in Monte-Carlo Algorithms

Reweighting in nonequilibrium simulations

Parallel Tempered Metadynamics: Overcoming potential barriers without surfing or tunneling

A critical study of the Monte Carlo replica method