Dunkl Paraproduct Operators

• Dunkl paraproduct operators are bilinear operators built via Littlewood–Paley decomposition using compactly supported Dunkl multipliers to split products into resonant and paraproduct components.
• They ensure Lᵖ boundedness and satisfy bilinear Calderón–Zygmund estimates, which are critical for deriving fractional Leibniz rules in the Dunkl analytic framework.
• The approach extends classical harmonic analysis to reflection-group geometries, offering practical tools for nonlinear analysis and the construction of Sobolev and Besov spaces.

Dunkl paraproduct operators are bilinear operators built from a Littlewood–Paley decomposition associated with the Dunkl transform, designed to decompose products $fg$ into interactions of comparable, low–high, and high–low Dunkl frequencies. In the formulation developed in "Dunkl paraproducts and fractional Leibniz rules for the Dunkl Laplacian" (Bui et al., 14 Jul 2025), they take the form

$$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$

where the kernels are defined through compactly supported Dunkl multipliers and the Dunkl convolution $*_k$. Their principal role is to provide the harmonic-analytic mechanism by which fractional Leibniz rules for the Dunkl Laplacian $\Delta_k$ are established without assuming $G$-invariance of either factor (Bui et al., 14 Jul 2025).

1. Dunkl-analytic setting

The theory is formulated on $\mathbb{R}^d$ relative to a finite root system $R\subset\mathbb{R}^d$, with $|\lambda|=\sqrt{2}$ for all $\lambda\in R$. Each $\lambda\in R$ determines the reflection

$$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$0

and the associated reflection group $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$1 is generated by these reflections. A multiplicity function is a $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$2-invariant map $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$3. Writing

$$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$4

the weighted density is

$$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$5

and the measure is

$$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$6

The parameter $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$7 acts as a homogeneous dimension in the volume-growth estimates for $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$8.

For $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$9, the Dunkl operator is the differential–difference operator

$*_k$0

and the Dunkl Laplacian is

$*_k$1

When $*_k$2, these reduce to the Euclidean derivative and Laplacian. A basic structural obstruction is that Dunkl operators do not satisfy the classical Leibniz rule in general; the identity

$*_k$3

holds only when at least one factor is $*_k$4-invariant.

The associated transform theory is organized around the Dunkl kernel $*_k$5, the unique analytic solution of

$*_k$6

For $*_k$7, the Dunkl transform is

$*_k$8

with $*_k$9 a homeomorphism, a Plancherel identity on $\Delta_k$0, and the intertwining relation

$\Delta_k$1

Dunkl translation is defined spectrally by

$\Delta_k$2

and Dunkl convolution by

$\Delta_k$3

A central technical difficulty is that, unlike Euclidean translation, the global $\Delta_k$4 boundedness of $\Delta_k$5 is known only for radial functions; the paper identifies this lack of full translation invariance as a key obstacle in the construction of a paraproduct theory (Bui et al., 14 Jul 2025).

2. Definition of Dunkl paraproducts and product decomposition

The Dunkl paraproduct construction begins with $\Delta_k$6 and Schwartz functions $\Delta_k$7 satisfying

$\Delta_k$8

For $\Delta_k$9, one sets

$G$0

so that $G$1, and similarly for $G$2. The operator $G$3 therefore functions as a Dunkl frequency projector at scale $G$4.

With these ingredients, the Dunkl paraproduct is

$G$5

Its heuristic content matches Bony’s Euclidean construction: both inputs are first localized in Dunkl frequency, then multiplied at matched or separated scales, and finally smoothed by $G$6.

The product decomposition is obtained from a partition of unity $G$7 with

$G$8

Writing $G$9 in Dunkl frequency variables and splitting the indices by relative scale yields

$\mathbb{R}^d$0

where $\mathbb{R}^d$1 is resonant, $\mathbb{R}^d$2 is low–high, and $\mathbb{R}^d$3 is high–low.

Term	Scale relation	Support pattern
$\mathbb{R}^d$4	$\mathbb{R}^d$5	both inputs at medium frequency
$\mathbb{R}^d$6	$\mathbb{R}^d$7	$\mathbb{R}^d$8 at scale $\mathbb{R}^d$9, $R\subset\mathbb{R}^d$0 low frequency
$R\subset\mathbb{R}^d$1	$R\subset\mathbb{R}^d$2	$R\subset\mathbb{R}^d$3 low frequency, $R\subset\mathbb{R}^d$4 at scale $R\subset\mathbb{R}^d$5

More precisely, Proposition 4.3 rewrites each $R\subset\mathbb{R}^d$6 as a paraproduct $R\subset\mathbb{R}^d$7 with explicit support conditions. For the resonant term,

$R\subset\mathbb{R}^d$8

For the low–high term,

$R\subset\mathbb{R}^d$9

and for the high–low term the roles of $|\lambda|=\sqrt{2}$0 and $|\lambda|=\sqrt{2}$1 are reversed. The support bookkeeping relies on a Dunkl convolution lemma stating that convolution of an annular cutoff with a much smaller ball-supported cutoff remains in a slightly enlarged annulus (Bui et al., 14 Jul 2025).

3. $|\lambda|=\sqrt{2}$2-boundedness and bilinear Calderón–Zygmund structure

The central boundedness theorem states that if $|\lambda|=\sqrt{2}$3 and the supports of at least two among $|\lambda|=\sqrt{2}$4 do not contain $|\lambda|=\sqrt{2}$5, then for

$|\lambda|=\sqrt{2}$6

the paraproduct maps

$|\lambda|=\sqrt{2}$7

boundedly, with estimate

$|\lambda|=\sqrt{2}$8

The proof proceeds in two stages. A preliminary proposition yields boundedness for exponents $|\lambda|=\sqrt{2}$9, $\lambda\in R$0, using Littlewood–Paley methods and vector-valued inequalities. The passage to the larger bilinear Calderón–Zygmund range is then obtained by showing that the integral kernel of $\lambda\in R$1 satisfies bilinear Dunkl–Calderón–Zygmund size and regularity conditions, after which a bilinear CZ theorem in the Dunkl setting is invoked.

The kernel representation has the form

$\lambda\in R$2

Its analysis is governed by the orbit distance

$\lambda\in R$3

rather than purely Euclidean separation. This reflects the geometry of the reflection group and the nonlocality of the Dunkl translation. The paper does not formulate a Coifman–Meyer theorem directly in Fourier variables; instead, the Coifman–Meyer-type content is encoded through the support and smoothness of $\lambda\in R$4 and the corresponding kernel estimates (Bui et al., 14 Jul 2025).

4. Fractional Leibniz rules generated by paraproducts

The principal motivation for introducing Dunkl paraproduct operators is the derivation of fractional Leibniz rules for the Dunkl Laplacian. Fractional powers are defined spectrally by

$\lambda\in R$5

The main estimate is: if

$\lambda\in R$6

then for $\lambda\in R$7,

$\lambda\in R$8

A mixed-order version is also proved: if $\lambda\in R$9 with $\lambda\in R$0, then

$\lambda\in R$1

The paraproduct method is structural. First, one decomposes $\lambda\in R$2. Second, one uses the support configuration of the cutoffs to move the multiplier $\lambda\in R$3 onto a single factor whenever possible. For the low–high and high–low terms this gives identities of the form

$\lambda\in R$4

and

$\lambda\in R$5

For the resonant term,

$\lambda\in R$6

where $\lambda\in R$7 is a variant of the standard paraproduct because $\lambda\in R$8 need not be Schwartz.

The low–high and high–low pieces are controlled by the paraproduct boundedness theorem. The resonant piece requires additional decay and maximal-function arguments. This is why, although the paraproduct theory itself reaches a larger bilinear range, the fractional Leibniz rules in the paper are established only for $\lambda\in R$9. When $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$00, these results recover the classical Kato–Ponce inequalities (Bui et al., 14 Jul 2025).

5. Auxiliary estimates underpinning the theory

Three auxiliary mechanisms support the paraproduct analysis.

First, the paper proves pointwise decay for fractional Dunkl Laplacians of Schwartz functions: $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$01 The proof uses the heat-semigroup representation

$$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$02

together with Gaussian bounds for the kernel $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$03 of $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$04: $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$05 A corollary states that for $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$06 and $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$07, one has $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$08 for all $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$09, and also for certain $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$10.

Second, an almost orthogonality estimate controls triple products of decaying factors in the Dunkl metric. If $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$11, then

$$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$12

is bounded by

$$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$13

This is the Dunkl counterpart of the almost orthogonality lemmas used in bilinear Euclidean Calderón–Zygmund theory.

Third, because the kernel representation of a paraproduct involves Dunkl translations of Schwartz functions, the paper establishes pointwise size and Hölder-type continuity bounds for $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$14. If $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$15 is smooth and supported in $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$16, $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$17, and $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$18, then for $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$19,

$$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$20

and when $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$21,

$$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$22

These three ingredients combine to verify the bilinear CZ kernel conditions and to control the modified resonant paraproduct by the Dunkl Hardy–Littlewood maximal operator $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$23, using the domination

$$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$24

in the resonant analysis (Bui et al., 14 Jul 2025).

6. Relation to classical theory, prior results, and scope

The structure of Dunkl paraproducts mirrors classical Bony decompositions: compactly supported multipliers define scale-localized pieces, the product is split into resonant and paraproduct components, and boundedness is obtained through Littlewood–Paley and bilinear singular-integral arguments. Their analytical role is also classical in form: they reduce fractional Leibniz estimates to bounds where the derivative falls on one factor at a time.

The differences are substantial. The measure $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$25 is not translation invariant, the relevant distance is $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$26 rather than Euclidean distance alone, and the Dunkl translation operator is not known to be globally bounded on all $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$27. For this reason, the Euclidean bilinear multiplier approach based on the symbol $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$28 is not transferred directly. The paper explicitly avoids treating $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$29 as a bilinear multiplier in $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$30 and instead insists on a paraproduct decomposition that always allows $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$31 to act on a single factor.

In the classical limit $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$32, one has $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$33, $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$34, and the Dunkl transform becomes the Fourier transform, so the theory reduces to the usual Kato–Ponce framework. Within Dunkl analysis, the paper generalizes a previous result of Wrȯbel, which treated a restricted case with $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$35 and one $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$36-invariant factor. The new work allows arbitrary finite reflection groups $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$37, removes the assumption that one factor is $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$38-invariant, and states that it corrects an error in Wrȯbel’s argument.

The paper also identifies broader consequences. Mixed-order Leibniz rules extend higher-order fractional product estimates to the Dunkl Laplacian, and the authors note that such inequalities are building blocks for Sobolev and Besov spaces defined via $$\Pi[\theta,\psi,\phi](f,g)(x) = \sum_{j\in\mathbb{Z}} \Theta_j *_k \Big((\Psi_j *_k f)(\cdot)\,(\Phi_j *_k g)(\cdot)\Big)(x),$$39, for nonlinear dispersive or parabolic equations with Dunkl operators, and for nonlinear functional calculus in Dunkl harmonic analysis. This suggests that Dunkl paraproduct operators provide a reusable framework for nonlinear analysis in reflection-group geometries, rather than a device limited to a single product estimate (Bui et al., 14 Jul 2025).