Mechanistic Investigation of the Catalytic Decomposition of Ammonia

Feb 18, 2014 - Catalytic decomposition of ammonia (NH3) is a promising chemical reaction in energy and environmental applications. Density functional ...
0 downloads 0 Views 4MB Size
Article pubs.acs.org/JPCC

Mechanistic Investigation of the Catalytic Decomposition of Ammonia (NH3) on an Fe(100) Surface: A DFT Study Sang Chul Yeo,† Sang Soo Han,‡ and Hyuck Mo Lee*,† †

Department of Materials Science and Engineering, KAIST, 291 Daehak-ro, Yuseong-gu, Daejeon 305-701, Republic of Korea Center for Computational Science, Korea Institute of Science and Technology (KIST), Hwarangno 14-gil 5, Seongbuk-Gu, Seoul 136-791, Republic of Korea



S Supporting Information *

ABSTRACT: Catalytic decomposition of ammonia (NH3) is a promising chemical reaction in energy and environmental applications. Density functional theory (DFT) calculations were performed to clarify the detailed catalytic mechanism of NH3 decomposition on an Fe(100) surface. Specifically, the elementary steps of the mechanism were calculated for the general dehydrogenation pathway of NH3. The adsorption of two types of ammonia dimers (2NH3), locally adsorbed NH3 and hydrogen-bonded NH3, were then compared, revealing that locally adsorbed NH3 is more stable than hydrogen-bonded NH3. By contrast, the dehydrogenation of dimeric NH3 results in a high energy barrier. Moreover, the catalytic characteristics of NH3 decomposition on a nitrogen (N)-covered Fe surface must be considered because the recombination of nitrogen (N2) and desorption have an extremely high energy barrier. Our results indicate that the catalytic characteristics of the NH3 decomposition reaction are altered by N coverage of the Fe surface. This study primarily focused on energetic and electronic analysis. Finally, we conclude that Fe is an alternative catalyst for the decomposition of NH3 in COx-free hydrogen production. NH3 decomposition follows a dehydrogenation mechanism.8 The rate-limiting step for NH3 decomposition on an Fe catalyst remains controversial. Lanzani and Laasonen suggested that the rate-limiting step is the first dehydrogenation step (NH3 → NH2 + H), which occurs in a nanosized Fe cluster.23 Lin et al.39 investigated the adsorption and dissociation of NH3 on Fe(111) and found that the third dehydrogenation reaction (NH → N + H) is the rate-limiting step. Conversely, the major rate-limiting step in catalytic NH3 decomposition is reported as the associative desorption of N2 on the surface.28 While these findings have provided insight into the decomposition of NH3 on an Fe surface, the intermediate species and a detailed catalytic mechanism have not yet been elucidated. A comprehensive DFT study would facilitate the determination of the complete reaction pathway of NH3 decomposition on Fe-based catalysts. Maier et al.40 reported a combined experimental (scanning tunneling microscopy)−theoretical (DFT) approach, which was used to demonstrate that one NH3 molecule (called αammonia) is chemisorbed to the surface via formation of an N−Ru chemical bond and another NH3 molecule (β-ammonia) is hydrogen bonded to the α-ammonia in the adsorption of an ammonia (NH3) dimer on Ru(0001). Similar dimer config-

1. INTRODUCTION Ammonia (NH3) is a toxic compound that can cause environmental pollution. However, NH3 can also provide COx-free hydrogen (H2) production through the decomposition of the NH3 molecule, which is a promising technology in a hydrogen economy.1−6 Thus, the catalytic performance and surface chemistry of NH3 decomposition have received considerable attention. A fundamental understanding of catalysts for the decomposition of NH3 has become very important.7−15 In general, Ru and Ir are the most active catalysts for the decomposition of NH3.16−22 However, they are expensive, and their supplies are limited. Less expensive alternatives (such as Ni-, Co-, and Fe-based catalysts) with good catalytic activity are needed.23−28 Consequently, NH3 decomposition on a variety of metal surfaces has been extensively studied. Several fundamental studies concerning NH3 decomposition have been performed using materials such as Ir29−32 and Ru.21,33−37 Recently, Duan et al.28 studied NH3 decomposition on inexpensive alternative catalysts, such as Ni(111), Co(111), and Fe(110) surfaces, using density functional theory (DFT) calculations. These authors compared the catalytic reactivity, including the adsorption of NH3 and the energy barrier for the N2 recombination reaction, and also experimentally determined that Fe nanoparticles on carbon nanofibers are highly active and stable due to the supporting effect.38 © 2014 American Chemical Society

Received: November 6, 2013 Revised: February 15, 2014 Published: February 18, 2014 5309

dx.doi.org/10.1021/jp410947d | J. Phys. Chem. C 2014, 118, 5309−5316

The Journal of Physical Chemistry C

Article

urations were also determined for Ni(111)41 and Pt(111).42 These findings suggest that the adsorption of dimeric NH3 on an Fe surface will be similar, although this reaction has not been studied. In addition, according to Li et al.,43 N atoms that are strongly adsorbed on the Fe surface inhibit the NH 3 decomposition reaction. These authors used a plasma-driven catalyst to increase the NH3 conversion from 7.8% to 99.9%. These findings indicate the importance of considering the effects of an N-covered Fe surface on NH3 decomposition reactions. In this study, we investigated the catalytic performance of an Fe(100) surface in the decomposition reaction of monomeric NH3. To clarify the exact mechanism, we calculated the elementary steps, including the adsorption energy (Ead), energy barrier (Eb), and reaction energy (ΔE). The results indicated that the catalytic properties of monomeric NH3 dehydrogenation on Fe are similar to those on the most active catalysts, such as Ru and Ir.29,30,32,44 Subsequently, we investigated the adsorption of two types of dimeric ammonia (2NH3) on the Fe surface: locally adsorbed NH3 and hydrogen-bonded NH3. Locally adsorbed NH3 was more stable than hydrogen-bonded NH3. Thus, we confirmed that the catalytic activity in the decomposition reaction of 2NH3 is altered due to locally adsorbed NH3. We also performed a quantitative analysis of the effect of surface coverage using energetic and electronic analyses, thereby reconsidering NH3 decomposition as a function of the surface coverage of N on the Fe surface. This study contributes to the fundamental understanding of the structural, energetic, and catalytic properties of Fe catalysts for NH3 decomposition and explains why Fe-based catalysts present major disadvantages for application in hydrogen (H2) production.

To determine the transition states and energy barriers for NH3 decomposition reactions on the Fe surface, we employed the climbing image nudged elastic band method (CI-NEB) and accordingly constructed the minimum energy path.56 Moreover, we calculated the density of state to investigate electronic structures for each geometry.57 To estimate the catalytic activity, the adsorption energy (Ead), energy barrier (Eb), and reaction energy (ΔE) were calculated for each reaction. In addition, to analyze the electronic structure, we evaluated the energy of the d-band center, εd, for the 4d electronic state. These equations are provided in the Supporting Information.

3. RESULTS AND DISCUSSION 3.1. Determination of Adsorption Intermediates. To investigate NH3 decomposition reactions, we considered an Fe(100) surface, which is the stable Fe surface.28,39,58 For the Fe(100) surface, we employed a (3 × 3) supercell containing a five-layer slab in which the bottom two layers were fixed and all the other atoms are allowed to relax freely. The vacuum layer was up to 15 Å, as shown in Figure 1. Hollow (h), bridge (b), and top (t) sites were considered as adsorption sites for NH3, NH2, NH, N, and H, as summarized in Table 1.

2. COMPUTATIONAL METHODS We performed DFT calculations within a plane wave basis set using the Vienna Ab-initio Software Package (VASP) to investigate the atomic and electronic structures of the Fe surfaces, in which the project augmented wave method45−48 with a frozen-core approximation was used. The spin-polarized DFT calculation was performed using the generalized gradient approximation of the Perdew−Burke−Ernzerhof (PBE) exchange-correlation functional.49 We studied adsorption and dissociation of NH3 on Fe(100) surface. In a previous study,50 the Norskov group reported that the PBE functional well describes the experimental data of N adsorption on Fe(100) surface in DFT calculation. And the type of PBE functional is believed to give more accurate energies than the type of PW91, especially in the case of adsorption of atoms and molecules, energetically.51 As well, the Cater group reported that the physical properties of Fe system computed with PBE functional is in good agreement with the experiment values.52−54 Therefore, the PBE functional is able to correctly describe our system. Moreover, an energy cutoff of 400 eV was used, and a (5 × 5 × 1) Monkhorst−Pack grid was considered for k-point sampling. The geometries were optimized using a conjugate gradient method until the maximum force in every degree of freedom was less than 0.03 eV/Å. For the relaxation of the electronic degrees of freedom, we achieved a total energy convergence of up to 10−6 eV. This DFT calculation provides a lattice parameter of 2.83 Å, a bulk modulus of 1.76 Mbar, a cohesive energy of 5.16 eV, and a magnetic moment of 2.21 μB for the bcc Fe; this value of the magnetic moment is close to reported values.55

Figure 1. Structural configuration of a 3 × 3 supercell containing five layers of Fe(100). The h in a black circle denotes a hollow site, b denotes a bridge site, and t denotes a top site. These sites are the preferred adsorption sites for (a) NH3, (b) NH2, (c) NH, (d) N, and (e) H. The purple, blue, and white spheres represent Fe, N, and H atoms, respectively.

An NH3 species is preferentially adsorbed at the top sites with a binding energy of −0.92 eV through the formation of an Fe−N chemical bond with C3v symmetry because a lone pair electron from the N atom interacts with the d-orbital of an Fe atom. Here, the N−H and Fe−N bond distances are 1.03 and 2.17 Å, respectively. Upon adsorption, the H−N−H bonding angle is 107.7°, similar to the value obtained in a previous theoretical study of NH3 adsorption on Fe(110).28 Grunze et al.7 experimentally reported that the adsorption energies of NH3 on Fe(100) are in the range of −0.43 to −0.74 eV, in reasonable agreement with our calculated values. An NH2 5310

dx.doi.org/10.1021/jp410947d | J. Phys. Chem. C 2014, 118, 5309−5316

The Journal of Physical Chemistry C

Article

Table 1. Adsorption Properties of NHx (with x = 0−3) Intermediates on Fe(100)a species NH3 NH2 NH N H NH3 (molecule)

adsorption site

Ead (eV)

dN−H (Ǻ )

dN−Fe (Ǻ )

∠HNH (deg)

top bridge hollow hollow bridge

−0.92 −3.01 −5.22 −7.35 −0.64

1.026 1.024 1.037

2.173 1.515 0.857 0.432

107.7 107.8

1.021

106.5

a

The energy was calculated with respect to the corresponding species in the gas phase and with respect to the NHx (with x = 0−3) species and the H atom.

species prefers to adsorb at the bridge site with an adsorption energy of −3.01 eV and N−H and Fe−N bond distances of 1.02 and 1.52 Å, respectively. Furthermore, an NH species is adsorbed at the hollow site with a binding energy of −5.22 eV and N−H and Fe−N bond distances of 1.04 and 0.86 Å, respectively. The preferential adsorption site for an N atom is a hollow site with a binding energy of −7.35 eV, while that for an H atom is a bridge site with a binding energy of −0.64 eV.55 Previously,59 we reported that an N2 molecule prefers a hollow site with an adsorption energy of −1.03 eV, with an N−N distance that is 1.243 Å longer than that of N2 in the gas phase. The adsorption energies of NHx (with x = 1−3) on Fe(100) exhibit a trend of NH > NH2 > NH3. This trend can be analyzed using projected density of states (PDOS) and d-band center theory, which indicates the trend in the reactivity of metals or the effects of alloying, surface configuration, and adsorbate−adsorbate interactions.60 The d-band center is calculated as shown in the Supporting Information. The DOS for the d-band center of the Fe atoms that are in contact with the NH3 species and the p-band of the nitrogen atom in the NHx (with x = 0−3) intermediate states were determined for the NH3, NH2, and NH species. The Fermi level was set to 0 (Figure S1). The black line denotes a clean Fe(100) surface prior to the adsorption of each species. By contrast, the red line indicates the adsorption conditions on the surface. We focused on the d-band center (εd) of each species (NH3, NH2, and NH) on the Fe(100) surface and determined εd = −0.87, −0.57, and −0.02, respectively. The d-band center of each species on the Fe(100) is close to the Fermi level, in the order NH3 < NH2 < NH. These features are similar to those describing the adsorption of NHx (with x = 0−3) species on Fe(100), as shown in Figure S1. 3.2. Monomeric Ammonia (NH3) Decomposition on Fe(100). To obtain a reaction pathway diagram for NH3 decomposition on Fe(100), we first investigated the relative elementary reaction steps or the NH3 dehydrogenation, as shown in Figure 2. The sequence of intermediate steps for NH3 decomposition can be written as follows: NH3* + * → NH 2* + H*

(1)

NH 2* + * → NH* + H*

(2)

NH* + * → N* + H*

(3)

Figure 2. Reaction pathway diagram of ammonia (NH3) decomposition on Fe(100). TS indicates the transition state.

reaction of NH3, which is important to estimate the catalytic activity. The initial and final configurations were determined at the lowest energy barrier from the calculations for the different reaction pathways. The upper three reaction steps for the NH3 decomposition reaction can also be considered geometrical changes, and they were defined as the initial, transition, and final states (IS, TS, and FS, respectively) of the elementary steps. Reaction 1 begins with the stable configuration at the top site, followed by dehydrogenation into NH2 at the bridge site and the presence of H at the bridge site. This step has an energy barrier of 0.95 eV and is exothermic by −0.26 eV. For the second reaction 2, NH2 begins with a stable configuration at the bridge site, followed by dehydrogenation into NH at a hollow site and H at the bridge site. This step has an energy barrier of 1.14 eV and is exothermic by −0.29 eV. For the final reaction 3, NH begins with a stable configuration at a hollow site, followed by dehydrogenation into N at a hollow site and H at the bridge site. This step has an energy barrier of 0.78 eV and is exothermic by −0.46 eV. On the basis of these results, we concluded that the rate-liming step for the decomposition of NH3 over an Fe(100) surface is the decomposition of NH2 to NH and H (reaction 2) because this step shows the largest energy barrier. To provide detailed information on the origin of the reaction energy barrier of the elementary reaction steps, we resolved the reaction energy barrier (Eb) by dividing it into the following six terms: def IS TS E b = ΔEsub + ΔEAB + EAB + E int − EATS − E BTS

(4)

where ΔEsub represents the energy of the substrate, which changes from its IS to the TS on the barrier (ΔEsub = ΔETS sub − def ΔEIS ), and ΔE is the deformation energy of the reactant sub AB (AB), which reflects the effect of the structural deformation of AB from the IS to the TS. A and B represent the adsorbates such as NHx (with x = 0−3) species and the H atom, respectively. IS and TS indicate initial state and transition state. This effect is also reflected in the way in which the N−H bonds are separated. EIS AB is called the binding energy of AB at the IS and is equal to the absolute adsorption energy of AB at the IS. ETS int represents the interaction energy between A and B at the TS. These four terms positively contribute to the energy TS barrier. Two other terms, ETS A and EB , negatively affect the

Here, the asterisk (*) denotes an empty site, and an NHx (with x = 0−3) intermediate species is adsorbed onto the Fe surface. We focused on the rate-limiting step for the decomposition 5311

dx.doi.org/10.1021/jp410947d | J. Phys. Chem. C 2014, 118, 5309−5316

The Journal of Physical Chemistry C

Article

Table 2. Energy Barriers and Contribution Factors (eV) for the Elementary Reactionsa reactions

ΔEsub

ΔEdef AB

EIS AB

ETS int

ETS A

ETS B

Eb

NH3 → NH2 + H NH2 → NH + H NH → N + H NH3 + NH3 → NH3 + NH2 + H NH3 + NH2 + H → NH3 + NH + 2H NH3 + NH + 2H → NH3 + N + 3H

0.03 0.05 0.19 0.09 0.12 0.26

3.56 2.21 1.28 3.65 2.73 1.55

0.92 3.01 5.22 0.48 2.36 3.34

1.26 1.98 2.57 1.61 2.23 2.92

2.38 3.69 5.86 2.41 3.68 4.63

2.44 2.42 2.62 2.45 2.43 2.65

0.95 1.14 0.78 0.97 1.33 0.79

a The energy was calculated with respect to the elementary reaction steps. A and B represent the adsorbates of the NHx (with x = 0−3) species and the H atom, respectively. IS and TS indicate initial state and transition state.

TS energy barrier. ETS A (EB ) represents the binding energy of the A (B) intermediate species in the TS without considering B (A). These values are shown in Table 2. To analyze the direct NH3 IS TS TS decomposition reaction, the four terms ΔEdef AB, EAB, Eint , and EA def should be studied in detail. The value of ΔEAB decreases, TS TS whereas those of EIS AB, Eint , and EA increase as the direct dehydrogenation process of NH3 proceeds. In particular, the relatively large interaction energy (ETS int ) at the TS leads to a higher energy barrier. The relatively small binding energies of TS EIS AB and EA lead to lower energy barriers. In addition, ΔEsub and ETS do not have any less effect on the energy barrier (Eb), B which is not as easily altered as the other terms. 3.3. Dimeric Ammonia (2NH3) Decomposition on Fe(100). For the dimeric ammonia (2NH3) decomposition, we considered the adsorption of two types of NH3 dimers. Preadsorbed NH3 leads to an additional locally adsorbed NH3 at the top site. Alternatively, an additional hydrogen-bonded NH3 is possible, as shown in Figure 3. To clarify the

NH3* + NH 2* + H* → NH3* + NH* + 2H*

(6)

NH3* + NH* + 2H* → NH3* + N* + 3H*

(7)

where the asterisk (*) denotes an empty site or an NHx (with x = 0−3) intermediate species adsorbed onto the Fe surface. In reaction 5, the locally adsorbed NH3 dimer has a binding energy of −0.48 eV as seen in Table 2, and the reaction begins with locally adsorbed NH3 on the top site of Fe. The locally adsorbed NH3 begins with the stable configuration at the top site, followed by dehydrogenation into NH2 at the bridge site and the presence of H at the bridge site. The energy barrier for the first decomposition reaction is 0.97 eV (see Table 2 or Figure 4), which is slightly higher than 0.95 eV for reaction 1 in

Figure 4. Reaction pathway of dimeric ammonia (2NH3) decomposition on Fe(100).

Figure 3. Configuration and energy difference between hydrogenbonded and locally adsorbed ammonia for dimeric ammonia (2NH3).

the mechanism of monomeric NH3 dehydrogenation. Moreover, the reaction is exothermic by −0.44 eV. For reaction 6, the first dehydrogenated NH2 begins with a stable configuration at the bridge site, followed by dehydrogenation into NH at a hollow site and H at the bridge site. This step has an energy barrier of 1.33 eV and is exothermic by −0.32 eV. In the final reaction step, NH begins with a stable configuration at a hollow site, followed by dehydrogenation into N at a hollow site and H at the bridge site. The energy barrier is 1.15 eV, and the reaction is exothermic by −0.43 eV. Overall, the reaction energy barrier of each step in the decomposition of dimeric NH3 is higher than for the decomposition of monomeric NH3, as shown in Figure 4. The rate-limiting step of dimeric NH3 dehydrogenation is therefore the second dehydrogenation reaction, which is identical to the dehydrogenation reaction of monomeric NH3.

decomposition of dimeric NH3, we determined the stability of the configuration between hydrogen-bonded and locally adsorbed ammonia. The energy state of locally adsorbed ammonia is lower than that of hydrogen-bonded ammonia as shown in Figure 3. Therefore, the decomposition of dimeric NH3 likely occurs through the dehydrogenation of locally adsorbed ammonia. To obtain a reaction pathway diagram for dimeric NH3 decomposition on Fe(100), we first investigated the relative elementary reaction steps. The sequence of elementary steps for the decomposition can be written as follows: NH3* + NH3* → NH3* + NH 2* + H*

(5) 5312

dx.doi.org/10.1021/jp410947d | J. Phys. Chem. C 2014, 118, 5309−5316

The Journal of Physical Chemistry C

Article

To investigate the energy barrier of the elementary steps in detail, we also considered the contribution of the energy barrier, as shown in Table 2. To analyze the locally adsorbed NH3 decomposition reaction, we analyzed the following terms IS TS def of ΔEdef AB, EAB, and EA . We observed that the value of ΔEAB increased for dimeric NH3 compared with monomeric NH3, whereas the values of EIS AB decreased as the NH3-hindered decomposition reaction proceeded. Moreover, the relatively strong binding energies of ETS int lead to higher energy barriers. By contrast, the higher energies of ETS A contribute to lower energy barriers. In addition, ΔEsub and ETS B do not have any less effect on the energy barrier (Eb), which is not as easily altered as the other terms. Consequently, the preadsorbed NH3 is constrained for the dehydrogenation of locally adsorbed NH3 at the top site. 3.4. Dependence of Ammonia Decomposition on NCovered Fe(100). The catalytic characteristics of metal surfaces are highly dependent on the adsorbates used to cover the surface. An energy barrier for the N2 recombination reaction, which is the highest energy barrier, of 2.17 eV has been reported by us.59 Moreover, the reaction is highly endothermic (1.03 eV), which is consistent with previous theoretical and experimental results.24,61,62 Therefore, nitrogen exists on the Fe surface and subsequently forms an N-covered surface. Thus, a quantitative analysis of the surface coverage effect via energetic and electronic analyses is required. We examined the adsorption energy and dissociation reaction to determine the effects of surface coverage on Fe(100). Initially, we optimized various N-precovered surfaces (θ = 0, 0.22, 0.56, 0.78, and 1.00 ML when θ is coverage and ML stands for monolayer) and analyzed the first dehydrogenation (NH3 → NH2 + H) elementary step of NH3. Quantitatively, we classified the precovered surfaces according to the number of Fe−N bonding events. As the N coverage increased, the number of Fe−N bonding events increased. Hammer investigated the N2 adsorption and dissociation at a Ru(0001) surface coated with N, O, or H using DFT63 and measured the N2 dissociation energy barrier and reaction energy for different levels of surface coverage of the adsorbate (N, O, or H). As discussed by this author, these energy changes are caused by the effects of N precoverage on the surface electronic structure, which affects the bonding strength in a manner suggested by previous theoretical and experimental results.64,65 To confirm the modification of the electronic structure due to surface precoverage, we calculated the d-band for the DOS of the Fe atom as the nitrogen surface coverage shifted. Moreover, the calculated d-band centers (εd) for an Fe atom with different numbers of Fe−N bonds, as shown in Figure 5, were as follows: −0.87 > −0.92 > −1.15 > −1.28 > −1.37 for θ = 0, 0.22, 0.56, 0.78, and 1.00 ML, respectively. Next, the adsorption energies (Ead) of NH3 on various N-precovered Fe(100) surfaces were calculated and are summarized in Table 3. The adsorption site of the NH3 molecule is located at the Fe top site, which establishes several Fe−N bonds, as shown in Figure 6. As the N precoverage increases, the NH3 adsorption energy decreases to −0.92, −0.72, −0.45, −0.41, and −0.39 eV. The first dehydrogenation reaction for the decomposition of NH3 at various levels of nitrogen surface coverage begins with a previously adsorbed NH3 at the Fe top site. The configurations of the initial, transition, and final states are shown in Figure 6. To determine the change in the catalytic reactivity as the N precoverage changes on the surface, the energy barrier (Eb) and the reaction energy (ΔE) were calculated using the CI-NEB method. The dissociation energy barrier (Eb) increased as the

Figure 5. Fe d-projected density of states for an Fe atom with different numbers of Fe−N bonds on Fe(100). The dotted line represents a Fermi level of 0.

Table 3. Coverage Dependence of the Adsorption and Dissociation Properties of Ammonia (NH3) on Fe(100) N coverage (ML)

no. of bonding Fe−Na

Ead (eV)

dN−H (Ǻ )

dN−Fe (Ǻ )

Eb (eV)

ΔE (eV)

0 0.22 0.56 0.78 1.00

0 1 2 3 4

−0.92 −0.72 −0.45 −0.41 −0.36

1.026 1.024 1.023 1.023 1.022

2.173 2.165 2.158 2.203 2.224

0.94 1.12 1.28 1.76 2.15

−0.44 0.46 0.78 1.31 1.94

a

Denotes the number of nitrogen atoms near the Fe atom.

Figure 6. Reaction pathway and configuration of the initial, transition, and final states for an ammonia-initiated dehydrogenation reaction on various nitrogen-precovered surfaces (θ = 0, 0.22, 0.56, 0.78, and 1.00 ML where θ is coverage and ML stands for monolayer).

surface nitrogen concentration increased, with values of 0.94, 1.12, 1.28, 1.76, and 2.15 eV. The variation in the energy barrier depends on the number of Fe−N bonds on the Fe surface, which depends on the electronic structure. Therefore, the electronic structure of the reaction is expected to correlate with 5313

dx.doi.org/10.1021/jp410947d | J. Phys. Chem. C 2014, 118, 5309−5316

The Journal of Physical Chemistry C

Article

inexpensive, Fe-based catalyst that can promote NH 3 decomposition for COx-free hydrogen production. Clarification of the exact mechanism and descriptions of the adsorption energy, the energy barrier, and the reaction energy of the NHx (with x = 0−3) intermediate species are important. We identified the exact catalytic decomposition reaction of NH3 on Fe(100). Subsequently, we investigated the adsorption of two types of ammonia dimers (2NH3) and concluded that decomposition probably occurs through locally adsorbed NH3 dimers. According to our results, the dehydrogenation of a hydrogen-bonded NH3 dimer leads to a higher energy barrier than that for monomeric NH 3 and that the second dehydrogenation reaction is the rate-limiting step. To analyze the energy barrier, we distributed the contribution factors of the elementary steps. Moreover, the investigation of the catalytic activity of an Fe surface precovered with N revealed a reduction in the catalytic activity as the N coverage on the Fe surface increased. In conclusion, the development of a method to remove nitrogen from the surface may facilitate the design of Fe catalysts with high activity for NH3 decomposition, thereby enabling the development of an Fe catalyst with similar performance as the most active Ru and Ir catalysts.

the adsorption energy and the energy barrier of the dissociation reaction. To confirm the effect of fully N-covered Fe(100) on the catalytic properties, we calculated the adsorption properties of NHx (with x = 0−3) intermediates on 1.00 ML N-covered Fe(100). Our results indicated that all species are adsorbed at the top site. The adsorption properties were worse than for clean Fe(100), as summarized in Table 4. These results indicate Table 4. Adsorption Properties of NHx (with x = 0−3) Intermediates on 1.00 ML N-Covered Fe(100)a species NH3 NH2 NH N H NH3 (molecule)

adsorption site

Ead (eV)

dN−H (Ǻ )

dN−Fe (Ǻ )

∠HNH (deg)

top top top top top

−0.36 −1.49 −1.73 −2.66

1.022 1.023 1.035

2.224 1.882 1.723 1.554

108.17 109.63

1.021

106.48

a

The energy was calculated with respect to the corresponding species in the gas phase and with respect to the NHx (with x = 0−3) species and the H atom.



ASSOCIATED CONTENT

S Supporting Information *

that the energy barrier of complete ammonia dehydrogenation was higher than for decomposition on a clean surface. In addition, the reaction energy (ΔE) indicates an endothermic reaction, as shown in Figure 7. Consequently, to improve the

The coordinates and thermodynamic data for the structural models reported in this paper. This material is available free of charge via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*E-mail [email protected]; Tel +82-42-350-3334; Fax +82-42350-3310 (H.M.L.). Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This research was supported by POSCO and a grant from the National Research Foundation of Korea (NRF) funded by the Korean government (MEST) (2011-0028612) and by the Future-based Technology Development Program (Nano Fields), which is funded through the NRF by the Ministry of Education, Science and Technology (2009-0082472). S.S.H. is grateful for the financial support from the Korea Institute of Science and Technology (Grant 2E24630).



Figure 7. Reaction pathway diagram of ammonia (NH3) decomposition on fully N-covered Fe(100).

REFERENCES

(1) Chein, R. Y.; Chen, Y. C.; Chang, C. S.; Chung, J. N. Numerical Modeling of Hydrogen Production from Ammonia Decomposition for Fuel Cell Applications. Int. J. Hydrogen Energy 2010, 35, 589−597. (2) Chellappa, A. S.; Fischer, C. M.; Thomson, W. J. Ammonia Decomposition Kinetics over Ni-Pt/Al2O3 for Pem Fuel Cell Applications. Appl. Catal., A 2002, 227, 231−240. (3) Lu, Y.; Wang, H.; Liu, Y.; Xue, Q. S.; Chen, L.; He, M. Y. Novel Microfibrous Composite Bed Reactor: High Efficiency H2 Production from NH3 with Potential for Portable Fuel Cell Power Supplies. Lab Chip 2007, 7, 133−140. (4) Metkemeijer, R.; Achard, P. Comparison of Ammonia and Methanol Applied Indirectly in a Hydrogen Fuel-Cell. Int. J. Hydrogen Energy 1994, 19, 535−542. (5) Appari, S.; Janardhanan, V. M.; Jayanti, S.; Maier, L.; Tischer, S.; Deutschmann, O. Micro-Kinetic Modeling of NH3 Decomposition on Ni and Its Application to Solid Oxide Fuel Cells. Chem. Eng. Sci. 2011, 66, 5184−5191.

catalytic performance of the Fe catalyst, a method to remove N surface coverage on the Fe surface should be considered. Several experiments7,24,62 examined adsorption and dissociation behaviors of the NH3 molecules and intermediates onto Fe(100), Fe(110), and Fe(111) surfaces and showed that as the concentration of chemisorbed nitrogen on Fe surfaces increases, the chemical activities of the Fe surfaces for NH3 dissociation are reduced, which is similar to our theoretical finding.

4. CONCLUSION In summary, we clarified the mechanism of the NH 3 decomposition reaction using DFT calculations to develop an 5314

dx.doi.org/10.1021/jp410947d | J. Phys. Chem. C 2014, 118, 5309−5316

The Journal of Physical Chemistry C

Article

Surface - Isolation and Identification of Adsorbed NH2 and NH. Sur. Sci. 1986, 175, 22−44. (28) Duan, X. Z.; Ji, J.; Qian, G.; Fan, C.; Zhu, Y.; Zhou, X. G.; Chen, D.; Yuan, W. K. Ammonia Decomposition on Fe(110), Co(111) and Ni(111) Surfaces: A Density Functional Theory Study. J. Mol. Catal. A: Chem. 2012, 357, 81−86. (29) Xiao, X. Z.; Cao, Y. L.; Cai, Y. Y. Decomposition of NH3 on Ir(110): A First-Principle Study. Surf. Sci. 2011, 605, 802−807. (30) Krekelberg, W. P.; Greeley, J.; Mavrikakis, M. Atomic and Molecular Adsorption on Ir(111). J. Phys. Chem. B 2004, 108, 987− 994. (31) Weststrate, C. J.; Bakker, J. W.; Rienks, E. D. L.; Martinez, J. R.; Vinod, C. P.; Lizzit, S.; Petaccia, L.; Baraldi, A.; Nieuwenhuys, B. E. Selective NH3 Oxidation on (110) and (111) Iridium Surfaces. J. Catal. 2005, 235, 92−102. (32) Huang, W. Y.; Lai, W. Z.; Xie, D. Q. First-Principles Study of Decomposition of NH3 on Ir(100). Surf. Sci. 2008, 602, 1288−1294. (33) Bradford, M. C. J.; Fanning, P. E.; Vannice, M. A. Kinetics of NH3 Decomposition over Well Dispersed Ru. J. Catal. 1997, 172, 479−484. (34) Dahl, S.; Sehested, J.; Jacobsen, C. J. H.; Tornqvist, E.; Chorkendorff, I. Surface Science Based Microkinetic Analysis of Ammonia Synthesis over Ruthenium Catalysts. J. Catal. 2000, 192, 391−399. (35) Logadottir, A.; Norskov, J. K. Ammonia Synthesis over a Ru(0001) Surface Studied by Density Functional Calculations. J. Catal. 2003, 220, 273−279. (36) Tsai, W.; Weinberg, W. H. Steady-State Decomposition of Ammonia on the Ru(001) Surface. J. Phys. Chem. 1987, 91, 5302− 5307. (37) Zhang, C. J.; Lynch, M.; Hu, P. A Density Functional Theory Study of Stepwise Addition Reactions in Ammonia Synthesis on Ru(0001). Surf. Sci. 2002, 496, 221−230. (38) Duan, X. Z.; Qian, G.; Zhou, X. G.; Sui, Z. J.; Chen, D.; Yuan, W. K. Tuning the Size and Shape of Fe Nanoparticles on Carbon Nanofibers for Catalytic Ammonia Decomposition. Appl. Catal., B 2011, 101, 189−196. (39) Lin, R. J.; Li, F. Y.; Chen, H. L. Computational Investigation on Adsorption and Dissociation of the NH3 Molecule on the Fe(111) Surface. J. Phys. Chem. C 2011, 115, 521−528. (40) Maier, S.; Stass, I.; Cerda, J. I.; Salmeron, M. Bonding of Ammonia and Its Dehydrogenated Fragments on Ru(0001). J. Phys. Chem. C 2012, 116, 25395−25400. (41) Diawara, B.; Joubert, L.; Costa, D.; Marcus, P.; Adamo, C. Ammonia on Ni(111) Surface Studied by First Principles: Bonding, Multilayers Structure and Comparison with Experimental IR and XPS Data. Surf. Sci. 2009, 603, 3025−3034. (42) Novell-Leruth, G.; Valcarcel, A.; Clotet, A.; Ricart, J. M.; PerezRamirez, J. Dft Characterization of Adsorbed NHx Species on Pt(100) and Pt(111) Surfaces. J. Phys. Chem. B 2005, 109, 18061−18069. (43) Wang, L.; Zhao, Y.; Liu, C. Y.; Gong, W. M.; Guo, H. C. Plasma Driven Ammonia Decomposition on a Fe-Catalyst: Eliminating Surface Nitrogen Poisoning. Chem. Commun. 2013, 49, 3787−3789. (44) Hellman, A.; Honkala, K.; Remediakis, I. N.; Logadottir, A.; Carlsson, A.; Dahl, S.; Christensen, C. H.; Norskov, J. K. Ammonia Synthesis and Decomposition on a Ru-Based Catalyst Modeled by First-Principles. Surf. Sci. 2009, 603, 1731−1739. (45) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B 1999, 59, 1758− 1775. (46) Kresse, G.; Furthmuller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169−11186. (47) Blochl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953−17979. (48) Kresse, G.; Hafner, J. Ab-Initio Molecular-Dynamics for OpenShell Transition-Metals. Phys. Rev. B 1993, 48, 13115−13118. (49) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868.

(6) Avci, A. K.; Trimm, D. L.; Onsan, Z. I. Quantitative Investigation of Catalytic Natural Gas Conversion for Hydrogen Fuel Cell Applications. Chem. Eng. J. 2002, 90, 77−87. (7) Grunze, M.; Bozso, F.; Ertl, G.; Weiss, M. Interaction of Ammonia with Fe(111) and Fe(100) Surfaces. Appl. Surf. Sci. 1978, 1, 241−265. (8) Hansgen, D. A.; Vlachos, D. G.; Chen, J. G. G. Using First Principles to Predict Bimetallic Catalysts for the Ammonia Decomposition Reaction. Nat. Chem. 2010, 2, 484−489. (9) Choudhary, T. V.; Sivadinarayana, C.; Goodman, D. W. Catalytic Ammonia Decomposition: COx-Free Hydrogen Production for Fuel Cell Applications. Catal. Lett. 2001, 72, 197−201. (10) Choudhary, T. V.; Santra, A. K.; Sivadinarayana, C.; Min, B. K.; Yi, C. W.; Davis, K.; Goodman, D. W. Ammonia Decomposition on Ir(100): From Ultrahigh Vacuum to Elevated Pressures. Catal. Lett. 2001, 77, 1−5. (11) Ganley, J. C.; Thomas, F. S.; Seebauer, E. G.; Masel, R. I. A Priori Catalytic Activity Correlations: The Difficult Case of Hydrogen Production from Ammonia. Catal. Lett. 2004, 96, 117−122. (12) Jacobsen, C. J. H.; Dahl, S.; Clausen, B. S.; Bahn, S.; Logadottir, A.; Norskov, J. K. Catalyst Design by Interpolation in the Periodic Table: Bimetallic Ammonia Synthesis Catalysts. J. Am. Chem. Soc. 2001, 123, 8404−8405. (13) Boisen, A.; Dahl, S.; Norskov, J. K.; Christensen, C. H. Why the Optimal Ammonia Synthesis Catalyst Is Not the Optimal Ammonia Decomposition Catalyst. J. Catal. 2005, 230, 309−312. (14) Erisman, J. W.; Sutton, M. A.; Galloway, J.; Klimont, Z.; Winiwarter, W. How a Century of Ammonia Synthesis Changed the World. Nat. Geosci. 2008, 1, 636−639. (15) Klerke, A.; Christensen, C. H.; Norskov, J. K.; Vegge, T. Ammonia for Hydrogen Storage: Challenges and Opportunities. J. Mater. Chem. 2008, 18, 2304−2310. (16) Garcia-Garcia, F. R.; Guerrero-Ruiz, A.; Rodriguez-Ramos, I. Role of B5-Type Sites in Ru Catalysts Used for the NH 3 Decomposition Reaction. Top. Catal. 2009, 52, 758−764. (17) Campbell, C. T. Finding the Rate-Determining Step in a Mechanism - Comparing Dedonder Relations with the “Degree of Rate Control”. J. Catal. 2001, 204, 520−524. (18) Dahl, S.; Logadottir, A.; Egeberg, R. C.; Larsen, J. H.; Chorkendorff, I.; Tornqvist, E.; Norskov, J. K. Role of Steps in N-2 Activation on Ru(0001). Phys. Rev. Lett. 1999, 83, 1814−1817. (19) Frechard, F.; van Santen, R. A.; Siokou, A.; Niemantsverdriet, J. W.; Hafner, J. Adsorption of Ammonia on the Rhodium (111), (100), and Stepped (100) Surfaces: An Ab Initio and Experimental Study. J. Chem. Phys. 1999, 111, 8124−8130. (20) Mortensen, H.; Diekhoner, L.; Baurichter, A.; Jensen, E.; Luntz, A. C. Dynamics of Ammonia Decomposition on Ru(0001). J. Chem. Phys. 2000, 113, 6882−6887. (21) Egawa, C.; Nishida, T.; Naito, S.; Tamaru, K. Ammonia Decomposition on (1 1 10) and (0 0 1) Surfaces of Ruthenium. J. Chem. Soc., Faraday Trans. I 1984, 80, 1595−1604. (22) Deshmukh, S. R.; Mhadeshwar, A. B.; Vlachos, D. G. Microreactor Modeling for Hydrogen Production from Ammonia Decomposition on Ruthenium. Ind. Eng. Chem. Res. 2004, 43, 2986− 2999. (23) Lanzani, G.; Laasonen, K. NH3 Adsorption and Dissociation on a Nanosized Iron Cluster. Int. J. Hydrogen Energy 2010, 35, 6571− 6577. (24) Ertl, G.; Huber, M. Mechanism and Kinetics of Ammonia Decomposition on Iron. J. Catal. 1980, 61, 537−539. (25) Seabury, C. W.; Rhodin, T. N.; Purtell, R. J.; Merrill, R. P. Chemisorption and Reaction of NH3 on Ni(111). Surf. Sci. 1980, 93, 117−126. (26) Stolbov, S.; Rahman, T. S. First-Principles Study of Some Factors Controlling the Rate of Ammonia Decomposition on Ni and Pd Surfaces. J. Chem. Phys. 2005, 123, 204716−1. (27) Bassignana, I. C.; Wagemann, K.; Kuppers, J.; Ertl, G. Adsorption and Thermal-Decomposition of Ammonia on a Ni(110) 5315

dx.doi.org/10.1021/jp410947d | J. Phys. Chem. C 2014, 118, 5309−5316

The Journal of Physical Chemistry C

Article

(50) Mortensen, J. J.; Ganduglia-Pirovano, M. V.; Hansen, L. B.; Hammer, B.; Stoltze, P.; Norskov, J. K. Nitrogen Adsorption on Fe(111), (100), and (110) Surfaces. Surf. Sci. 1999, 422, 8−16. (51) Logadottir, A.; Norskov, J. K. The Effect of Strain for N2 Dissociation on Fe Surfaces. Surf. Sci. 2001, 489, 135−143. (52) Jiang, D. E.; Carter, E. A. Diffusion of Interstitial Hydrogen into and through Bcc Fe from First Principles. Phys. Rev. B 2004, 70, 064102. (53) Jiang, D. E.; Carter, E. A. Carbon Atom Adsorption on and Diffusion into Fe(110) and Fe(100) from First Principles. Phys. Rev. B 2005, 71, 045402. (54) Jiang, D. E.; Carter, E. A. First Principles Study of H2S Adsorption and Dissociation on Fe(110). Surf. Sci. 2005, 583, 60−68. (55) Jiang, D. E.; Carter, E. A. Carbon Dissolution and Diffusion in Ferrite and Austenite from First Principles. Phys. Rev. B 2003, 67, 214103. (56) Henkelman, G.; Uberuaga, B. P.; Jonsson, H. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901−9904. (57) Tang, W.; Sanville, E.; Henkelman, G. A Grid-Based Bader Analysis Algorithm without Lattice Bias. J. Phys.: Condens. Matter 2009, 21, 084204. (58) Weiss, M.; Ertl, G.; Nitschke, F. Adsorption and Decomposition of Ammonia on Fe(110). Appl. Surf. Sci. 1979, 2, 614−635. (59) Yeo, S. C.; Han, S. S.; Lee, H. M. Adsorption, Dissociation, Penetration, and Diffusion of N2 on and in Bcc Fe: First-Principles Calculations. Phys. Chem. Chem. Phys. 2013, 15, 5186−5192. (60) Liu, Z. P.; Hu, P.; Alavi, A. Catalytic Role of Gold in Gold-Based Catalysts: A Density Functional Theory Study on the Co Oxidation on Gold. J. Am. Chem. Soc. 2002, 124, 14770−14779. (61) Shustorovich, E.; Bell, A. T. Synthesis and Decomposition of Ammonia on Transition-Metal Surfaces - Bond-Order-ConservationMorse-Potential Analysis. Surf. Sci. 1991, 259, L791−L796. (62) Ertl, G. Surface Science and Catalysis - Studies on the Mechanism of Ammonia-Synthesis - the Emmett, P.H. Award Address. Catal. Rev.Sci. Eng. 1980, 21, 201−223. (63) Hammer, B. Coverage Dependence of N2 Dissociation at an N, O, or H Precovered Ru(0001) Surface Investigated with Density Functional Theory. Phys. Rev. B 2001, 63, 205423. (64) Egawa, C.; Naito, S.; Tamaru, K. Adsorption and Decomposition of NO and NH3 on Ru(001) and Ru(1,1,10) Surfaces. Surf. Sci. 1984, 138, 279−291. (65) van Steen, E.; van Helden, P. A DFT Study of Hydrogen Dissociation on Co- and C-Precovered Fe(100) Surfaces. J. Phys. Chem. C 2010, 114, 5932−5940.

5316

dx.doi.org/10.1021/jp410947d | J. Phys. Chem. C 2014, 118, 5309−5316