Population Control Bias and Importance Sampling in Full Configuration Interaction Quantum Monte Carlo

Abstract

Population control is an essential component of any projector Monte Carlo algorithm. This control mechanism usually introduces a bias in the sampled quantities that is inversely proportional to the population size. In this paper, we investigate the population control bias in the full configuration interaction quantum Monte Carlo method. We identify the precise origin of this bias and quantify it in general. We show that it has different effects on different estimators and that the shift estimator is particularly susceptible. We derive a re-weighting technique, similar to the one used in diffusion Monte Carlo, for correcting this bias and apply it to the shift estimator. We also show that by using importance sampling, the bias can be reduced substantially. We demonstrate the necessity and the effectiveness of applying these techniques for sign-problem-free systems where this bias is especially notable. Specifically, we show results for large one-dimensional Hubbard models and the two-dimensional Heisenberg model, where corrected FCIQMC results are comparable to the other high-accuracy results.

arXiv comment8 pages, LaTeX; references and a note added

References in corpus (23)

• Solving the Quantum Many-Body Problem with Artificial Neural Networks
• The ALPS project release 1.3: open source software for strongly correlated systems
• The ALPS project release 2.0: Open source software for strongly correlated systems
• Computational Studies of Quantum Spin Systems
• Finite-Size Scaling of the Ground State Parameters of the Two-Dimensional Heisenberg Model
• The Loop Algorithm
• Deep autoregressive models for the efficient variational simulation of many-body quantum systems
• Algorithms for finite Projected Entangled Pair States
• Density matrix quantum Monte Carlo
• The sign problem and population dynamics in the full configuration interaction quantum Monte Carlo method
• NECI: N-Electron Configuration Interaction with emphasis on state-of-the-art stochastic methods
• Efficient Matrix Product State Method for periodic boundary conditions
• Unbiasing the initiator approximation in Full Configuration Interaction Quantum Monte Carlo
• Introduction to the Bethe ansatz II
• Compact numerical solutions to the two-dimensional repulsive Hubbard model obtained via non-unitary similarity transformations
• Gradient optimization of finite projected entangled pair states
• The Adaptive Shift Method in Full Configuration Interaction Quantum Monte Carlo: Development and Applications
• Projective quantum Monte Carlo simulations guided by unrestricted neural network states
• Minimising biases in Full Configuration Interaction Quantum Monte Carlo
• Understanding Quantum Tunneling using Diffusion Monte Carlo Simulations
• An Efficient Density Matrix Renormalization Group Algorithm for Chains with Periodic Boundary Condition
• Improved walker population control for full configuration interaction quantum Monte Carlo
• Stochastic differential equation approach to understanding the population control bias in full configuration interaction quantum Monte Carlo

Cited by in corpus (15)

• Combined unitary and symmetric group approach applied to low-dimensional spin systems
• Is Externally Corrected Coupled Cluster Always Better than the Underlying Truncated Configuration Interaction?
• Polaron-Depleton Transition in the Yrast Excitations of a One-Dimensional Bose Gas with a Mobile Impurity
• Stochastic differential equation approach to understanding the population control bias in full configuration interaction quantum Monte Carlo
• Quantum Monte Carlo study of the role of p-wave interactions in ultracold repulsive Fermi gases
• Lattice Bose polarons at strong coupling and quantum criticality
• Fixed and partial-node approximations in Slater determinant space for molecules
• Full configuration interaction quantum Monte Carlo for coupled electron--boson systems and infinite spaces
• Real-time Sign-Problem-Suppressed Quantum Monte Carlo Algorithm For Noisy Quantum Circuit Simulations
• Zero-temperature Monte Carlo simulations of two-dimensional quantum spin glasses guided by neural network states
• Spin-glass quantum phase transition in amorphous arrays of Rydberg atoms
• Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data
• Properties of fermionic systems with the Path-integral ground state method
• Quantum Monte Carlo study of thin parahydrogen films on graphite
• Bound excited states of Fröhlich polarons in one dimension