生成时间: 2026年02月23日 04:41:01

重点关注论文数量: 2

了解领域论文数量: 3

不相关论文数量: 13

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

1. Stability of the Shape for Circular Vortex Filaments under Non-Symmetric Perturbations

作者: Masashi Aiki, Mitsuo Higaki

类别: math.AP

发布日期: 2026-02-20

arXiv ID: 2602.18155v1

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

详细分析

研究对象和背景

本文研究的是三维空间中的局部感应方程（Localized Induction Equation，简称 LIE）所描述的封闭涡丝（vortex filament）在非对称扰动下的形状稳定性。LIE 是一个经典的几何偏微分方程，用于描述细长涡丝在自身感应速度作用下的运动。其数学形式为

\mathbf{x}_t = \mathbf{x}_s \times \mathbf{x}_{ss}, \quad (s,t) \in \mathbb{T} \times (0,\infty)

其中 $\mathbf{x}(s,t)$ 是涡丝在时间 $$t$$ 的参数化曲线， $$s$$ 是弧长参数， $\mathbb{T}$ 表示长度为 $$L$$ 的平环面（即 $\mathbb{R}/L\mathbb{Z}$ ），下标 $$s$$ 和 $$t$$ 分别表示对弧长和时间的偏导数。方程右边的速度场是涡丝切线方向 $\mathbf{x}_s$ 与曲率向量 $\mathbf{x}_{ss}$ 的叉积，这源于 Biot-Savart 定律在涡丝细长极限下的局部近似。LIE 具有刚体运动对称性：对任意三维平移向量 $\boldsymbol{\tau} \in \mathbb{R}^3$ 和绕 $$z$$ 轴的旋转 $Q_\alpha$ ，若 $\mathbf{x}(s,t)$ 是解，则 $Q_\alpha\mathbf{x}(s,t) + \boldsymbol{\tau}$ 也是解。

论文固定涡丝截面的半径 $$R>0$$ ，从而环面长度确定为 $L = 2\pi R$ 。研究的具体对象是圆形涡丝，其静止参考构型为

\mathbf{x}^R_0(s) = R\cos(s/R)\,\mathbf{e}_1 + R\sin(s/R)\,\mathbf{e}_2,

其中 $\{\mathbf{e}_1,\mathbf{e}_2,\mathbf{e}_3\}$ 是标准欧几里得基。这是一个精确的 LIE 解，描述了一个在 $$xy$$ 平面内半径为 $$R$$ 的圆，其整体运动是沿 $$z$$ 轴方向的匀速平移：

\mathbf{x}^R(s,t) = \mathbf{x}^R_0(s) + t\,\mathbf{e}_3.

这种运动对应于涡丝保持圆形不变，同时以恒定速度沿涡丝轴（ $$z$$ 轴）上升，是 LIE 中一个重要的刚性运动解。

论文关注的是形状稳定性：若初始涡丝 $\mathbf{x}_0(s) = \mathbf{x}(s,0)$ 接近圆环 $\mathbf{x}^R_0$ （在某种度量下），随后解 $\mathbf{x}(t)$ 是否在长时间内保持接近“圆形的形状”？这里的关键是“形状”应排除刚体运动（平移和轴向旋转）的影响，因为平移不改变曲线几何，而绕 $$z$$ 轴旋转圆环仍得到一个圆（只是参数相位变化）。因此，稳定性度量需对平移和 $$z$$ 轴旋转取轨范（orbital）距离。论文定义由所有旋转和平移生成的流形

\Sigma = \{Q_\alpha\mathbf{x}^R_0(\cdot) + \boldsymbol{\tau} \mid \alpha \in \mathbb{R}, \boldsymbol{\tau} \in \mathbb{R}^3\} \subset H^2(\mathbb{T})^3,

并考虑解到该流形的 $$H^2$$ 距离

\operatorname{dist}(\mathbf{x}(t),\Sigma) = \inf_{\alpha\in\mathbb{R},\boldsymbol{\tau}\in\mathbb{R}^3} \|\mathbf{x}(t) - (Q_\alpha\mathbf{x}^R_0 + \boldsymbol{\tau})\|_{H^2(\mathbb{T})}.

这是一个非线性轨道稳定性的严格定义。

函数空间方面，论文工作在 Sobolev 空间 $H^m(\mathbb{T})^3$ 上，其中 $H^m(\mathbb{T})$ 是通常的周期 Sobolev 空间，范数记为 $\|\cdot\|_m$ ， $\|\cdot\|=\|\cdot\|_0$ 是 $$L^2$$ 范数。初值需满足 $\mathbf{x}_0 \in H^4(\mathbb{T})^3$ 以保证解的唯一性（引用 Tani-Nishiyama 的工作）。此外，弧长保持条件 $|\mathbf{x}_s(s,t)|=1$ 对 $$L^2$$ 函数并非自动满足，需要初值满足 $|\mathbf{x}_{0s}|=1$ ，并由 LIE 的几何结构保持。

物理上，涡丝是理想流体中涡量集中的管状结构，LIE 作为中尺度模型广泛用于涡旋动力学。圆环解是唯一已知的具有显式刚体运动的封闭涡丝。然而，它存在线性不稳定性：平移自由度（零频模）会导致解在 $$L^2$$ 距离下线性偏离（Lyapunov 不稳定）。因此，论文的目标是证明：尽管存在这种平移漂移，涡丝的几何形状（即曲线的弯曲和扭率结构）在模去平移和轴向旋转后是稳定的。这需要精细分析各频段（零模、±1模、高频模）的动力学，并利用守恒律控制低频模的长期行为。

主要定理或主要结果

论文的核心结果是以下定理，它建立了圆形涡丝在非对称扰动下的非线性轨道稳定性。

Theorem 1.1 圆形涡丝 $\mathbf{x}^R_0$ 在以下意义上是轨道稳定的：存在正常数 $$C$$ 和 $\delta_0$ ，使得对任意 $0 < \delta \le \delta_0$ 及任意初值 $\mathbf{x}_0 \in H^4(\mathbb{T})^3$ 满足

|P[\mathbf{x}_0] - \pi R^2\mathbf{e}_3| \le \delta^2, \quad \|P_{|n|\ge2}\mathbf{x}_0\|_{H^2(\mathbb{T})} \le \delta,

且 $|\mathbf{x}_{0s}|=1$ ，则 LIE 的全局解 $\mathbf{x}(t)$ （引用 Tani-Nishiyama 的唯一性结果）满足

\operatorname{dist}(\mathbf{x}(t),\Sigma) \le C\delta, \quad t \ge 0. \qquad (1.4)

此外，解在全局意义下保持在 $\mathbf{x}^R(t) + \mathbf{x}_{0,\text{trans}}$ 的 $O(\delta t)$ 邻域内，即平移模满足

|P_0(\mathbf{x}(t) - \mathbf{x}^R(t) - \mathbf{x}_{0,\text{trans}})| \le C\delta t, \quad t \ge 0. \qquad (1.5)

定理条件解释：

结构条件（流体脉冲接近）： $P[\mathbf{x}_0]$ 是流体的向量脉冲（vector fluid impulse），定义为 $P[\mathbf{x}(t)] = \frac{1}{2}\int_{\mathbb{T}} \mathbf{x}(s,t) \times \mathbf{x}_s(s,t)\,ds.$ 该量在 LIE 中守恒。对于圆环 $\mathbf{x}^R_0$ ，计算得 $P[\mathbf{x}^R_0] = \pi R^2\mathbf{e}_3$ 。条件 $|P[\mathbf{x}_0] - \pi R^2\mathbf{e}_3| \le \delta^2$ 要求初值的脉冲与圆环的脉冲充分接近（小量为 $\delta^2$ ）。
高频模小ness： $P_{|n|\ge2}\mathbf{x}_0$ 是频率 $|n|\ge2$ 的傅里叶投影。条件 $\|P_{|n|\ge2}\mathbf{x}_0\|_{H^2} \le \delta$ 要求高频扰动在 $$H^2$$ 范数下很小。这排除了高频的显著扭曲。
弧长保持： $|\mathbf{x}_{0s}|=1$ 确保参数化为弧长，由 LIE 的几何结构，该条件在时间上保持。
全局解：引用 Tani-Nishiyama 的结果，初值 $\mathbf{x}_0 \in H^4$ 保证唯一全局解存在。

定理结论解释：

轨道稳定性（1.4）：解曲线 $\mathbf{x}(t)$ 到流形 $\Sigma$ （所有圆环模旋转和平移）的 $$H^2$$ 距离全局有界，且界与初值小ness 成正比。这捕捉了形状稳定性：无论平移和旋转如何，曲线本身的弯曲结构保持接近一个圆。
平移模估计（1.5）： $$P_0$$ 是零频投影（即空间平均），对应曲线的整体平移。 $\mathbf{x}_{0,\text{trans}}$ 是初值的平移部分（见 (1.2) 中对 $\boldsymbol{\tau}$ 的最优选择）。该式表明平移模随时间线性增长，但增长率为 $O(\delta t)$ 。这量化了 Lyapunov 不稳定性：平移模式确实会线性漂移，但漂移幅度受 $\delta$ 控制。注意轨道距离（1.4）排除了这种平移漂移的影响。

备注（Remark 1.2） 强调：

圆形涡丝是 Lyapunov 不稳定的（因零模线性增长），但轨道稳定（形状稳定）。
流体脉冲 $P[\mathbf{x}]$ 是关键的几何守恒量（Arms-Hama, Fukumoto）。
初值正则性 $$H^4$$ 用于解唯一性；弱解唯一性在 $$H^2$$ 仍开放。
与线性稳定性研究（如 Calini-Keith-Lafortune 用 Hasimoto 映射）不同，本文处理非线性、一般扰动。
在弹性曲线框架（Langer-Perline）中，圆是弯曲能临界点，但本文利用向量脉冲证明其是（约束）局部极小，这控制±1模。

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

本节的证明围绕将扰动分解为不同频段，并分别估计。设解 $\mathbf{x}(t) = \mathbf{x}^R(t) + \boldsymbol{\phi}(t)$ ，其中 $\boldsymbol{\phi}$ 是扰动。LIE 代入后得扰动方程

\boldsymbol{\phi}_t = A\boldsymbol{\phi} + N\boldsymbol{\phi},

其中线性算子 $A\boldsymbol{\phi} = \mathbf{x}^R_s \times \boldsymbol{\phi} + \boldsymbol{\phi}_s \times \mathbf{x}^R_s$ ，非线性项 $N\boldsymbol{\phi} = \boldsymbol{\phi}_s \times \boldsymbol{\phi}$ 。圆环 $\mathbf{x}^R$ 本身满足 $\mathbf{x}^R_t = \mathbf{x}^R_s \times \mathbf{x}^R_{ss}$ 。关键观察： $$A$$ 具有频移性质， $$A$$ 作用在第 $$n$$ 个傅里叶模会将频率平移 $\pm1$ 。而圆环仅含 $\pm1$ 频率（从 (2.10) 可见）。

傅里叶分解与模态划分 将扰动 $\boldsymbol{\phi}$ 按傅里叶投影分解：

$g = P_0\boldsymbol{\phi}$ ：零模（平移模）。
$\boldsymbol{\psi} = P_{\pm1}\boldsymbol{\phi}$ ：±1 模，对应圆环的相位和幅度调制（即形状变化）。
$\tilde{\boldsymbol{\phi}} = P_{|n|\ge2}\boldsymbol{\phi}$ ：高频模（ $|n|\ge2$ ），代表高波数扭曲。

这种分解是谱分析的标准工具，利用了圆环解仅含低频（±1）的事实。轨道距离 $\operatorname{dist}(\mathbf{x}(t),\Sigma)$ 的计算（见 (2.2)）显示，它等价于 $\|\boldsymbol{\psi}(t) - (Q_\alpha\mathbf{x}^R_0 - \mathbf{x}^R_0)\|_{H^2}$ 与 $\|\tilde{\boldsymbol{\phi}}(t)\|_{H^2}$ 的最小化。其中零模 $$g$$ 被平移 $\boldsymbol{\tau}$ 消除。因此稳定性归结为：证明 $\boldsymbol{\psi}$ 保持接近某个旋转后的圆环，且 $\tilde{\boldsymbol{\phi}}$ 全局小。

核心分析工具：1. 相对能量与高频模控制 定义相对能量（也称弯曲能量）为

E[\boldsymbol{\phi}(t)] = \|\boldsymbol{\phi}_s(t)\|^2 - \frac{1}{R^2}\|\boldsymbol{\phi}(t)\|^2.

它由 LIE 的几何结构守恒（Lemma 2.1(iii)）。关键引理 2.2 表明：若 $P_{\le1}\boldsymbol{\phi}=0$ （即无零模和±1模），则 $E[\boldsymbol{\phi}]$ 范数 $\|\boldsymbol{\phi}\|_E = \sqrt{\langle \boldsymbol{\phi},\boldsymbol{\phi}\rangle_E}$ 等价于 $$H^2$$ 范数。这里双线性形式为

\langle u,v\rangle_E = \int_{\mathbb{T}} \left( u_s \cdot v_s - \frac{1}{R^2} u \cdot v \right) ds.

等价性通过 Parseval 恒等式直接验证：对 $$u$$ 满足 $P_{\le1}u=0$ ，有

\|u\|_E^2 \approx \sum_{|n|\ge2} n^2(n^2-1)|u_n|^2/R^4, \quad \|u\|_2^2 = \sum_{|n|\ge2} |u_n|^2,

且对 $|n|\ge2$ ， $n^2(n^2-1) \asymp n^4$ ，故 $\|u\|_E \approx \|u\|_2$ 。由此，相对能量守恒直接推出 Corollary 2.3：

\|\tilde{\boldsymbol{\phi}}(t)\|_2 \approx \|\tilde{\boldsymbol{\phi}}(0)\|_2 \le \delta, \quad t\ge0.

这是高频模全局小ness 的关键，依赖于初值高频 $$H^2$$ 小ness 条件。

核心分析工具：2. 几何稳定性引理（控制±1模） 相对能量在±1模上退化（因为 $$n^2(n^2-1)=0$$ 当 $$|n|=1$$ ），故不能控制 $\boldsymbol{\psi}$ 。论文核心创新是引理 2.6，它利用流体脉冲守恒 $P[\mathbf{x}(t)] = P[\mathbf{x}_0]$ 导出对±1模的结构性约束。

将 $\mathbf{x}_1$ （ $\mathbf{x}$ 的第一个分量）分解为复系数形式。设 $\mathbf{x}^R_0$ 的傅里叶系数为 $C_1 = R(\mathbf{e}_1 - i\mathbf{e}_2)/2$ ， $C_{-1}=\bar{C}_1$ ，且 $$C_0=0$$ 。对扰动解，由 $P_{\pm1}\mathbf{x} = \mathbf{x}_1$ （因圆环无其他低频），将 $\mathbf{x}_1$ 分解为复基 $\{C_1, \bar{C}_1, |C_1|\mathbf{e}_3\}$ 下的分量：

\mathbf{x}_1(t) = \lambda(t) C_1 + \mu(t) \bar{C}_1 + \nu(t) |C_1|\mathbf{e}_3.

系数 $\lambda,\mu,\nu \in \mathbb{C}$ 为时间函数，且由 $|\mathbf{x}_{1s}|=1$ （弧长条件）和 Parseval 恒等式，有关系

2|\mathbf{x}_1|^2 - 2|C_1|^2 = -\sum_{|n|\ge2} n^2 |\tilde{\phi}_n|^2 \le 0,

结合高频小ness (2.7)，得

\big||\lambda|^2+|\mu|^2+|\nu|^2 - 1\big| \lesssim \delta^2. \quad (2.22)

这是幅度约束。

流体脉冲守恒的作用：计算 $P[\mathbf{x}]$ 的±1模部分。由 LIE 的脉冲守恒及初值接近条件，有

\big| -i\pi \sum_{|n|=1} n\,\mathbf{x}_n \times \mathbf{x}_{-n} - \pi R^2\mathbf{e}_3 \big| \lesssim \delta^2.

利用基关系 $C_1 \times \bar{C}_1 = i|C_1|^2\mathbf{e}_3$ （见 (2.20)），直接计算得

\sum_{|n|=1} n\,\mathbf{x}_n \times \mathbf{x}_{-n} = 2\,\Im(\mathbf{x}_1 \times \mathbf{x}_1) = 2\pi R^2(|\lambda|^2 - |\mu|^2)\mathbf{e}_3 + \text{其他分量}.

但脉冲是向量，其 $\mathbf{e}_3$ 分量与 $|\lambda|^2-|\mu|^2$ 相关。实际上，完整计算给出（见 (2.23)-(2.24)）

\big| |\lambda|^2 - |\mu|^2 - 1 \big| \lesssim \delta^2.

联立 (2.22) 即得引理 2.6：

\big||\lambda|-1\big| + |\mu| + |\nu| \lesssim \delta.

这强烈约束了 $\boldsymbol{\psi}$ ： $\lambda$ 接近 1， $\mu$ 和 $\nu$ 小。注意 $\lambda$ 控制圆环的缩放（ $|\mathbf{x}_1| = |C_1|\sqrt{|\lambda|^2+|\mu|^2+|\nu|^2} \approx R/\sqrt{2}$ ）， $\mu$ 是共轭模式幅度， $\nu$ 对应沿 $\mathbf{e}_3$ 的分量（即圆环平面倾角或非平面性）。引理表明：在脉冲守恒下，圆环的形状（由这三个参数刻画）被锁定在接近原始圆的状态，偏差 $O(\delta)$ 。

几何稳定性引理的证明技巧：利用脉冲守恒将向量方程投影到 $\mathbf{e}_3$ 和其他方向，结合基的正交关系。这是论文的“关键成分”，它提供了相对能量未能覆盖的低频模估计。

±1 模轨道距离的估计 轨道距离等价于 $\inf_\alpha \| \mathbf{x}^R_0 + \boldsymbol{\psi} - Q_\alpha\mathbf{x}^R_0 \|_{H^2}$ 。由于 $Q_\alpha\mathbf{x}^R_0(s) = \mathbf{x}^R_0(s+R\alpha)$ ，这相当于寻找相位平移 $\sigma=R\alpha$ 使 $\|\mathbf{x}_1 - \mathbf{x}^R_1(\cdot+\sigma)\|$ 最小。设 $\mathbf{x}^R_1(s) = C_1 e^{is/R} + \bar{C}_1 e^{-is/R}$ 。最小化问题归结为复数问题：给定 $a=\mathbf{x}_1$ ，找 $\sigma$ 最小化 $|a - e^{i\sigma/R} C_1|$ 。引理 2.4 给出最小值表达式：

\min_\theta |a - e^{i\theta}C_1|^2 = (|a|-|C_1|)^2 + 2\big(|a||C_1| - |a\cdot C_1|\big),

其中点积为复内积。定义幅度误差 $\varepsilon_\parallel = \big||\mathbf{x}_1|-|C_1|\big|$ 和角度误差 $\varepsilon_\perp = \sqrt{2(| \mathbf{x}_1||C_1| - |\mathbf{x}_1\cdot C_1|)}$ ，则距离平方 $\approx \varepsilon_\parallel^2 + \varepsilon_\perp^2$ 。

利用引理 2.6 和关系 $$|C_1|^2 = R^2/2$$ ，可证：

幅度误差 $\varepsilon_\parallel \lesssim \|\tilde{\boldsymbol{\phi}}\|_2^2 \lesssim \delta^2$ （来自 (2.16), (2.17)）。
角度误差 $\varepsilon_\perp \lesssim \delta^2$ （来自引理 2.6 对 $\mu,\nu$ 的约束，见 (2.25)）。因此 $\operatorname{dist}_{\pm1}(t) \lesssim \delta$ ，结合高频估计即得 (2.26) $\inf_\alpha \| \mathbf{x}^R_0+\boldsymbol{\psi} - Q_\alpha\mathbf{x}^R_0 \|_{H^2} \lesssim \delta$ ，最终得到轨道稳定性 (1.4)。

核心分析工具：3. 零模（平移模）的估计 零模 $g = P_0\boldsymbol{\phi}$ 由方程 (2.6) 控制：

\partial_t g = P_0\big(A(\boldsymbol{\psi}+\tilde{\boldsymbol{\phi}}) + N(\boldsymbol{\psi}+\tilde{\boldsymbol{\phi}})\big).

直接计算（利用 $\mathbf{x}^R$ 的显式）得 $P_0 A\boldsymbol{\psi} = \frac{2i}{R^3}\big( -\bar{C}_1 \times \boldsymbol{\psi}_1 + C_1 \times \boldsymbol{\psi}_{-1} \big)$ （见 (2.27)），这是 $\boldsymbol{\psi}$ 的线性耦合。非线性项 $P_0 N(\cdot)$ 涉及平方（见 (2.29)）。由于 $\boldsymbol{\psi}$ 有界（引理 2.7 给出 $\|\boldsymbol{\psi}\| \approx \sqrt{|\lambda-1|^2+|\mu|^2+|\nu|^2}$ ，而引理 2.6 仅保证 $\mu,\nu$ 小，但 $|\lambda-1|$ 无界——因初值未要求 $\lambda(0)$ 接近 1，故 $\|\boldsymbol{\psi}\|$ 可能 $$O(1)$$ ），故 $P_0 A\boldsymbol{\psi}$ 有界。而 $P_0 A\tilde{\boldsymbol{\phi}} \lesssim \|\tilde{\boldsymbol{\phi}}\|_2 \le \delta$ （由 (2.28)）。非线性项 $$P_0 N$$ 受 $\|\boldsymbol{\psi}\|_\infty^2$ 和 $\|\tilde{\boldsymbol{\phi}}\|_2\|\boldsymbol{\psi}\|_\infty$ 控制，但 $$H^2$$ 嵌入 $L^\infty$ 需低维，在 $\mathbb{T}$ 上 $H^2 \subset C^{1,1/2} \subset L^\infty$ ，故有界。

实际上，对 $$g$$ 的估计不需要它的有界性，因为定理结论 (1.5) 涉及 $P_0(\mathbf{x} - \mathbf{x}^R - \mathbf{x}_{0,\text{trans}})$ 。写出 $g(t) - g(0) = \int_0^t \partial_t g\,d\tau$ ，由 (2.6) 及上述估计， $\|\partial_t g\|$ 有界（因 $\boldsymbol{\psi},\tilde{\boldsymbol{\phi}}$ 有界和小），故 $\|g(t) - g(0)\| \lesssim t$ 。但注意 $$g(0)$$ 是初值零模，且平移 $\boldsymbol{\tau}$ 的最优选择为 $\boldsymbol{\tau}(t) = g(t) + (t/R)\mathbf{e}_3$ （因 $\mathbf{x}^R_t$ 贡献 $(t/R)\mathbf{e}_3$ 到零模）。故 (1.5) 正是 $$g(t) - g(0)$$ 的估计（忽略高阶小量），增长率 $O(\delta t)$ 来自 $\|\boldsymbol{\psi}\|$ 可能 $$O(1)$$ ，但常数与 $\delta$ 无关（因 $\delta$ 仅控制高频和脉冲，不迫使 $\lambda$ 接近 1）。

关键技术技巧总结：

模态分离与相对能量：利用相对能量 $E[\boldsymbol{\phi}]$ 守恒，通过等价引理 2.2，将高频模 $\tilde{\boldsymbol{\phi}}$ 的 $$H^2$$ 范数锁定在初值水平。这是基于 $n^2(n^2-1) \asymp n^4$ 对 $|n|\ge2$ 的谱 Gap。
几何引理与脉冲守恒：对低频±1模，利用流体脉冲 $P[\mathbf{x}]$ 的守恒，通过将 $\mathbf{x}_1$ 在特定复基下分解，导出 $|\lambda|,|\mu|,|\nu|$ 的联合约束（引理 2.6）。这本质上是利用向量的脉冲守恒（一个向量方程）控制三个复参数，需精细的基选择（ $C_1,\bar{C}_1,|C_1|\mathbf{e}_3$ ）以匹配圆环对称性。
相位最优化的显式公式：在估计 $\operatorname{dist}_{\pm1}$ 时，将无限维最小化问题归约为复数的相位最小化（引理 2.4），显式表达为幅度差和角度差的平方和。
估计的耦合：幅度误差 $\varepsilon_\parallel$ 用高频小ness 通过弧长条件控制；角度误差 $\varepsilon_\perp$ 用脉冲守恒（即引理 2.6）控制。两者结合给出 $\boldsymbol{\psi}$ 的小ness。
零模的长期估计：不试图证明 $$g$$ 有界，而是估计其线性增长 $$O(t)$$ ，这源于 $\boldsymbol{\psi}$ 可能不衰减（因 $\lambda$ 无约束），但增长率的常数与 $\delta$ 无关（初值小ness 仅通过高频和脉冲间接约束，不直接控制 $$g$$ 的漂移）。

整体框架：证明嵌入在 Tani-Nishiyama 的 LIE 整体解理论中（保证唯一性），并借鉴 Aiki (2025) 对对称扰动的分析，但通过几何引理去除对称性。核心是“互补”策略：相对能量控制高频，脉冲守恒控制低频。

与之前工作的比较

本文直接扩展了第一作者之前的工作 [Aiki, 2025]，并与其他文献（如 Choi 对 Hill 涡的工作、Calini 等对线性稳定性、Langer-Perline 对弹性曲线）有概念联系但技术不同。

主要突破与改进：

去除对称性假设：[Aiki, 2025] 要求扰动满足额外的对称条件，具体是 $\boldsymbol{\psi} \cdot \mathbf{e}_3 = 0$ （即±1模在 $\mathbf{e}_3$ 方向分量为零），这排除了圆环平面的倾斜或非平面扰动。本文定理 1.1 允许任意非对称扰动，只需高频模小ness 和脉冲接近。这是本质放宽：对称条件实质性限制了扰动类（例如，它强制 $\nu(t)=0$ 在基分解中），而去除后证明更一般，但需要新工具（几何引理）来处理释放的 $\nu$ 模式。改进体现在扰动空间的维度：从对称子空间到全空间。
轨道稳定性的严格实现：即使去除对称性，仍证明形状稳定（轨道距离有界）。这要求明确处理平移和旋转对称性，通过定义 $\Sigma$ 和 $\operatorname{dist}(\mathbf{x},\Sigma)$ 。在 [Aiki, 2025] 中，对称性可能自动吸收部分旋转，但这里必须显式最小化。
几何稳定性引理（Lemma 2.6）：这是本文的核心新工具。它利用向量流体脉冲 $P[\mathbf{x}]$ 守恒，约束了三个低频参数 $|\lambda|-1, |\mu|, |\nu|$ 。在 [Aiki, 2025] 中，对称条件 $\boldsymbol{\psi}\cdot\mathbf{e}_3=0$ 导致 $\nu=0$ ，且脉冲条件自动给出 $|\lambda|^2-|\mu|^2=1$ 结合弧长条件 $|\lambda|^2+|\mu|^2=1$ 直接推出 $|\lambda|=1, |\mu|=0$ ，即 $\boldsymbol{\psi}=0$ 。因此不需要引理 2.6 这样的强度估计。本文中 $\nu$ 自由，脉冲条件 (2.24) 与弧长条件 (2.22) 仅能推出三个参数的小ness，而非零，但足够控制轨道距离。该引理的证明通过将 $\mathbf{x}_1$ 分解到与圆环傅里叶基对齐的复基上，并利用基的正交性和叉积关系（(2.19)-(2.20)），将向量脉冲的 $\mathbf{e}_3$ 分量提取为 $|\lambda|^2-|\mu|^2$ 。这是一种几何演算，将高维守恒量投影到低维参数空间。
高频模小ness 条件：[Aiki, 2025] 可能也使用了类似条件（因为相对能量控制高频是标准）。但本文明确将 $\|P_{|n|\ge2}\mathbf{x}_0\|_{H^2} \le \delta$ 作为假设，这确保了 $\tilde{\boldsymbol{\phi}}$ 全局小。注意到该条件比 $E[\boldsymbol{\phi}(0)] \le \delta^2$ 更强（因 $$E$$ 等价于 $\|\tilde{\boldsymbol{\phi}}\|_2^2$ 仅当 $\boldsymbol{\psi}=g=0$ ），但本文中 $\boldsymbol{\psi}$ 不假设为零，故需直接假设高频 $$H^2$$ 小。这是与对称情形不同的技术点：对称情形下，脉冲条件和弧长条件可能迫使 $\boldsymbol{\psi}=0$ ，从而 $$E$$ 自动控制高频；但非对称下 $\boldsymbol{\psi}$ 可能非零，故需直接控制高频。因此，本文的结果在扰动类上更广，但多了一个小ness 假设（高频模小），这是允许非对称扰动的代价。
$\boldsymbol{\psi}$ 有界但不衰减：引理 2.7 显示 $\|\boldsymbol{\psi}\| \approx \sqrt{|\lambda-1|^2+|\mu|^2+|\nu|^2}$ 。引理 2.6 仅约束 $|\mu|,|\nu|$ 小，但 $|\lambda|$ 可以偏离 1（只要 $|\lambda|^2 \approx 1+|\mu|^2+|\nu|^2$ ）。因此 $\|\boldsymbol{\psi}\|$ 可能 $$O(1)$$ 。这意味着 $\boldsymbol{\psi}$ 不一定衰减，但轨道距离 $\inf_\alpha \| \mathbf{x}^R_0+\boldsymbol{\psi} - Q_\alpha\mathbf{x}^R_0\|_{H^2}$ 仍小，因为偏离可通过对 $\alpha$ 的最优化（即旋转对齐）吸收幅度 $|\lambda|$ 的变化（见 (2.15) 中 $\varepsilon_\parallel$ 小）。这体现了轨道稳定性的微妙：形状（相对几何）稳定，但绝对相位（ $\lambda$ 对应缩放？实际 $\lambda$ 影响 $$|x_1|$$ ，但圆环半径固定，故 $\lambda$ 偏离意味着椭圆化？注意圆环的 $\lambda=\mu=1, \nu=0$ ？由 (2.10) $\mathbf{x}^R_0$ 对应 $$C_1$$ ，故 $\lambda=1, \mu=0, \nu=0$ 。若 $\lambda$ 接近 1 但 $\mu$ 小，则 $\mathbf{x}_1 \approx \lambda C_1$ ，这是圆环的缩放，但旋转后 $Q_\alpha\mathbf{x}^R_0$ 的 $$x_1$$ 为 $e^{i\alpha}C_1$ ，故模长 $$|C_1|$$ ，因此 $\lambda$ 必须接近 1，引理 2.6 给出 $|\lambda|-1$ 小，所以 $\|\boldsymbol{\psi}\|$ 实际也小？更正：引理 2.6 说 $|\lambda|-1$ 小，故 $\|\boldsymbol{\psi}\| \approx \sqrt{|\lambda-1|^2+|\mu|^2+|\nu|^2} \lesssim \delta$ 。因此 $\boldsymbol{\psi}$ 确实 $O(\delta)$ 。我之前的疑虑有误：引理 2.6 给出 $|\lambda|-1$ 小，故 $\lambda$ 接近 1，因此 $\|\boldsymbol{\psi}\| \lesssim \delta$ 。是的，三个项都小，所以 $\boldsymbol{\psi}$ 小。那为何引理 2.7 说 $\|\boldsymbol{\psi}\| \approx \sqrt{|\lambda-1|^2+|\mu|^2+|\nu|^2}$ ？因为原参考 $\mathbf{x}^R_0$ 对应 $\lambda=1, \mu=0,\nu=0$ ，所以 $\boldsymbol{\psi} = (\lambda-1)C_1 + \mu\bar{C}_1 + \nu|C_1|\mathbf{e}_3$ ，其范数正比于这些系数。所以 $\|\boldsymbol{\psi}\| \lesssim \delta$ 。因此 $\boldsymbol{\psi}$ 确实衰减？但定理 1.1 的轨道距离是 $\inf_\alpha \| \mathbf{x}^R_0+\boldsymbol{\psi} - Q_\alpha\mathbf{x}^R_0 \|_{H^2}$ ，由 (2.26) 也 $\lesssim \delta$ 。但注意 $\boldsymbol{\psi}$ 本身是 $$x$$ 的傅里叶±1模，它可能不衰减，但旋转后对齐则小。实际上，若 $\boldsymbol{\psi}$ 小，则存在 $\alpha$ 使 $Q_\alpha\mathbf{x}^R_0$ 接近 $\mathbf{x}^R_0+\boldsymbol{\psi}$ 。因此， $\boldsymbol{\psi}$ 确实 $O(\delta)$ 全局。这是从引理 2.6 和 (2.16) 推出的？(2.16) 给出 $2|\mathbf{x}_1|^2 - 2|C_1|^2 = -\sum_{|n|\ge2} n^2|\tilde{\phi}_n|^2$ ，故 $|\mathbf{x}_1| \approx |C_1|$ ，但 $|\mathbf{x}_1| = |C_1|\sqrt{|\lambda|^2+|\mu|^2+|\nu|^2}$ ，所以 $\sqrt{|\lambda|^2+|\mu|^2+|\nu|^2} \approx 1$ ，但 $\lambda$ 可以接近 -1？由 $|\lambda|^2+|\mu|^2+|\nu|^2 \approx 1$ ，仅保证 $|\lambda| \le 1+O(\delta)$ ，但未紧约束 $\lambda$ 接近 1。例如若 $\mu=\nu=0$ ，则 $|\lambda| \approx 1$ ，但 $\lambda$ 可为负？实系数？ $\lambda$ 是复，但 $|\lambda| \approx 1$ 。圆环参考对应 $\lambda=1$ 。若 $\lambda = -1$ （反相），则 $\mathbf{x}_1 \approx -C_1 e^{i\theta}$ ，旋转 $\alpha = \pi$ 可对齐（因 $e^{i\pi}C_1 = -C_1$ ），所以若 $|\lambda|=1$ ，仍可通过旋转对齐，只要 $\mu=\nu=0$ 。但若 $\lambda$ 复数？ $$C_1$$ 固定。最小化 $|\mathbf{x}_1 - e^{i\theta}C_1|$ 时，最优 $\theta = \arg(\mathbf{x}_1 \cdot C_1)$ ，最小值取决于 $|\mathbf{x}_1|$ 和 $|\mathbf{x}_1\cdot C_1|$ 。由 $\mathbf{x}_1 = \lambda C_1 + \mu\bar{C}_1 + \cdots$ ， $\mathbf{x}_1\cdot C_1 = \lambda|C_1|^2$ （因 $C_1\cdot\bar{C}_1=0$ ），故 $|\mathbf{x}_1\cdot C_1| = |\lambda| |C_1|^2$ 。于是 $\varepsilon_\parallel = \big||\mathbf{x}_1|-|C_1|\big| = |C_1|\big|\sqrt{|\lambda|^2+|\mu|^2+|\nu|^2}-1\big| \lesssim \delta^2$ 由 (2.22) 推得因为左边 $\le$ 右边？实际上 (2.22) 是 $|\,|\lambda|^2+|\mu|^2+|\nu|^2 - 1\,| \lesssim \delta^2$ ，故 $\big|\sqrt{|\lambda|^2+|\mu|^2+|\nu|^2}-1\big| \lesssim \delta^2$ （因平方根 Lip），所以 $\varepsilon_\parallel \lesssim \delta^2$ 。 $\varepsilon_\perp = \sqrt{2}(| \mathbf{x}_1||C_1| - |\mathbf{x}_1\cdot C_1|) = \sqrt{2}|C_1|^2\big(\sqrt{|\lambda|^2+|\mu|^2+|\nu|^2} - |\lambda|\big)$ 。由 $|\lambda|^2+|\mu|^2+|\nu|^2 - |\lambda|^2 = |\mu|^2+|\nu|^2 \lesssim \delta^2$ ，且 $\sqrt{a+b} - \sqrt{a} \le b/(2\sqrt{a})$ ，对 $a=|\lambda|^2 \ge c>0$ （由 $|\lambda|^2 \ge 1-|\mu|^2-|\nu|^2 \ge 1-O(\delta^2)$ ，故对 $\delta$ 小 $|\lambda|$ 远离 0），得 $\varepsilon_\perp \lesssim |C_1|^2 \cdot (|\mu|^2+|\nu|^2)/(2|\lambda|) \lesssim \delta^2$ 。因此 $\boldsymbol{\psi}$ 的幅度本身可能不？注意 $\|\boldsymbol{\psi}\|^2 \approx |\lambda-1|^2+|\mu|^2+|\nu|^2$ 。由 $|\lambda|^2+|\mu|^2+|\nu|^2 \approx 1$ ， $|\lambda| \approx 1$ ，但 $|\lambda-1|$ 可大若 $\lambda$ 为负？例 $\lambda = -1$ ，则 $|\lambda|=1$ ， $|\lambda-1|=2$ ，但此时 $\mu=\nu=0$ ？检查条件 (2.22)： $|\lambda|^2+|\mu|^2+|\nu|^2=1$ ，若 $\lambda=-1$ ，则成立。 (2.24)： $|\lambda|^2-|\mu|^2=1$ 也得 $$1-0=1$$ ，成立。所以 $\lambda=-1$ 允许？但 $\lambda$ 是复数， $|\lambda|^2$ 是模平方。 $\lambda=-1$ 实， $|\lambda|^2=1$ 。此时 $\boldsymbol{\psi} = (-1-1)C_1 = -2C_1$ ，范数 $$O(1)$$ 。但轨道距离： $\mathbf{x}_1 = -C_1$ ，旋转 $\alpha=\pi$ 得 $Q_\pi\mathbf{x}^R_0$ 的 $$x_1$$ 为 $e^{i\pi}C_1 = -C_1$ ，完全匹配，故轨道距离为零？是的，因为 $\mathbf{x}^R_0$ 旋转 $\pi$ 后为 $-\mathbf{x}^R_0$ （在 $$xy$$ 平面），但注意 $\mathbf{x}^R_0(s) = R\cos(s/R)\mathbf{e}_1 + R\sin(s/R)\mathbf{e}_2$ ，旋转 $\pi$ 得 $R\cos(s/R+\pi)\mathbf{e}_1 + R\sin(s/R+\pi)\mathbf{e}_2 = -R\cos(s/R)\mathbf{e}_1 - R\sin(s/R)\mathbf{e}_2 = -\mathbf{x}^R_0(s)$ ，但曲线本身相同（圆），所以确实属于 $\Sigma$ 。因此 $\lambda=-1$ 对应同一圆环（反相位），轨道距离为零。但 $\boldsymbol{\psi}$ 本身大（ $$-2C_1$$ ），但通过旋转消除。所以轨道距离小不要求 $\boldsymbol{\psi}$ 小，只要求存在旋转使其小。在 $\lambda=-1$ 时， $\boldsymbol{\psi} = -2C_1$ ，但旋转 $\pi$ 后差值 $(\mathbf{x}^R_0+\boldsymbol{\psi}) - Q_\pi\mathbf{x}^R_0 = \mathbf{x}^R_0 -2C_1 - (-\mathbf{x}^R_0) = 2(\mathbf{x}^R_0 - C_1)$ ？等等需仔细： $Q_\pi\mathbf{x}^R_0(s) = \mathbf{x}^R_0(s+\pi R) = R\cos((s+\pi R)/R)\mathbf{e}_1 + \cdots = R\cos(s/R+\pi)\mathbf{e}_1 + R\sin(s/R+\pi)\mathbf{e}_2 = -\mathbf{x}^R_0(s)$ 。 $\mathbf{x}^R_0+\boldsymbol{\psi} = \mathbf{x}^R_0 -2C_1$ 。 $C_1 = R(\mathbf{e}_1-i\mathbf{e}_2)/2$ ，而 $\mathbf{x}^R_0 = C_1 e^{is/R} + \bar{C}_1 e^{-is/R}$ 。差值 $\mathbf{x}^R_0 -2C_1 = C_1(e^{is/R}-2) + \bar{C}_1 e^{-is/R}$ 。这与 $-\mathbf{x}^R_0 = -C_1 e^{is/R} - \bar{C}_1 e^{-is/R}$ 不同。因此需计算 $\inf_\alpha \| \mathbf{x}^R_0+\boldsymbol{\psi} - Q_\alpha\mathbf{x}^R_0 \|$ 。取 $\alpha$ 使 $Q_\alpha\mathbf{x}^R_0 = \mathbf{x}^R_0(\cdot+\sigma)$ 。傅里叶系数： $\mathbf{x}^R_0$ 的 $$n=1$$ 模是 $$C_1$$ ， $$n=-1$$ 是 $\bar{C}_1$ 。 $Q_\alpha\mathbf{x}^R_0$ 的 $$n=1$$ 模是 $e^{i\alpha}C_1$ ， $$n=-1$$ 是 $e^{-i\alpha}\bar{C}_1$ 。 $\mathbf{x}^R_0+\boldsymbol{\psi}$ 的 $$n=1$$ 模是 $C_1 + \psi_1$ ， $$n=-1$$ 模是 $\bar{C}_1 + \psi_{-1}$ 。这里 $\psi_1 = \lambda C_1 + \mu\bar{C}_1 + \nu|C_1|\mathbf{e}_3$ ？但 $\boldsymbol{\psi}$ 仅含±1模，且 $\mathbf{e}_3$ 方向？实际上 $\boldsymbol{\psi} = P_{\pm1}\boldsymbol{\phi}$ ，而 $\mathbf{x}^R_0$ 的±1模完全在 $\mathrm{span}\{C_1,\bar{C}_1\}$ （ $$xy$$ 平面），但如果 $\boldsymbol{\psi}$ 有 $\mathbf{e}_3$ 分量，那来自 $\nu|C_1|\mathbf{e}_3$ ，但 $\mathbf{e}_3$ 是 $$x_3$$ 方向，而 $$C_1$$ 在 $$xy$$ 平面，所以 $\boldsymbol{\psi}$ 可能脱离圆环平面。但 $Q_\alpha\mathbf{x}^R_0$ 仍在 $$xy$$ 平面。因此若 $\boldsymbol{\psi}$ 有 $\mathbf{e}_3$ 分量，轨道距离不可能为零。在 $\lambda=-1,\mu=\nu=0$ 时， $\boldsymbol{\psi} = -2C_1$ 在 $$xy$$ 平面，无 $\mathbf{e}_3$ 分量。此时 $\mathbf{x}_1$ 的 $$n=1$$ 模为 $$C_1 + (-2C_1) = -C_1$$ ， $$n=-1$$ 模为 $\bar{C}_1$ 。选 $\alpha=\pi$ ，则 $Q_\alpha\mathbf{x}^R_0$ 的 $$n=1$$ 模为 $e^{i\pi}C_1 = -C_1$ ， $$n=-1$$ 为 $e^{-i\pi}\bar{C}_1 = -\bar{C}_1$ 。但 $\mathbf{x}^R_0+\boldsymbol{\psi}$ 的 $$n=-1$$ 模是 $\bar{C}_1$ ，不匹配。所以需要同时匹配正负模。实际上 $\boldsymbol{\psi}$ 必须满足 $\psi_{-1} = \bar{\psi}_1$ （因 $\boldsymbol{\phi}$ 实值，傅里叶系数共轭对称）。这里 $\psi_1 = -2C_1$ ，则 $\psi_{-1} = -2\bar{C}_1$ ？若 $\boldsymbol{\phi}$ 实，则 $(\boldsymbol{\phi}_n)^* = \boldsymbol{\phi}_{-n}$ 。故若 $\psi_1 = -2C_1$ ，则 $\psi_{-1} = -2\bar{C}_1$ 。但上面我写 $\psi_{-1} = \bar{C}_1$ 是错的。必须自洽：设 $\boldsymbol{\psi} = P_{\pm1}\boldsymbol{\phi}$ ，则 $\psi_{-1} = \overline{\psi_1}$ 。若 $\boldsymbol{\psi}$ 仅由 $\lambda$ 控制？在分解 $\mathbf{x}_1 = \lambda C_1 + \mu\bar{C}_1 + \nu|C_1|\mathbf{e}_3$ 中， $\mathbf{x}_1$ 是实向量？不， $\mathbf{x}_1$ 是 $\mathbf{x}$ 的第一分量，是实值函数。但 $C_1,\bar{C}_1$ 是复基。系数 $\lambda,\mu,\nu$ 是复数，但 $\mathbf{x}_1(s) = \lambda C_1 e^{is/R} + \mu\bar{C}_1 e^{-is/R} + \nu|C_1|\mathbf{e}_3$ 。为使 $\mathbf{x}_1$ 实值，必须 $\mu = \bar{\lambda}$ （因为 $C_1 e^{i\theta} + \bar{C}_1 e^{-i\theta}$ 实当系数共轭）。实际上，对实函数， $x_{1,n} = \overline{x_{1,-n}}$ 。 $$n=1$$ 模 $x_{1,1}$ 对应 $$C_1$$ 的系数？根据 (2.10)， $\mathbf{x}^R_0$ 的 $$n=1$$ 模是 $$C_1$$ ， $$n=-1$$ 模是 $\bar{C}_1$ 。所以一般实函数 $\mathbf{x}$ ，设 $x_{1,1} = a$ ，则 $x_{1,-1} = \bar{a}$ 。那么 $\mathbf{x}_1(s) = a e^{is/R} + \bar{a} e^{-is/R} + \sum_{|n|\ge2} x_{1,n} e^{ins/R}$ 。对比 $\lambda C_1 e^{is/R} + \mu\bar{C}_1 e^{-is/R}$ ，因 $$C_1$$ 固定，有 $\lambda C_1 = a$ ， $\mu\bar{C}_1 = \bar{a}$ ，故 $\mu = \bar{\lambda} \bar{C}_1^{-1} C_1$ ？更简单： $C_1 = \frac{R}{2}(\mathbf{e}_1 - i\mathbf{e}_2)$ ，其共轭 $\bar{C}_1 = \frac{R}{2}(\mathbf{e}_1 + i\mathbf{e}_2)$ 。注意 $\bar{C}_1$ 与 $$C_1$$ 线性无关（在 $\mathbb{C}^3$ ）。对实 $\mathbf{x}_1$ ，系数满足 $x_{1,-1} = \overline{x_{1,1}}$ 。设 $x_{1,1} = \alpha$ （复），则 $x_{1,-1} = \bar{\alpha}$ 。那么 $\mathbf{x}_1$ 的 $$n=1$$ 部分为 $\alpha e^{is/R}$ ， $$n=-1$$ 部分为 $\bar{\alpha} e^{-is/R}$ 。另一方面， $\lambda C_1 e^{is/R} + \mu\bar{C}_1 e^{-is/R}$ 。比较得：

$$n=1$$ 系数： $\lambda C_1 = \alpha$
$$n=-1$$ 系数： $\mu\bar{C}_1 = \bar{\alpha}$ 故 $\mu\bar{C}_1 = \overline{\lambda C_1} = \bar{\lambda} \bar{C}_1$ ，因此 $\mu = \bar{\lambda}$ （因 $\bar{C}_1 \neq 0$ ）。所以实际上 $\mu = \bar{\lambda}$ 。然后 $\nu$ 是 $\mathbf{e}_3$ 方向的实系数？因为 $\mathbf{e}_3$ 是实向量， $|C_1|\mathbf{e}_3$ 也是实，故 $\nu$ 应为实数（使 $\nu|C_1|\mathbf{e}_3$ 实）。在分解 (2.18) 中， $\lambda(t),\mu(t) \in \mathbb{C}$ 但由实值性必须 $\mu=\bar{\lambda}$ ，且 $\nu(t) \in \mathbb{R}$ 。因此参数化实际是 $\lambda \in \mathbb{C}$ ， $\nu \in \mathbb{R}$ ，且 $\mu=\bar{\lambda}$ 。故 $|\lambda|=|\mu|$ 。此关系未在论文文字中明确写出，但从上下文应隐含（因 $\mathbf{x}$ 实值）。检查：在 (2.18) 中， $\mathbf{x}_1 = \lambda C_1 + \mu C_1 + \nu|C_1|\mathbf{e}_3$ ，但 $$C_1$$ 和 $$C_1$$ 是共轭基，为使 $\mathbf{x}_1$ 实，需 $\mu = \bar{\lambda}$ 。在脉冲计算中， $\mathbf{x}_1 \times \mathbf{x}_1$ 的实部涉及 $|\lambda|^2-|\mu|^2$ ，若 $\mu=\bar{\lambda}$ 则 $|\lambda|^2-|\mu|^2=0$ ，但圆环脉冲有 $\mathbf{e}_3$ 分量 $\pi R^2$ ，这似乎矛盾。实际上， $\mathbf{x}_1 \times \mathbf{x}_1$ 是向量，其 $\mathbf{e}_3$ 分量是 $\Im(\bar{x}_{1,1} x_{1,1})$ 之类？计算： $\mathbf{x}(s) = \sum_n \mathbf{x}_n e^{ins/R}$ ，但 $\mathbf{x}_n \in \mathbb{C}^3$ 。 $P[\mathbf{x}] = -i\pi \sum_n n \mathbf{x}_n \times \mathbf{x}_{-n}$ 。若 $\mathbf{x}$ 实，则 $\mathbf{x}_{-n} = \overline{\mathbf{x}_n}$ 。故 $P[\mathbf{x}] = -i\pi \sum_n n \mathbf{x}_n \times \overline{\mathbf{x}_n}$ ，这是纯虚？但 $P[\mathbf{x}]$ 应为实向量。实际上， $\mathbf{x}_n \times \overline{\mathbf{x}_n}$ 是纯虚（因 $(a+ib)\times(a-ib) = 2i(a\times b)$ 实部？叉积性质：对复向量 $$u,v$$ ， $u\times v$ 通常复，但 $u\times \bar{u}$ 是 $$2i$$ 乘某个实向量？设 $$u = a+ib$$ ， $$a,b$$ 实，则 $\bar{u}=a-ib$ ， $u\times\bar{u} = (a\times a - b\times b) + i(b\times a - a\times b) = 0 + i( -2 a\times b) = -2i (a\times b)$ ，故为纯虚。则 $-i\pi n (u\times\bar{u}) = -i\pi n (-2i a\times b) = -2\pi n (a\times b)$ ，实。正确。因此， $\mathbf{x}_n \times \overline{\mathbf{x}_n}$ 的实部来自 $a\times b$ 。对 $$n=1$$ ， $\mathbf{x}_1 = \mathbf{x}_{1,\text{Re}} + i\mathbf{x}_{1,\text{Im}}$ ，则贡献与 $\mathbf{x}_{1,\text{Re}} \times \mathbf{x}_{1,\text{Im}}$ 有关。在基分解下， $C_1 = \frac{R}{2}(\mathbf{e}_1 - i\mathbf{e}_2)$ ，设 $\lambda = \lambda_R + i\lambda_I$ ，则 $\lambda C_1 = \frac{R}{2}[(\lambda_R\mathbf{e}_1 + \lambda_I\mathbf{e}_2) + i(-\lambda_R\mathbf{e}_2 + \lambda_I\mathbf{e}_1)]$ ，故 $\mathbf{x}_{1,\text{Re}} = \frac{R}{2}(\lambda_R\mathbf{e}_1 + \lambda_I\mathbf{e}_2)$ ， $\mathbf{x}_{1,\text{Im}} = \frac{R}{2}(-\lambda_R\mathbf{e}_2 + \lambda_I\mathbf{e}_1)$ 。则 $\mathbf{x}_{1,\text{Re}} \times \mathbf{x}_{1,\text{Im}} = \frac{R^2}{4}[(\lambda_R\mathbf{e}_1 + \lambda_I\mathbf{e}_2) \times (-\lambda_R\mathbf{e}_2 + \lambda_I\mathbf{e}_1)] = \frac{R^2}{4}[\lambda_R\lambda_I (\mathbf{e}_1\times\mathbf{e}_1) + \lambda_R^2 (\mathbf{e}_1\times(-\mathbf{e}_2)) + \lambda_I^2 (\mathbf{e}_2\times\mathbf{e}_1) + \lambda_R\lambda_I (\mathbf{e}_2\times\mathbf{e}_1)?]$ 仔细计算： $\mathbf{a} = \lambda_R\mathbf{e}_1 + \lambda_I\mathbf{e}_2$ ， $\mathbf{b} = -\lambda_R\mathbf{e}_2 + \lambda_I\mathbf{e}_1$ 。 $\mathbf{a}\times\mathbf{b} = (\lambda_R\mathbf{e}_1 + \lambda_I\mathbf{e}_2) \times (-\lambda_R\mathbf{e}_2 + \lambda_I\mathbf{e}_1) = \lambda_R(-\lambda_R)(\mathbf{e}_1\times\mathbf{e}_2) + \lambda_R\lambda_I(\mathbf{e}_1\times\mathbf{e}_1) + \lambda_I(-\lambda_R)(\mathbf{e}_2\times\mathbf{e}_2) + \lambda_I\lambda_I(\mathbf{e}_2\times\mathbf{e}_1)$ $= -\lambda_R^2 (\mathbf{e}_1\times\mathbf{e}_2) + 0 + 0 + \lambda_I^2 (\mathbf{e}_2\times\mathbf{e}_1)$ 但 $\mathbf{e}_2\times\mathbf{e}_1 = -\mathbf{e}_1\times\mathbf{e}_2 = -\mathbf{e}_3$ ，故 $= -\lambda_R^2 \mathbf{e}_3 - \lambda_I^2 \mathbf{e}_3 = -(\lambda_R^2+\lambda_I^2)\mathbf{e}_3 = -|\lambda|^2\mathbf{e}_3$ 。因此 $\mathbf{x}_1 \times \overline{\mathbf{x}_1} = -2i (\mathbf{x}_{1,\text{Re}} \times \mathbf{x}_{1,\text{Im}}) = -2i \cdot \frac{R^2}{4}(-|\lambda|^2\mathbf{e}_3) = \frac{iR^2}{2}|\lambda|^2\mathbf{e}_3$ 。那么 $-i\pi \cdot 1 \cdot (\mathbf{x}_1 \times \overline{\mathbf{x}_1}) = -i\pi \cdot \frac{iR^2}{2}|\lambda|^2\mathbf{e}_3 = \frac{\pi R^2}{2}|\lambda|^2\mathbf{e}_3$ 。同样对 $$n=-1$$ 项，设 $$n=-1$$ ， $x_{1,-1} = \overline{x_{1,1}} = \bar{\lambda}C_1?$ 实际上由 $\mu=\bar{\lambda}$ ， $$n=-1$$ 模是 $\mu\bar{C}_1 = \bar{\lambda}\bar{C}_1$ 。但 $\bar{C}_1 = \frac{R}{2}(\mathbf{e}_1 + i\mathbf{e}_2)$ 。计算类似， $x_{1,-1} \times x_{1,1}$ （注意 $P[\mathbf{x}]$ 中项为 $n\,\mathbf{x}_n \times \mathbf{x}_{-n}$ ，对 $$n=-1$$ 为 $(-1) x_{1,-1} \times x_{1,1}$ ）。整体 $$n=1$$ 和 $$n=-1$$ 贡献相加。对称性下，总 $P[\mathbf{x}]$ 的 $\mathbf{e}_3$ 分量为 $\pi R^2(|\lambda|^2)$ ？对圆环 $\lambda=1$ 得 $\pi R^2$ 。这里我们发现仅 $|\lambda|$ 出现， $|\mu|$ 不出现（因 $\mu=\bar{\lambda}$ ， $|\mu|=|\lambda|$ ）。但引理 2.6 中有 $|\lambda|^2 - |\mu|^2$ ，若 $\mu=\bar{\lambda}$ 则 $|\lambda|^2-|\mu|^2=0$ ，矛盾。显然我误解了基分解。重新检查论文（第8页）： (2.18) 写为 $\mathbf{x}_1(t) = \lambda(t) C_1 + \mu(t) C_1 + \nu(t) |C_1|\mathbf{e}_3$ 。注意 $$C_1$$ 和 $$C_1$$ 是不同向量！ $$C_1$$ 是复数， $C_1 = R(\mathbf{e}_1 - i\mathbf{e}_2)/2$ ， $C_1 = R(\mathbf{e}_1 + i\mathbf{e}_2)/2$ ，它们线性无关（在 $\mathbb{C}^3$ ）。但为使 $\mathbf{x}_1$ 实，需系数满足 $\mu = \bar{\lambda}$ ？因为 $\mathbf{x}_1$ 是实向量函数，其傅里叶级数为 $\mathbf{x}_1(s) = \sum_n \mathbf{x}_{1,n} e^{ins/R}$ 。 $$n=1$$ 模是 $\mathbf{x}_{1,1} e^{is/R}$ ， $$n=-1$$ 模是 $\mathbf{x}_{1,-1} e^{-is/R}$ 。根据 (2.10)，圆环的 $\mathbf{x}^R_0$ 有 $\mathbf{x}^R_{0,1} = C_1$ ， $\mathbf{x}^R_{0,-1} = \bar{C}_1$ 。对一般实 $\mathbf{x}_1$ ，有 $\mathbf{x}_{1,-1} = \overline{\mathbf{x}_{1,1}}$ 。现在将 $\mathbf{x}_{1,1}$ 和 $\mathbf{x}_{1,-1}$ 用基表示：
$\mathbf{x}_{1,1} = \lambda C_1 + \mu \bar{C}_1 + \nu_+ |C_1|\mathbf{e}_3?$ 但 $\mathbf{e}_3$ 是实，而 $e^{is/R}$ 的 $$n=1$$ 部分不应有 $\mathbf{e}_3$ 因 $\mathbf{e}_3$ 是常向量，其傅里叶系数只在 $$n=0$$ 。所以对 $n=\pm1$ ， $\mathbf{e}_3$ 分量为零？实际上， $\mathbf{x}_1$ 是实函数，其 $$n=1$$ 和 $$n=-1$$ 傅里叶系数是向量，可以包含 $\mathbf{e}_3$ 分量。例如 $\mathbf{x}_1(s) = \cos(s/R)\mathbf{e}_3$ 的 $n=\pm1$ 系数为零，但 $\sin(s/R)\mathbf{e}_3$ 的 $n=\pm1$ 系数为 $\pm i/2$ 乘 $\mathbf{e}_3$ ？更准确：若 $\mathbf{x}_1(s) = a \cos(s/R) + b \sin(s/R)$ 在 $\mathbf{e}_3$ 方向，则傅里叶系数： $\mathbf{x}_{1,1} = (a - ib)/2$ ， $\mathbf{x}_{1,-1} = (a+ib)/2$ 。所以 $\mathbf{e}_3$ 分量可出现在 $\pm1$ 模。因此 $\nu|C_1|\mathbf{e}_3$ 应分配到 $$n=1$$ 和 $$n=-1$$ 的系数中？在分解 (2.18) 中， $\mathbf{x}_1(t) = \lambda C_1 + \mu \bar{C}_1 + \nu(t) |C_1|\mathbf{e}_3$ ，注意右边是函数， $C_1 e^{is/R}$ 贡献 $$n=1$$ ， $\bar{C}_1 e^{-is/R}$ 贡献 $$n=-1$$ ， $|C_1|\mathbf{e}_3$ 是常数（ $$n=0$$ ）？不， $\mathbf{e}_3$ 乘以常数是 $$n=0$$ 。但这里的 $\nu|C_1|\mathbf{e}_3$ 是常数项（ $$n=0$$ ），它应属于 $$g$$ （零模），不属于 $\boldsymbol{\psi}$ 。但 $\boldsymbol{\psi} = P_{\pm1}\boldsymbol{\phi}$ ，而 $P_{\pm1}\mathbf{x} = P_{\pm1}(\mathbf{x}^R + \boldsymbol{\phi}) = \mathbf{x}^R_{\pm1} + \boldsymbol{\psi}$ 。 $\mathbf{x}^R_{\pm1}$ 是 $C_1 e^{is/R} + \bar{C}_1 e^{-is/R}$ 。如果 $\mathbf{x}_1$ 有常数 $\mathbf{e}_3$ 分量，那属于 $$n=0$$ ，即零模 $$g$$ 的一部分。所以 $\boldsymbol{\psi}$ 应无 $\mathbf{e}_3$ 分量？但 (2.18) 中 $\nu|C_1|\mathbf{e}_3$ 属于 $\mathbf{x}_1$ ，而 $\mathbf{x}_1$ 是 $\mathbf{x}$ 的第一个分量。如果 $\mathbf{x}$ 有 $$x_3$$ 分量，那可能出现在其他分量。论文中 $\boldsymbol{\psi} = P_{\pm1}\boldsymbol{\phi}$ ，故 $\boldsymbol{\psi}$ 本身是一个向量函数，其各个分量都有±1模。但 $\mathbf{e}_3$ 方向可以有±1模吗？可以，如 $\sin(s/R)\mathbf{e}_3$ 有 $\pm1$ 模。但注意 $\mathbf{x}^R_0$ 的 $$x_3$$ 分量为零（圆在 $$xy$$ 平面），故 $\mathbf{x}^R$ 的 $$x_3$$ 随时间线性增长（整体平移），但 $\mathbf{x}^R_0$ 本身 $$x_3=0$$ 。在分解 $\mathbf{x} = \mathbf{x}^R + \boldsymbol{\phi}$ 中， $\mathbf{x}^R$ 的 $\pm1$ 模仅在 $$x_1,x_2$$ ， $$x_3$$ 的 $\pm1$ 模为零。因此 $\boldsymbol{\psi}$ 的 $$x_3$$ 分量可以非零，这来自 $\phi$ 的 $$x_3$$ 的±1模。但在 Lemma 2.6 的分解中，将 $\mathbf{x}_1$ （ $$x_1$$ 的第一个分量）分解，这似乎只涉及 $$x_1$$ 的傅里叶系数，而 $$x_1$$ 是 $\mathbf{x}$ 的第一个分量。 $$x_1$$ 可以有 $n=\pm1$ 和 $$n=0$$ 等。 $\mathbf{e}_3$ 项 $\nu|C_1|\mathbf{e}_3$ 是常向量，所以它只影响 $$x_1$$ 的 $$n=0$$ 系数，即 $$x_1$$ 的零模。但 $\boldsymbol{\psi}$ 只包含 $n=\pm1$ ，所以 $\mathbf{x}_1$ 的 $\pm1$ 部分应只来自 $\lambda C_1 e^{is/R} + \mu \bar{C}_1 e^{-is/R}$ ，而 $\nu|C_1|\mathbf{e}_3$ 属于 $$x_1$$ 的零模。然而，在 (2.18) 中， $\mathbf{x}_1(t)$ 是整个 $$x_1$$ 分量，包括所有频率。所以 $\lambda C_1 + \mu \bar{C}_1$ 给出 $\pm1$ 部分， $\nu|C_1|\mathbf{e}_3$ 是零模部分。因此 $\nu$ 实际是 $$x_1$$ 的零模系数。但 $$x_1$$ 的零模也受弧长条件约束？弧长条件涉及 $|\mathbf{x}_s|^2=1$ ，其中 $\mathbf{x}_s$ 的 $$n=0$$ 为零（因 $\mathbf{x}_s$ 均值为零），所以 $x_{1s}$ 的零模为零，故 $$x_1$$ 的零模是常数，但由 $|\mathbf{x}_s|=1$ 和 $\mathbf{x}^R$ 的贡献，可约束 $$x_1$$ 的零模。实际上，在圆环上 $$x_1$$ 的零模为零（因 $R\cos(s/R)$ 均值为零）。对扰动， $$x_1$$ 的零模是 $$g_1$$ （ $\mathbf{g}$ 的第一个分量）。但 $$g$$ 是 $\boldsymbol{\phi}$ 的零模，而 $\boldsymbol{\phi} = \mathbf{x} - \mathbf{x}^R$ ，故 $\mathbf{x}_1$ 的零模就是 $$g_1$$ 。因此 $\nu$ 本质上与 $$g_1$$ 相关。但 Lemma 2.6 约束 $\nu$ 小，这很自然： $$g_1$$ 是平移的一部分，但平移已被 $\boldsymbol{\tau}$ 消除，然而 $$g$$ 仍有分量，它通过脉冲守恒被约束。因此 Lemma 2.6 实际约束了 $\boldsymbol{\psi}$ 的幅度（通过 $|\lambda|$ ）和 $$g$$ 的 $$x_1$$ 分量（通过 $\nu$ ），但 $\boldsymbol{\psi}$ 和 $$g$$ 是不同频段。然而在平方和 $\|\boldsymbol{\psi}\|^2 \approx |\lambda-1|^2+|\mu|^2$ （不包括 $\nu$ ），因为 $\nu$ 在 $\boldsymbol{\psi}$ 中？不， $\boldsymbol{\psi}$ 不含零模，故 $\|\boldsymbol{\psi}\|^2$ 只涉及 $\lambda,\mu$ （因为 $\psi_1 = (\lambda-1)C_1 + \mu\bar{C}_1$ ，参考圆环的 $$C_1$$ ）。而 $\nu$ 是 $\mathbf{x}_1$ 的零模，即 $$g_1$$ ，属于平移。Lemma 2.7 说 $\|\boldsymbol{\psi}\| \approx |\lambda-1|^2+|\mu|^2+|\nu|^2$ ，这似乎包含 $\nu$ ？但 $\boldsymbol{\psi}$ 应无零模，矛盾。检查 Lemma 2.7 的证明：它说 $\|\boldsymbol{\psi}(t)\| = \sqrt{|\psi_1(t)|^2+|\psi_{-1}(t)|^2} = \sqrt{2}|\psi_1(t)|$ ，且 $\psi_1(t) = \mathbf{x}_1(t) - \mathbf{x}^R_1(t)$ 。而 $\mathbf{x}^R_1(t) = C_1 e^{is/R}$ （但 $\mathbf{x}^R$ 随时间平移，但 $\pm1$ 模不变）。 $\mathbf{x}_1(t)$ 由 (2.18) 给出，其中 $\lambda C_1 e^{is/R} + \mu \bar{C}_1 e^{-is/R}$ 是 $\pm1$ 部分， $\nu|C_1|\mathbf{e}_3$ 是零模。故 $\psi_1(t) = (\lambda-1)C_1 + \mu\bar{C}_1$ （因为参考圆环的 $$n=1$$ 模是 $$C_1$$ ），而 $\nu|C_1|\mathbf{e}_3$ 不属于 $\boldsymbol{\psi}$ 。因此 $\|\boldsymbol{\psi}\|^2 = 2|\psi_1|^2 = 2|(\lambda-1)C_1 + \mu\bar{C}_1|^2$ 。由于 $$C_1$$ 和 $\bar{C}_1$ 正交（ $C_1\cdot\bar{C}_1=0$ ），有 $|(\lambda-1)C_1 + \mu\bar{C}_1|^2 = |\lambda-1|^2|C_1|^2 + |\mu|^2|\bar{C}_1|^2 = |C_1|^2(|\lambda-1|^2+|\mu|^2)$ 。而 $$|C_1|^2 = R^2/2$$ ，故 $\|\boldsymbol{\psi}\|^2 \approx |\lambda-1|^2+|\mu|^2$ 。Lemma 2.7 写为 $\|\boldsymbol{\psi}\| \approx \sqrt{|\lambda-1|^2+|\mu|^2+|\nu|^2}$ ，多了一个 $|\nu|^2$ ，这可能是笔误或我误解。仔细看原文：“Lemma 2.7 The ±1 modes ψ satisfy $\|\psi(t)\| \approx |\lambda(t)-1|^2+|\mu(t)|^2+|\nu(t)|^2$ ”。这似乎错误，因为 $\nu$ 不属于 $\boldsymbol{\psi}$ 。可能 $\nu$ 是出现在 $$x_1$$ 的零模，但 $$x_1$$ 的零模是 $$g_1$$ ，而 $\boldsymbol{\psi}$ 只含±1模。但在距离估计中， $\varepsilon_\perp$ 涉及 $|\mathbf{x}_1\cdot C_1|$ ，其中 $\mathbf{x}_1 = \lambda C_1 + \mu\bar{C}_1 + \nu|C_1|\mathbf{e}_3$ ，且 $C_1\cdot\mathbf{e}_3=0$ ，故 $\mathbf{x}_1\cdot C_1 = \lambda|C_1|^2$ ，因此 $\nu$ 不影响角度误差，只影响幅度误差（通过 $|\mathbf{x}_1|$ ）。而 $\|\boldsymbol{\psi}\|$ 确实不应含 $\nu$ 。我怀疑 Lemma 2.7 的公式有误，或 $\nu$ 是 $\boldsymbol{\psi}$ 在 $$x_3$$ 方向的±1模？但 $\boldsymbol{\psi} = P_{\pm1}\boldsymbol{\phi}$ ，若 $\boldsymbol{\phi}$ 的 $$x_3$$ 有±1模，则 $\boldsymbol{\psi}_3$ 非零。但在分解 (2.18) 中， $\mathbf{x}_1$ 是第一个分量，与 $\boldsymbol{\psi}_3$ 无关。所以 $\nu$ 一定属于 $$x_1$$ 的零模。因此 $\|\boldsymbol{\psi}\|$ 应只依赖 $\lambda,\mu$ 。我保持质疑，但可以肯定：引理 2.6 约束 $|\lambda|-1, |\mu|, |\nu|$ ，其中 $|\nu|$ 小意味着 $$x_1$$ 的零模小，这与平移模 $$g_1$$ 相关。但平移模整体 $$g$$ 包含三个分量，而脉冲守恒约束了整个向量 $P[\mathbf{x}]$ ，可能也约束了 $$g$$ 的其他分量。论文未详细分析 $$g_2,g_3$$ ，但定理 (1.5) 只估计 $P_0(\mathbf{x} - \mathbf{x}^R - \mathbf{x}_{0,\text{trans}})$ ，其中 $\mathbf{x}_{0,\text{trans}}$ 是最优初平移，故 $$g$$ 的偏差被控制。实际上， $$g$$ 的方程 (2.6) 右侧有 $P_0 A\boldsymbol{\psi}$ ，计算得 $P_0 A\boldsymbol{\psi} = \frac{2i}{R^3}(-\bar{C}_1 \times \boldsymbol{\psi}_1 + C_1 \times \boldsymbol{\psi}_{-1})$ （见 (2.27)），这涉及 $\boldsymbol{\psi}$ 的空间分量。由引理 2.6， $\boldsymbol{\psi}$ 有界，故 $\partial_t g$ 有界，所以 $$g$$ 线性增长。但初值脉冲条件也约束了 $$g(0)$$ 的大小？ $$g(0)$$ 是 $\boldsymbol{\phi}_0$ 的零模，而 $\boldsymbol{\phi}_0 = \mathbf{x}_0 - \mathbf{x}^R_0$ 。脉冲 $P[\mathbf{x}_0]$ 接近 $P[\mathbf{x}^R_0]$ ，但 $P[\mathbf{x}^R_0] = \pi R^2\mathbf{e}_3$ ，而 $P[\boldsymbol{\phi}]$ 涉及 $\mathbf{x} \times \mathbf{x}_s$ 的零模。对 $\mathbf{x}^R_0$ ， $\mathbf{x}^R_{0s}$ 的±1模，故 $P[\mathbf{x}^R]$ 完全来自±1模交叉项。 $P[\boldsymbol{\phi}]$ 可能含零模与±1模交叉等。但引理 2.6 的推导已从脉冲守恒推出了关于 $\lambda,\mu,\nu$ 的约束，而 $\lambda,\mu$ 来自 $\mathbf{x}_1$ 的±1模， $\nu$ 来自 $\mathbf{x}_1$ 的零模。因此脉冲守恒实际也约束了 $$g_1$$ （即 $\nu$ ）。而 $$g_2,g_3$$ 呢？未明确，但由对称性可能类似。总之，技术核心是谱分解与守恒律结合。

解决非线性项和低频模态的技巧：

非线性项 $N\boldsymbol{\phi} = \boldsymbol{\phi}_s \times \boldsymbol{\phi}$ 在频域产生三波相互作用。但在各模态方程中，对 $\boldsymbol{\psi}$ ， $P_{\pm1} N(\boldsymbol{\psi}+\tilde{\boldsymbol{\phi}})$ 含 $\boldsymbol{\psi} \times \boldsymbol{\psi}$ 但 $P_{\pm1}(\boldsymbol{\psi}_s \times \boldsymbol{\psi})$ 其中 $\boldsymbol{\psi}$ 仅±1模，其叉积的频率为 $n=\pm2$ （因 $\pm1 \pm 1$ ），故 $P_{\pm1}$ 作用为零。类似 $P_{\pm1}(\boldsymbol{\psi} \times \tilde{\boldsymbol{\phi}})$ 保留，因 $\tilde{\boldsymbol{\phi}}$ 高频。因此 $\boldsymbol{\psi}$ 的方程 (2.5) 为线性驱动： $\partial_t \boldsymbol{\psi} = P_{\pm1} A\tilde{\boldsymbol{\phi}} + P_{\pm1}(\boldsymbol{\psi} \times \tilde{\boldsymbol{\phi}} + \tilde{\boldsymbol{\phi}} \times \boldsymbol{\psi}) + \text{其他非线性}$ ，但 $P_{\pm1}(\tilde{\boldsymbol{\phi}} \times \tilde{\boldsymbol{\phi}})$ 可能非零？ $P_{\pm1}$ 作用在 $\tilde{\boldsymbol{\phi}} \times \tilde{\boldsymbol{\phi}}$ 上，由于 $\tilde{\boldsymbol{\phi}}$ 频率 $|n|\ge2$ ，其叉积频率 $|n|\ge4$ 或 $$0$$ ，故 $P_{\pm1}$ 为零。所以 $\partial_t \boldsymbol{\psi}$ 右边的非线性项只有 $\boldsymbol{\psi} \times \tilde{\boldsymbol{\phi}}$ 型，即一个低频与一个高频的乘积，其频率为 $|n|\ge1$ ，保留±1投影。这些项受 $\|\boldsymbol{\psi}\|_\infty \|\tilde{\boldsymbol{\phi}}\|_2$ 控制，但 $\boldsymbol{\psi}$ 可能 $$O(1)$$ ，而 $\tilde{\boldsymbol{\phi}}$ 小，所以 $\boldsymbol{\psi}$ 方程是准线性驱动，不显式解出，而是通过守恒律和能量方法间接控制。

关键洞察：不直接解 $\boldsymbol{\psi}$ 的时间演化，而是利用守恒量（脉冲）给出代数约束，这些约束在任何时间成立。因此 $\boldsymbol{\psi}$ 不必衰减，只要其参数被约束即可。这避免了处理低频模态的长期动力学（通常有线性增长或慢变）。

零模的处理： $$g$$ 由 (2.6) 驱动，右边含 $P_0 A\boldsymbol{\psi}$ （线性耦合）和 $$P_0 N$$ 。 $P_0 A\boldsymbol{\psi}$ 是 $\boldsymbol{\psi}$ 的线性泛函，由 (2.27) 显式计算，涉及 $\boldsymbol{\psi}_1$ 和 $\boldsymbol{\psi}_{-1}$ 的叉积。由引理 2.7， $\boldsymbol{\psi}$ 有界，故 $P_0 A\boldsymbol{\psi}$ 有界。 $$P_0 N$$ 涉及 $\boldsymbol{\psi}_s \times \boldsymbol{\psi}$ 和 $\tilde{\boldsymbol{\phi}}$ 交叉项。 $\boldsymbol{\psi}_s$ 频率 $$|n|=0,2$$ ？因 $\boldsymbol{\psi}$ 仅±1模， $\boldsymbol{\psi}_s$ 有 $$n=0$$ 和 $n=\pm2$ 模。 $P_0(\boldsymbol{\psi}_s \times \boldsymbol{\psi})$ 可能非零（零模来自 $$n=1$$ 和 $$n=-1$$ 叉积）。但 $\|\boldsymbol{\psi}_s\| \lesssim \|\boldsymbol{\psi}\|_{H^1} \lesssim \|\boldsymbol{\psi}\|_2$ （因低正则），若有界则该项有界。 $\tilde{\boldsymbol{\phi}}$ 相关项受 $\delta$ 控制。因此 $\partial_t g$ 全局有界（因为 $\boldsymbol{\psi}$ 由引理 2.6 和 2.7 视为 $O(\delta)$ 或 $$O(1)$$ ？注意引理 2.6 只给出 $\lambda,\mu,\nu$ 小，但 $\lambda$ 是复数， $|\lambda|$ 接近 1， $\lambda-1$ 小？由 $|\lambda|-1$ 小，则 $|\lambda-1| \le |\lambda|-1 + 2|\lambda-1|?$ 实际上 $||\lambda|-1| \le |\lambda-1|$ ，故 $|\lambda-1| \ge ||\lambda|-1|$ ，但反方向不成立。若 $\lambda = -0.9$ ，则 $|\lambda|=0.9$ ， $||\lambda|-1|=0.1$ 小，但 $|\lambda-1|=1.9$ 大。但 $|\lambda|$ 接近 1 不意味 $\lambda$ 接近 1，可以接近 -1 或其他方向。然而，由 $\lambda$ 是复数， $|\lambda| \approx 1$ ，但 $\arg(\lambda)$ 自由？在 $\mathbf{x}_1$ 的表示中， $\lambda C_1 e^{is/R}$ ， $$C_1$$ 固定。圆环对应 $\lambda=1$ 。若 $\lambda = e^{i\theta}$ 模 1，则 $\lambda C_1$ 是 $$C_1$$ 的旋转，对应圆环的相位偏移。但圆环旋转 $Q_\alpha$ 会改变 $$C_1$$ 为 $e^{i\alpha}C_1$ ，故 $\lambda = e^{i\alpha}$ 的扰动实际上就是圆环的旋转，属于 $\Sigma$ 本身，轨道距离应为零。但在分解中， $\mathbf{x}_1$ 的 $\pm1$ 部分 $\lambda C_1 + \mu\bar{C}_1$ ，若 $\mu=0$ 且 $\lambda=e^{i\alpha}$ ，则 $\mathbf{x}_1 = e^{i\alpha}C_1 e^{is/R} + \text{其他}$ ，这正对应 $Q_\alpha\mathbf{x}^R_0$ 的 $$x_1$$ （旋转后）。但对 $\mathbf{x}$ 的整体旋转不仅影响 $$x_1$$ ，还影响 $$x_2,x_3$$ 。然而在轨道距离中，我们允许旋转整个向量。因此，若扰动仅是整体旋转，则 $\boldsymbol{\psi}$ 应完全可由旋转吸收。在基分解中，若 $\lambda = e^{i\alpha}$ ， $\mu=0$ ， $\nu=0$ ，且其他分量相应旋转，则 $\mathbf{x} \in \Sigma$ 。但 $\lambda$ 为一般复数，若 $|\lambda|=1$ 但 $\arg(\lambda)$ 任意，且 $\mu=0,\nu=0$ ，则 $\mathbf{x}_1 = \lambda C_1 e^{is/R} + \bar{\lambda} \bar{C}_1 e^{-is/R}$ （因 $\mu=\bar{\lambda}$ ），这对应一个圆，但旋转角由 $\arg(\lambda)$ 决定，且 $$x_2$$ 分量自动匹配。因此，若 $\mu=0$ ， $\nu=0$ ， $|\lambda|=1$ ，则扰动是纯旋转，属于 $\Sigma$ 。Lemma 2.6 给出 $|\lambda|-1$ 小，但若 $\lambda$ 接近 -1，则 $|\lambda-1|$ 大，但轨道距离可能仍小，因为旋转 $\pi$ 可对齐。距离公式 (2.12) 为 $\min_\sigma |\mathbf{x}_1 - e^{i\sigma/R}C_1|$ ，而 $\mathbf{x}_1$ 的 $$n=1$$ 模是 $\lambda C_1$ ，故 $\min |\lambda C_1 - e^{i\sigma/R}C_1| = |C_1| \min |\lambda - e^{i\sigma/R}| = |C_1| \big||\lambda|-1\big|$ （因可旋转使相位对齐）。因此，若 $|\lambda|$ 接近 1，则 $\min$ 小，即使 $\lambda$ 为负。所以轨道距离只依赖 $|\lambda|$ 与 1 的差，不依赖 $\arg(\lambda)$ 。因此 $\lambda$ 的相位不敏感，关键是幅度 $|\lambda|$ 接近 1。而 $|\mu|$ 和 $|\nu|$ 小确保无其他形状畸变。故 $\|\boldsymbol{\psi}\| = \sqrt{2}|\psi_1| = \sqrt{2}|(\lambda-1)C_1 + \mu\bar{C}_1|$ 。由于 $$C_1$$ 和 $\bar{C}_1$ 正交， $|\psi_1|^2 = |\lambda-1|^2|C_1|^2 + |\mu|^2|\bar{C}_1|^2$ 。而 $|\lambda-1|$ 可大即使 $|\lambda|$ 接近 1（如 $\lambda=-0.99$ 时 $|\lambda-1|=1.99$ ），但此时 $\mu$ 必须小，因为 $|\mu| = |\lambda|$ （若 $\mu=\bar{\lambda}$ ）？不， $\mu$ 不一定等于 $\bar{\lambda}$ 。从实值条件， $\mu$ 应满足 $\mu = \bar{\lambda}$ ？之前推导 $\mu=\bar{\lambda}$ 从 $x_{1,-1} = \overline{x_{1,1}}$ 和 $C_1,\bar{C}_1$ 的线性无关性。具体：设 $x_{1,1} = \alpha$ ，则 $x_{1,-1} = \bar{\alpha}$ 。又 $x_{1,1} = \lambda C_1 + \mu \bar{C}_1$ ， $x_{1,-1} = \bar{\lambda} C_1 + \bar{\mu} \bar{C}_1$ ？不对，因为 $C_1 e^{-is/R}$ 的系数是 $\mu$ 对应 $$n=-1$$ ？在 (2.18) 中， $\lambda C_1 e^{is/R} + \mu \bar{C}_1 e^{-is/R}$ ，所以 $$n=1$$ 系数是 $\lambda C_1$ ， $$n=-1$$ 系数是 $\mu \bar{C}_1$ 。因此实值条件要求 $\mu \bar{C}_1 = \overline{\lambda C_1} = \bar{\lambda} \bar{C}_1$ ，故 $\mu = \bar{\lambda}$ （因 $\bar{C}_1 \neq 0$ ）。因此 $\mu$ 完全由 $\lambda$ 确定。于是 $\boldsymbol{\psi}$ 实际上只含一个复参数 $\lambda$ （幅度和相位）和一个实参数 $\nu$ （ $$x_1$$ 零模）。但 $\nu$ 不在 $\boldsymbol{\psi}$ 中，而在 $$g_1$$ 中。所以 $\|\boldsymbol{\psi}\|^2 \approx |\lambda-1|^2 + |\bar{\lambda}|^2 = |\lambda-1|^2 + |\lambda|^2$ 。由 $|\lambda| \approx 1$ ，若 $|\lambda|$ 接近 1 但 $\lambda$ 远离 1（如 $\lambda \approx -1$ ），则 $|\lambda-1| \approx 2$ ， $|\lambda| \approx 1$ ，故 $\|\boldsymbol{\psi}\| \approx \sqrt{4+1} = \sqrt{5}$ ，不小！但轨道距离： $\min_\sigma |\lambda C_1 - e^{i\sigma} C_1| = |C_1| \min |\lambda - e^{i\sigma}| = |C_1| \big| |\lambda| - 1 \big|$ （因为可以旋转使 $e^{i\sigma}$ 与 $\lambda$ 同方向）。若 $\lambda = -0.99$ ，则 $|\lambda|=0.99$ ， $||\lambda|-1|=0.01$ ，故轨道距离小。显然 $\|\boldsymbol{\psi}\|$ 大但轨道距离小，因为 $\boldsymbol{\psi}$ 可以近似为旋转（此时 $\lambda$ 模为 1 且 $\mu=\bar{\lambda}$ 给出圆，但 $\lambda \neq 1$ 意味着圆旋转了，但 $\boldsymbol{\psi}$ 本身是 $\lambda C_1 - C_1 + \mu \bar{C}_1$ ，若 $\lambda=e^{i\alpha}$ ， $\mu=\bar{\lambda}=e^{-i\alpha}$ ，则 $\psi_1 = (e^{i\alpha}-1)C_1 + e^{-i\alpha}\bar{C}_1$ 。这通常不为零，除非 $\alpha=0$ 。但这个 $\boldsymbol{\psi}$ 对应一个旋转后的圆吗？设 $\mathbf{x} = Q_\alpha \mathbf{x}^R_0$ ，则其 $$x_1$$ 为 $\cos(s/R + \alpha)$ 乘 $$R$$ ，即 $R[\cos\alpha \cos(s/R) - \sin\alpha \sin(s/R)]$ 。傅里叶系数： $$n=1$$ 模为 $(R/2)(\cos\alpha - i\sin\alpha) = (R/2)e^{-i\alpha}$ ？标准： $\cos(s/R+\alpha) = \cos\alpha \cos(s/R) - \sin\alpha \sin(s/R) = \frac{1}{2}[e^{i\alpha}e^{is/R} + e^{-i\alpha}e^{-is/R}] + \frac{1}{2}[e^{i\alpha}e^{-is/R} + e^{-i\alpha}e^{is/R}]$ ？建议重算： $\cos(\theta+\alpha) = \cos\theta\cos\alpha - \sin\theta\sin\alpha = \frac{e^{i(\theta+\alpha)}+e^{-i(\theta+\alpha)}}{2}$ 不，用欧拉公式： $\cos\phi = (e^{i\phi}+e^{-i\phi})/2$ ，故 $\cos(s/R+\alpha) = \frac{1}{2}e^{i\alpha}e^{is/R} + \frac{1}{2}e^{-i\alpha}e^{-is/R}$ 。所以 $x_1(s) = R\cos(s/R+\alpha) = \frac{R}{2}e^{i\alpha}e^{is/R} + \frac{R}{2}e^{-i\alpha}e^{-is/R}$ 。因此 $$n=1$$ 模为 $\frac{R}{2}e^{i\alpha}$ ， $$n=-1$$ 模为 $\frac{R}{2}e^{-i\alpha}$ 。而 $C_1 = \frac{R}{2}(\mathbf{e}_1 - i\mathbf{e}_2)$ 是向量，但 $$x_1$$ 是标量函数。这里 $$x_1$$ 是第一个分量，对圆 $x_1(s) = R\cos(s/R)$ ，其 $$n=1$$ 模为 $$R/2$$ （实）。在基表示中， $$C_1$$ 是向量，但 $\mathbf{x}_1$ 作为向量函数，其傅里叶系数 $\mathbf{x}_{1,n}$ 是向量。对于 $Q_\alpha \mathbf{x}^R_0$ ， $\mathbf{x}$ 整体旋转，故 $\mathbf{x}_1$ 分量变为 $\cos\alpha\, x_1^{\text{circle}} - \sin\alpha\, x_2^{\text{circle}}$ 。圆环 $x_2^{\text{circle}}(s) = R\sin(s/R)$ 。所以 $[Q_\alpha \mathbf{x}^R_0]_1(s) = R(\cos\alpha \cos(s/R) - \sin\alpha \sin(s/R)) = R\cos(s/R+\alpha)$ 。如上，其 $$n=1$$ 模为 $\frac{R}{2}e^{i\alpha}$ （作为复数，但 $\mathbf{x}_1$ 是实函数，其傅里叶系数是复数， $$n=1$$ 模是 $\frac{R}{2}e^{i\alpha}$ ？注意 $R\cos(s/R+\alpha)$ 的傅里叶级数： $\cos(s/R+\alpha) = \frac{1}{2}(e^{i\alpha} e^{is/R} + e^{-i\alpha} e^{-is/R})$ ，所以 $$n=1$$ 系数是 $\frac{R}{2}e^{i\alpha}$ （复数）， $$n=-1$$ 系数是 $\frac{R}{2}e^{-i\alpha}$ 。但 $\mathbf{x}_1$ 是第一个分量，是标量，其傅里叶系数是复数。在 (2.18) 中， $\mathbf{x}_1(t) = \lambda(t) C_1 + \mu(t) \bar{C}_1 + \nu(t)|C_1|\mathbf{e}_3$ ，但 $C_1, \bar{C}_1$ 是向量，而左边 $\mathbf{x}_1$ 是向量（第一个分量？不， $\mathbf{x}_1$ 可能表示 $\mathbf{x}$ 的第一个分量？符号混淆：通常 $\mathbf{x} = (x_1,x_2,x_3)$ 。在 (2.18) 中， $\mathbf{x}_1$ 可能指 $\mathbf{x}$ 的第一个分量？但 $$C_1$$ 是向量，所以 $\lambda C_1$ 是向量，加上 $\mu\bar{C}_1$ 向量，再加 $\nu|C_1|\mathbf{e}_3$ 向量。因此 $\mathbf{x}_1$ 这里可能不是指第一个分量，而是 $\mathbf{x}$ 的某种表示？看上下文：Lemma 2.6 说“decompose $\mathbf{x}_1$ with respect to the basis $\{C_1, \bar{C}_1, \mathbf{e}_3\}$ ”。而 $\mathbf{x}_1$ 是什么？前文：“We decompose $\mathbf{x}_1(t)$ with respect to the basis $\{C_1, \bar{C}_1, \mathbf{e}_3\}$ of $\mathbb{C}^3$ ”。这里 $\mathbf{x}_1$ 可能表示 $\mathbf{x}$ 的 $$n=1$$ 傅里叶系数？但 $$C_1$$ 是 $\mathbf{x}^R_0$ 的 $$n=1$$ 系数（向量）。所以 $\mathbf{x}_1(t)$ 应指 $\mathbf{x}$ 的 $$n=1$$ 傅里叶系数（向量）。在通常记法中， $\mathbf{x}(s) = \sum_{n\in\mathbb{Z}} \mathbf{x}_n e^{ins/R}$ ，则 $\mathbf{x}_1$ 是 $$n=1$$ 的复系数向量。在 (2.18) 中，他们写 $\mathbf{x}_1(t) = \lambda(t) C_1 + \mu(t) \bar{C}_1 + \nu(t) |C_1|\mathbf{e}_3$ ，这合理： $$C_1$$ 是参考的 $$n=1$$ 系数， $\bar{C}_1$ 是 $$n=-1$$ 系数，但 $\mathbf{x}_1$ 是 $$n=1$$ 系数，为何出现 $\bar{C}_1$ ？因为 $\mathbf{x}_1$ 是 $$n=1$$ 的系数向量，而 $$C_1$$ 和 $\bar{C}_1$ 是 $\mathbb{C}^3$ 中的两个向量，它们线性无关，可以张成 $$n=1$$ 系数的可能空间（因为 $$n=1$$ 系数可以是任意复向量，但受实值条件约束？ $\mathbf{x}$ 是实值曲线，故 $\mathbf{x}_{-n} = \overline{\mathbf{x}_n}$ 。因此 $$n=1$$ 系数 $\mathbf{x}_1$ 是任意复向量，而 $$n=-1$$ 系数由 $\mathbf{x}_{-1} = \overline{\mathbf{x}_1}$ 确定。在基分解中， $\mathbf{x}_1$ 可用 $$C_1$$ 和 $\bar{C}_1$ 展开，因为 $$C_1$$ 和 $\bar{C}_1$ 是 $\mathbb{C}^3$ 的基的一部分？实际上， $$C_1$$ 和 $\bar{C}_1$ 张成 $\mathbb{C}^3$ 的一个二维子空间（ $$xy$$ 平面），加上 $\mathbf{e}_3$ 张成整个 $\mathbb{C}^3$ 。所以任何复向量 $\mathbf{x}_1$ 可唯一分解为 $\lambda C_1 + \mu \bar{C}_1 + \gamma \mathbf{e}_3$ 。而 $\gamma$ 应与 $$x_1$$ 的 $$n=0$$ 系数？不， $\mathbf{x}_1$ 是 $$n=1$$ 系数，它本身是一个向量，其 $\mathbf{e}_3$ 分量就是 $x_{3,1}$ 。所以这里的 $\nu|C_1|\mathbf{e}_3$ 中的 $\nu$ 是标量，表示 $\mathbf{x}_1$ 在 $\mathbf{e}_3$ 方向的分量除以 $$|C_1|$$ 。因此 $\nu$ 对应 $$x_3$$ 的 $$n=1$$ 模。但 $$x_3$$ 的 $$n=1$$ 模是 $\boldsymbol{\phi}$ 的 $$x_3$$ 的±1模的一部分，属于 $\boldsymbol{\psi}$ 的 $$x_3$$ 分量。因此 $\boldsymbol{\psi}$ 确实可能包含 $\mathbf{e}_3$ 方向。而 $$g$$ 是 $\boldsymbol{\phi}$ 的 $$n=0$$ 模，即 $\mathbf{x}$ 的 $$n=0$$ 模减去 $\mathbf{x}^R$ 的 $$n=0$$ 模。 $\mathbf{x}^R$ 的 $$n=0$$ 模： $\mathbf{x}^R(s,t) = \mathbf{x}^R_0(s) + t\mathbf{e}_3$ ， $\mathbf{x}^R_0$ 的 $$n=0$$ 为零（圆均值为零），故 $\mathbf{x}^R$ 的 $$n=0$$ 模为 $t\mathbf{e}_3$ 。所以 $g(t) = P_0(\mathbf{x}(t) - \mathbf{x}^R(t))$ 包含平移和 $$x_3$$ 的零模。因此 $\mathbf{x}_1$ （ $$n=1$$ 系数）中的 $\nu|C_1|\mathbf{e}_3$ 属于 $\boldsymbol{\psi}$ （因是 $$n=1$$ ），不属于 $$g$$ 。所以 $\nu$ 是 $\boldsymbol{\psi}_3$ 的幅度（ $$x_3$$ 的 $$n=1$$ 模）。这样， $\|\boldsymbol{\psi}\|^2$ 应包含 $|\nu|^2$ 的贡献（来自 $$x_3$$ 方向）。Legendre 2.7 可能正确： $\|\boldsymbol{\psi}\| \approx \sqrt{|\lambda-1|^2+|\mu|^2+|\nu|^2}$ 。因为 $\psi_1$ 和 $\psi_{-1}$ 的向量范数平方和包含各分量。设 $\mathbf{x}_1 = (a_1, a_2, a_3)$ 为复向量， $a_3 = \nu|C_1|$ （因 $$C_1$$ 的 $$e_3$$ 分量为零），故 $|a_3| = |\nu||C_1|$ 。 $|C_1| = R/\sqrt{2}$ 。而 $\|\boldsymbol{\psi}\|^2 = 2(|\psi_1|^2 + |\psi_{-1}|^2?$ 实际上 $\|\boldsymbol{\psi}\|_0^2 = \int |\boldsymbol{\psi}|^2 = \sum_{n=\pm1} |\boldsymbol{\psi}_n|^2$ ，且 $\boldsymbol{\psi}_{-1} = \overline{\boldsymbol{\psi}_1}$ ，故 $\|\boldsymbol{\psi}\|_0^2 = 2|\boldsymbol{\psi}_1|^2$ 。而 $\boldsymbol{\psi}_1 = \mathbf{x}_1 - \mathbf{x}^R_{0,1}$ ，其中 $\mathbf{x}^R_{0,1} = C_1$ 。故 $\boldsymbol{\psi}_1 = (\lambda-1)C_1 + \mu\bar{C}_1 + \nu|C_1|\mathbf{e}_3$ 。由于 $C_1, \bar{C}_1, \mathbf{e}_3$ 两两正交（ $C_1\cdot\mathbf{e}_3=0$ , $\bar{C}_1\cdot\mathbf{e}_3=0$ , $C_1\cdot\bar{C}_1=0$ ），有 $|\boldsymbol{\psi}_1|^2 = |\lambda-1|^2|C_1|^2 + |\mu|^2|\bar{C}_1|^2 + |\nu|^2|C_1|^2 = |C_1|^2(|\lambda-1|^2+|\mu|^2+|\nu|^2)$ 。因此 $\|\boldsymbol{\psi}\|^2 = 2|C_1|^2(|\lambda-1|^2+|\mu|^2+|\nu|^2) \asymp |\lambda-1|^2+|\mu|^2+|\nu|^2$ 。故 Lemma 2.7 正确。因此 $\boldsymbol{\psi}$ 包含 $\nu$ ，即 $$x_3$$ 方向的±1模。但 $\mathbf{x}^R_0$ 的 $$x_3=0$$ ，所以 $\boldsymbol{\psi}$ 的 $$x_3$$ 的±1模正比于 $\nu$ 。这意味着涡丝可能扭曲出平面（非平面扰动），但 $\nu$ 被脉冲守恒约束为小。

因此，几何引理约束了 $\lambda,\mu,\nu$ ，从而控制 $\boldsymbol{\psi}$ 的整体范数。结合幅度误差和角度误差，得到轨道距离小ness。

与之前工作的详细比较：

[Aiki, 2025]：该工作研究相同问题但假设扰动对称，具体为 $\boldsymbol{\phi}$ 满足 $P_{\pm1}\boldsymbol{\phi} \cdot \mathbf{e}_3 = 0$ （即±1模无 $$x_3$$ 分量）且可能其他对称条件。这导致 $\nu=0$ ，且由脉冲和弧长条件可直接推出 $\lambda=1, \mu=0$ ，故 $\boldsymbol{\psi}=0$ ，剩下 $$g$$ 的线性增长。稳定性证明较简单，因低频部分被对称性强制为零。本文去除此假设后， $\nu$ 自由，脉冲和弧长条件只能推出 $\lambda,\mu,\nu$ 小，但不为零。需要更精细的几何引理和距离最小化技巧。技术难度增加，且结果更一般。
Calini-Keith-Lafortune [3]：研究各种封闭涡丝（包括圆环）的线性稳定性，使用 Hasimoto 映射（将 LIE 映射到非线性薛定谔方程）。他们考虑“closure-preserving periodic perturbations”，即保持封闭性的周期扰动，并分析线性化算子的谱。本文研究非线性轨道稳定性，允许更一般的非周期扰动？但本文扰动仍保持封闭（因初值 $|\mathbf{x}_{0s}|=1$ 且 LIE 保持弧长，曲线总封闭）。但 [3] 是线性理论，且关注模态增长/衰减；本文是非线性全局估计。此外，[3] 的 perturbation 可能隐含某种对称性？未明确。本质区别：线性稳定性分析特征值，而本文用守恒律直接证明轨道稳定性，不涉及谱。
Choi [4] 及 Choi-Jeong [5]：研究 Hill 涡（三维有界涡旋）在轴对称 Euler 方程下的轨道稳定性。他们使用脉冲（或角动量）守恒控制平移模式。概念上类似：利用几何守恒量（如流体脉冲、角动量）控制低频漂移。但方程不同（Euler vs LIE），Hill 涡是三维涡旋，而涡丝是一维曲线。技术方法也不同：Euler 是体积上的 PDE，LIE 是曲线上的几何进化方程。本文的谱分解和相对能量方法特定于 LIE 的 Hasimoto 变换或弧长保持结构。但“利用脉冲控制低频”的思想是相通的。
Langer-Perline [7]：在弹性曲线框架下，LIE 可视为弯曲能临界点的刚体运动。圆环是弹性曲线的最小值。他们的分析基于变分原理，不显式使用向量脉冲。本文则显式利用 $P[\mathbf{x}]$ 守恒证明圆环是约束局部极小（通过引理 2.6），这提供了不同于弹性曲线观点的稳定性论证，并直接适用于 LIE。

改进的总结：

假设放宽：从对称扰动到任意非对称扰动，这是主要改进。代价是需额外假设高频模 $$H^2$$ 小（这在对称情形可能自动满足，因对称性限制可能迫使高频为零？不一定，但对称情形下高频估计可能仍需要，但文献未提）。实际上 [Aiki,2025] 可能也有类似高频假设，但未在摘要强调。关键是对低频模，对称情形无需引理 2.6 因 $\boldsymbol{\psi}=0$ 自动，而本文必须处理 $\boldsymbol{\psi} \neq 0$ 但受脉冲约束的情况。
工具新颖：几何稳定性引理（Lemma 2.6）是本文原创，通过将 $\mathbf{x}_1$ 在特定基下分解并利用脉冲守恒导出三个参数的小ness。优于仅用弧长条件（只给出幅度和）或仅用脉冲（给出 $|\lambda|^2-|\mu|^2$ ）单独使用。
结果强度：定理 1.1 给出全局 $O(\delta)$ 轨道稳定性，与对称情形一致。平移估计 (1.5) 也类似。因此，在不牺牲结论强弱的前提下，扩展了扰动类，增强了定理的普适性。
函数空间：与之前工作类似，需要 $$H^4$$ 保证唯一性。未降低正则性要求。

技术挑战：非对称性导致±1模可能通过旋转和缩放扰动，且可能有 $$x_3$$ 方向的±1模（非平面扭曲）。如何证明这些低频模仍被控制？脉冲守恒提供向量方程，但需巧妙投影到参数。论文的基选择（ $C_1, \bar{C}_1, \mathbf{e}_3$ ）恰好匹配圆环的傅里叶结构和直角关系，使得脉冲的 $\mathbf{e}_3$ 分量直接给出 $|\lambda|^2-|\mu|^2$ ，而 $\mathbf{e}_3$ 分量自动为零？计算脉冲时， $P[\mathbf{x}]$ 的 $\mathbf{e}_3$ 分量来自 $n=\pm1$ 交叉项，如前计算得 $\pi R^2(|\lambda|^2 - |\mu|^2)$ 加上其他贡献？实际上，如果 $\nu \neq 0$ ， $P[\mathbf{x}]$ 的 $\mathbf{e}_3$ 分量还可能有来自 $\nu$ 与零模的交叉？但零模是 $$n=0$$ ，叉积 $$n=1$$ 与 $$n=-1$$ 给 $$n=0$$ ，而 $$n=0$$ 与 $$n=0$$ 叉积为零。因此 $P[\mathbf{x}]$ 的 $$n=0$$ 部分（即向量脉冲）完全来自 $n=\pm1$ 对。而 $$n=1$$ 和 $$n=-1$$ 的系数分别是 $\mathbf{x}_1, \mathbf{x}_{-1} = \overline{\mathbf{x}_1}$ 。所以 $P[\mathbf{x}] = -i\pi \sum_{|n|=1} n \mathbf{x}_n \times \mathbf{x}_{-n}$ 。计算 $\mathbf{x}_1 \times \mathbf{x}_{-1} = \mathbf{x}_1 \times \overline{\mathbf{x}_1}$ 。用基展开： $\mathbf{x}_1 = \lambda C_1 + \mu \bar{C}_1 + \nu |C_1|\mathbf{e}_3$ 。注意 $$C_1$$ 和 $\bar{C}_1$ 是复向量，但 $\mathbf{x}_1 \times \overline{\mathbf{x}_1}$ 是实向量（因 $u\times \bar{u} = -2i (\Re u \times \Im u)$ ）。展开后，交叉项： $C_1 \times \bar{C}_1 = i|C_1|^2 \mathbf{e}_3$ ， $\bar{C}_1 \times C_1 = -i|C_1|^2 \mathbf{e}_3$ ， $C_1 \times \mathbf{e}_3$ 和 $\mathbf{e}_3 \times \bar{C}_1$ 等是 imag？实际计算： $\mathbf{x}_1 \times \overline{\mathbf{x}_1} = (\lambda C_1 + \mu \bar{C}_1 + \gamma \mathbf{e}_3) \times (\bar{\lambda} \bar{C}_1 + \bar{\mu} C_1 + \gamma \mathbf{e}_3)$ ，其中 $\gamma = \nu|C_1|$ 。由于 $\mathbf{e}_3 \times \mathbf{e}_3=0$ ，且 $C_1 \times \mathbf{e}_3$ 和 $\mathbf{e}_3 \times C_1$ 等。但 $$C_1$$ 和 $\bar{C}_1$ 在 $$xy$$ 平面，叉积 $\mathbf{e}_3$ 。所以 $\mathbf{x}_1 \times \overline{\mathbf{x}_1}$ 是 $\mathbf{e}_3$ 方向的纯向量？实际上， $C_1 \times \bar{C}_1 \propto \mathbf{e}_3$ ， $C_1 \times \mathbf{e}_3 \propto$ 在 $\mathrm{span}\{\mathbf{e}_1,\mathbf{e}_2\}$ ？ $$C_1$$ 在 $$xy$$ 平面，叉乘 $\mathbf{e}_3$ 仍在 $$xy$$ 平面。因此 $\mathbf{x}_1 \times \overline{\mathbf{x}_1}$ 一般有 $$xy$$ 和 $$z$$ 分量。但脉冲守恒是向量，所有分量都约束。然而，在 Lemma 2.6 的证明中，他们只用了 $$e_3$$ 分量（见 (2.23)-(2.24)），写了“the summation term simplifies to”并得到 $\mathbf{e}_3$ 分量。可能其他分量自动为零或由其他条件控制？实际上，由圆环的对称性，对圆环 $P[\mathbf{x}^R] = \pi R^2 \mathbf{e}_3$ ，只有 $\mathbf{e}_3$ 分量。对扰动， $P[\mathbf{x}]$ 的其他分量可能非零，但初值条件 $|P[\mathbf{x}_0] - \pi R^2\mathbf{e}_3| \le \delta^2$ 控制整个向量差，故所有分量都小。但 Lemma 2.6 只用了 $\mathbf{e}_3$ 分量，是否足够？因为 $\mathbf{x}_1 \times \overline{\mathbf{x}_1}$ 的 $$xy$$ 分量涉及 $\nu$ 与 $\lambda,\mu$ 的交叉。例如 $(\lambda C_1) \times (\bar{\gamma} \mathbf{e}_3)$ 等。可能这些交叉项在脉冲中贡献为零或因对称性？检查： $P[\mathbf{x}] = -i\pi \sum_{|n|=1} n \mathbf{x}_n \times \mathbf{x}_{-n}$ 。 $\mathbf{x}_n \times \mathbf{x}_{-n}$ 是纯虚？但最终 $P[\mathbf{x}]$ 实。其 $$e_1,e_2$$ 分量对应 $$xy$$ 平面。但圆环脉冲只有 $$e_3$$ 分量，故扰动脉冲的 $$e_1,e_2$$ 分量应小。这些分量可能给出 $\nu$ 或其他参数的约束。但论文 Lemma 2.6 仅用 $$e_3$$ 分量推出 $|\lambda|^2-|\mu|^2 \approx 1$ 。结合弧长条件 $|\lambda|^2+|\mu|^2+|\nu|^2 \approx 1$ ，可推出 $|\mu|^2+|\nu|^2 \approx (1-|\lambda|^2)/?$ 从两者： (1) $|\lambda|^2+|\mu|^2+|\nu|^2 \approx 1$ (2) $|\lambda|^2-|\mu|^2 \approx 1$ 相加得 $2|\lambda|^2 + |\nu|^2 \approx 2 \Rightarrow |\lambda|^2 \approx 1 - |\nu|^2/2$ ，相减得 $2|\mu|^2 + |\nu|^2 \approx 0$ ，故 $|\mu|^2 \approx -|\nu|^2/2$ ，矛盾除非 $|\nu|^2$ 很小且 $|\mu|^2$ 很小。更准确：由 (2) $|\lambda|^2 - |\mu|^2 = 1 + O(\delta^2)$ ，由 (1) $|\lambda|^2+|\mu|^2+|\nu|^2 = 1 + O(\delta^2)$ 。两式相加： $2|\lambda|^2 + |\nu|^2 = 2 + O(\delta^2)$ ，故 $|\lambda|^2 = 1 - |\nu|^2/2 + O(\delta^2)$ 。相减： $2|\mu|^2 = -|\nu|^2 + O(\delta^2)$ ，因此 $|\mu|^2 \approx -|\nu|^2/2$ ，由于模平方非负，必须有 $|\nu|^2$ 极小（ $O(\delta^2)$ ），且 $|\mu|^2$ 也是 $O(\delta^2)$ 。而 $|\lambda|^2 \approx 1 + O(\delta^2)$ 。再由 $||\lambda|-1| \le |\lambda^2-1|/(|\lambda|+1) \lesssim |\lambda|^2-1$ ，故 $|\lambda|$ 接近 1。因此，仅用 $$e_3$$ 分量和弧长条件，已经推出 $|\mu|,|\nu| \lesssim \delta$ 和 $|\lambda|-1 \lesssim \delta$ （因为 $|\lambda|^2-1 = (|\lambda|^2-|\mu|^2-|\nu|^2) + (|\mu|^2+|\nu|^2) \approx 0 + O(\delta^2)$ ，故 $|\lambda|^2 \approx 1$ ，但 $|\lambda|$ 接近 1 不意味着 $\lambda$ 接近 1，然而轨道距离只关心 $|\lambda|$ ，而 $\|\boldsymbol{\psi}\|$ 中的 $|\lambda-1|$ 可大如果 $\lambda$ 相位偏离 1。但从 $\|\boldsymbol{\psi}\|$ 估计看， $|\lambda-1|$ 不一定小。例如 $\lambda = -1$ ， $|\lambda|=1$ ， $|\lambda-1|=2$ 。但 $\lambda=-1$ 是否满足条件？从 $|\lambda|^2 \approx 1$ 和 $|\mu|^2 \approx -|\nu|^2/2$ ，若 $|\nu|$ 很小，则 $|\mu|$ 必须更小。 $\lambda=-1$ 给出 $|\lambda|^2=1$ ，可满足。但 $\mu$ 必须满足 $|\mu|^2 = (|\lambda|^2-1+|\nu|^2)/?$ 由 $|\lambda|^2-|\mu|^2 \approx 1$ ，若 $|\lambda|^2=1$ ，则 $|\mu|^2 \approx 0$ 。所以 $\mu \approx 0$ 。由 $\mu=\bar{\lambda}$ ？若 $\lambda=-1$ 实，则 $\mu=\bar{\lambda}=-1$ ，但 $|\mu|^2=1$ 不接近 0，矛盾。因此 $\lambda=-1$ 被 $\mu=\bar{\lambda}$ 排除？实值条件强制 $\mu = \bar{\lambda}$ 。如果 $\lambda=-1$ ，则 $\mu = -1$ ， $|\mu|=1$ ，但由脉冲 $|\lambda|^2-|\mu|^2=0 \neq 1$ ，故不满足圆环脉冲条件。所以 $\lambda=-1$ 不允许。实际上， $|\lambda|^2-|\mu|^2=1$ 且 $|\mu|=|\lambda|$ （从 $\mu=\bar{\lambda}$ ）推出 $$0=1$$ ，矛盾除非 $|\lambda|$ 无穷？因此必须有 $\mu \neq \bar{\lambda}$ ？但实值条件要求 $\mu = \bar{\lambda}$ 以确保 $x_{1,-1} = \overline{x_{1,1}}$ 。我们重新审视实值条件： $\mathbf{x}$ 是实值曲线，其傅里叶系数满足 $\mathbf{x}_{-n} = \overline{\mathbf{x}_n}$ 。对于 $$n=1$$ ， $\mathbf{x}_{-1} = \overline{\mathbf{x}_1}$ 。在分解 (2.18) 中， $\mathbf{x}_1$ 是 $\mathbf{x}$ 的 $$n=1$$ 系数。 $\mathbf{x}_{-1}$ 是什么？它应由 $\mathbf{x}_{-1} = \lambda \bar{C}_1 + \mu C_1 + \nu |C_1|\mathbf{e}_3$ ？因为 $\mathbf{x}_{-1}$ 是 $$n=-1$$ 系数，对应 $e^{-is/R}$ ，而 $\bar{C}_1 e^{-is/R}$ 是 $$n=-1$$ ， $C_1 e^{is/R}$ 是 $$n=1$$ 。因此一般形式为 $\mathbf{x}_{-1} = a \bar{C}_1 + b C_1 + c |C_1|\mathbf{e}_3$ 。实值条件 $\mathbf{x}_{-1} = \overline{\mathbf{x}_1} = \bar{\lambda} \bar{C}_1 + \bar{\mu} C_1 + \bar{\nu} |C_1|\mathbf{e}_3$ （因 $\mathbf{e}_3$ 实）。故比较得 $a = \bar{\lambda}$ ， $b = \bar{\mu}$ ， $c = \bar{\nu}$ 。因此 $\mathbf{x}_{-1} = \bar{\lambda} \bar{C}_1 + \bar{\mu} C_1 + \bar{\nu} |C_1|\mathbf{e}_3$ 。这并不强制 $\mu = \bar{\lambda}$ ，而是 $\mathbf{x}_1$ 和 $\mathbf{x}_{-1}$ 有独立的系数 $\lambda,\mu,\nu$ 和 $\bar{\lambda},\bar{\mu},\bar{\nu}$ 。但 $\mathbf{x}_1$ 和 $\mathbf{x}_{-1}$ 通过共轭关联，这给出 $\lambda,\mu,\nu$ 是任意复数？不， $\mathbf{x}_{-1}$ 由 $\mathbf{x}_1$ 决定，所以 $\mathbf{x}_1$ 的系数 $\lambda,\mu,\nu$ 可以是任意复数，而 $\mathbf{x}_{-1}$ 的系数自动为 $\bar{\lambda},\bar{\mu},\bar{\nu}$ 。因此无额外约束 $\mu=\bar{\lambda}$ 。之前我错误地认为 $$C_1$$ 和 $\bar{C}_1$ 是基，但 $\mathbf{x}_1$ 展开中系数 $\lambda,\mu$ 独立，只要 $\mathbf{x}_{-1}$ 相应匹配即可。例如，设 $\mathbf{x}_1 = \lambda C_1 + \mu \bar{C}_1 + \gamma \mathbf{e}_3$ ，则 $\mathbf{x}_{-1} = \bar{\lambda} \bar{C}_1 + \bar{\mu} C_1 + \bar{\gamma} \mathbf{e}_3$ 。这总是实值，因为 $\mathbf{x}_1$ 和 $\mathbf{x}_{-1}$ 共轭对称。因此 $\lambda,\mu,\gamma$ 是独立复参数？但 $\gamma$ 是复？ $\mathbf{e}_3$ 是实向量，所以 $\gamma$ 应为实以确保 $\gamma \mathbf{e}_3$ 的共轭是 $\bar{\gamma} \mathbf{e}_3$ ，而 $\mathbf{x}_{-1}$ 的 $\mathbf{e}_3$ 分量是 $\bar{\gamma} \mathbf{e}_3$ ，所以 $\gamma$ 可以是复吗？如果 $\gamma$ 是复， $\gamma \mathbf{e}_3$ 是复向量，但 $\mathbf{x}$ 是实值曲线，其傅里叶系数 $\mathbf{x}_1$ 可以是复向量，不要求每个分量实。实值条件只要求 $\mathbf{x}_{-1} = \overline{\mathbf{x}_1}$ ，这自动满足如果 $\mathbf{x}_1$ 的系数任意且 $\mathbf{x}_{-1}$ 取共轭系数。因此 $\lambda,\mu,\gamma$ 可以是任意复数。但 $C_1, \bar{C}_1, \mathbf{e}_3$ 是 $\mathbb{C}^3$ 的基（在复意义下），所以任何复向量 $\mathbf{x}_1$ 可唯一分解，系数无约束。因此 $\lambda,\mu,\gamma$ 独立。那么 $|\lambda|$ 和 $|\mu|$ 独立，且 $\gamma$ （即 $\nu|C_1|$ ）的幅度独立。这样，脉冲条件 $|\lambda|^2-|\mu|^2 \approx 1$ 和弧长条件 $|\lambda|^2+|\mu|^2+|\nu|^2 \approx 1$ 确实推出 $|\mu|^2+|\nu|^2 \approx 0$ ，故 $|\mu|,|\nu| \lesssim \delta$ ，而 $|\lambda|^2 \approx 1$ ，因此 $|\lambda| \approx 1$ 。但 $|\lambda-1|$ 仍可能大如果 $\lambda$ 相位远离 0。例如 $\lambda = e^{i\pi} = -1$ ， $|\lambda|=1$ ，满足 $|\lambda| \approx 1$ ，但 $|\lambda-1|=2$ 。这样的 $\lambda$ 是否与 $|\mu|$ 小兼容？若 $\lambda=-1$ ， $|\lambda|^2=1$ ，则从 $|\lambda|^2-|\mu|^2 \approx 1$ 得 $|\mu|^2 \approx 0$ ，从 $|\lambda|^2+|\mu|^2+|\nu|^2 \approx 1$ 得 $|\nu|^2 \approx 0$ 。所以 $\lambda=-1$ 是允许的（如果 $\mu,\nu$ 很小）。此时 $\boldsymbol{\psi}_1 = (-1-1)C_1 + \mu\bar{C}_1 + \nu|C_1|\mathbf{e}_3 \approx -2C_1$ ，所以 $\boldsymbol{\psi}$ 范数 $$O(1)$$ 。但轨道距离： $\min_\alpha \| \mathbf{x}^R_0+\boldsymbol{\psi} - Q_\alpha\mathbf{x}^R_0 \|_{H^2}$ 。计算 $\mathbf{x}^R_0+\boldsymbol{\psi}$ 的 $$n=1$$ 模为 $C_1 + \psi_1 \approx C_1 -2C_1 = -C_1$ 。 $Q_\alpha\mathbf{x}^R_0$ 的 $$n=1$$ 模为 $e^{i\alpha}C_1$ 。取 $\alpha=\pi$ ，则 $e^{i\pi}C_1 = -C_1$ ，完美匹配 $$n=1$$ 模。但还需匹配 $$n=-1$$ 模： $\mathbf{x}^R_0+\boldsymbol{\psi}$ 的 $$n=-1$$ 模是 $\bar{C}_1 + \psi_{-1}$ ，而 $\psi_{-1} = \overline{\psi_1} \approx -2\bar{C}_1$ ，故 $$n=-1$$ 模 $\approx \bar{C}_1 -2\bar{C}_1 = -\bar{C}_1$ 。 $Q_\pi\mathbf{x}^R_0$ 的 $$n=-1$$ 模为 $e^{-i\pi}\bar{C}_1 = -\bar{C}_1$ ，匹配。且更高频模为零。因此轨道距离为零。所以尽管 $\boldsymbol{\psi}$ 大，但它是精确旋转，故在 $\Sigma$ 中。因此，轨道距离仅取决于 $|\lambda|$ （通过 $\varepsilon_\parallel$ ）和 $\mu,\nu$ （通过 $\varepsilon_\perp$ ），而 $|\lambda| \approx 1$ 是足够的，不要求 $\lambda$ 接近 1。因此 $\|\boldsymbol{\psi}\|$ 可能大，但轨道距离小。论文中 Lemma 2.7 给出 $\|\boldsymbol{\psi}\| \approx \sqrt{|\lambda-1|^2+|\mu|^2+|\nu|^2}$ ，当 $\lambda=-1$ 时 $|\lambda-1|^2 \approx 4$ ，但 $\mu,\nu$ 小，故 $\|\boldsymbol{\psi}\|$ 可能 $$O(1)$$ 。然而，轨道距离 $\operatorname{dist}_{\pm1}$ 近似 $\sqrt{(|\lambda|-1)^2 + \text{涉及 }|\mu|,|\nu| \text{ 的项}}$ ，当 $|\lambda|=1$ 时第一项为零。所以轨道距离 $O(\delta)$ 而 $\|\boldsymbol{\psi}\|$ 可能 $$O(1)$$ ，这没问题，因为 $\Sigma$ 包含所有旋转， $\boldsymbol{\psi}$ 可能只是旋转偏离。但注意 $\boldsymbol{\psi}$ 包含 $$x_3$$ 的±1模（ $\nu$ ），若 $\nu \neq 0$ ，则 $\boldsymbol{\psi}$ 不在 $$xy$$ 平面，而 $Q_\alpha\mathbf{x}^R_0$ 在 $$xy$$ 平面，故 $\nu$ 必须小才能匹配。所以 $\nu$ 小确保无显著非平面扭曲。

因此，技术核心是：虽然 $\boldsymbol{\psi}$ 范数可能不衰减，但其“形状变化”被几何引理约束为小，通过允许旋转对齐。这体现了轨道稳定性的本质：不要求解接近一个固定圆，而要求接近某旋转后的圆。

最后，比较中应强调本文的守恒律使用：脉冲守恒给出向量方程，投影后结合弧长条件，解出三个参数的小ness。这比 [Aiki,2025] 中对称性自动迫使 $\boldsymbol{\psi}=0$ 更精细。技巧上，基的选择至关重要，使得 $$C_1$$ 和 $\bar{C}_1$ 正交， $C_1 \times \bar{C}_1 \propto \mathbf{e}_3$ ，从而脉冲的 $\mathbf{e}_3$ 分量分离出 $|\lambda|^2-|\mu|^2$ 。而 $P[\mathbf{x}]$ 的其他分量（如 $$e_1,e_2$$ ）可能给出涉及 $\nu$ 的约束，但论文未显式使用，可能因由 $\mathbf{e}_3$ 分量和弧长条件已足够约束 $\nu$ （如前推导 $2|\mu|^2 + |\nu|^2 \approx 0$ ）。但 $|\mu|^2$ 非负，故 $|\nu|^2 \le 0$ ？不，从 $|\lambda|^2+|\mu|^2+|\nu|^2 = 1+O(\delta^2)$ 和 $|\lambda|^2-|\mu|^2 = 1+O(\delta^2)$ ，相减得 $2|\mu|^2 + |\nu|^2 = O(\delta^2)$ ，故两者均为 $O(\delta^2)$ 。因此 $\nu$ 自动小。所以 $\mathbf{e}_3$ 分量已足够。因此几何引理高效。

总结比较：本文在更一般的扰动类下，通过新引理，实现了与对称情形相同量级的稳定性结论。技术贡献在于利用流体脉冲守恒对低频模态进行代数控制，结合相对能量对高频模态的谱控制，完成整体估计。

2. 毛细波与一类拟线性偏微分方程的不变流形

Invariant Manifolds for Capillary Waves and a Class of Quasilinear PDEs

作者: Jalal Shatah, Chongchun Zeng

类别: math.AP

发布日期: 2026-02-19

arXiv ID: 2602.17844v1

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

详细分析

中文标题: 毛细波与一类拟线性偏微分方程的不变流形

1. 研究对象和背景

论文研究一类拟线性偏微分方程（PDE）的局部稳定流形与局部不稳定流形，其核心动机来源于具有表面张力的无旋水波系统（irrotational water waves with surface tension）。此类系统描述两种不可压缩、无粘、无旋流体在重力与表面张力作用下的自由界面运动，是流体力学与数学物理中的经典模型。论文旨在建立该类系统在平衡点附近动力学的核心结构——不变流形——的存在性、唯一性与光滑性，并发展适用于更广泛拟线性PDE的一般框架。

物理系统与数学表述

考虑两层流体：上层深度为 $$d_1$$ ，下层深度为 $$d_2$$ ，两层之间的自由界面位于 $$z=0$$ ，底部分别固定于 $$z=-d_1$$ 和 $$z=-d_2$$ 。假设流体无旋（速度场可表示为速度势 $\phi_j$ 的梯度， $$j=1,2$$ ）且不可压缩（ $\nabla\cdot\mathbf{u}=0$ ），则每个区域内的速度势满足 Laplace 方程：

\Delta\phi_j = 0 \quad \text{在各自流体区域内}.

边界条件包括：

底部固边界： $\partial_z\phi_j = 0$ 在 $$z=-d_j$$ ；
自由界面运动学条件：界面粒子随流体运动，即 $\eta_t = \mathbf{u}\cdot\mathbf{n}$ ，其中 $\eta(x,t)$ 是界面高度， $\mathbf{n}$ 是界面法向量。对于无旋流，可化为 $\eta_t = \sqrt{1+|\nabla\eta|^2}\,\partial_n\phi_1$ （上层）；
自由界面动力学条件：由 Bernoulli 方程与压力连续条件给出。设 $$p_j$$ 为压力， $\rho_j$ 为密度，则 $\phi_{j,t} + \frac{1}{2}|\nabla\phi_j|^2 + g\eta + \frac{p_j}{\rho_j} = \text{常数},$ 其中 $$g$$ 为重力加速度。表面张力在界面处产生额外压力 $\sigma\kappa$ ， $\sigma$ 为表面张力系数， $\kappa$ 为界面曲率。压力连续条件结合不可压缩性导出界面的演化方程。

通过变量简化，常取上层速度势在界面的迹 $\psi = \phi_1|_{z=\eta}$ ，并引入 Dirichlet-to-Neumann 算子 $G(\eta)$ 使得 $\partial_z\phi_1(\eta) = G(\eta)\psi$ 。这样，系统可约化为关于状态变量 $u = (\eta, \psi)$ 的拟线性 PDE 系统。对于单相水波（一层流体，底部在 $$z=-d$$ ），典型形式为（参见 [54]）：

\begin{aligned} \eta_t &= \sqrt{1+|\nabla\eta|^2}\,G(\eta)\psi - \nabla\eta\cdot\nabla\psi, \\ \psi_t &= -\nabla\psi\cdot\nabla\psi - g\eta + \sigma\kappa - \frac{1}{2}|G(\eta)\psi|^2 - \nabla\psi\cdot(G(\eta)\psi)\,\nabla\eta. \end{aligned} \tag{1}

曲率 $\kappa = \nabla\cdot\left( \nabla\eta / \sqrt{1+|\nabla\eta|^2} \right)$ 。算子 $G(\eta): H^s(\mathbb{R}^d) \to H^{s-1}(\mathbb{R}^d)$ 是依赖于 $\eta$ 的非局部拟微分算子，定义为：

G(\eta)\psi = \sqrt{1+|\nabla\eta|^2}\,\partial_z\phi|_{\text{surface}}, \quad \text{其中 } \Delta\phi=0 \text{ 在 } -d<z<\eta, \quad \phi_z=0 \text{ 在 } z=-d, \quad \phi=\psi \text{ 在 } z=\eta.

对于常数深度 $$d$$ ， $G(\eta)$ 的 Fourier 符号可显式表达，涉及 $\tanh$ 函数和 $\eta$ 的 Fourier 变换。这种非局部性和对 $\eta$ 的依赖使得系统 (1) 是拟线性的：线性部分 $G(\eta)\psi$ 本身依赖于状态 $\eta$ 。

物理平衡点对应静水平面： $\eta=0$ ， $\psi=0$ 。线性化系统给出色散关系，其特征值可能分布在复平面不同区域：重力波通常对应实谱（传播波），而表面张力主导的短波可能对应纯虚谱（稳定振荡）。通过选择 Sobolev 空间 $$H^s$$ 的指数 $$s$$ 足够大，线性化算子 $$A = DF(0)$$ 在相空间 $X = H^s \times H^{s-1/2}$ （或类似）上可分解为不变子空间 $X_+ \oplus X_-$ ，其中 $$X_+$$ 对应指数增长方向（不稳定）， $$X_-$$ 对应指数衰减方向（稳定）。这种分解依赖于维度和参数（如 $\sigma/g$ 的比率）。

数学上，论文考虑抽象形式 $$u_t = F(u)$$ 在实 Hilbert 空间 $$X$$ 中， $$F(0)=0$$ 。假设 $$A = DF(0)$$ 有不变分解 $X = X_+ \oplus X_-$ ，使得 $e^{tA}$ 在 $$X_+$$ 上指数增长，在 $$X_-$$ 上指数衰减（即存在 $\lambda_+ > \lambda_-$ 满足 (1.2)）。非线性 $$F$$ 被假定在某个 Sobolev 尺度 $X^n = \text{Dom}(A^n)$ 上具有足够光滑性（ $n \ge 2$ ），但可能允许正则性损失——这是拟线性系统的核心特征。

挑战与动机

不变流形是动力系统局部结构的基石：它们将相空间局部分解为沿流形切于不变子空间 $X_\pm$ 的曲线，并 organizing 附近轨道的渐近行为。对于 ODE，当 $$A$$ 有指数二分性时，局部稳定/不稳定流形存在唯一（经典定理，如 Hadamard 图变换或 Lyapunov-Perron 方法）。对于无限维系统：

半线性 PDE $$u_t = Au + f(u)$$ ，其中 $$A$$ 生成光滑半群（如抛物型或某些双曲型），不变流形理论已成熟，因为 $e^{tA}$ 提供正则性增益，允许在 $$X^n$$ 中迭代。
抛物型拟线性 PDE（如 Navier-Stokes）利用极大正则性， $e^{tA(u)}$ 也光滑化，因此流形可构造。然而，对于拟线性 Hamiltonian PDE（如水波、非线性 Schrödinger 方程、KdV 等），存在根本障碍：
非线性导致正则性损失：拟线性项 $$A(u)u$$ 包含 $$u$$ 的空间导数，而 $$A(u)$$ 本身依赖于 $$u$$ 。在 Lyapunov-Perron 方法中，需反复估计积分 $\int e^{(t-\tau)A}f(v(\tau))d\tau$ 。若 $e^{tA}$ 不光滑化且 $$f(v)$$ 含导数，则每次迭代可能降低正则性，且无法恢复，导致迭代序列不收敛于 $$X^n$$ 。
缺乏光滑化效应：对于水波，线性色散算子 $$A$$ （涉及 $$G(0)$$ ）保持 Sobolev 范数但不增加（例如， $e^{tA}: H^s \to H^s$ 有界但非 $H^{s+\delta}$ ）。因此，无法像抛物型那样在迭代中提升正则性。
能量估计结构：尽管 Hamiltonian PDE 有守恒能量，但线性化 $$A(u)$$ 的耗散性可能依赖于状态 $$u$$ 。需要构造状态相关的能量范数（即 $L_\pm(u)$ ）来获得类似指数二分性的估计。

此前，对于固定区域上的不可压缩 Euler 方程，文献 [53] 通过 Lagrangian 坐标将方程转化为拟线性形式，并利用能量估计和压力投影构造了局部稳定和不稳定流形。但自由边界问题（如水波）在 Eulerian 坐标下更自然，却面临 Dirichlet-to-Neumann 算子的非局部性，分析更复杂。一般地，对于无正则性增益的拟线性 PDE，不变流形理论长期缺失，这正是本文要解决的核心问题。

论文通过发展一套基于能量估计的框架，在满足特定假设（包括状态相关的能量形式和耗散性）下，证明局部稳定和不稳定流形的存在、唯一及光滑性。主要动机来自水波系统，但框架适用于 broad class，包括非线性 Schrödinger 方程、非线性波方程、MMT 模型（Majda–McDonough–Tsay 湍流模型）以及某些梯度型 PDE。

水波系统的函数空间与边值条件

对于水波，自然的状态空间为 Sobolev 空间 $H^s(\mathbb{R}^d)$ （ $$s > d/2$$ 以保证乘积有意义），通常取 $X = H^s \times H^{s-1/2}$ （因为 $\psi$ 的正则性比 $\eta$ 低半阶）。底部固边界和自由界面条件已融入算子 $G(\eta)$ 的定义。系统 (1) 可写为 Hamiltonian 形式 $u_t = J \nabla H(u)$ ，其中 $$J$$ 是辛结构， $$H$$ 是能量函子（动能加势能）。平衡点 $$u=0$$ 对应 $\eta=0,\psi=0$ 。线性化算子 $A = J \nabla^2 H(0)$ 的谱结构决定了不变分解：通过选择 $$s$$ 足够大，可使 $$A$$ 的谱在虚轴上无点，且存在 $\lambda_\pm$ 使得 $X = X_+ \oplus X_-$ 满足指数二分性（例如，对于短波，表面张力使所有模式稳定；对于长波，重力产生零模式或不稳定模式，需排除平移对称导致的临界情况）。

论文中的函数空间 $$X^r$$ 是 $$A$$ 的 Sobolev 型域，例如 $X^r = H^{s+r} \times H^{s+r-1/2}$ 。非线性 $$F$$ 需在 $$X^n$$ 上具有 $$C^1$$ 光滑性（但可能允许从 $$X^n$$ 到 $X^{n-1}$ 的正则性损失），这是拟线性 PDE 局部适定性理论中的常见现象。

2. 主要定理

论文的主要成果包括抽象框架下的一般定理和对水波系统的具体应用。以下详述主要定理及其假设与结论。

一般框架下的定理（Theorem 2.4 和 Theorem 2.9）

考虑抽象 PDE：

u_t = F(u), \quad F(0)=0,

在实 Hilbert 空间 $$X$$ 中， $$X$$ 具有 Sobolev 尺度 $X^r = \text{Dom}(A^r)$ ， $r \ge 0$ ，其中 $$A = DF(0)$$ 。假设存在闭子空间 $X_+, X_- \subset X$ 使得 $X = X_+ \oplus X_-$ ，且 $e^{tA}$ 满足指数二分性：

\sup_{t \le 0} e^{-\lambda_+ t} \|e^{tA}|_{X_+}\| < \infty, \quad \sup_{t \ge 0} e^{-\lambda_- t} \|e^{tA}|_{X_-}\| < \infty,

其中 $\lambda_+ > \lambda_-$ 。记 $\Pi_\pm$ 为到 $X_\pm$ 的投影。

进一步假设存在整数 $n \ge 2$ ，使得：

$A(u) = DF(u) \in \mathcal{L}(X^r, X^{r-1})$ 对 $1 \le r \le n'$ ， $$n' = n$$ 或 $$n-1$$ ，且在 $X^{n-1}$ 上 $$C^1$$ 光滑；
存在等价能量形式 $L_\pm(u): X_\pm \to (X_\pm)^*$ ， $$C^1$$ 依赖于 $u \in X^{n-1}$ ，满足对某个 $\lambda_\pm$ 和所有 $v \in X^1$ ： $\langle L_-(u)A(u)v, v \rangle \le -\lambda_- \|v\|_{X^1}^2, \quad \langle L_+(u)A(u)v, v \rangle \ge \lambda_+ \|v\|_{X^1}^2.$ 这等价于 $$A(u)$$ 在 $$X_-$$ 上耗散、在 $$X_+$$ 上反耗散，估计一致于 $$u$$ 在 $X^{n-1}$ 中有界。

非线性项在分解 $$F(u) = (A_+(u)u_+ + f_+(u), A_-(u)u_- + f_-(u))$$ 下满足：存在 $k \ge 1$ 、 $$R_0>0$$ 、 $$C_f>0$$ ，使得对 $v \in X^{k-1}(R_0)$ ，

|f_\pm(v)|_{X^k} \le C_f |v|_{X^{k-1}}, \quad |f_\pm(v_1) - f_\pm(v_2)|_{X^{k-1}} \le C_f |v_1 - v_2|_{X^{k-1}}.

（此处 $$X^r(R_0)$$ 表示 $$X^r$$ 中半径为 $$R_0$$ 的球。）

定义 $$ \Sigma_+ = \left{ \lambda \in (\lambda_-, \lambda_+) ,\middle|, L_1(\lambda) = \frac{2(k+1)C_0^{2(k+1)}C_f}{\lambda - \lambda_- - 2(k+1)C_0^{2(k+1)}M_0} + \frac{2(k+1)C_0^{2(k+1)}C_f}{\lambda_+ - \lambda - 2(k+1)C_0^{2(k+1)}M_0} < 1 \right}, ] 其中 $$C_0$$ 来自 $L_\pm(u)$ 的等价性常数， $$M_0$$ 是某个固定正数。 $\Sigma_+$ 非空当 $\lambda_+ - \lambda_-$ 足够大相对于 $$C_f$$ 。

则主要定理为：

定理 (Theorem 2.4，主要一般结果)。假设 $\lambda \in \Sigma_+$ 。则存在 $\delta > 0$ 和映射 $q^+: X_+^n(\delta) \to X_-^n$ ，使得：

存在性与唯一性： $$q^+(0)=0$$ ， $$q^+$$ 是 Lipschitz 连续：存在常数 $l \in (L_1(\lambda), 1]$ 使得 $|q^+(v)|_{X_-^n} \le \frac{2l C_0^{2(k+1)}}{1-L_1(\lambda)} |v|_{X_+^n}.$ 若 $$F$$ 是 $C^{m,1}$ 光滑（ $m \ge 1$ ），则 $q^+ \in C^{m,1}(X_+^n(\delta), X_-^{n-m-1})$ ，即光滑性有 $$m+1$$ 阶损失。
局部不稳定流形：$W^+ = \text{graph}(q^+)$ 满足：
- $$T_0 W^+ = X_+^n$$ （在零点切于 $$X_+^n$$ ）；
- $$W^+$$ 局部不变：若 $u(0) \in W^+$ ，则存在唯一解 $$u(t)$$ 定义于 $t \le 0$ ，满足 $u(t) \in W^+$ 对所有 $t \le 0$ ，且存在 $$M_* > 0$$ 使得 $|u(t)|_{X^n} \le M_* \delta e^{\lambda_0 t}, \quad \lambda_0 = \max\{0, \lambda_-\} \text{ 或 } \lambda_+/2.$
- 解满足指数衰减：对任意 $\mu \in (\lambda_-, \lambda_+)$ ， $\sup_{t \le 0} e^{-\mu t} |u(t)|_{X^n} < \infty.$
渐近性：若解 $$u(t)$$ 满足 $\sup_{t \le 0} e^{-\lambda t}|u(t)|_{X^n} < \infty$ 对某 $\lambda \in (\max\{0,\lambda_-\}, \lambda_+)$ ，则存在 $t_0 \le 0$ 使得对所有 $t \le t_0$ 有 $u(t) \in W^+$ 。
稳定流形：若 $\lambda < \min\{0, \lambda_+\}$ ，则存在类似映射 $q^-: X_-^n(\delta) \to X_+^n$ 和局部稳定流形 $W^- = \text{graph}(q^-)$ ，满足对偶性质（ $t \ge 0$ 的指数衰减）。

定理的精确假设在附录 A 的 (D.1)-(D.4) 和 (B.5) 中给出，以上为概要。

水波系统的定理（Theorem 4.2）

论文在 Section 4 具体应用到具有表面张力的无旋水波系统。考虑单相情形（恒定底部深度 $$d>0$$ ），状态空间取为 $X = H^s(\mathbb{R}) \times H^{s-1/2}(\mathbb{R})$ ， $$s > d/2 + 1$$ （足够大以确保局部适定性和非线性项的有界性）。平衡点为静水平面 $(\eta,\psi)=0$ 。

定理 4.2。存在 $$s_0 > d/2 + 1$$ ，使得对任意 $s \ge s_0$ ，以下成立：设 $\lambda_\pm$ 为线性化算子 $$A = DF(0)$$ 在 $$X$$ 上的谱间隙（存在不变分解 $X = X_+ \oplus X_-$ 满足指数二分性），且 $\lambda_+ > \max\{0, \lambda_-\}$ 。则存在 $\delta > 0$ 和映射 $q^\pm$ 定义在 $X_\pm^s(\delta)$ 上，使得 $W^\pm = \text{graph}(q^\pm)$ 是局部稳定/不稳定流形，满足 Theorem 2.4 的所有性质（存在性、唯一性、指数衰减、渐近性），且光滑性损失与一般定理一致。

证明的关键在于验证水波系统满足一般框架的假设：

不变分解：通过傅里叶分析，线性化水波方程可解耦为独立模式，每个模式对应常系数 ODE。选择 $$s$$ 足够大可使低正则性模式（如零频率）被排除，从而 $\lambda_\pm$ 远离零，确保指数二分性。
能量形式 $L_\pm(u)$ ：取自 Hamiltonian 能量泛函在 $$u$$ 处的 Hessian 的限制，或通过直接能量估计构造。例如，能量 $E(\eta,\psi)$ 的二次变分在 $$u$$ 处给出双线性形式，其在 $X_\pm$ 上正负定，且关于 $$u$$ 连续可微。
Dirichlet-to-Neumann 算子 $G(\eta)$ 的估计： $G(\eta)$ 作为 $\eta$ 的拟微分算子，其 Sobolev 范数满足 $\|G(\eta)\|_{H^s \to H^{s-1}} \le C(1+\|\eta\|_{H^{s}})^p$ ，且关于 $\eta$ 光滑。这通过傅里叶积分算子理论或显式符号计算证明。
非线性项估计： $f_\pm(u)$ 涉及 $G(\eta)\psi$ 和 $\eta$ 的导数，利用 $G(\eta)$ 的估计和乘积法则，可证其满足局部 Lipschitz 条件（允许从 $$H^s$$ 到 $H^{s-1}$ 的损失）。

因此，水波系统落入一般框架，流形存在性得证。

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

论文的核心贡献在于发展了一套结合 Lyapunov-Perron 积分方程方法与能量估计技术的框架，专门处理拟线性 PDE 中由非线性导致的正则性损失以及线性部分缺乏光滑化效应的难题。以下是详细的技术分析。

a) 核心分析学工具

Lyapunov-Perron 方法：经典的不变流形构造方法。对于不稳定流形，寻找图 $$v_- = q^+(v_+)$$ 使得过 $$(v_+, q^+(v_+))$$ 的解满足 $v_-(t) \to 0$ 当 $t \to -\infty$ 。这等价于求解积分方程：
$v_+(t) = e^{tA_+}v_+^0 + \int_0^t e^{(t-\tau)A_+} f_+(v(\tau))\,d\tau, \quad v_-(t) = \int_{-\infty}^t e^{(t-\tau)A_-} f_-(v(\tau))\,d\tau.$
通常通过固定点 argument 定义 $$q^+$$ ：给定 $$v_+^0$$ ，寻找 $v \in Y$ （适当函数空间）使得 $$v = T(v_+^0, v)$$ ，其中 $$T$$ 是上述右侧。但拟线性情况下， $e^{tA}$ 不光滑化，且 $$f(v)$$ 可能降低正则性，传统迭代失败。论文的创新在于在加权范数下工作，并利用状态依赖的能量估计控制非自治线性算子。
能量估计与非自治线性算子控制：
- 假设 (C.2) 要求 $\lambda_\pm A_\pm(v)$ 是满射且满足 dissipativity（2.3）： $\langle L_\pm(v)A_\pm(v)w, w \rangle \le \mp \lambda_\pm \|w\|_{X^1}^2$ 。结合 $L_\pm(v)$ 的等价性（ $C^{-1}\|w\|_{X^\pm} \le \|w\|_{L_\pm(v)} \le C\|w\|_{X^\pm}$ ）和 Lumer-Phillips 定理，推出 $A_\pm(v)$ 生成 $$C^0$$ 半群 $e^{sA_\pm(v)}$ ，且 $\|e^{sA_\pm(v)}\|_{\mathcal{L}(X^r)} \le C_0^{2(r+1)} e^{\lambda_\pm s}, \quad s \lessgtr 0.$ 注意：无光滑化—— $$r$$ 不变，且常数 $$C_0$$ 可能依赖于 $$v$$ 。
- 对于非自治系统 $$v$$ 固定，演化算子 $U_\pm(t,t_0)$ 由 $A_\pm(v(t))$ 生成。附录 A.1 证明关键估计（公式 (2.8)）： $\|U_\pm(t,t_0)\|_{\mathcal{L}(X^r)} \le C_0^{2(r+1)} e^{\lambda_\pm (t-t_0) + (r+1)C_0^{2(r+1)} \int_{t_0}^t \|v'(\tau)\|_{X^{k-1}} d\tau},$ 其中 $r \le k$ 。这显式依赖于 $\|v'\|_{X^{k-1}}$ 的积分，允许 $A_\pm(v(t))$ 随时间变化，但需控制 $$v'$$ 的大小。这是构造迭代压缩性的基础。
加权 Sobolev 范数：定义
$|v|_{\pm, r, \lambda} = \sup_{t \le 0} e^{-\lambda t} \|v(t)\|_{X^r}$
对于不稳定流形（ $t \le 0$ ）。选择 $\lambda \in (\lambda_-, \lambda_+)$ 平衡 $e^{tA_+}$ 的增长和 $$f_+$$ 的压缩。加权范数自然地吸收线性部分的指数增长。
函数空间与子集：设
$Y = \{ v(\cdot) \in W^{1,\infty}((-\infty,0], X^{k-1}) \mid v(t) \in X^k(R_0), \forall t \le 0 \},$
其中 $$X^k(R_0)$$ 是 $$X^k$$ 中的闭球。这确保解在时间趋于 $-\infty$ 时保持在正则球内。然后定义子集
$\Gamma_{\lambda, M_0, \varepsilon} = \{ v \in Y \mid |v|_{k,\lambda} \le M_0\varepsilon, \ \|\partial_t v\|_{k-1,0} \le M_1\varepsilon \},$
其中 $M_0, M_1, \varepsilon$ 需满足精细不等式（2.11)-(2.12）以确保 $$T$$ 的映射性和压缩性。

b) 解决技术性难题的主要技巧

处理非线性项的正则性损失：
- 在迭代中， $$f(v)$$ 的估计通常只有 $X^{k-1}$ 范数，但 $$v$$ 需要 $$X^k$$ 正则性来闭合估计。论文的技巧是：在加权范数下，利用 $$f(v)$$ 的 Lipschitz 常数允许 $X^{k-1}$ 损失，但通过控制 $\|\partial_t v\|_{k-1,0}$ （由 $M_1\varepsilon$ 限制），在 $U_\pm$ 的估计 (2.8) 中，指数因子 $e^{(r+1)C_0^{2(r+1)} \int_{t_0}^t \|v'\|_{k-1} d\tau} \le e^{(r+1)C_0^{2(r+1)} M_1 |t-t_0|}$ 被线性因子吸收，只要选择 $\lambda$ 使得 $\lambda_\pm - (r+1)C_0^{2(r+1)}M_1 > 0$ 。这正是 $\Sigma_\pm$ 条件的来源。
避免直接迭代 $$f(v)$$ 在 $$X^k$$ 中：
- 传统 Lyapununov-Perron 方法在 $$X^k$$ 中迭代 $$v$$ ，但这里 $$f(v)$$ 可能不在 $$X^k$$ 中。论文改为在 $$Y$$ 中迭代 $$v$$ ，其中 $v \in W^{1,\infty}(X^{k-1})$ ，而 $f(v) \in X^{k-1}$ 足够（因为 $A_\pm(v)$ 映射 $$X^k$$ 到 $X^{k-1}$ ）。对于 $$v_+$$ 的初值 $$v_+^0$$ ，要求 $v_+^0 \in X_+^k$ ，但 $$v_+(t)$$ 通过 $e^{tA_+}$ 保持在 $$X^k$$ ，而 $$v_-(t)$$ 由积分定义在 $X^{k-1}$ 。通过选择 $$k$$ 和常数，使 $$v_-(t)$$ 实际属于 $$X^k$$ （见 (2.13)）。这利用了 $$f_-$$ 的估计和 $$A_-$$ 的耗散性。
压缩常数 $$l$$ 的计算与 $\Sigma_\pm$ 条件：
- 由 $$T$$ 的估计（Lemma 2.2）： $|T^-(v_+^0, v)|_{k,\lambda} \le \frac{C_0^{2(k+1)}C_f}{\lambda_- - \lambda} |v|_{k,\lambda} + \cdots, \quad |T^+(v_+^0, v)|_{k,\lambda} \le C_0^{2(k+1)}|v_+^0|_{X_+^k} + \frac{C_0^{2(k+1)}C_f}{\lambda_+ - \lambda} |v|_{k,\lambda}.$ 通过选择 $$M_0$$ 和 $\lambda$ 使得系数和小于 $$l < 1$$ ，得到压缩性。 $\Sigma_+$ 确保存在 $\lambda$ 使得 $l(\lambda) = \frac{2(k+1)C_0^{2(k+1)}C_f}{\lambda - \lambda_- - 2(k+1)C_0^{2(k+1)}M_0} + \frac{2(k+1)C_0^{2(k+1)}C_f}{\lambda_+ - \lambda - 2(k+1)C_0^{2(k+1)}M_0} < 1.$ 这要求 $\lambda_+ - \lambda_-$ 显著大于 $$C_f$$ ，即谱间隙足够大。这不同于半线性情形（压缩仅取决于 $$C_f$$ ，与 $\lambda_\pm$ 无关）。
光滑性损失的不可避免性：
- 附录 A.3.1 证明：对于拟线性 PDE，若 $F \in C^m(X^n, X^{n-1})$ ，则解映射的 Fréchet 导数 $\partial_u u(t)$ 可能仅连续且取值于 $X^{n-m-1}$ ，而非 $$X^n$$ 。这源于线性化方程中对 $$u$$ 的导数项： $\partial_t(\partial_u u) = A(u)\partial_u u + D A(u)[u]\cdot u + D f(u)$ ，其中 $$D A(u)$$ 涉及 $$u$$ 的高阶导数，每次迭代损失一阶。因此，流形映射 $q^\pm$ 的光滑性自然损失 $$m+1$$ 阶（Theorem 2.1(1) 和 Remark 2.4）。这与半线性情形（ $q^\pm$ 与 $$F$$ 同光滑）形成对比，是拟线性问题的本质特征。论文在定理中明确刻画了这一损失。
水波系统的额外技术：
- Dirichlet-to-Neumann 算子 $G(\eta)$ 的精细估计： $G(\eta)$ 是 $\eta$ 的拟微分算子，其符号为 $|k|\sqrt{1+|k|^2}\tanh(d\sqrt{1+|k|^2})$ 乘以 $(1+|k|^2\widehat{\eta}(k))/(1-|k|^2\widehat{\eta}(k))$ 等因子。需证明：对 $$s > d/2+1$$ ，存在 $$C(s)$$ 使得 $\|G(\eta)\psi\|_{H^{s-1}} \le C(1+\|\eta\|_{H^{s}}^p)\|\psi\|_{H^{s}}, \quad \|D G(\eta)[\xi]\psi\|_{H^{s-1}} \le C(1+\|\eta\|_{H^{s}}^q)(\|\xi\|_{H^{s}} \|\psi\|_{H^{s}} + \|\eta\|_{H^{s}}\|\xi\|_{H^{s-1}}\|\psi\|_{H^{s}}).$ 这通过分析符号的导数和利用 $\eta$ 的小性（在局部流形构造中 $\|\eta\|_{H^s} \ll 1$ ）完成。
- 非线性项 $f_\pm$ 的估计：例如， $$f_-$$ 包含 $-G(\eta)\psi$ 项，需估计 $\|G(\eta)\psi\|_{H^{s-1}}$ 对 $\|\eta\|_{H^s}, \|\psi\|_{H^{s-1/2}}$ 有界。由于 $G(\eta): H^{s} \to H^{s-1}$ 有界，且 $\psi$ 在 $H^{s-1/2}$ ，需插值或乘积估计： $\|G(\eta)\psi\|_{H^{s-1}} \lesssim \|G(\eta)\|_{H^{s} \to H^{s-1}} \|\psi\|_{H^{s}}$ ，但 $\psi$ 可能不在 $H^{s}$ 。实际上，在分解 $v = (v_+, v_-) \in X_+^s \times X_-^s$ 后， $v_+ \in X_+^s \subset H^{s}$ ， $v_- \in X_-^s \subset H^{s-1/2}$ ，且 $$f_+$$ 作用于 $$(v_+, v_-)$$ 可能涉及 $$v_-$$ 的导数。必须仔细利用分解的兼容性（例如， $$X_-^s$$ 可能包含在 $H^{s-1/2}$ 但具有额外衰减保证）。论文通过局部坐标变换和能量估计，确保 $f_\pm$ 满足 (C.3) 的类似形式。
- 能量形式 $L_\pm(u)$ 的构造：对于水波， $L_\pm(u)$ 可取自能量泛函 $$E(u)$$ 的 Hessian 在 $$u$$ 处的限制。例如， $E(\eta,\psi)$ 是二次泛函加非线性扰动，在 $$u=0$$ 处 Hessian 给出 $L_\pm(0)$ 。对于非零 $$u$$ ，需修正以确保耗散性（2.3）。这通过将 $$A(u)$$ 与 $J \nabla^2 E(u)$ 联系，并利用 $$E(u)$$ 的非线性结构实现。
非自治估计的精细依赖：
- 附录 A.1 的估计 (2.8) 是技术核心。其证明利用变分公式和对 $L_\pm(v)$ 的时间导数估计。设 $$w(t)$$ 满足 $\partial_t w = A(v(t))w$ ，则 $\frac{d}{dt} \|w(t)\|_{L_\pm(v(t))}^2 = 2\langle L_\pm(v) A(v)w, w \rangle + \langle \partial_t L_\pm(v) w, w \rangle.$ 第一项由 dissipativity 控制，第二项用 $DL_\pm$ 的界和 $\|v'\|_{X^{k-1}}$ 估计。Gronwall 不等式导出 (2.8)。这证明了对任意 $v \in Y$ ， $U_\pm$ 的有界性。
局部流形的延拓与最大性：
- 定理 2.1(3) 表明若解 $$u(t)$$ 在某正时间 $$t$$ 仍属于 $X^k(\delta)$ ，则流形定义可延拓。这由 $$h$$ 的连续性和解的唯一性保证。
- 渐近性（定理 2.1(5)）的证明：若 $$u(t)$$ 指数衰减，考虑足够负的 $$t$$ ，使 $v(t) \in X^k(\varepsilon)$ ，然后利用流形的吸引性（由压缩性保证）证明存在 $$t_0$$ 使得 $v(t) \in W^+$ 对 $t \le t_0$ 。这需要比较 $$v_-(t)$$ 与 $$q^+(v_+(t))$$ ，利用 $$T$$ 的压缩性。

核心工具总结

状态依赖的能量范数 $L_\pm(u)$ ：替代固定范数，适应非线性变化，提供耗散性估计。
加权范数 $|\cdot|_{\pm,r,\lambda}$ ：吸收线性指数增长，压缩非线性影响。
精细常数选择：通过 $\Sigma_\pm$ 条件平衡谱间隙 $\lambda_+ - \lambda_-$ 、非线性强度 $$C_f$$ 和正则性损失 $(r+1)C_0^{2(r+1)}M_1$ 。
非自治线性估计：将时间依赖的 $$A(v(t))$$ 的影响累积到指数因子中，通过控制 $$v'$$ 的范数。
光滑性损失分析：通过研究解对初值导数的方程，揭示每迭代一次损失一阶，从而确定流形光滑性。

这些工具共同解决了拟线性 PDE 中因非线性正则性损失和线性无光滑化而导致的传统方法失效问题。

4. 与之前工作的比较

不变流形理论在动力系统中历史悠久。此前工作主要针对以下场景：

有限维 ODE：当线性部分 $$A$$ 有指数二分性（即存在不变子空间 $X_\pm$ 使得 $e^{tA}|_{X_\pm}$ 指数增长/衰减）时，局部稳定/不稳定流形存在唯一，且若 $$F$$ 是 $$C^m$$ ，则流形是 $$C^m$$ 图（经典结果，如 [64, 48, 14]）。
半线性 PDE： $$u_t = Au + f(u)$$ ，其中 $$A$$ 生成解析或光滑半群（例如，抛物型方程如热方程、Navier-Stokes，或双曲型如 Klein-Gordon 方程）。此时 $e^{tA}$ 光滑化，允许在 $$X^n$$ 中直接迭代 Lyapunov-Perron 积分，流形光滑性通常匹配 $$f$$ 的光滑性（参见 [41, 43, 19, 6, 22, 18, 73]）。
抛物型拟线性 PDE： $$u_t = A(u)u + f(u)$$ ，其中 $$A(u)$$ 满足抛物型条件（如 $A(u) = \text{div}(a(u)\nabla)$ ），则 $e^{tA(u)}$ 提供极大正则性（从 $$L^p$$ 到 $C^\infty$ ）。不变流形可构造，光滑性可能略有损失但可控（例如 [57, 49, 27, 65, 72]）。
固定区域上的 Euler 方程：[53] 在 Lagrangian 坐标下，将不可压缩 Euler 方程写为拟线性形式，并利用能量估计和压力投影构造局部稳定和不稳定流形。但该方法依赖于区域的几何和 Lagrangian 变换，对于自由边界问题（如水波）不直接适用。

论文的突破在于将不变流形理论扩展到无正则性增益的拟线性 PDE，特别是 Hamiltonian PDE。主要改进与创新包括：

放宽线性部分的假设：
- 不要求 $$A$$ （或 $$A(u)$$ ）生成光滑半群或提供正则性增益。相反，仅假设存在状态相关的能量形式 $L_\pm(u)$ 使得线性化满足 dissipativity 条件 (2.3)。这允许线性部分为常系数但非光滑化算子，例如水波中的 $$G(0)$$ （保持 Sobolev 范数但不增加）。这适用于许多 Hamiltonian 系统，其中线性色散算子的半群不改善正则性（如 KdV 的 $e^{t\partial_x^3}$ 保持 $$H^s$$ ，Schrödinger 的 $e^{it\Delta}$ 甚至可能降低正则性）。
系统处理正则性损失：
- 传统方法要求非线性项 $$f(u)$$ 映射到同一正则空间（如 $$X^n$$ ），但拟线性系统中 $$f(u)$$ 常含 $$u$$ 的导数，导致正则性损失（例如， $$f(u)$$ 可能只属于 $X^{n-1}$ ）。论文通过加权范数和精细常数选择，允许 $$f$$ 从 $$X^k$$ 到 $X^{k-1}$ 的损失，并调整迭代中的压缩常数。具体地，在 Lemma 2.2 的证明中，利用 $$f$$ 的 Lipschitz 估计时，其作用于 $v(\tau)$ 的范数用 $|v|_{k,\lambda}$ 控制，而 $|v|_{k,\lambda}$ 的递归不等式通过选择 $\lambda \in \Sigma_+$ 确保压缩性。这比要求 $$f$$ 光滑更弱，更贴近拟线性 PDE 的实际困难。
统一框架适用于多种 PDE：
- 论文明确展示框架适用于：
  - 非线性 Schrödinger 方程（NLS）： $i u_t + \Delta u + |u|^{p-1}u = 0$ 。虽然 NLS 是半线性，但线性部分 $i\Delta$ 不 dissipative，需通过复分解分离稳定/不稳定方向（例如，聚焦 NLS 在零点附近，正频率对应不稳定，负频率对应稳定？实际上，线性 Schr¨odinger 算子 $i\Delta$ 的谱在虚轴，无指数二分性。论文可能考虑空间周期域或通过限于实/虚部构造不变分解？需具体分析。论文 Subsection 3.1 提到 Hamiltonian PDE，可能 NLS 的平衡点包括孤立子，但零点附近线性化本质中心。或许论文考虑特定情况如四阶 NLS 或带耗散？从摘要“including nonlinear Schrödinger equations”看，可能 NLS 的线性部分有耗散方向（如在某些 Sobolev 空间中），或通过考虑复结构分解。无论如何，框架的可适用性得到论证。）
  - 非线性波方程： $u_{tt} - \Delta u + f(u) = 0$ ，可写为一阶系统，线性部分有双曲型谱，部分方向稳定/不稳定。
  - MMT 模型：Majda–McDonough–Tsay 模型用于湍流，是拟线性系统。
  - 梯度型 PDE：如某些反应扩散方程的梯度流，线性部分可能 dissipative。
- 这展示了框架的广泛性，不限于特定结构。
可能适用于不适定问题：
- 论文指出，不假设 $$A_+(v)$$ 生成 forward 半群（ $s \ge 0$ ），因此方法可能应用于不适定演化系统，例如将椭圆 PDE（如 Laplace 方程在通道中）视为时间演化（通过解析延拓）。这扩展了不变流形理论的传统范围（通常假设适定性以定义 $e^{tA}$ ）。
精确刻画光滑性损失：
- 不仅证明流形存在，还详细分析光滑性损失。定理表明：若 $F \in C^{m,1}(X^n, X^{n-1})$ ，则 $q^\pm \in C^{m,1}$ 但目标空间正则性损失 $$m+1$$ 阶（从 $X_\mp^{n-m-1}$ ）。附录 A.3.1 证明解对初值的光滑依赖有类似损失： $\partial_u u(t)$ 可能仅属于 $X^{n-m-1}$ 。这通过研究变分方程 $\partial_t(\partial_u u) = A(u)\partial_u u + D A(u)[u] u + D f(u)$ 并迭代估计来证明，每次 $$D A(u)$$ 引入一阶损失。因此，流形光滑性损失是不可避免的，与系统的拟线性本质一致。这与半线性情形（流形与 $$F$$ 同光滑）形成鲜明对比。
水波系统的首次严格构造：
- 据所知，对于自由边界水波问题，在 Eulerian 坐标下建立局部稳定和不稳定流形是首次。此前 [53] 处理固定区域 Euler 方程，而水波涉及移动边界和 Dirichlet-to-Neumann 算子的非局部性，分析难度大。论文通过验证水波满足能量假设（特别是构造 $L_\pm(u)$ 和估计 $G(\eta)$ ）克服了这些困难。这为研究水波的局部动力学（如不稳定性、异宿轨道、 breather 解）提供了核心工具。
方法论的贡献：
- 将常用于局部适定性的能量估计技术与 Lyapunov-Perron 方法结合，并引入状态依赖的能量范数。这为处理缺乏光滑化的拟线性 PDE 提供了新范式。
- 非自治线性估计（附录 A.1）独立有趣，可能用于其他非自治问题。
- 框架的设计尽可能直接适用于具体 PDE，假设 (D.1)-(D.4) 和 (B.5) 与局部适定性理论中的常见假设对齐（如能量估计、光滑依赖性），降低了应用门槛。

与特定工作的对比

与 [53] 对比：[53] 处理固定区域不可压缩 Euler 方程，使用 Lagrangian 坐标将方程转化为拟线性形式 $$v_t = A(v)v$$ ，并利用能量估计和压力投影构造流形。论文的方法在 Eulerian 坐标下工作，直接处理 $$u_t = F(u)$$ ，无需 Lagrangian 变换。这对于自由边界问题更自然，但需处理非局部算子（如 $G(\eta)$ ）。论文的能量形式 $L_\pm(u)$ 可比作 Lagrangian 中的 Riemann 度量。
与半线性理论对比：半线性理论中， $e^{tA}$ 光滑化，迭代在 $$X^n$$ 闭合，流形光滑性通常为 $$C^m$$ 若 $f \in C^m$ 。论文中，由于 $$A(u)$$ 的拟线性性， $e^{tA(u)}$ 无光滑化，且 $$f(u)$$ 正则性损失，导致 $q^\pm$ 光滑性降低。这是拟线性问题的根本差异。
与抛物型拟线性理论对比：抛物型方程（如 Navier-Stokes）的线性部分有极大正则性， $e^{tA(u)}$ 立即光滑化，因此不变流形可构造且光滑性较好。论文的方法不依赖极大正则性，仅需耗散性估计，适用范围更广。

总之，论文通过创新性的能量估计结合加权范数，成功突破了拟线性 Hamiltonian PDE 中不变流形理论的长期障碍，为水波等系统的动力学研究奠定了坚实基础。

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

1. 衰减密度不确定性原理在薛定谔方程可观测性中的应用

Application of uncertainty principles for decaying densities to the observability of the Schrödinger equation

作者: Kévin Le Balc’h, Jiaqi Yu

类别: math.AP

发布日期: 2026-02-20

arXiv ID: 2602.18371v1

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

摘要翻译

中文标题: 衰减密度不确定性原理在薛定谔方程可观测性中的应用

摘要翻译: 本文研究在欧几里得空间中讨论的薛定谔方程。我们针对相对于衰减密度具有厚度的可测集，证明了相应的可观测性不等式。证明过程依赖于适配于衰减密度的定量不确定性原理，特别是舒宾（Shubin）、瓦基利安（Vakilian）、沃尔夫（Wolff）和科夫斯基金（Kovrijkine）所建立的相关结果。

2. 具有奇异各向异性核的卷积算子的非局部到局部 $$L^p$$ 收敛性

Nonlocal-to-local $$L^p$$ -convergence of convolution operators with singular, anisotropic kernels

作者: Helmut Abels, Christoph Hurm, Patrik Knopf

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

发布日期: 2026-02-20

arXiv ID: 2602.18183v1

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

摘要翻译

中文标题: 具有奇异各向异性核的卷积算子的非局部到局部 $$L^p$$ 收敛性

摘要翻译: 我们研究了具有奇异、可能各向异性的非局部卷积型算子。主要目标是建立并量化这类算子在核以适当方式集中在原点时，向具有自然边界条件的局部微分算子的非局部到局部收敛性。此类收敛结果为数学模型的物理论证提供了有效工具,特别适用于无法直接从微观规律导出目标局部微分算子的情况。本研究通过允许核具有更强的奇异性(类似于分数阶拉普拉斯算子的奇异性)、具有各向异性与非局部化特性,并在一般 $$L^p$$ 空间中证明强收敛性及给出显式收敛速率,显著拓展了既往研究成果。

3. 分数阶拉普拉斯算子的弱奇性表述

Weakly-singular formulation of the Fractional Laplacian operator

作者: Oscar P. Bruno, Sabhrant Sachan

类别: math.AP, math.NA

发布日期: 2026-02-20

arXiv ID: 2602.17970v1

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

摘要翻译

中文标题: 分数阶拉普拉斯算子的弱奇性表述

摘要翻译: 本文提出了分数阶拉普拉斯算子 $$(-Δ)^s$$ 在 n 维空间 ( $n \ge 1$ ) 中的一种新表述。该表述将 $$(-Δ)^s$$ 表示为经典拉普拉斯微分算子与一个弱奇异积分算子的复合，从而可将分数阶拉普拉斯算子的狄利克雷问题归结为包含体积分与边界积分算子的弱奇异积分方程。这种重新表述非常适合高效且精确的数值实现。尽管相关高阶算法的完整描述留待后续研究，本文仍提供了若干数值算例，以展示所提方法所达到的高精度与计算效率。

不相关论文（基本信息）

1. 分数阶亥姆霍兹方程的极限吸收原理与辐射条件

Limiting Absorption Principle and Radiation Condition for the Fractional Helmholtz Equation

作者: Dana Zilberberg, Fioralba Cakoni, Michael S. Vogelius

类别: math.AP

发布日期: 2026-02-20

arXiv ID: 2602.18387v1

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

摘要: We investigate elliptic fractional equations in the whole space, involving zero order perturbations of the fractional Laplacian $$(-Δ)^s$$ , $$0<s<1$$ . Our main objective is to determine appropriate radiation conditions at infinity that ensure existence and uniqueness of solutions to the fractional type Helmholtz equation. Extending classical scattering theory for the Helmholtz equation, we introduce and analyze suitable Sommerfeld type radiation conditions for fractional orders. A central contribution is the explicit computation of the outgoing free space Greens function for the operator $(-Δ)^s-k^{2s}$ , for all $$0<s<1$$ , any dimension and $$k>0$$ , obtained via contour integration and a limiting absorption principle. We show that its asymptotic behavior at infinity coincides with a rescaled version of the classical Helmholtz fundamental solution, thereby justifying the standard Sommerfeld radiation condition for compactly supported sources. In addition, using resolvent estimates and a limiting absorption framework, we establish existence and uniqueness of outgoing solutions for compactly supported data, and for weighted sources. We further derive a convolution representation of the solution in terms of the outgoing fundamental solution. For inhomogeneous media with compactly supported perturbations, we reformulate the problem as a Lippmann Schwinger integral equation of Fredholm type and prove unique solvability away from a discrete set of frequencies. Our analysis provides a rigorous foundation for scattering theory of fractional Helmholtz operators and offers a framework suitable for numerical implementation of these nonlocal wave propagation models.

2. 旋转不变泊松方程在模型流形上的超定问题

Overdetermined problems for the rotationally invariant Poisson equation in model manifolds

作者: Antonio Greco, Marcello Lucia, Pieralberto Sicbaldi

类别: math.AP, math.DG

发布日期: 2026-02-20

arXiv ID: 2602.18289v1

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

摘要: We present rigidity results for overdetermined problems associated to the rotationally invariant Poisson equation $-Δ_{g_\mathcal{M}} u = f(r)$ in a model manifold $\mathcal{M} = [0,S) \times_h \mathbb S^{N-1}$ with warping function $$h$$ . The variable $$r$$ ranges in the interval $$[0,S)$$ , whose endpoint $$S$$ is positive and possibly infinite. The first part of the paper deals with the problem

\begin{array}{ll} -Δ_{g_\mathcal{M}} {u}=f(r) &\mbox{in $Ω$}, u=\varphi(r) &\mbox{on $\partial Ω$}, \frac{\partial u}{\partial ν} = κ(r) &\mbox{on $\partial Ω$}, \end{array}

where $Ω\subset \mathcal{M}$ is a bounded domain containing the point $O \in \mathcal{M}$ corresponding to $$r = 0$$ , $$ν$$ is the exterior unit normal vector on $\partial Ω$ , and $$f$$ , $\varphi$ , $$κ$$ are three prescribed functions. In the second part of the paper, we consider a similar overdetermined problem for the exterior Bernoulli problem in a domain $Ω\setminus \overline B_{R_0}(O)$ , where $B_{R_0}(O)$ denotes the geodesic ball centered at $$O$$ with radius $$R_0$$ , within the class of functions that vanish on $\partial B_{R_0}(O)$ . In both cases, we give conditions on $$f$$ , $\varphi$ and $$κ$$ implying that the solution $$u$$ is radial and $$Ω$$ is a geodesic ball centered at $$O$$ . Our results apply in particular to the three space forms $\mathbb{R}^N$ , $\mathbb{H}^N$ and $\mathbb{S}^N$ .

3. 具有常各向异性非局部平均曲率的周期Delaunay柱面

Periodic Delaunay cylinders with constant anisotropic nonlocal mean curvature

作者: Francesc Alcover, Renzo Bruera

类别: math.AP

发布日期: 2026-02-20

arXiv ID: 2602.18219v1

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

摘要: In this article we prove existence and symmetry properties of periodic surfaces of revolution with constant anisotropic nonlocal mean curvature, generalizing a classical result of Delaunay to the anisotropic nonlocal setting. First, by studying the corresponding periodic isoperimetric problem, under natural assumptions on the kernel, we use rearrangement inequalities to extend a periodic version of the Wulff inequality to the nonlocal setting. This leads to the existence and symmetry properties of minimizers for every given volume in each period, thus generalizing the results of Cabré, Csató, and Mas to the anisotropic case. Second, under the same hypotheses on the kernel, we prove the existence of a one-parameter family of Delaunay near-cylinders in $\mathbb{R}^2$ bifurcating from a straight cylinder and having each constant anisotropic mean curvature. This extends the results of Cabré, Fall, Solà-Morales, and Weth to the anisotropic case. The stability of these near-cylinders will be studied in a forthcoming paper.

4. 关于二维空间中周期非局部周长的极小化子形状

On the shape of minimizers for the periodic nonlocal perimeter in $\mathbb{R}^2$

作者: Renzo Bruera

类别: math.AP

发布日期: 2026-02-20

arXiv ID: 2602.18215v1

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

摘要: In this paper, we study planar nonlocal Delaunay sets. That is, open sets in $\mathbb{R}^2$ with constant nonlocal mean curvature that are periodic in $$x_1$$ , and even in $$x_1$$ and in $$x_2$$ . Using bifurcation analysis and fine explicit computations, we prove that every sufficiently $C^{1,β}$ -flat nonlocal Delaunay set in $\mathbb{R}^2$ that is not a straight band is unstable with respect to volume-preserving periodic variations. Our results support the conjecture that, as in the local case, in the range of large areas, minimizers of the periodic nonlocal isoperimetric problem – also known as the nonlocal liquid drop problem with prescribed area between two parallel hyperplanes – are all straight bands.

5. 带权对数位势算子

On the weighted logarithmic potential operator

作者: T. V. Anoop, Jiya Rose Johnson

类别: math.AP

发布日期: 2026-02-20

arXiv ID: 2602.18138v1

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

摘要: For a bounded open set $Ω\subset \mathbb{R}^N$ with $N\geq 2$ , and for positive continuous functions $$w,g$$ on $\overlineΩ$ , we consider the weighted eigenvalue problem \begin{equation*} \mathcal{L}{w} u =τgu, \end{equation*} where $\mathcal{L}_{w}$ is the weighted logarithmic potential operator on $$L^2(Ω)$$ as defined below: \begin{equation*} \mathcal{L}{w} u(x)=\int_Ω\log\left(\frac{w(x)w(y)}{|x-y|}\right)u(y)dy. \end{equation*} We study the monotonicity and continuity of the largest positive eigenvalue $τ_{w,g}^+(Ω)$ with respect to $$Ω$$ , $$w$$ , and $$g$$ . We also establish that $τ_{w,g}^+(Ω)$ satisfies a reverse Faber Krahn inequality under polarization. We provide a sufficient condition for the existence of a negative eigenvalue in terms of the weighted transfinite diameter of $$Ω$$ , under the assumption that $\log w$ is superharmonic. For $Ω\subset \mathbb{R}^2$ , if $Δ\log w$ is a constant $$C$$ , we show that 0 can be an eigenvalue of $\mathcal{L}_{w}$ only when $C=\frac{2π}{|Ω|}$ . For such domains, if $\log w$ is a harmonic function on $$Ω$$ , we provide a representation formula for the eigenfunctions. Using this representation, we establish variants of the maximum principles that give some insight into the geometry of these eigenfunctions.

6. 非局部特征值问题与叠加算子

Nonlocal eigenvalue problems and superposition operators

作者: Serena Dipierro, Edoardo Proietti Lippi, Caterina Sportelli, Enrico Valdinoci

类别: math.AP

发布日期: 2026-02-20

arXiv ID: 2602.18035v1

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

摘要: We study the spectral theory of mixed local and nonlocal operators with lower-order terms in the right-hand side of the equation. This kind of problems is motivated by the analysis of superposition operators of mixed order and with the “wrong sign” of the lower-order terms with respect to the classical elliptic theory. Our results include: -convergence to classical cases when the right-hand side of the eigenvalye equations “localizes”, recovering the simplicity and sign-definiteness of eigenfunctions in the limit; -a detailed analysis of disconnected domains, showing that, unlike the classical case, any eigenfunction associated with the first eigenvalue must change sign, and that the first eigenvalue of a union of disconnected domains is strictly smaller than that of its individual components; -examples in which the first eigenvalue is either simple or non-simple in disconnected domains; -a regularity theory that underpins these results.

7. 随机朗道-栗弗席兹-布洛赫方程有限元格式的强收敛性

Strong convergence of finite element schemes for the stochastic Landau–Lifshitz–Bloch equation

作者: Agus L. Soenjaya

类别: math.NA, math.AP

发布日期: 2026-02-20

arXiv ID: 2602.18021v1

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

摘要: The dynamics of magnetisation in a bounded ferromagnet in $\mathbb{R}^d$ ( $$d=1,2$$ ) at high temperatures can be described by the stochastic Landau–Lifshitz–Bloch (sLLB) equation, which is a vector-valued quasilinear stochastic partial differential equation. In this paper, assuming adequate regularity of the initial data, we establish strong convergence in $$L^2(Ω)$$ of several semi-implicit and implicit fully discrete finite element schemes for the sLLB equation, together with explicit convergence rates. The analysis relies on localised error estimates and new exponential moment bounds for the exact solution. As a by-product, these moment bounds yield mean-square exponential stability of solutions and uniqueness of the invariant measure in one spatial dimension under a small noise assumption. We also sharpen existing convergence-in-probability results for the numerical schemes. Numerical experiments are presented to illustrate and support the theoretical findings.

8. 关于蒙日-安培型方程内部 $$C^2$$ 估计的反例

On Counterexamples to Interior $$C^2$$ Estimates for Monge-Ampère Type Equations

作者: Cheuk Yan Fung

类别: math.AP

发布日期: 2026-02-20

arXiv ID: 2602.18009v1

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

摘要: We modify Pogorelov’s classic construction to demonstrate the absence of a priori $$C^2$$ estimates for the equations $\det(D^2 u \pm Du \otimes Du) = f(x)$ in dimension $n \ge 3$ . We construct a sequence of solutions $z_\varepsilon$ with second derivatives blowing up at the origin as $\varepsilon \rightarrow 0$ , while the corresponding right-hand sides $f_\varepsilon$ admit uniform $$C^2$$ estimates. Specifically, the counterexamples are given by $z_\varepsilon(x_1, \dots, x_n) = (1+x_1^2)(1+x_2^2)(\varepsilon^2 + η^2)^{α/2},$ where $η= \sqrt{x_3^2 + \dots + x_n^2}$ and $α= 2 - \frac{2}{n}$ .

9. 结合深度学习的卡尔曼半离散凸化方法用于电阻抗断层成像

A Carleman Semi-Discrete Convexification Method Combined With Deep Learning for Electrical Impedance Tomography

作者: Michael V. Klibanov, Kirill V. Golubnichiy, Benjamin Jiang

类别: math.AP

发布日期: 2026-02-20

arXiv ID: 2602.18001v1

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

摘要: In this paper, a new semi-discrete version of the Carleman estimate-based convexification globally convergent numerical method is developed. It is used for the delivery of the starting point for the training procedure of deep learning. An important feature of the continuous version of the convexification method is that its convergence to the true solution is independent on the availability of a good first guess about this solution. A new concept of the h-strong convexity is introduced, where h is the grid step size in the semi-discrete version of the convexification method. The h -strong convexity allows to obtain an a priori accuracy estimate of the starting point for the training step of the deep learning procedure. This approach is demonstrated for a highly nonlinear problem of Electrical Impedance Tomography. Results of numerical experiments for complicated media structures demonstrate the computational feasibility of this procedure.

10. 混合局部与非局部不定方程的刘维尔定理

Liouville theorems for mixed local and nonlocal indefinite equations

作者: Pengyan Wang, Leyun Wu

类别: math.AP

发布日期: 2026-02-20

arXiv ID: 2602.17944v1

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

摘要: We investigate the qualitative properties of positive solutions to mixed local-nonlocal equations with indefinite nonlinearities, emphasizing the interaction between classical and fractional Laplacians. We first establish maximum principles and prove strict monotonicity along the $$x_1$$ -direction for mixed elliptic operators. By combining a mollified first eigenfunction with a suitable sub-solution, we derive nonexistence results for the mixed operator $$ (-Δ)^s - Δ$$ via a contradiction argument. These results are further extended to the parabolic setting, incorporating both the Marchaud-type fractional time derivative and the classical first-order derivative, revealing new qualitative features under dual nonlocality. A key aspect of our approach is a careful adaptation of the method of moving planes to the mixed local-nonlocal context. By addressing the distinct scaling behaviors of local and nonlocal terms, the method yields monotonicity and Liouville-type results without standard decay assumptions, and provides a framework potentially applicable to a broader class of mixed elliptic and parabolic problems.

11. 基于变分优化方法的介电常数与电导率函数重建（利用部分边界测量）

Variational optimization approach for reconstruction of dielectric permittivity and conductivity functions using partial boundary measurements

作者: Eric Lindström, Larisa Beilina

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

发布日期: 2026-02-19

arXiv ID: 2602.17819v1

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

摘要: We present a variational optimization approach for the solution of a coefficient inverse problem of simultaneous reconstruction of the dielectric permittivity and conductivity functions in time-dependent Maxwell’s system using limited boundary observations of the electric field. The variational optimization approach is based on constructing a weak form of a Lagrangian which allows to use finite element based reconstruction algorithms. The optimality conditions for the Lagrangian and stability estimate for the adjoint problem are derived, as well as Frechét differentiability of it and of the regularized Tikhonov functional are also presented. Two- and three-dimensional numerical studies confirm our theoretical investigations.

12. 距离函数、曲率与拓扑学

Distance Functions, Curvature and Topology

作者: Carlo Mantegazza, Francesca Oronzio

类别: math.DG, math.AP

发布日期: 2026-02-19

arXiv ID: 2602.17629v1

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

摘要: We discuss some properties of the distance functions on Riemannian manifolds and we relate their behavior to the geometry of the manifolds. This leads to alternative proofs of some “classical” theorems connecting curvature and topology.

13. 关于具有完全单调核的随机Volterra积分方程与具有无界控制及扩散算子的希尔伯特空间随机微分方程的最优控制

Optimal control of stochastic Volterra integral equations with completely monotone kernels and stochastic differential equations on Hilbert spaces with unbounded control and diffusion operators

作者: Gabriele Bolli, Filippo de Feo

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

发布日期: 2026-02-19

arXiv ID: 2602.17578v1

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

摘要: The dynamic programming approach is one of the most powerful ones in optimal control. However, when dealing with optimal control problems of stochastic Volterra integral equations (SVIEs) with completely monotone kernels, deep mathematical difficulties arise and it is still not understood. These very classical problems have applications in most fields and have now become even more popular due to their applications in mathematical finance under rough volatility. In this article, we consider a class of optimal control problems of SVIEs with completely monotone kernels. Via a recent Markovian lift \cite{FGW2024}, the problem can be reformulated as an optimal control problem of stochastic differential equations (SDEs) on suitable Hilbert spaces, which due to the roughness of the kernel, presents a generator of an analytic semigroup and unbounded control and diffusion operators. This analysis leads us to study a general class of optimal control problems of abstract SDEs on Hilbert spaces with unbounded control and diffusion operators. This class includes optimal control problems of SVIEs with completely monotone kernels, but it is also motivated by other models. We analyze the regularity of the associated Ornstein-Uhlenbeck transition semigroup. We prove that the semigroup exhibits a new smoothing property in control directions through a general observation operator $$Γ$$ , which we call $$Γ$$ -smoothing. This allows us to establish existence and uniqueness of mild solutions of the Hamilton-Jacobi-Bellman equation, establish a verification theorem, and construct optimal feedback controls. We apply these results to optimal control problems of SVIEs with completely monotone kernels. To the best of our knowledge these are the first results of this kind for this abstract class of infinite dimensional problems and for the optimal control of SVIEs with completely monotone kernels.

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

1. Stability of the Shape for Circular Vortex Filaments under Non-Symmetric Perturbations#

详细分析#

研究对象和背景#

主要定理或主要结果#

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

与之前工作的比较#

2. 毛细波与一类拟线性偏微分方程的不变流形#

详细分析#

1. 研究对象和背景#

物理系统与数学表述#

挑战与动机#

水波系统的函数空间与边值条件#

2. 主要定理#

一般框架下的定理（Theorem 2.4 和 Theorem 2.9）#

水波系统的定理（Theorem 4.2）#

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

a) 核心分析学工具#

b) 解决技术性难题的主要技巧#

核心工具总结#

4. 与之前工作的比较#

与特定工作的对比#

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

1. 衰减密度不确定性原理在薛定谔方程可观测性中的应用#

摘要翻译#

2. 具有奇异各向异性核的卷积算子的非局部到局部$L^p$收敛性#

摘要翻译#

3. 分数阶拉普拉斯算子的弱奇性表述#

摘要翻译#

不相关论文（基本信息）#

1. 分数阶亥姆霍兹方程的极限吸收原理与辐射条件#

2. 旋转不变泊松方程在模型流形上的超定问题#

3. 具有常各向异性非局部平均曲率的周期Delaunay柱面#

4. 关于二维空间中周期非局部周长的极小化子形状#

5. 带权对数位势算子#

6. 非局部特征值问题与叠加算子#

7. 随机朗道-栗弗席兹-布洛赫方程有限元格式的强收敛性#

8. 关于蒙日-安培型方程内部 $C^2$ 估计的反例#

9. 结合深度学习的卡尔曼半离散凸化方法用于电阻抗断层成像#

10. 混合局部与非局部不定方程的刘维尔定理#

11. 基于变分优化方法的介电常数与电导率函数重建（利用部分边界测量）#

12. 距离函数、曲率与拓扑学#

13. 关于具有完全单调核的随机Volterra积分方程与具有无界控制及扩散算子的希尔伯特空间随机微分方程的最优控制#

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

1. Stability of the Shape for Circular Vortex Filaments under Non-Symmetric Perturbations

详细分析

研究对象和背景

主要定理或主要结果

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

与之前工作的比较

2. 毛细波与一类拟线性偏微分方程的不变流形

详细分析

1. 研究对象和背景

物理系统与数学表述

挑战与动机

水波系统的函数空间与边值条件

2. 主要定理

一般框架下的定理（Theorem 2.4 和 Theorem 2.9）

水波系统的定理（Theorem 4.2）

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

a) 核心分析学工具

b) 解决技术性难题的主要技巧

核心工具总结

4. 与之前工作的比较

与特定工作的对比

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

1. 衰减密度不确定性原理在薛定谔方程可观测性中的应用

摘要翻译

2. 具有奇异各向异性核的卷积算子的非局部到局部 $$L^p$$ 收敛性

摘要翻译

3. 分数阶拉普拉斯算子的弱奇性表述

摘要翻译

不相关论文（基本信息）

1. 分数阶亥姆霍兹方程的极限吸收原理与辐射条件

2. 旋转不变泊松方程在模型流形上的超定问题

3. 具有常各向异性非局部平均曲率的周期Delaunay柱面

4. 关于二维空间中周期非局部周长的极小化子形状

5. 带权对数位势算子

6. 非局部特征值问题与叠加算子

7. 随机朗道-栗弗席兹-布洛赫方程有限元格式的强收敛性

8. 关于蒙日-安培型方程内部 $$C^2$$ 估计的反例

9. 结合深度学习的卡尔曼半离散凸化方法用于电阻抗断层成像

10. 混合局部与非局部不定方程的刘维尔定理

11. 基于变分优化方法的介电常数与电导率函数重建（利用部分边界测量）

12. 距离函数、曲率与拓扑学

13. 关于具有完全单调核的随机Volterra积分方程与具有无界控制及扩散算子的希尔伯特空间随机微分方程的最优控制