Uncovering the Mechanism of the Hydrogen Poisoning on Ru Nanoparticles via Density Functional Theory Calculations

: Although hydrogen plays a crucial role in ammonia synthesis, very little is known about its poisoning of Ru catalysts. In this study, density functional theory calculations of H 2 and N 2 dissociations, and H atom binding on Ru 153 were performed to provide a fundamental understanding of hydrogen poisoning. Because of the kinetic dominance of the H 2 dissociation over N 2 (vertically or horizontally adsorbed) splitting, the dissociated H atoms block the active sites required for horizontal (less energetically demanding dissociation) N 2 adsorption to occur either from the gas phase or after its geometrical transformation from being adsorbed vertically. Additionally, the dissociated H atoms withdraw electrons from the surface, which reduces the ability of the neighboring Ru atoms to donate electrons for N 2 activation, hindering its dissociation and suppressing ammonia synthesis.


Introduction
Liquescent NH 3 is both a carbon-free fuel and hydrogen carrier, spurring the development of highly efficient ammonia synthesis (AS) catalysts. Most NH 3 is synthesized by the Haber-Bosch process, an energy-intensive process, as the AS reaction over iron-based catalysts takes brute force; more than 20 MPa of pressure and temperatures higher than 450 • C are required to achieve a practically interesting conversion [1,2]. Because of its high activity under mild conditions, the Ru catalyst is more suitable for green AS than Fe-based [1,[3][4][5]. Thus, one of the strategies to decrease the energetic requirements of the Haber-Bosch process is replacing the Fe-based catalyst with Ru [1,3,4]. However, in the absence of basic promoters, the catalytic activity of Ru for AS is very little [6][7][8][9][10]. Generally, alkali or alkali earth promoters enhance the activity of the Ru catalyst via electron injection and, thus, lower the energy barrier for N 2 dissociation, the rate-determining step in AS [5,11]. Unfortunately, some of these promoters also stabilize H 2 adsorption; thus, the dissociated H and N atoms compete for the same sites leading to high surface coverage by H atoms, a process also known as hydrogen poisoning [11]. The negative H 2 reaction order for Ru nanoparticles suggests that the catalysts should operate at nonstoichiometric conditions, which are less thermodynamically favorable for the AS reaction [12,13]. Thus, a supported Ru catalyst that is less inhibited by hydrogen poisoning is highly desirable. Promoters such as Ba and La help against the inhibition effect of H 2 on the activity of Ru/MgO catalysts by a faster intrinsic rate of N 2 dissociation and increasing the competition for active sites with adsorbed H atoms [5].
Similarly, the synergy of C12A7e − to reversibly store/release hydrogen and donate electrons to the Ru particles suppressed hydrogen poisoning and accelerated the cleavage of the N-N triple bond [14]. Another practical strategy to increase the catalyst's lifetime against H poisoning is alloying Ru, a catalyst that interacts weakly/strongly with N/H, and Fe, which interacts weakly/strongly with H/N [15]. The FeRu solid solution exhibited a positive H 2 reaction order, attributed to a reduced overlap of Ru d-states with the s-states from H compared to pristine Ru(0001) [15]. Theoretically, the H atom binding [15][16][17] and its diffusion [18,19] on Ru surfaces and the H 2 dissociative adsorption on Ru nanoparticles [20][21][22] were studied extensively. Modulation-excitation infrared spectroscopy (MEIRS) measurements together with density functional theory (DFT) calculations revealed that the dissociation of activated H 2 molecules is more probable than their adsorption, leading to the hydrogenation of the surface of the Ru catalyst [23]. Given the critical role of hydrogen in the AS reaction, little is known about the H poisoning mechanism. Theoretically, to model the vast intricateness of the effects, conditions, and their synergy during the reaction conditions, such as the adsorbate's coverage, the effect of the support, impurities, defects, the surface roughness of the support and nanoparticles, the size, shape, symmetry, and dispersion of the nanoparticles, the pre-treatment effect, the effect of pressure and temperature, etc., is impractical. Thus, theoretical researchers must find structures or models that are representative of the experimental conditions. Homogeneous surfaces or small clusters fail to explain the adsorption properties of real nanoparticles because, unlike ideal flat surfaces, nanoparticles feature a high degree of site heterogeneity, which depends on nanoparticle size and shape [23,24]. The Ru(0001) is the most stable surface orientation of the hexagonal close-packed Ru and has been widely used to study the H binding [15,16,[24][25][26][27][28], N 2 adsorption [15,26,[28][29][30][31], and NH 3 formation paths [15,16,32]. However, in the case of the AS, the activation of N 2 occurred only on Ru steps, which confirmed that the reaction is extremely sensitive to the structure of the Ru catalyst [16,[25][26][27]. MEIRS measurements and DFT calculations exemplified the Ru(0001) limitation for the N 2 adsorption [23]. Under conditions similar to AS, the IR signal of N 2 adsorbed on Ru nanoparticles supported on MgO was deconvoluted into two major and several minor peaks [23]. DFT calculations of N 2 adsorption on Ru(0001) were only able to assign one of the largest peaks with considerable deviations from the experimental wavenumber [23]. Additionally, for the N 2 dissociation from its vertical interaction on flat [11] and stepped [13] Ru(0001), the horizontal N 2 adsorption appeared as a metastable molecularly adsorbed state. On the other hand, DFT calculations using a nanoparticle model containing 153 Ru atoms, Ru 153 , were able to assign the multiplicity of N-N stretching modes to different vertical and horizontal N 2 adsorption configurations [23]. While some deviations were observed, the DFT-calculated and experimental wavenumbers were in good qualitative and semi-quantitative agreement to explain the multiple N 2 /H 2 adsorption states and the source of the molecules' activation on Ru/MgO catalysts [23]. Moreover, DFT calculations of the H-H and H-Ru stretching modes from the H 2 and H atom interaction on Ru 153 , respectively, showed excellent agreement with the experimental vibration frequencies [23]. For the reasons mentioned above, herein, DFT calculations of the H 2 and N 2 dissociations and H atom interactions on Ru 153 were used to explain the mechanism of H poisoning on Ru catalysts. Our results show that H poisoning results from the barrierless H 2 dissociation. The dissociated H atoms physically and electronically hinder N 2 adsorption, activation, and dissociation, by which the AS is inhibited. These results are intended to illuminate the characteristics of H poisoning and help disrupt its effect on the activity of next-generation AS catalysts that are vital to establishing a viable carbon-free society.

N 2 Dissociation
The most stable configurations of N 2 adsorption on Ru 153 were reported to be vertically on top of a vertex atom with an adsorption energy of −1.031 eV, and at a 4-fold (4F)coordinated site with an adsorption energy of −1.021 eV [23]. From these two most stable configurations, the dissociation of the N 2 molecule has been studied, and all the stable structures are summarized in Figure S1 (top) and Figure S2 (4F). The N 2 dissociation energies from the top site range from −0.468 to −1.854 eV. The most stable configuration corresponds to each N atom binding to a 3-fold (3F)-coordinated site while sharing the same vertex atom. Both N atoms are charged negatively, −0.433 and −0.441 e, with average Ru-N bond orders of 0.970 and 0.938, respectively. Theoretically, electron donation from C12A7:O 2− , C12A7:O 2− /Ru 2 , C12A7:e − , and C12A7:e − /Ru 2 to N atoms has also been reported [33], and it was hypothesized that electronic effects such as depletion of surface charge gives rise to a new hydrogen adsorption site and weakens the H 2 chemisorption state that exists on the bare surface [34]. The electron density difference for the dissociated N 2 on Ru 153 is shown in Figure S3a. The energy barrier required to dissociate N 2 from its adsorption on top of a vertex atom to each N atom binding to a 3F site is approximately 1.96 eV, which is higher than our calculated energy barrier for the N 2 dissociation on Ru(0001), 1.49, and higher than the calculated energy required to split N 2 (0.66 eV) on Ru(1011) [35]. The transition state (TS) structure corresponding to the dissociation from the top site of Ru 153 is one of the N atoms on top of the vertex atom, while the remaining N atom is on a bridge/3F site with a distance of 2.057 Å between the N atoms. The N 2 dissociation pathways from its most stable adsorption on the top site of Ru 153 , vertically adsorbed on the top of a Ru atom of Ru(0001), and Ru(1011) [35] are represented in Figure 1a. The dissociation energies of N 2 from the 4F site range from −0.080 to −1.627 eV. After the N 2 adsorbs on a 4F site, the most stable dissociation leads to one of the N atoms binding to a 3F site and the other interacting on a bridge site. The N atoms for the most stable dissociated configuration are charged negatively, −0.404 and −0.392 e, with average Ru-N bond orders of 0.955 and 1.273 for the N atom binding on the 3F and bridge site, respectively. The electron density difference for the most stable N 2 dissociated from the 4F site of Ru 153 is shown in Figure S3b. The dissociation of the N 2 horizontally adsorbed on the 4F has an energy barrier of 0.94 eV, which is higher than our calculated activation energy (0.81 eV) for the N 2 horizontally adsorbed on a 3F site of Ru(0001), and higher than the energy required (0.79 eV) to split the N 2 adsorbed horizontally on a 4F site of Ru(1121) [27]. For the N 2 dissociation from the most stable 4F site of Ru 153 , the N-N distance at the TS is 1.902 Å, where each N atom is located at a bridge site. The difference in the energy barriers (1.02 eV) required to dissociate the N 2 from a top and 4F site arises from more electron donation from the four Ru atoms at the 4F site than a single atom at the top site [23]. Such charge transfer leads to a buildup of electron density at the N-N bond, destabilizing the N-N bond and activating the N 2 molecule [23]. These results prove the importance of the "side on" or horizontal adsorption of the N 2 to facilitate its dissociation (rate-determining step in AS) and speed up the reaction leading to the formation of NH 3 . The dissociation reaction from N 2 interacting on the most stable 4F site of Ru 153 is compared with the dissociation pathways of the N 2 horizontally adsorbed on Ru(0001), and Ru(1121) in Figure 1b.

H 2 Dissociation
The most stable configuration of the H 2 adsorption on Ru 153 is horizontally on top of the same vertex atom as for the vertical N 2 interaction, with an adsorption energy of −0.655 eV [23]. The dissociation energies of the H 2 from such vertex sites range from −0.516 to −1.348 eV, as pictured in Figure S4. The configuration corresponding to the most stable H 2 dissociation is with each H atom interacting on a bridge site between a vertex and edge atom while sharing the same vertex atom. Similar to the negatively charged N atoms, the dissociated H atoms exhibit charges of −0.131 and −0.124 e and averaged Ru-H bond orders of 0.540 and 0.560, respectively. The electron density difference for the most stable dissociated H 2 structure on Ru 153 is depicted in Figure S5. Note that the most stable configurations of the H 2 dissociation are less exothermic than the N 2 dissociation, both from the N 2 vertically and horizontally adsorbed. Notice that the dissociated H atoms ( Figure S4) interact with some of the same atoms as the dissociated N atoms shown in Figure S1. The energy barrier for the H 2 dissociation from the H 2 adsorbed on top of a vertex atom is 0.03 eV, which is the same for the H 2 dissociation on Ru(0001) and agrees with experimental observations of the spontaneous H 2 dissociation [36] and with MEIRS observations, suggesting that the dissociation of H 2 occurs immediately after the molecule's adsorption [23]. At the TS, the H atoms located on top of the vertex atom of Ru 153 are separated 1.214 Å. The dissociation reaction from H 2 adsorbed on top of the most stable vertex atom is compared with the molecule's dissociation on Ru(0001) in Figure 2.  [27,35], respectively. The configurations corresponding to the adsorbed, TS, and dissociated N2 on the Ru(1011) and Ru(1121) can be found in references [27,35], respectively. The TS and energy barriers for the N 2 dissociation are also shown.

H Atom Binding
From exploring all the possible binding sites of the H atom on Ru 153 , 25 symmetric adsorption sites were uncovered. The H atom interaction on the top of a vertex atom ( Figure S6) has a binding energy of −2.512 eV. Figure S7 summarizes the H atom binding on the bridge sites of Ru 153 with binding energies ranging from −2.499 to −2.973 eV, from which the interaction of the H atom at the bridge site between two edge atoms is the most stable. The binding energies of the H atom interaction on the Ru 153 s 3F sites range from −2.751 to −2.913 eV, as shown in Figure S8. From all the stable configurations summarized in Figure S8, the most stable is for the H atom interacting with two edge atoms and one {1011} facet atom. The binding energies of the H atom interaction on the 4F sites of Ru 153 range from −2.554 to −2.770 eV, from which the most stable one is the H atom binding with two edge and two atoms of the {1011} facet, as presented in Figure S9. Comparable to the H atoms from the H 2 dissociation, the Ru atoms donate electrons to the interacting H atom. It was reported that the H atom interacting on the threefold hollow site of Ru(0001) gains 0.051 e [37], which is in good agreement with our calculated value, 0.049 e ( Figure S8 Figure 3 using the projected crystal orbital Hamilton population (pCOHP) and the density of states (DOS) profiles.
Note that the DOS profiles of the Ru atoms at the adsorption site are averaged by the number of Ru atoms involved in the adsorption. Some dissimilarities exist between the DOS profiles of the individual Ru atoms of the most stable H binding on the 3F and 4F sites, and their averaged DOS profiles. These dissimilarities arise from the high degree of site heterogeneity arising from the unique features and properties (coordination number, charge, d-band center, etc.) of the symmetric atoms in Ru 153 . The heterogeneity of the Ru 153 surface sites is represented as heat maps of each Ru atom coordination number, charge, and d-band center elsewhere [23]. For the most stable structures of the H atom binding on the bridge, 3F, and 4F sites, the individual DOS profiles are compared to the averaged DOS in Figure S10. In this work, we use the averaged DOS profiles because it helps us explain the differences between the most stable structures of the H atom binding on top, bridge, 3F, and 4F sites. From the most stable H binding on the top, bridge, 3F, and 4F sites, the stability decreases in the following order: H atom interacting on the bridge, 3F, 4F, and top site. While some deviations arise from the difference in the size of the nanoparticle and its symmetry, the H atom binding energies reported here are in semi-quantitative and good qualitative agreement with DFT calculations of the H atom interaction on Ru 55 [22]. On the other hand, theoretical [16,28] and experimental [38] works reported that the H atom thermodynamically favors the threefold hollow sites of Ru(0001), with binding energies of −2.98 [16] and −2.96 eV [28]. From the pCOHP and the DOS analyses, the most stable interaction of the H atom binding on the bridge site is dominated by the coupling of the Ru atoms' 5s and e g orbitals with the 1s orbital from the H atom. Compared to the bridge interaction, the H atom binding on the top site is dominated by the combination of the Ru atom's 5p y , 4d xy , and the e g orbitals, which are more upshifted compared to all the other adsorption sites. Additionally, the antibonding states are the most downshifted among all the most stable adsorption sites considered, helping to explain the weak binding on the top site. The difference in stability between the H interaction on 3F and 4F sites arises from the lower coupling of the Ru atoms' 4d xy , 4d yz , and 4d x 2 −y 2 orbitals at the 4F site with the 1s from H atom. In Figure S11, the orbital combinations involved in the H-Ru interaction are decomposed for a more straightforward understanding.

Discussion
From the reported H 2 /N 2 adsorption states on Ru/MgO [23], there is a tendency that the N 2 adsorption is more stable than the H 2 when the adsorption site involves the same atoms (except for the vertical N 2 at the bridge site and N 2 horizontally on the top site). Although the difference in adsorption energies is sensitive to each of the adsorption sites in the nanoparticle, the more stable N 2 adsorption would be translated into N 2 molecules covering a large portion of the Ru nanoparticle. Additionally, DFT calculations [16] showed the more stable N atom interaction, −5.70 eV, compared to the binding energy of H atom, −2.98 eV on Ru(0001). For Ru(0001) pre-covered with N atoms, the presence of the N atoms inhibits H 2 adsorption [16], and it was hypothesized that as the concentration of N atoms on the surface increases, the depletion of the surface charge leads to a new hydrogen adsorption site [34]. Thus, H 2 must dissociate quickly before N 2 adsorption for H atom poisoning to occur. Such affirmation agrees with our calculated activation energy of 0.03 eV and MEIRS observations where H 2 and H atom peaks were observed on Ru/MgO [23]. It is even probable that the dissociated H atoms could hinder further H 2 adsorption and its dissociation as suggested by the time-domain IR spectrum where the peak corresponding to H 2 adsorption on Ru/MgO decreased with time and the one corresponding to H atom binding did not vary significantly, suggesting that the hydrogenation of the Ru nanoparticle proceeds with time [23]. When gas phase H 2 approaches the surface of the Ru 153 to interact vertically with the nanoparticle's adsorption sites, the molecule either adsorbs physically at distances between 3.313-3.652 Å or rotates to chemisorb horizontally as previously reported [23]. These results suggest that the H 2 orientation is important and that horizontal interaction is required for the activation of H 2 before its non-activated dissociation. Although the dissociated H atoms and the N 2 molecules compete for the same adsorption sites, some N 2 molecules could adsorb vertically on top of the Ru atoms as the dissociated H atoms prefer to bind to bridge, 3F, and 4F sites before the top sites, as shown in Figure S4 and from the interaction of an H atom on Ru 153 (Figures S6-S9). The effect of H 2 on the isotopic equilibration of N 2 on various Ru catalysts showed that the fraction of adsorption sites occupied by H atoms was ca. 5:1 compared to adsorbed N atoms [39]. A kinetic model revealed that the H atoms coverage is 46 times larger than the coverage of the N atoms [5]. While adsorbed N and H atoms retard N 2 activation over different Ru catalysts, the main retarding species are the adsorbed H atoms [40].
Similarly, there is a lower probability of horizontal N 2 adsorption either from the gas phase or after its rotation and diffusion from its vertical interaction because of the stable H atom binding preferences. To corroborate this statement, Figure S12 summarizes the possible H 2 dissociated configurations where the H atoms block one of the adsorption sites required for horizontal N 2 adsorption at the 4F site. Obstructing N 2 horizontal adsorption translates to a lower NH 3 production rate because of the higher energy barrier required to split the N 2 vertically adsorbed (1.96 eV) compared to the horizontal configuration (0.94 eV). Figures S1 and S2 summarize the stable configurations of the N 2 dissociation from the top and 4F sites, respectively, and Figures S4 and S12 show the stable configurations of dissociated H atoms that are occupying at least one of the same binding sites required for the N 2 dissociation. These results can explain why the activity of Ru catalysts in the absence of basic promoters is minimal [6][7][8][9][10], and the necessity of such promoters to improve the activity of the Ru catalyst via a trade-off effect between a faster N 2 dissociation and greater competition for the active sites with H atoms [5].
Besides the physical obstruction of the active sites for N 2 adsorption, activation, and dissociation, from the electron density difference and charge analyses ( Figure S5), the dissociated H atoms withdraw and localize electron density from the Ru nanoparticle. It has been shown that N 2 activation occurs from electron density accumulation between the N-N bonds due to electron donation from the Ru atoms to the N 2 molecule [23]. Thus, as the Ru atoms donate electrons to stabilize the dissociated H atoms, the surface Ru atoms will be less capable of donating electrons for the activation of N 2 , further dissociation, and stabilizing the dissociated N atoms that also withdraw and localize electron density. Following the hypothesis that the increased concentration of chemisorbed N atoms leads to a depletion of the surface charge [16], then H atoms present on the catalyst's surface would decrease the electron donation ability of the surface Ru atoms since both N and H atoms withdraw electrons. However, more calculations are required to understand the co-adsorption of N 2 and dissociated H atoms and their hindering effect on N 2 adsorption and activation. Additionally, the presence of preabsorbed atoms or molecules that can exhibit repulsive forces with the incoming H 2 molecules to prevent the spontaneous H 2 dissociation should also be addressed to avoid H atom poisoning.

Computational Methods
All calculations were based on the plane-wave DFT method implemented in the Vienna ab initio simulation package (VASP 5.6.2) [41][42][43]. Perdew−Burke−Ernzerhof parametrization under the generalized gradient approximation was employed as the exchange-correlation functional together with the projector-augmented wave method [44]. Spin-polarized calculations were performed throughout the study with a plane-wave cutoff energy of 600 eV. The convergence criterion for the geometry optimization calculations was set as the point at which the difference in the total energy between the two ionic steps was less than 10 −5 eV/atom. Such energy optimization criterion led to atomic forces of up to 0.003 eV/Å for the optimization of the nanoparticle. For the self-consistent-field iterations, a criterion of 10 −6 eV/atom was considered. The optimization of the HCP structure of the 1.58 nm diameter Ru nanoparticle, Ru 153 , was performed at the Γ point in reciprocal space owing to the significant spatial extent of the system, where all Ru atoms were allowed to relax. The distance between neighboring image Ru nanoparticles was set as 15 Å to avoid interaction between periodic images, as shown in Figure 4a. The nanoparticle model comprises two facets, {0001} and {1011}, 18 atoms located at vertices, and 42 atoms located at ridges, as illustrated in Figure 4b. We have previously investigated H 2 and N 2 adsorption on the 12 symmetric adsorption sites and their possible combinations on Ru 153 [23]. In this work, the dissociation energies of the most stable adsorbed N 2 and H 2 configurations are estimated. Additionally, the H atom binding to Ru 153 is described to complement the results of the H atom interaction only on the 12 symmetric sites of Ru 153 [23]. The N 2 and H 2 dissociation reactions on Ru(0001) were calculated on a 3 × 3 supercell containing seven atomic Ru layers. As shown in Figure 4a, the vacuum between the slabs was set to span a range of 15 Å to ensure no significant interaction between the slabs. The optimization was performed with a 4 × 4 × 1 Monkhorst−Pack k-point mesh for the Brillouin zone integration where the atoms of the bottom three layers were kept fixed, and those of the top four layers were allowed to relax. The convergence criterion for the geometry optimization calculations was set as the point at which the atomic forces were lower than 0.005 eV/Å. For the self-consistent-field iterations, a criterion of 10 −5 eV/atom was considered. The energies of N 2 and H 2 adsorption on Ru(0001) (E ads ), the energies for the dissociation of the molecules on Ru 153 and Ru(0001) (E dis ), and the H binding on Ru nanoparticles (E bin ) were calculated as: where E Ru n/s +adsorbate , E Ru n/s , and E adsorbate denote the total energies of N 2 , H 2 , or H atoms interacting with the Ru nanoparticle or N 2 , H 2 interacting with Ru(0001), the total energy of an isolated Ru nanoparticle or pristine Ru(0001), and the total energies of isolated N 2 , H 2 , and an H atom, respectively. According to this definition, negative adsorption/dissociation/binding energies denote more stable adsorbate-Ru interactions. While all atoms were relaxed during the optimization of the H 2 , N 2 dissociations, and the H atom binding to Ru 153 , only the adsorbate and the Ru atoms in the top four layers were allowed to relax for H 2 and N 2 adsorption/dissociation on Ru(0001). After the geometrical optimization of the adsorbate interaction on Ru 153 , the atomic forces were no greater than 0.020, 0.009, and 0.009 eV/Å for the N 2 and H 2 dissociation and the H atom interaction on Ru 153 , respectively. Because of the significant size of the system, no dipole correction was considered for the adsorbate-induced dipole. Bond order and charge analyses were performed using the sixth-generation density-derived electrostatic and chemical method [45].

Conclusions
To elucidate the mechanism of hydrogen poisoning of Ru catalysts, the N 2 /H 2 dissociations and the H atom interaction were investigated on a Ru nanoparticle model and compared to the N 2 /H 2 dissociation pathways on Ru(0001) using DFT calculations. Additionally, the energetics of the N 2 dissociation were compared to two stepped surfaces from the literature, Ru(1011) [35] and Ru(1121) [27]. Except for the horizontal N 2 adsorption on Ru(0001), the vertical and horizontal N 2 adsorption and dissociation are more thermodynamically stable on Ru 153 , followed by the planar Ru(0001), and the stepped surfaces. The opposite trend is observed for the energy barrier required to dissociate the N 2 ; stepped Ru surfaces, followed by the Ru(0001), and then Ru 153 . H 2 adsorbs and dissociates more stably on Ru 153 than on Ru(0001). However, for both Ru 153 and Ru(0001), the energy barrier to dissociate H 2 is almost zero. Because of the barrierless H 2 dissociation, the poisoning of the Ru catalyst can occur via (i) physical blocking of the surface sites, where the dissociated H atoms block the surface Ru atoms needed for N 2 vertical and horizontal adsorption. Compared to vertical N 2 adsorption, the dissociation of horizontal N 2 is less energetically demanding by 1.02 eV. However, the presence of dissociated H atoms will hinder horizontal N 2 adsorption from the gas phase or, after its geometrical transformation, from being adsorbed vertically. Additionally, from the dissociated N 2 and H 2 configurations, both N and H atoms prefer the same sites; hence, N 2 dissociation is stalled. (ii) Reducing the ability of neighboring Ru atoms to donate electrons to activate the adsorbed N 2 molecules in the vicinity of the dissociated H atoms. In this way, the catalytic activity of Ru toward N 2 dissociation and, henceforth, the AS, is suppressed.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/catal12030331/s1, Figure S1: Dissociated configurations from the most stable vertical N2 adsorption on a top site, Figure S2: Dissociated configurations from the most stable horizontal N2 adsorption on a 4F site, Figure S3: Electron density difference for the dissociated N2, Figure S4: Dissociated configurations from the most stable horizontal H2 adsorption on top of a vertex site, Figure S5: Electron density difference for the dissociated H2, Figure S6: H atom binding on the top site of Ru153, Figure S7: H atom binding on the bridge sites of Ru153, Figure S8: H atom binding on the 3F sites of Ru153, Figure S9: H atom binding on the 4F sites of Ru153, Figure S10: Comparison between the DOS profiles of the individual Ru atoms and their average DOS for the most stable structures of the H atom binding on bridge, 3F, and 4F sites, Figure S11: Atomic orbital decomposition of the H atom binding on the most stable top, bridge, 3F, and 4F sites of Ru153, Figure S12: Dissociated H atoms on the same sites required for the N2 adsorption and dissociation from the 4F site, and Figure  S13: Interatomic Ru-Ru distances of surface atoms of Ru153 and HCP bulk Ru.