Sharp Lifespan Estimates and Fujita Phenomena for Fractional Hardy-Hénon Type Parabolic Equations

技术摘要：分数阶 Hardy–Hénon 型抛物方程的精确寿命估计与 Fujita 现象

问题陈述
本文研究了具有 Hardy–Hénon 型权重的分数阶半线性抛物方程的 Cauchy 问题：
$$\begin{cases} u_t + (-\Delta)^s u = |x|^{-\gamma} |u|^p, & (t, x) \in (0, \infty) \times \mathbb{R}^N, \\ u(0, x) = \varepsilon u_0(x), & x \in \mathbb{R}^N, \end{cases}$$
其中 $0 < s < 1$，$0 \leq \gamma < \min(2s, N)$，$p > 1$，且初始数据 $u_0 \in L^1(\mathbb{R}^N) \cap L^\infty(\mathbb{R}^N)$ 满足 $\int_{\mathbb{R}^N} u_0 \, dx > 0$。参数 $\varepsilon > 0$ 为小量。主要目标是确定解的寿命 $T_\varepsilon$（最大存在时间），并刻画爆破行为，特别是识别临界 Fujita 指数并建立亚临界、临界及超临界状态下的精确寿命估计。

方法论
作者结合了半群理论、热核估计以及针对非局部算子改进的精细测试函数法。

1. 下界（Bootstrap 论证）：

   • 寿命的下界是通过基于辅助泛函 $M(t)$ 的 Bootstrap 论证推导出的，该泛函结合了 $L^1$ 和加权 $L^\infty$ 范数。
   • 一个关键工具是加权插值不等式（引理 2.2），它通过在 $L^1$ 和 $L^\infty$ 范数之间进行插值来控制非线性项 $\||x|^{-\gamma}|u|^p\|_{L^1}$。
   • 证明过程利用了分数阶热半群 $\{e^{-t(-\Delta)^s}\}$ 精确的 $L^p$–$L^q$ 光滑估计以及分数阶热核 $G_s(t, x)$ 特有的衰减性质。

2. 上界（测试函数法）：

   • 上界是通过利用后向分数阶热核 $G_s(b-t, \cdot)$ 对方程进行测试来建立的。
   • 至关重要的一点是，该方法避免了使用紧支集截断函数，因为对于非局部算子 $(-\Delta)^s$ 而言，紧支集函数是不兼容的（紧支集函数的算子作用结果永远不是紧支集的）。
   • 通过使用全局定义的、多项式衰减的测试函数，由于 $(-\Delta)^s$ 的自伴随性，线性项被完全抵消。这使得问题简化为关于泛函 $g_b(t) = \int u(x,t) G_s(b-t, x) \, dx$ 的标量微分不等式。
   • 对于临界情况（$p = p_F$），作者采用了尺度叠加技术（在对数尺度带上进行积分）来恢复对数形式的寿命行为。

3. 爆破率与剖面：

   • 利用非局部极大值原理不等式和缩放论证分析爆破率。
   • 通过自相似变量研究渐近行为，从而导出一个涉及分数阶拉普拉斯算子和 Ornstein–Uhlenbeck 漂移项的平稳问题。

核心贡献与结果

1. 精确寿命估计（定理 1.1）：
   本文建立了小 $\varepsilon$ 情况下的精确寿命 $T_\varepsilon$：

   • 亚临界 ($1 < p < p_F$)： $T_\varepsilon \approx \varepsilon^{-1/\beta}$, 其中 $\beta = \frac{(2s-\gamma) - N(p-1)}{2s(p-1)}$。
   • 临界 ($p = p_F$)： $T_\varepsilon \approx \exp(C \varepsilon^{-(p-1)})$。
   • 超临界 ($p > p_F$)： $T_\varepsilon = +\infty$（全局存在）。
     此处，临界 Fujita 指数确定为 $p_F = 1 + \frac{2s-\gamma}{N}$。指数 $\beta$ 被证明是精确的。

2. 大数据寿命（定理 6.1）：
   对于非衰减初始数据（$u_0 \geq c_0 > 0$），本文导出了不同的寿命律：$T_\varepsilon \approx \varepsilon^{-\frac{2s(p-1)}{2s-\gamma}}$。这与小数据情形形成了对比，凸显了寿命对初始数据空间衰减特性的依赖性。

3. 爆破率（定理 7.1）：

   • 普适下界： 无条件地，任何爆破非负解都满足 $N(t) \geq C(T-t)^{-\vartheta}$, 其中 $\vartheta = \frac{2s-\gamma}{2s(p-1)}$。这是自相似速率。
   • 内部下界： 如果爆破发生在远离原点处（$|x| \geq \delta$），则速率提升至 ODE 速率 $N(t) \geq C(T-t)^{-1/(p-1)}$。
   • 有条件上界： 若解满足特定的先验自相似界，则成立 Type-I 上界 $N(t) \leq C(T-t)^{-\vartheta}$。该条件与排除 Type-II 爆破相关。

4. 渐近剖面（定理 8.1）：
   在假设满足 Type-I 界的情况下，解（沿子序列）收敛于如下平稳方程的有界分布解：
   $$(-\Delta)^s V + \frac{1}{2s} y \cdot \nabla V + \vartheta V = |y|^{-\gamma} V^p.$$
   该剖面的唯一性取决于该平稳问题的 Liouville 型分类，作者指出该分类目前仅得到部分理解。

意义与主张
作者声称，即使在 $\gamma = 0$ 的情况下（即分数阶 Fujita 方程），这些结果也是全新的，提供了关于分数阶 Lee–Ni 指数的首次精确双侧寿命估计。

• 修正缩放直觉： 本文明确反驳了“自然猜想”，即认为可以通过简单地将 $2s$ 替换为 $2s-\gamma$ 来获得临界指数或寿命速率。作者证明了 Hardy 权重是通过加性方式（通过插值不等式）进入非线性的，从而导致指数 $\beta$ 中出现特定的加性修正，而非对扩散参数进行乘性修改。
• 方法论创新： 使用后向分数阶热核作为测试函数被强调为一项显著进展。它规避了经典测试函数法（依赖紧支集）在处理非局部算子时的局限性，实现了无需人工截断即可统一处理上界的方案。
• 状态区分： 本工作澄清了原点处的爆破（受奇异权重和较慢的自相似速率 $\vartheta$ 控制）与内部爆破（受局部 ODE 动力学和较快速率 $1/(p-1)$ 控制）之间的区别。

文章最后指出，虽然下界和寿命估计是无条件的，但匹配的上界以及爆破剖面的唯一性仍取决于与该平稳问题相关的分数阶 Liouville 分类的结果，这被确定为一个开放性问题。