Molecular Electron Affinities Using the Generalized Kohn-Sham Semicanonical Projected Random Phase Approximation

Taming Negative Ion Resonances Using Nonlocal Exchange-Correlation Functionals

The characterization of negative ion resonances poses a fundamental challenge to density functional methods due to the unbound nature of resonances. To overcome this challenge, we propose one-particle nonlocal exchange-correlation (xc) potentials combining the exact-exchange (EXX) and the random phase approximation (RPA) correlation potentials. The negative ion resonances are identified by perturbing the real Hermitian nonlocal xc potentials using complex absorbing local potentials. Our studies show that the nonlocal EXX+RPA potential significantly enhances the description of positions and widths of negative ion resonance states compared to potentials that exclude dynamic polarization in RPA or include only EXX. The use of low-scaling algorithms simplifies the computation of the RPA potential, thereby providing a practical solution for resonance-state characterization within the density functional framework. A theoretical framework and the underlying assumptions required for combining real Hermitian nonlocal xc potentials with complex local potentials are discussed.

Single-Pole Polarization Models: Rapid Evaluation of Electron Affinities of Solvated-Electron and Superatomic Molecular Anionic States

We propose a single-parameter effective one-particle potential, termed the single-pole exchange-correlation (1p-XC), to rapidly evaluate electron affinities (EAs) of nonvalence electronic states of molecular clusters and nanoassemblies. The model combines exact-exchange and the random phase approximation (RPA) correlation potential with a single-pole approximation to model the frequency-dependent polarization function. It captures long-range static and dynamic-frequency effects in the correlation potential, with mean absolute errors of 0.06 eV for EAs of hydrated- and ammoniated-electron clusters with EA values in the range 0.24-1.77 eV. The 1p-XC approximation enables EA estimation with a computational wall-time similar to that of hybrid functionals. The model also provides a compressed-basis, which significantly reduces the rank of higher-level parameter-free one-particle Hamiltonians and further simplifies the computation of EAs. The compressed-basis approach is used to model the hybridization of superatomic molecular states of (C60)2- and (C60)3-, thereby verifying previous model Hamiltonian studies.

Chemically accurate singlet-triplet gaps of organic chromophores and linear acenes by the random phase approximation and σ-functionals

Predicting the energy differences between different spin-states is challenging for many widely used ab initio electronic structure methods. We here assess the ability of the direct random phase approximation (dRPA), dRPA plus two different screened second-order exchange (SOX) corrections, and σ-functionals to predict adiabatic singlet-triplet gaps. With mean absolute deviations of below 0.1 eV to experimental reference values, independent of the Kohn-Sham starting point, dRPA and σ-functionals accurately predict singlet-triplet gaps of 18 organic chromophores. The addition of SOX corrections to dRPA considerably worsens agreement with experiment, adding to the mounting evidence that dRPA+SOX methods are not generally applicable beyond-RPA methods. Also for a series of linear acene chains with up to ten fused rings, dRPA, and σ-functionals are in excellent agreement with coupled-cluster single double triple reference data. In agreement with advanced multi-reference methods, dRPA@PBE and σ-functional@PBE predict a singlet ground state for all chain lengths, while dRPA@PBE0 and σ-functional@PBE0 predict a triplet ground state for longer acenes. Our work shows dRPA and σ-functionals to be reliable methods for calculating singlet-triplet gaps in aromatic molecules.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Accelerating Core-Level GW Calculations by Combining the Contour Deformation Approach with the Analytic Continuation of W

In recent years, the GW method has emerged as a reliable tool for computing core-level binding energies. The contour deformation (CD) technique has been established as an efficient, scalable, and numerically stable approach to compute the GW self-energy for deep core excitations. However, core-level GW calculations with CD face the challenge of higher scaling with respect to system size N compared to the conventional quartic scaling in valence-state algorithms. In this work, we present the CD-WAC method [CD with W analytic continuation (AC)], which reduces the scaling of CD applied to the inner shells from O(N5) to O(N4) by employing an AC of the screened Coulomb interaction W. Our proposed method retains the numerical accuracy of CD for the computationally challenging deep core case, yielding mean absolute errors <5 meV for well-established benchmark sets, such as CORE65, for single-shot GW calculations. More extensive testing for different GW flavors proves the reliability of the method. We have confirmed the theoretical scaling by performing scaling experiments on large acene chains and amorphous carbon clusters, achieving speedups of up to 10× for structures of only 116 atoms. This improvement in computational efficiency paves the way for more accurate and efficient core-level GW calculations on larger and more complex systems.

Quantifying the Impact of Halogenation on Intermolecular Interactions and Binding Modes of Aromatic Molecules

Halogenation of aromatic molecules is frequently used to modulate intermolecular interactions with ramifications for optoelectronic and mechanical properties. In this work, we accurately quantify and understand the nature of intermolecular interactions in perhalogenated benzene (PHB) clusters. Using benchmark binding energies from the fixed-node diffusion Monte Carlo (FN-DMC) method, we show that generalized Kohn-Sham semicanonical projected random phase approximation (GKS-spRPA) plus approximate exchange kernel (AKX) provides reliable interaction energies with mean absolute error (MAE) of 0.23 kcal/mol. Using the GKS-spRPA+AXK method, we quantify the interaction energies of several binding modes of PHB clusters ((C₆X₆)_n; X = F, Cl, Br, I; n = 2, 3). For a given binding mode, the interaction energies increase 3-4 times from X = F to X = I; the X-X binding modes have energies in the range of 2-4 kcal/mol, while the π-π binding mode has interaction energies in the range of 4-12 kcal/mol. SAPT-DFT-based energy decomposition analysis is then used to show that the equilibrium geometries are dictated primarily by the dispersion and exchange interactions. Finally, we test the accuracy of several dispersion-corrected density functional approximations and show that only the r2SCAN-D4 method has a low MAE and correct long-range behavior, which makes it suitable for large-scale simulations and for developing structure-function relationships of halogenated aromatic systems.

Modeling Nonresonant X-ray Emission of Second- and Third-Period Elements without Core-Hole Reference States and Empirical Parameters

Nonresonant X-ray emission (XE) energies and oscillator strengths are obtained using the effective potential of the generalized Kohn-Sham semi-canonical projected random phase approximation (GKS-spRPA) method. XE energies are estimated as a difference between the valence and core ionization eigenvalues, while the oscillator strengths are obtained within a frozen orbital approximation. This straightforward approach provides accurate XE energies without any need for core-hole reference states, empirical shifting parameters, or tuning of density functionals. To account for relativistic corrections to the core orbitals, we have formulated a scalar relativistic (sr) GKS-spRPA approach based on the spin-free X2C one-electron Hamiltonian. The sr-GKS-spRPA method provides highly reliable XE energies using uncontracted basis-sets on atoms where the core-hole is created prior to emission. For the largest basis-sets used in our study, using completely uncontracted polarized core-valence Dunning basis-sets, the mean absolute errors (MAEs) are within 0.7 eV compared to experimental reference values for a test-set consisting of 27 valence-to-core XE energies of molecules with second- and third-period elements. Considering a balance of accuracy and computational effort, we recommend the use of s-uncontracted def2-TZVP for second-period and all-uncontracted def2-TZVP for third-period elements. For this recommended basis-set, the MAE is 0.2 eV. The analytically continued sr-GKS-spRPA approach, with an O(N4) computational cost, enables efficient computation of XE spectra of molecules such as S₈ and C₆₀ with several core-hole states.

Sensing the ortho Positions in C6Cl6 and C6H4Cl2 from Cl2− Formation upon Molecular Reduction

The geometrical effect of chlorine atom positions in polyatomic molecules after capturing a low-energy electron is shown to be a prevalent mechanism yielding Cl₂⁻. In this work, we investigated hexachlorobenzene reduction in electron transfer experiments to determine the role of chlorine atom positions around the aromatic ring, and compared our results with those using ortho-, meta- and para-dichlorobenzene molecules. This was achieved by combining gas-phase experiments to determine the reaction threshold by means of mass spectrometry together with quantum chemical calculations. We also observed that Cl₂⁻ formation can only occur in 1,2-C₆H₄Cl₂, where the two closest C–Cl bonds are cleaved while the chlorine atoms are brought together within the ring framework due to excess energy dissipation. These results show that a strong coupling between electronic and C–Cl bending motion is responsible for a positional isomeric effect, where molecular recognition is a determining factor in chlorine anion formation.

Hexachlorobenzene-negative ion formation in electron attachment experiments

In this contribution, we report a comprehensive study on the formation of hexachlorobenzene (C₆Cl₆) negative ions probed by low-energy electron interactions from 0 up to 12 eV in a gas-phase crossed beam experiment. The anionic yields as a function of the electron energy reveal a rich fragmentation pattern of the dissociative electron attachment process, yet the most intense ion has been assigned to the non-dissociated parent anion that survives long enough within the detection time window. Other less intense fragment anions have been assigned as Cl^-, Cl₂^-, C₆Cl₄^-, and C₆Cl₅^-. The experimental results are accompanied by quantum chemical calculations at various levels of accuracy, providing an insight into the electronic structure, thermochemical thresholds, electron affinities and structures of neutral and anionic molecular species. The electron attachment process induces a considerable geometry change in the temporary-negative ion relative to the neutral molecule, where the most intense fragment anion assigned to Cl^- can be formed solely through a curve crossing involving a π*/σ* coupling. The yield of chlorine anions shows a signature of vibrational excitation reminiscent of a Jahn-Teller distortion.

Dissociative electron attachment to hexachlorobenzene

Gas phase molecules of hexachlorobenzene (C 6 Cl 6 ) were investigated by means of Dissociative Electron Attachment Spectroscopy (DEAS). Three channels of molecular negative ions decay have been recorded: Cl- and Cl 2 - as well as electron detachment (t a ~250 ms at 343 K). All of these decay channels exhibit temperature dependence. The electron affinity was estimated within the framework of the simple Arrhenius model (EA=1.6-1.9 eV), as well as quantum-chemically at various levels (EA a =0.9-1.0 eV). The reasons for the discrepancy between the obtained estimates of the EA value and the results of previous works (EA from 1 to 2 eV) are discussed.

Anionic states of C6Cl6 probed in electron transfer experiments

This is the first comprehensive investigation on the anionic species formed during collisions of fast neutral potassium (K) atoms with neutral hexachlorobenzene (C₆Cl₆) molecules in the laboratory frame range from 10 up to 100 eV. In such ion-pair formation experiments we also report a novel K⁺ energy loss spectrum obtained in the forward scattering giving evidence of the most accessible electronic states. The vertical electron affinity of (-3.76 ± 0.20) eV has been obtained and assigned to a purely repulsive transition from the C₆Cl₆ ground state to a state of the temporary negative ion yielding Cl^- formation. These experimental findings are also supported by state-of-the art theoretical calculations on the electronic structure of C₆Cl₆ in the presence of a potassium atom and are used for analysing the lowest unoccupied molecular orbitals participating in the collision process. From the time-of-flight mass spectra recorded in the wide collision energy range, more than 80% of the total anion yield is due to the undissociated parent anion C₆Cl₆^-, C₆Cl₅^- and Cl^- formation. Other fragment anions such as C₆Cl₄^-, C₃Cl₂^-, C₂Cl^- and Cl₂^- that undergo complex internal reactions with the temporary negative ion formed after electron transfer account for less than 20% of the total yield. The joint experimental and theoretical methodologies employed in these electron transfer studies provide the most comprehensive and unique assignments of the hexachlorobenzene anionic species and the role of C₆Cl₆ electronic states in collision induced dissociation to date.

Selfconsistent random phase approximation methods

This Perspective reviews recent efforts toward selfconsistent calculations of ground-state energies within the random phase approximation (RPA) in the (generalized) Kohn-Sham (KS) density functional theory context. Since the RPA correlation energy explicitly depends on the non-interacting KS potential, an additional condition to determine the energy as a functional of the density is necessary. This observation leads to the concept of functional selfconsistency (FSC), which requires that the KS density equals the interacting density defined as the functional derivative of the ground-state energy with respect to the external potential. While all existing selfconsistent RPA schemes violate FSC, the recent generalized KS semicanonical projected RPA (GKS-spRPA) method takes a step toward satisfying it. This leads to systematic improvements in densities, binding energy curves, reference state stability, and molecular properties compared to non-selfconsistent RPA as well as optimized effective potential RPA. GKS-spRPA orbital energies accurately approximate valence and core ionization potentials, and even electron affinities of non-valence bound anions. The computational cost and performance of GKS-spRPA are compared to those of related selfconsistent schemes, including GW and orbital optimization methods, and limitations are discussed. Large differences between KS and interacting densities observed in the absence of FSC and the well-rounded performance of GKS-spRPA suggest that the KS potential as a density functional should be defined via the FSC condition for explicitly potential-dependent density functionals.