纳米机械谐振器耦合量子比特非厄米哈密顿量诱导的声子阻塞

廖庆洪 邓伟灿 文健 周南润 刘念华

1) (南昌大学电子信息工程系,南昌 330031)

2) (清华大学物理系,低维量子物理国家重点实验室,北京 100084)

3) (南昌大学高等研究院,南昌 330031)

1 引 言

超导量子电路作为人工设计的量子力学系统,显示了宏观的量子相干性,可以用来研究量子光学和原子物理的现象[1].在超导量子电路中,超导传输线谐振器代替传统的光学腔,超导量子比特扮演自然原子的角色.超导量子比特作为宏观的人工原子,其他超导量子器件,例如超导传输线、机械纳米振子等的强耦合,可以实现一种崭新的全固态的腔量子电动力学(cavity quantum electrodynamics,C-QED)结构,而且这些人工的原子和光场,具有以前微腔QED没有的可控性和可集成化的优点.包含超导量子比特和纳米机械振子的混合量子电路是另一种新兴的混合量子系统[2].这种混合量子电路的有趣的理论研究和有实质性进展的实验结果,证明了该系统可以作为未来量子器件,并拥有广泛的应用[2−4].

2010年,加州大学O'Connell等[5]发现机械谐振器在超低温(超导)条件下的振动可以激发量子效应,产生单个量子(声子).随后,清华大学Liu等[6]根据这个实验,研究了纳米机械谐振器(nanomechanical resonator,NAMR)中的声子阻塞效应.声子阻塞,顾名思义,类似于光子阻塞[7],是一种纯量子现象.这种现象是在非线性机械振荡器[8]中,第一个声子被外部驱动场激发,第二个声子被第一个声子阻挡导致不能传输.光子阻塞效应的一项重要应用是可以用于实现理想的单光子源[9,10],而声子阻塞效应同样可以作为实现单声子源的机制之一[11,12].类比于产生光子阻塞的两种机制,声子阻塞同样地对应两种发生机制.声子阻塞的实现有两种主要的物理机理: 基于失谐机制的声子反聚束需要较大的非线性来改变系统的能级结构,这被称为传统声子阻塞(conventional phonon blockade,CPNB)[13−16].除了第一种物理机理外,还可以获得强声子反聚束,即利用双声子激发的不同路径产生的破坏性干涉,即使在非线性远小于腔模的衰减率(弱非线性)的条件下也能实现,称为非传统声子阻塞(unconventional phonon blockade,UPNB)[16−19].

自发现声子阻塞效应起,国内外人员对这种物理现象进行了大量的研究.例如双光力系统中,光子和声子不同耦合路径的干涉产生的光谱差别[20]、光学模式与机械模式之间产生的二阶非线性相互作用导致的强声子反聚束现象[21]、机械谐振器与两能级系统耦合的光力系统的双色二阶边带的生成[22].多模光力系统[23]、有限温度下耦合非线性微米/纳米机电系统的谐振器的声子反聚束效应并与热噪声的关系[24]、机械振动位移的平方耦合的声子阻塞现象研究[25]、光声阻塞现象[26]、俘获多个离子的库仑晶体中声子阻塞研究[27]等.

最近,关于非厄米哈密顿量的系统特性研究得到大量关注[28−44].其中包含了量子纠缠[28]、兰姆位移[31]、光子统计特性[44]等.本文研究 NAMR-量子比特系统中的声子统计性质,其中的系统是由一个纳米机械谐振器耦合一个超导量子比特(以下简称量子比特)组成.这里利用了非厄米哈密顿量的零本征值方法,并假设κ<0,研究了 NAMR-量子比特系统中的声子统计特性.同时给出探测声子阻塞或声子统计特性的方法.

2 NAMR-量子比特系统的物理模型

图1研究的是由一个纳米机械谐振器(右边)耦合一个量子比特(左边)组成的物理系统.纳米机械谐振器受到一个强度为ε的外加场驱动,量子比特受到一个强度为Ω的外加场驱动.整个系统的哈密顿量为[16]

其中 ℏ=1,σ+和σ−是量子比特的上升和下降算符,其跃迁频率是ω0;a和a+表示 NAMR 中声子的湮灭和产生算符,其频率是ωm;g是量子比特与NAMR的耦合强度; 量子比特的外加驱动场振幅是Ω,频率是ωq; NAMR的外加驱动场振幅是ε,频率是ωa;ϕ是两个外加驱动场的相位差.

假设两个外加驱动场的频率相同且量子比特与 NAMR 的频率相同,即ωq=ωa=ωd和ω0=ωm,在以入射光频率旋转的框架中,得出哈密顿量:

这里失谐量∆≡ω0−ωd=ωm−ωd.

图1 NAMR-量子比特系统的物理模型图Fig.1.Schematic diagram of the nanomechanical resonator coupled with a qubit.

3 非厄米哈密顿量作用下的声子阻塞效应的零本征值

3.1 声子阻塞的条件分析

考虑机械模衰减率κ和量子比特自发辐射速率γ.所以,整个系统可以用非厄米哈密顿量描述为

考虑系统总激发数为n+1,在基矢和下,可得到非厄米哈密顿量的矩阵形式:

考虑弱驱动的条件,即系统参数Ω,ε≪|κ|.

分析(4)式中的矩阵,可得知其两个特征值:

在弱驱动条件下,考虑系统激发数为1,即n=0.而且要使得(5)式特征值取到实数,则满足条件

此时,两个特征值可以化简为

当∆=∓g时,这一对特征值相等并且取得最小值0.此时可以知道特征值为零的条件为:

由图2中的系统能级可知,系统参数条件满足(8)式或(9)式时,可出现声子阻塞现象.

在我们的方案中,γ设为正数,κ设为负数.事实上,在光学腔中,腔模衰减率κC<0 已被广泛地研究[45,46],而且已经有实验实现了这个参数[47−51].本文基于此,做出了相应的假设,即声子场衰减系数κ<0.实际上已有实验通过相干驱动声子技术[5,52−54]、引入声子激光[55−57]和直接驱动机械振子[58,59]等方式来提供机械增益.这里简要讨论直接驱动纳米机械谐振器来产生力学增益的可能性[60].外场驱动在可分辨边带条件下,对应于纳米机械谐振器中声子冷却的斯托克斯散射过程被极大抑制,给纳米机械谐振器提供能量的反斯托克斯过程得到了增强并占主要优势.

3.2 数值模拟方法

通过计算系统的主方程来研究二阶关联函数.系统的密度矩阵为:

只需关注稳态时的情况,即∂ρ/∂t=0.通过二阶关联函数g(2)(0) 来研究声子的统计特性.

3.3 结果分析

下面通过计算机数值模拟分析研究声子阻塞效应,比较和分析(8)式和(9)式.系统中腔模的希尔伯特空间被截取到5维,以耗散率γ作为基本单位来设定其他参数.图3给出了不同条件下(g/γ=0.5,2.0,4.0)g(2)(0) 与∆/γ的关系图,可以得出出现声子阻塞的特点.

1)g/γ=0.5,当出现g(2)(0) 的波谷,即g(2)(0)≪1时,对应于强声子反聚束现象.此时∆/γ=0,与 (8)式相对应.

2)g/γ=2,当出现g(2)(0)≪1 时,此时∆/γ=±2,对应于 (9) 式.也就是说,在条件g/γ=0.5 下比在条件g/γ=2 下出现的反聚束现象更强.这是因为这里只取了的一个近似值.只要随着g的增大,近似效果越明显,反聚束效应就会越强,例如图3中的g/γ=4.

图3 在不同的 g /γ 关系下,二阶关联函数 g(2)(0) 与 ∆/γ的关系图 (a) κ /γ=−1,ε /γ=0.1,Ω=0,ε=0 ;(b) κ /γ=−1,Ω /γ=0.1,ε=0,ϕ=0Fig.3.Logarithmic plot (of base e) of the zero-delay-time second-order correlation functions g(2)(0) as a function of∆/γ for different g /γ with κ /γ= −1 : (a) ε /γ=0.1,Ω=0,ε=0 ; (b) Ω /γ=0.1,ε=0,ϕ=0.

3.4 与CPNB和UPNB比较

现在讨论在非厄米特哈密顿量作用下产生声子阻塞效应在零本征值下的物理机理.

图4描述了基于不同路径干涉的情况下的声子阻塞现象的能级情况.绿色箭头描述的是第一条路径蓝色箭头是第二条路径其中分别是衰 减到衰减到的衰减率.

物理原理上分析,能级分裂决定了声子阻塞效应的类型是UPNB还是CPNB.对于(8)式或者(9) 式的情况,驱动场频率ωd如果等于那就满足了单激发数共振条件,同时产生单声子的可能性就大大增加.在这里,由于失谐,产生双声子的可能性并没有像单声子那样大.

图4 UPNB 能级图Fig.4.Schematic energy-level diagram of unconventional phonon blockade.

在少声子数近似下,系统的态矢量可写作[16]:

在弱驱动的条件下,二阶关联函数

图5比较了CPNB,UPNB和本文方案的区别.不管是在强非线性还是弱非线性的限制下,本文方案都可以实现声子的反聚束现象.实现该方案的物理机理和CPNB相同,都是由于失谐所造成的声子的反聚束.实现CPNB的最佳条件如下:

这些最佳条件同时也分析了κ与γ之间的定量关系.这是因为仅仅考虑了(2)式中的特征频率,而没有引入κ与γ.这就导致了在∆=|g| 的条件下,不能得到 |g|=γ/2 这个条件(见图6),而只能注意到∆=|g|=0.本文方案最显著的特点就是加入了κ/γ=−1这个条件.事实上,如果(9)式中的条件变为κ>0 和γ>0,那么就变成了 (13) 式,也就是说变为了CPNB.

图5 本文方案与 CPNB 和 UPNB 的区别Fig.5.Different parameter range for the CPNB,the UPNB,and the present scheme.

图6 二阶关联函数 g(2)(0) 和平均声子数的对数图像 (a),(b) ε /γ=0.1,Ω=0 ; (c),(d) Ω /γ=0.1,ε=0,ϕ=0.其他参数为 κ /γ=−1.Fig.6.The logarithmic plot (of base e) of the zero-delay-time second-order correlation functions g(2)(0) and average photon numberN: (a) and (b) ε /γ=0.1,Ω=0 ; (c) and (d) Ω /γ=0.1,ε=0,ϕ=0.The shared parameters: κ /γ=−1.

图7 二阶关联函数 g(2)(0) (对数) 与 ∆ /γ 的关系图,其中 ε /γ=0.1,Ω=0 (a) g /γ=0.5,κ /γ=−1 (本文方案),κ/γ=1(UPNB); (b) g /γ=5,κ /γ= −1 (本文方案),κ /γ=1 (CPNB)Fig.7.The logarithmic plot (of base e) of the zero-delay-time second-order correlation functions g(2)(0) as a function of ∆ /γ.The shared parameters: ε /γ=0.1,Ω=0.Other parameters: (a) g /γ=0.5,κ /γ=−1 (our scheme),κ /γ=1 (UPNB); (b) g /γ=5,κ/γ=−1 (our scheme),κ /γ=1 (CPNB).

图8 二阶关联函数 g(2)(0) (对数) 与 ∆ /γ 的关系图,Ω=0,ε=0,ϕ=0 (a) Ω /γ=0.1,g /γ=0.5,κ /γ=−1 (本文方案),κ/γ=1 (UPNB); (b) Ω /γ=2,g /γ=10,κ /γ=−1 (本文方案),κ /γ=1 (CPNB)Fig.8.The logarithmic plot (of base e) of the zero-delay-time second-order correlation functions g(2)(0) as a function of ∆ /γ :Ω=0,ε=0,ϕ=0.Other parameters: (a) Ω /γ=0.1,g /γ=0.5,κ /γ=−1 (our scheme),κ /γ=1 (UPNB); (b) Ω /γ=2,g/γ=10,κ /γ= −1 (our scheme),κ /γ=1 (CPNB).

图9 声子统计特性探测模型Fig.9.The detecting model of the statistical properties of the phonons.

图7和图8比较了本文方案和CPNB,UPNB的不同.结果表明,在设定相应的参数下,本文方案得出的声子阻塞效果比其他两种方案要好.这是因为非厄米哈密顿量的作用方式更贴合主方程.这就导致了在其他条件相同时,本文方案会比直接使用(2)式的更精确.CPNB的触发是由于声子场频率与驱动场频率不同步(失谐)而导致的,而UPNB的触发是通过两个不同路径的量子干涉产生的.因此,当作用方式是量子干涉时,声子的反聚束会对干涉中的相位差极其敏感,所以干涉路径的改变会使反聚束的区域发生剧烈变化.然而,在本文方案中,驱动场的改变导致量子干涉路径的改变并不会对反聚束的物理环境产生影响.所以本方案中的强声子反聚束与UPNB方案的物理机理不同.

4 声子阻塞效应的探测

通过测量光学腔场的输出光子的相关性,可以间接观测到NAMR中的声子阻塞现象.对于NAMR-量子比特系统,我们参照了文献中的测量方法[16].在最近的实验中,一个量子比特耦合一个NAMR以后,NAMR再耦合一个超导腔组成了一个三体系统[61].假设这个超导光腔通过辐射压相互作用耦合到NAMR,就可以通过研究超导光腔中的光子统计特性,来研究NAMR中的声子统计特性,从而观察NAMR中的声子阻塞现象.声子统计特性探测模型见图9.

假设超导腔的驱动频率为ωb,在以入射光频率ωb旋转的框架中,整个系统的总哈密顿量为

其中b和b+表示超导腔中光子湮灭和产生算符,超导腔的共振频率为ωc;J是超导腔与NAMR的耦合强度;Ωc(假定为实数)是超导腔的外加驱动场振幅;∆b满足共振条件∆b≡ ωc− ωb=ωm.

超导腔的腔模算符可以写成量子波动算符和稳态平均场之和:b→b+α.这里α在频率为ωb驱动光场中的稳态平均场.在强驱动的条件下,|α|≫1.同样地,可以将Hom再在旋波近似下,哈密顿量变为其中光力系统有效的耦合强度G=gα,同样假定|G|≪ωm,在相互作用哈密顿量的演化下,声子场和腔场之间形成纠缠态,探测腔场光子统计特性可以实现对声子阻塞或声子统计特性的探测[61].

5 结 论

研究了在非厄米哈密顿量作用下的NAMR-量子比特系统的零本征值条件下的声子统计特性.充分利用了零本征值方法,通过分析和计算,得到了强声子反聚束的最佳条件,同时引入二阶关联函数讨论了声子阻塞效应.通过主方程的稳态求解,分析了各种条件对声子阻塞的影响,结果发现都能很好地匹配声子阻塞现象.比较了本方案与CPNB和UPNB的区别,发现本文方案在强非线性和弱非线性条件下都能实现强声子反聚束现象.阐明了这种声子阻塞的物理机理是不同于UPNB的失谐原理.最后,简要分析了利用外加超导腔去间接地探测NAMR中的声子统计特性的方法.本文的研究结果将在量子信息处理以及精密测量等研究领域具有重要的意义和应用前景.