{"id":"e9a594c9-4aec-48fa-8b41-43328ab8d6a7","arxiv_id":"2411.10266","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The paper extracts universal non-equilibrium scaling functions for the order parameter and its cumulants up to fourth order in a Z2 scalar field theory with Model A dynamics in two and three dimensions.","lead":"This paper uses large-scale computer simulations of a simple magnet-like field theory to extract universal curves describing how order parameter fluctuations, up to fourth order, evolve when the system is quenched through a critical point. These curves are intended as building blocks for models of the QCD critical point search in heavy-ion collision experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Universal-function extraction lacks a quantitative corrections-to-scaling test; reported Pade fits absorb unquantified window and z systematics, so the 'fully describe' claim is stronger than the evidence.","rationale":"The reader's weakest assumption identifies the scaling ansatz (19) holding with negligible corrections as the most vulnerable point, and I agree. My stress-test sharpens this: the paper provides no quantitative measure of collapse quality, no propagation of the uncertainty in z into the reported universal functions, and no test of subleading scaling corrections. The authors themselves note in Appendix B that the dynamic scaling region is selected by eye and that the d=3 z estimate is low, which underscores that the extracted functions may carry systematic errors beyond the quoted statistical uncertainties. The finite-size scaling functions are given only as spline wireframes, so the two-variable claim is even less quantitatively supported. These concerns do not overturn the central result, since the visual collapse in Figs. 3 and 4 is a nontrivial test of the scaling hypothesis, but they justify the CONDITIONAL verdict: the strength of the 'fully describe' statement exceeds what the current evidence demonstrates. My concrete test—a fast/slow subset split for d=3 f_chi—would directly reveal whether the reported function is stable against the choice of scaling window, and repeating with the authors' own z would quantify the exponent sensitivity. If both checks pass, the universality claim is substantially strengthened; if they fail, the reported fits must be reinterpreted as effective functions with the caveats stated. The verdict should remain CONDITIONAL pending these checks, matching the reader's assessment.","tokens_in":18810,"tokens_out":12119,"duration_ms":119486,"concrete_test":"Split the quench rates used in the Fig. 3/4 collapse for d=3 into two disjoint subsets (fast half vs slow half) within the nominal scaling window. Using the same literature exponents and amplitudes, extract f_chi(x) independently from each subset and compare over the overlapping x-range, with bootstrap errors. If the two estimates differ by more than the combined statistical uncertainty, corrections to scaling are present and the reported f_chi is an effective window-dependent fit rather than the asymptotic universal function. Repeat the same extraction using the authors' own z (1.949) to test sensitivity to the input exponent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the Pade fits in Tables A.3/A.4 are universal functions that fully describe the non-equilibrium evolution near the critical point. This requires the scaling ansatz (19) to hold in the probed (r_J, L) window with negligible corrections to scaling. The paper does not establish this quantitatively. Appendix B documents Kibble-Zurek breakdowns at slow and fast quenches and shows that the dynamic scaling region for the z extraction is identified by eye from local slope plateaus (Fig. B.7). No chi-square or residual analysis is reported for the collapses in Figs. 3 and 4; the fits quote only statistical errors. Two systematics are therefore unquantified: (i) the uncertainty in the input dynamic exponent z, where the authors' own d=3 estimate (1.949 +/- 0.054) is about 1.4 sigma below the literature value (2.0245 +/- 0.0015) used for the collapse, and (ii) corrections to scaling from the finite range of quench rates and lattice sizes. Since x = J/J_KZ depends on z through 1/(1+nu_c z), a shift in z or rate-dependent corrections changes the extracted f_A. If the effective z in the simulation window differs from the asymptotic value, or if subleading corrections are non-negligible, the reported functions are effective, window-dependent fits rather than universal ones. The finite-size extension (Sec. 4.3) is presented only as spline wireframes in Fig. 5 without quantitative collapse metrics, so the two-variable claim is even less constrained.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports classical-statistical lattice simulations of a Z2-symmetric scalar field theory with Langevin dynamics (Model A) in d=2 and d=3, driven through the critical point by a linearly varying external field J(t) = -r_J t at T = T_c. Using Kibble-Zurek scaling, the authors rescale the magnetization, susceptibility, skewness, and kurtosis and observe data collapse onto universal scaling functions f_M, f_chi, f_kappa3, f_kappa4, which are parametrized by Padé approximants given in Tables A.3 and A.4. The analysis is extended to a two-variable scaling function f_A(x, xi_KZ/L) for finite system sizes in Sec. 4.3, and the dynamic critical exponent z is extracted from the scaling of the zero crossing of the magnetization in Appendix B. The abstract claims that, together with the static and dynamic critical exponents, these scaling functions fully describe the universal non-equilibrium evolution of the system near the critical point.","tokens_in":19150,"tokens_out":6662,"duration_ms":61639,"significance":"If the extracted scaling functions are truly universal, the paper is a valuable contribution: it provides closed-form expressions for the out-of-equilibrium evolution of order-parameter cumulants up to fourth order in a canonical dynamic universality class, with direct relevance to phenomenological modeling near the QCD critical point. The dataset is substantial (about 10^4 and 5 x 10^3 independent runs per quench rate in 2D and 3D), statistical uncertainties are bootstrap-estimated, and the Padé parametrizations make the results usable by other groups. The paper also usefully catalogs regimes where Kibble-Zurek scaling breaks down. However, the central universality claim rests on visually assessed collapse; quantitative collapse metrics, a corrections-to-scaling analysis, and a full propagation of the z uncertainty are missing. The finite-size extension is presented only as a spline wireframe without quantitative validation. These gaps currently limit the claim to effective scaling functions valid in the probed window rather than fully established universal functions.","major_comments":[{"comment":"The collapse of the rescaled data is assessed visually only; no chi-square, residual analysis, or other quantitative metric is reported for the quality of the collapse or for the Padé fits. Given the central claim that the Padé functions in Tables A.3 and A.4 fully describe the universal non-equilibrium evolution, please add a quantitative goodness-of-fit measure (e.g., chi^2 per degree of freedom against the collapsed data), report the exact x-intervals used for each fit, and test the stability of the fits when the fitting range or the Padé order is varied. Without this, the reported functions may be effective, window-dependent fits that absorb corrections to scaling and systematics from the input exponent z.","section":"Sec. 4.1, Figs. 3 and 4, Tables A.3 and A.4"},{"comment":"The rescaling in Sec. 4.1 uses the literature values z = 2.1667(5) (d=2) and z = 2.0245(15) (d=3), while the authors' own estimates from the same data are z = 2.142(49) and z = 1.949(54), with the d=3 value about 1.4 sigma below the literature value. Since the scaling variable x = r_J^{-1/(1+nu_c z)} J depends on z, the extracted f_A inherit this uncertainty. Please quantify how the collapsed curves and Padé coefficients change when z is varied within the literature uncertainty, for example by recomputing the collapse at the 1-sigma bounds of z. Alternatively, treat z as a free parameter in a joint fit of the collapse and report the resulting f_A and confidence intervals. This is essential for assessing whether the reported functions are universal rather than specific to the chosen input z.","section":"Appendix B and Sec. 4.1"},{"comment":"The two-variable finite-size scaling function f_A(x, xi_KZ/L) is presented only as a polyharmonic spline wireframe, with no quantitative test that the rescaled data for different lattice sizes and quench rates collapse onto a single surface. Please provide a collapse metric for fixed values of xi_KZ/L (e.g., overlay slices obtained from different (L, r_J) combinations that yield the same xi_KZ/L and report the spread), and specify the reference scale L0 that defines the dimensionless size L = L/L0 in Eq. (31), including how L0 is determined from equilibrium finite-size scaling and its numerical value. Without this, the claim of a universal finite-size scaling function is not quantitatively established.","section":"Sec. 4.3 and Fig. 5"}],"minor_comments":[{"comment":"The statement that the scaling functions 'fully describe the universal non-equilibrium evolution' is stronger than what is demonstrated, given the documented breakdowns in Appendix B and the absence of a corrections-to-scaling analysis; please soften the wording or explicitly restrict the claim to the dynamic scaling region and the fitted range of x.","section":"Abstract and Sec. 5"},{"comment":"The expression for the time derivative of the relaxation time appears garbled: '(-t)^{-1}nu_c z' is likely intended in Eq. (15); please correct the typography.","section":"Sec. 3.1, Eqs. (14)-(15)"},{"comment":"The normalization condition f_A(1,0) = f_A(0,1) = 1 is stated in Eq. (20), but only f_M(0) = 1 is actually used to fix r_J0; please clarify whether the normalization conditions for chi, kappa3, and kappa4 are also enforced or whether their overall amplitudes are fixed solely by the literature amplitudes A0.","section":"Sec. 4.1, Eq. (20)"},{"comment":"The dynamic scaling region is identified 'by eye' from local slope plateaus; please state an objective criterion for the plateau selection and list the data points and quench-rate ranges used in the fit.","section":"Appendix B, Fig. B.7"},{"comment":"The conclusions refer to '2+1D' and '3+1D' while the simulations are in two and three spatial dimensions; please use consistent notation throughout the paper.","section":"Sec. 5"},{"comment":"The code and data are described but no repository or data availability statement is provided; please consider making the code and data publicly available to facilitate reproducibility.","section":"Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of Nuclear Physics B and builds on the authors' earlier work [41]. The main gap is between the strong 'fully describe' claim and the evidence: the collapse is visually convincing but not quantitatively characterized, the d=3 dynamic exponent extraction deviates from the literature value used for the collapse, and the finite-size scaling section lacks quantitative collapse metrics. With the requested additions (collapse goodness-of-fit, z-uncertainty propagation, and a finite-size collapse test), the paper would meet the journal's standards. I see no reason to doubt the honesty or technical quality of the simulations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth your time. It gives the first high-statistics extraction of universal non-equilibrium scaling functions for M, chi, kappa3, and kappa4 in the 2D and 3D Model A universality class, using classical-statistical lattice simulations where critical slowing down emerges naturally rather than being imposed. The data collapse across quench rates is visually convincing, and the two-variable finite-size extension is a genuinely useful addition. The Padé fits in Appendix A are practical tools for phenomenology. The authors also deserve credit for an unusually honest Appendix B: they document both breakdown regimes, report their own z estimate, and explicitly note the arbitrariness in defining the dynamic scaling region.\n\nThe soft spots are real but not fatal. The abstract says the scaling functions 'fully describe' the universal non-equilibrium evolution, which is stronger than what is actually shown. There is no residual analysis or chi-square for the collapses, no quantitative corrections-to-scaling test, and the Padé fits inherit unquantified systematics from the chosen z and the by-eye selection of the scaling region. The d=3 z estimate sits about 1.4 sigma below the high-precision Monte Carlo value; the authors attribute this to finite-size and rate limitations, which is plausible but not demonstrated. The finite-size scaling surface is only presented as spline wireframes without a quantitative collapse metric. And there is no released code or data, which would materially help independent verification. The stress-test note is largely on target here, though I would not call the central claim unsupported: within the probed window, the collapse genuinely demonstrates Kibble-Zurek scaling, and the limitations are acknowledged rather than hidden.\n\nThis is a benchmark paper for Model A quenches, directly useful for people connecting critical dynamics to heavy-ion phenomenology or testing similar methods in other universality classes. The conditional verdict from your reader is the right one: strong and worth publishing, but it should be tightened with code/data release, residual analysis, and more careful wording. Send it to peer review.","headline":"Solid numerical extraction of universal cumulant scaling functions for Model A; the claims are slightly ahead of the evidence, but the central result holds and deserves a referee.","tokens_in":19700,"tokens_out":1620,"would_cite":true,"duration_ms":17866,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that universal Kibble–Zurek scaling functions, fitted here with Padé approximants, completely describe the non-equilibrium evolution of the order parameter and its first four cumulants in linear magnetic quenches through…","keywords":["dynamic critical phenomena","non-equilibrium phase transitions","classical-statistical simulations","Kibble-Zurek scaling","Model A universality class","cumulants","finite-size scaling","lattice field theory"],"falsifier":"One decisive test is to run the same Model A dynamics at a quench rate inside the identified dynamic scaling region but with a different lattice coupling or dissipation strength, and check whether the rescaled curves still collapse onto the Padé fits of Tables A.3 and A.4; systematic deviation would show the functions are not fully universal. Alternatively, measure the local slope of the magnetization zero-crossing over a wider range of quench rates and verify that it plateaus at $1/(1+\\nu_c z)$ with the same value used to extract $z$ in this paper.","tokens_in":18599,"feed_emoji":"⚛️","tokens_out":14647,"duration_ms":119393,"temperature":0.7,"pith_summary":"This paper aims to establish that the full out-of-equilibrium response of a $\\mathbb{Z}_2$-symmetric scalar field theory driven linearly through its critical point is a universal function of a single scaling variable. Using classical-statistical lattice simulations of Model A dynamics in $d=2$ and $d=3$, the authors show that the magnetization, susceptibility, skewness, and kurtosis all collapse onto rate-independent curves $f_M(x)$, $f_\\chi(x)$, $f_{\\kappa_3}(x)$, $f_{\\kappa_4}(x)$ when plotted against $x = J/J_{\\mathrm{KZ}}$, where $J_{\\mathrm{KZ}}$ is the Kibble–Zurek field scale set by the quench rate. The claim matters because searches for the QCD critical point in heavy-ion collisions need quantitative predictions for how critical fluctuations evolve out of equilibrium; these functions, together with the static exponents and the dynamic exponent $z$, would supply such predictions without further dynamical simulation. The paper further claims that adding a second scaling variable, the Kibble–Zurek length divided by the system size $\\xi_{\\mathrm{KZ}}/L$, extends the description to finite systems and yields two-dimensional universal scaling functions.","feed_headline":"One scaling variable predicts quench dynamics at a critical point","feed_subtitle":"Lattice simulations fit universal curves for magnetization, susceptibility, skewness and kurtosis in 2D and 3D.","key_machinery":"The load-bearing object is the Kibble–Zurek finite-time scaling ansatz, Eq. (19): every observable $A$ with scaling dimension $\\Delta_A$ obeys $A(J,r_J) = A_0\\, s^{\\Delta_A} f_A(s^{1/\\nu_c} J,\\, s^{z+1/\\nu_c} r_J)$, where $s$ is a length-rescaling parameter and $r_J$ is the dimensionless quench rate. Choosing $s = r_J^{-\\nu_c/(1+\\nu_c z)}$ eliminates the quench-rate argument and leaves a function of the single variable $x = J/J_{\\mathrm{KZ}}$ with $J_{\\mathrm{KZ}} \\sim r_J^{1/(1+\\nu_c z)}$. The cumulants enter only through their equilibrium scaling dimensions $\\Delta_n = -\\beta/\\nu + (n-1)/\\nu_c$, with $\\nu_c=\\nu/(\\beta\\delta)$, which sets how each observable is rescaled before the collapse is tested. The finite-size extension, Eq. (31), adds a third argument $s L^{-1}$ and produces the two-variable functions $f_A(x,\\xi_{\\mathrm{KZ}}/L)$ that govern the crossover from quench-limited to system-size-limited behavior.","core_discovery":"For a scalar field theory with $\\mathbb{Z}_2$ symmetry in the relaxational Model A universality class, a constant-rate quench of the symmetry-breaking external field through the critical point produces non-equilibrium evolution of the order parameter and its cumulants that is described by universal Kibble–Zurek scaling functions. In $d=2$ and $d=3$, the rescaled data for the average magnetization, susceptibility, skewness, and kurtosis collapse onto single curves $f_M(x)$, $f_\\chi(x)$, $f_{\\kappa_3}(x)$, and $f_{\\kappa_4}(x)$ with $x=J/J_{\\mathrm{KZ}}$, and the paper provides closed-form Padé fits for these functions in Tables A.3 and A.4. The same scaling framework yields an independent estimate of the dynamic critical exponent from the quench-rate scaling of the magnetization zero-crossing, namely $z=2.142(49)$ in $d=2$ and $z=1.949(54)$ in $d=3$, consistent with Monte Carlo values in $d=2$ and slightly lower in $d=3$. With the system size included as an additional scaling variable, the data collapse onto two-variable functions that interpolate between infinite-volume non-equilibrium scaling and equilibrium finite-size scaling, confirming that the description remains valid when $\\xi_{\\mathrm{KZ}}/L$ is no longer small.","pith_inferences":["One could test the universality claim more stringently by changing the lattice coupling or the dissipation strength and checking whether the same Padé fits describe the rescaled data; the paper reports only a brief check of the coupling and found no improvement.","The same finite-time scaling machinery should apply to other control-parameter paths, such as temperature quenches or cis-critical protocols, and comparing the extracted functions across paths would show which features of the universal curves are genuinely universality-class properties rather than artifacts of the linear-field protocol.","The marked asymmetry of the susceptibility during the quench—a gradual departure from equilibrium before the critical point and a sharp return after—is a dynamical signature that an effective kinetic description would have to reproduce; it is not fixed by equilibrium critical exponents alone.","Mapping the 3D Ising scaling functions onto QCD phase-diagram trajectories would give concrete predictions for skewness and kurtosis at finite beam energies; if the QCD critical point is in Model H rather than Model A, the two sets of functions would differ in ways that could be used to identify the correct dynamical universality class."],"forward_implications":["Given the static critical exponents and the dynamic exponent $z$, the published Padé fits predict the time evolution of the magnetization, susceptibility, skewness and kurtosis for any linear magnetic quench through the critical point in $d=2$ and $d=3$ Model A dynamics, with no further simulation needed.","For finite systems, the two-variable scaling functions show that corrections become visible once $\\xi_{\\mathrm{KZ}}/L \\gtrsim 0.5$: the transition softens, cumulant peaks shrink, and the crossing shifts toward $J=0$, eventually reproducing equilibrium finite-size scaling in the adiabatic limit.","The Kibble–Zurek-based estimates of the dynamic exponent provide an independent cross-check of $z$ that is consistent with Monte Carlo results in $d=2$ and slightly lower in $d=3$, signaling the systematic limitations of finite-size and fast-quench simulations.","Because the scaling functions are universal, they provide a concrete benchmark for comparing more realistic dynamical frameworks, such as Model H stochastic fluid dynamics in heavy-ion phenomenology, against the predicted cumulant evolution."],"supporting_citations":[{"why":"It measures the non-universal critical amplitudes used to rescale the raw simulation data into universal form.","marker":"[41]"},{"why":"It defines the Kibble–Zurek scaling limit and the protocol classification on which the scaling ansatz rests.","marker":"[53]"},{"why":"It provides the exact static critical exponents for two dimensions used in the rescaling.","marker":"[55]"},{"why":"They supply the high-precision static critical exponents for three dimensions used in the rescaling.","marker":"[56, 57]"},{"why":"It supplies the two-dimensional dynamic critical exponent input used for the data collapse.","marker":"[60]"},{"why":"It supplies the three-dimensional dynamic critical exponent input used for the data collapse.","marker":"[61]"},{"why":"It is the earlier Model A Langevin study whose scaling functions, built with an imposed dynamic exponent, are contrasted with the emergent-z results in this paper.","marker":"[27]"},{"why":"It established universal off-equilibrium scaling of higher-order cumulants, the observable target extended here.","marker":"[34]"},{"why":"It documents the breakdown of Kibble–Zurek scaling in fast quenches, delimiting the window in which the present scaling functions are claimed to hold.","marker":"[74]"}],"fun_headline_variants":["Universal quench scaling for cumulants up to fourth order","Cumulant curves collapse across 2D and 3D quenches","Kibble-Zurek scaling predicts all order-parameter cumulants","Finite-size effects fit into universal quench curves","One scaling variable describes non-equilibrium criticality"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that in the chosen window of quench rates and lattice sizes the Kibble–Zurek scaling form holds with negligible corrections, so that every observable depends on the quench rate only through the combination $J/J_{\\mathrm{KZ}}$ (and, for finite systems, through $\\xi_{\\mathrm{KZ}}/L$); the paper's Appendix B documents how very slow and very fast quenches break this behavior.","fun_headline_variants_meta":{"raw":{"variants":["Universal quench scaling for cumulants up to fourth order","Cumulant curves collapse across 2D and 3D quenches","Kibble-Zurek scaling predicts all order-parameter cumulants","Finite-size effects fit into universal quench curves","One scaling variable describes non-equilibrium criticality"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000912,"raw_usage":{"total_tokens":3935,"prompt_tokens":976,"completion_tokens":2959,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":592,"completion_tokens_details":{"reasoning_tokens":2873}},"tokens_in":592,"tokens_out":2959,"duration_ms":22227,"temperature":1.0,"reasoning_tokens":2873,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:47:31.685352+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"One decisive test is to run the same Model A dynamics at a quench rate inside the identified dynamic scaling region but with a different lattice coupling or dissipation strength, and check whether the rescaled curves still collapse onto the Padé fits of Tables A.3 and A.4; systematic deviation would show the functions are not fully universal. Alternatively, measure the local slope of the magnetization zero-crossing over a wider range of quench rates and verify that it plateaus at $1/(1+\\nu_c z)$ with the same value used to extract $z$ in this paper.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It defines the Kibble–Zurek scaling limit and the protocol classification on which the scaling ansatz rests."},{"cited_title":"Schweitzer, Dynamic Critical Phenomena in the Classical Approximation on a Lattice, Doctoral dissertation, Justus-Liebig-University Giessen (2021)","cited_arxiv_id":null,"evidence_quote":"It provides the exact static critical exponents for two dimensions used in the rescaling."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the two-dimensional dynamic critical exponent input used for the data collapse."}],"review_version":1}