约瑟夫森结非线性特性的计算机模拟.doc

 物 理 实 验 第 25 卷 第 2 期 2005 年 2 月 Vol . 25 No . 2 Feb. ,2005 P H YSICS EXP ER IM EN TA T ION 约瑟夫森结非线性特性的计算机模拟 戴 岭 ,于 瑶 ,江洪建 ,蒋永兴 (南京大学 物理系 ,江苏 南京 210093) 摘 要 :通过计算机对约瑟夫森结方程求解 ,并模拟约瑟夫森结在直流和交流条件下 V 2I 特性曲线 . 通过约瑟夫森 结的相图与庞加莱截面图直观的表现了倍周期分岔以及混沌现象 ,并与用约瑟夫森结电子模拟器所做的实验进行比较 . 关键词 :约瑟夫森结 ;V2I 曲线 ;计算机模拟 ;分频 中图分类号 :O415 文献标识码 : A 文章编号 :100524642 (2005) 0220025203 dφ 2 e d t = V ( 2) 引 言 1 h- 式中 is 是流经约瑟夫森结的超流 , I0 是流经约瑟 夫森结的最大超流 , V 是约瑟夫森结两端的电压 , φ为 2 块超导体之间的量子相位差 . 图 1 等效电路的方程可以写成 约瑟夫森结是一种典型的弱连接超导体 , 具 有高度的非线性特性 , 如非线性的 V 2I 关系 、频率 变换 、混沌特性等 . 对于其特性的研究可以用数 学方程进行物理实现 , 如可以用相对应和比较的 方法 :大孔径角重摆 , 模拟电路[ 1 ] 等 . 其电路实现 比较直观和易于 控制 的 特点 可以 应 用于 多个 领 + V dV I ( t) = I si n φ + C d t 0 R 2 e I0 t , 对上式做变换 , 令 i (τ) = I ( t) ,τ= ω t = 域[ 2 ] . 但是实验参量的选择和实验的实现本身都 0 hC I0 受到设计 、材料以及可推广性的限制. 随着计算 机速度的不断提高 , 程序设计的简单和表达方式 的多样性和直观性的特点 , 使得计算机模拟实验 逐渐表现出明显的优势 . h G = . 将方程无量纲化后得到 2 e I0 C R2 d2φ dφ i (τ) = dτ2 + G dτ+ si n φ ( 3) 如果同时给约瑟夫森结注入交流电流 , 则 I ( t) = Id + Ia co s ωt , ( 3) 式又可以表示为 实验原理 2 点接触型的约瑟夫森结器件的等效电路如图 1 所示 . i (τ) = Id + Ia co s ωa = i τ + i co s ωτ ( 4) d a ω0 I0 I0 ( 4) 式 第 1 项 为 直 流 分 量 , 第 2 项 为 交 流 分 量. ω= ωa 为约化的外加交流角频率. ω0 微分方程 ( 3) 直接求解比较困难 , 为了得到 V 2I 特性曲线和系统的相图 , 虽然可以通过单摆 模拟法或电路模拟法来实现 , 但是这两种模拟方 法从设计到实现不可避免的要受到系统的参量调 整的限制 , 而采用计算机模拟可以比较方便的解 决上述问题 ( 见图 2) . 图 1 约瑟夫森结的等效电路 理想的约瑟夫森结遵循方程 is = I0 si n φ [ 3 ] 本 文 采 用 Mat la b 软 件 来 实 现 计 算 机 模 拟 , 程序主要采用了 4 阶龙格 2库塔法. 首先将方 ( 1) “第三届全国高等学校物理实验教学研讨会”论文 收稿日期 :2004206201 作者简介 :戴 岭 (1983 - ) ,男 ,江苏沭阳人 ,南京大学物理系 2001 级本科生 . 图 2 用计算机模拟输出 V2I 特性曲线的流程图 程降为 2 个一阶微分方程 , 则 ′y = z ( 5) ′z = - Gz - si n x + id + ia co s ωτ 图 4 结电容对 V2I 特性曲线的影响 初始条件的确定是以上一次的末态作为每次计算 的初始条 件 , 最初 点 的 初 始 条 件 取 为 ( 0 , 0 ) . 由 降为一阶 , 方程将存在唯一的解 . 2) 约瑟夫森结交流 V 2I 特性曲线 当同时向约瑟夫森结注入直流电流 id 和交 流电流 ia , 交流的 V 2I 特性曲线上出现很多的台 阶 , 数学计算的结果证明了这些台阶出现在 - ω0 ( 2) 式可知约瑟夫森结的结电压 V = h 〈 〉, 因 dφ 2 e dτ 此在用计算机模拟 V 2I 特性曲线时可用〈dφ ～ id 〉 dτ 曲线代替 . n h V n = 2 e f ( n = 0 , 1 , 2 , ) 3 计算机模拟约瑟夫 森结的非线性特性 曲线 1) 约瑟夫森结直流 V 2I 特性