1. 引言

量子纠缠的产生是量子信息处理(QIP) [1] 的一项基本任务，已在许多量子系统中得以实现，如在光子 [2] [3] [4] 、核磁共振 [5] [6] [7] 、金刚石氮空位中心 [8] [9] [10] [11] [12] 和腔量子电动力学(QED)中 [13] 。由于原子与光子之间一旦实现强耦合便可以对原子或光子进行快速QIP，所以人们对原子与单模光场组成的腔QED系统进行了大量研究。人工量子系统因为其具有大规模集成的良好潜力所以在模拟腔QED方面引起了广泛关注。

电路QED [14] [15] 作为一种人工量子系统，由人工原子超导量子比特(SQ)和超导谐振器(SR)组成，基于SQ或SR的通用量子计算研究也已经取得了很大进展。SQ方面，SQ一般由电容、电感和约瑟夫森结构成，其能量弛豫时间已提高到50 µs [16] 左右，近些年人们正在试图实现基于其的容错量子计算。例如，Chow等人 [17] 在2014年实现了一种可扩展的容错量子计算结构。Martinis等人 [18] 在2015年通过在一个由9个SQ组成的线性阵列的超导量子电路中，通过重复误差检测实现了量子状态的保持。SR方面，SR是基于超导量子电路的重要器件，其在存储和引导微波光子方面有重要作用。除此之外，SR还可以在超导人造原子间传递信息，充当量子数据线的作用。近些年，SR的品质因子已被提高到10⁶[19] ，这意味着该SR可以作为一种良好的量子信息载体。已经有大量实验和理论工作在研究基于SR的QIP，例如，Wang等人 [20] 在2011年实现了两个SR的NOON态纠缠。Taylor等人 [21] 在2013年提出了一种利用电路QED实现非线性光学量子计算的方案。与此同时，基于SR的通用量子门构造 [22] - [27] 、纠缠产生 [28] - [33] 和测量 [34] [35] [36] 等量子计算的基本任务也有很多研究工作。

基于电路量子电动力学的量子处理器的主要集成方式包括：一些SQ和一个SR [37] 的耦合，一些SR与一个SR总线 [38] 或一些SR与一个SQ [39] [40] [41] 的耦合，以及一些电路QED通过SR总线 [42] 的耦合等。在本文中，我们关注另一种集成多个SR或SQ的芯片类型——分布式量子计算(DQC) [43] ，它需要一条长传输线(TL)来连接许多本地处理器，以形成一个量子网络。在DQC中，量子器件之间的串扰可以进一步减少。为了实现远距离量子器件上的纠缠操作，人们需要在TL中使用真实的光子传播或暗光子的绝热演化。在电路QED中，为了连接多个处理器，需要对超导TL进行复杂的设计，从而降低TL的品质因子，提高了其中微波光子的泄漏率。因此，基于电路QED的分布式量子计算需要考虑超导TL中的暗光子的绝热演化过程。

本文是基于非局域性的超导量子谐振器的纠缠研究。在一个长超导TL连接两个SR，且这两个SR又分别与两个承载量子信息的SR耦合的系统中，提出了一种基于两个非局域超导量子SR的贝尔态产生方案。通过让所有SR和超导TL的频率相等，本方案可以使超导TL不被布局微波光子从而克服TL泄漏的影响。通过考虑可行的实验参数，贝尔态的保真度可达到99%左右。

2. 系统模型及其哈密顿量

为了构建基于两个远程超导谐振器的贝尔态纠缠产生方案。我们考虑系统如图1所示，该系统由长超导传输线(r_f)连接两个远距离一维超导谐振器(R_a1和R_a2)，且R_a1和R_a2又分别于两个承载量子信息的超导谐振器(b₁和b₂)耦合组成。

考虑到系统中不同原器件之间的相互作用，系统的哈密顿量在相互作用绘景下可以写为：

$\begin{array}{l}{\hat{H}}_1=g_1^{a1}\left({\hat{a}}_1^{+}{\hat{b}}_1{\text{e}}^{-i{\delta }_1^{a1}t}+{\hat{a}}_1{\hat{b}}_1^{+}{\text{e}}^{i{\delta }_1^{a1}t}\right)+g_2^{a2}\left({\hat{a}}_2^{+}{\hat{b}}_2{\text{e}}^{-i{\delta }_2^{a2}t}+{\hat{a}}_2{\hat{b}}_2^{+}{\text{e}}^{i{\delta }_2^{a2}t}\right)\\ +\underset{j=1}{\overset{\infty }{{\displaystyle \sum }}}{g}_{f,j}^I\left[{\hat{f}}_j\left({\hat{a}}_1^{+}+{\left(-1\right)}^j{\text{e}}^{i\varphi }{\hat{a}}_2^{+}\right)+H.c.\right]\end{array}$(1)

在公式(1)中， $g_1^{\hat{a}1}$是超导量子谐振器b₁与一维超导谐振器R_a1之间的耦合强度； $g_2^{a2}$是超导量子谐振器b₂与一维超导谐振器R_a2之间的耦合强度。其中频率失谐量 ${\delta }_{1\left(2\right)}^{a1\left(2\right)}={\omega }_{1\left(2\right)}-{\omega }_{a1\left(a2\right)}$， ${\omega }_{1\left(2\right)}$( ${\omega }_{a1\left(a2\right)}$)为谐振器b₁₍₂₎(R_a1(a2))的频率。 ${\hat{a}}_1^{+}$( ${\hat{a}}_2^{+}$)和 ${\hat{b}}_1^{+}$( ${\hat{b}}_2^{+}$)分别是谐振器R_a1(R_a2)和远程谐振器b₁(b₂)的产生算符。 $g_{f,j}^{\text{I}}$( $\text{I}=\text{a}、\text{b}$)是谐振器R_aI与传输线r_f在j模式下的耦合强度。 $\varphi$为传输场通过长度为L的传输线r_f所产生的相位，且满足公式： $\varphi =2\pi \omega l/c$，其中c为光速。利用精心设计的谐振器(R_a1)以及让共振操作时间远长与光子在TL中的往返时间，可以得到哈密顿量 ${\hat{H}}_1$的最后一项。

在选取 $2L{k}_f^{a1\left(a2\right)}/2\pi c\le 1$情形下TL的短极限时，我们只需考虑r_f的一个谐振模态f与R_a1和R_a2的谐振器相互作用。L是r_f的长度， $k_f^{a1\left(a2\right)}$是谐振器R_a1和R_a2的衰减率。这时，哈密顿量 ${\hat{H}}_1$可以约化为：

$\begin{array}{l}{\hat{H}}_2=g_1^{a_1}\left({\hat{a}}_1^{+}{\hat{b}}_1{\text{e}}^{-i{\delta }_1^{a1}t}+{\hat{a}}_1{\hat{b}}_1^{+}{\text{e}}^{i{\delta }_1^{a1}t}\right)+g_2^{a_2}\left({\hat{a}}_2^{+}{\hat{b}}_2{\text{e}}^{-i{\delta }_2^{a2}t}+{\hat{a}}_2{\hat{b}}_2^{+}{\text{e}}^{i{\delta }_2^{a2}t}\right)\\ +g_f^{a1}\left({\hat{f}}^{+}{\hat{a}}_1+\hat{f}{\hat{a}}_1^{+}\right)+g_f^{a2}\left({\hat{f}}^{+}{\hat{a}}_2+\hat{f}{\hat{a}}_2^{+}\right)\end{array}$(2)

当我们可取： ${\omega }_1/2\pi ={\omega }_2/2\pi ={\omega }_{a1}/2\pi ={\omega }_{a2}/2\pi ={\omega }_f/2\pi =\omega$，则系统哈密顿量 ${\hat{H}}_2$可以简化为：

${\hat{H}}_2=g_1^{a1}\left({\hat{a}}_1^{+}{\hat{b}}_1+{\hat{a}}_1{\hat{b}}_1^{+}\right)+g_2^{a2}\left({\hat{a}}_2^{+}{\hat{b}}_2+{\hat{a}}_2{\hat{b}}_2^{+}\right)+g_f^{a1}\left({\hat{f}}^{+}{\hat{a}}_1+\hat{f}{\hat{a}}_1^{+}\right)+g_f^{a2}\left({\hat{f}}^{+}{\hat{a}}_2+\hat{f}{\hat{a}}_2^{+}\right)$(3)

哈密顿量 ${\hat{H}}_2$在薛定谔绘景下可以表示成：

$\begin{array}{l}{\hat{H}}_3=\omega {\hat{a}}_1^{+}{\hat{a}}_1+\omega {\hat{a}}_2^{+}{\hat{a}}_2+\omega f^{+}f+\omega {\hat{b}}_1^{+}{\hat{b}}_1+\omega {\hat{b}}_2^{+}{\hat{b}}_2+g_1^{a_1}\left({\hat{a}}_1^{+}{\hat{b}}_1+{\hat{a}}_1{\hat{b}}_1^{+}\right)\\ +g_2^{a_2}\left({\hat{a}}_2^{+}{\hat{b}}_2+{\hat{a}}_2{\hat{b}}_2^{+}\right)+g_f^{a_1}\left({\hat{f}}^{+}{\hat{a}}_1+\hat{f}{\hat{a}}_1^{+}\right)+g_f^{a_2}\left({\hat{f}}^{+}{\hat{a}}_2+\hat{f}{\hat{a}}_2^{+}\right)\end{array}$(4)

通过取 $g_f^{a1}=g_f^{a2}=g$和规范变换 ${\hat{C}}_{\pm }=\frac{1}{2}\left({\hat{a}}_1+{\hat{a}}_2\pm \sqrt{2}\hat{f}\right)$和 $\hat{C}=\frac{\sqrt{2}}{2}\left({\hat{a}}_1-{\hat{a}}_2\right)$，哈密顿量 $H_3$可表示为：

$\begin{array}{l}{\hat{H}}_4=\omega {\hat{b}}_1^{+}{\hat{b}}_1+\omega {\hat{b}}_2^{+}{\hat{b}}_2+\omega {\hat{C}}^{+}\hat{C}+\left(\omega +\sqrt{2}g\right){\hat{C}}_{+}{\hat{C}}_{+}^{+}+\left(\omega -\sqrt{2}g\right){\hat{C}}_{-}{\hat{C}}_{-}^{+}\\ +\frac{1}{2}[g_1^{a1}\left({\hat{C}}_{+}+{\hat{C}}_{-}+\sqrt{2}\hat{C}\right){\hat{b}}_1^{+}+g_1^{a1}\left({\hat{C}}_{+}^{+}+{\hat{C}}_{-}^{+}+\sqrt{2}{\hat{C}}^{+}\right){\hat{b}}_1\\ +g_2^{a2}\left({\hat{C}}_{+}+{\hat{C}}_{-}-\sqrt{2}\hat{C}\right){\hat{b}}_2^{+}+g_2^{a2}\left({\hat{C}}_{+}^{+}+{\hat{C}}_{-}^{+}-\sqrt{2}{\hat{C}}^{+}\right){\hat{b}}_2]\end{array}$(5)

公式(5)中， $\hat{C}$和 ${\hat{C}}_{\pm }$是三种相互不耦合的玻色子模式。当R_a1、R_a2和r_f的频率相等时， ${\hat{C}}_{\pm }$和 $\hat{C}$的频率被分成三个不同的部分，频率分别为 ${\omega }_{C_{+}}=\omega +\sqrt{2}g$、 ${\omega }_{C_{-}}=\omega -\sqrt{2}g$和 ${\omega }_C=\omega$。当取 $g\gg \left\{g_1^{a1},g_2^{a2}\right\}$时，当 $\omega \pm \sqrt{2}g$与模式C、b₁和b₂的频率大失谐时，模式 $C_{\pm }$可以被忽略。因此，系统的哈密顿量可以简化为(相互作用绘景中)：

${\hat{H}}_{eff}=\frac{1}{\sqrt{2}}\left[g_1^{a1}\left(\hat{C}{\hat{b}}_1^{+}+{\hat{C}}^{+}{\hat{b}}_1\right)-g_2^{a2}\left(\hat{C}{\hat{b}}_2^{+}+{\hat{C}}^{+}{\hat{b}}_2\right)\right]$(6)

这里只剩下模式 $\hat{C}=\frac{\sqrt{2}}{2}\left({\hat{a}}_1-{\hat{a}}_2\right)$，这意味着在超导传输线r_f不参与整个演化过程，不会被布局微波光子。

假设哈密顿量 $H_{eff}$所描述的有效系统的初始状态为 $|{\varphi }_1={|1\rangle }_1{|0\rangle }_2{|0\rangle }_c$( ${|0}_c\equiv {|0\rangle }_a{|0\rangle }_b{|0\rangle }_f$)，系统的演化可以表示为：

$\begin{array}{l}|\varphi \left(t\right)={\text{e}}^{i{\hat{H}}_{eff}t}|{\varphi }_1\rangle \\ =\frac{1}{G}\left[{\left(g_2^{a2}\right)}^2+{\left(g_1^{a1}\right)}^2\cos \left(\sqrt{\frac{G}{2}}t\right)\right]|{\varphi }_1\rangle -\frac{g_1^{a1}{g}_2^{a2}}{G}\left[\cos \left(\sqrt{\frac{G}{2}}t\right)-1\right]|{\varphi }_2\rangle -\frac{i{g}_1^{a1}}{G}\sin \left(\sqrt{\frac{G}{2}}t\right)|{\varphi }_c\rangle \end{array}$(7)

其中， $G={\left(g_1^{a1}\right)}^2+{\left(g_2^{a2}\right)}^2$， $|{\varphi }_2\rangle ={|0\rangle }_1{|1\rangle }_2{|0\rangle }_c$和 $|{\varphi }_c\rangle ={|0\rangle }_1{|0\rangle }_2{|1\rangle }_c$。显然，如果取 $g_1^{a1}/g_2^{a2}=\sqrt{2}+1$，有效系统的状态可以演化为：

$|{\varphi }_{-}\rangle =\frac{1}{\sqrt{2}}\left({|0\rangle }_1{|1\rangle }_2-{|1\rangle }_1{|0\rangle }_2\right)\otimes {|0\rangle }_c$(8)

其操作时间为 $\sqrt{\frac{G}{2}}t=k\pi$( $k=1,2,3,\cdots$)。此外，如果取 $g_1^{a1}/g_2^{a2}=\sqrt{2}-1$，就可以得到另一个贝尔纠缠态：

$|{\varphi }_{+}\rangle =\frac{1}{\sqrt{2}}\left({|0\rangle }_1{|1\rangle }_2+{|1\rangle }_1{|0\rangle }_2\right)\otimes {|0\rangle }_c$(9)

其操作时间为 $\sqrt{\frac{G}{2}}t=k\pi$( $k=1,2,3,\cdots$)。

3. 数值模拟

为了研究系统在演化过程中每个器件的能量布局的变化，我们数值模拟了r_f中的能量布局随时间变化(如图2所示)、R_a1和R_a2中能量布局随时间变化(分别如图3(a)和图3(b)所示)、b₁和b₂中的能量布局随时间变化(分别如图3(c)和图3(d)所示)。从图2中可以看出，随着 $g/g_1^{a1}$越大，传输线中能量布局概率越小，即传输线的泄露率对整个演化过程的影响越小。

进一步，我们绘制了谐振器R_a1、R_a2、b₁和b₂的能量布局随时间的变化图3。

由图3可以发现谐振器b₁、b₂的能量布局图在100~200 ns区间内有一处b₁与b₂的峰值同时处于50%的位置，也就是说：在此时刻，谐振器b₁中50%的概率有一个光子，同时在b₂腔内也有50%的概率有一个光子。联合来看，当微波光子在b₁和b₂中均有一半概率存在时，在R_a1、R_a2和r_f中没有微波光子，那么此时谐振器b₁与b₂处于贝尔纠缠态 $|{\varphi }_{-}\rangle$。这里使用的能量在b₁、b₂、R_a1、R_a2和r_f中的布局 $P_m$( $m=1,2,\text{a}1,\text{a}2,\text{f}$)的定义为：

$P_m={\left|\langle {\varphi }_m|{\text{e}}^{-i{\hat{H}}_{1}t}|{\varphi }_1\rangle \right|}^2$(10)

其中 ${\varphi }_1={|1\rangle }_1{|0\rangle }_2{|0\rangle }_{a_1}{|0\rangle }_{a_2}{|0\rangle }_f$、 ${\varphi }_2={|0\rangle }_1{|1\rangle }_2{|0\rangle }_{a_1}{|0\rangle }_{a_2}{|0\rangle }_f$、 ${\varphi }_{a_1}={|0\rangle }_1{|0\rangle }_2{|1\rangle }_{a_1}{|0\rangle }_{a_2}{|0\rangle }_f$、 ${\varphi }_{a_2}={|0\rangle }_1{|0\rangle }_2{|0\rangle }_{a_1}{|1\rangle }_{a_2}{|0\rangle }_f$和 ${\varphi }_f={|0\rangle }_1{|0\rangle }_2{|0\rangle }_{a_1}{|0\rangle }_{a_2}{|1\rangle }_f$。在这里，我们取 $g_f^{a_1\left(a_2\right)}=20g_1^{a1\left(a2\right)}$可以完美的抑制 $C_{\pm }$(如图3(c)所示)。

最后我们需要计算我们的贝尔态的保真度，利用计算贝尔态的保真度的定义式：

$F={\left|\langle {\varphi }_m|{\text{e}}^{-i{\hat{H}}_{1}t}|{\varphi }_1\rangle \right|}^2$(11)

基于(11)式我们给出了态 $|{\varphi }_{-}\rangle$和 $|{\varphi }_{+}\rangle$的保真度随时间变化，分别如图4(a)和图4(b)所示。其中， $|{\varphi }_{-}\rangle$的保真度在131.7 ns可以达到99.8%， $|{\varphi }_{+}\rangle$的保真度在55.4 ns达到99.4%。

纠缠态是量子计算中的基础资源，无论是构造量子比特还是实现量子门操作都需要以其为基础。而在超导量子计算中，大量超导量子比特的集成会带来严重的串扰问题。为了规避串扰问题，我们可以考虑将分布式量子计算。而分布式计算中我们需要传输线连接各个超导谐振器，并且希望传输线中的能量布局概率越低越好，为此我们进行了一些数值模拟以保证我们的设计的贝尔态生成方案的可行性。

此外，无论是比特还是量子比特都依赖于超导量子电路或者半导体晶体管这些真实存在的元件。所以我们必须确保其正常运行的概率足够高，即量子态的保真度要足够高。所以，我们根据态保真度的定义计算我们方案中所制备的贝尔态的保真度，结果表明， $|{\varphi }_{-}\rangle$的保真度在131.7 ns可以达到99.8%， $|{\varphi }_{+}\rangle$的保真度在55.4 ns达到99.4%，结果较好，证明了我们的贝尔态生成方案的可行性。

4. 总结

我们基于非局域超导量子谐振器的设计处理贝尔态纠缠产生方案。我们首先构建了适当的系统，并根据系统中各器件及其相互作用关系确定了系统的原始哈密顿量，再引入特定参数关系，从原始哈密顿量获得有效哈密顿量。紧接着我们研究了系统在特定初态下的演化过程，并成功生成了基于两个非局域超导量子谐振器的贝尔纠缠态。而为了为验证我们方案的实验可行性，我们利用Python和QuTiP程序包进行了数值模拟。在采用真实实验参数的前提下，数值模拟的结果表明，通过我们的方案制备的贝尔态的保真度均在99%以上。鉴于超导谐振器的泄露时间约为50微秒，而贝尔态制备时间在50~100纳秒范围内，所以即使考虑谐振器泄露率的影响，我们的贝尔态保真度仍将高于99%，结果较好，理论上可以基于此方案进一步构建量子比特和量子门。我们的工作为分布式超导量子计算提供了一定理论支持。

参考文献