Do Dyson Orbitals resemble canonical Hartree–Fock orbitals?

Why Does the GW Approximation Give Accurate Quasiparticle Energies? The Cancellation of Vertex Corrections Quantified

Hedin's GW approximation to the electronic self-energy has been impressively successful in calculating quasiparticle energies, such as ionization potentials, electron affinities, or electronic band structures. The success of this fairly simple approximation has been ascribed to the cancellation of the so-called vertex corrections that go beyond the GW approximation. This claim is mostly based on past calculations using vertex corrections within the crude local-density approximation. Here, we explore a wide variety of nonlocal vertex corrections in the polarizability and the self-energy, using first-order approximations or infinite summations to all orders. In particular, we use vertices based on statically screened interactions like in the Bethe-Salpeter equation. We demonstrate on realistic molecular systems that the two vertices in Hedin's equation essentially compensate. We further show that consistency between the two vertices is crucial for obtaining realistic electronic properties. We finally consider increasingly large clusters and extrapolate that our conclusions about the compensation of the two vertices would hold for extended systems.

Time-dependent ab initio molecular-orbital decomposition for high-harmonic generation spectroscopy

We propose a real-time time-dependent ab initio approach within a configuration-interaction-singles ansatz to decompose the high-harmonic generation (HHG) signal of molecules in terms of individual molecular-orbital (MO) contributions. Calculations have been performed by propagating the time-dependent Schrödinger equation with complex energies, in order to account for ionization of the system, and by using tailored Gaussian basis sets for high-energy and continuum states. We have studied the strong-field electron dynamics and the HHG spectra in aligned CO2 and H2O molecules. Contribution from MOs in the strong-field dynamics depends on the interplay between the MO ionization energy and the coupling between the MO and the laser-pulse symmetries. Such contributions characterize different portions of the HHG spectrum, indicating that the orbital decomposition encodes nontrivial information on the modulation of the strong-field dynamics. Our results correctly reproduce the MO contributions to HHG for CO2 as described in the literature experimental and theoretical data and lead to an original analysis of the role of the highest occupied molecular orbitals HOMO, HOMO-1, and HOMO-2 of H2O according to the polarization direction of the laser pulse.

Temporary Anions of Benzoxazole in Charge-Transfer Cluster Photodetachment

The photoelectron spectra of cluster anions of superoxide (O2-) solvated by one molecule of benzoxazole (BzOx) reveal two competing photodetachment mechanisms: a direct photoemission from the solvated cluster core and an indirect pathway involving temporary anion states of benzoxazole accessed via the O2-·BzOx → O2·BzOx- charge-transfer transitions. Benzoxazole is a bicyclic unsaturated organic molecule that does not form permanent anions. However, its low-lying vacant π* orbitals permit a resonant capture of the electron emitted from the O2- cluster core. The non-Hermitian theory using a complex absorbing potential predicts the existence of two BzOx- π* resonances within the experimental energy range: resonance A (π1*) at 0.891 eV and resonance B (π2*) at 1.76 eV, relative to the onset of the BzOx + e- continuum at the ground-state geometry of neutral BzOx. Within the clusters, the O2·BzOx- charge-transfer states are partially stabilized relative to the free-electron limit by interactions with the O2 molecule. These interactions depend on the electronic states of both species. The theory predicts that at the O2-·BzOx cluster geometry, the O2(X3Σg-)·BzOx-(A) and O2(a1Δg)·BzOx-(A) states lie at 0.56 and 0.47 eV (vertically) above the respective neutral states. The O2(3Σg-)·BzOx-(B) resonance is found 1.43 eV (vertically) above O2(X3Σg-)·BzOx. Intense signatures of both BzOx- resonances and the three above-mentioned charge-transfer cluster states, O2(X3Σg-)·BzOx-(A), O2(a1Δg)·BzOx-(A), and O2(3Σg-)·BzOx-(B) are observed in the 355 nm (3.495 eV) and 532 nm (2.330 eV) photoelectron spectra of the O2-·BzOx cluster anion.

Electron Binding Energies of Open-Shell Species from Diagonal Electron-Propagator Self-Energies with Unrestricted Hartree-Fock Spin-Orbitals

For closed-shell molecules, valence electron binding energies may be calculated accurately and efficiently with ab initio electron-propagator methods that have surpassed their predecessors. Advantageous combinations of accuracy and efficiency range from cubically scaling methods with mean errors of 0.2 eV to quintically scaling methods with mean errors of 0.1 eV or less. The diagonal self-energy approximation in the canonical Hartree-Fock basis is responsible for the enhanced efficiency of several methods. This work examines the predictive capabilities of diagonal self-energy approximations when they are generalized to the canonical spin-orbital basis of unrestricted Hartree-Fock (UHF) theory. Experimental data on atomic electron binding energies and high-level, correlated calculations in a fixed basis for a set of open-shell molecules constitute standards of comparison. A review of the underlying theory and analysis of numerical errors lead to several recommendations for the calculation of electron binding energies: (1) In calculations that employ the diagonal self-energy approximation, Koopmans's identity for UHF must be qualitatively correct. (2) Closed-shell reference states are preferable to open-shell reference states in calculations of molecular ionization energies and electron affinities. (3) For molecular electron binding energies between doublets and triplets, calculations of electron detachment energies are more accurate and efficient than calculations of electron attachment energies. When these recommendations are followed, mean absolute errors increase by approximately 0.05 eV with respect to their counterparts obtained with closed-shell reference orbitals.

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.

Anomalous π-backbonding in complexes between B(SiR3)3 and N2: catalytic activation and breaking of scaling relations

Chemical transformations of molecular nitrogen (N2), including the nitrogen reduction reaction (NRR), are difficult to catalyze because of the weak Lewis basicity of N2. In this study, it is shown that Lewis acids of the types B(SiR3)3 and B(GeR3)3 bind N2 and CO with anomalously short and strong B-N or B-C bonds. B(SiH3)3·N2 has a B-N bond length of 1.48 Å and a complexation enthalpy of -15.9 kcal mol-1 at the M06-2X/jun-cc-pVTZ level. The selective binding enhancement of N2 and CO is due to π-backbonding from Lewis acid to Lewis base, as demonstrated by orbital analysis and density difference plots. The π-backbonding is found to be a consequence of constructive orbital interactions between the diffuse and highly polarizable B-Si and B-Ge bond regions and the π and π* orbitals of N2. This interaction is strengthened by electron donating substituents on Si or Ge. The π-backbonding interaction is predicted to activate N2 for chemical transformation and reduction, as it decreases the electron density and increases the length of the N-N bond. The binding of N2 and CO by the B(SiR3)3 and B(GeR3)3 types of Lewis acids also has a strong σ-bonding contribution. The relatively high σ-bond strength is connected to the highly positive surface electrostatic potential [VS(r)] above the B atom in the tetragonal binding conformation, but the σ-bonding also has a significant coordinate covalent (dative) contribution. Electron withdrawing substituents increase the potential and the σ-bond strength, but favor the binding of regular Lewis acids, such as NH3 and F-, more strongly than binding of N2 and CO. Molecules of the types B(SiR3)3 and B(GeR3)3 are chemically labile and difficult to synthesize. Heterogenous catalysts with the wanted B(Si-)3 or B(Ge-)3 bonding motif may be prepared by boron doping of nanostructured silicon or germanium compounds. B-doped and hydrogenated silicene is found to have promising properties as catalyst for the electrochemical NRR.

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.

Dipole effects in the photoelectron angular distributions of the sulfur monoxide anion

Photoelectron angular distributions (PADs) in SO^- photodetachment using linearly polarized 355 nm (3.49 eV), 532 nm (2.33 eV), and 611 nm (2.03 eV) light were investigated via photoelectron imaging spectroscopy. The measurements at 532 and 611 nm access the X³Σ^- and a¹Δ electronic states of SO, whereas the measurements at 355 nm also access the b¹Σ⁺ state. In aggregate, the photoelectron anisotropy parameter values follow the general trend with respect to electron kinetic energy (eKE) expected for π*-orbital photodetachment. The trend is similar to O₂^-, but the minimum of the SO^- curve is shifted to smaller eKE. This shift is mainly attributed to the exit-channel interactions of the departing electron with the dipole moment of the neutral SO core, rather than the differing shapes of the SO^- and O₂^- molecular orbitals. Of the several ab initio models considered, two approaches yield good agreement with the experiment: one representing the departing electron as a superposition of eigenfunctions of a point dipole-field Hamiltonian, and another describing the outgoing electron in terms of Coulomb waves originating from two separated charge centers, with a partial positive charge on the sulfur and an equal negative charge on the oxygen. These fundamentally related approaches support the conclusion that electron-dipole interactions in the exit channel of SO^- photodetachment play an important role in shaping the PADs. While a similar conclusion was previously reached for photodetachment from σ orbitals of CN^- (Hart, Lyle, Spellberg, Krylov, Mabbs, J. Phys. Chem. Lett., 2021, 12, 10086-10092), the present work includes the first extension of the dipole-field model to detachment from π* orbitals.

Assessing the performance of ΔSCF and the diagonal second-order self-energy approximation for calculating vertical core excitation energies

Vertical core excitation energies are obtained using a combination of the ΔSCF method and the diagonal second-order (D2) self-energy approximation. These methods are applied to a set of neutral molecules and their anionic forms. An assessment of the results with the inclusion of relativistic effects is presented. For core excitations involving delocalized symmetry orbitals, the applied composite method improves upon the overestimation of ΔSCF by providing approximate values close to experimental K-shell transition energies. The importance of both correlation and relaxation contributions to the vertical core-excited state energies, the concept of local and non-local core orbitals, and the consequences of breaking symmetry are discussed.

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...

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.

Exploring Dyson's Orbitals and Their Electron Binding Energies for Conceptualizing Excited States from Response Methodology

The molecular orbital (MO) concept is a useful tool, which relates the molecular ground-state energy with the energies (and occupations) of the individual orbitals. However, analysis of the excited states from linear response computations is performed in terms of the initial state MOs or some other forms of orbitals, e.g., natural or natural transition orbitals. Because these orbitals lack the respective energies, they do not allow developing a consistent orbital picture of the excited states. Herein, we argue that Dyson's orbitals enable description of the response states compatible with the concepts of molecular orbital theory. The Dyson orbitals and their energies obtained by mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) for the response ground state are remarkably similar to the canonical MOs obtained by the usual DFT calculation. For excited states, the Dyson orbitals provide a chemically sensible picture of the electronic transitions, thus bridging the chasm between orbital theory and response computations.

Electron binding energies and Dyson orbitals of OnH2n+1 +,0,- clusters: Double Rydberg anions, Rydberg radicals, and micro-solvated hydronium cations

Ab initio electron propagator methods are employed to predict the vertical electron attachment energies (VEAEs) of OH₃ ⁺(H₂O)_n clusters. The VEAEs decrease with increasing n, and the corresponding Dyson orbitals are diffused over exterior, non-hydrogen bonded protons. Clusters formed from OH₃ ^- double Rydberg anions (DRAs) and stabilized by hydrogen bonding or electrostatic interactions between ions and polar molecules are studied through calculations on OH₃ ^-(H₂O)_n complexes and are compared with more stable H^-(H₂O)_n+1 isomers. Remarkable changes in the geometry of the anionic hydronium-water clusters with respect to their cationic counterparts occur. Rydberg electrons in the uncharged and anionic clusters are held near the exterior protons of the water network. For all values of n, the anion-water complex H^-(H₂O)_n+1 is always the most stable, with large vertical electron detachment energies (VEDEs). OH₃ ^-(H₂O)_n DRA isomers have well separated VEDEs and may be visible in anion photoelectron spectra. Corresponding Dyson orbitals occupy regions beyond the peripheral O-H bonds and differ significantly from those obtained for the VEAEs of the cations.

Models of necessity

The way chemists represent chemical structures as two-dimensional sketches made up of atoms and bonds, simplifying the complex three-dimensional molecules comprising nuclei and electrons of the quantum mechanical description, is the everyday language of chemistry. This language uses models, particularly of bonding, that are not contained in the quantum mechanical description of chemical systems, but has been used to derive machine-readable formats for storing and manipulating chemical structures in digital computers. This language is fuzzy and varies from chemist to chemist but has been astonishingly successful and perhaps contributes with its fuzziness to the success of chemistry. It is this creative imagination of chemical structures that has been fundamental to the cognition of chemistry and has allowed thought experiments to take place. Within the everyday language, the model nature of these concepts is not always clear to practicing chemists, so that controversial discussions about the merits of alternative models often arise. However, the extensive use of artificial intelligence (AI) and machine learning (ML) in chemistry, with the aim of being able to make reliable predictions, will require that these models be extended to cover all relevant properties and characteristics of chemical systems. This, in turn, imposes conditions such as completeness, compactness, computational efficiency and non-redundancy on the extensions to the almost universal Lewis and VSEPR bonding models. Thus, AI and ML are likely to be important in rationalizing, extending and standardizing chemical bonding models. This will not affect the everyday language of chemistry but may help to understand the unique basis of chemical language.

Asymmetric Density Fitting with Modified Cholesky Decomposition Applied to Second-Order Electron Propagator

Computation of molecular orbital electron repulsion integrals (MO-ERIs) as a transformation from atomic orbital ERIs (AO-ERIs) is the bottleneck of second-order electron propagator calculations when a single orbital is studied. In this contribution, asymmetric density fitting is combined with modified Cholesky decomposition to generate efficiently the required MO-ERIs. The key point of the presented algorithms is to keep track of integrals through partial contractions performed on three-center AO-ERIs; these contractions are stored in RAM instead of the AO-ERIs. Two implementations are provided, an in-core, which reduces the arithmetic and memory scaling factors as compared to the four-center AO-ERIs contraction method, and a semidirect, which overcomes memory limitations by evaluating antisymmetrized MO-ERIs in batches. On the numerical side, the proposed approach is fast and stable. The bad effects due to ill conditioning, namely, several negative and close to zero eigenvalues due to machine round off errors of the matrix associated with the density fitting process, are effectively controlled by means of a modified Cholesky factorization that avoids the matrix inversion needed to perform the asymmetrical density fitting implementation. The numerical experience presented shows that the in-core implementation is highly competitive to perform calculations on medium and large basis sets, while the semidirect implementation has small variations in time by changes in the available memory. The general applicability is illustrated on a set of selected relatively large-size molecules.

Dyson orbitals within the fc-CVS-EOM-CCSD framework: theory and application to X-ray photoelectron spectroscopy of ground and excited states

Ionization energies and Dyson orbitals within frozen-core core–valence separated equation-of-motion coupled cluster singles and doubles (fc-CVS-EOM-CCSD) enable efficient and reliable calculations of standard XPS and of UV-pump/XPS probe spectra.

Dyson Orbitals and Double Rydberg Anions: Methylated, Annulated, and Paramagnetic

A double Rydberg anion (DRA) consists of a saturated, closed-shell, molecular cation and two electrons that occupy diffuse orbitals. Techniques of ab initio electron propagator theory (EPT) predict the existence and spectra of three new classes of DRAs. The first, with the formula NH_4-n(CH₃)_n^-, has vertical electron detachment energies (VEDEs) that vary between 0.24 and 0.39 eV and corresponding Dyson orbitals that accumulate near the periphery of N-H bonds. An internal hydrogen bond that forms a ring with five members occurs in the second class. In paramagnetic DRA isomers, electrons are assigned to two, diffuse, triplet-coupled spin-orbitals that localize outside the N-H bonds of a cationic, tetrahedral center or outside bonds on a nearby amide or methyl group. Effects of delocalization, dispersion, and radial correlation between diffuse electrons on VEDEs are described in terms of Dyson orbitals and their pole strengths. These concepts of EPT connect ground-state and spectral properties to each other and provide a rigorous, systematic, and insightful approach to predicting and characterizing novel patterns of chemical bonding and molecular electronic structure.