摘要:从微观角度上讲, 单个极端异常波事件可视为可积模型方程的时空局域有理函数解. 本文主要讨论了三类典型的可积谐振相互作用模型(即长波短波谐振方程, 三波谐振相互作用方程, 非线性薛定谔和麦克斯韦-布洛赫方程)的基阶Peregrine异常波解及其相关研究进展; 明确指出了这些基阶异常波解形式具有普适性, 可推广应用到多分量或更高阶的可积模型中; 借助数值模拟, 还展示了共存异常波、互补异常波、以及自感应透明Peregrine孤子等新颖动力学.
关键词: 异常波/
长波短波谐振/
三波谐振相互作用/
麦克斯韦–布洛赫方程

English Abstract

全文HTML

--> --> -->

疯狗浪(rogue wave, RW), 又称畸波(freak wave)、怪波(monster wave)、极端波(extreme wave)、异常波(abnormal wave)、杀人波(killer wave)、水墙(walls of water)等, 最初专指海洋上发生的一类稀有的且破坏力极大的暂态巨浪事件^[1-3]. 它不期而至出现在某个地方, 掀起滔天巨浪, 然后消失得无影无踪, 因此获得了“深海怪兽”的绰号(the real monsters of the deep)^[4]. 二三十年前, 这种神秘的疯狗浪事件还只是老水手之间的传说. 现在, 实时的卫星观测, 辅以先进的理论和实验研究, 已经确切证实了这种怪浪真的存在, 而且经常发生. 最为典型的例子是1995年元旦在北海德劳普纳(Draupner)海上平台检测到的时称“万年一遇”巨浪, 后来称之为“新年浪”, 其掀起的水墙就有26 m之高^[1].
学术界对疯狗浪的研究是落后的, 到今天, 连各方都能接受的严格定义都还没有^[3,5-7], 即使在称呼上, 也是五花八门. 鉴于此类极端波事件研究已扩展至多个学科, 如非线性光学^[8]、流体动力学^[9]、等离子体物理^[10,11]、声学^[12]、玻色-爱因斯坦凝聚^[13-16]、甚至金融学^[17], 再沿用“疯狗浪”的称呼将变得不合时宜. 因此, 本文将采用更中性的术语“异常波”(RW)来统称这类极端波事件.
在自然界, 哪些事件可以称之为异常波呢? 经过半个世纪特别是近十年的研究, 科学家们概括出了异常波所普遍具有的基本特征^[5], 其中包括: 1)具有巨大的峰振幅, 通常是有效波高的2倍以上(在海洋学里面, 有效波高是指某片海域内最高的三分之一海浪的平均高度)^[9,18]; 2)不可预测性, 从某种意义讲, 异常波似乎来无影去无踪^[19]; 3)具有L型的波振幅统计规律, 即异常波的出现概率比常规高斯或瑞利统计所预测的要更频繁^[20,21]. 上述三个基本特征目前已广泛用于异常波现象的学术定义、分析和讨论. 这也决定了异常波的科学研究存在“微观”和“宏观”两种方式.
从微观角度来看, 异常波可视为可积非线性偏微分方程的确定有理函数解, 其形式上表现为一类时间和空间双局域的波包, 以反映其不可预测特性^[9,22,23]. 特别地, 一个偏微分方程的解具有如下特征时可以视为基阶(或一阶)异常波^[23,24]: 1)该解的模可用有理函数来表示, 且在时间和空间上都是局域的; 2)该解在时空上通常表现为一个3倍背景振幅的波峰, 并伴有两个侧洞形成. 典型的例子就是Peregrine孤子, 其由Peregrine教授^[25]于1983年首次提出, 当时作为标准非线性薛定谔(nonlinear Schr?dinger, NLS)方程的基阶有理函数解. 由于该概念能很好地解释实验上观测到的单个极端波事件, 因此Peregrine孤子被普遍认为是最简单的异常波原型^[26], 并相继在非线性光纤^[27]、水波池^[18]、等离子体^[28]以及不规则的海洋波态^[29]中观测到. 但另一方面, 在很多情形下, 各种色散波会相互作用, 产生极端复杂的湍流行为, 继而产生异常波^[30,31]. 很自然, 由于湍流背景涉及大量的自由度, 这时人们需要用统计(即宏观)的方法来研究其中的异常波行为^[32]. 在探索极端波事件本质特征上, 这两种研究方式相辅相成, 相得益彰^[8,9].
本文将关注前者, 从可积模型视角上审视异常波事件的产生和演化动力学. 事实上, 在非线性科学中, 现实的实验设置和传播介质的多样性需要形式更广的可积模型^[23], 而不仅仅是简单的NLS方程. 下面列举几个光学方面的理由. 首先, 在很多光学晶体中, 主导非线性效应的是二次非线性系数, 而不是三次(或Kerr)非线性项 ^[33,34]. 其次, 对于飞秒激光脉冲, 基本方程模型还必须考虑进高阶色散和高阶非线性项, 以便能更准确地描述其传播动力学^[35]. 对于色散管理或非线性管理光学系统(如色散渐减光纤), 人们还需要考虑变系数非线性模型^[36,37]. 此外, 相互作用光场的矢量特性与传播介质的多维度特性也要求我们考虑一些耦合的或者高维的数学模型, 例如, 当考虑双折射光纤的脉冲传播^[38,39]或者光学晶体的光波相互作用^[40]时. 最后, 对于一些耗散系统如锁模光纤激光器^[41,42], 还需要额外考虑进耗散项(如增益^[43]), 尽管此时该模型方程(通常为复数Ginzburg-Landau方程^[44,45])将不再可积. 所有这些现实要求造就了众多的可积(或近可积)数学模型, 并持续推动理论工作者从事此类模型的求解工作^[46-49].
据我们所知, 绝大多数可积模型都允许有理数式的异常波解存在, 除了少数实场的波方程(例如经典Korteweg-de Vries (KdV)方程^[23], Kadomtsev-Petviashvili II 方程^[50]等)外. 这其中比较著名的可积模型包含有NLS方程^[51]、Hirota方程^[52]、Sasa-Satsuma方程^[53,54]、Chen-Lee-Liu方程^[55,56]、Fokas-Lenells方程^[57-59]、复数modified KdV方程^[60-62]、广义NLS方程^[63-66]等, 它们的基阶和高阶异常波解均已获得. 相应地, 耦合或矢量可积模型也得到了广泛的关注, 如Manakov系统^[67-69]及其推广^[70-72]、Davey-Stewartson方程^[73,74]等. 数学上, 对于可积模型而言, 这些异常波解可以借助逆散射变换^[75]、达布(Darboux)变换^[76]、Hirota双线性方法^[77]、Riemann-Hilbert方法^[78,79]、或同宿波尝试法(homoclinic test method)^[80]来求解得到. 至于非可积模型, 人们可以用微扰论^[81,82]或者变分法^[83]来近似求解.
异常波是一个快速发展的领域, 其话题的多样性可参看最新的综述论文^[84,85]. 特别是最近几年, 各种异常波新现象、新特性被预测或观测到. 例如, 继2016年观测到光学暗异常波后^[86], Baronio等^[87]又成功在通信光纤上观测到了暗三姊妹异常波, 其演化动力学与Chen等^[38]在2014年所做的解析预测完全一致. 2017年, Peregrine孤子的普适性得到了实验验证^[88], 次年, 周期背景Peregrine孤子^[89]以及反常Peregrine孤子^[90]等概念又相继提出. 此外, 在异常波机制阐释上, 调制不稳定性(modulation instability, MI)和可积湍流是两个长期存在的话题, 得到了持续的关注^[8,91]. 值得一提的是, 2016年Soto-Crespo等^[32]揭示了孤子湍流在异常波产生上扮演了一个常被忽视的重要角色. 所有这些新发现均加深了人们对异常波本质的理解.
本文将综述几类典型的可积谐振相互作用模型的异常波求解及相关研究进展. 它们分别是长波短波(long-wave short-wave, LWSW)谐振方程^[92], 三波谐振相互作用(three-wave resonant interaction, TWRI)方程^[40,93], 和非线性薛定谔和麦克斯韦-布洛赫(NLS and Maxwell-Bloch, NLS-MB)方程^[94,95]. 前两个可积模型可描述光波之间的谐振相互作用, 后一个模型则描述光波与谐振介质之间的相互作用. 对于每一个可积谐振系统, 文中将给出其Lax对、达布变换、基阶Peregrine RW解, 以及数值模拟验证, 并展示其新颖的异常波动力学. 当然, 还有其他可积谐振方程, 如AB模型^[96]、MTM(massive Thirring model)方程^[97,98]等, 它们也存在异常波解, 但由于篇幅关系, 这里不一一叙述.
本文的结构组织如下. 第1节为引言部分, 给出了异常波的背景知识和最新研究动态. 第2节讨论了LWSW方程的基阶异常波解, 数值展示了共存异常波概念. 第3节讨论了TWRI方程的基阶异常波解, 分析了互补型异常波的动力学. 第4节讨论了NLS-MB方程的基阶异常波解, 演示了自感应透明的光学Peregrine孤子和时空互补的物质波异常波动力学. 第5节为全文总结.

在耦合波动力学中, 长波短波(LWSW)谐振是一个典型的参量谐振过程, 通常发生于高频短波的群速度与低频长波的相速度匹配之时. 早在1972年, Zakharov^[99]在研究等离子体物理中的非线性朗缪尔波坍缩时就引入了LWSW谐振机制 , 五年后, Benney^[100]给出了LWSW相互作用方程. 作为一个重要的物理机制, LWSW谐振现象迅速在物理的各个分支得到了广泛的研究, 如流体动力学中重力波与毛细波之间的相互作用^[101], 负折射率光学介质中的简并基频波与差频波的三波混频过程^[102]等. 特别地, 2010年, Shats等^[103]在实验上直接观测到了毛细波异常波事件, 这一事实推动了人们对LWSW谐振介质中各类异常波事件的研究^[104,105].
LWSW谐振方程是一个基本的可积模型, 其归一化形式可写为^[102]

$\begin{split} & {\rm i}u_z+\frac{1}{2}u_{tt}+\phi{u} = 0, \\ & \phi_z-(\vert{u}\vert^2)_t = 0, \end{split}$

式中$ u(z, t) $和$ \phi(z, t) $分别表示短波和长波分量, z和t为空间和时间变量, 下标则表示偏导数. 该可积方程可等价转变为一个$ 3\times3 $线性本征值问题:

$\begin{split} & {{{R}}_t} = {\rm{}}(\lambda {{{U}}_{\rm{0}}} + {{{U}}_{\rm{1}}}){{R}},\\& {{{R}}_z} = {\rm{}}({\lambda ^{\rm{2}}}{{{V}}_{\rm{0}}} + \lambda {{{V}}_{\rm{1}}} + {{{V}}_{\rm{2}}}){{R}}, \end{split}$

式中$ \lambda $为任意谱参数, ${{{U}}_{\rm{0}}} = {\rm{diag}}({\rm{2}}{\rm i}, {\rm{0}}, {\rm{ - 2}}{\rm i})$, ${{{V}}_{\rm{0}}} = $${\rm{diag}}({\rm{ - 2}}{\rm i}{\rm{/3}}, {\rm{4}}{\rm i}{\rm{/3}}, {\rm{ - 2}}{\rm i}{\rm{/3}})$,

${{R}}(\lambda ) = {\rm{}}\left[ {\begin{aligned}{r(\lambda )}\\{s(\lambda )}\\{w(\lambda )}\end{aligned}} \right],\;\;\;\;{{{U}}_{\rm{1}}} = \left[\!\!\!{\begin{array}{*{20}{c}}{\rm{0}}&{\rm{u}}&{{\rm{ - }}{\rm i}\phi }\\{\rm{0}}&{\rm{0}}&{{\rm{ - }}{u^{\rm{*}}}}\\{{\rm{ - 2}}{\rm i}}&{\rm{0}}&{\rm{0}}\end{array}}\!\!\!\right],$

${{{V}}_{\rm{1}}} = \left[\!\!{\begin{array}{*{20}{c}}0 & -u &0 \\0 &0 & -{u^*}\\0 &0 &0 \end{array}}\!\!\right]\!,\;\;{{{V}}_2} =\! \left[\!\!\!{\begin{array}{*{20}{c}}0 &{\dfrac{{{\rm i}{u_t}}}{2}}& -\dfrac{{{\rm i} |u{ |^2}}}{2}\\ -{u^*}&0 &{\dfrac{{{\rm i}u_t^*}}{2}}\\0 &u&0 \end{array}}\!\!\!\right].$

为简洁起见, 这里所有函数的$ (z, t) $变量已省略掉. 据此可以构造出正确的达布变换式^[106]:

$u = {u_0} + \frac{{2{\rm i}(\lambda - {\lambda ^*})\lambda (\alpha - {\beta ^*})}}{\Delta }{s^*}w,$

$\phi = {\phi _0} + \frac{{8(\lambda - {\lambda ^*})}}{\Delta }\left[ {\alpha |\lambda {|^2}|w{|^2} + {\rm i}\lambda (\alpha - {\beta ^*}){\rm{Im}}({r^*}w)} \right],$

式中$ (u_{0}, \; \phi_{0}) $和$ (u, \; \phi) $分别表示方程(1)的种子解和新解, $ \rm{Im} $表示取虚部,

$\alpha = - \frac{{\lambda - {\lambda ^*}}}{{\lambda + {\lambda ^*}}}{R}{(\lambda )^\dagger }{\sigma _{\rm{1}}}{R}(\lambda ) = {\rm{ - }}{\alpha ^{\rm{*}}},$

$\beta = - {R}{(\lambda )^\dagger }{\sigma _{\rm{1}}}{R}({\rm{ - }}\lambda ),$

$\gamma = - \frac{{\lambda - {\lambda ^*}}}{{\lambda + {\lambda ^*}}}{R}{({\rm{ - }}\lambda )^\dagger }{\sigma _{\rm{1}}}{R}({\rm{ - }}\lambda ) = {\rm{ - }}{\gamma ^{\rm{*}}},$

$\varDelta = \alpha \gamma - |\beta {|^2}.$

在上述公式中, $ \lambda $为某个给定的谱参数值, ${\sigma _1}$为第一个泡利算符的$ 3\times3 $版形式, 星号*表示复共轭, 剑号$ \dagger $表示复共轭转置. 这里已经假设$ s(-\lambda) = s(\lambda) $和$ w(-\lambda) = w(\lambda) $, 但$ r(-\lambda) = r(\lambda)+ $$2\lambda w(\lambda) $. 因此, 容易发现, $ \alpha $, $ \beta $和$ \gamma $满足$ \alpha+\beta = \gamma+\beta^{*} $和$ \lambda(\alpha-\beta^{*})+\lambda^{*}(\alpha+\beta) = 0 $^[106].
利用Lax对(2)式以及达布变换(3)式和(4)式, 人们可以求出LWSW方程(1)的各类RW解. 具体推导过程可参看文献[92,106,107]. 这里仅提供其基阶RW解, 如下:

$u \!=\! {u_0}\!\left\{\!{1 - \frac{{2{\rm i}{n^2}z + 2{\rm i}(m - \omega )(t - mz) + 1}}{{[{n^2} \!+\! {{(m \!-\! \omega )}^2}]\! \left[ {{{(t \!-\! mz)}^2} \!+\! {n^2}{z^2} \!+\! \dfrac{1}{{4{n^2}}}} \right]}}} \!\right\},$

$\phi = b + 2\frac{{{n^2}{z^2} - {{(t - mz)}^2} + \dfrac{1}{{4{n^2}}}}}{{{{\left[ {{{(t - mz)}^2} + {n^2}{z^2} + \dfrac{1}{{4{n^2}}}} \right]}^2}}},$

式中$ u_{0}(z, t) = a\exp\left({\rm i}\omega t-{\rm i}k z\right) $为初始的平面波解, $ k =\omega^{2}/2-b $为其色散关系式, 实数b为长波分量的背景高度(通常可设其值为0). 这里实数m和n由下面两个代数方程确定:

$m + \frac{{{a^2}[{n^2} - {{(m - \omega )}^2}]}}{{{{[{n^2} + {{(m - \omega )}^2}]}^2}}} = 0,$

$\frac{1}{2} + \frac{{{a^2}(m - \omega )}}{{{{[{n^2} + {{(m - \omega )}^2}]}^2}}} = 0.$

异常波解(5)式和(6)式展示了丰富的异常波动力学, 如黑异常波^[92,107] (注: 所谓“黑异常波”是指振幅中心下陷至零的一类有限背景双局域波包结构). 特别地, 在文献[92]中, 作者借助解析预测和数值模拟, 分析了LWSW模型背景场的调制不稳定性(MI), 指出了上述基阶RW波的参数存在区间, 即$ \omega\leqslant3\omega_n/2 =3(2 a^2)^{1/3} /2$, 其恰好位于MI图的基带(baseband)区, 参看文献[92]图5. 此后, Baronio等^[68,108]明确提出了MI的基带理论, 即异常波只能存在于MI的基带区, 而不是通频带(passband)区. 现在已证明, 这一理论可用于解析确定可积甚至非可积非线性系统的异常波解参数存在区间^[23].
值得注意的是, 基阶RW解(5)式和(6)式是可以推广到多分量LWSW方程中去的. 例如, 对于两短波分量LWSW方程^[109]:

$\begin{split}& {\rm i}{u_z} + \frac{1}{2}{u_{tt}} + u\phi = 0,\\& {\rm i}{v_z} + \frac{1}{2}{v_{tt}} + v\phi = 0,\\& {\phi _z} - {(|u{|^2} + |v{|^2})_t} = 0,\end{split}$

其基阶RW解中的u和$ \phi $分量依然用(5)式和(6)式来表示, v短波分量可类似表示为

$v \!=\! {v_0}\left\{\!{1 \!-\! \frac{{2{\rm i}{n^2}z + 2{\rm i}(m - \varOmega )(t - mz) + 1}}{{[{n^2} \!+ \!{{(m \!-\! \varOmega )}^2}]\!\left[ {{{(t \!-\! mz)}^2} \!+\! {n^2}{z^2} \!+\! \dfrac{1}{{4{n^2}}}}\!\right]}}} \right\},$

式中$ v_{0}(z, t) = A\exp\left({\rm i}\varOmega t-{\rm i}K z\right) $表示v分量的初始平面波解 (相应色散关系为$K = \varOmega^{2}/2-b $). 同时, 确定实数m和n的两个代数方程变为:

$m + \frac{{{a^2}[{n^2} - {{(m - \omega )}^2}]}}{{{{[{n^2} + {{(m - \omega )}^2}]}^2}}} + \frac{{{A^2}[{n^2} - {{(m - \varOmega )}^2}]}}{{{{[{n^2} + {{(m - \varOmega )}^2}]}^2}}} = 0,$

$\frac{1}{2} + \frac{{{a^2}(m - \omega )}}{{{{[{n^2} + {{(m - \omega )}^2}]}^2}}} + \frac{{{A^2}(m - \varOmega )}}{{{{[{n^2} + {{(m - \varOmega )}^2}]}^2}}} = 0.$

很显然, 当v场消失时(即$ A = 0 $), 两短波分量LWSW方程(9)约化为基本LWSW方程(1). 相应地, 其基阶解将变回(5)式和(6)式的形式, 同时代数方程(11)和(12)也变回方程(7)和(8)的形式.
相比简单的LWSW谐振系统, 多分量LWSW可积系统将存在有趣的共存异常波现象. 所谓共存异常波指的是在相同的背景场中同时发生两类不同结构的确定性异常波事件, 且它们均是可积模型的有理函数解. 例如, 对于给定的$ a = A = 1 $, $ \omega = 0 $, 和$ \varOmega \!=\! -1.2469 $时, 方程(11)和(12)将给出两组有效值$ (m, n) \!=\! (-1.3514, 0.7803) $和(–0.4287,0.6442). 很显然, 把每组$ (m, n) $值代入(5)式, (6)式和(10)式中将得到基于相同背景参数的两种基阶RW解. 图1左列和中列给出了这两类RW解在初始白噪声微扰情况下的数值模拟结果. 很显然, 在足够的距离里, 它们各自均能完美展开, 尽管这时候MI已然发挥作用.
图 1 数值模拟验证初始白噪声微扰下的基阶RW解(5)式, (6)式和 (10)式的稳定性, 左列图对应$(m, n) =(-1.3514, ~0.7803)$, 中列图对应$(m, n) =(- 0.4287, ~0.6442)$. 右列图显示这两类RW结构在同一背景场中的数值激发. 图改编自文献[109]
Figure1. Simulations confirm the stability of the fundamental RW solutions (5), (6), and (10) against initial white noise perturbations. Left column: $(m, n) =(-1.3514, ~0.7803)$; Middle column: $(m, n) =(- 0.4287, ~0.6442)$. The right column shows the numerical excitation of such two rogue wave families from the same background field. Figure adapted from Ref. [109].

另一方面, 人们也许会问, 这两类RW解能否在现实条件下共存呢? 为此, Chen等^[109]做了广泛的数值激发实验. 具体上, 选取$ z = 0 $的平面波解作为u和v的初始条件, 并令$ \phi = 0.4\cos(2{\text{π}} t/40)\times $${\rm sech}[(t-2)/8] $, 然后利用分步傅里叶算法数值积分模型方程(9), 计算结果如图1右列所示(图中已去掉前五个单位距离的波演化轮廓图, 因为初始阶段MI引起的波结构不是很明显). 可以清晰地看到, 在$ z = 10 $附近, 同一背景中同时出现了两种明显不同的异常波结构, 它们均和方程(9)的解析解完全一致. 随后, 由于MI不断增强, 背景场中将出现多异常波动力学, 但所有这些动力学都是上述两种异常波类型的组合.

众所周知, 三波谐振相互作用(TWRI)在非线性科学中扮演着重要的角色^[110]. 例如, 在非线性光学中, TWRI可描述不同的光学过程, 如参量放大^[111]、频率转换^[112]、受激拉曼散射^[113]、受激发布里渊散射^[114]等. 此外, TWRI还可用来实现光脉冲的群速度控制^[115]、超短脉冲列产生^[116]、参量三波孤子产生^[117]以及激光-等离子体相互作用^[118]等. 早在20世纪70年代, 人们就建立了TWRI控制方程的可积性, 并给出了其孤子解^[110]. 不同于大家所熟悉的二次孤子情形^[119], 此类相干孤子主要产生于由非线性效应引起的能量交换与由群速度不匹配所引起的对流之间的平衡^[40]. 它们最终能以一个共同的(锁定的)速度传播, 尽管其三个波分量在相互捕获之前的群速度可以互不相同^[120]. 这种特性使得这类TWRI孤子在应用中相当诱人, 这是因为由群速度不匹配而引起的走散效应(walk-off)此时能被非线性耦合抵消掉.
在弱色散近似下, TWRI方程的基本形式可表示为^[110]:

$\begin{split}& {u_{1t}} + {V_1}{u_{1z}} = u_2^*u_3^*,\\& {u_{2t}} + {V_2}{u_{2z}} = - u_1^*u_3^*,\\& {u_{3t}} + {V_3}{u_{3z}} = u_1^*u_2^*,\end{split}$

式中$ u_{1, 2, 3}(z, t) $是三个光场的慢变复包络函数, $ V_{1, 2, 3} $为其相应的群速度常数. 下面不妨假定$ V_{3} = 0 $, 此即意味着模型(13)是建立在随$ u_{3} $场运动的参照系上. 一般地, 当$ V_1>V_2 $时, 该TWRI模型将允许孤子交换(soliton exchange)动力学(在非线性光学语境中, 又称为参量三波混频过程)得以存在, 但当$ V_1<V_2 $时, 其将展示有趣的受激背散射(stimulated backscattering)动力学^[110].
该方程具备完全可积性^[110], 因此具有下面的Lax对形式:

${{{R}}_t} = {\rm{}}(\lambda {{{U}}_{\rm{0}}} + {{{U}}_{\rm{1}}}){{R}},\;\;\;\;\;{{{R}}_z} = {\rm{}}({\rm{ - 3}}{\rm i}\lambda {{{V}}_{\rm{0}}} + {{{V}}_{\rm{1}}}){{R}},$

式中$ \lambda $为任意谱参数, ${{{U}}_{\rm{0}}} = {\rm{diag}}({\rm{ - 2}}{\rm i, i, i})$, $ {{{V}}_{\rm{0}}} = $ ${\rm{diag}}({\rm{0}}, {\rm{1/}}{V_{\rm{1}}}, {\rm{1/}}{V_{\rm{2}}})$,

${\bf{R}} = \left[ {\begin{aligned}r\\s\\w\end{aligned}} \right],\;\;\;\;{{{U}}_{\rm{1}}} = \left[\!\!{\begin{array}{*{20}{c}}{\rm{0}}&{\dfrac{{{V_{\rm{1}}}u_{\rm{1}}^{\rm{*}}}}{{{V_{\rm{2}}}{\rm{ - }}{V_{\rm{1}}}}}}&{\dfrac{{{V_{\rm{2}}}{u_{\rm{2}}}}}{{{V_{\rm{2}}}{\rm{ - }}{V_{\rm{1}}}}}}\\{{u_{\rm{1}}}}&{\rm{0}}&{\rm{0}}\\{u_{\rm{2}}^{\rm{*}}}&{\rm{0}}&{\rm{0}}\end{array}}\!\!\!\right],$

${{{V}}_{\rm{1}}} = \left[ {\begin{array}{*{20}{c}}{\rm{0}}&{\dfrac{{u_{\rm{1}}^{\rm{*}}}}{{{V_{\rm{1}}}-{V_{\rm{2}}}}}}&{\dfrac{{{u_{\rm{2}}}}}{{{V_{\rm{1}}}-{V_{\rm{2}}}}}}\\{-\dfrac{{{u_{\rm{1}}}}}{{{V_{\rm{1}}}}}}&{\rm{0}}&{\dfrac{{u_{\rm{3}}^{\rm{*}}}}{{{V_{\rm{1}}}}}}\\{-\dfrac{{u_{\rm{2}}^{\rm{*}}}}{{{V_{\rm{2}}}}}}&{-\dfrac{{{u_{\rm{3}}}}}{{{V_{\rm{2}}}}}}&{\rm{0}}\end{array}} \right].$

相应地, 其达布变换式可以表示为^[121]:

$\begin{split}& {u_1} = {u_{10}} + \frac{{3{\rm i}(\lambda - {\lambda ^*}){r^*}s}}{{{{R}^\dagger }{\sigma _{\rm{3}}}{R}}},\\& {u_2} = {u_{20}} + \frac{{3{\rm i}(\lambda - {\lambda ^*})r{w^*}}}{{{{R}^\dagger }{\sigma _{\rm{3}}}{R}}},\\& {u_3} = {u_{30}} + \frac{{3{\rm i}(\lambda - {\lambda ^*}){s^*}w}}{{{{R}^\dagger }{\sigma _{\rm{3}}}{R}}},\end{split}$

式中$ u_{n0} $和${u_n}\left( {n = 1, 2, 3} \right)$分别表示方程(13)的种子解和新解, $ \lambda $为某个任意给定的谱参数值, $ \sigma_{3} = \mathrm{diag}(1, \varGamma_1, \varGamma_2) $, ${\varGamma _j} = {V_j}/({V_1} - {V_2})\left( {j = 1, 2} \right)$.
考虑到三波谐振相互作用条件(即动量和能量守恒), 初始的平面波种子解可以表示为:

$\begin{split}& {u_{10}}(z,t) = {a_1}\exp \left[ { - {\rm i}({k_1}z - {\omega _1}t)} \right],\\& {u_{20}}(z,t) = {a_2}\exp \left[ {{\rm i}({k_2}z - {\omega _2}t)} \right],\\& {u_{30}}(z,t) = i{a_3}\exp \left[ {{\rm i}({k_1} - {k_2})z - {\rm i}({\omega _1} - {\omega _2})t} \right],\end{split}$

式中$ k_{1} = \dfrac{\omega_{1}}{V_{1}}+\dfrac{a_{2}^{2}}{\delta V_{1}} $, $ k_{2} = \dfrac{\omega_{2}}{V_{2}}+\dfrac{a_{1}^{2}}{\delta V_{2}} $, $ a_{3} = \dfrac{a_{1}a_{2}}{\delta} $. 这里$ a_{1, 2, 3} $表示三波各自的背景高度, $ \delta = \omega_{1}-\omega_{2} $ (不失一般性, 下面假设$ \delta>0 $).
利用上面的Lax对(14)式和达布变换(15)式, 很容易得到TWRI方程的基阶RW解^[93](借助达布变换求解基阶或高阶异常波解的详细过程可参看我们近期发表的论文[64,67,72]):

$\begin{split}& {u_1} = {u_{10}}\left[ {1 + \frac{{3({\lambda _0} - \lambda _0^*)\xi \theta _1^*/\alpha _1^*}}{{|\xi {|^2} + {\varGamma _1}a_1^2|{\theta _1}/{\alpha _1}{|^2} + {\varGamma _2}a_2^2|{\theta _2}/{\alpha _2}{|^2}}}} \right],\\& {u_2} = {u_{20}}\left[ {1 - \frac{{3({\lambda _0} - \lambda _0^*){\xi ^*}{\theta _2}/{\alpha _2}}}{{|\xi {|^2} + {\varGamma _1}a_1^2|{\theta _1}/{\alpha _1}{|^2} + {\varGamma _2}a_2^2|{\theta _2}/{\alpha _2}{|^2}}}} \right],\\& {u_3} = {u_{30}}\left[ {1 + \frac{{3({\lambda _0} - \lambda _0^*)\delta {\theta _1}\theta _2^*/({\alpha _1}\alpha _2^*)}}{{|\xi {|^2} + {\varGamma _1}a_1^2|{\theta _1}/{\alpha _1}{|^2} + {\varGamma _2}a_2^2|{\theta _2}/{\alpha _2}{|^2}}}} \right],\end{split}$

式中$ \xi = t+\dfrac{z}{V_{1}-V_{2}}\left (\dfrac{a_{1}^{2}}{\alpha_{1}^{2}}+\dfrac{a_{2}^{2}}{\alpha_{2}^{2}}\right) $, $ \alpha_{j} = \mu_{0}+\lambda_{0}- $$(-1)^{j}\delta/2 $, ${\theta _j} = \xi -{\rm i}/{\alpha _j}\left( {j = 1, 2} \right)$. 这里, $ \mu_{0} $是立方代数方程$ \mu^3-3\sigma\mu+2\rho = 0 $的一个等根(可实可复), $ \lambda_{0} $是等根判别式$ \varDelta = \sigma^{3}-\rho^{2} = 0 $的一个复根(虚部不为0), 式中$ \sigma \!=\! \lambda_{0}^2+\delta^{2}/12\!+\!\dfrac{1}{3}( \varGamma_{1}a_{1}^{2}\!+\!\varGamma_{2}a_{2}^{2}) $, $ \rho = \dfrac{1}{4}[\delta^{2}-2(\varGamma_{1}a_{1}^{2}+\varGamma_{2}a_{2}^{2})]\lambda_{0}\!-\!\lambda_{0}^3+\dfrac{1}{4}\delta (\varGamma_{1}a_{1}^{2}\!-\!\varGamma_{2}a_{2}^{2}) $. 可以证明, 不管是三波混频($ V_1>V_2 $)还是受激背散射($ V_1<V_2 $)过程, 在特定参数条件下, 解(17)式均能展示出标准的Peregrine孤子形式^[93].
毫无疑问, 人们会问, TWRI方程的简并情形(此时$ V_1 = V_2 = V $)是否存在RW解呢? 这似乎很难从上面的达布变换式寻求答案, 因为在这种情形下, Lax对(14)式将变得无意义. 但显然, 简并TWRI方程依然是可积的, 因此其RW解是可能存在的. 最近, 文献[122]采用另一个策略来回答这个问题, 也就是, 对广义解(17)式取如下极限: $ V_2\rightarrow V_1 = V $, 成功得到了简并TWRI方程的基阶RW解:

$\begin{split}& {u_1} = {u_{10}}\left[ {1 + \frac{{2{\rm i}\delta V(\theta - 2a_2^2z) - {\delta ^2}A{V^2}/a_1^2}}{{{\theta ^2} + 4{\delta ^2}a_3^2{z^2} + {A^2}{V^2}/(4a_3^2)}}} \right],\\& {u_2} = {u_{20}}\left[ {1 + \frac{{2{\rm i}\delta V(\theta + 2a_1^2z) - {\delta ^2}A{V^2}/a_2^2}}{{{\theta ^2} + 4{\delta ^2}a_3^2{z^2} + {A^2}{V^2}/(4a_3^2)}}} \right],\\& {u_3} = {u_{30}}\left[ {1 - \frac{{4{\rm i}{\delta ^3}V(Vt - z) + {A^2}{V^2}/a_3^2}}{{{\theta ^2} + 4{\delta ^2}a_3^2{z^2} + {A^2}{V^2}/(4a_3^2)}}} \right],\end{split}$

式中$ A = a_{1}^{2}+a_{2}^{2} $, $ \theta = \delta^{2}(Vt-z) -(a_{1}^{2}-a_{2}^{2})z $. 容易证明:

$|{u_1}{|^2} + |{u_2}{|^2} = a_1^2 + a_2^2 = A,$

此意味着$ u_1 $分量和$ u_2 $分量的强度和总是不变的, 不管它们各自的时空结构如何变化. 换句话说, 它们是时空互补的, 因此称为互补异常波^[122].
文献[122]给出了互补型异常波的数值模拟结果, 见图2. 具体来说, 作者将三个场分量的实部和虚部分别乘以因子$ [1+ \varepsilon r_i(z)] $ (这里, $ i = 1, \cdots, 6 $, $ r_i $为均值为0、方差为1的随机分布函数, $ \varepsilon $为噪声强度参数), 然后采用分步Fourier算法对TWRI模型进行数值积分. 图2左列为未添加白噪声($ \varepsilon = 0 $)的模拟结果, 右列为添加了白噪声($ \varepsilon = $10^–8)的模拟结果. 可以看到, 在微小的扰动下, 所有三个异常波分量依然可以在相当长的时间内稳定传播, 直到连续波背景的MI显著增长为止.
图 2 互补型基阶RW解(18)式的数值模拟结果. 左列图: 未微扰情形; 右列图: 白噪声微扰情形. 图摘自文献[122]
Figure2. Simulation results of the complementary fundamental rogue wave solutions (18). Left column: unperturbed; Right column: perturbed by initial white noises. Figure adapted from Ref. [122].

实验上, 我们预期这类亮-暗型的双色互补型异常波有可能在双模光纤上得以实现. 这是因为, 光波在双模光纤里传播时, 很容易发生前向受激布里渊散射过程^[123]; 在这种情况下, 泵浦光和斯托克斯(Stokes)光的群速度几乎相同, 而相比之下, 声波的速度可近似为零, 三者满足简并TWRI方程所需要的模间耦合条件, 从而出现上面提到的双色互补型异常波动力学^[122].

另一方面, 光波与非线性谐振介质的相互作用也是一个经久不衰的光学研究话题(注意区分上面讲的LWSW, TWRI过程, 其通常发生在非谐振介质的两波或三波谐振相互作用)^[124]. 而这其中最有名的莫过于光脉冲与两能级原子或离子的相互作用, 即所谓的麦克斯韦-布拉赫(MB)耦合^[125]. 正是由于这种谐振相互作用, 原本不透明(或吸收)的介质在超短激光脉冲照射下将变得完全透明, 此即大家所熟知的自感应透明现象^[126]. 除了自感应透明现象外, MB方程还可以产生面积为$ 2{\text{π}} $的基本sech型孤子^[124]. 后来, 一些****对MB系统做了推广, 研究了掺铒光纤中的脉冲传播动力学^[127]. 在慢变包络近似和旋波近似下, 该系统可用非线性薛定谔-麦克斯韦布洛赫(NLS-MB)方程来描述. 下面就来讨论这一有趣的可积模型及其基阶RW解.
为方便讨论, 把NLS-MB方程写成如下无量纲形式^[94,95]:

$ \begin{split} & {u_z} = \frac{\rm i}{2}{u_{tt}} + {\rm i}|u{|^2}u + M,\\ & {M_t} = {\rm i}\phi M - 2uF,\\ & {F_t} = {M^*}u + M{u^*}, \end{split} $

式中$ u(z, t) $表示光波复振幅, M和F表示物质波的特征函数(具体讲, M由谐振介质密度矩阵的离对角元$ \rho_{12} $来确定, 为复函数; F表示上下能级的布居差, 为实函数, 二者满足$ |M|^2+F^2 = 1 $^[95]), $ \phi $为去谐频率常数.
同样, 该可积方程可等价转变为一个$ 2\times2 $线性本征值问题(Lax对):

${{{R}}_t} = {{UR}},\;\;\;\;\;{{{R}}_z} = {{VR}},$

式中各矩阵形式表示如下:

${\bf{R}} = \left[ {\begin{aligned}r\\s\end{aligned}} \right],\;\;\;{{U}} = \left[\!\!\!{\begin{array}{*{20}{c}}{{\rm i}\lambda }&u\\{{\rm{ - }}{u^{\rm{*}}}}&{{\rm{ - }}{\rm i}\lambda }\end{array}}\!\!\right],$

${{V}} = \left[\!\!\!{\begin{array}{*{20}{c}}{{\rm{ - }}{\rm i}{\lambda ^{\rm{2}}} + \dfrac{{{\rm i}{\rm{|}}u{{\rm{|}}^{\rm{2}}}}}{{\rm{2}}} + \dfrac{{{\rm i}F}}{{{\rm{2}}\lambda {\rm{ - }}\phi }}}&{{\rm{ - }}\lambda u + \dfrac{{{\rm i}{u_t}}}{{\rm{2}}} + \dfrac{{{\rm i}M}}{{{\rm{2}}\lambda {\rm{ - }}\phi }}}\\{\lambda {u^{\rm{*}}} + \dfrac{{{\rm i}u_t^{\rm{*}}}}{{\rm{2}}} + \dfrac{{{\rm i}{M^{\rm{*}}}}}{{{\rm{2}}\lambda {\rm{ - }}\phi }}}&{{\rm i}{\lambda ^{\rm{2}}}{\rm{ - }}\dfrac{{{\rm i}{\rm{|}}u{{\rm{|}}^{\rm{2}}}}}{{\rm{2}}}{\rm{ - }}\dfrac{{{\rm i}F}}{{{\rm{2}}\lambda {\rm{ - }}\phi }}}\end{array}}\!\!\!\right].$

据此, 可以构造正确的达布变换式^[128]:

$\begin{split}& u \!=\! {u_0} + \frac{{2{\rm i}(\lambda - {\lambda ^*})r{s^*}}}{{|r{|^2} + |s{|^2}}},\\& M \!=\! \frac{{{\varDelta ^2}{M_0} \!-\! 4{{(\lambda \! - \!{\lambda ^*})}^2}{{(r{s^*})}^2}M_0^* \!-\! 4(\lambda \!-\! {\lambda ^*})r{s^*}\varDelta {F_0}}}{{|2\lambda - \phi {|^2}{{(|r{|^2} + |s{|^2})}^2}}},\\& F \!=\! \frac{{\left[ {|\varDelta {|^2}\! +\! 4{{(\lambda \!- \!{\lambda ^*})}^2}|rs{|^2}} \right]\!{F_0}\! - \!4{\rm i}(\lambda \! -\! {\lambda ^*}){\rm{Im}}({r^*}s\varDelta {M_0})}}{{|2\lambda - \phi {|^2}{{(|r{|^2} + |s{|^2})}^2}}},\end{split}$

式中$ \varDelta = ( 2\lambda-\phi) |r|^2+(2\lambda^{*}-\phi)|s|^2 $, $ \lambda $为某个给定的谱参数值, $ (u_0, M_0, F_0) $和$ (u, M, F) $分别表示NLS-MB方程(20)的旧解和新解.
利用上面的Lax对和达布变换式, 并参照文献[64,67,72]的做法, 可以很容易推得 NLS-MB方程(20)的基阶RW解^[95]:

$\begin{split}u =\, & {u_0}\left( {1 - \frac{{16{\rm i}\eta z + 4}}{{16{\eta ^2}{z^2} + 4{a^2}{\theta ^2} + 1}}} \right),\\M =\, & {M_0}\left\{ {1 - \frac{{16{\rm i}{\kappa ^2}\eta z + 4{\kappa ^2}}}{{({\kappa ^2} + 4{a^2})(16{\eta ^2}{z^2} + 4{a^2}{\theta ^2} + 1)}}} \right.\\ &\left. {- \frac{{64{\rm i}(4{a^2}{\theta ^2} - {\vartheta ^2}){a^2}\eta z - 32{\rm i}{a^2}\kappa \vartheta \theta }}{{({\kappa ^2} + 4{a^2}){{(16{\eta ^2}{z^2} + 4{a^2}{\theta ^2} + 1)}^2}}}} \right\},\\F =\, & \frac{{\kappa b}}{{2a}} - \frac{{8ab[(16{\eta ^2}{z^2} - 4{a^2}{\theta ^2} + 1)\kappa - 32{a^2}\eta z\theta ]}}{{({\kappa ^2} + 4{a^2}){{(16{\eta ^2}{z^2} + 4{a^2}{\theta ^2} + 1)}^2}}},\end{split}$

式中$ \kappa = \omega- \phi $, $ \vartheta = 4 {\rm i}\eta z+1 $, $ \theta = t-\chi z $, 以及

$\eta = - \frac{{ab}}{{{\kappa ^2} + 4{a^2}}} + \frac{{{a^2}}}{2},\;\;\;\;\chi = \frac{{\kappa b}}{{a({\kappa ^2} + 4{a^2})}} + \omega .$

这里初始的平面波解$ u_{0}(z, t) $, $ M_{0}(z, t) $和$ F_{0}(z, t) $可定义为

$ \begin{split} &{u_0} = a\exp \left[ {{\rm i}\left( {kz + \omega t} \right)} \right],~ {M_0} = \frac{{\rm i}b}{a}{u_0},~{F_0} = \frac{{\kappa b}}{{2a}}, \end{split} $

$k = \frac{b}{a} + {a^2} - \frac{{{\omega ^2}}}{2}.$

另外, 很容易证明, F和M的时空分布满足$F{(z, t)^2} + |M(z, t){|^2} = \left( {\dfrac{{{\kappa ^2}}}{{4{a^2}}} + 1} \right){b^2}$. 这表明物质波组分具有与简并TWRI系统相似的互补异常波特性^[122]. 尽管上述RW解适用于任意b值, 但考虑到概率守恒条件$ F^2 + |M|^2 = 1 $, 可令 b = $2 a/\sqrt{\kappa^{2}+4 a^{2}} $, 其值取决于参数a和$ \kappa $的取值.
图3(a)展示了这些解析解的时空分布图, 其初始参数为$ a = 1.5 $, $ \phi = 1/2 $, $ \omega = 0 $. 可以看出, 光场分量显示为一个标准的Peregrine孤子结构; 比较而言, 物质波分量M和F展示出更复杂的时空结构, 但二者满足$ F^2 + |M|^2 = 1 $, 即时空互补性. 此外, 为了评估异常波动力学及其稳定性, Chen等^[95]采用指数时间差分Crank-Nicolson (ETDCN)算法对模型方程(20)执行了数值模拟, 结果如图3(b)所示. 这里使用$ u(z = -1, t) $作为光场分量的初始条件, 而让物质波分量M和F分别取$ M_0(z = -1, t) $和$ F_0 $作为其初始值. 这些初始条件对应于对解析解(23)式的强烈扰动. 结果表明, 尽管自发性MI的发生呈指数增长, 并且倾向于干扰局域解的尾部部分, 但是异常波结构, 特别是对于光场分量, 均可以在相当长的距离上展开而不失真. 进一步地, 为了验证这类RW结构能否在真实条件下产生, 他们把光场的初始条件也换成平面波输入, 考察此类RW解的数值激发情况. 图3(c)为其数值模拟结果, 清晰显示了这些典型的RW结构是可以激发产生的, 见图中黑圈标出部分. 这些数值结果充分预示了实验观测的可行性, 具体实验方案可参看文献[95].
图 3 NLS–MB方程的基阶RW解(23)的时空演化, 其中(a)列图对应解析解的3D曲面和轮廓图; (b)列图为数值模拟结果, 初始条件已文中给出; (c)列图显示这类异常波结构在背景场中的数值激发产生, 已黑线圈出. 图改编自文献[95]
Figure3. Spatiotemporal evolution of the fundamental rogue wave solutions (23) of the NLS–MB equation. Column (a): Analytical solutions, given by 3D surface and contour plots; Column (b) the numerical results, with initial conditions being specified in the text; The column (c) shows the numerical excitation of the rogue waves, indicated by the black circles, from the background field. Figure adapted from Ref. [95].

最后指出, 这里给出的基阶RW解(23)式具有普适性, 可以推广应用到高阶NLS-MB方程中去. 例如, 对于下面的Hirota-MB 方程^[129]:

$\begin{aligned}& {{u_z} = {\rm i}g\left( {\frac{1}{2}{u_{tt}} + |u{|^2}u} \right) + h\left( {\frac{1}{6}{u_{ttt}} + |u{|^2}{u_t}} \right) + M,}\\& {{M_t} = {\rm i}\phi M - 2uF,}\\& {F_t} = {M^*}u + M{u^*},\\[-15pt]\end{aligned}$

(23)式依然可以作为其基阶RW解, 只需要把(24)式和(26)式定义的k, $ \eta $, 和$ \chi $参数换成如下形式即可:

$ \begin{aligned} & {k = \frac{b}{a} + g\left( {{a^2} - \frac{{{\omega ^2}}}{2}} \right) + \frac{h\omega (6{a^2} - {\omega ^2})}{6},}\\ & {\eta = - \frac{{ab}}{{{\kappa ^2} + 4{a^2}}} + \frac{{g{a^2}}}{2} + \frac{{h\omega {a^2}}}{2},}\\ & {\chi = \frac{{\kappa b}}{{a({\kappa ^2} + 4{a^2})}} + g\omega + \frac{{h({\omega ^2} - 2{a^2})}}{2}.} \end{aligned} $

这里g和h为任意实常数. 不难看出, 当$ h = 0 $和$ g = 1 $时, 上述 Hirota-MB方程及其RW解就回到了NLS-MB情形.

本文系统综述了LWSW谐振方程, TWRI方程和NLS-MB方程这三类典型可积谐振模型的新颖异常波动力学及其研究进展. 首先, 对于LWSW谐振系统, 提供了其基阶RW解的一般形式, 指出长波和短波之间的谐振相互作用能导致黑异常波的产生. 特别地, 该基阶RW解推广到多分量LWSW系统时, 可以产生有趣的共存异常波现象, 并得到了数值模拟验证. 其次, 在TWRI系统中, 给出了该可积模型的基阶RW解, 指出该解可以适用三波混频($ V_1>V_2 $)和受激背散射($ V_1<V_2 $)两种情形. 特别地, 对于简并TWRI情形(V₁ = V₂ = V), 也提供了其精确基阶RW解, 并数值展示了其互补异常波动力学. 最后, 在NLS-MB方程框架内, 讨论了光场与低温光纤中的两能级掺杂离子的谐振相互作用, 给出了光学和物质波的基阶RW解一般形式, 数值展示了它们的新颖异常波动力学. 特别地, 指出上述基阶RW解可以推广应用到更高阶的NLS-MB耦合系统(如Hirota-MB系统), 展示了该解的普适性.
这里附带提一下. 本文讨论的简并TWRI模型在一定条件下可以转换成著名的sine-Gordon方程^[130], 而后者与其他可积谐振模型如AB模型、MTM方程也有类似的转换关系^[131,132]. 因此可以预测, 本文呈现的基阶Peregrine RW解对其他可积谐振系统的异常波动力学研究也将有参考借鉴作用.
感谢Philippe Grelu教授、Jose M. Soto-Crespo教授、Dumitru Mihalache教授、Stefan Wabnitz教授、上海理工大学刘一教授富有成效的合作和深入的讨论.