Before and beyond the mixing time: New approximations for additive functionals of stationary Gauss-Markov processes
Pith reviewed 2026-05-12 04:20 UTC · model grok-4.3
The pith
Additive functionals of stationary Gauss-Markov processes exhibit a phase transition in their scaling limits when mixing time is order n.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Quantitative and functional limit theorems hold for additive functionals along triangular arrays of stationary Gaussian Markov processes when t_mix scales sublinearly, linearly, or superlinearly with n. A phase transition occurs at t_mix ≍ n, and the emerging new limit processes are identified along with their interrelation properties in the regimes at and before the mixing time.
What carries the argument
Triangular arrays of stationary Gaussian Markov processes whose mixing times scale relative to the number of observations n; the phase transition this induces in the scaling limits of their additive functionals.
If this is right
- Standard Brownian or other ergodic limits apply when t_mix is much smaller than n.
- Distinct new limit processes govern the additive functionals exactly when t_mix is order n.
- The limits obtained before the mixing time are interrelated with the critical-regime limits.
- The theorems supply approximations that remain valid uniformly across all three scaling regimes.
Where Pith is reading between the lines
- The phase-transition structure could guide adaptive block-length choices in resampling procedures for strongly dependent time series.
- Similar triangular-array constructions may allow extension of the limit identification to non-Gaussian Markov chains.
- In practice, an estimate of t_mix from data could select which limiting approximation to use for inference.
Load-bearing premise
The processes are stationary Gaussian Markov chains arranged in triangular arrays, with mixing times scaling in the stated sub-, super-, and proportionate regimes relative to n.
What would settle it
Simulate sequences with mixing time set to n/2, n, and 2n; the normalized additive functional should converge in law to three distinct objects with the middle regime distinct from the other two. Absence of this change in the empirical distributions would refute the phase transition.
Figures
read the original abstract
Whereas classical invariance principles for ergodic Markov chains address the situation in which the time horizon of observations is much larger than the mixing time, the quality of approximation is questionable when this is not the case anymore -- even when starting the Markov chain in the invariant law. In this article, we prove quantitative and functional limit theorems for additive functionals along triangular arrays of stationary Gaussian Markov processes when the mixing time $t_{\text{mix}}$ scales sub-, super- and proportionately to the number of observations $n$. Our major finding is a phase-transition at $t_{\text{mix}}\asymp n$, together with the identification and interrelation properties of the emerging new limit processes at and before the mixing time.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proves quantitative and functional limit theorems for additive functionals of stationary Gaussian Markov processes in triangular arrays, considering regimes where the mixing time t_mix scales sublinearly, superlinearly, or linearly with the number of observations n. The central result identifies a phase transition at t_mix ≍ n and characterizes the new limiting processes (including their interrelations) that emerge at and before the mixing time, leveraging the explicit exponential covariance structure for direct computations of finite-dimensional distributions and tightness.
Significance. If the derivations hold, the work fills a gap in invariance principles for Markov processes by providing approximations valid when the observation horizon is comparable to or smaller than the mixing time, extending beyond classical ergodic settings. The explicit covariance calculations enable parameter-free derivations of the limits in each regime and yield falsifiable predictions for the interrelation properties, which strengthens the contribution for applications in time-series analysis and stochastic modeling.
major comments (2)
- §3.2, Theorem 3.3 (proportionate regime): the tightness argument for the functional convergence uses moment bounds derived from the covariance kernel; however, the uniformity of these bounds as t_mix/n → c > 0 is not explicitly verified for all c in [0,∞), which is load-bearing for the phase-transition claim across the full range of scalings.
- §4.1, Eq. (4.5): the identification of the pre-mixing limit process as a specific Gaussian process with covariance min(s,t) exp(-|s-t|/α) assumes the triangular-array stationarity holds uniformly in the scaling parameter; a counter-example check or explicit verification for the boundary case t_mix/n → 0 would strengthen the interrelation properties stated in Theorem 4.4.
minor comments (3)
- Notation for the mixing time t_mix is introduced in the abstract but first defined only in §2.1; an early global definition would improve readability.
- Figure 1 (schematic of the three regimes) has axis labels that are too small for print; enlarging them and adding a legend for the limit processes would aid clarity.
- The proof of the quantitative bound in Proposition 2.4 cites a general result from [reference] but does not spell out the constant dependence on the Gaussian variance parameter; adding one line of explicit computation would help.
Simulated Author's Rebuttal
We thank the referee for the positive evaluation and the insightful comments, which have helped clarify the presentation of our results on the phase transition for additive functionals of stationary Gauss-Markov processes. We address the major comments point by point below.
read point-by-point responses
-
Referee: §3.2, Theorem 3.3 (proportionate regime): the tightness argument for the functional convergence uses moment bounds derived from the covariance kernel; however, the uniformity of these bounds as t_mix/n → c > 0 is not explicitly verified for all c in [0,∞), which is load-bearing for the phase-transition claim across the full range of scalings.
Authors: We thank the referee for highlighting this point. The moment bounds used in the tightness argument are derived directly from the explicit exponential covariance structure of the Gauss-Markov process. These bounds depend continuously on the scaling parameter c = t_mix/n, and are uniform on any compact interval [0, C]. For large c, the rapid decay of correlations reduces the process to one with nearly independent increments, for which standard tightness criteria apply with constants independent of c. In the revised manuscript, we have added a short lemma (new Lemma 3.4) that makes this uniformity explicit for all c ∈ [0, ∞), thereby supporting the phase-transition statement across all regimes. revision: yes
-
Referee: §4.1, Eq. (4.5): the identification of the pre-mixing limit process as a specific Gaussian process with covariance min(s,t) exp(-|s-t|/α) assumes the triangular-array stationarity holds uniformly in the scaling parameter; a counter-example check or explicit verification for the boundary case t_mix/n → 0 would strengthen the interrelation properties stated in Theorem 4.4.
Authors: We agree that an explicit boundary check strengthens the interrelation claims. By construction, each process in the triangular array is started in its invariant distribution, so stationarity holds for every n and every scaling. As t_mix/n → 0 the covariance kernel min(s,t) exp(-|s-t|/(t_mix/n)) converges pointwise to min(s,t), and the finite-dimensional distributions converge to those of standard Brownian motion. This is consistent with the classical invariance principle under fast mixing. In the revision we have inserted a remark following Theorem 4.4 that records this limit and the continuity argument, confirming the interrelation properties without requiring a separate counter-example search. revision: yes
Circularity Check
No significant circularity; derivation self-contained from covariance structure
full rationale
The paper establishes quantitative and functional limit theorems for additive functionals of stationary Gaussian Markov processes in triangular arrays by direct computation of finite-dimensional distributions and tightness criteria from the explicit exponential covariance structure. The phase transition at t_mix ≍ n and identification of emerging limit processes follow from these calculations in the sub-, super-, and proportionate mixing regimes without any reduction to fitted parameters, self-definitional constructions, or load-bearing self-citations. The central claims are independent of the inputs and rest on standard probabilistic tools applied to the given process assumptions.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Processes are stationary Gaussian Markov chains
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Our major finding is a phase-transition at t_mix ≍ n, together with the identification and interrelation properties of the emerging new limit processes at and before the mixing time.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
W_{q,γ}(t) = K_q^{(γ)} ∫_{R^q} … (spectral / time-domain multiple Wiener-Itô integral)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Bai, S. and Taqqu, M. S. (2013). Multivariate limit theorems in the context of long-range depen- dence.J. Time Series Anal., 34(6):717–743. Bhattacharya, R. N. and Lee, O. (1988). Asymptotics of a class of markov processes which are not in general irreducible.The Annals of Probability, 16(3):1333 –
work page 2013
-
[2]
Bhattacharyya, B. B., Richardson, G. D., and Franklin, L. A. (1997). Asymptotic inference for near unit roots in spatial autoregression.The Annals of Statistics, 25(4):1709 –
work page 1997
-
[3]
(1999).Convergence of probability measures
Billingsley, P. (1999).Convergence of probability measures. Wiley Series in Probability and Statis- tics: Probability and Statistics. John Wiley & Sons, Inc., New York, second edition. A Wiley- Interscience Publication. Buchmann, B. and Chan, N. H. (2007). Asymptotic theory of least squares estimators for nearly unstable processes under strong dependence....
work page 1999
-
[4]
Chan, N. H. and Wei, C. Z. (1987). Asymptotic inference for nearly nonstationary AR(1) processes. The Annals of Statistics, 15(3):1050 –
work page 1987
-
[5]
Chen, X. (1999). How often does a Harris recurrent Markov chain recur?The Annals of Probability, 27(3):1324 –
work page 1999
-
[6]
Chen, X., Ma, Z.-M., and Peng, X. (2019). Distribution flows associated with positivity preserving coercive forms.The Annals of Probability, 47(5):2894 –
work page 2019
-
[7]
J., Kuwae, K., and Zhang, T.-S
Chen, Z.-Q., Fitzsimmons, P. J., Kuwae, K., and Zhang, T.-S. (2008). Stochastic calculus for symmetric Markov processes.The Annals of Probability, 36(3):931 –
work page 2008
-
[8]
Cox, D. D. (1991). Gaussian likelihood estimation for nearly nonstationary AR(1) processes.The Annals of Statistics, 19(3):1129 –
work page 1991
-
[9]
Cox, D. D. and Llatas, I. (1991). Maximum likelihood type estimation for nearly nonstationary autoregressive time series.The Annals of Statistics, 19(3):1109 –
work page 1991
-
[10]
Davis, R. A. (1983). Stable limits for partial sums of dependent random variables.The Annals of Probability, 11(2):262 –
work page 1983
-
[11]
Dehling, H., Denker, M., and Philipp, W. (1986). Central limit theorems for mixing sequences of random variables under minimal conditions.The Annals of Probability, 14(4):1359 –
work page 1986
-
[12]
Deuschel, J.-D., Orenshtein, T., and Perkowski, N. (2021). Additive functionals as rough paths. The Annals of Probability, 49(3):1450 –
work page 2021
-
[13]
Dobrushin, R. L. and Major, P. (1979). Non-central limit theorems for non-linear functionals of Gaussian fields.Zeitschrift f ¨ur Wahrscheinlichkeitstheorie und Verwandte Gebiete, 50:27 –
work page 1979
-
[14]
Durrett, R. and Resnick, S. I. (1978). Functional limit theorems for dependent variables.The Annals of Probability, 6(5):829 –
work page 1978
-
[15]
Dynkin, E. B. (1977). Markov systems and their additive functionals.The Annals of Probability, 5(5):653 –
work page 1977
-
[16]
Dynkin, E. B. (1991). Branching particle systems and superprocesses.The Annals of Probability, 19(3):1157 –
work page 1991
-
[17]
29 Getoor, R. K. and Sharpe, M. J. (1973). Last exit times and additive functionals.The Annals of Probability, 1(4):550 –
work page 1973
-
[18]
Glover, J. (1980). Compactifications for dual processes.The Annals of Probability, 8(6):1119 –
work page 1980
-
[19]
Hairer, M. (2021). Introduction to Malliavin calculus. Lecture notes, available athttps://www. hairer.org/notes/Malliavin.pdf. Hardy, G. H., Littlewood, J. E., and P´ olya, G. (1988).Inequalities. Cambridge Mathematical Library. Cambridge University Press, Cambridge. Reprint of the 1952 edition. Hong, M., Liu, H., and Xu, F. (2024). Limit theorems for add...
work page 2021
-
[20]
Jain, N. C. (1990). Large deviation lower bounds for additive functionals of Markov processes. The Annals of Probability, 18(3):1071 –
work page 1990
-
[21]
Jara, M., Komorowski, T., and Olla, S. (2009). Limit theorems for additive functionals of a Markov chain.The Annals of Applied Probability, 19(6):2270 –
work page 2009
-
[22]
(2021).Foundations of modern probability, volume 99 ofProbability Theory and Stochastic Modelling
Kallenberg, O. (2021).Foundations of modern probability, volume 99 ofProbability Theory and Stochastic Modelling. Springer, Cham, third edition. Kaspi, H. (1988). Random time changes for processes with random birth and death.The Annals of Probability, 16(2):586 –
work page 2021
-
[23]
(2006).Introduction to stochastic integration
Kuo, H.-H. (2006).Introduction to stochastic integration. Universitext. Springer, New York. Kuwae, K. (2010). Stochastic calculus over symmetric Markov processes without time reversal. The Annals of Probability, 38(4):1532 –
work page 2006
-
[24]
Maisonneuve, B. (1975). Exit systems.The Annals of Probability, 3(3):399 –
work page 1975
-
[25]
(2014).Multiple Wiener-Itˆ o integrals, volume 849 ofLecture Notes in Mathematics
Major, P. (2014).Multiple Wiener-Itˆ o integrals, volume 849 ofLecture Notes in Mathematics. Springer, Cham, second edition. With applications to limit theorems. Marcus, M. B. and Rosen, J. (1996). Gaussian chaos and sample path properties of additive functionals of symmetric Markov processes.The Annals of Probability, 24(3):1130 –
work page 2014
-
[26]
Maxwell, M. and Woodroofe, M. (2000). Central limit theorems for additive functionals of Markov chains.The Annals of Probability, 28(2):713 –
work page 2000
-
[27]
Nelson, E. (1973). The free Markoff field.Journal of Functional Analysis, 12(2):211 –
work page 1973
-
[28]
Nevison, C. H. (1976). Potentials of Markov processes without duality.The Annals of Probability, 4(3):497 –
work page 1976
-
[29]
Ney, P. and de Acosta, A. (1998). Large deviation lower bounds for arbitrary additive functionals of a Markov chain.The Annals of Probability, 26(4):1660 –
work page 1998
-
[30]
Nourdin, I. (2012). Lectures on Gaussian approximations with Malliavin calculus. 30 Nourdin, I. and Peccati, G. (2009). Stein’s method on Wiener chaos.Probability Theory and Related Fields, 145(1):75 –
work page 2012
-
[31]
Nourdin, I., Peccati, G., and Reinert, G. (2010). Invariance principles for homogeneous sums: universality of Gaussian Wiener chaos.The Annals of Probability, 38(5):1947 –
work page 2010
-
[32]
Nourdin, I. and Poly, G. (2013). Convergence in total variation on Wiener chaos.Stochastic Processes and their Applications, 123(2):651–674. Nualart, D. and Peccati, G. (2005). Central limit theorems for sequences of multiple stochastic integrals.The Annals of Probability, 33(1):177 –
work page 2013
-
[33]
Pacheco, A. and Prabhu, N. U. (1996). A Markovian storage model.The Annals of Applied Probability, 6(1):76 –
work page 1996
-
[34]
Peccati, G. and Taqqu, M. S. (2011).Wiener chaos: moments, cumulants and diagrams, volume 1 ofBocconi & Springer Series. Springer, Milan; Bocconi University Press, Milan. A survey with computer implementation, Supplementary material available online. Peccati, G. and Tudor, C. A. (2005). Gaussian limits for vector-valued multiple stochastic integrals. InS´...
work page 2011
-
[35]
Peligrad, M., Utev, S., and Wu, W. B. (2007). A maximalL p-inequality for stationary sequences and its applications.Proceedings of the American Mathematical Society, 135(2):541 –
work page 2007
-
[36]
Pipiras, V. and Taqqu, M. S. (2010). Regularization and integral representations of Hermite processes.Statistics & Probability Letters, 80(23–24):2014 –
work page 2010
-
[37]
Sabzikar, F. (2015). Tempered Hermite process.Modern Stochastics: Theory and Applications, 2(4):327 –
work page 2015
-
[38]
Tsybakov, A. B. (2009).Introduction to nonparametric estimation. Springer Series in Statis- tics. Springer, New York. Revised and extended from the 2004 French original, Translated by Vladimir Zaiats. 31 Appendix A Discrete-time stationary Gauss-Markov process The following lemma is well-known, yet we did not find an appropriate discrete-time refer- ence ...
work page 2009
-
[39]
Ifα= 1, the process is constant in time whereasα=−1 contradicts stationarity
= Var(X 0). Ifα= 1, the process is constant in time whereasα=−1 contradicts stationarity. Thus, it remains to study|α|<1. Proposi- tion 14.7 in Kallenberg (2021) reveals Cov(X0, Xk+1)Var(X1) = Cov(X0, X1)Cov(X1, Xk+1). By stationarity, Var(X1) = Var(X0) and Cov(X1, Xk+1) = Cov(X0, Xk), hence Cov(X0, Xk+1) = Cov(X0, X1) Var(X0) Cov(X0, Xk) =αCov(X 0, Xk).(...
work page 2021
-
[40]
As concernsl < k, we writeε l =X l −µ−α(X l−1 −µ)
= Var(X0) +α 2Var(X0)−2α 2Var(X0) = Var(X0)(1−α 2). As concernsl < k, we writeε l =X l −µ−α(X l−1 −µ). Thenε l is represented as a function ofX l andX l−1, both in the past of time-stepk. Due to previous observations, εk ⊥ ⊥εl and more generallyε k ⊥ ⊥(ε1, . . . , εk−1). By induction, the sequence (ε k)k∈N is i.i.d. Finally, isolatingX k in the defining e...
work page 2009
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.