Deconvoluting Nonlinear Catalyst-Substrate Effects in the Intramolecular Dirhodium-Catalyzed C-H Insertion of Donor/Donor Carbenes Using Data Science Tools.

Interactions between catalysts and substrates can be highly complex and dynamic, often complicating the development of models to either predict or understand such processes. A dirhodium(II)-catalyzed C-H insertion of donor/donor carbenes into 2-alkoxybenzophenone substrates to form benzodihydrofurans was selected as a model system to explore nonlinear methods to achieve a mechanistic understanding. We found that the application of traditional methods of multivariate linear regression (MLR) correlating DFT-derived descriptors of catalysts and substrates leads to poorly performing models. This inspired the introduction of nonlinear descriptor relationships into modeling by applying the sure independence screening and sparsifying operator (SISSO) algorithm. Based on SISSO-generated descriptors, a high-performing MLR model was identified that predicts external validation points well. Mechanistic interpretation was aided by the deconstruction of feature relationships using chemical space maps, decision trees, and linear descriptors. Substrates were found to have a strong dependence on steric effects for determining their innate cyclization selectivity preferences. Catalyst reactive site features can then be matched to product features to tune or override the resultant diastereoselectivity within the substrate-dictated ranges. This case study presents a method for understanding complex interactions often encountered in catalysis by using nonlinear modeling methods and linear deconvolution by pattern recognition.

Leveraging Regio- and Stereoselective C(sp3)-H Functionalization of Silyl Ethers to Train a Logistic Regression Classification Model for Predicting Site-Selectivity Bias.

The C-H functionalization of silyl ethers via carbene-induced C-H insertion represents an efficient synthetic disconnection strategy. In this work, site- and stereoselective C(sp3)-H functionalization at α, γ, δ, and even more distal positions to the siloxy group has been achieved using donor/acceptor carbene intermediates. By exploiting the predilections of Rh2(R-TCPTAD)4 and Rh2(S-2-Cl-5-BrTPCP)4 catalysts to target either more electronically activated or more spatially accessible C-H sites, respectively, divergent desired products can be formed with good diastereocontrol and enantiocontrol. Notably, the reaction can also be extended to enable desymmetrization of meso silyl ethers. Leveraging the broad substrate scope examined in this study, we have trained a machine learning classification model using logistic regression to predict the major C-H functionalization site based on intrinsic substrate reactivity and catalyst propensity for overriding it. This model enables prediction of the major product when applying these C-H functionalization methods to a new substrate of interest. Applying this model broadly, we have demonstrated its utility for guiding late-stage functionalization in complex settings and developed an intuitive visualization tool to assist synthetic chemists in such endeavors.

Toggling the Z-type interaction off-on in nickel-boron dihydrogen and anionic hydride complexes.

Completing a series of nickel-group 13 complexes, a coordinatively unsaturated nickel-boron complex and its derivatives with a H2, N2, or hydride ligand were synthesized and characterized. The toggling "on" of a Ni(0)-B(III) inverse-dative bond enabled the stabilization of a nickel-bound anionic hydride with a remarkably low thermodynamic hydricity of kcal mol-1 in THF. The flexible topology of the boron metalloligand confers both favorable hydrogen binding affinity and strong hydride donicity, albeit at the cost of high H2 basicity during deprotonation to form the hydride.

Mechanistically Guided Workflow for Relating Complex Reactive Site Topologies to Catalyst Performance in C-H Functionalization Reactions.

Leveraging congested catalyst scaffolds has emerged as a key strategy for altering innate substrate site-selectivity profiles in C-H functionalization reactions. Similar to enzyme active sites, optimal small molecule catalysts often feature reactive cavities tailored for controlling substrate approach trajectories. However, relating three-dimensional catalyst shape to reaction output remains a formidable challenge, in part due to the lack of molecular features capable of succinctly describing complex reactive site topologies in terms of numerical inputs for machine learning applications. Herein, we present a new set of descriptors, "Spatial Molding for Approachable Rigid Targets" (SMART), which we have applied to quantify reactive site spatial constraints for an expansive library of dirhodium catalysts and to predict site-selectivity for C-H functionalization of 1-bromo-4-pentylbenzene via donor/acceptor carbene intermediates. Optimal site-selectivity for the terminal methylene position was obtained with Rh2(S-2-Cl-5-MesTPCP)4 (30.9:1 rr, 14:1 dr, 87% ee), while C-H functionalization at the electronically activated benzylic site was increasingly favored for Rh2(TPCP)4 catalysts lacking an ortho-Cl, Rh2(S-PTAD)4, and Rh2(S-TCPTAD)4, respectively. Intuitive global site-selectivity models for 25 disparate dirhodium catalysts were developed via multivariate linear regression to explicitly assess the contributing roles of steric congestion and dirhodium-carbene electrophilicity in controlling the site of C-H functionalization. The workflow utilizes spatial classification to extract descriptors only for reactive catalyst conformers, a nuance that may be widely applicable for establishing close correspondence between ground-state model systems and transition states. Broader still, SMART descriptors are amenable for delineating salient reactive site features to predict reactivity in other chemical and biological contexts.

Bioinspired Nickel Complexes Supported by an Iron Metalloligand.

Nature utilizes multimetallic sites in metalloenzymes to enable multielectron chemical transformations at ambient conditions and low overpotentials. One such example of multimetallic cooperativity can be found in the C-cluster of Ni-carbon monoxide dehydrogenase (CODH), which interconverts CO and CO2. Toward a potential functional model of the C-cluster, a family of Ni-Fe bimetallic complexes was synthesized that contain direct metal-metal bonding interactions. The complexes were characterized by X-ray crystallography, various spectroscopies (NMR, EPR, UV-vis, Mössbauer), and theoretical calculations. The Ni-Fe bimetallic system has a reversible Fe(III)/Fe(II) redox couple at -2.10 V (vs Fc+/Fc). The Fe-based "redox switch" can turn on CO2 reactivity at the Ni(0) center by leveraging the Ni→Fe dative interaction to attenuate the Ni(0) electron density. The reduced Ni(0)Fe(II) species mediated the formal two-electron reduction of CO2 to CO, providing a Ni-CO adduct and CO32- as products. During the reaction, an intermediate was observed that is proposed to be a Ni-CO2 species.

Thermodynamic and kinetic studies of H2 and N2 binding to bimetallic nickel-group 13 complexes and neutron structure of a Ni(η2-H2) adduct.

Understanding H2 binding and activation is important in the context of designing transition metal catalysts for many processes, including hydrogenation and the interconversion of H2 with protons and electrons. This work reports the first thermodynamic and kinetic H2 binding studies for an isostructural series of first-row metal complexes: NiML, where M = Al (1), Ga (2), and In (3), and L = [N(o-(NCH2PiPr2)C6H4)3]3-. Thermodynamic free energies (ΔG°) and free energies of activation (ΔG ‡) for binding equilibria were obtained via variable-temperature 31P NMR studies and lineshape analysis. The supporting metal exerts a large influence on the thermodynamic favorability of both H2 and N2 binding to Ni, with ΔG° values for H2 binding found to span nearly the entire range of previous reports. The non-classical H2 adduct, (η2-H2)NiInL (3-H2), was structurally characterized by single-crystal neutron diffraction-the first such study for a Ni(η2-H2) complex or any d10 M(η2-H2) complex. UV-Vis studies and TD-DFT calculations identified specific electronic structure perturbations of the supporting metal which poise NiML complexes for small-molecule binding. ETS-NOCV calculations indicate that H2 binding primarily occurs via H-H σ-donation to the Ni 4p z -based LUMO, which is proposed to become energetically accessible as the Ni(0)→M(iii) dative interaction increases for the larger M(iii) ions. Linear free-energy relationships are discussed, with the activation barrier for H2 binding (ΔG ‡) found to decrease proportionally for more thermodynamically favorable equilibria. The ΔG° values for H2 and N2 binding to NiML complexes were also found to be more exergonic for the larger M(iii) ions.

Formal Nickelate(-I) Complexes Supported by Group†13 Ions.

Formal nickelate(-I) complexes bearing Group 13 metalloligands (M=Al and Ga) were isolated. These 17 e- complexes were synthesized by one-electron reduction of the corresponding Ni0 →MIII precursors, and were investigated by single-crystal X-ray diffraction, EPR spectroscopy, and quantum chemical calculations. Collectively, the experimental and computational data support: 1) the strengthening of the Ni-M bond upon one-electron reduction, and 2) the delocalization of the unpaired spin across the Ni and M atoms. An intriguing electronic configuration is revealed where three valence electrons occupy two σ-type bonding interactions: Ni(3dz2 )2 →M and σ-(Ni-M)1 . The latter is an unusual Ni-M σ-bonding molecular orbital that comprises primarily the Ni 4pz and M npz /ns atomic orbitals.

A Bimetallic Nickel-Gallium Complex Catalyzes CO2 Hydrogenation via the Intermediacy of an Anionic d10 Nickel Hydride.

Large-scale CO2 hydrogenation could offer a renewable stream of industrially important C1 chemicals while reducing CO2 emissions. Critical to this opportunity is the requirement for inexpensive catalysts based on earth-abundant metals instead of precious metals. We report a nickel-gallium complex featuring a Ni(0)→Ga(III) bond that shows remarkable catalytic activity for hydrogenating CO2 to formate at ambient temperature (3150 turnovers, turnover frequency = 9700 h-1), compared with prior homogeneous Ni-centered catalysts. The Lewis acidic Ga(III) ion plays a pivotal role in stabilizing catalytic intermediates, including a rare anionic d10 Ni hydride. Structural and in situ characterization of this reactive intermediate support a terminal Ni-H moiety, for which the thermodynamic hydride donor strength rivals those of precious metal hydrides. Collectively, our experimental and computational results demonstrate that modulating a transition metal center via a direct interaction with a Lewis acidic support can be a powerful strategy for promoting new reactivity paradigms in base-metal catalysis.

Heterobimetallic Complexes That Bond Vanadium to Iron, Cobalt, and Nickel.

Zero-valent iron, cobalt, and nickel were installed into the metalloligand V[N(o-(NCH2P((i)Pr)2)C6H4)3] (1, VL), generating the heterobimetallic trio FeVL (2), CoVL (3), and NiVL (4), respectively. In addition, the one-electron-oxidized analogues [FeVL]X ([2(ox)]X, where X(-) = BPh4 or PF6) and [CoVL]BPh4 ([3(ox)]BPh4) were prepared. The complexes were characterized by a host of physical methods, including cyclic voltammetry, X-ray crystallography, magnetic susceptibility, electronic absorption, NMR, electron paramagnetic resonance (EPR), and Mössbauer spectroscopies. The CoV and FeV heterobimetallic compounds have short M-V bond lengths that are consistent with M-M multiple bonding. As revealed by theoretical calculations, the M-V bond is triple in 2, 2(ox), and 3(ox), double in 3, and dative (Ni → V) in 4. The (d-d)(10) species, 2 and 3(ox), are diamagnetic and exhibit large diamagnetic anisotropies of -4700 × 10(-36) m(3)/molecule. Complexes 2 and 3(ox) are also characterized by intense visible bands at 760 and 610 nm (ε > 1000 M(-1) cm(-1)), respectively, which correspond to an intermetal (M → V) charge-transfer transition. Magnetic susceptibility measurements and EPR characterization establish S = (1)/2 ground states for (d-d)(9) 2(ox) and (d-d)(11) 3, while (d-d)(12) 4 is S = 1 based on Evans' method.

Tuning Nickel with Lewis Acidic Group 13 Metalloligands for Catalytic Olefin Hydrogenation.

A series of bimetallic complexes pairing zero-valent nickel with group 13 M(III) ions is reported. Stronger Ni→M(III) dative bonds that render Ni more electron-deficient are seen for larger ions (In > Ga > Al). The larger Ga and In ions stabilize rare, nonclassical Ni-H2 adducts that catalyze olefin hydrogenation. In contrast, neither the Ni-Al complex nor a single nickel center enables H2 binding or olefin hydrogenation. By comparison of the structures, redox properties, and catalytic activities of the Ni-M series, the electronic and steric effects of the supporting metal ion are elucidated.