生成时间: 2026年08月20日 03:30:52

运行状态: 完成

计划处理论文数量: 31

已完成论文数量: 31

未完成论文数量: 0

重点关注论文数量: 9

了解领域论文数量: 14

不相关论文数量: 8

重点关注论文（完整分析）

1. 关于热扩散系统解的一致界

On the uniform bound of solutions to a thermo-diffusive system

作者: Kirill Leontev, Lenya Ryzhik

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.19111v1

摘要翻译: 我们证明了具有 $L^ ext{ extbackslash{}infty}$ 初值的热扩散反应 - 扩散方程组解的全局时间一致界。

Comment: 无

链接: http://arxiv.org/abs/2608.19111v1

1. 研究对象和背景

本文研究的是全空间 $\mathbb{R}^d$ 上的热扩散（thermo-diffusive）反应扩散方程组，其数学模型由以下耦合抛物型方程组描述：

$$ \begin{cases} \partial_t u = \nu \Delta u - u g(v), \\ \partial_t v = \kappa \Delta v + u g(v), \end{cases} $$

其中 $u(t, x)$ 表示燃料浓度，$v(t, x)$ 表示燃料温度，$\nu > 0$ 和 $\kappa > 0$ 分别为燃料和温度的扩散系数。非线性反应项 $g(v)$ 满足 $g(0)=0$，且在 $v>0$ 时非负。初始条件假设为有界函数：存在常数 $K \ge 1$ 使得 $0 \le u_0(x) \le K$ 且 $0 \le v_0(x) \le K$。该模型在燃烧理论中具有重要地位，描述了燃料消耗与热量释放之间的相互作用。

问题的核心难点在于：当扩散系数满足 $\kappa \ge \nu$（即温度扩散快于或等于燃料扩散）时，传统的能量估计或简单的比较原理难以直接导出解的全局一致有界性。对于单一的非线性热方程 $\partial_t v = \Delta v + v^p$，解可能在有限时间内发生爆破（blow-up）。虽然当 $\nu \ge \kappa$ 时，通过 Martin-Pierre 技巧可以证明 $v(t,x)$ 的有界性，但在 $\kappa \ge \nu$ 的情况下，由于热源项 $u g(v)$ 的正反馈机制，温度 $v$ 的爆炸风险显著增加。尽管已有文献（如 Collet & Xin [6]）针对多项式增长的非线性给出了对数级的时间依赖上界，而 Herrero et al. [13] 证明了指数增长情形下的全局存在性但未给出显式一致界，Chen & Qin [5] 及后续工作 [15] 则在特定假设（如 Fisher-KPP 型非线性、特定的初始条件或点火温度效应）下建立了统一界。然而，这些结果要么依赖于特殊的非线性结构（如 $g(v)=v$），要么要求初始条件满足“触发效应”（hair trigger effect）的强条件（如 $u_0+v_0 \ge K^{-1}$），或者无法处理燃烧理论中常见的点火非线性（ignition nonlinearity，即 $g(v)$ 在低温区间为零）。本文旨在填补这一空白：在仅假设初始条件有界（无衰减要求）、非线性 $g(v)$ 满足一般的增长条件（包括点火型非线性）以及 $\kappa \ge \nu$ 的条件下，严格证明温度 $v(t, x)$ 存在一个与时间无关的全局一致上界。

2. 主要定理或主要结果

本文的主要结论是定理 1.1，其断言如下：

设 $u(t, x)$ 和 $v(t, x)$ 是方程组 (1.1) 的非负解，初始条件 $u_0, v_0$ 满足 (1.2)，非线性项 $g(v)$ 满足 (1.3)-(1.5)。则存在常数 $C_0 > 0$（仅依赖于维度 $d$、扩散系数 $\nu, \kappa$、初始界 $K$ 以及非线性 $g$ 的具体形式），使得对所有 $t > 0$ 和 $x \in \mathbb{R}^d$，成立：

$$ 0 < v(t, x) \le C_0. $$

该结论的关键技术细节包括：

  1. 适用范围：适用于任意维数 $d \ge 1$，且不需要初始条件在无穷远处衰减。
  2. 非线性假设：涵盖了多项式增长和指数增长的非线性，特别是包含了点火型非线性（$g(v)=0$ for $v \le \theta$），这是此前基于“触发效应”的理论所无法处理的。
  3. 扩散系数关系：重点解决了 $\kappa \ge \nu$ 这一困难情形，此时温度扩散快于燃料扩散，传统技巧失效。
  4. 一致性：所得界 $C_0$ 不随时间 $t$ 增长，彻底消除了先前工作中出现的 $\log t$ 或 $\log \log t$ 的时间依赖性。

证明的核心步骤归结为引理 2.5，即证明热半群平均量 $[S_\kappa(\tau)v(t, \cdot)](x)$ 对所有 $t, x$ 一致有界：

$$ [S_\kappa(\tau)v(t, \cdot)](x) \le C_\tau, $$

结合抛物 BMO 范数估计（Proposition 2.2）和 John-Nirenberg 不等式（Proposition 2.3），这一平均量的有界性直接蕴含了点态一致界。

3. 研究方法、关键技术和核心工具

本文的证明框架极具创新性，主要克服了两大技术障碍：一是 $\kappa \ge \nu$ 时无法直接利用热核大小关系控制源项；二是缺乏“触发效应”导致无法保证燃料的快速耗尽。作者引入了三个关键的分析工具：

  1. 修正的 Martin-Pierre 技巧（Modified Martin-Pierre Trick）： 针对 $\kappa \ge \nu$ 的情形，原技巧失效。作者在 Lemma 2.1 中构造了一个非局部的版本。通过引入时间步长 $\tau$，将温度方程中的强迫项 $q = ug(v)$ 的热半群平均与燃料浓度的热半群平均联系起来。具体地，证明了：

    $$ [S_\kappa(\tau)\tilde{W}(t, t-\tau, \cdot)](x) \le M_{\kappa,\nu} [S_\kappa(2\tau)u(t-\tau, \cdot)](x), $$

    其中 $M_{\kappa,\nu} = ((2\kappa-\nu)/\kappa)^{d/2}$。这一估计利用了热核的缩放性质 $S_a(t)f \le (a/b)^{d/2} S_b(t)f$，并通过插值算子 $S_1(a_s)$ 将不同扩散系数的热核进行比较。这使得即使温度扩散更快，也能通过燃料浓度的平均量来约束温度方程的源项积分。

  2. 燃料耗竭的前向时间估计（Forward-in-time Fuel Decay Estimate）： 为了闭合论证，必须证明在高温区域燃料浓度会迅速下降。Lemma 3.1 指出，如果在某个抛物圆柱体 $Q^-_{\ell_\eta}(t, x)$ 内，温度超过点火阈值 $\vartheta$ 的区域占比足够大（即大部分区域被“点燃”），那么经过时间 $\tau_\eta$ 后，燃料浓度将按因子 $\eta$ 衰减：

    $$ u(t, x) \le \eta [S_\kappa(\tau_\eta)u(t-\tau_\eta, \cdot)](x). $$

    这一引理替代了以往文献中依赖的“触发效应”，它不要求全局的小扰动，而是依赖于局部高温区域的平均效应，从而能够处理点火非线性。

  3. 离散时间 Duhamel 展开与冷热分解（Discrete Time Duhamel Expansion and Hot-Cold Decomposition）： 这是证明中最精妙的部分。作者将时间离散化为步长 $\tau_\eta$，对 Duhamel 公式进行迭代展开。定义“热集” $U_t = \{x : (v)_{Q^-_{\ell_\eta}(t,x)} > Z_\eta\}$，并将演化算子分解为作用在热集上的算子 $T_n$ 和冷集上的算子 $F_n$。 在热集上，利用 Lemma 3.1 获得燃料的指数级衰减，从而抑制源项；在冷集上，利用抛物 BMO 理论和 John-Nirenberg 不等式获得温度的先验有界性（$v \le M_\eta$）。通过反复迭代，利用算子 $T_n$ 和 $F_n$ 的性质以及 Duhamel 展开中的关键抵消项（cancellation in the Duhamel expansion），最终证明了平均量 $V(t, x)$ 的一致有界性。这种冷热交替的迭代策略有效地平衡了局部高温导致的燃料消耗与整体温度控制的矛盾。

4. 与之前工作的比较

本文与既有文献相比，在多个维度上取得了实质性推进：

  1. 扩散系数条件的放宽： 早期的经典结果（如 Martin-Pierre [16]）仅在 $\nu \ge \kappa$（燃料扩散快于温度）时能直接给出一致界。Collet & Xin [6] 和 Herrero et al. [13] 虽然处理了 $\kappa \ge \nu$ 的情况，但给出的界依赖于时间（$\log \log t$ 或 $\log t$）。本文首次在不限制 $\nu \ge \kappa$ 的前提下，证明了 $\kappa \ge \nu$ 时的全局一致界，填补了扩散系数相对大小关系的理论缺口。

  2. 非线性假设的普适性： Chen & Qin [5] 的结果高度依赖于具体的线性非线性 $g(v)=v$ 以及特定的初始条件 $u_0 \equiv 1$。随后 Liouville 类型的工作 [15] 放宽了非线性到 $g'(0)>0$，但仍要求初始条件满足 $u_0+v_0 \ge K^{-1}$ 以触发“hair trigger effect”。本文完全去除了对初始条件下界的限制，允许 $u_0, v_0$ 仅为有界函数，甚至可以在某些区域为零。更重要的是，本文明确处理了点火型非线性（ignition nonlinearity），即 $g(v)=0$ 当 $v \le \theta$，这类非线性在物理上对应真实的燃烧过程，且因缺乏自动触发机制，使得 [15] 中的方法完全失效。本文通过新的燃料耗竭估计成功克服了这一障碍。

  3. 时间尺度的改进： 前人工作 [6, 13] 得到的界随时间增长，意味着随着时间推移，温度可能变得任意大（尽管增长缓慢）。本文证明的界 $C_0$ 是绝对常数，与时间 $t$ 无关，这在数学分析上是一个质的飞跃，表明系统具有更强的稳定性。

  4. 技术代价与局限： 本文的代价在于证明过程的复杂性显著增加。不再依赖简单的代数不等式或直接的比较原理，而是需要构建精细的迭代框架、引入热半群的平均估计以及复杂的冷热集合分解。此外，虽然结论是一致界，但常数 $C_0$ 的具体显式表达式并未给出，它依赖于非线性 $g$ 的复杂性质（如点火温度 $\theta$ 和增长率）。尽管如此，本文的方法论为解决更广泛的耦合反应扩散系统的正则性问题提供了强有力的新范式。

2. 无穷远处的完全刚性性与Levinson空腔的存在性

Complete Rigidity at infinity and Existence of the Levinson Cavity

作者: Yan Li, Sitan Lin, Georg S. Weiss, Chunjing Xie

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18913v1

摘要翻译: 我们提出一种位势论方法，将自由表面渐近形状的分析简化为对由该化简过程导出的一个精确常微分方程的分析。尽管该方法主要依赖于偏微分方程算子的主部以允许表示公式的导出，因此并不局限于椭圆型问题，但我们以三维轴对称稳态不可压缩空腔流动这一清晰范例展示该方法；此类问题属于Neumann型Bernoulli自由边界问题，其频率公式尚属未知，即便存在也不足以得出此处所证明的极为精确的渐近行为。

Comment: 无

链接: http://arxiv.org/abs/2608.18913v1

1. 研究对象和背景

本文研究的是三维轴对称定常不可压缩空腔流动问题，这是一个经典的自由边界问题。在物理上，当流体相对于刚体的高速运动时，在低压力区域会形成空腔，其边界（自由面）的形状是未知的，必须作为解的一部分确定。对于无粘、无旋的流动，速度势 $\phi$ 满足如下 Neumann 型 Bernoulli 自由边界问题：

$$ \begin{cases} \Delta \phi = 0 & \text{在流域 } D \subset \mathbb{R}^3 \text{ 内}, \\ \partial_\nu \phi = 0 & \text{在固定边界 } N \cup \Gamma \text{ 上}, \\ |\nabla \phi| = 1 & \text{在自由边界 } \Gamma \text{ 上}, \end{cases} $$

其中 $N$ 为固定边界，$\Gamma$ 为自由边界，$\nu$ 和 $s$ 分别为 $\partial D := N \cup \Gamma$ 上的单位外法向量和单位切向量。采用柱坐标 $(t, \rho, \lambda)$，其中 $t$ 为轴向坐标，$\rho$ 为径向距离，$\lambda$ 为方位角。假设固定边界和自由边界由连续函数 $g: [1, \infty) \to [0, \infty)$ 参数化：

$$ N := \{(t, g(t), \lambda) : t \in [1, a], \lambda \in [0, 2\pi)\}, \quad \Gamma := \{(t, g(t), \lambda) : t \in (a, \infty), \lambda \in [0, 2\pi)\}. $$

该问题的核心难点在于自由边界 $\Gamma$ 的渐近形状。早在 1946 年，Norman Levinson 通过幂律试探解加上缓慢变化的修正项，形式推导出了空腔形状的精确渐近公式：

$$ g(r) \sim C r^{1/2} (\log r)^{-1/4} \left[ 1 - \frac{1}{8} \frac{\log \log r}{\log r} + O\left(\frac{1}{\log r}\right) \right] \quad \text{当 } r \to \infty. $$

然而，Levinson 的推导依赖于极强的假设（如 $g(r)$ 在无穷远处凹性、特定的幂律形式 $r^k h(r)$ 以及导数比值的极限条件），且部分积分在 $k \in (3/4, 1)$ 范围内不收敛，导致其严格性存疑。尽管 Garabedian-Lewy-Schiffer [12] 利用变分法证明了此类轴对称无限空腔解的存在性，Caffarelli-Friedman [5] 进一步证明了空腔具有次线性增长速率，但关于这些已知存在的解是否真的具备 Levinson 所预测的精确渐近行为，或者至少是否存在一个满足该渐近行为的解，这一悬而未决的问题已持续了 80 年。

此外，Neumann 型 Bernoulli 问题在数学分析中 notoriously hard（极其困难），即使是在轴对称设定下，传统的频率公式（frequency formulas）或单调性公式似乎并不存在或不足以导出如此精确的渐近行为。Dirichlet 型 Bernoulli 问题（自由面为解的水平集）已有大量研究，但 Neumann 型问题缺乏“平坦性蕴含正则性”等关键工具。本文旨在填补这一空白：在温和的自然假设下，证明任何轴对称解在无穷远处都严格收敛于 Levinson 剖面，从而确立该剖面的“完全刚性”，并基于此结合存在性定理证明具有 Levinson 渐近行为的解确实存在。

2. 主要定理或主要结果

本文的主要结果建立在以下三个自然假设之上：

  1. 次线性增长：函数 $g$ 在无穷远处满足 $\lim_{t \to \infty} \frac{g(t)}{t} = 0$。
  2. 流速限制：固定边界 $N$ 上的流速不超过自由边界 $\Gamma$ 上的流速，即 $|\nabla \phi| \le 1$ on $N$。
  3. 渐近水平速度：当 $|X| \to \infty$ 时，流速趋于水平，即 $\lim_{|X| \to \infty} \nabla \phi(X) = (1, 0, 0)$。

在此假设框架下，作者证明了以下两个主要定理：

定理 1 (上界与下界)：假设上述 (1)-(3) 成立，则对于任意 $\epsilon > 0$，有

$$ \limsup_{r \to \infty} \frac{g(r)}{r^{1/2}(\log r)^{-1/4+\epsilon}} < \infty. $$

如果进一步假设固定边界 $N$ 是凹的（朝向空腔区域），即 $g'' \le 0$ on $(1, a)$，则对于任意 $\epsilon > 0$，有

$$ \liminf_{r \to \infty} \frac{g(r)}{r^{1/2}(\log r)^{-1/4-\epsilon}} > 0. $$

这意味着在凹性条件下，解的渐近行为被夹逼在 $r^{1/2}(\log r)^{-1/4}$ 的量级附近。

定理 2 (存在性)：结合 Garabedian-Lewy-Schiffer [12] 和 Caffarelli-Friedman [5] 的存在性结果与定理 1，证明了存在三维轴对称空腔流动，其自由边界满足精确的渐近不等式：

$$ r^{1/2}(\log r)^{-1/4-\epsilon} < g(r) < r^{1/2}(\log r)^{-1/4+\epsilon} $$

对每个 $\epsilon > 0$ 及充分大的 $r$ 成立。这直接确认了具有 Levinson 渐近行为的解的存在性，解决了长达 80 年的开放问题。

值得注意的是，定理 1 的证明不依赖能量最小化或有限 Morse 指标假设，仅基于 PDE 本身的性质和几何约束。

3. 研究方法、关键技术和核心工具

本文的核心创新在于提出了一种**位势理论约化至常微分方程（Potential Theoretic Reduction to ODE）**的分析框架，成功绕过了传统频率公式缺失的障碍。该方法并非泛泛的能量估计，而是通过精细的积分表示和迭代分析，将自由边界的渐近形状分析转化为一个精确的二阶常微分方程的定性分析。

1. 广义位势核与积分恒等式的构造 传统的格林函数在无穷远处发散，无法直接用于渐近分析。作者构造了一个混合类型的广义位势核（既非广义牛顿位势，也非单层或双层位势），并利用格林定理将速度势 $V$（定义为速度势与其 blow-down 极限之差）在自由边界点 $Y=(r, g(r), 0)$ 处进行积分表示。这导出了一个精确的积分恒等式，连接了速度势与自由边界轮廓 $g$：

$$ V(Y) = \int_{\partial D} G(X, Y) \dots dS_X, $$

其中 $G$ 是精心选择的广义核。这一步骤的关键在于选择正确的奇点邻域（transition window）大小。作者发现简单的 $\eta$-邻域不够，最终确定了过渡窗口为 $(r/2, 2r)$，并将积分分解为过渡窗口积分 $J_g$、尾部积分 $T_g$ 和头部积分 $H_g$。

2. 迭代估计与积分不等式的推导 证明过程具有显著的“两阶段/迭代”特征：

  • 第一阶段（初始上界）：通过对各分项积分进行粗略估计，作者首先建立了 $g(r) \le r^{1/2+\epsilon}$ 的初始上界。虽然系数 $1/9$ 不够 sharp，但这为后续细化提供了基础。
  • 第二阶段（精细估计）：利用初始上界控制尾部积分 $T_g$，并对过渡窗口积分 $J_g$ 进行更精细的拆分。特别是 Lemma 5.4 中，通过 Hardy 不等式和细致的局部估计，获得了 $g(r)g'(r)\log(r/g(r))$ 项的 sharp 系数，消除了非局部项 $V(r)$ 的干扰（Corollary 5.7）。这一过程克服了非局部算子带来的技术障碍。

3. 扰动积分等式与 Tauberian 步骤 经过反复迭代，作者导出了关键的扰动积分等式（Proposition 6.8）：

$$ V(r) + (1+o(1))g(r)g'(r)\log\frac{r}{g(r)} - \frac{1}{2}\int_a^r \frac{g(t)g'(t)}{t} dt = C + o(1), \quad r \to \infty. $$

该等式形式上对应于二阶 ODE：

$$ (g(r)g''(r) + g'(r)^2)\log\frac{r}{g(r)} + \frac{g(r)g'(r)}{2r} + \sqrt{1+g'(r)^2} - 1 - g'(r)^2 = 0. $$

为了从积分等式严格推导出 ODE 的渐近行为，作者引入了Tauberian 步骤。这是证明中的关键一环，用于弥合扰动积分等式与未扰动的 ODE 之间的差距，排除了另一个临界对数指数 $-1/2$ 的可能性，最终锁定唯一的渐近解为 $g(r) \sim r^{1/2}(\log r)^{-1/4}$。

整个证明框架避免了频率公式，转而依赖位势展开的奇点分析和 ODE 的动力学机制，展示了在处理无单调性公式的自由边界问题时，位势理论与积分不等式结合的强大约束力。

4. 与之前工作的比较

本文与既有文献相比，在多个维度上实现了重大推进：

1. 解决长期开放问题 vs. 形式推导 Levinson [20] 在 1946 年通过幂律试探解和渐近平衡分析给出了精确公式，但其推导依赖于极强的假设（如 $g(r)=r^k h(r)$ 且 $rh'/h \to 0$），且部分积分在关键参数范围内不收敛，因此其结论在 80 年间未被严格证明。Garabedian-Lewy-Schiffer [12] 和 Caffarelli-Friedman [5] 虽然证明了存在性，但仅能给出次线性增长的上界，无法触及具体的对数修正项。本文首次严格证明了 Levinson 公式的正确性，并证明了在温和假设下所有解都具有该渐近行为，填补了从“存在性”到“精确渐近性”的巨大鸿沟。

2. 方法论突破：绕过频率公式 在 Dirichlet 型 Bernoulli 问题（如障碍物问题）中，研究者常利用频率公式（frequency formula）或单调性公式来分析奇异点的 blow-down 极限（如 Eberle-Figalli-Weiss [7]）。然而，本文明确指出，对于 Neumann 型 Bernoulli 空腔问题，频率公式似乎不存在或不适用。Eberle-Figalli-Weiss [7] 的方法依赖于广义牛顿位势在对称轴上的展开，只能得到上界估计，难以处理下界。本文创新性地使用了混合位势核，并在自由边界点上展开，通过迭代积分估计和 Tauberian 理论，成功构建了无需频率公式的完整刚性证明。这种方法论不仅适用于椭圆型问题，作者还推测可推广至抛物型（如 Stefan 问题）甚至混合型方程。

3. 假设条件的放宽与适用范围 既有工作往往需要解满足特定的能量极小化或有限 Morse 指标假设才能讨论渐近行为。本文的刚性定理（Theorem 1）完全不依赖这些假设，仅要求自然的物理条件（次线性增长、流速限制、渐近水平性）和固定边界的凹性。这使得结论适用于更广泛的物理场景。此外，本文明确区分了“严格证明的结论”与“启发式解释”，确保了数学严谨性，而之前的数值模拟（如 Garabedian [11]）虽与 Levinson 公式吻合，但缺乏理论支撑。

4. 代价与局限 本文方法的代价是计算极其繁琐，涉及多层积分分解和非局部项的精细抵消，且目前仅限于轴对称情形。作者承认，将其推广至非轴对称情况（需考虑内外轴对称包络或平均值）以及带有涡度的流动（需满足特定结构假设）将是未来极具挑战性但可能的工作。尽管如此，本文提供的“位势约化至 ODE”框架为处理其他缺乏单调性公式的自由边界问题提供了全新的范式。

3. 一类拟线性抛物系统的局部零可控性及相关数值实验

Local null controllability of a quasi-linear system and related numerical experiments

作者: Enrique Fernández-Cara, Juan Límaco, Yuri Thamsten, Denilson Menezes

类别: math.OC, math.AP, math.NA

发布日期: 2026-08-19

arXiv ID: 2608.19023v1

摘要翻译: 本文研究扩散系数依赖于状态变量梯度的拟线性抛物系统的零控制问题。在我们的主要理论结果中，在关于扩散系数的正则性和增长性以及初始数据正则性的若干假设下，我们证明了局部零可控性成立。为此，我们考虑了线性化系统的零可控性问题，导出了关于控制量和状态的新估计，然后应用局部逆定理。我们还构造了一种准牛顿型的迭代算法，用于计算一个零控制量及其关联的状态。我们将此方法应用于该问题的若干数值逼近，并通过多个实验展示了结果。

Comment: 31 pages, 16 figures, 9 tables

链接: http://arxiv.org/abs/2608.19023v1

1. 研究对象和背景

本文研究一类扩散系数依赖于状态变量梯度的拟线性抛物系统的内部零控制问题。具体而言，考虑定义在 $d$ 维空间区域 $\Omega \subset \mathbb{R}^d$ ($1 \le d \le 3$) 上的非线性热传导方程：

$$ \begin{cases} y_t - \nabla \cdot [a(\nabla y)\nabla y] = \chi_\omega v & \text{in } Q := ]0, T[ \times \Omega, \\ y = 0 & \text{on } \Sigma := ]0, T[ \times \partial \Omega, \\ y|_{t=0} = y_0 & \text{in } \Omega. \end{cases} $$

其中 $y(t,x)$ 表示温度分布，$v$ 为作用在子区域 $\omega \Subset \Omega$ 上的外部热源（控制项），$\chi_\omega$ 是 $\omega$ 的特征函数。物理上，该模型描述了热导率 $a$ 随温度梯度变化的非牛顿流体或具有变参数固体材料的热传导过程。与经典的线性热方程不同，此处扩散项 $\nabla \cdot [a(\nabla y)\nabla y]$ 引入了强烈的非线性耦合，使得标准的线性控制理论无法直接适用。

在控制理论领域，线性及半线性抛物方程的零控性已得到充分发展（如 Fursikov-Imanuvilov 方法）。然而，对于扩散系数显式依赖于梯度 $\nabla y$ 的拟线性系统，尤其是当需要证明局部零控性时，现有的理论框架面临巨大挑战。主要难点在于：线性化后的系统虽然形式上是线性的，但其系数依赖于原非线性项的线性化结果，且为了应用不动点定理或反函数定理，必须建立关于线性化系统解及其控制的高阶正则性估计。特别是，由于非线性项 $a(\nabla y)$ 的存在，解的正则性传播受到限制，这要求初始数据 $y_0$ 具有极高的光滑度（$H^5$），且控制 $v$ 需满足特定的加权能量衰减条件。

本文旨在填补这一空白：首先针对线性化系统建立带权重的 Carleman 估计和高阶正则性先验估计；其次，通过精心构造合适的 Hilbert 空间 $Y$ 和 $Z$，将零控性问题转化为算子方程 $H(y, v) = (0, y_0)$，并证明映射 $H$ 在原点处的严格可微性及导数的满射性；最后，利用局部逆定理证明在初始数据足够小的前提下，系统具有局部零控性。此外，文章还探讨了边界控制的推广以及数值算法的实现。

2. 主要定理或主要结果

本文的核心成果建立在严格的假设 H1 和 H2 之上，证明了拟线性系统在特定正则性条件下的局部零控性。具体结论如下：

  1. 假设条件：

    • 扩散系数 $a: \mathbb{R}^d \to \mathbb{R}$ 属于 $C^4(\mathbb{R}^d)$。
    • 存在常数 $K_0, a_0 > 0$ 和 $r \ge 1$，使得对任意 $q \in \mathbb{R}^d$ 及多重指标 $|\beta| \le 4$，满足椭圆性条件 $a(q) \ge a_0 > 0$ 和增长条件 $|D^\beta a(q)| \le K_0 (1 + |q|^{(r-|\beta|)+})$。
    • 空间维度限制为 $1 \le d \le 3$。

  2. 线性化系统的高阶估计（定理 2.3 及其推论）： 对于线性化系统 $y_t - \sigma \Delta y = \chi_\omega v + f$，作者建立了带指数权重 $\rho_k(t,x) = e^{s\alpha(x)/\ell(t)^{k/2}}$ 的能量不等式。若初始数据 $y_0 \in H^5(\Omega) \cap H^1_0(\Omega)$ 且满足相容性条件 $\Delta y_0, \Delta^2 y_0 \in H^1_0(\Omega)$，则存在控制 $v$ 使得状态在时刻 $T$ 归零，且满足以下高阶正则性估计：

    $$ \int_Q \rho^2 |y|^2 + \int_{]0,T[\times\omega} \rho_3^2 |v|^2 \le C \left( \|y_0\|^2 + \int_Q \rho_3^2 |f|^2 \right). $$

    特别地，对于 $H^5$ 正则性的初值，控制 $v$ 满足：

    $$ \rho_{15} v \in L^2(0, T; H^4(\omega)), \quad \rho_{15} v_t \in L^2(0, T; H^2(\omega)), \quad \rho_{15} v_{tt} \in L^2(Q), $$

    且相应的状态 $y$ 满足 $\rho_{13} y \in L^\infty(0, T; H^3(\Omega))$。

  3. 局部零控性定理（定理 1.3）： 在希尔伯特空间 $Y$（包含高正则性状态和控制）和 $Z$（包含源项和初值）的定义下，非线性映射 $H(y, v) = (y_t - \nabla \cdot [a(\nabla y)\nabla y] - \chi_\omega v, y|_{t=0})$ 是严格可微的，且其导数 $H'(0,0)$ 是从 $Y$ 到 $Z$ 的有界满射算子。由此得出： 对于任意 $T > 0$，系统 (1.1) 是局部零可控的。 即存在 $\eta > 0$，当 $\|y_0\|_{H^5} \le \eta$ 且满足相容性条件时，存在控制 $v \in L^2(]0,T[\times\omega)$ 使得 $y(T, \cdot) = 0$。

  4. 边界控制推论（推论 1.4）： 通过域扩张技巧，证明了同样的局部零控性结论适用于边界控制情形：存在 $\eta > 0$，当 $\|y_0\|_{H^5} \le \eta$ 时，存在边界控制 $h \in L^2(\Sigma)$ 使得解归零，且支持集可限制在边界子集 $\gamma \subset \partial \Omega$ 上。

  5. 数值算法收敛性： 提出了基于准牛顿法的迭代算法 ALG 1：$(y^{n+1}, v^{n+1}) = (y^n, v^n) - H'(0,0)^{-1}(H(y^n, v^n) - (0, y_0))$。文中指出该算法在空间 $Y$ 中收敛，并给出了数值实验验证其在 $d=1, 2$ 维下的有效性。

3. 研究方法、关键技术和核心工具

本文证明的核心障碍在于拟线性项 $\nabla \cdot [a(\nabla y)\nabla y]$ 导致的标准线性化理论失效，特别是难以直接获得控制 $v$ 和状态 $y$ 的一致高阶正则性估计，从而无法满足应用局部逆定理所需的拓扑结构。

  1. 带权重的 Carleman 估计与正则性提升（Bootstrapping 策略）： 作者没有停留在基础的 $L^2$ 能控性估计上，而是深入研究了线性化系统 $y_t - \sigma \Delta y = \chi_\omega v + f$ 的高阶性质。关键工具是 Fursikov-Imanuvilov 的 Carleman 不等式，但作者对其进行了精细的加权处理。通过引入指数权重 $\rho_k(t,x) = e^{s\alpha(x)/\ell(t)^{k/2}}$，其中 $\alpha(x)$ 是与观测区域 $\omega_0$ 相关的几何函数，$\ell(t)$ 是时间截断函数，作者建立了形如 $\int_Q \rho_k^2 |D^m y|^2$ 的高阶能量估计。 这种加权技术的作用在于：在 $t \to T$ 时，权重 $\rho_k$ 趋于无穷大，迫使解及其导数在终端时刻迅速衰减至零，从而保证 $y(T)=0$ 的同时，赋予解在时间区间内部极高的正则性。作者通过归纳法（boot-strapping argument），从低阶估计出发，逐步提高对初值 $y_0$ 的正则性要求（最终达到 $H^5$），换取控制 $v$ 在 $H^4$ 甚至更高阶 Sobolev 空间中的正则性。例如，定理 2.3 和推论 2.4 详细列出了从 $y_0 \in H^1$ 到 $y_0 \in H^5$ 过程中，控制 $v$ 及其时间导数所满足的 $\rho_k$ 加权范数界限。这是证明后续映射可微性的基石。

  2. 泛函分析框架的重构与严格可微性证明： 为了应用逆定理，作者并未直接使用标准的 $L^2$ 空间，而是构造了高度定制化的 Hilbert 空间 $Y$ 和 $Z$。

    • 空间 $Y$ 定义为 $(y, v)$ 的乘积空间，其中 $v$ 不仅要求 $L^2$，还要求其加权的时空导数（如 $\rho_{15} v_{tt}$）属于 $L^2$，且 $y$ 需满足初值的 $H^5$ 正则性和相容性条件。
    • 空间 $Z$ 同样由加权空间和初值空间构成。 在这一框架下，作者定义了非线性映射 $H(y, v) = (y_t - \nabla \cdot [a(\nabla y)\nabla y] - \chi_\omega v, y|_{t=0})$。证明的关键步骤是验证 $H$ 在 $(0,0)$ 处的严格可微性（Strict Differentiability）。这需要证明余项 $R(y,v) = H(y,v) - H(0,0) - H'(0,0)(y,v)$ 满足 $\|R(y,v)\|_Z = o(\|(y,v)\|_Y)$。由于 $a(\nabla y)$ 的非线性，这一估计高度依赖于前述的高阶加权能量估计。作者通过精确计算 $a(\nabla y)\nabla y$ 在 $H^5$ 空间中的展开，结合 Sobolev 嵌入定理（仅在 $d \le 3$ 时成立），证明了 $H'(0,0)$ 不仅是连续的，而且是有界的满射算子，并建立了下界估计 $\|H'(0,0)(y,v)\|_Z \ge M \|(y,v)\|_Y$。

  3. 局部逆定理的应用： 基于上述准备，作者引用了文献 [2] 中的抽象局部逆定理（Theorem 3.2）。该定理保证了在 $H'(0,0)$ 满射的前提下，存在邻域内的右逆映射 $W$，使得 $H(W(z)) = z$。通过将目标状态设为 $(0, y_0)$，直接导出了局部零控性的存在性。这一路径避免了全局不动点论证中可能遇到的无界性或解的破裂问题，但也带来了“局部”性质的局限。

  4. 数值算法设计： 在数值部分，作者利用了线性化算子的逆 $H'(0,0)^{-1}$ 作为预条件器，构建了准牛顿迭代格式。这种方法本质上是在每一步迭代中求解一个带有相同线性主部（$y_t - a(0)\Delta y$）的线性控制问题，利用前面建立的 Carleman 估计对应的最小范数控制求解器来更新步长。

4. 与之前工作的比较

本文与既有偏微分方程控制理论文献相比，在以下几个方面取得了显著推进，同时也揭示了新的局限性：

  1. 扩散系数依赖性的突破： 早期关于抛物方程控制的工作主要集中在常系数或仅依赖于状态 $y$ 而非梯度 $\nabla y$ 的半线性情形（如文献 [19, 27, 28]）。对于扩散系数依赖于梯度的拟线性系统，尽管在存在性理论（De Giorgi-Moser-Nash）上已有深厚积累，但在控制理论方面进展缓慢。文献 [15] 和本文被作者明确指出是首批处理此类问题的论文之一，特别是允许扩散系数为一般幂律形式 $a(p) = a_0 + a_1|p|^r$。此前工作 [13, 21] 处理的是非局部依赖，而 [22, 26, 38] 涉及拟线性方程但可能未完全解决梯度依赖带来的高阶正则性难题。本文明确解决了 $a(\nabla y)$ 这一最困难的情形。

  2. 正则性要求的代价与权衡： 本文的一个核心特征是引入了极高的正则性要求。为了证明局部零控性，初值 $y_0$ 必须属于 $H^5(\Omega)$ 且满足相容性条件，控制 $v$ 也需要在加权空间中具备 $H^4$ 以上的正则性。相比之下，许多线性或半线性系统仅需 $L^2$ 初值即可实现零控性。这种高昂的正则性代价是作者为了克服 $a(\nabla y)$ 导致的非线性项在 $L^2$ 框架下不可微、甚至不连续而不得不付出的。作者承认，这种设定限制了结果的适用范围，使其仅适用于非常光滑的初始状态，但这在物理上对应于平滑的温度场，具有一定的合理性。

  3. 方法论的深化：从能控性到逆定理： 与文献 [22] 类似，本文采用了“线性化 + 逆定理”的策略，但极大地深化了线性化部分的分析。之前的工作可能止步于证明线性化系统的 $L^2$ 能控性，而本文通过引入复杂的加权 Carleman 估计，建立了控制 $v$ 的高阶正则性（$\rho_{15} v \in H^4$ 等）。这一技术细节是连接线性化系统与拟线性系统的关键桥梁，使得映射 $H$ 在 $H^5$ 空间中的严格可微性成为可能。这是本文相对于前人在技术深度上的主要贡献。

  4. 全局零控性的缺失与开放问题： 作者明确指出，目前的方法无法推广到全局零控性。尝试使用全局逆定理（如文献 [35]）需要证明 $H'(\bar{y}, \bar{v})$ 对所有解一致可逆且有界，这归结为求解形如 $y_t - \nabla \cdot [a(\nabla \bar{y})\nabla y + (a'(\nabla \bar{y}) \cdot \nabla y)\nabla \bar{y}] = \chi_\omega v$ 的线性化方程并获得与控制无关的一致估计。作者认为这是一个“开放且可能困难”的问题。这表明本文的结果虽然在理论上完善了拟线性梯度依赖系统的局部控制理论，但在处理大振幅或非光滑初值方面仍存在明显局限，与线性系统中常见的全局能控性形成鲜明对比。

  5. 数值实现的互补性： 不同于纯理论推导，本文提供了具体的准牛顿迭代算法及数值实验。这在控制理论的数值计算文献中（如 [6-8, 20]）较为常见，但将其应用于梯度依赖的拟线性系统尚属首次。数值实验不仅验证了理论的正确性，还展示了算法在实际计算中的收敛行为，弥补了纯理论分析在直观性上的不足。

4. ρ-凸域上球面动力学Sobolev空间的一个迹定理

A trace theorem for spherical kinetic Sobolev spaces on $ρ$-convex domains

作者: Herbert Egger, Matthias Schlottbom

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.19027v1

摘要翻译: 迹定理是分析动力学方程不可或缺的工具。Cessenat及其合作者已在辐射传输及相关应用的广泛情形下建立了这些定理。针对动力学Fokker-Planck方程和Kolmogorov方程，已建立了一系列迹估计，但通常要求底层的域具有光滑性；参见Niebel与Valentini的最新综述。在本工作中，我们证明了关于ρ-凸域（其边界通常可能非光滑）上球面动力学Sobolev空间的一个新迹估计。类似于Cessenat的工作，我们利用特征线方法获得了加权迹空间中带有显式常数的迹估计。为完备起见，我们还给出了Lipschitz域上相应函数空间的一个密度结果。

Comment: 无

链接: http://arxiv.org/abs/2608.19027v1

1. 研究对象和背景

在偏微分方程与分析理论，特别是输运理论与动力学方程的研究中，边界条件的严格数学处理是构建变分框架和证明解的存在唯一性的基石。本文聚焦于球面动力学索博列夫空间 $H^1_{\text{hyp}}(R; S^{d-1})$ 的迹理论问题，该空间定义为：

$$ H^1_{\text{hyp}}(R; S^{d-1}) = \{ v \in L^2(R; H^1(S^{d-1})) : s \cdot \nabla_r v \in L^2(R; H^{-1}(S^{d-1})) \}. $$

其中 $R \subset \mathbb{R}^d$ 是有界区域，$S^{d-1}$ 是 $d \ge 2$ 维单位球面，$H^1(S^{d-1})$ 为其上的索博列夫空间，$H^{-1}(S^{d-1})$ 是其对偶空间。该空间自然地出现在动能福克 - 普朗克（Kinetic Fokker-Planck）方程的研究中，这类方程作为玻尔兹曼输运的前向峰化近似，在粒子输运、放射治疗剂量计算及 MRI 引导下的放疗等应用中具有核心地位。

对于此类方程，基于 Lions 表示定理的存在性证明通常依赖于解在相空间边界上的迹（trace）存在性。然而，要确立唯一的变分解或进行严格的边界值问题处理，必须确保光滑函数在出流边界上为零时在相应的图范数空间中是稠密的。这一性质直接依赖于迹算子的有界性以及光滑函数的稠密性。尽管 Cessenat 及其合作者在辐射输运领域建立了广泛的迹定理，且 Niebel & Valentini 近期的综述总结了针对动能福克 - 普朗克方程和科尔莫戈罗夫方程的一系列迹估计，但这些结果通常要求底层的空间区域 $R$ 具有光滑边界（如 $C^{1,1}$ 正则性）。例如，对于形如

$$ \int_{\partial R \times S^{d-1}} |v(r, s)|^2 |s \cdot n(r)|^p d(r, s) \le C \|v\|^2_{\text{hyp}}, \quad p \ge 1, $$

的估计，已知在 $C^{1,1/(1+p)}$ 正则边界上成立，但对于非光滑边界，反例表明此类估计一般失效。

本文试图填补的空白在于：在区域 $R$ 仅具有“$\rho$-凸性”（$\rho$-convexity），即可能具有非光滑边界（如多面体或有限个球的交集）的情况下，建立适用于球面动力学索博列夫空间的精确迹估计。$\rho$-凸域定义为半径为 $\rho > 0$ 的若干开球的交集：

$$ R = \bigcap_{a \in A} B(a, \rho). $$

此类区域虽然可能是非光滑的，但保持了足够的几何结构（如局部星形性），使得特征线方法中的关键量——逃逸时间（exit time）$\tau(r, s)$ 具有良好的正则性。本文旨在证明在此类广义凸域上，加权迹估计依然成立，并补充证明光滑函数在对应 Lipschitz 域上的稠密性，从而为动力学方程在非光滑几何上的变分分析提供完整的理论基础。

2. 主要定理或主要结果

本文的核心成果包含两个主要定理，分别关于迹估计和稠密性。

定理 1 (球面动力学索博列夫空间的迹定理)：设 $R \subset \mathbb{R}^d$ 为一个 $\rho$-凸域。则对于所有 $v \in C^1(R \times S^{d-1})$，以下不等式成立：

$$ \int_{\partial R \times S^{d-1}} |v(r, s)|^2 \tau(r, s) |s \cdot n(r)| \, d(r, s) \le 2(1 + 4\rho) \|v\|^2_{\text{hyp}}, $$

其中 $\tau(r, s) = \tau_+(r, s) + \tau_-(r, s)$ 是穿过点 $r$ 沿方向 $s$ 的线段在 $R$ 内的长度，$n(r)$ 是边界上的单位外法向量（几乎处处存在）。这里的常数 $2(1+4\rho)$ 是显式的，依赖于区域的曲率半径参数 $\rho$。

通过连续性和密度延拓，迹映射 $\gamma: C^1(R \times S^{d-1}) \to L^2(\partial R \times S^{d-1}; \tau|s \cdot n|)$ 可以延拓为 $H^1_{\text{hyp}}(R; S^{d-1})$ 到加权迹空间的线性连续算子，且上述不等式保持有效。

定理 2 (Lipschitz 域上的稠密性)：设 $R \subset \mathbb{R}^d$ 是一个有界的 Lipschitz 域。则光滑函数空间 $C^\infty(R \times S^{d-1})$ 在 $H^1_{\text{hyp}}(R; S^{d-1})$ 中是稠密的。

由于 $\rho$-凸域是有界且 Lipschitz 的，定理 2 覆盖了定理 1 所需的稠密性条件，从而保证了定理 1 中延拓的合法性。此外，文中还给出了一个关键的比较结果（Remark 2.3）：对于 $\rho$-凸域，逃逸时间与法向量的关系满足 $\tau(r, s)|s \cdot n(r)| \le 2\rho |s \cdot n(r)|^p$ 对所有 $1 \le p \le 2$ 成立。这意味着本文建立的加权迹范数比文献 [26] 中要求的 $|s \cdot n|^p$ ($p \ge 1$) 范数更弱（当 $p \le 2$ 时），从而放宽了对边界的正则性要求。

3. 研究方法、关键技术和核心工具

本文证明的核心技术障碍在于如何在边界不光滑（缺乏 $C^1$ 或 $C^{1,1}$ 正则性）的情况下，控制沿特征线的积分行为。传统的迹估计往往依赖于边界法向量的连续性或曲率的有界性，而在 $\rho$-凸域中，这些性质可能仅在测度意义下存在或不一致。作者的创新点在于利用 $\rho$-凸域的几何结构，证明了逃逸时间 $\tau_\pm(r, s)$ 关于速度方向 $s$ 的 Lipschitz 连续性，并以此为基础重构了特征线方法。

具体而言，证明框架包含以下关键技术步骤：

  1. 逃逸时间的正则性分析 (Lemma 2.1)：这是全文的基石。作者首先指出，对于 $\rho$-凸域，逃逸时间 $\tau_\pm(r, s)$ 可以表示为一系列球体逃逸时间的下确界。利用球体逃逸时间的显式公式 $\tau_{+,a}(r, s) = \sqrt{\rho^2 - |r-a|^2 + \langle r-a, s \rangle^2} - \langle r-a, s \rangle$，通过形式微分和 Rademacher 定理，证明了 $\tau_\pm$ 关于 $s$ 的 Lipschitz 估计：

    $$ |\tau_\pm(r, s) - \tau_\pm(r, s')| \le 2\rho |s - s'|. $$

    这一估计表明，即使边界本身不光滑，只要区域由半径为 $\rho$ 的球交集构成，沿不同方向的特征线长度变化是受控的。这允许作者在积分变换中使用 $s$ 方向的导数估计。

  2. 带权重的特征线恒等式 (Lemma 2.2)：基于上述 Lipschitz 性质，作者构造了一个辅助函数 $\phi(t) = v(r+ts, s)^2 \tau_+(r+ts, s)$。通过对 $\phi(t)$ 求导并利用 $\frac{d}{dt}\tau_+ = -1$ 的性质，结合微积分基本定理，将边界上的积分转化为体积积分：

    $$ \int_{\Gamma_{\text{in}}} |v|^2 |s \cdot n| \tau_+ d(r,s) = \int_R (v^2 - 2v(s \cdot \nabla_r v)\tau_+) dr ds. $$

    这里的关键在于利用 $\tau_+$ 的有界性（$\tau_+ \le 2\rho$）以及 Hölder 不等式和 Young 不等式来处理交叉项 $v(s \cdot \nabla_r v)\tau_+$。最终得到：

    $$ \int_{\partial R \times S^{d-1}} |v|^2 |s \cdot n| \tau_+ \le (1 + 4\rho) \|v\|^2_{\text{hyp}}. $$

    这一过程避免了直接使用边界法向量的高阶导数，而是通过逃逸时间的几何性质来吸收正则性损失。

  3. 稠密性证明的几何缩放策略 (Proposition 3.1)：为了证明光滑函数的稠密性，作者没有直接使用标准的磨光卷积，因为标准磨光可能会破坏定义域。相反，作者引入了全局缩放映射 $T_\lambda(r) = a_s + \lambda(r - a_s)$，其中 $B(a_s, \rho_s)$ 是域内包含的一个球。利用 $\rho$-凸域（或更一般的关于球星形域）的性质，证明了缩放后的像 $T_\lambda(R)$ 严格包含在原域内，且与边界的距离至少为 $(1-\lambda)\rho_s$。这使得可以在 $T_\lambda(R)$ 上安全地进行欧几里得磨光和球面热半群平滑，最后通过取对角序列 $\nu_{\lambda_n, \epsilon_n, \delta_n}$ 逼近原函数。这种方法巧妙地利用了域的星形几何结构，将高维动力学的稠密性问题分解为空间缩放、空间磨光和速度磨光三个独立且可控的步骤。

整个证明框架不依赖边界的光滑性，而是完全建立在 $\rho$-凸性带来的逃逸时间 Lipschitz 性质之上，这是一种从几何结构出发解决分析问题的典型范例。

4. 与之前工作的比较

本文与既有文献在以下几个维度进行了明确的区分与推进：

  1. 边界正则性假设的显著放宽： 之前的工作，如 Niebel & Valentini 的综述 [26] 及相关文献，通常要求区域 $R$ 具有 $C^{1,1}$ 或更高阶的正则性才能建立形如 $\int |v|^2 |s \cdot n|^p d\sigma \le C \|v\|^2$ 的迹估计（见公式 (5)）。这种假设排除了许多在实际应用（如医学物理中的复杂几何模型）中常见的多面体或非光滑区域。本文证明，对于 $\rho$-凸域（包括 Reuleaux 三角形等具有尖角但仍满足特定凸性条件的区域），迹估计依然成立，尽管权函数从单纯的 $|s \cdot n|^p$ 变为 $\tau(r, s)|s \cdot n|$。这是一个实质性的推广，因为它不再要求边界切平面的连续旋转，仅需区域可由固定半径的球交集生成。

  2. 权函数的选择与等价性分析： 本文明确指出了权函数 $\tau(r, s)|s \cdot n|$ 与经典权函数 $|s \cdot n|^p$ 之间的关系。在 Remark 2.3 中，作者推导出不等式 $\tau(r, s)|s \cdot n| \le 2\rho |s \cdot n|^p$ 对于 $1 \le p \le 2$ 成立。这意味着本文建立的迹空间实际上比传统文献中考虑的 $L^2(|s \cdot n|^p)$ 空间更“大”（或者说条件更宽松），因为 $\tau$ 在靠近切向入射点（$s \cdot n \to 0$）时会趋于零，其衰减速度与 $|s \cdot n|$ 相当甚至更快。相比之下，Cessenat 等人的工作 [14] 虽然也使用了类似特征线的方法，但主要针对辐射输运方程，且未深入探讨此类非光滑几何下的显式常数与权函数优化。本文的结果表明，在 $\rho$-凸域上，无法简单地将权函数替换为 $|s \cdot n|^p$ ($p>2$)，这暗示了非光滑边界上迹理论的内在复杂性。

  3. 稠密性结果的完备性： 尽管文献 [2] 和 [12] 曾讨论过光滑函数在动力学索博列夫空间中的稠密性，但它们大多局限于 $C^1$ 或 $C^{1,1}$ 域。本文通过 Proposition 3.1 和 Theorem 2，明确证明了在有界 Lipschitz 域（涵盖所有 $\rho$-凸域）上，$C^\infty(R \times S^{d-1})$ 的稠密性。这一结果是连接迹定理与变分 formulations 的关键缺失环节。之前的研究往往隐含地假设了这一性质，或者仅在光滑域上给出证明。本文提供了一个自包含的、基于几何缩放的证明，填补了非光滑域上变分理论的基础空白。

  4. 代价与局限性： 本文的推进并非没有代价。首先，显式常数 $2(1+4\rho)$ 依赖于 $\rho$，当区域接近退化（即 $\rho \to 0$）时，常数会增大，这可能影响数值稳定性或误差估计的精度。其次，权函数 $\tau(r, s)$ 的引入使得迹空间的结构变得更为复杂，不同于经典的 $L^2$ 边界空间。作者也明确指出（Section 4），对于一般的非均匀凸域或多面体域（不完全满足 $\rho$-凸性），逃逸时间的 Lipschitz 性质可能失效，此时需要引入其他权函数 $w(r, s)$ 来建立估计，这将是未来工作的方向。因此，本文的贡献主要集中在 $\rho$-凸这一特定几何类，而非全类的 Lipschitz 域，尽管稠密性部分已扩展到 Lipschitz 域。

5. 一致局部Sobolev空间中拟线性模式形成系统的Ginzburg-Landau逼近定理

A Ginzburg-Landau approximation theorem for quasilinear pattern-forming systems in uniformly local Sobolev spaces

作者: Théo Belin, Guido Schneider

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.19035v1

摘要翻译: 我们证明了Ginzburg-Landau方程能够正确预测在第一不稳定性附近的拟线性或完全非线性模式形成系统的动力学。我们在一致局部Sobolev空间中建立了一种逼近理论，该理论被应用于一个拟线性Swift-Hohenberg模型、一个具有速度相关粘度的Bénard问题的拟线性版本，以及抽象的拟线性反应-扩散-平流系统。

Comment: 29 pages

链接: http://arxiv.org/abs/2608.19035v1

1. 研究对象和背景

在非线性偏微分方程与分析理论中，描述系统在第一不稳定性附近的模式形成（Pattern Formation）行为是一个核心问题。当系统参数接近临界值时，平凡解失稳，涌现出具有特定空间频率的结构化模式。对于半线性（semilinear）反应扩散系统，Ginzburg-Landau (GL) 方程作为通用的振幅方程已被广泛证明能正确预测原系统的动力学行为，其理论基础建立在多尺度摄动分析之上。

本文的研究对象是**拟线性（quasilinear）**甚至完全非线性的模式形成系统。这类系统通常出现在流体力学和物理应用中，其扩散项或非线性项依赖于未知函数的高阶导数，导致线性化算子的系数随解的变化而变化，从而极大地增加了正则性分析和误差估计的难度。文中具体考虑了以下三类模型：

  1. 拟线性 Swift-Hohenberg 模型：

    $$ \partial_t u = -(1 + \partial_x^2)^2 u + \alpha_{SH} u + g(u, \dots, \partial_x^4 u), $$

    其中 $g$ 是不含常数和线性项的光滑函数，例如 $g = -\partial_x^4(u^3)$。该模型在半线性情形下 ($g=0$ 或低阶非线性) 是研究 Turing 不稳定性和模式形成的标准范例。

  2. 具有速度依赖粘度的拟线性 B´enard 问题： 这是经典热对流问题的推广，其中粘度不再是常数，而是依赖于流速的函数，导致方程呈现拟线性结构。

  3. 抽象拟线性反应 - 扩散 - 输运系统： 涵盖更广泛的拟线性抛物型方程组。

问题设定与空白： 在这些系统中， bifurcating solutions（分岔解）通常在无穷远处不衰减，因此必须在一致局部索博列夫空间（uniformly local Sobolev spaces，记为 $H^r_{l,u}$）中进行讨论。该空间的范数定义为：

$$ \|u\|_{H^r_{l,u}} = \sup_{y\in\mathbb{R}} \left( \sum_{j=0}^r \int_y^{y+1} |\partial_x^j u(x)|^2 dx \right)^{1/2}. $$

尽管对于半线性问题，基于此类空间的 GL 逼近理论已经相当成熟（如 [SU17] 所述），但对于拟线性系统，此前仅有极少数结果（如 [Zim14] 针对 Marangoni 问题，[BS26] 针对特定反应扩散系统）。主要障碍在于拟线性项破坏了标准的谱分解和能量估计的简单结构，使得处理长时标 $O(\epsilon^{-2})$ 下的误差累积变得极其困难。本文旨在填补这一空白，建立适用于拟线性系统的通用 GL 逼近理论，并给出严格的误差界。

2. 主要定理或主要结果

本文的核心成果是在一致局部索博列夫空间中建立了拟线性模式形成系统的 Ginzburg-Landau 逼近定理。作者证明了在原系统参数处于弱不稳定性区域（即 $\alpha_{SH} = \epsilon^2$）时，原系统的解可以由 Ginzburg-Landau 方程的解通过特定的缩放和修正来近似，且误差具有明确的时间尺度和阶数。

主要定理（以拟线性 Swift-Hohenberg 模型为例，见 Theorem 2.1）： 设 $A(X, T) \in C([0, T_0], H^{r_0+2}_{l,u})$ 是实 Ginzburg-Landau 方程

$$ \partial_T A = 4\partial_X^2 A + A - 3A|A|^2 $$

的一个解，且满足一致有界性条件：

$$ \sup_{T\in[0,T_0]} \|A(\cdot, T)\|_{H^{r_0+2}_{l,u}} \leq C_{GL}. $$

则存在 $\epsilon_0 > 0$ 和 $C_2 > 0$（仅依赖于 $C_{GL}, T_0, r_0 \geq 5$），使得对于所有 $\epsilon \in (0, \epsilon_0)$，拟线性 Swift-Hohenberg 模型 (1) 存在解 $u(x,t)$，满足如下误差估计：

$$ \sup_{t\in[0, T_0/\epsilon^2]} \| u(\cdot, t) - \epsilon \Psi_{GL}(\cdot, t) \|_{H^{r_0}_{l,u}} \leq C_2 \epsilon^{3/2}. $$

其中 $\Psi_{GL}$ 是通过多重尺度展开得到的近似解：

$$ \Psi_{GL}(x,t) = A(\epsilon x, \epsilon^2 t) e^{ix} + c.c. $$

技术细节补充：

  1. 时间尺度：逼近成立的时间区间为 $t \in [0, T_0/\epsilon^2]$，这对应于 GL 方程的时间变量 $T \in [0, T_0]$，覆盖了模式形成的早期演化阶段。
  2. 误差阶：在 $H^{r_0}_{l,u}$ 范数下，误差为 $O(\epsilon^{3/2})$。通过索博列夫嵌入定理 $H^{r_0}_{l,u} \hookrightarrow C^0_{b,unif}$（当 $r_0 \geq 1$），该结果直接蕴含了点态估计：
    $$ \sup_{t\in[0, T_0/\epsilon^2]} \sup_{x\in\mathbb{R}} |u(x,t) - \epsilon \Psi_{GL}(x,t)| \leq C_2 \epsilon^{3/2}. $$
  3. 改进近似：为了达到 $O(\epsilon^{7/2})$ 的残差精度以支撑 $O(\epsilon^{3/2})$ 的误差界，作者引入了高阶修正项，构造了改进近似 $\epsilon \Psi = \epsilon \Psi_{GL} + \epsilon^3 A_3 e^{3ix} + c.c.$，其中 $A_3$ 由代数方程确定。
  4. 正则性要求：定理要求初始振幅 $A$ 具有 $H^{r_0+2}_{l,u}$ 的正则性，且 $r_0 \geq 5$，以确保在处理拟非线性项时所需的乘积规则和复合函数估计有效。

3. 研究方法、关键技术和核心工具

本文的技术核心在于解决拟线性项带来的**最优正则性估计（optimal regularity estimates）**问题，以及如何在长时标下控制误差增长。传统的半线性方法在此失效，因为拟线性项会导致导数损失或需要更高阶的时空正则性。作者构建了一套严密的分析框架，具体包含以下关键步骤和创新点：

  1. 时空加权空间与 Hölder 连续性： 为了处理长时标 $O(\epsilon^{-2})$ 上的解，作者引入了带有指数衰减权重的时空 Hölder 空间。 定义空间 $X^{r,\alpha}_{\eta, t_0} = C^\alpha([0, t_0], X)$，其中 $X = H^r_{l,u}$，并赋予范数：

    $$ \|R\|_{X^{r,\alpha}_{\eta, t_0}} = \sup_{t\in[0,t_0]} \|e^{-\eta t} R(t)\|_{C^\alpha([0,t_0], X)}. $$

    同时定义解空间 $Y^{r+4,\alpha}_{\eta, t_0} = C^\alpha([0, t_0], D(\Lambda)) \cap C^{1+\alpha}([0, t_0], X)$，其中 $D(\Lambda) = H^{r+4}_{l,u}$。 这些空间的选择至关重要：它们不仅捕捉了解在时间上的 Hölder 连续性（这对于处理拟非线性项的复合结构是必须的），而且通过权重 $e^{-\eta t}$ 允许在长时间积分中控制解的增长。

  2. 模态滤波与谱分解（Mode Filtering & Spectral Decomposition）： 为了分离临界模态（critical modes）和稳定模态（stable modes），作者引入了傅里叶空间中的模态滤波器 $E_c, E_\sigma$。 将误差方程 $\partial_t R = \Lambda R + f$ 分解为临界部分 $R_c$ 和稳定部分 $R_\sigma$：

    $$ \partial_t R_c = \Lambda R_c + E_c f, \quad \partial_t R_\sigma = \Lambda R_\sigma + E_\sigma f. $$

    对于稳定部分 $R_\sigma$，利用解析半群的指数衰减性质（exponential damping），结合最优正则性引理进行估计。 对于临界部分 $R_c$，由于其特征值位于虚轴附近，不能直接利用衰减，必须利用紧支集乘子（compact support multipliers）的性质，将其转化为半线性系统的标准形式进行处理。

  3. 最优正则性引理（Optimal Regularity Lemma）： 这是证明中最具技术含量的部分（Lemma 2.7）。作者引用并适配了 Lunardi (1995) 关于抛物方程最优正则性的理论，但在一致局部索博列夫空间和拟线性背景下进行了推广。 引理表明，对于方程 $\partial_t R = \Lambda R + f$，若 $f \in Z_0$（即 $C^\alpha$ 连续），则解 $R$ 属于 $Z_1$（即 $C^\alpha$ 连续到 $D(\Lambda)$，且 $C^{1+\alpha}$ 连续到 $X$），并给出了具体的范数估计：

    $$ \|R\|_{Z_1} \leq C ((\bar{\eta}-h)^{-1}\epsilon^{-2} + 1) (\|f\|_{Z_0} + \|R_0\|_{D(\Lambda)} + \|\Lambda R_0 + f|_{t=0}\|_{D_\Lambda(\alpha,\infty)}). $$

    这里的 $D_\Lambda(\alpha,\infty)$ 是半群生成的插值空间，用于处理初值与源项在 $t=0$ 时刻的不兼容性。这一估计确保了即使存在拟非线性项导致的复杂结构，解的正则性依然能够保持在可控范围内。

  4. 不动点论证（Fixed Point Argument）： 最终，作者将误差方程写为 $R = K(f(R))$ 的形式，其中 $K$ 是线性算子的逆映射。通过在空间 $Z_1$ 的球内构造压缩映射，证明了存在唯一的固定点。关键在于选择合适的参数 $\beta = 3/2$ 和足够大的权重 $\eta$，使得非线性项 $f(R)$ 中的高阶项（如 $\|R\|^2, \|R\|^3$）在长时间尺度下被 $\epsilon$ 的小量所压制，从而保证映射的收缩性。

  5. 残差估计的精细处理： 为了获得 $O(\epsilon^{3/2})$ 的误差界，必须将残差 $Res(\epsilon \Psi)$ 控制在 $O(\epsilon^{7/2})$。这要求不仅空间正则性要高，时间正则性也要达到 $C^\alpha$ 级别（Lemma 2.5），这是为了适应上述加权空间的要求，也是区别于以往仅关注 $L^2$ 或 $H^r$ 估计的关键创新。

4. 与之前工作的比较

本文与既有文献相比，在以下几个方面实现了显著的推进和突破：

  1. 从半线性到拟线性的跨越： 现有的 GL 逼近理论主要集中在半线性系统（如经典的 Swift-Hohenberg 方程、Gray-Scott 系统等）。虽然 [SU17] 教科书详细综述了半线性情况，但拟线性系统（Quasilinear systems）由于其扩散系数或非线性项依赖于高阶导数，导致线性化算子不再具有简单的谱性质，且正则性传播更为复杂。此前已知的拟线性结果非常有限：

    • [Zim14] 仅针对 Marangoni 问题在 Sobolev 空间 $H^r$ 中建立了结果，未涉及一致局部空间及更广泛的拟线性结构。
    • [BS26] 提供了一个针对特定拟线性反应 - 扩散 - 输运系统的“易用”逼近结果，但其适用范围较窄，且未深入探讨最优正则性估计的普适框架。 本文首次系统地建立了适用于一般拟线性模式形成系统（包括 Swift-Hohenberg 型、B´enard 型及抽象反应扩散型）的 GL 逼近理论，并统一在一致局部索博列夫空间 $H^r_{l,u}$ 的框架下。

  2. 函数空间的推广： 大多数早期工作假设解在 $L^2$ 或标准 Sobolev 空间 $H^r(\mathbb{R})$ 中衰减至零。然而，模式形成系统中的分岔解通常表现为行波或周期性结构，在无穷远处不消失。本文明确指出并严格处理了这一问题，采用了一致局部索博列夫空间 $H^r_{l,u}$。该空间允许解在局部区域具有任意高的正则性，同时在整体上不要求衰减，完美契合了物理上观察到的无限长条纹或图案的形成机制。这一选择使得理论能够覆盖更广泛的物理场景。

  3. 时间尺度与误差精度的优化： 本文的证明框架成功地将逼近的有效时间尺度推到了 $O(\epsilon^{-2})$，这与 GL 方程的自然时间尺度一致。更重要的是，通过引入时间 Hölder 连续性和最优正则性引理，作者克服了拟线性项在长时间积分中可能引发的正则性损失问题。相比之下，之前的拟线性工作往往受限于较短的时间尺度或较低的误差精度（例如仅 $O(\epsilon)$ 或依赖更强的正则性假设）。 本文获得的误差阶为 $O(\epsilon^{3/2})$，这是在 $H^r_{l,u}$ 范数下对拟线性系统的一个显著进步，且该结果是基于严格的解析推导而非启发式论证。

  4. 技术工具的深化： 与以往工作相比，本文最大的技术贡献在于对最优正则性估计的处理。在拟线性情形下，线性化算子的系数随解变化，传统的能量估计（Energy Estimates）往往不足以控制高阶导数。本文巧妙地结合了模态滤波（分离临界与稳定模态）、带权重的时空 Hölder 空间以及半群理论的插值空间技巧，构建了一个自洽的分析闭环。特别是 Lemma 2.7 的证明，展示了如何处理 $t=0$ 时刻的不兼容项以及如何利用指数权重控制长时行为，这些技术在处理其他类型的拟线性 PDE 时也具有借鉴意义。

  5. 局限性与代价： 当然，这种推广也付出了代价。首先，对初始数据的正则性要求较高（需 $r_0 \geq 5$ 以保证复合函数的可微性）。其次，证明过程高度依赖于线性化算子的解析半群性质，对于某些强耦合或非耗散型的拟线性系统可能需要额外的修改。此外，目前的框架主要针对第一不稳定性附近的弱非线性区域，对于远离临界点的强非线性动力学（如混沌或湍流）尚不适用。尽管如此，本文的工作为理解拟线性流体和材料科学中的模式形成提供了坚实的数学基础。

6. 具有Carathéodory强迫项的非线性耗散演化方程的存在性与定性理论

Existence and qualitative theory for nonlinear accretive evolutions with Carathéodory forcing

作者: Dieter Bothe

类别: math.FA, math.AP

发布日期: 2026-08-18

arXiv ID: 2608.18338v1

摘要翻译: 设$A$是实Banach空间$X$中的$m$-耗散算子。给定闭集$K(t)

Comment: 无

链接: http://arxiv.org/abs/2608.18338v1

1. 研究对象和背景

本文研究实Banach空间$X$中一类受时变约束控制的非线性吸积演化方程的解的存在性及其定性性质。具体而言，设$A: D(A) \subset X \to 2^X$为$m$-吸积算子（$m$-accretive operator），时间区间为$J=[0, a]$，并给定一族闭集$K(t) \subset X$表示随时间演化的约束集合。定义有效管（effective tube）为$K_A(t) := K(t) \cap D(A)$，其图像记为$\text{gr}(K_A)$。研究的动力学系统由以下微分包含式描述：

$$ u'(t) + A u(t) \ni f(t, u(t)), \quad t \in J, $$

其中强迫项$f: \text{gr}(K_A) \to X$在有效管上定义为Carathéodory映射（即对每个$x$，$f(\cdot, x)$强可测；对几乎每个$t$，$f(t, \cdot)$连续）。解$u(t)$必须满足时变约束条件$u(t) \in K_A(t)$对所有$t \in [t_0, a]$成立。

该问题的核心难点在于强迫项$f$的定义域是移动的、非圆柱形的集合$\text{gr}(K_A)$，而非固定的柱体$J \times X$。这导致$f$与解路径$u(\cdot)$的复合函数$f(\cdot, u(\cdot))$的可测性无法直接由Carathéodory性质保证，除非附加联合可测性假设或存在向固定柱体的Carathéodory延拓。此外，约束集合$K(t)$随时间变化，要求解轨迹始终保持在移动的有效管内部，这对应于可行性（viability）理论中的流不变性（flow invariance）问题。

已有理论背景主要集中在两个方面：一是Crandall-Liggett定理生成的非线性半群理论，用于处理无约束或静态约束下的$m$-吸积演化方程；二是基于Nagumo切向量定理的可行性理论，处理常微分方程和微分包含在时变集合上的不变性。对于无限维情形，Pavel、Pierre和Vrabie等人发展了适应半群理论的切性条件。然而，现有文献多集中于强迫项连续且约束静态的情形，或者强迫项连续但约束移动的情形。对于强迫项仅为Carathéodory正则性（允许不连续的时间依赖）且约束严格时变的情况，特别是如何处理定义域限制导致的可测性障碍和逼近过程中的约束保持问题，尚缺乏系统的存在性理论。

本文试图填补的空白是：在强迫项$f$仅具备Carathéodory正则性（甚至仅在有效管上定义）且约束$K(t)$随时间剧烈变化的情况下，建立$m$-吸积演化方程强解（mild solution）的存在性理论。作者不仅解决了可测性难题，还通过引入特定的逼近策略和子空间构造，证明了在满足线性增长条件和次切向（subtangential）条件的情况下，存在满足时变约束的解。这一理论框架进一步推广了比较原理、Lyapunov对理论以及周期解的存在性分析，适用于抽象反应扩散系统等应用模型。

2. 主要定理或主要结果

本文在实Banach空间$X$中，针对$m$-吸积算子$A$和时变闭集族$K(t)$，建立了以下主要结果：

  1. 存在性定理（Theorem on Existence of Viable Solutions）： 假设强迫项$f: \text{gr}(K_A) \to X$满足以下条件：

    • (i) $f$在$\text{gr}(K_A)$上是Carathéodory映射；
    • (ii) $f$在移动图像上联合可测，或者存在向固定柱体$J \times X$的Carathéodory延拓；
    • (iii) 线性增长条件：存在$c \in L^1_+(J)$使得对所有$(t, x) \in \text{gr}(K_A)$，有
      $$ \|f(t, x)\| \le c(t)(1 + \|x\|); $$
    • (iv) 次切向条件（Subtangential Conditions）： 对于几乎所有$t \in J$和所有$x \in K_A(t)$，有$f(t, x) \in T^A_K(t, x)$，其中$T^A_K(t, x)$是$A$-次切锥，定义为
      $$ T^A_K(t, x) := \left\{ z \in X : \liminf_{h \downarrow 0} \frac{d(S_z(h)x, K_A(t+h))}{h} = 0 \right\}, $$
      这里$S_z(h)$是由$-Az$生成的收缩半群。此外，在“例外时间”（exceptional times）上需满足相应的控制条件。

    在此假设下，对于任意初始值$(t_0, x_0) \in \text{gr}(K_A)$，存在至少一个强解（mild solution）$u: [t_0, a] \to D(A)$满足：

    $$ u(t) \in K_A(t), \quad \forall t \in [t_0, a], $$

    且$u$是微分包含$u' + Au \ni f(\cdot, u)$的解。

  2. 局部Lipschitz强迫项下的适定性： 若强迫项$f$关于状态变量局部Lipschitz连续，则在上述条件下解是唯一的，从而得到适定性（well-posedness）结果。

  3. 定性性质的推论： 基于上述存在性理论，文章推导了以下定性结论：

    • 比较原理（Comparison Principles）：若$A$是$m$-$T$-吸积算子（即在Banach格中保持序），且$f$满足拟单调条件，则解保持序关系。
    • 非自守Lyapunov对（Nonautonomous Lyapunov Pairs）：通过构造移动的上包络管（moving epigraphs），导出了非自守Lyapunov不等式，用于刻画耗散性。
    • 周期解（Periodic Solutions）：在满足返回条件（return conditions，即$K(T) \subset K(0)$）和紧性假设下，存在周期解。
    • 时变边界估计（Time-dependent Bounds）：对于抽象反应扩散系统，利用移动矩形管$K_y(t) = \{u \in L^p(\Omega; \mathbb{R}^m) : 0 \le u \le y(t)\}$，证明了正性和有界性。

  4. 关键的技术公式： 证明中使用的核心不等式包括解的稳定性估计：

    $$ \|u(t) - v(t)\| \le \|u(s) - v(s)\| + \int_s^t [u(\tau) - v(\tau), f(\tau) - g(\tau)]_+ d\tau, $$

    以及次切向条件的等价形式：

    $$ \lim_{h \downarrow 0} h^{-1} d(S_{f(t,x)}(h)x, K_A(t+h)) = 0 \quad \text{a.e. } t. $$

3. 研究方法、关键技术和核心工具

本文的证明面临两个核心技术障碍：一是Carathéodory强迫项在非圆柱形定义域上的可测性问题，二是时变约束下逼近序列的可行性（viability）保持问题。作者通过以下创新方法和分析框架克服了这些困难：

  1. 可测性障碍的克服：Scorza-Dragoni性质与子空间约化 由于$f$定义在移动图像$\text{gr}(K_A)$上，直接处理$f(\cdot, u(\cdot))$的可测性极其困难。作者首先利用线性增长条件将问题约化到一个有界的有效管上，并构造了一个特殊的闭可分子空间$Y \subset X$。该子空间的关键性质包括：

    • Resolvent-invariant：对任意$\lambda > 0$，有$J_\lambda Y \subset Y$，即$Y$对$A$的预解式保持不变。
    • 距离保持：$Y$保留了到移动集合$K_A(t)$的距离信息，即$d(x, K_A(t))$在$Y$中与在$X$中一致。
    • Scorza-Dragoni性质：通过将问题限制在$Y$上，并利用$Y$的可分性，作者证明了受限强迫项具有Scorza-Dragoni性质（即除去任意小测度的时间集后，映射在剩余集上连续）。这使得原本仅Carathéodory正则的$f$在约化空间中表现出近乎连续性，从而保证了复合函数$f(\cdot, u(\cdot))$的强可测性。 这一步骤是将无限维、非圆柱域上的不可测问题转化为有限维或可分空间上可处理问题的关键。

  2. 可行性保持的核心：两级逼近策略（Two-level Approximation） 这是本文最核心的技术创新。传统的Euler折线法在处理时变约束时，往往因为离散点落在约束外而导致失效。作者设计了一种包含三个分量的逼近结构：

    • 辅助路径$v_\varepsilon$：由向后Euler格式生成的阶梯函数近似解。
    • 当前时间选择器$x_\varepsilon(t)$：对于几乎每个$t$，选取$x_\varepsilon(t) \in K_A(t)$，并令强迫项取值为$w_\varepsilon(t) = f(t, x_\varepsilon(t))$。这一步确保了强迫项的取值严格位于定义域内，解决了定义域限制问题。
    • 滞后节点$\sigma_\varepsilon(t)$：对于每个$t$，存在$\sigma_\varepsilon(t) \in [(t-\varepsilon)^+, t]$使得$v_\varepsilon(\sigma_\varepsilon(t)) \in K_A(\sigma_\varepsilon(t))$。这个滞后节点机制保证了逼近序列在历史时刻满足约束，从而利用半群的收缩性和次切向条件，确保极限解在每一时刻都满足$u(t) \in K_A(t)$。

    这种“当前选择器识别极限强迫，滞后节点保证可行性”的双重机制，巧妙地分离了可测性问题和几何约束问题。作者通过精细的估计证明了当$\varepsilon \to 0$时，$v_\varepsilon$收敛到$m$-吸积方程的强解，且该解始终位于有效管内。

  3. 分析框架的推进 作者没有直接处理全局无界问题，而是先通过线性增长条件和局部极大值函数（local majorant）将问题限制在有界区域。随后，利用$m$-吸积算子的生成理论（Crandall-Liggett定理）和半群估计，结合上述逼近策略，建立了存在性。整个证明框架从“全局无界”到“局部有界”，再到“可分子空间约化”，最后通过“两级逼近”完成极限过程，逻辑严密且针对性强。

4. 与之前工作的比较

本文与既有文献相比，在以下几个方面实现了显著的推进和差异：

  1. 强迫项正则性的放宽： 既往工作如[17]（Bothe等）主要处理强迫项连续（continuous forcing）的情形，而[21]处理了静态约束下的Carathéodory强迫项。本文首次系统地处理了时变约束下Carathéodory强迫项的存在性问题。这意味着强迫项可以在时间上完全不连续，仅需满足Carathéodory条件（逐点连续、几乎处处可测），这大大扩展了理论的应用范围，使其能够涵盖更多物理和工程中的非光滑源项。

  2. 定义域处理的创新： 前人工作通常假设强迫项定义在固定的柱体$J \times X$上，或者假设约束是静态的。本文明确处理了强迫项仅定义在移动图像$\text{gr}(K_A)$上的情况。作者引入了“联合可测性”或“Carathéodory延拓”作为额外假设，并通过Scorza-Dragoni性质和子空间约化技术，成功解决了定义域限制带来的可测性障碍。这是前人未涉及的理论难点。

  3. 时变约束的处理机制： 在时变约束下，经典的Nagumo切向量条件虽然必要，但不足以保证逼近序列的可行性。本文提出的两级逼近策略（特别是滞后节点$\sigma_\varepsilon(t)$的使用）是对传统Euler离散法的重大改进。相比于仅依靠紧性论证或测度非紧性（measure of noncompactness）的方法，本文的方法在不需要强紧性假设的情况下，通过构造特定的选择器和滞后节点，直接保证了极限解的可行性。这使得理论适用于更广泛的非紧情形。

  4. 适用范围与代价：

    • 推进：本文将可行性理论从静态约束推广到时变约束，从连续强迫推广到Carathéodory强迫，并且无需假设解的唯一性即可证明存在性。同时，文章建立了完整的定性理论框架，包括比较原理、Lyapunov对和周期解，这些结果在之前的文献中多为孤立出现或仅针对特定模型。
    • 代价/局限：为了处理Carathéodory强迫项和时变约束，本文引入了额外的技术假设，如线性增长条件和严格的次切向条件（在例外时间上也需满足控制）。此外，证明过程依赖于构造特殊的子空间和复杂的逼近序列，这在计算实现上可能较为复杂。与[22]相比，本文未显式讨论范围条件（range conditions）的推广，而是专注于存在性本身。如果强迫项不具备局部Lipschitz性质，唯一性将丢失，此时需要借助紧半群或测度非紧性假设来替代，这是本文指出的另一种存在机制，但在主定理中未展开详细讨论。

总体而言，本文在方法论上填补了时变约束下非光滑强迫项演化方程的理论空白，为抽象反应扩散系统及其他无限维动力系统的定性分析提供了强有力的工具。

7. $Γ$-收敛：有界区域上弱型非局部泛函

$Γ$-Convergence of Weak-Type Nonlocal Functionals on Bounded Domains

作者: Xiaosheng Lin, Dachun Yang, Sibei Yang, Wen Yuan, Yangyang Zhang

类别: math.CA, math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18414v1

摘要翻译: 设 $N\ge1$, $p\in[1,\infty)$, $γ\in(0,\infty)$，且当 $N=1$ 时 $Ω\subset\mathbb R^N$ 为有界开区间，当 $N\ge2$ 时为有界 Lipschitz 区域。对于任意 $λ\in(0,\infty)$ 和任意可测函数 $u$，考虑弱型非局部泛函 \begin{align*} G_{λ,p,γ}(u;Ω) :=λ\iint_{Ω\timesΩ} \mathbf 1_{\left{(x,y)\inΩ\timesΩ:\ x\neq y,\ \frac{|u(x)-u(y)|^p}{|x-y|^{p+γ}}\geqλ\right}} |x-y|^{γ-N},dx,dy. \end{align*} 在本文中，我们证明了当 $λ\to\infty$ 时，族 $G_{λ,p,γ}$ 在 $L^p(Ω)$ 中按 $Γ$-收敛意义收敛于泛函 \begin{align*} Ψ_{p,γ}^{\mathrm{cell}}(u;Ω):= \begin{cases} C_{N,p,γ}^{\mathrm{cell}}\displaystyle\int_Ω|\nabla u|^p,dx, &p\in(1,\infty)\ \hbox{and}\ u\in W^{1,p}(Ω),\[2mm] C_{N,1,γ}^{\mathrm{cell}}|Du|(Ω), &p=1\ \hbox{and}\ u\in BV(Ω),\[1mm] \infty,&\hbox{otherwise}, \end{cases} \end{align*} 其中正常数 $C_{N,p,γ}^{\mathrm{cell}}$ 与 $Ω$ 无关，并由单元公式刻画。这给出了 Brezis [Open Problem~9.3, Rend. Lincei Mat. Appl. 2023] 所提问题的肯定回答。

Comment: 42 pages

链接: http://arxiv.org/abs/2608.18414v1

1. 研究对象和背景

本文研究的核心对象是定义在有界域 $\Omega \subset \mathbb{R}^N$ 上的弱型非局部泛函族 $G_{\lambda,p,\gamma}$。具体而言，给定参数 $N \ge 1$, $p \in [1, \infty)$, $\gamma \in (0, \infty)$，对于任意可测函数 $u: \Omega \to \mathbb{R}$ 和尺度参数 $\lambda \in (0, \infty)$，该泛函定义为：

$$ G_{\lambda,p,\gamma}(u;\Omega) := \lambda \iint_{\Omega \times \Omega} \mathbf{1}_{\left\{(x,y)\in\Omega\times\Omega:\ x\neq y,\ \frac{|u(x)-u(y)|^p}{|x-y|^{p+\gamma}}\geq\lambda\right\}} |x-y|^{\gamma-N}\,dx\,dy. $$

这一研究对象源于对 Sobolev 空间范数非局部逼近理论的深入探索。经典的 Bourgain–Brezis–Mironescu (BBM) 公式表明，当 $s \to 1-$ 时，分数阶 Sobolev 能量 $(1-s) \iint \frac{|u(x)-u(y)|^p}{|x-y|^{N+sp}} dx dy$ 收敛于 $C_{N,p} \int |\nabla u|^p dx$。然而，在端点情形或不同参数化下，存在多种非局部能量形式。特别地，Brezis 等人在研究强型（strong-type）与弱型（weak-type）非局部泛函的关系时指出，对于 $F_{\lambda,p,\gamma}$ 类泛函（基于超水平集测度），在全空间 $\mathbb{R}^N$ 上已建立了渐近估计。但 Brezis 在《Open Problem 9.3》中明确提出了一个悬而未决的问题：在有界域 $\Omega$ 上，当 $\gamma = N$ 时，泛函族 $G_{\lambda,p,N}$ 是否在 $L^p(\Omega)$ 拓扑下发生 $\Gamma$-收敛？其极限是否为形如 $c_p \int_\Omega |\nabla u|^p dx$ 的 Sobolev 能量（当 $p>1$）或 $c_1 |Du|(\Omega)$（当 $p=1$）？

这一问题之所以重要且困难，在于边界效应和非局部相互作用的复杂性。在非局部泛函中，点 $x$ 处的能量贡献不仅依赖于 $x$ 附近的梯度信息，还依赖于 $x$ 与域外点的相互作用（若延拓到全空间）或与域内远端点的耦合。在有界域上，这种“长程”相互作用导致泛函值受到边界几何形状的显著影响。此外，当 $p=1$ 时，目标空间从 Sobolev 空间 $W^{1,1}$ 转变为有界变差空间 $BV$，此时函数的正则性结构更为复杂（包含跳跃间断），且点态极限可能不存在（如 Lahti 的反例所示），这使得直接通过点态收敛来推导 $\Gamma$-极限变得不可行，必须依赖更精细的 $\Gamma$-收敛理论框架。本文旨在填补这一空白，证明对于任意 $\gamma \in (0, \infty)$（而不仅是 $\gamma=N$），该泛函族确实在 $L^p(\Omega)$ 中 $\Gamma$-收敛到上述形式的局部泛函，并精确刻画了比例常数。

已有理论背景主要涉及 BBM 公式的推广、Sobolev 空间的非局部刻画以及图像处理和变分法中的相关模型。Nguyen, Gobbino, Picenni 等学者曾在全空间情形或部分参数范围内取得进展，例如 Nguyen 证明了当 $\gamma \in (-1, 0)$ 时极限为零泛函，而 Lahti 和 Li 在特定收敛意义下研究了 $p=1$ 的情形。然而，针对 Brezis 提出的强 $L^p(\Omega)$ 拓扑下的 $\Gamma$-收敛问题，特别是涉及一般 $\gamma > 0$ 且有界域情形的严格证明，此前尚未完成。本文试图解决这一核心难题，并建立完整的 $\Gamma$-收敛理论框架。

2. 主要定理或主要结果

本文的主要成果是证明了弱型非局部泛函族 $G_{\lambda,p,\gamma}$ 在有界 Lipschitz 域上的 $\Gamma$-收敛性，并给出了精确的极限泛函表达式及常数特征。

设 $\Omega \subset \mathbb{R}^N$ 为有界开区间（$N=1$）或有界 Lipschitz 域（$N \ge 2$）。令 $\{\lambda_j\}_{j\in\mathbb{N}} \subset (0, \infty)$ 为满足 $\lambda_j \to \infty$ 的任意序列。定义极限泛函 $\Psi_{p,\gamma}^{\mathrm{cell}}: L^p(\Omega) \to [0, \infty]$ 如下：

$$ \Psi_{p,\gamma}^{\mathrm{cell}}(u;\Omega):= \begin{cases} C_{N,p,\gamma}^{\mathrm{cell}} \displaystyle\int_\Omega |\nabla u|^p\,dx, & p\in(1,\infty)\ \hbox{and}\ u\in W^{1,p}(\Omega),\\ C_{N,1,\gamma}^{\mathrm{cell}}|Du|(\Omega), & p=1\ \hbox{and}\ u\in BV(\Omega),\\ \infty,&\hbox{otherwise}, \end{cases} $$

其中 $C_{N,p,\gamma}^{\mathrm{cell}}$ 是一个仅依赖于维数 $N$、指数 $p$ 和参数 $\gamma$ 的正实数常数，由所谓的“胞腔公式”（cell formula）确定：

$$ C_{N,p,\gamma}^{\mathrm{cell}} := \inf \left\{ \liminf_{j\to\infty} G_{\lambda_j,p,\gamma}(v_j; Q_0) : \lambda_j \to \infty,\ v_j \to \ell \text{ in } L^p(Q_0) \right\}, $$

这里 $Q_0 = (-\frac{1}{2}, \frac{1}{2})^N$ 是单位立方体，$\ell(x) = x_N$ 是线性函数。

定理 1.1 ($\Gamma$-收敛)：在上述设定下，泛函族 $G_{\lambda_j,p,\gamma}(\cdot; \Omega)$ 在 $L^p(\Omega)$ 拓扑下 $\Gamma$-收敛到 $\Psi_{p,\gamma}^{\mathrm{cell}}$。具体包含两个不等式：

  1. $\Gamma$-下极限不等式：对于任意 $u \in L^p(\Omega)$ 及任意满足 $u_j \to u$ 在 $L^p(\Omega)$ 中的序列 $\{u_j\}$，有
    $$ \liminf_{j\to\infty} G_{\lambda_j,p,\gamma}(u_j; \Omega) \ge \Psi_{p,\gamma}^{\mathrm{cell}}(u; \Omega). $$
  2. $\Gamma$-上极限不等式（恢复序列）：对于任意 $u \in L^p(\Omega)$，存在序列 $\{u_j\} \subset L^p(\Omega)$ 使得 $u_j \to u$ 在 $L^p(\Omega)$ 中，且
    $$ \limsup_{j\to\infty} G_{\lambda_j,p,\gamma}(u_j; \Omega) \le \Psi_{p,\gamma}^{\mathrm{cell}}(u; \Omega). $$
    称这样的序列 $\{u_j\}$ 为 $u$ 的恢复序列（recovery sequence）。

此外，文章还确立了常数 $C_{N,p,\gamma}^{\mathrm{cell}}$ 的上下界估计：

$$ \frac{K_{N,p}}{\log 2 \cdot \frac{2^{\gamma+1}-1}{\gamma}} \le C_{N,p,\gamma}^{\mathrm{cell}} \le \frac{K_{N,p}}{\gamma}, $$

其中 $K_{N,p} = \int_{S^{N-1}} |\omega \cdot e_N|^p d\mathcal{H}^{N-1}(\omega)$ 是 BBM 公式中的经典常数。这表明极限常数与 BBM 常数密切相关，但在有界域和弱型设定下具有独特的“胞腔”修正因子。

3. 研究方法、关键技术和核心工具

证明的核心障碍在于处理有界域上的边界效应以及非局部算子在不同参数范围（特别是 $p=1$ 时）下的奇异性。作者并未采用传统的二进分割（dyadic discretization）或切片论证（slicing arguments，如 Gobbino 和 Picenni 所用），而是发展了一套基于“胞腔常数”（cell constant）的直接比较分析和精细的构造方法。

1. 下极限不等式的证明策略（The Gamma-Liminf Inequality）

对于 $p \in (1, \infty)$ 的情况，作者利用 Sobolev 函数的性质。关键在于证明：在几乎每个梯度非零的点处，函数的归一化缩放（normalized blow-ups）收敛到一个仿射函数。具体而言，对于 $u \in W^{1,p}(\Omega)$，在足够小的立方体 $Q$ 上，通过 Morse 类型的覆盖引理（Morse-type covering argument），将全局能量下界转化为局部能量下界。由于在小尺度下 $u$ 近似于仿射函数，其能量可由预先定义的胞腔常数 $C_{N,p,\gamma}^{\mathrm{cell}}$ 控制。即：

$$ G_{\lambda,p,\gamma}(u; Q) \gtrsim C_{N,p,\gamma}^{\mathrm{cell}} \int_Q |\nabla u|^p dx. $$

通过累加所有小立方体的贡献，并利用紧性论证，得到全局下界。

对于 $p=1$ 的情况，难度显著增加。因为 $BV$ 函数的归一化缩放未必收敛到仿射函数，它们可能收敛到一般的单调递增的一维 $BV$ 函数。这是证明的最大技术难点之一。作者必须证明：对于任意此类一维单调 $BV$ 函数 $v$，其胞腔能量满足

$$ G_{\lambda,1,\gamma}(v; Q_0) \ge C_{N,1,\gamma}^{\mathrm{cell}} |Dv|(Q_0). $$

为此，作者在 Section 4.2 和 4.3 中详细研究了缩放后 $BV$ 函数的一维极限行为，并建立了一个关于单调 $BV$ 函数的下界引理。这一部分避免了使用切片法，而是直接通过分析 $BV$ 函数的精细结构（fine structure）来处理跳跃集和绝对连续部分的能量贡献，确保了即使在没有仿射极限的情况下，能量依然被总变差乘以正确的常数所控制。

2. 上极限不等式的证明策略（The Gamma-Limsup Inequality）

构造恢复序列是另一个关键挑战。由于泛函的非局部性，直接将不同小立方体上的最小化序列拼接在一起会导致在立方体边界处产生额外的非局部相互作用能量（interaction energy）。为了消除这种边界效应，作者采用了受 Nguyen [25] 启发的“边界修改”（boundary modification）技术。

具体步骤如下：

  • Step 1: 首先，对于实现胞腔常数的极小化序列，构造特殊的修正序列 $V_j$。这些序列在靠近立方体边界 $\partial Q_0$ 的一个薄层（厚度为 $\rho_j$）内被强制修改为与仿射函数 $\ell$ 完全一致。这保证了在边界处没有不连续性或剧烈变化。
  • Step 2: 利用截断函数和光滑过渡函数（如文中定义的 $\theta$ 和 $\eta_j$），确保修改后的序列 $V_j$ 满足特定的衰减条件：
    $$ \lim_{j\to\infty} \Lambda_j \rho_j^\gamma = \infty \quad \text{且} \quad \lim_{j\to\infty} \frac{\|V_j - \ell\|_{L^\infty}^p}{\Lambda_j \rho_j^{p+\gamma}} = 0. $$
    这些条件确保了边界层内的额外能量贡献随着 $\lambda_j \to \infty$ 而消失。
  • Step 3: 对于一般的连续分段仿射函数，将其分解为若干个小立方体，在每个立方体上放置上述修正后的缩放序列。由于边界处已匹配仿射函数，不同块之间的非局部相互作用被有效抑制。
  • Step 4: 最后，利用稠密性论证，将结果推广到一般的 $W^{1,p}$ 和 $BV$ 函数。对于 $BV$ 函数，需要结合 $L^1$ 收敛性和总变差的半连续性进行逼近。

3. 核心工具总结

  • 胞腔公式（Cell Formula）：将全局极限常数局部化，转化为单位立方体上的极小化问题。
  • Morse 覆盖引理：用于处理 $p>1$ 时的局部化估计。
  • BV 函数的一维极限分析：专门针对 $p=1$ 时非仿射极限的处理。
  • 边界修正技术（Boundary Modification）：解决非局部泛函拼接时的界面能量问题。
  • 凸性不等式与特征函数估计：在处理 $G_{\lambda,p,\gamma}$ 的定义式时，利用 $|t|^p$ 的凸性和特征函数的集合包含关系进行精细的能量拆分。

4. 与之前工作的比较

本文与既有文献相比，在多个维度上实现了显著的推进和突破：

1. 问题解决范围的扩展：从全空间到有界域 Brezis 等人在早期工作（如 [9], [10]）中主要关注全空间 $\mathbb{R}^N$ 上的渐近行为。虽然他们在 $\mathbb{R}^N$ 上建立了点态极限和部分 $\Gamma$-收敛结果，但对于有界域 $\Omega$ 上的 $\Gamma$-收敛问题，特别是强 $L^p(\Omega)$ 拓扑下的收敛，一直是一个开放问题（Brezis Open Problem 9.3）。Nguyen [24, 25] 和 Gobbino-Picenni [20] 的工作主要集中在 $\mathbb{R}^N$ 或特定参数（如 $\gamma < 0$ 或 $\gamma \in (-1, 0)$）的情形。Lahti 和 Li [23] 虽然在 $p=1$ 时研究了 $\Gamma$-极限，但其收敛拓扑是面积严格收敛（area-strict convergence），而非本文关注的强 $L^p$ 拓扑。本文首次在有界 Lipschitz 域上，针对 $\gamma \in (0, \infty)$ 的所有情形，完整解决了强 $L^p(\Omega)$ 拓扑下的 $\Gamma$-收敛问题，给出了肯定的回答。

2. 参数范围的推广与常数精度的提升 既往工作多集中在 $\gamma = N$ 或特定的 $\gamma$ 值。本文不仅涵盖了 $\gamma = N$ 的特例，还将结论推广到了任意 $\gamma \in (0, \infty)$。更重要的是，本文精确刻画了极限泛函中的比例常数 $C_{N,p,\gamma}^{\mathrm{cell}}$。之前的文献（如 Nguyen [25]）虽然指出了常数存在且由隐式胞腔公式给出，但未能在有界域上明确其与 BBM 常数的关系及具体的上下界。本文通过 Lemma 2.5 给出了显式的上下界估计：

$$ \frac{K_{N,p}}{\log 2 \cdot \frac{2^{\gamma+1}-1}{\gamma}} \le C_{N,p,\gamma}^{\mathrm{cell}} \le \frac{K_{N,p}}{\gamma}. $$

这一结果澄清了极限常数并非简单的 BBM 常数，而是包含了域几何和边界效应的修正因子。

3. 处理 $p=1$ 时 BV 空间的独特性 在 $p=1$ 的情形下，点态极限可能不存在（Lahti [22] 的反例），这使得直接推导变得极其困难。Nguyen [26] 曾指出对于某些 $\gamma$ 范围极限为零，而 Lahti 和 Li [23] 仅在面积严格收敛意义下获得结果。本文突破了这一局限，证明了即使在点态极限缺失的情况下，$\Gamma$-极限依然存在且具有良好的结构（正比于总变差）。这一结论依赖于作者对 $BV$ 函数缩放极限的深刻分析——即证明了即使缩放后收敛到非仿射的单调 $BV$ 函数，其能量下界依然由相同的胞腔常数控制。这是对 $p=1$ 情形最本质的推进，弥补了之前文献在强拓扑下 $BV$ 空间分析的空白。

4. 方法论的创新：避免切片法与改进边界处理 与 Gobbino 和 Picenni [20] 使用的二进分割和切片论证不同，本文采用了更直接的“胞腔比较”方法和“边界修正”技术。这种方法论创新使得证明更加紧凑，且能更好地适应有界域的几何限制。特别是在处理非局部相互作用导致的边界能量溢出问题上，本文提出的构造修复序列的方法（在边界附近强制匹配仿射函数）比传统的全局延拓或截断更为有效，成功克服了非局部算子在拼接时的能量损失问题。

代价与局限 本文的证明依赖于域 $\Omega$ 为 Lipschitz 域（或 $N=1$ 时的开区间），以保证边界的良好行为和 Sobolev/BV 函数的迹性质。对于更不规则的域（如分形边界），结论是否成立尚需进一步探讨。此外，虽然给出了常数的上下界，但 $C_{N,p,\gamma}^{\mathrm{cell}}$ 的精确闭式表达（closed-form expression）对于一般 $N, p, \gamma$ 仍未给出，这仍是未来研究的一个潜在方向。总体而言，本文在假设条件和结论强度之间取得了很好的平衡，是目前该领域最完备的结果。

8. 超定奇异问题与分数扭转

Overdetermined Singular Problems and Fractional Torsion

作者: Daniel Baratta

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18838v1

摘要翻译: 我们研究有界开集上带有孤立不可去内部奇点的分数拉普拉斯算子的Serrin型超定问题。对于方程 $(- ext{Delta})^s u=f(u)$ 在 $ ext{Omega}ackslash ext{0}$ 中、边界条件 $u=0$ 在 $ ext{R}^N ext{Omega}$ 上以及法向导数条件 $( ext{partial}_ ext{eta})^s u=c$ 在 $ ext{partial Omega}$ 上的正弱解，其中 $0<s<1$, $N>2s$ 且 $c<0$，我们分离出边界商 $u/ ext{delta}^s$ 的一阶切向抵消，该抵消足以闭合分数角点论证所需的一阶抵消，且弱于要求在完整边界邻域内 $u/ ext{delta}^s ext{ in } C^1$。在此条件下，$ ext{Omega}$ 是以奇异点为中心的球，且 $u$ 关于径向变量是径向的且严格递减。极点处无需假设逐点爆破速率。若 $ ext{Omega}$ 属于 $C^{2, ext{alpha}}$ 类，则当 $s>1/2$ 时该抵消自动成立；当 $f$ 在零附近为常数时，对任意 $s ext{ in } (0,1)$ 均成立。特别地，在扭转情形 $f=1$ 下，$u(x)= ext{tau}_ ext{R}(x)+kG_ ext{R}(x,0)$，其中 $ ext{tau}_ ext{R}$ 和 $G_ ext{R}$ 分别是球的分数扭转函数和Green函数。

Comment: 33 pages

链接: http://arxiv.org/abs/2608.18838v1

1. 研究对象和背景

本文研究的是分数拉普拉斯算子 $(-\Delta)^s$ 在具有孤立不可去内奇点的有界开集上的 Serrin 型过定问题。具体而言，考虑如下边值问题：

$$ \begin{cases} (-\Delta)^s u = f(u) & \text{在 } \Omega \setminus \{0\} \text{ 中}, \\ u > 0 & \text{在 } \Omega \setminus \{0\} \text{ 中}, \\ u = 0 & \text{在 } \mathbb{R}^N \setminus \Omega \text{ 中}, \\ (\partial_\eta)^s u = c & \text{在 } \partial\Omega \text{ 上}, \end{cases} $$

其中 $0 < s < 1$, $N > 2s$, $c < 0$，且 $\Omega \subset \mathbb{R}^N$ 是包含原点 $0$ 的 $C^2$ 类有界开集。非线性项 $f \in \mathrm{Lip}_{\mathrm{loc}}(\mathbb{R}_+)$。

该问题的核心在于处理非局部算子下的刚性（rigidity）现象。经典的 Serrin 定理 [29] 指出，若光滑有界域 $\Omega$ 上存在满足 $-\Delta u = 1$、$u|_{\partial\Omega}=0$ 且法向导数 $\partial_\nu u$ 为常数的正解，则 $\Omega$ 必为球体。这一结果依赖于移动平面法（moving planes method）、极大值原理、Hopf 引理及角点引理。然而，当引入内部奇点（即解在原点处发散）时，传统的移动平面论证面临严峻挑战：反射过程中奇点的镜像会破坏测试函数的正则性，且非局部算子的积分定义使得在奇点附近的能量估计变得复杂。

近年来，Fall 和 Jarohs [20] 将 Serrin 定理推广至分数拉普拉斯算子，建立了无奇点情形下的刚性理论，其证明依赖于反对称函数的极大值原理、分数 Hopf 引理及分数角点引理。随后，Esposito, Sciunzi 和 Soave [19] 研究了局部椭圆方程中的奇异对应物，证明了在存在不可去极点且满足过定条件时，域必须是以极点为中心的球体。然而，针对分数阶情形，目前尚无文献处理“解本身具有不可去孤立内奇点”且“在边界上给定常数分数法向导数”的有界 Dirichlet 问题。本文旨在填补这一空白，建立非局部框架下的奇异刚性定理。

主要难点在于：

  1. 奇点处理：由于 $N > 2s$，单点在分数 $s$-容量意义下为零容量，这允许使用零容量截断函数移除原始奇点及其反射像，但需要精确控制由此产生的非局部交换子误差。
  2. 边界行为：在移动平面法的最后一步（排除正交接触），需要利用分数角点引理。该引理给出的下界阶数为 $t^{1+s}$。为了闭合论证，必须确保边界商 $u/\delta^s$ 在切向方向上具有特定的一阶抵消性质。传统的 $C^1$ 正则性假设过于严格，本文试图找到更弱的充分条件。
  3. 解的结构：在确定域为球体后，需对解进行显式分类，特别是在扭转情形（$f \equiv 1$）下，解应能分解为分数扭转函数与 Green 函数的线性组合。

本文的研究不仅推进了分数 PDE 的对称性理论，还揭示了非局部算子在奇异扰动下的精细边界行为，为理解分数扭转问题的物理模型（如反常扩散中的应力集中）提供了严格的数学基础。

2. 主要定理或主要结果

本文的核心成果建立在以下三个主要定理之上，均针对满足上述问题设定 (1.1) 的正弱解 $u$。

定理 1.1 (奇异刚性定理): 设 $0 < s < 1$ 且 $N > 2s$。令 $\Omega \subset \mathbb{R}^N$ 为包含原点的 $C^2$ 类有界开集。设 $f \in \mathrm{Lip}_{\mathrm{loc}}(\mathbb{R}_+)$，且 $u$ 是问题 (1.1) 的一个正弱解，其在原点处具有不可去奇点（定义为 Definition 2.7）。假设解满足正则性条件：

$$ u \in W^{s,2}_{\mathrm{loc}}(\mathbb{R}^N \setminus \{0\}) \cap L^1_s(\mathbb{R}^N) \cap C(\mathbb{R}^N \setminus \{0\}). $$

此外，假设边界商 $\psi := u/\delta^s$ 满足一阶切向抵消条件 (TC)：对于任意 $Q \in \partial\Omega$ 和单位切向量 $\tau \in T_Q\partial\Omega$，有

$$ \psi(Q + t(\nu_Q + \tau)) - \psi(Q + t(\nu_Q - \tau)) = o(t) \quad \text{当 } t \to 0^+ \text{ 时}. $$

在此条件下，结论成立：存在 $R > 0$ 使得 $\Omega = B_R(0)$，且 $u$ 关于径向变量 $|x|$ 是径向对称且严格递减的。 注：该定理未要求 $\Omega$ 连通，连通性由移动平面法导出的凸性自动推出。

推论 1.2 (自动正则性范围): 在定理 1.1 的假设下，若额外假设 $\Omega \in C^{2,\alpha}$ ($\alpha \in (0,1)$)，则切向抵消条件 (TC) 在以下两种情形下自动成立，从而无需作为独立假设：

  1. 半空间情形：当 $s > 1/2$ 时；
  2. 局部常数非线性：当存在 $a > 0$ 使得 $f$ 在 $[0, a]$ 上为常数时。 在这两种情况下，同样可推出 $\Omega$ 为以原点为中心的球体，且 $u$ 径向对称严格递减。

推论 1.3 (扭转情形的显式分类): 设 $0 < s < 1$, $N > 2s$，且 $\Omega \in C^{2,\alpha}$。考虑扭转情形 $f \equiv 1$。若 $u$ 是问题 (1.1) 的解且具有不可去奇点，则无需任何额外的切向抵消假设，必有 $\Omega = B_R(0)$，且解具有如下显式分解形式：

$$ u(x) = \tau_R(x) + k G_R(x, 0), $$

其中：

  • $\tau_R(x) = \gamma_{N,s}(R^2 - |x|^2)_+^s$ 是球体 $B_R$ 上的分数扭转函数；
  • $G_R(x, 0)$ 是 $B_R$ 上分数拉普拉斯算子的 Green 函数（基于 Boggio 公式）；
  • 系数 $k$ 由边界常数 $c$ 唯一确定为：
    $$ k = \frac{s R^{N-s}}{c_{N,s} 2^s} (-c - \gamma_{N,s} (2R)^s)^{-1} > 0, $$
    其中 $\gamma_{N,s}$ 和 $c_{N,s}$ 分别为扭转函数和 Green 函数表示中的正常数。 此外，解在原点附近的渐近行为由下式给出：
    $$ \lim_{x \to 0} |x|^{N-2s} u(x) = k c_{N,s} B(s, N/2 - s) = k c_{N,s} \frac{\Gamma(s)\Gamma(N/2-s)}{\Gamma(N/2)}. $$
    特别地，若允许奇点可去（即 $k=0$），则 $u = \tau_R$，此时 $-c = \gamma_{N,s}(2R)^s$。

这些结果明确给出了在何种几何和解析条件下，分数奇异过定问题强制域为球体，并给出了具体的解结构公式。

3. 研究方法、关键技术和核心工具

本文的证明框架建立在 Fall-Jarohs [20] 的非局部移动平面法之上，但针对奇异情形进行了关键性的修正与创新。主要技术障碍在于如何处理移动平面过程中产生的反射奇点，以及如何克服非局部算子在角点处的低阶估计限制。作者的创新点主要体现在以下三个方面：

1. 零容量截断策略与能量估计的精细控制 在标准的移动平面法中，比较函数通常直接取 $w_\lambda = u_\lambda - u$。但在本题中，$u$ 在原点 $0$ 处有奇点，其反射像 $0_\lambda$ 也可能落在移动区域 $\Sigma_\lambda$ 内。直接测试会导致能量泛函发散或定义不明确。 作者采用了 Montoro-Punzo-Sciunzi [26] 提出的零容量截断策略。由于 $N > 2s$，单点的分数 $s$-容量为零，即 $\mathrm{Cap}_s(\{0\}) = 0$。作者构造了一族反射不变的截断函数 $\psi_{a,\epsilon} = 1 - \eta_{a,\epsilon}$，其中 $\eta_{a,\epsilon}$ 是围绕奇点 $a$ 的光滑截断核。通过选取合适的 $\epsilon$，构造测试函数 $\phi = w_\lambda^- \cdot \psi_{0,\epsilon} \cdot \psi_{0_\lambda,\epsilon}$。 关键在于估计截断带来的非局部交换子误差。作者详细计算了形如 $E(w_\lambda, \phi)$ 的双线性型中涉及 $\psi$ 差分的项，证明了随着 $\epsilon \to 0$，这些误差项趋于零的速度快于主项，且该收敛性独立于 $u$ 在奇点附近的具体爆破速率。这一技巧成功地将原本在奇点处无定义的比较问题转化为在全局能量空间中有意义的测试问题，从而能够应用反对称函数的极大值原理。

2. 一阶切向抵消条件 (TC) 的引入与角点引理的匹配 这是本文最核心的分析创新。在移动平面法的最后阶段，需要排除平面与边界正交接触的可能性。根据 Fall-Jarohs 的分数角点引理，若在正交接触点附近，解的某种差值增长阶数低于 $t^{1+s}$，则导出矛盾。 具体而言，分数角点引理给出的下界是 $O(t^{1+s})$。由于 $u = \delta^s \psi$，而 $\delta \sim t$，因此 $u$ 的边界贡献为 $O(t^s)$。为了使反射后的差值 $w_\lambda$ 能够产生严格小于 $t^{1+s}$ 的上界（从而触发矛盾），边界商 $\psi = u/\delta^s$ 必须在切向方向上表现出比 $O(1)$ 更高的精度。 作者定义了条件 (TC)：

$$ \psi(Q + t(\nu_Q + \tau)) - \psi(Q + t(\nu_Q - \tau)) = o(t). $$

这一条件意味着 $\psi$ 在切向的一阶导数在边界上是“抵消”的。如果仅假设 $\psi \in C^\alpha$ ($\alpha < 1$)，则差值为 $O(t^\alpha)$，无法保证 $o(t)$ 的阶数。如果假设 $\psi \in C^1$，虽然满足 (TC)，但过于苛刻。Appendix A 中构造的反例表明，即使对于线性非线性 $f(u)=-u$，也存在解使得 $u/\delta^s \notin C^1$ 但仍满足 (TC)。因此，(TC) 是关闭角点论证所需的最小正则性条件，它比全 $C^1$ 正则性弱，但足以匹配 $t^{1+s}$ 的角点阶数。

3. 奇异部分的分布论分类与显式构造 一旦通过对称性证明 $\Omega$ 为球体，问题转化为在球体上求解带奇异源的分数方程。作者利用了 Li-Wu-Xu [25] 关于分数 Bôcher 定理的分布论版本。该定理允许在分布意义下描述解在原点附近的奇异性。 结合 Boggio 公式给出的球体 Green 函数 $G_R$ 和分数扭转函数 $\tau_R$ 的显式表达，作者证明了任何满足条件的解 $u$ 都可以唯一分解为 $u = \tau_R + k G_R$。这里的关键是利用了 $f \equiv 1$ 时的叠加原理（在分布意义下），以及 Green 函数在奇点处的精确渐近行为 $|x|^{2s-N}$。通过匹配边界上的常数法向导数 $c$，作者推导出了系数 $k$ 的显式公式，完成了从定性对称到定量结构的跨越。

整个证明框架逻辑严密：首先通过截断技术解决奇点带来的定义困难，然后利用 (TC) 条件完成对称性证明，最后利用分数势论工具完成解的分类。这种方法不仅适用于分数 Laplacian，也为其他非局部算子的奇异问题提供了范式。

4. 与之前工作的比较

本文与既有文献相比，在以下几个方面取得了显著的推进和突破：

1. 从非奇异到奇异情形的跨越 Fall 和 Jarohs [20] 建立的分数 Serrin 定理是针对无奇点情形（$u \in H^s_0(\Omega)$）的。他们证明了若存在正解满足过定条件，则域为球体。本文首次将这一理论扩展至具有不可去孤立内奇点的情形。在局部椭圆方程中，Esposito, Sciunzi 和 Soave [19] 已经处理了此类奇异问题，证明了刚性。本文是非局部算子领域对该结果的直接推广。主要区别在于：局部情形下，奇点附近的正则性相对容易处理，且移动平面法中的截断较为简单；而在分数情形下，由于非局部性，奇点的影响会通过长程相互作用传播到整个域，必须引入复杂的零容量截断技术和精细的能量估计来隔离奇点影响。本文填补了非局部奇异刚性理论的空白。

2. 边界正则性假设的弱化 在 Fall-Jarohs [20] 的后续讨论及修订版中，为了确保角点引理的有效性，往往隐含或明确要求边界商 $u/\delta^s$ 具有 $C^1$ 正则性（见 [20, Remark 1.3]）。本文通过引入条件 (TC)，显著降低了这一要求。(TC) 仅要求 $u/\delta^s$ 在切向的一阶差分是 $o(t)$，而非整体 $C^1$。Appendix A 中的反例明确构造了一个解，其 $u/\delta^s$ 在边界处不连续于 $C^1$ 空间（表现为 $A - B\delta^{2s}$ 形式的展开，其中 $2s < 1$），但仍满足 (TC) 并导致对称性。这意味着本文的结果适用于更广泛的解类，特别是那些在边界处具有较弱正则性的奇异解。这一改进对于处理更一般的非线性 $f$ 至关重要。

3. 适用范围与自动正则性的揭示 本文不仅提出了 (TC) 条件，还进一步分析了该条件何时自动成立。Corollary 1.2 指出，当 $s > 1/2$ 或 $f$ 在零点附近为常数时，(TC) 自动满足。这一发现具有重要的分类意义：

  • 对于 $s > 1/2$ 的高阶分数情况，边界正则性自然足够强，无需额外假设；
  • 对于 $f \equiv 1$ 的扭转问题，无论 $s \in (0,1)$ 如何，刚性都自动成立。 相比之下，之前的文献（如 Fall-Jarohs）在处理一般 $s$ 和一般 $f$ 时，往往需要更强的先验假设或无法处理奇异情形。本文明确了不同参数区间下的理论完备性，展示了分数阶参数 $s$ 对边界行为的重要影响。

4. 解的显式结构与常数依赖 在局部奇异 Serrin 问题中，解的结构通常被描述为正则部分加奇异部分。本文在推论 1.3 中给出了分数扭转情形下解的显式分解公式 $u = \tau_R + k G_R$，并精确计算了系数 $k$ 与边界数据 $c$ 之间的代数关系。这种显式刻画不仅验证了理论的一致性，还为数值模拟和物理应用提供了精确的基准解。此前的非局部文献多集中于存在性和对称性，较少给出如此详细的显式分类。

局限性与代价 本文的代价在于引入了 (TC) 条件作为一般情形的必要假设。虽然该条件比 $C^1$ 弱，但它仍然是一个额外的解析约束，并非对所有可能的 $f$ 和 $s$ 都能自动满足（例如 $s < 1/2$ 且 $f$ 非线性变化较快时）。此外，证明过程高度依赖于 $N > 2s$ 这一容量条件，这使得结果无法直接推广到临界或超临界情形（$N \le 2s$），因为在那些情形下单点具有非零容量，截断策略失效。尽管如此，本文在现有框架下实现了最优的理论推进，为非局部奇异 PDE 的研究树立了新的标准。

9. $Φ^*[I]$的正则性与非退化性蕴含固定边界$ackslash partialΩ$的正则性

Regularity and non-degeneracy of $Φ^*[I]$ implies regularity of the fixed boundary $\partialΩ$

作者: Henrik Garde, Michael S. Vogelius

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18400v1

摘要翻译: 对于定义域$ackslash overlineΩ$到自身的一个微分同胚$Φ$，若其在$ackslash partialΩ$上恒等，我们证明：在满足某个非退化条件的前提下，推前像$Φ^*[I]$的局部正则性蕴含$ackslash partialΩ$的局部正则性。具体而言，若$ν$是点$P$处$ackslash partialΩ$的法向量，则条件$(Φ^*I-I)ν

Comment: 8 pages

链接: http://arxiv.org/abs/2608.18400v1

1. 研究对象和背景

本文研究的核心对象是定义在 $\mathbb{R}^n$ ($n \geq 2$) 中某个有界区域 $\Omega$ 上的微分同胚映射 $\Phi: \overline{\Omega} \to \overline{\Omega}$，该映射满足边界恒等条件 $\Phi(x) = x$ 对所有 $x \in \partial\Omega$ 成立，且雅可比行列式 $\det D\Phi > 0$。在此设定下，作者考察了单位矩阵 $I$ 的推前向（push-forward）张量场，其定义为：

$$ \Phi^*[I](y) = (D\Phi D\Phi^T \circ \Phi^{-1})(y) / \det(D\Phi \circ \Phi^{-1})(y), \quad y \in \Omega. $$

更一般地，对于任意光滑对称正定矩阵值系数 $A$，其推前向定义为 $\Phi^*[A] = D\Phi A D\Phi^T \circ \Phi^{-1} / \det(D\Phi \circ \Phi^{-1})$。这一构造在电磁成像理论中具有关键地位，具体体现在三个方面：

  1. 电导层析成像（Electrical Impedance Tomography, EIT）：介质 $\Phi^*[A]$ 产生的边界测量数据（Cauchy 数据）与原始介质 $A$ 完全相同。在二维情形下，这是具有该性质的唯一介质类。
  2. 亥姆霍兹散射（Helmholtz Scattering）：若背景电导率为 $A$，则 $\Phi^*[A]$ 代表了一类对任意入射波均不产生散射的非散射介质（non-scattering inhomogeneities）。
  3. 透射特征值问题（Transmission Eigenvalue Problem）：对于背景电导率 $A$，$\Phi^*[A]$ 使得所有数值均为透射特征值，这通常违反了标准适定性所需的“覆盖条件”。

尽管上述性质在物理上看似与一般理论矛盾（例如各向同性介质的唯一性或谱的离散性），但通过李乌维尔定理或特定的边界条件破坏可以解释这些兼容性。然而，一个核心且尚未完全解决的问题是：在非散射背景下，如果介质 $\Phi^*[I]$ 具有某种正则性，是否意味着散射体（即区域 $\Omega$）的边界 $\partial\Omega$ 也必须具有相应的正则性？

已有的文献（如 Cakoni, Vogelius, Xiao [3]）证明了：若 $A$ 属于 $C^{k+1,\alpha}$ 类，且存在入射波梯度 $\xi$ 使得 $\xi^T(A-I)\nu \neq 0$（其中 $\nu$ 为边界法向量），则 $\partial\Omega$ 必须是 $C^{k+1,\alpha}$ 类的。这一结果后来被推广至仅假设 $\Omega$ 为 Lipschitz 区域的情形 [6]。为了理解 $\Phi^*[I]$ 作为非散射介质的相容性，必须阐明：为何当 $\Phi^*[I]$ 本身足够光滑时，$\partial\Omega$ 不能出现奇点？

本文试图填补的空白正是这一逻辑链条中的关键环节：在 $\Phi$ 仅为 $C^{1,\alpha}$ 微分同胚且 $\partial\Omega$ 初始假设为 $C^{1,\alpha}$ 的前提下，证明 $\Phi^*[I]$ 的局部 $C^{k+1,\alpha}$ 正则性结合适当的非退化条件，足以推出 $\partial\Omega$ 的 $C^{k+1,\alpha}$ 正则性。这直接回答了非散射反演问题中关于几何形状正则性的逆问题，并澄清了为何非散射介质无法在具有尖角或低正则性边界的区域上存在（除非满足特定的退化条件）。

2. 主要定理或主要结果

本文的主要成果建立在两个非退化假设之上，分别对应不同的技术路径和适用范围。设 $P \in \partial\Omega$，$\nu$ 为该点处的单位外法向量。

主要定理（Theorem 2.2）： 假设在点 $P$ 处满足非退化条件 Assumption 2.1：

$$ (\Phi^*[I](P) - I)\nu \neq 0, $$

等价于 $(D\Phi(P) - I)\nu \neq 0$。若进一步假设 $\Phi^*[I]$ 在 $P$ 附近的 $\Omega$ 内属于 $C^{k+1,\alpha}$ 类（其中 $k \geq 1$），则边界 $\partial\Omega$ 在 $P$ 附近也属于 $C^{k+1,\alpha}$ 类。

推论（Corollary 2.3）： 在相同的非退化条件下，若 $D\Phi D\Phi^T \circ \Phi^{-1}$ 在 $P$ 附近属于 $C^{k+1,\alpha}$ 类，同样可推出 $\partial\Omega$ 为 $C^{k+1,\alpha}$ 类。对于 $n \geq 3$，由于 $\det \Phi^*[I] = (\det D\Phi)^{2-n} \circ \Phi^{-1}$，该推论与主定理等价。

辅助定理（Theorem 2.5）： 引入更强的非退化条件 Assumption 2.4：存在切向量 $\xi$ 使得 $\xi^T(D\Phi(P) - I)\nu \neq 0$。在此条件下，证明更为直接。设局部坐标系下 $P=0$，边界由 $x_n = \psi(\tilde{x})$ 描述，$\psi \in C^{1,\alpha}$。结论如下：

  1. 若映射 $\tilde{x} \mapsto D\Phi D\Phi^T(\tilde{x}, \psi(\tilde{x}))$ 在 $0$ 附近属于 $C^{k,\alpha}$ 类，则 $\partial\Omega$ 属于 $C^{k+1,\alpha}$ 类。
  2. 若映射 $\tilde{x} \mapsto \Phi^*[I](\tilde{x}, \psi(\tilde{x}))$ 在 $0$ 附近属于 $C^{k,\alpha}$ 类，则 $\partial\Omega$ 属于 $C^{k+1,\alpha}$ 类。

技术细节与尺度关系：

  • 正则性提升阶数：从 $\Phi^*[I]$ 的 $C^{k+1,\alpha}$ 推导至 $\partial\Omega$ 的 $C^{k+1,\alpha}$，本质上是通过求解一个非线性椭圆型自由边界问题实现的。证明过程依赖于 $\psi$ 的一阶导数 $\nabla\psi$ 与 $D\Phi$ 分量之间的代数关系，以及后续利用 Hodograph 变换将自由边界转化为固定边界问题。
  • 关键不等式：在证明中，线性化算子的椭圆性依赖于估计式 $\min_{|\xi|=1} \xi \cdot A_0 \xi \geq c |\xi|^2$，其中 $c > 0$ 是由 $\Phi^*[I]$ 的正定性和非退化条件保证的常数。这确保了 Agmon-Douglis-Nirenberg 边界正则性理论的可应用性。
  • 递归论证：对于 $k \geq 1$ 的情况，证明采用递归策略。首先利用 $\Phi^*[I] \in C^{2,\alpha}$ 导出 $\partial\Omega \in C^{2,\alpha}$，进而将 $\Phi^*[I]$ 视为更高阶光滑系数的函数，逐步提升至 $C^{k+1,\alpha}$。

3. 研究方法、关键技术和核心工具

本文的证明面临的核心障碍在于：$\Phi$ 本身未必比 $\Phi^*[I]$ 更光滑，甚至 $\Phi^*[I]$ 的光滑性并不直接蕴含 $\Phi$ 的光滑性（因为 $\Phi^*[I]$ 仅包含 $D\Phi$ 的二次型信息，丢失了旋转部分的信息）。因此，不能简单地通过反解 $\Phi$ 来获得边界的正则性。作者根据非退化条件的强弱，采用了两种截然不同的分析框架。

第一部分：基于代数反演的初等方法（针对 Assumption 2.4） 当满足 Assumption 2.4（即存在切向量 $\xi$ 使得 $\xi^T(D\Phi - I)\nu \neq 0$）时，证明主要依赖于直接的代数计算和隐函数定理。作者首先在边界 $\partial\Omega$ 上建立 $D\Phi$ 的分量与边界法线方向 $\nabla\psi$ 之间的显式关系。在局部坐标系下，利用 $\Phi$ 在边界上的恒等性质 $\Phi(x) = x$，导出 $D\Phi$ 的切向导数与 $\psi'$ 的关系。具体而言，作者构造了一个映射 $F: \mathbb{R}^{n-1} \times \mathbb{R}^{n-1} \times \mathbb{R} \to \mathbb{R}^{(n-1)^2+n}$，将边界斜率 $\tilde{\eta} = \nabla\psi$、$D\Phi$ 的法向导数分量 $\tilde{\xi} = \partial_n e\Phi$ 和 $\xi_n = \partial_n \phi_n$ 映射到已知的 $D\Phi D\Phi^T$ 的边界迹上。

关键在于验证该映射 $F$ 的线性化算子 $DF$ 在特定点处的单射性。通过计算发现，$DF$ 的行列式非零当且仅当 $\tilde{\xi} \neq 0$ 且 $\xi_n \neq 0$，这正是 Assumption 2.4 所保证的。一旦确认 $DF$ 可逆，利用隐函数定理即可从 $D\Phi D\Phi^T$ 的光滑性反解出 $\nabla\psi$ 的光滑性，从而直接得到 $\partial\Omega$ 的 $C^{k+1,\alpha}$ 正则性。这种方法避免了复杂的 PDE 理论，仅利用了微积分中的隐函数定理和矩阵代数。

第二部分：基于 Hodograph 变换的自由边界问题（针对 Assumption 2.1 但非 2.4） 当 Assumption 2.1 成立但 Assumption 2.4 不成立时，意味着 $\nu$ 是 $D\Phi(P)$ 的特征向量且特征值不为 1。此时，上述代数方法失效，因为 $D\Phi$ 在切方向上的信息不足以唯一确定边界曲率。作者转而采用类似于 Kinderlehrer 和 Nirenberg [4] 以及 Cakoni 等人 [3] 的策略，引入 Hodograph 变换将自由边界问题转化为固定边界问题。

具体步骤如下：

  1. 构造辅助函数：定义 $v(y) = (\Phi^{-1}(y))_n - y_n$（在二维情形下简化为 $x_2 - y_2$）。该函数在 $\partial\Omega$ 上为零，且满足一个带有非齐次项的椭圆方程 $\nabla_y \cdot (\Phi^*[I] \nabla_y v) = -\nabla_y \cdot (\Phi^*[I] e_2)$。
  2. 边界条件分析：利用 Assumption 2.1，作者在边界上验证了 $v$ 的法向导数 $\partial_2 v(P) \neq 0$。这是应用 Hodograph 变换的关键前提，保证了变换的局部双射性。
  3. 坐标变换：引入变换 $H: (y_1, y_2) \mapsto (z_1, z_2) = (y_1, v(y_1, y_2))$，将 $\Omega$ 的邻域映射到上半平面 $z_2 > 0$。原边界 $\partial\Omega$ 变为新坐标系下的平面 $z_2 = 0$。
  4. 导出新的 PDE 系统：在新变量下，原问题转化为关于未知函数 $w(z)$（即 $\partial\Omega$ 的反函数图像）的非线性椭圆型偏微分方程组。该方程组的形式为 $\partial_{z_1} a_1 + \partial_{z_2} a_2 + a_0 = 0$，并伴随一个斜导数边界条件（oblique derivative condition）。
  5. 线性化与正则性提升：对系统进行线性化，证明线性化算子在 $C^{k+1,\alpha}$ 空间上是一致椭圆且满足 Agmon-Douglis-Nirenberg 边界条件的。利用已知结果（[1, Theorem 11.2]），从系数 $\Phi^*[I]$ 的光滑性推导出解 $w$ 的光滑性，进而恢复 $\partial\Omega$ 的 $C^{k+1,\alpha}$ 正则性。

创新点总结：本文的创新之处在于系统地分类讨论了非退化条件的不同情形，特别是处理了 $D\Phi$ 沿法向拉伸但不沿切向剪切（即 Assumption 2.4 失效）的特殊情况。通过巧妙结合代数反演和 Hodograph 变换，作者在不假设 $\Phi$ 本身高阶光滑的情况下，成功建立了推前向张量正则性与几何边界正则性之间的严格联系。

4. 与之前工作的比较

本文工作与既有文献（特别是 [3] 和 [6]）在以下维度存在显著差异与推进：

1. 问题设定的互补性： 文献 [3] 主要关注的是非散射介质 $A$ 的逆问题，即给定 $A$ 的光滑性和非散射条件，推导边界 $\partial\Omega$ 的正则性。而本文专门针对 $\Phi^*[I]$ 这一特殊形式的介质（即由恒等映射生成的非散射介质）。虽然 [3] 的结果在形式上涵盖了 $\Phi^*[I]$，但其证明依赖于入射波梯度 $\xi$ 的具体结构。本文则直接从 $\Phi^*[I]$ 的几何性质出发，证明了只要满足 $(\Phi^*[I]-I)\nu \neq 0$，边界正则性即成立。这使得结论更加纯粹，不依赖于外部入射波的特定方向，而是内蕴于介质本身的几何结构中。

2. 假设条件的细化与推广： 文献 [6] 将正则性结果推广到了 Lipschitz 边界，但并未详细讨论 $\Phi^*[I]$ 在 $C^{1,\alpha}$ 边界上的具体行为机制。本文在假设 $\Omega$ 初始为 $C^{1,\alpha}$ 的基础上，精确刻画了 $\Phi^*[I]$ 的正则性如何“传递”回边界。特别是，本文区分了两种非退化情形：

  • 情形 A（Assumption 2.4）：存在切向扰动。此时边界正则性可通过简单的代数反演获得，这与 [3] 中的某些直觉相符，但本文给出了更直接的证明。
  • 情形 B（仅 Assumption 2.1）：无切向扰动，仅有法向拉伸。这是 [3] 未明确处理的边缘情况。本文通过 Hodograph 变换填补了这一空白，证明了即使在没有切向剪切的情况下，只要法向拉伸非退化，边界依然保持光滑。

3. 适用范围与代价： 本文的结论要求 $\Phi$ 至少是 $C^{1,\alpha}$ 的微分同胚，且 $\Phi^*[I]$ 需达到 $C^{k+1,\alpha}$。相比之下，[3] 的处理可能允许更弱的初始正则性假设，但在处理 $\Phi^*[I]$ 这种特殊结构时，本文的方法更为精确和自洽。代价是本文无法处理 $\Phi^*[I]$ 退化的情形（即 $(\Phi^*[I]-I)\nu = 0$），此时边界可以是奇异的（例如尖角），这与文献中关于非散射介质可以在角域存在的结论一致。本文的工作实际上解释了为什么在非散射背景下，介质 $\Phi^*[I]$ 的高阶光滑性“强迫”边界必须光滑，除非发生特定的几何退化。

4. 理论贡献的深度： 本文不仅是一个存在性/正则性定理，更揭示了 $\Phi^*[I]$ 与边界几何之间深刻的内在联系。它表明，$\Phi^*[I]$ 的光滑性不仅仅是系数函数的属性，更是边界几何属性的反映。这对于理解电磁成像中的不可识别性（non-uniqueness）和非散射介质的构造具有重要意义。如果边界不够光滑，那么对应的 $\Phi^*[I]$ 必然会在边界附近失去高阶光滑性，或者必须满足特定的退化条件。这一结论完善了非散射反演理论的数学基础，为后续研究提供了严格的先验约束。

了解领域论文（摘要翻译）

1. 奇异摄动Korteweg–de Vries方程的相移纳米波解

Phase-Shifted Nanopteron Solutions to a Singularly Perturbed Korteweg–de Vries Equation

作者: Timothy E. Faver

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.19097v1

链接: http://arxiv.org/abs/2608.19097v1

摘要翻译

中文标题: 奇异摄动Korteweg–de Vries方程的相移纳米波解

摘要翻译: 我们构造了微分方程 $\varepsilon^2u^{(4)}+u''-u+u^2 =0$ 的纳米波（nanopteron）解，该方程是 Korteweg–de Vries (KdV) 方程行波问题 $u''-u+u^2=0$ 的奇异摄动。这些纳米波解是 sech$^2$ 型 KdV 行波轮廓、一个小局部误差和一个周期性"涟漪"的叠加。该涟漪由其振幅和相移参数化。在构造纳米波的过程中，我们先选定这两个参数中的一个，然后求解另一个，并对振幅与相移的求解过程如何不同提供详细评论。我们的工作结合关于格点微分方程中纳米波的现代发展，重新审视了此前针对该奇异摄动 KdV 方程的若干纳米波解构造。在这些当代方法的背景下，奇异摄动 KdV 方程的简单结构使其能够作为更复杂问题中纳米波构造的一个易入门模板。

2. 自发聚集模型的梯度流结构

Gradient Flow Structure of the Spontaneous Aggregation Model

作者: Jan Haskovec

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.19005v1

链接: http://arxiv.org/abs/2608.19005v1

摘要翻译

中文标题: 自发聚集模型的梯度流结构

摘要翻译: 我们识别出在自发粒子聚集模型中出现的 Fokker–Planck 方程先前未被注意到的梯度流结构。对于指数响应函数，该方程是具有非局部迁移率的吸引型 McKean–Vlasov 自由能的广义 Wasserstein 梯度流。我们表明，在一个自然的局部熵和对称相互作用能类中，这种结构本质上唯一地确定了指数响应。我们讨论了其对稳态、全局极小值、一维单簇结构的影响，并形式化地利用 Wasserstein–Łojasiewicz 不等式概述了向平衡态的收敛。

3. 二次 Aharonov–Bohm NLS 的长程渐近行为

Long Range Asymptotics for the Quadratic Aharonov–Bohm NLS

作者: Piero D’Ancona, Tohru Ozawa

类别: math.AP, math-ph

发布日期: 2026-08-19

arXiv ID: 2608.18960v1

链接: http://arxiv.org/abs/2608.18960v1

摘要翻译

中文标题: 二次 Aharonov–Bohm NLS 的长程渐近行为

摘要翻译: 我们研究了方程 $i\partial_tu=H_αu+λ|u|u$ 在 $\mathbb R^2$ 上解的长程行为，其中 $H_α$ 是带有单个极点的 Aharonov–Bohm Hamiltonian 的 Friedrichs 实现。长程 ansatz 的对数相位可能会将初值轮廓移出 $H_α$ 的定义域。我们通过边界迹在 0 处消失且阶数 $\le \frac 12$ 来刻画那些保持在算子定义域内的轮廓；在半通量 $α=1/2$ 时，没有非零的迹幸存。然而，对于 $H_α$ 定义域中的每一个具有小 $L^\infty$ 振幅的轮廓，都确定了一个具有修正终态的唯一全局解，其余项速率为 $t^{-b}$，对所有 $0<b<1/2+ν_α$ 成立，其中 $ν_α=\min\{α,1-α\}$。对于满足消失迹条件的轮廓，该速率改进为所有 $0<b<1$。这一结果在最优意义下是 sharp 的：如果 $α\neq \frac 12$，我们可以构造误差量为 $t^{-1/2-ν_α}\log t$ 的轮廓，从而排除所有更快的速率。上界来源于一个不属于 $L^2$ 的残差所对应的推迟 Strichartz 估计；下界则是通过 Hankel 变换进行的显式计算。对于更光滑的轮廓，我们还计算了第一阶修正项，这给出了具有 $1<b<2$ 的余项速率。

4. 具有阈值奇点的二维四阶薛定谔算子的时间衰减估计

Time-Decay Estimates for Two-Dimensional Fourth-Order Schrödinger Operators with Threshold Singularities

作者: Zijun Wan, Xiaohua Yao

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18916v1

链接: http://arxiv.org/abs/2608.18916v1

摘要翻译

中文标题: 具有阈值奇点的二维四阶薛定谔算子的时间衰减估计

摘要翻译: 我们建立了具有实值衰减势$V$的二维四阶薛定谔算子$H=Δ^2+V$的时间衰减估计，涵盖了所有可能的零能量阈值障碍。当零是正则点或第一类共振时，我们证明了 \[ \left| H^{\fracα{4}}e^{-itH}P_{\mathrm{ac}}(H) \right|{L^1\to L^\infty} \lesssim |t|^{-\frac{2+α}{4}}, \qquad -2<α\leq2, \]这与自由算子在$α$的整个范围内的最佳衰减速率一致。对于第二类共振，对于每个$-2<α\leq2$，衰减速率为$|t|^{-(2+α)/4}(\log(2+|t|))^2$，仅存在对数损失。对于更强的阈值奇点，我们表明长时间行为由$d$波共振的存在所支配。如果零是第三类共振，或者是伴随有$d$波共振的特征值，我们获得了$α=0$时的最佳衰减$(\log|t|)^{-1}$以及$0<α\leq2$时的$|t|^{-α/4}(\log|t|)^{-2}$。如果零是无$d$波共振的特征值，则在$-2<α\leq2$范围内恢复了第二类估计。此外，在正则和第一类共振情形下，对于每个$2<α\leq2$和$s>0$，我们获得了改进的对数加权估计： \[ \left| ω^{-s} H^{\fracα{4}}e^{-itH}P{\mathrm{ac}}(H)ω^{-s} \right|_{L^1\to L^\infty} \lesssim \frac{1}{|t|^{\frac{2+α}{4}}(\log|t|)^s}, \qquad |t|\geq2, \]其中$ω(x)=\log(2+|x|)$。相比之下，对于自由算子$Δ^2$，零是第二类共振，且自由演化不具有此类对数增益。因此，在正则和第一类情形中，势改变了自由算子的零能量谱结构，而这种改变伴随着改进的加权衰减。

5. 一类非线性热传导模型的最优控制

Optimal control of a class of nonlinear heat conduction models

作者: Denilson Menezes, Juan Límaco

类别: math.OC, math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18896v1

链接: http://arxiv.org/abs/2608.18896v1

摘要翻译

中文标题: 一类非线性热传导模型的最优控制

摘要翻译: 我们研究了一类描述金属热处理过程中非线性热传导的拟线性抛物方程的最优控制问题。在一般假设下，我们证明了最优解的存在性，并在一个具有实际相关性的框架中刻画了相应的最优控制-状态对。主要困难在于二次梯度项，我们通过能量估计和控制到状态映射的可微分性质来处理该项。

6. 受迫二维Navier-Stokes方程的有限时间爆破

Finite-time blow-up for the forced 2D Navier-Stokes equations

作者: Filippo Gazzola

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18802v1

链接: http://arxiv.org/abs/2608.18802v1

摘要翻译

中文标题: 受迫二维Navier-Stokes方程的有限时间爆破

摘要翻译: 我们提供了显式示例，表明文献中关于力的现有假设对于二维Navier-Stokes方程强解的所述性质是尖锐的。如果这些假设不满足，则可能出现不同形式的有限时间爆破。

7. Dirichlet-to-Neumann 算子的比较不等式

Comparison inequalities for Dirichlet-to-Neumann maps

作者: Denis S. Grebenkov, Michael Levitin, Karl-Mikael Perfekt, Iosif Polterovich

类别: math.SP, math.AP, math.FA

发布日期: 2026-08-19

arXiv ID: 2608.18678v1

链接: http://arxiv.org/abs/2608.18678v1

摘要翻译

中文标题: Dirichlet-to-Neumann 算子的比较不等式

摘要翻译: 我们证明了对应于不同非正亥姆霍兹参数的 Dirichlet-to-Neumann 算子的比较不等式。对于凸区域，我们的界限是尖锐的（sharp），所得出的特征值不等式部分证实了早期的一个猜想，但我们表明该猜想在完全一般情形下不成立。此外，我们还获得了适用于任意足够正则区域的依赖于几何的版本，并将其推广至带边界的紧致黎曼流形。我们还讨论了度量图上的类似问题。

8. 从磁化库仑量子动力学到磁化流体

From magnetized Coulombic quantum dynamics to magnetized fluids

作者: Immanuel Ben Porat

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18476v1

链接: http://arxiv.org/abs/2608.18476v1

摘要翻译

中文标题: 从磁化库仑量子动力学到磁化流体

摘要翻译: 我们扩展了文献 [17] 中发展的量子调制能量，旨在分别从磁化 Schrödinger-Poisson 方程和 von-Neumann 方程导出磁化无压 Euler-Poisson 方程的半经典极限与平均场半经典极限。我们还讨论了底层单速 PDE 的局部适定性。在这两种极限下，磁场均为外场且可能在空间上非均匀。我们的结果属于磁化量子动力学的半经典极限与量子平均场极限这一更广泛的范畴。

9. 具有有限秩代数几何背景的积分聚焦与反聚焦 Wadati-Konno-Ichikawa 方程的长时间渐近行为

Long-time asymptotics of the integrable defocusing Wadati-Konno-Ichikawa equation with a finite-genus algebro-geometric background

作者: Taohua Luo, Zhenya Yan, Guoqiang Zhang

类别: nlin.SI, math-ph, math.AP, nlin.PS, physics.optics

发布日期: 2026-08-19

arXiv ID: 2608.18630v1

链接: http://arxiv.org/abs/2608.18630v1

摘要翻译

中文标题: 具有有限秩代数几何背景的积分聚焦与反聚焦 Wadati-Konno-Ichikawa 方程的长时间渐近行为

摘要翻译: 我们研究了具有饱和非线性的 Wadati-Konno-Ichikawa (WKI) 方程的有限秩代数几何解，及其短程扰动的长时间渐近行为。首先，针对聚焦与反聚焦约化情形，我们将有限秩 Baker-Akhiezer 函数表述为复谱平面上可显式求解的矩阵 Riemann-Hilbert (RH) 问题，并获得 theta 函数表示以及 WKI 场和倒数坐标的重构公式。随后，我们考虑在有限秩代数几何背景下的反聚焦 WKI 方程的 Cauchy 问题。我们构造了散射数据和 RH 问题，并执行 Deift-Zhou 非线性最速下降分析。时空平面被划分为两个过渡区域：Zakharov-Manakov (ZM) 区域和快衰减区域。主导项为一个相位移动的有限秩 WKI 解。过渡修正由 Painlevé-XXXIV 模型控制，而 ZM 辐射则由抛物柱面函数描述。同时获得了倒数坐标的渐近性。

10. 度量空间中Eikonal方程Monge解的稳定性与Perron方法

Stability and Perron’s method for Monge solutions of eikonal equations in metric spaces

作者: Qing Liu, Made Benny Prasetya Wiranata

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18518v1

链接: http://arxiv.org/abs/2608.18518v1

摘要翻译

中文标题: 度量空间中Eikonal方程Monge解的稳定性与Perron方法

摘要翻译: 本文的主要结果涉及度量空间中Eikonal方程Monge解的稳定性。尽管标准定义的Monge解中使用的子斜率（subslope）在一致收敛下通常不具备稳定性，但我们通过引入基于球面子斜率条件的等价非局部刻画克服了这一困难。该方法补充了现有的基于度量粘性解的方法，使我们能够在不使用测试函数的情况下，直接在一致收敛以及定义域和度量的Gromov-Hausdorff扰动下建立稳定性。该方法通过光学度量扩展到非齐次方程，并推广到更广泛的Hamilton-Jacobi方程类。我们的稳定性论证进一步使得我们在一般完备长度空间上发展Monge解的Perron方法成为可能。这给出了一个直接获得、不依赖于粘性技术的度量空间中新的Perron型存在性结果。

11. 无质量共振的四维质量临界二次非线性薛定谔系统的有限时间爆破

Finite-time blow-up for the four-dimensional mass-critical quadratic nonlinear Schrödinger system without mass resonance

作者: Ngoc Uyen Cong Nguyen, Van Duong Dinh

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18453v1

链接: http://arxiv.org/abs/2608.18453v1

摘要翻译

中文标题: 无质量共振的四维质量临界二次非线性薛定谔系统的有限时间爆破

摘要翻译: 我们研究聚焦型二次非线性薛定谔系统：

$$ \begin{cases} i\partial_t u+Δu=-2v\overline{u}, \\ i\partial_t v+κΔv=-u^2, \end{cases} \qquad (t,x)\in I\times\mathbb R^4, $$

其中 $κ>0$。在非质量共振情形 $κ\neq \frac12$ 下，Inui–Kishimoto–Nishimura 与 Dinh–Forcella 的先前工作表明，具有负能量的径向解要么在有限时间内爆破，要么全局存在但其 $H^1$-范数无界增长。 在本文中，我们证明了每一个具有负能量的径向 $H^1\times H^1$ 解在正向和反向时间上均在有限时间内爆破。无需任何有限方差假设。 主要工具是基于有界指数权重的局域化位力论证：

$$ \nablaφ_R(x)=2x e^{-|x|^2/R^2}. $$

一个径向加权插值估计使得我们能够由相应的加权动能项控制非线性误差项，其误差为仅依赖于守恒质量的 $O(R^{-2})$。此外，局域化位力量本身可以直接由相同的加权动能缺陷进行界定。结合这些估计导出了一个超线性 Riccati 型微分不等式，该不等式无法对所有时间成立，从而迫使解在有限时间内爆破。

12. 软势Landau方程行波Maxwellian的线性稳定性

Linear stability of traveling Maxwellians for Landau equation with very soft potentials

作者: Sanchit Chaturvedi, Jonathan Luk

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18407v1

链接: http://arxiv.org/abs/2608.18407v1

摘要翻译

中文标题: 软势Landau方程行波Maxwellian的线性稳定性

摘要翻译: 考虑在库仑势下围绕行波Maxwellian线性化的Landau方程。我们证明了线性解的微观部分以多项式逆率衰减，而宏观部分保持有界但通常不衰减。特别地，其长时间动力学不同于以下两种情形的线性化：(1) 针对中等软势或硬势的行波Maxwellian；(2) 针对任意势的全局（非行波）Maxwellian。

13. $\mathbb{R}^N$上具有多孔介质扩散、线性产生和逻辑源的趋化模型弱解的整体存在性

Global existence of weak solutions to chemotaxis models with porous medium diffusion, linear production, and logistic source on $\mathbb{R}^N$

作者: Zulaihat Hassan

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18416v1

链接: http://arxiv.org/abs/2608.18416v1

摘要翻译

中文标题: $\mathbb{R}^N$上具有多孔介质扩散、线性产生和逻辑源的趋化模型弱解的整体存在性

摘要翻译: 本文研究了趋化系统 egin{equation*} \begin{cases} u_t = Δu^m - χ∇ · (u ∇v) + u(a - b u), & \text{in } (0,∞)×\mathbb{R}^N, \ τv_t = Δv - λv + μu, & \text{in } (0,∞)×\mathbb{R}^N, \end{cases} \end{equation*} 的弱解的整体适定性、有界性和唯一性，其中 $m>1$，$τ∈\{0,1\}$，$λ,μ,a,b>0$，且 $χ∈\mathbb{R}$。对于任意 $m>1$，我们建立了非必然可积初值下弱解的存在性，尽管此类解一般未必有界。随后，在抛物-抛物情形 $(τ=1)$ 中，当 $m>\frac{2N}{N+2}$ 时，我们证明了整体有界弱解的存在性；此外，当 $1<m≤\frac{2N}{N+2}$ 且逻辑阻尼系数 $b$ 足够大时，同样证明了整体有界弱解的存在性。在抛物-椭圆情形 $(τ=0)$ 中，当 $m>2-\frac{2}{N}$ 时，我们证明了整体有界弱解的存在性；此外，当 $1<m≤2-\frac{2}{N}$ 且 $b$ 足够大时，也证明了整体有界弱解的存在性。最后，对于 $1<m≤3$，我们证明了直至初始时刻 Hölder 连续的弱解的唯一性。

14. 具有非局部Kelvin–Voigt阻尼的五次波动方程的吸引子与奇异极限

Attractors and Singular Limits for a Quintic Wave Equation with Nonlocal Kelvin–Voigt Damping

作者: Yue Sun, Marcelo M. Cavalcanti, Vando Narciso

类别: math.AP, math.DS

发布日期: 2026-08-18

arXiv ID: 2608.18266v1

链接: http://arxiv.org/abs/2608.18266v1

摘要翻译

中文标题: 具有非局部Kelvin–Voigt阻尼的五次波动方程的吸引子与奇异极限

摘要翻译: 本文考虑定义在有界区域$Ω\subset\mathbb{R}^3$上的能量临界五次波动方程，其非线性且非局部的Kelvin–Voigt阻尼项形式为$-\|\nabla u_t\|_{L^2(Ω)}^αΔu_t$，其中$α\in\mathbb{R}_+=[0,\infty)$。在关于五次源项的适当假设下，我们建立了该问题的适定性，并在自然能量空间$\mathcal H=H_0^1(Ω)\times L^2(Ω)$中研究了其长时间动力学行为。对于任意$α\in\mathbb{R}_+$，我们证明了相应的动力系统$(\mathcal H,S^α(t))$是梯度耗散的，并导出一个稳定化估计，该估计给出了渐近光滑性，进而保证了紧致全局吸引子$\mathcal A_α$的存在性。同一估计提供了$\mathcal A_α$的Kolmogorov $\varepsilon$-熵的上界；在极限情形$α=0$下，该估计退化为拟稳定性不等式，由此表明$\mathcal A_0$具有有限分形维数。此外，我们证明了族$\{\mathcal A_α\}_{α\in\mathbb{R}_+}$在高正则性空间$\mathcal H_1=(H^2(Ω)\cap H_0^1(Ω))\times H_0^1(Ω)$中一致有界。最后，我们确立了$\{\mathcal A_α\}_{α\in\mathbb{R}_+}$在$α=0$处的上半连续性，表明具有非线性且非局部Kelvin–Voigt阻尼的吸引子收敛于具有经典线性Kelvin–Voigt阻尼的极限问题的全局吸引子。

不相关论文（基本信息）

1. Hessian 商方程的内估计

Interior estimates for the Hessian quotient equations

作者: Weisong Dong, Ruijia Zhang

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.19087v1

链接: http://arxiv.org/abs/2608.19087v1

摘要: In this paper, we establish interior $C^2$ estimates for admissible semiconvex solutions to the general Hessian quotient equation $ \frac{σ_k}{σ_l}(D^2u)=f(x,u), $ for the cases $l=k-1$ and $l=k-2$, where $f$ is a positive $C^{1,1}$ function. Such estimates are known to fail in general for $k-l\geq 3$, even for convex solutions, as shown by counterexamples due to Lu. The main ingredient is a quantitative concavity inequality for the Hessian quotient operator under the semiconvex condition. Our result provides a unified argument to such general Hessian quotient for $2\leq k\leq n-1$ in arbitrary dimensions.

2. 正交不变域中 Robin 函数的 Hessian 的特征分解

Eigendecomposition of the Hessian of the Robin function in orthogonally invariant domains

作者: Alejandro Ortega

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.19169v1

链接: http://arxiv.org/abs/2608.19169v1

摘要: In this work we analyze the eigendecomposition of the Hessian matrix of the Robin function $\mathcal{R}(x)$ for the spectral fractional Laplacian in orthogonally invariant domains. We prove that, if $Ω$ a smooth bounded convex domain invariant under the action of an orthogonal transformation $\mathcal{O}$ then, for $\overline{t}\in\{a\inΩ:\mathcal{O}(a)=a\}$, the gradient vector $\nabla\mathcal{R}(\overline{t})$ is an eigenvector of the Jacobian $D\mathcal{O}$ associated to the eigenvalue $1$. Moreover, if $Ω$ is a domain invariant under the reflection about a hyperplane $π_{v}=\{x\in\mathbb{R}^N:x\cdot v=0\}$, there exists $η>0$ such that $\mathbb{H}(\overline{t})v=ηv$ where $\mathbb{H}$ denotes the Hessian matrix of $\mathcal{R}(x)$. Consequently, if $Ω$ is invariant under the reflection about the hyperplanes $π_{v_i}$ for a linearly independent set $\{v_1,\ldots,v_N\}$, then the origin is a non degenerate critical point of $\mathcal{R}(x)$. A short proof of the Brezis-Peletier-like formulas for $\nabla \mathcal{R}$ is also provided which allows us to prove the former results for any $0<s<1$.

3. 分析极小模型纲领中的 Kähler-Ricci 切线流

Kähler-Ricci Tangent Flows in the Analytic Minimal Model Program

作者: Longteng Chen, Max Hallgren, Lucas Lavoyer

类别: math.DG, math.AG, math.AP, math.CV

发布日期: 2026-08-19

arXiv ID: 2608.19152v1

链接: http://arxiv.org/abs/2608.19152v1

摘要: We describe certain finite-time singularities of the Kähler-Ricci flow arising in the analytic minimal model program. Assuming that convergence to an asymptotically conical Kähler-Ricci shrinker is realized by holomorphic maps, we prove that, in a fixed holomorphic gauge, the nearby flow is modeled on the shrinker at the level of Kähler potentials. Consequently, every noncollapsed Kähler-Ricci flow through singularities in complex dimension two is modeled on a shrinker-cone-expander transition, confirming a strong form of Song’s conjectural picture. We also show analogous results in higher dimensions under the Calabi ansatz, and improve known results in the compact shrinker case. These give the first compact Ricci flows through conical singularities whose small-scale behavior is fully described.

4. Maxwell 透射问题中表面局域化特征函数的 Weyl 型规律

A Weyl-type law for surface-localized eigenfunctions of the Maxwell transmission problem

作者: Yan Jiang, Hongyu Liu, Kai Zhang, Haoran Zheng

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18872v1

链接: http://arxiv.org/abs/2608.18872v1

摘要: This work studies surface-localized transmission eigenfunctions for time-harmonic electromagnetic scattering. Prior results on this phenomenon are mostly qualitative, proving existence without quantifying distribution. Here we provide a quantitative analysis for the Maxwell transmission eigenvalue problem, covering both TE and TM modes. We first establish a Weyl asymptotic law for the full counting function, with cubic growth with respect to the radius $R$. We then prove Weyl-type upper and lower bounds for the counting functions restricted to surface-localized modes, and demonstrate that they share the same growth order three. To our knowledge, this is the first quantitative result of its kind, demonstrating that surface-localized eigenfunctions form a non-negligible fraction of the high-frequency spectrum.

5. 从局部 Neumann-to-Dirichlet 映射稳定确定复各向异性导纳及其导数

Stable boundary determination of a complex anisotropic admittivity and its derivatives from a local Neumann-to-Dirichlet map

作者: Jessica Crosse, Romina Gaburro

类别: math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18674v1

链接: http://arxiv.org/abs/2608.18674v1

摘要: We study the classical anisotropic Calderón problem associated to the elliptic equation $\text{div}(σ\nabla u)=0$, where the complex admittivity $σ$ is of the form $σ=A(\cdot,a(\cdot))$ in a domain $Ω\subset\mathbb{R}^n$, $n\ge3$. We establish boundary stability estimates for $σ$ and its derivatives of arbitrary order from a local Neumann-to-Dirichlet map. Our results extend those of Comm. Partial Differential Equations, 34 (2009) from the real-valued conductivity setting to complex anisotropic admittivities.

6. 有界域上二阶线性椭圆 Dirichlet 边值问题的解的对称性与区域缩减有限元方法

Symmetry of Solutions and Domain-Reduction Finite Element Method for Second-Order Linear Elliptic Dirichlet Boundary Value Problems on Bounded Domains

作者: Xianlong Pan, Wei Jiang

类别: math.NA, math.AP

发布日期: 2026-08-19

arXiv ID: 2608.18549v1

链接: http://arxiv.org/abs/2608.18549v1

摘要: Combining classical group theory and partial differential equation theory, this paper investigates the symmetry group $\operatorname{Sym}(u)$ of the unique solution $u$ to the second-order linear elliptic boundary value problem on an $n$-dimensional bounded domain $Ω$ $-\sum_{i,j=1}^{n} a_{ij}(x)u_{x_ix_j} + \sum_{i=1}^{n} b_i(x)u_{x_i} + c(x)u = f(x), x\in Ω$, $u(x) = h(x), x\in \partial Ω$, The following symmetry groups are defined and characterized respectively: the symmetry group $\operatorname{Sym}(A)$ of the second-order coefficient matrix function $A(x)=(a_{ij}(x))_{n\times n}$; the symmetry group $\operatorname{Sym}(b)$ of the first-order coefficient column vector function $b(x)=(b_{1}(x),b_{2}(x),\cdots,b_{n}(x))^{T}$; the symmetry group $\operatorname{Sym}(c)$ of the zero-order coefficient function $c(x)$; the symmetry group $\operatorname{Sym}(f)$ of the internal source function $f(x)$; and the symmetry group $\operatorname{Sym}(h)$ of the boundary source function $h(x)$. This paper rigorously proves that the common symmetry group $\operatorname{Sym}(A)\cap\operatorname{Sym}(b) \cap\operatorname{Sym}(c) \cap\operatorname{Sym}(f) \cap\operatorname{Sym}(h)$ is a subgroup of $\operatorname{Sym}(u)$. In addition, if the common symmetry group contains several mirror symmetry elements, the original second-order linear elliptic boundary value problem on the entire domain $Ω$ can be reduced to the corresponding boundary value problem on a certain subdomain. It is strictly proven in this paper that the new boundary condition imposed on the boundary of the subdomain is the homogeneous generalized Neumann boundary condition. The linear finite element method is used to numerically solve the second-order linear elliptic boundary value problem on the subdomain, thereby achieving domain reduction and significantly reducing the computational cost.

7. 有限博弈方最优停止博弈：随机化、$α$-势性与学习

Finite-player Optimal Stopping Games: Randomization, $α$-potentiality, and Learning

作者: Xin Guo, Mehdi Talbi, Qinxin Yan

类别: math.OC, math.AP

发布日期: 2026-08-18

arXiv ID: 2608.18355v1

链接: http://arxiv.org/abs/2608.18355v1

摘要: Finite-player nonzero-sum optimal stopping games typically lead to coupled equilibrium systems whose complexity grows rapidly with the number of players. We introduce an independently randomized formulation in which each stopping rule is represented by an adapted, nondecreasing cumulative stopping process. The canonical embedding preserves pure-profile payoffs, and a pure profile is a Nash equilibrium of the original game if and only if its embedding is a Nash equilibrium of the randomized game. We adopt the $α$-potential approach to construct an $α_N$-potential function, with the error $α_N=O(N^{-1})$ under weak-interaction. We also identify an exact-potential subclass with a closed-form threshold equilibrium. For local stopped-status interactions, randomized payoffs admit a local stopped-mass representation, and potential maximization can be formulated as a multidimensional singular-control problem with local gradient constraints and a nonlocal condition for finite jumps. Under suitable regularity assumptions, we study the associated Hamilton-Jacobi-Bellman quasi-variational inequality and its regularity properties. For unknown model coefficients, we propose a bounded-intensity Potential-CT-DDPG learning algorithm. Numerical experiments closely match the analytical benchmark and yield estimated best-response improvements consistent with $N^{-1}$ scaling.

8. 关于完全非线性一致椭圆方程的奇点集

On Singular Sets of Fully Nonlinear Uniformly Elliptic Equations

作者: BaoZhi Chu

类别: math.AP, math.DG

发布日期: 2026-08-18

arXiv ID: 2608.18217v1

链接: http://arxiv.org/abs/2608.18217v1

摘要: For a continuous viscosity solution $u$ of a fully nonlinear uniformly elliptic equation $F(D^2 u)=0$, the work of Nadirashvili–Tkachev–Vlăduţ shows that $u$ need not be $C^2$ in dimension $n\ge 5$. It is therefore natural to ask how large the singular set $\operatorname{Sing}(u)$, consisting of points having no $C^2$ neighborhood of $u$, can be. Under an additional $C^1$ assumption on $F$, Armstrong–Silvestre–Smart proved that $\operatorname{Sing}(u)$ has Hausdorff measure at most $n-ε$ for some $ε>0$. In this paper, we show that for general uniformly elliptic $F$, modulo a countable set, $\operatorname{Sing}(u)$ can be any nonempty compact nowhere-dense set in $\mathbb{R}^5$. In particular, without further assumptions on $F$, the Hausdorff dimension of $\operatorname{Sing}(u)$ can be any number between $0$ and $n$ in dimension $n\geq 5$; moreover, the singular set can even have positive Lebesgue measure.