约瑟夫森结I-V特性及非线性Matlab模拟Word文档下载推荐.doc

1、在其中，因具有各种独特性（量子干涉、特殊的I-V特性和高度的非线性等），约瑟夫森结得到广泛的研究和应用，并成为超导电子器件的核心部件。实际使用中的约瑟夫森结总处于某一电路之中，因此，利用等效电路理论来研究和分析约瑟夫森结的物理行为是一种很有效的方法。在各模型中，其物理行为均可用微分方程来描述，但这些方程大多不易直接求解析解，因而发展了很多间接解法。其中，利用电路模拟（RCSJ模型和RSJ模型等等），如图1、图2所示，并用数值计算来研究约瑟夫森结的方法最直接，简易。图1：RCSJ模型等效电路 图2：RSJ模型等效电路Resistively Capacitance Shunted Junction

2、 Resistively Shunted Junction 2 约瑟夫森效应及约瑟夫森结简介1962年，约瑟夫森（B. D. Josephson）提出：两块用绝缘薄层隔开且紧密地接近的超导体间，甚至在没有电势差的情况下，电子仍能够穿过绝缘薄层（隧道现象）。在不到一年的时间内，安德森（P. W. Anderson）和罗韦尔（J. M. Rowell）等人从实验上观察到这一现象，证实了约瑟夫森的预言，此即约瑟夫森效应（Josephson Effection），又称隧道效应。两块超导体通过一绝缘薄层（厚度为10埃左右）连接起来的组合称S-I-S超导隧道结或约瑟夫森结。绝缘层对电子来说是一势垒，一块超

3、导体中的电子能穿过势垒进入另一超导体中，这是约瑟夫森结中特有的量子隧道效应。绝缘层太厚时，隧道效应不明显，太薄时，两块超导体实际连成一体，这两种情况下均不易观察到约瑟夫森效应，只有当绝缘层不太厚也不太薄，即弱耦合时，才能观察到显著的约瑟夫森效应。3 电路模型及计算机模拟理想约瑟夫森结遵循如下方程： （1） （2）式中，是结中心量子波函数的相位差,是总临界电流,与外加磁场有关（文献3中有详细介绍）。等式（1）、（2）所描述的仅仅是电子对携带的电流，从理论上讲，结间还有位移电流、准粒子隧道电流和绝缘体漏电流，研究时应该把这些电流都考虑进去。研究约瑟夫森结时， 一般都把一个理想的约瑟夫森结与结电阻R

4、和结电容C相并联，即用电容模拟位移电流，用电阻模拟准粒子隧道电流和绝缘体漏电流。这样的模型称为电阻、电容分路结模型，简称RCSJ模型（Resistively Capacitance Shunted Junction） ，如图1所示。对于小面积隧道结、超导微桥以及点接触结，当其结电容很小时，可以不考虑电容的影响。于是可用电阻分路结模型（Resistively Shunted Junction）来描述：一个理想的约瑟夫森结与一个电阻相并联，简称RSJ模型，如图2所示。采用RCSJ模型，并同时加上直流电流源和交流电流源后，可用如下方程描述： （3）将（2）式代入上式并用除，再令， ， 即可得到无量纲

5、化后的方程： （4）在上式中，为阻尼参数，为归一化的外加交流电流圆频率。上述无量纲化的非线性微分方程没有确切的解析解，因此，我们将利用Matlab求其数值解。求解中采用Matlab自带的四阶龙格-库塔（Runge-Kutta）算法，此算法可以在保证高精度的前提下，快速求解微分方程的数值解，其Matlab函数为ode45（）。理论上已证明，在RCSJ模型中接入直流电源，当=时，约瑟夫森结间电压，且电容的存在会使电压的产生滞后。在RCSJ模型中接入交流电源时，会得到台阶状的I-V特性曲线，且台阶是按一定的规律分布的。用上述的方法建立数值模型，借助Matlab进行计算、绘图，可以精确、直观地从数值模

6、拟实验上证明和看到上述现象。此外，RCSJ的数值模型还是一个很好的非线性混沌现象研究工具，因为我们不仅可以在保证精度的前提下，随意的改变电路的各个参数来研究各参量之间的关系和系统进入混沌的过程，而且可以通过绘制相图、庞加莱截面、时序图和功率谱等对混沌产生机理进行分析。4 Matlab模拟结果及分析a. 约瑟夫森结的直流I-V特性由（2）式可知约瑟夫森结平均电压为： （5）因此，计算机模拟约瑟夫森结伏安特性时可用曲线代替I-V特性曲线。计算时用函数mean（）近似地计算电压平均值，可以得到比较精确的结果，绘出的I-V特性曲线如图3、图4所示。图3：约瑟夫森结直流I-V特性曲线 图4：约瑟夫森结直

7、流I-V特性曲线Josephson Junction DC I-V characteristic curve Josephson Junction DC I-V characteristic curve 分析：从图中可以看出，约瑟夫森结注入直流电流后（），当电流 时，U0，此时约瑟夫森结中的电流由正常态电流和超导电流组成。电容的存在并没有改变约瑟夫森结的性质，但会使得电压产生滞后，且随着电容的增大，滞后现象会增强，如图4所示。b. 约瑟夫森结的交流I-V特性图5：约瑟夫森结交流I-V特性曲线 Josephson Junction AC I-V characteristic curve当同时向约

8、瑟夫森结注入直流电流和交流电流，交流的I-V特性曲线上出现台阶现象，如图5所示。数学计算的结果证明了这些台阶出现在： n=0，1，2， （6）对于无量纲化的结果来说，交流台阶出现在： n=0，1，2， （7）c. 约瑟夫森结的高度非线性特性从上述无量纲化过程中，可以知道归一化电压u的表达式为： （8）因此，画分岔图、相图和庞加莱截面时纵轴就是归一化电压u。下面给出了通过Matlab模拟计算得到的非线性混沌现象的分岔图、相图、庞加莱截面及功率谱，也给出了一些简要的分析方法。i. 通向混沌的道路图6：系统u随参数的分岔图 图7：系统u随参数的分岔图The bifurcation diagram o

9、f U with the parameter The bifurcation diagram of U with the parameter 图8：系统u随参数的分岔图 图9：系统u随参数的分岔图 The bifurcation diagram of U with the parameter The bifurcation diagram of U with the parameter 通向混沌的道路有很多种，如倍周期分岔、阵发性方式、茹厄勒塔肯斯道路等。它们是具有普适性的，很多动力学系统随着参数的改变会以倍周期分岔或阵发性的方式进入混沌状态，在混沌区中，由于阵发性的作用，又会出现一些具有周期

10、性的窗口。如图6，图7，图8，图9所示，它们分别是约瑟夫森结电压u随参数、的分岔图。我们先简要分析图6和图7，然后在ii中用相图、时序图、庞加莱截面和功率谱着重分析图9。图6是以为控制参数（0.82.0），其他参数为（=1，）时，得到的u随变化的分岔图。当0.81.15时，系统处于1周期状态；当等于1.15时，系统进入2周期状态，接着是4周期状态、6周期状态、4周期状态；当1.241.43时，系统又处于2周期状态；当1.431.6时，系统处于内部有周期窗口的混沌状态；当1.61.68时，系统是3周期状态；当1.682.0时，系统处于混沌状态，且内部有周期窗口出现。从图6我们看到了不完全级联的过

11、程，即随直流电流的增加，系统由1周期解（单周期的运动）进入到2，4，6周期解（多周期的运动），然后返回到4，2周期解，之后才进入混沌状态，而没有直接进入混沌，这称为不完全级联；在选择的参数范围内，系统随的增加以阵发方式进入混沌状态，且内部有周期窗口。图7，图8（局部放大图）是以为控制参数（01.5），固定其他参数（，）时，得到的u随变化的分岔图。从图中很容易看出：随着参数的减小，周期运动状态和混沌状态会不断交替出现，且混沌区中含有较狭窄的周期窗口，最后完全进入混沌状态；陆续出现周期1，周期3，周期5，周期7，的周期状态区 ，且周期状态区域逐渐变窄，直到消失；在选择的参数范围内，系统随的减小以阵

12、发性方式进入混沌状态；随着参数的减小，系统由周期状态进入到混沌状态的过程是级联过程。 ii. 混沌现象的分析图9是约瑟夫森结电压u随参数的分岔图（）。由图中可以看出，它与前面两种情况不同，该条件下系统是通过倍周期分岔到混沌状态，再由阵法性方式到达周期状态的。由于单用一种方法对混沌进行研究时，在一些情况下并不能达到很好的效果。比如相轨道很密时不易分出是几个周期；画庞加莱截面时两点靠得太近而误以为是单周期，等等。因此，我们用四种方法相互比较对混沌进行研究，分别是相图、庞加莱截面、功率谱和时序图。 周期状态区、混沌状态区和周期窗口在分岔图中显而易见，在此就不列举了。我们截取有代表性的几个值进行研究，得到如下图像。图10：系统相图（Bc=0.5,w=0.66,ia=0.8） Phase Diagram of