Publication date: 14 July 2026
Source: Defect and Diffusion Forum Vol. 453
Author(s): Olga Kulitckaya, Vladislav Kulitckii
The vacancy-carbon interactions control defect kinetics and the properties of ferritic iron. Here we use spin‑polarized density‑functional theory on 3×3×3 bcc‑Fe supercells with one to four carbon atoms in octahedral interstitial sites to quantify how carbon modifies vacancy energetics. Two limiting families are considered: dilute configurations with C far from the vacancy and compact VCₙ clusters (n=1-4) with C placed in the nearest octahedral shell. For the dilute case, the vacancy formation energy remains close to that of pure Fe. In contrast, for compact clusters the effective formation energy of a vacancy bound to carbon, Ef_VC, decreases markedly with increasing n, while the total binding energy increases and then saturates. The incremental stabilization Eadd stays positive up to n = 3 and turns negative for n=4. Etrap is small and positive for n=1-2, but becomes negative for n ≥ 3, consistent with carbon‑rich microenvironments biasing vacancies into VCₙ states. Increasing n moderately reduces N(EF), particularly in the minority-spin channel, which is consistent with strengthened Fe-C hybridization and the larger binding energies of the VCn complexes. Finally, using a Seydel thermodynamic trapping model parameterized by our ab initio Ebind (VCn) we predict effective vacancy diffusivities Dv_eff(T,Ctot) that exhibit trends in line with prior analyses, while reflecting the stronger trapping implied by our energetics. Consistent with our DOS analysis, increasing carbon content strengthens Fe-C hybridisation, shifts Fe d states to lower energies and reduces N(EF), indicating a gradual transition from metallic Fe-Fe bonding to a more covalent Fe-C bond character. This work closes an important gap between electronic-structure data and mesoscale modelling by providing a consistent set of vacancy-carbon energetics and effective vacancy diffusivities for dilute C in α-Fe, which serves as a model system for ferritic Fe-based alloys.
The vacancy-carbon interactions control defect kinetics and the properties of ferritic iron. Here we use spin‑polarized density‑functional theory on 3×3×3 bcc‑Fe supercells with one to four carbon atoms in octahedral interstitial sites to quantify how carbon modifies vacancy energetics. Two limiting families are considered: dilute configurations with C far from the vacancy and compact VCₙ clusters (n=1-4) with C placed in the nearest octahedral shell. For the dilute case, the vacancy formation energy remains close to that of pure Fe. In contrast, for compact clusters the effective formation energy of a vacancy bound to carbon, Ef_VC, decreases markedly with increasing n, while the total binding energy increases and then saturates. The incremental stabilization Eadd stays positive up to n = 3 and turns negative for n=4. Etrap is small and positive for n=1-2, but becomes negative for n ≥ 3, consistent with carbon‑rich microenvironments biasing vacancies into VCₙ states. Increasing n moderately reduces N(EF), particularly in the minority-spin channel, which is consistent with strengthened Fe-C hybridization and the larger binding energies of the VCn complexes. Finally, using a Seydel thermodynamic trapping model parameterized by our ab initio Ebind (VCn) we predict effective vacancy diffusivities Dv_eff(T,Ctot) that exhibit trends in line with prior analyses, while reflecting the stronger trapping implied by our energetics. Consistent with our DOS analysis, increasing carbon content strengthens Fe-C hybridisation, shifts Fe d states to lower energies and reduces N(EF), indicating a gradual transition from metallic Fe-Fe bonding to a more covalent Fe-C bond character. This work closes an important gap between electronic-structure data and mesoscale modelling by providing a consistent set of vacancy-carbon energetics and effective vacancy diffusivities for dilute C in α-Fe, which serves as a model system for ferritic Fe-based alloys.
You might also be interested in these eBooks
Periodical:
Pages:
75-81
DOI:
Citation:
Online since:
July 2026
Authors:
Keywords:
Export:
Price:
Permissions CCC:
Permissions PLS:
Copyright:
© 2026 Trans Tech Publications Ltd. All Rights Reserved
Share:
Citation:
* - Corresponding Author
[1] C.C. Fu, E. Meslin, A. Barbu, F. Willaime and V. Oison: Solid State Phenom. Vol. 139 (2008), p.157–164.
[2] H. Ohtsuka and K. Tsuzaki: ISIJ Int. Vol. 61 (2021), p.2677–2686.
[3] E. Hristova, R. Janisch, R. Drautz and A. Hartmaier: Comput. Mater. Sci. Vol. 50 (2011), p.1088–1096.
[4] A.T. Paxton and C. Elsässer: Phys. Rev. B Vol. 87 (2013), p.224110.
[5] O. Seydel, G. Frohberg and H. Wever: Phys. Status Solidi (a) Vol. 144 (1994), p.69–79.
[6] L. Shao: Materials Vol. 17 (2024), p.3375.
[7] E. Hayward and C.-C. Fu: Phys. Rev. B Vol. 87 (2013), p.174103.
[8] P. Giannozzi et al.: J. Chem. Phys. Vol. 152 (2020), p.154105.
[9] C. Barouh, T. Schuler, C.-C. Fu and M. Nastar: Phys. Rev. B Vol. 90 (2014), p.054112.
[10] C. Domain, C.S. Becquart and J. Foct: Phys. Rev. B Vol. 69 (2004), p.144112.
[11] L. Wei, Z. Zhao, Y. Li, Q. Zheng, C. Zhang, J. Li, G. Zhao, B. Da and Z. Zeng: arXiv:2205.11050 (2022).
[12] V. Kulitckii, A. Schneider, O. Lukianova, G. Wilde, C.-C. Fu and S. Divinski: Acta Mater. Vol. 251 (2023), p.118883.