---
title: Sparse Domination Principles Overview
url: https://www.emergentmind.com/topics/sparse-domination-principles
type: topic
---

# Sparse Domination Principles Overview

Sparse domination principles are a family of estimates that majorize nonlocal operators by positive, local expressions indexed by sparse collections of cubes, balls, rectangles, or stopping regions. In their most standard form, they control either an operator pointwise by a sparse operator or its bilinear pairing by a sparse form; in more specialized settings they involve $\ell^r$-aggregation, square-function averages, cancellative percentiles, stopping times, or modified forms adapted to complexity and nondoubling geometry. Across one-parameter Calderón–Zygmund theory, metric spaces of homogeneous type, multiscale operators, vector-valued extensions, and martingale frameworks, sparse domination has become a unifying mechanism for deriving quantitative weighted inequalities, endpoint bounds, and structural necessary conditions. At the same time, the theory has sharp limitations: in product geometry, classical $(1,1)$-type sparse domination can fail outright for rectangles [1610.01958; 2009.00336; 1811.01243].

## 1. Basic objects and formulations

A sparse family of cubes $S$ in $\mathbb{R}^d$ is typically defined by the existence, for each $Q\in S$, of a measurable set $E_Q\subset Q$ with $|E_Q|\ge \eta |Q|$ such that the sets $\{E_Q\}_{Q\in S}$ are pairwise disjoint. The associated positive sparse bilinear form is
\[
\Lambda_S(f,g):=\sum_{Q\in S}|Q|\,\langle f\rangle_Q\,\langle g\rangle_Q,
\qquad 
\langle f\rangle_Q:=\frac{1}{|Q|}\int_Q |f(x)|\,dx,
\]
and the corresponding sparse operator is
\[
\mathcal A_S f(x):=\sum_{Q\in S}\langle |f|\rangle_Q\,\mathbf 1_Q(x).
\]
More generally, one uses $L^p$-averages,
\[
\langle f\rangle_{p,Q}:=\Big(\frac{1}{|Q|}\int_Q |f|^p\Big)^{1/p},
\]
and forms $\Lambda_S^{p,q}(f,g)=\sum_{Q\in S}|Q|\,\langle f\rangle_{p,Q}\langle g\rangle_{q,Q}$ [1610.01958].

In spaces of homogeneous type $(X,d,\mu)$, the same paradigm is stated for balls:
\[
\Lambda_{\mathcal S}^{p,q}(f,g)=\sum_{B\in\mathcal S}\mu(B)\,\langle |f|\rangle_{p,B}\,\langle |g|\rangle_{q,B},
\qquad
\langle |f|\rangle_{p,B}:=\Big(\frac{1}{\mu(B)}\int_B |f|^p\,d\mu\Big)^{1/p}.
\]
The sparse family is now a collection of measurable sets with pairwise disjoint major subsets of comparable measure [2009.00336].

Several variants arise when positivity alone is insufficient. In the $\ell^r$ framework on spaces of homogeneous type, the sparse operator becomes
\[
A_{S,a,p_0}^{(r)}f(x):=\Big(\sum_{Q\in S}\langle |f|\rangle_{p_0,aQ}^r\,\mathbf 1_Q(x)\Big)^{1/r},
\]
with a dilation parameter $a\ge 1$ [1907.00690]. In cancellative martingale settings one replaces local averages by conditional percentiles and uses
\[
C_S f(x):=\sum_{k\ge 0}\mathbf 1_{S_k}P_k^r(f)
\]
for an $\eta$-sparse adapted sequence $S=\{S_k\}_{k\ge 0}$ [2603.08446]. In balanced nondoubling filtrations, classical sparse forms are supplemented by a modified term
\[
C^{(N)}_{\mathcal S}(f,g)
\]
that records interactions between distinct sparse atoms at dyadic distance $\le N+2$ [2309.13943]. In bi-parameter dyadic analysis, a natural form is
\[
\Lambda_{\mathcal S}(f,g):=\sum_{R\in\mathcal S}|R|\,(Sf)_R\,(Sg)_R,
\]
with $R$ dyadic rectangles and $S$ the bi-parameter square function [1709.05009].

These formulations separate two basic notions. Pointwise sparse domination asserts an inequality such as
\[
|Tf(x)|\le C\,\mathcal A_S f(x),
\]
while bilinear sparse domination asserts
\[
|\langle Tf,g\rangle|\le C\,\Lambda_S(f,g).
\]
Much of the modern theory moves between these two forms, but the distinction becomes decisive in product and nondoubling settings.

## 2. Principal proof architectures

A central proof scheme is the Calderón–Zygmund decomposition combined with a recursive stopping-time argument. For dyadic shifts, one selects maximal subcubes $I\subset Q$ satisfying
\[
\langle f_j\rangle_I>2^8\langle f_j\rangle_Q
\]
for at least one input, obtaining the packing estimate
\[
\sum_{I\in\mathbb I_Q}|I|\le 2^{-7}|Q|.
\]
The key iterative lemma then bounds the localized form on $Q$ by a main sparse term plus recursively localized descendants; the dependence on the shift complexity $\rho$ is linear [1610.01958].

A closely related but more explicit kernel decomposition appears for Dini-smooth Calderón–Zygmund operators. One writes
\[
K_s(x,y)=\psi\!\left(\frac{x-y}{2^s}\right)K(x,y),
\]
decomposes $T$ into single-scale pieces, constructs an exceptional set by maximal functions, performs a Whitney decomposition, and estimates good–good, good–bad, and bad–bad interactions. The Dini integral
\[
\int_0^1\frac{\omega(t)}{t}\,dt
\]
enters precisely in summing oscillations across scales [2509.07659].

An alternative architecture dispenses with dyadic cubes altogether. In spaces of homogeneous type, sparse domination can be derived from a single-scale $L^p$-improving hypothesis. The assumptions are localization
\[
\mathrm{supp}\big(T_\sigma^s[f\,\mathbf 1_L]\big)\subset c_0L,
\]
uniform truncation bounds
\[
\sup_{\sigma<\tau}\|T_\sigma^\tau\|_{L^p\to L^p}<\infty,
\]
and the two single-scale estimates
\[
\left\langle T(s)(f\,\mathbf 1_L)\right\rangle_{p_2',\gamma_2L}
\le I_{p_1,p_2}\,\langle f\rangle_{p_1,L},
\]
\[
|\langle T(s)(f\,\mathbf 1_L),b\rangle|
\le \omega\!\left(\frac r{2^s}\right)\mu(L)\,\langle f\rangle_{p_1,L}\,\langle b\rangle_{p_2,\gamma_2L},
\]
for $(p_2,r)$-atoms $b$. Whitney coverings of balls replace dyadic systems, and stopping forms are summed by Carleson packing [2009.00336].

Other proof mechanisms localize oscillation rather than kernels. An improved pointwise principle assumes a localized $W_q$ property for $T$ and weak type for the oscillation operator
\[
M_a^\# f(x):=\sup_{Q\ni x}\operatorname*{ess\,sup}_{x',x''\in Q}
\big|T(f\mathbf 1_{\mathbb R^n\setminus aQ})(x')-T(f\mathbf 1_{\mathbb R^n\setminus aQ})(x'')\big|,
\]
which is weaker than assuming weak type for Lerner’s earlier grand maximal truncated operator $\mathcal M_T$ [1901.00195]. An even more abstract version is operator-free: one works with arbitrary families $f_Q,f_{P,Q}$ satisfying an $l^r$-condition and a quasi-locality estimate, and produces sparse domination directly for the family rather than for a specific operator [2106.16202].

## 3. Canonical positive theorems

Uniform sparse domination for dyadic shifts provides one of the basic one-parameter theorems. For every dyadic lattice $\mathcal D$ and every compactly supported $f_1,f_2\in L^1(\mathbb R^d)$, there exists a sparse collection $S_{\mathcal D}$ such that for all complexities $\rho\ge 1$,
\[
\sup_{\mathbb S^\rho}|\mathbb S^\rho(f_1,f_2)|
\le C\rho \sum_{Q\in S_{\mathcal D}}|Q|\langle f_1\rangle_Q\langle f_2\rangle_Q.
\]
Combined with Hytönen’s dyadic representation, this yields sparse domination for the class $\mathcal U$ of singular integrals satisfying the assumptions of the classical $T(1)$-theorem [1610.01958].

Rough singular integrals and Bochner–Riesz means admit bilinear sparse domination of $(1,p)$ type. For the rough homogeneous singular integrals $T_\Omega$ with $\Omega\in L^q(S^{d-1})$, $1<q\le \infty$, and mean zero, one has
\[
|\langle T_\Omega f_1,f_2\rangle|
\le \frac{Cp}{p-1}\sup_S \Lambda_S^{(1,p)}(f_1,f_2),
\]
with dependence on $\|\Omega\|_{L^{q,1}\log L(S^{d-1})}$ when $1<q<\infty$ and on $\|\Omega\|_{L^\infty(S^{d-1})}$ when $q=\infty$. The same scheme covers Bochner–Riesz means at the critical index [1612.09201].

Square-function operators also fit the paradigm. For the Marcinkiewicz integral $\mu_\Omega$ with $\Omega\in L^\infty(S^{n-1})$, homogeneous of degree zero and mean value zero, there exists a $1/3$-sparse family $\mathcal S$ such that for $r\in(1,2]$, $s\in[1,2]$,
\[
|\langle \mu_\Omega f,g\rangle|
\le C_n\|\Omega\|_{L^\infty(S^{n-1})}\,r'\,
\Lambda_{\mathcal S}^{r,s}(f,g).
\]
Here the square-function nature of $\mu_\Omega$ dictates the restriction $s\le 2$ [1905.10755].

The metric-space theorem based on single-scale improving recovers sparse domination for Dini-continuous Calderón–Zygmund kernels on spaces of homogeneous type and extends further to maximal functions associated to convolutions with measures exhibiting Fourier decay, as well as to Radon transforms along polynomial submanifolds of $\mathbb R^n$ [2009.00336]. A complementary multiscale theorem covers sums $\sum_j T_j$, maximal operators, square functions, and variation-norm operators under support near the diagonal, uniform weak-type $(p,p)$ and restricted strong-type $(q,q)$ bounds for partial sums, a single-scale $L^p\to L^q$ improving estimate, and an $\epsilon$-regularity condition; these hypotheses are also shown to be necessary in the stated sense [2009.00227].

A common misconception is that sparse domination is tied to smooth kernels. The rough-kernel results above, together with sparse bounds under $L^r$-Hörmander regularity and weak endpoint assumptions, show that the method survives substantial loss of smoothness, provided one has the correct local oscillation or improving estimates.

## 4. Vector-valued, multilinear, and operator-free extensions

For multilinear singular integrals, sparse domination can be proved under the multilinear $L^r$-Hörmander condition. If $T$ is bounded from $L^r\times\cdots\times L^r$ to weak $L^{r/m}$ and its kernel satisfies the $m$-linear $L^r$-Hörmander condition, then for compactly supported $f_j\in L^r(\mathbb R^n)$ there exists a sparse family $\mathcal S$ such that
\[
|T(f_1,\dots,f_m)(x)|
\le C\sum_{Q\in\mathcal S}\Big(\prod_{j=1}^m \langle |f_j|\rangle_{r,Q}\Big)\mathbf 1_Q(x)
\]
almost everywhere [1606.03925].

Scalar sparse domination can also be lifted systematically to vector-valued sparse domination. For tuples of quasi-Banach function spaces one introduces the multilinear analogue of the UMD condition, denoted $\mathrm{UMD}_{\vec r,s}$, characterized by the boundedness of the multisublinear Hardy–Littlewood maximal operator. Under the assumption $X^q\in \mathrm{UMD}_{\vec r,s}$, a scalar sparse bound of the form
\[
\|T(f)\cdot g\|_{L^q(\mathbb R^d)}
\le C_T\|\mathcal M_{(\vec r,q)}(f;g)\|_{L^q(\mathbb R^d)}
\]
implies
\[
\|\|T(f)(\cdot,\cdot)\|_X\,|g(\cdot)|\|_{L^q(\mathbb R^d)}
\le C_XC_T\|\mathcal M_{(\vec r,q)}(\|f\|_X,|g|)\|_{L^q(\mathbb R^d)}.
\]
The framework covers iterated Lebesgue, Lorentz, and Orlicz spaces and yields sharp vector-valued weighted bounds directly from scalar sparse domination, without Rubio de Francia extrapolation [2003.02233].

A related Banach-space-valued development replaces the usual $\ell^1$ sparse operator by an $\ell^r$ sparse operator dictated by the operator’s structure and the geometry of the underlying Banach spaces. On spaces of homogeneous type, this produces pointwise $\ell^r$-sparse domination for vector-valued operators, an $A_2$-theorem for vector-valued Calderón–Zygmund operators, an anisotropic mixed norm Mihlin multiplier theorem, and quantitative weighted inequalities for the Rademacher maximal operator [1907.00690].

Sparse domination via the helicoidal method extends the theory to multilinear Fourier multipliers whose symbols are singular along a $k$-dimensional subspace of
\[
\Gamma=\{\xi_1+\cdots+\xi_{n+1}=0\},
\]
with $k<\frac{n+1}{2}$, and to the variational Carleson operator. The method yields sparse estimates for multiple vector-valued extensions and Fefferman–Stein type inequalities in $L^q(v^q)$ for all $0<q<\infty$ [1707.05484].

The operator-free framework goes further still. For arbitrary families $f_Q,f_{P,Q}$ satisfying an $l^r$-condition and quasi-locality, one obtains sparse domination results that apply to generalized Poincaré–Sobolev inequalities, tent spaces, general dyadic sums, and vector-valued square functions that are not localizable in the sense of earlier operator-based schemes [2106.16202]. *This suggests that sparse domination is not merely an operator theorem but a structural principle for local oscillation and aggregation.*

## 5. Cancellative, martingale, and product-geometry variants

In continuous time, sparse domination can be formulated through stopping times rather than cubes. If $X$ and $Y$ are uniformly integrable càdlàg Hilbert-space-valued martingales and $Y$ is differentially subordinate to $X$, then there exists a sparse increasing sequence of stopping times $\{T_j\}_{j\ge 0}$ such that
\[
Y^*\le 8\sum_{j=0}^\infty |X|_{T_j}\,\mathbf 1_{E_j}
\quad\text{a.s.}
\]
This pathwise domination leads to weighted martingale bounds in terms of the martingale $A_p$ characteristic $Q_p(w)$, and in the special case $Y=X$ one gets the optimal exponent $Q_p(w)^{1/(p-1)}$ [1607.06319].

A distinct development is cancellative sparse domination. Instead of dominating by noncancellative averages, one precomposes sparse operators with percentile maximal functions. In the martingale setting,
\[
M_S f(x):=\sup_{k\ge 0}\mathbf 1_{S_k}(x)|E_kf(x)|,
\qquad
A_S f(x):=\sum_{k\ge 0}\mathbf 1_{S_k}(x)|E_kf(x)|,
\]
and the paper proves the sparse characterization
\[
\|f\|_{H_D^1}\eqsim \sup_{S\text{ sparse}}\|M_Sf\|_{L^1}
\eqsim \sup_{S\text{ sparse}}\|A_Sf\|_{L^1}.
\]
It also yields cancellative sparse domination for martingale transforms, square functions, dyadic Haar shifts, and smooth Calderón–Zygmund operators, and gives the bi-parameter bound
\[
M_{D\times D}f(x)\lesssim P_{D\times D}^{1/2}(M_Sf)(x)
\]
for the dyadic strong maximal function [2603.08446].

Product geometry requires further modifications. For bi-parameter martingale transforms, cancellative dyadic shifts, and paraproduct-free Journé operators, one has sparse domination not by averages of $|f|$ but by averages of square functions:
\[
|\langle Tf,g\rangle|
\lesssim \sum_{R\in\mathcal S}|R|\,(Sf)_R\,(Sg)_R.
\]
This reflects the Chang–Fefferman product structure and the role of Córdoba–Fefferman selection in place of one-parameter stopping-time trees [1709.05009].

In general filtered measure spaces with balanced measures, classical sparse domination may fail and must be replaced by a sparse domination-like inequality
\[
|\langle T_{s,t}f,g\rangle|
\le C\big(\Lambda_{\mathcal S}(f,g)+C_{\mathcal S}^{(N)}(f,g)\big),
\qquad N=s+t,
\]
where $C_{\mathcal S}^{(N)}$ records interactions up to dyadic distance $N+2$. The associated admissible weight classes depend on the complexity $N$ [2309.13943]. This is a sharp indication that cancellation, filtration regularity, and product geometry can force sparse principles away from the classical one-parameter model.

## 6. Weighted inequalities and quantitative sharpness

Sparse domination is especially effective because positive sparse forms admit precise weighted estimates. In the dyadic-shift/T(1) framework, sparse domination recovers Hytönen’s sharp $A_p$ bounds in the scalar setting and yields, for matrix weights,
\[
\|T\otimes \mathrm{Id}_{\mathbb R^n}\|_{L^2(W;\mathbb R^n)\to L^2(W;\mathbb R^n)}
\le C [W]_{A_2}^{3/2},
\]
uniformly over $T\in\mathcal U$ [1610.01958].

In spaces of homogeneous type, if $T$ satisfies the single-scale improving hypotheses with exponents $p_1<p<p_2'$, then for weights
\[
w\in A_{p/p_1}\cap RH_{p_2'/p}
\]
one has
\[
\|T:L^p(w)\to L^p(w)\|
\lesssim
\Big([w]_{A_{p/p_1}}\,[w]_{RH_{(\frac{p_2'}{p})'}}\Big)^{\max\!\left(\frac{1}{p-p_1},\,\frac{p_2'-1}{p_2'-p}\right)}.
\]
The same machinery yields vector-valued inequalities and endpoint consequences where appropriate [2009.00336].

For rough homogeneous singular integrals $T_\Omega$, sparse domination gives quantitative weighted estimates. If $1<q<\infty$ and $\Omega\in L^{q,1}\log L(S^{d-1})$, then for $q'<t<\infty$,
\[
\|T_\Omega f\|_{L^t(w)}
\le C(t,q)\|\Omega\|_{L^{q,1}\log L(S^{d-1})}\,[w]_{A_{t/q'}}^{\alpha(t,q)}\|f\|_{L^t(w)},
\]
with
\[
\alpha(t,q)=\max\{1,1/(t-q')\}.
\]
When $\Omega\in L^\infty(S^{d-1})$ one recovers the exponent
\[
\beta(t)=\max\{t,2\}/(t-1)
\]
for both $T_\Omega$ and critical Bochner–Riesz means [1612.09201].

For the Marcinkiewicz integral $\mu_\Omega$ with rough $L^\infty$ kernel, the bilinear sparse domination implies
\[
\|\mu_\Omega f\|_{L^p(w)}
\le C_{n,p}\|\Omega\|_{L^\infty(S^{n-1})}
[w]_{A_p}^{\max\{1/2,\,1/(p-1)\}}
\big([w]_{A_\infty}+[w^{1-p'}]_{A_\infty}\big)\|f\|_{L^p(w)},
\]
as well as the endpoint estimate
\[
\|\mu_\Omega f\|_{L^{1,\infty}(w)}
\le C_n\|\Omega\|_{L^\infty(S^{n-1})}
[w]_{A_1}[w]_{A_\infty}\log(e+[w]_{A_\infty})\,\|f\|_{L^1(w)}.
\]
A Coifman–Fefferman type inequality
\[
\|\mu_\Omega f\|_{L^p(w)}
\le C_{n,p}\|\Omega\|_{L^\infty(S^{n-1})}[w]_{A_\infty}\|Mf\|_{L^p(w)}
\]
also follows [1905.10755].

For variational Carleson, the sparse bound yields
\[
\|C^{var,r}\|_{L^p(w)\to L^p(w)}
\lesssim
[w]_{A_{p/(r'+\delta)}}^{\max\left(\frac{1}{p-r'-\delta},1\right)},
\]
for $\delta>0$ and $r'+\delta<p<\infty$ [1707.05484]. In the cancellative Hardy-space setting, smooth Calderón–Zygmund operators satisfy
\[
\|Tf\|_{L^p(w)}
\lesssim_{T,p,q,s}
[w]_{A_q}^{1/p}[w]_{A_\infty}^{(1-\frac1p)_+}\,\|M^sf\|_{L^p(w)},
\]
and if $s\ge d(\frac qp-1)$ then $T$ is bounded $H^p(w)\to L^p(w)$ [2603.08446].

## 7. Obstructions, modified forms, and open directions

Sparse domination is not universal in its classical $(1,1)$ form. For the dyadic strong maximal function $\mathcal M_S$ on axis-parallel rectangles in $\mathbb R^2$, the main theorem states that for every $0<\eta<1$ and every $C>0$ there exist compactly supported integrable $f,g$ such that for every $\eta$-sparse collection $\mathcal S$ of dyadic rectangles,
\[
|\langle \mathcal M_S f,g\rangle|
\ge C\sum_{R\in\mathcal S}|R|\,\langle |f|\rangle_R\,\langle |g|\rangle_R.
\]
Quantitatively, the construction produces discrete measures $\mu,\nu$ with
\[
\langle \mathcal M_S(\mu),\nu\rangle \gtrsim 2^k,
\qquad
\sum_{R\in\mathcal S}|R|\,\langle \mu\rangle_R\,\langle \nu\rangle_R
\lesssim \Lambda\,k\Bigl(1+\frac{2^{2k}}{m}\Bigr),
\]
and the choice $m=k2^{2k}$ forces any Carleson constant $\Lambda$ to blow up [1811.01243].

The obstruction is geometric. The proof uses an anisotropic “distance”
\[
\mathrm{dist}(p,q)
=
\inf\Bigl\{|R|^{1/2}: R\in\mathcal D\times\mathcal D,\ p,q\in R\Bigr\},
\]
a maximally separated set $\mathcal P$, and a near-far cloud $\mathcal Z$ with unique associations and controlled clustering. The same mechanism rules out $(1,1)$-type sparse domination for bi-parameter martingale transforms, excludes $(\phi,1)$ sparse bounds for
\[
\phi(x)=x[\log(e+x)]^\alpha,\qquad 0<\alpha<1/2,
\]
and propagates the failure to all dimensions $n\ge 2$ by tensor reduction [1811.01243].

A second obstruction occurs in balanced nondoubling filtrations. There, classical sparse domination fails even for positive complexity Haar shifts, and the correct weighted classes depend on the complexity $N=s+t$ through characteristics such as
\[
[w]_{A_p^{(N)}(\mu)}.
\]
The modified sparse form $C_{\mathcal S}^{(N)}$ and the maximal operator $M_N$ are not artefacts of proof but necessary features of the theory [2309.13943].

In bi-parameter harmonic analysis this does not mean that all sparse control disappears. Instead, one often has to replace averages of $|f|$ by averages of square functions, or replace classical sparse forms by cancellative or complexity-sensitive versions [1709.05009; 2603.08446]. *A plausible implication is that “sparse domination” in product and nondoubling settings is better understood as a family of domination mechanisms rather than a single canonical formula.*

Several open directions are explicit in the literature. For the strong maximal function and related bi-parameter operators, it remains open whether $(\phi,1)$-type sparse bounds hold for $\phi(x)=x[\log(e+x)]^\alpha$ with $\alpha\ge 1/2$, whether mixed sparse forms or restricted classes of rectangles can recover useful domination, and which other multiparameter operators admit salvageable sparse frameworks [1811.01243]. In the metric Whitney-covering approach, extending beyond Dini moduli and broadening maximal Radon applications remain open [2009.00336]. In the Calderón–Zygmund decomposition approach to Dini kernels, transferring the scheme to nondoubling measures is obstructed by the lack of classical Whitney coverings and good dyadic geometry for arbitrary measures [2509.07659].

Source: https://www.emergentmind.com/topics/sparse-domination-principles