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

Quantum Monte Carlo with coupled-cluster wave functions

Alessandro Roggero, Abhishek Mukherjee, Francesco Pederiva

Published Sep 23, 2013• 33 citations
arXivHTMLPDF

Abstract

We introduce a many-body method that combines two powerful many-body techniques, viz., quantum Monte Carlo and coupled cluster theory. Coupled cluster wave functions are introduced as importance functions in a Monte Carlo method designed for the configuration interaction framework to provide rigorous upper bounds to the ground-state energy. We benchmark our method on the homogeneous electron gas in momentum space. The importance function used is the coupled cluster doubles wave function. We show that the computational resources required in our method scale polynomially with system size. Our energy upper bounds are in very good agreement with previous calculations of similar accuracy, and they can be systematically improved by including higher order excitations in the coupled cluster wave function.

DOI
10.1103/physrevb.88.115138

Related papers

Coupled Cluster Theory in Materials Science

Igor Ying Zhang, Andreas Grüneis

The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.

arXiv:2004.06424•Frontiers in Materials•Jun 5, 2019•138 citations

Auxiliary-Field Quantum Monte Carlo Calculations of Excited States and Strongly Correlated Systems

S. Zhang

to develop quantum Monte Carlo (QMC) approaches for excited states and for more systematic and accurate treatment of ground states of strongly correlated sytems and heavier elements. We use the auxiliary-field quantum Monte Carlo (AFQMC) framework, which is a many-body total energy method that looks like many coupled mean-field calculations. I will discuss this method’s relations and dierences with traditional density-functional and quantum chemistry methods, and with other forms of QMC. The AFQMC method uses importance-sampled random walks in the space of nonorthogonal Slater determinants to project out the many-body ground state, or excited state. The random walks are constrained by a trial wave function (TWF) to remove the so-called sign problem, giving approximate (but non-perturbative) results. With singledeterminant TWFs from Hartree-Fock or density-functional theory (DFT), this approach can be thought of as a post-DFT method with no parameter. In this form, the accuracy is typically comparable to the quantum chemistry coupled-cluster CCSD(T) method near equilibrium geometries, and better than CCSD(T) when bonds are stretched. The computational cost of AFQMC scales as N 3 with characteristic size of the system. Its accuracy can be further improved at the cost of increased computing (releasing the constraint) or better TWF. With AFQMC, we can systematically target selected orbitals to correlate. This allows many-body calculations in solids with a frozen core, with down-folded Hamiltonians, or possibly embeded in a DFT calculations to extend scales. Results from a few on-going applications will be presented, including band structures in simple solids, magnetic phases in MnO, and Cobalt adsoption on graphene.

•Jan 1, 2012•0 citations

Electronic structure quantum Monte Carlo

Michal Bajdich, Luboš Mitáš

Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.

arXiv:1008.2369•Acta Physica Slovaca Reviews and Tutorials•Apr 1, 2009•33 citations

Variational block-correlated coupled-cluster Monte Carlo

Hunter Brodie, Alexander Shilcusky, Adrian E. Feiguin

Motivated by the recent success of quantum neural states, we present a variational implementation of the coupled-cluster method suitable for quantum Monte Carlo simulations. We introduce a Young tableaux representation of the wave-function coefficients and their derivatives that substantially reduces the computational complexity of the variational optimization. While the algorithm is generic, we illustrate its implementation with an application to the Heisenberg model in one and two dimensions, using the block-correlated coupled-cluster Ansatz with doubles and triples. Convergence is greatly improved by means of importance sampling based on a fictitious temperature. The resulting approach achieves competitive accuracy while providing a scalable framework for variational coupled-cluster calculations.

arXiv:2610.26646•The Journal of Chemical Physics•Oct 19, 2026•0 citations

Optimization of Ground- and Excited-State Wave Functions and van der Waals Clusters

M. P. Nightingale, V. Melik-Alaverdian

A quantum Monte Carlo method is introduced to optimize excited-state trial wave functions. The method is applied in a correlation function Monte Carlo calculation to compute ground- and excited-state energies of bosonic van der Waals clusters of up to seven particles. The calculations are performed using trial wave functions with general three-body correlations.

arXiv:physics/0010066•Physical Review Letters•Jul 9, 2001•104 citations

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
We introduce a many-body method that combines two powerful many-body techniques, viz., quantum Monte Carlo and coupled cluster theory. Coupled cluster wave functions are introduced as importance functions in a Monte Carlo method designed for the configuration interaction framework to provide rigorous upper bounds to the ground-state energy. We benchmark our method on the homogeneous electron gas in momentum space. The importance function used is the coupled cluster doubles wave function. We show that the computational resources required in our method scale polynomially with system size. Our energy upper bounds are in very good agreement with previous calculations of similar accuracy, and they can be systematically improved by including higher order excitations in the coupled cluster wave function.
The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.
to develop quantum Monte Carlo (QMC) approaches for excited states and for more systematic and accurate treatment of ground states of strongly correlated sytems and heavier elements. We use the auxiliary-field quantum Monte Carlo (AFQMC) framework, which is a many-body total energy method that looks like many coupled mean-field calculations. I will discuss this method’s relations and dierences with traditional density-functional and quantum chemistry methods, and with other forms of QMC. The AFQMC method uses importance-sampled random walks in the space of nonorthogonal Slater determinants to project out the many-body ground state, or excited state. The random walks are constrained by a trial wave function (TWF) to remove the so-called sign problem, giving approximate (but non-perturbative) results. With singledeterminant TWFs from Hartree-Fock or density-functional theory (DFT), this approach can be thought of as a post-DFT method with no parameter. In this form, the accuracy is typically comparable to the quantum chemistry coupled-cluster CCSD(T) method near equilibrium geometries, and better than CCSD(T) when bonds are stretched. The computational cost of AFQMC scales as N 3 with characteristic size of the system. Its accuracy can be further improved at the cost of increased computing (releasing the constraint) or better TWF. With AFQMC, we can systematically target selected orbitals to correlate. This allows many-body calculations in solids with a frozen core, with down-folded Hamiltonians, or possibly embeded in a DFT calculations to extend scales. Results from a few on-going applications will be presented, including band structures in simple solids, magnetic phases in MnO, and Cobalt adsoption on graphene.
Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.
Motivated by the recent success of quantum neural states, we present a variational implementation of the coupled-cluster method suitable for quantum Monte Carlo simulations. We introduce a Young tableaux representation of the wave-function coefficients and their derivatives that substantially reduces the computational complexity of the variational optimization. While the algorithm is generic, we illustrate its implementation with an application to the Heisenberg model in one and two dimensions, using the block-correlated coupled-cluster Ansatz with doubles and triples. Convergence is greatly improved by means of importance sampling based on a fictitious temperature. The resulting approach achieves competitive accuracy while providing a scalable framework for variational coupled-cluster calculations.
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.