基于超绝热捷径技术快速制备超导三量子比特Greenberger-Horne-Zeilinger态*

于宛让 计新

(延边大学理学院物理系, 延吉 133002)

(2018 年 10 月 29 日收到; 2018 年 12 月 3 日收到修改稿)

本文提出了一个基于超绝热捷径技术快速制备超导三量子比特Greenberger-Horne-Zeilinger态的理论方案. 该方案首先在量子Zeno动力学的帮助下得到系统的有效哈密顿量, 之后通过引入与有效哈密顿量具有相同形式的反向导热哈密顿量来构建绝热捷径, 加速了整个系统的演化过程. 该方案不需要初态和目标态之间的直接耦合, 在实验上也更容易实现. 数值模拟结果表明该方案对超导量子比特的自发辐射、波导谐振腔的泄漏以及超导量子比特的退相位是鲁棒的.

1 引 言

随着量子信息学的发展, 纠缠态已经应用在许多前沿领域, 如量子隐形传送、量子密钥分配、量子安全直接通信等[1−3]. 其中一种著名的纠缠态就是 Greenberger-Horne-Zeilinger (GHZ)态, 它提供了一种不使用贝尔不等式就可以验证量子力学违背局域隐变量理论的可能性[4,5]. 迄今为止, 人们在不同的系统中制备了GHZ态, 这些系统包括离子系统[6]、光子系统[7]、原子系统[8]、超导系统[9]等.相比较而言, 超导系统的优势在于它的可操控性,即通过调节电磁信号以及能级结构可以让超导量子比特之间的耦合强度达到理想数值[10−12].

近几年来, 在超导系统中人们应用了多种方法制备 GHZ 态. 2016 年, Wu 等[13]利用受激拉曼绝热过程的方法制备了n量子比特GHZ态, 但是,受激拉曼绝热过程的方法需要较长的演化时间来抑制非绝热过程, 而长的演化时间势必会对整个系统的消相干造成较大的影响. 为了解决绝热过程演化缓慢的问题, 2017年Zhang等和Wu等分别用无跃迁量子驱动的方法[14]以及Lewis-Riesenfeld不变量的方法[15]在超导系统中制备了GHZ态. 然而, 这两种方法在实验上实现起来分别遇到了不同的挑战. 无跃迁量子驱动的方法需要初态和目标态之间的直接耦合, 这在实际操作中是很难实现的;在Lewis-Riesenfeld不变量的方法中需要使用短时间截取的两个简谐脉冲, 在现有的实验技术条件下是很难获得这样的脉冲的. 为了克服上述两种方案中的不足之处, 本文提出基于超绝热迭代相互作用的方法来制备GHZ态. 超绝热迭代相互作用作为传统绝热近似的延伸第一次由Berry[16]介绍, 这里通过两次迭代相互作用获得超绝热态作为演化路径, 构建的反向导热哈密顿量与利用量子Zeno动力学化简后的哈密顿量具有相同的形式.与之前的两种方案相比, 该方案的优点在于不需要初态和目标态之间的直接耦合, 并且该方案中所使用的驱动脉冲可以在曲线拟合的帮助下用两个高斯型脉冲所替代, 这两个优点极大地增加了该方案在实验上实施的可行性. 数值模拟结果显示制备的纠缠态保真度较高, 并且对于波导谐振腔的耗散、超导量子比特的自发辐射以及超导量子比特的退相位是鲁棒的.

2 物理模型和有效动力学

如图1所示, 三个超导量子比特与两个共面的波导谐振腔(CPWRs)通过电容器相互耦合,三个超导量子比特具有相同的能级结构, 都具有一个激发态|e〉j, 两个基态|L〉j和|R〉j(j=1,2,3 ).对于超导量子比特SQ1和超导量子比特SQ3,的跃迁分别由拉比频率为的经典激光驱动.的跃 迁 与 CPWR1共 振 耦 合, 耦 合 强 度 为λL,的跃迁与CPWR2共振耦合, 耦合强度为λR. 在相互作用绘景下, 系统的总哈密顿量为 ( ℏ =1 )

图1 制备超导三量子比特 GHZ 态的装置图Fig.1. Setup for generating GHZ state of three superconducting qubits.

根据量子Zeno动力学理论, 得到系统的有效哈密顿量为

3 基于超绝热捷径方法快速制备超导三量子比特GHZ态

Heff(t)的瞬时本征态也称作绝热本征态, 其形式 如下:

4 数值模拟

式中Ω0是高斯脉冲的振幅,tf是操作时间, 选定参数通过这两个脉冲, 可以确定θ0(t) 和θ1(t) . 图 2 给出了θ0(t) 和θ1(t) 随时间变化的图像, 可以看出对于任意的Ω0, 均满足边界条件但 是Ω0数 值 的 大 小对边界条件的影响很大. 为了找到合适大小的Ω0, 图3绘制了保真度随Ω0变化的图像, 可以看到当Ω0的取值在10—50之内时,均可获得相对较高的目标态保真度, 但是考虑到量子 Zeno 条件又不能过大, 因为一旦不满足量子Zeno条件, 系统在演化的过程中会经历更多的耗散. 因此, 为了保证该方案的鲁棒性,在接下来的讨论中取为了进一步说明该方案的可行性, 在图4中根据脉冲绘 制 出七 个 状 态 的 布居转移, 可以看出在时可以近乎完美地获得GHZ态, 并且不包含在GHZ态中的其他状态在整个演化过程中布居占有数很小. 因此, 本文提出的方案是有效的.

图2 θ 0(t) 和 θ 1(t) 随时间的变化情况Fig.2. Dependence on t of θ 0(t) and θ 1(t) .

接下来, 考虑超导量子比特的自发辐射、波导谐振腔的泄漏以及超导量子比特的退相位对整个系统的影响. 当考虑这些消相干的影响时, 整个系统的主方程可以表示为

图3 保真度随 Ω0 的变化情况Fig.3. Dependence on Ω 0 of fidelity ( λ =100tf-1 ).

图4 量子态布局随时间的变化Fig.4. Population versus

图5(a)—(c)分别给出了保真度随γ和κ的变化情况, 保真度随γ和γϕ的变化情况, 以及保真度随κ和γϕ的变化情况. 从图 5 可以看出, 超导量子比特的自发辐射对系统的保真度几乎没有影响, 这主要是因为如图4 所展现的, 在该方案中, 激发态在整个演化过程中的布居数几乎为0; 与自发辐射相比, 波导谐振腔的泄漏对系统保真度的影响稍大, 但从图5(a)和图5(c)可以看出, 当κ从 0 变化到时, 系统的保真度变化很小, 所以该方案对波导谐振腔的泄漏也是鲁棒的; 对比图5(a)、图5(b)、图5(c)三个图可以看出,超导量子比特的退相位对系统的保真度影响最大,但即使当波导谐振腔的泄漏率且超导量子比特的退相位率时, 目标态保真度仍高达99.76%, 而在前面提到的通过无跃迁量子驱动的方法[14]和通过Lewis-Riesenfeld不变量的方法[15]制备 GHZ态中, 当时, 保真度分别为 0.922 和 0.87, 所以可以看出本文方案与之前的方案相比对于超导量子比特的退相位是更加鲁棒的.

图5 (a) γ 和 κ 对目标态保真度的影响; (b) γ ϕ 和 γ 对目标态保真度的影响; (c) κ 和 γ ϕ 对目标态保真度的影响Fig.5. (a) Denpendence on γ and κ of fidelity; (b) denpendence on γ and γϕ of fidelity; (c) dependence on κ and γϕ of fidelity.

5 结 论

在量子Zeno动力学的帮助下, 提出了一个实验上可行的理论方案, 通过超绝热捷径技术快速制备超导三量子比特GHZ态. 与无跃迁量子驱动的方法和Lewis-Riesenfeld不变量的方法相比, 该方案的优点在于以超绝热态作为演化路径, 构建的反向导热哈密顿量与利用量子Zeno动力学化简后的哈密顿量具有相同的形式, 不需要初态与目标态之间的直接耦合, 并且所使用的驱动脉冲在曲线拟合的帮助下在实验上可以用两个高斯型脉冲所替代,这些优点极大地保证了该方案在实验上实施的可行性. 数值模拟结果表明该方案对于超导量子比特的能量弛豫、波导谐振腔的泄漏以及超导量子比特的退相位是鲁棒的.