激光与光电子学进展, 2019, 56 (5): 050601, 网络出版: 2019-07-31  

非线性薛定谔方程的呼吸子解及其怪波极限 下载: 1341次

Breather Solutions and Their Rouge Wave Limits of NonlinearSchrödinger Equation
作者单位
山西大学物理电子工程学院, 山西 太原 030006
摘要
采用达布变换法得到了标准非线性薛定谔方程的一阶呼吸子解及其怪波极限,研究了一阶呼吸子解的动力学特性。借助达布变换的递推关系得到了非线性薛定谔方程的高阶呼吸子解,并分别研究了碰撞叠加、分离、简并和并行传输模式。当各呼吸子的频率趋于零时,得到非线性薛定谔方程怪波极限。研究结果表明,怪波幅值、凸起数以及怪波分裂后中心波峰的阶数和周围的波峰个数均与怪波阶数有关。
Abstract
Based on the standard nonlinear Schr dinger equation, the first-order breather solution and its rouge wave limit are obtained with Darboux transform method, and the dynamic characteristics of first-order breather solution are studied. High-order breather solutions of nonlinear Schr dinger equation are obtained by means of recurrence relation of Darboux transformation. And their collision superposition, separation, degeneracy and parallel transmission modes are studied, respectively. Nonlinear Schr dinger equation's rouge wave limit can be obtained when each breather frequency tends to zero. Research results show that the rouge wave's amplitude, number of bumps, order of center peaks and number of surrounding peaks after splitting are related to rouge wave's order.

1 引言

在现代科学与技术的发展过程中,系统的非线性特性一直受到广泛关注,非线性特性是理解很多自然现象的基础,在流体力学、化学、生物和非线性光学等领域均可以观察到,其重要性不容置疑。一般情况下,非线性现象利用线性和非线性效应的非线性演化方程进行建模,标准非线性薛定谔方程(NLSE)就是描述非线性现象的核心模型之一[1-2]

自20世纪60年代以来,非线性薛定谔方程已成为广泛研究的偏微分方程,一维自聚焦的非线性薛定谔方程被认为是描述海洋波的动力学模型[3-10]。非线性薛定谔方程之所以被广泛研究,是因为其在海洋动力学之外的其他物理学分支中均有重要应用,如非线性光学[11-12]、Bose-Einstein凝聚[13]以及其他非线性物理学领域。NLSE的精确解包括孤子[1]、Peregrine怪波(RW)[14-15]、Akhmediev呼吸子(AB)[16-17]、Kuznetsov-Ma孤子(KMS)[18-19]以及这些解的不同组合[20]

目前,非线性薛定谔方程的呼吸子解被认为是在海洋以及其他物理领域中怪波的潜在原型[21-23],当呼吸子周期趋于无穷大时,呼吸子解将成为有理分式的怪波解,因此,有关呼吸子的结果也将为怪波研究提供理论参考。1991年,呼吸子解首次被Kichenassamy[24]求得;随后几十年,呼吸子解分别在非线性薛定谔方程[25]、Korteweg-de Vries方程[26-27]、二维Fermi-Pasta-Ulam晶格[28]、向列相液晶盒[29]和掺杂光纤[14]等系统中获得;呼吸子由调制不稳定性(MI)[14,30-31]导致的周期性扰动所产生,这些扰动通常是混沌的,可能在频谱中包含多个频率成份,因此,如何构建高阶解成为研究者非常关心的问题,Akhmediev小组认为高阶解可以是几个低阶解的非线性叠加[32-33]

在自然界中,怪波是一种突然出现但短时间内又突然消失并具有极强破坏力的波[34],然而这种破坏性至今无法完全解释,航海家们一致认为线性理论无法解释怪波[35],只有非线性理论才能解释这种能量集中且峰值幅度通常为有效波高两倍以上的现象。在最近的研究中,除NLSE方程之外,Hirota方程[36]、Fokas-Leneels方程[37]、Hirota Maxwell-Bloch(MB)方程[38]、Sasa-Satsuma方程[39]以及耗散系统[40]中均已证明怪波的存在。经过研究,科学家们一致认为调制不稳定性是形成怪波最根本的动力学原因[41],而高阶怪波是由具有相同或特定频率的多个呼吸子之间的相互作用形成[42],即当两个或多个呼吸子碰撞时,会产生时空均局域且幅值较大的怪波[43]。另外,对怪波形成机制的理解非常有用,在海洋中可以有效预防怪波的产生,从而避免造成危害;在光学领域的光纤系统中产生怪波,可以用于产生高能量脉冲。

本文基于标准非线性薛定谔方程,采用达布变换方法推导其一阶和高阶呼吸子的精确表达形式,并详细讨论怪波极限及其他动力学特性。

2 NLSE的一呼吸子解及其怪波极限

标准的无量纲NLSE的表达式为

iψx+122ψt2+ψ2ψ=0,(1)

式中 ψ(x,t)为脉冲包络,在光纤光学中[21],x为沿光纤的传输距离,t为延迟时间。采用达布变换可以将NLSE的各阶呼吸解的递推表达式写成[44]

ψn=ψn-1+2ln*-lnsn1rn1*/rn12+sn12,(2)

式中n代表呼吸子解的阶数,ln为复特征值,rn1sn1为线性变化。

n=1时,一阶呼吸子解为

ψ1=ψ0+2l1*-l1s11r11*/r112+s112,(3)

式中ψ0=exp(ix)为种子解,复特征值l1=a1+ib1,a1为实部,b1为虚部,调制频率与特征值的关系为κj=2 1+lj2=κjr+iκji,下标r和i分别代表参数的实部与虚部,下标j表示n阶方程中第j个参数分量(1≤jn)。线性变化r11s11的表达式见文献[ 44-45]附录。当a1≠0且b1≠0时,一呼吸子解如图1所示。特征值l1的实部a1控制呼吸子解的速度,虚部b1控制呼吸子解的周期,当a1=-0.5,b1=0.8时,呼吸子的传输方向与t轴正向的夹角为锐角,在坐标平面内存在5个基本呼吸子周期,如图1(a)所示;当a1=0.5,b1=1时,呼吸子与t轴正向的夹角为钝角,且呼吸子的周期比图1(a)小,如图1(b)所示,图中色度条表示脉冲幅值 ψ的大小(下文不再赘述)。

图 1. NLSE方程一呼吸子解。(a) a1=-0.5,b1=0.8;(b) a1=0.5,b1=1

Fig. 1. One-breather solution of NLSE. (a) a1=-0.5, b1=0.8; (b) a1=0.5, b1=1

下载图片 查看所有图片

为了在坐标平面内得到时间和空间上都局域的波——Peregrine怪波,NLSE方程的一阶呼吸子解的周期应趋于无穷,即当a1=0,b1→1时,κ1≈0。分析可知,当a1=0时,波峰的山脊方向与x轴平行,即没有速度,如图2所示。图2(a)、图2(b)分别给出了坐标偏移量x1=0,t1=0时中心位置位于(0,0)和x1=3,t1=3时中心位置位于(3,3)处的Peregrine怪波解。

图 2. Peregrine孤子解。(a) x1=0,t1=0;(b) x1=3,t1=3

Fig. 2. Peregrine soliton solution. (a) x1=0, t1=0; (b) x1=3, t1=3

下载图片 查看所有图片

3 NLSE的双呼吸子解及其怪波极限

为了得到NLSE方程的双呼吸子解,需要两个复特征值,即lj=aj+ibj,j=1,2,基于NLSE方程的一阶呼吸子解,当n=2时,由(2)式可得NLSE方程的双呼吸子解为

ψ2=ψ1+2l2*-l2s21r21*/r212+s212(4)

n>1时,NLSE方程的呼吸子解需要更高阶次的rs,此时

rnp=ln-1*-ln-1sn-1,1*rn-1,1sn-1,p+1+lp+n-1-ln-1  rn-1,12rn-1,p+1+lp+n-1-ln-1*sn-1,12rn-1,p+1/rn-1,12+sn-1,12snp=ln-1*-ln-1rn-1,1*sn-1,1rn-1,p+1+lp+n-1-ln-1  sn-1,12sn-1,p+1+lp+n-1-ln-1*rn-1,12sn-1,p+1/rn-1,12+sn-1,12,(5)

式中下标p仅用于列举,并无特定物理意义,当n=2,p=1时,NLSE的双呼吸子传输如图3所示。

图3(a)可知,当l1=0.1+0.8i,l2=-0.1+0.8i,x1=x2=t1=t2=0时,两列呼吸子在平面内传输,并在坐标原点发生碰撞,产生一个幅值较高的二阶怪波,其幅值如色度条所示;碰撞之后,两列呼吸子保持原来的速度、大小和方向向前传播。当bj→1时,即l1=0.1+0.99i,l2=-0.1+0.99i时,碰撞产生的二阶怪波幅值与阶数n之间的关系为 ψ=2n+1,如图3(b)所示。二阶怪波的幅值 ψ=5,对比图3(a)与图3(b)可知,当特征值虚部增大时,二阶怪波的幅值也逐渐增大;当l1=0.1+0.9i,l2=-0.1+0.9i,且x1=t1=-1,x2=t2=1时,两列呼吸子没有发生碰撞,有效避免了二阶怪波的产生,如图3(c)所示。

图 3. 双呼吸子的碰撞。(a) b1=0.8;(b) b1→1;(c)双呼吸子碰撞的避免

Fig. 3. Collision of double breather. (a) b1=0.8; (b) b1→1; (c) avoid the collision of double breather

下载图片 查看所有图片

aj=0,坐标偏移量取不同值时,得到特殊双呼吸子的碰撞与分离,如图4所示。由图4(a)可知,当bj<1时,κj为实数,此时呼吸子表现为纵向局域横向周期呼吸的Akhmediev呼吸子,即当l1=0.95i,l2=0.8i,且x1=x2=t1=t2=0时,两列AB呼吸子在坐标平面内发生碰撞,在坐标平面的中心位置产生一个幅值较高的二阶怪波;当l1l2不变,x1=3,x2=-3时,两列碰撞的AB呼吸子在平面内发生分离,未产生二阶怪波,如图4(b)所示;当bj>1时,κj为虚数,此时呼吸子表现为横向局域纵向周期性呼吸的Kuznetsov-Ma孤子,即当l1=0.9i,l2=1.2i时,κ1为实数,κ2为虚数,AB呼吸子与KM孤子在平面内发生碰撞,产生一个幅值较高的二阶怪波,如图4(c)所示。

图 4. 呼吸子的碰撞与分离。(a)双AB碰撞;(b)双AB分离;(c) AB与KM孤子碰撞

Fig. 4. Collision and separation of breathers. (a) Double Akhmediev breather collision; (b) double Akhmediev breather separation; (c) collision between Akhmediev breather and Kuznetsov-Ma soliton

下载图片 查看所有图片

图 5. 二阶怪波的产生与分裂。(a)二阶怪波;(b)二维图;(c)分裂

Fig. 5. Generation and separation of second-order rouge waves. (a) Second-order rouge waves; (b) two-dimensional diagram; (c) separation

下载图片 查看所有图片

aj=0,bj→1时,坐标偏移量取不同值时,可分别获得Peregrine孤子的叠加和二阶怪波的分离,如图5所示。由图5(a)可知,当l1=0.99i,l2=0.9901i时,(4)式变为二阶有理怪波解,在坐标平面内产生一个具有唯一最大值的二阶怪波;以视线垂直 x,ψ平面方向为三维图的主视图方向,图5(b)为侧视图,则二阶怪波的幅值接近5,即当n=2时, ψ=5;当偏移量不为零,且x1=0.01,x2=-0.01时,一个二阶怪波分裂为三个一阶怪波,如图5(c)所示,且三个怪波的中心位置在x-t平面内呈等边三角形分布。

由于κj与呼吸子周期有关,当κ1κ2同比值不同时,将会得到不同的碰撞模式。令aj=0,图6分别给出了当κ1κ2=2κκ,κ1κ2=3κ∶2κ,κ1κ2=4κ∶3κ时,双呼吸子的几种碰撞模式,其中常数κ=0.4。由图6(a)可知,当κ1κ2=2κκ时,在坐标平面中心位置处沿t轴方向产生一个周期的二阶怪波;当κ1κ2=3κ∶2κ时,在坐标平面的中心位置处产生二阶怪波,且二阶怪波与呼吸子对交替出现,除中心位置的二阶怪波外,其他位置二阶怪波的能量都衰减,如图6(b)所示;当κ1κ2=4κ∶3κ时,二阶怪波和两对呼吸子对形成一个周期,如图6(c)所示。

图 6. 双呼吸子碰撞。(a) κ1∶κ2=2κ∶κ; (b) κ1∶κ2=3κ∶2κ; (c) κ1∶κ2=4κ∶3κ

Fig. 6. Double breather collision. (a) κ1∶κ2=2κ∶κ; (b) κ1∶κ2=3κ∶2κ; (c) κ1∶κ2=4κ∶3κ

下载图片 查看所有图片

图6中可观察到一个有趣的现象:在每个周期序列中,一阶怪波分量按照特定规则进行组织。首先,调制频率κ1κ2比值的整数定义了每个周期中一阶怪波的总数,即κ1κ2=2κκ时,怪波的总数为3;其次,三个一阶怪波分量总是合并形成如图5(a)所示的一个二阶怪波,其余怪波关于x=0对称。在图6(b)中,κ1κ2=3κ∶2κ,即在一个周期中,一阶怪波的总数为5,其中三个一阶怪波组合成一个二阶怪波,剩余两个一阶怪波关于x=0对称。

双呼吸子的碰撞特性具有极强的丰富性,本文只讨论几种具有代表性的碰撞关系。已经讨论了一般的双呼吸子的碰撞,接下来讨论双呼吸子碰撞的另一种特殊情况,即两个呼吸子并行传播时的碰撞特性[46]。在并行传输模式下,双呼吸子传播需要满足条件:

VH1/κ1i=VH2/κ2i(6)

根据文献[ 47]中对各参数的定义可以得到:

VH1=a1κ1i-b1κ1r/2VH2=a2κ2i-b2κ2r/2,(7)

式中VH1控制脉冲在x方向上的位置,且VHj/κji与NLSE方程中孤子解的群速度相对应,此时(6)式没有解析解,但可以计算得到数值解。假设l2已知,令l2=-0.5+0.4i,当(6)式成立时,可解得l1=-1.1705+1.4i,此时双呼吸子的并行传输如图7所示。从图7(a)可以看出,当x1=0,x2=0时,两列呼吸子周期性并行传输,且在传输过程中碰撞叠加,图7(b)为图7(a)的局部放大图;当x1=-2,x2=2时,并行传输的呼吸子序列分离为两列平行传输的呼吸子序列,并保持各自的周期并行传输,在坐标平面内不会发生碰撞叠加。

图 7. l1=-1.1705+1.4i,l2=-0.5+0.4i时双呼吸子的并行叠加。(a) x1=0,x2=0;(b)局部放大图;(c) x1=-2,x2=2

Fig. 7. Two breathers parallel superimposed with l1=-1.1705+1.4i and l2=-0.5+0.4i. (a) x1=0, x2=0; (b) partial enlargement; (c) x1=-2, x2=2

下载图片 查看所有图片

4 NLSE方程的三阶呼吸子解及其怪波极限

由(2)式可知,当n=3时,NLSE方程的三阶呼吸子解为

ψ3=ψ2+2l3*-l3s31r31*/r312+s312,(8)

式中lj=aj+ibj(j=1,2,3)为三个复特征值,参数s31r31由(5)式给出。当特征值分别取不同值时,可以得到呼吸子的双重碰撞、三重碰撞和简并态等,如图8所示。从图8(a)中可以看出,当l1=0.5+0.98i,x1=0,t1=10,l2=1.2i,x2=0,t2=-10,l3=-0.5+0.98i,x3=0,t3=10时,三阶呼吸子解在坐标平面内发生两两相互碰撞,且碰撞点分布呈三角形,当选取不同的特征值时,双重碰撞所呈现的图形也不相同;当l1=0.5+0.95i,l2=0.95i,l3=-0.5+0.95i,且xj=tj=0时,三阶呼吸子在坐标平面内相交于一点,在碰撞点处产生一个幅值较大的三阶怪波,如图8(b)所示;仔细比较图8(a)与(b)可知,合理选取参数可以有效避免高阶怪波的产生,一般地,各特征值不能相等,否则解析解将没有意义,为了克服该难题,可以考虑ljl的情况,即特征值近似相等的情况,令l=0.95i,特征值分别为l1=0.9501i,l2=0.9502i,l3=0.9503i时,三阶呼吸子的传输如图8(c)所示,三阶呼吸子解表现为简并模式,三列呼吸子先保持近似平行的向前传播,在坐标原点处发生碰撞叠加,产生三阶怪波,碰撞之后,三列呼吸子继续维持各自的速度、振幅和方向向前传播,且传播过程中方向仍保持近似平行。

aj=0,bj→1时,κj≈0,呼吸子的周期趋于无穷,三阶呼吸子解转变为三阶有理怪波解,如图9所示。由图9(a)可知,当l1=0.99i,l2=0.9901i,l3=0.9902i,且xj=tj=0时,在坐标原点处产生一个孤立的三阶怪波,幅值接近7,即幅值与阶数仍满足 ψ=2n+1的关系;图9(b)为图9(a)的等高图,颜色相同表示次峰的幅值相同,曲线从外到内颜色越明亮表示脉冲幅值越大,若从峰值功率最大点出发垂直于x轴做一条分界线,将三阶怪波分为左右两部分,假定每一部分最外层的次峰为一级次峰,级数由外向内逐渐递增,如左侧最外层的一级次峰为三瓣结构,而内层的二级次峰为双瓣结构,此时三阶怪波每一部分的次峰个数N均为5,即N=[n(n+1)/2]-1;当lj保持不变,坐标偏移不全为零,如x1=x3=10-5,x2=0时,此时三阶怪波发生分裂,分裂后的三阶怪波在坐标平面内呈五边形分布,且位于中心位置处怪波的阶数为n-2,即中心位置处为一阶Peregrine孤子,而包围中心位置处怪波的旁峰数目M与阶数n的关系为M=2n-1,旁峰的个数为5,如图9(c)所示。

同样地,当κ1κ2κ3比值不同时,将会得到不同的三阶呼吸子的叠加模式。为简单起见,令aj=0,图10给出了κ1κ2κ3满足不同比值时,三阶呼吸子的相互作用。

图 8. 三阶呼吸子的碰撞。(a)双重碰撞;(b)三重碰撞;(c)简并态

Fig. 8. Collisions of third-order breathers. (a) Double collision; (b) triple collision; (c) degeneration

下载图片 查看所有图片

图 9. 三阶怪波的产生与分裂。(a)三阶怪波;(b)等高图;(c)分裂

Fig. 9. Generation and separation of third-order rouge waves. (a) Third-order rouge waves; (b) contours; (c) separation

下载图片 查看所有图片

图 10. 三阶呼吸子的碰撞。(a) κ1∶κ2∶κ3=2∶3∶5且κ1=0.8;(b) κ1∶κ2∶κ3=2∶3∶5且κ1=0.6;(c) κ1∶κ2∶κ3=3∶2∶1且κ1=1.8

Fig. 10. Third-order breather collision. (a) κ1∶κ2∶κ3=2∶3∶5 and κ1=0.8; (b) κ1∶κ2∶κ3=2∶3∶5 and κ1=0.6; (c) κ1∶κ2∶κ3=3∶2∶1 and κ1=1.8

下载图片 查看所有图片

图10(a)可知,当κ1κ2κ3=2∶3∶5,且κ1=0.8时,三阶呼吸子相互叠加,在坐标平面内呈现怪波与呼吸子交替出现的周期性结构;当κ1κ2κ3=2∶3∶5,但κ1=0.6时,三阶呼吸子的叠加模式发生改变,即在坐标原点处产生一个高阶怪波,碰撞前后三个呼吸子均保持各自的传播特性,不具有周期性,如图10(b)所示;当κ1κ2κ3=3∶2∶1且κ1=1.8时,此时三阶呼吸子在x=0处产生怪波序列,怪波的幅度变化趋势呈现抛物线形,在t=0时峰值达到最大,如图10(c)所示。

由二阶呼吸子解及怪波极限的分析可知,三阶呼吸子并行传输需要满足条件:

VH1/κ1i=VH2/κ2i=VH3/κ3i(9)

假设l2已知,令l2=0.3+0.7i,当(9)式成立时,其数值解为l1=0.1935+0.7982i,l3=1.0046+1.2018i,此时三阶呼吸子的传输特性如图11所示。由图11(a)可知,当三阶呼吸子特征值满足(9)式,且偏移量均为零时,传输模式为并行传输模式,三列呼吸子先沿同一方向向前传输,传输一定距离后,在坐标平面内发生碰撞并形成三阶怪波,在坐标原点处三阶怪波的峰值功率大于其他碰撞位置的峰值功率,此时三阶呼吸子的并行传输具有周期性;图11(b)为图11(a)的局部放大图;当三列呼吸子的坐标偏移量不为零时,三列呼吸子的并行传输发生分离,如图11(c)所示,分离后三列呼吸子的传输方向平行,在传输过程中不会发生碰撞或叠加。

图 11. 三阶呼吸子的并行传输。(a) x1=0,x2=0,x3=0;(b)局部放大图;(c) x1=-4,x2=0,x3=4

Fig. 11. Third-order breather parallel transmission. (a) x1=0, x2=0, x3=0; (b) partial enlargement; (c) x1=-4, x2=0, x3=4

下载图片 查看所有图片

5 NLSE的四阶呼吸子解及其怪波极限

由(2)式可知,当n=4时,NLSE方程的四阶呼吸子解为

ψ4=ψ3+2l4*-l4s41r41*/r412+s412,(10)

式中s41r41由(5)式可得。当特征值取不同值时,四阶呼吸子解的传输特性如图12所示。由图12(a)可知,当特征值分别为l1=0.95i,l2=0.9501i,l3=1.05i,l4=1.0501i,x1=-2,x2=2,t3=-4,t4=4时,四阶呼吸子呈井字型分布,且在传输过程中没有发生碰撞叠加,即未产生高阶怪波;由图12(b)可知,当l1=0.9i,l2=1.2i,l3=-0.5+0.95i,l4=0.5+0.95i,且坐标偏移量均为零时,呼吸子在坐标原点处发生四重碰撞,产生四阶怪波,碰撞后,每列呼吸子仍保持其原有的速度、方向和大小向前传播,呼吸子呈米字型分布。

图 12. 四阶呼吸子的相互作用。(a)双重碰撞;(b)四重碰撞

Fig. 12. Fourth-order breather interaction. (a) Double collision; (b) fourfold collision

下载图片 查看所有图片

特别地,当l1=0.95i,l2=0.9501i时,特征频率κ1κ2为实数,且xj=tj=0时,可以获得AB呼吸子的简并模式,两列AB呼吸子在坐标原点处发生碰撞,产生一个二阶怪波;当l3=1.05i,l4=1.0501i时,特征频率κ3κ4为虚数,可以获得双KM孤子的简并模式,两列KM孤子在坐标原点处发生碰撞,产生一个二阶怪波。当双AB呼吸子的简并模式与双KM孤子的简并模式在坐标原点位置发生碰撞时,会产生一个四阶怪波,其幅值为9,如图13(a)所示;当l1=0.95i,l2=0.9501i,l3=0.9502i,l4=0.9503i,且xj=tj=0时,呼吸子的传输模式为四阶AB呼吸子的简并模式,即四列AB呼吸子以近似平行的方式传播,并在坐标原点碰撞产生一个四阶怪波,其幅值小于图13(a)中的幅值,碰撞后仍近似平行向前传播,如图13(b)所示。

aj=0,bj→1,κj→0时,呼吸子的周期趋于无穷,四阶呼吸子解转变为四阶怪波解。当l1=0.99i,l2=0.991i,l3=0.992i,l4=0.993i,且xj=tj=0时,在坐标平面内形成四阶怪波,如图14(a)所示;四阶怪波的幅度 ψ=9,与阶数n的关系满足 ψ=2n+1,如图14(b)所示;图14(c)为图14(a)的等高图,分界线左侧一级次峰为四瓣结构,二级次峰为三瓣结构,三级次峰为双瓣结构,即分界线左侧的次峰总数N=9,与阶数n的关系满足N=[n(n+1)/2]-1。

当四阶怪波解具有较小的坐标偏移时,四阶怪波将发生分裂,为了方便分析,将特征频率以及坐标偏移量分别改写为

κj=,xj=Xjnκ2(n-1),tj=Tjnκ2(n-1),(11)

式中参数κXjnTjn均为常数,适当取值使得坐标漂移量分别为x1=0.56,x2=x3=x4=0时,四阶怪波的分裂情况如图15(a)所示。可以看出,四阶怪波在一次分裂后产生一个圆形的怪波集群,其中心为一个二阶怪波,即中心怪波的阶数为n-2,分裂后外层旁峰的数目M与阶数n满足M=2n-1,即M=7(图中可见7个旁峰)。特别地,当x1=x2=0.56,x3=x4=0时,四阶怪波发生两次分裂,即在一次分裂的基础上,中心二阶怪波再次分裂,产生三个一阶Peregrine怪波,四阶怪波最多只能进行两次分裂,如图15(b)所示。

图 13. 呼吸子的简并。(a) AB简并与KM孤子简并叠加;(b) AB的简并

Fig. 13. Degenerate breather. (a) Superposition of Akhmediev breather and Kuznetsov-Ma soliton degenerate state; (b) degenerate Akhmediev breather

下载图片 查看所有图片

图 14. 四阶怪波的产生。(a)四阶怪波;(b)二维图;(c)等高图

Fig. 14. Generation of fourth-order rouge wave. (a) Fourth-order rouge waves; (b) two-dimensional diagram; (c) contour map

下载图片 查看所有图片

图 15. 四阶怪波的分裂。(a)一次分裂;(b)两次分裂

Fig. 15. Separation of fourth-order rouge wave. (a) Once separation; (b) twice separation

下载图片 查看所有图片

同样地,由于κj对呼吸子周期具有较大影响,当κj满足不同比值时,将会得到不同的四阶呼吸子叠加模式。为简单起见,令aj=0。图16给出了κ1κ2κ3κ4满足不同比值时,四阶呼吸子的相互作用。当κ1κ2κ3κ4=1∶3∶5∶7且κ1=0.2和κ1κ2κ3κ4=2∶4∶6∶8且κ1=0.4时,四阶呼吸子在中心位置处会叠加产生四阶怪波,其他位置表现为完全不同的呼吸子叠加行为,但脉冲的分布均关于t=0和x=0对称,等高图中曲线颜色越明亮说明脉冲幅值越大。

四阶呼吸子并行传输需要满足条件:

VH1/κ1i=VH2/κ2i=VH3/κ3i=VH4/κ4i(12)

同样地,假设l2已知,令l2=0.3+0.7i,当(12)式成立时,其数值解为l1=0.1935+0.7982i,l3=1.0046+1.2018i,l4=0.5991+0.5468i,此时四阶呼吸子的并行传输如图17所示。

图17(a)可知,当四阶呼吸子的特征值满足(12)式,且坐标偏移均为零时,四阶呼吸子呈周期性并行传输。但在图17(a)中,只能观察到三列明显的呼吸子序列,这是因为其中两列呼吸子传输方向完全重合,形成一条传输方向不变的中心呼吸子序列,该序列幅度高于两侧呼吸子序列。在坐标原点,四阶呼吸子相互碰撞产生四阶怪波,图17(b)为图17(a)的局部放大图;当坐标偏移量不为零时,四阶呼吸子的并行传输发生分离,表现为四列相互平行的呼吸子序列,呼吸子在传输过程中不会产生叠加,如图17(c)所示。

图 16. 四阶呼吸子相互作用的等高图。(a) κ1∶κ2∶κ3∶κ4=1∶3∶5∶7;(b) κ1∶κ2∶κ3∶κ4=2∶4∶6∶8

Fig. 16. Contour map of fourth-order breather interaction. (a) κ1∶κ2∶κ3∶κ4=1∶3∶5∶7; (b) κ1∶κ2∶κ3∶κ4=2∶4∶6∶8

下载图片 查看所有图片

图 17. 四阶呼吸子的并行传输。(a) x1=0,x2=0,x3=0,x4=0;(b)局部放大图;(c) x1=-6,x2=-2,x3=2,x4=6

Fig. 17. Third-order breather parallel transmission. (a) x1=0, x2=0, x3=0, x4=0; (b) partial enlargement; (c) x1=-6, x2=-2, x3=2, x4=6

下载图片 查看所有图片

6 结论

基于标准无量纲的非线性薛定谔方程,采用达布变换方法得到了NLSE方程各阶呼吸子解的一般表达形式,并详细研究了怪波极限及其他动力学特性。当呼吸子在坐标平面内发生有效碰撞或周期趋于无穷时,会产生一个对应阶数的怪波。

然而,高能量脉冲在光纤光学中有用,海洋中高能量脉冲具有极强的破坏力,应当有效预防和避免。研究结果表明,合理设定参数既可以有效控制高阶怪波的产生,还可以选择怪波出现的位置。此外,还研究了高阶呼吸子的并行传输模式和频率比对传输特性的影响,以及高阶怪波的特性和一次分裂、二次分裂。本文只研究了标准的NLSE方程,后续工作将针对其他NLSE方程进行,期望详尽了解呼吸子以及怪波的动力学特性。

参考文献

[1] Zakharov V E, Shabat A B. Exact theory of two-dimensional self-focussing and one-dimensional self-modulating waves in nonlinear media[J]. Journal of Mathematical Physics, 2015, 34(15): 62-69.

[2] AkhmedievN, AnkiewiczA. Solitons: nonlinear pulses and beams[M]. London: Chapman and Hall, 1997.

[3] Kharif C, Pelinovsky E. Physical mechanisms of the rogue wave phenomenon[J]. European Journal of Mechanics, 2003, 22(6): 603-634.

[4] . Nonlinear four-wave interactions and freak waves[J]. Journal of Physical Oceanography, 2003, 33(4): 863-884.

[5] Onorato M, Osborne A R, Serio M, et al. Freak waves in random oceanic sea states[J]. Physical Review Letters, 2001, 86(25): 5831-5834.

[6] Dyachenko A I, Zakharov V E. Spatial equation for water waves[J]. JETP Letters, 2016, 103(3): 181-184.

[7] Shrira V I, Geogjaev V V. What makes the Peregrine soliton so special as a prototype of freak waves?[J]. Journal of Engineering Mathematics, 2010, 67(1/2): 11-22.

[8] Osborne A R, Onorato M, Serio M. The nonlinear dynamics of rogue waves and holes in deep-water gravity wave trains[J]. Physics Letters A, 2000, 275(5/6): 386-393.

[9] Osborne A R. The random and deterministic dynamics of ‘rogue waves’ in unidirectional, deep-water wave trains[J]. Marine Structures, 2001, 14(3): 275-293.

[10] PerrieW. Nonlinear ocean waves[M]. Southampton: Computational Mechanics Publications, 2006.

[11] Hasegawa A, Tappert F. Transmission of stationary nonlinear optical pulses in dispersive dielectric fibers. I. Anomalous dispersion[J]. Applied Physics Letters, 1973, 23(3): 142-144.

[12] Agrawal GP. Nonlinear fiber optics[M]. 4 th ed. Amsterdam: Elsevier , 2006.

[13] Andreev P A. First principles derivation of NLS equation for BEC with cubic and quintic nonlinearities at nonzero temperature: dispersion of linear waves[J]. International Journal of Modern Physics B, 2013, 27(6): 1350017.

[14] Slunyaev A V. A high-order nonlinear envelope equation for gravity waves in finite-depth water[J]. Journal of Experimental and Theoretical Physics, 2005, 101(5): 926-941.

[15] 武达, 王娟芬, 石佳, 等. 掺杂光纤中Peregrine孤子的产生和传输[J]. 光学学报, 2017, 37(4): 0406002.

    Wu D, Wang J F, Shi J, et al. Generation and transmission of Peregrine solitons in doped fiber[J]. Acta Optica Sinica, 2017, 37(4): 0406002.

[16] Vishnu P N, Senthilvelan M, Lakshmanan M. Akhmediev breathers, Ma solitons, and general breathers from rogue waves: a case study in the Manakov system[J]. Physical Review E, 2013, 88(2): 022918.

[17] Frisquet B, Kibler B, Millot G. Collision of akhmediev breathers in nonlinear fiber optics[J]. Physical Review X, 2013, 3(4): 041032.

[18] Kedziora D J, Ankiewicz A, Akhmediev N. Second-order nonlinear Schrödinger equation breather solutions in the degenerate and rogue wave limits[J]. Physical Review E, 2012, 85(6): 066601.

[19] 杨光晔, 李禄, 田晋平. 基于谱过滤方法的Kuznetsov-Ma孤子向准基态孤子转化研究[J]. 光学学报, 2016, 36(6): 0619002.

    Yang G Y, Li L, Tian J P. Study on Transformation of Kuznetsov-Ma soliton to quasi-fundamental soliton based on spectral-filtering method[J]. Acta Optica Sinica, 2016, 36(6): 0619002.

[20] Bergé L. Wave collapse in physics: principles and applications to light and plasma waves[J]. Physics Reports, 1998, 303(5/6): 259-370.

[21] Efimov V B, Ganshin A N, Kolmakov G V, et al. Rogue waves in superfluid helium[J]. The European Physical Journal Special Topics, 2010, 185(1): 181-193.

[22] Shats M, Punzmann H, Xia H. Capillary rogue waves[J]. Physical Review Letters, 2010, 104(10): 104503.

[23] Kibler B, Fatome J, Finot C, et al. The Peregrine soliton in nonlinear fibre optics[J]. Nature Physics, 2010, 6(10): 790-795.

[24] Kichenassamy S. Breather solutions of the nonlinear wave equation[J]. Communications on Pure and Applied Mathematics, 1991, 44(7): 789-818.

[25] Tajiri M, Watanabe Y. Breather solutions to the focusing nonlinear Schrödinger equation[J]. Physical Review E, 1998, 57(3): 3510-3519.

[26] Kovalyov M. Modulating properties of harmonic breather solutions of KdV[J]. Journal of Physics A: Mathematical and General, 1998, 31(22): 5117-5128.

[27] Wang C J, Dai Z D, Lin S Q, et al. Breather-type soliton and two-soliton solutions for modified Korteweg-de Vries equation[J]. Applied Mathematics and Computation, 2010, 216(1): 341-343.

[28] Sarkar R, Dey B. Exact compact breather-like solutions of two-dimensional Fermi-Pasta-Ulam lattice[J]. Journal of Physics A: Mathematical and General, 2006, 39(4): L99-L104.

[29] 朱叶青, 胡巍. 无外置偏压的向列相液晶盒中的呼吸子传输[J]. 光学学报, 2015, 35(9): 0919001.

    Zhu Y Q, Hu W. Propagation of breathers in the nematic liquid crystal cell without bias voltage[J]. Acta Optica Sinica, 2015, 35(9): 0919001.

[30] Benjamin T B, Feir J E. The disintegration of wave trains on deep water. Part 1. Theory[J]. Journal of Fluid Mechanics, 1967, 27(3): 417-430.

[31] 卢洵, 王东升. 色散管理孤子系统的调制不稳定性[J]. 激光技术, 2012, 36(4): 557-561.

    Lu X, Wang D S. Modulation instability in dispersion managed soliton systems[J]. Laser Technology, 2012, 36(4): 557-561.

[32] Akhmediev N N, Korneev V I, Mitskevich N V. Modulation instability in an optical fiber induced by cross-phase modulation[J]. Radiophysics and Quantum Electronics, 1991, 34(1): 73-77.

[33] Akhmediev N. Soto-Crespo J M, Ankiewicz A. Extreme waves that appear from nowhere: on the nature of rogue waves[J]. Physics Letters A, 2009, 373(25): 2137-2145.

[34] Akhmediev N, Ankiewicz A, Taki M. Waves that appear from nowhere and disappear without a trace[J]. Physics Letters A, 2009, 373(6): 675-678.

[35] Shukla P K, Kourakis I, Eliasson B, et al. Instability and evolution of nonlinearly interacting water waves[J]. Physical Review Letters, 2006, 97(9): 094501.

[36] Tao Y S, He J S. Multisolitons, breathers, and rogue waves for the Hirota equation generated by the Darboux transformation[J]. Physical Review E, 2012, 85(2): 026601.

[37] He J S, Xu S W, Porsezian K. Rogue waves of the Fokas-Lenells equation[J]. Journal of the Physical Society of Japan, 2012, 81(12): 124007.

[38] Li C, He J, Porsezian K, et al. Rogue waves of the Hirota and the Maxwell-Bloch equations[J]. Physical Review E, 2013, 87(1): 012913.

[39] Bandelow U, Akhmediev N. Sasa-Satsuma equation: soliton on a background and its limiting cases[J]. Physical Review E, 2012, 86(2): 026606.

[40] Soto-Crespo J M, Grelu P, Akhmediev N. Dissipative rogue waves: extreme pulses generated by passively mode-locked lasers[J]. Physical Review E, 2011, 84: 016604.

[41] Bespalov V I, Talanov V I. Filamentary structure of light beams in nonlinear liquids[J]. Zhetf Pisma Redaktsiiu, 1966, 3(12): 471.

[42] Voronovich V V, Shrira V I, Thomas G. Can bottom friction suppress 'freak wave' formation?[J]. Journal of Fluid Mechanics, 2008, 604: 263-296.

[43] He J S, Zhang H R, Wang L H, et al. Generating mechanism for higher-order rogue waves[J]. Physical Review E, 2013, 87(5): 052914.

[44] Kedziora D J, Ankiewicz A, Akhmediev N. Circular rogue wave clusters[J]. Physical Review E, 2011, 84(5): 056611.

[45] Kedziora D J, Ankiewicz A, Akhmediev N. Classifying the hierarchy of nonlinear-Schrödinger-equation rogue-wave solutions[J]. Physical Review E, 2013, 88: 013207.

[46] Chowdury A, Kedziora D J, Ankiewicz A, et al. Breather solutions of the integrable quintic nonlinear Schrödinger equation and their interactions[J]. Physical Review E, 2015, 91(2): 022919.

[47] Chowdury A, Krolikowski W. Breather-to-soliton transformation rules in the hierarchy of nonlinear Schrödinger equations[J]. Physical Review E, 2017, 95(6): 062226.

杜志峰, 宋丽军, 王艳. 非线性薛定谔方程的呼吸子解及其怪波极限[J]. 激光与光电子学进展, 2019, 56(5): 050601. Zhifeng Du, Lijun Song, Yan Wang. Breather Solutions and Their Rouge Wave Limits of NonlinearSchrödinger Equation[J]. Laser & Optoelectronics Progress, 2019, 56(5): 050601.

引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!