Dianmin Zhang1, Xiaofang He1, Tao Yang1, Song Liu1,2. 1. Chongqing Key Laboratory of Environmental Materials and Remediation Technologies College of Chemistry and Environmental Engineering, Chongqing University of Arts and Sciences, Chongqing, 402160, China. 2. School of Chemistry and Chemical Engineering, Chongqing University, Chongqing, 400030, China.
Abstract
A computational study has been performed to investigate the mechanism of RhIII -catalyzed C-H bond activation using sulfoxonium ylides as a carbene precursor. The stepwise and concerted activation modes for sulfoxonium ylides were investigated. Detailed theoretical results showed that the favored stepwise pathway involves C-H bond activation, carbonization, carbene insertion, and protonation. The free energy profiles for dialkylation of 2-phenylpyridine were also calculated to account for the low yield of this reaction. Furthermore, the substituent effect was elucidated by comparing the energy barriers for the protonation of meta- and para-substituted sulfoxonium ylides calculated by density functional theory.
A computational study has been performed to investigate the mechanism of RhIII -catalyzed C-H bond activation using sulfoxonium ylides as a carbene precursor. The stepwise and concerted activation modes for sulfoxonium ylides were investigated. Detailed theoretical results showed that the favored stepwise pathway involves C-H bond activation, carbonization, carbene insertion, and protonation. The free energy profiles for dialkylation of 2-phenylpyridine were also calculated to account for the low yield of this reaction. Furthermore, the substituent effect was elucidated by comparing the energy barriers for the protonation of meta- and para-substituted sulfoxonium ylides calculated by density functional theory.
Construction of C−C bonds is a major area of both academic and industrial chemistry. The ubiquity and low cost of hydrocarbons make C−H bond functionalization an attractive alternative to classical transition‐metal‐catalyzed cross‐coupling reactions with functionalized organic compounds.
In recent years, transition‐metal‐catalyzed C−H bond activation has developed into an efficient, atom‐ and step‐economical strategy for C−C bond construction.
In particular, metal‐carbene insertion into C(sp2)−H bonds has emerged as a common approach for constructing C(sp2)−C(sp3) bonds.
This strategy is mainly considered to follow a pathway comprising C(sp2)−H metalation, metal carbene formation, carbene migratory insertion, and protonation (Scheme 1).
Scheme 1
C−H bond functionalization with metal carbenes.
C−H bond functionalization with metal carbenes.In recent years, Cp*RhIII (Cp*=1,2,3,4,5‐pentamethylcyclopentadienyl) complexes have shown significant catalytic activity in C−H functionalization.
A wide range of carbene precursors, including α‐diazo carbonyls, hydrazones, triazo compounds, cyclopropenes, and ketone‐functionalized enynes, have been explored as coupling partners in Cp*RhIII‐catalyzed direct aryl C−H activation to construct C(sp2)−C(sp3) bonds. Direct dediazoniation of α‐diazo carbonyls in the presence of Cp*RhIII generates Cp*RhIII−Carbene complexes (Scheme 2a).
The Cp*RhIII−Carbene complex results from dediazoniation of the diazo substrate generated in situ from a hydrazone in the presence of a base (Scheme 2b).
Additionally, 1,2,3‐triazoles share the characteristics of diazo compounds because of the equilibrium between the closed and open forms (Scheme 2c).
Cyclopropenes can also be used as vinyl Cp*RhIII−Carbene precursors through ring opening (Scheme 2d).
The Lewis‐acidic RhIII center can readily initiate a 5‐exo‐dig cyclization in ketone‐functionalized enynes to afford Cp*RhIII−Carbene complexes (Scheme 2e).
Furthermore, coordination of iodonium ylide and subsequent elimination of PhI can also give Cp*RhIII−Carbene (Scheme 2f).
Scheme 2
Activation mode of carbene precursors in Cp*RhIII‐catalyzed C−H bond functionalization.
Activation mode of carbene precursors in Cp*RhIII‐catalyzed C−H bond functionalization.In 2017, Aϊssa
and Li
independently reported Cp*RhIII‐catalyzed arene C−H bond activation and acylmethylation using sulfoxonium ylides
as a carbene precursor. The products were isolated in good to high yields from a range of 2‐sulfoxonium ylides bearing aryl, heteroaryl, cyclopropyl, and adamantyl ketone groups (Scheme 3). The previously reported plausible catalytic cycle, based on experimental observations, features C−H activation, migratory insertion of the ylide into the carbon−metal bond, and protodemetalation, but the detailed activation mode for the α‐carbonyl sulfoxonium ylides remains unclear. Whether migratory insertion of α‐carbonyl sulfoxonium ylides into C(aryl)‐RhIII bonds occurs through a concerted or stepwise mechanism is disputed.
Scheme 3
Cp*RhIII‐catalyzed C−H bond activation using a sulfoxonium ylide as a carbene precursor. [a] Cp*=1,2,3,4,5‐pentamethylcyclopentadienyl. [b] HFIP=1,1,1,3,3,3‐hexafluoro‐2‐propanol.
Cp*RhIII‐catalyzed C−H bond activation using a sulfoxonium ylide as a carbene precursor. [a] Cp*=1,2,3,4,5‐pentamethylcyclopentadienyl. [b] HFIP=1,1,1,3,3,3‐hexafluoro‐2‐propanol.Herein, we report a computational study of the mechanism of the reaction reported by Aϊssa.
The detailed reaction mechanism of this Cp*RhIII‐catalyzed arene C−H bond acylmethylation reaction using sulfoxonium ylides as a carbene precursor was elucidated. The potential energy surface for dialkylation of the 2‐phenylpyridine substrate was also calculated to account for the low yield of the dialkylation product reported by Li.
Furthermore, Hammett plots with some representative substituents on α‐carbonyl sulfoxonium ylides were applied to study the electronic effect of the rate‐determining step.
Computational Methods
All of the density functional theory (DFT) calculations were performed with using the Gaussian 09 series of programs.
The B3‐LYP
density functional with the standard 6‐31G(d) basis set (SDD basis set for Rh and I) was used for the geometry optimizations. Harmonic frequency calculations were performed for all of the stationary points to confirm whether they are as local minima or transition structures, and to derive the thermochemical corrections for the enthalpies and free energies. All of the minima had zero imaginary frequency and all of the transition states had only one imaginary frequency. The M06 functional with larger basis set 6‐311+G(d) (SDD basis set for Rh and I) was used to calculate the single point energies with 1,1,1,3,3,3‐hexafluoro‐2‐propanol as the solvent using a continuum solvation model (SMD).
For the solvent data not included in Gaussian 09, the key parameters were set in the input file manually (for 1,1,1,3,3,3‐hexafluoro‐2‐propanol, eps=16.70, while the other parameters were the same as for 2‐propanol). The reported energies are the M06‐calculated Gibbs free energies in 1,1,1,3,3,3‐hexafluoro‐2‐propanol solvent. The optimized structures were visualized using CYLview.
Results and Discussion
The proposed catalytic cycles for Cp*RhIII‐catalyzed C−H bond activation using sulfoxonium ylides as a carbene precursor are shown in Scheme 4. Two possible pathways, with cationic Cp*RhIII complex A set as the active catalytic complex, were considered during computational modeling of the reaction. Coordination of the 2‐phenylpyridine substrate to cationic Cp*RhIII complex A and subsequent electrophilic deprotonation affords five‐membered rhodacyclic intermediate B. Subsequent coordination of sulfoxonium ylide 2 gives alkyl‐RhIII species C, which then undergoes β‐elimination of dimethyl‐sulfoxide (DMSO) to afford reactive carbene‐RhIII species D (pathway A). Migratory insertion of the activated carbene into the aryl‐Rh bond generates alkyl‐RhIII species E. Finally, acylmethylation product 3 is generated by protonation of intermediate E by acetic acid, and the active catalyst, cationic Cp*RhIII complex A, is regenerated to complete the catalytic cycle. In pathway B, an SN2‐type displacement of DMSO by the aryl group in intermediate C gives alkyl‐RhIII species E. Both pathways were investigated by DFT calculations.
Scheme 4
Proposed catalytic cycle for Cp*RhIII‐catalyzed C−H bond activation using a sulfoxonium ylide as a carbene precursor.
Proposed catalytic cycle for Cp*RhIII‐catalyzed C−H bond activation using a sulfoxonium ylide as a carbene precursor.The calculated Gibbs energy profile for pathway A is shown in Figure 1, in which cationic Cp*RhIII species CP1 is set as the relative zero point. Coordination of 2‐phenylpyridine substrate 1 to the RhIII center gives intermediate CP2, which is 5.7 kcal mol−1 exergonic. Subsequent pyridyl‐directed electrophilic deprotonation occurs via transition state TS1 with an energy barrier of 17.0 kcal mol−1, from which five‐membered rhodacyclic intermediate CP3 is generated. In transition state TS1, the lengths of the breaking C−H and forming O−H bonds are 1.32 and 1.33 Å, respectively, while the length of the forming C−Rh bond is 2.23 Å. The free energy of CP3 is 5.2 kcal mol−1 higher than that of CP2. The calculated results indicate that the first C−H bond activation is reversible, which is consistent with the results of H/D exchange experiments. Coordination of sulfoxonium ylide 2 to the RhIII center in CP3 gives alkyl‐RhIII intermediate CP4, which is 10.3 kcal mol−1 exergonic.
Figure 1
Gibbs energy profile and structural information for pathway A of Cp*RhIII‐catalyzed C−H bond activation using a sulfoxonium ylide as a carbene precursor. Values for bond lengths are given in angstroms (Å).
Gibbs energy profile and structural information for pathway A of Cp*RhIII‐catalyzed C−H bond activation using a sulfoxonium ylide as a carbene precursor. Values for bond lengths are given in angstroms (Å).Subsequent β‐elimination of DMSO in CP4 occurs via transition state TS2 to afford reactive carbene‐RhIII intermediate CP5. The structural information regarding TS2 shows that the length of the breaking C−S bond is 2.45 Å. Migratory insertion of the carbene into the aryl‐Rh bond generates six‐membered alkyl‐RhIII species CP6 via transition state TS3 with an energy barrier of only 3.6 kcal mol−1. The structural information of TS3 shows that the lengths of the forming C−C and the breaking C−Rh bonds are 2.39 and 2.10 Å, respectively. Coordination of acetic acid to the RhIII center in CP6 gives CP7, which is 6.3 kcal mol−1 endergonic. Subsequent protonation of the alkyl‐RhIII bond by acetic acid occurs via transition state TS4 with an overall activation free energy of 28.4 kcal mol−1 to give CP8. In TS4, the lengths of the forming C−H and breaking O−H bonds are 1.43 and 1.26 Å, respectively. The calculated results show that the protonation process is the rate‐determining step in the Cp*RhIII‐catalyzed C−H bond acylmethylation reaction, which is consistent with experimental observations. Final ligand exchange of 2‐phenylpyridine substrate 1 with acylmethylation product 3 gives active catalytic intermediate CP2, completing the catalytic cycle.The calculated free energy profiles of the SN2‐type displacement pathway are shown in Figures 2 and 3. As shown in Figure 2, coordination of the carbonyl group in sulfoxonium ylides 2 to the RhIII center in CP3 gives intermediate CP9, which is 6.2 kcal mol−1 exergonic. Subsequent SN2‐type displacement of DMSO by the aryl group in intermediate CP9 occurs via concerted transition state TS5 to generate enol‐RhIII species CP10.The energy barrier for this step is 51.5 kcal mol−1, which indicates that this pathway is unfavorable. The geometric information of TS5 shows that the lengths of the forming C−H and breaking C−S bonds are 1.67 and 1.88 Å, respectively, while the length of the breaking C−Rh bond is 2.25 Å. The structural information of TS5 indicates a concerted process.
Figure 2
Gibbs energy profile and structural information for pathway B of Cp*RhIII‐catalyzed C−H bond activation using a sulfoxonium ylide as a carbene precursor. Values for bond lengths are given in angstroms (Å).
Figure 3
Gibbs energy profile and structural information for pathway B of Cp*RhIII‐catalyzed C−H bond activation using a sulfoxonium ylide as a carbene precursor. Values for bond lengths are given in angstroms (Å).
Gibbs energy profile and structural information for pathway B of Cp*RhIII‐catalyzed C−H bond activation using a sulfoxonium ylide as a carbene precursor. Values for bond lengths are given in angstroms (Å).Gibbs energy profile and structural information for pathway B of Cp*RhIII‐catalyzed C−H bond activation using a sulfoxonium ylide as a carbene precursor. Values for bond lengths are given in angstroms (Å).As shown in Figure 3, the coordination of sulfoxonium ylide 2 to the RhIII center by the sulfoxide group gives complex CP11, which is 1.0 kcal mol−1 endergonic. The subsequent concerted SN2‐type displacement process occurs via transition state TS6 with an energy barrier of 37.9 kcal mol−1 to generate CP6. The free energy of TS6 is 33.9 kcal mol−1 higher than that of TS2, which indicates that this concerted SN2‐type displacement pathway is also unfavorable.In Li's work,
when 2‐phenylpyridine was used as the substrate, Zn(OAc)2 as additive, both mono‐ and dialkylation products were obtained in a 2 : 1 ratio (Scheme 5). DFT calculations are performed to account for this observation. The calculated potential energy surface for dialkylation of the 2‐phenylpyridine substrate is shown in Figure S1 (see the Supporting Information for details). Pyridyl‐directed electrophilic deprotonation of CP8 occurs via transition state TS7 with an energy barrier of 21.3 kcal mol−1 to give five‐membered rhodacyclic intermediate CP12. The free energy of TS7 is 5.6 kcal mol−1 higher than that of TS1, which results in the low yield of the dialkylation product.
Scheme 5
RhIII‐catalyzed C−H mono‐ and dialkylation of arenes using a sulfoxonium ylide as a carbene precursor.
RhIII‐catalyzed C−H mono‐ and dialkylation of arenes using a sulfoxonium ylide as a carbene precursor.Hammett plots
with some representative substituents (NH2, CH3, OMe, F, Cl, Br, I, CN, and NO2) were then applied to study the electronic effect of the rate‐determining step in the Cp*RhIII‐catalyzed C−H bond acylmethylation reaction. As shown in Figure 4, for both the meta‐ and para‐positions of the aryl ketones in the sulfoxonium ylide, a negative Hammett effect was found. The calculated results indicate that an electron‐donating substituent on the aryl ketone in sulfoxonium ylides would increase the ionic character of the C(alkyl)‐RhIII bond, leading to a lower protonation barrier, which is consistent with the experimental results showing that electron‐donating substituents are favorable in terms of the reaction yield. The slope of the calculated Hammett plot for meta‐position substitution of the aryl ketone in sulfoxonium ylides was 0.15, while that for the para position was 0.09. The greater conjugate effect of meta‐position substitution than para‐position substitution on the ionic character of the C−Rh bond in the protonation transition state leads to the lower slope. Furthermore, meta‐position substitution showed a better correlation than para‐position substitution.
Figure 4
Hammett plots for the (a) meta‐position and (b) para‐position in the protonation process.
Hammett plots for the (a) meta‐position and (b) para‐position in the protonation process.
Conclusion
In summary, the M06 DFT method has been used to examine the mechanism of Cp*RhIII‐catalyzed C−H bond activation using sulfoxonium ylide as a carbene precursor. The theoretical calculations indicated that this reaction proceeds through C−H bond activation, carbonization, carbene insertion, and protonation. The rate‐determining step was found to be protonation, which is consistent with the experimental results. The calculated results showed that activation of the sulfoxonium ylide occurs through a stepwise pathway rather than a concerted pathway. The high activation energy of the second C−H bond activation leads to the low experimental yield of the dialkylation product. Furthermore, Hammett plot calculations indicated that the rate‐determining protonation step was more efficient if an electron‐donating substituent was present on the aryl ketone in the sulfoxonium ylide. In addition, meta‐position substitution of the aryl ketone had a greater effect on the rate of the protonation step than para‐position substitution.
Conflict of interest
The authors declare no conflict of interest.As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re‐organized for online delivery, but are not copy‐edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors.Supporting InformationClick here for additional data file.
Authors: Sharon R Neufeldt; Gonzalo Jiménez-Osés; John R Huckins; Oliver R Thiel; K N Houk Journal: J Am Chem Soc Date: 2015-08-03 Impact factor: 15.419