Next Article in Journal
Valorization of Glycerol through the Enzymatic Synthesis of Acylglycerides with High Nutritional Value
Next Article in Special Issue
Dry Reforming of Methane Using Ce-modified Ni Supported on 8%PO4 + ZrO2 Catalysts
Previous Article in Journal
The Composition and Structure of Ultra-Dispersed Mixed Oxide (II, III) Particles and Their Influence on In-Situ Conversion of Heavy Oil
Previous Article in Special Issue
Highly-Dispersed Ni-NiO Nanoparticles Anchored on an SiO2 Support for an Enhanced CO Methanation Performance
 
 
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 1
10.3390/catal10010115
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

Selective Hydrogenation of Acetylene Catalysed by a B12N12 Cluster Doped with a Single Nickel Atom: A DFT Study

by
Yun Wang
and
Lihua Kang
*
School of Chemistry & Chemical Engineering, Shihezi University, Shihezi 832003, Xinjiang, China
*
Author to whom correspondence should be addressed.
Catalysts 2020, 10(1), 115; https://doi.org/10.3390/catal10010115
Submission received: 9 December 2019 / Revised: 2 January 2020 / Accepted: 9 January 2020 / Published: 13 January 2020
(This article belongs to the Special Issue Recent Advances in Nickel-Based Catalysts)
Browse Figures
Versions Notes

Abstract

:
To obtain a catalyst based on a non-precious metal that can replace traditional palladium-based selective catalysts of acetylene hydrogenation, the catalytic performances of two different configurations of a B12N12 cluster doped with a single nickel atom were studied by a density functional theory computational approach. After analysing the effect that the adsorption of reactants onto the clusters has on the reaction path, we determined the lowest energy path for the acetylene double hydrogenation. Comparing the acetylene hydrogenation activities and ethylene product selectivities of the B11N12Ni and B12N11Ni clusters, which have different doping sites, we determined the activities of these two catalysts to be similar to each other; however, the B11N12Ni cluster was calculated to have higher selectivity for ethylene as a product. This difference may be related to the moderate adsorption of hydrogen and acetylene on the B11N12Ni cluster. As a new type of nickel-based single-atom catalyst, B11N12Ni clusters may have research value in the selective hydrogenation of acetylene.

Graphical Abstract

1. Introduction

Currently, the selective hydrogenation of acetylene to produce ethylene is still mainly realised through use of catalysts comprising palladium or other precious metals as the active material [1,2,3,4], in the presence of a carrier and a small amount of inhibitor such as Zn [5], Ag [6], Cu [7] and Si [8]. Palladium-based catalysts are favoured by researchers as a result of their high acetylene hydrogenation activity and their selectivity for ethylene as product [2,9]. On the other hand, palladium is an expensive heavy metal that is quite harmful to the environment [10,11]. Therefore, it is imperative to improve the efficiency and accessibility of the described reaction by finding selective acetylene hydrogenation catalysts of increased activity and decreased cost and environmental impact.
In 2011, Zhang et al. were the first to propose the concept of single-atom catalyst (SAC) [12], which has attracted much research attention, because SACs greatly improve the utilisation rate of the catalysts’ active components. SACs have special geometric and electronic characteristics through the coordination of the active centre of a single metal atom and the carrier material, and thus show a completely different catalytic ability than alloys and the pure metal [13]. A series of palladium-based acetylene selective hydrogenation SACs have been introduced, leading researchers to achieve a significant improvement in the catalytic performance [14,15,16,17]. The transition metal nickel can also catalyse the hydrogenation of acetylene to ethylene. However, pure metallic nickel displays poor catalytic selectivity, so, in the presence of nickel-based catalysts, the described hydrogenation reaction needs to be conducted at high temperature [18,19,20]. We thus tried to design a nickel-based selective acetylene hydrogenation SAC of improved efficacy with respect to the known ones. In such SAC, nickel is the active centre, and Ni atoms at high dispersion are used to dope the catalyst support surface, which is characterised by high surface area, based on the assumption that such a feature may improve the catalyst’s activity and selectivity [21].
Since the successful synthesis of the smallest boron nitride fullerene B12N12 cluster in 1998 [22], numerous studies have been conducted focusing on this material. Strout’s research proved the structure of the B12N12 cluster to be very stable [23]. Oku et al. performed molecular orbital calculations to demonstrate the possibility of hydrogen storage in boron nitride clusters [24,25,26]. Kuo et al. modified the B12N12 cluster by incorporating a carbon atom within it to improve this boron nitride fullerene’s hydrogen storage performance [27]. The results of Ayub et al.’s research indicate that hydrogen adsorption on the surface of Ni-modified B12N12 (Ni-BN) nanoclusters significantly increased with respect to their unmodified counterparts [28]. Therefore, we turned our attention to boron nitride fullerene B12N12. Previously, our group has reported that a Pd-doped BN nanocluster can be used as a selective catalyst for acetylene hydrogenation [17].
In the present study, a B12N12 cluster is combined with the transition metal nickel as part of the effort to develop an effective and selective catalyst for the hydrogenation of acetylene. Notably, the coordination of B or N atoms has been shown to enhance the ethylene product selectivity of metallic nickel catalysts, and the BN fullerene carrier has been observed to help improve the utilisation of metal atoms to 100%, while preventing excessive adsorption of hydrogen and acetylene onto the catalyst, thus ensuring the catalyst’s selectivity. This study has certain reference value in the search for non-precious metals catalysts that may replace traditional palladium-based catalysts of the selective hydrogenation of acetylene to ethylene.

2. Results and Discussion

2.1. Structure Optimisation of the Reactants

In Figure 1 are reported the results of the optimisation calculations performed on all reactants and catalysts. The “undoped” B12N12 fullerene is characterised by a Th-symmetrical structure with both B and N atoms. The value of Mulliken charge for the B atom is 0.57 and for the N atom is −0.57. As can be evinced from Figure 1a and Table S1, the B12N12 cluster consists of four- and six-membered rings. Two kinds of B-N bonds are present: the length of the B-N bond within the four-membered ring is 1.487 Å, whereas that within the six-membered ring is 1.439 Å. Following the introduction of a Ni atom in the B12N12 cluster, two different configurations are obtained, depending on which element is replaced by Ni: the B11N12Ni cluster (Figure 1b) and the B12N11Ni cluster (Figure 1c).
In the B11N12Ni cluster, the lengths of the bonds linking the metal centre to the two adjacent nitrogen atoms N(1) and N(2) are 1.893 Å and 1.927 Å, respectively, with the relevant bond angle measuring 108.7°. By comparison with the undoped B12N12 cluster, Ni doping results in bond length increases and bond angle decreases; furthermore, following the doping process, the Ni atom is convex significantly outward as shown in Figure 1b and Table S2. In the B11N12Ni cluster, the value of the Mulliken charge is −0.43 for N(1) and −0.34 for N(2) as a result of a charge transfer from the N atoms to the Ni centre.
In the B12N11Ni cluster, the lengths of the bonds linking the metal centre to the two adjacent boron atoms B(1) and B(2) are 1.960 Å and 1.977 Å, respectively; these bonds are thus longer than the Ni-N bonds in B11N12Ni, and the relevant bond angle measures 92.01°. The Ni atom is more convex outward as shown in Figure 1c and Table S3. In the B12N11Ni cluster, the value of the Mulliken charge is 1.26 for B(1) and 1.12 for B(2), as a result of a charge transfer from the B atoms to the Ni atom.
To render more intuitive the charge distribution in the BN fullerene before and after its doping with Ni, we analysed the molecular electrostatic potential of the B12N12, B11N12Ni and B12N11Ni clusters on the molecular van der Waals surface (electron density: 0.001). In the structures reported in Figure 2, the more intense the blue colour, the more positive the charge, and the more intense the red colour, the more negative the charge. As can be evinced from the mentioned Figure 2a, in undoped BN fullerene, the B atoms are positively charged and the N atoms are negatively charged. In Figure 2b, the electrostatic potential minimum is −0.95 eV above N(2) and the maximum is 2.23 eV above the Ni atom. In Figure 2c, the electrostatic potential minimum is −0.79 eV above N(1) and the maximum is 2.08 eV above the Ni atom. In the B11N12Ni and B12N11Ni clusters, the Ni atom is highest electropositive, so it has the potential to act as a nucleophilic reaction site, and the N atoms around the metal centre exhibit stronger electronegativity than their counterparts in undoped BN fullerene, which may in turn favour electrophilic reactions.
To deepen understanding of the characteristic of the BN fullerene before and after its doping with Ni, we performed the frontier molecular orbital analysis on the B12N12, B11N12Ni and B12N11Ni cluster. The energies of the frontier orbitals are listed in Table 1. The value of the LUMO − HOMO energy gap of the B12N12 cluster is 6.72 eV, which is consistent with that reported in the literature [29]. Furthermore, the energy of the HOMO of each of the two BN clusters doped with the nickel atom is higher than that of the HOMO of the undoped cluster and the energy of the LUMO of each of the two BN clusters doped with the nickel atom is lower than that of the LUMO of the undoped cluster. The lower the energy of LUMO, the more conducive to the filling of electrons, and the higher the energy level of HOMO, the more conducive to the loss of electrons. Therefore, the energy gap between HOMO and LUMO is reduced as a result of the doping process, the change that favours in the next adsorption of the reactants, thus promoting the catalytic reaction activity.
To better understand the molecular orbital energy changes associated with the doping of the B12N12 cluster with a Ni atom and the contribution of each component to the total density of states, we plotted the total density of states and the density of partial states of the undoped B12N12, B11N12Ni and B12N11Ni clusters. These plots are reported in Figure 3, where each cluster’s frontier orbitals are also indicated. As can be evinced from Figure 3a, the value of LUMO − HOMO gap is large, thus the B12N12 cluster displays the features of an insulator. The p orbitals of the N atoms are the main contributors of the HOMO and of the lower energy occupied molecular orbitals; on the other hand, the p orbitals of the B atoms are the main contributors of the LUMO and of the higher energy unoccupied molecular orbitals. The features of the plot reported in Figure 3b suggest that, after replacing a B atom with a Ni atom, a new unoccupied molecular orbital is added in the original forbidden band, and this new orbital becomes the LUMO of the B11N12Ni cluster. Notably, the original HOMO increases in energy and shifts to the right. Therefore, upon doping the B12N12 cluster with a nickel atom to form the B11N12Ni cluster, the LUMO − HOMO energy gap decreases. In the B11N12Ni cluster; moreover, the p orbitals of the N atoms are the main contributors to the HOMO, LUMO and the lower energy occupied molecular orbitals, with the contribution to these orbitals of the Ni atom being negligible. The HOMO and LUMO population maps of the B11N12Ni cluster are also reported in Figure 3b; in this case, the electrons are concentrated from N atoms to around the Ni atom compared to the undoped B12N12 cluster (see Figure 3a). As can be evinced from Figure 3c, after replacing an N atom with a Ni atom, a new occupied molecular orbital is added to the original forbidden band, and this new orbital becomes the HOMO of the B12N11Ni cluster. Notably, the energy of the original LUMO decreases and moves to the left. As a result of these changes, the LUMO − HOMO energy gap of the B12N11Ni cluster is smaller than that of the B12N12 cluster. Additionally, in the B12N11Ni cluster, the p orbitals of the B atoms are the main contributors of the HOMO and of the higher energy unoccupied molecular orbitals; on the other hand, the p orbitals of the N atoms are the main contributors of the lower energy occupied molecular orbitals. The s and p orbitals of the B atoms are major contributors of the LUMO, and the orbitals of the Ni atom make some contribution to the HOMO and LUMO. Based on the HOMO and LUMO population maps of the B12N11Ni cluster, we can infer that electrons accumulate near the Ni atom and the B atoms in the metal centre’s immediate vicinity. We speculate that Ni and its surrounding B atoms have strong activity.

2.2. Adsorption of C2H2 and H2 onto BN Fullerene or the BN-Ni Cluster

Implementing our computational approach, we obtained the stable adsorption configuration and the value for the adsorption energy of acetylene and hydrogen onto BN fullerene and nickel-doped BN clusters. The stable adsorption configurations are reported in Figure 4, whereas the values for the adsorption energies are reported in Table 2.
In Figure 4a and Table S4 are depicted the optimised configuration of acetylene adsorbed on the B12N12 cluster. According to our data, almost no differences exist between the structures of acetylene and the B12N12 clusters before and after adsorption. Acetylene and the cluster are quite far from each other, and the adsorption of acetylene on the cluster is physical in nature, with the value of the adsorption energy calculated to be −2.83 kcal/mol.
In Figure 4b and Table S5 are depicted the stable configuration of acetylene adsorbed on the B11N12Ni cluster. In this case, some changes in the configuration of acetylene and the B11N12Ni cluster are observed as a result of the adsorption process. Acetylene is close to the B11N12Ni cluster and is 2.146 Å apart. The C(2) of acetylene and the Ni atom of the cluster are close to each other, and the C≡C bond is slightly elongated with respect to ‘free’ acetylene (to a value of 1.215 Å). Following acetylene adsorption, the Ni atom is relatively more convex, with the length of the Ni–N(2) bond increasing to a value 1.942 Å and that of the Ni–N(3) bond increasing to 1.952 Å. The value of the N(2) –Ni–N(3) bond angle is reduced to 101.8°. Acetylene is physically adsorbed onto the B11N12Ni cluster, with the value for the adsorption energy calculated to be −11.92 kcal/mol. The Mulliken charge for C(1) has a value of −0.094 and for C(2) is −0.160; in this case, the Mulliken charge for Ni has a value of −0.163; consequently, the charge transfer from acetylene to the Ni centre.
In Figure 4c and Table S6 a stable configuration of acetylene adsorbed onto the B12N11Ni cluster is depicted. As can be evinced from this figure, the configurations of acetylene and of the B12N11Ni cluster are significantly different from their “free” counterparts, with acetylene undergoing an obvious deformation as a result of the adsorption process. Acetylene’s C(1) is adsorbed onto the Ni atom of the cluster, with a value for the length of the Ni–C(1) bond of 1.894 Å. Acetylene’s C(2) is adsorbed onto the B(1) atom of the cluster, and the C(2)–B(1) bond length has a value of 1.857 Å. The C–C bond is elongated to a value of 1.344 Å, and the Ni atom in the B12N11Ni cluster is convex to C(1). Following acetylene adsorption, the length of Ni–B(1) bond increases to 2.399Å and that of the Ni-B(2) bond increases to 2.035Å. Furthermore, the value for the B(1)–Ni–B(2) bond angle decreases to 88.19°. Acetylene is chemisorbed on the B12N11Ni cluster, and the adsorption energy is calculated to have a value of −40.40 kcal/mol. The Mulliken charge for C(2) has a value of −0.061 and that for B(1) is 0.584, as a result of the charge transfer from C(2) to B(1). The Mulliken charge for C(1) has a value of −0.497 and that for the Ni atom is −0.515, with the charge being transferred from the Ni atom to C(1); in this case, the Mulliken charge for B(2) has a value of 0.922, and charge is transferred from B(2) to the Ni atom.
In Figure 4d and Table S7 are depicted the optimised configuration of hydrogen adsorbed on the B12N12 cluster. Following this adsorption process, the configurations of hydrogen and cluster remain almost unchanged. Hydrogen and the B12N12 cluster are quite far from each other. Hydrogen adsorption on the B12N12 cluster is physical in nature, and the adsorption energy is calculated to have a value of −0.05 kcal/mol.
In Figure 4e and Table S8, the stable configuration of hydrogen adsorbed on the B11N12Ni cluster is depicted. As can be evinced from this figure, the configurations of hydrogen and the cluster change significantly following adsorption. Hydrogen is completely dissociated on the B11N12Ni cluster, and the adsorption energy is −1.90 kcal/mol. The distance between the two hydrogen atoms is 2.793 Å. H1 is adsorbed on the N(2) atom, and the length of the H–N bond is 1.019 Å. H(2) is adsorbed on the Ni atom, and the length of the H–Ni bond is 1.409 Å. The Ni atom protrudes away from the N(1) atom, and the value for the N(2)–Ni–N(3) bond angle is reduced to 98.28°. The Mulliken charge for H(2) has a value of 0.137, as a result of the charge transfer from H(2) to the Ni centre. The Mulliken charge for N(1) has a value of −0.834 and that for H(1) is 0.459, with the charge being transferred from H(1) to N(1).
In Figure 4f and Table S9, the stable configuration of hydrogen adsorbed on the B12N11Ni cluster is depicted. As can be evinced from this figure, the configurations of hydrogen and the B12N11Ni cluster change very significantly following adsorption. Hydrogen is completely dissociated on the B12N11Ni cluster, and the adsorption energy is calculated to have a value of −4.28 kcal/mol. The distance between the two H2-derived hydrogen atoms is 2.379 Å. H(1) is adsorbed on B(1), and the length of the H(1)–B(1) bond is 1.217 Å. H(2) is adsorbed on the Ni atom, and the length of the H(2)–Ni bond is 1.557 Å. The Ni atom protrudes away from B(1), and the value for the N(2)–Ni–N(3) bond angle is reduced to 80.62°. The value of the Mulliken charge for H(2) is −0.089, as a result of the charge transfer from H(2) to the Ni atom. The Mulliken charge for H(1) is 0.079 and that for B(1) is 1.633; in this case, charge is transferred from H(1) to B(1).

2.3. Acetylene Selective Hydrogenation Reaction Mechanism

The addition reaction between acetylene and hydrogen taking place on the catalyst’s surface starts with the stable co-adsorption (CO) of the two reactants on the said surface; it then passes through several transition states (TS) and intermediates (IM), to finally yield the reaction product (P).

2.3.1. Effect of Co-Adsorption Structure on Hydrogenation of Acetylene to Ethylene

To study the effect that different co-adsorbed structures have on acetylene hydrogenation, we investigated the mechanism of acetylene hydrogenation under two different co-adsorptions on the surface of the B12N11Ni cluster. In Figure 5a and Tables S10–S16, the configuration changes associated with each step of reaction pathway for the hydrogenation of acetylene to ethylene taking place on the surface of the B12N11Ni cluster are depicted. In Figure 5b, the energy changes associated with each step of the pathway is reported.
The reaction pathway starts with co-adsorption configuration CO involving the B12N11Ni cluster; initially, the two carbon atoms of acetylene are adsorbed on the Ni centre and boron atom B(2); on the other hand, dihydrogen is adsorbed around the Ni centre, with the two H atoms partially dissociating from each other. The distance between the two H2-derived hydrogen atoms is slightly larger than the isolation state. The adsorption energy of the co-adsorbed structure is −42.41 kcal/mol. H(2) vibrates in the C(1) direction, and H(1) vibrates in the B(3) direction to produce the transition state TS1 (characterised by a single imaginary frequency of −1248.70 cm−1), then reaction intermediate IM1 is produced. In it, H(1) has been adsorbed on B(3), H(2) has been added to C(1) and the distance between C(1) and Ni has increased with respect to the case of the co-adsorbed structure CO. Notably, the value for the energy barrier separating TS1 from CO is 6.54 kcal/mol. Next, H(1) attacks C(1) to produce transition state TS2, the species characterised by a single imaginary frequency (−356.31 cm−1); notably, the value for the energy barrier separating IM1 from TS2 is 26.07 kcal/mol. Intermediate IM2 is subsequently produced, whereby the CH3CH group has been adsorbed on the B12N11Ni cluster. Finally, H(3) migrates from C(1) to C(2), and ethylene is produced by way of transition state TS3 (characterised by a single imaginary frequency of −1028.14 cm−1). The process whereby IM2 turns into TS3 is the rate controlling step of the reaction pathway, and the overall pathway’s activation energy is calculated as 41.79 kcal/mol.
In another reaction route, the reaction pathway starts with co-adsorption configuration CO2. In Figure 6a and Tables S17–S30 are depicted the configuration changes associated with the various steps of the pathway and the energy changes are reported in Figure 6b. The adsorption energy of CO2 is calculated as −51.87kcal/mol. The two acetylene carbon atoms get adsorbed on the Ni atom, and the two H2-derived hydrogen atoms completely dissociate and get individually adsorbed onto the two B atoms bound to the nickel centre. Subsequently, the transition state TS1 is produced, which corresponds to the structure characterised by only one imaginary frequency (−1193.64 cm−1). In this structure, the B–H bond is elongated to a value of 1.736 Å, and H(2) attacks C(2) to produce the reaction intermediate IM1, which is characterised by an energy barrier of 19.51 kcal/mol. The H(1) atom, which is bound to B(1), then approaches C(1), and, after passing through the transition state TS2, which is the structure characterised by only one imaginary frequency (−1182.85 cm−1), the intermediate IM2 is produced. The process just described, whereby TS2 turns into IM2, is the rate-determining step in acetylene hydrogenation to ethylene. The energy barrier of the entire reaction is calculated to be 28.17 kcal/mol, whereas the desorption energy of ethylene is 17.15 kcal/mol.
Our data indicate that the activation energy of the second reaction pathway is lower than the first reaction pathway, so the second is the optimal pathway for acetylene hydrogenation to ethylene via catalysis by the B12N11Ni cluster. The starting points of these reaction pathways consist in different co-adsorbed structures. By comparing the data we collected on the two reaction pathways, we determined that the co-adsorbed structures would greatly affect the reaction.
In CO2, hydrogen is completely dissociated, and the H2-derived H atoms are adsorbed on the B atoms where the B atoms are nearby both sides of the two acetylenic carbons. The adsorption positions of the hydrogen atoms render favourable these atoms’ attack on different carbon atoms. In CO, as acetylene occupies the adsorption site on the catalyst, the two H2-derived hydrogen atoms can only be adsorbed around the Ni atom. After H(1) is added to C(1), the vinyl group rotates away from the Ni centre, as a result of an increase in steric hindrance, and H(2) can only be added to the C(1) atom, which is the closest to it, to produce the CH3CH moiety. For the migration of H atoms between two carbon atoms to take place, an extremely high energy barrier needs to be overcome, which accounts for the higher activation energy of the first pathway. What is why a suitable co-adsorbing structure can facilitate the reaction, and an unsuited co-adsorbing structure can be adverse.

2.3.2. Effect of Different Ni-Doped Catalysts on the Acetylene Hydrogenation Activity

To compare the effect that different Ni doping sites on the B12N12 clusters have on the acetylene hydrogenation activity, we also explored the mechanism of acetylene hydrogenation to ethylene on the B11N12Ni cluster. In Figure 7a and Tables S31–S42 are depicted the configuration changes associated with the various steps of the optimal pathway on the B11N12Ni cluster and the energy changes are reported in Figure 7b. Herein, acetylene is adsorbed around the Ni atom, and the two H2-derived hydrogen atoms are completely dissociated, with each of them adsorbed on either the Ni atom or one of the nitrogen atoms bound to the metal centre. The value for the co-adsorption energy is calculated to be −17.60 kcal/mol. The only imaginary frequency of the first transition state TS1 is −237.87 cm−1. In this state, H(1) attacks C(1) to produce the intermediate IM1, whose formation is characterised by an energy barrier of 2.13 kcal/mol. H(2), which is attached to N(1), then attacks C(2), and, after passing through transition state TS2 characterised by a single imaginary frequency of −1890.66 cm−1, the intermediate IM2 is produced. The process whereby TS2 is transformed into IM2 is the rate determining step in the hydrogenation of acetylene to ethylene; the energy barrier of reaction pathway is calculated to be 26.54 kcal/mol, whereas the desorption energy of ethylene is 9.75 kcal/mol.
Our data indicate that the activation energy on the B11N12Ni cluster is lower than the B12N11Ni cluster, so the hydrogenation of acetylene to ethylene on B11N12Ni clusters has better reaction activity. This shows that compared with the structure in which nickel is coordinated with three boron atoms, nickel is more active in coordination with three nitrogen atoms. Similarly, the desorption energy of ethylene on the B11N12Ni cluster is lower than the B12N11Ni cluster. It demonstrates that ethylene is more easily desorbed from the B11N12Ni cluster, thereby reducing the possibility of deep hydrogenation. We surmise that the B11N12Ni cluster is more ethylene-selective than the B12N11Ni cluster.

2.3.3. Selectivity of Acetylene Hydrogenation to Ethylene

For the purpose of determining whether the B11N12Ni and B12N11Ni clusters have ethylene selectivity, deep hydrogenation was continued to produce ethane after getting ethylene and compared its activation energy with that of ethylene.
In Figure 7, we hydrogenated intermediate IM2 further, so as to produce intermediate IM3, which is characterised by ethylene and hydrogen co-adsorbed on the B11N12Ni cluster. In this intermediate, hydrogen is completely dissociated: H(5) is adsorbed on Ni and H(6) is adsorbed on N(3); the distance between Ni and N(3) has increased from 1.952 Å to 3.099 Å with respect to IM2. After transition state TS3, which is characterised by having only one imaginary frequency of −185.86 cm−1, is produced, H(5) is added to C(1) to produce intermediate IM4. At this time, the distance between H(6) and C(2) is large. Following the ensuing formation of transition state TS4 (characterised a single imaginary frequency of −428.73 cm−1), H(6) rotates in the C(2) direction so that intermediate IM5 forms, wherein Ni and N(3) are close in space to each other. After the ensuing formation of transition state TS5, which is characterised by a single imaginary frequency of −1535.98 cm−1, H(6) is added to C(2) to form ethane. The energy required for the described process to take place is calculated to be 34.04 kcal/mol, and this process is the rate-controlling step in the generation of ethane. As the activation energy required to form ethane is greater than that required for ethylene production (26.54 kcal/mol), the hydrogenation of acetylene occurring on the B11N12Ni cluster is a selective process, with ethylene being the favoured product, with a calculated value for the selectivity of 24.29.
To compare the effect that different Ni doping sites on the B12N12 clusters have on the ethylene selectivity of the acetylene hydrogenation, we also explored the mechanism of acetylene hydrogenation to ethane after the formation of ethylene on the B12N11Ni cluster. As can be inferred from the data reported in Figure 6, intermediate IM2 is further hydrogenated to form intermediate IM3. In this intermediate, dihydrogen has dissociated completely, and the two resulting hydrogen atoms are adsorbed on different boron atoms located in the vicinity of the nickel centre. Ethylene, at this stage, is not obviously deformed and adsorbed on the top of the Ni atom. The distance between Ni and B(2) increases with respect to the IM2 case to 2.026 Å. After the transition state TS3, which is characterised by a single imaginary frequency of −154.31 cm−1, is produced, the C–C bond of ethylene increases in length, and C(1) approaches B(2), leading to the formation of intermediate IM4. The ensuing formation of intermediate IM5 is achieved by the preliminary formation of transition state TS4 (characterised by a single imaginary frequency of −1151.81 cm−1), whereby H(5) is added to ethylic carbon C(2). Notably, in IM5 C(2) rotates away from the N atom. This rotation may be due to the increase in steric hindrance associated with the attachment of three hydrogen atoms to C(2). After transition state TS5, characterised by a single imaginary frequency of −302.29 cm−1, has been produced, H(6) migrates from B(1) to the vicinity of B(2), the atom to which the ethyl group is attached, and the intermediate IM6 is produced. The energy required for this process to occur is calculated to be 34.65 kcal/mol, and this process is the rate-limiting step in ethane formation. Then transition state TS6 (characterised by a single imaginary frequency of −1063.04 cm−1) is produced, and finally H(6) is added to C(1) to form ethane. As the activation energy required for ethane formation is greater than that required for ethylene formation (28.17 kcal/mol), the acetylene hydrogenation reaction taking place on the B12N11Ni cluster is selective for ethylene as the product and the value for the selectivity of ethylene 17.49.

3. Computational Methods

In the present study, the Gaussian 09 software package was used [30], and all calculations were processed implementing the B3LYP hybrid density functional method [31,32,33,34,35]. The 6-31+G(d,p) basis set [36,37] was used for the B, C, N and H atoms, and the LANL2DZ base group [38,39,40] was used for the Ni atom. The counterpoise method was used to correct the base set superposition error [41]. The long-range dispersion correction was performed by the Grimme empirical dispersion correction method (DFT-D3) [42]. Frequency calculations were implemented to optimise the structures of the catalysts, reactants, intermediates, transition states and product in the absence of symmetry constraints. Notably, the transition state is the only species characterised by only one imaginary frequency. At the same time, intrinsic reaction coordinate calculations were performed in order to verify whether the transition state can be connected to the two intermediates in the reaction pathway [43,44]. All the energies involved in the calculation are taken as the zero-point corrected electronic energies, taking the Basis Set Superposition Error into consideration. Wave function analysis was performed using the Multiwfn program [45,46]. The possible reactive sites, adsorption configurations and charge changes of the reactants were explored through analysis of the frontier molecular orbitals, density of states, molecular electrostatic potentials and atomic charges of BNNi clusters.
The adsorption energy (Ead) was calculated employing Equation (1):
Ead = Eadsorption state − (Ehydrogen/acetylene + Ecatalyst),
The value of the catalyst’s selectivity can be obtained employing Equation (2):
Selectivity = Eethane barrier − Eethylene desorption,

4. Conclusions

The results of our investigation indicate that the acetylene hydrogenation reactions catalysed by the above mentioned single-nickel-doped clusters conform to the Horiuti–Polanyi mechanism [47]. The activity of the B11N12Ni and B12N11Ni clusters differs by only 1.63 kcal/mol, but the B11N12Ni cluster is expected to display higher ethylene product selectivity than the B12N11Ni cluster. This difference may be the result of the moderate degree of acetylene and hydrogen adsorption on the B11N12Ni cluster based on the Sabatier Principle, a different scenario to those of the undoped B12N12 and the B12N11Ni clusters. The low intensity of the described interaction is beneficial to the progress of the hydrogenation reaction, and it does not cause ethylene desorption to become difficult. Compared with the catalytic systems investigated in previously published studies, whose results are reported in Table 3 [17,48,49,50], the B11N12Ni cluster has moderate acetylene hydrogenation activity and increased ethylene selectivity. In particular, the hydrogenation activity of the B11N12Ni cluster is almost equivalent to that of the B11N12Pd cluster with similar structure, but the ethylene selectivity is much better. The B11N12Ni cluster has certain development prospects in catalytic the selective hydrogenation of acetylene by non-precious metals.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4344/10/1/115/s1, Table S1. The Cartesian coordinates of B12N12 cluster, Table S2. The Cartesian coordinates of B11N12Ni cluster, Table S3. The Cartesian coordinates of B12N11Ni cluster, Table S4. The Cartesian coordinates of acetylene adsorption on B12N12 cluster, Table S5. The Cartesian coordinates of acetylene adsorption on B11N12Ni cluster, Table S6. The Cartesian coordinates of acetylene adsorption on B12N11Ni cluster, Table S7. The Cartesian coordinates of hydrogen adsorption on B12N12 cluster, Table S8. The Cartesian coordinates of hydrogen adsorption on B11N12Ni cluster, Table S9. The Cartesian coordinates of hydrogen adsorption on B12N11Ni cluster, Table S10. The Cartesian coordinates of CO on B12N11Ni cluster, Table S11. The Cartesian coordinates of TS1 on B12N11Ni cluster, Table S12. The Cartesian coordinates of IM1 on B12N11Ni cluster, Table S13. The Cartesian coordinates of TS2 on B12N11Ni cluster, Table S14. The Cartesian coordinates of IM2 on B12N11Ni cluster, Table S15. The Cartesian coordinates of TS3 on B12N11Ni cluster, Table S16. The Cartesian coordinates of IM3 on B12N11Ni cluster, Table S17. The Cartesian coordinates of CO2 on B12N11Ni cluster, Table S18. The Cartesian coordinates of TS1 on B12N11Ni cluster, Table S19. The Cartesian coordinates of IM1 on B12N11Ni cluster, Table S20. The Cartesian coordinates of TS2 on B12N11Ni cluster, Table S21. The Cartesian coordinates of IM2 on B12N11Ni cluster, Table S22. The Cartesian coordinates of IM3 on B12N11Ni cluster, Table S23. The Cartesian coordinates of TS3 on B12N11Ni cluster, Table S24. The Cartesian coordinates of IM4 on B12N11Ni cluster, Table S25. The Cartesian coordinates of TS4 on B12N11Ni cluster, Table S26. The Cartesian coordinates of IM5 on B12N11Ni cluster, Table S27. The Cartesian coordinates of TS5 on B12N11Ni cluster, Table S28. The Cartesian coordinates of IM6 on B12N11Ni cluster, Table S29. The Cartesian coordinates of TS6 on B12N11Ni cluster, Table S30. The Cartesian coordinates of IM7 on B12N11Ni cluster, Table S31. The Cartesian coordinates of CO on B11N12Ni cluster, Table S32. The Cartesian coordinates of TS1 on B11N12Ni cluster, Table S33. The Cartesian coordinates of IM1 on B11N12Ni cluster, Table S34. The Cartesian coordinates of TS2 on B11N12Ni cluster, Table S35. The Cartesian coordinates of IM2 on B11N12Ni cluster, Table S36. The Cartesian coordinates of IM3 on B11N12Ni cluster, Table S37. The Cartesian coordinates of TS3 on B11N12Ni cluster, Table S38. The Cartesian coordinates of IM4 on B11N12Ni cluster, Table S39. The Cartesian coordinates of TS4 on B11N12Ni cluster, Table S40. The Cartesian coordinates of IM5 on B11N12Ni cluster, Table S41. The Cartesian coordinates of TS5 on B11N12Ni cluster, Table S42. The Cartesian coordinates of IM6 on B11N12Ni cluster.

Author Contributions

Writing—original draft, Y.W.; Writing—review & editing, L.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lindlar, V.H. Lindlar Ein neuer Katalysator für selektive Hydrierungen. Helv. Chim. Acta 1952, 35, 446–450. [Google Scholar] [CrossRef]
  2. Borodziński, A.; Geoffrey, C.B. Selective hydrogenation of ethyne in ethene-rich streams on palladium catalysts. Part 1. Effect of changes to the catalyst during reaction. Catal. Rev. 2006, 48, 91–144. [Google Scholar] [CrossRef]
  3. Zhao, Y.; Zhu, M.; Kang, L. The DFT study of Si-doped Pd6Si clusters for selective acetylene hydrogenation reaction. J. Mol. Graph. Model. 2018, 83, 129–137. [Google Scholar] [CrossRef] [PubMed]
  4. Zhao, Y.; Zhao, F.; Kang, L. Catalysis of the acetylene hydrochlorination reaction by Si-doped Au clusters: A DFT study. J. Mol. Model. 2018, 24, 61. [Google Scholar] [CrossRef] [PubMed]
  5. Zhou, H.; Yang, X.; Li, L.; Liu, X.; Huang, Y.; Pan, X.; Wang, A.; Li, J.; Zhang, T. PdZn intermetallic nanostructure with Pd–Zn–Pd ensembles for highly active and chemoselective semi-hydrogenation of acetylene. ACS Catal. 2016, 6, 1054–1061. [Google Scholar] [CrossRef]
  6. Zhang, Q.; Li, J.; Liu, X.; Zhu, Q. Synergetic effect of Pd and Ag dispersed on Al2O3 in the selective hydrogenation of acetylene. Appl. Catal. A Gen. 2000, 197, 221–228. [Google Scholar] [CrossRef]
  7. Studt, F.; Frank, A.P.; Thomas, B.; Rasmus, Z.S.; Claus, H.C.; Jens, K.N. On the role of surface modifications of palladium catalysts in the selective hydrogenation of acetylene. Angew. Chem. Int. Ed. 2008, 47, 9299–9302. [Google Scholar] [CrossRef]
  8. Shin, E.W.; Chang, H.C.; Ki, S.C.; Yoon, H.N.; Sang, H.M. Properties of Si-modified Pd catalyst for selective hydrogenation of acetylene. Catal. Today 1998, 44, 137–143. [Google Scholar] [CrossRef]
  9. McCue, A.J.; James, A.A. Recent advances in selective acetylene hydrogenation using palladium containing catalysts. Front. Chem. Sci. Eng. 2015, 9, 142–153. [Google Scholar] [CrossRef]
  10. Hageluken, C. Markets for the catalyst metals platinum, palladium and rhodium. Metall 2006, 60, 31. [Google Scholar]
  11. Janet, K.; Christine, M.; Detlef, K.; Mangelsdorf, I. Palladium—A review of exposure and effects to human health. Int. J. Hyg. Environ. Health 2002, 205, 417–432. [Google Scholar] [CrossRef]
  12. Qiao, B.; Wang, A.; Yang, X.; Allard, L.F.; Jiang, Z.; Cui, Y.; Liu, J.; Li, J.; Zhang, T. Single-atom catalysis of CO oxidation using Pt 1/FeO x. Nat. Chem. 2011, 3, 634. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, A.; Li, J.; Zhang, T. Heterogeneous single-atom catalysis. Nat. Rev. Chem. 2018, 2, 65–81. [Google Scholar] [CrossRef]
  14. Huang, X.; Xia, Y.; Cao, Y.; Zheng, X.; Pan, H.; Zhu, J.; Ma, C.; Wang, H.; Li, J.; You, R.; et al. Enhancing both selectivity and coking-resistance of a single-atom Pd1/C3N4 catalyst for acetylene hydrogenation. Nano Res. 2017, 10, 1302–1312. [Google Scholar] [CrossRef]
  15. Yu, Z.; Mingyuan, Z.; Lihua, K. The DFT Study of Single-Atom Pd 1/gC 3 N 4 Catalyst for Selective Acetylene Hydrogenation Reaction. Catal. Lett. 2018, 148, 2992–3002. [Google Scholar] [CrossRef]
  16. Zhou, S.; Shang, L.; Zhao, Y.; Shi, R.; Waterhouse, G.I.; Huang, Y.C.; Zheng, L.; Zhang, T. Pd Single-Atom Catalysts on Nitrogen-Doped Graphene for the Highly Selective Photothermal Hydrogenation of Acetylene to Ethylene. Adv. Mater. 2019, 31, 1900509. [Google Scholar] [CrossRef]
  17. Gong, W.; Kang, L. BNPd single-atom catalysts for selective hydrogenation of acetylene to ethylene: A density functional theory study. R. Soc. Open Sci. 2018, 5, 171598. [Google Scholar] [CrossRef] [Green Version]
  18. Ronnie, V.T.; Karoliina, H.; Søren, D.; Ebbe, K.V.; Joachim, S.; Erik, L.; Bjerne, S.C.; Jens, K.N.; Flemming, B. Controlling the catalytic bond-breaking selectivity of Ni surfaces by step blocking. Nat. Mater. 2005, 4, 160. [Google Scholar] [CrossRef]
  19. Ronnie, V.T.; Karoliina, H.; Søren, D.; Ebbe, K.V.; Joachim, S.; Erik, L.; Bjerne, S.C.; Jens, K.N.; Flemming, B. Ethylene dissociation on flat and stepped Ni (1 1 1): A combined STM and DFT study. Surf. Sci. 2006, 600, 66–77. [Google Scholar] [CrossRef]
  20. Mueller, J.E.; Adri, C.V.D.; William, A.G. Structures, energetics, and reaction barriers for CH x bound to the nickel (111) surface. J. Phys. Chem. C 2009, 113, 20290–20306. [Google Scholar] [CrossRef] [Green Version]
  21. Ledoux, M.J.; Pham-Huu, C. Carbon nanostructures with macroscopic shaping for catalytic applications. Catal. Today 2005, 102, 2–14. [Google Scholar] [CrossRef]
  22. Stephan, O.; Bando, Y.; Loiseau, A.; Willaime, F.; Shramchenko, N.; Tamiya, T.; Sato, T. Formation of small single-layer and nested BN cages under electron irradiation of nanotubes and bulk material. Appl. Phys. A Mater. 1998, 67. [Google Scholar] [CrossRef]
  23. Strout, D.L. Structure and stability of boron nitrides: Isomers of B12N12. J. Phys. Chem. A 2000, 104, 3364–3366. [Google Scholar] [CrossRef]
  24. Koi, N.; Takeo, O. Hydrogen storage in boron nitride and carbon clusters studied by molecular orbital calculations. Solid State Commun. 2004, 131, 121–124. [Google Scholar] [CrossRef]
  25. Takeo, O.; Atsushi, N.; Ichihito, N. Formation and atomic structure of B12N12 nanocage clusters studied by mass spectrometry and cluster calculation. Sci. Technol. Adv. Mater. 2004, 5, 635. [Google Scholar] [CrossRef]
  26. Takeo, O.; Ichihito, N.; Atsushi, N. Synthesis, atomic structures, and electronic states of boron nitride nanocage clusters and nanotubes. Mater. Manuf. Process. 2004, 19, 1215–1239. [Google Scholar] [CrossRef]
  27. Hongyu, W.; Xiaofeng, F.; Jer-Lai, K. Metal free hydrogenation reaction on carbon doped boron nitride fullerene: A DFT study on the kinetic issue. Int. J. Hydrogen Energy 2012, 37, 14336–14342. [Google Scholar] [CrossRef]
  28. Rad, A.S.; Khurshid, A. Enhancement in hydrogen molecule adsorption on B12N12 nano-cluster by decoration of nickel. Int. J. Hydrogen Energy 2016, 41, 22182–22191. [Google Scholar] [CrossRef]
  29. Zhao, F.; Wang, Y.; Zhu, M.; Kang, L. C-doped boron nitride fullerene as a novel catalyst for acetylene hydrochlorination: A DFT study. RSC Adv. 2015, 5, 56348–56355. [Google Scholar] [CrossRef]
  30. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09; Revision, A.1; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  31. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098. [Google Scholar] [CrossRef]
  32. Chengteh, L.; Weitao, Y.; Robert, G.P. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785. [Google Scholar] [CrossRef] [Green Version]
  33. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648. [Google Scholar] [CrossRef] [Green Version]
  34. Stephens, P.J.; Devlin, F.J.; Chabalowski, C.F.N.; Frisch, M.J. Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields. J. Phys. Chem. 1994, 98, 11623–11627. [Google Scholar] [CrossRef]
  35. Vosko, S.H.; Leslie, W.; Marwan, N. Accurate spin-dependent electron liquid correlation energies for local spin density calculations: A critical analysis. Can. J. Phys. 1980, 58, 1200–1211. [Google Scholar] [CrossRef] [Green Version]
  36. Timothy, C.; Jayaraman, C.; Günther, W.S. Paul Von Ragué Schleyer Efficient diffuse function-augmented basis-sets for anion calculations. 3. The 3-21+G basis set for 1st-row elements, Li-F. J. Comp. Chem. 1983, 4, 294–301. [Google Scholar] [CrossRef]
  37. Petersson, A.; Bennett, A.; Tensfeldt, T.G.; Al-Laham, M.A.; Shirley, W.A.; Mantzaris, J. A complete basis set model chemistry. I. The total energies of closed-shell atoms and hydrides of the first-row elements. J. Chem. Phys. 1988, 89, 2193–2218. [Google Scholar] [CrossRef]
  38. Hay, P.J.; Willard, R.W. Ab initio effective core potentials for molecular calculations. Potentials for the transition metal atoms Sc to Hg. J. Chem. Phys. 1985, 82, 270–283. [Google Scholar] [CrossRef]
  39. Willard, R.W.; Hay, P.J. Ab initio effective core potentials for molecular calculations. Potentials for main group elements Na to Bi. J. Chem. Phys. 1985, 82, 284–298. [Google Scholar] [CrossRef]
  40. Hay, P.J.; Willard, R.W. Ab initio effective core potentials for molecular calculations. Potentials for K to Au including the outermost core orbitals. J. Chem. Phys. 1985, 82, 299–310. [Google Scholar] [CrossRef]
  41. van Duijneveldt, F.B.; van Duijneveldt-van de Rijdt, J.G.; van Lenthe, J.H. State of the art in counterpoise theory. Chem. Rev. 1994, 94, 1873–1885. [Google Scholar] [CrossRef]
  42. Stefan, G.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [Green Version]
  43. Carlos, G.; Bernhard, H.S. An improved algorithm for reaction path following. J. Chem. Phys. 1989, 90, 2154–2161. [Google Scholar] [CrossRef]
  44. Carlos, G.; Bernhard, H.S. Reaction path following in mass-weighted internal coordinates. J. Phys. Chem. 1990, 94, 5523–5527. [Google Scholar] [CrossRef]
  45. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  46. Lu, T.; Chen, F. Calculation of Molecular Orbital Composition. Acta Chim. Sin. 2011, 69, 2393–2406. [Google Scholar]
  47. Iurô, H.; Polanyi, M. Exchange reactions of hydrogen on metallic catalysts. Trans. Faraday Soc. 1934, 30, 1164–1172. [Google Scholar] [CrossRef]
  48. Rao, D.M.; Sun, T.; Yang, Y.S.; Yin, P.; Pu, M.; Yan, H.; Wei, M. Theoretical study on the reaction mechanism and selectivity of acetylene semi-hydrogenation on Ni–Sn intermetallic catalysts. Phys. Chem. Chem. Phys. 2019, 21, 1384–1392. [Google Scholar] [CrossRef]
  49. Rao, D.M.; Zhang, S.T.; Li, C.M.; Chen, Y.D.; Pu, M.; Yan, H.; Wei, M. The reaction mechanism and selectivity of acetylene hydrogenation over Ni–Ga intermetallic compound catalysts: A density functional theory study. Dalton Trans. 2018, 47, 4198–4208. [Google Scholar] [CrossRef]
  50. Abdollahi, T.; Davood, F. Activity and selectivity of Ni nanoclusters in the selective hydrogenation of acetylene: A computational investigation. C. R. Chim. 2018, 21, 484–493. [Google Scholar] [CrossRef]
Catalysts 10 00115 g001 550
Figure 1. The optimised structures of undoped B12N12 cluster (a), B11N12Ni (b) and B12N11Ni (c). H atom is white; B atom is pink; C atom is grey; N atom is mazarine; Ni atom is blue.
Figure 1. The optimised structures of undoped B12N12 cluster (a), B11N12Ni (b) and B12N11Ni (c). H atom is white; B atom is pink; C atom is grey; N atom is mazarine; Ni atom is blue.
Catalysts 10 00115 g001
Catalysts 10 00115 g002 550
Figure 2. The molecular electrostatic potential of the B12N12 (a), B11N12Ni (b) and B12N11Ni (c) clusters on the molecular van der Waals surface (electron density: 0.001, energies in a.u.).
Figure 2. The molecular electrostatic potential of the B12N12 (a), B11N12Ni (b) and B12N11Ni (c) clusters on the molecular van der Waals surface (electron density: 0.001, energies in a.u.).
Catalysts 10 00115 g002
Catalysts 10 00115 g003 550
Figure 3. The TDOS and PDOS of B12N12 (a), B11N12Ni (b) and B12N11Ni (c) and their calculated HOMO and LUMO.
Figure 3. The TDOS and PDOS of B12N12 (a), B11N12Ni (b) and B12N11Ni (c) and their calculated HOMO and LUMO.
Catalysts 10 00115 g003
Catalysts 10 00115 g004 550
Figure 4. The optimised structures of C2H2 absorbed on B12N12 (a), B11N12Ni (b) and B12N11Ni (c); the optimised structures of H2 absorbed on B12N12 (d), B11N12Ni (e) and B12N11Ni (f). H atom is white; B atom is pink; C atom is grey; N atom is mazarine; Ni atom is blue.
Figure 4. The optimised structures of C2H2 absorbed on B12N12 (a), B11N12Ni (b) and B12N11Ni (c); the optimised structures of H2 absorbed on B12N12 (d), B11N12Ni (e) and B12N11Ni (f). H atom is white; B atom is pink; C atom is grey; N atom is mazarine; Ni atom is blue.
Catalysts 10 00115 g004
Catalysts 10 00115 g005 550
Figure 5. The reaction process and potential energy change for the different reaction pathways R3 of catalytic hydrogenation of acetylene to ethylene on B12N11Ni. (a) The structures of calculated hydrogenation reaction diagram for the B12N11Ni. (b) The reaction process and potential energy change of catalytic hydrogenation of ethylene to ethylene on the B12N11Ni. H atom is white; B atom is pink; C atom is grey; N atom is mazarine; Ni atom is blue.
Figure 5. The reaction process and potential energy change for the different reaction pathways R3 of catalytic hydrogenation of acetylene to ethylene on B12N11Ni. (a) The structures of calculated hydrogenation reaction diagram for the B12N11Ni. (b) The reaction process and potential energy change of catalytic hydrogenation of ethylene to ethylene on the B12N11Ni. H atom is white; B atom is pink; C atom is grey; N atom is mazarine; Ni atom is blue.
Catalysts 10 00115 g005
Catalysts 10 00115 g006 550
Figure 6. (a) The structures of calculated hydrogenation reaction diagram for the B12N11Ni. (b) The reaction process and potential energy change of catalytic hydrogenation of ethylene to ethane on the B12N11Ni. H atom is white; B atom is pink; C atom is grey; N atom is mazarine; Ni atom is blue.
Figure 6. (a) The structures of calculated hydrogenation reaction diagram for the B12N11Ni. (b) The reaction process and potential energy change of catalytic hydrogenation of ethylene to ethane on the B12N11Ni. H atom is white; B atom is pink; C atom is grey; N atom is mazarine; Ni atom is blue.
Catalysts 10 00115 g006
Catalysts 10 00115 g007 550
Figure 7. (a) The structures of calculated hydrogenation reaction diagram for the B11N12Ni. (b) The reaction process and potential energy change of catalytic hydrogenation of ethylene to ethane on the B12N11Ni. H atom is white; B atom is pink; C atom is grey; N atom is mazarine; Ni atom is blue.
Figure 7. (a) The structures of calculated hydrogenation reaction diagram for the B11N12Ni. (b) The reaction process and potential energy change of catalytic hydrogenation of ethylene to ethane on the B12N11Ni. H atom is white; B atom is pink; C atom is grey; N atom is mazarine; Ni atom is blue.
Catalysts 10 00115 g007
Table
Table 1. The orbital energies on the HOMO and LUMO of B12N12, B11N12Ni and B12N11Ni, and their energy gaps (energies in eV).
Table 1. The orbital energies on the HOMO and LUMO of B12N12, B11N12Ni and B12N11Ni, and their energy gaps (energies in eV).
HOMOLUMOΔEg
B12N12−7.95−1.246.72
B11N12Ni−7.08−5.111.97
B12N11Ni−5.18−2.522.66
Table
Table 2. The adsorption energies of C2H2 and H2 separately adsorb on B12N12, B11N12Ni and B12N11Ni (energies in kcal/mol).
Table 2. The adsorption energies of C2H2 and H2 separately adsorb on B12N12, B11N12Ni and B12N11Ni (energies in kcal/mol).
B12N12B11N12NiB12N11Ni
C2H2−2.83−11.92−40.40
H2−0.05−1.90−4.28
Table
Table 3. Comparison of the energy barriers for the desorption energy of ethylene and the selectivity with different catalysts (energies in kcal/mol).
Table 3. Comparison of the energy barriers for the desorption energy of ethylene and the selectivity with different catalysts (energies in kcal/mol).
CatalystEethylene barrierSelectivityB12N11Ni
B11N12Ni26.5424.29This work
B12N11Ni28.1717.49This work
Ni3Sn42.4312.91[29]
Ni3Sn312.223.23[29]
NiGa(111)24.4414.53[30]
Ni622.6015.35[31]
B11N12Pd26.5512.80[12]

Share and Cite

MDPI and ACS Style

Wang, Y.; Kang, L. Selective Hydrogenation of Acetylene Catalysed by a B12N12 Cluster Doped with a Single Nickel Atom: A DFT Study. Catalysts 2020, 10, 115. https://doi.org/10.3390/catal10010115

AMA Style

Wang Y, Kang L. Selective Hydrogenation of Acetylene Catalysed by a B12N12 Cluster Doped with a Single Nickel Atom: A DFT Study. Catalysts. 2020; 10(1):115. https://doi.org/10.3390/catal10010115

Chicago/Turabian Style

Wang, Yun, and Lihua Kang. 2020. "Selective Hydrogenation of Acetylene Catalysed by a B12N12 Cluster Doped with a Single Nickel Atom: A DFT Study" Catalysts 10, no. 1: 115. https://doi.org/10.3390/catal10010115

APA Style

Wang, Y., & Kang, L. (2020). Selective Hydrogenation of Acetylene Catalysed by a B12N12 Cluster Doped with a Single Nickel Atom: A DFT Study. Catalysts, 10(1), 115. https://doi.org/10.3390/catal10010115

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