Jacobi Heat Kernel Bounds

• Jacobi heat kernel bounds are sharp two-sided Gaussian estimates describing diffusion processes with variable coefficients and endpoint singularities.
• They leverage properties of Jacobi polynomials, explicit integral estimates, and weight-shifting techniques to track precise parameter dependencies.
• These bounds extend to diverse geometries such as symmetric spaces, balls, simplices, and conic domains, informing optimal analysis in these settings.

The Jacobi heat kernel describes the fundamental solution to the heat equation associated with the Jacobi differential operator, reflecting the geometry and analysis of weighted intervals, spheres, and related symmetric spaces. Its bounds encode the precise behavior of diffusion processes governed by variable coefficients with endpoint singularities and manifest in sharp Gaussian estimates, boundary singularities, and explicit dependence on parameters. Recent breakthroughs provide sharp, fully explicit two‐sided estimates across the full parameter range, unifying the heat kernel theory for broad families of orthogonal polynomial expansions and geometric settings.

1. Definition and Spectral Structure

The Jacobi heat kernel is defined for parameters $\alpha,\beta > -1$ on $[-1,1]$ via the weighted measure

$$d\rho_{\alpha,\beta}(x) = (1-x)^{\alpha}(1+x)^{\beta}\,dx,$$

and employs the family of Jacobi polynomials $\{P_n^{(\alpha,\beta)}\}_{n=0}^\infty$, which are eigenfunctions of the Jacobi operator

$$L_{\alpha,\beta}f(x) = -(1-x^2)f''(x) + [\beta-\alpha-(\alpha+\beta+2)x]f'(x)$$

with eigenvalues $\lambda_n = n(n+\alpha+\beta+1)$. The associated heat semigroup $e^{-tL_{\alpha,\beta}}$ admits the integral kernel

$$K_t^{(\alpha,\beta)}(x,y) = \sum_{n=0}^\infty e^{-t\lambda_n}\frac{P_n^{(\alpha,\beta)}(x)P_n^{(\alpha,\beta)}(y)}{h_n^{(\alpha,\beta)}},$$

where $h_n^{(\alpha,\beta)}$ is the $L^2$-norm squared of $P_n^{(\alpha,\beta)}$ with respect to $d\rho_{\alpha,\beta}$ (Nowak et al., 2019, Nowak et al., 2011, Nowak et al., 2024, Nowak et al., 18 Dec 2025).

2. Genuinely Sharp Two-Sided Bounds

The heat kernel exhibits fundamentally sharp two-sided Gaussian-type bounds incorporating both interior behavior and boundary singularities. For $\theta = \arccos x$, $\varphi = \arccos y$ and $0 < t \le T$, the standard two-sided bound is

$$C_{1}\,\mathcal{M}(t;\theta,\varphi) \leq K_t^{(\alpha,\beta)}(\cos\theta,\cos\varphi) \leq C_{2}\,\mathcal{M}(t;\theta,\varphi),$$

with

$$\mathcal{M}(t;\theta,\varphi) = [\sin(\tfrac{\theta}{2})\sin(\tfrac{\varphi}{2}) + t]^{-\alpha - \frac{1}{2}} [\cos(\tfrac{\theta}{2})\cos(\tfrac{\varphi}{2}) + t]^{-\beta - \frac{1}{2}}\, t^{-\frac{1}{2}} \exp\left(-\frac{(\theta-\varphi)^2}{4t}\right).$$

This structure captures:

• Gaussian off-diagonal decay at rate $e^{-d(x,y)^2/4t}$ with $d(x,y) = |\arccos x - \arccos y|$.
• Precise blow-up near boundaries, controlled by $$(\sin(\tfrac{\theta}{2})\sin(\tfrac{\varphi}{2}) + t)^{-\alpha-\frac12}$$ as $\theta,\varphi \to 0$ and $$(\cos(\tfrac{\theta}{2})\cos(\tfrac{\varphi}{2}) + t)^{-\beta-\frac12}$$ as $\theta,\varphi \to \pi$.
• Explicit dependency on $\alpha,\beta$ in all constants (Nowak et al., 2024, Nowak et al., 18 Dec 2025, Nowak et al., 2019).

For $\alpha, \beta \ge -\frac12$ and $K_t^{(\alpha,\beta)}(x, y)$ as above, alternative but essentially equivalent prefactor representations are:

$$t^{-(\alpha+\beta+2)/2} [(1-x)(1-y)]^{-\alpha/2} [(1+x)(1+y)]^{-\beta/2} \exp\left(-c\frac{d(x,y)^2}{t}\right),$$

with constants uniform for $t$ in compact intervals away from $0$ (Nowak et al., 2011). For $t \gtrsim 1$, the kernel is uniformly bounded.

3. Parametric Regimes and Explicit Constants

Recent work provides explicit tracking of all multiplicative constants in the bounds, with formulas in terms of Gamma functions, elementary functions of $\alpha,\beta$, and threshold time $T$:

• Lower constant: $$c_{\alpha, \beta, T} = [2^{\Lambda+1}\Gamma(\Lambda+1)]^{-1}$$, $\Lambda = \alpha + \tfrac12$.
• Upper constant: $C_{\alpha, \beta, T}$ combines $\Gamma(\Lambda+1)^{-1}$, powers of $2$, $\pi$, $e$, and numerical factors such as $w_0, w_0'$.

For special regimes (e.g., $\alpha, \beta \to \infty$), the ratio $C/c$ grows at most exponentially in $\alpha+\beta$. In the ultraspherical case $\alpha = \beta = \lambda - 1/2$, the bounds reduce exactly to the zonal spherical heat kernel, where constants behave as $2^d$ for $d=2\lambda$ (Nowak et al., 18 Dec 2025).

The extended product formula enables precise control for all $\alpha, \beta > -1$, not only for parameters above $-\frac12$. For $(\alpha, \beta) \in (-1, -\frac12)$, boundary terms in the product formula are controlled by analytic continuation and vanish or are quantitatively negligible (Nowak et al., 2024).

4. Extensions to Symmetric Spaces, Balls, Simplices, and Conic Domains

The Jacobi heat kernel bounds directly transfer, via transplantation/reduction, to a variety of geometric spaces:

• Compact rank-one symmetric spaces: The heat kernel satisfies

$$K_M(t;x,y) \simeq [t+\pi-d_M(x,y)]^{-(d-d_{ap}-1)/2} t^{-d/2} \exp\left(-\frac{d_M(x,y)^2}{4t}\right),$$

where $d_M(x,y)$ is the intrinsic geodesic distance (Nowak et al., 2019).

• Unit Euclidean ball: With weight $(1-|x|^2)^{\mu-\frac12}$, the kernel obeys

$$(1-|x|)^{-\mu/2}(1-|y|)^{-\mu/2} [t+ \pi - d_B(x,y)]^{-1/2} [t + d_B(x,y)]^{-\mu - 1/2} t^{-d/2} \exp\left(-\frac{d_B(x,y)^2}{4t}\right).$$

• Standard simplex: For multi-index $K$, the kernel bound is

$$\prod_{j=1}^{d+1}[t+\pi-d_j(x,y)]^{-K_j-1/2} t^{-d/2} \exp\left(-\frac{d_V(x,y)^2}{4t}\right).$$

• Conic domains: On $\mathbb{V}^{d+1}$ (solid cone) and its boundary, genuinely sharp bounds incorporate the explicit cone–geodesic distance, parameter-dependent power terms, and off-diagonal Gaussian suppression. These are among the few sharp bounds beyond classical settings (Hanrahan et al., 2022).
• Damek–Ricci spaces and noncompact rank-one symmetric spaces: The radial Jacobi heat kernel is naturally the model, with sharp two-sided estimates valid for all time and precise regime-specific asymptotics for derivatives in space and time (Bruno et al., 2022).

5. Key Methodologies and Proof Ingredients

All genuinely sharp bounds for the Jacobi heat kernel rely on a family of interconnected techniques:

• Reduction to Gegenbauer/ultraspherical kernel by the Dijksma–Koornwinder product formula; now extended via analytic continuation to all decisive parameter ranges (including when $\alpha$ or $\beta < -\frac12$) (Nowak et al., 2024, Nowak et al., 2019).
• Comparison/weight-shifting principles allow extension from dyadic/half-integer parameter sets to the full range via monotonicity and analytic iteration (Nowak et al., 2019, Nowak et al., 2011).
• Explicit integral estimates for oscillatory and singular-weighted integrals deliver uniform optimal powers and handle boundary singularities (Nowak et al., 2011, Nowak et al., 2024).
• Iteration and gluing via the semigroup property provide bounds at arbitrary time scales, beginning from sharp small-time control (Nowak et al., 18 Dec 2025).
• Spectral sum manipulation for conic and higher-dimensional cases, often via repeated reduction to the one-dimensional Jacobi kernel (Hanrahan et al., 2022).

The hierarchical reduction from general orthogonal polynomials or geometric settings to the classical Jacobi heat kernel is a central structural paradigm.

6. Applications and Theoretical Significance

Sharp Jacobi heat kernel bounds have broad consequences:

• Maximal operator bounds: For multi-dimensional settings, the associated maximal operator $T_*$ enjoys sharp weak-type $(1,1)$ norm control, uniform in $\alpha, \beta \ge -\frac12$ (Nowak et al., 2011).
• Sharp spectral multiplier and Riesz transform theory: Optimal kernel estimates underlie $L^p$ bounds for spectral multipliers, endpoint inequalities, and Littlewood-Paley theory (Bruno et al., 2022).
• Boundary phenomena, norm convergence, and transference: Quantification of singularity rates at the endpoints establishes the regularity and norm-convergence of Jacobi polynomial expansions (Nowak et al., 2019).
• Transfer to geometric PDE: Through transplantation, results on balls, spheres, simplices, and cones inform the fundamental understanding of heat flows, eigenfunction expansions, and boundary regularity in a wide class of compact or weighted manifolds and domains.

7. Recent Advances and Open Problems

Advancements since 2019 include:

• Extension of genuinely sharp two-sided bounds to all $(\alpha,\beta)>-1$ (Nowak et al., 2024).
• Complete tracking and explicit characterization of all multiplicative constants, including their dimension and parameter dependence (Nowak et al., 18 Dec 2025).
• Unified proof strategies applicable to spherical settings and associated spaces, refining and extending bounds to encompass geometric transitions, endpoint behaviors, and full parameter regimes.
• For conic domains and higher-dimensional analogues, genuinely sharp Gaussian bounds have been established, combining radial and angular features (Hanrahan et al., 2022).

A plausible implication is the stabilization of the theory: all classical settings with rank-one, symmetric, or polynomial structure now possess explicit kernel bounds. Further subtlety may arise in multidimensional generalizations and for non-polynomial structures, but for all classical settings the optimal form is known and quantitatively precise.