约瑟夫森结i-v特性及非线性matlab模拟.doc

江苏科技大学数理学院 约瑟夫森结 I-V 特性及非线性的数值模拟_______________________________________________________________________________________________作者：彭加福 0640502112 （jingyujiafu@ ） 共 16 页 - 1 -约瑟夫森结 I-V 特性及非线性的数值模拟彭加福（江苏科技大学数理学院，应用物理，0640502112）摘要：本文基于 Matlab 对约瑟夫森结（Josephson Junction）RCSJ 模型的交直流 I-V 特性及非线性混沌现象进行数值模拟通过计算机数值模拟得到该模型的非线性微分方程数值解，研究了 RCSJ 模型中各参量对约瑟夫森结的影响，进而简要分析其 I-V 特性和非线性混沌现象的产生机理，绘制出约瑟夫森结的交直流 I-V 特性曲线、非线性微分方程的相图及因其高度非线性而引起的通过倍周期分岔和阵发性原理进入混沌状态的分岔图。

关键词：超导器件 隧道效应 约瑟夫森结 弱耦合 倍周期分岔 庞加莱截面 混沌 1． 引 言自 1911 年荷兰科学家昂纳斯（H. K. Onnes）发现汞的超导现象以来，人们对超导进行了大量开拓性的研究，使超导理论 日趋成熟，与此同时，超导技术也在各个领域得到深入]1[而广泛的应用 约瑟夫森效应的发现开拓了超导量子干涉仪（SQUID）在弱电方面的应用]2[人们在对约瑟夫森效应进行研究的过程中发明了各种超导器件及应用电路 ，促使超导技]3[术应用的新领域——超导电子学逐渐发展起来在其中，因具有各种独特性（量子干涉、特殊的 I-V 特性和高度的非线性等） ，约瑟夫森结得到广泛的研究和应用，并成为超导电子器件的核心部件实际使用中的约瑟夫森结总处于某一电路之中，因此，利用等效电路理论来研究和分析约瑟夫森结的物理行为是一种很有效的方法在各模型中，其物理行为均可用微分方程来描述，但这些方程大多不易直接求解析解，因而发展了很多间接解法 其中，]5[,4利用电路模拟（RCSJ 模型和 RSJ 模型等等），如图 1、图 2 所示，并用数值计算来研究约瑟夫森结的方法最直接，简易图 1： RCSJ 模型等效电路 图 2： RSJ 模型等效电路Resistively Capacitance Shunted Junction Resistively Shunted Junction 江苏科技大学数理学院 约瑟夫森结 I-V 特性及非线性的数值模拟_______________________________________________________________________________________________作者：彭加福 0640502112 （jingyujiafu@ ） 共 16 页 - 2 -2．约瑟夫森效应及约瑟夫森结简介1962 年，约瑟夫森（B. D. Josephson）提出：两块用绝缘薄层隔开且紧密地接近的超导体间，甚至在没有电势差的情况下，电子仍能够穿过绝缘薄层（隧道现象） 。

在不到一年的时间内，安德森（P. W. Anderson）和罗韦尔（J. M. Rowell）等人从实验上观察到这一现象，证实了约 瑟 夫 森 的 预 言 ， 此 即 约瑟夫森效应 （ Josephson Effection） ，又称隧道]6[效应两块超导体通过一绝缘薄层（厚度为 10 埃左右）连接起来的组合称 S-I-S 超导隧道结或约瑟夫森结绝缘层对电子来说是一势垒，一块超导体中的电子能穿过势垒进入另一超导体中，这是约瑟夫森结中特有的量子隧道效应 绝缘层太厚时，隧道效应不明显，太薄时，]7[两块超导体实际连成一体，这两种情况下均不易观察到约瑟夫森效应，只有当 绝缘层不太厚也不太薄，即弱耦合时，才能观察到显著的约瑟夫森效应3．电路模型及计算机模拟理想约瑟夫森结遵循如下方程：（1）sincIVedth2（2）式中， 是结中心量子波函数的相位差, 是总临界电流,与外加磁场有关（文献3中有cI详细介绍）等式（1）、（2）所描述的仅仅是电子对携带的电流，从理论上讲，结间还有位移电流、准粒子隧道电流和绝缘体漏电流，研究时应该把这些电流都考虑进去研究约瑟夫森结时， 一般都把一个理想的约瑟夫森结与结电阻R和结电容C相并联，即用电容模拟位移电流，用电阻模拟准粒子隧道电流和绝缘体漏电流。

这样的模型称为电阻、电容分路结模型，简称RCSJ模型（Resistively Capacitance Shunted Junction） ，如图1所示对于小面积隧道结、超导微桥以及点接触结，当其结电容很小时，可以不考虑电容的影响于是可用电阻分路结模型（Resistively Shunted Junction）来描述：一个理想的约瑟夫森结与一个电阻相并联，简称RSJ模型，如图2所示采用RCSJ模型，并同时加上直流电流源和交流电流源后，可用如下方程描述：江苏科技大学数理学院 约瑟夫森结 I-V 特性及非线性的数值模拟_______________________________________________________________________________________________作者：彭加福 0640502112 （jingyujiafu@ ） 共 16 页 - 3 -dtVCRItIcad )sin()si(（3）将（2）式代入上式并用 除，再令 ， ， ， ，cIcdIicaIihek2tc即可得到无量纲化后的方程：CeIcchsin)sin(2di cad（4）在上式中， ， ， 为阻尼参数， 为归一化的外加交流电ca2CReIchc流圆频率。

上述无量纲化的非线性微分方程没有确切的解析解，因此，我们将利用Matlab求其数值解求解中采用Matlab自带的四阶龙格----库塔（Runge-Kutta ）算法 ，此算法可以在]8[保证高精度的前提下，快速求解微分方程的数值解，其Matlab 函数为ode45()理论上已证明，在RCSJ模型中接入直流电源，当 = 时，约瑟夫森结间电压 ，dIc 0V且电容的存在会使电压的产生滞后 在RCSJ模型中接入交流电源时，会得到台阶状的 I-V]9[特性曲线，且台阶是按一定的规律 分布的用上述的方法建立数值模型，借助 Matlab进]10[行计算、绘图，可以精确、直观地从数值模拟实验上证明和看到上述现象此外，RCSJ的数值模型还是一个很好的非线性混沌现象研究工具，因为我们不仅可以在保证精度的前提下，随意的改变电路的各个参数来研究各参量之间的关系和系统进入混沌的过程，而且可以通过绘制相图、庞加莱截面 、时序图和功率谱等对混沌产生机理进行分析]1[4．Matlab 模拟结果及分析a. 约瑟夫森结的直流 I-V 特性由（2）式可知约瑟夫森结平均电压为：江苏科技大学数理学院 约瑟夫森结 I-V 特性及非线性的数值模拟_______________________________________________________________________________________________作者：彭加福 0640502112 （jingyujiafu@ ） 共 16 页 - 4 -deUc/2h（5）因此，计算机模拟约瑟夫森结伏安特性时可用 — 曲线代替 I-V 特性曲线。

计dit/算时用函数 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 分析：从图中可以看出，约瑟夫森结注入直流电流后（ ） ，当电流 时，U 0，dIc c此时约瑟夫森结中的电流由正常态电流和超导电流组成电容的存在并没有改变约瑟夫森结的性质，但会使得电压产生滞后，且随着电容的增大，滞后现象会增强，如图 4 所示b. 约瑟夫森结的交流 I-V 特性图 5： 约瑟夫森结交流 I-V 特性曲线 Josephson Junction AC I-V characteristic curve江苏科技大学数理学院 约瑟夫森结 I-V 特性及非线性的数值模拟_______________________________________________________________________________________________作者：彭加福 0640502112 （jingyujiafu@ ） 共 16 页 - 5 -当同时向约瑟夫森结注入直流电流 和交流电流 ，交流的 I-V 特性曲线上出现台阶现diai象，如图 5 所示。

数学计算的结果证明了这些台阶出现在：n=0，1，2，… fenU2h（6）对于无量纲化的结果来说，交流台阶出现在：n=0，1，2，… nd（7）c. 约瑟夫森结的高度非线性特性从上述无量纲化过程中，可以知道归一化电压 u 的表达式为：cIeCVduh2（8）因此，画分岔图、相图和庞加莱截面时纵轴就是归一化电压 u下面给出了通过 Matlab模拟计算得到的非线性混沌现象的分岔图、相图、庞加莱截面及功率谱，也给出了一些简要的分析方法i. 通向混沌的道路图 6： 系统 u 随参数 的分岔图 图 7： 系统 u 随参数 的分岔图di cThe bifurcation diagram of U with the parameter The bifurcation diagram of U with the parameter di c江苏科技大学数理学院 约瑟夫森结 I-V 特性及非线性的数值模拟_______________________________________________________________________________________________作者：彭加福 0640502112 （jingyujiafu@ ） 共 16 页 - 6 -图 8： 系统 u 随参数 。