周艷珍，張素英

（山西大學 理論物理研究所，山西 太原 030006）

＊兩組分玻色－愛因斯坦凝聚體中矢量孤子的動力學性質

周艷珍，張素英＊

（山西大學 理論物理研究所，山西 太原 030006）

通過數(shù)值求解異種兩組分玻色愛因斯坦凝聚體在弱囚禁勢中的運動方程來討論其矢量孤子解的動力學性質.研究表明，種內和種間相互作用強度滿足不同的條件時，會形成亮亮孤子、亮暗孤子和暗暗孤子等不同的矢量孤子解.其中亮亮孤子和亮暗孤子是穩(wěn)定的，而暗暗孤子很不穩(wěn)定.適當改變種間相互作用強度，亮、暗孤子之間能夠相互轉換.

兩組分玻色愛因斯坦凝聚體；矢量孤子；穩(wěn)定性

玻色－愛因斯坦凝聚（簡稱BEC）中的暗、亮物質波孤子的成功觀測及其潛在的應用前景，使BEC中的物質波孤子成為當前低溫物理和凝聚態(tài)物理研究領域的研究熱點之一.孤子是非線性物質波的存在形式，能夠很好地保持空間輪廓，穩(wěn)定地傳播，它分為亮孤子和暗孤子兩類.目前，多組分BEC由于其具有復雜的量子相位［1］和在單組分BEC中所沒有的性質已成為研究暗、亮孤子的焦點之一.而且，理論［2－5］和實驗［6－7］的研究都表明異種兩組分BEC系統(tǒng)中種間相互作用強度起著很重要的作用［8－9］.本文通過數(shù)值模擬的方法研究了兩組分玻色－愛因斯坦凝聚體（g11＜0且g22＞0）中矢量孤子在弱囚禁勢中的動力學性質，并通過調節(jié)種間相互作用強度（實驗上采用Feshbach共振技術實現(xiàn)［10－11］）研究了孤子之間的相互轉換問題.

1 理論模型

兩組分玻色－愛因斯坦凝聚體的平均場動力學控制方程如下［2－4］：

2 不同的矢量孤子解

當無束縛外勢時，即λ1＝0，λ2＝0，方程（1）和（2）存在以下四種情況的矢量孤子解［12］：

為了得到不同的矢量孤子解，我們定義下面兩個量：

圖1 種間相互作用強度g12取不同值時凝聚體內形成不同的矢量孤子所對應的分段圖，其中κg 11＜g 22／κ，g 11＜0，g 22＜0Fig.1 Chart in stages of vector solitons in Bose－Einstein Condensates with different values of the interspecies interaction strength g 12，withκg 11＜g 22／κ，g 11＜0 and g 22＜0

第二類情況——g11＜0，g22＞0，g12可以調節(jié)，我們可以劃分出下面三個區(qū)域（如圖2，P666）：

f.g12∈（－∞，κg11］，此時Q1＞0，Q2＜0，形成亮 －亮孤子；

g.g12∈（κg11，g22／κ］，此時Q1＞0，Q2＞0，形成亮 －暗孤子；

h.g12∈（g22／κ，＋∞），此時Q1＜0，Q2＞0，形成暗 －暗孤子；

圖2 種間相互作用強度g12取不同值時凝聚體內形成不同的矢量孤子所對應的分段圖，其中κg 11＜g 22／κ，g 11＜0，g 22＞0Fig.2 Chart in stages of vector solitons in Bose－Einstein Condensates with different values of the interspecies interaction strength g 12，withκg 11＜g 22／κ，g 11＜0 and g 22＞0

圖3 種間相互作用強度g 12取不同值時凝聚體內形成不同的矢量孤子所對應的分段圖，其中κg 11＜g 22／κ，g 11＞0，g 22＞0Fig.3 Chart in stages of vector solitons in Bose－Einstein Condensates with different values of the interspecies interaction strength g 12，withκg 11＜g 22／κ，g 11＞0 and g 22＞0

關于亮－亮孤子的性質和暗－亮孤子的性質在異種兩組分玻色－愛因斯坦凝聚體（g11＞0且g22＞0）中已經(jīng)有大量的研究［13－14］.我們采用數(shù)值模擬的方法主要討論了在異種兩組分玻色－愛因斯坦凝聚體（g11＜0且g22＞0）中亮－亮孤子、亮－暗孤子和暗－暗孤子的穩(wěn)定性及其相互轉換.

3 不同孤子的穩(wěn)定性

兩種類玻色－愛因斯坦凝聚體中，我們選定第一種類為7Li原子，第二種類為23Na原子，并且種內相互作用強度滿足g11＜0，g22＞0，種間相互作用強度g12可以調節(jié).此條件滿足上面所列三類情況中的第二類情況，g12滿足不同的條件時可以形成亮－亮、亮－暗和暗－暗三種不同孤子對.下面我們通過數(shù)值模擬實驗來研究這三種孤子的穩(wěn)定性.初始態(tài)（t＝0）分別選用精確解（3）－（8），代入方程（1）和（2）進行動力學演化.結果表明，亮－亮孤子和亮－暗孤子在弱囚禁勢中都能夠穩(wěn)定地傳播，如圖4和圖5（P667）所示.而暗－暗孤子很不穩(wěn)定，迅速地被破壞，且由于種間強大的排斥作用，彼此遠離，如圖6（P667）所示.

另外，可調相互作用的兩組分玻色－愛因斯坦凝聚體中不同的孤子之間可以相互轉換，如圖7所示.首先，我們選擇一對速度為零的亮 －暗孤子且種間相互作用強度滿足g12∈（κg11，g22／κ］，代入方程（1）和（2）進行動力學演化，然后在t＝5時將種間相互作用強度g12瞬間變?yōu)樾∮讦蔳11的值，繼續(xù)動力學演化.由于g12＜κg11，滿足Q1＞0，Q2＜0，所以凝聚體內亮－暗孤子在t＝5時轉換為亮－亮孤子，且不同種類凝聚體中形成的亮孤子經(jīng)過幾次碰撞后合并為一個新的亮孤子.

4 結論

圖4 在弱囚禁勢中亮－亮孤子的動力學演化圖其中g 11＝－3.47，g 22＝5.89，g 12＝－3，κ＝7／23，λ1＝λ2＝0.12，η＝1.5，v＝0.8Fig.4 Dynamic evolution of bright－bright solitons in weak longitudinal trap with g 11＝－3.47，g 22＝5.89，g 12＝－3，κ＝7／23，λ1＝λ2＝0.12，η＝1.5 and v＝0.8

圖5 在弱囚禁勢中亮－暗孤子的動力學演化圖其中g11＝－3.47，g22＝5.89，g12＝3，κ＝7／23，λ1＝λ2＝0.12，η＝1.5，v＝0.4Fig.5 Dynamic evolution of bright－dark solitons in weak longitudinal trap with g 11＝－3.47，g 22＝5.89，g 12＝3，κ＝7／23，λ1＝λ2＝0.12，η＝1.5 and v＝0.4

圖6 在弱囚禁勢中暗－暗孤子的動力學演化圖其中g 11＝－3.47，g 22＝5.89，g 12＝30，κ＝7／23，λ1＝λ2＝0.12，η＝1.5，v＝0Fig.6 Dynamic evolution of dark－dark solitons in weak longitudinal trap with g 11＝－3.47，g 22＝5.89，g 12＝30，κ＝7／23，λ1＝λ2＝0.12，η＝1.5 and v＝0

圖7 在弱囚禁勢中亮－暗孤子相互轉換的動力學演化圖.初始態(tài)與圖5相同取g 11＝－3.47，g 22＝5.89，g 12＝3，在t＝5后，種間相互作用強度變?yōu)間12＝－20，其它參數(shù)不變，其中κ＝7／23，λ1＝λ2＝0.12，η＝1.5，v＝0Fig.7 Dynamic evolution of the mutual transforming bright and dark solitons in weak longitudinal trap withg 11＝－3.47，g22＝5.89，g12＝3 initially，and the interspecies interaction strength is changed tog 12＝－20 at t＝5，other parameters are changeless，andκ＝7／23，λ1＝λ2＝0.12，η＝1.5，v＝0

本文對可調相互作用的異種兩組分玻色－愛因斯坦凝聚體中孤子的研究表明，種內和種間相互作用強度滿足一定的條件時，會形成不同的矢量孤子解.通過數(shù)值模擬方法，調節(jié)相互作用強度，對所形成的亮－亮孤子、亮－暗孤子和暗－暗孤子在弱囚禁勢中進行動力學演化，結果顯示亮－亮孤子和亮－暗孤子是穩(wěn)定的，而暗－暗孤子很不穩(wěn)定.此外，當改變組分間相互作用強度時，亮、暗孤子之間可以相互轉換.這對進一步研究矢量孤子在多種類BEC中的性質有重要的意義.

［1］ Alon O E，Streltsov A I，Cederbaum L S.Zoo of Quantum Phases and Excitations of Cold Bosonic Atoms in Optical Lattices［J］.PhysRevLett，2005，95：030405.

［2］ Ho T L，Shenoy V B.Binary Mixtures of Bose Condensates of Alkali Atoms［J］.PhysRevLett，1996，77：3276－3279.

［3］ Esry B D，Greene C H，Burke J P，etal.Hartree－Fock Theory for Double Condensates［J］.PhysRevLett，1997，78：3594－3597.

［4］ Pu H，Bigelow N P.Properties of Two－Species Bose Condensates［J］.PhysRevLett，1998，80：1130－1133.

［5］ Pu H，Bigelow N P.Collective Excitations，Metastability，and Nonlinear Response of a Trapped Two－Species Bose－Einstein Condensate［J］.PhysRevLett，1998，80：1134－1137.

［6］ Myatt C J，Burt E A，Ghrist R W，etal.Production of Two Overlapping Bose－Einstein Condensates by Sympathetic Cooling［J］.PhysRevLett，1997，78：586－589.

［7］ Modugno G，F(xiàn)errari G，Roati G，etal.Bose－Einstein Condensation of Pota－ssium Atoms by Sympathetic Cooling［J］.Science，2001，294：1320.

［8］ Thalhammer G，Barontini G，Sarlo L De，etal.Double Species Bose－Einstein Condensate with Tunable Interspecies Interactions［J］.PhysRevLett，2008，100：210402.

［9］ Papp S B，Pino J M，Wieman C E.Tunable Miscibility in a Dual－Species Bose－Einstein Condensate［J］.PhysRevLett，2008，101：040402.

［10］ Roberts J L，Claussen N R，Burke J P，etal.Resonant Magnetic Field Control of Elastic Scattering in Cold85Rb ［J］.PhysRevLett，1998，81：5109－5l12.

［11］ Stenger J，Inouye S，Andrews M R，etal.Strongly Enhanced Inelastic Collisions in a Bose－Einstein Condensate near Feshbach Resonances［J］.PhysRevLett，1999，82：2422－2425.

［12］ Liu X X，Pu H，Xiong B，etal.Formation and Transformation of Vector Solitons in Two－Species Bose－Einstein Condensates with a Tunable Interaction［J］.PhysRevA，2009，79：013423.

［13］ Pérez－García Víctor M，Beitia Juan Belmonte.Symbiotic Solitons in Heteronuclear Multi－component Bose－Einstein Condensates［J］.PhysRevA，2005，72：033620.

［14］ Busch Th，Anglin J R.Dark－Bright Solitons in Inhomogeneous Bose－Einstein Condensates［J］.PhysRevLett，2001，87：010401.

Dynamic Properties of Vector Solitons in Two－species Bose－Einstein Condensates

ZHOU Yan－zhen，ZHANG Su－ying
（InstituteofTheoreticalPhysics，ShanxiUniversity，Taiyuan030006，China）

We discuss the dynamic properties of the vector solitons of two－species Bose－Einstein condensates in weak trapping potential by numerically solving its corresponding equations of motion.The results show that under some condition，different intraspecies and interspecies interaction coefficients result in different types of vector solitons，such as bright－bright，bright－dark and dark－dark vector solitons.Moreover，the bright－bright and bright－dark vector soliton are dynamically stable，while the dark－dark vector soliton is dynamically unstable.In addition，the different types of vector solitons can be transformed by tuning the interspecies interaction via Feshbach resonance.

two－species Bose－Einstein condensates；vector soliton；stability

O469

A

2011－09－27；

2012－02－24

國家自然科學基金（10972125）；高等學校博士學科點專項科研基金（20111401110004）；山西省自然科學基金（2010011001－2）；山西省留學回國人員科研基金

周艷珍（1988－），女，山西忻州人，研究方向：凝聚態(tài)物理.E－mail：Zhouyanzhen1988＠163.com.＊張素英（1967－），女，山西人，教授，博士生導師，研究方向：數(shù)學物理、計標物理.E－mail：zhangsu＠sxu.edu.cn

0253－2395（2012）04－0664－05

book=668,ebook=358