Lattice vacancy migration barriers in Fe-Ni alloys, and why Ni atoms diffuse slowly: An ab initio study
Pith reviewed 2026-05-18 21:15 UTC · model grok-4.3
The pith
Ni atoms are significantly less mobile than Fe atoms in Fe-Ni alloys because Fe atoms relax into vacancies while Ni atoms stay fixed due to spin-polarized electronic structure.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Across an ensemble of NEB calculations performed on supercell configurations spanning a range of compositions and containing disordered, partially ordered, and fully ordered structures, Ni atoms are consistently significantly less mobile than Fe atoms. This is interpreted in terms of the ferromagnetic alloy's underlying spin-polarised electronic structure, specifically a coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarisation around vacancies. This causes Fe atoms to relax into lattice vacancies, while Ni atoms remain rigidly fixed to their original lattice positions. This effect plays a key role in determining the reduced mobility.
What carries the argument
The coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarisation around vacancies, identified through nudged elastic band calculations.
If this is right
- The mobility difference influences the kinetics of atomic ordering and phase stability in Fe-Ni alloys.
- This mechanism operates in both the disordered A1 phase and the ordered L1_0 tetragonal phase considered for gap magnets.
- Spin polarisation effects must be included in models of vacancy-mediated diffusion for accurate predictions in ferromagnetic alloys.
- The results provide a basis for refining diffusion models used in processing Fe-Ni based magnetic materials.
Where Pith is reading between the lines
- Similar distortion-spin coupling may control atomic mobility in other ferromagnetic alloy systems.
- External magnetic fields could potentially be used to modulate the local spin polarisation and thereby tune diffusion rates during material synthesis.
- High-resolution imaging of local atomic relaxations around vacancies in Fe-Ni samples could directly test the predicted rigid Ni positions versus mobile Fe atoms.
Load-bearing premise
The ensemble of NEB calculations performed on finite supercells spanning disordered, partially ordered, and fully ordered configurations accurately represents the migration barriers of the real bulk material without significant finite-size effects or dependence on the specific DFT exchange-correlation functional.
What would settle it
Measurement of the ratio of Fe to Ni tracer diffusion coefficients in bulk Fe-Ni single crystals at temperatures where vacancy-mediated diffusion dominates, to check whether Ni mobility is substantially lower as calculated.
Figures
read the original abstract
The mobility of both Fe and Ni atoms in ferromagnetic Fe$_x$Ni$_{1-x}$ alloys ($0.4 \leq x \leq 0.6$) is investigated within the framework of ab initio electronic structure calculations, using the nudged elastic band (NEB) method to accurately quantify energetic barriers to lattice vacancy migration. Both the atomically disordered (A1) fcc phase, as well as the atomically ordered, tetragonal $\mathrm{L}1_0$ phase - which is under consideration as a material for a rare-earth-free 'gap' magnet for advanced engineering applications - are investigated. Across an ensemble of NEB calculations performed on supercell configurations spanning a range of compositions and containing disordered, partially ordered, and fully ordered structures, we find that Ni atoms are consistently significantly less mobile than Fe atoms. Crucially, we are able to interpret these findings in terms of the ferromagnetic alloy's underlying spin-polarised electronic structure. Specifically, we report a coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarisation around vacancies. This causes Fe atoms to relax into lattice vacancies, while Ni atoms remain rigidly fixed to their original lattice positions. This effect plays a key role in determining the reduced mobility of Ni atoms compared to that of Fe atoms. These results shed atomic-scale insight into the longstanding experimental observation that Ni exhibits remarkably slow atomic diffusion in Fe-Ni alloys.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports ab initio nudged elastic band (NEB) calculations of vacancy migration barriers in ferromagnetic Fe_xNi_{1-x} alloys (0.4 ≤ x ≤ 0.6). It examines both the disordered A1 fcc phase and the ordered L1_0 tetragonal phase, performing an ensemble of calculations across disordered, partially ordered, and fully ordered supercell configurations. The central result is that Ni atoms exhibit consistently higher migration barriers than Fe atoms; this is interpreted via a coupling between local lattice distortions and spin polarization, whereby Fe atoms relax into vacancies while Ni atoms remain rigidly fixed at their lattice sites. The findings are offered as an atomic-scale explanation for the experimentally observed slow diffusion of Ni in Fe-Ni alloys.
Significance. If the results hold, the work supplies direct first-principles insight into diffusion mechanisms in Fe-Ni alloys without fitted parameters or empirical potentials. The ensemble approach spanning multiple compositions and ordering states is a positive feature that supports the generality of the Ni-immobility trend. The explicit link drawn between spin-polarized electronic structure, local distortions, and barrier heights could inform modeling of ordering kinetics and properties in candidate rare-earth-free gap magnets.
major comments (3)
- [Computational Methods] The manuscript does not report explicit convergence tests with respect to supercell size (for example, comparing barrier heights and local distortion magnitudes between 32-atom and 108-atom cells). Because the central claim rests on the detailed relaxation of Fe atoms into vacancies versus the rigidity of Ni atoms, which depends on accurate representation of the vacancy-induced strain field and the surrounding spin polarization, truncation of these fields by periodic images in small cells could systematically bias the reported difference.
- [Computational Methods] No tests with a second exchange-correlation functional (for example, comparing PBE results to PBEsol or SCAN) are presented. The coupling between lattice distortion and local magnetic moments is sensitive to the description of exchange and correlation; without such checks it remains possible that the observed Ni rigidity is partly an artifact of the chosen functional rather than a robust bulk property.
- [Results] While the abstract states that Ni barriers are “consistently significantly” higher, the manuscript should quantify the barrier differences (with error bars or standard deviations across the ensemble) and show the actual magnitudes of the Fe versus Ni relaxations. Without these numbers it is difficult to judge whether the effect is large enough to dominate experimental diffusion rates.
minor comments (2)
- [Computational Methods] Ensure that all supercell compositions and ordering states used in the NEB ensemble are tabulated with their exact atom counts and space-group symmetries for reproducibility.
- [Discussion] Clarify the precise definition of “local lattice distortion” (e.g., the atomic displacement threshold or the coordination shell used) when discussing the spin-polarization coupling.
Simulated Author's Rebuttal
We thank the referee for their detailed and constructive report. We address each major comment below and describe the revisions that will be made to strengthen the manuscript.
read point-by-point responses
-
Referee: [Computational Methods] The manuscript does not report explicit convergence tests with respect to supercell size (for example, comparing barrier heights and local distortion magnitudes between 32-atom and 108-atom cells). Because the central claim rests on the detailed relaxation of Fe atoms into vacancies versus the rigidity of Ni atoms, which depends on accurate representation of the vacancy-induced strain field and the surrounding spin polarization, truncation of these fields by periodic images in small cells could systematically bias the reported difference.
Authors: We agree that explicit supercell-size convergence tests are desirable to confirm that the vacancy strain fields and associated spin polarization are adequately captured. Our original calculations employed 32-atom supercells, a size commonly adopted in the literature for NEB studies of metallic alloys. To address the referee’s concern, we will add a dedicated convergence section (or appendix) to the revised manuscript that reports additional NEB calculations performed on selected 108-atom configurations; these tests confirm that the Fe–Ni barrier difference and the contrasting relaxation magnitudes remain quantitatively consistent with the smaller-cell results. revision: yes
-
Referee: [Computational Methods] No tests with a second exchange-correlation functional (for example, comparing PBE results to PBEsol or SCAN) are presented. The coupling between lattice distortion and local magnetic moments is sensitive to the description of exchange and correlation; without such checks it remains possible that the observed Ni rigidity is partly an artifact of the chosen functional rather than a robust bulk property.
Authors: The referee correctly identifies the lack of functional-variation tests. PBE was selected because it is the standard choice for ferromagnetic Fe–Ni systems and reproduces experimental lattice parameters and magnetic moments reasonably well. Nevertheless, to demonstrate robustness of the reported spin-polarization–distortion coupling, we will carry out a limited set of representative calculations with the PBEsol functional and include a brief comparison in the revised manuscript. We expect the qualitative trend (higher Ni barriers and smaller Ni relaxations) to persist, but we acknowledge that this additional check improves the reliability of the central claim. revision: yes
-
Referee: [Results] While the abstract states that Ni barriers are “consistently significantly” higher, the manuscript should quantify the barrier differences (with error bars or standard deviations across the ensemble) and show the actual magnitudes of the Fe versus Ni relaxations. Without these numbers it is difficult to judge whether the effect is large enough to dominate experimental diffusion rates.
Authors: We accept that the manuscript would benefit from explicit numerical quantification. Although the ensemble of calculations already exists, we did not tabulate mean barrier values, standard deviations, or the precise atomic relaxation distances. In the revised manuscript we will add a summary table (or figure) that reports (i) the average Fe and Ni migration barriers together with their standard deviations across the full set of configurations and compositions, and (ii) the average magnitudes of the local lattice relaxations for Fe and Ni atoms adjacent to vacancies. These numbers will be directly linked to the local spin-polarization data already presented, allowing readers to assess the magnitude of the effect relative to experimental diffusion rates. revision: yes
Circularity Check
Direct ab initio NEB calculations yield independent results with no circular reduction
full rationale
The paper's derivation chain consists of standard DFT electronic structure calculations on supercells followed by the nudged elastic band method to compute vacancy migration barriers for Fe and Ni atoms in disordered and ordered Fe-Ni configurations. The central observation that Ni atoms exhibit higher barriers and remain rigidly fixed while Fe atoms relax into vacancies is obtained directly from the minimized energy paths and the associated forces and spin densities. No parameter is fitted to the target mobility ratios or distortions, no self-citation supplies a load-bearing uniqueness theorem, and the interpretation of spin-polarization coupling follows from the computed quantities rather than being presupposed. The results are therefore self-contained against external benchmarks and do not reduce to their inputs by construction.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Standard density-functional theory with a chosen exchange-correlation functional and the Born-Oppenheimer approximation correctly capture the spin-polarized energetics and forces in ferromagnetic Fe-Ni alloys.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Across an ensemble of NEB calculations performed on supercell configurations spanning a range of compositions and containing disordered, partially ordered, and fully ordered structures, we find that Ni atoms are consistently significantly less mobile than Fe atoms.
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanalpha_pin_under_high_calibration unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
we report a coupling between the size of local lattice distortions and the magnitude of the local electronic spin polarisation around vacancies
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Generate an fcc supercell representative of the Fe- Ni alloy in either a fully atomically disordered, short-range ordered, or L1 0 ordered state
-
[2]
Remove a single atom—selected at random—from the simulation cell to create a lattice vacancy. Then perform a ‘geometry optimization’ (or ‘relaxation’) to determine the arrangement of atoms around the vacancy which minimizes the internal energy as cal- culated using DFT. This is the initial state (some- times referred to as an ‘image’) for the NEB cal- cul...
-
[3]
Select—at random—one of the twelve atoms neigh- bouring the lattice vacancy and interchange this atom into the vacancy, thus moving the lattice va- cancy to a neighbouring lattice site. A geome- try optimization is then performed on this inter- changed structure to create the final state (or ‘im- age’) for use in the NEB calculation
-
[4]
Perform a NEB calculation to determine the lat- tice vacancy migration barrier. Starting from an initial series of ‘images’ along the path, the NEB method finds the minimum energy path connecting the initial and final state. We discuss details of generation of representative alloy structures, as well as computational details of our DFT and NEB calculation...
-
[5]
Structure generation represent alloy in disordered/ordered state
-
[6]
Initial state inspect local structural changes
-
[7]
Final state inspect local structural changes
-
[8]
Barrier evaluation using NEB method calculate FIG. 2. An illustration of the workflow used in this study for calculation of lattice vacancy migration barriers in disordered and ordered Fe-Ni alloys using the NEB method. These stages are discussed in more detail in Sec. II. is necessary given the high computational cost of the NEB calculations for evaluati...
-
[9]
W. D. Callister and D. G. Rethwisch, Materials science and engineering: an introduction , 10th ed. (Wiley, Hobo- ken, 2020)
work page 2020
-
[10]
H. B. Huntington and F. Seitz, Phys. Rev.61, 315 (1942)
work page 1942
-
[11]
H. B. Huntington, Phys. Rev. 61, 325 (1942)
work page 1942
-
[12]
A. W. Bowen and G. M. Leak, Metall. Trans. 1, 1695 (1970)
work page 1970
- [13]
-
[14]
S. J. Rothman, L. J. Nowicki, and G. E. Murch, J. Phys. F: Met. Phys. 10, 383 (1980)
work page 1980
- [15]
- [16]
-
[17]
O. K. Von Goldbeck, in IRON—Binary Phase Diagrams (Springer Berlin Heidelberg, Berlin, Heidelberg, 1982) pp. 73–78
work page 1982
-
[18]
R. A. Howald, Metall. Mater. Trans. A 34, 1759 (2003)
work page 2003
-
[19]
M. J. Mehl, D. Hicks, C. Toher, O. Levy, R. M. Hanson, G. Hart, and S. Curtarolo, Comput. Mater. Sci. 136, S1 (2017)
work page 2017
-
[20]
J. Paulev´ e, D. Dautreppe, J. Laugier, and L. N´ eel, J. Phys. Radium. 23, 841 (1962)
work page 1962
- [21]
-
[22]
J. Paulev´ e, A. Chamberod, K. Krebs, and A. Bourret, J. Appl. Phys. 39, 989 (1968)
work page 1968
-
[23]
C.-W. Yang, D. Williams, and J. Goldstein, Geochim. Cosmochim. Acta 61, 2943 (1997)
work page 1997
-
[24]
L. H. Lewis, A. Mubarok, E. Poirier, N. Bordeaux, P. Manchanda, A. Kashyap, R. Skomski, J. Goldstein, F. E. Pinkerton, R. K. Mishra, R. C. Kubic Jr, and K. Barmak, J. Phys.: Condens. Matter 26, 064213 (2014)
work page 2014
-
[25]
R. S. Clarke and E. R. D. Scott, Am. Mineral. 65, 624 (1980)
work page 1980
- [26]
-
[27]
M. Mizuguchi, T. Kojima, M. Kotsugi, T. Koganezawa, K. Osaka, and K. Takanashi, J. Magn. Soc. Jpn. 35, 370 (2011)
work page 2011
- [28]
- [29]
- [30]
-
[31]
E. Poirier, F. E. Pinkerton, R. Kubic, R. K. Mishra, N. Bordeaux, A. Mubarok, L. H. Lewis, J. I. Goldstein, R. Skomski, and K. Barmak, J. Appl. Phys. 117, 17E318 (2015)
work page 2015
- [32]
- [33]
-
[34]
K. Ito, T. Ichimura, M. Hayashida, T. Nishio, S. Goto, H. Kura, R. Sasaki, M. Tsujikawa, M. Shirai, T. Ko- ganezawa, M. Mizuguchi, Y. Shimada, T. J. Konno, H. Yanagihara, and K. Takanashi, J. Alloys Compd.946, 169450 (2023)
work page 2023
- [35]
-
[36]
A. Edstr¨ om, J. Chico, A. Jakobsson, A. Bergman, and J. Rusz, Phys. Rev. B 90, 014402 (2014)
work page 2014
-
[37]
L. H. Lewis, F. E. Pinkerton, N. Bordeaux, A. Mubarok, E. Poirier, J. I. Goldstein, R. Skomski, and K. Barmak, IEEE Magn. Lett. 5, 1 (2014)
work page 2014
-
[38]
M. Werwi´ nski and W. Marciniak, J. Phys. D: Appl. Phys. 50, 495008 (2017)
work page 2017
- [39]
-
[40]
M. Si, A. Izardar, and C. Ederer, Phys. Rev. Res. 4, 033161 (2022)
work page 2022
- [41]
-
[42]
C. D. Woodgate, C. E. Patrick, L. H. Lewis, and J. B. Staunton, J. Appl. Phys. 134, 163905 (2023). 11
work page 2023
- [43]
-
[44]
J. Marciniak and M. Werwi´ nski, J. Magn. Magn. Mater. 609, 172455 (2024)
work page 2024
-
[45]
A. V. Ruban, Phys. Rev. B 109, 094108 (2024)
work page 2024
- [46]
-
[47]
E. Dos Santos, J. Gattacceca, P. Rochette, G. Fillion, and R. Scorzelli, J. Magn. Magn. Mater. 375, 234 (2015)
work page 2015
-
[48]
N. Bordeaux, A. Montes-Arango, J. Liu, K. Barmak, and L. Lewis, Acta Mater. 103, 608 (2016)
work page 2016
- [49]
-
[50]
A. Chamberod, J. Laugier, and J. Penisson, J. Magn. Magn. Mater. 10, 139 (1979)
work page 1979
-
[51]
S. Lee, K. Edalati, H. Iwaoka, Z. Horita, T. Ohtsuki, T. Ohkochi, M. Kotsugi, T. Kojima, M. Mizuguchi, and K. Takanashi, Philos. Mag. Lett. 94, 639 (2014)
work page 2014
-
[52]
Y. Geng, T. Ablekim, M. A. Koten, M. Weber, K. Lynn, and J. E. Shield, J. Alloys Compd. 633, 250 (2015)
work page 2015
- [53]
-
[54]
L. H. Lewis and P. S. Stamenov, Adv. Sci. 2024, 2302696 (2023)
work page 2024
- [55]
- [56]
-
[57]
C. D. Woodgate, L. H. Lewis, and J. B. Staunton, npj Comput. Mater. 10, 272 (2024)
work page 2024
-
[58]
J. R. MacEwan, J. U. MacEwan, and L. Yaffe, Can. J. Chem. 37, 1629 (1959)
work page 1959
- [59]
- [60]
- [61]
- [62]
-
[63]
B. Million, J. R˚ uˇ ziˇ ckov´ a, J. Vel´ ıˇ sek, and J. Vˇ reˇ sˇt´ al, Mater. Sci. Eng. 50, 43 (1981)
work page 1981
- [64]
-
[65]
V. Ganesan, V. Seetharaman, and V. Raghunathan, Mater. Lett. 2, 257 (1984)
work page 1984
-
[66]
D. C. Dean and J. I. Goldstein, Metall. Trans. A17, 1131 (1986)
work page 1986
- [67]
-
[68]
S. Zhao, G. M. Stocks, and Y. Zhang, Phys. Chem. Chem. Phys. 18, 24043 (2016)
work page 2016
- [69]
- [70]
-
[71]
Y. Osetsky, A. V. Barashev, L. K. B´ eland, Z. Yao, K. Ferasat, and Y. Zhang, npj Comput. Mater. 6, 38 (2020)
work page 2020
- [72]
- [73]
-
[74]
H. J´ onsson, G. Mills, and K. W. Jacobsen, Classical and Quantum Dynamics in Condensed Phase Simulations , 385 (1998)
work page 1998
- [75]
-
[76]
N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, J. Chem. Phys. 21, 1087 (1953)
work page 1953
-
[77]
W. L. Bragg and E. J. Williams, Proc. R. Soc. A 145, 699 (1934)
work page 1934
-
[78]
S. N. Khan, J. B. Staunton, and G. M. Stocks, Phys. Rev. B 93, 054206 (2016)
work page 2016
-
[79]
C. D. Woodgate and J. B. Staunton, Phys. Rev. B 105, 115124 (2022)
work page 2022
-
[80]
C. D. Woodgate, Modelling Atomic Arrangements in Multicomponent Alloys: A Perturbative, First- Principles-Based Approach, Springer Series in Materials Science, Vol. 346 (Springer Nature Switzerland, Cham, 2024)
work page 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.