Optimized Broad-Histogram Ensembles for the Simulation of Quantum Systems
Abstract
The efficiency of statistical sampling in broad-histogram Monte Carlo simulations can be considerably improved by optimizing the simulated extended ensemble for fastest equilibration. Here we describe how a recently developed feedback algorithm can be generalized to find optimized sampling distributions for the simulation of quantum systems in the context of the stochastic series expansion (SSE) when defining an extended ensemble in the expansion order. If the chosen update method is efficient, such as non-local updates for systems undergoing a second-order phase transition, the optimized ensemble is characterized by a flat histogram in the expansion order if a variable-length formulation of the SSE is used. Whenever the update method suffers from slowdown, such as at a first-order phase transition, the feedback algorithm shifts weight towards the expansion orders in the transition region, resulting in a non-uniform histogram.
Cite
@article{arxiv.0710.0024,
title = {Optimized Broad-Histogram Ensembles for the Simulation of Quantum Systems},
author = {Stefan Wessel and Norbert Stoop and Emanuel Gull and Simon Trebst and Matthias Troyer},
journal= {arXiv preprint arXiv:0710.0024},
year = {2007}
}
Comments
15 pages, 13 figures