Evolution of Pd species over γ-Al2O3 and defective Al2O3-x during CH4 oxidationThe synthesis of defective alumina was realized by pyrolysis of aluminum sulfate dissolved in a deep eutectic solvent composed of urea and thiourea (Fig. 1a, see Supplementary information, Section 1 for details). Pd/Al2O3 and Pd/Al2O3-x catalysts with a nominal Pd loading of 1 wt%, were prepared by the incipient wetness impregnation (IWI) method followed by calcination at 400 °C for 2 h (see Supplementary information, Sections 1.3 and 4 for details). The 27Al solid-state nuclear magnetic resonance (ssNMR) results revealed the presence of three prominent peaks at around 5, 38 and 70 ppm over Al2O3-x (Fig. 1b), which can be ascribed to Al3+ species in tetrahedral (AlO4 at 2–8 ppm), pentahedral (AlO5 at 36–40 ppm) and octahedral (AlO6 at 68-70 ppm) coordination, respectively19,21. To quantify the different Al species, the 27Al ssNMR spectra were fitted by three components using a simple Czjzek model22,23. After loading of Pd, a decrease in the mole percent content of both AlV (from 40 to 34%) and AlIV (from 26 to 24%) on the Al2O3-x support was observed, accompanied by an elevation in the content of AlVI. While previous research has evidenced the role of AlV in the anchoring of isolated noble metals19, the observed decrease in the relative content of AlIV species here may also indicate binding of Pd atoms on AlIV. One possible explanation is the occurrence of localized surface reconstruction caused by the interaction between distorted AlIV species and Pd through oxygen bridging (Fig. 1a), leading to the conversion of tetrahedrally coordinated Al to an octahedral coordination structure. Furthermore, the total quantity of consumed AlV and AlIV sites was eight times greater than the amount of loaded Pd, which may be attributed to: (a) the fixation of Pd atoms requiring the consumption of multiple AlV and AlIV species24 (Fig. 1a and Supplementary Information, Section 8.1 for details); (b) the introduction of Pd causing surface restructuring of amorphous alumina due to strong metal-support interaction19,21; and (c) the presence of trace amounts of water in the impregnation solution resulting in hydroxylation of AlV and distorted AlIV sites. Notably, the remaining AlV was still sufficient to immobilize an equal or even greater amount of Pd species or to restrict the mobility of Pd during the reaction25. In contrast, 27Al NMR spectroscopy did not detect AlV sites on γ-Al2O3 at room temperature, regardless of the presence of Pd species (Fig. 1c). The proportion between AlIV and AlVI also remained invariant on γ-Al2O3 before and after the loading process.Fig. 1: Stabilizing active palladium species on defective alumina sites.a Schematic illustration of the synthesis of defective alumina (Al2O3-x) and Pd/Al2O3-x. 27Al solid state nuclear magnetic resonance (ssNMR) spectra and respective fitting for Al2O3-x (b bottom), Pd/Al2O3-x (b top), γ-Al2O3 (c bottom) and Pd/Al2O3 (c top) The deconvoluted peaks, denoted by blue, yellow, and red colors, respectively represent the tetrahedral, pentahedral, and octahedral coordination Al species. d Light-off curves for catalytic methane  (CH4) oxidation over fresh and aged Pd/Al2O3 and Pd/Al2O3-x catalysts under reaction conditions (500 °C, 1.5 h). Reaction conditions: CH4:O2:He = 2:8:90, 50 mg catalyst, GHSV = 60,000 h−1. e Arrhenius plots of the fresh and aged Pd catalysts. The high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) images of Pd/Al2O3-fresh (f, g); Pd/Al2O3-aged (h–j); Pd/Al2O3-x-fresh (k) and Pd/Al2O3-x-aged (l–n). The red and blue dashed boxes represent respectively the regions where PdOx clusters/nanoparticles (NPs) and single atoms are present on the catalysts. Pd clusters are marked with the red rings, while Pd single atoms are marked with blue rings.The catalytic activity of Pd/Al2O3 and Pd/Al2O3-x for the complete oxidation of CH4 was evaluated using a fixed-bed downflow reactor. The light-off curves (Fig. 1d) for CH4 conversion show that the fresh Pd/Al2O3-x exhibited higher catalytic activity (reaching 90% methane conversion at T90 = 367 °C) compared to Pd/Al2O3 (T90 = 429 °C). After ageing for 1.5 h at 500 °C with GHSV of 120,000 h−1, Pd/Al2O3-x-aged lost some of its low-temperature activity, but the catalytic activity was improved at relatively high temperatures (Fig. 1d). In contrast, Pd/Al2O3-aged showed noticeably lower activity compared to Pd/Al2O3-fresh across all temperatures. The corresponding Arrhenius plots (Fig. 1e) demonstrated that Pd/Al2O3-aged possesses a higher apparent activation energy (125 kJ/mol) in comparison to Pd/Al2O3-fresh (72 kJ/mol), Pd/Al2O3-x-fresh (64 kJ/mol) and Pd/Al2O3-x-aged (66 kJ/mol). Based on the apparent activation energy values, it can be inferred that the active centers over Pd/Al2O3-x exhibit higher intrinsic activity than those on Pd/Al2O3, and that they undergo negligible changes in the aging process. On the other hand, the ageing process significantly changed the active centers of Pd/Al2O3, generating intrinsically less active Pd species. The intercepts of the fitted Arrhenius plots, in first approximation related to the quantity of active sites, provide additional evidence that fresh Pd/Al2O3-x possesses a higher quantity of accessible active centers than Pd/Al2O3, whereas ageing causes the loss of active species on both catalysts26. Furthermore, during stability testing at 500 °C, Pd/Al2O3 gradually deactivated over time, consistent with previous studies10, while Pd/Al2O3-x remained relatively active for 15 h (Supplementary Information, Section 5 for details). Altogether, the results obtained suggest that Pd species on Al2O3-x were stabilized against thermally-induced deactivation processes during CH4 oxidation.We performed HAADF-STEM to obtain atomically resolved images of Pd species present on Pd/Al2O3 and Pd/Al2O3-x, for both fresh and aged catalysts. As shown in Fig. 1f, g, fresh Pd/Al2O3 consists mainly of small Pd NPs of 1–2 nm. After 1.5 h of ageing under reaction conditions, a fraction of the Pd species transformed into isolated Pd ions (Fig. 1j), potentially anchored on the generated AlV sites during reaction19. Meanwhile the remaining Pd species sintered and agglomerated into large PdO NPs (Fig. 1h, i). In contrast, the Pd on the fresh Pd/Al2O3-x catalyst consisted mainly of Pd single-atom and low-nuclear PdOx clusters (Fig. 1k). After ageing (Fig. 1l), a portion of the Pd species maintained their single-atom feature (Fig. 1m), while the other portion evolved into higher-nuclear Pd clusters (Fig. 1n), and no large NPs were observed (Supplementary Fig. 7), which we attribute to Pd anchoring to the more abundant AlV sites on the defective alumina support.In order to study the configuration of Pd species and their transformations during catalytic CH4 oxidation, we employed operando QEXAFS with a time resolution of 1s27,28,29. The catalysts were heated in reaction conditions and subsequently tested at 500 °C for 1.5 h, allowing us to track alterations in the Pd structure and compare these with catalyst deactivation. The normalized Pd K-edge X-ray absorption near-edge spectroscopy (XANES) and Fourier transformed EXAFS binned at intervals of 5 s to enhance the signal-to-noise ratio, are shown in Fig. 2a, b. Pd was predominantly in an oxidized state under reaction conditions. As the white line and the first minimum of the XANES are sensitive to a change in structure and oxidation state of the metal species, we plotted the normalized X-ray absorption of the white line (denoted as A for Pd/Al2O3 and C for Pd/Al2O3-x) and the first minimum (denoted as B for Pd/Al2O3 and D for Pd/Al2O3-x) against time on stream (Fig. 2c). This allowed us to compare the subtle changes in oxidation states of Pd species during the aging process. The white line absorption values of Pd species in Pd/Al2O3-x did not change appreciably, while for Pd/Al2O3 they increased slightly over time, suggesting a change in oxidation state of Pd species during the aging process. The R-space EXAFS indicated that the amount of metallic Pd species (at approximately 2.3 Å) on both catalysts remained roughly unchanged during this process (Fig. 2b), while the Pd-O contribution increased slightly. This can be attributed to growth of PdO NPs, consequently leading to an increase in Pd-O coordination number. It is worthwhile to mention that once steady state is reached under reaction conditions, the energies of the Pd K-edge absorption edge, the structures of the white line peak, and the absence of Pd-Pd scattering in the first shell of the EXAFS, all suggest the absence of metallic Pd as a composition of active centers. The on-line mass spectrometry (MS) results revealed a progressive decline in CO2 production over Pd/Al2O3 in the aging process (Fig. 2d). Consistent with the catalytic test results, these findings suggest a correlation between deactivation and the subsequent oxidation or sintering of Pd species. In contrast, the defective Al2O3-x support significantly suppressed the sintering/deactivation of Pd.Fig. 2: Correlation of Pd redox state and deactivation on different alumina supports.Operando X-ray absorption near-edge spectroscopy (XANES) (a) and R-space extended X-ray absorption fine structure (EXAFS) (b) of Pd/Al2O3 (top) and Pd/Al2O3-x (bottom) with 5 s time resolution obtained during aging in reaction conditions (CH4 : O2 : He = 2:8:90, reaction temperature: 500 °C, and GHSV = 9,554 h−1). c Plot of normalized X-ray absorption at the white line and first minimum (marked in Fig. 2a), as a function of time on stream. The corresponding on-line CO2 mass spectrometry (MS) signal collected during the operando X-ray absorption spectroscopy (XAS) measurements (d).Redox dynamics of Pd clusters and nanoparticles under CH4/O2 cyclesIn typical three-way catalytic converters, the exhaust composition undergoes fast perturbations between reducing (also known as fuel-rich) and oxidizing (lean) conditions, close the stoichiometric point30. This is known to cause changes in Pd speciation and can cause rapid aging of the catalyst. To mimic such rapid changes in redox conditions and study their effect on Pd species and catalytic activity, we have performed operando QEXAFS with 1 s time-resolution during alternating pulses (i.e., modulation excitation, in short ME) of CH4 and O2. While the spectra were recorded, the products were measured by online MS at 2.5 s time resolution. (See Supplementary Information, Section 7.4 for details). Time-resolved Pd K-edge XANES of Pd/Al2O3 and Pd/Al2O3-x at the white line region (24360–24400 eV) in a representative ME cycle are shown in Fig. 3a. Because Pd species are considered to follow the Mars-van Krevelen mechanism in CH4 oxidation reactions, i.e., the lattice O in PdO may be taken away or replenished in redox reactions, the change in the absorption intensity of the white line can quantitatively describe the change in the average valence of the Pd species8. Thus, from the white line feature of the time-resolved XANES, we can infer that the valence of the Pd species on Pd/Al2O3 changed more significantly than that on Pd/Al2O3-x in the ME cycle, and the Pd species on Pd/Al2O3-x maintain a higher oxidation state in CH4 atmosphere. In the ME experiments, the observed change in Pd species is mainly an interconversion between PdO and metallic Pd. To quantitatively characterize the changes in Pd oxidation state, we performed least squares linear combination fitting (LSLCF) of the Pd XANES, using PdO and metallic Pd as references. Considering that the XANES of atomically dispersed Pd species and PdO are very close (Supplementary Fig. 23), the fitted PdO fraction is the cumulative presence of single-atom Pd, PdOx clusters and NPs. Hence, we designate these components as Pd2+. The content of the Pd2+ and Pd0 species as a function of time in the ME experiment are shown in Fig. 3b. Interestingly, in the initial switch from O2 to CH4 atmosphere on Pd/Al2O3, around 90% of the Pd species underwent a transformation from Pd2+ to Pd0. However, in subsequent periodic pulsations alternating between CH4 and O2, a proportion of approximately 40% Pd0 and 10% Pd2+ was consistently maintained, while the remaining 50% of Pd species was oxidized or reduced depending on reaction atmosphere. The presence of irreducible Pd2+ species can be attributed to interfacial and single-atom Pd species, while Pd0 species that cannot be oxidized may possibly arise from Pd atoms located within the interior of the NPs. These atoms are generated during the initial reduction process but remain metallic throughout the oxidation process probably due to the significant energy barrier associated with the diffusion of O atoms from the surface to the interior31.Fig. 3: Redox dynamics of palladium species during simulated lean/rich cycles.a White line region of operando X-ray absorption near-edge spectroscopy (XANES) recorded at 1 s time resolution obtained upon switching between CH4 and O2. b Plots of Pd2+ and Pd0 relative contents over Pd/Al2O3 (top) and Pd/Al2O3-x (bottom) catalysts as a function of time, in the alternating CH4/O2 pulses. (Reaction temperature: 400 °C, and GHSV = 7,643 h−1). c Plot of online CO2 mass spectrometry (MS) signal as a function of time in a representative modulation excitation (ME) cycle. The variation in the content of Pd0 species over Pd/Al2O3 (d) and Pd/Al2O3-x (e) across the first to ninth reduction/oxidation cycles of the ME experiment, as a function of the relative CH4/O2 composition. The inserted bar chart illustrates the variation in the integrated area of the loops formed by the increase (during the reduction process) and decrease (during the oxidation process) of Pd0 species in each oxidation/reduction cycle, as a function of the cycle number. Error bars indicate the 95% confidence interval.In the case of Pd/Al2O3-x, on the other hand, approximately 40% Pd2+ and 10% Pd0 persisted unaffected by the oxidation/reduction atmosphere. This finding further suggests a higher proportion of interfacial/single-atom Pd species and smaller-sized PdOx clusters/nanoparticles over Pd/Al2O3-x compared to Pd/Al2O3. However, it is worth noting that the total quantity of Pd atoms capable of engaging in oxidation/reduction processes remained relatively similar between the two catalysts, accounting for approximately 50%. We propose that the remainder, non-redox-active species are located at the interface with alumina. Despite their comparable quantity, the redox-active Pd species on Pd/Al2O3 and Pd/Al2O3-x can be anticipated to exhibit different catalytic performance, because of their nature (more/less oxidic and dispersed).Figure 3c illustrates the plots of the online CO2 MS signal as a function of time in a representative ME cycle. Upon switching to a CH4 atmosphere, the CO2 signal on Pd/Al2O3-x exhibited an earlier onset and attained a higher peak compared to Pd/Al2O3. This potentially indicates a greater reducibility of PdOx clusters than NPs. Furthermore, on Pd/Al2O3-x there was a more pronounced decline in the CO2 signal subsequent to reaching its peak in comparison with Pd/Al2O3. This observation may suggest either a lower CO2 desorption energy on PdOx clusters, or a higher reduction energy barrier for the sub-surface Pd2+ on the PdOx NPs. The latter possibility can be ruled out as it does not align with the reduction rate of Pd2+ species indicated in Fig. 3b. Upon switching to the O2 atmosphere, CO2 MS signals were observed as well, albeit with an integrated area of approximately 30% compared to the previous process. The release of CO2 at this stage was attributed to the oxidation of residual CH4 from the previous gas pulse and/or to the oxidation of carbonaceous intermediates on the catalyst surface. The integration of CO2 MS signals during reduction in CH4 and oxidation in O2 reveals that the amount of CO2 produced on Pd/Al2O3-x is 1.45 and 2.92 times higher than on Pd/Al2O3, respectively. The potential reasons for the higher production of CO2 on Pd/Al2O3-x can may include: (a) the PdOx clusters are more active than NPs to catalyze CH4 oxidation; (b) the higher dispersion of PdOx clusters than NPs allows for the formation of a larger number of carbonaceous intermediates.Considering gas mixing and reaction during alternating pulses of CH4 and O2, the ratio of Pd0/Pd2+ species may exhibit dependence on the gas composition during the formed transients. Time-resolved QEXAFS and MS spectra can be combined in tracking the dependency of Pd species on the gaseous environment in the ME experiment. Since the combined content of Pd2+ and Pd0, as determined from the unconstrained LCF results, adds up to 1 with less than 0.06 error, our focus is specifically on plotting the changes in the Pd0 species content with respect to the relative concentrations of CH4 and O2 obtained from online MS, to gain a better understanding of the distinct responses of Pd atoms in clusters and NPs to the gas mixture interface during the reduction-oxidation cycles.The reduction of Pd2+ on Pd/Al2O3 proceeded slowly at lower methane concentrations, necessitating a specific threshold of the CH4/O2 ratio in the environment to expedite the reduction process. Conversely, during the oxidation process, Pd0 were oxidized upon the introduction of a small quantity of O2, with the oxidation rate gradually diminishing over time. In a single reduction/oxidation cycle, the reduction and oxidation curves of the Pd species interlace. As depicted in Fig. 3d, we denoted the two intertwined loops formed by the variations in the Pd0 species’ content with respect to CH4 concentration during the reduction-oxidation cycle as regions I and II. Region I is located at the higher relative CH4 concentrations, while region II at the lower relative methane concentrations. A larger area in region I indicates that a greater number of Pd2+ species can be reduced at lower CH4 concentrations, and more Pd0 species can be oxidized at lower relative O2 concentrations. Thus, the ratio of the areas I/II can serve as a descriptor of the reactivity of the Pd species. We integrated the areas of regions I and II in each ME cycle, and the changes in these integrated areas with cycle number are displayed in the inset of Fig. 3d. It is worth noting that the integrated area of region I increases with the number of cycles. This suggests that the ME cycling enhances the redox performance of Pd in Pd/Al2O3.On the other hand, Pd in Pd/Al2O3-x exhibited different responses to the alternating reaction gas pulses in the ME experiment. Specifically, the reducible Pd2+ species underwent immediate and complete reduction to Pd0 upon contact with CH4. Similarly, during the oxidation process, the oxidizable Pd0 species rapidly oxidize to Pd2+ upon exposure to O2. Considering that the relative concentrations of CH4 and O2 were obtained from online MS, it can be inferred that once Pd species on Pd/Al2O3-x were subjected to alternating CH4/O2 pulses, they preferentially consumed all CH4 or O2 present in the gas mixture interface through the reduction/oxidation of Pd species, or the catalytic CH4 oxidation. This process occurs at such a rapid rate that the apparent composition of Pd species appears to be independent of the relative concentrations of CH4/O2 in the gas mixture interface. Pd0 content as a function of CH4 concentration in the CH4/O2 pulse cycles forms a singular loop, region I in Fig. 3d, serving as a descriptor of the CH4 oxidation reactivity. Notably, the integrated area of region III shows minimal sensitivity to the number of ME cycles. As discussed earlier, the reducible/oxidizable Pd species on Pd/Al2O3 and Pd/Al2O3-x originate from non-interface and non-core Pd atoms located on PdOx NPs and PdOx clusters, respectively. The above analysis of their responses to CH4 and O2 pulses suggests that Pd atoms within PdOx clusters exhibit higher catalytic reactivity towards CH4 oxidation compared to PdOx NPs, with a lesser dependence on the relative concentration of CH4 or O2. With an increasing number of ME cycles, the activity of Pd atoms contained within PdOx NPs gradually improved, although it consistently remained lower than that of PdOx clusters. One possible explanation for this heightened activity is an increase in surface disorder or boundaries31,32,33. The gradual reduction in the amount of non-reducible Pd2+ species on Pd/Al2O3-x over time indicates the sintering of Pd in the ME experiment. However, the relatively stable area of region I of Fig. 3e implies that sintering does not significantly affect the activity of Pd atoms on the cluster surface.Unraveling the reaction mechanism of CH4 oxidation on Pd/Al2O3-x
Operando FT-IR spectroscopy was employed to further investigate the reaction mechanism of catalytic CH4 oxidation. In all cases, the FT-IR spectra of the catalysts in a He atmosphere at 400 °C were subtracted (see Supplementary Information, Section 6.1 for details), to visualize the adsorbed organic intermediates. After the aged Pd/Al2O3 and Pd/Al2O3-x catalysts were subjected to reaction conditions (400 °C, CH4: O2: He = 2:4:90, GHSV = 60,000 h−1) the FT-IR spectra were collected until reaching the steady states (Supplementary Information, Section 6.1 and Fig. 4a for details). The assignment of observed FT-IR absorption peaks to organic species is listed in Supplementary Table 6. A higher CH4 absorption band at the wavenumber of 1305 cm-1 was observed on Pd/Al2O3 catalyst, because of the lower catalytic activity of this catalyst. Nevertheless, the gradual accumulation of formate, carbonate, and bicarbonate species on the surface of Pd/Al2O3 was observed. The accumulation of carbonates and bicarbonates as reaction products on Pd/Al2O3 implies that under reaction conditions, Pd species may be poisoned by these species, leading to reduced activity. This is consistent with the observation of a lower rate of CO2 release over Pd/Al2O3 during ME (Fig. 3c). In contrast, when Pd/Al2O3-x is switched into the reaction atmosphere, the catalyst could maintain a relatively cleaner surface, with a certain amount of adsorbed bicarbonate and water, despite higher CH4 conversion.Fig. 4: Surface chemistry of palladium clusters and nanoparticles during CH4 oxidation.a Operando Fourier transform-infrared (FT-IR) spectra of Pd/Al2O3 (Pd NPs) and Pd/Al2O3-x (Pd clusters) at the steady states under reaction condition. CH4 : O2 : He = 2:8:90, reaction temperature = 400 °C, GHSV = 60,000 mL·gcat−1·h−1. In-situ FT-IR spectra of Pd/Al2O3 (b) and Pd/Al2O3-x (c) in the successive O2-to-CH4 and CH4-to-O2 switching processes, reaction temperature = 400 °C. Plots of the normalized integrated area of the IR absorption of each surface species on Pd/Al2O3 (d) and Pd/Al2O3-x (e) as a function of time during gas switching experiments.To record the build-up and degradation of the intermediates on the surface after excitation by different reactants, we performed operando FT-IR spectroscopy experiments with a time resolution of 1 min during alternating pulses of CH4 and O2 at 400 °C. The FT-IR spectra obtained over Pd/Al2O3 and Pd/Al2O3-x are shown in the Fig. 4b, c, respectively. To trace the occupation and exchange of different organic species with time, the FT-IR absorption peaks of all the observed species were integrated (Fig. 4d, e).When switching from O2 to CH4 atmosphere, two intermediates, *CO and *HCOO, gradually accumulated on the surface of Pd/Al2O3 and Pd/Al2O3-x. The vibrational frequencies of the adsorbed CO on both catalysts are less than 2100 cm−1, indicating that surface Pd species are metallic under these conditions34,35, which agrees with the operando XAS results (Fig. 3b). Furthermore, it is interesting to note that *CO and *HCOO feature different adsorption modes on Pd/Al2O3 and Pd/Al2O3-x, because of their different surface structures. Specifically, the absence of infrared absorption peaks in the wavenumber range of 1800–2000 cm−1 (bridge- and hollow- CO absorption bands) on Pd/Al2O3-x suggest that the active centers of Pd/Al2O3-x are composed of single atoms and small clusters during reaction, because only in that case would geometrical effects prevent the bridge- or hollow-adsorption of CO36. Besides, the linearly adsorbed CO on Pd/Al2O3-x, exhibits a higher vibrational frequency than that on Pd/Al2O3, implying that the strong interaction between Pd and Al2O3-x affords the atomically dispersed Pd species a partial positive charge37,38. This is consistent with the observation that the C-H vibration of the adsorbed formate species exhibits a slightly lower frequency on Pd/Al2O3-x (1364 cm−1) than on Pd/Al2O3 (1368 cm−1)39,40. From a time-dependence perspective (Fig. 3c), once switched from O2 to CH4 atmosphere, the hollow- and bridge-adsorbed CO covers the surface of Pd/Al2O3 preferentially amongst all intermediates. This can suggest the crucial role of the CO pathway on Pd/Al2O3. The coverage of the three CO intermediates: hollow-, bridge- and linear-CO, begin to decrease at minutes 8, 11 and 14, respectively. This is either due to the desorption of CO caused by a competition with other adsorbates, or because the CO species are further oxidized by subsurface oxygen. Once the CO intermediates are formed, carbonates begin to accumulate on the surface of Pd/Al2O3. After ten minutes of CH4 flow, when the coverage of carbonate reaches a specific threshold, the infrared absorption peaks for formate, bicarbonate and water become evident on Pd/Al2O3. Importantly, the accumulation of carbonate and bicarbonate on the surface of Pd NPs may induce alterations in the reaction pathways, favouring the formation of *HCOO. In contrast, although the *CO intermediate on Pd/Al2O3-x also appears early in the atmosphere of CH4, its coverage remains relatively low. Interestingly, the coverage of *HCOO on Pd/Al2O3-x increases rapidly from 6 min onwards and became the dominant surface species.On switching from CH4 to O2 atmosphere, *CO are immediately consumed on both Pd/Al2O3 and Pd/Al2O3-x, while the coverage of *HCOO increased abruptly. This may indicate that the oxidation of the CO intermediate is kinetically more favorable. More importantly, from the degradation rates of *HCOO we can infer that Pd clusters or single atoms are significantly more active than Pd NPs in catalyzing the oxidation of formate intermediates. As the (OCO)s vibrational frequency of the adsorbed *HCOO is correlated with the oxidation states of the Pd species, it can be deduced that the Pd-O coordination number of the surface Pd over the clusters or single atoms is lower than that over PdOx NPs in the atmosphere of O2, because Pd/Al2O3-x exhibits a higher *HCOO vibrational frequency (1358 cm−1) than Pd/Al2O3 (1349 cm−1). This may account for the higher activity of Pd/Al2O3-x towards *HCOO oxidation.In the previous section, we proposed the hypothesis that Pd single-atoms may also participate in CH4 oxidation without yielding Pd0 species. That is, the Pd single atoms may accomplish a catalytic cycle through a combination of oxidative addition, reductive elimination, and transfer insertion, which is similar to the extensively reported catalytic mechanism of a homogeneous catalyst41,42. In light of the difficulty to corroborate this hypothesis experimentally, ab initio DFT calculations may afford us an insight into CH4 oxidation pathways over the Pd single atoms.In spite of the considerable theoretical effort being expended on the mechanism of the CH4 oxidation over the specific facets of Pd or PdO, comparable work has rarely been done on the single-atom Pd species32,43,44,45,46,47,48,49. To implement this knowledge and examine our hypothesis, we have built a Pd1/Al2O3 slab in which an isolated Pd atom was anchored on the (100) facet of γ-Al2O3 via O bridges19 (Supplementary Information, Section 8.1 for details). Referring to the individual CH4 oxidation mechanisms postulated by previous studies, we proposed all CH4 oxidation pathways over a single-atom Pd, and subsequently screened and optimized these pathways with DFT calculations. The proposed mechanism network is illustrated in Fig. 5, where a-j represent sketches of the structures of the catalyst (a) and different intermediates (b-j). The numbers following the letters are the calculated Gibbs free energies of the corresponding intermediates relative to the initial catalyst, at 400 °C. The arrows symbolize the conversion between individual intermediates, and TS1-14 refer to the transition states present in these conversion processes, whose relative Gibbs free energies are also indicated below. The configuration of the individual intermediates and transition state species is shown in Supplementary information, Sections 8.2 and 8.4. For comparison purposes, an energy profile of the entire reaction network was made (Fig. 6).Fig. 5: A model of CH4 oxidation pathways on single-atom Pd.Reaction pathway network for CH4 oxidation on single-atom Pd obtained from density functional theory (DFT) calculations (top). The values are the corresponding Gibbs free energies of the intermediates and transition states, in kJ/mol. The blue, yellow, pink, and green arrows respectively represent the formate branch of *OCH3 pathway, the *CO branch of the *OCH3 pathway, the *OCH2 branch of the *CH3 pathway, and the interconversions between these pathways. The green shaded region represents the elementary steps for CH4 activation, leading to either *CH3 or *OCH3 pathways; while the yellow shaded region is the elemental steps where branching occurs between *CO and formate paths on the *OCH3 pathway.Fig. 6: Gibbs free-energy profile of CH4 oxidation on single-atom Pd.The letters and numbers annotated on the energy profile represent the intermediates and transition states along the various reaction pathways, corresponding to those depicted in Fig. 5. The green shaded region represents the elementary steps for CH4 activation, leading to either *CH3 or *OCH3 pathways; while the yellow shaded region is the elemental steps where branching occurs between *CO and formate paths on the *OCH3 pathway.It is worth noting that the well-known carbide pathway is highly unlikely on the single-atom Pd, as the adsorption of carbene, carbyl and *C species on the single-atom Pd is thermodynamically unfavorable (Supplementary information, Section 8.2 for details). The activation of CH4, highlighted in green shading in Fig. 5, is identified as a crucial rate-determining step in the overall CH4 oxidation reaction. Compared to the *CH3 pathway (marked by pink arrows), the *OCHx pathway (marked by blue arrows) is more favorable in terms of methane activation, as the *OCHx generation has a significantly lower activation energy than the *CH3 production. The subsequent oxidation of the carbonous intermediates all exhibits small energy barriers regardless of the pathways. The interconversions (marked by the green arrows) amongst the two pathways are all in favor of OCHx pathways from the microkinetic perspective. These DFT results are corroborating the dominance of the OCHx pathway in the catalytic mechanism.According to the proposed model, the *CO and *HCOO intermediates observed in the operando FT-IR spectra may have originated from the oxidation of *OCH2 by the coexistent active O (highlighted with yellow shadowing). Compared to *CO, the generation of *HCOO is thermodynamically more favorable while kinetically unfavorable, but the energy barrier difference (23.8 kJ/mol) between these two intermediates is much smaller than the thermodynamic Gibbs free energy difference (87.6 kJ/mol). This may explain the preferential appearance of *HCOO adsorbates. Notably, possibly because a reductive elimination leads to large changes in the valence of the single atom Pd, a high energy barrier for the desorption of water from the single-atom Pd was observed, which was higher than the activation energy of the CH4 molecule. This is also consistent with the operando FT-IR spectroscopy results, implying that water can poison the single-atom Pd centers.In comparison with the reaction pathways over a crystalline PdO surface unveiled by previous research43, single-atom Pd manifests 20.3 kJ/mol higher CH4 activation energy but lower energy barriers for further oxidation of carbon-containing intermediates. This explains the relatively cleaner surface of Pd/Al2O3-x under the reaction conditions as shown in the operando FT-IR spectra. The other rate-determining step to the single-atom Pd, the desorption of water, is less demanding for the crystalline PdO thermodynamically43. From a practical point of view, the low tolerance of Pd single atoms to water vapor, just as the poor resistance of PdOx NPs to carbonates, may be detrimental to their activity, regardless of their high dispersion or intrinsic activity. These are conceivably important reasons for the growth of apparent activation energy during the deactivation of commercial Pd/Al2O3 catalysts and highlight the potential of Pd clusters as more stable active centers for CH4 activation.To sum up, a highly dispersed Pd catalyst was synthesized using defective Al2O3-x as a support and exhibited higher catalytic activity and stability in CH4 oxidation compared to a counterpart Pd/Al2O3 catalyst. The evolution of the Pd active centers was significantly influenced by the type of support under the reaction conditions: the Pd clusters on γ-Al2O3 partially sintered and agglomerated into PdOx nanoparticles (NPs), meanwhile being partially anchored as single atoms. In contrast, only a small fraction of the initial Pd single atoms in Pd/Al2O3-x were fused into clusters due to the constraint effect of Alv species. The evolution of active centers under reaction conditions not only changes the quantity of Pd atoms available for participating in catalytic reaction, but also alters the intrinsic activities of the catalysts and their resistance to different poisoning species by changing the relative proportions of Pd NPs, clusters, and single atoms. The combined operando FT-IR spectroscopy, operando QEXAFS and DFT calculations revealed the mechanism of CH4 oxidation reactions over Pd NPs, clusters and single atoms. The higher dispersion and lower adsorbate coverage of PdOx clusters made them the most active in catalytic CH4 oxidation. This paper showcases a method to ensure the abundance of Pd clusters by imposing restrictions on the evolution of Pd with defective alumina. We foresee that similar methods may be applied to other supports and metals, paving the way to more thermally stable and highly dispersed metal catalysts.
Anchoring PdOx clusters on defective alumina for improved catalytic methane oxidation
