生成时间: 2026年06月05日 06:53:30
运行状态: 完成
计划处理论文数量: 21
已完成论文数量: 21
未完成论文数量: 0
重点关注论文数量: 2
了解领域论文数量: 12
不相关论文数量: 7
重点关注论文(完整分析)
1. 关于在全局吸引子上演化之非线性动力系统的统计推断
On statistical inference for non-linear dynamical systems evolving in their global attractor
作者: Dimitri Konen, Richard Nickl
类别: math.ST, math.AP, math.DS
发布日期: 2026-06-04
arXiv ID: 2606.06018v1
摘要翻译: 我们在反应函数满足自然条件且初始条件为 $θ$ 的情形下,考虑一个二维周期反应-扩散系统。我们证明:在所得动力系统 $(u_θ(t):t>0)$ 的全局吸引子 $\mathcal A$ 上,反向 Poincaré 不等式成立;从而对任意固定的 $t>0$,映射 $θ\mapsto u_θ(t)$ 在 $\mathcal A$ 上满足 $L^2$-Lipschitz 稳定性估计。随后我们证明:从系统在‘快速’离散时间点上的观测数据出发,可在参数收敛速率附近实现对吸引子 $\mathcal A$ 中初始条件 $θ$ 的统计重构,以及对状态 $u_θ$ 的预测。
Comment: 无
链接: http://arxiv.org/abs/2606.06018v1
1. 研究对象和背景
本文研究二维环面 $\Omega = (0,1]^2$ 上的周期反应-扩散方程
其中 $f \in C^3(\mathbb{R})$ 满足双侧增长与有界导数条件:存在 $p > 2$、$\alpha_1 > \alpha_2 > 0$、$k \geq 0$、$L \in \mathbb{R}$,使得
该模型属于典型的耗散非线性抛物型 PDE,其长期行为由全局吸引子 $\mathcal{A} \subset L^2(\Omega)$ 刻画,定义为
其中 $S(t)\theta = u_\theta(t)$ 是对应于 (4) 的解半群,$B \subset L^2(\Omega)$ 是任一包含所有有界初值轨道最终进入的紧吸收集。由 Babin & Vishik (1992)、Temam (1997)、Robinson (2001) 等经典理论可知,在 (5) 下,系统存在唯一紧不变全局吸引子 $\mathcal{A}$,且 $\mathcal{A}$ 有界于 $H^2(\Omega)$,具有有限正 Hausdorff 维数 $D < \infty$;进一步,$\mathcal{A}$ 被包含于一个有限维惯性流形 $\mathcal{M} = \{v + \Phi(v) : v \in P_n[L^2(\Omega)]\}$ 中,其中 $P_n$ 是 Laplacian $-\Delta$ 前 $n$ 个特征函数 $e_j$(对应特征值 $\lambda_j \asymp j$)张成的投影,$\Phi : P_n[L^2] \to Q_n[L^2]$ 是 Lipschitz 映射,满足 $\|\Phi(v)-\Phi(w)\|_{H^2} \leq C \|v-w\|_{L^2}$。问题的核心在于:当系统已进入平衡态——即初值 $\theta$ 已知属于 $\mathcal{A}$,且观测数据仅来自固定时间 $t > 0$ 的离散采样 $Y_i = u_\theta(t, X_i) + \varepsilon_i$($X_i \sim \mathrm{Unif}(\Omega)$,$\varepsilon_i \sim \mathcal{N}(0,1)$ 独立同分布)——此时经典的初始条件反演问题 $\theta \mapsto u_\theta(t)$ 极度不适定。若不限制 $\theta$ 在 $\mathcal{A}$ 内,而仅假设 $\theta \in H^\alpha$(如 Nickl & Titi (2024) 对 2D 周期 Navier-Stokes 所证),则最优极小极大收敛速率为 $1/\log N$,属对数慢速;即便对热方程扰动情形,此速率亦不可改进。本文旨在填补的关键空白是:在 $\theta \in \mathcal{A}$ 这一物理上自然的“平衡态先验”下,能否突破对数速率瓶颈?这不仅关乎统计推断效率,更触及非线性动力系统长期行为的可识别性本质:全局吸引子虽为无限维空间 $L^2(\Omega)$ 的子集,但其几何结构(有限维覆盖数、惯性流形参数化)是否足以支撑近参数速率?本文将严格证明:正是 $\mathcal{A}$ 的内在几何——特别是其上成立的反向 Poincaré 不等式——彻底改变了不适定性格局,使统计恢复从 $1/\log N$ 加速至 $\sqrt{\log N / N}$。
2. 主要定理或主要结果
定理 1.1(统计预测速率):设 $\theta_0 \in \mathcal{A}$,数据 $(Y_i, X_i)_{i=1}^N$ 依模型 (3) 生成,其中 $u_\theta$ 满足 PDE (4) 且 $f$ 满足 (5)。则存在估计量 $\hat{\theta}_N = \hat{\theta}_N((Y_i,X_i)_{i=1}^N)$ 取值于 $\mathcal{A}$,使得对任意 $\tau > 0$,有
定理 1.2 / 2.4(Lipschitz 稳定性):对任意 $t \geq 0$,存在仅依赖于 $f$ 的常数 $c = c(f) > 0$,使得
该估计的精确常数见定理 2.4,且结论对 $\theta, \vartheta \in \mathcal{M}$(惯性流形)亦成立。定理 2.1 与 2.2(后验收缩):对 $\delta_N = \sqrt{\log N / N}$,分别采用 $\varepsilon$-网先验 $\Pi = \sum_{j\geq1} c_j \Pi_{\varepsilon_j}$(其中 $\varepsilon_j = j^{-\tau}$,$c_j \propto j^{-(1+\tau)}$,$\Pi_{\varepsilon_j}$ 在 $\mathcal{A}$ 的 $\varepsilon_j$-覆盖上均匀分布)或惯性流形先验 $\Pi = \mathrm{Law}(F(X))$($X$ 在 $n$-维椭球 $D(\Phi) \subset P_n[L^2]$ 上具正密度),对应的后验分布满足
定理 2.5(初值估计):上述任一先验下,后验分布以 $P^{N}_{\theta_0}$-概率收缩于 $\theta_0$,即
推论 2.6(点估计存在性):存在 $\mathcal{A}$-取值估计量 $\hat{\theta}_N$,满足
此外,命题 2.3(反向 Poincaré 不等式):存在常数 $C > 0$,使得
3. 研究方法、关键技术和核心工具
核心障碍在于:标准反演理论(如 Nickl & Titi (2024))依赖于初值空间的线性结构与局部渐近正态性(LAN),而全局吸引子 $\mathcal{A}$ 是非线性的紧集,缺乏向量空间结构,导致 Fisher 信息等经典工具失效。本文的真正创新点在于将动力系统几何与贝叶斯非参数统计深度耦合,其技术支柱有三:
第一,建立并利用 $\mathcal{A}$ 上的反向 Poincaré 不等式(命题 2.3)。该不等式并非直接来自谱理论(因 $\mathcal{A}$ 非特征子空间),而是通过惯性流形 $\mathcal{M}$ 的精细构造获得:由定理 3.7,$\mathcal{A}$ 可唯一分解为 $u = u_n + \Phi(u_n)$,其中 $u_n = P_n[u] \in \mathrm{span}\{e_0,\dots,e_{n-1}\}$,$\Phi$ 是 $L^2$-Lipschitz 映射到正交补空间。由此导出
其中关键使用了 $\lambda_j \asymp j$ 的谱渐近与 $\|\cdot\|_{H^2} \asymp \sum_j (1+\lambda_j^2) |\langle \cdot, e_j \rangle|^2$ 的等价范数。
第二,将反向 Poincaré 不等式升级为解映射的 $L^2$-Lipschitz 稳定性(定理 2.4)。考虑差 $w(t) = u_\theta(t) - u_\vartheta(t)$,其满足
利用 $\mathcal{A}$ 的不变性($u_\theta(t), u_\vartheta(t) \in \mathcal{A} \subset \mathcal{M}$)及 $\mathcal{M}$ 的 $L^\infty$ 有界性($\|u\|_{L^\infty} \leq R$),结合 $f'$ 在 $[-R,R]$ 上的 Lipschitz 性,得 $|f(u_\theta)-f(u_\vartheta)| \leq c_f |w|$。再代入命题 2.3 的 $\|\nabla w(t)\|_{L^2} \leq c_P \|w(t)\|_{L^2}$,得到微分不等式
积分即得 Grönwall 型下界 $\|w(t)\|_{L^2}^2 \geq \|w(0)\|_{L^2}^2 e^{-(c_P^2 + c_f)t}$,从而证得 (1.2)。
第三,设计适配 $\mathcal{A}$ 几何的先验与后验分析框架。一方面,$\varepsilon$-网先验(2.1.1)直接利用 $\mathcal{A}$ 的有限覆盖数 $N(\mathcal{A}, \|\cdot\|_{L^2}, \varepsilon) \lesssim \varepsilon^{-D}$,通过加权 $\varepsilon_j$-覆盖构造支持在 $\mathcal{A}$ 上的先验,其小球概率 $\Pi(B'_{N}) \gtrsim \delta_N^{D+1}$ 满足贝叶斯收缩定理(Nickl (2023) 定理 1.3.2)所需的指数下界 $e^{-A N \delta_N^2}$。另一方面,惯性流形先验(2.1.2)将无限维问题降维:$\mathcal{M} = F(D(\Phi))$ 是 $n$-维椭球在 Lipschitz 映射下的像,故其覆盖数 $N(\mathcal{M}, \|\cdot\|_{L^2}, \varepsilon) \lesssim \varepsilon^{-n}$,且 $\theta_0 = F(\bar{\theta}_0)$ 的原像 $\bar{\theta}_0$ 位于 $D(\Phi)$ 内部,从而小球概率 $\Pi(B'_N) \gtrsim \delta_N^n$,同样满足收缩条件。两种先验均确保后验与估计量取值于 $\mathcal{A}$ 或 $\mathcal{M}$,保持了物理一致性。
4. 与之前工作的比较
与既有工作相比,本文在多个维度实现实质性推进:
首先,相对于 Nickl & Titi (2024) 对 2D 周期 Navier-Stokes 方程的研究,本文将 Lipschitz 稳定性从假设性的反向 Poincaré 不等式(需额外假设 $\|\nabla(\theta-\vartheta)\|_{L^2} \lesssim \|\theta-\vartheta\|_{L^2}$)转变为在反应-扩散系统中可被严格证明的定理(命题 2.3),其证明完全依赖于 $\mathcal{A}$ 的惯性流形结构与二维周期 Laplacian 的谱特性($\lambda_j \asymp j$),而非 Navier-Stokes 的复杂涡量结构。代价是目前结果局限于二维反应-扩散,作者明确指出“证明其他维数或 Navier-Stokes 的类似结果仍是巨大挑战”。
其次,相对于 Nickl (2024a) 和 Konen & Nickl (2026) 在 $t \approx 0$ 附近获得的代数速率 $N^{-\alpha}$,本文在 $t > 0$ 固定(即无零时刻数据)的苛刻条件下,仅凭 $\theta \in \mathcal{A}$ 的先验知识,就将速率从对数慢速 $1/\log N$ 提升至 $\sqrt{\log N / N}$。这一提升的关键在于:$\mathcal{A}$ 的有限维几何(覆盖维数 $D$ 或惯性维数 $n$)替代了时间局部性,使统计复杂度从无限维退化为有限维。然而,该速率仍慢于理想参数速率 $1/\sqrt{N}$,作者明确解释其根源在于 $\mathcal{A}$ 缺乏线性结构,导致无法定义 Fisher 信息与 LAN 展开,故不能直接套用 van der Vaart (1998) 的经典理论。
再次,相对于一般 Sobolev 类先验(如 $\theta \in H^\alpha$)下的贝叶斯分析(Stuart (2010), Cotter et al. (2009)),本文提出的 $\mathcal{A}$-支撑先验($\varepsilon$-网或惯性流形)具有根本性物理意义:它强制先验与后验严格位于系统实际可达的平衡态集合内,避免了传统高斯随机场先验将质量分配到非物理状态(如高能量、非吸引子轨道)的缺陷。这种先验的构造本身即是一项贡献,特别是惯性流形先验(11)将无限维推断问题转化为 $n$-维参数空间上的贝叶斯推断,其计算可行性远高于抽象 $\mathcal{A}$-支撑先验。
最后,在方法论上,本文拓展了 Monard et al. (2021) 和 Nickl (2023) 将贝叶斯非参数理论应用于 PDE 反问题的框架,首次将动力系统不变流形理论(Babin & Vishik, Temam)作为统计先验设计的严格数学基础。前人工作(如 Abraham & Nickl (2019) 对 Calderón 问题)处理的是线性化后的有限维参数,而本文处理的是非线性、非凸、无限维嵌入的紧集,其覆盖数估计(6)、小球概率下界(定理 2.1 证明)和 Lipschitz 参数化(13)均需针对 $\mathcal{A}$ 的特殊几何定制。从提供内容中无法确认本文方法是否可推广至带边界的有界域(如 $\Omega = (0,1)^2$),因全文基于周期边界与 Laplacian 的完整谱展开。
2. 具有水平耗散与小水平扰动的三维各向异性MHD方程组的整体存在性
Global Existence for 3D Anisotropic MHD system with Horizontal Dissipation and Small Horizontal Variations
作者: Qiliang Lin, Chenyin Qian, Daoyao Zhou
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.05764v1
摘要翻译: 本文建立了具有部分耗散的三维各向异性MHD方程组的整体适定性:速度场耗散项为 $Δ_\mathrm{h}u$,磁场耗散项为 $\partial_1^2b$,初值位于背景磁场 $(0,1,0)$ 附近。关键在于,仅需水平分量 $(u^\mathrm{h}_0,b^\mathrm{h}_0)$ 在 $H^2(\R^3)$ 中充分小,而垂直分量 $(u^3_0,b^3_0)$ 可任意大。我们的分析发展了若干新技巧,包括分量解耦型能量估计,以及利用背景磁场结构对危险非线性项进行迭代控制。该结果首次实现了各向异性MHD方程组在允许大垂直初值情形下的整体存在性,突破了以往工作所要求的全分量小初值条件。
Comment: 无
链接: http://arxiv.org/abs/2606.05764v1
1. 研究对象和背景
本文研究三维不可压各向异性磁流体动力学(MHD)系统,其速度场 $u = (u^\mathrm{h}, u^3) \in \mathbb{R}^3$ 与磁场 $b = (b^\mathrm{h}, b^3) \in \mathbb{R}^3$ 满足如下以背景磁场 $B^* = (0,1,0)$ 为基准的扰动形式:
其中 $\Delta_\mathrm{h} = \partial_1^2 + \partial_2^2$ 为水平拉普拉斯算子,$\nabla_\mathrm{h} = (\partial_1, \partial_2)$,$\mathrm{div}_\mathrm{h} = \partial_1 + \partial_2$;线性项 $\partial_2 b$ 与 $\partial_2 u$ 源于背景场 $B^*$ 的 Alfvén 波结构,提供弱稳定化效应。该模型属于混合部分耗散型 MHD 系统:速度仅具水平耗散 $\Delta_\mathrm{h} u$,而磁场仅在 $x_1$ 方向具二阶耗散 $\partial_1^2 b$,其余方向(如 $\partial_2^2$, $\partial_3^2$)无耗散。物理上,此类各向异性源于强背景磁场主导下的等离子体输运各向异性,广泛见于托卡马克装置、太阳日冕与行星发电机中。
核心数学挑战在于:标准能量方法要求所有分量均小以控制高阶非线性耦合,但物理现实(如对流胞、喷流不稳定性)允许垂直运动 $u^3_0$, $b^3_0$ 幅度极大。已有工作(如 [7, 13, 26, 27])虽处理同类各向异性耗散,却一律要求全初值 $(u_0,b_0) \in H^2(\mathbb{R}^3)$ 小,即 $\|(u_0,b_0)\|_{H^2} \leq \varepsilon$,从而将 $u^3_0$, $b^3_0$ 同样限制在小尺度内,严重脱离物理场景。本文旨在填补这一关键空白:是否可仅假设水平扰动 $(u^\mathrm{h}_0, b^\mathrm{h}_0)$ 小,而让垂直分量 $(u^3_0, b^3_0)$ 在 $H^2$ 中任意大?该问题本质是探究背景磁场 $B^* = (0,1,0)$ 与耗散模式 $(\Delta_\mathrm{h} u, \partial_1^2 b)$ 的协同稳定机制能否将控制权局部化至水平自由度。特别地,$\partial_1^2 b$ 耗散与 $B^*$ 的第二分量对齐,使得危险非线性项(如含 $\partial_3^2 b^\mathrm{h}$ 的项)可通过背景结构迭代控制,而非依赖 $b^3$ 的小性——这与 $B^* = (1,0,0)$ 搭配 $(\Delta_\mathrm{h} u, \partial_3^2 b)$ 的配置形成根本区别(后者因 $b^3$ 完全无耗散,强制 $b^3_0$ 小)。因此,本文试图建立首个允许垂直大初值的三维各向异性 MHD 整体适定性结果,突破“全小性”范式,为模拟强垂直扰动共存于受控水平波动的物理态提供严格数学基础。
2. 主要定理或主要结果
在初值 $(u_0,b_0) \in H^2(\mathbb{R}^3)$ 且满足 $\mathrm{div}\, u_0 = \mathrm{div}\, b_0 = 0$ 的前提下,存在绝对常数 $\varepsilon > 0$,使得若水平初值满足
则系统 (1.2) 存在唯一全局解
该结论由定理 1.1 给出,其核心在于条件 (1.3) 允许 $\|(u^3_0, b^3_0)\|_{H^2}$ 任意大:只要水平初值足够小(其容许大小随垂直初值增长呈双指数衰减),即可保证整体存在性。进一步,由命题 2.1 与 2.2 的精细估计可知,解的能量演化满足闭合不等式:对任意 $t \in [0,T^{**}]$,有
$$ \frac{1}{2} \frac{d}{dt} |(u^3,b^3)|{H^2}^2 + \frac{4}{25} |(\nabla\mathrm{h} u^3, \partial_1 b^3)|{H^2}^2 + \frac{3}{100} |(\partial_2 b^\mathrm{h}, \partial_2 b^3)|{H^1}^2 \leq C_1 |(\nabla_\mathrm{h} u^\mathrm{h}, \partial_1 b^\mathrm{h})|{H^2}^2 |(u^3,b^3)|{H^2}^2 + C_1 |(\nabla_\mathrm{h} u^\mathrm{h}, \partial_1 b^\mathrm{h})|{H^2}^2 \left(1 + |(u^\mathrm{h},b^\mathrm{h})|{H^2}^2 \right)
- \frac{1}{5} \frac{d}{dt} \int_{\mathbb{R}^3} \left( u^\mathrm{h} \partial_2 b^\mathrm{h} + u^3 \partial_2 b^3 + |\nabla u^\mathrm{h} \nabla \partial_2 b^\mathrm{h}| + |\nabla u^3 \nabla \partial_2 b^3| \right) dx. $$
该不等式通过时间积分与 Gronwall 型处理(利用 $\|(\nabla_\mathrm{h} u^\mathrm{h}, \partial_1 b^\mathrm{h})\|_{L^2_t H^2_x}^2 \leq \eta$ 的先验控制),结合 $Y(t) - Y(0)$ 的边界项吸收,最终导出 $\|(u^3,b^3)(t)\|_{H^2}^2$ 的一致有界性,从而完成整体性证明。
3. 研究方法、关键技术和核心工具
证明的核心障碍源于三重张力:(1)耗散各向异性导致标准 Sobolev 嵌入与能量估计失效,尤其 $b^\mathrm{h}$ 缺乏 $\partial_2^2$、$\partial_3^2$ 耗散,使 $\partial_3^2 b^\mathrm{h}$ 成为“危险因子”;(2)初值小性仅施加于水平分量,而垂直分量 $u^3$, $b^3$ 可任意大,传统统一能量法无法闭合;(3)压力项 $-\nabla p$ 在分量分解后失去 $\int u \cdot \nabla p$ 的自动抵消,其水平与垂直贡献互为对偶,需兼容估计。
作者的真正创新点在于构建一套“分量解耦—迭代控制—结构利用”的三级分析框架:
第一,分量解耦能量法(Stage 1)。定义独立能量泛函 $E_\mathrm{v}(t) = \frac{1}{2} \|(u^3,b^3)\|_{H^2}^2$ 与 $E_\mathrm{h}(t) = \frac{1}{2} \|(u^\mathrm{h},b^\mathrm{h})\|_{H^2}^2$,并分别推导其演化方程。关键在于利用不可压条件 $\partial_3 u^3 = -\mathrm{div}_\mathrm{h} u^\mathrm{h}$、$\partial_3 b^3 = -\mathrm{div}_\mathrm{h} b^\mathrm{h}$,将垂直高阶导数(如 $\partial_3^2 u^3$)转化为水平导数作用于小量:$\partial_3^2 u^3 = -\partial_3 \mathrm{div}_\mathrm{h} u^\mathrm{h}$。此“导数转移机制”系统应用于所有垂直高阶项,例如将临界积分
重写为仅含至多一个高阶垂直导数的组合,最终获得可控界
从而将 $E_\mathrm{v}$ 的耗散项 $\|(\nabla_\mathrm{h} u^3, \partial_1 b^3)\|_{H^2}^2$ 与水平小量耦合。
第二,压力项对偶估计(Stage 2)。由 $-\Delta p = \mathrm{div}\, \mathrm{div}\, (u \otimes u - b \otimes b)$ 解出 $p$,其难点在于形如 $\langle (-\Delta)^{-1} \partial_3^2 (u^3)^2, \mathrm{div}_\mathrm{h} u^\mathrm{h} \rangle_{L^2}$ 的项含过多 $u^3$。作者发现该压力项在水平与垂直能量平衡中实为同一量的不同表现:$-(\nabla_\mathrm{h} p | u^\mathrm{h})_{L^2} = (p | \mathrm{div}_\mathrm{h} u^\mathrm{h})_{L^2} = -(p | \partial_3 u^3)_{L^2} = (\partial_3 p | u^3)_{L^2}$。为此,发展两种组件特异性估计:对水平能量,利用水平 Gagliardo–Nirenberg 不等式 $\|f\|_{L^4_h} \lesssim \|f\|_{L^2_h}^{1/2} \|\nabla_\mathrm{h} f\|_{L^2_h}^{1/2}$ 与 Riesz 算子在混合范数 $L^p_h(L^q_v)$ 中的有界性,将 $u^3$ 的高阶范数转化为 $u^\mathrm{h}$ 的低阶范数;对垂直能量,则直接控制 $\|\nabla_\mathrm{h} u^3\|_{H^2}$ 与 $\|\nabla_\mathrm{h} u^\mathrm{h}\|_{H^2}$ 的乘积。
第三,迭代非线性分析(Stage 3)。针对 $b^\mathrm{h}$ 的 $H^2$ 估计中最危险项 $\langle \partial_i^2 (u \cdot \nabla b^\mathrm{h}), \partial_i^2 b^\mathrm{h} \rangle_{L^2}$,作者引入突破性迭代程序:利用磁场方程 $\partial_t b + u \cdot \nabla b - \partial_1^2 b = b \cdot \nabla u + \partial_2 u$,将 $\partial_2 u$ 替换为 $\partial_t b + u \cdot \nabla b - \partial_1^2 b - b \cdot \nabla u$,从而将危险因子 $\partial_3^2 b^\mathrm{h}$ 上的 $\partial_2$ 导数转移至时间导数或 $\partial_1^2$ 耗散项。反复应用此替换,生成含 $\partial_t \partial_3^2 b^\mathrm{h}$ 的项,再用方程二次替换,最终揭示隐藏对称性:经分部积分,高阶项 $-\int u \cdot \nabla b^i \, \partial_3^2 u^j \, \partial_3^2 b^k dx$ 与 $\int u \cdot \nabla b^i \, \partial_3^2 u^k \, \partial_3^2 b^j dx$ 相消,剩余项均为低阶可控形式。此过程由引理 1.2 的各向异性插值不等式精确支撑,例如
确保 $\partial_1$(唯一有耗散的方向)被最优分配至最危险因子。
4. 与之前工作的比较
本文与既有工作的比较可从四个维度展开:
(1)与全耗散 MHD 的对比:全耗散情形($\mu_i, \nu_i > 0$ 对所有 $i$)下,标准能量法可直接应用,小初值整体适定性属常规推广 [6,21];近期 [2] 虽处理大初值,但仍依赖全耗散以保证高频阻尼。本文仅具混合部分耗散,且无任何 $\partial_2^2$、$\partial_3^2$ 项作用于 $b$,故技术难度远超全耗散框架,其成功依赖于背景场结构而非额外耗散。
(2)与单一部耗散 MHD 的对比:Lin 等 [11] 及后续工作研究仅有速度耗散或仅有磁场耗散的情形,但均要求全初值小。本文耗散更弱(仅 $\Delta_\mathrm{h} u$ 与 $\partial_1^2 b$),却实现更强结论(仅水平小),表明背景场 $B^* = (0,1,0)$ 的定向稳定效应可部分替代缺失耗散。
(3)与混合耗散但全小性要求的工作对比:文献 [7,13,26,27] 均处理混合耗散,但无一例外要求 $\|(u_0,b_0)\|_{H^2} \leq \varepsilon$。例如 [26] 考察 $B^* = (1,0,0)$ 与 $(\Delta_\mathrm{h} u, \partial_3^2 b)$,其危险项 $\int \partial_1 u_1 \, \partial_2^2 b^3 \, \partial_2^2 b^3 dx$ 含 $\partial_2^2 b^3$,而 $b^3$ 完全无耗散,故强制 $b^3_0$ 小。本文配置 $(\Delta_\mathrm{h} u, \partial_1^2 b)$ 与 $B^* = (0,1,0)$ 则使同类危险项变为 $\int \partial_2 u_2 \, \partial_3^2 b^\mathrm{h} \, \partial_3^2 b^\mathrm{h} dx$,其中 $\partial_3^2 b^\mathrm{h}$ 可通过迭代程序(利用 $\partial_1^2 b$ 耗散与 $\partial_2 u$ 的背景结构)控制,无需 $b^3$ 小。此结构性差异是本文突破的根源。
(4)与各向异性 Navier-Stokes 的对比:当 $b \equiv 0$,系统退化为水平耗散 NS 方程,Chemin 等 [5] 已证水平小性结果。但 MHD 的强耦合($u$ 与 $b$ 相互驱动、Alfvén 波传播)使 NS 技术完全失效:NS 中压力与非线性项可独立处理,而 MHD 中 $b \cdot \nabla u$ 与 $u \cdot \nabla b$ 形成闭环,且 $\partial_2 u$、$\partial_2 b$ 的线性项引入新耦合通道。本文的分量解耦法、压力对偶估计及迭代非线性分析均为 MHD 特有设计,无法从 NS 结果直接类比。代价在于:条件 (1.3) 中水平小性需随垂直初值呈双指数衰减,较全小性假设更苛刻;且解的空间正则性受限于耗散模式,如 $\partial_2 b$ 仅在 $L^2_t H^1_x$ 中有界,无法获得更高阶垂直正则性。从提供内容中无法确认该结果是否可推广至 $H^s$($s>2$)或 Besov 空间。
了解领域论文(摘要翻译)
1. 椭圆系统在角域中算子铅笔谱界结果的补遗
Addendum to Spectral bounds for the operator pencil of an elliptic system in an angle
作者: Michael Tsopanopoulos
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.06264v1
链接: http://arxiv.org/abs/2606.06264v1
摘要翻译
中文标题: 椭圆系统在角域中算子铅笔谱界结果的补遗
摘要翻译: 本补遗将原始论文中建立的谱下界结果推广至 Neumann 问题。在椭圆组元的形式正定性假设下,我们证明:与 Dirichlet 情形相同,这些谱下界对 Neumann 边界条件亦成立。
2. Vlasov-Poisson 系统的同胚型修正波算子
Homeomorphic modified wave operators for the Vlasov-Poisson system
作者: Léo Bigorgne
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.06488v1
链接: http://arxiv.org/abs/2606.06488v1
摘要翻译
中文标题: Vlasov-Poisson 系统的同胚型修正波算子
摘要翻译: 我们在一个函数框架下证明了 Vlasov-Poisson 系统小初值解的修正散射性,其中初值、散射态及渐近收敛性均在同一拓扑下度量。此外,我们证明相应的波算子在初值空间与散射态空间之间定义了同胚映射,并且在更弱范数下具有局部 Lipschitz 连续性。作为推论,在排斥情形下,大球对称解是渐近稳定的。该证明特别依赖于引入一组适配于渐近非线性流的动态坐标系。
3. 能量空间中非聚焦 $k$-色散广义 Benjamin-Ono 方程的大初值散射
Large data scattering for the defocusing $k$-dispersion generalized Benjamin-Ono equation in the energy space
作者: Luccas Campos, Felipe Linares, Thyago S. R. Santos
类别: math.AP, math.CA
发布日期: 2026-06-04
arXiv ID: 2606.06463v1
链接: http://arxiv.org/abs/2606.06463v1
摘要翻译
中文标题: 能量空间中非聚焦 $k$-色散广义 Benjamin-Ono 方程的大初值散射
摘要翻译: 我们研究非聚焦 $k$-色散广义 Benjamin-Ono 方程。对每个偶整数 $k\geq 4$,我们证明:初始数据属于能量空间 $H^{\frac{\alpha}{2}}$ 的解整体存在且散射。证明结合了 Kenig 与 Merle 的集中紧致性-刚性方法,以及基于 Caffarelli-Silvestre 扩展和 Tao 单调性公式的技术——该单调性公式已适配于分数阶色散情形。
4. 从外部测量中恢复稳定核
Recovering stable kernels from exterior measurements
作者: Yi-Hsuan Lin
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.06427v1
链接: http://arxiv.org/abs/2606.06427v1
摘要翻译
中文标题: 从外部测量中恢复稳定核
摘要翻译: 我们研究一类平移不变对称稳定算子的反问题,其形式为 \begin{equation*} L_a u(x)=\mathrm{P.V.}\int_{\mathbb R^n}(u(x)-u(y))\frac{a((x-y)/|x-y|)}{|x-y|^{n+2s}},dy, \quad 0<s<1, \end{equation*} 其中未知量为球面 $\mathbb S^{n-1}$ 上的偶角密度 $a$。对有界开集 $Ω\subset\mathbb R^n$,记 $Ω_e=\mathbb R^n\setminus\overlineΩ$,我们考虑受限的外部 Dirichlet-to-Neumann 映射 $Λ_a^{W_1,W_2}$,其中外部数据支集于 $W_1\SubsetΩ_e$,而非局部 Neumann 数据在 $W_2\SubsetΩ_e$ 上被观测。我们证明关于主阶角密度的三项恢复结果:在重叠情形 $W_1\cap W_2\ne\emptyset$ 下,外部对角奇性唯一确定每个光滑椭圆角密度;在分离情形 $\overline W_1\cap\overline W_2=\emptyset$(此时该奇性消失)下,我们通过稳定符号的精确因式分解,证明了有限调和角类中的唯一性;此外,当源集与观测集均位于无界外部连通分支中时,我们利用非对角 Dirichlet-to-Neumann 核的解析延拓及远场渐近论证,证明了实解析角密度的分离数据唯一性。
5. 超椭圆阻尼下波动方程的最优衰减
Optimal decay for waves damped by superellipses
作者: B. Achammer, Perry Kleinhenz
类别: math.AP, math.SP
发布日期: 2026-06-04
arXiv ID: 2606.06166v1
链接: http://arxiv.org/abs/2606.06166v1
摘要翻译
中文标题: 超椭圆阻尼下波动方程的最优衰减
摘要翻译: 环面上阻尼波动方程解的能量衰减率已知依赖于阻尼集的几何结构以及阻尼函数的增长性质。本文针对一类阻尼函数建立能量衰减率的下界估计:该阻尼函数在超椭圆上为正,且以到超椭圆边界的距离的多项式形式增长,即形如 $\mathrm{dist}(x, \partial \mathcal{E}_p)^\gamma$,其中 $\mathcal{E}_p = \{(x,y) \in \mathbb{T}^2 : |x|^p + |y|^p \leq 1\}$ 为指数为 $p > 2$ 的超椭圆,$\gamma > 0$。我们所得的能量衰减率显式依赖于超椭圆定义中的指数 $p$ 和多项式幂次 $\gamma$;进一步,我们证明这些衰减率在某些参数范围内是最优的。证明通过简化标准的法型(normal form)论证,将源于 $y$-不变阻尼的拟模态(quasimodes)适配到当前设定中。
6. 一个离散时间格式对退化抛物-双曲系统的收敛性
Convergence of a discrete-in-time Approximation to a Degenerate Parabolic-Hyperbolic System
作者: Julia Hauser, Hideki Murakawa, Markus Schmidtchen
类别: math.AP, math.NA
发布日期: 2026-06-04
arXiv ID: 2606.05879v1
链接: http://arxiv.org/abs/2606.05879v1
摘要翻译
中文标题: 一个离散时间格式对退化抛物-双曲系统的收敛性
摘要翻译: 本文考虑一类退化反应-交叉扩散系统的隐式半离散逼近。由于抛物部分具有对称性,该系统已知保持密度的分离性——初始时刻互不重叠、属于不同物种的密度在所有时刻均保持分离,从而在不同物种之间形成内部层。我们证明:当时间步长趋于零时,时间离散逼近存在并收敛于一个弱解。
7. 临界 $p$-拉普拉斯情形下条件 Harnack 不等式的 Pohozaev 型细颈证明
A Pohozaev-type neck proof of a conditional Harnack inequality in the critical $p$-Laplacian setting
作者: Guolin Qin, Yi Ru-Ya Zhang
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.05990v1
链接: http://arxiv.org/abs/2606.05990v1
摘要翻译
中文标题: 临界 $p$-拉普拉斯情形下条件 Harnack 不等式的 Pohozaev 型细颈证明
摘要翻译: 我们在全局临界 Sobolev 增长假设及单调性条件 $s^{-(p^*-1)}g(s)$ 非增的前提下,证明了临界 $p$-拉普拉斯方程
的正弱解的条件 Schoen 型 Harnack 不等式。该结果依赖于两个前提:一是将有界正整体爆破极限分类为 Aubin–Talenti $p$-气泡;二是对归一化细颈给出初步的奇性速率上界控制。在上述两个假设下,定义于 $B_{3R}$ 中的解满足
核心论证是一种 Pohozaev-细颈方法,它将初步奇性衰减率 $|x|^{-(n-p)/p}$ 提升至尖锐的 $p$-调和基本解速率 $|x|^{-(n-p)/(p-1)}$。该论证替代了共形不变半线性情形 $p=2$ 下可用的 Kelvin 变换与移动球面法,而后者在一般 $p$-拉普拉斯情形下不可用。
8. 聚焦质量临界NLS在 $\Bbb R^d\times\Bbb T$ 上的Dancer型解的尖锐质量阈值
Sharp mass-threshold for Dancer-type solutions of the focusing mass-critical NLS on $\Bbb R^d\times\Bbb T$
作者: Yongming Luo
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.05707v1
链接: http://arxiv.org/abs/2606.05707v1
摘要翻译
中文标题: 聚焦质量临界NLS在 $\Bbb R^d\times\Bbb T$ 上的Dancer型解的尖锐质量阈值
摘要翻译: 欧氏空间 $\R^d$ 上的质量临界NLS展现出强质量刚性:所有正基态均由单一剖面生成,且具有相同的基态质量 $\widehat{M}(Q)$。Dancer在其著名论文 \cite{DancerSolution} 中借助分歧方法构造了相应方程在 $\R^d\times\T$ 上的解,这些解在非紧方向上衰减,并在一个方向上非平凡地周期振荡。然而,此类分歧方法并未给出解的能量刻画,尤其未能阐明其与欧氏基态之间的关系。本文通过引入一种新的关于预设质量能量水平的严格单调性机制,并结合作者近期工作中发展的半维里尔-消没几何框架,证明:对任意质量 $c\in(0,2π\widehat{M}(Q))$,半维里尔-消没变分问题 $m_c$ 存在一个归一化的Dancer型极小元,该极小元亦满足 $\R^d\times\T$ 上的聚焦质量临界NLS。该结果亦为我们在前期工作 \cite{Luo_LegendreFenchel} 中利用Legendre-Fenchel对偶所导出的存在性结果提供了尖锐补遗。
9. 混合空间阶的界面问题
Interface problems of mixed spatial order
作者: Dionyssios Mantzavinos, Ravindra Pethiyagoda, Dave Smith
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.05593v1
链接: http://arxiv.org/abs/2606.05593v1
摘要翻译
中文标题: 混合空间阶的界面问题
摘要翻译: 我们在直线上求解各类常系数线性发展型偏微分方程之间的界面问题。我们的典型例子包括热方程、线性薛定谔方程、Airy 方程、线性化 Korteweg–de Vries 方程以及双谐薛定谔方程。在每个问题中,所列方程之一被定义在一个空间半直线 $(-\infty,0)$ 上,另一个则定义在另一半直线 $(0,\infty)$ 上,并配有适当的界面条件。这些问题通过 Fokas 统一变换方法的新颖推广得以解决,且显式解公式利用 Filon 数值积分法进行计算。
10. 具有Hartree型非局部非线性项叠加的临界方程的解
Solutions to critical equations with a superposition of nonlocal Hartree-type nonlinearities
作者: Artur Jorge Marinho, Kanishka Perera
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.05592v1
链接: http://arxiv.org/abs/2606.05592v1
摘要翻译
中文标题: 具有Hartree型非局部非线性项叠加的临界方程的解
摘要翻译: 我们研究一类定义在 $\mathbb{R}^N$ 上的非线性非局部椭圆方程,其非线性项为 Hartree 型非线性项的叠加。受 Schrödinger-Poisson-Slater 系统启发,这类方程自然地推广了仅含单一非局部相互作用项的问题。更精确地,我们考虑由一族以正 Borel 测度加权的 Riesz 位势所驱动的方程,从而产生一个叠加的非局部算子。为变分地处理该问题,我们引入合适的作用空间,即叠加 Coulomb 空间及其相应的叠加 Coulomb-Sobolev 空间,并研究其主要性质。结合变分方法与近期发展的基于缩放的临界点理论,我们证明了径向解的存在性与多重性结果。我们还研究了一类 Brezis-Nirenberg 型问题,并在相关非线性特征值问题的特征值附近获得了多重性结果。我们的结果推广了此前关于单 Hartree 型方程的工作,并为处理 Hartree 型非局部相互作用的叠加提供了一个统一框架。
11. 通过切比雪夫型积分不等式实现的鲁棒插值不等式
Robust interpolation inequalities via Chebyshev-type integral inequalities
作者: Guy Foghem
类别: math.AP
发布日期: 2026-06-03
arXiv ID: 2606.05477v1
链接: http://arxiv.org/abs/2606.05477v1
摘要翻译
中文标题: 通过切比雪夫型积分不等式实现的鲁棒插值不等式
摘要翻译: 我们在Gagliardo半范数尺度下建立了鲁棒的对数凸插值不等式。为此,我们推导了适用于一般非同步函数的若干切比雪夫型积分不等式。建立这些鲁棒插值不等式的主要动机源于对区域分数阶 $p$-拉普拉斯算子相关边界Dirichlet问题弱解的渐近非局部到局部稳定性的研究。更准确地说,若 $u_s \in W^{s,p}(Ω)$ 弱满足 $(-Δ)_{p, Ω}^s u_s = f_s $ 在 $Ω$ 中,且 $ γ^s_0(u_s) = g_s$ 在 $\partialΩ$ 上,其中 $\frac{1}{p} < s \leq 1$,且 $Ω\subset \mathbb{R}^d$ 是有界Lipschitz区域,则在数据 $f_s$ 与 $g_s$ 当 $s \to 1^-$ 时适当收敛的条件下,我们证明了 $\| u_s - u_1 \|_{W^{η,p}(Ω)} \xrightarrow{s \to 1^-} 0 $ 对所有 $0 \leq η< 1$ 成立。
12. 薄区域、约化与慢流形
Thin Domains, Reduction, and Slow Manifolds
作者: Christian Kuehn, Jan-Eric Sulzbach
类别: math.AP
发布日期: 2026-06-03
arXiv ID: 2606.05398v1
链接: http://arxiv.org/abs/2606.05398v1
摘要翻译
中文标题: 薄区域、约化与慢流形
摘要翻译: 我们为定义在薄区域上的偏微分方程提出一个统一的维数约化框架。我们的方法融合了三个互补要素:细致的边界条件分析、适用于一般薄几何结构的基于平均化的分裂方法,以及对所得快慢系统所采用的慢流形观点。薄面(thin faces)上的齐次Neumann条件作为最相关且物理上自然的设定而浮现,因其保持横向零模(transverse zero mode),从而导出真正低维的约化方程。对于一般的薄区域,我们推导出平均化后的快慢系统,并分离出由几何引起的修正项。接着,当慢变量存在合适的谱分解时,我们构造一种基于分裂的Lyapunov-Perron 方法以得到精确慢流形;同时,我们亦直接从不变性方程出发,通过渐近展开构造近似慢流形及经修正的约化动力学。此外,我们提出Schnakenberg 反应-扩散系统作为一个典型测试问题,用于比较完整薄区域动力学、平均化模型以及流形修正后的约化动力学。最后,我们还将该框架推广至薄管状区域,并在弯曲几何中导出相应的重标度快慢系统。
不相关论文(基本信息)
1. 一个完全非线性椭圆方程的反源问题
An inverse source problem for a fully nonlinear elliptic equation
作者: Ching-Lung Lin, Yi-Hsuan Lin, Jenn-Nan Wang
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.06431v1
链接: http://arxiv.org/abs/2606.06431v1
摘要: We study an inverse source problem for fully nonlinear elliptic equations of the form
The question is whether the source term can be recovered from the Dirichlet-to-Neumann map. In two dimensions, the first linearization does not immediately give uniqueness: it leaves a natural conformal ambiguity in the linearized coefficients. For homogeneous nonlinearities $F$ with injective differential $DF$, we show that this ambiguity has a precise meaning at the level of the equation itself, namely that the source is determined up to an explicit scalar factor. The main point of the paper is to show how this remaining factor can be removed. We use the second linearization to extract information which is invisible at first order, and combine it with an algebraic nondegeneracy condition on the nonlinearity. Under this condition, the residual ambiguity is forced to be trivial, and the Dirichlet-to-Neumann map uniquely determines the source. The result applies, in particular, to homogeneous admissible Hessian equations of Monge–Ampère type and related examples.
2. 含指数项与梯度反应的椭圆方程解的局部与整体性质
Local and global properties of solutions of an elliptic equation involving exponential and gradient reaction
作者: Marie-Françoise Bidaut-Véron, Marta Garcia-Huidobro, Laurent Véron
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.06206v1
链接: http://arxiv.org/abs/2606.06206v1
摘要: We study some local and global properties of solutions of $-Δu- m\abs{\nabla u}^q-e^{u}=0$ in a punctured domain $Ω\setminus\{0\}$, or in an exterior domain of $R^N$, $N\geq 2$, where $m$ is a positive parameter and $q>1$. We study particularly the local behaviour of solutions with an isolated singularity or the asymptotic behaviour for solutions defined in an exterior domain, and also the existence of solutions with the behaviours previously described. These behaviours change drastically according $q$ is smaller or larger than $2$. Many results are obtained by introducing various dynamical systems associated to the equation.
3. Canham-Helfrich 能量下近圆柱形管的微小形变及其在生物膜中的应用
Small deformations of a near cylindrical tube for the Canham-Helfrich Energy with applications to biological membranes
作者: Charles M. Elliott, Carsten Gräser, Philip J. Herbert
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.06172v1
链接: http://arxiv.org/abs/2606.06172v1
摘要: In this article we develop a quadratic energy which approximates the Canham-Helfrich energy for a tube-like surface with clamped boundary and area constraint. The energy is suited to the study of small deformations of biological membranes where the deformations are induced by point forces or point constraints due to the cytoskeleton or a phase dependent spontaneous curvature. Since the deformations we consider are small, we may assume that the surface of interest is a graph over a fixed, undeformed surface. A Lagrangian and the associated Euler-Lagrange equations for the graph are derived. Well-posedness of the Euler-Lagrange equations in suitable spaces is shown. Finally, we provide some illustrative numerical examples.
4. 为变分相场断裂模型赋予定制强度准则
Endowing variational phase-field fracture models with custom strength criteria
作者: Roberto Alessi, Matteo Brunetti, Roshan Udaram Patil, Jacinto Ulloa
类别: cond-mat.mtrl-sci, math.AP
发布日期: 2026-06-04
arXiv ID: 2606.06169v1
链接: http://arxiv.org/abs/2606.06169v1
摘要: By now, several approaches have been proposed to endow phase-field fracture models with the ability to describe crack nucleation under multiaxial stress states. These include techniques for splitting the free energy, direct modifications of the phase-field driving or resisting forces that sacrifice the variational structure of the problem, and the introduction of additional internal variables, such as plastic strains or other nonlinear strains. In this paper, we propose a fundamentally different strategy for incorporating arbitrary elastic domains into phase-field fracture models, formulated within the variational framework of generalized standard materials. The proposed approach relies on letting the dissipation potential depend on the current state of the material. In this way, the variational structure of the problem is preserved, while elastic degradation and the strength criterion remain two distinct and independently controllable aspects of the material response. Simple yet representative models are presented and thoroughly discussed to demonstrate the effectiveness of the proposed methodology. The resulting evolution of the elastic domain is investigated in both strain and stress spaces. Moreover, numerical simulations demonstrate a range of crack nucleation processes under multiaxial loading conditions for various analytical strength surfaces. This work paves the way for future developments and applications in several directions.
5. 宽两层神经网络中不确定性的量化:关于极限涨落过程的律
Quantifying Uncertainty In Wide Two-Layer Neural Networks: On The Law Of The Limiting Fluctuation Process
作者: Arnaud Descours, Arnaud Guillin, Geoffrey Lacour, Manon Michel, Boris Nectoux, Paul Stos
类别: cs.NE, math.AP, math.PR
发布日期: 2026-06-04
arXiv ID: 2606.05982v1
链接: http://arxiv.org/abs/2606.05982v1
摘要: Uncertainty quantification in neural networks prediction is a main issue for usual applications. Our approach seeks at reducing computation costs by directly evaluating uncertainty using PDE’s information on the asymptotic variance, rather than the deep ensemble method which may be seen as a Monte Carlo estimation of the prediction, requiring the training of multiple networks. We thus study the law of the limiting process describing the random fluctuations around the mean-field limit of wide two-layer neural networks trained by stochastic gradient descent in a weak-noise regime. Building on a recent trajectorial central limit theorem, in which this limit is characterized as the weak solution of a linear stochastic evolution equation, we identify its law explicitly. More precisely, we show that it is a centered Gaussian process in the dual of a weighted Sobolev space, and we derive a closed covariance representation for the finite-dimensional distributions obtained by testing it against smooth functions. This covariance is expressed through the solution of a backward transport equation with a nonlocal source term, whose coefficients are driven by the mean-field trajectory. As a consequence, by testing against the activation function at a fixed input, we obtain an expression for the limiting variance of the corresponding network-output fluctuations. We illustrate this result numerically on a one-dimensional regression example.
6. 通过Keller-Segel系统中信号生成的局部各向异性防止 $L^p$ 爆破:适用于强不同扩散速率情形
Preventing $L^p$ blow-up by local anisotropy of signal production in the Keller-Segel system with strongly differing diffusion rates
作者: Youshan Tao, Michael Winkler
类别: math.AP
发布日期: 2026-06-04
arXiv ID: 2606.06009v1
链接: http://arxiv.org/abs/2606.06009v1
摘要: In a smoothly bounded domain $Ω\subset R^n$, $n\le 5$, the manuscript considers the variant of the Keller-Segel system given by
which involves an additional contribution $\nabla \cdot (u\nabla v)$ to the chemoattractant evolution, in line with refined modeling literature reflecting an anisotropic correction to the isotropic signal production term $+u$ in the classical Keller-Segel model. It is shown that for arbitrary $D>0$ and $d>0$ and any nonnegative intial data from $W^{1,\infty}(Ω)\times W^{1, \infty}(Ω)$, an associated Neumann problem admits a global weak solution $(u,v)$ which, inter alia, satisfies
with some $α>0$ and some null set $N\subset (0,\infty)$.
7. 稳定性边界上的梯度下降:双层网络的自由能模型与动力学描述
Gradient descent at the Edge of Stability: free energy model and kinetic description of the two-layer network
作者: Antonin Chodron de Courcel
类别: math.OC, cs.AI, cs.LG, math-ph, math.AP
发布日期: 2026-06-03
arXiv ID: 2606.05326v1
链接: http://arxiv.org/abs/2606.05326v1
摘要: We study the dynamics of gradient descent in the Edge of Stability regime, where the learning rate is large enough to induce persistent oscillations in the loss and the sharpness. We propose a continuous-time effective model that tracks the evolution of the average trajectory coupled with the time-averaged covariance of its fast oscillations. Our analysis reveals that the natural quantity to monitor in such unstable regimes is an effective free energy, which combines the original risk functional with a curvature-related “entropic” term. Our model allows us to track the envelope of the oscillations even in situations where its dynamics evolve on similar timescales as the averaged weights. Otherwise stated, we can track the spikes that occur during the training of some neural network architectures. For wide two-layer neural networks optimized under stable non-vanishing oscillations, we derive a mean-field limit that results in a novel kinetic equation describing the joint distribution of weights and their fluctuations. We show that this equation can be interpreted as a Wasserstein-2 gradient flow of a macroscopic free energy. Finally, we provide numerical evidence on matrix factorization and deep learning tasks (CIFAR-10) to demonstrate the model’s accuracy in capturing the envelope of the oscillations and the predictive power of the effective free energy.