Hydrogen clustering in bcc metals: atomic origin and strong stress anisotropy

Hydrogen (H) induced damage in metals has been a long-standing woe for many industrial applications. One form of such damage is linked to H clustering, for which the atomic origin remains contended, particularly for non-hydride forming metals. In this work, we systematically studied H clustering behavior in bcc metals represented by W, Fe, Mo, and Cr, combining first-principles calculations, atomistic and Monte Carlo simulations. H clustering has been shown to be energetically favorable, and can be strongly facilitated by anisotropic stress field, dominated by the tensile component along one of the <001> crystalline directions. We showed that the stress effect can be well predicted by the continuum model based on H formation volume tensor, and that H clustering is thermodynamically possible at edge dislocations, evidenced by nanohydride formation at rather low levels of H concentration. Moreover, anisotropy in the stress effect is well reflected in nanohydride morphology around dislocations, with nanohydride growth occurring in the form of thin platelet structures that maximize one <001> tension. In particular, the <001> type edge dislocation, with the <001> tensile component maximized, has been shown to be highly effective in facilitating H aggregation, thus expected to play an important role in H clustering in bcc metals, in close agreement with recent experimental observations. This work explicitly and quantitatively clarifies the anisotropic nature of stress effect on H energetics and H clustering behaviors, offering mechanistic insights critical towards understanding H-induced damages in metals.

Accurate prediction of nanovoid structures and energetics in bcc metals

Knowledge on structures and energetics of nanovoids is fundamental to understand defect evolution in metals. Yet there remain no reliable methods able to determine essential structural details or to provide accurate assessment of energetics for general nanovoids. Here, we performed systematic first-principles investigations to examine stable structures and energetics of nanovoids in bcc metals, explicitly demonstrated the stable structures can be precisely determined by minimizing their Wigner-Seitz area, and revealed a linear relationship between formation energy and Wigner-Seitz area of nanovoids. We further developed a new physics-based model to accurately predict stable structures and energetics for arbitrary-sized nanovoids. This model was well validated by first-principles calculations and recent nanovoid annealing experiments, and showed distinct advantages over the widely used spherical approximation. The present work offers mechanistic insights that crucial for understanding nanovoid formation and evolution, being a critical step towards predictive control and prevention of nanovoid related damage processes in structural metals.

Predictive Model of Hydrogen Trapping and Bubbling in Nanovoids in BCC Metals

Interplay between hydrogen and nanovoids, despite long-recognized as a central aspect in hydrogen-induced damages in structural materials, remains poorly understood. Focusing on tungsten as a model BCC system, the present study, for the first time, explicitly demonstrated sequential adsorption of hydrogen adatoms on Wigner-Seitz squares of nanovoids with distinct energy levels. Interaction between hydrogen adatoms on the nanovoid surface is shown to be dominated by pairwise power law repulsion. A predictive model was established for quantitative prediction of configurations and energetics of hydrogen adatoms in nanovoids. This model, further combined with equation of states of hydrogen gas, enables prediction of hydrogen molecule formation in nanovoids. Multiscale simulations based on the predictive model were performed, showing excellent agreement with experiments. This work clarifies fundamental physics and provides full-scale predictive model for hydrogen trapping and bubbling in nanovoids, offering long-sought mechanistic insights crucial for understanding hydrogen-induced damages in structural materials.

Hydrogen bubble nucleation by self-clustering: Density Functional Theory and statistical models studies using tungsten as a model system

Low-energy hydrogen irradiation is known to induce bubble formation in tungsten, while its atomistic mechanisms remain little understood. Using first-principles calculations and statistical models, we studied the self-clustering behavior of hydrogen in tungsten. Unlike previous speculations that hydrogen self-clusters are energetically unstable owing to the general repulsion between two hydrogens, we demonstrated that hydrogen self-cluster becomes more favorable as the cluster size increases. We found that hydrogen atoms would form two-dimensional platelet-like structures along {100} planes. These hydrogen self-clustering behaviors can be quantitative understood by the competition between long-ranged elastic attraction and local electronic repulsion. Further statistical analysis showed that there exists a critical hydrogen concentration above which hydrogen self-clusters are thermodynamically stable and kinetically feasible. Based on this critical hydrogen concentration, the plasma loading conditions under which hydrogen self-clusters form were predicted. Our predictions showed excellent agreement with experimental results of hydrogen bubble formation in tungsten exposed to low-energy hydrogen irradiation. Finally, we proposed a possible mechanism for the hydrogen bubble nucleation via hydrogen self-clustering. This work provides mechanistic insights and quantitative models towards understanding of plasma-induced hydrogen bubble formation in plasma-facing tungsten.

First principle study of intrinsic defects in hexagonal tungsten carbide

The characteristics of intrinsic defects are important for the understanding of self-diffusion processes, mechanical strength, brittleness, and plasticity of tungsten carbide, which present in the divertor of fusion reactors. Here, we use first-principles calculations to investigate the stability of point defects and their complexes in WC. Our calculation results confirm that the formation energies of carbon defects are much lower than that of tungsten defects. The outward relaxations around vacancy are found. Both interstitial carbon and interstitial tungsten atom prefer to occupy the carbon basal plane projection of octahedral interstitial site. The results of isolated carbon defect diffusion show that the carbon vacancy stay for a wide range of temperature because of extremely high diffusion barriers, while carbon interstitial migration is activated at lower temperatures for its considerable lower activation energy. These results provide evidence for the presumption that the 800K stage is attributed by the annealing out of carbon vacancies by long-range migration.