Coupled cluster Green function: Model involving single and double excitations
TL;DR: A key aspect of this work is the determination of the frequency dependent self-energy, Σ(ω), which reveals a diagonally dominate block structure where as the non-local correlation increases, the block size increases proportionally.
Abstract: In this paper, we report on the development of a parallel implementation of the coupled-cluster (CC) Green function formulation (GFCC) employing single and double excitations in the cluster operator (GFCCSD). A key aspect of this work is the determination of the frequency dependent self-energy, Σ(ω). The detailed description of the underlying algorithm is provided, including approximations used that preserve the pole structure of the full GFCCSD method, thereby reducing the computational costs while maintaining an accurate character of methodology. Furthermore, for systems with strong local correlation, our formulation reveals a diagonally dominate block structure where as the non-local correlation increases, the block size increases proportionally. To demonstrate the accuracy of our approach, several examples including calculations of ionization potentials for benchmark systems are presented and compared against experiment.
Citations
More filters
Journal Article•
TL;DR: Vibrational progressions for the S(0)-D(0) vibronic bands computed within double-harmonic approximation with Duschinsky rotations are compared with previously reported experimental photoelectron spectra and differentiated PIE curves.
Abstract: Electronic structure and spectroscopy of nucleic acid bases: Ionization energies, ionization-induced structural changes, and photoelectron spectra Ksenia B. Bravaya a , Oleg Kostko b , Stanislav Dolgikh a , Arie Landau a , Musahid Ahmed b , and Anna I. Krylov a a Department of Chemistry, University of Southern California, Los Angeles, CA 90089-0482, USA b Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA We report high-level ab initio calculations and single-photon ionization mass spec- trometry study of ionization of adenine (A), thymine (T), cytosine (C) and guanine (G). For thymine and adenine, only the lowest-energy tautomers were considered, whereas for cytosine and guanine we characterized five lowest-energy tautomeric forms. The first adiabatic and several vertical ionization energies were computed us- ing equation-of-motion coupled-cluster method for ionization potentials with single and double substitutions. Equilibrium structures of the cationic ground states were characterized by DFT with the ωB97X-D functional. The ionization-induced geome- try changes of the bases are consistent with the shapes of the corresponding molecular orbitals. For the lowest-energy tautomers, the magnitude of the structural relaxation decreases in the following series G>C>A>T, the respective relaxation energies being 0.41, 0.32, 0.25 and 0.20 eV. The computed adiabatic ionization energies (8.13, 8.89, 8.51-8.67 and 7.75-7.87 eV for A,T,C and G, respectively) agree well with the on- sets of the photoionization efficiency (PIE) curves (8.20±0.05, 8.95±0.05, 8.60±0.05 and 7.75±0.05 eV). Vibrational progressions for the S 0 -D 0 vibronic bands computed within double-harmonic approximation with Duschinsky rotations are compared with previously reported experimental photoelectron spectra. I. INTRODUCTION Ionization of DNA bases is one of the processes involved in DNA radiation- and photo- damage resulting in dangerous mutations with potential risk for cancer and neurodegenera-
79 citations
TL;DR: It is argued that (in molecules) exchange is as important as screening, advocating for a Hartree-Fock reference and second-order exchange in the self-energy.
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 Green’s function, and we show that the Green’s 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 former’s unbalanced treatment of time-ordering. However, the EOM-CCSD Green’s 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 Green’s 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 classifie...
68 citations
TL;DR: Hedin's scheme is solved with the inclusion of the vertex function (GWΓ) for a set of small molecules, and very good agreement with expensive wave function-based methods is obtained for both approximations.
Abstract: Hedin's scheme is solved with the inclusion of the vertex function (GWΓ) for a set of small molecules. The computational scheme allows for the consistent inclusion of the vertex both at the polarizability level and in the self-energy. A diagrammatic analysis shows that the self-energy formed with this four-point vertex does not lead to double counting of diagrams, that can be classied as direct "bubbles" and exchange diagrams. By removing the exchange diagrams from the self-energy, a simpler approximation is obtained, called GWtc-tc. Very good agreement with expensive wavefunction-based methods is obtained for both approximations.
55 citations
TL;DR: The diagrammatic linkedness and thus size-consistency of the one-particle Green's function and self-energy are demonstrated at any perturbation order on the basis of the algebraic recursions in an entirely time-independent (frequency-domain) framework.
Abstract: A thorough analytical and numerical characterization of the whole perturbation series of one-particle many-body Green’s function (MBGF) theory is presented in a pedagogical manner. Three distinct but equivalent algebraic (first-quantized) recursive definitions of the perturbation series of the Green’s function are derived, which can be combined with the well-known recursion for the self-energy. Six general-order algorithms of MBGF are developed, each implementing one of the three recursions, the ΔMPn method (where n is the perturbation order) [S. Hirata et al., J. Chem. Theory Comput. 11, 1595 (2015)], the automatic generation and interpretation of diagrams, or the numerical differentiation of the exact Green’s function with a perturbation-scaled Hamiltonian. They all display the identical, nondivergent perturbation series except ΔMPn, which agrees with MBGF in the diagonal and frequency-independent approximations at 1≤n≤3 but converges at the full-configuration-interaction (FCI) limit at n=∞ (unless it diverges). Numerical data of the perturbation series are presented for Koopmans and non-Koopmans states to quantify the rate of convergence towards the FCI limit and the impact of the diagonal, frequency-independent, or ΔMPn approximation. The diagrammatic linkedness and thus size-consistency of the one-particle Green’s function and self-energy are demonstrated at any perturbation order on the basis of the algebraic recursions in an entirely time-independent (frequency-domain) framework. The trimming of external lines in a one-particle Green’s function to expose a self-energy diagram and the removal of reducible diagrams are also justified mathematically using the factorization theorem of Frantz and Mills. Equivalence of ΔMPn and MBGF in the diagonal and frequency-independent approximations at 1≤n≤3 is algebraically proven, also ascribing the differences at n = 4 to the so-called semi-reducible and linked-disconnected diagrams.
55 citations
TL;DR: The unusual intrinsic electronic structure of the [B12X12]2- MCAs provides the basis for a molecular level understanding of their observed unique physical and chemical properties and a new paradigm for understanding the properties of these MCAs with not well-separated charges that departs from the prevailing model used for spatially separated charges.
Abstract: The stability and electron loss process of numerous multiply charged anions (MCAs) have been traditionally explained in terms of the classical Coulomb interaction between spatially separated charged groups. An understanding of these processes in MCAs with not well-separated excess charges is still lacking. We report the surprising properties and physical behavior of [B12X12]2–, X = F, Cl, Br, I, At, which are MCAs with not well-separated excess charges and cannot be described by the prevailing classical picture. In this series of MCAs, comprising a “boron core” surrounded by a “halogen shell”, the sign of the total charge in these two regions changes along the halogen series from X = F–At. With the aid of experimental photoelectron spectroscopy and highly correlated ab initio electronic structure calculations, we demonstrate that the trend in the electronic stability of these MCAs is determined by the interplay between the Coulomb (de)stabilization originating from the “boron core” and “halogen shell” and...
51 citations
References
More filters
TL;DR: In this paper, a detailed study of correlation effects in the oxygen atom was conducted, and it was shown that primitive basis sets of primitive Gaussian functions effectively and efficiently describe correlation effects.
Abstract: In the past, basis sets for use in correlated molecular calculations have largely been taken from single configuration calculations. Recently, Almlof, Taylor, and co‐workers have found that basis sets of natural orbitals derived from correlated atomic calculations (ANOs) provide an excellent description of molecular correlation effects. We report here a careful study of correlation effects in the oxygen atom, establishing that compact sets of primitive Gaussian functions effectively and efficiently describe correlation effects i f the exponents of the functions are optimized in atomic correlated calculations, although the primitive (s p) functions for describing correlation effects can be taken from atomic Hartree–Fock calculations i f the appropriate primitive set is used. Test calculations on oxygen‐containing molecules indicate that these primitive basis sets describe molecular correlation effects as well as the ANO sets of Almlof and Taylor. Guided by the calculations on oxygen, basis sets for use in correlated atomic and molecular calculations were developed for all of the first row atoms from boron through neon and for hydrogen. As in the oxygen atom calculations, it was found that the incremental energy lowerings due to the addition of correlating functions fall into distinct groups. This leads to the concept of c o r r e l a t i o n c o n s i s t e n t b a s i s s e t s, i.e., sets which include all functions in a given group as well as all functions in any higher groups. Correlation consistent sets are given for all of the atoms considered. The most accurate sets determined in this way, [5s4p3d2f1g], consistently yield 99% of the correlation energy obtained with the corresponding ANO sets, even though the latter contains 50% more primitive functions and twice as many primitive polarization functions. It is estimated that this set yields 94%–97% of the total (HF+1+2) correlation energy for the atoms neon through boron.
26,705 citations
TL;DR: This review gives a general introduction to the materials, production techniques, working principles, critical parameters, and stability of the organic solar cells, and discusses the alternative approaches such as polymer/polymer solar cells and organic/inorganic hybrid solar cells.
Abstract: The need to develop inexpensive renewable energy sources stimulates scientific research for efficient, low-cost photovoltaic devices.1 The organic, polymer-based photovoltaic elements have introduced at least the potential of obtaining cheap and easy methods to produce energy from light.2 The possibility of chemically manipulating the material properties of polymers (plastics) combined with a variety of easy and cheap processing techniques has made polymer-based materials present in almost every aspect of modern society.3 Organic semiconductors have several advantages: (a) lowcost synthesis, and (b) easy manufacture of thin film devices by vacuum evaporation/sublimation or solution cast or printing technologies. Furthermore, organic semiconductor thin films may show high absorption coefficients4 exceeding 105 cm-1, which makes them good chromophores for optoelectronic applications. The electronic band gap of organic semiconductors can be engineered by chemical synthesis for simple color changing of light emitting diodes (LEDs).5 Charge carrier mobilities as high as 10 cm2/V‚s6 made them competitive with amorphous silicon.7 This review is organized as follows. In the first part, we will give a general introduction to the materials, production techniques, working principles, critical parameters, and stability of the organic solar cells. In the second part, we will focus on conjugated polymer/fullerene bulk heterojunction solar cells, mainly on polyphenylenevinylene (PPV) derivatives/(1-(3-methoxycarbonyl) propyl-1-phenyl[6,6]C61) (PCBM) fullerene derivatives and poly(3-hexylthiophene) (P3HT)/PCBM systems. In the third part, we will discuss the alternative approaches such as polymer/polymer solar cells and organic/inorganic hybrid solar cells. In the fourth part, we will suggest possible routes for further improvements and finish with some conclusions. The different papers mentioned in the text have been chosen for didactical purposes and cannot reflect the chronology of the research field nor have a claim of completeness. The further interested reader is referred to the vast amount of quality papers published in this field during the past decade.
6,059 citations
TL;DR: The coupled cluster singles and doubles model (CCSD) as discussed by the authors is derived algebraically, presenting the full set of equations for a general reference function explicitly in spin-orbital form, and the computational implementation of the CCSD model, which involves cubic and quartic terms, is discussed and results are compared with full CI calculations for H2O and BeH2.
Abstract: The coupled‐cluster singles and doubles model (CCSD) is derived algebraically, presenting the full set of equations for a general reference function explicitly in spin–orbital form. The computational implementation of the CCSD model, which involves cubic and quartic terms, is discussed and results are reported and compared with full CI calculations for H2O and BeH2. We demonstrate that the CCSD exponential ansatz sums higher‐order correlation effects efficiently even for BeH2, near its transition state geometry where quasidegeneracy efforts are quite large, recovering 98% of the full CI correlation energy. For H2O, CCSD plus the fourth‐order triple excitation correction agrees with the full CI energy to 0.5 kcal/mol. Comparisons with low‐order models provide estimates of the effect of the higher‐order terms T1T2, T21T2, T31, and T41 on the correlation energy.
5,603 citations
TL;DR: An overview of NWChem is provided focusing primarily on the core theoretical modules provided by the code and their parallel performance, as well as Scalable parallel implementations and modular software design enable efficient utilization of current computational architectures.
4,666 citations