Miguel
Rivero-Crespo‡
*a,
Judit
Oliver-Meseguer
a,
Klaudia
Kapłońska
b,
Piotr
Kuśtrowski
b,
Emilio
Pardo
c,
José Pedro
Cerón-Carrasco
*d and
Antonio
Leyva-Pérez
*a
aInstituto de Tecnología Química, Universitat Politècnica de València–Consejo Superior de Investigaciones Científicas, Avda. de los Naranjos s/n, 46022, Valencia, Spain. E-mail: miguel_angel.rivero_crespo@org.chem.ethz.ch; anleyva@itq.upv.es; Fax: +34963877809; Tel: +34963877800
bFaculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Krakow, Poland
cInstituto de Ciencia Molecular (ICMol), Universidad de Valencia, 46980 Paterna, Valencia, Spain
dReconocimiento y Encapsulación Molecular (REM), Universidad Católica de Murcia (UCAM), Spain. E-mail: jpceron@ucam.edu
First published on 6th July 2020
The Pt-catalysed addition of silanes to functional groups such as alkenes, alkynes, carbonyls and alcohols, i.e. the hydrosilylation reaction, is a fundamental transformation in industrial and academic chemistry, often claimed as the most important application of Pt catalysts in solution. However, the exact nature of the Pt active species and its mechanism of action is not well understood yet, particularly regarding regioselectivity. Here, experimental and computational studies together with an ad hoc graphical method show that the hydroaddition of alkynes proceeds through Pt–Si–H clusters of 3–5 atoms (metal(oid) association) in parts per million amounts (ppm), which decrease the energy of the transition state and direct the regioselectivity of the reaction. Based on these findings, new extremely-active (ppm) microporous solid catalysts for the hydrosilylation of alkynes, alkenes and alcohols have been developed, paving the way for more environmentally-benign industrial applications.
Heterogeneous catalysis is considered a suitable solution when soluble catalysts are toxic and not recyclable, and the more direct and simple approach to heterogenise a soluble catalyst is by supporting it on a simple solid.4 However, examples of truly heterogeneous Pt-supported solid catalysts for hydrosilylation reactions, which mimic the extremely high catalytic activity of Pt in solution, are difficult to find in the literature.5 The main reason for this shortage of reported examples is that the exact nature of the catalytically active Pt species formed in solution during reactions, as well as the structure of key intermediates during the different hydrosilylation reactions, remains mainly unknown, beyond the postulated Chalk–Harrod intermediate shown in Fig. 1.6 The extremely low amount of Pt required to catalyse the hydrosilylation of alkenes and alkynes in solution is well below the detection limits for many experimental techniques, which has precluded conclusive studies on Pt active species. If these catalytically active Pt reaction intermediates are properly detected, the design of new solid catalysts would be easier and more efficient.
Here, we show experimental and computational studies, together with the development of an ad hoc graphical method, which support the association of silane molecules with Pt atoms in solution, prior to or during the rate determining step (rds) of the reaction, to generate 3–5 atom cyclic metal(oid) clusters that facilitate the reaction by providing an alternative low-energy pathway to the classical Chalk–Harrod mechanism. It is also shown here how these cyclic metal(oid) clusters control the electronics and sterics during the alkyne hydrosilylation and govern the final chemo-, regio- and stereoselectivity of the reaction, and with this information in hand, new and extremely active soluble and solid-supported Pt catalysts can be designed not only for the hydrosilylation of alkynes but also for the hydrosilylation of alkenes and alcohols.
The results in Fig. 2 [see also Fig. S1–S3 in the ESI†] show the rate equations for the hydrosilylation of phenylacetylene 1 with triethylsilane 2 catalysed by either Pt1 or Pt3, obtained from the initial kinetic points, which are r0 = kexp[Pt1][1]−1[2]2 and r0 = kexp[Pt3][1]−1[2], respectively. In other words, the hydrosilylation reaction rate is, in both cases, linearly proportional to the amount of Pt and inversely proportional to the amount of alkyne but, unexpectedly, it is quadratically proportional to the amount of silane for Pt1 but not for Pt3. For the latter, saturation kinetics is observed at high silane concentrations ([2] > 0.50 M). It is noteworthy that the main product for Pt1 is the typical β-(E) alkenylsilane,15 with >95% selectivity, while for Pt3 is the less common α-alkenylsilane,9,16 with >90% selectivity. However, here the formation of both products is equally considered to assess the catalytic activity of the different Pt catalysts.
A Hammett plot for different para-substituted dimethylphenylsilanes prepared by a reported method17 gives a similar negative slope for both catalysts (ρPt1 = −0.50, ρPt3 = −0.60, Fig. S4†), and the kinetic isotope effect (KIE) values for the hydrosilylation of 1 with H(D)SiEt3 are 2.0 and 1.7 ± 0.2 for Pt1 and Pt3 respectively (see further below for a more precise value through the graphical method), which confirm that a positive charge is being built and that the Si–H(D) bond is breaking before or during the rate determining step of the hydrosilylation reaction, for both Pt1 and Pt3 catalysts respectively. However, this catalytic similarity between Pt1 and Pt3 is clearly broken when considering the enthalpy and entropy values in the transition state (ΔH‡ and ΔS‡). Table 1 shows the ΔH‡ and ΔS‡ values for the hydrosilylation reaction of 1 with 2, calculated experimentally with the Eyring–Polanyi equation, and the results show that the Pt1 catalyst significantly decreases both the enthalpy and entropy values (∼50 kJ mol−1 and −190 J K−1 mol−1, respectively, entry 1) relative to Pt3 (∼110 kJ mol−1 and −25 J K−1 mol−1, entry 2), which raises the apparent contradiction that the transition state is more associated for Pt1 than for Pt3. To assess the influence of the alkyne and the silane on the ΔH‡ and ΔS‡ values, the study was expanded to 1-octyne 3 and 1,1,3,3-tetramethyldisiloxane (TMDS) 4, and the results (entries 3–6) show that a different alkyne does not modify the energetics of the transition state while a different silane indeed does. These results fit with the different reaction order found with respect to the Pt atomicity for the silane but not for the alkyne and suggest that not only the nature but also the number of silane molecules is key in the transition state of the hydrosilylation reaction, at least with Pt1.
Entry | Alkyne/silane | Pt | ΔH‡ (kJ mol−1) | ΔS‡ (J K−1 mol−1) |
---|---|---|---|---|
1 | 1/2 | Pt1 | 53 ± 4 (52 ± 2) | −175 ± 10 (−156 ± 6) |
2 | Pt3 | 114 ± 6 (109 ± 2) | −2 ± 16 (5 ± 6) | |
3 | 3/2 | Pt1 | 52 ± 3 (54 ± 3) | −171 ± 8 (−168 ± 15) |
4 | Pt3 | 103 ± 12 (123 ± 4) | 22 ± 33 (40 ± 17) | |
5 | 1/4 | Pt1 | 24 ± 3 (38 ± 2) | −252 ± 9 (−210 ± 20) |
6 | Pt3 | 28 ± 3 (41 ± 2) | −242 ± 6 (−188 ± 5) |
In order to confirm these energetic differences, a more accurate graphical method was developed, since the initial rate approximation encounters some limitations in reactions with induction times, catalyst deactivation, late steady state attainment and very fast rates,18 which can be the case here. In this context, Reaction Progress Kinetic Analysis (RPKA)19 and Variable Time Normalization Analysis (VTNA)20 have emerged in recent years as powerful tools to measure reaction orders from catalysed reactions with high precision in few experiments.21 However, both RPKA and VTNA do not cover other kinetic and thermodynamic parameters such as the transition state enthalpy and entropy, linear free energy relationships (LFER) and KIE (Fig. S5†).22 Despite temperature scanning methodologies having been recently developed to calculate energetic activation parameters,22d this methodology is much more complex in terms of instrumentation and data treatment. For all these reasons, and inspired by the simple graphical methodologies mentioned above, a graphical method to calculate the ΔH‡, ΔS‡, Hammett parameter (ρ) and KIE values, using the tools available in our laboratory, was developed here. This new methodology does not only allow the calculation of the energetic values with high precision in few experiments, but also proposes a way to calculate the errors associated with the measurements (see the ESI† for discussion).
Fig. 3 illustrates the basis of our approach for a simulated bimolecular reaction A + B → P, i.e. the hydrosilylation reaction, with the calculation of the transition state enthalpy and entropy values (ΔH‡ and ΔS‡, see the ESI† for the mathematical model development). The final methodology is very simple and involves representing the product concentration vs. the normalized time-scale f or g(t, T, ΔH‡ or ΔS‡). In this way, the transition state enthalpy ΔH‡ or entropy value ΔS‡ is easily obtained by iterations with any available user graphic software (Excel®, Origin®, …) until all kinetic profiles overlay at different temperatures, which can be done either by just visually, adjusting to a polynomic equation until obtaining the best fit in terms of R2 or by minimizing the Euclidean distance between the data points. The latter approach was used in this work (see the ESI†). This methodology gives the activation parameters directly from raw concentration data, avoiding the propagation of errors associated with the initial rates, which can be as high as ±8% at R2 = 0.997 for the classical method (Fig. S6 and Tables S1–S5†), but only up to ±4% at nearly perfect fittings (R2 > 0.999) for the new graphical method. Noteworthily, the graphical method overrides the experimental distortions arising from induction times, since these induction times are included in the whole kinetic profile, which was validated for a simulated (Fig. S7, S8, Tables S6 and S7†) and a reported related reaction23 (Fig. S9†).
The results of ΔH‡ and ΔS‡ obtained by the graphical method are shown in Table 1 (see also Fig. S10–S15†) and reinforce the conclusions obtained with the traditional kinetic values, i.e. the transition stateis less energetic (lower ΔH‡) and more associated (lower ΔS‡) for Pt1 than for Pt3. Fig. 4 shows that the graphical method can also be applied to find the Hammett parameter (ρ, Fig. S16 and Tables S8–S10†) and KIE values (see the ESI† for discussion), which gives ρ = 0.50 for the Pt1 catalyst and KIE = 1.91(5) for the Pt3 catalyst, exactly the same value as that for Pt1 but with much more precision (see above and Table S11†). These results confirm the energetic and isotopic kinetic data obtained for the hydrosilylation of alkynes with Pt1 and Pt3, and support the hypothesis of metal(oid) Pt–silane cluster formation during the transition state of the reaction.
The occupancy of the free coordinating sites in TS-1 by alkynes is supported by the negative dependence of the hydrosilylation rate on the alkyne concentration, for both Pt1 and Pt3 catalysts, which indicates that the saturation of the coordinative sites of Pt with alkyne molecules occurs at the very beginning of the catalytic cycle, in accordance with the alkynophilic nature of the metal. This hypothesis is supported by the enhanced reaction rate found for electron-rich alkynes (worse π-acceptors), following the order: electron-deficient phenylacetylenes < electron-rich phenylacetylenes < alkyl alkynes, assessed by the negative slope in the Hammett plot (Fig. S18†), and the decrease in the reaction rate when phenylacetylene 1 is added during the hydrosilylation of 1-octyne 3, but not the reverse, in accordance with the stronger π-accepting ability of 1 (Fig. S19†), and the lack of the kinetic isotope effect (KIE) found for deuterated phenylacetylene (1.1 ± 0.1 for both Pt1 and Pt3), which indicates that σ-Pt alkynylide species do not participate in the reaction. Thus, the incoming silane must displace coordinating alkyne molecules in the inner sphere of Pt through σ-interactions,27 to give the key cyclic metalloid cluster.
A four-membered metal(oid) cyclic transition state may also be involved for Pt3 clusters, in this case with just one silane molecule, shown in Fig. 6 as TS-2 (see also Fig. 1). The linearity of the rate equation with [2] (see above) supports this intermediate. Furthermore, taking advantage of the different regioselectivity found for Pt3 and the fact that the regioselectivity comes from the electronic effect produced by alkyne polarization, the reactivity of different internal alkynes where the two carbon atoms are electronically similar, was tested. Fig. 6 also shows that no α/β differentiation occurs between Pt1 and Pt3 catalysts while, in contrast, an extremely polarized internal alkyne such as trimethylsilylphenylacetylene 7, where the charge difference between the two carbon atoms is comparable to a terminal alkyne (Table S12†), produces a change in the α/β isomer ratio from 0.37 (Pt1) to 0.60 (Pt3). These results support that the electronics of the alkyne is determinant for the regioselectivity during the Pt3-catalysed hydrosilylation reaction, but not with the Pt1 catalyst.28 Regarding sterics, electronically similar ortho- and para-substituted phenylacetylenes show that the former react less with Pt1 clusters than with Pt3 (Fig. 6, S20 and S21†), and these results make sense considering that the cyclic metal(oid) transition state around Pt1 is more congested than around Pt3, thus enhancing steric effects during the reaction (Fig. S21†).29 The so-called ortho effect, initially described for the hydrostannylation of alkynes30 and associated with higher differences in 13C NMR for internal and external carbon atoms, cannot be invoked here since it is caused not only by electronic but also by steric effects.31 Two-dimensional nuclear magnetic resonance nuclear Overhauser effect spectroscopy (NOESY) of the α-vinylsilane product, deuterated or not, shows that both the anti- and Markovnikov hydrosilylation reactions with Pt1 and Pt3 proceed completely in a syn-fashion (Fig. S22–S24†), which is in accordance with the observed products and previous studies on silylplatination and reductive elimination steps.32
According to the calculations, the proposed TS-1 metal(oid) cluster decreases the barrier for the hydride addition (intermediates [A]) by ca. 10 kcal mol−1 (black and blue lines) with respect to the standard Chalk–Harrod intermediate TS-3 (red and green lines), regardless of the regiochemistry of the product. This stabilization comes from the significantly less constrained environment around the Pt atom when two silane molecules bind and cooperate, compared to the canonical Chalk–Harrod configuration. Furthermore, the smallest energy barrier for the final silylation step (intermediates [B]) corresponds to the formation of the β-isomer through the TS-1 intermediate (black line), and thus the DFT calculations predict the right product after a silane self-assisted hydroplatination process. All these accumulated theoretical outcomes back up the hypothesis of cyclic metal(oid) cluster formation.
![]() | ||
Fig. 8 (a) Chalk–Harrod mechanism (n = 1) and the alternative mechanism found in this work (n = 1 and 3). (b) Transition states for the oxidative addition proposed in this study. |
Taking into consideration the elementary steps of the mechanism proposed in Fig. 8, validated not only by experimental results but also by computational results and the graphical method, a theoretical kinetic equation for the Markovnikov hydrosilylation of alkynes catalysed by Pt1 can be obtained, using the steady state approximation (Fig. S25† for mathematical development). With the appropriate simplifications, the resulting calculated equation resembles the experimentally obtained equation (vide supra), with the quadratic [silane]2 dependence:
Therefore, the mechanism in Fig. 8 is supported not only by kinetic, isotopic, thermodynamic and computational data, but also by the theoretical rate equations.
Both Pt1/NaY and Pt3/NaY, in ppm, give the α-vinylsilane product, but while Pt1/NaY showed a long induction time (∼25 minutes) and no leaching of Pt into solution, as assessed by the hot filtration test (Fig. 9A) and ICP analysis (<1% of supported Pt into solution), Pt3/NaY showed, in contrast, a shorter induction and leaching of inactive Pt species into the solution (∼10% measured by ICP), which evolve to active species after 50 minutes (Fig. 9B). In a similar way to Pt3 clusters being formed in situ in solution from simple homogeneous precursors,9,16 we reasoned here that Pt atoms inside Pt1/NaY could also evolve to form Pt3 clusters under reaction conditions, remaining within the zeolitic framework, while the pre-formed Pt3 clusters inside the zeolite Pt3/NaY are less stabilized by the framework and tend to leach out, which converts the Pt3/NaY just to a simple reservoir of Pt.
Fig. 9A also shows that the Pt1/NaY catalyst could be reused without any induction time (up to five times), without loss of activity and with just a small fraction of the confined Pt3 aggregates being converted into nanoparticles (<2 nm),40 according to high resolution transmission electron microscopy (HR-TEM) measurements (Fig. 9, bottom). These results strongly support that the active Pt3 species are formed in situ and stabilized within the zeolite framework, ready to be further reused without leaching. In accordance with the leaching results, the preformed Pt3/NaY sample loses activity and still shows some induction time (Fig. 9B).
Fig. 10 shows the results for the hydrosilylation reaction of different alkynes and silanes catalysed by Pt1/NaY, and it can be seen that just 1 mg of the zeolite per mmol of alkyne (50 ppm of Pt) is enough to achieve full conversion and good selectivity towards the corresponding α-vinylsilanes, including aromatic (products 8a–c–10 and 12–13) and alkyl alkynes (products 11a,b and 14–15), aromatic (products 8b,c, 11b, 12 and 15) and alkyl silanes (products 8a–11a and 14–14), and also an internal alkyne (product 14). Monosusbtituted phenylsilane PhSiH3 reacts selectively once with an aromatic alkyne to give the corresponding silane 8c in reasonable 51% yield (the rest is β isomer and not polyalkylated silanes), and it does not react with an alkyl alkyne (product 11c). The silane HSi(OEt)3 did not yield the desired product either, in any case, but rather decomposed. The reaction can also be run at room temperature, provided that more solid catalyst (20 mg per mmol of alkyne) is used, and this result indicates that the active and zeolite-stabilized Pt3 species can be formed inside the zeolite without thermal requirements. To our knowledge, this is the first efficient and recyclable solid catalyst for the Markovnikov hydrosilylation of alkynes.41 Previous tests9 showed that Pt/C does not really work as a catalyst to obtain the alpha product but just generate species in solution with low alpha selectivity. Besides, not only the selectivity but also the activity is considerably higher in Pt/zeolites than in Pt/C. At this stage, and in order to gain more information about the Pt species on the solid during the reaction, a 195Pt solid NMR experiment was envisaged; however, the acquisition time to have a reliable result was estimated to be unacceptably long.
Fig. 11 shows the single crystal X-ray diffractogram of Pt1/MOF, where the single Pt atoms are marked with an arrow.38b As can be seen, the MOF is constituted by ∼1 nm channels with homogeneously distributed isolated Pt atoms stabilized on the framework walls by electronic interactions.35,38b,d This structure resembles that expected in the cavities and channels of the zeolite and, thus, it is not surprising that the Pt1-supported microporous MOF catalyses the hydrosilylation of alkenes similarly or even better than Pt1/NaY does. Fig. 11 also shows the scope for the hydrosilylation of terminal alkenes catalysed by Pt1/MOF, and not only aromatic (products 18–23) but also alkyl alkenes (products 24–25) engage well in the reaction. Different electron poor silanes are reactive, including PhSiH3 as a selective monosilylation partner; however, triethylsilane 2 was not suitable in this case.
A hot filtration test for the hydrosilylation of styrene 16 with the Pt1/MOF catalyst shows that the ∼30% of the catalytic activity comes from active species in solution (Fig. S27,† left); however, most of the catalytic activity lies in the solid under reaction conditions. Furthermore, when one looks at the kinetic profile closely, a short but clear induction time of 5 min can be seen (Fig. S28,† right), an induction time that is not present at all when the hydrosilylation of 16 is catalyzed even by tiny amounts of Karstedt's catalyst in solution, according to our experiments (Fig. S29†). The induction time with Pt1/MOF may indicate that initial Pt1 on the MOF is not the active species and evolves to active Pt1 during the reaction. Coulombic forces are behind the stabilization of Pt on the negatively-charge walls of the microporous solids. Thus subtle electronic modifications of the Pt1 atom by donation of electron density from the solid framework to the metal atom could modulate and ultimately increase the catalytic activity.45 Indeed, by just heating Pt1/NaY in an oven at 200 °C for 1 hour before the reaction, the hydrosilylation of 16 proceeds with full conversion and with the highest TOF measured (entry 9 in Table 2). X-ray photoelectron spectroscopy (XPS) in combination with CO-probe low temperature infrared (IR) measurements of the Pt1/NaY sample calcined at 200 °C show that the amount of Pt(0) increases at least twice with respect to Pt(II) in non-calcined Pt1/NaY and Pt1/MOF (Fig. S29†). These results support the formation of catalytically active, reduced Pt1 species within the microporous catalyst. It should be noted that the Pt1 annotation used here does not necessarily presuppose hardly reduced species, since discussion focuses on Pt atomicity rather than on the exact valence state.
Entry | Pt catalyst | Yield (%) | TOF0 (h−1) |
---|---|---|---|
1 | Karstedt | 96 | 126 |
2 | Pt(acac)2 | 26 | 30 |
3 | Pt(NH3)4Cl2·xH2O | 33 | 37 |
4 | Pt(COD)(Me)2 | 38 | 43 |
5 | Pt/C | 87 | 99 |
6 | PtCl2/SiO2 | 37 | 42 |
7 | Pt1/MOF | 35 | 42 |
8 | Pt1/NaY | 93 | 106 |
9 | Pt 1 /NaY (calc. at 200 °C) | 99 | 118 |
10 | Pt1/NaY (calc. at 300 °C) | 64 | 76 |
11 | Pt1/HY | 61 | 72 |
12 | Pt1/CsY | 19 | 22 |
13 | Pt3/NaY | 22 | 26 |
Karstedt's catalyst is the only efficient soluble Pt catalyst from those tested (entries 1–4), and considering that Karstedt's catalyst only contains Pt(0) atoms, and that Pt/C is also quite effective (entry 5), it seems that Pt(0) could be the active redox species for the dehydrogenative silylation of alcohols. This hypothesis fits well with the studies above on alkynes and alkenes and also with the lower activity shown by PtCl2/SiO2 and Pt1/MOF (entries 6 and 7) with respect to Pt1/NaY (entry 8), which is able to easily modulate the Pt site by the negatively-charged zeolite framework. Indeed, Pt1/NaY calcined at 200 °C, which tends to reduce the Pt atoms further,38a gives a 99% yield of 27 and the highest initial TOF (entry 9). Calcination at higher temperatures to provoke metal aggregation (entry 10),52 and decreasing or increasing the electronic density of the zeolite with a different charge compensating cation to Na+ (H+ and Cs+, entries 11 and 12 respectively),38a do not improve the catalytic activity and supports the need of fine tuning the Pt site on truly heterogeneous microporous solid catalysts for the hydrosilylation reaction, when leaching does not occur (Fig. S30†). These facts point again to an in situ formation of the Pt catalyst, and having in mind the H-atom bridges found in the metal(oid) cluster intermediates and the decisive effect of the zeolite confined space, it is possible to propose a tentative mechanism for the dehydrogenative silylation of alcohols based on metal(oid) cycles (Fig. S31†).53 Although speculative, this mechanism should be taking into consideration for future catalyst design based on Pt for the dehydrogenative hydrosilylation of alcohols.
Fig. 12 shows that the Pt1/NaY solid calcined at 200 °C for 1 hour, catalyses the dehydrogenative hydrosilylation of alkyl alcohols (products 28–32, 34 and 36), benzyl alcohols (product 33) and phenol derivatives (products 34 and 37–38) in good yields and with just 100 ppm of Pt, and with different silanes. Interestingly, the disilane H2SiPh2 yields only the monosilylated product without further reaction, and the silane HSi(OEt)3 was also reactive in this case, to give products 30 and 37, respectively, in very good yields. To our knowledge, this is the first Pt-based efficient solid catalyst for the dehydrogenative hydrosilylation of alcohols.
Footnotes |
† Electronic supplementary information (ESI) available: General procedures, additional tables, figures and movies, and compound characterization. See DOI: 10.1039/d0sc02391d |
‡ Current address: Department of Chemistry and Applied Biosciences, Swiss Federal Institute of Technology (ETH) Zürich, Vladimir-Prelog-Weg 1-5/10, 8093 Zürich, Switzerland. |
This journal is © The Royal Society of Chemistry 2020 |