Computational Studies of Synthetically Relevant Homogeneous Organometallic Catalysis Involving Ni, Pd, Ir, and Rh: An Overview of Commonly Employed DFT Methods and Mechanistic Insights

An effective and versatile strategy for the synthesis of structurally diverse heteroarylsilanes via Ir(iii)-catalyzed C-H silylation

A versatile silylation of heteroaryl C–H bonds is accomplished under the catalysis of a well-defined spirocyclic NHC Ir(iii) complex (SNIr), generating a variety of heteroarylsilanes. A significant advantage of this catalytic system is that multiple types of intermolecular C–H silylation can be achieved using one catalytic system at α, β, γ, or δ positions of heteroatoms with excellent regioselectivities. Mechanistic experiments and DFT calculations indicate that the polycyclic ligand of SNIr can form an isolable cyclometalated intermediate, which leaves a phenyl dentate free and provides a hemi-open space for activating substrates. In general, favorable silylations occur at γ or δ positions of chelating heteroatoms, forming 5- or 6-membered C–Ir–N cyclic intermediates. If such an activation mode is prohibited sterically, silylations would take place at the α or β positions. The mechanistic studies would be helpful for further explaining the reactivity of the SNIr system.

A versatile silylation of heteroaryl C–H bonds is accomplished under the catalysis of a well-defined spirocyclic NHC Ir(iii) complex (SNIr), generating a variety of heteroarylsilanes.

Five-Membered Cyclic Carbonates: Versatility for Applications in Organic Synthesis, Pharmaceutical, and Materials Sciences

This review presents the recent advances involving several applications of five-membered cyclic carbonates and derivatives. With more than 150 references, it covers the period from 2012 to 2020, with special emphasis on the use of five-membered cyclic carbonates as building blocks for organic synthesis and material elaboration. We demonstrate the application of cyclic carbonates in several important chemical transformations, such as decarboxylation, hydrogenation, and transesterification reactions, among others. The presence of cyclic carbonates in molecules with high biological potential is also displayed, together with the importance of these compounds in the preparation of materials such as urethanes, polyurethanes, and flame retardants.

Mechanism and Selectivity Control in Ni- and Pd-Catalyzed Cross-Couplings Involving Carbon-Oxygen Bond Activation

Transition-metal-catalyzed C-O bond activation provides a useful strategy for utilizing alcohol- and phenol-derived electrophiles in cross-coupling reactions, which has become a research field of active and growing interest in organic chemistry. The synergy between computation and experiment elucidated the mechanistic model and controlling factors of selectivities in these transformations, leading to advances in innovative C-O bond activation and functionalization methods.Toward the rational design of C-O bond activation, our collaborations with the Jarvo group bridged the mechanistic models of C(sp²)-O and C(sp³)-O bond activations. We found that the nickel catalyst cleaves the benzylic and allylic C(sp³)-O bonds via two general mechanisms: the stereoinvertive S_N2 back-side attack model and the stereoretentive chelation-assisted model. These two models control the stereochemistry in a wide array of stereospecific Ni-catalyzed cross-coupling reactions with benzylic or allylic alcohol derivatives. Because of the catalyst distortion, the ligands can differentiate the competing stereospecific C(sp³)-O bond activations. The PCy₃ ligand interacts with nickel mainly through σ-donation, and the Ni(PCy₃) catalyst can undergo facile bending of the substrate-nickel-ligand angle, which favors the stereoretentive benzylic C-O bond activation. The N-heterocyclic carbene SIMes ligand has additional d(metal)-p(ligand) back-donation with nickel, which leads to an extra energy penalty for the same angle bending. This results in the preference of stereoinvertive benzylic C-O bond activation under Ni/SIMes catalysis. In addition to ligand control, a Lewis acid can increase the selectivity for stereoinvertive C(sp³)-O activation by stabilizing the S_N2 back-side attack transition state. The oxygen leaving group complexes with the MgI₂ Lewis acid in the stereoinvertive activation, leading to the exclusive stereoinvertive Kumada coupling of benzylic ethers. We also identified that the competing C(sp³)-O bond activation models have noticeable differences in charge separation. This leads to the solvent polarity control of the stereospecificity in C(sp³)-O activations. Low-polarity solvents favor the neutral stereoretentive C-O bond activation, while high-polarity solvents favor the zwitterionic stereoinvertive cleavage.In sharp contrast to the nickel catalysts, the C(sp²)-O bond activation under palladium catalysis mainly proceeds via the classic three-membered ring oxidative addition mechanism instead of the chelation-assisted mechanism. This is due to the lower oxophilicity of palladium, which disfavors the oxygen coordination in the chelation-assisted-type activation. The three-membered ring activation model selectively cleaves the weak C-O bond, resulting in the exclusive chemoselectivity of acyl C-O bond activation in Pd-catalyzed cross-coupling reactions with aryl carboxylic acid derivatives. This explains the overall acylation in the Pd-catalyzed Suzuki-Miyaura coupling with aryl esters. In collaboration with the Szostak group, we revealed that the three-membered ring model applies in the Pd-catalyzed C-O bond activation of carboxylic acid anhydride, which stimulated the development of a series of Pd-catalyzed decarbonylative functionalizations of aryl carboxylic acids.

Reaction-based machine learning representations for predicting the enantioselectivity of organocatalysts

Hundreds of catalytic methods are developed each year to meet the demand for high-purity chiral compounds. The computational design of enantioselective organocatalysts remains a significant challenge, as catalysts are typically discovered through experimental screening. Recent advances in combining quantum chemical computations and machine learning (ML) hold great potential to propel the next leap forward in asymmetric catalysis. Within the context of quantum chemical machine learning (QML, or atomistic ML), the ML representations used to encode the three-dimensional structure of molecules and evaluate their similarity cannot easily capture the subtle energy differences that govern enantioselectivity. Here, we present a general strategy for improving molecular representations within an atomistic machine learning model to predict the DFT-computed enantiomeric excess of asymmetric propargylation organocatalysts solely from the structure of catalytic cycle intermediates. Mean absolute errors as low as 0.25 kcal mol-1 were achieved in predictions of the activation energy with respect to DFT computations. By virtue of its design, this strategy is generalisable to other ML models, to experimental data and to any catalytic asymmetric reaction, enabling the rapid screening of structurally diverse organocatalysts from available structural information.

Base-Free Cross-Couplings of Aryl Diazonium Salts in Methanol: PdII -Alkoxy as Reactivity-Controlling Intermediate

Pd-catalyzed cross-coupling reactions of aryl diazonium salts are generally assumed to proceed via cationic Pd^II intermediates which in turn would be highly reactive in the subsequent transmetalation step. Contrary to this belief, we herein report our observation and rationalization of opposing reactivities of ArN₂ ⁺ in Suzuki (=effective) and Stille (=ineffective) cross-couplings in MeOH. Our systematic experimental and computational studies on the roles of transmetalating agent, solvent, base and the likely involvement of in situ formed diazoether derivatives challenge the currently accepted mechanism. Our data suggest that the observed solvent dichotomy is primarily due to Pd^II -methoxy intermediates being formed, which are unreactive with arylstannanes, but highly reactive with arylboronic acids, complementing the Suzuki "Pd-oxy" mechanism with the direct demonstration of transmetalation of a Pd^II -alkoxy complex. Lewis acids were found to circumvent this reactivity divergence, promoting efficient couplings regardless of the employed conditions or coupling partners.

Graphene-supported single-atom catalysts and applications in electrocatalysis

Supported metal nanostructures are the most extensively studied heterogeneous catalysts, benefiting from easy separation, regeneration and affordable cost. The size of the supported metal species is one of the decisive factors in determining the activity of heterogeneous catalysts. Particularly, the unsaturated coordination environment of metal atoms preferably act as the active centers, minimizing these metal species can significantly boost the specific activity of every single metal atom. Single-atom catalysts/catalysis (SACs), containing isolated metals atomically dispersed on or coordinated with the surface of a support material, represent the ultimate utilization of supported metals and maximize metal usage efficiency. Graphene, a two-dimensional star material, exhibiting extraordinary physical and chemical properties, has been approved as an excellent platform for constructing SACs. When atomically dispersed metal atoms are strongly anchored on the graphene surface, featuring ultra-high surface area and excellent electronic properties, SACs offer a great potential to significantly innovate the conventional heterogeneous catalysis, especially in the field of electrocatalysis. In this review, a detailed discussion of graphene-supported SACs, including preparation approaches, characterization techniques and applications on typical electrocatalytic reactions is provided. The advantages and unique features of graphene-supported SACs as efficient electrocatalysts and the upcoming challenges for improving their performance and further practical applications are also highlighted.

Asymmetric hydrogenation catalyzed by first-row transition metal complexes

This review focuses on asymmetric direct and transfer hydrogenation with first-row transition metal complexes. The reaction mechanisms and the models of enantiomeric induction were summarized and emphasized.

Theoretical insight into the origins of chemo- and diastereo-selectivity in the palladium-catalysed (3 + 2) cyclisation of 5-alkenyl thiazolones

The mechanism of the palladium-catalysed (3 + 2) cyclisation of 5-alkenyl thiazolones and VECs has been investigated from a computational perspective, and the origins of unique chemoselectivity and excellent diastereoselectivity have been disclosed.

Comparative study of aryl halides in Pd-mediated reactions: key factors beyond the oxidative addition step

The comparative experimental study of Ar–X (X = Cl, Br, I) reactivity and analysis reported herein suggest that oxidative addition cannot be considered the sole reason of the observed low reactivity of aryl chlorides.

Reactions of nickel(0) with organochlorides, organobromides, and organoiodides: mechanisms and structure/reactivity relationships

The reactions of nickel(0) complexes with phosphine, bipyridine-type, and N-heterocyclic carbene ligands with aryl, vinyl, and alkyl halides is reviewed.

Oxidative Addition of Water, Alcohols, and Amines in Palladium Catalysis

The homolytic cleavage of O-H and N-H or weak C-H bonds is a key elementary step in redox catalysis, but is thought to be unfeasible for palladium. In stark contrast, reported here is the room temperature and reversible oxidative addition of water, isopropanol, hexafluoroisopropanol, phenol, and aniline to a palladium(0) complex with a cyclic (alkyl)(amino)carbene (CAAC) and a labile pyridino ligand, as is also the case in popular N-heterocyclic carbene (NHC) palladium(II) precatalysts. The oxidative addition of protic solvents or adventitious water switches the chemoselectivity in catalysis with alkynes through activation of the terminal C-H bond. Most salient, the homolytic activation of alcohols and amines allows atom-efficient, additive-free cross-coupling and transfer hydrogenation under mild reaction conditions with usually unreactive, yet desirable reagents, including esters and bis(pinacolato)diboron.

Understanding the unique reactivity patterns of nickel/JoSPOphos manifold in the nickel-catalyzed enantioselective C-H cyclization of imidazoles

The 3d transition metal-catalyzed enantioselective C–H functionalization provides a sustainable strategy for the construction of chiral molecules. A better understanding of the catalytic nature of the reactions and the factors controlling the enantioselectivity is important for rational design of more efficient systems. Herein, the mechanisms of Ni-catalyzed enantioselective C–H cyclization of imidazoles are investigated by density functional theory (DFT) calculations. Both the π-allyl nickel(ii)-promoted σ-complex-assisted metathesis (σ-CAM) and the nickel(0)-catalyzed oxidative addition (OA) mechanisms are disfavored. In addition to the typically proposed ligand-to-ligand hydrogen transfer (LLHT) mechanism, the reaction can also proceed via an unconventional σ-CAM mechanism that involves hydrogen transfer from the JoSPOphos ligand to the alkene through P–H oxidative addition/migratory insertion, C(sp²)–H activation via σ-CAM, and C–C reductive elimination. Importantly, computational results based on this new mechanism can indeed reproduce the experimentally observed enantioselectivities. Further, the catalytic activity of the π-allyl nickel(ii) complex can be rationalized by the regeneration of the active nickel(0) catalyst via a stepwise hydrogen transfer, which was confirmed by experimental studies. The calculations reveal several significant roles of the secondary phosphine oxide (SPO) unit in JoSPOphos during the reaction. The improved mechanistic understanding will enable design of novel enantioselective C–H transformations.

Several unique reactivity patterns of the Ni/JoSPOphos manifold, including facile hydrogen transfer via the two-step oxidative addition/migratory insertion and C(sp²)–H activation via an unconventional σ-CAM mechanism, were disclosed in this work.

Glycerol Boosted Rh-Catalyzed Hydroaminomethylation Reaction: A Mechanistic Insight

We report a Rh-catalyzed hydroaminomethylation reaction of terminal alkenes in glycerol that proceeds efficiently under mild conditions to produce the corresponding amines in relatively high selectivity towards linear amines, moderate to excellent yields by using a low catalyst loading (1 mol % [Rh], 2 mol % phosphine) and relative low pressure (H₂ /CO, 1:1, total pressure 10 bar). This work sheds light on the importance of glycerol in enabling enamine reduction via hydrogen transfer. Moreover, evidence for the crucial role of Rh as chemoselective catalyst in the condensation step has been obtained for the first time in the frame of the hydroaminomethylation reaction by precluding deleterious aldol condensation reactions. The hydroaminomethylation proceeds under a molecular regime; the outcome of catalytically active species into metal-based nanoparticles renders the catalytic system inactive.

Insights into the Regioselectivity of Hydroheteroarylation of Allylbenzene with Pyridine Catalyzed by Ni/AlMe3 with N-Heterocyclic Carbene: The Concerted Hydrogen Transfer Mechanism

The hydroheteroarylation of allylbenzene with pyridine as catalyzed by Ni/AlMe₃ and a N-heterocyclic carbene ligand has recently been established. Density functional calculations revealed that the common stepwise pathway, which involves the C-H oxidative addition of pyridine-AlMe₃ before the migratory insertion of allylbenzene, is unlikely as the migratory insertion needs to overcome a prohibitively high energy barrier. In contrast, the ligand-to-ligand hydrogen transfer pathway is more favorable in which the hydrogen is transferred directly from the para-position of pyridine-AlMe₃ to C2 of allylbenzene. Our distortion-interaction analysis and natural bond orbital analysis indicate that the interaction energy is strongly correlated with the extent of the charge transfer from the alkene (hydrogen acceptor) to the pyridine-AlMe₃ (hydrogen donor), which dictates the selectivity of the H-transfer to the C2 position of allylbenzene. Then, the subsequent C-C reductive elimination of the regioselective linear product is facilitated by the steric hindrance of the IPr ligand. Understanding these key factors affecting the product regioselectivity is important to the development of catalysts for hydroheteroarylation of alkenes.

Nickel(II)-Catalyzed Borylation of Alkenyl Methyl Ethers via C-O Bond Cleavage

A new protocol has been developed for the borylation of conjugated alkenyl methyl ethers using B₂Pin₂ via C-O bond cleavage catalyzed by Ni(II). In this cross-coupling reaction, both E/Z isomers of alkenyl ethers are converted into (E)-alkenyl boronic esters with good reactivity. This transformation exhibits high chemoselectivity in the presence of competitive C-O bonds such as aryl ether, ester, amide, and thioether groups, thus providing a new method for the construction of various alkenyl boronates.

Rationalizing the AlI -Promoted Oxidative Addition of C-C Versus C-H Bonds in Arenes

The factors controlling the oxidative addition of C-C and C-H bonds in arenes mediated by Al^I have been computationally explored by means of Density Functional Theory calculations. To this end, we compared the processes involving benzene, naphthalene and anthracene which are promoted by a recently prepared anionic Al^I -carbenoid. It is found that this species exhibits a strong tendency to oxidatively activate C-H bonds over C-C bonds, with the notable exception of benzene, where the C-C bond activation is feasible but only under kinetic control reaction conditions. State-of-the-art computational methods based on the combination of the Activation Strain Model of reactivity and the Energy Decomposition Analysis have been used to rationalize the competition between both bond activation reactions as well as to quantitatively analyze in detail the ultimate factors controlling these transformations.

Sulfonate Versus Sulfonate: Nickel and Palladium Multimetallic Cross-Electrophile Coupling of Aryl Triflates with Aryl Tosylates

While phenols are frequent and convenient aryl sources in cross-coupling, typically as sulfonate esters, the direct cross-Ullmann coupling of two different sulfonate esters is unknown. We report here a general solution to this challenge catalyzed by a combination of Ni and Pd with Zn reductant and LiBr as an additive. The reaction has broad scope, as demonstrated in 33 examples (65% ± 11% average yield). Mechanistic studies show that Pd strongly prefers the aryl triflate, the Ni catalyst has a small preference for the aryl tosylate, aryl transfer between catalysts is mediated by Zn, and Pd improves yields by consuming arylzinc intermediates.

Carboxylate breaks the arene C-H bond via a hydrogen-atom-transfer mechanism in electrochemical cobalt catalysis

Combined computational and experimental studies elucidated the distinctive mechanistic features of electrochemical cobalt-catalyzed C-H oxygenation. A sequential electrochemical-chemical (EC) process was identified for the formation of an amidylcobalt(iii) intermediate. The synthesis, characterization, cyclic voltammetry studies, and stoichiometric reactions of the related amidylcobalt(iii) intermediate suggested that a second on-cycle electro-oxidation occurs on the amidylcobalt(iii) species, which leads to a formal Co(iv) intermediate. This amidylcobalt(iv) intermediate is essentially a cobalt(iii) complex with one additional single electron distributed on the coordinating heteroatoms. The radical nature of the coordinating pivalate allows the formal Co(iv) intermediate to undergo a novel carboxylate-assisted HAT mechanism to cleave the arene C-H bond, and a CMD mechanism could be excluded for a Co(iii/i) catalytic scenario. The mechanistic understanding of electrochemical cobalt-catalyzed C-H bond activation highlights the multi-tasking electro-oxidation and the underexplored reaction channels in electrochemical transition metal catalysis.

Machine learning dihydrogen activation in the chemical space surrounding Vaska's complex

Homogeneous catalysis using transition metal complexes is ubiquitously used for organic synthesis, as well as technologically relevant in applications such as water splitting and CO2 reduction. The key steps underlying homogeneous catalysis require a specific combination of electronic and steric effects from the ligands bound to the metal center. Finding the optimal combination of ligands is a challenging task due to the exceedingly large number of possibilities and the non-trivial ligand-ligand interactions. The classic example of Vaska's complex, trans-[Ir(PPh3)2(CO)(Cl)], illustrates this scenario. The ligands of this species activate iridium for the oxidative addition of hydrogen, yielding the dihydride cis-[Ir(H)2(PPh3)2(CO)(Cl)] complex. Despite the simplicity of this system, thousands of derivatives can be formulated for the activation of H2, with a limited number of ligands belonging to the same general categories found in the original complex. In this work, we show how DFT and machine learning (ML) methods can be combined to enable the prediction of reactivity within large chemical spaces containing thousands of complexes. In a space of 2574 species derived from Vaska's complex, data from DFT calculations are used to train and test ML models that predict the H2-activation barrier. In contrast to experiments and calculations requiring several days to be completed, the ML models were trained and used on a laptop on a time-scale of minutes. As a first approach, we combined Bayesian-optimized artificial neural networks (ANN) with features derived from autocorrelation and deltametric functions. The resulting ANNs achieved high accuracies, with mean absolute errors (MAE) between 1 and 2 kcal mol-1, depending on the size of the training set. By using a Gaussian process (GP) model trained with a set of selected features, including fingerprints, accuracy was further enhanced. Remarkably, this GP model minimized the MAE below 1 kcal mol-1, by using only 20% or less of the data available for training. The gradient boosting (GB) method was also used to assess the relevance of the features, which was used for both feature selection and model interpretation purposes. Features accounting for chemical composition, atom size and electronegativity were found to be the most determinant in the predictions. Further, the ligand fragments with the strongest influence on the H2-activation barrier were identified.

Selective ortho-Functionalization of Adamantylarenes Enabled by Dispersion and an Air-Stable Palladium(I) Dimer

Contrary to the general belief that Pd-catalyzed cross-coupling at sites of severe steric hindrance are disfavored, we herein show that the oxidative addition to C-Br ortho to an adamantyl group is as favored as the corresponding adamantyl-free system due to attractive dispersion forces. This enabled the development of a fully selective arylation and alkylation of C-Br ortho to an adamantyl group, even if challenged with competing non-hindered C-OTf or C-Cl sites. The method makes use of an air-stable PdI dimer and enables straightforward access to diversely substituted therapeutically important adamantylarenes in 5-30 min.

On the Nature of C(sp3)-C(sp2) Bond Formation in Nickel-Catalyzed Tertiary Radical Cross-Couplings: A Case Study of Ni/Photoredox Catalytic Cross-Coupling of Alkyl Radicals and Aryl Halides

The merger of photoredox and nickel catalysis has enabled the construction of quaternary centers. However, the mechanism, role of the ligand, and effect of the spin state for this transformation and related Ni-catalyzed cross-couplings involving tertiary alkyl radicals in combination with bipyridine and diketonate ligands remain unknown. Several mechanisms have been proposed, all invoking a key Ni(III) species prior to undergoing irreversible inner-sphere reductive elimination. In this work, we have used open-shell dispersion-corrected DFT calculations, quasi-classical dynamics calculations, and experiments to study in detail the mechanism of carbon-carbon bond formation in Ni bipyridine- and diketonate-based catalytic systems. These calculations revealed that access to high spin states (e.g., triplet spin state tetrahedral Ni(II) species) is critical for effective radical cross-coupling of tertiary alkyl radicals. Further, these calculations revealed a disparate mechanism for the C-C bond formation. Specifically, contrary to the neutral Ni-bipyridyl system, diketonate ligands lead directly to the corresponding tertiary radical cross-coupling products via an outer-sphere reductive elimination step via triplet spin state from the Ni(III) intermediates. Implications to related Ni-catalyzed radical cross-couplings and the design of new transformations are discussed.