Petrus C. M.
Laan
a,
Martijn J.
Mekkering
a,
Felix J.
de Zwart
a,
Alessandro
Troglia
b,
Roland
Bliem
b,
Kai
Zhao
c,
Norbert J.
Geels
a,
Bas
de Bruin
a,
Gadi
Rothenberg
*a,
Joost N. H.
Reek
a and
Ning
Yan
*ac
aVan't Hoff Institute for Molecular Sciences, University of Amsterdam, Science Park 904, 1098XH Amsterdam, The Netherlands. E-mail: g.rothenberg@uva.nl
bAdvanced Research Center for Nanolithography (ARCNL), Science Park 106, 1098XG Amsterdam, The Netherlands
cKey Laboratory of Artificial Micro- and Nano-Structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, 430072, China. E-mail: ning.yan@whu.edu.cn
First published on 14th August 2024
The local coordination environment of single atom catalysts (SACs) often determines their catalytic performance. To understand these metal–support interactions, we prepared Pt SACs on cerium dioxide (CeO2) cubes, octahedra and rods, with well-structured exposed crystal facets. The CeO2 crystals were characterized by SEM, TEM, pXRD, and N2 sorption, confirming the shape-selective synthesis, identical bulk structure, and variations in specific surface area, respectively. EPR, XPS, TEM and XANES measurements showed differences in the oxygen vacancy density following the trend rods > octahedra > cubes. AC-HAADF-STEM, XPS and CO-DRIFTS measurements confirmed the presence of only single Pt2+ sites, with different surface platinum surface concentrations. We then compared the performance of the three catalysts in ammonia borane hydrolysis. Precise monitoring of reaction kinetics between 30–80 °C gave Arrhenius plots with hundreds of data points. All plots showed a clear inflection point, the temperature of which (rods > octahedra > cubes) correlates to the energy barrier of ammonia borane diffusion to the Pt sites. These activity differences reflect variations in the – facet dependent – degree of stabilization of intermediates by surface oxygen lone pairs and surface–metal binding strength. Our results show how choosing the right macroscopic support shape can give control over single atom catalysed reactions on the microscopic scale.
Controlling these EMSIs is crucial. Cerium dioxide (CeO2) is known for having strong EMSIs promoting high metal dispersity as well as tailoring catalytic performances for nanoparticle-based catalysts.5,6 These strong interactions are mainly controlled by surface oxygen vacancies which differ per crystal facet. CeO2 crystals can be synthesized as cubes, octahedra or rods, and each shape features a unique terminal facet: (100), (111) and (110), respectively. Controlling the shape of the CeO2 support is thus a straightforward way for tailoring these EMSIs and has led to developments in sensing, solar cells and catalysis.7 There are indeed various reports in the literature showing the shape (facet) effects of ceria supported catalysts. Yet, understanding such effects remains challenging for typical heterogeneous catalysts that consists out of supported transition metal nanoparticles. On one hand, only the metal atoms close to the support are subject to strong EMSIs whereas reactions eventually occur across the entire surface of the nanoparticle. On the other hand, the control of facet effects is often accompanied by the change of other physico-chemical parameters of ceria (e.g., defect density, surface area) which will affect the catalytic performance as well. Therefore, disentangling the real facet effects from various influential factors remains ambiguous.8,9
In this context, we prepared platinum SACs on CeO2 cubes, octahedra and rods with distinct surface facets and studied their catalytic performance in the solvent phase, an environment which is rarely probed on the study of EMSIs using SACs.10–13 Using aberration corrected scanning transmission electron microscopy (AC-STEM) and diffuse reflectance infrared Fourier transform spectroscopy of adsorbed CO (CO-DRIFTS), we examined the surface defect structure of ceria on various facets and confirmed the atomic distribution of Pt. A deeper understanding of the defect chemistry was obtained by investigating X-ray absorption near edge structure (XANES) and electron paramagnetic resonance spectroscopy (EPR). In the ammonia borane hydrolysis carried out at various temperatures, all catalysts showed a kink in their Arrhenius plots at different temperatures which reflects the surface travel of ammonia borane to the active sites. Our results show that the shape of the support particles influences the stability of Pt sites and reaction intermediates, and thereby the catalytic activity.
Most importantly, we probed the density of oxygen vacancies by CW X-band electron paramagnetic resonance (EPR) spectroscopy at room temperature (Fig. 1c).18,19 The typical Ce3+–O−–Ce4+ (g = 1.963) defect was observed in all samples as indicated in yellow.20–23 What differs among the samples is the signal intensity, which follows the order rod > octa > cube. While exact quantification remains difficult, this shows that the density of oxygen defects is directly correlated to the terminal crystal facet type.24
Next, we ran aberration-corrected HAADF-STEM measurements to study the structure and dispersion of Pt on the three supports (Fig. 2). These show the known interplanar spacings for CeO2 cubes (100), octahedra (111) and rods (110) of 0.27, 0.31 and 0.19 nm, respectively, confirming the expected terminal facets.14–16 Moreover, they show the atomic dispersion of Pt on the surfaces of all three catalysts without the formation of any clusters or nanoparticles.
We then studied the oxidation state and surface concentration of platinum using X-ray photoelectron spectroscopy (XPS) and diffuse reflectance infrared Fourier transform spectroscopy of adsorbed CO (CO-DRIFTS, Fig. 3). The XP spectra of the Pt 4f core levels were all dominated by one component with minor core level shifts among the samples (<0.6 eV), indicating a single and highly similar type of platinum species close to a 2+ oxidation state for all catalysts (Fig. 3a).26,27 The CO-DRIFTS was used for evaluating the electronic structure and dispersion of Pt (Fig. 3b).28 All samples show a symmetrical vibration band around 2105 cm−1 which is characteristic of linear bound CO on isolated Pt2+ atoms.29–31 The lack of an additional band around 2030 cm−1 indicates the absence of linearly bound CO on Pt0, confirming that all of the platinum on the catalyst surface is Pt2+, in agreement with the XPS data. More importantly, the absence of additional bands around 1950 cm−1 (CO adsorbed on the interface between Pt clusters and support) and 1860 cm−1 (bridge bound CO on two Pt atoms) strongly suggests that Pt was atomically dispersed on all the samples and not agglomerated into larger clusters, in line with AC-HAADF-STEM-EDX data. The identical platinum oxidation state of all three samples (Pt2+) allows us to study purely the facet effect.32 Next, we determined the amount of Pt atoms present at the surface by XPS (Table 1).
The relative metal content was proportional to their corresponding CO-DRIFTS peak areas (Fig. S9–S11†). These amounts do depend on the crystals' shape, following the trend cube > octa > rod.33 This cross-correlation between two independent measurement techniques is important, as it confirms that both techniques provide information solely on the surface Pt species. Based on these values and the respective SSAs of the supports, we estimated the average distances between two Pt atoms (Table 1). The Pt1–Pt1 distance is significantly shorter for Pt1@cube and Pt1@octa compared to Pt1@rod. However, it is unlikely that their catalytic behavior will be influenced by their neighbouring SACs – this is generally observed at metal-to-metal distances below 1.2 nm only.34 This is further supported by the fact that we only see Pt2+ species at the surface in all samples (vide supra). Thus, the catalysts differ in only two aspects: their density of oxygen defects and their density of Pt2+ sites.
NH3BH3 + 4H2O → NH4+ + B(OH)4− + 3H2 | (1) |
The rate of the reaction over this temperature range was determined by quantifying the amount of hydrogen evolved in a single experiment by applying a temperature ramp of 2 °C using a home-built bubble counter.36 This simple and inexpensive setup measures gas evolution with volume steps of ca. 10 μL, generating highly precise Arrhenius plots based on hundreds of data points. Therefore, we can study even subtle physico-chemical changes of the catalysts. The supports themselves were not active in the reaction (Table S3†) and rate data is corrected for the amount of Pt present at the surface (see Arrhenius plots in Fig. 4).
All three catalysts show low- and high-temperature regions separated by a clear inflection point. Before this kink, the activation energies and pre-exponential factors are relatively low, indicating moderate rates. Thereafter the rates increase significantly (Fig. S12†). Such a kink in the Arrhenius plot often reflects a change in reaction mechanism. However, when we ran additional isothermal experiments at temperatures before and after the kink we observed similar kinetic profiles, only at different rates (Fig. S13 and S14†). These reaction profiles are characteristic of a zero-order rate in ammonia borane, indicating that O–H bond cleavage in water is the rate-determining step.37 Interestingly, this kink is not observed for ceria-supported Pt nanoparticles (cf. Fig. S15†), affirming that our catalysts are SACs rather than nanoparticles.
Based on this, we maintain that water activation is the rate-determining step between 30–80 °C. The kinks in the Arrhenius plots are not caused by a change in mechanism, but rather by temperature-dependent surface diffusion of ammonia borane to the Pt active sites. Because water is both reactant and solvent, the active sites are nearly always occupied by water molecules (the [H2O]:
[NH3BH3] ratio is ca. 550
:
1). Yet as both substrates are required for a successful reaction, the rate depends on the availability of ammonia borane at the active site, and therefore on the travel of ammonia borane to this site.37 This surface travel depends on the respective surface–adsorbate binding constants and on the distance to be travelled (i.e. the active site density and distribution). Surface oxygen defects are known to stabilize such adsorbates, increasing the energy barrier for surface travel.38 Given the trends in oxygen defect density (rod > octa > cube) and active site density (rod < octa < cube), the energy barrier for ammonia borane traveling to the active site should follow the trend rod > octa > cube. Indeed, we see this in our Arrhenius plots: with increasing temperatures, the cubes show the kink first (43 °C) followed by the octahedra (52 °C) and the rods (60 °C).
To further test this hypothesis, we increased the NH3BH3:
catalyst ratio to see if this would indeed enhance the ammonia borane diffusion to the Pt sites. We did this by running a non-isothermal ammonia borane experiment at 50% of the catalyst loading while keeping the ammonia borane concentration identical and using Pt1@octa as an example. We see that the kink temperature decreased from 52 °C to 38 °C (Fig. S16†). This indicates that ammonia borane diffusion to the Pt sites has a lower barrier at higher relative substrate concentrations which is in line with the observed trend for the different facets.
Besides the fact that the Arrhenius plots showed kinks, the slopes also differ significantly (Fig. 4). This means that the activation energies for ammonia borane hydrolysis on Pt SACs depend on the facet they are supported on. This is partially caused by the differences in support–ammonia borane binding (vide supra), yet there are more parameters governing catalyst performance which are facet dependent. Firstly, Pt SACs bind stronger to CeO2 supports with more oxygen defects.10,11 Thus, based on the EPR, XPS and XANES data, the surface–metal atom binding strength should be rod > octa > cube (Fig. 1c and S6†). This means that Pt atoms supported on cubes are available for catalysis at lower temperatures compared to octahedra and rods (Fig. 4 and S12†). Secondly, reaction intermediates can be stabilized by surface oxygen atom lone pairs on ceria based catalysts.39–41 Reactions are therefore faster when the active site is surrounded by oxygen atoms rather than by surface defects – the most common surface defect is an oxygen vacancy in the surface lattice.7 This further supports the rod < octa < cube activity trend.
The recyclability of the Pt1@cube catalyst in ammonia borane hydrolysis was also tested at 65 °C for three consecutive runs (Fig. S17†). The maximum conversion and initial rate decrease in each run, likely due to catalyst poisoning by metaborate salts (commonly seen in borohydride hydrolysis42). Deactivation due to Pt leaching is less likely,43 as one would expect activity of the leached species. This was affirmed by SEM-imaging showing plate-like aggregates on the surface of the catalyst material (Fig. S18,†cf. with Fig. S1a† that shows the pristine cube-shaped crystals).
N2 adsorption–desorption isotherms were measured on a Thermo Scientific Surfer instrument at 77 K, using vacuum-dried samples. More specific, around 100 mg of each sample was dried at 100 °C for 16 h on a Belprep-vacIII prior to analysis. The specific surface area was determined based on the adsorption branch and the BET analysis was performed according to the Rouquerol consistency criteria (Fig. S2–S4†).44,45
Powder X-ray diffraction (pXRD) patterns were obtained with a Rigaku MiniFlex II diffractometer (Tokyo, Japan) using Ni-filtered CuKα radiation (λ = 1.541874 Å) at 30 kV and 15 mA. For each measurement, the sample was ground and loaded on a monocrystalline silicon sample holder with an 8 mm wide and 0.2 mm deep cavity. The powdered sample was pressed firmly in the cavity to make a uniform flat sample area. Residual sample outside the sample cavity was removed to minimize background scattering. Diffraction patterns were collected between the 2θ range of 20° and 90° using a rotation speed of 2° min−1, a step size of 0.05° and 1 s dwell time.
CW X-band electron paramagnetic resonance (EPR) spectra of the samples were measured in EPR quartz tubes on a Bruker EMX-plus CW X-band spectrometer at room temperature. The effective g values were defined as the magnetic field strength at the maximum microwave absorption according to eqn (S1) in which geff is the effective g-value, h is Planck's constant being 4.135 × 10−15 eV s, v is the microwave frequency of the spectrometer being 9.643 GHz, μB is the Bohr magneton being 5.788 × 10−5 eV T−1 and B is the applied magnetic field at the maximum microwave absorption maximum in T (1 Gauss = 1 × 10−4 T).
![]() | (S1) |
Metal loadings were determined based on inductively coupled plasma optical emission spectrometry (ICP-OES) analysis performed by the Mikroanalytisches Laboratorium Kolbe, Oberhausen, Germany. Samples were prepared using microwave digestion and then analysed with a Spectro Arcos analyzer of Spectro capable of maintaining a standard error of ±1.5 ppm.
Aberration-corrected high angle annular dark-field scanning transmission electron microscopy (AC-HAADF-STEM) measurements, elemental mappings, and energy-dispersive X-ray (EDX) spectroscopy measurements were taken on a JEOL JEM-ARM300F2 GRAND ARM™ 2 instrument coupled with a high-angle annular dark field (HAADF) detector and an energy-dispersive X-ray spectroscopy detector. The set-up was operated at 300 kV and delivers a spatial resolution of ≤60 pm in both TEM and STEM resulting in atomic-resolution imaging of the samples. Intensity profiles and interplanar spacings were determined using the Gwyddion 2.61 software package.
Ce M4,5-edge and O K-edge X-ray absorption near edge structure (XANES) spectra were obtained at the X-ray Magnetic Circular Dichroism (XMCD) beamline of Hefei Light Source (HLS). After baseline correction, the Ce M4,5-edge data was fitted by standard Gaussian curves using peak fitting in Origin 2018.
X-ray photoelectron spectroscopy (XPS) was performed in ultra-high vacuum (base pressure below 2 × 10−9 mbar) using a Scienta Omicron HiPP-3 analyzer with a 1 mm entrance slit operated in Swift Acceleration mode and a monochromatic Al Kα source. Survey spectra were acquired at a pass energy of 500 eV and the high-resolution spectra were acquired at pass energies of 100 eV (Pt1@cube and Pt1@octa) or 300 eV (Pt1@rod). Prior to data processing, the binding energies were calibrated using that of adventitious carbon (C 1s at 284.8 eV). The Ce 3d and O 1s regions are fitted using peak positions and ratios from literature.46 XPS peak fitting was performed using KolXPD (Kolibrik), employing a (ranged) Shirley background and Voigt functions for the individual components.
Diffuse reflectance infrared Fourier-transform spectroscopy (DRIFTS) measurements were performed on a Nicolet iS50 FTIR (Thermo Fisher, United States) spectrometer equipped with a MCTA detector and a KBr beam splitter. The diffuse reflection accessory used was a Praying Mantis™ (DRP, Harrick, United States) in combination with a high temperature reaction chamber with ZnSe windows (HVC-DRM-5, Harrick, United States) and a temperature control unit (ATK-024-4, Harrick, United States). Typically, ∼50 mg KBr was put in the reaction chamber which was topped with 10 mg of finely grounded sample. This was placed in the in situ cell and pre-treated at 200 °C under vacuum (typically 5 × 10−4 mbar, Pfeiffer HiCube) for 30 minutes followed by a pure O2 (99.999%) flow for 30 minutes to remove any other gases, adsorbed water and other impurities. After cooling to room temperature, the cell was vacuumed again, and a background spectrum was collected. Then, the cell was saturated with CO (99.9%). After 1 min, the cell was vacuumized (5 × 10−4 mbar) and the spectra were recorded in the range between 4000 and 800 cm−1 with a spectral resolution of 4.0 cm−1 until no changes between the collected spectra were visible. Each spectrum was recorded by averaging 32 scans. The peak areas were determined by peak fitting using a standard Gaussian curve using Origin 2018 (Fig. S8–S10†).
For non-isothermal ammonia borane hydrolysis, a screwcap vial (10.0 mL) was charged with a stirring bean (8.0 × 3.0 mm), catalyst (7.0 mg) and water (8.0 mL). The catalyst was homogeneously suspended with the aid of an ultrasonic bath for ten minutes at room temperature. The screwcap vial was mounted on the reactor head. A freshly prepared ammonia borane solution (0.4 mL, 2.0 M) was loaded into syringe (1.0 mL) equipped with a glass capillary (∅ = 0.32 mm, l = 15 cm). The capillary was inserted through one of the syringe ports into the reaction mixture and purged with nitrogen (5.0 mL min−1) while the reaction mixture was cooled down to 5.0 °C by an external ice bath. When this temperature was reached, the reactor was placed in the heating mantle under stirring and the other three remaining syringe ports were closed off. After five seconds, the ammonia borane solution was injected and a ramp of 2.0 °C min−1 was initiated achieving a controlled temperature increase between 30 and 80 °C. The sample was held at 80.0 °C until no gas production was observed anymore.
For isothermal ammonia borane hydrolysis, a screwcap vial (10.0 mL) was charged with a stirring bean (8.0 × 3.0 mm), catalyst (7.0 mg) and water (8.0 mL). The catalyst was homogeneously suspended with the aid of an ultrasonic bath for ten minutes at room temperature. The screwcap vial was mounted on the reactor head. A freshly prepared ammonia borane solution (0.4 mL, 0.2 M) was loaded into syringe (1.0 mL) equipped with a glass capillary (∅ = 0.32 mm, l = 15 cm). The capillary was inserted through one of the syringe ports into the reaction mixture and purged with nitrogen (5.0 mL min−1) while pre-heating the reaction mixture to the desired temperature under stirring by the use of a heating mantle. When stabilized at the desired temperature, the other three remaining syringe ports were closed off and the ammonia borane solution was injected. Gas production was monitored until reaction completion keeping the temperature of the reaction mixture at the desired temperature.
For the recyclability tests, the reaction mixture was centrifuged at 4000 rpm for 10 min after reaction completion in its original screwcap vial, giving a white residue and a colorless supernatant. The supernatant was removed with a needle and replaced by fresh water to a total volume of 8.0 mL. Thereafter, the purified specimen was resuspended by ultrasonication for 15 min. This washing procedure was repeated three times before a subsequent catalytic reaction was started. After the last run, the supernatant was removed, and the remaining white powder was further dried under vacuum (5 mbar at 50 °C) for 16 h. The remaining catalyst was weighed, confirming that there were no losses during the washing procedures.
Footnote |
† Electronic supplementary information (ESI) available: SEM imaging, nitrogen physisorption studies, additional calculations, XPS spectra and fittings, XANES spectra and fittings, CO-DRIFTS spectra and fittings and additional kinetics studies. See DOI: https://doi.org/10.1039/d4cy00484a |
This journal is © The Royal Society of Chemistry 2024 |