李 倩, 舒 級(jí), 楊 袁, 王云肖, 汪春江

(四川師范大學(xué) 數(shù)學(xué)與軟件科學(xué)學(xué)院, 四川 成都 610066)

眾所周知,非線性偏微分方程是用來描述非線性科學(xué)問題的重要模型,而演化方程的孤子解的多樣性反應(yīng)出物理世界各種形式的時(shí)空結(jié)構(gòu),因此求解非線性偏微分方程具有非常重要的理論和應(yīng)用價(jià)值.非線性演化方程在某些特殊情況下才能得到其顯示表達(dá)式.這些年來,科學(xué)家們發(fā)現(xiàn)了許多求解非線性偏微分方程的方法,如Painleve分析法[1]、tanh函數(shù)法[2]、齊次平衡法[3]、Hirota方法[4]、不變子空間法[5]、Wronskin技巧[6]、達(dá)布變換方法[7-12]、貝克隆變換法[13]和反散射方法[14]等.

研究發(fā)現(xiàn),怪波現(xiàn)象可以用非線性Schr?dinger方程的解來描述,并且這一性質(zhì)已被非線性光纖中的怪波實(shí)驗(yàn)證實(shí)[15-16].怪波最初是描述海洋上出現(xiàn)的一種奇怪的水波,以其出現(xiàn)的突然性和異常的波高得名.這里的突然性是指,它出現(xiàn)時(shí)無任何征兆,而后又很快地消失.歷史上,有記載的怪波事件已有很多,比如它對(duì)在海洋航行的各類船只、海上油井等的致命性破壞.因?yàn)榻?jīng)歷過怪波災(zāi)難的人們很少能有機(jī)會(huì)生還,長期以來,大家都以為是海怪造成了這些災(zāi)難,直到1995年初,人們?cè)诒焙Ｖ苯犹綔y(cè)到了被稱為“新年波”的怪波,才使得大家相信這是一種海洋現(xiàn)象而不是所謂的海怪所為[17].由于玻色凝聚體的動(dòng)力學(xué)方程類似于非線性光纖中的動(dòng)力學(xué)方程,可以通過研究玻色體系中的怪波動(dòng)力學(xué)來獲取對(duì)怪波的一般性認(rèn)識(shí)[18-19].研究表明,Schr?dinger方程具有很多非線性波解,包括亮暗孤子解、呼吸子解以及有理形式解,它們可以用來描述諸多豐富的物理現(xiàn)象,比如玻色凝聚體中的孤子性質(zhì)和非線性光纖、波導(dǎo)管等中的孤子傳輸可以用它的孤子解來描述[20-21].

本文研究具有阻尼項(xiàng)的Gross-Pitaevskii方程

iqt+qxx+2|q|2q-

(αx-β2x2)q+iβ|q|2q=0,

(1)

其中,α、β是常數(shù),q(t,x)是復(fù)值波函數(shù).借鑒文獻(xiàn)[22-24]的方法,應(yīng)用達(dá)布變換研究方程(1)的孤子解.

1 Gross-Pitaevskii方程的Lax對(duì)和達(dá)布變換

將構(gòu)造方程(1)的Lax對(duì)和達(dá)布變換.由于方程(1)是可積的,可以用AKNS方法構(gòu)造其Lax對(duì)

φx=Mφ, φt=Nφ,

(2)

其中

這里*表示復(fù)共軛,λ是復(fù)值譜參數(shù),

A=-2iλ2+2iλβx+i|Q|2-

B=2λQ+iQx-2βxQ,

D=2iλ2-2iλβx-i|Q|2+

下面考慮上述譜問題的一個(gè)規(guī)范變換

φ[1]=Tφ,

(3)

將線性問題(2)轉(zhuǎn)化為

(4)

并且M、N與M[1]、N[1]具有相同的形式,除了將M、N中的Q、Q*換成M[1]、N[1]中的Q[1]、Q[1]*.將(4)式代入(3)式中,知道T滿足

M[1]T=Tx+TM,

(5)

N[1]T=Tt+TN,

(6)

其中T是λ的多項(xiàng)式形式變換,即

其中a1、b1、c1、d1、a、b、c、d是關(guān)于x和t的實(shí)函數(shù).由(5)和(6)式知

對(duì)比上式λk(k=0,1,2)的系數(shù),得到:

當(dāng)k=2時(shí),

b1=c1=0.

(7)

當(dāng)k=1時(shí),

a1x=d1x=0,

ax=Q[1]c+Q*b,

bx=Q[1]d+Qa,

-2ib+Q[1]d1-Qa1=0,

2ic+Q[1]*a1-Q*d1=0.

(8)

當(dāng)k=0時(shí),

cx=-Q[1]*a+Q*d,

dx=-Q[1]*b+Qc.

(9)

由(8)式的第一式知a1、d1是常數(shù),不失一般性,取a1=d1=1,因此方程(3)的達(dá)布變換可以寫成下列形式[25]

φ[1]=Tφ=(λI-S)φ,

(10)

其中,λ是復(fù)的譜參數(shù),I是2×2的單位矩陣,S是非奇異矩陣.

將M、M[1]和T代入(5)式,比較譜參數(shù)λ的系數(shù)可得

其中

經(jīng)過一次達(dá)布變換[26],新的特征函數(shù)與原來的特征函數(shù)有如下關(guān)系:

Q[1]=Q-2is12,

-Q[1]*=-Q*+2is21,

(11)

且滿足限制條件

(12)

為了得到矩陣S的表達(dá)式,可以通過Lax對(duì)的解來定義矩陣

S=HΛH-1,

(13)

而

又由(13)式可知

其中Δ=|f1|2+|f2|2,這時(shí)可以驗(yàn)證矩陣S滿足限制條件(12)式.

由(11)和(12)式可以得到方程的一次達(dá)布變換

(14)

(15)

其中

最后,方程的n次達(dá)布變換的行列式為

(16)

其中

2 方程的孤子解

(17)

(18)

(19)

(20)

可以解得

(21)

(22)

其中

Q[1]=

顯然得到方程的解為

q[1]=

當(dāng)α=1、β=1、n=1和c=-0.5i時(shí),孤子解q[1]如圖1.

圖 1 孤子解

討論一類具非線性阻尼項(xiàng)的GP方程的達(dá)布變換和孤子解,該方程在玻色-愛因斯坦凝聚中有重要意義.首先通過AKNS方法構(gòu)造Lax對(duì)并推導(dǎo)出達(dá)布變換公式,再應(yīng)用此公式求得方程在零種子解情形下的孤子解.下一步,將從方程的非零種子解出發(fā)進(jìn)行求解,根據(jù)線性偏微分方程的疊加原理,將特征函數(shù)線性疊加組成新的特征函數(shù),從而得到方程的呼吸子解,并對(duì)此呼吸子解進(jìn)行泰勒級(jí)數(shù)展開,最后得到怪波解.

[1] 陳美同. Painleve方法構(gòu)造非線性偏微分方程的精確解[D]. 錦州:渤海大學(xué),2015.

[2] 羅琳,徐國進(jìn). 對(duì)tanh函數(shù)法的推廣及其在非線性方程中的應(yīng)用[J]. 湖北工程學(xué)院學(xué)報(bào),2004,24(3):58-62.

[3] WANG M, ZHOU Y, LI Z. Application of a homogeneous balance method to exact solutions of nonlinear equations in Mathematical Physics[J]. Physics Letters,1996,A216(1/5):67-75.

[4] ANKIEWICZ A, SOTO-CRESPO J M, AKHMEDIEV N. Rogue waves and rational solutions of the Hirota equation[J]. Physical Review Statistical Nonlinear and Soft Matter Physics,2010,E81(4):387-395.

[5] WU J P. A new wronskian condition for a (3+1)-dimensional nonlinear evolution equation[J]. Chinese Physics Letters,2011,28(5):50501-50503(3).

[5] 姜丙利,柳銀萍. 不變子空間方法及一個(gè)非線性演化方程的精確解[J]. 浙江師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,36(2):155-160.

[7] 谷超豪,胡和生,周子翔. 孤立子理論中的達(dá)布變換及其幾何應(yīng)用[M]. 2版. 上海:上海科學(xué)技術(shù)出版社,2005.

[8] GUO B, LING L, LIU Q P. Nonlinear schrodinger equation:generalized darboux transformation and rogue wave solutions[J]. Physical Review Statistical Nonlinear and Soft Matter Physics,2012,E85(2):317-344.

[9] 張金順,李華夏. 2+1維Levi孤子方程的Darboux變換[J]. 鄭州大學(xué)學(xué)報(bào)(理學(xué)版),2001,23(3):13-17.

[10] ZAKHAROV V E, SHABAT A B. A scheme for integrating the nonlinear equations of mathematical physics by the method of the inverse scattering problem. I[J]. Functional Analysis and Its Applications,1974,8(3):226-235.

[11] ABLOWITZ M J, KAUP D J, NEWELL A C, et al. Nonlinear-evolution equations of physical significance[J]. Physical Review Letters,1973,31(2):125-127.

[12] MATVEEV V B, SALLE M A. Darboux transformations and solitons[J]. J Neurochemistry,1991,42(6):1667-1676.

[13] 李翊神. 規(guī)范變換,貝克隆變換與非線性疊加公式[J]. 數(shù)學(xué)進(jìn)展,1989,5(3):356-372.

[14] 広田良吾,紅艷. 孤子理論中的直接方法[M]. 北京:清華大學(xué)出版社,2008.

[15] PERKINS S. Dashing rogues:freak ocean waves pose threat to ships, deep-sea oil platforms[J]. Science News,2006,170(21):328-329.

[16] KIBLER B, FATOME J, FINOT C, et al. The Peregrine soliton in nonlinear fibre optics[J]. Nature Physics,2010,6(10):790-795.

[17] HAVER S. A possible freak wave event measured at the Draupner Jacket january 1 1995[J]. Actes De Colloques-IFREMER,2004,56(2):24-35.

[18] DENSCHLAG J, SIMSARIAN J E, FEDER D L, et al. Generating solitons by phase engineering of a Bose-Einstein condensate[J]. Science,2000,287(5450):97-101.

[19] SEAMAN B T, CARR L D, HOLLAND M J. Effect of a potential step or impurity on the Bose-Einstein condensate mean field[J]. Physical Review,2004,A71(3):309-315.

[20] LOSSEVA T V, POPEL S I, GOLUB A P, et al. Weakly dissipative dust-ion-acoustic solitons in complex plasmas and the effect of electromagnetic radiation[J]. Physics of Plasmas,2012,19(1):133-144.

[21] 郭柏靈. 孤立子[M]. 北京:科學(xué)出版社,1987.

[22] 李德孝. 一個(gè)非線性方程的延拓結(jié)構(gòu)[J]. 青海師范大學(xué)學(xué)報(bào)(自科版),2012,28(3):16-18.

[23] 秦振云. 高階矩陣譜問題非線性化與孤立子方程的擬周期解[D]. 上海:復(fù)旦大學(xué),2006.

[24] 田守富. 非線性微分方程的若干解析解方法與可積系統(tǒng)[D]. 大連:大連理工大學(xué),2012.

[25] 陶勇勝. Hirota型方程的孤立子解和怪波解[D]. 寧波:寧波大學(xué),2012.

[26] 郭玉翠. 非線性偏微分方程引論[M]. 北京:清華大學(xué)出版社,2008.

[27] 孫業(yè)朋. 可積系統(tǒng)與非等譜孤子方程的求解[D]. 上海:上海大學(xué),2006.