Effect of Vacancies on Hydrogen Mobility and Trapping in Elemental Fe and Cr: A DFT and kMC Study
Pith reviewed 2026-06-29 11:39 UTC · model grok-4.3
The pith
Vacancy defects reduce hydrogen mobility and raise diffusion activation energy more in Cr than in Fe
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Vacancy defects significantly reduce hydrogen mobility and increase the effective activation energy, with a more pronounced effect observed in Cr due to stronger trapping.
What carries the argument
Multiscale DFT-kMC framework where DFT supplies binding energies and migration barriers for hydrogen-vacancy complexes and kMC evolves diffusion over extended scales
If this is right
- Hydrogen diffusion rates drop and activation energies rise when vacancies are present
- The reduction in mobility is larger in chromium than in iron
- Detrapping barriers fall as more hydrogen atoms occupy a vacancy except for the sixth atom in iron
- The combined approach supplies atomic-scale detail on trapping and detrapping that controls embrittlement
Where Pith is reading between the lines
- Similar vacancy effects may appear in other BCC metals and could be used to tune hydrogen transport in alloys
- The stronger trapping in chromium suggests chromium-rich regions in steels might locally retain hydrogen
- The finite barrier found for the sixth hydrogen in iron could be checked by targeted DFT calculations with different functionals
Load-bearing premise
The DFT binding energies and migration barriers accurately represent real material behavior at the temperatures and concentrations relevant to kMC simulations without significant errors from functional choice or finite-size effects.
What would settle it
Direct experimental measurement of hydrogen diffusion coefficients and activation energies in iron and chromium samples prepared with controlled vacancy concentrations.
read the original abstract
Hydrogen-vacancy interactions play an important role in governing hydrogen transport and embrittlement in body-centered cubic (BCC) metals. In this study, a multiscale approach combining density functional theory (DFT) and kinetic Monte Carlo (kMC) simulations is employed to investigate hydrogen behavior in BCC Fe and Cr. The DFT-calculated binding energies and Bader charge analysis indicate stronger hydrogen trapping in Cr than in Fe. Migration and detrapping energy barriers are determined using the climbing-image nudged elastic band method, showing that the detrapping energy generally decreases with increasing hydrogen occupancy. However, the sixth hydrogen atom in Fe exhibits a finite barrier, contrary to some previous reports. kMC simulations are then used to evaluate hydrogen diffusion over extended time and length scales. The results demonstrate that vacancy defects significantly reduce hydrogen mobility and increase the effective activation energy, with a more pronounced effect observed in Cr due to stronger trapping. The combined DFT-kMC framework provides detailed insight into the mechanisms of hydrogen trapping, detrapping, and diffusion in BCC metals, offering important implications for understanding hydrogen embrittlement in structural materials.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript employs DFT calculations (including CI-NEB for migration and detrapping barriers and Bader charge analysis) together with kMC simulations to examine hydrogen-vacancy interactions in BCC Fe and Cr. It reports stronger hydrogen trapping in Cr than in Fe, a general decrease in detrapping barriers with increasing hydrogen occupancy, a finite barrier for the sixth hydrogen atom in Fe (contrary to some prior reports), and kMC results showing that vacancies reduce hydrogen mobility while raising the effective activation energy, with a more pronounced effect in Cr.
Significance. If the DFT-derived binding energies and barriers are transferable to the kMC model, the multiscale framework supplies quantitative mechanistic insight into hydrogen trapping, detrapping, and long-time diffusion in BCC metals, with direct relevance to hydrogen embrittlement in structural alloys. The combination of atomistic DFT data with extended-scale kMC is a methodological strength.
minor comments (1)
- [Abstract] The abstract states that the sixth hydrogen atom in Fe exhibits a finite barrier 'contrary to some previous reports' but does not identify or cite those reports; this citation should be supplied in the main text.
Simulated Author's Rebuttal
We thank the referee for their summary of our work and for recognizing the methodological value of the DFT-kMC multiscale approach in providing mechanistic insight into hydrogen trapping and diffusion in BCC Fe and Cr. The recommendation is listed as uncertain, but no specific major comments are provided in the report. We therefore have no individual points to rebut or revise at this stage. We remain available to address any additional questions or clarifications the referee may wish to raise.
Circularity Check
No significant circularity identified
full rationale
The derivation chain consists of independent DFT computations (binding energies via Bader analysis, NEB migration/detrapping barriers) that supply numerical inputs to separate kMC rate tables; the kMC outputs (mobility reduction, effective activation energies) are therefore computed consequences rather than redefinitions or statistical fits of those inputs. No self-definitional equations, fitted-input predictions, load-bearing self-citations, or ansatz smuggling appear in the described workflow. The approach is a standard, externally falsifiable multiscale pipeline whose central claims remain logically independent of the final results.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption DFT with chosen functional and pseudopotentials yields binding and migration energies accurate enough for kMC timescales
- domain assumption Bader charge analysis correctly partitions electron density to quantify trapping strength
Reference graph
Works this paper leans on
-
[1]
Dadfarnia, M., P. Novak, D. C. Ahn, J. B. Liu, P. Sofronis, D. D. Johnson, and I. M. Robertson., ‘Recent Advances in the Study of Structural Materials Compatibility with Hydrogen’, Adv. Mater., vol. 22, no. 10, pp. 1128–1135, 2010, doi: 10.1002/adma.200904354
-
[2]
Song and W
J. Song and W. A. Curtin, ‘Atomic mechanism and prediction of hydrogen embrittlement in iron’, Nat. Mater., vol. 12, no. 2, pp. 145–151, 2013
2013
-
[3]
Moore, Emilio Martínez-Pañeda, Andrej Atrens, and Julie M
Chen, Yi-Sheng, Chao Huang, Pang-Yu Liu, Hung-Wei Yen, Ranming Niu, Patrick Burr, Katie L. Moore, Emilio Martínez-Pañeda, Andrej Atrens, and Julie M. Cairney, ‘Hydrogen trapping and embrittlement in metals–A review’, Int. J. Hydrog. Energy, vol. 136, pp. 789–821, 2025
2025
-
[4]
Kappes, M
M. Kappes, M. Iannuzzi, and R. M. Carranz a, ‘Hydrogen embrittlement of magnesium and magnesium alloys: a review’, J. Electrochem. Soc., vol. 160, no. 4, pp. C168–C178, 2013
2013
-
[5]
Bloch and M
J. Bloch and M. H. Mintz, ‘Kinetics and mechanisms of metal hydrides formation—a review’, J. Alloys Compd., vol. 253, pp. 529–541, 1997
1997
-
[6]
Winzer, A
N. Winzer, A. Atrens, W. Dietzel, G. Song, and K. U. Kainer, ‘Evaluation of the delayed hydride cracking mechanism for transgranular stress corrosion cracking of magnesium alloys’, Mater. Sci. Eng. A, vol. 466, no. 1–2, pp. 18–31, 2007
2007
-
[7]
Jiang and E
D.-E. Jiang and E. A. Carter, ‘First principles assessment of ideal fracture energies of materials with mobile impurities: implications for hydrogen embrittlement of metals’, Acta Mater., vol. 52, no. 16, pp. 4801–4807, 2004
2004
-
[8]
A. Van Der Ven and G. Ceder, ‘Impurity-induced van der Waals transition during decohesion’, Phys. Rev. B, vol. 67, no. 6, p. 060101, 2003, doi: 10.1103/PhysRevB.67.060101
-
[9]
Martínez -Pañeda, C
E. Martínez -Pañeda, C. F. Niordson, and R. P. Gangloff, ‘Strain gradient plasticity -based modeling of hydrogen environment assisted cracking’, Acta Mater., vol. 117, pp. 321–332, 2016
2016
-
[10]
P. Gong, J. Nutter, P. E. J. Rivera -Diaz-Del-Castillo, and W. M. Rainforth, ‘Hydrogen embrittlement through the formation of low -energy dislocation nanostructures i n nanoprecipitation-strengthened steels’, Sci. Adv. , v ol. 6, no. 46, p. eabb6152, 2020, doi: 10.1126/sciadv.abb6152
-
[11]
H. K. Birnbaum and P. Sofronis, ‘Hydrogen -enhanced localized plasticity —a mechanism for hydrogen-related fracture’, Mater. Sci. Eng. A, vol. 176, no. 1–2, pp. 191–202, 1994
1994
-
[12]
Mater., vol
Huang, Longchao, Dengke Chen, Degang Xie, Suzhi Li, Yin Zhang, Ting Zhu, Dierk Raabe, En Ma, Ju Li, and Zhiwei Shan, ‘Quantitative tests revealing hydrogen-enhanced dislocation motion in α-iron’, Nat. Mater., vol. 22, no. 6, pp. 710–716, 2023
2023
-
[13]
Nagumo and K
M. Nagumo and K. Takai, ‘The predominant role of strain -induced vacancies in hydrogen embrittlement of steels: Overview’, Acta Mater., vol. 165, pp. 722–733, 2019
2019
-
[14]
Wesson, Tokamaks, International Series of Monographs on Physics, vol
J. Wesson, Tokamaks, International Series of Monographs on Physics, vol. 149. Oxford university press, 2011
2011
-
[15]
Serra, G
E. Serra, G. Benamati, and O. V. Ogorodnikova, ‘Hydrogen isotopes transport parameters in fusion reactor materials’, J. Nucl. Mater., vol. 255, no. 2–3, pp. 105–115, 1998
1998
-
[16]
Chen, Yi -Sheng, Hongzhou Lu, Jiangtao Liang, Alexander Rosenthal, Hongwei Liu, Glenn Sneddon, Ingrid McCarroll, Zhengzhi Zhao, Wei Li, Aimin Guo, Julie M Cairney, ‘Observation of hydrogen trapping at dislocations, grain boundaries, and precipitates’, Science, vol. 367 , no. 6474, pp. 171–175, Jan. 2020, doi: 10.1126/science.aaz0122
-
[17]
, ‘Hydrogen-enhanced grain boundary vacancy stockpiling causes transgranular to intergranular fracture transition’, Acta Mater
Ding, Yu, Haiyang Yu, Meichao Lin, Kai Zhao, Senbo Xiao, Alexey Vinogradov, Lijie Qiao, Michael Ortiz, Jianying He, and Zhiliang Zhang. , ‘Hydrogen-enhanced grain boundary vacancy stockpiling causes transgranular to intergranular fracture transition’, Acta Mater. , vol. 239, p. 118279, 2022
2022
-
[18]
Y. He, Y. Su, H. Yu, and C. Chen, ‘First -principles study of hydrogen trapping and diffusion at grain boundaries in γ-Fe’, Int. J. Hydrog. Energy, vol. 46, no. 10, pp. 7589–7600, 2021
2021
-
[19]
A. Boda, S. Bajania, S. M. Ali, K. T. Shenoy, and S. Mohan, ‘Chemisorption, diffusion and permeation of hydrog en isotopes in bcc bulk Cr and C r (100) surface: First -principles dft simulations’, J. Nucl. Mater., vol. 543, p. 152538, 2021
2021
-
[20]
T. C. Cui, P. F. Liu, and C. H. Gu, ‘Finite element analysis of hydrogen diffusion/plasticity coupled behaviors of low-alloy ferritic steel at large strain’, Int. J. Hydrog. Energy, vol. 42, no. 31, pp. 20324–20335, 2017
2017
-
[21]
A. Khosravi, J. Song, and N. Mousseau, ‘Kinetics of hydrogen and vacancy diffusion in iron: A kinetic activation relaxation technique (k -ART) study’, Phys. Rev. Mater. , vol. 7, no. 12, p. 123605, Dec. 2023, doi: 10.1103/PhysRevMaterials.7.123605
-
[22]
Kaur, Sanjeet, Daniel R. Mason, Prashanth Srinivasan, Stephen Dixon, Sid Mungale, Teresa Orr, Mikhail Yu Lavrentiev, and Duc Nguyen-Manh, ‘Effect of multi-occupancy traps on the diffusion and retention of multiple hydrogen isotopes in irradiated tungsten an d vanadium’, Phys. Rev. Mater., vol. 9, no. 12, p. 125404, 2025, doi: 10.1103/nbwm-bs8m
-
[23]
Wu, Meidie, Yiliang Liu, Luyao Hai, Xuan Liu, Siqi Jiang, Yawen Hua, and Shuoxue Jin, ‘Effect of vacancy and chromium on H behavior in Fe and Fe9Cr: a DFT study’, Int. J. Hydrog. Energy, vol. 188, p. 152066, 2025
2025
-
[24]
G. Kresse and J. Hafner, ‘Ab initio molecular dynamics for liquid metals’, Phys. Rev. B, vol. 47, no. 1, pp. 558–561, 1993, doi: 10.1103/PhysRevB.47.558
-
[25]
G. Kresse and J. Hafner, ‘ Ab initio molecular-dynamics simulation of the liquid -metal– amorphous-semiconductor transition in germanium’, Phys. Rev. B , vol. 4 9, no. 20, pp. 14251 – 14269,1994, doi: 10.1103/PhysRevB.49.14251
-
[26]
G. Kresse and J. Furthmüller, ‘Efficient iterative schemes for ab initio total-energy calculations using a plane -wave basis set’, Phys. Rev. B , vol. 5 4, no. 16, pp. 11169 –11186, 1996, doi: 10.1103/PhysRevB.54.11169
-
[27]
J. P. Perdew, K. Burke, and M. Ernzerhof, ‘Generalized Gradient Approximation Made Simple’, Phys. Rev. Lett., vol. 77, no. 18, pp. 3865–3868, 1996, doi: 10.1103/PhysRevLett.77.3865
-
[28]
G. Kresse and D. Joubert, ‘From ultrasoft pseudopotentials to the projector augmented -wave method’, Phys. Rev. B, vol. 59, no. 3, pp. 1758–1775, Jan. 1999, doi: 10.1103/PhysRevB.59.1758
-
[29]
H. J. Monkhorst and J. D. Pack, ‘Special points for Brillouin-zone integrations’, Phys. Rev. B, vol. 13, no. 12, pp. 5188–5192, Jun. 1976, doi: 10.1103/PhysRevB.13.5188
-
[30]
Henkelman, B
G. Henkelman, B. P. Uberuaga, and H. Jónsson, ‘A climbing image nudged elastic band method for finding saddle points and minimum energy paths’, J. Chem. Phys., vol. 113, no. 22, pp. 9901– 9904, 2000
2000
-
[31]
D. F. Johnson and E. A. Carter, ‘Hydrogen in tungsten: Absorption, diffusion, vacancy trapping, and decohesion’, J. Mater. Res., vol. 25, no. 2, pp. 315–327, 2010
2010
-
[32]
Wang, Austin J
Shang, Shun-Li, Bi-Cheng Zhou, William Y. Wang, Austin J. Ross, Xuan L. Liu, Yong -Jie Hu, Hua-Zhi Fang, Yi Wang, and Zi -Kui Liu, ‘A comprehensive first -principles study of pure elements: Vacancy formation and migration energies and self-diffusion coefficients’, Acta Mater., vol. 109, pp. 128–141, 2016
2016
-
[33]
K. Ohsawa, K. Eguchi, H. Watanabe, M. Yamaguchi, and M. Yagi, ‘Configuration and binding energy of multiple hydrogen atoms trapped in monovacancy in bcc transition metals’, Phys. Rev. B, vol. 85, no. 9, p. 094102, 2012, doi: 10.1103/PhysRevB.85.094102
-
[34]
Xing, X.-Q
W. Xing, X.-Q. Chen, Q. Xie, G. Lu, D. Li, and Y. Li, ‘Unified mechanism for hydrogen trapping at metal vacancies’, Int. J. Hydrog. Energy, vol. 39, no. 21, pp. 11321–11327, 2014
2014
-
[35]
Y. He, Y. Li, C. Chen, and H. Yu, ‘Diffusion coefficient of hydrogen interstitial atom in α-Fe, γ- Fe and ε-Fe crystals by first -principle calculations’, Int. J. Hydrog. Energy, vol. 42, no. 44, pp. 27438–27445, 2017
2017
-
[36]
D. A. Mirzaev, A. A. Mirzoev, K. Yu. Okishev, and A. V. Verkhovykh, ‘Hydrogen –vacancy interaction in bcc iron: ab initio calculations and thermodynamics’, Mol. Phys., vol. 112, no. 13, pp. 1745–1754, Jul. 2014, doi: 10.1080/00268976.2013.861087
-
[37]
Matsumoto, S
R. Matsumoto, S. Taketomi, and N. Miyazaki, ‘Atomistic study of hydrogen distributions around lattice defects and defect energies under hydrogen environment’, presented at the 12th International C onference on Fracture 2009, 2009, p. 6109. [Online]. Availabl e: https://www.proceedings.com/content/007/007799webtoc.pdf
2009
-
[38]
K. Ohsawa, J. Goto, M. Yamakami, M. Yamaguchi, and M. Yagi, ‘Trapping of multiple hydrogen atoms in a tungsten monovacancy from first principles’, Phys. Rev. B, vol. 82, no. 18, p. 184117, 2010, doi: 10.1103/PhysRevB.82.184117
-
[39]
Jia, Dian, William Yi Wang, Peixuan Li, Ying Zhang, Weijia Gong, Biao Chen, Jun Wang, and Jinshan Li, ‘The hydrogen-helium-vacancy interaction and hydrogen-vacancy clusters formation mechanisms in chromium: A first-principles study’, J. Nucl. Mater., vol. 601, p. 155338, 2024
2024
-
[40]
A. J. Samin, D. A. Andersson, E. F. Holby, and B. P. Uberuaga, ‘First-principles localized cluster expansion study of the kinetics of hydrogen diffusion in homogeneous and heterogeneous Fe-Cr alloys’, Phys. Rev. B, vol. 99, no. 1, p. 014110, 2019, doi: 10.1103/PhysRevB.99.014110
-
[41]
G. Lv, M. Zhang, H. Zhang, and Y. Su, ‘Hydrogen diffusion and vacancy clusterization in iron’, Int. J. Hydrog. Energy, vol. 43, no. 32, pp. 15378–15385, 2018
2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.