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

General Correlated Geminal Ansatz for Electronic Structure Calculations: Exploiting Pfaffians in Place of Determinants

Claudio Genovese, Tomonori Shirakawa, Kousuke Nakano, Sandro Sorella

Published Aug 17, 2020• 21 citations
arXivHTMLPDF

Abstract

High Resolution Image Download MS PowerPoint Slide We propose here a single Pfaffian correlated variational ansatz that dramatically improves the accuracy with respect to the single determinant one, while remaining at a similar computational cost. A much larger correlation energy is indeed determined by the most general two electron pairing function, including both singlet and triplet channels, combined with a many-body Jastrow factor, including all possible spin–spin, spin–density, and density–density terms. The main technical ingredient to exploit this accuracy is the use of the Pfaffian for antisymmetrizing a highly correlated pairing function, thus recovering the Fermi statistics for electrons with an affordable computational cost. Moreover, the application of the diffusion Monte Carlo, within the fixed node approximation, allows us to obtain very accurate binding energies for the first preliminary calculations reported in this study: C 2, N 2, and O 2 and the benzene molecule. This is promising and remarkable, considering that they represent extremely difficult molecules even for computationally demanding multideterminant approaches, and opens therefore the way for realistic and accurate electronic simulations with an algorithm scaling at most as the fourth power of the number of electrons.

DOI
10.1021/acs.jctc.0c00165

Related papers

Toward Chemical Accuracy Using the Jastrow Correlated Antisymmetrized Geminal Power Ansatz

Abhishek Raghav, Ryo Maezono, Kenta Hongo +2 more

High Resolution Image Download MS PowerPoint Slide Herein, we report accurate atomization energy calculations for 55 molecules in the Gaussian-2 (G2) set using lattice regularized diffusion Monte Carlo (LRDMC). We compare the Jastrow–Slater determinant ansatz with a more flexible JsAGPs (Jastrow correlated antisymmetrized geminal power with singlet correlation) ansatz . AGPs is built from pairing functions, which explicitly include pairwise correlations among electrons, and hence, this ansatz is expected to be more efficient in recovering the correlation energy. The AGPs wave functions are first optimized at the variational Monte Carlo (VMC) level, which includes both the Jastrow factor and the nodal surface optimization. This is followed by the LRDMC projection of the ansatz . Remarkably, for many molecules, the LRDMC atomization energies obtained using the JsAGPs ansatz reach chemical accuracy (∼1 kcal/mol), and for most other molecules, the atomization energies are accurate within ∼5 kcal/mol. We obtained a mean absolute deviation of 1.6 kcal/mol with JsAGPs and 3.2 kcal/mol with JDFT ( J astrow factor + Slater determinant with DFT orbitals) ansatzes . This work shows the effectiveness of the flexible AGPs ansatz for atomization energy calculations and electronic structure simulations in general.

arXiv:2211.05447•Journal of Chemical Theory and Computation•Apr 4, 2023•12 citations

Pfaffian Pairing Wave Functions in Electronic-Structure Quantum Monte Carlo Simulations

Michal Bajdich, Luboš Mitáš, G. Drobný +2 more

We investigate the accuracy of trial wave functions for quantum Monte Carlo based on Pfaffian functional form with singlet and triplet pairing. Using a set of first row atoms and molecules we find that these wave functions provide very consistent and systematic behavior in recovering the correlation energies on the level of 95%. In order to get beyond this limit we explore the possibilities of multi-Pfaffian pairing wave functions. We show that a small number of Pfaffians recovers another large fraction of the missing correlation energy comparable to the larger-scale configuration interaction wave functions. We also find that Pfaffians lead to substantial improvements in fermion nodes when compared to Hartree-Fock wave functions.

arXiv:cond-mat/0512327•Physical Review Letters•Apr 5, 2006•145 citations

Efficient calculation of unbiased atomic forces in ab initio variational Monte Carlo

Kousuke Nakano, Michele Casula, Giacomo Tenti

quantum Monte Carlo (QMC) is a state-of-the-art numerical approach for evaluating accurate expectation values of many-body wave functions. However, one of the major drawbacks that still hinders widespread QMC applications is the lack of an affordable scheme to compute unbiased atomic forces. In this study, we propose an efficient method to obtain unbiased atomic forces and pressures in the variational Monte Carlo (VMC) framework with the Jastrow-correlated Slater determinant ansatz or the Jastrow antisymmetrized geminal power ansatz, exploiting the gauge-invariant and locality properties of their geminal representation. We demonstrate the effectiveness of our method for H 2 and Cl 2 molecules and for the cubic boron nitride crystal. Our framework has a better algorithmic scaling with the system size than the traditional finite-difference method and, in practical applications, is as efficient as single-point VMC calculations. Thus, it paves the way to study dynamical properties of materials, such as phonons, and is beneficial for pursuing more reliable machine-learning interatomic potentials based on unbiased VMC forces. Published by the American Physical Society 2024

arXiv:2312.17608•Physical review. B./Physical review. B•May 28, 2024•9 citations

Basis-free neural-network geminal and Jastrow factors for variational Monte Carlo

Jan Kessler, Thomas D. Kühne

Neural-network quantum states offer a flexible route to compact many-electron wave functions, but their practical accuracy depends strongly on how fermionic antisymmetry, electron correlation, and optimization noise are treated. Here we combine an antisymmetrized geminal power (AGP) determinant with feed-forward neural networks that replace conventional basis-set expansions in the geminal and in two Jastrow-factor constructions. The resulting basis-free Jastrow--AGP ansatz is optimized by variational Monte Carlo and is designed to separate two tasks: the AGP part defines the nodal surface, while the neural-network Jastrow factor recovers dynamical correlation at fixed nodes. This separation makes it possible to distinguish errors associated with dynamical correlation from those caused by static, multireference correlation. Applications to the hydrogen molecule and the rectangular hydrogen tetramer show sub-millihartree accuracy when the AGP nodes are adequate, and expose the residual nodal limitation near the large-radius square geometry of the hydrogen tetramer. These results clarify where neural-network building blocks can improve a compact geminal ansatz and where additional nodal flexibility is required.

arXiv:2605.17083•May 16, 2026•0 citations

Influence of single particle orbital sets and configuration selection on multideterminant wavefunctions in quantum Monte Carlo

Raymond C. Clay, Miguel A. Morales

Multideterminant wavefunctions, while having a long history in quantum chemistry, are increasingly being used in highly accurate quantum Monte Carlo calculations. Since the accuracy of QMC is ultimately limited by the quality of the trial wavefunction, multi-Slater determinants wavefunctions offer an attractive alternative to Slater-Jastrow and more sophisticated wavefunction ansatz for several reasons. They can be efficiently calculated, straightforwardly optimized, and systematically improved by increasing the number of included determinants. In spite of their potential, however, the convergence properties of multi-Slater determinant wavefunctions with respect to orbital set choice and excited determinant selection are poorly understood, which hinders the application of these wavefunctions to large systems and solids. In this paper, by performing QMC calculations on the equilibrium and stretched carbon dimer, we find that convergence of the recovered correlation energy with respect to number of determinants can depend quite strongly on basis set and determinant selection methods, especially where there is strong correlation. We demonstrate that properly chosen orbital sets and determinant selection techniques from quantum chemistry methods can dramatically reduce the required number of determinants (and thus the computational cost) to reach a given accuracy, which we argue shows clear need for an automatic QMC-only method for selecting determinants and generating optimal orbital sets.

The Journal of Chemical Physics•Jun 16, 2015•23 citations

Beyond Single-Reference Fixed-Node Approximation in Ab Initio Diffusion Monte Carlo Using Antisymmetrized Geminal Power Applied to Systems with Hundreds of Electrons

Kousuke Nakano, Sandro Sorella, Dario Alfè +1 more

High Resolution Image Download MS PowerPoint Slide Diffusion Monte Carlo (DMC) is an exact technique to project out the ground state (GS) of a Hamiltonian. Since the GS is always bosonic, in Fermionic systems, the projection needs to be carried out while imposing antisymmetric constraints, which is a nondeterministic polynomial hard problem. In practice, therefore, the application of DMC on electronic structure problems is made by employing the fixed-node (FN) approximation, consisting of performing DMC with the constraint of having a fixed, predefined nodal surface. How do we get the nodal surface? The typical approach, applied in systems having up to hundreds or even thousands of electrons, is to obtain the nodal surface from a preliminary mean-field approach (typically, a density functional theory calculation) used to obtain a single Slater determinant. This is known as single reference. In this paper, we propose a new approach, applicable to systems as large as the C 60 fullerene, which improves the nodes by going beyond the single reference. In practice, we employ an implicitly multireference ansatz (antisymmetrized geminal power wave function constraint with molecular orbitals), initialized on the preliminary mean-field approach, which is relaxed by optimizing a few parameters of the wave function determining the nodal surface by minimizing the FN-DMC energy. We highlight the improvements of the proposed approach over the standard single-reference method on several examples and, where feasible, the computational gain over the standard multireference ansatz, which makes the methods applicable to large systems. We also show that physical properties relying on relative energies, such as binding energies, are affordable and reliable within the proposed scheme.

arXiv:2402.01458•Journal of Chemical Theory and Computation•May 24, 2024•6 citations
High Resolution Image Download MS PowerPoint Slide We propose here a single Pfaffian correlated variational ansatz that dramatically improves the accuracy with respect to the single determinant one, while remaining at a similar computational cost. A much larger correlation energy is indeed determined by the most general two electron pairing function, including both singlet and triplet channels, combined with a many-body Jastrow factor, including all possible spin–spin, spin–density, and density–density terms. The main technical ingredient to exploit this accuracy is the use of the Pfaffian for antisymmetrizing a highly correlated pairing function, thus recovering the Fermi statistics for electrons with an affordable computational cost. Moreover, the application of the diffusion Monte Carlo, within the fixed node approximation, allows us to obtain very accurate binding energies for the first preliminary calculations reported in this study: C 2, N 2, and O 2 and the benzene molecule. This is promising and remarkable, considering that they represent extremely difficult molecules even for computationally demanding multideterminant approaches, and opens therefore the way for realistic and accurate electronic simulations with an algorithm scaling at most as the fourth power of the number of electrons.
High Resolution Image Download MS PowerPoint Slide Herein, we report accurate atomization energy calculations for 55 molecules in the Gaussian-2 (G2) set using lattice regularized diffusion Monte Carlo (LRDMC). We compare the Jastrow–Slater determinant ansatz with a more flexible JsAGPs (Jastrow correlated antisymmetrized geminal power with singlet correlation) ansatz . AGPs is built from pairing functions, which explicitly include pairwise correlations among electrons, and hence, this ansatz is expected to be more efficient in recovering the correlation energy. The AGPs wave functions are first optimized at the variational Monte Carlo (VMC) level, which includes both the Jastrow factor and the nodal surface optimization. This is followed by the LRDMC projection of the ansatz . Remarkably, for many molecules, the LRDMC atomization energies obtained using the JsAGPs ansatz reach chemical accuracy (∼1 kcal/mol), and for most other molecules, the atomization energies are accurate within ∼5 kcal/mol. We obtained a mean absolute deviation of 1.6 kcal/mol with JsAGPs and 3.2 kcal/mol with JDFT ( J astrow factor + Slater determinant with DFT orbitals) ansatzes . This work shows the effectiveness of the flexible AGPs ansatz for atomization energy calculations and electronic structure simulations in general.
We investigate the accuracy of trial wave functions for quantum Monte Carlo based on Pfaffian functional form with singlet and triplet pairing. Using a set of first row atoms and molecules we find that these wave functions provide very consistent and systematic behavior in recovering the correlation energies on the level of 95%. In order to get beyond this limit we explore the possibilities of multi-Pfaffian pairing wave functions. We show that a small number of Pfaffians recovers another large fraction of the missing correlation energy comparable to the larger-scale configuration interaction wave functions. We also find that Pfaffians lead to substantial improvements in fermion nodes when compared to Hartree-Fock wave functions.
quantum Monte Carlo (QMC) is a state-of-the-art numerical approach for evaluating accurate expectation values of many-body wave functions. However, one of the major drawbacks that still hinders widespread QMC applications is the lack of an affordable scheme to compute unbiased atomic forces. In this study, we propose an efficient method to obtain unbiased atomic forces and pressures in the variational Monte Carlo (VMC) framework with the Jastrow-correlated Slater determinant ansatz or the Jastrow antisymmetrized geminal power ansatz, exploiting the gauge-invariant and locality properties of their geminal representation. We demonstrate the effectiveness of our method for H 2 and Cl 2 molecules and for the cubic boron nitride crystal. Our framework has a better algorithmic scaling with the system size than the traditional finite-difference method and, in practical applications, is as efficient as single-point VMC calculations. Thus, it paves the way to study dynamical properties of materials, such as phonons, and is beneficial for pursuing more reliable machine-learning interatomic potentials based on unbiased VMC forces. Published by the American Physical Society 2024
Neural-network quantum states offer a flexible route to compact many-electron wave functions, but their practical accuracy depends strongly on how fermionic antisymmetry, electron correlation, and optimization noise are treated. Here we combine an antisymmetrized geminal power (AGP) determinant with feed-forward neural networks that replace conventional basis-set expansions in the geminal and in two Jastrow-factor constructions. The resulting basis-free Jastrow--AGP ansatz is optimized by variational Monte Carlo and is designed to separate two tasks: the AGP part defines the nodal surface, while the neural-network Jastrow factor recovers dynamical correlation at fixed nodes. This separation makes it possible to distinguish errors associated with dynamical correlation from those caused by static, multireference correlation. Applications to the hydrogen molecule and the rectangular hydrogen tetramer show sub-millihartree accuracy when the AGP nodes are adequate, and expose the residual nodal limitation near the large-radius square geometry of the hydrogen tetramer. These results clarify where neural-network building blocks can improve a compact geminal ansatz and where additional nodal flexibility is required.
Multideterminant wavefunctions, while having a long history in quantum chemistry, are increasingly being used in highly accurate quantum Monte Carlo calculations. Since the accuracy of QMC is ultimately limited by the quality of the trial wavefunction, multi-Slater determinants wavefunctions offer an attractive alternative to Slater-Jastrow and more sophisticated wavefunction ansatz for several reasons. They can be efficiently calculated, straightforwardly optimized, and systematically improved by increasing the number of included determinants. In spite of their potential, however, the convergence properties of multi-Slater determinant wavefunctions with respect to orbital set choice and excited determinant selection are poorly understood, which hinders the application of these wavefunctions to large systems and solids. In this paper, by performing QMC calculations on the equilibrium and stretched carbon dimer, we find that convergence of the recovered correlation energy with respect to number of determinants can depend quite strongly on basis set and determinant selection methods, especially where there is strong correlation. We demonstrate that properly chosen orbital sets and determinant selection techniques from quantum chemistry methods can dramatically reduce the required number of determinants (and thus the computational cost) to reach a given accuracy, which we argue shows clear need for an automatic QMC-only method for selecting determinants and generating optimal orbital sets.
High Resolution Image Download MS PowerPoint Slide Diffusion Monte Carlo (DMC) is an exact technique to project out the ground state (GS) of a Hamiltonian. Since the GS is always bosonic, in Fermionic systems, the projection needs to be carried out while imposing antisymmetric constraints, which is a nondeterministic polynomial hard problem. In practice, therefore, the application of DMC on electronic structure problems is made by employing the fixed-node (FN) approximation, consisting of performing DMC with the constraint of having a fixed, predefined nodal surface. How do we get the nodal surface? The typical approach, applied in systems having up to hundreds or even thousands of electrons, is to obtain the nodal surface from a preliminary mean-field approach (typically, a density functional theory calculation) used to obtain a single Slater determinant. This is known as single reference. In this paper, we propose a new approach, applicable to systems as large as the C 60 fullerene, which improves the nodes by going beyond the single reference. In practice, we employ an implicitly multireference ansatz (antisymmetrized geminal power wave function constraint with molecular orbitals), initialized on the preliminary mean-field approach, which is relaxed by optimizing a few parameters of the wave function determining the nodal surface by minimizing the FN-DMC energy. We highlight the improvements of the proposed approach over the standard single-reference method on several examples and, where feasible, the computational gain over the standard multireference ansatz, which makes the methods applicable to large systems. We also show that physical properties relying on relative energies, such as binding energies, are affordable and reliable within the proposed scheme.