Pith. sign in

REVIEW 2 major objections 4 minor 1 cited by

Derivation of non-polynomial fractional diffusions from the generalized exclusion with a slow barrier

T0 review · 2 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash

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

Pith's one-line read This paper proves that a generalized exclusion process with long jumps and a slow barrier converges, in the hydrodynamic limit, to the unique weak solution of ∂_tρ = L^γ_κ F(ρ), where F is an absolutely convergent power series and the…

desk verdict Genuinely new hydrodynamic-limit result, but the proof as written has a real gap: the Replacement Lemma is not proved for the constant-profile entropy bound that Theorem 2.27 actually uses. read the letter →

arxiv 2412.12414 v1 pith:JCDZGRAH submitted 2024-12-16 math.PR

Pedro Cardoso , Patrícia Gonçalves , Gabriel Nahum This is my paper · ORCID
classification math.PR MSC 60K3535R1135S15
keywords hydrodynamiclimitgeneralizedexclusionprocessfractionalLaplacianregionalslowbarrierporousmediumequationpowerseriesnonlinearityreplacementlemma
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Hydrodynamic limits of interacting particle systems usually produce linear or polynomial equations. This paper constructs a particle system whose limiting density obeys a fractional diffusion equation with a genuinely non-polynomial nonlinearity: ∂_tρ = L^γ_κ F(ρ), where L^γ_κ is the distorted fractional Laplacian $κΔ^{{γ/2}}$ + (1-κ)$Δ^{{γ/2}}$_⋆ and F is an absolutely convergent power series of the form (2.8), so F can be an exponential rather than just x or $x^{2}$. The microscopic model is a symmetric generalized exclusion process with long jumps, up to N_e particles per site, and a slow barrier between two half-spaces; the barrier strength α_n becomes the parameter κ in the limiting operator, and the boundary condition at the barrier is fixed by the limit α = lim α_n and by whether α_n r_n^γ vanishes or stays positive. The result matters because it provides a constructive dictionary: choose the coefficients of F, build the jump rates from products of occupation variables, and the hydrodynamic limit returns the chosen PDE. The proof is an adaptation of the entropy method to infinite volume, with a replacement lemma that holds under entropy bounds relative to non-constant profiles.

What carries the argument

The argument is carried by a long-range gradient property of the rates. Identity (3.5) writes the current $c^{{(k),j}}$_{x,y}(η)[η(y)-η(x)] as $P^{{(k),j}}$(τ_y η) - $P^{{(k),j}}$(τ_x η) - ∇_j $A^{{(k),j}}$_{x,y}(η), with A antisymmetric; after summation by parts this turns the generator's integral term into a discrete version of the distorted fractional Laplacian acting on the local function F_n(η). The other essential mechanism is the Replacement Lemma (Lemma 3.15), proved by one-block and two-block estimates, which allows products of occupation variables to be replaced by products of empirical box averages under entropy bounds relative to non-constant profiles; this is what extracts the nonlinear function F from the microscopic products.

What would settle it

Take a case covered by the theorem, e.g. F(ρ)=e^ρ (b_k^+ = 1/k!, b_k^- = 0), d=1, γ=3/2, N_e=1, α_n → 1/2, and initial density a step function. Simulate the accelerated process at two large system sizes and compare the empirical density with a high-resolution numerical solution of (2.32); agreement at the barrier would confirm the boundary-condition classification, while disagreement in the bulk or across the barrier beyond discretization error would falsify the hydrodynamic-limit claim.

Watch

Extended reading notes

Core claim

Under Hypotheses 2.25, the n^γ-time-accelerated process converges, in probability and with respect to the empirical measure, to the unique weak solution of ∂_tρ = L^γ_κ F(ρ), where L^γ_κ = $κΔ^{{γ/2}}$ + (1-κ)$Δ^{{γ/2}}$_⋆ and F is the power series (2.8). The limiting equation is (2.32) when α = lim α_n lies in R_+\setminus{1}, (2.33) when α=1, and (2.34) or (2.35) when α=0, depending on whether α_n r_n^γ tends to 0 or has a positive limit, with r_n^γ defined in (2.42). The nonlinearity is realized microscopically by jump constraints that are linear combinations of products of occupation variables; the novel step is a replacement lemma that, under an entropy bound with non-constant profiles, replaces products of occupancies by products of local density averages, so the martingale's integral term converges to ∫ L^γ_κ F(ρ) G. Uniqueness is proved in the regimes covered by the theorem; in the α=0, γ∈(1,2) case with lim α_n r_n^γ ∈ (0,∞], the paper proves tightness and that all limit points solve (2.35) but leaves uniqueness open.

Load-bearing premise

The load-bearing premise is that the coefficient-weighted derivative sum f'_∞ = Σ_{k≥1} k(|b_k^+|+|b_k^-|) $N_e^{{k-1}}$ is finite, which keeps the microscopic rates bounded and the macroscopic F regular enough for energy estimates and uniqueness; this fails for fast-diffusion nonlinearities such as F(ρ)=ρ^m with m∈(0,1).

Editorial extensions

If this is right

  • Any absolutely convergent F of the form (2.8), including F(ρ)=e^ρ, can be realized as the hydrodynamic nonlinearity of a symmetric long-range exclusion model, so rigorously derivable fractional PDEs are no longer limited to polynomials.
  • The slow barrier's limiting strength α selects the operator: α=1 gives the standard fractional Laplacian, α=0 gives the regional fractional Laplacian, and other values give a distorted mixture of the two.
  • The speed at which α_n approaches its limit selects the boundary condition: α_n r_n^γ → 0 yields vanishing fractional flux through the barrier, while a positive limit yields equality of the two one-sided fractional derivatives.
  • For fast diffusion F(ρ)=ρ^m with m∈(0,1), the weaker Hypotheses 2.28–2.29 still give tightness and the integral equations, even though full uniqueness uses the stronger slow-growth condition (2.39).
  • The replacement lemma works under the weaker entropy bound (3.29) with non-constant profiles, so the hydrodynamic limit does not require initial measures close to a constant-density equilibrium.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The dictionary F ↔ jump rates is likely reversible: any sufficiently regular increasing diffusion coefficient with an absolutely convergent expansion of the type (2.8) should be realizable by some choice of b_k^±, which would extend the construction to nonlinearities such as D(ρ)=1/(1+ρ) or D(ρ)=ρ^m with m>2.
  • The threshold r_n^γ = 1_{γ<1} + log n 1_{γ=1} + n^{γ-1} 1_{γ>1} probably governs slow bonds in other long-range systems: the barrier is effectively insulating when α_n r_n^γ → 0 and transparent when the product has a positive limit, a classification one could test in zero-range or inclusion processes.
  • A practical check of the replacement mechanism: choose F with moderate coefficients so that (2.39) is close to failing, for example b_k^+ ≈ 1/k^{1+δ} with δ small, and look for finite-n corrections near the barrier; the proof predicts these corrections vanish only after the two limits n→∞ then ǫ→0, so simulations should show slow convergence in ǫ.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 4 minor

Summary. The paper derives, by means of the entropy method, the hydrodynamic limit of a generalized exclusion process with long jumps and a slow barrier. The limiting density is shown to solve a fractional nonlinear equation ∂tρ = L^γ_κ F(ρ), where F is a power series (2.8), L^γ_κ = κΔ^{γ/2} + (1−κ)Δ^{γ/2}_⋆, and the boundary conditions depend on α = lim α_n through Definitions 2.16–2.19. The main result, Theorem 2.27, covers polynomial and non-polynomial F under Hypotheses 2.25 and the entropy bound (2.41), with a partial statement (tightness only rather than full convergence) in the regime α = 0, γ ∈ (1,2), lim α_n r_n^γ ∈ (0,∞] because uniqueness of weak solutions for (2.35) is left open. The proof is organized as tightness (Section 4), energy estimates (Section 5), characterization of limit points via Dynkin's martingale (Section 6), and a replacement lemma (Section 7), with substantial appendices on fractional operators, discrete convergences, and measure-theoretic tools.

Significance. If completed, the paper would be a significant contribution: it gives rigorous microscopic derivations of non-polynomial fractional porous-medium-type equations, extends the slow-barrier analysis to a general scaling sequence α_n, and treats arbitrary dimension d ≥ 1 with occupation number Ne ≥ 1. The paper is carefully structured, with the solution notions, test-function spaces, and hypotheses set out precisely, and with detailed appendix material on the infinite-volume technicalities. A particular strength is the explicit separation of the hypotheses that are needed for tightness, energy estimates, uniqueness, and the replacement lemma, rather than hiding them in vague regularity conditions. The identified gap concerns the proof of the replacement lemma under one of its stated hypotheses; it is localized and likely repairable, but it currently affects the proof of the main theorem for nonlinear F.

major comments (2)
  1. [Section 7 / Lemma 3.15] Theorem 2.27 assumes the relative entropy bound (2.41), which together with Hypothesis 2.25(ii) is exactly Hypothesis 3.12(i). The Replacement Lemma 3.15 is stated under Hypothesis 3.12, but Section 7 proves it only for Hypothesis 3.12(ii): after Proposition 7.5 the text says “We leave the case where Hypothesis 3.12 (i) holds for the reader”. This is not a cosmetic omission, because Lemma 3.15 is used in the nonlinear case of Proposition 5.18 (Step III, for the energy estimate (5.3)) and in Proposition 6.5 (characterization of limit points). Consequently Theorem 2.27 is not fully proved for nonlinear F under its own hypotheses. Since h ≡ θ is an admissible element of Ref and (5.26) gives an exact Dirichlet-form identity for constant profiles, I expect the missing case to be repairable, but the proof as written should not rely on an unproved case of a central lemma.
  2. [Proposition 5.18] Proposition 5.18 is stated under (2.39) and (2.41) (plus the relevant choice of O), but its proof for nonlinear F applies Lemma 3.15. Under Hypothesis 3.12(i), Lemma 3.15 additionally requires Hypothesis 2.25(ii); this hypothesis is not listed among the assumptions of Proposition 5.18. The gap disappears inside Theorem 2.27 because Hypothesis 2.25(ii) is assumed there, but as a standalone proposition the statement is missing an assumption. Please add Hypothesis 2.25(ii) to Proposition 5.18, or replace the reference to Lemma 3.15 by a version whose hypotheses are met.
minor comments (4)
  1. [Section 1.1] The phrase “by means of a of a partial differential equation” contains a duplicated article and should be corrected.
  2. [Proposition 3.8] The statement writes “Y^{n,G}_F ≲ …” while the definitions in (3.26) and (3.27) use the notation Y^{n,γ}_F(G) and Y^{n,γ}_S(G); please unify the notation.
  3. [Equation (7.37)] In the definition of Ω_w, the second case appears to repeat “if k⋆ = k+” where it should read “if k⋆ = k−”.
  4. [Throughout] The notation Λ^γ_κ is introduced in (2.5) but the text sometimes writes Λ^γ_α, Λ^γ_{n,α}, and Λ^γ_κ in adjacent formulas; a notational index or table would improve readability.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the hydrodynamic limit is a genuine derivation whose limit PDE is not an input; the one admitted proof gap (constant-profile Replacement Lemma) is a completeness issue, not a circular reduction.

full rationale

The limiting equation is not an input: F is fixed by the coefficients (b±_k) in (2.8), the operator L^γ_κ by (2.31), and weak solutions by Definitions 2.16–2.19, independently of the convergence argument; no parameter is fitted to any data. The same coefficients appear in the generator (2.15), but this is model design—the microscopic rates are specified first and the limit flux is computed via the gradient identity (Prop. 3.2), the entropy bounds, and the Replacement Lemma, which convert the rate functional into powers ρ^k and (N_e−ρ)^k under the local equilibrium measures; the theorem's content is the convergence proof, not a data-driven prediction of F. Uniqueness (Lemma 2.20) is proved in-appendix (D.1) under 'F increasing' (which follows from Hypothesis 2.25(ii)) and is not imported from the authors' prior papers; [7,8,10,11,17] are invoked for techniques (Oleinik's method, the interpolating-model ansatz, Skorohod/entropy toolkit) and each has independent, external content. The claimed novelty, non-polynomial F, is explicitly an ansatz ('we will consider F in (2.6) given by (2.8)'), so Theorem 2.27 states convergence for that class. Two proof-completeness flags, weighed here but NOT treated as circularity: (1) Section 7 proves the Replacement Lemma only for the non-constant-profile setting and states: 'We leave the case where Hypothesis 3.12 (i) holds for the reader, but we observe that the arguments are totally analogous and, in fact, simpler.' Since Theorem 2.27 relies on Hypothesis 3.12(i)/(2.41) (via Props. 5.18 and 6.5), the stated proof is incomplete for the case actually needed; the omitted case is a special case of the proved one, so this is a gap, not a reduction to inputs. (2) Prop. 5.18 is stated under (2.39)+(2.41) only, yet its nonlinear Step III applies Lemma 3.15, whose Hypothesis 3.12(i) also requires Hypothesis 2.25(ii); under the theorem's own hypotheses that assumption is available, so the mismatch is fillable. Neither step exhibits a fit-renamed-as-prediction or a self-citation that forces the result.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No fitted parameters; the theorem is universal over all coefficient sequences satisfying the explicit hypotheses. The axioms are the model-defining assumptions and the standard mathematical background. No new physical entities are introduced.

assumptions (4)
  • domain assumption The generator L^γ_{n,α} in (2.16) defines a Markov process on Ω={0,...,Ne}^{Z^d}; the rates c^n_{x,y} are nonnegative (Hypothesis 2.25(ii)).
    Defines the microscopic model. Nonnegativity is imposed explicitly so that the generator is well-defined.
  • domain assumption Initial measures µn are associated to the profile g and satisfy H(µn|νθ) ≤ Cθ n^d for some θ∈(0,1) (Eq. (2.41)).
    Controls initial entropy and is the basis of the energy estimates in Section 5; standard in hydrodynamic limit proofs.
  • domain assumption Slow-growth condition f'_∞ = Σ k(|b+_k|+|b−_k|)N_e^{k-1} < ∞ (Hypothesis 2.25(i), Eq. (2.39)).
    Ensures truncation of rates is uniformly bounded and F is regular enough for uniqueness by Oleinik's method; also used in the replacement lemma under non-constant profile entropy. Excludes fast-diffusion nonlinearities.
  • standard math Standard analytic and probabilistic background: martingale theory, Skorokhod topology, Riesz representation, Portmanteau theorem, fractional Sobolev inequalities, Lebesgue differentiation theorem.
    Used throughout Sections 3-7 and the appendices; accepted background results assumed without proof.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Derivation of non-polynomial fractional diffusions from the generalized exclusion with a slow barrier." pith.science (2026). https://pith.science/paper/JCDZGRAH

@misc{pith2026241212414,
  author       = {Pith},
  title        = {Pith review of: Derivation of non-polynomial fractional diffusions from the generalized exclusion with a slow barrier},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JCDZGRAH}},
  note         = {Machine review of arXiv:2412.12414}
}
read the original abstract

In this article we derive in the hydrodynamic limit a generalized fractional porous medium equation, in the sense that the regional fractional Laplacian is applied to a function of the density given in terms of a power series, instead of a polynomial. The hydrodynamic limit is obtained considering a microscopic dynamics of random particles with long range interactions, but the jump rate highly depends on the occupancy near the sites where the interactions take place. This system is also studied in the presence of a "slow barrier" that hinders the flow of mass between two half-spaces.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. A gradient model for the Bernstein polynomial basis

    math.PR 2024-11 conditional novelty 7.0 of 10

    A symmetric exclusion process realizing Bernstein-polynomial diffusivities is constructed and claimed to satisfy the gradient condition, generalizing the Porous Media Model.

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages · cited by 1 Pith paper

  1. [17]

    From Exclusion to Slow and Fast Diffusion

    P . Gonçalves, G. Nahum, and M. Simon. From exclusion to slow a nd fast diffusion. arXiv preprint arXiv:2301.06585, 2023. NON-POL YNOMIAL FRACTIONAL DIFFUSIONS 81

    work page Pith review arXiv 2023

  2. [1]

    R. Adams. Sobolev spaces. Pure and Applied Mathematics, V ol. 65. Academic Pre ss [Harcourt Brace Jovanovich, Publishers], New Y ork-London, 1975

    work page 1975

  3. [2]

    Baldasso, O

    R. Baldasso, O. Menezes, A. Neumann, and R. Souza. Exclusion process with slow boundary . J. Stat. Phys., 167(5):1112–1142, 2017

    work page 2017

  4. [3]

    H. Bauer . Measure and integration theory, volume 26 of De Gruyter Studies in Mathematics . Walter de Gruyter & Co., Berlin, 2001. Translated from the German by Robert B. B urckel

    work page 2001

  5. [4]

    Billingsley .Convergence of probability measures

    P . Billingsley .Convergence of probability measures. John Wiley & Sons, Inc., New Y ork-London-Sydney , 1968

    work page 1968

  6. [5]

    Bonorino, R

    L. Bonorino, R. de Paula, P . Gonçalves, and A. Neumann. Hydro dynamics of porous medium model with slow reservoirs. J. Stat. Phys., 179(3):748–788, 2020

    work page 2020

  7. [6]

    H. Brezis. Functional analysis, Sobolev spaces and partial differential equations. Springer Science & Business Media, 2010

    work page 2010

  8. [7]

    Cardoso, R

    P . Cardoso, R. De Paula, and P . Gonçalves. Derivation of the fractional porous medium equation from a microscopic dynamics. Nonlinearity, 36(3):1840–1872, 2023

    work page 2023

Show all 27 references
  1. [8]

    Cardoso and P

    P . Cardoso and P . Gonçalves. Linear and nonlinear fractional pdes from interacting particle systems. to appear in Proceedings, 2024+

    2024

  2. [9]

    Cardoso, P

    P . Cardoso, P . Gonçalves, and B. Jiménez-Oviedo. Hydrodynamic behavior of long-range symmetric exclusion with a slow barrier: diffusive regime. to appear in AIHP, 2023+

    2023

  3. [10]

    Cardoso, P

    P . Cardoso, P . Gonçalves, and B. Jiménez-Oviedo. Hydrodyna mic behavior of long-range symmetric exclusion with a slow barrier: superdiffusive regime. to appear in Annali della Scuola Normale Superiore di Pisa, Classe di Scienze, 2023+

    2023

  4. [11]

    Cardoso, P

    P . Cardoso, P . Gonçalves, and B. Jiménez-Oviedo. Hydrodyna mics of a d-dimensionnal long jumps symmetric exclusion with a slow barrier . arXiv preprint arXiv:2304.01152, 2023

    2023 arXiv

  5. [12]

    Di Nezza, G

    E. Di Nezza, G. Palatucci, and E. V aldinoci. Hitchhiker’s gu ide to the fractional Sobolev spaces. Bull. Sci. Math., 136(5):521–573, 2012

    2012

  6. [13]

    Franceschini, P

    C. Franceschini, P . Gonçalves, G. Nahum, and B. Salvador . Hydrodynamical behavior for the generalized symmetric exclusion with open boundary . arXiv preprint arXiv:2201.10241, 2022

    2022 arXiv

  7. [14]

    Franco, P

    T . Franco, P . Gonçalves, and A. Neumann. Hydrodynamical beh avior of symmetric exclusion with slow bonds. Ann. Inst. Henri Poincaré Probab. Stat., 49(2):402–427, 2013

    2013

  8. [15]

    Franco and M

    T . Franco and M. T avares. Hydrodynamic limit for the SSEP wit h a slow membrane. J. Stat. Phys., 175(2):233– 268, 2019

    2019

  9. [16]

    Gonçalves, C

    P . Gonçalves, C. Landim, and C. T oninelli. Hydrodynamic limit for a particle system with degenerate rates. Annales de l’I.H.P .Probabilités et Statistiques, 45(4):887–909, 2009

    2009

  10. [18]

    Grisvard.Elliptic problems in nonsmooth domains, volume 69 of Classics in Applied Mathematics

    P . Grisvard.Elliptic problems in nonsmooth domains, volume 69 of Classics in Applied Mathematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, P A, 2011. Reprint of the 1985 original [ MR0775683], With a foreword by Susanne C. Brenner

    2011

  11. [19]

    Guan and Z

    Q. Guan and Z. Ma. Reflected symmetric α-stable processes and regional fractional Laplacian. Probab. Theory Related Fields, 134(4):649–694, 2006

    2006

  12. [20]

    M. Z. Guo, G. C. Papanicolaou, and S. R. S. V aradhan. Nonlinear diffusion limit for a system with nearest neighbor interactions. Comm. Math. Phys., 118(1):31–59, 1988

    1988

  13. [21]

    M. Jara. Hydrodynamic limit of particle systems with long ju mps. arXiv preprint arXiv:0805.1326v2, 2009

    2009 arXiv

  14. [22]

    Bernstein polynomials for functions of two variables of class c (k)

    Edward H Kingsley . Bernstein polynomials for functions of two variables of class c (k). Proceedings of the American Mathematical Society, 2(1):64–71, 1951

    1951

  15. [23]

    Kipnis and C

    C. Kipnis and C. Landim. Scaling limits of interacting particle systems, volume 320 of Grundlehren der mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Springer-V erlag, Berlin, 1999

    1999

  16. [24]

    O. Menezes. Limite hidrodinâmico para um processo de exclus ão: Pontos limites das medidas empíricas. Master’s thesis, 7 2013

    2013

  17. [25]

    Sethuraman and D

    S. Sethuraman and D. Shahar . Hydrodynamic limits for long-r ange asymmetric interacting particle systems. Electron. J. Probab., 23:Paper No. 130, 54, 2018

    2018

  18. [26]

    H. Spohn. Large Scale Dynamics of Interacting Particles . Springer-V erlag, 1991

    1991

  19. [27]

    Zeidler .Nonlinear functional analysis and its applications

    E. Zeidler .Nonlinear functional analysis and its applications. II/A. Springer-V erlag, New Y ork, 1990. Linear mono- tone operators, Translated from the German by the author and Leo F . Boron. Pedro Cardoso, I NSTITUTE FOR APPLIED MATHEMATICS UNIVERSITY OF BONN ENDENICHER ALL...

    1990

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.