REVIEW 1 major objections 3 minor 1 cited by
Poisson statistics, vanishing correlations, and extremal particle limits for symmetric exclusion in d > 1
T0 review · 1 major / 3 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read For symmetric exclusion in d>1, extremes match independent-particle statistics when the expected count converges, with an explicit Gumbel limit for polynomial initial shapes.
desk verdict Solid extension of the d=1 Gumbel work to d>1 with detailed proofs; the abstract overstates the d=2,3 scope slightly but the polynomial-shape theorem holds. 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 argument runs through the stirring representation of SSEP, self-duality, and negative association. The load-bearing identity is the strong Rayleigh criterion (Lemma 3.8): for each t the sum of occupation variables in any finite region has the same law as a sum of independent Bernoulli variables, so N_t ⇒ Poisson(λ) iff E[N_t] → λ and E[N_t] − Var(N_t) → 0. The correlation term E[N_t] − Var(N_t) = S_t + C_t is then controlled by Lemma 2.1, which bounds C_t by an integral over a two-particle displacement kernel, and by Propositions 4.5–4.6, which reduce the bounds to expectations of the one-dimensional random walk ζ_{t/d}. The dimension enters through γ_d(t), the expected number of returns to 0 of a (d−1)-dimensional random walk: √t for d = 2, log t for d = 3, and 1 for d ≥ 4; this drives the sharpness of correlation decay. For polynomial profiles, Gaussian large-deviation asymptotics (Lemma 9.4) evaluate the mean E[N_t] and the rates of S_t and C_t.
What would settle it
In d = 2 or 3, construct a shape function g_i satisfying Conditions (A)–(C) with sup_t E[N_t] < ∞ and lim E[N_t] = λ > 0 but with E_0[G(ζ_{t/d} − z)1(ζ_{t/d} > z)] decaying slower than the rates in (4.2); then simulate or estimate S_t + C_t to see whether it stays bounded away from zero. If it does, the Poisson convergence in Corollary 4.3 fails for that profile, disproving the sufficiency of the mean-convergence condition without the extra hypothesis. A direct numerical check for polynomial shapes in d = 2 and 3, comparing the empirical distribution of N_t to Poisson(M $e^{{−βx}}$) at finite t and verifying the Gumbel tail, would also settle whether the stated rates are sharp.
Extended reading notes
Core claim
The paper establishes that for SSEP in d ≥ 2 with step initial conditions satisfying Conditions (A)–(C), the criterion for Poisson convergence of N_t is the same as for independent Bernoulli arrays: lim E[N_t] = λ and lim Var(N_t) = λ, equivalently E[N_t] − Var(N_t) → 0. Using the strong Rayleigh property, this reduces to showing that the quantities S_t (sum of squared single-particle occupation probabilities) and C_t (sum of negative covariances) tend to zero. The main theorem proves S_t + C_t → 0 whenever sup_t E[N_t] < ∞ in d ≥ 4, and in d = 2, 3 under condition (4.2), which bounds E_0[G(ζ_{t/d} − z)1(ζ_{t/d} > z)] by explicit o-rates. For polynomial shapes g_i(u) = c_i $u^{{α_i}}$ + r_i, the paper identifies z = b_t(x + a_t) with b_t = (βt/(d log t))^{1/2} and a_t = log(t / ((2π)^{1/(2β)} (log t)^{(β+1)/(2β)})), β = 1 + Σ α_i, and proves N_t ⇒ Poisson(M $e^{{−βx}}$) with the constant M in (5.4). Consequently, the rescaled maximum satisfies a Gumbel law and each order statistic $X_t^{{(m)}}$ has limit (5.3).
Load-bearing premise
In d = 2 and 3, the conclusion that correlations vanish and N_t becomes Poisson rests on the extra geometric condition (4.2), which the paper verifies for polynomial shapes but does not prove is necessary for all allowed profiles; if this condition fails for some subexponential shape with sup_t E[N_t] < ∞, the main theorem no longer applies in low dimensions.
Editorial extensions
If this is right
- If the central claim is right, then in d ≥ 4 the Poisson limit for the number of particles beyond a moving frontier follows directly from convergence of its expectation, for every shape satisfying the mild Conditions (A)–(C), with no extra hypothesis.
- For polynomial profiles, the explicit scaling gives the maximum position X_t with mean of order √((β/d) t log t) and variance of order (π²/(6β)) t / log t, so extremes move faster than the diffusive scale √t.
- The limiting constant M depends only on the leading-order coefficients c_i and the exponents α_i, not on the intercepts r_i, so adding or subtracting finitely many particles (or periodically sparse particles) does not change the extreme-value asymptotics.
- The order-statistics limit (5.3) gives the full joint asymptotic distribution of the top particles: for each fixed m, P(X_t^{(m)} ≤ b_t(x+a_t)) converges to the Poisson tail sum, meaning the top particles behave like the record points of an independent-particle system.
- In periodic product-measure initial conditions, the same extremes result holds with the mean scaled by the average density ρ̄, so the conclusion is robust to deterministic initial configurations beyond pure step profiles.
Reading between the lines
- The result suggests a universality principle for d > 1: the extreme-value statistics of the exclusion process coincide with those of non-interacting random walks whenever the expected frontier count stabilizes, so interaction effects are asymptotically invisible in the far tail even though they shape the bulk.
- The extra condition (4.2) in d = 2, 3 may well be necessary for subexponential but non-polynomial shapes, because the slower decorrelation in low dimensions could leave a residual correlation that no choice of level z can remove; a counterexample would settle the open question the paper raises.
- Since the paper's machinery reduces everything to one-dimensional random-walk functionals, one could numerically test convergence rates for finite t by simulating the random walk expectations E_0[G(ζ_{t/d} − z)1(ζ_{t/d} > z)] and comparing S_t + C_t bounds to the predicted Poisson behavior, even without running the full exclusion process.
- The explicit Gaussian large-deviation route suggests that the correct frontier for shapes with growth between polynomial and exponential may interpolate between √(t log t) and linear scale t, and the sharp transition point could be identified by analyzing when the tail P(ζ_{t/d} > z) becomes so small that expectations of G(ζ_{t/d} − z) cease to dominate.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper studies the symmetric simple exclusion process (SSEP) on Z^d, d ≥ 2, starting from half-space step initial conditions whose cross-section is determined by shape functions g_2,...,g_d. It defines N_t as the number of particles that have moved beyond a level z = z(t) in the first coordinate, and asks whether convergence of E[N_t] suffices for a Poisson limit, as in the corresponding system of independent particles. The main theorem (Theorem 4.1) shows that, under conditions (A)-(C) on the shape functions and uniform boundedness of E[N_t], the correlation quantity E_t = S_t + C_t tends to 0: unconditionally for d ≥ 4, and in d = 2,3 under an additional decay condition (4.2). Corollary 4.3 converts this into a Poisson limit for N_t when lim E[N_t] exists. For polynomial shapes g_i(u) = c_i u^{α_i} + r_i, Theorem 5.1 identifies the scaling z = b_t(x+a_t), proves N_t converges to a Poisson law with mean M e^{-βx}, and consequently obtains a Gumbel limit for the maximum and explicit finite-sum limits for all order statistics. The proofs are built on the stirring construction, self-duality, negative association and the strong Rayleigh property, and are carried out in Sections 6-9 with auxiliary lemmas stated in detail.
Significance. If correct, this is a substantial contribution: it provides the first extreme-value and Poisson-statistics limits for the exclusion process in dimensions d > 1, showing that the statistics of the particle front match those of independent particles for a nontrivial class of interacting initial conditions. The d ≥ 4 result is an unconditional statement under natural hypotheses, and the polynomial-shape theorem gives explicit scaling and dimension-dependent constants that are directly checkable. A particular strength is that the proofs are largely self-contained and include quantitative rates of decay for S_t and C_t in Lemma 5.4; no free parameters are fitted and the main predictions are sharp and falsifiable. The d = 2,3 part is correctly proved only under the additional condition (4.2), whose necessity is explicitly left open; this qualification must be reflected in the abstract and introductory claims.
major comments (1)
- [Abstract, Theorem 4.1, Eq. (4.2), Section 2.4] The abstract states without qualification that 'when lim_{t→∞} E[N_t] exists, correlations between particles beyond z vanish so as to allow convergence of N_t to the same Poisson distribution' as for independent particles. This is stronger than what is proved. In d = 2,3, Theorem 4.1 and Corollary 4.3 require the additional condition (4.2) on E_0[G(ζ_{t/d}-z)1(ζ_{t/d}>z)], and the necessity of (4.2) is left open in Section 2.4 and Remark 4.4(a). The paper should qualify the abstract and the corresponding sentences in the Introduction and Section 2.2, either by stating that in d = 2,3 the result is conditional on (4.2) or by restricting the unqualified statement to d ≥ 4. This is a scope correction rather than a mathematical error: Theorem 5.1 is unaffected because Lemma 5.4 verifies (4.2) for polynomial shapes, but the advertised generality of the main theorem needs to be corrected before publication.
minor comments (3)
- [Section 7, Eq. (7.1)] The display preceding the proof of Theorem 4.1 reads γ_d(t)(E_0[G(ζ_{t/d}-z)1(ζ_{t/d}>z)])^2 = 0 for d ∈ {2,3}; the intended statement is that this quantity tends to 0 as t → ∞. As written, the equality makes the subsequent argument in (7.2) meaningless. Please correct '= 0' to '→ 0'.
- [Section 2.2, Eq. (2.4)] The displayed Gumbel rate appears to have (d-1)! ∏_i(2c_i) d^{d+1/2} e^{-dx} in the exponential. Theorem 5.1 and Remark 5.2(c) give the same rate with d^{d+1/2} in the denominator, i.e. (d-1)!∏_i(2c_i) d^{-(d+1/2)} e^{-dx}. If the inline formula is not a typesetting artifact, the exponent should be corrected.
- [Section 8, Lemma 8.2 and its proof] In the definition of \anH^i_{k,s}(m), the product over l ≠ i contains ζ_{u/d} with an undefined variable u; it should be ζ_{s/d} to agree with the rest of the lemma. A similar u/s typo appears in the proof of (8.17), where one expectation uses ζ_{u/d} instead of ζ_{s/d}.
Circularity Check
No load-bearing circularity: the Poisson and Gumbel limits are derived from new correlation bounds and explicit random-walk asymptotics; self-citations to [4] are proof templates, not circular inputs.
full rationale
The central derivation is self-contained against the stated model. Theorem 4.1 and Corollary 4.3 rest on new bounds on S_t and C_t (Propositions 4.5, 4.6, and Corollary 4.7), which are proved from self-duality, negative association, and random-walk estimates; the Poisson criterion (Lemma 3.8) is imported from known strong-Rayleigh theory, not from the paper's target result. Condition (4.2) in d=2,3 is an explicit extra hypothesis: the paper proves (4.1) under (4.2), and for polynomial shapes Lemma 5.4 verifies it by direct bounds. The abstract's phrase 'when lim E[N_t] exists' is broader than Theorem 4.1 in d=2,3, but this is a scope/overstatement concern, not circularity. For Theorem 5.1, the Gumbel and order-statistic limits follow from a genuine Poisson convergence proof plus the explicit computation E[N_t] -> M e^{-beta x} in Lemma 5.3; the centering constants a_t and b_t are derived from that computation, not fitted to the limit law. Self-citations to [4] occur as proof templates (Lemma 2.1 and Lemma 3.9) and as a moderate-deviation estimate; the load-bearing arguments are reproduced in the text or cited to external sources, so they are not circular. No equation assumes the target result, and no fitted parameter is renamed as a prediction. Therefore no circularity steps are identified; the score 2 only reflects minor, non-load-bearing self-citations.
Assumptions & free parameters
assumptions (6)
- domain assumption Self-duality of SSEP and the stirring construction, equations (2.5) and (3.6)-(3.8).
- domain assumption Strong Rayleigh property and negative association of SSEP occupation variables, Lemmas 3.6 and 3.8.
- domain assumption Conditions (A)-(C) on the shape functions g_i: nondecreasing, derivative regularity (B), and subexponential growth g'/g to 0 (C).
- ad hoc to paper Condition (4.2) in d=2,3: E0[G(zeta_{t/d}-z) 1(zeta_{t/d}>z)] = o((log t)^{-1/2}) for d=3 and o(t^{-1/4}) for d=2.
- standard math Random walk large-deviation comparison, Lemma 9.4, requiring z = o(t^{2/3}).
- domain assumption Infinite system with finite escape counts: sup_t E[N_t] < infinity and t^{-1/2} z tends to infinity.
Cite this review
Pith. "Pith review of Poisson statistics, vanishing correlations, and extremal particle limits for symmetric exclusion in d > 1." pith.science (2026). https://pith.science/paper/TCGEWFUT
@misc{pith2026250110522,
author = {Pith},
title = {Pith review of: Poisson statistics, vanishing correlations, and extremal particle limits for symmetric exclusion in d > 1},
year = {2026},
howpublished = {\url{https://pith.science/paper/TCGEWFUT}},
note = {Machine review of arXiv:2501.10522}
}
abstract
We consider the symmetric simple exclusion system on $\mathbb{Z}^d$, $d \ge 2$, starting from a class of ``step'' initial conditions in which particles are constrained within a half-space. One may count the number $N_t$ of particles that have moved beyond a distance $z = z(t)$ into the initially-empty half of $\mathbb{Z}^d$ at time $t$. We show in large generality that when $\lim_{t\to\infty} E[N_t]$ exists, correlations between particles beyond $z$ vanish as $t \to \infty$ so as to allow convergence of $N_t$ to the same Poisson distribution one would get were the particles allowed to move independently. When the initial condition constrains a region of polynomial growth, we identify $z(t)$ and the limit of $E[N_t]$ explicitly. As a consequence of the limit, we obtain a Gumbel limit distribution for the extremal particle position, as well as the limiting distributions of all order statistics.
Figures
Forward citations
Cited by 1 Pith paper
-
Point process convergence of extremes in $K$-symmetric exclusion
For K-symmetric exclusion from a step profile, the rescaled point process of extreme particles converges to a Poisson random measure with intensity proportional to e^{-x} dx.
Reference graph
Works this paper leans on
-
[4]
M. Conroy and S. Sethuraman. Gumbel laws in the symmetric exclu sion process. Com- mun. Math. Phys. 402: 723–764 (2023). 41
work page 2023
-
[1]
R. Arratia. The motion of a tagged particle in the simple symmetric e xclusion system on Z 1. Ann. Probab. 11(2): 362–373 (1983)
work page 1983
-
[2]
T. Bodineau and B. Dagallier. Large deviations for out of equilibrium correlations in the symmetric simple exclusion process. Electron. J. Probab. 29: 1–96 (2024)
work page 2024
- [3]
-
[5]
J. B. Hass, A. N. Caroll-Godfrey, I. Corwin, E. I. Corwin. Anoma lous fluctuations of extremes in many-particle diffusion. Phys. Rev. E 107 L022101 (2023)
work page 2023
- [6]
-
[7]
W. Feller. An Introduction to Probability Theory , vol. 2, ed. 2. Wiley, New York (1971)
work page 1971
-
[8]
B. Gess and V. Konarovskyi. A quantitative central limit theorem for the simple sym- metric exclusion process. arXiv:2408.01238 (2024)
arXiv 2024
Show all 19 references
-
[9]
Jansen and N
S. Jansen and N. Kurt. On the notion(s) of duality for Markov pr ocesses. Probab. Surveys 11: 59–120 (2014)
2014
-
[10]
Kipnis and C
C. Kipnis and C. Landim. Scaling limits of interacting particle systems. Grundlehren der Mathematischen Wissenschaften 320. Springer-Verlag, Berlin (1991)
1991
-
[11]
G. F. Lawler and V. Limic. Random Walk: A Modern Introduction . Cambridge Studies in Advanced Mathematics, vol. 123. Cambridge University Press, Ca mbridge (2010)
2010
-
[12]
T. M. Liggett. Interacting Particle Systems . Springer, Berlin (2005)
2005
-
[13]
T. M. Liggett. Stochastic Interacting Systems: Contact, Voter, Exclusio n Processes. Springer-Verlag, Berlin (1999)
1999
-
[14]
T. M. Liggett. Distributional limits for the symmetric exclusion pr ocess. Stoch. Process. Appl. 119: 1–15 (2009)
2009
-
[15]
Mikosch and J
T. Mikosch and J. Yslas. Gumbel and Fr´ echet convergence ofthe maxima of independent random walks. Adv. Appl. Prob. 52: 213–236 (2020)
2020
-
[16]
Redig and E
F. Redig and E. Saada. Particle systems with sources and sinks. Markov. Proc. Rel. Fields; arXiv:2403.01893 (2025)
2025 arXiv
-
[17]
H. Spohn. Large Scale Stochastic Dynamics. Springer-Verlag, Berlin (1991)
1991
-
[18]
Vandenberg-Rodes
A. Vandenberg-Rodes. A limit theorem for particle current in th e symmetric exclusion process. Electron. Commun. Probab. 15: 240–252 (2010)
2010
-
[19]
X. F. Xue and L. J. Zhao. Moderate deviations for the SSEP with a slow bond. J. Stat. Phys. 182: 48 (2021). 42
2021
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.