生成时间: 2026年05月27日 06:49:29

运行状态: 完成

计划处理论文数量: 19

已完成论文数量: 19

未完成论文数量: 0

重点关注论文数量: 4

了解领域论文数量: 13

不相关论文数量: 2

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

1. 任意余维数嵌入子流形上的实用张量微积分

Practical tensor calculus on embedded submanifolds of arbitrary codimension

作者: Vladimir Yushutin

类别: math.DG, math-ph, math.AP

发布日期: 2026-05-26

arXiv ID: 2605.27260v1

摘要翻译: 我们提出了一种完全外在的、无需参数化的张量微积分变体，适用于任意维度与余维数、带边界的嵌入子流形（可能随时间演化）。所提出的框架是分量无关的，对一般秩张量覆盖了若干基本概念，例如切向投影、外在导数与协变导数、外在Stokes公式以及Laplace-Beltrami算子。该发展框架的显著特征在于其算法递归性，以及由张量的特殊行表示法所提供的透明性——该表示法使人联想到完备树的递归数据结构。因此，所建议的张量微积分既适用于具体计算，也适用于理论分析；后者通过三个独立的应用在一般维度与余维数情形下得以展示：第一，我们推导出一个新的外在守恒律，即黎曼流形上不可压Euler流的外在动量消失原理（principle of vanishing extrinsic momentum）；第二，我们重新审视正余维数嵌入子流形上的Cauchy应力概念，并论证：仅当应力张量被限制为作用于切向取向上时，角动量守恒才蕴含应力张量的切向性与对称性；第三，对演化子流形，我们以一种外在方式引入一般秩张量场的物质导数，并推导出相应张量型Dirichlet能量变化率的表达式。本文提供了一套实用的记号与工具，可立即用于数学建模、几何感知偏微分方程（geometry-aware PDEs）的分析，以及嵌入子流形上的数值方法。

Comment: 无

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

1. 研究对象和背景

本文研究的是嵌入欧氏空间 $\mathbb{R}^n$ 中、具有任意维度 $d = n - m$ 与任意余维数 $m$（$1 \leq m < n$）的定向闭子流形 $M$ 上的张量场分析。核心对象是定义在 $M$ 上的 $(0,q)$ 型张量场 $T \in \mathcal{T}_q(M)$，其在每点 $x \in M$ 处取值为一个多重线性映射

$$ T(x) : (\mathbb{R}^n)^q \to \mathbb{R}, \quad (v_1,\dots,v_q) \mapsto T(x)(v_1,\dots,v_q), $$

并非依赖于局部参数化或切丛截面的内在描述，而是完全基于 $\mathbb{R}^n$ 的笛卡尔结构与子流形的外在几何——即由 $m$ 个 $C^2$ 水平集函数 $d_1,\dots,d_m \in C^2(\Omega_\delta)$ 同时为零所定义的嵌入：

$$ M = \{ x \in \mathbb{R}^n : d_1(x) = \cdots = d_m(x) = 0 \}, \quad \Omega_\delta = \left\{ x \in \mathbb{R}^n : \sum_{i=1}^m |d_i(x)| < \delta \right\}. $$

对应单位法向量场 $n_i(x) = \nabla d_i / \|\nabla d_i\| \in C^1(\Omega_\delta)^n$，投影矩阵 $N = \sum_{i=1}^m n_i n_i^T$ 与 $P = I - N$ 分别给出法向与切向正交投影。经典子流形微分算子如切向梯度、散度被定义为

$$ \nabla_M f = (\nabla f) P, \quad \nabla_M u = (\nabla u) P, \quad \mathrm{div}_M u = \mathrm{tr}(\nabla_M u), $$

其中 $\nabla f$ 是行向量，$\nabla u$ 是 $n \times n$ 矩阵，作用于列向量 $u$。关键几何量是平均曲率向量 $\kappa \in C^1(\Omega_\delta)^n$，定义为

$$ \kappa := -\mathrm{div}_M(\nabla_M r), \quad r(x) = x, $$

满足 $P\kappa = 0$，且在余维 $m=1$ 时退化为经典标量平均曲率 $\kappa = \mathrm{tr}(\nabla_M n)$。本文试图填补的关键空白在于：现有文献（如 [Bur+18]、[Bac+24]、[FS20]、[HL14] 等）对嵌入子流形上张量微积分的处理严重受限于余维数（常限于 $m=1$）或张量阶数（常限于 $q \leq 2$），且高度依赖坐标参数化或 Ricci 指标记号，导致计算不透明、难以推广至高阶张量与演化流形。本文旨在构建一个完全外在、无参数化、无分量、递归可计算的张量微积分框架，覆盖任意 $d$, $m$, $q$，并使其能直接服务于几何感知型偏微分方程（如曲面上的 Euler 方程、应力平衡方程、演化能量泛函变分）的建模、分析与数值实现。

2. 主要定理或主要结果

本文在任意维度 $d = n - m$、任意余维数 $m$、任意张量阶数 $q$ 的嵌入子流形 $M \subset \mathbb{R}^n$ 上，建立了如下严格证明的核心结论：（1）【递归张量表示与切向投影】定义了 $\mathcal{T}_q$ 上的“行表示”（row representation）：对 $T \in \mathcal{T}_q$，有 $T = \{T_*\} = \{T_x, T_y, \dots, T_z\}$，其中 $T_* \in \mathcal{T}_{q-1}$ 是其行分量；切向张量子空间 $\mathcal{PT}_q \subset \mathcal{T}_q$ 定义为满足 $T(v_1,\dots,Nv_k,\dots,v_q) = 0$（$\forall k$）的张量集合；切向投影 $P : \mathcal{T}_q \to \mathcal{PT}_q$ 由递归两步算法实现：先计算 $\tilde{T} = \{P T_*\}$，再令 $P T = \tilde{T}_n := \tilde{T} - \sum_{i=1}^m n_i^T \otimes \tilde{T}(n_i)$（定理 3.15）。（2）【外在 Stokes 公式推广】对任意 $q$ 阶张量场 $T \in \mathcal{T}_q(M)$ 及其切向部分 $P T$，有

$$ \int_M \mathrm{div}_M T = \int_{\partial M} (P T) \cdot t + \int_M (N T) \cdot \kappa, $$

其中 $t$ 是边界 $\partial M$ 的外余法向量场，$\kappa$ 是平均曲率向量（命题 2.2 推广至高阶，见第 4 节隐含推导）。（3）【外在动量守恒律】对黎曼流形 $M$ 上的不可压 Euler 流 $u \in \mathcal{PT}_1(M)$，证明存在全局外在动量守恒律：

$$ \frac{d}{dt} \int_M u = \int_{\partial M} (u \otimes u) \cdot t + \int_M p \kappa - \int_M (u \cdot \nabla_M u) = 0, $$

即“外在动量消失原理”，该结论不依赖于任何对称性假设，且适用于任意 $d,m$（定理 5.1）。（4）【Cauchy 应力的本构约束】对余维 $m > 0$ 的嵌入子流形，证明 Cauchy 应力张量 $\sigma \in \mathcal{T}_2(M)$ 若满足力与力矩平衡方程

$$ \mathrm{div}_M \sigma + f = 0, \quad \sigma^{[ij]} = \frac{1}{2}(\sigma^{ij} - \sigma^{ji}) = 0 \text{ 仅当 } \sigma \text{ 限制在切向输入上}, $$

则其切向性 $P \sigma P = \sigma$ 与对称性 $\sigma^T = \sigma$ 并非牛顿力学的必然推论，而取决于对“允许切割方向”的建模选择；更一般地，$\sigma$ 可为非切向、非对称，只要其法向-切向分量 $N \sigma P$ 在平衡中自动消去（定理 6.3）。（5）【材料导数与 Dirichlet 能量变分】对演化子流形 $M_t$，定义外在材料导数 $\frac{D}{Dt} T := \partial_t T + v \cdot \nabla_M T + \sum_{k=1}^q T(\dots, \nabla_M v \cdot v_k, \dots)$（定义 7.1），并证明任意 $q$ 阶张量场 $T$ 的张量 Dirichlet 能量 $E(t) = \frac{1}{2} \int_{M_t} |\nabla_M T|^2$ 的时间变化率为

$$ \frac{d}{dt} E(t) = \int_{M_t} \left( \nabla_M T \right) : \left( \frac{D}{Dt} \nabla_M T \right) + \int_{M_t} |\nabla_M T|^2 \, \mathrm{div}_M v + \int_{\partial M_t} |\nabla_M T|^2 \, v \cdot t, $$

其中 $:$ 表示 Frobenius 内积（定理 7.2）。所有结果均在 $C^2$ 正则性假设下成立，且未引入额外小参数或渐近展开。

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

本文的核心技术障碍在于：如何在完全摒弃局部坐标参数化与指标记号的前提下，为任意余维、任意阶张量建立一套既保持几何直观又具备严格代数结构与计算可行性的微积分框架。作者的真正创新点并非发明新算子，而是重构整个张量语言的基础语法——通过“递归行表示”（recursive row representation）将张量视为一种数据结构而非抽象多线性映射。具体而言：（1）【行表示的树状结构】将 $q$ 阶张量 $T \in \mathcal{T}_q$ 视为一棵深度为 $q$ 的完全 $n$ 元树，其叶节点存储 $n^q$ 个标量 $T(e_{i_1},\dots,e_{i_q})$，而内部节点对应行分量 $T_* \in \mathcal{T}_{q-1}$（定义 3.1，引理 3.2）。这一结构使所有张量运算（加法、数乘、外积、插入、收缩）均可通过递归应用笛卡尔内积 $u \cdot v$ 实现，彻底避免矩阵索引与求和约定。例如，左插入 $T(v)$ 定义为 $\sum_* v_* T_*$，右插入 $T \cdot v$ 定义为 $\{T_* \cdot v\}$（命题 3.9），而 Frobenius 内积 $S \odot T$ 定义为 $\sum_* S_* \odot T_*$（定义 3.4）。这种“计算即定义”的范式，使得张量微积分从符号游戏转变为可编程的递归算法。（2）【切向投影的递归消除法】传统方法需解约束优化或使用 Christoffel 符号，而本文提出纯代数的两步递归投影算法（命题 3.15）：第一步对每个行分量 $T_*$ 施加切向投影 $P T_*$；第二步通过减去法向分量 $\sum_i n_i^T \otimes T(n_i)$ 来强制 $T$ 在所有位置上对法向输入为零。该算法的正确性依赖于关键引理 3.14，它证明了若 $T$ 仅在第一个变量上不满足切向性，则 $T_n = T - \sum_i n_i^T \otimes T(n_i)$ 必然切向；而递归性保证了此性质对所有 $q$ 个变量成立。（3）【外在微分算子的统一构造】所有微分算子（梯度、散度、协变导数、Laplace-Beltrami）均被定义为笛卡尔导数 $\nabla$ 与切向投影 $P$ 的组合，并通过行表示递归实现。例如，$q$ 阶张量的切向梯度 $\nabla_M T$ 被定义为一个 $(0,q+1)$ 阶张量，其行分量为 $(\nabla_M T)_* = \nabla_M (T_*)$（第 4 节隐含）。这种构造天然兼容外在 Stokes 公式（命题 2.2），因其证明仅依赖于笛卡尔导数的迹恒等式 $\mathrm{tr}(\nabla(Nu)P) = u \cdot \kappa$ 与投影算子的幂等性 $P^2 = P$，无需任何曲率张量或联络系数。（4）【应用中的几何-物理耦合机制】三个应用均依赖同一核心工具：将物理定律（动量守恒、角动量守恒、能量守恒）在 $\mathbb{R}^n$ 中表述，再通过行表示与切向投影将其“挤压”回 $M$。例如，外在动量守恒律的推导，本质是将 $\mathbb{R}^n$ 中的欧拉方程 $\partial_t u + u \cdot \nabla u + \nabla p = 0$ 与 $\mathrm{div}_M u = 0$ 结合，利用 $\int_M \nabla p = \int_{\partial M} p t + \int_M p \kappa$（推论 2.3）及 $\int_M u \cdot \nabla_M u = -\int_M (u \otimes u) : \nabla_M P$，最终得到仅含 $M$ 上几何量的闭合表达式。这表明，行表示不仅是记号革新，更是连接外在欧氏分析与内在几何物理的结构性桥梁。

4. 与之前工作的比较

本文与既有工作在多个根本维度上形成鲜明对比：（1）【适用范围与通用性】先前文献普遍受限于特定余维或张量阶。[Bur+18] 仅处理标量情形（$q=1$）下的椭圆 PDE；[Bac+24] 虽涉及抛物方程与高阶张量，但显式计算仅限 $q \leq 2$；[FS20]、[HL14]、[SF21] 等均聚焦于 $m=1$ 或 $q=2$ 的物理模型（膜、壳、梁）。本文首次将理论框架推广至任意 $d$, $m$, $q$，且所有定义与定理均以统一递归形式陈述，无特例分支。（2）【方法论基础】主流方法（如 [Bet86]、[FT09]、[NV22]）严重依赖局部坐标卡与 Ricci 记号，导致计算复杂、几何意义模糊；[Gri+13]、[PAM25] 引入弱形式但仍未摆脱参数化。本文则彻底转向完全外在、无参数化、无分量范式，所有运算均基于 $\mathbb{R}^n$ 的标准基与笛卡尔导数，其代价是牺牲了内在微分几何的某些优雅性（如无显式联络），但换取了无可比拟的计算透明性与跨领域适用性。（3）【物理结论的普适性与修正性】在 Euler 方程应用中，[Mar+17]、[GV23] 所述的动量守恒依赖于 Killing 场对称性，而本文的“外在动量消失原理”是全局、无对称性假设、适用于任意 $m$ 的新守恒律，推广了 [CF88]、[MB01] 中关于 $\mathbb{R}^2/\mathbb{R}^3$ 中固定区域压力合力为零的经典事实。在 Cauchy 应力分析中，[GI75] 被广泛引用为切向性与对称性的数学依据，但本文指出其结论仅在预设切向切割（即应力张量仅接受切向法向量输入）下成立；本文通过允许任意方向切割，证明应力张量可为非切向、非对称，从而修正了 [Moh+13]、[JOR18]、[PSW21] 等文献中对该假设的隐含依赖。在演化能量分析中，[DE13b] 给出了标量 Dirichlet 能量的变分公式，本文则将其完整推广至任意 $q$ 阶张量，且材料导数定义明确包含 $\nabla_M v$ 对张量各分量的耦合作用，超越了 [BGN08]、[ES10] 等仅考虑标量或向量的情形。（4）【技术代价与局限】本文框架的代价在于：其几何对象（如平均曲率向量 $\kappa$）虽外在，但仍是 $M$ 的内蕴量；对高度振荡或奇异性较强的流形（如尖角、自交），水平集函数 $d_i$ 的 $C^2$ 正则性假设可能成为瓶颈；且行表示虽利于计算，但对高阶张量（$q \gg 1$）可能导致 $n^q$ 维度爆炸，此时需结合稀疏性或低秩近似。作者未显式展开与 [OR17]、[ORG09] 等弱形式工作的比较，亦未讨论与现代机器学习中流形张量网络的潜在联系。

2. Dyson 方程解的正则性及其应用

Regularity of solutions of the Dyson equation and applications

作者: Valentin Pesce

类别: math.AP, math.PR

发布日期: 2026-05-26

arXiv ID: 2605.26968v1

摘要翻译: 本文旨在研究 Dyson 方程解的正则性。在 Bertucci 等人 [3, 4, 5] 的工作中，借助粘性解理论引入了 Dyson 方程的一种新解的概念，并证明了解在 $L^\infty$ 中具有正则化性质。根据 Biane [6] 的工作，我们预期解在 $C^{1/3}$ 中具有正则化性质（且不能更好）。受 Bertucci 等人工作的启发，我们将利用偏微分方程方法证明：对几乎所有 $t \ge 0$，解确如预期属于 $C^{1/3}$。我们的方法还可将该结果推广至含漂移项的 Dyson 方程情形，此时半显式解未必已知。此外，我们将给出该结果的一个应用：证明周期 Dyson 方程的解在长时间下按 $L^\infty$ 范数收敛于圆周上的均匀分布，这在 [5] 中是一个尚未解决的问题。

Comment: 无

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

1. 研究对象和背景

本文研究一维实轴上的Dyson方程

$$ \partial_t u + \partial_x(u H[u]) = 0 \quad \text{in } (0, +\infty) \times \mathbb{R}, $$

其中 $H[u]$ 是Hilbert变换，定义为

$$ H[u](x) = \frac{1}{\pi} \mathrm{p.v.}\int_{\mathbb{R}} \frac{u(y) - u(x)}{x - y}\, dy. $$

初始数据 $u_0 \in \mathcal{P}(\mathbb{R})$ 为概率测度，且在主要分析中进一步假设 $u_0 \in \mathcal{P}_2(\mathbb{R})$（即具有有限二阶矩）。该方程源于随机矩阵理论：它是Dyson布朗运动特征值经验分布的宏观极限PDE，亦被称为自由概率论中热方程的“自由类比”。从物理与分析角度看，它是一个非局部连续性方程，其输运速度 $H[u]$ 依赖于解的整体结构，而非局部梯度；这种非局部性导致标准双曲或抛物正则化理论（如Kružkov熵解或标准扩散估计）不直接适用。关键未知量是解 $u(t,x)$ 在时间 $t>0$ 后的空间正则性：已知由Bertucci等[3,4,5]基于粘性解理论证明了 $L^\infty$ 正则化，即存在普适常数 $C>0$ 使得

$$ \|u(t,\cdot)\|_{L^\infty} \le \frac{C}{\sqrt{t}}, \quad \forall t > 0, $$

但该估计仅给出有界性，无法刻画Hölder连续性。而Biane[6]利用自由概率中的复分析方法（通过自由卷积将解表示为初值与半圆律的自由和）证明：对任意 $t>0$，解 $u(t,\cdot)$ 至少属于 $C^{1/3}(\mathbb{R})$，且该指数不可改进（例如对一般初值，不存在 $\alpha > 1/3$ 使得 $u(t,\cdot) \in C^\alpha$ 对所有 $t>0$ 成立）。这一 $C^{1/3}$ 正则性具有根本重要性：一方面，它标志着非局部输运方程所能达到的最优Hölder正则性尺度，反映了Hilbert变换的 $1/2$ 阶奇异积分本质（因 $H$ 将 $L^2$ 映至自身，但将 $C^\alpha$ 映至 $C^{\alpha-1/2}$，故 $\partial_x H[u]$ 的正则性损失为 $1/2$）；另一方面，它在长时行为分析中构成关键缺口——Bertucci等[5]已证周期Dyson方程解在 $L^p$（$1 \le p < \infty$）中收敛于圆周上的均匀分布，但 $L^\infty$ 收敛性悬而未决，其障碍正在于缺乏一致 $C^{1/3}$ 控制以支撑紧性论证。因此，本文试图填补的核心空白是：在纯PDE框架下，不依赖复分析或自由概率工具，严格证明Dyson方程解的 $C^{1/3}$ 正则性，并将其推广至含漂移项的更一般情形，从而为长时渐近提供可控制的函数空间框架。

2. 主要定理或主要结果

在初始数据 $u_0 \in \mathcal{P}_2(\mathbb{R})$ 的假设下，本文证明了以下主要结果：（i）解 $u$ 满足先验估计

$$ \int_{\mathbb{R}} u(t,x) \log u(t,x)\,dx + \int_0^t \|u(s,\cdot)\|_{\dot{H}^{1/2}}^2\,ds \le \int_{\mathbb{R}} u_0(x) \log u_0(x)\,dx, \quad \forall t > 0, $$

以及

$$ \frac{1}{2} \|u(h,\cdot)\|_{\dot{H}^{1/2}}^2 + \frac{4}{9} \int_t^h \|u^{3/2}(s,\cdot)\|_{\dot{H}^1}^2\,ds \le \frac{1}{2} \|u(t,\cdot)\|_{\dot{H}^{1/2}}^2, \quad \forall h \ge t > 0; $$

（ii）由此导出 $\int_t^T \|u(s,\cdot)\|_{C^{1/3}}^3\,ds < +\infty$ 对任意 $T \ge t > 0$ 成立，进而推出 $u(t,\cdot) \in C^{1/3}(\mathbb{R})$ 对几乎每个 $t > 0$ 成立（定理2.7）；（iii）对含漂移项的Dyson方程

$$ \partial_t u + \partial_x(u(H[u] + b)) = 0, $$

在漂移 $b \in L^\infty_t(W^{1,\infty}_x)$ 的假设下，证明了修正估计

$$ \frac{1}{2} \|u(h,\cdot)\|_{\dot{H}^{1/2}}^2 + \frac{4}{9} \int_t^h \|u^{3/2}(s,\cdot)\|_{\dot{H}^1}^2\,ds \le \frac{1}{2} \|u(t,\cdot)\|_{\dot{H}^{1/2}}^2 + C_1 (h-t) + C_2 \int_t^h \|u(s,\cdot)\|_{\dot{H}^{1/2}}^2\,ds, $$

其中 $C_1 = \sup_{s\ge0} \|\partial_x b(s,\cdot)\|_{L^\infty}^2$, $C_2 = 2 \sup_{s\ge0} \|\partial_x b(s,\cdot)\|_{L^\infty}$，并由此得到 $u(t,\cdot) \in C^{1/3}(\mathbb{R})$ 对几乎每个 $t > 0$ 仍成立；（iv）作为应用，证明了周期Dyson方程的解 $u_{\mathrm{per}}$ 满足 $\|u_{\mathrm{per}}(t,\cdot) - \frac{1}{2\pi}\|_{L^\infty(\mathbb{T})} \to 0$ 当 $t \to +\infty$，彻底解决了文献[5]中遗留的 $L^\infty$ 长时收敛问题。

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

本文的核心技术障碍在于：Dyson方程既非标准抛物型（无二阶耗散项），亦非经典双曲型（因 $H[u]$ 非局部且不满足Lipschitz连续性），传统De Giorgi-Nash-Moser方法在维度 $n=1$ 下存在严重技术困难[9]，而Biane的复分析方法难以推广至含漂移情形。作者的真正创新点在于构建了一套基于分数阶Sobolev范数与非线性变换的PDE能量估计框架。具体而言：（1）第一层估计（引理2.1）利用熵泛函 $\int u \log u$ 与 $\dot{H}^{1/2}$ 范数的耦合，通过分部积分与恒等式 $(−Δ)^{1/2} u = \partial_x H[u]$ 推出 $\frac{d}{dt} \int u \log u + \|u\|_{\dot{H}^{1/2}}^2 = 0$，此式揭示了 $\dot{H}^{1/2}$ 能量的耗散机制，但 $\dot{H}^{1/2} \not\hookrightarrow C^\alpha$ 故不足以得Hölder正则性。（2）第二层估计（引理2.2）是突破性步骤：作者发现对 $u^{3/2}$ 施加 $\dot{H}^1$ 范数能捕获 $C^{1/3}$ 结构，通过精细计算导出

$$ \frac{1}{2} \frac{d}{dt} \|u\|_{\dot{H}^{1/2}}^2 + \frac{4}{9} \|u^{3/2}\|_{\dot{H}^1}^2 + \frac{1}{2} \int ((−Δ)^{1/2} u)^2 u\,dx = 0. $$

此处的关键技巧是两次运用Cotlar恒等式 $H[f]^2 − f^2 = 2H[f H[f]]$ 及 $\partial_x H[f] = (−Δ)^{1/2} f$，将含 $\partial_x u \cdot H[u] \cdot (−Δ)^{1/2} u$ 的棘手项转化为 $\frac{1}{2} \int (((−Δ)^{1/2} u)^2 − (\partial_x u)^2) u\,dx$，再将 $(\partial_x u)^2 u$ 重写为 $\frac{4}{9} (\partial_x (u^{3/2}))^2$，从而自然引出 $u^{3/2} \in \dot{H}^1$。（3）正则性传递引理（引理2.3）建立了非线性映射 $f \mapsto f^{3/2}$ 的Hölder指数缩放关系：若 $f^{3/2} \in C^{1/2}$，则 $f \in C^{1/3}$，其证明依赖于幂函数 $x \mapsto x^{2/3}$ 的 $2/3$-Hölder连续性，即 $\|f\|_{C^{1/3}} \le C \|f^{3/2}\|_{C^{1/2}}^{2/3}$。（4）对于含漂移情形（命题2.8），核心障碍是 $\partial_x (u b)$ 项破坏能量耗散。作者引入Kato-Ponce交换子估计处理 $[(−Δ)^{1/4}, b] \partial_x u$，结合Gronwall不等式控制 $\|u\|_{\dot{H}^{1/2}}$ 的指数增长，确保 $\|u^{3/2}\|_{\dot{H}^1}^2$ 的积分控制依然有效。整个框架的鲁棒性体现在：所有估计均通过光滑逼近 $u^\varepsilon$（添加 $−\varepsilon \Delta u^\varepsilon$ 正则化）建立，再利用Aubin-Lions-Simon引理获得强紧性，最终通过相对熵技巧（以标准正态密度 $\gamma$ 为参考）克服 $\int u \log u$ 的非凸性，使极限过程严格成立。

4. 与之前工作的比较

本文与既有工作存在多维度实质性推进：（1）与Bertucci等[3,4,5]相比：后者基于粘性解理论确立了 $L^\infty$ 正则化 $\|u(t,\cdot)\|_{L^\infty} \le C / \sqrt{t}$，但未触及Hölder正则性；本文首次在纯PDE框架下证明 $C^{1/3}$ 正则性，且将时间集从“所有 $t>0$”（Biane[6]的复分析结论）弱化为“几乎每个 $t>0$”，但代价是获得了可显式控制的积分型估计 $\int_t^T \|u\|_{C^{1/3}}^3 < \infty$，这对长时分析至关重要。（2）与Biane[6]相比：Biane的方法依赖自由概率中的Stieltjes变换与Nevanlinna函数理论，其结论 $u(t,\cdot) \in C^{1/3}$ 对所有 $t>0$ 成立且指数最优，但高度特化于无漂移Dyson方程，无法处理 $b \neq 0$ 的情形；本文虽仅得“几乎处处”结论，却提供了完全PDE化的、可推广的证明路径，并成功将 $C^{1/3}$ 正则性延拓至含任意 $W^{1,\infty}$ 漂移的方程，这是自由概率方法目前无法企及的。（3）与Caffarelli等[9,10]的De Giorgi-Nash-Moser型工作相比：后者在 $\mathbb{R}^n$（$n \ge 2$）中证明了类似方程的Hölder正则性，但明确指出一维情形存在技术障碍（如非局部算子的临界尺度问题）；本文专为 $n=1$ 设计，绕过De Giorgi迭代，转而利用 $u^{3/2}$ 的 $\dot{H}^1$ 嵌入 $C^{1/2}$ 这一更直接的Sobolev嵌入路径，显著简化了证明。（4）与长时行为文献[5]相比：此前仅知周期解在 $L^p$（$p<\infty$）中收敛，$L^\infty$ 收敛性因缺乏一致Hölder模而悬置；本文提供的 $C^{1/3}$ 估计（特别是 $\int_t^\infty \|u\|_{C^{1/3}}^3 < \infty$）蕴含了解在 $C^{1/3}$ 中的紧性，结合周期性与唯一平衡态，最终确立 $L^\infty$ 收敛。代价在于：本文结果要求初值具二阶矩（$\mathcal{P}_2$），而[3,4,5]仅需 $\mathcal{P}(\mathbb{R})$；此外，含漂移情形的估计中出现了依赖 $\|\partial_x b\|_{L^\infty}$ 的指数增长因子，这反映了漂移对正则性传播的潜在破坏作用，是当前方法的固有局限。

3. KMP 过程在维数 $d\ge 2$ 下的病态大偏差

Pathological Large Deviations of the KMP Process in Dimension $d\ge 2$

作者: Daniel Heydecker

类别: math.PR, math.AP

发布日期: 2026-05-26

arXiv ID: 2605.26652v1

摘要翻译: 我们在离散环面 $\mathbb{T}_N^d$ 上研究 Kipnis–Marchioro–Presutti 过程的动力学大偏差。通过将候选速率函数重写为具有粗糙漂移项的线性双曲-抛物型方程，我们证明：在任意维数 $d\ge 2$ 下，病态轨迹以有限速率出现在大偏差中。在此过程中，我们严格验证了 Bertini-Gabrielli-Lebowitz 在任意维数下所导出的下界。

Comment: 无

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

1. 研究对象和背景

本文研究Kipnis–Marchioro–Presutti（KMP）热传导模型在离散环面 $\mathbb{T}_N^d = (N^{-1}\mathbb{Z}/\mathbb{Z})^d$ 上的动力学大偏差行为。该模型描述非负能量变量 $\xi_N(t,x) \in [0,\infty)$ 在格点上的随机热输运：对每条无向边 $\{x,y\}$，以速率 $N^2$ 触发泊松钟；当钟响时，按均匀分布 $P \sim \mathrm{Uniform}[0,1]$ 重分配总能量 $\xi(x)+\xi(y)$，即

$$ \xi(x) \mapsto P(\xi(x)+\xi(y)), \quad \xi(y) \mapsto (1-P)(\xi(x)+\xi(y)). $$

其经验测度定义为 $\pi^N(\xi_N(t)) := N^{-d} \sum_{x \in \mathbb{T}_N^d} \xi_N(t,x) \delta_x \in \mathcal{M}_+(\mathbb{T}^d)$，总能量守恒：$\langle 1, \pi^N(\xi_N(t)) \rangle = N^{-d} \sum_x \xi_N(t,x) =: e_N$ 恒定。在平衡初值 $\xi_N(0) \sim \nu^N_\rho$（各站点独立服从均值为 $\rho$ 的指数分布）下，宏观尺度极限满足热方程

$$ \partial_t u = \frac{1}{2} \Delta u, \quad x \in \mathbb{T}^d,\ t \in [0,T], \tag{1.1} $$

这是KMP模型水动力极限的核心结果，亦是宏观涨落理论（MFT）的出发点。然而，MFT预言的动力学大偏差率函数 $I_\rho$ 并非直接由 (1.1) 的线性扩散结构导出，而是源于一个形式上更复杂的Fokker–Planck型方程：

$$ \partial_t u = \frac{1}{2} \Delta u - \nabla \cdot (u^2 \nabla H), \quad H \in C^{1,2}_{t,x}, \tag{1.2} $$

其中 $H$ 是控制场（control field），而 $u^2 \nabla H$ 构成非线性漂移项。与排它过程或零范围过程（ZRP）不同，KMP的扩散系数 $D(u) = 1/2$ 是常数（线性），而迁移率 $\chi(u) = u^2$ 是二次的——这一根本不对称性导致大偏差尺度下涨落与扩散效应之间存在非平凡竞争：当 $u$ 局部极小或趋于零时，$u^2$ 的衰减快于线性扩散项，使得漂移主导局部演化，可能诱发奇异性。问题的关键在于：MFT所依赖的形式梯度流结构（即以 $J_{\mathrm{FP}}(u) := \frac{1}{2} \int_0^T \int_{\mathbb{T}^d} u(t,x)^2 |\nabla H(t,x)|^2 \, dx dt$ 为动态代价）是否完备？抑或存在不满足 (1.2) 弱解正则性要求、却仍以有限速率出现的大偏差轨迹？已有静态大偏差结果（如 [47, Theorem 1]）已确立熵泛函

$$ S_\rho(\xi_0) = \left\langle \frac{1}{\rho}, \xi_0 \right\rangle - \int_{\mathbb{T}^d} \log\left( \frac{u_0(x)}{\rho} \right) dx - 1, \tag{1.3} $$

其中 $u_0$ 是 $\xi_0$ 绝对连续部分的密度，但动力学部分仍悬而未决。Bertini–Gabrielli–Lebowitz 在 [12] 中曾形式推导出 $J_{\mathrm{FP}}$，但其严格下界证明依赖于超指数替换引理（superexponential replacement lemma），而该引理在KMP中因二次迁移率失效：现有方法（如 [40, 38]）要求可观测量 $F$ 有界且Lipschitz，但KMP大偏差所需的 $F$ 增长性为二次，且不变测度 $\nu^N_\rho$ 的指数尾无法提供足够截断。因此，核心空白在于：是否存在真正‘病态’（pathological）的轨迹——即不具连续时间代表、非弱解、甚至不满足 $L^2_{t,x}$ 可积性——却仍以 $e^{-N^d I}$ 量级概率出现？若存在，其奇异性如何随维度 $d$ 演化？这直接挑战MFT作为普适框架的完备性。

2. 主要定理或主要结果

本文通过构造性论证，严格确立了病态大偏差轨迹的存在性及其维度依赖性，主结果如下：

定理1（受限下界）：设 $\xi^N_t$ 为初值 $\xi^N_0 \sim \nu^N_\rho$ 的KMP过程，$R := X \cap \{ u \in C^\infty([0,T] \times \mathbb{T}^d) : \inf_{(t,x)} u(t,x) > 0 \}$ 为正则路径集，则对任意 $u \in R$ 及其开邻域 $U \ni u$，有

$$ \liminf_{N \to \infty} \frac{1}{N^d} \log \mathbb{P}(\pi^N(\xi^N) \in U) \ge -I_\rho(u), \quad \text{其中 } I_\rho(u) = S_\rho(u_0) + J(u), $$

且在 $R$ 上 $J_{\mathrm{FP}} = J$，其中 $J(u) := \frac{1}{2} \inf \{ \|g\|_{L^2_{t,x}}^2 : u \text{ 弱解 } \partial_t u = \frac{1}{2} \Delta u - \nabla \cdot (u g),\ g \in L^2_{t,x} \}$。

定理2（病态轨迹作为大偏差）：在相同设定下，存在以有限速率出现的病态轨迹：

(i) 一维情形 ($d=1$)：对任意 $t_0 \in [0,T]$，存在 $\Xi \in X$，其唯一连续代表 $\xi_t$ 满足：当 $t \neq t_0$ 时 $\xi_t$ 具光滑密度 $u_t(x)$；而 $\xi_{t_0}$ 关于Lebesgue测度奇异（如Dirac质量），且

$$ \inf_{U \ni \Xi} \liminf_{N \to \infty} \frac{1}{N} \log \mathbb{P}(\pi^N(\xi^N) \in U) > -\infty. $$

(ii) 二维情形 ($d=2$)：对任意满足 $\sum_k \epsilon_k^\gamma < \infty$（$0<\gamma<1$）的正数列 $\epsilon_k>0$ 及互异时刻 $t_k \in (0,T)$，存在 $\Xi \in X$，其唯一càdlàg代表 $\xi_t$ 满足

$$ \|\xi_{t_k} - \xi_{t_k^-}\|_{\mathrm{TV}} = 2\epsilon_k, \quad \forall k, $$

即在每个 $t_k$ 处发生总变差为 $2\epsilon_k$ 的跳跃，且

$$ \inf_{U \ni \Xi} \liminf_{N \to \infty} \frac{1}{N^2} \log \mathbb{P}(\pi^N(\xi^N) \in U) > -\infty. $$

(iii) 高维情形 ($d \ge 3$)：对任意 $c,M \in (0,\infty)$，存在 $J = J(\rho,c,M) < \infty$，使得对任意满足 $\Xi(dt,dx) \ge c \, dt dx$ 且总能量 $e(\Xi) \le M$ 的松弛测度 $\Xi$，有

$$ \inf_{U \ni \Xi} \liminf_{N \to \infty} \frac{1}{N^d} \log \mathbb{P}(\pi^N(\xi^N) \in U) \ge -J(\rho,c,M). $$

此表明在 $d \ge 3$ 时，任意满足下界和能量约束的松弛测度均可被病态轨迹逼近，且代价一致有界。

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

证明的核心障碍源于KMP大偏差率函数的双重脆弱性：一是分析层面，骨架方程 $\partial_t u = \frac{1}{2} \Delta u - \nabla \cdot (u g)$（记为(Skg)）中，仅 $g \in L^2_{t,x}$ 受控，而 $u$ 的正则性完全缺失；二是概率层面，超指数替换引理失效，无法将粗糙可观测量平滑化。作者的创新在于绕过全局上界，转而从率函数的下半连续包络（lower semicontinuous envelope） 出发，利用维度相关的Sobolev嵌入失效来构造病态序列。

关键工具一：对数熵耗散的稳定性估计。对光滑 $\psi$ 形式计算 $\frac{d}{dt} \int_{\mathbb{T}^d} \psi(u) dx$，结合Cauchy-Schwarz，得到闭合估计的充要条件是 $(\psi''(u))^2 u^2 \lesssim \psi''(u)$，解得 $\psi''(u) = u^{-2}$，即 $\psi(u) = -\log u$。由此导出关键不等式：

$$ S_\rho(u_T) + \frac{1}{2} \|\nabla \log u\|_{L^2_{t,x}}^2 \le S_\rho(u_0) + \frac{1}{2} \|g\|_{L^2_{t,x}}^2, \tag{1.11} $$

其中 $S_\rho$ 由(1.3)定义。此式表明，有限率 $J(u) < \infty$ 蕴含 $\|\nabla \log u\|_{L^2_{t,x}} < \infty$，但 $S_\rho(u)$ 在 $u$ 奇异时仍可有限（因 $S_\rho$ 渐近线性于 $u$），故 $u$ 可坍缩为奇异测度。

关键工具二：维度驱动的紧性失效。病态构造依赖 $\|\nabla \log u\|_{L^2_{t,x}}$ 所能提供的紧性强度随 $d$ 递减：

  • $d=1$：Sobolev嵌入 $W^{1,1}_x \hookrightarrow L^\infty_x$ 结合能量守恒，保证 $u \in L^2_{t,x}$，故(Skg)可弱理解；病态体现为单点奇异（如 $\delta$-质量），由集中序列 $u_n(t,x) = n \mathbf{1}_{[0,1/n]}(x)$ 实现。
  • $d=2$：Moser-Trudinger不等式给出 $\{ u \in L^1(\mathbb{T}^2) : \|\nabla \log u\|_{L^2_x} \le M, \langle u,1 \rangle \le M \}$ 在 $L^1_x$ 弱拓扑下紧，故 $u_t \in L^1_x$ 对 $dt$-a.e. $t$ 成立；但嵌入不足以保证 $u \in L^2_{t,x}$，导致乘积 $u g$ 未定义。作者构造 $u_n \in R$, $g_n$ 解(Skg)，使在渐近消失的时间片 $I_n$ 上 $\int_{I_n} \int_{\mathbb{T}^2} |g_n|^2 \to 0$，但 $\int_{I_n} \int_{\mathbb{T}^2} u_n g_n \not\to 0$，从而在极限中产生瞬时跳跃。
  • $d \ge 3$：$\|\nabla \log u\|_{L^2_x}$ 对 $L^1_x$ 弱拓扑甚至不再下半连续。作者构造 $u_n$ 为大量尖峰的叠加（如 $u_n = \sum_j a_{n,j} \phi_n(x-x_{n,j})$），满足 $S_\rho(u_n) \le C$, $\|\nabla \log u_n\|_{L^2_x} \to 0$；再设计 $g_n$ 快速移动尖峰，使动态成本 $\frac{1}{2} \|g_n\|_{L^2_{t,x}}^2 \to 0$，从而逼近任意满足 $\Xi \ge c \, dt dx$ 的松弛测度 $\Xi$。

4. 与之前工作的比较

本文与既有工作存在本质差异，主要体现在以下三方面：

1. 与经典大偏差框架（[40, 38, 49]）的对比：传统方法（如排它过程、ZRP）依赖超指数替换引理，要求可观测量 $F$ 有界Lipschitz，从而确保局部平均 $\bar{\xi}^N_\epsilon$ 与真实值接近。KMP因二次迁移率 $\chi(u)=u^2$，需处理 $F$ 二次增长的情形，而 $\nu^N_\rho$ 的指数尾无法提供 $\mathbb{E}_{\nu^N_\rho}[e^{\lambda \xi^2}] < \infty$，故该引理失效。本文不尝试修复此引理，而是首创‘受限下界+下半连续包络构造’双轨策略：定理1在 $R$ 上验证 $I_\rho$，而定理2则证明 $I_\rho$ 的下半连续包络远大于其原始定义域，揭示MFT框架的固有缺陷。此前文献（如 [12, 47]）仅形式推导 $J_{\mathrm{FP}}$，未触及病态性。

2. 与同类热传导模型（[14, 18, 32, 47]）的对比：BEP($m$)、Brownian Energy Process等模型虽也具非线性迁移率，但其迁移率 $\chi(u) = u^{1+2/m}$ 在 $m \to \infty$ 时退化为KMP的 $u^2$。本文首次严格证明，仅当 $\chi(u)$ 为二次时，维度 $d \ge 2$ 下才会出现跳跃与任意松弛测度逼近等强病态。例如，[18] 对非平衡KMP导出了准势（quasipotential），但假设轨迹连续；本文则证明该假设在 $d \ge 2$ 下不成立。此外，[47] 的形式推导未区分维度，而本文明确指出病态严重性随 $d$ 升级：$d=1$ 仅单点奇异，$d=2$ 出现可数跳跃，$d \ge 3$ 则完全丧失正则性控制。

3. 与粗糙漂移PDE理论（[20, 36, 43, 46]）的对比：本文骨架方程(Skg)要求 $g \in L^2_{t,x}$，远粗于文献中常见假设（如 [36] 要求 $g \in L^1_t W^{1,1}_x \cap L^\infty_t C^0_x$）。这种极端粗糙性使得标准解理论（如重整化解、熵解）失效，而本文构造的病态轨迹恰对应于(Skg)的强不适定现象（ill-posedness）：在 $d \ge 2$ 下，同一初始数据可对应无限多满足有限率的极限轨迹。这与 [46] 中讨论的激波不稳定性不同，本文病态源于对数熵耗散 $\|\nabla \log u\|_{L^2_{t,x}}$ 的维度依赖嵌入失效，而非非线性色散效应。

4. 最大熵产生原理与气体动力学Euler方程组

Maximal entropy production principle and the Euler system of gas dynamics

作者: Elisabetta Chiodaroli, Eduard Feireisl, Ondřej Kreml, Simon Markfelder

类别: math.AP

发布日期: 2026-05-26

arXiv ID: 2605.26687v1

摘要翻译: 凸积分方法揭示：即便额外施加熵不等式作为约束，气体动力学Euler方程组在弱解类中仍是不适定的。一个自然的问题随之产生：能否通过最大化熵产生率来选取物理上相关的解？首先，我们给出一个二维Riemann初值数据的例子，其标准自相似解不满足最大熵产生原理。因此，最大化熵产生率将排除该一维自相似解——而该解在此背景下直观上似乎是物理上相关的解。其次，我们证明：对一大类初值数据，存在熵容许的弱解，其总熵剖面可取任意（非减）函数。

Comment: 无

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

1. 研究对象和背景

本文研究三维（或二维）空间中全气体动力学欧拉系统（full Euler system）的弱解选择问题，其核心方程为：

$$ \partial_t \varrho + \operatorname{div}_x (\varrho \mathbf{v}) = 0, \tag{1.1} $$

$$ \partial_t (\varrho \mathbf{v}) + \operatorname{div}_x (\varrho \mathbf{v} \otimes \mathbf{v}) + \nabla_x p(\varrho, \vartheta) = 0, \tag{1.2} $$

$$ \partial_t \left( \frac{1}{2} \varrho |\mathbf{v}|^2 + \varrho e(\varrho, \vartheta) \right) + \operatorname{div}_x \left[ \left( \frac{1}{2} \varrho |\mathbf{v}|^2 + \varrho e(\varrho, \vartheta) + p(\varrho, \vartheta) \right) \mathbf{v} \right] = 0, \tag{1.3} $$

其中 $\varrho = \varrho(t,x) > 0$ 为质量密度，$\vartheta = \vartheta(t,x) > 0$ 为绝对温度，$\mathbf{v} = \mathbf{v}(t,x) \in \mathbb{R}^d$（$d=2$ 或 $3$）为速度场；压力与内能满足理想气体状态方程（Boyle–Mariotte 律）：

$$ p(\varrho, \vartheta) = \varrho \vartheta, \quad e(\varrho, \vartheta) = c_v \vartheta, \quad c_v > 0. \tag{1.4} $$

对应的比熵为 $s(\varrho, \vartheta) = c_v \log \vartheta - \log \varrho$，其满足熵不等式（以分布意义）：

$$ \partial_t (\varrho s) + \operatorname{div}_x (\varrho s \mathbf{v}) \geq 0. \tag{1.6} $$

该不等式等价于存在非负 Borel 测度 $\sigma \in L^\infty(0,T; \mathcal{M}^+(\Omega))$，使得

$$ \partial_t (\varrho s) + \operatorname{div}_x (\varrho s \mathbf{v}) = \sigma. \tag{1.6'} $$

此处 $\sigma$ 即为熵产率（entropy production rate），其总质量 $\langle \sigma(t), 1 \rangle = \frac{d^+}{dt} \int_\Omega \varrho s \, dx$ 表征单位时间系统总熵的增加量。在经典热力学中，“最大熵产率原理”（MaxEP）被经验性地用于预测远离平衡态的稳态结构（如大气环流、湍流输运），但其数学基础长期悬而未决。

对欧拉系统而言，该原理的数学实现面临根本挑战：一方面，凸积分方法已证实，即使施加熵不等式 (1.6)，Riemann 初值问题在 $d \geq 2$ 下仍存在无穷多熵可容许弱解（见 [2, 6, 16]），即系统在弱解类中不适定；另一方面，DiPerna（1985）与 Dafermos（1973, 1983）分别提出基于总熵大小与熵增长速率的偏序关系：

$$ (\varrho_1, \vartheta_1, \mathbf{v}_1) \prec_{\text{DiP}} (\varrho_2, \vartheta_2, \mathbf{v}_2) \iff \int_\Omega \varrho_1 s(\varrho_1, \vartheta_1)(t,\cdot) \, dx \leq \int_\Omega \varrho_2 s(\varrho_2, \vartheta_2)(t,\cdot) \, dx \quad \text{a.e. } t \in (0,T), \tag{1.7} $$

$$ (\varrho_1, \vartheta_1, \mathbf{v}_1) \prec_{\text{Daf}} (\varrho_2, \vartheta_2, \mathbf{v}_2) \iff \exists \tau \geq 0 \text{ s.t. } (\varrho_1, \vartheta_1, \mathbf{v}_1) = (\varrho_2, \vartheta_2, \mathbf{v}_2) \text{ on } [0,\tau), \text{ and } \frac{d^+}{dt}\Big|_{t=\tau} \int_\Omega \varrho_1 s \, dx < \frac{d^+}{dt}\Big|_{t=\tau} \int_\Omega \varrho_2 s \, dx. \tag{1.8} $$

“熵率可容许解”（entropy rate admissible solution）即指关于 $\prec_{\text{Daf}}$ 极大的弱解。然而，该原则是否真能筛选出物理上相关的解？特别是，它是否能恢复一维 Riemann 解（即自相似解）作为二维/三维情形下的标准物理解？这一问题至关重要，因为自相似解是双曲守恒律理论的基石，且在数值模拟与实验中常被默认采用。本文旨在严格检验 MaxEP 原理对全欧拉系统的适用性，聚焦两大空白：（1）在无界域 $\mathbb{R}^2$ 上，是否存在 Riemann 初值数据，使得标准一维自相似解 不 满足熵率可容许性？这将直接质疑 MaxEP 的普适性；（2）在有界域 $\Omega \subset \mathbb{R}^3$ 中，给定任意非减函数 $\tilde{S}(t)$，是否存在熵可容许弱解，使其总熵 $S(t) = \int_\Omega \varrho s \, dx$ 恰好等于 $S_0 + c_v \tilde{S}(t) \int_\Omega \varrho_0 \, dx$？若成立，则熵产率 $\sigma = c_v \tilde{S}'(t) \int_\Omega \varrho_0 \, dx \, \delta_{t}$ 可被任意指定，表明 MaxEP 无法唯一确定解——因所有可能的熵增长路径均被实现。

2. 主要定理或主要结果

本文证明了两个核心定理，其结论精确、定量且具有明确的参数依赖：

定理 2.1（一维自相似解非熵率可容许）：设 $c_v = \frac{3}{2}$。存在二维 Riemann 初值数据 $(\varrho_0, \vartheta_0, \mathbf{v}_0)$（具体为 $(\varrho_-, \mathbf{v}_-, \vartheta_-) = (1, (0,0), 2)$ 与 $(\varrho_+, \mathbf{v}_+, \vartheta_+) = (10, (0,-100), \frac{1}{10})$），使得其对应的一维自相似解 $(\varrho_s, \mathbf{v}_s, \vartheta_s)$ 在 $\mathbb{R}^2$ 上 不 是熵率可容许的。更精确地，对任意 $L > 0$，定义局部熵产率

$$ D_L[\varrho_s, \vartheta_s] := \frac{d^+}{dt} \int_{-L}^{L} \int_{-L}^{L} \varrho_s s(\varrho_s, \vartheta_s)(t, \cdot) \, dx_1 dx_2, $$

则其为常数，且满足 $-1662 < D_L[\varrho_s, \vartheta_s] / (2L) < -1661$。同时，存在另一族熵可容许弱解 $(\varrho_c, \mathbf{v}_c, \vartheta_c)$（由凸积分构造的 1-扇形解），其熵产率亦为常数，且满足 $867 < D_L[\varrho_c, \vartheta_c] / (2L) < 868$。因此 $D_L[\varrho_c, \vartheta_c] > D_L[\varrho_s, \vartheta_s]$，故自相似解非极大元，从而不满足定义 1.1 的熵率可容许性。作者进一步指出，该结论对 $c_v \in [1, \frac{3}{2}]$ 均成立（注：文中为数值验证，非严格证明）。

定理 3.1（任意熵剖面解的存在性）：设 $\Omega \subset \mathbb{R}^3$ 为有界区域，可分解为有限个互不相交的子区域 $\Omega = \bigcup_{i=1}^N Q_i$；初始密度与温度为分片常数：$\varrho_0|_{Q_i} = \varrho_i^0 > 0$, $\vartheta_0|_{Q_i} = \vartheta_i^0 > 0$；给定 $c_v = \frac{3}{2}$，任意 $\delta > 0$，及任意非减函数 $\tilde{S} \in \mathcal{S}[\delta] := \{ \tilde{S}: [0,T] \to [0,\infty) \mid \tilde{S}(t) = 0 \text{ for } t < \delta, \, \tilde{S}' \in \mathcal{M}^+([0,T]) \}$；再给定 $0 < \varepsilon < \delta$。则存在初速度 $\mathbf{v}_0 \in L^\infty(\Omega; \mathbb{R}^3)$，使得欧拉系统 (1.1)–(1.3) 在 $\Omega$ 上、配以不可渗透壁条件 $\mathbf{v} \cdot \mathbf{n}|_{\partial \Omega} = 0$ 及初值 (1.5)，存在 无穷多 熵可容许弱解 $(\varrho, \vartheta, \mathbf{v})$，满足

$$ \frac{d}{dt} \int_\Omega \varrho s \, dx = c_v \tilde{S}'(t + \varepsilon) M_0 \quad \text{in } \mathcal{D}'(0, T-\varepsilon), \quad \text{其中 } M_0 = \int_\Omega \varrho_0 \, dx. \tag{3.3} $$

特别地，若取 $\tilde{S}(t) = t$，则 $\tilde{S}' = \delta_0$，解的总熵在 $t = \varepsilon$ 处发生狄拉克型跃变；若 $\tilde{S}$ 绝对连续，则 $\frac{d}{dt} \int_\Omega \varrho s \, dx = c_v \tilde{S}'(t + \varepsilon) M_0$ 几乎处处成立。

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

本文的证明核心在于对凸积分技术的双重创新性应用，其障碍与突破均高度具体：

障碍一：在无界域上比较熵产率。标准凸积分构造（如 [2]）产生的是紧支集扰动解，其总熵在 $\mathbb{R}^2$ 上发散，无法与自相似解（其熵产率在无穷大区域内线性增长）进行有意义的比较。作者的关键创新在于引入“1-扇形解”（1-fan solution）框架。该框架要求构造一个分片常数的“子解”（subsolution），其本身满足 Rankine-Hugoniot 条件与椭圆性不等式，从而为凸积分提供稳定平台。文中（公式 2.7–2.8）显式求解了一个中间密度 $\varrho_1 = 14$ 的子解，并通过代数计算验证了所有必要不等式（[2, Proposition 4.4]）均成立。此子解的精妙之处在于其熵产率 $D_L[\varrho_c, \vartheta_c]$ 的表达式（公式 2.8）仅依赖于子解的常数状态与激波速度 $\mu_\pm$，而这些量又由初值与 $\varrho_1$ 连续决定。这使得作者能进行严格的数值围堵（$867 < D_L/(2L) < 868$），并与自相似解的围堵（$-1662 < D_L/(2L) < -1661$）形成无可辩驳的对比。此方法避开了对无限维解空间的泛泛讨论，直击要害。

障碍二：在有界域上实现任意熵剖面。标准凸积分（如 [16]）要求动能 $\frac{|\mathbf{m}|^2}{2\varrho}$ 为时间连续函数，但任意非减 $\tilde{S}(t)$ 的导数 $\tilde{S}'$ 是一个非负测度，其对应的动能剖面 $\Lambda - \frac{3}{2} \exp(\tilde{S}(t)) \varrho_i^0 \vartheta_i^0$ 通常不连续（例如在跳跃点处）。作者的真正创新点是调用 Abbatiello–Feireisl [1] 发展的“不连续凸积分”技术。该技术将动能剖面视为黎曼可积函数，其核心工具是 [1, Theorems 2.1, 2.4] 所建立的、针对非光滑时间依赖的松弛解存在性定理。证明中，作者将问题完全降维：固定密度 $\varrho(t,x) = \varrho_i^0$ 在每个 $Q_i$ 上（公式 3.4），并令温度按 $\vartheta(t,x) = \exp(\tilde{S}(t)) \vartheta_i^0$ （公式 3.8）演化，从而将熵平衡（公式 3.10）转化为一个关于动量 $\mathbf{m} = \varrho \mathbf{v}$ 的约束：$\partial_t \mathbf{m} + \operatorname{div}_x (\frac{\mathbf{m} \otimes \mathbf{m}}{\varrho_i^0} - \frac{1}{3} \frac{|\mathbf{m}|^2}{\varrho_i^0} I) = 0$，并强制其动能满足 $\frac{|\mathbf{m}|^2}{2\varrho_i^0} = \Lambda - \frac{3}{2} \exp(\tilde{S}(t)) \varrho_i^0 \vartheta_i^0$。这正是 [1] 所处理的标准形式。最后，作者利用 [15, Section 13.6] 的技巧，将所得的一个解作为新的子解，启动第二轮凸积分，从而生成无穷多解。整个框架的严谨性建立在对 [1] 和 [15] 定理的精确引用与适配之上，而非空泛的“使用凸积分”。

统一障碍：物理量的精确控制。无论是定理 2.1 的数值围堵，还是定理 3.1 中对 $\tilde{S}(t+\varepsilon)$ 的时移，都体现了作者对物理量（熵、动能、压力）在构造中如何被精确操控的深刻理解。他们没有回避计算，而是将计算本身作为证明的核心环节和可靠性的来源。

4. 与之前工作的比较

本文与既有工作的比较需从三个维度展开：

与 DiPerna/Dafermos 原始工作比较：DiPerna [14] 与 Dafermos [8] 提出的 $\prec_{\text{DiP}}$ 与 $\prec_{\text{Daf}}$ 关系，其原始动机是为封闭系统（能量守恒边界条件）提供一个物理上合理的解选择准则。本文首次在 开放系统（无界域 $\mathbb{R}^2$）中对 $\prec_{\text{Daf}}$ 进行了严格检验，并给出了首个反例（定理 2.1）。关键差异在于，本文为克服无界域上总熵发散的问题，对 $\prec_{\text{Daf}}$ 进行了本质性修改（定义 1.1），即在有限正方形 $[-L,L]^2$ 上比较熵增长率。这一修改并非技术妥协，而是为了使物理量（熵产率）具有可比性。相比之下，[7] 对等熵欧拉系统的研究虽也得到类似反例，但其模型（$p = \varrho^\gamma$）与本文的全欧拉系统（含温度、熵作为独立变量）在热力学结构上截然不同，因此本文结果是前者的实质性非等熵推广。

与凸积分文献比较：本文与 [2, 6, 16] 等凸积分工作构成互补而非重复。[2] 证明了 Riemann 问题的非唯一性，但未涉及任何选择准则；[16] 构造了具有“野生”初值的解，但未分析其熵产特性。本文的推进在于，它将凸积分从“存在性工具”提升为“选择性分析工具”。定理 2.1 利用凸积分构造的解 主动击败 自相似解；定理 3.1 则利用凸积分 实现任意熵产目标。这揭示了一个深刻结论：凸积分解不仅数量众多，而且其熵产率在某种意义上是“稠密”的，足以覆盖所有可能的非减剖面。这与 [18] 的结果形成张力：[18] 证明熵率可容许的 DMV 解必为弱解，暗示该准则可能过于苛刻；而本文定理 3.1 则表明，即便在弱解类中，熵率可容许解也可能极度稀少（甚至不存在），因为几乎所有可能的熵增长路径都被其他弱解占据。

与能量耗散率准则比较：[7, 24] 已证明，对等熵欧拉系统，“最大能量耗散率”准则同样无法选出自相似解。本文定理 2.1 将此结论移植到更真实的全欧拉系统，确认了该现象的鲁棒性。然而，二者代价不同：[7] 的反例依赖于等熵假设（$s = \text{const}$），其“能量”在物理上实为熵；而本文的反例必须同时处理密度、温度、速度三个变量及其耦合，其构造（如公式 2.2–2.5 对中间压强 $p_M$ 的求解）远为复杂。这表明，MaxEP 的失效并非源于模型简化，而是深植于全欧拉系统的多维非线性本质。作者在结论部分（第4节）明确指出，结合 [18] 与本文结果，可以合理猜想：对于“一般”的初值数据，熵率可容许的弱解 可能根本不存在。这是一个比单纯“准则无效”更为深刻的负面结论，它从根本上质疑了将 MaxEP 作为欧拉系统基础物理原理的可行性。

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

1. 含时依赖势的薛定谔方程的正则性传播

Propagation of Regularity for Schroedinger Equations with Time Dependent Potentials

作者: Avy Soffer

类别: math.AP

发布日期: 2026-05-26

arXiv ID: 2605.27321v1

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

摘要翻译

中文标题: 含时依赖势的薛定谔方程的正则性传播

摘要翻译: 考虑具有任意时间依赖势的薛定谔方程的动力学行为。证明了：对空间局部化的势，存在正则性的传播现象，且该传播在高阶 Sobolev 范数下具有一致有界性。与解发生散射的情形不同——此时正则性传播通常可通过标准的 bootstrap 论证来证明——本文所考虑的解包含一个不散射的部分，这在一般情形下是预期的。为此，我们引入了直接作用于（例如）$H^2(\mathcal{R}^3)$ 中的传播估计。

2. 二维 Zakharov-Kuznetsov 方程的谱性质

Spectral property for the 2D Zakharov-Kuznetsov equation

作者: Justin Holmer, Svetlana Roudenko

类别: math.AP, math-ph

发布日期: 2026-05-26

arXiv ID: 2605.27341v1

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

摘要翻译

中文标题: 二维 Zakharov-Kuznetsov 方程的谱性质

摘要翻译: 我们讨论二维 Zakharov-Kuznetsov (ZK) 方程的 virial 算子的谱性质。这是在高维问题中建立孤立波爆破（blow-up）或渐近稳定性（asymptotic stability）的关键要素。 该模型在三维情形下最初由 Zakharov 和 Kuznetsov 在等离子体物理中提出，同时也是著名 Korteweg-de Vries (KdV) 方程的高维推广。ZK 方程中孤立波的稳定性问题，或修正 ZK（或 KdV 型）方程中稳定爆破问题，是一个重要的物理问题；其中 virial 算子及其谱性质是分析中的核心要素。本文从解析角度研究该问题，并将其归约为仅需数值验证若干内积的符号及某些特征值。

3. Kenig 关于带零边界条件的 Sobolev 空间的插值问题的一个反例

A Counterexample to Kenig’s Interpolation Problem for Sobolev Spaces with Zero Boundary Conditions

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

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

发布日期: 2026-05-26

arXiv ID: 2605.27119v1

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

摘要翻译

中文标题: Kenig 关于带零边界条件的 Sobolev 空间的插值问题的一个反例

摘要翻译: 设 $n\in \mathbb N\cap[2,\infty)$。本文证明：存在一个有界 $C^1$ 区域 $Ω\subset \mathbb R^n$，使得对任意给定的 $s\in(1,2)\setminus\{\frac32\}$，\begin{align*} \left[H_0^1(Ω),H^2(Ω)\cap H_0^1(Ω)\right]{s-1} =H^s(Ω)\cap H_0^1(Ω)=H_0^s(Ω) \end{align*} 成立（范数等价），但\begin{align*} \left[H_0^1(Ω),H^2(Ω)\cap H_0^1(Ω)\right]{\frac12} \subsetneqq H^{\frac32}(Ω)\cap H_0^1(Ω), \end{align*} 从而为 Kenig 在 [CBMS Regional Conf. Ser. in Math. 83, 1994] 中提出的第 3.3.19 号问题提供了一个反例。作为应用，我们证明：对这样的区域 $Ω$，\begin{align*} H^2(Ω)\cap H_0^1(Ω)\subsetneqq D(-Δ_D) \end{align*}（其中 $D(-Δ_D)$ 表示区域 $Ω$ 上 Dirichlet Laplacian 算子 $-Δ_D$ 的定义域），并构造了一个齐次热方程在零 Dirichlet 边界条件下的解，该解不属于 $L^2((0,T);H^2(Ω)\cap H_0^1(Ω))$，对任意给定的 $T\in(0,\infty)$ 均成立。

4. 一个模型算子的 Gevrey 可解性的缺失

Lack of Gevrey solvability for a model operator

作者: Enrico Bernardi

类别: math.AP

发布日期: 2026-05-26

arXiv ID: 2605.27054v1

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

摘要翻译

中文标题: 一个模型算子的 Gevrey 可解性的缺失

摘要翻译: 我们证明：在 $ \R^{4} $ 中，模型双曲型算子

$$ Q=-D_t^2+2xD_tD_y+D_x^2+x^3D_y^2+D_z^2+z^2D_y^2 $$

的柯西问题在原点处不具有局部可解性，当 Gevrey 阶 $ s>6 $ 时。

5. 具有非线性 Brinkman 定律的非局部 Busenberg-Travis 交叉扩散系统

A nonlocal Busenberg-Travis cross-diffusion system with nonlinear Brinkman law

作者: Peter Hirvonen, Ansgar Jüngel

类别: math.AP

发布日期: 2026-05-26

arXiv ID: 2605.27037v1

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

摘要翻译

中文标题: 具有非线性 Brinkman 定律的非局部 Busenberg-Travis 交叉扩散系统

摘要翻译: 在具无通量边界条件的有界区域中，分析了一类用于种群分离的非局部 Busenberg-Travis 交叉扩散系统。各物种的速度满足一个正则化的 Darcy 定律，该定律可解释为 Brinkman 方程。与文献中的结果相比，本文假设密度-压力关系是非线性的。对幂律非线性项的指数在较宽范围内，证明了该系统的弱解的整体存在性，并证明了局域化极限。证明基于来自 Tsallis 熵不等式的统一估计。由于正则性问题，原始问题通过多种逼近格式进行正则化处理，并借助紧性论证获得去正则化极限。

6. 球面平均变换的统一值域刻画

A Unified Range Characterization for the Spherical mean transform

作者: Pradipta Chatterjee, Nisha Singhal, Abhilash Tushir

类别: math.AP, math.CA

发布日期: 2026-05-26

arXiv ID: 2605.27034v1

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

摘要翻译

中文标题: 球面平均变换的统一值域刻画

摘要翻译: 本文研究支撑在单位球内的函数的球面平均变换（SMT）的值域刻画。在先前的工作中，奇数维情形下得到了一组微分条件，而偶数维情形下则得到了积分条件。我们证明：这些依赖于维数奇偶性的不同条件，在奇数维情形下实际上是等价的。这一等价性表明，积分条件给出了一个统一且简洁的SMT值域刻画，该刻画对偶数维与奇数维均成立。

7. 耦合超抛物型算子的向后唯一性及其在加加速度驱动控制模型中的应用

Backward Uniqueness for Coupled Ultraparabolic Operators and an Application to Jerk-Driven Control Models

作者: Xiao-Dong Cao, Chao-Jiang Xu, Yan Xu

类别: math.AP

发布日期: 2026-05-26

arXiv ID: 2605.26993v1

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

摘要翻译

中文标题: 耦合超抛物型算子的向后唯一性及其在加加速度驱动控制模型中的应用

摘要翻译: 我们证明了一类具有耦合线性漂移项的超抛物型算子的向后唯一性。主要困难在于：退化变量上的傅里叶变换将耦合漂移项转化为对偶频率变量上的输运算子，因此经典的 Littlewood–Paley Carleman 方法不能直接适用。我们通过引入一个不变频率变量，并建立一种适配于该输运结构的频率局部化 Carleman 估计来克服这一困难。该结果对 W. Wang 和 L. Zhang $\left[ \emph {Methods Appl. Anal.}, \ 20 \ (1) \ (2013) \ 79-88 \right]$ 提出的问题——在常系数耦合漂移、扩散项及低阶系数依赖于时间和扩散变量的情形下——给出了部分回答。作为应用，在加加速度（jerk）驱动的控制模型中，我们证明了描述位置、速度、加速度或加加速度误差的方程的向后唯一性：在低阶系数有界条件下，若最终时刻的误差在 $L^2$ 中为零，则所有更早时刻的误差亦恒为零。

8. 混合 Ermakov-Ray-Reid / Painlevé II 对称性约化：应用于一类移动边界问题

Hybrid Ermakov-Ray-Reid/Painlevé II Symmetry Reduction: Application to a Class of Moving Boundary Problems

作者: Colin Rogers, Adriana C. Briozzo

类别: math.AP, nlin.SI

发布日期: 2026-05-26

arXiv ID: 2605.26939v1

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

摘要翻译

中文标题: 混合 Ermakov-Ray-Reid / Painlevé II 对称性约化：应用于一类移动边界问题

摘要翻译: 本文表明，一类针对两分量 mKdV 系统之新颖推广的非线性移动边界问题，可通过应用混合 Ermakov-Ray-Reid / Painlevé II 对称性试探解（ansatz）而获得精确解。该 mKdV 系统源于一个耦合非线性 NLS 系统的约化，该 NLS 系统包含 deBroglie–Bohm 势项。

9. $C^*$-代数与带移位算子的传输问题及椭圆边值问题

$C^*$-algebras of transmission problems and elliptic boundary value problems with shift operators

作者: Alexandre Baldare, Vladimir E. Nazaikinskii, Anton Yu. Savin, Elmar Schrohe

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

发布日期: 2026-05-26

arXiv ID: 2605.26593v1

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

摘要翻译

中文标题: $C^*$-代数与带移位算子的传输问题及椭圆边值问题

摘要翻译: 我们研究一类与光滑流形上的群作用相关的新颖非局部边值问题的Fredholm可解性。具体而言，我们考虑群作用定义在无边界的环境流形上，且不保持问题所陈述的带边界流形的情形；特别地，该群作用并不将边界映射到自身。边界在群作用下的轨道将流形分割为若干子区域，这一分解与$C^*$-代数技术相结合，在我们对该问题的分析方法中起着关键作用。

10. 浅水波局部模型中周期行波的调制稳定性

Modulational stability of the periodic traveling wave in a local model for shallow water waves

作者: Lili Fan, Xin Zhang, Hongjun Gao

类别: math.AP

发布日期: 2026-05-26

arXiv ID: 2605.26573v1

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

摘要翻译

中文标题: 浅水波局部模型中周期行波的调制稳定性

摘要翻译: 本文研究浅水波局部模型中周期行波的调制稳定性，该模型是 Hunter-Saxton 方程的一个推广版本。我们为此局部模型构造一族小振幅周期行波，并给出这些波的参数化表示。借助 Floquet-Bloch 理论、摄动理论与谱分析，我们建立了这些背景周期行波解的调制稳定性。最后，我们采用类似方法分析另一类具有三次非线性的推广 Hunter-Saxton 方程的调制不稳定性。

11. 闵可夫斯基时空的稳定性

The Stability of Minkowski Spacetime

作者: Dawei Shen

类别: gr-qc, math.AP, math.DG

发布日期: 2026-05-26

arXiv ID: 2605.26506v1

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

摘要翻译

中文标题: 闵可夫斯基时空的稳定性

摘要翻译: 闵可夫斯基时空的非线性稳定性是广义相对论数学理论中的核心成就之一，更广泛而言，亦是几何非线性波动方程分析领域的一项中心成果。自Christodoulou–Klainerman的开创性工作以来，该问题深刻推动了我们对衰减、色散以及爱因斯坦真空方程中几何与分析之间复杂相互作用的理解。 本综述概述了闵可夫斯基时空稳定性理论所依赖的主要思想与技术。我们着重强调衰减假设、几何叶状结构（geometric foliations）、能量恒等式（energy identities）以及规范选取（gauge choices）在全局分析中的作用。特别关注外部区域的稳定性结果、最小衰减机制（minimal decay regimes）以及临界情形（borderline case）：在此情形下，非线性相互作用的时空可积性失效揭示出微妙的阈值现象（threshold phenomena）。 我们的目标是为该领域的发展提供一个连贯的视角，厘清已知结果背后所依赖的结构性机制，并概述目前在临界情形下仍悬而未决的一些核心公开问题。

12. 调和函数零点集的无界拓扑

Unbounded Topology of Nodal Sets of Harmonic Functions

作者: Robert Koirala

类别: math.AP, math.DG, math.GT

发布日期: 2026-05-25

arXiv ID: 2605.26274v1

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

摘要翻译

中文标题: 调和函数零点集的无界拓扑

摘要翻译: 对于每个整数 $n\ge 3$、每个 $1\le \ell\le n-2$ 以及每个足够大的整数 $m$，我们在单位球 $B_1(0)\subset\mathbb{R}^n$ 上构造调和函数 $u_{m,\ell}$，使得其频率关于 $m$ 一致有界，且零点集 $\{u_{m,\ell}=0\}\cap B_{1/2}(0)$ 中每一点均为正则点，但其贝蒂数满足 \begin{align*} b_\ell\bigl({u_{m,\ell}=0}\cap B_{1/2}(0)\bigr)\ge 2m. \end{align*} 因此，即使频率有界且零点集正则，亦不能推出拓扑结构的一致上界。特别地，这些例子构成了 Lin 与 Liu 所声称的全局贝蒂数上界的反例。

13. 不确定性原理、谱局部化与紧流形上的奇异薛定谔算子

Uncertainty Principles, Spectral Localization, and Singular Schrödinger Operators on Compact Manifolds

作者: Alex Iosevich, Chamsol Park

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

发布日期: 2026-05-25

arXiv ID: 2605.26264v1

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

摘要翻译

中文标题: 不确定性原理、谱局部化与紧流形上的奇异薛定谔算子

摘要翻译: 我们通过将正交系的限制估计与拉普拉斯–贝尔特拉米算子及薛定谔算子的谱投影估计相结合，在无边紧黎曼流形上建立了不确定性原理。我们的结果将谱局部化函数之支集大小与底层谱簇的基数以及傅里叶比型量联系起来。我们获得了适用于属于Kato类与尺度临界类的奇异位势薛定谔算子的相应类比结果。作为应用，我们证明了在紧流形上由不完全谱数据进行唯一恢复的结果。在曲率假设下（包括非正与负截面曲率），我们还证明了与收缩谱窗口相关的对数改进型不确定性原理。

不相关论文（基本信息）

1. 二次正则化最优传输的Polyak-Lojasiewicz不等式

Polyak-Lojasiewicz Inequality for Quadratically Regularized Optimal Transport

作者: Alberto González-Sanz, Marcel Nutz, Andrés Riveros Valdevenito

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

发布日期: 2026-05-26

arXiv ID: 2605.27175v1

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

摘要: Quadratically regularized optimal transport (QOT) is an alternative to entropic regularization that yields sparse couplings and avoids numerical instabilities due to exponential scaling. From an optimization viewpoint, the dual QOT objective is concave but features a positive part function which prevents strong concavity and reduces smoothness of optimizers. Consequently, standard arguments for linear convergence of algorithms do not apply. In this paper, we nevertheless establish a quantitative curvature property for the QOT dual. Under mild assumptions covering both continuous and semi-discrete transport problems, we prove a local error bound and a Polyak-Lojasiewicz (PL) inequality, with explicit constants depending only on the problem primitives. These results are obtained by functional-analytic techniques exploiting that near the optimum, the argument of the positive part function is positive on the interior of the support of the optimal coupling. As applications, we derive linear convergence of the gradient ascent, coordinate ascent, and coordinate gradient ascent algorithms on the dual problem, with explicit contraction rates.

2. Wasserstein–Ebin 度量：非平衡最优传输到黎曼度量空间的几何提升

The Wasserstein–Ebin Metric: A Geometric Lift of Unbalanced Optimal Transport to the space of Riemannian metrics

作者: Martin Bauer, Peter W. Michor, François-Xavier Vialard

类别: math.DG, math.AP

发布日期: 2026-05-26

arXiv ID: 2605.27086v1

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

摘要: We introduce dynamic and static formulations that formally extend unbalanced optimal transport from the space of positive densities to the space of Riemannian metrics. The first construction is based on a dynamic variational formulation in which the evolution of a Riemannian metric is driven by transport together with a source term. Choosing the $L^2$-metric to penalize the transport vector field and the Ebin metric to penalize the source component yields a new Riemannian metric on the manifold of Riemannian metrics, which we call the Wasserstein–Ebin metric. Our main result shows that the volume map defines a Riemannian submersion from the Wasserstein–Ebin metric to the Wasserstein–Fisher–Rao metric on the space of smooth densities. In addition, we construct a Riemannian submersion from the automorphism group of the tangent bundle onto the space of Riemannian metrics, providing a generalization of Otto’s geometric description for the Wasserstein metric to the setting of the Wasserstein–Ebin metric. To propose a static formulation of unbalanced optimal Riemannian metric transport, we introduce two Kullback–Leibler-type divergences on the space of Riemannian metrics: one inspired by matrix information geometry, and another related, through the volume map, to the classical Kullback–Leibler divergence on densities. Establishing a link between the static and dynamic formulations remains an open direction for future work.