劉新建，李衛(wèi)東

(山西大學(xué)理論物理研究所，中國 太原 030006)

?

弱聯(lián)接玻色愛因斯坦凝聚體中勢壘寬度對非線性耦合及其動力學(xué)的影響

劉新建，李衛(wèi)東

(山西大學(xué)理論物理研究所，中國 太原 030006)

利用解析與數(shù)值方法，對處于對稱雙勢阱中的玻色愛因斯坦凝聚體中，勢壘寬度對系統(tǒng)非線性耦合及其動力學(xué)的影響進行了研究．研究發(fā)現(xiàn)當(dāng)勢壘寬度較大時，系統(tǒng)的線性耦合強度可迅速減??；在勢壘寬度大于0.3且非線性強度較大時，線性耦合強度遠小于非線性耦合項，此時玻色約瑟夫森結(jié)模型的動力學(xué)特性由非線性耦合強度來決定．同時對勢壘寬度對BEC約瑟夫森振蕩的周期和發(fā)生宏觀量子自俘獲時的非線性臨界值進行了詳細的研究．

對稱雙方勢阱；雙模近似； 玻色約瑟夫森結(jié)； 宏觀量子自俘獲

自玻色愛因斯坦凝聚體(BECs)實現(xiàn)以來，理論[1-9]和實驗[10-11]都對其進行了廣泛的研究．雙勢阱模型作為一個簡單的物理模型，主要研究量子干涉效應(yīng)和約瑟夫森效應(yīng)的典型結(jié)構(gòu)．早在1997年Smerzi等人就利用玻色約瑟夫森模型對雙勢阱中排斥相互作用原子的隧穿動力學(xué)進行了研究[1]，并在理論上得出了一種非線性現(xiàn)象：宏觀量子自俘獲(MQST)[2-4]．十年后在光晶格中人們觀察到了約瑟夫森振蕩[10]，在雙勢阱中觀察到了約瑟夫森振蕩和宏觀量子自俘獲[11]．這些工作都是基于用一束失諧的激光把一個磁諧振子勢從中一分為二，實現(xiàn)雙勢阱，這樣囚禁在其中的BECs就一分為二，適當(dāng)調(diào)節(jié)雙勢阱之間的勢壘高度就形成了兩團弱耦合BECs．李衛(wèi)東等人指出在修正的玻色約瑟夫森結(jié)模型中不僅存在線性耦合項而且存在非線性耦合項，并且發(fā)現(xiàn)在弱耦合，即滿足雙模近似的條件下，勢壘高度和寬度不發(fā)生變化的情況下隨著非線性增加非線性耦合項將達到甚至超過線性耦合項，從而影響雙勢阱中BECs的動力學(xué)性質(zhì)[7-8]．我們知道在實驗中可以通過調(diào)節(jié)Feshbach共振改變非線性強度[10-11]．同樣可以改變激光脈沖的強度來調(diào)節(jié)勢壘的高度和寬度，從而影響兩勢阱的耦合強度．本文在對稱雙方勢阱這一模型中，利用文獻[8]的方法,研究了勢壘寬度對對稱雙方勢阱中BECs非線性耦合強度及其動力學(xué)影響．

1 非線性雙模近似模型

零溫下對稱雙方勢阱中弱相互作用的BECs的波函數(shù)滿足一維GP方程：

(1)

其中V(x)是雙勢阱函數(shù)形式為：

(2)

(3)

將雙模近似代入GP方程可推導(dǎo)出修正的非線性雙模近似動力學(xué)方程：

(4)

其中R(..)表示括號中式子的實部．其中κ為線性耦合項，為非線性耦合項．

模糊層次分析法FAHP(Fuzzy Analytic Hierarchy Process)是在傳統(tǒng)的層次分析法基礎(chǔ)上,考慮到人們對復(fù)雜事物判斷的模糊性,引入模糊一致性矩陣的決策方法。模糊層次分析法很好地解決了AHP判斷矩陣的一致性問題。本文在具體采用模糊層次分析法確定各個路由度量權(quán)重系數(shù)之前,首先作如下假設(shè):

(5)

定義兩個勢阱中的粒子布局數(shù)差z(t)≡n1(t)-n2(t)和相位差φ(t)=θ2(t)-θ1(t)，則非線性雙模近似的動力學(xué)方程可寫為：

(6)

(7)

表1 線性耦合強度與非線性耦合強度隨勢壘寬度和非線性強度的變化表

取勢壘高度V0=200, 通過計算表(1)發(fā)現(xiàn)當(dāng)勢壘寬度b≥0.30，即兩勢阱間的耦合強度較弱，同時非線性參數(shù)非線性較強時，會出現(xiàn)κ?χ 即非線性耦合項占主導(dǎo)作用．此時方程(6)可以忽略κ簡化為：

(8)

系統(tǒng)的總哈密頓量為：

(9)

此式比方程(6)更簡潔．但此式必須在弱耦合和強非線性同時滿足的條件下才能滿足．對于雙勢阱模型所采用的雙模近似在耦合強度越弱的情況下，近似程度越高；同時非線性強度可以通過調(diào)節(jié)Feshbach共振改變．所以上式能很好地描述雙模近似模型下強相互作用的量子氣體．

2 勢壘寬度對動力學(xué)的影響

2.1 約瑟夫森振蕩

利用(6)式取勢壘高度V0=200,勢壘寬度分別取b=0.15,0.30，0.45．粒子布局數(shù)差初始值z(0)=0.01，初始相位差為φ(0)=0計算出其粒子布局數(shù)差隨時間變化圖像，及其對應(yīng)的相位差圖像．如圖1可知，粒子數(shù)布局數(shù)差平均值為零，所以BECs進行約瑟夫森振蕩．勢壘寬度增加，勢阱間耦合強度減弱，其他條件相同的情況下，勢壘寬度越大，雙勢阱系統(tǒng)中BECs做約瑟夫森振蕩的周期越長；相位差的變化周期與粒子數(shù)布局差周期規(guī)律一致．

圖1 不同勢壘寬度和η=5,40時，粒子布局數(shù)差和相位差隨時間演化Fig.1 The time evolution of particle population imbalance and phase difference for various barrier width and η=5,40

2.2 宏觀量子自俘獲

同樣利用公式(6)取勢壘高度V0=200，初始相位差為φ(0)=0，且勢壘寬度分別取b=0.15，0.30，0.45．與小幅振蕩不同的是這里選取粒子布局數(shù)差的初始值z(0)=0.1．計算出其粒子布局數(shù)差隨時間變化圖像及其對應(yīng)的相位差圖像．由圖2可以看出不同的勢壘寬度在其他條件一致的情況下出現(xiàn)的振蕩形式不同；勢壘寬度越大的越容易出現(xiàn)粒子數(shù)布局數(shù)差平均值為非零，即宏觀量子自俘獲，關(guān)于這一點文章后續(xù)還有說明．從圖2還可以看出當(dāng)粒子布局數(shù)差進入宏觀量子自俘獲相時，左右勢阱中波函數(shù)相位差變成了相位的持續(xù)增加．

圖2 不同勢壘寬度和η=5,40時，粒子布局數(shù)差和相位差隨時間演化Fig.2 The time evolution of particle population imbalance and phase difference for various barrier width and η=5,40

2.3 弱耦合強相互作用下的動力學(xué)

當(dāng)勢壘的寬度較寬且非線性較強時，對稱雙勢阱中BECs的隧穿動力學(xué)可以用(8)式描述，選取勢壘寬度b=0.30，非線性強度為η=80，來比較保留線性耦合項和忽略線性耦合項對其動力學(xué)的影響，驗證簡化式(8)的合理性．從圖3可以看出在上述條件下忽略線性耦合項對BECs的動力學(xué)影響較?。虼丝芍?8)式能較好地刻畫強相互作用下BECs的動力學(xué)行為．

圖3 φ(0)=0，z(0)分別為0.05,0.10時保留線性耦合項和忽略線性耦合項時，粒子布局數(shù)差和相位差隨時間演化對比Fig.3 The contrast time evolution images of particle population imbalance and phase difference between keeping the linear coupling term and ignoring it when φ(0)=0, z(0)=0.05 and 0.10

2.4 勢壘寬度對宏觀量子自俘獲的非線性臨界值的影響

宏觀量子自俘獲現(xiàn)象是玻色約瑟夫森結(jié)模型中的一種新奇的量子現(xiàn)象．對于給定的模型中當(dāng)初始相位差φ(0)=0，對應(yīng)的初始粒子布局數(shù)差臨界值為：

(10)

圖4 不同勢壘寬度下初始粒子數(shù)差臨界值隨非線性強度變化圖像Fig.4 The threshold of initial population imbalance varies with the change of nonlinear intensity under different barrier widths

因此改變勢壘寬度和非線性強度都會改變臨界值z(0)．由圖4可知，隨著勢壘寬度越大，在相同的粒子數(shù)差初始情況下，動力學(xué)從約瑟夫森振動向宏觀量子自俘獲轉(zhuǎn)變的非線性值越?。瑫r發(fā)現(xiàn)勢壘較寬的模型出現(xiàn)自俘獲的區(qū)間較大．

3 結(jié)論

在對稱雙方勢阱中，利用修正的約瑟夫森模型，系統(tǒng)地研究了勢壘寬度對線性耦合項及其動力學(xué)的影響．發(fā)現(xiàn)在b≥0.30弱耦合模型，強相互作用下，非線性耦合項將變成主要耦合項，修正的約瑟夫森雙模動力學(xué)模型(6)式可以簡化為(8)式．同時發(fā)現(xiàn)在相同非線性強度下，勢壘越寬，約瑟夫森振蕩的周期越大，宏觀量子自俘獲的初始臨界值越?。ㄟ^計算發(fā)現(xiàn)在弱耦合強相互作用下簡化的約瑟夫森模型能夠很好地描述雙方勢阱中BECs的動力學(xué)．

[1] SMERZI A, FANTONI S, GIOVANAZZI S,etal. Quantum coherent atomic tunnneling between two trapped Bose-Einstein condensates[J]. Phys Rev Lett, 1998,79(25):4950-4953.

[2] RAGHAVAN S, SMERZI A, FANTONI S,etal. Transitions in coherent oscillations between two trapped Bose-Einstein condensates[J]. Phys Rev A, 1999,60(3):620-624.

[3] THIAGO F VOSCONDI, FURUYA K. Dynamics of a Bose-Einstein condensate in a symmetric triple-well trap[J]. J Phys A: Math Theor, 2011,44(17):175301-175303.

[4] MAHMUD K W, KUTZ J N, REINHARDT W P. Bose-Einstein condensates in a one-dimensional double square well: Analytical solutions of the nonlinear Schr?dinger equation[J]. Phys Rev A, 2002,66(6):063607.

[5] LIN B, FU L B, YANG S P,etal. Josephson oscillation and transition to self-trapping for Bose- Einstein condensates in a triple-well trap[J]. Phys Rev A, 2007,75(3):033601.

[6] CHENG J, JING H, YAN Y J. Spin-mixing dynamics in a spin-1 atomic condensate coupled with a molecular[J]. Phys Rev A, 2008,77(R)(6):061604.

[7] LI W D. Stationary solutions of Gross-Pitaevskii equations in a double square well[J]. Phys Rev A, 2006,74(6):063612.

[8] JIA X Y, LI W D, LIANG J Q. Nonliner correction to the boson Josephson-junction model[J]. Phys Rev A, 2008,78(2):023613.

[9] XIE Q T, HAI W H. Coherent control of self-trapping of two weakly coupled Bose-Einstein condensates[J]. Phys Rev A, 2007,75(1):015603.

[10] MAHMUD K W, PENG H, REINHARDT W P. Quantum phase-space picture of Bose-Einstein condensates in a double well[J]. Phys Rev A, 2005,71(2):023615.

[11] ALBIEZ M, GATI R, FOLLING J,etal. Direct observation of tunneling and nonlinear self-trapping in a single bosonic josephson junction[J]. Phys Rev Lett, 2005,95(1):010402.

(編輯 陳笑梅)

The Effect of the Barrier Width on the Nonlinear Coupling and Dynamics of Weakly Coupling Bose-Einstein Condensates

LIUXin-jian*,LIWei-dong

(Institute of Theoretical Physics, Shanxi University, Taiyuan 030006, China)

By using analytical and numerical methods, the effects of barrier width on the strength of nonlinear coupling and the dynamics of Bose-Einstein Condensates with symmetrical double square well were investigated. The results show that the linear coupling could be very small once the barrier width is large. When the width is larger than 0.3 and nonlinear parameters larger, the linear coupling can be safely neglected. In this case, the dynamics character of Bose-Jsoephson model is determined only by nonlinear coupling strength. Furthermore, the effects of this width on the period of the nonlinear Josephson oscillation and the nonlinear threshold value in macroscopic quantum self-trapping were presented in details.

symmetric double square well potential; two mode approximation; Bose-Josephson junction; macroscopic quantum self-trapping

2014-02-24

國家自然科學(xué)基金資助項目(11374197；11074155；10934004)

O469

A

1000-2537(2014)03-0053-05

*通訊作者,E-mail：15364839097@163.com