REVIEW 3 major objections 5 minor 3 cited by
A sequential, online particle system approximating the McKean–Vlasov diffusion achieves the same N^{−1/2} statistical accuracy as the exchangeable system, with a factor-2 feedback in its Gaussian fluctuation limit.
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 →
T0 review · deepseek-v4-flash
2026-08-03 05:36 UTC pith:RGZ5TH4I
load-bearing objection Theorems 1–2 on incremental entropy and N^{-1/2} empirical convergence look solid; the Gaussian fluctuation theorem has a genuine gap at the a_j→a identification step, so the paper needs a serious referee rather than acceptance as is. the 3 major comments →
Quantitative Convergence for Sequential Interacting Diffusions via Incremental Relative Entropy
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the incremental relative entropy, which measures how far particle i's path law deviates from an independent McKean–Vlasov copy given the predecessors, decays as (i−1)^{−1}. This rate is optimal: even with i.i.d. inputs drawn from the limit law, the sampling variance of the interaction field forces exactly this scale. Summing the increments yields a global path-space entropy of order log N, and the same mechanism yields N^{−1/2} convergence of the empirical measure in H^{−β} and a Gaussian limit whose linear feedback term carries a factor 2, originating from the weighted sum ∑_{i=2}^N 1/√(i−1) ~ 2√N.
What carries the argument
The argument is carried by the incremental path-space relative entropy R_i(t), computed through a Girsanov identity as the quadratic energy of the drift mismatch between the conditional law of particle i and an independent limit copy. The predecessor empirical measure is replaced by an average of conditional laws, producing a martingale-difference term (handled by sub-Gaussian concentration under a decoupled reference law) and a predictable bias term (handled by Pinsker's inequality). An upper-envelope closure and a discrete Hardy inequality yield the sharp decay; a negative Sobolev energy estimate then converts the entropy control into the N^{−1/2} empirical-measure rate.
Load-bearing premise
The Gaussian fluctuation limit depends on the limit SPDE being well-posed in the weighted Sobolev space and on the asserted convergence a_j→a in the identification step; neither is proved in this paper.
What would settle it
Compute directly, for a one-dimensional linear example with explicit solvable kernel and initial density, the limiting variance of ∫φ dη_N at a fixed time; if the prefactor differs from 2, or if the asserted convergence a_j→a fails for a non-smooth initial density, the identification of the SPDE limit collapses. Alternatively, simulate the sequential system with i.i.d. initial data and compare the empirical covariance with the prediction of the factor-2 SPDE.
If this is right
- Sequential particle simulation with online refinement achieves the same N^{−1/2} empirical-measure accuracy as fully-coupled mean-field Monte Carlo, despite its O(N) marginal cost for adding particles.
- Tail blocks of m particles are chaotic with total-variation error O(√(m/N)); in particular, the terminal particle's path law is within O(N^{−1/2}) of the McKean–Vlasov law.
- The factor-2 feedback in the limiting SPDE is a testable signature of the directed architecture, distinguishing it from exchangeable and independent-copy systems.
- For general step-size sequences, the method predicts incremental entropy decay governed by the effective sample size θ_{i−1}^{−1}, with a phase transition at the divergence of ∑α_i.
Where Pith is reading between the lines
- The factor-2 accumulation mechanism likely persists in streaming and random-batch variants of sequential samplers, wherever predecessor averages contribute at the 1/√(i−1) fluctuation scale; a testable extension is to measure the feedback coefficient in such variants.
- The sharp incremental entropy bound suggests an analogous 1/(i−1) decay for time-marginal Fisher information, which could open a route to uniform-in-time and singular-kernel sequential propagation of chaos, a direction the paper sketches but does not fully develop.
- The Gaussian limit is the least settled part of the paper: the identification step relies on an unproved convergence a_j→a and on imported well-posedness of the limit SPDE, so a direct proof of uniqueness in the weighted Sobolev space would close the main remaining gap.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the sequential (lower-triangular) interacting diffusion system (1.4), in which particle i interacts only with the empirical measure of its predecessors. It introduces an incremental path-space relative entropy R_i(t) and proves Theorem 1: R_i(t) ≤ C_T (i−1)^{-1}(1+M_{0,i}), with the consequences S_N(T) ≲ log N and quantitative tail propagation of chaos. Theorem 2 establishes the canonical N^{-1/2} rate for the empirical measure in negative Sobolev norms, E[sup_t ||μ^N_t−ρ_t||_{H^{-β}}^2] + E∫||μ^N_t−ρ_t||_{H^{-β+1}}^2 dt ≤ C_T/N, under b≡0 and a regularity condition on K⋆ρ. Theorem 3 states a Gaussian fluctuation limit in C([0,T];H^{-β',-m}) for η^N=√N(μ^N−ρ), with the limiting SPDE (1.10) containing a factor-2 feedback term. Section 6 extends the incremental entropy method to general sequential weights and states a threshold phenomenon. The main technical machinery is a Girsanov identity, a martingale-difference replacement of predecessor empirical measures by averaged conditional measures, an upper-envelope closure, and a negative Sobolev energy estimate.
Significance. Theorems 1 and 2 are substantial and, as far as I can verify, correctly proven. The sharp (i−1)^{-1} decay for a directed non-exchangeable system and the logarithmic global entropy are genuinely new and well motivated by the online/streaming interpretation. The negative-Sobolev N^{-1/2} bound is also an honest improvement over what global entropy arguments would give. The factor-2 feedback in the fluctuation SPDE is an interesting and potentially important structural effect of the sequential architecture. However, the identification step in Theorem 3 currently contains a load-bearing gap, and the regularity assumptions for Theorem 3 are not sufficient as stated. The paper's proofs are largely transparent and the estimates are explicit, which is a strength; with a complete proof of the fluctuation limit the contribution would be significant.
major comments (3)
- [§5.1, Step 2] The identification of the factor-2 term rests on the assertion 'We have a_j→a as j→∞', where a_j = ∫_s^t ⟨ρ_r,(K*η^j_r)·∇φ⟩ dr and a = ∫_s^t ⟨ρ_r,(K*η_r)·∇φ⟩ dr. This is not proved and does not follow from the preceding tightness argument. Tightness gives at most subsequential convergence of the global fields η^N to some η; the prefix fields η^j = √j(μ^j−ρ) are not shown to converge to the same η, and no estimate controlling ||η^j−η|| is supplied. The Stolz–Cesàro passage requires convergence of the entire sequence a_j, not of a subsequence. This is the load-bearing step for the factor-2 term in (1.10). Without a proof of a_j→a, the martingale-problem identification of Theorem 3 is incomplete.
- [§5.1, first paragraph] The claim that ρ_0∈S(R^d) and 'standard parabolic bootstrapping' imply ρ∈C([0,T];S(R^d)) is not justified for bounded measurable K. In general K⋆ρ_t is only bounded measurable; it need not lie in W^{m,∞}, so high-order Sobolev bounds on ρ_t do not follow solely from smoothness of ρ_0. Consequently the well-definedness of the product ρ_t(K⋆η_t) in (1.10) and the estimates such as ||L_r φ||_{H^β} ≲ ||φ||_{H^{β'}} are not established. Theorem 3 should either impose an explicit regularity condition on K⋆ρ_t, as Theorem 2 does in (1.8), or assume smoothness of K. This is a second load-bearing gap for Theorem 3.
- [§6, Theorem 4] The general-weight result is stated as a theorem, but the proof is only a 'Sketch of proof.' The weighted upper-envelope closure is not written out, and the threshold phenomenon for r>1 relies on unproved assertions about θ_{i−1} and non-convergence of R_i. If Theorem 4 is to be part of the paper's formal results, it needs a complete proof; otherwise it should be presented as a conjecture or heuristic extension.
minor comments (5)
- [§2 and §3] There are frequent typos: 'Propsition', 'postive integer', 'inequlaity', 'concerntration', 'defintions', 'sup|b|<∞' in Section 6, etc. The paper would benefit from a careful proofreading pass.
- [§4, Theorem 2] Theorem 2 and the whole negative-Sobolev section assume b≡0, while the introduction and Theorem 1 allow a bounded drift b. This restriction is fine if stated clearly at the start of Section 4, but the presentation could make it more prominent.
- [§5.1, Step 3] Uniqueness is asserted via a one-line energy estimate and then attributed to [31, Section 4]. Please state precisely which theorem or proposition of [31] is being invoked, since the factor-2 drift and the weighted spaces are not literally the exchangeable setting of [31].
- [§5.1, Step 1] The tightness argument imports the martingale modulus bound from [31, Lemmas 3.1–3.2] without checking that the sequential drift decompositions verify the same hypotheses. At least the relevant verifications should be indicated.
- [§1.2, Remark 1.7] The comparison of the three SPDE limits is useful, but the 'add–subtract trick' displaying the decomposition of the sequential interaction term contains a typo: the index in the last sum should be 1/√(i−1)δ_{X^i}(K⋆η^{i−1}), which is later written correctly in Section 5.
Circularity Check
No circular derivation: Theorems 1 and 2 are proved from the dynamics; Theorem 3 has an unproved convergence assertion (a_j→a) that is a correctness gap, not a constructional circularity.
full rationale
The paper's main derivation chain is not circular. Theorem 1 follows from the Girsanov identity (3.2), which is not an identity in the target quantities alone: the drift mismatch Δ^i_s is estimated through conditional-measure replacement, martingale-difference sub-Gaussian bounds, and an upper-envelope closure. No fitted parameter is renamed as a prediction, and the 1/(i−1) bound is an honest estimate rather than a restatement of the definition of R_i. Corollaries 1.1–1.2 and Theorem 2 are derived consequences of these bounds together with the energy identity (4.4) and the negative-Sobolev estimates; they do not presuppose the conclusions. The self-citations to [30] and [31] are methodological and technical borrowings from prior published work, and the central quantitative estimates are derived here from stated assumptions rather than imported. The only load-bearing issue I see is in Theorem 3, Step 2: the text asserts "We have a_j → a as j→∞" for the prefix fluctuation fields η^j without proof and without an estimate relating η^j to the subsequential limit η; this convergence is needed for the Stolz–Cesàro factor 2. This is a genuine omitted justification / correctness gap, but it is not a definitional equivalence or a fitted-input-called-prediction reduction, so it does not make the derivation circular in the sense of this review. The score 2 reflects minor self-citation and reliance on prior machinery, not circularity of the central proofs.
Axiom & Free-Parameter Ledger
axioms (5)
- domain assumption σ is constant invertible; b and K are bounded measurable.
- domain assumption V_T = sup_{t≤T} ||K⋆ρ_t||_{W^{m,∞}} < ∞ with m=⌈β⌉.
- ad hoc to paper ρ_0 is a Schwartz function.
- domain assumption Well-posedness (uniqueness) of the limit SPDE (1.10) in H^{-β',-m}.
- standard math Girsanov theorem, Itô formula, Azuma–Hoeffding inequality, discrete Hardy inequality, Sobolev embedding.
read the original abstract
We study a lower-triangular system of interacting diffusions in which particle \(i\) interacts only with its predecessors through the empirical measure \(\mu^{i-1}_t\). This gives a directed, non-exchangeable approximation of the same McKean--Vlasov diffusion as the classical exchangeable particle system. We introduce an incremental path-space relative entropy adapted to the causal structure, $$ R_i(T) = H \left(P^{1:i}_{[0,T]}\,\middle|\,P^{1:i-1}_{[0,T]}\otimes \bar P_{[0,T]}\right), $$ and prove the sharp estimate \(R_i(T)\lesssim (i-1)^{-1}\). Furthermore, we obtain convergence of the empirical measure to the McKean--Vlasov law at the canonical \(N^{-1/2}\) scale in negative Sobolev norms. The proof combines a Girsanov representation, a martingale-difference replacement of predecessor empirical measures by averaged conditional measures, an upper-envelope closure, and a negative Sobolev energy estimate.
Figures
Forward citations
Cited by 3 Pith papers
-
A Logarithmic Fluctuation Hierarchy for Sequential Interacting Diffusions
For lower-triangular interacting diffusions, the N^{-1/2} fluctuation field converges to an infinite logarithmic hierarchy that couples Y^n to Y^{n+1}, not to the closed fluctuation SPDE of exchangeable mean-field systems.
-
Quantitative Propagation of Chaos and Fluctuations for Kinetic McKean--Vlasov SDEs with Singular Interaction Kernels
Kinetic McKean–Vlasov systems with singular Kato-class kernels enjoy path-space entropy chaos at rate k/N and a Gaussian fluctuation CLT with N^{-1/6} Berry–Esseen projections.
-
Universal Central Limit Theorem for non-exchangeable interacting diffusions
The global fluctuation field of non-exchangeable interacting diffusions converges to the universal mean field SPDE limit under deterministic denseness conditions on the interaction matrix, with a sharp n^{-1/2} threshold.
Reference graph
Works this paper leans on
-
[1]
John Wiley & Sons, 2016
Noga Alon and Joel H Spencer.The probabilistic method. John Wiley & Sons, 2016
2016
-
[2]
New particle repre- sentations for ergodic McKean–Vlasov SDEs.ESAIM: Proceedings and Surveys, 65:145–185, 2019
Houssam AlRachid, Mireille Bossy, Cristiano Ricci, and Lukasz Szpruch. New particle repre- sentations for ergodic McKean–Vlasov SDEs.ESAIM: Proceedings and Surveys, 65:145–185, 2019
2019
-
[3]
Springer Science & Business Media, 2008
Luigi Ambrosio, Nicola Gigli, and Giuseppe Savar´ e.Gradient flows: in metric spaces and in the space of probability measures. Springer Science & Business Media, 2008
2008
-
[4]
Graphon mean field systems.The Annals of Applied Probability, 33(5):3587–3619, 2023
Erhan Bayraktar, Suman Chakraborty, and Ruoyu Wu. Graphon mean field systems.The Annals of Applied Probability, 33(5):3587–3619, 2023
2023
-
[5]
Self-interacting diffusions.Probability Theory and Related Fields, 122(1):1–41, 2002
Michel Bena ¨ ım, Michel Ledoux, and Olivier Raimond. Self-interacting diffusions.Probability Theory and Related Fields, 122(1):1–41, 2002. 46
2002
-
[6]
A stochastic particle method for the McKean–Vlasov and the Burgers equation.Mathematics of Computation, 66(217):157–192, 1997
Mireille Bossy and Denis Talay. A stochastic particle method for the McKean–Vlasov and the Burgers equation.Mathematics of Computation, 66(217):157–192, 1997
1997
-
[7]
Large-scale machine learning with stochastic gradient descent
L´ eon Bottou. Large-scale machine learning with stochastic gradient descent. InProceedings of COMPSTAT’2010: 19th International Conference on Computational StatisticsParis France, August 22-27, 2010 Keynote, Invited and Contributed Papers, pages 177–186. Springer, 2010
2010
-
[8]
Mean field limit and quantitative estimates with singular attractive kernels.Duke Mathematical Journal, 172(13):2591–2641, 2023
Didier Bresch, Pierre-Emmanuel Jabin, and Zhenfu Wang. Mean field limit and quantitative estimates with singular attractive kernels.Duke Mathematical Journal, 172(13):2591–2641, 2023
2023
-
[9]
Sequentially interacting Markov chain Monte Carlo methods.The Annals of Statistics, 38(6):3387–3411, 2010
Anthony Brockwell, Pierre Del Moral, and Arnaud Doucet. Sequentially interacting Markov chain Monte Carlo methods.The Annals of Statistics, 38(6):3387–3411, 2010
2010
-
[10]
Interacting Markov chain Monte Carlo methods for solving nonlinear measure-valued equations
Pierre Del Moral and Arnaud Doucet. Interacting Markov chain Monte Carlo methods for solving nonlinear measure-valued equations. 2010
2010
-
[11]
Sequential propagation of chaos.arXiv preprint arXiv:2301.09913, 2023
Kai Du, Yifan Jiang, and Xiaochen Li. Sequential propagation of chaos.arXiv preprint arXiv:2301.09913, 2023
Pith/arXiv arXiv 2023
-
[12]
John Wiley & Sons, 2011
Paul Dupuis and Richard S Ellis.A weak convergence approach to the theory of large deviations, volume 902. John Wiley & Sons, 2011
2011
-
[13]
Propagation of chaos for the 2d viscous vortex model.Journal of the European Mathematical Society, 16(7):1423–1466, 2014
Nicolas Fournier, Maxime Hauray, and St´ ephane Mischler. Propagation of chaos for the 2d viscous vortex model.Journal of the European Mathematical Society, 16(7):1423–1466, 2014
2014
-
[14]
Uniform in time propagation of chaos for the 2d vortex model and other singular stochastic systems.Journal of the European Mathematical Society, pages 1–28, 2024
Arnaud Guillin, Pierre Le Bris, and Pierre Monmarch´ e. Uniform in time propagation of chaos for the 2d vortex model and other singular stochastic systems.Journal of the European Mathematical Society, pages 1–28, 2024
2024
-
[15]
Strong convergence of propagation of chaos for mckean–vlasov sdes with singular interactions.SIAM Journal on Mathematical Analysis, 56(2):2661–2713, 2024
Zimo Hao, Michael R¨ ockner, and Xicheng Zhang. Strong convergence of propagation of chaos for mckean–vlasov sdes with singular interactions.SIAM Journal on Mathematical Analysis, 56(2):2661–2713, 2024
2024
-
[16]
Cambridge university press, 1952
Godfrey Harold Hardy, John Edensor Littlewood, and George P´ olya.Inequalities. Cambridge university press, 1952
1952
-
[17]
A review of the mean field limits for vlasov equations.Kinetic & Related Models, 7(4):661, 2014
Pierre-Emmanuel Jabin. A review of the mean field limits for vlasov equations.Kinetic & Related Models, 7(4):661, 2014
2014
-
[18]
Quantitative estimates of propagation of chaos for stochastic systems withw −1,∞ kernels.Inventiones mathematicae, 214(1):523–591, 2018
Pierre-Emmanuel Jabin and Zhenfu Wang. Quantitative estimates of propagation of chaos for stochastic systems withw −1,∞ kernels.Inventiones mathematicae, 214(1):523–591, 2018
2018
-
[19]
The mean-field limit of sparse networks of integrate- and-fire neurons.Annales de l’Institut Henri Poincar´ e C, 2025
Pierre-Emmanuel Jabin and Datong Zhou. The mean-field limit of sparse networks of integrate- and-fire neurons.Annales de l’Institut Henri Poincar´ e C, 2025
2025
-
[20]
Random Batch Methods (RBM) for interacting particle systems.Journal of Computational Physics, 400:108877, 2020
Shi Jin, Lei Li, and Jian-Guo Liu. Random Batch Methods (RBM) for interacting particle systems.Journal of Computational Physics, 400:108877, 2020
2020
-
[21]
Strong solutions of stochastic equations with singular time dependent drift.Probability theory and related fields, 131(2):154–196, 2005
Nicolai V Krylov and Michael R¨ ockner. Strong solutions of stochastic equations with singular time dependent drift.Probability theory and related fields, 131(2):154–196, 2005. 47
2005
-
[22]
Springer, 2003
Harold J Kushner and G George Yin.Stochastic approximation and recursive algorithms and applications. Springer, 2003
2003
-
[23]
Hierarchies, entropy, and quantitative propagation of chaos for mean field diffusions.Probability and mathematical physics, 4(2):377–432, 2023
Daniel Lacker. Hierarchies, entropy, and quantitative propagation of chaos for mean field diffusions.Probability and mathematical physics, 4(2):377–432, 2023
2023
-
[24]
Daniel Lacker, Lane Chun Yeung, and Fuzhong Zhou. Quantitative propagation of chaos for non-exchangeable diffusions via first-passage percolation.arXiv preprint arXiv:2409.08882, 2024
Pith/arXiv arXiv 2024
-
[25]
Henry P McKean. Propagation of chaos for a class of non-linear parabolic equations.Stochastic Differential Equations (Lecture Series in Differential Equations, Session 7, Catholic Univ., 1967), pages 41–57, 1967
1967
-
[26]
A stochastic approximation method.The annals of mathematical statistics, pages 400–407, 1951
Herbert Robbins and Sutton Monro. A stochastic approximation method.The annals of mathematical statistics, pages 400–407, 1951
1951
-
[27]
Mean field limit for coulomb-type flows.Duke Mathematical Journal, 169(15):2887–2935, 2020
Sylvia Serfaty. Mean field limit for coulomb-type flows.Duke Mathematical Journal, 169(15):2887–2935, 2020
2020
-
[28]
Topics in propagation of chaos
Alain-Sol Sznitman. Topics in propagation of chaos. InEcole d’´ et´ e de probabilit´ es de Saint- Flour XIX—1989, pages 165–251. Springer, 1991
1989
-
[29]
Theory of function spaces
Hans Triebel. Theory of function spaces. iii,.Monographs in Mathematics, 2006
2006
-
[30]
Zhenfu Wang, Xianliang Zhao, and Rongchan Zhu. Mean-field limit of non-exchangeable interacting diffusions with singular kernels.To appear in SIAM J. Math. Analysis. arXiv preprint arXiv:2209.14002, 2022
arXiv 2022
-
[31]
Gaussian fluctuations for interacting par- ticle systems with singular kernels.Archive for Rational Mechanics and Analysis, 247(5):101, 2023
Zhenfu Wang, Xianliang Zhao, and Rongchan Zhu. Gaussian fluctuations for interacting par- ticle systems with singular kernels.Archive for Rational Mechanics and Analysis, 247(5):101, 2023
2023
-
[32]
Datong Zhou. Non-exchangeable mean-field theory for adaptive weights: propagation of chaos and graphon sampling lemma.arXiv preprint arXiv:2506.13587, 2025. 48
arXiv 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.