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

Combining the Transcorrelated Method with Full Configuration Interaction Quantum Monte Carlo: Application to the Homogeneous Electron Gas.

Hongjun Luo, Ali Alavi

Published Dec 20, 2017• 0 citations
arXivHTMLPDF

Abstract

We suggest an efficient method to resolve electronic cusps in electronic structure calculations by using an effective transcorrelated Hamiltonian. This effective Hamiltonian takes a simple form for plane wave bases, containing up to two-body operators only, and its use incurs almost no additional computational overhead compared to that of the original Hamiltonian. We apply this method in combination with the full configuration interaction quantum Monte Carlo (FCIQMC) method to the homogeneous electron gas. As a projection technique, the non-Hermitian nature of the transcorrelated Hamiltonian does not cause complications or numerical difficulties for FCIQMC. The rate of convergence of the total energy to the complete basis set limit is improved from ${\cal O}(M^{-1})$ to ${\cal O}\left({M^{-5/3}}\right)$, where $M$ is the total number of orbital basis functions.

DOI
10.17863/cam.32443

Related papers

Combining the Transcorrelated Method with Full Configuration Interaction Quantum Monte Carlo: Application to the Homogeneous Electron Gas

Hongjun Luo, Ali Alavi

We suggest an efficient method to resolve electronic cusps in electronic structure calculations through the use of an effective transcorrelated Hamiltonian. This effective Hamiltonian takes a simple form for plane wave bases, containing up to two-body operators only, and its use incurs almost no additional computational overhead compared to that of the original Hamiltonian. We apply this method in combination with the full configuration interaction quantum Monte Carlo (FCIQMC) method to the homogeneous electron gas. As a projection technique, the non-Hermitian nature of the transcorrelated Hamiltonian does not cause complications or numerical difficulties for FCIQMC. The rate of convergence of the total energy to the complete basis set limit is improved from [Formula: see text] to [Formula: see text], where M is the total number of orbital basis functions.

Journal of Chemical Theory and Computation•Feb 12, 2018•83 citations

Full configuration interaction perspective on the homogeneous electron gas

James J. Shepherd, George H. Booth, Andreas Grüneis +1 more

Highly accurate results for the homogeneous electron gas (HEG) have only been achieved to date within a diffusion Monte Carlo (DMC) framework. Here, we introduce a recently developed stochastic technique, full configuration interaction quantum Monte Carlo (FCIQMC), which samples the exact wave function expanded in plane-wave Slater determinants. Despite the introduction of a basis-set incompleteness error, we obtain a finite-basis energy, which is significantly and variationally lower than any previously published work for the 54-electron HEG at ${r}_{s}$ $=$ 0.5 a.u., in a Hilbert space of ${10}^{108}$ Slater determinants. At this value of ${r}_{s}$, as well as of 1.0 a.u., we remove the remaining basis-set incompleteness error by extrapolation, yielding results comparable to state-of-the-art DMC backflow energies. In doing so, we demonstrate that it is possible to yield highly accurate results with the FCIQMC method in periodic systems.

arXiv:1109.2635•Physical Review B•Feb 7, 2012•147 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

Accelerating the convergence of exact diagonalization with the transcorrelated method: Quantum gas in one dimension with contact interactions

Péter Jeszenszki, Hongjun Luo, Ali Alavi +1 more

Exact diagonalization expansions of Bose or Fermi gases with contact interactions converge very slowly due to a nonanalytic cusp in the wave function. Here we develop a transcorrelated approach where the cusp is treated exactly and folded into the many-body Hamiltonian with a similarity transformation that removes the leading-order singularity. The resulting transcorrelated Hamiltonian is not Hermitian but can be treated numerically with a standard projection approach. The smoothness of the wave function improves by at least one order and thus the convergence rate for the ground-state energy improves. By numerical investigation of a one-dimensional gas of spin-$\frac{1}{2}$ fermions we find the error in the transcorrelated energy to scale as ${M}^{\ensuremath{-}3}$ with a single-particle basis of $M$ plane waves compared to ${M}^{\ensuremath{-}1}$ for the expansion of the original Hamiltonian and ${M}^{\ensuremath{-}2}$ using conventional lattice renormalization.

arXiv:1806.11268•Physical Review A•Nov 26, 2018•27 citations

Investigation of the full configuration interaction quantum Monte Carlo method using homogeneous electron gas models

James J. Shepherd, George H. Booth, Ali Alavi

Using the homogeneous electron gas (HEG) as a model, we investigate the sources of error in the "initiator" adaptation to full configuration interaction quantum Monte Carlo (i-FCIQMC), with a view to accelerating convergence. In particular, we find that the fixed-shift phase, where the walker number is allowed to grow slowly, can be used to effectively assess stochastic and initiator error. Using this approach we provide simple explanations for the internal parameters of an i-FCIQMC simulation. We exploit the consistent basis sets and adjustable correlation strength of the HEG to analyze properties of the algorithm, and present finite basis benchmark energies for N = 14 over a range of densities 0.5 ≤ r(s) ≤ 5.0 a.u. A single-point extrapolation scheme is introduced to produce complete basis energies for 14, 38, and 54 electrons. It is empirically found that, in the weakly correlated regime, the computational cost scales linearly with the plane wave basis set size, which is justifiable on physical grounds. We expect the fixed-shift strategy to reduce the computational cost of many i-FCIQMC calculations of weakly correlated systems. In addition, we provide benchmarks for the electron gas, to be used by other quantum chemical methods in exploring periodic solid state systems.

arXiv:1201.4691•The Journal of Chemical Physics•Jun 22, 2012•110 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
We suggest an efficient method to resolve electronic cusps in electronic structure calculations by using an effective transcorrelated Hamiltonian. This effective Hamiltonian takes a simple form for plane wave bases, containing up to two-body operators only, and its use incurs almost no additional computational overhead compared to that of the original Hamiltonian. We apply this method in combination with the full configuration interaction quantum Monte Carlo (FCIQMC) method to the homogeneous electron gas. As a projection technique, the non-Hermitian nature of the transcorrelated Hamiltonian does not cause complications or numerical difficulties for FCIQMC. The rate of convergence of the total energy to the complete basis set limit is improved from O(M−1){\cal O}(M^{-1})O(M−1) to O(M−5/3){\cal O}\left({M^{-5/3}}\right)O(M−5/3), where MMM is the total number of orbital basis functions.
We suggest an efficient method to resolve electronic cusps in electronic structure calculations through the use of an effective transcorrelated Hamiltonian. This effective Hamiltonian takes a simple form for plane wave bases, containing up to two-body operators only, and its use incurs almost no additional computational overhead compared to that of the original Hamiltonian. We apply this method in combination with the full configuration interaction quantum Monte Carlo (FCIQMC) method to the homogeneous electron gas. As a projection technique, the non-Hermitian nature of the transcorrelated Hamiltonian does not cause complications or numerical difficulties for FCIQMC. The rate of convergence of the total energy to the complete basis set limit is improved from [Formula: see text] to [Formula: see text], where M is the total number of orbital basis functions.
Highly accurate results for the homogeneous electron gas (HEG) have only been achieved to date within a diffusion Monte Carlo (DMC) framework. Here, we introduce a recently developed stochastic technique, full configuration interaction quantum Monte Carlo (FCIQMC), which samples the exact wave function expanded in plane-wave Slater determinants. Despite the introduction of a basis-set incompleteness error, we obtain a finite-basis energy, which is significantly and variationally lower than any previously published work for the 54-electron HEG at rs{r}_{s}rs​ === 0.5 a.u., in a Hilbert space of 10108{10}^{108}10108 Slater determinants. At this value of rs{r}_{s}rs​, as well as of 1.0 a.u., we remove the remaining basis-set incompleteness error by extrapolation, yielding results comparable to state-of-the-art DMC backflow energies. In doing so, we demonstrate that it is possible to yield highly accurate results with the FCIQMC method in periodic systems.
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.
Exact diagonalization expansions of Bose or Fermi gases with contact interactions converge very slowly due to a nonanalytic cusp in the wave function. Here we develop a transcorrelated approach where the cusp is treated exactly and folded into the many-body Hamiltonian with a similarity transformation that removes the leading-order singularity. The resulting transcorrelated Hamiltonian is not Hermitian but can be treated numerically with a standard projection approach. The smoothness of the wave function improves by at least one order and thus the convergence rate for the ground-state energy improves. By numerical investigation of a one-dimensional gas of spin-12\frac{1}{2}21​ fermions we find the error in the transcorrelated energy to scale as ${M}^{\ensuremath{-}3}$ with a single-particle basis of MMM plane waves compared to ${M}^{\ensuremath{-}1}$ for the expansion of the original Hamiltonian and ${M}^{\ensuremath{-}2}$ using conventional lattice renormalization.
Using the homogeneous electron gas (HEG) as a model, we investigate the sources of error in the "initiator" adaptation to full configuration interaction quantum Monte Carlo (i-FCIQMC), with a view to accelerating convergence. In particular, we find that the fixed-shift phase, where the walker number is allowed to grow slowly, can be used to effectively assess stochastic and initiator error. Using this approach we provide simple explanations for the internal parameters of an i-FCIQMC simulation. We exploit the consistent basis sets and adjustable correlation strength of the HEG to analyze properties of the algorithm, and present finite basis benchmark energies for N = 14 over a range of densities 0.5 ≤ r(s) ≤ 5.0 a.u. A single-point extrapolation scheme is introduced to produce complete basis energies for 14, 38, and 54 electrons. It is empirically found that, in the weakly correlated regime, the computational cost scales linearly with the plane wave basis set size, which is justifiable on physical grounds. We expect the fixed-shift strategy to reduce the computational cost of many i-FCIQMC calculations of weakly correlated systems. In addition, we provide benchmarks for the electron gas, to be used by other quantum chemical methods in exploring periodic solid state systems.
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).