生成时间: 2026年05月04日 06:03:14
重点关注论文数量: 5
了解领域论文数量: 7
不相关论文数量: 4
重点关注论文（完整分析）
1. 光滑且稳定的欧拉坍缩解
Smooth and stable Euler implosions
作者: Jiajie Chen, Steve Shkoller, Vlad Vicol
类别: math.AP
发布日期: 2026-05-01
arXiv ID: 2605.00808v1
摘要翻译: 我们为多维可压欧拉方程构造了一类新的自相似坍缩型剖面（implosion profiles）。这些剖面是光滑的、真正非等熵的、径向/球对称的，且具有显式（闭式）的相似指数。我们证明：对应于基态坍缩剖面的精确欧拉解，在全非线性可压欧拉方程意义下，关于径向对称扰动是稳定的，前提是初始数据满足一个一维相容性条件。对于单原子或双原子气体对应的基态坍缩剖面之欧拉解的扰动（不作任何对称性假设），我们完全刻画了导致非线性稳定的初值集合。
Comment: 153 pages, 6 figures
链接: http://arxiv.org/abs/2605.00808v1
1. 研究对象和背景
本研究聚焦多维可压缩欧拉方程的内爆解存在性与稳定性问题。研究对象为全（非等熵）可压缩欧拉方程系统：
其中压力 $p$ 满足 $p = (\gamma - 1)\left(E - \frac{1}{2}\rho|u|^2\right)$，$\gamma > 1$ 为绝热指数（物理上典型值为 $\gamma = 5/3$（三维单原子气体）和 $\gamma = 7/5$（三维双原子气体））。在经典解范畴内（密度严格正且流动变量 $C^1$ 光滑），系统等价于主要流动变量的演化方程：
作者引入重缩放声速 $\sigma$、流体速度 $u$ 和伪熵平方根 $b$ 作为基本变量（$\alpha = \frac{\gamma - 1}{2} > 0$）：
其中 $b = \sqrt{\gamma p / \rho^\gamma} = e^{s/2}$（$s$ 为比熵），$\sigma = \frac{1}{\alpha} \rho^\alpha b$。内爆解的数学定义（Definition 1.3）要求：(i) 在坍缩时空点 $(x_*, t_*)$ 处至少一个主要热力学变量（密度或压力）无界，即 $\forall \varepsilon \in (0,1]$ 满足
(ii) 流动具有向内聚焦特性：$\exists r(t) \to 0$ 使得 $(u(x,t) - u(x_*,t)) \cdot (x - x_*) \leq 0$ 对 $|x - x_*| \leq r(t)$ 成立；(iii) 远离 $(x_*, t_*)$ 的时空点上主要流动变量有界。
已有内爆解主要分两类：(1) 冲击波内爆（Guderley解）：含光滑余维1冲击前沿，相似指数需数值求解声学点非线性特征值问题（Lazarus[45]），不能从全局 $C^2$ 光滑初始数据产生，且假设零压力静止核心，与Navier-Stokes方程零粘性极限不一致；(2) 光滑等熵内爆（Merle等[48,7,60]）：无真空区域、$C^\infty$ 光滑、无冲击，但相似指数隐式确定且数值分析[3]表明在径向对称下可能不稳定。两类解的局限性在于：相似指数非显式、稳定性未严格证明、物理适用性受限。
本文试图填补的关键空白包括：(1) 构造具有显式闭式相似指数的光滑、真正非等熵自相似内爆剖面；(2) 严格证明基态剖面在径向及非径向扰动下的非线性稳定性；(3) 建立适用于数值模拟的解析基准解，解决冲击波内爆的数值实现困难和光滑等熵内爆的稳定性缺失问题。特别地，作者关注压力仅在原点为零（动态保持的单点相容性条件）的严格正密度解，避免Guderley解中零压力区域的数值困难，并突破等熵假设以实现更丰富的爆破行为（密度必爆破但压力未必）。
2. 主要定理或主要结果
在空间维度 $d \in \{1,2,3\}$ 和绝热指数 $\gamma \in (1, 2d+1]$ 下，作者证明以下核心结论：
自相似剖面构造（Section 2）：构造由整数 $N \geq 1$ 索引的光滑自相似内爆剖面序列。相似指数 $\bar{c}_r$、$\bar{c}_u$、$\bar{c}_b$ 具有显式闭式表达式（Definition 2.5）：
$$ \bar{c}_r = \frac{d + 2N - 1}{d + 2N + 1}, \quad \bar{c}_u = \bar{c}_r - 1, \quad \bar{c}_b = \frac{2N}{d + 2N + 1}. $$当 $N=1$ 时称为基态剖面（相似指数最大）。自相似变量满足径向对称欧拉方程（1.7）的稳态解：
$$ \begin{aligned} (\bar{c}_r R - U) \partial_R \Sigma + \alpha \Sigma \left( \partial_R U + \frac{d-1}{R} U \right) &= 0, \\ (\bar{c}_r R - U) \partial_R U + \alpha \Sigma \partial_R \Sigma - \frac{\alpha}{\gamma} \Sigma^2 \frac{\partial_R B}{B} &= 0, \\ (\bar{c}_r R - U) \partial_R B &= 0, \end{aligned} $$其中 $U(R)$、$\Sigma(R)$、$B(R)$ 为实解析函数，密度 $\rho(R) = (\alpha \Sigma / B)^{1/\alpha} > 0$ 处处严格正，压力 $p(R) = \alpha^2 \rho \Sigma^2 / \gamma$ 仅在 $R=0$ 处为零。当 $\gamma(\bar{c}_r - 1) \geq \bar{c}_b$ 时压力有界（如 $d=2$、$\gamma=2$、$N=1$），否则压力爆破（如 $d=3$、$\gamma=5/3$、$N=1$）；速度幅值爆破条件为 $\bar{c}_r < 1$（Remark 1.6）。
径向对称稳定性（Theorem 3.15）：对任意 $d$、$\gamma$，基态剖面（$N=1$）在径向对称扰动下渐近稳定。设初始扰动满足 $\|\delta \sigma(\cdot, -1)\|_{C^M} + \|\delta u(\cdot, -1)\|_{C^M} + \|\delta b(\cdot, -1)\|_{C^M} \leq \varepsilon$（$M$ 足够大）且压力在原点为零的相容性条件，存在 $\varepsilon_0 > 0$ 使得当 $\varepsilon < \varepsilon_0$ 时，扰动解收敛到精确内爆解，误差满足：
$$ \sup_{\tau \geq 0} e^{\kappa \tau} \left( \|\delta U(\cdot, \tau)\|_{H^k} + \|\delta \Sigma(\cdot, \tau)\|_{H^k} + \|\delta B(\cdot, \tau)\|_{H^k} \right) \leq C \varepsilon, $$其中 $\kappa > 0$ 为衰减速率，$k \geq 3$。稳定性模去一维相容性条件（初始数据需满足 $\partial_R^{2N} B(0,-1) = \text{const}$）。
非径向稳定性（Theorems 4.4 & 5.1）：对单原子气体（$\gamma = \frac{d+2}{d}$）和双原子气体（$\gamma = \frac{2d+1}{2d-1}$），基态剖面在非径向扰动下的稳定性完全刻画。设扰动初始数据满足：
$$ \sum_{|\beta| \leq M} \sup_{x \in \mathbb{R}^d} \left| \partial^\beta_x (\rho_0 - \rho_{\text{exact}})(x) \right| (1 + |x|^{2N}) \leq \varepsilon, $$其中 $\rho_{\text{exact}}$ 为精确解，$M$ 足够大。不稳定+中心流形维度为 $\binom{d+1}{2}$（对应平移和旋转对称性），稳定性条件为：
$$ \int_{\mathbb{R}^d} x^\alpha (\rho_0 - \rho_{\text{exact}})(x) dx = 0, \quad |\alpha| \leq 2N, \quad \alpha \neq (2N,0,\dots,0), \dots, (0,\dots,0,2N), $$即初始扰动需满足有限个（$O(N^d)$）显式相容性条件。满足条件时，扰动解满足：
$$ \sup_{t \in [-1,t_*)} (t_* - t)^{-\bar{c}_r} \| \rho(\cdot,t) - \rho_{\text{exact}}(\cdot,t) \|_{L^\infty} \leq C \varepsilon. $$爆破行为精确刻画（Section 1.6.1）：对基态剖面（$N=1$），密度在 $r \to 0$ 时满足 $\rho(r,t) \sim (t_* - t)^{-(d+1)/(d+3)} r^{-(2d+2)/(d+3)}$，压力行为分两种情形：(i) $d=3$、$\gamma=5/3$ 时 $p(r,t) \sim (t_* - t)^{-10/9} r^{-8/3}$（爆破）；(ii) $d=2$、$\gamma=2$ 时 $p(r,t) \to p_\infty(r)$ 有界（Figure 4）。速度幅值 $|u|$ 在 $\bar{c}_r \geq 1$ 时有界（如 $d=3$、$\gamma=5/3$、$N=1$ 时 $\bar{c}_r = 3/5 < 1$ 故爆破；$N$ 较大时 $\bar{c}_r \geq 1$ 则有界）。
3. 研究方法、关键技术和核心工具
证明中的核心障碍与创新方法如下：
1. 自相似剖面的显式构造（Section 2） 传统内爆解（Guderley或Merle等）需数值求解声学点非线性特征值问题确定相似指数，作者创新性地将问题转化为求解常微分方程组的显式解。关键步骤包括：(i) 在自相似坐标 $(R,U,\Sigma,B)$ 下将欧拉方程化为：
其中分母 $D(R) = (\bar{c}_r R - U)^2 - \alpha \Sigma^2$ 在声学点 $R=1$ 处为零（$D(1)=0$），要求分子同时为零以保证解光滑；(ii) 通过设定 $\bar{c}_r = \frac{d+2N-1}{d+2N+1}$ 和 $\bar{c}_b = \frac{2N}{d+2N+1}$，使声学点条件 $D(1)=0$ 与分子为零条件相容，导出相似指数的闭式解；(iii) 证明解在 $R=0$ 和 $R \to \infty$ 处的渐近行为：$R \to 0$ 时 $\Sigma(R) \sim R^{2N}$，$U(R) \sim -R^{2N+1}$，$B(R) \sim \text{const}$；$R \to \infty$ 时 $\Sigma(R) \sim R^{-d}$，$U(R) \sim -\bar{c}_r R$，$B(R) \sim R^{-\bar{c}_b / \bar{c}_r}$。关键引理2.8通过归纳法证明了泰勒系数的显式表达式，避免了数值 shooting 方法。
2. 径向稳定性分析（Section 3） 主要障碍是处理自相似坐标下的非线性演化方程：
其中 $\mathcal{N}$ 为非线性项。作者创新性地构造加权能量泛函：
权重 $w_j(R) = e^{-\beta R^2} R^{2j} (1 + R^{2m})$（$\beta > 0$，$m$ 大）用于控制 $R \to 0$ 和 $R \to \infty$ 的行为。关键估计在于处理边界项：在 $R=0$ 处利用 $\delta \Sigma \sim R^{2N}$ 的相容性条件消除奇性；在 $R \to \infty$ 处利用 $w_j(R)$ 的指数衰减控制线性项。引理3.12证明了能量耗散不等式：
其中 $\kappa > 0$ 由线性化算子的谱隙保证。通过Gronwall不等式导出指数稳定性，模去一维相容性条件（对应 $\partial_R^{2N} B(0)$ 的自由度）。
3. 非径向稳定性分析（Sections 4-5） 突破径向对称的关键障碍是处理高维扰动空间的谱结构。作者采用球谐展开将扰动分解：
其中 $Y_{lm}$ 为球谐函数。对 $l=0$（径向模式），稳定性由Section 3保证；对 $l \geq 1$，作者证明：(i) $l=1$ 模式对应平移对称性，其不稳定性由初始数据的质心条件 $\int x^\alpha (\rho_0 - \rho_{\text{exact}}) dx = 0$（$|\alpha|=1$）消除；(ii) $l \geq 2$ 模式的线性化算子具有正谱隙。关键工具是构造 $l$ 模式的加权能量：
其中 $\lambda_l = l(l+d-2) > 0$ 为球拉普拉斯特征值。引理4.7证明了对 $l \geq 2$ 有 $\frac{d}{d\tau} \mathcal{E}_l \leq -( \kappa + c_l) \mathcal{E}_l$（$c_l > 0$），而 $l=1$ 模式需满足相容性条件 $\int_0^\infty \delta \Sigma_1(R,0) R^d dR = 0$ 才稳定。定理5.1通过高阶泰勒展开精确刻画了不稳定流形：初始数据需满足 $\partial^{2N} \delta B(0) = 0$（$N$ 阶相容性条件），其维度由 $2N$ 阶导数的自由度决定。
4. 水动力学极限应用（Section 1.6.2） 创新性地将内爆解应用于Landau方程的流体极限。关键发现是：非等熵剖面的 $\bar{c}_b > 1/3$（$d=3$，$\gamma=5/3$）使水动力学极限条件 $\gamma_{kin} \bar{c}_u < 1/2$ 对所有 $\gamma_{kin} \in [-3,2]$ 成立（对比Merle等人的等熵解要求 $\gamma_{kin} > \sqrt{3}$）。通过选取 $N$ 使 $\bar{c}_u = \frac{2}{d+2N+1}$ 足够小，可满足任意 $\gamma_{kin}$ 的极限条件，特别地 $N=1$ 时覆盖 $\gamma_{kin} \in \{-3,0,1\}$（Coulomb、Maxwell分子、硬球）。
4. 与之前工作的比较
本文与既有工作的比较如下：
1. 与Guderley冲击波内爆的比较（Section 1.8）
- 相似指数确定性：Guderley解的相似指数需数值求解声学点非线性特征值问题（Lazarus[45]报告精度 $10^{-8}$ 需数小时计算），而本文剖面的 $\bar{c}_r = \frac{d+2N-1}{d+2N+1}$ 为显式闭式，计算复杂度 $O(1)$。例如 $d=3$、$\gamma=7/5$ 时，Guderley指数需迭代求解 $\lambda \approx 0.825$，而本文 $N=1$ 时 $\bar{c}_r = 3/5 = 0.6$ 直接可得。
- 光滑性与物理一致性：Guderley解含冲击前沿（密度有界但压力爆破），且假设零压力静止核心（与Navier-Stokes零粘性极限矛盾）；本文解 $C^\infty$ 光滑（无冲击），压力仅在原点为零（单点相容性条件），严格满足 $\rho > 0$ 且与Boltzmann方程流体极限兼容。数值实现上，Guderley解需冲击追踪和零压力区域正则化，而本文解可直接用标准有限差分法模拟（Remark 4.8）。
- 稳定性理论：Guderley解仅有线性稳定性分析[51,15]（支持径向稳定但非径向不稳定[6,68]），无非线性稳定性证明；本文首次给出非线性稳定性严格证明（Theorem 3.15），且对非径向扰动提供完整刻画（Theorems 4.4,5.1）。代价是本文解仅适用于 $\gamma \leq 2d+1$（虽覆盖物理相关值）。
2. 与Merle等光滑等熵内爆的比较（Section 1.8）
- 熵条件：Merle等解是等熵的（$\bar{c}_b = 0$），压力必与密度同步爆破；本文解真正非等熵（$\bar{c}_b > 0$），允许压力有界而密度爆破（$d=2$、$\gamma=2$ 时）。
- 相似指数显式性：Merle等解的相似指数需隐式求解（[48]中通过零点确定），而本文 $\bar{c}_r$、$\bar{c}_b$ 为显式闭式，便于数值实现。
- 稳定性：[3]的数值分析表明Merle等解在径向扰动下不稳定，而本文基态剖面在径向对称下严格稳定（Theorem 3.15），且对非径向扰动提供相容性条件（Theorems 4.4,5.1）。特别地，$N=1$ 时稳定性最优。
- 水动力学极限：Merle等解要求 $\gamma_{kin} > \sqrt{3}$ 才能满足自相似流体极限条件，而本文解通过非等熵项使条件 $\gamma_{kin} \bar{c}_u < 1/2$ 对所有 $\gamma_{kin} \in [-3,2]$ 成立，覆盖Coulomb碰撞（$\gamma_{kin} = -3$）等重要情形。
3. 与数值模拟应用的比较
- 径向对称数值代码：本文解作为基准优于Guderley解：(i) 相似指数显式（免去shooting方法）；(ii) 无冲击（免去前端追踪）；(iii) 压力仅在单点为零（避免零压力区域的数值困难）。
- 非对称数值代码：本文提供初始数据的显式相容性条件（如 $\partial^{2N} B(0) = \text{const}$），确保非径向稳定内爆。例如 $d=3$、$\gamma=5/3$ 时，初始数据需满足 $O(N^3)$ 个相容性条件（$N=1$ 时仅需6个），即可产生含非零涡度和熵变化的稳定内爆（Remark 4.8）。
4. 代价与局限
- 参数范围限制：非径向稳定性仅对单/双原子气体（$\gamma = \frac{d+2}{d}$, $\frac{2d+1}{2d-1}$）完全刻画，其他 $\gamma$ 值需进一步分析。
- 自相似性假设：解仅适用于自相似内爆，无法处理远场为非真空常数的情形（需渐近自相似分析[49,7]）。
- 高正则性要求：稳定性证明需初始扰动在 $C^M$ 范数下小（$M$ 足够大），可能排除某些物理相关但正则性较低的扰动。
2. 最优控制理论中不稳定的自由边界问题：存在性与正则性
Unstable free boundary problems in optimal control theory: existence and regularity
作者: Lorenzo Ferreri, Idriss Mazari-Fouquer, Raphaël Prunier
类别: math.AP, math.OC
发布日期: 2026-05-01
arXiv ID: 2605.00694v1
摘要翻译: 我们建立了源于数学物理与数学生物学的带约束最优控制问题的首个一般性正则性结果。具体而言，我们证明：对一大类形如“在约束 $0\leq m\leq 1$ a.e. 下最大化 $\int ψ(Θ_m)-c\int m$，其中 $-ΔΘ_m=mΘ_m+B(x,Θ_m)$”的问题，其解 $m^*$ 是 bang-bang 的，即 $m^*=χ_{E^*}$，且 $\partial E^*$ 在一个 $(d-2)$-维子集之外是光滑的。此外，我们还证明：对体积约束问题“在约束 $0\leq m\leq 1$ a.e. 与 $\int m=m_0$ 下最大化 $\int ψ(Θ_m)$，其中 $-ΔΘ_m=mΘ_m+B(x,Θ_m)$”，其解亦为 bang-bang 形式，即 $m^*=χ_{E^*}$，且在二维情形下，$\partial E^*$ 是有限条光滑曲线的并集。该结论通过将问题约化为一个不稳定的自由边界问题而获得；此类问题的正则性分析最早由 Monneau & Weiss 以及 Chanillo, Kenig & To 开创。在我们的设定中，自由边界并非极小化型，且状态函数的拉普拉斯算子 $ΔΘ_m$ 变号，这带来了显著困难，尤其体现在爆破解（blow-ups）的非退化性上。为此，我们发展了一种新方法，融合了最优控制理论、自由边界理论与测度论的工具，以建立该自由边界的正则性。
Comment: 61 Pages
链接: http://arxiv.org/abs/2605.00694v1
1. 研究对象和背景
本文研究一类源于数学物理和数学生物学的双线性最优控制问题，其核心是分析以下两类问题的解的正则性：约束问题 $\max_{m \in \mathcal{M}} \int_{\mathbb{T}^d} \psi(x, \Theta_m)$ 和惩罚问题 $\max_{m \in \mathcal{M}_0} \int_{\mathbb{T}^d} \psi(x, \Theta_m) - c \int_{\mathbb{T}^d} m$，其中状态方程为
且 $\Theta_m \geq 0$, $\Theta_m \not\equiv 0$。这里 $\mathcal{M}_0(\mathbb{T}^d) := \{m \in L^\infty(\mathbb{T}^d) : 0 \leq m \leq 1 \text{ a.e.}\}$，$\mathcal{M}(\mathbb{T}^d) := \{m \in L^\infty(\mathbb{T}^d) : 0 \leq m \leq 1 \text{ a.e.}, \int_{\mathbb{T}^d} m = m_0\}$，非线性项 $B$ 满足假设 (A1)：$B \in L^\infty(\mathbb{T}^d; C^2(\mathbb{R}))$，$B(x, 0) = 0$，$u \mapsto B(x,u)/u$ 在 $(0, +\infty)$ 上递减，且 $\lim_{u\to\infty}\sup_{x\in\mathbb{T}^d} B(x,u)/u = -\infty$。通过变换 $\theta_m := \ln(\Theta_m)$，状态方程可重写为 Hamilton-Jacobi 方程
其中 $Q$ 满足假设 (H1)：$\forall u \in \mathbb{R}$, $\max_{x\in\mathbb{T}^d} \partial_u Q(x, u) < 0$，$\lim_{u\to+\infty}\sup_{x\in\mathbb{T}^d} Q(x, u) = -\infty$，$\lim_{u\to-\infty}\sup_{x\in\mathbb{T}^d} |Q(x, u)| = 0$。通过引入切换函数 $\eta_{m^*}$，最优控制问题可转化为不稳定自由边界问题
其中 $m^* = \mathbf{1}_{E^*}$ 是 bang-bang 控制，$E^*$ 是最优集。这类问题与已有文献中的三类情形相关：(1) $g \equiv 0$（不稳定障碍问题）；(2) $g < 0$（复合膜问题）；(3) $g > 0$（不稳定两相膜问题）。然而，本文的关键困难在于无法保证 $\eta_{m^*}$ 是能量极小化解，且 $g$ 的符号不固定，这导致标准的非退化性分析失效。已有文献（如 Monneau & Weiss、Chanillo & Kenig 等）依赖于能量最小化性质或 $g$ 的固定符号来建立正则性，而本文需要发展新的方法来处理一般情形。特别地，在生物学应用中（如资源分布优化），bang-bang 性质和自由边界的正则性对理解种群分布的几何结构至关重要，但现有理论无法处理 $g$ 变号的情形。本文试图填补这一理论空白，建立最优控制问题解的 bang-bang 性质和自由边界的系统正则性理论。
2. 主要定理或主要结果
在假设 (A1)-(A2)-(A3) 或 (H1)-(H2)-(H3) 下，本文证明了以下主要结果：(1) 约束问题正则性：(Qconstr) 的任意解 $m^*$ 具有 bang-bang 性质，即 $m^* = \mathbf{1}_{E^*}$；进一步，$E^*$ 的测度论边界 $\partial_e E^*$（定义为 $\partial_e A := \{x \in \mathbb{T}^d : \forall r > 0, |A \cap B(x; r)| \cdot |(\mathbb{T}^d \backslash A) \cap B(x; r)| > 0\}$）在 $d=2$ 时是有限个 $C^{2,\gamma}$ 曲线的并集（定理 1(1) 和定理 2(1)）。(2) 惩罚问题正则性：对任意 $c > 0$，(Qpenal) 的任意解 $m^*$ 满足 $m^* = \mathbf{1}_{E^*}$；且 $\partial_e E^*$ 满足：(a) 当 $d=2$ 时，$\partial_e E^*$ 是有限个 $C^{2,\gamma}$ 曲线的并集；(b) 当 $d \geq 3$ 时，$\partial_e E^*$ 除 Hausdorff $(d-2)$ 维集合外是 $C^{2,\gamma}$ 的（定理 1(2) 和定理 2(2)）。(3) 奇异点分类与排除：在二维情况下，自由边界可分解为 $\partial E^* = \mathcal{R} \sqcup S^e_h \sqcup S^e_{nd} \sqcup S^e_d \sqcup S_d$，其中 $\mathcal{R}$ 为正则点集（$|\nabla \eta_{m^*}| \neq 0$），$S^e_h$ 为调和奇异点集（$\Psi(0+) = -\infty$，吹胀极限为非退化 2 次齐次调和函数），$S^e_{nd}$ 为非退化奇异点集（$\Psi(0+) > -\infty$，吹胀极限为非退化 2 次齐次函数），$S^e_d$ 为退化奇异点集（$\Psi(0+) = 0$，吹胀极限恒为零），$S_d$ 为退化部分。本文证明 $S^e_h = S^e_{nd,s} = S^e_d \cup S^e_{nd,p} = \emptyset$（定理 8 和命题 4.2.2、4.3.3），从而 $\partial_e E^* = \mathcal{R}$。(4) 中间密度点非退化性：对任意中间密度点 $x_0 \in \partial_e E^*$（满足 $\liminf_{r \to 0} \frac{|E^* \cap B(x_0; r)|}{|B(x_0; r)|} > 0$ 且 $\limsup_{r \to 0} \frac{|E^* \cap B(x_0; r)|}{|B(x_0; r)|} < 1$），存在常数 $C > 0$ 使得 $\sup_{B(x_0; r)} |\eta_{m^*}| \geq C r^2$（定理 10），这保证了吹胀极限的非退化性。
3. 研究方法、关键技术和核心工具
本文的核心技术障碍在于：(1) 无法依赖能量最小化性质或 $g$ 的固定符号，导致标准的非退化性分析失效；(2) 对于约束问题，允许的扰动集合较小，二阶最优性条件较弱。针对这些问题，作者发展了以下创新方法：首先，通过最优控制理论将原问题转化为不稳定自由边界问题。关键步骤是证明 bang-bang 性质（$m^* = \mathbf{1}_{E^*}$）和建立切换函数 $\eta_{m^*}$ 满足的方程
进而重写为自由边界问题 $- abla \eta_{m^*} = f \mathbf{1}_{E^*} - g \mathbf{1}_{\mathbb{T}^d \backslash E^*}$（引理 2.1.2）。其次，引入 Weiss 单调性公式 $\Psi(r) = \frac{1}{r^{d+2}} \int_{B_r(x_0)} |\nabla \eta_{m^*}|^2 - \frac{2}{r^{d+3}} \int_{\partial B_r(x_0)} \eta_{m^*}^2$，并证明其拟单调性：存在常数 $C > 0$ 使得 $r \mapsto \Psi(r) + C r$ 非减（定理 6），从而保证吹胀极限 $\eta_0(x) = \lim_{r_k \to 0} r_k^{-2} \eta_{m^*}(x_0 + r_k x)$ 的存在性（定理 7）。第三，对中间密度点发展新的非退化性分析方法：通过弱稳定性条件（命题 2.1.11）
结合傅里叶分析，证明若吹胀极限退化则导致矛盾（定理 10）。具体而言，将吹胀极限在极坐标下展开为 $\eta_0(r,\theta) = \sum_{k \in \mathbb{Z}} a_k r^{|k|} e^{ik\theta}$，利用弱稳定性条件导出系数 $a_k$ 的递推关系，结合 $f + g > 0$ 的假设证明 $a_k \neq 0$ 对某些 $k$ 成立，从而排除退化情形。第四，利用测度论方法连接中间密度点与奇异点：证明中间密度点在 $\partial_e E^*$ 上稠密（命题 4.2.1），且每个中间密度点附近自由边界为简单闭曲线，其对强稳定性条件的贡献有正下界（引理 4.3.1）。若存在奇异点 $x_0 \in S^e_d \cup S^e_{nd,p}$，则存在中间密度点序列 $x_{0,i} \to x_0$，生成无限多条互不相交的自由边界曲线 $\Gamma_i \subset \partial_e E^*$，与有限性矛盾（命题 4.3.3）。最后，对高维惩罚问题，利用维度约简论证和 Monneau & Weiss 的方法，证明奇异点集的 Hausdorff 维数至多为 $(d-2)$（定理 9）。这些技术的创新性在于将最优控制理论、自由边界理论和测度论有机结合，突破了传统方法对能量最小化或固定符号的依赖。
4. 与之前工作的比较
本文与已有工作的主要比较如下：首先，与 Monneau & Weiss [54] 的工作相比，他们研究 $g \equiv 0$ 且 $\eta_{m^*}$ 为能量最小化解的情形，证明了二维情况下自由边界的局部解析性及高维下除 $(d-2)$ 维集合外的光滑性。本文的关键推进在于不假设 $\eta_{m^*}$ 是最小化解，且允许 $g$ 变号，通过最优控制的二阶条件替代能量最小化来建立正则性。特别地，当 $g$ 变号时，Weiss 单调性公式的标准证明失效（因 $g$ 的符号影响导数符号），本文通过引入弱稳定性条件和傅里叶分析克服了这一障碍。其次，与 Chanillo & Kenig [19] 和 Chanillo, Kenig & To [20] 的复合膜问题研究相比，他们要求 $g > 0$ 以保证解的严格超调和性，从而获得非退化性。本文的突破在于处理 $g$ 变号的情形，代价是仅在二维约束问题和高维惩罚问题中获得完整正则性；对于高维约束问题，目前仍无法排除退化奇异点。值得注意的是，当 $g > 0$ 时，$-\frac{1}{2}x_2^2$ 在 $d=2$ 时是解且具有整条奇异点线，但本文通过最优控制的二阶条件排除了此类情形。第三，与 Soave & Terracini [62] 的不稳定两相膜问题研究相比，他们假设 $g > 0$ 以保证 Weiss 单调性公式的几乎单调性，而本文在 $g$ 变号时无法直接应用该公式，转而发展了基于中间密度点的新方法。本文的自由边界正则性结果在适用范围上更广（覆盖所有三种情形），但技术路径截然不同：他们依赖 $g$ 的固定符号导出吹胀极限的分类，而本文通过最优控制条件约束吹胀极限的结构（如命题 4.1.1 证明符号变化的吹胀极限必须具有特定星形结构）。第四，与 Andersson, Shahgholian & Weiss [1-4] 对非最小化解的研究相比，他们证明奇异点集的 Hausdorff 维数至多为 $(d-2)$，但未给出自由边界的正则性。本文的创新在于结合最优控制的几何约束（bang-bang 性质），在二维情况下证明了自由边界的全局 $C^{2,\gamma}$ 正则性。特别地，他们构造的十字形奇异自由边界（如 $\eta(x,y) = x^2 - y^2$）在本文框架下被排除，因为此类解不满足最优控制的二阶条件（命题 4.2.2）。最后，本文的局限性在于：(1) 对约束问题，目前仅在 $d=2$ 获得完整结果，高维情形下中间密度点的稠密性论证失效；(2) 需要较强的二阶最优性条件（定理 4），这在某些应用场景中可能不满足；(3) 无法处理 Hamilton-Jacobi 方程中非二次 Hamiltonian 的情形（见备注 1.1.5）。相比之下，惩罚问题因具有更强的局部最优性条件，在高维下也能获得除 $(d-2)$ 维集合外的正则性（定理 2(2b)），这反映了两类问题的本质差异：约束问题的扰动空间较小，导致非退化性分析更困难。
3. 具有临界数据的拟线性偏微分方程弱解的存在性与唯一性
Existence and uniqueness of weak solutions to quasilinear PDEs with critical data
作者: Sebastian Bechtel, Pascal Auscher
类别: math.AP, math.CA, math.FA
发布日期: 2026-05-01
arXiv ID: 2605.00439v1
摘要翻译: 我们建立了具有有界、一致连续初值数据的拟线性偏微分方程的整体、有界弱解的存在性与唯一性，并研究了其性质。此外，当初始数据仅为有界时，我们亦建立了有界弱解的存在性。
Comment: 19 pages. Comments welcome
链接: http://arxiv.org/abs/2605.00439v1
1. 研究对象和背景
本文研究拟线性抛物型偏微分方程：
其中系数函数 $a: [0, \infty) \times \mathbb{R}^n \times O \to \mathbb{R}^{n \times n}$ 满足假设A，$O \subseteq \mathbb{R}$ 为开集，初始数据 $u_0 \in \text{BUC}(\mathbb{R}^n)$ 满足 $\overline{\text{ran}(u_0)} \subseteq O$。假设A包含四个关键条件：(1) 局部椭圆性：对任意紧集 $K \subseteq O$，存在 $\lambda(K) > 0$ 使得 $a(t, x, y)\xi \cdot \xi \geq \lambda(K)|\xi|^2$；(2) 局部Lipschitz条件：对任意 $T > 0$ 和紧集 $K \subseteq O$，存在非递减常数 $C_L(T, K)$ 满足 $|a(t, x, y) - a(t, x, y')| \leq C_L(T, K)|y - y'|$；(3) 有界平衡：对任意 $T > 0$，存在非递减常数 $C_E(T)$ 使得 $|a(t, x, 0)| \leq C_E(T)$；(4) 局部一致连续性：对任意 $T > 0$ 和紧集 $K \subseteq O$，$a$ 在 $(0, T) \times \mathbb{R}^n \times K$ 上一致连续。$a$ 可能在 $y$ 上无界且退化，例如 $a(t, x, y) = y^m \cdot \text{Id}$（$m > 0$）对应方程 $\partial_t u - \text{div}(u^m\nabla u) = 0$，此时要求 $u_0 \in \text{BUC}(\mathbb{R}^n)$ 且 $\inf_{x \in \mathbb{R}^n} u_0(x) > 0$。有界弱解定义为：$u: (0, T) \times \mathbb{R}^n \to O$ 本质有界且 $\text{ran}(u) \subseteq O$，$\nabla u \in L^2_{\text{loc}}((0, T) \times \mathbb{R}^n)$，满足弱形式方程：
以及分布意义下的初值条件 $\langle u(t), \psi \rangle \to \langle u_0, \psi \rangle$（当 $t \to 0$）。该问题的重要性在于：有界性在抛物型尺度变换下是尺度不变的，因此研究临界正则性情形（即初始数据仅满足有界性而非更高正则性）具有根本意义。已有工作（如Ladyzhenskaya等[5]）主要处理高正则性情形，而本文聚焦于尺度临界最小正则性设置。此外，$\mathbb{R}^n$ 的无界性导致无法使用经典紧性论证，且初始数据质量通常无限，排除了 $L^1$ 技巧的应用。本文填补的空白是：在局部椭圆性、局部Lipschitz条件等最弱假设下，建立全局有界弱解的存在唯一性及正则性理论，特别处理了初始数据仅为有界（而非一致连续）的情形。关键创新在于引入加权Z空间处理低正则性，避免了对解的逐点控制，转而采用符合尺度分析的局部平均控制。
2. 主要定理或主要结果
定理1.4（拟线性方程适定性）：设 $a$ 满足假设A，$u_0 \in \text{BUC}(\mathbb{R}^n)$ 且 $\text{ran}(u_0) \subseteq O$，则存在唯一全局有界弱解 $u: (0, \infty) \times \mathbb{R}^n \to O$ 满足：(i) 正则性：$u$ 在 $[0, \infty) \times \mathbb{R}^n$ 上一致连续，且对任意 $T > 0$ 和 $q \in (1, \infty)$ 有
(ii) 稳定性：$\text{ran}(u) = \text{ran}(u_0)$ 且 $\|u(t) - u_0\|_{L^\infty} \to 0$（当 $t \to 0$）；若 $\lim_{|x| \to \infty} u_0(x) = c$ 存在，则 $\|u(t) - c\|_{L^\infty} \to 0$（当 $t \to \infty$）。推论1.5（有界数据存在性）：当 $u_0 \in L^\infty(\mathbb{R}^n)$ 且 $\text{ran}(u_0) \subseteq O$ 时，存在全局有界弱解，但失去唯一性、梯度正则性及 $t \to 0$ 时的稳定性。此时对 $a$ 的要求可弱化为：除局部椭圆性外，仅需 $a(t, x, \cdot)$ 在 $O$ 上连续，且 $a \in L^\infty([0, T] \times \mathbb{R}^n \times K)$（对任意 $T > 0$ 和紧集 $K \subseteq O$）。推论1.7（Hölder正则性）：对推论1.5中的解 $u$，由Nash局部正则性[8]得 $u \in C^\alpha((\delta', \infty) \times \mathbb{R}^n)$（$\alpha \in (0,1)$，$\delta' > 0$）；进一步由假设A(4)及经典抛物理论[5]，对任意 $\delta \in (\delta', T)$ 有 $u \in C^\alpha((\delta, T) \times \mathbb{R}^n)$（$\alpha \in (0,1)$）。
3. 研究方法、关键技术和核心工具
主要技术障碍包括：$\mathbb{R}^n$ 无界性导致无法使用经典紧性论证；初始数据质量无限使 $L^1$ 技巧失效；低正则性下 $\sqrt{t}\nabla u$ 的逐点控制过于苛刻。作者的创新点在于：1. 加权Z空间的构造与应用：定义 $Z^{\infty,q}_{-1/2}(T)$ 空间，其范数为
该空间通过局部平均控制替代逐点控制，符合抛物尺度 $t \sim |x|^2$ 的分析需求。关键性质是嵌套性：对 $q_0 \geq q_1$ 有 $Z^{\infty,q_0}_{-1/2}(T) \subseteq Z^{\infty,q_1}_{-1/2}(T)$，且 $L^\infty_{-1/2}(T) \subseteq Z^{\infty,q}_{-1/2}(T)$（$\forall q > 1$）。2. 非自治线性PDE的先验估计：命题3.4(e)证明当 $A$ 一致连续时，$\nabla u \in Z^{\infty,q}_{-1/2}(T')$ 且满足
其中 $C(T')$ 依赖于 $A$ 的模连续性和系数界。证明核心是迭代不等式：对 $r \in (1^*, \infty)$ 和 $0 < c' < c < 1 < d < d' < \infty$，有
通过有限次迭代降至 $p \in (1,2]$，再结合Caccioppoli不等式 $\|\nabla v\|_{L^p} \lesssim \|v\|_{L^2}$ 完成证明。3. 自治线性问题的表示公式：引理4.6给出关键表示
其中 $E_{A_0}$ 为自治方程解，$R_{A_0}$ 由引理4.3定义。该公式将一般问题转化为热方程的扰动，利用 $\nabla E_{\text{Id}}(u_0) \in L^\infty_{-1/2}$（引理4.4）及 $R_{A_0}$ 的有界性建立估计。4. 不动点论证框架：对拟线性问题，定义映射 $\Phi(v) = E_{a(\cdot,\cdot,v)}(u_0)$，利用Z空间中的梯度估计证明 $\Phi$ 在 $\text{BUC} \cap \{ \nabla u \in Z^{\infty,q}_{-1/2} \}$ 上是压缩映射。关键步骤是：当 $v \in \text{BUC}$ 时，$a(t,x,v)$ 满足非自治线性方程的系数条件，从而 $\nabla \Phi(v) \in Z^{\infty,q}_{-1/2}$ 且满足一致估计。唯一性由命题3.4(b)（$u(0)=0$ 蕴含 $u=0$）保证。5. 范围保持性质：通过Aronson基本解的Gauss界 $C_0(t-s)^{-n/2}e^{-c_0|x-y|^2/(t-s)} \leq \Gamma(t,x,s,y) \leq C_1(t-s)^{-n/2}e^{-c_1|x-y|^2/(t-s)}$ 证明 $\text{conv}(\text{ran}(u)) = \text{conv}(\text{ran}(u_0))$，结合 $u_0 \in \text{BUC}$ 的连通性得 $\text{ran}(u) = \text{ran}(u_0)$（备注3.5）。这为全局解的存在性提供关键支撑。
4. 与之前工作的比较
与既有工作相比，本文的推进主要体现在以下维度：1. 正则性假设的放宽：Ladyzhenskaya等[5, Ch. V, Thm. 8.1]要求 $a$ 和 $u_0$ 具有较高正则性（如 $C^{2+\alpha}$），而本文仅需 $u_0 \in \text{BUC}(\mathbb{R}^n)$ 和 $a$ 满足局部椭圆性等弱条件。特别地，定理1.4将初始数据正则性降至尺度临界情形（$\text{BUC}$ 是抛物尺度下的临界空间），这是文献中未解决的问题。推论1.5进一步将 $u_0$ 放宽至 $L^\infty(\mathbb{R}^n)$，仅需额外正则化 $a$，改进了[5, Ch. V, Rem. 8.1]的结果。2. 解空间的创新：传统方法（如[3]）依赖 $L^\infty$ 或 $C^\alpha$ 范数控制 $\sqrt{t}\nabla u$，但在 $\mathbb{R}^n$ 上对临界数据失效。本文引入的Z空间 $Z^{\infty,q}_{-1/2}$ 通过局部平均控制 $\sqrt{t}\nabla u$，既保留尺度不变性又适应低正则性。例如，对 $q > n+2$，$Z^{\infty,q}_{-1/2}$ 包含 $L^\infty_{-1/2}$ 但允许局部奇性，而经典方法要求 $\sqrt{t}\nabla u \in L^\infty$。3. 全局解的构造：在无界域 $\mathbb{R}^n$ 上，文献常依赖紧性论证（如[9]在有界域的工作），但本文通过证明解的先验一致连续性（命题3.4(d)）避免了紧性，直接建立全局解。关键突破是：当 $u_0 \in \text{BUC}$ 时，$u$ 的模连续性仅依赖于 $u_0$ 的模连续性、$\|u_0\|_\infty$ 及 $A$ 的系数界（命题3.4(d)），这保证了时间延拓的可能性。4. 唯一性机制的革新：传统 $L^1$ 技巧（如[7]）在 $\mathbb{R}^n$ 上对无限质量数据失效。本文利用命题3.4(b)（$u(0)=0$ 蕴含 $u=0$）和不动点框架建立唯一性，其核心是Z空间中的梯度估计：对 $v_1, v_2 \in \text{BUC}$，有 $\|\Phi(v_1) - \Phi(v_2)\|_{Z^{\infty,q}_{-1/2}} \leq C\|v_1 - v_2\|_\infty$，当 $C < 1$ 时压缩映射原理适用。5. 代价与局限：放宽 $u_0$ 正则性至 $L^\infty$ 时（推论1.5），失去唯一性、$t \to 0$ 的稳定性和梯度正则性。例如，对 $u_0 \in L^\infty \setminus \text{BUC}$，无法保证 $u(t) \to u_0$ 于 $L^\infty$ 中，且 $\nabla u$ 可能不属于 $Z^{\infty,q}_{-1/2}$。此外，系统情形（开放问题1）尚未解决，因命题4.1虽适用于系统，但拟线性系统的局部Lipschitz条件需矩阵范数控制，技术障碍更大。6. 与Nash理论的互补：推论1.7显示，当 $u_0 \in L^\infty$ 时，Nash正则性[8]仅给出 $t > 0$ 后的 $C^\alpha$ 正则性，而本文定理1.4在 $t=0$ 时提供 $\text{BUC}$ 稳定性。这表明：对于临界数据，全局正则性需初始数据满足 $\text{BUC}$，而 $L^\infty$ 数据仅能获得后验正则性。7. 与Aronson工作的区别：在证明范围性质时，本文不能直接使用[1]的唯一性结果（因其需解的符号信息），而是通过命题3.4(b)和Gauss界构造比较函数 $v = u - a$，避免了符号假设（备注3.6）。这使范围性质在更弱条件下成立，适用于变系数情形。
4. 一维可压缩欧拉方程光滑真空解的梯度爆破
Gradient blowup of smooth vacuum solutions to 1D compressible Euler equations
作者: Juhi Jang, Jiaqi Liu, Nader Masmoudi
类别: math.AP
发布日期: 2026-05-01
arXiv ID: 2605.00407v1
摘要翻译: 我们考虑定义在半直线上的等熵可压缩欧拉方程，该方程描述了与静止真空边界接触的气体流体的运动。我们构造了一大类初始时刻光滑且平方可积的解；这些解在有限时间内，在边界附近退化为 $C^{1-μ}$ 正则性（其中 $μ\in [1/2,1)$），从而导致边界处的梯度爆破。该构造基于作者近期所构造的自相似等待时间解 \cite{JLN2025} 的稳定性分析。
Comment: 无
链接: http://arxiv.org/abs/2605.00407v1
1. 研究对象和背景
本文研究半直线上等熵可压缩欧拉方程在静止真空边界条件下的解行为：
其中密度 $\rho \geq 0$、速度 $u$ 和压力 $p \geq 0$ 满足 $\gamma$-律 $p = \frac{\gamma-1}{\gamma}\rho^\gamma$，$\gamma \in (1,3)$。声速 $c = \sqrt{\frac{dp}{d\rho}}$ 定义为 $c = \sqrt{(\gamma-1)\rho^{\frac{\gamma-1}{2}}}$，在非真空区域系统等价于：
考虑空间域 $\Omega = \{0 < x < \infty\}$，左边界 $x = 0$ 为静止真空边界，满足：
初始数据 $c(t_0,x) = c_0(x)$, $u(t_0,x) = u_0(x)$ 满足 $c_0(x) > 0$ ($x > 0$), $c_0(0) = 0$, $u_0(0) = 0$。问题的关键在于真空边界附近 $c_0(x)$ 的正则性：当 $c_0(x) \sim x$ 时（即 $\alpha = 1$），Liu-Yang [35] 证明了 Sobolev 空间中的短时间适定性；当 $c_0(x) \sim x^{1/2}$ 时（即 $\alpha = 1/2$），对应物理真空边界条件 $(c^2)_x \sim 1$，此时边界会移动，需考虑真空自由边界问题。[25] 指出当 $\alpha = 1$ 时边界行为通常不能长期保持，可能经过等待时间后转变为物理真空边界（即等待时间解），而其他 $\alpha$ 值可能导致边界行为的瞬时变化。[27] 首次构造了从 $\alpha = 1$ 到 $\alpha = 1/2$ 的自相似等待时间解。本文试图填补的关键空白是：是否存在非自相似的光滑、平方可积初始数据也能产生类似奇性现象（如梯度爆破或静态边界演变为移动边界），以及这类解的普遍性如何。这一问题在浅水方程 ($\gamma = 2$) 和真空溃坝问题中也有重要意义。
2. 主要定理或主要结果
设 $\gamma \in (1,3)$，$\mu \in [\frac{1}{2}, 1)$，$\beta = \frac{1}{1-\mu} \geq 2$。令 $t_0 < 0$ 且 $|t_0|$ 充分小，则存在有限余维集的 $C^{\lfloor\beta\rfloor,\beta-\lfloor\beta\rfloor}$ 初始数据（在 $t = t_0$ 处），使得它们生成真空边界欧拉系统 (1.3)-(1.5) 在 $[t_0, 0)$ 上的强解，并在 $t = 0$ 时在 $x = 0$ 附近转变为 $C^{1-\mu}$-Hölder 连续，首次梯度爆破发生在 $(t,x) = (0,0)$。具体地，初始数据在边界附近满足：
其中 $a_0 > 0$，$b_i = b_i(a_2,\ldots,a_i)$。对 $0 < x \ll 1$，当 $t \to 0^-$ 时解满足：
其中 $\tilde{a}, \tilde{b} > 0$。速度与声速在每个 $t \in [t_0, 0)$ 的 $L^2$ 中有界，密度 $\rho > 0$ 在 $[t_0, 0) \times (0, \infty)$ 上保持正。精确表述见定理 2.7：对 $\mu \in [\frac{1}{2}, 1)$，当 $\mu \notin S = \{\mu \in [\frac{1}{2}, 1) \mid \beta - \lfloor\beta\rfloor = 0\}$ 时，初始数据集由条件 $a_1 = b_1 = 0$、$b_2 = b_2(a_2)$ 及 $a_0 > 0$ 确定；当 $\mu \in S$ 且 $\beta$ 为半整数时，需额外添加 $w_2 = 0$ 的约束 (2.40)。正则性指标 $n_\mu = \lfloor\beta + \frac{1}{2}\rfloor + 1 = \lfloor\frac{3-\mu}{2(1-\mu)}\rfloor + 1$ 确保了线性部分的严格耗散性。
3. 研究方法、关键技术和核心工具
核心障碍在于处理非自相似解的梯度爆破，需克服自相似剖面在远场无界、非整数 $\beta$ 引起的正则性问题，以及低频模态增长与高频稳定性之间的平衡。作者的创新点在于：首先引入黎曼不变量 $z = u - \frac{2}{\gamma-1}c$，$w = u + \frac{2}{\gamma-1}c$，其满足特征方程 $\partial_t z + \lambda_1 \partial_x z = 0$，$\partial_t w + \lambda_2 \partial_x w = 0$，其中特征速度 $\lambda_1 = \frac{\gamma+1}{4}z + \frac{3-\gamma}{4}w$，$\lambda_2 = \frac{3-\gamma}{4}z + \frac{\gamma+1}{4}w$。通过尺度对称性 $(u,c)(t,x) \mapsto \nu^{\delta-1}(u,c)(t/\nu, x/\nu^\delta)$ 定义自相似变量 $\tau = -\log(-t)$，$y = x/(-t)^\delta$（$\mu = 1/\delta$），系统转化为稳态形式：
关键引理 2.1 给出稳态解 $\bar{z}(y)$ 的渐近行为：当 $y \to 0^+$ 时 $\bar{z}(y) = -\frac{4\mu}{\gamma+1}y + c_1 y^\beta + O(y^{2\beta-1})$；当 $y \to \infty$ 时 $\bar{z}(y) = -c_2 y^{1-\mu} + o(y^{1-\mu})$。为构造平方可积解，作者通过输运方程 $\partial_\tau \chi + (y + \frac{\gamma+1}{4}\bar{z})\partial_y \chi = 0$ 定义 $\tau$-依赖截断函数 $\chi$，构造局部化剖面 $\mathring{z} = \chi \bar{z}$。扰动方程 (2.34) 的两个输运速度均为外向，对高阶导数产生阻尼效应。作者将解分解为：
其中 $\phi$ 为紧支截断函数，$z_i,w_i$ 为原点处泰勒系数，$(Z,W)$ 为余项。引理 3.1 通过高阶方程在原点的消失条件导出模态 ODE，确定导致时间衰减的初始数据集。对余项 $(Z,W)$，作者采用带奇异权的能量估计：分析系统 $(Z/y^{n_\mu}, W/y^{n_\mu})$ 及 $(\partial_y^{n_\mu} Z, \partial_y^{n_\mu} W)$，其线性部分在 $n_\mu = \lfloor\beta + 1/2\rfloor + 1$ 下严格耗散。非线性分析中，当 $\beta$ 为半整数时，$w_2 = 0$ 的额外约束 (2.40) 消除了低阶模态与剖面间的非可积交互项。技术关键在于 Hardy 不等式与插值不等式的精细应用，以及通过 bootstrap 论证调节不同范数的衰减/增长条件。
4. 与之前工作的比较
与既有工作相比，本文的推进主要体现在三方面：首先，[35] 仅处理 $c_0(x) \sim x$ ($\alpha = 1$) 的短时间适定性，而本文证明了此类初始数据可演化至 $C^{1-\mu}$ 正则性并发生梯度爆破；其次，[27] 构造的自相似等待时间解无界于远场，本文通过局部化技术首次构造了平方可积的非自相似解；第三，[30] 仅讨论自相似解的边界行为（瞬时移动或静止），本文揭示了导致梯度爆破的精确有限余维初始数据集。具体差异体现在：(1) 问题设定上，既有真空自由边界工作 [7,13,18-20,22,24,26] 均假设物理真空条件 $(c^2)_x \sim 1$，而本文处理的是初始满足 $c_0(x) \sim x$ ($\alpha = 1$) 的静止边界，证明其会经历等待时间后转变为梯度爆破；(2) 正则性范围上，[41] 仅研究 $\beta \in (1,2]$ 对应的 $C^{1/\beta}$ 奇性，本文将 $\mu$ 扩展至 $[1/2,1)$（即 $\beta \geq 2$），覆盖更广的初始正则性（$C^{\lfloor\beta\rfloor,\beta-\lfloor\beta\rfloor}$）；(3) 稳定性机制上，Burgers 自相似解 [8,11] 的 $\beta = 2k+1$ ($k \in \mathbb{N}$) 最稳定，而本文因边界约束 $\beta \geq 2$ 且 $\mu = 1/2$ ($\beta = 2$) 时最稳定，此时初始条件 $a_1 = b_1 = 0$、$b_2 = b_2(a_2)$ 均源于系统对称性（Remark 2.10）。代价与局限在于：(1) 仅考虑 $t < 0$ 的爆破前行为，$t > 0$ 的延续问题留待未来研究；(2) 当 $\beta$ 为半整数时需额外约束 $w_2 = 0$ 以消除非可积项；(3) 有限余维集意味着需精确调节初始数据，但作者证明这是必要条件（通过模态 ODE 的稳定流形）。值得注意的是，$\mu = 1/2$ 时解在爆破点附近呈现 $C^{1/2}$ 正则性，与 Burgers 方程 $\beta = 3$ 的 $C^{1/3}$ 奇性形成对比，凸显了真空边界对爆破结构的调制作用。
5. 守恒律粘性逼近的局部渐近模式
Local Asymptotic Patterns for Viscous Approximations of Conservation Laws
作者: Alberto Bressan, Laura Caravenna, Wen Shen
类别: math.AP
发布日期: 2026-04-30
arXiv ID: 2605.00189v1
摘要翻译: 双曲守恒律的解可通过多种方式逼近：粘性消失法、松弛法、离散或半离散数值格式、具有非局部通量的逼近等……对其中某些方法，已存在一般的 ${\bf L}^1$ 收敛性结果。本文旨在理解这些逼近解在双曲解具有奇性的点附近的局部行为，具体而言：沿激波的某点、两道激波相互作用的点、或新激波生成的点。
Comment: 无
链接: http://arxiv.org/abs/2605.00189v1
1. 研究对象和背景
本文研究标量守恒律方程 $u_t + f(u)_x = 0$ 的粘性逼近解在结构稳定奇异点附近的局部渐近行为，其中通量函数 $f$ 满足严格凸性条件 $0 < c_1 \leq f''(u) \leq c_2$ 对所有 $u \in \mathbb{R}$ 成立。对于一般 $C^3$ 初值 $\bar{u} \in \mathcal{U} \subset C^3(\mathbb{R})$（$\mathcal{U}$ 为 $C^3(\mathbb{R})$ 中可数个开稠密集的交），熵弱解是分片光滑的，具有有限多条激波曲线，且结构稳定的奇异点仅限于三类：(I) 沿单个激波的点；(II) 两个激波相互作用的点；(III) 新激波形成的点。研究对象为粘性逼近方程
关注当 $\epsilon \to 0$ 时，在上述三类奇异点附近，粘性解的局部行为是否可通过坐标缩放收敛到特定的“永恒解”（即在整个时空 $\mathbb{R}^2$ 上定义的解）。具体而言，给定奇异点 $(\tau, \xi)$，考虑缩放形式
问题在于是否存在 $\alpha, \beta, \gamma > 0$ 及点 $(\tau_\epsilon, \xi_\epsilon) \to (\tau, \xi)$，使得 $U_\epsilon$ 收敛到满足 $\epsilon = 1$ 时同类方程的全局解。这一研究填补了粘性逼近理论中局部行为刻画的空白：既有 $L^1$ 收敛结果（如文献 [4]）仅关注整体收敛性，而本文聚焦于奇异点附近的精细结构，为理解数值方法在激波附近的局部行为及不同逼近方法的关联提供理论基础。问题的重要性在于，结构稳定奇异点（由开稠密集初值生成）在物理情境中更可能被观测到，其局部渐近模式的刻画对数值模拟和理论分析均具有关键意义。
2. 主要定理或主要结果
定理 1.1：设假设 (A1)-(A2) 成立，$u = u(t, x)$ 为守恒律方程 (1.3)-(1.4) 的熵弱解，在点 $P = (\tau, \xi)$ 附近分片 Lipschitz 连续，$\Omega$ 为 $P$ 的后向域且 $\bar{u} \in C^1(\Omega_0)$。令 $(u_\epsilon)_{\epsilon > 0}$ 为粘性方程 (1.10) 的解族，初值满足
则当 $\epsilon \to 0$ 时：(i) 若 $(\tau, \xi)$ 为单激波点，存在 $(\tau_\epsilon, \xi_\epsilon) \to (\tau, \xi)$ 使得缩放函数
在 $\mathbb{R}^2$ 有界子集上一致收敛到粘性行波解 $U(t, x) = S(x - \lambda t)$，其中 $S$ 满足
且归一化条件 $S(0) = \frac{u_- + u_+}{2}$。(ii) 若 $(\tau, \xi)$ 为两激波相互作用点，存在 $(\tau_\epsilon, \xi_\epsilon) \to (\tau, \xi)$ 使得 $U_\epsilon$ 收敛到全局粘性解 $W(t, x)$，满足对任意 $\lambda_2 < \lambda < \lambda_1$，
其中 $S_1, S_2$ 为连接 $(u_-, u_*)$ 与 $(u_*, u_+)$ 的粘性行波剖面，$\lambda_1 > \lambda_2$ 为 Rankine-Hugoniot 速度。定理 1.2：设 $P = (\tau, \xi)$ 为新激波形成点，$\bar{u} \in C^4(\Omega_0)$，初值满足
则存在 $(\tau_\epsilon, \xi_\epsilon) \to (\tau, \xi)$ 及常数 $c, \sigma, \lambda$，使得缩放函数
在有界集上一致收敛到 $\Phi(t, x) = c Z(\sigma t, x - \lambda t)$，其中 $Z$ 为粘性 Burgers 方程的全局解
$z$ 满足无粘 Burgers 方程 $z_t + z z_x = 0$，$z(0, x) = -x^{1/3}$。
3. 研究方法、关键技术和核心工具
核心障碍在于精确刻画粘性解在奇异点附近的局部行为，需解决坐标缩放选择、全局解构造及收敛性证明三大挑战。首先，针对新激波形成点，作者通过变量缩放分析确定关键指数：将 $t = \epsilon^\alpha \tau$, $x = \epsilon^\beta y$, $u = \epsilon^\gamma v$ 代入守恒律方程，利用 $f'(u) = u + \psi_1(u)u^2$ 得
为使 $\epsilon \to 0$ 时退化为 Burgers 方程 $v_\tau + vv_y = 0$，需满足 $\gamma - \alpha = 2\gamma - \beta$ 与 $\gamma = \beta - \alpha > 0$。结合解 $x = -bu^3 + tu$ 的缩放不变性 $\epsilon^\beta y = -b\epsilon^{3\gamma}v^3 + \epsilon^{\alpha+\gamma}\tau v$，导出 $\beta = 3\gamma = \alpha + \gamma$，即 $\alpha : \beta : \gamma = 2 : 3 : 1$。对于粘性逼近，进一步推导出新激波情形的缩放关系 $\alpha = 1/2$, $\beta = 3/4$, $\gamma = 1/4$，即 $U_\epsilon(t, x) = \epsilon^{-1/4}[u_\epsilon(\cdots)]$。
其次，针对两激波相互作用情形，作者构造了全局“永恒解” $W(t,x)$。关键创新在于证明存在唯一解 $W$ 满足当 $\tau \to -\infty$ 时，
其中 $\lambda_2 < \lambda < \lambda_1$。此构造依赖于特征线方法与精细渐近分析：通过分析粘性激波剖面 $S_1, S_2$ 的渐近行为，结合比较原理与能量估计，证明 $W$ 的存在唯一性。特别地，作者利用 $S_1, S_2$ 的指数衰减性质（$|S_i(x) - u_\pm| \leq C e^{-\mu |x|}$）控制 $W$ 在 $t \to -\infty$ 时的误差。
第三，对于新激波形成情形，作者构造了粘性 Burgers 方程的全局解 $Z(t,x)$。关键步骤包括：(1) 证明无粘解 $z(-n,x)$ 是 $Z^{(n)}$ 的上下解，利用最大值原理得 $\|Z^{(n)}\|_{L^\infty} \leq C$；(2) 通过 Schauder 估计建立高阶导数一致界 $\|Z^{(n)}\|_{C^{2+\alpha}} \leq C$；(3) 应用 Arzelà-Ascoli 定理证明 $Z^{(n)} \to Z$ 在紧集上一致收敛。此构造的难点在于 $t \to -\infty$ 时 $z(-n,x)$ 的奇性（$z_x \sim x^{-2/3}$），作者通过精细的 $C^0$ 估计与加权 Sobolev 范数克服此障碍。
此外，作者发展了初值修正技术：第 6 节引理证明可在有界区间外修改初值，而不影响缩放解在 $(\tau,\xi)$ 附近的渐近性质。此技术将问题简化为特定初值情形（如新激波形成时取 $\bar{u}(x) = -x^{1/3}$），大幅降低分析复杂度。收敛性证明中，作者结合能量方法（如 $\int |u_\epsilon - U_\epsilon|^2 dx$ 的微分不等式）与紧性论证，验证极限函数满足目标方程。特别在定理 1.2 中，需利用 $C^4$ 初值条件控制缩放后的四阶导数项，确保 $\epsilon^{-1/4}$ 缩放下的稳定性。
4. 与之前工作的比较
本文与既有工作的比较可从五个维度展开。首先，相较于整体收敛结果（如文献 [4] 的 $L^1$ 收敛定理），本文聚焦局部行为刻画：前者仅证明 $\|u_\epsilon - u\|_{L^1} \to 0$，而本文精确描述了三类奇异点附近的缩放极限，揭示了粘性逼近在激波附近的精细结构。例如，定理 1.1(i) 证明单激波点处的缩放解收敛到粘性行波 $S(x - \lambda t)$，其剖面宽度为 $O(\epsilon)$，而整体收敛结果无法捕捉此尺度。
其次，与匹配渐近展开方法相比（如文献 [20] 在孤立激波附近、[1, 13] 在新激波形成点的构造），本文采用“反问题”视角：假设粘性逼近已给定，分析其渐近行为，而非主动构造逼近。关键差异在于时间范围——文献 [1, 13] 仅处理奇点形成前（$t < \tau$）的行为，而本文通过全局“永恒解” $Z(t,x)$ 描述跨越 $t = \tau$ 的完整动态。例如，定理 1.2 中的 $Z(t,x)$ 定义于 $t \in \mathbb{R}$，包含 $t \to -\infty$ 时的无粘解与 $t \to +\infty$ 时的粘性激波，而 [1, 13] 仅涉及 $t < \tau$ 的局部展开。
第三，初值光滑性要求体现推进深度：定理 1.2 需 $\bar{u} \in C^4(\Omega_0)$，而文献 [1, 13] 对新激波形成点通常要求 $C^3$ 光滑性。此差异源于本文需控制缩放后的高阶导数（如 $\epsilon^{-1/4}$ 缩放下 $u_{xxx}$ 的贡献），作者在 Remark 1.3 中明确指出，若初值收敛速度不够快（如仅 $C^0$ 收敛），则 $\|u^n(t, \cdot) - u(t, \cdot)\|_{C^3} \to 0$ 可能失效，这解释了数值模拟中网格尺寸与时间步长的耦合需求。
第四，适用范围的扩展与代价：本文结果原则上可推广至其他逼近方法（如松弛逼近、非局部通量逼近），但仅严格处理了粘性情形。相较之下，文献 [26, 29, 33, 34] 研究的松弛系统
其局部行为可能需不同缩放指数（如 $\alpha, \beta$ 依赖松弛时间尺度 $\epsilon$）。本文的代价是限于标量方程与结构稳定奇异点，未处理多维情形或三个以上激波相互作用的复杂奇点。
最后，与分类奇异点文献的联系：作者在引言中引用 [22, 6, 10, 28] 对 PDE 奇异点的分类工作，将此类分类从无粘解推广至粘性逼近层面。例如，[22] 证明标量守恒律的结构稳定奇点仅三类，本文则进一步刻画了每类奇点对应的粘性逼近局部模式：单激波点对应粘性行波（TYPE 1），激波相互作用点对应两行波合并的永恒解（TYPE 2），新激波形成点对应 Burgers 全局解（TYPE 3）。此推进揭示了粘性项如何“平滑”无粘奇点：无粘解在 CASE 3 中的缩放极限为分片光滑解（含间断），而粘性解的缩放极限 $Z(t,x)$ 处处光滑，体现了粘性对奇点的正则化效应。
了解领域论文（摘要翻译）
1. Klein-Gordon-Zakharov 系统中偶极子解的存在性
Existence of dipoles of Klein-Gordon-Zakharov system
作者: Vicente Alvarez, Amin Esfahani
类别: math.AP
发布日期: 2026-05-01
arXiv ID: 2605.00805v1
链接: http://arxiv.org/abs/2605.00805v1
摘要翻译
中文标题: Klein-Gordon-Zakharov 系统中偶极子解的存在性
摘要翻译: 本文研究 Klein-Gordon-Zakharov 系统解的长时间行为。我们证明存在一类具有特殊结构的解，称之为偶极子解，即存在一个解 $\vec{u}$，使得
其中 $\vec{R}_{k}$ 对每个 $k$ 表示一个孤立波，其位置由平移 $z_k$ 刻画，且满足
我们的方法首先聚焦于该系统相关 Hamiltonian 算子的谱分析；随后，我们旨在建立一个强制性估计（coercivity estimate），从而导出保证该解存在的条件。需特别指出的是，在本问题中，我们的目标是通过求解终值问题（final data problem）来构造近似解；再借助一致估计与紧性结果，利用密度论证推导出所期望的结论。
2. 加权Sobolev迹算子的紧性与非线性Steklov问题
Compactness of weighted Sobolev trace operators and non-linear Steklov problems
作者: Alexander Menovschikov, Alexander Ukhlov
类别: math.AP
发布日期: 2026-05-01
arXiv ID: 2605.00795v1
链接: http://arxiv.org/abs/2605.00795v1
摘要翻译
中文标题: 加权Sobolev迹算子的紧性与非线性Steklov问题
摘要翻译: 我们利用Sobolev空间上的复合算子，证明了向外尖角区域中加权Sobolev迹算子的紧性。该结果使得我们能在恰当的函数框架下对向外尖角区域中的非线性Steklov问题进行严格表述，并建立其非平凡解的存在性。
3. 带BMO项与分数阶拉普拉斯算子的Gagliardo-Nirenberg型不等式
Gagliardo-Nirenberg type inequalities with a BMO term and fractional Laplacians
作者: Dung Le
类别: math.AP
发布日期: 2026-05-01
arXiv ID: 2605.00701v1
链接: http://arxiv.org/abs/2605.00701v1
摘要翻译
中文标题: 带BMO项与分数阶拉普拉斯算子的Gagliardo-Nirenberg型不等式
摘要翻译: 通过将局部导数替换为分数阶拉普拉斯算子 $(-\Delta)^{\alpha/2}$，改进了带BMO项的全局（强）Gagliardo-Nirenberg不等式；同时也给出了相应的局部版本。
4. 标量守恒律熵解的非局部逼近原理
Nonlocal Approximation Principle for Entropy Solutions of Scalar Conservation Laws
作者: Alexander Keimer, Lukas Pflug
类别: math.AP
发布日期: 2026-05-01
arXiv ID: 2605.00635v1
链接: http://arxiv.org/abs/2605.00635v1
摘要翻译
中文标题: 标量守恒律熵解的非局部逼近原理
摘要翻译: 我们在 $\mathbb{R}$ 上建立了标量守恒律熵解的一般非局部逼近原理。更准确地说，我们证明：对非负初值数据，其熵解可表示为相应标量非局部守恒律解的弱星极限。该非局部守恒律的通量函数依赖于密度的适当空间平均。证明基于 Hamilton–Jacobi 层面的重构：通过对原函数（primitives）进行操作，我们利用粘性解的稳定性性质识别出极限；随后借助 Hamilton–Jacobi 方程与标量守恒律之间的经典关系，恢复出熵解。我们进一步证明：在作适当平移后，该逼近可延拓至变号初值数据；并对凸通量情形，以非局部核的一阶矩为基准，证明了定量收敛估计。该结果使得我们能够借助满足质量传播有限速度这一双曲守恒律关键特征的非局部逼近来定义一般通量下的熵解。
5. 亚临界情形下具两个变系数的半线性梁方程柯西问题的渐近轮廓
Asymptotic profiles for the Cauchy problem of semilinear beam equation with two variable coefficients in the subcritical case
作者: Mohamed Ali Hamza, Yuta Wakasugi, Shuji Yoshikawa
类别: math.AP
发布日期: 2026-05-01
arXiv ID: 2605.00550v1
链接: http://arxiv.org/abs/2605.00550v1
摘要翻译
中文标题: 亚临界情形下具两个变系数的半线性梁方程柯西问题的渐近轮廓
摘要翻译: 本文研究亚临界情形下半线性梁方程的柯西问题。我们证明了相应抛物型问题的自相似解的渐近稳定性结果。该结果的证明基于在加权空间中建立精细的能量估计，并将这些估计重写为抛物型自相似变量形式。
6. KdV方程的几乎全局大偏差原理
Almost global large deviations principle for the KdV equation
作者: Riccardo Berforini D’Aquino, Ricardo Grande
类别: math.AP, math-ph, math.PR
发布日期: 2026-05-01
arXiv ID: 2605.00484v1
链接: http://arxiv.org/abs/2605.00484v1
摘要翻译
中文标题: KdV方程的几乎全局大偏差原理
摘要翻译: 我们研究环面上Korteweg-de Vries方程在平均尺度为 $ε$ 的随机初值下的极端波形成现象。我们建立了该方程解在其上确界上的大偏差原理，时间尺度可取任意长的多项式形式 $t \leq ε^{-n}$，其中 $n$ 为任意固定的自然数。该原理刻画了观测到异常大幅值解的概率的首阶渐近行为。在此可积情形下，动力学演化于不变环面上，其傅里叶模长几乎守恒，从而排除了基于共振能量交换的极端波形成机制。因此，大幅值仅可能通过相干结构或色散聚焦产生，而后者对应于大量相位的拟同步化。我们证明，在弱非线性情形下，色散聚焦起主导作用。我们的方法结合了Birkhoff正规形分析与概率论证，利用可积动力学的稳定性来控制长时间尺度下相位拟同步化的概率。
7. 高维 $n\geq 3$ 抛物-椭圆型 Keller-Segel 系统中 Dirac 奇异性的形成与行为
Formation and Behavior of Dirac Singularities in the Parabolic-Elliptic Keller-Segel System in Dimensions $n\geq 3$
作者: Gregor Flüchter
类别: math.AP
发布日期: 2026-04-30
arXiv ID: 2605.00110v1
链接: http://arxiv.org/abs/2605.00110v1
摘要翻译
中文标题: 高维 $n\geq 3$ 抛物-椭圆型 Keller-Segel 系统中 Dirac 奇异性的形成与行为
摘要翻译: 我们考虑在 $\mathbb{R}^n$ 中单位球上、满足齐次 Neumann 边界条件的非负径向对称解，其满足抛物-椭圆型 Keller-Segel 系统 \begin{align*} \left\lbrace \begin{array}{r@{}l@{\quad}l} &u_t=Δu-\nabla \cdot \big(u\nabla v\big),\ &0=Δv -μ+ u , \ \end{array}\right. \end{align*} 其中 $μ$ 是 $u$ 的空间平均值，且 $n\geq 3$。在二维情形下，已有结论表明：有限时间内爆破的解在模糊拓扑（vague topology）下收敛于 Dirac 剖面。相比之下，当 $n\geq 3$ 时，具有有限存在时间的爆破解似乎并不呈现此类集中行为。我们通过将解推广至对应于 $u$ 的累积密度的测度值解，将分析延拓至爆破时刻之后。在此框架下，我们证明了极小解的存在性
其中 $ρ$ 可积，$θ$ 单调递增且右连续。进一步，我们构造了一类初始数据，使得对某个 $t_0>0$ 有 $θ(t_0)>0$，从而确立了原点处 Dirac 质量的形成。与 $n=2$ 情形不同，该奇异质量并非瞬时跃升至正水平；相反，$θ$ 连续地变为正值。此外，$θ$ 在 $[t_0,\infty)$ 上严格递增，且全部质量渐近地被原点吸收。
不相关论文（基本信息）
1. 等离子体光子学中定量的对称性破缺与非线性谐波生成
Quantitative symmetry-breaking and nonlinear harmonic generation in plasmonics
作者: Hongyu Liu, Zhi-Qiang Miao, Jingfeng Yao, Chengxun Yuan, Guang-Hui Zheng
类别: math.AP
发布日期: 2026-05-01
arXiv ID: 2605.00683v1
链接: http://arxiv.org/abs/2605.00683v1
摘要: We develop a quantitative mathematical theory that offers new perspectives on nonlinear harmonic generation in plasmonic structures arising from symmetry breaking. Focusing on second harmonic generation–the most fundamental process and the most extensively studied owing to its practical significance–we establish a theoretical framework that can be readily extended to higher-order harmonics. We investigate the plasmonic system in the static regime using a columnar nanowire with $n$-fold rotational symmetry ($n \in \mathbb{N}$) and construct a phenomenological model in which the second harmonic response originates from nonlinear sources confined to a selvedge region near the surface. By introducing a notion of symmetry degree grounded in group theory, we precisely quantify the second harmonic generation in terms of multipolar contributions. Our theory complements existing physical descriptions of this practically important phenomenon and provides a rigorous account of how nonlinear optical efficiency depends on shape, size, symmetry, and defects in plasmonic structures.
2. Robin 零度与临界双曲悬链面的渐近几何
Robin nullity and asymptotic geometry of the critical hyperbolic catenoid
作者: Alexander Pigazzini
类别: math.AP, math.DG, math.SP
发布日期: 2026-05-01
arXiv ID: 2605.00617v1
链接: http://arxiv.org/abs/2605.00617v1
摘要: For each parameter $a>1$, the critical hyperbolic catenoid $Σ_a$ is a rotationally symmetric, free boundary minimal annulus in a geodesic ball $B^3(r(a))\subset\mathbb{H}^3$. The Morse index of $Σ_a$ is at least $4$ by Medvedev [7], who conjectures equality. In this paper we identify a new geometric and spectral phenomenon for the family $\{Σ_a\}_{a>1}$, which we call “parameter-criticality”, and study its consequences for the Robin spectrum. Specifically, we prove two main results: (I) Parameter-criticality (Theorem 1.5). The boundary radius $r(a)$ is non-monotone on $(1,\infty)$: it satisfies $r'(1^+)<0$ and $r(a)=\frac{3}{2}\log a+d_\infty+o(1)$ as $a\to\infty$ with $d_\infty=\log[Γ(1/4)/Γ(3/4)]-\frac{1}{2}\log(2π)$ (Theorem 1.4). Hence there exists a parameter-critical value $a^\sharp\in(1,\infty)$ with $r'(a^\sharp)=0$. (II) Robin nullity jump (Theorem 1.6). At every such $a^\sharp$, the Robin nullity of $Σ_{a^\sharp}$ satisfies $\text{nul}(L_{Σ_{a^\sharp}})\geq 3$, with an additional kernel element in mode $k=0$ generated by the parametric variation field $j_a=\langle\partial_aΦ_a,ν\rangle_L|_{a=a^\sharp}$, which we show is non-vanishing at the catenoid neck via the closed-form $j_a(0)=1/(2\sqrt{a^2-1})$. The argument requires the limit $r_0:=\lim_{a\to 1^+}r(a)$ characterized as the unique positive solution of the transcendental equation $\tanh(r_0)\,\tanh(2r_0/\sqrt{3})=\sqrt{3}/2$ (Theorem 1.3), giving a clean parametrization of the degeneration $Σ_a\toΣ_1$. The Robin nullity of $Σ_a$ in mode $|k|=1$ is shown to equal $2$ (Proposition 1.1); this extends to the hyperbolic setting the mode-by-mode Fourier decomposition technique of Devyver [2] for the Euclidean critical catenoid, and is used in the proof of (II) to identify the extra kernel as a mode-$k=0$ phenomenon.
3. Pauli-Fierz 动力学的经典极限
Classical limit of the Pauli-Fierz dynamics
作者: Lucas Jougla, Nikolai Leopold
类别: math-ph, math.AP
发布日期: 2026-04-30
arXiv ID: 2605.00228v1
链接: http://arxiv.org/abs/2605.00228v1
摘要: We study the Schrödinger evolution generated by the Pauli-Fierz Hamiltonian, a model for nonrelativistic quantum electrodynamics, in the classical limit $\hbar \rightarrow 0$. In this regime, we rigorously derive the Newton-Maxwell equations of classical electrodynamics as effective dynamics approximating the time evolution. Our result complements prior work by an alternative derivation that provides explicit estimates on the rate of convergence, justifying the validity of the approximation for a special class of initial data.
4. 保持结构的麦克斯韦方程组最优控制及其在源隐身中的应用
Structure-Preserving Optimal Control of Maxwell’s Equations with Applications to Source Cloaking
作者: Harbir Antil, Yaw Owusu-Agyemang, Rohit Khandelwal, Jimmie Adriazola, Denis Ridzal
类别: math.OC, math.AP, math.NA, physics.comp-ph
发布日期: 2026-04-30
arXiv ID: 2605.00212v1
链接: http://arxiv.org/abs/2605.00212v1
摘要: We develop a structure-preserving solution framework for the optimal control of the time-dependent Maxwell’s equations. Building on a well-posedness theory for a weak form of the forward problem, we first analyze a forward solver that couples Nédélec and Raviart–Thomas finite elements with Crank–Nicolson time stepping. The solver preserves the de~Rham structure, enforces a discrete Gauss law, exactly satisfies a per-time-step energy balance, and converges to the weak solution under low regularity assumptions on the problem data, which are dictated by the optimal control setting. To control the Maxwell system, we add the curl of a space-time current density as a source to Ampére’s law. The curl form yields charge conservation without auxiliary constraints. We prove the well-posedness and continuity of the control-to-state map, derive the adjoint system and a gradient representation for a tracking-type objective functional, and formulate a discrete optimization scheme that inherits structure preservation from the forward solver. Our discrete stationarity conditions are consistent with their continuous counterparts, and the discrete optimal controls converge, with mesh and time refinements, to the continuous optima. We demonstrate the merits of our optimal control formulation and the theoretical developments by numerically solving a series of source-cloaking model problems.