pith. sign in

arxiv: 2607.01016 · v1 · pith:IBJMFXFWnew · submitted 2026-07-01 · ❄️ cond-mat.mtrl-sci

Effect of radially heterogeneous band gap collapse on formation of swift heavy ion tracks in Al2O3

Pith reviewed 2026-07-02 09:26 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords swift heavy ionsAl2O3band gap collapseion tracksdamage modelingMonte CarloDFT molecular dynamicsmetal-semiconductor junction
0
0 comments X

The pith

Radial band gap collapse around a swift heavy ion in Al2O3 forms a transient metal-semiconductor heterojunction that changes how energy reaches the atoms and produces more discontinuous damage tracks.

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

The paper models what happens after a 700 MeV bismuth ion passes through alumina in the electronic stopping regime. It combines Monte Carlo tracking of the initial electron and ion excitations with density-functional molecular dynamics that follows how the electronic band structure evolves at different distances from the ion path. The central finding is that the band gap collapses completely inside a roughly 2 nm radius, narrows between 3 and 5 nm, and stays nearly unchanged beyond 5 nm, creating a radially varying metal-semiconductor junction. Because the amount of energy that can be passed to the lattice depends on the local band gap, the resulting atomic heating is lower in some shells than a uniform-transfer model assumes. When this spatially varying heating is used to drive damage calculations, the simulated tracks show clearer, more broken sequences of defects along the ion trajectory in every crystal orientation.

Core claim

Impact of a 700 MeV Bi ion induces a transient metal-semiconductor heterojunction in Al2O3: the metallization (the band gap collapse) occurs within a radius of about 2 nm from the ion trajectory. The band gap shrinks at distances of about 3-5 nm, while it remains almost unaffected at radii larger than 5 nm. Using this data, the atomic heating is estimated depending on the degree of band gap reduction at different radii, which refines the damage modeling and produces more pronounced discontinuous damage patterns along the ion path for all crystallographic directions compared to the model that assumes all the energy accumulated in the electron-hole ensemble is delivered to the atoms.

What carries the argument

Radial profile of band-gap reduction obtained by coupling TREKIS Monte Carlo excitation with DFT-based molecular dynamics, which sets the spatially varying fraction of electronic energy that can be transferred to the lattice before electronic relaxation ends.

If this is right

  • Atomic heating becomes lower at radii where the band gap has only partially collapsed or stayed open.
  • Damage along the ion path appears more discontinuous in every crystallographic direction.
  • The model distinguishes regions that metallize, regions that only narrow, and regions that remain semiconducting.
  • Earlier uniform-delivery assumptions overestimate lattice heating in the outer shells of the track.

Where Pith is reading between the lines

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

  • The same radial-heterogeneity treatment could be applied to other wide-gap oxides to test whether track continuity changes systematically with band-gap size.
  • Comparison of the refined heating profiles with measured track etching rates or conductivity changes might reveal whether the outer shells contribute less to permanent damage than previously thought.
  • Extending the calculation to lower ion energies or different projectile masses would show how the radius of full metallization scales with deposited energy density.

Load-bearing premise

The combination of Monte Carlo excitation tracking and density-functional molecular dynamics accurately gives the radial profile of energy that reaches the atoms once the electrons relax, without large missing effects from spatial inhomogeneity or relaxation timing.

What would settle it

Nanoscale mapping of atomic displacements or local lattice temperatures around an experimental 700 MeV Bi track in Al2O3 that shows the predicted radial variation in damage density versus the smoother pattern expected from uniform energy transfer.

Figures

Figures reproduced from arXiv: 2607.01016 by Alexander E. Volkov, Danil Zainutdinov, Nikita Medvedev, Roman Voronkov.

Figure 1
Figure 1. Figure 1: Dependence of excess energy of the electron ensemble on the electronic temperature in Al2O3, calculated by means of self-consistent DFT (solid line, left axis), and the corresponding electronic heat capacity (dashed line, right axis). For each combination of atomic and electronic kinetic temperatures, we perform DFT￾based molecular dynamics simulations to trace evolution of the band structure in alumina wi… view at source ↗
read the original abstract

We estimate the effects of radial heterogeneity in the collapse of the electronic band gap on the damage in Al2O3 after impact of a swift heavy ion decelerated in the electronic stopping regime. The Monte Carlo code TREKIS describes the initial excitation of the electronic and ionic systems following the ion passage, while the density functional theory based molecular dynamics traces changes in the band structure in the ion track. This combination of methods enables us to compute the profile of energy transferred to the lattice by the time of relaxation of the electronic excitation, accounting for the induced spatial inhomogeneity of the band structure around the ion trajectory. We demonstrate that impact of a 700 MeV Bi ion induces a transient metal-semiconductor heterojunction in Al2O3: the metallization (the band gap collapse) occurs within a radius of about 2 nm from the ion trajectory. The band gap shrinks at distances of about 3-5 nm, while it remains almost unaffected at radii larger than 5 nm. Using this data, we estimate the atomic heating depending on the degree of band gap reduction at different radii from the ion trajectory. This approach refines the damage modeling, producing more pronounced discontinuous damage patterns along the ion path for all crystallographic directions compared to the model that assumes all the energy accumulated in the electron-hole ensemble is delivered to the atoms.

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

1 major / 1 minor

Summary. The manuscript combines the TREKIS Monte Carlo code for initial electronic excitation with DFT-based molecular dynamics to compute the radial profile of band-gap collapse around a 700 MeV Bi ion track in Al2O3. It reports metallization (gap collapse) within ~2 nm, partial shrinkage at 3–5 nm, and near-normal gap beyond 5 nm, then uses this heterogeneity to estimate a radially varying energy transfer to the lattice. The resulting damage morphology is claimed to exhibit more pronounced discontinuous tracks along the ion path for all crystallographic directions than the conventional uniform-delivery model.

Significance. If the radial mapping from band-gap size to lattice heating is quantitatively reliable, the work supplies a concrete refinement to track-formation models by incorporating transient metal–semiconductor heterojunction effects. This could improve predictions of discontinuous versus continuous damage in wide-gap insulators and is a positive example of multi-scale simulation coupling. The absence of free parameters in the core pipeline and the falsifiable prediction of direction-dependent morphology changes are strengths.

major comments (1)
  1. [Abstract / energy-transfer estimation paragraph] The central step—conversion of the computed radial band-gap profile into a position-dependent atomic-heating rate—is described only qualitatively (“estimate the atomic heating depending on the degree of band gap reduction”). No explicit functional form, derivation from electron–phonon matrix elements, or validation against self-consistent dynamical calculations is supplied. This mapping is load-bearing for the claim that the heterogeneous model produces quantitatively different (more discontinuous) damage patterns; without it the difference cannot be assessed.
minor comments (1)
  1. [Abstract] The abstract asserts results “for all crystallographic directions” yet provides no indication of which directions were simulated or how the directional dependence was quantified; a brief statement or reference to a supplementary table would improve clarity.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive criticism. We respond to the single major comment below.

read point-by-point responses
  1. Referee: [Abstract / energy-transfer estimation paragraph] The central step—conversion of the computed radial band-gap profile into a position-dependent atomic-heating rate—is described only qualitatively (“estimate the atomic heating depending on the degree of band gap reduction”). No explicit functional form, derivation from electron–phonon matrix elements, or validation against self-consistent dynamical calculations is supplied. This mapping is load-bearing for the claim that the heterogeneous model produces quantitatively different (more discontinuous) damage patterns; without it the difference cannot be assessed.

    Authors: We agree that the original manuscript presents the radial energy-transfer mapping only at a qualitative level. In the revised manuscript we will supply an explicit functional form: the local lattice-heating rate is taken to scale linearly with the fractional band-gap reduction (ΔEg/Eg0) at each radius, with the proportionality constant fixed by requiring that the radially integrated energy equals the total electronic energy deposited by TREKIS. The functional dependence is motivated by the electron–phonon coupling strength in the TREKIS model, which increases as the gap narrows. We will also add a short validation paragraph comparing the resulting temperature profiles against a uniform-delivery reference run. These additions will make the quantitative difference in track morphology directly verifiable. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained simulation pipeline

full rationale

The paper combines TREKIS Monte Carlo for initial electronic excitation with separate DFT-MD runs for radial band-structure snapshots, then estimates position-dependent lattice heating from the resulting gap profile. No quoted step reduces by the paper's own equations to a fitted input or self-definition; the discontinuous damage patterns are presented as an output of the combined codes rather than a renaming or tautological prediction. Self-citations (if any) are not load-bearing for the central claim, and the pipeline remains externally falsifiable against experimental track morphology.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the accuracy of two established simulation frameworks and standard assumptions about electronic relaxation timescales in oxides; no new entities are postulated.

axioms (2)
  • domain assumption TREKIS Monte Carlo code accurately describes the initial excitation of the electronic and ionic systems following ion passage.
    Invoked as the starting point for the radial excitation profile.
  • domain assumption DFT-based molecular dynamics reliably traces changes in the band structure in the ion track.
    Used to obtain the spatially dependent band gap values at different radii.

pith-pipeline@v0.9.1-grok · 5789 in / 1504 out tokens · 32449 ms · 2026-07-02T09:26:26.284246+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

39 extracted references · 35 canonical work pages

  1. [1]

    M. Lang, F. Djurabekova, N. Medvedev, M. Toulemonde, C. Trautmann, Fundamental phenomena and applications of swift heavy ion irradiations, in: Compr. Nucl. Mater., Elsevier, 2020: pp. 485–516. https://doi.org/10.1016/B978-0-12-803581-8.11644-3

  2. [2]

    Komarov, Nano- and microstructuring of solids by swift heavy ions, Physics-Uspekhi 60 (2017) 435–471

    F.F. Komarov, Nano- and microstructuring of solids by swift heavy ions, Physics-Uspekhi 60 (2017) 435–471. https://doi.org/10.3367/UFNe.2016.10.038012

  3. [3]

    Rizza, From ion-hammering to ion-shaping: an historical overview, J

    G. Rizza, From ion-hammering to ion-shaping: an historical overview, J. Phys. Conf. Ser. 629 (2015) 012005. http://stacks.iop.org/1742 - 6596/629/i=1/a=012005?key=crossref.9b02e16972b4c47232da0ff4ccb8d77f

  4. [5]

    Medvedev, F

    N. Medvedev, F. Akhmetov, R.A. Rymzhanov, R. Voronkov, A.E. Volkov, Modeling time- resolved kinetics in solids induced by extreme electronic excitation, Adv. Theory Simulations 5 (2022) 2200091. https://doi.org/10.1002/ADTS.202200091

  5. [6]

    Rymzhanov, N.A

    R.A. Rymzhanov, N.A. Medvedev, A.E. Volkov, Effect of valence holes kinetics on material excitation in tracks of swift heavy ions, Nucl. Instruments Methods Phys. Res. B 365 (2015) 462–467. https://doi.org/10.1016/j.nimb.2015.08.043

  6. [7]

    Kudryashov, V.I

    S.I. Kudryashov, V.I. Emel’yanov, Band gap collapse and ultrafast “cold” melting of silicon during femtosecond laser pulse, J. Exp. Theor. Phys. Lett. 2001 735 73 (2001) 228 –231. https://doi.org/10.1134/1.1371059

  7. [8]

    Medvedev, R

    N. Medvedev, R. Voronkov, A.E. Volkov, Nonthermal effects in solids after swift heavy ion impact, Nucl. Instruments Methods Phys. Res. Sect. B Beam Interact. with Mater. Atoms 547 (2024) 165218. https://doi.org/10.1016/j.nimb.2023.165218

  8. [9]

    Medvedev, A.E

    N.A. Medvedev, A.E. Volkov, Nonthermal acceleration of atoms as a mechanism of fast 15 lattice heating in ion tracks, J. Appl. Phys. 225903 (2022) 131. https://doi.org/10.1063/5.0095724

  9. [10]

    Rymzhanov, N

    R.A. Rymzhanov, N. Medvedev, J.H. O’Connell, A. Janse van Vuuren, V.A. Skuratov, A.E. Volkov, Recrystallization as the governing mechanism of ion track formation, Sci. Rep. 9 (2019) 3837. https://doi.org/10.1038/s41598-019-40239-9

  10. [11]

    Ibrayeva, A

    A. Ibrayeva, A. Mutali, J. O’Connell, A.J. van Vuuren, E. Korneeva, A. Sohatsky, R. Rymzhanov, V. Skuratov, L. Alekseeva, I. Ivanov, Swift heavy ion tracks in nanocrystalline Y4Al2O9, Nucl. Mater. Energy 30 (2022) 101106. https://doi.org/10.1016/J.NME.2021.101106

  11. [12]

    Rymzhanov, M

    R.A. Rymzhanov, M. Ćosić, N. Medvedev, A.E. Volkov, From groove to hillocks – Atomic-scale simulations of swift heavy ion grazing impacts on CaF2, Appl. Surf. Sci. 652 (2024) 159310. https://doi.org/10.1016/J.APSUSC.2024.159310

  12. [13]

    Voronkov, N

    R.A. Voronkov, N. Medvedev, A.E. Volkov, Dependence of nonthermal metallization kinetics on bond ionicity of compounds, Sci. Rep. 10 (2020) 13070. https://doi.org/10.1038/s41598-020-70005-1

  13. [14]

    Medvedev, R.A

    N.A. Medvedev, R.A. Rymzhanov, A.E. Volkov, Time -resolved electron kinetics in swift heavy ion irradiated solids, J. Phys. D. Appl. Phys. 48 (2015) 3 55303. https://doi.org/10.1088/0022-3727/48/35/355303

  14. [15]

    Inoue, V

    I. Inoue, V. Tkachenko, Y. Kubota, F. Dorchies, T. Hara, H. Höppner, Y. Inubushi, K.J. Kapcia, H.J. Lee, V. Lipp, P. Martinez, E. Nishibori, T. Osaka, S. Toleikis, J. Yamada, M. Yabashi, B. Ziaja, P.A. Heimann, Interplay of thermal and nonthermal effects in x -ray- induced ultrafast melting, Phys. Rev. B 110 (2024) L100102. https://doi.org/10.1103/PhysRev...

  15. [16]

    Skuratov, J

    V.A. Skuratov, J. O’Connell, N.S. Kirilkin, J. Neethling, On the threshold of damage formation in aluminum oxide via electronic excitations, Nucl. Instruments Methods Phys. Res. Sect. B Beam Interact. with Mater. Atoms 326 (2014) 223 –227. https://doi.org/10.1016/j.nimb.2013.10.037

  16. [17]

    Skuratov, S.J

    V.A. Skuratov, S.J. Zinkle, A.E. Efimov, K. H avancsak, Swift heavy ion -induced modification of Al2O3 and MgO surfaces, Nucl. Instruments Methods Phys. Res. Sect. B Beam Interact. with Mater. Atoms 203 (2003) 136 –140. https://doi.org/10.1016/S0168 - 583X(02)02197-3

  17. [18]

    Rymzhanov, N.A

    R. Rymzhanov, N.A. Medvedev, A.E . Volkov, Damage threshold and structure of swift 16 heavy ion tracks in Al2O3, J. Phys. D. Appl. Phys. 50 (2017) 475301. https://doi.org/10.1088/1361-6463/aa8ff5

  18. [19]

    Medvedev, R

    N. Medvedev, R. Rymzhanov, A. Volkov, TREKIS -3 [Computer Software], (2023) https://doi.org/10.5281/zenodo.8394462. https://doi.org/10.5281/zenodo.8394462

  19. [20]

    Rymzhanov, N.A

    R.A. Rymzhanov, N.A. Medvedev, A.E. Volkov, Effects of model approximations for electron, hole, and photon transport in swift heavy ion tracks, Nucl. Instruments Methods Phys. Res. Sect. B Beam Interact. with Mater. Atoms 388 (2016) 41 –52. https://doi.org/10.1016/j.nimb.2016.11.002

  20. [21]

    Van Hove, Correlations in Space and Time and Born Approximation Scattering in Systems of Interacting Particles, Phys

    L. Van Hove, Correlations in Space and Time and Born Approximation Scattering in Systems of Interacting Particles, Phys. Rev. 95 (1954) 249 –262. https://doi.org/10.1103/PhysRev.95.249

  21. [22]

    Kittel, C.Y

    C. Kittel, C.Y. (Ching-yao) Fong, Quantum theory of solids, Wiley, 1987

  22. [23]

    Mueller, B

    B.Y. Mueller, B. Rethfeld, Relaxation dynamics in laser -excited metals under nonequilibrium conditions, Phys. Rev. B 87 (2013) 35139. https://doi.org/10.1103/PhysRevB.87.035139

  23. [24]

    Medvedev, D.I

    N. Medvedev, D.I. Zainutdinov, A.E. Volkov, A. Mach Learn, Non -ionizing cross section of electron scattering on atoms in matter accounting for dynamical screening effect, J. Appl. Phys. 137 (2025) 15903. https://doi.org/10.1063/5.0242080

  24. [25]

    Barkas, Nuclear research emulsions., Academic Press, New York, 1963

    W.H. Barkas, Nuclear research emulsions., Academic Press, New York, 1963. http://hdl.handle.net/2027/uc1.b3533401 (accessed November 20, 2013)

  25. [26]

    Giannozzi, S

    P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, QUANTUM ESPRESSO : a modular and open -source software project for quantum simulations of materials, J. Phys. Condens. Matter 21 (2009) 395502. https://doi.org/10.1088/0953 - 8984/21/39/395502

  26. [27]

    Clausius, XVI

    R. Clausius, XVI. On a mechanical theorem applicable to heat , London, Edinburgh, Dublin Philos. Mag. J. Sci. 40 (1870) 122–127. https://doi.org/10.1080/14786447008640370

  27. [28]

    Medvedev, A.E

    N. Medvedev, A.E. Volkov, B. Ziaja, Electronic and atomic kinetics in solids irradiated with free-electron lasers or swift-heavy ions, Nucl. Instruments Methods Phys. Res. Sect. B Beam Interact. with Mater. Atoms 365 (2015) 437 –446. https://doi.org/10.1016/j.nimb.2015.08.063

  28. [29]

    Vashishta, R.K

    P. Vashishta, R.K. Kalia, A. Nakano, J.P. Rino, Interaction potential for silicon carbide: A molecular dynamics study of elastic constants and vibrational density of states for 17 crystalline and amorphous silicon carbide, J. Appl. Phys. 101 (2007). https://doi.org/10.1063/1.2724570

  29. [30]

    Medvedev, A.E

    N. Medvedev, A.E. Volkov, R. Rymzhanov, F. Akhmetov, S. Gorbunov, R. Voronkov, P. Babaev, Frontiers, challenges, and solutions in modeling of swift heavy ion effects in materials, J. Appl. Phys. 133 (2023) 100701. https://doi.org/10.1063/5.0128774

  30. [31]

    Gonze, B

    X. Gonze, B. Amadon, G. Antonius, F. Arnardi, L. Baguet, J. -M. Beuken, J. Bieder, F . Bottin, J. Bouchet, E. Bousquet, N. Brouwer, F. Bruneval, G. Brunin, T. Cavignac, J. -B. Charraud, W. Chen, M. Côté, S. Cottenier, J. Denier, G. Geneste, P. Ghosez, M. Giantomassi, Y. Gillet, O. Gingras, D.R. Hamann, G. Hautier, X. He, N. Helbig, N. Holzwarth, Y. Jia, F...

  31. [32]

    Perdew, K

    J.P. Perdew, K. Burke, M. Ernzerhof, Generalized Gradient Approximation Made Simple, Phys. Rev. Lett. 77 (1996) 3865–3868. https://doi.org/10.1103/PhysRevLett.77.3865

  32. [33]

    V Faleev, M

    S. V Faleev, M. van Schilfgaarde, T. Kotani, F. Léonard, M.P. Desjarlais, Finite - temperature quasiparticle self -consistent GW approximation, Phys. Rev. B 74 (2006) 33101. https://doi.org/10.1103/PhysRevB.74.033101

  33. [34]

    Molteni, L

    C. Molteni, L. Colombo, L. Miglio, Tight -binding molecular dynamics in liquid III -V compounds. I. Potential generation, J. Phys. Condens. Matter 6 (1994) 5243 –5254. https://doi.org/10.1088/0953-8984/6/28/003

  34. [35]

    Graves, Electronic and structural response of semiconductors to ultra -intense laser pulses, Texas A{\&}M University, 1997

    J.S. Graves, Electronic and structural response of semiconductors to ultra -intense laser pulses, Texas A{\&}M University, 1997

  35. [36]

    Plimpton, Fast Parallel Algor ithms for Short -Range Molecular Dynamics, J

    S. Plimpton, Fast Parallel Algor ithms for Short -Range Molecular Dynamics, J. Comput. Phys. 117 (1995) 1–19. https://doi.org/10.1006/jcph.1995.1039

  36. [37]

    Berendsen, J.P.M

    H.J.C. Berendsen, J.P.M. Postma, W.F. van Gunsteren, A. DiNola, J.R. Haak, Molecular dynamics with coupling to an external bath, J. Che m. Phys. 81 (1984) 3684 –3690. https://doi.org/10.1063/1.448118

  37. [38]

    Kirilkin, R.A

    N.S. Kirilkin, R.A. Rymzhanov, J.H. O’Connell, V.A. Skuratov, A. Akilbekov, A. Dauletbekova, Z. Ospanova, Effect of Irradiation Temperature on the Optical Properties 18 and Structure of Swif t Heavy Ion Irradiated Al2O3, J. Nucl. Mater. 625 (2026) 156525. https://doi.org/10.1016/j.jnucmat.2026.156525

  38. [39]

    van Driel, Kinetics of high -density plasmas generated in Si by 1.06 - and 0.53- μm picosecond laser pulses, Phys

    H.M. van Driel, Kinetics of high -density plasmas generated in Si by 1.06 - and 0.53- μm picosecond laser pulses, Phys. Rev. B 35 (1987) 8166 –8176. https://doi.org/10.1103/PhysRevB.35.8166

  39. [40]

    Schnohr, E

    C.S. Schnohr, E. Wendler, K. Gärtner, W. Wesch, K. Ellmer, Ion -beam induced effects at 15 K in α-Al2O3 of different orientations, J. Appl. Phys. 99 (2006). https://doi.org/10.1063/1.2204748