Abstract
The large-scale conversion of N2 and H2 into NH3 (refs. 1,2) over Fe and Ru catalysts3 for fertilizer production occurs through the Haber–Bosch process, which has been considered the most important scientific invention of the twentieth century4. The active component of the catalyst enabling the conversion was variously considered to be the oxide5, nitride2, metallic phase or surface nitride6, and the rate-limiting step has been associated with N2 dissociation7,8,9, reaction of the adsorbed nitrogen10 and also NH3 desorption11. This range of views reflects that the Haber–Bosch process operates at high temperatures and pressures, whereas surface-sensitive techniques that might differentiate between different mechanistic proposals require vacuum conditions. Mechanistic studies have accordingly long been limited to theoretical calculations12. Here we use X-ray photoelectron spectroscopy—capable of revealing the chemical state of catalytic surfaces and recently adapted to operando investigations13 of methanol14 and Fischer–Tropsch synthesis15—to determine the surface composition of Fe and Ru catalysts during NH3 production at pressures up to 1 bar and temperatures as high as 723 K. We find that, although flat and stepped Fe surfaces and Ru single-crystal surfaces all remain metallic, the latter are almost adsorbate free, whereas Fe catalysts retain a small amount of adsorbed N and develop at lower temperatures high amine (NHx) coverages on the stepped surfaces. These observations indicate that the rate-limiting step on Ru is always N2 dissociation. On Fe catalysts, by contrast and as predicted by theory16, hydrogenation of adsorbed N atoms is less efficient to the extent that the rate-limiting step switches following temperature lowering from N2 dissociation to the hydrogenation of surface species.
Similar content being viewed by others
Main
Figure 1a shows how surface-sensitive operando X-ray photoelectron spectroscopy (XPS) is measured together with reaction-product detection during the Haber–Bosch process in the POLARIS instrument13. XPS is a powerful technique for investigating the chemical state of catalytic surfaces through core-level shifts that traditionally required vacuum conditions, but operando studies can be conducted using a differential pumping scheme17. The Fe and Ru single-crystal surfaces are mounted in front of the electron spectrometer with a gap of 30 µm and gases are fed through the front cone of the electron lens, creating a localized virtual catalytic reactor of elevated pressure with a rapid gas flow13. The typical operational pressure for ammonia synthesis is 50–200 bar (ref. 18), at which the gas-phase equilibrium is strongly shifted towards the product, giving a high final conversion to ammonia. However, during the initial phase of the Haber–Bosch process, when not much ammonia has yet been produced, the reaction also proceeds with a high rate at our operational pressures of up to 1 bar (refs. 19,20).
The incoming X-rays were set to an energy of 4,600 eV and the incidence at an angle below total reflection, allowing for high surface sensitivity despite high kinetic energy electron detection. The emitted photoelectrons will pass into the spectrometer through orifices in the front cone and be detected in a hemispherical analyser. The inset in Fig. 1a shows an example of an N1s spectrum of 1:3 N2:H2 gases at 1 bar at 673 K, indicating NH3 (blue), NH2 (purple), NH (red), surface N (green) and nitride surface (yellow) components. The measurements were conducted at a photon flux at which no detectable X-ray-beam-induced changes could be seen during the Haber–Bosch process (see Methods for further details).
To track the production of NH3, masses 15 and 16 were monitored in the mass spectrometer (see Methods), as shown in Fig. 1b. The relative chemical reactivities shown in Fig. 1c were determined by measuring the mass spectrometer ammonia signal with respect to the signal of all constituents to compute the number of ammonia molecules formed per second per surface site, which is then further normalized to the highest activity shown by any surface at any temperature (see Methods for further details). The reaction rate increases with increasing temperature and is higher for the stepped Fe(210) than the flat Fe(110) surface, in agreement with previous high-pressure-reactor studies9. The highest rate is seen for the \({\rm{Ru}}(10\bar{1}3)\) surface, as expected based on polycrystalline studies showing that Ru has higher activity than Fe (ref. 21). The maximum rate for Ru is not at the highest temperature of 723 K, as for the Fe surfaces, but at 623 K, also in accordance with catalytic-reactor studies22.
On exposure to pure N2 gas at 150 mbar, the two Fe surfaces have a delayed but eventually rapid increase in the N1s intensity, showing bulk nitride formation (Fig. 2a,b). On the basis of the binding-energy position of the N1s peaks in the spectra, this corresponds to the formation of γ′-nitride and ε-nitride plus some small amount of chemisorbed N atoms on the bare Fe surface (see Extended Data Table 1). The nitride formation is more rapid on the Fe(210) surface, specifically the γ′-nitride, whereas on the Fe(110) surface, there is an equal amount of the two nitrides and slower growth. The thicknesses of the nitride layers are greater than ten monolayers; exact quantification depends on the reaction time, as the surface continues to evolve even after hours of observation (see Methods for details on monolayer calculations). We attribute the faster growth on the Fe(210) facet to the higher probability of N2 dissociation on the stepped surface23. At temperatures below 523 K, no nitride formation is observed.
The \({\rm{Ru}}(10\bar{1}3)\) reacts completely differently. Almost instantaneously after N2 exposure, the N1s intensity saturates and remains constant, corresponding to a coverage of 5% of a monolayer, and there is no bulk nitride formation at 623 K (Fig. 2c). The coverage is comparable with previous work, which predicts 17% of a monolayer at 500 K and a pressure of 100 mbar (ref. 23). The small amount of N2 on the Ru surface indicates a much weaker N–metal interaction than on Fe, as expected from theoretical predictions16. The two components are at 397.4 eV and 397.9 eV, and we tentatively assign these to N adsorbed on terraces and steps, respectively (Extended Data Fig. 1). It is interesting that a weak, broad feature is seen at approximately 399–400 eV, with a binding energy consistent with adsorbed N2 (ref. 24); see Extended Data Fig. 1.
When the pure N2 gas is replaced by 1:1 N2:H2 at 300 mbar, a marked change on the two Fe surfaces occurs within the first spectral sweep (90 s), shown at the bottom of Fig. 2a,b. The nitrides instantaneously disappear and only a small amount of adsorbed N atoms with a coverage of 2% of a monolayer on Fe(110) and 5% on Fe(210) remains. At the same time as the gas mixture is introduced, NH3 is detected by the mass spectrometer. The rapid removal of the nitrides shows the strong reduction ability of the H2. The slow growth of nitrides (10–15 min) compared with the fast reduction (<100 ms) shows the difference in rates of N2 and H2 dissociation. The adsorbed N atom coverage is also substantially lowered on the \({\rm{Ru}}(10\bar{1}3)\) surface following the introduction of the 1:1 N2:H2 mixture at 300 mbar and decreases from 5% to <0.05% of a monolayer as NH3 is produced.
Next, we address the question of oxides potentially not being reduced on Fe under operando conditions owing to trace contaminations of water or CO2 in the gas phase5. Iron is known to oxidize in trace amounts of water or CO2 at room temperature, yet iron oxide is not readily reduced below 500 K and, as a result, even under pure hydrogen, iron will oxidize with high flows (see Methods for a detailed description). Figure 3 shows data collected at 500 mbar, 1:3 N2:H2 and various temperatures. The Fe 2p2/3 peaks in Fig. 3a from metallic iron at 706.5 eV and 707.4 eV are split owing to exchange interactions with the ferromagnetic valence electrons, and there is a broad Fe oxide peak at 710.8–709.8 eV, indicated by the grey rectangle. The Fe(110) sample is fully reduced as the temperature reaches 523 K at 500 mbar and the Fe(210) surface requires a higher temperature of 573 K, as seen in Fig. 3b. Fe(210) needs a higher temperature because of the stronger binding of oxygen on a stepped surface. Ru is metallic at all conditions. All surfaces are in a metallic state during the Haber–Bosch process, as expected because of the high concentration of adsorbed hydrogen (Fig. 3c). Note that these measurements were gathered simultaneously with the data in Fig. 4.
The adsorbed nitrogen species can be measured operando as NH3 is produced. First, focusing on the two Fe single-crystal surfaces (Fig. 4a,b), we observe only adsorbed N atoms on the surface at a binding energy of 397.4 eV, consistent with previous surface-science vacuum experiments once the recoil effect of the emitted atoms is considered (see Extended Data Table 1). Adsorbed molecular N2 could not be detected and would have been observed at 399.0, 401.2 or 405.9 eV (Extended Data Table 1), depending on the adsorption site and bonding type. The coverage of adsorbed N is 1.3% at 200 mbar and 0.6% at 500 mbar on the Fe(110) surface and increases on the Fe(210) surface to 5.0% and 1.5%, respectively. The higher coverage on the stepped surface is related to availability and stronger bonding of undercoordinated sites16. What is most surprising is that the coverage is not increasing at higher pressures; on the contrary, the coverage decreases slightly with increased pressure. Inspecting the N1s spectra in Fig. 4d, measured at 1 bar and 673 K, the peak is barely distinguishable from the noise, implying an even lower coverage. It would be tempting to expect an increase in N coverage with increasing pressure because the impinging rate of N2 molecules increases, but obviously also does the rate of H adsorption. Although we cannot determine the H coverage with XPS, our data suggest that the hydrogenation ability of the surface increases with the total pressure; this would explain a more efficient further reaction of the adsorbed N atoms. Extrapolating to much higher pressures, we predict that the Fe surface is an almost pristine metal under realistic conditions. The fact that no amines (NH or NH2) or NH3 are observed at the reaction temperature of 673 K indicates that the rate-limiting step after N2 dissociation is the hydrogenation of adsorbed N, and the rates of the other hydrogenation steps of NH and NH2 as well as NH3 desorption are much faster. At high temperatures, the Ru surface (Fig. 4g) has adsorbed N at 397.4 eV and the adsorbate coverage is almost negligible, with <0.1% of a monolayer of both NH and NH2 species, independent of pressure within the noise limit. Here the surface is almost entirely clean of any species at conditions of high reaction rate.
At 523 K, for which the reaction proceeds very slowly, the population of the adsorbates changes. There is a slight increase of the adsorbed N on Fe(110) at 500 mbar to 2.3% of a monolayer (Fig. 4d). The Fe(210) surface shows large differences compared with the higher-temperature spectra (Fig. 4e). Further peaks at 398.0 eV, 398.9 eV and 400.2 eV formed, corresponding to NH, NH2 and NH3, as determined by previous XPS vacuum studies9,25,26 and calculated relative peak positions (Extended Data Table 1). Note that the peak at 399 eV is not related to adsorbed N2 because ex situ XPS studies observed the peak when the Fe catalyst was cooled down to room temperature in the reaction mixture and moved to a vacuum, in which all molecular N2 would desorb. We observe a relatively high coverage of NH2 (24.8%), adsorbed N (4.3%), NH (6.7%) and NH3 (5.2%) at 200 mbar. There is a slight pressure dependence, for which—in particular—the NH2 decreases to 9.3%. Clearly, there exist conditions in which the adsorbed N and NHx species are strongly adsorbed on step sites owing to a substantially lower hydrogenation rate. Decreasing the temperature further to 423 K, adsorbed NHx and NH3 become visible on the Fe(110) surface. These trends are seen across 423 to 623 K (Extended Data Fig. 2).
On Ru at 523 K at 500 mbar (Fig. 4c), we still see very low coverages, although the coverage of adsorbed N at steps has increased to 0.5%, as well as adsorbed NH2 to 0.1% and adsorbed NH3 to 0.1% at around 400 eV. The NH signal increases with pressure, but the nitrogen coverage quantification of these results is nearly within the margin of error. If there is an increase in coverage with pressure for Ru, it may indicate that the H2–metal interaction for Ru is weaker than for Fe, possibly leading to higher coverages at operational pressures. The adsorbed N species is much more reactive on Ru than Fe, supporting previous theoretical predictions16.
We can discriminate the various proposed hypotheses and put forward ideas consistent with the data on the chemical state of the catalysts and reaction mechanism in terms of rate-limiting steps. We have shown that nitride formation is far slower than nitride reduction and that the surface states are all metallic with low coverages of atomic nitrogen. There is no evidence for interstitial nitrogen, oxides or high coverage of any species of nitrogen, especially over the most active catalysts. It is interesting to compare the hydrogenation reactions of CO and N2, which are isoelectronic molecules. In the case of the Fischer–Tropsch reaction on Fe(110), a thick carbide is formed15, whereas in the Haber–Bosch process, on the same surface, only a pristine metallic phase is generated. Clearly, the difference in the bond breaking of the CO molecule with respect to N2 and the strength of the adsorbed C and N play an essential role.
The different reaction steps in NH3 synthesis have been proposed as the following10:
in which * means surface species and θ* indicate empty sites available for bonding.
The simplest case is the \({\rm{Ru}}(10\bar{1}3)\) surface, for which we can directly explain that steps 3–6 are extremely rapid with no build-up of intermediates, pointing to 1 and 2 as the rate-limiting steps. We observe that the population of adsorbed N2 is extremely low at high temperatures. The adsorbed molecular state is indeed observed at the low reaction temperature of 523 K, at which its dissociation limits the reaction. We conclude that the rate-limiting step of NH3 production is the dissociation of the adsorbed N2, fully in line with theoretical estimations12. Even at low temperatures, the surface is mostly adsorbate free, with little adsorbed NHx seen, because of the strong bonding to step sites in comparison with terrace atoms23. Although we have not observed definitive pressure dependence in the population of adsorbed N, it is plausible that the step sites will become more populated but are expected to remain well below a monolayer.
On Fe it is well established that the rate-limited steps is the molecular dissociation7,8,9, supported by the correlation between the NH3 production rate and the N2 dissociative sticking coefficient for the different single-crystal surface facets9,27. However, the results here show that, at all temperatures, a factor of around 100 times higher population of adsorbates is observed in comparison with the stepped Ru surface and we can no longer postulate that the reaction proceeds with a high rate after the molecular dissociative steps. Furthermore, there are no signs of molecularly adsorbed N2 even at the lowest temperatures, indicative of a much higher rate of step 1b. Above 573 K, we observe adsorbed N that is more populated on the stepped crystal, indicating that the hydrogenation step 3 also partly controls the rate12.
The coverage of N species on the Fe surfaces decreases with increasing total pressure at a constant N2:H2 ratio, implying that the N2 dissociation step is slower than the hydrogenation step10. Most likely, the coverage of adsorbed H increases with pressure, resulting in faster hydrogenation. Because the coverage of H2 at the reaction temperatures is expected to be low, we can assume that there is no inhibition of N2 dissociation caused by the adsorbed hydrogen27.
The population of intermediates shows that, as the reaction temperature lowers, the rate-limiting step switches to become hydrogenation of N, NH and NH2 as well as NH3 desorption (steps 3–6), demonstrating differences in the bonding at different high and low coordinated Fe sites. This agrees with earlier observations of the activation energy for hydrogenation being much higher than for N2 dissociation10 and the difference in the barriers of these two steps thus becoming prominent at low temperatures: although the N2 dissociation rate at high temperatures is low owing to a low sticking coefficient that limits N2 adsorption10, we see a large population of amines NHx and NH3 on Fe at low temperatures. This trend, not seen with Ru, points to the hydrogenation steps affecting the overall rate on Fe. At higher pressures at which more N2 is converted and the NH3 content is higher, the back reaction may become important. Indeed, for Ru, it has been theoretically predicted that the coverage of nitrogen species may become substantially higher28.
In closing, we note that, although concerns over the environmental impact of ammonia synthesis have spurned interest in low-pressure alternatives and these might indeed be feasible29, the Haber–Bosch process looks set to remain the primary method of ammonia production for many years to come. A better understanding of the mechanism at play might help to further improve the efficiency and, thereby, lower the environmental impact of this important industrial process. We anticipate that our approach to operando studies will contribute to this endeavour, by making it possible to explore the surface chemistry associated with ammonia formation in the presence of promotors and by making it possible, once measurements at higher pressures and with a higher NH3 content are feasible, to explore the impact of the ammonia decomposition back reaction.
Methods
Ambient pressure XPS
All XPS data are collected with total external reflection X-rays and normalized to core levels of the substrate. The emitted photoelectrons and gases pass into the spectrometer through orifices in the front cone to be detected in a hemispherical analyser. The overall resolution in the measurement was 0.2 eV, all spectra were normalized with respect to the Fe 2p3/2 or Ru 3d5/2 core levels unless stated otherwise and spectra are presented in counts per second (c.p.s.).
POLARIS is an ambient pressure X-ray photoelectron spectrometer that operates with several key differences from typical ambient pressure XPS systems. The sample is approached to 30 µm from a set of roughly 20-µm-diameter apertures that lead to the analyser; the X-rays used are in the tender range 4.6 keV for all data collected. Most importantly, though, is that the gas is delivered through the front cone directly to the sample, making a virtual pressure cell in which only the sample and aperture to the analyser are pressurized. To achieve surface sensitivity, grazing incidence X-rays are used within the total external reflection range 0.3° for iron and 0.45° for ruthenium. This geometry allows for surface sensitivity despite high kinetic energy electron detection; the probe depths are 15.5 Å and 12.6 Å for Fe and Ru, respectively30. The electron spectrometer is a HiPP-2 hemispherical analyser manufactured by Scienta Omicron; see ref. 13 for more details. The single-crystal samples (Surface Preparation Laboratory, 99.99% purity) are mounted in a steel sample holder and heated from the back side with a resistive heater. The temperature of the sample is measured with a type C thermocouple pressed between the sample and the heater. The separation between the sample and the apertures is held constant by PID feedback based on the pressure over the sample31. A Si(311) double-crystal monochromator was used, yielding a photon-energy bandwidth of approximately 130 meV, a 0.8 mm curved entry slit and 100 pass energy was used in the electron analyser.
Extended Data Fig. 1a shows an example of an N1s spectrum of 1:3 N2:H2 gases at 200 mbar and the sample at 423 K, indicating adsorbed N atoms on the surface. The measurements were conducted at a photon flux at which no detectable X-ray-beam-induced changes could be seen. Individual spectra were gathered for 30 to 300 min with no decreeable spectral changes when hydrogen was present in the gas phase. Extended Data Fig. 3 shows an example time interval of 2 h over Fe(210) at 423 K and 500 mbar in 1:3 N2:H2 gas mixture. Extended Data Fig. 3a,b shows the data over this time for mass fragments 15, 16, 17 and 18, with and without processing. Extended Data Fig. 3c shows the XPS spectra evolution with time and Extended Data Fig. 3d shows the time-averaged results. From these, it is clear that the only change observed with time is the decrease in water signal owing to the slow improvement of vacuum conditions under constant hydrogen conditions.
XPS data processing
All presented spectra are scaled by the number of sweeps and dwell time per data point. Further scaling is done based on the relative cross-section of the materials, as mentioned in the main text. To fit the spectra, CasaXPS was used with linear or Shirley backgrounds as needed. Peaks were fitted with modified Voigt function (LA) line shapes, which allows for asymmetry. Asymmetry was tuned for each component.
Sample preparation
Sample cleaning was performed by ion sputtering with 3 keV Ar+ for iron and 1 keV for ruthenium. The samples were annealed to 900 K for Fe and 1,100 K for Ru. Chemical cleaning was performed as needed by exposing the sample to either hydrogen or oxygen at elevated temperatures to remove oxygen or carbon, respectively. Small contaminations of sulfur and silicon were present, but the atomic composition was maintained at or below 1%.
Coverage
To model the coverage of the surfaces, the Fe4N and RuN nitrides were used as the physical representation of the surface species. Although not a perfect model, reference data of commercial nitrides verify that the surface constituents are similar in atomic bonding. We used the method previously established30 and typical XPS coverage formulation32 to calculate the coverages. Elemental cross-section data was taken from ref. 33.
To calculate the probe depth, the X-ray and electron mean free path need to be combined; this is done by calculating the X-ray field in the material34 at a given angle and using the TTP2M electron mean free path35 to determine the electron-signal intensity as a function of depth within the sample. Then the integral is evaluated over all depths to determine the effective probe depth. Once the probe depth is determined, the coverage is then calculated on the basis of the ratio of substrate to surface species intensity-weighted by the cross-sections and atomic densities32.
Mass spectrometry
To determine the amount of NH3 formed by the catalyst, a differentially pumped mass spectrometer (Hiden HAL/3F RC 301 PIC system) was attached to the first differential pumping stage of the XPS analyser. By leaking a small amount of gas from the pumping stage to the mass spectrometer, the composition of the gas over the sample was determined. To ascertain the gas composition, mass fragments of all relevant peaks were monitored. Impurities in the N2 and H2 gases were predominantly H2O and CO2. These contributions to the NH3 fragments were subtracted on the basis of the measured ratio of pure gas to contaminate. Owing to the marked overlap of water and NH3 ionization patterns, m/z = 15 and 16 were used as the markers for NH3. Further smoothing is done with a third-order Savitzky–Golay filter over a window of 1 s. The result of this analysis is shown in Fig. 1. As the mass spectrometer is highly sensitive, there is signal before any experiment from the chamber at all masses, including masses 15 and 16, most likely because of hydrocarbons. With the high-flow conditions required to establish the pressure for the XPS measurements, the amount of ammonia in the gas stream into the mass spectrometer is small and the signal becomes noisy. Therefore, to make a more accurate measurement of the NH3 production, time integration was done between the background level in pure N2 and that of pure H2. The background was subtracted from the time integration during ammonia production. Extended Data Fig. 4a shows an example mass spectrometer time trace in which there is negligible NH3 production, showing how the background change with gas flow is within the noise of the measurement and therefore requires time integration. Owing to the specific design of the high-flow virtual cell, unwanted gas molecules originating from reactions of the sample holder or heater cannot reach the single-crystal surface area that is examined by the opening into the electron spectrometer. Thereby, all measurement conditions are constant. The increase in ammonia production at higher temperatures is as expected according to refs. 9,21,22, providing further confidence that ammonia is produced.
The relative chemical activity (RCA) was calculated using the following equation. Time-averaged NH3% was calculated from the amount of signal from ammonia as described above per total signal from the mass spectrometer. Volume (V) of gas is the total volume of gas used during the measurement, pressure (P) over the sample, temperature (T) of the sample, the gas constant (R), time is the duration of the time when ammonia could have been produced, An is Avogadro’s number and sites is the number of active sites under the high-pressure area. Finally, the highest activity on any surface is a normalization to the maximum of any surface. The normalization is to account for systematic errors, such as the fact that most of the volume of gas used does not pass over the sample or the fact that not all sites in the high-pressure region under the front cone would be examined by the mass spectrometer.
With an instrument exposed to many gases over the years, there are signals at all masses, including masses 15 and 16, before any ammonia-synthesis experiment is performed, owing to desorption from the chamber walls. This desorption in the first differential pumping stage most likely comes from hydrocarbons. In the mass spectrometer, it is possible for crosstalk between channels or other instrumental errors to affect the signal. This is particularly true when the signal is very near the noise level, as in the work presented herein. Extended Data Fig. 3 shows the masses 15, 16, 17 and 18. Mass 17, corresponding to ammonia, is strongly affected by water production from H2 interaction with the chamber walls and the mass spectrometer itself, making quantitative analysis impossible. To decrease the possibility of the ammonia signal originating from instrumental errors, both masses 16 and 15 are included in the signal of ammonia. As discussed above, mass 17 is not included because of the large water signal. Extended Data Fig. 3a,b shows the effect of both processing and long acquisition times. The sample is Fe(210), 423 K, 500 mbar, 1:3 N2:H2 ratio. Here we can see that the atomic mass units of both 15 and 16 are constantly above the background signal from hydrocarbons or water. Meanwhile, mass 17 is not, owing to the strong overlap of OH and NH3 masses. Note that the background signal removed at this point in the processing does not account for all of the background signals. As described above, to determine the relative chemical activity, the signal of ammonia (masses 15 and 16) above the background signal in either pure N2 or pure H2 is taken. That ammonia signal is then compared with the total signal in the mass spectrometer over the same time period. By this method, the plotted data do not remove all of the background signals, yet when the data are processed for relative chemical activity, the entirety of the background is removed.
The error of the relative chemical activity is estimated on the basis of the signal-to-noise ratio of the background of the ammonia signal. Part of the calculation is to subtract the background, seen in Fig. 1b, between times 15 and 23 min; the fluctuations in the background can have a notable effect on the calculation of ammonia content. To ascertain the estimated error, the 95% confidence interval of the noise average and standard deviation over the collected time were introduced as an error source in the equation for relative chemical activity. Because the background signal and noise are similar for all experiments, the estimated error introduced is also similar. The relative chemical activity is meant to be a semiquantitative description of the abundance of ammonia, only as a comparative description of these similar systems, and to demonstrate that the trends follow previous more absolute activity measurements. Extended Data Fig. 4a shows an example of when there is no ammonia production at the lowest temperature with the least active catalyst, Fe(110) at 523 K (300 mbar, 1:1 ratio). Here we can see that extremely small ammonia production occurs and this is most likely the background level. Extended Data Fig. 4b, by contrast, shows the same surface and experiment at the higher temperature of 673 K and clearly shows that ammonia production increases with increasing temperature.
Beam effects
To determine the effect of the X-ray beam intensity on the observed species, two types of beam damage test were performed, attenuation and dark spectra, both carried out at 423 K, 500 mbar and a 3:1 mixture of H2 and N2. The dark spectra were performed by aligning the sample and gathering highly attenuated spectra, then closing off the light, cleaning the sample in hydrogen, then opening the shutter and obtaining a new spectrum. The result showed no change owing to the amount of time the sample was exposed to X-rays. The second test was done by gradually increasing the intensity of the beam to determine whether any beam damage would accumulate; here no change in the spectra once normalized to the attenuation factor was observed.
Nitrogen spectra gathered outside reaction conditions (that is, without hydrogen), such as those shown in Fig. 2, show not only chemical changes inherent to the reaction but also accumulation of beam-induced effects. The build-up of beam-induced nitride formation is slower than chemical activity but is not possible to fully avoid. For this reason, no attempt is made to quantify the formation rate of the various nitrides. The main finding of the paper is the distinct lack of nitrogen on the surface during the reaction and the slow nitride formation compared with fast reduction caused by hydrogen. The beam-induced nitride formation only serves to increase nitrogen formation. Therefore, the beam effects do not alter any conclusions of the activity of N2 compared with H2 but rather strengthen the finding.
To determine whether the beam has any effect on the mass spectrometer findings, two experiments of gas switching (in which the sample was exposed to pure N2 then 1:1 N2:H2 then pure H2) were performed with X-ray light and without. Extended Data Fig. 4b,c shows the mass fragments of all relevant species for these experiments. For these data, the sample was Fe(110) at 673 K. Extended Data Fig. 4b shows the measurement with beam, whereas the measurement in Extended Data Fig. 4c was collected without beam. Although there are some minor differences in the ammonia signal between the two datasets, none of the changes are what would be expected from beam effects. It is expected that beam effects in the mass spectrometer would strictly cause an increase or decrease in ammonia signal. The change in the relative chemical activity between the test done with and without beam is approximately 2% and well within the error of the experiment. From this, it is clear that the X-ray beam does not have any effect on the mass spectrometer findings.
Extended Data Fig. 5a shows the effect of photon flux on the nitrogen content. Note that the line of best fit shown in grey has a forced intercept to zero. Extended Data Fig. 5b shows the effect of flux on nitrogen speciation, showing no change to the components of the N1s spectra. These spectra were gathered at 500 mbar in a 1:3 gas ratio at 423 K over Ru, the conditions expected to be the most sensitive to beam damage. Equivalent tests were performed for all catalysts.
Trace contaminations
Although highly pure gas (5N for nitrogen and hydrogen) was used with in-line chemical purifiers (N2, model no. MC45-804; H2 gas, model no. MC45-904, SAES Group), trace impurities are still present; on the basis of the mass spectrometer data, approximately 6 ppm of water and 3 ppm of CO2. Owing to the high flows used, these small contaminants can react and build on the surface instantly. Furthermore, molecules will readily react with iron to form iron oxides; the same is not true for ruthenium. As a result, the iron surface at low temperatures will form a partial oxide but, as the temperature increases, the reduction by hydrogen outpaces the oxidation of the contaminants, yielding a metallic surface at relevant conditions.
Nitrogen species binding energies
Extended Data Table 1 presents previously published and computed binding energies for various amine and nitrogen species over iron and ruthenium. Note that, for comparison with studies at low photon energy, the recoil effect, in which the energy of the electron causes the nucleus to recoil, thereby decreasing the kinetic energy of the emitted electrons, needs to be considered36.
Data availability
Experimental data were generated at the PETRA III facility at the DESY Research Centre of the Helmholtz Association. Raw datasets are available from the corresponding authors on reasonable request.
References
Haber, F. & le Rossignol, R. Über die technische Darstellung von Ammoniak aus den Elementen. Z. Elektrochem. Angew. Phys. Chem. 19, 53–72 (1913).
Ertl, G. The arduous way to the Haber–Bosch process. Z. Anorg. Allg. Chem. 638, 487–489 (2012).
Mittasch, A. & Frankenburg, W. Early studies of multicomponent catalysts. Adv. Catal. 2, 81–104 (1950).
Smil, V. Detonator of the population explosion. Nature 400, 415–415 (1999).
Arabczyk, W., Jasińska, I. & Lubkowski, K. The surface properties of iron catalyst for ammonia synthesis. React. Kinet. Catal. Lett. 83, 385–392 (2004).
Kandemir, T., Schuster, M. E., Senyshyn, A., Behrens, M. & Schlögl, R. The Haber–Bosch process revisited: on the real structure and stability of “ammonia iron” under working conditions. Angew. Chem. Int. Ed. 52, 12723–12726 (2013).
Stoltze, P. & Nørskov, J. K. Bridging the ‘pressure gap’ between ultrahigh-vacuum surface physics and high-pressure catalysis. Phys. Rev. Lett. 55, 2502 (1985).
Ertl, G., Lee, S. B. & Weiss, M. Adsorption of nitrogen on potassium promoted Fe(111) and (100) surfaces. Surf. Sci. 114, 527–545 (1982).
Somorjai, G. A. & Materer, N. Surface structures in ammonia synthesis. Top. Catal. 1, 215–231 (1994).
Ertl, G. Surface science and catalysis—studies on the mechanism of ammonia synthesis: the P. H. Emmett Award address. Catal. Rev. 21, 201–223 (2006).
Spencer, M. S. On the rate-determining step and the role of potassium in the catalytic synthesis of ammonia. Catal. Lett. 13, 45–53 (1992).
Honkala, K. et al. Ammonia synthesis from first-principles calculations. Science 307, 555–558 (2005).
Amann, P. et al. A high-pressure x-ray photoelectron spectroscopy instrument for studies of industrially relevant catalytic reactions at pressures of several bars. Rev. Sci. Instrum. 90, 103102 (2019).
Amann, P. et al. The state of zinc in methanol synthesis over a Zn/ZnO/Cu(211) model catalyst. Science 376, 603–608 (2022).
Shipilin, M. et al. In situ surface-sensitive investigation of multiple carbon phases on Fe(110) in the Fischer–Tropsch synthesis. ACS Catal. 12, 7609–7621 (2022).
Medford, A. J. et al. Assessing the reliability of calculated catalytic ammonia synthesis rates. Science 345, 197–200 (2014).
Salmeron, M. & Schlögl, R. Ambient pressure photoelectron spectroscopy: a new tool for surface science and nanotechnology. Surf. Sci. Rep. 63, 169–199 (2008).
Humphreys, J., Lan, R. & Tao, S. Development and recent progress on ammonia synthesis catalysts for Haber–Bosch process. Adv. Energy Sustain. Res. 2, 2000043 (2021).
Ertl, G. & Thiele, N. XPS studies with ammonia synthesis catalysts. Appl. Surf. Sci. 3, 99–112 (1979).
Dahl, S., Logadottir, A., Jacobsen, C. J. H. & Norskov, J. K. Electronic factors in catalysis: the volcano curve and the effect of promotion in catalytic ammonia synthesis. Appl. Catal. A Gen. 222, 19–29 (2001).
Saadatjou, N., Jafari, A. & Sahebdelfar, S. Ruthenium nanocatalysts for ammonia synthesis: a review. Chem. Eng. Commun. 202, 420–448 (2014).
Pan, C., Li, Y., Jiang, W. & Liu, H. Effects of reaction conditions on performance of Ru catalyst and iron catalyst for ammonia synthesis. Chin. J. Chem. Eng. 19, 273–277 (2011).
Dahl, S. et al. Role of steps in N2 activation on Ru(0001). Phys. Rev. Lett. 83, 1814 (1999).
Umbach, E. Satellite structures in photoemission spectra from different types of absorbates. Surf. Sci. 117, 482–502 (1982).
Grunze, M. et al. π-bonded N2 on Fe(111): the precursor for dissociation. Phys. Rev. Lett. 53, 850 (1984).
Kishi, K. & Roberts, M. W. Adsorption of nitrogen and ammonia by polycrystalline iron surfaces in the temperature range 80–290 K studied by electron spectroscopy. Surf. Sci. 62, 252–266 (1977).
Ertl, G. Primary steps in catalytic synthesis of ammonia. J. Vac. Sci. Technol. A 1, 1247–1253 (1983).
Logadóttir, Á. & Nørskov, J. K. Ammonia synthesis over a Ru(0001) surface studied by density functional calculations. J. Catal. 220, 273–279 (2003).
Vojvodic, A. et al. Exploring the limits: a low-pressure, low-temperature Haber–Bosch process. Chem. Phys. Lett. 598, 108–112 (2014).
Goodwin, C. M. et al. The structure of the active Pd state during catalytic carbon monoxide oxidization. J. Phys. Chem. Lett. 12, 4461–4465 (2021).
Goodwin, C. M. et al. A novel method to maintain the sample position and pressure in differentially pumped systems below the resolution limit of optical microscopy techniques. Appl. Spectrosc. 75, 137–144 (2021).
Shard, A. G. Practical guides for x-ray photoelectron spectroscopy: quantitative XPS. J. Vac. Sci. Technol. A 38, 041201 (2020).
Trzhaskovskaya, M. B. & Yarzhemsky, V. G. Dirac–Fock photoionization parameters for HAXPES applications. At. Data Nucl. Data Tables 119, 99–174 (2018).
Stepanov, S. A. X-ray server: an online resource for simulations of x-ray diffraction and scattering. Proc. SPIE 5536, 16–26 (2004).
Shinotsuka, H., Tanuma, S., Powell, C. J. & Penn, D. R. Calculations of electron inelastic mean free paths. X. Data for 41 elemental solids over the 50 eV to 200 keV range with the relativistic full Penn algorithm. Surf. Interface Anal. 47, 871–888 (2015).
Takata, Y. et al. Recoil effect of photoelectrons in the Fermi edge of simple metals. Phys. Rev. Lett. 101, 137601 (2008).
Jiang, Y., Zhang, X., al Mehedi, A., Yang, M. & Wang, J. P. A method to evaluate α″-Fe16N2 volume ratio in FeN bulk material by XPS. Mater. Res. Express 2, 116103 (2015).
Moreno-Armenta, M. G., Diaz, J., Martinez-Ruiz, A. & Soto, G. Synthesis of cubic ruthenium nitride by reactive pulsed laser ablation. J. Phys. Chem. Solids 68, 1989–1994 (2007).
Gröning, P., Nowak, S. & Schlapbach, L. Surface modifications of nitrogen-plasma-treated stainless steels. Appl. Surf. Sci. 64, 265–273 (1993).
Ogata, Y., Aika, K. I. & Onishi, T. Isotopic equilibration reaction of dinitrogen over Raney ruthenium: importance of the structural factor. J. Catal. 112, 469–477 (1988).
Grunze, M. The interaction of hydrazine with an Fe(111) surface. Surf. Sci. 81, 603–625 (1979).
Egawa, C., Naito, S. & Tamaru, K. Adsorption and decomposition of NO and NH3 on Ru(001) and Ru(1,1,10) surfaces. Surf. Sci. 138, 279–291 (1984).
Acknowledgements
This research was financed by the Swedish Research Council (Vetenskapsrådet, VR) under award numbers VR 2013-8823 and 2017-00559, as well as the Knut and Alice Wallenberg Foundation under award number KAW 2016.0042. We acknowledge Deutsches Elektronen-Synchrotron DESY (Hamburg, Germany), a member of the Helmholtz Association (HGF), for the provision of experimental facilities. Parts of this research were carried out at PETRA III using beamline P22. Beamtime was allocated for proposals I-20200291 EC, I-20200292 EC and II-20211048 EC. J.K.M. is grateful for financial support from VILLUM FONDEN (research grant 41388).
Funding
Open access funding provided by Stockholm University.
Author information
Authors and Affiliations
Contributions
A.N. and C.M.G. proposed the study. C.S. developed the beamline. C.M.G., A.N., D.D. and P.A. developed the experimental set-up. C.M.G., P.L., D.D., M.S., P.A., F.G.-M., S.K., B.D., J.K.M. and R.R. performed the experiments. C.M.G. analysed the data. G.L.S.R. performed the binding-energy calculations. A.N., C.M.G. and B.D. wrote the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature thanks Jorge Boscoboinik and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Comparison of N1s over ruthenium under pure N2 and low-pressure condition.
The spectra of N1s on \({\rm{Ru}}(10\bar{1}3)\) under 200 mbar 1:3 N2:H2 gas mixture at 423 K and the thickest nitride film made on the same surface at 623 K. The black lines show the offset between the NH species in red in the top spectra to N on terraces and light green in the bottom spectra cumulating to 150 meV.
Extended Data Fig. 2 N1s data at every temperature and pressure.
a–c, 200 mbar 1:3 N2:H2 nitrogen spectra over Fe(110), Fe(210) and \({\rm{Ru}}(10\bar{1}3)\), respectively. d–f, 500 mbar 1:3 N2:H2 nitrogen spectra over Fe(110), Fe(210) and \({\rm{Ru}}(10\bar{1}3)\), respectively. Temperature increases from top to bottom, from 423 K to 723 K.
Extended Data Fig. 3 Stability of long acquisition times.
a, The set of mass spectrometer data before any processing gathered over a 2-h window. The trend in masses 16, 17 and 18 are because of the slow improvement of the vacuum conditions. b, The same data as in a with processing as described in Methods. c, XPS spectra gathered simultaneously normalized to background signal. d, The sum of the spectra with the fitted peaks shown according to the colour scheme used throughout the paper.
Extended Data Fig. 4 Mass spectrometer time trace with and without ammonia production.
a, The equivalent experiment as for the data shown in Fig. 1b but at 523 K over Fe(110) with the complete set of mass fragments shown before any data processing, as well as the ammonia signal with all of the data processing. b, The equivalent experiment but at 673 K over Fe(110). Note here that the timescale is far longer than in Fig. 1b to ensure sufficient statistics, ensuring that the relative reactivity is representative of the lack of ammonia production. Note that, in a, there is a glitch in the mass spectrometer background as a response to the removal of N2, but it is clear that mass 15 is decreasing in pure H2 after the switchover. c, The equivalent experiment as b but performed without any X-rays on the sample.
Extended Data Fig. 5 Beam flux effect on surface species.
a, The trend of N1s intensity over Ru at 500 mbar, 423 K, 1:3 N2:H2 gas ratio as beam flux increased, with a line of best fit in grey (intercept set to zero) showing that the beam does not have an appreciable effect on total coverage. b, The same experiment showing that the N1s spectra gathered sequentially from lowest flux to full flux, showing no change in speciation with beam flux.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Goodwin, C.M., Lömker, P., Degerman, D. et al. Operando probing of the surface chemistry during the Haber–Bosch process. Nature 625, 282–286 (2024). https://doi.org/10.1038/s41586-023-06844-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-023-06844-5
This article is cited by
-
Reverse water gas-shift reaction product driven dynamic activation of molybdenum nitride catalyst surface
Nature Communications (2024)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.