Search

Article

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Simulation of random photon loss in boson sampling of different optical networks

Ji Yang Chen Mei-Ling Huang Xun Wu Yong-Zheng Lan Bing

Ji Yang, Chen Mei-Ling, Huang Xun, Wu Yong-Zheng, Lan Bing. Simulation of random photon loss in boson sampling of different optical networks. Acta Phys. Sin., 2022, 71(19): 190301. doi: 10.7498/aps.71.20220331
Citation: Ji Yang, Chen Mei-Ling, Huang Xun, Wu Yong-Zheng, Lan Bing. Simulation of random photon loss in boson sampling of different optical networks. Acta Phys. Sin., 2022, 71(19): 190301. doi: 10.7498/aps.71.20220331

Simulation of random photon loss in boson sampling of different optical networks

Ji Yang, Chen Mei-Ling, Huang Xun, Wu Yong-Zheng, Lan Bing
Article Text (iFLYTEK Translation)
PDF
HTML
Get Citation
  • Boson sampling is a candidate for quantum protocols to truly realize the quantum computation advantage and to be used in advanced fields where complex computations are needed, such as quantum chemistry. However, this proposal is hard to achieve due to the existence of noise sources such as photon losses. In order to quantificationally analyze the influences of photon losses in optical networks, boson sampling is classically simulated based on the equivalent beam splitter mechanism, where the photon loss happening in optical units is equivalent to the photon transmission into the environmental paths through a virtual beam splitter. In our simulation, networks corresponding to random unitary matrices are made up, considering both the Reck structure and the Clements structure. The photon loss probability in an optical unit is well controlled by adjusting the parameters of the virtual beam splitter. Therefore, to simulate boson sampling with photon losses in optical networks is actually to simulate ideal boson sampling with more modes. It is found that when the photon loss probability is constant, boson sampling with Clements structures distinctly performs much better than that with Reck structures. Furthermore, the photon loss probability is also set to follow the normal distribution, which is thought to be closer to the situation in reality. It is found that when the mean value of photon loss probability is constant, for both network structures, errors of outputs become more obvious with the increase of standard deviation. It can be inferred that the increase of error rate can be explained by the network depth and the conclusion is suitable for larger-scale boson sampling. Finally, the number of output photons is taken into consideration, which is directly related to the classical computation complexity. It is found that with the photon loss probability, the ratio of output combinations without photon losses decreases sharply, implying that photon losses can obviously affect the quantum computation advantage of boson sampling. Our results indicate that photon losses can result in serious errors for boson sampling, even with a stable network structure such as that of Clements. This work is helpful for boson sampling experiments in reality and it is desired to develop a better protocol, for example, a well-designed network or excellent optical units, to well suppress photon losses.
      PACS:
      03.67.Ac(Quantum algorithms, protocols, and simulations)
      42.50.-p(Quantum optics)
      Corresponding author: Chen Mei-Ling, cml13262283846@163.com
    • Funds: Project supported by the sub-project of Science and Technology Major Project of the Ministry of Science and Technology of Shanghai (Grant No. 2019SHZDZX01-ZX03), the project of “Quantum Computing and Simulation Technology Research” of the development department of China Electronics Technology Group Corporation (Grant No. GQ201173-00), and the project of “Practical Application Research of Quantum Computing” of the 32nd Research Institute of China Electronics Technology Group Corporation (Grant No. GY200906-00).

    量子计算机在解决某些复杂问题, 如大数分解问题[1]、搜索问题[2]时相比于经典计算机具有指数级加速的潜力, 因此吸引了越来越多的关注. 但受限于目前的技术手段, 具有大规模以及高保真度的通用量子计算机仍然难以在短时间内实现. 基于此, 研究专用量子计算机具有很强的现实意义. 这主要是因为, 专用量子计算机可以率先实现真正意义上的量子优势, 即量子计算机在解决某些高复杂度问题上相比传统计算机具有明显的加速效果.

    玻色采样机作为一种专用量子计算机, 由于其具有结构简单、类型多样[3-6]、可以通过多种手段实现[3,7]、不需要量子逻辑门[8]、经典验证方法成熟[9-11]等优势, 是近年来量子计算领域的研究热点之一. 玻色采样概念的提出是为了在计算复杂度问题上推翻广义“丘奇-图灵”论题[8], 这一目标随着2020年的76光子100模式高斯玻色采样实验的开展已基本实现[4]. 除此之外, 在超导量子线路上进行的玻色采样也成功模拟了分子振动谱[12], 这说明玻色采样在量子化学领域也具有很好的应用前景.

    但是, 玻色采样的超大规模实现始终面临着一些困难, 主要有光子非全同和光子损失. 当光子非全同和光子损失所引起的噪声足够明显时, 玻色采样的样本可以由经典算法在多项式时间内模拟得到, 此时玻色采样机失去了量子优势[13]. 因此, 对光子非全同和光子损失的研究具有重要意义. 光子非全同指的是光子在某一个内在量子态上存在差异, 例如光子的频率或者位置. 人们通过一系列理论研究, 得到了光子非全同的玻色采样输出组合概率及对应的高效经典估算方法[14,15]. 在实验方面, 相比于传统的自发参量下转换, 后来发展的光学微腔-半导体量子点体系制备出来的光子全同度更高, 更加适合大规模实现典型玻色采样[16].

    同时, 研究者们也对光子损失玻色采样的理论和实验开展了深入的研究. 光子损失是指光子在产生、传输、探测过程中发生丢失的情况, 其与光子源、光学网络、探测器的质量密切相关. 与光子非全同一样, 光子损失会明显影响玻色采样的量子优势. 对于光子损失对玻色采样量子优势影响的分析, 主要基于两种理论模型. 一种模型是对输出端损失光子的个数进行统计分析. 例如, 对于输入光子数为n, 输出光子数为(nl): 若l为一个常数, 则玻色采样仍然具有量子优势; 若(nl)<O(n), 则玻色采样失去量子优势[17]. 更一般地, 随着输入光子数的增加: 若损失光子数与输入光子数之比趋向于0, 则玻色采样仍然具有量子优势; 若损失光子数与输入光子数之比趋向于1, 则玻色采样失去量子优势[13]. 第二种模型则是考虑光子损失有可能发生在每一个光学元件中, 其更加符合实际实验模型[18]. 这种模型将光子损失等效为光子通过一个虚构的分束器进入到环境光路中. 基于此, 对光子损失的玻色采样模拟实际上等效于对更大规模光学网络的玻色采样模拟. 这种模型揭示了光子损失玻色采样的量子优势与光学网络的深度密切相关[17].

    本工作基于光子损失情况下玻色采样的第二种理论模型展开, 对不同结构的光学网络模型进行了定量的分析, 给出了直观的比较结果. 同时本文的模拟结果由于考虑了光子损失概率的随机性, 更加符合真实物理场景.

    目前比较常用的两种光学线性网络, 即Reck结构和Clements结构的光学网络, 如图1所示, 图中每一个交点都表示一个光学单元.

    图 1 光学线性网络示意图: (a) Reck结构; (b) Clements结构\r\nFig. 1. Sketches of optical linear networks: (a) the Reck structure; (b) the Clements structure.
    图 1  光学线性网络示意图: (a) Reck结构; (b) Clements结构
    Fig. 1.  Sketches of optical linear networks: (a) the Reck structure; (b) the Clements structure.

    其中, Reck结构为三角形非对称结构, 而Clements结构为四边形对称结构[19]. 直接对任意m维幺正矩阵进行分解, 使其可以表示为一系列矩阵相乘的形式. 在该系列矩阵中, 第一个矩阵是对角矩阵, 对应光学网络末端的m个相移器; 其余矩阵均对应一个光学单元, 每个光学单元由一个相移器和一个分束器组成. 每个光学单元对应的幺正矩阵M可以表示为分束器矩阵MBS和相移器矩阵MPS的乘积[19]:

    M=MBSMPS=[10001cosωsinωsinωcosω10001][10001eiφ110001]
    (1)

    其中分束器的反射率是cosω, 相移器改变的相位是φ, 分束器矩阵含参元素所在行(列)数为分束器所在的两个通道数, 相移器矩阵含参元素所在行(列)数为相移器所在的通道数.

    对于任意的幺正矩阵U, 将其分解成一系列矩阵相乘的形式, 该系列矩阵对应的光学单元的排列符合Reck结构, 其分解方法如算法1所示[20].

    算法1 任意幺正矩阵符合Reck结构的分解算法[20]
    输入 任意幺正矩阵U
    输出 符合Reck结构的矩阵列表M_list
    1: M_list ← ϕ //初始化M_list为空集
    2: for k ← 1 to (m – 1) do //mU的维数
    3: for t ← 0 to (k – 1) do
    4:  ω ← arctan(|um–t, k–t/um–t, k–t+1|) //计算分束器参数ω, uU中元素, 下标分别表示u所在的行数和列数
    5:  temp ← (um–t, k–t/um–t, k–t+1)cotω //计算临时变量temp
    6:  φ ← arctan(Imag(temp)/Real(temp)) //计算相移器参数φ, Imag和Real分别表示temp的虚部和实部
    7:  if |um–t, k–t+1sinω – eiφum–t, k–tcosω| ≠ 0 then
    8:   ω –ω //修正ω符号, 使得步骤10的消元能够顺利进行
    9:  M ← OpticalUnit(m, k – t, k – t + 1, ω, φ)    //利用(1)式计算m维光学单元矩阵, 其中分束器所在通道数为(kt)和(kt + 1), 相移器所在通道数为(kt)
    10: M_list ← (M_list, M), U UM–1  //对矩阵U进行消元
    11: end for
    12: end for
    13: M_list (M_list, U) //将消元后得到的对角矩阵放入M_list
    14: M_list←Reverse(M_list) //反向排列M_list, 使得M_list中所有元素乘积为待分解矩阵
    下载: 导出CSV 
    | 显示表格

    对于Clements结构, U的分解方法如算法2所示[19]. Clements结构的光学网络对应的矩阵构造方法要比Reck结构复杂得多. 这主要是因为, U需要交替右乘和左乘一系列矩阵来实现对角化, 即MMMUM1M1M1=D, 其中D为对角矩阵, M表示一系列参数(通道数, 反射率以及相位)不尽相同的光学单元矩阵, 则U可以表示为形如M1M1M1DMMM的形式. 在此基础上, 可以初始化一个对角矩阵D'以及一个光学单元矩阵M', 使得DM=M1D. 经过计算可以发现, D'M'的参数是唯一的, 则U可以表示为形如M–1M–1···M–1D'M'MM···M的形式, 此时该形式中对角矩阵左乘的矩阵减少1个, 而右乘的矩阵增加1个. 通过重复构造D'M', 最终U可以表示为形如D'M'M'···M'MM···M的形式, 实现基于Clements结构对U的分解.

    算法2 任意幺正矩阵符合Clements结构的分解算法[19]
    输入 任意幺正矩阵U
    输出 符合Clements结构的矩阵列表M_list
    1: M_list ← M_list1 M_list2M_list3ω_list←ϕ   //初始化M_list为空集, 并初始化辅助集合M_list1, M_list2, M_list3以及ω_list为空集
    2: for k ← 1 to (m – 1) do //mU的维数
    3: for t← 0 to (k – 1) do
    4:  if k mod 2 ≠ 0 then
    5:   Compute(ω, φ) //根据算法1计算ω, φ
    6:   M ← OpticalUnit(m, kt, kt + 1, ω, φ)
    7:   M_list1 ← (M_list1, M), U UM–1  //若k为奇数, 则利用M–1右乘U, 对U进行消元
    8:  else
    9:    ω ← –arctan(|umk+t+1, t+1/umk+t, t+1|)
    10:    temp ← –(umk+t+1, t+1/umk+t, t+1)cotω
    11:   φ ← arctan(Imag(temp)/Real(temp))
    12:    if |umk+t+1, t+1cosω + eiφumk+t, t+1sinω| ≠ 0 then
    13:     ω ← –ω //修正ω符号, 使得步骤15的消元能够顺利进行
    14:   M ← OpticalUnit(m, mk + t, mk + t + 1, ω, φ)
    15:   ω_list (ω_list, ω), M_list2 (M_list2, M), U MU //若k为偶数, 则利用M左乘U, 对U进行消元
    16: end for
    17: end for
    18: ω_list←Reverse(ω_list), M_list2←Reverse(M_list2), D U, p ← 1  //U经过消元变换为对角矩阵D, 待分解矩阵可以表示为形如M–1M–1···M–1DMM···M的形式
    19: for k ← (m – 1) to 1 by –1 do
    20: for t ← (k – 1) to 0 by –1 do
    21:  if k mod 2 = 0 then
    22:   ω ← |ωp| //计算M'的分束器参数ω, ωpω_list的第p个元素
    23:   temp ← –tanωcotωpdmk+t/dmk+t+1 //d表示D的对角元素, 下标表示d所在行(列)数
    24:   φ ← arctan(Imag(temp)/Real(temp)) //计算M'的相移器参数φ
    25:   if |dmk+tsinωp + eiφdmk+t+1cosωp| ≠ 0 then
    26:     ω ← –ω //修正ω符号
    27:     Compute(φ) //在得到修正的ω后, 根据步骤23和24计算φ, 使得DM=M1pD, 其中MpM_list2的第p个元素
    28:   M M' ← OpticalUnit(m, mk + t, mk + t + 1, ω, φ) //计算M'并赋值给M
    29:   M_list3(M_list3, M), DDM1pDM1 , pp + 1 //计算D'并赋值给 D
    30: end for
    31: end for
    32: M_list3←(M_list3, D) //将对角矩阵放入M_list3
    33: M_list1←Reverse(M_list1), M_list3←Reverse(M_list3)
    34: M_list←(M_list3, M_list1) //M_list中所有元素乘积为待分解矩阵
    下载: 导出CSV 
    | 显示表格

    在得到对应任意幺正矩阵的光学网络后, 为了研究在光学网络中发生的光子损失对玻色采样的影响, 基于等效分束器原理对光学网络进行重构. 2个虚构的分束器和2条环境光路被加入每个光学单元中. 虚构的分束器参数ω和单条实际光路在单个光学单元中的光子损失概率Ploss满足cosω =1 – (1 – Ploss)1/2关系[18]. 当光子在该光学单元中发生损失, 则其等效于光子通过虚构的分束器进入环境光路中. 由于两种光学网络都具有m(m – 1)/2个光学单元, 在重构的光学网络中, 一共附加了m(m – 1)个虚构分束器和相同数目的环境光路, 则重构的光学网络的模式数变为m2, 因此原始m维矩阵M需要变换为m2维矩阵, 即MEm(m–1), 其中Em(m – 1)m(m – 1)维单位矩阵. 相应地, 虚构的分束器矩阵也是m2维矩阵. 在得到重构后的光学网络后, 则可以计算考虑光子损失的玻色采样输出组合概率, 其方法如算法3所示, 在计算输出组合概率的过程中, 需要计算矩阵的积和式, 可以采用Ryser或者Glynn方法, 并利用格雷码以简化计算过程[21]. 计算了4光子16模式的玻色采样输出组合概率, 其中16维输入组合S = [1, 1, 1, 1, 0, 0, ···, 0]. 一共随机产生了100个16维随机幺正矩阵, 并对最终计算结果取平均值.

    算法3 考虑光子损失的玻色采样输出组合概率算法
    输入 输入组合S, 输出组合T, m维光学网络矩阵U, 单条实际光路在单个光学单元处的光子损失概率Ploss
    输出 输出组合概率Pout
    1: M_list ← Decompose(U) //按照Reck或Clements结构分解矩阵U
    2: A M1Em(m 1), M1M_list //将对角矩阵的维度扩展为m2
    3: for k ← 2 to (m(m – 1)/2 + 1) do
    4: Mk MkEm(m – 1), MkM_list  //将光学单元矩阵的维度扩展为m2
    5: ω ← arccos(1 – (1 – Ploss)1/2) //计算虚构分束器参数ω
    6: Compute(ch) //确定原始光学单元所在通道数ch和(ch + 1)
    7: B ← OpticalUnit(m2, ch, m + 2k – 3, ω, 0), B' ← OpticalUnit(m2, ch + 1, m + 2k – 2, ω, 0)   //构造单个光学单元中的2个虚构分束器矩阵
    8: A AMkBB' //构造考虑光子损失的等效光学网络矩阵
    9: end for
    10: S ← (S, [0]m(m–1)) //将输入组合维度变为m2
    11: Pout ← ∑T' | Perm(AS, (T, T' ))|2, T' ∈{0, 1}m(m–1) //Perm()表示矩阵积和式函数, T' 中所有元素之和为损失光子数
    下载: 导出CSV 
    | 显示表格

    基于玻色子生日悖论, 在只考虑输出组合中每个通道只有0或1个光子的情况时, 可以通过总变差距离来评估玻色采样的输出结果与理想情况的偏差程度[21]. 总变差距离可以表示为T412|PoutPideal|, 其中T4表示没有发生光子损失的任意16维输出组合, Pideal表示理想情况下归一化的输出组合概率. 为了便于与理想情况对比, Pout也需要进行归一化处理, 使得对于所有的T4, Pout之和为1. 如图2(a)所示, 我们发现随着光子损失概率的提高, 总变差距离以准线性的方式逐步变大, 这表明光子损失会对玻色采样的精度产生显著影响. 我们可以直观地发现, 对于固定的Ploss, Clements结构比Reck结构对应的总变差距离更小, 表明Clements结构在光子损失情况下稳定性更好. 这主要是因为Clements结构受到平衡与非平衡光子损失的影响更小[19]. 在光子损失情况下, 一方面, 对于平衡光子损失, 光学网络深度越大, 光子损失程度越高. 对于较大规模光学网络, Clements结构的深度比Reck结构小很多, 因此受到平衡光子损失的影响更小; 另一方面, 对于非平衡光子损失, 光子损失程度与通道中分束器的光子损失概率密切相关. 在Reck结构中, 由于每个通道经过的分束器数目不同, 导致通道中光子损失程度不同, 这极大地影响了多光子相干干涉结果. 而在Clements结构中, 每个通道经过的分束器数目大致相同, 通道中光子损失程度大致相当, 因此在Clements结构中非平衡光子损失对多光子相干干涉结果的影响相比于Reck结构较小.

    图 2 (a) 考虑光子损失的玻色采样输出组合概率的总变差距离随光子损失概率Ploss的变化关系; (b) Ploss为0.1时, 两种光学网络结构随着模式数增加对应玻色采样结果总变差距离的变化关系. 图中空心图形代表实验结果平均值, 误差棒代表标准差\r\nFig. 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.
    图 2  (a) 考虑光子损失的玻色采样输出组合概率的总变差距离随光子损失概率Ploss的变化关系; (b) Ploss为0.1时, 两种光学网络结构随着模式数增加对应玻色采样结果总变差距离的变化关系. 图中空心图形代表实验结果平均值, 误差棒代表标准差
    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.

    为了验证上述结论, 我们对不同规模的玻色采样进行了仿真, 如图2(b)所示, 所采用玻色采样规模分别为1光子3模式、2光子6模式、3光子9模式、3光子12模式和4光子16模式. 在模式数为3时, 两种结构具有同样的形式, 因此总变差距离相同. 随着模式数的增加, 光学网络深度增加, 通道经过的最大分束器数目增加, 导致两种结构对应的总变差距离也逐渐增加, 表明光子损失(包括平衡与非平衡光子损失)对玻色采样结果的影响逐渐加大. 而Reck结构对应的总变差距离增幅明显高于Clements结构, 表明Clements结构受到平衡与非平衡光子损失的影响相对于Reck结构更小.

    进一步假设每个光学单元对应的Ploss不是一个定值, 而是符合一定的随机分布, 即正态分布. 可认为这样的假设比较能够反映出真实的物理实验情况. 分两种情况进行了模拟: 第1种情况是, 固定Ploss的标准差不变, 改变其平均值. 对于这种情况, 发现总变差距离的变化情况与图2(a)类似; 第2种情况是, 固定Ploss的平均值不变, 改变其标准差. 研究发现, 对于两种光学网络结构, 随着Ploss标准差的增大, 输出组合概率的总变差距离呈现变大的趋势, 如图3所示, 表明光子损失概率的离散程度也会明显影响玻色采样的精度. 以深度为3的光学网络为例, 对这种情况进行简单的分析. 此时光学网络的每一层只含有一个光学单元. 在Ploss为一个定值P时, 每一层的光子透过率为T = 1 – P, 因此总光子透过率是T 3. 若Ploss符合平均值为P的正态分布, 则每一层的光子透过率也符合平均值为T正态分布, 总光子透过率的期望值仍为T 3. 为了满足这个条件, 随着Ploss标准差的增加, 总光子透过率小于T 3的概率需要逐渐增加. 受较低的总光子透过率的影响, 输出组合概率的总变差距离相应增加, 表明光子损失对玻色采样精度的影响变得更为显著.

    图 3 考虑光子损失的玻色采样输出组合概率的总变差距离平均值随光子损失概率的标准差σPloss的变化关系\r\nFig. 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.
    图 3  考虑光子损失的玻色采样输出组合概率的总变差距离平均值随光子损失概率的标准差σPloss的变化关系
    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.

    最后, 统计了固定光子个数的输出组合的概率之和随着光子损失概率变化的情况, 如图4所示. 当Ploss = 0.01时, Clements结构相比于Reck结构, 4光子输出的概率明显较高. 由于输出光子数越高, 其对应的经典计算的时间复杂度也越高, 我们的结果表明在光子损失情况下, Clements结构的玻色采样对于经典计算来说具有更高的时间复杂度. 而随着Ploss的提高, 4光子输出概率急剧减小, 3光子输出占据优势. 进一步提高Ploss, 则2光子输出变成主导. 这表明, 发生在光学网络中的光子损失能够显著地影响玻色采样的量子优势. 当低光子数的输出组合占据主导时, 玻色采样变得易于被经典算法模拟, 如MCMC算法[22]、Clifford-Clifford算法[23]等.

    图 4 玻色采样固定光子个数的输出组合概率之和随着Ploss的变化关系. (a) Reck结构; (b) Clements结构光学网络. 图中空心图形代表实验结果平均值, 误差棒代表标准差\r\nFig. 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.
    图 4  玻色采样固定光子个数的输出组合概率之和随着Ploss的变化关系. (a) Reck结构; (b) Clements结构光学网络. 图中空心图形代表实验结果平均值, 误差棒代表标准差
    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.

    基于对光子损失情况下玻色采样的经典仿真结果, 发现光子损失能够显著影响玻色采样结果的精度, 进而影响玻色采样量子优势的可实现性. 如何抑制光子损失对玻色采样机等光量子器件的影响一直是研究热点之一. 在理论方面, 系统研究了输出光子个数[24]、光学网络深度[17]等对玻色采样对应经典计算时间复杂度的影响. 在工程技术方面, 减反射涂层[16]、光路拓展[25]、三维光学网络[26]的应用提供了抑制光子损失影响的可能性. 值得一提的是, 在容错全光量子中继器方案中, 利用量子奇偶码有效抑制了光子损失的影响, 这给未来光量子计算器件抵抗光子损失带来了启示[27].

    系统研究了两种结构的光学网络中光子损失对玻色采样的影响. 可以直观地发现, 当存在光子损失时, Clements结构相比于Reck结构具有更好的稳定性, 其输出高光子个数的概率更高, 对应经典计算的时间复杂度更高. 随着光子损失概率的提高, 玻色采样的精度明显降低, 且输出高光子个数的概率急剧降低, 表明光子损失对玻色采样的经典计算时间复杂度具有十分显著的影响. 研究还发现, 若光学单元中光子损失概率不是一个定值, 而是服从正态分布, 则光子损失概率的离散程度也能够影响玻色采样的精度. 结果表明, 发生在光学网络中的光子损失能够明显影响玻色采样的精度和对应经典计算的时间复杂度, 这对真正意义上具有量子优势的玻色采样的物理实现提出了又一个挑战.

    [1]

    Shor P W 1994 Proceedings 35th Annual Symposium on Foundations of Computer Science Santa Fe, NM, USA, IEEE p124

    [2]

    Grover L K 1997 Phys. Rev. Lett. 79 325Google Scholar

    [3]

    Spring J B, Metcalf B J, Humphreys P C, Kolthammer W S, Jin X M, Barbieri M, Datta A., Thomas N, Langford N K, Kundys D, Gates J C, Smith B J, Smith P G, Walmsley I A 2013 Science 339 798Google Scholar

    [4]

    Zhong H S, Wang H, Deng Y H, Chen M C, Peng L C, Luo Y H, Qin J, Wu D, Ding X, Hu Y H, Yang X Y, Zhang W J, Li H, Li Y X, Jiang X, Gan L, Yang G G, You L X, Wang Z, Li L, Liu N L, Lu C Y, Pan J W 2020 Science 370 1460Google Scholar

    [5]

    Zhong H S, Li Y, Li W, Peng L C, Su Z E, Hu Y, He Y M, Ding X, Zhang W J, Li H, Zhang L, Wang Z, You L X, Wang X L, Jiang X, Li L, Chen Y A, Liu N L, Lu C Y, Pan J W 2018 Phys. Rev. Lett. 121 250505Google Scholar

    [6]

    Wang H, Li W, Jiang X, He Y M, Li Y H, Ding X, Chen M C, Qin J, Peng C Z, Schneider C, Kamp M, Zhang W J, Li H, You LX, Wang Z, Dowling J P, Höfling S, Lu C Y, Pan J W 2018 Phys. Rev. Lett. 120 230502Google Scholar

    [7]

    Peropadre B, Guerreschi G G, Huh J, Aspuru-Guzik A 2016 Phys. Rev. Lett. 117 140505Google Scholar

    [8]

    Aaronson S, Arkhipov A 2011 Proceedings of the forty-third annual ACM symposium on Theory of computing San Jose, California, June 6–8, 2011 p333

    [9]

    Spagnolo N, Vitelli C, Bentivegna M, Brod D J, Crespi A, Flamini F, Giacomini S, Milani G, Ramponi R, Mataloni P, Osellame R, Galvão E F, Sciarrino F 2014 Nat. Photonics 8 615Google Scholar

    [10]

    Bentivegna M, Spagnolo N, Vitelli C, Brod D J, Crespi A, Flamini F, Ramponi R, Mataloni P, Osellame R, Galvão E F, Sciarrino F 2015 Int. J. Quantum. Inf. 12 1560028Google Scholar

    [11]

    Agresti I, Viggianiello N, Flamini F, Spagnolo N, Crespi A, Osellame R, Wiebe N, Sciarrino F 2017 Phys. Rev. X 9 2160Google Scholar

    [12]

    Huh J, Guerreschi G G, Peropadre B, Mcclean J R, Aspuru-Guzik A 2015 Nat. Photonics 9 615Google Scholar

    [13]

    Renema1 J J, Shchesnovich V, Garcia Patron R 2019 arXiv: 1809.01953 v2 [quant-ph]

    [14]

    Tichy C M 2015 Phys. Rev. A 91 022316Google Scholar

    [15]

    Renema J J, Menssen A, Clements W R, Triginer G, Kolthammer W S, Walmsley I A 2018 Phys. Rev. Lett. 120 220502Google Scholar

    [16]

    Wang H, He Y, Li Y H, Su Z E, Li B, Huang H L, Ding X, Chen M C, Liu C, Qin J, Li J P, He Y M, Schneider C, Kamp M, Peng C Z, Höfling S, Lu C Y, Pan J W 2017 Nat. Photonics 11 361Google Scholar

    [17]

    García-Patrón R, Renema J J, Shchesnovich V 2019 Quantum 3 169Google Scholar

    [18]

    Oszmaniec M, Brod D J 2018 New J. Phys. 20 092002Google Scholar

    [19]

    Clements W R, Humphreys P C, Metcalf B J, Steven K W, Walsmley I A 2016 Optica 3 1460Google Scholar

    [20]

    Reck M, Zeilinger A, Bernstein H J, Bertani P 1994 Phys. Rev. Lett. 73 58Google Scholar

    [21]

    Brod D J, Galvão E F, Crespi A, Osellame R, Spagnolo N, Sciarrino F 2019 Adv. Photon. 1 034001Google Scholar

    [22]

    Neville A, Sparrow C, Clifford R, Johnston E, Birchall P M, Montanaro A, Laing A 2017 Nat. Phys. 13 1153Google Scholar

    [23]

    Clifford P, Clifford R 2018 Proceedings of the Twenty-Ninth Annual ACM-SIAM Symposium on Discrete Algorithms New Orleans Louisiana January 7–10, 2018 p146

    [24]

    Aaronson S, Brod D J 2015 Phys. Rev. A 93 012335Google Scholar

    [25]

    Miller D A B 2015 Optica 2 8Google Scholar

    [26]

    Wang H, Qin J, Ding X, Chen M C, Chen S, You X, He Y M, Jiang X, You L, Wang Z, Schneider C, Renema J J, Höfling S, Lu C Y, Pan J W 2019 Phys. Rev. Lett. 123 250503Google Scholar

    [27]

    Zhang R, Liu L Z, Li Z D, Fei Y Y, Yin X F, Li L, Liu N L, Mao Y, Chen Y A, Pan J W 2022 Optica 9 2Google Scholar

  • 图 1  光学线性网络示意图: (a) Reck结构; (b) Clements结构

    Figure 1.  Sketches of optical linear networks: (a) the Reck structure; (b) the Clements structure.

    图 2  (a) 考虑光子损失的玻色采样输出组合概率的总变差距离随光子损失概率Ploss的变化关系; (b) Ploss为0.1时, 两种光学网络结构随着模式数增加对应玻色采样结果总变差距离的变化关系. 图中空心图形代表实验结果平均值, 误差棒代表标准差

    Figure 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.

    图 3  考虑光子损失的玻色采样输出组合概率的总变差距离平均值随光子损失概率的标准差σPloss的变化关系

    Figure 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.

    图 4  玻色采样固定光子个数的输出组合概率之和随着Ploss的变化关系. (a) Reck结构; (b) Clements结构光学网络. 图中空心图形代表实验结果平均值, 误差棒代表标准差

    Figure 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.

    算法1 任意幺正矩阵符合Reck结构的分解算法[20]
    输入 任意幺正矩阵U
    输出 符合Reck结构的矩阵列表M_list
    1: M_list ← ϕ //初始化M_list为空集
    2: for k ← 1 to (m – 1) do //mU的维数
    3: for t ← 0 to (k – 1) do
    4:  ω ← arctan(|um–t, k–t/um–t, k–t+1|) //计算分束器参数ω, uU中元素, 下标分别表示u所在的行数和列数
    5:  temp ← (um–t, k–t/um–t, k–t+1)cotω //计算临时变量temp
    6:  φ ← arctan(Imag(temp)/Real(temp)) //计算相移器参数φ, Imag和Real分别表示temp的虚部和实部
    7:  if |um–t, k–t+1sinω – eiφum–t, k–tcosω| ≠ 0 then
    8:   ω –ω //修正ω符号, 使得步骤10的消元能够顺利进行
    9:  M ← OpticalUnit(m, k – t, k – t + 1, ω, φ)    //利用(1)式计算m维光学单元矩阵, 其中分束器所在通道数为(kt)和(kt + 1), 相移器所在通道数为(kt)
    10: M_list ← (M_list, M), U UM–1  //对矩阵U进行消元
    11: end for
    12: end for
    13: M_list (M_list, U) //将消元后得到的对角矩阵放入M_list
    14: M_list←Reverse(M_list) //反向排列M_list, 使得M_list中所有元素乘积为待分解矩阵
    DownLoad: CSV
    算法2 任意幺正矩阵符合Clements结构的分解算法[19]
    输入 任意幺正矩阵U
    输出 符合Clements结构的矩阵列表M_list
    1: M_list ← M_list1 M_list2M_list3ω_list←ϕ   //初始化M_list为空集, 并初始化辅助集合M_list1, M_list2, M_list3以及ω_list为空集
    2: for k ← 1 to (m – 1) do //mU的维数
    3: for t← 0 to (k – 1) do
    4:  if k mod 2 ≠ 0 then
    5:   Compute(ω, φ) //根据算法1计算ω, φ
    6:   M ← OpticalUnit(m, kt, kt + 1, ω, φ)
    7:   M_list1 ← (M_list1, M), U UM–1  //若k为奇数, 则利用M–1右乘U, 对U进行消元
    8:  else
    9:    ω ← –arctan(|umk+t+1, t+1/umk+t, t+1|)
    10:    temp ← –(umk+t+1, t+1/umk+t, t+1)cotω
    11:   φ ← arctan(Imag(temp)/Real(temp))
    12:    if |umk+t+1, t+1cosω + eiφumk+t, t+1sinω| ≠ 0 then
    13:     ω ← –ω //修正ω符号, 使得步骤15的消元能够顺利进行
    14:   M ← OpticalUnit(m, mk + t, mk + t + 1, ω, φ)
    15:   ω_list (ω_list, ω), M_list2 (M_list2, M), U MU //若k为偶数, 则利用M左乘U, 对U进行消元
    16: end for
    17: end for
    18: ω_list←Reverse(ω_list), M_list2←Reverse(M_list2), D U, p ← 1  //U经过消元变换为对角矩阵D, 待分解矩阵可以表示为形如M–1M–1···M–1DMM···M的形式
    19: for k ← (m – 1) to 1 by –1 do
    20: for t ← (k – 1) to 0 by –1 do
    21:  if k mod 2 = 0 then
    22:   ω ← |ωp| //计算M'的分束器参数ω, ωpω_list的第p个元素
    23:   temp ← –tanωcotωpdmk+t/dmk+t+1 //d表示D的对角元素, 下标表示d所在行(列)数
    24:   φ ← arctan(Imag(temp)/Real(temp)) //计算M'的相移器参数φ
    25:   if |dmk+tsinωp + eiφdmk+t+1cosωp| ≠ 0 then
    26:     ω ← –ω //修正ω符号
    27:     Compute(φ) //在得到修正的ω后, 根据步骤23和24计算φ, 使得DM=M1pD, 其中MpM_list2的第p个元素
    28:   M M' ← OpticalUnit(m, mk + t, mk + t + 1, ω, φ) //计算M'并赋值给M
    29:   M_list3(M_list3, M), DDM1pDM1 , pp + 1 //计算D'并赋值给 D
    30: end for
    31: end for
    32: M_list3←(M_list3, D) //将对角矩阵放入M_list3
    33: M_list1←Reverse(M_list1), M_list3←Reverse(M_list3)
    34: M_list←(M_list3, M_list1) //M_list中所有元素乘积为待分解矩阵
    DownLoad: CSV
    算法3 考虑光子损失的玻色采样输出组合概率算法
    输入 输入组合S, 输出组合T, m维光学网络矩阵U, 单条实际光路在单个光学单元处的光子损失概率Ploss
    输出 输出组合概率Pout
    1: M_list ← Decompose(U) //按照Reck或Clements结构分解矩阵U
    2: A M1Em(m 1), M1M_list //将对角矩阵的维度扩展为m2
    3: for k ← 2 to (m(m – 1)/2 + 1) do
    4: Mk MkEm(m – 1), MkM_list  //将光学单元矩阵的维度扩展为m2
    5: ω ← arccos(1 – (1 – Ploss)1/2) //计算虚构分束器参数ω
    6: Compute(ch) //确定原始光学单元所在通道数ch和(ch + 1)
    7: B ← OpticalUnit(m2, ch, m + 2k – 3, ω, 0), B' ← OpticalUnit(m2, ch + 1, m + 2k – 2, ω, 0)   //构造单个光学单元中的2个虚构分束器矩阵
    8: A AMkBB' //构造考虑光子损失的等效光学网络矩阵
    9: end for
    10: S ← (S, [0]m(m–1)) //将输入组合维度变为m2
    11: Pout ← ∑T' | Perm(AS, (T, T' ))|2, T' ∈{0, 1}m(m–1) //Perm()表示矩阵积和式函数, T' 中所有元素之和为损失光子数
    DownLoad: CSV
  • [1]

    Shor P W 1994 Proceedings 35th Annual Symposium on Foundations of Computer Science Santa Fe, NM, USA, IEEE p124

    [2]

    Grover L K 1997 Phys. Rev. Lett. 79 325Google Scholar

    [3]

    Spring J B, Metcalf B J, Humphreys P C, Kolthammer W S, Jin X M, Barbieri M, Datta A., Thomas N, Langford N K, Kundys D, Gates J C, Smith B J, Smith P G, Walmsley I A 2013 Science 339 798Google Scholar

    [4]

    Zhong H S, Wang H, Deng Y H, Chen M C, Peng L C, Luo Y H, Qin J, Wu D, Ding X, Hu Y H, Yang X Y, Zhang W J, Li H, Li Y X, Jiang X, Gan L, Yang G G, You L X, Wang Z, Li L, Liu N L, Lu C Y, Pan J W 2020 Science 370 1460Google Scholar

    [5]

    Zhong H S, Li Y, Li W, Peng L C, Su Z E, Hu Y, He Y M, Ding X, Zhang W J, Li H, Zhang L, Wang Z, You L X, Wang X L, Jiang X, Li L, Chen Y A, Liu N L, Lu C Y, Pan J W 2018 Phys. Rev. Lett. 121 250505Google Scholar

    [6]

    Wang H, Li W, Jiang X, He Y M, Li Y H, Ding X, Chen M C, Qin J, Peng C Z, Schneider C, Kamp M, Zhang W J, Li H, You LX, Wang Z, Dowling J P, Höfling S, Lu C Y, Pan J W 2018 Phys. Rev. Lett. 120 230502Google Scholar

    [7]

    Peropadre B, Guerreschi G G, Huh J, Aspuru-Guzik A 2016 Phys. Rev. Lett. 117 140505Google Scholar

    [8]

    Aaronson S, Arkhipov A 2011 Proceedings of the forty-third annual ACM symposium on Theory of computing San Jose, California, June 6–8, 2011 p333

    [9]

    Spagnolo N, Vitelli C, Bentivegna M, Brod D J, Crespi A, Flamini F, Giacomini S, Milani G, Ramponi R, Mataloni P, Osellame R, Galvão E F, Sciarrino F 2014 Nat. Photonics 8 615Google Scholar

    [10]

    Bentivegna M, Spagnolo N, Vitelli C, Brod D J, Crespi A, Flamini F, Ramponi R, Mataloni P, Osellame R, Galvão E F, Sciarrino F 2015 Int. J. Quantum. Inf. 12 1560028Google Scholar

    [11]

    Agresti I, Viggianiello N, Flamini F, Spagnolo N, Crespi A, Osellame R, Wiebe N, Sciarrino F 2017 Phys. Rev. X 9 2160Google Scholar

    [12]

    Huh J, Guerreschi G G, Peropadre B, Mcclean J R, Aspuru-Guzik A 2015 Nat. Photonics 9 615Google Scholar

    [13]

    Renema1 J J, Shchesnovich V, Garcia Patron R 2019 arXiv: 1809.01953 v2 [quant-ph]

    [14]

    Tichy C M 2015 Phys. Rev. A 91 022316Google Scholar

    [15]

    Renema J J, Menssen A, Clements W R, Triginer G, Kolthammer W S, Walmsley I A 2018 Phys. Rev. Lett. 120 220502Google Scholar

    [16]

    Wang H, He Y, Li Y H, Su Z E, Li B, Huang H L, Ding X, Chen M C, Liu C, Qin J, Li J P, He Y M, Schneider C, Kamp M, Peng C Z, Höfling S, Lu C Y, Pan J W 2017 Nat. Photonics 11 361Google Scholar

    [17]

    García-Patrón R, Renema J J, Shchesnovich V 2019 Quantum 3 169Google Scholar

    [18]

    Oszmaniec M, Brod D J 2018 New J. Phys. 20 092002Google Scholar

    [19]

    Clements W R, Humphreys P C, Metcalf B J, Steven K W, Walsmley I A 2016 Optica 3 1460Google Scholar

    [20]

    Reck M, Zeilinger A, Bernstein H J, Bertani P 1994 Phys. Rev. Lett. 73 58Google Scholar

    [21]

    Brod D J, Galvão E F, Crespi A, Osellame R, Spagnolo N, Sciarrino F 2019 Adv. Photon. 1 034001Google Scholar

    [22]

    Neville A, Sparrow C, Clifford R, Johnston E, Birchall P M, Montanaro A, Laing A 2017 Nat. Phys. 13 1153Google Scholar

    [23]

    Clifford P, Clifford R 2018 Proceedings of the Twenty-Ninth Annual ACM-SIAM Symposium on Discrete Algorithms New Orleans Louisiana January 7–10, 2018 p146

    [24]

    Aaronson S, Brod D J 2015 Phys. Rev. A 93 012335Google Scholar

    [25]

    Miller D A B 2015 Optica 2 8Google Scholar

    [26]

    Wang H, Qin J, Ding X, Chen M C, Chen S, You X, He Y M, Jiang X, You L, Wang Z, Schneider C, Renema J J, Höfling S, Lu C Y, Pan J W 2019 Phys. Rev. Lett. 123 250503Google Scholar

    [27]

    Zhang R, Liu L Z, Li Z D, Fei Y Y, Yin X F, Li L, Liu N L, Mao Y, Chen Y A, Pan J W 2022 Optica 9 2Google Scholar

Metrics
  • Abstract views:  4850
  • PDF Downloads:  57
Publishing process
  • Received Date:  25 February 2022
  • Accepted Date:  06 June 2022
  • Available Online:  16 September 2022
  • Published Online:  05 October 2022

/

返回文章
返回