{"id":"50e08966-151e-494f-b7a7-88030bf296b6","arxiv_id":"2602.13614","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Simulations show swift heavy ion tracks in beta-Ga2O3 form anisotropic core-shell structures, where recrystallization along [010] is strongest due to high elastic stiffness.","lead":"This paper uses computer simulations to show that when heavy ions hit beta-Ga2O3 crystals, the resulting damage tracks have different shapes and sizes depending on the crystal direction, even though the initial energy deposit is the same. This matters because it explains a key mechanism of radiation damage in an important semiconductor material.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Verify TTM-to-MD energy injection is truly isotropic; otherwise the observed track anisotropy may be an artifact of the coupling method.","rationale":"The reader's weakest_assumption focused on TTM parameters and ML-potential accuracy, with a secondary mention of the asserted 'identical lattice energy profile.' I believe the more fundamental and load-bearing issue is the TTM-to-MD coupling, because it directly determines whether the initial damage is actually isotropic. If the injection is anisotropic, the separation between primary damage and recovery—the core claim—collapses, and the stiffness explanation becomes moot. This is a concrete, falsifiable concern, testable by a controlled change in the injection protocol. The paper's statement that the combined core+shell radius is orientation-independent is suggestive but not conclusive; it does not establish that the per-atom energy landscape (and hence the local melting/amorphization probability) is isotropic. The reported high-Ek fraction for the (100) plane is a red flag that the injection may be per-atom rather than per-volume, which would make the primary damage anisotropic. I therefore recommend the verdict remain CONDITIONAL, with the added condition that the authors demonstrate isotropic energy injection, either by describing their coupling scheme in detail and showing the initial damage is circular, or by performing the test I propose. If the test fails, the central claim should be rejected.","tokens_in":13214,"tokens_out":8780,"duration_ms":83348,"concrete_test":"Rerun the MD simulations for the four orientations using a strictly isotropic energy-deposition protocol: from the TTM temperature profile, compute the kinetic energy per unit volume and assign each atom a velocity such that the kinetic energy density is isotropic (i.e., per-atom energy = E_vol(r)/n_local, where n_local is the local atomic density). Compare the initial kinetic-energy distribution, the transient amorphous region, and the final track shapes against the current results. If the anisotropy disappears, the paper's central claim is an artifact of the injection method; if it persists, the claim is supported. Also report the fraction of high-Ek atoms for each orientation under this protocol to check consistency with Table SX.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central conclusion—that the observed crystallographic anisotropy in ion-track shape is caused by orientation-dependent recovery (and ultimately elastic stiffness), not by primary damage—requires that the initial lattice energy deposition be isotropic. The paper asserts 'an identical lattice energy profile was imposed for all orientations' but never specifies how the continuum TTM energy profile is projected onto discrete atoms. Two natural choices exist: (i) assign each atom the same kinetic energy at a given radius (per-atom profile), which makes volumetric energy density proportional to atomic density and hence anisotropic; or (ii) assign the same volumetric energy density, which makes per-atom energy inversely proportional to density. The reported observation that the higher-density (100) plane has the highest fraction of atoms above 0.5 eV is difficult to reconcile with (ii) and suggests (i), which would violate the isotropy premise. The invariance of the combined (core+γ) radius across orientations is taken as evidence of isotropic primary damage, but it is a weak test: different volumetric energy densities could still yield similar combined radii if the recrystallization threshold shifts. Without a demonstrated isotropic injection, the final track anisotropy could be a direct artifact of orienting the same per-atom energy profile in different crystal directions.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":13376,"tokens_out":8320,"duration_ms":70488,"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":[{"comment":"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.","section":"Results, 'identical lattice energy profile'; Methods, 'Molecular dynamics simulation'"},{"comment":"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.","section":"Results, comparison with experiments; Table II"},{"comment":"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.","section":"Methods, 'Molecular dynamics simulation'; Table II"},{"comment":"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.","section":"Results, 'entropy-based criterion'; Methods, 'Molecular dynamics simulation'"},{"comment":"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.","section":"Methods, 'Two-temperature model'"}],"minor_comments":[{"comment":"Typo: 'Morte Carlo' should be 'Monte Carlo'.","section":"Results, first paragraph"},{"comment":"The supercell atom count reads '115,2000'—likely a typo for '1,152,000'. Please correct.","section":"Methods, 'Molecular dynamics simulation'"},{"comment":"The reference 'Table SX' is a placeholder; the actual supplementary table number should be inserted.","section":"Results, subsection on low-Se response"},{"comment":"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.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The review is based solely on the main text, as the supplementary information is not part of the arXiv posting. The authors should be asked to provide the SI containing TTM parameters, the energy-injection method, and the entropy-threshold definition. The paper's qualitative phenomenology is plausible and interesting, but the isotropy premise and quantitative agreement need to be demonstrated more rigorously before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi [Colleague],\n\nQuick take: this paper is worth reading, but don't trust the 'excellent agreement' headline. The core result—orientation-dependent core-shell (amorphous + gamma) tracks in beta-Ga2O3 under SHI, with recovery anisotropy—is new and plausible. The multiscale method is standard, and the sequence from full recovery -> gamma-phase -> core-shell with increasing Se is coherent.\n\nWhat it does well: it is the first atomistic simulation, as far as I know, to resolve the core-shell structure and metastable gamma phase from electronic excitation alone in this material. The qualitative match with Tracy and Ai on track sizes is real: the simulated diameters at 18 and 44 keV/nm bracket the measured ones, and the combined core+gamma radius is impressively isotropic across four orientations, which supports the claim that the primary damage is roughly isotropic and the final anisotropy comes from recovery.\n\nNow the soft spots. The quantitative case is thinner than claimed. Three runs per condition and no reported dispersion means you cannot assess the noise. The entropy-based definition of the track radius uses an arbitrary plateau threshold; there is no sensitivity check. The TTM parameters are in a missing SI, which is a reproducibility problem. And the 'elastic stiffness controls recovery' link is asserted rather than tested—you cannot tell if the [010] recovery advantage is actually causal or correlated.\n\nThere is also a real red flag in the low-Se discussion. The paper says the (100) orientation has the highest fraction of high-kinetic-energy atoms and attributes it to the high atomic planar density of (100). That logic is backwards: if the energy is deposited per unit volume, a denser plane means the same energy spread over more atoms, lowering per-atom energy. If instead each atom gets the same energy at a given radius, the planar density is irrelevant. The paper never specifies how the TTM energy profile is projected onto atoms. This is exactly the point the stress-test raises. The combined-radius invariance is good evidence that the injection is not grossly anisotropic, but the authors need to clearly state and justify the mapping. As written, the planar density explanation is confused and makes me question the 'identical energy profile' claim.\n\nBottom line: I agree with the reader's conditional verdict. The core idea survives my reading, but it needs major revision to be convincing: report parameters, show dispersion and sensitivity, and clean up the energy-deposition description. I would send it to peer review—it is a novel result in a tech-relevant material—but I would expect substantial changes before acceptance.\n\nRecommendation: engage with it; it is worth referee time.","headline":"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.","tokens_in":14033,"tokens_out":8994,"would_cite":true,"duration_ms":85972,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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","keywords":["swift heavy ion","ion track","beta-Ga2O3","core-shell structure","anisotropic recovery","two-temperature model","machine-learned potential","recrystallization"],"falsifier":"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.","tokens_in":12972,"feed_emoji":"⚛️","tokens_out":3311,"duration_ms":27894,"temperature":0.7,"pith_summary":"The paper argues that the common assumption of radially symmetric ion tracks fails for low-symmetry crystals. Using multiscale simulations, it shows that beta-Ga2O3 responds to swift heavy ions with a sequence: full recovery at low energy loss, recrystallization into a metastable gamma phase at intermediate loss, and core-shell amorphous/gamma tracks at high loss. Even though energy is deposited essentially isotropically, final track shapes are anisotropic: the amorphous core is consistently narrowed along the stiff [010] direction. The simulated track sizes quantitatively match experimental measurements across a wide range of energy losses. If right, this makes elastic anisotropy a central predictor of radiation damage morphology in oxides.","feed_headline":"Elastic stiffness, not energy spread, sets ion-track shape","feed_subtitle":"Simulations show beta-Ga2O3 tracks stay narrow along the stiffest crystal direction, matching experiment.","key_machinery":"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.","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Track shape in Ga2O3 set by stiffness, not energy spread","Stiff crystal directions heal faster, shrinking ion tracks","Monoclinic Ga2O3: recovery dynamics dictate track shape","Fast healing along stiff axis narrows ion track in beta-Ga2O3","Simulations: elastic modulus, not heat spike, controls ion-track shape"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Track shape in Ga2O3 set by stiffness, not energy spread","Stiff crystal directions heal faster, shrinking ion tracks","Monoclinic Ga2O3: recovery dynamics dictate track shape","Fast healing along stiff axis narrows ion track in beta-Ga2O3","Simulations: elastic modulus, not heat spike, controls ion-track shape"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000564,"raw_usage":{"total_tokens":2536,"prompt_tokens":796,"completion_tokens":1740,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":540,"completion_tokens_details":{"reasoning_tokens":1646}},"tokens_in":540,"tokens_out":1740,"duration_ms":11044,"temperature":1.0,"reasoning_tokens":1646,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T23:27:07.347837+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}