Decay rates to equilibrium in a nonlinear subdiffusion equation with two counteracting terms

技术摘要：具有两个相互抵消项的非线性亚扩散方程向平衡态的衰减率

问题陈述
本文研究了具有两个相互抵消项的非线性亚扩散方程解的长时渐近行为。该模型由以下初边值问题给出：
$$\partial_t^\alpha u - Lu = q(x)u - p(x)f(u) + r(t,x)$$
其中 $\partial_t^\alpha$ 表示阶数为 $\alpha \in (0, 1]$ 的 Djirbashian-Caputo 分数阶导数（当 $\alpha=1$ 时代表标准抛物型情形），$L$ 是一个椭圆算子，$p, q$ 是空间依赖系数。项 $q(x)u$ 和 $-p(x)f(u)$ 起着相反的作用，这一特征在诸如种群动力学、化学反应和相分离等应用中十分常见。

主要目标是量化解 $u(t)$ 随时间 $t \to \infty$ 向稳态 $u_\infty$ 收敛的过程。具体而言，本文旨在建立差值 $u - u_\infty$ 及其时间导数的衰减率（对于 $\alpha=1$ 为指数衰减，对于 $\alpha \in (0,1)$ 为幂律衰减），其条件允许系数 $p$ 和 $q$ 具有低正则性（$p \in L^r, q \in L^s$）且具有空间变化性。这一分析受到通过间接观测重建系数 $p(x)$ 和 $q(x)$ 所使用的不动点方案收敛性证明需求的驱动。

方法论
作者结合使用了能量估计、分数阶微积分工具和渐近分析技术：

1. 线性化与差值公式化： 差值 $u_\sim = u - u_\infty$ 被表征为另一个新的初边值问题的解。非线性项使用带有积分余项的泰勒公式进行展开，使得方程可以写成包含势函数 $p_\infty(x) = p(x)f'(u_\infty) - q(x)$ 的线性化算子的形式。
2. 能量估计： 分析的核心涉及用 $u_\sim$ 测试差分方程，以及在更高正则性要求下用 $\mu u_\sim - L u_\sim$ 进行测试。这些测试利用了分数阶导数的强制性（具体而言是估计式 $v \partial_t^\alpha v \geq \frac{1}{2} \partial_t^\alpha [v^2]$）以及算子 $L$ 的椭圆性。
3. 处理非线性和低正则性： 为了适应低正则性系数（$p \in L^r$），作者利用 Sobolev 嵌入和广义 Hölder 不等式仔细估计了非线性余项。这对于非线性 $f$ 的增长（受 $f''$ 控制）相对于维度 $d$ 和 $p$ 的可积性所施加的约束提出了特定要求。
4. Gronwall 型与 Tauberian 估计： 推导出的微分不等式通过涉及 Mittag-Leffler 函数的分数阶 Gronwall 型引理进行分析。对于渐近衰减率，作者利用 Tauberian 定理将解的拉普拉斯变换在零点附近的性质与时间域中的长时行为联系起来。
5. 屏障方法（Barrier Methods）： 为了处理能量估计中由于非线性产生的负项，采用了屏障方法。这确保了如果初始数据足够接近稳态（或者满足全局符号条件），解将保持在由线性稳定性占主导地位的机制内。

主要贡献与结果

• 低正则性系数的衰减率： 本文确立了对于 $\alpha \in (0,1)$ 的幂律衰减（$O(t^{-\alpha})$）以及对于 $\alpha=1$ 的指数衰减，针对差值 $\|u(t) - u_\infty\|_{H^s}$，其中 $s \in \{0, 1, 2\}$。至关重要的是，这些结果对于处于 $L^r$ 空间中的系数 $p, q$ 依然成立，而此前依赖于高正则性或常系数的结论并未涵盖这一设定。
• 时间导数的衰减： 作者证明了时间导数 $\partial_t^\alpha u$（以及 $\alpha=1$ 时的 $u_t$）在适当的范数下也趋于零。这是前言中所提应用的批判性结果，因为它是证明用于逆问题的固定点算子收缩性的必要条件。
• 局部稳定性与全局稳定性：
  • 局部稳定性（定理 2.2）： 在稳态处线性化势函数满足符号条件（$p(x)f'(u_\infty) - q(x) \geq c > 0$）且初始数据较小时，证明了衰减性。
  • 全局稳定性（定理 2.6）： 通过要求在整个轨迹上（或在瞬态与平衡态之间）满足势函数的符号条件而非仅仅在稳态处，去除了对初始数据大小的限制，从而得到了全局衰减结果。
• 正则性要求： 论文规定了在维度 $d=1, 2, 3$ 下，为了使估计成立，关于 $f$ 的增长和 $p$ 的可积性所需的精确条件。例如，在 3D 中，对于三次非线性（Allen-Cahn 型），$p \in L^6$ 即已足够。

意义与声明
本文声称填补了文献中关于具有竞争项和空间变化系数的非线性亚扩散方程向稳态收敛量化研究的空白。虽然之前的研究（如 [18]）处理过类似问题，但它们通常要求 $p \equiv 1$ 或 $q \leq 0$（非抵消项），或者依赖于需要系数具有更高空间正则性的变分公式。

作者明确指出，他们基于能量估计而非变分公式的方法，使得他们能够在处理逆问题所需的低正则性系数设定下开展工作。他们也承认了一些局限性，例如由于缺乏适用于 $\alpha < 1$ 的分数阶导数特定强制性估计（无法使用奇异函数进行测试），导致无法在衰减估计中移除范数上的平方。此外，他们指出对于 $\alpha < 1$，初始时间导数 $\|u_t(0)\|_{L^2}$ 的有限性无法仅从 PDE 数据中通过提升（bootstrapping）获得，必须作为先验假设。

综上所述，这项工作为具有相互抵消项的非线性亚扩散方程的稳定性与衰减提供了一个严谨的分析框架，将其扩展到了与实际逆问题相关的更一般的低正则性系数设定中。