Abstract:Helium diffusion in metals is the basic requirement of nucleation and growth of bubble, which gives rise to adverse degradation effects on mechanical properties of structural materials in reactors under irradiation. Multi-scale modeling scheme has been developed to study effects of helium on the long-term microstructural evolution. However, the implementation of Arrhenius law based on the quasi-equilibrium reaction process is not appropriate to predict the migration behavior of helium in metals due to low-energy barrier. A coarse-grained formula is required to incorporate the non-equilibrium nature, e.g., the dissipative friction coefficient $\gamma$. In this paper, we derive an analytical expression for $\gamma$ based on a coarse-grained model of Brownian motion upon a periodic potential, in terms of dissipative feature of the thermal excitations in the many-body dynamical system by constructing an adiabatic relaxation process, which is then confirmed by a numerical example of vacancy migration in BCC W. Then the many-body dynamics simulations are performed for helium migration in BCC W and Fe, where the classical mobility are obtained and in good agreement with the data from experiments and other calculations. Finally, we propose a coarse-grained formula for the helium migration in BCC W and Fe, i.e., Eq.(47) in the context, using the calculated parameters from the adiabatic relaxation simulations. This work would help to develop a new multi-scale modeling scheme for effects of helium in metals, as well as the atomistic reactions with low-energy pathways in materials science.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **How to accurately describe and predict the non - equilibrium migration behavior of helium atoms in body - centered cubic (BCC) tungsten (W) and iron (Fe), especially the migration under low - energy - barrier conditions**. Specifically, the traditional Arrhenius - law - based method has limitations in describing low - energy - barrier reactions, so a new coarse - grained formula is required to better describe the migration behavior of helium atoms in metals.
### Problem Background
Helium (He) is one of the common products in nuclear reactor structural materials, especially in fission and fusion reactors. Due to the insolubility of helium atoms in metals, they are easily trapped and form helium bubbles, which will lead to the degradation of the mechanical properties of materials, especially the helium embrittlement phenomenon that occurs at high temperatures. Therefore, understanding the nucleation and growth kinetics of helium bubbles is crucial for nuclear materials science and engineering.
### Limitations of Existing Methods
Traditional multi - scale modeling schemes usually use the Arrhenius law to describe the atomic activation process. However, the Arrhenius law is only applicable to quasi - equilibrium processes, that is, the effective reaction temperature \(T \ll 1\). For low - energy - barrier reactions (such as the migration of helium in metals), when the effective temperature \(T \gg 1\), the reaction is non - equilibrium, and the traditional Arrhenius law cannot accurately describe this behavior. Therefore, a dissipative friction coefficient \(\gamma\) needs to be introduced to describe the kinetic characteristics of non - equilibrium reactions.
### Main Contributions of the Paper
To solve the above problems, this paper proposes a new coarse - grained formula for describing the classical mobility of helium in BCC W and Fe. The specific steps are as follows:
1. **Derivation of the dissipative friction coefficient \(\gamma\)**:
- Based on the coarse - grained model of Brownian motion, the analytical expression of \(\gamma\) is derived by constructing an adiabatic relaxation process.
- The migration behavior of helium in BCC W and Fe is calculated by multi - body dynamics simulation, and the numerical value of \(\gamma\) is obtained.
2. **Application to helium migration simulation**:
- Using the obtained parameters, a practical coarse - grained formula is proposed to describe the migration behavior of helium in BCC W and Fe.
- The validity of this formula is verified by molecular dynamics simulation and compared with experimental and other computational results.
### Conclusion
This work provides a new multi - scale modeling scheme that can more accurately describe the low - energy - barrier migration behavior of helium in metals, which is of great significance for studying the helium - induced microstructure evolution. In addition, this method can also be applied to other atomic reactions with low - energy - barrier paths, further promoting the development of materials science.
### Formula Summary
- Arrhenius law:
\[
\nu=\nu_0 e^{-\frac{G}{k_B T}}=\nu_0 e^{\frac{S}{k_B}} e^{-\frac{\Phi_m}{k_B T}}
\]
where \(\nu_0\) is the attempt frequency, \(G = \Phi_m - TS\) is the free - energy barrier, \(k_B\) is the Boltzmann constant, and \(T\) is the absolute temperature.
- Total energy dissipation rate in the adiabatic relaxation process:
\[
\frac{dU(t)}{dt}=- 2\gamma K(t)+\eta(t)
\]
where \(U(t)=\Phi(X)+K(t)\) is the total energy, \(K(t)=\frac{m^* \dot{X}^2}{2}\) is the kinetic energy, and \(\eta(t)=f(t)\dot{X}\) is the instantaneous power of the random force.
- Characteristic dissipation time \(\tau\):
\[
\tau = \gamma^{-1}=\int_0^\infty \frac{U(t)}{U_0} dt
\]
- Classical mobility \(\mu\) of helium:
\[
\mu=\frac{1}{k_B T} \lim_{t \to \infty} \left\l