两个超导比特纠缠动力学的反转动波效应

收稿日期：基金项目：国家自然科学基金(11174254)和浙江省自然科学基金（Z7080203）和浙江师范大学创新研究团队项目资助的课题．作者简介：钱雷(1987—)，男，硕士研究生，主要从事量子光学研究.*通信作者，Email：[e-mail].1两个超导比特纠缠动力学的反转动波效应*钱雷1贺树2段立伟2陈庆虎1,2*（1.浙江师范大学,海峡两岸统计物理与凝聚态理论研究中心,金华321004;2.浙江大学物理系,杭州310027）摘要：本文研究了两个无相互作用的与各自量子振子耦合的超导量子比特的纠缠动力学。数值严格结果证明在强耦合下,以往的转动波近似不再适用，非转动波效应必须考虑。基于著名的推广的转动波近似的解析结果,在目前实验上可实现的强耦合区域与严格数值结果有明显差别。从我们近年来发展的转动波近似的一级校正，我们得到的结果更接近于精确数值解。该理论结果可激发基于近来可以实现超导比特的相关实验。关键词：纠缠动力学;量子比特;比特-振子耦合;反转动波.QIANLei1,HEShu2,DUANLiwei2,CHENQinghu1,2*（1.CenterforStatisticalandTheoreticalCondensedMatterPhysics,ZhejiangNormalUniversity,Jinhua321004,China;2.DepartmentofPhysics,ZhejiangUniversity,Hangzhou310027,China)Title:EffectofCounter-RotatingtermontheentanglementdynamicsoftwosuperconductingqubitscouplingtoquantumoscillatorsAbstract:Inthispaper,weinvestigateentanglementdynamicsoftwonon-interactingsuperconductingqubitscoupledwiththeirownquantumoscillatorsbeyondtherotating-waveapproximation(RWA).Itisshowninthenumericallyexactstudiesthatinthestrongcouplingregime,theRWAisnotvalidandtheeffectofthecounter-rotatingwaveshouldbetakenintoaccount.Theanalyticresultsforentanglementbasedonthewell-knowngeneralizedRWAdeviatesfromthenumericaloneobviouslyinthepresent-dayexperimentallyaccessiblecouplingregime.Inthis收稿日期：基金项目：国家自然科学基金(11174254)和浙江省自然科学基金（Z7080203）和浙江师范大学创新研究团队项目资助的课题．作者简介：钱雷(1987—)，男，硕士研究生，主要从事量子光学研究.*通信作者，Email：[e-mail].2paper,basedonourfirst-ordercorrectionstotheRWA,theanalyticresultsaremoveclosetothenumericalones.Itmaymotivatetherelevantexperimentsbasedonthesuperconductingqubitsrealizedrecently.KeyWords:Entanglementdynamics;qubits;qubit-oscillatorcoupling;Counter-rotatingwave.1引言超导比特作为量子信息科学的基础和量子计算机的重要元件,在基础研究和实际应用中受到了广泛的关注。近年来，二能级系统与LC谐振子耦合系统的动力学研究在实验上取得重要进展，对量子计算，量子相干性和量子逻辑门等领域的发展有较大的促进。量子纠缠是一种典型的量子关联,描述了量子力学的本质的非局域性质,根本性地区分了量子和经典物理学。纠缠作为一种重要的量子信息资源,受到越来越多的重视，而由于与环境耦合导致退相干和纠缠突然死亡（纠缠在一段时间消失）等问题是目前量子信息领域内研究的热点课题。Jaynes-Cummings（JC）模型【1】是光学腔量子电动力学(CQED)系统中最简单的一类模型,描述腔场中单个二能级原子和单模量子化的电磁场的相互作用问题。在转动波近似下,该模型可以严格解析求解。超出转动波近似,闭合形式的解析解至今缺乏,各种近似的解析解相继提出,其中最为著名的是推广的转动波近似【6】.Yu等人研究了在转动波近似下两个JC原子的量子纠缠动力学【2，3】,并给出了相应的解析表达式。陈等人利用数值方法研究了无转动波近似下的纠缠行为【4】,发现反转动波对纠缠有压制效应。至今，未见到反转动波下的纠缠动力学的解析研究或近似解析收稿日期：基金项目：国家自然科学基金(11174254)和浙江省自然科学基金（Z7080203）和浙江师范大学创新研究团队项目资助的课题．作者简介：钱雷(1987—)，男，硕士研究生，主要从事量子光学研究.*通信作者，Email：[e-mail].3研究。近期的超导量子比特与振子耦合系统的实验表明,当其耦合强度达到g=0.1（这比量子光学系统里的二能级原子和腔耦合系统的耦合强度大3-6个数量级），转动波近似已被证明实效.因此,纠缠动力学的反转动波的效应的解析研究是十分必须的和值得的.在本文中，我们基于推广转动波近似和转动波近似一级修正两种方法解析研究双JC原子的纠缠演化。通过与数值精确解【4】的比较，我们发现转动波近似一级修正相对于推广转动波近似更加合适用来描述强耦合区间的行为。2理论模型与计算方法两个独立的超导量子比特与各自LC振子耦合系统可用下面哈密顿量表示††,,,()2kkzkkkkxtotalkABkABkABHaagaa(1)其中是二态的能级劈裂，kz是第k个量子比特的z方向的泡利算符，†ka是从第k个LC振子的光子的产生算符，ka是相应的湮灭算符，g是原子和LC振子的耦合强度,ω为光子的频率。采用自然单位1。为了研究在两个无相关作用的JC系统的纠缠动力学,我们只需要给出单个二能级原子和光场耦合系统的本征波函数和能量。推广的转动波近似方法(简称GRWA)【6】在很大耦合范围内给出的JC原子能谱与数值精确解符合很好，引起很多关注。我们简要回顾下GRWA求解的方法【6】。首先把单个JC哈密顿量绕y轴旋转2，得到††'()2xzHaagaa(2)然后对上式做幺正变换†exp[()]zUgaa.仅保留第一阶激发项，可以得到GRWA下的哈密顿量，††20[()]()2GRWAzGnHaagRraa(3)这里的系数0()Gn为收稿日期：基金项目：国家自然科学基金(11174254)和浙江省自然科学基金（Z7080203）和浙江师范大学创新研究团队项目资助的课题．作者简介：钱雷(1987—)，男，硕士研究生，主要从事量子光学研究.*通信作者，Email：[e-mail].4†222202()|cosh[()]|exp(2/)(4/)ngGnnaangLg(4)其中nL为拉盖尔多项式，0(0)G,()2RrFn。由于系统总粒子数在哈密顿量(3)下守恒，我们可以在基矢空间下|,n和|,1n下对其对角化，得到它的本征能量222222242222221221111[(4)(4)]{[(4)(4)]}[(4)]2424(1)gggnnnnnneegeEngLgLgLgLgLgn(6)同时其相应的本征波函数为|cos||sin|1|22nnn(7)|sin||cos|1|22nnn(8)基态为|,0,基态能为22GSEg(9)这样我们通过GRWA方法获得了单个JC模型的本征波函数和本征值。最近,我们提出了另一种更简洁的哈密顿量(3)的本征解.称为RWA的一级校正,简称CRWA[9].相应的本征能量为2311(34)m124321mmEmgmggm(15)23121(m2)(32)m124321mmEmgmggm(16)其中,-,+对应奇偶宇称，对应的本征波函数为11|||1|1mmmmcmcmcm(17)21||2||1mmmmcmcmcm(18)为了研究通常的纠缠演化,初态通常选为光子真空态、关联或反关联的Bell态,这些演化只会涉及基态和低激发态.收稿日期：基金项目：国家自然科学基金(11174254)和浙江省自然科学基金（Z7080203）和浙江师范大学创新研究团队项目资助的课题．作者简介：钱雷(1987—)，男，硕士研究生，主要从事量子光学研究.*通信作者，Email：[e-mail].53纠缠演化根据两种近似解析方法给出的本征波函数和本征能量，我们就可以求出任意时刻的两个无相互作用JC原子的约化密度矩阵。对于两个独立的JC原子，本征波函数在两种方法下分别为12|||GRWAGRWAGRWA(23)12|||CRWACRWACRWA(24)所以整个系统的在任意时刻的波函数可以写成(i)(t)1||iiEtGRWAGRWAiGe(25)(i)(t)1||iiEtCRWACRWAiCe(26)其中(i)(i),GC是由体系初始态确定的系数。这里我们选取原子的具有代表性地反关联贝尔1态和关联贝尔2态为初始态，(1)|cos|sin|Bell(27)(2)|cos|sin|Bell(28)而初始光子态为真空态.系数(i)(i),GC由(j)(0)||,(j1,2)Bell确定的。这样就很容易获得约化密度矩阵(GRWA)和(CRWA)。两个比特态的纠缠可用并发(Concurrence)[10]来度量.根据并发的定义,先计算赝密度矩阵()()yyyy，可以得到并发的表达式为(t)1234max[0,]ABC其中1234,,,是的本征值且以递减顺序排列。收稿日期：基金项目：国家自然科学基金(11174254)和浙江省自然科学基金（Z7080203）和浙江师范大学创新研究团队项目资助的课题．作者简介：钱雷(1987—)，男，硕士研究生，主要从事量子光学研究.*通信作者，Email：[e-mail].6图1系统初态(1)(0)||Bell下/4纠缠随时间演化，黑色线是数值严格解，蓝色线是推广的转动波近似，红色是转动波近似的一级修正，绿色是转动波近似.其中1,11,1Fig.1Concurrenceforatom-atomentanglementwiththeinitialatomicstate(1)(0)||Bellfor/4，theblackcurvesisthenumeriallyexactresults,thebluecurvesisGRWA,theredcurvesisthefirstcorrelationsofRWA,thegreencurvesisRWA,and1,1收稿日期：基金项目：国家自然科学基金(11174254)和浙江省自然科学基金（Z70802