江苏科技大学数理学院 约瑟夫森结I-V特性及非线性的数值模拟
约瑟夫森结I-V特性及非线性的数值模拟
彭加福
（江苏科技大学数理学院，应用物理，0640502112）
摘要：本文基于Matlab对约瑟夫森结（Josephson Junction）RCSJ模型的交直流I-V特性
及非线性混沌现象进行数值模拟。通过计算机数值模拟得到该模型的非线性微分方程数值解，研究了RCSJ模型中各参量对约瑟夫森结的影响，进而简要分析其I-V特性和非线性混沌现象的产生机理，绘制出约瑟夫森结的交直流I-V特性曲线、非线性微分方程的相图及因其高度非线性而引起的通过倍周期分岔和阵发性原理进入混沌状态的分岔图。
关键词：超导器件 隧道效应 约瑟夫森结 弱耦合 倍周期分岔 庞加莱截面 混沌 1． 引 言
自1911年荷兰科学家昂纳斯（H. K. Onnes）发现汞的超导现象以来，人们对超导进行了大量开拓性的研究，使超导理论[1]日趋成熟，与此同时，超导技术也在各个领域得到深入而广泛的应用[2]。约瑟夫森效应的发现开拓了超导量子干涉仪（SQUID）在弱电方面的应用。人们在对约瑟夫森效应进行研究的过程中发明了各种超导器件及应用电路[3]，促使超导技术应用的新领域——超导电子学逐渐发展起来。在其中，因具有各种独特性（量子干涉、特殊的I-V特性和高度的非线性等），约瑟夫森结得到广泛的研究和应用，并成为超导电子器件的核心部件。实际使用中的约瑟夫森结总处于某一电路之中，因此，利用等效电路理论来研究和分析约瑟夫森结的物理行为是一种很有效的方法。在各模型中，其物理行为均可用微分方程来描述，但这些方程大多不易直接求解析解，因而发展了很多间接解法[4],[5]。其中，利用电路模拟（RCSJ模型和RSJ模型等等），如图1、图2所示，并用数值计算来研究约瑟夫森结的方法最直接，简易。
图1：RCSJ模型等效电路 图2：RSJ模型等效电路
Resistively Capacitance Shunted Junction Resistively Shunted Junction
_______________________________________________________________________________________________作者：彭加福 0640502112 （jingyujiafu@163.com） 共16页 - 1 -
江苏科技大学数理学院 约瑟夫森结I-V特性及非线性的数值模拟
2． 约瑟夫森效应及约瑟夫森结简介
1962年，约瑟夫森（B. D. Josephson）提出：两块用绝缘薄层隔开且紧密地接近的超导体间，甚至在没有电势差的情况下，电子仍能够穿过绝缘薄层（隧道现象）。在不到一年的时间内，安德森（P. W. Anderson）和罗韦尔（J. M. Rowell）等人从实验上观察到这一现象，证实了约瑟夫森的预言，此即约瑟夫森效应[6]（Josephson Effection），又称隧道效应。
两块超导体通过一绝缘薄层（厚度为10埃左右）连接起来的组合称S-I-S超导隧道结或约瑟夫森结。绝缘层对电子来说是一势垒，一块超导体中的电子能穿过势垒进入另一超导体中，这是约瑟夫森结中特有的量子隧道效应[7]。绝缘层太厚时，隧道效应不明显，太薄时，两块超导体实际连成一体，这两种情况下均不易观察到约瑟夫森效应，只有当绝缘层不太厚也不太薄，即弱耦合时，才能观察到显著的约瑟夫森效应。
3． 电路模型及计算机模拟
理想约瑟夫森结遵循如下方程：
I?Icsin? （1）
d?2e?V （2） dt?式中，?是结中心量子波函数的相位差,Ic是总临界电流,与外加磁场有关（文献3中有详细介绍）。
等式（1）、（2）所描述的仅仅是电子对携带的电流，从理论上讲，结间还有位移电流、准粒子隧道电流和绝缘体漏电流，研究时应该把这些电流都考虑进去。研究约瑟夫森结时， 一般都把一个理想的约瑟夫森结与结电阻R和结电容C相并联，即用电容模拟位移电流，用电阻模拟准粒子隧道电流和绝缘体漏电流。这样的模型称为电阻、电容分路结模型，简称RCSJ模型（Resistively Capacitance Shunted Junction） ，如图1所示。对于小面积隧道结、超导微桥以及点接触结，当其结电容很小时，可以不考虑电容的影响。于是可用电阻分路结模型（Resistively Shunted Junction）来描述：一个理想的约瑟夫森结与一个电阻相并联，简称RSJ模型，如图2所示。
采用RCSJ模型，并同时加上直流电流源和交流电流源后，可用如下方程描述：
Id?Iasin(?at)?Icsin(?)?VdV?C （3） Rdt_______________________________________________________________________________________________
作者：彭加福 0640502112 （jingyujiafu@163.com） 共16页 - 2 -
江苏科技大学数理学院 约瑟夫森结I-V特性及非线性的数值模拟
将（2）式代入上式并用Ic除，再令id?即可得到无量纲化后的方程：
IdI2e?c?k?，，ia?a ，???ct，
?IcIc2eIc ?Cd2?d?id?iasin(??)????sin? （4） cd?d?2在上式中，??圆频率。
上述无量纲化的非线性微分方程没有确切的解析解，因此，我们将利用Matlab求其数值解。求解中采用Matlab自带的四阶龙格----库塔（Runge-Kutta）算法
[8]?a?，?c?，?c为阻尼参数，?为归一化的外加交流电流
2eIcCR2?c，此算法可以在保
证高精度的前提下，快速求解微分方程的数值解，其Matlab函数为ode45()。
理论上已证明，在RCSJ模型中接入直流电源，当Id=Ic时，约瑟夫森结间电压V?0，且电容的存在会使电压的产生滞后
[9]。在RCSJ模型中接入交流电源时，会得到台阶状的I-V
特性曲线，且台阶是按一定的规律[10]分布的。用上述的方法建立数值模型，借助Matlab进行计算、绘图，可以精确、直观地从数值模拟实验上证明和看到上述现象。此外，RCSJ的数值模型还是一个很好的非线性混沌现象研究工具，因为我们不仅可以在保证精度的前提下，随意的改变电路的各个参数来研究各参量之间的关系和系统进入混沌的过程，而且可以通过绘制相图、庞加莱截面
[11]
、时序图和功率谱等对混沌产生机理进行分析。
4． Matlab模拟结果及分析
a. 约瑟夫森结的直流I-V特性
由（2）式可知约瑟夫森结平均电压为：
U???cd?/d? （5） 2e因此，计算机模拟约瑟夫森结伏安特性时可用id—d?/dt曲线代替I-V特性曲线。计算时用函数mean()近似地计算电压平均值，可以得到比较精确的结果，绘出的I-V特性曲线如图3、图4所示。
_______________________________________________________________________________________________作者：彭加福 0640502112 （jingyujiafu@163.com） 共16页 - 3 -
约瑟夫森结I-V特性及非线性Matlab模拟
