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

X marks the spot: accurate energies from intersecting extrapolations of continuum quantum Monte Carlo data

Seyed Mohammadreza Hosseini, Ali Alavi, Pablo López Ríos

Published Mar 1, 2024• 0 citations
arXivHTMLPDF

Abstract

We explore the application of an extrapolative method that yields very accurate total and relative energies from variational and diffusion quantum Monte Carlo (VMC and DMC) results. For a trial wave function consisting of a small configuration interaction (CI) wave function obtained from full CI quantum Monte Carlo and reoptimized in the presence of a Jastrow factor and an optional backflow transformation, we find that the VMC and DMC energies are smooth functions of the sum of the squared coefficients of the initial CI wave function, and that quadratic extrapolations of the non-backflow VMC and backflow DMC energies intersect within uncertainty of the exact total energy. With adequate statistical treatment of quasi-random fluctuations, the extrapolate and intersect with polynomials of order two (XSPOT) method is shown to yield results in agreement with benchmark-quality total and relative energies for the C2, N2, CO2, and H2O molecules, as well as for the C2 molecule in its first electronic singlet excited state, using only small CI expansion sizes.

DOI
10.48550/arxiv.2403.00649

Related papers

X marks the spot: Accurate energies from intersecting extrapolations of continuum quantum Monte Carlo data

Seyed Mohammadreza Hosseini, Ali Alavi, Pablo López Ríos

We explore the application of an extrapolative method that yields very accurate total and relative energies from variational and diffusion quantum Monte Carlo (VMC and DMC) results. For a trial wave function consisting of a small configuration interaction (CI) wave function obtained from full CI quantum Monte Carlo and reoptimized in the presence of a Jastrow factor and an optional backflow transformation, we find that the VMC and DMC energies are smooth functions of the sum of the squared coefficients of the initial CI wave function and that quadratic extrapolations of the non-backflow VMC and backflow DMC energies intersect within uncertainty of the exact total energy. With adequate statistical treatment of quasi-random fluctuations, the extrapolate and intersect with polynomials of order two method is shown to yield results in agreement with benchmark-quality total and relative energies for the C2, N2, CO2, and H2O molecules, as well as for the C2 molecule in its first electronic singlet excited state, using only small CI expansion sizes.

The Journal of Chemical Physics•May 15, 2024•1 citation

Toward Compact Selected Configuration Interaction Wave Functions with Quantum Monte Carlo─A Case Study of C2

Jil Ludovicy, Robin Dahl, Arne Lüchow

is investigated using truncated CIPSI-Jastrow CSF wave functions with Hartree-Fock orbitals within the framework of variational and diffusion quantum Monte Carlo. The truncation is performed based on the absolute value of the CI coefficients, and the Jastrow, molecular orbitals, and CI parameters are either partially or fully reoptimized with respect to the variational energy. Excellent absolute as well as bond dissociation energies are obtained at DMC level with very compact, fully optimized wave functions. By studying the expansions in more detail, we observe a change in the CI picture when reoptimizing the antisymmetric part of the CIPSI-Jastrow wave functions. Furthermore, we demonstrate that a decrease in the VMC energy as well as an improvement of the nodal surface quality can be achieved─with the same expansion size─if the CSFs are selected in the presence of a Jastrow correlation function, laying the foundation for a Jastrow selected CI scheme with quantum Monte Carlo.

Journal of Chemical Theory and Computation•May 2, 2023•7 citations

Variational Quantum Monte Carlo with Inclusion of Orbital Correlations

Shigenori Tanaka

Variational quantum Monte Carlo (VMC) method applied to interacting electron systems has played a pivotal role for the studies of electronic states of atoms, molecules and condensed matters. First, the many-body wave functions optimized by the VMC method provide good guiding functions in more accurate, diffusion quantum Monte Carlo (DMC) calculations for identical systems. Second, the VMC calculations themselves can give a good upper bound for the exact ground-state energy of many-electron systems when an appropriate form for the variational wave function is employed. In addition to a conventional form for the variational wave function that is expressed by the product of the Jastrow factor and the single Slater determinant, recent investigations have found that the expression in terms of many Slater determinants can provide a very accurate description of electronic states, especially for intrinsically multi-reference systems, which corresponds to the configuration interaction (CI) expansion method in conventional quantum chemistry. Although the CI method needs an extremely large number of expansion coefficients as variational parameters, the number of parameters can be reduced efficiently by using an idea of (molecular) orbital pair correlations. The purpose of the present paper is to examine the feasibility and usefulness of this idea within the framework of VMC method with the Jastrow–Slater wave function, where the ground state of a water molecule is investigated as a test case. In this study we rely on a variational wave function composed of a product of the Jastrow and Slater parts as 1⁄4 J S; ð1Þ

Journal of the Physical Society of Japan•Jun 10, 2013•3 citations

Chemical accuracy from quantum Monte Carlo for the benzene dimer

Sam Azadi, R. E. Cohen

We report an accurate study of interactions between benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory using different van der Waals functionals. In our quantum Monte Carlo (QMC) calculations, we use accurate correlated trial wave functions including three-body Jastrow factors and backflow transformations. We consider two benzene molecules in the parallel displaced geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the coupled-cluster theory through perturbative triplets/complete basis set limit is -2.65(2) kcal/mol [Miliordos et al., J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.

arXiv:1509.00995•The Journal of Chemical Physics•Sep 8, 2015•25 citations

Quantum Monte Carlo calculations of energy gaps from first principles

R. J. Hunt, Marcin Szyniszewski, Genki I. Prayogo +2 more

We review the use of continuum quantum Monte Carlo (QMC) methods for the calculation of energy gaps from first principles, and present a broad set of excited-state calculations carried out with the variational and fixed-node diffusion QMC methods on atoms, molecules, and solids. We propose a finite-size-error correction scheme for bulk energy gaps calculated in finite cells subject to periodic boundary conditions. We show that finite-size effects are qualitatively different in two-dimensional materials, demonstrating the effect in a QMC calculation of the band gap and exciton binding energy of monolayer phosphorene. We investigate the fixed-node errors in diffusion Monte Carlo gaps evaluated with Slater-Jastrow trial wave functions by examining the effects of backflow transformations, and also by considering the formation of restricted multideterminant expansions for excited-state wave functions. For several molecules, we examine the importance of structural relaxation in the excited state in determining excited-state energies. We study the feasibility of using variational Monte Carlo with backflow correlations to obtain accurate excited-state energies at reduced computational cost, finding that this approach can be valid. We find that diffusion Monte Carlo gap calculations can be performed with much larger time steps than are typically required to converge the total energy, at significantly diminished computational expense, but that in order to alleviate fixed-node errors in calculations on solids the inclusion of backflow correlations is sometimes necessary.

arXiv:1806.04750•Physical review. B./Physical review. B•Aug 13, 2018•50 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 explore the application of an extrapolative method that yields very accurate total and relative energies from variational and diffusion quantum Monte Carlo (VMC and DMC) results. For a trial wave function consisting of a small configuration interaction (CI) wave function obtained from full CI quantum Monte Carlo and reoptimized in the presence of a Jastrow factor and an optional backflow transformation, we find that the VMC and DMC energies are smooth functions of the sum of the squared coefficients of the initial CI wave function, and that quadratic extrapolations of the non-backflow VMC and backflow DMC energies intersect within uncertainty of the exact total energy. With adequate statistical treatment of quasi-random fluctuations, the extrapolate and intersect with polynomials of order two (XSPOT) method is shown to yield results in agreement with benchmark-quality total and relative energies for the C2, N2, CO2, and H2O molecules, as well as for the C2 molecule in its first electronic singlet excited state, using only small CI expansion sizes.
We explore the application of an extrapolative method that yields very accurate total and relative energies from variational and diffusion quantum Monte Carlo (VMC and DMC) results. For a trial wave function consisting of a small configuration interaction (CI) wave function obtained from full CI quantum Monte Carlo and reoptimized in the presence of a Jastrow factor and an optional backflow transformation, we find that the VMC and DMC energies are smooth functions of the sum of the squared coefficients of the initial CI wave function and that quadratic extrapolations of the non-backflow VMC and backflow DMC energies intersect within uncertainty of the exact total energy. With adequate statistical treatment of quasi-random fluctuations, the extrapolate and intersect with polynomials of order two method is shown to yield results in agreement with benchmark-quality total and relative energies for the C2, N2, CO2, and H2O molecules, as well as for the C2 molecule in its first electronic singlet excited state, using only small CI expansion sizes.
is investigated using truncated CIPSI-Jastrow CSF wave functions with Hartree-Fock orbitals within the framework of variational and diffusion quantum Monte Carlo. The truncation is performed based on the absolute value of the CI coefficients, and the Jastrow, molecular orbitals, and CI parameters are either partially or fully reoptimized with respect to the variational energy. Excellent absolute as well as bond dissociation energies are obtained at DMC level with very compact, fully optimized wave functions. By studying the expansions in more detail, we observe a change in the CI picture when reoptimizing the antisymmetric part of the CIPSI-Jastrow wave functions. Furthermore, we demonstrate that a decrease in the VMC energy as well as an improvement of the nodal surface quality can be achieved─with the same expansion size─if the CSFs are selected in the presence of a Jastrow correlation function, laying the foundation for a Jastrow selected CI scheme with quantum Monte Carlo.
Variational quantum Monte Carlo (VMC) method applied to interacting electron systems has played a pivotal role for the studies of electronic states of atoms, molecules and condensed matters. First, the many-body wave functions optimized by the VMC method provide good guiding functions in more accurate, diffusion quantum Monte Carlo (DMC) calculations for identical systems. Second, the VMC calculations themselves can give a good upper bound for the exact ground-state energy of many-electron systems when an appropriate form for the variational wave function is employed. In addition to a conventional form for the variational wave function that is expressed by the product of the Jastrow factor and the single Slater determinant, recent investigations have found that the expression in terms of many Slater determinants can provide a very accurate description of electronic states, especially for intrinsically multi-reference systems, which corresponds to the configuration interaction (CI) expansion method in conventional quantum chemistry. Although the CI method needs an extremely large number of expansion coefficients as variational parameters, the number of parameters can be reduced efficiently by using an idea of (molecular) orbital pair correlations. The purpose of the present paper is to examine the feasibility and usefulness of this idea within the framework of VMC method with the Jastrow–Slater wave function, where the ground state of a water molecule is investigated as a test case. In this study we rely on a variational wave function composed of a product of the Jastrow and Slater parts as 1⁄4 J S; ð1Þ
We report an accurate study of interactions between benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory using different van der Waals functionals. In our quantum Monte Carlo (QMC) calculations, we use accurate correlated trial wave functions including three-body Jastrow factors and backflow transformations. We consider two benzene molecules in the parallel displaced geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the coupled-cluster theory through perturbative triplets/complete basis set limit is -2.65(2) kcal/mol [Miliordos et al., J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods.
We review the use of continuum quantum Monte Carlo (QMC) methods for the calculation of energy gaps from first principles, and present a broad set of excited-state calculations carried out with the variational and fixed-node diffusion QMC methods on atoms, molecules, and solids. We propose a finite-size-error correction scheme for bulk energy gaps calculated in finite cells subject to periodic boundary conditions. We show that finite-size effects are qualitatively different in two-dimensional materials, demonstrating the effect in a QMC calculation of the band gap and exciton binding energy of monolayer phosphorene. We investigate the fixed-node errors in diffusion Monte Carlo gaps evaluated with Slater-Jastrow trial wave functions by examining the effects of backflow transformations, and also by considering the formation of restricted multideterminant expansions for excited-state wave functions. For several molecules, we examine the importance of structural relaxation in the excited state in determining excited-state energies. We study the feasibility of using variational Monte Carlo with backflow correlations to obtain accurate excited-state energies at reduced computational cost, finding that this approach can be valid. We find that diffusion Monte Carlo gap calculations can be performed with much larger time steps than are typically required to converge the total energy, at significantly diminished computational expense, but that in order to alleviate fixed-node errors in calculations on solids the inclusion of backflow correlations is sometimes necessary.
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).