Takuya
Suguro
,
Fuminao
Kishimoto
*,
Sota
Kuramoto
,
William J.
Movick
and
Kazuhiro
Takanabe
*
Department of Chemical System Engineering, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, Japan. E-mail: kfuminao@chemsys.t.u-tokyo.ac.jp; takanabe@chemsys.t.u-tokyo.ac.jp
First published on 5th July 2024
This study investigated the ammonia decomposition mechanism over Ru/CeO2. Isotopic tests using ND3 revealed that the rate-determining step involves adsorbed nitrogen atoms on Ru. Moreover, an inverse kinetic isotope effect where ND3 decomposition was faster than NH3 was clearly observed. The origin of the inverse effect was explained by the lower D coverage on the catalyst surface compared to H coverage for mitigating the inhibition of ND3 activation.
NH3 can be decomposed using various metals.6–11 Ru shows outstanding performance for NH3 decomposition under high NH3 concentrations and shows optimal dissociative N2 adsorption energy.8–12 The various types of support materials were investigated for enhancing the NH3 decomposition performance over the Ru catalyst, e.g. CeO210,11 (with PrOx13), MgAl2O4,14 Al2O3,15 carbon nanotubes,16–18 carbon nanofiber,19 SiC,20 C12A7:e− (ref. 21), and MgO8,9 with various exposed crystal facets.22
The differences in ammonia decomposition kinetics in terms of the exposed facets of Ru23–25 have been investigated. Egawa et al. demonstrated NH3 and ND3 isotope tests at ultra-low ammonia pressure (∼10−5 kPa) over single-crystal surfaces of Ru (1110) and Ru (001), which showed the normal kinetic isotope effect (NH3 decomposition occurred faster than ND3 decomposition).23 The rate-determining step was N–H bond cleavage under such low ammonia pressure conditions. Tsai et al. demonstrated a temperature dependence of the reaction rate of Ru (0001) and showed that the associative desorption of N2 was the rate-determining step below 377 °C. The rate-determining step changed to cleavage of the N–H bonding at higher temperature.24 A study using density functional theory predicted that the formation of the molecular nitrogen would be the rate-determining step with Ru (0001) and Ru (111).25 The Temkin–Pyzhev model is widely accepted to express the ammonia decomposition kinetics and assumes that associative N2 desorption step is the rate-determining step.15,26 This model fits well with conditions where reaction inhibition by hydrogen is observed at low temperatures or high pressures.27 However, several other models of reaction kinetics show different rate-determining steps.28–30 Armenise et al. reported a microkinetic analysis of ammonia decomposition with a Langmuir–Hinshelwood approach and modeled the reaction to predict an integral reactor by assuming the rate-determining step to be N2 desorption and strong hydrogen and ammonia adsorption.31
Although isotope experiments are a powerful tool for examining rate-determining steps, to our best knowledge, there are hardly any reports on isotope tests in the pressure range where reaction inhibition by adsorbed molecules is observed. Antunes et al. demonstrated the normal kinetic isotope effect over supported Ru/Al2O3 at an ammonia partial pressure of 1.5 Pa.32 However, different ammonia partial pressures should result in different catalyst surface coverage and thus different kinetics.12
Therefore, this paper reports the kinetic isotope effect of ammonia decomposition reaction using ND3 at 1 kPa over Ru/CeO2 (TEM image is shown in Fig. S1, ESI†). Fig. 1 shows the experimental setup used for NH3 decomposition. The catalyst bed was placed in a quartz tube with an internal diameter of 4 mm and heated with an electric furnace. The temperature of the cylindrical catalyst bed was measured using a thermocouple located outside the quartz tube at the catalyst bed. The catalytic reaction was performed with a closed circulation configuration (gradient-less batch reactor). The quantity of the accumulated NH3 gas was estimated from the infrared absorption spectra passing through the in-line gas cell.
Fig. 1 Schematic diagram of the closed recirculation batch system for catalytic ammonia decomposition. |
Fig. S2(a) (ESI†) shows typical results of the NH3 decomposition at various temperatures over a supported Ru/CeO2 catalyst. Deposition of colloidal Ru nanoparticles (mean particle size ∼4 nm) on the CeO2 support was employed to eliminate particle-size effects among the different samples. For tests at higher temperature, the Ru-deposited CeO2 catalyst was diluted with pristine CeO2 to control the reaction rate. As shown in Fig. S2(b) (ESI†), the NH3 decomposition rates were obtained from the slope of the line connecting the several plotted points in Fig. S2(a) (ESI†).
Fig. 2 shows the partial pressure dependence of ammonia decomposition over the supported Ru/CeO2 catalyst. The apparent kinetic order with respect to NH3 partial pressure was 0.23 at 200 °C, which gradually increased with the temperature, reaching 0.62 at 400 °C. Next, N2 was introduced into the initial gas composition to study its partial pressure dependence on the reaction rate. The N2 partial pressure was insensitive to the NH3 decomposition rate at both 200 and 400 °C (Fig. 2(b)).
Fig. 2 NH3 (red square) and N2 (blue circle) pressure dependence of ammonia decomposition rates (a) over 0.1 mg of 5 wt% Ru/CeO2 diluted with 9.9 mg of CeO2 within pellets at 200 °C and (b) over 1 mg pellets of 10× further dilution with CeO2 within the pellets at 400 °C. Initial conditions: ∼0.5 kPa NH3 balanced with Ar, total 101 kPa. Raw data are shown in Fig. S3 and S4 (ESI†). |
H2 introduction caused a decrease of the ammonia decomposition rate at 200–400 °C. Substantially negative order of −0.66 on H2 partial pressure was observed at 200 °C (Fig. 3). The inhibition effects by H2 became less severe with increasing temperature, reaching −0.35 at 400 °C.
Fig. 3 H2 pressure dependence of ammonia decomposition rates over 0.1 mg of 5 wt% Ru/CeO2 diluted with 9.9 mg of CeO2 within the pellets at 200 (blue square), 250 (light blue circle) and 300 °C (green diamond) and over 1 mg pellets of 10× further dilution with CeO2 within pellets at 350 (yellow pentagon) and 400 °C (orange pentagon) (5 kPa NH3, balanced by Ar). Raw data are shown in Fig. S5 (ESI†). |
Elementary reaction steps were examined for detailed analysis of the kinetics as follows:
(1) |
(2) |
(3) |
(4) |
(5) |
(6) |
(7) |
Fig. 4(a) shows the Arrhenius plot of NH3 and ND3 decomposition over the supported Ru/CeO2 catalyst. Notably, an inverse kinetic isotope effect was observed at all temperatures, where the decomposition rate of ND3 was greater than that of NH3. The apparent activation energy for NH3 and ND3 decomposition was almost identical (109 and 108 kJ mol−1). If N–H bond breaking (steps (2–4)) is involved in the rate-determining step, the activation energies should obviously be different. Therefore, the identical activation energy suggested that the rate-determining step should be steps (1), (5), (6) and (7), and most likely for associative N–N bonding formation (step (5)).
Fig. 4 (a) Arrhenius plot of NH3 (black square) and ND3 (red circle) decomposition over 0.1 mg of 5 wt% Ru/CeO2 diluted with 9.9 mg of CeO2 within pellets below 300 °C and over 1 mg pellets of 10× further dilution with CeO2 within pellets above 400 °C. The initial gas ratio was Ar:NH3 = 100:1 kPa or Ar/ND3 = 100:1 kPa. The raw data are shown in Fig. S6–S10 (ESI†). Error bars: 34.1% (1σ of Gaussian distribution). (b) Ratio of NH3 decomposition rate, rNH3, to ND3 decomposition rate, rND3, as a function of reaction temperature. Averaged rNH3/rND3 was 0.414 with a standard error of 0.062 (standard deviation = 0.139). |
Fig. 4(b) shows the ratio of NH3 and ND3 decomposition rates (rNH3/rND3) in the Arrhenius expression. The degree of kinetic isotope effect was independent of the reaction temperature and was almost constant (∼0.5). Aika and Ozaki studied the inverse kinetic isotope effect on NH3 synthesis (i.e., reverse reaction of NH3 decomposition). Strong inverse kinetic isotope effects were observed in the reaction of H2 and D2 with N2 over iron,33,34 molybdenum-based,35 or even Ru-based36,37 catalysts. The inverse kinetic isotope effects were explained by the larger equilibrium constant for ND3 decomposition relative to NH3 decomposition, which increased the empty sites for N2 activation.
Based on elementary steps (1)–(7), the overall reaction rate can be expressed as follows with an assumption of nitrogen triple bond formation in step (5) as the rate-determining step (adopted from Stoltze38).
(8) |
(9) |
(10) |
Assuming NH3 as the MARI, the order of the partial pressure dependence on H2 is still minus three. To be consistent with the zero order of N2 and the slightly negative order of H2 demonstrated by experiments, we derived the rate expression by assuming NH or NH2 species as the MARI, respectively:
(11) |
(12) |
From the experimental results, both the ammonia and hydrogen orders increased with increasing temperature. Typically, the surface reaction intermediates decreased as the temperature increases, meaning that the coverage term, θ*, becomes closer to zero at higher temperatures. Therefore, based on eqn (8), the dependence on ammonia partial pressure increases with increasing temperature under ammonia partial pressure. On the other hand, the effect of reaction inhibition by intentionally introduced hydrogen weakened as the temperature increased. This change can be explained by the change in MARI with increasing temperature. At high temperatures, MARI approaches N, and the apparent hydrogen partial pressure dependence is expected to approach zero, as shown in the equation below.
−rNH3 = k5+ | (13) |
From the rate expression as eqn (11) or (12), one possibility for explanation of the inverse kinetic isotope effect is the difference of the equilibrium constant of H2 or D2 adsorption/desorption (K7). Kellner and Bell studied H2/D2 isotope effects on Fischer–Tropsch synthesis over supported Ru catalysts.39 They calculated the equilibrium constant of H2 and D2 adsorption on Ru from vibrational frequencies for molecular and atomically adsorbed H2 and D2 (derivation is shown in the ESI†). The ratio of the equilibrium constants of D2 adsorption and that of H2 adsorption (K7,D/K7,H) was less than 1 between 180 °C and 270 °C, or the equilibrium shifted more toward desorption for D2 than for H2. Based on rate expression (11) or (12), lower K7 leads to a higher overall reaction rate due to more active surface species for the forward reaction (less inhibition). Therefore, the difference in the adsorption equilibrium constants of H2 and D2 can be one possible explanation for the inverse kinetic isotope effect in the NH3 decomposition reaction over the Ru/CeO2 catalyst observed in this study, suggesting that the secondary kinetic isotope effect appears inversely.
In conclusion, kinetic analysis using Ru/CeO2 revealed that comparable apparent activation energies between NH3 and ND3 decomposition reactions (∼108 kJ mol−1) were obtained with an inverse kinetic isotope effect (rNH3/rND3 ∼ 0.5). The kinetic orders were zero with respect to N2 partial pressure and negative with respect to H2 partial pressure on the NH3 decomposition rate. The results suggest that the rate-determining step does not involve N–H or N–D-bond dissociation, but it likely involves N–N triple bond formation, of which the rate is strongly perturbed by the surface H(D) adsorption equilibrium from H2 and D2 to the Ru surface. The adsorption equilibrium constant for D2 is smaller than that for H2, resulting in lower D coverage than H coverage, which is advantageous for high N* coverage. It can be concluded that the observed isotope effect was the secondary kinetic isotope effect inversely influenced by H(D) surface coverage. A more comprehensive understanding of the kinetics can be attained with further studies on isotope effects under diverse conditions.
This work was financially supported by the Mohammed bin Salman Center for Future Science and Technology for Saudi-Japan Vision 2030 at The University of Tokyo (MbSC2030), and the Science and Technology Research Partnership for Sustainable Development (SATREPS) in collaboration between the Japan Science and Technology Agency (JST, JPMJSA2104) and Japan International Cooperation Agency (JICA).
Footnote |
† Electronic supplementary information (ESI) available: Experimental methods, partial pressure dependence, and fitting result. See DOI: https://doi.org/10.1039/d4cc02372b |
This journal is © The Royal Society of Chemistry 2024 |