Mechanistic Insight into the Hydrogenation of Acetylene on the Pd2/g-C3N4 Catalyst: Effect of Pd Clustering on the Barrier Energy and Selectivity

Mechanistic Insight into the Hydrogenation of Acetylene on the Pd2/g-C3N4 Catalyst: Effect of Pd Clustering on the Barrier Energy and Selectivity | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Mechanistic Insight into the Hydrogenation of Acetylene on the Pd2/g-C3N4 Catalyst: Effect of Pd Clustering on the Barrier Energy and Selectivity Sayyed Mahdi Hosseini, Mehran Ghiaci, Hossein Farrokhpour This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-259183/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 19 Apr, 2021 Read the published version in Structural Chemistry → Version 1 posted 3 You are reading this latest preprint version Abstract In this work, the hydrogenation of acetylene on the Pd 2 /g-C 3 N 4 catalyst is investigated by the Density Functional Theory (DFT) and Quantum Theory of Atoms in Molecules (QTAIM) calculations. The Pre-reactant (R), transition states (TSs), and the intermediates (IMs), involved in the hydrogenation process, are characterized from the point of view of energy and structure. The calculated energy barrier for the hydrogen transfer to the acetylene and ethylene are 6.77 and 12.28 kcal/mol, respectively which shows that the Pd 2 /g-C 3 N 4 catalyst has good selectivity for the conversion of acetylene to ethylene rather than ethane. Comparing the values of these energy barriers with those of the hydrogenation of acetylene on the Pd/g-C 3 N 4 catalyst (21.53 and 38.88 kcal/mol, respectively) shows that the increase in the number of the Pd atoms decreases the energy barriers of the hydrogenation reaction and increases the selectivity of the catalyst for the ethylene production. General Biochemistry Acetylene Hydrogenation Pd cluster g-C3N4 Selectivity Ethylene Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Figure 9 Introduction Ethylene is an important raw material with wide applications in petrochemical industries [1]. Generally, the production of ethylene is performed via the steam cracking of hydrocarbons. The main disadvantage of this method is the catalyst poisoning leading to the reduction of catalyst lifetime.[2]. So, the catalytic hydrogenation of acetylene followed by the classical Horiuti–Polanyi mechanism has been considered as the best method for ethylene production on large scales [3]. An efficient catalyst for acetylene hydrogenation should have high activity and selectivity for ethylene production [4]. In this regard, the noble transition metals, especially, palladium (Pd)-based catalysts are widely used for this purpose. Because they show good selectivity for ethylene production because of their superior hydrogen solubility and good ability for dissociation of H 2 molecules [5,6]. There are several experimental works aiming production of ethylene via acetylene hydrogenation in a selective manner in which the structure of the Pd-based catalyst plays a critical role in achieving a good selectivity [7-12]. Huang et al. prepared the dispersed Pd atoms on the g-C 3 N 4 using atomic layer deposition and performed the hydrogenation of acetylene on the Pd 1 /g-C 3 N 4 catalyst [7]. Boucher et al. performed the hydrogenation of styrene and acetylene on the catalyst composed of the Pd atoms binding to the Cu surface [8]. Panpranot et al. investigated the catalytic efficiency of Pd catalyst supported on nanocrystalline TiO 2 (Pd/TiO 2 ) for the selective hydrogenation of acetylene in the presence of excess ethylene [11]. The same group also examined the effect of the phase composition of the TiO 2 crystalline support on the physicochemical and catalytic properties of the Pd/TiO 2 catalyst on the selective hydrogenation of acetylene [12]. There are several theoretical studies of the clarification of the hydrogenation mechanism of acetylene on the Pd-based catalysts, and the interaction between the different forms of Pd and the substrates [13-16]. Meng et al. studied the hydrogenation of acetylene on the TiO 2 supported Pd a Ag b ( a + b =4) catalyst using the density functional theory (DFT)+U method [13]. The hydrogenation of acetylene using the Pd 4 cluster on the anatase TiO 2 (101) support was studied by Yang et al. [14]. Ma et al. performed the DFT study of the selective hydrogenation of acetylene on the P-doped Cu (111) surface [15]. The effect of the size of Pd cluster (Pd n ) and its geometry on the hydrogenation of acetylene has been investigated [17-20]. Li et al. investigated the effect of the size of Pd clusters on acetylene hydrogenation in the absence of any support [17]. Pu et al. investigated the effect of Pd n clusters with different sizes (n=2-8) on acetylene hydrogenation by the DFT calculations. They proposed two pathways for acetylene hydrogenation including ( a ) ethene formation as the main product via vinyl intermediate and ( b ) the vinylidene, ethylidene, and ethylidene as the intermediates participating in the hydrogenation process. Also, They showed that the reaction pathway depends strongly on the size of Pd n clusters so that for n <4, the hydrogenation reaction proceeded through the vinyl intermediate while the other pathway via the vinylidene intermediate was the main pathway for the larger Pd clusters [17]. Due to the high activity of Pd catalysts, the over-hydrogenation of ethylene is always a deep concern. Therefore, designing a selective catalyst for ethylene production to eliminate extra purification steps in its synthesis and prevent ethane production is considered a challenging task in acetylene hydrogenation [21-22]. Various strategies have been so far reported for the promotion of the performance of Pd-based catalysts for the acetylene hydrogenation which one of the important of them is the doping of the Pd cluster with other noble metals such as Au [23], Pb [24,25], Cu [26], Ag [27,28,29], and Ga [30,31]. The other important strategy is using single-atom alloys (SAA) and single-atom catalysts (SAC), immobilized on solid supports such as ZnO [32], fiberglass [33], graphene [34], Al 2 O 3 [35], and graphitic carbon nitride (g-C 3 N 4 ) [6-7] which has shown interesting results in acetylene hydrogenation. Among the solid supports, g-C 3 N 4 is always considered a suitable candidate for the deposition of Pd active sites [36-38]. It has some outstanding features such as low-cost fabrication, short time and easy synthesis, adjustable surface area, good dispersion potential for active sites, and high coordination ability for metal clusters and nanoparticles. Huang et.al showed that Pd single-atom immobilized on the surface of g-C 3 N 4 has a better selectivity in acetylene hydrogenation compared to nanoparticle ones [7]. Zhao et al.’s theoretical study aiming to demonstrate the acetylene hydrogenation mechanism over the Pd 1 /g-C 3 N 4 catalyst showed that the Pd atom tended to adsorb on the top of the sixth fold cavity of g-C 3 N 4 as the most accessible site. Moreover, the Pd 1 /g-C 3 N 4 catalyst showed better selectivity to ethylene production in comparison to ethane production as expected from the experimental studies [39]. Kang et al. investigated the catalytic activity of the Pd 1 /g-C 3 N 4 catalyst for acetylene hydrogenation by replacing the nitrogen and carbon atoms of g-C 3 N 4 with the sulfur atom (S) [40]. It was observed that the Pd 1 /g-C 3 N 4 containing sulfur-doped carbon atoms have higher activity for the dissociation of H 2 molecules while better selectivity is achieved for sulfur-doped nitrogen ones. To the best of our knowledge, there is no theoretical mechanistic study of acetylene hydrogenation over Pd n /g-C 3 N 4 structures with n>1 in the literature. Our recent theoretical study on the adsorption of small-size Pd n clusters (n=2, 4, 6, and 8), immobilized on the g-C 3 N 4 quantum dot, reveals that the Pd 2 cluster has the highest adsorption energy among the other selected Pd clusters [41]. In this work, the DFT calculations are performed to understand the exact mechanism of acetylene hydrogenation on the Pd 2 /g-C 3 N 4 catalyst and see the effect of the added Pd atom on the activity of catalyst compared to that of Pd 1 /g-C 3 N 4 catalyst reported in references 39 and 40. Our calculated results show less activation energy for acetylene hydrogenation using the Pd 2 /g-C 3 N 4 compared to the Pd 1 /g-C 3 N 4 catalyst and better selectivity for ethylene production which is in good accordance with the experimental data in the literature [39,45-47]. Computational details The structure of all reactants (Pd 2 /g-C 3 N 4 , C 2 H 2 , and H 2 ), pre-reactant (R), transition states (TSs), and intermediates (IMs) were optimized in the gas phase using the DFT method employing B3PW91 functional [42]. The 6-311+G(d,p) and Los Alamos effective core pseudo-potential (ECP) (LANL2DZ) basis set were used for the g-C 3 N 4 and Pd 2 , respectively. The frequency calculations were performed on all optimized structures to check their correct positions on their potential energy surfaces in the hydrogenation pathway and determine their electronic energies, enthalpies, Gibbs free energies, and zero-point energies. The suitability of the selected DFT functional and basis set for the considered systems have been investigated by Zhao et al [39]. The intrinsic reaction coordinate (IRC) calculations were performed to make sure that each TS structure was related to its corresponding reactants and products in the reaction mechanism. The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of the structures were calculated to obtain the bandgap energies and donor-acceptor properties of structures. The Quantum Theory of Atoms in Molecules (QTAIM) calculations were performed using Multiwfn 3.4 software [43] for better structural characterization of the TSs and confirming the proposed mechanism of the reaction. All DFT calculations were performed using the Gaussian 09 Quantum Chemistry Package [44]. Results And Discussion The optimized structure of the Pd 2 /g-C 3 N 4 catalyst and its charge distribution A quantum dot of g-C 3 N 4 is selected as a model of substrate for the deposition of the Pd 2 cluster. The g-C 3 N 4 has two main structural forms based on the s-triazine and tri-s-triazine building blocks so that the latter form is always considered as the most stable form of g-C 3 N 4 [39-41]. A layer of g-C 3 N 4 composed of three tri-s-triazine units was selected and optimized (see Figure 1). The nitrogen atoms in this model are classified into two classes: ( a ) saturated nitrogen atoms which are bonded to three carbon atoms (N1) and ( b ) pyridine-like nitrogen atoms which are bonded to two carbon atoms and are located at the corners of the sixth fold cavity (N2). Since 2p orbitals of N2 atoms have the main contribution in the valance band of g-C 3 N 4 , their lone pairs electrons have the main role in the chemical activity of g-C 3 N 4 compared to the N1 atoms [41]. On the other hand, 2p orbitals of the N1 and C atoms orbitals have a major contribution to the conduction band of g-C 3 N 4 [41]. The Pd 2 cluster is optimized, separately and placed on the top of the sixth fold cavity of g-C 3 N 4 aiming to achieve the maximum interaction with the lone pairs of N2 atoms. The comparison of the structure of g-C 3 N 4 in the Pd 2 /g-C 3 N 4 complex with the optimized structure of isolated g-C 3 N 4 shows that the presence of the Pd 2 creates more wrinkles and waves in the structure of g-C 3 N 4 (Figure 1). The bond length of the Pd 2 cluster in the optimized structure of the Pd 2 /g-C 3 N 4 is 2.74 Å which shows the elongation of the Pd-Pd bond length compared to the isolated Pd 2 cluster due to the strong interaction between the Pd atom, directed to the support, and the carbon atom of one of the tri-s-triazine units. The calculated partial atomic charges of these Pd and C atoms in the Pd 2 /g-C 3 N 4 is 0.37e and -1.458e, respectively. Comparison of the charge of this C atom in the complex with that in the isolated g-C 3 N 4 (0.297e) shows a strong charge transfer from the Pd 2 to g-C 3 N 4 confirming a strong interaction between the Pd and carbon atom (see Figure 2). The angle between the Pd 2 axis and the g-C 3 N 4 plane is 46.0°. HOMO and LUMO of reactants Table 1 reports the energies of the HOMO and LUMO of C 2 H 2 , H 2 , g-C 3 N 4 , Pd 2 /g-C 3 N 4 , and C 2 H 2 .Pd 2 /g-C 3 N 4 structures along with their energy gaps. The shapes of the HOMO and LUMO of g-C 3 N 4 , Pd 2 /g-C 3 N 4 , and C 2 H 2 .Pd 2 /g-C 3 N 4 structures have been shown in Figure 3. According to our previous work [41], immobilization of the Pd 2 cluster on the surface of g-C 3 N 4 causes the decrease of the energy gap of the Pd 2 /g-C 3 N 4 complex compared to the isolated g-C 3 N 4 because of creating new electronic states between the HOMO and LUMO of g-C 3 N 4 [41]. These new electronic states cause the increase of the photocatalytic activity of the Pd 2 /g-C 3 N 4 complex compared to the bare g-C 3 N 4 . Comparison of the difference between the energy of the HOMO of Pd 2 /g-C 3 N 4 and that of LUMO of C 2 H 2 (-4.25 eV) with the energy difference between the HOMO of C 2 H 2 and LUMO of Pd 2 /g-C 3 N 4 (-5.24 eV) shows that the C 2 H 2 is an electron acceptor species. A similar comparison between the Pd 2 /g-C 3 N 4 and H 2 shows that the H 2 is also an electron acceptor. The smaller value of the energy difference between the HOMO of C 2 H 2 and LUMO of Pd 2 /g-C 3 N 4 compared to the corresponding value for the H 2 and Pd 2 /g-C 3 N 4 indicates the better adsorption of the C 2 H 2 on the Pd 2 /g-C 3 N 4 compared to the H 2 molecule. Adsorption of H 2 and C 2 H 2 on Pd 2 /g-C 3 N 4 Essentially for a heterogeneous catalyst containing the Pd cluster, it is important to use the maximum active site of the cluster after immobilization. Therefore, it is necessary to find the most stable site for the adsorption of H 2 and C 2 H 2 molecules on the surface of the Pd 2 /g-C 3 N 4 catalyst. For this purpose, the H 2 molecule approaches the catalyst in the dissociated and bonded forms, separately (see Figures S1(a) and (b)). As seen in Figure S1(a), the H 2 molecule is located near one Pd atom which is farther from the surface of g-C 3 N 4 . In Figure S1(b), it is assumed that the H 2 molecule has been dissociated on the Pd 2 /g-C 3 N 4 catalyst and each of its H atoms is interacting with one of the Pd atoms. The two initial structures lead to the same structure after the optimization so that the two hydrogen atoms of the H 2 molecule is interacting with one Pd atom (H 2 .Pd 2 /g-C 3 N 4 ; see Figure S1(c)). The reason for the interaction of the H 2 molecule with this Pd atom could be attributed to its negative Mulliken charge which provides a good potential for the participation in electron donation to reactants. In the next step, the C 2 H 2 is approached to the H 2 .Pd 2 /g-C 3 N 4 complex to form the pre-reactant (C 2 H 2 .H 2 /g-C 3 N 4 ) for the start of the reaction. Figure S2 shows two proposed structures for the pre-reactant (R1 and R2). The best selected initial orientation of the C 2 H 2 relative to the Pd 2 in the C 2 H 2 .H 2 .Pd 2 /g-C 3 N 4 complex is such that its molecular axis is in a vertical position relative to the axis of the Pd 2 (R1 in Figure S2). In this case, the simultaneous interaction of each C atom of C 2 H 2 with both Pd atoms is possible. It will be shown that the interaction of the C atoms of C 2 H 2 with the Pd atom, located far away from the surface of g-C 3 N 4 , is stronger than the interaction with the other Pd atom. As previously shown in Figure 2b, the Pd atom located further away from the g-C 3 N 4 has a more negative Mulliken charge compared to the other Pd atom. Therefore, this Pd atom has good potential for electron donation to the reactants (H 2 and C 2 H 2 ). Therefore, the C 2 H 2 and H 2 interact mainly with the upper Pd atom in the optimized structure of C 2 H 2 .H 2 .Pd 2 /g-C 3 N 4 complex. Due to the adsorption, the H-H bond length increases to 0.93 Å compared to that of an isolated H 2 molecule (0.74 Å). Moreover, the increase of the C≡C bond length is seen for the acetylene in the C 2 H 2 .H 2 .Pd 2 /g-C 3 N 4 complex (1.28 Å) compared to isolated acetylene molecule (1.19 Å). These observations confirm the activation of the H 2 and C 2 H 2 molecules by the catalyst and indicate that the initial complex (R1 in Figure S2) is suitable for starting the hydrogenation process. Considering the adsorption energies of reactants, it can be concluded that the adsorption of acetylene on the catalyst (-28.84 kcal/mol) is stronger than that of the H 2 molecule (-13.23 kcal/mol) which is in agreement with the frontier molecular orbital analysis performed in the previous section. The energy difference between the HOMO of C 2 H 2 .Pd 2 /g-C 3 N 4 and the LUMO of H 2 (-5.93 eV) is smaller than the energy difference between the HOMO of H 2 and the LUMO of C 2 H 2 .Pd 2 /g-C 3 N 4 (-9.07 eV). This shows a charge transfer from the HOMO of C 2 H 2 .Pd 2 /g-C 3 N 4 complex to the LUMO of H 2 molecule and leads to a strong interaction between these reactants and can better describe the co-adsorption configuration of the H 2 .C 2 H 2 .Pd 2 /g-C 3 N 4 complex.The calculated co-adsorption energy of the H 2 .C 2 H 2 .Pd 2 /g-C 3 N 4 complex is -45.3 kcal/mol which is greater than the value related to the adsorption of isolated C 2 H 2 (-28.84 kcal/mol) and H 2 (-13.23 kcal/mol) molecules on the Pd 2 /g-C 3 N 4 catalyst. Eq. (1) is used for the calculation of the co-adsorption energy of H 2 .C 2 H 2 .Pd 2 /g-C 3 N 4 complex. where E complex is the electronic energy of H 2 .C 2 H 2 .Pd 2 /g-C 3 N 4 complex. The E acetylene , E catalyst , and E H2 are the electronic energies of isolated optimized acetylene, Pd 2 /g-C 3 N 4, and H 2 molecule, respectively Mechanism of acetylene hydrogenation on the Pd 2 /g-C 3 N 4 catalyst As previously mentioned, two initial structures, R1 and R2 (Figure S2), were considered as the starting structures for the hydrogenation mechanism. Zhao et al. demonstrated that the hydrogenation pathway via the R1 complex in the presence of Pd 1 /g-C 3 N 4 single-atom catalyst requires less activation energy compared to the R2 ones [39]. Also, they reported that the hydrogenation pathway via the R1 complex had more selectivity to the ethylene production rather than ethane production [39]. Therefore, the R1 structure is selected as the initial structure for the hydrogenation mechanism on the Pd 2 /g-C 3 N 4 catalyst. Considering the R1 as a starting complex, two consecutive pathways are proposed for the hydrogenation mechanism. In the first pathway, the acetylene is converted to ethylene (IM3) via the three transition states (TS1, TS2, and TS3) (see Figure 4). In the second pathway, ethylene (IM3) is hydrogenated more and converted to ethane (P) via two transition states (TS4 and TS5). To obtain the activation energy of each step, the E co-adsorption of each complex involve in each step is calculated using Eq. (1). The energy difference between the E co-adsorption of the reactant and TS of each step is considered as the activation energy (Figure 4). Figures 5-9 show the optimized structures of the complexes (R1, IMs, and TSs) involved in the pathways shown in Figure 4. Also, the important interaction paths for each complex, obtained from the QTAIM calculations, are shown in the Figures. The position of the bond critical points (BCPs) (dots) and their calculated electron density ( r ) have been identified in the Figures. R1 ® TS1 ® IM1. The calculated value of the activation energy for the conversion of the R1 complex to IM1 via TS1 is 6.77 kcal/mol (see Figure 5). The optimized structure of TS1 is similar to the R1 complex except that the H 2 bond length increases to 2.77 Å in the structure of TS1 indicating the complete dissociation of the H 2 molecule on the catalyst. The calculated vibrational mode of TS1 with the imaginary vibrational frequency (n=-470 i cm -1 ) is shown in Figure S3. The direction of the displacement vector of one of the H atoms of the dissociated H 2 molecule is towards the C atom of the acetylene indicating the H transfer from the Pd to C atom. Comparison of the arrangement of the BCPs in the R1 complex with that of TS1 shows that the BCP No.1 (BCP1; r =0.1732 a.u), related to the interaction path between two H atoms of H 2 molecule in the R1 complex, disappears. Also, a new BCP (BCP4) is formed between one of the H atoms of the H 2 molecule and the C atom of acetylene in the TS1 complex. As seen, this H atom has simultaneous interaction with the Pd and C atoms via the interaction paths related to the BCP2 ( r =0.1189 a.u) and BCP4 (( r =0.1444 a.u) in the structure of TS1. Due to the interaction of this H atom with the C atom of acetylene, the C-H bond length increases from 1.08 Å in the R1 complex to 1.43 Å in the IM1 complex. These observations confirm the transfer of the H atom from the Pd site to the C atom of acetylene for the formation of the vinyl compound. It is important to note that there is no BCP between two Pd atoms in the structure of the complexes shown in Figure 5 and each Pd atom is interacting with the C atoms of acetylene independent of the other one. In the R1 complex, each of the Pd atoms is interacting with one C atom of acetylene while one of the Pd atoms is in the direct interaction with two C atoms, simultaneously in the TS1 and IM1 complexes via the two BCPs (No. 5 and 6 in TS1 and No. 3 and 4 in IM1). The increase of the value of r for the BCP of Pd-H bond (BCP2, BCP1, and BCP1 in the R1, TS1, and IM1 complexes, respectively) shows that this bond becomes stronger due to the H transfer. Also, the strength of the Pd-N bond between the Pd atom and g-C 3 N 4 becomes weaker in the IM1 (BCP No. 5; r =0.7159 a.u.) complex compared to the R1 complex (BCP No. 5; r =0.7015 a.u.) due to the decrease in its r value. IM1 ® TS2 ® IM2. Before the transfer of the second H atom in the structure of IM1 to the C atom of vinyl for forming IM3 (see Figure 7), a step is necessary (IM1®TS2®IM2). The second H atom is interacting with the Pd atom via the BCP1 in the structure of IM1. Notably, this H transfer becomes easier when the orientation of the complex above the g-C 3 N 4 changes so that this Pd atom can directly interact with the N atoms of the substrate (see Figure 6; IM2). Due to this interaction, the IM2 complex is more stable than the IM1 complex and the orientation of the Pd-H bond is more suitable for the second H transfer to the C atom of vinyl compared to what is observed in the IM1. The bond angle of C-Pd-H in the structure of IM1 is 160.8° and reaches 83.9° in the structure of IM2 so that the distance between the second H atom and the C atom of vinyl decreases (see Figure 6). As seen, the orientation of the complex on the g-C 3 N 4 changes so that the two Pd atoms are in direct interaction with the N atoms of g-C 3 N 4 . The strength of the H-Pd bond does not change due to this configuration change because of the small change of its r value. The calculated value of r of the H-Pd bond in the IM1 (BCP1) is 0.1484 a.u. and the corresponding value in the IM2 is 0.1502 a.u. (BCP5). Moreover, a decrease is seen in the bond length of the formed C-H bond in step 1 in the IM2 complex (2.36 Å) compared to that in the IM1 complex (2.60 Å) which shows that this bond is stronger. Therefore, the transfer of the second H atom from the Pd atom to the C atom in the IM2 needs less activation energy compared to that of this process in the IM1. A considerable decrease in the r values of the C-Pd interaction paths of the TS2 structure compared to those of the C-Pd interaction paths of the IM1 structure is seen. For example, the r of the BCPs No.3 and 4 (0.9617 and 0.9885 a.u, respectively), related to the simultaneous interaction of one Pd atom with two C atoms in the structure of IM1, reduces to 0.3508 and 0.4902 a.u., respectively in the structure of TS2 showing the decrease of the strength of these bonds (see Figure 6). Also, the arrangement of the interaction paths in the IM2 shows that each Pd atom is mainly interacting with one C atom unlike what is observed for the IM1 and TS2 structures in which one of the Pd atoms is in the direct interaction with two C atoms, simultaneously. One important point related to the TS2 structure is that the orientation of the complex on the g-C 3 N 4 is so that its interaction with the substrate takes place mainly via the H atoms (see BCPs No.4, 7, and 8 in Figure 6). The imaginary vibrational mode of the TS2 is shown in Figure S3. As seen, the H atom which is interacting with the g-C 3 N 4 via the BCPs No. 7 and 8, moves out of the substrate plane. IM2 ® TS3 ® IM3. Figure 7 shows the transfer of the second H atom from the Pd atom to the C atom of vinyl via the TS3 structure and form of the IM3. The comparison of the structure of TS3 with IM2 shows that the bond angle of C-Pd-H related to BCPs No. 2 and 5 in IM2 decreases in the TS3 structure (see BCPs No. 2 and 4). Comparison of the value of r of the BCP5 in the IM2 (0.1502 a.u.) with that of the corresponding BCP in the structure of TS3 (BCP4 ; r = 0.1341 a.u. ) shows that the decrease of the strength of the Pd-H bond due to the decrease of the C-Pd-H bond angle. The imaginary vibrational mode of the TS3 is shown in Figure S4. The displacement vector of the H atom shows its transfer from the Pd atom to the C atom. Notably, a BCP ( r =0.5594 a.u.) is seen between two Pd atoms in the structure of IM3 showing the interaction of two Pd atoms with each other which is absent in the structure of IM2 and TS3. IM4 ® TS4 ® IM5. Figure 8 shows the first step of the hydrogenation of ethylene to form ethane via TS4. A new H 2 molecule is added to the IM3 for forming the IM4 so that the H 2 molecule is located near the Pd atoms to feasible the breaking of the H 2 bond. As seen in the figure, the dissociation of the H 2 molecule takes place in the TS4 structure so that one of the H atoms of H 2 is connected to one of the Pd atoms via the BCP2, and the other one is pulled toward the carbon atom of ethylene so that it is in the direct interaction with three atoms via the three BCPs (BCP4, BCP5, and BCP3), simultaneously. The imaginary vibrational frequency of the TS4 is shown in Figure S4. As seen, the displacement vector of one of the H atoms shows its tendency for forming the bond with the C atom of ethylene. The imaginary frequency of this vibrational mode and the barrier energy for this hydrogen transfer is 409 i cm -1 and 12.28 kcal/mol, respectively. The comparison of the structure of IM4 with TS4 shows the decrease of the Pd-Pd distance so that a new BCP appears between two Pd atoms in the IM4 structure. The strength of the Pd-Pd interaction in IM5 is lower than that in IM4 because of the decrease of the r of BCP3 from IM4 to IM5. Comparison of the r value of C-H bond in TS4 (BCP4 ; 0.1557 a.u.) with the corresponding value in IM5 (BCP1 ; 0.2384 a.u.) shows the complete formation of the new C-H bond in IM5 compared to TS4. IM5 ® TS5 ® P. Figure 9 shows the last step of the hydrogenation mechanism related to the second H transfer from the Pd atom to the C atom of CH 3 CH 2 via the TS5 structure to form ethane. As seen, this H atom interacts with the Pd atom via the BCP5 ( r = 0.1960 a.u.) in the structure of IM5. Also, the increase of the number of BCPs between the Pd atom and g-C 3 N 4 is seen in the TS5 structure (BCPs No. 3 and 4) compared to the IM5 structure (BCP No. 6) which helps to the hydrogen transfer by the weakening of the Pd-H bond as well as the presence of the BCP1 between the H and the C atom in this structure. Comparison of the calculated value of r of the BCP2 (0.1484 a.u.) in the structure of TS5 with that of BCP5 (0.1960 a.u.) in the structure of IM5 shows the decrease of the strength of the Pd-H interaction and confirms the H transfer from the Pd atom to C atom. The calculated value of vibrational frequency of the imaginary vibrational mode of TS5 (see Figure S4) related to the H transfer is 535 i cm -1 which is in agreement with the literature [39]. After the transfer of the last H atom to the CH 3 CH 2 and the formation of ethane (structure P in Figure 9), one of the Pd atoms becomes close to the substrate so which creates a folding in the structure of g-C 3 N 4 . A comparison has been made between the imaginary vibrational frequencies of the TSs of the hydrogenation of acetylene using the Pd 2 /g-C 3 N 4 catalyst with that reported in the literature [39] for Pd 1 /g-C 3 N 4 (see Table 2). As seen, the calculated vibrational frequencies for the R1®TS1®IM1, IM4®TS4®IM5, and IM5®TS5®P steps of our work are smaller than the value reported by Zhao et al. [39]. The value of the imaginary vibrational frequency of the TS of IM1®TS2®IM2 is nearly equal to that reported by Zhao et al. for the Pd 1 /g-C 3 N 4 catalyst. The imaginary frequency of the transition state of IM2®TS3®IM3 is greater than that reported for the Pd 1 /g-C 3 N 4 catalyst. Selectivity and reactivity As previously mentioned, the effective performance of a heterogeneous catalyst in the hydrogenation of acetylene is the production of ethylene with a smaller value of activation energy compared to that of ethane production [39-40]. According to Figure 4, the activation energy for the first H transfer reaction (R1-TS1-IM1) is equal to 6.77 kcal/mol and it is smaller than the values reported in the literature for the other heterogeneous catalyst (see Table 3). Also, the comparison of the activation energy of the R1-TS1-IM1 step obtained in this work using the Pd 2 /g-C 3 N 4 catalyst (6.77 kcal/mol) with the corresponding value reported by Zhao et al. [39] using the Pd 1 /g-C 3 N 4 catalyst (21.53 kcal/mol) shows that the increase of the number of Pd atom has a considerable effect on the decrease of activation energy of this step. The selectivity of a heterogeneous catalyst in the acetylene hydrogenation, for the production of ethylene, increases by the decrease of the activation energy of the ethylene production step compared to that of the ethane production step. The selectivity (D E ) is considered as the difference between the hydrogenation barrier ( E hydrogenation ) and desorption barrier ( E desorption ) of ethylene [39-40]. The higher value of ∆ E leads to the higher selectivity of the catalyst for ethylene production. The calculated value of ∆ E obtained in this work is 21.9 kcal/mol (see Table3) which is more than that reported by Zhao et al. for the Pd 1 /g-C 3 N 4 catalyst. This shows that the increase in the number of Pd atoms of the catalyst is accompanied by the increase of the selectivity of the catalyst for ethylene production. Conclusion In this work, a theoretical mechanistic study of the hydrogenation of acetylene on the Pd 2 /g-C 3 N 4 heterogeneous catalyst was performed using the DFT and QTAIM calculations. The selectivity of this catalyst for the ethylene production was compared with those of the other catalysts reported in the literature (Pd(111) [45], Al 13 CoO 4 [46], AlPd [47], Pd1/g-C 3 N 4 [39]). It was observed that the Pd 2 /g-C 3 N 4 has more selectivity and activity for ethylene production compared to the mentioned catalyst. Also, the results of this work show that the clustering of the Pd atoms on the g-C 3 N 4 increases the selectivity and activity of the catalyst for the acetylene hydrogenation. This work provides strong insights for further research on the role of different palladium clusters deposited on graphitic carbon nitride material in acetylene hydrogenation. Declarations Acknowledgment The authors thank Isfahan University of Technology (IUT) for its supports. The authors gratefully acknowledge the Shaikh Bahaie National High-Performance Computing Center (SBNHPCC) for providing a computing Facility. Funding The authors thank Isfahan University of Technology (IUT) for its supports. The contribution of Authors All authors conceived and designed the calculations; analyzed and interpreted the data; contributed materials, analysis tools, or data and software; and wrote the paper. Conflict of interest The authors declare that they have no conflict of interest Availability of data and material All data generated or analyzed during this study are included in this published article [and its supplementary information files]. Code availability N/A References [1] Pradier C M, Mazina M, Berthier Y, Oudar J (1994) Hydrogenation of acetylene on palladium. J Mol Catal 89: 211-220. [2] Sarkany A, Geszti O, Safran G (2008) Preparation of Pd shell–Au core/SiO 2 catalyst and catalytic activity for acetylene hydrogenation. Appl Catal A: Gen 350: 157–163. [3] Wang P, Yang B Influence of surface strain on activity and selectivity of Pd‐based catalysts for the hydrogenation of acetylene: A DFT study (2018) Chinese J Catal 39: 1493–1499. [4] Lopez N, Vargas-Fuentes C (2012) Promoters in the hydrogenation of alkynes in mixtures: insights from density functional theory. Chem Commun 48: 1379–1391. [5] Studt F, Abild‐Pedersen F, Bligaard T, Sørensen R Z, Christensen C H, Norskov J K (2008) On the role of surface modifications of palladium catalysts in the selective hydrogenation of acetylene. Angew Chem Int Ed 47: 9299–9302. [6] Vile G, Albani D, Nachtegaal M, Chen Z, Dontsova D, Antonietti M (2015) A stable single-site palladium catalyst for hydrogenations. Angew Chem 54: 11265. [7] Huang X, Xia Y, Cao Y, Zheng X, Pan H, Zhu J (2017) Enhancing both selectivity and coking-resistance of a single-atom Pd 1 /C 3 N 4 catalyst for acetylene hydrogenation. Nano Res 10: 1302–1312. [8] Kyriakou G, Boucher M B, Jewell A D, Lewis E A, Lawton T J, Baber A E (2012) Isolated metal atom geometries as a strategy for selective heterogeneous hydrogenations. Science 2012 335: 1209. [9] Liu J, Bunes B R, Zang L, Wang C (2013) Supported single-atom catalysts: synthesis, characterization, properties, and applications. Environ Chem Lett 2017: 1–29. [10] Kang L, Cheng B, Zhu M (2019) Pd/MCM-41 catalyst for acetylene hydrogenation to ethylene. R Soc Open Sci 6: 191155-191166. [11] Panpranot J, Kontapakdee K, Praserthdam P (2006) Selective hydrogenation of acetylene in excess ethylene on micron-sized and nanocrystalline TiO 2 supported Pd catalysts. Appl Cat A Gen 314: 128–133. [12] Panpranot J, Kontapakdee K, Praserthdam P (2006) Effect of TiO 2 crystalline phase composition on the physicochemical and catalytic properties of Pd/TiO 2 in selective acetylene hydrogenation. J Phys Chem B 110: 8019–8024. [13] Meng L D, Wang G C (2014) A DFT + U study of acetylene selective hydrogenation over anatase supported Pd a Ag b (a + b = 4) cluster. Phys Chem Chem Phys 16: 17541-17551. [14] Yang J, Cao L X, Wang G C (2012) Acetylene hydrogenation on anatase TiO 2 (101) supported Pd 4 cluster: oxygen deficiency effect. J Mol Model 18: 3329–3339. [15] Ma L L, Lv C Q, Wang G C (2017) A DFT study and micro-kinetic analysis of acetylene selective hydrogenation on Pd-doped Cu (111) surfaces. Appl Surf Sci 410: 154–165. [16] Yang J, Lv C Q, Guo Y, Wang G C (2012) A DFT+U study of acetylene selective hydrogenation on oxygen defective anatase (101) and rutile (110) TiO 2 supported Pd 4 cluster. J Chem Phys 136: 104107-104122. [17] Li J N, Pu M, Ma C C, Tian Y, He J, Evans D J (2012) The effect of palladium clusters (Pdn, n = 2–8) on mechanisms of acetylene hydrogenation: A DFT study. J Mol Catal A: Chem 359: 14–20. [18] Boyukata M, Belchior J C (2008) Molecular dynamics study of palladium clusters: size dependent analysis of structural stabilities and energetics of Pdn (n<40) via a lennard-jones type potential. Croatica Chem Act 81: 289-297. [19] Moc J, Musaev D J, Morukama D (2003) Activation and adsorption of multiple H 2 molecules on a Pd 5 cluster: a density functional study. J Phys Chem A 107: 4929-4939. [20] Zhou C, Yao S, Wu J, Forrey R C, Chen. L, Tachibana A, Chen H (2008) Hydrogen dissociative chemisorption and desorption on saturated sub nano palladium clusters (Pdn, n = 2–9). Phys Chem Chem Phys 10: 5445–5451. [21] Cao Y, Sui Z, Zhu Y, Zhou X, Chen D (2017) Selective hydrogenation of acetylene over Pd-In/Al 2 O 3 catalyst: promotional effect of indium and composition-dependent performance. ACS Catal 7: 7835–7846. [22] Albani D, Shahrokhi M, Chen Z, Mitchell S, Hauert R, Lopez N, Pérez-Ramírez J (2018) Selective ensembles in supported palladium sulfide nanoparticles for alkyne semi-hydrogenation. J Nat Commun 9: 2634-2644. [23] Lopez N, Vargas-Fuentes C (2012) Promoters in the hydrogenation of alkynes in mixtures: insights from density functional theory. Chem Commun 48: 1379–1391. [24] Volpe M A, Rodriguez P, Gigola C E (1999) Preparation of Pd–Pb/α-Al 2 O 3 catalysts for selective hydrogenation using PbBu 4 : the role of metal-support boundary atoms and the formation of a stable surface complex. Catal Lett 61: 27–32. [25] Anderson J A, Mellor J, Wells R P K (2009) Pd catalyzed hexyne hydrogenation modified by Bi and by Pb. J Catal 261: 208–216. [26] Kim S K, Lee L H, Ahn I K, Kim W J, Moon S H (2011) Performance of Cu-promoted Pd catalysts prepared by adding Cu using a surface redox method in acetylene hydrogenation. ApplCatalA Gen 401: 12–19. [27] Gonzalez S, Neyman K M, Shaikhutdinov S, Freund H J, Illas F (2007) On the promoting role of Ag in selective hydrogenation reactions over Pd–Ag bimetallic catalysts: a theoretical study. J Phys Chem C 111: 6852–6856. [28] Mei D, Neurock M, Smith C M (2009) Hydrogenation of acetylene–ethylene mixtures over Pd and Pd–Ag alloys: first-principle based kinetic Monte Carlo simulations. J Catal 268: 181–195. [29] Sheth P A, Neurock M, Smith C M (2009) First-principles analysis of the effects of alloying Pd with Ag for the catalytic hydrogenation of acetylene–ethylene mixtures. J Phys Chem B 109: 12449–12466. [30] Kumar N, Ghosh P (2016) Selectivity and reactivity of Pd-rich Pd-Ga surfaces toward selective hydrogenation of acetylene: interplay of surface roughness and ensemble effect. J Phys Chem C 120: 28654–28663. [31] Osswald J, Giedigkeit R, Jentoft R E, Armbrüster M, Girgsdies F, Kovnir K (2008) Palladium–gallium intermetallic compounds for the selective hydrogenation of acetylene. Part I. Preparation and structural investigation under reaction conditions. J Catal 258: 210–218. [32] Zhou H, Yang X, Wang A, Miao S, Liu X, Pan X (2016) Pd/ZnO catalysts with different origins for high chemoselectivity in acetylene semi-hydrogenation. Chin J Catal 37: 692–699. [33] Gulyaeva Y K, Kaichev V V, Zaikovskii V I, Kovalyov E V, Suknev A P, Balzhinimaev B S (2015) Selective hydrogenation of acetylene over novel Pd/fiberglass catalysts. Catal Today 245: 139–146. [34] Yan H, Cheng H, Yi H, Lin Y, Yao T, Wang T (2015) Single-atom Pd 1 /Graphene catalyst achieved by atomic layer deposition: remarkable performance in selective hydrogenation of 1,3-butadiene. J Am Chem Soc 137: 10484–10487. [35] Komhom S, Mekasuwandumrong O, Praserthdam P, Panpranot J (2008) Improvement of Pd/Al 2 O 3 catalyst performance in selective acetylene hydrogenation using mixed phases Al 2 O 3 support. Catal Commun 10: 86–91. [36] Hosseini S M, Ghiaci M, Kulinich S A, Wunderlich W, Farrokhpour H, Saraji M, Shahvar A (2018) Au-Pd@g-C 3 N 4 as an efficient photocatalyst for visible-light oxidation of benzene to phenol: experimental and mechanistic study. J Phys Chem C 122: 27477–27485. [37] Shahzeydi A, Ghiaci M, Farrokhpour H, Shahvar A, Sun M, Saraji M (2019) Facile and green synthesis of copper nanoparticles loaded on the amorphous carbon nitride for the oxidation of cyclohexane. Chem Eng J 370: 1310-1321. [38] Nilforoushan S, Ghiaci M, Hosseini S M, Laurent S, Muller R N (2019) Selective liquid phase oxidation of ethyl benzene to acetophenone by palladium nanoparticles immobilized on a g-C 3 N 4 –rGO composite as a recyclable catalyst. New J Chem 46: 6921-6931. [39] Zhao Y, Zhu M, Kang L (2018) The DFT study of single-atom Pd 1 /g-C 3 N 4 catalyst for selective acetylene hydrogenation reaction. Cattal Lett 148: 2992–3002. [40] Kang L, Zhu M, Zhao Y (2019) A DFT study of acetylene hydrogenation catalyzed by S-doped Pd 1 /g-C 3 N 4 . Catalysts 9: 887-899. [41] Hosseini S M, Ghiaci M, Farrokhpour H (2019) The adsorption of small size Pd clusters on a g-C 3 N 4 quantum dot: DFT and TD-DFT study. Mater Res Express 6: 105079-105093. [42] Singh N B, Sarkar U (2014) Structure, vibrational, and optical properties of platinum cluster: a density functional theory approach. J Mol Model 20: 2537–2548. [43] Lu T, Chen F (2012) Multiwfn: a multifunctional wavefunction analyzer. J Comput Chem 33: 580–592. [44] Frisch M J (2013) Gaussian 09, Revision D.01, Wallingford, CT: Gaussian, Inc. [45] Mei D, Sheth P A, Neurock M, Smith C M (2006) First-principles based kinetic Monte Carlo simulation of the selective hydrogenation of acetylene over Pd (111). J Catal 242: 1–15. [46] Krajci M, Hafner J (2011) Complex intermetallic compounds as selective hydrogenation catalysts—a case study for the (100) surface of Al 13 Co 4 . J Catal 278: 200–207. [47] Krajci M, Hafner J (2012) Intermetallic compound AlPd as a selective hydrogenation catalyst: a DFT study. J Phys Chem C 116: 6307–6319. Table Table 1. The energies of the HOMO and LUMO of C 2 H 2 , H 2 , g-C 3 N 4 , Pd 2 /g-C 3 N 4 , and Pd 2 /g-C 3 N 4 ·C 2 H 2 compounds and their corresponding energy gaps HOMO LUMO ∆E g H 2 -11.74 (eV) 1.84 -9.89 C 2 H 2 -8.25 0.45 -7.79 g-C 3 N 4 -6.53 -2.62 -3.90 Pd 2 /g-C 3 N 4 -4.71 -2.71 -2.00 Pd 2 /g-C 3 N 4 ·C 2 H 2 -4.09 -2.67 -1.41 Table 2. The imaginary frequencies of the transition states (TSs) involved in the different steps of the acetylene hydrogenation mechanism taken from Zhao’s group work [39] for the Pd1/g-C 3 N 4 with those obtained in this work. Imaginary Frequency (cm -1 ) Catalyst IM5-TS5-P IM4-TS4-IM5 IM2-TS3-IM3 IM1-TS2-IM2 R1-TS1-IM1 760.9 i 685.08 i 641.95 i 724.3 i 750.62 i Pd 1 /g-C 3 N 4 535 i 409 i 711 i 723.1 i 476 i Pd 2 /g-C 3 N 4 Table 3. Comparison of the ctivation energies for the first hydrogen transfer to the acetylene and ethylene reported in literature with those obtained in this work. Reference Selectivity (∆E) Energy barrier (kcal/mol) Catalyst C 2 H 4 + H 2 → C 2 H 6 C 2 H 2 + H 2 → C 2 H 4 45 2.45 17.61 18.1 Pd(111) 46 6.34 19.57 15.41 Al 13 CoO 4 47 11.25 13.7 18.59 AlPd 39 15.4 38.88 21.53 Pd 1 /g-C 3 N 4 This Study 21.9 12.28 6.77 Pd 2 /g-C 3 N 4 Supplementary Files SupportingInformation22.docx Cite Share Download PDF Status: Published Journal Publication published 19 Apr, 2021 Read the published version in Structural Chemistry → Version 1 posted Reviewers invited by journal 03 Mar, 2021 Editor assigned by journal 02 Mar, 2021 First submitted to journal 19 Feb, 2021 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-259183","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":15011385,"identity":"5ecfc081-b551-4811-8f7e-3e89c3df23de","order_by":0,"name":"Sayyed Mahdi Hosseini","email":"","orcid":"","institution":"Isfahan University of Technology","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Sayyed","middleName":"Mahdi","lastName":"Hosseini","suffix":""},{"id":15011386,"identity":"f9f0df94-efe9-4271-8fa2-d91480ee1953","order_by":1,"name":"Mehran Ghiaci","email":"","orcid":"","institution":"Isfahan University of Technology","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Mehran","middleName":"","lastName":"Ghiaci","suffix":""},{"id":15011387,"identity":"f0f50ff4-acba-4348-bd82-e4cd20ad0e04","order_by":2,"name":"Hossein Farrokhpour","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABBklEQVRIiWNgGAWjYBAC9gYQeeAAgwGYa8AgB+Y/wKOF5wCaFmMwP4F4LQwMiWBr8WphP/vsw4czd/LM2dsff/hQcCd9ftjhh0Bb7OR0G3Bo4Uk3njnjxrNiy54zZpIzDJ7lbrydZgDUkmxsdgC7FnuGNGZmng+HEzfcyGFj5jE4nLtxdgJIy4HEbTi08PA/g2lJf/z5j8HhdMPZ6R/wa5EA2XIDpCXBQJrB4HCCvHQOAVsknjEzzjhzuNjgDNAvPQaHDTdI5xQcSDDA7Rce/jRmhg/HDucZHAeG2I8/h+XlZ6dv/vChwk4OlxYYSICzDMAqDXCow6pFvoGw6lEwCkbBKBhZAAAXFGuHg1ogPAAAAABJRU5ErkJggg==","orcid":"","institution":"Chemistry Department-Isfahan University of Technology","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Hossein","middleName":"","lastName":"Farrokhpour","suffix":""}],"badges":[],"createdAt":"2021-02-20 16:09:55","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-259183/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-259183/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s11224-021-01781-3","type":"published","date":"2021-04-19T19:03:36+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":6824289,"identity":"3ade4601-0c0f-44b9-b849-e52242d9225e","added_by":"auto","created_at":"2021-03-11 01:14:10","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":120685,"visible":true,"origin":"","legend":"The optimized structures of the g-C3N4 quantum dot and the Pd2/g-C3N4 catalyst shown in two different views. ","description":"","filename":"1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-259183/v1/f8b22f591d9c207ed85c3646.jpg"},{"id":6824291,"identity":"cc3f7c73-d79c-4be6-a9c4-78ac42373b11","added_by":"auto","created_at":"2021-03-11 01:14:11","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":79788,"visible":true,"origin":"","legend":"Mulliken partial atomic charges distribution of (a) g-C3N4 quantum dot and (b) Pd2/g-C3N4 complex.","description":"","filename":"2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-259183/v1/753365f2abe6c8d66894b967.jpg"},{"id":6824292,"identity":"562623c1-f194-4327-a21c-1760a4bcf3e6","added_by":"auto","created_at":"2021-03-11 01:14:11","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":118977,"visible":true,"origin":"","legend":"Frontier molecular orbitals of the g-C3N4, Pd2/g-C3N4, and C2H2.Pd2/g-C3N4 structures ","description":"","filename":"3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-259183/v1/4f0ad2daf0edc752cc4aae4b.jpg"},{"id":6823830,"identity":"6126478c-74fe-4189-a8c0-640a0434a228","added_by":"auto","created_at":"2021-03-11 01:11:10","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":47829,"visible":true,"origin":"","legend":"Energy diagram for the acetylene hydrogenation over the Pd2/g-C3N4 catalyst","description":"","filename":"4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-259183/v1/6024aa662b126b368aeed7a6.jpg"},{"id":6824290,"identity":"93e3f5e3-d9fe-47eb-8a90-93b57761265a","added_by":"auto","created_at":"2021-03-11 01:14:10","extension":"jpg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":86547,"visible":true,"origin":"","legend":"The optimized structure of the R1, TS1, and IM1 of the first step of hydrogenation of acetylene on the Pd2/g-C3N4 catalyst. The red points show the position of the bond critical points (BCPs) obtained from the QTAIM calculations. The value of electron density () for each BCP has been given.","description":"","filename":"5.jpg","url":"https://assets-eu.researchsquare.com/files/rs-259183/v1/c62687889ff80bbad2ab9c9e.jpg"},{"id":6823837,"identity":"33e17eb6-655b-40b3-bd57-fb6af22362e9","added_by":"auto","created_at":"2021-03-11 01:11:11","extension":"jpg","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":83800,"visible":true,"origin":"","legend":"The optimized structure of the IM1, TS2, and IM2 of the second step of hydrogenation of acetylene on the Pd2/g-C3N4 catalyst. The red points show the position of the bond critical points (BCPs) obtained from the QTAIM calculations. The value of electron density () for each BCP has been given.","description":"","filename":"6.jpg","url":"https://assets-eu.researchsquare.com/files/rs-259183/v1/429fd735292a6eb820fe2698.jpg"},{"id":6823834,"identity":"be6aeedb-c999-4a27-8caa-46f4feb17d00","added_by":"auto","created_at":"2021-03-11 01:11:11","extension":"jpg","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":91776,"visible":true,"origin":"","legend":"The optimized structure of the IM2, TS3, and IM3 of the third step of hydrogenation of acetylene on the Pd2/g-C3N4 catalyst. The red points show the position of the bond critical points (BCPs) obtained from the QTAIM calculations. The value of electron density () for each BCP has been given.","description":"","filename":"7.jpg","url":"https://assets-eu.researchsquare.com/files/rs-259183/v1/8b416b24b5839a37ac30ef97.jpg"},{"id":6824293,"identity":"a5a2817c-10f3-45a5-9120-2ab9c0002dc4","added_by":"auto","created_at":"2021-03-11 01:14:11","extension":"jpg","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":83565,"visible":true,"origin":"","legend":"The optimized structure of the IM4, TS4, and IM5 of the first step of hydrogenation of ethylene on the Pd2/g-C3N4 catalyst. The red points show the position of the bond critical points (BCPs) obtained from the QTAIM calculations. The value of electron density () for each BCP has been given.","description":"","filename":"8.jpg","url":"https://assets-eu.researchsquare.com/files/rs-259183/v1/4a5bdf45734ee1fb2faadf4c.jpg"},{"id":6823836,"identity":"297234d5-fe07-49aa-8f7c-cb159bd89b45","added_by":"auto","created_at":"2021-03-11 01:11:11","extension":"jpg","order_by":9,"title":"Figure 9","display":"","copyAsset":false,"role":"figure","size":73714,"visible":true,"origin":"","legend":"The optimized structure of the IM5, TS5, and P of the second step of hydrogenation of ethylene on the Pd2/g-C3N4 catalyst. The red points show the position of the bond critical points (BCPs) obtained from the QTAIM calculations. The value of electron density () for each BCP has been given.","description":"","filename":"9.jpg","url":"https://assets-eu.researchsquare.com/files/rs-259183/v1/96f8f7e66f56dd4456a2161f.jpg"},{"id":13677276,"identity":"2bcc9d56-00f5-4fcc-8839-37d2151a516a","added_by":"auto","created_at":"2021-09-17 11:34:27","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1058996,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-259183/v1/6d07b411-4c12-47ae-b93c-8ef6e2800cf0.pdf"},{"id":6823839,"identity":"f4126c68-f760-40e7-8e47-4d54238c807f","added_by":"auto","created_at":"2021-03-11 01:11:11","extension":"docx","order_by":7,"title":"","display":"","copyAsset":false,"role":"supplement","size":516945,"visible":true,"origin":"","legend":"","description":"","filename":"SupportingInformation22.docx","url":"https://assets-eu.researchsquare.com/files/rs-259183/v1/73201c645741dd6292f0f7ca.docx"}],"financialInterests":"","formattedTitle":"\u003cp\u003eMechanistic Insight into the Hydrogenation of Acetylene on the Pd2/g-C3N4 Catalyst: Effect of Pd Clustering on the Barrier Energy and Selectivity\u003c/p\u003e","fulltext":[{"header":"Introduction","content":"\u003cp\u003eEthylene is an important raw material with wide applications in petrochemical industries [1]. Generally, the production of ethylene is performed via the steam cracking of hydrocarbons. The main disadvantage of this method is the catalyst poisoning leading to the reduction of catalyst lifetime.[2]. So, the catalytic hydrogenation of acetylene followed by the classical Horiuti\u0026ndash;Polanyi mechanism has been considered as the best method for ethylene production on large scales [3]. An efficient catalyst for acetylene hydrogenation should have high activity and selectivity for ethylene production [4]. In this regard, the noble transition metals, especially, palladium (Pd)-based catalysts are widely used for this purpose. Because they show good selectivity for ethylene production because of their superior hydrogen solubility and good ability for dissociation of H\u003csub\u003e2\u003c/sub\u003e molecules [5,6]. There are several experimental works aiming production of ethylene via acetylene hydrogenation in a selective manner in which the structure of the Pd-based catalyst plays a critical role in achieving a good selectivity [7-12].\u003c/p\u003e\n\u003cp\u003eHuang et al. prepared the dispersed Pd atoms on the g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e using atomic layer deposition and performed the hydrogenation of acetylene on the Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecatalyst [7]. Boucher et al. performed the hydrogenation of styrene and acetylene on the catalyst composed of the Pd atoms binding to the Cu surface [8]. Panpranot et al. investigated the catalytic efficiency of Pd catalyst supported on nanocrystalline TiO\u003csub\u003e2\u003c/sub\u003e (Pd/TiO\u003csub\u003e2\u003c/sub\u003e) for the selective hydrogenation of acetylene in the presence of excess ethylene [11]. The same group also examined the effect of the phase composition of the TiO\u003csub\u003e2\u003c/sub\u003e crystalline support on the physicochemical and catalytic properties of the Pd/TiO\u003csub\u003e2\u003c/sub\u003e catalyst on the selective hydrogenation of acetylene [12].\u003c/p\u003e\n\u003cp\u003eThere are several theoretical studies of the clarification of the hydrogenation mechanism of acetylene on the Pd-based catalysts, and the interaction between the different forms of Pd and the substrates [13-16]. Meng et al. studied the hydrogenation of acetylene on the TiO\u003csub\u003e2\u003c/sub\u003e supported Pd\u003cem\u003e\u003csub\u003ea\u003c/sub\u003e\u003c/em\u003eAg\u003cem\u003e\u003csub\u003eb\u003c/sub\u003e\u003c/em\u003e (\u003cem\u003ea\u003c/em\u003e+\u003cem\u003eb\u003c/em\u003e=4) catalyst using the density functional theory (DFT)+U method [13]. The hydrogenation of acetylene using the Pd\u003csub\u003e4\u003c/sub\u003e cluster on the anatase TiO\u003csub\u003e2 \u003c/sub\u003e(101) support was studied by Yang et al. [14]. Ma et al. performed the DFT study of the selective hydrogenation of acetylene on the P-doped Cu (111) surface [15]. The effect of the size of Pd cluster (Pd\u003cem\u003e\u003csub\u003en\u003c/sub\u003e\u003c/em\u003e) and its geometry on the hydrogenation of acetylene has been investigated [17-20]. Li et al. investigated the effect of the size of Pd clusters on acetylene hydrogenation in the absence of any support [17]. Pu et al. investigated the effect of Pd\u003cem\u003e\u003csub\u003en\u003c/sub\u003e\u003c/em\u003e clusters with different sizes (n=2-8) on acetylene hydrogenation by the DFT calculations. They proposed two pathways for acetylene hydrogenation including (\u003cem\u003ea\u003c/em\u003e) ethene formation as the main product via vinyl intermediate and (\u003cem\u003eb\u003c/em\u003e) the vinylidene, ethylidene, and ethylidene as the intermediates participating in the hydrogenation process. Also, They showed that the reaction pathway depends strongly on the size of Pd\u003cem\u003e\u003csub\u003en\u003c/sub\u003e\u003c/em\u003e clusters so that for \u003cem\u003en\u003c/em\u003e\u0026lt;4, the hydrogenation reaction proceeded through the vinyl intermediate while the other pathway via the vinylidene intermediate was the main pathway for the larger Pd clusters [17].\u003c/p\u003e\n\u003cp\u003eDue to the high activity of Pd catalysts, the over-hydrogenation of ethylene is always a deep concern. Therefore, designing a selective catalyst for ethylene production to eliminate extra purification steps in its synthesis and prevent ethane production is considered a challenging task in acetylene hydrogenation [21-22]. Various strategies have been so far reported for the promotion of the performance of Pd-based catalysts for the acetylene hydrogenation which one of the important of them is the doping of the Pd cluster with other noble metals such as Au [23], Pb [24,25], Cu [26], Ag [27,28,29], and Ga [30,31]. The other important strategy is using single-atom alloys (SAA) and single-atom catalysts (SAC), immobilized on solid supports such as ZnO [32], fiberglass [33], graphene [34], Al\u003csub\u003e2\u003c/sub\u003eO\u003csub\u003e3 \u003c/sub\u003e[35], and graphitic carbon nitride (g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e) [6-7] which has shown interesting results in acetylene hydrogenation. Among the solid supports, g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003eis always considered a suitable candidate for the deposition of Pd active sites [36-38]. It has some outstanding features such as low-cost fabrication, short time and easy synthesis, adjustable surface area, good dispersion potential for active sites, and high coordination ability for metal clusters and nanoparticles. Huang et.al showed that Pd single-atom immobilized on the surface of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e has a better selectivity in acetylene hydrogenation compared to nanoparticle ones [7]. Zhao et al.\u0026rsquo;s theoretical study aiming to demonstrate the acetylene hydrogenation mechanism over the Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst showed that the Pd atom tended to adsorb on the top of the sixth fold cavity of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e as the most accessible site. Moreover, the Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst showed better selectivity to ethylene production in comparison to ethane production as expected from the experimental studies [39]. Kang et al. investigated the catalytic activity of the Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst for acetylene hydrogenation by replacing the nitrogen and carbon atoms of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ewith the sulfur atom (S) [40]. It was observed that the Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003econtaining sulfur-doped carbon atoms have higher activity for the dissociation of H\u003csub\u003e2\u003c/sub\u003e molecules while better selectivity is achieved for sulfur-doped nitrogen ones.\u003c/p\u003e\n\u003cp\u003eTo the best of our knowledge, there is no theoretical mechanistic study of acetylene hydrogenation over Pd\u003csub\u003en\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e structures with n\u0026gt;1 in the literature. Our recent theoretical study on the adsorption of small-size Pd\u003csub\u003en\u003c/sub\u003e clusters (n=2, 4, 6, and 8), immobilized on the g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e quantum dot, reveals that the Pd\u003csub\u003e2\u003c/sub\u003e cluster has the highest adsorption energy among the other selected Pd clusters [41]. In this work, the DFT calculations are performed to understand the exact mechanism of acetylene hydrogenation on the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst and see the effect of the added Pd atom on the activity of catalyst compared to that of Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst reported in references 39 and 40. Our calculated results show less activation energy for acetylene hydrogenation using the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e compared to the Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecatalyst and better selectivity for ethylene production which is in good accordance with the experimental data in the literature [39,45-47].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eComputational details\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe structure of all reactants (Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e, C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e, and H\u003csub\u003e2\u003c/sub\u003e), pre-reactant (R), transition states (TSs), and intermediates (IMs) were optimized in the gas phase using the DFT method employing B3PW91 functional [42]. The 6-311+G(d,p) and Los Alamos effective core pseudo-potential (ECP) (LANL2DZ) basis set were used for the g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e and Pd\u003csub\u003e2\u003c/sub\u003e, respectively. The frequency calculations were performed on all optimized structures to check their correct positions on their potential energy surfaces in the hydrogenation pathway and determine their electronic energies, enthalpies, Gibbs free energies, and zero-point energies. The suitability of the selected DFT functional and basis set for the considered systems have been investigated by Zhao et al [39]. The intrinsic reaction coordinate (IRC) calculations were performed to make sure that each TS structure was related to its corresponding reactants and products in the reaction mechanism. The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of the structures were calculated to obtain the bandgap energies and donor-acceptor properties of structures. The Quantum Theory of Atoms in Molecules (QTAIM) calculations were performed using Multiwfn 3.4 software [43] for better structural characterization of the TSs and confirming the proposed mechanism of the reaction. All DFT calculations were performed using the Gaussian 09 Quantum Chemistry Package [44].\u003c/p\u003e"},{"header":"Results And Discussion","content":"\u003cp\u003e\u003cstrong\u003eThe optimized structure of the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecatalyst and its charge distribution\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eA quantum dot of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e is selected as a model of substrate for the deposition of the Pd\u003csub\u003e2\u003c/sub\u003e cluster. The g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e has two main structural forms based on the s-triazine and tri-s-triazine building blocks so that the latter form is always considered as the most stable form of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003e[39-41]. A layer of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e composed of three tri-s-triazine units was selected and optimized (see Figure 1). The nitrogen atoms in this model are classified into two classes: (\u003cem\u003ea\u003c/em\u003e) saturated nitrogen atoms which are bonded to three carbon atoms (N1) and (\u003cem\u003eb\u003c/em\u003e) pyridine-like nitrogen atoms which are bonded to two carbon atoms and are located at the corners of the sixth fold cavity (N2). Since 2p orbitals of N2 atoms have the main contribution in the valance band of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e, their lone pairs electrons have the main role in the chemical activity of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecompared to the N1 atoms [41]. On the other hand, 2p orbitals of the N1 and C atoms orbitals have a major contribution to the conduction band of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e [41].\u003c/p\u003e\n\u003cp\u003eThe Pd\u003csub\u003e2\u003c/sub\u003e cluster is optimized, separately and placed on the top of the sixth fold cavity of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003eaiming to achieve the maximum interaction with the lone pairs of N2 atoms. The comparison of the structure of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e in the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecomplex with the optimized structure of isolated g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e shows that the presence of the Pd\u003csub\u003e2\u003c/sub\u003e creates more wrinkles and waves in the structure of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003e(Figure 1). The bond length of the Pd\u003csub\u003e2\u003c/sub\u003e cluster in the optimized structure of the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003eis 2.74 \u0026Aring; which shows the elongation of the Pd-Pd bond length compared to the isolated Pd\u003csub\u003e2 \u003c/sub\u003ecluster due to the strong interaction between the Pd atom, directed to the support, and the carbon atom of one of the tri-s-triazine units. The calculated partial atomic charges of these Pd and C atoms in the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e is 0.37e and -1.458e, respectively. Comparison of the charge of this C atom in the complex with that in the isolated g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e (0.297e) shows a strong charge transfer from the Pd\u003csub\u003e2 \u003c/sub\u003eto g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003econfirming a strong interaction between the Pd and carbon atom (see Figure 2). The angle between the Pd\u003csub\u003e2\u003c/sub\u003e axis and the g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003eplane is 46.0\u0026deg;.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eHOMO and LUMO of reactants\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTable 1 reports the energies of the HOMO and LUMO of C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e, H\u003csub\u003e2\u003c/sub\u003e, g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e, Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e, and C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003estructures along with their energy gaps. The shapes of the HOMO and LUMO of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e, Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e, and C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003estructures have been shown in Figure 3. According to our previous work [41], immobilization of the Pd\u003csub\u003e2 \u003c/sub\u003ecluster on the surface of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecauses the decrease of the energy gap of the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecomplex compared to the isolated g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e because of creating new electronic states between the HOMO and LUMO of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e [41]. These new electronic states cause the increase of the photocatalytic activity of the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecomplex compared to the bare g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e. Comparison of the difference between the energy of the HOMO of Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003eand that of LUMO of C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e (-4.25 eV) with the energy difference between the HOMO of C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e and LUMO of Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e (-5.24 eV) shows that the C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e is an electron acceptor species. A similar comparison between the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e and H\u003csub\u003e2\u003c/sub\u003e shows that the H\u003csub\u003e2\u003c/sub\u003e is also an electron acceptor. The smaller value of the energy difference between the HOMO of C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2 \u003c/sub\u003eand LUMO of Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e compared to the corresponding value for the H\u003csub\u003e2\u003c/sub\u003e and Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e indicates the better adsorption of the C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e on the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e compared to the H\u003csub\u003e2 \u003c/sub\u003emolecule.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAdsorption of H\u003csub\u003e2\u003c/sub\u003e and C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e on Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eEssentially for a heterogeneous catalyst containing the Pd cluster, it is important to use the maximum active site of the cluster after immobilization. Therefore, it is necessary to find the most stable site for the adsorption of H\u003csub\u003e2\u003c/sub\u003e and C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e molecules on the surface of the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst. For this purpose, the H\u003csub\u003e2\u003c/sub\u003e molecule approaches the catalyst in the dissociated and bonded forms, separately (see Figures S1(a) and (b)). As seen in Figure S1(a), the H\u003csub\u003e2 \u003c/sub\u003emolecule is located near one Pd atom which is farther from the surface of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e. In Figure S1(b), it is assumed that the H\u003csub\u003e2\u003c/sub\u003e molecule has been dissociated on the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst and each of its H atoms is interacting with one of the Pd atoms. The two initial structures lead to the same structure after the optimization so that the two hydrogen atoms of the H\u003csub\u003e2\u003c/sub\u003e molecule is interacting with one Pd atom (H\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e; see Figure S1(c)). The reason for the interaction of the H\u003csub\u003e2\u003c/sub\u003e molecule with this Pd atom could be attributed to its negative Mulliken charge which provides a good potential for the participation in electron donation to reactants.\u003c/p\u003e\n\u003cp\u003eIn the next step, the C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2 \u003c/sub\u003eis approached to the H\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecomplex to form the pre-reactant (C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.H\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e) for the start of the reaction. Figure S2 shows two proposed structures for the pre-reactant (R1 and R2). The best selected initial orientation of the C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2 \u003c/sub\u003erelative to the Pd\u003csub\u003e2\u003c/sub\u003e in the C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.H\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e complex is such that its molecular axis is in a vertical position relative to the axis of the Pd\u003csub\u003e2\u003c/sub\u003e (R1 in Figure S2). In this case, the simultaneous interaction of each C atom of C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2 \u003c/sub\u003ewith both Pd atoms is possible. It will be shown that the interaction of the C atoms of C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e with the Pd atom, located far away from the surface of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e, is stronger than the interaction with the other Pd atom. As previously shown in Figure 2b, the Pd atom located further away from the g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ehas a more negative Mulliken charge compared to the other Pd atom. Therefore, this Pd atom has good potential for electron donation to the reactants (H\u003csub\u003e2\u003c/sub\u003e and C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e). Therefore, the C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e and H\u003csub\u003e2 \u003c/sub\u003einteract mainly with the upper Pd atom in the optimized structure of C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.H\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e complex. Due to the adsorption, the H-H bond length increases to 0.93 \u0026Aring; compared to that of an isolated H\u003csub\u003e2\u003c/sub\u003e molecule (0.74 \u0026Aring;). Moreover, the increase of the C\u0026equiv;C bond length is seen for the acetylene in the C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.H\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e complex (1.28 \u0026Aring;) compared to isolated acetylene molecule (1.19 \u0026Aring;). These observations confirm the activation of the H\u003csub\u003e2\u003c/sub\u003e and C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2 \u003c/sub\u003emolecules by the catalyst and indicate that the initial complex (R1 in Figure S2) is suitable for starting the hydrogenation process. Considering the adsorption energies of reactants, it can be concluded that the adsorption of acetylene on the catalyst (-28.84 kcal/mol) is stronger than that of the H\u003csub\u003e2\u003c/sub\u003e molecule (-13.23 kcal/mol) which is in agreement with the frontier molecular orbital analysis performed in the previous section.\u003c/p\u003e\n\u003cp\u003eThe energy difference between the HOMO of C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e and the LUMO of H\u003csub\u003e2\u003c/sub\u003e (-5.93 eV) is smaller than the energy difference between the HOMO of H\u003csub\u003e2 \u003c/sub\u003eand the LUMO of C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e (-9.07 eV). This shows a charge transfer from the HOMO of C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2 \u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecomplex to the LUMO of H\u003csub\u003e2\u003c/sub\u003e molecule and leads to a strong interaction between these reactants and can better describe the co-adsorption configuration of the H\u003csub\u003e2\u003c/sub\u003e.C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecomplex.The calculated co-adsorption energy of the H\u003csub\u003e2\u003c/sub\u003e.C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecomplex is -45.3 kcal/mol which is greater than the value related to the adsorption of isolated C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2 \u003c/sub\u003e(-28.84 kcal/mol) and H\u003csub\u003e2 \u003c/sub\u003e(-13.23 kcal/mol) molecules on the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecatalyst. Eq. (1) is used for the calculation of the co-adsorption energy of H\u003csub\u003e2\u003c/sub\u003e.C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecomplex.\u003c/p\u003e\n\u003cp\u003e\u003cimg src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAA6AAAAA7CAYAAABljOesAAAXuklEQVR4Ae2cgZXjtg5Ft4XUkBbSQ0pIDWkhHaSDdJAKUkEaSAPpID3MPzf/378YLClL9kgjy4/neGVTJAhcgCApe/bLW0oIhEAIhEAIhEAIhEAIhEAIhEAIHEDgywFjZIgQCIEQCIEQCIEQCIEQCIEQCIEQeMsBNEEQAiEQAiEQAiEQAiEQAiEQAiFwCIEcQA/BnEFCIARCIARCIARCIARCIARCIARyAE0MhEAIhEAIhEAIhEAIhEAIhEAIHEIgB9BDMGeQEAiBEAiBEAiBEAiBEAiBEAiBHEATAyEQAiEQAiEQAiEQAiEQAiEQAocQyAH0EMwZJARCIARCIARCIARCIARCIARCIAfQxEAIhEAIhEAIhEAIhEAIhEAIhMAhBHIAPQRzBgmBEAiBEAiBEAiBEAiBEAiBEMgBNDEQAiEQAiEQAiEQAiEQAiEQAiFwCIEcQA/BnEFCIARCIARCIARCIARCIARCIARyAE0MhEAIhEAIhEAIhEAIhEAIhEAIHEIgB9BDMGeQEAiBEAiBEAiBEAiBEAiBEAiBHEATAyEQAiEQAiEQAiEQAiEQAiEQAocQyAH0EMwZJARCIARCIARCIARCIARCIARCIAfQxEAIhEAIhEAIhEAIhEAIhEAIhMAhBHIAPQRzBgmBEAiBEAiBEAiBEAiBEAiBEMgBNDEQAiEQAiEQAiEQAiEQAiEQAiFwCIEcQA/BnEFCIARCIARCIARCIARCIARC4HMI/P3332+//PLL219//fWQAr///vvbb7/99pCM3Q6gX758edvy+u677x4yJJ1D4AoE/vzzz03zhjn2008/XcH02LCRwJb8Stvk2I2AL9D8xx9/3JxPHt2YXADb26+//rqZ26ObsStwiw0hEALnJfDHH3+8sSZwCF0q7ENpRx5cKhxCaffPP/8sNZve2+0AioHff//9/5P4SAOU/uGHH/5tk030iFDqXpFAPYQyuUeFTSIHCg4W2fiMCF2/Ljn2+j5+1ELWWNZWH1aQW0bl559//rdNHlJ8pVMPobONGBsw2ebg/pVd3oVACJyLALmK89bSYdE25rRZ3quW0Yez3q1Dbe3j+90OoAzA5hlDMHpWTPLZRM8Ipf4VCZgA+KnErDi/svGZEbp+vTGQHHt9X99roWssOWW2+fChVx4Ef6UsE7jxzcGscD8H9xmd1IdACHw2AfIXeWppr0i+42Vb2q85gGIb68atw+2Iwa4HUL8B5enqrJjkl8DM+qY+BK5IgCdJTH5ePF2aFZJDNj4zOq9Rnxz7Gn5+xEq/3SRWlgr5Jg+CvxIi95qHl57u8xAoB/ev3PIuBELgPAR46Mg+cekc1rXlMEnuW3sAdYyteXC3AygKmbyzqHX35nMIzAnUJ1B5MDPn9Op3kmNfPQLW2e9mYuvmYJ3067bi1yfsYfKQ77o+jmUhcHUCPoBc+hVHZ+Avq9YeQOnvOHypuLbsdgCtm+gtCq1VPO1C4KoE3Piw+UkJgRmB5NgZmdRLoD6k2LKZsP8rXz24sxlLCYEQCIFnI+Cv6bY+RLvnAOp+ZEu+3G2Hu7SJ9qctz+bM6BsCRxBY2vjwlGnLBD9C34zxOQSSYz+H+zON6qaAh1n9QbCbk17/TPbtpeutgzs/Z86Bfi/6kRsCIfARBNwjbP31yz0H0Joz1/5yb7cDqAawme6Fe9lEdyr5fAYCxq0/H7/n+mhsOybJoxYmOE+ysvGpVF73vbGaHPu6MXDLcnKF+YT8UYv3al3e/5cAh3K5cYivxXs5uFcqeR8CIXA2AuwXyWNb94zuLbb28/+k6HvXGZfdDqAm7/6Hry56axWcKZ76ENiDgBPP+L3n+sgB1M0N49b/gIjNI0+xqO8boj04ROb5CRibybHn99VnaWg+6/8BkXlm9PDis3Q907juU5hj9T8g4r2brH6gP5P+0SUEQuC1CfAtpHuEupdcQ8V1Y+sB1H5r15VdDqAubho/umYTvSYM0ubVCNSNz2jeUJeNz6tFxbf2Jsd+yyQ13xKY5RDr8yD4W2bUuJGSU7+u3WCNpac2BEIgBPYl4J86krvYL2wp5r+tB1C/JGHMNWVdqzWSShv+11sTdn16SBMV7PWle96GwMsScH70byw4dFLX618W1AkMX3MINA8uXbcuDpieHHuCAFipgov5UgzcuoeMrWXpCbjxs/XJ+FYdHmm/5mHcLW7cv6f407X+6wKZ9vp7xkifEAiBENiLQM2fW/cYrllbD6B1zDV/B3pfdr5BbLaJphsKbv0fmW4Ml9tPRIANQfw/d9hs40MPkgJzK+UcBD7zAJoce44YWKOFi/maw9KszT0HUA+ZyOwPfI3dNZuENTbu0aZuZmZc1tRv1c1DJrJHB3TqYZsSAiEQAmclUPPnZxxA14y5ywF0aRN9VmcdpZcLvwvnK/2cUttziBpH262Nz7jX9WuNm1ecMzPvvlqOrYtpfv44i4r39UsPKd63HH+qzJl79Vu/Piev5JOlg/uY1GvV1gcq+Qn3a/k+1j4PgZq/yddbinMcGVvK1jE//ADKk1Y3inlKOHadjK60aHOQZlP87AuSE88Yvud6z7cVREr9zf6Zv5kYR/W+tR7OrzRn7iVm/iA2z5pjmQMf7SsPBkfnGA9bz/b/Fvif5TzywA/WxNlIhv5gnCs9SOWgjc2sZ89e9opdH24825x4dn9G/xBYS2DrYbDKdR/8yAF0zR72ww+gazfRwqlGv8p7kjYLXH2i/Oy2689nX5CcePjn3te9B9C1Gx/k3zvGs8aZG6mjDx9n5HX2HLuXrz4rx5gTnumQtfYhhb7iOirmpFFeNw6vNifXHNx94Lp1gzZivGfdXrG7l9w9WUR2CLwSAfMz+9hZfp/xcH5vzW8+sGTMNWVdqzWS/tfGBWtJAZM3Rs4KCyiy/KkZ8kbtezsWj63QZjr0evRGtguUY6HX7JuIkX46lwDphfY+XcRm7R5tAEayu+3oLEdY+lQCedrBfUodm28v6EtBpn7omw1soK261qv2ucnxnnL/FV7+oV21HbndbvS3jToiD73UsTMoQ5z6rf4YxbmKy3LJxs4RLt1vyHskfoj3ru8tP9T7s1jE9lF8+G2LMSUPr/SBSY3FbjNtjB10r+MgF514nb18do6Fsf43n8CM9+Y257pXc0xl2+MU3+BDfDMq+o64tXSf1nu0cXyu1d/c0wbbMC5jWIgfY9w2Xvs49jnLlXhWVzjPiv6a3Xc+dXa0hw9j9Bw9k/UM9fhVbks5lnu0m7GFl7FMO+KIz/Am7kal96Ef/Oscsx/+1Xe043310ZrYRa7+pX8t2IVcdaj3eE+enM1T7t/Sj5hxbHMD7J3jyF+KK3MH7dARXUacut75HAKvQoD54Bye5ZwZC3PLUg4c9bUfc3tN+fADqAt2T2hVGReumXEuniQjkhKJFbndKACTgKh3Q4BMoM9kVz22vEcHxmE8F53qYHTuZaQfgWBQqLP9GKPbY/seQCPZI9vRGd1cUGgDf5I746EL92WMXPnTBz9qr4tDXxiU0/2jXVxlNWvjmCxGyOPleI4PLxcabUUXddQGbOL9MxVsMy6WYhd+tJNJt9HDCTJgIHd8Xsuj8WM86qOZH+qY+Ak/a6uxqCxjveuKDMfpc4Z72EJc8JKL7WusUsf4tIGhnI096pbyVrXlM99j5y1d98yx+gkdRv4wRpfmoMyd79XH2Dcq5MZ+D3/hf3VyM2t/47/3w/fob+xhB3ozRi3ED+1GMVnbne29eQDdZ0XbZjGP/+jf2SlvjZ9t+yxX4xK7zSVdd2KFOJmxhZuxpAzjkD6M0Qv3XfvtYw6r84j3+Iu2xq6x32NU//Z6xsYG5CDPcbpeyp3NqV6P3DX6oRd8bAsT7Ecf9JIVbUbF3MYVGbfaj2SkLgRegQC5m/nlXmetzczFe/q5Joxyzmjs+eo0an2jzo3dTHGShQs/bUy0VazJBAC18LkmPBcBACO3FmTPFs3abst7HdJ1th59akEnFgleVT8ZUd+L9/pCgFPruFtsd2xlI4sFYFRsq4+wDX9YrO/BrOzqH/t4XWrjQslCWMuozy0dZ/6ocs/4XrazeYHP3SjQZlTcMHT/0L7G1KPxU2WpN3FVY0U/VD31nX6lT5VlPX17cS71emQy17lf56Ax5SaNfo7Pe/oY3/SlHX2WYriP/RmfZYRPu5+1UZ/MYmltjoWX3Cs7xkH2LMcu3aOvNnQ/Wz/ygTr3HKFeXEfj2q/GGTowdtcf/3edZNn7I+OsRb/Bo9ujztiDb2kziiPaOYdG/pA3c/hKBV4w4TUqxCg2L7FVBm1rsb7mKe7D0hzGewt+6THqGoBvLPqiz41bsetYMz+7nnR9ra+5VV3W6ufYMiHGrEMWfLvt1Dt23+DiE+I5JQRC4CsBH9Ywz9YW12HnYJ//S3Low6vuBRfbL93cco/E4IKmEreuI/kmpFtGk7CQP9oYOO5I/j11OmS0EDPWKFHOkr+y+mKBXm6WSKZL9t9ju4l7NG5novy+wGgT11qsH/nCdtyDVZfJffhxr9ssq9HkcaHruny079V/z6u81f3WdcSDxZv5N4rFrrvjjfzl2L3PLH70Q/crcmYbAmOhx+LM386L3h4djb0eB+o1S4Teh9eIQ7f/DJ9hfGSOlW3no5+Io16W7tl2luNn49HP+OtxpkyuHgzqRpb2ozmh/5fkIXOmax33TO/hbz51Lt+60mdU3Lws9e8HgZGcZ6gjZoyfJXvrvZ5zsHMp/uk7isWluK/slN3z/2xurI1dbEc37K+FOTKa486dnlu36sdY5rMagzM5rnH0qXN8aX2o9uR9CLwaAfbUzO3ZXqzycN7VHFff17aj9/bv+WnU1rrxYz7vHnw1kfREOFLDxFUTke2E5meuwvHerWvta8KtTx25r3O538tMv9liYX8PB/QfLXC0m8nmnnYpz6ty+6Lh/Xp1Ie5sXSi7XmsWutn4Sz7XZ6OAdpNVdVTWqH2174rvjas1G8KPjB/8gLwtfnBj2+fTzN/a1uMOPxoHzEVeyDQe6TcrxvJSm1nfZ653jjySY32AMGIn135oldnS+LMcQV9zMDEyK8ZVbUN8jHRBD+cB8cLnUSGfIuMVi/NoxEY/9zlMHoA3vqQ/L+Kk5ocrszROO5d79gqdk/MDvvAkzq0b7UG2xK5rvmOqL9denDe9Xl3W6ucYPRcRL+je8z1yqWccCvbTBn2QMdK165jPIfBqBJyXPSd9NIfR+ntrjFMdQGeJpxvhJobFrReSEkmq37Oee2teVa4Jty+iJsS+EVO/nliR6QKFPrNCoLix7jKU3e1DljaO7ilvNqb12Djix33kcq9vSODDa6nMxl/yuZucvhC5cHU7lbXmELak6zPeWxNX2PWR8TOLFf3Q/SZXYpo46kV/9/lkAu1zxjio85mYwP89RvtYyuxj9XZX+3zLN9prnPT8w33ZjRgbh6N79HX80RwlR8zyyCx/qG+VrU/JzyP97UP8agsx1BdoGYw298q48hUmM3+4FvRNv/XMZearc7rn6qtygxfctu4V1vBxDpjviG3it8ctbLfGbp+3yB3NUXPuSN8t+qGj+6e+TqhLt6vOVRgwHm2R03lfNb5iVwhsJcDcIC+N5uxWWbP2jrF1rfx2Fzgb4YB6F6uekPrQLGwkoBFQk5SbkN73ns8m/N7XhNs3xjP9TN7IW1NMxFX+TDbyZrYTHDNeXQ/l98VH3bG5Fhe6kS9sZ99RmyWfy7dvZlmYsKfrKK++cKnHla+whUmNlZG9+nfki63xo6w+X2/pgp6jg4H+Jl5qsb7W8d7xR7b0tvUzdqoj71+pLM23ymHG1rmMD0fFDfjoHnWz8d2MjnzpmKN7dRx1ZgwXROpulVkO87Dc4/uWvCvcl+VsfuD/vhZgNz7qfZjrs3i5AqtqA3aObDWH9XiU863YZoyZ7Dq+77fGrvMSfXgxj5lDvThPR3Nii37IZf2mT1+vzSF9fBjRvjPsOuZzCITAewLmmb3mDvmAHNfn7Hstvv003kV82+6QmlnSxKiapGYbEjcS94BYMnCUWAXOvV50dl9UWJhp3+vpz71+aK6LgmPcY/vswKbMenVMrrXMDncjW/FDDUTbdJnIny1oxsLoiYoLF31rcZGvY9f7V36vf3qCwRe8LB8ZP8ZKnZuM4wbCMas/jIU+B/R3f6iADObM6MDqfO+y6ENsVLtpw9jMMzdXyCRmLFVP6652lXOfi9he/Tjz01IOg1XPlcitfjBm6vi08ZBS62WvLqPYsI1X8yvzYdSe+Mf/vdiv1o90RZdXKCPbtdv1ZJSb8TWMa2GO1XlW713tfY9/7IOla1O391YOq+17XvUeMvpaOPLfUuzWOYauNRc4Dlflju5v0Q9ZzvkaL65Po3w/W+PIH+iVEgIhMCdAjnDvM2+1/Q75B7lct5ZvT09bJXxg+5p8SCoUEh0JsW8m2FCS7E2qXIFwL4glM0x86GKyw5nqYF1NguqCTbzYuLnYkFzpg1wXDuqwk7YUrtQhRxbq6LhrbXfRcCzljK7KRjfG5QV7WFf77Fttog6d+gJGHf2Rg12Moe7Ix8Zqu5tkOTmWV3WUlfXIYRwKuvaYsd0Vr/gW24kzClzhCNfud/npA66w4zVKIrP4Uc6SH9Cnxo2yur8ZG3m9IBu7jAXk1c2Pc0S9aU/sUm8dMpHhy3rjGvthRL+rl84Te2c51pigD69ZDjOOkEUfOBN/9MEP1f+w5z6+pg199Sn1+IQ23kemsV3zB7JHhbhCB67I7wVbGafqZL6pdfSz3jzifJqN3cd61s9wgx+csLkX5gn31swXfVfnbJd3pc+ywV44ElMwMFdaV2PNmDUvGf91DsDIfIU8irKI9853a+wS0/gUWUt+Za7SjvGYu+ho2aIffZCD7bX0/IAu2saVPuhggQU6V57eyzUEQuA9Aeesueb93e2flHfvmniqAyjmk2Rc/EyGGNkLyZdERBsTGUmI+o8uOEud2FCpj4sr9/pCbWJEN5KoDjeB08fEir7UI5v2vLCLupFjt9ru4qcOS3wYF91cSNEFvaquvb8Lj7Z2ndFX29Cl68HnPp6LbB+Lz/Lp96o/XnFBwuY6H/aOn3v8oJ/Rreq65C/aMlafM/h/TewQf6P+xCkyuccYe+SOHqNn+Lw2x/Yc5rzWHzUXalftM7pPO/MF/uc93GuOIEYci/b13ih/ODZXc505ut7jPfW0MfbwPXqO8g3jKq/q2mVe7bM2w4ZXzf36zntL89Z8PGJ7NWbaQz4yp9T4lwX3+l6BPpU570fMiMdR3qS+l3tiF58S5yN5ykdX507PmVv0Yx4yHnO9l5pfauzRjniTL3rQtuaKLiufQyAE3hNgnpLHH503zM2l/P9+1PGn0x1Ax2qm9ggCLC4sCiT1lBDYgwCbBjYQKSHw0QTcHI827x89VuQtE8AHzPXZg4Dl3rl7NAH9xRxKCYEQCIEjCOQAegTlJxmDRYgDaH9C+yTqR82TE+CJ2+yp98lVj3onJEC+4tsiCjnLb0ROqOpLqeQ6ksPned3OQ2a+vfDbEPIyfksJgRAIgaMI5AB6FOknGMefvmQhegJnPaGKbkzzDfsTOu+EKrOBZuPsK3H1+U7yp5X9ISa+ybdrn+8fNag/+c3hUyq5hkAIHEkgB9AjaZ98LP5mxc3co7/tPrmpUe8TCNS/H/Obq09QI0NehIB/T0ze6n8rdhETn84M/z6P+e3LdeXRvzl6OhgnVphfC7DWM4fyYODEjopqIXBhAjmAXti5MS0EQiAEQiAEjiLAN6Cz11E6ZJwQCIEQCIHzE8gB9Pw+ioYhEAIhEAIhEAIhEAIhEAIhcAkCOYBewo0xIgRCIARCIARCIARCIARCIATOTyAH0PP7KBqGQAiEQAiEQAiEQAiEQAiEwCUI5AB6CTfGiBAIgRAIgRAIgRAIgRAIgRA4P4EcQM/vo2gYAiEQAiEQAiEQAiEQAiEQApcgkAPoJdwYI0IgBEIgBEIgBEIgBEIgBELg/ARyAD2/j6JhCIRACIRACIRACIRACIRACFyCQA6gl3BjjAiBEAiBEAiBEAiBEAiBEAiB8xPIAfT8PoqGIRACIRACIRACIRACIRACIXAJAv8BArKo8erqzCwAAAAASUVORK5CYII=\" alt=\"\" /\u003e\u003c/p\u003e\n\u003cp\u003ewhere \u003cem\u003eE\u003c/em\u003e\u003csub\u003ecomplex\u003c/sub\u003e is the electronic energy of H\u003csub\u003e2\u003c/sub\u003e.C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e.Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecomplex. The \u003cem\u003eE\u003c/em\u003e\u003csub\u003eacetylene\u003c/sub\u003e, \u003cem\u003eE\u003c/em\u003e\u003csub\u003ecatalyst\u003c/sub\u003e, and \u003cem\u003eE\u003c/em\u003e\u003csub\u003eH2 \u003c/sub\u003eare the electronic energies of isolated optimized acetylene, Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4,\u003c/sub\u003e and H\u003csub\u003e2\u003c/sub\u003e molecule, respectively\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMechanism of acetylene hydrogenation on the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003ecatalyst\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAs previously mentioned, two initial structures, R1 and R2 (Figure S2), were considered as the starting structures for the hydrogenation mechanism. Zhao et al. demonstrated that the hydrogenation pathway via the R1 complex in the presence of Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003esingle-atom catalyst requires less activation energy compared to the R2 ones [39]. Also, they reported that the hydrogenation pathway via the R1 complex had more selectivity to the ethylene production rather than ethane production [39]. Therefore, the R1 structure is selected as the initial structure for the hydrogenation mechanism on the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst.\u003c/p\u003e\n\u003cp\u003eConsidering the R1 as a starting complex, two consecutive pathways are proposed for the hydrogenation mechanism. In the first pathway, the acetylene is converted to ethylene (IM3) via the three transition states (TS1, TS2, and TS3) (see Figure 4). In the second pathway, ethylene (IM3) is hydrogenated more and converted to ethane (P) via two transition states (TS4 and TS5). To obtain the activation energy of each step, the \u003cem\u003eE\u003c/em\u003e\u003csub\u003eco-adsorption \u003c/sub\u003eof each complex involve in each step is calculated using Eq. (1). The energy difference between the \u003cem\u003eE\u003c/em\u003e\u003csub\u003eco-adsorption \u003c/sub\u003eof the reactant and TS of each step is considered as the activation energy (Figure 4). Figures 5-9 show the optimized structures of the complexes (R1, IMs, and TSs) involved in the pathways shown in Figure 4. Also, the important interaction paths for each complex, obtained from the QTAIM calculations, are shown in the Figures. The position of the bond critical points (BCPs) (dots) and their calculated electron density (\u003cem\u003er\u003c/em\u003e) have been identified in the Figures.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u003cem\u003eR1\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e\u0026reg;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e TS1\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e\u0026reg;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e IM1.\u003c/em\u003e\u003c/strong\u003e The calculated value of the activation energy for the conversion of the R1 complex to IM1 via TS1 is 6.77 kcal/mol (see Figure 5). The optimized structure of TS1 is similar to the R1 complex except that the H\u003csub\u003e2\u003c/sub\u003e bond length increases to 2.77 \u0026Aring; in the structure of TS1 indicating the complete dissociation of the H\u003csub\u003e2\u003c/sub\u003e molecule on the catalyst. The calculated vibrational mode of TS1 with the imaginary vibrational frequency (n=-470\u003cem\u003ei\u003c/em\u003e cm\u003csup\u003e-1\u003c/sup\u003e) is shown in Figure S3. The direction of the displacement vector of one of the H atoms of the dissociated H\u003csub\u003e2 \u003c/sub\u003emolecule is towards the C atom of the acetylene indicating the H transfer from the Pd to C atom. Comparison of the arrangement of the BCPs in the R1 complex with that of TS1 shows that the BCP No.1 (BCP1; \u003cem\u003er\u003c/em\u003e=0.1732 a.u), related to the interaction path between two H atoms of H\u003csub\u003e2\u003c/sub\u003e molecule in the R1 complex, disappears. Also, a new BCP (BCP4) is formed between one of the H atoms of the H\u003csub\u003e2\u003c/sub\u003e molecule and the C atom of acetylene in the TS1 complex. As seen, this H atom has simultaneous interaction with the Pd and C atoms via the interaction paths related to the BCP2 (\u003cem\u003er\u003c/em\u003e=0.1189 a.u) and BCP4 ((\u003cem\u003er\u003c/em\u003e=0.1444 a.u) in the structure of TS1. Due to the interaction of this H atom with the C atom of acetylene, the C-H bond length increases from 1.08 \u0026Aring; in the R1 complex to 1.43 \u0026Aring; in the IM1 complex. These observations confirm the transfer of the H atom from the Pd site to the C atom of acetylene for the formation of the vinyl compound. It is important to note that there is no BCP between two Pd atoms in the structure of the complexes shown in Figure 5 and each Pd atom is interacting with the C atoms of acetylene independent of the other one. In the R1 complex, each of the Pd atoms is interacting with one C atom of acetylene while one of the Pd atoms is in the direct interaction with two C atoms, simultaneously in the TS1 and IM1 complexes via the two BCPs (No. 5 and 6 in TS1 and No. 3 and 4 in IM1). The increase of the value of \u003cem\u003er\u003c/em\u003e for the BCP of Pd-H bond (BCP2, BCP1, and BCP1 in the R1, TS1, and IM1 complexes, respectively) shows that this bond becomes stronger due to the H transfer. Also, the strength of the Pd-N bond between the Pd atom and g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e becomes weaker in the IM1 (BCP No. 5; \u003cem\u003er\u003c/em\u003e=0.7159 a.u.) complex compared to the R1 complex (BCP No. 5; \u003cem\u003er\u003c/em\u003e=0.7015 a.u.) due to the decrease in its \u003cem\u003er\u003c/em\u003e value.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u003cem\u003eIM1\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e\u0026reg;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e TS2\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e\u0026reg;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e IM2.\u003c/em\u003e\u003c/strong\u003e Before the transfer of the second H atom in the structure of IM1 to the C atom of vinyl for forming IM3 (see Figure 7), a step is necessary (IM1\u0026reg;TS2\u0026reg;IM2). The second H atom is interacting with the Pd atom via the BCP1 in the structure of IM1. Notably, this H transfer becomes easier when the orientation of the complex above the g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e changes so that this Pd atom can directly interact with the N atoms of the substrate (see Figure 6; IM2). Due to this interaction, the IM2 complex is more stable than the IM1 complex and the orientation of the Pd-H bond is more suitable for the second H transfer to the C atom of vinyl compared to what is observed in the IM1. The bond angle of C-Pd-H in the structure of IM1 is 160.8\u0026deg; and reaches 83.9\u0026deg; in the structure of IM2 so that the distance between the second H atom and the C atom of vinyl decreases (see Figure 6). As seen, the orientation of the complex on the g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e changes so that the two Pd atoms are in direct interaction with the N atoms of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e. The strength of the H-Pd bond does not change due to this configuration change because of the small change of its \u003cem\u003er\u003c/em\u003e value. The calculated value of \u003cem\u003er\u003c/em\u003e of the H-Pd bond in the IM1 (BCP1) is 0.1484 a.u. and the corresponding value in the IM2 is 0.1502 a.u. (BCP5). Moreover, a decrease is seen in the bond length of the formed C-H bond in step 1 in the IM2 complex (2.36 \u0026Aring;) compared to that in the IM1 complex (2.60 \u0026Aring;) which shows that this bond is stronger. Therefore, the transfer of the second H atom from the Pd atom to the C atom in the IM2 needs less activation energy compared to that of this process in the IM1. A considerable decrease in the \u003cem\u003er\u003c/em\u003e values of the C-Pd interaction paths of the TS2 structure compared to those of the C-Pd interaction paths of the IM1 structure is seen. For example, the \u003cem\u003er\u003c/em\u003e of the BCPs No.3 and 4 (0.9617 and 0.9885 a.u, respectively), related to the simultaneous interaction of one Pd atom with two C atoms in the structure of IM1, reduces to 0.3508 and 0.4902 a.u., respectively in the structure of TS2 showing the decrease of the strength of these bonds (see Figure 6). Also, the arrangement of the interaction paths in the IM2 shows that each Pd atom is mainly interacting with one C atom unlike what is observed for the IM1 and TS2 structures in which one of the Pd atoms is in the direct interaction with two C atoms, simultaneously. One important point related to the TS2 structure is that the orientation of the complex on the g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e is so that its interaction with the substrate takes place mainly via the H atoms (see BCPs No.4, 7, and 8 in Figure 6). The imaginary vibrational mode of the TS2 is shown in Figure S3. As seen, the H atom which is interacting with the g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e via the BCPs No. 7 and 8, moves out of the substrate plane.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u003cem\u003eIM2\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e\u0026reg;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e TS3\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e\u0026reg;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e IM3.\u003c/em\u003e\u003c/strong\u003e Figure 7 shows the transfer of the second H atom from the Pd atom to the C atom of vinyl via the TS3 structure and form of the IM3. The comparison of the structure of TS3 with IM2 shows that the bond angle of C-Pd-H related to BCPs No. 2 and 5 in IM2 decreases in the TS3 structure (see BCPs No. 2 and 4). Comparison of the value of \u003cem\u003er\u003c/em\u003e of the BCP5 in the IM2 (0.1502 a.u.) with that of the corresponding BCP in the structure of TS3 (BCP4 ; \u003cem\u003er\u003c/em\u003e= 0.1341 a.u. ) shows that the decrease of the strength of the Pd-H bond due to the decrease of the C-Pd-H bond angle. The imaginary vibrational mode of the TS3 is shown in Figure S4. The displacement vector of the H atom shows its transfer from the Pd atom to the C atom. Notably, a BCP (\u003cem\u003er\u003c/em\u003e=0.5594 a.u.) is seen between two Pd atoms in the structure of IM3 showing the interaction of two Pd atoms with each other which is absent in the structure of IM2 and TS3.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u003cem\u003eIM4\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e\u0026reg;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e TS4\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e\u0026reg;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e IM5.\u003c/em\u003e\u003c/strong\u003e Figure 8 shows the first step of the hydrogenation of ethylene to form ethane via TS4. A new H\u003csub\u003e2\u003c/sub\u003e molecule is added to the IM3 for forming the IM4 so that the H\u003csub\u003e2\u003c/sub\u003e molecule is located near the Pd atoms to feasible the breaking of the H\u003csub\u003e2\u003c/sub\u003e bond. As seen in the figure, the dissociation of the H\u003csub\u003e2\u003c/sub\u003e molecule takes place in the TS4 structure so that one of the H atoms of H\u003csub\u003e2 \u003c/sub\u003eis connected to one of the Pd atoms via the BCP2, and the other one is pulled toward the carbon atom of ethylene so that it is in the direct interaction with three atoms via the three BCPs (BCP4, BCP5, and BCP3), simultaneously. The imaginary vibrational frequency of the TS4 is shown in Figure S4. As seen, the displacement vector of one of the H atoms shows its tendency for forming the bond with the C atom of ethylene. The imaginary frequency of this vibrational mode and the barrier energy for this hydrogen transfer is 409\u003cem\u003ei\u003c/em\u003e cm\u003csup\u003e-1\u003c/sup\u003e and 12.28 kcal/mol, respectively. The comparison of the structure of IM4 with TS4 shows the decrease of the Pd-Pd distance so that a new BCP appears between two Pd atoms in the IM4 structure. The strength of the Pd-Pd interaction in IM5 is lower than that in IM4 because of the decrease of the \u003cem\u003er\u003c/em\u003e of BCP3 from IM4 to IM5. Comparison of the \u003cem\u003er\u003c/em\u003e value of C-H bond in TS4 (BCP4 ; 0.1557 a.u.) with the corresponding value in IM5 (BCP1 ; 0.2384 a.u.) shows the complete formation of the new C-H bond in IM5 compared to TS4.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u003cem\u003eIM5\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e\u0026reg;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e TS5\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e\u0026reg;\u003c/em\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cem\u003e P.\u003c/em\u003e\u003c/strong\u003e Figure 9 shows the last step of the hydrogenation mechanism related to the second H transfer from the Pd atom to the C atom of CH\u003csub\u003e3\u003c/sub\u003eCH\u003csub\u003e2 \u003c/sub\u003evia the TS5 structure to form ethane. As seen, this H atom interacts with the Pd atom via the BCP5 ( \u003cem\u003er\u003c/em\u003e= 0.1960 a.u.) in the structure of IM5. Also, the increase of the number of BCPs between the Pd atom and g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e is seen in the TS5 structure (BCPs No. 3 and 4) compared to the IM5 structure (BCP No. 6) which helps to the hydrogen transfer by the weakening of the Pd-H bond as well as the presence of the BCP1 between the H and the C atom in this structure. Comparison of the calculated value of \u003cem\u003er\u003c/em\u003e of the BCP2 (0.1484 a.u.) in the structure of TS5 with that of BCP5 (0.1960 a.u.) in the structure of IM5 shows the decrease of the strength of the Pd-H interaction and confirms the H transfer from the Pd atom to C atom. The calculated value of vibrational frequency of the imaginary vibrational mode of TS5 (see Figure S4) related to the H transfer is 535\u003cem\u003ei \u003c/em\u003ecm\u003csup\u003e-1 \u003c/sup\u003ewhich is in agreement with the literature [39]. After the transfer of the last H atom to the CH\u003csub\u003e3\u003c/sub\u003eCH\u003csub\u003e2\u003c/sub\u003e and the formation of ethane (structure P in Figure 9), one of the Pd atoms becomes close to the substrate so which creates a folding in the structure of g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e.\u003c/p\u003e\n\u003cp\u003eA comparison has been made between the imaginary vibrational frequencies of the TSs of the hydrogenation of acetylene using the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst with that reported in the literature [39] for Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e (see Table 2). As seen, the calculated vibrational frequencies for the R1\u0026reg;TS1\u0026reg;IM1, IM4\u0026reg;TS4\u0026reg;IM5, and IM5\u0026reg;TS5\u0026reg;P steps of our work are smaller than the value reported by Zhao et al. [39]. The value of the imaginary vibrational frequency of the TS of IM1\u0026reg;TS2\u0026reg;IM2 is nearly equal to that reported by Zhao et al. for the Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst. The imaginary frequency of the transition state of IM2\u0026reg;TS3\u0026reg;IM3 is greater than that reported for the Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSelectivity and reactivity\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAs previously mentioned, the effective performance of a heterogeneous catalyst in the hydrogenation of acetylene is the production of ethylene with a smaller value of activation energy compared to that of ethane production [39-40]. According to Figure 4, the activation energy for the first H transfer reaction (R1-TS1-IM1) is equal to 6.77 kcal/mol and it is smaller than the values reported in the literature for the other heterogeneous catalyst (see Table 3). Also, the comparison of the activation energy of the R1-TS1-IM1 step obtained in this work using the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst (6.77 kcal/mol) with the corresponding value reported by Zhao et al. [39] using the Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst (21.53 kcal/mol) shows that the increase of the number of Pd atom has a considerable effect on the decrease of activation energy of this step.\u003c/p\u003e\n\u003cp\u003eThe selectivity of a heterogeneous catalyst in the acetylene hydrogenation, for the production of ethylene, increases by the decrease of the activation energy of the ethylene production step compared to that of the ethane production step. The selectivity (D\u003cem\u003eE\u003c/em\u003e) is considered as the difference between the hydrogenation barrier (\u003cem\u003eE\u003c/em\u003e\u003csub\u003ehydrogenation\u003c/sub\u003e) and desorption barrier (\u003cem\u003eE\u003c/em\u003e\u003csub\u003edesorption\u003c/sub\u003e) of ethylene [39-40]. The higher value of ∆\u003cem\u003eE \u003c/em\u003eleads to the higher selectivity of the catalyst for ethylene production. The calculated value of ∆\u003cem\u003eE\u003c/em\u003e obtained in this work is 21.9 kcal/mol (see Table3) which is more than that reported by Zhao et al. for the Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst. This shows that the increase in the number of Pd atoms of the catalyst is accompanied by the increase of the selectivity of the catalyst for ethylene production.\u003c/p\u003e"},{"header":"Conclusion","content":"\u003cp\u003eIn this work, a theoretical mechanistic study of the hydrogenation of acetylene on the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e heterogeneous catalyst was performed using the DFT and QTAIM calculations. The selectivity of this catalyst for the ethylene production was compared with those of the other catalysts reported in the literature (Pd(111) [45], Al\u003csub\u003e13\u003c/sub\u003eCoO\u003csub\u003e4\u003c/sub\u003e [46], AlPd [47], Pd1/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e [39]). It was observed that the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e has more selectivity and activity for ethylene production compared to the mentioned catalyst. Also, the results of this work show that the clustering of the Pd atoms on the g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4 \u003c/sub\u003eincreases the selectivity and activity of the catalyst for the acetylene hydrogenation. This work provides strong insights for further research on the role of different palladium clusters deposited on graphitic carbon nitride material in acetylene hydrogenation.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAcknowledgment\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors thank Isfahan University of Technology (IUT) for its supports. The authors gratefully acknowledge the Shaikh Bahaie National High-Performance Computing Center (SBNHPCC) for providing a computing Facility.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors thank Isfahan University of Technology (IUT) for its supports.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eThe contribution of Authors\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll authors conceived and designed the calculations; analyzed and interpreted the data; contributed materials, analysis tools, or data and software; and wrote the paper.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConflict of interest\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;The authors declare that they have no conflict of interest\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and material\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll data generated or analyzed during this study are included in this published article [and its supplementary information files].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCode availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eN/A\u003c/p\u003e"},{"header":"References","content":"\u003cp\u003e[1] Pradier C M, Mazina M, Berthier Y, Oudar J (1994) Hydrogenation of acetylene on palladium. J Mol Catal 89: 211-220.\u003c/p\u003e\n\u003cp\u003e[2] Sarkany A, Geszti O, Safran G (2008) Preparation of Pd shell\u0026ndash;Au core/SiO\u003csub\u003e2\u003c/sub\u003e catalyst and catalytic activity for acetylene hydrogenation. Appl Catal A: Gen 350: 157\u0026ndash;163.\u003c/p\u003e\n\u003cp\u003e[3] Wang P, Yang B Influence of surface strain on activity and selectivity of Pd‐based catalysts for the hydrogenation of acetylene: A DFT study (2018) Chinese J Catal 39: 1493\u0026ndash;1499.\u003c/p\u003e\n\u003cp\u003e[4] Lopez N, Vargas-Fuentes C (2012) Promoters in the hydrogenation of alkynes in mixtures: insights from density functional theory. Chem Commun 48: 1379\u0026ndash;1391.\u003c/p\u003e\n\u003cp\u003e[5] Studt F, Abild‐Pedersen F, Bligaard T, S\u0026oslash;rensen R Z, Christensen C H, Norskov J K (2008) On the role of surface modifications of palladium catalysts in the selective hydrogenation of acetylene. Angew Chem Int Ed 47: 9299\u0026ndash;9302.\u003c/p\u003e\n\u003cp\u003e[6] Vile G, Albani D, Nachtegaal M, Chen Z, Dontsova D, Antonietti M (2015) A stable single-site palladium catalyst for hydrogenations. Angew Chem 54: 11265.\u003c/p\u003e\n\u003cp\u003e[7] Huang X, Xia Y, Cao Y, Zheng X, Pan H, Zhu J (2017) Enhancing both selectivity and coking-resistance of a single-atom Pd\u003csub\u003e1\u003c/sub\u003e/C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst for acetylene hydrogenation. Nano Res 10: 1302\u0026ndash;1312.\u003c/p\u003e\n\u003cp\u003e[8] Kyriakou G, Boucher M B, Jewell A D, Lewis E A, Lawton T J, Baber A E (2012) Isolated metal atom geometries as a strategy for selective heterogeneous hydrogenations. Science 2012 335: 1209.\u003c/p\u003e\n\u003cp\u003e[9] Liu J, Bunes B R, Zang L, Wang C (2013) Supported single-atom catalysts: synthesis, characterization, properties, and applications. Environ Chem Lett 2017: 1\u0026ndash;29.\u003c/p\u003e\n\u003cp\u003e[10] Kang L, Cheng B, Zhu M (2019) Pd/MCM-41 catalyst for acetylene hydrogenation to ethylene. R Soc Open Sci 6: 191155-191166.\u003c/p\u003e\n\u003cp\u003e[11] Panpranot J, Kontapakdee K, Praserthdam P (2006) Selective hydrogenation of acetylene in excess ethylene on micron-sized and nanocrystalline TiO\u003csub\u003e2\u003c/sub\u003e supported Pd catalysts. Appl Cat A Gen 314: 128\u0026ndash;133.\u003c/p\u003e\n\u003cp\u003e[12] Panpranot J, Kontapakdee K, Praserthdam P (2006) Effect of TiO\u003csub\u003e2\u003c/sub\u003e crystalline phase composition on the physicochemical and catalytic properties of Pd/TiO\u003csub\u003e2\u003c/sub\u003e in selective acetylene hydrogenation. J Phys Chem B 110: 8019\u0026ndash;8024.\u003c/p\u003e\n\u003cp\u003e[13] Meng L D, Wang G C (2014) A DFT + U study of acetylene selective hydrogenation over anatase supported Pd\u003csub\u003ea\u003c/sub\u003eAg\u003csub\u003eb\u003c/sub\u003e (a + b = 4) cluster. Phys Chem Chem Phys 16: 17541-17551.\u003c/p\u003e\n\u003cp\u003e[14] Yang J, Cao L X, Wang G C (2012) Acetylene hydrogenation on anatase TiO\u003csub\u003e2 \u003c/sub\u003e(101) supported Pd\u003csub\u003e4\u003c/sub\u003e cluster: oxygen deficiency effect. J Mol Model 18: 3329\u0026ndash;3339.\u003c/p\u003e\n\u003cp\u003e[15] Ma L L, Lv C Q, Wang G C (2017) A DFT study and micro-kinetic analysis of acetylene selective hydrogenation on Pd-doped Cu (111) surfaces. Appl Surf Sci 410: 154\u0026ndash;165.\u003c/p\u003e\n\u003cp\u003e[16] Yang J, Lv C Q, Guo Y, Wang G C (2012) A DFT+U study of acetylene selective hydrogenation on oxygen defective anatase (101) and rutile (110) TiO\u003csub\u003e2\u003c/sub\u003e supported Pd\u003csub\u003e4\u003c/sub\u003e cluster. J Chem Phys 136: 104107-104122.\u003c/p\u003e\n\u003cp\u003e[17] Li J N, Pu M, Ma C C, Tian Y, He J, Evans D J (2012) The effect of palladium clusters (Pdn, n = 2\u0026ndash;8) on mechanisms of acetylene hydrogenation: A DFT study. J Mol Catal A: Chem 359: 14\u0026ndash;20.\u003c/p\u003e\n\u003cp\u003e[18] Boyukata M, Belchior J C (2008) Molecular dynamics study of palladium clusters: size dependent analysis of structural stabilities and energetics of Pdn (n\u0026lt;40) via a lennard-jones type potential. Croatica Chem Act 81: 289-297.\u003c/p\u003e\n\u003cp\u003e[19] Moc J, Musaev D J, Morukama D (2003) Activation and adsorption of multiple H\u003csub\u003e2\u003c/sub\u003emolecules on a Pd\u003csub\u003e5\u003c/sub\u003ecluster: a density functional study. J Phys Chem A 107: 4929-4939.\u003c/p\u003e\n\u003cp\u003e[20] Zhou C, Yao S, Wu J, Forrey R C, Chen. L, Tachibana A, Chen H (2008) Hydrogen dissociative chemisorption and desorption on saturated sub nano palladium clusters (Pdn, n = 2\u0026ndash;9). Phys Chem Chem Phys 10: 5445\u0026ndash;5451.\u003c/p\u003e\n\u003cp\u003e[21] Cao Y, Sui Z, Zhu Y, Zhou X, Chen D (2017) Selective hydrogenation of acetylene over Pd-In/Al\u003csub\u003e2\u003c/sub\u003eO\u003csub\u003e3\u003c/sub\u003ecatalyst: promotional effect of indium and composition-dependent performance. ACS Catal 7: 7835\u0026ndash;7846.\u003c/p\u003e\n\u003cp\u003e[22] Albani D, Shahrokhi M, Chen Z, Mitchell S, Hauert R, Lopez N, P\u0026eacute;rez-Ram\u0026iacute;rez J (2018) Selective ensembles in supported palladium sulfide nanoparticles for alkyne semi-hydrogenation. J Nat Commun 9: 2634-2644.\u003c/p\u003e\n\u003cp\u003e[23] Lopez N, Vargas-Fuentes C (2012) Promoters in the hydrogenation of alkynes in mixtures: insights from density functional theory. Chem Commun 48: 1379\u0026ndash;1391.\u003c/p\u003e\n\u003cp\u003e[24] Volpe M A, Rodriguez P, Gigola C E (1999) Preparation of Pd\u0026ndash;Pb/\u0026alpha;-Al\u003csub\u003e2\u003c/sub\u003eO\u003csub\u003e3\u003c/sub\u003e catalysts for selective hydrogenation using PbBu\u003csub\u003e4\u003c/sub\u003e: the role of metal-support boundary atoms and the formation of a stable surface complex. Catal Lett 61: 27\u0026ndash;32.\u003c/p\u003e\n\u003cp\u003e[25] Anderson J A, Mellor J, Wells R P K (2009) Pd catalyzed hexyne hydrogenation modified by Bi and by Pb. J Catal 261: 208\u0026ndash;216.\u003c/p\u003e\n\u003cp\u003e[26] Kim S K, Lee L H, Ahn I K, Kim W J, Moon S H (2011) Performance of Cu-promoted Pd catalysts prepared by adding Cu using a surface redox method in acetylene hydrogenation. ApplCatalA Gen 401: 12\u0026ndash;19.\u003c/p\u003e\n\u003cp\u003e[27] Gonzalez S, Neyman K M, Shaikhutdinov S, Freund H J, Illas F (2007) On the promoting role of Ag in selective hydrogenation reactions over Pd\u0026ndash;Ag bimetallic catalysts: a theoretical study. J Phys Chem C 111: 6852\u0026ndash;6856.\u003c/p\u003e\n\u003cp\u003e[28] Mei D, Neurock M, Smith C M (2009) Hydrogenation of acetylene\u0026ndash;ethylene mixtures over Pd and Pd\u0026ndash;Ag alloys: first-principle based kinetic Monte Carlo simulations. J Catal 268: 181\u0026ndash;195.\u003c/p\u003e\n\u003cp\u003e[29] Sheth P A, Neurock M, Smith C M (2009) First-principles analysis of the effects of alloying Pd with Ag for the catalytic hydrogenation of acetylene\u0026ndash;ethylene mixtures. J Phys Chem B 109: 12449\u0026ndash;12466.\u003c/p\u003e\n\u003cp\u003e[30] Kumar N, Ghosh P (2016) Selectivity and reactivity of Pd-rich Pd-Ga surfaces toward selective hydrogenation of acetylene: interplay of surface roughness and ensemble effect. J Phys Chem C 120: 28654\u0026ndash;28663.\u003c/p\u003e\n\u003cp\u003e[31] Osswald J, Giedigkeit R, Jentoft R E, Armbr\u0026uuml;ster M, Girgsdies F, Kovnir K (2008) Palladium\u0026ndash;gallium intermetallic compounds for the selective hydrogenation of acetylene. Part I. Preparation and structural investigation under reaction conditions. J Catal 258: 210\u0026ndash;218.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;[32] Zhou H, Yang X, Wang A, Miao S, Liu X, Pan X (2016) Pd/ZnO catalysts with different origins for high chemoselectivity in acetylene semi-hydrogenation. Chin J Catal 37: 692\u0026ndash;699.\u003c/p\u003e\n\u003cp\u003e[33] Gulyaeva Y K, Kaichev V V, Zaikovskii V I, Kovalyov E V, Suknev A P, Balzhinimaev B S (2015) Selective hydrogenation of acetylene over novel Pd/fiberglass catalysts. Catal Today 245: 139\u0026ndash;146.\u003c/p\u003e\n\u003cp\u003e[34] Yan H, Cheng H, Yi H, Lin Y, Yao T, Wang T (2015) Single-atom Pd\u003csub\u003e1\u003c/sub\u003e/Graphene catalyst achieved by atomic layer deposition: remarkable performance in selective hydrogenation of 1,3-butadiene. J Am Chem Soc 137: 10484\u0026ndash;10487.\u003c/p\u003e\n\u003cp\u003e[35] Komhom S, Mekasuwandumrong O, Praserthdam P, Panpranot J (2008) Improvement of Pd/Al\u003csub\u003e2\u003c/sub\u003eO\u003csub\u003e3\u003c/sub\u003e catalyst performance in selective acetylene hydrogenation using mixed phases Al\u003csub\u003e2\u003c/sub\u003eO\u003csub\u003e3\u003c/sub\u003e support. Catal Commun 10: 86\u0026ndash;91.\u003c/p\u003e\n\u003cp\u003e[36] Hosseini S M, Ghiaci M, Kulinich S A, Wunderlich W, Farrokhpour H, Saraji M, Shahvar A (2018) Au-Pd@g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e as an efficient photocatalyst for visible-light oxidation of benzene to phenol: experimental and mechanistic study. J Phys Chem C 122: 27477\u0026ndash;27485.\u003c/p\u003e\n\u003cp\u003e[37] Shahzeydi A, Ghiaci M, Farrokhpour H, Shahvar A, Sun M, Saraji M (2019) Facile and green synthesis of copper nanoparticles loaded on the amorphous carbon nitride for the oxidation of cyclohexane. Chem Eng J 370: 1310-1321.\u003c/p\u003e\n\u003cp\u003e[38] Nilforoushan S, Ghiaci M, Hosseini S M, Laurent S, Muller R N (2019) Selective liquid phase oxidation of ethyl benzene to acetophenone by palladium nanoparticles immobilized on a g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e\u0026ndash;rGO composite as a recyclable catalyst. New J Chem 46: 6921-6931.\u003c/p\u003e\n\u003cp\u003e[39] Zhao Y, Zhu M, Kang L (2018) The DFT study of single-atom Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003ecatalyst for selective acetylene hydrogenation reaction. Cattal Lett 148: 2992\u0026ndash;3002.\u003c/p\u003e\n\u003cp\u003e[40] Kang L, Zhu M, Zhao Y (2019) A DFT study of acetylene hydrogenation catalyzed by S-doped Pd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e. Catalysts 9: 887-899.\u003c/p\u003e\n\u003cp\u003e[41] Hosseini S M, Ghiaci M, Farrokhpour H (2019) The adsorption of small size Pd clusters on a g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e quantum dot: DFT and TD-DFT study. Mater Res Express 6: 105079-105093.\u003c/p\u003e\n\u003cp\u003e[42] Singh N B, Sarkar U (2014) Structure, vibrational, and optical properties of platinum cluster: a density functional theory approach. J Mol Model 20: 2537\u0026ndash;2548.\u003c/p\u003e\n\u003cp\u003e[43] Lu T, Chen F (2012) Multiwfn: a multifunctional wavefunction analyzer. J Comput Chem 33: 580\u0026ndash;592.\u003c/p\u003e\n\u003cp\u003e[44] Frisch M J (2013) Gaussian 09, Revision D.01, Wallingford, CT: Gaussian, Inc.\u003c/p\u003e\n\u003cp\u003e[45] Mei D, Sheth P A, Neurock M, Smith C M (2006) First-principles based kinetic Monte Carlo simulation of the selective hydrogenation of acetylene over Pd (111). J Catal 242: 1\u0026ndash;15.\u003c/p\u003e\n\u003cp\u003e[46] Krajci M, Hafner J (2011) Complex intermetallic compounds as selective hydrogenation catalysts\u0026mdash;a case study for the (100) surface of Al\u003csub\u003e13\u003c/sub\u003eCo\u003csub\u003e4\u003c/sub\u003e. J Catal 278: 200\u0026ndash;207.\u003c/p\u003e\n\u003cp\u003e[47] Krajci M, Hafner J (2012) Intermetallic compound AlPd as a selective hydrogenation catalyst: a DFT study. J Phys Chem C 116: 6307\u0026ndash;6319.\u003c/p\u003e"},{"header":"Table","content":"\u003cp\u003e\u003cstrong\u003eTable 1. The energies of the HOMO and LUMO of C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e, H\u003csub\u003e2\u003c/sub\u003e, g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e, Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e, and Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e\u0026middot;C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2 \u003c/sub\u003ecompounds and their corresponding energy gaps\u003c/strong\u003e\u003c/p\u003e\n\u003ctable border=\"1\" width=\"0\"\u003e\n\u003ctbody\u003e\n\u003ctr\u003e\n\u003ctd rowspan=\"2\" width=\"138\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd rowspan=\"2\" width=\"81\"\u003e\n\u003cp\u003eHOMO\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd rowspan=\"2\" width=\"58\"\u003e\n\u003cp\u003eLUMO\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd rowspan=\"2\" width=\"59\"\u003e\n\u003cp\u003e∆E\u003csub\u003eg\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"0\"\u003e\u0026nbsp;\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"0\"\u003e\u0026nbsp;\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"138\"\u003e\n\u003cp\u003eH\u003csub\u003e2\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"81\"\u003e\n\u003cp\u003e-11.74 (eV)\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"58\"\u003e\n\u003cp\u003e1.84\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"59\"\u003e\n\u003cp\u003e-9.89\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"0\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"138\"\u003e\n\u003cp\u003eC\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"81\"\u003e\n\u003cp\u003e-8.25\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"58\"\u003e\n\u003cp\u003e0.45\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"59\"\u003e\n\u003cp\u003e-7.79\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"0\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"138\"\u003e\n\u003cp\u003eg-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"81\"\u003e\n\u003cp\u003e-6.53\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"58\"\u003e\n\u003cp\u003e-2.62\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"59\"\u003e\n\u003cp\u003e-3.90\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"0\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"138\"\u003e\n\u003cp\u003ePd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"81\"\u003e\n\u003cp\u003e-4.71\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"58\"\u003e\n\u003cp\u003e-2.71\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"59\"\u003e\n\u003cp\u003e-2.00\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"0\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"138\"\u003e\n\u003cp\u003ePd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e\u0026middot;C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"81\"\u003e\n\u003cp\u003e-4.09\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"58\"\u003e\n\u003cp\u003e-2.67\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"59\"\u003e\n\u003cp\u003e-1.41\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"0\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003e\u003cstrong\u003eTable 2. The imaginary frequencies of the transition states (TSs) involved in the different steps of the acetylene hydrogenation mechanism taken from Zhao\u0026rsquo;s group work [39] for the Pd1/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e\u0026nbsp;\u0026nbsp; with those obtained in this work.\u003c/strong\u003e\u003c/p\u003e\n\u003ctable border=\"1\"\u003e\n\u003ctbody\u003e\n\u003ctr\u003e\n\u003ctd colspan=\"5\" width=\"520\"\u003e\n\u003cp\u003eImaginary Frequency (cm\u003csup\u003e-1\u003c/sup\u003e)\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd rowspan=\"2\" width=\"104\"\u003e\n\u003cp\u003eCatalyst\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"85\"\u003e\n\u003cp\u003eIM5-TS5-P\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"104\"\u003e\n\u003cp\u003eIM4-TS4-IM5\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"104\"\u003e\n\u003cp\u003eIM2-TS3-IM3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"104\"\u003e\n\u003cp\u003eIM1-TS2-IM2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"123\"\u003e\n\u003cp\u003eR1-TS1-IM1\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"85\"\u003e\n\u003cp\u003e760.9\u003cem\u003ei\u003c/em\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"104\"\u003e\n\u003cp\u003e685.08\u003cem\u003ei\u003c/em\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"104\"\u003e\n\u003cp\u003e641.95\u003cem\u003ei\u003c/em\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"104\"\u003e\n\u003cp\u003e724.3\u003cem\u003ei\u003c/em\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"123\"\u003e\n\u003cp\u003e750.62\u003cem\u003ei\u003c/em\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"104\"\u003e\n\u003cp\u003ePd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"85\"\u003e\n\u003cp\u003e535\u003cem\u003ei\u003c/em\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"104\"\u003e\n\u003cp\u003e409\u003cem\u003ei\u003c/em\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"104\"\u003e\n\u003cp\u003e711\u003cem\u003ei\u003c/em\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"104\"\u003e\n\u003cp\u003e723.1\u003cem\u003ei\u003c/em\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"123\"\u003e\n\u003cp\u003e476\u003cem\u003ei\u003c/em\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"104\"\u003e\n\u003cp\u003ePd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003e\u003cstrong\u003eTable 3. Comparison of the ctivation energies for the first hydrogen transfer to the acetylene and ethylene reported in literature with those obtained in this work.\u003c/strong\u003e\u003c/p\u003e\n\u003ctable border=\"1\"\u003e\n\u003ctbody\u003e\n\u003ctr\u003e\n\u003ctd rowspan=\"2\" width=\"105\"\u003e\n\u003cp\u003eReference\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd rowspan=\"2\" width=\"93\"\u003e\n\u003cp\u003eSelectivity (∆E)\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd colspan=\"2\" width=\"283\"\u003e\n\u003cp\u003eEnergy barrier (kcal/mol)\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd rowspan=\"2\" width=\"94\"\u003e\n\u003cp\u003eCatalyst\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"137\"\u003e\n\u003cp\u003eC\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e4\u003c/sub\u003e + H\u003csub\u003e2\u003c/sub\u003e \u0026rarr; C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e6\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"146\"\u003e\n\u003cp\u003eC\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e2\u003c/sub\u003e + H\u003csub\u003e2\u003c/sub\u003e \u0026rarr; C\u003csub\u003e2\u003c/sub\u003eH\u003csub\u003e4\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e45\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"93\"\u003e\n\u003cp\u003e2.45\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"137\"\u003e\n\u003cp\u003e17.61\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"146\"\u003e\n\u003cp\u003e18.1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"94\"\u003e\n\u003cp\u003ePd(111)\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e46\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"93\"\u003e\n\u003cp\u003e6.34\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"137\"\u003e\n\u003cp\u003e19.57\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"146\"\u003e\n\u003cp\u003e15.41\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"94\"\u003e\n\u003cp\u003eAl\u003csub\u003e13\u003c/sub\u003eCoO\u003csub\u003e4\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e47\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"93\"\u003e\n\u003cp\u003e11.25\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"137\"\u003e\n\u003cp\u003e13.7\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"146\"\u003e\n\u003cp\u003e18.59\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"94\"\u003e\n\u003cp\u003eAlPd\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e39\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"93\"\u003e\n\u003cp\u003e15.4\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"137\"\u003e\n\u003cp\u003e38.88\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"146\"\u003e\n\u003cp\u003e21.53\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"94\"\u003e\n\u003cp\u003ePd\u003csub\u003e1\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003eThis Study\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"93\"\u003e\n\u003cp\u003e21.9\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"137\"\u003e\n\u003cp\u003e12.28\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"146\"\u003e\n\u003cp\u003e6.77\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"94\"\u003e\n\u003cp\u003ePd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003c/tbody\u003e\n\u003c/table\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"structural-chemistry","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"stuc","sideBox":"Learn more about [Structural Chemistry](https://www.springer.com/journal/11224)","snPcode":"11224","submissionUrl":"https://submission.nature.com/new-submission/11224/3","title":"Structural Chemistry","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"Acetylene, Hydrogenation, Pd cluster, g-C3N4, Selectivity, Ethylene","lastPublishedDoi":"10.21203/rs.3.rs-259183/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-259183/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eIn this work, the hydrogenation of acetylene on the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst is investigated by the Density Functional Theory (DFT) and Quantum Theory of Atoms in Molecules (QTAIM) calculations. The Pre-reactant (R), transition states (TSs), and the intermediates (IMs), involved in the hydrogenation process, are characterized from the point of view of energy and structure. The calculated energy barrier for the hydrogen transfer to the acetylene and ethylene are 6.77 and 12.28 kcal/mol, respectively which shows that the Pd\u003csub\u003e2\u003c/sub\u003e/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst has good selectivity for the conversion of acetylene to ethylene rather than ethane. Comparing the values of these energy barriers with those of the hydrogenation of acetylene on the Pd/g-C\u003csub\u003e3\u003c/sub\u003eN\u003csub\u003e4\u003c/sub\u003e catalyst (21.53 and 38.88 kcal/mol, respectively) shows that the increase in the number of the Pd atoms decreases the energy barriers of the hydrogenation reaction and increases the selectivity of the catalyst for the ethylene production.\u003c/p\u003e","manuscriptTitle":"Mechanistic Insight into the Hydrogenation of Acetylene on the Pd2/g-C3N4 Catalyst: Effect of Pd Clustering on the Barrier Energy and Selectivity","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2021-03-11 01:11:08","doi":"10.21203/rs.3.rs-259183/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"reviewersInvited","content":"","date":"2021-03-04T00:00:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2021-03-03T00:00:00+00:00","index":"","fulltext":""},{"type":"submitted","content":"Structural Chemistry","date":"2021-02-19T07:29:19+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"structural-chemistry","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"stuc","sideBox":"Learn more about [Structural Chemistry](https://www.springer.com/journal/11224)","snPcode":"11224","submissionUrl":"https://submission.nature.com/new-submission/11224/3","title":"Structural Chemistry","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"87f06d06-9d98-43fc-afc8-573773e0215e","owner":[],"postedDate":"March 11th, 2021","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":2873289,"name":"General Biochemistry"}],"tags":[],"updatedAt":"2021-08-18T19:26:51+00:00","versionOfRecord":{"articleIdentity":"rs-259183","link":"https://doi.org/10.1007/s11224-021-01781-3","journal":{"identity":"structural-chemistry","isVorOnly":false,"title":"Structural Chemistry"},"publishedOn":"2021-04-19 19:03:36","publishedOnDateReadable":"April 19th, 2021"},"versionCreatedAt":"2021-03-11 01:11:08","video":"","vorDoi":"10.1007/s11224-021-01781-3","vorDoiUrl":"https://doi.org/10.1007/s11224-021-01781-3","workflowStages":[]},"version":"v1","identity":"rs-259183","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-259183","identity":"rs-259183","version":["v1"]},"buildId":"cBFmMYwuxLRRLfASyISRj","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}