The two-particle-hole Tamm-Dancoff approximation (2ph-TDA) equations for closed-shell atoms and molecules

Ab Initio Electron Propagators with an Hermitian, Intermediately Normalized Superoperator Metric Applied to Vertical Electron Affinities

New-generation ab initio electron propagator methods for calculating electron detachment energies of closed-shell molecules and anions have surpassed their predecessors' accuracy and computational efficiency. Derived from an Hermitian, intermediately normalized superoperator metric, these methods contain no adjustable parameters. To assess their versatility, a standard set (NIST-50-EA) of 50 vertical electron affinities of small closed-shell molecules based on NIST reference data has been created. Errors with respect to reference data on 23 large, conjugated organic photovoltaic (OPV23) molecules have also been analyzed. All final states are valence anions that correspond to electron affinities between 0.2 and 4.2 eV. For a given scaling of the arithmetic bottleneck, the new-generation methods obtain the lowest mean absolute errors (MAEs). The best methods with fifth-power arithmetic scaling realize MAEs below 0.1 eV. Composite models comprising cubically and quintically scaling calculations executed with large and small basis sets, respectively, produce OPV23 MAEs near 0.05 eV. The accuracy of quintically scaling methods executed with large basis sets is thereby procured with reduced computational effort. New-generation results obtained with and without the diagonal self-energy approximation in the canonical Hartree-Fock basis have been compared. These results indicate that Dyson orbitals closely resemble canonical Hartree-Fock orbitals multiplied by the square root of a probability factor above 0.85.

Reference Energies for Valence Ionizations and Satellite Transitions

Upon ionization of an atom or a molecule, another electron (or more) can be simultaneously excited. These concurrently generated states are called "satellites" (or shakeup transitions) as they appear in ionization spectra as higher-energy peaks with weaker intensity and larger width than the main peaks associated with single-particle ionizations. Satellites, which correspond to electronically excited states of the cationic species, are notoriously challenging to model using conventional single-reference methods due to their high excitation degree compared to the neutral reference state. This work reports 42 satellite transition energies and 58 valence ionization potentials (IPs) of full configuration interaction quality computed in small molecular systems. Following the protocol developed for the quest database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; and Loos, P.-F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], these reference energies are computed using the configuration interaction using a perturbative selection made iteratively (CIPSI) method. In addition, the accuracy of the well-known coupled-cluster (CC) hierarchy (CC2, CCSD, CC3, CCSDT, CC4, and CCSDTQ) is gauged against these new accurate references. The performances of various approximations based on many-body Green's functions (GW, GF2, and T-matrix) for IPs are also analyzed. Their limitations in correctly modeling satellite transitions are discussed.

New-Generation Electron-Propagator Methods for Molecular Electron-Binding Energies

A new generation of electron-propagator methods for the calculation of electron binding energies has surpassed its antecedents with respect to accuracy, efficiency, and interpretability. No adjustable parameters are introduced in these fully ab initio procedures. Numerical tests versus several databases of valence, vertical electron binding energies of closed-shell molecules and atoms have been performed. Easily interpreted self-energy approximations with cubic arithmetic scaling produce mean absolute errors (MAEs) of 0.2 and 0.3 eV for electron detachments and attachments, respectively. The most accurate explicitly renormalized methods with fifth-power arithmetic scaling yield MAEs below 0.1 eV for detachments and attachments. Approximate renormalization leads to more efficient fifth-power alternatives for electron detachments that achieve similar accuracy with fewer bottleneck operations. Composite protocols generate excellent predictions versus highly accurate basis-extrapolated standards and experiments. The validity of the diagonal self-energy approximation and the accuracy of the approximate renormalizations are confirmed. The success of these perturbative methods based on canonical Hartree-Fock orbitals rests on a Hermitized, intermediately normalized superoperator metric. The results of all of the new-generation calculations may be analyzed in terms of final-state orbital relaxation and differential correlation effects.

New-Generation Electron-Propagator Methods for Calculations of Electron Affinities and Ionization Energies: Tests on Organic Photovoltaic Molecules

A new generation of ab initio electron-propagator self-energies recently superseded its antecedents' accuracy and computational efficiency in calculating vertical ionization energies (VIEs) of closed-shell molecules. (See J. Chem. Phys. 2021, 155, 204107, J. Chem. Theory Comput. 2022, 18, 4927, J. Chem. Phys. 2023, 159, 124109.) No adjustable parameters were introduced in the generation of reference orbitals or in the construction of self-energies. The same approach has been extended in this work to vertical electron affinities (VEAs). Calculations were performed on 24 conjugated, organic photovoltaic molecules with diverse functional groups. These molecules are considerably larger than those studied in previous tests on VIEs. Several new-generation self-energies produce mean absolute errors (MAEs) below 0.1 eV versus ΔCCSD(T) (i.e., total energy differences from the coupled-cluster singles, doubles, and perturbative triples method) VIEs and VEAs obtained with identical basis sets. A composite model employs cubically and quintically scaling algorithms and power-law basis-set extrapolations based on augmented double-triple or triple-quadruple ζ data. Its MAEs are near 0.05 eV versus benchmark values, with 0.03 eV error bars for the lowest VIE and the highest VEA of each molecule. A more efficient and equally accurate composite model for calculating VIEs avoids full transformations of electron repulsion integrals to the molecular orbital basis. High probability factors support the diagonal self-energy approximation, wherein Dyson orbitals are proportional to canonical, Hartree-Fock orbitals.

A new generation of non-diagonal, renormalized self-energies for calculation of electron removal energies

A new generation of diagonal self-energies for the calculation of electron removal energies of molecules and molecular ions that has superseded its predecessors with respect to accuracy, efficiency, and interpretability is extended to include non-diagonal self-energies that permit Dyson orbitals to be expressed as linear combinations of canonical Hartree-Fock orbitals. In addition, an improved algorithm for renormalized methods eliminates the convergence difficulties encountered in the first studies of the new, diagonal self-energies. A dataset of outer-valence, vertical ionization energies with almost full-configuration-interaction quality serves as a standard of comparison in numerical tests. The new non-diagonal, renormalized methods are slightly more accurate than their diagonal counterparts, with mean absolute errors between 0.10 and 0.06 eV for outer-valence final states. This advantage is procured at the cost of an increase in the scaling of arithmetic bottlenecks that accompany the inclusion of non-diagonal self-energy terms. The new, non-diagonal, renormalized self-energies are also more accurate and efficient than their non-diagonal predecessors.

Electron Propagator Theory of Vertical Electron Detachment Energies of Anions: Benchmarks and Applications to Nucleotides

A new generation of ab initio electron-propagator self-energy approximations that are free of adjustable parameters is tested on a benchmark set of 55 vertical electron detachment energies of closed-shell anions. Comparisons with older self-energy approximations indicate that several new methods that make the diagonal self-energy approximation in the canonical Hartree-Fock orbital basis provide superior accuracy and computational efficiency. These methods and their acronyms, mean absolute errors (in eV), and arithmetic bottlenecks expressed in terms of occupied (O) and virtual (V) orbitals are the opposite-spin, non-Dyson, diagonal second-order method (os-nD-D2, 0.2, OV²), the approximately renormalized quasiparticle third-order method (Q3+, 0.15, O²V³) and the approximately renormalized, non-Dyson, linear, third-order method (nD-L3+, 0.1, OV⁴). The Brueckner doubles with triple field operators (BD-T1) nondiagonal electron-propagator method provides such close agreement with coupled-cluster single, double, and perturbative triple replacement total energy differences that it may be used as an alternative means of obtaining standard data. The new methods with diagonal self-energy matrices are the foundation of a composite procedure for estimating basis-set effects. This model produces accurate predictions and clear interpretations based on Dyson orbitals for the photoelectron spectra of the nucleotides found in DNA.

Electron Propagator Self-Energies versus Improved GW100 Vertical Ionization Energies

Ab initio electron propagator (EP) methods that are free of adjustable parameters in their self-energy formulae and in the generation of their orbital bases have been applied to the calculation of the lowest vertical ionization energies (VIEs) of the GW100 set. An improved set of standard results accompanied by irreducible representation assignments has been produced indirectly with coupled-cluster singles and doubles plus perturbative triples, i.e., CCSD(T), total energy differences at initial-state geometries reoptimized (in 28 cases) with the largest applicable point groups. The best compromises of accuracy and efficiency belong to a new generation of EP self-energies, several members of which may be derived from an intermediately normalized, Hermitized super-operator metric. The following diagonal self-energy methods are optimal: opposite-spin non-Dyson second order (os-nD-D2), approximately renormalized partial third order (P3+), approximately renormalized quasiparticle third order (Q3+), and non-Dyson approximately renormalized linear third order version B (nD-L3+B). Their mean absolute errors (MAEs) in electron volts and arithmetic scaling factors expressed in terms of occupied (O) and virtual (V) orbital dimensions are, respectively, (0.18, OV²), (0.14, O²V³), (0.15, O²V³), and (0.11, OV⁴). The 0.06 eV MAE for the non-diagonal, sixth-power (O²V⁴) Brueckner doubles, triple-field operator (BD-T1) EP method is exceeded by the 0.1 eV MAE with respect to experiments in seventh-power, ΔCCSD(T) calculations and indicates that BD-T1 may serve as a direct, spin-symmetry-conserving alternative in the generation of standard results for VIEs of larger, closed-shell molecules.

A new generation of diagonal self-energies for the calculation of electron removal energies

A new generation of diagonal self-energy approximations in ab initio electron propagator theory for the calculation of electron removal energies of molecules and molecular ions has been derived from an intermediately normalized, Hermitized super-operator metric. These methods and widely used antecedents such as the outer valence Green's function and the approximately renormalized partial third order method are tested with respect to a dataset of vertical ionization energies generated with a valence, triple-ζ, correlation-consistent basis set and a converged series of many-body calculations whose accuracy approaches that of full configuration interaction. Several modifications of the diagonal second-order self-energy, a version of G₀W₀ theory based on Tamm-Dancoff excitations and several non-diagonal self-energies are also included in the tests. All new methods employ canonical Hartree-Fock orbitals. No adjustable or empirical parameters appear. A hierarchy of methods with optimal accuracy for a given level of computational efficiency is established. Several widely used diagonal self-energy methods are rendered obsolete by the new hierarchy whose members, in order of increasing accuracy, are (1) the opposite-spin non-Dyson diagonal second-order or os-nD-D2, (2) the approximately renormalized third-order quasiparticle or Q3+, (3) the renormalized third-order quasiparticle or RQ3, (4) the approximately renormalized linear third-order or L3+, and (5) the renormalized linear third-order or RL3 self-energies.

A Fock space coupled cluster based probing of the single- and double-ionization profiles for the poly-cyclic aromatic hydrocarbons and conjugated polyenes

Sequential formation of a poly-cyclic aromatic hydrocarbon (PAH) dication in the H I regions of the interstellar medium (ISM) is proposed to be a function of internal energy of the doubly ionized PAHs, which, in turn, is dependent on the single- and double-ionization potentials of the system. This sets a limit on the single- and double-ionization energies of the system(s) that can further undergo sequential absorption of two photons, leading to a dication (PAH⁺²). Here, we report the single-ionization (I⁺¹) and double-ionization (I⁺²) energies and the I⁺²/I⁺¹ ratio for some selected PAHs and conjugated polyenes obtained using the Fock space coupled cluster technique, enabling simultaneous consideration of several electronic states of different characters. The I⁺² to I⁺¹ ratio bears a constant ratio, giving allowance to determine I⁺² from the knowledge of single-ionization (I⁺¹) and vice versa. Our observations are in good agreement with the established literature findings, confirming the reliability of our estimates. The measured single- and double-ionization energies further demonstrate that the sequential formation and fragmentation of a PAH dication in the H I regions of the ISM for systems such as benzene and conjugated polyenes such as ethylene and butadiene are quite unlikely because I⁺²-I⁺¹ for such system(s) is higher than the available photon energy in the H I regions of the ISM. Present findings may be useful to understand the formation and underlying decay mechanisms of multiply charged ions from PAHs and related compounds that may accentuate the exploration of the phenomenon of high-temperature superconductivity.

Full-frequency GW without frequency

Efficient computer implementations of the GW approximation must approximate a numerically challenging frequency integral; the integral can be performed analytically, but doing so leads to an expensive implementation whose computational cost scales as O(N6), where N is the size of the system. Here, we introduce a new formulation of the full-frequency GW approximation by exactly recasting it as an eigenvalue problem in an expanded space. This new formulation (1) avoids the use of time or frequency grids, (2) naturally obviates the need for the common "diagonal" approximation, (3) enables common iterative eigensolvers that reduce the canonical scaling to O(N5), and (4) enables a density-fitted implementation that reduces the scaling to O(N4). We numerically verify these scaling behaviors and test a variety of approximations that are motivated by this new formulation. The new formulation is found to be competitive with conventional O(N4) methods based on analytic continuation or contour deformation. In this new formulation, the relation of the GW approximation to configuration interaction, coupled-cluster theory, and the algebraic diagrammatic construction is made especially apparent, providing a new direction for improvements to the GW approximation.

The ionic states of cyclooctatetraene: Analysis of a new experimental photoelectron spectrum by ab initio and density functional methods

A synchrotron-based study of the photoelectron spectrum (PES) for cyclooctatetraene (COT) is reported, and this has been subjected to theoretical analysis in unprecedented detail. Weak vibrational structure was observed on the lowest ionization energy (IE₁), but the peaks generally show very broad features. Multiconfiguration self-consistent field study confirms that the adiabatic IE (AIE) sequence is 1²A₁ < 1²B₁ < 2²A₁ < 1²A₂ < 2²B₁ < 3²A₁. The Tamm-Dancoff approximation gives an acceptable interpretation of the PES below 20 eV. Vibrational analysis of the PES bands by Franck-Condon methods predicts well-defined vibrational structure for these ionic states. The principal contributors to the PES envelopes are the a₁ modes, and only a few are responsible for the overall shape of most bands. The high density of vibrational states, together with the known D_2d ⇆ D_2d interconversion process, where the C=C and C-C bonds interchange, is attributed to the lack of the observed structure. The transition state (TS) structures for the interconversion above of several ionic states of COT have been elucidated. The intrinsic reaction coordinate procedure gave a planar TS for the X¹A₁ (D_2d) ⇆ ¹A_1g (D_4h) ⇆ X¹A₁ (D_2d) process; this shows alternating C-C (1.4791 Å), C=C (1.3261 Å), and C-H (1.0780 Å) bond lengths. The planar TS is a very shallow maximum, with energy varying with the method used, Hartree-Fock 0.575 eV and second order Møller-Plesset 0.653 eV. A polynomial containing quadratic, quartic, and sextic terms gives an effectively exact fit to the surface. The lowest ionic state of COT (X²A₁) shows a similar TS for the process X²A₁ (D_2d) ⇆ ²B_2u (D_4h). This structure has alternating C-C (1.4366 Å), C=C (1.3572 Å), and C-H (1.0756 Å) bonds.

Green's function coupled cluster formulations utilizing extended inner excitations

In this paper, we analyze new approximations of the Green's function coupled cluster (GFCC) method where locations of poles are improved by extending the excitation level of inner auxiliary operators. These new GFCC approximations can be categorized as the GFCC-i(n, m) method, where the excitation level of the inner auxiliary operators (m) used to describe the ionization potential and electron affinity effects in the N - 1 and N + 1 particle spaces is higher than the excitation level (n) used to correlate the ground-state coupled cluster wave function for the N-electron system. Furthermore, we reveal the so-called "n + 1" rule in this category [or the GFCC-i(n, n + 1) method], which states that in order to maintain size-extensivity of the Green's function matrix elements, the excitation level of inner auxiliary operators X _p (ω) and Y _q (ω) cannot exceed n + 1. We also discuss the role of the moments of coupled cluster equations that in a natural way assures these properties. Our implementation in the present study is focused on the first approximation in this GFCC category, i.e., the GFCC-i(2,3) method. As our first practice, we use the GFCC-i(2,3) method to compute the spectral functions for the N₂ and CO molecules in the inner and outer valence regimes. In comparison with the Green's function coupled cluster singles, doubles results, the computed spectral functions from the GFCC-i(2,3) method exhibit better agreement with the experimental results and other theoretical results, particularly in terms of providing higher resolution of satellite peaks and more accurate relative positions of these satellite peaks with respect to the main peak positions.

Density functional theory study of the vertical ionization energies of the valence and core electrons of cyclopentadiene, pyrrole, furan, and thiophene

The procedure abbreviated as ΔPBE0(SAOP)/et-pVQZ, based on density functional theory, was developed recently for the calculation of vertical ionization energies of the valence electrons of organic and other small molecules and succeeded in giving results with an average absolute deviation of 0.21 eV from experiment for a collection of 115 reliable test cases of nonperhalo molecules. The objective of this work is to add a number of test cases to the benchmark database. We chose the set of molecules cyclo-C4H4X, with X = CH2, NH, O, and S, previously studied by many workers both experimentally and theoretically. The results show that the ΔPBE0(SAOP)/et-pVQZ procedure is not only as good as ab initio methods such as SAC-CI, OVGF, and ADC(3) in performance, but also handles inner valence ionized cations more efficiently. Although the core-electron binding energies of the titled molecules have not been as well investigated theoretically, we apply the methods we developed in recent years to calculate the binding energies of C1s, N1s, O1s, S1s, and S2p, which compare well with available experimental data.