Fundamental insight into electrochemical oxidation of methane towards methanol on transition metal oxides.

Electrochemical oxidation of CH4 is known to be inefficient in aqueous electrolytes. The lower activity of methane oxidation reaction (MOR) is primarily attributed to the dominant oxygen evolution reaction (OER) and the higher barrier for CH4 activation on transition metal oxides (TMOs). However, a satisfactory explanation for the origins of such lower activity of MOR on TMOs, along with the enabling strategies to partially oxidize CH4 to CH3OH, have not been developed yet. We report here the activation of CH4 is governed by a previously unrecognized consequence of electrostatic (or Madelung) potential of metal atom in TMOs. The measured binding energies of CH4 on 12 different TMOs scale linearly with the Madelung potentials of the metal in the TMOs. The MOR active TMOs are the ones with higher CH4 binding energy and lower Madelung potential. Out of 12 TMOs studied here, only TiO2, IrO2, PbO2, and PtO2 are active for MOR, where the stable active site is the O on top of the metal in TMOs. The reaction pathway for MOR proceeds primarily through *CH x intermediates at lower potentials and through *CH3OH intermediates at higher potentials. The key MOR intermediate *CH3OH is identified on TiO2 under operando conditions at higher potential using transient open-circuit potential measurement. To minimize the overoxidation of *CH3OH, a bimetallic Cu2O3 on TiO2 catalysts is developed, in which Cu reduces the barrier for the reaction of *CH3 and *OH and facilitates the desorption of *CH3OH. The highest faradaic efficiency of 6% is obtained using Cu-Ti bimetallic TMO.

Contrasting Mechanisms of Aromatic and Aryl-Methyl Substituent Hydroxylation by the Rieske Monooxygenase Salicylate 5-Hydroxylase.

The hydroxylase component (S5HH) of salicylate-5-hydroxylase catalyzes C5 ring hydroxylation of salicylate but switches to methyl hydroxylation when a C5 methyl substituent is present. The use of 18O2 reveals that both aromatic and aryl-methyl hydroxylations result from monooxygenase chemistry. The functional unit of S5HH comprises a nonheme Fe(II) site located 12 Šacross a subunit boundary from a one-electron reduced Rieske-type iron-sulfur cluster. Past studies determined that substrates bind near the Fe(II), followed by O2 binding to the iron to initiate catalysis. Stopped-flow-single-turnover reactions (STOs) demonstrated that the Rieske cluster transfers an electron to the iron site during catalysis. It is shown here that fluorine ring substituents decrease the rate constant for Rieske electron transfer, implying a prior reaction of an Fe(III)-superoxo intermediate with a substrate. We propose that the iron becomes fully oxidized in the resulting Fe(III)-peroxo-substrate-radical intermediate, allowing Rieske electron transfer to occur. STO using 5-CD3-salicylate-d8 occurs with an inverse kinetic isotope effect (KIE). In contrast, STO of a 1:1 mixture of unlabeled and 5-CD3-salicylate-d8 yields a normal product isotope effect. It is proposed that aromatic and aryl-methyl hydroxylation reactions both begin with the Fe(III)-superoxo reaction with a ring carbon, yielding the inverse KIE due to sp2 → sp3 carbon hybridization. After Rieske electron transfer, the resulting Fe(III)-peroxo-salicylate intermediate can continue to aromatic hydroxylation, whereas the equivalent aryl-methyl intermediate formation must be reversible to allow the substrate exchange necessary to yield a normal product isotope effect. The resulting Fe(III)-(hydro)peroxo intermediate may be reactive or evolve through a high-valent iron intermediate to complete the aryl-methyl hydroxylation.

Huzinaga projection embedding for efficient and accurate energies of systems with localized spin-densities.

We demonstrate the accuracy and efficiency of the restricted open-shell and unrestricted formulation of the absolutely localized Huzinaga projection operator embedding method. Restricted open-shell and unrestricted Huzinaga projection embedding in the full system basis is formally exact to restricted open-shell and unrestricted Kohn-Sham density functional theory, respectively. By utilizing the absolutely localized basis, we significantly improve the efficiency of the method while maintaining high accuracy. Furthermore, the absolutely localized basis allows for high accuracy open-shell wave function methods to be embedded into a closed-shell density functional theory environment. The open-shell embedding method is shown to calculate electronic energies of a variety of systems to within 1 kcal/mol accuracy of the full system wave function result. For certain highly localized reactions, such as spin transition energies on transition metals, we find that very few atoms are necessary to include in the wave function region in order to achieve the desired accuracy. This extension further broadens the applicability of our absolutely localized Huzinaga level-shift projection operator method to include open-shell species. Here, we apply our method to several representative examples, such as spin splitting energies, catalysis on transition metals, and radical reactions.

Hydrogen Evolving Activity of Dithiolene-Based Metal-Organic Frameworks with Mixed Cobalt and Iron Centers.

Electrocatalytic systems based on metal-organic frameworks (MOFs) have attracted great attention due to their potential application in commercially viable renewable energy-converting devices. We have recently shown that the cobalt 2,3,6,7,10,11-triphenylenehexathiolate (CoTHT) framework can catalyze the hydrogen evolution reaction (HER) in fully aqueous media with Tafel slopes as low as 71 mV/dec and near-unity Faradaic efficiency (FE). Taking advantage of the high synthetic tunability of MOFs, here, we synthesize a series of iron and mixed iron/cobalt THT-based MOFs. The incorporation of the iron and cobalt dithiolene moieties is verified by various spectroscopic techniques, and the integrity of the crystalline structure is maintained regardless of the stoichiometries of the two metals. The hydrogen evolving activity of the materials was explored in pH 1.3 aqueous electrolyte solutions. Unlike CoTHT, the FeTHT framework exhibits minimal activity due to a late catalytic onset [-0.440 V versus reversible hydrogen electrode (RHE)] and a large Tafel slope (210 mV/dec). The performance of the mixed-metal MOFs is adversely affected by the incorporation of Fe, where increasing Fe content results in MOFs with lower HER activity and diminished long-term stability and FE for H2 production. It is proposed that the FeTHT domains undergo alternative Faradaic processes under catalytic conditions, which alter its local structure and electrochemical behavior, eventually resulting in a material with diminished HER performance.

Cp2Ti(κ2-tBuNCNtBu): A Complex with an Unusual κ2 Coordination Mode of a Heterocumulene Featuring a Free Carbene.

Although there are myriad binding modes of heterocumulenes to metal centers, the monometallic κ2-ECE (E = O, S, NR) coordination mode has not been reported. Herein, the synthesis, isolation, and physical characterization of Cp2Ti(κ2-tBuNCNtBu) (3) (Cp = cyclopentadienyl, tBu = tert-butyl), a strained 4-membered metallacycle bearing a free carbene, is described. Computational (DFT, CASSCF, QT-AIM, ELF) and solid-state CP-MAS 13C NMR spectroscopic analysis indicate that 3 is best described as a free carbene with partial Ti-Cß bonding that results from Ti-N π-bonding mixing with N-C-N σ-bonding of the bent N-C-N framework. Reactivity studies of 3 corroborate its carbene-like nature: protonation with [LutH]I results in the formation of a Ti-formamidinate (4), while oxidation with S8 yields a Ti-thioureate (5). Additionally, a related bridged dititanamidocarbene, (Cp2Ti)2(µ-η1,η1-CyNCNCy) (10) (Cy = cyclohexyl) is reported. Taken together, this work suggests that the 2-electron reduction of heterocumulene moieties can allow access to unusual free carbene coordination geometries given the proper stabilizing coordination environment from the reducing transition metal fragment.

Mechanistic insights into electrochemical reduction of CO2 over Ag using density functional theory and transport models.

Electrochemical reduction of CO2 using renewable sources of electrical energy holds promise for converting CO2 to fuels and chemicals. Since this process is complex and involves a large number of species and physical phenomena, a comprehensive understanding of the factors controlling product distribution is required. While the most plausible reaction pathway is usually identified from quantum-chemical calculation of the lowest free-energy pathway, this approach can be misleading when coverages of adsorbed species determined for alternative mechanism differ significantly, since elementary reaction rates depend on the product of the rate coefficient and the coverage of species involved in the reaction. Moreover, cathode polarization can influence the kinetics of CO2 reduction. Here, we present a multiscale framework for ab initio simulation of the electrochemical reduction of CO2 over an Ag(110) surface. A continuum model for species transport is combined with a microkinetic model for the cathode reaction dynamics. Free energies of activation for all elementary reactions are determined from density functional theory calculations. Using this approach, three alternative mechanisms for CO2 reduction were examined. The rate-limiting step in each mechanism is **COOH formation at higher negative potentials. However, only via the multiscale simulation was it possible to identify the mechanism that leads to a dependence of the rate of CO formation on the partial pressure of CO2 that is consistent with experiments. Simulations based on this mechanism also describe the dependence of the H2 and CO current densities on cathode voltage that are in strikingly good agreement with experimental observation.

The Stability of the Complex and the Basicity of the Anion Impact the Selectivity and Affinity of Tripodal Gadolinium Complexes for Anions.

The affinities and selectivities of lanthanide complexes with open coordination sites for anions vary considerably with the chelate. In order to determine the effect of the stability of a lanthanide complex on its affinity for anions, five different complexes featuring different bidentate chelating moieties were synthesized, and their affinity for anions in water at neutral pH were evaluated by longitudinal relaxometry measurements. The chelates comprise both oxygen and nitrogen donors including maltol, 1,2-hydroxypyridinone, hydroxamic acid, pyridin-2-ylmethanol, and carbamoylmethylphosphonate diester. They were chosen to span a range of basicities all the while maintaining a similar tripodal tris-bidentate architecture, thereby allowing for a direct study of the role of the coordinating motif on the supramolecular recognition of anions by the corresponding GdIII complex. Overall, for ligands containing the same number of protonation steps, and therefore the same charge at neutral pH, the lower the acidity of the chelate (higher ∑pKa's), the less stable the corresponding GdIII complex, and the higher its affinity for anions. Regardless of the number of protonation steps, the more stable GdIII complexes form ternary or quaternary assemblies with coordinating anions. In contrast, the same anions readily displace the chelate of the least stable complexes, resulting instead in the formation of GdIII·anion precipitates. Irrespective of the chelate, in the absence of steric hindrance at the open coordination site, the affinity of GdIII complexes for anions follows the order phosphate > arsenate > bicarbonate > fluoride. Hence, the selectivity and affinity of GdIII complexes of tripodal tris-bidentate chelates for anions is a function of the stability of the GdIII complex and the basicity of the anion.

Mechanism of Ti-Catalyzed Oxidative Nitrene Transfer in [2 + 2 + 1] Pyrrole Synthesis from Alkynes and Azobenzene.

A combined computational and experimental study on the mechanism of Ti-catalyzed formal [2 + 2 + 1] pyrrole synthesis from alkynes and aryl diazenes is reported. This reaction proceeds through a formally TiII/TiIV redox catalytic cycle as determined by natural bond orbital (NBO) and intrinsic bond orbital (IBO) analysis. Kinetic analysis of the reaction of internal alkynes with azobenzene reveals a complex equilibrium involving Ti═NPh monomer/dimer equilibrium and Ti═NPh + alkyne [2 + 2] cycloaddition equilibrium along with azobenzene and pyridine inhibition equilibria prior to rate-determining second alkyne insertion. Computations support this kinetic analysis, provide insights into the structure of the active species in catalysis and the roles of solvent, and provide a new mechanism for regeneration of the Ti imido catalyst via disproportionation. Reductive elimination from a 6-membered azatitanacyclohexadiene species to generate pyrrole-bound TiII is surprisingly facile and occurs through a unique electrocyclic reductive elimination pathway similar to a Nazarov cyclization. The resulting TiII species are stabilized through backbonding into the π* of the pyrrole framework, although solvent effects also significantly stabilize free TiII species that are required for pyrrole loss and catalytic turnover. Further computational and kinetic analysis reveals that in complex reactions with unysmmetric alkynes the resulting pyrrole regioselectivity is driven primarily by steric effects for terminal alkynes and inductive effects for internal alkynes.

Accurate and systematically improvable density functional theory embedding for correlated wavefunctions.

We analyze the sources of error in quantum embedding calculations in which an active subsystem is treated using wavefunction methods, and the remainder using density functional theory. We show that the embedding potential felt by the electrons in the active subsystem makes only a small contribution to the error of the method, whereas the error in the nonadditive exchange-correlation energy dominates. We test an MP2 correction for this term and demonstrate that the corrected embedding scheme accurately reproduces wavefunction calculations for a series of chemical reactions. Our projector-based embedding method uses localized occupied orbitals to partition the system; as with other local correlation methods, abrupt changes in the character of the localized orbitals along a reaction coordinate can lead to discontinuities in the embedded energy, but we show that these discontinuities are small and can be systematically reduced by increasing the size of the active region. Convergence of reaction energies with respect to the size of the active subsystem is shown to be rapid for all cases where the density functional treatment is able to capture the polarization of the environment, even in conjugated systems, and even when the partition cuts across a double bond.

CatEmbed: A Machine-Learned Representation Obtained via Categorical Entity Embedding for Predicting Adsorption and Reaction Energies on Bimetallic Alloy Surfaces.

Machine-learning models for predicting adsorption energies on metallic surfaces often rely on basic elemental properties and electronic and geometric descriptors. Here, we apply categorical entity embedding, a featurization method inspired by natural language processing techniques, to predict adsorption energies on bimetallic alloy surfaces using categorical descriptors. Using this method, we develop a machine-learned representation from categorical descriptors (e.g., surface composition, adsorbate type, and site type) of the slab/adsorbate complex. By combining this representation with numerical features (e.g., slab metal stoichiometric ratios), we create the CatEmbed representation. Remarkably, decision tree models trained using CatEmbed, which includes no explicit geometric information, achieve a Mean Absolute Error (MAE) of 0.12 eV. Additionally, we extend this technique to predict reaction energies on bimetallic surfaces, creating the CatEmbed-React representation, which achieves an MAE of 0.08 eV. These findings highlight the effectiveness of categorical entity embedding for predicting adsorption and reaction energies on bimetallic alloy surfaces.

Accurate basis set truncation for wavefunction embedding.

Density functional theory (DFT) provides a formally exact framework for performing embedded subsystem electronic structure calculations, including DFT-in-DFT and wavefunction theory-in-DFT descriptions. In the interest of efficiency, it is desirable to truncate the atomic orbital basis set in which the subsystem calculation is performed, thus avoiding high-order scaling with respect to the size of the MO virtual space. In this study, we extend a recently introduced projection-based embedding method [F. R. Manby, M. Stella, J. D. Goodpaster, and T. F. Miller III, J. Chem. Theory Comput. 8, 2564 (2012)] to allow for the systematic and accurate truncation of the embedded subsystem basis set. The approach is applied to both covalently and non-covalently bound test cases, including water clusters and polypeptide chains, and it is demonstrated that errors associated with basis set truncation are controllable to well within chemical accuracy. Furthermore, we show that this approach allows for switching between accurate projection-based embedding and DFT embedding with approximate kinetic energy (KE) functionals; in this sense, the approach provides a means of systematically improving upon the use of approximate KE functionals in DFT embedding.

Density functional theory embedding for correlated wavefunctions: improved methods for open-shell systems and transition metal complexes.

Density functional theory (DFT) embedding provides a formally exact framework for interfacing correlated wave-function theory (WFT) methods with lower-level descriptions of electronic structure. Here, we report techniques to improve the accuracy and stability of WFT-in-DFT embedding calculations. In particular, we develop spin-dependent embedding potentials in both restricted and unrestricted orbital formulations to enable WFT-in-DFT embedding for open-shell systems, and develop an orbital-occupation-freezing technique to improve the convergence of optimized effective potential calculations that arise in the evaluation of the embedding potential. The new techniques are demonstrated in applications to the van-der-Waals-bound ethylene-propylene dimer and to the hexa-aquairon(II) transition-metal cation. Calculation of the dissociation curve for the ethylene-propylene dimer reveals that WFT-in-DFT embedding reproduces full CCSD(T) energies to within 0.1 kcal/mol at all distances, eliminating errors in the dispersion interactions due to conventional exchange-correlation (XC) functionals while simultaneously avoiding errors due to subsystem partitioning across covalent bonds. Application of WFT-in-DFT embedding to the calculation of the low-spin/high-spin splitting energy in the hexaaquairon(II) cation reveals that the majority of the dependence on the DFT XC functional can be eliminated by treating only the single transition-metal atom at the WFT level; furthermore, these calculations demonstrate the substantial effects of open-shell contributions to the embedding potential, and they suggest that restricted open-shell WFT-in-DFT embedding provides better accuracy than unrestricted open-shell WFT-in-DFT embedding due to the removal of spin contamination.

Embedded density functional theory for covalently bonded and strongly interacting subsystems.

Embedded density functional theory (e-DFT) is used to describe the electronic structure of strongly interacting molecular subsystems. We present a general implementation of the Exact Embedding (EE) method [J. Chem. Phys. 133, 084103 (2010)] to calculate the large contributions of the nonadditive kinetic potential (NAKP) in such applications. Potential energy curves are computed for the dissociation of Li(+)-Be, CH(3)-CF(3), and hydrogen-bonded water clusters, and e-DFT results obtained using the EE method are compared with those obtained using approximate kinetic energy functionals. In all cases, the EE method preserves excellent agreement with reference Kohn-Sham calculations, whereas the approximate functionals lead to qualitative failures in the calculated energies and equilibrium structures. We also demonstrate an accurate pairwise approximation to the NAKP that allows for efficient parallelization of the EE method in large systems; benchmark calculations on molecular crystals reveal ideal, size-independent scaling of wall-clock time with increasing system size.

Improving and Understanding the Hydrogen Evolving Activity of a Cobalt Dithiolene Metal-Organic Framework.

Despite the promising previous reports on the development of electrocatalytic dithiolene-based metal-organic frameworks (MOFs) for the hydrogen evolution reaction (HER), these materials often display poor reproducibility of the HER performance because of their poor bulk properties upon integration with electrode materials. We demonstrate here an in-depth investigation of the electrocatalytic HER activity of a cobalt 2,3,6,7,10,11-triphenylenehexathiolate (CoTHT) MOF. To enhance the durability and charge transport properties of the constructed CoTHT/electrode architecture, CoTHT is deposited as an ink composite (1) composed of Nafion and carbon black. We leverage here the well-established use of catalyst inks in the literature to increase adhesion of the catalyst to the electrode surface and to improve the overall electrical conductivity of the integrated catalyst/electrode. The utilization of the composite 1 leads to a significant improvement in the overpotential (η) to reach a current density of 10 mA/cm2 (η = 143 mV) compared to prior reports, resulting in the most active MOF-based electrocatalyst for the HER that contains only earth-abundant elements. Extensive density functional theory (DFT) calculations were applied to understand the structure of CoTHT and the mechanistic pathways of the HER. The computational results suggest that an AB stacking geometry is energetically favorable, where one layer is slipped by 1.6 Šrelative to the neighboring one along the a and b vectors. Additionally, the DFT calculations indicate that the catalytic cycle likely involves a Volmer discharge step to generate a cobalt hydride, followed by a Heyrovsky step to form a cobalt-H2 intermediate, and finally the dissociation of H2.

Magnetic ordering through itinerant ferromagnetism in a metal-organic framework.

Materials that combine magnetic order with other desirable physical attributes could find transformative applications in spintronics, quantum sensing, low-density magnets and gas separations. Among potential multifunctional magnetic materials, metal-organic frameworks, in particular, bear structures that offer intrinsic porosity, vast chemical and structural programmability, and the tunability of electronic properties. Nevertheless, magnetic order within metal-organic frameworks has generally been limited to low temperatures, owing largely to challenges in creating a strong magnetic exchange. Here we employ the phenomenon of itinerant ferromagnetism to realize magnetic ordering at TC = 225 K in a mixed-valence chromium(II/III) triazolate compound, which represents the highest ferromagnetic ordering temperature yet observed in a metal-organic framework. The itinerant ferromagnetism proceeds through a double-exchange mechanism, which results in a barrierless charge transport below the Curie temperature and a large negative magnetoresistance of 23% at 5 K. These observations suggest applications for double-exchange-based coordination solids in the emergent fields of magnetoelectrics and spintronics.

Exact nonadditive kinetic potentials for embedded density functional theory.

We describe an embedded density functional theory (DFT) protocol in which the nonadditive kinetic energy component of the embedding potential is treated exactly. At each iteration of the Kohn-Sham equations for constrained electron density, the Zhao-Morrison-Parr constrained search method for constructing Kohn-Sham orbitals is combined with the King-Handy expression for the exact kinetic potential. We use this formally exact embedding protocol to calculate ionization energies for a series of three- and four-electron atomic systems, and the results are compared to embedded DFT calculations that utilize the Thomas-Fermi (TF) and the Thomas-Fermi-von Weisacker approximations to the kinetic energy functional. These calculations illustrate the expected breakdown due to the TF approximation for the nonadditive kinetic potential, with errors of 30%-80% in the calculated ionization energies; by contrast, the exact protocol is found to be accurate and stable. To significantly improve the convergence of the new protocol, we introduce a density-based switching function to map between the exact nonadditive kinetic potential and the TF approximation in the region of the nuclear cusp, and we demonstrate that this approximation has little effect on the accuracy of the calculated ionization energies. Finally, we describe possible extensions of the exact protocol to perform accurate embedded DFT calculations in large systems with strongly overlapping subsystem densities.

Mechanistic Studies of Bioorthogonal ATP Analogues for Assessment of Histidine Kinase Autophosphorylation.

Phosphorylation is an essential protein modification and is most commonly associated with hydroxyl-containing amino acids via an adenosine triphosphate (ATP) substrate. The last decades have brought greater appreciation to the roles that phosphorylation of myriad amino acids plays in biological signaling, metabolism, and gene transcription. Histidine phosphorylation occurs in both eukaryotes and prokaryotes but has been shown to dominate signaling networks in the latter due to its role in microbial two-component systems. Methods to investigate histidine phosphorylation have lagged behind those to study serine, threonine, and tyrosine modifications due to its inherent instability and the historical view that this protein modification was rare. An important strategy to overcome the reactivity of phosphohistidine is the development of substrate-based probes with altered chemical properties that improve modification longevity but that do not suffer from poor recognition or transfer by the protein. Here, we present combined experimental and computational studies to better understand the molecular requirements for efficient histidine phosphorylation by comparison of the native kinase substrate, ATP, and alkylated ATP derivatives. While recognition of the substrates by the histidine kinases is an important parameter for the formation of phosphohistidine derivatives, reaction sterics also affect the outcome. In addition, we found that stability of the resulting phosphohistidine moieties correlates with the stability of their hydrolysis products, specifically with their free energy in solution. Interestingly, alkylation dramatically affects the stability of the phosphohistidine derivatives at very acidic pH values. These results provide critical mechanistic insights into histidine phosphorylation and will facilitate the design of future probes to study enzymatic histidine phosphorylation.

Robust, Accurate, and Efficient: Quantum Embedding Using the Huzinaga Level-Shift Projection Operator for Complex Systems.

Using wave function (WF) in density functional theory (DFT) embedding methods provides a framework for performing localized, high-accuracy WF calculations on a system, while not incurring the full computational cost of the WF calculation on the full system. To effectively partition a system into localized WF and DFT subsystems, we utilize the Huzinaga level-shift projection operator within an absolutely localized basis. In this work, we study the ability of the absolutely localized Huzinaga level-shift projection operator method to study complex WF and DFT partitions, including partitions between multiple covalent bonds, a double bond, and transition-metal-ligand bonds. We find that our methodology can accurately describe all of these complex partitions. Additionally, we study the robustness of this method with respect to the WF method, specifically where the embedded systems were described using a multiconfigurational WF method. We found that the method is systematically improvable with respect to both the number of atoms in the WF region and the size of the basis set used, with energy errors less than 1 kcal/mol. Additionally, we calculated the adsorption energy of H2 to a model of an iron metal-organic framework (Fe-MOF-74) to within 1 kcal/mol compared to CASPT2 calculations performed on the full model while incurring only a small fraction of the full computational cost. This work demonstrates that the absolutely localized Huzinaga level-shift projection operator method is applicable to very complex systems with difficult electronic structures.

Absolutely Localized Projection-Based Embedding for Excited States.

We present a quantum embedding method that allows for calculation of local excited states embedded in a Kohn-Sham density functional theory (DFT) environment. Projection-based quantum embedding methodologies provide a rigorous framework for performing DFT-in-DFT and wave function in DFT (WF-in-DFT) calculations. The use of absolute localization, where the density of each subsystem is expanded in only the basis functions associated with the atoms of that subsystem, provide improved computationally efficiency for WF-in-DFT calculations by reducing the number of orbitals in the WF calculation. In this work, we extend absolutely localized projection-based quantum embedding to study localized excited states using EOM-CCSD-in-DFT and TDDFT-in-DFT. The embedding results are highly accurate compared to the corresponding canonical EOM-CCSD and TDDFT results on the full system, with TDDFT-in-DFT frequently more accurate than canonical TDDFT. The absolute localization method is shown to eliminate the spurious low-lying excitation energies for charge-transfer states and prevent overdelocalization of excited states. Additionally, we attempt to recover the environment response caused by the electronic excitations in the high-level subsystem using different schemes and compare their accuracy. Finally, we apply this method to the calculation of the excited-state energy of green fluorescent protein and show that we systematically converge to the full system results. Here we demonstrate how this method can be useful in understanding excited states, specifically which chemical moieties polarize to the excitation. This work shows absolutely localized projection-based quantum embedding can treat local electronic excitations accurately and make computationally expensive WF methods applicable to systems beyond current computational limits.

Iterative Supervised Principal Component Analysis Driven Ligand Design for Regioselective Ti-Catalyzed Pyrrole Synthesis.

The rational design of catalysts remains a challenging endeavor within the broader chemical community owing to the myriad variables that can affect key bond-forming events. Designing selective catalysts for any reaction requires an efficient strategy for discovering predictive structure-activity relationships. Herein, we describe the use of iterative supervised principal component analysis (ISPCA) in de novo catalyst design. The regioselective synthesis of 2,5-dimethyl-1,3,4-triphenyl-1H-pyrrole (C) via a Ti-catalyzed formal [2 + 2 +1] cycloaddition of phenylpropyne and azobenzene was targeted as a proof of principle. The initial reaction conditions led to an unselective mixture of all possible pyrrole regioisomers. ISPCA was conducted on a training set of catalysts, and their performance was regressed against the scores from the top three principal components. Component loadings from this PCA space and k-means clustering were used to inform the design of new test catalysts. The selectivity of a prospective test set was predicted in silico using the ISPCA model, and optimal candidates were synthesized and tested experimentally. This data-driven predictive-modeling workflow was iterated, and after only three generations the catalytic selectivity was improved from 0.5 (statistical mixture of products) to over 11 (>90% C) by incorporating 2,6-dimethyl-4-(pyrrolidin-1-yl)pyridine as a ligand. The origin of catalyst selectivity was probed by examining ISPCA variable loadings in combination with DFT modeling, revealing that ligand lability plays an important role in selectivity. A parallel catalyst search using multivariate linear regression (MLR), a popular approach in catalysis informatics, was also conducted in order to compare these strategies in a hypothetical catalyst scouting campaign. ISPCA appears to be more robust and predictive than MLR when sparse training sets are used that are representative of the data available during the early search for an optimal catalyst. The successful development of a highly selective catalyst without resorting to long, stochastic screening processes demonstrates the inherent power of ISPCA in de novo catalyst design and should motivate the general use of ISPCA in reaction development.