二维超导薄膜中拓扑缺陷数目与冷却速度-KUBB相变.docx

二维超导薄膜中拓扑缺陷数目与冷却速度KUBB相变论文导读:：基于含时Ginzburg-Landau(TDGL)方程建立二维超导薄膜经历快速冷却后的拓扑缺陷数目与冷却速度的关系模型。采用差分法对此模型进行了数值模拟，计算出不同冷却速度下拓扑缺陷数目随着时间的变化关系。计算结果发现，在不同冷却速度下，缺陷数目随着时间的变化可以分为两个阶段：第一阶段缺陷数目较多，拓扑缺陷数目随时间急剧减少；在第二阶段，缺陷数目较少，缺陷数目缓慢减少。而且发现，较快的冷却速度有利于形成更多的拓扑缺陷。这种变化趋势与文献[8]报道的实验结果相一致。论文关键词：超导薄膜，TDGL，KUBB相变，拓扑缺陷，冷却速度1引言相变现象一直是凝聚态物理学中的研究热点之一。在文献[1]中，Kibble预言如果一个复序参量系统经历快速冷却后发生相变进入到有序状态，就会产生拓扑缺陷。Zurek[2]发展了Kibble的理论，他预言拓扑缺陷的初始密度是冷却速度的函数。并且KUBB相变，Zurek等人用数值方法研究了一维和二维对称性破缺发生时，拓扑缺陷的数目与冷却速度的关系[3，4]中国期刊全文数据库。以前关于这类现象的研究主要集中在超流体系[5-7]。近期，在实验上二维超导薄膜经历快速冷却后形成拓扑缺陷的现象也有了相关报道[8]。在超导体系中，形成的拓扑缺陷就是磁通量子。但二维超导薄膜经历快速冷缺后形成拓扑缺陷（磁通量子）的数目与冷却速度的关系没有相关的数值模拟。含时Ginzburg-Landau(TDGL)方程是研究超导磁通动力学的一个有力工具[9-15]。基于上述情况，我们基于TDGL方程，采用数值模拟的方法研究了二维超导薄膜在经历快速冷却条件下所形成的拓扑缺陷数目与冷却速度的关系。数值模拟结果发现，在不同冷却速度下KUBB相变，缺陷数目随着时间的变化情况不同，而且较快的冷却速度有利于形成更多的拓扑缺陷。这种变化趋势与文献[8]所报道的实验结果相一致。2基于TDGL方程的二维模型2.1模型的建立通过对TDGL方程的矢量和标量势的规范变换，TDGL方程可以进行数值积分并且用于二维超导薄膜经历快速冷却后的拓扑缺陷数目与冷却速度的关系研究中国期刊全文数据库。本文中对TDGL方程进行的归一化和规范变换如下：，，，，，其中是一个正的常量KUBB相变，是热力学临界磁场强度。和分别定义为和。经过上述变换后，TDGL方程可以改写为（1）其中波函数Ψ以0K时的平衡值为单位，温度T以超导临界温度Tc为单位，T随时间t的变化关系为，为冷却速度。经过上述变换后，建立的方程即为本文中用于二维超导薄膜经历快速冷却后的拓扑缺陷数目与冷却速度的关系模型。2.2数值模拟考虑在X-Y平面上的二维超导薄膜。首先对TDGL方程进行离散化为差分方程。我们取一块单位晶胞（600ξ×600ξ），和分别表示单位晶胞沿x方向和y方向的长度，包含（2400×2400）个晶格点KUBB相变，在此二维区域上离散TDGL方程，采用差分法中的九点差分格式离散空间变量。和分别表示X方向和y方向的网格间隔，取值0.25。在我们的模拟里，采用欧拉方法求解中国期刊全文数据库。其余相关变量的初值分别为=2，=1，=0.5，设置好初值后就开始进行模拟。模拟时的边界条件有自然边界条件和周期性边界条件两种。文中采用后者，因为周期性边界条件能够较好的给出拓扑缺陷的动力学行为。根据本文的TDGL模型KUBB相变，采用的周期性边界条件为：（2）在模拟计算过程中，逐渐增加晶格点数目，最后达到了满意结果。接近Tc时计算较长时间,达到平衡后再降温。2.3结果与讨论本文得到的模拟结果如图1和图2所示。图1显示了不同冷却速度下，缺陷数目随着时间的变化。曲线可以分为两个阶段，第一阶段缺陷数目较多，间距较小，涡旋方向相反的磁通量子由于近程作用而靠近湮灭，拓扑缺陷数目随时间急剧减少；在第二阶段KUBB相变，缺陷数目较少，间距较大，磁通量子间的相互作用较弱，缺陷数目缓慢减少。在这两个阶段的分界点处缺陷数目与冷却速度的关系见图2，较快的冷却速度有助于形成更多的拓扑缺陷中国期刊全文数据库。图1.不同冷却速度下缺陷数目随着时间的变化曲线Fig.1Thenumberofdefectsduringdifferentcoolingrate图2.第一二阶段分界处缺陷数目与冷却速度的关系。插图为文献[8]中的实验结果Fig.2Therelationofnumberofdefectsandcoolingspeedatthedemarcationpointbetweenthefirststageandthesecondstage.Theillustrationistheexperimentsresultsindocument[8].3结论建立了基于TDGL方程的二维超导薄膜在快速冷却条