{"id":"8693b6bb-dfa2-4aae-a2cc-3fc742df484d","arxiv_id":"2607.27108","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.5,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"Ab initio concentration-wave EPIs plus Monte Carlo recover three composition-dependent regimes of L12 Ni3Al ordering and Cu separation in Cux(Ni3/4Al1/4)1−x, matching experiment qualitatively.","lead":"A first-principles workflow maps how Cu–Ni–Al alloys form strengthening Ni3Al precipitates across composition, recovering three experimental regimes. It links those regimes to electronic p–d bonding and offers a parameter-light route to coherent precipitation design.","discovery_kind":"new_application","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The entire three-regime phase diagram is sampled with a pair-only Bragg–Williams Hamiltonian linearised about the disordered CPA reference state, yet the paper never checks that this Hamiltonian reproduces the energetics of the fully ordered/phase-separated configurations it is used to simulate.","rationale":"The reader identified the neglected vibrational free energy as the weakest assumption. I agree that is the largest quantitative error source, but it is only partially load-bearing for the stated central claim: the claim is explicitly qualitative, the authors flag the omission prominently (Secs. III B, III D, IV), and the cited ~30% reduction in ordering temperature for Ni3Al would shift boundaries in T without obviously reordering the sequence (L1_2 ordering is strongly favoured energetically and would still precede or coincide with Cu separation). The less-discussed and more structural assumption is that a pair Hamiltonian derived as a second derivative about the disordered state faithfully describes the strongly ordered and phase-separated states that constitute the claimed regimes. Everything in Fig. 7 — the regime boundaries, the two-peak structure at intermediate x, the single peak at high x — inherits this assumption. A secondary concern is finite-size/energy-resolution sensitivity of the small, broad low-T specific-heat peak in 2048-atom cells with 1024 energy bins, but the deposited dataset and published BraWl code make that straightforwardly checkable, and WL convergence criteria are stated. Because the paper claims only qualitative agreement, caveats its temperatures, and the method has prior independent validation in related alloys, my concern does not overturn the reader's ACCEPT; the proposed cross-check would either close the gap or, if it failed, convert the result into a conditional one.","tokens_in":35022,"tokens_out":2364,"duration_ms":54256,"concrete_test":"Using the deposited EPIs (Zenodo dataset), evaluate Eq. 6 for random and L1_2 configurations at x = 0 and compare ΔE to the directly calculated KKR-CPA value of 315 meV/atom (Sec. III E); repeat at x = 0.25 against the quoted 183 and 259 meV/atom values. If the pair Hamiltonian deviates by more than ~15–20%, the inferred transition temperatures and possibly the regime boundaries in Fig. 7 are unreliable; agreement within that band validates the truncation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim — three composition-dependent regimes of ordering and separation — is inferred entirely from Wang–Landau Monte Carlo on the Hamiltonian of Eq. (6), whose EPIs are V = −S^(2), i.e. second concentration derivatives of the KKR-CPA internal energy evaluated about the homogeneous disordered solid solution (Secs. II B–C). This is a harmonic expansion around the random alloy: it is reliable for infinitesimal fluctuations, but the configurations that define regimes (ii) and (iii) — fully L1_2-ordered Ni3Al and macroscopically Cu-separated two-phase states — are far from that reference. Higher-order (three-body and beyond) interactions and any concentration-derivative terms beyond second order are dropped, and the model is further truncated to six fcc shells. Within the paper itself there is an independent benchmark available: Sec. III E reports direct KKR-CPA enthalpy differences (315 meV/atom for A1→L1_2 at x = 0; 183 and 259 meV/atom at x = 0.25), but these are never cross-checked against what the pair Hamiltonian predicts for the same configuration pairs. If the pair model systematically mis-weights ordering vs. separation energetics, the relative placement of the two specific-heat peaks — hence the identification of the intermediate-x regime and the boundary near x ≈ 0.1 — could shift, and the claimed qualitative agreement with Semboshi et al. would be less secure than presented. Prior validation of the method in other Al-bearing systems (Refs. 58, 63) mitigates but does not settle this for a system combining strong ordering with weak, separation-driving Cu interactions.","agreement_with_reader":"partial"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The manuscript applies an established ab initio workflow — KKR-CPA electronic structure of the disordered solid solution, a multicomponent concentration-wave (S^(2)) analysis yielding a Landau-type stability criterion, extraction of real-space effective pair interactions (EPIs), and fixed-lattice Wang–Landau/Metropolis Monte Carlo — to the Cu_x(Ni_3/4Al_1/4)_1−x pseudobinary. The authors report three composition-dependent regimes: (i) at low x, Cu dissolves in L1_2 Ni3Al with a single (virtual) Ni–Al ordering transition; (ii) at intermediate x, high-temperature L1_2 ordering with Cu soluble is followed at lower temperature by Cu/Ni3Al phase separation; (iii) at high x, L1_2 Ni3Al precipitates directly from the solid solution with a single transition. The inferred diagram agrees qualitatively with the experimental pseudobinary diagram of Semboshi et al. (J. Alloys Compd. 921, 166124 (2022)), including a Cu solubility limit near x ≈ 0.1. Bloch spectral function calculations attribute the ordering to strong Al(p)–Ni(d) hybridisation and the Cu expulsion to the energetic offset of the Cu and Ni 3d complexes. Transition temperatures are acknowledged to be overestimated because vibrational entropy and melting are neglected.","tokens_in":35383,"tokens_out":3399,"duration_ms":66102,"significance":"If the results hold, the paper delivers a parameter-free, fully ab initio workflow — KKR-CPA electronic structure, multicomponent concentration-wave theory including charge rearrangement, and Wang–Landau Monte Carlo — that reproduces the experimentally observed pseudobinary phase behaviour of a technologically relevant precipitation-strengthened alloy, while also providing an electronic-structure rationalisation (Al(p)–Ni(d) hybridisation driving ordering; Cu–Ni d-band mismatch driving expulsion of Cu). Specific strengths worth naming: the method is free of empirically fitted parameters; the simulations resolve the L1_2 vs D0_22 ambiguity that the mean-field theory cannot; the comparison against the experimental diagram of Semboshi et al. is a genuinely falsifiable external check; and the work is reproducible, with an open dataset (Zenodo) and simulations performed with the published open-source BraWl package. The Cu–Ni–Al pseudobinary is also a demanding test case for any method that must capture ordering and phase separation within a single Hamiltonian, so this is a useful benchmark demonstration beyond the specific system.","major_comments":[{"comment":"The Hamiltonian of Eq. (6), with EPIs defined as V = −S^(2) about the homogeneous disordered CPA reference (Secs. II B–C) and truncated to six fcc shells, is used via Wang–Landau sampling to simulate configurations far from that reference: fully L1_2-ordered Ni3Al and macroscopically Cu-separated two-phase states. Three-body and higher terms are dropped by construction. The manuscript itself contains an available benchmark that is never used: Sec. III E reports direct KKR-CPA enthalpy differences (A1→L1_2 at x=0: 315 meV/atom; at x=0.25: 183 meV/atom for ordering with Cu dissolved, 259 meV/atom for full separation). The authors should evaluate Eq. (6) on these same configuration pairs and report the comparison. If the pair model mis-weights ordering vs. separation energetics by even a modest factor, the relative placement of the two specific-heat peaks — and hence the identification of a","section":"§II C, Eq. (6); §III D, Fig. 7"},{"comment":"All WL simulations use a single cell of 8×8×8 fcc unit cells (2048 atoms) with periodic boundaries. Phase separation on a fixed lattice in a small periodic cell is known to be sensitive to finite-size effects: miscibility-gap temperatures can be underestimated, and small-cell commensurability can artificially sharpen or merge transitions. The central evidence for the intermediate-x regime is the splitting of one specific-heat peak into two, with the lower-T peak described as 'smaller, broader' — precisely the signature most vulnerable to finite-size smearing. No convergence test in system size is reported. A repeat at a second size (e.g. 12^3 or 16^3 cells) at a few representative compositions (x = 0, 0.25, 0.5, 0.8) would establish whether the two-peak structure and the regime boundary near x ≈ 0.1 are robust, and is inexpensive relative to the calculations already performed.","section":"§III D and Appendix A.3"}],"minor_comments":[{"comment":"The manuscript switches from full-potential PBE (lattice parameters, bulk moduli) to LDA/ASA for the concentration wave analysis. The justification (prior insensitivity to lattice parameter) is reasonable, but a one-line quantitative statement of how much the x=0 EPIs or T_ord shift between the two settings would preempt the obvious question, since the PBE under-binding caveat is itself raised for the structural data.","section":"§III A and Appendix A.1a"},{"comment":"The neglect of vibrational free-energy differences is disclosed and a 30% reduction of T_ord for Ni3Al (Ref. 91) is cited, but the discussion frames the effect purely as a uniform shift of temperatures. Since ordered and separated phases differ in stiffness (Fig. 4, bottom), vibrational entropy could also shift the topology, not just the scale. A sentence acknowledging this — and noting which regime boundaries are most/least sensitive — would strengthen Sec. III D.","section":"§III D"},{"comment":"The ASRO axes in Figs. S17–S20 extend to very large negative values (down to −40) as Al becomes dilute; α_pq_n is unbounded below for c_q → 0, so these magnitudes mostly reflect dilution rather than physical ordering strength. A brief note in the SM caption would help readers interpret the high-x panels.","section":"Supplemental Material"},{"comment":"The non-magnetic treatment is justified by a DLM test at x=0 only (the highest Ni content), which is the right worst-case choice, but the text should state explicitly that local-moment collapse was checked at the optimised lattice parameter and whether any composition dependence was sampled.","section":"Appendix A.1b"},{"comment":"The SM URL placeholder 'URL will be inserted by publisher' should be replaced with the Zenodo DOI before publication; also the duplicated C_V panels (linear and zoomed) in each SM figure would benefit from axis labels distinguishing them.","section":"Supplemental Material / Ref. 94"},{"comment":"Ref. 66: 'agitaion' → 'agitation'. Also 'ferronmagnetism' in Ref. 95 is reproduced from the original title, which is fine, but worth a check against the journal's style.","section":"References"},{"comment":"The claim that the bulk-modulus trend supports 'modulus mismatch strengthening' is plausible but soft as stated; consider citing or briefly noting that coherent strengthening here is primarily misfit-strain and order-strengthening driven, with modulus mismatch a secondary contribution.","section":"§III A"}],"recommendation":"major_revision","confidential_remarks":"The citation pattern leans heavily on the group's own prior S^(2) methodology papers (Refs. 45–48, 58, 59, 63, 106), but this is expected given that they developed the formalism and it is not excessive relative to the method's lineage. Methodologically the work is incremental — an application of an established workflow to a new pseudobinary — but the head-to-head comparison with the Semboshi et al. experimental diagram and the demonstration of concurrent ordering/separation capture make it a useful contribution within the journal's scope. The authors should be aware that the requested Hamiltonian benchmark and finite-size check are the kind of additions a skeptical reader of the follow-on screening studies (proposed in Sec. IV) will demand anyway."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful takeaway is simple: they run the group’s ab initio concentration-wave → effective-pair → Wang–Landau pipeline on the Cuₓ(Ni₃/₄Al₁/₄)₁₋ₓ line and get three composition regimes—Cu soluble in L1₂ at low x, ordering then separation at intermediate x, direct precipitation at high x—in qualitative line with Semboshi. That is a real, usable result for people who care about coherent precipitation and precipitation strengthening, not a re-derivation of the method paper.\n\nWhat works. The internal chain is consistent: X-point instabilities and Ni–Al polarised eigenvectors, short-ranged EPIs dominated by Al–Al/Ni–Al, specific-heat peaks plus Warren–Cowley ASRO, and Metropolis snapshots that match the story. The electronic-structure section (Al–Ni p–d hybridisation; Cu as a diluent with lower-lying d states) actually explains why ordering comes first and Cu is expelled later, rather than just fitting a diagram. They flag the fixed-lattice, no-vibration limitation and the overestimated temperatures up front. Data and code are pointed at Zenodo. Citations to their own prior S^(2) work are methodological, not circular; the experimental comparison is external.\n\nSoft spot, in proportion. The stress-test note is fair: the Bragg–Williams Hamiltonian is second-order about the disordered CPA medium, truncated to six shells, and they never close the loop by comparing its energy for fully ordered or phase-separated cells against the direct KKR-CPA enthalpy differences they already quote in Sec. III E (315, 183, 259 meV/atom). If the pair model mis-weights ordering versus separation, the intermediate-x double-peak topology and the ~0.1 boundary could shift. Prior Al-bearing validations help, but a one-paragraph cross-check would have hardened the claim. Vibrational free energy is omitted by design; they do not pretend the temperatures are predictive. Neither issue breaks the qualitative central claim as stated.\n\nWho it is for: computational alloy design and anyone modelling coherent γ′ precipitation. A serious referee should see it. I would engage—read the SI ASRO plots, maybe ask for the pair-vs-CPA energy table in revision—and I would cite it if I were working on multicomponent precipitation or testing similar workflows. Send it to review.","headline":"Clean application of their established S^(2)+EPI+MC pipeline to Cu–Ni–Al that recovers the known three-regime topology with honest caveats; the pair-model energetics check is the one real soft spot.","tokens_in":36153,"tokens_out":607,"would_cite":true,"duration_ms":18034,"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":"First-principles modelling recovers three composition-dependent regimes of L12 Ni3Al precipitation in Cu–Ni–Al, matching experiment.","keywords":["Cu-Ni-Al alloys","L12 Ni3Al precipitation","concentration wave theory","effective pair interactions","Monte Carlo simulation","coherent precipitation","KKR-CPA","precipitation strengthening"],"falsifier":"Recompute the same compositions with vibrational entropy (or off-lattice relaxations) included and check whether the three-regime topology survives and whether the high-Cu precipitation temperature moves downward toward the experimental ~1020 K at x≈0.8.","tokens_in":35947,"feed_emoji":"⚙️","tokens_out":973,"duration_ms":25093,"temperature":0.7,"pith_summary":"This paper sets out to show that a parameter-free computational workflow—ab initio electronic structure, concentration-wave analysis, and lattice Monte Carlo—can capture both chemical ordering and coherent precipitation in the technologically important Cu–Ni–Al system. Along the pseudobinary line Cu_x(Ni3/4Al1/4)1−x it finds three distinct regimes: Cu dissolved in ordered Ni3Al at low copper content, sequential ordering then phase separation at intermediate content, and direct precipitation of L12 Ni3Al from the solid solution at high copper content. Those regimes line up qualitatively with the experimental phase diagram. The same calculations link the behaviour to electronic structure: strong Al–Ni p–d hybridisation drives ordering, while copper’s lower-lying d bands leave it relatively inert and eventually expelled. A sympathetic reader cares because the workflow offers a practical route to studying precipitation strengthening without empirical fitting, while still connecting microstructure to the underlying electronic mechanisms.","feed_headline":"Three regimes of Ni3Al precipitation recovered from first principles","feed_subtitle":"Ab initio pair interactions plus Monte Carlo match the Cu–Ni–Al phase diagram and link it to electronic structure","key_machinery":"Ab initio concentration-wave (S^(2)) analysis of the disordered CPA medium, which yields reciprocal-space second concentration derivatives of the internal energy that are inverted to real-space atom–atom effective pair interactions for fixed-lattice Monte Carlo (Wang–Landau and Metropolis) sampling.","core_discovery":"Across the pseudobinary Cu_x(Ni3/4Al1/4)1−x system the model produces three composition-dependent regimes of phase behaviour in qualitative agreement with experiment: at low Cu, a single Ni–Al L12 ordering transition with Cu soluble in the ordered phase; at intermediate Cu, high-temperature L12 ordering followed by lower-temperature separation of Cu from L12 Ni3Al; and at high Cu, direct precipitation of L12 Ni3Al from the ternary solid solution with no separate secondary transition.","pith_inferences":["Because the method already recovers both ordering and separation from one set of EPIs, it is a natural candidate for mapping how small quaternary additions (e.g., Cr, Fe, or Zn) shift the high-Cu precipitation boundary and precipitate volume fraction.","The large calculated enthalpy gain from Al–Ni ordering versus the smaller gain from Cu expulsion implies that processing windows that freeze in dissolved Cu inside L12 Ni3Al would degrade transport properties—an experimentally testable materials-design rule.","Extending the concentration-wave step to finite magnetic moments or to vibrational free-energy corrections would turn the present qualitative phase-diagram topology into a quantitative design tool for precipitation-strengthened conductors."],"forward_implications":["At copper-rich compositions the model predicts that L12 Ni3Al precipitates directly, leaving a near-pure Cu matrix that should retain good electrical and thermal conductivity while the precipitates harden the alloy.","The same workflow can be used to screen other multicomponent solid solutions for coherent precipitation without empirical interaction parameters.","Smooth composition dependence of the recovered pair interactions suggests that sparse sampling plus interpolation can map larger ternary or higher-order composition spaces at modest cost.","Electronic-structure diagnostics (p–d hybridisation versus inert late-transition-metal d bands) supply a physical criterion for when ordering and phase separation will compete."],"fun_headline_variants":["Three Cu-dependent regimes for L12 Ni3Al precipitation from ab initio model","Ab initio model recovers three regimes of Ni3Al precipitation in Cu-Ni-Al","Concentration waves plus Monte Carlo yield three Ni3Al phase regimes","Cu content sets three paths for L12 Ni3Al ordering and precipitation","Model links electronic structure to three Cu-Ni-Al precipitation regimes"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"The entire phase diagram is built from a fixed-lattice configurational model that ignores vibrational free-energy differences between ordered and disordered phases, so all transition temperatures are expected to be substantially too high.","fun_headline_variants_meta":{"raw":{"variants":["Three Cu-dependent regimes for L12 Ni3Al precipitation from ab initio model","Ab initio model recovers three regimes of Ni3Al precipitation in Cu-Ni-Al","Concentration waves plus Monte Carlo yield three Ni3Al phase regimes","Cu content sets three paths for L12 Ni3Al ordering and precipitation","Model links electronic structure to three Cu-Ni-Al precipitation regimes"]},"model":"grok-4.5","effort":"low","cost_usd":0.006045,"raw_usage":{"total_tokens":1718,"prompt_tokens":957,"num_sources_used":0,"completion_tokens":81,"cost_in_usd_ticks":60448000,"prompt_tokens_details":{"text_tokens":957,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":680,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":957,"tokens_out":81,"duration_ms":11771,"temperature":1.0,"reasoning_tokens":680,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-30T11:38:32.116343+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Recompute the same compositions with vibrational entropy (or off-lattice relaxations) included and check whether the three-regime topology survives and whether the high-Cu precipitation temperature moves downward toward the experimental ~1020 K at x≈0.8.","supporting_citations":[],"review_version":1}