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

Error propagation in the fully self-consistent stochastic second-order Green's function method

Blair A. Winograd, Emanuel Gull, Dominika Zgid

Published Sep 21, 2020• 0 citations
arXivHTMLPDF

Abstract

We present an implementation of a fully self-consistent finite temperature second order Green's function perturbation theory (GF2) within the diagrammatic Monte Carlo framework. In contrast to the previous implementations of stochastic GF2 ({\it J. Chem. Phys.},{\bf 151}, 044144 (2019)), the current self-consistent stochastic GF2 does not introduce a systematic bias of the resulting electronic energies. Instead, the introduced implementation accounts for the stochastic errors appearing during the solution of the Dyson equation. We present an extensive discussion of the error handling necessary in a self-consistent procedure resulting in dressed Green's function lines. We test our method on a series of simple molecular examples.

DOI
10.48550/arxiv.2009.10048

Related papers

Stochastic self-consistent Green's function second-order perturbation theory (sGF2)

Daniel Neuhauser, Roi Baer, Dominika Zgid

The second-order Green's function method (GF2) was shown recently to be an accurate self-consistent approach for electronic structure of correlated systems since the self-energy accounts for both the weak and some of the strong correlation. The numerical scaling of GF2 is quite steep however, $O({N^5})$ (where the pre-factor is often hundreds), effectively preventing its application to large systems. Here, we develop a stochastic approach to GF2 (sGF2) where the self-energy is evaluated by a random-vector decomposition of Green's functions so that the dominant part of the calculation scales quasi linearly with system size. A study of hydrogen chains shows that the resulting approach is numerically efficient and accurate, as the stochastic errors are very small, 0.05% of the correlation energy for large systems with only a moderate computational effort. The method also yields automatically efficient MP2 energies and is automatically temperature dependent.

arXiv:1603.04141•Mar 14, 2016•3 citations

Monte Carlo Second- and Third-Order Many-Body Green’s Function Methods with Frequency-Dependent, Nondiagonal Self-Energy

Alexander E. Doran, So Hirata

We fully develop the Monte Carlo many-body Green’s function (MC-GF) method with the following enhancements: (1) The truncation order of the perturbation expansion of the Dyson self-energy is raised from the second order (MC-GF2) to the third order (MC-GF3) with the aid of a computerized procedure to enumerate and transform all 84 third-order Goldstone diagrams into Monte Carlo integrable expressions and then into central processing unit (CPU)/graphical processing unit (GPU)-parallel computer codes. (2) An efficient algorithm is proposed that computes all off-diagonal and diagonal elements of the MC-GF2 and MC-GF3 self-energy matrices by common subexpression elimination. (3) The frequency-independent approximation is lifted by introducing a method that computes frequency derivatives of the MC-GF2 and MC-GF3 self-energies up to any arbitrarily high order at nearly no additional computational cost. (4) The imaginary-time integration in the Laplace-transformed expressions of the self-energy is carried out stochastically (instead of using a quadrature in the previous implementations), resulting in a 50- to 200-fold speedup. (5) The efficiency of the redundant-walker convergence acceleration scheme is analyzed numerically, and the guidelines are established to select an optimal number of walkers for maximal efficiency. When such an optimal number is used, the cost per sample is constant of molecular size on either many CPUs or many GPUs. (6) The computational cost to obtain a binding energy within a given statistical uncertainty is observed to increase as (tentatively) O( n 4 ) and O( n 5 ) of molecular size ( n is the number of electrons) with and without the diagonal approximation, respectively, as opposed to O( n 5 ) and O( n 6 ) of the corresponding deterministic algorithms. With this method applied to the electron binding energies of C 60, we show that the third-order corrections to the self-energies are much greater in electron binding energies than in ground-state energies. They display a sign of oscillatory convergence toward experimental results, not necessarily improving the agreement with increasing perturbation order, justifying MC-GF3 and motivating even higher-order methods.

Journal of Chemical Theory and Computation•Oct 3, 2019•15 citations

Communication: The description of strong correlation within self-consistent Green's function second-order perturbation theory

Jordan J. Phillips, Dominika Zgid

We report an implementation of self-consistent Green's function many-body theory within a second-order approximation (GF2) for application with molecular systems. This is done by iterative solution of the Dyson equation expressed in matrix form in an atomic orbital basis, where the Green's function and self-energy are built on the imaginary frequency and imaginary time domain, respectively, and fast Fourier transform is used to efficiently transform these quantities as needed. We apply this method to several archetypical examples of strong correlation, such as a H32 finite lattice that displays a highly multireference electronic ground state even at equilibrium lattice spacing. In all cases, GF2 gives a physically meaningful description of the metal to insulator transition in these systems, without resorting to spin-symmetry breaking. Our results show that self-consistent Green's function many-body theory offers a viable route to describing strong correlations while remaining within a computationally tractable single-particle formalism.

arXiv:1405.3273•The Journal of Chemical Physics•Jun 25, 2014•123 citations

Stochastic Self-Consistent Second-Order Green’s Function Method for Correlation Energies of Large Electronic Systems

Daniel Neuhauser, Roi Baer, Dominika Zgid

The second-order Matsubara Green's function method (GF2) is a robust temperature-dependent quantum chemistry approach, extending beyond the random-phase approximation. However, until now the scope of GF2 applications was quite limited as they require computer resources that rise steeply with system size. In each step of the self-consistent GF2 calculation there are two parts: estimating of the self-energy from the previous step's Green's function and updating the Green's function from the self-energy. The first part formally scales as the fifth power of the system size, while the second has a much gentler cubic scaling. Here, we develop a stochastic approach to GF2 (sGF2), which reduces the fifth power scaling of the first step to merely quadratic, leaving the overall sGF2 scaling as cubic. We apply the method to linear hydrogen chains with up to 1000 electrons, showing that the approach is numerically stable, efficient, and accurate. The stochastic errors are very small, on the order of 0.1% or less of the correlation energy for large systems, with only a moderate computational effort. The first iteration of GF2 is an MP2 calculation that is done in linear scaling; hence we obtain an extremely fast stochastic MP2 (sMP2) method as a byproduct. While here we consider finite systems with large band gaps where at low temperatures effects are negligible, the sGF2 formalism is temperature dependent and general and can be applied to finite or periodic systems with small gaps at finite temperatures.

arXiv:1707.08296•Journal of Chemical Theory and Computation•Sep 29, 2017•54 citations

Self-consistent second-order Green’s function perturbation theory for periodic systems

Alexander A. Rusakov, Dominika Zgid

Despite recent advances, systematic quantitative treatment of the electron correlation problem in extended systems remains a formidable task. Systematically improvable Green's function methods capable of quantitatively describing weak and at least qualitatively strong correlations appear as promising candidates for computational treatment of periodic systems. We present a periodic implementation of temperature-dependent self-consistent 2nd-order Green's function (GF2) method, where the self-energy is evaluated in the basis of atomic orbitals. Evaluating the real-space self-energy in atomic orbitals and solving the Dyson equation in k-space are the key components of a computationally feasible algorithm. We apply this technique to the one-dimensional hydrogen lattice--a prototypical crystalline system with a realistic Hamiltonian. By analyzing the behavior of the spectral functions, natural occupations, and self-energies, we claim that GF2 is able to recover metallic, band insulating, and at least qualitatively Mott regimes. We observe that the iterative nature of GF2 is essential to the emergence of the metallic and Mott phases.

arXiv:1511.03911•The Journal of Chemical Physics•Feb 3, 2016•92 citations

Green/WeakCoupling: Implementation of fully self-consistent finite-temperature many-body perturbation theory for molecules and solids

Sergei Iskakov, Chia-Nan Yeh, Pavel Pokhilko +10 more

The accurate ab initio simulation of molecules and periodic solids with diagrammatic perturbation theory is an important task in quantum chemistry, condensed matter physics, and materials science. In this article, we present the WeakCoupling module of the open-source software package Green, which implements fully self-consistent diagrammatic weak coupling simulations, capable of dealing with real materials in the finite-temperature formalism. The code is licensed under the permissive MIT license. We provide self-consistent GW (scGW) and self-consistent second-order Green's function perturbation theory (GF2) solvers, analysis tools, and post-processing methods. This paper summarizes the theoretical methods implemented and provides background, tutorials and practical instructions for running simulations.

arXiv:2406.18479•Jun 26, 2024•1 citation
We present an implementation of a fully self-consistent finite temperature second order Green's function perturbation theory (GF2) within the diagrammatic Monte Carlo framework. In contrast to the previous implementations of stochastic GF2 ({\it J. Chem. Phys.},{\bf 151}, 044144 (2019)), the current self-consistent stochastic GF2 does not introduce a systematic bias of the resulting electronic energies. Instead, the introduced implementation accounts for the stochastic errors appearing during the solution of the Dyson equation. We present an extensive discussion of the error handling necessary in a self-consistent procedure resulting in dressed Green's function lines. We test our method on a series of simple molecular examples.
The second-order Green's function method (GF2) was shown recently to be an accurate self-consistent approach for electronic structure of correlated systems since the self-energy accounts for both the weak and some of the strong correlation. The numerical scaling of GF2 is quite steep however, O(N5)O({N^5})O(N5) (where the pre-factor is often hundreds), effectively preventing its application to large systems. Here, we develop a stochastic approach to GF2 (sGF2) where the self-energy is evaluated by a random-vector decomposition of Green's functions so that the dominant part of the calculation scales quasi linearly with system size. A study of hydrogen chains shows that the resulting approach is numerically efficient and accurate, as the stochastic errors are very small, 0.05% of the correlation energy for large systems with only a moderate computational effort. The method also yields automatically efficient MP2 energies and is automatically temperature dependent.
We fully develop the Monte Carlo many-body Green’s function (MC-GF) method with the following enhancements: (1) The truncation order of the perturbation expansion of the Dyson self-energy is raised from the second order (MC-GF2) to the third order (MC-GF3) with the aid of a computerized procedure to enumerate and transform all 84 third-order Goldstone diagrams into Monte Carlo integrable expressions and then into central processing unit (CPU)/graphical processing unit (GPU)-parallel computer codes. (2) An efficient algorithm is proposed that computes all off-diagonal and diagonal elements of the MC-GF2 and MC-GF3 self-energy matrices by common subexpression elimination. (3) The frequency-independent approximation is lifted by introducing a method that computes frequency derivatives of the MC-GF2 and MC-GF3 self-energies up to any arbitrarily high order at nearly no additional computational cost. (4) The imaginary-time integration in the Laplace-transformed expressions of the self-energy is carried out stochastically (instead of using a quadrature in the previous implementations), resulting in a 50- to 200-fold speedup. (5) The efficiency of the redundant-walker convergence acceleration scheme is analyzed numerically, and the guidelines are established to select an optimal number of walkers for maximal efficiency. When such an optimal number is used, the cost per sample is constant of molecular size on either many CPUs or many GPUs. (6) The computational cost to obtain a binding energy within a given statistical uncertainty is observed to increase as (tentatively) O( n 4 ) and O( n 5 ) of molecular size ( n is the number of electrons) with and without the diagonal approximation, respectively, as opposed to O( n 5 ) and O( n 6 ) of the corresponding deterministic algorithms. With this method applied to the electron binding energies of C 60, we show that the third-order corrections to the self-energies are much greater in electron binding energies than in ground-state energies. They display a sign of oscillatory convergence toward experimental results, not necessarily improving the agreement with increasing perturbation order, justifying MC-GF3 and motivating even higher-order methods.
We report an implementation of self-consistent Green's function many-body theory within a second-order approximation (GF2) for application with molecular systems. This is done by iterative solution of the Dyson equation expressed in matrix form in an atomic orbital basis, where the Green's function and self-energy are built on the imaginary frequency and imaginary time domain, respectively, and fast Fourier transform is used to efficiently transform these quantities as needed. We apply this method to several archetypical examples of strong correlation, such as a H32 finite lattice that displays a highly multireference electronic ground state even at equilibrium lattice spacing. In all cases, GF2 gives a physically meaningful description of the metal to insulator transition in these systems, without resorting to spin-symmetry breaking. Our results show that self-consistent Green's function many-body theory offers a viable route to describing strong correlations while remaining within a computationally tractable single-particle formalism.
The second-order Matsubara Green's function method (GF2) is a robust temperature-dependent quantum chemistry approach, extending beyond the random-phase approximation. However, until now the scope of GF2 applications was quite limited as they require computer resources that rise steeply with system size. In each step of the self-consistent GF2 calculation there are two parts: estimating of the self-energy from the previous step's Green's function and updating the Green's function from the self-energy. The first part formally scales as the fifth power of the system size, while the second has a much gentler cubic scaling. Here, we develop a stochastic approach to GF2 (sGF2), which reduces the fifth power scaling of the first step to merely quadratic, leaving the overall sGF2 scaling as cubic. We apply the method to linear hydrogen chains with up to 1000 electrons, showing that the approach is numerically stable, efficient, and accurate. The stochastic errors are very small, on the order of 0.1% or less of the correlation energy for large systems, with only a moderate computational effort. The first iteration of GF2 is an MP2 calculation that is done in linear scaling; hence we obtain an extremely fast stochastic MP2 (sMP2) method as a byproduct. While here we consider finite systems with large band gaps where at low temperatures effects are negligible, the sGF2 formalism is temperature dependent and general and can be applied to finite or periodic systems with small gaps at finite temperatures.
Despite recent advances, systematic quantitative treatment of the electron correlation problem in extended systems remains a formidable task. Systematically improvable Green's function methods capable of quantitatively describing weak and at least qualitatively strong correlations appear as promising candidates for computational treatment of periodic systems. We present a periodic implementation of temperature-dependent self-consistent 2nd-order Green's function (GF2) method, where the self-energy is evaluated in the basis of atomic orbitals. Evaluating the real-space self-energy in atomic orbitals and solving the Dyson equation in k-space are the key components of a computationally feasible algorithm. We apply this technique to the one-dimensional hydrogen lattice--a prototypical crystalline system with a realistic Hamiltonian. By analyzing the behavior of the spectral functions, natural occupations, and self-energies, we claim that GF2 is able to recover metallic, band insulating, and at least qualitatively Mott regimes. We observe that the iterative nature of GF2 is essential to the emergence of the metallic and Mott phases.
The accurate ab initio simulation of molecules and periodic solids with diagrammatic perturbation theory is an important task in quantum chemistry, condensed matter physics, and materials science. In this article, we present the WeakCoupling module of the open-source software package Green, which implements fully self-consistent diagrammatic weak coupling simulations, capable of dealing with real materials in the finite-temperature formalism. The code is licensed under the permissive MIT license. We provide self-consistent GW (scGW) and self-consistent second-order Green's function perturbation theory (GF2) solvers, analysis tools, and post-processing methods. This paper summarizes the theoretical methods implemented and provides background, tutorials and practical instructions for running simulations.