一种通过约瑟夫森结非线性频率响应确定微波耗散的方法*

陈恒杰 薛航 李邵雄 王镇

1) (中国科学院上海微系统与信息技术研究所,信息功能材料国家重点实验室,上海 200050)

2) (中国科学院超导电子学卓越创新中心,上海 200050)

3) (中国科学院大学,北京 100049)

1 引 言

超导约瑟夫森结是超导电子学应用的关键元件之一.尤其是在最近二十年里,它在超导量子信息处理器件中有着至关重要的应用[1].随着结制备水平的进步、材料损耗研究的深入和电路设计的优化,以约瑟夫森结为基础的超导量子比特有了巨大的进步.为了不断提升约瑟夫森结器件水平,建立一种适用于约瑟夫森结微波器件,特别是量子比特器件,而且简单、可靠的结质量评估方法是十分有意义的.

通常测量约瑟夫森结的电流-电压(I-V)关系曲线是一个简单评估结质量的方法.例如对于在超导电子学中有广泛应用的Nb/AlOx/Nb约瑟夫森结,可定义一个量Vm来衡量结的品质.Vm为测得的I-V曲线中2 mV处的能隙内电阻与结临界电流的乘积[2].Vm值越大,表明得到的结质量越好.但是,Vm不直接反映结的高频性质.在结的质量比较好、没有额外漏电流的情况下,能隙内电阻值的大小取决于能隙内低能激发的准粒子数,因此Vm基本反映了结中准粒子的情况.在高频应用中,准粒子会带来能量耗散[3].而超导量子比特器件通常工作在数吉赫兹的频率,它在这一频域上的能量损耗是最受关注的性质之一.高频损耗会导致量子比特的能量弛豫,能量弛豫率1/T1的大小直接正比于能量耗散的大小[4,5].能量弛豫时间T1,相位退相干时间Tø是限制量子比特量子相干时间T2的两个主要因素 (1/T2=1/(2T1) + 1/Tø).因此,测量结的高频耗散,能有效地评估制备的约瑟夫森结在量子比特器件中的适用性.约瑟夫森结中高频耗散主要来自三个方面: 准粒子损耗、辐射损耗、介电损耗[6].其中准粒子损耗会随着温度降低而指数减少[7],而且从I-V曲线中得到的Vm反映的是热平衡下的准粒子平均数,不反映准粒子的涨落情况[8,9],因此Vm不能直接全面表征结的高频性质.高频耗散的准确确定通常需要复杂的动力学建模,不合理的系统建模会给出不可靠的结果[10,11].另外,直接测量能量弛豫时间,会有严格的实验温度要求(通常在mK范围),为此需要极其谨慎地实验排除所有强度可比于该温度涨落的噪声干扰[12,13].

基于结的非线性频率响应,本文提出一种温度限制相对宽松的测量结微波损耗的简单可靠方法.若外加微波频率在约瑟夫森结的线性等离子振荡频率附近,随着微波驱动能量的增强,结因非线性效应产生的共振频率偏离结的线性等离子振荡频率.这一偏离的幅度与振荡经历的耗散直接相关.这种非线性频率响应效应可以被用来确定微波耗散.

本文首先介绍电流偏置约瑟夫森结的动力学描述.采用结的等效势阱的四阶近似,得到约瑟夫森结的非线性响应解析结果.从得到的非线性响应方程中建立结的共振频率-振幅-微波品质因子关系.用数值方法模拟未做势阱四阶近似的电流偏置约瑟夫森结微波动力学行为.数值结果证实了非线性响应方程的定量关系可在约瑟夫森结中应用.最后讨论提出的方法在约瑟夫森结中的实验应用情况.

2 约瑟夫森结中的非线性响应

约瑟夫森结是超导量子比特器件中至关重要的非线性元件.电流偏置的约瑟夫森结在较强微波驱动下会表现出比较强的非线性响应效应.为了能够定量描述这一非线性响应,需要考虑约瑟夫森结动力学具体行为.电流偏置的约瑟夫森结的动力学行为,可以由电阻电容并联结(RCSJ)模型描述[14,15]:

其中j是通过约瑟夫森结的相位差,是约化普朗克常数,e是电荷电量,C是结电容,R是损耗带来的有效阻抗,Ic0是结临界电流,Idc是偏置电流,Irfcos(ωrft)为外加微波驱动.在固定偏置电流(

其中idc=Idc/Ic0是归一化的偏置电流,Q=ωp0RC是品质因子.运动方程(2)式中对应的非线性势能(洗衣板势)为

考虑虚拟粒子是在势阱的极小值附近振动,可将势函数U在极小值φ0=arcsinidc附近展开,取四阶近似,运动方程(2)式可简化为

发生共振时,振幅b达到最大值,

因此结线性等离子频率与非线性共振频率的差 ∆ω=γres− ωp与微波品质因子Q有一个直接简单的关系:

利用这个关系,可以通过测量不同强度微波驱动下的频率差值来确定结的微波品质因子.

3 数值模拟验证

由于(7)式是采用了势阱近似得到的,为了验证这一关系可以被定量地应用于约瑟夫森结中,我们对微波驱动的约瑟夫森结动力学行为进行了数值模拟.采用四阶Runge-Kutta方法数值求解方程(1),式中势能U未做前述四阶近似处理.数值计算步长为0.05.为获得虚拟相位粒子在相空间中完整的运动图像,总计算步长数为1800.选取常见的结面积为4 µm2的结,其他所有用到的结参数值见表1.表2列出了所有数值模拟用到的实验条件参数值.其中结的直流偏置选取在0.5左右是考虑到实验在热激发逃离势阱为主导的高温区进行,相较于量子隧穿逃离为主的低温区,高温区的势阱需要更深才能容许较大的微波驱动[17].由于ωp=(1−不同的直流偏置电流选择会导致工作频率的不同.但是我们的数值模拟计算中采用了归一化因子ωp0,实验上总是可以通过磁场调节结的临界电流值达到将工作频率调整到一个感兴趣范围的目的.

表1 数值模拟采用的结参数Table 1.Parameters of Josephson junctions used in numerical simulations.

表2 数值模拟采用的实验参数Table 2.Experiment settings used in numerical simulations.

图1为数值计算得到的一个典型的虚拟相位粒子在相空间的运动轨迹.由于模拟计算时采用的初始条件不一定处于相位粒子的稳定态,采用的总计算步长数必须保证粒子运动轨迹能够达到一个不变稳定态.稳定态的运动轨迹(例如图1中内圈轨迹被重复了很多遍因而更密集)可被认定是方程(1)的稳定解.运动轨迹在位置坐标上的变化幅度对应粒子的振动幅度.

图1 计算得到的相位粒子在相空间中的运动轨迹 纵坐标 v=dϕ/dτ,计算采用的实验参数为idc=0.473,irf=5.5×10−4,g=0.9306Fig.1.Calculated trajectory of phase particle with experiment parameters for junction 3: idc=0.473,irf=5.5×10–4,g=0.9306.Vertical axis v=dϕ/dτ.

微波频率响应通常需要变化外加微波频率.但是实验上测量约瑟夫森结的微波响应时,选取固定外加微波频率,改变偏置电流idc的方法(即改变ωp)可以保证耦合到结的实际微波功率的可定标性.为方便与实验结果比较,在不同的微波强度驱动条件下,在计算非线性响应时只改变idc,这对应于(5)式中固定g,改变ωp.图2给出了数值模拟得到的结非线性响应结果与(5)式计算结果的比较,可以看到两者定量符合.图2中,响应幅度b最大处对应约瑟夫森结发生共振吸收.发生共振时共振频率γres=γ,对应的ωp由最大振幅处的idc决定.从计算的结果可见当外加微波强度增大,约瑟夫森结非线性共振时对应的idc变小,从而相应的ωp变大,这意味着共振频率相对更小,这说明较强微波驱动下约瑟夫森结振子会软化.为进一步验证由非线性频率响应确定微波品质因子的准确性及在约瑟夫森结中的适用性,图3给出了表1中三个样品结模拟得到的结果与(7)式结果的比较.考虑到同样工艺制备的结可能被运用到不同的器件,遇到不同的电路环境,图4为结2在不同环境带来不同损耗情况下的模拟结果比较.同样,(7)式与模拟结果定量符合.从图3和图4的结果可知,只要实验确定了约瑟夫森结的共振峰对应的idc、结临界电流、结电容参数,结的微波品质便可定量确定(图3和图4中直线斜率大小对应微波品质因子).由于模拟选取了常见的约瑟夫森结参数,而且在模拟的品质因子变化范围内的约瑟夫森结中已有实验观测到了量子相干态[18],我们认为由非线性频率响应来确定微波品质因子的方法完全适用于约瑟夫森结系统,且可用来评估制备的约瑟夫森结在相关器件应用中的微波性能.

图2 结 (IC=8 µA,Q=515.7)在不同的微波驱动下,(5)式(曲线)和模拟(误差棒)得到的响应振幅b随归一直流偏置 idc 的函数变化关系,响应曲线按其最大振幅从小到大顺序对应的外加微波强度分别为irf=1.63×10−4,2.44×10−4,3.67×10−4 和 5 .50×10−4,点划线显示了最大振幅对应的 idc 随外加微波强度的变化Fig.2.Microwave response curves obtained by Eq.(5)(curves) and numerical simulation (error bars),for junction(IC=8 µA,Q=515.7) with applied microwave irf=1.63×10−4,2.44×10−4,3.67×10−4and5.50×10−4 for the curves with the maximum amplitude from small to large respectively.Dot-dash line shows the dependence of idc where corresponding to the maximum oscillation amplitude on the power of the applied microwave.

图3 表1 中三个样品结的随 irf的变化关系图,线 是 (7) 式的结果,误 差棒是 v˜ 的 数 值模拟结 果,图中直线斜率从小到大分别对应Q值163.0(结1),508.6(结2),1549.1(结 3)Fig.3. as a function of irf for different parameters of sample Josephson junctions in Table 1.Lines are results of Eq.(7).Error bars are numerical simulation results of v˜.The lines of slops from small to large corresponding to Q values: 163.0 (junction 1),508.6 (junction 2),1549.1 (junction 3).

图4 结 2 在不同的电路环境下,导致不同的品质因子 Q(51.6,257.8,515.7,1533.0,对 应 直 线 斜 率 从 小 到 大) 时,v˜随 irf 的变化关系,直线是(7)式的结果,误差棒是 v˜ 的数值模拟结果Fig.4.v˜ as a function of irf for different quality factors.Lines are Eq.(7)’s results.Error bars are numerical simulation results of v˜.The lines of slops from small to large corresponding toQvalues: 51.6,257.8,515.7,1533.0accounting for junction 2 with different environment influences.

4 讨 论

对于过阻尼的约瑟夫森结,其线性等离子振动频率(ωp)也受微波损耗影响很大.但是大多数约瑟夫森微波器件,尤其是量子比特器件,要求比较高的微波品质 (Q≫1).欠阻尼约瑟夫森结的等离子振动频率几乎不随微波品质变化.对于欠阻尼约瑟夫森结,将其驱动到较强的非线性区域,可由非线性频率响应来确定微波品质.虽然约瑟夫森结的线性微波响应也与微波损耗相关[19],但是与线性响应相比,非线性响应信号更强,不容易受额外的噪声干扰,可以在相对较高的温区准确测量.图3和图4的结果显示,结的微波非线性响应(图中直线斜率)随微波品质变化很大,已有的实验结果也已经显示微波品质因子严重影响非线性振动[20].Mao等[21]用数值模拟比对由较强微波驱动的、测量温度在4.2 K的共振激发实验数据得到的约瑟夫森结高频等效阻抗.因为约瑟夫森结的共振频率、临界电流、结电容等参数都是实验可测量的,我们提出的非线性频率响应模型方法不需要数值模拟便可直接实验确定微波品质因子.而且由(7)式可知微波品质因子的定量准确性也是主要由结参数—临界电流、结电容参数的准确性保证的[13,22].考虑到共振激发的实验测量可以适用于单个电流偏置的约瑟夫森结中[19,21],不需要额外的电路设计,数值模拟也证实了非线性响应模型在约瑟夫森结中的定量准确性,因此利用非线性频率响应来测量微波损耗的方法是方便可靠的.