非線性Schr?dinger方程的對稱約化和精確解

曹 瑞

(菏澤學院 數學系，山東 菏澤 274015)

非線性Schr?dinger方程的對稱約化和精確解

曹 瑞*

(菏澤學院 數學系，山東 菏澤 274015)

本文研究一類立方非線性Schr?dinger方程的對稱約化和精確解問題。 首先，利用直接對稱方法，得到非線性Schr?dinger方程的對稱；其次，根據求解相應的特征方程獲得非線性Schr?dinger 方程的相似約化；最后，結合輔助方程獲得非線性Schr?dinger方程的精確解。 這些解包括孤立波解、Jacobi橢圓函數解以及三角函數解。

非線性Schr?dinger方程的對稱；對稱約化；精確解；孤立子

作為描述復雜非線性現象的數學模型，非線性發展方程涉及了眾多自然科學領域，如物理學，化學，生物，工程等。 非線性發展方程的精確解在解釋復雜非線性現象中有著重要作用。 為了尋找非線性發展方程的精確解，許多專家和學者提出了一系列行之有效的求解非線性發展方程精確解的方法，例如齊次平衡法[1]，Painleve截尾展開法[2]，Hirota直接法[3]，sine- cosine方法[4]，雙曲函數法[5]，試探函數法[6]，Jacobi橢圓函數展開法[7]以及Jacobi橢圓函數展開方法一般化的F-展開法[8]，G'/G展開法[9，10]，對稱約化法[11，12]等。

利用這些方法得到了非線性發展方程許多豐富的精確解，其中包括孤立波解，周期波解，激波解等。

文獻[13]中研究了如下非線性Schr?dinger方程

其中p，q是非零實常數，ψ=ψ(x，t) 是關于變量x，t的復值函數。 N.Taghizadeh 等運用首次積分方法構造了非線性Schr?dinger方程的精確解。 當p=1，q=μ時，我們可以得到文獻[14]中的方程，Ma 和Chen 利用直接方法獲得了非線性Schr?dinger方程的精確解。 在本文中我們討論如下的立方非線性Schr?dinger方程

(1)

其中ψ=ψ(x，t)是復函數，i2=-1。 非線性Schr?dinger方程(1)與它的各種推廣形式在非線性光學、等離子物理以及玻色-愛因斯坦凝聚等領域有著廣泛應用[15-17]。 本文的目的是通過直接對稱方法尋找非線性Schr?dinger方程(1)的精確解。

本文第一部分，介紹直接對稱方法的基本思想；第二部分，利用直接對稱方法構造非線性Schr?dinger方程的對稱；然后，根據求解相應的特征方程組得到非線性Schr?dinger方程的相似約化；最后，通過一個四階輔助方程獲得了非線性Schr?dinger方程的精確解；第三部分，得出一些結論。

1 直接對稱方法

下面列出直接對稱方法的基本步驟: 對于給定的非線性偏微分方程

F(x，t，u，ut，ux，uxx，…)=0，

(2)

其中F是關于u及其導數的函數。 利用直接對稱方法求方程對稱的基本步驟為：

第一步：假設方程(2)有下列形式的對稱

σ=α(x，t)ut+β(x，t)ux+γ(x，t)u+δ(x，t)，

(3)

其中α，β，γ，δ都是x，t的待定函數，而函數u滿足方程(2)。

(4)

將對稱(3)代入(4)式，得到關于u的決定方程組，令u的各階導數系數為零，解的α，β，γ，δ，進而可求得相應的對稱σ。

第二步：通過對稱σ為零求解一階偏微分方程，可以得到群不變量。

第三步：利用第二步得到的不變量約化原方程。

2 立方非線性Schr?dinger方程的對稱約化和精確解

2.1 方程的對稱

為了求出方程(1)的對稱，我們假定

σ=a(x，t)ψt+b(x，t)ψx+c(x，t)ψ+d(x，t)，

(5)

這里a(x，t)，b(x，t)，c(x，t)，d(x，t) 是待定函數，其中σ滿足下面的方程

(6)

將方程(5)代入方程(6)，并考慮方程(1)以及如下方程

(7)

a(x，t)=2c4t+c2，b(x，t)=c4x-c3t+c1，

c(x，t)=ic3x+ic5+c4，d(x，t)=0。

(8)

其中c1，c2，c3，c4，c5是任意常數。 由此可得所求對稱

σ= (2c4t+c2)ψt+(c4x-c3t+c1)ψx+

(ic3x+ic5+c4)ψ，

(9)

其交換圖如圖1所示。

[σi，σj]σ1σ2σ3σ4σ5σ100σ5σ10σ200-σ12σ20σ3-σ5σ10-σ30σ4-σ1-2σ2σ300σ500000

圖1 交換圖

(10)

從向量場(10)，可得對應的五維李代數L={V1，V2，V3，V4，V5} ，有下列一組基

(11)

并且，可以得到與向量場(11)相對應的單參數變換群

g1:(x，t，ψ)→(x+ε，t，ψ)，

g2:(x，t，ψ)→(x，t+ε，ψ)，

g4:(x，t，ψ)→(eεx，e2εt，e-εψ)，

g5:(x，t，ψ)→(x，t，e-iεψ)。

(12)

這里ε群參數，g1是空間變換 ，g2是時間變換，g3是伽利略變換 ，g4是尺度變換，g5是標量變換。

根據上述單參數群可知，如果f(x，t)是方程(1)的解，那么ψi(i=1，2，3，4，5)也是方程(1)的解。

ψ1(x，t)=f(x-ε，t)，

ψ2(x，t)=f(x，t-ε)，

ψ4(x，t)=e-εf(e-εx，e-2εt)，

ψ5(x，t)=e-iεf(x，t)。

(13)

2.2 方程的對稱約化和精確解

本節利用方程σ=0的相容性求方程(1)的對稱約化及精確解。 方程(1)的對稱σ=(2c4t+c2)ψt+(c4x-c3t+c1)ψx+(ic3x+ic5+c4)ψ。

特征方程為

1. 當c1=0，c2=0，c3=0，c4≠0，c5=0時，求解特征方程

可得

此時，f(ξ)滿足方程

(14)

2.當 c1≠0，c2≠0，c3=0，c4=0，c5≠0時，求解特征方程

可得

ξ=x-vt，ψ=f(ξ)e-iωt。

此時，f(ξ)滿足方程

(15)

并且若f(ξ)是上式的解，那么ψ=f(ξ)e-iωt就是方程(1)的解。

下面我們求方程(1)的精確解。 假設

ψ(x，t)=v(x，t)exp(iη)，η=αx+βt，

(16)

其中v(x，t)是實值函數，且α，β是實數。 把(16)代入，方程(1)成為一個實系統

vt+αvx=0，vxx-(α2+2β)v+2μv3=0。

(17)

令v(x，t)=v(ξ)，ξ=k(x-αt)，k≠0。方程(17)可化為

k2v″-(α2+2β)v+2μv3=0，

(18)

上式乘以2v′積分一次得到

(19)

情形1 Jacobi橢圓函數解

(1.1)當a0=1，a2=-(1+m2)，a4=m2時，

ψ1=eiηsn(ξ)，ψ2=eiηcd(ξ)。

(1.2)當a0=1-m2，a2=2m2-1，a4=-m2時，v(ξ)=cn(ξ)， 此時

ψ3=eiηcn(ξ)。

(1.3)當a0=m2-1，a2=2-m2，a4=-1時，v(ξ)=dn(ξ)，此時

ψ4=eiηdn(ξ)。

(1.4)當a0=m2，a2=-(1+m2)，a4=1時，

ψ5=eiηns(ξ)，ψ6=eiηdc(ξ)。

(1.5)當a0=-m2，a2=2m2-1，a4=1-m2時，

ψ7=eiηnc(ξ)。

(1.6)當a0=-1，a2=2-m2，a4=m2-1時，v(ξ)=nd(ξ)， 此時

ψ8=eiηnd(ξ)。

(1.7)當a0=1，a2=2-m2，a4=1-m2時， v(ξ)=sc(ξ)，此時

ψ9=eiηsc(ξ)。

(1.8)當a0=1，a2=2m2-1，a4=-m2(1-m2)時，v(ξ)=sd(ξ)，此時

ψ10=eiηsd(ξ)。

(1.9)當a0=1-m2，a2=2-m2，a4=1時，v(ξ)=cs(ξ)，此時

ψ11=eiηcs(ξ)。

(1.10)當a0=-m2(1-m2)，a2=2m2-1，a4=1時，v(ξ)=ds(ξ)，此時

ψ12=eiηds(ξ)。

ψ13=eiη(ns(ξ)±cs(ξ))。

ψ14=eiη(nc(ξ)±sc(ξ))。

ψ15=eiη(ns(ξ)±ds(ξ))。

ψ16=eiη(sn(ξ)±icn(ξ))。

情形2 三角函數解

(2.3)當m→0時，情形1中的橢圓函數解ψ9退化三角函數解

ψ19=eiηtan(ξ)。

情形3 孤立子解

(3.3)當m→1時，情形1中的橢圓函數解ψ3退化亮孤子解

ψ22=eiηsech(ξ)。

(3.4)當m→1時，情形1中的橢圓函數解ψ1退化暗孤子解

ψ23=eiηtanh(ξ)。

注：(1)當模m→1時，Jacobi橢圓函數退化為雙曲函數。

snξ→tanhξ，cnξ→sechξ，scξ→sinhξ，

csξ→cschξ，ncξ→coshξ。

(2)當模m→0時，Jacobi橢圓函數退化為三角函數。

snξ→sinξ，cnξ→cosξ，scξ→tanξ，csξ→cotξ，

ncξ→secξ，nsξ→cscξ。

3 結語

本文把直接對稱方法應用于具有重要物理背景的立方非線性Schr?dinger方程，成功得到了立方非線性Schr?dinger方程的對稱。 并根據所得對稱相應的特征方程求出了方程的不變量和相似約化。 再通過行波約化，結合四階輔助方程獲得了立方非線性Schr?dinger方程的豐富的精確解，這些解包括孤立波解、Jacobi橢圓函數解以及三角函數解。當Jacobi橢圓函數解中參數取特殊值時可以得到雙曲函數解、孤立波解。 這些精確解豐富了立方非線性Schr?dinger方程精確解的解系，同時，由于該方程是光學中的重要方程，這些得到的精確解將有助于物理上對該方程的研究。

[1] WANG M L，ZHOU Y B，LI Z B. Application of a homogeneous balance method to exact solutions of nonlinear equations in mathematical physics[J]. Phys Lett A，1996，216:67-75.

[2] 張解放，陳芳躍.截斷展開方法和廣義變系數Kdv方程[J].物理學報，2001，50(9):1648-1650.

[3] Hirota R. The Direct Methods in Soliton Theory[M].Cambridge:Cambridge University，2004.

[4] Yan C T. A simple transformation for nonlinear waves[J].Phys Lett A，1996，224:77-84.

[5] 郭冠平，張解放. 關于雙曲函數方法求孤波解的注記[J]. 物理學報，2002，51(6):1159-1162.

[6] Cao Rui，Zhang Jian.Trial function method and exact solutions to the generalized nonlinear Schr?dinger equation with time- dependent coefficient[J].Chinese Physics B，2013，22(10):100507-1-100507-4.

[7] 向以華，石義霞. (2+1)維色散長波方程的擴展橢圓函數有理展開解法[J]. 數學雜志，2009(2):206-210.

[8] 曹瑞. 改進的F-展開方法和藕合Klein- Gordon方程的精確解[J]. 蘭州大學學報(自然科學版)，2007，43(6):112-116.

[9] BIAN Chunquan，PANG Jing，JIN Linghua，et al. Solving two fifth order strong nonlinear evolution equations by using the (G/G')- expansion method[J]. Communications in Nonlinear Science and Numerical Simulation，2010，15(9):2337-2343.

[10] 李幫慶，馬玉蘭. G/G'展開法和(2+1)維非對稱Nizhnik- Novikov- Veselov系統的新精確解[J].物理學報，2009，58(7):4373-4378.

[11] MIAO Q，XIN X P，CHEN Y. Nonlocal symmetries and explicit solutions of the AKNS system[J]. Applied Mathematics Letters，2014，28:713.

[12] CHEN J C，XIN X P，CHEN Y. Symmetry Analysis and Conservation Laws to the (2+1)- Dimensional Coupled Nonlinear Extension of the Reaction- Diffusion Equation[J]. Commun Theor Phys， 2014，62:173-182.

[13] Taghizadeh N，Mirzazadeh M，Farahrooz F. Exact solutions of the nonlinear Schr?dinger equation by the first integral method[J]. J Math Anal Appl， 2011，374:549-553.

[14] MA W X，CHEN M. Direct search for exact solutions to the nonlinear Schr?dinger equation[J]. Appl Math Comput， 2009，215:2835-2842.

[15] Zabusky N J，Kruskal M D. Interaction of "soliton" in a collisionless plasma and the recurrence of initial states[J]. Phys Rev Lett，1965，19:1095-1097.

[16] Ablowitz M J，Segur H. Solitons and the Inverse Scattering Transformation[M]. Philadelphia:SIAM，1981.

[17] Hasegawa A，Kodama Y. Solitons in Optical Communications[M].New York:Oxford University Press，1995.

[18] 劉式適，劉式達.物理學中的非線性方程[M].北京: 北京大學出版社，2000.

(責任編輯：周曉南)

Symmetry Reduction and Exact Solutions of Nonlinear Schr?dinger Equation

CAO Rui*

(Department of Mathematics, Heze University, Heze 274000, China)

Symmetry of nonlinear Schr?dinger equation was investigated. Applying direct symmetry method, symmetry of nonlinear Schr?dinger equation was derived. By solving corresponding characteristic equation associated with the symmetry equation, symmetry reduction of nonlinear Schr?dinger were obtained based on a fourth- order auxiliary equation, these solutions include solitary wave solutions, Jacobi elliptic function solutions and trigonometric function solutions.

nonlinear Schr?dinger equation; direct symmetry method; exact solution; soliton

1000-5269(2016)06-0001-04

10.15958/j.cnki.gdxbzrb.2016.06.01

2016-05-30

國家自然科學基金項目資助(11347102，11401409)；山東省自然科學基金項目資助(ZR2011AL018，ZR2011AQ008)；山東省高?？萍加媱濏椖抠Y助(J13LI02)；荷澤學院自然科學基金項目資助(xy14KJ04)

曹瑞(1979-)，女，副教授，研究方向: 孤立子和可積系統，Email:ruicao999@126.com.

曹瑞，Email:ruicao999@126.com.

O175.2

A