# Diffusion Behaviors of Hydrogen Isotopes in Incoloy 800H: A First-Principles Study

Authors: Hongyu Chen; Hong Li; Chuan Li; Jianzhu Cao; Chao Fang
Journal: Science and Technology of Nuclear Installations
 (2017)

---

## Abstract

Incoloy 800H is one of the main stainless steel materials used in steam generators with High Temperature Reactor Pebble-bed Modules (HTR-PM). In this study, the diffusion behaviors of hydrogen isotopes in Incoloy 800H were investigated with first-principle calculations. Numerical results reveal that the starting and ending positions of the diffusion process are the two adjacent and most stable octahedral sites surrounded by Fe atoms and Ni atoms, and the diffusion follows an indirect path via the metastable tetrahedral sites and octahedral sites surrounded by Fe atoms and Cr atoms. The diffusion activation energies of hydrogen (H), deuterium (D), and tritium (T) in Incoloy 800H are investigated by first-principles calculations with the same approximate value ofQ=0.757 eV; the diffusion coefficient frequency factors are also obtained with values of D0=1.56×10-6, 1.10×10-6, and 8.99×10-7 (m2/s) for H, D, and T, respectively. Furthermore, the theoretical results are compared with the experimental data, and it is found that both are in agreement with each other. These results are very helpful for understanding the diffusion behaviors of hydrogen isotopes in Incoloy 800H and can be used to guide the tritium source term analysis of secondary circuits in HTR-PM, which are first studied from a microperspective.


---

## Body

## 1. Introduction

Tritium is one of the most important radionuclides in reactors owing to its very harmfulβ-radiation and long half-life. Although the β-particle produced by the decay of tritium is easily absorbed and blocked by the skin, tritium oxide (HTO, DTO, or T2O) can be inhaled and can combine with organic matter, resulting in internal irradiation, which is very harmful to the human body. In the primary circuit of High Temperature Reactor Pebble-bed Modules (HTR-PM), tritium could be generated in two ways: ternary fission and neutron reactions [1]. Additionally, the experimental results have indicated that tritium diffusion is a serious issue for High Temperature Gas-cooled Reactors (HTGR) because tritium could easily traverse metal and permeate from the primary circuit into the secondary circuit under high temperature conditions [2]. It is reported that, at the end of the 20-year lifetime of a 10 MW High Temperature Gas-cooled Test Reactor (HTR-10), the radioactivity of 3H in the primary circuit will have accumulated to 5.7 × 109 Bq, and it is estimated that the concentration of 3H in the secondary circuit of a 1000 MWe HTGR could be 1.85 × 105 Bq/L in the equilibrium state [1, 3]. The steam generator pipes are the main boundary for tritium permeation, and Incoloy 800H is a major material of those pipes in HTR-PM, which makes it essential to understand the diffusion mechanism of tritium in Incoloy 800H.In recent years, a number of theoretical studies based on first-principles calculations of diffusion mechanisms and behaviors of gas atoms in metal have been published [4–9]. Wimmer et al.have predicted the temperature-dependent diffusion coefficients of interstitial hydrogen, deuterium, and tritium atoms in nickel [4]. Yang et al. studied the diffusion behaviors of C, N, and O atoms in V metal [5]. Liu et al. have researched the diffusion activation energies of C, N, and O in α-Fe [6]. Apart from these theoretical studies, many experimental measurements on diffusion have been performed, proving that the first-principle calculation is a powerful research tool in materials science [10–14].However, there are few studies on the diffusion of gas atoms in Incoloy 800H and other alloy materials. In this study, the diffusion mechanism and behavior of hydrogen isotopes (H, D, and T) in Incoloy 800H were investigated with first-principle calculations using density functional theory (DFT) [15]. In addition, the numerical results of the calculation, including the diffusion activation energies and the frequency factors of diffusion coefficients, were compared with the experimental results of hydrogen and deuterium obtained by our team. The similarity between these results shows that the study is reliable and sometimes, in nuclear industry, the first-principle calculation could take the place of complicated experiments in rigorous conditions, saving costs and resources.


## 2. Method of Calculation and Experiment

### 2.1. The Construction of Incoloy 800H Model

The calculation was performed based on DFT as implemented in the Vienna Ab initio Simulation Package (VSAP) with the projector augmented wave (PAW) potential method [16–19]. The exchange and correlation effects are described by the generalized gradient approximation (GGA) in the Perdew-Burke-Ernzerhof (PBE) form [20]. To carry out this study, the microstructure of Incoloy 800H has to be determined in the first step.Incoloy 800H is an austenitic alloy with a face-centered cubic structure. Its main components are Fe, Cr, and Ni with the mass fractions of approximately 46%, 32%, and 21%, respectively [21]. The numerical ratio of these three elements can be easily calculated as being approximately 6 : 4 : 3. This ratio is similar to that obtained on replacing two iron atoms by one chromium atom and one nickel atom in each ideal primitive cell of face-centered cubic iron (Figure 1(a)), and, therefore, the ratio can be simplified to 2 : 1 : 1. Finally, the system is constructed using a 2×2×2 supercell containing 16 Fe atoms, 8 Cr atoms, and 8 Ni atoms (Figure 1(b)). To check the convergence of the activation energies, we have also considered a 3×3×3 supercell containing 72 host atoms and one hydrogen atom. The lattice constant was optimized to 3.486 Å, the calculation was performed by using a 5 × 5 × 5k-point mesh, and the plane wave cutoff was 500 eV.Figure 1
(a) The cell of face-centered cubic iron with two Fe atoms replaced by a Cr atom and an Ni atom. (b) 2 × 2 × 2 face-centered cubic supercell. In all figures, green atoms represent Fe, yellow Ni atoms, and orange Cr atoms.
(a)
(b)


### 2.2. Diffusion Activation Energy and Diffusion Coefficient

The most common diffusion mechanism of hydrogen isotopes in crystalline solids is via the vacancies composed of lattice atoms [8]. There are two types of conventional vacancy trap sites in face-centered cubic cells: octahedral sites and tetrahedral sites. In addition, owing to the differences in vicinal atoms, the octahedral sites fall into three types, while there is only one type of tetrahedral site (Figure 2).Figure 2
(a–c) The octahedral site surrounded by different base atoms. (a) The adjacent atoms are Fe and Ni (O-FeNi). (b) The adjacent atoms are Fe and Cr(O-FeCr). (c) The adjacent atoms are Fe, Cr, and Ni (O-FeCrNi). (d) The single tetrahedral site (Te). In all figures, white atoms represent H impurities. To represent the vacancy site more clearly, only the adjacent atoms are shown in these figures.
(a)
(b)
(c)
(d)The general form of the diffusion coefficient is given by(1)D=D0×exp⁡-QkBT,where D0 is the temperature independent frequency factor, kB is the Boltzmann constant, andT is the temperature. Q is the diffusion activation energy, that is, the highest energy barrier for the impurity atoms to diffuse from one most stable vacancy to another [5].According to the transition state theory, the transition rate of atoms in solids is given by [22](2)υ=v·exp⁡-QkBT,where v is the vibrational frequency of the impurity. Therefore, the diffusion coefficient can be expressed as [23](3)D=l2vexp⁡-QkBT,where l is the length of every transition. Using Zener and Wert’s theory, the vibrational frequency can be defined as [24, 25](4)v=Q2ml2,where m is the mass of the impurity. Substituting (3) and (4) into (1), D0 and D can be written as(5)D0=Ql22m,(6)D=Ql22mexp⁡-QkBT.It is clear from the above derivation that the main challenge of the study is to determine the diffusion path and quantify the parametersQ and l.


## 2.1. The Construction of Incoloy 800H Model

The calculation was performed based on DFT as implemented in the Vienna Ab initio Simulation Package (VSAP) with the projector augmented wave (PAW) potential method [16–19]. The exchange and correlation effects are described by the generalized gradient approximation (GGA) in the Perdew-Burke-Ernzerhof (PBE) form [20]. To carry out this study, the microstructure of Incoloy 800H has to be determined in the first step.Incoloy 800H is an austenitic alloy with a face-centered cubic structure. Its main components are Fe, Cr, and Ni with the mass fractions of approximately 46%, 32%, and 21%, respectively [21]. The numerical ratio of these three elements can be easily calculated as being approximately 6 : 4 : 3. This ratio is similar to that obtained on replacing two iron atoms by one chromium atom and one nickel atom in each ideal primitive cell of face-centered cubic iron (Figure 1(a)), and, therefore, the ratio can be simplified to 2 : 1 : 1. Finally, the system is constructed using a 2×2×2 supercell containing 16 Fe atoms, 8 Cr atoms, and 8 Ni atoms (Figure 1(b)). To check the convergence of the activation energies, we have also considered a 3×3×3 supercell containing 72 host atoms and one hydrogen atom. The lattice constant was optimized to 3.486 Å, the calculation was performed by using a 5 × 5 × 5k-point mesh, and the plane wave cutoff was 500 eV.Figure 1
(a) The cell of face-centered cubic iron with two Fe atoms replaced by a Cr atom and an Ni atom. (b) 2 × 2 × 2 face-centered cubic supercell. In all figures, green atoms represent Fe, yellow Ni atoms, and orange Cr atoms.
(a)
(b)


## 2.2. Diffusion Activation Energy and Diffusion Coefficient

The most common diffusion mechanism of hydrogen isotopes in crystalline solids is via the vacancies composed of lattice atoms [8]. There are two types of conventional vacancy trap sites in face-centered cubic cells: octahedral sites and tetrahedral sites. In addition, owing to the differences in vicinal atoms, the octahedral sites fall into three types, while there is only one type of tetrahedral site (Figure 2).Figure 2
(a–c) The octahedral site surrounded by different base atoms. (a) The adjacent atoms are Fe and Ni (O-FeNi). (b) The adjacent atoms are Fe and Cr(O-FeCr). (c) The adjacent atoms are Fe, Cr, and Ni (O-FeCrNi). (d) The single tetrahedral site (Te). In all figures, white atoms represent H impurities. To represent the vacancy site more clearly, only the adjacent atoms are shown in these figures.
(a)
(b)
(c)
(d)The general form of the diffusion coefficient is given by(1)D=D0×exp⁡-QkBT,where D0 is the temperature independent frequency factor, kB is the Boltzmann constant, andT is the temperature. Q is the diffusion activation energy, that is, the highest energy barrier for the impurity atoms to diffuse from one most stable vacancy to another [5].According to the transition state theory, the transition rate of atoms in solids is given by [22](2)υ=v·exp⁡-QkBT,where v is the vibrational frequency of the impurity. Therefore, the diffusion coefficient can be expressed as [23](3)D=l2vexp⁡-QkBT,where l is the length of every transition. Using Zener and Wert’s theory, the vibrational frequency can be defined as [24, 25](4)v=Q2ml2,where m is the mass of the impurity. Substituting (3) and (4) into (1), D0 and D can be written as(5)D0=Ql22m,(6)D=Ql22mexp⁡-QkBT.It is clear from the above derivation that the main challenge of the study is to determine the diffusion path and quantify the parametersQ and l.


## 3. Results and Discussion

### 3.1. Diffusion Activation Energy

To confirm the physical diffusion path, the most stable vacancy site should be found in the first step. The system energyEs of H located in the three octahedral sites and the tetrahedral site was calculated (Table 1).Table 1
System energy of H atom located in the four sites for2×2×2 and 3×3×3 supercell models, respectively.
Vacancy sites
O-FeCr
O-FeNi
O-FeCrNi
Te
2 × 2 × 2
E
s (eV)
−251.937
−251.785
−251.641
−251.321
Δ
E
s (eV)
—
0.152
0.296
0.616
3 × 3 × 3
E
s (eV)
−777.936
−777.908
−777.790
−777.321
Δ
E
s (eV)
—
0.028
0.146
0.615In this calculation with VASP, the projector augmented wave (PAW) pseudopotential was employed to describe the electron-ion interaction while neutron effects and mass effects were not involved. Therefore, the energies of H, D, and T should be the same; this is confirmed by our calculation. This is usually well justified as the barrier for an atomic jump is determined by the electronic interaction, which is identical for two isotopes of the same element, and not by the masses of the nuclei [25].It is found that the most stable site is the O-FeCr site whereEs equals to −251.937 eV, lower than the Te site with a 2×2×2 supercell by 0.616 eV. The 3×3×3 supercell gives the same order of Es of the four sites and an essentially same energy difference of 0.615 eV. Therefore, the impurity atoms must transit through an O-FeNi site or an O-FeCrNi site when they complete a whole diffusion process from one O-FeCr site to another. Owing to the lower Es, the O-FeNi site is considered as the midpoint. As the starting, end, and midpoints were confirmed, there are two possible diffusion paths in the model, as illustrated in Figure 3. In path A, a H atom transits from O1 to O2 in a straight line; in path B, the impurity atom jumps by a broken line, passing by the Te position. To compare the two diffusion paths, the “Transition State Search” method, nudged elastic band calculation [26], was used to search for the transition state.Figure 3
Two possible paths of one impurity atom diffusing from an O-FeCr site O1 to another O-FeCr siteO1′. O2 is an O-FeNi site. Te denotes a tetrahedral interstitial site.As shown in Figure4, the highest energy barriers (diffusion barriers) of path A and B are 1.152 eV and 0.757 eV, respectively. It is obvious that impurity atoms tend to diffuse through path B owing to the lower energy barrier, and its diffusion activation energy Q is 0.757 eV. The resulting energy curves also show that the energy at the tetrahedral sites is a local minimum of path B. Thus, the diffusion atoms would pass through the metastable tetrahedral sites acting as a transition state.Figure 4
The curve for diffusion barriers of path A and path B. Owing to the symmetry of the present system, only half of the path is shown.
(a)
(b)An experiment testing hydrogen diffusion performance at medium-high temperatures by gas-phase hydrogen permeation methods was also conducted.As it is illustrated in Figure5, an Incoloy 800H sample tube was maintained under high vacuum and several experimental temperatures. Hydrogen flowed from its source to inside the tube and then diffused outwards driven by the concentration difference. A quadrupole mass spectrometer was used to measure hydrogen concentration.Figure 5
Schematic of experimental system.According to the measured hydrogen permeation curve, the temperature-dependent diffusion coefficientD(T) was calculated by the lag time method:(7)DT=d26t,wheredis the thickness of the sample tube andtis the characteristic time, that is, the time for the hydrogen permeation to reach 0.617 times the saturated permeability.The experimental data can be fitted with the logarithmic form of the Arrhenius relation, and the activation energy results are listed in Table2 with their corresponding theoretical values. Considering the parameter difference between experimental and theoretical approach (such as temperature, pressure, gas concentration, and the structural defects of experimental material), this deviation is acceptable. In addition, since the values of QH and QD are similar to each other in experiment, it is reasonable to infer that the value of QT would also be close to QH and QD.Table 2
Comparison between first-principles and experimental diffusion activation energyQ.
Studies
Q
H (eV)
Q
D (eV)
T (°C)
Calculations
0.757
0.757
—
Experimental [15]
0.493 ± 0.037
0.497 ± 0.036
300–550
η
53.55%
52.31%
—


### 3.2. Diffusion Coefficients

As is shown in Figure4, the length of every transition l is 2.584 Å (O1 to O2). Substituting Q and l to (5), D0 of H, D, and T are shown in Table 3.Table 3
The calculation and experimental results of temperature independent preexponential factorD0.
Studies
D
0
H (m2/s)
D
0
D (m2/s)
D
0
T (m2/s)
Calculation
1.56
×
10
-
6
1.10
×
10
-
6
8.99
×
10
-
7
Experimental
6.31
×
10
-
7
5.19
×
10
-
7
—
η
147.23%
115.76%
—Finally, the complete diffusion coefficients of H, D, and T can be written as(8)DH=1.56×10-6×exp⁡-0.757eVkBTm2·s-1DD=1.10×10-6×exp⁡-0.757eVkBTm2·s-1.DT=8.99×10-7×exp⁡-0.757eVkBTm2·s-1It can be seen that the diffusion coefficients of isotopes decrease with the increasing mass. This law is consistent with the experimental results (Figure6). The explanation of this behavior is that the lower vibrational frequency of heavier isotopes reduces the diffusion rate linearly as shown in (3).Figure 6
Computed and experimental diffusion coefficients of hydrogen isotopes.


## 3.1. Diffusion Activation Energy

To confirm the physical diffusion path, the most stable vacancy site should be found in the first step. The system energyEs of H located in the three octahedral sites and the tetrahedral site was calculated (Table 1).Table 1
System energy of H atom located in the four sites for2×2×2 and 3×3×3 supercell models, respectively.
Vacancy sites
O-FeCr
O-FeNi
O-FeCrNi
Te
2 × 2 × 2
E
s (eV)
−251.937
−251.785
−251.641
−251.321
Δ
E
s (eV)
—
0.152
0.296
0.616
3 × 3 × 3
E
s (eV)
−777.936
−777.908
−777.790
−777.321
Δ
E
s (eV)
—
0.028
0.146
0.615In this calculation with VASP, the projector augmented wave (PAW) pseudopotential was employed to describe the electron-ion interaction while neutron effects and mass effects were not involved. Therefore, the energies of H, D, and T should be the same; this is confirmed by our calculation. This is usually well justified as the barrier for an atomic jump is determined by the electronic interaction, which is identical for two isotopes of the same element, and not by the masses of the nuclei [25].It is found that the most stable site is the O-FeCr site whereEs equals to −251.937 eV, lower than the Te site with a 2×2×2 supercell by 0.616 eV. The 3×3×3 supercell gives the same order of Es of the four sites and an essentially same energy difference of 0.615 eV. Therefore, the impurity atoms must transit through an O-FeNi site or an O-FeCrNi site when they complete a whole diffusion process from one O-FeCr site to another. Owing to the lower Es, the O-FeNi site is considered as the midpoint. As the starting, end, and midpoints were confirmed, there are two possible diffusion paths in the model, as illustrated in Figure 3. In path A, a H atom transits from O1 to O2 in a straight line; in path B, the impurity atom jumps by a broken line, passing by the Te position. To compare the two diffusion paths, the “Transition State Search” method, nudged elastic band calculation [26], was used to search for the transition state.Figure 3
Two possible paths of one impurity atom diffusing from an O-FeCr site O1 to another O-FeCr siteO1′. O2 is an O-FeNi site. Te denotes a tetrahedral interstitial site.As shown in Figure4, the highest energy barriers (diffusion barriers) of path A and B are 1.152 eV and 0.757 eV, respectively. It is obvious that impurity atoms tend to diffuse through path B owing to the lower energy barrier, and its diffusion activation energy Q is 0.757 eV. The resulting energy curves also show that the energy at the tetrahedral sites is a local minimum of path B. Thus, the diffusion atoms would pass through the metastable tetrahedral sites acting as a transition state.Figure 4
The curve for diffusion barriers of path A and path B. Owing to the symmetry of the present system, only half of the path is shown.
(a)
(b)An experiment testing hydrogen diffusion performance at medium-high temperatures by gas-phase hydrogen permeation methods was also conducted.As it is illustrated in Figure5, an Incoloy 800H sample tube was maintained under high vacuum and several experimental temperatures. Hydrogen flowed from its source to inside the tube and then diffused outwards driven by the concentration difference. A quadrupole mass spectrometer was used to measure hydrogen concentration.Figure 5
Schematic of experimental system.According to the measured hydrogen permeation curve, the temperature-dependent diffusion coefficientD(T) was calculated by the lag time method:(7)DT=d26t,wheredis the thickness of the sample tube andtis the characteristic time, that is, the time for the hydrogen permeation to reach 0.617 times the saturated permeability.The experimental data can be fitted with the logarithmic form of the Arrhenius relation, and the activation energy results are listed in Table2 with their corresponding theoretical values. Considering the parameter difference between experimental and theoretical approach (such as temperature, pressure, gas concentration, and the structural defects of experimental material), this deviation is acceptable. In addition, since the values of QH and QD are similar to each other in experiment, it is reasonable to infer that the value of QT would also be close to QH and QD.Table 2
Comparison between first-principles and experimental diffusion activation energyQ.
Studies
Q
H (eV)
Q
D (eV)
T (°C)
Calculations
0.757
0.757
—
Experimental [15]
0.493 ± 0.037
0.497 ± 0.036
300–550
η
53.55%
52.31%
—


## 3.2. Diffusion Coefficients

As is shown in Figure4, the length of every transition l is 2.584 Å (O1 to O2). Substituting Q and l to (5), D0 of H, D, and T are shown in Table 3.Table 3
The calculation and experimental results of temperature independent preexponential factorD0.
Studies
D
0
H (m2/s)
D
0
D (m2/s)
D
0
T (m2/s)
Calculation
1.56
×
10
-
6
1.10
×
10
-
6
8.99
×
10
-
7
Experimental
6.31
×
10
-
7
5.19
×
10
-
7
—
η
147.23%
115.76%
—Finally, the complete diffusion coefficients of H, D, and T can be written as(8)DH=1.56×10-6×exp⁡-0.757eVkBTm2·s-1DD=1.10×10-6×exp⁡-0.757eVkBTm2·s-1.DT=8.99×10-7×exp⁡-0.757eVkBTm2·s-1It can be seen that the diffusion coefficients of isotopes decrease with the increasing mass. This law is consistent with the experimental results (Figure6). The explanation of this behavior is that the lower vibrational frequency of heavier isotopes reduces the diffusion rate linearly as shown in (3).Figure 6
Computed and experimental diffusion coefficients of hydrogen isotopes.


## 4. Conclusion

In this study, a first-principles calculation was performed to study the diffusion behavior of tritium and its isotopes in Incoloy 800H. Three octahedral sites and one tetrahedral site were considered, and the most stable vacancy was found to be the octahedral site surrounded by Fe atoms and Ni atoms. It was found that the diffusion of impurity atoms occurred through an indirect path via the metastable tetrahedral site and the octahedral site surrounded by Fe atoms and Cr atoms. The numerical results show that the diffusion activation energies of the three isotopes have the same value of approximately 0.757 eV, and the diffusion coefficients are also estimated; andD0=1.56×10-6, 1.10×10-6, and 8.99×10-7 (m2/s) for H, D, and T, respectively. Owing to the limitations related to experimental conditions, there are some deviations between the experimental and theoretical results.In addition, the diffusion coefficients of tritium in Incoloy 800H and other alloys are still very difficult to obtain experimentally. However, the first-principles calculation provides an effective way to investigate the diffusion mechanism and related parameters. These theoretical results can be used for material selection and nuclear safety evaluation in HTR-PM.



---

*Source: 4038021-2017-10-10.xml*