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:1309.0427 · Journal of Chemical Theory and Computation

Molecular Properties by Quantum Monte Carlo: An Investigation on the Role of the Wave Function Ansatz and the Basis Set in the Water Molecule

Andrea Zen, Ye Luo, Sandro Sorella, Leonardo Guidoni

Published Aug 28, 2013• 44 citations
arXivHTMLPDF

Abstract

Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely, the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudopotential, and the basis set for QMC calculations. We also introduce a new method for the computation of forces with finite variance on open systems and a new strategy for the definition of the atomic orbitals involved in the Jastrow-Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.

DOI
10.1021/ct400382m

Related papers

Molecular Properties by Quantum Monte Carlo: An Investigation on the Role of the Wave Function Ansatz and the Basis Set in the Water Molecule

Andrea Zen, Luo, Y, Sorella, S +1 more

Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely: the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudo potential, and the basis set for QMC calculations. We also introduce a new strategy for the definition of the atomic orbitals involved in the Jastrow - Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.

•Oct 1, 2013•37 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

High-accuracy electronic structure quantum Monte Carlo for molecular systems

Tommaso Torelli

Quantum Monte Carlo (QMC) is one of the most promising methods for solving quantum many-body problems. QMC combines known analytical properties of wave functions, results from other methods, and stochastic techniques, into a powerful tool for investigation of the electronic structure of real materials. Furthermore, the method has a favorable scaling ( O[N3] in the number of electrons), perfect scalability on parallel architectures, and can be applied to a wide range of systems. Our work has focused both on expanding the capabilities of QMC for calculation of quantities beyond energies, such as interatomic forces, and on application of the method to challenging problems in nanostructure materials research. For the calculation of interatomic forces we have implemented a finite difference correlated sampling method. The correlated sampling enables us to avoid the problem of statistical noise in evaluating the energy differences and allows us to obtain forces with an average 1% error. We have analyzed the effect of the method's main approximations on the estimate of the forces, both by theoretical arguments and by performing a variety of tests on small systems. Furthermore, we have used this newly developed method for estimating equilibrium geometries in a problem where other methods have failed. QMC enabled us to treat the electron correlation effects with high accuracy and to provide accurate predictions in a number of interesting applications. In particular, for carbon rings of C4N+2 stoichiometry we have resolved a long-standing problem regarding the competition between the aromaticity and second order Jahn-Teller effects. We have found that the dominant mechanism at small and intermediate sizes, which stabilizes the bond angle and bond length alternated geometries, is the secondorder Jahn-Teller effect while the aromatic isomer is always found to be a transition state. These high accuracy calculations which have involved multi-reference wave functions, extensive variational optimizations, and the use of massively parallel platforms demonstrate clearly the unique capabilities of the QMC approach.

•Jun 1, 2001•0 citations

Ab initio molecular dynamics of water by quantum Monte Carlo

Ye Luo, Sandro Sorella

The chapter 2, we deal with the challenge b). It focuses on the variational Monte Carlo (VMC) and the wavefunction optimization methods based on VMC. The performance of different methods are displayed through the op- timization of the Jastrow factor in our test case Beryllium dimer and the efficiency is improving surprisingly during the evolution of these methods. In chapter 3, we focus on the challenge a). It describes the wavefunc- tion ansatz used by our simulation. In this thesis, we introduce the atomic hybrid orbitals which significantly increase the compactness of our wavefunc- tion without hurting accuracy. This chapter also explain how to optimize the determinant in a way that the number of variational parameters scales only linearly with the system size. This further helps the efficiency of the wavefunction optimization. In chapters 4 and 5, the issue c) is explained in detail. In chapter 4, a second order Langevin dynamics (SLD) scheme is devised particularly for QMC and this thesis improves this scheme by developing a better integration method. Here, we also highlight the remarkable power of the force covari- ance matrix which can be defined only in QMC and is capable of accelerating the slow modes of a dynamics. In chapter 5, this SLD for QMC is validated through intensive benchmarking on the calculation of the vibrational frequen- cies of water and other small molecules. It is shown that many systematic biases in our MD scheme and QMC evaluation can be controlled so that we are confident to push forward this ab initio molecular dynamics for applica- tions on large systems. Finally in chapter 6, we perform the simulation of liquid water with all the preparation done in the previous chapters. The results are encouraging since we’ve closed the discrepancy of the peak positions of RDFs between experi- ments and ab initio simulations. The power of QMC is also demonstrated by the fact that the shapes of our RDFs are much less structured than previous DFT-based ab initio simulations even if the two water molecule interaction is dealt with the same level of accuracy as the DFT/BLYP calculation. In this chapter, we have also studied the features of hydrogen bonds in our simulation of liquid water. All our results indicate that it is important to consider the quantum nature of the ions for a faithful description of liquid water. This will be left for future studies, possible in principle even within the QMC approach.

•Oct 31, 2014•0 citations

Computing the energy of a water molecule using multideterminants: A simple, efficient algorithm

Bryan K. Clark, Miguel A. Morales, Jeremy McMinis +2 more

Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave functions are critical to ascertaining new physics. One such wave function is the multi-Slater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wave functions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally, we implement this method and use it to compute the ground state energy of a water molecule.

arXiv:1106.2456•The Journal of Chemical Physics•Dec 27, 2011•71 citations

Computing the energy of a water molecule using MultiDeterminants: A simple, efficient algorithm

Clark, B K, Morales, M A, McMinis, J +2 more

Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave functions are critical to ascertaining new physics. One such wave function is the multi-Slater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wave functions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally, we implement this method and use it to compute the ground state energy of a water molecule.

•May 3, 2011•60 citations
Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely, the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudopotential, and the basis set for QMC calculations. We also introduce a new method for the computation of forces with finite variance on open systems and a new strategy for the definition of the atomic orbitals involved in the Jastrow-Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.
Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely: the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudo potential, and the basis set for QMC calculations. We also introduce a new strategy for the definition of the atomic orbitals involved in the Jastrow - Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.
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Þ
Quantum Monte Carlo (QMC) is one of the most promising methods for solving quantum many-body problems. QMC combines known analytical properties of wave functions, results from other methods, and stochastic techniques, into a powerful tool for investigation of the electronic structure of real materials. Furthermore, the method has a favorable scaling ( O[N3] in the number of electrons), perfect scalability on parallel architectures, and can be applied to a wide range of systems. Our work has focused both on expanding the capabilities of QMC for calculation of quantities beyond energies, such as interatomic forces, and on application of the method to challenging problems in nanostructure materials research. For the calculation of interatomic forces we have implemented a finite difference correlated sampling method. The correlated sampling enables us to avoid the problem of statistical noise in evaluating the energy differences and allows us to obtain forces with an average 1% error. We have analyzed the effect of the method's main approximations on the estimate of the forces, both by theoretical arguments and by performing a variety of tests on small systems. Furthermore, we have used this newly developed method for estimating equilibrium geometries in a problem where other methods have failed. QMC enabled us to treat the electron correlation effects with high accuracy and to provide accurate predictions in a number of interesting applications. In particular, for carbon rings of C4N+2 stoichiometry we have resolved a long-standing problem regarding the competition between the aromaticity and second order Jahn-Teller effects. We have found that the dominant mechanism at small and intermediate sizes, which stabilizes the bond angle and bond length alternated geometries, is the secondorder Jahn-Teller effect while the aromatic isomer is always found to be a transition state. These high accuracy calculations which have involved multi-reference wave functions, extensive variational optimizations, and the use of massively parallel platforms demonstrate clearly the unique capabilities of the QMC approach.
The chapter 2, we deal with the challenge b). It focuses on the variational Monte Carlo (VMC) and the wavefunction optimization methods based on VMC. The performance of different methods are displayed through the op- timization of the Jastrow factor in our test case Beryllium dimer and the efficiency is improving surprisingly during the evolution of these methods. In chapter 3, we focus on the challenge a). It describes the wavefunc- tion ansatz used by our simulation. In this thesis, we introduce the atomic hybrid orbitals which significantly increase the compactness of our wavefunc- tion without hurting accuracy. This chapter also explain how to optimize the determinant in a way that the number of variational parameters scales only linearly with the system size. This further helps the efficiency of the wavefunction optimization. In chapters 4 and 5, the issue c) is explained in detail. In chapter 4, a second order Langevin dynamics (SLD) scheme is devised particularly for QMC and this thesis improves this scheme by developing a better integration method. Here, we also highlight the remarkable power of the force covari- ance matrix which can be defined only in QMC and is capable of accelerating the slow modes of a dynamics. In chapter 5, this SLD for QMC is validated through intensive benchmarking on the calculation of the vibrational frequen- cies of water and other small molecules. It is shown that many systematic biases in our MD scheme and QMC evaluation can be controlled so that we are confident to push forward this ab initio molecular dynamics for applica- tions on large systems. Finally in chapter 6, we perform the simulation of liquid water with all the preparation done in the previous chapters. The results are encouraging since we’ve closed the discrepancy of the peak positions of RDFs between experi- ments and ab initio simulations. The power of QMC is also demonstrated by the fact that the shapes of our RDFs are much less structured than previous DFT-based ab initio simulations even if the two water molecule interaction is dealt with the same level of accuracy as the DFT/BLYP calculation. In this chapter, we have also studied the features of hydrogen bonds in our simulation of liquid water. All our results indicate that it is important to consider the quantum nature of the ions for a faithful description of liquid water. This will be left for future studies, possible in principle even within the QMC approach.
Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave functions are critical to ascertaining new physics. One such wave function is the multi-Slater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wave functions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally, we implement this method and use it to compute the ground state energy of a water molecule.
Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave functions are critical to ascertaining new physics. One such wave function is the multi-Slater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wave functions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally, we implement this method and use it to compute the ground state energy of a water molecule.