Energetics of silicon in the bulk and near surfaces of tungsten: a first-principles study

Siliconization of the tokamak walls is a candidate method to improve plasma confinement in fusion tokamaks containing tungsten plasma facing components (W PFCs). To understand the interactions of silicon (Si) with W, the Si behavior in bulk W, and near three low-index W surfaces ((100), (110) and (111)) has been investigated using first-principles density functional theory. In bulk W, Si interstitial atoms have a low solution ability and high mobility, and Si atoms can be strongly trapped by W vacancies. The interaction between two Si adatoms is responsible for the stability of adatom superstructures on W surfaces, consistent with previous experimental observation (Tsong and Casanova 1981 Phys. Rev. Lett. 47 113). Although the coverage dependence of Si adsorption and diffusion energetics on surfaces is related to surface orientation, the W(110) surface has lower Si adsorption affinity and higher Si diffusivity than either the W(111) or W(100) surfaces. The most stable Si adatom superstructure on W surfaces is: square c(2 × 2) pattern on W(100) covered with 0.5 ML Si; rectangular c(4 × 2) pattern on W(110) with 0.25 ML Si; and rhombus p(1 × 1) pattern on W(111) with 1 ML Si. The coverage dependence of Si mobility on/toward W surfaces is generally related to the stability of the Si superstructures as a function of coverage on each surface. Interestingly, Si adatoms prefer to transport below the surface and into W subsurface by an exchange mechanism with W atoms, indicating the likelihood of epitaxial growth of W silicide layers on W surfaces during the operation of W PFCs.