{"id":"477a0248-cfde-4b53-8bc9-ccfa35378eca","arxiv_id":"2412.03541","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Specialized graph neural network potentials trained on small DFT datasets qualitatively reproduce the temperature dependence of CaTiO3's dielectric tensor and its orthorhombic-to-cubic phase transition in molecular dynamics.","lead":"The authors train two machine learning force fields on a few hundred DFT snapshots of perovskite crystals and run molecular dynamics to compute how the dielectric tensor of CaTiO3 changes with temperature. The result, a decreasing dielectric response with heating and an orthorhombic-to-cubic phase transition, suggests that small specialized machine-learned potentials could replace costly DFT for such studies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"High-temperature phase-transition and dielectric results rest on GNN extrapolation far outside the 300–500 K training window; without high-T DFT validation, the agreement could be coincidental.","rationale":"I read the paper as a demonstration that a small specialized equivariant-GNN force field can produce finite-T MD for CaTiO3 that qualitatively tracks the experimental dielectric temperature dependence and a structural phase transition. The reader's verdict CONDITIONAL is justified. I find the reader's weakest_assumption (ionic-only dielectric model, Eqs. 1-2) real but not the most load-bearing: the missing clamped-ion electronic contribution to the dielectric constant is nearly temperature-independent at the few-percent level, so it would shift values but not destroy the qualitative trend. A more decisive risk is the temperature extrapolation: the training set is confined to 300-500 K, yet the key results — the 1280 K transition and the high-T dielectric slope — are computed at temperatures far above it. The model is never validated on high-T configurations; the BTO soft-mode test is zero-temperature and on a different compound. If the PES is inaccurate in this extrapolation regime, the simulated phase transition and dielectric curve could be artifacts. I therefore recommend keeping the conditional verdict, with one additional explicit condition: validate the GNN predictions on high-T MD snapshots against DFT (e.g., force MAE vs T, or retrain-and-compare). I partially agree with the reader because their weakest_assumption and rationale already gesture at extrapolation, but the highlighted weakest point differs.","tokens_in":5877,"tokens_out":12961,"duration_ms":139051,"concrete_test":"Extract 50 snapshots from the production MD at T = 1200 K and 1500 K, compute PBEsol+U forces and Born charges with DFT, and compare against MACE and Equivar predictions. If the per-atom force MAE exceeds the 300–500 K cross-validation MAE by more than a factor of two, or if the Born-charge MAE worsens by more than 2×, the high-T phase-transition and dielectric-trend claims are not supported; if errors stay within 2×, the extrapolation concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The datasets (A–D) are ab initio MD snapshots at 300–500 K, so the MACE force field and Equivar Born-charge model are trained only in that narrow window. Production MD then runs to 1800 K for the phase-transition scan (Fig. 4) and up to at least ~1000 K for the dielectric tensor (Fig. 5). At these temperatures the system explores tilt patterns, bond lengths and coordination distortions that are not represented in the training set and are never checked against DFT. The only transferability evidence offered is a zero-temperature soft-phonon path in BaTiO3 (Fig. 2b), which does not cover thermal disorder in CaTiO3. Therefore the observed orthorhombic-to-cubic switch at ~1280 K (without the tetragonal phase) and the decreasing dielectric trend could be artifacts of the learned potential rather than properties of CaTiO3. This is load-bearing because the central claim is qualitative agreement with experiment; if the high-T PES is wrong, the agreement is not predictive. The paper's own caveat that the phase transition error 'could also be due to the approximated density functional' does not address the additional extrapolation error from the 300–500 K training window.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports machine-learned interatomic potentials based on the MACE and Equivar equivariant graph neural network architectures, trained on small DFT datasets (178-502 ab initio MD snapshots) for CaTiO3, BaTiO3, and CaZrO3, including energies, forces, stresses, and Born effective charge tensors. The potentials are used in finite-temperature molecular dynamics simulations with an external electric field coupled to Born charges (Eqs. 1-2), and the resulting polarization is used to extract a temperature-dependent dielectric tensor. In addition, NpT heating runs are used to study structural phase transitions in CaTiO3. The central claim is qualitative agreement with experiment for the temperature dependence of the dielectric tensor and for the orthorhombic-to-cubic structural transition.","tokens_in":6162,"tokens_out":7923,"duration_ms":86481,"significance":"If the claims hold, the paper would demonstrate that a specialized equivariant GNN trained on roughly 178 snapshots can capture subtle temperature-dependent dielectric and phase behavior in a perovskite oxide, and that the specialized model outperforms a general-purpose foundation model on this system. The paper includes useful cross-validation and direct comparisons with DFPT Born charges, which are strengths. The practical relevance for finite-temperature property prediction is clear. However, several gaps in the dielectric model, the lack of error estimates, and the unvalidated extrapolation to temperatures far above the training window currently prevent the paper from fully supporting its central qualitative claim.","major_comments":[{"comment":"The dielectric tensor shown in Fig. 5 is the ionic (lattice) contribution only, as the caption states, but the text refers to it as the 'dielectric tensor' and compares it with experimental total values. The zero-temperature DFPT points are conventionally total static dielectric constants; if so, the comparison mixes the ionic-only result with the total experimental/DFPT values. The omission of the electronic (clamped-ion) contribution is not merely a density-functional limitation but a model restriction. Please add the electronic contribution (e.g., from a DFPT epsilon_infinity) and compare consistently, or explicitly and consistently restrict all claims and comparisons to the ionic contribution.","section":"Methods, Eqs. (1)-(2); Fig. 5"},{"comment":"No statistical uncertainties are shown in Fig. 5, and the Methods do not state how many field strengths or independent MD runs were used for each dielectric component. If the tensor is obtained from the slope of P(E), multiple field points and an error estimate are required; if it is obtained from a single finite field, the linear-response assumption should be justified. Without this information, the claim that the simulations reproduce the decreasing trend and shape of the experimental curve is not quantitatively supported.","section":"Fig. 5 and Methods"},{"comment":"The training datasets are ab initio MD snapshots in the 300-500 K range, while the production runs reach 1000-1800 K. The phase-transition simulation switches from orthorhombic to cubic at approximately 1280 K and does not visit the tetragonal phase, and the dielectric simulations also extend beyond the training temperatures. No validation against DFT energies or forces for high-temperature structures is provided. As a result, the agreement with experiment could be an artifact of the learned potential's extrapolation rather than a property of CaTiO3. Please add direct validation (e.g., single-point DFT on snapshots at several temperatures above 500 K) or otherwise demonstrate that the relevant tilt instabilities and thermal disorder are correctly described at high temperature.","section":"Fig. 4 and dataset description"},{"comment":"In Eqs. (1)-(2), the Born charges Z* are predicted by Equivar as functions of the local environment and therefore change during the MD simulation. The force in Eq. (1) is not the negative gradient of the energy in Eq. (2) unless Z* is independent of the displacements, because the derivative of Z* with respect to the displacement contributes an additional term. The manuscript does not state whether Z* is held fixed during the field-coupled MD or whether this derivative term is neglected. Since the field-induced forces are central to the dielectric response, please clarify the implementation and, if the term is neglected, estimate its magnitude.","section":"Methods, Eqs. (1)-(2)"},{"comment":"The phase-transition claim in Fig. 4 is difficult to evaluate because no order parameter or structural metric is defined. The text says the system 'switched from the orthorhombic to cubic phase and back' before adopting the cubic structure, but without a precise definition of phase identity (e.g., octahedral tilt angles or space-group assignment) and without discussion of hysteresis or heating-rate effects, the reported transition temperature of about 1280 K is not well characterized. Please provide the order-parameter curves and address finite-size and heating-rate dependence.","section":"Fig. 4"}],"minor_comments":[{"comment":"The final sentence of the Conclusions says 'molecular dynamics simulations at finite temperatures and infinite electric fields are performed'; this should read 'finite electric fields.'","section":"Conclusions"},{"comment":"The sentence 'Hubbard U correction of 3 eV was applied to d electrons in transition metals, except for Ti, where U=0 was used' is confusing because Ti is a transition metal; please specify that the U=3 eV applies to Zr (or otherwise clarify the treatment for each element).","section":"Methods"},{"comment":"The caption refers to 'Dynamic charges'; for consistency with the text, please use 'Born effective charges' or 'dynamical charges' explicitly.","section":"Fig. 1 caption"},{"comment":"The statement that 'MP-0 FF underestimates the force constants' needs a definition of how force constants are extracted from the P(E) curves, and the temperature at which the P(E) data are computed should be stated.","section":"Fig. 3"},{"comment":"The text says the simulations reproduce the shape of the experimental curve, but the experimental data are not shown in Fig. 5; please add the experimental curve or specify the source and values used for the comparison.","section":"Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The novelty relative to the authors' prior work (refs. 11, 19, 20) should be clarified in the revision; the current text does not fully delineate what is new beyond applying the previously developed methods to CaTiO3. This is a scoping concern rather than a technical flaw. The manuscript fits the journal's scope, but the dielectric-model completeness and the extrapolation validation need to be addressed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a practical and mostly honest paper, but the central claims are softer than the abstract makes them sound. The genuinely new thing is the workflow: train a specialized MACE force field and an Equivar Born-charge predictor on a few hundred MD snapshots, then run finite-field MD to get a temperature-dependent dielectric tensor and a structural phase transition in CaTiO3. The small-data learning results are convincing, especially the comparison with the foundation model MACE-MP-0, and the authors are upfront about zero-point motion and DFT limitations. That part deserves credit.\n\nThe soft spots are real and they matter. First, the models are trained on 300–500 K ab initio MD snapshots, then run up to 1800 K for the phase transition and to roughly 1000 K for the dielectric tensor. There is no high-temperature DFT validation at all, so the orthorhombic-to-cubic switch at ~1280 K and the decreasing dielectric trend could be artifacts of the learned potential. The paper's caveat that the transition error might come from the density functional does not address this extrapolation gap. Second, the dielectric computation uses Eqs. 1–2, which couple the field only to Born charges and read out polarization from ionic displacements. That gives the ionic contribution, not the total dielectric response, and Fig. 5's caption says 'ionic part' while the text compares to experimental total values. That is a mismatch, not a minor wording issue. Third, there are no error bars on the MD-averaged dielectric tensor, so the claimed agreement with the experimental curve is hard to assess statistically. The abstract also overstates the phase-transition result: the simulation shows a single orthorhombic-to-cubic transition, while experiments show two transitions; the paper admits this, but the plural 'phase transitions' in the abstract is misleading.\n\nNone of this kills the paper. The qualitative trends might well be right, and the workflow is worth publishing so others can build on it. But the evidence as presented does not establish that the high-temperature behavior is predictive rather than coincidental. The authors should provide at least a few DFT checks at 800–1500 K, add error bars, and clarify exactly which dielectric quantity is being compared to experiment.\n\nThis deserves a serious referee, not a desk reject. I would send it out, with the expectation that the extrapolation issue and the dielectric definition get addressed before publication.","headline":"A useful small-data demonstration of GNN potentials for finite-T dielectric/phase behavior, but the headline results rest on unvalidated temperature extrapolation and an ionic-only dielectric model.","tokens_in":6674,"tokens_out":2026,"would_cite":true,"duration_ms":23987,"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":"A neural network trained on 178 snapshots captures CaTiO3's dielectric curve and phase change.","keywords":["equivariant graph neural networks","machine-learned force fields","Born effective charges","dielectric tensor","finite-temperature molecular dynamics","calcium titanate","perovskite oxides","structural phase transitions"],"falsifier":"At a representative temperature (for example, 1000 K), recompute the dielectric tensor from the same MD trajectories but with Born charges evaluated by density-functional perturbation theory for the sampled snapshots instead of by the equivariant model; if the decreasing-with-temperature trend disappears or reverses, the neural-network charge model, rather than the dynamics, is what reproduces the experimental curve.","tokens_in":5624,"feed_emoji":"⚛️","tokens_out":15609,"duration_ms":139939,"temperature":0.7,"pith_summary":"The authors aim to show that finite-temperature molecular dynamics with machine-learned potentials can capture two hard-to-compute properties of a real perovskite oxide: the temperature dependence of the dielectric tensor and the structural phase transition. They train equivariant graph neural networks on only about 178 density-functional-theory snapshots of CaTiO3, with one network supplying zero-field energies and forces and another supplying Born effective charges. Applying a small electric field and reading the polarization from the field-induced ionic displacements through the Born charges, they reproduce the experimentally observed decreasing, bell-shaped dielectric response of CaTiO3 with temperature and an orthorhombic-to-cubic transition near 1280 K on heating. Because such finite-temperature response properties are normally too expensive for direct ab initio simulation, a working small-data recipe would make routine studies of dielectric and phase behavior in perovskite oxides practical.","feed_headline":"Small neural network reproduces CaTiO3's dielectric trend","feed_subtitle":"Trained on 178 atomic snapshots, the potential reproduces both the dielectric drop and the phase switch of CaTiO3.","key_machinery":"The load-bearing mechanism is a combined machine-learned force field assembled from Eq. 1 and Eq. 2: zero-field forces from MACE plus an electric-field term $|e|\\sum_\\beta \\mathcal{E}_\\beta Z^*_{\\kappa,\\beta\\alpha}$ from Born effective charges predicted by Equivar. Born effective charges are atomic tensors that measure how the force on an atom responds to an applied electric field (equivalently, how polarization changes per atomic displacement), and in perovskite titanates they are anomalously large because Ti-O hybridization shifts charge along bonds. The polarization is then read out as $P_\\alpha = (e/\\Omega)\\sum_{\\kappa,\\beta} Z^*_{\\kappa,\\alpha\\beta} u_{\\kappa,\\beta}$, with the static dielectric tensor obtained from the slope of polarization versus field in finite-temperature MD. Equivariant graph neural networks are the enabler because their outputs transform correctly under rotations and translations, so a model trained on roughly 100-500 perturbed snapshots can carry both forces and tensorial Born charges with close-to-ab initio accuracy.","core_discovery":"The central claim is that a specialized equivariant graph neural network potential trained on a small, system-specific dataset can replace direct ab initio molecular dynamics for finite-temperature response properties. For CaTiO3, the authors combine a MACE force field for zero-field forces and energies with an Equivar model for Born effective charges, and compute the static dielectric tensor from the field-induced displacements via $P_\\alpha = (e/\\Omega)\\sum_{\\kappa,\\beta} Z^*_{\\kappa,\\alpha\\beta} u_{\\kappa,\\beta}$ under an applied field of $5 \\times 10^7$ V/m. The resulting tensor decreases with temperature and follows the shape of the experimental curve, although the absolute values are smaller than experiment because of the density-functional training data and the absence of a low-temperature zero-point plateau in classical dynamics. In the same simulations, heating switches CaTiO3 from the orthorhombic to the cubic phase at about 1280 K, qualitatively reproducing the experimental phase transition, albeit without the intermediate tetragonal phase and at a lower temperature. The authors take this as evidence that small-data equivariant potentials can capture the finite-temperature physics of realistic perovskite oxides.","pith_inferences":["Not stated in the paper: the missing clamped-ion electronic contribution could be estimated by adding a Born-charge-independent electronic dielectric constant from DFPT to the MD ionic values, and the remaining gap would test whether Eqs. 1-2 are complete.","Not stated in the paper: retraining on PBE rather than PBEsol+U and re-running the heating protocol would directly test whether the 1280 K transition temperature and the missing tetragonal phase are functional errors or model errors.","Not stated in the paper: substituting strain for the electric field in the same coupling scheme could yield finite-temperature piezoelectric or flexoelectric tensors with no new architecture."],"forward_implications":["Specialized equivariant GNNs trained on about 178 CaTiO3 snapshots predict forces, energies, and Born charges with errors more than an order of magnitude smaller than a generalist model trained on around 10^6 structures.","The dielectric tensor of CaTiO3 computed from finite-temperature MD decreases with increasing temperature and follows the shape of the experimental curve, so such simulations can be used to screen perovskite dielectrics without expensive ab initio MD.","Heating dynamics of CaTiO3 show an orthorhombic-to-cubic phase transition at roughly 1280 K, demonstrating that structural phase stability can be probed with small-data machine-learned potentials.","Because the classical MD lacks zero-point motion and the training functional suppresses the absolute dielectric value, the computed curve is a qualitative rather than quantitative match; improving the underlying DFT functional or adding quantum nuclear effects should be the next correction.","The same training and field-coupling recipe is presented as applicable to other perovskite oxides and to materials whose dielectric response is dominated by anomalous Born charges."],"supporting_citations":[{"why":"Supplies the MACE message-passing architecture that produces the zero-field forces, energies, and stresses used in the runs.","marker":"6"},{"why":"Introduces the Equivar equivariant GNN for predicting Born effective charges and provides the earlier BM1 baseline model used for comparison.","marker":"11"},{"why":"The generalist MACE-MP-0 model trained on roughly 10^6 structures is the large-data baseline that the specialized small-data models outperform.","marker":"12"},{"why":"Provides the large pretraining dataset behind the generalist baseline, establishing the data-size contrast emphasized by the paper.","marker":"13"},{"why":"Gives the microscopic theory of anomalous Born charges in ABO3 perovskites, the key tensorial quantity the models must reproduce.","marker":"16"},{"why":"Supplies reference Born effective charge values for CaTiO3 that validate the DFPT training data for the equivariant charge model.","marker":"22"},{"why":"Computes the static dielectric constant of CaTiO3 from first principles, providing the zero-temperature reference for the finite-temperature dielectric tensor.","marker":"23"},{"why":"Reports the experimental orthorhombic-to-tetragonal and tetragonal-to-cubic transition temperatures of CaTiO3 used to benchmark the heating simulations.","marker":"24"},{"why":"Measures the experimental temperature-dependent dielectric response of CaTiO3, the curve the MD simulations reproduce qualitatively.","marker":"26"}],"fun_headline_variants":["CaTiO3 dielectric trend captured by tiny ML potential","Small equivariant net matches CaTiO3 dielectric temp curve","Tiny training set, accurate perovskite dielectric simulation","Equivariant GNN model reproduces CaTiO3 phase and dielectric","Neural net potential predicts CaTiO3 finite-T dielectric tensor"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the applied electric field acts on the crystal only through the Born effective charges of individual ions, with no separate electronic (clamped-ion) contribution to the polarization, and that this purely ionic, linear, local response can be compared directly with the experimental total dielectric tensor.","fun_headline_variants_meta":{"raw":{"variants":["CaTiO3 dielectric trend captured by tiny ML potential","Small equivariant net matches CaTiO3 dielectric temp curve","Tiny training set, accurate perovskite dielectric simulation","Equivariant GNN model reproduces CaTiO3 phase and dielectric","Neural net potential predicts CaTiO3 finite-T dielectric tensor"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000396,"raw_usage":{"total_tokens":2046,"prompt_tokens":890,"completion_tokens":1156,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":506,"completion_tokens_details":{"reasoning_tokens":1074}},"tokens_in":506,"tokens_out":1156,"duration_ms":8868,"temperature":1.0,"reasoning_tokens":1074,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:17:46.691517+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"At a representative temperature (for example, 1000 K), recompute the dielectric tensor from the same MD trajectories but with Born charges evaluated by density-functional perturbation theory for the sampled snapshots instead of by the equivariant model; if the decreasing-with-temperature trend disappears or reverses, the neural-network charge model, rather than the dynamics, is what reproduces the experimental curve.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the MACE message-passing architecture that produces the zero-field forces, energies, and stresses used in the runs."}],"review_version":1}