{"id":"8d6aace2-3d34-4a00-9dc3-539001e0c109","arxiv_id":"2607.10775","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Shell-resolved bond statistics and PCA reveal slow L12-like V-sublattice ordering in CoNiV that RDF convergence alone misses.","lead":"Molecular-dynamics Monte Carlo simulations of equiatomic CoNiV show that chemical short-range order keeps evolving after the radial distribution function has stabilized. The analysis supplies a practical diagnostic so alloy simulations used for nuclear materials do not stop too early.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"Incomplete second-shell equilibration undercuts the L12-like mode claim that the paper treats as established.","rationale":"The Reader correctly flags the imported MTP as an unvalidated assumption and assigns CONDITIONAL. That concern is legitimate, yet the manuscript’s own figures already show a more immediate, internal soft spot: the second-shell order that supplies the dominant PCA mode never saturates. Because the paper’s strongest claim is precisely the existence and interpretability of that mode, incomplete equilibration of shell 2 is the single most load-bearing issue for the argument as written. Extending the trajectories and re-deriving the PCs is a direct, paper-internal check that does not require new potentials or experiments. If the late-window PCA remains stable, the Reader’s MTP caveat still keeps the verdict conditional; if it shifts, the mode analysis itself needs revision. Hence I keep CONDITIONAL but re-center the weakest link on equilibration rather than solely on potential fidelity.","tokens_in":7112,"tokens_out":604,"duration_ms":8297,"concrete_test":"Extend every replica by at least another 150–200 k MC/MD steps (or until the lag-JSD of P_pair plateaus and α_VV^(2) stops decreasing). Recompute the six-replica mean WC matrix for shell 2 and re-run the per-replica-centered PCA on the new late window. If |α_VV^(2)| changes by >0.05 or the shell-2 weight / sign pattern of PC1 rearranges, the L12-like mode interpretation is not yet robust.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that shell-resolved bond statistics plus PCA reveal a coherent V-sublattice (L12-like) ordering amplitude on shells 2 and 4 that continues after RDF convergence. The paper itself states that second-shell WC entries (especially α_VV^(2)) keep evolving through the full 170k-step trajectory and that reported second-shell magnitudes are only lower bounds (Results, Slow Chemical Ordering; Fig. 4). PCA is performed on the post-saturation window (steps ≥ 21 000) that was chosen because first-shell α had flattened; the leading PC is then interpreted as the L12-like mode because its loadings couple shells 2 and 4 with the expected V-pair signs. If the second-shell order parameter is still drifting, the loadings, the 29 % variance attributed to PC1, and the claim of a robust short-range L12-like motif are all conditioned on an incompletely equilibrated ensemble. The MTP-accuracy concern raised by the Reader is real but secondary: even a perfect potential would leave the same gap if the MC sampling has not converged the very shell that defines the mode.","agreement_with_reader":"partial"},"referee_report":{"model":"grok-4.5","summary":"The manuscript analyzes chemical short-range order (SRO) in equiatomic fcc CoNiV using hybrid MD/MC snapshots generated with a level-20 moment tensor potential. Shell boundaries are fixed from the total RDF; shell-resolved Warren–Cowley parameters and pair fractions are then tracked. The dominant pairwise signal is first-shell V–V avoidance and second-shell V–V enrichment, interpreted as a short-range L12-like motif, while shell 3 remains weak. Lagged Jensen–Shannon distances show that RDF statistics plateau while the pair-fraction vectors continue to evolve. PCA of per-replica-centered pair fractions in a post-saturation window (steps ≥ 21 000) yields three modes: a V-sublattice ordering amplitude (PC1, shells 2 and 4), a Ni–Co redistribution mode (PC2), and a mixed Co–V mode (PC3). The authors conclude that scalar RDF convergence can miss slow chemical relaxation and that shell-resolved bond statistics plus PCA provide a compact diagnostic for SRO development.","tokens_in":7402,"tokens_out":1166,"duration_ms":37609,"significance":"If the sampling and potential are adequate, the work supplies a practical, unsupervised protocol (lag-JSD on pair fractions plus PCA of shell-resolved bond statistics) that cleanly separates structural from chemical equilibration in multicomponent alloys. The explicit demonstration that first-shell WC parameters flatten while second-shell V–V order continues, and that this slow direction is recovered as the leading PC, is useful for the community that relies on hybrid MD/MC for SRO. Six independent replicas, between-replica error bars, and the frank statement that second-shell magnitudes are lower bounds are methodological strengths. The main limitation is that the L12-like mode claim rests on an incompletely equilibrated second shell and on an MTP whose ranking of pair energies is not validated inside this study.","major_comments":[{"comment":"Results, Slow Chemical Ordering and Fig. 4: the paper itself states that α_VV^(2) continues to decrease through the full 170 000-step trajectory and that reported second-shell magnitudes are only lower bounds. PCA is performed on the post-saturation window chosen because first-shell α had flattened; PC1 (29 % variance) is then interpreted as a coherent V-sublattice / L12-like mode because its loadings couple shells 2 and 4 with the expected V-pair signs. If the order parameter that defines the mode is still drifting, the loadings, the variance attribution, and the claim of a “robust short-range L12-like motif” are conditioned on an incompletely equilibrated ensemble. Either longer sampling until second-shell stationarity, or a quantitative demonstration that the PC1 direction is stable under further evolution, is needed before the mode can be treated as established.","section":"Results, Slow Chemical Ordering; Fig. 4; Collective Modes"},{"comment":"Methods, Hybrid MD/MC Simulations: all energetics and therefore all SRO are generated by a level-20 MTP imported from a concurrent arXiv preprint (ref. 9) without any independent validation against DFT or experiment inside this manuscript. Because the central physical claim is a specific L12-like pair preference (first-shell V–V avoidance, second-shell enrichment), the ranking of relative pair energies by the potential is load-bearing. At minimum, a short comparison of key WC parameters or formation energies against DFT (or against published experimental CSRO motifs for CoNiV) should be supplied so that readers can judge whether the observed modes are physical or potential-specific.","section":"Methods, Hybrid MD/MC Simulations"}],"minor_comments":[{"comment":"Fig. 3 caption and body text refer to “Fig. ??” / “as shown in Fig. ??”; the cross-reference is broken and should be fixed.","section":"Results, Slow Chemical Ordering; Fig. 3"},{"comment":"Equation numbering: the WC definition is labeled (1) in the Introduction and then the JSD is also labeled (2); the text later says “using Eq. (2)” for WC. Renumber consistently.","section":"Introduction; Methods"},{"comment":"Title and abstract use “CoNiV” / “NiCoV” interchangeably; pick one ordering and keep it throughout.","section":"Title; Abstract"},{"comment":"The shell-weight definition used to localize PCs (“sum of the squared PC-vector entries over the six pair features”) is clear but should be written as an explicit formula so that the 0.78 / 0.94 numbers are reproducible.","section":"Results, Collective Modes"},{"comment":"A brief statement of how many MC swaps are accepted per 100 attempts (or the acceptance rate) would help readers judge sampling efficiency.","section":"Methods"}],"recommendation":"major_revision","confidential_remarks":"The incomplete second-shell equilibration is the load-bearing issue; the MTP-validation gap is real but secondary and fixable with a short DFT benchmark. If the authors can either equilibrate further or show that the PC1 direction is robust under continued drift, the paper becomes a solid methods contribution. Scope is appropriate for a computational-materials or alloy-physics venue."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful takeaway is simple: scalar RDF convergence is not enough for chemical equilibration in multi-principal alloys, and shell-resolved bond probabilities plus a quick PCA give a compact way to track the slow part. That is new relative to the existing CoNiV SRO literature, which already had the first-shell WC signs and L11/L12-like motifs. Here they add the lagged JS comparison (Ppair keeps drifting while RDF is flat) and the three unsupervised modes from per-replica-centered pair fractions.\n\nWhat they do well is the bookkeeping. Six independent replicas, between-replica error bars, an explicit post-saturation window, and the frank statement that second-shell magnitudes are only lower bounds. The loadings make sense: PC1 is mostly shell-2 (and weakly shell-4) V-pair terms with the expected L12-like signs, PC2 is Ni-Co redistribution at fixed V, PC3 mixed Co-V. First-shell flattens fast; second-shell V-V keeps moving. That matches the lag-JSD and the time traces. Math and data handling look clean; citations cover the experimental and prior simulation work without padding.\n\nSoft spots are real but proportionate. The level-20 MTP is taken from a concurrent preprint with no DFT or experimental check inside this study, so the absolute pair preferences rest on that potential. Second-shell order is still evolving at 170k steps, so the 29 % variance on PC1 and the “robust short-range L12-like motif” language are conditioned on an incompletely equilibrated ensemble; they acknowledge the lower-bound issue, which keeps it honest rather than overstated. No code or data release. Free choices (window start after inspecting first-shell WC, fixed shell cuts from total g(r), swap schedule) are standard but not exhaustively tested.\n\nThis is for people who run hybrid MD/MC on concentrated alloys and need a practical stationarity check beyond energy or RDF. It deserves a serious referee; the central observation is directly visible in the figures and does not collapse under the incomplete-equilibration caveat. I would engage with it, cite the diagnostic, and push for longer sampling or potential validation in revision.","headline":"Solid methods note: lag-JSD and PCA cleanly show RDF can miss slow SRO, with a coherent V-sublattice mode in CoNiV, though second-shell sampling is incomplete and the MTP is imported unvalidated.","tokens_in":7987,"tokens_out":566,"would_cite":true,"duration_ms":22373,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Scalar RDF convergence can miss slow chemical ordering in CoNiV; shell-resolved bond statistics and their PCA modes track the real SRO.","keywords":["chemical short-range order","Warren-Cowley parameters","CoNiV","principal component analysis","hybrid MD/MC","Jensen-Shannon divergence","L12-like ordering","machine-learned potential"],"falsifier":"An independent long hybrid MD/MC trajectory of the same alloy driven by a different, validated potential (or by direct DFT-based Monte Carlo) that shows either no second-shell V–V enrichment or first-shell ordering that continues to evolve after the RDF has plateaued would overturn the claimed mode structure and the RDF-versus-SRO timescale separation.","tokens_in":8016,"feed_emoji":"🔬","tokens_out":664,"duration_ms":8321,"temperature":0.7,"pith_summary":"In equiatomic face-centered-cubic CoNiV, atoms are not randomly mixed: vanadium strongly avoids other vanadium atoms as nearest neighbors and piles up at second-neighbor positions, forming a short-range L12-like motif. The paper shows that ordinary radial distribution functions and total energy can look fully relaxed long before this chemical pattern has finished developing. By grouping neighbors into fixed coordination shells, computing the full set of pair probabilities, and running principal component analysis on those vectors, the authors extract three collective modes—the dominant one being coherent V-sublattice ordering on shells 2 and 4. Lagged Jensen–Shannon distances confirm that the bond statistics keep evolving after the RDF has plateaued. The practical payoff is a compact, shell-resolved diagnostic that tells simulators when chemical short-range order has truly converged and which pair correlations move together.","feed_headline":"RDF can look done while alloy chemistry still orders","feed_subtitle":"Shell-resolved bond modes catch the slow V-sublattice ordering that total g(r) misses in CoNiV","key_machinery":"Shell-resolved pair-probability vectors (24-dimensional for three elements and four FCC shells) fed to per-replica-centered PCA, which extracts collective modes whose loadings identify which shells and element pairs vary together; the first mode is the V-sublattice ordering amplitude.","core_discovery":"In hybrid MD/MC snapshots of equiatomic CoNiV at 300 K, first-shell V–V avoidance and second-shell V–V enrichment develop after the radial structure has already stabilized; principal-component analysis of the shell-resolved pair-probability vectors isolates this process as a coherent V-sublattice ordering mode on shells 2 and 4 (with shell 3 remaining weak), together with separate Ni–Co redistribution and Co–V exchange modes, proving that scalar RDF convergence alone is an insufficient equilibration metric for chemical short-range order.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["RDF stabilizes while CoNiV chemical SRO keeps evolving","V-V avoidance lags radial order in equiatomic CoNiV","Shell bond modes track slow SRO after g(r) converges","PCA isolates V-sublattice ordering missed by scalar RDF","Bond statistics relax slower than RDF in NiCoV alloys"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The machine-learned interatomic potential used for every energy evaluation correctly ranks the relative pair preferences that drive the observed ordering.","fun_headline_variants_meta":{"raw":{"variants":["RDF stabilizes while CoNiV chemical SRO keeps evolving","V-V avoidance lags radial order in equiatomic CoNiV","Shell bond modes track slow SRO after g(r) converges","PCA isolates V-sublattice ordering missed by scalar RDF","Bond statistics relax slower than RDF in NiCoV alloys"]},"model":"grok-4.5","effort":"low","cost_usd":0.008024,"raw_usage":{"total_tokens":1910,"prompt_tokens":763,"num_sources_used":0,"completion_tokens":89,"cost_in_usd_ticks":80240000,"prompt_tokens_details":{"text_tokens":763,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":1058,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":763,"tokens_out":89,"duration_ms":12240,"temperature":1.0,"reasoning_tokens":1058,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T09:20:22.340996+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"An independent long hybrid MD/MC trajectory of the same alloy driven by a different, validated potential (or by direct DFT-based Monte Carlo) that shows either no second-shell V–V enrichment or first-shell ordering that continues to evolve after the RDF has plateaued would overturn the claimed mode structure and the RDF-versus-SRO timescale separation.","supporting_citations":[],"review_version":1}