超导非绝热几何量子计算及其优化控制

徐靖; 薛正远

doi:10.6054/j.jscnun.2020087

超导非绝热几何量子计算及其优化控制

徐靖,
薛正远^,

华南师范大学物理与电信工程学院，广州 510006

基金项目:

国家自然科学基金项目 11874156

详细信息

通讯作者:
薛正远，教授，Email:zyxue@scnu.edu.cn

中图分类号: O413
计量
- 文章访问数: 572
- HTML全文浏览量: 188
- PDF下载量: 92
出版历程
- 收稿日期: 2020-06-04
- 网络出版日期: 2021-01-04
- 刊出日期: 2020-12-24

Superconducting Nonadiabatic Geometric Quantum Computation with Optimal Control

XU Jing,
XUE Zhengyuan^,

School of Physics and Telecommunication Engineering, South China Normal University, Guangzhou 510006, China

摘要

摘要: 为实现量子门的高保真度和强鲁棒性，提出基于超导量子电路体系的非绝热几何量子计算方案.仅通过对超导比特施加含时共振微波驱动的方式，可以在超导比特上实现任意的单比特几何量子门.同时，在2个电容耦合的超导比特体系中，非平庸的2比特几何量子门也可以类似地实现.结果表明：提出的非绝热几何量子计算方案不仅对几何量子操作具有较好的鲁棒性，还可以与优化控制技术兼容，进一步增强量子门的鲁棒性.该方案的提出使容错固态量子计算的研究与发展向前迈出了重要的一步.
- 量子计算 /
- 超导量子电路 /
- 非绝热演化 /
- 几何相位 /
- 优化控制
Abstract: In order to realize quantum gates with high fidelity and strong robustness, a scheme of nonadiabatic geometric quantum computation based on superconducting quantum circuits is proposed. Arbitrary single-qubit geometric quantum gates can be realized by applying a time-dependent resonant microwave field to a superconducting qubit. Meanwhile, nontrivial two-qubit geometric quantum gates can be realized similarly on two capacitively coupled qubits. The results show that the proposed scheme not only have good robustness of geometric operations but also are compatible with optimal control technology to further enhance the gate robustness. The study makes an important step towards fault-tolerant solid-state quantum computation.
- quantum computing /
- superconducting circuit /
- nonadiabatic evolution /
- geometric phase /
- optimal control

HTML全文

量子计算具有巨大的应用潜力，因此量子计算已成为当前最受瞩目的新兴科技.但是，量子计算的实验实现研究仍然面临着诸多挑战，其中最显著的因素包括量子系统与外界环境相互作用导致退相干，以及量子操控时存在系统误差等问题.因此，针对这些问题，快速高鲁棒性的量子逻辑操作就显得尤为重要.因为几何相位^[1-3]具有内禀的抗操作误差特性^[4-5]，所以它可被用于实现高鲁棒性的量子门.

在早期的研究中，几何量子计算通过绝热演化的方式实现.此时，由于绝热演化过程要求的操作时间很长，所以在此过程中，量子态受退相干的影响非常严重.为了克服这个问题，相关研究报道了基于非绝热阿贝尔^[6-9]和非阿贝尔几何相位^[10-11]的量子计算策略.然而，这些实现方案对系统的控制误差与动力学门一样敏感，因而几何量子门的鲁棒性优势被减弱^[12-13].为了解决这个问题，基于非阿贝尔几何相位，有研究从理论上^[14-15]提出并通过实验^[16]实现了可与优化控制技术^[17]相结合的量子计算方案.

基于非阿贝尔几何相位量子门方案的实验实现研究比较复杂，而基于二能级体系的阿贝尔几何量子门的实验更易实现.另外，由于超导电路系统^[18]易实现大规模集成且具有非常好的可操控性，因而该系统受到了广泛关注.所以，本文基于超导电路体系提出一种在实验上切实可行的非绝热几何量子计算方案.该方案可以很好地与优化控制技术相结合，可进一步提高量子门的鲁棒性.因此，本文的工作向着容错固态量子计算方向的发展迈出了重要的一步.

1. 研究方法

利用动力学不变量加速绝热演化过程，是避免因演化时间过长导致退相干对量子门影响过大的有效手段.在计算基矢S₁={|0〉, |1〉}中，对于哈密顿量为

$H\left( t \right) = \frac{1}{2}\left( {\begin{array}{*{20}{c}} 0&{\mathit{\Omega }(t){{\rm{e}}^{{\rm{ - i}}\xi \left( t \right)}}}\\ {\mathit{\Omega }(t){{\rm{e}}^{{\rm{i}}\xi \left( t \right)}}}&0 \end{array}} \right)$

(1)

的二能级系统，可以找到满足冯诺依曼关系式

$\frac{\partial }{{\partial t}}\mathit{\boldsymbol{I}} + {\rm{i}}\left[ {\mathit{\boldsymbol{H}}, \mathit{\boldsymbol{I}}} \right] = 0$

的对应不变量I.不变量的矩阵形式可写为

$\mathit{\boldsymbol{I}}\left( t \right) = \frac{{{\mathit{\Omega }_0}}}{2}\left( {\begin{array}{*{20}{c}} {{\rm{cos}}\chi \left( t \right)}&{{\rm{sin}}\chi (t){{\rm{e}}^{ - {\rm{i}}\beta \left( t \right)}}}\\ {{\rm{sin}}\chi (t){{\rm{e}}^{{\rm{i}}\beta \left( t \right)}}}&{ - {\rm{cos}}\chi \left( t \right)} \end{array}} \right){\rm{, }}$

(2)

令h=1，Ω₀为具有频率单位的常数.不变量的本征态为

$|\left\langle {{\mathit{\boldsymbol{ \mathit{\boldsymbol{\phi}} }}_ + }\left( t \right)} \right. = \cos \frac{\chi }{2}{{\rm{e}}^{ - {\rm{i}}\beta \left( t \right)}}|\left. 0 \right\rangle + \sin \frac{\chi }{2}{{\rm{e}}^{{\rm{i}}\beta /2}}|\left. 1 \right\rangle ,$

$|\left\langle {{\mathit{\boldsymbol{ \mathit{\boldsymbol{\phi}} }}_ - }\left( t \right)} \right. = - \sin \frac{\chi }{2}{{\rm{e}}^{{\rm{ - i}}\beta /2}}|\left. 0 \right\rangle + \cos \frac{\chi }{2}{{\rm{e}}^{{\rm{i}}\beta /2}}|\left. 1 \right\rangle ,$

其中, 含时变化的角度χ(t)和β(t)可作为反向设计H(t)的辅助参量.

根据不变量理论，可以找到H(t)相应薛定谔方程的2个正交完备解，即|ψ₊(t)〉=e^if₊(t)|ϕ₊(t)〉、|ψ_-(t)〉=e^if_-(t)|ϕ₊(t)〉，其中f₊和f_－为相应的相位，且有f₊=－f_－.为了更方便地得到1个纯几何相位，本文考虑哈密顿量驱动演化初态并进行循环演化的情况.设定在演化初始时刻满足f₊(0)=f_－(0)=0，此时薛定谔方程的2个解与不变量本征态重合，从而可以利用不变量的本征态作为获得几何相位的演化态.由于|ϕ₊(t)〉和|ϕ_-(t)〉的演化路径具有对称性，所以仅需设计其中一个态矢量(如|ϕ₊(t)〉的演化路径)即可反向确定微波驱动场的波形Ω(t)和相位ξ(t).在循环演化的终止时刻，令γ=f₊(τ)，可以将哈密顿量对应的演化算符写为

$\mathit{\boldsymbol{U}}\left( \tau \right) = {{\rm{e}}^{{\rm{i \mathit{ γ} }}}}{\rm{|}}{\mathit{\boldsymbol{\phi}} _{\rm{ + }}}\left. {\left( 0 \right)} \right\rangle \left\langle {{\mathit{\boldsymbol{\phi}} _{\rm{ + }}}\left( 0 \right)} \right.| + {{\rm{e}}^{{\rm{ - i \mathit{ γ} }}}}|{\mathit{\boldsymbol{\phi}} _ - }\left. {\left( 0 \right)} \right\rangle \left\langle {{\mathit{\boldsymbol{\phi}} _ - }\left( 0 \right)} \right.|,$

(3)

其中，相位γ包含动力学部分

$\begin{array}{l} {{\rm{ \mathit{ γ} }}_d} = - \int_0^\tau {\left\langle {\mathit{\boldsymbol{\phi}} \left( t \right)|\mathit{\boldsymbol{H}}\left( t \right)} \right.} \left. {|\mathit{\boldsymbol{\phi}} \left( t \right)} \right\rangle {\rm{d}}t= \\ \;\;\;\;\;\;\;\frac{1}{2}\int_0^\tau {\frac{{\dot \beta \left( t \right){{\sin }^2}\chi \left( t \right)}}{{\cos \chi \left( t \right)}} {\rm{d}}t} \end{array}$

(4)

和几何部分

${{\rm{ \mathit{ γ} }}_g} = {\rm{i}}\int_0^\tau {\left\langle {{\mathit{\boldsymbol{\phi}}} \left( t \right)} \right.} \left. {|\dot{\mathit{\boldsymbol{\phi}}} \left( t \right)} \right\rangle {\rm{d}}\mathit{t} = \frac{1}{2}\int_0^\tau {\dot \beta \left( t \right)\cos \chi \left( t \right){\rm{d}}\mathit{t}} .$

(5)

当消除动力学部分后，即可得到1个纯几何的相位.最后，将初始时刻的不变量本征态代入演化算符表达式中可得：

$\mathit{\boldsymbol{U}}\left( \tau \right) = \left( \begin{array}{l} \cos \gamma + {\rm{i}}\cos {\chi _0}\sin \gamma \;\;{\rm{i}}\sin \gamma \sin {\chi _0}{{\rm{e}}^{ - {\rm{i}}{\beta _0}}}\\ {\rm{i}}\sin \gamma \sin {\chi _0}{{\rm{e}}^{{\rm{i}}{\beta _0}}}\;\;\cos \gamma - {\rm{i}}\cos {\chi _0}\sin \gamma \end{array} \right) = {{\rm{e}}^{{\rm{i}}\gamma \mathit{\boldsymbol{\vec n}} \cdot \mathit{\boldsymbol{\vec \sigma }}}},$

(6)

其中, ${\chi _0} = \chi \left( 0 \right), {\beta _0} = \beta \left( 0 \right), \mathit{\boldsymbol{\vec n = }}\left( {\sin {\chi _0}\cos {\beta _0}, \sin {\chi _0} \times \sin {\beta _0}, \cos {\chi _0}} \right)$ .式(6)即绕任意轴 $\mathit{\boldsymbol{\vec n}}$ 旋转-2γ角度的几何操作，可以用来实现任意的单比特量子门.

另外，在2个电容耦合的超导比特体系中^[9]，通过对其中一个超导比特频率的含时调制，使2个比特之间的耦合达到等效共振.此时，在{|11〉_AB, |20〉_AB}子空间共振的情况下，2比特的相互作用形式与式(1)的形式完全相同，因此可以类比单比特的方式实现非平庸的2比特几何量子门.

2. 几何量子门的实现及其优化控制

根据上述研究方法，首先，利用非绝热循环演化的方式产生阿贝尔几何相位，从而在超导transmon比特^[19]中实现任意的单比特几何量子门; 然后给出绕X和Z轴旋转这2种门操作的具体实现方案(图 1)，并进行数值模拟; 最后采用优化控制技术进一步提升量子门的鲁棒性.

图 1 单量子比特门

Figure 1. The single-qubit gate

下载: 全尺寸图片幻灯片

2.1 量子门的构建

利用微波共振驱动超导transmon量子比特最低的2个能级可以实现单比特量子门(图 1A).但此时驱动微波场也会在超导比特高能态之间引起不必要的色散相互作用.这里先作简化处理，只考虑最低能态的2个能级，则系统的哈密顿量可写为式(1)的形式.根据不变量的冯诺依曼公式，可得到辅助参数χ(t)、β(t)与驱动频率Ω(t)、相位ξ(t)间的关系式：

$\dot \chi \left( t \right) = - \mathit{\Omega }\left( t \right){\rm{sin}}\left[ {\beta \left( t \right) + \xi \left( t \right)} \right]{\rm{, }}$

(7)

$\dot \beta \left( t \right) = - \mathit{\Omega }\left( t \right){\rm{cos}}\chi \left( t \right){\rm{cos}}[\beta \left( t \right) + \xi \left( t \right)].$

(8)

由此可以通过设计循环演化过程中的辅助参数χ(t)和β(t), 进行数值计算, 求出Ω(t)和ξ(t).同时为了简化表述，引入一个含时参数f(t)=－2f₊(t).联立(3)、(4)可以得到f(t)的导数函数:

$\dot f\left( t \right) = - \frac{{\dot \beta \left( t \right)}}{{{\rm{cos}}\chi \left( t \right)}}.$

(9)

在上式的基础上，只需要选定χ(t)和f(t)的函数形式，β(t)的形式就可以确定.接下来将具体实现绕X轴旋转的几何门和绕Z轴旋转的几何门.

2.1.1 绕X轴旋转操作的几何门

为实现绕X轴旋转的几何门，将循环过程分为4个相等的时间段，态矢量|ϕ₊(t)〉的极角和方位角函数可设计为

$\begin{array}{l} {\chi _1}\left( t \right) = {\rm{ \mathit{ π} }}\left[ {1 + {\rm{si}}{{\rm{n}}^2}(2{\rm{ \mathit{ π} }}t/\tau )} \right]/2, \\ \;\;\;\;{\beta _1}\left( 0 \right) = 0, t \in \left[ {0, \tau /4} \right]; \end{array}$

(10)

$\begin{array}{l} \;\;\;\;\;\;\;{\chi _2}\left( t \right) = {\rm{ \mathit{ π} }}\left[ {1 + {\rm{si}}{{\rm{n}}^2}(2{\rm{ \mathit{ π} }}t/\tau )} \right]/2, \\ {\beta _2}\left( {\tau /4} \right) = {\beta _1}\left( {\tau /4} \right) - \gamma , t \in \left[ {\tau /4, \tau /2} \right]; \end{array}$

(11)

$\begin{array}{l} \;\;\;\;\;\;\;{\chi _3}\left( t \right) = {\rm{ \mathit{ π} }}\left[ {1 - {\rm{si}}{{\rm{n}}^2}(2{\rm{ \mathit{ π} }}t/\tau )} \right]/2, \\ {\beta _3}\left( {\tau /2} \right) = {\beta _2}\left( {\tau /2} \right), t \in \left[ {\tau /2, 3\tau /4} \right]; \end{array}$

(12)

$\begin{array}{l} \;\;\;\;\;\;\;{\chi _4}\left( t \right) = {\rm{ \mathit{ π} }}\left[ {1 - {\rm{si}}{{\rm{n}}^2}(2{\rm{ \mathit{ π} }}t/\tau )} \right]/2, \\ {\beta _4}\left( {3\tau /4} \right) = {\beta _3}\left( {3\tau /4} \right) + \gamma , t \in \left[ {3\tau /4, \tau } \right]. \end{array}$

(13)

选取f_j(t)=cos2χ_j(t)/5，每段方位角用

${\beta _j}\left( t \right) = - \int {{{\dot f}_j}\left( t \right){\rm{cos}}{\chi _j}\left( t \right){\rm{d}}t}$

来计算.为了验证循环演化产生的相位为纯几何相位，首先计算演化过程中的动力学相位

${\gamma _d} = \sum\limits_{j = 1}^4 {\int_{\left( {j - 1} \right)\tau /4}^{j\tau /4} { - {{\dot f}_j}\left( t \right){\rm{si}}{{\rm{n}}^2}{\chi _j}\left( t \right)/2{\rm{d}}t = 0} } ,$

(14)

同时，由于方位角β(t)在t=τ /4和t=3τ/4这2个时刻发生突变，从而容易证明几何相位γ_g=γ.这样就实现了绕X轴旋转操作的纯几何门.

2.1.2 绕Z轴旋转操作的几何门

为了实现绕Z轴旋转的几何操作，将单圈演化成2段.这2段路径上的极角和方位角分别为

$\begin{array}{l} {\chi _1}\left( t \right) = \left. {{\rm{ \mathit{ π} si}}{{\rm{n}}^2}\left( {{\rm{ \mathit{ π} }}t/\tau } \right)} \right]\\ {\beta _1}\left( 0 \right) = 0, t \in \left[ {0, \tau /2} \right]{\rm{;}} \end{array}$

(15)

$\begin{array}{l} {\chi _2}\left( t \right) = \left. {{\rm{ \mathit{ π} si}}{{\rm{n}}^2}\left( {{\rm{ \mathit{ π} }}t/\tau } \right)} \right]\\ {\beta _2}\left( {\tau /2} \right) = {\beta _1}\left( {\tau /2} \right) - \gamma , \\ \;\;\;\;\;\;\;t \in \left[ {\tau /2, \tau } \right]. \end{array}$

(16)

选取f_j(t)=(2χ_j(t)－sin2χ_j(t))/5，和实现绕X轴旋转操作一样，每段方位角也用

${\beta _j}\left( t \right) = - \int {{{\dot f}_j}\left( t \right){\rm{cos}}{\chi _j}\left( t \right){\rm{d}}t}$

来计算.同样可计算演化过程的动力学相位

${\gamma _d} = \sum\limits_{j = 1}^2 {\int_{\left( {j - 1} \right)\tau /2}^{j\tau /2} { - {{\dot f}_j}\left( t \right){\rm{si}}{{\rm{n}}^2}{\chi _j}\left( t \right)/2{\rm{d}}t = 0} {\rm{. }}}$

(17)

2.2 量子门的性能

在量子门的物理实现过程中，退相干是影响保真度的一个重要因素.由于超导transmon量子比特的非谐性较小，本文采取有效的修正方案^[20]来抑制因量子信息泄露导致的误差.通常使用Lindblad主方程评估量子门的性能.主方程为

${\mathit{\boldsymbol{\dot \rho }}_1} = {\rm{i}}\left[ {{\mathit{\boldsymbol{\rho }}_1}, \mathit{\boldsymbol{H}} + {\mathit{\boldsymbol{H}}_{{\rm{leak}}}}} \right] + \left[ {{\mathit{\boldsymbol{ \boldsymbol{\varGamma} }}_1}\mathit{\boldsymbol{L}}\left( {{\mathit{\boldsymbol{\sigma }}_1}} \right) + {\mathit{\boldsymbol{ \boldsymbol{\varGamma} }}_2}\mathit{\boldsymbol{L}}\left( {{\mathit{\boldsymbol{\sigma }}_2}} \right)} \right],$

(18)

且有

${\mathit{\boldsymbol{H}}_{{\rm{leak}}}}\left( t \right) = - \alpha \left. {|2} \right\rangle \left\langle {2|} \right. + \left[ {\sqrt 2 \mathit{\Omega }(t){\rm{e}^{{\rm{i}}\xi \left( t \right)}}\left. {|1} \right\rangle \left\langle {2| + {\rm{H}}.{\rm{c}}.} \right.} \right],$

(19)

其中，L(A)=2Aρ₁A⁺－A⁺Aρ₁－ρ₁A⁺A, 是算符A的Lindblad算符， ${\mathit{\boldsymbol{\sigma }}_1} = |\left. 0 \right\rangle \left\langle {1|} \right. + \sqrt 2 \left. {|1} \right\rangle$ 和σ₂=|1〉〈1|+2|2〉〈2|. Γ₁和Γ₂分别代表transmon的衰减和退相速率.根据现有的实验技术条件^[18]，设Γ₁=Γ₁=2π×2 kHz，transmon比特的非谐为α=2π×300 MHz，微波驱动波形振幅最大值Ω_max=2π×16 MHz.这里通过选取χ₀=π/2、β₀=0、γ=π/2来实现非门，通过选取χ₀=0、β₀=0、γ=-π/8来实现相位门.二者循环演化周期分别为102 ns和125 ns.实现这2个门所需的驱动振幅Ω(t)波形和相位(t)波形如图 2A和图 2B所示. 1个量子态|Φ(t)〉经过单比特量子门演化以后的保真度F=〈 Φ (τ)|ρ₁|Φ(τ)〉.对于非门，假设演化初态| Φ(0)〉_N=|0〉，则其演化末态|Φ(τ)〉_N=|1〉.对于相位门，假设演化初态 $\mathit{\boldsymbol{| \boldsymbol{\varPhi} }}{\left. {\left( 0 \right)} \right\rangle _T} = \left( {|\left. 0 \right\rangle + |\left. 1 \right\rangle } \right)/\sqrt 2$ ，则其演化末态 $\mathit{\boldsymbol{| \boldsymbol{\varPhi} }}{\left. {\left( \tau \right)} \right\rangle _T} = \left( {|\left. 0 \right\rangle + {\rm{e}^{{\rm{i \mathit{ π} }}/4}}|\left. 1 \right\rangle } \right)/\sqrt 2$ .根据数值计算结果，非门和相位门的态保真度分别为F_N=99.87%和F_T=99.80%(图 2C、D).为了计算量子门的保真度，这里假设初态均为|Φ(0)〉=cosθ|0〉+sinθ|1〉，则演化以后，非门得到的末态为|Φ(τ)〉_N=sinθ|0〉+ cosθ|1〉，相位门得到的末态为|Φ(τ)〉_T=cosθ|0〉+e^iπ/4sinθ|1〉. θ在[0, 2π]间均匀取1 001个值，通过定义量子门的保真度

图 2 单量子比特几何门的实现及其性能

Figure 2. The implementation of geometric single-qubit gates and their performance

下载: 全尺寸图片幻灯片

${F^G} = 1/2{\rm{ \mathit{ π} }}\int_0^{2{\rm{ \mathit{ π} }}} {\left\langle {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}\left( \tau \right)} \right.\left| {{\rho _1}} \right|\left. {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}\left( \tau \right)} \right\rangle } {\rm{d}}\theta ,$

最终计算得到非门和相位门的保真度分别为F_N^G=99.87%和F_T^G=99.84%.

2.3 量子门的优化控制

基于前面实现的单比特量子门，接下来将优化控制技术^[17]应用到演化路径的设计上，以实现绕Z轴旋转的几何操作.假设存在一个静态的系统误差，即V(t)=εH(t).由于演化路径由对称的2段组成，所以只需要优化第一段演化路径.这里用|ψ₊(t)〉来表示演化过程，则演化初态为|ψ₊(0)〉，理想的目标态为|ψ₊(τ/2)〉.设目标函数

$P = |\left\langle {{\mathit{\boldsymbol{\psi }}_ + }} \right.\left( {\tau /2} \right)|\left. {\mathit{\boldsymbol{ \boldsymbol{\varPsi} }}\left( {\tau /2} \right)} \right\rangle {|^2} = 1 + {\tilde O_1} + {\tilde O_2} + \cdots ,$

(20)

其中，|Ψ(t)〉为存在系统误差时系统的量子态， ${\tilde O_n}$ 表示n阶微扰.这里忽略高阶项，只保留到二阶，即 $P = 1 + {\tilde O_1} + {\tilde O_2}.$ .通过计算可得

${\tilde O_1} = 0, {\tilde O_2} = - |\int_0^{\tau /2} {z\left( t \right){\rm{d}}t} {|^2}.$

(21)

其中,

$z\left( t \right) = \frac{\varepsilon }{2}(\frac{1}{2}{\rm{ }}\dot f(t){\rm{sin}}\left( {2\chi \left( t \right)} \right) - {\rm{i}}\dot \chi \left( t \right)){e^{{\rm{i}}f(t)}}.$

为了消除了二阶微扰的影响，必须要求z(t)在t=[0, τ/2]时间段积分为0.对应方程解的形式可以写为

$f\left( \chi \right) = 2\chi + {C_1}\sin \left( {2\chi } \right) + {C_2}\sin \left( {4\chi } \right) + \cdots ,$

(22)

通过数值计算，当取C₁=-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)相位门的保真度

Figure 3. The fidelity of phase gate without (η=0) and with (η=1) optimization

下载: 全尺寸图片幻灯片

3. 结论

基于超导量子电路，提出了构建快速通用的几何量子计算方案.通过施加共振耦合驱动，在单个超导比特上实现了任意的单量子比特门.同时提出的方案可以结合优化技术，提高量子门对静态系统误差的鲁棒性.在2比特超导电容耦合体系中可实现非平庸的双比特量子门.由于本文的方案可以直接扩展到二维比特阵列情形，因此研究结果为实现容错量子计算提供了一条可靠的道路.

图 1 单量子比特门

Figure 1. The single-qubit gate

下载: 全尺寸图片幻灯片

图 2 单量子比特几何门的实现及其性能

Figure 2. The implementation of geometric single-qubit gates and their performance

下载: 全尺寸图片幻灯片

图 3 未优化(η=0)与优化(η=1)相位门的保真度

Figure 3. The fidelity of phase gate without (η=0) and with (η=1) optimization

下载: 全尺寸图片幻灯片

参考文献(20)

[1]	BERRY M V. Quantal phase factors accompanying adiabatic changes[J]. Proceedings of the Royal Society A, 1984, 392:45-47.
[2]	WILCZEK F, ZEE A. Appearance of gauge structure in simple dynamical systems[J]. Physical Review Letters, 1984, 52(24):2111-2114.
[3]	AHARONOV Y, ANANDAN J. Phase change during a cyclic quantum evolution[J]. Physical Review Letters, 1987, 58(16):1593-1596.
[4]	ZHU S L, ZANARDI P. Geometric quantum gates that are robust against stochastic control errors[J]. Physical Review A, 2005, 72(2):020301/1-4. http://scitation.aip.org/getabs/servlet/GetabsServlet?prog=normal&id=PLRAAN000072000002020301000001&idtype=cvips&gifs=Yes
[5]	JOHANSSON M, ANDERSSON L M, ERICSSON M, et al. Robustness of nonadiabatic holonomic gates[J]. Physical Review A, 2012, 86(6):062322/1-10.
[6]	WANG X B, KEIJI M. Nonadiabatic conditional geometric phase shift with NMR[J]. Physical Review Letters, 2001, 87(9):097901/1-4. http://europepmc.org/abstract/med/11531598
[7]	ZHU S L, WANG Z D. Implementation of universal quantum gates based on nonadiabatic geometric phases[J]. Physical Review Letters, 2002, 89(9):097902/1-4. http://scitation.aip.org/getabs/servlet/GetabsServlet?prog=normal&id=VIRT01000006000008000071000001&idtype=cvips&gifs=Yes
[8]	ZHAO P Z, CUI X D, XU G F, et al. Rydberg-atom-based scheme of nonadiabatic geometric quantum computation[J]. Physical Review A, 2017, 96(5):052316/1-6. doi: 10.1103/PhysRevA.96.052316
[9]	CHEN T, XUE Z Y. Nonadiabatic geometric quantum computation with parametrically tunable coupling[J]. Physical Review Applied, 2018, 10(5):054051/1-13.
[10]	SJÖVIST E, TONG D M, ANDERSSON L M, et al. Non-adiabatic holonomic quantum computation[J]. New Journal of Physics, 2012, 14(10):103035/1-10.
[11]	XU G F, ZHANG J, TONG D M, et al. Nonadiabatic holonomic quantum computation in decoherence-free subspaces[J]. Physical Review Letters, 2012, 109(17):170501/1-5. http://www.ncbi.nlm.nih.gov/pubmed/23215167
[12]	ZHENG S B, YANG C P, NORI F. Comparison of the sensitivity to systematic errors between nonadiabatic non-Abelian geometric gates and their dynamical counterparts[J]. Physical Review A, 2016, 93(3):032313/1-5.
[13]	JING J, LAM C H, WU L A. Non-Abelian holonomic transformation in the presence of classical noise[J]. Phy-sical Review A, 2017, 95(1):012334/1-7.
[14]	LIU B J, SONG X K, XUE Z Y, et al. Plug-and-play approach to nonadiabatic geometric quantum gates[J]. Phy-sical Review Letters, 2019, 123(10):100501/1-6. http://www.ncbi.nlm.nih.gov/pubmed/31573289
[15]	LI S, CHEN T, XUE Z Y. Fast holonomic quantum computation on superconducting circuits with optimal control[J]. Advanced Quantum Technologies, 2020, 3(3):2000001/1-7.
[16]	YAN T, LIU B J, XU K, et al. Experimental realization of nonadiabatic shortcut to non-Abelian geometric gates[J]. Physical Review Letters, 2019, 122(8):080501/1-6.
[17]	DAEMS D, RUSCHHAUPT A, SUGNY D, et al. Robust quantum control by a single-shot shaped pulse[J]. Phy-sical Review Letters, 2013, 111(5):050404/1-5. http://europepmc.org/abstract/med/23952372
[18]	DEVORET M H, SCHOELKOPF R J. Superconducting circuits for quantum information:an outlook[J]. Science, 2013, 339:1169-1174.
[19]	KOCH J, YU T M, MAJER J, et al. Charge-insensitive qubit design derived from the Cooper pair box[J]. Physical Review A, 2007, 76(4):042319/1-19. http://scitation.aip.org/getabs/servlet/GetabsServlet?prog=normal&id=VIRT04000007000010000140000001&idtype=cvips&gifs=Yes
[20]	WANG T H, ZHANG Z X, XIANG L, et al. The experimental realization of high-fidelity ' shortcut-to-adiabati- city' quantum gates in a superconducting Xmon qubit[J]. New Journal of Physics, 2018, 20(6):065003/1-10.

施引文献

资源附件(0)

图(3)

计量

文章访问数: 572
HTML全文浏览量: 188
PDF下载量: 92
被引次数: 0

1. 研究方法
2. 几何量子门的实现及其优化控制
2.1 量子门的构建
2.1.1 绕X轴旋转操作的几何门
2.1.2 绕Z轴旋转操作的几何门
2.2 量子门的性能
2.3 量子门的优化控制
3. 结论

超导非绝热几何量子计算及其优化控制

通讯作者: 薛正远，教授，Email:zyxue@scnu.edu.cn

计量

出版历程