Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1908.04158 · The Journal of Chemical Physics

A hybrid approach to extending selected configuration interaction and full configuration interaction quantum Monte Carlo

Nick S. Blunt

Published Nov 4, 2019• 14 citations
arXivHTMLPDF

Abstract

We present an approach to combining selected configuration interaction (SCI) and initiator full configuration interaction quantum Monte Carlo (i-FCIQMC). In the current i-FCIQMC scheme, the space of initiators is chosen dynamically by a population threshold. Here, we instead choose initiators as the selected space (V) from a prior SCI calculation, allowing substantially larger initiator spaces for a given walker population. While SCI+PT2 adds a perturbative correction in the first-order interacting space beyond V, the approach presented here allows a variational calculation in the same space and a perturbative correction in the second-order interacting space. The use of a fixed initiator space reintroduces population plateaus into FCIQMC, but it is shown that the plateau height is typically only a small multiple of the size of V. Thus, for a comparable fundamental memory cost to SCI+PT2, a substantially larger space can be sampled. The resulting method can be seen as a complementary approach to SCI+PT2, which is more accurate but slower for a common selected/initiator space. More generally, our results show that approaches exist to significantly improve initiator energies in i-FCIQMC while still ameliorating the fermion sign problem relative to the original FCIQMC method.

DOI
10.1063/1.5123146

Related papers

Improved walker population control for full configuration interaction quantum Monte Carlo

Mingrui Yang, Elke Pahl, Joachim Brand

Full configuration interaction quantum Monte Carlo (FCIQMC) is a stochastic approach for finding the ground state of a quantum many-body Hamiltonian. It is based on the dynamical evolution of a walker population in Hilbert space, which samples the ground state configuration vector over many iterations. Here, we present a modification of the original protocol for walker population control of Booth et al. [J. Chem. Phys. 131, 054106 (2009)] in order to achieve equilibration at a pre-defined average walker number and to avoid walker number overshoots. The dynamics of the walker population is described by a noisy damped harmonic oscillator and controlled by two parameters responsible for damping and forcing, respectively, for which reasonable values are suggested. We further introduce a population growth witness that can be used to detect annihilation plateaus related to overcoming the FCIQMC sign problem. Features of the new population control procedure such as precise walker number control and fast equilibration are demonstrated. The standard error of the shift estimator for the ground state energy as well as the population control bias is found to be unaffected by the population control procedure or its parameters. The improved control of the walker number, and thereby memory consumption, is a desirable feature required for automating FCIQMC calculations and requires minimal modifications to the existing code.

arXiv:2008.01927•The Journal of Chemical Physics•Nov 2, 2020•20 citations

An explicitly correlated approach to basis set incompleteness in full configuration interaction quantum Monte Carlo

George H. Booth, Deidre M. Cleland, Ali Alavi +1 more

By performing a stochastic dynamic in a space of Slater determinants, the full configuration interaction quantum Monte Carlo (FCIQMC) method has been able to obtain energies which are essentially free from systematic error to the basis set correlation energy, within small and systematically improvable error bars. However, the weakly exponential scaling with basis size makes converging the energy with respect to basis set costly and in larger systems, impossible. To ameliorate these basis set issues, here we use perturbation theory to couple the FCIQMC wavefunction to an explicitly correlated strongly orthogonal basis of geminals, following the [2](R12) approach of Valeev et al. The required one- and two-particle density matrices are computed on-the-fly during the FCIQMC dynamic, using a sampling procedure which incurs relatively little additional computation expense. The F12 energy corrections are shown to converge rapidly as a function of sampling, both in imaginary time and number of walkers. Our pilot calculations on the binding curve for the carbon dimer, which exhibits strong correlation effects as well as substantial basis set dependence, demonstrate that the accuracy of the FCIQMC-F12 method surpasses that of all previous FCIQMC calculations, and that the F12 correction improves results equivalent to increasing the quality of the one-electron basis by two cardinal numbers.

arXiv:1208.0980•The Journal of Chemical Physics•Oct 28, 2012•60 citations

Preconditioning and Perturbative Estimators in Full Configuration Interaction Quantum Monte Carlo

Nick S. Blunt, Alex J. W. Thom, Charles J. C. Scott

We propose the use of preconditioning in FCIQMC which, in combination with perturbative estimators, greatly increases the efficiency of the algorithm. The use of preconditioning allows a time step close to unity to be used (without time-step errors), provided that multiple spawning attempts are made per walker. We show that this approach substantially reduces statistical noise on perturbative corrections to initiator error, which improve the accuracy of FCIQMC but which can suffer from significant noise in the original scheme. Therefore, the use of preconditioning and perturbatively corrected estimators in combination leads to a significantly more efficient algorithm. In addition, a simpler approach to sampling variational and perturbative estimators in FCIQMC is presented, which also allows the variance of the energy to be calculated. These developments are investigated and applied to benzene (30e, 108o), an example where accurate treatment is not possible with the original method.

arXiv:1901.06348•Journal of Chemical Theory and Computation•May 3, 2019•27 citations

Full configuration interaction quantum Monte Carlo for accurate $\textit{ab initio}$ nuclear structure calculations: algorithms and calculation details

Rongzhe Hu, Furong Xu, Baishan Hu +1 more

Full configuration interaction quantum Monte Carlo (FCIQMC) is a stochastic many-body solver that has been widely applied to electronic, molecular, and condensed-matter systems. In this work we apply FCIQMC to $\textit{ab initio}$ nuclear structure calculations using interactions derived from chiral effective field theory. We describe the algorithm in detail, including imaginary-time propagation, excitation generation, estimator choices, the initiator approximation with adaptive shift correction, and reduced-density-matrix (RDM) sampling. Benchmark calculations in small model spaces, where deterministic full configuration interaction (FCI) results are available, validate the stochastic calculation of energies, radii, and RDM-based pure estimators. For large model spaces, we analyze the residual finite-walker bias through systematic walker-number convergence and infinite-walker extrapolations. We also demonstrate that FCIQMC can be extended beyond ground-state calculations by computing the low-lying spectrum of $^6$Li.

arXiv:2607.05525•Jul 6, 2026•0 citations

Beyond Walkers in Stochastic Quantum Chemistry: Reducing Error Using Fast Randomized Iteration

Samuel M. Greene, Robert J. Webber, Jonathan Weare +1 more

We introduce a family of methods for the full configuration interaction problem in quantum chemistry, based on the fast randomized iteration (FRI) framework [Lim, L.-H.; Weare, J. SIAM Rev. 2017, 59, 547; DOI: 10.1137/15M1040827 ]. These methods, which we term “FCI-FRI”, stochastically impose sparsity during iterations of the power method and can be viewed as a generalization of full configuration interaction quantum Monte Carlo (FCIQMC) without walkers. In addition to the multinomial scheme commonly used to sample excitations in FCIQMC, we present a systematic scheme where excitations are not sampled independently. Performing ground-state calculations on five small molecules at fixed cost, we find that the systematic FCI-FRI scheme is 11–45 times more statistically efficient than the multinomial FCI-FRI scheme, which is in turn 1.4–178 times more statistically efficient than the original FCIQMC algorithm.

arXiv:1905.00995•Journal of Chemical Theory and Computation•Aug 7, 2019•5 citations

Emergence of Critical Phenomena in Full Configuration Interaction Quantum Monte Carlo

James J. Shepherd, Lauretta R. Schwarz, Robert E. Thomas +3 more

There has been recent literature discussion on the origin and severity of the `sign problem' in full configuration interaction quantum Monte Carlo (FCIQMC) and its `initiator' adaptation (i-FCIQMC), methods of interest and potential because they allow for exact (FCI) ground-state solutions to be obtained often at a much reduced computational cost. In this study we aim to use a simple order parameter, describing the `sign structure' of the stochastic wavefunction representation, to empirically characterise the fundamentally different collective behaviour of the walker population in both methods.

arXiv:1209.4023•Sep 18, 2012•0 citations
We present an approach to combining selected configuration interaction (SCI) and initiator full configuration interaction quantum Monte Carlo (i-FCIQMC). In the current i-FCIQMC scheme, the space of initiators is chosen dynamically by a population threshold. Here, we instead choose initiators as the selected space (V) from a prior SCI calculation, allowing substantially larger initiator spaces for a given walker population. While SCI+PT2 adds a perturbative correction in the first-order interacting space beyond V, the approach presented here allows a variational calculation in the same space and a perturbative correction in the second-order interacting space. The use of a fixed initiator space reintroduces population plateaus into FCIQMC, but it is shown that the plateau height is typically only a small multiple of the size of V. Thus, for a comparable fundamental memory cost to SCI+PT2, a substantially larger space can be sampled. The resulting method can be seen as a complementary approach to SCI+PT2, which is more accurate but slower for a common selected/initiator space. More generally, our results show that approaches exist to significantly improve initiator energies in i-FCIQMC while still ameliorating the fermion sign problem relative to the original FCIQMC method.
Full configuration interaction quantum Monte Carlo (FCIQMC) is a stochastic approach for finding the ground state of a quantum many-body Hamiltonian. It is based on the dynamical evolution of a walker population in Hilbert space, which samples the ground state configuration vector over many iterations. Here, we present a modification of the original protocol for walker population control of Booth et al. [J. Chem. Phys. 131, 054106 (2009)] in order to achieve equilibration at a pre-defined average walker number and to avoid walker number overshoots. The dynamics of the walker population is described by a noisy damped harmonic oscillator and controlled by two parameters responsible for damping and forcing, respectively, for which reasonable values are suggested. We further introduce a population growth witness that can be used to detect annihilation plateaus related to overcoming the FCIQMC sign problem. Features of the new population control procedure such as precise walker number control and fast equilibration are demonstrated. The standard error of the shift estimator for the ground state energy as well as the population control bias is found to be unaffected by the population control procedure or its parameters. The improved control of the walker number, and thereby memory consumption, is a desirable feature required for automating FCIQMC calculations and requires minimal modifications to the existing code.
By performing a stochastic dynamic in a space of Slater determinants, the full configuration interaction quantum Monte Carlo (FCIQMC) method has been able to obtain energies which are essentially free from systematic error to the basis set correlation energy, within small and systematically improvable error bars. However, the weakly exponential scaling with basis size makes converging the energy with respect to basis set costly and in larger systems, impossible. To ameliorate these basis set issues, here we use perturbation theory to couple the FCIQMC wavefunction to an explicitly correlated strongly orthogonal basis of geminals, following the [2](R12) approach of Valeev et al. The required one- and two-particle density matrices are computed on-the-fly during the FCIQMC dynamic, using a sampling procedure which incurs relatively little additional computation expense. The F12 energy corrections are shown to converge rapidly as a function of sampling, both in imaginary time and number of walkers. Our pilot calculations on the binding curve for the carbon dimer, which exhibits strong correlation effects as well as substantial basis set dependence, demonstrate that the accuracy of the FCIQMC-F12 method surpasses that of all previous FCIQMC calculations, and that the F12 correction improves results equivalent to increasing the quality of the one-electron basis by two cardinal numbers.
We propose the use of preconditioning in FCIQMC which, in combination with perturbative estimators, greatly increases the efficiency of the algorithm. The use of preconditioning allows a time step close to unity to be used (without time-step errors), provided that multiple spawning attempts are made per walker. We show that this approach substantially reduces statistical noise on perturbative corrections to initiator error, which improve the accuracy of FCIQMC but which can suffer from significant noise in the original scheme. Therefore, the use of preconditioning and perturbatively corrected estimators in combination leads to a significantly more efficient algorithm. In addition, a simpler approach to sampling variational and perturbative estimators in FCIQMC is presented, which also allows the variance of the energy to be calculated. These developments are investigated and applied to benzene (30e, 108o), an example where accurate treatment is not possible with the original method.
Full configuration interaction quantum Monte Carlo for accurate ab initio\textit{ab initio}ab initio nuclear structure calculations: algorithms and calculation details
Full configuration interaction quantum Monte Carlo (FCIQMC) is a stochastic many-body solver that has been widely applied to electronic, molecular, and condensed-matter systems. In this work we apply FCIQMC to ab initio\textit{ab initio}ab initio nuclear structure calculations using interactions derived from chiral effective field theory. We describe the algorithm in detail, including imaginary-time propagation, excitation generation, estimator choices, the initiator approximation with adaptive shift correction, and reduced-density-matrix (RDM) sampling. Benchmark calculations in small model spaces, where deterministic full configuration interaction (FCI) results are available, validate the stochastic calculation of energies, radii, and RDM-based pure estimators. For large model spaces, we analyze the residual finite-walker bias through systematic walker-number convergence and infinite-walker extrapolations. We also demonstrate that FCIQMC can be extended beyond ground-state calculations by computing the low-lying spectrum of 6^66Li.
We introduce a family of methods for the full configuration interaction problem in quantum chemistry, based on the fast randomized iteration (FRI) framework [Lim, L.-H.; Weare, J. SIAM Rev. 2017, 59, 547; DOI: 10.1137/15M1040827 ]. These methods, which we term “FCI-FRI”, stochastically impose sparsity during iterations of the power method and can be viewed as a generalization of full configuration interaction quantum Monte Carlo (FCIQMC) without walkers. In addition to the multinomial scheme commonly used to sample excitations in FCIQMC, we present a systematic scheme where excitations are not sampled independently. Performing ground-state calculations on five small molecules at fixed cost, we find that the systematic FCI-FRI scheme is 11–45 times more statistically efficient than the multinomial FCI-FRI scheme, which is in turn 1.4–178 times more statistically efficient than the original FCIQMC algorithm.
There has been recent literature discussion on the origin and severity of the `sign problem' in full configuration interaction quantum Monte Carlo (FCIQMC) and its `initiator' adaptation (i-FCIQMC), methods of interest and potential because they allow for exact (FCI) ground-state solutions to be obtained often at a much reduced computational cost. In this study we aim to use a simple order parameter, describing the `sign structure' of the stochastic wavefunction representation, to empirically characterise the fundamentally different collective behaviour of the walker population in both methods.