徐 靖,薛正远
(华南师范大学物理与电信工程学院,广州 510006)
量子计算具有巨大的应用潜力,因此量子计算已成为当前最受瞩目的新兴科技. 但是,量子计算的实验实现研究仍然面临着诸多挑战,其中最显著的因素包括量子系统与外界环境相互作用导致退相干,以及量子操控时存在系统误差等问题. 因此,针对这些问题,快速高鲁棒性的量子逻辑操作就显得尤为重要. 因为几何相位[1-3]具有内禀的抗操作误差特性[4-5],所以它可被用于实现高鲁棒性的量子门.
在早期的研究中,几何量子计算通过绝热演化的方式实现. 此时,由于绝热演化过程要求的操作时间很长,所以在此过程中,量子态受退相干的影响非常严重. 为了克服这个问题,相关研究报道了基于非绝热阿贝尔[6-9]和非阿贝尔几何相位[10-11]的量子计算策略. 然而,这些实现方案对系统的控制误差与动力学门一样敏感,因而几何量子门的鲁棒性优势被减弱[12-13]. 为了解决这个问题,基于非阿贝尔几何相位,有研究从理论上[14-15]提出并通过实验[16]实现了可与优化控制技术[17]相结合的量子计算方案.
基于非阿贝尔几何相位量子门方案的实验实现研究比较复杂,而基于二能级体系的阿贝尔几何量子门的实验更易实现. 另外,由于超导电路系统[18]易实现大规模集成且具有非常好的可操控性,因而该系统受到了广泛关注. 所以,本文基于超导电路体系提出一种在实验上切实可行的非绝热几何量子计算方案. 该方案可以很好地与优化控制技术相结合,可进一步提高量子门的鲁棒性. 因此,本文的工作向着容错固态量子计算方向的发展迈出了重要的一步.
利用动力学不变量加速绝热演化过程,是避免因演化时间过长导致退相干对量子门影响过大的有效手段. 在计算基矢S1={|0〉,|1〉}中,对于哈密顿量为
(1)
的二能级系统,可以找到满足冯诺依曼关系式
的对应不变量I. 不变量的矩阵形式可写为
(2)
令ћ=1,Ω0为具有频率单位的常数. 不变量的本征态为
其中,含时变化的角度χ(t)和β(t)可作为反向设计H(t)的辅助参量.
根据不变量理论,可以找到H(t)相应薛定谔方程的2个正交完备解,即|ψ+(t)〉=eif+(t)|φ+(t)〉、|ψ-(t)〉=eif-(t)|φ-(t)〉,其中f+和f-为相应的相位,且有f+=-f-. 为了更方便地得到1个纯几何相位,本文考虑哈密顿量驱动演化初态并进行循环演化的情况. 设定在演化初始时刻满足f+(0)=f-(0)=0,此时薛定谔方程的2个解与不变量本征态重合,从而可以利用不变量的本征态作为获得几何相位的演化态. 由于|φ+(t)〉和|φ-(t)〉的演化路径具有对称性,所以仅需设计其中一个态矢量(如|φ+(t)〉的演化路径)即可反向确定微波驱动场的波形Ω(t)和相位ξ(t). 在循环演化的终止时刻,令γ=f+(),可以将哈密顿量对应的演化算符写为
U()=eiγ|φ+(0)〉〈φ+(0)|+e-iγ|φ-(0)〉〈φ-(0)|,
(3)
其中,相位γ包含动力学部分
(4)
和几何部分
(5)
当消除动力学部分后,即可得到1个纯几何的相位. 最后,将初始时刻的不变量本征态代入演化算符表达式中可得:
U(
(6)
另外,在2个电容耦合的超导比特体系中[9],通过对其中一个超导比特频率的含时调制,使2个比特之间的耦合达到等效共振. 此时,在{|11〉AB,|20〉AB}子空间共振的情况下,2比特的相互作用形式与式(1)的形式完全相同,因此可以类比单比特的方式实现非平庸的2比特几何量子门.
根据上述研究方法,首先,利用非绝热循环演化的方式产生阿贝尔几何相位,从而在超导transmon比特[19]中实现任意的单比特几何量子门;然后给出绕X和Z轴旋转这2种门操作的具体实现方案(图1),并进行数值模拟;最后采用优化控制技术进一步提升量子门的鲁棒性.
图1 单量子比特门
利用微波共振驱动超导transmon量子比特最低的2个能级可以实现单比特量子门(图1A). 但此时驱动微波场也会在超导比特高能态之间引起不必要的色散相互作用. 这里先作简化处理,只考虑最低能态的2个能级,则系统的哈密顿量可写为式(1)的形式. 根据不变量的冯诺依曼公式,可得到辅助参数χ(t)、β(t)与驱动频率Ω(t)、相位ξ(t)间的关系式:
(7)
(8)
由此可以通过设计循环演化过程中的辅助参数χ(t)和β(t),进行数值计算,求出Ω(t)和ξ(t). 同时为了简化表述,引入一个含时参数f(t)=-2f+(t). 联立(3)、(4)可以得到f(t)的导数函数:
(9)
在上式的基础上,只需要选定χ(t)和f(t)的函数形式,β(t)的形式就可以确定. 接下来将具体实现绕X轴旋转的几何门和绕Z轴旋转的几何门.
2.1.1 绕X轴旋转操作的几何门 为实现绕X轴旋转的几何门,将循环过程分为4个相等的时间段,态矢量|φ+(t)〉的极角和方位角函数可设计为
χ1(t)=π[1+sin2(2πt/)]/2,
β1(0)=0,t[0,/4];
(10)
χ2(t)=π[1+sin2(2πt/)]/2,
β2(/4)=β1(/4)-γ,t[/4,/2];
(11)
χ3(t)=π[1-sin2(2πt/)]/2,
β3(/2)=β2(/2),t[/2,3/4];
(12)
χ4(t)=π[1-sin2(2πt/)]/2,
β4(3/4)=β3(3/4)+γ,t[3/4,].
(13)
选取fj(t)=cos 2χj(t)/5,每段方位角用
来计算. 为了验证循环演化产生的相位为纯几何相位,首先计算演化过程中的动力学相位
(14)
同时,由于方位角β(t)在t=/4和t=3/4这2个时刻发生突变,从而容易证明几何相位γg=γ. 这样就实现了绕X轴旋转操作的纯几何门.
2.1.2 绕Z轴旋转操作的几何门 为了实现绕Z轴旋转的几何操作,将单圈演化成2段. 这2段路径上的极角和方位角分别为
χ1(t)=πsin2(πt/)],
β1(0)=0,t[0,/2];
(15)
χ2(t)=πsin2(πt/)],
β2(/2)=β1(/2)-γ,
t[/2,].
(16)
选取fj(t)=(2χj(t)-sin 2χj(t))/5,和实现绕X轴旋转操作一样,每段方位角也用
来计算. 同样可计算演化过程的动力学相位
(17)
在量子门的物理实现过程中,退相干是影响保真度的一个重要因素. 由于超导transmon量子比特的非谐性较小,本文采取有效的修正方案[20]来抑制因量子信息泄露导致的误差. 通常使用Lindblad主方程评估量子门的性能. 主方程为
(18)
且有
(19)
图2 单量子比特几何门的实现及其性能
基于前面实现的单比特量子门,接下来将优化控制技术[17]应用到演化路径的设计上,以实现绕Z轴旋转的几何操作. 假设存在一个静态的系统误差,即V(t)=εH(t). 由于演化路径由对称的2段组成,所以只需要优化第一段演化路径. 这里用|ψ+(t)〉来表示演化过程,则演化初态为|ψ+(0)〉,理想的目标态为|ψ+(/2)〉. 设目标函数
(20)
(21)
其中,
为了消除了二阶微扰的影响,必须要求z(t)在t=[0,/2]时间段积分为0. 对应方程解的形式可以写为
f(χ)=2χ+C1sin(2χ)+C2sin(4χ)+…,
(22)
通过数值计算,当取C1=-1,其余系数均取0时,关于z(t)的积分为0.
为了评估本文优化方案的效果,给f(t)加系数f(t)=η[2χ(t)-sin(2χ(t))]. 结果表明:当η=0时,该方案就退回到优化前的方案[6-9];当η=1时,则对应于本文的优化方案.
在系统误差εΩmax的取值范围(2π×[-5,5]MHz)内,给出2个方案对系统误差的鲁棒性对比图(图3). 结果表明:当不考虑退相干影响时,优化方案对系统误差的鲁棒性效果要远优于未优化方案的效果;当考虑实际情况中退相干的影响时,优化方案的鲁棒性依然很出色.
图3 未优化(η=0)与优化(η=1)相位门的保真度
基于超导量子电路,提出了构建快速通用的几何量子计算方案. 通过施加共振耦合驱动,在单个超导比特上实现了任意的单量子比特门. 同时提出的方案可以结合优化技术,提高量子门对静态系统误差的鲁棒性. 在2比特超导电容耦合体系中可实现非平庸的双比特量子门. 由于本文的方案可以直接扩展到二维比特阵列情形,因此研究结果为实现容错量子计算提供了一条可靠的道路.