黃 芳 李海彬

(浙江工業(yè)大學(xué)理學(xué)院應(yīng)用物理系，杭州 310023)

雙勢阱中玻色-愛因斯坦凝聚的絕熱隧穿*

黃 芳 李海彬

(浙江工業(yè)大學(xué)理學(xué)院應(yīng)用物理系，杭州 310023)

(2010年2月9日收到;2010年5月16日收到修改稿)

研究了玻色-愛因斯坦凝聚體在雙勢阱中隨著能級差絕熱循環(huán)變化而發(fā)生的絕熱隧穿.發(fā)現(xiàn)當(dāng)相互作用較強且初態(tài)選擇為凝聚體全部置于較淺勢阱時，演化過程破壞絕熱定理，而演化結(jié)果有可能回到初態(tài)，也有可能不回到初態(tài)，取決于演化周期的選擇;另外還發(fā)現(xiàn)演化過程表現(xiàn)出對初態(tài)選擇的依賴，具有不對稱的特征.利用能級圖和相圖，對上述現(xiàn)象給出了解釋.

玻色-愛因斯坦凝聚，Landau-Zener模型，絕熱隧穿

PACS:03.65.-w，03.75.Lm，42.50.-p

1.引 言

1995年在實驗中實現(xiàn)了超冷稀薄原子氣體的玻色-愛因斯坦凝聚 (Bose-Einstein condensates，BEC)，激起了物理學(xué)界的研究熱潮，不論從實驗還是理論角度都取得了豐碩的成果.一方面，關(guān)于凝聚體的超流性質(zhì)，集體激發(fā)，相干性和量子相變等性質(zhì)的研究對于深入理解物理世界具有重要的意義［1］.另一方面，BEC的研究已經(jīng)顯著地滲透到物理學(xué)的其他各個領(lǐng)域，如量子通信、量子計算與量子信息等領(lǐng)域;而且，在原子鐘制造、原子干涉儀和重力常數(shù)精確測量等方面也有廣泛的應(yīng)用.

在關(guān)于凝聚體自身性質(zhì)的研究中，宏觀量子效應(yīng)是一個重要課題，例如宏觀量子隧穿現(xiàn)象，考慮放置在雙勢阱中的凝聚體，理論上首先發(fā)現(xiàn)了BEC在雙勢阱中的約瑟夫森效應(yīng)［2，3］，若考慮 BEC原子間的相互作用，還發(fā)現(xiàn)了自陷獲(Self-trap)等奇特現(xiàn)象［2，3］，這些都已在實驗中得到了證實［4，5］，并引起廣泛的研究興趣［6—11］.

實際上，雙勢阱中BEC的性質(zhì)在平均場近似和雙模近似下，可以用非線性 Landau-Zener模型［12］來描述.我們知道，Landau-Zener模型是一個典型的兩能級系統(tǒng)，在量子力學(xué)中，它是一個基本的物理模型，同時在量子化學(xué)、碰撞理論、鈉磁鐵的自旋隧穿和量子計算等領(lǐng)域有廣泛的應(yīng)用.在Landau-Zener模型基礎(chǔ)上，如果考慮相互作用等效應(yīng)，就可以得到非線性 Landau-Zener模型.可以在其中發(fā)現(xiàn)不同于線性模型的效應(yīng)，如量子絕熱定理的破壞等［12—14］.另外可以考慮對此模型的相關(guān)參數(shù)進行調(diào)制，研究對模型動力學(xué)的影響［15—20］，并由此研究光子晶格中的非線性 Landau-Zener隧穿［21，22］，非線性三能級體系的絕熱 Landau-Zener隧穿［23］，光晶格勢阱中 BEC的 Landau-Zener隧穿行為［24］等問題.

本文主要研究在雙勢阱中BEC的絕熱隧穿動力學(xué)特性.我們發(fā)現(xiàn)由于相互作用的存在，當(dāng)對兩勢阱凝聚體的能級差進行絕熱循環(huán)變化時，系統(tǒng)或者回到初態(tài)，或者不回到初態(tài)，而這取決于周期的選擇;另外，還發(fā)現(xiàn)隧穿結(jié)果對初態(tài)選擇的依賴，表現(xiàn)出不對稱的特性.

2.模 型

2.1.BEC雙勢阱模型

我們考慮在簡諧勢阱中形成的BEC，可以利用紅失諧的激光入射勢阱中間，從而形成雙勢阱，如圖1所示.此系統(tǒng)的動力學(xué)可以用無量綱的雙模非線性薛定諤方程

在本文中，我們將研究凝聚體在雙勢阱之間的絕熱隧穿，即當(dāng)以非常慢的速度(滿足量子絕熱定理)改變某個參數(shù)，凝聚體經(jīng)歷這個過程之后最終的狀態(tài)，特別是當(dāng)考慮原子間相互作用存在時，演化結(jié)果會有什么特征.在辨析系統(tǒng)動力學(xué)特征時，需要指出的是將雙勢阱BEC模型用二能級系統(tǒng)來表示，實際上存在兩種過程.首先，這里二能級實際上是兩個勢阱中的基態(tài)能量，而一般二能級系統(tǒng)中能級和位置無關(guān)，凝聚體在這兩能級上的重新分配，稱之為躍遷;其次，由于存在雙勢阱，凝聚體可以從一個勢阱轉(zhuǎn)移到另外一個勢阱，稱之為隧穿.從下面的研究結(jié)果可以看到，這兩種過程可以是一致的，也可以各自發(fā)生.

圖1 雙勢阱中的BEC 圖上部為分布示意圖，圖下部為概率.(a)為初始時刻狀態(tài)，BEC凝聚體處于左阱(淺勢阱)中，(b)當(dāng)c=0或cv時，系統(tǒng)演化結(jié)束后凝聚體的狀態(tài)

2.2.BEC的絕熱循環(huán)演化

首先考慮能級差隨時間單向線性變化的情況，即γ=αt，t從-∞向+∞演化.當(dāng)c=0，初始時刻凝聚體全部放置在深度較淺的勢阱里，即為圖1(a)中左勢阱，演化過程中左邊的勢阱逐漸降低，而右勢阱逐漸抬高，最后右勢阱成為較淺的勢阱.根據(jù)線性Landau-Zener模型的理論，隧穿概率可以解析求得為，如果參數(shù)變化很慢，即α→0，則p=0，這意味著原來在左勢阱的凝聚體全部隧穿至右勢阱，始終留在較淺勢阱里，即能級間的躍遷沒有發(fā)生，滿足量子絕熱定理，如圖1(b)所示.若考慮原子之間的相互作用即 c≠0，當(dāng)cv時，演化結(jié)束后，凝聚體狀態(tài)如圖1(c)的模擬，這意味著量子絕熱定理的破壞［12，13］，在這種情況下，能級間的躍遷和空間隧穿同時發(fā)生.

根據(jù)量子力學(xué)，一個孤立的量子系統(tǒng)，經(jīng)歷參數(shù)絕熱循環(huán)演化會回到初態(tài)，只是相位會增加一個幾何相(berry phase).本文中，我們將考慮當(dāng)原子間的相互作用不能忽略，雙勢阱中的BEC凝聚體經(jīng)歷能級差絕熱循環(huán)變化時的演化特征.選擇了兩種演化方式作為比較的對象.

由前面的介紹，能級差的變化方式一般采取隨時間線性變化的形式即γ=αt.現(xiàn)在依舊考慮采取能級差線性變化的方式，但是對γ進行分階段掃描，在所選參數(shù)下，可以用函數(shù)表示為 γ=-αt(ti≤t≤t0)，γ=-2αt0+αt(t0≤t≤tf)，其中ti，tf的選擇應(yīng)保證演化初態(tài)和末態(tài)近似為系統(tǒng)本征態(tài)，并且t0-ti=tf-t0.

圖2 能級差分段線性絕熱循環(huán)演化時BEC凝聚體分布隨時間的變化曲線 初態(tài)為模式.參數(shù)取值為v= 0.1，α=0.0001.(a)c=0.0和c=0.05，圖中兩個取值的曲線重合(虛線)，c=0.2(實線)，以上三種情況，時間為t0-ti=100000;(b)取c=0.2所得的結(jié)果，時間為t0-ti=100030

依舊選擇初態(tài)為凝聚體全部位于較淺勢阱(圖1(a)中左勢阱)，數(shù)值計算結(jié)果顯示于圖2.當(dāng)c=0或cv情況下，在正向掃描過程，當(dāng)γ演化跨越γ=0發(fā)生能級轉(zhuǎn)換后，p≈0.77，此時系統(tǒng)處于非線性本征態(tài)，這意味著凝聚體同時發(fā)生能級間躍遷和空間隧穿，如圖1(c)所示，之后系統(tǒng)保持在此狀態(tài)進入逆向掃描，當(dāng)再次跨越γ=0發(fā)生能級轉(zhuǎn)換后，我們發(fā)現(xiàn)當(dāng)演化周期的選擇不同時，會出現(xiàn)兩種結(jié)果，一個近似為0.0，也就是凝聚體全部重新回到初始勢阱，這種情況似乎滿足量子絕熱定理，但是，需要指出的是只是首尾兩點滿足，中間過程卻遭到了破壞.另一個結(jié)果近似為0.23，即保持在非線性本征態(tài)上，此時，凝聚體不能回到初始狀態(tài)，量子絕熱定理遭到破壞.

但是我們注意到上述演化函數(shù)存在拐點，而在此拐點處，量子絕熱定理不滿足，所以考慮采取能級差變化形式為這種情況下絕熱條件為以 γ=γ0為能級差初始值，凝聚體全部布局在淺勢阱中為初始態(tài)開始演化，且使γ0的選擇足夠大以保證此初始態(tài)近似為系統(tǒng)的本征態(tài)，計算結(jié)果見圖3.可以發(fā)現(xiàn)和分段線性演化情況(圖2)有相同的結(jié)果，即對于c=0和cv，則出現(xiàn)演化過程絕熱定理破壞，且最后演化結(jié)果同樣不唯一.當(dāng)對周期做改變時，或者得到圖3(a)(實線)所示，回不到初態(tài);或者如圖3(b)所示，回到初態(tài);而且我們同樣發(fā)現(xiàn)p的數(shù)值同樣近似為0.0和0.23，與對 γ分階段掃描情況一致.也就是說，對能級差分階段掃描過程中雖然有拐點出現(xiàn)，不符合絕熱定理，但是得到了相同的結(jié)果.

圖3 能級差以余弦形式絕熱循環(huán)演化時BEC凝聚體分布隨時間的變化曲線 初態(tài)為模式.圖中取值為γ0= 4.0，v=0.1，(a)取T=120000，c=0.0和0.05，這兩個取值的曲線重合(虛線)，c=0.2(實線);(b)取T=120020，c=0.2的概率

另外，我們還發(fā)現(xiàn)，在上述兩種演化過程中，演化結(jié)果還依賴于初始態(tài)的選擇，如果選擇將凝聚體全部布局在較深勢阱，即作為初態(tài)，則不論 c取多大，凝聚體都會經(jīng)歷滿足絕熱定理的演化過程而回到初態(tài)，如圖4所示，這和作為初態(tài)的情況完全不同，體現(xiàn)出一種不對稱的特征.那么，如何解釋這些奇特的行為呢?

圖4 BEC凝聚體分布隨時間的變化曲線 初態(tài)為模式.(a)能級差分階段線性變化，其中 v=0.1，α= 0.0001，時間為t0-ti=100000，c分別為0.0，0.05，2.0，所得到的曲線一致;(b)能級差余弦變化，其中，T= 120000，γ0=4.0，v=0.1，c分別為0.0，0.05，2.0，所得到的曲線一致

2.3.分析與討論

我們在圖5中給出了模型(1)式的能級結(jié)構(gòu)圖，其中c>v.可以看到(見圖5(a))，在高能級上出現(xiàn)了一個環(huán)狀結(jié)構(gòu).

圖5 非線性Landau-Zener模型的能級圖 圖中取值為v=0.1，c=0.2.(a)為正向掃描過程的能級圖，(b)為正向掃描結(jié)束后逆向掃描過程的能級圖

圖6中均取c>v的情況，γ所取的值對應(yīng)于能級圖5(a)上的虛線a，b，c，d，e.若選擇為初態(tài)，γ是從正值向負值變化.首先看圖6(a)，當(dāng)γ=0.4時，出現(xiàn)對應(yīng)本征態(tài)的兩個橢圓點p1，p2，當(dāng)勢阱中凝聚體從系統(tǒng)的本征態(tài)開始演化，隨 γ減小，演化到圖5(a)的環(huán)狀結(jié)構(gòu)邊緣，此時劈裂出雙曲線點 p3.當(dāng)γ=0時，兩勢阱能級發(fā)生轉(zhuǎn)換，凝聚體原子發(fā)生空間隧穿和能級間躍遷，較多部分原子隧穿到另一個勢阱中，即保留在高能級，其他較少部分留在原來勢阱.之后系統(tǒng)演化至圖6(d)中橢圓點p1和雙曲線點p3融合消失.圖6(e)為環(huán)結(jié)構(gòu)消失后的演化態(tài).如圖5(b)所示，此時初態(tài)處于非本征態(tài)，大部分原子處于高能級勢阱中，小部分原子處于低能級勢阱中，此狀態(tài)演化至γ=-γ0然后開始逆向演化，由于此時系統(tǒng)處于非本征態(tài)，p≈0.77，當(dāng)演化至 γ= 0，勢阱的能級位置再一次發(fā)生轉(zhuǎn)換時，一種可能是各勢阱中凝聚體原子不發(fā)生能級躍遷，依舊停留在這個非本征態(tài)上，此時p≈0.23，如圖3(a)(實線)所示，即在逆向掃描過程中遵從絕熱定理.另一種可能性為低能級勢阱中那小部分的原子全部躍遷到另一個能級上，得到圖3(b)p=0.0，即又回到了初態(tài).從計算中可以看出是否回到初態(tài)依賴于周期的選擇，這是因為在能級轉(zhuǎn)變的時候，模型實際上并沒有嚴格停留在一個狀態(tài)上，而是有小幅的振蕩，如圖3所示.當(dāng)周期不同，會導(dǎo)致在轉(zhuǎn)變點(γ= 0)系統(tǒng)的狀態(tài)，即振幅和相位，都有所不同，從而導(dǎo)致了演化結(jié)果的不同.

圖6 取c=0.2，v=0.1時Hamilton量在不同參數(shù)γ下的相空間結(jié)構(gòu)圖 (a)γ=0.4，(b)γ=0.045，(c)γ=0，(d)γ=-0.045，(e)γ=-0.4.箭頭代表不動點pi隨著γ減小的移動方向

3.結(jié) 論

通過緩慢改變雙勢阱中BEC凝聚體的能級差，使此系統(tǒng)產(chǎn)生循環(huán)演化來研究凝聚體的絕熱隧穿.發(fā)現(xiàn)當(dāng)原子相互作用達到一定強度時，初態(tài)全部放置在淺勢阱(高能級)中的凝聚體所經(jīng)歷的演化過程會破壞量子絕熱定理，發(fā)生了凝聚體的空間隧穿和能級躍遷，而且演化結(jié)束后，系統(tǒng)可能回到初態(tài)，也可能不回到初態(tài)，而這依賴于演化周期的選擇.另外，我們還發(fā)現(xiàn)，隧穿結(jié)果還依賴初態(tài)的選擇，當(dāng)從淺勢阱演化會出現(xiàn)絕熱隧穿的破壞，而從深勢阱(低能級)出發(fā)則不會，體現(xiàn)出不對稱特征.通過給出此模型的能級圖及相圖，我們發(fā)現(xiàn)能級中出現(xiàn)的環(huán)狀結(jié)構(gòu)改變了本征態(tài)的演化特征，從而導(dǎo)致了上述破壞絕熱定理的結(jié)果和初態(tài)依賴的不對稱特征.

［1］Dalfovo F，Giorgini S，Pitaevskii L P，Stringari S 1999 Rev. Mod.Phys.71 463

［2］Smerzi A，F(xiàn)antoni S，Giovanazzi S，Shenoy S R 1997 Phys. Rev.Lett.79 4950

［3］Raghavan S，Smerzi A，F(xiàn)antoni S，Shenoy S R 1999 Phys.Rev. A 59 620

［4］Albiez M，Gati R，F(xiàn)olling J，Hunsmann S，Cristiani M，Oberthaler M K 2005 Phys.Rev.Lett.95 010402

［5］Cataliotti F S，Burger S，F(xiàn)ort C，Maddaloni P，Minardi F，Trombettoni A.Smerrzi A，Inguscio M 2001 Science 293 843

［6］Shin Y，Saba M，Schirotzek A，Pasquini T A，Leanhardt A E，Pritchard D E，Ketterle W 2004 Phys.Rev.Lett.92 150401

［7］Zollner S，Meyer H D，Schmelcher P 2008 Phys.Rev.A 78 013621

［8］Fang Y C，Yang Z A，Yang L Y 2008 Acta Phys.Sin.57 661 (in Chinese)［房永翠、楊志安、楊麗云 2008物理學(xué)報 57 661］

［9］Adhikari S K，Lu H，Pu H 2009 Phys.Rev.A 80 063607

［10］Li G Q，Chen H J，Xue J K 2010 Acta Phys.Sin.59 1449(in Chinese)［李高清、陳海軍、薛具奎2010物理學(xué)報59 1449］

［11］Ma Y，F(xiàn)u L B，Yang Z A，Liu J 2006 Acta Phys.Sin.55 5623 (in Chinese)［馬 云、傅立斌、楊志安、劉 杰 2006物理學(xué)報55 5623］

［12］Wu B，Niu Q 2000 Phys.Rev.A 61 023402

［13］Liu J，F(xiàn)u L B，Ou B Y，Chen S G，Choi D，Wu B，Niu Q 2002 Phys.Rev.A 66 023404

［14］Liu J，Wu B，Niu Q 2003 Phys.Rev.Lett.90 170404

［15］Wang G F，Ye D F，F(xiàn)u L B，Chen X Z，Liu J 2006 Phys.Rev. A 74 033414

［16］Ishkhanyan A，Mackie M，Carmichael A，Gould P L，Javanainen J 2004 Phys.Rev.A 69 043612

［17］Duan S Q，F(xiàn)u L B，Liu J，Zhao X G 2005 Phys.Lett.A 346 315

［18］Ye D F，F(xiàn)u L B，Zhao H，Liu J 2007 Acta Phys.Sin.56 5071 (in Chinese)［葉地發(fā)、傅立斌、趙 鴻、劉 杰2007物理學(xué)報56 5071］

［19］Nesterenko V O，Novikov A N，Cherny A Y，Cruz F F S，Suraud E 2009 J.Phys.B:At.Mol.Opt.Phys.42 235303

［20］Wu L H，Duan W S 2009 Chin.Phys.B 18 4110

［21］Wang S，Yang Z A 2009 Acta Phys.Sin.58 3699(in Chinese)［王 沙、楊志安2009物理學(xué)報58 3699］

［22］Wang S，Yang Z A 2009 Acta Phys.Sin.58 729(in Chinese)［王 沙、楊志安2009物理學(xué)報58 729］

［23］Wang G F，Liu B，F(xiàn)u L B，Zhao H 2007 Acta Phys.Sin.56 3733(in Chinese)［王冠芳、劉 彬、傅立斌、趙 鴻2007物理學(xué)報56 3733］

［24］Xi Y D，Wang D L，Ding J W，She Y C，Wang F J 2010 Acta Phys.Sin.59 3720(in Chinese)［奚玉東、王登龍、丁建文、佘彥超、王鳳姣2010物理學(xué)報59 3720］

PACS:03.65.-w，03.75.Lm，42.50.-p

Adiabatic tunneling of Bose-Einstein condensate in double-well potential*

Huang Fang Li Hai-Bin
(Department of Applied Physics，Zhejiang University of Technology，Hangzhou 310023，China)

9 February 2010;revised manuscript

16 May 2010)

We studied the adiabatic tunneling of Bose-Einstein condensates(BEC)in double-well potential with circle adiabatic variation of energy level spacing.We found that the adiabatic theorem is violated during the evolution of BEC，when the initial state is chosen such that all condensate atoms are in the shallow well and the interaction between atoms is strong enough.Furthermore，we found that the condensates can return to the initial state or not，depending on the evolution periods.In addition，we also found the evolution is dependent on the choice of initial state，which indicates an asymmetric character.With the help of the graph of energy levels and the phase diagram，we explained the results obtained.

Bose-Einstein condensates，Landau-Zener model，adiabatic tunneling

*國家自然科學(xué)基金(批準(zhǔn)號:10605020)和浙江省自然科學(xué)基金(批準(zhǔn)號:R107377)資助的課題.

*Project supported by the National Natural Science Foundation of China(Grant No.10605020)and the Natural Science Foundation of Zhejiang Province，China(Grant No.R107377).