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

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.

New-generation electron-propagator methods for vertical electron detachment energies of molecular anions: benchmarks and applications to model green-fluorescent-protein chromophores

Ab initio electron-propagator calculations continue to be useful companions to experimental investigations of electronic structure in molecular anions. A new generation of electron-propagator methods recently has surpassed its antecedents' predictive accuracy and computational efficiency. Interpretive clarity has been conserved, for no adjustable parameters have been introduced in the preparation of molecular orbitals or in the formulation of approximate self-energies. These methods have employed the diagonal self-energy approximation wherein each Dyson orbital equals a canonical Hartree-Fock orbital times the square root of a probability factor. Numerical tests indicate that explicitly renormalized, diagonal self-energies are needed when Dyson orbitals have large valence nitrogen, oxygen or fluorine components. They also demonstrate that even greater accuracy can be realized with generalizations that do not employ the diagonal self-energy approximation in the canonical Hartree-Fock basis. Whereas the diagonal methods have fifth-power arithmetic scaling factors, the non-diagonal generalizations introduce only non-iterative sixth-power contractions. Composite models conserve the accuracy of the most demanding combinations of self-energy approximations and flexible basis sets with drastically reduced computational effort. Composite-model results on anions that resemble the chromophore of the green fluorescent protein illustrate the interpretive capabilities of explicitly renormalized self-energies. Accurate predictions on the lowest vertical electron detachment energy of each anion confirm experimental data and the utility of the diagonal self-energy approximation.

Achieving Accuracy and Economy for Calculating Vertical Detachment Energies of Molecular Anions: A Model Chemistry Composite Methods

The vertical detachment energy (VDE) is a vital factor for predicting the stability of anions that have important applications in the atom, molecule and cluster science. Due to the synthetic or characterization difficulty of anions, accurate and efficient predictions of VDE independent of laboratory data have always been an appealing task to remedy the experimental deficiencies. Unfortunately, the generally adopted CCSD(T) and electron propagator theory (EPT) methods have respectively been proven to be reliable but very cost-expensive, and cost-effective but sometimes problematic when Koopman's theorem is invalid. Here, we for the first time introduced and benchmarked a series of model chemistry composite methods (e. g., CBS-QB3, G4 and W1BD) on calculating VDE for 57 molecular anions. Notably, CBS-QB3 exceeds the accuracy of CCSD(T) while approaching the economy of EPT. Therefore, we highly recommend the composite method CBS-QB3 to compute VDEs for molecular anions in the attractive "killing two birds with one stone" manner.

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.

Fractional-Electron and Transition-Potential Methods for Core-to-Valence Excitation Energies Using Density Functional Theory

Methods for computing X-ray absorption spectra based on a constrained core hole (possibly containing a fractional electron) are examined. These methods are based on Slater's transition concept and its generalizations, wherein core-to-valence excitation energies are determined using Kohn-Sham orbital energies. Methods examined here avoid promoting electrons beyond the lowest unoccupied molecular orbital, facilitating robust convergence. Variants of these ideas are systematically tested, revealing a best-case accuracy of 0.3-0.4 eV (with respect to experiment) for K-edge transition energies. Absolute errors are much larger for higher-lying near-edge transitions but can be reduced below 1 eV by introducing an empirical shift based on a charge-neutral transition-potential method, in conjunction with functionals such as SCAN, SCAN0, or B3LYP. This procedure affords an entire excitation spectrum from a single fractional-electron calculation, at the cost of ground-state density functional theory and without the need for state-by-state calculations. This shifted transition-potential approach may be especially useful for simulating transient spectroscopies or in complex systems where excited-state Kohn-Sham calculations are challenging.

Vertical Ionization Potentials and Electron Affinities at the Double-Hybrid Density Functional Level

The double-hybrid (DH) time-dependent density functional theory is extended to vertical ionization potentials (VIPs) and electron affinities (VEAs). Utilizing the density fitting approximation, efficient implementations are presented for the genuine DH ansatz relying on the perturbative second-order correction, while an iterative analogue is also elaborated using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH approach. The favorable computational requirements of the present schemes are discussed in detail. The performance of the recently proposed spin-component-scaled and spin-opposite-scaled (SOS) range-separated (RS) and long-range corrected (LC) DH functionals is comprehensively assessed, while popular hybrid and global DH approaches are also discussed. For the benchmark calculations, up-to-date test sets are selected with high-level coupled-cluster references. Our results show that the ADC(2)-based SOS-RS-PBE-P86 approach is the most accurate and robust functional. This method consistently outperforms the excellent SOS-ADC(2) approach for VIPs, although the results are somewhat less satisfactory for VEAs. Among the genuine DH functionals, the SOS-ωPBEPP86 approach is also recommended for describing ionization processes, but its performance is even less reliable for electron-attached states. In addition, surprisingly good results are attained by the LC hybrid ωB97X-D functional, where the corresponding occupied (unoccupied) orbital energies are retrieved as VIPs (VEAs) within the present formalism.

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.

Accurate Prediction of Vertical Ionization Potentials and Electron Affinities from Spin-Component Scaled CC2 and ADC(2) Models

The CC2 and ADC(2) wave function models and their spin-component scaled modifications are adopted for predicting vertical ionization potentials (VIPs) and electron affinities (VEAs). The ionic solutions are obtained as electronic excitations in the continuum orbital formalism, making possible the use of existing, widespread quantum chemistry codes with minimal modifications, in full consistency with the treatment of charge transfer excitations. The performance of different variants is evaluated via benchmark calculations on various sets from previous works, containing small- and medium-sized systems, including the nucleobases. It is shown that with the spin-scaled approximate methods, in particular the scaled opposite-spin variant of the ADC(2) method, the accuracy of EOM-CCSD is achievable at a fraction of the computational cost, also outperforming many common electron propagator approaches.

Double Rydberg anions, Rydberg radicals and micro-solvated cations with ammonium-water kernels

Highly accurate ab initio electron-propagator and coupled-cluster methods are employed to predict the vertical electron attachment energies (VEAEs) of NH₄⁺(H₂O)_n (n = 1-4) cationic clusters. The VEAEs decrease with increasing n and the corresponding Dyson orbitals are diffused over peripheral, non-hydrogen bonded protons. Clusters formed from NH₄^- double Rydberg anions (DRAs) and stabilized by hydrogen bonding or electrostatic interactions are studied through calculations on NH₄^-(H₂O)_n complexes and are compared with more stable H^-(NH₃)(H₂O)_n isomers. Structures that have cationic and anionic congeners have notable changes in geometry. For all values of n, the hydride-molecule complex H^-(NH₃)(H₂O)_n is always the most stable, with large vertical electron detachment energies (VEDEs). NH₄^-(H₂O)_n DRA isomers are predicted to have VEDEs that correspond to energetically well-separated peaks in an anion photoelectron spectrum. Less stable DRA isomers display proton donation from the tetrahedral NH₄^- fragment to water molecules and VEDEs close to those of previously discovered DRAs. The most stable DRA isomers feature tetrahedral NH₄^- fragments without H bridges to water molecules and VEDEs that increase with n. Dyson orbitals of NH₄^-(H₂O)_n DRAs occupy regions beyond the exterior non-bridging O-H and N-H bonds. Thus, the Rydberg electrons in the uncharged Rydberg radicals and DRAs are held near the outer protons of the water and ammonia molecules. Several bound low-lying excited states of the doublet Rydberg radicals have single electrons occupying delocalized Dyson orbitals of s-like, p-like, d-like, or f-like nodal patterns with the following Aufbau principle: 1s, 1p, 1d, 2s, 2p, 1f.

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.

Multi-electron excitation contributions towards primary and satellite states in the photoelectron spectrum

The computation of Dyson orbitals and corresponding ionization energies has been implemented within the Equation of Motion Coupled Cluster Singles, Doubles and Perturbative Triples (EOM-CC3) method. Coupled to an accurate...