REVIEW 5 major objections 4 minor 33 references
Anisotropic Core-Shell Swift Heavy Ion Tracks in beta-Ga2O3
T0 review · 5 major / 4 minor · reviewed 2026-08-02 · deepseek-v4-flash
Pith's one-line read This paper shows that swift heavy ion tracks in beta-Ga2O3 are intrinsically anisotropic core-shell structures—an amorphous core wrapped in a gamma-phase shell—and that the anisotropy is set by orientation-dependent recovery controlled by e
desk verdict A plausible but under-supported simulation study of anisotropic core-shell SHI tracks in beta-Ga2O3; the new physics is interesting but the quantitative and methodological details need a major revision. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The argument is carried by a multiscale chain: Monte Carlo particle transport produces the spatial ionization profile; a two-temperature model converts it to a lattice energy deposition profile; and molecular dynamics with a machine-learned interatomic potential evolves the atomic structure. The central conceptual object is the core-shell track (amorphous core plus gamma-phase shell), characterized by local configurational entropy. The orientation-dependent recovery is tied to direction-dependent elastic stiffness, with the highest stiffness along [010] enabling faster recrystallization.
What would settle it
A direct test is to measure track cross-sections for irradiation perpendicular to (010), (001), and (201) planes at about 44 keV/nm and compare with the predicted shapes—especially the strongly asymmetric profile for (010). If the amorphous core is not consistently narrowest along [010], the elastic-stiffness mechanism fails. A complementary check is to rerun the two-temperature model with different published electron-phonon coupling constants and see whether the anisotropy and quantitative sizes survive.
Extended reading notes
Core claim
Swift heavy ion irradiation of monoclinic beta-Ga2O3 produces core-shell tracks whose morphology is governed by recovery dynamics, not by primary damage. The paper identifies a universal structural sequence with increasing electronic energy loss: complete lattice recovery; recrystallization into metastable gamma-Ga2O3; and finally an amorphous core with a gamma-phase shell. Despite isotropic energy deposition, the amorphous core is consistently more confined along [010], the direction of highest Young's modulus (~284 GPa), because stiff directions recover faster. The simulations quantitatively reproduce experimental track diameters for (100) irradiation from 18 to 44 keV/nm, challenging the
Load-bearing premise
The result leans on the specific two-temperature-model parameters (electron-phonon coupling and energy source) and on the accuracy of the machine-learned potential for amorphous and recrystallizing gallium oxide; if either misrepresents energy transfer from electrons to the lattice, the quantitative track sizes and anisotropy could shift.
Editorial extensions
If this is right
- Track diameters and shapes in beta-Ga2O3 can be predicted from elastic anisotropy rather than assumed cylindrical symmetry.
- The metastable gamma-phase shell is a stable product of electronic excitation alone, not only of nuclear collision cascades.
- Irradiation perpendicular to (100) gives smaller residual tracks at high energy loss, making orientation a lever for radiation-hardness engineering.
- Simulated track sizes match experiments, supporting use of this multiscale approach for other low-symmetry oxides.
Reading between the lines
- If elastic stiffness controls recovery, then strain engineering or alloying that changes directional moduli could tune track morphology and radiation tolerance.
- The same core-shell picture may apply to other anisotropic oxides, where experimental tracks might show similar orientation-dependent fine structure.
- A testable extension: measure track cross-sections for (010), (001), and (201) irradiations to check the predicted shapes, since the paper only compares (100) sizes with experiments.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a multiscale simulation study of swift heavy ion (SHI) track formation in monoclinic beta-Ga2O3, combining Geant4 Monte Carlo electronic energy deposition, a two-temperature model (TTM), and machine-learned interatomic potential molecular dynamics (tabGAP). The authors report a sequence of structural responses with increasing electronic energy loss (Se): full recovery, recrystallization into a metastable gamma-Ga2O3 phase, and core-shell track formation (amorphous core + gamma shell) at high Se. They claim that despite isotropic initial electronic energy deposition, the final track morphology is strongly anisotropic across orientations, governed by orientation-dependent recovery dynamics, with the highest stiffness along [010] promoting recrystallization. They also claim quantitative agreement with experimental track sizes over Se = 18–44 keV/nm.
Significance. If substantiated, the paper would provide an atomistic mechanism for anisotropic core-shell track formation in a low-symmetry oxide, establishing elastic anisotropy as a key factor controlling track recovery. The multiscale framework (Geant4 + TTM + ML-MD with a generally applicable tabGAP potential) is state-of-the-art, and the predicted orientation-dependent track shapes are concrete, experimentally testable. The paper explicitly separates primary damage from recovery, which is a valuable conceptual contribution. However, the quantitative and mechanistic claims currently rest on several incompletely specified procedures, as detailed below.
major comments (5)
- [Results, 'identical lattice energy profile'; Methods, 'Molecular dynamics simulation'] The central claim—that anisotropy arises from recovery, not primary damage—requires the initial volumetric energy deposition to be isotropic. The text states 'an identical lattice energy profile was imposed for all orientations' but does not specify whether this is a per-atom kinetic energy profile or a per-volume energy density. The observation that the (100) plane has the highest fraction of high-Ek atoms, attributed to its higher planar density, indicates a per-atom profile, which makes the total deposited energy per unit path length orientation-dependent. The invariance of the combined (core+gamma) radius is a weak test of isotropy. Please specify the injection method explicitly and demonstrate that the volumetric energy density is orientation-independent.
- [Results, comparison with experiments; Table II] The claim of 'excellent quantitative agreement' is not supported by the numbers. For (100) at 44 keV/nm, the simulated track is elliptical with diameters 8.55 nm and 11.61 nm. The cited experimental values are 7.8±0.9 nm and 8.3±0.4 nm. If the experiments measure the track cross-section in the (100) plane, the simulated major axis exceeds the Tracy et al. value by roughly 3.3 nm, far outside the stated uncertainty. The relevant simulated quantity (minor axis, area-equivalent diameter, or full ellipse) must be compared explicitly, with the experimental measurement geometry stated.
- [Methods, 'Molecular dynamics simulation'; Table II] Three independent runs are performed for each configuration, but Table II reports only single values for track size, gamma-phase size, and recovery ratio, with no standard deviation or range. Since quantitative agreement with experiment is a central claim, the run-to-run spread must be reported to assess the significance of the comparisons. Without error bars, the apparent agreement could be fortuitous.
- [Results, 'entropy-based criterion'; Methods, 'Molecular dynamics simulation'] The size criterion used to extract track and gamma-phase radii is not specified. The text says 'Based on the entropy-based criterion, the size ... are summarized in Table II,' but the exact entropy threshold(s) for the amorphous core (Omega in [-2.5,-2.0]) and the boundary of the gamma-phase shell are not defined, nor is the procedure for converting one-dimensional entropy profiles into radial extents. Without this, Table II is not reproducible. Please provide the explicit criterion and analysis procedure.
- [Methods, 'Two-temperature model'] The TTM parameters (electron-phonon coupling constant g, heat capacities Ce and Cl, thermal conductivities Ke and Kl) and the normalization of the source term Se(r,t) are not given in the main text; the text only refers to supporting information. The TTM output is the sole energy input to the MD simulations, so the quantitative track sizes depend directly on these values. Please provide the numerical values and their provenance, or ensure the SI is available with the manuscript.
minor comments (4)
- [Results, first paragraph] Typo: 'Morte Carlo' should be 'Monte Carlo'.
- [Methods, 'Molecular dynamics simulation'] The supercell atom count reads '115,2000'—likely a typo for '1,152,000'. Please correct.
- [Results, subsection on low-Se response] The reference 'Table SX' is a placeholder; the actual supplementary table number should be inserted.
- [General] Several formatting issues: missing spaces before 'beta-Ga2O3' in a few places (e.g., abstract), and 'byLAMMPSsoftware' in Methods should be 'by LAMMPS software'. The supplementary information is referenced but not included in the arXiv version; please ensure it is accessible for review.
Circularity Check
No significant circularity: simulated track-size predictions are independent MD outputs benchmarked against experiment, not fitted inputs.
full rationale
The central quantities—orientation-dependent track diameters, gamma-phase shell sizes, and recovery ratios—are outputs of the multiscale simulation chain (Geant4 → TTM → ML-MD with tabGAP), not parameters fitted to the experimental diameters quoted for comparison. The Se values are independently calculated inputs, the TTM equations are standard, and no coefficient is tuned to reproduce the measured track sizes. The paper's use of the authors' own tabGAP potential and earlier gamma-phase characterizations is methodological and interpretive; those citations do not contain the target result, and the final track-size agreement with independent experiments (Ai et al., Tracy et al.) provides external validation. The paper does not invoke a uniqueness theorem or rename a fitted parameter as a prediction. A transparency concern remains: TTM parameter values and the exact discrete projection of the continuum lattice-energy profile onto atoms are not given in the main text, so the physical isotropy of the injected energy is asserted rather than demonstrated. This is a reproducibility/validity risk, not circularity, because no output is reduced to its input by construction.
Assumptions & free parameters
free parameters (5)
- electron-phonon coupling g =
not stated
- TTM heat capacities and conductivities Ce, Cl, Ke, Kl =
not stated in main text
- Entropy-based track-size criterion (Omega plateau threshold) =
Omega in [-2.5,-2.0] defines the amorphous core; the gamma shell is defined by the sharp increase at 5.0-7.5 nm
- Kinetic-energy fraction threshold (0.5 eV) =
0.5 eV
- Young's modulus along [010] =
284 GPa
assumptions (5)
- domain assumption The two-temperature model with parameters Ce, Cl, Ke, Kl, g adequately describes the electron-lattice energy transfer in beta-Ga2O3 under SHI.
- domain assumption The tabGAP machine-learned interatomic potential accurately describes beta, gamma, amorphous Ga2O3, melting, recrystallization, and defect energies.
- domain assumption The isotropic, orientation-averaged electronic energy deposition profile can be imposed identically for the four crystal orientations by rotating the crystal, without modifying the TTM source term.
- domain assumption The 320 ps simulation time is sufficient to reach the final stable track morphology.
- domain assumption A supercell of roughly 150x300x300 A^3 and periodic boundary conditions in all directions do not significantly affect the track morphology.
Cite this review
Pith. "Pith review of Anisotropic Core-Shell Swift Heavy Ion Tracks in beta-Ga2O3." pith.science (2026). https://pith.science/paper/B7ST7L5Q
@misc{pith2026260213614,
author = {Pith},
title = {Pith review of: Anisotropic Core-Shell Swift Heavy Ion Tracks in beta-Ga2O3},
year = {2026},
howpublished = {\url{https://pith.science/paper/B7ST7L5Q}},
note = {Machine review of arXiv:2602.13614}
}
abstract
Swift heavy ion (SHI) irradiation generates nanoscale ion tracks through intense electronic excitation, yet the microscopic mechanisms governing their morphology and phase stability in low symmetry oxides remain poorly understood. Here, a multiscale atomistic simulation framework is used to investigate the formation and recovery of SHI-induced tracks in monoclinic $\beta$-Ga2O3 over a wide range of electronic energy losses (Se) and crystallographic orientations. A sequence of distinct structural responses is identified with increasing Se: (i) complete lattice recovery at low Se; (ii) recrystallization into a metastable $\gamma$-Ga2O3 phase at intermediate Se; and (iii) the formation of core-shell ion tracks at high Se, consisting of an amorphous core surrounded by a recrystallized $\gamma$-phase shell. Despite the essentially isotropic initial energy deposition, the final ion-track morphology exhibits pronounced crystallographic anisotropy, governed by orientation-dependent recovery dynamics. The superior recrystallization along the [010] direction is attributed to its exceptionally high elastic stiffness. Notably, SHI irradiation perpendicular to the (100) plane induces a more severe structural response at low Se ($\le$ 10 keV/nm), however, at higher Se, it yields a smaller residual ion track compared to the other orientations. The simulated ion-track sizes show excellent quantitative agreement with the available experimental measurements over a wide range of Se values. These findings establish a unified atomic-scale picture of core-shell track formation and anisotropic recovery in $\beta$-Ga2O3.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[10]
A. Abdullaev, J. Garcia Fernandez, C. Nozais, J. O’Connell, R. Tlegenov, K. Sekerbayev, A. Azarov, A. Leino, T. Fernandez Bouvier, J. Zhao, A. Ar- timez Pena, N. Medvedev, Z. Utegulov, Ø. Prytz, F. Djurabekova, and A. Kuznetsov, Ions leaving no tracks, arXiv 10.48550/arXiv.2509.07440 (2025), 2509.07440
work page Pith review arXiv doi:10.48550/arxiv.2509.07440 2025
-
[1]
M. C. Ridgway, R. Giulian, D. J. Sprouster, P. Kluth, L. L. Araujo, D. J. Llewellyn, A. P. Byrne, F. Kremer, P. F. P. Fichtner, G. Rizza, H. Amekura, and M. Toule- monde, Role of thermodynamics in the shape trans- formation of embedded metal nanoparticles induced by swift heavy-ion irradiation, Phys. Rev. Lett.106, 095505 (2011)
2011
-
[2]
Ochedowski, O
O. Ochedowski, O. Osmani, M. Schade, B. K. Bussmann, B. Ban-d’Etat, H. Lebius, and M. Schleberger, Graphitic nanostripes in silicon carbide surfaces created by swift heavy ion irradiation, Nat. Commun.5, 3913 (2014)
2014
-
[3]
Zhang, R
Y. Zhang, R. Sachan, O. H. Pakarinen, M. F. Chisholm, P. Liu, H. Xue, and W. J. Weber, Ionization-induced annealing of pre-existing defects in silicon carbide, Nat. Commun.6, 8049 (2015)
2015
-
[4]
M. Lang, F. Djurabekova, N. Medvedev, M. Toulemonde, and C. Trautmann, Fundamental phenomena and appli- cations of swift heavy ion irradiations, inComprehensive Nuclear Materials, edited by R. J. M. Konings and R. E. Stoller (Elsevier, Oxford, 2020) 2nd ed., pp. 485–516
2020
-
[5]
Medvedev, A
N. Medvedev, A. E. Volkov, R. Rymzhanov, F. Akhme- tov, S. Gorbunov, R. Voronkov, and P. Babaev, Frontiers, challenges, and solutions in modeling of swift heavy ion effects in materials, J. Appl. Phys.133, 100701 (2023)
2023
-
[6]
Amekura, A
H. Amekura, A. Chettah, K. Narumi, A. Chiba, Y. Hirano, K. Yamada, S. Yamamoto, A. A. Leino, F. Djurabekova, K. Nordlund, N. Ishikawa, N. Okubo, and Y. Saitoh, Latent ion tracks were finally observed in diamond, Nat. Commun.15, 1786 (2024)
2024
-
[7]
M. C. Ridgway, T. Bierschenk, R. Giulian, B. Afra, M. D. Rodriguez, L. L. Araujo, A. P. Byrne, N. Kirby, O. H. Pakarinen, F. Djurabekova, K. Nordlund, M. Schle- berger, O. Osmani, N. Medvedev, B. Rethfeld, and P. Kluth, Tracks and voids in amorphous Ge induced by swift heavy-ion irradiation, Phys. Rev. Lett.110, 245502 (2013)
2013
Show all 33 references
-
[8]
M. C. Sequeira, J.-G. Mattei, H. Vazquez, F. Djurabekova, K. Nordlund, I. Monnet, P. Mota- Santiago, P. Kluth, C. Grygiel, S. Zhang, E. Alves, and K. Lorenz, Unravelling the secrets of the resistance of GaN to strongly ionising radiation, Commun. Phys.4, 51 (2021)
2021
-
[9]
P. P. Hu, Z. N. Jin, L. J. Xu, S. X. Zhang, P. F. Zhai, J. Zeng, Z. Z. Li, X. Y. Yan, L. Cai, and J. Liu, Latent tracks induced by swift heavy ions in gallium nitride, Nucl. Instrum. Methods Phys. Res., Sect. B568, 165855 (2025)
2025
-
[11]
A. K. Gupta, E. Zarkadoula, B. L. Musico, J. A. Hachtel, M. A. Roldan, V. R. Paduri, C. Harris, R. Subedi, M. Ziatdinov, S. V. Kalinin, C. Trautmann, V. Kep- pens, J. Liu, Y. Zhang, W. J. Weber, and R. Sachan, Unveiling swift heavy ion track morphology in Sr-based high-entropy...
2025
-
[12]
Wierbik, H
J. Wierbik, H. Heimes, C. Notthoff, S. Dutt, T. Al- wadi, A. Kiy, P. Mota-Santiago, N. Kirby, and P. Kluth, Anisotropic fine structure of ion tracks in single crystals, Phys. Rev. B113, 035306 (2026)
2026
-
[13]
Kluth, C
P. Kluth, C. S. Schnohr, O. H. Pakarinen, F. Djurabekova, D. J. Sprouster, R. Giulian, M. C. Ridgway, A. P. Byrne, C. Trautmann, D. J. Cookson, K. Nordlund, and M. Toulemonde, Fine structure in swift heavy ion tracks in amorphous SiO 2, Phys. Rev. Lett.101, 175503 (2008)
2008
-
[14]
Bierschenk, A
T. Bierschenk, A. A. Leino, W. Wesch, B. Afra, M. D. Rodriguez, F. Djurabekova, L. Keller, O. H. Pakarinen, K. Nordlund, M. C. Ridgway, and P. Kluth, Forma- tion and self-organisation of nano-porosity in swift heavy ion irradiated amorphous Ge, Acta Mater.261, 119396 (2023)
2023
-
[15]
S. J. Pearton, J. Yang, P. H. Cary, F. Ren, J. H. Kim, M. J. Tadjer, and M. A. Mastro, A review of Ga 2O3 materials, processing, and devices, Appl. Phys. Rev.5, 011301 (2018)
2018
-
[16]
W.-S. Ai, J. Liu, Q. Feng, P.-F. Zhai, P.-P. Hu, J. Zeng, S.-X. Zhang, Z.-Z. Li, L. Liu, X.-Y. Yan, and Y.-M. Sun, Degradation ofβ-Ga 2O3 Schottky barrier diode under swift heavy ion irradiation, Chin. Phys. B30, 056110 (2021)
2021
-
[17]
S. J. Pearton, F. Ren, A. Y. Polyakov, E. B. Yakimov, L. Chernyak, and A. Haque, Perspective on comparative radiation hardness of Ga 2O3 polymorphs, J. Vac. Sci. Technol. A43, 038501 (2025)
2025
-
[18]
Manikanthababu, B
N. Manikanthababu, B. R. Tak, K. Prajna, S. Sarkar, K. Asokan, D. Kanjilal, S. R. Barman, R. Singh, and B. K. Panigrahi, Swift heavy ion irradiation-induced modifications in the electrical and surface properties of β-Ga2O3, Appl. Phys. Lett.117, 142105 (2020). 10
2020
-
[19]
H. Zhu, Y. Tang, A. Zhong, L. Wang, F. Liu, P. Zhao, J. Liu, L. Shu, Z. Wu, and B. Li, Swift heavy ion irradiation-driven energy band engineering and its pro- found influence on the photoresponse ofβ-Ga 2O3 ul- traviolet photodetectors, Appl. Phys. Lett.124, 153505 (2024)
2024
-
[20]
T. Ma, X. Li, S. Zhong, Y. Feng, Y. Wang, X. Zhang, C. Peng, H. Zhang, Z. Zhang, X. Liang, and Z. Lei, Single-event burnout inβ-Ga 2O3 Schottky barrier diode induced by heavy ion irradiation, Appl. Phys. Lett.127, 123501 (2025)
2025
-
[21]
C. L. Tracy, M. Lang, D. Severin, M. Bender, C. Traut- mann, and R. C. Ewing, Anisotropic expansion and amorphization of Ga 2O3 irradiated with 946 MeV Au ions, Nucl. Instrum. Methods Phys. Res., Sect. B374, 40 (2016)
2016
-
[22]
W. Ai, L. Xu, S. Nan, P. Zhai, W. Li, Z. Li, P. Hu, J. Zeng, S. Zhang, L. Liu, Y. Sun, and J. Liu, Radiation damage inβ-Ga 2O3 induced by swift heavy ions, Jpn. J. Appl. Phys.58, 120914 (2019)
2019
-
[23]
Iancu, E
D. Iancu, E. Zarkadoula, V. Leca, A. Hotnog, Y. Zhang, W. J. Weber, and G. Veli¸ sa, Intrinsic property of defec- tiveβ-Ga 2O3 to self-heal under ionizing irradiation, Scr. Mater.268, 116858 (2025)
2025
-
[24]
J. Zhao, J. Byggm¨ astar, H. He, K. Nordlund, F. Djurabekova, and M. Hua, Complex Ga 2O3 poly- morphs explored by accurate and general-purpose machine-learning interatomic potentials, npj Comput. Mater.9, 159 (2023)
2023
-
[25]
Azarov, J
A. Azarov, J. Fern´ andez, J. Zhao,et al., Universal ra- diation tolerant semiconductor, Nature Communications 14, 4855 (2023)
2023
-
[26]
J. Zhao, J. G. Fern´ andez, A. Azarov, R. He, Ø. Prytz, K. Nordlund, M. Hua, F. Djurabekova, and A. Kuznetsov, Crystallization instead of amorphization in collision cascades in gallium oxide, Phys. Rev. Lett. 134, 126101 (2025)
2025
-
[27]
Azarov, J
A. Azarov, J. G. Fern´ andez, J. Zhao, R. He, J.- H. Park, D.-W. Jeon, Ø. Prytz, F. Djurabekova, and A. Kuznetsov, Phase glides and self-organization of atom- ically abrupt interfaces out of stochastic disorder inα- Ga2O3, Nat. Commun.16, 3245 (2025)
2025
-
[28]
Agostinelli, J
S. Agostinelli, J. Allison, K. Amako, J. Apostolakis, H. Araujo, P. Arce, M. Asai, D. Axen, S. Banerjee, G. Barrand, F. Behner, L. Bellagamba, J. Boudreau, L. Broglia, A. Brunengo, H. Burkhardt, S. Chauvie, J. Chuma, R. Chytracek, G. Cooperman, G. Cosmo, P. Degtyarenko, A. Del...
2003
-
[29]
A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolin- tineanu, W. M. Brown, P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, LAMMPS - a flexible simulation tool for particle-based mat...
2022
-
[30]
H. He, J. Zhao, J. Byggm¨ astar, R. He, K. Nord- lund, C. He, and F. Djurabekova, Threshold displace- ment energy map of Frenkel pair generation inβ-Ga 2O3 from machine-learning-driven molecular dynamics simu- lations, Acta Mater.276, 120087 (2024)
2024
-
[31]
Zhang, J
J. Zhang, J. Zhao, J. Byggm¨ astar, E. J. Frankberg, and A. Kuronen, Large-scale atomistic study of plasticity in amorphous gallium oxide with ab-initio accuracy, Sci. Rep.15, 9492 (2025)
2025
-
[32]
Stukowski, Visualization and analysis of atomistic simulation data with OVITO–the open visualization tool, Modell
A. Stukowski, Visualization and analysis of atomistic simulation data with OVITO–the open visualization tool, Modell. Simul. Mater. Sci. Eng.18, 015012 (2010)
2010
-
[33]
T. Liu, T. Wang, Z. Li, E. Zhou, J. Zhao, J. Feng, X. Fei, Y. Guo, F. Djurabekova, S. Liu, and Z. Zhang, Anisotropic kinetics of ion-irradiation- induced phase transition in gallium oxide, arXiv 10.48550/arXiv.2601.02770 (2026), 2601.02770
2026 doi
Reviewed August 2, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.