两个超导比特纠缠动力学的反转动波效应*钱雷1贺树 2 段立伟2陈庆虎1,2*（1.浙江师范大学,海峡两岸统计物理与凝聚态理论研究中心,金华 321004;2.浙江大学物理系,杭州 310027）摘要：本文研究了两个无相互作用的与各自量子振子耦合的超导量子比特的纠缠动力学。
数值严格结果证明在强耦合下,以往的转动波近似不再适用，非转动波效应必须考虑。
基于著名的推广的转动波近似的解析结果,在目前实验上可实现的强耦合区域与严格数值结果有明显差别。
从我们近年来发展的转动波近似的一级校正，我们得到的结果更接近于精确数值解。
该理论结果可激发基于近来可以实现超导比特的相关实验。
关键词：纠缠动力学; 量子比特; 比特-振子耦合; 反转动波.QIAN Lei1 , HE Shu2, DUAN Liwei2, CHEN Qinghu1,2*（1.Center for Statistical and Theoretical Condensed Matter Physics, Zhejiang Normal University, Jinhua 321004, China; 2.Department of Physics, Zhejiang University, Hangzhou 310027, China)Title: Effect of Counter-Rotating term on the entanglement dynamics of two superconducting qubits coupling to quantum oscillatorsAbstract:In this paper, we investigate entanglement dynamics of two non-interacting superconducting qubits coupled with their own quantum oscillators beyond the rotating-wave approximation (RWA). It is shown in the numerically exact studies that in the strong coupling regime, the RWA is not valid and the effect of the counter-rotating wave should be taken into account. The analytic results for entanglement based on the well-known generalized RWA deviates from the numerical one obviously in the present-day experimentally accessible coupling regime. In this 收稿日期： 1paper, based on our first-order corrections to the RWA, the analytic results are move close to the numerical ones. It may motivate the relevant experiments based on the superconducting qubits realized recently.Key Words:Entanglement dynamics; qubits; qubit-oscillator coupling;Counter-rotating wave.1 引言超导比特作为量子信息科学的基础和量子计算机的重要元件, 在基础研究和实际应用中受到了广泛的关注。
近年来，二能级系统与LC谐振子耦合系统的动力学研究在实验上取得重要进展，对量子计算，量子相干性和量子逻辑门等领域的发展有较大的促进。
量子纠缠是一种典型的量子关联, 描述了量子力学的本质的非局域性质, 根本性地区分了量子和经典物理学。
纠缠作为一种重要的量子信息资源, 受到越来越多的重视，而由于与环境耦合导致退相干和纠缠突然死亡（纠缠在一段时间消失）等问题是目前量子信息领域内研究的热点课题。
Jaynes-Cummings（JC）模型【1】是光学腔量子电动力学(CQED)系统中最简单的一类模型, 描述腔场中单个二能级原子和单模量子化的电磁场的相互作用问题。
在转动波近似下, 该模型可以严格解析求解。
超出转动波近似, 闭合形式的解析解至今缺乏, 各种近似的解析解相继提出, 其中最为著名的是推广的转动波近似【6】.Yu 等人研究了在转动波近似下两个JC原子的量子纠缠动力学【2，3】,并给出了相应的解析表达式。
陈等人利用数值方法研究了无转动波近似下的纠缠行为【4】, 发现反转动波对纠缠有压制效应。
至今，未见到反转动波下的纠缠动力学的解析研究或近似解析收稿日期： 2收稿日期： 3研究。
近期的超导量子比特与振子耦合系统的实验表明, 当其耦合强度达到g=0.1 （这比量子光学系统里的二能级原子和腔耦合系统的耦合强度大 3-6个数量级）， 转动波近似已被证明实效. 因此,纠缠动力学的反转动波的效应的解析研究是十分必须的和值得的.在本文中，我们基于推广转动波近似和转动波近似一级修正两种方法解析研究双JC 原子的纠缠演化。
通过与数值精确解【4】的比较，我们发现转动波近似一级修正相对于推广转动波近似更加合适用来描述强耦合区间的行为。
2 理论模型与计算方法两个独立的超导量子比特与各自 LC 振子耦合系统可用下面哈密顿量表示 ††,,,()2k k z k k k k x total k A B k A B k A B H a a g a a σωσ===∆=+++∑∑∑(1) 其中∆ 是二态的能级劈裂，kz σ是第k 个量子比特的 z 方向的泡利算符，†k a 是从第 k 个LC 振子的光子的产生算符，k a 是相应的湮灭算符，g 是原子和LC 振子的耦合强度,ω为光子的频率。
采用自然单位 1=。
为了研究在两个无相关作用的JC 系统的纠缠动力学, 我们只需要给出单个二能级原子和光场耦合系统的本征波函数和能量。
推广的转动波近似方法(简称GRWA)【6】 在很大耦合范围内给出的JC 原子能谱与数值精确解符合很好，引起很多关注。
我们简要回顾下GRWA 求解的方法【6】。
首先把单个JC 哈密顿量绕y 轴旋转2π，得到 ††'()2x z H a a g a a σωσ∆=-+++ (2)然后对上式做幺正变换†exp[()]z U g a a σ=--. 仅保留第一阶激发项，可以得到GRWA 下的哈密顿量，††20[()]()2GRWA z G n H a a g Rr a a ηωσσσ-+∆-=+-++ (3)这里的系数0()G n 为收稿日期： 4 †222202()|cosh[()]|exp(2/)(4/)n gG n n a a n g L g ωωω=〈-〉=- (4)其中n L 为拉盖尔多项式，0(0)G η=,()2Rr F n ∆=。
由于系统总粒子数在哈密顿量(3)下守恒， 我们可以在基矢空间下 |,n ↑〉 和 |,1n ↓+〉 下对其对角化，得到它的本征能量222222242222221221111[(4)(4)]{[(4)(4)]}[(4)]2424(1)g g g n n n n n n e e g e E n g L g L g L g L g L g n ---±++∆∆∆=+-+-±-+++ (6) 同时其相应的本征波函数为 |cos||sin |1|22n n n θθϕ+〉=〉↑〉++〉↓〉 (7) |sin ||cos |1|22n n n θθϕ-〉=〉↑〉-+〉↓〉 (8)基态为 |,0↓〉 , 基态能为 22GS E g η∆=--(9) 这样我们通过GRWA 方法获得了单个JC 模型的本征波函数和本征值。
最近,我们提出了另一种更简洁的哈密顿量(3)的本征解. 称为RWA 的一级校正,简称CRW A [9]. 相应的本征能量为 2311(34)m 124321m m E m g m g g m -+=+-+++++ (15)23121(m 2)(32)m 124321m m E m g m g g m ++++=+++--++ (16)其中, -, + 对应奇偶宇称，对应的本征波函数为11|||1|1m m m m c m c m c m ϕ--+〉⎛⎫〉∝ ⎪-〉++〉⎝⎭ (17)21||2||1m m m m c m c m c m ϕ+++〉++〉⎛⎫〉∝ ⎪+〉⎝⎭ (18)为了研究通常的纠缠演化,初态通常选为光子真空态、关联或反关联的Bell 态,这些演化只会涉及基态和低激发态.收稿日期： 53 纠缠演化根据两种近似解析方法给出的本征波函数和本征能量，我们就可以求出任意时刻的两个无相互作用JC 原子的约化密度矩阵。
对于两个独立的JC 原子，本征波函数在两种方法下分别为12|||GRWA GRWA GRWA ψϕϕ〉=〉⊗〉 (23)12|||CRWA CRWACRWAψϕϕ〉=〉⊗〉 (24)所以整个系统的在任意时刻的波函数可以写成(i)(t)1||i iE t GRWA GRWA i G e ψψ-=〉=〉∑ (25)(i)(t)1||i iE t CRWA CRWA i C e ψψ-=〉=〉∑ (26) 其中(i)(i),G C 是由体系初始态确定的系数。
这里我们选取原子的具有代表性地反关联贝尔1态和关联贝尔2态为初始态，(1)|cos |sin |Bell ψαα〉=↑↓〉+↓↑〉 (27)(2)|cos |sin |Bell ψαα〉=↑↑〉+↓↓〉 (28) 而初始光子态为真空态. 系数(i)(i),G C 由(j)(0)||,(j 1,2)Bell ψψ〉=〉=确定的。
这样就很容易获得约化密度矩阵(GRWA)ρ和(CRWA)ρ。
两个比特态的纠缠可用并发(Concurrence)[10]来度量. 根据并发的定义, 先计算赝密度矩阵()()y y y y ρσσρσσ=⊗⊗，可以得到并发的表达式为(t)1234max[0,]AB C λλλλ=---其中1234,,,λλλλ是ρ的本征值且以递减顺序排列。
收稿日期：6图1 系统初态(1)(0)||Bellψψ〉=〉下/4απ=纠缠随时间演化，黑色线是数值严格解，蓝色线是 推广的转动波近似，红色是转动波近似的一级修正，绿色是转动波近似.其中1,1ω∆== Fig.1 Concurrence for atom-atom entanglement with the initial atomic state (1)(0)||Bell ψψ〉=〉for /4απ=，theblack curves is the numerially exact results,the blue curves is GRWA,the red curves is the first correlations ofRWA,the green curves is RWA,and 1,1ω∆==收稿日期： 7图2 系统初态(2)(0)||Bell ψψ〉=〉下/12απ=纠缠随时间演化，黑色线是数值严格解，蓝色线是推广的转动波近似，红色是转动波近似一级修正，绿色是转动波近似.其中1,1ω∆== Fig.2 Concurrence for atom-atom entanglement with the initial atomic state (2)(0)||Bell ψψ〉=〉 for /12απ=，the black curves is the numerially exact results,the blue curves is GRWA,the red curves is the first correlations ofRWA,the green curves is RWA,and 1,1ω∆==图1给出了由反关联Bell 态出发的纠缠演化曲线，其中黑线是数值严格解，蓝线是推广的转动波近似，红线是我们的转动波近似的一级修正，绿色是原先的转动波近似 。
![基于超导量子比特和里德堡原子的纠缠态制备方法及装置[发明专利]](https://img.taocdn.com/s1/m/bdb259da31126edb6e1a10a7.png)
