约瑟夫森结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 特性和高度的非线性等），约瑟夫森结得到广泛的研究和应用，并成为超导电子器件的核心部件。

实际使用中的约瑟夫森结总处于某一电路之中，因此，利用等效电路理论来研究和分析约瑟夫森结的物理行为是一种很有效的方法。

在各模型中，其物理行为均可用微分方程来描述，但这些方程大多不易直接求解析解，因而发展了很多间接解法]5[],4[。

其中，利用电路模拟（RCSJ 模型和RSJ 模型等等），如图1、图2所示，并用数值计算来研究约瑟夫森结的方法最直接，简易。

图1：RCSJ 模型等效电路 图2：RSJ 模型等效电路Resistively Capacitance Shunted Junction Resistively Shunted Junction2． 约瑟夫森效应及约瑟夫森结简介1962年，约瑟夫森（B. D. Josephson ）提出：两块用绝缘薄层隔开且紧密地接近的超导体间，甚至在没有电势差的情况下，电子仍能够穿过绝缘薄层（隧道现象）。

在不到一年的时间内，安德森（P. W. Anderson ）和罗韦尔（J. M. Rowell ）等人从实验上观察到这一现象，证实了约瑟夫森的预言，此即约瑟夫森效应]6[（Josephson Effection ），又称隧道效应。

两块超导体通过一绝缘薄层（厚度为10埃左右）连接起来的组合称S-I-S 超导隧道结或约瑟夫森结。

绝缘层对电子来说是一势垒，一块超导体中的电子能穿过势垒进入另一超导体中，这是约瑟夫森结中特有的量子隧道效应]7[。

绝缘层太厚时，隧道效应不明显，太薄时，两块超导体实际连成一体，这两种情况下均不易观察到约瑟夫森效应，只有当绝缘层不太厚也不太薄，即弱耦合时，才能观察到显著的约瑟夫森效应。

3． 电路模型及计算机模拟理想约瑟夫森结遵循如下方程：ϕsin c I I = （1）V e dt d2=ϕ （2） 式中，ϕ是结中心量子波函数的相位差,c I 是总临界电流,与外加磁场有关（文献3中有详细介绍）。

等式（1）、（2）所描述的仅仅是电子对携带的电流，从理论上讲，结间还有位移电流、准粒子隧道电流和绝缘体漏电流，研究时应该把这些电流都考虑进去。

研究约瑟夫森结时， 一般都把一个理想的约瑟夫森结与结电阻R 和结电容C 相并联，即用电容模拟位移电流，用电阻模拟准粒子隧道电流和绝缘体漏电流。

这样的模型称为电阻、电容分路结模型，简称RCSJ 模型（Resistively Capacitance Shunted Junction ） ，如图1所示。

对于小面积隧道结、超导微桥以及点接触结，当其结电容很小时，可以不考虑电容的影响。

于是可用电阻分路结模型（Resistively Shunted Junction ）来描述：一个理想的约瑟夫森结与一个电阻相并联，简称RSJ 模型，如图2所示。

采用RCSJ 模型，并同时加上直流电流源和交流电流源后，可用如下方程描述：dtdV C R V I t I I c a a d ++=+)sin()sin(ϕω （3）将（2）式代入上式并用c I 除，再令c d d I I i =，c a a I I i = ， e k 2=，t c ωτ=，CeI c c 2=ω 即可得到无量纲化后的方程： ϕτϕβτϕωτsin )sin(22++=+d d d d i i c a d （4） 在上式中，c a ωωω=，22CR eI c c =β，c β为阻尼参数，ω为归一化的外加交流电流圆频率。

上述无量纲化的非线性微分方程没有确切的解析解，因此，我们将利用Matlab 求其数值解。

求解中采用Matlab 自带的四阶龙格----库塔（Runge-Kutta ）算法]8[，此算法可以在保证高精度的前提下，快速求解微分方程的数值解，其Matlab 函数为ode45()。

理论上已证明，在RCSJ 模型中接入直流电源，当d I =c I 时，约瑟夫森结间电压0≠V ，且电容的存在会使电压的产生滞后]9[。

在RCSJ 模型中接入交流电源时，会得到台阶状的I-V 特性曲线，且台阶是按一定的规律]10[分布的。

用上述的方法建立数值模型，借助Matlab 进行计算、绘图，可以精确、直观地从数值模拟实验上证明和看到上述现象。

此外，RCSJ 的数值模型还是一个很好的非线性混沌现象研究工具，因为我们不仅可以在保证精度的前提下，随意的改变电路的各个参数来研究各参量之间的关系和系统进入混沌的过程，而且可以通过绘制相图、庞加莱截面]11[、时序图和功率谱等对混沌产生机理进行分析。

4． Matlab 模拟结果及分析a. 约瑟夫森结的直流I-V 特性由（2）式可知约瑟夫森结平均电压为：τϕωd d eU c /2 = （5） 因此，计算机模拟约瑟夫森结伏安特性时可用d i —dt d /ϕ曲线代替I-V 特性曲线。

计算时用函数mean()近似地计算电压平均值，可以得到比较精确的结果，绘出的I-V 特性曲线如图3、图4所示。

图3：约瑟夫森结直流I-V 特性曲线 图4：约瑟夫森结直流I-V 特性曲线Josephson Junction DC I-V characteristic curve Josephson Junction DC I-V characteristic curve分析：从图中可以看出，约瑟夫森结注入直流电流后（1=c I ），当电流d I < c I 时，U=0，此为约瑟夫森结的超导态电流；当d I =c I 时是最大临界电流；当d I > c I 时，U ≠0，此时约瑟夫森结中的电流由正常态电流和超导电流组成。

电容的存在并没有改变约瑟夫森结的性质，但会使得电压产生滞后，且随着电容的增大，滞后现象会增强，如图4所示。

b. 约瑟夫森结的交流I-V 特性图5：约瑟夫森结交流I-V 特性曲线Josephson Junction AC I-V characteristic curve当同时向约瑟夫森结注入直流电流d i 和交流电流a i ，交流的I-V 特性曲线上出现台阶现象，如图5所示。

数学计算的结果证明了这些台阶出现在：f en U n 2 = n=0，1，2，… （6） 对于无量纲化的结果来说，交流台阶出现在：ωτϕn d d n =〉〈 n=0，1，2，… （7） c. 约瑟夫森结的高度非线性特性从上述无量纲化过程中，可以知道归一化电压u 的表达式为：cI eC V d d u 2==τϕ （8） 因此，画分岔图、相图和庞加莱截面时纵轴就是归一化电压u 。

下面给出了通过Matlab 模拟计算得到的非线性混沌现象的分岔图、相图、庞加莱截面及功率谱，也给出了一些简要的分析方法。

i. 通向混沌的道路图6：系统u 随参数d i 的分岔图 图7：系统u 随参数c β的分岔图The bifurcation diagram of U with the parameter d i The bifurcation diagram of U with the parameter c β图8：系统u 随参数c β的分岔图 图9：系统u 随参数a i 的分岔图 The bifurcation diagram of U with the parameter c β The bifurcation diagram of U with the parameter a i通向混沌的道路有很多种]12[，如倍周期分岔、阵发性方式、茹厄勒——塔肯斯道路等。

它们是具有普适性的，很多动力学系统随着参数的改变会以倍周期分岔或阵发性的方式进入混沌状态，在混沌区中，由于阵发性的作用，又会出现一些具有周期性的窗口。

如图6，图7，图8，图9所示，它们分别是约瑟夫森结电压u 随参数d i 、c β、a i 的分岔图。

我们先简要分析图6和图7，然后在ii 中用相图、时序图、庞加莱截面和功率谱着重分析图9。

图6是以d i 为控制参数(0.8<d i <2.0)，其他参数为(c β=1，5.0=ω，65.0=a i )时，得到的u 随d i 变化的分岔图。

当0.8<d i <1.15时，系统处于1周期状态；当d i 等于1.15时，系统进入2周期状态，接着是4周期状态、6周期状态、4周期状态；当1.24<d i <1.43时，系统又处于2周期状态；当1.43<d i <1.6时，系统处于内部有周期窗口的混沌状态；当1.6<d i <1.68 时，系统是3周期状态；当1.68<d i <2.0时，系统处于混沌状态，且内部有周期窗口出现。

从图6我们看到了不完全级联的过程，即随直流电流的增加，系统由1周期解(单周期的运动)进入到2，4，6周期解(多周期的运动)，然后返回到4，2周期解，之后才进入混沌状态，而没有直接进入混沌，这称为不完全级联；在选择的参数范围内，系统随d i 的增加以阵发方式进入混沌状态，且内部有周期窗口。

图7，图8（局部放大图）是以c β为控制参数(0<c β<1.5)，固定其他参数(0=d i ，2.0=ω，50.1=a i )时，得到的u 随c β变化的分岔图。

从图中很容易看出：随着参数c β的减小，周期运动状态和混沌状态会不断交替出现，且混沌区中含有较狭窄的周期窗口，最后完全进入混沌状态；陆续出现周期1，周期3，周期5，周期7，…的周期状态区 ，且周期状态区域逐渐变窄，直到消失；在选择的参数范围内，系统随c β的减小以阵发性方式进入混沌状态；随着参数c β的减小，系统由周期状态进入到混沌状态的过程是级联过程。