The Nature of the Barrier to Phosphane Dissociation from Grubbs Olefin Metathesis Catalysts

Assessment of Density Functional Theory Methods for the Structural Prediction of Transition and Post-Transition Metal-Nucleic Acid Complexes

Understanding the structure of metal-nucleic acid systems is important for many applications such as the design of new pharmaceuticals, metal detection platforms, and nanomaterials. Herein, we explore the ability of 20 density functional theory (DFT) functionals to reproduce the crystal structure geometry of transition and post-transition metal-nucleic acid complexes identified in the Protein Data Bank and Cambridge Structural Database. The environmental extremes of the gas phase and implicit water were considered, and analysis focused on the global and inner coordination geometry, including the coordination distances. Although gas-phase calculations were unable to describe the structure of 12 out of the 53 complexes in our test set regardless of the DFT functional considered, accounting for the broader environment through implicit solvation or constraining the model truncation points to crystallographic coordinates generally afforded agreement with the experimental structure, suggesting that functional performance for these systems is likely due to the models rather than the methods. For the remaining 41 complexes, our results show that the reliability of functionals depends on the metal identity, with the magnitude of error varying across the periodic table. Furthermore, minimal changes in the geometries of these metal-nucleic acid complexes occur upon use of the Stuttgart-Dresden effective core potential and/or inclusion of an implicit water environment. The overall top three performing functionals are ωB97X-V, ωB97X-D3(BJ), and MN15, which reliably describe the structure of a broad range of metal-nucleic acid systems. Other suitable functionals include MN15-L, which is a cheaper alternative to MN15, and PBEh-3c, which is commonly used in QM/MM calculations of biomolecules. In fact, these five methods were the only functionals tested to reproduce the coordination sphere of Cu²⁺-containing complexes. For metal-nucleic acid systems that do not contain Cu²⁺, ωB97X and ωB97X-D are also suitable choices. These top-performing methods can be utilized in future investigations of diverse metal-nucleic acid complexes of relevance to biology and material science.

Solv: An Alternative Continuum Model Implementation Based on Fixed Atomic Charges, Scaled Particle Theory, and the Atom-Atom Potential Method

An alternative continuum model implementation is reported. The electrostatic contribution to the solvation Gibbs free energy utilizes the noniterative conductor-like screening model of Vyboishchikov and Voityuk (DOI: 10.1002/jcc.26531) based on the fixed partial atomic charges. The nonelectrostatic solute-solvent dispersion-repulsion energy is calculated through the Caillet-Claverie atom-atom potential method employing the grid-based approach. The nonelectrostatic cavitation energy is computed within the scaled particle theory (SPT) formalism with the solute hard-sphere radius obtained via the Pierotti-Claverie (PC) scheme, from the solute molecular surface (SPT-S) or volume (SPT-V). The solvent hard-sphere radius is derived through the fitting to the experimental total solvation free energies of 2530 neutral species in 92 solvents. Application of the model to reproduce both the absolute and relative (reaction net) solvation free energies indicates that the SPT-V approach based on the CM5 charges is the best performer. The method is suggested for the solvation free energy calculation in the nonaqueous solvents.

Formation of active species from ruthenium alkylidene catalysts-an insight from computational perspective

Ruthenium alkylidene complexes are commonly used as olefin metathesis catalysts. Initiation of the catalytic process requires formation of a 14-electron active ruthenium species via dissociation of a respective ligand. In the present work, this initiation step has been computationally studied for the Grubbs-type catalysts (H₂IMes)(PCy₃)(Cl)₂Ru=CHPh, (H₂IMes)(PCy₃)(Cl)₂Ru=CH-CH=CMe₂ and (H₂IMes)(3-Br-py)₂(Cl)₂Ru=CHPh, and the Hoveyda-Grubbs-type catalysts (H₂IMes)(Cl)₂Ru=CH(o-OiPrC₆H₄), (H₂IMes)(Cl)₂Ru=CH(5-NO₂-2-OiPrC₆H₃), and (H₂IMes)(Cl)₂Ru=CH(2-OiPr-3-PhC₆H₃), using density functional theory (DFT). Additionally, the extended-transition-state combined with the natural orbitals for the chemical valence (ETS-NOCV) and the interacting quantum atoms (IQA) energy decomposition methods were applied. The computationally determined activity order within both families of the catalysts and the activation parameters are in agreement with reported experimental data. The significance of solvent simulation and the basis set superposition error (BSSE) correction is discussed. ETS-NOCV demonstrates that the bond between the dissociating ligand and the Ru-based fragment is largely ionic followed by the charge delocalizations: σ(Ru-P) and π(Ru-P) and the secondary CH^…Cl, CH^…π, and CH^…HC interactions. In the case of transition state structures, the majority of stabilization stems from London dispersion forces exerted by the efficient CH^…Cl, CH^…π, and CH^…HC interactions. Interestingly, the height of the electronic dissociation barriers is, however, directly connected with the prevalent (unfavourable) changes in the electrostatic and orbital interaction contributions despite the favourable relief in Pauli repulsion and geometry reorganization terms during the activation process. According to the IQA results, the isopropoxy group in the Hoveyda-Grubbs-type catalysts is an efficient donor of intra-molecular interactions which are important for the activity of these catalysts.

Further understanding of the Ru-centered [2+2] cycloreversion/cycloaddition involved into the interconversion of ruthenacyclobutane using the Grubbs catalysts from a reaction force analysis

The chemical reactivity of the first- and second-generation Grubbs catalysts has always been a significant issue in olefin metathesis. In the present work, we study the [2+2] cycloreversion/cycloaddition and the alkylidene rotation involved into the interconversion of the ruthenacyclobutane intermediate, through the reaction force and reaction force constant analysis. It has been found that the structural contribution controls the barrier energy in the interconversion of ruthenacyclobutane via [2+2] cycloreversion/cycloaddition, which is slightly lower in the second generation of Grubbs catalysts while its electronic contribution is slightly higher, which unveils a major rigidity and donor/acceptor properties of the NHC. This finding explains a greater structural contribution in the rate constant. Moreover, on the basis of the reaction force constant, the process can be classified as "two-stage"-concerted reactions, noting a more asynchronous process when the first generation is used as a catalyst.Finally, a similar analysis into the alkylidene rotation was performed. It was determined that [2+2] cycloreversion and alkylidene rotations take place in a sequential manner, the energy barrier is again controlled by structural reorganization, and the pathway is less asynchronous.

Disentangling Ligand Effects on Metathesis Catalyst Activity: Experimental and Computational Studies of Ruthenium-Aminophosphine Complexes

Second-generation ruthenium olefin metathesis catalysts bearing aminophosphine ligands were investigated with systematic variation of the ligand structure. The rates of phosphine dissociation ( k₁; initiation rate) and relative phosphine reassociation ( k_-1) were determined for two series of catalysts bearing cyclohexyl(morpholino)phosphine and cyclohexyl(piperidino)phosphine ligands. In both cases, incorporating P-N bonds into the architecture of the dissociating phosphine accelerates catalyst initiation relative to the parent [Ru]-PCy₃ complex; however, this effect is muted for the tris(amino)phosphine-ligated complexes, which exhibit higher ligand binding constants in comparison to those with phosphines containing one or two cyclohexyl substituents. These results, along with X-ray crystallographic data and DFT calculations, were used to understand the influence of ligand structure on catalyst activity. Especially noteworthy is the application of phosphines containing incongruent substituents (PR₁R'₂); detailed analyses of factors affecting ligand dissociation, including steric effects, inductive effects, and ligand conformation, are presented. Computational studies of the reaction coordinate for ligand dissociation reveal that ligand conformational changes contribute to the rapid dissociation for the fastest-initiating catalyst of these series, which bears a cyclohexyl-bis(morpholino)phosphine ligand. Furthermore, the effect of amine incorporation was examined in the context of ring-opening metathesis polymerization, and reaction rates were found to correlate well with catalyst initiation rates. The combined experimental and computational studies presented in this report reveal important considerations for designing efficient ruthenium olefin metathesis catalysts.

Ground-State Gas-Phase Structures of Inorganic Molecules Predicted by Density Functional Theory Methods

We tested a battery of density functional theory (DFT) methods ranging from generalized gradient approximation (GGA) via meta-GGA to hybrid meta-GGA schemes as well as Møller-Plesset perturbation theory of the second order and a single and double excitation coupled-cluster (CCSD) theory for their ability to reproduce accurate gas-phase structures of di- and triatomic molecules derived from microwave spectroscopy. We obtained the most accurate molecular structures using the hybrid and hybrid meta-GGA approximations with B3PW91, APF, TPSSh, mPW1PW91, PBE0, mPW1PBE, B972, and B98 functionals, resulting in lowest errors. We recommend using these methods to predict accurate three-dimensional structures of inorganic molecules when intramolecular dispersion interactions play an insignificant role. The structures that the CCSD method predicts are of similar quality although at considerably larger computational cost. The structures that GGA and meta-GGA schemes predict are less accurate with the largest absolute errors detected with BLYP and M11-L, suggesting that these methods should not be used if accurate three-dimensional molecular structures are required. Because of numerical problems related to the integration of the exchange-correlation part of the functional and large scattering of errors, most of the Minnesota models tested, particularly MN12-L, M11, M06-L, SOGGA11, and VSXC, are also not recommended for geometry optimization. When maintaining a low computational budget is essential, the nonseparable gradient functional N12 might work within an acceptable range of error. As expected, the DFT-D3 dispersion correction had a negligible effect on the internuclear distances when combined with the functionals tested on nonweakly bonded di- and triatomic inorganic molecules. By contrast, the dispersion correction for the APF-D functional has been found to shorten the bonds significantly, up to 0.064 Å (AgI), in Ag halides, BaO, BaS, BaF, BaCl, Cu halides, and Li and Na halides and hydrides. These results do not agree well with very accurate structures derived from microwave spectroscopy; we therefore believe that the dispersion correction in the APF-D method should be reconsidered. Finally, we found that inaccurate structures can easily lead to errors of few kcal/mol in single-point energies.

Theoretical characterization of first and second generation Grubbs catalysts in styrene cross-metathesis reactions: insights from conceptual DFT

Reorganization energy and DFT-based reactivity descriptors revealed important issues regarding the performance of Grubbs catalysts.

The divergent effects of strong NHC donation in catalysis

The inverse relationship between NHC donicity and catalyst initiation.

Strong σ-donation from NHC ligands (NHC = N-heterocyclic carbene) is shown to have profoundly conflicting consequences for the reactivity of transition-metal catalysts. Such donation is regarded as central to high catalyst activity in many contexts, of which the second-generation Grubbs metathesis catalysts (RuCl₂(NHC)(PCy₃)( Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> CHPh), GII) offer an early, prominent example. Less widely recognized is the dramatically inhibiting impact of NHC ligation on initiation of GII, and on re-entry into the catalytic cycle from the resting-state methylidene species RuCl₂(NHC)(PCy₃)( Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> CH₂), GIIm. Both GII and the methylidene complexes are activated by dissociation of PCy₃. The impact of NHC donicity on the rate of PCy₃ loss is explored in a comparison of s-GIIm, vs.u-GIIm, in which the NHC ligand is saturated H₂IMes or unsaturated IMes, respectively. PCy₃ loss is nearly an order of magnitude slower for the IMes derivative (a difference that is replicated, albeit smaller, for the benzylidene precatalysts GII). Proposed as an overlooked contributor to these rate differences is an increase in the Ru–PCy₃ bond strength arising from π-back-donation onto the phosphine ligand. Strong σ-donation from the IMes ligand, coupled with the inability of this unsaturated NHC to participate in significant π-backbonding, amplifies Ru → PCy₃ π-back-donation. The resulting increase in Ru–P bond strength greatly inhibits entry into the active cycle. For s-GII, in contrast, the greater π-acceptor capacity of the NHC ligand enables competing Ru → H₂IMes back-donation (as confirmed by NOE experiments, which reveal restricted rotation about the Ru–NHC bond for H₂IMes, but not IMes). Ru → PCy₃ back-donation is thus attenuated in the H₂IMes complexes, accounting for the greater lability of the PCy₃ ligand in s-GIIm and s-GII. Similarly inhibited initiation is predicted for other metal–NHC catalysts in which a π-acceptor ligand L must be dissociated to permit substrate binding. Conversely, enhanced reactivity can be expected where such L ligands are pure σ-donors. These effects are expected to be particularly dramatic where the NHC ligand has minimal π-acceptor capacity (as in the unsaturated Arduengo carbenes), and in geometries that maximize NHC–M–L orbital interactions.

A comprehensive study of olefin metathesis catalyzed by Ru-based catalysts

During a Ru-catalyzed reaction of an olefin with an alkylidene moiety that leads to a metallacycle intermediate, the cis insertion of the olefin can occur from two different directions, namely side and bottom with respect to the phosphine or N-heterocyclic ligand (NHC), depending on the first or second generation Grubbs catalyst. Here, DFT calculations unravel to which extent the bottom coordination of olefins with respect is favored over the side coordination through screening a wide range of catalysts, including first and second generation Grubbs catalysts as well as the subsequent Hoveyda derivatives. The equilibrium between bottom and side coordination is influenced by sterics, electronics, and polarity of the solvent. The side attack is favored for sterically less demanding NHC and/or alkylidene ligands. Moreover the generation of a 14-electron species is also discussed, with either pyridine or phosphine ligands to dissociate.

Key processes in ruthenium-catalysed olefin metathesis

While the fundamental series of [2+2]cycloadditions and retro[2+2]cycloadditions that make up the pathways of ruthenium-catalysed metathesis reactions is well-established, the exploration of mechanistic aspects of alkene metathesis continues. In this Feature Article, modern mechanistic studies of the alkene metathesis reaction, catalysed by well-defined ruthenium complexes, are discussed. Broadly, these concern the processes of pre-catalyst initiation, propagation and decomposition, which all have a considerable impact on the overall efficiency of metathesis reactions.

Theory-assisted development of a robust and Z-selective olefin metathesis catalyst

DFT calculations have predicted a new, highly Z-selective ruthenium-based olefin metathesis catalyst that is considerably more robust than the recently reported (SIMes)(Cl)(RS)RuCH(o-OiPrC6H4) (3a, SIMes = 1,3-dimesityl-4,5-dihydroimidazol-2-ylidene, R = 2,4,6-triphenylbenzene) [J. Am. Chem. Soc., 2013, 135, 3331]. Replacing the chloride of 3a by an isocyanate ligand to give 5a was predicted to increase the stability of the complex considerably, at the same time moderately improving the Z-selectivity. Compound 5a is easily prepared in a two-step synthesis starting from the Hoveyda-Grubbs second-generation catalyst 3. In agreement with the calculations, the isocyanate-substituted 5a appears to be somewhat more Z-selective than the chloride analogue 3a. More importantly, 5a can be used in air, with unpurified and non-degassed substrates and solvents, and in the presence of acids. These are traits that are unprecedented among highly Z-selective olefin metathesis catalysts and also very promising with respect to applications of the new catalyst.

The influence of anionic ligands on stereoisomerism of Ru carbenes and their importance to efficiency and selectivity of catalytic olefin metathesis reactions

Investigations detailed herein provide insight regarding the mechanism of stereochemical inversion of stereogenic-at-Ru carbene complexes through a nonolefin metathesis-based polytopal rearrangement pathway. Computational analyses (DFT) reveal that there are two key factors that generate sufficient energy barriers that are responsible for the possibility of isolation and characterization of high-energy, but kinetically stable, intermediates: (1) donor-donor interactions that involve the anionic ligands and the strongly electron donating carbene groups and (2) dipolar effects arising from the syn relationship between the anionic groups (iodide and phenoxide). We demonstrate that a Brønsted acid lowers barriers to facilitate isomerization, and that the positive influence of a proton source is the result of its ability to diminish the repulsive electronic interactions originating from the anionic ligands. The implications of the present studies regarding a more sophisticated knowledge of the role of anionic units on the efficiency of Ru-catalyzed olefin metathesis reactions are discussed. The electronic basis for the increased facility with which allylic alcohols participate in olefin metathesis processes will be presented as well. Finally, we illustrate how a better understanding of the role of anionic ligands has served as the basis for successful design of Ru-based Z-selective catalysts for alkene metathesis.