Indolo[2,3-b]quinoxaline as a Low Reduction Potential and High Stability Anolyte Scaffold for Nonaqueous Redox Flow Batteries.

Redox flow batteries (RFBs) are a promising stationary energy storage technology for leveling power supply from intermittent renewable energy sources with demand. A central objective for the development of practical, scalable RFBs is to identify affordable and high-performance redox-active molecules as storage materials. Herein, we report the design, synthesis, and evaluation of a new organic scaffold, indolo[2,3-b]quinoxaline, for highly stable, low-reduction potential, and high-solubility anolytes for nonaqueous redox flow batteries (NARFBs). The mixture of 2- and 3-(tert-butyl)-6-(2-methoxyethyl)-6H-indolo[2,3-b]quinoxaline exhibits a low reduction potential (-2.01 V vs Fc/Fc+), high solubility (>2.7 M in acetonitrile), and remarkable stability (99.86% capacity retention over 49.5 h (202 cycles) of H-cell cycling). This anolyte was paired with N-(2-(2-methoxyethoxy)-ethyl)phenothiazine (MEEPT) to achieve a 2.3 V all-organic NARFB exhibiting 95.8% capacity retention over 75.1 h (120 cycles) of cycling.

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory.

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter (α) was set to recover the exact dissociation limit of H2 in a minimal basis set. Alternatively α can be viewed as a regularization parameter, where each value of α represents a valid variant of BW-s2, which we denote as BW-s2(α). In this work, we semi-empirically optimize α for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the ab initio parameterization of BW-s2 and MP2. We demonstrate that the optimal α parameter (α = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2(α) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2(α) is iterative O(N5), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond.

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within "chemical accuracy" (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob's ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 - N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

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.

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.

Multiple Stable Isoprene-Ozone Complexes Reveal Complex Entrance Channel Dynamics in the Isoprene + Ozone Reaction.

Accurately characterizing isoprene ozonolysis continues to challenge atmospheric chemists. The reaction is believed to be a spontaneous, concerted cycloaddition. However, little information is available about the entrance channel and isoprene-ozone complexes thought to define the long-range portion of the reaction coordinate. Our coupled cluster and auxiliary field quantum Monte Carlo calculations predict multiple stable isoprene-ozone van der Waals complexes for trans-isoprene in the gas phase with moderate association energies. These results indicate that long-range dynamics in the isoprene-ozone entrance channel can impact the overall reaction in the troposphere and provide the spectroscopic information necessary to extend the microwave characterization of isoprene ozonolysis to prereactive complexes. At the air-water interface, Born-Oppenheimer molecular dynamics simulations indicate that the cycloaddition reaction between ozone and trans-isoprene follows a stepwise mechanism, which is quite distinct from our proposed gas-phase mechanism and occurs on a femtosecond time scale. The stepwise nature of isoprene ozonolysis on the aqueous surface is more consistent with the DeMore mechanism than with the Criegee mechanism suggested by the gas-phase calculations, suggesting that the reaction media may play an important role in the reaction. Overall, these predictions aim to provide a missing fundamental piece of molecular insight into isoprene ozonolysis, which has broad tropospheric implications due to its critical role as a nighttime source of hydroxyl radicals.

Molecular Engineering of Chromophores to Enable Triplet-Triplet Annihilation Upconversion.

Triplet-triplet annihilation upconversion (TTA-UC) is an unconventional photophysical process that yields high-energy photons from low-energy incident light and offers pathways for innovation across many technologies, including solar energy harvesting, photochemistry, and optogenetics. Within aromatic organic chromophores, TTA-UC is achieved through several consecutive energy conversion events that ultimately fuse two triplet excitons into a singlet exciton. In chromophores where the singlet exciton is roughly isoergic with two triplet excitons, the limiting step is the triplet-triplet annihilation pathway, where the kinetics and yield depend sensitively on the energies of the lowest singlet and triplet excited states. Herein we report up to 40-fold improvements in upconversion quantum yields using molecular engineering to selectively tailor the relative energies of the lowest singlet and triplet excited states, enhancing the yield of triplet-triplet annihilation and promoting radiative decay of the resulting singlet exciton. Using this general and effective strategy, we obtain upconversion yields with red emission that are among the highest reported, with remarkable chemical stability under ambient conditions.

The role of high-order electron correlation effects in a model system for non-valence correlation-bound anions.

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H2O)4 cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree-Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger EBEs are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

Electrophotocatalysis with a Trisaminocyclopropenium Radical Dication.

Visible-light photocatalysis and electrocatalysis are two powerful strategies for the promotion of chemical reactions. Here, these two modalities are combined in an electrophotocatalytic oxidation platform. This chemistry employs a trisaminocyclopropenium (TAC) ion catalyst, which is electrochemically oxidized to form a cyclopropenium radical dication intermediate. The radical dication undergoes photoexcitation with visible light to produce an excited-state species with oxidizing power (3.33 V vs. SCE) sufficient to oxidize benzene and halogenated benzenes via single-electron transfer (SET), resulting in C-H/N-H coupling with azoles. A rationale for the strongly oxidizing behavior of the photoexcited species is provided, while the stability of the catalyst is rationalized by a particular conformation of the cis-2,6-dimethylpiperidine moieties.

Toward Benchmark-Quality Ab Initio Predictions for 3d Transition Metal Electrocatalysts: A Comparison of CCSD(T) and ph-AFQMC.

Generating accurate ab initio ionization energies for transition metal complexes is an important step toward the accurate computational description of their electrocatalytic reactions. Benchmark-quality data is required for testing existing theoretical methods and developing new ones but is complicated to obtain for many transition metal compounds due to the potential presence of both strong dynamical and static electron correlation. In this regime, it is questionable whether the so-called gold standard, coupled cluster with singles, doubles, and perturbative triples (CCSD(T)), provides the desired level of accuracy─roughly 1-3 kcal/mol. In this work, we compiled a test set of 28 3d metal-containing molecules relevant to homogeneous electrocatalysis (termed 3dTMV) and computed their vertical ionization energies (ionization potentials) with CCSD(T) and phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) in the def2-SVP basis set. A substantial effort has been made to converge away the phaseless bias in the ph-AFQMC reference values. We assess a wide variety of multireference diagnostics and find that spin-symmetry breaking of the CCSD wave function and the PBE0 density functional correlate well with our analysis of multiconfigurational wave functions. We propose quantitative criteria based on symmetry breaking to delineate correlation regimes inside of which appropriately performed CCSD(T) can produce mean absolute deviations from the ph-AFQMC reference values of roughly 2 kcal/mol or less and outside of which CCSD(T) is expected to fail. We also present a preliminary assessment of density functional theory (DFT) functionals on the 3dTMV set.

Bridging physical intuition and hardware efficiency for correlated electronic states: the local unitary cluster Jastrow ansatz for electronic structure.

A prominent goal in quantum chemistry is to solve the molecular electronic structure problem for ground state energy with high accuracy. While classical quantum chemistry is a relatively mature field, the accurate and scalable prediction of strongly correlated states found, e.g., in bond breaking and polynuclear transition metal compounds remains an open problem. Within the context of a variational quantum eigensolver, we propose a new family of ansatzes which provides a more physically appropriate description of strongly correlated electrons than a unitary coupled cluster with single and double excitations (qUCCSD), with vastly reduced quantum resource requirements. Specifically, we present a set of local approximations to the unitary cluster Jastrow wavefunction motivated by Hubbard physics. As in the case of qUCCSD, exactly computing the energy scales factorially with system size on classical computers but polynomially on quantum devices. The local unitary cluster Jastrow ansatz removes the need for SWAP gates, can be tailored to arbitrary qubit topologies (e.g., square, hex, and heavy-hex), and is well-suited to take advantage of continuous sets of quantum gates recently realized on superconducting devices with tunable couplers. The proposed family of ansatzes demonstrates that hardware efficiency and physical transparency are not mutually exclusive; indeed, chemical and physical intuition regarding electron correlation can illuminate a useful path towards hardware-friendly quantum circuits.

Caught in the act: real-time observation of the solvent response that promotes excited-state proton transfer in pyranine.

Photo-induced excited-state proton transfer (ESPT) reactions are of central importance in many biological and chemical processes. Identifying mechanistic details of the solvent reorganizations that facilitate proton transfer however, is challenging for current experimental and theoretical approaches. Using optical pump THz probe (OPTP) spectroscopy and molecular dynamics simulations, we were able to elucidate the ultrafast changes in the solvation environment for three derivatives of pyranine: the photoacid HPTS, the methoxy derivative MPTS, and the photobase OPTS. Experimentally, we find damped oscillations in the THz signal at short times and our simulations enable their assignment to vibrational energy transfer beatings between the photoexcited chromophore and nearby solvent molecules. The simulations of HPTS reveal strikingly efficient sub-ps energy transfer into a particular solvent mode, that is active near 4 THz, and which can provide the requisite energy required for solvent reorganization promoting proton transfer. Similar oscillations are present in the THz signal for all three derivatives, however the signal is damped rapidly for HPTS (within 0.4 ps) and more slowly for MPTS (within 1.4 ps) and OPTS (within 2.0 ps). For HPTS, we also characterize an additional phonon-like propagation of the proton into the bulk with a 140 ps period and an 83 ps damping time. Thermalization of the solvent occurs on a time scale exceeding 120 ps.

Revisiting the Orbital Energy-Dependent Regularization of Orbital-Optimized Second-Order Møller-Plesset Theory.

Optimizing orbitals in the presence of electron correlation, as in orbital-optimized second-order Møller-Plesset perturbation theory (OOMP2), can remove artifacts associated with mean-field orbitals such as spin contamination and artificial symmetry-breaking. However, OOMP2 is known to suffer from divergent correlation energies in regimes of small orbital energy gaps. To address this issue, several approaches to amplitude regularization have been explored, with those featuring energy-gap-dependent regularizers appearing to be most transferable and physically justifiable. For instance, κ-OOMP2 was shown to address the energy divergence issue in, for example, bond-breaking processes while offering a significant improvement in accuracy for the W4-11 thermochemistry data set, and a parameter of κ = 1.45 was recommended. A more recent investigation of regularized MP2 with Hartree-Fock orbitals revealed that stronger regularization (i.e., smaller values of κ) than what had previously been recommended for κ-OOMP2 may offer huge improvements in certain cases such as noncovalent interactions while retaining a high level of accuracy for main-group thermochemistry data sets. In this study, we investigate the transferability of those findings to κ-OOMP2 and assess the implications of stronger regularization on the ability of κ-OOMP2 to diagnose strong static correlation. We found similar results using κ-OOMP2 for several main-group thermochemistry, barrier height, and noncovalent interaction data sets including both closed shell and open shell species. However, stronger regularization yielded substantially higher accuracy for open-shell transition-metal (TM) thermochemistry and is necessary to provide qualitatively correct spin symmetry breaking behavior for several large and electrochemically relevant TM systems. We therefore find a single κ value insufficient to treat all systems using κ-OOMP2.

A Localized-Orbital Energy Evaluation for Auxiliary-Field Quantum Monte Carlo.

Phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) has recently emerged as a promising method for the production of benchmark-level simulations of medium- to large-sized molecules because of its accuracy and favorable polynomial scaling with system size. Unfortunately, the memory footprints of standard energy evaluation algorithms are nontrivial, which can significantly impact timings on graphical processing units (GPUs) where memory is limited. Previous attempts to reduce scaling by taking advantage of the low-rank structure of the Coulombic integrals have been successful but exhibit high prefactors, making their utility limited to very large systems. Here we present a complementary cubic-scaling route to reduce memory and computational scaling based on the low rank of the Coulombic interactions between localized orbitals, focusing on the application to ph-AFQMC. We show that the error due to this approximation, which we term localized-orbital AFQMC (LO-AFQMC), is systematic and controllable via a single variable and that the method is computationally favorable even for small systems. We present results demonstrating robust retention of accuracy versus both experiment and full ph-AFQMC for a variety of test cases chosen for their potential difficulty for localized-orbital-based methods, including the singlet-triplet gaps of the polyacenes benzene through pentacene, the heats of formation for a set of Platonic hydrocarbon cages, and the total energy of ferrocene, Fe(Cp)2. Finally, we reproduce our previous result for the gas-phase ionization energy of Ni(Cp)2, agreeing with full ph-AFQMC to within statistical error while using less than 1/15th of the computer time.

Calculation of Metallocene Ionization Potentials via Auxiliary Field Quantum Monte Carlo: Toward Benchmark Quantum Chemistry for Transition Metals.

The accurate ab initio prediction of ionization energies is essential to understanding the electrochemistry of transition metal complexes in both materials science and biological applications. However, such predictions have been complicated by the scarcity of gas phase experimental data, the relatively large size of the relevant molecules, and the presence of strong electron correlation effects. In this work, we apply all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multideterminant trial wave functions to six metallocene complexes to compare the computed adiabatic and vertical ionization energies with experimental results. We find that ph-AFQMC yields mean absolute errors (MAEs) of 1.69 ± 1.02 kcal/mol for the adiabatic energies and 2.85 ± 1.13 kcal/mol for the vertical energies. We also carry out density functional theory (DFT) calculations using a variety of functionals, which yields MAEs of 3.62-6.98 kcal/mol and 3.31-9.88 kcal/mol, as well as one variant of localized coupled cluster calculations (DLPNO-CCSD(T0) with moderate PNO cutoffs), which has MAEs of 4.96 and 6.08 kcal/mol, respectively. We also test the reliability of DLPNO-CCSD(T0) and DFT on acetylacetonate (acac) complexes for adiabatic energies measured in the same manner experimentally, and we find higher MAEs, ranging from 4.56 to 10.99 kcal/mol (with a different ordering) for DFT and 6.97 kcal/mol for DLPNO-CCSD(T0). Finally, by utilizing experimental solvation energies, we show that accurate reduction potentials in solution for the metallocene series can be obtained from the AFQMC gas phase results.

Concerted Electron-Nuclear Motion in Proton-Coupled Electron Transfer-Driven Grotthuss-Type Proton Translocation.

Photoinduced proton-coupled electron transfer and long-range two-proton transport via a Grotthuss-type mechanism are investigated in a biomimetic construct. The ultrafast, nonequilibrium dynamics are assessed via two-dimensional electronic vibrational spectroscopy, in concert with electrochemical and computational techniques. A low-frequency mode is identified experimentally and found to promote double proton and electron transfer, supported by recent theoretical simulations of a similar but abbreviated (non-photoactive) system. Excitation frequency peak evolution and center line slope dynamics show direct evidence of strongly coupled nuclear and electronic degrees of freedom, from which we can conclude that the double proton and electron transfer processes are concerted (up to an uncertainty of 24 fs). The nonequilibrium pathway from the photoexcited Franck-Condon region to the E2PT state is characterized by an ∼110 fs time scale. This study and the tools presented herein constitute a new window into hot charge transfer processes involving an electron and multiple protons.

Regularized Second-Order Møller-Plesset Theory: A More Accurate Alternative to Conventional MP2 for Noncovalent Interactions and Transition Metal Thermochemistry for the Same Computational Cost.

Second-order Møller-Plesset theory (MP2) notoriously breaks down for π-driven dispersion interactions and dative bonds in transition metal complexes. Herein, we investigate three physically justified forms of single-parameter, energy-gap dependent regularization which can yield high and transferable accuracy for a variety of noncovalent interactions (including S22, S66, and L7 test sets) and (mostly closed shell) transition metal thermochemistry. Regularization serves to damp overestimated pairwise additive contributions, renormalizing first-order amplitudes such that the effects of higher-order correlations are incorporated. The optimal parameter values for the noncovalent and transition metal sets are 1.1, 0.7, and 0.4 for κ, σ, and σ2 regularizers, respectively. However, such regularization slightly degrades the accuracy of conventional MP2 for some small-molecule test sets, most of which have relatively large average frontier energy gaps. Our results suggest that appropriately regularized MP2 models may improve double hybrid density functionals, at no additional cost over conventional MP2.

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile.

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

Solvent Mediated Excited State Proton Transfer in Indigo Carmine.

Excited state proton transfer (ESPT) is thought to be responsible for the photostability of biological molecules, including DNA and proteins, and natural dyes such as indigo. However, the mechanistic role of the solvent interaction in driving ESPT is not well understood. Here, the electronic excited state deactivation dynamics of indigo carmine (InC) is mapped by visible pump-infrared probe and two-dimensional electronic-vibrational (2DEV) spectroscopy and complemented by electronic structure calculations. The observed dynamics reveal notable differences between InC in a protic solvent, D2O, and an aprotic solvent, deuterated dimethyl sulfoxide (dDMSO). Notably, an acceleration in the excited state decay is observed in D2O (<10 ps) compared to dDMSO (130 ps). Our data reveals clear evidence for ESPT in D2O accompanied by a significant change in dipole moment, which is found not to occur in dDMSO. We conclude that the ability of protic solvents to form intermolecular H-bonds with InC enables ESPT, which facilitates a rapid nonradiative S1 → S0 transition via the monoenol intermediate.

Predicting Ligand-Dissociation Energies of 3d Coordination Complexes with Auxiliary-Field Quantum Monte Carlo.

Transition-metal complexes are ubiquitous in biology and chemical catalysis, yet they remain difficult to accurately describe with ab initio methods because of the presence of a large degree of dynamic electron correlation, and, in some cases, strong static correlation which results from a manifold of low-lying states. Progress has been hindered by a scarcity of high-quality gas-phase experimental data, while exact ab initio predictions are usually computationally unaffordable because of the large size of the relevant complexes. In this work, we present a data set of 34 tetrahedral, square planar, and octahedral 3d metal-containing complexes with gas-phase ligand-dissociation energies that have reported uncertainties of ≤2 kcal/mol. We perform all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) calculations utilizing multideterminant trial wave functions selected by a black box procedure. We compare the results with those from the density functional theory (DFT) with the B3LYP, B97, M06, PBE0, ωB97X-V, and DSD-PBEP86/2013 functionals and a localized orbital variant of the coupled cluster theory with single, double, and perturbative triple excitations (DLPNO-CCSD(T)). We find mean averaged errors of 1.07 ± 0.27 kcal/mol for our most sophisticated ph-AFQMC approach versus 2.81 kcal/mol for DLPNO-CCSD(T) and 1.49-3.78 kcal/mol for DFT. We find maximum errors of 2.96 ± 1.71 kcal/mol for our best ph-AFQMC method versus 9.15 kcal/mol for DLPNO-CCSD(T) and 5.98-13.69 kcal/mol for DFT. The reasonable performance of a number of DFT functionals is in stark contrast to the much poorer accuracy previously demonstrated for diatomic species, suggesting a moderation in electron correlation because of ligand coordination in most cases. However, the unpredictably large errors for a small subset of cases with both DFT and DLPNO-CCSD(T) methods leave cause for concern, especially in light of the unreliability of common multireference indicators. In contrast, the robust and, in principle, systematically improvable results of ph-AFQMC for these realistic complexes establish the method as a useful tool for elucidating the electronic structure of transition-metal-containing complexes and predicting their gas-phase properties.