Weighted $W^{1,p}$-estimates for Parabolic Equations of Fabes-Kenig-Seraponi singular-degenerate type

Abstract: We investigate Dirichlet boundary value problems for a class of second-order parabolic equations in divergence-form with coefficient matrices that exhibit singular and degenerate behaviors characterized by a Muckenhoupt weight class. This framework serves as the parabolic analogue to the singular-degenerate elliptic equations pioneered by Fabes, Kenig, and Seraponi. Under a smallness assumption on the partially weighted mean oscillation of the coefficients, we establish the existence, uniqueness, and local interior and boundary regularity estimates for weak solutions within appropriately defined weighted Sobolev spaces. The proofs rely on the freezing coefficient technique alongside the level-set method introduced by Caffarelli and Peral. Additionally, we develop the necessary weighted Sobolev space framework and related weighted inequalities. Finally, a compactness argument is utilized to demonstrate that solutions to these equations remain locally close, in the weighted Sobolev norm, to their frozen-coefficient counterparts.

• The paper introduces a quantitative Calderón–Zygmund framework for obtaining weighted W¹,ᵖ estimates in parabolic equations with singular-degenerate coefficients.
• It establishes existence, uniqueness, and optimal global/local regularity of weak solutions using small BMO conditions and precise Muckenhoupt weight assumptions.
• The work opens pathways for analyzing nonlinear and multi-scale degenerate systems in heterogeneous media with rigorous, quantifiable estimates.

Weighted $W^{1,p}$-Estimates for Parabolic Equations of Singular-Degenerate Type: An Expert Analysis

Introduction and Problem Formulation

The paper "Weighted $W^{1,p}$-estimates for Parabolic Equations of Fabes-Kenig-Seraponi singular-degenerate type" (2607.09895) presents a comprehensive study of divergence-form linear parabolic equations whose principal part features singular-degenerate coefficients governed by a Muckenhoupt weight. This extends the elliptic framework of Fabes-Kenig-Serapioni (FKS) to parabolic equations, where the coefficient matrix $A(x,t)$ satisfies an ellipticity condition with a spatial weight $\beta(x)$. Specifically,

$$\nu \beta(x) |\xi|^2 \leq \langle A(x,t)\xi,\xi \rangle \leq \nu^{-1}\beta(x)|\xi|^2$$

for almost all $(x,t)$ and all vectors $\xi$, with $\beta$ in an $A_{1+1/n_0}$ Muckenhoupt class, $n_0 = \max\{n-1,1\}$.

Consideration is given to Dirichlet problems:

$W^{1,p}$0

where $W^{1,p}$1 is a vector field. The spatial degeneracy/singularity can result in both vanishing and unbounded eigenvalues, with the framework admitting weights with multiple isolated singularities.

Analytical Framework and Main Results

The analytic contribution is the establishment of a robust quantitative Calderón-Zygmund theory in parabolic weighted Sobolev spaces associated with singular-degenerate coefficients. The main statement is as follows: under a small BMO norm condition measuring the partial mean oscillation of $W^{1,p}$2 with respect to $W^{1,p}$3, and boundedness of the Muckenhoupt constant $W^{1,p}$4, for each $W^{1,p}$5, the following hold:

• Existence and uniqueness of weak solutions $W^{1,p}$6 (weighted parabolic Sobolev space with trace vanishing on the parabolic boundary),
• Global estimate: There is $W^{1,p}$7 independent of $W^{1,p}$8 such that

$W^{1,p}$9

• Interior and boundary local regularity estimates: For sufficiently small local oscillation of the coefficients, solutions enjoy local higher Lebesgue norm bounds (interior and at the flat boundary), with constants depending polynomially on the parameters.

The sharpness of the partial BMO condition is discussed: unless $A(x,t)$0, counterexamples (cf. [CMP]) confirm that this is essentially necessary.

Methodological Advances

Key technical innovations and strategies employed include:

• Weighted Sobolev and Poincaré inequalities specialized to the $A(x,t)$1 Muckenhoupt class, including new scale-invariant, weighted parabolic cylinders reflecting the degeneracy structure.
• Level-set/Caffarelli-Peral perturbative argument for Calderón-Zygmund decay: Local solutions are approximated by solutions to associated frozen-coefficient equations where boundedness and regularity are under tighter control.
• Freezing-coefficient compactness argument: The regularity theory requires demonstrating uniform closeness (in the weighted Sobolev norm) between solutions to equations with close coefficients, which is obtained via a sophisticated limiting procedure and weighted compactness/duality.
• Explicit control of degeneracy: The proofs handle quantitative estimates in function spaces where the underlying measure is singular or vanishing, requiring intricate covering, localization, and comparison techniques that adapt classical arguments (e.g., Moser-De Giorgi) to the weighted context.

Additionally, the work provides technical results such as the derivation of weighted embedding, Poincaré, and reverse Hölder inequalities needed for the global theory.

Implications and Relation to Literature

The results unify and generalize several strands in the literature:

• Connection with classical FKS theory: The interior regularity, existence, and uniqueness theorems complete the analogous parabolic program that FKS developed for elliptic equations with $A(x,t)$2 weights. The BMO assumption aligns with the necessary hypothesis for weighted Calderón-Zygmund theory for degenerate/singular operators.
• Sharpness and flexibility: The admissible weight class includes, but is not limited to, density weights of the form $A(x,t)$3 (with small $A(x,t)$4), combinations of logarithmic singularities, or even more complex products/combinations, allowing for models with multi-point degeneracy and singularity.
• Sobolev regularity in the degenerate regime: The theoretical advance is foundational for further analysis of nonlinear and nonlocal parabolic problems in highly heterogeneous or stratified media.

The approaches are robust enough to allow extensions to:

• Operators with matrix weights or systems,
• Domains with minimal geometric regularity (e.g., Reifenberg flat),
• Variants with non-divergence structure (via analogous techniques),
• Related equations in stochastic, geometric, or biological models where media heterogeneity is extreme.

Numerical and Analytical Significance

Strong claims are made: For the considered class, optimal global and local weighted $A(x,t)$5-estimates are established for the full range $A(x,t)$6. The bounds depend only polynomially on the involved parameters (weight Muckenhoupt constant and local BMO oscillation), providing a practical a priori framework for further numerical or analytical studies of degenerate/singular parabolic systems.

The paper also incorporates a wide array of examples and connections to models arising in mathematical biology (e.g., nonlinear biofilm/porous media flows), materials science (phase transitions, composite media), and geometric PDEs.

Prospective Developments

Future research directions suggested by the present framework include:

• Nonlinear extensions: Extension of these techniques to fully nonlinear or quasilinear singular-degenerate parabolic problems, where the dependence on the weighted structure is essentially nonperturbative.
• Systems and coupled equations: Addressing vectorial or multi-physics coupled cases where the weighted degeneracy is both spatially and temporally multi-scale.
• Optimality of weight class: While $A(x,t)$7 is essential for the present method, further work may elucidate whether weaker assumptions, or localized variants, could suffice.
• Integration with nonlocal and fractional operators: The weighted theory developed here interacts naturally with the burgeoning theory of nonlocal degenerate evolution equations involving fractional structure and Muckenhoupt-type weights.
• Quantitative control in minimal regularity geometry: The techniques can be adapted to sharp regularity results on non-smooth (e.g., fractal or metric measure space) domains, subject only to minimal geometric assumptions.

Conclusion

This work provides the first fully quantitative, global and local weighted Sobolev theory for divergence-form parabolic PDEs with FKS-type singular-degenerate coefficients. It synthesizes weighted harmonic analysis, parabolic regularity, and perturbation theory, and offers powerful tools for both pure and applied analysis in heterogeneous media. The methods are expected to serve as a foundation for future investigations in degenerate/singular PDEs with spatial and temporal anisotropy (2607.09895).