Open Access Article
Masahiko Shimizu
*abc,
Yuta Inami
a,
Ryuichi Shimogawa
ad,
Takeshi Matsuoa,
Yu Fujikata
a,
Hajime Matsumotoab,
Kazutaka Mitsuishi
*e and
Ayako Hashimoto
*bc
aScience & Innovation Center, Mitsubishi Chemical Corporation, 1000 Kamoshida-cho, Aoba-ku, Yokohama, Kanagawa, Japan. E-mail: masahiko.shimizu.ma@mcgc.com
bResearch Center for Energy and Environmental Materials, National Institute for Materials Science, 1-2-1 Sengen, Tsukuba, Ibaraki, Japan. E-mail: hashimoto.ayako@nims.go.jp
cGraduate School of Science and Technology, University of Tsukuba, 1-2-1 Sengen, Tsukuba, Ibaraki, Japan
dDepartment of Materials Science and Chemical Engineering, Stony Brook University, Stony Brook, New York 11794, USA
eCenter for Basic Research on Materials, National Institute for Materials Science, 1-2-1 Sengen, Tsukuba, Ibaraki, Japan. E-mail: mitsuishi.kazutaka@nims.go.jp
First published on 27th April 2026
Rational design of high-performance bimetallic nanocatalysts requires an understanding of the unique atomic structures governing their performance. This study focused on the Re–Ge/TiO2 catalyst, which exhibits high performance for carboxylic acid hydrogenation. We analyzed structural and electronic state changes occurring during its multistep preparation involving calcination in air, hydrogen reduction, and oxidative stabilization to elucidate the key structural factors responsible for its high performance. To this end, a complementary approach combining in situ X-ray absorption fine structure analysis and ex situ scanning transmission electron microscopy (STEM) was employed. In particular, the STEM measurements utilized an air-free transfer holder to track the same individual nanoparticles throughout all preparation steps, thereby addressing two issues: the challenge of distinguishing true structural changes from particle-to-particle variations among different particles and electron beam damage associated with long exposure times. It was revealed that the highly active catalytic state after hydrogen reduction is associated with approximately 1 nm crystalline Re–Ge alloy nanoparticles. The crystal structure of the nanoparticles was a unique, low-energy, and face-centered cubic fragment, verified by density functional theory calculations. Furthermore, quantitative STEM image analysis demonstrated that the nanoparticles formed a Re–Ge random alloy. This atomic-level mixing, as detected by the spatially averaged X-ray absorption, is considered to stabilize the metallic Re(0) state, rendering Re electron-rich. These findings not only provide a new strategy for designing high-performance nanocatalysts but also establish a correlative methodology as a way to identify performance-determining factors in complex nanomaterials.
We recently developed a novel nanocatalyst by combining Re with germanium (Ge).16,17 It exhibits higher selectivity in direct hydrogenation of carboxylic acids compared with monometallic Re catalysts and has high activity and selectivity comparable to commercial Ru–Pt–Sn catalysts.18 Furthermore, its precious-metal-free composition and low-cost preparation make it a promising candidate for industrial applications. However, the atomic-level structural and electronic origins of this performance enhancement remain unclear. Understanding these factors is essential for further enhancements and for rational design of next-generation nanomaterials. This requires accurately tracking atomic-level structural changes throughout multistep thermal treatments during catalyst preparation, which is difficult to achieve with conventional methods. For example, in situ X-ray techniques provide spatially averaged information under reaction conditions18,19 but do not resolve the local structures of individual active sites.20,21 Atomic-resolution scanning transmission electron microscopy (STEM) enables local structural analysis.22–26 However, since structural changes associated with each preparation step are typically assessed by observing different samples, the effects of sample inhomogeneity—such as variations in particle size and shape—cannot be excluded. Therefore, it is challenging to distinguish whether observed structural differences originate from changes during the catalyst preparation or reflect simply particle-to-particle variations. Although in situ STEM enables dynamic imaging, prolonged electron beam exposure can cause significant beam damage to the sample.27,28
We addressed these challenges through a complementary multi-scale approach. Using in situ X-ray absorption fine structure (XAFS) analysis, we tracked the spatially averaged chemical state changes across the entire sample. We then used ex situ STEM with an air-free transfer holder to directly image the corresponding structural changes of individual nanoparticles at the atomic scale. This approach enabled us to elucidate the changes that occur in both the structural and electronic states occurring during each preparation step of Re–Ge catalysts and to identify the key structural factors responsible for the high performance. Unique crystalline nanoparticles of approximately 1 nm in diameter were observed in their highly active state after hydrogen reduction. Furthermore, XAFS analysis indicated that Re–Ge alloy formation stabilized the metallic Re(0) state, rendering Re electron-rich. These results not only confirm a new strategy for designing high-performance nanocatalysts but also demonstrate an effective analytical methodology for the structural analysis of various complex nanomaterials.
For the calcined and oxidized Re–Ge/TiO2 catalysts, powder XRD patterns were re-measured using a diffractometer (Empyrean, Malvern Panalytical) with Cu Kα radiation (45 kV, 40 mA) in reflection mode. Data were collected over a 2θ range of 10–80° with a step size of 0.008°. Each sample was measured by four repeated scans.
535 eV) and Ge K edge (11
103 eV) spectra were obtained in transmission mode. In situ measurements were performed using an electric furnace cell equipped with an automated gas mixing system. The Re–Ge/TiO2 or Re/TiO2 sample was loaded into the cell, and spectra were first recorded in air. The sample was then reduced under flowing hydrogen as the temperature was increased from room temperature to 500 °C at 10 °C min−1. XAFS spectra were recorded continuously during the temperature ramp. After holding it at 500 °C, the sample was then cooled to 50 °C under hydrogen. The gas atmosphere was then switched sequentially from H2 to N2, 5% O2/N2, and 20% O2/N2 at 50 °C to investigate the oxidation behavior. For comparison, Re/TiO2 (without Ge) was analyzed under identical conditions.
XAFS data processing was performed with the Larch software package.29 Energy calibration was conducted using Re and Ge metal references, where the edge energy (E0) was defined as the first zero-crossing point of the second derivative of the absorption spectrum. This was 10
535 eV for the Re L3 edge and 11
103 eV for the Ge K edge. Extended XAFS (EXAFS) fitting was performed in R-space with FEFF8.5L software bundled with Larch. Details of the fitting procedure are provided in the SI.
Fig. 1 shows the XAFS spectra of Re–Ge/TiO2 at different treatment stages, while Fig. S1 shows the temperature-dependent changes in the X-ray absorption near edge structure (XANES) spectra during the reduction step. Fig. S2 presents the XANES spectra in Fig. 1 overlaid with reference compounds: Re metal, ReO2, ReO3, and NH4ReO4 for the Re L3 edge and Ge(0) and hexagonal GeO2 for the Ge K edge. Before reduction, the Re L3 edge spectrum exhibited a peak position comparable with that of Re(IV), Re(VI), and Re(VII) references, confirming the presence of oxidized Re compounds (Fig. S2a). The corresponding Fourier transform spectrum (Fig. 1c) matches well with that of NH4ReO4 in R space. After hydrogen reduction at 500 °C, the white line intensity significantly decreases, and the edge position shifts to lower energy. However, the position remains higher than that of Re(0), indicating partial reduction toward the metallic state (Fig. S2a). The Fourier transform spectrum shows that the intensity of the peak corresponding to the first coordination shell decreases as the reduction proceeds (Fig. 1c).
Similarly, for the Ge K edge, the XANES spectrum before reduction matches well with that of GeO2, confirming that Ge exists in the Ge(IV) oxidized state (Fig. S2b). The Fourier transform spectrum (Fig. 1f) also shows good agreement with that of GeO2 in R space. After hydrogen reduction, the edge position shifts to lower energy, and the intensity of the first-shell peak decreases. However, the white line intensity remains higher than that of Ge(0), indicating partial reduction, as observed for Re.
The temperature-dependent XANES spectra (Fig. S1) provide insight into the reduction mechanism. For the Re L3 edge (Fig. S1a), the absence of an isosbestic point during reduction indicates that the transformation involves more than two chemical compounds. This suggests that the reduction of Re proceeds through multiple intermediate oxidation states, from Re(VII) and/or Re(VI) through Re(IV) to Re(0), rather than through a simple two-state conversion. Such reductions involving multiple intermediate species have also been reported in the Re/TiO2 system.10 In contrast, the Ge K edge spectra (Fig. S1b) exhibit behavior consistent with a two-component system, suggesting that the reduction of Ge proceeds directly from Ge(IV) to Ge(0) without detectable intermediate oxidation states.
To quantitatively analyze the structural evolution during each preparation step, EXAFS curve-fitting analyses were performed, as summarized in Tables 1 and 2. Fig. 2, Tables S2 and S3 show the temperature dependences of the structural and fitting parameters during hydrogen reduction. The individual fitting results in R-space are presented in Fig. S3–S5.
| Sample | Path | N | R (Å) | σ2 (×10−3 Å2) | ΔE0 (eV) | R-Factor | χ2v |
|---|---|---|---|---|---|---|---|
| Re–Ge/TiO2 before reduction | Re–O | 5.3 ± 2.7 | 1.775 ± 0.032 | 0.0 ± 1.9 | 7.5 ± 8.3 | 2.9% | 557.0 |
| Re–Ge/TiO2 H2 500 °C | Re–O | 1.2 ± 1.4 | 1.899 ± 0.074 | 0.1 ± 11.9 | −11.0 ± 12.3 | 4.9% | 135.1 |
| Re–Ge | 2.8 ± 3.3 | 2.448 ± 0.066 | 10.7 ± 8.8 | ||||
| Re–Ge/TiO2 O2 50 °C | Re–O | 1.3 ± 1.4 | 2.067 ± 0.055 | 0.1 ± 5.0 | 14.6 ± 10.0 | 5.4% | 81.7 |
| Re–Ge | 1.8 ± 2.3 | 2.554 ± 0.051 | 6.0 ± 10.3 |
| Sample | Path | N | R (Å) | σ2 (×10−3 Å2) | ΔE0 (eV) | R-Factor | χ2v |
|---|---|---|---|---|---|---|---|
| Re–Ge/TiO2 before reduction | Ge–O | 4.7 ± 2.2 | 1.729 ± 0.034 | 0.0 ± 5.1 | −1.0 ± 7.2 | 5.0% | 1092.0 |
| Re–Ge/TiO2 H2 500 °C | Ge–O | 1.5 ± 0.4 | 1.735 ± 0.018 | 0.0 ± 3.6 | −10.0 ± 4.0 | 1.8% | 29.8 |
| Ge–Re | 3.0 ± 2.1 | 2.547 ± 0.029 | 12.7 ± 7.0 | ||||
| Re–Ge/TiO2 O2 50 °C | Ge–O | 3.2 ± 1.4 | 1.738 ± 0.032 | 0.6 ± 5.0 | 0.7 ± 6.6 | 4.6% | 881.3 |
A significant difference was observed in the EXAFS spectra in R-space for Re–Ge/TiO2 and Re/TiO2 at 500 °C under hydrogen (Fig. 3b), indicating the presence of Re–Ge bonds in the bimetallic catalyst. A two-shell fit with Re–O and Re–Ge paths reproduced the EXAFS signals, yielding a Re–Ge bond distance of 2.448 ± 0.066 Å, which is significantly shorter than the Re–Re bond distance reported for Re/TiO2
10 (2.76 ± 0.01 Å). In contrast, a two-shell fit with Re–O and Re–Re paths did not adequately reproduce the experimental data owing to the longer bond distance of the Re–Re shell.
For the Re L3 edge (Fig. 2a and b), the Re–O coordination number, N(Re–O), and the Re–O bond distance, R(Re–O) (∼1.7 Å) remain consistent with the tetrahedral ReO4− structure of NH4ReO4 below 280 °C. Above this temperature, N(Re–O) begins to decrease, indicating the onset of Re reduction. At intermediate temperatures (approximately 335 °C), EXAFS fitting required both short and long Re–O paths (Fig. S4l, Table S2), indicating the coexistence of multiple Re–O coordination environments consistent with the multi-step reduction suggested by the XANES data. Above 335 °C, the model with a long Re–O bond provided physically reasonable parameters.
Above 335 °C, N(Re–Ge) begins to increase as N(Re–O) continues to decrease, indicating the formation of Re–Ge bonds concurrently with the loss of Re–O coordination. However, owing to the strong correlation between the coordination number and the Debye–Waller factor (σ2) in the EXAFS fitting, the error bars associated with N(Re–Ge) and σ2(Re–Ge) are too large to allow a quantitative comparison; therefore, only a trend of increasing N(Re–Ge) can be inferred. At 500 °C, the fitting results show N(Re–Ge) = 2.8 ± 3.3 with R(Re–Ge) = 2.448 ± 0.066 Å (Table 1).
For the Ge K edge (Fig. 2c and d), N(Ge–O) and R(Ge–O) (∼1.73 Å) remain consistent with those of GeO2 below 300 °C (Fig. S5). Above this temperature, N(Ge–O) decreases similarly to Re, suggesting concurrent reduction and alloy formation. It is noteworthy that even after 500 °C reduction under hydrogen, a significant Ge–O contribution remains with N(Ge–O) = 1.5 ± 0.4 (Table 2), indicating that a substantial fraction of Ge exists as oxides on the TiO2 support.
For the spectrum at 500 °C, fitting with a two-shell model including Ge–O and Ge–Re paths yielded N(Ge–Re) = 3.0 ± 2.1 and R(Ge–Re) = 2.547 ± 0.029 Å (Table 2). This spectrum could also be fitted with a two-shell model including Ge–O and Ge–Ge paths, yielding N(Ge–Ge) = 1.2 ± 1.3 and R(Ge–Ge) = 2.431 ± 0.033 Å (Table S1). This Ge–Ge bond distance is comparable to that in Ge bulk (2.45 Å).37 Although the Ge–Re model provided a marginally better fit (R-factor = 1.8% vs. 3.3%), both structural models remain plausible given the similar fit quality, and the Ge K-edge data alone do not allow us to exclude Ge–Ge bonding. The fitted Ge–Re bond distance (2.547 Å) is slightly longer than the Re–Ge distance obtained from the Re L3 edge analysis (2.448 Å). Nevertheless, these values are comparable within experimental uncertainty. Considering the clear Re–Ge contribution observed from the Re L3 edge, we conclude that Ge–Re bonds are present in the reduced catalyst. However, the Ge–Ge contribution may also be significant, suggesting that both Ge–Re and Ge–Ge coordination environments coexist in the reduced catalyst.
The XANES spectra of Re–Ge/TiO2 and Re/TiO2 were compared to investigate the effect of Ge addition on the electronic state of Re. Fig. S6 shows the Re L3 edge XAFS analysis of the Re/TiO2 catalyst at different treatment stages and Fig. 3 presents a direct comparison between Re–Ge/TiO2 and Re/TiO2. After hydrogen reduction at 500 °C, the Re L3 edge XANES spectrum of Re–Ge/TiO2 exhibits a lower white line intensity and a slight shift to lower energy compared with that of Re/TiO2 (Fig. 3a). This decrease in white line intensity indicates that Re in the bimetallic catalyst has a higher electron density, suggesting that the formation of the Re–Ge alloy stabilizes the metallic Re(0) state more effectively than in the monometallic catalyst. This electron-rich Re is closely related to the enhanced Re–Ge catalytic performance.
The structural changes were investigated to understand the origin of the decreased catalytic performance after oxidative stabilization. Upon exposure to 20% O2 at 50 °C, distinctly different behaviors were observed for Re–Ge/TiO2 and Re/TiO2 catalysts (Fig. 3). For Re/TiO2, significant changes in both the XANES spectrum and the Fourier transform magnitude were observed after oxidation (Fig. S6), indicating substantial re-oxidation of Re(0). In contrast, Re–Ge/TiO2 exhibited relatively minor changes in the Re L3 edge spectra after oxidative stabilization (Fig. 1a–c), with the white line intensity remaining lower than that of Re/TiO2 (Fig. 3a). This suggests that the Re–Ge alloy structure provides enhanced oxidation resistance when compared with Re/TiO2.
The EXAFS fitting results (Tables 1 and 2) also support this interpretation. For the Re L3 edge, the Re–Ge bond remains detectable after oxidative stabilization with N(Re–Ge) = 1.8 ± 2.3 and R(Re–Ge) = 2.554 ± 0.051 Å (Table 1). Although no definitive conclusion can be drawn regarding the coordination number owing to the strong N–σ2 correlation in the fit, the persistence of the Re–Ge contribution suggests that the Re–Ge alloy structure formed at high temperature is at least partially preserved, indicating that the Re–Ge alloy layer has higher resistance to oxidation. In contrast, at the Ge K edge, the coordination number N(Ge–O) increases from 1.5 ± 0.4 (H2, 500 °C) to 3.2 ± 1.4 (O2, 50 °C) after oxidation (Table 2), while the second shell peak corresponding to the Ge–Re (or Ge–Ge) bonds significantly decreases in intensity (Fig. 1f, comparing the green and red curves). This indicates that the reduced Ge compounds have been substantially re-oxidized.
These observations suggest that the addition of Ge leads to the formation of a stable Re–Ge alloy layer that confers oxidation resistance, thereby stabilizing the Re(0) state. Meanwhile, most of the Ge exists as Ge–Ge-rich regions or oxide species on the support and are readily converted to GeO2 under oxidizing conditions. The re-oxidation of Ge compounds within the alloy nanoparticles upon oxidative stabilization is expected to disrupt the atomic-level mixing of Re and Ge that characterizes the active state. This would potentially lead to structural degradation of the nanoparticles.
Fig. 4a–c show the ex situ HAADF-STEM images of the same field of view at each preparation step. After calcination in air (Fig. 4a), both Re and Ge species were found to be highly dispersed on the TiO2 support, primarily as sub-nanometer clusters and some individual atoms. In contrast, nanoparticles were formed after hydrogen reduction (Fig. 4b). The particle size distribution of approximately 700 particles, according to a reported method,38 yielded an average diameter of 1.3 ± 0.3 nm (see the histogram in Fig. S8). The nanoparticles remain clearly visible after the oxidative stabilization step (Fig. 4c), with no observable changes in their size or position. These observations directly demonstrate that approximately 1 nm nanoparticles are formed from the dispersed species during hydrogen reduction. Furthermore, their size and position on the support remain unchanged at the atomic scale during oxidation stabilization. Fig. 4d and e show the HAADF-STEM images of the same particle after hydrogen reduction and after oxidation stabilization, respectively. The highly active particles after hydrogen reduction exhibit lattice fringes. In contrast, the fringes disappear for the low-activity particles after oxidative stabilization. This loss of crystallinity was also confirmed by the disappearance of the diffraction spots from the nanoparticles in the corresponding fast Fourier transform (FFT) images (Fig. 4f and g and Fig. S9). Under our imaging conditions, control electron beam irradiation experiments confirmed that these crystallization/amorphization changes were not beam-induced but resulted from the respective preparation steps (Fig. S10).
The Re–Ge/TiO2 catalytic performance differs between the post-hydrogen-reduction state (high activity and selectivity) and the post-oxidative stabilization state (low activity and selectivity).36 The results shown in Fig. 4b and c indicate that this performance decline cannot be explained by the changes in particle size or dispersion. Therefore, the loss of nanoparticle crystallinity, revealed in Fig. 4d–g, is likely the primary cause of catalyst deactivation. This is consistent with the in situ XAFS results, which showed that Ge in the nanoparticles was re-oxidized upon oxygen exposure, as evidenced by the increased Ge–O coordination number and decreased Ge–Re contributions. This Ge oxidation disrupts the atomic-level Re–Ge mixing essential for stabilizing the metallic Re(0) state, thereby explaining both the loss of crystallinity and the accompanying decline in catalytic performance. Thus, the ability to clearly decouple multiple structural factors affecting catalytic performance was made possible by the complementary approach of combining spatially averaged XAFS with direct atomic-resolution STEM imaging.
First, the elemental composition of the nanoparticles was investigated using STEM-EDS. Fig. 5a shows an elemental map overlaying Re M (red) and Ge L (green) maps (the individual maps for Re and Ge are shown in Fig. S11). The nanoparticles contain Re and Ge, as shown for the representative particles highlighted by the frames. Quantitative analysis of approximately 50 particles indicated that both Re and Ge were detected in every particles, with an average Ge/Re atomic ratio of 0.53 ± 0.25 (see the histogram in Fig. S12). This is consistent with the overall ratio of 0.51 obtained from the entire mapped area. The observed Ge/Re ratio is significantly lower than the nominal Ge/Re loading ratio of approximately 1.5, indicating that only a small fraction of the total loaded Ge is incorporated into the Re-Ge nanoparticle alloys. The remaining Ge most likely exists as highly dispersed GeOx compounds on the TiO2 support that are below the STEM-EDS detection limit. This interpretation is consistent with the substantial Ge–O coordination number detected by XAFS after hydrogen reduction. This Re-rich composition, combined with the in situ XAFS finding that both elements are in the metallic state, indicates that the nanoparticles are Re-based Re–Ge alloys. The formation of this Re-rich alloy is inconsistent with the bulk Re-Ge equilibrium phase diagram, where only a Ge-rich intermediate phase is considered thermodynamically stable.39,40 Therefore, the nanoparticles exist in a nanoscale-specific metastable phase.1
The crystal structure of the nanoparticles was investigated based on the HAADF-STEM images. Fig. 5b shows an image acquired from an area nearly identical to that shown in Fig. 5a, while Fig. 5c is a magnified image of a representative nanoparticle selected for analysis. This particle reflects the average characteristics of the nanoparticle ensemble in terms of both size (1.5 nm in diameter) and composition (Ge/Re = 0.47). The FFT pattern obtained from the image of this particle (Fig. 5c, inset) shows hexagonal reflections. This characteristic structure was commonly observed among the nanoparticles of this catalyst (Fig. S13). However, the FFT pattern did not match those of the hexagonal close-packed (HCP) structure of bulk Re,41 the diamond structure of Ge,41 or the theoretically calculated icosahedral structure of Re nanoparticles.30 Based on the EDS analysis indicating that this nanoparticle is Re-rich, we focused on Re-based structural models and first tested whether a distorted HCP Re structure could reproduce the experimental FFT pattern. We therefore examined anisotropically distorted HCP Re models by varying the a-axis strain (εa) and c-axis strain (εc); however, none of them could reproduce the experimental FFT patterns or the real-space atomic arrangement (Table S4 and Fig. S14). To explore alternative Re-based structural models that could explain the experimental results, DFT calculations were performed. The calculations revealed that, in addition to the previously reported icosahedral structure,30 the face-centered cubic fragment (FCCf) structure30,42 is also identified as a low-energy candidate (Fig. S15). Therefore, a hemispherical atomic model based on this FCCf structure was constructed, reflecting the experimentally observed particle size. Fig. 5d shows the atomic model and the corresponding simulated diffraction. This model showed good agreement with the experimental results (Fig. 5c) in both the real-space atomic arrangement (HAADF-STEM image) and the reciprocal-space diffraction pattern. It should be noted that the present DFT calculations were performed for free-standing Re clusters without explicitly considering the TiO2 support. Although support interactions may affect cluster stability and structure, such an analysis is beyond the scope of this study.
To elucidate the atomic mixing pattern of Re and Ge within this FCCf structure, the Z-contrast in the HAADF-STEM images was quantitatively analyzed.43–45 The experimental intensity profiles were compared with those simulated based on the major atomic ordering models for bimetallic nanoparticles, including core–shell, subcluster-segregated, ordered alloy and random alloy structures1,2 (see Fig. S16). The results demonstrated that the random alloy model, in which Re and Ge are atomically mixed, as depicted in Fig. 5d, best reproduced the experimental image. This finding was consistent with the spatially averaged XAFS results, which indicated that Re and Ge are present within the first coordination spheres of each other. These results directly demonstrate, through complementary spectroscopic analysis and atomically resolved STEM imaging, that the Re–Ge nanoparticles form random alloys with a unique FCCf crystal structure. Although the present results of this study do not allow a clear distinction between the individual contributions of electronic effects (the electron-rich Re(0) state) and geometric effects (the shortened bond distances, the FCCf structure and random alloy mixing), both are likely to play important roles in the enhanced catalytic performance.
| This journal is © The Royal Society of Chemistry 2026 |