Calderón-Zygmund estimates for double phase problems with matrix weights

技术摘要：具有矩阵权重的双相问题的 Calderón-Zygmund 估计

问题陈述
本研究探讨了涉及退化/奇异矩阵值权重的非均匀椭圆双相问题的正则性理论。作者研究了如下在有界区域 $\Omega \subset \mathbb{R}^n$ ($n \geq 2$) 中的原型方程：
$$\text{div}\left( |MDu|^{p-2}M^2Du + a(x)|MDu|^{q-2}M^2Du \right) = \text{div}\left( |MF|^{p-2}M^2F + a(x)|MF|^{q-2}M^2F \right),$$
其中 $1 < p < q < \infty$，$a(\cdot) \in C^{0,\alpha}(\Omega)$ 是调节系数，且 $M: \mathbb{R}^n \to \mathbb{R}^{n \times n}$ 是一个对称、几乎处处正定的矩阵权重。矩阵 $M$ 满足定量各向异性界限 $|M(x)||M^{-1}(x)| \leq \Lambda$。

主要目标是建立该系统的最优 Calderón-Zygmund 理论。具体而言，作者旨在证明：如果数据项 $F$ 满足可积条件 $(|MF|^p + a(x)|MF|^q) \in L^\gamma_{\text{loc}}$（对于某个 $\gamma > 1$），那么解的梯度满足相应的估计 $(|MDu|^p + a(x)|MDu|^q) \in L^\gamma_{\text{loc}}$。该结果必须在对矩阵权重 $M$ 和非线性项进行最小正则性假设的前提下成立，同时解决潜在的“Lavrentiev 现象”（不同能量空间之间的间隙）。

方法论
证明过程采用了扰动法，并结合了多层比较策略和高级调和分析工具。该方法依赖于两个新颖的结构特征：

1. 冻结矩阵的对数：
   作者并非直接冻结矩阵 $M$，而是在每个球 $B$ 上通过其对数平均值 $(M)^{\log}_B = \exp\left(\fint_B \log M(x) dx\right)$ 来冻结 $\log M$。这种技术此前已用于单相问题，现被扩展到双相设定中。通过对数 BMO 估计控制偏差，$\log M$ 的 BMO 范数 $|\log M|_{\text{BMO}}$ 的微小性使得变量矩阵权重可以被常数矩阵近似。

2. 分数极大算子与 $A_{p,s}$ 权重：
   由 Hölder 连续系数 $a(x)$ 控制的 $q$ 相问题，需要使用加权 Lebesgue 空间中的分数极大算子。作者利用了分数 Muckenhoupt-Wheeden $A_{p,s}$ 类，其中参数满足 $1/s = 1/p - \alpha/(nq)$。该类对于控制分数极大算子的有界性以及防止在锐利阈值 $q/p \leq 1 + \alpha/n$ 处出现 Lavrentiev 现象至关重要。

证明策略与参考问题
核心论证涉及将原始解 $u$ 与定义在同心球上的四个参考问题进行嵌套比较：

1. 齐次问题： 去除源项后的原始方程。
2. 冻结矩阵问题： 将矩阵权重替换为冻结对数平均值 $(M)^{\log}_B$ 的方程。
3. 冻结系数问题： 去除非线性项的 $x$ 依赖性（在球上取平均）的方程。
4. 冻结高度问题： 在稍小的球上将调节系数 $a(x)$ 替换为常数 $a_0$ 的方程。

从变量权重 $M$ 到冻结权重的过渡由 $|\log M|_{\text{BMO}}$ 的微小性来控制。一旦权重被冻结，分析便按照类似于无权重双相问题的路径进行，确保估计与常数矩阵的大小无关。证明最终通过基于 Vitali 覆盖引理和能量密度 $H(x, MDu) = |MDu|^p + a(x)|MDu|^q$ 的退出时间估计导出的“好-$\lambda$”不等式而告终。通过迭代该不等式，可以得到所需的 $L^\gamma_{\text{loc}}$ 提升。

关键结果
主定理（定理 2.10）确立了：对于每一个 $\gamma > 1$，若满足以下条件，则蕴含关系
$$(|MF|^p + a(x)|MF|^q) \in L^\gamma_{\text{loc}} \implies (|MDu|^p + a(x)|MDu|^q) \in L^\gamma_{\text{loc}}$$
成立：

• 矩阵权重 $M$ 满足各向异性界限，且其对数的 BMO 范数足够小（$|\log M|_{\text{BMO}} \leq \delta$）。
• 非线性项满足 $\delta$-消失条件（在 $x$ 变量上的小振荡）。
• 结构参数满足锐利阈值 $q/p \leq 1 + \alpha/n$。

该结果包含了对常数 $\Lambda$、BMO 范数及数据的定量依赖关系。此外，本文证明了在此矩阵加权设定下不存在 Lavrentiev 现象，从而确保能量泛函在自然空间上的下确界与在光滑函数上的下确界一致。

意义与主张
作者声称这项工作提供了首个针对具有一般矩阵权重的双相问题的全面正则性理论。其意义在于弥合了经典双相问题（其中 $M \equiv I_n$）与由变量矩阵诱导的各向异性模型之间的鸿沟。

• 恢复经典结果： 当矩阵权重 $M$ 趋于单位矩阵（$M \to I_n$）时，结果恢复了先前工作中建立的经典（无权重）双相问题的 Calderón-Zygmund 理论（例如 [18, 20]）。
• 处理边界情况： 该方法成功处理了 $q/p = 1 + \alpha/n$ 这一边界情况——在这一机制下通常会出现 Lavrentiev 现象——通过利用分数 $A_{p,s}$ 条件与小 $\log\text{-BMO}$ 控制之间的特定相互作用。
• 泛化性： 本工作将现有的单相问题（具有变量矩阵）的框架推广到了完整的双相结构，捕捉了与强各向异性材料和均质化相关的更广泛的各向异性行为。

论文最后指出，这些技术为未来的研究方向提供了启发，包括向 Musielak-Orlicz 空间的扩展以及对边界正则性的研究，但并未详细阐述具体的实验提案或超出理论框架之外的应用。