The helium atom is the simplest many-body electronic system provided by nature. The exact solution to the Schrödinger equation is known for helium ground and excited states, and it represents a benchmark for any many-body methodology. Here, we check the ab initio many-body GW approximation and the Bethe-Salpeter equation (BSE) against the exact solution for helium. Starting from the Hartree-Fock method, we show that the GW and the BSE yield impressively accurate results on excitation energies and oscillator strength, systematically improving the time-dependent Hartree-Fock method. These findings suggest that the accuracy of the BSE and GW approximations is not significantly limited by self-interaction and self-screening problems even in this few electron limit. We further discuss our results in comparison to those obtained by time-dependent density-functional theory.
The helium atom is the simplest many-body electronic system provided by nature. The exact solution to the Schrödinger equation is known for helium ground and excited states, and it represents a benchmark for any many-body methodology. Here, we check the ab initio many-body GW approximation and the Bethe-Salpeter equation (BSE) against the exact solution for helium. Starting from the Hartree-Fock method, we show that the GW and the BSE yield impressively accurate results on excitation energies and oscillator strength, systematically improving the time-dependent Hartree-Fock method. These findings suggest that the accuracy of the BSE and GW approximations is not significantly limited by self-interaction and self-screening problems even in this few electron limit. We further discuss our results in comparison to those obtained by time-dependent density-functional theory.
Over time, many different theories and approaches have been developed to tackle the many-body problem in quantum chemistry, condensed-matter physics, and nuclear physics. Here we use the helium atom, a real system rather than a model, and we use the exact solution of its Schrödinger equation as a benchmark for comparison between methods. We present new results beyond the random-phase approximation (RPA) from a renormalized RPA (r-RPA) in the framework of the self-consistent RPA (SCRPA) originally developed in nuclear physics, and compare them with various other approaches like configuration interaction (CI), quantum Monte Carlo (QMC), time-dependent density-functional theory (TDDFT), and the Bethe-Salpeter equation on top of the GW approximation. Most of the calculations are consistently done on the same footing, e.g. using the same basis set, in an effort for a most faithful comparison between methods.
Over time, many different theories and approaches have been developed to tackle the many-body problem in quantum chemistry, condensed-matter physics, and nuclear physics. Here we use the helium atom, a real system rather than a model, and we use the exact solution of its Schrödinger equation as a benchmark for comparison between methods. We present new results beyond the random-phase approximation (RPA) from a renormalized RPA (r-RPA) in the framework of the self-consistent RPA (SCRPA) originally developed in nuclear physics, and compare them with various other approaches like configuration interaction (CI), quantum Monte Carlo (QMC), time-dependent density-functional theory (TDDFT), and the Bethe-Salpeter equation on top of the \boldsymbol{GW} 𝐆 𝐖 approximation. Most of the calculations are consistently done on the same footing, e.g. using the same basis set, in an effort for a most faithful comparison between methods.
The most of the experimental measurable quantities in many-body electron system, such as the optical absorption, luminescent spectra and exciton effect, are related to the excitation states description. The density functional theory with local density approximation (DFT-LDA) as a ground state theory, i.e. the solutions based on the Kohn-Sham equation is a powerful tool for studying the ground state properties in many-particle systems. However, the first principles computation of excited states is more complexity than ground-state calculations. A key problem is that the exchange-correlation interaction in excited states is differ from the ground-states. In recent years, nevertheless, several electronic-excitation theories have been developed. The most important theoretical and computational method include the many-body perturbation theory which is based on the quasi-particle concept and the Green function equation, the time-dependent DFT and the Bethe-Salpeter equation for describing the electron-hole interaction. Among them the central ingredient is an electron's self-energy operator Σ that describes the exchange correlation effect beyond Hartree approximation. The implementations of the above theories are unavoidable to introduce some of approximations. A good approximation for Σ is the GW approach by Hedin. It is shown from computation simulations for many real condensed matter systems that the GW approach is a successful method for electronic-excitation problems. Another effective method for excited-states is so-called screen-exchange local-density approximation (sX-LDA). It is a combination of LDA and Hartree-Fock (HF)theory, in which the intrinsic local screen-exchange interaction has been replaced by a non-local one and a generalized Kohn-Sham equation (GKS) is given. Based on a PWscf package, the sX-LDA-PWscf method has been tested and indicates that an available approach for treating band-gap problem in semiconductors. In this paper, the many-body theories for the excitation process, its development and signification are reviewed. The relation and differences consist in different theoretical methods are discussed. On the based of above discussions, the applications of many-body theory to the band-band transition (band-gap underestimation problem) and exciton effect in semiconductors are presented.
In a many-body perturbation treatment of electronic excitations, one has to solve the so-called Bethe–Salpeter equation (BSE) for the kernel (vertex) describing the interaction between electrons. In general, the BSE exhibits a non-local and frequency-dependent screened interaction and consequently it is extremely difficult to solve. We have developed a scheme that solves the BSE equation iteratively, including dynamically screening, by considering the change in the electron self-energy upon applying a time-dependent field. The BSE is obtained from the self-energy in the GW approximation (GWA) using the Schwinger functional technique. We apply our formalism to a simple model system and discuss briefly changes in the charge response, when dynamical screening is taken into account.
Here, we build on the works of Scuseria et al. [J. Chem. Phys. 129, 231101 (2008)] and Berkelbach [J. Chem. Phys. 149, 041103 (2018)] to show connections between the Bethe–Salpeter equation (BSE) formalism combined with the GW approximation from many-body perturbation theory and coupled-cluster (CC) theory at the ground- and excited-state levels. In particular, we show how to recast the GW and Bethe–Salpeter equations as non-linear CC-like equations. Similitudes between BSE@GW and the similarity-transformed equation-of-motion CC method are also put forward. The present work allows us to easily transfer key developments and the general knowledge gathered in CC theory to many-body perturbation theory. In particular, it may provide a path for the computation of ground- and excited-state properties (such as nuclear gradients) within the GW and BSE frameworks.
We calculate the excitation energies of finite 1D Hubbard chains with a variety of different site energies from two perspectives: (i) the physics-based Bethe-Salpeter equation (BSE) method and (ii) the chemistry-based configuration interaction (CI) approach. Results obtained from all methods are compared against the exact values for three classes of systems: metallic, impurity-doped, and molecular (semiconducting/insulating) systems. While in a previous study we showed that the GW method holds comparative advantages versus traditional quantum chemistry approaches for calculating the ionization potentials and electron affinities across a large range of Hamiltonians, we show now that the BSE method outperforms CI approaches only for metallic and semiconducting systems. For insulating molecular systems, CI approaches generate better results.