Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOn The Relation Between Equation-of-Motion Coupled-Cluster Theory and the GW Approximation
169
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We discuss the analytic and diagrammatic structure of ionization potential (IP) and electron affinity (EA) equation-of-motion coupled-cluster (EOM-CC) theory, in order to put it on equal footing with the prevalent $GW$ approximation. The comparison is most straightforward for the time-ordered one-particle Greens function, and we show that the Greens function calculated by EOM-CC with single and double excitations (EOM-CCSD) includes fewer ring diagrams at higher order than does the $GW$ approximation, due to the formers unbalanced treatment of time-ordering. However, the EOM-CCSD Greens function contains a large number of vertex corrections, including ladder diagrams, mixed ring-ladder diagrams, and exchange diagrams. By including triple excitations, the EOM-CCSDT Greens function includes all diagrams contained in the $GW$ approximation, along with many high-order vertex corrections. In the same language, we discuss a number of common approximations to the EOM-CCSD equations, many of which can be classified as elimination of diagrams. Finally, we present numerical results by calculating the principal charged excitations energies of the molecules contained in the so-called $GW$100 test set [J. Chem. Theory Comput. 2015, 11, 5665-5687]. We argue that (in molecules) exchange is as important as screening, advocating for a Hartree-Fock reference and second-order exchange in the self-energy.
rate research
Read More
We present ab initio absorption spectra of six three-dimensional semiconductors and insulators calculated using Gaussian-based periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). The spectra are calculated efficiently by solving a system of linear equations at each frequency, giving access to an energy range of tens of eV without explicit enumeration of excited states. We assess the impact of Brillouin zone sampling, for which it is hard to achieve convergence due to the cost of EOM-CCSD. Although our most converged spectra exhibit lineshapes that are in good agreement with experiment, they are uniformly shifted to higher energies by about 1 eV. We tentatively attribute this discrepancy to a combination of vibrational effects and the remaining electron correlation, i.e., triple excitations and above.
We study the decomposition of the Coulomb integrals of periodic systems into a tensor contraction of six matrices of which only two are distinct. We find that the Coulomb integrals can be well approximated in this form already with small matrices compared to the number of real space grid points. The cost of computing the matrices scales as O(N^4) using a regularized form of the alternating least squares algorithm. The studied factorization of the Coulomb integrals can be exploited to reduce the scaling of the computational cost of expensive tensor contractions appearing in the amplitude equations of coupled cluster methods with respect to system size. We apply the developed methodologies to calculate the adsorption energy of a single water molecule on a hexagonal boron nitride monolayer in a plane wave basis set and periodic boundary conditions.
Charged excitations of the oligoacene family of molecules, relevant for astrophysics and technological applications, are widely studied and therefore provide an excellent system for benchmarking theoretical methods. In this work, we evaluate the performance of many-body perturbation theory within the GW approximation relative to new high-quality CCSD(T) reference data for charged excitations of the acenes. We compare GW calculations with a number of hybrid density functional theory starting points and with eigenvalue self-consistency. Special focus is given to elucidating the trend of GW-predicted excitations with molecule length increasing from benzene to hexacene. We find that GW calculations with starting points based on an optimally tuned range-separated hybrid (OTRSH) density functional and eigenvalue self-consistency can yield quantitative ionization potentials for the acenes. However, for larger acenes, the predicted electron affinities can deviate considerably from reference values. Our work paves the way for predictive and cost-effective GW calculations of charged excitations of molecules and identifies certain limitations of current GW methods used in practice for larger molecules.
We describe a finite-field approach to compute density response functions, which allows for efficient $G_0W_0$ and $G_0W_0Gamma_0$ calculations beyond the random phase approximation. The method is easily applicable to density functional calculations performed with hybrid functionals. We present results for the electronic properties of molecules and solids and we discuss a general scheme to overcome slow convergence of quasiparticle energies obtained from $G_0W_0Gamma_0$ calculations, as a function of the basis set used to represent the dielectric matrix.
We present an equation of motion coupled cluster approach for calculating and understanding intrinsic inelastic losses in core level x-ray absorption spectra (XAS). The method is based on a factorization of the transition amplitude in the time-domain, which leads to a convolution of an effective one-body spectrum and the core-hole spectral function. The spectral function characterizes these losses in terms of shake-up excitations and satellites, and is calculated using a cumulant representation of the core-hole Greens function that includes non-linear corrections. The one-body spectrum also includes orthogonality corrections that enhance the XAS at the edge.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
