Carleson operators are maximally modulated singular integrals that characterize convergence phenomena in Fourier analysis and underpin modern time-frequency methods.
Generalizations include polynomial, oscillatory, and weighted variants that use specialized modulation symmetries and curvature to achieve L^p boundedness.
Recent developments extend Carleson theory to abstract metric spaces, bilinear and anisotropic settings, revealing new insights into oscillatory behavior and endpoint estimates.
Carleson operators are maximally modulated singular integral operators whose prototype is the classical operator
Cf(x)=N∈Rsup∫Rf(x−t)eiNttdt,
a central object in harmonic analysis because its boundedness is equivalent to almost-everywhere convergence phenomena for Fourier expansions and partial Fourier integrals (Ramos, 2020). In contemporary usage, the term also encompasses polynomial, oscillatory, variational, anisotropic, multilinear, and abstract metric-space analogues, all of which preserve the characteristic feature of taking a supremum over a modulation family while retaining a singular kernel. The subject now spans classical time-frequency analysis, oscillatory integral theory, weighted inequalities, vector-valued extensions, and non-Euclidean generalizations (Lie, 2011).
1. Classical operator and the convergence problem
The classical Carleson operator can be written either as the maximally modulated Hilbert transform
Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),
or in the oscillatory singular-integral form displayed above (Plinio et al., 2013). In one dimension, this operator governs maximal partial Fourier integrals and is the analytic core of Carleson’s theorem on almost-everywhere convergence of Fourier series.
A detailed modern metric-space presentation states the classical theorem in quantitative Egorov form: for a 2π-periodic uniformly continuous f with ∣f(x)∣≤1, and every 0<ϵ<1, there is a Borel set E⊂[0,2π] with ∣E∣≤ϵ and an integer N0 such that
∣f(x)−SNf(x)∣≤ϵ
for all Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),0 and all Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),1 (Becker et al., 2024). That work also makes explicit the reduction from Fourier partial sums to a real-line Carleson maximal operator built from the localized kernel
Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),2
thereby exhibiting the classical theorem as a special case of a broader Carleson framework on doubling metric measure spaces (Becker et al., 2024).
Classical boundedness is now understood in several complementary forms. Carleson’s original theorem gives the Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),3 result, while later developments yield Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),4-boundedness for every Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),5 (Ramos, 2020). In the weighted setting, Di Plinio and Lerner proved that the classical operator satisfies the same linear Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),6-dependence as the Hilbert transform for Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),7: Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),8
and also obtained the Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),9-bound
2π0
showing that weighted Carleson theory reflects the near-2π1 endpoint behavior of Fourier convergence (Plinio et al., 2013).
2. Polynomial and oscillatory generalizations
A fundamental nonlinear extension is the polynomial Carleson operator
2π2
whose one-dimensional 2π3-boundedness for 2π4 was proved by Victor Lie, resolving the one-dimensional case of Stein’s conjecture (Ramos, 2020). Lie’s formulation emphasizes the generalized modulation symmetries
2π5
and recasts the proof in a higher-order wave-packet framework adapted to polynomial phases rather than linear frequencies (Lie, 2011).
That higher-order framework replaces ordinary time-frequency rectangles by curved polynomial tiles. In Lie’s one-dimensional construction, a tile is a 2π6-tuple
2π7
encoding a localized neighborhood of a polynomial graph over 2π8. The operator is decomposed as
2π9
and the proof proceeds through trees, forests, density and geometric factors, and a mass decomposition that yields BMO control of counting functions and eliminates exceptional sets from the analysis (Lie, 2011). One consequence highlighted in that paper is a direct strong f0 argument for generalized Carleson operators, without passing through interpolation from weak-f1 (Lie, 2011).
A more restricted polynomial family fixes the nonlinear phase and maximizes only over a linear modulation: f2
Ramos proved that f3 is bounded on f4, f5, with bounds depending only on f6 and f7, and did so by reducing the problem to the classical linear Carleson operator together with a quadratic analysis, rather than using the full polynomial Carleson machinery (Ramos, 2020). The result is conceptually narrower than the full theorem of Lie and Zorin-Kranich, but it isolates a mechanism by which linear Carleson control forces boundedness for a nonlinear modulation family (Ramos, 2020).
Oscillatory Carleson operators with nonlinear monomial phase form another major branch: f8
For f9, ∣f(x)∣≤10, Ramos proved ∣f(x)∣≤11-boundedness for ∣f(x)∣≤12 together with the uniform estimate
∣f(x)∣≤13
thereby controlling the constants in a neighborhood of the delicate quadratic regime ∣f(x)∣≤14 (Ramos, 2020). Earlier proofs treated ∣f(x)∣≤15 and ∣f(x)∣≤16 by sharply different methods and lost control as ∣f(x)∣≤17; the novelty is the uniformity in ∣f(x)∣≤18 (Ramos, 2020). The same paper notes that the analogous odd-phase family with ∣f(x)∣≤19 admits the same local uniformity near 0<ϵ<10, while the regimes0<ϵ<11 and especially 0<ϵ<12 remain more delicate (Ramos, 2020).
3. Weighted, variational, and vector-valued forms
Weighted Carleson theory has developed into a precise extension of sharp Calderón–Zygmund inequalities. Di Plinio and Lerner treated the classical Carleson operator 0<ϵ<13, its lacunary version 0<ϵ<14, and the Walsh-Carleson operator 0<ϵ<15 using local mean oscillation and sparse domination rather than the original time-frequency decomposition (Plinio et al., 2013). Their abstract framework applies to maximally modulated Calderón–Zygmund operators
For the classical operator they proved the weak estimate
0<ϵ<19
while for the lacunary operator
E⊂[0,2π]0
and for the Walsh-Carleson operator
E⊂[0,2π]1
(Plinio et al., 2013). The resulting E⊂[0,2π]2-bounds exhibit a hierarchy: a doubly logarithmic loss for E⊂[0,2π]3, a single logarithm for E⊂[0,2π]4, and no logarithmic loss for E⊂[0,2π]5, mirroring the corresponding endpoint theories (Plinio et al., 2013).
A different refinement replaces maximal control by variation control. For Banach-space-valued functions E⊂[0,2π]6, partial Fourier integrals are
E⊂[0,2π]7
and the E⊂[0,2π]8-variational Carleson operator is
E⊂[0,2π]9
Hytönen, Lorist, and Veraar proved that if ∣E∣≤ϵ0 is an ∣E∣≤ϵ1-intermediate UMD space, then
∣E∣≤ϵ2
for all ∣E∣≤ϵ3 and ∣E∣≤ϵ4 (Amenta et al., 2020). This strengthens vector-valued Carleson–Hunt theory by controlling the oscillation of the truncation path ∣E∣≤ϵ5, not merely its supremum (Amenta et al., 2020).
Their method combines time-frequency analysis with outer Lebesgue spaces on the time-frequency-scale domain ∣E∣≤ϵ6, using a wave-packet embedding
∣E∣≤ϵ7
and a truncated embedding ∣E∣≤ϵ8 adapted to the linearized variation increments (Amenta et al., 2020). The result recovers the maximal vector-valued Carleson theorem for intermediate UMD spaces as the limiting case ∣E∣≤ϵ9, but its actual content is stricter: finite N00-variation yields almost-everywhere convergence with quantitative oscillation control (Amenta et al., 2020).
4. Bilinear, anisotropic, and geometric variants
Carleson-type maximal modulation also appears in bilinear settings. The non-resonant bilinear Hilbert–Carleson operators are
N01
with N02 interpreted as N03 or N04 for noninteger N05 (Benea et al., 2021). These operators interpolate between the bilinear Hilbert transform, recovered at N06, and nonlinear-phase Carleson operators, recovered when one input is frozen to N07 (Benea et al., 2021).
For every non-resonant exponent
N08
Guo, Roos, and Yung proved that N09 maps ∣f(x)−SNf(x)∣≤ϵ0 to ∣f(x)−SNf(x)∣≤ϵ1 whenever
∣f(x)−SNf(x)∣≤ϵ2
and they established decay for high-oscillation pieces: ∣f(x)−SNf(x)∣≤ϵ3
for some ∣f(x)−SNf(x)∣≤ϵ4 (Benea et al., 2021). The excluded exponents ∣f(x)−SNf(x)∣≤ϵ5 are resonant because additional modulation symmetries destroy the scale-decay mechanism (Benea et al., 2021). The paper describes the resulting operators as genuinely hybrid, combining zero-curvature features inherited from bilinear Hilbert transform modulation symmetry with nonzero-curvature oscillatory cancellation (Benea et al., 2021).
In higher-dimensional anisotropic Euclidean settings, Roos studied
∣f(x)−SNf(x)∣≤ϵ6
where ∣f(x)−SNf(x)∣≤ϵ7 is a bounded ∣f(x)−SNf(x)∣≤ϵ8 multiplier homogeneous of degree ∣f(x)−SNf(x)∣≤ϵ9 under anisotropic dilations
Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),00
(Roos, 2017). If Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),01, then
Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),02
recovering the classical weak Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),03 Carleson bound when Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),04, Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),05 (Roos, 2017). The proof extends the Lacey–Thiele tree method to anisotropic tiles, anisotropic wave packets, and an anisotropic cone decomposition (Roos, 2017).
A geometric Radon-transform analogue appears in polynomial Carleson operators along the monomial curve Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),06: Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),07
Guo and Pierce studied partial linearizations in which the stopping-time function depends on only one ambient variable: Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),08
and proved Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),09-boundedness for all Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),10 when Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),11 and Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),12, as well as for the special modulation-invariant cases Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),13 with Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),14 and Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),15 with Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),16 (Guo et al., 2016). Their analysis combines Stein–Wainger-style Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),17 estimates, including phases with fractional monomials, with vector-valued Carleson–Hunt bounds in the symmetric cases (Guo et al., 2016). The same paper observes that Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),18 bounds for certain partial curved operators imply the usual Carleson theorem or even the quadratic Carleson theorem, so these are not merely auxiliary models (Guo et al., 2016).
A related higher-dimensional geometric variant acts along the paraboloid in Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),19: Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),20
where the phase has no linear term and its quadratic homogeneous part is not a nonzero multiple of Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),21 (Pierce et al., 2015). For Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),22, Pierce proved
Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),23
using a Littlewood–Paley decomposition and square-function comparison with a smoother auxiliary operator, precisely because the usual stopping-time Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),24 method fails in the Radon setting (Pierce et al., 2015).
5. Abstract metric-space formulations
Recent work has detached Carleson theory from Euclidean algebra and rebuilt it on doubling metric measure spaces. In that formulation, the ambient singular kernel is a one-sided Calderón–Zygmund kernel Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),25 satisfying
Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),26
together with Hölder regularity in the second variable (Becker et al., 2024). The generalized Carleson operator is then
Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),27
where Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),28 is a compatible cancellative collection of phase functions equipped, for each ball Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),29, with a metric Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),30 measuring relative oscillation (Becker et al., 2024).
The metric-space theorem of Buss, Cladek, and others proves a restricted weak-type bound: if the associated non-tangential maximal truncation Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),31 satisfies
Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),32
then for Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),33, Borel sets Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),34, and Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),35,
Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),36
(Becker et al., 2024). The 2025 sequel sharpens and streamlines the abstract theory, formulating main results as Theorem 1.1 and Theorem 1.2, and states that these results have been computer verified in Lean/mathlib, with the 2024 paper serving as the detailed formalization-oriented blueprint (Becker et al., 7 Aug 2025).
The central structural innovation is the axiomatization of the phase family. A compatible collection Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),37 is required to satisfy oscillation control
Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),38
monotonicity and doubling-type properties for the local metrics Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),39, and a covering condition in phase space (Becker et al., 2024). A cancellative collection also satisfies an oscillatory integral estimate against Lipschitz test functions on each ball, abstracting Euclidean van der Corput or integration-by-parts arguments (Becker et al., 2024). This framework is explicitly designed to include the classical linear phases, polynomial phases, and related non-Euclidean modulation classes (Becker et al., 7 Aug 2025).
The proof architecture remains recognizably time-frequency-theoretic: Christ-type dyadic cubes, tiles Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),40, antichains, trees, forests, density functionals, and Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),41-style tile correlation estimates all reappear in metric form (Becker et al., 2024). What changes is that Euclidean frequency intervals are replaced by modulation cells inside Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),42, and derivative-based phase separation is replaced by the local metrics Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),43 (Becker et al., 7 Aug 2025).
6. Methods, structural themes, and terminological boundaries
Across the subject, three proof paradigms recur. The first is classical time-frequency analysis: tiles, trees, forests, and almost orthogonality in the style of Fefferman and Lacey–Thiele. Lie’s polynomial theorem extends this to higher-order wave packets and curved tiles adapted to polynomial phase graphs (Lie, 2011). The second is oscillatory integral analysis via Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),44, stationary phase, and van der Corput estimates, especially in Stein–Wainger-type settings, oscillatory Carleson families, and Radon-transform analogues (Ramos, 2020). The third is domination technology—sparse domination, local mean oscillation, and outer Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),45 embeddings—which has been particularly effective in weighted and vector-valued questions (Plinio et al., 2013).
Several papers emphasize that the main technical difficulty is not merely singularity, but the interaction between maximal modulation and underlying symmetry. In quadratic or near-quadratic regimes, modulation symmetry can survive strongly enough to force a time-frequency treatment, even when oscillatory curvature is present (Ramos, 2020). In bilinear Hilbert–Carleson theory, this interaction is explicit: non-resonant exponents permit curvature-driven decay, whereas the resonant cases Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),46 and Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),47 produce additional modulation invariances that obstruct the available mechanism (Benea et al., 2021). In the monomial-curve and paraboloid settings, curvature from the underlying geometry and oscillation from the phase are inseparable, and partial or smoothed models become necessary (Guo et al., 2016).
The available results also leave a clear map of unresolved regimes. Uniform bounds for oscillatory Carleson operators are established only locally near Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),48, not globally in Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),49 (Ramos, 2020). The resonant bilinear Hilbert–Carleson cases Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),50 remain open (Benea et al., 2021). Several partial operators along monomial curves, including Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),51, Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),52, and Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),53, are not covered by full Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),54 theory (Guo et al., 2016). In the variational vector-valued setting, no endpoint statements are proved at Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),55 or Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),56, and it remains open whether intermediate UMD can be replaced by general UMD (Amenta et al., 2020). Abstract metric-space theorems presently yield restricted weak type in the lower range Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),57, with strong Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),58 bounds for Cf(x)=ξ∈RsupH(Mξf)(x),Mξf(x)=e2πiξxf(x),59 obtained by interpolation, while broader endpoint and full-range extensions remain natural targets (Becker et al., 2024).
The term “Carleson” also appears outside Fourier-analytic operator theory. In elliptic PDE, “Carleson perturbations” refers to Carleson-measure-controlled perturbations of elliptic coefficients, not to maximally modulated singular integrals (Feneuil et al., 2020). That usage is historically connected through Carleson measure estimates but designates a different theory. Within harmonic analysis proper, however, Carleson operators remain the prototypical maximally modulated singular integrals, and their modern generalizations continue to organize a large part of time-frequency and oscillatory analysis.