Next Article in Journal
A Critical Overview of ASP and Future Perspectives of NASP in EOR of Hydrocarbon Reservoirs: Potential Application, Prospects, Challenges and Governing Mechanisms
Next Article in Special Issue
Two-Dimensional (2D) TM-Tetrahydroxyquinone Metal–Organic Framework for Selective CO2 Electrocatalysis: A DFT Investigation
Previous Article in Journal
Influence of Cocatalysts (Ni, Co, and Cu) and Synthesis Method on the Photocatalytic Activity of Exfoliated Graphitic Carbon Nitride for Hydrogen Production
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two-Dimensional Transition Metal-Hexaaminobenzene Monolayer Single-Atom Catalyst for Electrocatalytic Carbon Dioxide Reduction

1
Institute for Advanced Study, School of Physics and Materials Science, Nanchang University, Nanchang 330031, China
2
School of Chemistry and Chemical Engineering, Nanchang University, Nanchang 330031, China
3
School of Mechanical and Electrical Engineering, Xinyu University, Xinyu 338004, China
4
School of Mathematical Sciences and Physics, Jinggangshan University, Ji’an 343009, China
5
Materials Genome Institute, National Center for International Research on Photoelectric and Energy Materials, School of Materials and Energy, Yunnan University, Kunming 650091, China
6
Advanced Computing Center, Information Technology Center, Yunnan University, Kunming 650091, China
*
Authors to whom correspondence should be addressed.
Nanomaterials 2022, 12(22), 4005; https://doi.org/10.3390/nano12224005
Submission received: 15 October 2022 / Revised: 7 November 2022 / Accepted: 11 November 2022 / Published: 14 November 2022
(This article belongs to the Special Issue 1D and 2D Nanomaterials for Energy Storage and Conversion)

Abstract

:
Electrocatalytic reduction of CO2 to valuable fuels and chemicals can not only alleviate the energy crisis but also improve the atmospheric environment. The key is to develop electrocatalysts that are extremely stable, efficient, selective, and reasonably priced. In this study, spin-polarized density function theory (DFT) calculations were used to comprehensively examine the catalytic efficacy of transition metal-hexaaminobenzene (TM-HAB) monolayers as single-atom catalysts for the electroreduction of CO2. In the modified two-dimensional TM-HAB monolayer, our findings demonstrate that the binding of individual metal atoms to HAB can be strong enough for the atoms to be evenly disseminated and immobilized. In light of the conflicting hydrogen evolution processes, TM-HAB effectively inhibits hydrogen evolution. CH4 dominates the reduction byproducts of Sc, Ti, V, Cr, and Cu. HCOOH makes up the majority of Zn’s reduction products. Co’s primary reduction products are CH3OH and CH4, whereas Mn and Fe’s primary reduction products are HCHO, CH3OH, and CH4. Among these, the Ti-HAB reduction products have a 1.14 eV limiting potential and a 1.31 V overpotential. The other monolayers have relatively low overpotentials between 0.01 V and 0.7 V; therefore, we predict that TM-HAB monolayers will exhibit strong catalytic activity in the electrocatalytic reduction of CO2, making them promising electrocatalysts for CO2 reduction.

1. Introduction

In recent years, a large amount of CO2 has been released into the atmosphere, leading to a series of environmental and socio-economic problems. Converting CO2 into hydrocarbon fuels and valuable chemical raw materials such as methane, methanol, and formic acid is important for resource utilization and environmental protection, as well as for global sustainable development [1,2,3,4,5,6,7]. However, it is not easy to achieve efficient electrocatalytic CO2 reduction, which is a challenge because the carbon atoms in the CO2 molecule adopt sp hybridization to bond with two oxygen atoms to form C=O bonds, which endows the CO2 molecule with high chemical inertness and low solubility in aqueous solutions. Currently, CO2 conversion can be carried out using different catalytic approaches involving chemical [8], photochemical [9], electrochemical [10], and biological methods [11]. Among the various methods, electrochemical methods have been widely noted for their mild conditions, simple operation, safety and stability, and easy access to intermittent energy sources such as solar and wind power [12]. In the electrocatalytic CO2 reduction process, CO2 molecules are activated by the interaction with the catalyst surface under the action of external electrical energy, and the subsequent reaction steps can be carried out further under mild conditions with relatively low energy cost. However, CO2 activation by high-activity electrocatalysts, the products of CO2RR, generally occurs with low selectivity. Therefore, it is quite urgent to design and fabricate catalysts with high activity and distinct selectivity. It is well known that effective activation of carbon dioxide requires electron transfer from the substrate to the carbon dioxide molecule, which leads to appropriate accompanying structural deformation [13,14,15,16,17,18]. In the past few years, various heterogeneous catalysts for the activation of inert CO2 have been proposed, including pure metals [19], metal oxide interfaces [20], graphite-based materials [21], sulfides [22], and metal organic frameworks (MOFs) [23]. To meet the high activity and selectivity for specific products, researchers have combined the aforementioned heterogeneous catalysts with the emerging nanoscience for functional applications [24]. Supported metal nanoparticles have been a widespread approach in recent years, by which to reshape heterogeneous catalysts to meet different needs. The size of metal particles is an important factor affecting catalytic performance. As the size of metal particles decreases, their surface area/volume ratio and catalytic activity can be significantly improved. Rationally designed single-atom catalysts (SACs) for CO2 conversion show high activity and significant selectivity. The reduction of CO2 to CH3OH is catalyzed by β 12 boron monolayers supported by V atoms [25], as well as the conversion of CO2 to CH4 on the surface of SACs fabricated from porphyrin-like graphene supported by single Co, Rh, and Ir atoms [26]. In addition to the type of metal atoms supported on the surface of nanomaterials, organometallic complexes anchored on substances, as well as anchored on porous materials involving MOFs, zeolites, and other ion-exchange metals on the surface, can also be considered as SACs [27,28]. Thus, MOFs that include low-coordination metal centers individually and regularly dispersed on the surface are promising for a wide range of SAC applications in the near future [29,30,31,32,33,34,35,36,37,38,39,40,41,42].
Two-dimensional MOFs may be directly synthesized without the need for material modification to become highly useful SACs. Additionally, a variety of metals may be used for the metal centers of two-dimensional MOFs, which always serve as active sites and impact the performance of the catalyst. By changing the kind of metal center to suit our demands, we may make use of the tunability of 2D MOFs. Additionally, MOFs possess metallic characteristics that guarantee quick electron transfer during electrochemical catalysis, a crucial quality of effective electrocatalysts. The formation of TMN4 complexes, which have high catalytic characteristics comparable to those of noble metals in a variety of catalytic interactions such oxygen reduction, nitrogen fixation, and CO2 reduction, is noteworthy [43,44,45,46]. Recently, Louie’s group successfully fabricated a series of novel two-dimensional MOFs using a bottom-up technique in which the transition metal (TM) atoms serve as the metal center and hexamethylbenzene (HAB) serves as the ligand [47]. Low-coordination TM atoms as metal centers are regularly and stably immobilized to form MOFs with the molecular formula TM3(HAB) 2 . The low-coordinated TM atoms are stably anchored as metal centers in the MOF TM3(HAB) 2 , exhibiting the characteristics of single-atom catalysts (SACs). Therefore, they can be considered SACs with high practicality. The metal center in the two-dimensional TM-HAB plays the role of an active site with catalytic properties, and the type of central atom can be tuned to meet the catalytic requirements [29]. Notably, TM-HAB is formed by each TM atom with four surrounding N atoms. The TMN4 complex is an analogue of TMNx, which exhibits good catalytic properties similar to those of noble metals in various catalysts, such as in oxygen reduction, nitrogen fixation, and carbon dioxide reduction [43,44,45,46]. Park et al. [48] proposed a two-dimensional (2D) conductive metal-organic framework consisting of M-N4 units (M = Ni, Cu) and hexaaminobenzene (HAB) linkers as a catalyst for oxygen reduction reactions, and the results showed that the catalytic performance depends strongly on the metal species. However, the application of TM-HAB monolayers for CO2 reduction has been little reported so far. This motivated us to explore whether TM-HAB constructed with different metal species could be used as a prospective electrocatalyst for CO2 reduction. In this paper, the catalytic properties of the first periodic transition-metal series TM-HAB monolayers for CO2 were systematically investigated using spin-polarized density function theory (DFT). Our results show that most TM-HAB monolayers exhibit excellent catalytic activity and strong stability. Therefore, we predict that TM-HAB monolayers can contribute to the next generation of low-cost, high-stability electroreduction catalysts.

2. Computational Methods

The Dmol3 package’s spin-polarized density functional theory (DFT) [49] was used for all computations in this work. The Perdew-Burke-Ernzerhof (PBE) function of the generalized gradient approximation (GGA) was used to accomplish the exchange correlation of electrons [50]. Through Grimme’s technique, which has been employed for different gas adsorptions and gas phase catalysts, the van der Waals (vdW) force was used to deal with long-range dispersion correction [51,52,53,54]. For transition metals, the density functional semicore pseudopotential (DSPP) was used to properly calculate other electrons and substitute core electrons with an effective pseudopotential [55]. The double numerical (DN) function was selected as the basis set for other atoms. The periodic supercell of each MOF contains 12 hydrogen atoms, 12 nitrogen atoms, 12 carbon atoms, and 3 transition metal atoms, and a unit cell of 14.33 × 14.33 × 25 Å3 with a vacuum space of 25 Å along the Z axis, which is spacious enough to avoid nonphysical interaction with the periodic image. The Brillouin zone (BZ) k-point sample was a 5 × 5 × 1 grid with a Monkhorst–Pack design [56]. The convergence thresholds for energy, gradient, and displacement for all geometric optimizations were 1.0 × 10−5 Ha, 2.0 × 10−3 Ha/Å, and 5.0 × 10−3 Å, respectively. The conductor-like screening model (COSMO) was used to simulate the electrolyte, which is always an aqueous environment, with a dielectric constant of 78.54 for H2O as the solvent [57]. The sorption energy (Eads) of CO2 on the TM-HAB periodical unilayer is defined as:
E ads = E TM - HAB - CO 2 E TM - HAB E CO 2
where E TM - HAB - CO 2 , E TM - HAB , and E CO 2 are the total energy for the CO2 adsorbed on the TM-HAB monolayer, the pure TM-HAB monolayer, and the CO2 molecule, respectively. The computational hydrogen electrode (CHE) model proposed by Nrskov and coworkers to account for the energy of a proton–electron pair in aqueous solution was used to estimate the change in Gibbs free energy (ΔG) for each of the CO2RR steps [58,59,60]. Gibbs free energy is defined as:
Δ G = Δ E + Δ E ZEP T Δ S + Δ G pH + Δ G U
where ΔE is the response energy, which can be obtained directly from the DFT calculation. ΔEZPE is the zero-point change in energy, ΔS is the change in the entropy, and T is the systematic temperature (298.15 K). EZPE and S of the CO2RR interstitial are calculated from the vibrational frequency. Δ G pH = 2.303 k BT   pH is a modification of the free energy due to the variation in H+ concentration. pH is taken to be zero for acidic mediators in this article. ΔU = −neU, where n is the number of transferred electrons, e is the electronic charge and U is the applied voltage. The limiting potential (UL) of the CO2RR can be obtained from the free potential change ( Δ G Max ) by the relation U L = Δ G Max / ne . The overpotential ( η ) is evaluated as the difference between the balanced potential and the limiting potential. Therefore, the overpotential is defined as:
η = U equilibrium U L

3. Results and Discussion

3.1. Structural Features and Properties of the TM-HAB Monolayer

The optimized architectures of the TM-HAB materials are depicted in Figure 1. Table S1 lists specific lattice constant values. Each MOF has three transitional metal atoms, 12 hydrogen atoms, 12 nitrogen atoms, and 12 carbon atoms in each periodic cell. The lattice parameters of each cell are shown in Table S1. Figure 1 displays models of a periodic 2 × 2 supercell. Each transition metal atom binds to two HAB molecules in the TM-HAB monolayer. The atoms in each of the 10 transition metals that we took into consideration (from Sc to Zn) all reside in the same plane (Figure 1). The bond length decreases from Sc to Co, the metal-N bond length from 1.837–2.110 Å, and after that, the metal-N bond length increments to a Zn-N length of 2.032 Å. Two of the monolayers (V and Mn), according to our calculations, exhibit differing degrees of spin polarization in their ground states, which lowers their energy below that of the non-magnetic state. The magnetic moment with the highest value is 3.22 mB for Mn, as indicated in Table 1. Additionally, we studied the Hirshfeld charge of these monolayers. As can be seen in Table 1, the nitrogen atoms have a small negative charge, while the 10 metal atoms have a small positive charge. Ionic and covalent bonds exist between the metal atoms and the atoms in their vicinity.

3.2. Stabilization of TM-HAB Monolayer

The formation energy represents the complexity of the catalyst preparation process. If the formation energy is negative, the preparation process will be exothermic. As a result, the lower the formation energy, the easier the preparation and the more stable the material. Figure 2 demonstrates that the formation energies of these materials are all negative, ranging from 2.59 to 7.88 eV (for further information, see Table S2). This suggests that the ten TM-HAB single-atom catalysts we investigated should be quite simple to manufacture experimentally.
The stability of these materials is a key feature to consider when evaluating their catalytic performance. The higher the binding energy of the transition metal atom (TM) to the substrate (HAB) in TM-HAB catalysts, the better the catalytic stability. If the binding energy is too low, the metal atoms may agglomerate into clusters, reducing the catalytic effectiveness of the single-atom catalyst. According to Figure 2, the binding energies of the ten 3d transition metal elements with HAB are all negative, with Ti-HAB having the greatest negative binding energy of 13.31 eV and Zn-HAB having a lower negative binding energy of 5.38 eV. These values are also mentioned in Table S2.
Additionally, we looked at the cohesion energy of bulk metals. Table S2 demonstrates that the cohesion energy of bulk metals, which ranges from 1.05 to 6.58 eV for Sc to Zn, is negative. If the metal atoms are likely to form clusters or are stably implanted in the skeletal layer as single atoms, it may be determined by comparing the cohesion energy and binding energy. The ten 3d transition elements from Sc to Zn that are bound to HAB have binding energies that are all lower than the cohesion energy, as shown in Figure 2. This suggests that the metal atoms can be firmly buried in the HAB monolayer as active sites.

3.3. The First Hydrogenation Step: Selectivity for CO2RR vs. HER

Each stage of the whole electrocatalytic CO2 reduction process calls for the participation of a proton–electron pair ( H + + e ). Depending on the location of the H addition, two distinct intermediates are produced during the initial protonation of the CO2 reduction. Oxygen atoms may be transformed into the intermediate *COOH by adding H to them. However, the intermediate *OCHO will be created if H is added to the carbon atom. It is also possible to combine H with metal atoms to create adsorbed *H, which will cause the hydrogen evolution process. Due to the competing nature of the two CO2RR and HER responses, this final instance is not what we would want to observe for the CO2RR. Therefore, we must take into account the CO2RR catalysts’ ability to prevent hydrogen evolution. These steps are detailed in the equations below:
* + C O 2 + H + + e * C O O H * + C O 2 + H + + e * O C H O * + H + + e * H
Figure 3 compares the Gibbs free energy changes for the first protonation reaction steps for the formation of *OCHO, *COOH, and *H. The detailed values are shown in Table S3. As shown in Figure 3, for the 10 transition metals, the Gibbs free energy change for the formation of *OCHO is lower than that for the formation of *COOH intermediates for the other 9 metals except metal V, indicating that the further protonation of these 9 catalysts tends to produce *OCHO intermediates more after the activation of the adsorbed CO2 molecules. If Δ G [ * COOH ] or Δ G [ * OCHO ] are smaller than Δ G [ * H ] , it is easier to form *COOH or *OCHO than H*. Once the active site is occupied by *COOH or *OCHO, there are few remaining active sites to form *H, so the hydrogen evolution reaction is inhibited. Therefore, the catalysts above the dashed line in Figure 3 will be dominated by the hydrogen evolution reaction, while the catalysts below the dashed line will be dominated by the desired CO2 reduction reaction. As can be seen from Figure 3, the Gibbs free energy change for the formation of *COOH or *OCHO is lower than that for the formation of *H. All 10 metal values are below the dashed line, so the TM-HAB monolayer we studied will have strong hydrogen evolution reaction inhibition properties.

3.4. Possible Product Pathways and Adsorption Energy

Since TM-HAB electrocatalytic CO2 reduction is a single-atom catalytic process, it is generally believed that it is difficult to generate multi-carbon products because the single-atom catalytic process cannot achieve C-C coupling between intermediates. Therefore, theoretically, it is sufficient to predict the monoatomic catalytic CO2 reduction process by considering the C1 product. The most common C1 products of CO2 electrocatalytic reduction are CO, CH4, HCOOH, CH3OH, and HCHO. Figure 4 shows the scheme for the electrocatalytic reduction of CO2 to obtain the C1 product [61].
As can be seen in Figure 4, the reduction of CO2 produces CO and HCOOH as a 2e process. The reduction paths are *CO2 → *COOH → *CO → CO and *CO2 → *OCHO → *HCOOH → HCOOH. The generation of HCHO is a 4e process and the reduction path is *CO2 → *COOH → *CO → *CHO → *OCH2 → HCHO. The obtaining of CH3OH product is a 6e process and the reduction path is *CO2 → *COOH → *CO → *CHO → *OCH2 → *OCH3 → *OHCH3 → CH3OH. The most complicated is the obtaining of CH4 product, which is an 8e process, and there are three possible paths, which are (1) *CO2 → *COOH → *CO → *CHO → *C → *CH → *CH2 → *CH3 → CH4; (2) *CO2 → *COOH → *CO → *CHO → *OCH2 → *OCH3 → *OHCH3 → *O + CH4 → *OH + CH4 → * + H2O + CH4; and (3) *CO2 → *COOH → *CO → *CHO → *OCH2 → *OCH3 → *OHCH3 → *OH + CH4 → * + H2O + CH4.
Based on the complexity of the CO2 electrocatalytic reduction reaction pathway, in order to predict the most likely products for each catalyst, we first calculated the adsorption energy of the catalyst for the C1 products as shown in Table 2.
Ni-HAB is relatively large for both species’ C1 product energy absorption, thus the product is firmly adsorbed by the catalyst and cannot be desorbed during the catalytic process, and the whole process is poisoned and no product can be obtained; this process is considered to be catalytically inactive. Similarly, for the three single-atom catalysts, Sc, Ti, and V, they have a strong adsorption capacity for CO, HCOOH, HCHO, and CH3OH, resulting in the inability to desorb and obtain the product; fortunately, this is not strong for CH4 adsorption, thus making it possible to obtain the CH4 product. The possible products obtained from Cr single atoms are HCHO, CH3OH, and CH4. Mn-HAB is strong enough to adsorb only CO; thus, the possible products of Mn-HAB are HCOOH, HCOH, CH3OH, and CH4. Similarly, the possible products of Co-HAB are only HCOOH, CH3OH, and CH4. However, the single-atom catalysts of Fe, Cu, and Zn are weak to adsorb the products, and all possible products can be obtained.

3.5. Reaction Pathways for CO2 Electrochemical Reduction

3.5.1. HCOOH as the Main Catalytic Product

MOF electrocatalytic CO2 reduction to produce a single HCOOH product is a fascinating thing. Janire et al. prepared zirconium-based MFO applied to electrocatalytic CO2 reduction to produce a single-product formic acid in the liquid fraction [62]. We calculated the catalytic process step diagram and found that the main product of Zn-HAB for electrocatalytic CO2 reduction is HCOOH. The free-energy step curve is shown in Figure 5. According to the scheme in Figure 4, after CO2 adsorbs Zn-HAB, the first step of protonation occurs under the action of external potential to generate *COOH or *OCHO intermediates, and it can be seen in Figure 5 that the generation of *OCHO intermediates is a Gibbs-free-energy drop process, and the reaction is easily carried out. In addition, the generation of *COOH intermediate needs to cross a very high energy barrier of 0.934 eV (see Table S4), so the first step of protonation to generate *OCHO is dominant, and the second step of paper protonation to generate *OCHOH intermediate that occurs on this basis is also a free-energy drop process. Then, after the formation of *OCHOH, the intermediate has the possibility of a *OCHOH + H+ + e → *CHO/*OCH + H2O protonation reaction, but it needs to cross the energy barrier of 1.314/1.292 eV, relatively, and the direct desorption of *OCHOH to form the HCOOH process only needs to cross the energy barrier of 0.24 eV. Detailed values of the free energy are given in Table S4. Therefore, this step is more inclined to be terminated by HCOOH desorption. Since the tendency to form *CHO/*OCH intermediates to compete is weak, the occurrence of multi-electron steps such as 3e and 4e proceeds with difficulty, thus making it difficult to obtain CH3OH, HCHO, or even CH4 products. In summary, the main product of Zn-HAB for electrocatalytic CO2 reduction is HCOOH. The reaction path is * + CO2 → *OCHO → *OCHOH → * + HCOOH. The rate-determining step is *OCHOH → * + HCOOH with a limiting potential of 0.24 V.

3.5.2. CH3OH and CH4 Are Produced Simultaneously as the Main Reduction Products

In Table 2, it is shown that the adsorption energy of Co-HAB on CO and HCHO is too large to desorb and obtain the product. Here, we analyzed the theoretical calculation of HCOOH, CH3OH, and CH4 as products and found that the main products obtained simultaneously by this catalyst were CH3OH and CH4. Figure 6 shows the whole process.
After CO2 adsorption by Co-HAB, the steps of paper protonation to form *COOH or *OCHO are all free energy reduction processes, which are exothermic and occur easily. Table S5 gives the detailed protonation steps and the Gibbs free energy change of each step for the electrocatalytic reduction of CO2 by Co-HAB. The further protonation of *COOH to *CO is also a free energy reduction process; however, the further protonation of *OCHO to *OCHOH requires a high external energy supply. Table S5 shows that 2.63 eV is required, so this step and the subsequent pathway need not be considered. In fact, *CO is further protonated to generate *CHO/*COH, but the *COH generation step needs to overcome the energy barrier of 1.73 eV (see Table S5), and this process need not be considered, so Figure 6 only gives the *CHO generation step and considers its subsequent protonation step. It is clear that further protonation of *CHO to generate *OCH2 requires overcoming an energy barrier of 0.53 eV, and the subsequent steps, whether CH3OH or CH4 intermediates, are either free energy reduction processes or require overcoming energies below 0.53 eV. Thus, CH3OH and CH4 can both be generated if 0.53 eV energy is obtained from outside.
In the 6e conversion process, the formation of *CH3OH from *OCH3 intermediates is a free energy reduction process, while the formation of *O intermediates requires external energy; thus, the reaction pathway prefers the formation of *CH3OH. The subsequent CH3OH desorption process needs to overcome an energy barrier of 0.465 eV (see Table S5), which is lower than the 0.53 eV of the *OCH2 formation step, and therefore the reaction pathway at 0.53 eV driven by external energy can occur smoothly. In addition, further protonation of *CH3OH to form *OH + CH4 is a free energy reduction process, and the next 8e process to generate H2O needs to cross the energy barrier of 0.095 eV. In summary, the main products of Co-HAB electrocatalytic CO2 reduction are CH4 and CH3OH, and the reaction path is * + CO2 → *COOH → *CO → *CHO → *OCH2 → *OCH3 → *O/*CH3OH → (* + CH3OH)/(*OH + CH4) → * + H2O + CH4. The rate-determining steps are *CHO + H2O + H+ + e → *OCH2 + H2O with a limiting potential of 0.53 V.

3.5.3. HCHO, CH3OH, and CH4 Are Produced Simultaneously as the Main Reduction Products

From the adsorption calculations, it is clear that the adsorption energy of Mn-HAB is too large only for CO to obtain the C1 product. The adsorption energy of Fe-HAB is not large for all the five C1 products. Our calculations show that HCHO, CH3OH, and CH4 can be obtained with both Mn-HAB and Fe-HAB catalysts. Figure 7 shows the free energy steps of electrocatalytic CO2 reduction by Mn-HAB and Fe-HAB.
For the electrocatalytic CO2 process with both Mn-HAB and Fe-HAB catalysts, the first step of the protonation process after the activation of CO2 by adsorption is a free energy reduction process, whether or not *COOH or *OCHO intermediate is generated. In the 2e conversion process, the generation of *CO or *OCHOH from Mn-HAB is a free energy reduction process (see Figure 7a), while the generation of *OCHOH intermediate from Fe-HAB needs to overcome the energy barrier of 0.51 eV (see Figure 7b and Table S7), and thus its subsequent protonation step need not be considered. In the subsequent 3e process, the lowest energy barrier step for both catalysts is *CO + H2O + H+ + e → *CHO + H2O, with 0.27 eV for Mn-HAB and 0.1 eV for Fe-HAB. However, in the subsequent 4e electron process, Fe-HAB also needs to overcome an energy barrier of 0.27 eV to produce the *OCH2 intermediate. In the further protonation process, Mn-HAB is almost always a free energy step-down process, although the final steps of HCHO, CH3OH desorption, and CH4 generation step up and need to overcome the energy barriers, but their values are 0.23 eV, 0.27 eV, and 0.14 eV, respectively (Table S6). For 0.1 eV energy, all three products are accessible to both catalysts at an external energy of 0.27 eV. Fe-HAB is also a free energy reduction process in the subsequent protonation process, except that the HCHO and CH3OH desorption steps need to cross energy barriers of 0.102 eV and 0.098 eV (Table S7), respectively. In conclusion, the Mn-HAB and Fe-HAB electrocatalytic CO2 reduction reactions are similar with the following pathways: * + CO2 → *COOH → *CO → *CHO → *OCH2 → *OCH3/(* + HCHO) → *O/*CH3OH → (* + CH3OH)/(*OH + CH4) → * + H2O + CH4. The limiting potential is 0.27 V. The difference is that the rate-determining step of Mn-HAB is *CO + H2O + H+ + e → *CHO + H2O, while the rate-determining step of Fe-HAB is *CHO + H2O + H+ + e → *OCH2 + H2O.

3.5.4. CH4 as the Main Catalytic Product

Since Sc-HAB (Table S9), Ti-HAB, and V-HAB only have a small adsorption energy for CH4, the other four C1 products will be poisoned during the catalytic process and no product will be obtained. Thus, they can only obtain a single product of CH4 for electrocatalytic CO2 reduction. In fact, the CO, HCOOH, HCHO, and CH3OH desorption steps in the Cu-HAB catalytic process need to overcome 0.86 eV, 1.22 eV, 0.31 eV, and 0.7 eV energy barriers, respectively (see Table S11), all of which are higher than the energy barriers required for the CH4 generation step, and thus CH4 is considered the most likely product to be obtained. In terms of adsorption energy, HCHO, CH3OH, and CH4 can all be desorbed on the Cr-HAB surface to obtain the products, but in the actual catalytic process, the HCHO and CH3OH desorption steps need to overcome 0.56 eV and 0.80 eV energy barriers, respectively (see Table S12), which are higher than the energy barriers required for the CH4 generation step, as shown in Figure 8d. Thus, CH4 is also considered to be the most likely product to be obtained.
Figure 9 shows the step diagram of Ti-HAB electrocatalytic CO2 reduction along the possible pathways of reduction, and the reaction equations and Gibbs free energy changes for each protonation step are shown in Table S8. It is clear that in the 3e process, the protonation from *CO to *COH step needs to overcome the energy barrier. Table S8 shows the value of 1.12 eV for this process, but the re-protonation of *COH to *C needs to overcome the 2.12 eV energy barrier (see Table S8). Although the free energy of the *OCHO intermediate is reduced and can be easily generated, the energy barriers required for the further protonation of *OCHO are all higher. Therefore, the reaction path for Ti-HAB electrocatalytic CO2 reduction is *CO2 → *COOH → *CO → *CHO → *OCH2 → *OCH3 → *O/*CH3OH → *OH + CH4 → * + H2O + CH4. The reaction rate-determining step is *OH + CH4 + H2O + H+ + e → * + CH4 + 2H2O with a limiting potential of 1.14 V.
Figure 8b shows the step diagram of V-HAB electrocatalytic CO2 reduction along the possible pathway reduction, and the reaction equations of each protonation step and Gibbs free energy change are shown in Table S10. Both the first and second protonation steps are free energy reduction processes, and the *CO/*OCHOH intermediate is easily obtained. However, in the further protonation process, only the *OCHOH → *OCH step requires the lowest energy barrier to be crossed, which is 0.27 eV, and this energy barrier is the highest energy barrier for CH4 production and thus the limiting potential. The energy barriers to be crossed in the *CO → *COH/*CHO step are 1.523 eV and 0.453 eV, respectively (Table S10), and the energy barrier to be crossed in the *OCHOH → *CHO step is 1.205 eV (Table S10), which are all higher than 0.27 eV and are therefore not considered. Therefore, the pathway of V-HAB electrocatalytic CO2 reduction is *CO2 → *OCHO → *OCHOH → *OCH → *OCH2 → *OCH3 → *O/*CH3OH → *OH + CH4 → * + H2O + CH4. The rate-determining step is *OCHOH + H+ + e → *OCH + H2O, corresponding to a limiting potential of 0.27 V.
Sc-HAB, Cu-HAB, and Cr-HAB are similar to Ti-HAB in that the first step of the protonation process to generate *COOH/*OCHO after CO2 activation by adsorption is exothermic, but the further protonation of *OCHOH to generate *CHO/*OCH needs to overcome a higher energy barrier, while the step of *CO intermediate to generate *COH also needs to overcome a higher energy barrier. Figure 8a,c,d represent the free energy changes of intermediates of the catalytic processes of Sc-HAB, Cu-HAB, and Cr-HAB, respectively. These three catalysts have the same reaction path within the 3e process as *CO2 → *COOH → *CO → *CHO. However, for the Cu-HAB catalyst, the energy barrier of 1.17 eV (see Table S11) is required to overcome the *OCH3 + H2O + H+ + e → *O + CH4 + H2O step in the 6th protonation step, and thus this reaction is not considered. Otherwise, the reaction paths of these three catalysts can be expressed as *CO2 → *COOH → *CO → *CHO → *OCH2 → *OCH3 → *CH3OH/*O (except Cu) → *OH + CH4 → * + H2O + CH4. Among them, the rate-determining step of the Sc-HAB catalytic process is *OCH2 + H2O + H+ + e → *OCH3 + H2O with a limiting potential of 0.29 eV; the rate-determining step of the Cu-HAB catalytic process is *CH3OH + H2O + H++ e → *OH + CH4 + H2O with a limiting potential of 0.18 eV; and the rate-determining step of the Cr-HAB catalytic process is *OH + CH4 + H2O + H+ + e → * + CH4 + 2 H2O with a limiting potential of 0.49 eV. The rate-determining step of the Cr-HAB process is *OH + CH4 + H2O + H+ + e → * + CH4 + 2 H2O with a limiting potential of 0.49 eV.

3.6. Electronic Structure

In Section 3.5, after the Gibbs free energy change for each intermediate step, we discussed the rate-determining step, the limiting potential, and the corresponding major products for each catalyst as shown in Table 3. In summary, all catalysts except Ni-HAB are electrocatalytically active for CO2 reduction, and the main products of Sc, Ti, V, Cr, and Cu are CH4. The limiting potential and overpotential of Ti-HAB are the highest, at 1.14 V and 1.31 V, respectively, and the overpotentials of the other monolayer catalysts are in the range 0.01–0.7 V, which are comparable to those of the most active step surface Cu(211) ( η = 0.77 V) and the overpotential of the most active metal surface Pt(111) ( η = 0.46 V) [63]. In addition, our results were compared with experimentally prepared catalysts. An example is the work of Xu’s team, who synthesized non-peripheral octamethyl-substituted cobalt(II) phthalocyanine (N-CoMePc) catalysts, which achieved a Faraday efficiency of up to 94.1% for CO production at a low overpotential of 0.6 V [64]. Another typical example is the work of Ivan et al. [65], which involves the synthesis, description, and preliminary evaluation of bimetallic copper-based hollow fiber electrodes with a compact three-dimensional geometry to overcome mass transfer limitations and improve the electrochemical conversion of CO. It is noted that the generation of CO occurs in the range 1–1.5 V. Some of our predicted TM-HABs have comparable or even lower overpotentials than N-CoMePc and bimetallic Cu-based hollow fiber. Thus, our theoretical results suggest a very promising single-atom catalyst for electrocatalytic CO2 reduction.
The metal-ligand bonding theory of organometallic catalysts reveals that the interactions between catalysts and intermediates are mainly σ - and π -bonds. Figure 10 shows a clear overlap between the 3d orbital of the metal atom and the 2p orbital of the O atom or C atom in the decisive step intermediates (*OCH2, O*H, O*CHOH, *CO, *CHO, or *OHCH3), either spin-up or spin-down, which indicates a good interaction between the TM-HAB and the intermediate. However, the overlap effect of 3d and 2p orbitals in Figure 10b is better than the overlap of 3d and 2p orbitals in Figure 10a,c–i. This indicates that the interaction of Ti-HAB with the corresponding intermediates is stronger than that of the other catalysts. In addition, it can be seen from Table 3 that the limiting potential UL of −1.14 eV for the electrocatalytic CO2 reduction by Ti-HAB is larger than that of the other catalysts. This is in good agreement with the results of PDOS. The stronger the interaction, the more stable the adsorption intermediate system, and the higher energy barrier that needs to be overcome to ensure that the catalytic reaction occurs, which leads to a larger increase in the free energy of the Ti-HAB catalyzed CO2 reduction decisive step, and thus a more negative limiting potential for the reduction reaction.

4. Conclusions

In summary, we investigated the electrocatalytic CO2 reduction reaction of single-atom catalysts created from transition metal-hexaaminobenzene two-dimensional coordination network materials. Density functional theory calculations show that for 10 transition metal TM-HAB monolayers ranging from Sc to Zn, the binding energy of the metal atoms to the HAB is large enough for the metal atoms to be stably dispersed in the HAB monolayer. All of these materials inhibit the hydrogen evolution reaction. The reduction products of Sc, Ti, V, Cr, and Cu are mainly CH4. The reduction products of Zn are mainly HCOOH. The main reduction products of Co are CH3OH and CH4. The reduction products of Mn and Fe are mainly HCHO, CH3OH, and CH4. Among them, the limiting potential of Ti-HAB is 1.14 eV and the overpotential is 1.31 V. The overpotentials of the other monolayers were in the range 0.01–0.7 V. All of the values are relatively low; therefore, we predict that a TM-HAB monolayer will exhibit strong catalytic activity in the electrocatalytic reduction of CO2, making it a promising electrocatalyst for CO2 reduction.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano12224005/s1, Table S1: Lattice constant a and b of TM-HAB in Å, where TM are the metal atoms of the first transition metal series; Table S2: Ec is the cohesive energy of the bulk TM, Eb is the binding energy between the TM and the TM-HAB, and Ef is the formation energy of TM-HAB, where TM are the metal atoms of the first transition metal series; Table S3: Gibbs free energy change (ΔG/eV) of the first protonation step in the CO2 reduction reaction (CRR) and H2 evolution reaction (HER) on the TM-HAB; Table S4: Gibbs free energy change for each protonation step of Zn-HAB electrocatalytic CO2 reduction; Table S5: Gibbs free energy change for each protonation step of Co-HAB electrocatalytic CO2 reduction; Table S6: Gibbs free energy change for each protonation step of Mn-HAB electrocatalytic CO2 reduction; Table S7: Gibbs free energy change for each protonation step of Fe-HAB electrocatalytic CO2 reduction; Table S8: Gibbs free energy change for ench protonation step of Ti-HAB electrocatalytic CO2 reduction; Table S9: Gibbs free energy change for each protonation step of Sc-HAB electrocatalytic CO2 reduction; Table S10: Gibbs free energy change for each protonation step of V-HAB electrocatalytic CO2 reduction; Table S11: Gibbs free energy change for ench protonation step of Cu-HAB electrocatalytic CO2 reduction; Table S12: Gibbs free energy change for ench protonation step of Cr-HAB electrocatalytic CO2 reduction.

Author Contributions

Conceptualization, X.Z. and Y.W.; methodology, C.X.; software, K.X. and Y.W.; validation, X.Z., C.X. and Y.W.; formal analysis, Z.T.; investigation, Z.T. and Y.Y.; data curation, L.L. and C.X.; writing—original draft preparation, X.Z. and C.X.; writing—review and editing, Z.T. and Y.W.; visualization, X.Z., C.X. and Y.W.; project administration, C.X., L.L. and Y.W.; funding acquisition, C.X., L.L. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant numbers 12064019 and 51661013) and the Science and Technology Program of the Education Office of Jiangxi Province (grant numbers GJJ212326 and GJJ202321).

Data Availability Statement

The data presented in this study are available in Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. He, M.; Sun, Y.; Han, B. Green carbon science: Scientific basis for the integration of efficient carbon resource processing, utilization and recycling. Angew. Chem. Int. Ed. 2013, 52, 9620–9633. [Google Scholar] [CrossRef] [PubMed]
  2. D’Alessandro, D.M.; Smit, B.; Long, J.R. Carbon dioxide capture: Prospects for new materials. Angew. Chem. Int. Ed. 2010, 49, 6058–6082. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Pires, J.C.M.; Martins, F.G.; Alvim-Ferraz, M.C.M.; Simões, M. Recent developments on carbon capture and storage: An overview. Chem. Eng. Res. Des. 2011, 88, 1446–1460. [Google Scholar] [CrossRef]
  4. Whipple, D.; Kenis, P. Prospects of CO2 utilization via direct heterogeneous electrochemical reduction. J. Phys. Chem. Lett. 2010, 1, 3451–3458. [Google Scholar] [CrossRef]
  5. Kim, S.K.; Zhang, Y.J.; Bergstrom, H.; Michalsky, R.; Peterson, A. Understanding the Low-Overpotential Production of CH4 from CO2 on Mo2C Catalysts. ACS Catal. 2016, 6, 2003–2013. [Google Scholar] [CrossRef]
  6. Li, N.; Chen, X.; Ong, W.J.; MacFarlane, D.R.; Zhao, X.; Cheetham, A.K.; Sun, C. Understanding of Electrochemical Mechanisms for CO2 Capture and Conversion into Hydrocarbon Fuels in Transition-Metal Carbides (MXenes). ACS Nano 2017, 11, 10825–10833. [Google Scholar] [CrossRef]
  7. Liu, C.; Yang, B.; Tyo, E.; Seifert, S.; DeBartolo, J.; von Issendorff, B.; Curtiss, L.A. Carbon Dioxide Conversion to Methanol over Size-Selected Cu4 Clusters at Low Pressures. J. Am. Chem. Soc. 2015, 137, 8676–8679. [Google Scholar] [CrossRef]
  8. Zaman, S.F.; Pasupulety, N.; Al-Zahrani, A.A.; Daous, M.A.; Al-Shahrani, S.S.; Inokawa, H.; Driss, H. Efficient homogeneous catalysis in the reduction of CO2 to CO. J. Am. Chem. Soc. 2005, 127, 17196–17197. [Google Scholar] [CrossRef]
  9. Bai, S.; Jiang, J.; Zhang, Q.; Xiong, Y. Steering charge kinetics in photocatalysis: Intersection of materials syntheses, characterization techniques and theoretical simulations. Chem. Soc. Rev. 2015, 44, 2893–2939. [Google Scholar] [CrossRef]
  10. Costentin, C.; Robert, M.; Savéant, J.M. Catalysis of the electrochemical reduction of carbon dioxide. Chem. Soc. Rev. 2013, 42, 2423–2436. [Google Scholar] [CrossRef]
  11. Xu, S.W.; Lu, Y.; Li, J.; Jiang, Z.Y.; Wu, H. Efficient Conversion of CO2 to Methanol Catalyzed by Three Dehydrogenases Co-encapsulated in an Alginate–Silica (ALG–SiO2) Hybrid Gel. Ind. Eng. Chem. Res. 2006, 45, 4567–4573. [Google Scholar] [CrossRef]
  12. Lu, Q.; Rosen, J.; Zhou, Y.; Hutchings, G.S.; Kimmel, Y.C.; Chen, J.G.; Jiao, F. A selective and efficient electrocatalyst for carbon dioxide reduction. Nat. Commun. 2014, 5, 3242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sun, Q.; Li, Z.; Searles, D.J.; Chen, Y.; Lu, G.M.; Du, A. Charge-Controlled Switchable CO2 Capture on Boron Nitride Nanomaterials. J. Am. Chem. Soc. 2013, 135, 8246–8253. [Google Scholar] [CrossRef]
  14. Qin, G.Q.; Du, A.J.; Sun, Q. Charge- and Electric-Field-Controlled Switchable Carbon Dioxide Capture and Gas Separation on a C2N Monolayer. Energy Technol. 2018, 6, 205–212. [Google Scholar] [CrossRef]
  15. Sun, Q.; Wang, M.; Li, Z.; Du, A.; Searles, D.J. Carbon Dioxide Capture and Gas Separation on B80 Fullerene. J. Phys. Chem. C 2014, 118, 2170–2177. [Google Scholar] [CrossRef]
  16. Sun, Q.; Wang, M.; Li, Z.; Du, A.; Searles, D.J. A computational study of carbon dioxide adsorption on solid boron. Phys. Chem. Chem. Phys. 2014, 16, 12695–12702. [Google Scholar] [CrossRef] [Green Version]
  17. Sun, Q.; Wang, M.; Li, Z.; Ma, Y.; Du, A. CO2 capture and gas separation on boron carbon nanotubes. Chem. Phys. Lett. 2013, 575, 59–66. [Google Scholar] [CrossRef] [Green Version]
  18. Qin, G.; Cui, Q.; Wang, W.; Li, P.; Du, A.; Sun, Q. First-Principles Study of Electrocatalytically Reversible CO2 Capture on Graphene-like C3N. ChemPhysChem 2018, 19, 2788–2795. [Google Scholar] [CrossRef]
  19. Peterson, A.A.; Abild-Pedersen, F.; Studt, F.; Rossmeisl, J.; Nørskov, J.K. How copper catalyzes the electroreduction of carbon dioxide into hydrocarbon fuels. Energy Environ. Sci. 2010, 3, 1311–1315. [Google Scholar] [CrossRef]
  20. Gao, D.; Zhang, Y.; Zhou, Z.; Cai, F.; Zhao, X.; Huang, W.; Li, Y.; Zhu, J.; Liu, P.; Yang, F.; et al. Enhancing CO2 Electroreduction with the Metal–Oxide Interface. J. Am. Chem. Soc. 2017, 139, 5652–5655. [Google Scholar] [CrossRef]
  21. Back, S.; Lim, J.; Kim, N.Y.; Kim, Y.H.; Jung, Y. Single-atom catalysts for CO2 electroreduction with significant activity and selectivity improvements. Chem. Sci. 2017, 8, 1090–1096. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Liu, X.; Yang, H.; He, J.; Liu, H.; Song, L.; Li, L.; Luo, J. Highly Active, Durable Ultrathin MoTe2 Layers for the Electroreduction of CO2 to CH4. Small 2018, 14, 1704049. [Google Scholar] [CrossRef] [PubMed]
  23. Li, P.Z.; Wang, X.J.; Liu, J.; Phang, H.S.; Li, Y.; Zhao, Y. Highly Effective Carbon Fixation via Catalytic Conversion of CO2 by an Acylamide-Containing Metal–Organic Framework. Chem. Mater. 2017, 29, 9256–9261. [Google Scholar] [CrossRef]
  24. Zhang, W.; Zheng, W. Single Atom Excels as the Smallest Functional Material. Adv. Funct. Mater. 2016, 26, 2988–2993. [Google Scholar] [CrossRef]
  25. Ling, C.; Li, Q.; Du, A.; Wang, J. Computation-Aided Design of Single-Atom Catalysts for One-Pot CO2 Capture, Activation, and Conversion. ACS Appl. Mater. Interfaces 2018, 10, 36866–36872. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, Z.; Zhao, J.; Cai, Q. CO2 electroreduction performance of a single transition metal atom supported on porphyrin-like graphene: A computational study. Phys. Chem. Chem. Phys. 2017, 19, 23113–23121. [Google Scholar] [CrossRef]
  27. Qiao, B.; Wang, A.; Yang, X.; Allard, L.F.; Jiang, Z.; Cui, Y.; Zhang, T. Singleatom catalysis of CO oxidation using Pt1/FeOx. Nat. Chem. 2011, 3, 634–641. [Google Scholar] [CrossRef]
  28. Lin, J.; Wang, A.; Qiao, B.; Liu, X.; Yang, X.; Wang, X.; Liang, J.; Li, J.; Liu, J.; Zhang, T. Remarkable Performance of Ir1/FeOx Single-Atom Catalyst in Water Gas Shift Reaction. J. Am. Chem. Soc. 2013, 135, 15314–15317. [Google Scholar] [CrossRef]
  29. Yang, X.F.; Wang, A.; Qiao, B.; Li, J.; Liu, J.; Zhang, T. Single-Atom Catalysts: A New Frontier in Heterogeneous Catalysis. Accounts Chem. Res. 2013, 46, 1740–1748. [Google Scholar] [CrossRef]
  30. Liu, N.; Huang, W.; Zhang, X.; Tang, L.; Wang, L.; Wang, Y.; Wu, M. Ultrathin graphene oxide encapsulated in uniform MIL-88A(Fe) for enhanced visible light-driven photodegradation of RhB. Appl. Catal. B Environ. 2018, 221, 119–128. [Google Scholar] [CrossRef]
  31. Zhang, X.; Yang, Y.; Huang, W.; Yang, Y.; Wang, Y.; He, C.; Tang, L. g-C3N4/UiO-66 nanohybrids with enhanced photocatalytic activities for the oxidation of dye under visible light irradiation. Mater. Res. Bull. 2018, 99, 349–358. [Google Scholar] [CrossRef]
  32. Huang, W.; Liu, N.; Zhang, X.; Wu, M.; Tang, L. Metal organic framework g-C3N4/MIL-53(Fe) heterojunctions with enhanced photocatalytic activity for Cr(VI) reduction under visible light. Appl. Surf. Sci. 2017, 425, 107–116. [Google Scholar] [CrossRef]
  33. Li, J.R.; Kuppler, R.J.; Zhou, H.C. Selective gas adsorption and separation in metal–organic frameworks. Chem. Soc. Rev. 2009, 38, 1477–1504. [Google Scholar] [CrossRef]
  34. Zhang, X.; Yang, Y.; Song, L.; Chen, J.; Yang, Y.; Wang, Y. Enhanced adsorption performance of gaseous toluene on defective UiO-66 metal organic framework: Equilibrium and kinetic studies. J. Hazard. Mater. 2019, 365, 597–605. [Google Scholar] [CrossRef] [PubMed]
  35. Zhang, X.; Yang, Y.; Lv, X.; Wang, Y.; Liu, N.; Chen, D.; Cui, L. Adsorption/desorption kinetics and breakthrough of gaseous toluene for modified microporous-mesoporous UiO-66 metal organic framework. J. Hazard. Mater. 2019, 366, 140–150. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, X.; Lv, X.; Shi, X.; Yang, Y.; Yang, Y. Enhanced hydrophobic UiO-66 (University of Oslo 66) metal-organic framework with high capacity and selectivity for toluene capture from high humid air. J. Colloid Interface Sci. 2019, 539, 152–160. [Google Scholar] [CrossRef] [PubMed]
  37. Liu, N.; Huang, W.; Tang, M.; Yin, C.; Gao, B.; Li, Z.; Tang, L.; Lei, J.; Cui, L.; Zhang, X. In-situ fabrication of needle-shaped MIL-53(Fe) with 1T-MoS2 and study on its enhanced photocatalytic mechanism of ibuprofen. Chem. Eng. J. 2019, 359, 254–264. [Google Scholar] [CrossRef]
  38. Cui, Q.; Qin, G.; Wang, W.; Geethalakshmi, K.R.; Du, A.; Sun, Q. Mo-based 2D MOF as a highly efficient electrocatalyst for reduction of N2 to NH3: A density functional theory study. J. Mater. Chem. A 2019, 7, 14510–14518. [Google Scholar] [CrossRef]
  39. Zhang, X.; Hou, F.; Yang, Y.; Wang, Y.; Liu, N.; Chen, D.; Yang, Y. A facile synthesis for cauliflower like CeO2 catalysts from Ce-BTC precursor and their catalytic performance for CO oxidation. Appl. Surf. Sci. 2017, 423, 771–779. [Google Scholar] [CrossRef]
  40. Zhang, X.; Zhang, X.; Song, L.; Hou, F.; Yang, Y.; Wang, Y.; Liu, N. Enhanced catalytic performance for CO oxidation and preferential CO oxidation over CuO/CeO2 catalysts synthesized from metal organic framework: Effects of preparation methods. Int. J. Hydrog. Energy 2018, 43, 18279–18288. [Google Scholar] [CrossRef]
  41. Zhang, X.; Li, H.; Hou, F.; Yang, Y.; Dong, H.; Liu, N.; Wang, Y.; Cui, L. Synthesis of highly efficient Mn2O3 catalysts for CO oxidation derived from Mn-MIL-100. Appl. Surf. Sci. 2017, 411, 27–33. [Google Scholar] [CrossRef]
  42. Yang, Y.; Dong, H.; Wang, Y.; He, C.; Wang, Y.; Zhang, X. Synthesis of octahedral like Cu-BTC derivatives derived from MOF calcined under different atmosphere for application in CO oxidation. J. Solid State Chem. 2018, 258, 582–587. [Google Scholar] [CrossRef]
  43. Kattel, S.; Wang, G. Reaction Pathway for Oxygen Reduction on FeN4 Embedded Graphene. J. Phys. Chem. Lett. 2014, 5, 452–456. [Google Scholar] [CrossRef] [PubMed]
  44. Li, X.F.; Li, Q.K.; Cheng, J.; Liu, L.; Yan, Q.; Wu, Y.; Zhang, X.H.; Wang, Z.Y.; Qiu, Q.; Luo, Y. Conversion of Dinitrogen to Ammonia by FeN3-Embedded Graphene. J. Am. Chem. Soc. 2016, 138, 8706–8709. [Google Scholar] [CrossRef]
  45. Mahmood, J.; Li, F.; Kim, C.; Choi, H.J.; Gwon, O.; Jung, S.M.; Seo, J.M.; Cho, S.J.; Ju, Y.W.; Jeong, H.Y.; et al. Fe@C2N: A highly-efficient indirect-contact oxygen reduction catalyst. Nano Energy 2018, 44, 304–310. [Google Scholar] [CrossRef]
  46. Tian, Y.; Wang, Y.; Yan, L.; Zhao, J.; Su, Z. Electrochemical reduction of carbon dioxide on the two–dimensional M3(Hexaiminotriphenylene)2 sheet: A computational study. Appl. Surf. Sci. 2019, 467, 98–103. [Google Scholar] [CrossRef]
  47. Lahiri, N.; Lotfizadeh, N.; Tsuchikawa, R.; Deshpande, V.V.; Louie, J. Hexaaminobenzene as a building block for a Family of 2D Coordination Polymers. J. Am. Chem. Soc. 2017, 139, 19–22. [Google Scholar] [CrossRef]
  48. Park, J.; Chen, Z.; Flores, R.A.; Wallnerstrom, G.; Kulkarni, A.; Nørskov, J.K.; Bao, Z. Two-Dimensional Conductive Ni-HAB as a Catalyst for the Electrochemical Oxygen Reduction Reaction. ACS Appl. Mater. Interfaces 2020, 12, 39074–39081. [Google Scholar] [CrossRef]
  49. Delley, B. From molecules to solids with the DMol3 approach. J. Chem. Phys. 2000, 113, 7756–7764. [Google Scholar] [CrossRef]
  50. Perdew, J.P.; Ernzerhof, M. Rationale for mixing exact exchange with density functional approximations. J. Chem. Phys. 1996, 105, 9982–9985. [Google Scholar] [CrossRef]
  51. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef] [PubMed]
  52. Sun, Q.; Sun, C.; Du, A.; Dou, S.; Li, Z. In-plane graphene/boron-nitride heterostructures as an efficient metal-free electrocatalyst for the oxygen reduction reaction. Nanoscale 2016, 8, 14084–14091. [Google Scholar] [CrossRef] [PubMed]
  53. Qin, G.; Du, A.; Sun, Q. A theoretical insight into a feasible strategy for the fabrication of borophane. Phys. Chem. Chem. Phys. 2018, 20, 16216–16221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Qin, G.; Cui, Q.; Yun, B.; Sun, L.; Du, A.; Sun, Q. High capacity and reversible hydrogen storage on two dimensional C2N monolayer membrane. Int. J. Hydrogen Energy 2018, 43, 9895–9901. [Google Scholar] [CrossRef]
  55. Delley, B. Hardness conserving semilocal pseudopotentials. Phys. Rev. B 2002, 66, 155125. [Google Scholar] [CrossRef]
  56. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  57. Klamt, A.; Schüürmann, G. COSMO: A new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient. J. Chem. Soc. Perkin Trans. 1993, 2, 799–805. [Google Scholar] [CrossRef]
  58. Back, S.; Min, S.Y.; Jung, Y. Active Sites of Au and Ag Nanoparticle Catalysts for CO2 Electroreduction to CO. Am. Chem. Soc. 2015, 5, 5089–5096. [Google Scholar] [CrossRef] [Green Version]
  59. Rossmeisl, J.; Logadottir, A.; Nørskov, J.K. Electrolysis of water on (oxidized) metal surfaces. Chem. Phys. 2005, 319, 178–184. [Google Scholar] [CrossRef]
  60. Nørskov, J.K.; Rossmeisl, J.; Logadottir, A.; Lindqvist, L.; Kitchin, J.R.; Bligaard, T.; Jónsson, H. Origin of the Overpotential for Oxygen Reduction at a Fuel-Cell Cathode. J. Phys. Chem. B 2004, 108, 17886–17892. [Google Scholar] [CrossRef]
  61. Liu, J.H.; Yang, L.M.; Ganz, E. Electrochemical reduction of CO2 by single atom catalyst TM–TCNQ monolayers. J. Mater. Chem. A 2019, 7, 3805–3814. [Google Scholar] [CrossRef]
  62. Santos-Lorenzo, J.; San José-Velado, R.; Albo, J.; Beobide, G.; Castaño, P.; Castillo, O.; Luque, A.; Pérez-Yáñez, S. A straightforward route to obtain zirconium based metal-organic gels. Microporous Mesoporous Mater. 2019, 284, 128–132. [Google Scholar] [CrossRef]
  63. Shi, C.; Chan, K.; Yoo, J.S.; Nørskov, J.K. Barriers of Electrochemical CO2 Reduction on Transition Metals. Org. Process. Res. Dev. 2016, 20, 1424–1430. [Google Scholar] [CrossRef]
  64. Li, M.; Yan, C.; Ramachandran, R.; Lan, Y.; Dai, H.; Shan, H.; Meng, X.; Cui, D.; Wang, F.; Xu, Z.X. Non-peripheral octamethyl-substituted cobalt phthalocyanine nanorods supported on N-doped reduced graphene oxide achieve efficient electrocatalytic CO2 reduction to CO. Chem. Eng. J. 2022, 430, 133050. [Google Scholar] [CrossRef]
  65. Merino-Garcia, I.; Albo, J.; Krzywda, P.; Mul, G.; Irabien, A. Bimetallic Cu-based hollow fibre electrodes for CO2 electroreduction. Catal. Today 2020, 346, 34–39. [Google Scholar] [CrossRef]
Figure 1. The top (up) and side (down) views of the structures in a 2 × 2 supercell for ten TM-HAB monolayers.
Figure 1. The top (up) and side (down) views of the structures in a 2 × 2 supercell for ten TM-HAB monolayers.
Nanomaterials 12 04005 g001
Figure 2. Stability of different TM-HAB monolayers structurally. TM is the transition metal atom, Ec is the cohesive energy of the TM bulk, Eb is the binding energy between TM and HAB in TM-HAB, and Ef is the formation energy of TM-HAB.
Figure 2. Stability of different TM-HAB monolayers structurally. TM is the transition metal atom, Ec is the cohesive energy of the TM bulk, Eb is the binding energy between TM and HAB in TM-HAB, and Ef is the formation energy of TM-HAB.
Nanomaterials 12 04005 g002
Figure 3. The Gibbs free energy change of the first protonation step in the CO2 reduction reaction (CO2RR) and H2 evolution reaction (HER). CO2RR selective catalysts are those located below the dotted line.
Figure 3. The Gibbs free energy change of the first protonation step in the CO2 reduction reaction (CO2RR) and H2 evolution reaction (HER). CO2RR selective catalysts are those located below the dotted line.
Nanomaterials 12 04005 g003
Figure 4. Flow chart of electrocatalytic CO2 reduction to C1 product scheme; red is the final product [61].
Figure 4. Flow chart of electrocatalytic CO2 reduction to C1 product scheme; red is the final product [61].
Nanomaterials 12 04005 g004
Figure 5. At zero potential, Gibbs free energy profiles for the CO2RR along the most favorable routes for Zn-HAB. A CO2 molecule in the gas phase with a clean catalytic surface is assigned a free energy of zero.
Figure 5. At zero potential, Gibbs free energy profiles for the CO2RR along the most favorable routes for Zn-HAB. A CO2 molecule in the gas phase with a clean catalytic surface is assigned a free energy of zero.
Nanomaterials 12 04005 g005
Figure 6. Gibbs free energy profiles for the CO2RR along the most favorable pathways for Co–HAB at zero potential. The free energy zero is set to a CO2 molecule in the gas phase with a clean catalyst surface.
Figure 6. Gibbs free energy profiles for the CO2RR along the most favorable pathways for Co–HAB at zero potential. The free energy zero is set to a CO2 molecule in the gas phase with a clean catalyst surface.
Nanomaterials 12 04005 g006
Figure 7. Gibbs free energy curves for (a) Mn-HAB and (b) Fe-HAB, at zero potential, along the most favorable path of the CO2RR. The free energy zero point was set as the carbon dioxide molecule in the gas phase with a clean catalyst surface.
Figure 7. Gibbs free energy curves for (a) Mn-HAB and (b) Fe-HAB, at zero potential, along the most favorable path of the CO2RR. The free energy zero point was set as the carbon dioxide molecule in the gas phase with a clean catalyst surface.
Nanomaterials 12 04005 g007
Figure 8. Gibbs free energy profiles for the CRR along the most favorable pathways for (a) Sc–HAB, (b) V–HAB, (c) Cu–HAB, and (d) Cr–HAB at zero potential. The free energy zero is set to a CO2 molecule in the gas phase with a clean catalyst surface.
Figure 8. Gibbs free energy profiles for the CRR along the most favorable pathways for (a) Sc–HAB, (b) V–HAB, (c) Cu–HAB, and (d) Cr–HAB at zero potential. The free energy zero is set to a CO2 molecule in the gas phase with a clean catalyst surface.
Nanomaterials 12 04005 g008
Figure 9. Gibbs free energy curve of CO2RR along the most favored path of Ti-HAB at zero potential. The free energy zero has been set as the free energy of CO2 molecules in the gas phase with a clean catalyst surface.
Figure 9. Gibbs free energy curve of CO2RR along the most favored path of Ti-HAB at zero potential. The free energy zero has been set as the free energy of CO2 molecules in the gas phase with a clean catalyst surface.
Nanomaterials 12 04005 g009
Figure 10. The projected partial density of states of *OCH2 adsorbed on Sc, *OH adsorbed on Ti and Cr, O*CHOH adsorbed on V and Zn, *CO adsorbed on Mn and Fe, *CHO adsorbed on Co, and *OHCH3 adsorbed on Cu. The dotted lines denote the Fermi level. The red, green, and blue lines represent the 3d orbital of the metal atoms, the 2p orbital of the oxygen atoms, and the 2p orbital of the carbon atoms, respectively.
Figure 10. The projected partial density of states of *OCH2 adsorbed on Sc, *OH adsorbed on Ti and Cr, O*CHOH adsorbed on V and Zn, *CO adsorbed on Mn and Fe, *CHO adsorbed on Co, and *OHCH3 adsorbed on Cu. The dotted lines denote the Fermi level. The red, green, and blue lines represent the 3d orbital of the metal atoms, the 2p orbital of the oxygen atoms, and the 2p orbital of the carbon atoms, respectively.
Nanomaterials 12 04005 g010
Table 1. Hirshfeld charge on metal atoms (QTM) and on nitrogen atoms (QN) for the ten TM-HAB model systems. The Hirshfeld spin of the metal atoms and the M-N bond length (RM−N) are shown.
Table 1. Hirshfeld charge on metal atoms (QTM) and on nitrogen atoms (QN) for the ten TM-HAB model systems. The Hirshfeld spin of the metal atoms and the M-N bond length (RM−N) are shown.
TM-HABQTMSpin-TMQN/eRTM-N/Å
Sc0.7080.000−0.2592.110
Ti0.5880.000−0.2361.962
V0.349−1.935−0.2021.959
Cr0.4460.000−0.2221.929
Mn0.321−3.225−0.2031.876
Fe0.1320.000−0.1621.924
Co0.0490.000−0.1531.837
Ni0.0520.000−0.1471.838
Cu0.3240.000−0.2081.945
Zn0.3940.000−0.2202.032
Table 2. Adsorption energy (Eads/eV) of different CO2 reduction products.
Table 2. Adsorption energy (Eads/eV) of different CO2 reduction products.
TM-HABCOHCOOHHCHOCH3OHCH4
Sc-HAB−1.263−1.432−0.998−1.361−0.100
Ti-HAB−2.767−2.013−2.461−2.334−0.848
V-HAB−1.910−0.804−1.296−1.258−0.125
Cr-HAB−2.209−1.558−0.161−0.838−0.109
Mn-HAB−1.617−0.081−0.073−0.274−0.243
Fe-HAB−0.165−0.073−0.070−0.123−0.071
Co-HAB−1.794−0.895−1.008−0.781−0.525
Ni-HAB−2.309−1.994−1.825−2.275−2.044
Cu-HAB−0.5470.0530.068−0.473−0.553
Zn-HAB−0.218−0.109−0.151−0.185−0.106
Table 3. Calculated potential-determining steps (PDS) for Sc-Zn (except Ni) for nine materials, the limiting potential (UL) for the CO2RR reaction, and the corresponding major products.
Table 3. Calculated potential-determining steps (PDS) for Sc-Zn (except Ni) for nine materials, the limiting potential (UL) for the CO2RR reaction, and the corresponding major products.
TM-HABPDSUL/VMain Products and Corresponding Overpotentials ( η /V)
Sc-HAB*OCH2 + H2O + H+ + e → *OCH3 + H2O−0.29CH4(0.46)
Ti-HAB*OH + CH4 + H2O + H+ + e → * + CH4 + 2H2O−1.14CH4(1.31)
V-HAB*OCHOH + H+ + e → *OCH + H2O−0.27CH4(0.44)
Cr-HAB*OH + CH4 + H2O + H+ + e → * + CH4 + 2H2O−0.27CH4(0.66)
Mn-HAB*CO + H2O + H+ + e → *CHO + 2H2O−0.27HCHO(0.2),CH3OH(0.29),CH4(0.44)
Fe-HAB*CO + H2O + H+ + e → *CHO + 2H2O−0.27HCHO(0.2),CH3OH(0.29),CH4(0.44)
Co-HAB*CHO + H2O + H+ + e → *OCH2 + 2H2O−0.53CH3OH(0.55),CH4(0.70)
Cu-HAB*CH3OH + H2O + H+ + e → *OH + CH4 + H2O−0.18CH4(0.35)
Zn-HAB*OCHOH → * + HCOOH0.24HCOOH(0.01)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zeng, X.; Tu, Z.; Yuan, Y.; Liao, L.; Xiao, C.; Wen, Y.; Xiong, K. Two-Dimensional Transition Metal-Hexaaminobenzene Monolayer Single-Atom Catalyst for Electrocatalytic Carbon Dioxide Reduction. Nanomaterials 2022, 12, 4005. https://doi.org/10.3390/nano12224005

AMA Style

Zeng X, Tu Z, Yuan Y, Liao L, Xiao C, Wen Y, Xiong K. Two-Dimensional Transition Metal-Hexaaminobenzene Monolayer Single-Atom Catalyst for Electrocatalytic Carbon Dioxide Reduction. Nanomaterials. 2022; 12(22):4005. https://doi.org/10.3390/nano12224005

Chicago/Turabian Style

Zeng, Xianshi, Zongxing Tu, Yanli Yuan, Luliang Liao, Chuncai Xiao, Yufeng Wen, and Kai Xiong. 2022. "Two-Dimensional Transition Metal-Hexaaminobenzene Monolayer Single-Atom Catalyst for Electrocatalytic Carbon Dioxide Reduction" Nanomaterials 12, no. 22: 4005. https://doi.org/10.3390/nano12224005

APA Style

Zeng, X., Tu, Z., Yuan, Y., Liao, L., Xiao, C., Wen, Y., & Xiong, K. (2022). Two-Dimensional Transition Metal-Hexaaminobenzene Monolayer Single-Atom Catalyst for Electrocatalytic Carbon Dioxide Reduction. Nanomaterials, 12(22), 4005. https://doi.org/10.3390/nano12224005

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

Article Metrics

Back to TopTop