張道華，李華新，朱國良，石大慶，李金星

(1.湖北師范學院 物理與電子科學學院，湖北 黃石 435002；2.夏店鎮(zhèn)中學，湖北 大悟 432804)

0 引言

自20世紀60年代以來，非線性科學飛速發(fā)展，與此對應(yīng)，物理學中的非線性知識也日趨豐富，而作為自然界中普遍存在的孤子，是最早能夠觀察到并且可以在實驗室產(chǎn)生的非線性現(xiàn)象之一。隨著孤子問題研究的深入，對它的理論建立、研究方法和實際應(yīng)用的研究，已經(jīng)在物理學各研究領(lǐng)域中占據(jù)了越來越重要的地位。在討論固體的非線性物理問題時，一維非線性晶格常常因為模型簡單、物理圖像清晰，并且存在各種非線性元激發(fā)而被普遍重視和研究[1～3]。其中一維FPU晶格模型是現(xiàn)有描述固體晶格非線性模型中最經(jīng)典、最簡單的模型之一，對FPU晶格動力學的研究導(dǎo)致了內(nèi)稟局域模的誕生，促進了晶格中孤子理論的發(fā)展[4～6]。在FPU晶格模型中，原子鏈上的非線性元激發(fā)包含了包絡(luò)孤子、扭結(jié)包絡(luò)孤子、反扭包絡(luò)結(jié)孤子等，已被廣泛地應(yīng)用于一維非線性晶格局域模、熱傳導(dǎo)現(xiàn)象和混沌現(xiàn)象等的研究之中。對于眾多的非線性晶格模型來說，除了有限幾個的運動學方程可積外，大多數(shù)其它模型的運動學方程都是不可積的，F(xiàn)PU晶格模型就是其中的一種。為了求解這些不可積的運動學方程，人們發(fā)展了很多的方法。起初，主要是采用連續(xù)性近似。但是在特征波長與晶格常數(shù)相差不大時，忽略晶格的離散性就會帶來較大的誤差[2]。為了解決這個問題，許多求解非線性離散晶格體系的運動學方程的方法應(yīng)運而生，如格林函數(shù)方法、準連續(xù)性近似法、多重尺度法、旋轉(zhuǎn)波近似等。對于FPU晶格模型，在長波近似下采用多重尺度法可以很方便的找出其孤子解，在處理這一問題上，多重尺度法較其他方法有一定的優(yōu)勢。

本文主要是利用多重尺度法結(jié)合長波近似求解FPU晶格模型中的運動學方程，研究非線性晶格中的孤子。

1 模型和方法

對于由n個質(zhì)量為m的原子、分子或原子集團組成的FPU晶格體系，哈密頓量通?？杀硎緸閇2，6]：

(1)

其中un表示第n個原子相對平衡位置的位移，右邊第一項表示動能，第二項是相鄰粒子之間的相互作用勢，可表示為：

k2是線性系數(shù)，k3、k4為非線性系數(shù)。對于特殊情況，當k3=0，k4≠0時，體系對應(yīng)FPU-α模型；當k3≠0，k4=0時，體系對應(yīng)FPU-β模型[6]。近幾年對于FPU模型的研究主要集中于FPU-β模型，其數(shù)學處理過程往往比FPU-α模型簡單，甚至有些方法只對前者適用。本文也僅以前者為例。

對FPU-β模型應(yīng)用海森堡運動方程可得[2]：

(2)

做變換un→εφn，ε是小參量，(2)式可變?yōu)椋?/p>

(3)

假定孤子解的形式如下：

φn=φn(t1,t2,εna)ei(kna-ωt)+c.c

(4)

上式中a和k分別是晶格常量和波數(shù)，φn是孤子的慢變外在包絡(luò)函數(shù)，ei(kna-ωt)是快變載波。引入多尺度變量t1=εt，t2=ε2t,x1=εx， 由(4)我們可以得到：

(5)

(6)

(7)

應(yīng)用長波近似，并將(4)待入(3)式，這里我們僅保留到O(ε3) 項，其線性項可化為：

(8)

非線性項：

(9)

(10)

將(6)～(10)代入(3)式，并比較等式兩端參量ε系數(shù)得：

ε0:mω3φn=2k2(1-coska)φn

(11)

(12)

(13)

(14)

顯然有(11)可得色算關(guān)系：

(15)

因此我們很容易得出孤子的群速度為：

(16)

有(12)～(13)可得：

(17)

上式即為標準的非線性薛定諤方程。

2 結(jié)果推導(dǎo)

對比標準非線性薛定諤方程的解，我們可以得到(17)式晶格模型中的孤子解。方程(17)式可直接寫成以下非線性薛定諤方程[3]：

(18)

其中

(19)

設(shè)(18)式有包絡(luò)波解：

φ=u(ξ)ei(kx-ωt),ξ=x-vgt

(20)

將(20)代入(18)式，得到

(21)

因u(ξ) 為實函數(shù)，故要求du/dξ前的復(fù)系數(shù)為零，方程簡化為：

(22)

當ω=(k2+2)α時，方程(22)有暗孤子解[7]：

(23)

當ω=(k2-4)α時，方程有亮孤子解：

(24)

當ω=(k2-1)α時，方程有亮孤子解：

(25)

上面的孤子解表明，在適當?shù)臈l件下，一維FPU原子鏈中能夠形成亮孤子和暗孤子。

3 總結(jié)

本文在長波近似的條件下，運用多重尺度法對一維非線性振動體系進行了研究，對于FPU-β模型，很方便的將其運動方程變換成了標準的非線性薛定諤方程，根據(jù)非線性薛定諤方程的一般求解方法，得到了一維FPU原子鏈中的亮孤子和暗孤子解。這將有利于促進晶格中孤子理論的發(fā)展。

參考文獻：

[1]周凌云,王瑞麗,吳光敏,等.非線性物理理論及應(yīng)用[M]. 北京:科學出版社，2000.

[2]龐小峰. 孤子物理學[M]. 成都: 四川科學技術(shù)出版社, 2003.

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

[4]Chen W Z.Experimental observation of solitons in a 1D nonlinear lattice [J].Phys Rev B, 1994, 49(21): 15063～15066.

[5]Ford J.The Fermi-Pasta-Ulma Problem: Paradox Turns Discovery[J].Phys Rep,1992, (213): 271～310.

[6]Berman G P,Izraileva F M.The Fermi-Pasta-Ulam problem: Fifty years of progress[J].chaos, 2005, (15):1～18.

[7]Liu Ji-Bing, Xie Xiao-Tao,Wu Ying.Acoustic solitary wave and breather in molecular magnets via electromagnetically induced Transparency[J].EPL,2010 (89 ):1～6.