pith. sign in

arxiv: 2606.26782 · v1 · pith:PESJHW3Pnew · submitted 2026-06-25 · ❄️ cond-mat.mtrl-sci

Hydrogen segregation around a straight screw dislocation in bcc iron

Margot Lucas , Marie Landeiro Dos Reis , Sylvain Queyreau , Xavier Feaugas This is my paper

Pith reviewed 2026-06-26 03:45 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords hydrogen segregationscrew dislocationbcc ironneural network potentialhydrogen embrittlementeasy corethermodynamic modelsolubility limit
0
0 comments X

The pith

Easy-core screw dislocations in bcc iron dominate hydrogen trapping and explain measured solubility limits.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper maps hydrogen binding energies to straight screw dislocations in body-centered cubic iron across easy and hard core structures, many atomic sites, and varying concentrations using molecular simulations driven by a neural network interatomic potential. From these energies the authors build a thermodynamic model that ranks the statistical importance of each trapping configuration and thereby collapses the high-dimensional interaction space into a much smaller set of relevant sites. This ranking shows that the easy-core geometry accounts for the experimentally observed hydrogen solubility limits far better than models relying only on hard cores or elastic approximations alone. The same framework also identifies the range of conditions where a simpler elastic-dipole description remains valid and supplies binding data suitable for insertion into larger plasticity models.

Core claim

Simulations with a neural network potential establish that hydrogen atoms preferentially segregate to sites around the easy-core configuration of a straight screw dislocation in bcc iron; when these site energies are inserted into a thermodynamic weighting scheme the resulting effective solubility matches experimental limits, whereas hard-core or purely elastic models do not.

What carries the argument

Thermodynamic framework that assigns statistical weights to hydrogen trapping sites on the basis of binding energies computed for easy-core and hard-core dislocation structures.

If this is right

  • Hydrogen solubility in bcc iron is controlled primarily by easy-core trapping rather than hard-core or bulk interstitial sites.
  • An elastic-dipole model of hydrogen-dislocation interaction is quantitatively reliable only inside a limited range of distances and concentrations.
  • Binding-energy tables from the simulations can be used directly in continuum or discrete-dislocation models of hydrogen-assisted plasticity.
  • The same thermodynamic ranking procedure can be applied to other interstitial species or to curved or mixed-character dislocations.

Where Pith is reading between the lines

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

  • If the easy-core preference holds under applied stress, hydrogen segregation may alter the Peierls barrier and therefore the mobility of screw dislocations at low temperature.
  • Extending the thermodynamic framework to finite hydrogen concentrations could predict the onset of hydride-like atmospheres or pipe diffusion along dislocation lines.
  • The same site-ranking approach might be tested on edge dislocations or on other bcc metals where experimental solubility data also exist.

Load-bearing premise

The neural network potential produces binding energies that remain accurate for all explored hydrogen-dislocation configurations even though direct density-functional-theory checks exist for only a small subset of those configurations.

What would settle it

An atom-probe tomography or similar measurement that finds hydrogen concentrations around screw dislocations in bcc iron to be independent of whether the core adopts the easy or hard configuration would falsify the claim that the easy core is required to explain solubility limits.

Figures

Figures reproduced from arXiv: 2606.26782 by Margot Lucas, Marie Landeiro Dos Reis, Sylvain Queyreau, Xavier Feaugas.

Figure 1
Figure 1. Figure 1: FIG. 1. Simulation cell with a screw dislocation (in green) [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Interaction energy between a hydrogen atom and a [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Interaction energy between a hydrogen atom and a [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Comparison of the interaction energy in the easy-core [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figure 4
Figure 4. Figure 4: The second-order approximation slightly im [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: ); and (ii) a configuration where the whole line is taken to be in the hard-core structure, giving L = nHC b (see [PITH_FULL_IMAGE:figures/full_fig_p005_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Saturation ratio of the [PITH_FULL_IMAGE:figures/full_fig_p006_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Saturation ratio of the E0, H0, and H1 sites as a [PITH_FULL_IMAGE:figures/full_fig_p006_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Interaction energy between a hydrogen atom, at the [PITH_FULL_IMAGE:figures/full_fig_p007_9.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Correlation between H solubility measure experi [PITH_FULL_IMAGE:figures/full_fig_p008_11.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Interaction energy computed within the frame of [PITH_FULL_IMAGE:figures/full_fig_p008_10.png] view at source ↗
read the original abstract

The interaction between hydrogen and screw dislocations in bcc iron is central to understanding hydrogen embrittlement. A major challenge lies in the high-dimensional parametric landscape governing this interaction. In this work, we perform a comprehensive set of molecular simulations using a reliable neural network interatomic potential, systematically exploring hydrogen binding across dislocation core structures (easy and hard cores), site types, and concentrations. From these energetics, we construct a thermodynamic framework that quantifies the statistical relevance of the various trapping configurations, thereby significantly reducing the complexity of the problem. Our results show good agreement with the limited density functional theory data available in the literature. We further delineate the validity domain of an elastic dipole description of hydrogen-dislocation interactions, providing a simplified yet physically grounded modeling approach. Finally, we demonstrate that the easy-core configuration plays a key role in rationalizing experimental hydrogen solubility limits. These findings establish a consistent multiscale foundation for incorporating hydrogen-dislocation interactions into larger-scale models of plasticity and embrittlement.

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

2 major / 2 minor

Summary. The manuscript uses molecular simulations with a neural network interatomic potential (NNIP) to systematically map hydrogen binding energetics around straight screw dislocations in bcc Fe, distinguishing easy- and hard-core structures across site types and concentrations. From these energies a thermodynamic framework is constructed to quantify the statistical weight of different trapping configurations, reducing the parametric complexity. Results are compared to the limited available DFT data, the validity range of an elastic-dipole model is delineated, and the easy-core configuration is shown to be essential for rationalizing experimental hydrogen solubility limits.

Significance. If the NNIP binding energies are sufficiently accurate, the work supplies a concrete multiscale bridge from atomistic energetics to macroscopic solubility, offering a physically grounded simplification for larger-scale models of hydrogen embrittlement. The systematic exploration of core structures and concentrations together with the thermodynamic reduction constitute clear strengths.

major comments (2)
  1. [Abstract / Results (easy-core)] Abstract and the results section on easy-core configurations: the claim that the easy-core plays a key role in rationalizing experimental solubility limits rests on the statistical weighting derived from NNIP binding energies; the manuscript reports agreement only with the limited existing DFT data and does not supply extensive new DFT benchmarks across the systematically explored site types, concentrations, and core structures. A systematic bias of even 20–50 meV would alter the relative weights and weaken the link to solubility data.
  2. [Thermodynamic framework] Thermodynamic framework section: the framework is built directly from the NNIP-derived energies; it is unclear whether any implicit fitting or selection of data was performed to achieve agreement with solubility limits, which would affect the independence asserted in the abstract.
minor comments (2)
  1. [Figures] Figures showing binding-energy maps should include explicit error bars or convergence checks from the NNIP simulations.
  2. [Methods / Thermodynamic framework] All symbols appearing in the thermodynamic equations (e.g., site occupancies, partition functions) should be defined in a single notation table or paragraph.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and valuable comments. We address each major point below and have revised the manuscript to improve clarity and address concerns where possible.

read point-by-point responses

  1. Referee: [Abstract / Results (easy-core)] Abstract and the results section on easy-core configurations: the claim that the easy-core plays a key role in rationalizing experimental solubility limits rests on the statistical weighting derived from NNIP binding energies; the manuscript reports agreement only with the limited existing DFT data and does not supply extensive new DFT benchmarks across the systematically explored site types, concentrations, and core structures. A systematic bias of even 20–50 meV would alter the relative weights and weaken the link to solubility data.

    Authors: We acknowledge that the validation relies on the limited existing DFT data and that the NNIP, while extensively benchmarked on hydrogen-iron configurations with errors typically below 15 meV/H, does not include new DFT calculations for every site and concentration explored. Performing such benchmarks would be computationally intensive and outside the primary scope of mapping NNIP energetics. However, the NNIP training set includes dislocation-relevant environments, and trends match available DFT. To address potential bias concerns, we have added a sensitivity analysis in the revised results section showing that a uniform 30 meV shift does not change the dominance of easy-core sites in the solubility model. This supports the robustness of the conclusions. revision: partial

  2. Referee: [Thermodynamic framework] Thermodynamic framework section: the framework is built directly from the NNIP-derived energies; it is unclear whether any implicit fitting or selection of data was performed to achieve agreement with solubility limits, which would affect the independence asserted in the abstract.

    Authors: No fitting, parameter adjustment, or selective data inclusion was performed to match solubility limits. The thermodynamic framework uses the full set of NNIP binding energies directly, with statistical weights computed without any post-hoc tuning. The agreement with experimental solubility is an emergent outcome. We have revised the thermodynamic framework section and abstract to explicitly state that the model contains no parameters fitted to solubility data, thereby clarifying the independence of the prediction. revision: yes

Circularity Check

0 steps flagged

No significant circularity: derivation relies on independent NNIP simulations and external DFT benchmarks

full rationale

The paper derives binding energetics from molecular simulations with a neural network interatomic potential, constructs a thermodynamic framework directly from those computed energies, validates against independent literature DFT data, and applies the framework to interpret experimental solubility limits. No step reduces by construction to a fitted input, self-definition, or self-citation chain; the central claim follows from external simulation outputs without renaming or smuggling ansatzes.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Central claim rests on accuracy of the NN interatomic potential (domain assumption) and validity of the thermodynamic framework built from simulation energetics. No explicit free parameters or invented entities stated in abstract.

axioms (1)
  • domain assumption Neural network interatomic potential accurately models H-Fe interactions near dislocations
    All simulations and derived framework rely on this; validated only against limited DFT.

pith-pipeline@v0.9.1-grok · 5708 in / 1090 out tokens · 64468 ms · 2026-06-26T03:45:28.251841+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

47 extracted references · 2 canonical work pages

  1. [1]

    R. A. Oriani, A mechanistic theory of hydrogen em- brittlement of steels, Berichte der Bunsengesellschaft für physikalische Chemie 76, 848 (1972)

    1972

  2. [2]

    A. R. Troiano, The Role of Hydrogen and Other Intersti- tials in the Mechanical Behavior of Metals: (1959 Edward De Mille Campbell Memorial Lecture), Metallography, Microstructure, and Analysis 5, 557 (2016)

    1959

  3. [3]

    Ehlers, M

    F. Ehlers, M. Seydou, D. Tingaud, F. Maurel, Y. Charles, and S. Queyreau, H induced decohesion of an al grain boundary investigated with first principles: General con- ditions for instant breakage and local delayed fracture, Computational Materials Science 173, 109403 (2020)

    2020

  4. [4]

    embrittlement

    C. D. Beachem, A new model for hydrogen-assisted cracking (hydrogen “embrittlement”), Metallurgical and Materials Transactions B 3, 441 (1972)

    1972

  5. [5]

    Birnbaum and P

    H. Birnbaum and P. Sofronis, Hydrogen-enhanced local- ized plasticity—a mechanism for hydrogen-related frac- ture, Materials Science and Engineering: A 176, 191 (1994)

    1994

  6. [6]

    Lynch, Hydrogen embrittlement phenomena and mechanisms, Corrosion Reviews 30, 10.1515/corrrev- 2012-0502 (2012)

    S. Lynch, Hydrogen embrittlement phenomena and mechanisms, Corrosion Reviews 30, 10.1515/corrrev- 2012-0502 (2012)

    work page doi:10.1515/corrrev- 2012

  7. [7]

    Feaugas and C

    X. Feaugas and C. Scott, The relationship between me- chanical strength and hydrogen embrittlement, New Ad- vanced High Strength Steels: Optimizing Properties , 263 (2023)

    2023

  8. [8]

    Nagumo et al

    M. Nagumo et al. , Fundamentals of hydrogen embrittle- ment, Vol. 921 (Springer, 2016)

    2016

  9. [9]

    I. M. Robertson, P. Sofronis, A. Nagao, M. L. Martin, S. Wang, D. W. Gross, and K. E. Nygren, Hydrogen Embrittlement Understood, Metallurgical and Material Transaction B , 20 (2015)

    2015

  10. [10]

    Barrera, D

    O. Barrera, D. Bombac, Y. Chen, T. D. Daff, E. Galindo- Nava, P. Gong, D. Haley, R. Horton, I. Katzarov, J. R. Kermode, C. Liverani, M. Stopher, and F. Sweeney, Un- derstanding and mitigating hydrogen embrittlement of steels: a review of experimental, modelling and design progress from atomistic to continuum, Journal of Mate- rials Science 53, 6251 (2018)

    2018

  11. [11]

    H. Yu, A. Díaz, X. Lu, B. Sun, Y. Ding, M. Koyama, J. He, X. Zhou, A. Oudriss, X. Feaugas, and Z. Zhang, Hydrogen embrittlement as a conspicuous material chal- lenge￿comprehensive review and future directions, Chem- ical Reviews 124, 6271 (2024)

    2024

  12. [12]

    S. L. Frederiksen and K. W. Jacobsen, Density functional theory studies of screw dislocation core structures in bcc metals, Philosophical Magazine 83, 365 (2003)

    2003

  13. [14]

    Dezerald, L

    L. Dezerald, L. Ventelon, E. Clouet, C. Denoual, D. Rod- ney, and F. Willaime, Ab initio modeling of the two- dimensional energy landscape of screw dislocations in bcc transition metals, Physical Review B 89, 024104 (2014)

    2014

  14. [15]

    Clouet, B

    E. Clouet, B. Bienvenu, L. Dezerald, and D. Rod- ney, Screw dislocations in BCC transition metals: from ab initio modeling to yield criterion, Comptes Rendus. Physique 22, 83 (2021)

    2021

  15. [16]

    M. R. Gilbert, S. Queyreau, and J. Marian, Stress and temperature dependence of screw dislocation mobility in α-Fe by molecular dynamics, Physical Review B 84, 174103 (2011)

    2011

  16. [17]

    M. R. Gilbert, P. Schuck, B. Sadigh, and J. Mar- ian, Free Energy Generalization of the Peierls Poten- tial in Iron, Physical Review Letters 111, 10.1103/Phys- RevLett.111.095502 (2013)

    work page doi:10.1103/phys- 2013

  17. [18]

    Rodney and L

    D. Rodney and L. Proville, Stress-dependent Peierls po- tential: Influence on kink-pair activation, Physical Re- view B 79, 094108 (2009)

    2009

  18. [19]

    Freitas, M

    R. Freitas, M. Asta, and V. V. Bulatov, Quantum ef- fects on dislocation motion from ring-polymer molecular dynamics, npj Computational Materials 4, 55 (2018)

    2018

  19. [20]

    Itakura, H

    M. Itakura, H. Kaburaki, and M. Yamaguchi, First- principles study on the mobility of screw dislocations in bcc iron, Acta Materialia 60, 3698 (2012). 10

    2012

  20. [21]

    M. e. a. Itakura, The effect of hydrogen atoms on the screw dislocation mobility in bcc iron: A first-principles study, Acta Materialia 61, 6857 (2013)

    2013

  21. [22]

    E. L. Simpson and A. Paxton, Effect of applied strain on the interaction between hydrogen atoms and 1/2<111> screw dislocations in α-iron, International Journal of Hy- drogen Energy 45, 20069 (2020)

    2020

  22. [23]

    L. e. a. Huang, Quantitative tests revealing hydrogen- enhanced dislocation motion in α-iron, Nature Materials 22, 710 (2023)

    2023

  23. [24]

    Zheng, P

    Y. Zheng, P. Yu, and L. Zhang, Atomistic study of the in- teraction nature of h-dislocation and the validity of elas- ticity in bcc iron, Metals 13, 1267 (2023)

    2023

  24. [25]

    Moro, Fragilisation par l’hydrogène gazeux d’un acier ferrito-perlitique de grade API X80 , Ph.D

    I. Moro, Fragilisation par l’hydrogène gazeux d’un acier ferrito-perlitique de grade API X80 , Ph.D. thesis, Institut National Polytechnique de Toulouse-INPT (2009)

    2009

  25. [26]

    Meng and J.-P

    F.-S. Meng and J.-P. e. a. Du, General-purpose neural network interatomic potential for the α-iron and hydro- gen binary system, Physical Review Materials 5, 113606 (2021)

    2021

  26. [27]

    M. J. Gillan, The elastic dipole tensor for point defects in ionic crystals, Journal of Physics C: Solid State Physics 17, 1473 (1984)

    1984

  27. [28]

    R. e. a. Nazarov, First-principles calculation of the elastic dipole tensor of a point defect, Physical Review B 94, 241112 (2016)

    2016

  28. [29]

    Clouet, C

    E. Clouet, C. Varvenne, and T. Jourdan, Elastic model- ing of point-defects and their interaction, Computational Materials Science 147, 49 (2018)

    2018

  29. [30]

    A. P. Thompson, H. M. Aktulga, and R. e. a. Berger, LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales, Computer Physics Communications 271, 108171 (2022)

    2022

  30. [31]

    Lee and J.-W

    B.-J. Lee and J.-W. Jang, A modified embedded-atom method interatomic potential for the fe–h system, Acta Materialia 55, 6779 (2007)

    2007

  31. [32]

    Wen, A new interatomic potential describing fe-h and hh interactions in bcc iron, Computational Materials Sci- ence 197, 110640 (2021)

    M. Wen, A new interatomic potential describing fe-h and hh interactions in bcc iron, Computational Materials Sci- ence 197, 110640 (2021)

    2021

  32. [33]

    P. P. P. O. Borges, E. Clouet, and L. Ventelon, Ab ini- tio investigation of the screw dislocation-hydrogen inter- action in bcc tungsten and iron, Acta Materialia 234, 118048 (2022)

    2022

  33. [34]

    Luthi, Modelisation ab initio des interactions dislocation-solute dans les metaux cubiques centres, Ph.D

    B. Luthi, Modelisation ab initio des interactions dislocation-solute dans les metaux cubiques centres, Ph.D. thesis, Universite de Lyon (2017)

    2017

  34. [35]

    Luthi, L

    B. Luthi, L. Ventelon, D. Rodney, and F. Willaime, At- tractive interaction between interstitial solutes and screw dislocations in bcc iron from first principles, Computa- tional Materials Science 148, 21 (2018)

    2018

  35. [36]

    Guén̈olé and W

    J. Guén̈olé and W. G. e. a. Nöhring, Assessment and op- timization of the fast inertial relaxation engine (fire) for energy minimization in atomistic simulations and its im- plementation in lammps, Computational Materials Sci- ence 175, 109584 (2020)

    2020

  36. [37]

    A. Stukowski, Visualization and analysis of atomistic simulation data with OVITO–the Open Visualization Tool, Modelling and Simulation in Materials Science and Engineering 18, 015012 (2009)

    2009

  37. [38]

    Zhao and G

    Y. Zhao and G. Lu, Qm/mm study of dislocation— hydrogen/helium interactions in α-fe, Modelling and Simulation in Materials Science and Engineering 19, 065004 (2011)

    2011

  38. [39]

    J. H. Andrew, H. Lee, Bhat, U V, and Lloyd, H K, The Effect of Cold-Work on Steel, JOURNAL OF THE IRON AND STEEL INSTITUTE , 382 (1950)

    1950

  39. [40]

    L. S. Darken and R. P. Smith, Behavior of Hydrogen in Steel During and After Immersion in Acid, Corrosion 5, 1 (1949)

    1949

  40. [41]

    Tabor, The hardness of metals (Oxford university press, 2000)

    D. Tabor, The hardness of metals (Oxford university press, 2000)

    2000

  41. [42]

    Queyreau, G

    S. Queyreau, G. Monnet, and B. Devincre, Slip systems interactions in α-iron determined by dislocation dynam- ics simulations, International Journal of Plasticity 25, 361 (2009)

    2009

  42. [43]

    Kahloun, G

    C. Kahloun, G. Monnet, S. Queyreau, L. Le, and P. Fran- ciosi, A comparison of collective dislocation motion from single slip quantitative topographic analysis during in- situ afm room temperature tensile tests on cu and α- fe crystals, International Journal of Plasticity 84, 277 (2016)

    2016

  43. [44]

    Monnet, L

    G. Monnet, L. Vincent, and B. Devincre, Dislocation- dynamics based crystal plasticity law for the low- and high-temperature deformation regimes of bcc crystal, Acta Materialia 61, 6178 (2013)

    2013

  44. [45]

    Queyreau, Dislocation Based Mechanics: the various contributions of Dislocation Dynamics simulations (John Wiley & Sons, 2024) p

    S. Queyreau, Dislocation Based Mechanics: the various contributions of Dislocation Dynamics simulations (John Wiley & Sons, 2024) p. 1

    2024

  45. [46]

    Takeuchi, Work Hardening of Copper Single Crys- tals with Multiple Glide Orientations, Transactions of the Japan Institute of Metals 16, 629 (1975)

    T. Takeuchi, Work Hardening of Copper Single Crys- tals with Multiple Glide Orientations, Transactions of the Japan Institute of Metals 16, 629 (1975)

    1975

  46. [47]

    F. B. Pickering, Physical metallurgy and the design of steels (1978)

    1978

  47. [48]

    Spitzig and A

    W. Spitzig and A. Keh, The effect of orientation and temperature on the plastic flow properties of iron single crystals, Acta Metallurgica 18, 611 (1970)

    1970