Chapter 9: Fractional Discrete Linear and Nonlinear Models

1. 作者

Mario I. Molina（独立作者），智利大学（University of Chile）理学院物理系，圣地亚哥 Nunoa，Email: mmolina@uchile.cl。

Molina 是离散非线性晶格 + 分数阶扩展的专家，长期在 Chile 当地工作。他的工作连接 Ch1-Ch8 的抽象理论与"具体离散晶格实现"（紧束缚模型、传输线、Anderson 模型）。本章与 Ch4 (PT 对称 FNLS) 有部分重叠（PT 对称 + 离散），但视角更"实验"——以数值的、可观测的现象为重心。

工作由 Fondecyt Grant 1200120 资助（智利国家科学发展基金）。

2. 内容概述

本章是全书最"应用驱动"的章节之一。5 大主题：(1) 1D 分数阶 DNLS 方程 + PT 对称（§2.1）——紧束缚模型，核 $K_s(m) \sim |m|^{-1-2s}$，色散带宽 $4V$ (s=1) → $V$ (s=0)，PT 对称窗口；(2) 2D 分数阶 DNLS（§2.2）——$K_s(m)$ 用修正 Bessel 函数 $I_1, I_2$ 表达；(3) 分数阶电传输线（§3）——双电感 LC 电路，RWA 下非线性方程，PT 对称与分数阶的相互作用；(4) 分数阶非线性杂质（§4）——单杂质处的分数阶离散方程，Green 函数 + 束缚态 + 透射率 + 自陷阈值随 $s$ 减小而降低；(5) 分数阶 Anderson 局域化（§5）——分数阶对无序晶格的影响，意外的"反常迁移率"现象（在 $0.15 < s < 0.4$ 区间）。

关键发现 1（§2）：分数阶减小 → 带宽减小 → 模式兼并度增加；PT 对称窗口出现在 $s$ 接近 1 时。 关键发现 2（§4）：分数阶减小 → 自陷阈值降低（弱非线性即可"捕获"激发）。 关键发现 3（§5）：分数阶 Anderson 模型中，MSD 在 $0.15 < s < 0.4$ 出现"驼峰"——分数阶反而增强迁移率，与经典无序导致的"全局抑制"相反。

3. 核心方程与概念

3.1 1D 分数阶离散非线性 Schrödinger 方程（fDNLS）

标准 DNLS 方程（无 PT 对称）： $$i\frac{dC_n}{dt} + V(C_{n+1} + C_{n-1}) + \chi|C_n|^2 C_n = 0$$

分数阶离散 Laplacian（关键，eq. 7-8）： $$(\Delta_n)^s \phi_n = \sum_{m \ne n} K_s(n - m)(\phi_m - \phi_n), \quad K_s(m) = \frac{4^s \Gamma(s + 1/2)}{\sqrt\pi|\Gamma(-s)|\Gamma(|m| + 1 + s)}$$ 渐近行为 $K_s(m) \approx |m|^{-(1+2s)}$ ——长程耦合。 - $s = 1$：退化为 $(\Delta_n)\phi_n = \phi_{n+1} - 2\phi_n + \phi_{n-1}$； - $s \to 0$：$K_s(m) \to 1/|m|$，所有格点全局耦合（"$\alpha$ 稳定"）。

含 PT 对称 + 分数阶的 fDNLS（eq. 6, 9，odd/even 交替 $\gamma_n = \pm\gamma$）： $$i\frac{da_n}{dt} + i\gamma a_n + 2V b_n + V\sum_m K_s(n-m)(b_m - b_n) + \chi|a_n|^2 a_n = 0$$ $$i\frac{db_n}{dt} - i\gamma b_n + 2V a_n + V\sum_m K_s(n-m)(a_m - a_n) + \chi|b_n|^2 b_n = 0$$

色散关系（线性，$\chi = 0$，$a_n, b_n \sim e^{i(kx - \lambda t)}$，eq. 13-14）： $$\lambda_0(k) = 2V - \frac{4^s V \Gamma(s + 1/2)}{\sqrt\pi\Gamma(1+s)}[2(1 - s) e^{-ik} F(1, 1-s; 2+s; e^{-ik}) + c.c.]$$ 用正则化超几何函数 $F(a,b;c;z) = {}_2F_1(a,b;c;z)/\Gamma(c)$ 表达。

关键观察（图 1）：$s = 1^-$ 时带宽 $4V$（标准），$s \to 0$ 时带宽 $V$——分数阶减小 → 模式兼并度增加。

PT 对称窗口（图 3，eq. 14）：$\gamma$ 固定时，$s$ 接近 1 → 谱实部，$s$ 远小于 1 → 复本征值出现 → PT 对称破缺。

均方位移（MSD，eq. 16-18）： $$\langle n^2 \rangle = 2\sum_{m=1}^\infty (m K_s(m))^2 (Vt)^2 = 24^{s-1} s \Gamma(s+1/2)^2 \left[\frac{1}{\Gamma(1+s)^2} + \frac{8s(s-1)^2}{\Gamma(3+s)^2} {}_3F_2(\{3, 2-s, 2-s\}; \{3+s, 3+s\}; 1)\right](Vt)^2$$ 总是 ballistic 输运（$t^2$ 增长），但有效速度随 $s$ 减小而减小。

3.2 2D 分数阶 DNLS

2D 分数阶离散 Laplacian（[X1] Ciaurri-Roncal 2020，eq. 25-26）： $$(\Delta_n)^s C_n = \sum_{m \ne n} (C_m - C_n) K_s(n - m), \quad K_s(m) = \frac{1}{|\Gamma(-s)|}\int_0^\infty e^{-4t} I_{m_1}(2t) I_{m_2}(2t) t^{-1-s}dt$$ 用修正 Bessel 函数 $I_m(x)$ 表达。

色散关系（eq. 27）： $$\lambda(\mathbf{k}) = 4 + \sum_m (e^{i\mathbf{k}\cdot\mathbf{m}} - 1) K_s(\mathbf{m})$$ 2D 带宽:沿三个 $\mathbf{k}$ 方向都 $s$ 减小 → 带宽减小（图 6）。

MSD（eq. 28-29）： $$\langle n^2 \rangle = \frac{1}{4\pi^2}\int_{FBZ} |\nabla\lambda(\mathbf{k})|^2 d^2k \cdot t^2$$ 弹道速度 $\propto s$：$s = 0$ 时速度 $\to 0$（"冻结"），$s = 1$ 时速度 = 2。

3.3 分数阶电传输线

双电感 LC 电路（图 7）：第 $n$ 单元包含 $L_1$、$L_2$ 串联电感和 $C_n = C_0(1 + \chi U_n^2)$ 非线性电容（Kerr 介质）。

基尔霍夫定律（eq. 30）： $$\frac{d^2 Q_n}{dt^2} = \frac{1}{L_1}(U_{n+1} - 2U_n + U_{n-1}) - \frac{1}{L_2}U_n$$

RWA（旋转波近似）下的稳态方程（eq. 31）： $$-\Omega^2 V_n + \frac{3}{4}\chi V_n^3 + \sum_m K_s(n-m)(V_m - V_n) + \omega^2 V_n = 0$$

PT 对称版本（eq. 33-34，$R_n = (-1)^n R$ 交替电阻）： $$(-\Omega^2 + i\gamma\Omega + 2\omega^2)a_n - \omega^2\left(2b_n + \sum_m K_s(n-m)(b_m - b_n)\right) = 0$$ $$(-\Omega^2 - i\gamma\Omega + 2\omega^2)b_n - \omega^2\left(2a_n + \sum_m K_s(n-m)(a_m - a_n)\right) = 0$$

PT 对称 + 分数阶的色散（eq. 35-36）： $$\Omega_k^\pm = \frac{1}{2}\left(-\gamma + 4\omega^2 \pm \sqrt{\gamma^4 + 4H(k)^2 - 8\gamma^2\omega^2}\right), \quad H(k) = 2\omega^2 - 4\omega^2\sum_m K_s(m)\sin^2(mk/2)$$

稳定性（图 10）：PT 对称窗口仅在 $s$ 接近 1 时存在；$s$ 减小 → 谱更易变复。

3.4 分数阶非线性杂质

单非线性杂质（$n = 0$ 处，eq. 37）： $$i\frac{dC_n}{dt} + V(\Delta_n)^s C_n + \delta_{n,0}\chi|C_n|^\beta C_n = 0$$

Hamiltonian 框架（eq. 39-43）： $$H = H_0 + H_1, \quad H_0 = \sum_n \epsilon_n |n\rangle\langle n| + \sum_{n,m} V_{nm}|m\rangle\langle n|, \quad H_1 = \chi|\phi_0|^\beta |0\rangle\langle 0|$$ 其中 $\epsilon_n = 2V - V\sum_{m\ne n} K_s(n-m)$，$V_{nm} = K_s(n-m)$。

Green 函数方法（eq. 44-50）： $$G^{(0)}_{nm}(z) = \frac{1}{2\pi}\int_{-\pi}^\pi \frac{e^{ik(n-m)}}{z - \lambda(k)}dk$$ 束缚态能量 $z_b$ 满足 $1 = \chi|\phi_0|^\beta G^{(0)}_{00}(z_b)$；模在 $n$ 处的振幅 $|\phi_n|^2 = -G^{(0)2}_{n0}(z_b)/G^{(0)'}_{00}(z_b)$。

透射率（eq. 51-53）： $$t(z) = \frac{1}{1 + \chi^2 |G^+_{00}(z)|^2 t(z)^\beta}, \quad \beta = 2 \text{ 时为立方方程} bt^3 + t - 1 = 0$$ 分数阶效应：$s$ 减小 → 整体透射率增大（反直觉——分数阶"减少限制"）。

自陷阈值（图 12）：分数阶减小 → 自陷阈值降低。$s \to 0$ 时即使零非线性也有"完全捕获"——长程耦合导致任何局域激发都无法逃逸。

3.5 分数阶 Anderson 局域化

分数阶 Anderson 模型（eq. 54，$\epsilon_n \in [-W, W]$ 均匀分布）： $$-\lambda C_n + 2VC_n + \epsilon_n C_n + V\sum_m K_s(m-n)(C_m - C_n) = 0$$

态密度 + 局域化（图 13-16）： - $W = 0$（无序为 0）时：$s \to 0$ → 带宽 $\to 0$（平带），$s \approx 0.415$ 处两个最低本征值从主带脱离形成离散"舌头"； - $W$ 增大：舌头逐渐模糊、并入主带； - $s$ 固定，$W$ 增大 → 谱展宽；$W$ 固定，$s$ 减小 → 谱收窄——两者可相互抵消。

参与比（Participation Ratio, PR, eq. 56）： $$R = \frac{(\sum_n |\phi_n|^2)^2}{\sum_n |\phi_n|^4}, \quad R \to N \text{ 延展态}, \quad R \to 1 \text{ 完全局域}$$

关键观察（图 15）：$s$ 固定，$W$ 增大 → $\langle\langle R\rangle\rangle$ 减小（更局域化）；$W$ 固定，$s$ 减小 → $\langle\langle R\rangle\rangle$ 减小——分数阶也增强局域化。

MSD 反常"驼峰"（图 16，$W = 1$）： - $1 > s \gtrsim 0.4$：MSD 趋于饱和（局域化），且 $s$ 减小 → MSD 减小（更局域）； - $0.15 < s < 0.4$：驼峰——MSD 不饱和，$s$ 减小 → MSD 增大（反常增强迁移率）； - $0 < s < 0.15$：MSD 急剧下降到 0（强局域化，"冻结"）。

3 个 regime 的物理解释： - Regime 1（$s \gtrsim 0.4$）：分数阶使"跳跃长度"变长但仍能局部绕行——经典局域化； - Regime 2（$0.15 < s < 0.4$）：长程跳跃（$K_s(m) \to |m|^{-(1+2s)}$） 允许绕过无序"绕过路径"——存在"extended states population"，与 [X2] 长期跳跃 Anderson 模型一致； - Regime 3（$s < 0.15$）：跳跃太长"反而冻结"——波包被"钉在原地"。

4. 关键结论

分数阶减小 → 带宽减小 → 模式兼并度增加——这是 1D 和 2D DNLS 共有的现象。
PT 对称 + 分数阶的相互作用复杂——分数阶使 PT 对称更容易破缺（$\gamma$ 临界值减小），但小窗口仍存在（$s \approx 1$）。
分数阶减小 → 自陷阈值降低——长程耦合使激发"难以逃逸"非线性势阱。
分数阶 Anderson 局域化的 3 个 regime——尤其是"驼峰"现象 $0.15 < s < 0.4$ 揭示分数阶与无序的非单调相互作用：分数阶可以抑制也可以增强局域化。
MSD 总是 ballistic（$t^2$）（无 PT 情形）——分数阶只改变"有效速度"不改变 scaling 律。
分数阶离散 Laplacian（$K_s(m) \approx |m|^{-(1+2s)}$）是"长程耦合"——既可视为"分数阶算子"也可视为"带长程跳跃的紧束缚模型"，两种解释等价。

5. 挑战和开放性问题

分数阶与 PT 对称在非线性情形的相互作用——目前只在线性情形给出谱特征，非线性的自洽 PT 对称未深入。
$s \to 0$ 的极限（$K_s \to 1/|m|$）——所有格点"全连接"——是否对应一个"超均匀"或"平均场"极限？
Anderson 局域化的"驼峰"机制——为什么 $0.15 < s < 0.4$ 出现反常迁移率？需要从"重整化群"或"局域化长度"角度给出理论解释。
2D/3D 分数阶 DNLS的完整色散关系——本章只给出 2D，3D 未涉及。
分数阶 + 准晶（quasicrystal）——既非周期又非完全无序，分数阶 + PT 对称 + 准晶是开放方向。
实验实现——分数阶 DNLS 的"长程跳跃"在耦合波导阵列中如何实现？需要非局域耦合的工程方案。
分数阶对拓扑相的影响——分数阶紧束缚模型是否改变拓扑边界态的存在？文献中无系统研究。
分数阶 + 无序 + 非线性三者相互作用——同时引入分数阶、无序、非线性，局域化行为如何？

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

优点： 1. "现象学驱动"的典范——不是从抽象理论出发，而是从可观测现象出发：带宽、PT 对称窗口、透射率、自陷阈值、MSD、PR。 2. 分数阶 + 无序的"驼峰"（§5）是出乎意料的发现——揭示分数阶与经典无序理论的非平凡相互作用。 3. 电传输线（§3）提供实验可实现的物理系统——与 Ch1-8 的抽象模型形成对照。 4. 2D 分数阶离散 Laplacian用修正 Bessel 函数表达（[X1] 2018-2020 的近期工作）——严格的数学基础。 5. Green 函数方法（§4）对分数阶离散方程的束缚态分析——把"分数阶非局域"转化为"动量空间积分"。 6. 5 个主题统一（1D DNLS, 2D DNLS, 传输线, 非线性杂质, Anderson）的"分数阶效应"——分数阶 $s$ 减小 → 带宽减小 + 兼并增加 + 阈值降低 + 局域化增强。

缺点： 1. 5 个主题之间缺乏统一理论——每个主题独立处理，没有共同框架（如统一的有效 Hamiltonian 或统一的有效 Lagrangian）。 2. 数值方法只简单提及 Newton-Raphson、Newton-GMRES、旋转波近似——没有详细分析精度和效率。 3. 与 Ch8 的 IST 方法没有交叉引用——本章是"经验/数值"，Ch8 是"严格解析"，但两章应当互相支撑。 4. 自陷的"长程跳跃"图像过于直觉化——缺乏对"分数阶"与"自陷"对应关系的严格推导。 5. Anderson "驼峰"机制未给出物理解释——只说"长程跳跃可绕过无序"但这是含糊的。 6. PT 对称部分（§2.1.1, §3.2）相对简短——分数阶与 PT 对称的相互作用细节不够。 7. 与 Ch3 的"时间分数阶 A.add / A.lin 失稳/稳定反转"是平行现象——但本章未引用 Ch3，可形成"分数阶离散 vs 连续"的对照。

对本人研究（血管生物力学 / SMC 介导的血管重塑）的关联： - 分数阶离散晶格与血管壁的层状结构类比——内弹力板、中膜、外膜可用 1D 离散链建模，分数阶描述层间长程应力传递。 - Anderson 局域化的"驼峰"与血管壁细胞的无序分布相关——SMC 在血管壁中不是均匀分布（有"螺旋"或"螺旋状"排布），分数阶 + 无序可能产生意外的"长程应力绕过"——这与血管的"自适应重塑"机制相关。 - 分数阶自陷阈值降低与血管壁在高血压下的"锁定"现象（血管壁硬化、难以逆转）有类比——分数阶"长程跳跃"使血管壁在病理条件下"自我增强"。 - 2D 分数阶 DNLS与血管壁的"二维网络"（细胞 + 细胞外基质）对应——可作为血管壁生长-重塑的简化模型。 - PT 对称 + 分数阶的"窗口"现象与血管壁的"动态平衡"（生长 vs 退化）类比——PT 对称破缺 = 平衡失稳 = 病理。 - 电传输线与血管壁的"波动传播"（例如脉搏波）的物理图像同源——都是"波在 1D 离散系统中的传播"。

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

[X1] O. Ciaurri, L. Roncal, P.R. Stinga, J.L. Torrea, J.L. Varona. Adv. Math. 330, 688 (2018). —— 2D 分数阶离散 Laplacian 的严格定义。
[X2] D.J. Luitz, et al. —— 长期跳跃 Anderson 模型有 extended states 的存在性证明（[X48-X50] 出现在文中）。
[X3] N. Laskin. Phys. Rev. E 62, 3135 (2000). —— 分数阶量子力学。
[X4] A.A. Kilbas, H.M. Srivastava, J.J. Trujillo. —— 分数阶微积分。
[X5] R. Metzler, J. Klafter. Phys. Rep. 339, 1 (2000). —— 分数阶反常扩散。
[X6] B.B. Mandelbrot. The Fractal Geometry of Nature (Freeman, 1982). —— 分形几何。
[X7] R. Gorenflo, F. Mainardi. —— 分数阶积分的经典工作。
[X8] I. Podlubny. Fractional Differential Equations (Academic Press, 1998). —— Caputo 导数标准教材。
[X9] M. Caputo. Geophys. J. Int. 13, 529 (1967). —— Caputo 导数原始论文。
[X10] M. Riesz. Acta Math. 81, 1 (1949). —— Riesz 位势（Ch2, Ch7 也引）。
[X11] B.J. West. Rev. Mod. Phys. 86, 1169 (2014). —— 分数阶应用于自然科学。
[X12] S. Liu, Y. Zhang, B.A. Malomed, E. Karimi. Nat. Commun. 14, 222 (2023). —— 首次实验分数阶 GVD。
[X13] P. Li, B.A. Malomed, D. Mihalache. Opt. Lett. 46, 3267 (2021). —— PT 对称 FNLS 的 GSs（Ch4 也引）。
[X14] P.W. Anderson. Phys. Rev. 109, 1492 (1958). —— Anderson 局域化原始论文。
[X15] B.A. Malomed. Photonics 8, 353 (2021). —— 分数阶孤子综述。
[X16] M.J. Ablowitz, J.B. Been, L.D. Carr. Phys. Rev. Lett. 128, 184101 (2022). —— ABC 构造（Ch8 也引）。
[X17] C. Klein, N. Stoilov. —— 多区域谱方法（Ch6）。
[X18] M.J. Ablowitz, D.J. Kaup, A.C. Newell, H. Segur. Stud. Appl. Math. 53, 249 (1974). —— AKNS（Ch8 也引）。
[X19] P.G. Kevrekidis, J. Cuevas-Maraver (eds.). Fractional Dispersive Models and Applications (Springer, 2024). —— 本书。
[X20] A.M. Turing. Philos. Trans. Royal Soc. London, Ser. B 237, 37 (1952). —— Turing 不稳定（Ch2, Ch3, Ch5, Ch8 也引）。

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