Jose
Pinto
ab,
Andreas
Weilhard
*a,
Luke T.
Norman
a,
Rhys W.
Lodge
a,
David M.
Rogers
a,
Aitor
Gual
c,
Israel
Cano
d,
Andrei N.
Khlobystov
a,
Peter
Licence
*b and
Jesum
Alves Fernandes
*a
aSchool of Chemistry, University of Nottingham, NG7 2RD, Nottingham, UK. E-mail: Andreas.Weilhard1@nottingham.ac.uk; Jesum.AlvesFernandes@nottingham.ac.uk
bGSK Carbon Neutral Laboratories for Sustainable Chemistry, University of Nottingham, NG7 2TU, Nottingham, UK. E-mail: Peter.Licence@nottingham.ac.uk
cCTQC-Eurecat-UTQ, C/Marcel·lí Domingo 2, Building N5, 43007 Tarragona, Spain
dDepartment of Inorganic Chemistry, Faculty of Chemistry, Universidad Complutense de Madrid, 28040 Madrid, Spain
First published on 2nd June 2023
In this work, we demonstrate that the synergistic effect of PdAu nanoparticles (NPs) in hydrogenation reactions is not only related to high activity but also to their stability when compared to Pd mono-metallic NPs. To demonstrate this, a series of mono- and bi-metallic NPs: Pd, Pd0.75Au0.25, Pd0.5Au0.5, Pd0.25Au0.75 and Au in ionic liquid [C4C1Im][NTf2] have been fabricated via a magnetron sputtering process. Bi-metallic NPs possess external shells enriched with Pd atoms that interact with [NTf2]− of the ionic liquid resulting in enhanced catalytic performance in hydrogenation of cinnamaldehyde compared to their mono-metallic counterparts. This is ascribed to their higher stability over 24 h of reaction, whilst the catalytic activity and selectivity are comparable for both catalysts. Using a bespoke kinetic model for in situ catalyst deactivation investigations and electron microscopy imaging at the nanoscale, we have shown that PdAu has a deactivation rate constant of 0.13 h−1, compared to 0.33 h−1 for Pd NPs, leaving 60% and 40% of available sites after the reaction, respectively. Beyond that, the kinetic model demonstrates that the reaction product has a strong stabilizing factor for bimetallic NPs against coarsening and deactivation, which is not the case for Pd NPs. In summary, our kinetic model enables the evaluation of the catalyst performance over the entire chemical reaction space, probing the contribution of each individual component of the reaction mixture and allowing the design of high-performance catalysts.
An excellent model reaction for the investigation of catalyst stability, selectivity and activity is the selective hydrogenation of cinnamaldehyde (CAL) to hydrocinnamaldehyde (HCAL). HCAL is employed as an important building block for the synthesis of important pharmaceuticals;12 however, the hydrogenation of CAL can result in a wide range of products, including cinnamyl alcohol (COL) and hydrocinnamyl alcohol (HCOL).13 Due to the competing selectivity, extensive research has been undertaken to obtain an active and selective catalyst for this transformation. In general, Ni14 and Pd15,16 based catalysts are found to be more active and selective towards HCAL compared to Pt,17 Ru,18,19 and Au,16 although a substantial effect can also be attributed to the employed supports and ligands.20–23 Another important aspect is the synergistic effect observed in bi-metallic systems leading to an increase of the activity and selectivity for this reaction either by facilitating the adsorption and activation of reactants and therefore enhancing the reaction rate,24–26 or by forming electrophilic sites, thus tuning the selectivity.27,28 In this regard, PdAu nanoparticles (NPs) represent an attractive catalyst as the H2 activation can be facilitated due to the synergistic effects of Pd and Au species.24–26 Moreover, PdAu NPs have found widespread application not only in the selective hydrogenation of CAL,29–31 but also in other catalytic hydrogenations32–34 and oxidations.35–39 However, mechanisms of mono- and bi-metallic catalyst deactivation remain largely unknown due to the lack of in situ deactivation kinetic measurements combined with local-probe analysis (e.g. electron microscopy) of the catalyst evolution.
In this work we developed a kinetic model describing mono- and bimetallic catalyst deactivation, by correlating transmission electron microscopy (TEM) imaging of catalytic Pd and PdAu NPs with reaction kinetics measurements. Direct growth of Pd and PdAu NPs in [C4C1Im][NTf2] ionic liquid (IL) via the magnetron sputtering approach provided 2.0–2.8 nm highly active catalytic centres with a core–shell structure for the bi-metallic NPs where Pd atoms are on the surface due to interactions with [NTf2]−. In the reaction of hydrogenation of CAL, the activity and selectivity of Pd0.5Au0.5 and Pd NPs were found to be similar but with stability significantly higher in the bi-metallic case due to a slower decrease of their available sites linked to a slower rate of coarsening of Pd0.5Au0.5 NPs as compared to Pd NPs. Interactions of the reaction product (HCAL) with the catalyst surface has been identified as the key factor for bi-metallic catalyst stabilisation during the reaction.
Fig. 1 The magnetron sputtering process, and structural and chemical characterisation of Pd and PdAu NPs. (a) Scheme of Pd and Au magnetron sputtering co-deposition in [C4C1Im][NTf2]. (b) Transmission electron microscopy (TEM) imaging of Pd0.5Au0.5 NPs in [C4C1Im][NTf2] shows the homogeneous NP size distribution generated by the magnetron sputtering process. A complete set of TEM images is available in the ESI† (Fig. S5). (c) X-ray photoelectron spectroscopy (XPS) spectra of Pd, Au and PdAu NPs in [C4C1Im][NTf2]. (d) Pd K-edge and Au L-edge X-ray absorption near edge structure (XANES). (e) Pd K-edge Fourier transform of k2-weighted extended X-ray absorption fine structure (EXAFS) spectra of PdAu NPs. (f) Schematic showing the formation of PdAu NPs with Au-core and Pd-shell as the proportion of Pd increases. |
Bright-field TEM imaging revealed a similar average NP size for Pd and Pd0.75Au0.25 and Pd0.5Au0.5 NPs, from ca. 1.98 to 2.18 nm, with a slightly larger NP size increase for high Au NP content Pd0.25Au0.75 and Au NPs from ca. 2.35 to 2.84 nm. This behaviour can be associated with the weaker interaction of [NTf2]− with Au atoms compared to Pd atoms, thus leading to a larger Au NP average size.40–42
X-ray photoelectron spectroscopy (XPS) and X-ray absorption spectroscopy (XAS) provided insights into the interactions between Pd–Au atoms within NPs and metal atoms with ionic liquid. The XPS spectra of monometallic Pd and Au NPs in [C4C1Im][NTf2] show two distinct electronic environments: the metallic Pd 3d5/2 at 336.1 eV and Au 4f7/2 at 84.3 eV ascribed to atoms in a metallic environment in the NP core, and another set of peaks at 338.4 eV for Pd 3d5/2 and 85.7 eV for Au 4f7/2 corresponding to surface atoms interacting with ionic liquid (Fig. 1c).43 For PdAu NPs in [C4C1Im][NTf2], XPS detects significant changes in the electronic environment of the metals as the Pd:Au atomic ratio is varied (Fig. 1c). As the Pd concentration increases in the Au matrix (Au → Pd0.25Au0.75 → Pd0.5Au0.5), the intensity of the peak associated with surface Au atoms interacting with ionic liquid (ca. 85.7 eV) decreases and completely disappears for Pd0.75Au0.25. Additionally, the Au peak at 84.3 eV, attributed to Au NP core metal atoms, showed a significant shift towards a lower binding energy as the Pd concentration increased (ca. 0.5 eV shift for Pd0.75Au0.25), while no significant changes were observed for Pd peaks.44 These XPS measurements demonstrate the interaction between Pd and Au atoms within NPs, thus confirming the formation of bi-metallic PdAu NPs, as well as suggesting that the surface of NP is enriched with Pd atoms.
X-ray absorption near edge structure (XANES) analysis for Pd and Au NPs displayed a similar adsorption edge to metal foil but with a higher white line (Hw) intensity, which can be assigned to the Pd and Au surface atoms interacting with the IL (Fig. 1d).43 For PdAu NPs, the Hw intensity of the Pd K-edge is higher than that for Pd NPs, whereas for the Au L-edge, the Hw intensity peak of PdAu NPs is lower than that of Au NPs, which corroborates the XPS results.44 Extended X-ray absorption fine structure (EXAFS) measurements of Pd NPs show a peak at 1.5 Å associated with the Pd–IL interaction and another peak at 2.5 Å ascribed to the Pd–Pd bond (Fig. 1e). As the Au concentration increases in the Pd matrix, a peak for longer bond distances is observed at 2.7 Å, which is attributed to the Pd–Au bond.43 Additionally, in PdAu NPs the peak intensity associated with the Pd–Pd bond (ca. 2.5 Å) is dramatically decreased with respect to the Pd–IL peak (ca. 1.5 Å), which indicates the tendency of Pd atoms to be on the NP surface in the presence of Au. Interestingly, the formation of Pd shell/Au core NPs contrasts the thermodynamically most stable alloy configuration for these bimetallic nanoparticles.45 However, it has been shown that the interaction with strongly coordinating molecules can lead to the segregation of Pd from the nanoparticle surface.45 Similarly, the interaction with [NTf2]− could lead to the segregation of Pd atoms from the NP surface, especially considering that during the initial metal NP formation in IL by magnetron sputtering dynamic clusters are formed, as we recently demonstrated.43 The dynamic metal clusters are most likely to be thermodynamically unstable and kinetically reactive, thus small differences between Pd and Au in coordination strengths with [NTf2]− may lead to the segregation of Pd atoms to the NP surface. Furthermore, XPS and XAS data are consistent with a core–shell model in which Au atoms are preferentially located in the NP core and Pd atoms on the NP surface (Fig. 1f).46,47
Pd, Au and PdAu NPs in [C4C1Im][NTf2] were tested in the hydrogenation of CAL. Under the optimised reaction conditions, all NPs containing Pd showed high selectivity towards HCAL (>90%, Fig. 2a). Our results are comparable to PdAu NPs supported on high surface area mesoporous silica48 (however, it must be noted that the catalyst support and reaction conditions may have a significant influence on this reaction).29–31,49 Furthermore, Pd0.75Au0.25 and Pd0.5Au0.5 NPs exhibited higher conversions after 24 h than Pd and Pd0.25Au0.75 NPs (Fig. 2a, entries 3, 2 and 1, 4, 5).24–26 Au NPs were inactive under the assayed reaction conditions, indicating that Pd is crucial for the catalytic turnover in this transformation (Fig. 2a, entries 4, 5).
A more detailed analysis of the reaction profiles and rates provided insight into the superior catalytic performance of PdAu NPs compared to the Pd catalyst, both in terms of activity and stability (Fig. 2b and c). For example, for Pd NPs the reaction rate drops and thus the conversion has only a small increase after reaching ca. 40% conversion, whereas no such effect was observed for Pd0.5Au0.5 and Pd0.25Au0.75 NPs (Fig. 2b). Pd0.75Au0.25 NPs displayed the highest initial rate, achieving 70% conversion of CAL to HCAL after 8 h, but thereafter the conversion increased only slightly (Fig. 2b). Even though Pd0.75Au0.25 NPs exhibit much higher initial rates than Pd0.25Au0.75, Pd0.5Au0.5 and Pd NPs, the highest conversion after 24 h was found for Pd0.5Au0.5 NPs with a continuous increase in conversion over time (Fig. 2c). These results show that the correlation of the initial reaction rate (related to catalyst activity) with the overall conversion (related to catalyst activity and stability) is not always correct, thus the overall catalytic performance cannot be solely attributed to an enhancement of the initial reaction rates,24–26 and must include catalyst stability over an extended period of time as well as any changes in the reaction mechanism caused by the different catalysts.24 As these differences are most pronounced with Pd and Pd0.5Au0.5 NPs, these two catalysts were investigated in greater detail. TEM analysis of the catalysts after reaction was performed to determine the catalyst changes occurring during the reaction (Fig. 3).
TEM analysis of Pd NPs showed a significant increase in NP size after the reaction with the mean diameter changing from 2.0 to 5.46 nm, whereas the size distribution of Pd0.5Au0.5 NPs changed only slightly with the mean diameter before and after the reaction from 2.2 to 3.3 nm, respectively (Fig. 3). Significant coarsening of Pd NPs during the reaction leads to reduction of the surface area available for catalysis, whilst Pd0.5Au0.5 NPs changing only slightly under the same reaction conditions and retain much of their surface area. Indeed, a more detailed analysis of the size distribution shows that for Pd0.5Au0.5 NPs 61% of surface sites after the reaction remain available for catalysis, whilst in Pd NPs only 36% of initially active sites remain available. This suggests that the long-term catalyst stability (rather than activity) is the main factor discriminating the performance of Pd and Pd0.5Au0.5 NPs.
To reveal how the compounds of the reaction mixture (THF, H2, CAL, HCAL) affect the catalyst stability, Pd and Pd0.5Au0.5 NPs were exposed to the individual components for an incubation period (τ) prior to the reaction initiation (Fig. 4).
Fig. 4 Dependence of the relative concentration of total available sites [ST] on the incubation time τ for (a) Pd NPs and (b) Pd0.5Au0.5 NPs in the presence of only THF (black), in the presence of 0.5 M HCAL (green), in the presence of 0.5 M CAL (red), and under 5 bar H2 (blue). In all experiments, it was assumed that the catalyst deactivation takes place only during the incubation period, thus the initial rates are not affected by further catalyst deactivation. The relative concentration of available sites was estimated using the fraction of the initial rate after the incubation period and the initial rate at τ = 0 (see ESI† for further details). |
For the first set of experiments, Pd and Pd0.5Au0.5 NPs were subjected to incubation periods, without addition of reagents/products (only THF solvent), prior to initiating the reaction (Fig. 4). This led to an exponential decay of the initial reaction rates as a function of pre-reaction incubation time, suggesting a first order dependence for the catalyst deactivation (Fig. 4, S7 and S8†) measured as the total concentration of available sites [ST]. Comparing Pd and Pd0.5Au0.5 NPs, a much faster deactivation rate for Pd NPs was observed, which is consistent with extensive Pd NP coarsening shown by TEM analysis (Fig. 3). The presence of reagents or product(s) in the pre-reaction exposure of catalysts retarded deactivation (e.g. increased their stability) for both catalysts, which can be ascribed to catalyst surface stabilisation due to adsorption of reagents or products (Fig. 4). For example, the presence of H2 considerably slowed the deactivation of Pd NPs during the incubation period as compared to pure solvent, whereas the addition of CAL and HCAL had a negligible effect on Pd NPs. For Pd0.5Au0.5 NPs, under the same conditions, the presence of HCAL had the greatest effect on the catalyst during the incubation period, with effects of H2 and CAL being less pronounced. It is important to highlight that these results were obtained using the standard reaction conditions assayed in this work (Fig. 3a).
The experimental results presented in Fig. 4 strongly indicate a pseudo-first order mechanism for the catalyst deactivation due to the particle coarsening (hence leading to reduced surface area), which is retarded by both reagent and product adsorption on the catalysts. By considering the quasi-equilibria taking place during the incubation periods, eqn (1) can describe the concentration of available sites ([ST]τ) as a function of time τ (see ESI† for further details), where [ST]0 is the total concentration of initial available sites, k5 is the rate constant for the catalyst deactivation, and K1 and K4 are the equilibrium constants for the adsorption of CAL and HCAL on the catalyst surface, respectively. The adsorption equilibrium constant for H2 (K2) is multiplied by the Henry constant to give Kp which allows the use of the hydrogen pressure P(H2) in the equation:
(1) |
As expected, Pd0.5Au0.5 NPs exhibit a significantly lower rate constant for the catalyst deactivation than Pd NPs (Table 1). For Pd NPs, CAL showed a stronger effect on their stability compared to H2 (Table 1), which was an unexpected outcome based on the results shown in Fig. 4 explained by the much higher concentration of H2 than that of CAL in the reaction conditions (Fig. 4), whereas in the quasi-equilibrium constant calculation, the concentrations of reagents and the product were normalised (Table 1). Furthermore, these results were confirmed by DFT calculations (see details in the ESI†), which showed that the computed Gibbs free energy for adsorption of CAL and H2 adsorbed on a Pd19 cluster are −36.76 and −4.87 kcal mol−1, respectively, in line with the experimental quasi equilibrium constants (Table 1). For Pd0.5Au0.5 NPs, HCAL had the largest adsorption equilibrium constant (Table 1), following a similar trend as that observed under the reaction conditions (Fig. 4), which underpins the strong effect of HCAL on Pd0.5Au0.5 NP stabilisation.
Constant | Pd NPs | Pd0.5Au0.5 NPs |
---|---|---|
Catalyst deactivation rate constant (k5/h−1) | 0.33 | 0.13 |
CAL adsorption equilibrium constant (K1/M−1) | 2.75 | 0.54 |
H2 adsorption equilibrium constant (Kp/bar−1) | 1.30 | 0.70 |
HCAL adsorption equilibrium constant (K4/M−1) | 0.25 | 6.36 |
Using the kinetic parameters in Table 1, the relative number of catalyst available sites was calculated and correlated with TEM results (Fig. 4) providing deep insights into the mechanism of the catalyst deactivation process (Fig. 5).
Fig. 5 Calculated relative numbers of the catalyst available sites after the reaction using kinetic parameters and correlating them with TEM imaging. (a and b) Calculated kinetic traces of Pd and Pd0.5Au0.5 NPs for the loss of their available sites in different excess experiments, respectively (for original traces see Fig. S10†). (c) Comparison of Pd NPs and Pd0.5Au0.5 NPs for estimated available sites after the reaction determined by TEM and kinetic traces. |
Remarkably the loss of available sites calculated by the kinetic analysis (Fig. 5) for both Pd and Pd0.5Au0.5 NPs agrees with the TEM observations for the spent catalysts. For instance, Pd NPs increase in size by 3.5 nm during the reaction, whilst Pd0.5Au0.5 NPs increase only by 1.3 nm under the same conditions indicating 36% and 61% of remaining available sites, respectively, which correlates well with the kinetic data, emphasising the key difference in stability of mono- and bi-metallic catalysts (Fig. 5c).
In order to generate a complete kinetic model involving the catalytic transformation, the apparent rate orders were investigated using different excess experiments and analysed via the variable time normalization analysis (VTNA) method (detailed discussion in the ESI,† Fig. S10–S13).51–54 Therefore, with both the kinetic rate orders of reagents and product(s) (Fig. S11†) and the individual parameters influencing the catalyst deactivation in hand (Fig. 5), a complete picture of the catalytic transformation was established. As shown in Fig. 6a, the hydrogenation proceeds via L–H type transformation of CAL to HCAL, whilst the catalyst deactivates due to particle growth which is described as a pseudo first order. The relative number of active sites is utilised in conjunction with the rate orders determined from kinetics measurements (Fig. S10–S13†) to produce a complete normalisation of the kinetic traces (Fig. 6b and c). Note that the rate order for CAL changes for Pd NPs from 1 to −1, whilst Pd0.5Au0.5 NPs remaining +1 throughout the assayed reaction conditions. Furthermore, the unusual rate orders found for Pd NPs (Fig. S10–S13†) can be ascribed to the effects of the individual components (e.g. solvent, reactants and products) on the catalyst deactivation as only the order of +1 and −1 are present in the total normalisation of the kinetic traces (Fig. 6b). It is important to emphasise that it is the combination of complete normalisation of the kinetic traces (Fig. 6b and c) with TEM imaging of the catalyst that allows us to validate the kinetic model for in situ analysis of the catalyst deactivation process.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3cy00289f |
This journal is © The Royal Society of Chemistry 2023 |