不同光學(xué)網(wǎng)絡(luò)結(jié)構(gòu)玻色采樣發(fā)生隨機光子損失的模擬研究*

季陽 陳美玲 黃汛 吳永政 蘭冰

(中國電子科技集團公司第三十二研究所,上海 201808)

玻色采樣機是最有可能真正意義上實現(xiàn)量子優(yōu)勢的專用量子計算機之一,其在量子化學(xué)等領(lǐng)域也有著很好的應(yīng)用前景.然而,光子損失相關(guān)的噪聲會引起玻色采樣樣本的誤差.為了研究光子損失對玻色采樣的影響,基于等效分束器原理,對玻色采樣展開經(jīng)典的模擬研究.對于對應(yīng)任意幺正矩陣的兩種光學(xué)網(wǎng)絡(luò),當(dāng)在每一個光學(xué)單元中有一定概率發(fā)生光子損失時,發(fā)現(xiàn)具有Clements 結(jié)構(gòu)光學(xué)網(wǎng)絡(luò)的玻色采樣相比于Reck 結(jié)構(gòu)的樣本誤差更小.進一步地,當(dāng)光子損失的概率符合正態(tài)分布時,發(fā)現(xiàn)若光子損失概率平均值不變,標(biāo)準(zhǔn)差越大,玻色采樣的樣本誤差越大.最后,考慮輸出光子數(shù)的實驗結(jié)果表明,隨著光子損失概率的提高,無光子損失的輸出組合比例急劇下降,說明光子損失會明顯影響玻色采樣的量子優(yōu)勢.隨機光子損失的玻色采樣模擬研究有助于玻色采樣實驗的開展,為含噪聲量子計算研究提供思路.

1 引言

量子計算機在解決某些復(fù)雜問題,如大數(shù)分解問題[1]、搜索問題[2]時相比于經(jīng)典計算機具有指數(shù)級加速的潛力,因此吸引了越來越多的關(guān)注.但受限于目前的技術(shù)手段,具有大規(guī)模以及高保真度的通用量子計算機仍然難以在短時間內(nèi)實現(xiàn).基于此,研究專用量子計算機具有很強的現(xiàn)實意義.這主要是因為,專用量子計算機可以率先實現(xiàn)真正意義上的量子優(yōu)勢,即量子計算機在解決某些高復(fù)雜度問題上相比傳統(tǒng)計算機具有明顯的加速效果.

玻色采樣機作為一種專用量子計算機,由于其具有結(jié)構(gòu)簡單、類型多樣[3-6]、可以通過多種手段實現(xiàn)[3,7]、不需要量子邏輯門[8]、經(jīng)典驗證方法成熟[9-11]等優(yōu)勢,是近年來量子計算領(lǐng)域的研究熱點之一.玻色采樣概念的提出是為了在計算復(fù)雜度問題上推翻廣義“丘奇-圖靈”論題[8],這一目標(biāo)隨著2020年的76 光子100 模式高斯玻色采樣實驗的開展已基本實現(xiàn)[4].除此之外,在超導(dǎo)量子線路上進行的玻色采樣也成功模擬了分子振動譜[12],這說明玻色采樣在量子化學(xué)領(lǐng)域也具有很好的應(yīng)用前景.

但是,玻色采樣的超大規(guī)模實現(xiàn)始終面臨著一些困難,主要有光子非全同和光子損失.當(dāng)光子非全同和光子損失所引起的噪聲足夠明顯時,玻色采樣的樣本可以由經(jīng)典算法在多項式時間內(nèi)模擬得到,此時玻色采樣機失去了量子優(yōu)勢[13].因此,對光子非全同和光子損失的研究具有重要意義.光子非全同指的是光子在某一個內(nèi)在量子態(tài)上存在差異,例如光子的頻率或者位置.人們通過一系列理論研究,得到了光子非全同的玻色采樣輸出組合概率及對應(yīng)的高效經(jīng)典估算方法[14,15].在實驗方面,相比于傳統(tǒng)的自發(fā)參量下轉(zhuǎn)換,后來發(fā)展的光學(xué)微腔-半導(dǎo)體量子點體系制備出來的光子全同度更高,更加適合大規(guī)模實現(xiàn)典型玻色采樣[16].

同時,研究者們也對光子損失玻色采樣的理論和實驗開展了深入的研究.光子損失是指光子在產(chǎn)生、傳輸、探測過程中發(fā)生丟失的情況,其與光子源、光學(xué)網(wǎng)絡(luò)、探測器的質(zhì)量密切相關(guān).與光子非全同一樣,光子損失會明顯影響玻色采樣的量子優(yōu)勢.對于光子損失對玻色采樣量子優(yōu)勢影響的分析,主要基于兩種理論模型.一種模型是對輸出端損失光子的個數(shù)進行統(tǒng)計分析.例如,對于輸入光子數(shù)為n,輸出光子數(shù)為(n?l):若l為一個常數(shù),則玻色采樣仍然具有量子優(yōu)勢;若則玻色采樣失去量子優(yōu)勢[17].更一般地,隨著輸入光子數(shù)的增加: 若損失光子數(shù)與輸入光子數(shù)之比趨向于0,則玻色采樣仍然具有量子優(yōu)勢;若損失光子數(shù)與輸入光子數(shù)之比趨向于1,則玻色采樣失去量子優(yōu)勢[13].第二種模型則是考慮光子損失有可能發(fā)生在每一個光學(xué)元件中,其更加符合實際實驗?zāi)Ｐ蚚18].這種模型將光子損失等效為光子通過一個虛構(gòu)的分束器進入到環(huán)境光路中.基于此,對光子損失的玻色采樣模擬實際上等效于對更大規(guī)模光學(xué)網(wǎng)絡(luò)的玻色采樣模擬.這種模型揭示了光子損失玻色采樣的量子優(yōu)勢與光學(xué)網(wǎng)絡(luò)的深度密切相關(guān)[17].

本工作基于光子損失情況下玻色采樣的第二種理論模型展開,對不同結(jié)構(gòu)的光學(xué)網(wǎng)絡(luò)模型進行了定量的分析,給出了直觀的比較結(jié)果.同時本文的模擬結(jié)果由于考慮了光子損失概率的隨機性,更加符合真實物理場景.

2 模 型

2.1 理論模型

目前比較常用的兩種光學(xué)線性網(wǎng)絡(luò),即Reck結(jié)構(gòu)和Clements 結(jié)構(gòu)的光學(xué)網(wǎng)絡(luò),如圖1 所示,圖中每一個交點都表示一個光學(xué)單元.

圖1 光學(xué)線性網(wǎng)絡(luò)示意圖: (a) Reck 結(jié)構(gòu);(b) Clements 結(jié)構(gòu)Fig.1.Sketches of optical linear networks: (a) the Reck structure;(b) the Clements structure.

其中,Reck 結(jié)構(gòu)為三角形非對稱結(jié)構(gòu),而Clements 結(jié)構(gòu)為四邊形對稱結(jié)構(gòu)[19].直接對任意m維幺正矩陣進行分解,使其可以表示為一系列矩陣相乘的形式.在該系列矩陣中,第一個矩陣是對角矩陣,對應(yīng)光學(xué)網(wǎng)絡(luò)末端的m個相移器;其余矩陣均對應(yīng)一個光學(xué)單元,每個光學(xué)單元由一個相移器和一個分束器組成.每個光學(xué)單元對應(yīng)的幺正矩陣M可以表示為分束器矩陣MBS和相移器矩陣MPS的乘積[19]:其中分束器的反射率是cosω,相移器改變的相位是φ,分束器矩陣含參元素所在行(列)數(shù)為分束器所在的兩個通道數(shù),相移器矩陣含參元素所在行(列)數(shù)為相移器所在的通道數(shù).

2.2 實驗

對于任意的幺正矩陣U,將其分解成一系列矩陣相乘的形式,該系列矩陣對應(yīng)的光學(xué)單元的排列符合Reck 結(jié)構(gòu),其分解方法如算法1 所示[20].

對于Clements 結(jié)構(gòu),U的分解方法如算法2所示[19].Clements 結(jié)構(gòu)的光學(xué)網(wǎng)絡(luò)對應(yīng)的矩陣構(gòu)造方法要比Reck 結(jié)構(gòu)復(fù)雜得多.這主要是因為,U需要交替右乘和左乘一系列矩陣來實現(xiàn)對角化,即MM···MUM?1M?1···M?1=D,其 中D為對角矩陣,M表示一系列參數(shù)(通道數(shù),反射率以及相位)不盡相同的光學(xué)單元矩陣,則U可以表示為形如M?1M?1···M?1DMM···M的形式.在此基礎(chǔ)上,可以初始化一個對角矩陣D'以及一個光學(xué)單元矩陣M',使得D′M′=M?1D.經(jīng)過計算可以發(fā)現(xiàn),D'和M'的參數(shù)是唯一的,則U可以表示為形如M-1M-1···M-1D'M'MM···M的形式,此時該形式中對角矩陣左乘的矩陣減少1 個,而右乘的矩陣增加1 個.通過重復(fù)構(gòu)造D'和M',最終U可以表示為形如D'M'M'···M'MM···M的形式,實現(xiàn)基于Clements 結(jié)構(gòu)對U的分解.

在得到對應(yīng)任意幺正矩陣的光學(xué)網(wǎng)絡(luò)后,為了研究在光學(xué)網(wǎng)絡(luò)中發(fā)生的光子損失對玻色采樣的影響,基于等效分束器原理對光學(xué)網(wǎng)絡(luò)進行重構(gòu).2 個虛構(gòu)的分束器和2 條環(huán)境光路被加入每個光學(xué)單元中.虛構(gòu)的分束器參數(shù)ω和單條實際光路在單個光學(xué)單元中的光子損失概率Ploss滿足cosω=1 — (1 —Ploss)1/2關(guān)系[18].當(dāng)光子在該光學(xué)單元中發(fā)生損失,則其等效于光子通過虛構(gòu)的分束器進入環(huán)境光路中.由于兩種光學(xué)網(wǎng)絡(luò)都具有m(m— 1)/2個光學(xué)單元,在重構(gòu)的光學(xué)網(wǎng)絡(luò)中,一共附加了m(m— 1)個虛構(gòu)分束器和相同數(shù)目的環(huán)境光路,則重構(gòu)的光學(xué)網(wǎng)絡(luò)的模式數(shù)變?yōu)閙2,因此原始m維矩陣M需要變換為m2維矩陣,即M⊕Em(m—1),其中Em(m—1)為m(m— 1)維單位矩陣.相應(yīng)地,虛構(gòu)的分束器矩陣也是m2維矩陣.在得到重構(gòu)后的光學(xué)網(wǎng)絡(luò)后,則可以計算考慮光子損失的玻色采樣輸出組合概率,其方法如算法3 所示,在計算輸出組合概率的過程中,需要計算矩陣的積和式,可以采用Ryser 或者Glynn 方法,并利用格雷碼以簡化計算過程[21].計算了4 光子16 模式的玻色采樣輸出組合概率,其中16 維輸入組合S=[1,1,1,1,0,0,···,0].一共隨機產(chǎn)生了100 個16 維隨機幺正矩陣,并對最終計算結(jié)果取平均值.

3 實驗結(jié)果及討論

基于玻色子生日悖論,在只考慮輸出組合中每個通道只有0 或1 個光子的情況時,可以通過總變差距離來評估玻色采樣的輸出結(jié)果與理想情況的偏差程度[21].總變差距離可以表示為Pideal|,其中T4表示沒有發(fā)生光子損失的任意16維輸出組合,Pideal表示理想情況下歸一化的輸出組合概率.為了便于與理想情況對比,Pout也需要進行歸一化處理,使得對于所有的T4,Pout之和為1.如圖2(a)所示,我們發(fā)現(xiàn)隨著光子損失概率的提高,總變差距離以準(zhǔn)線性的方式逐步變大,這表明光子損失會對玻色采樣的精度產(chǎn)生顯著影響.我們可以直觀地發(fā)現(xiàn),對于固定的Ploss,Clements結(jié)構(gòu)比Reck 結(jié)構(gòu)對應(yīng)的總變差距離更小,表明Clements 結(jié)構(gòu)在光子損失情況下穩(wěn)定性更好.這主要是因為Clements 結(jié)構(gòu)受到平衡與非平衡光子損失的影響更小[19].在光子損失情況下,一方面,對于平衡光子損失,光學(xué)網(wǎng)絡(luò)深度越大,光子損失程度越高.對于較大規(guī)模光學(xué)網(wǎng)絡(luò),Clements 結(jié)構(gòu)的深度比Reck 結(jié)構(gòu)小很多,因此受到平衡光子損失的影響更小;另一方面,對于非平衡光子損失,光子損失程度與通道中分束器的光子損失概率密切相關(guān).在Reck 結(jié)構(gòu)中,由于每個通道經(jīng)過的分束器數(shù)目不同,導(dǎo)致通道中光子損失程度不同,這極大地影響了多光子相干干涉結(jié)果.而在Clements結(jié)構(gòu)中,每個通道經(jīng)過的分束器數(shù)目大致相同,通道中光子損失程度大致相當(dāng),因此在Clements 結(jié)構(gòu)中非平衡光子損失對多光子相干干涉結(jié)果的影響相比于Reck 結(jié)構(gòu)較小.

圖2 (a) 考慮光子損失的玻色采樣輸出組合概率的總變差距離隨光子損失概率Ploss的變化關(guān)系;(b) Ploss為0.1 時,兩種光學(xué)網(wǎng)絡(luò)結(jié)構(gòu)隨著模式數(shù)增加對應(yīng)玻色采樣結(jié)果總變差距離的變化關(guān)系.圖中空心圖形代表實驗結(jié)果平均值,誤差棒代表標(biāo)準(zhǔn)差Fig.2.Plots between the total variation distance of boson sampling output combination probabilities considering photon losses and(a) the photon loss probability Ploss or (b) mode numbers when Ploss is 0.1.Open shapes represent mean values and error bars represent standard deviations.

為了驗證上述結(jié)論,我們對不同規(guī)模的玻色采樣進行了仿真,如圖2(b)所示,所采用玻色采樣規(guī)模分別為1 光子3 模式、2 光子6 模式、3 光子9 模式、3 光子12 模式和4 光子16 模式.在模式數(shù)為3 時,兩種結(jié)構(gòu)具有同樣的形式,因此總變差距離相同.隨著模式數(shù)的增加,光學(xué)網(wǎng)絡(luò)深度增加,通道經(jīng)過的最大分束器數(shù)目增加,導(dǎo)致兩種結(jié)構(gòu)對應(yīng)的總變差距離也逐漸增加,表明光子損失(包括平衡與非平衡光子損失)對玻色采樣結(jié)果的影響逐漸加大.而Reck 結(jié)構(gòu)對應(yīng)的總變差距離增幅明顯高于Clements 結(jié)構(gòu),表明Clements 結(jié)構(gòu)受到平衡與非平衡光子損失的影響相對于Reck 結(jié)構(gòu)更小.

進一步假設(shè)每個光學(xué)單元對應(yīng)的Ploss不是一個定值,而是符合一定的隨機分布,即正態(tài)分布.可認(rèn)為這樣的假設(shè)比較能夠反映出真實的物理實驗情況.分兩種情況進行了模擬: 第1 種情況是,固定Ploss的標(biāo)準(zhǔn)差不變,改變其平均值.對于這種情況,發(fā)現(xiàn)總變差距離的變化情況與圖2(a)類似;第2 種情況是,固定Ploss的平均值不變,改變其標(biāo)準(zhǔn)差.研究發(fā)現(xiàn),對于兩種光學(xué)網(wǎng)絡(luò)結(jié)構(gòu),隨著Ploss標(biāo)準(zhǔn)差的增大,輸出組合概率的總變差距離呈現(xiàn)變大的趨勢,如圖3 所示,表明光子損失概率的離散程度也會明顯影響玻色采樣的精度.以深度為3的光學(xué)網(wǎng)絡(luò)為例,對這種情況進行簡單的分析.此時光學(xué)網(wǎng)絡(luò)的每一層只含有一個光學(xué)單元.在Ploss為一個定值P時,每一層的光子透過率為T=1 —P,因此總光子透過率是T3.若Ploss符合平均值為P的正態(tài)分布,則每一層的光子透過率也符合平均值為T的正態(tài)分布,總光子透過率的期望值仍為T3.為了滿足這個條件,隨著Ploss標(biāo)準(zhǔn)差的增加,總光子透過率小于T3的概率需要逐漸增加.受較低的總光子透過率的影響,輸出組合概率的總變差距離相應(yīng)增加,表明光子損失對玻色采樣精度的影響變得更為顯著.

圖3 考慮光子損失的玻色采樣輸出組合概率的總變差距離平均值隨光子損失概率的標(biāo)準(zhǔn)差σPloss的變化關(guān)系Fig.3.Plots between the mean values of total variation distance of boson sampling output combination probabilities considering photon losses and the standard deviation σPloss of photon loss probability.

最后,統(tǒng)計了固定光子個數(shù)的輸出組合的概率之和隨著光子損失概率變化的情況,如圖4 所示.當(dāng)Ploss=0.01 時,Clements 結(jié)構(gòu)相比于Reck 結(jié)構(gòu),4 光子輸出的概率明顯較高.由于輸出光子數(shù)越高,其對應(yīng)的經(jīng)典計算的時間復(fù)雜度也越高,我們的結(jié)果表明在光子損失情況下,Clements 結(jié)構(gòu)的玻色采樣對于經(jīng)典計算來說具有更高的時間復(fù)雜度.而隨著Ploss的提高,4 光子輸出概率急劇減小,3 光子輸出占據(jù)優(yōu)勢.進一步提高Ploss,則2 光子輸出變成主導(dǎo).這表明,發(fā)生在光學(xué)網(wǎng)絡(luò)中的光子損失能夠顯著地影響玻色采樣的量子優(yōu)勢.當(dāng)?shù)凸庾訑?shù)的輸出組合占據(jù)主導(dǎo)時,玻色采樣變得易于被經(jīng)典算法模擬,如MCMC 算法[22]、Clifford-Clifford 算法[23]等.

圖4 玻色采樣固定光子個數(shù)的輸出組合概率之和隨著Ploss的變化關(guān)系.(a) Reck結(jié)構(gòu);(b) Clements 結(jié)構(gòu)光 學(xué)網(wǎng)絡(luò).圖中空心圖形代表實驗結(jié)果平均值,誤差棒代表標(biāo)準(zhǔn)差Fig.4.Plots between the sum of boson sampling output combination probabilities with fixed photon numbers and Ploss: (a) the Reck structure;(b) the Clements structure.Open shapes represent mean values and error bars represent standard deviations.

基于對光子損失情況下玻色采樣的經(jīng)典仿真結(jié)果,發(fā)現(xiàn)光子損失能夠顯著影響玻色采樣結(jié)果的精度,進而影響玻色采樣量子優(yōu)勢的可實現(xiàn)性.如何抑制光子損失對玻色采樣機等光量子器件的影響一直是研究熱點之一.在理論方面,系統(tǒng)研究了輸出光子個數(shù)[24]、光學(xué)網(wǎng)絡(luò)深度[17]等對玻色采樣對應(yīng)經(jīng)典計算時間復(fù)雜度的影響.在工程技術(shù)方面,減反射涂層[16]、光路拓展[25]、三維光學(xué)網(wǎng)絡(luò)[26]的應(yīng)用提供了抑制光子損失影響的可能性.值得一提的是,在容錯全光量子中繼器方案中,利用量子奇偶碼有效抑制了光子損失的影響,這給未來光量子計算器件抵抗光子損失帶來了啟示[27].

4 結(jié)論

系統(tǒng)研究了兩種結(jié)構(gòu)的光學(xué)網(wǎng)絡(luò)中光子損失對玻色采樣的影響.可以直觀地發(fā)現(xiàn),當(dāng)存在光子損失時,Clements 結(jié)構(gòu)相比于Reck 結(jié)構(gòu)具有更好的穩(wěn)定性,其輸出高光子個數(shù)的概率更高,對應(yīng)經(jīng)典計算的時間復(fù)雜度更高.隨著光子損失概率的提高,玻色采樣的精度明顯降低,且輸出高光子個數(shù)的概率急劇降低,表明光子損失對玻色采樣的經(jīng)典計算時間復(fù)雜度具有十分顯著的影響.研究還發(fā)現(xiàn),若光學(xué)單元中光子損失概率不是一個定值,而是服從正態(tài)分布,則光子損失概率的離散程度也能夠影響玻色采樣的精度.結(jié)果表明,發(fā)生在光學(xué)網(wǎng)絡(luò)中的光子損失能夠明顯影響玻色采樣的精度和對應(yīng)經(jīng)典計算的時間復(fù)雜度,這對真正意義上具有量子優(yōu)勢的玻色采樣的物理實現(xiàn)提出了又一個挑戰(zhàn).