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:1204.1490 · Physical Review B

Partial node configuration-interaction Monte Carlo as applied to the Fermi polaron

Michael Kolodrubetz, Bryan K. Clark

Published Aug 6, 2012• 21 citations
arXivHTMLPDF

Abstract

Finding the ground state of a fermionic Hamiltonian using quantum Monte Carlo (QMC) is a very difficult problem, due to the Fermi sign problem. While still scaling exponentially, full configuration-interaction Monte Carlo (FCIQMC) mitigates some of the exponential variance by allowing annihilation of noise---whenever two walkers arrive at the same configuration with opposite signs, they are removed from the simulation. While FCIQMC has been quite successful for quantum chemistry problems, its application to problems in condensed systems has been limited. In this paper, we apply the FCIQMC algorithm to the Fermi polaron problem, which provides an ideal test bed for improving the algorithm, although since we restrict the number of many-body excitations, our algorithm is more precisely a QMC implementation of the configuration-interaction method (CIQMC). In its simplest form, CIQMC is unstable for even a fairly small system sizes. However, with a series of algorithmic improvements, we are able to significantly increase its effectiveness. We modify fixed node QMC to work in these systems, introduce a well chosen importance sampled trial wave function, a partial node approximation, and a variant of release node. Finally, we develop a way to perform CIQMC directly in the thermodynamic limit.

DOI
10.1103/physrevb.86.075109

Related papers

Approaching Strongly Correlated Systems Using Partial Node FCIQMC

Bryan K. Clark, Michael Kolodrubetz

problem, due to the Fermion sign problem. While still scaling exponentially, full configuration-interaction Monte Carlo (FCI-QMC) [1] mitigates some of the exponential variance by allowing annihilation of noise ‐ whenever two walkers arrive at the same configuration with opposite signs, they are removed from the simulation. While FCI-QMC has been quite successful for quantum chemistry problems, it has had less success with strongly correlated systems; for example, it is unstable for even small systems sizes when applied to the one of the simplest many-body systems: the Fermi-Polaron. In this talk, we explore a series of algorithmic improvements which significantly increase its eectiveness. These improvements include modifying fixed node to work in this method, developing a partial node approach which allows for an extrapolation to the exact answer and developing a variant of release node. Additionally, we find a way to work directly in the thermodynamic limit for some subset of important systems. We apply these new approaches to the Fermi-Polaron finding that it accurately reproduces known analytic and numerical results going beyond them in certain cases [2].

•Jan 1, 2012•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

APPLYING QUANTUM MONTE CARLO TO THE ELECTRONIC STRUCTURE PROBLEM

Andrew Powell, Richard Dawes

Two distinct types of Quantum Monte Carlo (QMC) calculations are applied to electronic structure problems such as calculating potential energy curves and producing benchmark values for reaction barriers. First, Variational and Diffusion Monte Carlo (VMC and DMC) methods using a trial wavefunction subject to the fixed node approximation were tested using the CASINO code.[1] Next, Full Configuration Interaction Quantum Monte Carlo (FCIQMC), along with its initiator extension (i-FCIQMC) were tested using the NECI code.[2] FCIQMC seeks the FCI energy for a specific basis set. At a reduced cost, the efficient i-FCIQMC method can be applied to systems in which the standard FCIQMC approach proves to be too costly. Since all of these methods are statistical approaches, uncertainties (error-bars) are introduced for each calculated energy. This study tests the performance of the methods relative to traditional quantum chemistry for some benchmark systems. References: [1] R. J. Needs et al., J. Phys.: Condensed Matter 22, 023201 (2010). [2] G. H. Booth et al., J. Chem. Phys. 131, 054106 (2009).

International Symposium on Molecular Spectroscopy (ISMS)•Jun 21, 2016•0 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

Electronic structure quantum Monte Carlo

Michal Bajdich, Luboš Mitáš

Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.

arXiv:1008.2369•Acta Physica Slovaca Reviews and Tutorials•Apr 1, 2009•33 citations

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
Finding the ground state of a fermionic Hamiltonian using quantum Monte Carlo (QMC) is a very difficult problem, due to the Fermi sign problem. While still scaling exponentially, full configuration-interaction Monte Carlo (FCIQMC) mitigates some of the exponential variance by allowing annihilation of noise---whenever two walkers arrive at the same configuration with opposite signs, they are removed from the simulation. While FCIQMC has been quite successful for quantum chemistry problems, its application to problems in condensed systems has been limited. In this paper, we apply the FCIQMC algorithm to the Fermi polaron problem, which provides an ideal test bed for improving the algorithm, although since we restrict the number of many-body excitations, our algorithm is more precisely a QMC implementation of the configuration-interaction method (CIQMC). In its simplest form, CIQMC is unstable for even a fairly small system sizes. However, with a series of algorithmic improvements, we are able to significantly increase its effectiveness. We modify fixed node QMC to work in these systems, introduce a well chosen importance sampled trial wave function, a partial node approximation, and a variant of release node. Finally, we develop a way to perform CIQMC directly in the thermodynamic limit.
problem, due to the Fermion sign problem. While still scaling exponentially, full configuration-interaction Monte Carlo (FCI-QMC) [1] mitigates some of the exponential variance by allowing annihilation of noise ‐ whenever two walkers arrive at the same configuration with opposite signs, they are removed from the simulation. While FCI-QMC has been quite successful for quantum chemistry problems, it has had less success with strongly correlated systems; for example, it is unstable for even small systems sizes when applied to the one of the simplest many-body systems: the Fermi-Polaron. In this talk, we explore a series of algorithmic improvements which significantly increase its eectiveness. These improvements include modifying fixed node to work in this method, developing a partial node approach which allows for an extrapolation to the exact answer and developing a variant of release node. Additionally, we find a way to work directly in the thermodynamic limit for some subset of important systems. We apply these new approaches to the Fermi-Polaron finding that it accurately reproduces known analytic and numerical results going beyond them in certain cases [2].
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.
Two distinct types of Quantum Monte Carlo (QMC) calculations are applied to electronic structure problems such as calculating potential energy curves and producing benchmark values for reaction barriers. First, Variational and Diffusion Monte Carlo (VMC and DMC) methods using a trial wavefunction subject to the fixed node approximation were tested using the CASINO code.[1] Next, Full Configuration Interaction Quantum Monte Carlo (FCIQMC), along with its initiator extension (i-FCIQMC) were tested using the NECI code.[2] FCIQMC seeks the FCI energy for a specific basis set. At a reduced cost, the efficient i-FCIQMC method can be applied to systems in which the standard FCIQMC approach proves to be too costly. Since all of these methods are statistical approaches, uncertainties (error-bars) are introduced for each calculated energy. This study tests the performance of the methods relative to traditional quantum chemistry for some benchmark systems. References: [1] R. J. Needs et al., J. Phys.: Condensed Matter 22, 023201 (2010). [2] G. H. Booth et al., J. Chem. Phys. 131, 054106 (2009).
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.
Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.
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.