Juan Liua,
Tongxiang Liang*b,
Chen Wang*a and
Wenzheng Laic
aState Key Laboratory of New Ceramics and Fine Processing, Institute of Nuclear and New Energy Technology, Tsinghua University, Beijing, 100084, China. E-mail: 9120160073@jxust.edu.cn
bSchool of Materials Science and Engineering, Jiangxi University of Science and Technology, Ganzhou, 341000, China
cSchool of Materials Science and Engineering, Tsinghua University, Beijing, 100084, China
First published on 13th January 2017
Oxygen adsorption and desorption of gasification products are two factors that influence graphite oxidation behavior. Whether the catalytic or inhibiting effect occurs after B and N doping still remains controversial. In this paper, the activation energies for the adsorption of oxygen and desorption of CO on a modified nuclear graphite surface, combining vacancy and dopants, are studied by density functional theory (DFT). Vacancy defected graphite is more sensitive to O2 adsorption because of the unpaired electron on the dangling carbon atom. However, substitution with B or N makes the adsorption more difficult, since substituting C by a B or N atom eliminates the effect of the unpaired electron, which means that these dopants can be applied to improve oxidation resistance of nuclear graphite or other carbon materials. Desorption calculation of adsorbed radical (C–O) indicates that substitutional B or N atoms in graphite facilitate CO desorption, in this respect B or N dopants play a role of catalyst for graphite oxidation. Introduction of dopants such as B, N results in weaker O2 adsorption capability on vacancy defected graphite while they help desorption of oxidation products such as CO. Both the inhibiting and catalytic effects of these dopants are an example of the compensation effect.
Defect investigations in graphite and related materials have attracted considerable attention by the desire to investigate irradiation damage taking place in a HTR.1 Vacancy defects are proved to be more prevalent in irradiated materials and are considered as the primary defects in irradiated graphite or nanotubes.2 The effect of neutron irradiation in nuclear graphite leads to the change of volume, mechanical and thermal physical properties, meanwhile to the happening of radiolytic oxidation, all these changes are associated with the defect of graphite. Reactive species such as oxygen can easily react with defective graphite surfaces thus resulting in enhanced carbon erosion yields. DFT study shows that adsorption energy of O2 on monovacancy defected graphite is approximately 10 times stronger than that on defect-free graphite surface.3 Theoretical and experimental investigations have proved that O2 tends to physically adsorb on perfect graphite surface.4,5 It is the extensive conjugate π bonding among carbon atoms that makes the graphite sheet chemically inert. However, the displacement of C atoms in graphite layers by neutron irradiation causes electronic redistribution and behaves as active sites that promote dissociation, transformation and activation of adsorbed gases.6,7 Besides oxygen, other gas molecule such as CO2 on monovacancy surface is also found to be more active than that on pure graphene layer.8 These studies suggest that graphite plane with vacancy defect is active enough to interact with reaction species.
Based on the reactivity of vacancy defected graphite, it is of significant importance to find a way to reduce oxidation. The use of boron as doping element for carbon-based materials has been proved to improve oxidation resistance. Previous studies have found that oxidation rate of graphite decreases 3–6 times owing to boron substitution.9,10 Boron is one electron less compared with carbon, and the substitution is expected to modify electronic and physicochemical properties.11,12 Boron substitution has also been applied in Li ion batteries to improve carbon's capacity.13 Furthermore, boron has also been introduced in carbon nanotubes and diamond electrodes to modify electronic structure as well as electrochemical performance.14,15 Other impurity atoms (such as N, S, Al) are also applied into graphitic materials to affect sensitivity for gas adsorption.18 Peng et al. has investigated the improvement of sensor devices by substitutional doping of impurity N atoms into pure single-walled carbon nanotubes (SWCNTs). It is found that CNx nanotube is more efficient to detect toxic species than pristine CNTs probably because of the existence of highly reactive pyridine-like composition in tube walls.16,17
Improving oxidation resistance of nuclear graphite is critical to security operation of HTR reactor. Even though the observation of catalytic oxidation19,20 of boron at low concentration makes the role of substitutional boron more mysterious, we expect to understand the mechanism of B or N doping on the behavior of O2 adsorption and CO desorption. In the present study, the effect of substitutional boron and nitrogen on the structure and electronic states of vacancy defected graphite surface is investigated by using DFT method. Adsorption of O2 molecule on impurity doped systems is compared with that on un-doped one. Meanwhile, desorption of adsorbed radical for CO formation is also studied by nudged elastic band method.
In our calculation basal plane is chosen because of its least reactive property than zigzag and armchair surfaces. A 4 × 2 × 1 supercell containing 32 atoms (without vacancy) is set within a single graphite layer. Periodic boundary condition is set in two spatial directions. The vacancy defected graphite is created by removing one carbon atom and then fully optimized. One layer model is used by the reason that van der Waals interactions between layers are relatively weaker than the planar bonding, and consequently the impact on the surface processes is not important.25,26 To further confirm the accuracy of one layer model, we have compared the vacancy formation energy, chemisorption energy of O2 and desorption energy barrier of CO on one layer and four layer vacancy defected graphite surfaces. The details are listed in Table 1. The energy values are almost the same for these two models. Therefore one layer model is selected in this work by computational resource limitation. The optimized lattice constants of the planar structure are a = b = 2.46 Å, that is agree with theoretical and experimental reports.27,28 Our supercell is working within a box in which three lattice vectors are perpendicular to each other. The graphite layers are separated to have a vacuum of ∼16 Å to avoid the influence of neighboring slabs although the fact that an 11 Å vacuum layer is large enough to reduce the periodic neighbor interaction.29,30 The same supercells are applied for the systems with O2 interaction as for the isolated graphite species. The k-point is set to 5 × 5 × 1 for structure optimization, which increases to 12 × 12 × 1 for more accurate electronic calculation including density of states (DOS) and charge density difference (CDD). All atoms in the supercell are relaxed to reach the equilibrium configurations, and the calculations of reaction energy barriers are performed by employing the nudged elastic band (NEB) method.
Type | Ef_v (eV) | Eads_O2 (eV) | Edes_CO (eV) |
---|---|---|---|
1 layer | 7.76 | −3.85 | 5.94 |
4 layer | 7.58 | −3.82 | 6.01 |
The adsorption energy during the interaction between oxygen molecule and graphite is defined by following formula:
Eads = Egra+O2 − Egra − EO2 | (1) |
The charge density difference (CDD) is defined as:
Δρ = ρ(graphite + O2) − ρ(graphite) − ρ(O2) | (2) |
Fig. 1 Relaxed structures of graphite with or without dopants. (a) Graphite with vacancy defect (VG), (b) boron doped and (c) nitrogen doped vacancy-defected graphite. |
Fig. 2 Density of state (DOS) analysis for (a) vacancy-defected graphite, (b) B-doped vacancy-defected graphite and (c) N-doped vacancy-defected graphite. |
For dopant defected graphite, a single impurity atom is inserted in the defect region. Three doping sites (C1, C2 and C4) are considered after taking into account the symmetry of the vacancy. It is found that C1 substitution exhibits lower energy and higher stability in comparison with C2 and C4 sites for B, N dopants. Therefore substitutions at C1 site are mainly considered in further study upon O2 adsorption. Fig. 1 shows the structures of B- and N-doped vacancy-defected graphite after optimization, named as BVG and NVG. For BVG, substitution of B causes small distortion of C2 and C3 atom and makes them move out of the original plane. The substitution in this case causes a larger B–C distance, which elongates from 1.37 Å in VG to 1.61 Å. Meanwhile, the distance between C2 and C3 is elongated from 2.11 Å in VG to 2.55 Å, while the distance between B–C2 is changed from 2.56 Å in VG to 2.16 Å. The structure modification is mainly contributed by the introduction of C–B bond. As we know, the covalent radius of boron (0.88 Å) is slightly larger than that of carbon (0.77 Å).33 A C–B bond cannot replace a C–C bond in equal length thus resulting in structure modification. In order to have a deeper insight into the effect of dopants in vacancy graphite (VG), the obtained DOS of impurity-doped graphite are also illustrated in Fig. 2. It is learned that the defected graphite types are metallic with a zero band gap which is different from perfect graphite.34 Substitution of B leads to more symmetrical DOS distribution near the Fermi level as performed in Fig. 2(b). This result is different from Zhou's work which gives the DOS for BVG with asymmetric disposition.35 In our opinion, the symmetrical DOS of BVG is more convincing. As we know, the existence of unpaired electron in VG causes asymmetrical spin-up and spin-down DOS. The valance electron configuration of B is 2s22p1, one electron less than C. Substituting a C atom located at the vacancy site by a B atom would eliminate the effect of the unpaired electron thus leading to symmetrical DOS. Faccio's work has also proved the symmetric DOS distribution by substitution of the dangling C atom.32 Large DOS overlap of C2 and C3 demonstrates the interaction of two carbon atoms. The BVG exhibited a slight magnetic moment of −0.002 μB.
In the case of N-doped VG (NVG), the bond length of N–C is shortened to 1.33 Å along with the distance between C2 and C3 (1.87 Å), which is shorter than that of VG (2.11 Å). The substitution of N makes all atoms locate in the graphene plane and no out-of-plane atomic displacement is discovered. According to Fig. 2(c), it shows that the DOS distribution for spin-up and spin-down are completely symmetric. Based on the valence bond theory, two N-2s electrons are pairing and the remaining three electrons named 2px, 2py and 2pz are unpaired. The newly formed C–C bond is attributed to the interaction of unpaired electrons in C2 and C3. Two unpaired electrons of N atom interact with neighbouring C atoms and the rest unpaired electron contributes to the π electron system, which results in the NVG presented non-magnetic state.36 In general, the dopants mostly affect the valence electrons around Fermi energy. Considering the symmetrical or asymmetrical distribution of DOS peaks near Fermi energy, we suspect the VG with B, N impurity are more inert than VG itself. The relationship between symmetrical property of DOS distribution and reactivity of graphite will be studied in detail in our next paper. In the present work, we anticipate that the introduction of dopants will improve oxidation resistance of graphite. In following sections, interaction of molecular oxygen with VG and dopant induced VG will be analysed.
Fig. 4 CDD for adsorption structures of O2 on graphite (a) NVG-A and (b) NVG-B. Orange area indicates a positive value while blue colour indicates the negative value of electron transfer. |
For a better understanding of the mechanism, DOS of the dopants before and after O2 interaction is calculated. The most stable adsorption configuration in each graphitic type (VG-C, BVG-C or NVG-C) is chosen for DOS calculation. As performed in Fig. 5, local density of states (LDOS) of the dangling atoms (C1, B and N) before and after O2 adsorption is depicted. It is found that the intensive peaks in valance states near the Fermi level largely decrease after bonding with O atom. Especially for C1 in VG (Fig. 5(a)), the adsorption produces the disappearance of electronic spin-up states above the Fermi energy as well as the spin-down states below the Fermi level. The change of electronic properties is due to the changing of distribution of unpaired electrons by orbital interaction. The charge distribution of systems is also examined by using Bader charge analysis. As summarized in Table 2, O2 molecule behaves as an accepter capturing charge from graphite. This is due to the stronger ability of O atom to grab charge, which is also observed by other studies.25,37 The charge transfer between oxygen and graphite surface is much larger in VG-C (−0.96e) than that of BVG-C (−0.78e) and NVG-C (−0.38e). The results indicate that the presence of dopants such as B or N inhibits the interaction between O2 molecule and vacancy defected graphite surface. Fig. 5(b) and (c) plot the LDOS of dopants (B, N) with and without adsorption of O2. It shows that interaction with B or N leads to unsymmetrical spin-up and spin-down state due to the orbital interaction changing the original structures of paired electrons. Meanwhile, spin polarization is introduced to the system during adsorption on BVG or NVG, with a magnetic moment of 0.84 μB for BVG-C and 1.09 μB for NVG-C. Furthermore, partial electronic density of states (LDOS) of the dangling atoms in these three graphite systems is analysed (Fig. 6). Compared to the un-doped systems, there is rather larger overlap between the doped atom and C3 which implies the existence of strong hybridization between the dopant and the C atoms around the vacancy. This is in accordance with Qin's results that the combination of defect-dopant affects the adsorption mechanism to some extent due to the orbital hybridization between C atoms and the dopants.38
Type | VG-C | BVG-C | NVG-C |
Ead (eV) | −3.85 | −2.57 | −0.13 |
d(Å) | C–O: 1.23 | B–O: 1.39 | N–O: 1.29 |
Q(e) | C: +0.91 | B: +0.33 | N: +0.84 |
O: −0.96 | O: −0.78 | O: −0.38 |
Fig. 6 Partial electronic density of states (LDOS) of the dangling atoms in (a) VG-C, (b) BVG-C and (c) NVG-C after O2 adsorption. |
As shown in Fig. 7, desorption paths for CO on VG and BVG are displayed. The diffusion path on NVG is similar to that on BVG and thus will not be presented here. The initial situations are those with the adsorbed oxygen atom bonding to the dangling C atom with C–O length of 1.23 Å, and the dopant is located at the position neighbor to the dangling carbon. The substitution of B leads to much longer C–C (1.50 Å) and B–C (1.56 Å) bond compared to the corresponding C–C (1.48 Å) bond in VG system. To some extent, the substitution weakens C–C bond. NEB calculation shows that desorption of CO takes place step by step on these graphitic surfaces. The C–O radical tends to break a C–C bond in the first step. During the second stage, the other C–C (C–B or C–N) bond is elongated for preparation of bond breaking. As the process continues, the C–O radical move off the surface and ultimately a CO molecule desorbs. The formation of CO on doped graphite surface takes place with much lower activation energies. For instance, the energy barrier for VG is 5.94 eV, which is 2.63 eV higher than BVG (3.31 eV) and 2.32 eV higher than NVG (3.62 eV). This implies that substitutional B or N facilitates CO formation and finally results in the elimination of carbon in graphite.
In fact, the lowering of activation energy barrier for oxidation of B-doped graphite has been observed by Allardice et al.20 Thus substitutional boron is considered to anticipate the catalytic effect even though this effect is masked by inhibitory action of B2O3 in their study. Based on experimental and theoretical evidence, Radovic provides an effect that the redistribution of π electrons in the presence of substitutional boron results not in weaker C–O bonds and stronger C–C bonds,11 but in exactly the opposite effect that as π electron acceptor, boron is to both reduce the π electron density and accept additional electron density into its vacant 2pz orbital. By accepting electrons from the graphene layer, the substitution ultimately results in the weakening of C–C bonds along with strengthening of C–O bonds.19
Our investigation shows that under conditions favouring desorption control boron appears to be a catalyst, while under conditions favouring adsorption control boron behaves as an inhibitor. This phenomenon indicates another example of the compensation effect, which is interpreted as a “macroscopic complementarity principle” by Radovic.19
In summary, as shown in Fig. 8, the effect of B or N doping on graphite oxidation is two-fold. Dopants of B and N induce a weaker O2 adsorption capability and they behave as an inhibitor to graphite oxidation. On the other hand, the removal of CO gasification product becomes easier after B or N doping. At this point, dopants behave as catalyst to graphite oxidation.
This journal is © The Royal Society of Chemistry 2017 |