{"id":"d32177cf-2854-4bd3-a9b7-4828cdfcb77d","arxiv_id":"2608.09791","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Simulations with universal machine-learned and polarizable force fields reproduce the known electrolyte dependence of SEI chemistry, with fluorinated solvents yielding LiF-rich interphases.","lead":"This paper combines two universal force fields to simulate how lithium metal and different electrolytes form solid-electrolyte interphases, and finds that fluorinated solvents produce lithium-fluoride-rich layers while carbonate solvents produce organic-rich layers. It is a step toward using computer simulation to screen electrolytes for next-generation batteries without years of trial-and-error experiments.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified beyond the reader's conditional; the central trend replication is credible, but the universal-prediction claim needs an out-of-distribution reactivity check.","rationale":"The reader's weakest_assumption is exactly the weakest point: the universal MLFF's reactivity beyond training distribution is unvalidated. The paper itself contains no direct evidence of this assumption's correctness, so the conditional verdict is appropriate. The paper has independent support in the form of replicating experimental trends, which counts as real evidence, but the lack of error bars, the single 200 ps trajectory per system, and the absence of code/data prevent a stronger verdict. The correct response is to keep the conditional verdict and request the specific transferability check. I recommend UNCHANGED because the reader's conditional captures the appropriate level of confidence; no additional concrete flaw emerged in the stress test.","tokens_in":8355,"tokens_out":1327,"duration_ms":10736,"concrete_test":"Run short (50-100 ps) reactive MD on the LiFSI/FSA 1:3.5 and LiFSI/EC/DMC/FEC interfaces using a second independent universal MLFF (e.g., a MACE-MP or CHGNet-based potential) or a QM/MM reference for the same interface. Compare (a) the LiF/organic mole-fraction ordering and (b) the dominant decomposition products to the UMA-s-1p2 results. If the qualitative trend is preserved, the central prediction is robust; if not, the claim of universal prediction is specific to one model and needs revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that combining ByteFF-Pol and UMA-s-1p2 enables universal predictive simulation of interphase chemistry. The abstract's strongest claim is the replication of the experimental LiF-rich vs organic-rich trend across four chemically distinct electrolytes. The trend is physically plausible and consistent with known electrolyte chemistry. The weakest assumption is that UMA-s-1p2, trained on OC25, reliably describes spontaneous bond-breaking/forming at the Li/electrolyte interface for chemistries beyond its training distribution. This is a genuine risk, but the paper's design partially mitigates it: OC25 includes solid-liquid interfaces, and the observed chemistry qualitatively matches experiment. The paper does not validate UMA-s-1p2 against quantum chemistry for the specific reaction pathways, and no error bars or multiple seeds are reported, so quantitative claims like SEI thicknesses and mole fractions are unsupported. Still, the load-bearing concern is not internal inconsistency; it is a strong but unproven extrapolation. A direct test would be to run the same reactions with a QM/MM or DFTB approach on small model systems to check whether the MLFF's predicted pathways and barrier heights are accurate. Another critical check is the sensitivity of the results to the MLFF's training set: does an older or different universal potential (e.g., MACE-MP or CHGNet) produce the same qualitative trend? This settles whether the protocol is universal or specific to UMA-s-1p2.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a molecular dynamics protocol that combines a commercially developed polarizable force field (ByteFF-Pol) for bulk electrolyte equilibration with the UMA-s-1p2 universal machine-learned interatomic potential for reactive interfacial dynamics. Six electrolyte compositions are simulated on a Li(001) slab for 200 ps each, and the authors report SEI thickness, molecular composition, LiF density, and dominant reaction pathways for LiFSI/FSA, LiFSI/F5DEE, LiFSI/SFL/TOL, and LiFSI/EC/DMC/FEC. The central claim is that the protocol reproduces the qualitative experimental trend that fluorinated solvents form LiF-rich interphases whereas conventional carbonate solvents form organic-rich interphases, establishing a route toward universal predictive SEI simulation.","tokens_in":8636,"tokens_out":4616,"duration_ms":42345,"significance":"If the claims hold, the protocol would be a valuable high-throughput screening tool for electrolyte engineering, because it uses pre-trained universal potentials that are independent of the target SEI outcome, avoiding the circularity often present in chemistry-specific fitted potentials. The qualitative reproduction of the LiF-rich versus organic-rich trend across four chemically distinct electrolyte classes is encouraging, and the explicit description of the construction, equilibration, and production protocols supports reproducibility. The main significance is conditional, however: the reliability of the reactive simulation rests entirely on the extrapolative accuracy of UMA-s-1p2 for bond-breaking and bond-forming at the Li/electrolyte interface, and the quantitative SEI properties are reported without uncertainty estimates. The paper would be substantially strengthened by adding direct validation of representative reaction pathways against quantum chemistry and by reporting multiple independent trajectories.","major_comments":[{"comment":"The load-bearing assumption that UMA-s-1p2 correctly describes spontaneous bond-breaking and bond-forming at the Li/electrolyte interface is not validated. The manuscript states that universal MLFFs 'maintain quantum level accuracy' and that OC25 was selected because it includes solid-liquid configurations, but no evidence is provided that the model is accurate for reaction pathways and barrier heights of the specific molecules studied (FSA, F5DEE, SFL, EC, DMC, FEC). Because the entire central claim depends on reactions emerging from the MLFF, the authors should add a direct test: compare the predicted dominant pathways (e.g., F5DEE defluorination order, FSA C-N bond stability, EC ring-opening) against DFT or QM/MM calculations on small model systems, or at least perform sensitivity tests with an alternative universal MLFF (e.g., MACE-MP or CHGNet) to see whether the same qualitative trend is reproduced. Without such a check, the simulated SEI compositions could reflect force-field artifacts rather than chemistry.","section":"Methods (Interfacial molecular dynamics) and Results (Molecular reaction pathways)"},{"comment":"All SEI thicknesses, mole fractions, and LiF densities are derived from a single 200 ps trajectory per composition, with no independent seeds, no block averaging, and no reported statistical uncertainty. The paper makes quantitative claims such as 'thickness increases from 23.39 Å at 1:10 to 28.33 Å at 1:6 and 30.85 Å at 1:3.5' and compares values like 7.63, 9.77, and 9.35 Å across electrolytes; these differences may be within noise given the small system size and short timescale. The authors should provide error bars from multiple independent trajectories or time-block decomposition, and should state whether the observed differences are statistically significant. This is essential to support the 'systematic increase' and 'lowest/highest' statements that the paper makes.","section":"Results, Table 4 and Figure 2"},{"comment":"The reaction mechanisms in Figure 4 are inferred from a single reactive trajectory, and the text presents them as the predominant pathways (e.g., 'the C–F bonds at the –CF2 end undergo reduction first, followed by those at the –CF3 end' for F5DEE; 'no C-N bond breaking is observed' for FSA). The generalizability of these mechanistic claims is unclear without either multiple independent trajectories showing reproducibility or quantum-chemical calculation of the relevant barriers. The paper should at minimum state the number of times each pathway was observed across the trajectory and how the 'predominant' pathway was selected, and ideally confirm the key pathways against DFT.","section":"Results (Molecular reaction pathways, Figure 4)"},{"comment":"The abstract claims that the combined force fields allow the authors to 'simulate interphasial chemistry across chemically diverse electrolyte formulations' and describes a 'route toward universal and high-throughput predictive simulation of interphases.' The conclusion uses 'could lead to a universal tool,' which is more measured. Given that only six compositions are tested, all of them LiFSI salts and all lithium-metal anodes, the evidence is too narrow to support the adjective 'universal' as applied to the simulation method itself. The authors should either expand the validation to include substantially different chemistries (e.g., different salts, anode materials, or polymer electrolytes) or explicitly temper the abstract's universal claim to reflect the demonstrated domain, so that the title and abstract do not overstate the current scope.","section":"Abstract and Conclusion"}],"minor_comments":[{"comment":"The caption lists panels (a), (b), (c), (e), (f), and (g) but omits panel (d) and then describes 'representative structures' in panel (d) and (h); the panel labels in the caption and the figure should be reconciled.","section":"Figure 3 caption"},{"comment":"The header 'Molecular ratio of reacted molecules' is used inconsistently: the first column of Table 4 lists the electrolyte name, then a ratio expression such as 'FSI:EC:DMC:FEC = 1:1.50:0.58:0.58'; the formatting should be unified and the units clarified, because the column mixes a ratio and a count-like quantity.","section":"Table 4"},{"comment":"The abbreviation 'UMA' is used in the Introduction but the full model name 'UMA-s-1p2' and its training data (OC25) are first defined only in Methods; defining it at first use would make the paper easier for a nonspecialist to follow.","section":"Introduction"},{"comment":"The paper states that bulk equilibration used 'constrained anisotropic isothermal-isobaric ensemble (NPT)' with lateral cell vectors fixed and the z length fluctuating; the text also says 'Monte Carlo anisotropic barostat at 1 atm, 298 K.' The description is understandable but the fixed-lateral-cell variant of NPT is unusual and should be explicitly justified as appropriate for later matching to the Li slab, which is fixed in x and y.","section":"Methods (Electrolyte construction)"}],"recommendation":"major_revision","confidential_remarks":"This manuscript falls within the scope of the journal and addresses an important problem. The primary concern is the unvalidated extrapolative reactivity of the universal MLFF, which is central to the paper's claims; the requested quantum-chemical validation or cross-MLFF sensitivity test is feasible within the project's scope and would address the concern. The lack of error bars is also fixable. If the authors are unwilling or unable to add those tests, the paper would be more appropriately positioned as a demonstration of a promising but unvalidated protocol rather than as a predictive simulation study. The manuscript's heavy reliance on references to tools from the authors' own groups (e.g., Bytesize, FAIR-Chem) is not itself a concern, but the novelty relative to prior SEI-simulation work is somewhat incremental. I do not see any evidence of misconduct or duplication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here is my read. The paper does something genuinely new: it pairs two general-purpose models, ByteFF-Pol for the liquid and UMA-s-1p2 for the reactive interface, and runs a single protocol across chemically different electrolytes to see what kind of SEI forms spontaneously. Reproducing the known LiF-rich versus organic-rich trend across FSA, F5DEE, carbonate, and SFL/TOL is real evidence that the approach has signal, not just noise. The reaction-pathway analysis (all five F on F5DEE ending as LiF, TOL staying inert) is chemically sensible and consistent with prior work. The protocol is described carefully enough to reproduce.\n\nThe soft spot is the scope of the claim. 'Universally predict' is supported by six formulations, one 200 ps trajectory per system, and no error bars or repeated seeds. The MLFF's ability to describe spontaneous bond-breaking at the Li interface for chemistries beyond its training distribution is assumed, not validated. There is no check against quantum chemistry for the specific pathways reported. That is the load-bearing assumption. Without it, readers cannot tell whether the LiF is real chemistry or a force-field artifact. No code or data accompanies the paper, which slows independent verification. These are fixable: run multiple seeds, compare selected pathways with QM/MM or DFTB, and test whether another universal potential gives the same qualitative ranking.\n\nThe stress-test note matches my view. The central trend is plausible; the universality claim is unproven. I would also note the SEI thickness values are likely sensitive to the 200 ps cutoff, since growth has not visibly plateaued for all systems. That does not destroy the qualitative ranking, but it does undercut quantitative statements.\n\nWho this is for: people working on electrolyte screening and MLFF applications to battery interfaces. It is a useful proof-of-concept, not yet a validated predictive tool. I would send it to a serious referee because it deserves expert time, but the revision should require reactivity validation and uncertainty quantification. I would not cite it as evidence of universal prediction, though I might cite it as an early example of combined universal force fields for SEI simulation. Reading group: maybe, as a discussion of what counts as prediction in ML-driven simulation.","headline":"A credible first pass at universal-MLFF SEI simulation that reproduces known trends, but the universality claim outruns the evidence; send to review with a demand for reactivity validation and uncertainty quantification.","tokens_in":714,"tokens_out":2160,"would_cite":false,"duration_ms":34694,"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 claims that pairing a polarizable force field with a universal machine-learning force field predicts the solid–electrolyte interphase chemistry on lithium metal from electrolyte composition alone.","keywords":["solid electrolyte interphase","lithium metal anode","machine learning force field","reactive molecular dynamics","electrolyte decomposition","LiF-rich interphase","force field simulation","battery electrolyte design"],"falsifier":"Take one electrolyte outside the force field's training distribution, run the same 200 ps interfacial protocol, and compare the predicted interphase composition and reaction products with depth-profiled X-ray photoelectron spectroscopy or cryogenic electron microscopy of a cell cycled under otherwise identical conditions. If the observed LiF fraction or organic-to-inorganic ratio disagrees with the simulation, the claim of universal predictive power fails.","tokens_in":8190,"feed_emoji":"🔋","tokens_out":5666,"duration_ms":45068,"temperature":0.7,"pith_summary":"This paper claims that the chemistry of the solid–electrolyte interphase that forms on lithium metal can be predicted from the electrolyte composition alone, using atomistic simulation instead of trial-and-error experiments. It combines a transferable polarizable force field with a universal machine-learning force field and runs reactive molecular dynamics at a lithium–electrolyte interface for four chemically distinct electrolyte families. The simulations reproduce the established experimental pattern: fluorinated solvents such as FSA and F5DEE yield LiF-rich interphases, while carbonate-based electrolytes yield organic-rich interphases. If the claim holds, electrolyte design for lithium-metal batteries could be screened computationally before any cell is built.","feed_headline":"Simulation predicts battery interphase chemistry from electrolyte alone","feed_subtitle":"Pairing two force fields reproduces LiF-rich interphases from fluorinated solvents and organic-rich ones from carbonates.","key_machinery":"The load-bearing instrument is a reactive molecular-dynamics protocol in which a graph-neural-network polarizable force field equilibrates the bulk electrolyte, and a universal machine-learning interatomic potential trained on solid–liquid interface configurations propagates the interfacial dynamics without predefined reaction templates. A 10x10 lithium slab with a fixed bottom layer is placed against the equilibrated electrolyte, with a harmonic wall preventing electrolyte atoms from escaping into vacuum; 200 ps trajectories are then analyzed by counting reacted molecules, locating LiF along the interface, and clustering decomposition products. Because bond breaking and bond formation emerge from the learned potential, the same protocol can be applied to an arbitrary electrolyte without reparametrizing a reactive force field for each chemistry.","core_discovery":"The central claim is that spontaneous interfacial decomposition can be captured by combining a polarizable force field for realistic bulk electrolyte packing with a universal machine-learning force field for reactive dynamics, requiring no chemistry-specific reactive potential. Across six formulations, the simulations generate distinct interphase morphologies and compositions in 200 ps of dynamics. The computed LiF fractions and thicknesses track the known experimental ordering: F5DEE gives the highest LiF fraction, FSA electrolytes form thick inorganic-rich interphases, and the carbonate system gives the thinnest, most organic-rich interphase. The simulation also resolves molecular pathways, including stepwise defluorination of F5DEE and dominance of FSA as a fluorine donor, and shows that higher salt concentration thickens the SEI by increasing anion participation while leaving LiF density nearly unchanged.","pith_inferences":["I infer that the authors intend this two-force-field protocol, not the specific 200 ps trajectories, as the transferable contribution; the quantitative SEI thicknesses should be read as trends rather than converged values.","An extension the paper leaves implicit is that the same combination could be pointed at other reactive interfaces, such as catalytic or electrosynthetic surfaces, by swapping the electrode slab and the electrolyte.","A testable consequence not stated in the paper is that repeating the simulation with different initial packings or longer trajectories should preserve the compositional ordering of the six electrolytes, even if absolute thicknesses drift."],"forward_implications":["If the protocol is correct, the solid–electrolyte interphase from a new electrolyte can be predicted by simulation before any electrochemical cell is assembled.","The method yields a mechanistic explanation for why fluorinated solvents, and higher salt concentrations, favor inorganic LiF-rich interphases rather than organic decomposition products.","The same simulation protocol could be applied to other electrode surfaces and interphase-forming processes without changing the force fields.","The computed reaction counts and LiF distributions provide direct molecular-level input for AI-driven electrolyte screening pipelines."],"supporting_citations":[{"why":"supplies the graph-neural-network polarizable force field used to equilibrate the bulk electrolytes","marker":"[12]"},{"why":"provides the harmonic-wall boundary treatment and the precedent for machine-learning-force-field molecular dynamics of SEI formation","marker":"[13]"},{"why":"supplies the universal machine-learning interatomic potential that carries the reactive dynamics","marker":"[16]"},{"why":"supplies the solid–liquid interface dataset whose training domain contains the reacting species","marker":"[24]"},{"why":"provides the experimental evidence that FSA acts as a strong fluorine donor, which the simulations reproduce","marker":"[29]"},{"why":"provides experimental characterization of F5DEE electrolyte interphases that the simulations are compared with","marker":"[30]"}],"fun_headline_variants":["Dual force fields predict Li metal interphase chemistry","Simulated SEI composition matches experiments without reactive potentials","Predicting interphase chemistry from electrolyte alone","Machine learning plus polarizable field predicts SEI chemistry","Spontaneous SEI reactions simulated for diverse electrolytes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method stands or falls on the assumption that the machine-learning force field correctly describes bonds breaking and forming at the lithium–electrolyte interface for these particular molecules, even though it was not trained specifically on these reaction pathways, and the paper does not validate that assumption against quantum-chemistry calculations for the claimed reactions.","fun_headline_variants_meta":{"raw":{"variants":["Dual force fields predict Li metal interphase chemistry","Simulated SEI composition matches experiments without reactive potentials","Predicting interphase chemistry from electrolyte alone","Machine learning plus polarizable field predicts SEI chemistry","Spontaneous SEI reactions simulated for diverse electrolytes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000204,"raw_usage":{"total_tokens":1368,"prompt_tokens":904,"completion_tokens":464,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":520,"completion_tokens_details":{"reasoning_tokens":391}},"tokens_in":520,"tokens_out":464,"duration_ms":5008,"temperature":1.0,"reasoning_tokens":391,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:40:39.113507+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take one electrolyte outside the force field's training distribution, run the same 200 ps interfacial protocol, and compare the predicted interphase composition and reaction products with depth-profiled X-ray photoelectron spectroscopy or cryogenic electron microscopy of a cell cycled under otherwise identical conditions. If the observed LiF fraction or organic-to-inorganic ratio disagrees with the simulation, the claim of universal predictive power fails.","supporting_citations":[{"cited_title":"UMA: A family of universal models for atoms","cited_arxiv_id":null,"evidence_quote":"supplies the universal machine-learning interatomic potential that carries the reactive dynamics"}],"review_version":1}