{"id":"942833d7-ac9d-47d2-825b-8b949b3e59da","arxiv_id":"1908.01880","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The 3D Ising model with bond disorder generated from critical pure Ising configurations shows a new universality class with νd = 1.13(5) and ηd = 0.48(3), distinct from the Weinrib-Halperin prediction.","lead":"This paper uses Monte Carlo simulations to test how slowly decaying power-law-correlated disorder changes the critical behavior of the three-dimensional Ising model. It reports new critical exponents that disagree with the standard Weinrib-Halperin prediction and shows that a finite-sample two-peak anomaly is a finite-size effect.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported exponents rest on a beta_c estimate validated with the same nu_d they are used to determine; a single-peak M=0 cross-check is needed before the new universality class is accepted.","rationale":"The reader's weakest assumption concerned the inheritance of the disorder correlation exponent a=d-2+eta_pure and the non-Gaussian nature of the bimodal bond distribution. Those are real caveats, but the spin-channel contribution to the bond-bond correlator does give the leading r^{-(d-2+eta)} decay, so the correlation exponent is likely correct to leading order; the non-Gaussianity weakens the comparison to WH but not the existence of a new universality class if the FSS is sound. The more load-bearing uncertainty is the numerical extraction of beta_c and the exponents in the presence of the two-peak structure, together with the circular validation of the beta_c extrapolation. The paper already possesses the M=0 data needed to settle this, but does not report exponents for it. The overall result is plausible and the simulations are extensive, so a conditional verdict remains appropriate; the proposed cross-check would either affirm or disprove the central exponent claim.","tokens_in":11531,"tokens_out":10710,"duration_ms":157704,"concrete_test":"Repeat the finite-size scaling analysis on the restricted M=0 disorder ensemble, which shows a single susceptibility peak for all sizes in Fig. 2. Extract beta_c, nu_d, and eta_d from g_T and chi using the same seven largest sizes, or preferably perform a global simultaneous fit of beta_c, nu, eta, and a correction-to-scaling term. If the M=0 exponents agree with nu_d=1.13(5) and eta_d=0.48(3) within error bars, the two-peak ambiguity and the circular beta_c validation are not biasing the result. If they disagree, the new universality class claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of a new universality class with nu_d=1.13(5), eta_d=0.48(3) depends on the finite-size scaling analysis at the thermodynamic critical temperature. The paper estimates beta_c by extrapolating the left peak of g_T and chi with a cubic polynomial in 1/L (Sec. III B), and then validates this extrapolation using scaling fits that assume nu_d=1.13, the very exponent obtained from the subsequent FSS. This is circular. The two-peak structure for the unrestricted disorder ensemble means the left peak is one of two competing features; for smaller sizes the right peak is comparable, and the identification of the true thermodynamic transition is not independently established. The paper shows in Fig. 2 that the restricted M=0 disorder ensemble has a single peak for all sizes and argues it has the same thermodynamic limit, but no exponent analysis is reported for that ensemble. Consequently, the reported exponents could be effective values reflecting the beta_c estimate and the disorder-averaging procedure rather than an asymptotic universality class. In addition, the Weinrib-Halperin comparison is explicitly outside its Gaussian disorder assumption, so the observation of a new fixed point does not by itself validate the WH-predicted fixed point or the quantitative exponent comparison.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the three-dimensional Ising model with quenched random bonds generated from critical equilibrium configurations of the pure 3D Ising model. The random couplings take two values and are intended to be power-law correlated with exponent a≈1.036. Using parallel tempering and population annealing, the authors perform large-scale Monte Carlo simulations for system sizes up to L=50 and analyze finite-size scaling of the Binder-ratio derivative and the susceptibility. They report a single phase transition in the thermodynamic limit with exponents νd=1.13(5) and ηd=0.48(3), which differ from both the pure Ising and uncorrelated-disorder universality classes and from the Weinrib-Halperin prediction. They also argue that the double-peak structure in disorder-averaged quantities for finite sizes is a finite-size effect, and they contrast this with a bilayer model that genuinely has two transitions.","tokens_in":11744,"tokens_out":12860,"duration_ms":140983,"significance":"If confirmed, the result is significant: it provides a simple numerical construction of power-law-correlated disorder with a small decay exponent, evidence for a universality class different from both pure and uncorrelated-disorder Ising models, and a concrete test of the Weinrib-Halperin criterion outside its Gaussian-disorder assumption. The strengths of the paper include the use of two independent Monte Carlo algorithms, large disorder samples (2000-5000 realizations), bootstrap error estimates, the data collapse in Fig. 4, and the direct comparison with pure, uncorrelated, and layer models. The main weaknesses are the circular validation of the βc extrapolation, the absence of an exponent analysis for the single-peaked M=0 ensemble, and the unproven inheritance of the spin correlation exponent by the composite bond disorder.","major_comments":[{"comment":"The thermodynamic transition temperature βc=0.1396(3) is obtained by a cubic-polynomial extrapolation of the left-peak positions in 1/L, yet the only stated validation of that extrapolation is that it was verified using scaling fits assuming νd=1.13 estimated below—the very exponent that is then extracted from finite-size scaling at the estimated βc. This validation is circular. Please replace it with an independent determination of βc, for example a simultaneous fit of βc(L)=βc + a L^{-1/νd} with νd free, or an analysis of the M=0 ensemble, whose susceptibility is single-peaked at all sizes (Fig. 2).","section":"Sec. III B (paragraph after Fig. 3)"},{"comment":"The restricted M=0 disorder ensemble is used only to demonstrate that the two-peak structure is a finite-size effect; no finite-size scaling or exponent estimates are reported for that ensemble. Because the reported νd and ηd come entirely from the unrestricted M≠0 ensemble, where the disorder-averaged χ and c have two competing peaks, a quantitative cross-check of the exponents in the M=0 ensemble is needed to rule out that νd=1.13(5) and ηd=0.48(3) are effective values tied to the peak-selection procedure.","section":"Sec. III A and Sec. III B, Figs. 2-4"},{"comment":"The statement that the random couplings inherit the spin correlation exponent a=d−2+η_pure≈1.036 is asserted without derivation. Since Jij=1+(Si+1)(Sj+1)/4 is a composite of neighboring spins, its covariance contains both spin-spin and energy-density contributions; the leading large-distance decay should be derived explicitly, and ideally confirmed by a direct measurement of the bond-bond correlation function, because the comparison with the Weinrib-Halperin criterion depends directly on the value of a.","section":"Sec. II A, second paragraph"},{"comment":"The exponents are obtained from power-law fits over the seven largest sizes, but the paper does not report the fit range, goodness-of-fit, or the stability of νd and ηd when the analysis window is varied. Since Fig. 4(a) shows visible curvature in χ for smaller sizes, evidence that the corrections to scaling have converged by L≈24 is necessary to support the claim that νd=1.13(5) and ηd=0.48(3) represent an asymptotic universality class rather than an effective finite-size description.","section":"Sec. III B, Fig. 4"}],"minor_comments":[{"comment":"The phrase 'Ideally one may would like' should read 'Ideally one would like'.","section":"Sec. III B"},{"comment":"The entry for the PT simulations is unclear: '16 − 6×10^6' should specify the number of replicas and the number of sweeps unambiguously.","section":"Table I"},{"comment":"The expression m0(βJdn/Jup) should be written as m0(β Jdn/Jup) or with an explicit multiplication symbol to avoid ambiguity.","section":"Eq. (11)"},{"comment":"The sentence 'This prediction agrees with our simulation results but the predicted exponents disagree' is potentially confusing; since the model explicitly violates the Gaussian-disorder assumption of WH, please clarify whether the agreement is only about the existence of a new fixed point, with the exponent comparison intended as a heuristic test.","section":"Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is in scope for the journal and the simulations are substantial. My main request is an independent beta_c and finite-size scaling analysis for the M=0 ensemble, or an equivalent non-circular determination of beta_c; without it, the central universality-class claim rests on a single, partly circular analysis path. The layer-model section is somewhat tangential but acceptable as presented. I do not see grounds for rejection if the requested cross-check is provided."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know before you read it. The disorder-generation scheme—equilibrium pure Ising configurations mapped onto two-valued bonds—is new for 3D random-bond Ising models, and the finite-size double-peak effect in disorder-averaged quantities is a real, interesting observation. The second thing is that the headline exponent claim (nu_d=1.13(5), eta_d=0.48(3)) is plausible but not yet nailed down: the beta_c used in the FSS fits is extracted by a cubic in 1/L extrapolation of the left peak, and the validation of that extrapolation uses the fitted nu_d from the same beta_c. That is circular, though mildly so.\n\nWhat the paper does well: large-scale simulations with PA and PT that agree, bootstrap error bars, thousands of disorder realizations, and a data collapse that works for the seven largest sizes. The double-peak structure being a finite-size effect is supported by the M=0 restricted ensemble having a single peak and by the two peaks converging with L. The layer model is a nice, simple contrast.\n\nSoft spots, in order of importance. First, the M=0 single-peak ensemble is displayed but never given its own FSS analysis. That would be an independent way to estimate beta_c and exponents, avoiding the circular validation. Second, the assumption that the bond disorder inherits the pure spin correlation exponent a = d - 2 + eta_pure is physically plausible—the bond operator is dominated by S_i S_j at criticality—but it is asserted rather than demonstrated. A direct measurement of the bond-bond correlation function would settle it. Third, the WH comparison is strictly outside its validity because the disorder is bimodal, not Gaussian; the paper acknowledges this, but it means the disagreement with WH is not a sharp test of the theory. Also, the L=4 data are dropped without much discussion; not a big deal.\n\nNone of these sink the paper. The central scenario—one transition in the thermodynamic limit with exponents intermediate between pure and WH—is credible, and the method opens a useful door for generating correlated disorder. The main fix is to repeat the FSS analysis on the M=0 ensemble. I'd send it to a good referee, with the request that the M=0 analysis be added before acceptance. The audience is people working on correlated disorder, random-bond Ising models, and FSS methodology.","headline":"Careful Monte Carlo study with a new disorder-generation scheme, but the new-universality-class claim rests on a beta_c extrapolation validated with the very exponent it produces.","tokens_in":12310,"tokens_out":2964,"would_cite":true,"duration_ms":31868,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B27","82B20","82B80"],"pacs":["05.50.+q","64.60.Fr","75.10.Hk","75.40.Mg"],"model":"deepseek-v4-flash","headline":"The 3D Ising model with power-law-correlated bond disorder from critical Ising configurations has a single transition in a new universality class, νd=1.13(5), ηd=0.48(3), clearly different from the Weinrib-Halperin predictions.","keywords":["3D Ising model","power-law-correlated disorder","quenched disorder","Weinrib-Halperin criterion","universality class","finite-size scaling","Monte Carlo simulation","population annealing"],"falsifier":"Compute the disorder-averaged bond correlation function C(r)=[(J_ij−J̄)(J_kl−J̄)] as a function of separation in large systems and extract its decay exponent; if it is much larger than 1.036, or if the two finite-size peaks move apart rather than merging when L exceeds 50, the claim of a single new universality class driven by slow power-law disorder is falsified.","tokens_in":11316,"feed_emoji":"🧲","tokens_out":6833,"duration_ms":60317,"temperature":0.7,"pith_summary":"This paper reports large-scale Monte Carlo simulations of the three-dimensional Ising model with quenched random bonds that are power-law correlated, generated by freezing critical equilibrium configurations of the pure Ising model. The authors aim to test the Weinrib-Halperin criterion, which predicts that sufficiently slow power-law disorder creates a new universality class with specific exponents. They find that a new universality class does emerge and that there is one phase transition in the thermodynamic limit, with exponents νd=1.13(5) and ηd=0.48(3). These values differ strongly from the Weinrib-Halperin predictions (νd≈1.93, ηd=0), so the theory's quantitative predictions need correction for this bimodal disorder distribution. Double peaks in disorder-averaged susceptibility and heat capacity at finite sizes are shown to be finite-size effects that merge in the thermodynamic limit.","feed_headline":"New universality class for 3D Ising with correlated disorder","feed_subtitle":"Monte Carlo finds νd=1.13(5), ηd=0.48(3), far from the predicted νd≈1.93.","key_machinery":"The central object is the disorder-generation map Jij=1+(Si+1)(Sj+1)/4, applied to equilibrium critical configurations of the pure 3D Ising model. It defines two-valued random bonds (1 and 2) whose spatial correlations inherit the pure spin correlation exponent a=d−2+η_pure≈1.036, yielding slowly decaying power-law disorder with a<d. This construction is what brings the Weinrib-Halperin criterion into play, since that criterion predicts a new fixed point for such slow decay, with νd=2/a≈1.93 and ηd=0. The analysis machinery is finite-size scaling of disorder-averaged susceptibility, χ∼$L^{{2−η}}$, and Binder-parameter derivative, gT∼$L^{{1/ν}}$, at the extrapolated transition temperature, with parallel tempering and population annealing Monte Carlo used for equilibration.","core_discovery":"The central discovery is that the 3D Ising model with quenched bond disorder constructed from critical pure Ising spin configurations—setting Jij=1+(Si+1)(Sj+1)/4 so that bonds inside spin-up clusters are strong (J=2) and all others weak (J=1)—belongs to a universality class distinct from both the pure Ising class and the Weinrib-Halperin prediction. Finite-size scaling of the Binder-parameter derivative and susceptibility at the extrapolated critical temperature gives νd=1.13(5) and ηd=0.48(3). The authors interpret the two-peak structure seen in disorder-averaged quantities for small systems as a finite-size artifact of the disorder-averaging and generation procedure, not as two transitions: each disorder realization has a single peak, the two peak temperatures approach each other as L grows, and restricting the generating configurations to zero magnetization eliminates the second peak.","pith_inferences":["If the reported exponents hold, the specific-heat exponent given by hyperscaling αd=2−dνd≈−1.39 is negative, so the transition should remain sharp at large sizes, with disorder rounding rather than a true double transition; this is a checkable prediction for future larger-scale runs.","The same disorder-generation recipe applied to other O(N) models in three dimensions would produce power-law-correlated bonds with a≈1+η_pure(N); comparing νd/(2/a) across models could reveal whether the deviation from the Weinrib-Halperin value seen here is a universal correction or specific to the bimodal bond distribution.","The double-peak diagnostic suggests a practical rule for simulations of correlated disorder: restrict the generating ensemble to fixed magnetization; if the second peak disappears, it is a finite-size artifact, while persistence would signal genuine multiple transitions."],"forward_implications":["The long-range correlations in the disorder are relevant: the measured exponents differ from both pure Ising (ν≈0.630, η≈0.036) and uncorrelated bond disorder (ν≈0.685, η≈0.036), confirming a distinct universality class.","There is a single phase transition in the thermodynamic limit despite finite-size two-peak structures; the double peaks arise from disorder realizations with a majority of strong or weak bonds and vanish when the generating ensemble is restricted to zero magnetization.","The Weinrib-Halperin criterion correctly predicts that a new universality class emerges for slowly decaying power-law disorder, but its quantitative exponents (νd=2/a, ηd=0) do not describe the bimodal bond distribution studied here; higher-order corrections or a different fixed point are needed.","A layered non-random model with strong and weak bond halves has two genuine transitions and can be tuned to design magnetization curves, for example m(β)=∑_i ω_i m0(J_i β) for multilayers.","Parallel tempering is more efficient than population annealing for studying phase transitions of pure or weakly disordered systems, where equilibrating directly near the transition is faster than a full temperature sweep."],"supporting_citations":[{"why":"Supplies the Weinrib-Halperin criterion and the predicted fixed-point exponents νd=2/a, ηd=0 that the paper tests.","marker":"[2]"},{"why":"Provides high-precision pure Ising exponents (νpure=0.629971(4), ηpure=0.036298(2)) used to set the disorder correlation exponent a=d−2+ηpure.","marker":"[22]"},{"why":"Provides the high-precision inverse critical temperature βc,pure of the pure 3D Ising model used to generate critical equilibrium configurations.","marker":"[21]"},{"why":"Supports the claim that restricted and unrestricted disorder ensembles have the same thermodynamic limit, used to argue the two peaks are finite-size effects.","marker":"[27]"},{"why":"Supplies the uncorrelated bond-disorder exponents used as a baseline showing correlations change the universality class.","marker":"[28]"},{"why":"Provides a line-defect example where measured exponents sit between pure and Weinrib-Halperin values, a pattern similar to the present result.","marker":"[4]"},{"why":"Reviews the controversy over the Weinrib-Halperin exponents and offers comparison for correlated disorder simulations.","marker":"[5]"},{"why":"Gives a 2D Potts model with a similar disorder construction, used as a comparison point for hyperscaling violations and Griffiths phases.","marker":"[10]"}],"fun_headline_variants":["3D Ising with critical disorder: new universality","Critical disorder defies Ising universality prediction","Unexpected ν,η from Ising with power-law disorder","New Ising universality class from critical bond disorder","Monte Carlo reveals new Ising exponents for correlated disorder"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis leans on the claim that the random bonds built from critical Ising configurations are power-law correlated with exponent a=d−2+η_pure≈1.036; if the actual bond-bond correlations decay faster or are not characterized by this exponent, the comparison with Weinrib-Halperin theory and the interpretation of νd and ηd would need revision.","fun_headline_variants_meta":{"raw":{"variants":["3D Ising with critical disorder: new universality","Critical disorder defies Ising universality prediction","Unexpected ν,η from Ising with power-law disorder","New Ising universality class from critical bond disorder","Monte Carlo reveals new Ising exponents for correlated disorder"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001009,"raw_usage":{"total_tokens":4265,"prompt_tokens":948,"completion_tokens":3317,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":3238}},"tokens_in":564,"tokens_out":3317,"duration_ms":26478,"temperature":1.0,"reasoning_tokens":3238,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:00:41.405747+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the disorder-averaged bond correlation function C(r)=[(J_ij−J̄)(J_kl−J̄)] as a function of separation in large systems and extract its decay exponent; if it is much larger than 1.036, or if the two finite-size peaks move apart rather than merging when L exceeds 50, the claim of a single new universality class driven by slow power-law disorder is falsified.","supporting_citations":[{"cited_title":"temperature journey","cited_arxiv_id":null,"evidence_quote":"Supplies the Weinrib-Halperin criterion and the predicted fixed-point exponents νd=2/a, ηd=0 that the paper tests."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides high-precision pure Ising exponents (νpure=0.629971(4), ηpure=0.036298(2)) used to set the disorder correlation exponent a=d−2+ηpure."},{"cited_title":"Machta, Population annealing with weighted averages: A Monte Carlo method for rough free-energy landscapes , Phys","cited_arxiv_id":null,"evidence_quote":"Provides the high-precision inverse critical temperature βc,pure of the pure 3D Ising model used to generate critical equilibrium configurations."},{"cited_title":"Wolff, Collective Monte Carlo updating for spin systems , Phys","cited_arxiv_id":null,"evidence_quote":"Supports the claim that restricted and unrestricted disorder ensembles have the same thermodynamic limit, used to argue the two peaks are finite-size effects."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the uncorrelated bond-disorder exponents used as a baseline showing correlations change the universality class."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides a line-defect example where measured exponents sit between pure and Weinrib-Halperin values, a pattern similar to the present result."},{"cited_title":"Three-dimensional universality class of Ising model with power-law-correlated critical disorder","cited_arxiv_id":"1908.01880","evidence_quote":"Reviews the controversy over the Weinrib-Halperin exponents and offers comparison for correlated disorder simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives a 2D Potts model with a similar disorder construction, used as a comparison point for hyperscaling violations and Griffiths phases."}],"review_version":1}