Abstract

Graphite and silicon carbide (SiC) are important materials of fuel elements in High Temperature Reactor-Pebble-bed Modules (HTR-PM) and it is essential to analyze the source term about the radioactive products adsorbed on graphite and SiC surface in HTR-PM. In this article, the adsorption behaviors of activation product Cobalt (Co) on graphite and SiC surface have been studied with the first-principle calculation, including the adsorption energy, charge density difference, density of states, and adsorption ratios. It shows that the adsorption behaviors of Co on graphite and SiC both belong to chemisorption, with an adsorption energy 2.971 eV located at the Hollow site and 6.677 eV located at the hcp-Hollow site, respectively. Combining the charge density difference and density of states, it indicates that the interaction of Co-SiC system is stronger than Co-graphite system. Furthermore, the variation of adsorption ratios of Co on different substrate is obtained by a model of grand canonical ensemble, and it is found that when the temperature is close to 650 K and 1700 K for graphite surface and SiC surface, respectively, the Co adatom on the substrate will desorb dramatically. These results show that SiC layer in fuel element could obstruct the diffusion of Co effectively in normal and accidental operation conditions, but the graphite may become a carrier of Co radioactivity nuclide in the primary circuit of HTR-PM.

1. Introduction

Cobalt-60 (60Co) is a kind of long half-life γ-ray radionuclide and it could be generated through activation reaction of impurities (59Co and 60Ni) in the material of fuel elements and metal/nonmetal reactor internals of High Temperature Reactor-Pebble-bed Modules (HTR-PM). The behavior of 60Co is important for the safety analysis and radiation protection design in HTR-PM and has attracted a lot of attentions. In the fuel element of HTR-PM, 60Co in porous pyrolytic carbon (buffer layer) and dense inner pyrolytic carbon (IPyC) of TRISO fuel particle will diffuse and interact with silicon carbide (SiC) of TRISO fuel particle ‎[1], which is considered to possibly influence the performance of SiC. In the primary circuit of HTR-PM, 60Co on the surfaces of fuel elements, graphite reflectors, graphite reactor internals, and metal reactor internals will adsorb on the graphite dust, which is generated through abrasion or corrosion effect when the fuel elements flow in the primary circuit and is playing a significant role in contributing to the source term of HTR-PM ‎[2]. It is reported that, in the end of lifetime of HTR-PM, the radioactivity of 60Co in the primary circuit will be accumulated to 8.6 1010 Bq ‎[3], while the activity concentration of other radionuclides is at least two orders of magnitude lower than that of 60Co.

It is known that the study of interaction between Co and reactor material is essential but it is difficult to obtain the experimental result. Fortunately, the first-principle calculation, a powerful tool to study on an atomic scale, provides a way to research the above issue ‎[47]. Based on the density functional theory (DFT) ‎[8], the first-principle calculation is implemented in the Vienna Ab initio Simulation Package (VSAP) by the group of Kresse et al. ‎[911], which could be used to study the atom-material interaction in microlevel ‎[1214]. A number of theoretical studies about the adsorption behaviors using first-principle calculations, especially for the graphite ‎[1519] and SiC nanotubes ‎[2027], have been published. Ancilotto and Toigo have performed first-principles total-energy calculations to study the adsorption of potassium on graphite ‎[15]. Electronic structure calculations based on spin-polarized DFT with the generalized gradient approximation (GGA) and ultrasoft pseudopotentials are used to investigate the interaction between H atoms and a graphite surface ‎[16]. The single Co atom adsorbed on some graphite materials also is discussed ‎[2831]. Wehling et al. have researched the orbitally controlled Kondo effect of Co adatoms on graphene ‎[28]. Rudenko et al. have researched the adsorption of Co on graphene and analysis of the electron correlation effects from a quantum chemical perspective ‎[29].

However, there is few works on the adsorption/desorption of Co on graphite/SiC and the mechanism of interaction is also not clear. In this work, the behavior of Co adsorbed on graphite and SiC surface with DFT will be studied, including the charge density difference (CDD) and the density of states (DOS). Furthermore, the mechanism of adsorption will also be discussed. At last, the variation of adsorption ratios of Co will be given by a model of grand canonical ensemble, which is significant for understanding the adsorption of Co on graphite and SiC macroscopically.

2. Method of Calculation

2.1. The Construction of Graphite and SiC Micromodel

The first-principle calculation is also known as ab initio calculation; it has been performed based on DFT as implemented in the VASP code and employed the projector augmented wave (PAW) pseudo-potential ‎[32] and GGA-PBE exchange-correlation functional [33] to describe the interaction of electron-ion. Since a large number of basis functions are usually required to describe the electronic wave functions appropriately, it is very demanding to use the first-principles pseudo-potential method to calculate the carbon properties ‎[15]. To simplify the calculations, we choose a single, isolated graphite monolayer with 6 6 primitive cells (Figure 1(a)), by using a -point mesh of 1 1 1, as the substrate to calculate the adsorption energy. According to the anisotropic character of bonding in graphite, the coupling between adjacent graphite layers is much weaker than the in-plane coupling between carbon atoms ‎[14, 34, 35]. This means that the electronic properties of a monolayer can usually provide a reasonable picture of the electronic properties of the infinite crystal ‎[36]. As shown in Figure 1(b), graphite is here represented by graphene; there are three different adsorption sites on the graphene including Top (T), Bridge (B), and Hollow (H). SiC used in TRISO fuel particle of fuel element in HTR-PM is -SiC and the corresponding structure is diamond cubic crystal. Different from pure diamond, there are four silicon atoms surrounding carbon in each SiC crystal. Therefore, the interaction between adjacent atoms in a crystal is very strong and makes the structure of SiC especially stable, which is conducive to SiC to block the activation products diffuse from the fuel particle. The supersurface of SiC(001) lattice with 3 3 1 primitive cells (Figure 2(a)) and the -point mesh as 2 2 1 are employed to calculate the adsorption energy. The four kinds of adsorptive sites are including Top (T), Bridge (B), hcp-Hollow (hH), and fcc-Hollow (fH) as shown in Figure 2(b).

2.2. Adsorption Energy

The adsorption energy of Co is defined aswhere , , , and are the energies of the adsorption, substrate, a single Co atom, and the Co-substrate system, respectively. , , and could be obtained by the first-principle calculations directly, and then could be obtained from (1). There are two influencing factors that need to be considered which may affect the value of . One is the adsorptive sites and the other is the height between the adatom and substrate. The adsorption energies at different heights of different adsorptive sites are calculated firstly and after obtaining the most suitable position which located at the height corresponding to the largest adsorption energy, the total system will be relaxed to obtain the value of and obtain the optimum adsorptive site and height. It should be mentioned that the “-” sign of represents the direction of energy change, and the adsorption energy is only related to the absolute value of . Meanwhile, the CDD and the DOS could be obtained and are employed to analysis. The CDD is defined as the difference between the total charge density and the atomic charge densities ‎[20]. It shows the visual image of electron transfer between nuclide and substrate directly. As suggested by Khantha et al. ‎[37], the redistribution of the charge density indicates that the interaction potential between the adatom and substrate is composed of two distinct contributions: a screened Coulomb interaction resulting from the charge transfer between the adatom and substrate atoms and a van der Waals type interaction ‎[38, 39]. The equilibrium position of an adatom above graphite/SiC is largely determined by the charge redistribution that takes place in this short-separation region. DOS shows the distributions of electrons in each energy level, which are helpful to understand the behaviors of atom absorbed/desorbed on material in microlevel. Fermi energy level could also be obtained in the abovementioned process. The relative position of the Fermi level in the DOS shows the occupation of the states and the features of bonding ‎[39, 40]. If atom is absorbed on the surface of substrate, the peaks of the state for both adatom and substrate compared to will be shifted. Supposing that the DOS is changed, the adsorption behavior between the adatom and substrate is chemisorption. Otherwise, it is physisorption. If the energy state of atom has a left shift, this means that the atom acquires electrons. Then its energy decreases, so that the peak width of state will become wider and more occupied. If the energy state shifted towards right, the atom will lose electrons.

What should be emphasized here is , CDD, DOS, and are obtained in the condition of 0 K with DFT method, but the temperature condition of the fuel element and primary circuit of HTR-PM is below 1000 K under normal operation and below 1850 K under accident operation. It is necessary to point out that the energy of electron is not sensitive to the temperature at the range of T ~103 K [41], so the results obtained by DFT can be applied to study the interaction between Co and graphite and SiC of HTR-PM.

2.3. Adsorption Ratio

The adsorbed nuclide atoms on substrate and separate nuclide atoms will form a dynamic equilibrium system. According to literatures ‎[14, 42], a model of grand canonical ensemble is considered to describe this system and calculate the balanced adsorption ratio of Co on the graphite and SiC. In this model, the substrate materials adsorbing Co atoms are considered as an open system, and the separate Co atoms in circumstance (i.e., the coolant in the primary circuit of HTR-PM and the gas interlayer between SiC and IPyC of TRISO fuel particle) are treated as particle and gas source. On the theory of the statistic physics, the relation between the number of adsorbed Co atoms and the total adsorptive centers is given by where is the chemical potential of Co atom, is Boltzmann’s constant, and is the temperature of the circumstance. So the adsorption ratio can be defined as

Because the concentration of radioactive products is very low ‎[42], the ideal gas model can be employed to deal with nuclide gas source directly. If an equilibrium state is achieved in the reactor, the chemical potentials of nuclide atoms in the open system will be identical to those of gas source. Therefore, according to the ideal gas model, the chemical potential of Co can be written as follows:where is the density of the nuclide in the circumstance, is the atomic mass of Co, and is the Planck constant.

Combined with (2), (3), and (4), the adsorption ratio is accessible to obtain asIt is seen from (5) that only three parameters of , , and are undetermined. The numerical value of and of circumvents in HTR-PM can be determined by either calculation or experimental measurement and can be obtained from (1).

3. Results and Discussion

3.1. Adsorption Energy

In order to obtain the most stable adsorptive site, the adsorption energy with different adsorptive sites has been calculated. For each adsorptive site, the adsorption energy depends on the height between the adatom and substrate and there exists a maximum value which could be considered as the exact adsorption energy of corresponding site. Figure 3 shows the variations of adsorption energies of Co on graphite and SiC at different sites with different height.

The adsorption energies of Co on graphite and SiC at different adsorptive sites are listed in Tables 1 and 2. To different substrate, the adsorption energies and the most stable adsorptive sites of Co are different. It shows that the most stable sites for Co-graphite and Co-SiC system are Hollow and hcp-Hollow with the adsorption energies 2.971 eV and 6.677 eV, respectively. These results indicate that the interaction of Co-SiC system is stronger than Co-graphite system.

3.2. Analysis of Electronic Structure

To identify the adsorption behavior of the Co adatom on different substrates, the electronic structure of the two systems should be analyzed. Figure 4 shows the optimized stable configurations and their corresponding CDD for the Co-graphite system and Co-SiC system. The blue areas represent the decreases of charge density and the yellow areas represent the increase of charge density. Even though there is no obvious electron transfer between Co and graphite in Figure 4(a), the adsorption behavior of the Co-graphite system could be considered as chemisorption due to its high adsorption energy in Table 1, while in Figure 4(b) there are a lot of yellow areas that appeared between Co and SiC surface, which suggests that a strong chemical bond is formed in this area. This result indicates that the Co has a stronger interaction with SiC than graphite. The result also coincides with the abovementioned adsorption energies of these two systems.

In order to understand the deep features of interactions between the adatom and substrates, the DOS for single Co atom, the substrate of graphite and SiC, and the Co-substrate system are calculated. It should be mentioned that the -point meshes to integrate the Brillouin zone have been changed to 2 2 2 and 4 4 1 for the DOS calculation of Co-graphite system and Co-SiC system, respectively. The results of DOS are shown in Figure 5. When the Co adsorbed on graphite and SiC surface, the positions of Co adatom are different. In the following discussion, different systems and different positions corresponding to the adsorption state are used in the calculation of isolated Co atom for Co-graphite and Co-SiC, respectively, which makes the DOS of Co have some differences in these two systems.

Figures 5(a) and 5(b) give a comparison of DOS for the Co-graphite system, the isolated graphite substrate, the isolated Co atom, and the adsorbed Co adatom. For isolated Co atom, the peaks of the -state nearly locate at and that of -state locate at the right of . When the Co atom is adsorbed on the graphite surface, the peaks of the -state have a left shift towards a lower energy than and have evident broadenings owing to adsorption interaction. It is also found that the peak of the -state shifts slightly right and the peak intensity becomes too weak to be identified when the Co atom is adsorbed on the graphite surface. Combining the above results and total DOS of graphite, it is feasible to know that the interaction between Co and graphite is very strong and there are electrons transferring on the interactive interface, which is the direct evidence of chemisorption. The situation of the Co-SiC system, as shown in Figures 5(c) and 5(d), is a little different. When the Co atom is adsorbed on the SiC surface, the peak intensities of the -state become too weak to recognize and there are only shifts and broadenings of the peaks of -state. Remarkably, although most of the peaks of -state move to the left, there is still right shift, which indicates that Co adatom not only gains partial electrons from SiC but also loses electrons to SiC. These results could be explained from the view of electrochemistry that the electronegativity of Co is close to Si and there exists resonance between the -state of Co adatom and SiC surface when the partial electron transferring turns out. That is to say, due to the orbital hybridization, the adsorptive behavior between Co and SiC is also chemisorption. Based on the electronic analysis, it presented that the interaction between Co and SiC is stronger than that between Co and graphite significantly.

3.3. Adsorption Ratio

The adsorption ratio, as shown in (4), is determined by three parameters , , and . Among them, can be obtained from (1) by the first-principle calculation, the temperature is variational, and the density of Co () in circumstance could be obtained by experiment or theoretical calculations. Considering that the adsorption conditions of 60Co on graphite and SiC surface are different in HTR-PM, the density of 60Co used for calculation is not the same. Based on the experimental data, the density of 60Co in the primary circuit could be set to be  m−3 [3]. The density of 60Co on the surface of the SiC was calculated according to the content of 59Co impurity in the graphite matrix, and the calculated result was  m−3. In order to investigate the trend of adsorption ratio varying with temperature T, the adsorption ratios of Co adatom on graphite and SiC surface calculated by (4) are shown in Figure 6. It can be easily noted that the adsorption ratio of the Co decreases exponentially with the increase of temperature.

According to Figure 6, there exists an inflexion point while the adsorption ratio varies with the temperature. When the temperature is close to 650 K and 1700 K for graphite surface and SiC surface, respectively, the adsorption ratio of Co decreases dramatically. The temperature at the inflexion point is suggested to be defined as the critical temperature, . If the temperature is lower than , most or even all of Co will be absorbed on the substrate surface. If the temperature is higher than , it stands for the idea that the Co on the substrate will have a lot of transitions from absorption to desorption.

In HTR-PM, the average temperature of the fuel element and the primary circuit under normal operation is usually less than 1000 K, and the maximum design value of temperature under accident conditions is about 1850 K. From these data, it shows that Co could be adsorbed well on the surface of SiC in the normal operating temperature range and only at the highest temperature in supposed accidents, Co will be desorbed and diffused from the SiC surface. For graphite, Co is susceptible to a large amount of desorption and diffusion from the graphite surface whether in normal operation or accident conditions. As a result, the Co absorbed on SiC surface performed a better adsorption behavior than on graphite surface. Therefore, SiC could block the diffusion of Co, while graphite could not if the temperature is higher than 650 K.

4. Conclusion

In this paper, the first-principles DFT theory was used to study the adsorption behavior of Co on graphite and SiC surface in HTR-PM. The adsorption energies corresponding to different adsorptive site and height were calculated, and the adsorption energy for Co on graphite and SiC surface is 2.971 eV and 6.677 eV, respectively. In order to know the underlying physical nature of the adsorption behavior of Co, we have analyzed the electronic structure of the adsorption system by CDD and DOS. The results suggested that the adsorptive behavior of the two systems belongs to chemisorption. All of these results could be used to explain the result of adsorption energy from different aspects; the source term from the microscopic mechanism has been analyzed; it helps us to understand the adsorption behavior of the Co more clearly. In addition, combining the adsorption energy obtained by the first-principle calculation, the adsorption ratio was derived with different temperature and the density of Co in HTR-PM. The critical temperature for graphite surface is 650 K, and for SiC surface it is 1700 K.

In a conclusion, the adsorption performance of Co on the SiC surface was better than that on the graphite surface. In the primary circuit of HTR-PM, the Co adatom seems to be easily desorbed from the graphite surface. Compared with graphite, SiC could suppress the diffusion of Co effectively for both the normal operation and the accident conditions in the fuel elements. The results of adsorption ratios can also be used for nuclear safety evaluation of Co release in the normal operation and the supposed accident of the HTR-PM.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China (Grant no. ZX06901).