On the Cholesky decomposition for electron propagator methods: General aspects and application on C60

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.

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.

How many bound valence states does the C₆₀⁻ anion have?

We report on unprecedentedly large coupled cluster calculations for the C60(-) anion, and on a heuristic model uncovering the valence states of C60(-) that allow the resolution of the headlined question. Our results convincingly demonstrate that C60(-) possesses as many as four bound valence states: (2)T1u, (2)T1g, (2)T2u and (2)Hg. Our findings reconcile previous controversies regarding the existence of the bound (2)T2u and (2)Hg states. For all bound states of C60(-) we present an analysis of the radial and angular distributions of the excess electron, which reveals some unique properties of the valence states. Some interesting features of the introduced model are analyzed and discussed.

Quantifying the relative molecular orbital alignment for molecular junctions with similar chemical linkage to electrodes

Estimating the relative alignment between the frontier molecular orbitals (MOs) that dominates the charge transport through single-molecule junctions represents a challenge for theory. This requires approaches beyond the widely employed framework provided by the density functional theory, wherein the Kohn-Sham 'orbitals' are treated as if they were real MOs, which is not the case. In this paper, we report results obtained by means of quantum chemical calculations, including the equation-of-motion coupled-cluster singles and doubles, which is the state-of-the-art of quantum chemistry for medium-size molecules like those considered here. These theoretical results are validated against data on the MO energy offset relative to the electrodes' Fermi energy extracted from experiments for junctions based on 4,4'-bipyridine and 1,4-dicyanobenzene.

A quantum chemical study from a molecular transport perspective: ionization and electron attachment energies for species often used to fabricate single-molecule junctions

The accurate determination of the lowest electron attachment (EA) and ionization (IP) energies for molecules embedded in molecular junctions is important for correctly estimating, for example, the magnitude of the currents (I) or the biases (V) where an I-V curve exhibits significant non-Ohmic behavior. Benchmark calculations for the lowest electron attachment and ionization energies of several typical molecules utilized to fabricate single-molecule junctions characterized by n-type conduction (4,4'-bipyridine, 1,4-dicyanobenzene and 4,4'-dicyano-1,1'-biphenyl) and p-type conduction (benzenedithiol, biphenyldithiol, hexanemonothiol and hexanedithiol) based on the EOM-CCSD (equation-of-motion coupled-cluster singles and doubles) state-of-the-art method of quantum chemistry are presented. They indicate significant differences from the results obtained within current approaches to molecular transport. The present study emphasizes that, in addition to a reliable quantum chemical method, basis sets much better than the ubiquitous double-zeta set employed for transport calculations are needed. The latter is a particularly critical issue for correctly determining EAs, which is impossible without including sufficient diffuse basis functions. The spatial distribution of the dominant molecular orbitals (MOs) is another important issue, on which the present study draws attention, because it sensitively affects the MO energy shifts Φ due to image charges formed in electrodes. The present results cannot substantiate the common assumption of a point-like MO midway between electrodes, which substantially affects the actual Φ-values.

Existence of a Correlation Bound s-Type Anion State of C60

It is established using high-level electronic structure calculations that C60 has an s-type correlation-bound anion state with an electron binding energy of about 120 meV. Examination of the "singly occupied" natural orbital of the anion reveals that about 9% of the charge density of the excess electron is localized inside, and about 91% is localized outside the C60 cage. Calculations were also carried out for the He@C60, Ne@C60, and H2O@C60 endohedral complexes. For each of these species, the s-type anion is predicted to be less strongly bound than for C60 itself.

Anions of xenon clusters bound by long-range electron correlations

In contrast with the single atom, atomic van der Waals clusters can form stable anions where the excess electron is bound due to long-range correlations with the electrons of the cluster. We report on extensive all-electron many-body ab initio studies on Xe clusters. Three-dimensional, planar, and linear structures of the clusters are investigated and compared. In particular, we find that the minimal number of Xe atoms in the cluster required to form a stable anion is 5 independently of the dimensionality of the cluster. We provide electron affinities for clusters made of 5, 6, and 7 atoms in all dimensions and find that the planar clusters form the most stable anions. The Dyson orbitals of the excess electrons are computed and analyzed.