{"id":"931f6d26-8336-4600-b7ea-13c291f41b67","arxiv_id":"1908.01702","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":5,"one_line_summary":"Monte Carlo simulations of the improved Blume-Capel model give the dynamic critical exponent of the 3D Ising universality class as z = 2.0245(15).","lead":"The paper measures a number called the dynamic critical exponent z for the 3D Ising model class, using large-scale Monte Carlo simulations of an improved lattice model. It finds z = 2.0245(15), ending a long disagreement between older simulations and field-theory estimates.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The L^{-2} single-correction ansatz (Eq. 27) may hide opposite-sign subleading corrections (2-η and ω_NR), leaving a systematic bias in z that the quoted error does not cover.","rationale":"The reader's weakest assumption, suppression of leading corrections at D = 0.655, is actually addressed directly in the paper. The synthetic-data check in Section V.B.2 multiplies τ by (1 ± 0.43/30 L^{-ω}) and changes z by only ±0.00045, well within the quoted error. The Binder-cumulant analysis in Appendix D also shows a vanishing leading amplitude at D = 0.655 within precision. So the leading-correction suppression is not the most load-bearing weakness.\n\nThe more serious, under-tested assumption is the form of the correction ansatz. The fits employ a single L^{-2} correction to cover two subleading corrections with exponents 2−η ≈ 1.9637 and ω_NR ≈ 2.0227. The paper itself acknowledges for the susceptibility that these corrections could have opposite signs and cancel over the fitted L-range; the same logic applies to the autocorrelation time. If realized, this cancellation would cause the single-term ansatz to miss a real correction, biasing z, and the quoted error, being an envelope of fits using the same ansatz, would not capture the bias. A two-term fit is a concrete and feasible check. Since the central value is independently supported by the quench analysis and field theory, this concern does not overturn the result, but it does reinforce the need for a conditional verdict and an explicit systematic-error statement.","tokens_in":21973,"tokens_out":15658,"duration_ms":154389,"concrete_test":"Fit the joint Metropolis+HB τ_int,χ data with the two-term ansatz τ = a L^z (1 + c_1 L^{-(2−η)} + c_2 L^{-ω_NR}), fixing 2−η = 1.9637022 and ω_NR = 2.022665, for L_min = 14, 16, 18. Check whether the resulting z is consistent with 2.0245(15) and whether c_1 and c_2 have roughly opposite signs. As a complementary analytical check, generate synthetic data from the same two-term ansatz with c_1 = −c_2 chosen so the corrections cancel near L ≈ 20–30, then fit with the single-term ansatz (27) and compare the recovered z with the input z.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section V.B.2 derives z from fits of τ_int,χ(L) with the single-correction ansatz (27), τ = a L^z (1 + c L^{-ε}) with ε = 2. This term is intended to represent two subleading corrections with nearly degenerate exponents: the analytic background (2−η = 1.9637) and the lattice rotational-symmetry breaking (ω_NR = 2.0227). In Section V.B.1 the author explicitly acknowledges for the susceptibility that 'we can not exclude that these two corrections have amplitudes with opposite sign and cancel to a large extent in the range of lattice sizes considered here.' The same ambiguity applies to τ. If such a cancellation occurs, the single-term ansatz absorbs the combined effect and can leave a residual L-dependent bias in the fitted z. The quoted error of 0.0015 is an envelope over fits that all use this same single-term ansatz, so it does not include the model uncertainty from a possible two-term cancellation. The synthetic-data check in Section V.B.2 only perturbs data with the leading correction amplitude (∝ L^{-ω}, ω = 0.82968), not with the subleading pair. Replacing ε by 2−η or ω_NR individually tests one exponent at a time, not a simultaneous two-term cancellation. Therefore the systematic error on z could exceed the reported 0.0015 if the cancellation scenario is realized.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a high-precision Monte Carlo determination of the dynamic critical exponent z of the three-dimensional Ising universality class for purely dissipative relaxational dynamics (model A). The author simulates the improved Blume-Capel model on the simple cubic lattice at D=0.655, β=0.387721735, using local heat-bath and Metropolis algorithms with checkerboard ordering. Equilibrium integrated autocorrelation times of the magnetic susceptibility are computed for lattice sizes up to L=72 and fit with finite-size scaling ansaetze τ=a L^z (1+c L^{-ε}) with ε=2, as well as without correction terms. The quoted final estimate is z=2.0245(15). Supplementary off-equilibrium quenches from a fully magnetized state to criticality give z=2.0245(10) (Metropolis) and z=2.0240(8) (heat bath), consistent with the equilibrium result. The paper also analyzes the four-loop epsilon expansion in Appendix C, obtaining z=2.0243, and includes simulations of the Ising model and of the Blume-Capel model at D=1.15 to quantify leading corrections to scaling. The central claim is that the new value reconciles Monte Carlo results with recent functional RG and field-theoretic estimates.","tokens_in":22198,"tokens_out":6702,"duration_ms":67343,"significance":"If the result holds, this is the most accurate Monte Carlo determination of the dynamic critical exponent of the three-dimensional Ising universality class for model A dynamics, and it resolves a long-standing tension between older lattice simulations (z≈2.03–2.08) and field-theoretic estimates (z≈2.024). The paper's strengths are the use of an improved Hamiltonian to suppress the leading correction to scaling, the cross-check between two independent local algorithms and between two dynamical protocols, the consistency of the extracted static exponent η with the conformal-bootstrap value, and the explicit synthetic-data tests for residual leading corrections. The appendix D analysis of correction amplitudes for χ, U4, and τ is a useful and nontrivial check of universality of amplitude ratios. The final value is also consistent with functional RG and four-loop epsilon-expansion estimates, which lends independent support.","major_comments":[{"comment":"The final error bar is not robust against the two-correction cancellation that the paper itself flags. In §V.B.1 the author notes for χ that one \"can not exclude that these two corrections have amplitudes with opposite sign and cancel to a large extent in the range of lattice sizes considered here\"; the two corrections are the analytic background (exponent 2−η) and the lattice rotational-symmetry breaking (exponent ω_NR). The same two corrections enter τ, and the fits of τ with ansatz (27) all replace them by a single L^{-2} term. Since every fit used to set z=2.0245(15) employs this same single-term ansatz, the quoted 0.0015 is an envelope over statistical and Lmin variations only; it does not include the model uncertainty from a possible opposite-sign cancellation. The synthetic-data check for residual leading corrections (multiplying by 1±[0.43/30]L^{-ω}) tests only the leading correction amplitude and cannot detect cancellation between the two subleading terms. I ask for a quantitative bound on this scenario, e.g., fits of τ with the two subleading exponents 2−η and ω_NR included simultaneously with independent amplitudes, or with amplitudes constrained by the corresponding χ and U4 analyses; the systematic error should be inflated by the observed spread.","section":"§V.B.2, eqs. (26)–(27)"},{"comment":"The error assignment is not fully documented. The text says the error bar covers fits with ansatz (27) and Lmin=14,16,18 and the Metropolis no-correction fits, but for the heat-bath no-correction fits \"at least the central values are covered\" for Lmin≥40, and that \"completely covering also the error bars of these fits ... seems too pessimistic.\" This is a subjective exclusion of fits that are not statistically rejected: for the heat-bath algorithm with ansatz (26), χ2/d.o.f. drops below one at Lmin=28. Please state the quantitative criterion by which these fits are excluded, or include a systematic term that captures the difference between the no-correction and correction fits; otherwise the quoted error is smaller than the spread of statistically acceptable analyses.","section":"§V.B.2, paragraph containing eq. (28)"}],"minor_comments":[{"comment":"In the definition of the heat-bath probabilities, p(0) is written twice; the third line should be p(+1)=exp(−D+βSx)/z.","section":"Appendix A, eq. (A3)"},{"comment":"\"to a large extend\" should read \"to a large extent\".","section":"§V.B.1"},{"comment":"The inequality symbols appear corrupted (\"t /greaterorapproxeql840\"); they should be rendered as \"≳\".","section":"§VI.A, text near Figure 3"},{"comment":"The two-dimensional boundary condition is enforced as z=2.167 exactly, although the quoted literature values are 2.1665(12) and 2.1667(5). The effect of this rounding on the extracted three-dimensional value is not stated; please give the resulting uncertainty.","section":"Appendix C, eqs. (C4) and (C6)"},{"comment":"The sentence \"Replacing the 2 by 2−η or ω_NR has only little effect on the results for z\" would be more informative if the numerical shifts in z were reported together with the fit ranges for which they were obtained.","section":"§V.B.2, first paragraph"}],"recommendation":"major_revision","confidential_remarks":"The central numerical result appears well supported, and the paper is a valuable contribution, but the error bar needs to be made robust against the subleading-correction ambiguity before publication. The author should be asked to provide a quantitative treatment of the two-correction scenario rather than relying on the single L^{-2} ansatz."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a high-quality numerical paper. Hasenbusch computes the dynamic critical exponent z for 3D Ising model A dynamics using the improved Blume-Capel model at D = 0.655, with two local algorithms, equilibrium autocorrelation times and off-equilibrium quenches, and gets z = 2.0245(15). The result is consistent with functional RG and 4-loop epsilon expansion, and roughly an order of magnitude more precise than the previous improved-model simulation (z = 2.020(8)). I believe the central value.\n\nWhat the paper does well: the simulation campaign is large (about 8.4 core-years total), the author checks algorithm dependence, lattice-size dependence, residual leading corrections via synthetic data, and explicitly analyzes subleading corrections (2 − η and ω_NR). The quench analysis uses variance reduction and gives consistent values with smaller statistical errors. Appendix C's field-theoretic analysis is a reasonable bonus. The paper is transparent about places where assumptions are made.\n\nSoft spots: (1) The quoted error is an envelope over fits that all use a single L^{-2} correction term. The two subleading exponents, 2 − η ≈ 1.964 and ω_NR ≈ 2.023, are within 3% of each other, and the author admits in the susceptibility analysis that opposite-sign amplitudes could cancel in the studied L range. The same possibility applies to the autocorrelation time. Replacing ε = 2 with either exponent individually does not test simultaneous cancellation, and the synthetic-data checks only perturb the leading correction amplitude, not this subleading pair. So the quoted ±0.0015 is a statement about fitting variability, not a rigorous bound on this model uncertainty. I don't think this is fatal — the quench results and field theory agree, so a large bias is unlikely — but it should be addressed explicitly. (2) Reproducibility: raw data are not tabulated and source code is not provided. For a precision metrology paper, that is a weakness, though not a reason to reject. (3) D* is taken from the author's earlier work rather than re-derived; that is fine since it is a consistency anchor, not an input to z.\n\nWho it is for: statistical mechanics people working on critical dynamics, and anyone using improved models for Monte Carlo. It deserves a serious referee. I would recommend acceptance after minor revision, with the author asked to provide the data or a more formal estimate of the subleading-correction systematic.","headline":"A careful, high-statistics Monte Carlo determination of z for 3D Ising model A dynamics that gives z = 2.0245(15) and reconciles simulation with field theory; the central value is credible, but the quoted error is an envelope over one-correction fits and may leave a small systematic uncovered.","tokens_in":22779,"tokens_out":3966,"would_cite":true,"duration_ms":37233,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B20","82B27","82B80"],"pacs":["05.10.Ln","64.60.Fr","75.10.Hk"],"model":"deepseek-v4-flash","headline":"Precise Monte Carlo measurement puts the 3D Ising dynamic critical exponent at z = 2.0245(15).","keywords":["dynamic critical exponent","Ising universality class","Blume-Capel model","critical slowing down","autocorrelation time","Monte Carlo simulation","model A dynamics","finite-size scaling"],"falsifier":"Measure $\\tau_{\\mathrm{int},\\chi}$ on larger lattices, say $L = 96$ and $128$, with errors on $z$ below $5\\times10^{-4}$, and refit with the same $L^{-2}$ correction form: if the exponent drifts by more than the quoted $1.5\\times10^{-4}$ as $L_{\\min}$ increases, the correction ansatz is incomplete. Alternatively, take a second improved coupling inside $D^*=0.656(20)$, such as $D=0.656$; if the two determinations of $z$ disagree beyond error bars, the assumed suppression of leading corrections is wrong.","tokens_in":21700,"feed_emoji":"⏱","tokens_out":8112,"duration_ms":72597,"temperature":0.7,"pith_summary":"The paper aims to pin down the dynamic critical exponent $z$ of the three-dimensional Ising universality class for purely dissipative relaxational dynamics (model A), the power law that governs how autocorrelation times grow near criticality. The author simulates the improved Blume-Capel model on the simple cubic lattice, tuning the coupling to $D = 0.655$ so that the dominant correction to scaling is suppressed by at least a factor of 30. Finite-size scaling of the equilibrium integrated autocorrelation time of the magnetic susceptibility yields $z = 2.0245(15)$, and sudden quenches of fully magnetized configurations to the critical point give consistent values. The result matters because most Monte Carlo estimates for the plain Ising model were not compatible with field-theoretic calculations, and this measurement closes that gap.","feed_headline":"3D Ising dynamics pinned at z = 2.0245(15)","feed_subtitle":"Improved Blume-Capel simulations reconcile Monte Carlo and field-theory values for critical slowing down.","key_machinery":"The carrying mechanism is the improved Blume-Capel model at the parameter $D^* = 0.656(20)$ (simulated at $D = 0.655$), where the amplitude of the leading correction to scaling is suppressed by at least a factor of 30 relative to the spin-$1/2$ Ising model, so the remaining corrections to $\\tau \\sim L^z$ can be represented by a single effective $L^{-2}$ term. The measured observable is the integrated autocorrelation time of the magnetic susceptibility, evaluated from autocorrelation functions with a self-consistent truncation and single-exponential tail correction, and then fitted with and without the $L^{-2}$ correction over lattice sizes up to $L = 72$. The heat bath and Metropolis algorithms give consistent exponents, and a combined fit yields the central value.","core_discovery":"The central claim is that, at the improved point $D = 0.655$, the integrated autocorrelation time of the magnetic susceptibility at criticality obeys $\\tau_{\\mathrm{int},\\chi} = a L^z (1 + c L^{-\\epsilon})$ with an effective correction exponent $\\epsilon = 2$, and fits to this form across lattice sizes $L \\ge 14$ give $z = 2.0245(15)$ for the model-A dynamic critical exponent. The same value is obtained from out-of-equilibrium quenches, where the magnetization decays as $m(t) = a (t - t_0)^{-\\beta/(\\nu z)}$. This estimate is fully consistent with functional renormalization group results and with the four-loop $\\epsilon$-expansion analyzed in the paper, and it attributes the earlier spread of Monte Carlo values to unsubtracted leading corrections to scaling.","pith_inferences":["Beyond the paper: the same improved-model strategy could be applied to other dynamics (e.g., conserved order parameter, model B) to test whether a single static universality class fixes the dynamic exponent once corrections are removed, or whether the conservation law changes it as expected.","Beyond the paper: applying this correction-free analysis to the two-dimensional Ising class might sharpen the accepted value $z = 2.1665(12)$ and test whether the one-term $L^{-2}$ ansatz remains adequate at higher precision.","Beyond the paper: the reported correction amplitudes predict that re-analyzing published Ising-model autocorrelation data with a leading $L^{-\\omega}$ term of amplitude $a_\\tau \\approx -0.44(3)$ should bring those older estimates down to $z \\approx 2.024$; this is directly checkable with existing data."],"forward_implications":["Earlier high Monte Carlo estimates of $z$ for the 3D Ising model, such as $z \\approx 2.04$ to $2.08$, can be explained as uncorrected leading corrections to scaling, whose amplitude is particularly large for autocorrelation times.","The dynamic critical exponent of the 3D Ising universality class for model A dynamics is now consistent across Monte Carlo, functional renormalization group, and resummed four-loop $\\epsilon$-expansion results, all near $z = 2.024$.","The universal ratio of correction amplitudes for the autocorrelation time relative to the Binder cumulant, $a_\\tau / a_U = 3.1(6)$, provides a quantitative target for other models in the same universality class.","Equilibrium and off-equilibrium determinations of $z$ agree, supporting the transferability of the exponent across different dynamical protocols within model A dynamics."],"supporting_citations":[{"why":"provides the improved coupling D*=0.656(20), the critical inverse temperature, and the factor-30 suppression of leading corrections that the whole analysis relies on","marker":"[37]"},{"why":"supplies the accurate static exponents nu, eta, omega and correction exponents used as inputs in the fitting ansaetze","marker":"[1]"},{"why":"gives the functional renormalization group estimate z≈2.025 that the simulation result is compared with","marker":"[13]"},{"why":"reports functional RG values z=2.024, 2.024, and 2.023 using frequency regulators, the field-theoretic target","marker":"[14]"},{"why":"provides the four-loop epsilon-expansion series whose reanalysis in appendix C yields z≈2.0243","marker":"[15]"},{"why":"is the previous improved Blume-Capel off-equilibrium simulation with z=2.020(8), the direct predecessor that this work extends","marker":"[32]"},{"why":"is an equilibrium Ising-model Monte Carlo estimate z=2.03(4) representing the older results that disagreed with field theory","marker":"[24]"}],"fun_headline_variants":["3D Ising z = 2.0245(15) via improved Blume-Capel","Monte Carlo matches field theory: 3D Ising z = 2.0245(15)","Improved Blume-Capel yields 3D Ising dynamic exponent 2.0245(15)","3D Ising critical slowing down z = 2.0245(15) from MC","Dynamic exponent z=2.0245(15) from improved 3D Ising model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that at $D = 0.655$ the leading correction to scaling is suppressed by at least a factor of 30, so all remaining corrections to $\\tau \\sim L^z$ can be absorbed into a single effective $L^{-2}$ term; the improved point itself is taken from the author's earlier work and is not re-derived in this paper.","fun_headline_variants_meta":{"raw":{"variants":["3D Ising z = 2.0245(15) via improved Blume-Capel","Monte Carlo matches field theory: 3D Ising z = 2.0245(15)","Improved Blume-Capel yields 3D Ising dynamic exponent 2.0245(15)","3D Ising critical slowing down z = 2.0245(15) from MC","Dynamic exponent z=2.0245(15) from improved 3D Ising model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00115,"raw_usage":{"total_tokens":4717,"prompt_tokens":842,"completion_tokens":3875,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":458,"completion_tokens_details":{"reasoning_tokens":3748}},"tokens_in":458,"tokens_out":3875,"duration_ms":23781,"temperature":1.0,"reasoning_tokens":3748,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:05:48.623364+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $\\tau_{\\mathrm{int},\\chi}$ on larger lattices, say $L = 96$ and $128$, with errors on $z$ below $5\\times10^{-4}$, and refit with the same $L^{-2}$ correction form: if the exponent drifts by more than the quoted $1.5\\times10^{-4}$ as $L_{\\min}$ increases, the correction ansatz is incomplete. Alternatively, take a second improved coupling inside $D^*=0.656(20)$, such as $D=0.656$; if the two determinations of $z$ disagree beyond error bars, the assumed suppression of leading corrections is wrong.","supporting_citations":[{"cited_title":"Wansleben and D","cited_arxiv_id":null,"evidence_quote":"provides the improved coupling D*=0.656(20), the critical inverse temperature, and the factor-30 suppression of leading corrections that the whole analysis relies on"},{"cited_title":"Below we analyze the merged re sults","cited_arxiv_id":null,"evidence_quote":"supplies the accurate static exponents nu, eta, omega and correction exponents used as inputs in the fitting ansaetze"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"reports functional RG values z=2.024, 2.024, and 2.023 using frequency regulators, the field-theoretic target"},{"cited_title":"Diagram Reduction in Problem of Critical Dynamics of Ferromagnets: 4-Loop Approximation","cited_arxiv_id":"1712.05917","evidence_quote":"provides the four-loop epsilon-expansion series whose reanalysis in appendix C yields z≈2.0243"},{"cited_title":"In particular all studies that are performed later than 1991 are not compatible within the quoted error bars with eq","cited_arxiv_id":null,"evidence_quote":"is the previous improved Blume-Capel off-equilibrium simulation with z=2.020(8), the direct predecessor that this work extends"},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"is an equilibrium Ising-model Monte Carlo estimate z=2.03(4) representing the older results that disagreed with field theory"}],"review_version":1}