Unbiasing the initiator approximation in Full Configuration Interaction Quantum Monte Carlo

Abstract

We identify and rectify a crucial source of bias in the initiator FCIQMC algorithm. Non-initiator determinants (i.e. determinants whose population is below the initiator threshold) are subject to a systematic {\em undersampling} bias, which in large systems leads to a bias in the energy when an insufficient number of walkers is used. We show that the acceptance probability ($p_{acc}$), that a non-initiator determinant has its spawns accepted, can be used to unbias the initiator bias, in a simple and accurate manner, by reducing the applied shift to the non-initiator proportionately to $p_{acc}$. This modification preserves the property that in the large walker limit, when $p_{acc}\to 1$, the unbiasing procedure disappears, and the initiator approximation becomes exact. We demonstrate that this algorithm shows rapid convergence to the FCI limit with respect to walker number, and furthermore largely removes the dependence of the algorithm on the initiator threshold, enabling highly accurate results to be obtained even with large values of the threshold. This is exemplified in the case of butadiene/ANO-L-pVDZ and benzene/cc-pVDZ, correlating 22 and 30 electrons in 82 and 108 orbitals respectively. In butadiene $5\times 10^7$ and in benzene $10^8$ walkers suffice to obtain an energy to within a milli-Hartree of the CCSDT(Q) result, in Hilbert spaces of $10^{26}$ and $10^{35}$ respectively. Essentially converged results require $\sim 10^8$ walkers for butadiene and $\sim 10^9$ walkers for benzene, and lie slightly lower than CCSDT(Q). Owing to large-scale parallelisability, these calculations can be executed in a matter of hours on a few hundred processors. The present method largely solves the initiator-bias problems that the initiator method suffered from when applied to medium-sized molecules.

12 pages, 4 figures

References in corpus (22)

• Heat-bath Configuration Interaction: An efficient selected CI algorithm inspired by heat-bath sampling
• Semistochastic Heat-bath Configuration Interaction method: selected configuration interaction with semistochastic perturbation theory
• A deterministic alternative to the full configuration interaction quantum Monte Carlo method
• An adaptive configuration interaction approach for strongly correlated electrons with tunable accuracy
• Semistochastic Projector Monte Carlo Method
• Hybrid stochastic-deterministic calculation of the second-order perturbative contribution of multireference perturbation theory
• Unbiased Reduced Density Matrices and Electronic Properties from Full Configuration Interaction Quantum Monte Carlo
• Excited States of Methylene, Polyenes, and Ozone from Heat-Bath Configuration Interaction
• Selected configuration interaction dressed by perturbation
• Adaptive multiconfigurational wave functions
• Linear-scaling and parallelizable algorithms for stochastic quantum chemistry
• An excited-state approach within full configuration interaction quantum Monte Carlo
• The sign problem and population dynamics in the full configuration interaction quantum Monte Carlo method
• Semi-stochastic full configuration interaction quantum Monte Carlo: developments and application
• Efficient Formulation of Full Configuration Interaction Quantum Monte Carlo in a Spin Eigenbasis via the Graphical Unitary Group Approach
• Virtual orbital many-body expansions: A possible route towards the full configuration interaction limit
• Compact numerical solutions to the two-dimensional repulsive Hubbard model obtained via non-unitary similarity transformations
• Similarity transformation of the electronic Schrödinger equation via Jastrow factorisation
• Stochastic multi-reference perturbation theory with application to linearized coupled cluster method
• An efficient and accurate perturbative correction to initiator full configuration interaction quantum Monte Carlo
• Multi-state effective Hamiltonian and size-consistency corrections in stochastic configuration interactions
• Time propagation and spectroscopy of Fermionic systems using a stochastic technique

Cited by in corpus (40)

• QUESTDB: a database of highly-accurate excitation energies for the electronic structure community
• The Ground State Electronic Energy of Benzene
• Variational and Diffusion Quantum Monte Carlo Calculations with the CASINO Code
• NECI: N-Electron Configuration Interaction with emphasis on state-of-the-art stochastic methods
• Prospects of Quantum Computing for Molecular Sciences
• The Shape of Full Configuration Interaction to Come
• The performance of CIPSI on the ground state electronic energy of benzene
• Electronic correlations and magnetic interactions in infinite-layer NdNiO${}_2$
• The Adaptive Shift Method in Full Configuration Interaction Quantum Monte Carlo: Development and Applications
• Accurate full configuration interaction correlation energy estimates for five- and six-membered rings
• Molecule-specific Uncertainty Quantification in Quantum Chemical Studies
• Population Control Bias and Importance Sampling in Full Configuration Interaction Quantum Monte Carlo
• Is Externally Corrected Coupled Cluster Always Better than the Underlying Truncated Configuration Interaction?
• Eliminating the wave function singularity for ultracold atoms by similarity transformation
• A new form of transcorrelated Hamiltonian inspired by range-separated DFT
• Towards near-exact solutions of molecular electronic structure: Full coupled-cluster reduction with a second-order perturbative correction
• Taming the fixed-node error in diffusion Monte Carlo via range separation
• High-Level Coupled-Cluster Energetics by Monte Carlo Sampling and Moment Expansions: Further Details and Comparisons
• Improved Optimization for the Neural-network Quantum States and Tests on the Chromium Dimer
• Benchmark study of Nagaoka ferromagnetism by spin-adapted full configuration interaction quantum Monte Carlo
• Improved walker population control for full configuration interaction quantum Monte Carlo
• A full configuration interaction quantum Monte Carlo study of ScO, TiO and VO molecules
• Rationale for the Extrapolation Procedure in Selected Configuration Interaction
• Stochastic differential equation approach to understanding the population control bias in full configuration interaction quantum Monte Carlo
• Phaseless auxiliary field quantum Monte Carlo with projector-augmented wave method for solids
• Multi-configurational nature of electron correlation within nitrogen vacancy centers in diamond
• Converging High-Level Coupled-Cluster Energetics via Adaptive Selection of Excitation Manifolds Driven by Moment Expansions
• Benchmarking the Semi-Stochastic CC(P;Q) Approach for Singlet-Triplet Gaps in Biradicals
• Accelerating convergence to the thermodynamic limit with twist angle selection applied to methods beyond many-body perturbation theory
• Fixed and partial-node approximations in Slater determinant space for molecules
• Full Configuration Interaction Excited-State Energies in Large Active Spaces from Subspace Iteration with Repeated Random Sparsification
• 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
• Electronic specific heat capacities and entropies from density matrix quantum Monte Carlo using Gaussian process regression to find gradients of noisy data
• Rapidly convergent quantum Monte Carlo using a Chebyshev projector
• The Singlet-Triplet Gap of Cyclobutadiene: The CIPSI-Driven CC($P$;$Q$) Study
• High-Level Coupled-Cluster Energetics by Merging Moment Expansions with Selected Configuration Interaction
• Graphic characterization and clustering configuration descriptors of determinant space for molecules
• Addressing Strong Correlation by Approximate Coupled-Pair Methods with Active-Space and Full Treatments of Three-Body Clusters
• Improved Fast Randomized Iteration Approach to Full Configuration Interaction