Comparison of NOR on various transition metals. A suite of transition metals including Ti, V, Cu, Ag, Au, Ni, Pd and Pt were prepared as electrodes (See details in Methods), and the catalytic activity was evaluated in a sulfuric acid-diluted K2SO4 (0.5 M, pH=3) aqueous solution under 1 atm of NO in a two-compartment gastight H-cell. Each transition-metal electrode was tested using chronoamperometry at multiple potentials. The gas products were quantified via an on-line gas chromatography (GC), and liquid products were confirmed off-line by three techniques, UV-vis, ion chromatography (IC) and nuclear magnetic resonance (NMR) spectroscopy, for accuracy (see Methods, Supplementary Figures S1-4). Here we should emphasize that the removal of oxygen in electrochemical system is essential to ensure the stability of NO before getting electrochemically reduced on cathode. Otherwise, the residual oxygen will chemically oxidize NO to NO2, NO2- and finally to NO3-, leading to an inflated current density and complex reaction kinetics (Supplementary Figure S5). Beyond product distribution (represented by Faradaic efficiency, FE), we also take emphasis on quantifying product yield rates. To accurately determine the intrinsic activity, the topological surface area (TSA) of electrode was quantified by measuring the topological surface roughness factor (see Methods, Supplementary Figure S6 and Table S1). Figure 1a describes the FEs and yield rates of NH3 at -0.6 V vs. RHE on different transition metals. Notably, Ni exhibited the highest FE (~94%) and highest yield rate (9.48 µmol cm-2TSA h-1) towards NH3 simultaneously. Cu shows a lower FE of 80% with a yield rate of 6.61 µmol cm-2 TSA h-1. Instead, the VIII group metal Pd and Pt are more favorable to generate H2, having a remarkable decreased NH3 selectivity of 59% and 51%, respectively. As another competing reaction pathway, NH2OH was observed with a decent selectivity on some metals, such as Ag, Ti and V. The comparison of activity and selectivity at multiple potentials was summarized in Supplementary Figure S7.
Concurrently, we employed DFT calculations to understand electrochemical NOR. It is confirmed that the adsorption free energy of N (∆GN*) is a good descriptor for NH3 production. Figure 1b demonstrates a volcano-shaped plot between calculated ΔGN* and NH3 production (shown by both FE and yield rate) on various transition metals. The ΔGN* of Ni is close to the optimal value, neither too strong nor too weak. Accordingly, Ni exhibited simultaneously the highest FE and highest yield rate towards NH3. As the adjacent transition-metal to Ni, Cu exhibits a weaker N adsorption, thus a lower FE and production rate. The poor performance of metal Pd and Pt is due to their remarkable capability to generate H2 [24–25], which is a major competing reaction to produce NH3. Moreover, Ti, V, Ag and Au show much lower selectivities and production rates because of large deviations from the optimal value of ΔGN*.
To further understand the activities and selectivities of these transition metal electrodes, the electrochemical process of NOR were theoretically investigated on the most stable facets of various transition metals (see Methods, Supplementary Figure S8). After screening possible NOR pathways, two feasible ones, O* and N* pathways, were proposed as listed in Supplementary Table S2. Rather than directly dissociating into atomic O* and N*, the NO molecule forms adsorbed NHO* (O* pathway) or NOH* (N* pathway) intermediate at the hollow site of metal atoms through an associative pathway [16]. Regardless of each pathway, the adsorption energies of intermediates involved in the rate-determining steps (RDS), NO→NHO* for O* pathway and NH*→NH2*for N* pathway, depend almost linearly on ΔGN* due to the scaling relation (Supplementary Figure S9), which can well explain the remarkable correlation between ΔGN* and NH3 production.
Design and fabrication of single-crystal Ni foils. To understand the active site motifs of Ni towards NOR, we rationally designed monocrystalline Ni with well-defined lattices as electrocatalysts. According to the Miller indices {hkl}, Ni facets can be classified into low-index facets including {100}, {111} and {110}, and high-index facets which constitute terraces, steps and/or kinks with variable orientations. A general microfacet notation for cubic structure has been proposed by Somorjai in 1980 [26], in which each high-index facet can be described by a combination of low-index facets with different constitution ratios. For example, Ni(210) is denoted as Ni(S)-[11(110)+11(100)], indicating that one-atom-wide (100) terraces are separated by one-atom-wide (110) steps. Likewise, Ni(310), denoted as Ni(S)-[11(110)+22(100)], can be visualized as a kinked surface with a repetitive unit having two-atom-wide (100) terraces and one-atom-wide (110) steps. The presence of low-index terraces/steps suggests a high degree of site heterogeneity on high-index facets, which can be clearly observed in the side view of Ni(210), where the kinks are shown to have atomic sites with six-(pink), nine-(yellow) and eleven-(blue) fold coordination (Figure 2a). To elucidate the structure-function correlation, we rationally prepared five single-crystal Ni facets with distinct coordination environments and site motifs, including two low-index facets, e.g., Ni(100) and Ni(111), and three high-index facets, e.g., Ni(210), Ni(310) and Ni(520). Those high-index facets can be denoted as different combinations of low-index facets, and are located on the inverse pole figure as shown in Supplementary Figure S10.
The growth of single-crystal Ni follows our recently developed “seeded abnormal grain growth” technique [27]. Figure 2b shows an optical image of a large-format single-crystal Ni foil with (111) facet, the crystal orientation of which is corroborated by the low energy electron diffraction (LEED) (Figure 2c). The uniform color contrast of five single-crystal Ni foils shown in electron-back scattered diffraction (EBSD) at the micro-meter scale strongly suggests the absence of in-plane misorientation on surface (Figure 2d). X-ray diffraction was engaged to track the crystallinity of bulk Ni foils. As shown in Figure 2e, XRD symmetric scans confirm the single crystallinity nature of fcc Ni foils textured in five out-of-plane directions, which are in good agreement with EBSD mappings. On the whole, these characterizations proved that the synthesis of single-crystal Ni is atomically controllable.
Understanding NO-to-NH 3 conversion on Ni. The electrochemical performance of single-crystal Ni was evaluated under the identical condition as we described above. Selectivity and yield rate on five single-crystal Ni facets were plotted as a function of applied potential, as summarized in Figure 3a and Supplementary Figure S11. Remarkably, NH3 is the major product on monocrystalline Ni at a more positive bias. The hydrogen evolution reaction (HER), as a major competing reaction, was gradually dominant as the applied potential increased negatively from -0.8 V to -1.0 V vs. RHE, due to the mass transfer limitation of NO (solubility ~1.94 mM at 25 ℃) [28] under a large electrochemical polarization. In the medium potential region, NOR kinetically dominates the reaction. To exclude the influence from mass transport, we focused on the medium potential region to understand NOR. As shown in Figure 3a, all high-index Ni facets explicitly demonstrate higher selectivity towards NH3. In particular, Ni(210) achieved almost a unique selectivity with a yield rate nearly two-folds higher than those of Ni(100) and Ni(111) (Supplementary Figure S11). It is worth to note that the facet orientation of all selected single-crystal Ni foils remained stable after electrochemical NOR, confirmed by EBSD characterization (Supplementary Figure S12).
We performed DFT calculations to understand the activities on various Ni facets, including both N* pathway (NO→NOH*→N*+H2O, N*→NH*→NH2*→NH3) and O* pathway (NO→NHO*→NH2O*→NH3+O*, O*→OH*→H2O)[18]. DFT results suggest that N* pathway with NH*→NH2* as the RDS dominants NOR on Ni(111), Ni(100), Ni(310) and Ni(520) facets (see Supplementary Figures S13-14 and Table S3 for adsorption energies of all intermediates), and the energy increase associated with RDS are 0.16, 0.33, 0.10, and 0.09 eV, respectively (Figure 3b). Since the stepped surface Ni(310) and Ni(520) demand less energy than that of flat surface Ni(100) and Ni(111), they would exhibit better NH3 production performances, which is in-line with the experimental results in Figure 3a. For Ni(210), both N* and O* pathway were plotted in Figure 3c. Although the N* pathway exhibits a fairly smaller barrier of 0.07 eV, the O* pathway is energetically preferred on Ni(210) as it shows a continuous downhill conversion, which can well explain its best performance towards NH3 production in experiment.
We also observed a switch of adsorption sites for intermediates in the RDS of N* pathway on all involved facets. DFT calculations show that the first three generated species, NOH*, N* and NH*, prefer to occupy a hollow active site, the NH2* intermediate, instead, is only stabilized at a bridge site. Therefore, the adsorption site has to switch from a hollow to a bridge site at the step of NH*→NH2*. It is presumable that the free energy demanded for such switch, denoted as ∆Gs, may determine the activity towards NH3 production. Figure 3d plotted the NH3 yield rate in log scale as a function of ∆Gs on Ni(111), Ni(310), Ni(210) and Ni(520), which demonstrates an roughly exponential correlation. As for Ni(210), although it favors O* pathway, the hydrogenation of O* to produce OH* similarly associates with a location switch from hollow to bridge site, which is exothermic considering the energetic downhill nature of O* pathway on Ni(210). In this regard, ∆Gs can be regarded as 0 eV on Ni(210). Suprisingly, the ∆Gs of O* pathway on Ni(210) appears along the trend line in Figure 3d, suggesting ∆Gs is a reasonable theoretical parameter to describe the yield of NH3. Such correlation can be well explained from a view of geometric structure. Ni(100) displays the largest ∆Gs as the coordination number of NH* has to decrease from four (hollow site) to two (bridge site) associated with its hydrogenation. Instead, the NH* coordination only need to drop from three (hollow site) to two (bridge site) at the closely-packed Ni(111) surface, which exhibits a smaller ∆Gs in relative to Ni(100). For high-index Ni facets, the adsorption capability of bridge site will be strongly elevated because of the abundant undercoordinated step sites, which effectively reduce the required ∆Gs and ultimately accelerate the yield of NH3.