Computational Investigation of Scandium-Based Catalysts for Olefin

3 days ago - ... for developing new olefin hydroaminoalkylation reactions by using simple ... Transition-Metal-Free C–C, C–O, and C–N Cross-Coup...
5 downloads 0 Views 3MB Size
Article Cite This: Organometallics XXXX, XXX, XXX−XXX

pubs.acs.org/Organometallics

Computational Investigation of Scandium-Based Catalysts for Olefin Hydroaminoalkylation and C−H Addition Gen Luo,‡ Fan Liu,‡ Yi Luo,*,‡ Guangli Zhou,‡ Xiaohui Kang,‡,§ Zhaomin Hou,*,∥ and Lun Luo*,‡,†

Downloaded via IDAHO STATE UNIV on April 16, 2019 at 12:12:41 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.



Hubei Key Laboratory of Wudang Local Chinese Medicine Research, School of Pharmaceutical Sciences, Hubei University of Medicine, Shiyan 442000, China ‡ State Key Laboratory of Fine Chemicals, School of Chemical Engineering, Dalian University of Technology, Dalian 116024, China § College of Pharmacy, Dalian Medical University, Dalian, Liaoning 116044, China ∥ Organometallic Chemistry Laboratory, RIKEN Cluster for Pioneering Research, and Advanced Catalysis Research Group, RIKEN Center for Sustainable Resource Science, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan S Supporting Information *

ABSTRACT: Great progress has been achieved in the olefin hydroaminoalkylation by using amines, which is an atom-efficient route for the synthesis of alkylated amine derivatives. However, success in the catalytic olefin hydroaminoalkylation with a simple tertiary amine is hitherto very limited. In this study, density functional theory was applied to investigate the hydroaminoalkylation of olefins with tertiary amines, catalyzed by a series of homoleptic tris(benzyl) scandium complexes. It is found that the catalytic performance can be improved via substitution of electron-withdrawing groups and modifying ligand frameworks to reduce their steric hindrance. In addition, the potential applications of scandium catalysts in the α-C(sp3)−H alkylation of various heteroatomcontaining (P, As, O, S, and Se) substrates were explored. The results suggest that alkyl sulfides and selenides are promising substrates to undergo α-C(sp3)−H addition to olefins. Importantly, the effects of ligand backbone and substituent on catalytic performance and the different reactivities of the heteroatom-containing substrates were elucidated by frontier orbital, natural charge, topographic steric map, and distortion−interaction analyses, which give considerable insight into catalytic systems. This work provides useful information for developing new olefin hydroaminoalkylation reactions by using simple tertiary amines and for the addition of α-C(sp3)−H bond of heteroatom-containing substrates to alkenes.



INTRODUCTION Amine motifs are important functional groups, which are involved in a large number of agrochemicals, pharmaceuticals, and fine chemicals. Therefore, the development of effective approaches toward amine functionalization has attracted considerable interests in the past few decades.1 Of the possible approaches, the catalytic olefin hydroaminoalkylation by using amines is an atom-efficient method to synthesize alkylated amine derivatives.2−4 Significant advances have been made in early-transition-metal-catalyzed hydroaminoalkylation.4−6 These reactions are generally thought to initially go through the formation of metal amido intermediates via deprotonation of an N−H bond, which limits the amine substrates to primary and secondary amines.4 Although olefin hydroaminoalkylation by using tertiary amines can be accomplished through latetransition-metal catalysts,7,8 the reactions are often limited to the amines having a directing group such as pyridyl.7 The development of new systems for the catalytic olefin hydroaminoalkylation with simple tertiary amines is therefore still of great interest and significance. Hou et al. recently reported that a homoleptic tris(benzyl) scandium complex was an excellent © XXXX American Chemical Society

catalyst working for the addition of C−H bond of various aliphatic tertiary amines into olefins (Scheme 1).9 Remarkably, Scheme 1. Sc-Catalyzed Hydroaminoalkylation9

that work demonstrated the first catalytic olefin hydroaminoalkylation using simple aliphatic tertiary amines without a directing group.9,10 Subsequent theoretical work revealed that such a reaction generally follows the sequential events: active species formation, olefin insertion, and C−H activation of another molecule of amine (Scheme 2).11,12 Furthermore, electronic factors were found to be responsible for the regioselectivity experimentally observed.11 To date, however, Received: December 15, 2018

A

DOI: 10.1021/acs.organomet.8b00906 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

series of proposed homoleptic tris(benzyl) scandium catalysts for the hydroaminoalkylation of olefins to systematically compare their catalytic performances. The effects of the substituents and the ligand framework of these Sc-based catalysts were elucidated, providing useful information for catalyst design. In addition, the α-C(sp3)−H alkylation reactivities of various model substrates, that is, PMe3, AsMe3, OMe2, SMe2, and SeMe2, were also investigated. It is theoretically predicted that the alkyl sulfides and selenides are promising substrates for achieving the addition of α-C(sp3)−H to alkenes.

Scheme 2. Mechanistic Synopsis of Sc-Catalyzed Hydroaminoalkylation11



COMPUTATIONAL METHODS

The Gaussian 09 program15 was applied in all calculations. The geometry optimizations were carried out by using B3PW91 functional.16 In these optimizations, the 6-31G* basis set was used for nonmetallic atoms and the Stuttgart/Dresden relativistic effective core potentials17 as well as the associated valence basis sets were used for Sc atom. All structures were subsequently characterized, at the same level theory, by harmonic vibrational frequency analysis. Each stationary point was characterized as a minimum (no imaginary frequency) or transition state (one imaginary frequency). Single-point calculations were carried out at the M06/6-311+G** level with the CPCM model18 to consider the solvation effect of toluene. The free energy in solution is reported here, which includes the Gibbs free-energy correction calculated in the gas phase. Considering the overestimation of the entropy contribution, the free energies were corrected by +2.6 (or −2.6) kcal/mol for one-to-two (or two-to-one) molecular conversions, based on the free-volume theory.19 This empirical value is on the basis of energy of activation experimentally determined for the collision frequency Z′AB in liquids. Similar corrections were also applied in many earlier computational studies.20 The theoretical method used here was successfully applied to a similar system in our previous studies.11 The three-dimensional pictures of optimized structures were generated by using CYLview.21

factors governing the catalytic performance of the homoleptic tris(benzyl) scandium complex in hydroaminoalkylation, such as the effects of substituents and ligand frameworks with various heteroatoms, are not well understood. Because of the importance of this reaction, continuous efforts toward the development of more efficient catalysts are still highly desired in this field. Theoretical calculations are an effective approach for assessing various factors of a catalyst, such as the metal center, substituents, and ligand backbone, at the molecular level. This tool can facilitate the design of new catalysts with better catalytic performances.13 To identify more efficient catalytic systems for olefin hydroaminoalkylation, a deep theoretical understanding of the effects of the metal center, substituents, and ligand framework would be helpful. In addition, the exploration of αC(sp3)−H alkylation of various heteroatom (e.g., phosphorus, arsenic, oxygen, sulfur, and selenium)-containing substrates is also of much interest and significance. Similar to the addition of α-C(sp3)−H bond of amines, the direct addition of α-C(sp3)− H of other heteroatom-containing substrates into simple olefins is also challenging and still underexplored.14 In this work, density functional theory calculations were utilized to study a



RESULTS AND DISCUSSION Overview of the Catalytic Cycle. Regarding rare-earthcatalyzed hydroaminoalkylation, some analogous complexes

Scheme 3. Detailed Mechanism for the Sc-Catalyzed n-Hexene Hydroaminoalkylation with N-Methylpiperidine11

B

DOI: 10.1021/acs.organomet.8b00906 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

another molecule of N-methylpiperidine coordinates to the metal center to form complex G.23 Subsequently, G undergoes σ-BM via TS3 to give complex H with the coordinating alkylation product P. Finally, the catalytic cycle completed via ligand exchange between a1 and H23 to yield the product P and regenerate D. It is obvious that the σ-BM (C−H activation) is the rate-determining step with an energy barrier of 28.3 kcal/ mol.11 Hereafter, the exploration of catalyst designs and αC(sp3)−H alkylation of the other heteroatom-containing substrates is discussed based on this mechanism. Effects of the Substituents. To the base scandium catalyst Sc-1,9 electron-donating group (EDG) NH2 and electronwithdrawing group (EWG) NO2 substitute the para position of NMe2, respectively (Sc-NH2 and Sc-NO2, Figure 1). The calculated energy profiles of the hydroaminoalkylation of nhexene with N-methylpiperidine catalyzed by Sc-NH2 or ScNO2 are indicated in Figure 2. The previously reported energy profile of Sc-1-catalyzed reaction is also shown in this figure for a comparison.11 The coordination complex BSc‑NO2 with NO2 group is lower in energy than BSc‑1 by 4.1 kcal/mol, whereas the complex BSc‑NH2 with the NH2 EDG is higher in energy than BSc‑1 by 1.7 kcal/mol. For the active species generation, all three systems have moderate energy barriers of 25−27 kcal/mol, which is kinetically accessible under typical experimental conditions. It is noted that the active species DSc‑NO2 is significantly more stable than the corresponding species DSc‑1 and DSc‑NH2, and its formation is exergonic by 22.2 kcal/mol. In the catalytic cycle, the olefin insertion energy barriers for all three cases are similar (ΔG2⧧ = 21−23 kcal/mol). The subsequent C−H activation becomes the rate-determining step of catalytic cycle and has energy barriers of 28.3, 27.8, and 27.5 kcal/mol for Sc-1, Sc-NH2, and Sc-NO2, respectively. It should be noted that although the newly designed Sc-NH2 has a slightly lower energy barrier than Sc-1 (27.8 vs 28.3 kcal/mol for the rate-determining steps), the stationary points involved in the Sc-NH2 system are obviously higher in energy than the corresponding stationary points in the Sc-1 system. Therefore, modifying the Sc catalyst by introducing EDGs is thermodynamically unfavorable for this reaction. In contrast, the reaction

Figure 1. Sc-based catalysts with different substituents and their corresponding cationic species A showing selected bond distances (Å).

(e.g., Y, Lu, Gd, and Sm) were demonstrated experimentally to be inactive;9 thus, various metal centers are not considered in this work. Before discussing the effects of the substituents and ligand backbone, the features of the catalytic cycle are described.11 As shown in Schemes 1 and 2, the combination of [Ph3C][B(C6F5)4] and a neutral complex [Sc] can generate a cationic species A by loss of one alkyl group.22 In addition, the hydroaminoalkylation reaction includes the generation of the active species and the subsequent catalytic cycle (Scheme 3). In the step of active species generation, N-methylpiperidine (a1) coordinates to A, yielding the more stable complex B.23 Then, B transforms through the σ-bond metathesis (σ-BM)24,25 transition state (TS1) to give C with a coordinating N,Ndimethyl-o-toluidine part. Subsequently, C could undergo a ligand exchange process23 to generate the N-methylpiperidinecoordinated η2-azametallacyclic species D, which is the catalytically active species. In the catalytic cycle, D undergoes ligand exchange23 to yield π-complex E. Thereafter, the coordinating n-hexene in E inserts to the Sc−C bond through a conventional four-center transition state TS2, producing the ring-expanded azametallocyclic species F. Once olefin inserted,

Figure 2. Calculated energy profiles for the hydroaminoalkylation of n-hexene with N-methylpiperidine catalyzed by Sc-1, Sc-NH2, and Sc-NO2. The energies (kcal/mol) in the profiles include the energies of all species involved in the corresponding reactions. C

DOI: 10.1021/acs.organomet.8b00906 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

Figure 3. Frontier molecular orbital analyses of the cationic scandium species A (ASc‑1, ASc‑NH2, and ASc‑NO2) and N-methylpiperidine a1 (orbital energies in eV).

Figure 4. Sc-based catalysts with different ligands and their corresponding cationic species A showing selected bond distances (Å).

Figure 5. Energy profiles of the hydroaminoalkylation of n-hexene with N-methylpiperidine catalyzed by Sc-1, Sc-P, Sc-O, and Sc-S. The energies (kcal/mol) shown in the profiles include those of all species involved in the corresponding reactions.

catalyzed by Sc-NO2 is thermodynamically more favorable than that catalyzed by Sc-1, and the energy barrier of C−H activation in the former system is slightly lower (27.7 vs 28.3 kcal/mol, Figure 2). Therefore, introducing an EWG into the ligand could improve the catalytic activity13i and stabilize the reaction intermediates. The small differences in the bond lengths of Sc−N (2.18− 2.20 Å) and Sc−C (2.26−2.27 Å) in the three cationic species A

(Figure 1) suggest that the substituents NH2 and NO2 have little influence on the steric environment around the metal center (see also the topographic steric map analysis26,27 in Figure S1), and thus, steric factors should not have an obvious effect on the performance of the proposed catalysts. Therefore, the electronic effects of the substituents were investigated. Accordingly, a natural charge analysis was performed. It is found that the EWG increases the Lewis acidity of the Sc center, whereas the EDG D

DOI: 10.1021/acs.organomet.8b00906 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

Figure 6. Topographic steric maps of the cationic scandium species A with different heteroatom-containing frameworks. % VBur denotes the percent of buried volume of Sc atom.

Scheme 4. α-C(sp3)−H Alkylation of Various HeteroatomContaining Substrates

In addition, the frontier orbitals of catalysts and Nmethylpiperidine (a1) were analyzed. During the reaction, the LUMO (lowest unoccupied molecular orbital) of the electrophilic Sc catalyst could interact with the HOMO (the highest occupied molecular orbital) of the nucleophilic amine substrate. As illustrated in Figure 3, the LUMO of A mostly consists of the Sc 3dZ2 orbital and the HOMO of N-methylpiperidine is dominated by a sp-hybrided lone pair orbital of the N atom. Compared to ASc‑1, NH2 substitution clearly increases the LUMO energy from −3.14 to −2.85 eV, whereas NO2 substitution reduces the LUMO energy from −3.14 to −4.19 eV. Therefore, introducing an electron-withdrawing substituent (e.g., NO2) could effectively reduce the energy gap of E(LUMO,A) − E(HOMO,a1). In short, the electron-withdrawing substituent increases the Lewis acidity of Sc center and decreases the energy gap of E(LUMO,A) − E(HOMO,a1) (Figure 3), resulting in tighter binding between the catalyst and the substrate/product and thus stabilizing the reaction intermediates in this catalytic system. To verify this finding, the binding energies (ΔE) of the catalyst−substrate (complex D) and catalyst−product (complex H) complexes involved in both the Sc-1 and Sc-NO2 systems were calculated. As expected, the results clearly indicate that the binding between the catalyst and the substrate/product is stronger in the Sc-NO2 system than that in the Sc-1 system. Furthermore, the binding (in H) between the catalyst and product is stronger than that (in D) between the catalyst and amine substrate (ΔED(Sc‑1) = −48.5 kcal/mol; ΔED(Sc‑NO2) = −50.0 kcal/mol; ΔEH(Sc‑1) = −50.1 kcal/mol; ΔEH(Sc‑NO2) = −52.4 kcal/mol). Effects of the Heteroatom-Containing Ligand Framework. In addition to the substituents, the heteroatom of the ligand framework probably affects the catalytic performance as well. In this study, the N atom of the ligand framework was replaced with several common heteroatoms (O, S, and P). As shown in Figure 4, all of the cationic scandium analogues were successfully optimized and have similar structural characteristics. The energy profiles of the hydroaminoalkylation of hexene with N-methylpiperidine catalyzed by Sc[CH2C6H4XMen-o]3 [X = P (Sc-P); X = O (Sc-O); and X = S (Sc-S)] are shown in Figure 5, which also includes the results for the reaction catalyzed by Sc-1 for comparison. It is found that during active

Figure 7. Energy profiles for the generation of catalytically active species via various heteroatom-containing substrates XMe3 (X = P and As) and XMe2 (X = O, S, and Se). The energies (kcal/mol) shown in the profiles include those of all species involved in the corresponding reactions.

decreases the Lewis acidity, as revealed by the calculated natural charges on the Sc center of 1.70, 1.67, and 1.74 in ASc‑1, ASc‑NH2, and ASc‑NO2, respectively. This result explains the relative stability of the corresponding stationary points in the three systems (Figure 2); the more positive metal center binds the ligand more strongly, resulting in its higher stability. E

DOI: 10.1021/acs.organomet.8b00906 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

Figure 8. Optimized geometrical structures of the coordination complexes and transition states in the C(sp3)−H activation of heteroatom-containing substrates by the cationic scandium species ASc‑1 and the distortion−interaction analysis for corresponding transition states (energies in kcal/mol).

Figure 9. Topographic steric maps of the TSs for α-C(sp3)−H activation of heteroatom-containing substrates. The % VBur denotes the percent of buried volume of the heteroatom.

for improving the catalytic performance of this type of hydroaminoalkylation reaction. To obtain better understanding of the higher thermodynamic stability of the stationary points in the Sc-P/Sc-O/Sc-S systems compared to those in the Sc-1 system, the electronic effects of these catalysts were analyzed. The natural charges on the Sc atom of A suggest that only the O atom substitution (charge of 1.80 on Sc of ASc‑O) increases the Lewis acidity of the Sc center, and the P and S ligands (charges of 1.38 and 1.43 on Sc of ASc‑P and ASc‑S, respectively) decrease the Lewis acidity compared to the catalyst with N ligands (1.70 on Sc of ASc‑1). This result does not explain the relative thermodynamic stabilities of the stationary points indicated in Figure 5. Besides, the frontier orbital analyses of these catalysts do not provide a rational explanation as well (see Figure S2). These results indicate that

species generation, the coordination complex B and the active species D in the Sc-P/Sc-O/Sc-S systems are significantly more stable than those in the Sc-1 system. Although the energy barriers of TS1 (ΔG⧧1 ) in the newly designed Sc-P/Sc-O/Sc-S catalyst systems are slightly higher than that in the Sc-1 system, the reaction should be kinetically accessible under experimental conditions in all of the systems. The energy barriers of olefin insertion and C−H activation in the three newly designed systems are moderate (less than 28 kcal/mol).28 It is noteworthy that all of the stationary points in the three systems (Sc-P, Sc-O, and Sc-S) are more stable than that in the Sc-1 system. The energy barriers are energetically accessible, especially for the ScS system having the lowest energy barrier (ca. 26 kcal/mol). On the basis of these comparisons, the S-containing complex is the most promising candidate of the three newly designed catalysts F

DOI: 10.1021/acs.organomet.8b00906 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

performed.29 In Figure 8, ΔEdist(cat) is the energy required to distort the cationic scandium species structure in coordination complex B to that in transition state TS1, whereas ΔEdist(XMen) is the energy needed to distort the coordinated XMen in coordination complex B to that in TS1. The difference in the interaction energies between the distorted Sc active species and XMen in complex B and transition state TS1 is denoted ΔEint. The activation energy ΔE⧧ is the sum of ΔEdist(cat), ΔE dist (XMe n ), and ΔE int . As shown in Figure 8, the distortion−interaction analysis reveals that the activation energies (ΔE⧧) of the five TS1 transition states are consistent with the trends in their corresponding free-energy barriers (ΔG⧧). Small variations in ΔEdist(cat) (+7.8 to +8.9 kcal/mol) and ΔEint (−12.2 to −15.4 kcal/mol) are clearly observed between all five systems. However, the values of ΔEdist(XMen) for the SMe2 (+31.7 kcal/mol) and SeMe2 (+31.9 kcal/mol) systems are significantly smaller than those for the other three systems (+38.2, +36.0, and +37.1 kcal/mol for PMe3, OMe2, and AsMe3, respectively). The differences between ΔEdist(XMen) of the SeMe2 system and those of the other three systems (PMe3, OMe2, and AsMe3) are nearly equivalent to the differences between related activation energies (ΔE⧧). Therefore, the differences in distortion energy of the XMen moiety in transition states are the origin of energy barrier differences of these systems. The relatively small values of ΔEdist(XMen) for SMe2 and SeMe2 suggest that the steric effects in TS1(SMe2) and TS1(SeMe2) are smaller than those in the other transition states. To confirm these steric effects, the topographic steric map analysis of the five TS1 was carried out (Figure 9).26,27 As expected, the percent of buried volume of the heteroatom center in TS1(PMe3) (% VBur = 75.3), TS1(AsMe3) (% VBur = 73.8), and TS1(OMe2) (% VBur = 71.7) is significantly larger than those in TS1(SMe2) (% VBur = 62.3) and TS1(SeMe2) (% VBur = 59.2), which is in line with the results of distortion−interaction analyses. The results of the distortion−interaction and topographic steric map analyses can elucidate the origin of the energy barrier differences of the reactions with various substrates. To further explore the feasibility of the catalytic cycle for the SMe2 and SeMe2 substrates, the reaction pathways of the C−H activation and olefin insertion were also calculated (Figure 10). It is found that the energy barriers of the olefin insertion are 21.4 and 23.1 kcal/mol for SMe2 and SeMe2, respectively. In addition, the energy barriers of corresponding C−H activation are 22.2 and 24.3 kcal/mol, respectively. Such moderate energy barriers are kinetically accessible under typical experimental conditions. In this sense, alkyl sulfides and selenides are promising candidates for α-C(sp3)−H alkylation.

Figure 10. Energy profiles of XMe2 (X = S and Se) in the catalytic cycle. The free energies are relative to A and the substrates and include the energies of all species involved in the corresponding reactions (energy in kcal/mol).

electronic factors might not be the origin of the differences in the catalytic performances of these catalysts. Because the differences in the sizes of heteroatoms (N, O, P, and S) and the number of methyl groups of the ligands might result in differences in the steric hindrance of these catalysts, the topographic steric maps of the four cationic scandium species (ASc‑1, ASc‑O, ASc‑P, and ASc‑S) were analyzed (Figure 6).26,27 The percent of buried volume (% VBur) of the Sc atom demonstrates the differences in the steric bulkiness of these four ligands; the N ligand in ASc‑1 (% VBur = 81.5) is clearly more bulky than the heteroatom-containing ligands in ASc‑P (% VBur = 74.3), ASc‑O (% VBur = 72.2), and ASc‑S (% VBur = 76.3). Thus, the stationary points involved in the former system are less thermodynamically favorable than those in the latter three systems. This steric effect correlates well with the differences in the thermodynamics of the stationary points in the systems shown in Figure 5. The results obtained indicate that the catalytic performance could be improved through modifying the catalyst with less sterically bulky ligands. Reactivities of the Substrates with Various Heteroatoms. Similarly, to hydroaminoalkylation, the direct αC(sp3)−H alkylation of various heteroatom-containing (e.g., P, As, O, S, and Se) substrates with simple olefins is also a highly desirable but underexplored methodology (Scheme 4).14 To computationally explore the application scope of the Sc-based catalyst, α-C(sp3)−H alkylation of model substrates XMe3 (X = P and As) and XMe2 (X = O, S, and Se) with n-hexene was investigated (Scheme 4 and Figure 7). As shown in Figure 7, the generation of the active species through the reaction of ASc‑1 with PMe3 or AsMe3 has energy barrier of up to 33 kcal/mol. Besides, the formation of active species D is endothermic, which is 15.4 and 20.0 kcal/mol higher in energy compared with B in PMe3 and AsMe3 systems, respectively. Therefore, the C−H activation of PMe3 or AsMe3 by ASc‑1 is difficult to achieve under typical experimental conditions. Regarding the group VIA substrates (XMe2, X = O, S, and Se), the reaction for the OMe2 substrate also shows high-energy barrier (32.4 kcal/mol), whereas moderate energy barriers are observed for SeMe2 (28.4 kcal/mol) and SMe2 (28.2 kcal/mol), suggesting that αC(sp3)−H bonds in SeMe2 and SMe2 can be activated by ASc‑1. Therefore, the alkyl sulfides and selenides are promising substrates for generating active species via α-C(sp3)−H activation under appropriate conditions. To explore the origin of the differences in free-energy barriers of transition states, distortion−interaction analyses were



CONCLUSIONS In summary, a series of homoleptic tris(benzyl) scandium complexes were computationally investigated as potential catalysts for olefin hydroaminoalkylation with simple tertiary amines to determine the effects of the substituents and heteroatoms. The key discoveries and insights are as follows. (1) The catalytic performance of the Sc-based catalyst could be improved by the introduction of EWGs because they increase the Lewis acidity of the metal atom and reduce the energy gap between the LUMO of the catalyst and the HOMO of the amine substrate. (2) Modifying the ligand framework by heteroatom substitution is also a promising approach for improving the catalytic performance. By these means, the steric hindrance G

DOI: 10.1021/acs.organomet.8b00906 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

transition-metal dialkylamides and alkoxides. Organometallics 1983, 2, 161−162. (4) For reviews, see (a) Roesky, P. W. Catalytic hydroaminoalkylation. Angew. Chem., Int. Ed. 2009, 48, 4892−4894. (b) Chong, E.; Garcia, P.; Schafer, L. L. Hydroaminoalkylation: Early-transition-metal-catalyzed α-alkylation of amines. Synthesis 2014, 46, 2884−2896. (c) Ryken, S. A.; Schafer, L. L. N,O-Chelating Four-Membered Metallacyclic Titanium(IV) Complexes for Atom-Economic Catalytic Reactions. Acc. Chem. Res. 2015, 48, 2576−2586. (d) Edwards, P. M.; Schafer, L. L. Early transition metal-catalyzed C−H alkylation: hydroaminoalkylation for Csp3−Csp3 bond formation in the synthesis of selectively substituted amines. Chem. Commun. 2018, 54, 12543−12560. (e) Hannedouche, J.; Schulz, E. Hydroamination and hydroaminoalkylation of alkenes by group 3−5 elements: Recent developments and comparison with late transition metals. Organometallics 2018, 37, 4313−4326. (5) Examples of group 4 metal catalysts, see: (a) Bielefeld, J.; Doye, S. Dimethylamine as a substrate in hydroaminoalkylation reactions. Angew. Chem., Int. Ed. 2017, 56, 15155−15158. (b) Manßen, M.; Lauterbach, N.; Dörfler, J.; Schmidtmann, M.; Saak, W.; Doye, S.; Beckhaus, R. Efficient access to titanaaziridines by C−H activation of N-methylanilines at ambient temperature. Angew. Chem., Int. Ed. 2015, 54, 4383−4387. (c) Dörfler, J.; Preuß, T.; Schischko, A.; Schmidtmann, M.; Doye, S. A 2,6-Bis(phenylamino)pyridinato titanium catalyst for the highly regioselective hydroaminoalkylation of styrenes and 1,3butadienes. Angew. Chem., Int. Ed. 2014, 53, 7918−7922. (d) Chong, E.; Schafer, L. L. 2-Pyridonate titanium complexes for chemoselectivity. Accessing intramolecular hydroaminoalkylation over hydroamination. Org. Lett. 2013, 15, 6002−6005. (e) Bexrud, J. A.; Eisenberger, P.; Leitch, D. C.; Payne, P. R.; Schafer, L. L. Selective C−H activation α to primary amines. Bridging metallaaziridines for catalytic, intramolecular α-alkylation. J. Am. Chem. Soc. 2009, 131, 2116−2118. (f) Preuß, T.; Saak, W.; Doye, S. Titanium-catalyzed intermolecular hydroaminoalkylation of conjugated dienes. Chem.Eur. J. 2013, 19, 3833−3837. (g) Dörfler, J.; Doye, S. Aminopyridinato titanium catalysts for the hydroaminoalkylation of alkenes and styrenes. Angew. Chem., Int. Ed. 2013, 52, 1806−1809. (h) Prochnow, I.; Zark, P.; Müller, T.; Doye, S. The mechanism of the titanium-catalyzed hydroaminoalkylation of alkenes. Angew. Chem., Int. Ed. 2011, 50, 6401−6405. (i) Jaspers, D.; Saak, W.; Doye, S. Dinuclear titanium complexes with sulfamide ligands as precatalysts for hydroaminoalkylation and hydroamination reactions. Synlett 2012, 23, 2098−2102. (j) Kubiak, R.; Prochnow, I.; Doye, S. Titanium-catalyzed hydroaminoalkylation of alkenes by C−H bond activation at sp3 centers in the α-position to a nitrogen atom. Angew. Chem., Int. Ed. 2009, 48, 1153−1156. (k) Kubiak, R.; Prochnow, I.; Doye, S. [Ind2TiMe2]: A catalyst for the hydroaminomethylation of alkenes and styrenes. Angew. Chem., Int. Ed. 2010, 49, 2626−2629. (l) Stelter, L.; Teusch, T.; Bielefeld, J.; Doye, S.; Klüner, T. Theoretical studies on the hydroaminoalkylation of alkenes with primary and secondary amines. Chem.Eur. J. 2018, 24, 12485−12489. (6) Examples of group 5 metal catalysts, see: (a) DiPucchio, R. C.; Roşca, S.-C.; Schafer, L. L. Catalytic and atom-economic Csp3-Csp3 bond formation: Alkyl tantalum ureates for hydroaminoalkylation. Angew. Chem., Int. Ed. 2018, 57, 3469−3472. (b) Chong, E.; Brandt, J. W.; Schafer, L. L. 2-Pyridonate tantalum complexes for the intermolecular hydroaminoalkylation of sterically demanding alkenes. J. Am. Chem. Soc. 2014, 136, 10898−10901. (c) Garcia, P.; Lau, Y. Y.; Perry, M. R.; Schafer, L. L. Phosphoramidate tantalum complexes for room-temperature C−H functionalization: hydroaminoalkylation catalysis. Angew. Chem., Int. Ed. 2013, 52, 9144−9148. (d) Reznichenko, A. L.; Hultzsch, K. C. The mechanism of hydroaminoalkylation catalyzed by group 5 metal binaphtholate complexes. J. Am. Chem. Soc. 2012, 134, 3300−3311. (e) Eisenberger, P.; Ayinla, R. O.; Lauzon, J. M. P.; Schafer, L. L. Tantalum-amidate complexes for the hydroaminoalkylation of secondary amines: Enhanced substrate scope and enantioselective chiral amine synthesis. Angew. Chem., Int. Ed. 2009, 48, 8361−8365. (f) Herzon, S. B.; Hartwig, J. F. Hydroaminoalkylation of unactivated olefins with dialkylamines. J. Am. Chem. Soc. 2008, 130, 14940−14941. (g) Herzon, S. B.; Hartwig, J. F. Direct, catalytic hydroaminoalkylation of unactivated olefins with N-alkyl arylamines. J.

surrounding the metal atom could be effectively decreased. Of the N-, P-, O-, and S-containing ligand frameworks, the Scontaining catalyst exhibits the best catalytic performance. (3) To explore the application scope of the catalysts, the α-C(sp3)− H addition of various heteroatom-containing substrates to olefins was calculated. It is found that P-, As-, and O-containing substrates are inert because of steric factors in the transition state of C−H activation. In contrast, alkyl sulfides and selenides are calculated to be promising candidates for α-C(sp3)−H alkylation. These results provide valuable information on catalyst design and the application of these catalysts in the alkylation of α-C(sp3)−H of heteroatom-containing compounds.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.organomet.8b00906. Topographic steric maps and frontier orbital analysis (PDF) Cartesian coordinates of all optimized structures and their single-point energies (XYZ)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected] (Y.L.). *E-mail: [email protected] (Z.H.). *E-mail: [email protected] (L.L.). ORCID

Gen Luo: 0000-0002-5297-6756 Yi Luo: 0000-0001-6390-8639 Xiaohui Kang: 0000-0003-2793-7111 Zhaomin Hou: 0000-0003-2841-5120 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS Financial supports from the NSFC (nos. 21429201, 21674014, and 21704011), the Fundamental Research Funds for the Central Universities (DUT18GJ201 for Y.L. and DUT18RC(3) 002 for G.L.), the State Key Laboratory of Fine Chemicals for a research grant (KF1713 for X.K.), and a Faculty Development Grant from Hubei University of Medicine (2018QDJZR13 for L.L.) are gratefully acknowledged. We also thank the RIKEN’s HOKUSAI system and the Network and Information Center of Dalian University of Technology for computational support.



REFERENCES

(1) (a) Salvatore, R. N.; Yoon, C. H.; Jung, K. W. Synthesis of secondary amines. Tetrahedron 2001, 57, 7785−7811. (b) Hager, A.; Vrielink, N.; Hager, D.; Lefranc, J.; Trauner, D. Synthetic approaches towards alkaloids bearing α-tertiary amines. Nat. Prod. Rep. 2016, 33, 491−522. (c) Li, B.; Sortais, J.-B.; Darcel, C. Amine synthesis via transition metal homogeneous catalysed hydrosilylation. RSC Adv. 2016, 6, 57603−57625. (2) Dong, Z.; Ren, Z.; Thompson, S. J.; Xu, Y.; Dong, G. Transitionmetal-catalyzed C−H alkylation using alkenes. Chem. Rev. 2017, 117, 9333−9403. (3) (a) Clerici, M. G.; Maspero, F. Catalytic C-alkylation of secondary amines with alkenes. Synthesis 1980, 305−306. (b) Nugent, W. A.; Ovenall, D. W.; Holmes, S. J. Catalytic C−H activation in early H

DOI: 10.1021/acs.organomet.8b00906 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics Am. Chem. Soc. 2007, 129, 6690−6691. (h) Reznichenko, A. L.; Emge, T. J.; Audörsch, S.; Klauber, E. G.; Hultzsch, K. C.; Schmidt, B. Group 5 metal binaphtholate complexes for catalytic asymmetric hydroaminoalkylation and hydroamination/cyclization. Organometallics 2011, 30, 921−924. (i) Lauzon, J. M.; Eisenberger, P.; Roşca, S.-C.; Schafer, L. L. Amidate complexes of tantalum and niobium for the hydroaminoalkylation of unactivated alkenes. ACS Catal. 2017, 7, 5921−5931. (7) Examples of late-transition metal catalysts, see: (a) Chatani, N.; Asaumi, T.; Yorimitsu, S.; Ikeda, T.; Kakiuchi, F.; Murai, S. Ru3(CO)12catalyzed coupling reaction of sp3 C−H bonds adjacent to a nitrogen atom in alkylamines with alkenes. J. Am. Chem. Soc. 2001, 123, 10935− 10941. (b) Bergman, S. D.; Storr, T. E.; Prokopcová, H.; Aelvoet, K.; Diels, G.; Meerpoel, L.; Maes, B. U. W. The role of the alcohol and carboxylic acid in directed ruthenium-catalyzed C(sp3)−H α-alkylation of cyclic amines. Chem.Eur. J. 2012, 18, 10393−10398. (c) Schinkel, M.; Wang, L.; Bielefeld, K.; Ackermann, L. Ruthenium(II)-catalyzed C(sp3)−H α-alkylation of pyrrolidines. Org. Lett. 2014, 16, 1876− 1879. (8) (a) Thullen, S. M.; Rovis, T. A mild hydroaminoalkylation of conjugated dienes using a unified cobalt and photoredox catalytic system. J. Am. Chem. Soc. 2017, 139, 15504−15508. (b) Perez, F.; Oda, S.; Geary, L. M.; Krische, M. J. Ruthenium-catalyzed transfer hydrogenation for C−C bond formation: Hydrohydroxyalkylation and hydroaminoalkylation via reactant redox pairs. Top. Curr. Chem. 2016, 374, 35. (c) Chen, T.-Y.; Tsutsumi, R.; Montgomery, T. P.; Volchkov, I.; Krische, M. J. Ruthenium-catalyzed C−C coupling of amino alcohols with dienes via transfer hydrogenation: Redox-triggered imine addition and related hydroaminoalkylations. J. Am. Chem. Soc. 2015, 137, 1798−1801. (d) Schmitt, D. C.; Lee, J.; Dechert-Schmitt, A.-M. R.; Yamaguchi, E.; Krische, M. J. Ruthenium catalyzed hydroaminoalkylation of isoprene via transfer hydrogenation: Byproduct-free prenylation of hydantoins. Chem. Commun. 2013, 49, 6096−6098. (e) Trowbridge, A.; Reich, D.; Gaunt, M. J. Multicomponent synthesis of tertiary alkylamines by photocatalytic olefinhydroaminoalkylation. Nature 2018, 561, 522−527. (9) Nako, A. E.; Oyamada, J.; Nishiura, M.; Hou, Z. Scandiumcatalysed intermolecular hydroaminoalkylation of olefins with aliphatic tertiary amines. Chem. Sci. 2016, 7, 6429−6434. (10) Gao, H.; Su, J.; Xu, P.; Xu, X. Scandium-catalyzed C(sp3)−H alkylation of N,N-dimethyl anilines with alkenes. Org. Chem. Front. 2018, 5, 59−63. (11) Liu, F.; Luo, G.; Hou, Z.; Luo, Y. Mechanistic insights into scandium-catalyzed hydroaminoalkylation of olefins with amines: Origin of regioselectivity and charge-based prediction model. Organometallics 2017, 36, 1557−1565. (12) (a) Luo, G.; Luo, Y.; Qu, J.; Hou, Z. Mechanistic investigation on scandium-catalyzed C−H addition of pyridines to olefins. Organometallics 2012, 31, 3930−3937. (b) Song, G.; Luo, G.; Oyamada, J.; Luo, Y.; Hou, Z. Ortho-selective C−H addition of N,N-dimethyl anilines to alkenes by a yttrium catalyst. Chem. Sci. 2016, 7, 5265−5270. (c) Zhou, G.; Luo, G.; Kang, X.; Hou, Z.; Luo, Y. Origin of product selectivity in yttrium-catalyzed benzylic C−H alkylations of alkyl pyridines with olefins: A DFT study. Organometallics 2018, 37, 2741− 2748. (13) (a) Wang, G.; Zhang, H.; Zhao, J.; Li, W.; Cao, J.; Zhu, C.; Li, S. Homolytic cleavage of a B−B bond by the cooperative catalysis of two Lewis bases: Computational design and experimental verification. Angew. Chem., Int. Ed. 2016, 55, 5985−5989. (b) Ge, H.; Chen, X.; Yang, X. Hydrogenation of carbon dioxide to methanol catalyzed by iron, cobalt, and manganese cyclopentadienone complexes: Mechanistic insights and computational design. Chem.Eur. J. 2017, 23, 8850−8856. (c) Ge, H.; Chen, X.; Yang, X. A mechanistic study and computational prediction of iron, cobalt and manganese cyclopentadienone complexes for hydrogenation of carbon dioxide. Chem. Commun. 2016, 52, 12422−12425. (d) Chen, X.; Yang, X. Bioinspired design and computational prediction of iron complexes with pendant amines for the production of methanol from CO2 and H2. J. Phys. Chem. Lett. 2016, 7, 1035−1041. (e) Zeng, G.; Maeda, S.; Taketsugu, T.;

Sakaki, S. Catalytic hydrogenation of carbon dioxide with ammoniaborane by pincer-type phosphorus compounds: Theoretical prediction. J. Am. Chem. Soc. 2016, 138, 13481−13484. (f) Zeng, G.; Maeda, S.; Taketsugu, T.; Sakaki, S. Theoretical study of hydrogenation catalysis of phosphorus compound and prediction of catalyst with high activity and wide application scope. ACS Catal. 2016, 6, 4859−4870. (g) Zhu, B.; Yan, L.-K.; Yao, L.-S.; Ren, H.; Li, R.-H.; Guan, W.; Su, Z.-M. Orthogonal reactivity of Ni(I)/Pd(0) dual catalysts for Ullmann C−C cross-coupling: Theoretical insight. Chem. Commun. 2018, 54, 7959− 7962. (h) Zhu, B.; Yan, L.-K.; Geng, Y.; Ren, H.; Guan, W.; Su, Z.-M. IrIII/NiII-metallaphotoredox catalysis: the oxidation state modulation mechanism versus the radical mechanism. Chem. Commun. 2018, 54, 5968−5971. (i) Gilmour, D. J.; Lauzon, J. M. P.; Clot, E.; Schafer, L. L. Ta-catalyzed hydroaminoalkylation of alkenes: Insights into ligandmodified reactivity using DFT. Organometallics 2018, 37, 4387−4394. (14) (a) Chen, J.; Kirchmeier, R. L.; Shreeve, J. n. M. Insertion of fluoroalkenes into activated C−H bonds for the preparation of polyfluorinated sulfanes, alcohols, and acyclic and cyclic ethers. Inorg. Chem. 1996, 35, 6676−6681. (b) Zhang, S.-Y.; Zhang, F.-M.; Tu, Y.-Q. Direct sp3α-C−H activation and functionalization of alcohol and ether. Chem. Soc. Rev. 2011, 40, 1937−1949. (c) Cao, H.; Liu, D.; Liu, C.; Hu, X.; Lei, A. Copper-catalyzed oxidative alkenylation of thioethers via Csp3−H functionalization. Org. Biomol. Chem. 2015, 13, 2264−2266. (d) Luo, Y.; Ma, Y.; Hou, Z. α-C−H alkylation of methyl sulfides with alkenes by a scandium catalyst. J. Am. Chem. Soc. 2018, 140, 114−117. (15) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, N. J.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö .; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09, Revision A.01; Gaussian, Inc.: Wallingford, CT, 2009. (16) (a) Becke, A. D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648−5652. (b) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti correlationenergy formula into a functional of the electron density. Phys. Rev. B: Condens. Matter Mater. Phys. 1988, 37, 785−789. (c) Perdew, J. P.; Burke, K.; Wang, Y. Generalized gradient approximation for the exchange-correlation hole of a many-electron system. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54, 16533−16539. (17) (a) Dolg, M.; Wedig, U.; Stoll, H.; Preuss, H. Energy-adjusted ab initio pseudopotentials for the first row transition elements. J. Chem. Phys. 1987, 86, 866−872. (b) Schwerdtfeger, P.; Dolg, M.; Schwarz, W. H. E.; Bowmaker, G. A.; Boyd, P. D. W. Relativistic effects in gold chemistry. I. Diatomic gold compounds. J. Chem. Phys. 1989, 91, 1762− 1774. (c) Dolg, M.; Stoll, H.; Savin, A.; Preuss, H. Energy-adjusted pseudopotentials for the rare earth elements. Theor. Chim. Acta 1989, 75, 173−194. (18) (a) Barone, V.; Cossi, M. Quantum calculation of molecular energies and energy gradients in solution by a conductor solvent model. J. Phys. Chem. A 1998, 102, 1995−2001. (b) Cossi, M.; Rega, N.; Scalmani, G.; Barone, V. Energies, structures, and electronic properties of molecules in solution with the C-PCM solvation model. J. Comput. Chem. 2003, 24, 669−681. (19) Benson, S. W. The Foundations of Chemical Kinetics; Krieger: Malabar, FL, 1982. (20) (a) Bai, W.; Wei, W.; Sung, H. H. Y.; Williams, I. D.; Lin, Z.; Jia, G. Syntheses of Re(V) alkylidyne complexes and ligand effect on the reactivity of Re(V) alkylidyne complexes toward alkynes. Organometallics 2018, 37, 559−569. (b) Yu, H.; Lu, Q.; Dang, Z.; Fu, Y. I

DOI: 10.1021/acs.organomet.8b00906 Organometallics XXXX, XXX, XXX−XXX

Article

Organometallics

(28) Considering the experimental condition (70 °C, 24−48 h) and the moderate yields (57−93%, see ref 9), the calculated energy barrier of ca. 26−28 kcal/mol could be reasonable. (29) For reviews, see: (a) Bickelhaupt, F. M.; Houk, K. N. Analyzing reaction rates with the distortion/interaction-activation strain model. Angew. Chem., Int. Ed. 2017, 56, 10070−10086. (b) Liu, F.; Liang, Y.; Houk, K. N. Bioorthogonal cycloadditions: Computational analysis with the distortion/interaction model and predictions of reactivities. Acc. Chem. Res. 2017, 50, 2297−2308. (c) Fernández, I.; Bickelhaupt, F. M. The activation strain model and molecular orbital theory: Understanding and designing chemical reactions. Chem. Soc. Rev. 2014, 43, 4953−4967.

Mechanistic study of the rhodium-catalyzed [3+2+2] carbocyclization of alkenylidenecyclopropanes with alkynes. Chem.Asian J. 2013, 8, 2262−2273. (c) Schoenebeck, F.; Houk, K. N. Ligand-controlled regioselectivity in palladium-catalyzed cross coupling reactions. J. Am. Chem. Soc. 2010, 132, 2496−2497. (d) Ardura, D.; López, R.; Sordo, T. L. Relative Gibbs energies in solution through continuum models: Effect of the loss of translational degrees of freedom in bimolecular reactions on Gibbs energy barriers. J. Phys. Chem. B 2005, 109, 23618− 23623. (e) Liu, Q.; Lan, Y.; Liu, J.; Li, G.; Wu, Y.-D.; Lei, A. Revealing a second transmetalation step in the Negishi coupling and its competition with reductive elimination: Improvement in the interpretation of the mechanism of biaryl syntheses. J. Am. Chem. Soc. 2009, 131, 10201− 10210. (f) Wang, B.; Luo, G.; Nishiura, M.; Hu, S.; Shima, T.; Luo, Y.; Hou, Z. Dinitrogen activation by dihydrogen and a PNP-ligated titanium complex. J. Am. Chem. Soc. 2017, 139, 1818−1821. (g) Kang, X.; Luo, G.; Luo, L.; Hu, S.; Luo, Y.; Hou, Z. Mechanistic insights into ring cleavage and contraction of benzene over a titanium hydride cluster. J. Am. Chem. Soc. 2016, 138, 11550−11559. (h) Bai, W.; Lee, K.-H.; Sung, H. H. Y.; Williams, I. D.; Lin, Z.; Jia, G. Alkyne metathesis reactions of rhenium(V) carbyne complexes. Organometallics 2016, 35, 3808−3815. (i) Xie, H.; Lin, Z. Understanding the reactivity difference of isocyanate and isothiocyanate toward a ruthenium silylene hydride complex. Organometallics 2014, 33, 892−897. (j) Yuan, R.; Lin, Z. Computational insight into the mechanism of nickel-catalyzed reductive carboxylation of styrenes using CO2. Organometallics 2014, 33, 7147−7156. (21) Legault, C. Y. CYLview, 1.0b; Université de Sherbrooke, 2009. http://www.cylview.org. (22) (a) Nishiura, M.; Guo, F.; Hou, Z. Half-sandwich rare-earthcatalyzed olefin polymerization, carbometalation, and hydroarylation. Acc. Chem. Res. 2015, 48, 2209−2220. (b) Nishiura, M.; Hou, Z. Novel polymerization catalysts and hydride clusters from rare-earth metal dialkyls. Nat. Chem. 2010, 2, 257−268. (23) (a) The transition states for the processes of ligand association (coordination), dissociation, and exchange were located for Sc-1 case (see Figures S3 and S4). It is found that these transition states have small imaginary frequencies (i30−i60 cm−1) and relatively low energy barriers, which could not affect the rate-determining barriers. Therefore, such transition states were not further considered for the remaining cases. (b) A method for estimation of coordination barrier, see: Hartwig, J. F.; Cook, K. S.; Hapke, M.; Incarvito, C. D.; Fan, Y.; Webster, C. E.; Hall, M. B. Rhodium boryl complexes in the catalytic, terminal functionalization of alkanes. J. Am. Chem. Soc. 2005, 127, 2538−2552. (24) Selected reviews, see: (a) Lin, Z. Current understanding of the C−H activation reactions of LnMR+R′−H→LnMR′ + R−H. Coord. Chem. Rev. 2007, 251, 2280−2291. (b) Balcells, D.; Clot, E.; Eisenstein, O. C−H bond activation in transition metal species from a computational perspective. Chem. Rev. 2010, 110, 749−823. (c) Kefalidis, C. E.; Castro, L.; Perrin, L.; Rosal, I. D.; Maron, L. New perspectives in organolanthanide chemistry from redox to bond metathesis: Insights from theory. Chem. Soc. Rev. 2016, 45, 2516−2543. (25) Examples of multimetal cooperating σ-bond metathesis E−H (E = C, N, P) activations, see: (a) Luo, G.; Luo, Y.; Zhang, W.; Qu, J.; Hou, Z. DFT studies on the methane elimination reaction of a trinuclear rareearth polymethyl complex: σ-Bond metathesis assisted by cooperation of multimetal sites. Organometallics 2014, 33, 1126−1134. (b) Luo, G.; Luo, Y.; Hou, Z. E−H (E = N and P) bond activation of PhEH2 by a trinuclear yttrium methylidene complex: Theoretical insights into mechanism and multimetal cooperation behavior. Organometallics 2017, 36, 4611−4619. (26) Falivene, L.; Credendino, R.; Poater, A.; Petta, A.; Serra, L.; Oliva, R.; Scarano, V.; Cavallo, L. SambVca 2. A web tool for analyzing catalytic pockets with topographic steric maps. Organometallics 2016, 35, 2286−2293. (27) The topographic steric map was produced by the web application SambVca 2.0. For the details of steric effects analysis, please refer to https://www.molnac.unisa.it/OMtools/sambvca2.0/. The default parameters embedded in the software were used in the current work. J

DOI: 10.1021/acs.organomet.8b00906 Organometallics XXXX, XXX, XXX−XXX