Effect of temperature on hydrogen diffusion mechanism in tungsten: A molecular dynamics simulation study

In this research, the diffusion of hydrogen atoms in the crystalline structure of tungsten was investigated using LAMMPS molecular dynamics simulation code. To describe the interatomic interactions of the W–H system, the EAM potential was used. Hydrogen atoms were placed in the tetrahedral sites of a perfect BCC tungsten lattice to simulate a realistic impurity distribution with a concentration of 2 %. Simulations were performed in a temperature range of 1400–2700K. After structure optimization, the Mean Squared Displacement (MSD) was calculated using the Einstein relation to determine the diffusion coefficients for each temperature. The results showed that with increasing temperature, the hydrogen diffusion coefficient increases exponentially and verifies the Arrhenius relationship. The effective activation energy parameter is calculated 1.48 eV, with a pre-exponential factor of 3.2×10 −6 m 2 /s. Physical analysis revealed three distinct diffusion regimes: at low temperatures, hydrogen mobility is limited by trapping effects; at intermediate temperatures, the TIS-TIS pathway is the dominant mechanism; and at high temperatures, the transition to TIS-OIS-TIS pathways is activated, leading to a sharp increase in the diffusion coefficient. The high value of the effective activation energy is attributed to the collective motion and interactions of the hydrogen atoms at this concentration. These results are applicable in predicting the behavior of tungsten under the high-temperature conditions of fusion reactors. • Molecular dynamics simulations were performed to investigate hydrogen diffusion in single-crystal tungsten. • Simulations were conducted at temperatures from 1400 K to 2700 K using an EAM1 potential. • Hydrogen atoms were uniformly distributed in the tungsten lattice without random clustering. • Diffusion coefficients were calculated from MSD data and fitted to the Arrhenius equation to obtain D 0 and E a . • Hydrogen diffusion rates increase with temperature; single atoms have higher mobility than clusters.

Authors

Institutions

Publication Details

Journal
Nuclear Analysis
Published
2025-09-01
DOI
https://doi.org/10.1016/j.nucana.2025.100186
Primary Topic
Advanced materials and composites
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Effect of temperature on hydrogen diffusion mechanism in tungsten: A molecular dynamics simulation study

Dariush Rostamifard, Ehsanollah Noori, Maryam Ebrahimi, Mir mohammadreza Seyedhabashi et al.
Nuclear Analysis
Advanced materials and composites
article

Effect of temperature on hydrogen diffusion mechanism in tungsten: A molecular dynamics simulation study

Dariush Rostamifard, Ehsanollah Noori, Maryam Ebrahimi, Mir mohammadreza Seyedhabashi, Ali Reza Asle Zaeem, Reza Goodarzi, Amir raeisdana
article en

Abstract

In this research, the diffusion of hydrogen atoms in the crystalline structure of tungsten was investigated using LAMMPS molecular dynamics simulation code. To describe the interatomic interactions of the W–H system, the EAM potential was used. Hydrogen atoms were placed in the tetrahedral sites of a perfect BCC tungsten lattice to simulate a realistic impurity distribution with a concentration of 2 %. Simulations were performed in a temperature range of 1400–2700K. After structure optimization, the Mean Squared Displacement (MSD) was calculated using the Einstein relation to determine the diffusion coefficients for each temperature. The results showed that with increasing temperature, the hydrogen diffusion coefficient increases exponentially and verifies the Arrhenius relationship. The effective activation energy parameter is calculated 1.48 eV, with a pre-exponential factor of 3.2×10 −6 m 2 /s. Physical analysis revealed three distinct diffusion regimes: at low temperatures, hydrogen mobility is limited by trapping effects; at intermediate temperatures, the TIS-TIS pathway is the dominant mechanism; and at high temperatures, the transition to TIS-OIS-TIS pathways is activated, leading to a sharp increase in the diffusion coefficient. The high value of the effective activation energy is attributed to the collective motion and interactions of the hydrogen atoms at this concentration. These results are applicable in predicting the behavior of tungsten under the high-temperature conditions of fusion reactors. • Molecular dynamics simulations were performed to investigate hydrogen diffusion in single-crystal tungsten. • Simulations were conducted at temperatures from 1400 K to 2700 K using an EAM1 potential. • Hydrogen atoms were uniformly distributed in the tungsten lattice without random clustering. • Diffusion coefficients were calculated from MSD data and fitted to the Arrhenius equation to obtain D 0 and E a . • Hydrogen diffusion rates increase with temperature; single atoms have higher mobility than clusters.

Nuclear AnalysisVol. 4(3)
Atomic Energy Organization of Iran (IR), Amirkabir University of Technology (IR), Iranian Center for Quantum Technologies (IR), Nuclear Science and Technology Research Institute
Openalex Percentile: Top 75%
Advanced materials and composites
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.