🔢 mathematics

Non-Uniqueness for Nonlinear Fokker--Planck Equations and Their Associated Distribution-Dependent SDEs

本文通过构造具有临界正则性的无散漂移项，从而证明了在 $\mathbb T^d$ 或 $\mathbb R^d$ ( $d\geq 2$ ) 上，依赖于分布的随机微分方程及其相关的非线性福克-普朗克方程存在平稳解非唯一性，且这些漂移项会产生无限多个或任意多个不同的平衡态。

原作者： Huaxiang Lü

发布于 2026-07-01

📖 1 分钟阅读🧠 深度阅读

CC BY 4.0

原作者： Huaxiang Lü

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一下，你正在观察一群在城市中移动的、隐形的巨型人群。在物理学和数学中，我们经常使用方程来预测这群人的行为。通常情况下，我们假设如果我们知道每个人从哪里开始，以及风（或外部力量）是如何吹拂的，那么人群就会以一种特定的、可预测的方式移动。这被称为“适定性”（well-posedness）。

由吕华祥（Huaxiang Lü）撰写的这篇论文挑战了这一假设。研究表明，对于某些类型的群体（具体而言，是那些每个人的运动都取决于整个群体密度的群体），预测规则可能会失效。即使你从一个完美的、静止的群体开始，并施加一个微小的、几乎察觉不到的推力，人群也可能会分裂成许多种不同的可能未来，而不是遵循单一路径。

以下是使用简单类比对该论文主要思想的拆解：

1. 设置：人群与风

该论文研究了一种被称为**分布依赖随机微分方程（DDSDE）**的数学模型。

类比： 想象一群人在公园里散步。每个人的运动都受到两件事的影响：
1. 风：一个推动他们的外部力量（比如微风）。
2. 人群： 他们周围有多拥挤。如果很拥挤，他们的移动方式会与空旷时不同。
方程： 这是“非线性福克-普朗克方程”（Nonlinear Fokker-Planck Equation）。这是一个复杂的配方，它告诉我们人群的密度随时间如何变化。

2. 问题：“光滑”与“粗糙”

在现实世界中，风通常是平滑的。在数学中，我们更喜欢“光滑”的力量，因为它们使方程易于求解，并且能得到唯一的答案（一个特定的未来）。

然而，作者提出了一个问题：如果风是“粗糙”的呢？

隐喻： 想象风不是一阵温柔的微风，而是一阵混乱、锯齿状的狂风，瞬间且反复无常地改变方向。在数学术语中，这是“低正则性”或“粗糙”的漂移。
发现： 论文证明，如果这种粗糙的风仅仅是稍微有一点点“锯齿感”（具体来说是在数学上的 $L^d$ 类别中），系统就会失去预测单一未来的能力。

3. 主要结果：一个起点，多种结局

论文最大的主张是非唯一性（Non-Uniqueness）。

场景： 想象人群正静止不动（一个“平稳态”）。风很平静。
转折： 作者引入了一种微小的、隐形的、粗糙的风，它如此之小，以至于你几乎无法测量它。
结果： 数学表明，从这个完全相同的起点出发，无数种不同的情景都可能发生。
- 在一种情景下，人群向左旋转。
- 在另一种情景下，他们向右旋转。
- 在第三种情景下，他们分裂成两组。
- 所有这些都是该方程在相同初始条件下合法的数学解。

4. 方法：“凸集成”（建筑师的技巧）

作者是如何证明这一点的？他们使用了一种名为**凸集成（Convex Integration）**的技术。

类比： 想象你是一位建筑师，试图建造一座从远处看像直线、但近看其实是由成千上万个微小的、锯齿状折线组成的塔楼。
过程：
1. 从一个简单的、光滑的解开始（直线）。
2. 向风中加入一个微小的、高频的“抖动”（粗糙扰动）。
3. 这个抖动会在预测中产生一个误差。
4. 作者随后加入另一个抖动来抵消这个误差，但这个新抖动又会产生一个新的微小误差。
5. 他们无数次地重复这个过程。
魔力： 通过仔细控制这些抖动的大小和频率，他们可以使最终结果看起来像是一个有效的解，但由于他们添加了如此多不同的“抖动”，他们可以从同一个起点构建出许多不同的有效解。

5. 两项具体发现

论文展示了两类主要的混沌现象：

发现 A：从平稳态出发的无限演化
即使人群开始时完全静止，微小的粗糙风也会导致系统演化成无数种不同的模式。这发生在“临界阈值”处——这意味着风足够粗糙以至于打破了规则，但又没粗糙到让数学完全崩溃。
发现 B：多重稳态（Multistability）
论文还表明，一个系统可以拥有多种不同的“静止状态”。
- 类比： 想象一个景观中的球。通常，球会滚向单个山谷的底部（一个稳定状态）。
- 结果： 作者展示了通过适当的粗糙风，你可以创造出一个拥有许多山谷的景观。球可以落在其中任何一个山谷中，而系统无法“知道”它应该选择哪一个。这模拟了物理学中的相变（例如水变成冰或蒸汽），即系统可以存在于多种稳定形式之中。

6. 应用领域

作者在两种类型的“城市”上测试了这一点：

环面（Torus，甜甜圈形状）： 一个有限的、循环的空间（就像一个视频游戏世界，你从左侧走出，会从右侧出现）。
全空间（The Whole Space，无限平面）： 一个开放的、无限的世界。

在这两种情况下，结果都成立：微小的、粗糙的扰动破坏了解的唯一性。

总结

简单来说，这篇论文说：“如果你有一群相互作用的粒子，并用一种非常特殊的、粗糙且锯齿状的力量去推动它们，你就无法预测它们的未来。即使它们都从完全相同的地方出发，它们也可能走向无数个不同的方向，或者沉淀为许多不同的静止位置。”

这是关于复杂系统预测极限的一项基础性发现，它表明“光滑性”往往是维持我们的数学模型不崩塌为混沌的唯一保障。

技术摘要：非线性 Fokker–Planck 方程及其相关依赖分布的 SDE 的非唯一性

问题陈述
本文研究了依赖分布的随机微分方程（DDSDEs）及其相关的非线性 Fokker–Planck 方程（nFPEs）的适定性问题。具体而言，本文探讨了非唯一性现象，即在这些系统中，当漂移项（drift）具有某种分布时，解的测度（law）是否存在非唯一性。虽然经典理论已确立了对于具有足够正则系数（例如漂移项属于 $L^r_t L^p_x$ 且满足 $\frac{d}{p} + \frac{2}{p} < 1$ ）的 DDSDEs 的适定性，但本研究探索了漂移项较为粗糙的临界（critical）与超临界（supercritical）情形。核心问题在于：一个任意小且足够粗糙的无散（divergence-free）扰动漂移场，是否会破坏从平稳态开始的时间演化过程以及平稳态本身的唯一性。

研究考虑了两种主要的设定：

周期性设定： 定义域为环面 $\mathbb{T}^d$ （ $d \ge 2$ ）。
全空间设定： 定义域为 $\mathbb{R}^d$ （ $d \ge 2$ ），侧重于具有限制势（confining potentials）和非线性扩散的方程。

所考虑的方程具有如下一般形式：
$dX_t = v(t, X_t) dt + (V * \rho_t)(X_t) dt + \sqrt{2}\sigma(t, X_t, \rho_t) dW_t,$
其中 $\rho_t$ 是 $X_t$ 的概率密度， $v$ 是外部漂移， $V$ 代表相互作用核， $\sigma$ 是依赖于分布的扩散系数。相关的 nFPE 为：
$\partial_t \rho - \text{div}\text{div}((\sigma\sigma^T)(\rho)\rho) + \text{div}(v\rho) + \text{div}((V*\rho)\rho) = 0.$

方法论
主要采用的工具是凸集成法（convex integration method），该方法最初由流体力学（如 Euler 和 Navier–Stokes 方程）的发展而产生，用于构造“野蛮解”（wild solutions）。作者将该技术应用于非线性、依赖分布的抛物型方程。

关键的技术创新与改进包括：

处理非线性： 不同于线性 Fokker–Planck 方程中解的凸组合仍为解，nFPE 是真正的非线性方程。作者无法简单地叠加解，而是通过迭代消除耦合方程组中的应力项（stress terms），从而构造出无穷多不同的解。
无穷系统构造： 为了生成无穷多个解（定理 1.3），作者构造了一个递增的整数序列 $\{N_q\}$ ，并在第 $q+1$ 步仅对前 $N_{q+1}$ 个方程应用凸集成方案。他们引入了一个衰减因子 $2^{-i}$ 来处理索引 $i > N_{q+1}$ 的应力项，以确保总扰动收敛，同时保持活跃方程所需的估计。
构建模块： 构造过程利用了特定的“间歇性时空喷流”（intermittent spatial-time jets，针对时间相关情况）和“Mikado 型流体”（Mikado-type flows，针对平稳情况）。这些构建模块旨在具有不相交的支持集和特定的振荡性质，以抵消误差项（振荡、非线性、线性及交换子误差）并控制漂移的正则性。
正则性阈值： 构造目标是属于 $C_t L^{d-}$ 空间的漂移项（在时间上连续，在空间上属于 $L^p$ 其中 $p < d$ ），这处于自然的临界正则性阈值。论文证明了在低正则性拓扑下，任何平滑漂移项的邻域内都存在非唯一性。
全空间适配： 对于 $\mathbb{R}^d$ 的情况，作者利用 Bogovskii 算子 来处理逆散度问题，同时保持紧支撑特性。这使得构造可以被局部化到一个固定区域，确保扰动不会影响系统在无穷远处的行为或限制势。

主要贡献与结果

时间演化的非唯一性（定理 1.3 & 1.6）：
- 作者证明了对于 $d \ge 2$ ，存在无散漂移场 $v \in C_t L^{d-}$ （在低正则性范数下任意接近平滑漂移），使得 nFPE 从同一个平稳初始密度（例如环面上的均匀测度或 $\mathbb{R}^d$ 中的特定平稳态）开始，具有无穷多个不同的概率密度解。
- 根据叠加原理，这意味着对应的 DDSDE 存在无穷多个不同的鞅解（martingale solutions）。
- 该结果被证明是本质上尖锐的，因为对于 $C_t L^{d+}$ 中的漂移项，预计是适定的。
平稳态的非唯一性（定理 1.5 & 1.7）：
- 对于 $d \ge 3$ ，论文构造了无散漂移场，使得平稳 nFPE 对于任何给定的 $N \in \mathbb{N}$ ，至少拥有 $N$ 个不同的平稳概率密度解。
- 这种平衡态的多样性类似于物理系统中的多稳态（multistability）和相变现象。
- 该结果适用于正则相互作用核以及更奇异的情况（例如 $V \sim |x|^{1-d}$ 的格林函数相互作用）。
向奇异相互作用的扩展（定理 1.4）：
- 该非唯一性机制被扩展到具有更奇异相互作用核的情况（例如 $V = \text{div } b$ ，其中 $b \in L^{\frac{d}{d-\gamma}, \infty}$ ）。
- 在此情形下，作者获得了具有改进的时间正则性估计（ $C_t$ 边界）的 $N$ 个不同解，以证明涉及奇异核的非线性相互作用项的合理性。
特定物理模型：
该框架被证明适用于广泛的类别，包括：
- 具有粗糙漂移的独立于分布的 SDE。
- 卷积型 McKean–Vlasov 方程。
- 类似孔隙介质（porous-medium）类型的方程。
- Keller–Segel 方程（趋化作用）。
- 二维涡量 Navier–Stokes 方程。
- Landau 型方程。

意义与主张
本文声称建立了第一类针对广泛的非线性且依赖分布的 Fokker–Planck 方程的非唯一性结果。作者强调，构造无穷多个解的过程与线性情况有着本质区别，因为在非线性设定下，解的凸组合通常不再是解。

研究结果表明，即使在漂移项是无散且在低正则性范数下任意小时，适定性也是脆弱的。这些发现表明，这些系统中的唯一性临界正则性阈值是尖锐的。此外，构造多个平稳态为具有粗糙外部场的平均场系统中的多稳态提供了数学机制，将对相变的理解扩展到了超越经典常数扩散模型的领域。

这项工作并非提出新的实验应用，而是为一类广泛的相互作用粒子系统极限中的唯一性猜想提供了严谨的理论反例。