Gaussian estimates for fundamental solutions of higher-order parabolic equations with time-independent coefficients

1. Problem Statement

The paper addresses the regularity and pointwise behavior of fundamental solutions for higher-order parabolic operators in divergence form with time-independent, complex-valued, and measurable coefficients.

Specifically, the author studies the operator $H$ defined on $\mathbb{R}^{n+1}$:
$$H := \partial_t + L := \partial_t + \sum_{|\alpha|=|\beta|=m} (-1)^{|\alpha|} \partial^\alpha (a_{\alpha,\beta}(x) \partial^\beta)$$
where $m \geq 1$ is the order of the operator (so the total order is $2m$). The coefficients $a_{\alpha,\beta}(x)$ are:

• Bounded and measurable.
• Complex-valued.
• Time-independent (depend only on $x$).
• Uniformly elliptic in the sense of the Gårding inequality.

The central goal is to establish Gaussian upper bounds and Hölder regularity estimates for the fundamental solution $\Gamma(t, x, s, y)$ of $H$, particularly in the regime where the order of the operator satisfies $2m \geq n$.

2. Methodology

The paper employs a multi-step strategy combining elliptic regularity theory, perturbation arguments, and semigroup methods:

• Transfer from Elliptic to Parabolic: The core methodology relies on the observation that if the associated elliptic operator $L$ satisfies specific Hölder continuity properties (De Giorgi-Moser-Nash estimates), then the parabolic operator $H$ (with time-independent coefficients) inherits similar properties.
• Technical Lemmas and Energy Estimates: The author constructs a series of technical lemmas (Section 3) involving:
  • Caccioppoli-type inequalities for higher-order parabolic equations.
  • Iteration arguments to derive uniform estimates on time slices.
  • Poincaré-Sobolev inequalities for higher-order elliptic settings.
  • Morrey space estimates to handle inhomogeneous terms.
• Perturbation Theory: The paper utilizes a small $L^\infty$ perturbation argument. It establishes that if a "base" operator satisfies Hölder estimates, a sufficiently small perturbation of its coefficients also satisfies these estimates. This allows the extension of results from operators with constant or VMO (Vanishing Mean Oscillation) coefficients to general bounded measurable coefficients.
• Semigroup and Resolvent Analysis: To prove Gaussian bounds without assuming a priori Hölder continuity, the author uses the analytic semigroup generated by $L$. By analyzing the resolvent $(z+L)^{-1}$ and utilizing Sobolev embeddings (specifically when $2m \geq n$), the author derives pointwise bounds for the heat kernel $K_t(x,y)$ of $L$.
• Fundamental Solution Construction: The fundamental solution $\Gamma$ for the parabolic equation is constructed via the heat kernel of the elliptic operator: $\Gamma(t, x, s, y) = K_{t-s}(x, y)$ for $t > s$.

3. Key Contributions

A. Extension of De Giorgi-Moser-Nash Estimates

The paper proves that if the elliptic operator $L$ satisfies the Hölder property $H(\mu, \infty)$ (a higher-order version of De Giorgi-Moser-Nash estimates), then the corresponding parabolic operator $H$ with time-independent coefficients satisfies a similar parabolic Hölder property $H(\mu_0, \infty)$ for some $\mu_0 < \mu$.

• Theorem 4.1: Establishes the transfer of Hölder regularity from elliptic to parabolic equations.
• Theorem 4.3: Provides a perturbation theorem, showing that the Hölder property is stable under small $L^\infty$ perturbations of the coefficients. This generalizes previous results by Auscher (for second-order) to higher-order systems.

B. Gaussian Upper Bounds for Fundamental Solutions

The most significant contribution is the derivation of Gaussian upper bounds for the fundamental solution $\Gamma(t, x, s, y)$ under the condition $2m \geq n$.

• Theorem 5.5: Proves that for any parabolic operator $H$ with time-independent coefficients in the class $E_{2m}$, if $2m \geq n$, the fundamental solution satisfies Gaussian estimates without requiring the coefficients to be a small perturbation of real or constant coefficients.
• The estimates include pointwise bounds on derivatives $D_x^\gamma \Gamma$ and $D_y^\gamma \Gamma$, as well as Hölder continuity in the spatial variables.

C. Handling Complex Coefficients

The results hold for complex-valued coefficients, which is a significant generalization over classical real-coefficient theories. The paper demonstrates that the principal part of the operator dictates the kernel's behavior, and lower-order terms (or complex perturbations) do not destroy the Gaussian bounds provided the Gårding inequality holds.

4. Main Results

1. Hölder Regularity Transfer: Weak solutions to the parabolic equation $H u = 0$ possess Hölder continuity if the associated elliptic operator $L$ possesses it. Specifically, if $L$ satisfies $H(\mu, \infty)$, then $H$ satisfies $H(\mu_0, \infty)$.
2. Stability under Perturbation: The class of operators satisfying Hölder estimates is open with respect to the $L^\infty$ norm of the coefficients. This implies that operators with VMO coefficients or small perturbations of complex constants satisfy these estimates.
3. Gaussian Estimates (Theorem 5.5): For $2m \geq n$, the fundamental solution $\Gamma(t, x, s, y)$ satisfies:
   $$|D_x^\gamma \Gamma(t, x, s, y)| + |D_y^\gamma \Gamma(t, x, s, y)| \leq \frac{C}{(t-s)^{\frac{n+|\gamma|}{2m}}} g_{c}\left( \frac{|x-y|}{(t-s)^{1/2m}} \right)$$
   where $g_c(s) = e^{-c s^{\frac{2m}{2m-1}}}$.
   Furthermore, Hölder continuity estimates for the derivatives are established:
   $$|D_x^\gamma \Gamma(t, x+h, s, y) - D_x^\gamma \Gamma(t, x, s, y)| \leq \frac{C |h|^\nu}{(t-s)^{\frac{n+|\gamma|+\nu}{2m}}} g_c(\dots)$$
   for $|\gamma| = l$ (where $l$ is an integer related to the regularity).

5. Significance and Impact

• Generalization of Second-Order Theory: The paper successfully extends the celebrated results of Auscher, Hofmann, Kim, and others (originally for second-order parabolic equations) to higher-order systems ($2m \geq 4$).
• Removal of Small Perturbation Assumption: Unlike earlier works that required coefficients to be small perturbations of real or constant coefficients, this paper proves Gaussian bounds for general time-independent complex coefficients provided $2m \geq n$. This is a crucial step toward understanding the behavior of fundamental solutions in the most general measurable setting.
• Foundation for Boundary Value Problems: The author notes that these Gaussian estimates are essential for developing higher-order layer potential theory. This theory is the primary tool for solving boundary value problems (Dirichlet, Neumann, Regularity problems) for higher-order parabolic operators in non-smooth domains.
• Complex Coefficients: By handling complex coefficients, the results apply to a broader range of physical and mathematical models, including those involving damping or non-self-adjoint operators.

In summary, Guoming Zhang's work provides a rigorous analytical framework establishing that the fundamental solutions of higher-order parabolic equations with time-independent complex coefficients exhibit Gaussian decay and Hölder regularity, provided the order of the equation is sufficiently high relative to the spatial dimension ($2m \geq n$). This bridges a critical gap between elliptic regularity theory and parabolic fundamental solution estimates.