{"id":"e30593b0-9d00-49f4-b710-44aa9891b9ed","arxiv_id":"2607.26597","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Polarization switching in Hf0.5Zr0.5O2 proceeds through cooperative exchange of 3- and 4-coordinated oxygen atoms, with a transient detour and internal volume compensation.","lead":"Using a machine-learned force field, the authors simulate electric-field-driven switching in hafnium-zirconium oxide and find that oxygen atoms swap their 3- and 4-coordinated positions instead of shifting simply. The proposed mechanism could explain why HZO films stay ferroelectric when squeezed into ultra-thin memory devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"MLFF validation does not cover the transient O3c/O4c exchange pathway; the central switching mechanism could be an extrapolation artifact.","rationale":"The reader's weakest-assumption analysis identifies the MLFF's out-of-distribution extrapolation to transient switching configurations as the core risk. My reading of the manuscript supports this: the training set, validation metrics, phonon comparisons, and temperature-dependent trends all sample equilibrium or near-equilibrium phase space, while the central mechanistic claim concerns a short-lived, bond-breaking, strongly field-driven pathway. The ClNEB comparison does not resolve this because it relies on the same model. The proposed DFT single-point check and independent DFT-NEB calculation directly target the gap between what is validated and what is claimed. This does not introduce a new concern beyond the reader's, so the appropriate verdict remains CONDITIONAL: the mechanism is plausible and internally consistent, but the current evidence does not rule out a force-field artifact. I therefore recommend no change to the reader's verdict.","tokens_in":15891,"tokens_out":2997,"duration_ms":32548,"concrete_test":"Extract 50–100 configurations from the EF-MD trajectory, deliberately oversampling the switching windows and the detour region, and recompute their energies, atomic forces, and Born effective charges with VASP PBEsol using the same ENCUT/EDIFF/k-point settings as §2.1. Compare these DFT references with MACEField predictions; if force MAEs on these transient configurations are substantially larger than the reported 12.1 meV/Å validation MAE (e.g., >50 meV/Å), the pathway is not reliably predicted. Additionally, recompute one representative O3c→O4c exchange minimum-energy path with DFT-NEB (or with an independently trained MLFF); if the DFT path does not reproduce the detour, the mechanism should be treated as unconfirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that polarization switching in HZO proceeds via cooperative O3c/O4c exchange with a transient detour, rather than via conventional S:N/S:T displacements. For this to hold, the fine-tuned MACEField model must remain quantitatively accurate on transient, strongly field-driven, bond-breaking configurations visited during switching. The training set (§2.1) contains only 2,257 SQS structures from four equilibrium crystalline phases (Pca21, P21/c, P42/nmc, Pbca), sampled via heating MD. Validation (§3.1, Figures S1–S3) is performed on a held-out test split from the same phases and on phonon dispersions; none of these tests include configurations on the O3c/O4c exchange pathway, where coordination numbers change and new cation–oxygen bonds form. The ClNEB comparison in §3.3 uses the same force field, so it does not provide independent validation of the pathway. If the MLFF is biased in this out-of-distribution region, the detour trajectory, the O3c/O4c classification, and the internal self-compensation volumes could all be simulation artifacts, even if the P–E loop and equilibrium properties are reproduced well. The specificity of the claimed detour—present in O3c→O4c but absent in O4c→O3c—makes the conclusion especially sensitive to small force errors near the newly forming fourth bond.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents electric-field-induced molecular dynamics simulations of Hf0.5Zr0.5O2 using a fine-tuned MACEField machine-learning force field that predicts energies, forces, and Born effective charges. The authors reproduce a P–E hysteresis loop and report that polarization switching proceeds via cooperative O3c↔O4c oxygen exchange along a 'C:N34ex' pathway, with an asymmetric 'detour' in the O3c→O4c direction that they attribute to transient formation of a fourth cation–oxygen bond. They further propose an 'internal self-compensation' mechanism in which local tetrahedral volume expansion and contraction cancel, explaining the small macroscopic strain change during switching and HZO's ultra-thin-film stability. The central claim is that this dynamic exchange, rather than conventional S:N/S:T displacement, is the intrinsic switching mechanism.","tokens_in":16227,"tokens_out":7746,"duration_ms":77358,"significance":"If substantiated, the work would provide a finite-temperature, dynamic atomistic mechanism for polarization switching in hafnia-based ferroelectrics, going beyond static NEB analyses and offering a microscopic rationalization of HZO's strain tolerance. The paper has clear strengths: a DFT-validated MLFF with BEC capability, explicit P–E hysteresis simulation, reproduction of experimental trends (temperature-dependent coercive field and permittivity, d111 spacing), supercell-size convergence checks, and a concrete, falsifiable pathway. However, the central mechanism rests on MLFF extrapolation to transient, bond-breaking configurations and on a single MD trajectory; these issues must be addressed before the claims are fully supported.","major_comments":[{"comment":"The MLFF is trained on 2,257 structures from four equilibrium phases (Pca21, P21/c, P42/nmc, Pbca) and validated on a held-out test split from the same phases plus phonon dispersions. The switching pathway, however, involves transient O3c/O4c exchange with bond breaking and formation under a strong electric field. No validation structures are drawn from this region, and the ClNEB comparison in §3.3 cannot supply independent confirmation because it uses the same MLFF. The paper should provide direct DFT (VASP) energies, forces, and BECs for representative frames along the MD switching path (e.g., 20–50 configurations before/after switching), or a DFT-NEB path for at least one O3c→O4c event. Without this, the detour and exchange mechanism may be an artifact of force-field extrapolation.","section":"§2.1, §3.1"},{"comment":"The central observation is based on a single MD trajectory: one SQS arrangement, one set of initial velocities, and one field-sweep rate of 5.0 × 10⁻⁵ MV/(cm·fs). This sweep rate is extremely fast compared with experimental switching conditions, and the pathway may be rate-dependent. No replicate SQS cells, repeated runs with different random seeds, or error bars are reported. The statement in §3.2 that the results 'definitively establish' the mechanism is therefore too strong. At a minimum, repeat the switching simulation with 2–3 independent SQS cells (or different initial velocities) and report whether the O3c/O4c exchange, detour asymmetry, and self-compensation are reproduced; ideally, test at least one slower sweep rate.","section":"§2.4, §3.2"},{"comment":"The random baselines P<100> = 0.7008, P<110> = 0.8354, and P<111> = 0.7577 are inconsistent with the stated definition of P as <cos²α> for a fixed reference axis: for uniformly random 3D vectors, this expectation is 1/3 for every axis. The reported values appear instead to correspond to a metric based on the nearest angle among all equivalent directions in the <uvw> family. If that is the intended definition, Eq. (1) and the surrounding text must be revised to state this explicitly and to derive the baselines. Since the orientation sequence <111> → <100> → random → <111> in Figure 2 is the main evidence for the exchange mechanism, this ambiguity is load-bearing and must be resolved.","section":"§2.3, Eq. (1)"},{"comment":"The 'internal self-compensation mechanism' is inferred from the qualitative observation of simultaneous expansion and contraction of tetrahedral volumes. No quantitative correlation, net volume balance, or error analysis is given; only 20 structures are plotted. The claim that local volume changes are 'perfectly offset' is not demonstrated. Additionally, the text overstates the macroscopic effect: §3.4 says 'without altering the macroscopic lattice constant at all', but earlier in the same section volume changes of ~1.7% and a-axis changes of ~1% are reported. Please provide statistical distributions/correlations and reconcile the wording.","section":"§3.4, Figure 7"}],"minor_comments":[{"comment":"The O3c/O4c coordination-number assignment is never explicitly defined. Specify the cutoff distance or neighbor criterion used; the entire classification depends on it.","section":"§2.1/§3.2"},{"comment":"Eq. (1) appears with a typographical issue ('cos2α' rather than a superscripted square). Also clarify whether α is measured to a single fixed axis or to the closest direction in the family.","section":"§2.3"},{"comment":"The phrase 'a microscopic physical origin of for HZO's exceptional ability' contains a grammatical error and should be corrected.","section":"Abstract"},{"comment":"The trained MLFF models and trajectories are only available 'upon reasonable request.' For reproducibility, please deposit the model checkpoints, training data, and representative trajectories in a public repository.","section":"Data Availability"},{"comment":"The distinction from Ref. [16] (Mu et al.) should be sharpened: the paper initially describes that work as proposing a 'theoretical possibility', but then the C:N34ex pathway of [16] is directly invoked. State explicitly what is new beyond Ref. [16] (e.g., the dynamic observation, the detour, and the self-compensation) to avoid ambiguity about novelty.","section":"Introduction / §3.3"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely and important question and the methodology is a natural fit for the journal. The main risk is that the central mechanism is drawn from a single MLFF trajectory whose force field is not validated on the transient exchange pathway. If the authors can supply DFT validation on switching-path configurations and demonstrate reproducibility across SQS arrangements, the work could become a strong contribution. The heavy self-citation pattern is understandable in context and does not appear inappropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: this is a useful dynamic complement to Mu et al.'s static NEB picture of HZO switching, and the detour plus internal volume compensation are genuinely new observations. But the central mechanism is supported by one MLFF trajectory, and the MLFF has not been validated on the very configurations that define the mechanism. The claim as stated is overstrong.\n\nWhat the paper does well: it fine-tunes MACEField on a decent DFPT dataset (2,257 structures, four phases, energies/forces/BECs), and the reported test errors and phonon agreement are good. The P-E loop for the defect-free ideal cell is in a reasonable place (Ec 7.59 MV/cm, Pr 48 μC/cm²), and the authors explicitly argue why that differs from experiment. The O3c→O4c detour—present in one direction and absent in the reverse—is a specific, falsifiable shape of claim that can't be dismissed as a generic displacement artifact. The clamped-cell result and the local tetrahedral volume compensation give a coherent explanation for clamping resistance.\n\nThe soft spots are real. The training set samples four equilibrium phases; the switching path is out of that distribution. The ClNEB 'confirmation' uses the same force field, so it is an internal consistency check, not independent validation. There is one SQS arrangement, one trajectory, no error bars, no replicates; the coordination-number cutoff is not stated; data and code are only available 'on request.' The field sweep rate is far from experimental, and the ideal lattice gives intrinsic behavior, so 'endurance and scalability' language goes beyond the evidence.\n\nIs this a deal-breaker? Not obviously. The mechanism is not fitted to the target result; it emerges from the MD. The asymmetry of the detour and the averaged-path reproducibility make a pure noise artifact less likely. But an out-of-distribution bias in the MLFF could produce exactly this kind of specific artifact, so the central claim should be treated as a strong hypothesis from a single simulation, not a definitive mechanism.\n\nRecommendation: send to peer review. A good referee should push for uncertainty quantification, a second SQS arrangement, and ideally a few DFT or independent MLFF checks on the exchange path—or at least release the models and trajectories. The paper deserves referee time, but it needs revision before it supports the strength of its conclusions.","headline":"Useful dynamic complement to the static NEB picture, but the central mechanism rests on one MLFF trajectory and needs validation before it can be sold as definitive.","tokens_in":16735,"tokens_out":3522,"would_cite":true,"duration_ms":39116,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["77.80.-e","77.80.Fm","71.15.Pd"],"model":"deepseek-v4-flash","headline":"Polarization switching in ferroelectric HZO proceeds not by simple oxygen shifts but by a cooperative exchange between 3- and 4-coordinated oxygen, with a transient detour that keeps the lattice nearly strain-free; the paper establishes thi","keywords":["HZO","ferroelectric switching","O3c/O4c exchange","machine learning force field","electric-field-induced MD","polarization hysteresis","strain-free switching","ultra-thin film stability"],"falsifier":"A direct density-functional-theory calculation (static or AIMD) of the O3c-to-O4c transition under an applied electric field that shows a straight minimum-energy path without the detour, or that shows the machine-learned force field's forces deviate substantially (e.g., >100 meV/Å) from DFT forces on the detour configurations, would falsify the proposed dynamic mechanism.","tokens_in":15759,"feed_emoji":"🔁","tokens_out":2586,"duration_ms":31267,"temperature":0.7,"pith_summary":"The paper claims that polarization switching in ferroelectric HZO is driven by the dynamic mutual exchange of 3-coordinated (O3c) and 4-coordinated (O4c) oxygen atoms, rather than by the conventional simple-displacement S:N/S:T models. Using electric-field-induced molecular dynamics with a machine-learned force field fine-tuned on first-principles data, the authors reproduce the P-E hysteresis loop and directly observe the oxygen exchange pathway. During the O3c-to-O4c transition, oxygen atoms follow a curved 'detour' trajectory caused by transient cation-oxygen bond formation, while the reverse transition is more direct. Crucially, the paper identifies an 'internal self-compensation mechanism' in which local volumetric expansion and contraction from the two simultaneous transitions offset each other within the cell, keeping the macroscopic lattice nearly unchanged. This mechanism provides a microscopic origin for HZO's exceptional strain-free switching and its ability to retain ferroelectricity in ultra-thin films, and suggests that preserving cooperative O3c/O4c exchange pathways, rather than minimizing individual atomic displacements, is the key design principle for endurance and scalability.","feed_headline":"Oxygen swapping, not shifting, drives HZO polarization switching","feed_subtitle":"Field-driven MD shows cooperative O3c/O4c exchange that keeps the lattice nearly strain-free, explaining ultra-thin-film stability.","key_machinery":"The central machinery is a fine-tuned machine-learned force field that predicts energies, forces, and Born effective charges on-the-fly, combined with electric-field-induced molecular dynamics under an applied field (swept 0 -> +10 -> -10 MV/cm) in an NPT ensemble. The field-induced force on each atom uses the predicted Born effective charge tensor. To quantify oxygen motion, the paper introduces the 'degree of orientation' P<uvw>, which measures how aligned oxygen displacement vectors are with crystallographic directions, with known random baselines. The trajectory analysis uses the cation cage (the rectangular cuboid of 4-coordinated cations) and the distance from its center of mass to ide","core_discovery":"The central claim is that polarization switching in Hf0.5Zr0.5O2 proceeds through the C:N34ex pathway: oxygen atoms remain inside the 4-coordinated cation cage without crossing the cation plane, while exchanging their coordination number between 3 and 4. The forward O3c-to-O4c transition involves a distinct detour in the x-direction, traced to the oxygen being pulled by a newly forming cation-oxygen bond. The forward and reverse transitions are not time-reversals: the initial O3c state sits close to the cation center, while the final O3c state rests farther away, explaining the asymmetry. Throughout the cycle, local tetrahedral volume expansions from O4c-to-O3c transitions are offset by cont","pith_inferences":["A testable extension: if O3c/O4c exchange is the rate-limiting step, then isotopic tracer or time-resolved X-ray absorption experiments that track oxygen coordination changes during field cycling should observe the detour signature and the asymmetry between forward and reverse switching.","The authors' mechanism suggests that oxygen vacancy formation—which would remove O3c or O4c sites—could disrupt the cooperative exchange and be a primary cause of wake-up or fatigue; this is an editorial inference beyond the paper's defect-free simulations, but it offers a concrete avenue for future study.","The same machine-learning MD approach could probe whether other fluorite ferroelectrics (e.g., doped HfO2 or ZrO2) share the O3c/O4c exchange mechanism, or whether their switching follows different transient pathways; the paper does not make this comparison, but it is a natural next step.","The detour trajectory's dependence on field sweep rate is not discussed in the paper; if the exchange mechanism is cooperative, one would predict that faster sweeps alter the detour shape or even suppress the exchange in favor of a more direct path, a prediction that could be checked with the same simulation framework."],"forward_implications":["If the O3c/O4c exchange is the true switching mechanism, then static minimum-energy-path analyses alone are insufficient: cooperative, transient bond formation under an applied field must be modeled to capture ferroelectric dynamics.","The internal self-compensation mechanism implies that HZO's switching is intrinsically nearly strain-free, explaining why ferroelectricity survives under strong epitaxial clamping and in ultra-thin films without external strain relief.","The design principle for endurance and scalability shifts from minimizing individual atomic displacements to preserving the integrity of the cooperative O3c/O4c exchange pathways, e.g., by avoiding defects that disrupt coordination changes.","The calculated coercive field (7.59 MV/cm) and remnant polarization (48 μC/cm2) in a defect-free ideal lattice exceed typical experimental values, supporting the view that defects such as oxygen vacancies are what reduce coercivity in real devices.","Because the forward and reverse transitions are geometrically asymmetric (detour vs. direct), the switching dynamics may be intrinsically rate-dependent and direction-dependent, which could matter for device speed and reliability modeling."],"fun_headline_variants":["HZO polarization switch runs on cooperative oxygen exchange","Oxygen swap, not displacement, flips HZO polarization","Transient detour in oxygen exchange mediates HZO switching","HZO's strain-free switching comes from O3c/O4c exchange","Cooperative oxygen detour keeps HZO lattice stable during switching"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The machine-learned force field, trained on 2,257 structures sampled from four equilibrium crystalline phases, extrapolates accurately to the transient, bond-breaking, strongly field-driven configurations visited along the O3c/O4c exchange pathway, which were not explicitly included in the training set.","fun_headline_variants_meta":{"raw":{"variants":["HZO polarization switch runs on cooperative oxygen exchange","Oxygen swap, not displacement, flips HZO polarization","Transient detour in oxygen exchange mediates HZO switching","HZO's strain-free switching comes from O3c/O4c exchange","Cooperative oxygen detour keeps HZO lattice stable during switching"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000635,"raw_usage":{"total_tokens":2826,"prompt_tokens":865,"completion_tokens":1961,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":609,"completion_tokens_details":{"reasoning_tokens":1876}},"tokens_in":609,"tokens_out":1961,"duration_ms":13805,"temperature":1.0,"reasoning_tokens":1876,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T12:35:27.925288+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct density-functional-theory calculation (static or AIMD) of the O3c-to-O4c transition under an applied electric field that shows a straight minimum-energy path without the detour, or that shows the machine-learned force field's forces deviate substantially (e.g., >100 meV/Å) from DFT forces on the detour configurations, would falsify the proposed dynamic mechanism.","supporting_citations":[],"review_version":1}