Abstract
Helium (He) accumulation in tungsten-widely used as a plasma-facing material in fusion reactors-can lead to clustering, trap mutation, and eventual formation of helium bubbles, critically impacting material performance. To clarify the atomic-scale mechanisms governing this process, we conducted systematic molecular statics and molecular dynamics simulations across a wide range of vacancy cluster sizes (n = 1-27) and temperatures (500-2000 K). We identified the onset of trap mutation through abrupt increases in tungsten atomic displacement. At 0 K, the critical helium-to-vacancy (He/V) ratio required to trigger mutation was found to scale inversely with cluster size, converging to ~5.6 for large clusters. At elevated temperatures, thermal activation lowered the mutation threshold and introduced a distinct He/V stability window. Below this window, clusters tend to dissociate; above it, trap mutation occurs with near certainty. This critical He/V ratio exhibits a linear dependence on temperature and can be described by a size- and temperature-dependent empirical relation. Our results provide a quantitative framework for predicting trap mutation behavior in tungsten, offering key input for multiscale models and informing the design of radiation-resistant materials for fusion applications.