Bijan
Mirabi
abc,
Mark
Lautens
c and
Mu-Hyun
Baik
*ab
aDepartment of Chemistry, Korean Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea. E-mail: mbaik2805@kaist.ac.kr
bCenter for Catalytic Hydrocarbon Functionalizations, Institute for Basic Science (IBS), Daejeon 34141, Republic of Korea
cDepartment of Chemistry, Davenport Chemical Laboratories, University of Toronto, Toronto, Ontario M5S 3H6, Canada
First published on 6th March 2025
Metal-catalyzed β-eliminations are elementary reaction mechanisms commonly leveraged in organometallic processes, including the renowned Mizoroki–Heck reaction. Although β-hydride elimination has traditionally been the focus of study, β-heteroatom elimination, in particular β-fluoride elimination, has seen a significant rise in contemporary organic methodologies. Intriguingly, rhodium(I) and palladium(II), which are isoelectronic, display opposite chemoselectivity for β-hydride vs. β-fluoride elimination. We investigated the origin of preferential β-fluoride over β-hydride elimination under rhodium(I) catalysis using density functional theory (DFT) calculations. Our modelling indicates that the kinetic preference is to undergo β-hydride elimination, but the observed chemoselectivity arises due to the reversible nature of the reaction. Additional modelling reveals that a Curtin–Hammett scenario enabled by reversibility of β-hydride elimination under the employed reaction conditions accounts for the enantioselectivity observed experimentally.
In medicinal chemistry, the incorporation of fluorine atoms provide benefits aside from improved metabolic stability and lipophilicity, influencing conformation8 and attenuating basicity.5c Fluorine atoms have been proposed to act as hydrogen bond acceptors in biological settings and to enable π-stacking interactions, playing a key role in increasing the potency of drug molecules compared to their non-fluorinated counterparts (Scheme 1a).9 Many fluorine-containing functional groups are bioisosteres of commonly utilized motifs in drug discovery (Scheme 1b).10 The earliest, and perhaps simplest, example of bioisosterism involving fluorine is the substitution of hydrogen atoms for fluorine atoms, again, typically to improve metabolic stability. Contemporary bioisosteres are more complex, containing varying topographies and interesting physicochemical properties, and have been used to replace more intricate functional groups. As an example, the gem-difluoroolefin has been leveraged as a carbonyl bioisostere and is present in various pesticides and bioactive compounds, such as phenstatin or artemisinin derivatives.11gem-Difluoroolefins also serve as reliable platforms for functionalization to other value-added molecules, often with high levels of selectivity.12 As such, methods enabling the synthesis of gem-difluoroolefins have seen remarkable advances and have been of increasing focus in the synthetic organic community.13
One method for the synthesis of this functionality is from the corresponding trifluoromethyl compounds. Classical conditions relied on SN2′-like reactions of strongly nucleophilic reagents with trifluoromethylolefins,14 Wittig-like,15 or Julia–Kocienski-like16 reactions. A milder method to access this valuable structural motif is via a transition metal-catalyzed β-fluoride elimination. Numerous transition metals have displayed the ability to perform this elementary step, wherein a metal alkyl species can undergo elimination to generate an alkene and the corresponding metal–fluoride.17
A meaningful understanding of this process can be obtained by comparison with the related β-hydride elimination reaction. The mechanistic underpinnings of the two elementary processes are similar. The β-fluoride elimination reaction is thought to occur via a syn co-planar conformation of the M–C–C–F dihedral angle, highly reminiscent of the orientation required in β-hydride elimination.17c,d,18 The presence of a vacant coordination site at the metal center enables the formation of a 3-center, 2-electron agostic interaction consisting of σC–F → d donation.19 Mechanistically, these features exactly mirror those for syn stereospecific β-hydride elimination. However, the key differences lie in the nature of the C–H bond compared to the C–F bond and the formation of the resulting products. First, the C–F σ-orbital is significantly polarized toward the fluorine atom due to its greater electronegativity and is a poor donor.20 Thus, the agostic interaction arising from a C–F bond is not as favourable as the agostic interaction from a C–H bond. It would follow that the C–F bond cleavage is less favoured from a thermodynamic perspective since it is stronger than a C–H bond. However, a general claim of this type cannot be made since the thermodynamics of a β-fluoride, in comparison to the corresponding β-hydride, elimination must also account for the other bonds forming in the reaction. If we assume that the strength of the metal–alkyl and metal–olefin bonds in both cases is approximately equal, then the thermodynamics will depend on the strength of the C–H and C–F bonds, as well as that of the metal–hydride and metal–fluoride bonds in the products. Although the strength of the C–F bond is greater than that of the C–H bond, the strength of the metal–fluoride bond21 is typically also greater than that of the corresponding metal–hydride bond,22 which could bias the thermodynamics in favour of β-fluoride elimination.23
Our interest in this realm lies specifically in the observed chemoselectivity in β-hydride versus β-fluoride elimination in isoelectronic palladium(II) and rhodium(I) catalyst systems. In particular, palladium(II)–hydrides have been invoked as intermediates in countless mechanisms involving β-hydride elimination, most commonly in Mizoroki–Heck reactions.24 Concurrently, β-fluoride elimination as an elementary process in catalysis involving palladium(II) has been established in the 90s25 and has since been proven to be a facile process enabling syntheses of fluorinated alkenes.26 Early computational work by Lin on model systems determined that there is a kinetic preference for palladium(II) to undergo β-hydride elimination over β-fluoride elimination.27 Experimental and further computational evidence supporting these findings was later obtained by the Altman and Cheong groups in a copper- and palladium-catalyzed arylation reaction of gem-difluoroolefins using arylsulfonyl chlorides (Scheme 2a, top).28 Following arylation of the gem-difluoroolefin, the resulting alkylpalladium complex could either undergo β-hydride or β-fluoride elimination, yet perfect selectivity was observed for the hydride elimination. The Morandi group similarly observed exclusive formation of olefin products arising from β-hydride elimination when competitive β-fluoride elimination was viable (Scheme 2a, bottom).29
In 2008, Murakami reported the first example of a rhodium-catalyzed arylative defluorination reaction between aryl boron reagents and α-trifluoromethylstyrenes.30 Interestingly, the reaction led to the near exclusive formation of the gem-difluoroolefin arising from β-fluoride elimination, despite competing β-hydride elimination available as an alternative pathway. Although the authors highlight that the reaction proceeds with high chemoselectivity for the β-fluoride elimination product, no rationale or explanation was provided. The exploitation of β-fluoride elimination reactions following a rhodium(I)-catalyzed C–C bond forming migratory insertion reaction of electron-poor alkenes was not explored further until 2016, when Hayashi reported the asymmetric variant.31 As in the case of Murakami, the reaction displayed remarkable chemoselectivity for the β-fluoride elimination pathway. The Lautens group later reported a rhodium-catalyzed defluorinative arylation reaction of secondary amides.32 Utilizing aryl potassium trifluoroborate salts and β-trifluoromethylacrylamides as electronically-biased Michael acceptors, the rhodium-catalyzed enantioselective arylation reaction occurs with inverted regiochemistry via polarity inversion of the α-position.33 Interestingly, the alkylrhodium(I) intermediate could undergo a β-hydride elimination to afford the trisubstituted acrylamide, however products arising from this mechanistic pathway were never observed.
Palladium and rhodium are arguably two of the most privileged metal catalysts in organic synthesis, and β-eliminations are commonly observed for both metals. While the experimental trends observed in the literature show a clear preference for palladium to undergo β-hydride elimination and for rhodium to undergo β-fluoride elimination, there exists a fundamental gap accounting for this reactivity difference. We were interested in determining the origin of chemoselectivity for β-fluoride elimination over β-hydride elimination in rhodium(I) catalysis (Scheme 2b). To the best of our knowledge, the methods disclosed by Murakami, Hayashi, and our group are the only known examples of β-fluoride elimination occurring at rhodium(I) where, in theory, competing β-hydride elimination can also take place. Furthermore, all three reported systems are significantly different with respect to catalyst structure, although all three utilize diene ligands, and substrate choice. As described earlier, computational work by Lin and experimental work by Altman and Cheong, and Morandi both support that palladium catalysts prefer β-hydride elimination. Considering this preference of palladium to perform β-hydride elimination when β-fluoride elimination is also possible, we were curious why the same inclination was not also observed with rhodium(I) catalysts. We were particularly intrigued in the orthogonal behaviour between rhodium(I) and palladium(II) since these complexes are found on the same row of the periodic table, precluding differences due to relativistic effects34 or atomic radii, and are isoelectronic d8-metal complexes and should display similar behaviour.
The significant preference for β-fluoride elimination in these systems renders experimental investigation difficult. In the studies by Hayashi and Murakami, products arising from β-hydride elimination pathways are formed as minor products in yields of less than 6%. The report by Lautens occurred with even greater chemoselectivity and no products arising from this mechanistic manifold were observed. Experiments deducing why one of these products is not formed are not simple to design and significant perturbation to the system would not provide adequate evidence for the observed reactivity. Thus, we reasoned that a computational approach would be more appropriate and would enable investigation into the apparent selectivity of rhodium(I) to undergo β-fluoride elimination.
The energy components were computed with the following protocol. The free energy in solution, Gsol, has been calculated as follows with T = 298.15 K, T = 373.15 K, T = 333.15 K, or T = 303.15 K to match the experimental conditions for each system under study.
Gsol = Ggas + Gsolv | (1) |
Ggas = Hgas + TSgas | (2) |
Hgas = ESCF + ZPE | (3) |
ΔESCF = ∑ ESCF, products – ∑ ESCF, reactants | (4) |
ΔGsol = ∑ Gsol, products – ∑ Gsol, reactants | (5) |
Our calculations indicate that the β-hydride elimination path traversing 2-TS (ΔG‡ = 6.5 kcal mol−1) is more favoured than the β-fluoride elimination occurring via3-TS (ΔG‡ = 9.4 kcal mol−1, ΔΔG‡ = 2.9 kcal mol−1). This somewhat surprising result, considering the experimental observations, suggests that the rhodium(I) complex 1 has a kinetic preference for β-hydride elimination for the model system. Further examination of the reaction coordinate reveals that the formation of 4via β-hydride elimination is endergonic by 5.5 kcal mol−1, relative to 2. However, the formation of 5 from 2 is endergonic by only 2.8 kcal mol−1. This suggests that the formation of the rhodium(I)–fluoride species 5via β-fluoride elimination could be irreversible, while the formation of the rhodium(I)–hydride 4 can be reversible. Although not evaluated computationally, the resulting rhodium(I)–fluoride complex can react with an equivalent of base or an organoboron species via transmetalation in an irreversible fashion, driving the reaction. On the other hand, rhodium(I)–hydride complexes do not engage in this type of behaviour and, in the absence of an additional driving force, would terminate catalysis if re-insertion to the bound olefin was not viable. Taken together, the initial computational results suggest that the observed chemoselectivity for β-fluoride elimination under rhodium(I) catalysis may not be a kinetic preference, but an artifact of a thermodynamically-driven capture of a kinetically unfavorable intermediate. If this is the case, the observed product ratio between 4 and 5 would be approximately 1:
35 at 100 °C based on the computed ΔΔG = −2.6 kcal mol−1, which is more consistent with the experimental observation of real systems.
Intermediates 7 and 8 serve as common intermediates for both β-elimination pathways and were used as starting points for examination of the β-hydride elimination (Scheme 4b). In this case, two possible pathways for β-hydride elimination exist, one forming the (E)-olefin and one forming the (Z)-olefin. Experimentally, only the (E)-olefin was observed but both pathways were calculated. Intermediate 11 features a C–H agostic interaction that enforces a trans relationship between the two phenyl rings and is exergonic by 3.6 kcal mol−1 relative to 8. β-Hydride elimination via11-TS has an activation barrier of 17.6 kcal mol−1 affording the (Z)-olefin–rhodium(I) complex 12, a process that is endergonic by 6.8 kcal mol−1. In contrast, the formation of the C–H agostic interaction that places the two phenyl rings cis relative to each other 13 is slightly less exergonic (−2.8 kcal mol−1 relative to 8). In this case, the β-hydride elimination transition state 13-TS has an activation barrier of only 3.2 kcal mol−1, affording the (E)-olefin–rhodium(I) product 14 in a process that is endergonic by 2.7 kcal mol−1. These results correlate well with experiments and suggest that the preferred product arising from β-hydride elimination should indeed be the (E)-olefin.
Analysis of 9-TS for β-fluoride elimination (Scheme 4a) and 13-TS for β-hydride elimination (Scheme 4b) again reveals that the rhodium(I) catalyst has a kinetic preference for the unobserved β-hydride elimination pathway. Under the reaction conditions, which employ an organoboron reagent that is activated by an equivalent of MeMgCl, the rhodium(I)–fluoride complex 10 can undergo a transmetalation reaction in an irreversible fashion, which serves to drive the reaction toward the observed difluoroolefin. On the other hand, both of the rhodium(I)–hydride complexes 12 and 14 are formed in an apparent reversible fashion and cannot be funneled in a manner that would liberate the olefin and enable catalysis.
At the onset of our investigation, we were hoping to find a unified explanation that could be applied for all three of the known experimental systems. Our analysis thus far has revealed that the observed selectivity for β-fluoride elimination over β-hydride elimination is likely due to the reactions occurring under thermodynamic control. While this explanation is understandable for the model system and Murakami's achiral reaction,30 the observed product selectivity being a result of thermodynamics is not tractable for the method reported by Hayashi31 or Lautens.32 More specifically, these reactions are highly enantioselective. Thus, if the reaction is driven by thermodynamics under equilibrating conditions, then observation of a racemic mixture must necessarily occur. Therefore, we surmised that an alternative explanation must exist for the enantioselective rhodium-catalyzed defluorinative reactions.
![]() | ||
Scheme 5 (a) Enantioselective rhodium-catalyzed arylation–defluorination reaction of trifluoromethyl olefins and (b) β-hydride elimination pathway accounting for side product formation. (c) Reaction coordinate for the enantioselective defluorinative arylation developed by Hayashi in 2016.31 |
We began our investigation of the reaction by examining the rhodium–trifluoromethylolefin complex 15 (Scheme 5c). The enantioselectivity of the reaction is determined during the C–C bond forming migratory insertion of the olefin into the rhodium(I)–aryl bond. The transition state leading to the product with (S)-absolute configuration, 15-TS-(S), was found to have an activation barrier of 9.6 kcal mol−1. By comparison, 15-TS-(R) leading to the enantiomeric product had a much larger activation energy of 22.4 kcal mol−1. The alkyl rhodium(I) complex 16 features η2-coordination of the aryl group. Displacement of the aryl ring by the C–H σ-bond leading to agostic complex 17 was endergonic by 11 kcal mol−1, which undergoes a β-hydride elimination traversing 17-TS, leading to olefin-bound rhodium species 18. Notably, this process has an activation barrier of 12.5 kcal mol−1 from 16. By comparison, formation of the β-fluoride agostic complex 19 is endergonic by 15 kcal mol–1 from 16. The β-fluoride elimination occurs via19-TS leading to the difluoroolefin-bound rhodium–fluoride 20, with an activation barrier of 30.0 kcal mol−1.
Although these results qualitatively support the experimental conditions, which require heat to promote the reaction, there are two troubling observations. Firstly, the quantum chemical calculations once again support that the β-hydride elimination pathway should be favourable over the observed β-fluoride elimination pathway. Secondly, the calculations predict that the enantiomer that should form during the course of the reaction is opposite to that experimentally observed. We will first examine the defluorinative reaction disclosed by Lautens before returning to this alarming finding.
The reaction developed by Lautens utilizes a C1-symmetric catalyst and β-trifluoromethyl acrylamides. The conjugated olefin could react with the C1-symmetric catalyst in either the s–cis or s–trans conformation at either enantiotopic face, leading to 8 possible encounter complexes for the initial migratory insertion. Only the lowest energy transition state leading to each enantiomer will be considered for this discussion (see the ESI† for more details). The lowest energy encounter complex 25 was chosen as the reference point (Scheme 6c).
The rhodium–olefin complex 26 at an energy of 6.1 kcal mol−1 was found to lead to the lowest energy migratory insertion transition state 26-TS at 10.4 kcal mol−1. 26-TS leads to the observed enantiomer and is 3.9 kcal mol−1 lower in energy than the lowest energy enantiomeric transition state (see the ESI† for more details), in excellent agreement with the experimental observations. The direct product arising from 26-TS is the alkyl–rhodium(I) complex coordinated by the arene in an η2-fashion 27, found at −14.5 kcal mol−1. Formation of the β-C–F agostic interaction 28 is endergonic by 17.5 kcal mol−1, and subsequent β-fluoride elimination 28-TS was higher in energy by another 4.7 kcal mol−1. This process forming the rhodium(I) complex bound to the product difluoroolefin 29 was exergonic by 4.6 kcal mol−1 relative to 28. The β-fluoride elimination occurs with an activation barrier of 22.2 kcal mol−1 relative to the η2-coordinated species 27.
Instead of forming the C–F agostic interaction, intermediate 27 could also undergo ligand substitution of the bound arene and form the C–H agostic interaction. Since we are using 27 as the common intermediate for these two competitive pathways, the activation energy of the β-hydride elimination transition state 27-TS was calculated relative to 27 with an activation energy of 17.7 kcal mol−1, which is lower in energy than the β-fluoride elimination transition state 28-TS (Scheme 7a). This result is unsurprising based on our findings thus far, but indicates that the preferred pathway involves a stereoablative β-elimination leading to the achiral intermediate 30 in an endergonic fashion. An olefin facial swap occurring via decoordination and recoordination leads to 31 and is slightly exergonic. A hydrorhodation traversing 31-TS with an activation energy of 15.5 kcal mol−1 leads to the C–H agostic intermediate 32. An endergonic ligand exchange to the C–F agostic complex 33 sets the stage for another β-fluoride elimination occurring via33-TS (ΔG‡ = 7.4 kcal mol−1 relative to 33) affording the rhodium–fluoride complex 34.
![]() | ||
Scheme 7 (a) Potential energy surface for competing β-hydride elimination and β-fluoride elimination, and (b) the Curtin–Hammett scenario explaining selectivity. |
At first glance, this profile appears to suggest that the absolute configuration of the final product should be opposite to that which is observed since the β-hydride elimination pathway offers a low energy pathway that removes the chiral information imparted in the initial phase of the reaction. However, careful examination of the β-hydride elimination and reverse hydrorhodation steps suggests, again, that these processes should be reversible under the conditions of the reaction. Thus, if the two diastereomeric β-C–F agostic complexes 33 and 28 exist in equilibrium via these intermediary rhodium–olefin and rhodium–hydride complexes, then the observed product distribution should be governed by Curtin–Hammett kinetics (Scheme 7b).47 In other words, the difference in activation energy between 28-TS and 33-TS, measured from the respective equilibrating β-C–F agostic complexes, will determine the observed selectivity. The computed ΔΔG‡ of the reaction is 2.7 kcal mol−1 in favour of the enantiomer with (S)-absolute configuration, in alignment with experimental observations. Interestingly, the ablation of the stereocenter installed in the migratory insertion step during the β-hydride elimination process suggests that the enantio-determining step is actually the β-fluoride elimination and not the C–C bond forming reaction.
Based on these results, we re-examined the Hayashi reaction discussed earlier. We discovered that our quantum chemical calculations were unable to account for the observed enantioselectivity based on the fact that the preferred C–C bond forming migratory insertion was favoured to give the experimentally unobserved isomer. In all of the cases examined in this work, β-hydride elimination has consistently been the kinetically favoured pathway, providing the catalyst an alternative pathway to enantioselectivity. We re-examined the potential energy surface of the defluorinative arylation of trifluoromethyl alkenes (Scheme 8).
![]() | ||
Scheme 8 Curtin–Hammett scenario accounting for enantioselectivity in the defluorinative arylation reaction reported by Hayashi in 2016.31 |
As discussed earlier, the β-hydride elimination pathway occurring via17-TS is significantly favoured over β-fluoride elimination 19-TS. Facial swapping of the bound olefin converts rhodium–hydride 18 to its isomer 35. Hydrorhodation via35-TS leads to the β-hydride agostic complex 36. Ligand substitution to the β-fluoride agostic complex 37 is endergonic by 1.9 kcal mol−1. Subsequent β-fluoride elimination traversing 37-TS leads to the difluoroolefin-bound rhodium–fluoride complex 38. Notably, dissociation of the product from 38 would lead to the observed enantiomer of the product. If we again consider that the β-hydride elimination/re-insertion pathway acts as a means of equilibration between diastereomeric rhodium complexes 19 and 37 containing β-fluoride agostic interactions, then the enantioselectivity would be determined by the difference in activation energy between 19-TS and 37-TS, and not the C–C bond forming migratory insertion (15-TS-(R) and 15-TS-(S)). Based on this, the favoured product is the one that forms through 37-TS (ΔΔG‡ = 6.8 kcal mol−1), which is indeed the enantiomer observed by Hayashi.
![]() | ||
Scheme 9 Explanation for the different reactivity between (a) rhodium(I) and (b) palladium(II) when comparing β-hydride and β-fluoride elimination. |
This work discloses a comparison on rhodium(I) and palladium(II) species, which suggests that chemoselectivity can be obtained in β-hydride vs. β-fluoride elimination by the choice of metal catalyst – palladium(II) would lead to preferential β-hydride elimination, whereas rhodium(I) would lead to preferential β-fluoride elimination. Importantly, this work revealed that the chemoselectivity in this reaction class can be due to mechanistic pathways available to the metal catalyst and may not necessarily be affected by ligand choice. Thus, future methodologies studying chemoselective β-hydride and β-X elimination should consider which mechanistic pathways are available to the metal and not overly rely on ligand discovery, if the metal's intrinsic reactivity cannot be overridden or diverted. Broadly speaking, we anticipate that future work can utilize the disclosed results in designing new synthetic methodologies involving β-elimination with rhodium, palladium, and other organometallic species. We believe that this work provides fundamental understanding for elementary organometallic reaction mechanisms and will serve as a framework for methodology development leveraging chemoselective β-elimination as a crucial step in catalysis.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4cy01495b |
This journal is © The Royal Society of Chemistry 2025 |