买欣蕾 李 威

(北京化工大学 数理学院， 北京 100029)

引 言

非线性薛定谔方程(NLSE)是一类重要的非线性演化方程，在非线性光学、量子力学、等离子体理论等领域具有广泛的应用，因此研究该类方程的精确解具有重要的理论意义。近年来，人们提出许多求解NLSE精确解的有效方法，如齐次平衡法[1]、tanh函数展开法及其扩展[2]、正弦- 余弦法[2]、exp-函数法[2]、试验函数法及其扩展[3]、扩展辅助方程法[4]、改进的Fan-子方程法[5]、(G′/G)-展开法[6]、exp(-φ(ξ))-展开法[6]等。众所周知，经典的NLSE描述了脉冲宽度在100 fs以上的光脉冲在光纤中的传播[7]。随着脉冲宽度逐渐变窄，如接近50 fs甚至低于10 fs时，系统必须考虑非线性项和高阶色散(如三阶和四阶色散)。

本文考虑描述超短光脉冲在光纤中传播的具有高阶色散和立方- 五次非线性项的NLSE[8]

(1)

利用动力系统的分岔理论不仅能够得到方程局部解的表达式，而且还可以探查方程解的整体结构[10-11]，因此本文先利用代入法将方程(1)转化为动力系统，再运用分岔理论来分析系统的动力特性。通过分析相图的有界轨道不仅得到了方程有界行波解的显式表达式，而且还可直接获知行波解的分类结构。

1 动力系统分岔分析

对方程(1)采用行波变换E(z,t)=φ(ξ)eiη，其中ξ=v0z-vt,η=ω0z-ωt，并分离实部、虚部可得

γ1φ3(ξ)-γ2φ5(ξ)=0

(2)

式中，

当l1=l3=0时文献[5]采用改进的Fan-子方程法获得对应的行波解；当θ=l2l3-l1l4≠0时将式(2)一式微分并代入式(2)二式中得

φ″(ξ)=c2φ(ξ)+2c4φ3(ξ)+3c6φ5(ξ)

(3)

式中，c2=-l0l3/θ,2c4=-γ1l3/θ,3c6=γ2l3/θ。式(3)等价于如下所示的二维自治系统

(4)

这是一个具有Hamilton函数

H(φ,y)=y2-c2φ2-c4φ4-c6φ6=h

图1 平衡点个数及其类型的分布区域Fig.1 The distribution regions of the number and type of equilibrium points

图2 不同区域所对应的相图Fig.2 Phase diagrams for different regions

2 NLSE的精确解

对于一个固定的h∈R，曲线

Ch={(φ,y)∈R×R:H(φ,y)=h}

被称为具有h能量水平的能量曲线[12]。显然，H(φ,y)=h的每个轨道都是一条能量曲线，因此可以研究其有界轨道与能量水平h之间的关系。设

Fh(φ)=h+c2φ2+c4φ4+c6φ6

显然能量曲线Ch相当于y2=Fh(φ)所定义的曲线，标记H(Q±)=h1=-F0(φ1±)，H(R±)=h2=-F0(φ2±)。下面将通过分析情况③～④所对应能量h的变化来得到系统精确解的分类及分布。

2.1 区域Ⅴ对应系统的精确解

根据能量水平h范围的不同，分为如下3种情况。

1)h=0时，这是一条连接鞍点O的稳定流形和不稳定流形的同宿轨

E1±(z,t)=

(5)

其分别对应于左半平面的同宿轨道族L-(O,O)和右半平面的同宿轨道族L+(O,O)，波形图如图4(a)所示。

2)h∈(h1,0)时，这是位于同宿轨内部，分别围绕中心点Q±的两个周期轨

(6)

(7)

(8)

图时能量曲线与相图的对应关系Fig.3 Correspondence between the energy curve and phase diagram at

图4 孤立波解E1±和周期波解E2±、E3±、E4±的三维波形图Fig.4 The three-dimensional waveform of the solitary wave solution E1± and periodic wave solutions E2±,E3±,E4±

2.2 区域Ⅵ对应系统的精确解

图时能量曲线与相图对应关系Fig.5 Correspondence between the energy curve and phase diagram at

根据能量水平h范围的不同，分为如下两种情况。

1)h=h2时，这是一条连接鞍点R±的稳定流形和不稳定流形的异宿轨

E5±(z,t)=

(9)

分别对应于异宿轨道族L-(R+,R-)和L+(R-,R+)，波形图如图6(a)所示。

2)h∈(0,h2)时，这是位于异宿轨内部，包围中心点O的一个周期轨

(10)

图6 孤立波解E5±和周期波解E6±的三维波形图Fig.6 The three-dimensional waveforms of the solitary wave solution E5±and periodic wave solution E6±

2.3 区域Ⅶ中情况ⓐ对应系统的精确解

图时能量曲线与相图对应关系Fig.7 Correspondence between the energy curve and phase diagram at

根据能量水平h范围的不同，分为如下两种情况。

1)h=h2=0时，这是一条连接鞍点O和R±的稳定流形和不稳定流形的异宿轨

因此得到扭结波解为

(11)

分别对应于异宿轨道族L(O→R+→O)和L(O→R-→O)，波形图如图8所示。

图8 扭结波解E7±的三维波形图Fig.8 The three-dimensional waveform of the kink wave solution E7±

2)h∈(h1,0)时，这是位于同宿轨内部分别围绕中心点Q±的两个周期轨。因此周期波解为

(12)

本文对图2中区域Ⅶ的情况ⓑ、ⓒ及区域Ⅷ分别对应的3个相图也进行了与前面类似的讨论，得到对应的周期波解和孤立波解，并且发现其对应的波形无太大差别，在此不再赘述。

3 结论

本文研究了描述超短光脉冲在光纤中传播的具有高阶色散和立方- 五次非线性项的薛定谔方程的精确解。采用代入法和动力系统的分岔理论得到方程的10个相图，对其中比较典型的3个相图进行理论分析，找到了系统的孤立波解、扭结波解和周期波解等精确解，并进行了数值验证，给出了波形的三维图。孤立波解的存在意味着非线性和色散之间的完美平衡，我们给出了孤立波解的存在条件，从而也证实了该系统所描述的光纤在一定条件下能够稳定地进行长时传播。同时本文对系统动态分析的结果体现了行波结构的多样性和完整性。