Next Article in Journal
Hydrogels Synthesized by Electron Beam Irradiation for Heavy Metal Adsorption
Previous Article in Journal
The Antifungal Activity of Functionalized Chitin Nanocrystals in Poly (Lactid Acid) Films
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Materials
Volume 10
Issue 5
10.3390/ma10050549
materials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

DFT Study of the Oxygen Reduction Reaction Activity on Fe−N4-Patched Carbon Nanotubes: The Influence of the Diameter and Length

by
Xin Chen
*,
Rui Hu
and
Fan Bai
The Center of New Energy Materials and Technology, College of Chemistry and Chemical Engineering, Southwest Petroleum University, Chengdu 610500, China
*
Author to whom correspondence should be addressed.
Materials 2017, 10(5), 549; https://doi.org/10.3390/ma10050549
Submission received: 6 April 2017 / Revised: 15 May 2017 / Accepted: 15 May 2017 / Published: 18 May 2017
(This article belongs to the Section Energy Materials)
Browse Figures
Versions Notes

Abstract

:
The influences of diameter and length of the Fe−N4-patched carbon nanotubes (Fe−N4/CNTs) on oxygen reduction reaction (ORR) activity were investigated by density functional theory method using the BLYP/DZP basis set. The results indicate that the stability of the Fe−N4 catalytic site in Fe−N4/CNTs will be enhanced with a larger tube diameter, but reduced with shorter tube length. A tube with too small a diameter makes a Fe−N4 site unstable in acid medium since Fe−N and C−N bonds must be significantly bent at smaller diameters due to hoop strain. The adsorption energy of the ORR intermediates, especially of the OH group, becomes weaker with the increase of the tube diameter. The OH adsorption energy of Fe−N4/CNT with the largest tube diameter is close to that on Pt(111) surface, indicating that its catalytic property is similar to Pt. Electronic structure analysis shows that the OH adsorption energy is mainly controlled by the energy levels of Fe 3d orbital. The calculation results uncover that Fe−N4/CNTs with larger tube diameters and shorter lengths will exhibit better ORR activity and stability.

1. Introduction

In recent years, the development of nonprecious metal catalysts with high oxygen reduction reaction (ORR) activity and stability has been a major focus of fuel cell research. As early as 1964, Jasinski discovered that metal−N4 macrocyclic compounds are able to catalyze the ORR [1]. Nevertheless, these metal macrocycle-based electrodes (such as Fe-phthalocyanine and Fe-porphyrin) have not shown yet any acceptable ORR efficiency and stability in practical use [2]. The main reason is that previous methods employed simple mixing of metal-macrocycle with carbon black, which resulted in poor mechanical contacts with carbon support and thus inferior electrical conductivity [3]. Subsequent studies have shown that pyrolysis of transition metal, carbon, and some nitrogen containing materials—such as organic complexes, nitrogen containing salts and even gaseous nitrogen precursors—can also produce active ORR sites [4]. It was proposed that the active sites are the Fe ions coordinated by four nitrogen atoms, which are incorporated to graphitic carbon and have the same local structure of Fe-porphyrin catalysts [5]. These Fe−N4/C catalytic sites in graphitic carbon matrix were demonstrated to have excellent ORR catalytic performance [2,6,7,8,9]. Recently, graphene- and carbon nanotube (CNT)-based ORR catalysts containing transition metal and nitrogen attracted much research interest due to their excellent ORR activities [10,11,12,13]. The related research indicates that the activity of metal−Nx/C catalytic sites in graphene matrix is mainly affected by metal types, the N-coordination number and matrix size [13]. However, for Fe−Nx catalytic sites in CNT matrix, little is known about the relationship between the activity and the catalytic structures. The ORR activity and stability of CNTs are highly related to the tube diameter and length, which can obviously influence the geometrical and electronic structures of CNTs [14,15,16,17]. However, there is little related research reported so far. A systemic investigation of this issue will be helpful to build a reasonable Fe−Nx/CNT model and to understand its catalytic activity towards ORR.
In this paper, we applied the density functional theory (DFT) method to study how the tube diameter and length influence the ORR activity and stability of Fe−N4/CNTs nano-catalysts. The theoretical calculations indicate that Fe−N4/CNT catalysts with large tube diameter and short length should have good ORR activity and stability.

2. Computational Methods and Models

The DFT calculations carried out in this study were based on GGA/BLYP [18] using the Amsterdam Density Functional (ADF) program package [19]. ADF calculations were run using the DZP basis set keeping the 1s core of C, N, and O, and 1s–3p core for iron atoms frozen. Note that we did not consider the d-band correction for Fe atom. For all stationary states, spin multiplicity was allowed to relax: possible geometries with varying spin states were carefully checked and the ground state was determined to be the one with the lowest electronic energy.
The models used in this work are armchair nanotubes terminated with C−H bonds, and they are denoted as Fe−N4 (N, N)-L. Here, (N, N) represents the diameter of the CNT and L indicates the length (N = 2–8, L = 9.8, 12.3, 14.7, 17.3, 19.7, and 22.1 Å, respectively). From Fe−N4 (2, 2)-L to Fe−N4 (8, 8)-L, the diameters are 2.8, 4.3, 5.6, 6.8, 8.2, 9.7, and 11 Å, respectively. From L = 9.8 to 22.1 Å, the length of the CNT was increased by adding a cell containing several six-membered rings (note that the entire findings in this work could only apply to the studied armchair nanotubes).
The adsorption energy of the ORR species is an important reference point for determining the activity and stability of an electrocatalyst. In this paper, the adsorption energy (AE) of the ORR species is calculated using the equation: AE(molecule) = E(catalyst−molecule system) − E(catalyst) − E(molecule). Therefore, a negative adsorption energy suggests that the ORR species would be energetically favorable to be adducted to the Fe−N4 (N, N)-L catalysts.

3. Results and Discussion

3.1. The Stability of Fe−N4 Catalytic Sites in Acid Medium

The optimized configurations of Fe−N4/CNTs with different tube diameters (the tube length is 9.8 Å) and lengths (taken Fe−N4 (4, 4)-L as an example) are shown in Figure 1. In order to evaluate the stability of these Fe−N4 catalytic sites in acid medium, we calculated the energies required (ΔE) to remove the Fe ion from the bulk Fe−N4/CNT structures using the following equation
Fe−N4 (N, N)-L + 2[H3O+(H2O)2] → 2H−N4 (N, N)-L + [Fe(H2O)6]2+
Obviously, a larger ΔE denotes a more stable structure of the Fe−N4 site in the catalytic process (note that this is only a thermodynamic approximation and does not address the issue on kinetics of bond breaking which are also important for stability, as well as that we only consider the stability of the Fe−N4 site, not the corrosion process of bulk carbons).
The calculated ΔE values are shown in Figure 2. At a given length (such as 9.8 Å), as the tube diameter increases, the stability of Fe−N4 (N, N)-9.8 (N = 2–8) increases, as clearly shown in Figure 2a. However, as seen from the decreased ΔE values, the stability decreases with the increasing CNT length, as shown in Figure 2b. The ΔE is only +0.12 eV for Fe−N4 (2, 2)-9.8, indicating that it is not very stable in acid medium. This is mainly because the tube diameter is so small that the resulting hoop strain could significantly bend the Fe−N and C−N bonds. Therefore, the catalytic Fe site would be partially exposed to the external acid environment, leading to its instability. Except for Fe−N4 (2, 2)-9.8, other catalytic structures with tube length of 9.8 Å have higher ΔE values, suggesting they are stable in acid solution. Furthermore, for Fe−N4 (N, N)-L with a larger tube diameter (for example, the Fe−N4 (5, 5)-L), although the stability decreases with the increasing of tube length, the ΔE values are always above +2.2 eV, meaning that they are still stable in acid medium. It should be noted that the ΔE values in Figure 2b do not change monotonically with the increasing tube length. The possible reason for this is that the local geometrical structure (such as the average bond length of four Fe−N bonds) and surface electronic structure of Fe−N4/CNTs might be slightly affected by the tube length. Furthermore, although the curves are not completely monotonic, a general conclusion can also be made from Figure 2b that the ΔE values of longer structures are lower than those of shorter ones.

3.2. Adsorption of ORR Species

The adsorption energy (AE) of the ORR species is an important criterion in assessing the activity of a catalyst. Platinum-based materials are well known for their high catalytic activities for ORR. We can evaluate whether or not Fe−N4 (N, N)-L based catalysts have superior catalytic properties by comparing their adsorption energies with Pt.
All the adsorption energies of the ORR species on Fe−N4 (N, N)-9.8 (N = 2–8) are shown in Figure 3, and only the adsorption configurations on Fe−N4 (4, 4)-9.8 are shown in Figure 4 for the sake of clarity. The AE(O2) on Fe−N4 (2, 2)-9.8 and Fe−N4 (3, 3)-9.8 are respectively −2.43 and −1.00 eV, while others are in the range of −0.63 to −0.76 eV. The experimentally determined low-coverage adsorption energy of O2 on Pt(111) is −0.3 to −0.5 eV [20,21,22], and theoretical values range between −0.41 and −1.04 eV depending on the methods and models used [23,24,25,26]. In general, the O2 adsorption energy on an ideal catalytic material should be as small as possible, but large enough to prevent O2 from drifting away or desorbing from the catalytic center [27]. The adsorption energy of O2 on Fe−N4 (2, 2)-9.8 is −2.43 eV, which is much stronger than that on the Pt(111) surface. As is well known, during the operation of the fuel cell, the energy loss of this non-electron-transfer step is unavailable for electrical work. The waste reaction energy at this step will lead to the reduction of the energy at those electron transfer steps. In this case, a large overpotential is inevitable. Furthermore, considering its instable nature in acid solution, a conclusion can be made that the Fe−N4 (2, 2)-9.8 is not a good catalyst for ORR. Conversely, other tubes with a length of 9.8 Å have suitable O2 adsorption energy, and are suitable for the ORR.
The AE(OOH) on Fe−N4 (2, 2)-9.8, as shown in Figure 3, is −2.56 eV, which is very close to the value of O2 adsorption energy. Generally, for an effective electrocatalyst, the heats of formation of the OOH should be higher than the heats of formation of the adsorbed O2 on the surface, otherwise the ORR will not be an energetically favorable process [27]. Therefore, in view of this, the Fe−N4 (2, 2)-9.8 nanomaterial is again proven to be an ineffective electrocatalyst for ORR. Thus, the ORR catalyzed by Fe−N4 (2, 2)-9.8 is no longer considered in the following sections except when being specially pointed out. Interestingly, the AE(OOH) on other tubes almost has no change with the increasing tube diameter. The −1.60 eV adsorption energy seems to be unfavorable for ORR since it is slightly higher than the reported theoretical value of −1.06 to −1.16 eV for OOH adsorbed on Pt(111) surface [23,28]. However, as is well known, Pt is not a perfect catalyst for the ORR because it bonds OOH too weakly while it bonds O and OH too strongly [29,30]. Therefore, a higher OOH adsorption energy than Pt will lead to a larger heat loss during the first ORR step, which may facilitate the electron transfer. This conclusion is further confirmed by the calculated reaction energy for each ORR step, which we will show later.
The adsorption energy of atomic O, similar to AE(OOH), does not change significantly with the increasing of tube diameter. However, the situation is quite different for OH adsorption. As is clearly shown in Figure 3, the AE(OH) gradually becomes weaker with the increasing tube diameter. Let us compare these OH adsorption results with those on Pt(111) surface. Previous DFT calculations of OH adsorption on Pt(111) surface led to AE(OH) = −2.26 to −2.45 eV [23,31]. DFT calculations performed on a 35 atom Pt cluster imitating the Pt(111) surface led to AE(OH) = −2.06 eV [32]. Therefore, the AE(OH) on the tubes with smaller diameter are relatively too strong, while on those tubes with larger diameters—for example, the Fe−N4 (8, 8)-9.8—are too close to those on Pt(111) surface.
The variation of AE(OH) as a function of the tube length is shown in Figure 5. Unlike the case shown in Figure 3, it seems that there are no systematic changes with the increasing of tube length. The AE(OH) on Fe−N4 (4, 4)-L and Fe−N4 (5, 5)-L does not change much in the studied range of tube length, while on Fe−N4 (3, 3)-L, the change becomes obvious. The weakest AE(OH) on Fe−N4 (3, 3)-L is −2.87 eV, which is still much stronger than that on Pt(111) surface. The strong OH adsorption, as is well known, will increase the difficulty of removing OH from the catalyst’s surface and therefore lead to a high overpotential [33,34].
The calculated AE(H2O) on the Fe−N4 (N, N)-9.8 (N = 3–8) is in the range of −0.28 to −0.36 eV. These results are very close to the experimental values of −0.43 to −0.65 eV [35] and theoretical values of −0.22 to −0.60 eV on Pt(111) surface [23,36]. In general, the calculated adsorption energies of the ORR species on the above catalysts, especially of O and OH, are weaker than that on Fe−N4 embedded in the graphene [37] (note that this Fe−N4 was embedded in the bulk of graphene, not the edge). This may be attributed to the curvature effect of CNT [38]. Besides, we found that for all the studied Fe−N4 (N, N)-L (N = 3–8) nano-materials, the AE(H2O) are much weaker than their corresponding AE(O2). Generally, for an effective ORR catalyst, the H2O molecule should be adsorbed more weakly than O2, so that the catalytic cycle could repeat easily [39]. Thus, in principle, all the studied Fe−N4 (N, N)-L (N = 3–8) nano-materials have catalytic activities for ORR. Nevertheless, just as discussed above, only the tubes with larger diameters have AE(OH) close to the one on the Pt(111) surface. Since the OH group is regarded as a poison on the catalyst’s surface, stronger adsorption of OH on the tubes with smaller diameters will lead to higher overpotential for the ORR. Therefore, based on the above analysis, the Fe−N4 (N, N)-L based catalysts with large tube diameter and short length should have higher catalytic ORR activities.
In order to further support these conclusions, the relative energies for each electron transfer step are calculated, as shown in Table 1. The values on the Pt(111) are cited from the published literature [40]. The rate-determining step (RDS) is considered as the one with the smallest energy change. The higher the energy difference of the RDS is, the better the catalyst. Obviously, the RDS for all the structures is the reduction of OH group. From Fe−N4 (2, 2)-9.8 to Fe−N4 (4, 4)-9.8, the energy change values of RDS are positive, indicating these steps are energetically unfavorable. From Fe−N4 (5, 5)-9.8, due to the decreased OH adsorption, all the energy changes of the RDS become exothermic, suggesting the spontaneous nature of the process. As is expected, Fe−N4 (8, 8)-9.8 possesses the largest energy change of RDS, indicating it possesses the best ORR activity among all the screened catalysts. This conclusion is consistent with the analysis that the decreased OH adsorption is favorable for ORR, and the catalysts with larger tube diameters should have higher catalytic ORR activities.

3.3. Electronic Structural Effect on the ORR Activity

It is well known that the ORR activity is largely controlled by the electronic structure of the catalysts. DFT method is a good theoretical approach to investigate the relationship between the electronic structure and catalytic activity. Previous studies on periodic metal nanostructures, such as Pt-based materials [41,42] and other noble metals [43], show that the ORR activity is highly determined by the d-band center of the catalysts. For non-periodic molecular catalysts, such as various metal−Nx/C materials, the occupied metal 3d orbital is proven to be of great significance in describing the catalytic activity [44]. In this paper, we calculated the energy levels of Fe 3d HOMO (the highest occupied molecular orbital) for both spin-up and spin-down electrons to investigate how orbital levels affect the ORR activity. The obtained results are shown in Figure 6. The energy levels of Fe 3d HOMO, especially for spin-up electrons, have been decreased with the increasing of tube diameter. Since the change of electronic structure is directly related to the variation of chemical reactivity, the adsorption energies of some intermediates will be affected. As shown in Figure 3, the AE(OH) is most affected by the electronic structure. Higher energy levels of Fe 3d HOMO could lead to stronger adsorption interactions. On the other hand, an appropriate HOMO orbital can induce an appropriate adsorption energy, which is favorable for ORR.

4. Conclusions

In this paper, the influences of diameter (2.8–11 Å) and length (9.8–22.1 Å) of carbon nanotube on Fe−N4/CNTs catalytic activities for ORR were investigated by theoretical calculations at the BLYP/DZP level of theory. The results indicate that the stability of the Fe−N4 site increases with the increase of tube diameter, but reduces with the increase of tube length. Besides, the tube with a too small diameter is not very stable in acid medium due to the fact that, at smaller diameters, the Fe−N and C−N bonds must be significantly bent due to hoop strain. The adsorption energy of the ORR species, especially of the OH group, becomes weaker with the increase of tube diameter. The tube with the largest diameter has the AE(OH) closest to that on the Pt(111) surface, suggesting a catalytic property similar to Pt. Electronic structure analysis shows that the AE(OH) is mainly controlled by the energy levels of Fe 3d orbital. It suggests that the Fe−N4 (N, N)-L nano-catalysts with large tube diameters and short lengths should have good ORR activity and stability.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (No. 51602270). The authors thank the Supercomputing Center of the Chinese Academy of Sciences and Beijing University of Technology for providing the computational resources and software.

Author Contributions

Xin Chen conceived and designed the study; Xin Chen and Rui Hu performed the calculations; Rui Hu and Fan Bai analyzed the data; Xin Chen wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jasinski, R. A new fuel cell cathode catalyst. Nature 1964, 201, 1212–1213. [Google Scholar] [CrossRef]
  2. Lee, D.H.; Lee, W.J.; Lee, W.J.; Kim, S.O.; Kim, Y.H. Theory, synthesis, and oxygen reduction catalysis of Fe-porphyrin-like carbon nanotube. Phys. Rev. Lett. 2011, 106, 175502. [Google Scholar] [CrossRef] [PubMed]
  3. Chen, R.; Li, H.; Chu, D.; Wang, G. Unraveling oxygen reduction reaction mechanisms on carbon-supported Fe-phthalocyanine and Co-phthalocyanine catalysts in alkaline solutions. J. Phys. Chem. C 2009, 113, 20689–20697. [Google Scholar] [CrossRef]
  4. Chen, Z.; Higgins, D.; Yu, A.; Zhang, L.; Zhang, J. A review on non-precious metal electrocatalysts for PEM fuel cells. Energy Environ. Sci. 2011, 4, 3167–3192. [Google Scholar] [CrossRef]
  5. Orellana, W. Catalytic properties of transition metal-N4 moieties in graphene for the oxygen reduction reaction: Evidence of spin-dependent mechanisms. J. Phys. Chem. C 2013, 117, 9812–9818. [Google Scholar] [CrossRef]
  6. Lefèvre, M.; Proietti, E.; Jaouen, F.; Dodelet, J.P. Iron-based catalysts with improved oxygen reduction activity in polymer electrolyte fuel cells. Science 2009, 324, 71–74. [Google Scholar] [CrossRef] [PubMed]
  7. Wu, G.; More, K.L.; Johnston, C.M.; Zelenay, P. High-performance electrocatalysts for oxygen reduction derived from polyaniline, iron, and cobalt. Science 2011, 332, 443–447. [Google Scholar] [CrossRef] [PubMed]
  8. Liu, R.; Malotki, C.; Arnold, L.; Koshino, N.; Higashimura, H.; Baumgarten, M.; Müllen, K. Triangular trinuclear metal-N4 complexes with high electrocatalytic activity for oxygen reduction. J. Am. Chem. Soc. 2011, 133, 10372–10375. [Google Scholar] [CrossRef] [PubMed]
  9. Byon, H.R.; Suntivich, J.; Shao-Horn, Y. Graphene-based non-noble-metal catalysts for oxygen reduction reaction in acid. Chem. Mater. 2011, 23, 3421–3428. [Google Scholar] [CrossRef]
  10. Li, Y.; Zhou, W.; Wang, H.; Xie, L.; Liang, Y.; Wei, F.; Idrobo, J.C.; Pennycook, S.J.; Dai, H. An oxygen reduction electrocatalyst based on carbon nanotube-graphene complexes. Nat. Nanotechnol. 2012, 7, 394–400. [Google Scholar] [CrossRef] [PubMed]
  11. Kattel, S.; Atanassov, P.; Kiefer, B. Density functional theory study of Ni−Nx/C electrocatalyst for oxygen reduction in alkaline and acidic medium. J. Phys. Chem. C 2012, 116, 17378–17383. [Google Scholar] [CrossRef]
  12. Kattel, S.; Atanassov, P.; Kiefer, B. Catalytic activity of Co−Nx/C electrocatalysts for oxygen reduction reaction: A density functional theory study. Phys. Chem. Chem. Phys. 2013, 15, 148–153. [Google Scholar] [CrossRef] [PubMed]
  13. Chen, X.; Li, F.; Zhang, N.; An, L.; Xia, D. Mechanism of oxygen reduction reaction catalyzed by Fe (Co)−Nx/C. Phys. Chem. Chem. Phys. 2013, 15, 19330–19336. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, Y.F.; Liu, Z.F. Oxidation of zigzag carbon nanotubes by singlet O2: Dependence on the tube diameter and the electronic structure. J. Phys. Chem. B 2004, 108, 11435–11441. [Google Scholar] [CrossRef]
  15. Zurek, E.; Autschbach, J. Density functional calculations of the 13C NMR chemical shifts in (9, 0) single-walled carbon nanotubes. J. Am. Chem. Soc. 2004, 126, 13079–13088. [Google Scholar] [CrossRef] [PubMed]
  16. Zhao, J.; Balbuena, P.B. Structural and reactivity properties of finite length cap-ended single-wall carbon nanotubes. J. Phys. Chem. A 2006, 110, 2771–2775. [Google Scholar] [CrossRef] [PubMed]
  17. Hu, X.; Liu, C.; Wu, Y.; Zhang, Z. Density functional theory study on nitrogen-doped carbon nanotubes with and without oxygen adsorption: The influence of length and diameter. New J. Chem. 2011, 35, 2601–2606. [Google Scholar] [CrossRef]
  18. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef]
  19. Velde, G.; Bickelhaupt, F.M.; Baerends, E.J.; Fonseca Guerra, C.; van Gisbergen, S.J.A.; Snijders, J.G.; Ziegler, T. Chemistry with ADF. J. Comput. Chem. 2001, 22, 931–967. [Google Scholar] [CrossRef]
  20. Gland, J.L.; Sexton, B.A.; Fisher, G.B. Oxygen interactions with the Pt(111) surface. Surf. Sci. 1980, 95, 587–602. [Google Scholar] [CrossRef]
  21. Campbell, C.T.; Ertl, G.; Kuipers, H.; Segner, J. A molecular beam study of the adsorption and desorption of oxygen from a Pt(111) surface. Surf. Sci. 1981, 107, 220–236. [Google Scholar] [CrossRef]
  22. Nolan, P.D.; Lutz, B.R.; Tanaka, P.L.; Davis, J.E.; Mullins, C.B. Molecularly chemisorbed intermediates to oxygen adsorption on Pt (111): A molecular beam and electron energy-loss spectroscopy study. J. Chem. Phys. 1999, 111, 3696–3704. [Google Scholar] [CrossRef]
  23. Sha, Y.; Yu, T.H.; Liu, Y.; Merinov, B.V.; Goddard, W.A. Theoretical study of solvent effects on the platinum-catalyzed oxygen reduction reaction. J. Phys. Chem. Lett. 2010, 1, 856–861. [Google Scholar] [CrossRef]
  24. Qi, L.; Qian, X.; Li, J. Near neutrality of an oxygen molecule adsorbed on a Pt (111) surface. Phys. Rev. Lett. 2008, 101, 146101. [Google Scholar] [CrossRef] [PubMed]
  25. Hyman, M.P.; Medlin, J.W. Effects of electronic structure modifications on the adsorption of oxygen reduction reaction intermediates on model Pt (111)-alloy surfaces. J. Phys. Chem. C 2007, 111, 17052–17060. [Google Scholar] [CrossRef]
  26. Feng, Y.; Li, F.; Hu, Z.; Luo, X.; Zhang, L.; Zhou, X.F.; Wang, H.T.; Xu, J.J.; Wang, E.G. Tuning the catalytic property of nitrogen-doped graphene for cathode oxygen reduction reaction. Phys. Rev. B 2012, 85, 155454. [Google Scholar] [CrossRef]
  27. Lyalin, A.; Nakayama, A.; Uosaki, K.; Taketsugu, T. Theoretical predictions for hexagonal BN based nanomaterials as electrocatalysts for the oxygen reduction reaction. Phys. Chem. Chem. Phys. 2013, 15, 2809–2820. [Google Scholar] [CrossRef] [PubMed]
  28. Roques, J.; Anderson, A.B. Pt3Cr (111) alloy effect on the reversible potential of OOH (ads) formation from O2 (ads) relative to Pt (111). J. Fuel Cell Sci. Technol. 2005, 2, 86–93. [Google Scholar] [CrossRef]
  29. Chen, X.; Chen, S.; Wang, J. Screening of catalytic oxygen reduction reaction activity of metal-doped graphene by density functional theory. Appl. Surf. Sci. 2016, 379, 291–295. [Google Scholar] [CrossRef]
  30. Chen, X.; Li, M.; Yu, Z.; Ke, Q. A comparative DFT study of oxygen reduction reaction on mononuclear and binuclear cobalt and iron phthalocyanines. Russ. J. Phys. Chem. A 2016, 90, 2413–2417. [Google Scholar] [CrossRef]
  31. Han, B.C.; Miranda, C.R.; Ceder, G. Effect of particle size and surface structure on adsorption of O and OH on platinum nanoparticles: A first-principles study. Phys. Rev. B 2008, 77, 075410. [Google Scholar] [CrossRef]
  32. Jacob, T.; Goddard, W.A. Water formation on Pt and Pt-based alloys: A theoretical description of a catalytic reaction. ChemPhysChem 2006, 7, 992–1005. [Google Scholar] [CrossRef] [PubMed]
  33. Chen, X.; Chang, J.; Yan, H.; Xia, D. Boron nitride nanocages as high activity electrocatalysts for oxygen reduction reaction: Synergistic catalysis by dual active sites. J. Phys. Chem. C 2016, 120, 28912–28916. [Google Scholar] [CrossRef]
  34. Chen, X. Graphyne nanotubes as electrocatalysts for oxygen reduction reaction: The effect of doping elements on the catalytic mechanisms. Phys. Chem. Chem. Phys. 2015, 17, 29340–29343. [Google Scholar] [CrossRef] [PubMed]
  35. Fisher, G.B.; Gland, J.L. The interaction of water with the Pt(111) surface. Surf. Sci. 1980, 94, 446–455. [Google Scholar] [CrossRef]
  36. Keith, J.A.; Jacob, T. Theoretical studies of potential-dependent and competing mechanisms of the electrocatalytic oxygen reduction reaction on Pt (111). Angew. Chem. Int. Ed. 2010, 49, 9521–9525. [Google Scholar] [CrossRef] [PubMed]
  37. 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]
  38. Chai, G.L.; Hou, Z.; Shu, D.J.; Ikeda, T.; Terakura, K. Active sites and mechanisms for oxygen reduction reaction on nitrogen-doped carbon alloy catalysts: Stone−Wales defect and curvature effect. J. Am. Chem. Soc. 2014, 136, 13629–13640. [Google Scholar] [CrossRef] [PubMed]
  39. Vayner, E.; Anderson, A.B. Theoretical predictions concerning oxygen reduction on nitrided graphite edges and a cobalt center bonded to them. J. Phys. Chem. C 2007, 111, 9330–9336. [Google Scholar] [CrossRef]
  40. Chen, X.; Sun, F.; Chang, J. Cobalt or nickel doped SiC nanocages as efficient electrocatalyst for oxygen reduction reaction: A computational prediction. J. Electrochem. Soc. 2017, 164, F616–F619. [Google Scholar] [CrossRef]
  41. Stamenkovic, V.; Mun, B.S.; Mayrhofer, K.J.J.; Ross, P.N.; Markovic, N.M.; Rossmeisl, J.; Greeley, J.; Norskov, J.K. Changing the activity of electrocatalysts for oxygen reduction by tuning the surface electronic structure. Angew. Chem. 2006, 118, 2963–2967. [Google Scholar] [CrossRef]
  42. Stamenkovic, V.R.; Fowler, B.; Mun, B.S.; Wang, G.; Ross, P.N.; Lucas, C.A.; Markovic, N.M. Improved oxygen reduction activity on Pt3Ni(111) via increased surface site availability. Science 2007, 315, 493–497. [Google Scholar] [CrossRef] [PubMed]
  43. Lima, F.H.B.; Zhang, J.; Shao, M.H.; Sasaki, K.; Vukmirovic, M.B.; Ticianelli, E.A.; Adzic, R.R. Catalytic activity−d-band center correlation for the O2 reduction reaction on platinum in alkaline solutions. J. Phys. Chem. C 2007, 111, 404–410. [Google Scholar] [CrossRef]
  44. Chen, X. Oxygen reduction reaction on cobalt–(n)pyrrole clusters from DFT studies. RSC Adv. 2016, 6, 5535–5540. [Google Scholar] [CrossRef]
Materials 10 00549 g001 550
Figure 1. Optimized structures of Fe−N4 site in CNTs.
Figure 1. Optimized structures of Fe−N4 site in CNTs.
Materials 10 00549 g001
Materials 10 00549 g002 550
Figure 2. Reaction energies required for (a) Fe−N4 (N, N)-9.8 with different tube diameters (N = 2–8, the CNT length is 9.8 Å); and (b) Fe−N4 (N, N)-L with different tube lengths (N = 3–5, L = 9.8–22.1 Å).
Figure 2. Reaction energies required for (a) Fe−N4 (N, N)-9.8 with different tube diameters (N = 2–8, the CNT length is 9.8 Å); and (b) Fe−N4 (N, N)-L with different tube lengths (N = 3–5, L = 9.8–22.1 Å).
Materials 10 00549 g002
Materials 10 00549 g003 550
Figure 3. Calculated adsorption energies of the ORR species.
Figure 3. Calculated adsorption energies of the ORR species.
Materials 10 00549 g003
Materials 10 00549 g004 550
Figure 4. Calculated adsorption configurations of ORR species on Fe−N4 (4, 4)-9.8.
Figure 4. Calculated adsorption configurations of ORR species on Fe−N4 (4, 4)-9.8.
Materials 10 00549 g004
Materials 10 00549 g005 550
Figure 5. Calculated adsorption energies of OH with the increasing of tube length.
Figure 5. Calculated adsorption energies of OH with the increasing of tube length.
Materials 10 00549 g005
Materials 10 00549 g006 550
Figure 6. Calculated energy levels of Fe 3d HOMO with the increasing of tube diameter.
Figure 6. Calculated energy levels of Fe 3d HOMO with the increasing of tube diameter.
Materials 10 00549 g006
Table
Table 1. Calculated reaction energy changes (unit: eV) for each electron transfer step on Fe−N4 (N, N)-9.8 (N = 3–8) (* adsorbed species).
Table 1. Calculated reaction energy changes (unit: eV) for each electron transfer step on Fe−N4 (N, N)-9.8 (N = 3–8) (* adsorbed species).
Reaction Step(2, 2)(3, 3)(4, 4)(5, 5)(6, 6)(7, 7)(8, 8)Pt(111)
O2 + H+ + e− → * OOH−2.50−1.59−1.55−1.55−1.55−1.55−1.54−1.02
* OOH + H+ + e− → * O + H2O−1.89−1.90−1.85−1.86−1.87−1.88−1.88−2.01
* O + H+ + e− → * OH−1.26−1.07−1.09−0.99−0.93−0.88−0.85−0.77
* OH + H+ + e− → H2O+0.88+0.21+0.28−0.37−0.42−0.46−0.50−0.88

Share and Cite

MDPI and ACS Style

Chen, X.; Hu, R.; Bai, F. DFT Study of the Oxygen Reduction Reaction Activity on Fe−N4-Patched Carbon Nanotubes: The Influence of the Diameter and Length. Materials 2017, 10, 549. https://doi.org/10.3390/ma10050549

AMA Style

Chen X, Hu R, Bai F. DFT Study of the Oxygen Reduction Reaction Activity on Fe−N4-Patched Carbon Nanotubes: The Influence of the Diameter and Length. Materials. 2017; 10(5):549. https://doi.org/10.3390/ma10050549

Chicago/Turabian Style

Chen, Xin, Rui Hu, and Fan Bai. 2017. "DFT Study of the Oxygen Reduction Reaction Activity on Fe−N4-Patched Carbon Nanotubes: The Influence of the Diameter and Length" Materials 10, no. 5: 549. https://doi.org/10.3390/ma10050549

APA Style

Chen, X., Hu, R., & Bai, F. (2017). DFT Study of the Oxygen Reduction Reaction Activity on Fe−N4-Patched Carbon Nanotubes: The Influence of the Diameter and Length. Materials, 10(5), 549. https://doi.org/10.3390/ma10050549

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