Chapter 4: Symmetry Breaking in Fractional NLSE and FCGL Soliton Dynamics

1. 作者

Pengfei Li（通信作者），太原师范学院物理系，Email: lpf281888@gmail.com
Boris A. Malomed，以色列特拉维夫大学（与 Ch1 同作者）
Dumitru Mihalache，罗马尼亚 Horia Hulubei 核物理与工程研究所

Li 和 Malomed 是分数阶光孤子领域最活跃的合作组之一，特别在 FNLS 对称破缺和 PT 对称系统方面贡献了大量工作。Mihalache 长期在耗散孤子、非线性波传播方面有重要贡献。三人从 2020 年起在分数阶光孤子领域发表了系列工作。

2. 内容概述

本章是 Ch1（光学模拟分数阶）的逻辑延伸，聚焦分数阶非线性薛定谔方程（FNLS）和分数阶复 Ginzburg-Landau 方程（FCGL）中对称破缺（SSB）与孤子动力学。内容分 4 部分：(1) 1D/2D FNLS 的 SSB：双势阱（double-well potential）下的对称/反对称/不对称孤子家族，及其分岔图（pitchfork、twisting loops）；(2) 2D FNLS 的涡旋孤子：用立方-五次（cubic-quintic, CQ）非线性抑制临界塌缩后的稳定涡旋；(3) PT 对称 FNLS 的对称破缺和鬼态（Ghost States, GSs）——SSB 产生复共轭传播常数的解对（$\alpha < 2$ 时独有）；(4) FCGL 方程的耗散孤子传播、相空间分区、同相孤子合并距离对 Lévy 指数 $\alpha$ 的依赖。

3. 核心方程与概念

3.1 1D FNLS 方程与双势阱

1D FNLS 方程（归一化形式，Riesz 分数阶导数）： $$i\frac{\partial \Psi}{\partial \zeta} - \frac{1}{2}\left(-\frac{\partial^2}{\partial \xi^2}\right)^{\alpha/2}\Psi + V(\xi)\Psi + \sigma|\Psi|^2\Psi = 0$$ $\zeta$ 是传播方向，$\xi$ 是横向坐标，$\alpha \in (1,2]$ 是 Lévy 指数，$\sigma = \pm 1$ 表自聚焦/自散焦。

Riesz 分数阶导数（Fourier 乘子形式）： $$\mathcal{F}\left[-\frac{d^2}{d\xi^2}\right]^{\alpha/2}\psi = |k|^\alpha \mathcal{F}\psi(k)$$

双势阱（沿 $\xi$ 方向）： $$-V(\xi) = -V_0\left[\exp\left(-\frac{(\xi + \xi_0)^2}{W_0^2}\right) + \exp\left(-\frac{(\xi - \xi_0)^2}{W_0^2}\right)\right]$$ 典型值 $V_0 = 2$，$\xi_0 = 1.5$，$W_0 = 1.4$。

稳态孤子 $\Psi(\zeta, \xi) = \psi(\xi)e^{i\beta\zeta}$，$\psi$ 满足 $$-\frac{1}{2}(-\partial_{\xi\xi})^{\alpha/2}\psi + V(\xi)\psi + \sigma|\psi|^2\psi - \beta\psi = 0$$ 功率（norm）$P(\beta) = \int|\Psi|^2 d\xi$，可拆为左右两半 $P = P_L + P_R$。

不对称度参数（SSB 的量化）： $$\Theta = \frac{P_L - P_R}{P_L + P_R}$$

稳定性分析（BoG 微扰）：$\Psi = e^{i\beta\zeta}[\psi + u e^{\delta\zeta} + v^* e^{\delta^*\zeta}]$，$\delta$ 为扰动增长率；$\text{Re}(\delta) > 0$ 表明失稳。

数值方法——Newton-conjugate-gradient： $$\mathcal{L}_0\psi(\xi) = 0, \quad \mathcal{L}_0 = -\frac{1}{2}(-\partial_{\xi\xi})^{\alpha/2} + V + \sigma|\psi|^2 - \beta$$ $$\psi_{n+1} = \psi_n + \Delta\psi_n, \quad \mathcal{L}_1^n \Delta\psi_n = -\mathcal{L}_0(\psi_n)$$ 其中 $\mathcal{L}_1^n$ 是 $\mathcal{L}_1$ 在 $\psi_n$ 处的 Fréchet 线性化；用 FFT 高效计算分数阶导数。

3.2 1D 双势阱中的 SSB（核心结果）

自聚焦 Kerr + 1D 双势阱（$\sigma = +1$，$\alpha = 1.1$）—— pitchfork 分岔（图 1）： - 对称孤子：$\text{Re}(\delta) > 0$ 失稳，超临界 pitchfork 处被稳定不对称孤子替代； - 反对称孤子：保持稳定； - 临界功率 $P_{th}$ 随 $\alpha$ 减小而降低（图 4）—— 分数阶衍射增强 SSB。

自散焦 Kerr + 1D 双势阱（$\sigma = -1$，$\alpha = 1.1$）—— 反对称破缺（图 2）： - 对称孤子保持稳定； - 反对称孤子失稳，被稳定不对称孤子替代； - 同样有临界功率随 $\alpha$ 下降的趋势（图 5）。

$\alpha < 1$ 的例外（图 3）：即使 $\alpha = 0.8$（无势阱时将出现超临界塌缩），双势阱的"trap"机制仍能稳定对称/反对称/不对称孤子——这是分数阶对称破缺的"势阱拯救"现象。

立方-五次（cubic-quintic, CQ）非线性 + 双 hump 势（[X1]）： $$i\partial_\zeta\Psi - \frac{1}{2}(-\partial_{\xi\xi})^{\alpha/2}\Psi + V(\xi)\Psi + |\Psi|^2\Psi - |\Psi|^4\Psi = 0$$ $$-V(\xi) = -V_0[\text{sech}^2((\xi + \xi_0)/W_0) + \text{sech}^2((\xi - \xi_0)/W_0)]$$ 两种分岔并存： - 对称孤子：出现双环分岔（double loops）—— twisting 分支与基分支交叉； - 反对称孤子：仍为超临界 pitchfork。

Lévy 指数影响（图 7, 8）： - LI 增大 → twisting 分岔区域扩大（基分支与扭转分支交叉点更宽）； - LI 增大 → pitchfork 区域收缩； - 稳定性区间也对应变化—— LI 增大时 twisting loops 失稳区扩大，pitchfork 稳定区扩大； - 结论：分数阶衍射对对称和反对称孤子有相反效应。

3.3 2D FNLS 中的 SSB 与涡旋

2D FNLS + 饱和非线性（避免立方临界塌缩）： $$i\partial_\zeta\Psi - (-\nabla_\perp^2)^{\alpha/2}\Psi + V(\xi, \eta)\Psi + \sigma\frac{|\Psi|^2\Psi}{1 + S|\Psi|^2} = 0$$ 饱和参数 $S > 0$。

对称 2D 势： $$-V(\xi, \eta) = -V_0\exp\left[-\frac{(\xi + \xi_0)^2 + \eta^2}{\chi_0^2}\right] - V_0\exp\left[-\frac{(\xi - \xi_0)^2 + \eta^2}{\chi_0^2}\right]$$

关键发现（图 10, 11）： - 自聚焦情形：pitchfork 普适存在，$P_{th}$ 随 $\alpha$ 减小而下降； - 自散焦情形：$\alpha \ge 1.9$ 时 SSB 完全被抑制；$\alpha = 1.9$ 残留，$\alpha = 1.5, 1.1$ 出现分岔环。

2D 涡旋孤子（自由空间 + CQ 非线性，无势阱）（[X2]）： $$i\partial_\zeta\Psi - (-\nabla_\perp^2)^{\alpha/2}\Psi + |\Psi|^2\Psi - |\Psi|^4\Psi = 0$$ 涡旋 ansatz：$\Psi(\zeta, r, \theta) = U(r)e^{is\theta}e^{i\beta\zeta}$，$s$ 是拓扑荷（涡量）。相位 $\phi(r,\theta) = s\theta$。

涡旋稳定性（图 12）： - VK 判据 $d\beta/dP > 0$（$\alpha = 2$ 标准 NLS 时充要）在分数阶情形下仅为必要条件； - 上稳定分支的极限 $\beta_{max} = 3/16$ 是 CQ 非线性的特征（[X3]）； - 失稳模式：角向分裂（angular instability）—— 高涡量 $s = 2, 3$ 涡旋分裂为多个碎片（图 15）。

3.4 PT 对称 FNLS 的 SSB 与鬼态（Ghost States, GSs）

PT 对称势（特殊类）： $$V(\xi) = g^2(\xi) + hg(\xi) + i\frac{dg}{d\xi}$$ $g(\xi)$ 实偶函数，$h$ 实常数。取 $g(\xi) = 2[\text{sech}(\xi + \xi_0) + \text{sech}(\xi - \xi_0)]$，$h = 0$。

FNLS 方程（Kerr 或 CQ 非线性）： $$i\partial_\zeta\Psi - (-\partial_{\xi\xi})^{\alpha/2}\Psi + V(\xi)\Psi + \sigma|\Psi|^2\Psi - \gamma|\Psi|^4\Psi = 0$$

关键现象——Ghost States（GSs, 鬼态）： - 稳态解：$\Psi = \psi e^{i\beta\zeta}$，$\beta$ 可以是复数 $\beta = \beta_R + i\beta_I$； - SSB 产生两支互为复共轭的 GSs——$\beta$ 的实部对称，虚部反号（图 17e, f）； - GSs 只在 $\alpha < 2$（分数阶）下存在；$\alpha = 2$ 时退化为实传播常数的非 PT 对称孤子（图 18b, d）； - 这是"分数阶衍射诱导的 SSB"——与标准 PT 对称系统（$\alpha = 2$）的 SSB 本质不同。

Kerr 非线性：单一 SSB 临界点，GSs 单点对； CQ 非线性：双 SSB 临界点（双环），GSs 出现在双环之间（图 17c, f）。

3.5 分数阶复 Ginzburg-Landau 方程（FCGL）

1D FCGL 方程（带耗散项）： $$i\partial_z\Psi - \frac{1}{2}(-\partial_x^2)^{\alpha/2}\Psi + |\Psi|^2\Psi - \gamma|\Psi|^4\Psi = iR[\Psi]$$ $$R[\Psi] = -\delta\Psi + \varepsilon|\Psi|^2\Psi - \mu|\Psi|^4\Psi + \check D\Psi$$ $$\check D\Psi = -\rho(-\partial_x^2)^{\alpha/2}\Psi$$ 其中 $\delta > 0$ 线性损耗，$\varepsilon > 0$ 立方增益，$\mu > 0$ 五次损耗，$\rho > 0$ 分数阶扩散强度。

输入高斯 $\Psi = A_0\exp(-x^2/2x_0^2)$ 的传播相图（$(\delta, \alpha)$ 平面，图 19）： - 区域 B：欠阻尼，高斯扩张； - 区域 C：不稳定演化 → 自束缚稳定孤子； - 区域 D：快速形成稳定孤子； - 区域 E：强损耗下快速衰减。

$(\alpha, \varepsilon)$ 相图（图 20）：分数阶衍射增强 + 立方增益减小 → 形成稳定孤子所需的"增益阈值"降低。

同相孤子合并距离（图 21）： - 两个同相高斯在 FCGL 中吸引合并； - 合并距离 $z_{merge}$ vs LI $\alpha$：先减后增，呈"V 形"； - LI 调节等价于微调"等效非线性"。

4. 关键结论

分数阶衍射 $\alpha$ 是 SSB 的关键调节器：$\alpha$ 减小 → 双势阱中临界功率 $P_{th}$ 降低 → SSB 更容易发生。这一规律在 1D 和 2D 中一致。
对称/反对称破缺是 $\alpha$ 的"对偶现象"：自聚焦引起对称破缺，自散焦引起反对称破缺——同一双势阱下两个"分岔"可分别发生。
CQ 非线性引入"双环分岔"——比 Kerr 的 pitchfork 更复杂，且分数阶对两者的影响相反：twisting 区域随 LI 增大而扩，pitchfork 区域随 LI 增大而缩。
2D 涡旋孤子在 CQ 非线性下可在 $\alpha \in (1, 2)$ 范围稳定；高涡量 $s \ge 2$ 涡旋失稳后沿角向分裂为多个碎片。
PT 对称分数阶孤子的 Ghost States（$\alpha < 2$ 独有）——是分数阶衍射诱导的新型 SSB 现象，标准 PT 对称系统（$\alpha = 2$）中不存在。GSs 的复共轭传播常数是分数阶非局部算子的"内禀非厄米性"的几何表现。
FCGL 耗散孤子的稳定性边界（$\delta, \alpha, \varepsilon$ 三参数）形成 4 块相图：分数阶衍射调节与立方增益"等效"，$\alpha$ 减小等价于 $\varepsilon$ 增大。
同相孤子合并距离与 $\alpha$ 的"V 形"关系意味着LI 可作为合并距离的"调节旋钮"——是光通信系统中可调参数的物理基础。

5. 挑战和开放性问题

3D FNLS 的 SSB（Ch4 §5 明确点出"a challenging direction"）：2D 已成功，3D + 分数阶的 SSB 全貌仍待揭示。
暗孤子（dark solitons）在分数阶自散焦介质中的稳定性（Ch4 §5 末提）——分数阶扩散下"暗孤子"的等价性（如 $u_0 - U$ 形式）需要新定义。
PT 对称 GSs 的物理解释——复传播常数意味着孤子指数衰减/增长，是"非厄米分数阶算子"的本征现象。其是否可观测？需要实验验证。
FCGL + 分数阶 + 边界条件的实际模拟仍以 1D 为主；2D/3D + 复杂域 + 分数阶耗散的参数空间巨大但缺乏系统性扫描。
同相孤子合并的"V 形"曲线的物理原因——为何 LI 增大先减后增？是群速度与非线性色散的竞争还是其他？需要解析论证。
Lévy 指数管理（LI management, LIM）作为"控制参数"的实验实现——2023 年的 [X4] Liu et al. 首次实现了分数阶 GVD，但 LIM 用于孤子控制仍是开放。
分数阶 GSs 在 $\alpha \to 2$ 极限下的连续性——非平凡，因为 $\alpha = 2$ 时 GSs 消失。该极限下 GSs 的"消失方式"是连续还是突变？
PT 对称双 hump 势的实验可行性——$V(\xi) = g^2 + ig'$ 在光学中需要复折射率/增益损耗分布，技术实现门槛高。

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

优点： 1. SSB 与 GSs 的现象学非常丰富——pitchfork、twisting、双环、ghost state 是真正的"新现象"，不是分数阶算子简单替代 Laplacian 的延伸。 2. 数值工作非常扎实：图 1-21 共 21 张图，每张图对应一个明确物理结论，参数扫描细致。Newton-conjugate-gradient 方法对高维非线性 eigenvalue 求解效率高。 3. 与 Ch1 的衔接自然——同样的 $\alpha$ 参数、同样的 Riesz 算子，但 Ch1 关注"立方/二次非线性 + 单势阱"，Ch4 关注"双势阱 + SSB"——是 Ch1 内容的自然纵深。 4. PT 对称 + 分数阶的组合是首创性工作——标准 PT 对称系统已研究很多，但分数阶 PT 对称的 GSs 现象是 Li et al. 2021 年首次报告。 5. 涡旋分裂的动力学（图 15）直观展示分数阶对角向对称性破缺的影响。

缺点： 1. 大量重复自己的图和参数——许多图仅 $\alpha$ 不同但物理结论完全一样。章节本可以精简到 60% 的体量。 2. 数学分析不深入——通篇是数值结果，没有用 Ch1 那种变分近似（VA）给出解析公式。Ch1 的 VA 在 SSB 问题中本可发挥更大作用。 3. PT 对称 GSs 的物理解释含糊。复传播常数 → 指数增长/衰减，这与"稳定孤子"的概念矛盾—— GSs 是真的"稳态解"还是"短时瞬态"？需要更多讨论。 4. 稳定性分析的覆盖主要在 $\alpha = 1.5, 1.7$ 等少数值；完整的 $\alpha$ 参数稳定性图谱缺失。 5. FCGL 部分相对简略——只有 5 张图，缺少对耗散孤子相互作用、矢量孤子、暗孤子的全面分析。 6. 与 Ch1 重复较多——FNLS 方程、SSB 概念、VA 工具、双核耦合器等内容与 Ch1 重叠。建议 Ch4 应集中在 Ch1 未涉及的内容（PT 对称、GSs、FCGL）。 7. 没有与其他"非 PT 对称 SSB"系统的对比——例如 Ch1 的双核耦合器 SSB（[X5]）与 Ch4 的双势阱 SSB 在机制上的差异未被讨论。

对本人研究（血管生物力学 / SMC 介导的血管重塑）的关联： - SSB 的"临界功率"是分数阶算子可调性的直接体现。类比到血管：双势阱可类比为血管壁的双层（内弹力板两侧）介质异质性——不同的分数阶阶数表示不同层间的"长程相互作用强度"，SSB 类比"层间不对称性出现"。 - 2D 涡旋的角向分裂直接对应生物流体力学中的涡旋-反涡旋对形成。Lévy 指数调节与血液涡量扩散的非局部性相关。 - PT 对称 GSs 的"复传播常数"思想可用于血管壁的"失稳与稳定"临界态——例如动脉瘤前兆期的"复耗散率"。 - FCGL 的"分数阶"调节与血管壁粘弹性的分数阶 Maxwell 模型是同构的——分数阶阶数 $\alpha$ 越大，对应 Maxwell 模型的"短时记忆"占比越高。 - 同相孤子合并距离的"V 形"曲线与血管壁的剪切波传播类比：分数阶 $\alpha$ 影响剪切波的色散关系。

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

[X1] P.F. Li, C.Q. Dai. Ann. Phys. 532, 2000048 (2020). —— CQ 非线性 + 双 hump 势的 twisting/pitchfork 分岔。
[X2] P.F. Li, B.A. Malomed, D. Mihalache. Chaos Solitons Fractals 137, 109783 (2020). —— 2D 涡旋孤子。
[X3] V.A. Stephanovich, W. Olchawa, E.V. Kirichenko, V.K. Dugaev. Sci. Rep. 112, 15031 (2022). —— CQ 非线性 $\beta_{max} = 3/16$。
[X4] S.L. Liu, Y.W. Zhang, B.A. Malomed, E. Karimi. Nat. Commun. 14, 222 (2023). —— 首次实验实现分数阶 GVD（与 Ch1 同）。
[X5] D.V. Strunin, B.A. Malomed. Phys. Rev. E 107, 064203 (2023). —— 分数阶双核耦合器 SSB（Ch1 详细讨论）。
[X6] P.F. Li, B.A. Malomed, D. Mihalache. Opt. Lett. 46, 3267 (2021). —— PT 对称 FNLS 的 GSs（本章核心）。
[X7] P.F. Li, B.A. Malomed, D. Mihalache. Chaos Solitons Fractals 132, 109602 (2020). —— 双势阱中 1D SSB。
[X8] P.F. Li, R.J. Li, C.Q. Dai. Opt. Express 29, 3193 (2021). —— 2D SSB + 饱和非线性。
[X9] P.F. Li, H. Sakaguchi, L.W. Zeng, X. Zhu, D. Mihalache, B.A. Malomed. Chaos Solitons Fractals 173, 113701 (2023). —— 分数阶 SHG 孤子（Ch1 也引）。
[X10] M. Zhong, L. Wang, P.F. Li, Z.Y. Yan. Chaos 33, 013106 (2023). —— 饱和非线性 + PT 对称的扩展。
[X11] M. Zhong, Z.Y. Yan. Commun. Phys. 6, 92 (2023). —— 2D PT 对称 GSs。
[X12] S. Longhi. Opt. Lett. 40, 1117 (2015). —— 4f 腔分数阶衍射（Ch1 也引）。
[X13] B.A. Malomed. Photonics 8, 353 (2021). —— 分数阶孤子综述。
[X14] N. Laskin. Phys. Rev. E 62, 3135 (2000). —— 分数阶 Schrödinger 方程的奠基论文。
[X15] N. Laskin. Fractional Quantum Mechanics (World Scientific, 2018). —— 分数阶量子力学专著。
[X16] I.S. Aranson, L. Kramer. Rev. Mod. Phys. 74, 99 (2002). —— CGL 方程综述。
[X17] N.N. Rosanov. Spatial Hysteresis and Optical Patterns (Springer, 2002). —— CGL 在光学中的应用。
[X18] M.C. Cross, P.C. Hohenberg. Rev. Mod. Phys. 65, 851 (1993). —— 模式形成经典综述。
[X19] E.M. Wright, G.I. Stegeman, S. Wabnitz. Phys. Rev. A 40, 4455 (1989). —— 标准（非分数）耦合器 SSB 早期工作。
[X20] V. Skarka, N.B. Aleksić, W. Krolikowski, D.N. Christodoulides, S. Rakotoarimalala, B.N. Aleksić, M.R. Belić. Opt. Express 25, 10090 (2017). —— 耗散孤子最新综述。

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