pith. sign in

arxiv: 2605.27983 · v1 · pith:NLOPUVORnew · submitted 2026-05-27 · ❄️ cond-mat.mtrl-sci

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

classification ❄️ cond-mat.mtrl-sci
keywords hydrogen trappingvacanciesBCC ironBCC chromiumDFTkinetic Monte Carlohydrogen diffusionhydrogen embrittlement
0
0 comments X

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.

The study combines density functional theory calculations with kinetic Monte Carlo simulations to track hydrogen atoms around vacancies in body-centered cubic iron and chromium. DFT results show hydrogen binds more tightly to vacancies in chromium than in iron, and migration barriers are extracted for different occupancies. Kinetic Monte Carlo runs then demonstrate that these vacancies slow long-range hydrogen movement and increase the effective activation energy, with the slowdown stronger in chromium. The work focuses on the mechanisms that govern how defects control hydrogen transport in these metals.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 1 minor

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)
  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

0 responses · 0 unresolved

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

0 steps flagged

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

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard assumptions of DFT accuracy for defect energetics and the validity of kMC rate models derived from those energies. No new entities are postulated.

axioms (2)
  • domain assumption DFT with chosen functional and pseudopotentials yields binding and migration energies accurate enough for kMC timescales
    Invoked implicitly when using DFT results as input for kMC diffusion simulations
  • domain assumption Bader charge analysis correctly partitions electron density to quantify trapping strength
    Used to support stronger trapping claim in Cr vs Fe

pith-pipeline@v0.9.1-grok · 5735 in / 1321 out tokens · 20028 ms · 2026-06-29T11:39:54.939196+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

41 extracted references · 16 canonical work pages

  1. [1]

    Novak, D

    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. [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

  3. [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

  4. [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

  5. [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

  6. [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

  7. [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

  8. [8]

    Van Der Ven and G

    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. [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

  10. [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. [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

  12. [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

  13. [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

  14. [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

  15. [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

  16. [16]

    367 , no

    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. [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

  18. [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

  19. [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

  20. [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

  21. [21]

    Khosravi, J

    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. [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. [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

  24. [24]

    Kresse and J

    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. [25]

    Kresse and J

    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. [26]

    Kresse and J

    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. [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. [28]

    Kresse and D

    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. [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. [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

  31. [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

  32. [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

  33. [33]

    Ohsawa, K

    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. [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

  35. [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

  36. [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. [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

  38. [38]

    Ohsawa, J

    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. [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

  40. [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. [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