1. 引言
2001年日本科学家秋光纯等报道了MgB2具有超导电性,其临界转变温度
高达40 K [1],该发现迅速激起了全世界范围内关于该超导材料的研究热潮。MgB2是由六方紧密排列的镁原子层和类石墨结构的硼原子层交替组成,其中硼原子对二硼化镁超导电性起着关键作用。研究发现MgB2属于常规超导体,符合传统的电声耦合机制 [2],其费米面是由两个二维的σ带和两个三维的π带组成 [3] [4]。随后,有关该超导体两个能隙的图像被比热、核磁共振、角分辨光电子谱和电子隧道等一系列实验广泛证实 [5] [6] [7] [8] [9]。
MgB2超导薄膜对于电子、高场和频射腔等应用具有重要意义 [10] [11] [12],而制备出高质量的超导薄膜是研究其物理特性及实际应用的重要前提。目前人们已对MgB2超导薄膜的沉积技术进行了大量研究,并积累了非常丰富的制备经验,其中混合物理化学气相沉积法(HPCVD)是制备MgB2超导薄膜较为有效的技术 [13] [14] [15]。例如,冯庆荣团队利用HPCVD技术在SiC衬底上制备出一系列厚度为8 μm到8 nm的MgB2超导薄膜 [16] [17],发现当薄膜厚度从30 nm减小到8 nm时,临界温度下降了13%。潘杰云等人在MgO衬底上生长出厚度在40 nm到10 nm的MgB2薄膜,并研究了超导临界温度和厚度间的关系 [18],发现在10 nm的薄膜中,
急剧下降到34 K。于是人们提出了各种不同的机制来解释MgB2超导薄膜中临界温度被抑制的问题,如超导序参量的量子相位波动 [19] 和无序引发的库仑相互作用增强等 [20]。但是到目前为止,对于上述实验数据的解释仍然没有达成共识。
在本文中,我们将针对两带超导体–绝缘体(或真空)之间的边界效应展开研究,并试图探究两带超导体边界效应对MgB2薄膜临界温度的影响。我们首先在超导体界面为超导序参量引入合适的边界条件,然后基于两带Bogoliubov-de Gennes理论得到了MgB2薄膜表面项的特征尺度。然后我们根据两带Ginzburg-Landau (GL)理论计算了MgB2薄膜临界温度与厚度之间的依赖关系,并将理论计算结果与实验数据进行比较,发现理论分析与实验数据一致。我们的工作表明边界效应是理解MgB2薄膜超导电性的一个重要因素。
本文主要由四部分组成,在接下来的一部分中,我们回顾了两带Bogoliubov-de Gennes理论和GL方程,对自由能中引入的新表面项进行讨论。第三部分是MgB2临界温度的理论计算。第四部分是文章的结论。
2. 两带Ginzburg-Landau理论和超导体–绝缘体界面的边界条件
根据Zhitomirsky和Dao的理论工作 [21],我们可以写出两带超导体的哈密顿量为:
(1)
方程中
指能隙序数,
是自旋态。
表示正常金属单粒子哈密顿量,
为散射振幅且
。
为简化哈密顿量的形式,引入能隙函数
(2)
在平均场近似下,我们由此可以得到有效哈密顿量的形式为
(3)
通过引入描述准粒子激发的产生湮灭算符
,我们可以利用Bogoliubov变换对有效的哈密顿量形式进行对角化
(4)
(5)
方程中
代表波矢量。对角化后的哈密顿量可写为
(6)
是基态的能量,
是激发态的能量。
利用对易关系
和公式(4)~(6),我们可以推导出两带Bogoliubov-de Gennes方程
(7)
由公式(6),我们同时也能得到
,其中
。结合公式(2),可得能隙函数的自洽方程为
(8)
与单带类似 [22],当
为小量时,我们可以通过公式(7)和(8)得到能隙自洽方程的线性表达式
(9)
其中
(10)
是正常金属电子本征函数,
。
代表频率,
是整数。
我们结合块体材料
的精确表达形式并在能隙方程中加入非线性项,由公式(9)可以得到两带GL方程 [18]
(11)
(12)
方程中参数取值如下
,
(13)
,
(14)
其中
,
为两个能带各自费米能级态密度。
以及
为
矩阵的最大特征值。
是块体材料的临界温度,
表示每个能带的平均费米速度。
在各向同性的均匀体系中,我们可以略去两带GL方程的梯度项。当
时,由(11)、(12)式得能隙自洽方程
(15)
与此同时,根据公式(11)和(12)我们可以写出弱耦合的两带GL自由能表达式
(16)
其中
(17)
以及
(18)
考虑超导体–绝缘体(或真空)间存在界面,我们添加新表面项
(19)
该积分在超导体表面S进行,其中
为常数。
于是总的自由能可表示为
(20)
对
做变分,得
(21)
当满足条件
时,我们不仅可以得到两带GL方程(11)和(12),而且可得边界条件
(22)
对于通常应用于超导体–绝缘体(或真空)界面的Neumann边界条件,我们有
,但在此条件下超导临界温度与薄膜厚度的依赖关系无法得到合理解释。我们基于de Gennes [22] 关于单带超导体边界条件的理论计算,推广得到两带超导体边界条件中
的表达式为
(23)
我们认为由此表面项引起的边界效应将为解释
被抑制的行为提供一种可能。
3. MgB2超导薄膜临界温度的计算
在本节中,我们将根据边界效应解释临界温度
与MgB2薄膜厚度之间的依赖关系。设薄膜厚度为d,其边界位置由
和
决定。
根据公式(11)和(12),两带GL方程可写为
(24)
其中
(25)
和
(26)
由于临界温度附近序参量较小,我们忽略了包含
的高阶项贡献。
与单带超导体类似 [23],我们设超导薄膜能隙方程的解为
(27)
为常数。由边界条件
(28)
可得
(29)
由
(30)
得波矢k满足方程
(31)
同时,从公式(24)和(27)出发,两带GL方程可改写为
(32)
其中
(33)
和
(34)
在
时,决定两带超导薄膜临界温度的条件为
(35)
对于典型的两带超导体MgB2,已知
[1],平均晶格常数
[24]。
带费米面的态密度是
,
带费米面的态密度是
[2]。由比热数据的理论拟合得
,比值
[21],据此我们可以计算出
,
,
。结合费米平均速度
,
[2],公式(14)给出
,
。将所得参数带入公式(23),最终得特征尺度
,
,
,
。基于公式(31)和(35),我们给出了超导临界温度随薄膜厚度的变化趋势,如图1所示。从图中可以看到临界温度随薄膜厚度的减小而逐渐降低,当膜厚度降低到10 nm时,临界温度
。值得注意的是,该理论结果与不同衬底生长的MgB2实验数据均相符。
Figure 1.The critical temperature as a function of MgB2 film thickness
图1. 超导临界温度随MgB2薄膜厚度的变化曲线
4. 结论
在本文中,我们对两带超导体–绝缘体(或真空)之间的边界效应展开细致研究,发现对于典型的两带超导体二硼化镁Neumann边界条件无法解释薄膜厚度与临界温度的依赖关系。于是我们在两带GL理论中引入了新表面项,并给出表面项参数的微观表达式,该理论能够很好地解释当薄膜厚度减少时临界温度被显著抑制的实验现象。研究结果表明,在MgB2等多带超导薄膜的研究中由新表面项引起的边界效应可能产生不可忽略的物理效应。