REVIEW 2 major objections 5 minor 28 references
Functional Central Limit Theorem and SPDE for epidemic model with memory of the last infection and waning immunity
T0 review · 2 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read This paper proves that in a stochastic epidemic model with memory of the last infection and waning immunity, the rescaled fluctuations of the empirical infection measure converge to a unique Gaussian-driven stochastic limit, and derives…
desk verdict Useful FCLT/SPDE result for a structured epidemic model with memory, but Proposition 3.13 has a genuine topology gap in identifying the Gaussian limit. 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 is carried by the fluctuation empirical measure $\hat\mu^N$, by a per-individual reinfection counter $A_k^N(t)$ used to compare the interacting system with independent copies of the deterministic mean-field trajectory, and by the operator $R$ that encodes the expected jump effect of a reinfection through the memory kernel $K$. The counter comparison gives the $O(1/\sqrt{N})$ coupling estimate (3.5) under Assumption 2.3; this feeds the moment inequalities that yield tightness in the weighted Sobolev dual space $W^{-(m_d+1),\alpha}_0$. The martingale term $W^N$ has quadratic variation converging to the covariance (2.14) of the limiting Gaussian process, and a standard tightness criterion for Hilbert-space-valued processes turns these bounds into convergence.
What would settle it
Simulate the age-and-trait epidemic for growing $N$ with a memory kernel that obeys the model's normalization $\int_\Theta K(\theta,\tilde\theta)\,\nu(d\tilde\theta)=1$ for every $\theta$ but violates $\int_\Theta \sup_{\theta\in\Theta} K(\theta,\tilde\theta)\,\nu(d\tilde\theta)<\infty$, for example a heavy-tailed trait distribution; the theorem predicts tightness and the Gaussian limit (2.15), so observing divergence or subsequence-dependent limits would refute it.
Extended reading notes
Core claim
The central result is Theorem 2.6: under Assumptions 2.1–2.3, as $N\to\infty$, $\hat\mu^N=\sqrt{N}(\mu^N-\mu)$ converges in distribution in $D(\mathbb{R}_+, W^{-(m_d+1),\alpha}_0)$ to the unique continuous solution $\hat\mu$ of the stochastic integral equation (2.15), where $W$ is a centered Gaussian process with covariance (2.14). When $\hat\mu_t$ has a density, the density solves the SPDE (2.16); equivalently, the pair $(\hat F,\hat S)$—the fluctuations of the average force of infection and average susceptibility—solves the Volterra system (2.18)–(2.19). The paper thus claims that memory of previous infections changes the fluctuation limit only through the kernel $K$, and that the no-memory FCLT is recovered when $K(\theta,\tilde\theta)=K(\tilde\theta)$.
Load-bearing premise
The proof leans on the technical condition that the memory kernel $K$ satisfies $\int_\Theta \sup_{\theta\in\Theta} K(\theta,\tilde\theta)\,\nu(d\tilde\theta)<\infty$; the coupling estimate that drives the whole tightness argument collapses if this integrability bound fails.
Editorial extensions
If this is right
- Fluctuations of the empirical infection measure around the FLLN limit are of order $1/\sqrt{N}$ and converge to a Gaussian process, so forecast bands around mean-field epidemic curves become mathematically justified.
- The limiting covariance (2.14) is explicit in terms of the deterministic solution, the infectivity and susceptibility functions, and the memory kernel, so it can be computed without simulating the full stochastic system.
- When the limit measure has a density, the SPDE (2.16) gives a closed evolution equation for the fluctuation density, including the boundary condition at infection age zero.
- The Volterra representation (2.18)–(2.19) shows that all fluctuations are captured by the pair $(\hat F,\hat S)$, and setting $K(\theta,\tilde\theta)=K(\tilde\theta)$ recovers the earlier no-memory FCLT.
Reading between the lines
- The sup-integrability condition on $K$ is likely stronger than needed; one could test whether the same Gaussian limit emerges under weaker $L^1$-type kernels, which would extend the theorem to heavier-tailed trait distributions.
- The explicit covariance (2.14) makes the limit directly usable for uncertainty quantification: compare the Gaussian fluctuations with Monte Carlo simulations of the $N$-individual process and use the mismatch to estimate the memory kernel.
- The same Hilbert-space tightness route should transfer to other mean-field models with individual memory and unbounded initial ages, such as age-structured Hawkes or branching models, where bounded-age assumptions are unrealistic.
- The SPDE (2.16) suggests a numerical discretization in age and trait; a particle method against the Gaussian limit would be a direct, implementable test of the theorem.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper establishes a Functional Central Limit Theorem for a stochastic epidemic model with infection-age structure, individual traits, and memory of the last infection. The fluctuation process \hat\mu^N = \sqrt{N}(\mu^N - \mu) is shown to converge, in weighted Sobolev spaces, to a continuous Gaussian process solving a stochastic integral equation; the paper also derives an SPDE formulation and an alternative representation through fluctuations of the force of infection and average susceptibility. The proof follows the standard route of tightness via Aldous's criterion, identification of limits via martingale and quadratic-variation arguments, and uniqueness via Gronwall's inequality, with many technical estimates delegated to previous papers by the author and coauthors.
Significance. If the main theorem is valid, the paper makes a useful contribution to the non-Markovian epidemic FCLT literature: it incorporates memory of previous infections through a trait kernel, avoids a bounded-initial-age assumption, and uses weighted Sobolev spaces to handle the fluctuations. The derived SPDE and the comparison with the memoryless case in [28] are valuable. However, the proof as written contains load-bearing gaps in the identification of the limiting covariance and in the well-posedness/uniqueness step, so the central claim is not yet fully established.
major comments (2)
- [Proposition 3.13] The convergence of the Doob-Meyer process H_φ(µ^N) to H_φ(µ) is not justified. The proof shows only that H_φ is Lipschitz with respect to the sup-norm total variation distance and then invokes Theorem 2.1 via the continuous mapping theorem. But Theorem 2.1 gives convergence in the Skorohod weak topology, not in total variation; when µ_0 has a density, µ_t has a density while µ^N_t is atomic, so ‖µ^N_t − µ_t‖_TV = 2 for every N. Moreover, λ and R^(2)φ are only bounded measurable under Assumptions 2.1–2.2, so H_φ is not continuous in the weak topology used in Theorem 2.1. Since this convergence identifies the covariance of the limiting Gaussian process in (2.14), the argument is load-bearing. A direct proof using the coupling estimates of Lemma 3.1 and Proposition 3.2 to obtain pointwise convergence of ⟨µ^N_s, f⟩ for bounded measurable f, followed by dominated convergence in s, would likely repair the gap, but the manuscript does not supply it.
- [Proposition 3.14 / Theorem 2.6] The uniqueness argument requires that Rφ_{t−s} belongs to the Sobolev space W^{md+2,α}_0 so that ⟨\hatµ_s, Rφ_{t−s}⟩ is a valid dual pairing for \hatµ_s ∈ W^{−(md+2),α}_0. The proof of (3.29) only gives the pointwise bound |Rφ(a,θ)| ≤ 2‖φ‖_{C_b^{md+2}}; it does not control derivatives in θ or decay in a. Indeed, for φ with compact support in a and γ not decaying, Rφ(a,θ) → γ(a,θ)∫ φ(0,θ̃)K(θ,θ̃)ν(dθ̃) as a→∞, so Rφ need not lie in the weighted Sobolev space at all. Thus equation (2.15) may not even be well defined for the stated class of test functions, and the Gronwall argument does not establish uniqueness. The authors need either stronger regularity/decay assumptions on K and γ, a different solution space, or a different uniqueness argument.
minor comments (5)
- [Definition 2.5] Definition 2.5 says W takes values in W^{md+1,α}, but the convergence result is in W^{−(md+1),α}; W should be a W^{−(md+1),α}-valued process, or the text should clarify the intended embedding.
- [Definition 2.5] In the covariance formula, the right-hand side uses φ_s and ψ_s while the left-hand side uses φ_t and ψ_{t′}; the notation should be reconciled, and there is a typo where ϕ_s appears in the definition of eR.
- [Assumption 2.3] The condition ∫_Θ sup_{θ∈Θ} K(θ,θ̃)ν(dθ̃) < ∞ is used crucially in Lemma 3.1 and Proposition 3.2, but the paper gives no examples or discussion of natural memory kernels satisfying it, nor does it discuss whether it can be relaxed. A short remark would improve applicability.
- [Theorem 2.6] The theorem says the test function φ is bounded measurable and C^1 in the first two variables, and also φ_t ∈ W^{md+2,α}_0; such a φ need not belong to W^{md+2,α}_0 from the C^1 assumption alone. The regularity assumptions on φ should be stated precisely.
- [Proposition 3.13] The statement says the convergence takes place in C(R_+, W^{−(md+1),α}_0), while W^N was defined and shown tight in W^{−md,α}; the embedding W^{−md,α} ↪ W^{−(md+1),α} should be explicitly invoked.
Circularity Check
The FCLT derivation is conditional on the companion FLLN [13] and on the author's prior moment estimates [28], but these are independent inputs; no step reduces the target conclusion to its own assumptions.
full rationale
The claimed FCLT is not circular. The covariance (2.14) of the limiting Gaussian process is not postulated; it is obtained as the limit of the quadratic variation (3.4) of the compensated-noise martingale W^N. The proof of convergence of the fluctuation process is via Aldous tightness (Prop. 3.9) and a Gronwall uniqueness argument (Prop. 3.14), and the identification of the limit's Doob-Meyer process (Prop. 3.13) uses the FLLN of the companion paper [13] as an external input. Heavy self-citation occurs — Theorem 2.1 is from the authors' companion paper, and Prop. 3.2 adapts [28, Prop. 5.4] — but these are independent prior theorems, not restatements of the target result; the equivalence with [28] in Prop. 2.9 is a consistency check. One genuine technical gap, not a circularity: in Prop. 3.13 the paper proves H_phi is Lipschitz in total variation and invokes Theorem 2.1, although Theorem 2.1 supplies only weak convergence in D(R_+, P(R_+ × Theta)); since mu^N_t is atomic while mu_t typically has a density, ||mu^N_t - mu_t||_TV = 2 for all N. This leaves the convergence of the quadratic variation unproved as written, but the missing argument (pointwise convergence of <mu^N_s, R^(2)phi> via the coupling estimates of Lemma 3.1/Prop. 3.2 plus dominated convergence) would still not make the result equivalent to its input. Thus the circularity score is low, reflecting only the heavy but non-circular reliance on the author's own prior work.
Assumptions & free parameters
assumptions (4)
- domain assumption The FLLN from Guerin and Zotsa-Ngoufack [13, Theorem 2.1]: mu^N converges to the unique solution mu of (2.5).
- domain assumption nu is absolutely continuous with respect to Lebesgue measure.
- standard math Sobolev embedding theorems (2.10)-(2.13), including Hilbert-Schmidt embeddings.
- standard math Aldous tightness criterion for Hilbert space-valued processes (Definition A.1) and the martingale characterization of Gaussian limits.
Cite this review
Pith. "Pith review of Functional Central Limit Theorem and SPDE for epidemic model with memory of the last infection and waning immunity." pith.science (2026). https://pith.science/paper/TQY2CVAQ
@misc{pith2026250515617,
author = {Pith},
title = {Pith review of: Functional Central Limit Theorem and SPDE for epidemic model with memory of the last infection and waning immunity},
year = {2026},
howpublished = {\url{https://pith.science/paper/TQY2CVAQ}},
note = {Machine review of arXiv:2505.15617}
}
read the original abstract
We study the fluctuations of a stochastic epidemic model with memory of the last infections, varying infectivity, and waning immunity, as introduced in Guerin and Zotsa-Ngoufack:arXiv preprint arXiv:2505.00601. The dynamics of the epidemic model are described by a measure-valued process with respect to infection age and individual traits. The Functional Law of Large Numbers (FLLN) is formulated as an integral equation, which is solved by a deterministic measure. In this article, we establish the Functional Central Limit Theorem (FCLT), capturing the fluctuations of the stochastic model around its deterministic limit. The limit of the FCLT is given by a nonlinear stochastic integral equation which is solved by a random signed-measure. We further derive the weak solution in the form of a stochastic partial differential equation (SPDE) and propose an alternative representation of the FCLT, as fluctuations in the average total force of infection and average susceptibility.
Reference graph
Works this paper leans on
-
[4]
Julien Chevallier. Fluctuations for mean-field interacting age-dependent Hawkes processes.Electronic Journal of Probability, 22, 2017
work page 2017
-
[9]
Stochastic epidemic models with varying infectivity and waning immunity
Rapha¨ el Forien, Guodong Pang,´Etienne Pardoux, and Arsene-Brice Zotsa-Ngoufack. Stochastic epidemic models with varying infectivity and susceptibility.arXiv preprint arXiv:2210.04667, 2022
work page Pith review arXiv 2022
-
[28]
Arsene Brice Zotsa Ngoufack. Functional central limit theorems for epidemic models with varying infectivity and waning immunity.ESAIM: Probability and Statistics, 29:45–112, 2025. AppendixA. A.1.Definition of tightness. Definition A.1.Let H be a separable Hilbert space. A sequence of stochastic process( X N )N ∈ D(R+, H)is tight if the following condition...
work page 2025
-
[13]
A stochastic epidemic model with memory of the last infection and waning immunity
H´ el` ene Gu´ erin and Arsene Brice Zotsa-Ngoufack. A stochastic epidemic model with memory of the last infection and waning immunity.arXiv preprint arXiv:2505.00601, 2025
work page Pith review arXiv 2025
-
[1]
Elsevier, 2003
Robert A Adams and John JF Fournier.Sobolev spaces. Elsevier, 2003
2003
-
[2]
John Wiley & Sons, 1999
Patrick Billingsley.Convergence of probability measures. John Wiley & Sons, 1999
1999
-
[3]
Tom Britton, Etienne Pardoux, Franck Ball, Catherine Laredo, David Sirl, and Viet Chi Tran.Stochastic epidemic models with inference, volume 132. Springer, 2019
work page 2019
-
[5]
A general age-dependent branching process
Kenny S Crump and Charles J Mode. A general age-dependent branching process. i.Journal of Mathematical Analysis and Applications, 24(3):494–508, 1968
work page 1968
Show all 28 references
-
[6]
John Wiley & Sons, 2009
Stewart N Ethier and Thomas G Kurtz.Markov processes: characterization and convergence. John Wiley & Sons, 2009
2009
-
[7]
A hilbertian approach for fluctuations on the mckean-vlasov model
Begona Fernandez and Sylvie M´ el´ eard. A hilbertian approach for fluctuations on the mckean-vlasov model. Stochastic processes and their applications, 71(1):33–53, 1997
1997
-
[8]
Epidemic models with varying infectivity.SIAM Journal on Applied Mathematics, 81(5):1893–1930, 2021
Rapha¨ el Forien, Guodong Pang, and´Etienne Pardoux. Epidemic models with varying infectivity.SIAM Journal on Applied Mathematics, 81(5):1893–1930, 2021
1930
-
[10]
From individual-based epidemic models to McKendrick-von Foerster PDEs: A guide to modeling and inferring COVID-19 dynamics.J
F´ elix Foutel-Rodier, Fran¸ cois Blanquart, Philibert Courau, Peter Czuppon, Jean-Jil Duchamps, Jasmine Gamblin, ´Elise Kerdoncuff, Rob Kulathinal, L´ eo R´ egnier, Laura Vuduc, Amaury Lambert, and Emmanuel Schertzer. From individual-based epidemic models to McKendrick-von Fo...
2022
-
[11]
Optimal vaccination policy to prevent endemicity: a stochastic model.Journal of Mathematical Biology, 90(1):1–55, 2025
F´ elix Foutel-Rodier, Arthur Charpentier, and H´ el` ene Gu´ erin. Optimal vaccination policy to prevent endemicity: a stochastic model.Journal of Mathematical Biology, 90(1):1–55, 2025
2025
-
[12]
Gu´ erin and A.-B
H. Gu´ erin and A.-B. Zotsa-Ngoufack. Stochastic epidemic model with memory of previous infection and waning immunity.Work in progress, 2025
2025
-
[14]
Weak convergence of sequences of semimartingales with applications to multitype branching processes.Advances in Applied Probability, 18(1):20–65, 1986
Anatole Joffe and Michel M´ etivier. Weak convergence of sequences of semimartingales with applications to multitype branching processes.Advances in Applied Probability, 18(1):20–65, 1986
1986
-
[15]
Spatial sir epidemic model with varying infectivity without movement of individuals: Law of large numbers.arXiv preprint arXiv:2412.01673, 2024
Armand Kanga and Etienne Pardoux. Spatial sir epidemic model with varying infectivity without movement of individuals: Law of large numbers.arXiv preprint arXiv:2412.01673, 2024
2024
-
[16]
McKendrick
William Ogilvy Kermack and Anderson G. McKendrick. A contribution to the mathematical theory of epidemics. Proceedings of the Royal Society of London. Series A., 115(772):700–721, 1927
1927
-
[17]
McKendrick
William Ogilvy Kermack and Anderson G. McKendrick. Contributions to the mathematical theory of epidemics. II.—The problem of endemicity.Proceedings of the Royal Society of London. Series A, containing papers of a mathematical and physical character, 138(834):55–83, 1932
1932
-
[18]
McKendrick
William Ogilvy Kermack and Anderson G. McKendrick. Contributions to the mathematical theory of epidemics–III. Further studies of the problem of endemicity. 1933.Proceedings of the Royal Society of London. Series A, containing papers of a mathematical and physical character, 14...
1933
-
[19]
Springer, 2006
David Nualart.The Malliavin calculus and related topics, volume 1995. Springer, 2006
1995
-
[20]
Functional central limit theorems for epidemic models with varying infectivity.Stochastics, pages 1–48, 2022
Guodong Pang and ´Etienne Pardoux. Functional central limit theorems for epidemic models with varying infectivity.Stochastics, pages 1–48, 2022
2022
-
[21]
Functional limit theorems for non-markovian epidemic models.The Annals of Applied Probability, 32(3):1615–1665, 2022
Guodong Pang and ´Etienne Pardoux. Functional limit theorems for non-markovian epidemic models.The Annals of Applied Probability, 32(3):1615–1665, 2022
2022
-
[22]
Functional law of large numbers and pdes for epidemic models with infection-age dependent infectivity.Applied Mathematics & Optimization, 87(3):50, 2023
Guodong Pang and ´Etienne Pardoux. Functional law of large numbers and pdes for epidemic models with infection-age dependent infectivity.Applied Mathematics & Optimization, 87(3):50, 2023
2023
-
[23]
Spde for stochastic sir epidemic models with infection-age dependent infectivity.arXiv preprint arXiv:2410.17200, 2024
Guodong Pang and Etienne Pardoux. Spde for stochastic sir epidemic models with infection-age dependent infectivity.arXiv preprint arXiv:2410.17200, 2024
2024 arXiv
-
[24]
Stochastic sir model with individual heterogeneity and infection-age dependent infectivity on large non-homogeneous random graphs.arXiv preprint arXiv:2502.04225, 2025
Guodong Pang, ´Etienne Pardoux, and Aur´ elien Velleret. Stochastic sir model with individual heterogeneity and infection-age dependent infectivity on large non-homogeneous random graphs.arXiv preprint arXiv:2502.04225, 2025. 24 ARSENE BRICE ZOTSA NGOUF ACK
2025 arXiv
-
[25]
PhD thesis, Universit´ e de Nanterre-Paris X, 2006
Viet Chi Tran.Mod` eles particulaires stochastiques pour des probl` emes d’´ evolution adaptative et pour l’approximation de solutions statistiques. PhD thesis, Universit´ e de Nanterre-Paris X, 2006
2006
-
[26]
Springer Science & Business Media, 2012
Kˆ osaku Yosida.Functional analysis, volume 123. Springer Science & Business Media, 2012
2012
-
[27]
Weighted norm inequality in the variable Lebesgue spaces for Bergman Projector on the unit ball ofC n.PhD thesis, Aix-Marseille Universit´ e, 2024
Arsene Brice Zotsa Ngoufack.Stochastic Epidemic model with varying infectivity and waning immunity: Law of large numbers and Central limit theorem. Weighted norm inequality in the variable Lebesgue spaces for Bergman Projector on the unit ball ofC n.PhD thesis, Aix-Marseille U...
2024
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.