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.

Generalized perturbative singles corrections to the random phase approximation method: Impact on noncovalent interaction energies of closed- and open-shell dimers

The post-Kohn-Sham (KS) random phase approximation (RPA) method may provide a poor description of interaction energies of weakly bonded molecules due to inherent density errors in approximate KS functionals. To overcome these errors, we develop a generalized formalism to incorporate perturbative singles (pS) corrections to the RPA method using orbital rotations as a perturbation parameter. The pS schemes differ in the choice of orbital-rotation gradient and Hessian. We propose a pS scheme termed RPA singles (RPAS)[Hartree-Fock (HF)] that uses the RPA orbital-rotation gradient and time-dependent HF Hessian. This correction reduces the errors in noncovalent interaction energies of closed- and open-shell dimers. For the open-shell dimers, the RPAS(HF) method leads to a consistent error reduction by 50% or more compared to the RPA method for the cases of hydrogen-bonding, metal-solvent, carbene-solvent, and dispersion interactions. We also find that the pS corrections are more important in error reduction compared to higher-order exchange corrections to the RPA method. Overall, for open shells, the RPAS(HF)-corrected RPA method provides chemical accuracy for noncovalent interactions and is more reliable than other perturbative schemes and dispersion-corrected density functional approximations, highlighting its importance as a reliable beyond-RPA correction.

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.

Excited state dynamics of a spontaneously generated TTF radical cation inside a water-soluble nanocage

Supramolecular cavities have been traditionally used to stabilize reactive redox intermediates. Recently with the success of multiple new photoredox catalytic strategies that use supramolecular cages, there is a growing demand for photogeneration strategies of diverse reactive intermediates inside confined spaces, which will drive enzyme-like catalysis in real time. Here we report the excited state dynamics of a redox-active TTF radical cation and its corresponding dimethyl-derivative DiMeTTF inside a confined supramolecular cavity. We prepare the radical cation by spontaneous oxidation of neutral TTF upon incarceration inside a water-soluble nanocage Pd₆L₄¹²⁺, and characterize it with a combination of resonance Raman and electron paramagnetic resonance spectroscopy. Using broadband transient absorption spectroscopy, we demonstrate that the confined native TTF radical cation and its dimethyl derivative upon photoexcitation rapidly de-excite to form the hot ground state, thereby inhibiting further oxidation to a TTF⁺² dication. We discuss our results in the context of excited state crossings of the radical cation potentials as well as modifying the cage energetics to generate a stable dication. Our work has important implications for the usage of such radical cations for photoactivated catalysis.

Exploring the Solvation of Acetic Acid in Water Using Liquid Jet X-ray Photoelectron Spectroscopy and Core Level Electron Binding Energy Calculations

Liquid jet X-ray photoelectron spectroscopy was used to investigate changes in the local electronic structure of acetic acid in the bulk of aqueous solutions induced by solvation effects. These effects manifest themselves as shifts in the difference in the carbon 1s binding energy (ΔBE) between the methyl and carboxyl carbons of acetic acid. Furthermore, molecular dynamics simulations, coupled with correlated electronic structure calculations of the first solvation sphere, provide insight into the number of water molecules directly interacting with the carboxyl group that are required to match the ΔBE from the photoelectron spectroscopy experiments. This comparison shows that a single water molecule in the first solvation shell describes the photoelectron ΔBE of acetic acid while at least 20 water molecules are required for the conjugate base, acetate, in aqueous solutions.

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

The suitability of one-particle energies from the generalized Kohn-Sham semicanonical projected random phase approximation (GKS-spRPA) method for electron affinities of molecules is investigated. It is shown that the GKS-spPRA effective potential includes exact exchange and polarizability-dependent correlation terms that are necessary for the correct description of anionic systems. An O(N4) implementation that enables fast computation of electron affinities is presented. For model systems, I show that the GKS-spRPA approach is applicable for valence and nonvalence type anions with a maximum error of 0.13 eV for valence anionic states and 0.03 eV for nonvalence anionic states compared to equation of motion coupled cluster methods. For a series of perhalobenzene molecules, C₆X₆ (X= F, Cl, Br, and I), GKS-spRPA predicts that the ground-state character changes from a nonvalence-σ* type in C₆F₆^- to valence-π* in C₆Cl₆^- and valence-σ* in C₆Br₆^- and C₆I₆^-. Experimental implications of these findings are discussed.

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

Effective one-particle energies from generalized Kohn-Sham random phase approximation: A direct approach for computing and analyzing core ionization energies

Generalized-Kohn-Sham (GKS) orbital energies obtained self-consistently from the random phase approximation energy functional with a semicanonical projection (spRPA) were recently shown to rival the accuracy of GW quasiparticle energies for valence ionization potentials. Here, we extend the scope of GKS-spRPA correlated one-particle energies from frontier-orbital ionization to core orbital ionization energies, which are notoriously difficult for GW and other response methods due to strong orbital relaxation effects. For a benchmark consisting of 23 1s core electron binding energies (CEBEs) of second-row elements, chemical shifts estimated from GKS-spRPA one-particle energies yield mean absolute deviations from experiment of 0.2 eV, which are significantly more accurate than the standard GW and comparable to Δ self-consistent field theory without semiempirical adjustment of the energy functional. For small ammonia clusters and cytosine tautomers, GKS-spRPA based chemical shifts capture subtle variations in covalent and noncovalent bonding environments; GKS-spRPA 1s CEBEs for these systems agree with equation-of-motion coupled cluster singles and doubles and ADC(4) results within 0.2-0.3 eV. Two perturbative approximations to GKS-spRPA orbital energies, which reduce the scaling from O(N⁶) to O(N⁵) and O(N⁴), are introduced and tested. We illustrate the application of GKS-spRPA orbital energies to larger systems by using oxygen 1s CEBEs to probe solvation and packing effects in condensed phases of water. GKS-spRPA predicts a lowering of the oxygen 1s CEBE of approximately 1.6-1.7 eV in solid and liquid phases, consistent with liquid-jet X-ray photoelectron spectroscopy and gas phase cluster experiments. The results are rationalized by partitioning GKS-spRPA electron binding energies into static, relaxation, and correlation parts.

Understanding the role of intermolecular interactions between lissoclimides and the eukaryotic ribosome

Natural products that target the eukaryotic ribosome are promising therapeutics to treat a variety of cancers. It is therefore essential to determine their molecular mechanism of action to fully understand their mode of interaction with the target and to inform the development of new synthetic compounds with improved potency and reduced cytotoxicity. Toward this goal, we have previously established a short synthesis pathway that grants access to multiple congeners of the lissoclimide family. Here we present the X-ray co-crystal structure at 3.1 Å resolution of C45, a potent congener with two A-ring chlorine-bearing stereogenic centers with ‘unnatural’ configurations, with the yeast 80S ribosome, intermolecular interaction energies of the C45/ribosome complex, and single-molecule FRET data quantifying the impact of C45 on both human and yeast ribosomes. Together, these data provide new insights into the role of unusual non-covalent halogen bonding interactions involved in the binding of this synthetic compound to the 80S ribosome.

Comparisons of lanthanide/actinide +2 ions in a tris(aryloxide)arene coordination environment

Nd, like U, prefers a f⁴ configuration with the tris(aryloxide)arene ligand rather than the 4f³5d¹ configuration found in tris(cyclopentadienyl) complexes.

A new series of Ln³⁺ and Ln²⁺ complexes has been synthesized using the tris(aryloxide)arene ligand system, ((^Ad,MeArO)₃mes)^3–, recently used to isolate a complex of U²⁺. The triphenol precursor, (^Ad,MeArOH)₃mes, reacts with the Ln³⁺ amides, Ln(NR₂)₃ (R = SiMe₃), to form a series of [((^Ad,MeArO)₃mes)Ln] complexes, 1-Ln. Crystallographic characterization was achieved for Ln = Nd, Gd, Dy, and Er. The complexes 1-Ln can be reduced with potassium graphite in the presence of 2.2.2-cryptand (crypt) to form highly absorbing solutions with properties consistent with Ln²⁺ complexes, [K(crypt)][((^Ad,MeArO)₃mes)Ln], 2-Ln. The synthesis of the Nd²⁺ complex [K(crypt)][((^Ad,MeArO)₃mes)Nd], 2-Nd, was unambiguously confirmed by X-ray crystallography. In the case of the other lanthanides, crystals were found to contain mixtures of 2-Ln co-crystallized with either a Ln³⁺ hydride complex, [K(crypt)][((^Ad,MeArO)₃mes)LnH], 3-Ln, for Ln = Gd, Dy, and Er, or a hydroxide complex, [K(crypt)][((^Ad,MeArO)₃mes)Ln(OH)], 4-Ln, for Ln = Dy. A Dy²⁺ complex with 18-crown-6 as the potassium chelator, [K(18-crown-6)(THF)₂][((^Ad,MeArO)₃mes)Dy], 5-Dy, was isolated as a co-crystallized mixture with the Dy³⁺ hydride complex, [K(18-crown-6)(THF)₂][((^Ad,MeArO)₃mes)DyH], 6-Dy. Structural comparisons of 1-Ln and 2-Ln are presented with respect to their uranium analogs and correlated with density functional theory calculations on their electronic structures.

Synthesis facilitates an understanding of the structural basis for translation inhibition by the lissoclimides

The lissoclimides are unusual succinimide-containing labdane diterpenoids that were reported to be potent cytotoxins. Our short semisynthesis and analogue-oriented synthesis approaches provide a series of lissoclimide natural products and analogues that expand the structure-activity relationships (SARs) in this family. The semisynthesis approach yielded significant quantities of chlorolissoclimide (CL) to permit an evaluation against the National Cancer Institute's 60-cell line panel and allowed us to obtain an X-ray co-crystal structure of the synthetic secondary metabolite with the eukaryotic 80S ribosome. Although it shares a binding site with other imide-based natural product translation inhibitors, CL engages in a particularly interesting and novel face-on halogen-π interaction between the ligand's alkyl chloride and a guanine residue. Our analogue-oriented synthesis provides many more lissoclimide compounds, which were tested against aggressive human cancer cell lines and for protein synthesis inhibitory activity. Finally, computational modelling was used to explain the SARs of certain key compounds and set the stage for the structure-guided design of better translation inhibitors.

Nonvalence Correlation-Bound Anion States of Polycyclic Aromatic Hydrocarbons

In this work, we characterize the nonvalence correlation-bound anion states of several polycyclic aromatic hydrocarbon (PAH) molecules. Unlike the analogous image potential states of graphene that localize the charge density of the excess electron above and below the plane of the sheet, we find that for PAHs, much of the charge distribution of the excess electron is localized around the periphery of the molecule. This is a consequence of the electrostatic interaction of the electron with the polar CH groups. By replacing the H atoms by F atoms or the CH groups by N atoms, the charge density of the excess electron shifts from the periphery to above and below the plane of the ring systems.

Nonvalence correlation-bound anion states of spherical fullerenes

We present a one-electron model Hamiltonian for characterizing nonvalence correlation-bound anion states of fullerene molecules. These states are the finite system analogs of image potential states of metallic surfaces. The model potential accounts for both atomic and charge-flow polarization and is used to characterize the nonvalence correlation-bound anion states of the C60, (C60)2, C240, and C60@C240 fullerene systems. Although C60 is found to have a single (s-type) nonvalence correlation-bound anion state, the larger fullerenes are demonstrated to have multiple nonvalence correlation-bound anion states.

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.