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:1209.4023

Emergence of Critical Phenomena in Full Configuration Interaction Quantum Monte Carlo

James J. Shepherd, Lauretta R. Schwarz, Robert E. Thomas, Booth, George. H., Daan Frenkel, Ali Alavi

Published Sep 18, 2012• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.1209.4023

Related papers

Sign problem in full configuration interaction quantum Monte Carlo: Linear and sublinear representation regimes for the exact wave function

James J. Shepherd, Gustavo E. Scuseria, James S. Spencer

We investigate the sign problem for full configuration interaction quantum Monte Carlo (FCIQMC), a stochastic algorithm for finding the ground-state solution of the Schr\"odinger equation with substantially reduced computational cost compared with exact diagonalization. We find $k$-space Hubbard models for which the solution is yielded with storage that grows sublinearly in the size of the many-body Hilbert space, in spite of using a wave function that is simply a linear combination of states. The FCIQMC algorithm is able to find this sublinear scaling regime without bias and with only a choice of the Hamiltonian basis. By means of a demonstration we solve for the energy of a 70-site half-filled system (with a space of ${10}^{38}$ determinants) in 250 core hours, substantially quicker than the $\ensuremath{\sim}{10}^{36}$ core hours that would be required by exact diagonalization. This is the largest space that has been sampled in an unbiased fashion. The challenge for the recently developed FCIQMC method is made clear: Expand the sublinear scaling regime while retaining exact-on-average accuracy. We comment upon the relationship between this and the scaling law previously observed in the initiator adaptation (i-FCIQMC). We argue that our results change the landscape for the development of FCIQMC and related methods.

arXiv:1407.4800•Physical Review B•Oct 27, 2014•35 citations

The sign problem and population dynamics in the full configuration interaction quantum Monte Carlo method

J. S. Spencer, Nick S. Blunt, W. M. C. Foulkes

The recently proposed full configuration interaction quantum Monte Carlo method allows access to essentially exact ground-state energies of systems of interacting fermions substantially larger than previously tractable without knowledge of the nodal structure of the ground-state wave function. We investigate the nature of the sign problem in this method and how its severity depends on the system studied. We explain how cancellation of the positive and negative particles sampling the wave function ensures convergence to a stochastic representation of the many-fermion ground state and accounts for the characteristic population dynamics observed in simulations.

arXiv:1110.5479•The Journal of Chemical Physics•Feb 6, 2012•108 citations

The effect of quantization on the full configuration interaction quantum Monte Carlo sign problem

Michael Kolodrubetz, James S. Spencer, Bryan K. Clark +1 more

The sign problem in full configuration interaction quantum Monte Carlo (FCIQMC) without annihilation can be understood as an instability of the psi-particle population to the ground state of the matrix obtained by making all off-diagonal elements of the Hamiltonian negative. Such a matrix, and hence the sign problem, is basis dependent. In this paper, we discuss the properties of a physically important basis choice: first versus second quantization. For a given choice of single-particle orbitals, we identify the conditions under which the fermion sign problem in the second quantized basis of antisymmetric Slater determinants is identical to the sign problem in the first quantized basis of unsymmetrized Hartree products. We also show that, when the two differ, the fermion sign problem is always less severe in the second quantized basis. This supports the idea that FCIQMC, even in the absence of annihilation, improves the sign problem relative to first quantized methods. Finally, we point out some theoretically interesting classes of Hamiltonians where first and second quantized sign problems differ, and others where they do not.

arXiv:1209.3044•The Journal of Chemical Physics•Jan 11, 2013•30 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

Full configuration interaction quantum Monte Carlo and its application to the lattice Bose polaron at quantum criticality

Matija Čufar

This thesis provides an investigation of the application of the full configuration interaction quantum Monte Carlo (FCIQMC) algorithm to ultracold atom problems. This includes working on the algorithm, understanding the sign problem, and applying it to the study of strongly-correlated lattice polaron physics. The thesis conceptually consists of two parts. In the first, we discuss the FCIQMC algorithm, which is a modern variant of projector Monte Carlo and can be used for calculating the ground states or low-lying excited states of quantum Hamiltonians. We present our implementation of the algorithm in the Julia programming language and discuss the details of the algorithm and data structures employed. Next, we turn to a treatment of the FCIQMC sign problem, which arises when the number of walkers employed by FCIQMC is insufficient to capture the sign structure of the ground state eigenvector of the Hamiltonian. We present a new scalar model that accurately reconstructs the dynamic behaviour of FCIQMC and shows us that as the number of walkers increases, the walker population switches from a sign-incoherent regime to a sign-coherent one in way that is reminiscent of a continuous phase transition. In the second part, we focus on the application of FCIQMC to the Bose-Hubbard model (BHM) with an impurity near quantum criticality. We use FCIQMC to guide and validate a diagrammatic theory based on the quantum Gutzwiller method, which describes the polaron by treating interactions between the impurity and the elementary excitations in the BHM. This theory shows that at strong interactions, the polaron reflects the state of the host system and as such, polarons could be used as quantum sensors in experiments. Lastly, we show numerical evidence that the polaron energy --- the energy shift caused by introducing an impurity to the BHM --- is scale invariant near the phase transition and that it exhibits critical phenomena not yet explained by theory. Furthermore, we show that the polaron energy is linked to the impurity-boson density-density correlation function and that it exhibits a flattening near the phase transition. This flattening indicates that in the thermodynamic limit, the length scales associated with the polaron diverge.

•Jan 1, 2026•0 citations

Semi-stochastic full configuration interaction quantum Monte Carlo: Developments and application

Nick S. Blunt, Simon D. Smart, J. A. F. Kersten +3 more

We expand upon the recent semi-stochastic adaptation to full configuration interaction quantum Monte Carlo (FCIQMC). We present an alternate method for generating the deterministic space without a priori knowledge of the wave function and present stochastic efficiencies for a variety of both molecular and lattice systems. The algorithmic details of an efficient semi-stochastic implementation are presented, with particular consideration given to the effect that the adaptation has on parallel performance in FCIQMC. We further demonstrate the benefit for calculation of reduced density matrices in FCIQMC through replica sampling, where the semi-stochastic adaptation seems to have even larger efficiency gains. We then combine these ideas to produce explicitly correlated corrected FCIQMC energies for the beryllium dimer, for which stochastic errors on the order of wavenumber accuracy are achievable.

arXiv:1502.04847•The Journal of Chemical Physics•May 14, 2015•100 citations
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.
We investigate the sign problem for full configuration interaction quantum Monte Carlo (FCIQMC), a stochastic algorithm for finding the ground-state solution of the Schr\"odinger equation with substantially reduced computational cost compared with exact diagonalization. We find kkk-space Hubbard models for which the solution is yielded with storage that grows sublinearly in the size of the many-body Hilbert space, in spite of using a wave function that is simply a linear combination of states. The FCIQMC algorithm is able to find this sublinear scaling regime without bias and with only a choice of the Hamiltonian basis. By means of a demonstration we solve for the energy of a 70-site half-filled system (with a space of 1038{10}^{38}1038 determinants) in 250 core hours, substantially quicker than the $\ensuremath{\sim}{10}^{36}$ core hours that would be required by exact diagonalization. This is the largest space that has been sampled in an unbiased fashion. The challenge for the recently developed FCIQMC method is made clear: Expand the sublinear scaling regime while retaining exact-on-average accuracy. We comment upon the relationship between this and the scaling law previously observed in the initiator adaptation (i-FCIQMC). We argue that our results change the landscape for the development of FCIQMC and related methods.
The recently proposed full configuration interaction quantum Monte Carlo method allows access to essentially exact ground-state energies of systems of interacting fermions substantially larger than previously tractable without knowledge of the nodal structure of the ground-state wave function. We investigate the nature of the sign problem in this method and how its severity depends on the system studied. We explain how cancellation of the positive and negative particles sampling the wave function ensures convergence to a stochastic representation of the many-fermion ground state and accounts for the characteristic population dynamics observed in simulations.
The sign problem in full configuration interaction quantum Monte Carlo (FCIQMC) without annihilation can be understood as an instability of the psi-particle population to the ground state of the matrix obtained by making all off-diagonal elements of the Hamiltonian negative. Such a matrix, and hence the sign problem, is basis dependent. In this paper, we discuss the properties of a physically important basis choice: first versus second quantization. For a given choice of single-particle orbitals, we identify the conditions under which the fermion sign problem in the second quantized basis of antisymmetric Slater determinants is identical to the sign problem in the first quantized basis of unsymmetrized Hartree products. We also show that, when the two differ, the fermion sign problem is always less severe in the second quantized basis. This supports the idea that FCIQMC, even in the absence of annihilation, improves the sign problem relative to first quantized methods. Finally, we point out some theoretically interesting classes of Hamiltonians where first and second quantized sign problems differ, and others where they do not.
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.
This thesis provides an investigation of the application of the full configuration interaction quantum Monte Carlo (FCIQMC) algorithm to ultracold atom problems. This includes working on the algorithm, understanding the sign problem, and applying it to the study of strongly-correlated lattice polaron physics. The thesis conceptually consists of two parts. In the first, we discuss the FCIQMC algorithm, which is a modern variant of projector Monte Carlo and can be used for calculating the ground states or low-lying excited states of quantum Hamiltonians. We present our implementation of the algorithm in the Julia programming language and discuss the details of the algorithm and data structures employed. Next, we turn to a treatment of the FCIQMC sign problem, which arises when the number of walkers employed by FCIQMC is insufficient to capture the sign structure of the ground state eigenvector of the Hamiltonian. We present a new scalar model that accurately reconstructs the dynamic behaviour of FCIQMC and shows us that as the number of walkers increases, the walker population switches from a sign-incoherent regime to a sign-coherent one in way that is reminiscent of a continuous phase transition. In the second part, we focus on the application of FCIQMC to the Bose-Hubbard model (BHM) with an impurity near quantum criticality. We use FCIQMC to guide and validate a diagrammatic theory based on the quantum Gutzwiller method, which describes the polaron by treating interactions between the impurity and the elementary excitations in the BHM. This theory shows that at strong interactions, the polaron reflects the state of the host system and as such, polarons could be used as quantum sensors in experiments. Lastly, we show numerical evidence that the polaron energy --- the energy shift caused by introducing an impurity to the BHM --- is scale invariant near the phase transition and that it exhibits critical phenomena not yet explained by theory. Furthermore, we show that the polaron energy is linked to the impurity-boson density-density correlation function and that it exhibits a flattening near the phase transition. This flattening indicates that in the thermodynamic limit, the length scales associated with the polaron diverge.
We expand upon the recent semi-stochastic adaptation to full configuration interaction quantum Monte Carlo (FCIQMC). We present an alternate method for generating the deterministic space without a priori knowledge of the wave function and present stochastic efficiencies for a variety of both molecular and lattice systems. The algorithmic details of an efficient semi-stochastic implementation are presented, with particular consideration given to the effect that the adaptation has on parallel performance in FCIQMC. We further demonstrate the benefit for calculation of reduced density matrices in FCIQMC through replica sampling, where the semi-stochastic adaptation seems to have even larger efficiency gains. We then combine these ideas to produce explicitly correlated corrected FCIQMC energies for the beryllium dimer, for which stochastic errors on the order of wavenumber accuracy are achievable.