arXiv paper· probability theory (applied probability / operations research)· method paper (pure theory)

Gaussian Approximation for Multivariate Martingale Sums from Uniformly Ergodic Markov Chains

Yixuan Zhang, Qiaomin Xie

math.PR · cs.LG · stat.ML

Quick verdict

Claims first optimal O(n^-1/2) W_p bounds for Markov-chain martingale sums; proof looks coherent though dense, no numerical verification included.

Summary

The paper proves nonasymptotic Gaussian approximation bounds in Wasserstein distance W_p (p ≥ 2) for sums of multivariate martingale differences generated by a uniformly ergodic Markov chain. Starting from the Poisson-equation martingale decomposition of additive functionals, the authors obtain, under an L^{(2+η)p} moment condition (η > 0) and covariance normalization, the explicit bound C(η,C0,ρ0)·(p³‖A‖₄² + p·d^{1/4}‖A‖₂^{1/2}‖A‖₄²), where A collects the L^{(2+η)p}-sizes of the weighted increments. In the balanced-increment regime (increments of order n^{-1/2}) this yields the first optimal O(n^{-1/2}) W_p rate for fixed p and d, and Corollary 1 transfers it to additive functionals of Markov chains, improving the O(n^{-1/4}) available under weaker geometric-ergodicity assumptions. Methodologically the paper contributes (i) a tensor-level restatement of the Ornstein–Uhlenbeck relative-score argument of Fang and Koike (2023) for antisymmetric Stein couplings, which avoids 4p-moment requirements, and (ii) a refresh-then-maximal coupling combining an independent first transition with symmetric blockwise maximal coupling, yielding the exact Stein identity, conditional exchangeability, geometric meeting times, and prefix measurability. Limitations: the p-dependence (p³) is weaker than the near-linear scaling in the independent case; constants are unspecified; there are no numerical experiments.

Review complete

Assessment at a glance

probability theory (applied probability / operations research)

3

risk flags identified

● high 0

● medium 2

● low 1

Reproducibility readinessmedium
Code & data availablenot stated
Review confidencehigh
AI text signal0/20 passages flagged
Deep research viableCode-reproducible

4 strengths · 5 weaknesses · 5 author questions

AI text signal = Binoculars zero-shot detector score over 20 sampled passages (Fast-DetectGPT co-reported) — absence does not prove human authorship, and a flag does not imply misconduct.

Levels reflect verification complexity — not a quality score. Evolve never rates your paper.

Key claims under review

  • Under uniform ergodicity, stationary initialization, L^{(2+η)p} moments and Cov(S_n)=I_d, W_p(L(S_n), N(0,I_d)) ≤ C(η,C0,ρ0)·(p³‖A‖₄² + p·d^{1/4}‖A‖₂^{1/2}‖A‖₄²) (Theorem 1).
  • In the balanced-increment regime this is the first optimal O(n^{-1/2}) W_p rate (fixed p, d) for multivariate Markov-chain-induced martingale sums.
  • Corollary 1 gives the first optimal O(n^{-1/2}) W_p bound for multivariate additive functionals of uniformly ergodic chains under an L^{(2+η)p} condition, improving the O(n^{-1/4}) of Zhang–Xie (2026).
  • Formulating the OU relative-score bound over antisymmetric Stein couplings with conditional tensors avoids the 4p-moment requirement of applying Fang–Koike (2023, Thm 7.1) directly.

Strengths

  • Optimal rate achieved

    Theorem 1's bound O(p³‖A‖₄² + p·d^{1/4}‖A‖₂^{1/2}‖A‖₄²) yields O(n^{-1/2}) in the balanced regime; Corollary 1 gives the first optimal rate for additive functionals under uniform ergodicity, vs the O(n^{-1/4}) of Zhang–Xie (2026) under weaker geometric conditions.

  • Moment-efficient reformulation of the OU score argument

    Proposition 1 retains conditional tensors from an antisymmetric Stein coupling so contributions cancel before norms are taken; Section 4.1 argues this drops the 4p-moment requirement of a direct Fang–Koike (2023, Thm 7.1) application to (2+η)p.

  • Genuinely new coupling construction

  • Honest scoping of limitations

Weaknesses

  • No numerical or empirical sanity check

  • p-dependence far from benchmark

  • Headline moment condition stronger than needed for the n-rate

  • Novelty is partly reformulation

  • Restrictive ergodicity assumption

Full Review Report

$1.99 · one-time

Unlock the complete preliminary analysis Evolve produced for this paper.

Report preview

Scope & Claims

The paper targets quantitative Gaussian approximation in higher-order Wasserstein distance W_p (p ≥ 2) for S_n = Σ Y_i, where Y_i = M_i(g(X_i) − Pg(X_{i−1})) are martingale differences induced through the Poisson equation g − Pg = h of a stationary, uniformly ergodic Markov chain (Assumption 1: sup_x ‖P^m(x,·) − π‖_TV ≤ C0ρ0^m). Under an L^{(2+η)p} moment condition on h (q = (2+η)p) and normalization Cov(S_n) = I_d, Theorem 1 proves W_p(L(S_n), N(0,I_d)) ≤ C(η,C0,ρ0)·(p³‖A‖₄² + p·d^{1/4}‖A‖₂^{1/2}‖A‖₄²), with A = (‖Y_i‖_{L^q}). Corollary 1 (Appendix B) transfers the bound to normalized additive functionals Σh(X_i)/√n under the covariance limit (9). The headline claims are the first optimal O(n^{-1/2}) W_p rate for fixed p and d in the balanced-increment regime — for both the martingale sums and the additive functionals — improving the O(n^{-1/4}) obtainable from the authors' own prior geometric-ergodicity result (Zhang & Xie 2026).

The remaining sections cover Methodology Assessment, Statistical Evidence, Results & Reproducibility, Limitations & Risks, and Verdict & Recommendations.

  • All 4 strengths & 5 weaknesses, analyzed in detail
  • 3 risk flags with severity, section & reviewer questions
  • Methods & statistical evidence assessment
  • 5 questions a peer reviewer would ask
  • Full six-dimension evaluation report & recommendations

One-time purchase · this paper only

One-time purchase, this paper only. AI-generated preliminary analysis — not peer review.

Frequently asked questions

What is the main result of this paper?

It proves a quantitative bound on how far the sum of n multivariate martingale differences, generated by a uniformly ergodic Markov chain, is from a matching Gaussian in p-th order Wasserstein distance (any p ≥ 2). The bound is O(p³‖A‖₄² + p·d^{1/4}‖A‖₂^{1/2}‖A‖₄²), where the entries of A are the L^{(2+η)p}-sizes of the individual increments. When all increments have size about n^{-1/2} (the 'balanced' case), the distance shrinks at the optimal O(n^{-1/2}) rate for fixed p and dimension d.

How does this improve on previous Markov-chain Gaussian approximation results?

Earlier W_p (p ≥ 2) results for Markov-chain additive functionals, notably Zhang and Xie (2026), achieved only an O(n^{-1/4}) rate under weaker geometric-ergodicity assumptions; this paper reaches the optimal O(n^{-1/2}) for uniformly ergodic chains under an L^{(2+η)p} moment condition. Prior multivariate martingale work covered mainly W_1 (Srikant 2025; Wu et al. 2025), and the scalar higher-order result of Guo (2026) degrades to O(n^{-1/4}) in W₂ and O(n^{-1/6}) in W₃ in this regime, per the authors' comparison.

Is there code or numerical experiments accompanying the paper?

No. The paper is purely theoretical: it presents theorems and full proofs (Sections 4–6, Appendices A–B) with no simulation study and no software repository. A reader could nonetheless verify the predicted O(n^{-1/2}) decay on a simple finite-state Markov chain without much effort.

What is the 'refresh-then-maximal coupling'?

A two-stage coupling of two copies of the Markov chain: the first transition of the second copy is resampled independently (the 'refresh'), which makes an exact Stein identity provable, and from then on the two trajectories are tied with a symmetric blockwise maximal coupling that makes them meet in geometrically few blocks. This combination supplies both the Stein identity and the small coupling increments that the paper's score bound requires, while keeping conditional expectations measurable with respect to short prefixes of the observed chain.

About this review

This is an AI-generated preliminary review by Evolve.Pub, produced from the openly available text of the paper. It is not peer review and does not represent the authors' views. The canonical record of the paper is at the publisher link above.

  • PDF extraction jumbles sub/superscripts (e.g. '∥A∥2 4' stands for ‖A‖₄²); formulas in this review were reconstructed from context and the abstract, which matches the body.
  • The extraction runs from the title through references and Appendices A–B, ending at page 28 — no truncation detected at the far end.

Go deeper on this paper

Wrote a paper? Stress-test it before reviewers do.

Upload your manuscript for the same structured scrutiny — reasoning, statistics, evidence, figures, and reproducibility.

Scan your paper
Gaussian Approximation for Multivariate Martingale Sums from Uniformly Ergodic Markov Chains | Evolve.Pub