汪野 張靜寧 金奇奐

(清華大學(xué)交叉信息研究院量子信息中心, 北京 100084)

(2018 年 9 月 19 日收到; 2018 年 10 月 18 日收到修改稿)

能夠長時間儲存量子信息的量子存儲設(shè)備是實現(xiàn)大規(guī)模量子計算和量子通信的基本要素. 與其他量子計算平臺相比, 囚禁離子系統(tǒng)的優(yōu)勢之一在于具有很長的相干時間. 此前, 基于囚禁離子的單量子比特相干時間不到1 min. 研究發(fā)現(xiàn), 在囚禁離子系統(tǒng)中, 限制量子比特相干時間的主要因素是運(yùn)動能級加熱和環(huán)境噪聲, 其中后者包含環(huán)境磁場漲落和微波相位噪聲. 在同時囚禁171Yb+離子和138Ba+離子的混合囚禁系統(tǒng)中, 通過實施協(xié)同冷卻和動力學(xué)解耦, 可以實現(xiàn)相干時間超過10 min的單離子量子比特. 這一技術(shù)有望用于實現(xiàn)量子密碼學(xué)和搭建混合量子計算平臺.

1 引 言

量子計算和量子模擬[1]是近年來飛速發(fā)展的理論與實驗并重的熱點(diǎn)研究領(lǐng)域之一. 它以量子力學(xué)和量子信息學(xué)為理論基礎(chǔ), 以囚禁離子、超導(dǎo)電路等實驗技術(shù)的進(jìn)步為前提, 致力于實現(xiàn)實用的量子計算和量子通信[2,3], 具有廣闊的應(yīng)用前景. 量子計算平臺的基本元素是量子比特, 它是完全可控的量子二能級系統(tǒng), 人們可以對其進(jìn)行初始化、量子邏輯門和測量等基本操作. 量子力學(xué)中的不可克隆原理告訴我們, 儲存在量子比特上的未知量子信息不能被精確地復(fù)制[4,5]. 這項量子信息獨(dú)有的性質(zhì)只能夠通過對單個量子比特的相干操作和探測來驗證. 如果不能完全控制單個量子比特, 量子信息的儲存、操作和讀出從根本上是難以令人信服的. 長期以來, 人們一直致力于發(fā)展在單量子比特層面的相干操作和探測技術(shù).

囚禁離子量子比特的長相干時間將成為實際量子計算和量子通信的基本元素之一. 囚禁離子系統(tǒng)是實現(xiàn)大規(guī)模量子計算機(jī)的領(lǐng)先物理平臺之一.利用離子阱技術(shù)構(gòu)建大規(guī)模量子網(wǎng)絡(luò)的規(guī)?；桨竅6]包括光子鏈接[3]和離子穿梭[7]. 光子鏈接方案結(jié)合了遠(yuǎn)程量子通信和局域量子計算, 為實現(xiàn)量子網(wǎng)絡(luò)提供了一個競爭性很強(qiáng)的物理平臺. 人們已經(jīng)演示了離子光子糾纏[8]以及以光子為媒介的遠(yuǎn)距離離子糾纏[9], 這是構(gòu)建離子-光子網(wǎng)絡(luò)的基本單元. 離子穿梭方案是將囚禁離子系統(tǒng)劃分為操作區(qū)域和存儲區(qū)域, 通過離子穿梭技術(shù)在區(qū)域間移動離子.在這一架構(gòu)下, 人們已通過實驗演示了操作區(qū)域的基本結(jié)構(gòu)和操作[10]. 隨著系統(tǒng)規(guī)模的增大, 對存儲區(qū)域量子比特儲存時間的需求將相應(yīng)增加. 為了將量子比特誤碼率保持在容錯量子計算的閾值之下,延長量子比特的相干時間是十分關(guān)鍵的. 隨著系統(tǒng)規(guī)模的增大, 這兩種方案對量子比特的相干時間的要求都會相應(yīng)地提高.

二十多年前, 人們在離子系綜里利用磁場不敏感的內(nèi)態(tài)實現(xiàn)了大約10 min的相干時間[11,12]. 然而, 經(jīng)過十多年的發(fā)展后, 單量子比特系統(tǒng)的相干時間仍然遠(yuǎn)低于系綜相干時間, 只有不到1 min[13?15]. 對于囚禁離子量子比特, 特別是基于超精細(xì)能級的量子比特, 由自發(fā)輻射導(dǎo)致的固有弛豫時間T1遠(yuǎn)長于由磁場漲落和本振源的相位噪聲導(dǎo)致的退相干時間T2*. 目前在離子阱系統(tǒng)中限制相干時間的主要是環(huán)境磁場漲落和本振源相位噪聲導(dǎo)致的非相干演化, 同時我們還注意到, 加熱導(dǎo)致的測量效率降低也極大地限制了離子阱系統(tǒng)中相干時間的測量[16,17]. 具體來說, 由于運(yùn)動能級加熱, 離子的波包隨時間不斷彌散, 收集到的熒光光子數(shù)減少, 這使得區(qū)分離子內(nèi)態(tài)的測量效率低下.在硬件方面, 通過清潔離子阱表面[18]或?qū)㈦x子阱置入低溫環(huán)境中[19]可以顯著地抑制運(yùn)動能級加熱效應(yīng). 在同時囚禁不同種類離子的混合囚禁系統(tǒng)中,還可以通過協(xié)同冷卻來消除運(yùn)動能級加熱, 從而保證在整個實驗過程中探測效率不降低. 此外, 動力學(xué)解耦技術(shù)能夠消除來自環(huán)境磁場和本振源的噪聲, 從而延長量子系統(tǒng)的相干時間[20?22]. 最近有文獻(xiàn)報道, 在離子化供體的核自旋系綜中實現(xiàn)了小時量級的相干時間[23,24].

2 協(xié)同冷卻

令人奇怪的是, 二十年前關(guān)于離子系綜的實驗中[11,12]測得的相干時間遠(yuǎn)比在最近的單離子實驗中測得的相干時間來得長, 前者大約是10 min, 而后者最長不超過1 min[13?15]. 現(xiàn)在已弄清了導(dǎo)致這一現(xiàn)象的原因. 原因之一是能夠操作單量子比特的離子阱系統(tǒng)復(fù)雜度遠(yuǎn)高于之前的離子阱系綜系統(tǒng),導(dǎo)致很難在系統(tǒng)上實施磁屏蔽以降低磁場噪聲; 原因之二是離子阱系統(tǒng)的背景加熱使得同種離子的探測效率隨著時間的增加而大量降低. 之前離子阱系綜的系統(tǒng)使用了非常大型的離子阱結(jié)構(gòu), 這種結(jié)構(gòu)的電極與離子云距離較遠(yuǎn), 背景加熱較低. 但是為了能夠?qū)崿F(xiàn)速度較為理想的量子比特門, 量子信息領(lǐng)域中使用的離子阱結(jié)構(gòu)都更加微小, 導(dǎo)致背景加熱效應(yīng)更加顯著. 因此我們相信, 只要找到合適的方法, 能夠在激光冷卻停止后減輕運(yùn)動加熱效應(yīng), 單離子實驗應(yīng)該可以展示出和離子系綜可比擬的長相干時間.

對光學(xué)性質(zhì)不同的兩種離子的混合離子晶體執(zhí)行協(xié)同冷卻操作可以壓制運(yùn)動加熱效應(yīng), 同時不影響編碼在其中一種離子的內(nèi)態(tài)上的量子比特中的量子信息的相干性. 在基于囚禁離子的量子信息領(lǐng)域典型的Paul離子阱中, 通過同時囚禁兩種類型的離子, 這里選擇的是171Yb+和138Ba+, 我們可以實施協(xié)同冷卻. 其中, 單個171Yb+離子作為量子比特, 單個138Ba+作為制冷離子.171Yb+離子具有的基態(tài)電子組態(tài)具有兩個磁場不敏感的超精細(xì)能級, 這使得編碼于其上的量子比特本身具有較長的內(nèi)稟相干時間. 同時, 選用138Ba+離子是因為這兩種離子的質(zhì)量相對接近, 使得協(xié)同冷卻效率更高.

2.1 實驗裝置

協(xié)同冷卻的原理是通過對一種離子施加冷卻激光, 從而冷卻由兩種離子組成的混合離子鏈的運(yùn)動模式, 因此需要在離子阱中同時囚禁兩種類型的離子. 圖1(a)所示為一個標(biāo)準(zhǔn)的Paul離子阱同時囚禁了171Yb+離子和138Ba+離子. 圖中, 兩種離子所對應(yīng)的操作激光均覆蓋整個離子鏈, 這是因為離子鏈中的離子會在熱噪聲的作用下交換位置. 同時, 由171Yb+離子和138Ba+離子的能級圖可知, 兩種離子的光學(xué)特性差異很大, 對138Ba+離子進(jìn)行激光冷卻不會影響編碼在171Yb+離子上的量子比特的相干性.

171Yb+離子2S1/2電子組態(tài)的兩個超精細(xì)能級被定義為量子比特的兩個能級. 這兩個能級的特征頻率之差為 1 2642812118+310.8B2Hz, 其中B是以高斯為單位的磁場強(qiáng)度的數(shù)值. 在實驗中, 可以通過標(biāo)準(zhǔn)的光泵技術(shù)將量子比特制備到 |↓〉 態(tài), 并通過熒光測量技術(shù)來區(qū)分量子比特的兩個能級.

圖1(b)是用于產(chǎn)生激光和微波信號的控制系統(tǒng)的示意圖. 其中微波信號由微波振蕩器的信號和直接數(shù)字合成器的信號混合而成, 并受到可編程邏輯門陣列的控制. 該控制系統(tǒng)能夠在 1 00 ns 內(nèi)改變微波信號的相位. 所有微波源相位均以銣原子鐘為標(biāo)準(zhǔn). 為了處理波長為369 nm的激光的泄漏問題,設(shè)計了由聲光調(diào)制器、電光脈沖選擇器和機(jī)械快門組成的三段開關(guān). 該激光用于171Yb+離子的多普勒冷卻、熒光探測和光泵操作. 因為該激光頻率與量子比特的躍遷頻率接近, 稍有泄漏即會導(dǎo)致巨大的退相干效應(yīng).

2.2 冷卻效果

衡量協(xié)同冷卻效果的直接途徑是測量171Yb+離子在協(xié)同冷卻過程中達(dá)到的平衡態(tài)溫度. 為測量平衡態(tài)溫度, 需要利用雙光子拉曼躍遷來耦合內(nèi)態(tài)與運(yùn)動能級, 從而測量運(yùn)動能級布居數(shù)分布. 然而,這需要在實驗系統(tǒng)中增加額外的激光和相應(yīng)控制系統(tǒng).

圖1 實驗裝置和控制系統(tǒng) (a) 同時囚禁 1 71Yb+ 離子和 1 38Ba+ 離子的混合囚禁系統(tǒng)及相關(guān)能譜圖; (b) 微波和激光信號的控制系統(tǒng)Fig.1. Experimental setup and control system: (a) Hybrid trapping system that traps 1 71Yb+ and 1 38Ba+ simultaneously;(b) control system for generating laser and microwave signals.

另一方面, 協(xié)同冷卻的效果可以從離子鏈的穩(wěn)定性來側(cè)面反映. 在協(xié)同冷卻的過程中觀測到, 兩個離子大約每5 min交換一次位置. 除此之外,離子從未自阱中逃逸, 甚至可以在阱中穩(wěn)定存在超過一星期, 這使得超長相干時間的測量成為可能.

由于在實驗過程中離子的位置會交換, 保證離子鏈中兩個離子所處的環(huán)境完全一致是十分重要的, 比如熒光探測和光抽運(yùn)的效率、磁場和微波場的強(qiáng)度等. 在不進(jìn)行動力學(xué)解耦的情況下, 該系統(tǒng)中單量子比特相干時間大約為800 ms.

3 動力學(xué)解耦

在利用協(xié)同冷卻技術(shù)抑制了運(yùn)動能級加熱效應(yīng)后, 剩下的退相干的主要因素是磁場漲落和用于操作量子比特的微波的相位漲落, 它們將導(dǎo)致量子比特的相位隨機(jī)化. 利用動力學(xué)解耦技術(shù), 可以抑制來自環(huán)境磁場漲落的外界噪聲和本振源的相位噪聲[20,25?29].

量子比特的演化由如下隨機(jī)哈密頓量決定,

其中,σ?z是泡利自旋算符的z分量,ω0是兩個能級的頻率差,是隨機(jī)相位噪聲,βB和βLO分別代表環(huán)境磁場漲落和操作微波相位漲 落 的 影 響 .β(t) 滿 足〈β(t)〉=0 和

為延長量子系統(tǒng)的相干時間, 動力學(xué)解耦技術(shù)[30]在量子信息領(lǐng)域被廣泛應(yīng)用, 是用來消除隨機(jī)相位噪聲影響的一項標(biāo)準(zhǔn)技術(shù). 動力學(xué)解耦最簡單的實施方案是Hanh回波[31], 又被稱作自旋回波.在經(jīng)過一段時間τ的退相干后, 該技術(shù)通過一個翻轉(zhuǎn)脈沖( π 脈沖)來實現(xiàn)自旋系統(tǒng)的重對焦, 從而在時間為 2τ時得到增強(qiáng)的回波信號. 然而, Hahn 回波無法消除時間尺度比τ更小的系統(tǒng)-環(huán)境相互作用擾動. 為解決這一問題, 人們試圖使用一系列等距 π 脈沖來消除高頻擾動, 即 Carr-Purcell (CP) 序列[32]和 Carr-Purcell-Meiboom-Gill (CPMG) 序列[33]. 可以形象地把這一類型的動力學(xué)解耦序列對系統(tǒng)的作用理解為將不需要的相位漲落通過時間平均的方式消除掉. 這類方案原則上要求 π 脈沖之間的間隔τ盡可能小, 然而在實際實驗中τ的取值受到硬件和系統(tǒng)最大可承受功率的限制.

另一方面, 人們也嘗試尋找利用有限數(shù)量的脈沖獲得最優(yōu)性能的動力學(xué)解耦方案. 對純退相干的系統(tǒng)-環(huán)境相互作用, 動力學(xué)解耦序列可看做系統(tǒng)演化的濾波器, 濾波函數(shù)與環(huán)境噪聲譜的重疊決定了量子態(tài)的衰減率. 比如, 前述CPMG序列可看做帶寬為 2 π/τ的高通濾波器. 因此, 如果環(huán)境噪聲譜已知, 可以設(shè)計合適的濾波函數(shù)來得到最優(yōu)的衰減率, 這就是基于非等距脈沖序列的動力學(xué)解耦方案[34]的核心思想.

此外, 動力學(xué)解耦序列的實際性能還受到脈沖誤差的限制. 在許多情況下, 脈沖誤差主要由系統(tǒng)誤差構(gòu)成, 比如經(jīng)典控制系統(tǒng)的校準(zhǔn)誤差. 源自核磁共振系統(tǒng)的組合脈沖方法[35]可以將系統(tǒng)誤差減少到O(εn) , 其中ε是裸的系統(tǒng)誤差率. 完全補(bǔ)償?shù)慕M合脈沖可以用于任意初態(tài), 替換任意單個脈沖而不影響其他脈沖. 將組合脈沖與動力學(xué)解耦結(jié)合起來的方案是 KDD?方案[21], 其中下標(biāo)?是 Knill組合脈沖[36]的參數(shù), 決定了 π 脈沖在x-y平面上的角度.

由于動力學(xué)解耦方案的性能對噪音環(huán)境的特性敏感, 我們測量了系統(tǒng)的噪聲譜, 并據(jù)此為動力學(xué)解耦序列選擇合適的參數(shù). 具體來說, 采用CPMG方案來測量系統(tǒng)噪聲譜, 采用 K DDxy方案來延長量子比特的相干時間, 圖2所示即CPMG方案及 K DDxy方案對應(yīng)的脈沖序列.

圖2 動力學(xué)解耦脈沖序列 (a) CPMG 方案; (b) K DDxy 方案Fig.2. Pulse sequence for dyanmical decoupling: (a) CPMG protocol; (b) K DDxy protocol.

3.1 濾波函數(shù)

為了刻畫動力學(xué)解耦序列對量子態(tài)演化的影響, 需要計算動力學(xué)解耦序列對應(yīng)的濾波函數(shù). 任意動力學(xué)解耦序列可以看作在T時間內(nèi)執(zhí)行n個π脈沖操作, 第i個 π 脈沖對應(yīng)的方位角和時刻分別為?i和τi, 并約定

其中與動力學(xué)解耦序列和系統(tǒng)在噪聲環(huán)境中的演化對應(yīng)的演化算符可寫為

其中單量子比特旋轉(zhuǎn)

考慮到γi=π , 上述演化算符可化簡為

其中Θ(·) 是Heaviside階梯函數(shù).

在拉姆齊條紋實驗中, 執(zhí)行動力學(xué)解耦序列之后, 測量的信號是自旋算符z分量的期待值, 即

推導(dǎo)(5)式時, 使用了對高斯型隨機(jī)變量X成立的等式定 義時域濾波函數(shù)的傅里葉變換

利用卷積定理, 信號的衰減函數(shù)χ(T) 可寫成如下形式:

3.2 環(huán)境噪聲

如前所述, 系統(tǒng)的隨機(jī)相位噪聲β(t) 包含來自隨機(jī)磁場漲落和操作微波相位漲落兩方面的影響,分別用βB(t) 和βLO(t) 表示. 操作微波相位漲落βLO與本地振蕩器的品質(zhì)有關(guān), 在使用高品質(zhì)本地振蕩器、銣原子鐘和直接數(shù)字合成器等設(shè)備后, 由頻譜分析儀測得的典型相位漲落如圖3所示.

圖3 頻譜分析儀測得的典型相位漲落噪聲譜Fig.3. Typical phase noise measured by spectrum analyzer.

其中ωHF是超精細(xì)能級劈裂,Bi是磁場在i方向的平均值,bi(t) 是i方向的磁場漲落 假設(shè)降低磁場平均值不影響漲落的大小, 則可以通過減小磁場的方式來壓制βB(t) , 從而得到更長的相干時間. 由通量計測得因此使磁場沿x方向. 在研究系統(tǒng)噪聲譜時, 采用在實現(xiàn)長相干時間量子存儲時, 采用這是能夠?qū)﹄x子的內(nèi)態(tài)進(jìn)行有效熒光測量的最小磁場.

系統(tǒng)中最強(qiáng)的噪聲組分是來自于電源線的50 Hz振蕩及其諧波, 因此頻域噪聲譜可以用分立模型來模擬[37,38],其中δ(·)代表狄拉克δ函數(shù). 將β?(ω) 代入 (5) 式, 得到最終 信號為

其中 J0(·) 是0階Bessel函數(shù). 實驗上測得末態(tài)平均值的平均脈沖間距變化的數(shù)據(jù)后, 可通過數(shù)值擬合得到分立噪聲譜強(qiáng)度βk.

接著, 通過連續(xù)噪聲模型來進(jìn)一步研究系統(tǒng)噪聲譜. 對任意噪聲譜拉姆齊條紋的對比度可表示為, 其中

4 基于微波的量子門操作

量子比特躍遷可通過微波來驅(qū)動, 微波信號是通過混合來自安捷倫微波振蕩器的12.442812 GHz的信號和來自直接數(shù)字合成器的 2 00 MHz 的信號產(chǎn)生的. 經(jīng)過一次放大, 微波信號通過放置在視窗附近的微波喇叭施加于系統(tǒng). 由于動力學(xué)解耦序列通常包含數(shù)千乃至數(shù)萬個單比特邏輯門, 邏輯門誤差的累積會導(dǎo)致量子比特退相干. 因此, 微波驅(qū)動的單比特量子邏輯門的保真度是實現(xiàn)長相干時間量子存儲的最重要的因素之一.

4.1 隨機(jī)化基準(zhǔn)測試

隨機(jī)化基準(zhǔn)測試[39]是用來測量極低的單比特量子邏輯門的保真度的一種標(biāo)準(zhǔn)方法. 傳統(tǒng)上, 量子操作的誤差行為可以通過量子過程層析來標(biāo)定.然而, 雖然量子過程層析可以提供關(guān)于量子操作的許多信息, 但卻無法保證該量子操作嵌入一長串其他量子操作中的表現(xiàn). 另一方面, 規(guī)?；孔佑嬎阋髥蝹€量子操作的誤差率在 1 0-4量級或更低.通過傳統(tǒng)的量子過程層析在實驗上測定這樣低的誤差率是很困難的. 此外, 量子過程層析無法排除初態(tài)制備和測量過程引入的誤差. 因此, 人們提出了一種可以克服以上問題的對量子操作的誤差行為進(jìn)行實驗研究的新方法, 即隨機(jī)化基準(zhǔn)方法.

圖4 (a) 執(zhí)行包含 N 個脈沖的動力學(xué)解耦序列, 不同的總演化時間 T 對應(yīng)的條紋對比度; (b) 通過分析圖(a)中的數(shù)據(jù)得到的環(huán)境噪聲譜Fig.4. (a) Ramsey contrasts depending on the total evolution time T for various numbers of pulses N in the dynamical decoupling sequence; (b) the noise spectra analized from the measured data in Fig.(a).

隨機(jī)化方法最初是在量子噪聲領(lǐng)域提出的, 其具體過程是對系統(tǒng)實施隨機(jī)幺正操作及其逆操作假如噪聲模型可表示為施加在之間的獨(dú)立的量子操作, 保真度的下降程度即代表了系統(tǒng)的噪聲強(qiáng)度. 隨機(jī)化基準(zhǔn)方法是隨機(jī)化方法的簡化版本, 表現(xiàn)在前者將隨機(jī)幺正操作限制為Clifford門操作.

隨機(jī)化基準(zhǔn)方法對標(biāo)準(zhǔn)初態(tài)執(zhí)行不同長度的隨機(jī)的量子邏輯門序列, 每一序列均以隨機(jī)化測量結(jié)束, 該測量決定系統(tǒng)是否達(dá)到正確的末態(tài). 根據(jù)最終測量的錯誤率隨序列長度的增加可以確定序列中每一個門操作平均的計算相干的誤差率. 序列中的隨機(jī)門操作出自Clifford群, 該群是用e-iσπ/4的 π /2 旋轉(zhuǎn)生成的, 其中σ是作用于不同量子比特的泡利算符的直積.

對于單量子比特系統(tǒng), 隨機(jī)化基準(zhǔn)方法包含大量實驗, 每一實驗的流程均包含初態(tài)制備、執(zhí)行不同的量子邏輯門序列和測量, 其中量子邏輯門序列包含泡利隨機(jī)化和門計算兩部分. 泡利隨機(jī)化執(zhí)行具有形式的幺正操作, 其中隨機(jī)變量是恒等算符. 門計算執(zhí)行形式為的幺正操作,u∈{x,y}. 除了最后一個 π /2 脈沖, 上述符號、隨機(jī)變量b和μ都由均勻隨機(jī)采樣得到. 選擇系統(tǒng)初態(tài)為|0〉, 是由于序列中的操作均出自Clifford 群, 理想情況下系統(tǒng)的狀態(tài)總是泡利算符的本征態(tài). 最后一個 π /2 脈沖將系統(tǒng)狀態(tài)旋轉(zhuǎn)為的本征態(tài), 最終測量是關(guān)于的馮諾依曼測量. 在理想情況下, 最終測量的結(jié)果是已知的. 在誤差環(huán)境下, 由于隨機(jī)化, 序列的誤差與組成序列的單個或數(shù)個操作無關(guān). 序列的長度l是序列中 π /2 脈沖的個數(shù). 可以測量不同長度l的序列對應(yīng)的誤差率pl, 由pl與l的關(guān)系可以估算門脈沖的平均誤差率,

其中d是隨機(jī)化計算門(一個 π /2 脈沖和一個泡利脈沖的隨機(jī)組合)的平均誤差率,dif為包含初態(tài)制備和測量過程的誤差率.

圖5展示了不同序列的保真度, 對于序列長度l的每一個不同取值都存在32個數(shù)據(jù)點(diǎn), 由此可以估算單量子比特邏輯門保真度為99.994%±0.002%.

5 測量單量子比特相干時間

圖5 序列保真度隨序列長度的變化Fig.5. Sequence fidelity pl as a function of the sequence length l.

結(jié)合協(xié)同冷卻和動力學(xué)解耦, 最近我們實現(xiàn)了具有超長相干時間的單比特量子存儲, 并測得超過10 min的相干時間[40]. 我們在標(biāo)準(zhǔn)的Paul阱中同時囚禁171Yb+ 和138Ba+ 兩種離子, 其中171Yb+離子2S1/2電子組態(tài)中的兩個超精細(xì)能級被定義為量 子 比 特mF=0〉. 在實驗過程中, 持續(xù)對138Ba+離子進(jìn)行激光冷卻. 在測量了系統(tǒng)的環(huán)境噪聲譜后, 通過合適的動力學(xué)解耦方案, 即脈沖平均間隔的 K DDxy方案, 將環(huán)境噪聲漲落對編碼于138Yb+離子之上的量子比特的影響降到最低.

在測量單量子比特相干時間的過程中, 將磁場設(shè)定為 3.5 G 來最大限度地減小磁場漲落的影響.選取六個不同的初態(tài), 分別是對于每一個初態(tài)執(zhí)行動力學(xué)解耦序列, 并測量拉姆齊條紋的對比度. 以為例, 實驗的脈沖序列如圖6所示. 其中初始化包括多普勒冷卻、激光抽運(yùn)和 π /2 脈沖. 隨后, 執(zhí)行一定時間的KDDxy序列, 再執(zhí)行 π /2 脈沖. 在理想情況下, 最后一個π/2脈沖將系統(tǒng)狀態(tài)旋轉(zhuǎn)到|↑〉. 最后, 測量結(jié)果為0或1. 為得到代表拉姆齊條紋對比度的每個實驗被重復(fù)30—50次.

圖6 測量單量子比特相干時間的脈沖序列Fig.6. Pulse sequences for measuring the single-qubit coherent time.

圖7 單量子比特六個不同初態(tài)的相干時間, 其中 |↑〉 和 |↓〉 ,對應(yīng)的相干時間是(4740 ± 1760) s; 其他四個初態(tài)對應(yīng)的相干時間為(667 ± 17) s; 圖中的誤差線代表標(biāo)準(zhǔn)差Fig.7. Single-qubit coherece time for six different initial states. For |↑〉 and |↓〉 , the coherence time is (4740±1760) s. For the other four initial states, the coherence time is (667 ± 17) s. The error bars are the standard deviation.

對每一個初態(tài), 均測量演化時間T后的拉姆齊條紋對比度, 其中600 s}.圖7展示了實驗上測得的六個不同初態(tài)在演化時間T后的拉姆齊條紋對比度. 其中,|↑〉和|↓〉的條紋對比度在測量時間內(nèi)有略微下降, 導(dǎo)致該下降的主要原因是量子邏輯門操作誤差的累積, 對應(yīng)的相干時間是 ( 4740±1760)s . 其余四個初態(tài)對應(yīng)的相干時間是 ( 667±17)s , 導(dǎo)致退相干的主要原因是動力學(xué)解耦的帶通濾波函數(shù)中心頻率附近的殘留噪聲譜強(qiáng)度. 上述相干時間是通過對拉姆齊條紋對比度的數(shù)據(jù)進(jìn)行指數(shù)擬合得到的.

6 總結(jié)與展望

在同時囚禁171Yb+ 離子和138Ba+ 離子的混合囚禁離子系統(tǒng)中, 結(jié)合協(xié)同冷卻和動力學(xué)解耦, 本文實現(xiàn)了相干時間超過10 min的單離子量子比特. 雖然超過10 min的單量子比特相干時間已經(jīng)刷新了此前該領(lǐng)域的世界紀(jì)錄, 但我們注意到該相干時間仍然可以進(jìn)一步延長, 即實現(xiàn)更長的單比特相干時間并沒有基礎(chǔ)原理上的限制. 具體來說, 可以通過以下幾種方式來進(jìn)一步延長量子比特相干時間: 1) 安裝磁場屏蔽; 2) 使用磁場不敏感的量子比特; 3) 使用相干性更好、相位更穩(wěn)定的微波源.

如果囚禁更多的用于存儲量子信息的量子比特, 輔以協(xié)同冷卻技術(shù), 該系統(tǒng)有望實現(xiàn)量子密碼學(xué)領(lǐng)域的理論方案, 包括量子錢幣[5]等. 本文所實現(xiàn)的具有超長相干時間的量子存儲將刺激混合量子計算系統(tǒng)[41,42]的發(fā)展, 比如將具有有限相干時間的量子系統(tǒng)中的量子信息映射并儲存到囚禁離子量子比特中.