First-Principles Study on Various Point Defects Formed by Hydrogen and Helium Atoms in Tungsten.
Tungsten (W) is widely considered as the most potential plasma facing material (PFM)  and the first wall material (FWM) in fusion reactor due to its excellent properties, such as high melting point, high mechanical strength at high temperatures, good thermal conductivity, low thermal expansion coefficient, and high sputtering threshold energy . Studies on tungsten receive an increasing amount of attention from scientists and engineers. Tungsten is subjected to high particle fluxes in fusion reactors, not only the plasma background ions (H, D, and T) but also the intrinsic impurities (He, Be, and C). Because hydrogen (H) is one of the major plasma background ions, hydrogen irradiation may result in the change of the mechanical properties of tungsten. Helium (He) irradiation leads to the blister formation and the subsequent degradation of the mechanical properties of metals. Thus, it is very important to investigate the trapping and blister formation of hydrogen or helium in tungsten and the behaviors of hydrogen and helium in tungsten.
A number of experimental and theoretical studies were carried out to understand the interaction between the hydrogen or helium impurity and the tungsten host lattice . For example, Arkhipov et al. found that hydrogen could form ablister on the tungsten surface . You presented the experimental results in a divertor plasma simulator NAGDIS-II and observed the formation of hydrogen bubbles on the surface of tungsten . Grigorev et al. observed the nucleation and growth of hydrogen bubbles on dislocations in tungsten under high-flux, low energy plasma exposure . Liu et al. studied the behavior of hydrogen in tungsten through first-principles calculation . Wang et al. studied the behavior of helium on tungsten surfaces and the hydrogen in tungsten by using first-principles calculation . Li et al. studied the behavior of hydrogen or helium in tungsten by using the molecular dynamics simulation method . With the rapid development of the computer technique, the first-principles calculation method plays an important role in the materials research; it can investigate the mechanical properties, electronic structure, and optical properties of some new kinds of materials [9-11].
There are a lot of works focused on defects in tungsten. For example, research work on the interaction between two hydrogen atoms that are located in tetrahedral interstitial sites in tungsten has been completed . Some researchers investigated the interaction between hydrogen or helium atoms and the vacancy . However, a systematic study of the interaction between two hydrogen or helium atoms in tungsten without other defects is still lacked. In this paper, by using first-principles calculation method based on density functional theory (DFT), we have studied the behavior of hydrogen atoms and helium atoms in body centered cubic (BCC) tungsten. We calculated the formation energy of single hydrogen atom and two hydrogen atoms in tungsten and the formation energy of single helium atom and two helium atoms in tungsten. Using the charge density difference and atomic population, we analyzed the reason why there is a repulsive interaction between two hydrogen atoms and an attractive interaction between two helium atoms. Our computational results might shed some light on the development of tungsten-based materials that are used in ITER and other fusion rectors. For example, Cheng's research group investigated the blistering behavior and deuterium retention in tungsten exposed to high-flux deuterium-neon mixed plasmas by experiment research  (the temperature of target surface is 10 K) after a lot of first-principles calculations were completed; the results of Nguyen-Manh and Dudarev's  first-principles calculation are in good agreement with experimental data derived from thermal desorption spectroscopy  (annealing-temperature is about 2150 K).
2. Models and Computational Details
Tungsten belongs to the BCC crystal system, its space group is Im-3m. The simulation supercell composed of 24 lattice points (2x3x2) was used in this paper. Initially, we considered three possible sites for a single H atom and He atom entering the lattice of BCC W. According to their positions, we marked the H and He atoms as [H.sub.Sub], [H.sub.Tet], [H.sub.Oct], [He.sub.Sub], [He.sub.Tet], and [He.sub.Oct] when the H atom and He atom were located in the substitutional site, the tetrahedral interstitial site, and the octahedral interstitial site, respectively. The calculation models that contain a single H atom or a He atom are shown in Figure 1. To explore the interaction between two H atoms or two He atoms, we build a series of two H atoms or two He atoms models. The symbols [H.sup.21.sub.Sub], [H.sup.22.sub.Sub], and [H.sup.23.sub.Sub] stand for the first-nearest-neighbor, second-nearest-neighbor, and third-nearest-neighbor H atom around the preexisting H atom that is located in the substitutional site, and we use the symbols [H.sup.21.sub.Tet], [H.sup.22.sub.Tet], [H.sup.23.sub.Tet], [H.sup.22.sub.Oct], [H.sup.23.sub.Oct], and HOC standing for two H atoms that are located in the different tetrahedral interstitial sites and octahedral interstitial sites, respectively. We use the symbols [He.sup.21.sub.Sub], [He.sup.22.sub.Sub], [He.sup.23.sub.Sub], [He.sup.21.sub.Tet], [He.sup.22.sub.Tet], [He.sup.23.sub.Tet], [He.sup.21.sub.Oct], [He.sup.22.sub.Oct], and [He.sup.23.sub.Oct] standing for the first-nearest-neighbor, second-nearest- neighbor, and third-nearest-neighbor He around the preexisting He atom that is located in substitutional, tetrahedral, and octahedral sites, respectively. These models contain two H atoms or two He atoms and are shown in Figure 2.
First-principles calculations were performed by using the DFT  and the plane wave pseudo-potential technique, as implemented in the Cambridge Sequential Total Energy Package (CASTEP) . The generalized gradient approximation (GGA)  with Perdew-Burke-Ernzerhof (PBE)  functional was used to describe the exchange-correlation interaction among the electrons, and ultrasoft pseudopotentials were employed for ion-electron interaction. For W, the four 5d electrons are considered to be the valence electrons together with two 6s electrons (the reference state is 5[d.sup.4]6[s.sup.2]). During the calculation, the supercell size, shape, and atomic position are relaxed to equilibrium. The plane wave energy cutoff is 400 eV for all calculations. The basic parameters were chosen as follows: the space representation is reciprocal, SCF tolerance equals 1.0 x [10.sup.-6] eV/atom, and k sampling with 5 x 5 x 5 k-point mesh in the Brillouin zone was used, after the convergence test was completed. The total energy of the supercell deceases as the plane wave cutoff and k-point mesh increasing; the structure also becomes more stable. Considering the optimization of computing resource, we set the plane wave cutoff and k-point mesh as mentioned above. Our parameters setting is in agreement with that of Liu's group ; our calculated formation energy of a H or He atom in W is in good agreement with Liu et al.'s results. Moreover, the increase of these two calculation parameters does not significantly change the total energy of the supercell. The optimal atomic positions were determined by satisfying the following conditions: (1) the maximum interatomic forces were smaller than 0.05 eV/nm; (2) the maximum change of energy per atom was smaller than 1.0 x [10.sup.-5] eV; (3) the maximum displacement was smaller than 0.001 A; and (4) the maximum stress of the crystal was smaller than 0.02 GPa. All the properties that we calculated based on the crystal structure have been optimized.
3. Results and Discussion
3.1. Optimization of the Crystal Structure. In order to study various defects formed by H or He atoms in W, first of all, we optimized the crystal structure of the W supercell with and without the various defects formed by H or He atoms. Table 1 shows the calculated results, lattice parameters, elastic constants, and elastic modulus of the BCC W. The result shows that our results are in good agreement with the theoretical simulation results  and the experimental data . Figure 3 shows the change of the volume of BCC W when the H or He atoms were introduced in W. From Figure 3, we can find that the volume of W decreases when the H or He atoms are located in the substitutional sites while the volume increases when the H or He atoms are located in the tetrahedral and octahedral sites. The reason for the volume decrease is that the atomic radius of H or He is much smaller than that of W, and the volume increases due to the extra atoms that enter the W supercell.
3.2. Formation Energy. The formation energy of various defect configurations, [E.sup.f.sub.defect], was defined as follows:
[E.sup.f.sub.defect] = [E.sup.tot.sub.system] - [E.sup.24.sub.W] (N/24) - [n.sub.He][E.sub.He] - [n.sub.H][E.sub.H], (1)
where [E.sup.tot.sub.system] is the total energy of the system with defects, N, [n.sub.He], and [n.sub.H] are the numbers of W, He, and H atoms, respectively, [E.sup.24.sub.W] is the total energy of a perfect 2 x 3 x 2 supercell of BCC W with 24 atoms, [E.sub.He] is the energy of an isolated He atom, and [E.sub.H] is the total energy of an isolated H atom. The formation energy of various single point defects in W is summarized in Table 2, and the results are in good agreement with other calculation results [1, 20, 21].
According to Table 2, we plotted the formation energies of the single point defect in Figure 4. For He, the substitutional site is the most favorable over the tetrahedral and octahedral sites. For H, the tetrahedral site is the most preferred site over the substitutional and octahedral sites. Interestingly, the formation energies of [H.sub.Tet] and [H.sub.Oct] are negative, which implies that defects are more stable than other defects, and the formation process of the defects, [H.sub.Tet] and [H.sub.Oct], is an exothermic process. Figure 5 shows the formation energy of two H and He atoms in W, and these two atoms are located in different sites. Figure 5(a) shows the formation energy of two H atoms that are located in different substitution, tetrahedral, and octahedral sites. According to the distance between two H atoms, they are marked as the first-, second-, and third-nearest-neighbor atoms. Figure 5(b) shows the formation energy of two He atoms located in different substitutional, tetrahedral, and octahedral sites. For two H atoms, the highest formation energy is the H atoms located in different substitutional sites, and the lowest formation energy is closely related to the sites of the two H atoms. The formation energies of the two H atoms located in different octahedral sites are the lowest when the two H atoms are the first- or third-nearest-neighbor atoms, the formation energy of the two H atoms that located in different tetrahedral sites is the lowest, and the two H atoms are the second-nearest-neighbor atoms. For He atoms, the lowest formation energy can be got when the two He atoms located in different substitutional sites. The formation energy of two He atoms located in different tetrahedral sites is the highest when the two He atoms are the first- or second-nearest-neighbor atoms. When two He atoms are the third-nearest-neighbor atoms, the highest formation energy can be got when the two He atoms located in different octahedral sites. The formation energy of the two foreign atoms (H or He) is not equal to twice value of the single atom, and the formation energy of the two foreign atoms is two times smaller than the value of the formation energy of a single atom. In short, the formation energy of two H and He atoms are determined by their located sites.
3.3. Interaction of H-H and He-He. The binding energy can reflect the interaction between two H or He atoms in W. The binding energy of two H or He atoms is calculated as follows:
[E.sub.b] ([X.sub.2]) = 2E(X)-E([X.sub.2])-[E.sup.24.sub.W], (2)
where E(X) is the total energy of the W supercell with a H or He atom, E([X.sub.2]) is the total energy of the supercell with two H or He atoms, and [E.sup.24.sub.W] is the total energy of 24 W atoms. A positive binding energy indicates an attraction exists between two atoms, while a negative binding energy shows a repulsion interaction exists between two atoms. The binding energies of two H or He atoms are shown in Figure 6. The most binding energies existing between two H atoms are near the zero while the most binding energies existing between two He atoms are in a positive value. The results show that there is an attraction interaction among the He atoms, and the interaction between two H atoms is very weak.
To further study the interaction between two H or He atoms, we recorded the distance between two H and He atoms before and after relaxation. Figure 7 shows the change of the distance between two H or He atoms after relaxation, and the atoms are located in different sites. The distance between two H atoms becomes longer after relaxation than the values before relaxation. This result suggests that a repulsive interaction exists between H atoms when the H atoms are close to each other, because each H atom in W lattice carries certain amount of negative charge. The distance between two H atoms is larger than 0.75 [Angstrom] (the distance between two H atoms in a [H.sub.2] molecule) after relaxation; these results indicate there is no [H.sub.2] molecule formed. The result of our study is in good agreement with the results of Liu et al. . The distance between two He atoms becomes shorter after relaxation. This result indicates that there is an attraction interaction between the two He atoms when they are close to each other.
3.4. Analysis of the Mulliken Population and Charge Density. The electronic structure can help us understand the root causes of the interaction between two H and He atoms, so the Mulliken population and the charge density were analyzed in this paper. Table 3 shows the atomic population of the H and He atoms in W, and two H and He atoms are the first-nearest-neighbor atoms. Mulliken population analysis shows that certain electrons will transfer from W to H, and H atoms can stand in a stable state. The H atoms will accept 0.10 electrons from the surrounding W atoms when two H atoms are located in substitutional sites. The H atoms accept 0.24 electrons when two H atoms are located in the tetrahedral sites, and they accept 0.26 electrons when they are located in the octahedral sites. The amount of charge transfer is insensitive to the interstitial H and He atoms' positions. The H atoms have a negative charge. This implies that there is a Coulomb repulsion between two H atoms. On the contrary, W will get some electrons from the He atoms. W atoms will accept 0.05 electrons from a He atom when the He atoms are located in the substitution sites; W atoms will get 0.11 electrons from a He atom when He atoms are located in the tetrahedral and octahedral sites. Thus, the Van der Waals force between two He atoms becomes stronger when they are close to each other. The more electrons the W atoms get, the more stable the structure is. Interestingly, we found that the H-H repulsive interaction fits well to the screened Coulomb potential. To explain the differences between the interactions between two H atoms and He atoms, charge density distribution was used. We selected a single atom and two atoms located in the octahedral site. Figure 8 shows the charge density of one and two atoms. Figures 8(a) and 8(d) show the charge density of single H and He atom in W, respectively. Figures 8(b) and 8(c) show two H atoms in W, and Figures 8(e) and 8(f) also show two He atoms in W from several different perspectives. The bonds between H and W contain covalent components. The bond is stronger than that between He and W, and the charge density difference between H and W is bigger than that between He and W. The H atoms are more likely to be attracted by W atoms than by each other while the He atoms are likely to attract each other. There is only one electron in the 1s orbital of an H atom, and the H atom obtains the electrons from the W atoms to maintain a relatively stable state. He atom is a closed shell atom. It prefers to occupy the low electron density zone (interstitial sites).
Using first-principles calculation, we studied the BCC W with single and two point defects formed by H or He atoms, respectively. The results of our calculations are as follows:
(1) For a single H atom, the most preferred site is the tetrahedral interstitial site. For He, the substitutional site is more favorable than the tetrahedral and octahedral sites.
(2) For two H and He atoms, the formation energy is closely related to the sites of the two atoms. The formation energy of two H atoms that were located in different substitutional sites is the highest, while the formation energy of He atoms is lowest when the two He atoms are located in different substitutional sites.
(3) After relaxation, two H atoms are away from each other because the H atom accepts some charge from the surrounding W atoms and they form a covalent bond. The He atoms are close to each other because the W atoms get some electrons from the He atoms, and the more electrons the W atoms get, the more stable the crystal structure is.
Our results are good in agreement with other simulation results and experimental data. This paper might shed some light on the development of tungsten materials used in ITER and other fusion reactors.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This work was supported by the Fundamental Research Funds for the Central Universities under Grant no. 2017MS079.
 S. C. Lee, J. H. Choi, and J. G. Lee, "Energetics of He and H atoms with vacancies in tungsten: First-principles approach," Journal of Nuclear Materials, vol. 383, no. 3, pp. 244-246, 2009.
 C. H. Skinner, A. A. Haasz, V. K. H. Almov et al., "Recent advances on hydrogen retention in iter's plasma-facing materials: Beryllium, carbon. And tungsten," Fusion Science and Technology, vol. 54, no. 4, pp. 891-945, 2008.
 1.1. Arkhipov, S. L. Kanashenko, V M. Sharapov, R. K. Zalavutdinov, and A. E. Gorodetsky, "Deuterium trapping in iondamaged tungsten single crystal," Journal of Nuclear Materials, vol. 363-365, no. 1-3, pp. 1168-1172, 2007.
 Y. M. You, "Effects of low energy and high flux helium/hydrogen plasma irradiation on tungsten as plasma facing material," Plasma Science and Technology, vol. 7, no. 3, pp. 2828-2834, 2005.
 P. Grigorev, D. Terentyev, V. Dubinko et al., "Nucleation and growth of hydrogen bubbles on dislocations in tungsten under high flux low energy plasma exposure," Nuclear Instruments and Methods in Physics Research, Section B: Beam Interactions with Materials and Atoms, vol. 352, pp. 96-99, 2015.
 Y. L. Liu, W. Lu, A. Y. Gao, L. J. Gui, and Y. Zhang, "First-principles investigation of diffusion behaviours of H isotopes: from W(110) surface into bulk and in bulk W," Chinese Physics B, vol. 21, no. 12, 2012.
 X. X. Wang, Y. Zhang, H. B. Zhou, and J. L. Wang, "Effects of niobium on helium behaviors in tungsten: a first-principles investigation," Acta Physica Sinica, vol. 63, no. 4, article 046103, 2014.
 M. Li, J. Cui, J. Wang, and Q. Hou, "Molecular dynamics simulations of cumulative helium bombardments on tungsten surfaces," Nuclear Instruments and Methods in Physics Research, Section B: Beam Interactions with Materials and Atoms, vol. 337, pp. 45-54, 2014.
 Q. Zhao, Z. Zhang, Y. Li, and X. P. Ouyang, "The mechanical and thermodynamic properties of [beta]-[si.sub.1-x][c.sub.x]," RSC Advances, vol. 7, no. 45, pp. 28499-28505, 2017.
 D. Y. Jiang, Z. Zhang, R. X. Liang et al., "Electronic structure and optical properties of the scintillation material wurtzite ZnS(Ag)," Nuclear Science and Techniques, vol. 28, no. 3, pp. 35-41, 2017.
 Z. Zhang, Q. Zhao, Y. Li, and X. P. Ouyang, "Electronic structure and optical properties of CsI, CsI(Ag), and CsI(Tl)," Journal of the Korean Physical Society, vol. 68, no. 9, pp. 1069-1074, 2016.
 J. Xu and J. Zhao, "First-principles study of hydrogen in perfect tungsten crystal," Nuclear Instruments and Methods in Physics Research, Section B, vol. 267, no. 18, pp. 3170-3174, 2009.
 L. Cheng, G. De Temmerman, T. W. Morgan et al., "Mitigated blistering and deuterium retention in tungsten exposed to high-flux deuterium-neon mixed plasmas," Nuclear Fusion, vol. 57, no. 4, article 046028, 2017.
 D. Nguyen-Manh and S. L. Dudarev, "Trapping of He clusters by inert-gas impurities in tungsten: First-principles predictions and experimental validation," Nuclear Instruments and Methods in Physics Research, Section B, vol. 352, pp. 86-91, 2015.
 E. V. Kornelsen and A. A. Van Gorkum, "Study of bubble nucleation in tungsten using thermal desorption spectrometry: clusters of 2 to 100 helium atoms," Journal of Nuclear Materials, vol. 92, no. 1, pp. 79-88, 1980.
 P. Hohenberg and W. Kohn, "Inhomogeneous electron gas," Physical Review, vol. 136, no. 3B, pp. B864-B871,1964.
 S. Clark, M. Segall, C. Pickard, P. Hasnip, M. Probert, and K. Refson, "First principles methods using castep," Zeitschrift fur Kristallographie, vol. 13, pp. B864-B871, 1964.
 J. P. Perdew and W. Yue, "Accurate and simple density functional for the electronic exchange energy: generalized gradient approximation," Physical Review B, vol. 33, no. 12, pp. 8800-8802, 1986.
 J. P. Perdew, K. Burke, and M. Ernzerhof, "Generalized gradient approximation made simple," Physical Review Letters, vol. 77, no. 18, pp. 3865-3868, 1996.
 Y. L. Liu, Y. Zhang, G. N. Luo, and G. H. Lu, "Structure, stability and diffusion of hydrogen in tungsten: a first-principles study," Journal of Nuclear Materials, vol. 390-391, no. 1, pp. 1032-1034, 2009.
 C. S. Becquart and C. Domain, "Ab initio calculations about intrinsic point defects and He in W," Nuclear Instruments and Methods in Physics Research, Section B, vol. 255, no. 1, pp. 2326, 2007.
 P. Soderlind, O. Eriksson, J. M. Wills, and A. M. Boring, "Theory of elastic constants of cubic transition metals and alloys," Physical Review B, vol. 48, no. 9, pp. 5844-5851, 1993.
Qiang Zhao, (1) Zheng Zhang, (1) Yang Li, (1) and Xiaoping Ouyang (1,2,3)
(1) Beijing Key Laboratory of Passive Safety Technology for Nuclear Energy, School of Nuclear Science and Engineering, North China Electric Power University, Beijing 102206, China
(2) Northwest Institute of Nuclear Technology, Xi'an 710024, China
(3) School of Materials Science and Engineering, Xiangtan University, Xiangtan 411105, China Correspondence should be addressed to Qiang Zhao; email@example.com
Received 2 April 2017; Revised 10 July 2017; Accepted 13 July 2017; Published 7 September 2017 Academic Editor: Eugenijus Uspuras
Caption: FIGURE 1: The crystal structure of BCC W. The blue balls stand for the tungsten atoms, and the red ball stands for the H atom or He atom. A pure W supercell (a) and a W supercell contain a H atom or He atom in the substitutional site (b), tetrahedral site (c), and octahedral site (d).
Caption: FIGURE 2: The crystal structure of BCC W contains two H atoms or He atoms. The blue balls and red balls stand for W atoms and the H atoms or He atoms, respectively. (a), (b), and (c) show two atoms located in the substitutional site. (d), (e), and (f) show two atoms located in the tetrahedral site, and (g), (h), and (i) show two atoms located in the octahedral site. The numbers in the figure stand for the distances (in A) between the two foreign atoms.
Caption: FIGURE 3: The volume (in A3) of the supercell (2 x 3 x 2) of the W supercell with single and two point defects formed by H and He atoms. The dashed line represents the volume of the perfect BCC W supercell.
Caption: FIGURE 4: The formation energy of the single defects in W.
Caption: FIGURE 5: The formation energy of two H (a) and He (b) atoms located in the different sites.
Caption: FIGURE 6: The binding energy of two H (a) and He (b) atoms in W when they are located in different positions and the distance of two H and He atoms is increased.
Caption: FIGURE 7: The distance (in A) between the two H (a) or He (b) atoms before the relaxation and after relaxation.
Caption: FIGURE 8: Charge density distribution. (a), (b), and (c) show the single H atom and two H atoms located in octahedral sites; (d), (e), and (f) show the single He atom and two He atoms located in different octahedral sites. The red color represents high charge density, and the blue color represents low valence charge density.
TABLE 1: Calculated and experimental lattice constants (a in [Angstrom]), elastic constants ([C.sub.ij] in GPa), bulk modulus (B in GPa), Young's modulus (E in GPa), shear modulus (G in GPa), Poisson's ratio (v), and G/B for BCC W. Parameters a [C.sub.11] [C.sub.12] [C.sub.44] B This paper 3.18 528.06 196.73 176.28 307.17 Cal.  3.17 536.32 202.25 138.70 313.61 Exp.  3.16 533 205 163 314.33 Parameters E G v G/B This paper 421.26 171.95 0.27 0.560 Cal.  386.83 149.42 0.29 0.478 Exp.  417.80 163.40 0.28 0.54 TABLE 2: Summary of the formation energy of the single defect in W. Configuration [H.sub.Sub] [H.sub.Tet] [H.sub.Oct] [He.sub.Sub] This paper 0.965 -2.327 -2.034 5.059 Cal.  0.92 -2.47 -2.07 5.00 Cal. [20, 21] 0.78 -2.44 -2.06 4.70 Configuration [He.sub.Tet] [He.sub.Oct] This paper 6.365 6.583 Cal.  6.23 6.48 Cal. [20, 21] 6.16 6.38 TABLE 3: Atomic populations (Mulliken) of H and He atoms in W. Sites Species Ion s Charge (e) [H.sup.21.sub.Sub] H 1 1.10 -0.10 H 2 1.10 -0.10 [H.sup.21.sub.Sub] He 1 1.95 0.05 He 2 1.95 0.05 [H.sup.21.sub.Tet] H 1 1.24 -0.24 H 2 1.10 -0.24 [H.sup.21.sub.Tet] He 1 1.89 0.11 He 2 1.89 0.11 [H.sup.21.sub.Oct] H 1 1.26 -0.26 H 2 1.26 -0.26 [H.sup.21.sub.Oct] He 1 1.89 0.11 He 2 1.89 0.11
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Zhao, Qiang; Zhang, Zheng; Li, Yang; Ouyang, Xiaoping|
|Publication:||Science and Technology of Nuclear Installations|
|Date:||Jan 1, 2017|
|Previous Article:||Dynamic Modeling and Control Characteristics of the Two-Modular HTR-PM Nuclear Plant.|
|Next Article:||Development and Application of a New High-Efficiency Sparse Linear System Solver in the Thermal-Hydraulic System Analysis Code.|