Open Access Article
Lydia J. Bachmann
ab,
Jagrati Dwivedi
a,
Dmitry Lapkin†
c,
Bihan Wangc,
Jan-Christian Schober
ab,
Gerard N. Hinsley
c,
Sarah Bernart‡
d,
Kuan Hoon Ngoi
c,
Rustam Rysov§
c,
Arti Dangwal Pandey
a,
Thomas F. Keller
ab,
Ivan A. Vartanyants
c and
Andreas Stierle
*ab
aCentre for X-ray and Nano Science CXNS, Deutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607 Hamburg, Germany. E-mail: andreas.stierle@desy.de
bDepartment of Physics, University of Hamburg, Notkestraße 9-11 Hamburg, 22607 Hamburg, Germany
cDeutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607 Hamburg, Germany
dInstitute of Catalysis Research and Technology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
First published on 20th March 2026
The high global warming potential of methane requires the development of advanced catalysts for conversion with increased activity and a longer lifetime. Here we utilize in situ Bragg coherent diffraction imaging (BCDI) to follow the evolution of two PtPd model catalyst nanoparticles during methane oxidation with different initial arrangements, one core–shell and one partially alloyed nanoparticle. While heating the nanoparticles under reaction conditions we observe the mixing and segregation of Pt and Pd, indicating a partial alloying and de-alloying, and a preferred segregation of Pt towards {111} oriented surface facets with potential use as bifunctional catalyst. Our in situ BCDI data also indicate a complete core–shell inversion at elevated temperatures with Pd moving into the nanoparticle. This suggests that a maximum operation temperature should not be exceeded, since Pd oxide is the more active phase. In spite of such limitations, binary metal core–shell nanoparticles are promising candidates for novel catalysts, as the elemental arrangement at their surface may be better controlled – a pivotal step for improving future catalysts.
Common catalyst nanoparticles (NPs) for methane oxidation consist of Pd, due to the highest activity for complete methane conversion of its oxide, or Pt as it is an active methane oxidation catalyst especially at lower pressures.2–4 Combining Pd with Pt allows to build NPs which are more stable.5–7 The surface chemical composition is crucial for the catalyst, therefore segregation and mixing processes need to be taken into account that may alter the catalyst during operation. It is known that PtPd alloy NPs preferentially form Pd oxide at the surface under oxidizing conditions.8 This Pd oxide formation can lead to Pd segregation to the surface, resulting in a metal Pt enriched core and a Pd oxide shell.4,9,10 Such segregation can significantly impact catalyst properties, as demonstrated by the increased methane oxidation activity associated with Pd oxide formation.11 The amount of the migrated Pd and the formed Pd oxide depends on the Pd to Pt ratio of the initial alloy, for very high Pt concentration less Pd oxide is formed.4,12
Another very important influence on the segregation and mixing process is the temperature. The maximum Pd concentration at the surface of PtPd under pure O2 is reached at about 500 °C (thermodynamic equilibrium), at higher temperatures the Pd concentration at the surface is decreasing.12 To determine the composition inside the NPs, strain measurements can give insights, since both elements have different lattice constants with a mismatch of 0.8% at room temperature (RT).13–16 A powerful method to investigate the strain on the nanometer scale is Bragg coherent diffraction imaging (BCDI), which enables to measure the shape and the strain inside single NPs in situ with a spatial resolution of about 5 nm.17–21
Apart from the increased stability, the combination of two materials can influence the activity. It has been shown, that core–shell systems can have different strain and electronic properties compared to pure or alloyed materials, which might boost the catalytic activity.22–25 In the present work we compare the mixing and segregation behavior of an initial core–shell NP with an initially partially alloyed NP under methane oxidation conditions at various temperatures between RT and 680 °C. We influenced the initial state by increasing the temperature of the second sample during the overgrowth with Pd, resulting in Pd diffusion into the NP. While the initial core–shell NP showed a facet orientation dependent segregation of Pt towards the {111} oriented facets between 470 °C and 560 °C during methane oxidation, the initially partially alloyed NP did not show such segregation behavior. At the final temperature (620 °C) at which the core–shell sample was analyzed, it was observed that Pd and Pt started to mix. Heating the initially partially alloyed NP to even higher temperatures (680 °C) revealed that the Pt is segregating to the surface, leading to a Pt shell and Pd rich core. These compositional changes at the surface are expected to have a high influence on the catalytic activity of the NPs.
The treatment after the EBL procedure was performed at the DESY NanoLab.28 Sample-150 was annealed in air at 950 °C for 2 h. Then Sample-150 was brought into UHV and annealed for 1 h at 410 °C, cleaned in 1 × 10−6 mbar O2 for 30 min at 410 °C, annealed again for 30 min at 410 °C and finally overgrown with Pd at 150 °C for 13 min. This procedure resulted in a nominal Pd thickness around (1.7 ± 0.7) nm as the calculations in SI section S2 show and a diffusion length of around 1.5 × 10−12 nm shown in SI section S3. Sample-670 on the other hand, was annealed in UHV by very slowly ramping to 1100 °C and keeping it at this temperature for 15 min. The sample was taken out of UHV conditions, to test the quality of such EBL samples for BCDI experiments at beamline P10 under N2 flow at RT. Afterwards, Sample-670 was annealed in UHV for 90 min at 670 °C and then overgrown with Pd at 670 °C for 1 h. This sample has a higher nominal thickness of (4.9 ± 0.2) nm and a much larger nominal diffusion length of around 0.4 nm as it is a function of the temperature and time, as described in SI sections S2 and S3.
To retrieve the crystalline real space shape from each dataset, the collected images were stacked and cut to individual region of interests, centered around the (111) Bragg peak as displayed in SI section S6. Then, an iterative phase retrieval algorithm was applied via the software PHASOR.30 Thereby, different algorithm sequences were used. On Sample-670 either four cycles of 300 iterations difference map (DM) followed by 300 iterations error reduction (ER) or three cycles of 300 iterations hybrid input–output (HIO) followed by 300 iterations ER or three cycles of 300 iterations DM followed by 300 iterations ER were used, all with a shrink-wrap threshold of 0.2. For Sample-150 the reconstruction was more difficult, since the high strain of the core–shell structure made it more challenging for the algorithm to converge, so additionally to the same sequences as on Sample-670 also the shrink-wrap threshold was adapted between 0.07 and 0.2 for each dataset. The sequences for each dataset are listed in SI section S6. Afterwards, the results of the phase retrieval algorithm were transformed into laboratory coordinates and the ten best reconstructions were selected and averaged. The voxel size in each reconstruction for Sample-670 is approximately 5 × 5 × 5 nm3, while for Sample-150 it is approximately 7 × 3 × 3 nm3. Due to the different ROI and angular ranges, the voxel size in each dataset differs slightly. The final values are listed in SI section S6.
The strain εzz is defined along the (111) direction, so the scattering vector g is parallel to the z-axis. With εzz(r) = ∂uz(r)/∂z the strain can be calculated from the displacement field u(r). This displacement field is connected to the reconstructed phase of the real space object by Φ(r) = −u(r)·g, so it follows u(r) = −Φ(r)/|g| = −Φ(r)/(2π/dhkl) = −dhklΦ(r)/(2π) with the distance dhkl between the lattice planes in (hkl) direction. So for the 111 Bragg peak
with aNP, determined from the Bragg peak position as described in SI section S7. A discussion of the strain error can be found in SI section S12.
For comparability the limits of the color scales in this manuscript were set to ±0.334% for all images displaying the strain. This limit was chosen by taking the average strain of a dataset with high strain values (dataset 22 on Sample-670) and adding the standard deviation to it.
Throughout the experiment a total pressure of 0.1 bar with a total flow of around 50 mL min−1 was applied. Three different gas mixtures were flown to the reactor: either 50 mL min−1 Ar (hereafter referred to ‘Ar atmosphere’) or 48 mL min−1 Ar + 2 mL min−1 H2 (‘H2 treatment’) or 28.5 mL min−1 Ar + 2.6 mL min−1 CH4 + 20 mL min−1 O2 (‘methane oxidation conditions’). Detailed gas and temperature conditions can be found in SI section S9 for both samples.
Comparing the SEM images on Sample-150 before and after the Pd overgrowth (Fig. 1a and b), reveals small Pd NPs on the substrate and the Pt NP side facets. In contrast, an SEM image of Sample-670 after Pd overgrowth (Fig. 1e) does not show similar Pd NPs on the substrate or the Pt NP, indicating that the higher overgrowth temperature either produced a smooth Pd film or caused Pd diffusion into Pt. Furthermore, the {100} facets on Sample-150 are larger than the {111} facets resulting in a more triangular shape, while Sample-670 has a more hexagonal shape and the facet orientation is not identifiable in the SEM image. Regarding particle dimensions, Sample-150 has a diameter of 275–283 nm and is therefore around 100 nm larger than Sample-670 with a diameter of 179–195 nm. The heights of the particles were determined by AFM, shown in (Fig. 1c and 1f). The line profiles extracted from these AFM images, shown in (Fig. 1d and 1e), identify that Sample-150 has a height of (91 ± 1) nm, and is approximately only 13 nm taller than Sample-670 which has a height of (78 ± 1) nm. The measured height and diameter together with the nominal thickness of the overgrown Pd, lead to an estimated nominal volume fraction of 2 vol% Pd and 98 vol% Pt on Sample-150 assuming a complete coating and a simplified NP shape, as shown in SI section S4. Similarly, Sample-670 has an estimated nominal volume fraction of 6 vol% Pd and 94 vol% Pt.
The initial state of the NPs was also characterized by BCDI. Thereby, the intensity distribution around the 111 Bragg peak was measured by scanning the incident angle in a range of ±0.4° to ±0.6°, as explained in the Experimental section. The height of the crystalline part of the NPs can be calculated from the fringes of this intensity distribution as discussed in SI section S5, resulting in heights of (97 ± 7) nm for Sample-150 and (77 ± 3) nm for Sample-670. These are in good agreement with the heights measured by AFM. Applying phase retrieval algorithms to the measured intensity distributions allows to reconstruct the shape and strain of the crystalline part of the NPs, as described in the Experimental section. The reconstructed shapes of both initial NPs are shown in (Fig. 2a and 3a), for Sample-150 and Sample-670, respectively, which agrees well with the shapes measured by SEM and AFM. A list with the specific applied algorithms for each dataset can be found in SI section S6.
The initial strain distribution at RT and under Ar atmosphere (50 mL min−1) is shown in Fig. 2 for Sample-150. The strain was calculated as described in the Experimental section with the lattice constant aNP = (3.920 ± 0.001) Å calculated from the observed Bragg peak position. The most distinctive feature is the high tensile strain (positive sign) inside the NP and the compressive strain (negative sign) towards the surface. This is clearly visible by the two peaks in the histogram (Fig. 2f) of the strain values of the cut shown in (Fig. 2e), which can be fitted by two Gaussian functions as described in SI section S8. Such bimodal strain distributions can be observed when investigating core–shell materials.13–15 In this case, we associate compressive strain in Pd rich and tensile strain in Pt rich regions, because aPd = 3.890 Å < aNP < aPt = 3.923 Å. We calculate the strain from the lattice constant aPt of Pt relative to the lattice constant aNP from the measured Bragg peak position, using the strain definition shown in SI section S7. For the initial state, this results in a tensile strain of (0.08 ± 0.02)%. This is the same order as the mean of (0.165 ± 0.006)% fitted in the histogram which indicates that the areas with tensile strain are regions of pure Pt. For a region with 100% Pd, the compositional strain would be (−0.77 ± 0.02)% compared to aNP. This strain value is significantly lower than the mean of (−0.032 ± 0.009)% fitted for the second peak in the strain histogram and even lower than the minimum strain value of around −0.25%. This discrepancy is related to the difference between the observed thickness of the compressive strain ring of (27 ± 9) nm, and the much smaller nominal thickness of the evaporated Pd of (1.7 ± 0.7) nm. Both is an indication that the observed ring of compressive strain is not only compositional strain from the Pd shell, but primarily mechanical strain in the Pt whose lattice is adapting to the Pd shell. Similar lattice changes of the core material to the shell material has been shown in previous studies.15,32,33 This effect of the overgrown shell material on the core lattice depends on several parameters. First, the lattice mismatch between both metal materials cannot be larger than 5%. Otherwise, epitaxic growth of the shell material on the core can be excluded.32 The lattice mismatch between Pt and Pd of 0.8% is much smaller than this limit. Second, the shell material should be able to grow epitaxially on the core material if the temperatures during and after the overgrowth are low enough to suppress any inter-diffusion. For Sample-150 the overgrowth temperature was 150 °C, which does not imply any diffusion. Third, the thickness of the overgrown layer must be below a critical thickness, otherwise the shell lattice will incorporate dislocations rather than expanding or compressing the core lattice. The critical thickness hc for coherent growth of Pd on Pt is calculated to be around 56 nm at RT as described in SI section S13. The nominal thickness of Pd on Sample-150 is (1.7 ± 0.7) nm, far below this critical thickness. Simulations of Pd on Au have shown that, if all three conditions are fulfilled, an evaporated layer of 4 nm Pd can have an influence on the lattice constant of an Au core up to 20 nm inside the NP, which agrees well with our observation.32
The initial strain distribution of Sample-670 is shown in Fig. 3 and the strain was calculated with the lattice constant aNP = (3.920 ± 0.001) Å from the observed Bragg peak position, which is identical to the lattice constant of the initial state of Sample-150. In comparison to Sample-150 the strain in Sample-670 is more mixed, with positive and negative strain values in the center and near the surface. The origin of this non core–shell structure in Sample-670 lies in the higher overgrowth temperature, which resulted in Pd diffusion into the Pt NP. Although the calculated diffusion length on Sample-670 is still smaller than the nominal Pd thickness of (4.9 ± 0.2) nm, it was sufficient to reduce the strain induced by the overgrown Pd. Simultaneously, also Pt surface segregation can occur.34
Fig. 4 displays the reconstructed shape of a selection of datasets on Sample-150, showing some changes throughout the in situ experiment. Foremost, when comparing (Fig. 4a) to (Fig. 2), the strain at the top facet has an opposite strain. This change of sign appears after elevating the temperature to 470 °C and exposing the sample to methane oxidation condition for the first time. It was not possible to measure under this condition due to experimental issues, instead the dataset displayed in (Fig. 4a) was measured at RT and under Ar atmosphere after this first exposure. The tensile strain at the top facet indicates Pt segregation through the Pd shell to the surface of the top facet. A H2 treatment at 320 °C was leading to Pt segregation towards other facets as (Fig. 4b) shows. This is in line with previous findings, showing that H2 treatment at elevated temperatures promotes the movement of Pt atoms to the surface in PtPd NPs.34 As the black arrows in (Fig. 4b) indicate, the Pt only segregated to the {111} facets. Similar facet orientation dependent segregation behavior towards {111} oriented facets was observed for a PtRh NP under CO conditions.16 In our case, Pt stayed at the {111} facets, after switching back to methane oxidation conditions and ramping the temperature to 520 °C. Only after increasing the temperature further to 560 °C, Pt disappeared from the surfaces, as (Fig. 4c) shows. This could be related to preferred Pd segregation and surface oxide formation at this condition.12 The strain inside the NP stays the same, unaffected by the Pt segregation to the {111} facets and back, as the horizontal (−20 nm) and the vertical (‘cf’) cuts in Fig. 4 show.
![]() | ||
| Fig. 4 Cuts through selected datasets on Sample-150 (a) At RT under Ar atmosphere after first exposure to methane oxidation conditions at 470 °C. (b) At 320 °C under Ar atmosphere, after H2 treatment during the ramp up. The black arrows are pointing towards the {111} oriented facets, where we observed the tensile strain. (c) At 560 °C under methane oxidation conditions. (d) At 620 °C under methane oxidation conditions, the strain is decreasing and the ring structure seems to vanish. (e) Histogram of the strain values from the horizontal cut at 20 nm of the dataset shown in (c). The histogram was fitted using two Gaussians, with the peak positions indicated by the dashed lines at −(0.163 ± 0.018)% and (0.044 ± 0.005)%. (f) Histogram of the strain values from the horizontal cut at 20 nm of the dataset shown in (d). The histogram was fitted using one Gaussian, with the peak position indicated by the dashed line at (0.015 ± 0.006)%. All datasets are cut as indicated in Fig. 2 and the NP facets are indicated in the same manner. The numbers marked by # denotes the number of the datasets. | ||
When increasing the temperature further to 620 °C, the strain at the surface and inside the NP decreases, as seen in the vertical and horizontal cuts in (Fig. 4d). This decrease of the strain could be the first step towards the disappearance of the core–shell structure, as the ring structure in the horizontal cut nearly vanished. The fits of the histograms support this assumption, since one Gaussian in (Fig. 4f), instead of two in (Fig. 4e), is sufficient to fit the histogram and results in decreasing mean values from −(0.163 ± 0.018)% and (0.044 ± 0.005)% at 560 °C to (0.015 ± 0.006)% at 620 °C. Additionally, the minimal values in the cuts are getting higher, from −0.33% to −0.2%. Additional values and histograms are shown in SI section S8. Overall, Pt and Pd seem to start alloying at this temperature.
As one can see in Fig. 5 Sample-670 does not show such kind of facet-dependent Pt segregation even at temperatures up to 680 °C (see SI section S9), which might originate from the different initial state.
![]() | ||
| Fig. 5 Cuts through selected datasets on Sample-670 (a) At RT under Ar atmosphere. (b) At 320 °C under methane oxidation conditions. (c) At 680 °C under methane oxidation conditions. (d) At 320 °C under methane oxidation conditions, after switching to Ar atmosphere and then to a H2 treatment at this temperature, as shown in Supplementary section S9. (e) At 680 °C under methane oxidation conditions. (f) Histogram of the strain values from the horizontal cut at −20 nm shown in (b). The histogram was fitted using one Gaussian, with the peak position indicated by the dashed line at −(0.014 ± 0.057)%. (g) Histogram of the strain values from the horizontal cut at −20 nm shown in (c). The histogram was fitted using two Gaussians, with the peak positions indicated by the dashed lines at −(0.113 ± 0.004)% and (0.066 ± 0.005)%. All datasets are cut as indicated in Fig. 3 and the NP sides are indicated in the same manner. The gaps in the reconstructed shape do not result from real gaps in the NP, but rather from failure in the reconstruction at those points due to high strain values. | ||
After applying several cycles of methane oxidation condition with maximal 550 °C and including H2 treatments (see SI section S9), Sample-670 does not change (Fig. 5a and b) in comparison to the initial state in Fig. 3. Upon increasing the temperature further to 680 °C, the strain distribution is however changing drastically. As shown in (Fig. 5c) a clear ring structure is surprisingly present now, with high tensile strain at the surface (max. 0.35% in horizontal cut) and high compressive strain at the inside (min. −0.37% in horizontal cut, see SI section S8 for the values). This range of the compressive strain values agrees with the calculated compositional strain of (−0.46 ± 0.2)% for 100% Pd compared to the calculated lattice constant from the Bragg peak center of this dataset at 680 °C (see SI section S7 for the calculation). The presence of a Pd rich core Pt shell structure is further confirmed by comparing observed features with literature on de-alloying NPs with different materials.35,36 Especially, the splitting of the strain histogram into two peaks (from Fig. 5f to g) and the increase of the asymmetry in the line profile through the fringes of the intensity distribution around the Bragg peak (shown in SI section S10) are typical for de-alloying.37,38 Furthermore, the radius r = (40 ± 5) nm of the region with compressive strain agrees with the radius r = 42 nm of a half-sphere filled by the volume the of evaporated Pd, see SI section S11. Altogether one can conclude a Pd rich core, Pt shell structure was formed, thus an inverted core–shell structure compared to the initial state of Sample-150. The Pd rich core, Pt shell structure was not changed by decreasing the temperature and applying an H2 treatment (see Fig. 5d) and when increasing the temperature to 680 °C for the second time (Fig. 5e).
Apart from the temperature, the biggest influence is the different initial state. Sample-150 has a Pd shell at the beginning of the in situ experiment, due to the lower overgrowth temperature. Thus more Pd has to be transported further into the NP and more Pt to the surface to form a Pt shell Pd rich core structure. Furthermore, we speculate that the higher availability of Pd at the surface has led to more oxidized Pd in Sample-150, which needs to be reduced before Pd can segregate into the bulk. Taking this into account, we assume that the observed mixing of the Pt and Pd on Sample-150 is a first step towards a Pt shell, Pd core like in Sample-670, inline with a core–shell inversion.
Next, we discuss how the initial state and its transformation through segregation and mixing during the in situ experiment may influence the catalytic activity. The different initial states of both samples are expected to have a significant impact on the performance. It has been shown in several experiments and simulations, that the core–shell structure of NPs can boost their catalytic activity compared to alloyed NPs of the same materials, since the core–shell structure has a high impact on the strain.22–25 Such increased strain of the core–shell system, compared to the partially alloyed system was also observed in our experiment. Additionally, the initial core–shell NP might be more stable compared to PtPd alloy NPs, because no mass transport is needed to form Pd oxide. The partially alloyed state on the other hand is not preferential, since the high Pt concentration prevents the formation of Pd oxide, so the catalytic performance may be close to the performance of pure Pt which is lower than the activity of Pd oxide.4–6
Furthermore, the segregation and mixing is influencing the catalytic performance of both samples as well. For the initial Pt core, Pd shell Sample-150 it was observed that Pt segregates to {111} oriented facets under methane oxidation conditions at temperatures up to around 520 °C. The existence of Pt segregated to the Pd rich surface might be beneficial for the activity for methane oxidation for partial oxygen pressures up to 25 mbar, above this pressure Pd oxide facets are more preferred sites for methane oxidation.4 On the other hand, such facet dependent surface composition could be utilized to prepare bifunctional catalyst NPs. In case of the initially partially alloyed Sample-670, we observed Pt segregation to the surface at 680 °C, resulting in a Pd rich core, Pt shell structure. Other than the numerously observed Pd segregation to the surface,4,8–10 is Pt segregation is not driven by oxide formation but it thermodynamically favored. In contrary, the Pt segregation was promoted because Pd oxide gets reduced at such high temperatures.4,44 Since Pd is more active than Pt,4–6,45 this compositional change might have a negative effect on the activity of the NP. Thus, such PtPd alloyed catalyst NPs should only be operated for temperatures below 680 °C for methane oxidation.
Supplementary information (SI): additional calculations, images and further information. See DOI: https://doi.org/10.1039/d5nr05321h.
:
Pt ratio of Pd/Pt bimetallic catalyst on CH4 oxidation, Appl. Catal., B, 2022, 316, 121623 CrossRef CAS.Footnotes |
| † Current address: Institute of Applied Physics, University of Tübingen, 72076 Tübingen, Germany. |
| ‡ Current address: Inorganic Materials and Catalysis, Eindhoven University of Technology, 5600 MB Eindhoven, Netherlands. |
| § Current address: European X-ray Free-Electron Laser Facility, Holzkoppel 4, 22869 Schenefeld, Germany. |
| This journal is © The Royal Society of Chemistry 2026 |