Next Article in Journal
Selective Catalytic Reduction of NO by NH3 Using a Combination of Non-Thermal Plasma and Mn-Cu/ZSM5 Catalyst
Next Article in Special Issue
CeZrOx Promoted Water-Gas Shift Reaction under Steam–Methane Reforming Conditions on Ni-HTASO5
Previous Article in Journal
Synthesis and Characterization of Zinc Peroxide Nanoparticles for the Photodegradation of Nitrobenzene Assisted by UV-Light
Previous Article in Special Issue
A Review on Catalysts Development for Steam Reforming of Biodiesel Derived Glycerol; Promoters and Supports
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Catalysts
Volume 10
Issue 9
10.3390/catal10091043
catalysts-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanistic Insights for Dry Reforming of Methane on Cu/Ni Bimetallic Catalysts: DFT-Assisted Microkinetic Analysis for Coke Resistance

by
Ahmed Omran
1,
Sun Hee Yoon
1,
Murtaza Khan
1,
Minhaj Ghouri
1,
Anjaneyulu Chatla
1 and
Nimir Elbashir
1,2,*
1
Chemical Engineering Program, Texas A&M University at Qatar, Doha 23874, Qatar
2
TEES Gas & Fuels Research Center, Texas A&M University, College Station, TX 77843, USA
*
Author to whom correspondence should be addressed.
Catalysts 2020, 10(9), 1043; https://0-doi-org.brum.beds.ac.uk/10.3390/catal10091043
Submission received: 22 August 2020 / Revised: 7 September 2020 / Accepted: 8 September 2020 / Published: 10 September 2020
(This article belongs to the Special Issue Catalytic Steam Reforming)
Browse Figures
Versions Notes

Abstract

:
Density functional theory (DFT) calculations have been utilized to evaluate the complete reaction mechanism of methane dry reforming (DRM) over Ni2Cu (111) bimetallic catalyst. The detailed catalytic cycle on Ni2Cu (111) catalyst demonstrated superior coke resistance compared to pure Ni (111) and Ni2Fe (111) reported in the literature. Doping Cu in the Ni–Ni network enhanced the competitive CH oxidation by both atomic O and OH species with the latter having only 0.02 eV higher than the 1.06 eV energy barrier required for CH oxidation by atomic O. Among the C/CH oxidation pathways, C* + O* → CO (g) was the most favorable with an energy barrier of 0.72 eV. This was almost half of the energy barrier required for the rate-limiting step of CH decomposition (1.40 eV) and indicated enhanced coke deposition removal. Finally, we investigated the effect of temperature (800~1000 K) on the carbon deposition and elimination mechanism over Ni2Cu (111) catalyst. Under those realistic DRM conditions, the calculations showed a periodic cycle of simultaneous carbon deposition and elimination resulting in improved catalyst stability.

1. Introduction

Dry reforming of methane (DRM) has been regarded as a process to produce syngas (CO and H2), while providing a prospective solution for mitigation of greenhouse gases (CH4 and CO2). The produced syngas is a valuable precursor for many other downstream applications [1,2]. It can be used directly as fuel or converted to other hydrocarbons to produce valuable ultra-clean fuels and products such as methanol, diesel, gasoline, kerosene and naphtha via Fischer-Tropsch synthesis, and other reactions [3,4]. In addition to syngas, there are ongoing efforts to make the process more economical by producing valuable byproducts such as solid carbon as added value [5,6]. Despite its numerous advantages, DRM is a highly endothermic process (Equation (1)) that requires high temperatures of around 700 °C to achieve equilibrium conversions of CH4 and CO2 [7]. Therefore, the reaction proceeds typically with the use of a catalyst.
CH4 + CO2 ⇌ 2CO + 2H2    ∆H289 = 247.3 kJ mol−1
Catalyst deactivation, commonly caused by a rapid accumulation of carbon or “coking,” is one of DRM’s major drawbacks. Metal catalysts such as platinum (Pt) [8], rhodium (Rh) [9], ruthenium (Ru) [10], iridium (Ir) [11], and cobalt (Co) [12] have been used to overcome the drawbacks of coking, as well as the unfavorable reactivity of the formed carbon towards oxidation reactions. Noble metals, such as Ru and Rh, have improved catalytic activity, stability, and coke resistance [3]. However, those metals are expensive and less abundant, making their usage economically unattractive for large industrial applications [13]. On the other hand, nickel catalysts have proved catalytic activity comparable to that of noble metals and can be a cheaper alternative with greater availability. However, the use of Ni catalysts in DRM is still limited, due to their high thermodynamic potential to coke formation [14,15].
The recent developments of bimetallic Ni-transition metal catalysts [16,17] have attracted attention as an effective strategy to overcome the drawbacks of pure Ni catalysts [18]. Xu et al. [19] conducted a DFT study on DRM for Ni2Fe overlayer of Ni (111) surface. They concluded that the surface carbon reacted with the lattice oxygen in FeO and produced CO with a low energy barrier of 0.16 eV via Mars-van Krevelen (MvK) mechanism. Their conclusion is consistent with the experimental results from Theofanidis et al. [20,21].
Cu, another example of the alloyed transition metals, acted as an active phase promoter in CH4/CO2 reforming. The incorporation of Cu into the Ni–Ni surface has been demonstrated to improve stability, activity, and coke resistance in DRM. According to Faro et al. [22], the partial substitution of Ni by Cu interrupts the Ni-Ni network and reduces coke formation by altering their crystal and surface features. Rahemi et al. [23] synthesized Ni–Cu/Al2O3 nanocatalysts to explain the effect of Cu. They found that the size of the NiO crystals affected the stability of the nanocatalysts and their activities are highly dependent on the morphology and synergies between nickel and copper. This synergistic effect is very closely related to the amount of Cu loading [23]. For instance, Wu et al. found that Cu/Ni showed a higher activity at a 1:3 loading ratio, but the higher the Cu content, the significantly lower the activation of methane. This may be due to the competitive adsorption of CH4 and slow dissociation kinetics at the copper atoms [24]. These outcomes are supported by Wang et al., who investigated the effect of coinage metals doped on the Ni-Ni network for DRM and concluded that the adsorption energies of CHx (x = 0–3) are reduced on the alloy surface [18]. Moreover, they also agreed with Liu et al. [25], who performed DFT study on CH4 dissociation over a Ni–Cu (111) catalyst surface. Their study outcomes clearly showed that the Cu-rich surface induced better coke resistance.
Understanding the CH4 dissociation pathway is very important to evaluate both the catalyst activity and coke resistance. An et al. [26] found that the activation energy barrier for CH → C + H and CH4 → CH3 + H on NiCu (111) in the DRM mechanism were 1.8 times and 1.3 times higher than Ni (111), respectively. Similarly, Wei et al. [7] studied methane dehydrogenation on the bimetallic Cu/Ni (111) surface, and observed that coke formation was reduced by increasing the CH energy barrier. Very recently, our research team (Chatla et al. [27]) conducted a DFT study to investigate the effect of Cu addition to Ni/Al2O3 catalysts on DRM mechanism. Our findings showed that the addition of Cu increases the activation energy barrier for CHx dehydrogenation to carbon nearly 1.5 times, but considerably decreases the activation energy barrier for carbon elimination through carbon gasification of Ni–Cu alloy system. However, the stability of Cu/Ni alloy cannot be sufficiently explained by just the ensemble effect of Cu on the Ni metal, which can suppress the rate of carbon formation. More comprehensively, the incorporation effect of Cu into Ni on the oxidation pathway of C/CH should be considered. This factor would be more important in the stability control of Cu–Ni alloy.
In this study, density functional theory (DFT) calculations are employed to consider how the Cu-alloyed Ni (111) catalyst affects the complete catalytic cycle of DRM, focusing on the essential factors that can control the Ni-Cu catalyst ability to resist coke formation. We investigated all possible coke deposition and elimination pathways in realistic reaction conditions over Ni2Cu (111) bimetallic system. In particular, both the effects of temperature and Cu doping on coke deposition resistance are explained by calculating the reaction rate constant over the temperature range of 800–1000 K through complete vibration analysis of the elementary reactions.

2. Results and Discussion

2.1. Adsorption Geometries and Energies on Ni2Cu Overlayer of Ni (111) Surface

The electronic properties of Ni2Cu overlayer of Ni (111) surface described in Figure 1 were investigated, and three different adsorption sites were identified: two top sites (TopNi and TopCu), two bridge sites (B2Ni and BNiCu), and two hollow sites (Hexagonal Close Packed, HCP and Face Centered Cubic, FCC).
First, optimal adsorption sites for all DRM reaction species were selected after examining all possible adsorption structures. Then, the adsorption energies of all DRM reaction species from the most stable adsorption configurations were calculated, as shown in Table 1. CH4 species was physically adsorbed at the HCP site with a weak adsorption energy of −0.02 eV.
On the other hand, CH3, CH2, and CH are strongly adsorbed via C atom on the FCC site with corresponding adsorption energies of −2.17 eV, −4.99 eV, and −5.43 eV, respectively. Similarly, the H atom tends to be adsorbed to the FCC site. These results are in good agreement with Ni2Fe overlayer [20,21] with same diagonal arrangement and pure Ni (111) surfaces [28]. In the case of CO, the adsorption is preferable on the TopNi site with Eads of −1.52 eV, while the O atom is favorably bound to the FCC site with an Eads of −6.34 eV. It indicates that the Cu dopant in Ni (111) network improved atomic O adsorption while reducing CO adsorption compared to both pure Ni (111) and Ni2Fe (111) surfaces, respectively [19]. As a descriptor of coke resistance, the adsorption energy of the C atom indicates that the more strongly the carbon is adsorbed on a specific surface, the lower the coke resistance. The C atom preferentially adsorbed to HCP site with an Eads of −7.04 eV and is less than the Eads of −7.50 eV in a monometallic Ni catalyst [27,29,30]. This means that the Cu dopant on Ni (111) reduces coke deposition, and, as a result, it can increase the coke resistance of the catalyst surface. The adsorption of CO2 on the Ni2Cu (111) surface has been investigated at all possible sites. After geometry optimization, it has been found that CO2 is physically adsorbed towards a vacuum with weak adsorption energy of −0.03 eV in parallel to the surface with C–Cu. This result agrees with Heil and [31], who used both fast ion beam and fast molecular beam techniques to study CO2 dissociation, and found no chemisorbed CO2 on Ni (111). Moreover, it agrees with the computational studies performed on Ni (111) and on Ni2Fe (111) surfaces [19,32]. CO molecules prefer to adsorb at the B2Ni site with an Eads of −1.52 eV. The COH* intermediate produced by the combination of C* and OH* is strongly adsorbed at the B2Ni site with an Eads of −4.07 eV. In the case of H2O, it is physically adsorbed to the surface with Eads of −0.06 eV. The adsorption sites and energies of the other DRM intermediates are detailed in Table 1.

2.2. DRM Reaction Mechanism

DRM turnover rate is known to be limited only by C–H bond activation, regardless of the concentration of reactants [7]. Consequently, it is expected that the CH4 is consumed via dissociative adsorption regardless of CO2 concentration. The resulting CHx species are then oxidized by atomic O and OH species to form CHxO and CHxOH species, respectively, or they undergo further dehydrogenation.

2.2.1. CHx Dissociation (x = 1–4)

The successive steps of CH4 dissociation start with the physical adsorption of CH4 (g) on the Ni2Cu (111) surface to form CH4*(* denotes an adsorption state). Then, the activation energy in the first step of CH4 dissociation is found to be 1.30 eV, as shown in Figure 2. The first transition state shows the detached H atom located at the HCP site [33], and the remaining CH3 fragment settled down at 2.08 Å to achieve maximum C–H–Ni three-center bonding [34,35]
After the dissociative adsorption of CH4 on the Ni2Cu (111) surface, the two subsequent dehydrogenation steps of CH3* and CH2* species are more facile with energy barriers of 0.75 eV and 0.49 eV, respectively. However, CH3* dehydrogenation showed endothermicity of 0.33 eV, while CH2* exhibited exothermicity by −0.06 eV. On the other hand, oxidations of CH3* and CH2* by O* or OH* are considerably difficult, due to relatively higher energy barriers than the dehydrogenation pathway mentioned above.
Figure 3 shows the activation energy barriers of all possible pathways of CH3*. It confirms that subsequent dissociation to CH2* is dominant over oxidation pathways.CH2* and CH3* oxidation can produce CH3OH* through CH3*, CH3O*, or CH2OH* intermediates. Of these intermediates, CH2OH* was found to be more favorable with a lower activation energy (Ea,r) of 1.10 eV, as shown in Table 2. In addition, Table 2 shows that the activation energy of the forward reaction (Ea,f) and the activation energy for the reverse reaction (Ea,r). Both of these reaction rates are calculated from Climbing-NEB method considering the enthalpy (Hf) as the energy difference between the final state (F.S) and the initial state (I.S). Next, for the CH dissociation, the activation energy (Ea,f) increases significantly to 1.40 eV, which is higher than that of CH* oxidation by OH* (1.08 eV) and O* species (1.06 eV), as shown in Table 2. The conclusion of this result is that CH dissociation is the rate-limiting step.
Consequently, the energy barriers for CHx (x = 1–4) dissociation increase in order of CH2 < CH3 < CH4 < CH [36], which is in line with CH4 decomposition on pure Ni (111) and Ni2Fe (111) systems. In addition, the energy barriers of those subsequent steps are much higher in the case of Ni2Cu (111) than Ni (111) and Ni2Fe (111) systems, which reduce carbon deposition from CH radical dehydrogenation [19,28,37].

2.2.2. Two Reaction Pathways for CO2 Dissociation

According to the previous experimental studies [7,38], there are two possible routes for CO2 dissociation. First, adsorbed CO2* is directly dissociated to form CO* and atomic O*. The latter acts as an oxidant for CHx (x = 1–4), ultimately producing CO (g). In the second pathway, CH4* is first dehydrogenated to produce atomic H, which subsequently activates the CO2* to produce a COOH* intermediate. Then, the decomposition of COOH* produces CO* and an OH* oxidant in the following scheme.
CO2* + H* → COOH* → CO* + OH*
As described above, H* atoms released from the CH4 dissociation pathway can either contribute towards or activate CO2 dissociation. The potential energy profile and structure associated with CO2 hydrogenation (or H-induced CO2 dissociation) are shown in Figure 4. At the first step, CO2 and H atoms are co-adsorbed on Ni2Cu (111) surface to form a COOH* intermediate with an Ea,f of 0.93 eV, and this reaction is exothermic by −0.11 eV. In the final state, CO and OH are co-adsorbed at the B2Ni site after the breakage of COOH*. For adsorbed CO*, the C-Ni bond length is 1.89 Å, and in OH*, the O–Ni bond length is 1.98 Å. The activation energy (Ea,f) of COOH* dissociation is 0.43 eV, and the reaction is exothermic by 0.59 eV. Previous studies of DRM on pure Ni (111) surfaces have shown that direct CO2 dissociation has a low energy barrier of 0.67 eV compared to 1.13 eV for the H-induced CO2 dissociation pathway [28]. Meanwhile, the Ni2Cu (111) surface has a lower (Ea,f) of 0.93 eV for H-induced CO2 dissociation. Thus, doping Cu in Ni (111) network slightly favors the CO2 hydrogenation process, producing an OH oxidant via COOH* intermediate dissociation. Comparing the direct CO2 dissociation and the H-induced CO2 dissociation pathways, the latter is more favorable reaction in Ni2Cu (111) bimetallic system than in the case of Ni (111), where both atomic O* and OH* species participate as oxidants of CHx (x = 1–3) intermediates. Similar results have been reported for the Ni2Fe (111) surface, where both pathways were found to be competitive [19].

2.2.3. Oxidation of C and CH

When CH4 dissociates into C and H, CHx (x = 0–3) can combine with atomic O or OH species to produce CHO (H). Table 2 shows that the CH3 and CH2 dehydrogenation reactions are more favorable with relatively lower Ea,f of 0.75 eV and 0.49 eV, respectively. In contrast, CH dehydrogenation has a relatively higher Ea,f of 1.40 eV, while the reverse reaction is energetically favorable with Ea,r of 0.71 eV. Thus, coke formation due to CH decomposition is less desirable for Ni2Cu (111) compared to pure Ni (111) and Ni2Fe (111) surfaces with energy barriers of 1.33 eV [19] and 1.36 eV [28], respectively. Moreover, only CH* and C* are the most stable major intermediates remaining on the surface than the kinetically unstable CH3* and CH2* intermediates, as discussed in Table 1. It turns out that from all possible oxidation reactions of CHx (x = 0–3), C* + O* oxidation is the most favorable reaction.

2.2.4. Carbon elimination by C + O and C + OH

The C* + O* oxidation pathway was found to be more favorable at an Ea,f of 0.72 eV compared to 1.59 eV and 1.58 eV for Ni (111) and Ni2Fe (111), respectively, as shown in Figure 5.
The C* + OH* reaction pathway shows an Ea,f of 1.13 eV compared to 1.46 eV for Ni (111) and 1.42 eV for Ni2Fe (111) surfaces. The reverse reaction of CO* → C* + O* has an Ea,r of 2.89 eV, which is clearly in favor of carbon elimination. These results indicate that Cu-rich surfaces can accelerate specific oxidation reactions and change surface properties to improve carbon elimination by both atomic O and even OH species. As will be explained later, temperature plays a crucial role in enhancing these oxidation steps. These results conclude that Cu doped in Ni (111) reduced coke formation by increasing the activation energy barrier of CO* and CH* dissociation, but improved carbon elimination by both atomic O* and OH* species [18,39].
This study shows that the addition of Cu can overcome the limitation of Ni (111) and Ni2Fe (111) in terms of coke resistance, which is consistent with the results of previous calculations [26] and experimental studies [37].

2.2.5. CH + O and CH + OH Reactions

Oxidation of CH* by atomic O* generated from direct CO2 dissociation begins through co-adsorption of CH on the FCC site and an O atom on an HCP site via one Cu atom. In the transition state, atomic O moves to the BNiCu site at a Ni–O distance of 1.97 Å and a Cu–O distance of 2.01 Å. CH* settles at the B2Ni site at a Ni–C distance of 1.86 Å and a Ni-H distance of 1.81 Å. At this step, the energy barrier is 1.06 eV, which is much higher than 0.47 eV for pure Ni (111) and 0.78 eV for the Ni2Fe (111) surface [19,28].
Figure 6 shows that the CHOH* intermediate formed by co-adsorption of CH* and OH* at adjacent FCC and B2Ni sites easily overcomes the energy barrier of 1.08 eV compared to 1.48 eV and 1.42 eV in pure Ni (111) and Ni2Fe (111), respectively [19,28].

2.2.6. CHOH and COH Decomposition

A stable CHOH* configuration was found at the B2Ni site with the same C-Ni distance of 1.95 Å. Figure 7 shows two possible pathways for subsequent CHOH*. The first pathway involves CHOH* dissociation (CHOH* → CHO* + H*) with an energy barrier of 0.60 eV. Then, the produced CHO* has a low energy barrier of 0.18 eV, which readily forms CO* and atomic H*. The second pathway follows CHOH* dissociation, which forms COH* and H* with an energy barrier of 0.20 eV, and then the COH* dissociates to CO* and atomic H* with an energy barrier of 1.09 eV.
From the results above, it can be concluded that the CH oxidation by atomic O (1.06 eV) and OH oxidants (1.08 eV) in the Ni2Cu (111) system is very competitive, due to the small difference of 0.02 eV. Moreover, both oxidation pathways have a lower energy barrier than that of Ni (111) and Ni2Fe (111) [19,28].
The energy barrier of CH oxidation by atomic O was slightly lower than that of CH oxidation by OH. Conversely, the oxidation of CH by OH showed a higher reaction rate constant in the temperature range of 800~1000K, which will be described in detail in Section 2.3.

2.2.7. H2 and H2O Formation

The reaction of H* + H* → H2 (g) requires an Ea,f of 0.73 eV, and is an endothermic reaction at a Hf of 0.57 eV, as shown in Table 2. The H2 (g) is physically adsorbed on the catalyst surface with an Eads of −0.016 eV, which prefers the product formation in a higher reaction temperature. Hence, H2 (g) formation in Ni2Cu (111) catalyst is more facile than in pure Ni (111) and Ni2Fe (111), which have relatively higher energy barriers of 0.92 eV and 1.18 eV, respectively. In the first step of H2O* formation, the adsorbed H* and O* react to produce OH* with an Ea,f of 0.94 eV and exothermic reaction energy (Hf) of −0.24 eV, as shown in Table 2. Figure 8, which describes the next step, shows that the formed OH* undergoes an exothermic reaction with H* with an energy barrier of 0.91 eV. These two steps are energetically preferred over pure Ni (111) and Ni2Fe (111) surfaces [19,28]. At high temperatures in DRM, H2O (or steam) facilitates the physical removal of coke deposition from the catalyst surface.

2.3. Effect of Temperature on Carbon Deposition Resistance

One drawback of the DFT calculations used to investigate the DRM reaction is that all atoms are localized to the “minimum energy” position at zero temperature (0 K). A complete vibrational analysis is performed on the elementary reaction to overcome this problem. The analysis has two advantages; the first is that the transition state configuration as well as the activation energy barrier can be verified with a single imaginary frequency. The second advantage is that the rate constants can be calculated through the Vineyard formula [38], a crucial step in developing microkinetic models and identifying atomic models [40,41].
In this regard, the calculated vibrational frequencies are used to evaluate the rate constant at the confirmed transition states, focusing on the effect of Cu doping and temperature on coke deposition resistance. The calculations of the vibrational frequency and rate constant are described in more detail in the supplementary information. Coke deposition can come from two primary sources: CH4 and CO2.Inthe CH4 dissociation pathway, the ratio of the rate constant of the CH oxidation pathway to the C oxidation pathway determines the carbon deposition. According to a previous study, increasing this ratio reduces carbon formation [29]. Coke can be also formed during the dissociation of CO into C* and O* in the CO2 reaction pathway. Table 3 shows the rate constants k (s−1) of carbon elimination and deposition on the Ni2Cu (111) surface over various temperature ranges. We found that k of C* + O* (or + OH*) or CH* + O* (or + OH*) oxidation reactions increases with the increasing temperature, and the carbon deposition from CO* is negligible due to the very small k value. In addition, the rate constant of CH dissociation was lower than that of carbon oxidation (kC(O) and kC(OH)) by atomic O or OH species, respectively. As the temperature increases, kCH/kC(O) increases while kCH/kC(OH) decreases.
Figure 9 shows the temperature effect on carbon elimination reactions. The reaction of C* + O* has the highest rate constant in all carbon elimination reactions, followed by CH* + OH*. C* + OH* and CH* + O* reactions have similarly low rate constants as a function of temperature. OH* is the primary oxidizing agent for CH*, whereas atomic O is dominant in coke removal at different temperatures. It indicates that coke deposition by atomic O is easily removed at higher temperature ranges. On the other hand, it also reveals the importance of OH species in bypassing the coke formation step from CH dissociation where kCH/kC(OH) >> kCH/kC(O), as shown in Figure 10. Comparing these results with previous studies [42,43], we can conclude that carbon elimination via C* + O* is enhanced in Ni2Cu (111) system over the Ni2Fe (111) and pure Ni (111) systems [19,25,28]. It suggests that the incorporation of Cu to Ni–Ni surface can simultaneously improve carbon elimination by C* + O* reaction, while reducing carbon deposition via CH dissociation. The higher the Cu/Ni ratio, the more it is expected to suppress catalyst activity when occupying the Ni active site. This result is consistent with the experimental results [22,24].
Since the carbon elimination rate by both OH and atomic O is much higher than that of carbon deposition through CH dissociation, the Ni2Cu (111) surface is expected to exhibit higher catalytic activity. In addition, this periodic cycle of carbon deposition and elimination will result in catalyst stability.

2.4. Dominant Reaction Pathways and Rate-Limiting Steps

After determining the initial and final state of each elementary reaction of DRM, transition states are identified and confirmed by vibrational analysis. Then, zero-point energy corrections and rate constants are evaluated. As a result, the dominant reaction pathway and rate-limiting step are determined. CH4 undergoes successive dehydrogenation steps, of which CH3 and CH2 dissociation more readily occur with lower energy barriers of 0.75 eV and 0.49 eV, respectively. However, CH3 and CH2 oxidation by atomic O or OH are energetically unfavorable due to the high energy barriers. To illustrate, CH3 requires energy barriers of 1.37 eV and 1.6 eV while CH2 needs 0.76 eV and 0.78 eV for oxidation by O and OH, respectively.
Consequently, CH4 is dissociated to CH, C, and H intermediates. The main intermediate of CH4 dissociation, CH, follows three primary pathways: (a) CH* → C* + H*, (b) CH* + O* → CHO* → CO* + H*, and (c) CH* + OH* → CHOH*. The CH dissociation, as a significant source of coke formation, is energetically and kinetically unfavorable. Therefore, CH dissociation is the rate-limiting step of CH4 dissociation and determines the overall reaction rate. In the case of CH* oxidation pathways, both (b) and (c) are very competitive, with only a 0.02 eV difference and (b) slightly preferred with a lower energy barrier of 1.06eV. However, the kinetic data shown in Table 3 indicate that CHOH* formation is more favorable reaction than CHO* at the temperature range 800~1000 K. Concerning CHOH*, the most preferred pathway goes through COH* formation with an energy barrier of 0.20 eV, which further dissociates into CO* and H*. From the standpoint of coke elimination, C oxidation by atomic O is a more favorable reaction with an energy barrier of 0.72 eV, which is lower than all other C and CH oxidation pathways. As for H* intermediate, H* + H* → H2* is the most preferred pathway. In the CO2 activation pathway, we investigated the direct CO2 dissociation to form CO and O, and the H-induced CO2 formation pathway to form CO* and OH*. Unlike pure Ni (111) system, the H-induced CO2 formation pathway is the more favorable than the direct CO2 dissociation with a lower energy barrier of 0.93 eV. Eventually, this result suggests that the doped Cu atoms promote the H-induced CO2 formation pathway and form O and OH oxidants that remove CH and C from the Ni2Cu (111) surface, as shown in Figure 10.

3. Computational Details

The first-principles calculations performed in this study are based on spin-polarized DFT implemented in the Vienna Ab initio Simulation Package (VASP) [28,44]. Blöchl’s all-electron-like projector augmented wave (PAW) method [45] was employed to explain the interactions between the ion cores and the valence electrons. The exchange-correlation of the Kohn–Sham theory are treated by the generalized gradient approximation (GGA) [46] with Revised Perdew-Burke-Ernzerhof (RPBE), which gives more accurate chemisorption energies than the PW91 and PBE functional [47,48]. A plane-wave basis set with a kinetic energy cut-off of 400 eV is used for DFT calculations. All geometries were optimized using a force-based conjugate gradient algorithm until the forces acting on each atom converge better than 0.01 eV/Å. The optimized lattice constant for Ni is 3.523 Å, which is consistent with the experimental value of 3.524 Å [19]. Brillouin zone samplings are defined using Monkhorst-Pack scheme [39] and electronic occupancies are determined in a Methfessel-Paxton smearing method [49] with smearing energy of 0.2 eV.
Previous DFT studies on Ni catalysts considered the spin polarization of the surface magnetism to obtain more accurate total energies [50]. The Ni2Cu (111) surface is represented as a four-layer slab with a p(3 × 3) supercell where only the bottom layer of the slab is constrained. The neighboring slabs are separated by a vacuum region of 12 Å. The first Brillouin zone of the supercell is sampled with a 3 × 3 × 1 k-point grid, which is sufficient for this cell structure. The climbing nudged elastic band [51,52] and dimer methods [53] are used to locate the transition state (TS). Total energy and band structure energy are converged to obtain accurate forces within 1 × 10−7 eV/atom during the electronic optimization. The adsorption energy (Eads) is calculated as follows:
Eads = E(adsorbate/slab) − Eslab − Eadsorbate
where E(adsorbate/slab) is the total energy of the slab model with an adsorbate, Eslab is the total energy of the bare slab model, Eadsorbate is the total energy of the isolated adsorbate. As agreed in the literature [54], the more negative adsorption energy indicates that the gas species is more strongly adsorbed on the surface.
The Ni structure is crystallized in a face-centered cubic (FCC) crystal lattice. We selected the (111) orientation of this lattice due to their abundance in nickel-based alloys [50]. The Cu/Ni (111) surface alloy was then constructed by replacing three nickel atoms in the topmost layer of Ni (111) by copper atoms (i.e., at the dopant coverage of 1/3 ML). The ratio (1/3 ML) was chosen from the experimental results [27] as having better carbon gasification and suitable activation energies for the CH4 dehydrogenations, compared to high Cu loading (1.0 ML) and low Cu loading (1/9ML). The slab model surface is shown in Figure 1. Unless otherwise specified, the bottom-most layer in the slab models was constrained at the bulk position, and the top three layers, as well as the adsorbates, were allowed to relax during geometry optimization and transition state (TS) search [55]. The TS and minimum energy paths are obtained using the climbing-image nudged elastic band (Climbing-NEB) method, and confirmed by a single imaginary frequency from the vibrational analysis [56].

4. Conclusions

In this study, density functional theory (DFT) calculations were used to understand the mechanism of DRM over a Ni2Cu (111) catalyst. In addition, carbon deposition and elimination for have been investigated at the molecular level and under realistic reaction conditions. All elementary reaction steps, including CH4 dehydrogenation, CO2 activation, CH and C oxidations, were comprehensively described. The results showed a significant improvement in inhibiting carbon deposition compared to similar reactions on the Ni (111) and Ni2Fe (111) surfaces. First, the energy barriers in CH4 dissociation steps were found to be relatively increased to Ni (111) and Ni2Fe (111) catalysts. Specifically, the CH dissociation was found to be the rate-limiting step with the highest energy barrier of 1.41 eV when compared to 1.33 eV for Ni (111) and 1.36 eV for Ni2Fe (111) catalysts. Moreover, the source of coke had been confirmed to be from CH4 dissociation rather than CO2, which agrees with previous studies on Ni-based catalysts. The coke elimination through oxidation is remarkably improved by both OH and atomic O species as competitive oxidants. From the C/CH oxidation pathways, C* + O* → CO (g) was found to be the most favorable reaction with an energy barrier of 0.72 eV. This value is less than half of similar oxidation reaction on Ni (111) and Ni2Fe (111) catalysts. In addition, coke deposition is physically prevented by H2O (g), which is easily formed in a Ni2Cu (111) system. Combining the above results showed that the Ni-Ni surface with Cu atoms could improve carbon elimination (by C* + O* reaction), while simultaneously reducing carbon deposition from CH dissociation.
The vibrational frequencies were used to evaluate the rate constant k of coke deposition and elimination reactions at typical DRM reaction temperatures. The rate constant of CH dissociation was lower than those of carbon oxidation (kC(O)and kC(OH)) by atomic O or OH species, respectively. Their calculated values revealed the importance of OH species in bypassing the coke formation step from CH dissociation. Last, when carbon deposition occurs on the Ni2Cu (111) surface, coke could be instantaneously removed by various oxidation reactions. In conclusion, this efficient periodic cycle of carbon deposition and carbon elimination results in a more stable catalytic performance and excellent carbon deposition resistance in DRM.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4344/10/9/1043/s1.

Author Contributions

A.O. “investigation; writing original draft”; S.H.Y. “formal analysis; writing review and editing”; M.K. “review and editing”; M.G. “conceptualization; A.C. “writing review and editing”; N.E. “conceptualization, supervision, and project administration”. All authors have read and agreed to the published version of the manuscript.

Funding

This paper was made possible by an NPRP exceptional grant award [NPRP-EP X-100-2-024] from the Qatar National Research Fund (a member of the Qatar Foundation). The statements made herein are solely the responsibility of the authors.

Acknowledgments

All the simulations were performed on the supercomputing clusters RAAD & RAAD2. The hardware & software support from Texas Advanced Scientific Computing Center, Texas A&M University at Qatar, is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Uykun Mangaloğlu, D.; Baranak, M.; Ataç, Ö.; Atakül, H. Effect of the promoter presence in catalysts on the compositions of Fischer–Tropsch synthesis products. J. Ind. Eng. Chem. 2018, 66, 298–310. [Google Scholar] [CrossRef]
  2. Arora, S.; Prasad, R. An overview on dry reforming of methane: Strategies to reduce carbonaceous deactivation of catalysts. RSC Adv. 2016, 6, 108668–108688. [Google Scholar] [CrossRef]
  3. Challiwala, M.S.; Wilhite, B.A.; Ghouri, M.M.; Elbashir, N.O. Multidimensional modeling of a microfibrous entrapped cobalt catalyst Fischer-Tropsch reactor bed. AIChE J. 2018, 64, 1723–1731. [Google Scholar] [CrossRef]
  4. Gu, B.; Khodakov, A.Y.; Ordomsky, V.V. Selectivity shift from paraffins to α-olefins in low temperature Fischer–Tropsch synthesis in the presence of carboxylic acids. Chem. Commun. 2018, 54, 2345–2348. [Google Scholar] [CrossRef] [PubMed]
  5. Singh, R.; Dhir, A.; Mohapatra, S.K.; Mahla, S.K. Dry reforming of methane using various catalysts in the process: Review. Biomass Convers. Biorefin. 2020, 10, 567–587. [Google Scholar] [CrossRef]
  6. Challiwala, M.S.; Afzal, S.; Choudhury, H.A.; Sengupa, D.; El-Halwagi, M.M.; Elbashir, N.O. Alternative via Pathways Dry Reforming for CO of 2 Utilization Methane. In Advances in Carbon Management Technologies Volume 1: Carbon Removal, Renewable and Nuclear Energy; CRC Press: Boca Raton, FL, USA, 2020; p. 253. Available online: https://www.taylorfrancis.com/books/9780429243608/chapters/10.1201/9780429243608-14 (accessed on 10 September 2020).
  7. Wei, J.; Iglesia, E. Isotopic and kinetic assessment of the mechanism of reactions of CH4 with CO2 or H2O to form synthesis gas and carbon on nickel catalysts. J. Catal. 2004, 224, 370–383. [Google Scholar] [CrossRef]
  8. Niu, J.; Ran, J.; Chen, D. Understanding the mechanism of CO2 reforming of methane to syngas on Ni@Pt surface compared with Ni(111) and Pt(111). Appl. Surf. Sci. 2020, 513, 145840. [Google Scholar] [CrossRef]
  9. Takami, D.; Yamamoto, A.; Yoshida, H. Dry reforming of methane over alumina-supported rhodium catalysts at low temperatures under visible and near-infrared light. Catal. Sci. Technol. 2020, 10, 5811–5814. [Google Scholar] [CrossRef]
  10. Zhu, Y.; Chen, K.; Barat, R.; Mitra, S. Dry reforming of methane over ruthenium/carbon nanotube catalyst. Chem. Eng. 2020, 4, 16. [Google Scholar] [CrossRef] [Green Version]
  11. Maina, S.C.P.; Ballarini, A.D.; Vilella, J.I.; de Miguel, S.R. Study of the performance and stability in the dry reforming of methane of doped alumina supported iridium catalysts. Catal. Today 2020, 344, 129–142. [Google Scholar] [CrossRef]
  12. Leba, A.; Yıldırım, R. Determining most effective structural form of nickel-cobalt catalysts for dry reforming of methane. Int. J. Hydro. Energy 2020, 45, 4268–4283. [Google Scholar] [CrossRef]
  13. Pakhare, D.; Spivey, J. A review of dry (CO2) reforming of methane over noble metal catalysts. Chem. Soc. Rev. 2014, 43, 7813–7837. [Google Scholar] [CrossRef]
  14. Kawi, S.; Kathiraser, Y.; Ni, J.; Oemar, U.; Li, Z.W.; Saw, E.T. Progress in synthesis of highly active and stable nickel-based catalysts for carbon dioxide reforming of methane. Chemsuschem. 2015, 8, 3556–3575. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, Y.; Yao, L.; Wang, S.H.; Mao, D.H.; Hu, C.W. Low-temperature catalytic CO2 dry reforming of methane on Ni-based catalysts: A review. Fuel Process. Technol. 2018, 169, 199–206. [Google Scholar] [CrossRef]
  16. Ren, J.; Qin, X.; Yang, J.Z.; Qin, Z.F.; Guo, H.L.; Lin, J.Y.; Li, Z. Methanation of carbon dioxide over Ni-M/ZrO2 (M = Fe, Co, Cu) catalysts: Effect of addition of a second metal. Fuel Process. Technol. 2015, 137, 204–211. [Google Scholar] [CrossRef]
  17. Bian, Z.F.; Das, S.; Wai, M.H.; Hongmanorom, P.; Kawi, S. A review on bimetallic nickel-based catalysts for CO2 reforming of methane. Chemphyschem 2017, 18, 3117–3134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Wang, M.Y.; Fu, Z.M.; Yang, Z.X. The carbon-tolerance mechanism of Ni-based alloy with coinage metals. Phys. Lett. A 2013, 377, 2189–2194. [Google Scholar] [CrossRef]
  19. Xu, L.-L.; Wen, H.; Jin, X.; Bing, Q.-M.; Liu, J.-Y. DFT study on dry reforming of methane over Ni2Fe overlayer of Ni (1 1 1) surface. Appl. Surf. Sci. 2018, 443, 515–524. [Google Scholar] [CrossRef]
  20. Theofanidis, S.A.; Batchu, R.; Galvita, V.V.; Poelman, H.; Marin, G.B. Carbon gasification from Fe–Ni catalysts after methane dry reforming. Appl. Catal. B Environ. 2016, 185, 42–55. [Google Scholar] [CrossRef]
  21. Theofanidis, S.A.; Galvita, V.V.; Poelman, H.; Marin, G.B. Enhanced carbon-resistant dry reforming Fe–Ni catalyst: Role of Fe. ACS Catal. 2015, 5, 3028–3039. [Google Scholar] [CrossRef]
  22. Lo Faro, M.; Frontera, P.; Antonucci, P.; Arico, A.S. Ni-Cu based catalysts prepared by two different methods and their catalytic activity toward the ATR of methane. Chem. Eng. Res. Des. 2015, 93, 269–277. [Google Scholar] [CrossRef]
  23. Rahemi, N.; Haghighi, M.; Babaluo, A.A.; Allahyari, S.; Jafari, M.F. Syngas production from reforming of greenhouse gases CH4/CO2 over Ni–Cu/Al2O3 nanocatalyst: Impregnated vs. plasma-treated catalyst. Energ. Convers. Manag. 2014, 84, 50–59. [Google Scholar] [CrossRef]
  24. Wu, T.; Cai, W.Y.; Zhang, P.; Song, X.F.; Gao, L. Cu-Ni@ SiO2 alloy nanocomposites for methane dry reforming catalysis. RSC Adv. 2013, 3, 23976–23979. [Google Scholar] [CrossRef]
  25. Liu, H.Y.; Zhang, R.G.; Yan, R.X.; Li, J.R.; Wang, B.J.; Xie, K.C. Insight into CH4 dissociation on NiCu catalyst: A first-principles study. Appl. Surf. Sci. 2012, 258, 8177–8184. [Google Scholar] [CrossRef]
  26. An, W.; Zeng, X.C.; Turner, C.H. First-principles study of methane dehydrogenation on a bimetallic Cu/Ni (111) surface. J. Chem. Phys. 2009, 131, 174702. [Google Scholar] [CrossRef] [PubMed]
  27. Chatla, A.; Ghouri, M.M.; El Hassan, O.W.; Mohamed, N.; Prakash, A.V.; Elbashir, N.O. An experimental and first principles DFT investigation on the effect of Cu addition to Ni/Al2O3 catalyst for the dry reforming of methane. Appl. Catal. A-Gen. 2020, 602, 117699. [Google Scholar] [CrossRef]
  28. Kresse, G.; Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15–50. [Google Scholar] [CrossRef]
  29. Wang, Z.Y.; Cao, X.M.; Zhu, J.H.; Hu, P. Activity and coke formation of nickel and nickel carbide in dry reforming: A deactivation scheme from density functional theory. J. Catal. 2014, 311, 469–480. [Google Scholar] [CrossRef]
  30. Arevalo, R.L.; Aspera, S.M.; Escano, M.C.S.; Nakanishi, H.; Kasai, H. Tuning methane decomposition on stepped Ni surface: The role of subsurface atoms in catalyst design. Sci. Rep. 2017, 7, 13963. [Google Scholar] [CrossRef] [Green Version]
  31. Heiland, W. Charge-exchange processes and surface-chemistry. Surf. Sci. 1991, 251, 942–946. [Google Scholar] [CrossRef]
  32. Jang, W.J.; Shim, J.O.; Kim, H.M.; Yoo, S.Y.; Roh, H.S. A review on dry reforming of methane in aspect of catalytic properties. Catal. Today 2019, 324, 15–26. [Google Scholar] [CrossRef]
  33. Kim, S.M.; Abdala, P.M.; Margossian, T.; Hosseini, D.; Foppa, L.; Armutlulu, A.; van Beek, W.; Comas-Vives, A.; Copéret, C.; Müller, C. Cooperativity and dynamics increase the performance of NiFe dry reforming catalysts. J. Am. Chem. Soc. 2017, 139, 1937–1949. [Google Scholar] [CrossRef] [PubMed]
  34. Fan, X.; Liu, Z.; Zhu, Y.-A.; Tong, G.; Zhang, J.; Engelbrekt, C.; Ulstrup, J.; Zhu, K.; Zhou, X. Tuning the composition of metastable CoxNiyMg100−x−y(OH)(OCH3) nanoplates for optimizing robust methane dry reforming catalyst. J. Catal. 2015, 330, 106–119. [Google Scholar] [CrossRef]
  35. Michaelides, A.; Hu, P. Methyl chemisorption on Ni (111) and C-H-M multicentre bonding: A density functional theory study. Surf. Sci. 1999, 437, 362–376. [Google Scholar] [CrossRef]
  36. Zhu, Y.-A.; Dai, Y.-C.; Chen, D.; Yuan, W.-K. First-principles calculations of CH4 dissociation on Ni (100) surface along different reaction pathways. J. Mol. Catal. A Chem. 2007, 264, 299–308. [Google Scholar] [CrossRef]
  37. Bradford, M.C.J.; Vannice, M.A. Catalytic reforming of methane with carbon dioxide over nickel catalysts II. Reaction kinetics. Appl. Catal. A-Gen. 1996, 142, 97–122. [Google Scholar] [CrossRef]
  38. Maslov, M.M.; Openov, L.A.; Podlivaev, A.I. On the vineyard formula for the pre-exponential factor in the Arrhenius law. Phys. Solid State 2014, 56, 1239–1244. [Google Scholar] [CrossRef] [Green Version]
  39. Monkhorst, H.J.; Pack, J.D. Special points for brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  40. Jørgensen, M. Microkinetic Modeling of Nanoparticl Catalysis Using Density Functional Theory. Ph.D. Thesis, Chalmers University of Technology, Gothenburg, Sweden, 2017. [Google Scholar]
  41. Jørgensen, M.; Grönbeck, H. First-principles microkinetic modeling of methane oxidation over Pd (100) and Pd (111). ACS Catal. 2016, 6, 6730–6738. [Google Scholar] [CrossRef]
  42. Ray, D.; Ghosh, S.; Tiwari, A.K. Controlling heterogeneous catalysis of water dissociation using Cu-Ni bimetallic alloy surfaces: A quantum dynamics study. J. Phys. Chem. A 2018, 122, 5698–5709. [Google Scholar] [CrossRef]
  43. Zhang, R.G.; Guo, X.Q.; Wang, B.J.; Ling, L.M. Insight Into the effect of CuNi (111) and FeNi (111) surface structure and second metal composition on surface carbon elimination by O or OH: A comparison study with Ni (111) surface. J. Phys. Chem. C 2015, 119, 14135–14144. [Google Scholar] [CrossRef]
  44. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef] [PubMed]
  45. Blochl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [Green Version]
  47. Hammer, B.; Hansen, L.B.; Norskov, J.K. Improved adsorption energetics within density-functional theory using revised Perdew-Burke-Ernzerhof functionals. Phys. Rev. B 1999, 59, 7413–7421. [Google Scholar] [CrossRef] [Green Version]
  48. Hammer, B.; Norskov, J.K. Theoretical surface science and catalysis—Calculations and concepts. Adv. Catal. 2000, 45, 71–129. [Google Scholar] [CrossRef]
  49. Methfessel, M.; Paxton, A.T. High-precision sampling for brillouin-zone integration in metals. Phys. Rev. B 1989, 40, 3616–3621. [Google Scholar] [CrossRef] [Green Version]
  50. Zhu, Y.A.; Chen, D.; Zhou, X.G.; Yuan, W.K. DFT studies of dry reforming of methane on Ni catalyst. Catal. Today 2009, 148, 260–267. [Google Scholar] [CrossRef]
  51. Henkelman, G.; Jonsson, H. Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points. J. Chem. Phys. 2000, 113, 9978–9985. [Google Scholar] [CrossRef] [Green Version]
  52. Yoon, S.H.; Kang, U.; Park, H.; Abdel-Wahab, A.; Han, D.S. Computational density functional theory study on the selective conversion of CO2 to formate on homogeneously and heterogeneously mixed CuFeO2 and CuO surfaces. Catal. Today 2019, 335, 345–353. [Google Scholar] [CrossRef]
  53. Henkelman, G.; Jonsson, H. A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives. J. Chem. Phys. 1999, 111, 7010–7022. [Google Scholar] [CrossRef]
  54. Yoon, S.H.; Yu, C.; Han, A.; Park, H.; Elbashir, N.; Han, D.S. Computational characterization of nitrogen-doped carbon nanotube functionalized by Fe adatom and Fe substituent for oxygen reduction reaction. Appl. Surf. Sci. 2019, 485, 342–352. [Google Scholar] [CrossRef]
  55. Henkelman, G.; Arnaldsson, A.; Jónsson, H. Theoretical calculations of CH4 and H2 associative desorption from Ni (111): Could subsurface hydrogen play an important role? J. Chem. Phys. 2006, 124, 044706. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Fan, C.; Zhu, Y.A.; Yang, M.L.; Sui, Z.J.; Zhou, X.G.; Chen, D. Density functional theory-assisted microkinetic analysis of methane dry reforming on Ni catalyst. Ind. Eng. Chem. Res. 2015, 54, 5901–5913. [Google Scholar] [CrossRef]
Catalysts 10 01043 g001 550
Figure 1. Schematic illustration of Ni2Cu overlayer of Ni (111) surface. FCC, HCP, TopNi, TopCu, B2Ni, and BNiCu represent a face centered cubic, hexagonal close packed, Ni at the top site, Cu at the top site, Ni–Ni at the bridge site and Ni–Cu at the bridge site, respectively. Green and brown spheres represent Ni and Cu atoms, respectively.
Figure 1. Schematic illustration of Ni2Cu overlayer of Ni (111) surface. FCC, HCP, TopNi, TopCu, B2Ni, and BNiCu represent a face centered cubic, hexagonal close packed, Ni at the top site, Cu at the top site, Ni–Ni at the bridge site and Ni–Cu at the bridge site, respectively. Green and brown spheres represent Ni and Cu atoms, respectively.
Catalysts 10 01043 g001
Catalysts 10 01043 g002 550
Figure 2. Activation energy profile for successive CH4 dissociation steps. Green, brown, blue, and grey spheres represent Ni, Cu, H, and C atoms, respectively.
Figure 2. Activation energy profile for successive CH4 dissociation steps. Green, brown, blue, and grey spheres represent Ni, Cu, H, and C atoms, respectively.
Catalysts 10 01043 g002
Catalysts 10 01043 g003 550
Figure 3. Activation energy profile for possible CH3 reaction pathways. Green, brown, blue, grey, and red spheres represent Ni, Cu, H, C, and O atoms, respectively.
Figure 3. Activation energy profile for possible CH3 reaction pathways. Green, brown, blue, grey, and red spheres represent Ni, Cu, H, C, and O atoms, respectively.
Catalysts 10 01043 g003
Catalysts 10 01043 g004 550
Figure 4. Activation energy profile for H-induced CO2 dissociation (hydrogenation).
Figure 4. Activation energy profile for H-induced CO2 dissociation (hydrogenation).
Catalysts 10 01043 g004
Catalysts 10 01043 g005 550
Figure 5. Comparison of activation energy profiles of C* + O* oxidation reaction for Ni2Cu (111) (red) and Ni2Fe (111) (green) surfaces.
Figure 5. Comparison of activation energy profiles of C* + O* oxidation reaction for Ni2Cu (111) (red) and Ni2Fe (111) (green) surfaces.
Catalysts 10 01043 g005
Catalysts 10 01043 g006 550
Figure 6. Activation energy profile of CH oxidation by OH. The OH is an effective oxidant that helps to skip CH dissociation. The transition state (TS) configuration shows where OH is adsorbed on the TopCu site while CH is adsorbed on the TopNi site.
Figure 6. Activation energy profile of CH oxidation by OH. The OH is an effective oxidant that helps to skip CH dissociation. The transition state (TS) configuration shows where OH is adsorbed on the TopCu site while CH is adsorbed on the TopNi site.
Catalysts 10 01043 g006
Catalysts 10 01043 g007 550
Figure 7. Activation energy profiles for two different CHOH* dissociation pathways.
Figure 7. Activation energy profiles for two different CHOH* dissociation pathways.
Catalysts 10 01043 g007
Catalysts 10 01043 g008 550
Figure 8. The energy of H2O formation via OH* + H* reaction and configuration of transition state. The weakly adsorbed H2O physically removes coke from the catalyst surface at high temperatures (steam).
Figure 8. The energy of H2O formation via OH* + H* reaction and configuration of transition state. The weakly adsorbed H2O physically removes coke from the catalyst surface at high temperatures (steam).
Catalysts 10 01043 g008
Catalysts 10 01043 g009 550
Figure 9. Rate constants of the coke removal reaction at different temperatures.
Figure 9. Rate constants of the coke removal reaction at different temperatures.
Catalysts 10 01043 g009
Catalysts 10 01043 g010 550
Figure 10. A schematic diagram of the DRM dominant reaction pathway on Ni2Cu (111).
Figure 10. A schematic diagram of the DRM dominant reaction pathway on Ni2Cu (111).
Catalysts 10 01043 g010
Table
Table 1. Calculated adsorption energies of all adsorbed species involved in DRM on Ni2Cu overlayer of Ni (111) surface.
Table 1. Calculated adsorption energies of all adsorbed species involved in DRM on Ni2Cu overlayer of Ni (111) surface.
SpeciesEads (eV)Favored Adsorption Site
1CH4*−0.02N/A
2CH3*−2.17FCC
3CH2*−4.99FCC
4CH*−5.43FCC
5C*−7.04HCP
6O*−6.34FCC
7OH*−3.28B2Ni
8CO*−1.52B2Ni
9COH*−4.07B2Ni
10CHOH*−2.40B2Ni
11CO2*−0.03N/A
12H*−3.60FCC
13H2*−0.06N/A
14CH3O*−2.43B2Ni
15CH2O−0.04N/A
16CH2OH*−1.87TopNi
17CH3OH*−0.05N/A
18CHO*a−2.05TopNi
19COH*a−4.07B2Ni
20H2O*−0.06N/A
21COOH*−2.28TopNi
a CHO is “H–C=O*” and COH* is “C–O–H*”.
Table
Table 2. Calculated activation energies of forward (Ea,f), reverse (Ea,r), and enthalpy (Hf) in an elementary reaction of dry reforming of methane (DRM).
Table 2. Calculated activation energies of forward (Ea,f), reverse (Ea,r), and enthalpy (Hf) in an elementary reaction of dry reforming of methane (DRM).
ReactionEa,f+ (eV)Hf+ (eV)Ea,r+ (eV)
CH4* → CH3* + H*1.300.780.52
CH3* → CH2* + H*0.750.330.42
CH2* → CH* + H*0.49−0.060.55
CH* → C* + H*1.400.690.71
C* + O* → CO*0.72−2.172.89
CH* + O* → CHO*1.06−0.301.36
CHO* → CO* + H*0.18−1.171.35
C* + OH* → COH*1.13−0.952.07
H* + H* → H20.730.570.16
CH* + OH* → CHOH*1.080.240.84
CHOH* → CHO* + H*0.60−0.310.91
CO2* + H* → COOH*0.93−0.121.05
COOH* → CO* + OH*0.43−0.591.02
O* + H* → OH*0.94−0.241.18
H* + OH* → H2O*0.91−0.151.06
CO2* → CO* + O*1.690.481.21
CH3* + O* → CH3O*1.37−0.221.59
CH2* + O* → CH2O*0.76−0.431.19
CHOH* → COH* + H*0.20−0.530.73
COH* → CO* + H*1.01−1.092.10
CH3* + OH* → CH3OH*1.60−0.141.74
CH2* + OH* → CH2OH*0.78−0.100.88
CH3O* → CH2O* + H*0.960.570.39
CH2OH* → CH2O* + H*0.75−0.070.82
CH2OH* → CHOH* + H*0.920.420.50
CH3OH* → CH3O* + H*4.290.014.28
CH3OH* → CH2OH* + H*2.461.391.10
+ Zero–point energy (ZPE) corrections are not included.
Table
Table 3. The rate constant k (s−1) of carbon elimination and deposition on Ni2Cu (111) surface at different temperature ranges.
Table 3. The rate constant k (s−1) of carbon elimination and deposition on Ni2Cu (111) surface at different temperature ranges.
Reactions800 K850 K900 K950 K1000 K
CO* → C* + O*4.36 × 10−75.11 × 10−64.56 × 10−53.23 × 10−41.88 × 10−3
CH* → C* + H*1.18 × 1053.90 × 1051.13 × 1062.92 × 1066.88 × 106
C* + O* → CO*2.43 × 1084.50 × 1087.74 × 1081.25 × 1091.95 × 109
C* + OH* → COH*5.75 × 1061.50 × 1073.53 × 1077.59 × 1071.51 × 108
CH* + O* → CHO*6.01 × 1061.49 × 1073.32 × 1076.83 × 1081.31 × 108
CH* + OH* → CHOH*3.28 × 1078.27 × 1071.88 × 1083.93 × 1087.62 × 108
kCH/kC(O)0.020.030.040.050.07
kCH/kC(OH)5.705.505.335.185.05

Share and Cite

MDPI and ACS Style

Omran, A.; Yoon, S.H.; Khan, M.; Ghouri, M.; Chatla, A.; Elbashir, N. Mechanistic Insights for Dry Reforming of Methane on Cu/Ni Bimetallic Catalysts: DFT-Assisted Microkinetic Analysis for Coke Resistance. Catalysts 2020, 10, 1043. https://0-doi-org.brum.beds.ac.uk/10.3390/catal10091043

AMA Style

Omran A, Yoon SH, Khan M, Ghouri M, Chatla A, Elbashir N. Mechanistic Insights for Dry Reforming of Methane on Cu/Ni Bimetallic Catalysts: DFT-Assisted Microkinetic Analysis for Coke Resistance. Catalysts. 2020; 10(9):1043. https://0-doi-org.brum.beds.ac.uk/10.3390/catal10091043

Chicago/Turabian Style

Omran, Ahmed, Sun Hee Yoon, Murtaza Khan, Minhaj Ghouri, Anjaneyulu Chatla, and Nimir Elbashir. 2020. "Mechanistic Insights for Dry Reforming of Methane on Cu/Ni Bimetallic Catalysts: DFT-Assisted Microkinetic Analysis for Coke Resistance" Catalysts 10, no. 9: 1043. https://0-doi-org.brum.beds.ac.uk/10.3390/catal10091043

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop