Stochastic differential equation approach to understanding the population control bias in full configuration interaction quantum Monte Carlo
arXiv:2103.07800 · doi:10.1103/PhysRevB.105.235144
Abstract
We investigate a systematic statistical bias found in full configuration quantum Monte Carlo (FCIQMC) that originates from controlling a walker population with a fluctuating shift parameter. This bias can become the dominant error when the sign problem is absent, e.g. in bosonic systems. FCIQMC is a powerful statistical method for obtaining information about the ground state of a sparse and abstract matrix. We show that, when the sign problem is absent, the shift estimator has the nice property of providing an upper bound for the exact ground state energy and all projected energy estimators, while a variational estimator is still an upper bound to the exact energy with substantially reduced bias. A scalar model of the general FCIQMC population dynamics leads to an exactly solvable Itô stochastic differential equation. It provides further insights into the nature of the bias and gives accurate analytical predictions for delayed cross-covariance and auto-covariance functions of the shift energy estimator and the walker number. The model provides a toe-hold on finding a cure for the population control bias. We provide evidence for non-universal power-law scaling of the population control bias with walker number in the Bose-Hubbard model for various estimators of the ground state energy based on the shift or on projected energies. For the specific case of the non-interacting Bose-Hubbard Hamiltonian we obtain a full analytical prediction for the bias of the shift energy estimator.
24 pages, 14 figures; fixed typo in Eq. (72)
References in corpus (15)
- Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms
- Computational complexity and fundamental limitations to fermionic quantum Monte Carlo simulations
- Unbiased Reduced Density Matrices and Electronic Properties from Full Configuration Interaction Quantum Monte Carlo
- Phase diagram of the extended Bose Hubbard model
- Unbiasing the initiator approximation in Full Configuration Interaction Quantum Monte Carlo
- Population size bias in Diffusion Monte Carlo
- The Adaptive Shift Method in Full Configuration Interaction Quantum Monte Carlo: Development and Applications
- The sign problem in full configuration interaction quantum Monte Carlo: Linear and sub-linear representation regimes for the exact wave function
- The effect of quantization on the FCIQMC sign problem
- Population Control Bias and Importance Sampling in Full Configuration Interaction Quantum Monte Carlo
- Minimising biases in Full Configuration Interaction Quantum Monte Carlo
- Eliminating the wave function singularity for ultracold atoms by similarity transformation
- Improved walker population control for full configuration interaction quantum Monte Carlo
- Signatures of the BCS-BEC crossover in the yrast spectra of Fermi quantum rings
- Polaron-Depleton Transition in the Yrast Excitations of a One-Dimensional Bose Gas with a Mobile Impurity