{"id":"05225e01-ed3e-46eb-927b-2a4d7265865a","arxiv_id":"2603.04152","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A machine-learned atomic model for Ti-carbide MXenes reproduces DFT accuracy and yields the first large-scale molecular-dynamics statistics for ion irradiation damage and implantation.","lead":"MXenes are ultra-thin metal-carbide sheets used in batteries and electronics. This paper builds a fast computer model of their atoms, trained on quantum mechanics, and uses it to map what happens when ions strike the sheets.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Non-spin-polarized DFT reference is the critical soft spot: tabGAP is fitted to a constrained nonmagnetic PES, while spin-polarized Ti2C is 33 meV/atom lower and has a 0.05 Å larger lattice constant; the potential cannot represent this.","rationale":"The reader's weakest assumption identifies exactly the same soft spot: the entire near-equilibrium reference is non-spin-polarized DFT despite a known 33 meV/atom spin-polarization energy and 0.05 Å lattice-constant shift. My stress test found no fatal error in the ML methodology itself; the tabGAP appears to reproduce its target reference faithfully, and the training/validation strategy is sound given that reference. However, the claim that the potential captures the DFT energy surface and provides reliable irradiation guidelines requires that the chosen reference be adequate for the intended physical regime. Low-energy He implantation (15–40 eV) and near-threshold sputtering are precisely the observables most likely to be affected by a systematic ~1.6% lattice-constant offset. The authors are transparent about the limitation, which is why this is a conditional concern rather than a rejection, but the paper does not quantify the downstream effect on the reported statistics. The proposed test—re-running the irradiation series at the spin-polarized lattice constant—would directly determine whether the concern actually lands. If the results are unchanged, the spin issue is a caveat; if they shift, the quantitative guidelines need revision. I therefore keep the reader's CONDITIONAL verdict unchanged.","tokens_in":12944,"tokens_out":8177,"duration_ms":91225,"concrete_test":"Re-run the published He and Ti irradiation matrix for Ti2C (10,000 impacts per energy, 15 eV to 4 keV) using tabGAP but with the MXene sheet pre-expanded to the spin-polarized lattice constant a = 3.08 Å instead of the non-spin 3.03 Å, keeping all other settings identical. If sputtering yields, reflection/implantation probabilities, or coordination-defect counts shift by more than the sampling uncertainty of 10,000 impacts, the quantitative conclusions depend on the spin treatment and the paper must re-scope its accuracy claim. If they do not shift, the concern is bounded and the conditional verdict can be accepted with only a clarifying caveat.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing concern is the choice of DFT reference, not the ML methodology. Methods 4.1 and Supplemental S2 state that all training data are non-spin-polarized PBE (except isolated-atom energies), even though spin-polarized Ti2C has a ferromagnetic state 33 meV/atom lower in energy and a relaxed lattice constant of 3.08 Å versus 3.03 Å. tabGAP descriptors (two-body, three-body, EAM density) contain no spin/magnetization degrees of freedom, so the trained potential cannot represent the magnetic state or the strain-induced collapse of the moment documented in S2. The close agreement in Table 1 is therefore agreement with a constrained nonmagnetic reference, not necessarily with the physical PBE ground state. Because the irradiation simulations are initialized by relaxing at the tabGAP (non-spin) lattice constant, a systematic 0.05 Å (~1.6%) in-plane strain is built into every trajectory if the spin-polarized state is the relevant ground state. Sputtering thresholds, reflection/implantation probabilities, and damage counts are sensitive to such strains at impact energies near threshold. The authors disclose the limitation and partially correct isolated-atom limits, but a 33 meV/atom systematic reference bias is large enough to affect the quantitative 'guidelines for defect engineering' that constitute the applied half of the central claim. What is missing is any bound on how much the non-spin reference shifts the reported irradiation statistics or near-equilibrium properties.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a tabulated Gaussian approximation potential (tabGAP) for bare Ti_{n+1}C_n MXenes, trained on a diverse set of DFT structures including MXenes, non-MXene Ti–C phases, defects, and iteratively sampled disordered/recrystallized configurations. The potential is validated against DFT for lattice constants, monovacancy formation energies, shear stacking-fault profiles, phonon dispersions, and quasi-static atom drag. The authors then apply the potential to large-scale MD irradiation simulations of Ti2C and Ti3C2 sheets, reporting He and Ti ion reflection, implantation, pass-through, sputtering yields, and coordination defects as functions of impact energy from ~15 eV to 100 keV, and conclude with guidelines for defect engineering by ion irradiation. The main claims are that the ML potential is accurate and robust across a wide range of bond environments and that the irradiation statistics provide first-of-their-kind MD-based quantitative guidance.","tokens_in":13303,"tokens_out":4619,"duration_ms":49731,"significance":"If the claims hold, the paper addresses a genuine gap: despite the technological importance of MXenes, there are very few ML interatomic potentials for them, and no dedicated MD potential for Ti_{n+1}C_n irradiation studies. The approach is physically sensible, the training database is reasonably diverse, and the validation shows good agreement with DFT for several properties (lattice constants, vacancy energies, shear profiles, phonons, drag). The use of tabGAP is a strength: it is fast, transferable, and the repulsive ZBL-corrected short-range part is appropriate for collision cascades. The irradiation study is extensive — about one million individual impact simulations — and the qualitative finding of preferential Ti sputtering and the resilience/healing of the sheets is consistent with available experiments. The authors are also transparent about the non-spin-polarized DFT reference and the isolated-atom correction in Figure 5. However, the central applied claim — quantitative defect-engineering guidelines — depends on the adequacy of a constrained nonmagnetic reference, which is not quantitatively bounded. The reproducibility claim is also weakened by the missing open-reposito","major_comments":[{"comment":"The training reference is non-spin-polarized PBE except for isolated-atom energies. The manuscript itself reports that spin-polarized Ti2C is 33 meV/atom lower in energy and has a lattice constant of 3.08 Å instead of 3.03 Å. Since the tabGAP descriptors (two-body, three-body, EAM density) contain no spin/magnetization degrees of freedom, the potential cannot represent the magnetic state. Consequently, Table 1 and Figures 2–5 validate agreement against a constrained nonmagnetic PES, not necessarily the physical PBE ground state. The irradiation simulations are initialized at the tabGAP/non-spin lattice constant, so if the magnetic state is the relevant ground state, every trajectory carries ~1.6% in-plane strain. Sputtering thresholds, reflection/implantation probabilities, and damage counts are sensitive to such strain near thresholds. The authors disclose the limitation and correct iso","section":"§4.1, Supplemental S2, Table 1, Figure 5"},{"comment":"The text states 'All input parameters and training data are available in an open repository [47]', but reference [47] is 'Link to be added later.' This is load-bearing for a paper whose central product is a trained ML potential: without the exact training database, descriptor parameters, and regression settings, the potential cannot be reproduced or independently retrained. The repository link or a permanent DOI must be provided; this is a necessary condition for the reproducibility of the method.","section":"§4.2, Ref. [47]"},{"comment":"Most validation quantities are drawn from categories that are explicitly part of the training database: strained/sheared MXenes, vacancies, AIMD thermal frames, and defect configurations. The parity plots in Figure 1 show training errors, not holdout test errors. Thus the agreement in Table 1 and Figures 3–5 is largely interpolative and does not by itself demonstrate robustness outside the training distribution. The genuinely out-of-sample part is the irradiation application, but no DFT check is provided for the conditions used there (e.g., threshold displacement energies, recoil damage, or implanted configurations). Please report a holdout test RMSE and benchmark at least one irradiation-relevant quantity against DFT (for example, the displacement threshold energy along a low-index direction, or the energy of a Frenkel pair) to substantiate the extrapolation claim.","section":"§2.1 vs §2.2, Figure 1"}],"minor_comments":[{"comment":"Typo: 'with with both light and heavier ions' — delete the second 'with'.","section":"§2.3"},{"comment":"The sentence 'for Ti3C2 sheets they are much higher (Figure 5(b))' should refer to Figure 6(b), not Figure 5(b). Figure 5 shows atom drag, not irradiation probabilities.","section":"§2.3, last paragraph"},{"comment":"The caption says the Tersoff data are shifted in energy to match DFT minima. This is fine for comparing the shape of E(a), but it means the absolute Tersoff energy scale is not directly compared; the text should state this more explicitly to avoid misleading absolute-energy comparisons.","section":"Figure 2 caption"},{"comment":"The speed is quoted as '70–135 katom-steps/s' on 'a single CPU node' with a specific AMD processor. Please specify whether this is one core or multiple cores, and clarify the units (kilo atom-steps per second) so that the comparison with Tersoff is reproducible.","section":"§2.2, speed benchmark"},{"comment":"The labels 'spin' in the figure are ambiguous. It would be clearer to label the curves as 'DFT no spin-pol.', 'DFT spin-pol. isolated atom', and 'spin-corrected DFT'.","section":"Figure 5"}],"recommendation":"major_revision","confidential_remarks":"The core concern is not the ML methodology but the choice of DFT reference. The authors are transparent about the non-spin-polarized training data, but because the irradiation application is quantitative and the magnetic state is significantly lower in energy near equilibrium, a sensitivity analysis is needed to determine whether the reported defect-engineering guidelines survive. The missing data repository link is also a clear reproducibility problem. If the authors can provide the missing link and a quantitative bound on the spin-reference effect, the paper would be a strong candidate for acceptance. I do not see evidence of circularity in the fitting; the irradiation predictions are not fitted to any benchmark."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper delivers the first machine-learned interatomic potential for bare Ti_{n+1}C_n MXenes, and that alone makes it worth a serious look. It builds a 1522-structure training set with an iterative self-learning loop, validates against DFT for lattice constants, vacancy energies, shear profiles, phonons, and atom drag, and then runs roughly one million irradiation impacts to map sputtering, reflection, implantation, and defect production as functions of ion energy. The irradiation statistics are genuinely new — nothing of that scope exists for MXenes. The author is also transparent about the main limitation: all training data are non-spin-polarized PBE, and the paper itself reports that spin-polarized Ti2C sits 33 meV/atom lower and relaxes to a 0.05 Å larger lattice constant. That is not hidden; it is in Methods 4.1 and Figure S3.\n\nThe soft spots are real but proportionally manageable. The spin issue is the most consequential. The potential cannot represent the magnetic state, so it is a potential for a constrained nonmagnetic PES, not for the DFT ground state of Ti2C. The 0.05 Å lattice constant shift means every irradiation trajectory starts from a 1.6% in-plane strain if the spin-polarized state is the physical one. Sputtering thresholds near threshold energies could shift. The author discloses this, but does not bound the effect. A few test runs at the spin-polarized lattice constant would settle it. That is the main thing I would ask for. The validation overlap with training categories is a lesser concern — the melt-quench tests provide some evidence of extrapolative robustness. And the data repository link is still a placeholder, which blocks reproducibility right now.\n\nOn balance, the central claim holds up with caveats. This is a useful tool and the first MD-level guidance for defect engineering of Ti MXenes, but the accuracy statement should be read as \"accurate to the nonmagnetic DFT reference.\" The irradiation probabilities are probably in the right ballpark, but the near-quantitative guidelines need the sensitivity check.\n\nWho is this for? Anyone doing MD simulations of MXene mechanics or radiation effects, especially experimental groups planning ion-irradiation defect engineering. It deserves a serious referee — I would accept it with major revision requests: supply the data, and add a sensitivity analysis on the spin-polarized reference. Not a takedown; a good paper that needs a tighter qualification on its headline claim.","headline":"First ML potential for bare Ti MXenes, solid and useful, but the non-spin-polarized DFT reference leaves a real, unbounded caveat in the near-equilibrium regime.","tokens_in":724,"tokens_out":766,"would_cite":true,"duration_ms":28942,"reading_group":"yes","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 develops a machine-learned interatomic potential for titanium carbide MXenes that is fast and accurate enough for molecular dynamics, and uses it to produce the first quantitative ion-irradiation statistics—sputtering yields, ref","keywords":["MXenes","machine-learned interatomic potential","tabGAP","ion irradiation","sputtering","defect engineering","titanium carbide","molecular dynamics"],"falsifier":"Measure the equilibrium in-plane lattice constant of a bare single-layer Ti2C MXene using spin-polarized DFT or a high-level method such as quantum Monte Carlo: the paper's potential gives 3.03 Å, while spin-polarized DFT gives 3.08 Å. If a reliable reference confirms 3.08 Å (or shows the potential's relaxed geometry is off by more than ~0.03 Å), the non-spin reference is inadequate and the potential's near-equilibrium description—and all irradiation simulations starting from its relaxed sheets—would be biased.","tokens_in":12826,"feed_emoji":"⚛️","tokens_out":4045,"duration_ms":41327,"temperature":0.7,"pith_summary":"The paper aims to show that a machine-learned interatomic potential of the tabGAP type can be trained on a diverse set of density functional theory data to accurately describe bare Ti_{n+1}C_n MXenes across a wide range of bond distances and atomic environments. If true, this would remove a major bottleneck: the lack of reliable interatomic potentials has previously forced researchers to use either expensive ab initio dynamics or inaccurate classical potentials for molecular dynamics studies of MXenes, especially for ion irradiation. The author applies the trained potential to simulate roughly one million individual He and Ti ion impacts on single-layer Ti2C and Ti3C2 sheets, extracting statistics for sputtering, ion reflection, implantation, and defect production as functions of impact energy. These results provide concrete, simulation-based guidelines for experimental defect engineering of MXenes by ion irradiation, and they confirm experimentally observed preferential sputtering of Ti at higher energies while revealing a low-energy regime where C sputters more easily under light-ion bombardment.","feed_headline":"ML potential opens Ti-carbide MXenes to ion-irradiation simulation","feed_subtitle":"Trained on 1,522 DFT structures, the tabGAP potential predicts sputtering and defect stats for He and Ti impacts.","key_machinery":"The central object is the tabulated Gaussian approximation potential (tabGAP), a low-dimensional machine-learned interatomic potential that combines a two-body distance descriptor, a three-body permutation-invariant descriptor, and a scalar embedded-atom-method density descriptor, with predictions tabulated on 1D and 3D spline grids for computational efficiency. It is supplemented by a Ziegler–Biersack–Littmark screened-Coulomb repulsive pair potential fitted to all-electron DFT data, which handles close-range collisions during ion irradiation. The tabGAP formalism carries the argument by delivering near-DFT accuracy at a computational cost only 20–40% slower than an analytical Tersoff poten","core_discovery":"The central claim is that a tabulated Gaussian approximation potential (tabGAP) trained on 1,522 DFT-generated structures containing 24,150 atoms—covering strained/sheared MXenes, defective and disordered sheets, non-MXene Ti–C phases, and iterative-learning structures—reproduces the DFT energy surface for Ti_{n+1}C_n MXenes with root-mean-square errors of 4.5 meV/atom and 0.19 eV/Å in forces on MXene structures. Validation against DFT shows accurate lattice constants (within ~0.01 Å), vacancy formation energies (within ~0.2 eV), shear stacking-fault profiles, phonon dispersions, and quasi-static atom-drag curves. Using this potential, the paper reports that for He irradiation of Ti2C and Ti","pith_inferences":["Because all training data except isolated-atom references were generated without spin polarization, the potential's near-equilibrium description of bare MXenes may be biased: spin-polarized DFT lowers Ti2C energy by 33 meV/atom and shifts the lattice constant from 3.03 Å to 3.08 Å, so if spin polarization is physically important for bare sheets, the potential's relaxed geometry—and any irradiation","The use of purely repulsive He–Ti and He–C pair potentials likely underestimates low-energy He implantation probabilities; adding an attractive interaction (as the paper acknowledges) would raise the chance of He sticking to the sheet, which is directly testable by comparing to low-energy ion-scattering experiments.","The iterative 'learn from its mistakes' loop that eliminated unphysical dense carbon clusters in high-temperature MD suggests a general, inexpensive way to detect extrapolation failures in ML potentials for 2D materials: run high-temperature or damage simulations and scrutinize any energetically over-stable but structurally anomalous configurations, then retrain.","The reported double-peak in C sputtering yield for He on Ti3C2, attributed to correlated two-atom sputtering events, implies that standard single-sputtering models might poorly describe light-ion damage in multilayer MXenes; this could be probed by analyzing impact-parameter-resolved sputtering events in the simulations."],"forward_implications":["Provides the first molecular-dynamics-based quantitative guidelines for ion-beam defect engineering of Ti_{n+1}C_n MXenes, including sputtering yields, ion reflection, implantation probabilities, and damage extent as functions of impact energy.","Confirms experimentally observed preferential sputtering of Ti under heavy-ion and high-energy light-ion irradiation, and predicts a low-energy light-ion regime where C sputtering is enhanced due to more efficient energy transfer to the lighter C atoms.","Shows that single-layer Ti2C and Ti3C2 sheets remain intact and self-heal after heavy Ti impacts at energies up to 100 keV, consistent with experimental stability observations and suggesting that irradiation-induced defect engineering is a viable route for MXenes.","Demonstrates a repeatable training strategy—mixing manually constructed, ab initio molecular dynamics, and iterative-learning structures—that the author argues can be transferred to develop ML potentials for other MXene compositions and surface terminations.","Quantifies low but non-negligible He implantation probabilities (0.0001–0.0005 for 15–40 eV He on Ti2C), relevant for understanding ion-beam doping of MXenes."],"fun_headline_variants":["ML potential simulates ion irradiation damage in Ti-carbide MXenes","Fast tabGAP force field opens MXene irradiation studies","Ion impacts on MXenes modeled with machine-learned potential","Ti2C MXene sputtering stats from ML interatomic potential","Machine-learned potential predicts MXene defect engineering"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire reference dataset for training (except isolated atoms) was computed with non-spin-polarized density functional theory, yet near equilibrium spin-polarized Ti2C is lower in energy by 33 meV/atom and has a ~0.05 Å larger lattice constant; if non-spin DFT is not an adequate reference for near-equilibrium MXene energetics, the potential's accuracy fails in precisely the regime used to relax the sheets before irradiation.","fun_headline_variants_meta":{"raw":{"variants":["ML potential simulates ion irradiation damage in Ti-carbide MXenes","Fast tabGAP force field opens MXene irradiation studies","Ion impacts on MXenes modeled with machine-learned potential","Ti2C MXene sputtering stats from ML interatomic potential","Machine-learned potential predicts MXene defect engineering"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000217,"raw_usage":{"total_tokens":1252,"prompt_tokens":705,"completion_tokens":547,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":449,"completion_tokens_details":{"reasoning_tokens":473}},"tokens_in":449,"tokens_out":547,"duration_ms":5244,"temperature":1.0,"reasoning_tokens":473,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T18:53:23.739247+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the equilibrium in-plane lattice constant of a bare single-layer Ti2C MXene using spin-polarized DFT or a high-level method such as quantum Monte Carlo: the paper's potential gives 3.03 Å, while spin-polarized DFT gives 3.08 Å. If a reliable reference confirms 3.08 Å (or shows the potential's relaxed geometry is off by more than ~0.03 Å), the non-spin reference is inadequate and the potential's near-equilibrium description—and all irradiation simulations starting from its relaxed sheets—would be biased.","supporting_citations":[],"review_version":1}