Mixed modulation spaces are function spaces defined by measuring the short-time Fourier transform in mixed phase-space norms to capture time-frequency concentration.
They include variants such as permutation-dependent, quasi-Banach, and α-modulation spaces that refine classical modulation frameworks through anisotropic and order-sensitive norms.
These spaces facilitate discrete decompositions via Gabor frames and underpin sharp kernel theorems and operator characterizations in advanced time-frequency analysis.
Mixed modulation spaces are time–frequency function spaces defined by measuring the short-time Fourier transform (STFT) in mixed phase-space norms. In the standard setting, for a nonzero window g, the STFT is
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,
and the classical modulation space Mp,q(Rd) is obtained by requiring Vgf∈Lp,q(R2d). In broader usage, the subject includes permutation-dependent spaces M(c)p1,…,p2d, multi-index mixed quasi-norm spaces Mwp, mixed-norm α-modulation spaces, and modulation spaces associated to tensor products of amalgam spaces; all of them retain the same basic principle that time–frequency concentration is quantified by anisotropic or order-sensitive mixed norms on phase space (Cordero et al., 2017, Toft, 2014, Han et al., 2011, Feichtinger et al., 2020).
1. Foundational definitions and the STFT framework
The core object throughout the theory is the STFT, which identifies a distribution f with a phase-space function Vgf(x,ξ). In the classical mixed-norm case,
∥f∥Mp,q=(∫Rd(∫Rd∣Vgf(x,ξ)∣pdx)q/pdξ)1/q,
with the usual modifications when Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,0 or Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,1. The window Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,2 can be changed without changing the space, provided Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,3 is chosen in the admissible class specified in the relevant theory; this window-independence is a basic structural property in both the Banach and quasi-Banach settings (Cordero et al., 2017, Cordero et al., 2023, Toft, 2014).
A more flexible construction replaces the two exponents Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,4 by a full multi-index and, in some formulations, a permutation of the phase-space coordinates. In the permutation-based setting of kernel theorems, one fixes a permutation Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,5 of Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,6, lets Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,7 reorder the coordinates, and defines
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,8
When Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,9 is the identity and all exponents coincide, one recovers the standard space Mp,q(Rd)0; when the exponents split into spatial and frequency blocks, one recovers Mp,q(Rd)1. A distinguished permutation Mp,q(Rd)2 satisfies
Mp,q(Rd)3
so the mixed modulation formalism also contains Wiener amalgam spaces (Cordero et al., 2017).
The quasi-Banach generalization replaces Mp,q(Rd)4 by iterated mixed quasi-norms Mp,q(Rd)5, with Mp,q(Rd)6, possibly after a prescribed permutation of coordinates. In that setting,
Mp,q(Rd)7
and the analysis is formulated in Gelfand–Shilov classes when the weights are more general than polynomially moderate. This extends the classical Mp,q(Rd)8 scale to a genuinely multi-parameter mixed-norm theory (Toft, 2014).
2. Variants of mixed modulation spaces
Across the literature, mixed modulation spaces appear in several structurally distinct but closely related forms. One major variant is the Mp,q(Rd)9-modulation scale Vgf∈Lp,q(R2d)0, Vgf∈Lp,q(R2d)1, built from an Vgf∈Lp,q(R2d)2-covering of frequency space. With associated projections Vgf∈Lp,q(R2d)3, the norm is
Vgf∈Lp,q(R2d)4
For Vgf∈Lp,q(R2d)5 one recovers modulation-space geometry, while Vgf∈Lp,q(R2d)6 is treated by convention as the Besov space Vgf∈Lp,q(R2d)7. The space is therefore a mixed decomposition scale interpolating geometrically between uniform and dyadic frequency tilings (Han et al., 2011).
A further anisotropic refinement is the mixed-norm Vgf∈Lp,q(R2d)8-modulation space
Vgf∈Lp,q(R2d)9
where M(c)p1,…,p2d0. Here the mixedness is spatially anisotropic: each coordinate direction may carry a different integrability exponent, while the global M(c)p1,…,p2d1-summation remains tied to the M(c)p1,…,p2d2-dependent frequency covering (Nielsen, 2023).
Another variant starts from a Banach space M(c)p1,…,p2d3 on phase space and defines
M(c)p1,…,p2d4
When M(c)p1,…,p2d5, one recovers M(c)p1,…,p2d6. When M(c)p1,…,p2d7 is a completed projective or injective tensor product of Wiener amalgam spaces, the resulting modulation space can itself be identified with a Wiener amalgam space. For example, one of the paper’s principal special cases is
M(c)p1,…,p2d8
together with the injective analogue
M(c)p1,…,p2d9
This identifies a broad family of apparently different mixed modulation constructions with concrete amalgam spaces (Feichtinger et al., 2020).
The special affine Fourier transform generates yet another realization. The SAFT-based spaces Mwp0 are defined through Mwp1-convolution and Mwp2-modulation, but the key structural theorem is that they are chirp-conjugated classical modulation spaces: Mwp3
This shows that the SAFT theory does not introduce an unrelated scale; it re-encodes mixed modulation spaces in a transform-adapted form (Biswas et al., 2022).
3. Discretization, Gabor frames, and molecular decompositions
Discretization by Gabor analysis is one of the defining technical mechanisms of the subject. If Mwp4 is a Gabor frame with dual window Mwp5, then the coefficient and synthesis maps
Mwp6
satisfy Mwp7, and the resulting coefficient sequence yields an equivalent discrete description of the modulation norm. For mixed modulation spaces Mwp8, the Gabor coefficients characterize the norm in mixed Mwp9 spaces, with unconditional reconstruction when all exponents are finite (Cordero et al., 2017).
This discretization persists in the broad quasi-Banach setting. For α0, the analysis operator α1 and synthesis operator α2 extend continuously between modulation spaces and the corresponding mixed sequence spaces α3, and dual Gabor frames provide the reconstruction formula
α4
with unconditional convergence when α5 and weak-α6 convergence otherwise. This result places mixed-norm modulation spaces within a fully reconstructible frame theory rather than a merely abstract norm definition (Toft, 2014).
For mixed-norm α7-modulation spaces, the analogous tool is the α8-transform, which produces a discrete sequence space α9 and a tight frame adapted to the f0-covering. The boundedness of the analysis map f1, the boundedness of the synthesis map f2, and the identity f3 provide a discrete decomposition theory parallel to Frazier–Jawerth constructions (Nielsen, 2023).
The molecular formulation refines this further. A family f4 is a system of f5-molecules if it obeys simultaneous spatial and frequency localization estimates of the form
f6
f7
The change-of-frame matrix between two such molecular systems is almost diagonal, and the associated almost diagonal matrices form an algebra under composition. This algebraic closure is the basis for perturbation theory, compactly supported frame constructions, and multiplier estimates in mixed-norm f8-modulation spaces (Nielsen, 2023).
4. Kernel theorems and operator-theoretic characterizations
The sharpest operator-theoretic use of mixed modulation spaces is the kernel theorem for operators on f9. Let Vgf(x,ξ)0 be a linear continuous operator on Vgf(x,ξ)1 with distribution kernel Vgf(x,ξ)2. For Vgf(x,ξ)3,
Vgf(x,ξ)4
and
Vgf(x,ξ)5
From these endpoint statements one deduces
Vgf(x,ξ)6
The result completely characterizes bounded operators on Vgf(x,ξ)7 by placing the kernel in carefully permuted mixed modulation spaces (Cordero et al., 2017).
The theorem is quantitative: the operator norm is equivalent to the corresponding kernel norm. In particular,
Vgf(x,ξ)8
The proof identifies the kernel with a Gabor matrix
Vgf(x,ξ)9
and translates operator boundedness into mixed ∥f∥Mp,q=(∫Rd(∫Rd∣Vgf(x,ξ)∣pdx)q/pdξ)1/q,0 estimates for this matrix. Feichtinger’s kernel theorem appears as the special case
∥f∥Mp,q=(∫Rd(∫Rd∣Vgf(x,ξ)∣pdx)q/pdξ)1/q,1
so the mixed formulation is a genuine extension of the earlier theorem (Cordero et al., 2017).
A persistent point of interpretation is that the clean iff-characterization is specific to the diagonal scale ∥f∥Mp,q=(∫Rd(∫Rd∣Vgf(x,ξ)∣pdx)q/pdξ)1/q,2. For operators on ∥f∥Mp,q=(∫Rd(∫Rd∣Vgf(x,ξ)∣pdx)q/pdξ)1/q,3, the same paper states that a full kernel characterization analogous to the ∥f∥Mp,q=(∫Rd(∫Rd∣Vgf(x,ξ)∣pdx)q/pdξ)1/q,4 case is not expected. Instead, it proves sufficient conditions such as
∥f∥Mp,q=(∫Rd(∫Rd∣Vgf(x,ξ)∣pdx)q/pdξ)1/q,5
∥f∥Mp,q=(∫Rd(∫Rd∣Vgf(x,ξ)∣pdx)q/pdξ)1/q,6
and, by interpolation,
∥f∥Mp,q=(∫Rd(∫Rd∣Vgf(x,ξ)∣pdx)q/pdξ)1/q,7
implies boundedness on ∥f∥Mp,q=(∫Rd(∫Rd∣Vgf(x,ξ)∣pdx)q/pdξ)1/q,8 for all ∥f∥Mp,q=(∫Rd(∫Rd∣Vgf(x,ξ)∣pdx)q/pdξ)1/q,9. This distinction is central: mixed modulation spaces provide robust kernel criteria beyond the classical kernel theorem, but the diagonal and off-diagonal cases behave differently (Cordero et al., 2017).
The same kernel perspective underlies recent boundedness theory for Fourier integral operators. In the study of Schrödinger-type propagators, the kernel
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,00
is estimated in permutation-dependent mixed modulation spaces such as Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,01, while the phase derivatives are measured in working spaces of Wiener amalgam type. The resulting theorems give sharp Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,02-boundedness thresholds in low, mild, critical, and high growth regimes for the phase, with necessity shown by discrete embedding counterexamples. This suggests that mixed modulation spaces have become a primary language for sharp time–frequency kernel estimates of oscillatory operators (Guo et al., 7 Jul 2025).
5. Embeddings, symmetries, and comparison with neighboring scales
Embedding theory makes explicit how mixed norms interact with smoothness and frequency geometry. For classical modulation spaces Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,03, one may write
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,04
which isolates two different mechanisms: local Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,05-control on uniform frequency boxes and global Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,06-summation over the box index. The sharp comparison with Besov, Triebel–Lizorkin, Sobolev, and Fourier Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,07 spaces is governed by the thresholds
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,08
These determine the exact smoothness needed for embeddings into modulation spaces and the maximal smoothness compatible with embeddings out of modulation spaces (Lu, 2021).
Within the Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,09-modulation scale, the corresponding sharp embedding criterion is
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,10
where
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,11
The same paper proves exact complex interpolation and duality, including the formula
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,12
It also shows that the heuristic identification of Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,13-modulation spaces as interpolation spaces between modulation and Besov scales fails in general. This is a notable correction to an oversimplified view of the subject (Han et al., 2011).
Symmetry properties are equally sensitive to mixed structure. For a metaplectic operator Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,14 projecting to Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,15, the covariance formula
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,16
reduces boundedness on Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,17 to boundedness of the pullback Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,18 on Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,19. The precise criterion is that Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,20 is bounded if and only if either Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,21, or Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,22 and Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,23 is upper triangular. The mixed-norm geometry is therefore not invariant under arbitrary symplectic changes of variables when Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,24; only the diagonal case Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,25 enjoys the full metaplectic flexibility stated in that theorem (Cordero et al., 2023).
The same metaplectic viewpoint yields alternative norm characterizations. Shift-invertible metaplectic Wigner distributions Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,26 satisfy
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,27
under the upper-triangularity hypothesis on the relevant linear map, with the triangularity condition dropping in the case Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,28. This places STFTs, Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,29-Wigner distributions, and related phase-space representations within a common mixed-norm invariance theory (Cordero et al., 2023).
6. Operator classes, compactness, and analytical applications
The mixed-norm viewpoint has extensive operator-theoretic consequences. On the functional-analytic side, weighted mixed-norm Lebesgue spaces Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,30 have the metric approximation property under the factorization hypothesis
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,31
and this transfers to weighted modulation spaces Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,32 via Gabor frame representations. As a consequence, Grothendieck’s nuclearity theory applies, and an operator Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,33 is Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,34-nuclear precisely when its kernel admits a decomposition
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,35
with
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,36
For Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,37, the associated trace formula is
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,38
This ties mixed modulation norms directly to trace and spectral theory (Delgado et al., 2014).
Compactness questions can also be phrased in mixed modulation terms. For localization operators Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,39 acting on weighted spaces Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,40 of Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,41-tempered distributions, boundedness follows from Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,42, while compactness follows from a vanishing-at-infinity condition on the STFT of the symbol. In particular, if Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,43, then
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,44
is compact for Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,45. The proof passes through a kernel expansion using a tight Gabor frame and then approximates general symbols by test symbols in Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,46 (Boiti et al., 2019).
Spectral multiplier theory provides another major application. For the Hermite operator
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,47
the multiplier Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,48 is bounded on Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,49 under a localized Sobolev condition of order Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,50 when
Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,51
and on the diagonal spaces Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,52 under the sharper condition Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,53. The same work states that these mixed-norm conditions are not sharp in the off-diagonal case and uses Heisenberg-group transference and torus multiplier theory to prove the diagonal improvement. It also deduces boundedness of Hermite Riesz transforms in the same mixed range and shows that the wave and Schrödinger propagators for Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,54 preserve modulation-space regularity (Bhimani et al., 2017).
Nonlinear dispersive theory likewise exploits mixed modulation spaces. For the mixed fractional Hartree equation and the Hartree equation with harmonic potential, local and global well-posedness are established in Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,55 and in Fourier amalgam spaces Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,56, based on trilinear Hartree estimates and Strichartz estimates. The paper emphasizes that this yields Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,57- and Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,58-regularity for all Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,59 in the stated regimes, extending earlier Sobolev-based results and admitting low-regularity data (Bhimani et al., 2023).
A final analytical direction is variational. The uncertainty-principle framework based on modulation spaces formulates minimization problems on constraint sets defined by mixed-norm modulation spaces Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,60, proves a mixed-norm extension of Lieb’s inequality for ambiguity functions, and identifies the optimal constant with an extremal eigenvalue of an inverse compact localization operator in the Hilbert case. This connects mixed modulation norms to compact embeddings, localization operators, and Euler–Lagrange equations, including the harmonic-oscillator ground-state equation in the classical one-dimensional case (Dias et al., 2022).
Taken together, these developments show that mixed modulation spaces are not merely a notational extension of Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,61. They form a family of STFT-based spaces in which anisotropy, coordinate ordering, frequency geometry, and operator structure can be encoded with high precision. The literature also makes clear that this precision comes with sharp constraints: full kernel theorems are currently specific to certain scales, metaplectic invariance is restricted in the genuinely mixed case Vgf(x,ξ)=∫Rdf(t)g(t−x)e−2πiξ⋅tdt,62, and interpolation heuristics can fail outside special parameter regimes (Cordero et al., 2017, Cordero et al., 2023, Han et al., 2011).