Formation of Active Species from Ruthenium Alkylidene Catalysts-An
Total Page:16
File Type:pdf, Size:1020Kb
Journal of Molecular Modeling (2019) 25: 331 https://doi.org/10.1007/s00894-019-4202-5 ORIGINAL PAPER Formation of active species from ruthenium alkylidene catalysts—an insight from computational perspective PawełŚliwa1 & Mariusz P. Mitoraj2 & Filip Sagan2 & Jarosław Handzlik1 Received: 27 July 2019 /Accepted: 3 September 2019 /Published online: 7 November 2019 # The Author(s) 2019 Abstract 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 (H2IMes)(PCy3)(Cl)2Ru=CHPh, (H2IMes)(PCy3)(Cl)2Ru=CH- CH=CMe2 and (H2IMes)(3-Br-py)2(Cl)2Ru=CHPh, and the Hoveyda-Grubbs-type catalysts (H2IMes)(Cl)2Ru=CH(o- OiPrC6H4), (H2IMes)(Cl)2Ru=CH(5-NO2–2-OiPrC6H3), and (H2IMes)(Cl)2Ru=CH(2-OiPr-3-PhC6H3), using density function- al 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 dem- onstrates 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…π,andCH…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. Keywords Density functional theory . Natural orbitals for chemical valence . Extended-transition-state . Grubbs catalysts . Hoveyda-Grubbs catalysts . Olefin metathesis This paper is dedicated to Professor Zdzisław Latajka, honouring his significant contribution to quantum chemistry. Introduction This paper belongs to the Topical Collection Zdzisław Latajka 70th Birthday Festschrift Well-defined ruthenium alkylidene complexes (Fig. 1)are Electronic supplementary material The online version of this article commonly used as highly efficient catalysts for olefin metath- (https://doi.org/10.1007/s00894-019-4202-5) contains supplementary esis. These systems show a high tolerance towards most of the material, which is available to authorized users. functional groups in the substrates, air and moisture insensi- tivity, and thermal stability [1–4]. * Mariusz P. Mitoraj The key step during the olefin metathesis process is initial [email protected] formation of a 14-electron active complex [5–14]. In the case * Jarosław Handzlik of the Grubbs-type catalysts (1–3), it occurs by dissociation of [email protected] ligands L, whereas for the Hoveyda-Grubbs catalysts (4–6), it requires decoordination of the ether moiety (Fig. 2). Electron 1 Faculty of Chemical Engineering and Technology, Cracow University of Technology, ul. Warszawska 24, and steric properties of these ligands significantly affect the 31-155 Kraków, Poland rate of the catalyst activation [6–8, 10–12, 14, 15]. In addition 2 Department of Theoretical Chemistry, Faculty of Chemistry, to this often postulated dissociative mechanism, where the Jagiellonian University, ul. Gronostajowa 2, 30-387 Kraków, Poland dissociation/decoordination is the first step, associative and 331 Page 2 of 15 J Mol Model (2019) 25: 331 Fig. 1 Example ruthenium alkylidene catalysts interchange mechanisms for the activation of the second and rate-determining step [12]. For the Grubbs-type catalysts, the third generation ruthenium alkylidene catalysts were also pro- dissociative mechanism was rather computationally predicted posed [8, 9, 12, 14, 16, 17]. In the case of the associative as the kinetically favoured one, with the dissociation step be- mechanism, the alkene molecule coordinates to the metal cen- ing rate determining [12, 17]. Calculations also showed that tre before dissociation of the appropriate ligand, whereas for the Gibbs energy barrier for phosphine dissociation from the the interchange mechanism, both steps take place simulta- Grubbs complex is solvent- and entropy-related [18]. neously. Experimental kinetic studies for the Hoveyda- In this work, the formation of the propagating 14-electron Grubbs catalysts indicated that two mechanisms, dissociative alkylidene species from the Grubbs-type (1–3) and Hoveyda- and interchange, are in work, depending on the steric and Grubbs-type (4–6) initiators has been investigated using den- electronic properties of the ruthenium complex and the sub- sity functional theory (DFT). Our main objectives are to com- strate employed [8]. Accordingly, the dissociative pathway pare the activation parameters for catalysts with different li- was theoretically shown to be preferred over the interchange gands within a given family of compounds and to analyse in one in the case of the Hoveyda-Grubbs complexes bearing a detail the nature of the interactions between the dissociating large chelating alkoxy ligand, although the Hoveyda ligand ligand L and the rest of the system. Bonding between the decoordination during the cross-metathesis reaction is the ligand L and the remaining Ru-based fragment is R NN N N 1 Mes Mes -L Mes Mes Cl Cl R Ru Ru 1 Cl R +L Cl R L R N N Mes N N Mes Mes Mes Cl Cl Ru Ru Cl Cl R1 R1 N N Mes Mes Mes N N Mes R' Cl Cl Ru Ru O Cl Cl R R 1 O 2 R2 R R Fig. 2 Mechanism of olefin metathesis catalysed by ruthenium alkylidene complexes [7] J Mol Model (2019) 25: 331 Page 3 of 15 331 characterized for catalysts 1–3 and the corresponding transi- It is well-established that the ground state for 16- tion state structures, using the extended-transition-state meth- electron and 14-electron ruthenium alkylidene complexes, od combined with the natural orbitals of the chemical valence studied in this work, is singlet [18, 42]. According to the (ETS-NOCV) [19–22]. The Hoveyda-Grubbs catalysts have PBE0/def2-SVP calculations, the triplet and quintet states been additionally characterized by means of the interacting of the 16-electron complex 1 are less stable than the sin- quantum atoms (IQA) [23] energy decomposition. This is glet state, by 93 and 191 kJ mol−1, respectively. In the the first time when these complementary theoretical ap- case of the 14-electron species 1_dis, formed after disso- proaches (ETS-NOCVand IQA) are used to get deeper insight ciation of 1, the corresponding differences are 75 and into formation of active species from olefin metathesis 171 kJ mol−1. Hence, the spin-restricted formalism is used catalysts. for other computations. Furthermore, the DFT approach was proved sufficiently accurate in studying the reactivity of ruthenium alkylidene complexes [12, 13, 18, 28–33]. However, it shall be noticed, that for more difficult cases Computational methods (e.g., open-shell transition metal complexes), single- determinant DFT shall be applied with special caution Geometry optimization of all structures was carried out using due to the approximated character of electron correlation, the hybrid GGA PBE0 functional [24] combined with the as well as a possible multireference character of split-valence def2-SVP basis set [25]. The 28 innermost elec- wavefunction. trons of Ru were replaced by the Stuttgart effective core po- All the above calculations were done with the Gaussian 09 tential [26]. Further single point energy calculations were per- software [43]. formed with the hybrid meta-GGA M06 functional [27]and For the Grubbs type catalysts, the charge and energy de- the triple-zeta valence def2-TZVPP basis set [25, 26]. The composition ETS-NOCV [19, 20] analysis was performed. M06 method is recommended for main-group and transition DFT (BLYP [44, 45]) calculations were carried out with the metals thermochemistry, kinetics, and studies of noncovalent Amsterdam Density Functional (ADF) program [21, 22], in interactions [27]. A good performance of the M06 functional which the ETS-NOCV scheme was implemented. The pre- in predicting energies of reactions involving ruthenium sented energies include DFT-D3(BJ) dispersion corrections alkylidene complexes [28–33], including the dissociation en- [46, 47]. The standard triple-zeta STO basis set containing ergy for the Grubbs catalysts [28–31], was proved and it was one set of polarization functions (TZP) was used for the ru- often used for investigations of real ruthenium systems thenium atom, whereas the double-zeta basis set (DZP) was [28–36]. On the other hand, the PBE0 method was shown to employed for other elements as implemented in the ADF [21, be accurate in reproducing experimental geometry of the sec- 22]. The contours of deformation densities were plotted based ond generation Grubbs catalyst (1)[28, 37], being even better on the ADF-GUI