Functional Breuer–Major Theorem

• The functional Breuer–Major theorem is a path-space extension of the classical Gaussian central limit theorem, defined via the Hermite rank of nonlinear functionals of stationary processes.
• It employs Malliavin calculus and tightness criteria to rigorously establish convergence in function spaces, including Hilbert-valued and rough-path topologies.
• Recent advances provide quantitative convergence rates and broaden the theorem’s applicability to non-Gaussian, Poisson, and infinite-dimensional settings.

Searching arXiv for recent and foundational papers on the functional Breuer–Major theorem and related quantitative/Hilbert-valued extensions. The functional Breuer–Major theorem is the path-space extension of the Breuer–Major central limit theorem for nonlinear functionals of stationary Gaussian inputs. In its discrete Gaussian form, one considers normalized partial sums

$$Y_n(t)=\frac{1}{\sqrt n}\sum_{i=0}^{\lfloor nt\rfloor-1}\varphi(X_i), \qquad t\in[0,1],$$

where $X=\{X_n\}_{n\in\mathbb Z}$ is a zero-mean stationary Gaussian sequence with covariance $\rho$, and $\varphi$ has Hermite rank $d$. Under the short-memory condition $\sum_{r\in\mathbb Z}|\rho(r)|^d<\infty$, the finite-dimensional distributions converge to those of a scaled Brownian motion, and, under an additional moment assumption $E[|\varphi(X_0)|^p]<\infty$ for some $p>2$, one obtains weak convergence in path space. Modern formulations treat convergence in $D([0,1])$, $C([0,\infty))$, Hilbert-valued path spaces, Hölder topologies, and rough-path topologies, and also provide quantitative rates via Malliavin–Stein, Dirichlet-structure, and second-order Gaussian Poincaré methods (Nourdin et al., 2018, Bourguin et al., 2019).

1. Discrete Gaussian formulation

Let $X=\{X_n\}_{n\in\mathbb Z}$0 be a zero-mean stationary Gaussian sequence with

$X=\{X_n\}_{n\in\mathbb Z}$1

and let $X=\{X_n\}_{n\in\mathbb Z}$2 belong to $X=\{X_n\}_{n\in\mathbb Z}$3, where $X=\{X_n\}_{n\in\mathbb Z}$4 is the standard Gaussian law. Writing the Hermite expansion

$X=\{X_n\}_{n\in\mathbb Z}$5

the integer $X=\{X_n\}_{n\in\mathbb Z}$6 is the Hermite rank, meaning $X=\{X_n\}_{n\in\mathbb Z}$7 and $X=\{X_n\}_{n\in\mathbb Z}$8. The normalized partial-sum process is

$X=\{X_n\}_{n\in\mathbb Z}$9

The classical Breuer–Major theorem yields convergence of the finite-dimensional distributions of $\rho$0 to those of $\rho$1, where $\rho$2 is a standard Brownian motion, under the summability condition

$\rho$3

The functional theorem strengthens this by establishing convergence in $\rho$4 endowed with the Skorohod topology, provided one assumes in addition that $\rho$5 for some $\rho$6 (Nourdin et al., 2018).

A Hilbert-valued formulation replaces path-space arguments by embedding the process into a separable Hilbert space $\rho$7 containing $\rho$8 or $\rho$9 isometrically, for example $\varphi$0. In that setting one views $\varphi$1 as a $\varphi$2-valued random variable and obtains

$\varphi$3

in distribution in $\varphi$4, hence in $\varphi$5, again under $\varphi$6 for Hermite rank $\varphi$7 (Bourguin et al., 2019).

2. Variance structure and the role of Hermite rank

The limiting variance is given by the Breuer–Major formula

$\varphi$8

Under $\varphi$9, this series converges because $d$0 (Nourdin et al., 2018).

The Hermite rank is the decisive threshold in the dependence condition. It is the smallest chaos order that contributes nontrivially, and therefore the minimal power of the covariance that must be summable. In the continuous stationary version,

$d$1

the corresponding condition is

$d$2

and the limiting variance becomes

$d$3

In that formulation the condition $d$4 is both necessary and sufficient for nondegeneracy of $d$5 (Campese et al., 2018).

The same structural principle persists in examples. For fractional Gaussian noise $d$6 of Hurst index $d$7, the functional theorem applies when $d$8. For $d$9, one enters the non-Gaussian Hermite regime of Taqqu, so the Brownian conclusion is no longer the relevant limit statement (Nourdin et al., 2018).

3. Tightness, Malliavin calculus, and proof architecture

A central difficulty is that finite-dimensional convergence does not imply convergence in $\sum_{r\in\mathbb Z}|\rho(r)|^d<\infty$0 or $\sum_{r\in\mathbb Z}|\rho(r)|^d<\infty$1. The functional Breuer–Major theorem is therefore a tightness theorem as much as a central limit theorem. In the discrete Gaussian case, the proof proceeds by combining the classical Breuer–Major finite-dimensional convergence with a Billingsley-type tightness criterion. It suffices to prove that for some $\sum_{r\in\mathbb Z}|\rho(r)|^d<\infty$2,

$\sum_{r\in\mathbb Z}|\rho(r)|^d<\infty$3

This yields the required two-parameter moment condition and hence tightness in $\sum_{r\in\mathbb Z}|\rho(r)|^d<\infty$4 (Nourdin et al., 2018).

The key analytical input is a Malliavin representation of each summand as an iterated divergence. If $\sum_{r\in\mathbb Z}|\rho(r)|^d<\infty$5 is the Hermite rank and $\sum_{r\in\mathbb Z}|\rho(r)|^d<\infty$6 is the isonormal kernel of $\sum_{r\in\mathbb Z}|\rho(r)|^d<\infty$7, then

$\sum_{r\in\mathbb Z}|\rho(r)|^d<\infty$8

where

$\sum_{r\in\mathbb Z}|\rho(r)|^d<\infty$9

Meyer inequalities imply boundedness of $E[|\varphi(X_0)|^p]<\infty$0 from a Sobolev-type norm into $E[|\varphi(X_0)|^p]<\infty$1, and hypercontractivity yields uniform $E[|\varphi(X_0)|^p]<\infty$2-bounds on the relevant chaos components. Together with stationarity and $E[|\varphi(X_0)|^p]<\infty$3, one recovers the increment estimate needed for tightness (Nourdin et al., 2018).

The same scheme appears in continuous and vector-valued settings. In the continuous stationary theorem one writes

$E[|\varphi(X_0)|^p]<\infty$4

and then applies Meyer inequalities and Kolmogorov–Chentsov to obtain tightness in $E[|\varphi(X_0)|^p]<\infty$5 (Campese et al., 2018). For Gaussian vector fields $E[|\varphi(X_0)|^p]<\infty$6, one proves a $E[|\varphi(X_0)|^p]<\infty$7-fold divergence representation

$E[|\varphi(X_0)|^p]<\infty$8

with $E[|\varphi(X_0)|^p]<\infty$9 involving shift operators $p>2$0, and again derives

$p>2$1

which yields tightness in $p>2$2 (Nualart et al., 2019).

4. Quantitative functional CLTs and Hilbert-valued formulations

A major development is the transition from qualitative weak convergence to explicit rates in functional metrics. Bourguin–Campese formulate the theorem in a separable Hilbert space $p>2$3 and work with the distance

$p>2$4

where $p>2$5 consists of twice-Fréchet differentiable test functions with

$p>2$6

This metric metrizes weak convergence in $p>2$7 (Bourguin et al., 2019).

Their framework combines an infinite-dimensional version of Stein’s method, due to Shih, with a diffusive Dirichlet structure and $p>2$8-calculus. For a centered Gaussian $p>2$9 on $D([0,1])$0 with covariance $D([0,1])$1, and a $D([0,1])$2-valued random variable $D([0,1])$3, the carré-du-champ bound takes the form

$D([0,1])$4

Expanding

$D([0,1])$5

into Wiener–Itô chaos then leads to fourth-moment and contraction estimates in terms of $D([0,1])$6 (Bourguin et al., 2019).

Applied to the partial-sum process

$D([0,1])$7

this yields a fully quantitative functional Breuer–Major theorem. If

$D([0,1])$8

with $D([0,1])$9 slowly varying at $C([0,\infty))$0, then there exists $C([0,\infty))$1 such that

$C([0,\infty))$2

with more refined cases when $C([0,\infty))$3 (Bourguin et al., 2019).

For $C([0,\infty))$4 and $C([0,\infty))$5 with $C([0,\infty))$6, this recovers the rates of Nourdin-Peccati–Podolskij:

• if $C([0,\infty))$7, then $C([0,\infty))$8;
• if $C([0,\infty))$9, then $X=\{X_n\}_{n\in\mathbb Z}$00;
• if $X=\{X_n\}_{n\in\mathbb Z}$01, then $X=\{X_n\}_{n\in\mathbb Z}$02 (Bourguin et al., 2019).

Related quantitative work includes contraction-based bounds for the one-dimensional Breuer–Major CLT on Wiener space, where the rate depends on the Hermite rank and the chaotic gap, and convex-distance bounds for the finite-dimensional distributions of the functional theorem (Fissler et al., 2016, Nourdin et al., 2020).

5. Extensions across settings and topologies

The functional Breuer–Major paradigm now spans several ambient spaces and noise structures. In the continuous stationary setting,

$X=\{X_n\}_{n\in\mathbb Z}$03

converges in $X=\{X_n\}_{n\in\mathbb Z}$04 under $X=\{X_n\}_{n\in\mathbb Z}$05 and $X=\{X_n\}_{n\in\mathbb Z}$06 for some $X=\{X_n\}_{n\in\mathbb Z}$07 (Campese et al., 2018). For Gaussian vector fields $X=\{X_n\}_{n\in\mathbb Z}$08, one obtains

$X=\{X_n\}_{n\in\mathbb Z}$09

assuming $X=\{X_n\}_{n\in\mathbb Z}$10 for some $X=\{X_n\}_{n\in\mathbb Z}$11 and $X=\{X_n\}_{n\in\mathbb Z}$12 (Nualart et al., 2019).

A genuinely infinite-dimensional version appears in the Hilbert space-valued theorem for stationary Gaussian processes $X=\{X_n\}_{n\in\mathbb Z}$13 with values in $X=\{X_n\}_{n\in\mathbb Z}$14 and measurable $X=\{X_n\}_{n\in\mathbb Z}$15. If $X=\{X_n\}_{n\in\mathbb Z}$16 and

$X=\{X_n\}_{n\in\mathbb Z}$17

then

$X=\{X_n\}_{n\in\mathbb Z}$18

and the continuous-time version

$X=\{X_n\}_{n\in\mathbb Z}$19

converges in $X=\{X_n\}_{n\in\mathbb Z}$20 to $X=\{X_n\}_{n\in\mathbb Z}$21 (Düker et al., 2024).

A second-order Gaussian Poincaré approach gives another Hilbert-valued functional theorem. With

$X=\{X_n\}_{n\in\mathbb Z}$22

and an envelope condition on the kernels $X=\{X_n\}_{n\in\mathbb Z}$23, one has

$X=\{X_n\}_{n\in\mathbb Z}$24

where $X=\{X_n\}_{n\in\mathbb Z}$25 is a centered Gaussian random element in $X=\{X_n\}_{n\in\mathbb Z}$26 with the same covariance operator as $X=\{X_n\}_{n\in\mathbb Z}$27. Moreover $X=\{X_n\}_{n\in\mathbb Z}$28 converges in $X=\{X_n\}_{n\in\mathbb Z}$29 to $X=\{X_n\}_{n\in\mathbb Z}$30, so $X=\{X_n\}_{n\in\mathbb Z}$31 (Vidotto et al., 16 Jun 2025).

The theorem also extends beyond Gaussian input. In the Poisson setting, if $X=\{X_n\}_{n\in\mathbb Z}$32 is a Poisson point process on $X=\{X_n\}_{n\in\mathbb Z}$33, $X=\{X_n\}_{n\in\mathbb Z}$34, $X=\{X_n\}_{n\in\mathbb Z}$35 with $X=\{X_n\}_{n\in\mathbb Z}$36, and $X=\{X_n\}_{n\in\mathbb Z}$37 with $X=\{X_n\}_{n\in\mathbb Z}$38, then the $X=\{X_n\}_{n\in\mathbb Z}$39-th-and-higher-order chaos truncation $X=\{X_n\}_{n\in\mathbb Z}$40 converges in law in $X=\{X_n\}_{n\in\mathbb Z}$41 to $X=\{X_n\}_{n\in\mathbb Z}$42. Tightness is proved by an $X=\{X_n\}_{n\in\mathbb Z}$43-spectral-gap inequality for Poisson functionals (Kong et al., 30 Oct 2025).

Recent work strengthens the topology. In the summable-covariance regime, interpolated partial sums

$X=\{X_n\}_{n\in\mathbb Z}$44

converge to $X=\{X_n\}_{n\in\mathbb Z}$45 in $X=\{X_n\}_{n\in\mathbb Z}$46 for every $X=\{X_n\}_{n\in\mathbb Z}$47, assuming $X=\{X_n\}_{n\in\mathbb Z}$48 for some $X=\{X_n\}_{n\in\mathbb Z}$49 and $X=\{X_n\}_{n\in\mathbb Z}$50 (Benning et al., 2 Jun 2026). A rough-path version lifts the theorem to the càdlàg rough-path space $X=\{X_n\}_{n\in\mathbb Z}$51, with a Brownian rough-path limit carrying explicit symmetric and antisymmetric area corrections (Altman et al., 18 Feb 2026).

6. Optimality, misconceptions, and non-Gaussian boundaries

A persistent misconception is that the classical Breuer–Major theorem already implies path-space convergence. It does not. The classical theorem gives finite-dimensional convergence under $X=\{X_n\}_{n\in\mathbb Z}$52-integrability and covariance summability, but functional convergence requires tightness, and tightness cannot in general be obtained from mere $X=\{X_n\}_{n\in\mathbb Z}$53-bounds because one needs control of higher moments of increments (Nourdin et al., 2018, Campese et al., 2018).

The additional assumption $X=\{X_n\}_{n\in\mathbb Z}$54 or $X=\{X_n\}_{n\in\mathbb Z}$55 for some $X=\{X_n\}_{n\in\mathbb Z}$56 is therefore not a cosmetic strengthening. In the discrete functional theorem it is described as a “sufficient (and almost necessary) natural condition,” and in the detailed exposition it is said to be minimal in the sense that it exactly produces the moment power $X=\{X_n\}_{n\in\mathbb Z}$57 needed in Billingsley’s criterion (Nourdin et al., 2018). Earlier criteria, such as those of Ben Hariz and Chambers–Slud, required rather unnatural summability conditions on all Hermite coefficients; the Malliavin–Meyer method reduces the problem to a single extra $X=\{X_n\}_{n\in\mathbb Z}$58-assumption (Nourdin et al., 2018).

Another conceptual boundary concerns the dependence regime. The Brownian limit is tied to the summable-covariance regime $X=\{X_n\}_{n\in\mathbb Z}$59 or $X=\{X_n\}_{n\in\mathbb Z}$60. When this fails, the limit need not remain Gaussian. In the continuous non-stationary self-similar setting, the non-central case $X=\{X_n\}_{n\in\mathbb Z}$61 leads to a Hermite process of order $X=\{X_n\}_{n\in\mathbb Z}$62; in the critical case $X=\{X_n\}_{n\in\mathbb Z}$63, the Gaussian limit survives only after logarithmic renormalization (Campese et al., 2018). Similarly, for fractional Gaussian noise with $X=\{X_n\}_{n\in\mathbb Z}$64, one leaves the Brownian regime and enters the non-Gaussian Hermite regime (Nourdin et al., 2018).

Finally, the modern literature clarifies that the functional theorem is not tied to the Skorohod topology. Hölder convergence is strictly stronger than Skorohod $X=\{X_n\}_{n\in\mathbb Z}$65 convergence, and rough Hölder or rough-path convergence is stronger still. This matters because many operations on paths, including Young integration and solution maps of differential equations, are continuous in Hölder or rough Hölder topologies but not in Skorohod topology (Benning et al., 2 Jun 2026). The functional Breuer–Major theorem has thus evolved from a pathwise strengthening of a classical CLT into a family of invariance principles across Gaussian, Poisson, Hilbert-valued, and rough-path frameworks.