Quanguo Jiang
Jianfeng Zhang1*
Zhimin Ao
Haiyan He- 1College of Mechanics and Materials, Hohai University, Nanjing, China
- 2Guangzhou Key Laboratory of Environmental Catalysis and Pollution Control, Institute of Environmental Health and Pollution Control, School of Environmental Science and Engineering, Guangdong University of Technology, Guangzhou, China
The CO oxidation mechanism on graphene with divacancy (DG) embedded with transition metal from Sc to Zn has been studied by using first principles calculations. The results indicate that O2 molecule is preferentially adsorbed on Sc, Ti, V, Cr, Mn, and Fe-DG, which can avoid the CO poisoning problem that many catalysts facing and is beneficial to the CO oxidation progress. Further study indicates that Mn-DG shows the best catalytic properties for CO oxidation with consideration of both Langmuir-Hinshelwood (LH) and Eley-Rideal (ER) oxidation mechanisms. Along the ER mechanism, the reaction energy barrier for the first step (CO free + O2 pre-adsorbed → OOCO) is 0.96 eV. Along the LH mechanism, the energy barrier for the rate limiting step (CO adsorbed + O2 adsorbed → OOCO) is only 0.41 eV, indicating that the CO oxidation on Mn-DG will occur along LH mechanism. The Hirshfeld charge distributions of O2 and CO molecules is tuned by the embedded Mn atom, and the charge transfer from the embedded Mn atom to the adsorbed molecules plays an important role for the CO oxidation. The result shows that the Mn-embedded divacancy graphene is a noble-metal free and efficient catalyst for CO oxidation at low temperature.
Introduction
Carbon monoxide (CO) is colorless, tasteless and toxic in air, while oxidation of CO is an efficient way to eliminate of the air pollutant (Xie et al., 2009). CO oxidation has important applications in atmosphere purification for hydrogen gas fuel in fuel cells as well (Qiao et al., 2015; Saavedra et al., 2016). Noble metals(Bleakley and Hu, 1999; Zhang and Hu, 2001; Liu et al., 2002; Gong et al., 2004; Dupont et al., 2006; Zhang et al., 2006; Liu, 2007) are common catalysts for the CO oxidation, where the rate limiting energy barriers are 0.46 eV (Liu et al., 2002) for Au(221), 0.91 eV (Gong et al., 2004) and 0.93 eV (Zhang and Hu, 2001) for Pd(111), 0.79 eV (Gong et al., 2004) and 0.82 eV (Dupont et al., 2006) for Pt(111), 1.17 eV (Gong et al., 2004) and 1.01 eV (Zhang et al., 2006) for Rh(111), 1.00 eV (Liu, 2007) for Rh(100). Due to the high cost and high reaction temperature of these noble metals, it is desirable to develop noble-metal-free catalysts for CO oxidation at low temperature. Noble metal clusters on supports are further studied to decrease the reaction barriers for CO oxidation (Tang et al., 2015; Ma et al., 2016; Wang et al., 2016; Ali et al., 2017; Chen et al., 2017). Furthermore, single atom catalyst decorated on appropriate matrix is attracted a lot of interests due to the excellent catalytic performance (Dvorák et al., 2016; Jones et al., 2016). The outstanding physical properties (Novoselov et al., 2005; Balandin et al., 2008; Lee et al., 2008) and the large surface-to-volume ratio make graphene (Novoselov et al., 2004) to be a promising substrate to realize high performance single atom catalysis. However, the inert nature of pristine graphene usually causes clustering problems for the adsorbed metal atoms (Liu and Huang, 2017; Liu et al., 2017). The interactions between the pristine graphene and the adsorbed atoms can be enhanced by introducing different carbon vacancies. Experimentally, single carbon vacancy (SV) and double carbon vacancies (DV) are common point defects on graphene. Vacancies can be introduced into graphene by exposing it to the focussed electron beam, and the vacancy defects in graphene can be tailored by controlling the exposing time, where the DV is more often observed than SV due to the high energy of the SV under electron beam irradiation (Robertson et al., 2012). In addition, the concentration of DV changes from 0.1 to 0.5 nm2 by controlling the total electron beam dose (Robertson et al., 2012). In addition, it is reported that the metal atoms resident on the SV and DV positions in graphene surfaces are stable in comparison to that on graphene edge (Robertson et al., 2013). Hence, it is important to evaluate the catalytic performance of a single atom supported on graphene with both SV and DV defects.
Many works about the CO oxidation have been reported for single atom decorated at the single carbon vacancy on graphene. Theoretically, Au- (Lu et al., 2009), Fe- (Li et al., 2010), Cu- (Song et al., 2011), Pt- (Tang et al., 2012), Si- (Zhao et al., 2012), and Al-embedded (Jiang et al., 2014) SV graphenes show high activity for the CO oxidation. Although the SV graphene with decorating metal atoms has high catalytic activity, decorating metal atom on the controllable carbon vacancies in graphene to realize different single atom catalyst is still highly desirable. The decorated metal atom is three-bond coordinated on SV graphene, while it is four-bond coordinated on DV graphene, indicating that the chemical activity of the decorated metal atom on different carbon vacancies should be different due to the different chemical environment. Furthermore, as mentioned above divacancy is commonly present in graphene obtained through chemical synthesis. However, only the catalytic activity of Fe-decorated DV graphene has been studied (Tang et al., 2016; Liu et al., 2017). Therefore, further research about the catalytic performance of DV graphene decorated with metal atom for CO oxidation is needed.
For the CO oxidation, Langmuir-Hinshelwood (LH) and Eley-Rideal (ER) are two mainly mechanisms. Along the ER mechanism, the activated O2 molecule directly reacts with the free CO molecule, where the activation of O2 is the rate-limiting step (Lu et al., 2009). Alone the LH mechanism, the CO and O2 molecules are first coadsorbed, and then react to form OOCO intermediate, which is the rate-limiting step for the oxidation progress (Lu et al., 2009). In general, the reaction energy barriers are proportional to the adsorption energy of adsorbed molecules on supported catalyst (Gong et al., 2004), indicating that the adsorption energy of adsorbed CO and O2 molecules could be a benchmark for the catalytic performance of graphene for CO oxidation. In addition, a larger adsorption energy for O2 molecule than that of CO is desired during the CO oxidation progress, because the preferential adsorption of CO will block the active sites and prevent the continuous oxidation reaction (Jiang et al., 2014; Tang et al., 2016). Therefore, the adsorption energy for O2 and CO on transition metals doped graphene is first calculated and the alternative decorating atom is chosen based on this rule.
In this work, by using first principles calculations, we will systematically study the CO oxidation mechanism on DV graphene decorated with transition metals from Sc to Zn, which are non-noble metals and commonly used to decorate two dimensional materials. The reaction barriers for each step are analyzed and the corresponding reaction mechanisms are discussed through analyzing the electronic property of the graphene systems.
Calculation Methods
The density functional theory (DFT) calculations are carried out by using Dmol3 package (Delley, 2000). Exchange-correlation functions are taken as generalized gradient approximation (GGA) with Perdew-Burke-Ernzerhof (PBE) (Perdew et al., 1996). The selection of exchange–correlation functional has evidential effect on the result of adsorption energies, while has much smaller effect on the reaction energy barriers (Roldán et al., 2009). DFT semicore pseudopotentials (DSPPs) core treatment is implemented for relativistic effects. Double numerical plus polarization (DNP) is employed as the basis set. The convergence tolerance of energy of 10−5 Hartree is taken (1 Hartree = 27.21 eV), and the maximal allowed force and displacement are 0.002 Hartree/Å and 0.005 Å, respectively. Linear synchronous transit/quadratic synchronous transit (LST/QST) (Halgren and Lipscomb, 1977) and nudged elastic band (NEB) (Henkelman and Jonsson, 2000) tools in Dmol3 module are used to investigate the minimum energy pathway for CO oxidation on graphene, which have been well validated to determine the transition state. Three-dimensional periodic boundary conditions are taken in the simulation. The simulation cell consists of a 4 × 4 graphene supercell with a vacuum width of 20 Å above the graphene layer to minimize the interlayer interaction. The k-point is set to 5 × 5 × 1, and all atoms are allowed to relax according to previous reports (Jiang et al., 2014). After structure relaxations, the density of states (DOS) are calculated with a finer k-point grid of 15 × 15 × 1. The DFT+D method within the Grimme scheme (Grimme, 2006) is used in all calculations to consider the van der Waals forces. The electron orbits of the free and adsorbed molecules are calculated with CASTEP code (Segall et al., 2002), where the ultrasoft pseudopotentials, GGA-PBE functional, an energy cutoff of 340 eV and 5 × 5 × 1 k-point meshes are used. We have compared the total energy of the graphene system with different spin state, and choose the proper spin state for the graphene system with the smallest total energy.
The adsorption energy Ead of molecules on graphene is determined by,
where Emolecules/graphene, Egraphene, and Emolecule are total energies of the adsorbed graphene system, the isolate graphene and molecules respectively.
Results and Discussion
Adsorption of O2 and CO on TM-DG
Based on the literature results in the introduction section, it is expected that transition-metal-embedded divacancy graphene (TM-DG) could also exhibit excellent catalytic behaviors for CO oxidation, similar to the cases of graphene with single vacancy systems. Herein, we present systematic DFT calculations on the adsorption energy of O2 and CO molecules on the divacancy graphene embedded with transition metal (from Sc to Zn). To comprehensively understanding the adsorption behaviors of O2 and CO molecules on transition metal decorated divacancy graphene, the adsorption energies of O2 and CO, as well as the co-adsorption energy of O2 and CO molecules on TM-DG are shown in Figure 1, where we can see that the adsorption energy of O2 molecule is larger than that of CO molecule on TM-DG (from Sc to Fe), while the adsorption energy of CO molecule is larger than that of O2 molecule on TM-DG (from Co to Zn). This result indicates that the O2 molecule has priority during the adsorption progress and avoids the CO poisoning at the active sites on TM-DG (from Sc to Fe) during the CO oxidation. In addition, the co-adsorption energy of O2 and CO molecules on TM-DG (from Sc to Fe) is also shown in Figure 1, where a local adsorption energy minimum on Mn-DG is found, which indicates that Mn-DG can facilitate the CO oxidation better, due to the fact that the energy barrier is proportional to the adsorption energy of molecules (Gong et al., 2004). Therefore, CO oxidation on Mn-DG is mainly studied in the following, and the CO oxidation on other TM-DG (Sc, Ti, V, Cr and Fe) is also studied for comparison purpose.
Figure 1. The adsorption energies of O2 and CO as well as the co-adsorption energy of O2 and CO on TM-DG, where the light and dark areas indicate the preferential adsorption of O2 and CO molecules, respectively.
Electronic Properties and Stability of Mn-DG
Before further study the catalytic properties of Mn-DG, the electronic properties and structure stability of Mn-DG are first studied in Figure 2. After embedding Mn atom into the divancancy graphene through substituting two carbon atoms, the reconstructed structure of graphene is shown in Figure 2A, where four chemical Mn-C bonds are formed in graphene with bond length lMn−C = 1.99 Å. Mn atom is out of the plane of graphene with a distance of 0.68 Å due to the larger atomic radius of 1.79 Å compared with that of C 0.91 Å. The adsorption energy of the Mn atom on the divacancy graphene is −6.81 eV. The Hirshfeld charge distributions near the dopant are also given in Figure 2A, where the electron-deficiency position is formed for the Mn atom with losing electrons of 0.284 e, which promotes the adsorption of O2 and CO molecules. The differential charge density along C-Mn-C bonds for Mn-DG is further studied and is shown in Figure 2B, where the blue and red isosurfaces correspond to the increase in the number of electrons and the depletion zone, respectively. It shows that electrons accumulate near the doped Mn atom, indicating the high chemical active area.
Figure 2. (A) The atomic structure of divacancy graphene embedded with Mn atom (Mn-DG), where the gray and pink atoms are C and Mn atoms, respectively, in this and following figures. The atomic charge obtained by Hirshfeld analysis near the Mn dopant is also given. (B) The differential charge density along C-Mn-C bonds for Mn-DG is shown, where the blue and red isosurfaces correspond to the increase in the number of electrons and the depletion zone, respectively. (C) The diffusion pathway of the Mn atom on Mn-DG, where IS, TS, and FS represent initial, transition and final structures, respectively, in this and following figures. (D) PDOS of Mn atom and C atom on the Mn-DG, where the vertical line indicates the Fermi level. Inset is the positive and negative spin density of Mn-DG shown in green and red, respectively.
The aggregation problems for the adsorbed metal atoms on substrate are significant for the catalytic performance, especially when the concentration of metal atom is high (Ao and Peeters, 2010a). To determine the possibility of aggregation for Mn atoms on divacancy graphene, the diffusion pathway of Mn atom to its neighboring positions is investigated based on DFT calculations (see Figure 2C), where the corresponding diffusion energy barrier for the decorated Mn atom is 5.91 eV. It is claimed that a surface reaction will occur when the reaction barrier is smaller than the critical value of Ecbar = 0.91 eV (Young, 2001), the decorated Mn atom on divacancy graphene is thus stable. In addition, the adsorption energy of Mn atom on divacancy graphene is −6.81 eV and it is much larger than the cohesive energy −2.92 eV/atom for Mn element (http://www.knowledgedoor.com/). Therefore, the Mn decorated divacancy graphene is quite stable without aggregation problems. Partial density of states (PDOS) are further analyzed to confirm the enhanced interactions between the Mn atom and graphene (Figure 2D), where the energy bands between the decorated Mn atom and carbon atoms overlap significantly. The electron-deficiency character is confirmed by the Fermi level crossing the valence band, and it is consistent with the Hirshfeld charge distributions. As shown in the insert of Figure 2D, the positive and negative spin density of Mn-DG is shown in green and red, respectively. Therefore, the magnetic moment of Mn-DG is 3.053 μB, which is mainly contributed by the Mn atom. The remaining unsaturated d orbital of Mn atom is reactive, which can adsorb small molecules and promote the subsequent reactions.
Adsorption of Molecules on Mn-DG
To investigate the oxidation of CO on Mn-embedded divacancy graphene, the adsorptions of O2 and CO on Mn-DG are studied carefully. Figure 3A shows the most stable configuration for O2 molecule adsorbed on Mn-DG (Ead = −1.96 eV), where the O-O bond is parallel to the graphene sheet, and ~0.18 e is transferred from Mn-DG to O2 molecule based on Hirshfeld charge analysis. To assess the stability of the adsorbed O2 molecule, the dissociative adsorption of O2 molecule on Mn-DG is then studied in Figure 3B. After NEB calculations, the dissociation reaction barrier for an O2 molecule on the Mn-DG is 1.18 eV > Ecbar = 0.91 eV (Young, 2001), which indicates that the adsorbed O2 molecule prefers to stay on Mn-DG in molecular state at room temperature. The PDOS of free O2 molecule is shown in Figure 3C, where the 2π* anti-bond orbital is half filled (Honkala and Laasonen, 2000). Then, the PDOS of Mn atom and the adsorbed O2 molecule is shown in Figure 3D, where all orbitals of the adsorbed O2 molecule are also labeled to understand the interaction between the adsorbed O2 and Mn-DG. About 0.18 e is transferred from Mn-DG to the adsorbed O2 molecule based on Hirshfeld method, which occupies the O2-2π* orbital above the Fermi level (see Figure 3C) and is confirmed by the new peaks for O2-2π* orbital below the Fermi level (see Figure 3D). This charge transfer elongates the O-O bond from 1.23 Å in free O2 to 1.40 Å in adsorbed O2 molecule. The activated O2 molecule with longer O-O bond will be beneficial for the subsequent CO oxidation. The O2-2π* orbitals and Mn atom is strong hybridized near the Fermi level (see Figure 3D), which mainly responses for the chemical adsorption of O2 molecule on Mn-DG.
Figure 3. (A) The most stable structure of O2 adsorbed on Mn-DG. (B) The reaction pathway for the dissociative adsorption of O2 molecule on Mn-DG. (C) PDOS and orbitals of the free O2 molecule, where the number of electrons for each orbital is also shown by arrows. (D) PDOS of the adsorbed O2 molecule and Mn atom for Mn-DG, inset is the charge density of 4s, 5s, 1p, and 2p* orbitals of the adsorbed O2 molecule. The vertical lines indicate the Fermi level.
The adsorption configuration of a CO molecule on Mn-DG is shown in Figure 4A, where the CO molecule is vertically adsorbed on the top of the decorated Mn atom. CO is chemically adsorbed on Mn-DG, which is confirmed by the chemical bond between Mn atom and the carbon atom of CO molecule. The adsorption energy of a CO molecule on Mn-DG is Ead = −1.70 eV, and CO molecule loses 0.006 e to the Mn-DG. The binding energy Eb of C-O bond is 11.57 eV, which is much larger than Eb = 6.36 eV for O-O bond based on DFT calculations. Thus CO should be more difficult to be dissociated, which is confirmed by the fact that the dissociative energy barrier for CO molecule on Mn-DG is 5.55 eV based on DFT calculation as shown in Figure 4B. The PDOS of free CO molecule is shown in Figure 4C, where all orbitals are labeled. Figure 4D shows the PDOS of the adsorbed CO molecule on Mn-embedded divacancy graphene, where the orbitals of adsorbed CO molecule are labeled and displayed. The 5σ peak of the CO molecule adsorbed on Mn-DG is significantly depressed than the free CO molecule due to the charge transfer. Although the 2π* anti-bond orbital far above Fermi level for free CO molecule is fully empty, Mn atom transfers some electrons to CO-2π* orbital due to the fact that CO-2π* is close to the Fermi level at absorbed state, which slightly elongates the C-O bond from 1.14 Å for free CO molecule to 1.16 Å for adsorbed CO molecule. In addition, CO molecules will act as a donor with the carbon atom near the graphene surface, dur to fact that the CO-5σ orbital locates on the carbon atom (Leenaerts et al., 2008). This agrees with the above Hirshfeld analysis result. The above discussions indicate that O2 and CO have strong interactions with Mn-DG (corresponding adsorption energies are −1.96 and −1.70 eV, respectively), but the adsorption of O2 is much stronger. The O2 molecule is activated on Mn-DG, which will facilitate the CO oxidation progress on graphene.
Figure 4. (A) The most stable structure of CO adsorbed on Mn-DG. (B) The reaction pathway for the dissociative adsorption of CO molecule on Mn-DG. (C) PDOS and orbitals of the free CO molecule, where the number of electrons for each orbital is also shown by arrows. (D) PDOS of the adsorbed CO molecule and Mn atom for the Mn-DG, inset is the charge density of 4σ, 1π, 5σ, and 2π* orbitals of the adsorbed CO molecule. The vertical lines indicate the Fermi level.
CO Oxidation on Mn-DG
Two reaction mechanisms have been established for the oxidation of CO molecule: Langmuir-Hinshelwood (LH) mechanism and Eley-Rideal (ER) mechanism (Molina and Hammer, 2005; An et al., 2008; Lu et al., 2009; Li et al., 2010; Song et al., 2011; Tang et al., 2012; Zhao et al., 2012; Jiang et al., 2014). For the ER mechanism, the O2 molecule is first adsorbed and activated by the Mn-DG, then a free CO molecule approaches the substrate to form an intermediate product. Along the LH mechanism, the O2 and CO molecules first co-adsorb on the Mn-DG, and then form an intermediate product. Since O2 has a larger adsorption energy (−1.96 eV) on Mn-DG than that of CO (−1.70 eV), the adsorption of O2 on Mn-DG has higher priority, thus the ER mechanism for the CO oxidation seems to be favorable. However, the lower co-adsorption energy for O2 and CO molecules (−2.06 eV) indicates that O2 and CO may co-adsorb on Mn atom as discussed in literatures (Lu et al., 2009; Song et al., 2011; Tang et al., 2012; Jiang et al., 2014). Therefore, both mechanisms for CO oxidation will be discussed in the following.
In order to search the preferred reaction path for CO oxidation on Mn-DG, the first reaction step along both LH and ER mechanisms is studied in Figure 5. If the LH mechanism is more favorable, the free CO molecule will co-adsorb on Mn atom with pre-adsorbed O2 molecule after overcoming a small reaction barrier. If the ER mechanism is more favorable, the free CO molecule will react with the O2 molecule to form an intermediate (OOCO). The structure of physisorbed CO molecule on the Mn-DG with pre-adsorbed O2 molecule was selected as the IS after studying all adsorption configurations as shown in Figure 5A. The reaction profile for the co-adsorption of CO and O2 molecules on Mn atom is shown in Figure 5A, where the CO and O2 are titled toward the graphene surface at TS state. After overcoming an energy barrier of 0.19 eV, the CO and O2 co-adsorb on Mn-DG with releasing energy of 0.02 eV, as shown in Figure 5A.
Figure 5. (A) The reaction pathway for the adsorption of CO molecule on Mn-embedded graphene with pre-adsorbed O2 molecule. (B) The primary reaction step for the oxidation of CO molecule on Mn-embedded graphene for ER mechanism.
Figure 5B shows the reaction profile for intermediate product along the ER mechanism. When the free CO approaches the Mn atom, the O-O bond is broken (see TS in Figure 5B). After overcoming a relative large energy barrier of 0.96 eV, the C atom from the CO binds with two O atoms from the dissociation of the O2 molecule to form a OOCO complex over the Mn atom (see FS in Figure 5B). This exothermic process releases an energy of 3.39 eV.
The effect of temperature on the energy barriers is considered. The free energy change (ΔG) between the reactant and transition state is considered as the temperature-dependent energy barrier E′bar and ΔG = ΔH- TΔS, where ΔH is the enthalpy change, ΔS is the entropy change, and T is the room temperature (298.15 K). In addition, ΔH = (ΔU + PΔV), ΔU = (ΔEtot + ΔEvib + ΔEtrans + ΔErot) and ΔS = ΔSvib + ΔStrans + ΔSrot, where ΔU is the internal energy change, ΔEtot is the change of total electronic energy, the vib, trans and rot indicate vibration, translation and rotation, respectively, which can be obtained through calculations of vibrational frequency. The temperature-dependent energy barriers at 298.15 K in Figure 5 are: E′bar = 0.33 eV in Figure 5A and E′bar = 1.15 eV in Figure 5B. This indicates that the free CO molecule will desorb from substrate with increasing temperature, and the reaction is more difficult to occur. The co-adsorption of CO and O2 on Mn-DG in Figure 5A will happen at 298.15 K due to the lower reaction barrier than the critical barrier of Ecbar = 0.91 eV (Young, 2001), thus the LH mechanism is mainly studied for the CO oxidation in the following.
For oxidation of CO on the Mn-DG along LH mechanism, several steps and also intermediate products (MS) for the oxidation procedure exist (Lu et al., 2009; Song et al., 2011; Tang et al., 2012; Jiang et al., 2014). For each step, e.g., from initial state to intermediate state in Figure 6, a transition state also exists. The configuration of co-adsorbed CO and O2 molecules on Mn-DG is taken as the reactant (IS in Figure 6A) based on the above discussions. At transition state, one O-Mn bond changes from 1.93 to 2.10 Å and a C-O bond between CO and O2 is formed. After overcoming an energy barrier of 0.41 eV, an OOCO intermediate (MS in Figure 6A) is formed and the elongated O-Mn bond is cleaved. Then the first CO2 molecule (MS2) is formed on Mn-DG after overcoming an energy barrier of 0.15 eV (see Figure 6A), where the CO2 is physically adsorbed on Mn-DG and its adsorption energy is −0.20 eV. The reaction for this step can release a heat of 3.81 eV, which can sufficiently overcome the adsorption energy of CO2, and the first produced CO2 molecule would desorb from the Mn-DG efficiently. The O-O bond for OOCO configuration in MS is broken at TS2, and a CO2 molecule with bond angle of 126.6° is formed. Figure 6B shows the formation of the second CO2 molecule on Mn-embedded divacancy graphene. The subsequent CO molecule will react with the remaining O atom to produce CO2 after surmounting an energy barrier of 0.17 eV (see Figure 6B), similar to the Au-embedded graphene (Lu et al., 2009). The reaction for this step releases a heat of 0.90 eV, which can also surmount the adsorption energy of CO2 (−0.27 eV) in Mn-DG, and the second CO2 molecule will desorb from Mn-DG efficiently. The reaction profile for the production of the two CO2 molecules is shown in Figure 6, where the rate limiting step is the formation of OOCO intermediate and the energy barrier of 0.41 eV is quite small. Therefore, the CO oxidation reaction along the LH mechanism releases 0.21 eV, 3.81 eV, and 0.90 eV for step 1, step 2, and step 3, respectively, which indicates the favorable thermodynamics for the CO oxidation on Mn-DG. In addition, the small energy barrier of 0.41 eV for the rate limiting step indicates the favorable kinetics for the CO oxidation on Mn-DG.
Figure 6. The reaction pathway of CO oxidation on Mn-embedded graphene for LH mechanism: production of the first CO2 molecule (A), and the second CO2 moelcule (B). The energy profile for the formation of CO2 molecule is also shown. The unit of E is eV, where E is the energy barrier. The Hirshfeld charge near the adsorbate is also shown.
The energy barrier of each step for the CO oxidation along LH mechanism on divacancy graphene decorated with transition metals (from Sc to Fe) are also studied as shown in Table 1, where the rate limiting energy barriers are 1.15, 0.95, 0.85, 1.17, 0.41, and 0.88 eV for Sc-, Ti-, V-, Cr-, Mn-, and Fe-DG, respectively. The reaction barrier is the smallest on Mn-DG, which confirms the fact that the reaction barriers of molecules are proportional to the adsorption energy on supported catalyst (Gong et al., 2004), and the co-adsorption energy of adsorbed CO and O2 molecules can be a benchmark for the catalytic performance of the graphene for CO oxidation. Therefore, CO can be oxidized easily on Mn-DG at low temperature, indicating that Mn-DG is an excellent candidate during the catalysts for CO oxidation.
Table 1. The reaction energy barriers of each step for the CO oxidation on divacancy graphene decorated with transition metal (from Sc to Fe) along the LH mechanism, where Ebar1, Ebar2, and Ebar3 are the energy barriers for step 1, step 2, and step 3 during CO oxidation, respectively.
The energy barriers at 298.15 K for each step of CO oxidation on Mn-DG in Figure 6 are: E′bar = 0.49 eV for step 1, E′bar = 0.10 eV for step 2, and E′bar = 0.25 eV for step 3, where the energy barriers for step 1 and step 3 slightly increase than those at 0 K. The reaction time for each step in Figure 6 is calculated by Arrhenius equation (Pauling, 1988)
where KB is the Boltzmann constant, ν is in order of 1012 Hz, and T = 298.15 K. τ1 = 1.9 × 10−4 s for step 1, τ2 = 4.9 × 10−11 s for step 2, and τ3 = 1.7 × 10−8 s for step 3, respectively. Therefore, the CO oxidation on Mn-DG along LH mechanism has fast kinetics.
The possibility for the reversing reaction of CO oxidation, i.e., for the dissociation of the first and the second CO2 molecule into CO molecule and O atom on Mn-DG are further considered as shown in Figure 6. The reversing energy barriers for the step 1, step 2, and step 3 based on DFT calculations are 0.62, 3.96, and 1.07 eV, respectively. After considering the effect of temperature, the reversing energy barriers at 298.15 K are: E′bar = 0.69 eV for step 1, E′bar = 4.05 eV for step 2, E′bar = 1.01 eV for step 3, while the corresponding reaction time is: τ1 = 0.44 s for step 1, τ2 = 2.3 × 1056 s for step 2, and τ3 = 1.1 × 105 s for step 3. The long reaction time indicates that the reversing reactions for the CO oxidation on Mn-DG are hardly to occur and the CO oxidation can be finished thoroughly.
To consider the effect of interactions between vacancies on the oxidation of CO, the graphene with more vacancies in the 4 × 4 supercell is studied as shown in Figure 7. The atomic structure of graphene with two divacancies embedded with Mn atoms is shown in Figure 7A, and the corresponding reaction pathway of CO oxidation for LH mechanism is shown in Figure 7B, where the rate limiting energy barrier is 0.65 eV for the formation of OOCO complex. The atomic structure of graphene with one divacancy and one single vacancy embedded with Mn atoms is shown in Figure 7C, and the corresponding reaction pathway of CO oxidation for LH mechanism is shown in Figure 7D, where the rate limiting energy barrier is 0.88 eV for the formation of OOCO complex. Therefore, the reaction barriers on graphene with two double vacancies or one single and one double vacancies in the supercell in Figure 7 are both large than that on graphene with one double vacancy in Figure 6. Considering that the distances between the nearest Mn dopants in Figures 7A,B are 4.295 Å and 5.062 Å, while that in the periodic supercell with one double vacancy in Figure 2A is 9.908 Å with Ebar = 0.41 eV, thus the relative larger distance between carbon vacancies is beneficial for the CO oxidation on Mn-DG.
Figure 7. The atomic structure of graphene with two divacancies embedded with Mn atoms (A), and the corresponding reaction pathway of CO oxidation along LH mechanism (B). The atomic structure of graphene with one divacancy and one single vacancy embedded with Mn atoms (C), and the corresponding reaction pathway of CO oxidation along LH mechanism (D).
Origin of the High Activity of Mn-DG
To further understand the high activity of the Mn-DG for CO oxidation, the Hirshfeld charge near the adsorbed molecules along the LH reaction path is shown in Figure 6, and the corresponding PDOS near the adsorbed molecules is shown in Figure 8. O2 obtains 0.115 e while CO loses 0.076 e for the initial state in Figure 6A. The empty component of the O2-2π* and CO-2π* orbitals are partially filled (see Figure 8A) due to the electron transformation during the adsorption, which causes the elongation of the O-O bond and C-O bond to 1.34 Å and 1.15 Å, respectively. At transition state, the charge of adsorbed O2 is −0.136 e while that of adsorbed CO is −0.006 e in Figure 6A. This indicates that the O2-2π* and CO-2π* orbitals are more filled as shown in Figure 8B, which is confirmed by the elongation of the O-O bond and C-O bond to 1.46 and 1.17 Å, respectively. After surmounting the energy barrier, more charge transfers from graphene to the adsorbed O2 and CO, where the charge of O2 is −0.154 e while that of CO is −0.058 e for MS in Figure 6A. The O2-2p* and CO-2π* orbitals are more filled as shown in Figure 8C. Therefore, the Mn-DG can tune the charge distributions of the adsorbed O2 and CO, and the charge transfer from Mn-DG to O2 and CO molecules plays an important role for the OOCO intermediate formation.
Figure 8. PDOS of adsorbed O2, CO, and Mn atom on Mn-DG for the IS (A), TS (B), and MS (C). PDOS of adsorbed O2, CO and Fe atom on Fe-DG for the TS state is also shown in (D). Red, blue and black curves represent PDOS of adsorbed O2, CO, and Mn (or Fe) atom, respectively. The orbitals of the adsorbed O2 and CO molecules are roughly labeled. Inset is the charge density of orbitals near the Fermi level. The vertical lines indicate the Fermi level.
Due to the fact that Fe is a commonly used dopant and the energy barrier for the first step on Fe-DG is 0.88 eV for CO oxidation, which is larger than that of 0.41 eV on Mn-DG, the PDOS of the TS configurations on Fe-DG is analyzed to further understood the mechanism for depressing the formation energy barrier of OOCO intermediate, as shown in Figure 8D. It can be seen that the overlapping area between O2-2π* and Mn-3d orbitals at Fermi level on Mn-DG (see Figure 8B) becomes much weaker compared with that on Fe-DG (see Figure 8D), indicating that the interaction of Mn-O bond is significant weakened on Mn-DG. It is reported that weaker interaction nearby the Fermi level generally causes smaller reaction barrier (Arellano et al., 2000; Ao and Peeters, 2010b; Jiang et al., 2016). In addition, the interaction from −16 to −5 eV is strengthened on Mn-DG (see Figures 8B,D). Therefore, the enhanced interactions in the low energy range could lead to smaller reaction barrier for the CO oxidation on Mn-DG.
Conclusion
The oxidation of CO molecule on transition metals decorated graphene with divacancy (DG) has been studied by using DFT calculations. We found that Mn-DG has the best performance for the CO oxidation, while the LH mechanism is preferred, where the rate limiting energy barrier is only 0.41 eV, indicating the efficient oxidation process. The charge transfer from Mn-DG to the O2-2p* and CO-2π* orbitals through the Mn atom along the LH mechanism plays a key role for depressing the energy barrier of the CO oxidation. The results indicate that Mn-DG can be a noble-metal-free and efficient catalyst for CO oxidation.
Author Contributions
QJ, JZ, and ZA designed the research and wrote the paper. QJ carried out the simulation. HJH, HYH, and YW entered the discussion. All authors commented on the manuscript.
Funding
We acknowledge supports by the Fundamental Research Funds for National Natural Science Foundation of China (Grant No. 21703052), the Central Universities (Grant Nos. 2017B12914 and 2015B01914), China Postdoctoral Science Foundation (2015M571652), Natural Science Foundation of Jiangsu Province (BK20161506) and National 973 Plan Project (2015CB057803). ZA acknowledges the financial supports from National Natural Science Foundation of China (Grant No. 21607029, 21777033), Science and Technology Program of Guangdong Province (2017B020216003), Science and Technology Program of Guangzhou City (201707010359), and 1000 Plan for Young Professionals Program of China, 100 Talents Program of Guangdong University of Technology.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Ali, S., Liu, T. F., Lian, Z., Li, B., and Su, D. S. (2017). The tunable effect of nitrogen and boron dopants on a single walled carbon nanotube support on the catalytic properties of a single gold atom catalyst: a first principles study of CO oxidation. J. Mater. Chem. A 5, 16653–16662. doi: 10.1039/C7TA03551A
An, W., Pei, Y., and Zeng, X. C. (2008). CO oxidation catalyzed by single-walled helical gold nanotube. Nano Lett. 8, 195–202. doi: 10.1021/nl072409t
Ao, Z. M., and Peeters, F. M. (2010a). High-capacity hydrogen storage in Al-adsorbed graphene. Phys. Rev. B 81:205406. doi: 10.1103/PhysRevB.81.205406
Ao, Z. M., and Peeters, F. M. (2010b). Electric field: A catalyst for hydrogenation of graphene. Appl. Phys. Lett. 96:253106. doi: 10.1063/1.3456384
Arellano, J. S., Molina, L. M., Rubio, A., and Alonso, J. A. (2000). Density functional study of adsorption of molecular hydrogen on graphene layers. J. Chem. Phys. 112:8114. doi: 10.1063/1.481411
Balandin, A. A., Ghosh, S., Bao, W., Calizo, I., Teweldebrhan, D., Miao, F., et al. (2008). Superior thermal conductivity of single-layer graphene. Nano Lett. 8, 902–907. doi: 10.1021/nl0731872
Bleakley, K., and Hu, P. (1999). A density functional theory study of the interaction between CO and O on a Pt Surface: CO/Pt (111), O/Pt (111), and CO/O/Pt (111). J. Am. Chem. Soc. 121, 7644–7652. doi: 10.1021/ja983363w
Chen, Z. W., Wen, Z., and Jiang, Q. (2017). Rational design of Ag38 cluster supported by graphdiyne for catalytic CO oxidation. J. Phys. Chem. C 121, 3463–3468. doi: 10.1021/acs.jpcc.6b12434
Delley, B. (2000). From molecules to solids with the Dmol3 approach. J. Chem. Phys. 113:7756. doi: 10.1063/1.1316015
Dupont, C., Jugnet, Y., and Loffreda, D. (2006).Theoretical evidence of PtSn alloy efficiency for CO oxidation. J. Am. Chem. Soc. 128, 9129–9136. doi: 10.1021/ja061303h
Dvorák, F., Camellone, M. F., Tovt, A., Tran, N. D., Negreiros, F. R., Vorokhta, M., et al. (2016). Creating single-atom Pt-ceria catalysts by surface step decoration. Nat. Commun. 7:10801. doi: 10.1038/ncomms10801
Gong, X. Q., Liu, Z. P., Raval, R., and Hu, P. (2004). A systematic study of CO oxidation on metals and metal oxides: density functional theory calculations. J. Am. Chem. Soc. 126, 8–9. doi: 10.1021/ja030392k
Grimme, S. (2006). Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 27, 1787–1799. doi: 10.1002/jcc.20495
Halgren, T. A., and Lipscomb, W. N. (1977). The synchronous-transit method for determining reaction pathways and locating molecular transition states. Chem. Phys. Lett. 49, 225–232. doi: 10.1016/0009-2614(77)80574-5
Henkelman, G., and Jonsson, H. (2000). A climbing image nudged elastic band method for finding saddle points and minimum energy paths. J. Chem. Phys. 113:9978. doi: 10.1063/1.1323224
Honkala, K., and Laasonen, K. (2000). Oxygen molecule dissociation on the Al (111) surface. Phys. Rev. Lett. 84:705. doi: 10.1103/PhysRevLett.84.705
Jiang, Q. G., Ao, Z. M., Li, S., and Wen, Z. (2014). Density functional theory calculations on the CO catalytic oxidation on Al-embedded graphene. RSC Adv. 4, 20290–20296. doi: 10.1039/C4RA01908C
Jiang, Q. G., Wu, W. C., Zhang, J. F., Ao, Z. M., Wu, Y. P., and Huang, H. J. (2016). Defections induced hydrogenation of silicene: a density functional theory calculation study. RSC Adv. 6, 69861–69868. doi: 10.1039/C6RA11885B
Jones, J., Xiong, H. F., Delariva, A. T., Peterson, E. J., Pham, H., Challa, S. R., et al. (2016). Thermally stable single-atom platinum-on-ceria catalysts via atom trapping. Science 353, 150–154. doi: 10.1126/science.aaf8800
Lee, C. G., Wei, X. D., Kysar, J. W., and Hone, J. (2008). Measurement of the elastic properties and intrinsic strength of monolayer graphene. Science 321, 385–388. doi: 10.1126/science.1157996
Leenaerts, O., Partoens, B., and Peeters, F. M. (2008). Adsorption of H2O, NH3, CO, NO2, and NO on graphene: a first-principles study. Phys. Rev. B 77:125416. doi: 10.1103/PhysRevB.77.125416
Li, Y. F., Zhou, Z., Yu, G. T., Chen, W., and Chen, Z. F. (2010). CO catalytic oxidation on iron-embedded graphene: computational quest for low-cost nanocatalysts. J. Phys. Chem. C 114, 6250–6254. doi: 10.1021/jp911535v
Liu, D. J (2007). CO oxidation on Rh (100): multisite atomistic lattice-gas modelling. J. Phys. Chem. C 111, 14698–14706. doi: 10.1021/jp071944e
Liu, S. G., and Huang, S. P. (2017). Theoretical insights into the activation of O2 by Pt single atom and Pt4 nanocluster on functionalized graphene support: critical role of Pt positive polarized charges. Carbon 115, 11–17. doi: 10.1016/j.carbon.2016.12.094
Liu, Z. P., Hu, P., and Alavi, A. (2002). Catalytic role of gold in gold-based catalysts: a density functional theory study on the CO oxidation on gold. J. Am. Chem. Soc. 124, 14770–14779. doi: 10.1021/ja0205885
Liu, Z. Y., He, T. W., Liu, K. K., Chen, W. G., and Tang, Y. N. (2017). Structural, electronic and catalytic performances of single-atom Fe stabilized by divacancy-nitrogen-doped graphene. RSC Adv. 7, 7920–7928. doi: 10.1039/C6RA28387J
Lu, Y. H., Zhou, M., Zhang, C., and Feng, Y. P. (2009). Metal-embedded graphene: a possible catalyst with high activity. J. Phys. Chem. C 113, 20156–20160. doi: 10.1021/jp908829m
Ma, L., Melander, M., Weckman, T., Laasonen, K., and Akola, J. (2016). CO oxidation on the Au15Cu15 cluster and the role of vacancies in the MgO(100) support. J. Phys. Chem. C 120, 26747–26758. doi: 10.1021/acs.jpcc.6b06876
Molina, L. M., and Hammer, B. (2005). The activity of the tetrahedral Au20 cluster: charging and impurity effects. J. Catal. 233, 399–404. doi: 10.1016/j.jcat.2005.04.037
Novoselov, K. S., Geim, A. K., Morozov, S. V., Jiang, D., Katsnelson, M. I., Grigorieva, I. V., et al. (2005). Two-dimensional gas of massless dirac fermions in graphene. Nature 438, 197–200. doi: 10.1038/nature04233
Novoselov, K. S., Geim, A. K., Morozov, S. V., Jiang, D., Zhang, Y., Dubonos, S. V., et al. (2004). Electric field effect in atomically thin carbon films. Science 306, 666–669. doi: 10.1126/science.1102896
Perdew, J. P., Burke, K., and Ernzerhof, M. (1996). Generalized gradient approximation made simple. Phys. Rev. Lett. 77:3865. doi: 10.1103/PhysRevLett.77.3865
Qiao, B. T., Liu, J. X., Wang, Y. G., Lin, Q. Q., Liu, X. Y., Wang, A. Q., et al. (2015). Highly efficient catalysis of preferential oxidation of CO in H2-rich stream by gold single-atom catalysts. ACS Catal. 5, 6249–6254. doi: 10.1021/acscatal.5b01114
Robertson, A. W., Allen, C. S., Wu, Y. A., He, K., Olivier, J., Neethling, J., et al. (2012). Spatial control of defect creation in graphene at the nanoscale. Nat. Commun. 3:1144. doi: 10.1038/ncomms2141
Robertson, A. W., Montanari, B., He, K., Kim, J., Allen, C. S., Wu, Y. A., et al. (2013). Dynamics of single Fe atoms in graphene vacancies. Nano Lett. 13, 1468–1475. doi: 10.1021/nl304495v
Roldán, A., Ricart, J. M., and Illas, F. (2009). Influence of the exchange-correlation potential on the description of the molecular mechanism of oxygen dissociation by Au nanoparticles. Theor. Chem. Acc. 123, 119–126. doi: 10.1007/s00214-009-0540-1
Saavedra, J., Whittaker, T., Chen, Z. F., Pursell, C. J., Rioux, R. M., and Chandler, B. D. (2016). Controlling activity and selectivity using water in the Au-catalysed preferential oxidation of CO in H2. Nat. Chem. 8, 584–589. doi: 10.1038/nchem.2494
Segall, M. D., Lindan, P. L. D., Probert, M. J., Pickard, C. J., Hasnip, P. J., Clark, S. J., et al. (2002). First-principles simulation: ideas, illustrations and the CASTEP code. J. Phys. 14, 2717–2744. doi: 10.1088/0953-8984/14/11/301
Song, E. H., Wen, Z., and Jiang, Q. (2011). CO catalytic oxidation on copper-embedded graphene. J. Phys. Chem. C 115, 3678–3683. doi: 10.1021/jp108978c
Tang, Y. A., Zhou, J. C., Shen, Z. G., Chen, W. G., Li, C. G., and Dai, X. Q. (2016). High catalytic activity for CO oxidation on single Fe atom stabilized in graphene vacancies. RSC Adv. 6, 93985–93996. doi: 10.1039/C6RA14476D
Tang, Y. N., Lu, Z. S., Chen, W. G., Li, W., and Dai, X. Q. (2015). Geometric stability and reaction activity of Pt clusters adsorbed graphene substrates for catalytic CO oxidation. Phys. Chem. Chem. Phys. 17, 11598–11608. doi: 10.1039/C5CP00052A
Tang, Y. N., Yang, Z. X., and Dai, X. Q. (2012). A theoretical simulation on the catalytic oxidation of CO on Pt/graphene. Phys. Chem. Chem. Phys. 14, 16566–16572. doi: 10.1039/c2cp41441d
Wang, Y. G., Cantu, D. C., Lee, M. S., Li, J., Glezakou, V. A., and Rousseau, R. (2016). CO oxidation on Au/TiO2: condition-dependent active sites and mechanistic pathways. J. Am. Chem. Soc. 138, 10467–10476. doi: 10.1021/jacs.6b04187
Xie, X. W., Li, Y., Liu, Z. Q., Haruta, M., and Shen, W. J. (2009). Low-temperature oxidation of CO catalysed by Co3O4 nanorods. Nature 458, 746–749. doi: 10.1038/nature07877
Young, D. C. (2001). Computational Chemistry: A Practical Guide for Applying Techniques to Real World Problems. New York, NY: Wiley.
Zhang, A. H., Zhu, J., and Duan, W. H. (2006). Oxidation of carbon monoxide on Rh (111): a density functional theory study. J. Chem. Phys. 124:234703. doi: 10.1063/1.2209239
Zhang, C. J., and Hu, P. (2001). CO oxidation on Pd (100) and Pd (111): a comparative study of reaction pathways and reactivity at low and medium coverages. J. Am. Chem. Soc. 123, 1166–1172. doi: 10.1021/ja002432f
Keywords: graphene, divacancy, Mn-embedded, CO oxidation, first principles calculations
Citation: Jiang Q, Zhang J, Ao Z, Huang H, He H and Wu Y (2018) First Principles Study on the CO Oxidation on Mn-Embedded Divacancy Graphene. Front. Chem. 6:187. doi: 10.3389/fchem.2018.00187
Received: 06 January 2018; Accepted: 07 May 2018;
Published: 29 May 2018.
Edited by:
Philippe Vernoux, UMR5256 Institut de Recherches sur la Catalyse et l'Environnement de Lyon (IRCELYON), FranceReviewed by:
Xingang Li, Tianjin University, ChinaJose Luis Valverde, Universidad de Castilla-La Mancha, Spain
Copyright © 2018 Jiang, Zhang, Ao, Huang, He and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianfeng Zhang, jfzhang@hhu.edu.cn
Zhimin Ao, zhimin.ao@gdut.edu.cn