REVIEW 1 major objections 3 minor 1 cited by
Hydrodynamic limit for some gradient and attractive spin models
T0 review · 1 major / 3 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Under diffusive scaling, three unbounded spin models converge to the heat equation, with model-dependent diffusion coefficients.
desk verdict A genuine advance on hydrodynamics for unbounded spin models, but the absolute-continuity proof has a wrong LLN constant that must be fixed before this is cited. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the combination of the gradient identity and attractiveness inside the entropy method, a standard route that proves hydrodynamic limits by tightness plus characterization of limit points. The gradient identity $W^{\sigma}_{x,x+1}=D_\sigma(\eta_{x+1}-\eta_x)$ turns the generator's action on occupation variables into $N^2 L_N^\sigma \eta_x = D_\sigma \Delta_N \eta_x$, which produces the diffusion term in the limiting equation. Attractiveness—preservation of the pointwise partial order between coupled configurations—is proved for generalized KMP by a basic Beta coupling and for discrete KMP and Harmonic models by checking the rate inequalities of Theorem 2.9 of [10]; it lets every expectation under the unknown initial measure be bounded by an expectation under the invariant measure $\nu^{\sigma}_{\hat\rho}$. A model-independent bound on the carré du champ, $D_\sigma(\eta_x-\eta_{x+1})^2 - L^{\sigma}_{x,x+1}(\eta_x\eta_{x+1}) \le D_\sigma(\eta_x^2+\eta_{x+1}^2)$, is the tightness input that controls the martingale oscillations uniformly in $N$.
What would settle it
A direct numerical check would evolve each model on a large torus from an initial profile satisfying the domination assumption and compare $N^{-1}\sum_x \eta^\sigma_x(N^2t)\delta_{x/N}$ with the heat-kernel solution with diffusivity $D_\sigma$; a mismatch that persists as $N$ grows would contradict Theorem 2.8. For the harmonic spin restriction, an algebraic search for ordered configurations with $s<1/2$ violating inequalities (3.1)-(3.2) would either confirm that the restriction is necessary or suggest it can be removed.
Extended reading notes
Core claim
The central discovery is Theorem 2.8: for each of the three models, if the initial measures are associated to a bounded density profile $\rho_0$ and are stochastically dominated by the invariant product measure $\nu^{\sigma}_{\hat\rho}$ for some $\hat\rho>0$ (and if the harmonic spin satisfies $s\ge 1/2$), then the empirical measure $N^{-1}\sum_x \eta_x^\sigma(N^2 t)\delta_{x/N}$ converges in probability to the unique weak solution of $\partial_t \rho = D_\sigma \Delta \rho$, with $D_{\mathrm{gKMP}}=D_{\mathrm{dKMP}}=1/2$ and $D_{\mathrm{Harm}}=1/(2s)$. The proof's contribution is to verify the two structural properties that make the entropy method run: the gradient identity, which writes the instantaneous current as a discrete gradient with constant $D_\sigma$, and attractiveness, which gives monotonicity under coupling. The harmonic models require $s\ge 1/2$ because that is where the attractiveness proof works.
Load-bearing premise
The load-bearing premise is that the initial random configuration is statistically no larger, site by site, than an equilibrium configuration at some positive density; all moment estimates for the unbounded occupation variables use this comparison, so without it the proof has no control.
Editorial extensions
If this is right
- At macroscopic scale the evolution is deterministic diffusion: the random empirical measure converges in probability to the heat-kernel solution, so fluctuations vanish in the $N\to\infty$ limit.
- The diffusion coefficient is density-independent, so the hydrodynamic equation is linear even though the microscopic jump rates are nonlinear functions of the occupation variables.
- The result extends without conceptual change to any dimension and to the infinite lattice, where the same heat equation with the same $D_\sigma$ is predicted.
- Because the weak solution of the heat equation is unique, the whole sequence of empirical measures converges—subsequence extraction is only an intermediate step.
Reading between the lines
- The restriction to dominated initial measures is probably not essential: approximating arbitrary bounded profiles by dominated ones and using the linearity of the heat equation should give the same limit; if true, the theorem would cover all physically natural initial states.
- The three models share constant diffusivity and quadratic mobility, so the next-order fluctuations around the heat profile should be Gaussian with a covariance fixed by the mobility; this is a concrete prediction that could be checked by computing the associated non-equilibrium fluctuation field.
- The attractiveness criterion used for the harmonic family suggests an easy numerical test for $s<1/2$: if ordered configurations can be found whose jump rates violate inequalities (3.1)-(3.2), that pinpoints why the spin restriction appears; if no violation exists, the restriction may be removable.
- The continuous harmonic models mentioned in the paper have formally identical transport coefficients; once their well-posedness is settled, the same heat equation should emerge.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proves hydrodynamic limits for three one-dimensional conservative spin models with unbounded occupation variables: the generalized KMP model, the discrete KMP model, and the harmonic models. Under diffusive scaling and for initial laws satisfying Assumption 2.6 (association to a profile plus stochastic domination by an invariant measure), Theorem 2.8 states that the empirical measure converges in probability to the unique weak solution of the heat equation, with diffusion coefficients D=1/2 for the two KMP-type models and D=1/(2s) for the harmonic models. The strategy combines the gradient identity (2.13) with a proof of attractiveness (Section 3) and the entropy method; tightness uses a common carré du champ bound (Lemma 4.1), and the limit points are characterized via Dynkin martingales. The spin restriction s≥1/2 for the harmonic models and the initial-measure restriction are stated explicitly by the authors.
Significance. The paper is a solid contribution to the hydrodynamic limit theory for unbounded spin systems. Its strengths are explicit: the diffusion coefficients are computed directly from the generators in (2.12), the gradient identities are given in detail, and the attractiveness of all three models is proved in full, including a novel verification for the harmonic family via the Gobron-Saada criterion. The common carré du champ estimate (4.6) reduces the tightness proof to a uniform second-moment bound and is transparently adapted to each model. If the identified gap in Section 4.2.2 is repaired, the main theorem provides a unified treatment of three models that previous methods handled only under stronger assumptions. The authors are also honest about the limitations imposed by Assumption 2.6 and by the harmonic spin range.
major comments (1)
- [Section 4.2.2, Eq. (4.8)] The inequality (4.8) is asserted with the wrong constant. For the generalized KMP model, (2.3) with m=1 gives E_{ν_ρ}[η_0]=2sρ, and for the harmonic model, (2.9) with m=1 gives the same value 2sρ; the discrete KMP model has mean ρ. Therefore, under the dominating measure ν_{ρhat}, the law of large numbers gives <π_t,|G|> → 2s ρhat ||G||_{L^1} for gKMP and Harm, not ρhat ||G||_{L^1}. For s>1/2, which is allowed for both gKMP and Harm, the event used in (4.8) has probability tending to 1 rather than 0; for example μ_N=ν_{ρhat} satisfies Assumption 2.6 and produces a limit violating (4.8). This is a genuine error in a load-bearing step, because (4.8) is the only input that yields absolute continuity of the limit measures. The repair is local: replace ρhat everywhere in (4.8) and in the subsequent supremum argument by C=(2s∨1)ρhat. Since Proposition B.1 only requires a finite constant, the rest of the characterization, the martingale estimates, and the uniqueness argument are unaffected.
minor comments (3)
- [Theorem 2.8] The set {gKMP, dKPM, Harm} should read {gKMP, dKMP, Harm}; the same typo appears in the paragraph immediately before Theorem 2.8.
- [Appendix A.1] The displayed equality for D_gKMP(η_x−η_{x+1})^2−L(η_xη_{x+1}) omits the cross term −2Ips η_x η_{x+1}; the subsequent inequality remains valid because this omitted term is negative, but the equality as written is false.
- [Section 4.1 and Proposition 4.2] The notation ||ΔG||_∞ in Proposition 4.2 refers to the discrete Laplacian Δ_N G defined there, but the same symbol Δ is also used for the continuum Laplacian in Definition 2.7; a footnote or a repeated definition would remove the ambiguity.
Circularity Check
No circularity: diffusion coefficients are computed from the generators, attractiveness is proved directly, and weak-solution uniqueness is standard.
full rationale
The derivation chain of Theorem 2.8 is self-contained and does not reduce to its inputs. The diffusion coefficients in (2.12) are obtained by explicit generator computations (2.11)-(2.13) and Appendix A; they are not fitted to the heat equation. The characterization of limit points uses Dynkin martingales, tightness, the carré du champ bound, and the standard uniqueness of weak solutions cited from an external textbook [12, Theorem A.2.4.4]. No step assumes the heat equation in order to prove it. Attractiveness is established in Section 3 either directly via the basic coupling (Theorem 3.2) or by verifying the externally stated Gobron-Saada criterion (Proposition 3.3); the criterion is not a self-citation. Assumption 2.6 is an explicitly stated restriction on initial measures, and the paper clearly marks this limitation in Sections 1.3 and 1.6; it is not a hidden use of the target PDE. The few self-citations that occur, such as [6] in Remark 2.2 and previous related work in the introduction, are for model provenance and context, not load-bearing for Theorem 2.8. The separate technical concern about the law-of-large-numbers constant in Section 4.2.2 is a correctness issue, not a circularity issue: for gKMP and Harmonic models (2.3) and (2.9) give mean 2sρ rather than ρ under ν_ρ, but the absolute-continuity argument only needs a finite constant, so the claimed theorem is not forced by a circular fit or by definition.
Assumptions & free parameters
assumptions (5)
- standard math The entropy method of Guo-Papanicolaou-Varadhan and Kipnis-Landim is applicable and sufficient for the hydrodynamic limit.
- standard math The Gobron-Saada attractiveness criterion (Theorem 2.9 of [10]) correctly characterizes attractiveness for the particle models.
- ad hoc to paper Initial measures are stochastically dominated by some invariant measure (Assumption 2.6).
- ad hoc to paper The spin parameter satisfies s >= 1/2 for the Harmonic models.
- domain assumption The single-site invariant measures are Gamma, geometric, and negative binomial, with the stated moments.
Cite this review
Pith. "Pith review of Hydrodynamic limit for some gradient and attractive spin models." pith.science (2026). https://pith.science/paper/EZNV6P3W
@misc{pith2026250511092,
author = {Pith},
title = {Pith review of: Hydrodynamic limit for some gradient and attractive spin models},
year = {2026},
howpublished = {\url{https://pith.science/paper/EZNV6P3W}},
note = {Machine review of arXiv:2505.11092}
}
read the original abstract
We study the hydrodynamic limit for three gradient spin models: generalized Kipnis-Marchioro-Presutti (KMP), its discrete version and a family of harmonic models, under symmetric and nearest-neighbor interactions. These three models share some universal properties: occupation variables are unbounded, all these processes are of gradient type, their invariant measures are product with spatially homogeneous weights, and, notably, they are all attractive, meaning that the process preserves the partial order of measures along the dynamics. In view of hydrodynamics of large-scale interacting systems, dealing with processes taking values in unbounded configuration spaces is known to be a challenging problem. In the present paper, we show the hydrodynamic limit for all three models listed above in a comprehensive way, and show as a main result, that, under the diffusive time scaling, the hydrodynamic equation is given by the heat equation with model-dependent diffusion coefficient. Our novelty is showing the attractiveness for each model, which is crucial for the proof of hydrodynamics.
Forward citations
Cited by 1 Pith paper
-
Convergence of the KMP model to the KPZ equation
The KMP heat transport process converges, in a t^{3/4} scaling window, to the multiplicative-noise stochastic heat equation (the exponential of KPZ) with noise coefficient 1/(2√α).
Reference graph
Works this paper leans on
-
[1]
Billingsley.Convergence of probability measures
P. Billingsley.Convergence of probability measures. John Wiley & Sons, 2013
2013
-
[2]
M. Capanna, D. Gabrielli, and D. Tsagkarogiannis. On a class of solvable stationary non equilibrium states for mass exchange models.Journal of Statistical Physics, 191(2):25, 2024
work page 2024
-
[3]
G. Carinci, C. Giardinà, C. Giberti, and F. Redig. Duality for stochastic models of transport. Journal of Statistical Physics, 152:657–697, 2013
work page 2013
-
[4]
L. Fajfrová, T. Gobron, and E. Saada. Invariant measures of mass migration processes. Electronic Journal of Probability, 21(60):1–52, 2016
work page 2016
-
[5]
S. Feng, I. Iscoe, and T. Seppäläinen. A microscopic mechanism for the porous medium equation. Stochastic Processes and their Applications, 66(2):147–182, 1997
work page 1997
-
[6]
C. Franceschini, R. Frassek, and C. Giardinà. Integrable heat conduction model.Journal of Mathematical Physics, 64(4), 2023
work page 2023
-
[7]
C. Franceschini, P. Gonçalves, and F. Sau. Symmetric inclusion process with slow bound- ary: Hydrodynamics and hydrostatics.Bernoulli, 28(2):1340–1381, 2022
work page 2022
-
[8]
R.Frassek, C.Giardinà, andJ.Kurchan. Non-compactquantumspinchainsasintegrable stochastic particle processes.Journal of Statistical Physics, 180(1):135–171, 2020
work page 2020
Show all 17 references
-
[9]
Giardinà and F
C. Giardinà and F. Redig.Duality for Markov processes: a Lie-algebraic approach. In preparation, 2025
2025
-
[10]
Gobron and E
T. Gobron and E. Saada. Couplings, attractiveness and hydrodynamics for conservative particle systems. Annales de l’IHP Probabilités et statistiques, 46(4):1132–1177, 2010
2010
-
[11]
M. Z. Guo, G. C. Papanicolaou, and S. S. Varadhan. Nonlinear diffusion limit for a system with nearest neighbor interactions.Communications in Mathematical Physics, 118(1):31–59, 1988
1988
-
[12]
Kipnis and C
C. Kipnis and C. Landim. Scaling limits of interacting particle systems, volume 320. Springer Science & Business Media, 1999
1999
-
[13]
Kipnis, C
C. Kipnis, C. Marchioro, and E. Presutti. Heat flow in an exactly solvable model.Journal of Statistical Physics, 27:65–74, 1982
1982
-
[14]
T. M. Liggett.Interacting particle systems. Springer, 1985
1985
-
[15]
L. Lipatov. High energy asymptotics of multi–colour QCD and exactly solvable lattice models. Journal of Experimental and Theoretical Physics Letters, 59:596–599, 1994
1994
-
[16]
Rezakhanlou
F. Rezakhanlou. Hydrodynamic limit for attractive particle systems onZd. Communi- cations in mathematical physics, 140(3):417–448, 1991
1991
-
[17]
Hydrodynamiclimitforaspinsystemonamultidimensional lattice
Y.SuzukiandK.Uchiyama. Hydrodynamiclimitforaspinsystemonamultidimensional lattice. Probability theory and related fields, 95:47–74, 1993. 24
1993
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.