Pith. sign in

REVIEW 4 major objections 5 minor 12 references

This paper reports the first numerical simulation of model H and measures the dynamic critical exponent z ≈ 3.01 at the critical point.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

First simulations of model H stochastic fluid dynamics near a critical point give a dynamic critical exponent z ≈ 3.01 and viscosity renormalization.

T0 review reviewed 2026-08-05 challenge →

load-bearing objection A clear, honest summary of the authors' own prior model H simulations, but the quantitative evidence lives elsewhere—treat it as a briefing, not a self-contained proof. the 4 major comments →

arxiv 2509.00545 v1 pith:4TPNA7LG submitted 2025-08-30 nucl-th

Simulating stochastic fluid dynamics

Chandrodoy Chattopadhyay , Josh Ott , Thomas Schaefer , Vladimir Skokov This is my paper
classification nucl-th
keywords model Hstochastic fluid dynamicsdynamic critical exponentQCD critical pointcritical slowing downMetropolis algorithmviscosity renormalizationfinite-size scaling
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Model H is the stochastic hydrodynamic theory expected to govern long-wavelength fluctuations near a liquid-gas-type critical endpoint, including a possible critical point in the QCD phase diagram. Until now it had never been simulated numerically, because fluctuations on all scales make renormalization and fluctuation-dissipation balance hard to preserve on a lattice. This paper presents a discretization scheme built on a conserving Metropolis update with transverse projection and skew-symmetrized derivatives, and simulates the real-time dynamics of the order parameter and fluid momentum. At the critical point the measured dynamic critical exponent is z ≈ 3.01, clearly different from mean-field behavior and consistent with epsilon-expansion predictions. The simulations also reproduce the predicted renormalization of shear viscosity, the 'stickiness of sound,' and show a finite-size crossover from z ≈ 4 to z ≈ 3.

Core claim

Model H—stochastic hydrodynamics for an order parameter coupled to transverse momentum—has not previously been simulated, and this paper reports simulations of it in two and three dimensions. The order parameter is the specific entropy, the momentum is projected transverse, and the dynamics are generated by a conserving Metropolis update whose first and second moments reproduce the diffusion and noise terms of the continuum equations. Tuning the mass parameter to the critical point reproduces the static 3d Ising exponents; the finite-size scaling of the order-parameter correlation time then gives a dynamical critical exponent z ≈ 3.01, consistent with epsilon-expansion predictions and clearl

What carries the argument

The conserving Metropolis update is the central mechanism: a trial update moves a charge between neighboring lattice sites, accepted with probability min(1, e^{-ΔH/T}). Its first moment realizes the diffusion equation and its second moment reproduces the noise correlation, so fluctuation-dissipation balance is built in. Transverse projection removes the longitudinal component of momentum, and skew-symmetrized derivatives keep the advection step conservative despite large sub-lattice fluctuations.

Load-bearing premise

The lattice discretization is assumed to faithfully represent the continuum model H equations, preserving conservation laws and fluctuation-dissipation balance; the paper does not show a convergence or conservation-law validation.

What would settle it

Run the finite-size scaling measurement at several different lattice spacings: if the extracted z shifts systematically away from ≈3.01, or if global momentum conservation is violated as resolution increases, the central claim fails.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • If z ≈ 3.01 is correct, the relaxation time near a QCD critical point grows as ξ^3 rather than ξ^2, changing predicted fluctuation observables in heavy-ion collisions.
  • The algorithm lifts a longstanding barrier: model H can now be studied quantitatively in real time, not just through analytic approximations.
  • The observed finite-size crossover from model-B-like z ≈ 4 to model-H z ≈ 3 implies that in a finite fireball the effective exponent depends on viscosity and system size.
  • Reproducing the 'stickiness of sound' viscosity renormalization means the lattice scheme captures non-perturbative momentum transport, not just order-parameter diffusion.
  • Embedding the framework into an expanding relativistic fluid is the stated next step; the current results provide the benchmark dynamics to embed.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Editorial extension: the Metropolis-plus-projection construction should generalize to other stochastic hydrodynamic theories with different conserved charges, but the paper does not claim this.
  • Editorial extension: because the paper flags possible violations of ideal conservation laws and does not show a lattice-spacing convergence study, the continuum limit of z ≈ 3.01 remains an open check.
  • Editorial extension: once the advertised mapping to QCD parameters is completed, z ≈ 3 implies critical slowing-down stronger than mean-field estimates in a heavy-ion context; this translation is not made in the paper.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper reports simulations of real-time stochastic fluid dynamics in the universality class of model H, motivated by QCD critical-point phenomenology. The authors describe a numerical scheme combining conserving Metropolis updates for the dissipative and noise parts of the equations of motion with a skew-symmetrized derivative discretization for the advective terms. They report two main results: (i) a UV-sensitive renormalization of the shear viscosity, dominated by self-advection and consistent with the 'stickiness of sound' effect, and (ii) a dynamic critical exponent z ≃ 3.01 extracted from finite-volume scaling directly at the critical point, which they claim is clearly different from mean-field behavior and consistent with the epsilon expansion. They also observe a crossover from z ≃ 4 to z ≃ 3 as the physical viscosity is reduced and evidence for universality between model H and a truncated 'model H0'.

Significance. If the numerical method can indeed simulate the stochastic equations of model H and reproduce the expected universal dynamic scaling, this is a substantial step toward first-principles modeling of critical fluctuations in heavy-ion collisions and other fluid systems near a critical point. The paper explicitly identifies a challenging combination of issues—large thermal fluctuations on all scales, renormalization of UV divergences, and the need to preserve conservation laws and fluctuation-dissipation relations—and the proposed conserving Metropolis update is an interesting and potentially powerful idea. However, the quantitative claims made here (z ≈ 3.01, viscosity renormalization, crossover, Ising static exponents) are not actually demonstrated within this manuscript; the text repeatedly refers to the authors' previous papers [4,5] for the supporting analysis. As a standalone contribution, the paper therefore lacks the evidence needed to substantiate its central numerical results, even though those results may be correct.

major comments (4)
  1. [§4, Fig. 2] The central quantitative claim is 'We observe a value z≃3.01 which is clearly different from mean-field behavior', but no error bars, scaling-fit details, system sizes, or lattice-spacing dependence are provided. The right panel of Fig. 2 shows a curve labeled 'model H' but with no uncertainty and no description of how z_eff is extracted from the order-parameter correlation function. Since the entire paper's significance rests on this exponent, the result is under-determined in this manuscript. The pointer to Refs. [4,5] is not a substitute, because the present paper makes a standalone claim; please either include the scaling analysis or explicitly state that this is a summary and rephrase the claim as a report of previous work.
  2. [§3] The discretization of the advective terms is the step on which the model-H universality claim depends, but the manuscript gives no validation that the skew-symmetrized derivative scheme combined with the Metropolis update and transverse projection preserves conservation laws and fluctuation-dissipation balance at finite lattice spacing. The text itself acknowledges 'we may encounter large violations of conservation laws that follow from the ideal equations of motion in the continuum limit'. No convergence study with lattice spacing, no comparison between different advection discretizations, and no conservation-law residual plots are shown. Without such tests, the measured relaxation times and the extracted z could be biased by lattice-scale violations, so the identification of the simulated model with continuum model H is not established.
  3. [§4] The statement 'We have verified that m^2 can be tuned to a critical point at which the simulation reproduces the static critical exponents of the 3d Ising Model' is made without any supporting data, method, or reference to which figure shows this verification. The static critical point provides the physical condition under which the dynamic exponent is measured; if the static tuning is not shown, the dynamic scaling claim lacks its foundation. Please provide the static-exponent determination or specify that it is taken from Ref. [4] and give the relevant quantities.
  4. [§4, Fig. 2] The crossover from z_eff ≈ 4 to z_eff ≈ 3 as a function of physical viscosity is described only qualitatively. The right panel of Fig. 2 is unclear: the definition of z_eff, the observable from which it is extracted, the system size, and the number of independent runs are not given. The curve labeled 'Kawasaki approximation' is not defined or referenced. Because this crossover is used to argue for universality, it needs a precise definition and an uncertainty estimate, or it should be explicitly deferred to a forthcoming/previous detailed paper.
minor comments (5)
  1. [§2, Eq. (2)] The index structure of the last term in Eq. (2) is difficult to parse: the expression 'πT_j δH/δπT_k' is not explicitly summed over j and k in the displayed formula, and the overall index of the term is not clear. Please clarify the notation.
  2. [§3] The sentence 'We follow the same procedure for ⃗π. We perform a conserving Metropolis update for all components of the momentum density, and then apply a transverse projection operator.' leaves open the question of whether the transverse projection is part of the Markov chain or a separate step applied to the output; this affects the conservation properties. Please state the order of operations and the effect on the momentum components.
  3. [Fig. 2] 'Model H0' is used in the figure caption and the text but is never defined in the body. Please define it explicitly, stating which terms of Eq. (1)-(2) are dropped.
  4. [References] Reference [3] is a arXiv preprint without a journal or DOI; if possible, update to the published version or include the arXiv identifier consistently. Also, the reference to 'Kawasaki approximation' in Fig. 2 has no citation; either add the appropriate reference or remove the label.
  5. [§4] The sentence 'We have demonstrated that we can also measure more complicated observables, such as the relaxation rate of higher moments of the order parameter.' is not accompanied by a figure, equation, or data. Either present an example or delete the sentence, as it currently serves only as a placeholder.

Circularity Check

0 steps flagged

No significant circularity; central results are simulation outputs compared with external benchmarks.

full rationale

The paper is primarily a companion summary that points to the authors' earlier papers [4,5] for detailed numerical evidence. That self-citation is not a circular derivation: [4,5] are prior published reports containing the simulations, and the present text presents measured outcomes (z approximately 3.01, viscosity renormalization, static Ising exponents) that are compared with external expectations such as the epsilon expansion, the 'stickiness of sound' effect, and 3d Ising static critical exponents. The only parameter explicitly tuned is m^2, which is adjusted to the critical point; the resulting exponents and the dynamic critical exponent are outputs of the simulation, not inputs constrained to match the claimed conclusion. The discretization method (Metropolis updates plus skew-symmetrized derivatives) is adopted from the literature [8,9] and is a numerical implementation choice, not a device that defines z. Potential conservation-law violations are a validation/correctness concern, not evidence that a prediction was built from the result it claims. No equation or fitted parameter is shown to equal z or the renormalized viscosity by construction.

Axiom & Free-Parameter Ledger

1 free parameters · 4 axioms · 0 invented entities

The main assumptions are physical (model H universality, order parameter identification) and algorithmic (Metropolis discretization preserves conservation laws and fluctuation-dissipation). No new entities are introduced. The only fitted quantity is the bare mass squared used to tune to the critical point.

free parameters (1)
  • bare squared mass m^2 = not given
    Tuned in the simulations to locate the critical point where static exponents match the 3d Ising model, as stated in Section 4. This is a parameter adjusted to make the simulation sit at criticality.
axioms (4)
  • domain assumption Model H equations describe the universal critical dynamics near a QCD critical endpoint.
    Invoked in the Introduction with citations to Hohenberg-Halperin and Son-Stephanov; the relevance of the entire simulation program rests on this physical assumption, not proven in the paper.
  • domain assumption The order parameter phi can be identified with the specific entropy s/n.
    Used in Section 2 following Ref. [6]; this mapping is needed to connect the abstract model H field to QCD variables.
  • ad hoc to paper The conserving Metropolis update realizes the diffusion equation and noise correlations, and skew-symmetrized derivatives preserve conservation laws.
    Central to the numerical method in Section 3. The paper cites Refs. [8,9] for these properties but provides no derivation or validation within this text; the correctness of the simulation depends on this unproved algorithmic premise.
  • domain assumption Finite-volume scaling tau ~ L^z is a valid way to extract the dynamic critical exponent z.
    Used in Section 4 to measure z by varying the system size L; a standard, but unstated in detail, assumption of the scaling analysis.

reviewed 2026-08-05 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Simulating stochastic fluid dynamics." pith.science (2026). https://pith.science/paper/4TPNA7LG

@misc{pith2026250900545,
  author       = {Pith},
  title        = {Pith review of: Simulating stochastic fluid dynamics},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4TPNA7LG}},
  note         = {Machine review of arXiv:2509.00545}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

We present simulations of the real time dynamics of a fluid in the vicinity of a critical endpoint in the phase diagram. The relevant hydrodynamic theory is known as model H, and it is expected to describe the long-distance dynamics of QCD matter in the vicinity of a possible critical endpoint in the QCD phase diagram.

Figures

Figures reproduced from arXiv: 2509.00545 by Chandrodoy Chattopadhyay, Josh Ott, Thomas Schaefer, Vladimir Skokov.

Figure 1
Figure 1. Figure 1: Left panel: Order parameter ϕ (color coded) and fluid momentum ⃗π (arrows) configuration in a stochastic two-dimensional fluid. Right panel: Order parameter configuration in a three dimensional fluid. 2 Model H Model H describes the interaction of an order parameter density ϕ with the momentum density ⃗π of the fluid. This equations of motion are given by [1] ∂tϕ = Γ ∇ 2 [PITH_FULL_IMAGE:figures/full_fig_… view at source ↗
Figure 2
Figure 2. Figure 2: Left panel: Physical shear viscosity as a function of the bare viscosity in a non-critical fluid. Several truncations are shown: Model H0 (no self-advection), self-advection only, and pure diffusion (no interactions). Right panel: Dynamical critical exponent z extracted from a scaling analysis of the order parameter correlation function in a finite volume as a function of the physical viscosity. 3 Numerica… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

12 extracted references · 2 canonical work pages · 1 internal anchor

  1. [1]

    Hohenberg, B.I

    P.C. Hohenberg, B.I. Halperin, Theory of Dynamic Critical Phenomena , Rev. Mod. Phys. 49, 435 (1977). 10.1103/RevModPhys.49.435

    doi:10.1103/revmodphys.49.435 1977

  2. [2]

    Son, M.A

    D.T. Son, M.A. Stephanov, Dynamic universality class of the QCD critical point , Phys. Rev. D 70, 056001 (2004), hep-ph/0401052

    Pith/arXiv arXiv 2004

  3. [3]

    Basar, Recent developments in relativistic hydrodynamic fluctuations (2024), 2410.02866

    G. Basar, Recent developments in relativistic hydrodynamic fluctuations (2024), 2410.02866

    Pith/arXiv arXiv 2024

  4. [4]

    Chattopadhyay, J

    C. Chattopadhyay, J. Ott, T. Sch \"a fer, V.V. Skokov, Simulations of Stochastic Fluid Dynamics near a Critical Point in the Phase Diagram , Phys. Rev. Lett. 133, 032301 (2024), 2403.10608 . 10.1103/PhysRevLett.133.032301

    Pith/arXiv arXiv 2024

  5. [5]

    Chattopadhyay, J

    C. Chattopadhyay, J. Ott, T. Sch \"a fer, V.V. Skokov, Critical fluid dynamics in two and three dimensions , Phys. Rev. D 111, 034026 (2025), 2411.15994 . 10.1103/PhysRevD.111.034026

    Pith/arXiv arXiv 2025

  6. [6]

    Transits of the QCD Critical Point

    Y. Akamatsu, D. Teaney, F. Yan, Y. Yin, Transits of the QCD critical point , Phys. Rev. C 100, 044901 (2019), 1811.05081 . 10.1103/PhysRevC.100.044901

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  7. [7]

    Parotto, M

    P. Parotto, M. Bluhm, D. Mroczek, M. Nahrgang, J. Noronha-Hostler, K. Rajagopal, C. Ratti, T. Sch\"afer, M. Stephanov, QCD equation of state matched to lattice data and exhibiting a critical point singularity , Phys. Rev. C 101, 034901 (2020), 1805.05249 . 10.1103/PhysRevC.101.034901

    Pith/arXiv arXiv 2020

  8. [8]

    Florio, E

    A. Florio, E. Grossi, A. Soloviev, D. Teaney, Dynamics of the O(4) critical point in QCD , Phys. Rev. D 105, 054512 (2022), 2111.03640

    Pith/arXiv arXiv 2022

  9. [9]

    Morinishi, T

    Y. Morinishi, T. Lund, O. Vasilyev, P. Moin, Fully conservative higher order finite difference schemes for incompressible flow, Journal of Computational Physics 143, 90 (1998)

    work page 1998

  10. [10]

    Kovtun, G.D

    P. Kovtun, G.D. Moore, P. Romatschke, The stickiness of sound: An absolute lower limit on viscosity and the breakdown of second order relativistic hydrodynamics , Phys. Rev. D 84, 025006 (2011), 1104.1586

    Pith/arXiv arXiv 2011

  11. [11]

    , " * write output.state after.block = add.period write newline

    ENTRY address archive author booktitle chapter collaboration edition editor eid eprint howpublished institution isbn issn journal key month note number numpages organization pages publisher school series title type url doi volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION ...

  12. [12]

    write FUNCTION output.bibitem newline

    " write FUNCTION output.bibitem newline " " " " * write cite write " " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max su...

This paper was first reviewed by deepseek-v4-flash on August 5, 2026.