高帥 PREST Matthew ILO-OKEKE Ebubechukwu O. KONDAPPAN ARISTIZABAL-ZULUAGA Juan E. IVANNIKOV Valentin BYRNES Tim

摘要： 著重研究了一種在玻色-愛因斯坦凝聚態(tài)（Bose-Einstein Condensate， BEC）之間產(chǎn)生光學(xué)介導(dǎo)糾纏的方案. 該方案使用量子非破壞性哈密頓量 ， 并將 BEC 置于馬赫-曾德（Mach-Zehnder）構(gòu)型中. 結(jié)果表明 ， 通過對光進行測量 ， 可以誘導(dǎo)產(chǎn)生糾纏態(tài). 還特別分析了糾纏態(tài)在退相干作用下的效應(yīng). 研究表明 ， 該糾纏態(tài)的行為表現(xiàn)對于原子與光的作用時間較為敏感：當(dāng)相互作用時間?時 ， 該糾纏態(tài)相對穩(wěn)定;當(dāng)相互作用時間>時 ， 該糾纏態(tài)則相對脆弱.

關(guān)鍵詞：玻色-愛因斯坦凝聚態(tài);? 量子非破壞性測量;? 雙自旋糾纏態(tài);? 退相干

中圖分類號： O431.2??? 文獻標(biāo)志碼： A??? DOI：10.3969/j.issn.1000-5641.2022.02.011

Optically mediated entanglement between Bose-Einstein condensates

GAO Shuai1 ，? PREST Matthew1，2 ，? ILO-OKEKE Ebubechukwu O.2，3 ，? KONDAPPAN

Manikandan4???? ARISTIZABAL-ZULUAGA Juan E.1，2 ，? IVANNIKOV Valentin1，5，6，7 ，

BYRNES Tim1，5，6，7

（1. State Key Laboratory of Precision Spectroscopy， East China Normal University， Shanghai? 200241，China;2. Department of Physics， New York University Shanghai， Shanghai? 200122， China;3. Departmentof Physics， School of Science， Federal University of Technology， Imo State? 460001， Nigeria;4. Departmentof Precision Science and Natural Sciences， Institute of Physics， Universidad de Antioquia UdeA， Medellín05001000， Colombia;5. NYU-ECNU Institute of Physics at NYU Shanghai， Shanghai? 200062， China;6. National Institute of Informatics， Tokyo? 101-8430， Japan;7. Department of Physics， New YorkUniversity， New York， NY? 10003， USA）

Abstract： This paper explores a method for generating optically mediated entanglement between Bose- Einstein condensates （BECs）. Using a quantum nondemolition Hamiltonian with BECs placed in a Mach- Zehnder configuration， it is shown that entangled states can be induced by performing measurement on light. In particular， the effects of the entangled state in the presence of decoherence were analyzed. The behavior of the entangled state was found to be sensitive to the atom-light interaction time. The entangledstate is relatively stable when the dimensionless interaction time? ??? and relatively fragile when the time is greater.

Keywords： Bose-Einstein condensate;?? quantum nondemolition measurement;?? two-spin entangled state; decoherence

0? 引言

量子糾纏[1-2]作為一種重要的物理資源 ， 被廣泛地應(yīng)用于量子信息[3-5]和量子計算[6-7]等任務(wù)中.糾纏是量子力學(xué)的核心內(nèi)容之一 ， 因此如何產(chǎn)生多粒子糾纏態(tài)是一項重要的任務(wù). 而這其中的關(guān)鍵是多粒子糾纏態(tài)必須能被視作純單粒子能級 ， 玻色-愛因斯坦凝聚態(tài)（BEC）則恰好滿足這樣的條件. 最近人們越發(fā)對研究宏觀系統(tǒng)中的量子機理感興趣. BEC 是一種典型的介觀系綜 ， 它內(nèi)部包含103個或更多原子 ， 并且具有高度的相干性. 鑒于以上優(yōu)點 ， BEC 提供了一種在宏觀尺度上觀測糾纏的可能性. 最近有3 個研究小組[8-10]演示了 BEC 中的糾纏 ， 并證明了單個 BEC 的2 個空間區(qū)域之間存在糾纏. 但是 ， 目前還沒有 2個獨立的 BEC 之間糾纏的演示. 而本文所使用的光學(xué)介導(dǎo)方案則可以在 2個 BEC 之間產(chǎn)生糾纏 ， 因此本文方案具有重要的意義.

量子計量是原子系綜的主要應(yīng)用之一.使用自旋壓縮態(tài)可以提高測量靈敏度， 甚至可以超過標(biāo)準(zhǔn)量子極限.事實上 ， 單原子系綜中的量子相干態(tài)可以轉(zhuǎn)換成自旋壓縮態(tài)[11-12]. 根據(jù)海森堡不確定原理[13-15] ， 自旋壓縮態(tài)可以減小特定測量方向的不確定度， 反過來又可以增加另一個反方向的不確定度.實驗中， 已經(jīng)成功地將這種方法用在光學(xué)系統(tǒng)中[15-17]. 類似的 ， 這種方法也可以用在本文考慮的原子系統(tǒng)中.在測量內(nèi)部自旋能級時 ， 有許多不同的相互作用可以降低噪聲[12，18-20]. 在 BEC 中， 實驗中得到論證的有效壓縮作用有單軸和雙軸反扭轉(zhuǎn)哈密頓量[21-23] ， 這樣的相互作用廣泛應(yīng)用于量子計算 [24]、計量 [25]、密碼學(xué)[26] ， 以及引力波探測[27].

一般地 ， 真實的量子系統(tǒng)不可避免地與周圍環(huán)境產(chǎn)生相互作用 ， 此時量子系統(tǒng)會變得非常脆弱 ， 而且非常容易失去相干性 ， 這就是常說的退相干[28]. 這是實現(xiàn)量子信息處理和量子計算的重大障礙. 因此研究量子系統(tǒng)的退相干機制 ， 以及如何在最大程度上抑制退相干是量子信息學(xué)中難以忽視的問題 .基于上述原因 ， 本文研究了糾纏態(tài)在退相干機制下的魯棒性 ， 這將給相關(guān)量子邏輯門的設(shè)計提供可靠的量子比特來源 ， 并促進相應(yīng)技術(shù)的發(fā)展.

1? 物理模型

1.1? 理論模型

簡化的2 個 BEC 之間產(chǎn)生糾纏的物理實現(xiàn)裝置見圖 1. 圖1中 ， Coherent light 代表的是相干光源 ， Beam splitter 代表分束器 ， 光子數(shù)探測器則用 Photodetector 進行表示.該方案包含2個放置在馬赫-曾德干涉儀（置于相干光路中） 的不同臂的、空間分離的 BEC.光與 BEC 通過量子非破壞性 （Quantum Nondemolition， QND）哈密頓量相互作用. 隨后 ， 使用光電探測器對相應(yīng)的光子數(shù)進行探測 ， 最終得到一個結(jié)構(gòu)豐富的糾纏態(tài). 這里考慮的 BEC 位于2個分離良好的阱中.此處的阱可能是原子芯片上分離的2 個磁阱 ， 也可能是光學(xué)偶極阱[29-31]. 每個 BEC 都有2 個可以布居原子的內(nèi)部能級 ， 相應(yīng)的玻色湮滅算子為 gj ， ej ， 其中 j ∈{1， 2}. 常見的內(nèi)部能級態(tài)是超精細基態(tài) ， 例如87Rb 的 F =1， mF = 1和 F =2， mF = 2時鐘態(tài)[8，32-33].

為了推導(dǎo)的方便 ， 采用施溫格（Schwinger）玻色算符定義有效自旋 ， 即

其中 ， j ∈{1， 2}用來標(biāo)記 BEC.自旋算符遵循對易關(guān)系[Sl ， Sm]= 2i?lmnSn ， 這里 l， m， n ∈{x， y， z} ， x， y， z 表示3 個坐標(biāo)軸方向 ， ?lmn 是完全反對稱 Levi-Civita 張量.這里我們考慮每個 BEC 中的原子數(shù)是常量的情況 ， 因此有

（2）

這里原子數(shù)算符是?？耍‵ock）態(tài)的本征值 ， 可以記作

其中， Fock 態(tài)可以表示為

（4）

公式（4）中， |vac?是沒有原子和光子的真空態(tài). 通常用自旋相干態(tài)來描述 BEC， 其形式為

（5）

其中， 0? θ? π， ?π? φ? π是布洛赫（Bloch）球上的2 個任意球坐標(biāo)角. 具體推導(dǎo)過程如下.

2種不同光學(xué)模式的光束和各自的 BEC 相互作用.模式 b 是激光最初照射光束的玻色湮滅算符.該模式進入分束器 ， 進而轉(zhuǎn)換成

通過2 個 BEC 的光束的玻色湮滅算符是 a1 ， a2. 光學(xué)模式的斯托克斯（Stokes）算符可以表示成

（7）

其中 np 是光子數(shù)算符.在光和 BEC 干涉后 ， 第二個分束器進一步對光學(xué)模式進行轉(zhuǎn)換 ， 最后得到的相應(yīng)光學(xué)模式的玻色湮滅算符 ， 用 c， d 進行表示. 該轉(zhuǎn)換關(guān)系為

（8）

光子數(shù)算符的本征態(tài) ， 可以表示為

（9）

其中 b 是相應(yīng)光束的光學(xué)模式的湮滅算符.

光的初始態(tài)可以用光學(xué)相干態(tài)表征 ， 為

（10）

原子和光通過量子非破壞性（QND）哈密頓量相互作用 ， 該哈密頓量形式為

其中 ， Jz是描述光源的斯托克斯算符 ， S 和 S 是描述原子系統(tǒng)的施溫格玻色算符.注意到[H，Jz]= 0 ， 這說明Jz沒有時間依賴 ， 該哈密頓量確實是非破壞性的. 最后對光子數(shù)進行測量.

1.2? QND 糾纏方程

對于 QND 誘導(dǎo)糾纏的分析 ， 通常采用 Holstein-Primakoff （HP）近似

本文采用精確推導(dǎo)的方法 ， 這將允許用更長的時間尺度研究其行為模式. 這里考慮 BEC 的通用初始態(tài)

一種合適的初始態(tài)是2個偏振在Sx 方向的 BEC， 將其分別用自旋相干態(tài)表示 ， 易得

（14）

將其在 Fock 基中進行展開 ， 可以得到

（15）

激光將光學(xué)模式 b 映射成光學(xué)相干態(tài) ， 因此光學(xué)初始態(tài)可表示為|α?b .在經(jīng)過第一個分束器后， 初始波函數(shù)變?yōu)?/p>

其中 ， a1，a2是光學(xué)模式 b 經(jīng)過分束器后所產(chǎn)生的新光學(xué)模式.隨后光和 BEC 產(chǎn)生相互作用 ， 假設(shè)該相互作用時間為 t， 則相互作用后的初始態(tài)演化為

|?（τ）?= e ?? |?（0）?= exp [?? （S ? S ）Jz]|?（0）?.???? （17）

這里定義無量綱時間τ≡ ?t .將式（13）帶入式（16）， 得到

（18）

下一步應(yīng)用公式（8）的分束算符， 初始態(tài)演變?yōu)?/p>

最后， 分別探測在光學(xué)模式c， d 下的光子數(shù) ， 這等價于將式（19）用測量算符進行投影 ， 有

（20）

忽略不相關(guān)的全局相位 ， 未歸一化態(tài)為

（21）

其中， Fock 基下的系數(shù)Cnc nd （χ）具有

該式可由

（23）

進行歸一化.

公式（21）給出的形式表明 ， 該光介導(dǎo)方案的效應(yīng)是非常簡單的 ， 最初的波函數(shù)被一個額外因子 Cnc nd? [（ k1 ? k2）τ] 所調(diào)制. 在測量光子數(shù)后 ， 歸一化糾纏態(tài)的波函數(shù)為

（24）

這里歸一化因子等價于光子探測概率

（25）

公式（24）是我們最終求得的2 個 BEC 之間的波函數(shù). 顯然 ， 公式（19）的形式表明 ， 初始光學(xué)態(tài)的取值可以是任意的. 精確波函數(shù)表達式（21）適用于任意時間 ， 從而大大擴展了模型的適用范圍.

1.3? 近似分析

在獲得了波函數(shù)精確表達式后 ， 為了直觀研究該糾纏態(tài)的性質(zhì) ， 我們對其進行了近似分析. 這里研究較短的相互作用時間τ?? ， 因為該區(qū)域是產(chǎn)生 BEC 間糾纏的最相關(guān)區(qū)域 ， 且適合近似分析.

首先， 當(dāng) N ?1時 ， 可以將二項式系數(shù)近似為

（26）

可以看出 ， 只有在 ?? ? k1， k2 ?? +范圍內(nèi)的 Fock 項 ， 其系數(shù)才有明顯的權(quán)重 ， 這也導(dǎo)致了波函數(shù)（21）系數(shù)產(chǎn)生了不均勻疊加.

如果總光子數(shù) Np ? 1 ， 可以進一步對式（21）做斯特林（Stirling）近似[34] ， 可以得到高斯（Gauss）型函數(shù)

（27）

其中nc 是探測到的 c 光學(xué)模式下的光子數(shù). 在較短時間內(nèi) ， 事實上 ， 光子數(shù)大多集中在模式 c 上 ， 所以我們選取參數(shù)nc = Np， nd = 0. 從公式（27）可以明顯得到 ， 最大值發(fā)生在

（28）

由于余弦函數(shù)的振蕩性質(zhì) ， 它通常具有k1 ， k2的多個解.

事實上 ， 還可以直接對余弦函數(shù)進行指數(shù)近似.此時不需要限定 Np ? 1 ， 只需保證時間較小即可 ， 便得到這樣的近似

（29）

綜合公式（26）和式（29）， 最終公式（21）的波函數(shù)近似為

（30）

該近似對于 N ?1， |τ| ?? ， |α|? 1有效.這里 ， N 是歸一化系數(shù) ， 為

此近似的概率分布為

其中， k1， k2是 2個 BEC 在 Fock 態(tài)下的表示模式.從公式（32）可以看到， 糾纏脈沖在k1， k2兩個模式間產(chǎn)生了關(guān)聯(lián). 公式（32）的第一項是個對稱的非關(guān)聯(lián)高斯分布 ， 平均值?k1， k2?=? ， 標(biāo)準(zhǔn)差為; 第二項則產(chǎn)生了關(guān)聯(lián)并導(dǎo)致了最后的糾纏. 也可以看到 ， k1 = k2的情況概率分布最明顯 ， 且關(guān)聯(lián)的強度和無量綱的時間τ以及光子數(shù)nc 有關(guān). 為了產(chǎn)生強關(guān)聯(lián) ， 式（30）表明， 相互作用時間τ應(yīng)該滿足

（33）

聯(lián)想到時間τ?? ， 可以得到關(guān)系式

所以， 這說明使用明亮的相干光 ， 可以幫助產(chǎn)生 BEC 之間的更強關(guān)聯(lián).

2? 結(jié)果與分析

在得到雙自旋糾纏態(tài)的波函數(shù)表達式（21）后 ， 對其概率密度分布進行研究 ， 可以對其性質(zhì)有更好的認識. 另一方面 ， 通過施加馬爾可夫退相干作用 ， 研究該糾纏態(tài)在不同時間尺度下的行為特征及魯棒性， 這將有助于其在量子信息以及未來的量子計算方面的應(yīng)用.

2.1? 概率密度分布

圖 2給出了不同時間下的雙自旋糾纏態(tài)（公式（21））的概率密度分布.橫坐標(biāo)k1是第一個 BEC 在 Fock 態(tài)下的表示模式 ， 縱坐標(biāo)k2是第二個 BEC 在 Fock 態(tài)下的表示模式.通過概率分布 ， 可以直觀地分析2 個 BEC 之間的關(guān)聯(lián)狀態(tài). 所采用的參數(shù)：原子數(shù) N =20; 光子數(shù)nc = 20; 光子數(shù)nd = 0; α= ?? ?=? .

從圖 2（a）可以看到 ， 在τ= 0時 ， 概率分布呈現(xiàn)高斯分布 ， 且中心處于k1 = k2 =? .這是非常容易理解的：考慮到當(dāng)τ= 0時 ， 沒有施加相互作用 ， 此時就是兩個自旋相干態(tài)（式（14））.該糾纏態(tài)繼續(xù)演化 ， 直到時間在圖 2（b）的τ=? 尺度時 ， 明顯地觀察到壓縮現(xiàn)象 ， 方向沿著k1 = k2 ， 這完全符合公式（32）所表達的內(nèi)容.事實上 ， 圖 2（b）的概率分布和雙模壓縮態(tài)的概率分布是一致的 ， 這說明初始態(tài)經(jīng)過τ=? 的演變后 ， 已演變成雙模壓縮態(tài). 繼續(xù)增加時間τ ， 當(dāng)τ= 時 （圖 2（c））， 此時觀察到更大程度的壓縮和更強的關(guān)聯(lián) ， 并在主關(guān)聯(lián)峰旁邊出現(xiàn)了兩個副關(guān)聯(lián)峰. 當(dāng)時間來到圖 2（d）的τ= 時 ， 出現(xiàn)了更多的多關(guān)聯(lián)峰 ， 這意味著波函數(shù)出現(xiàn)了非高斯行為模式 ， 同時也意味著 Holstein-Primakoff （HP）近似的瓦解. 這足以體現(xiàn)本文采用的精確推導(dǎo)的方法的優(yōu)越性. 因為我們可以在更長的時間下研究態(tài)的演化. 全面對比圖 2的4 幅圖 ， 不難看出 ， 單關(guān)聯(lián)峰和多關(guān)聯(lián)峰分布的分界線是在τ=? ， 即τ～級別上.

事實上 ， 之所以會出現(xiàn)多關(guān)聯(lián)峰 ， 這是由于公式（28）中三角函數(shù)的周期性引起的：由于 n 取值的多樣性（n =0， ±1， ±2， ±3， ·· ·）， 因此可以有多組 k1 ? k2滿足條件. 而對于時間分界線的出現(xiàn) ， 這可以從公式（26）得到解釋：由于只有在 ?? ? k1， k2 ?? +范圍內(nèi)的 Fock 項 ， 其系數(shù)才有明顯的權(quán)重.所以此時便有k1 ? k2 ～? .從公式（28）容易看出， 只有當(dāng)時間τ～? ， 才會出現(xiàn)2τ（k1 ? k2）～nπ ， 此時多關(guān)聯(lián)峰效應(yīng)才開始顯著. 而在此之前 ， 則由于時間尺度太小 ， 超出了有效 Fock 態(tài)范圍 ， 因此呈現(xiàn)單關(guān)聯(lián)峰分布.

2.2? 馬爾可夫退相干

近似態(tài)（式（30））直觀展示出了2 個 BEC 之間的糾纏 ， 然而這是沒有考慮任何耗散的理想情況. 接下來 ， 將要探究馬爾可夫退相干下糾纏態(tài)的演化 ， 該效應(yīng)會顯著影響原子系綜里糾纏的大小.抑制退相干影響的能力 ， 是衡量糾纏方案價值的有力標(biāo)準(zhǔn).在本文所述系統(tǒng)中 ， 退相干的物理來源是交流斯塔克（AC-Stark）效應(yīng)和文獻[35]所述的原子阱電流擾動. 退相干存在于任何系統(tǒng)中 ， 初始態(tài)密度矩陣可以寫成

（35）

作為一個通用的退相干模型 ， 對每個自旋態(tài)施加馬爾可夫退相干

（36）

這里考慮 j = z， x 的效應(yīng). Sz 退相干和 Sx 退相干是2種不同的退相干機制. 事實上 ， 由于林德布拉德（Lindblad）算符Sz 和 QND 哈密頓量對易 ， 我們可以在相互作用產(chǎn)生糾纏態(tài)后 ， 直接施加退相干效應(yīng) ， 這將仍然得到精確表達式.因此 ， 將馬爾可夫退相干直接應(yīng)用到波函數(shù)表達式（34）， 最終獲得Sz 退相干精確密度矩陣 ， 其表達式為

其中， N1 =? nc nd （τ） | nc nd （τ）?.

而對于Sx 退相干 ， 精準(zhǔn)的求解方式是解公式（36）.但是在較短的時間內(nèi) ， 可以做一種近似分析 ， 即在獲得波函數(shù)精確表達式（21）后 ， 在布洛赫（Bloch）球上對其進行旋轉(zhuǎn)操作 ， 然后施加Sz退相干作用.

事實上 ， 從一個（ S ? S ）哈密頓量開始計量糾纏 ， 最后施加Sx 退相干 ， 和從（ S ? S ）哈密頓量開始然后施加Sz退相干的方式是等價的.在Sy 基中 ， 對于公式（21）做一個角度的旋轉(zhuǎn) ， 最后乘上退相干表達式 e ?2Γτ [（k1?k ）2+（k2?k ）2]便可以得到 Sx退相干下的波函數(shù) ， 也即

（38）

其中， N2 =??k ，k? （τ） |?k ，k? （τ）?是歸一化因子. ?k|e ?iSy?? |k′?由

（39）

給出.

最后， 得到Sx退相干下的密度矩陣 ， 其表達式為

（40）

圖 3 給出了施加馬爾可夫退相干作用后， 不同退相干因子下， 歸一化糾纏值與時間的關(guān)系 ， 其中 Emax 由公式（42）給出. 所采用參數(shù)： 原子數(shù) N =20; 光子數(shù) nc = 20; 光子數(shù) nd = 0; α=? =? .

一般來說 ， 退相干會急速地破壞現(xiàn)有的糾纏. 對于純態(tài)來說 ， 馮·諾依曼熵（von Neumann entropy）可以很好地量化糾纏. 由于退相干的存在 ， 原有的糾纏態(tài)演變成混合態(tài) ， 更好的量化方式是對數(shù)負性[36-37]. 糾纏值的計算表達公式為

其中， ρT2是第二個 BEC 的部分轉(zhuǎn)置. 可將糾纏值用最大糾纏值進行歸一化 ， 即

（42）

在圖3中 ， 在Γ =0 時 ， 糾纏展示出了一種魔鬼裂縫（devil’s crevasse）， 這和采用馮·諾依曼熵量化的糾纏所表現(xiàn)的行為相似.在文獻[38]中 ， 糾纏展示出了分形擾動 ， 且糾纏結(jié)構(gòu)圖的周期為 .考慮到哈密頓量的相似性 ， 這可以非常容易理解.注意到一個典型的時間τ=? ， 在該時間點之前 ， 糾纏線性增加 ， 在該時間點后 ， 糾纏展現(xiàn)出復(fù)雜的谷狀結(jié)構(gòu). 另一個特點是糾纏天花板 ， 這是由于二項式項不均勻疊加產(chǎn)生的 ， 它只包括中心在 k1， k2 =? 附近～ O（） 范圍的 Fock 項. 因此限制了糾纏值并產(chǎn)生了糾纏天花板. 此外， 由于退相干效應(yīng) ， 歸一化糾纏值隨著退相干因子增加而變小.

為了研究糾纏的特性 ， 這里選擇了3 個典型的時間τ=? ，? ，?? . 圖4給出了Sz 和Sx 退相干方式下 ， 不同時間尺度糾纏值 E （公式 （41））與原子數(shù)的關(guān)系 .所采用參數(shù)：最大原子數(shù) Nmax =20;光子數(shù)nc = 20; 光子數(shù)nd = 0; α=? =? ;Γz，x = 0.1. 圖 5則給出了Sz 和Sx 退相干方式下 ， 不同時間尺度糾纏比 R （公式（43））與原子數(shù)的關(guān)系圖. 所采用參數(shù)：最大原子數(shù)Nmax =20; 光子數(shù) nc = 20; 光子數(shù)nd = 0; α=? =? .結(jié)合圖 4和圖 5， 很明顯地可以觀察到 ， 3種時間下糾纏行為差異十分顯著.

在圖4中 ， 圖4（a）代表Sz 退相干方式 ， 而圖4（b）代表Sx 退相干方式. 對于τ=? ， 糾纏值 E 與Γj = 0的情況對于原子數(shù) N具有相似的依賴. 并且在 N →∞ 時 ， 糾纏比

（43）

趨向于1個常數(shù)（圖5（b））.隨著 N 增大， 由于相互作用門減小時間∝? .因此， 退相干能夠作用的時間也隨之減小 ， 兩者得以抵消.

在τ=? 時 ， 圖 4中的糾纏逐漸增加.在圖5中 ， 對于Sz 退相干（ 圖5（a））， 糾纏比基本上是恒定的， 而對于Sx 退相干（ 圖5（b））， 在 N →∞ 時 ， 糾纏比緩慢衰減到 R =0 .

關(guān)于時間τ= 的Sx 退相干 ， 糾纏迅速以～ e ?N2依賴衰減（ 圖4 （b））.對于Sx 退相干 ， 糾纏比也以指數(shù)形式衰減（圖5（b））.這種高斯依賴是薛定諤貓態(tài)在退相干下的典型特征. 事實上 ， 在τ= 時 ， 波函數(shù)可以寫成

（44）

值得注意的是 ， 對于時間τ= 的Sz 退相干 ， 隨著原子數(shù) N 的增加 ， 糾纏比 R 先保持常量后增加.考慮到主方程（36）的退相干自由子空間 ， 可以很容易理解[39]. 這里退相只改變量子位的局部相位 ， 隨機相位并不會影響糾纏 ， 因為量子位之間的糾纏早已建立.

總結(jié)上面的分析：對于所有的3 個相互作用時間 ， 在施加Sz退相干效應(yīng)時 ， R 和 N 有這樣的關(guān)系其中， 0? q <1 ;而對于Sx 退相干， 3種時間尺度表現(xiàn)出不同的行為模式

其中 ， 0< q <1 ， a =0.95569， b =0.790911， c =2.18937， d =0.293331， h =4.26138. 這表明如果退相干和哈密頓量在對易基上相互作用（Sz 退相干）， 2個 BEC 之間的糾纏可以足夠穩(wěn)定 ， 即使 N 非常大 ， 依舊如此. 綜合兩種退相干方式 ， 對于較短的時間尺度τ～? ， 糾纏對于退相干表現(xiàn)得足夠穩(wěn)定. 對于更長的時間尺度τ～ O （1） ， 由于薛定諤貓態(tài)的出現(xiàn) ， 糾纏態(tài)會變的非常脆弱.這兩種截然不同行為的分界時間線是τ～? .這說明 ， 當(dāng)τ?時 ， 糾纏態(tài)足夠穩(wěn)定并且適合實驗觀測.

可以看到 ， 對于Sz 退相干和Sx 退相干 ， 其行為表現(xiàn)存在著顯著差異. 這的確是2 種不同的退相干方式.上述已經(jīng)介紹了2 種退相干方式密度矩陣來源的不同.這里有必要其物理來源做出進一步的解釋 .事實上 ， 這是源于不同的自旋算符對于 Fock 態(tài)具有不同的作用形式 ， 即

（47）

此處 Fock 態(tài)是 Sz算符的本征態(tài). 容易看到 ， Sz算符不會對 Fock 態(tài)有移位作用 ， 而Sx退相干則會對 Fock 態(tài)產(chǎn)生移位作用（|k +1? ， |k ?1?）， 因此必然會造成退相干密度矩陣的不同. 這也最終導(dǎo)致了 Sx 退相干和Sz 退相干下， 該雙自旋糾纏態(tài)的不同行為表現(xiàn).

3? 結(jié)論

本文主要研究了玻色-愛因斯坦凝聚態(tài)間的光學(xué)介導(dǎo)糾纏 ， 在馬赫曾德干涉儀光路下， 利用光和原子的干涉 ， 通過精確推導(dǎo) ， 得到了具體的波函數(shù). 并進一步研究了其在不同時間下的概率分布.發(fā)現(xiàn)在較短的時間范圍內(nèi) ， 該糾纏態(tài)呈現(xiàn)明顯的壓縮效應(yīng) ， 而當(dāng)時間τ= 時 ， 出現(xiàn)了很多的關(guān)聯(lián)峰 ， 這表明非高斯行為的出現(xiàn). 這種現(xiàn)象出現(xiàn)的原因 ， 是由于公式（28）中三角函數(shù)的周期性導(dǎo)致的.且出現(xiàn)單關(guān)聯(lián)峰和多關(guān)聯(lián)峰分布的分界線是在 τ= 即 τ～級別上 .這是由于二項式系數(shù)限制了有效Fock 態(tài)范圍 ? ?k1 ， k2 ?? + ， 同時在公式（28）的多關(guān)聯(lián)峰公式的雙重作用下 ， 導(dǎo)致了只有τ ～1/才開始觀察到多關(guān)聯(lián)峰. 本文還研究了在馬爾可夫退相干作用下 ， 糾纏態(tài)在不同時間下的魯棒性. 發(fā)現(xiàn)對于較短的時間尺度τ～ 1/N ， 退相干對于糾纏態(tài)的影響較弱. 對于更長的時間尺度τ ～ O （1） ， 由于薛定諤貓態(tài)的出現(xiàn) ， 糾纏態(tài)會變得非常脆弱. 這兩種截然不同行為的分界時間線是τ ～? .對于Sz 退相干和Sx 退相干 ， 其行為模式差異顯著 ， 這是由不同的自旋算符對 Fock 態(tài)具有不同的作用形式而導(dǎo)致的.本文工作可為量子測量和量子計算提供一種可靠的糾纏來源 ， 并通過對退相干的探究， 促進相應(yīng)量子技術(shù)的發(fā)展.

[參考文獻]

[1] NIELSEN M A， CHUANG I L. Quantum Computation and Quantum Information [M]. 10th Anniversary ed. [S.l.]： CambridgeUniversity Press， 2011.

[2] HORODECKI R， HORODECKI P， HORODECKI M， et al. Quantum entanglement [J]. Reviews of Modern Physics， 2009， 81（2）：865.

[3] BENNETT C H， BRASSARD G， CRéPEAU C， et al. Teleporting an unknown quantum state via dual classical and Einstein-Podolsky-Rosen channels [J]. Physical Review Letters， 1993， 70（13）：1895-1899.

[4] EKERT A K. Quantum cryptography based on Bell’s theorem [J]. Physical Review Letters， 1991， 67（6）：661.

[5] BYRNES T， ROSSEAU D， KHOSLA M， et al. Macroscopic quantum information processing using spin coherent states [J]. OpticsCommunications， 2015， 337：102-109.

[6] BYRNES T， WEN K， YAMAMOTO Y. Macroscopic quantum computation using Bose-Einstein condensates [J]. Physical Review A，2012， 85（4）：4233-4237.

[7] GROVER L K. Quantum mechanics helps in searching for a needle in a haystack [J]. Physical Review Letters， 1997， 79（2）：325.

[8] FADEL M， ZIBOLD T， DéCAMPS B， et al. Spatial entanglement patterns and Einstein-Podolsky-Rosen steering in Bose-Einsteincondensates [J]. Science， 2018， 360（6387）：409-413.

[9] KUNKEL P， PRüFER M， STROBEL H， et al. Spatially distributed multipartite entanglement enables EPR steering of atomic clouds[J]. Science， 2018， 360（6387）：413-416.

[10] LANGE K， PEISE J， LüCKE B， et al. Entanglement between two spatially separated atomic modes [J]. Science， 2018， 360（6387）：416-418.

[11] BYRNES T， ILO-OKEKE E O. Quantum Atom Optics： Theory and Applications to Quantum Technology [M]. [S.l.]： CambridgeUniversity Press， 2021.

[12] GROSS C. Spin squeezing， entanglement and quantum metrology with Bose-Einstein condensates [J]. Journal of Physics B， 2012，45（10）：103001.

[13] S?RENSEN A， DUAN L M， CIRAC J I， et al. Many-particle entanglement with Bose–Einstein condensates [J]. Nature， 2001， 409（6816）：63-66.

[14] MACHIDA S， YAMAMOTO Y， ITAYA Y. Observation of amplitude squeezing in a constant-current –driven semiconductor laser [J].Physical Review Letters， 1987， 58（10）：1000-1003.

[15] WU L A， KIMBLE H J， HALL J L， et al. Generation of squeezed states by parametric down conversion [J]. Physical Review Letters，1986， 57（20）：2520-2523.

[16] SLUSHER R E， HOLLBERG L W， YURKE B， et al. Observation of squeezed states generated by four-wave mixing in an opticalcavity [J]. Physical Review Letters， 1985， 55（22）：2409-2412.

[17] BREITENBACH G， SCHILLER S， MLYNEK J. Measurement of the quantum states of squeezed light [J]. Nature， 1997， 387（6632）：471-475.

[18] MACOMBER J D， LYNCH R. Squeezed spin states [J]. The Journal of Chemical Physics， 1985， 83（12）：6514-6519.

[19] WALLS D F， ZOLLER P. Reduced quantum fluctuations in resonance fluorescence [J]. Physical Review Letters， 1981， 47（10）：709-711.

[20] WODKIEWICZ K， EBERLY J H. Coherent states， squeezed fluctuations， and the? SU（2） am? SU（1，1） groups in quantum-opticsapplications [J]. Journal of the Optical Society of America B， 1985， 2（3）：458-466.

[21] KITAGAwA M， UEDA M. Squeezed spin states [J]. Physical Review A， 1993， 47（6）：5138-5143.

[22] MUESSEL W， STROBEL H， LINNEMANN D， et al. Scalable spin squeezing for quantum-enhanced magnetometry with Bose-Einsteincondensates [J]. Physical Review Letters， 2014， 113（10）：103004.

[23] HALD J， S?RENSEN J L， SCHORI C， et al. Spin squeezed atoms： A macroscopic entangled ensemble created by light [J]. PhysicalReview Letters， 1999， 83（7）：1319-1322.

[24] NAVARRETE-BENLLOCH C . Quantum information with continuous variables [M]// An Introduction to the Formalism of QuantumInformation with Continuous Variables.[S.l]：[s.n.]， 2015.

[25] MOXLEY F， DOWLING J， DAI W， et al. Sagnac interferometry with coherent vortex superposition states in exciton-polaritoncondensates [J]. Physical Review A， 2016， 93：053603.

[26] HILLERY M. Quantum cryptography with squeezed states [J]. Physical Review A， 2000， 61（2）：022309.

[27] BONDURANT R S， SHAPIRO J H. Squeezed states in phase-sensing interferometers [J]. Physical Review D Particles & Fields， 1984，30（12）：2548-2556.

[28] BREUER H P， PETRUCCIONE F. The Theory of Open Quantum Systems [M].[S.l.]： Oxford University Press， 2002.

[29] REICHEL J， VULETIC V. Atom Chips [M].[S.l.]： John Wiley & Sons， 2011.

[30] WHITLOCK S， GERRITSMA R， FERNHOLZ T， et al. Two-dimensional array of microtraps with atomic shift register on a chip [J].New Journal of Physics， 2009， 11（2）：023021.

[31] ABDELRAHMAN A， MUKAI T， H?FFNER H， et al. Coherent all-optical control of ultracold atoms arrays in permanent magnetictraps [J]. Optics Express， 2014， 22（3）：3501-3513.

[32] B?HI P， RIEDEL M F， HOFFROGGE J， et al. Coherent manipulation of Bose –Einstein condensates with state-dependent microwavepotentials on an atom chip [J]. Nature Physics， 2009， 5（8）：592-597.

[33] RIEDEL M F， B?HI P， LI Y， et al. Atom-chip-based generation of entanglement for quantum metrology [J]. Nature， 2010， 464（7292）：1170-1173.

[34] ILO-OKEKE E O， BYRNES T. Information and backaction due to phase-contrast-imaging measurements of cold atomic gases： BeyondGaussian states [J]. Physical Review A， 2016， 94（1）：013617.

[35] LONE M Q， BYRNES T. Suppression of the ac-Stark-shift scattering rate due to non-Markovian behavior [J]. Physical Review A，2015， 92（1）：011401.

[36] VIDAL G， WERNER R F. Computable measure of entanglement [J]. Physical Review A， 2002， 65（3）：032314.

[37] PLENIO M B. Logarithmic negativity： A full entanglement monotone that is not convex [J]. Physical Review Letters， 2005， 95（9）：090503.

[38] BYRNES T. Fractality and macroscopic entanglement in two-component Bose-Einstein condensates [J]. Physical Review A， 2013，88（2）：023609.

[39] LIDAR D A， WHALEY K B. Decoherence-free subspaces and subsystems [M]// Irreversible Quantum Dynamics， Lecture Notes inPhysics， vol 622. Berlin： Springer， Berlin， 2003：83-120.

（責(zé)任編輯：李藝）