Hölder regularity of harmonic functions on metric measure spaces

Technical Summary: Hölder Regularity of Harmonic Functions on Metric Measure Spaces

Problem Statement
The paper addresses the regularity of harmonic functions and heat kernels on unbounded metric measure Dirichlet (MMD) spaces, specifically focusing on spaces arising from strongly recurrent fractals (e.g., Sierpiński gasket, Sierpiński carpet). While classical De Giorgi-Nash-Moser theory establishes the local Hölder continuity of weak solutions to elliptic equations on Riemannian manifolds and Euclidean domains, it typically yields an implicit Hölder exponent dependent on structural constants rather than an explicit value.

On fractals, the behavior of the heat kernel is governed by sub-Gaussian estimates involving a "walk dimension" $\beta$ strictly greater than the Hausdorff dimension $\alpha$. Previous work by Devyver, Russ, and Yang (2023) established gradient estimates for heat kernels on certain cable systems (graphs with edges) under a "generalized reverse Hölder inequality" (GRH). However, this approach relied on harmonic extension algorithms intrinsic to p.c.f. (post-critically finite) self-similar sets, failing to apply to the Sierpiński carpet, and left the validity of GRH on the Sierpiński carpet cable system as an open problem. Furthermore, existing methods often required explicit gradient-type conditions or probabilistic frameworks (fractional diffusions) to derive regularity.

Methodology
The authors introduce a new Hölder regularity condition, denoted $HR(\Phi, \Psi)$, tailored for metric measure spaces where explicit exponents can be derived. The methodology relies on the following framework:

1. Geometric and Analytic Assumptions: The paper assumes a volume regular condition $V(\Phi)$ where volume scales as $\Phi(r) \approx r^\alpha$ (with different exponents $\alpha_1, \alpha_2$ for small and large scales) and an upper heat kernel estimate $UHK(\Psi)$ with a time-distance scaling function $\Psi(r) \approx r^\beta$. Crucially, the authors assume a strongly recurrent condition ($\alpha_i < \beta_i$), which distinguishes these spaces from transient ones and ensures the existence of a resistance form.
2. Equivalence of Conditions: The core technical achievement is proving the equivalence of five distinct conditions under the assumptions of $V(\Phi)$ and $UHK(\Psi)$:
   • The Hölder regularity condition $HR(\Phi, \Psi)$ for harmonic functions.
   • The weak Bakry-Émery non-negative curvature condition $wBE(\Phi, \Psi)$.
   • Hölder continuity of the heat kernel $HHK(\Phi, \Psi)$.
   • Hölder continuity of the heat kernel with exponential terms $HHK_{exp}(\Phi, \Psi)$.
   • The near-diagonal lower bound for the heat kernel $NLE(\Psi)$.
3. Proof Techniques:
   • Poisson Equation Approach: To prove $HR \Rightarrow wBE$ and $HR \Rightarrow HHK_{exp}$, the authors treat the heat equation $\Delta(P_t f) = \partial_t(P_t f)$ as a Poisson equation. They utilize existence and uniqueness results for Poisson equations on unbounded MMD spaces, combined with pointwise estimates derived from the heat kernel bounds.
   • Resistance Form Theory: To prove the reverse implication $NLE \Rightarrow HR$, the authors leverage the theory of resistance forms (Kigami). They establish that under the strongly recurrent condition, the Dirichlet form induces a resistance form. Using effective resistance estimates and a maximum principle for harmonic functions in the extended Dirichlet space, they derive the Hölder regularity without assuming the existence of a gradient operator (carré du champ).
   • Analytic vs. Probabilistic: The authors explicitly avoid the probabilistic framework of fractional diffusions used in prior works (e.g., Barlow), opting for a concise analytic approach.

Key Contributions and Results

1. Equivalence Theorem (Theorem 2.1): The paper establishes that for unbounded MMD spaces satisfying volume regularity and upper heat kernel bounds, the Hölder regularity of harmonic functions is equivalent to the near-diagonal lower bound of the heat kernel. This generalizes previous gradient-type equivalences (e.g., Coulhon et al.) to a Hölder-type setting without requiring a gradient operator.
2. Resolution of the Sierpiński Carpet Open Problem: By proving that two-sided heat kernel estimates ($HK(\Psi)$) imply the Hölder regularity condition $HR(\Phi, \Psi)$, and subsequently the generalized reverse Hölder inequality $GRH(\Phi, \Psi)$, the authors resolve the open problem left by Devyver, Russ, and Yang regarding the validity of $GRH$ on the Sierpiński carpet cable system.
3. Gradient Estimates without Gradient Assumptions: The paper demonstrates that on strongly recurrent fractal-like cable systems, pointwise two-sided heat kernel estimates alone are sufficient to derive gradient estimates for the heat kernel. This improves upon previous results that required assuming $GRH$ or specific gradient-type conditions.
4. Explicit Hölder Exponents: Unlike classical iteration techniques that provide implicit exponents, this framework yields explicit Hölder exponents for harmonic functions and heat kernels: $\beta_1 - \alpha_1$ at small scales and $\beta_2 - \alpha_2$ at large scales.
5. Generalization to Unbounded Spaces and Blowups: The results extend to unbounded spaces with potentially different scaling behaviors at small and large scales. The authors apply these findings to "blowup" constructions, where one fractal is scaled by another (e.g., blowing up the Vicsek set by the Sierpiński carpet), showing that the resulting spaces exhibit scale-dependent gradient behaviors (sub-Gaussian at small scales, Gaussian or sub-Gaussian at large scales).

Significance
The paper claims significance in three primary areas:

• Theoretical Unification: It provides a unified analytic framework linking heat kernel estimates, curvature conditions (in the weak Bakry-Émery sense), and Hölder regularity for a broad class of fractal spaces, extending classical Riemannian manifold results (Li-Yau, Cheeger-Yau) to the metric measure Dirichlet setting.
• Stability: The authors argue that their Hölder-type conditions are quasi-isometry stable, unlike many gradient-type conditions which often require curvature assumptions that are not stable under quasi-isometry. This stability stems from the strongly recurrent nature of the spaces.
• Resolution of Open Problems: It definitively settles the validity of the generalized reverse Hölder inequality on the Sierpiński carpet cable system, a case where previous methods relying on harmonic extension algorithms failed. Furthermore, it removes the need for a priori gradient assumptions to obtain heat kernel gradient estimates, relying solely on heat kernel bounds and volume regularity.

The work is presented as a rigorous extension of regularity theory to non-smooth, fractal-like geometries, offering explicit regularity exponents and resolving specific technical gaps in the analysis of cable systems associated with complex fractals.