Distinguishing homolytic vs heterolytic bond dissociation of phenylsulfonium cations with localized active space methods.

Modeling chemical reactions with quantum chemical methods is challenging when the electronic structure varies significantly throughout the reaction and when electronic excited states are involved. Multireference methods, such as complete active space self-consistent field (CASSCF), can handle these multiconfigurational situations. However, even if the size of the needed active space is affordable, in many cases, the active space does not change consistently from reactant to product, causing discontinuities in the potential energy surface. The localized active space SCF (LASSCF) is a cheaper alternative to CASSCF for strongly correlated systems with weakly correlated fragments. The method is used for the first time to study a chemical reaction, namely the bond dissociation of a mono-, di-, and triphenylsulfonium cation. LASSCF calculations generate smooth potential energy scans more easily than the corresponding, more computationally expensive CASSCF calculations while predicting similar bond dissociation energies. Our calculations suggest a homolytic bond cleavage for di- and triphenylsulfonium and a heterolytic pathway for monophenylsulfonium.

Competitive Valerate Binding Enables RuO2-Mediated Butene Electrosynthesis in Water.

The (non)-Kolbe oxidation of valeric acid, sourced from a hydrolysis product of cellulose, provides a sustainable synthetic route to access value-added products, such as butene. An essential mechanistic step preceding product formation involves the oxidative and decarboxylative cleavage of a C-C bond. Yet, the role of the electrode surface in mediating this oxidative step remains an open question: the electron transfer can occur either via an inner-sphere or outer-sphere mechanism. Here, we report the electrochemical, in situ spectroscopic, computational, and reactivity studies of RuO2-mediated oxidative decarboxylation of valeric acid to butene in aqueous electrolytes. We find that carboxylates bind to RuO2 anode surfaces at potential values where decarboxylation products are observed. Our results are consistent with a reaction scheme where the competitive and catalytic oxygen evolution reaction (OER) is impeded by these bound carboxylate species while these species are inert toward butene formation. Our results implicate an outer-sphere electron transfer mechanism for decarboxylation where the surface chemistry of the RuO2 electrode serves to enable higher non-Kolbe reaction selectivity by suppressing the parasitic OER. Our findings delineate interfacial design principles for selective electrochemical systems that utilize water as the ultimate oxidant for sustainable decarboxylation.

Core Binding Energy Calculations: A Scalable Approach with the Quantum Embedding-Based Equation-of-Motion Coupled-Cluster Method.

We investigated the use of density matrix embedding theory to facilitate the computation of core ionization energies (IPs) of large molecules at the equation-of-motion coupled-cluster singles doubles with perturbative triples (EOM-CCSD*) level in combination with the core-valence separation (CVS) approximation. The unembedded IP-CVS-EOM-CCSD* method with a triple-ζ basis set produced ionization energies within 1 eV of experiment with a standard deviation of ∼0.2 eV for the core65 data set. The embedded variant contributed very little systematic error relative to the unembedded method, with a mean unsigned error of 0.07 eV and a standard deviation of ∼0.1 eV, in exchange for accelerating the calculations by many orders of magnitude. By employing embedded EOM-CC methods, we computed the core ionization energies of the uracil hexamer, doped fullerene, and chlorophyll molecule, utilizing up to ∼4000 basis functions within 1 eV from experimental values. Such calculations are not currently possible with the unembedded EOM-CC method.

A US perspective on closing the carbon cycle to defossilize difficult-to-electrify segments of our economy.

Electrification to reduce or eliminate greenhouse gas emissions is essential to mitigate climate change. However, a substantial portion of our manufacturing and transportation infrastructure will be difficult to electrify and/or will continue to use carbon as a key component, including areas in aviation, heavy-duty and marine transportation, and the chemical industry. In this Roadmap, we explore how multidisciplinary approaches will enable us to close the carbon cycle and create a circular economy by defossilizing these difficult-to-electrify areas and those that will continue to need carbon. We discuss two approaches for this: developing carbon alternatives and improving our ability to reuse carbon, enabled by separations. Furthermore, we posit that co-design and use-driven fundamental science are essential to reach aggressive greenhouse gas reduction targets.

Automatic State Interaction with Large Localized Active Spaces for Multimetallic Systems.

The localized active space self-consistent field method factorizes a complete active space wave function into an antisymmetrized product of localized active space wave function fragments. Correlation between fragments is then reintroduced through localized active space state interaction (LASSI), in which the Hamiltonian is diagonalized in a model space of LAS states. However, the optimal procedure for defining the LAS fragments and LASSI model space is unknown. We here present an automated framework to explore systematically convergent sets of model spaces, which we call LASSI[r, q]. This method requires the user to select only r, the number of electron hops from one fragment to another, and q, the number of fragment basis functions per Hilbert space, which converges to CASCI in the limit of r, q → ∞. Numerical tests of this method on the trimetal oxo-centered complexes [Fe(III)Al(III)Fe(II)(µ3-O)(HCOO)6] and [Fe(III)2Fe(II)(µ3-O)(HCOO)6] show efficient convergence to the CASCI limit with 4-10 orders of magnitude fewer states than CASCI.

Aliovalent Substitution Tunes Physical Properties in a Conductive Bis(dithiolene) Two-Dimensional Metal-Organic Framework.

Two-dimensional conductive metal-organic frameworks have emerged as promising electronic materials for applications in (opto)electronic, thermoelectric, magnetic, electrocatalytic, and energy storage devices. Many bottom-up or postsynthetic protocols have been developed to isolate these materials or further modulate their electronic properties. However, some methodologies commonly used in classic semiconductors, notably, aliovalent substitution, are conspicuously absent. Here, we demonstrate how aliovalent Fe(III) to Ni(II) substitution enables the isolation of a Ni bis(dithiolene) material from a previously reported Fe analogue. Detailed characterization supports the idea that aliovalent substitution of Fe(III) to Ni(II) results in an in situ oxidation of the organic dithiolene linker. This substitution-induced redox tuning modulates the electronic properties in the system, leading to higher electrical conductivity and Hall mobility but slightly lower carrier densities and weaker antiferromagnetic interactions. Moreover, this aliovalent substitution improves the material's electrochemical stability and thus enables pseudocapacitive behavior in the Ni material. These results demonstrate how classic aliovalent substitution strategies in semiconductors can also be leveraged in conductive MOFs and add further support to this class of compounds as emerging electronic materials.

Catalytic, Spectroscopic, and Theoretical Studies of Fe4S4-Based Coordination Polymers as Heterogenous Coupled Proton-Electron Transfer Mediators for Electrocatalysis.

Iron-sulfur clusters play essential roles in biological systems, and thus synthetic [Fe4S4] clusters have been an area of active research. Recent studies have demonstrated that soluble [Fe4S4] clusters can serve as net H atom transfer mediators, improving the activity and selectivity of a homogeneous Mn CO2 reduction catalyst. Here, we demonstrate that incorporating these [Fe4S4] clusters into a coordination polymer enables heterogeneous H atom transfer from an electrode surface to a Mn complex dissolved in solution. A previously reported solution-processable Fe4S4-based coordination polymer was successfully deposited on the surfaces of different electrodes. The coated electrodes serve as H atom transfer mediators to a soluble Mn CO2 reduction catalyst displaying good product selectivity for formic acid. Furthermore, these electrodes are recyclable with a minimal decrease in activity after multiple catalytic cycles. The heterogenization of the mediator also enables the characterization of solution-phase and electrode surface species separately. Surface enhanced infrared absorption spectroscopy (SEIRAS) reveals spectroscopic signatures for an in situ generated active Mn-H species, providing a more complete mechanistic picture for this system. The active species, reaction mechanism, and the protonation sites on the [Fe4S4] clusters were further confirmed by density functional theory calculations. The observed H atom transfer reactivity of these coordination polymer-coated electrodes motivates additional applications of this composite material in reductive H atom transfer electrocatalysis.

State Preparation in Quantum Algorithms for Fragment-Based Quantum Chemistry.

State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space-unitary coupled cluster (LAS-UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. In this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C-C double bond breaking in transbutadiene and the spin ladder in a bimetallic system. We analyze the impact of QPE parameters, such as the number of ancilla qubits and Trotter steps, on the prepared state. We find a trade-off between the methods, where DI requires fewer resources for smaller fragments, while QPE is more efficient for larger fragments. Our resource estimates highlight the benefits of system fragmentation in state preparation for subsequent quantum chemical calculations. These findings have broad applications for preparing multireference quantum chemical wave functions on quantum circuits that can be used for realistic chemical applications.

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory.

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear-electronic coupling in addition to accurately predicting Franck-Condon vertical excitations. In this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

Organic Reactivity Made Easy and Accurate with Automated Multireference Calculations.

In organic reactivity studies, quantum chemical calculations play a pivotal role as the foundation of understanding and machine learning model development. While prevalent black-box methods like density functional theory (DFT) and coupled-cluster theory (e.g., CCSD(T)) have significantly advanced our understanding of chemical reactivity, they frequently fall short in describing multiconfigurational transition states and intermediates. Achieving a more accurate description necessitates the use of multireference methods. However, these methods have not been used at scale due to their often-faulty predictions without expert input. Here, we overcome this deficiency with automated multiconfigurational pair-density functional theory (MC-PDFT) calculations. We apply this method to 908 automatically generated organic reactions. We find 68% of these reactions present significant multiconfigurational character in which the automated multiconfigurational approach often provides a more accurate and/or efficient description than DFT and CCSD(T). This work presents the first high-throughput application of automated multiconfigurational methods to reactivity, enabled by automated active space selection algorithms and the computation of electronic correlation with MC-PDFT on-top functionals. This approach can be used in a black-box fashion, avoiding significant active space inconsistency error in both single- and multireference cases and providing accurate multiconfigurational descriptions when needed.

Minimum-Energy Conical Intersections by Compressed Multistate Pair-Density Functional Theory.

Compressed multistate pair-density functional theory (CMS-PDFT) is a multistate version of multiconfiguration pair-density functional theory that can capture the correct topology of coupled potential energy surfaces (PESs) around conical intersections. In this work, we develop interstate coupling vectors (ISCs) for CMS-PDFT in the OpenMolcas and PySCF/mrh electronic structure packages. Yet, the main focus of this work is using ISCs to calculate minimum-energy conical intersections (MECIs) by CMS-PDFT. This is performed using the projected constrained optimization method in OpenMolcas, which uses ISCs to restrain the iterations to the conical intersection seam. We optimize the S1/S0 MECIs for ethylene, butadiene, and benzene and show that CMS-PDFT gives smooth PESs in the vicinities of the MECIs. Furthermore, the CMS-PDFT MECIs are in good agreement with the MECI calculated by the more expensive XMS-CASPT2 method.

Harvesting Water from Air with High-Capacity, Stable Furan-Based Metal-Organic Frameworks.

We synthesized two isoreticular furan-based metal-organic frameworks (MOFs), MOF-LA2-1(furan) and MOF-LA2-2(furan) with rod-like secondary building units (SBUs) featuring 1D channels, as sorbents for atmospheric water harvesting (LA = long arm). These aluminum-based MOFs demonstrated a combination of high water uptake and stability, exhibiting working capacities of 0.41 and 0.48 gwater/gMOF (under isobaric conditions of 1.70 kPa), respectively. Remarkably, both MOFs showed a negligible loss in water uptake after 165 adsorption-desorption cycles. These working capacities rival that of MOF-LA2-1(pyrazole), which has a working capacity of 0.55 gwater/gMOF. The current MOFs stand out for their high water stability, as evidenced by 165 cycles of water uptake and release. MOF-LA2-2(furan) is the first aluminum MOF to employ a double 'long arm' extension strategy, which is confirmed through single-crystal X-ray diffraction (SCXRD). The MOFs were synthesized by using a straightforward synthesis route. This study offers valuable insights into the design of durable, water-stable MOFs and underscores their potential for efficient water harvesting.

Divergent Bimetallic Mechanisms in Copper(II)-Mediated C-C, N-N, and O-O Oxidative Coupling Reactions.

Copper-catalyzed aerobic oxidative coupling of diaryl imines provides a route for conversion of ammonia to hydrazine. The present study uses experimental and density functional theory computational methods to investigate the mechanism of N-N bond formation, and the data support a mechanism involving bimolecular coupling of Cu-coordinated iminyl radicals. Computational analysis is extended to CuII-mediated C-C, N-N, and O-O coupling reactions involved in the formation of cyanogen (NC-CN) from HCN, 1,3-butadiyne from ethyne (i.e., Glaser coupling), hydrazine from ammonia, and hydrogen peroxide from water. The results reveal two different mechanistic pathways. Heteroatom ligands with an uncoordinated lone pair (iminyl, NH2, OH) undergo charge transfer to CuII, generating ligand-centered radicals that undergo facile bimolecular radical-radical coupling. Ligands lacking a lone pair (CN and CCH) form bridged binuclear diamond-core structures that undergo C-C coupling. This mechanistic bifurcation is rationalized by analysis of spin densities in key intermediates and transition states, as well as multiconfigurational calculations. Radical-radical coupling is especially favorable for N-N coupling owing to energetically favorable charge transfer in the intermediate and thermodynamically favorable product formation.