---
title: Hardy Martingales
url: https://www.emergentmind.com/topics/hardy-martingales
type: topic
---

# Hardy Martingales

Hardy martingales are martingales on the product torus \( \mathbb{T}^{\mathbb{N}} \) whose increments are analytic in the last variable, equivalently martingales for which each conditioned difference belongs to \( H_0^1(\mathbb{T}) \). Introduced and developed by Garling, Maurey, and Bourgain, they occupy a central position at the interface of martingale theory, Hardy space methods, Banach space geometry, and harmonic analysis; the subject includes sharp square-function estimates, analyticity-preserving Davis-type decompositions, dyadic perturbation theory, product decompositions with BMO, and vector-valued as well as noncommutative extensions [1009.3629][1504.06513][2509.07616].

## 1. Definition and analytic model

In the standard model, one works on \( \mathbb{T}^{\mathbb{N}} \) with product Haar measure \( \mathbb{P} \) and the natural filtration \( \mathcal{F}_k \) generated by the first \( k \) coordinates. A martingale \( F=(F_k) \) is a Hardy martingale if, for each \( k \), the martingale difference
\[
\Delta F_k = F_k - F_{k-1}
\]
defines, conditionally on \( \mathcal{F}_{k-1} \), a function in \( H_0^1(\mathbb{T}) \); in the formulation used for \( H^1(\mathbb{T}^{\mathbb{N}}) \), this means that
\[
t \mapsto \Delta_n f(x_1,\ldots,x_{n-1},t)
\]
belongs to \( H^1(\mathbb{T}) \) for almost every fixed history \((x_1,\ldots,x_{n-1})\) [1009.3629][2509.07616].

The same analytic constraint can be expressed spectrally. For functions in the Hardy martingale space on \( \mathbb{T}^{\mathbb{N}} \), the Fourier transform is supported in the positive cone for the lexicographic order, corresponding to analytic directions in each variable. At the \(H^1\)-level, the norm is equivalent to the square-function norm
\[
\|f\|_{H^1_{\mathrm{last}}} \simeq \left\| \left( \sum_k |\Delta_k f|^2 \right)^{1/2} \right\|_{L^1(\mathbb{T}^{\mathbb{N}})},
\]
which already indicates that analyticity and martingale square functions are inseparable in the subject [2509.07616].

Classical \(L^1\)-bounds persist in this analytic setting. For complex-valued Hardy martingales one has
\[
\| \sup_k |F_k| \|_1 \leq e \sup_k \|F_k\|_1,
\]
and the conditional square function
\[
\Bigl( \sum_k \mathbb{E}_{k-1} |\Delta G_k|^2 \Bigr)^{1/2}
\]
defines the previsible norm
\[
\|G\|_{\mathcal P}
=
\Big\|
\Bigl( \sum_k \mathbb{E}_{k-1} |\Delta G_k|^2 \Bigr)^{1/2}
\Big\|_{L^1},
\]
which plays a structural role in transform estimates, decomposition theorems, and quotient-space arguments [1611.02653].

## 2. Square functions, thin-thick splitting, and analytic Davis-Garsia theory

A foundational theorem due to Garling states that for every Hardy martingale \(F\) there are constants \(c,C>0\) such that
\[
c \,\mathbb{E}|F_n|
\leq
\mathbb{E}\left(\sum_{k=1}^n |\Delta F_k|^2\right)^{1/2}
\leq
C\,\mathbb{E}|F_n|.
\]
This square-function characterization already distinguishes Hardy martingales from general martingales by replacing a purely probabilistic control with one adapted to analytic increments [1009.3629].

A decisive refinement is the thin-thick decomposition. Every Hardy martingale \(F=(F_k)_{k=1}^n\) can be written as
\[
F=G+B,
\]
where \(G\) and \(B\) are again Hardy martingales and satisfy
\[
\mathbb{E} \left( \sum_{k = 1 }^n \mathbb{E}_{k-1} |\Delta G_k |^2 \right)^{1/2}
+
\mathbb{E} \left( \sum_{k = 1 }^n  |\Delta B_k | \right)
\leq
C_1\mathbb{E} | F_n |,
\]
together with the previsible estimate
\[
|\Delta G_k| \leq A_0 |F_{k-1}|, \qquad \forall k \leq n.
\]
The right-hand side involves only \( \mathbb{E}|F_n| \), not the square function, and the control of \( \Delta G_k \) depends only on the immediately previous value rather than on the full past maximum [1009.3629].

The analytic input for this splitting is Bourgain’s complex convexity inequality:
\[
\int_{\mathbb{T}} \big(|z|^2 + \alpha_0^2 |h|^2\big)^{1/2} dm \leq \int_{\mathbb{T}} |z + h| dm
\]
for all \(z\in\mathbb C\) and \(h\in H_0^1(\mathbb T)\). Combined with Brownian motion and a stopping time
\[
\rho = \inf\{t < \tau: |h(B_t)| > 2 \alpha_0^{-1}|z|\},
\]
this yields a “pointwise” analytic splitting \(h=g+b\) with \(g\in H_0^\infty\), \(b\in H_0^1\), and estimates that globalize to the martingale decomposition \(F=G+B\) [1009.3629].

More recently, the Davis-Garsia decomposition has also been identified as an intrinsic decomposition inside the Hardy martingale space itself:
\[
\mathbb{E} \left(\sum_{n=0}^\infty \left|\Delta_n f\right|^2\right)^\frac{1}{2}
\simeq
\inf_{f=g+h}
\mathbb{E}\sum_{n=0}^\infty \left|\Delta_n g\right|
+
\mathbb{E}\left(\sum_{n=0}^\infty \mathbb{E}\left(\left|\Delta_n h\right|^2\mid \mathcal{F}_{n-1}\right)\right)^\frac{1}{2},
\]
with \(g,h\) remaining in the Hardy martingale space. That analyticity-preserving form is the key input in the extension of Fefferman’s multiplier theorem from \(H^1(\mathbb T)\) to \(H^1(\mathbb T^{\mathbb N})\) [2509.07616].

## 3. Vector-valued Hardy martingales and Banach-space geometry

The vector-valued theory addresses a basic obstruction: the classical Davis decomposition does not respect the analytic structure when applied to Hardy martingales. For an \(X\)-valued Hardy martingale \(F=(F_k)_{k=1}^n\), where \(X\) is a complex Banach space, one has a Davis-type decomposition
\[
F=G+B
\]
into \(X\)-valued Hardy martingales such that
\[
\|\Delta G_k\|_X \le 2 \max_{m \le k-1} \|F_m\|_X
\]
and
\[
\sum_{k=1}^n \mathbb{E}\|\Delta B_k\|_X \le C\,\mathbb{E}\|F\|_X.
\]
A stronger version sharpens the predictable part to
\[
\|\Delta G_k\|_X \le C\|F_{k-1}\|_X,
\qquad
\sum_{k=1}^\infty \mathbb{E}\|\Delta B_k\|_X \le C\,\mathbb{E}\|F\|_X.
\]
These decompositions preserve the Hardy structure and mirror the scalar theory at the level of predictable control [1504.06513].

The proofs combine Brownian motion, stopping times, and Doob’s projection with nonlinear telescoping and maximal estimates; the refined decomposition also uses Havin’s Lemma, outer-function methods, and a vector-valued version of Bourgain’s functional inequalities for Hardy spaces. The analytic content is therefore not decorative: it is the mechanism that permits predictable control without leaving the class of Hardy martingales [1504.06513].

Under additional geometric assumptions on \(X\), the decomposition yields Davis-Garsia inequalities of \(L^q\)-type. If \(X\) satisfies hypothesis \((q)\), meaning that for any \(M\) there is \( \delta(M)>0 \) such that for \(\|x\|=1\) and \(g\in H_0^\infty(\mathbb T,X)\) with \(\|g\|_\infty\le M\),
\[
\int_\mathbb T \| z + g \|_X dm
\ge
\left( 1 + \delta(M) \int_\mathbb T \|g\|_X^q dm \right)^{1/q},
\]
then
\[
\left( \sum_{k=1}^\infty \mathbb{E}_{k-1} \|\Delta G_k\|_X^q \right)^{1/q}
+
\sum_{k=1}^\infty \mathbb{E}\|\Delta B_k\|_X
\le
A_q\,\mathbb{E}\|F\|_X.
\]
For \(q=2\) and \(X=\mathbb C\), this recovers the classical Garsia/Davis inequalities for scalar Hardy martingales [1504.06513].

Two applications are especially prominent. First, if a martingale multiplier \(T_\varepsilon\) has an \(L^2\) bound on square-integrable Hardy martingales, the decomposition extrapolates this to \(L^1\):
\[
\mathbb{E}\|T_\varepsilon F_k\|_X \le C \mathbb{E}\|F_k\|_X.
\]
Second, any \(X\)-valued Hardy martingale can be embedded as a subsequence into another Hardy martingale with small, predictable increments, providing an alternative route to extrapolation arguments in the style of Burkholder or Maurey [1504.06513].

## 4. Sine-cosine decomposition, dyadic perturbation, and quotient embeddings

For a martingale \(G=(G_k)\), the sine-cosine decomposition isolates the even and odd parts in the last variable. Writing \(x=(x_1,\ldots,x_{k-1})\) and \(y\in\mathbb T\),
\[
U_k(x,y)=\frac{1}{2}\bigl(G_k(x,y)+G_k(x,\overline y)\bigr)
\]
defines the cosine part \(U=(U_k)\), while \(V_k=G_k-U_k\) defines the sine part. Then \(G=U+V\), \(U_k(x,y)=U_k(x,\overline y)\), and \(\Delta V_k(x,y)=-\Delta V_k(x,\overline y)\). In previsible norm, one has the transform identity
\[
\|U\|_{\mathcal P}=\|T_W(G)\|_{\mathcal P},
\qquad
T_W(G)=\Im \sum w_{k-1}\Delta_k G,
\]
for adapted unimodular \(W=(w_k)\) [1611.02653].

This identity remains stable under dyadic perturbation. If \( \mathcal D \) is the dyadic sigma-algebra generated by \( \sigma_k(x)=\operatorname{sign}(\Re x_k) \), then
\[
\|U - \mathbb{E}(U|\mathcal D)\|_{\mathcal P}
\leq
C\, \|T_W(G - \mathbb{E}(G | \mathcal D))\|_{\mathcal P}^{1/2} \|G\|_{\mathcal P}^{1/2}.
\]
The estimate shows that the \( \mathcal P \)-norm relation between a Hardy martingale and its cosine part survives the subtraction of coarse dyadic information [1611.02653].

A related dyadic-perturbation theory for Davis-Garsia inequalities starts with a Hardy martingale \(F\) and a dyadic martingale \(D\). There exists a Hardy martingale \(G\) and a remainder \(B=F-G\) such that
\[
|\Delta G_k| \leq C_0 |F_{k-1} - D_{k-1}|,
\qquad
\mathbb{E} \sum |\Delta B_k| \leq C \mathbb{E}|F_n - D_n|,
\]
and
\[
\| T(G - D) \|_{\mathcal{P}}
\leq
C \|F - D\|_{L^1}^{1/2} \|F - D\|_{H^1}^{1/2},
\]
where
\[
T(H)=\Im \left[ \sum w_{k-1}\Delta H_k \right],
\qquad
w_{k-1}=\frac{\overline{F_{k-1}-D_{k-1}}}{|F_{k-1}-D_{k-1}|}.
\]
For \(D=\mathbb E_{\mathcal D}F\), this gives
\[
\|B\|_{\mathcal A}\leq C\|F-\mathbb E_{\mathcal D}F\|_{L^1}.
\]
These are analyticity-preserving perturbative versions of the Davis-Garsia mechanism [1209.3964].

The same circle of ideas leads to geometric consequences in quotient spaces. For dyadic martingales \(D\),
\[
\inf_{F\in \text{Hardy}} \|D-F\|_{L^1} \geq \delta \|D\|_{L^1},
\]
and Bourgain’s embedding \(L^1\hookrightarrow L^1/H^1\) is revisited through this lower bound and the dyadic stability of transform estimates. Both the dyadic perturbation theory and the sine-cosine stability estimate are used in the proof that \(L^1\) embeds into \(L^1/H^1\) [1209.3964][1611.02653].

## 5. Products with dual spaces, paraproducts, and endpoint commutators

The product of a Hardy-space martingale with an element of its dual does not generally remain in the same Hardy class, and its correct range is described by paraproduct decompositions. For \(f\in H^1(\Omega)\) and \(g\in \mathrm{BMO}(\Omega)\), there exist continuous bilinear operators
\[
\Pi_1: H^1(\Omega)\times \mathrm{BMO}(\Omega)\to L^1(\Omega),\qquad
\Pi_2: H^1(\Omega)\times \mathrm{BMO}(\Omega)\to H^1(\Omega),
\]
\[
\Pi_3: H^1(\Omega)\times \mathrm{BMO}(\Omega)\to H_\Phi^1(\Omega),
\qquad
\Phi(t)=\frac{t}{\log(e+t)},
\]
such that
\[
f\cdot g=\Pi_1(f,g)+\Pi_2(f,g)+\Pi_3(f,g).
\]
For finite martingale expansions,
\[
\Pi_1(f,g)=\sum_k d_kf\,d_kg,\qquad
\Pi_2(f,g)=\sum_k f_{k-1}d_kg,\qquad
\Pi_3(f,g)=\sum_k g_{k-1}d_kf,
\]
and analogous decompositions hold for \(0<p<1\) with the martingale Lipschitz/Campanato dual space \(\Lambda_{1,\alpha_p}\), where \(\alpha_p=1/p-1\) [2301.08723].

A sharper endpoint decomposition isolates bounded variation. For \(f\in H_1\) and \(g\in \mathrm{BMO}\),
\[
f\cdot g = G(f,g)+L(f,g),
\]
where
\[
L(f,g)_n:=\sum_{k=0}^n d_k(f)d_k(g)
\]
is a process with bounded variation and limit in \(L^1\), while \(G(f,g)\) belongs to the martingale Hardy-Orlicz space \(H_{\log}\) associated with
\[
\Phi(t):=\frac{t}{\log(e+t)}.
\]
The range \(L^1+H_{\log}\) is sharp: for particular martingales, it cannot be replaced by a smaller space having a larger dual [2301.09392].

This endpoint structure controls commutators. For \(b\in \mathrm{BMO}\), the largest subspace \(H_1^b\subset H_1\) for which the commutators
\[
[T,b](f):=T(bf)-bT(f)
\]
are bounded from \(H_1^b\) to \(L^1\) is characterized by
\[
\|f\|_{H_1^b}
=
\|f\|_{H_1}
+
\sup_{n\ge 0}\|[E_n,b](f)\|_{L^1},
\qquad
[E_n,b](f)=E_n(bf)-bE_n(f),
\]
and equivalently by
\[
f\in H_1^b
\iff
L(f,b)\in H_1.
\]
For operators in the class \( \mathcal K_1 \), one has
\[
\|f\|_{H_1^b}\sim \|f\|_{H_1}+\|L(f,b)\|_{H_1}\sim \|f\|_{H_1}+\|[T,b](f)\|_{L^1}.
\]
The relevant operator class includes martingale transforms, martingale fractional integrals, the Doob maximal operator, the martingale square function, the dyadic Hilbert transform, and the maximal operator of Césaro means for Walsh–Fourier series [2301.09392].

## 6. Extensions, related frameworks, and recent directions

Several adjacent theories extend the Hardy-martingale toolkit beyond the analytic torus model. Some of these concern martingale Hardy spaces rather than Hardy martingales in the strict sense, but they preserve the same structural themes: square functions, atomic decompositions, interpolation, perturbation stability, and endpoint duality.

| Framework | Representative result | Source |
|---|---|---|
| Non-homogeneous martingale Hardy spaces | \( \|Mf\|_1 \approx \|Sf\|_1 \), proved via the relation between Carleson sequences and balayage | [1411.5407] |
| Variable and mixed exponents | \(H^s_{p(\cdot)}\) admits a \((1,p(\cdot),\infty)\)-atomic decomposition; five mixed spaces \(H_{\mathbf p}^s,H_{\mathbf p}^S,H_{\mathbf p}^M,\mathcal P_{\mathbf p},\mathcal Q_{\mathbf p}\) have atomic decompositions and satisfy a mixed Burkholder-Davis-Gundy inequality | [1404.2395][1911.05356] |
| Noncommutative martingale Hardy spaces | Constructive atomic decomposition for \(\mathrm h_p\), \(0<p<2\); atomic block decomposition for arbitrary filtrations; K-closedness and real interpolation for \(H_p(\mathcal M)\); square-function inequalities for weakly differentially subordinate martingales | [2001.08775][1407.5451][2407.12335][1901.08752] |
| Fractional integration and irregular filtrations | For arbitrary filtrations, \( \mathfrak I_\alpha:L_p\to L_q \) is bounded when \( \frac1p-\frac1q=\alpha \), and \( \tilde{\mathcal I}_\alpha \) shares the same \(L_p\to L_q\) boundedness | [2009.05293] |
| Multipliers, interpolation, and partial sums | Fefferman’s multiplier theorem extends to \(H^1(\mathbb T^{\mathbb N})\); for holomorphic random variables, new truncation formulae yield Marcinkiewicz decompositions and \((H^1(\Omega),H^\infty(\Omega))_{1/2}=H^2(\Omega)\); on Vilenkin systems, subsequences of partial sums on \(H_p\), \(0<p<1\), satisfy sharp boundedness and convergence criteria | [2509.07616][1609.07364][1802.07669] |

Within noncommutative martingale Hardy spaces, several structural facts parallel the classical analytic theory. For \(0<p<2\), \(\mathrm h_p\) admits an explicit constructive decomposition into algebraic atoms, and for \(0<p\le 1\) the \((p,q)_c\)-atomic spaces coincide for all \(2\le q\le \infty\) [2001.08775]. For arbitrary filtrations, \(H_1(\Omega)\) and its noncommutative analogue admit atomic block decompositions,
\[
H_1(\Omega)\simeq H_{\mathrm{atb},p}(\Omega),
\qquad
H_1(\mathcal M)\simeq H_{\mathrm{atb},p}^{1,c}(\mathcal M),
\]
which replace classical atomic decompositions that fail outside the regular setting [1407.5451]. Real interpolation also persists: for \(1\le p,q\le\infty\), \((H_p(\mathcal M),H_q(\mathcal M))\) is K-complemented in \((L_p(\mathcal M,\ell_2),L_q(\mathcal M,\ell_2))\), and
\[
(H_{p_0}(\mathcal M),H_{p_1}(\mathcal M))_{\theta,p}=H_p(\mathcal M)
\]
with equivalent norms [2407.12335].

This suggests that the defining analytic increment condition of Hardy martingales is compatible with a broad transfer apparatus: square-function equivalences, atomic and block decompositions, interpolation, multiplier theory, and endpoint product estimates all survive, but typically only after the theory is reformulated so that analyticity, predictability, or noncommutative support conditions are preserved at every stage.

Source: https://www.emergentmind.com/topics/hardy-martingales