Yang
Li
and
Zhenyang
Lin
*
Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China. E-mail: chzlin@ust.hk; Fax: (+852) 2358 1594
First published on 26th September 2014
Density functional theory calculations were carried out to study the Pd-catalysed oxidation reactions of alcohols using phenyl chloride as an oxidant. Our calculations supported that the mechanism mainly involves oxidative addition, β-hydride elimination, reductive elimination, and finally ligand substitution. Through our calculations, we have explained why oxidation of secondary alcohol was experimentally observed but not that of primary alcohol. The oxidation products (ketones) of secondary alcohols bind much more weakly than the oxidation products (aldehydes) of primary alcohols, contributing to the reactivity difference. The mechanism of the Pd-catalysed domino oxidation–arylation reactions of secondary alcohols was also studied. We have explained the experimental observation that α-arylation of the oxidation products (ketones) of secondary alcohols occurs only when the temperature was raised to 80 °C from below 40 °C.
Development of catalyst systems that can control the ability to selectively oxidize an alcohol in the presence of other groups susceptible to oxidation has been identified as one of the key challenges in alcohol oxidations.11 In 2009, Navarro et al. found that NHC-Pd complexes were able to catalyse oxidation of secondary alcohols at mild temperatures when phenyl chloride was used as an oxidant (eqn (1)). Interestingly, the NHC-Pd complexes did not catalyse the oxidation of primary alcohols under the same reaction conditions. Although a reaction mechanism involving oxidative addition, β-hydride elimination, and reductive elimination has been proposed, it cannot account for the experimental observation that secondary alcohols are far more reactive than primary alcohols.12 On the basis of their work, the same research group reported in 2011 the domino oxidation–arylation reactions of secondary alcohols to produce α-arylated products under the same NHC-Pd catalytic system but at a higher temperature (80 °C) and in the presence of 2.1 equivalent of the oxidant phenyl chloride (eqn (2)).13
(1) |
(2) |
Theoretical studies reported so far on metal-catalysed oxidation of alcohols are limited. Baerends et al. theoretically investigated the aerobic oxidation of primary alcohols to aldehydes, using [Cu(bipy)]2+, the TEMPO radical, and a strong base as the catalytic system. The most important conclusion of this work is that oxidation of alcohols to aldehydes takes place by electrophilic attack of uncoordinated TEMPO+ on one C–H(α) bond of a primary alcohol.14 Goddard et al. theoretically studied the reactivity of (NHC)Pd(OAc)2 and emphasized the role played by the strong trans influence ligand NHC. They found that β-hydride elimination, in which the β-hydrogen of a palladium-bound alkoxide was transferred directly to the free oxygen of the bound carboxylate, provides the lowest-energy route.15
In this work, with the aid of DFT calculations, we will investigate the reactivity difference between primary and secondary alcohols in the Pd-catalysed oxidation reactions (eqn (1)) by studying the detailed reaction mechanisms. We hope to understand the factors affecting the reactivity difference. We will also study the mechanism related to the Pd-catalysed domino oxidation–arylation reaction shown in eqn (2) and try to understand how higher temperature promotes the α-arylation immediately after the alcohol oxidation. An understanding of the reaction mechanism is expected to lead to more efficient synthetic strategies and more efficient catalysts for the oxidation of alcohols to aldehydes or ketones.
Worth mentioning here is that in recent years, environment-friendly alcohol dehydrogenation that does not require an oxidant has been developed to synthesize amides and pyrroles.16 Wang et al. investigated the mechanism of Ru-catalysed reactions of amines with primary alcohols to produce amides.17 In these reactions, the dehydrogenation reactions normally proceed via the bifunctional double hydrogen transfer mechanism rather than the β-H elimination mechanism.
To reduce the computational cost, we used 1,3-dimethylimidazol-2-ylidine as the model NHC ligand, in which the substituents at N in the experimentally used NHC ligand (IPr = 1,3-bis-(2,6-diisopropylphenyl)imidazol-2-ylidene) were replaced by CH3. Our early work established the validity of using 1,3-dimethylimidazol-2-ylidine as a model for IPr.26 To provide further support for the use of the model NHC ligand, we also calculated the reaction free energies for 1E + PhCH2OH → 1A + PhCHO and 2E + PhEtCHOH → 2A + PhEtCO, which are important for our conclusions, with the full NHC ligand used in the experimental reactions. The results show that with the model NHC ligand the reaction free energies are 8.8 and 3.2 kcal mol−1 (Fig. 1 and 3). With the full NHC model, the reaction free energies are 8.7 and 3.9 kcal mol−1. We also calculated the free energy changes for 1A + PhCl → 1TS(A–B) and 2A + PhCl → 2TS(A–B) with the full NHC ligand. The free energy changes were calculated to be 17.0 and 17.7 kcal mol−1 (Fig. S1 and S2†), compared with 20.6 and 20.8 kcal mol−1 from the model NHC results. Additional comments on the relative barriers using different models are given in the ESI.†
Fig. 1 Energy profile calculated for the Pd-catalysed oxidation of PhCH2OH (PhCH2OH + PhCl + tBuOK → PhCHO + PhH + tBuOH + KCl). The relative free energies are given in kcal mol−1. |
Oxidative addition of phenyl chloride to complex A gives a four-coordinate, square planar Pd(II) complex (B).28 The complex B reacts with tBuOK to give tBuOH, KCl and the intermediate C containing an agostic interaction. The intermediate C then undergoes β-hydride elimination to give complex D, which is a precursor complex ready to undergo reductive elimination to give the co-product benzene molecule, and complex E in which the carbonyl product molecule acts as a ligand. Finally, a ligand substitution of substrate (alcohol) for product (aldehyde/ketone) occurs to regenerate the active species A and complete the catalytic cycle.
The energy profile shown in Fig. 1 indicates that the overall reaction is highly exergonic with a reaction free energy of −46.4 kcal mol−1 (PhCH2OH + PhCl + tBuOK → PhCHO + PhH + tBuOH + KCl). The highest energy structure in the energy profile corresponds to the transition state 1TS(A–B) which lies 20.6 kcal mol−1 higher in energy than the energy reference point (1A + PhCH2OH + PhCl + tBuOK). From the energy profile shown in Fig. 1, one may rush to conclude that the overall barrier for the catalytic reaction corresponds to the barrier calculated for the oxidative addition of PhCl to 1A, which is 20.6 kcal mol−1. This is not the case because the lowest energy point in the catalytic cycle corresponds to the intermediate 1E, not the end point of the catalytic cycle. We can see that 1E + PhCH2OH is more stable than 1A + PhCHO by 8.9 kcal mol−1 (Fig. 1). In fact, the overall barrier for the catalytic reaction should be 20.6 + 8.9 = 29.5 (kcal mol−1) according the energetic span model developed by Kozuch et al.29
To further elaborate the application of the energetic span model to the case here, we have depicted the schematic energy profile diagram in Fig. 2 by extending the energy profile to an additional catalytic cycle. In Fig. 2, RDTS and RDI stand for rate-determining transition state and intermediate, respectively. ΔGR is the reaction free energy for the catalytic reaction under consideration. From Fig. 2, we can also view the first RDI to the second RDI as a catalytic cycle. In this view, we can easily come to the conclusion that the overall barrier is ΔG1 + ΔG2, not a simple addition of the two relevant barriers (ΔG1 + ΔG3). In Fig. 1, we have that ΔG1 = 20.6 kcal mol−1 and ΔG2 = 8.9 kcal mol−1.
Fig. 2 Schematic illustration of the energy profile given in Fig. 1, which includes the energy profile for an additional catalytic cycle. |
Similar to what we found for the reaction of the primary alcohol PhCH2OH, the intermediate 2E having the product molecule PhEtCO as a ligand also corresponds to the lowest energy point in the energy profile. However, 2E + PhEtCHOH is more stable than 2A + PhEtCO by 3.2 kcal mol−1 only. In the reaction of the primary alcohol PhCH2OH, 1E + PhCH2OH is more stable than 1A + PhCHO by 8.9 kcal mol−1. Thus, the overall barrier for the catalytic reaction of the secondary alcohol PhEtCHOH is 20.8 + 3.2 = 24.0 (kcal mol−1). The overall reaction (PhEtCHOH + PhCl + tBuOK → PhEtCO + PhH + tBuOH + KCl) is exergonic with a reaction free energy of −47.6 kcal mol−1.
Fig. 3 Energy profile calculated for the Pd-catalysed oxidation of PhEtCHOH (PhEtCHOH + PhCl + tBuOK → PhEtCO + PhH + tBuOH + KCl). The relative free energies are given in kcal mol−1. |
The analysis given above indicates that the overall barrier for the catalysed oxidation of a given alcohol substrate molecule is the sum of (1) the barrier (ΔG1) for oxidative addition of PhCl to the active species LPd(0)(substrate) (1A or 2A) and (2) the reaction free energy (ΔG2) for the ligand substitution of an alcohol substrate molecule for an aldehyde or ketone product molecule in LPd(0)(product) to regenerate the active species.
Our calculations show that the barriers (ΔG1) for the oxidative additions of phenyl chloride to 1A and 2A are approximately the same. Examining the structures of 1A, 2A and the corresponding oxidative addition transition states 1TS(A–B) and 2TS(A–B) (Fig. 4), we see that the optimized structures of 1A and 2A as well as 1TS(A–B) and 2TS(A–B) are similar to each other. The results are consistent with the fact that approximately the same barriers were calculated for the two oxidative additions.
As discussed above, the ligand substitution of an alcohol substrate molecule for an aldehyde or a ketone product molecule is another step that contributed to the overall barrier. The energetics (ΔG2) associated with the ligand substitution step reflects how relatively easily an aldehyde or a ketone product molecule, once it is formed, is released from the catalyst to which it is coordinated. The calculation results show that 1E + PhCH2OH → 1A + PhCHO is more endergonic than 2E + PhEtCHOH → 2A + PhEtCO by 5.5 kcal mol−1. Here, we can see that it is the difference in (ΔG2) that contributes to the significant reactivity difference of the two substrate molecules.
Scheme 2 compares the binding energies of the alcohol substrate molecules and the corresponding aldehyde/ketone product molecules with the Pd(0) metal centre. A relatively small difference in the binding energy was calculated between the two alcohol substrate molecules. However, a noticeably significant difference in the binding energy can be found between the two product molecules. Clearly, the product molecule PhCHO derived from the former alcohol substrate molecule binds much more strongly to the Pd(0) metal centre than the product molecule PhEtCO derived from the latter alcohol substrate molecule. The stronger η2-binding of PhCHO versus PhEtCO with Pd(0) can be attributed to both steric and electronic effects. PhCHO is less sterically hindered than PhEtCO when it is η2-coordinated to the metal centre. Electronically, PhCHO has a lower lying CO π* orbital when compared with PhEtCO (Fig. 5), giving rise to greater Pd(0)-to-CO(π*) back-bonding interaction. In the calculated structures shown in Fig. 6, 1E gives a shorter Pd–C(O) distance than 2E.
Fig. 6 Selected bond distances (Å) calculated for the ligand substitution intermediates 1E and 2E. Hydrogen atoms have been omitted for the purpose of clarity. |
Scheme 2 The binding energies of the alcohol substrate molecules and the corresponding aldehyde/ketone product molecules with the Pd(0) metal centre. |
A plausible explanation for the anionic species [LPd(OtBu)]− not being involved as the active species is given as follows. Experimentally, tertbutoxide has very poor solubility in toluene, which is the solvent used in the reactions. Alcohol and the catalyst are much more soluble in toluene than tertbutoxide. Therefore, it is reasonable to assume that formation of LPd(0)(alcohol) is much faster than that of [LPd(OtBu)]−.
With LPd(0)(alcohol) as the active species, oxidative addition of PhCl has a reasonably small barrier and is expected to have a fast step also. Once LPd(II)Ph(Cl)(alcohol) is formed from the oxidative addition, the coordinated alcohol in this Pd(II) complex is much more acidic (when compared with the coordinated alcohol in LPd(0)(alcohol)) and easily deprotonated by tertbutoxide as proposed in Scheme 1.
Here, we use the product molecule PhEtCO (discussed above) as a substrate ketone to react with phenyl chloride (the α-arylation reaction). Fig. 3 shows that the last step 2E + PhEtCHOH → 2A + PhEtCO, which is a ligand substitution process, is reversible. Therefore, instead of undergoing the ligand substitution, 2E can further react with PhCl via oxidative addition to give a four-coordinate, square planar Pd(II) species 2F (Scheme 3). From the species 2F, a ligand substitution of tBuO− for Cl− gives the intermediate 2G. Subsequently, deprotonation at the α-carbon by the tBuO− ligand produces the η3 intermediate 2H and tBuOH. Then, reductive elimination in the 2H takes place to give the intermediate 2I in which the α-arylated product molecule acts as a ligand. Finally, a ligand substitution of a ketone molecule for the α-arylated product molecule regenerates the active catalytic species 2E and completes the α-arylation reaction.
Fig. 7 shows the energy profile calculated for the Pd-catalysed α-arylation reaction of PhEtCO based on the mechanism proposed in Scheme 3. The first step is the oxidative addition of phenyl chloride to the two-coordinate complex 2E to give the intermediate 2Fvia the transition state 2TS(E–F). The free energy barrier calculated for this step is 26.7 kcal mol−1. From 2F, a ligand substitution of tBuO− for Cl− occurs to give the intermediate 2G, in which the tBuO− ligand coordinates to the metal centre. The deprotonation from 2G proceeds to 2G2via the transition state 2TSG, which involves migration of one proton at the coordinated-PhEtCO α-carbon to the tBuO− ligand. In the 2G2 intermediate, tBuOH is easily dissociated to form the more stable intermediate 2H. Reductive elimination in 2H occurs to give the two-coordinated intermediate 2Ivia the transition state 2TS(H–I), in which the α-arylated product molecule is coordinated as a ligand to the metal centre in an η2 coordination mode. In the final step, a ligand substitution of PhEtCO for the α-arylation product regenerates the active catalytic species 2E and completes the α-arylation reaction.
Fig. 7 Energy profile calculated for the Pd-catalysed α-arylation reaction of PhEtCO. The relative free energies are given in kcal mol−1. |
The energy profile shown in Fig. 7 indicates that the overall reaction PhEtCO + PhCl + tBuOK → CH3CH(Ph)C(O)Ph + tBuOH + KCl is exergonic with a reaction free energy of −43.5 kcal mol−1. The highest energy structure in the energy profile corresponds to the transition state 2TS(E–F) which lies 26.7 kcal mol−1 higher in energy than the energy reference point (2E + PhEtCO + PhCl + tBuOK). The lowest energy structure corresponds to the intermediate 2I. 2I + PhEtCO is more stable than 2E + CH3CH(Ph)C(O)Ph by 3.2 kcal mol−1. As mentioned above, the overall barrier for the α-arylation reaction is 29.9 kcal mol−1, the sum of 26.7 kcal mol−1 and 3.2 kcal mol−1 according to what we discussed above (Fig. 2).
The calculation results indicate that the α-arylation of the oxidation product PhEtCO has a higher overall free energy barrier than the oxidation of PhEtCHOH. Therefore, it is expected that at a lower temperature, the reaction stops at the oxidation reaction while at a higher temperature the reaction can continue to give the α-arylated product. These results well explain that the oxidation reaction occurs below 40 °C and further reaction (α-arylation reaction) of the oxidation product occurs when the temperature was raised to 80 °C. We would like to point out here that the ability of a ketone to be enolized (formation of 2H), but not benzaldehyde, contributes to the fact that the secondary alcohols can further react with one more equivalent of PhCl leading to α-arylation.
Experimentally, selective oxidation of secondary alcohols versus primary alcohols has been observed. Consistent with this experimental observation, the oxidation reaction of the secondary alcohol PhEtCHOH was calculated to be kinetically more favorable than that of the primary alcohol PhCH2OH. Our detailed mechanistic analysis indicates that the overall reaction barrier for the oxidation of a given alcohol substrate molecule is the sum of (1) the barrier for oxidative addition of phenyl chloride to the active species LPd(0)(alcohol) and (2) the reaction free energy for the ligand substitution of the alcohol substrate molecule for the aldehyde or ketone product molecule in LPd(0)(aldehyde) or LPd(0)(ketone) to regenerate the active species LPd(0)(alcohol). We found that a ketone product molecule binds much more weakly to the Pd(0) metal centre than an aldehyde product molecule, explaining that oxidation of secondary alcohols, but not primary alcohols, occurs under the catalytic conditions.
The Pd-catalysed domino oxidation–arylation reactions of secondary alcohols have also been investigated. Our calculation results support that the secondary alcohols were first oxidized by phenyl chloride to their corresponding ketones according the mechanism summarized above and then further reacted with the second equivalent of the oxidant phenyl chloride at a higher temperature to finally produce the α-arylated products. In the α-arylation reaction of the ketone product molecule, we have established that the mechanism involves (1) oxidative addition of PhCl to LPd(0)(ketone), (2) ligand substitution of tBuO− for Cl−, (3) deprotonation of an α-H of the coordinated ketone molecule by the tBuO− ligand, (4) reductive elimination to give a species which contains the α-arylated product molecule as a ligand, and (5) finally a ligand substitution of a ketone molecule for the α-arylated product molecule to regenerate the active species LPd(0)(ketone). The overall barrier calculated for the α-arylation of the oxidation product PhEtCO was found to be higher than the oxidation of the secondary alcohol PhEtCHOH. These results explain the experimental observation that α-ketone arylation occurs when the temperature was raised to 80 °C from below 40 °C.
As a final note to readers, we would like to point out that the conclusions summarized above are based on the assumption that LPd(0)(alcohol) is the active species for the catalytic cycle. In the presence of tBuOK, one would expect that the anionic species [LPd(0)(OtBu)]− is likely the active species. However, our calculations show that the anionic species [LPd(0)(OtBu)]− gives an overall reaction barrier of 17.3 kcal mol−1 for the catalysed oxidation of PhCH2OH, which is too small to account for the oxidation reactions observed experimentally (slow and with elevated temperature). We argue that the poor solubility of the tertbutoxide and the solvent used (toluene) prevents the fast formation of the anionic species. We also examined whether decarbonylation is responsible for the oxidation reaction of primary alcohol PhCH2OH not being observed. Our calculation results do not support this. We hope that this work will invite more research on these reactions to gain more insights into the reaction mechanism.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c4qo00214h |
This journal is © the Partner Organisations 2014 |