Prediction of photodynamics of 200 nm excited cyclobutanone with linear response electronic structure and ab initio multiple spawning.

Simulations of photochemical reaction dynamics have been a challenge to the theoretical chemistry community for some time. In an effort to determine the predictive character of current approaches, we predict the results of an upcoming ultrafast diffraction experiment on the photodynamics of cyclobutanone after excitation to the lowest lying Rydberg state (S2). A picosecond of nonadiabatic dynamics is described with ab initio multiple spawning. We use both time dependent density functional theory (TDDFT) and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) theory for the underlying electronic structure theory. We find that the lifetime of the S2 state is more than a picosecond (with both TDDFT and EOM-CCSD). The predicted ultrafast electron diffraction spectrum exhibits numerous structural features, but weak time dependence over the course of the simulations.

Ultrafast X-ray Spectroscopy of Intersystem Crossing in Hexafluoroacetylacetone: Chromophore Photophysics and Spectral Changes in the Face of Electron-Withdrawing Groups.

Intersystem crossings between singlet and triplet states represent a crucial relaxation pathway in photochemical processes. Herein, we probe the intersystem crossing in hexafluoro-acetylacetone with ultrafast X-ray transient absorption spectroscopy at the carbon K-edge. We observe the excited state dynamics following excitation with 266 nm UV light to the 1ππ* (S2) state with element and site-specificity using a broadband soft X-ray pulse produced by high harmonic generation. These results are compared to X-ray spectra computed from orbital optimized density functional theory methods. It is found that the electron-withdrawing fluorine atoms decongest the X-ray absorption spectrum by enhancing separation between features originating from different carbon atoms. This facilitates the elucidation of structural and electronic dynamics at the chromophore. The evolution of the core-to-valence resonances at the carbon K-edge reveals an ultrafast population transfer between the 1nπ* (S1) and 3ππ* (T1) states on a 1.6 ± 0.4 ps time scale, which is similar to the 1.5 ps time scale earlier observed for acetylacetone [ J. Am. Chem. Soc. 2017, 139, 16576-16583, DOI: 10.1021/jacs.7b07532]. It therefore appears that terminal fluorination has little influence on the intersystem crossing rate of the acetylacetone chromophore. In addition, the significant role of hydrogen-bond opened and twisted rotational isomers is elucidated in the excited state dynamics by comparison of the experimental transient X-ray spectra with theory.

Probing C-I bond fission in the UV photochemistry of 2-iodothiophene with core-to-valence transient absorption spectroscopy.

The UV photochemistry of small heteroaromatic molecules serves as a testbed for understanding fundamental photo-induced chemical transformations in moderately complex compounds, including isomerization, ring-opening, and molecular dissociation. Here, a combined experimental-theoretical study of 268 nm UV light-induced dynamics in 2-iodothiophene (C4H3IS) is performed. The dynamics are experimentally monitored with a femtosecond extreme ultraviolet (XUV) probe that measures iodine N-edge 4d core-to-valence transitions. Experiments are complemented by density functional theory calculations of both the pump-pulse induced valence excitations and the XUV probe-induced core-to-valence transitions. Possible intramolecular relaxation dynamics are investigated by ab initio molecular dynamics simulations. Gradual absorption changes up to ∼0.5 to 1 ps after excitation are observed for both the parent molecular species and emerging iodine fragments, with the latter appearing with a characteristic rise time of 160 ± 30 fs. Comparison of spectral intensities and energies with the calculations identifies an iodine dissociation pathway initiated by a predominant π → π* excitation. In contrast, initial excitation to a nearby n⟂ → σ* state appears unlikely based on a significantly smaller oscillator strength and the absence of any corresponding XUV absorption signatures. Excitation to the π → π* state is followed by contraction of the C-I bond, enabling a nonadiabatic transition to a dissociative π→σC-I* state. For the subsequent fragmentation, a relatively narrow bond-length region along the C-I stretch coordinate between 230 and 280 pm is identified, where the transition between the parent molecule and the thienyl radical + iodine atom products becomes prominent in the XUV spectrum due to rapid localization of two singly occupied molecular orbitals on the two fragments.

When Is a Bond Broken? The Polarizability Perspective.

The question of when a chemical bond can be said to be broken is of fundamental chemical interest but has not been widely studied. Herein we propose that the maxima of static polarizability along bond dissociation coordinates naturally define cutoff points for bond rupture, as they represent the onset of localization of shared electron density into constituent fragments. Examples of computed polarizability maxima over the course of bond cleavage in main-group and transition metal compounds are provided, across covalent, dative and charge-shift bonds. The behavior along reaction paths is also considered. Overall, the static polarizability is found to be a sensitive reporter of electronic structure reorganization associated with bond stretching, and thus can serve as a metric for describing bond cleavage (or diagnose the absence of a chemical bond).

Computing x-ray absorption spectra from linear-response particles atop optimized holes.

State specific orbital optimized density functional theory (OO-DFT) methods, such as restricted open-shell Kohn-Sham (ROKS), can attain semiquantitative accuracy for predicting x-ray absorption spectra of closed-shell molecules. OO-DFT methods, however, require that each state be individually optimized. In this Communication, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, which is capable of giving reasonable accuracy without requiring state-specific optimization. This is achieved by fully optimizing the core-hole through the core-ionized state, followed by the use of electron-addition configuration interaction singles to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method and can attain ∼0.6 eV rms error for the K-edges of C-F through the use of local functionals, such as PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment and can thus help reduce the computational cost of obtaining OO-DFT quality spectra. ROKS(STEX), therefore, appears to be a useful technique for the efficient prediction of x-ray absorption spectra.

Electron-Nuclear Dynamics Accompanying Proton-Coupled Electron Transfer.

Although photoinduced proton-coupled electron transfer (PCET) plays an essential role in photosynthesis, a full understanding of the mechanism is still lacking due to the complex nonequilibrium dynamics arising from the strongly coupled electronic and nuclear degrees of freedom. Here we report the photoinduced PCET dynamics of a biomimetic model system investigated by means of transient IR and two-dimensional electronic-vibrational (2DEV) spectroscopies, IR spectroelectrochemistry (IRSEC), and calculations utilizing long-range-corrected hybrid density functionals. This collective experimental and theoretical effort provides a nuanced picture of the complicated dynamics and synergistic motions involved in photoinduced PCET. In particular, the evolution of the 2DEV line shape, which is highly sensitive to the mixing of vibronic states, is interpreted by accurate computational modeling of the charge separated state and is shown to represent a gradual change in electron density distribution associated with a dihedral twist that occurs on a 120 fs time scale.

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems.

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

Revealing the nature of electron correlation in transition metal complexes with symmetry breaking and chemical intuition.

In this work, we provide a nuanced view of electron correlation in the context of transition metal complexes, reconciling computational characterization via spin and spatial symmetry breaking in single-reference methods with qualitative concepts from ligand-field and molecular orbital theories. These insights provide the tools to reliably diagnose the multi-reference character, and our analysis reveals that while strong (i.e., static) correlation can be found in linear molecules (e.g., diatomics) and weakly bound and antiferromagnetically coupled (monometal-noninnocent ligand or multi-metal) complexes, it is rarely found in the ground-states of mono-transition-metal complexes. This leads to a picture of static correlation that is no more complex for transition metals than it is, e.g., for organic biradicaloids. In contrast, the ability of organometallic species to form more complex interactions, involving both ligand-to-metal σ-donation and metal-to-ligand π-backdonation, places a larger burden on a theory's treatment of dynamic correlation. We hypothesize that chemical bonds in which inter-electron pair correlation is non-negligible cannot be adequately described by theories using MP2 correlation energies and indeed find large errors vs experiment for carbonyl-dissociation energies from double-hybrid density functionals. A theory's description of dynamic correlation (and to a less important extent, delocalization error), which affects relative spin-state energetics and thus spin symmetry breaking, is found to govern the efficacy of its use to diagnose static correlation.

Exploring spin symmetry-breaking effects for static field ionization of atoms: Is there an analog to the Coulson-Fischer point in bond dissociation?

Löwdin's symmetry dilemma is an ubiquitous issue in approximate quantum chemistry. In the context of Hartree-Fock (HF) theory, the use of Slater determinants with some imposed constraints to preserve symmetries of the exact problem may lead to physically unreasonable potential energy surfaces. On the other hand, lifting these constraints leads to the so-called broken symmetry solutions that usually provide better energetics, at the cost of losing information about good quantum numbers that describe the state of the system. This behavior has previously been extensively studied in the context of bond dissociation. This paper studies the behavior of different classes of HF spin polarized solutions (restricted, unrestricted, and generalized) in the context of ionization by strong static electric fields. We find that, for simple two electron systems, unrestricted Hartree-Fock (UHF) is able to provide a qualitatively good description of states involved during the ionization process (neutral, singly ionized, and doubly ionized states), whereas RHF fails to describe the singly ionized state. For more complex systems, even though UHF is able to capture some of the expected characteristics of the ionized states, it is constrained to a single Ms (diabatic) manifold in the energy surface as a function of field intensity. In this case, a better qualitative picture can be painted by using generalized Hartree-Fock as it is able to explore different spin manifolds and follow the lowest solution due to lack of collinearity constraints on the spin quantization axis.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package.

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Generalized single excitation configuration interaction: an investigation into the impact of the inclusion of non-orthogonality on the calculation of core-excited states.

In this paper, we investigate different non-orthogonal generalizations of the configuration interaction with single substitutions (CIS) method for the calculation of core-excited states. Fully non-orthogonal CIS (NOCIS) has been described previously for species with singlet and doublet ground states, and this paper reports the extension to molecules in their triplet ground state. In addition to NOCIS, we present a novel method, one-center NOCIS (1C-NOCIS), for open-shell molecules which is intermediate between NOCIS and the computationally less demanding static exchange approximation (STEX). We explore this hierarchy of spin-pure methods for core excitations of molecules with singlet, doublet, and triplet ground states. We conclude that, while NOCIS provides the best results and preserves the spatial symmetry of the wavefunction, 1C-NOCIS retains much of the accuracy of NOCIS at a dramatically reduced cost. For molecules with closed-shell ground states, STEX and 1C-NOCIS are identical.

Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations.

State-specific orbital optimized approaches are more accurate at predicting core-level spectra than traditional linear-response protocols, but their utility had been restricted due to the risk of "variational collapse" down to the ground state. We employ the recently developed square gradient minimization [D. Hait and M. Head-Gordon, J. Chem. Theory Comput. 16, 1699 (2020)] algorithm to reliably avoid variational collapse and study the effectiveness of orbital optimized density functional theory (DFT) at predicting second period element 1s core-level spectra of open-shell systems. Several density functionals (including SCAN, B3LYP, and ωB97X-D3) are found to predict excitation energies from the core to singly occupied levels with high accuracy (≤0.3 eV RMS error) against available experimental data. Higher excited states are, however, more challenging by virtue of being intrinsically multiconfigurational. We thus present a configuration interaction inspired route to self-consistently recouple single determinant mixed configurations obtained from DFT, in order to obtain approximate doublet states. This recoupling scheme is used to predict the C K-edge spectra of the allyl radical, the O K-edge spectra of CO+, and the N K-edge of NO2 with high accuracy relative to experiment, indicating substantial promise in using this approach for the computation of core-level spectra for doublet species [vs more traditional time dependent DFT, equation of motion coupled cluster singles and doubles (EOM-CCSD), or using unrecoupled mixed configurations]. We also present general guidelines for computing core-excited states from orbital optimized DFT.

Beyond the Coulson-Fischer point: characterizing single excitation CI and TDDFT for excited states in single bond dissociations.

Linear response time dependent density functional theory (TDDFT), which builds upon configuration interaction singles (CIS) and TD-Hartree-Fock (TDHF), is the most widely used class of excited state quantum chemistry methods and is often employed to study photochemical processes. This paper studies the behavior of the resulting excited state potential energy surfaces beyond the Coulson-Fischer (CF) point in single bond dissociations, when the optimal reference determinant is spin-polarized. Many excited states exhibit sharp kinks at the CF point, and connect to different dissociation limits via a zone of unphysical concave curvature. In particular, the unrestricted MS = 0 lowest triplet T1 state changes character, and does not dissociate into ground state fragments. The unrestricted MS = ±1 T1 CIS states better approximate the physical dissociation limit, but their degeneracy is broken beyond the CF point for most single bond dissociations. On the other hand, the MS = ±1 T1 TDHF states reach the asymptote too soon, by merging with the ground state from the CF point onwards. Use of local exchange-correlation functionals causes MS = ±1 T1 TDDFT states to resemble their unphysical MS = 0 counterpart. The 2 orbital, 2-electron model system of minimal basis H2 is analytically treated to understand the origin of these issues, revealing that the lack of double excitations is at the root of these remarkable observations. The behavior of excited state surfaces is also numerically examined for species like H2, NH3, C2H6 and LiH in extended basis sets.

Well-behaved versus ill-behaved density functionals for single bond dissociation: Separating success from disaster functional by functional for stretched H2.

Unrestricted density functional theory (DFT) methods are typically expected to describe the homolytic dissociation of nonpolar single bonds in neutral species with qualitative accuracy, due to the lack of significant delocalization error. We however find that many widely used density functional approximations fail to describe features along the dissociation curve of the simple H2 molecule. This is not a universal failure of DFT in the sense that many classic functionals like PBE and B3LYP give very reasonable results, as do some more modern methods like MS2. However, some other widely used functionals like B97-D (empirically fitted) and TPSS (non-empirically constrained) predict qualitatively wrong static polarizabilities, force constants, and some even introduce an artificial barrier against association of independent H atoms to form H2. The polarizability and force constant prediction failures appear to stem from incomplete spin localization into individual H atoms beyond the Coulson-Fischer point, resulting in "fractionally bonded" species where the ionic contributions to the Slater determinant are not completely eliminated, unlike the case of unrestricted Hartree-Fock. These errors therefore appear to be a consequence of poor self-consistent density prediction by the problematic functional. The same reasons could potentially lead to spurious barriers toward H atom association, indirectly also leading to incorrect forces. These unphysicalities suggest that the use of problematic functionals is probably unwise in ab initio dynamics calculations, especially if strong electrostatic interactions are possible.

Chemoenzymatic Platform for Synthesis of Chiral Organofluorines Based on Type†II Aldolases.

Aldolases are C-C bond forming enzymes that have become prominent tools for sustainable synthesis of complex synthons. However, enzymatic methods of fluorine incorporation into such compounds are lacking due to the rarity of fluorine in nature. Recently, the use of fluoropyruvate as a non-native aldolase substrate has arisen as a solution. Here, we report that the type II HpcH aldolases efficiently catalyze fluoropyruvate addition to diverse aldehydes, with exclusive (3S)-selectivity at fluorine that is rationalized by DFT calculations on a mechanistic model. We also measure the kinetic parameters of aldol addition and demonstrate engineering of the hydroxyl group stereoselectivity. Our aldolase collection is then employed in the chemoenzymatic synthesis of novel fluoroacids and ester derivatives in high stereopurity (d.r. 80-98 %). The compounds made available by this method serve as precursors to fluorinated analogs of sugars, amino acids, and other valuable chiral building blocks.

How accurate are static polarizability predictions from density functional theory? An assessment over 132 species at equilibrium geometry.

Static polarizabilities are the first response of the electron density to electric fields, and are therefore important for predicting intermolecular and molecule-field interactions. They also offer a global measure of the accuracy of the treatment of excited states by density functionals in a formally exact manner. We have developed a database of benchmark static polarizabilities for 132 small species at equilibrium geometry, using coupled cluster theory through triple excitations (extrapolated to the complete basis set limit), for the purpose of developing and assessing density functionals. The performance of 60 popular and recent functionals are also assessed, which indicates that double hybrid functionals perform the best, having RMS relative errors in the range of 2.5-3.8%. Many hybrid functionals also give quite reasonable estimates with 4-5% RMS relative error. A few meta-GGAs like mBEEF and MVS yield performance comparable to hybrids, indicating potential for improved excited state predictions relative to typical local functionals. Some recent functionals however are found to be prone to catastrophic failure (possibly as a consequence of overparameterization), indicating a need for caution in applying these.

Communication: xDH double hybrid functionals can be qualitatively incorrect for non-equilibrium geometries: Dipole moment inversion and barriers to radical-radical association using XYG3 and XYGJ-OS.

Double hybrid (DH) density functionals are amongst the most accurate density functional approximations developed so far, largely due to the incorporation of correlation effects from unoccupied orbitals via second order perturbation theory (PT2). The xDH family of DH functionals calculate energy directly from orbitals optimized by a lower level approach like B3LYP, without self-consistent optimization. XYG3 and XYGJ-OS are two widely used xDH functionals that are known to be quite accurate at equilibrium geometries. Here, we show that the XYG3 and XYGJ-OS functionals can be ill behaved for stretched bonds well beyond the Coulson-Fischer point, predicting unphysical dipole moments and humps in potential energy curves for some simple systems like the hydrogen fluoride molecule. Numerical experiments and analysis show that these failures are not due to PT2. Instead, a large mismatch at stretched bond-lengths between the reference B3LYP orbitals and the optimized orbitals associated with the non-PT2 part of XYG3 leads to an unphysically large non-Hellman-Feynman contribution to first order properties like forces and electron densities.

A hybrid memory kernel approach for condensed phase non-adiabatic dynamics.

The spin-boson model is a simplified Hamiltonian often used to study non-adiabatic dynamics in large condensed phase systems, even though it has not been solved in a fully analytic fashion. Herein, we present an exact analytic expression for the dynamics of the spin-boson model in the infinitely slow-bath limit and generalize it to approximate dynamics for faster baths. We achieve the latter by developing a hybrid approach that combines the exact slow-bath result with the popular non-interacting blip approximation (NIBA) method to generate a memory kernel that is formally exact to second-order in the diabatic coupling but also contains higher-order contributions approximated from the second-order term alone. This kernel has the same computational complexity as the NIBA, but is found to yield dramatically superior dynamics in regimes where the NIBA breaks down-such as systems with large diabatic coupling or energy bias. This indicates that this hybrid approach could be used to cheaply incorporate higher-order effects into second-order methods and could potentially be generalized to develop alternate kernel resummation schemes.

Condensed phase electron transfer beyond the Condon approximation.

Condensed phase electron transfer problems are often simplified by making the Condon approximation: the approximation that the coupling connecting two charge-transfer diabatic states is a constant. Unfortunately, the Condon approximation does not predict the existence of conical intersections, which are ubiquitous in both gas-phase and condensed-phase photochemical dynamics. In this paper, we develop a formalism to treat condensed-phase dynamics beyond the Condon approximation. We show that even for an extremely simple test system, hexaaquairon(ii)/hexaaquairon(iii) self-exchange in water, the electronic coupling is expected to fluctuate rapidly and non-Condon effects must be considered to obtain quantitatively accurate ultrafast nonequilibrium dynamics. As diabatic couplings are expected to fluctuate substantially in many condensed-phase electron transfer systems, non-Condon effects may be essential to quantitatively capture accurate short-time dynamics.

Predicting the X-ray Absorption Spectrum of Ozone with Single Configuration State Functions.

X-ray absorption spectra (XAS) of biradicaloid species are often thought to represent a challenge to theoretical methods. This has led to the testing of recently developed multireference techniques on the XAS of ozone, but reproduction of the experimental spectral profile has proven difficult. We utilize a minimal model consisting of a single configuration state function (CSF) per excited state to model core-level excitations of ozone, with the orbitals of each CSF optimized using the restricted open-shell Kohn-Sham (ROKS) method. This protocol leads to semiquantitative agreement with experimental XAS. In fact, we find that low-lying core-hole excited states in biradicaloids can be approximated with individual CSFs, despite the presence of multireference character in the ground state. We also report that the 1s → π* and 1s → σ* transitions have quite distinct widths for O3. This reveals the importance of sampling over a representative range of geometries from the vibrational ground state for properly assessing the accuracy of electronic structure methods against experiments instead of the popular procedure of uniformly broadening stick spectra at the equilibrium geometry.