The calculated properties of dioxygen O2 adsorption and reduction on the N-doped short CNTs. In the table, O-C is the shortest O-C distance; O-O is the O-O bond length, and MM is the dioxygen magnetic moment in µB. In the gas state, O2 has a magnetic moment of 2µB and a bond length of 1.22 Å from our first-principles DFT calculation.
The electronic and materials properties of carbon nanotubes (CNTs), like those of recent discovered graphene and earlier found fullerenes, attracted lots of research interests due to their appealing applications1 in the field of molecular electronics or for the refinement of materials, such as antistatic paints and shieldings, sensor, and catalytic functionality in fuel cells, etc. Like its two/three dimensional siblings graphene and graphite, CNTs consist
In a 2009 publication in
To further understand the properties of nitrogen doped carbon nanotubes as a metal-free electrocatalytic electrode for dioxygen reduction, Yang
2. Computational method
The first principles DFT calculations with the projector augmented wave (PAW) method12,13 were performed. The relativistic effect was included in the calculations. The Vienna Ab-initio Simulation Package (VASP)14,15,16,17 was used in the simulations. The exchange-correlation interaction potentials of the many electron system both in local density approximation (LDA) and in the generalized gradient
3. Optimized carbon nanotube structures
The synthesized nitrogen-containing carbon nanotubes exhibit a bamboo-like structure.4,5,7 Distinctively, a great part of the bamboo-like nano-structures consists of open ends of relatively short CNTs.18,19,20Individual short CNTs are weakly stacked one on top of the other to create a long nano-fiber. The observed bamboo-like structure along the nano-fiber length can be attributed to the integration of nitrogen into the graphitic structure, altering the nanotube surface from straight cylinder geometry. In order to simulate the open structure of the short carbon nanotubes, a short piece (10, 0) CNT with 16 atomic rings that have 10 carbon atoms per ring was used. There are a total of 160 carbon atoms in the short piece CNT that are included in the full self-consistent DFT calculations. The CNT has a diameter of 7.84 Å, through the center of the carbon atoms, and a total length of about 15.5 Å after the optimization of the structure. Different from the previous calculations, which utilized hydrogen atoms to artificially edge-saturate those CNTs to facility the convergence of the DFT computations,7 the open-edge of the short piece CNT was not saturated using hydrogen or other atoms. Wei
In this work, nitrogen-doped CNTs either at open-edge sites or in CNT walls was simulated. In addition, the open-edges of the short piece CNTs play an important role in the nitrogen doping, as demonstrated by the previous experimental work, 4,5 as well as by the first-principles calculations as discussed late. In the first-principles self-consistent DFT calculations, all atomic positions are allowed to be fully relaxed. The calculated results show that the carbon atoms on the open-edge of the short CNT have a substantial relaxation, due to the existence of some dangling bonds of the carbon atoms at the open-edge sites. The calculated C-C bond length of the carbon atoms on the open-edge to the atoms of the second atomic ring away from the open-edge of the CNT is 1.39 Å, which is shorter than the C-C bond length of 1.41 to 1.42 Å in the inner atomic rings of the short CNT. This reduction of the C-C bond length for the carbon atoms at the open-edge sites is mainly due to the absence of atoms on the empty side. The calculated bond length of the carbon atoms from the second to third atomic ring is 1.44 Å, which is slightly larger than those in the inner rings (such as those from the third to fourth atomic rings), away from the open-edge of the short CNT. The results indicate that the effect of the relaxation of the third atomic ring from the open-edge of the CNTs and other inner rings would be insignificant.
The spin-polarized DFT calculations for the short piece CNT were also performed. It is interesting to see that the carbon atoms on the open-edge of the CNT possess a magnetic moment of about 0.59 µB/atom. Similar spin polarization effect was also shown on a recent study of Möbius graphene nanoribbons.22 Other carbon atoms away from the open-edge and in the inner wall of the CNT do not have a noticeable magnetic moment from our computational results. This property of the carbon atoms on the open-edge with a magnetic moment is attributed to the existence of some dangling bonds. In the next sub-sections, it is demonstrate that the absorption of other atoms to the open-edge of the CNT will reduce the magnetic moment of the carbon atoms, because the dangling bond is reduced.
The nitrogen substitutionally doped carbon nanotubes, i. e., a nitrogen atom replacing a carbon atom in the short CNTs were then calculated. Four types of doping configurations were considered. They are,
In the following discussions, the properties of the nitrogen doping on the open-edges of the short CNTs are introduced. The nitrogen atom on the open-edge of the CNT bonds to two carbon atoms of the second atomic ring. The two N-C bond lengths are 1.31 Å and 1.34 Å, respectively. Because we only included one nitrogen atom substituting one carbon atom on the open-edge, the short CNT in the simulation did not have a cylindrical symmetry. The cross-section of the CNT is distorted and formed to nearly triangular shape by the N doping that can be seen from Figure 1(b). The distortion of the CNT cross-section is also partly attributed to the formation of new C-C bond of about 1.49 to 1.55 Å on the open-edge of the CNT, reducing the dangling bonds of the carbon atoms on the open-edge and forming in more stable
The effective charge (Bader charge) and the charge transfer of the N and C atoms in the doped CNT were also calculated. The nitrogen atom on the open-edge of the CNT obtained about 3.0 electrons from the nearest neighboring carbon atoms, i.e., turning to N-3.0. The C atom with a N-C bond length of 1.31 Å loses about 1.6 electrons to N atom; and the C atom with a N-C bond length of 1.34 Å loses about 1.3 electrons to the N atom. The calculated effective charge may have a computational uncertainty within an estimated range of about 0.1 electrons. The remaining of the charge transfer of about 0.1 electrons to the N atom can be attributed to the second nearest neighboring C atoms or to the computational uncertainty.
The spin-polarized DFT calculation for the short piece CNT with a doping N atom on the open-edge was then calculated. The results show that the nitrogen atom did not have a noticeable magnetic moment. In addition, some of the carbon atoms around the distorted locations of the open-edge of the short CNT have a strong relaxation and lose their magnetic moments. Only the carbon atoms on the open-edges of CNTs, which still have a dangling bond, possess a magnetic moment from 0.1 to 0.5 µB/atom, depending on their local environments.
4. Dioxygen adsorption on one nitrogen open edge doped CNT(10, 0)
The first-principles DFT calculations to study dioxygen O2 adsorption and reduction on the nitrogen edge-doped carbon nanotubes were performed. The calculated results show that the O2 can chemisorbed on a site close to the nitrogen-carbon complex with a tilted slant away from the nitrogen atom and more close to a carbon atom (Fig. 2a and b, Pauling adsorption model). The O-O bond length of the dioxygen on the end-on Pauling site increased to 1.33 Å, which is longer than the calculated bond-length of 1.22 Å in free gas O2 molecule state, see Table 1. The O(1)-C and O(2)-N distances for the Pauling adsorption are 1.31 Å and 2.64 Å, respectively. The dioxygen absorbed on the Pauling site obtained partial electrons from carbon-nitrogen complex. The calculated effective charges of two oxygen are O(1)-0.99 and O(2)-0.24, respectively. The effective charge of the nitrogen atom becomes N-2.9 after the dioxygen adsorption, which is not much different from that before the dioxygen adsorption. The electrons obtained by the dioxygen are mainly transferred from the neighboring carbon atoms. The adsorption energy of the dioxygen on the Pauling site is about 1.53 eV/atom. The magnetic momentum of the dioxygen at the Pauling site is greatly reduced from the gas state magnetic momentum of 2 µB to 0.027 µB. The carbon atom that bonded to O(1) loses its magnetic moment to nearly zero, from about 0.5 µB before the dioxygen adsorption. The nitrogen atom still has no noticeable magnetic moment.
The dioxygen adsorption and reduction at two other locations on the open-edge of the nitrogen doped carbon nanotubes was studied. The first case is the dioxygen adsorption on the opposite side of the CNT open-section away from nitrogen, which is illustrated in Figure 2 (c) and (d). The dioxygen is stabilized in a long bridge site on the open-edge of the CNT. The O-O bond length of the dioxygen on the long bridge site is 1.57 Å, which is even larger than the O-O bond length of the dioxygen at the end-on Pauling site as discussed above. The O(1)-C(1) and O(2)-C(2) distances for the bridge site adsorption are 1.34 Å and 1.38 Å, respectively. The calculated effective charges of two oxygen are O(1)-0.91 and O(2)-0.85, respectively. The electrons are mainly transferred from the neighboring carbon atoms. The effective charge of nitrogen did not change. The chemisorption energy of the dioxygen on this bridge site is 2.57 eV/atom. The results of the spin-polarized DFT calculations showed that the dioxygen loses its magnetic moment after the chemisorption on the long bridge site to nearly zero. The carbon atoms that adsorb the dioxygen also do not have noticeable magnetic moments, which are different from those of the carbon atoms on the open-edge of CNT with dangling bonds.
The second case of dioxygen adsorption in a short bridge-site on the open-edge of the CNT was also studied. In this case, near the dioxygen adsorption site, the open-edge of the CNT has a noticeable distortion that is attributed to the formation of new C-C bond. The dioxygen adsorption energy at this bridge site is about 0.015 eV/atom, which is much lower than that of the first case of the dioxygen adsorption discussed above. The calculated charges transfer to the two oxygen atom are O(1)-0.03 and O(2)-0.01, which is much smaller than that of the fist case of the bridge-site adsorption discussed earlier. The O-O bond length is about 1.23 Å, which is nearly the same as that in the free O2 molecule state. The O(1)-C(1) and O(2)-C(2) distances for this short bridge site adsorption are 3.65 Å and 3.76 Å, respectively. The dioxygen on this site maintained a magnetic momentum of 2 µB, which is the same as that of the free gas state. All the data clearly show that the dioxygen adsorption on this short bridge-site is a physisorption, which is attributed to the loss of the dangling bonds of the carbon atoms on the open-edge of the CNT. To save space, this typical physisorption site graph is omitted here.
From Figure 2 (c) and the Bader charge analysis in Table I, it can be seen that the oxygen-carbon form a weak chemical bond (sp hybridization) at long bridge site, since the two carbon atoms at long bridge site would otherwise possess some dangling bonds, while there is strong C-C bond at the short bridge site on the open-edge of the short CNT, thus the interaction between O2 and carbon atoms at short bridge site is weak physical adsorption.
|O2 Site||O-C (Å)||O-O (Å)||MM(µB)||Bader charge change|
|Pauling site (End-on C-N complex, Chemisorption)||1.31||1.33||0.027||O(1)-0.99 , O(2)-0.24||1.53|
|Long Bridge site|
|1.34||1.57||0.001||O(1)-0.91 , O(2)-0.85||2.57|
|Short Bridge site (Physisorption)||3.71||1.23||2.0||O(1)-0.03, O(2)-0.01||0.015|
5. Dioxygen dissociation on one nitrogen open edge doped CNT(10, 0)
To study the dioxygen dissociation energy barrier, the NEB method was used to study the minimum energy path from the Pauling site to the dissociated N-O and C-O state. This section will introduce the major results of the simulations and the detailed results will be organized and publish late.27 To save the computation time, only top ring atoms are relaxed since the bottom carbon atoms are relatively stable and the relaxation could be safely ignored. It should be mentioned that dioxygen dissociative adsorption minimum energy path from the NEB method is not necessary the only possible reaction path for the ORR or the most important step. However the relative energy barrier comparison can be used as a measure of reactivity.
The calculated initial Pauling site is shown in Figure 3 (a), the stable dissociated N-O(1) and C-O(2) state is shown in Figure 3 (c), while the transition state is shown in Figure 3 (b). The calculated minimum energy barrier is 0.55 eV which is reduced tremendously than the gas phase dissociation energy 5.1 eV. In the transition state (Figure 3(b)), the dioxygen bond length is elongated and the dangling O(1) is moved toward the nitrogen atom. The top carbon ring atoms have a slight re-adjustment to balance the residue forces. It should be emphasized here that in real PEMFC application, the Pauling site dioxygen may be dissociated through other channels even before it is dissociated into N-O(1) and C-O(2) considering the existence of the hydrogen atoms and OH. At the transition state, there is significant charge redistribution on nitrogen, O(1)-O(2), and bonded carbon atoms. It is the local electric field that weakens the O(1)-O(2) bond and dramatically decreases the transition barrier energy.
In this review report, first principles spin polarized DFT simulations of nitrogen substitutionally doped (10, 0) carbon nanotube (CNT) for dioxygen adsorption and dissociation are performed. The calculated results show that nitrogen prefers to stay at the open-edge of the short CNTs. Two O2 chemisorption sites are found, the carbon-nitrogen complex (Pauling site) and carbon-carbon long bridge (long bridge) sites. Dioxygen O2 can be chemisorbed on and reduced on the carbon-nitrogen complex at the open-edge of the CNT and on the open carbon-carbon sites. The carbon atoms on the open-edge of the short CNT can possess a magnetic moment of 0.59 μB/atom, which is due to the existence of the dangling bonds of these C atoms. The chemisorption of dioxygen O2 onto both Pauling site and long bridge sites at the open-edge of the short CNTs will reduce the magnetic moments of the carbon atoms to nearly zero. Further spin polarized NEB method minimum energy barrier simulations show that the Pauling site is the possible O2 dissociation site with a reaction barrier 0.55 eV. The unique open-edge structure and charge redistribution are crucial to the novel doped CNT catalyst design.
The work is funded in part by the NSF-LASiGMA program (grant number #EPS-1003897), the NASA/LEQSF (2009-2012)-Phase3-03, the LaSPACE/NASA (grant number NNG05GH22H), DOE award No. DE-FE0004734 and DE-FE0003693, and the LONI institute. GLZ acknowledged the funding supports of the National Science Foundation (Award No CBET-0754821) and the Air Force Office of Scientific Research (Award No FA9550-09-1-0367).