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

Wave functions for quantum Monte Carlo calculations in solids: Orbitals from density functional theory with hybrid exchange-correlation functionals

Jindřich Kolorenč, Shuming Hu, Luboš Mitáš

Published Sep 8, 2010• 64 citations
arXivHTMLPDF

Abstract

We investigate how the fixed-node diffusion Monte Carlo energy of solids depends on single-particle orbitals used in Slater-Jastrow wave functions. We demonstrate that the dependence can be significant, in particular in the case of $3d$ transition-metal compounds, which we adopt as examples. We illustrate how exchange-correlation functionals with variable exact-exchange component can be exploited to reduce the fixed-node errors. On the basis of these results we argue that the fixed-node quantum Monte Carlo provides a variational approach for optimization of effective single-particle Hamiltonians with parameters.

DOI
10.1103/physrevb.82.115108

Related papers

Orbital-dependent backflow wave functions for real-space quantum Monte Carlo

Markus Holzmann, Saverio Moroni

We present and motivate an efficient way to include orbital-dependent many-body correlations in trial wave functions of real-space quantum Monte Carlo methods for use in electronic structure calculations. We apply our new orbital-dependent backflow wave function to calculate ground state energies of the first row atoms using variational and diffusion Monte Carlo methods. The systematic overall gain of correlation energy with respect to single determinant Jastrow-Slater wave functions is competitive with the best single determinant trial wave functions currently available. The computational cost per Monte Carlo step is comparable to that of simple backflow calculations.

arXiv:1910.07167•Physical review. B./Physical review. B•Feb 12, 2019•18 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

Quantum Monte Carlo studies of density functional theory

Maziar Nekovee, W. M. C. Foulkes, Richard J. Needs

Abstract It is difficult to imagine that density functional theory (DFT) could be useful in improving quantum Monte Carlo methods, but it is easy to imagine that QMC could be useful in improving density functional methods. Since QMC can produce accurate wavefunctions and since wavefunctions lead to all other properties one should be able, in principle, to use QMC to predict the exchangecorrelation energy functional Exc[n(r)], the key (unknown) quantity required in DFT. In practice, that is not so easy, but this paper shows how QMC can take steps in that direction by providing insights gained from comparisons of results from QMC and DFT. The authors considered the case of the strongly inhomogeneous electron gas using variational Q:tvIC to obtain wavefunctions 1/J>. and energies E[7/J>.] for the ground state of the Hamiltonian iI>- = T +AVee+V >. associated with the coupling constant A in the range 0 to 1. Optimization of a Slater-Jastrow many-body wavefunction was made by minimizing a penalty function given by a weighted sum of the variance in local energies and a sum of squares of electron density differences. The variational parameters in 7/J>. and v>. were repeatedly varied to reach convergence followed by additional calculations to obtain the local values of the exchange-correlation energies e c and the exchange-correlation holes n c from which Exe and nxc could be determined.

•Jun 18, 2007•0 citations

Starting-point-independent quantum Monte Carlo calculations of iron oxide

Joshua Townsend, Sergio D. Pineda Flores, Raymond C. Clay +5 more

Quantum Monte Carlo (QMC) methods are useful for studies of strongly correlated materials because they are many body in nature and use the physical Hamiltonian. Typical calculations assume as a starting point a wave function constructed from single-particle orbitals obtained from one-body methods, e.g., density functional theory. However, mean-field-derived wave functions can sometimes lead to systematic QMC biases if the mean-field result poorly describes the true ground state. Here, we study the accuracy and flexibility of QMC trial wave functions using variational and fixed-node diffusion QMC estimates of the total spin density and lattice distortion of antiferromagnetic iron oxide (FeO) in the ground state $B1$ crystal structure. We found that for relatively simple wave functions the predicted lattice distortion was controlled by the choice of single-particle orbitals used to construct the wave function, rather than by subsequent wave function optimization techniques within QMC. By optimizing the orbitals with QMC, we then demonstrate starting-point independence of the trial wave function with respect to the method by which the orbitals were constructed by demonstrating convergence of the energy, spin density, and predicted lattice distortion for two qualitatively different sets of orbitals. The results suggest that orbital optimization is a promising method for accurate many-body calculations of strongly correlated condensed phases.

Physical review. B./Physical review. B•Oct 30, 2020•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

Performance of the Diffusion Quantum Monte Carlo Method with a Single-Slater-Jastrow Trial Wavefunction Using Natural Orbitals and Density Functional Theory Orbitals on Atomization Energies of the Gaussian-2 Set

Ting Wang, Xiaojun Zhou, Fan Wang

Performance of the fixed-node diffusion quantum Monte Carlo method (FN-DMC) with a single Slater-Jastrow trial wavefunction using natural orbitals (NOs) from MP2, CCSD, CCSD(T), and CASSCF as well as density functional theory orbitals on atomization energies of the molecules in the Gaussian-2 set is investigated in this work. The effects of spin contamination and pseudopotentials (PPs) are also studied. Our results show that the DMC energy with NOs from MP2 or CCSD(T) is the lowest on average, whereas that from CASSCF is only lower than the DMC energy using Hartree-Fock orbitals. Atomization energies are generally underestimated with DMC, and mean absolute deviations of DMC atomization energies are about 2.8 kcal/mol with NOs from MP2 or CCSD(T) and about 2.7 kcal/mol with NOs from CASSCF or B3LYP orbitals. The better performance of the latter two orbitals is due to a more effective error cancellation. The accuracy of the present FN-DMC is similar to CCSD(T)/aug-cc-pVQZ on atomization energies. In addition, the error of DMC atomization energies tends to be larger for molecules with multiple bonds. DMC energies of the open-shell systems with spin-restricted orbitals are generally higher than those with spin-unrestricted orbitals if spin contamination is not serious and the atomization energies are improved to some extent. Furthermore, our results indicate that the error of the PPs is rather small at the CCSD(T) level, but could be more pronounced in DMC calculations and DMC atomization energies are improved with some newly developed PPs.

The Journal of Physical Chemistry A•Apr 5, 2019•18 citations
We investigate how the fixed-node diffusion Monte Carlo energy of solids depends on single-particle orbitals used in Slater-Jastrow wave functions. We demonstrate that the dependence can be significant, in particular in the case of 3d3d3d transition-metal compounds, which we adopt as examples. We illustrate how exchange-correlation functionals with variable exact-exchange component can be exploited to reduce the fixed-node errors. On the basis of these results we argue that the fixed-node quantum Monte Carlo provides a variational approach for optimization of effective single-particle Hamiltonians with parameters.
We present and motivate an efficient way to include orbital-dependent many-body correlations in trial wave functions of real-space quantum Monte Carlo methods for use in electronic structure calculations. We apply our new orbital-dependent backflow wave function to calculate ground state energies of the first row atoms using variational and diffusion Monte Carlo methods. The systematic overall gain of correlation energy with respect to single determinant Jastrow-Slater wave functions is competitive with the best single determinant trial wave functions currently available. The computational cost per Monte Carlo step is comparable to that of simple backflow calculations.
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.
Abstract It is difficult to imagine that density functional theory (DFT) could be useful in improving quantum Monte Carlo methods, but it is easy to imagine that QMC could be useful in improving density functional methods. Since QMC can produce accurate wavefunctions and since wavefunctions lead to all other properties one should be able, in principle, to use QMC to predict the exchangecorrelation energy functional Exc[n(r)], the key (unknown) quantity required in DFT. In practice, that is not so easy, but this paper shows how QMC can take steps in that direction by providing insights gained from comparisons of results from QMC and DFT. The authors considered the case of the strongly inhomogeneous electron gas using variational Q:tvIC to obtain wavefunctions 1/J>. and energies E[7/J>.] for the ground state of the Hamiltonian iI>- = T +AVee+V >. associated with the coupling constant A in the range 0 to 1. Optimization of a Slater-Jastrow many-body wavefunction was made by minimizing a penalty function given by a weighted sum of the variance in local energies and a sum of squares of electron density differences. The variational parameters in 7/J>. and v>. were repeatedly varied to reach convergence followed by additional calculations to obtain the local values of the exchange-correlation energies e c and the exchange-correlation holes n c from which Exe and nxc could be determined.
Quantum Monte Carlo (QMC) methods are useful for studies of strongly correlated materials because they are many body in nature and use the physical Hamiltonian. Typical calculations assume as a starting point a wave function constructed from single-particle orbitals obtained from one-body methods, e.g., density functional theory. However, mean-field-derived wave functions can sometimes lead to systematic QMC biases if the mean-field result poorly describes the true ground state. Here, we study the accuracy and flexibility of QMC trial wave functions using variational and fixed-node diffusion QMC estimates of the total spin density and lattice distortion of antiferromagnetic iron oxide (FeO) in the ground state B1B1B1 crystal structure. We found that for relatively simple wave functions the predicted lattice distortion was controlled by the choice of single-particle orbitals used to construct the wave function, rather than by subsequent wave function optimization techniques within QMC. By optimizing the orbitals with QMC, we then demonstrate starting-point independence of the trial wave function with respect to the method by which the orbitals were constructed by demonstrating convergence of the energy, spin density, and predicted lattice distortion for two qualitatively different sets of orbitals. The results suggest that orbital optimization is a promising method for accurate many-body calculations of strongly correlated condensed phases.
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Þ
Performance of the fixed-node diffusion quantum Monte Carlo method (FN-DMC) with a single Slater-Jastrow trial wavefunction using natural orbitals (NOs) from MP2, CCSD, CCSD(T), and CASSCF as well as density functional theory orbitals on atomization energies of the molecules in the Gaussian-2 set is investigated in this work. The effects of spin contamination and pseudopotentials (PPs) are also studied. Our results show that the DMC energy with NOs from MP2 or CCSD(T) is the lowest on average, whereas that from CASSCF is only lower than the DMC energy using Hartree-Fock orbitals. Atomization energies are generally underestimated with DMC, and mean absolute deviations of DMC atomization energies are about 2.8 kcal/mol with NOs from MP2 or CCSD(T) and about 2.7 kcal/mol with NOs from CASSCF or B3LYP orbitals. The better performance of the latter two orbitals is due to a more effective error cancellation. The accuracy of the present FN-DMC is similar to CCSD(T)/aug-cc-pVQZ on atomization energies. In addition, the error of DMC atomization energies tends to be larger for molecules with multiple bonds. DMC energies of the open-shell systems with spin-restricted orbitals are generally higher than those with spin-unrestricted orbitals if spin contamination is not serious and the atomization energies are improved to some extent. Furthermore, our results indicate that the error of the PPs is rather small at the CCSD(T) level, but could be more pronounced in DMC calculations and DMC atomization energies are improved with some newly developed PPs.