Chapter 10: Breathers in the Fractional Frenkel-Kontorova Model

1. 作者

J. Catarecha（通信作者），西班牙塞维利亚大学（Universidad de Sevilla）应用物理系 + 数学研究所（IMUS），Email: catarechajorge@gmail.com
J. Cuevas-Maraver，同上机构，Email: jcuevas@us.es
P. G. Kevrekidis（系列书编者），美国马萨诸塞大学阿默斯特分校（UMass Amherst）数学与统计系，Email: kevrekid@umass.edu

Kevrekidis 是该书主编，Cuevas-Maraver 是 sine-Gordon 模型专家，Catarecha 是博士生。三人合作聚焦于分数阶 Frenkel-Kontorova (FK) 模型——经典 FK 模型的长程相互作用推广。工作由 NSF PHY-2110030、DMS-2204702 和西班牙 PID2020-112620GB-I00、PID2022-143120OB-I00 资助。

2. 内容概述

本章是全书最"具体"的章节——把分数阶离散思想应用到Frenkel-Kontorova 模型（即离散 sine-Gordon），研究离散呼吸子（breather）——时间周期、空间局域的离散行波。3 大主题：(1) 模型与色散（§2）——分数阶 FK 方程 $\ddot u_n + \sin u_n + C(-\Delta)^\alpha u_n = 0$，$K_\alpha(m) \sim |m|^{-1-2\alpha}$ 长程耦合，色散关系 $\Omega^2(q) = 1 + 2^\alpha C[1 - \cos q]^\alpha$；(2) 数值结果（§3）——单格点（onsite）、同相双格点（in-phase inter-site）、反相双格点（out-of-phase）三种呼吸子的代数尾衰减 $|n|^{-2}$（onsite/in-phase）和 $|n|^{-3}$（out-of-phase），"双衰减"——近中心指数、远尾代数；Floquet 稳定性分析，4 种分岔（$C_1$ Hopf、$C_2$ 稳定性交换、$C_3$ phantom、$C_4$ 重稳定化）；(3) 动力学——多稳态、nanopteron、kick 激发移动呼吸子。

3. 核心方程与概念

3.1 分数阶 FK 模型

控制方程（eq. 1）： $$\ddot u_n + \sin u_n + C(-\Delta)^\alpha u_n = 0$$ 分数阶离散 Laplacian（eq. 2-3）： $$(-\Delta)^\alpha u_n = \sum_{m \ne n} K_\alpha(n-m)(u_n - u_m), \quad K_\alpha(m) = \frac{4^\alpha \Gamma(\alpha+1/2)}{\sqrt\pi|\Gamma(-\alpha)|\Gamma(|m|+1+\alpha)}$$ - $\alpha = 1$：$K_{\alpha=1}(m) = \delta_{|m|,1}$，退化为标准离散 Laplacian； - $\alpha \to 0$：$K_\alpha(m) \to 1/|m|$，全格点长程耦合。

等价形式（eq. 4）： $$\ddot u_n + \sin u_n + C\sum_{m>0} K_\alpha(m)(2u_n - u_{n+m} + u_{n-m}) = 0$$

Hamiltonian（eq. 5）： $$H = \sum_n E_n, \quad E_n = \frac{\dot u_n^2}{2} + (1 - \cos u_n) + \frac{C}{4}\sum_{m>0} K_\alpha(m)[(u_n - u_{n+m})^2 + (u_n - u_{n-m})^2]$$

色散关系（声子谱，eq. 10-11）： $$\Omega^2(q) = 1 + 4C\sum_{m \ge 1} K_\alpha(m)\sin^2(mq/2) = 1 + 2^\alpha C[1 - \cos q]^\alpha$$ 声子带 $\Omega_{ph}^2 \in [1, 1 + 2^\alpha C \cdot 4^\alpha = 1 + 4^\alpha C]$。

关键性质： - $\alpha \to 0$ → $\Omega(q) \to 1$（平带，无色散）； - $0.5 < \alpha < 1$：$\Omega^2$ 在 $q = 0$ 处三阶导数不连续； - $\alpha \le 0.5$：二阶导数不连续。

3.2 离散呼吸子的代数衰减

时间周期解（eq. 6，时间可逆的余弦 Fourier 截断）： $$u_n(t) = z_{0,n} + 2\sum_{k=1}^{k_m} z_{k,n}\cos(k\omega t)$$ 求解 $z_{k,n}$ 的非线性代数方程用 Newton-Raphson（[X1]）。

3 种基本呼吸子： 1. 单格点（onsite）——经典单点呼吸子 $z_{0,0} \ne 0$，其他 $z_{k,0} = 0$； 2. 同相双格点（in-phase inter-site）——相邻两点同相振动； 3. 反相双格点（out-of-phase, twisted）——相邻两点反相振动。

双衰减（关键发现，图 3）： - 近中心：指数衰减 $\sim e^{-\kappa n}$； - 远尾（$n > n_c$）：代数衰减 $\sim |n|^{-p}$，其中 $n_c$ 随 $\alpha$ 减小而减小。

衰减指数（[X2] Flach 1998 的工作预测 $K(n) = |n|^{-s}$，$s = 1+2\alpha$）： - 单格点 + 同相双格点：$p = s + 1 = 2\alpha + 2$，对 $\alpha = 0.5$ → $p = 2$（验证）； - 反相双格点：$p = s + 1 + 1 = 2\alpha + 3$？实际 $p = 3$（与 $\alpha$ 无关的反常现象）——两个 onsite 反相叠加导致 $\sim |n|^{-s} - |n+1|^{-s} \approx |n|^{-(s+1)}$，加上对 onsite 的 $|n|^{-s}$ 修正，总为 $|n|^{-3}$。

3.3 Floquet 稳定性

Floquet 算子（eq. 8）： $$\Xi(T) = F\Xi(0), \quad \Xi(t) = [\xi_n(t), \dot\xi_n(t)]$$ Floquet 乘子 $\lambda = e^{i\theta}$： - 实本征值：$\lambda, 1/\lambda$ 对； - 复本征值：$\lambda, \lambda^*, 1/\lambda, 1/\lambda^*$ 四重态； - Hamilton 系统 + 周期解 → $\lambda = 1$（$\theta = 0$）总是存在。

Krein 签名（eq. 9）： $$\kappa(\lambda) = \text{sgn}\sum_n \text{Im}\xi_n^{(\lambda)}(0)\dot\xi_n^{(\lambda)}(0)$$ 决定双稳定性——当 $\kappa$ 反号的乘子对碰撞时会出现 Hopf 型失稳。

3.4 4 类分岔

$C_1$（Hopf 型，弱耦合失稳）： $$C_1 = \frac{9\omega^2 - 4}{4^{1+\alpha}}$$ 对 $\alpha = 0.95, \omega = 0.8$：$C_1 = 0.12$；$\alpha = 0.5$：$C_1 = 0.23$。声子弧碰撞引起——声子带边界 $\theta_{ph} = 2\pi\Omega_{ph}/\omega$ 抵达 $\pi$。

$C_2$（稳定性交换，指数型）： $\alpha = 0.95$：$C_2 = 0.83$；$\alpha = 0.5$：$C_2 = 1.10$。反对称内模变实——单格点 → 同相双格点 + 非对称双格点（"stability exchange bifurcation"）。

$C_3$（phantom 呼吸子）： $$C_3 = \frac{9\omega^2 - 1}{4^\alpha}$$ $\alpha = 0.95$：$C_3 = 1.28$；$\alpha = 0.5$：$C_3 = 2.39$。呼吸子频率的奇数倍与声子带共振。

$C_4$（重稳定化）：$C_2$ 失稳的乘子"回入"单位圆。$\alpha = 0.5$ 时 $C_4 = 1.81$。

4 类分岔的总体模式（图 4-5）： - $C_1$ 几乎与 $\alpha$ 无关（声子带边缘仅微弱受 $\alpha$ 影响）； - $C_2, C_3$ 都随 $\alpha$ 减小而增大（声子带变窄 → 临界值右移）； - 失稳增长率随 $\alpha$ 减小而减小——长程耦合"软化"了失稳。

3.5 Floquet 乘子的解析近似

小耦合下乘子（eq. 15）： $$\theta = \pm\frac{2\pi i}{\omega}\sqrt{2C K_\alpha(1)} \frac{\omega^2\sqrt{3(4\omega^2 - 1)} - (2\omega^2 - 1)\sqrt{4\omega^2 - 1}}{16\omega^4 - 7\omega^2 + 1}$$ 对 onsite 呼吸子：$\text{Im}(\theta) = \ln|\lambda|$ → 失稳。对反相双格点：$\theta$ 实数 → 谱稳定。对同相双格点（in-phase inter-site）：$\text{Im}(\theta)$ 与 $C$ 线性——指数失稳。对反相双格点（out-of-phase）：$\theta$ 实数 → 谱稳定。 反相双格点不能继续到连续极限——纯离散态。

3.6 能量-频率关系与多稳态

能量-频率单调性判据（[X3] 2017 定理）： $$H'(\omega) > 0 \Rightarrow \text{局部失稳} \Leftrightarrow \text{呼吸子在该频段不稳定}$$ 对 $\alpha = 0.5$, $C = 0.5$ 时，$\omega \in [0.909, 0.962]$ 区间内 $H(\omega)$ 出现"局部极大" → 失稳区间（图 8-9）。这一现象在 $\alpha = 1$（标准 FK）下不存在——分数阶"长程相互作用"诱发了多稳态。

3.7 动力学：移动呼吸子与 Nanopteron

Kick 激发（沿不稳定 Floquet 乘子的本征方向扰动）： - 在 $[C_2, C_4]$ 区间：单格点呼吸子因反对称内模失稳 → 加速度 → 移动呼吸子； - $\alpha = 0.95$：移动呼吸子稳定传播（图 12 左）； - $\alpha = 0.5$：移动呼吸子减速 + 辐射耗散（图 12 右）——长程耦合使尾巴振幅大 → kick 后辐射损失更多。

Nanopteron——失稳 + 静息 → 双频束缚结构（呼吸子上叠加高频振动）——长寿命且具有"呼吸子叠加在平台上"的形状。

4. 关键结论

分数阶 FK 模型保持经典 FK 的大部分特征——离散呼吸子、$C = 0$ 极限的延续、双稳态等，但尾衰减从指数变为代数。
"双衰减"现象（近中心指数、远尾代数）——$n_c$ 随 $\alpha$ 减小而减小；这是分数阶与最近邻耦合的本质差异。
3 类呼吸子的不同尾衰减指数：
单格点/同相双格点：$|n|^{-2}$（$\alpha = 0.5$）；
反相双格点：$|n|^{-3}$（与 $\alpha$ 无关的反常现象）。
4 类 Floquet 分岔——$C_1$ Hopf、$C_2$ 稳定性交换、$C_3$ phantom、$C_4$ 重稳定化——临界值 $C_2, C_3$ 随 $\alpha$ 减小而右移（长程耦合"软化"失稳）。
能量-频率单调性判据（[X3] 2017）在分数阶情形仍然有效——但多稳态区间仅在 $\alpha < 1$（分数阶）下出现。
Nanopteron + 移动呼吸子是分数阶 FK 模型的特征动力学态——长程耦合使失稳后的辐射更多，移动呼吸子难以稳定。
$\alpha = 0.5$ 临界值——文章主要研究 $0.5 \le \alpha \le 1$；$\alpha < 0.5$ 时代数衰减太慢，数值上需要极大晶格，Ch6 的多区域方法可作处理。

5. 挑战和开放性问题

$\alpha < 0.5$ 的强分数阶情形——代数衰减 $|n|^{-1-2\alpha}$ 太慢，需要 Ch6 的多区域方法处理。
2D 分数阶 FK 模型——平面呼吸子、带涡度的呼吸子（$2\pi$ 相位缠绕，[X4] 已有相关研究）。
移动呼吸子的严格存在性——$C \in [C_2, C_4]$ 时有移动呼吸子但稳定性如何？
Nanopteron 的长寿性——长程耦合下，nanopteron 是稳定还是缓慢衰减？
能量-频率单调性与稳定性的更深联系——多稳态区间的宽度随 $\alpha, C$ 如何变化？
与连续分数阶 FK 模型的衔接——$C \to \infty$ 时的连续极限行为？分数阶 sine-Gordon 方程是否也有连续呼吸子？
反相多格点呼吸子的详细分析——本文只给出 $\alpha = 0.5$ 时的稳定态，3-site, 4-site 等待研究。
时间-不可逆呼吸子（multibreather 破缺时间对称性）——本文方法基于余弦 Fourier 截断，限制为可逆态；非可逆态需新方法 [X5]。
离散分数阶 sine-Gordon 的实验实现——长程耦合 torsion pendulum 实验 [X6] 是可行的。

6. 个人反思与批判性分析

优点： 1. "统一现象学"——在 1 个方程框架内统一处理单格点/同相/反相呼吸子、4 类分岔、能量-频率单调性——读者获得"全景图"。 2. "双衰减"是真正的发现——单格点/同相双格点的 $|n|^{-2}$ vs 反相双格点的 $|n|^{-3}$ 反常差异，给出清晰的解析论证。 3. Floquet 分析严格执行——计算 monodromy 矩阵、特征值、Krein 签名——与 Ch3 亚扩散系统的谱分析一脉相承。 4. 能量-频率判据（[X3] 2017）的分数阶推广——把"分数阶长程耦合"与"稳定性的能量判据"联系起来。 5. "标准 Laplacian vs 分数阶"对照——$\alpha = 0.95$（接近 1）vs $\alpha = 0.5$（远离 1）——明确展示分数阶效应何时"出现"。 6. Nanopteron + 移动呼吸子是"活的"动力学示例——理论（Floquet）+ 计算（直接模拟）双管齐下。

缺点： 1. $\alpha$ 范围限制为 $[0.5, 1]$——题目"fractional FK"但实际只覆盖"弱分数阶"情形。 2. 没有反相多格点呼吸子的详细分析（只给 2-site 结果）。 3. 与 Ch3 亚扩散的 A.lin（"分数阶使稳定"） 是平行现象——但本章未引用。 4. 数值方法细节不足——只说 Newton-Raphson，未给出实际精度测试或效率比较。 5. 无连续极限衔接——$C \to \infty$ 时的行为（连续分数阶 sine-Gordon 是否有连续呼吸子？）未讨论。 6. "双衰减"的解析论证只到"$\sim |n|^{-s} - |n+1|^{-s}$"——更精细的"内层指数 + 外层代数 + 过渡区 $n_c$"的解析公式未给出。 7. 与 Ch8 的分数阶可积方程无任何交叉引用——FK 模型是否"近似可积"？是否存在分数阶 sine-Gordon 的反散射变换？

对本人研究（血管生物力学 / SMC 介导的血管重塑）的关联： - 分数阶 FK 模型与血管壁的"层状晶格"类比——内弹力板 / 中膜 / 外膜的弹性 + 主动收缩构成"离散弹性链"，分数阶描述长程应力传递。 - "双衰减"与血管壁中"SMC 应力的非局部传播"类比——近 SMC 局部应力指数衰减、远距离代数衰减——这与 Goriely 生长-重塑理论中的"非局部效应"同源。 - 稳定性交换分岔与血管壁的"相变"类比——从"稳态"到"不稳定"的相变对应病理重塑（如动脉瘤前兆）。 - 移动呼吸子与脉搏波类比——脉搏波是"波"在血管壁的传播，分数阶 FK 的移动呼吸子是理论模型。 - 多稳态（nanopteron）与血管壁的"病态稳态"类比——动脉瘤形成前的血管壁可能是"多稳态"——一个稳定的、但脆弱的态。

7. 重要参考文献（按出现顺序编号）

[X1] J.L. Marín, S. Aubry. Nonlinearity 9, 1501 (1996). —— Newton-Raphson 方法用于离散呼吸子。
[X2] S. Flach. Phys. Rev. E 58, R4116 (1998). —— 长程耦合晶格的呼吸子代数衰减理论。
[X3] P.G. Kevrekidis, J. Cuevas-Maraver, D.E. Pelinovsky. Phys. Rev. Lett. 119, 094101 (2016). —— 能量-频率单调性稳定性判据。
[X4] C. Mejía-Cortés, M.I. Molina. Opt. Lett. 46, 2256 (2021). —— 带涡度的离散呼吸子。
[X5] V. Koukouloyannis, P.G. Kevrekidis, J. Cuevas, V. Rothos. Phys. D 242, 16 (2013). —— 时间-不可逆多呼吸子。
[X6] J. Cuevas, L.Q. English, P.G. Kevrekidis, M. Anderson. Phys. Rev. Lett. 102, 224101 (2009). —— 长程耦合 torsion pendulum 实验。
[X7] S. Longhi. Opt. Lett. 40, 1117 (2015). —— 4f 腔分数阶衍射（Ch1, Ch8, Ch9 也引）。
[X8] S. Liu, Y.W. Zhang, B.A. Malomed, E. Karimi. Nat. Commun. 14, 222 (2023). —— 首次实验分数阶 GVD（Ch1, Ch4, Ch8, Ch9 也引）。
[X9] M.I. Molina. Phys. Lett. A 384, 126180 (2020). —— 分数阶离散 NLS（Ch9 也引）。
[X10] M.I. Molina. Phys. Lett. A 384, 126835 (2020). —— 2D 分数阶离散 NLS（Ch9 也引）。
[X11] O.M. Braun, Yu.S. Kivshar. Phys. Rep. 306, 1 (1998). —— FK 模型综述。
[X12] O.M. Braun, Yu.S. Kivshar. The Frenkel-Kontorova Model (Springer, 2004). —— FK 模型专著。
[X13] A.C. Scott. Am. J. Phys. 37, 52 (1969). —— 扭转摆链实验。
[X14] L.V. Yakushevich. Nonlinear Physics of DNA (Wiley, 1998). —— DNA 与 sine-Gordon 模型。
[X15] S. Flach, A.V. Gorbach. Phys. Rep. 467, 1 (2008). —— 离散呼吸子综述。
[X16] S. Aubry. Phys. D 103, 201 (1997). —— 离散呼吸子 Aubry 理论。
[X17] O. Ciaurri, L. Roncal, P.R. Stinga, J.L. Torrea, J.L. Varona. Adv. Math. 330, 688 (2018). —— 分数阶离散 Laplacian 定义（Ch9 也引）。
[X18] A.A. Kilbas, H.M. Srivastava, J.J. Trujillo. —— 分数阶微积分。
[X19] A.M. Morgante, M. Johansson, S. Aubry, G.A. Kopidakis. J. Phys. A 35, 4999 (2002). —— Phantom 呼吸子。
[X20] A.M. Turing. Philos. Trans. Royal Soc. London, Ser. B 237, 37 (1952). —— Turing 不稳定（Ch2, Ch3, Ch5, Ch8, Ch9 也引）。

本章由主 agent 亲自从 PDF 提取并撰写，无子任务委托。
撰写日期：2026-06-01