{"id":"c5aad8ba-d683-498c-86ac-d2fc41be9e4d","arxiv_id":"2411.15994","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Direct numerical simulations of stochastic model H give a dynamic critical exponent z ≈ 3 in 3D and z ≈ 2 in 2D, with a crossover from mean-field z = 4 controlled by the renormalized shear viscosity.","lead":"Using a Metropolis-based lattice algorithm for stochastic model H hydrodynamics, the paper reports direct simulations of critical fluctuations and extracts the dynamic critical exponent z. The results confirm z around 3 in three dimensions and around 2 in two dimensions, giving numerical backing for the critical dynamics expected near a possible QCD critical endpoint.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"2D model H0 may not have the same universal dynamics as model H, so the central 2D claim rests on an untested universality assumption.","rationale":"The reader identified exactly the same weakest point: the 2D exponent is measured on model H0, not model H, based on an assumed shared universality class. I agree that this is the most load-bearing assumption of the paper. The 3D result is substantially safer: Fig. 9 shows model H and model H0 data and demonstrates that the effective exponent collapses as a function of the renormalized viscosity, and the extracted z = 3.013 ± 0.058 is consistent with the established model H exponent. In 2D, no such cross-check exists. The theoretical argument for dropping self-advection is parametric at the level of linear damping rates, but near the critical point with z = 2 the two damping rates are parametrically comparable, so the truncation is not guaranteed to be irrelevant. There is also a subtlety specific to d=2: the shear viscosity itself has a logarithmic infrared issue, and the ε-expansion and FRG disagree about xη, so the expectation that model H0 and model H share the same static/dynamic universality class is less secure than in d=3. I would therefore keep the verdict as CONDITIONAL rather than upgrading to a clean accept: the 2D headline is contingent on a checkable, unperformed simulation. I do not see grounds for REJECT; the algorithm, static checks, 3D result, and data deposition all support the paper's main methodological claims, and the 2D concern is a well-posed empirical question rather than an internal inconsistency.","tokens_in":22420,"tokens_out":1909,"duration_ms":15382,"concrete_test":"Run the d=2 full model H simulation (including the self-advection term, Eq. 42) at η = 10^-2 with L = 40 and 48, and extract z by the same data-collapse procedure used in Fig. 11. If the full model H z is consistent with 2.11 ± 0.015, the concern is settled; if it shifts toward 2.18-2.2 (the ε-expansion value) or moves by more than ~0.05, the model H0 truncation does change the universality class and the 2D headline must be revised. A cheaper supporting check: re-extract z from model H0 using the lowest non-zero momentum mode (n=1) and a third volume (L=64) to quantify finite-size and mode-selection systematics.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's headline result for d=2 is z ≈ 2. This is extracted exclusively from model H0, defined as model H minus the π^T self-advection term. Section V.D does not report a full model H simulation in d=2; the only d=2 dynamic collapse is shown for model H0 (Fig. 11). The justification for replacing model H by model H0 is the statement in Section II.B that the self-coupling is negligible because shear modes decay as (η/ρ)k^2 while order-parameter fluctuations decay as Γk^4. However, in two dimensions shear fluctuations are enhanced at low momentum, and even non-critically the d=2 shear mode has subtle infrared behavior (the zero-frequency viscosity has a log divergence; Refs. [9,74,75]). For z ≈ 2, the ratio of damping rates (Γk^4)/( (η/ρ)k^2 ) behaves as k^{z-2} = k^0 at the critical fixed point, so the two mode-coupling terms are marginal relative to each other precisely in d=2. Therefore the truncation is not obviously innocuous at the renormalization-group level, and a z measured on model H0 could differ from model H. Moreover, the quoted 2D uncertainty (z = 2.11 ± 0.015) is extremely small and is derived from only two volumes (L=40, 48) and one momentum mode (second non-trivial mode); the statistical error does not include systematics of the H→H0 truncation or of finite-L scaling. This makes the 2D claim substantially weaker than the 3D claim, where a full model H curve exists and agrees with model H0 once plotted against renormalized viscosity.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a Metropolis-based lattice algorithm for stochastic model H hydrodynamics, in which the dissipative update and noise are combined so that the equilibrium distribution is the Gibbs distribution of the microscopic Hamiltonian. The method is applied to a conserved Ising-type order parameter coupled to transverse momentum density in d=3 and d=2. Static checks (magnetization histograms, subvolume scaling, correlation functions) are presented, the momentum correlation function is used to extract a renormalized shear viscosity, and the order-parameter correlation function is used to extract effective dynamic critical exponents via finite-size dynamic scaling with L=40 and L=48. In three dimensions the authors report z_eff=3.013±0.058 at ηR=10^-2, with a crossover from z≈4 to z≈3 that they compare with the Kawasaki approximation. In two dimensions, using only model H0 (model H without the transverse momentum self-advection), they report z=2.11±0.015 at η=10^-2 and argue for an asymptotic value z≈2.","tokens_in":22836,"tokens_out":5155,"duration_ms":52257,"significance":"If the central claims hold, the paper provides a direct numerical test of model H dynamic scaling and demonstrates a practical method for stochastic fluid dynamics near a critical point, which is relevant for QCD critical-endpoint phenomenology. The Metropolis construction is attractive because fluctuation-dissipation relations are satisfied by construction, the static equilibrium behavior can be checked directly, and the code is deposited on Zenodo, which is a concrete reproducibility asset. The three-dimensional result, including the agreement between model H and model H0 when plotted against the renormalized viscosity, is a genuine measurement rather than an input. However, the two-dimensional conclusion is substantially weaker: it rests entirely on model H0, a universality assumption that is not verified in d=2, and on a two-volume, one-mode finite-size scaling analysis. Because the 2D claim appears in the abstract as a headline result, this limitation is load-bearing.","major_comments":[{"comment":"The two-dimensional dynamic exponent is extracted exclusively from model H0, not from the full model H. The justification in §II.B that the π^T self-advection can be neglected relies on the disparity between shear damping (η/ρ)k^2 and order-parameter damping Γk^4; at z≈2 in d=2 that ratio behaves as k^{z-2}=k^0, so the two mode-coupling terms are marginal relative to each other and the truncation is not obviously innocuous. Since §V.D states that the renormalized viscosity was not studied in 2D and no full model H simulation is reported, the quoted z≈2 is presently a measurement of model H0 rather than of model H. The claim should be supported by a 2D model H run or explicitly qualified as a model H0 result with an assessment of the universality risk.","section":"§V.D, Fig. 11"},{"comment":"The quoted exponents are derived from a single pair of volumes (L=40 and L=48), the second non-trivial momentum mode, and the collapse window C(t)>0.15, with no documented scan over the fitting window, the momentum mode, or additional lattice sizes. The statistical errors (for example z=2.11±0.015 in 2D) therefore do not include finite-size scaling or analysis-parameter systematics, and the extrapolation from z=2.11 to the asymptotic value z≈2 is made without an estimate of those systematics. The authors should report stability of the extracted exponents under variation of the window, mode number, and volume range, or state a corresponding systematic uncertainty.","section":"§V.C–V.D, Figs. 8, 9, 11"},{"comment":"The manuscript itself notes a shift in the critical bare mass for model H, δm2_c=-0.03 in 3D, attributed to the advection step not exactly conserving the potential-energy part of the Hamiltonian. This means the discretized dynamics does not sample exactly the target Hamiltonian, and the shift must be determined empirically. The text does not state explicitly which critical mass value is used in the 2D dynamic scaling runs (the model A/B value m2_c=-3.8240±0.0003 or the shifted model H value m2_c≈-3.859), and the uncertainty in this empirical shift is not propagated into the extracted exponent. Please state the value used and estimate the resulting systematic uncertainty.","section":"§V.A"}],"minor_comments":[{"comment":"There is a typo: 'centererd derivative kinetic energy' should read 'centered derivative kinetic energy'.","section":"Footnote 80"},{"comment":"References [78] and [79] list the same Zenodo DOI (10.5281/zenodo.14706997); if the 2D and 3D data sets have separate DOIs, the second DOI should be corrected.","section":"References [78,79]"},{"comment":"The caption states that self-advection results are offset horizontally for readability, but the offset is not quantified; adding the offset value to the caption would make the figure self-contained.","section":"Fig. 7"},{"comment":"The Kawasaki comparison uses only the broad estimate ξ∈[L/(2π),L/2] for the finite-volume correlation length; this yields a wide band and makes the comparison semi-quantitative. A direct measurement of the finite-volume correlation length from the static correlator would strengthen the crossover test.","section":"Eq. (67) and Fig. 9"}],"recommendation":"major_revision","confidential_remarks":"The 3D part of the paper is the strongest contribution and appears to support the model H dynamic-scaling claim. I would advise the editor that the 2D claim needs either additional simulation evidence from the full model H or a much more cautious presentation, since the current abstract-level claim goes beyond what the reported data establish. The paper's fit to the journal scope is good, and the availability of the code is a positive point."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the detailed companion to the model H PRL, and most of the value is in the method rather than a new exponent. The 3D result z = 3.013 ± 0.058 at ηR = 10^-2 is consistent with the epsilon expansion and with the companion PRL; that is a genuine confirmation. What is actually new in this manuscript is the lattice implementation of the Metropolis update with skew-symmetric advection, the viscosity renormalization study (Fig. 7), and the first 2D dynamic scaling result.\n\nThe algorithm deserves credit. The skew-symmetric discretization preserves energy conservation for the advection steps, the checkerboard Metropolis update is local and conserving, and the code and data are deposited. Static tests pass, though they are partly self-consistency checks because the Metropolis step enforces the equilibrium distribution. The dynamic exponent extraction is a measurement, not an input, and the crossover from z = 4 to z = 3 as a function of ηR is well presented and semi-quantitatively matches Kawasaki.\n\nNow the soft spots. The 2D claim z ≈ 2 is extracted exclusively from model H0, the truncation that drops the π^T self-advection. The paper justifies this by the relative damping rates, but in d = 2 that argument is marginal: at z = 2 the ratio (Γk^4)/( (η/ρ)k^2 ) behaves as k^0, so the two mode-couplings are marginal relative to each other at the fixed point. The stress-test note has this right. The authors do not simulate full model H in 2D, and the existing 2D viscosity literature already has log divergences. So the reported 2D exponent may be the exponent of a different theory. At minimum the claim should be softened to “model H0 in 2D” and the truncation tested, for example by running full model H at one small η.\n\nSecond, the finite-size extraction is thin. Two volumes, one momentum mode, and a window C(t) > 0.15. The quoted statistical error, z = 2.11 ± 0.015, is suspiciously small given that the two-volume collapse has no systematic scan over modes or fitting windows. Missing error bars in several key figures make it hard to judge the scatter. This is more a presentation and robustness issue than a fatal flaw.\n\nThird, trivial but embarrassing: refs [78] and [79] give the same Zenodo DOI. That should be corrected.\n\nWho is this for? People building stochastic hydrodynamic simulations for the QCD critical endpoint, and anyone interested in nonperturbative methods for critical dynamics. It deserves serious refereeing; the 3D part is solid, and the 2D result is worth publishing after the truncation concern is addressed head-on. I would accept for review with the expectation of a major revision on the 2D section.","headline":"Solid 3D model H result, but the new 2D claim rests on an untested truncation and a thin finite-size extraction.","tokens_in":23337,"tokens_out":3309,"would_cite":true,"duration_ms":27540,"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 Metropolis lattice algorithm simulates stochastic fluid dynamics near a critical point and measures the dynamic critical exponent $z\\simeq3$ in three dimensions and $z\\simeq2$ in two dimensions.","keywords":["stochastic fluid dynamics","model H","dynamic critical exponent","Metropolis algorithm","critical point","Ising universality class","QCD critical endpoint","finite-size scaling"],"falsifier":"Run full model H, with the self-advection term included, in two dimensions at bare viscosity $\\eta=10^{-2}$ on $L=40$ and $L=48$ lattices, and check whether the order-parameter correlation functions collapse with $z\\approx2$; if a clearly different exponent is required, the reported two-dimensional asymptotic value is not the exponent of the full model H.","tokens_in":22220,"feed_emoji":"🌀","tokens_out":10722,"duration_ms":96835,"temperature":0.7,"pith_summary":"This paper introduces a Metropolis accept-reject update for stochastic fluid dynamics in the Ising universality class, the theory known as model H, and uses it to measure the dynamic critical exponent $z$. Because the same step produces both diffusion and thermal noise, fluctuation-dissipation balance is automatic and the equilibrium distribution is governed solely by the microscopic free energy. In three dimensions the authors report $z\\simeq3$, with a best value $z_{\\rm eff}=3.013\\pm0.058$, and in two dimensions they find $z\\simeq2$, with a measured value $2.11\\pm0.015$ at the smallest viscosity simulated. They also observe a finite-size crossover from the mean-field value $z=4$ to the critical value, controlled by the correlation length and the renormalized shear viscosity. This matters because model H is expected to describe the non-equilibrium dynamics of the quark-gluon plasma near a possible critical endpoint, where existing analytic methods are low-order and hard to extend to expanding flows.","feed_headline":"Metropolis simulation gives z≈3 in 3D, z≈2 in 2D critical fluids","feed_subtitle":"A stochastic fluid-dynamics algorithm confirms model-H critical scaling and tracks the z=4-to-critical crossover.","key_machinery":"The object that carries the argument is a local, conserving Metropolis update: a random amount of the conserved field, either $\\phi$ or the momentum density, is moved between neighboring lattice sites and the move is accepted with probability $\\min(1,e^{-\\Delta H/T})$. The first moment of this step reproduces the dissipative current and the second moment supplies the noise, so fluctuation-dissipation balance is automatic and equilibrium is governed by the Hamiltonian alone. The advective mode-coupling terms are integrated with a strong-stability-preserving Runge-Kutta scheme and a skew-symmetric spatial discretization that conserves kinetic energy, and after each sweep the momentum field is projected onto transverse modes in Fourier space. The diagnostic that fixes $z$ is the dynamic scaling collapse of the order-parameter correlation function, $C_\\phi(t,k;L)=\\tilde C_\\phi(t/L^z,kL)$, compared between two volumes.","core_discovery":"The central claim is that the Metropolis algorithm correctly simulates model H and that the dynamic critical exponent takes its standard values: $z\\simeq3$ in three dimensions and $z\\simeq2$ in two dimensions. The best three-dimensional estimate, $z_{\\rm eff}=3.013\\pm0.058$ at renormalized viscosity $\\eta_R=10^{-2}$, is obtained by demanding data collapse of the order-parameter correlation function $C_\\phi(t,k;L)$ under the rescaling $t\\to(40/48)^z t$ for lattices with $L=40$ and $L=48$. In two dimensions the measured value at bare viscosity $\\eta=10^{-2}$ is $z=2.11\\pm0.015$, and the paper argues that the asymptotic value is $z\\simeq2$. The two-dimensional extraction is performed in model H0, the truncation of model H without the self-advection of the transverse momentum, which the authors expect to lie in the same dynamical universality class as model H. A crossover from the mean-field value $z=4$ to the critical value is observed and is semi-quantitatively described by the mode-coupling relaxation-rate formula used in the paper, with the crossover governed by the correlation length and the renormalized shear viscosity.","pith_inferences":["If the model-H0 and model-H universality assumption holds, the two-dimensional result is a first direct numerical value for the $d=2$ critical fluid; simulating the full two-dimensional model H at the same parameters would test this identification.","The paper's own estimate that the critical enhancement of $\\eta_R$ between $L=40$ and $L=48$ is only about one percent suggests that isolating the universal part of the shear viscosity will require substantially larger volumes or a different estimator.","The Metropolis framework is naturally suited to non-equilibrium and expanding flows, so the main obstacle to QCD-relevant simulations is likely to be the renormalization of the equation of state rather than the noise implementation."],"forward_implications":["The measured $z\\simeq3$ in $d=3$ agrees with two-loop epsilon-expansion values and provides direct numerical confirmation of model-H dynamic scaling.","The observed crossover from $z=4$ to the critical exponent means finite systems will generically show effective exponents between these limits, so size and viscosity must be controlled before comparing to experiments.","Because equilibrium is independent of the transport coefficients, the same Metropolis update can be reused at different viscosities and densities without re-tuning the noise.","The two-dimensional measurement gives $z\\simeq2.11\\pm0.015$ at the smallest viscosity simulated, with asymptotic $z\\simeq2$, supplying a numerical reference in a dimension where analytic methods are least controlled."],"supporting_citations":[{"why":"Defines model H and the dynamic-scaling framework that the paper's simulations test.","marker":"[12]"},{"why":"Companion paper reporting the first numerical model-H relaxation rate and exponent with this Metropolis method, which this work details and extends.","marker":"[49]"},{"why":"Earlier conserving Metropolis algorithm for model B from which the momentum-density update is adapted.","marker":"[47]"},{"why":"Introduces the Metropolis approach to stochastic hydrodynamics on which the algorithm is built.","marker":"[44]"},{"why":"Mathematical literature establishing the Metropolis method for stochastic partial differential equations.","marker":"[50]"},{"why":"Provides the mode-coupling relaxation-rate formula used to predict the crossover from $z=4$ to the critical value.","marker":"[67]"},{"why":"Supplies the two-loop epsilon-expansion values $z\\simeq3.0712$ in $d=3$ and $z=2.179$ in $d=2$ used as benchmarks.","marker":"[68]"},{"why":"Supplies the skew-symmetric discretization that makes the advection terms conserve kinetic energy on the lattice.","marker":"[62]"},{"why":"Provides the strong-stability-preserving Runge-Kutta scheme used to integrate the advective terms.","marker":"[63]"},{"why":"Gives the critical Binder-cumulant value used to locate the critical point in three dimensions.","marker":"[70]"}],"fun_headline_variants":["Metropolis fluid sim yields z~3 (3D), z~2 (2D)","Critical fluid dynamics: Metropolis method extracts z~3 and z~2","Metropolis model H simulation nails critical exponents z~3, z~2","Simulating critical fluid dynamics: z~3 (3D), z~2 (2D)"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The two-dimensional result is extracted from model H0, a truncation that drops the self-advection of the transverse momentum, and the analysis assumes both that model H0 lies in the same dynamical universality class as model H and that two volumes and one momentum mode are enough to determine the infinite-volume exponent.","fun_headline_variants_meta":{"raw":{"variants":["Metropolis fluid sim yields z~3 (3D), z~2 (2D)","Critical fluid dynamics: Metropolis method extracts z~3 and z~2","Metropolis model H simulation nails critical exponents z~3, z~2","Simulating critical fluid dynamics: z~3 (3D), z~2 (2D)"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00076,"raw_usage":{"total_tokens":3407,"prompt_tokens":1009,"completion_tokens":2398,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":625,"completion_tokens_details":{"reasoning_tokens":2308}},"tokens_in":625,"tokens_out":2398,"duration_ms":16320,"temperature":1.0,"reasoning_tokens":2308,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:39:50.363050+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run full model H, with the self-advection term included, in two dimensions at bare viscosity $\\eta=10^{-2}$ on $L=40$ and $L=48$ lattices, and check whether the order-parameter correlation functions collapse with $z\\approx2$; if a clearly different exponent is required, the reported two-dimensional asymptotic value is not the exponent of the full model H.","supporting_citations":[{"cited_title":"Vasil’ev, The Field Theoretic Renormalization Group in Critical Behavior Theory and Stochastic Dynamics (Chapman & Hall/CRC, 2004)","cited_arxiv_id":null,"evidence_quote":"Provides the strong-stability-preserving Runge-Kutta scheme used to integrate the advective terms."},{"cited_title":"Morinishi, T","cited_arxiv_id":null,"evidence_quote":"Provides the mode-coupling relaxation-rate formula used to predict the crossover from $z=4$ to the critical value."},{"cited_title":"Shu and S","cited_arxiv_id":null,"evidence_quote":"Supplies the two-loop epsilon-expansion values $z\\simeq3.0712$ in $d=3$ and $z=2.179$ in $d=2$ used as benchmarks."},{"cited_title":"Onuki, Phase Transition Dynamics(Cambridge University Press, 2002)","cited_arxiv_id":null,"evidence_quote":"Supplies the skew-symmetric discretization that makes the advection terms conserve kinetic energy on the lattice."}],"review_version":1}