超导约瑟夫森结物理参数的实验推算

韩金舸; 欧阳鹏辉; 李恩平; 王轶文; 韦联福

doi:10.7498/aps.70.20210393

摘要

Abstract

作者及机构信息

Authors and contacts

文章全文

1. 引　言
2. 约瑟夫森结物理参数的推算方法
2.1 结电容和结临界电流的推算方法
2.2 结阻尼参数的推算方法
3. 约瑟夫森结的制备及其物理参数推算
3.1 约瑟夫森结的制备
3.2 约瑟夫森结I-V曲线和跳变电流的测量
3.3 约瑟夫森结物理参数的推算
4. 结　论

参考文献

施引文献

摘要

超导约瑟夫森结是实现超导量子计算和微波单光子探测的核心器件, 其物理参数很难直接测定. 与之前常用的测量结微波激励效应估计方法不同, 本文通过实验测量低频电流驱动下的约瑟夫森结I-V曲线及其跳变电流统计分布, 并与基于标准电阻电容分路结模型数值模拟进行比对, 推算出了约瑟夫森结的临界电流

$I_{\rm c}$ 、电容C、电阻R及阻尼参数

$\beta_{\rm c}$ 等物理参数. 结果表明, 所推算的参数值与基于微观理论推导所得到的Ambgaokar-Baratoff公式基本符合, 可供约瑟夫森结的器件参数按需设计和制备工艺的参数设置等参考.

关键词:

Superconducting Josephson junctions are the key devices for superconducting quantum computation and microwave single photon detection. It is important to fabricate the Josephson junctions with designable parameters. Different from the typical methods to calibrate the parameters of the Josephson junctions,, e.g., by using the microwave drivings and measuring the ratio of hysteresis current to critical one, in this paper we achieve the calibrations with the low frequency current biases. First, we measure the I-V characteristic curves of the fabricated Al/AlO_x/Al junctions. Second, we measure the statistical distributions of the jump currents of the Josephson junction samples driven by the low frequency (@71.3 Hz) biased currents at an extremely low temperature of 50 mK. These two sets of experimental data are utilized to estimate the typical parameters of the Josephson junction, i.e., junction capacitance, critical current, and the damping coefficient, which are difficult to be directly measured in the usual experiments. The critical current and capacitance of the Josephson junction are estimated by fitting the statistical distribution of the measured jump currents with the relevant theoretical model of the "particle" escape from the potential driven by the thermal excitations and quantum tunnelings. With the calibrated critical current of the junction, the relation between

$I/I_{\rm{c}}$ and

${\rm{d}}\varphi/{\rm{d}}\tau,\,\tau=\omega_{\rm{c}}t$ (with

$\omega_{\rm{c}}$ being the plasmon frequency) is obtained from the measured

$I\text-V$ curve. Using the standard resistively capacitance shunted junction model to fit such a relation, the damping coefficient of the junction can be estimated. With the estimated critical current, capacitance, and damping coefficient, the resistance

$R_n$ of the junction at the working temperature is calibrated consequently. It is shown that our estimated results are in good agreement with that predicted by the famous Ambgaokar-Baratoff formula. Physically, the method demonstrated here possesses two advantages. First, it is relatively insensitive to the noise during the measurement of the junction's I-V characteristic curve, compared with the usual method to calibrate damping coefficient by measuring the ratio of hysteresis current to critical current. Second, only the low frequency driving is required to measure the jump current of the junction for estimating the damping coefficient. The microwave driving is not necessary. Hopefully, the present work is useful for the on-demand designs of the Josephson junctions for various applications.

Keywords:

Josephson junction /
resistively capacitance shunted junction model /
jump current /
parameter estimation

PACS:

03.75.Lm	(Tunneling, Josephson effect, Bose-Einstein condensates in periodic potentials, solitons, vortices, and topological excitations)
06.20.Dk	(Measurement and error theory)
06.20.fb	(Standards and calibration)

作者及机构信息

1.
西南交通大学信息科学与技术学院, 信息量子技术实验室, 成都　610031

2.
西南交通大学物理科学与技术学院, 量子光电实验室, 成都　610031

3.
东华大学理学院, 光子学实验室, 上海　201620

通信作者: 韦联福, lfwei@dhu.edu.cn

基金项目: 国家自然科学基金(批准号: 11974290, 61871333)资助的课题

Authors and contacts

1.
School of Information Science and Technology, Southwest Jiaotong University, Chengdu 610031, China

2.
School of Physics Science and Technology, Southwest Jiaotong University, Chengdu 610031, China

3.
School of Science, Donghua University, Shanghai 201620, China

Corresponding author: Wei Lian-Fu, lfwei@dhu.edu.cn

Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 11974290, 61871333)

文章全文

1. 引　言

目前经典计算机的计算能力渐趋极限, 量子计算机的超强计算能力得到了广泛关注^[1,2]. 如谷歌构建的53位的超导量子“霸权”机和中国科技大学实现的76个光子的“九章”量子计算原型机都展示了经典电子计算机所不可比拟的计算功能^[3,4]. 作为一种广受关注的量子计算实现模式, 超导量子计算机因其特有的易于集成、芯片设计和加工技术相对成熟等方面的优势, 成为量子计算领域的研究热点. 超导量子计算的基本物理单元比特是超导量子比特, 它的三种基本的构型——电荷量子比特、磁通量子比特和相位量子比特等的实现, 取决于核心器件——约瑟夫森结的物理参数. 这些参数决定了结的电容充电能 $E_{\rm c}$ 和约瑟夫森隧穿效应能量 $E_{\rm J}$ 之间的比值, 进而决定了超导量子比特的构型. 此外, 约瑟夫森结也是实现微波单光子探测的核心器件^[5,6].

$E_{\rm c}$ 的大小与约瑟夫森结的电容有关, 而 $E_{\rm J}$ 的取值则由约瑟夫森结的临界电流决定^[7-9]. 这些参数是量子比特设计与制备的重要依据, 但实际上它们很难在实验上直接测量得到, 通常需要根据其他的某些相关参量的测量来推算得到. 之前采用的方法有测量微波驱动下的约瑟夫森结跳变电流统计分布, 通过与理论模型拟合来推算约瑟夫森结各物理参数. 显然, 这一方法有两点不足: 第一, 通过微波激励来进行的参数标定, 实际上是一种动态标定, 有一定的不确定性; 其二, 微波激励本身需要考虑阻抗匹配实现问题. 此外, 还可采用回滞电流与临界电流的比值来标定阻尼参数, 此方法也存在一定的局限性. 当外界干扰造成约瑟夫森结I-V特性曲线零点上下涨落时, 回滞电流与临界电流比值对原点上下涨落比较敏感, 会给标定结果带来较大误差^[10,11]. 为尽可能精确地测定结的各种静态物理参数, 本文提出一种基于标准RCSJ (resistively capacitance shunted junction)模型拟合、仅需要低频电流偏置下的结跳变电流统计分布测量和结I-V特征曲线测量的更为简便的方法, 实现约瑟夫森结各静态物理参数的推算. 当然, 实际应用时可能还需要通过其他方法来标定这些参数, 但本推算方法在器件按需设计和制备工艺相关参数设置中仍具有一定的参考意义.

2. 约瑟夫森结物理参数的推算方法

2.1 结电容和结临界电流的推算方法

电流偏置下的约瑟夫森结可以看作是一个质量为 $m = C\left[\varPhi_{0} / (2 \pi)\right]^{2}$ 的“粒子”在搓衣板势 $U(\varphi) =$ $-E_{\rm J}\left[\left(I_{\rm b} / I_{\rm c}\right) \varphi+\cos (\varphi)\right]$ 中运动的动力学模型^[12], 如图1所示. 其中, $I_{\rm b}$ 是外加偏置电流, C为结电容, $\varPhi_{0} = h/(2 e)$ 为磁通量子, $\varphi$ 为约瑟夫森结两端的位相差, 结的约瑟夫森能为 $E_{\rm J} = I_{\rm c}\varPhi_{0}/(2\pi)$ , $I_{\rm c}$ 为约瑟夫森结的临界电流.

$图 1 电流偏置下的约瑟夫森结势能曲线\r\nFig. 1. Potential of a current-biased Josephson junction.$

图 1 电流偏置下的约瑟夫森结势能曲线

Fig. 1. Potential of a current-biased Josephson junction.

下载: 全尺寸图片幻灯片

势阱中“粒子”的运动由方程

$m \ddot{\varphi}+\frac{1}{R C} \dot{\varphi}+\frac{{\rm d} U(\varphi)}{{\rm d} \varphi} = 0$

(1)

描述, 其中R为结电阻. 势阱的高度(即“粒子”跑出势阱的束缚需要克服的势垒高度)由下式给出^[13]:

$\Delta U = 2 E_{\rm {J}}\bigg[\sqrt {\left(1-\frac{I_{\rm {b}}}{I_{\rm {c}}}\right)}-\frac{I_{\rm{ b}}}{I_{\rm {c}}} \arccos (\varphi)\bigg].$

(2)

显然, 当外加偏置电流远小于约瑟夫森结临界电流时, “粒子”被很好地囚禁在势阱中. 这时, 结处于超导状态因而结两端的电压为零. 当偏置电流缓慢增加到接近结的临界电流 $I_{\rm c}$ 时, “粒子”可以通过热激发或者量子隧穿效应逃逸出势阱, 在宏观上就表现为约瑟夫森结两端可测量的电压信号, 这时所偏置的电流就称为跳变电流. 当然, 由于热激发效应和量子隧穿效应都是随机事件, 因此测量中记录下来的与电压信号对应的跳变电流也就是一种随机分布.

对应于“粒子”逃逸出或囚禁于势阱的两种不同状态. 假设“粒子”逃逸出势阱的分布函数由P 来描述, N(t)为粒子在势阱中的数目. 则:

$N(t) = 1-\int_{0}^{t} P\left(t^{\prime}\right) {\rm d} t^{\prime}.$

(3)

因此, 在 $\Delta t$ 时间内“粒子”从势阱中逃逸的数目为

${\rm d} N = -N(t)\varGamma(t){\rm d}t = -N\varGamma(I)\left(\frac{{\rm d} I}{{\rm d} t}\right)^{-1} {\rm d} I,$

(4)

其中, $\varGamma(t)$ 为“粒子”的逃逸概率, I为实验上可测量的跳变电流大小. 在 $N(0) = 1$ 的初始条件下, (4)式积分得: $N(I) = \exp \Big[-({\rm d} I / {\rm d} t)^{-1} \displaystyle\int_{0}^{I} \varGamma(I^{\prime}){\rm d} I^{\prime}\Big]$ . 可得跳变电流的分布函数为^[14]

$\begin{split} P(I) = \;&-\frac{{\rm d} N}{{\rm d} I} = \varGamma(I)\left(\frac{{\rm d} I}{{\rm d} t}\right)^{-1} \\ &\times\exp \left[-\left(\frac{{\rm d} I}{{\rm d} t}\right)^{-1} \int_{0}^{I} \varGamma(I^{\prime}) {\rm d} I^{\prime}\right]. \end{split}$

(5)

当温度高于某个临界温度 $T^*$ 时, 热激发占主导地位, “粒子”逃逸的速度由Kramers公式^[15]给出:

$\varGamma_{t} = \frac{\omega_{\rm p}}{2 \pi} a_{\rm t} \exp \left(-\frac{\Delta U}{k_{\rm B} T}\right),$

(6)

其中, $a_{\rm t}$ 是阻尼系数, 其取值范围为 $0 < a_{\rm t}\leqslant 1$ , 一般将 $a_{\rm t}$ 近似为 $1$ . (6)式中, $\omega_{\rm p} = \omega_{\rm p 0}\big[1-\left(I_{\rm b} / I_{\rm c}\right)^{2}\big]^{1 / 4},$ $\omega_{\rm p 0} = 2 \pi I_{\rm c} / (\varPhi_{0} C)$ . 当温度小于 $T^*$ 时, “粒子”逃逸出势阱主要是由量子隧穿效应导致, 其速度可通过WKB方法计算^[16]:

$\varGamma_{\rm q} = \frac{\omega_{\rm p}}{2 \pi} a_{\rm q} \exp \left[-\frac{7.2 \Delta U}{h \omega_{\rm p}}\left(1+\frac{0.87}{Q}\right)\right].$

(7)

它几乎与温度无关. 其中 $a_{\rm q} \!\approx\![120 \pi(7.2 \Delta U / \hbar \omega_{\rm p})]^{1 / 2}$ . 以上的理论结果被广泛应用于早期的约瑟夫森电子学特性研究, 下面将其应用于推算结的物理参数, 即实验测量跳变电流的统计分布, 并通过与理论模型的拟合, 来推算约瑟夫森结的临界电流 $I_{\rm c}$ 和电容参数C.

2.2 结阻尼参数的推算方法

直流偏置下的约瑟夫森结RCSJ模型由如下方程描述^[17]:

$I_{\rm b} = I_{\rm c} \sin (\varphi)+\frac{V}{R}+C \frac{{\rm d} V}{{\rm d} t},$

(8)

其中 $I_{\rm b}$ 为结的直流偏置, $V = [\hbar/( 2 e)] {\rm d} \varphi / {\rm d} t$ 为结两端的电压. 令 $\tau = \omega_{\rm c} t$ , 其中 $\omega_{\rm c} = \sqrt{2 e I_{\rm c} / (\hbar C)}$ , 则方程(8)可化为如下形式:

$i_{\rm b} = \frac{{\rm d}^{2} \varphi}{{\rm d} \tau^{2}}+\frac{1}{\sqrt{\beta_{\rm c}}} \frac{{\rm d} \varphi}{{\rm d} \tau}+\sin \varphi,$

(9)

这里, $\beta_{\rm c} = 2 e I_{\rm c} C R^{2} / \hbar$ 为阻尼参数, $i_{\rm b} = I_{\rm b} / I_{\rm c}$ . 因此, 对 $i_{\rm b} > 1$ 的非超导态, $i_{\rm b}$ 应该是 ${{\rm d} \varphi}/{{\rm d} \tau}$ 的一个线性函数. 不失一般性, 假设:

$i_{\rm b} = K\left(\frac{{\rm d} \varphi}{{\rm d} \tau}\right)+i_{\rm b0},$

(10)

这里, $i_{\rm b0}$ 为常数, 斜率K可由实验测得的结I-V曲线提取. 由此, 通过与(9)式数值解的比较即可建立阻尼参数 $\beta_{\rm c}$ 与斜率K的联系, 就可实现阻尼参数 $\beta_{\rm c}$ 的推算, 进而计算出低温条件下的结电阻R.

3. 约瑟夫森结的制备及其物理参数推算

3.1 约瑟夫森结的制备

采用悬空掩膜斜角度蒸发镀膜的工艺来制备SIS (superconductor-insulator-superconductor)约瑟夫森结^[18]. 所采用的超导材料为Al, 对Al氧化生成的氧化铝做势垒层, 制备的主要步骤分为匀胶、光刻、显影、斜角度蒸发镀膜和去胶等. 在清洗好的硅片上先后匀上LOR10B负胶和S1805正胶, 匀胶后分别进行烘胶, 溶解掉其中的有机溶剂, 使其与硅片贴牢固. 用深紫外曝光机和设计好的掩模版对其进行曝光; 曝光后的硅片使用显影液进行显影, 便可得到所需要的悬空掩膜结构. 进而, 将其放入电子束蒸发镀膜仪完成斜角度镀膜、氧化、和再一次斜角度镀膜三个步骤, 制备出约瑟夫森结样品. 整个结样品的制备流程如图2所示.

$图 2 约瑟夫森结制备流程图\r\nFig. 2. Preparation process of a Josephson junction.$

图 2 约瑟夫森结制备流程图

Fig. 2. Preparation process of a Josephson junction.

下载: 全尺寸图片幻灯片

最后将剩余的正、负光刻胶分别使用丙酮和N-甲基吡咯烷酮进行去除, 就可完成图3所示的约瑟夫森结样品的制备.

$图 3 约瑟夫森结样品\r\nFig. 3. Josephson junction sample.$

图 3 约瑟夫森结样品

Fig. 3. Josephson junction sample.

下载: 全尺寸图片幻灯片

3.2 约瑟夫森结I-V曲线和跳变电流的测量

约瑟夫森结I-V曲线一般采用四端子法进行测量, 一条线路通过一个1000倍的放大器测量结两端的电压信号, 另一条线路测流过结的电流信号. 任意波形发生器将电压信号通过可调电阻加到结的两端, 实现结的电流偏置. 这里, 可调电阻阻值远大于约瑟夫森结电阻, 以逐渐增大结的低频偏置电流. 流过结的电流和结两端可能产生的电压信号由数据采集卡采集. 将数据采集卡采集的数据在PC端进行处理, 便可得到图4所示的过阻尼约瑟夫森结标准I-V特性曲线.

$图 4 约瑟夫森结I-V特性曲线\r\nFig. 4. Measured I-V characteristic curve of the fabricated Josephson junction.$

图 4 约瑟夫森结I-V特性曲线

Fig. 4. Measured I-V characteristic curve of the fabricated Josephson junction.

下载: 全尺寸图片幻灯片

下面介绍在50 mK的极低环境温度下对约瑟夫森结跳变电流的实验测量. 为此, 需要给约瑟夫森结施加一个大小可调的直流偏置(直流值以 ${\rm d}I /{\rm d}t$ 的速率逐渐增大). 在0时刻触发同步信号, 开始计时; 随后调节偏置电流使之逐渐增大, 同时监测结两端的电压; 当结电压从零跳变到一阈值时, 计时结束, 记录此时的偏置电流大小作为跳变电流. 测量信号的时序设置如图5所示^[19].

$图 5 约瑟夫森结跳变电流实验测量的时序图\r\nFig. 5. Time sequence diagram for the junction jump current measurements.$

图 5 约瑟夫森结跳变电流实验测量的时序图

Fig. 5. Time sequence diagram for the junction jump current measurements.

下载: 全尺寸图片幻灯片

重复以上测量过程 $10^4$ 次, 并对采集的测量数据在PC端进行数据处理, 得到如图6所示的跳变电流大小测量值的统计分布图, M为测得对应跳变电流值的次数. 从图6中的数据可以看出, 在偏置电流小于 $4.0\times 10^{-7}$ A时, 没有测到任何跳变事件(即没有非零的电压信号); 在偏置电流大于 $4.8\times 10^{-7}$ A时, 结两端的电压几乎不再发生跳变, 而是逐渐增大. 因此, 可以预计所测量结的临界电流应该大于 $4.8\times 10^{-7}$ A.

$图 6 约瑟夫森结跳变电流及其次数统计\r\nFig. 6. Josephson junction jump currents and their statistical distributions.$

图 6 约瑟夫森结跳变电流及其次数统计

Fig. 6. Josephson junction jump currents and their statistical distributions.

下载: 全尺寸图片幻灯片

3.3 约瑟夫森结物理参数的推算

首先, 根据测量得到的跳变电流统计分布特性, 来推算结的临界电流 $I_{\rm c}$ 和结电容C. 对实验测得的结跳变电流的分布做归一化处理, 得到图7中的点状线跳变电流的统计分布. 利用跳变电流的统计分布公式(5)式, 通过调整约瑟夫森结临界电流 $I_{\rm c}$ 和电容C的预定值, 得到如图7所示的实验测量数据拟合. 在拟合实验数据时, 理论模型的参数设定为 $I_{\rm c} = 5.56\times10^{-7}$ A, $C = 23.3$ fF. 由于跳变电流统计分布理论模型中对结电阻的取值很不敏感, 因此在理论拟合中暂时取为其室温测量值 $R = 598$ $\Omega$ .

$图 7 结跳变电流的归一化统计分布: 理论拟合(红实线)与实验数据(点状线)\r\nFig. 7. Statistical distributions of the junction jump currents: theoretical simulations (red soild line), and measurement data (dot line).$

图 7 结跳变电流的归一化统计分布: 理论拟合(红实线)与实验数据(点状线)

Fig. 7. Statistical distributions of the junction jump currents: theoretical simulations (red soild line), and measurement data (dot line).

下载: 全尺寸图片幻灯片

为证实理论拟合中预定的参数是最优的, 具体分析了临近参数值对实验数据拟合的偏差程度, 如表1所列. 综合表1的数据可见, 图7中实验数据的理论拟合是最优的, 因此所测量约瑟夫森结的临界电流可推算为 $I_{\rm c} = 5.56\times10^{-7}$ A, 结电容可估计为 $C = 23.3$ fF.

表 1 不同参数取值对实验数据拟合的偏差度分析

Table 1. Deviations from the data simulated by using the different theoretical parameters.

$I_{\rm c}/(10^{-7}{\rm A})$	C/fF	峰值位置偏差(误差/实验值)	峰值高度偏差(误差/实验值)	半高宽偏差(误差/实验值)
$5.56$	$23.3$	0.06%	9.75%	3.52%
$5.57$	$23.3$	0.27%	9.75%	3.52%
$5.55$	$23.3$	0.15%	9.75%	3.52%
$5.56$	$23.4$	0.17%	9.74%	3.53%
$5.56$	$23.2$	0.31%	9.76%	3.52%

下载: 导出CSV

| 显示表格

下面推算结的阻尼参数值和电阻值. 由于约瑟夫森结的跳变电流分布的理论模型对电阻的变化很不敏感, 所以以上拟合时所取的室温结电阻测量值与极低温工作环境下的结电阻的实际数值有所偏离. 为此, 采用与归一化RCSJ方程(9)的数值解进行拟合的方法, 数值解求解采用四阶龙格库塔算法. 先实现约瑟夫森结阻尼参数 $\beta_{\rm c}$ 的推算, 进而计算出极低温工作时结电阻R的取值. 为此, 将实验测得的结I-V曲线(图4)中的I变量值除以 $I_{\rm c}$ 、电压变量除以 $\hbar\omega_{\rm c}/(2 e)$ 得到图8中黑色点状线所示的等效I-V曲线; 进而, 对非超导态工作区的结I-V特征曲线进行线性拟合, 其均方根误差RMSE (root mean squared error)为0.197. 得到如图8中红实线所示的拟合直线. 由此, 可提取方程(10)中的斜率参数为 $K = 0.295$ .

$图 8 约瑟夫森结的等效I-V曲线\r\nFig. 8. Effective I-V curve of the measured Josephson junction$

图 8 约瑟夫森结的等效I-V曲线

Fig. 8. Effective I-V curve of the measured Josephson junction

下载: 全尺寸图片幻灯片

接下来, 在非超导态区域(即取 $i_{\rm b} > 1$ )对方程(9)采用四阶龙格库塔算法进行数值求解, 得出阻尼参数 $\beta_{\rm c}$ 和斜率参数K的数值解, 并进行拟合. 图9中, 点状线表示由方程(9)数值解所得的阻尼参数 $\beta_{\rm c}$ 和斜率参数K的对应关系, 红实线是利用以下公式对点状线变化规律的拟合:

$图 9 斜率K随阻尼参数变化$\beta_{\rm c}$关系\r\nFig. 9. Relationship between the parameters $\beta_{\rm c}$ and K.$

图 9 斜率K随阻尼参数变化

$\beta_{\rm c}$

关系

Fig. 9. Relationship between the parameters

$\beta_{\rm c}$

and K.

下载: 全尺寸图片幻灯片

$K = 0.9679 \times \beta_{\rm c}^{-0.4721},$

(11)

均方根RMSE为0.006. 因此, 根据以上对所测试约瑟夫森结的K参数推算值 $K = 0.295$ , 再根据(11)式便可推算出该约瑟夫森结的阻尼系数为 $\beta_{\rm c} = 12.3$ , 由此说明所制备的约瑟夫森结是一个过阻尼结, 进而可推算极低温工作环境下结的正常态电阻值为 $R = 558.99$ $\Omega$ .

至此完成了对所制备的约瑟夫森结样品的物理参数推算. 为了检验根据低频电流偏置下进行实验测量所推算的约瑟夫森结物理参数的准确度, 将根据实验推算所得到的参数代入计算得到 $I_{\rm c}R =$ $5.56\times 10^{-7}\times 558.99 = 3.108\times 10^{-4}$ ${\rm{A}}\cdot\Omega$ , 与由著名的Ambgaokar-Baratoff公式^[20]

$I_{\rm c}R = \frac{\pi}{2e}\varDelta\tanh\left(\frac{\varDelta}{2k_{\rm B}T}\right)$

(12)

在 $T = 50$ mK工作温度下代入Al的超导能隙值 $\varDelta = 2.906\times 10^{-23} \;{\rm{J}}$ 计算所得到的值 $I^{\rm {AB}}_{\rm c}R_{\rm n} =$ $2.849\times$ $10^{-4}$ ${\rm{A}}\cdot \Omega$ 进行比较, 两者相差 $9.10\%$ . 需要说明的是, 由于实验测量是在远离超导转变温区进行的, 所以这里超导Al膜超导能隙 $\varDelta$ 可看成是不随温度变化的常数. 其中存在误差的主要原因是约瑟夫森结跳变电流的测量误差, 在与理论模型拟合时对临界电流和电容的推算误差导致了后续推算结果与超导能隙理论存在一定误差. 因此, 通过跳变电流分布实验测量所推算得到的 $I_{\rm c}, R$ 基本符合Ambgaokar-Baratoff公式, 这说明通过测量低频电流偏置下跳变电流的统计分布响应, 来实现约瑟夫森结物理参数推算的方法是可行的.

4. 结　论

约瑟夫森结作为超导量子计算芯片和微波单光子探测芯片的核心器件, 其物理参数很难直接测定但对器件的按需设计和制备却极其重要. 不同于通常的高频微波激励响应测量, 本文提出了一种基于低频电流偏置下跳变电流分布测量和I-V特性曲线的测量, 结合结RCSJ模型和跳变电流统计分布的理论模型, 实现了约瑟夫森结物理参数——结的临界电流 $I_{\rm c}$ 、结电容C、结电阻R 以及结阻尼参数 $\beta_{\rm c}$ 的推算方法. 这些根据实验数据推算出来的参数, 基本符合著名的Ambgaokar-Baratoff公式, 所以该实验推算方法是可行的. 本文所提出的方法, 与测量微波驱动下约瑟夫森结跳变电流统计分布的方案相比, 在实验上更加简便; 与采用回滞电流与临界电流比值相比, 不受I-V特性曲线原点上下涨落的影响, 具有独有的优势.

当然, 本文所实现的基于跳变电流分布实验测量的约瑟夫森结物理参数推算, 其精度有待于进一步提高. 首先, 跳变电流测量的样本数可以更多, 从而减少统计误差; 其次, 根据电压信号的跳变记录到的跳变电流, 通过精密地控制信号同步使其测量精度还有很大的提升空间, 等等. 不过, 受限于数值拟合的理论公式大多数也只是某种近似处理下的结果, 通过基于近似理论模型来拟合实验测量数据从中提取结的物理参数值, 本身也是影响参数推算的准确度. 所以, 约瑟夫森结的物理参数值还需要在器件后续应用(比如超导量子计算和微波单光子探测等)的研究中, 应用其他方法(如器件的能谱测量等)来进行进一步的测定. 本文的工作只是提供了可应用于约瑟夫森结器件设计和制备中, 如何实现物理参数实验推算的简单而有效方法, 为按需设计约瑟夫森结器件参数和设置器件制备工艺参数提供参考.

参考文献

[1]	Barends R, Kelly J, Megrant A, Veitia A, Sank D, Jeffrey E, White T C, Mutus J, Fowler A G, Campbell B, Chen Y, Chen Z, Chiaro B, Dunsworth A, Neill C, Malley P O, Roushan P, Vainsencher A, Wenner J, Korotkov A N, Cleland A N, Martinis J M 2014 Nature 508 500 Google Scholar
[2]	李春光, 王佳, 吴云, 王旭, 孙亮, 董慧, 高波, 李浩, 尤立星, 林志荣, 任浩, 李婧, 张文, 贺青, 王轶文, 韦联福, 孙汉聪, 王华兵, 李劲劲, 屈继峰 2021 物理学报 70 018501 Google Scholar Li C G, Wang J, Wu Y, Wang X, Sun L, Dong H, Gao B, Li H, You L X, Lin Z R, Ren H, Li J, Zhang W, He Q, Wang Y W, Wei L F, Sun H C, Wang H B, Li J J, Qu J F 2021 Acta Phys. Sin. 70 018501 Google Scholar
[3]	Arute F, Arya K, Babbush R, Bacon D, Bardin J, Barends R, Biswas R, Boixo S, Brandao F, Buell D, Burkett B, Chen Y, Chen Z J, Chiaro B, Collins R, Courtney W, Dunsworth A, Farhi E, Foxen B, Fowler A, Gidney C, Giustina M, Graff B, Guerin K, Habegger S, Harrigan M, Hartmann M, Ho A, Hoffmann M, Huang T, Humble T, Isakov S, Jeffrey E, Zhang J, Kafri D, Kechedzhi K, Kelly J, Klimov P, Knysh S, Korotkov A, Kostritsa F, Landhuis D, Lindmark M, Lucero E, Lyakh D, Mandrà S, McClean J, McEwen M, Megrant A, Mi X, Michielsen K, Mohseni M, Mutus J, Naaman O, Neeley M, Neill C, Niu M Y, Ostby E, Petukhov A, Platt J, Quintana C, Rieffel E, Roushan P, Rubin N, Sank D, Satzinger K, Smelyanskiy V, Sung K, Trevithick M, Vainsencher A, Villalonga B, Yao T J, Yeh P, Zalcman A, Neven H, Martinis J 2019 Nature 574 505 Google Scholar
[4]	Zhong H S, Wang H, Deng Y H, Chen M C, Peng L C, Luo Y H, Qin J, Wu D, Ding X, Hu Y, Hu P, Yang X Y, Zhang W J, Li H, Li Y X, Jiang X, Gan L, Yang G W, You L X, Wang Z, Li L, Liu N L, Lu C Y, Pan J W 2020 Science 370 1460
[5]	Sathyamoorthy S R, Stace T M, Johansson G 2016 C. R. Phys. 17 756 Google Scholar
[6]	Chen Y F, Hover D, Sendelbach S, Maurer L, Merkel S T, Pritchett E J, Wilhelm F K, McDermott R 2011 Phys. Rev. Lett. 370 1460
[7]	张裕恒, 李玉芝 2009 超导物理 (卷3) (合肥: 中国科学技术大学出版社) 第342−484页 Zhang Y H, Li Y Z 2009 Superconductor Physics (Vol. 3) (Hefei: China University of science and Technology Press) pp342−484 (in Chinese)
[8]	郑东宁 2021 物理学报 70 018502 Zheng D L 2021 Acta Phys. Sin. 70 018502
[9]	Clarke J, Wilhelm F K 2008 Nature 453 1031 Google Scholar
[10]	Sun G Z, Wang Y W, Cao J Y, Chen J, Ji Z M, Kang L, Xu W W, Yu Y, Han S Y, Wu P H 2008 Phys. Rev. B 77 104531 Google Scholar
[11]	Stewart W C 1968 Appl. Phys. Lett. 12 277 Google Scholar
[12]	Makhlin Y, Sch?n G, Shnirman A 2001 Rev. Mod. Phys. 73 357 Google Scholar
[13]	Stoutimore M J A, Rossolenko A N, Bolginov V V, Oboznov V A, Rusanov A Y, Baranov D S, Pugach N, Frolov S M, Ryazanov V V, Van Harlingen D J 2018 Phys. Rev. Lett. 121 17
[14]	崔大健, 林德华, 于海峰, 彭智慧, 朱晓波, 郑东宁, 景秀年, 吕力, 赵士平 2008 物理学报 09 5933 Google Scholar Cui D J, Lin D H, Yu H F, Peng Z H, Zhu X B, Zheng D N, Jing X N, Lu L, Zhao S P 2008 Acta Phys. Sin. 09 5933 Google Scholar
[15]	Kramers H A 1940 Physica 7 284 Google Scholar
[16]	Caldeira A O, Leggett A J 1981 Phys. Rev. Lett. 46 211 Google Scholar
[17]	Trees B R, Saranathan V, Stroud D 2005 Phys. Rev. E 71 016215 Google Scholar
[18]	陈伟 2019 博士学位论文 (南京: 南京大学) Chen W 2019 Ph. D. Dissertation (Nanjing: Nanjing University) (in Chinese)
[19]	曹俊宇, 孙国柱, 王轶文, 陈健, 吴培亨 2007 低温物理学报 03 196 Cao J Y, Sun G Z, Wang Y W, Chen J, Wu P H 2007 Chin. J. Low Temp. 2007 03 196 (in Chinese)
[20]	Ambegaokar V, Baratoff A 1963 Phys. Rev. Lett. 10 486 Google Scholar

施引文献

期刊类型引用(2)

1.	Aleksey G.Vostretsov，Svetlana G.Filatova. Algorithms for automatic measurement of SIS-type hysteretic underdamped Josephson junction's parameters by current-voltage characteristics. Journal of Electronic Science and Technology. 2023(04): 62-76 . 必应学术
2.	贺青，刘剑，韦联福. 微弱电磁信号的物理极限检测：单光子探测器及其研究进展. 广西师范大学学报(自然科学版). 2022(05): 1-23 . 百度学术

其他类型引用(1)

图 1 电流偏置下的约瑟夫森结势能曲线

Fig. 1. Potential of a current-biased Josephson junction.

下载: 全尺寸图片幻灯片

图 2 约瑟夫森结制备流程图

Fig. 2. Preparation process of a Josephson junction.

下载: 全尺寸图片幻灯片

图 3 约瑟夫森结样品

Fig. 3. Josephson junction sample.

下载: 全尺寸图片幻灯片

图 4 约瑟夫森结I-V特性曲线

Fig. 4. Measured I-V characteristic curve of the fabricated Josephson junction.

下载: 全尺寸图片幻灯片

图 5 约瑟夫森结跳变电流实验测量的时序图

Fig. 5. Time sequence diagram for the junction jump current measurements.

下载: 全尺寸图片幻灯片

图 6 约瑟夫森结跳变电流及其次数统计

Fig. 6. Josephson junction jump currents and their statistical distributions.

下载: 全尺寸图片幻灯片

图 7 结跳变电流的归一化统计分布: 理论拟合(红实线)与实验数据(点状线)

Fig. 7. Statistical distributions of the junction jump currents: theoretical simulations (red soild line), and measurement data (dot line).

下载: 全尺寸图片幻灯片

图 8 约瑟夫森结的等效I-V曲线

Fig. 8. Effective I-V curve of the measured Josephson junction

下载: 全尺寸图片幻灯片

图 9 斜率K随阻尼参数变化 $\beta_{\rm c}$ 关系

Fig. 9. Relationship between the parameters $\beta_{\rm c}$ and K.

下载: 全尺寸图片幻灯片

表 1 不同参数取值对实验数据拟合的偏差度分析

Table 1. Deviations from the data simulated by using the different theoretical parameters.

$I_{\rm c}/(10^{-7}{\rm A})$	C/fF	峰值位置偏差(误差/实验值)	峰值高度偏差(误差/实验值)	半高宽偏差(误差/实验值)
$5.56$	$23.3$	0.06%	9.75%	3.52%
$5.57$	$23.3$	0.27%	9.75%	3.52%
$5.55$	$23.3$	0.15%	9.75%	3.52%
$5.56$	$23.4$	0.17%	9.74%	3.53%
$5.56$	$23.2$	0.31%	9.76%	3.52%

下载: 导出CSV

[1]	Barends R, Kelly J, Megrant A, Veitia A, Sank D, Jeffrey E, White T C, Mutus J, Fowler A G, Campbell B, Chen Y, Chen Z, Chiaro B, Dunsworth A, Neill C, Malley P O, Roushan P, Vainsencher A, Wenner J, Korotkov A N, Cleland A N, Martinis J M 2014 Nature 508 500 Google Scholar
[2]	李春光, 王佳, 吴云, 王旭, 孙亮, 董慧, 高波, 李浩, 尤立星, 林志荣, 任浩, 李婧, 张文, 贺青, 王轶文, 韦联福, 孙汉聪, 王华兵, 李劲劲, 屈继峰 2021 物理学报 70 018501 Google Scholar Li C G, Wang J, Wu Y, Wang X, Sun L, Dong H, Gao B, Li H, You L X, Lin Z R, Ren H, Li J, Zhang W, He Q, Wang Y W, Wei L F, Sun H C, Wang H B, Li J J, Qu J F 2021 Acta Phys. Sin. 70 018501 Google Scholar
[3]	Arute F, Arya K, Babbush R, Bacon D, Bardin J, Barends R, Biswas R, Boixo S, Brandao F, Buell D, Burkett B, Chen Y, Chen Z J, Chiaro B, Collins R, Courtney W, Dunsworth A, Farhi E, Foxen B, Fowler A, Gidney C, Giustina M, Graff B, Guerin K, Habegger S, Harrigan M, Hartmann M, Ho A, Hoffmann M, Huang T, Humble T, Isakov S, Jeffrey E, Zhang J, Kafri D, Kechedzhi K, Kelly J, Klimov P, Knysh S, Korotkov A, Kostritsa F, Landhuis D, Lindmark M, Lucero E, Lyakh D, Mandrà S, McClean J, McEwen M, Megrant A, Mi X, Michielsen K, Mohseni M, Mutus J, Naaman O, Neeley M, Neill C, Niu M Y, Ostby E, Petukhov A, Platt J, Quintana C, Rieffel E, Roushan P, Rubin N, Sank D, Satzinger K, Smelyanskiy V, Sung K, Trevithick M, Vainsencher A, Villalonga B, Yao T J, Yeh P, Zalcman A, Neven H, Martinis J 2019 Nature 574 505 Google Scholar
[4]	Zhong H S, Wang H, Deng Y H, Chen M C, Peng L C, Luo Y H, Qin J, Wu D, Ding X, Hu Y, Hu P, Yang X Y, Zhang W J, Li H, Li Y X, Jiang X, Gan L, Yang G W, You L X, Wang Z, Li L, Liu N L, Lu C Y, Pan J W 2020 Science 370 1460
[5]	Sathyamoorthy S R, Stace T M, Johansson G 2016 C. R. Phys. 17 756 Google Scholar
[6]	Chen Y F, Hover D, Sendelbach S, Maurer L, Merkel S T, Pritchett E J, Wilhelm F K, McDermott R 2011 Phys. Rev. Lett. 370 1460
[7]	张裕恒, 李玉芝 2009 超导物理 (卷3) (合肥: 中国科学技术大学出版社) 第342−484页 Zhang Y H, Li Y Z 2009 Superconductor Physics (Vol. 3) (Hefei: China University of science and Technology Press) pp342−484 (in Chinese)
[8]	郑东宁 2021 物理学报 70 018502 Zheng D L 2021 Acta Phys. Sin. 70 018502
[9]	Clarke J, Wilhelm F K 2008 Nature 453 1031 Google Scholar
[10]	Sun G Z, Wang Y W, Cao J Y, Chen J, Ji Z M, Kang L, Xu W W, Yu Y, Han S Y, Wu P H 2008 Phys. Rev. B 77 104531 Google Scholar
[11]	Stewart W C 1968 Appl. Phys. Lett. 12 277 Google Scholar
[12]	Makhlin Y, Sch?n G, Shnirman A 2001 Rev. Mod. Phys. 73 357 Google Scholar
[13]	Stoutimore M J A, Rossolenko A N, Bolginov V V, Oboznov V A, Rusanov A Y, Baranov D S, Pugach N, Frolov S M, Ryazanov V V, Van Harlingen D J 2018 Phys. Rev. Lett. 121 17
[14]	崔大健, 林德华, 于海峰, 彭智慧, 朱晓波, 郑东宁, 景秀年, 吕力, 赵士平 2008 物理学报 09 5933 Google Scholar Cui D J, Lin D H, Yu H F, Peng Z H, Zhu X B, Zheng D N, Jing X N, Lu L, Zhao S P 2008 Acta Phys. Sin. 09 5933 Google Scholar
[15]	Kramers H A 1940 Physica 7 284 Google Scholar
[16]	Caldeira A O, Leggett A J 1981 Phys. Rev. Lett. 46 211 Google Scholar
[17]	Trees B R, Saranathan V, Stroud D 2005 Phys. Rev. E 71 016215 Google Scholar
[18]	陈伟 2019 博士学位论文 (南京: 南京大学) Chen W 2019 Ph. D. Dissertation (Nanjing: Nanjing University) (in Chinese)
[19]	曹俊宇, 孙国柱, 王轶文, 陈健, 吴培亨 2007 低温物理学报 03 196 Cao J Y, Sun G Z, Wang Y W, Chen J, Wu P H 2007 Chin. J. Low Temp. 2007 03 196 (in Chinese)
[20]	Ambegaokar V, Baratoff A 1963 Phys. Rev. Lett. 10 486 Google Scholar

[1]	张定, 朱玉莹, 汪恒, 薛其坤. 转角铜氧化物中的约瑟夫森效应. 物理学报, 2023, 72(23): 237402. doi: 10.7498/aps.72.20231815
[2]	李中祥, 王淑亚, 黄自强, 王晨, 穆清. 原子级控制的约瑟夫森结中Al₂O₃势垒层制备工艺. 物理学报, 2022, 71(21): 218102. doi: 10.7498/aps.71.20220820
[3]	陈恒杰, 薛航, 李邵雄, 王镇. 一种通过约瑟夫森结非线性频率响应确定微波耗散的方法. 物理学报, 2019, 68(11): 118501. doi: 10.7498/aps.68.20190167
[4]	王兰若, 钟源, 李劲劲, 屈继峰, 钟青, 曹文会, 王雪深, 周志强, 付凯, 石勇. 用于精密测量玻尔兹曼常数的量子电压噪声源芯片研制. 物理学报, 2018, 67(10): 108501. doi: 10.7498/aps.67.20172643
[5]	王松, 王星云, 周章渝, 杨发顺, 杨健, 傅兴华. 硼膜制备工艺、微观结构及其在硼化镁超导约瑟夫森结中的应用. 物理学报, 2016, 65(1): 017401. doi: 10.7498/aps.65.017401
[6]	陈钊, 何根芳, 张青雅, 刘建设, 李铁夫, 陈炜. 具有Washer型输入线圈的超导量子干涉放大器的制备与表征. 物理学报, 2015, 64(12): 128501. doi: 10.7498/aps.64.128501
[7]	曹文会, 李劲劲, 钟青, 郭小玮, 贺青, 迟宗涛. 用于电压基准的Nb/NbxSi1-x/Nb约瑟夫森单结的研制. 物理学报, 2012, 61(17): 170304. doi: 10.7498/aps.61.170304
[8]	张立森, 蔡理, 冯朝文. 线性延时反馈Josephson结的Hopf分岔和混沌化. 物理学报, 2011, 60(6): 060306. doi: 10.7498/aps.60.060306
[9]	张立森, 蔡理, 冯朝文. 约瑟夫森结中周期解及其稳定性的解析分析. 物理学报, 2011, 60(3): 030308. doi: 10.7498/aps.60.030308
[10]	岳宏卫, 阎少林, 周铁戈, 谢清连, 游峰, 王争, 何明, 赵新杰, 方兰, 杨扬, 王福音, 陶薇薇. 嵌入Fabry-Perot谐振腔的高温超导双晶约瑟夫森结的毫米波辐照特性研究. 物理学报, 2010, 59(2): 1282-1287. doi: 10.7498/aps.59.1282
[11]	岳宏卫, 王争, 樊彬, 宋凤斌, 游峰, 赵新杰, 何明, 方兰, 阎少林. 高温超导双晶约瑟夫森结阵列毫米波相干辐射. 物理学报, 2010, 59(8): 5755-5758. doi: 10.7498/aps.59.5755
[12]	雷佑铭, 徐伟. 一类约瑟夫森结混沌系统的谐和共振控制. 物理学报, 2008, 57(6): 3342-3352. doi: 10.7498/aps.57.3342
[13]	崔大健, 林德华, 于海峰, 彭智慧, 朱晓波, 郑东宁, 景秀年, 吕力, 赵士平. 本征约瑟夫森结跳变电流分布的量子修正. 物理学报, 2008, 57(9): 5933-5936. doi: 10.7498/aps.57.5933
[14]	吴炳国, 赵志刚, 尤育新, 刘楣. 二维约瑟夫森结阵列中的相变及噪声频谱研究. 物理学报, 2007, 56(3): 1680-1685. doi: 10.7498/aps.56.1680
[15]	周铁戈, 宋凤斌, 左涛, 顾静, 夏侯海, 胡雅婷, 赵新杰, 方兰, 阎少林. 本征约瑟夫森结阵列的PSpice模型及混沌行为研究. 物理学报, 2007, 56(11): 6307-6314. doi: 10.7498/aps.56.6307
[16]	李照鑫, 邹健, 蔡金芳, 邵彬. 电荷量子比特与量子化光场之间的纠缠. 物理学报, 2006, 55(4): 1580-1584. doi: 10.7498/aps.55.1580
[17]	肖宇飞, 王登龙, 王凤姣, 颜晓红. 非对称的玻色-爱因斯坦凝聚中的约瑟夫森结的动力学性质. 物理学报, 2006, 55(2): 547-550. doi: 10.7498/aps.55.547
[18]	王震宇, 廖红印, 周世平. 直流偏置的与RLC谐振器耦合的约瑟夫森结动力学行为的数值模拟. 物理学报, 2001, 50(10): 1996-2000. doi: 10.7498/aps.50.1996
[19]	姚希贤. 关于电阻分路约瑟夫森结的解析解. 物理学报, 1978, 27(5): 559-568. doi: 10.7498/aps.27.559
[20]	刘东, 曹淑兰, 桑亚男, 王慧芝, 刘福绥, 万玉珍. 约瑟夫森隧道结微波感应效应的实验和理论研究. 物理学报, 1976, 25(2): 97-104. doi: 10.7498/aps.25.97

期刊类型引用(2)

1.	Aleksey G.Vostretsov，Svetlana G.Filatova. Algorithms for automatic measurement of SIS-type hysteretic underdamped Josephson junction's parameters by current-voltage characteristics. Journal of Electronic Science and Technology. 2023(04): 62-76 . 必应学术
2.	贺青，刘剑，韦联福. 微弱电磁信号的物理极限检测：单光子探测器及其研究进展. 广西师范大学学报(自然科学版). 2022(05): 1-23 . 百度学术

其他类型引用(1)

计量

文章访问数: 8095
PDF下载量: 285
被引次数: 3

姓名
邮箱
手机号码
标题
留言内容
验证码

搜索

留言板