本文是「固体物理学」系列的第 34 篇(共 34 篇),整理自 Neil W. Ashcroft 与 N. David Mermin 所著 Solid State Physics 第 34 章的中文译稿。
第三十二章指出,独立电子近似不能充分描述大多数磁有序固体。在许多没有磁有序的金属中,独立电子近似在极低温下突然发生更为显著的失效,此时形成另一种电子有序态,称为超导态(superconducting state)。超导性并非少数金属所特有。超过 20 种金属元素可以成为超导体(表 34.1)。在适当条件下,甚至某些半导体也可以成为超导体1;已测量超导性质的合金则多达数千种。2
表 34.1 超导元素a
| H | He | ||||||||||||||||
| Li | Be● | B | C | N | O | F | Ne | ||||||||||
| Na | Mg | Al | Si● | P | S | Cl | Ar | ||||||||||
| K | Ca | Sc | Ti | V | Cr | Mn | Fe | Co | Ni | Cu | Zn | Ga | Ge● | As | Se● | Br | Kr |
| Rb | Sr | Y | Zr | Nb | Mo | Tc | Ru | Rh | Pd | Ag | Cd | In | Sn | Sb● | Te● | I | Xe |
| Cs● | Ba● | Lu | Hf | Ta | W | Re | Os | Ir | Pt | Au | Hg | Tl | Pb | Bi● | Po | At | Rn |
| Fr | Ra | ||||||||||||||||
| La | Ce● | Pr | Nd | Pm | Sm | Eu | Gd | Tb | Dy | Ho | Er | Tm | Yb | ||||
| Ac | Th | Pa | U | Np | Pu |
a 仅在特殊条件下超导的元素另作标记。注意超导有序与磁有序的不相容性。据 G. Gladstone 等,Parks,op. cit.,脚注 6。
图例:
| Al 超导 | B 非金属元素 |
| Si● 在高压或薄膜条件下超导 | Fe 具有磁有序的元素 |
| Li 金属性,但尚未发现超导 |
从独立电子近似的角度看,金属在超导态中的特征性质显得十分反常。超导体最显著的特征是:
- 超导体可以表现得如同没有可测量的直流(DC)电阻率。在超导体中已建立电流,即使没有任何驱动场,只要人们耐心持续观测,始终未发现可辨别的衰减。3
- 超导体可以表现为完全抗磁体(perfect diamagnet)。样品在外加磁场中处于热平衡时,只要磁场不过强,就存在表面电流。这些电流产生附加磁场,在超导体内部恰好抵消外加磁场。
- 超导体通常表现得如同在允许的单电子能级集合中,存在一个以费米能为中心、宽度为 $2\Delta$ 的能隙。4因此,只有当 $\mathcal E-\mathcal E_F$(或 $\mathcal E_F-\mathcal E$)超过 $\Delta$ 时,能量为 $\mathcal E$ 的电子才能被超导体接纳(或从中取出)。5能隙 $\Delta$ 随温度下降而增大,在极低温下趋于最大值 $\Delta(0)$。
超导理论内容十分广泛,且高度专门化。与本书其他地方描述的理论一样,它以电子和离子的非相对论量子力学为基础,但除此之外,它与此前考察的其他模型和理论的相似性迅速减少。超导的微观理论不能用独立电子近似的语言描述。即使相对初等的超导体微观计算,也依赖于形式化技术(场论方法,field theoretic methods)。这些方法在概念上并不比普通量子力学方法更复杂,但必须积累相当多的经验和实践,才能有把握地理解和使用。
因此,与其他章节相比,这里将更严格地把超导理论的概述限制为对某些主要概念的定性描述,以及对少数较简单预言的陈述。即使只希望获得本领域初步实用知识的读者,也必须参考现有的众多著作之一。6
本章安排如下:
- 概述超导的基本实验事实。
- 描述唯象伦敦方程及其与完全抗磁性的关系。
- 定性描述巴丁、库珀和施里弗的微观理论。
- 概括微观理论的一些基本平衡性质预言,并与实验比较。
- 定性讨论微观理论、“序参量”(order parameter)概念和超导体输运性质之间的关系。
- 描述 B. D. Josephson 预言的超导体之间的奇特隧穿现象。
临界温度
块体样品向超导态的转变是突变的。在临界温度7 $T_c$ 以上,金属的性质完全正常;在 $T_c$ 以下则表现出超导性质,其中最显著的是不存在任何可测量的直流(DC)电阻。测得的临界温度从几毫开尔文8到略高于 20 K 不等。相应的热能 $k_BT_c$ 从约 $10^{-7}$ eV 到几千分之一电子伏特不等。与固体中通常认为重要的能量相比,这极其微小。9超导元素的转变温度列于表 34.2。
持续电流
图 34.1 显示超导金属跨过临界温度 $T_c$ 时电阻率随温度的变化。在 $T_c$ 以上,电阻率具有正常金属的特征形式 $\rho(T)=\rho_0+BT^5$:常数项来自杂质10和缺陷散射,$T^5$ 项来自声子散射。在 $T_c$ 以下,这些机制不再使电流衰减,电阻率骤降为零。电流可以在超导体中流动,而没有可辨别的能量耗散。11不过,存在一些限制:
- 施加足够大的磁场会破坏超导性(见下文)。
表 34.2 超导元素的 $T_c$ 和 $H_c$ 值a
| 元素 | $T_c$ (K) | $H_c$(GAUSS)b | |
|---|---|---|---|
| Al | 1.196 | 99 | |
| Cd | 0.56 | 30 | |
| Ga | 1.091 | 51 | |
| Hf | 0.09 | — | |
| Hg | $\alpha$ (rhomb) | 4.15 | 411 |
| $\beta$ | 3.95 | 339 | |
| In | 3.40 | 293 | |
| Ir | 0.14 | 19 | |
| La | $\alpha$ (hcp) | 4.9 | 798 |
| $\beta$ (fcc) | 6.06 | 1096 | |
| Mo | 0.92 | 98 | |
| Nb | 9.26 | 1980 | |
| Os | 0.655 | 65 | |
| Pa | 1.4 | — | |
| Pb | 7.19 | 803 | |
| Re | 1.698 | 198 | |
| Ru | 0.49 | 66 | |
| Sn | 3.72 | 305 | |
| Ta | 4.48 | 830 | |
| Tc | 7.77 | 1410 | |
| Th | 1.368 | 162 | |
| Ti | 0.39 | 100 | |
| Tl | 2.39 | 171 | |
| U | $\alpha$ | 0.68 | — |
| $\gamma$ | 1.80 | — | |
| V | 5.30 | 1020 | |
| W | 0.012 | 1 | |
| Zn | 0.875 | 53 | |
| Zr | 0.65 | 47 | |
a 对于第二类超导体(type II superconductors),所列零温临界磁场通过等面积构造得到:将低场($H\lt H_{c1}$)磁化强度线性外推至某个磁场 $H_c$,选择该磁场,使围成的面积等于实际磁化曲线下的面积。
b 在 $T=0$(K)时。
来源:B. W. Roberts,Progr. Cryog. 4,161(1964);G. Gladstone、M. A. Jensen 和 J. R. Schrieffer,Superconductivity,R. D. Parks 编,Dekker,New York,1969;Handbook of Chemistry and Physics,第 55 版,Chemical Rubber Publishing Co.,Cleveland,1974–1975。
图 34.1 (a)含非磁性杂质的正常金属的低温电阻率($\rho(T)=\rho_0+BT^5$)。(b)含非磁性杂质的超导体在零磁场下的低温电阻率。在 $T_c$ 处,$\rho$ 骤降为零。
- 如果电流超过“临界电流”(critical current),超导态会被破坏(西尔斯比效应,Silsbee effect)。临界电流的大小(在 1-mm 导线中可高达 100 amp)取决于样品的性质和几何形状,并与电流产生的磁场是否超过超导体表面的临界磁场有关。12
- 远低于转变温度的超导体,在频率不太高时,对交流(AC)电场也会作出无耗散响应。从无耗散响应转变为正常响应,发生在 $\Delta/\hbar$ 量级的频率 $\omega$ 处,其中 $\Delta$ 为能隙。
热电性质
在独立电子近似中,良好的电导体也是良好的热导体,因为传导电子不仅输运电荷,也输运熵。13与此相反,超导体是热的不良导体(图 34.2)。14它们也不表现出珀尔帖效应(Peltier effect);即在温度均匀的超导体中,电流不像正常金属中的电流那样伴随热流。没有珀尔帖效应表明,参与持续电流的电子不携带熵。较差的热导率表明,即使超导体中没有电流,也只有一部分传导电子能够输运熵。15
图 34.2 铅的热导率。在 $T_c$ 以下,下方曲线给出超导态热导率,上方曲线给出正常态热导率。通过施加磁场,在 $T_c$ 以下获得正常态样品;假定磁场除此之外对热导率没有显著影响。(经加拿大国家研究委员会许可,转载自 J. H. P. Watson 和 G. M. Graham,Can. J. Phys. 41,1738(1963)。)
磁性质:完全抗磁性
磁场只要不过强,就不能穿透超导体内部。迈斯纳–奥克森费尔德效应(Meissner-Ochsenfeld effect)最显著地说明这一点:如果将磁场中的正常金属16冷却至超导转变温度以下,磁通会被突然排出。因此,在磁场中发生的转变伴随着表面电流的出现,这些电流恰好抵消样品内部磁场。
注意,单凭完全导电性(perfect conductivity,即 $\sigma=\infty$)并不能推出这一结果,尽管完全导电性确实意味着某种相关性质:如果将最初处于零磁场的完全导体移入非零磁场区域(或开启磁场),法拉第感应定律会产生涡流,抵消内部磁场。然而,如果完全导体内部原已建立磁场,排出该磁场也同样受到阻碍。若将样品移入无磁场区域(或关闭外加磁场),会感应出涡流以维持磁场。因此,完全导电性意味着内部磁场与时间无关,却不规定该磁场必须具有何值。在超导体中,磁场不仅与时间无关,而且为零。下面讨论伦敦方程时,将更定量地考察完全导电性与迈斯纳效应之间的关系。
磁性质:临界磁场
考虑温度 $T$ 低于临界温度 $T_c$ 的超导体。开启磁场 $H$ 时,需消耗一定能量来建立屏蔽电流(screening currents)的磁场,以抵消超导体内部磁场。如果外加磁场足够大,样品恢复正常态并允许磁场穿透,就会在能量上变得有利。因为虽然在零磁场下、$T_c$ 以下,正常态的自由能高于超导态,但磁场足够强时,屏蔽电流消失、磁场得以进入样品所引起的磁场能降低,会超过正常态自由能增加。
磁场强度增大时,穿透方式一般取决于样品的几何形状。不过,对最简单的几何形状,即轴线平行于外加磁场的细长圆柱形样品,存在两种可以明确区分的行为:
第一类(Type I) 在临界磁场 $H_c(T)$ 以下,磁通不穿透;随着 $T$ 降至 $T_c$ 以下,$H_c(T)$ 增大。当外加磁场超过 $H_c(T)$ 时,整个样品恢复正常态,磁场完全穿透。17相应的 $H$-$T$ 平面相图见图 34.3。18通常将宏观抗磁磁化密度 $M$ 对外加磁场 $H$ 作图,描述这种磁场穿透(图 34.4a)。
图 34.3 第一类超导体在 $H$-$T$ 平面中超导态与正常态之间的相界。相界由曲线 $H_c(T)$ 给出。
第二类(Type II) 在下临界磁场(lower critical field)$H_{c1}(T)$ 以下,磁通不穿透;当外加磁场超过上临界磁场(upper critical field)$H_{c2}(T)\gt H_{c1}(T)$ 时,整个样品恢复正常态,磁场完全穿透。当外加磁场强度介于 $H_{c1}(T)$ 和 $H_{c2}(T)$ 之间时,磁通部分穿透,且样品形成正常区域与超导区域共存的相当复杂的微观结构,称为混合态(mixed state)。19第二类行为对应的磁化曲线见图 34.4b。
图 34.4 (a)第一类超导体的磁化曲线。在 $H_c$ 以下,磁场不穿透:$B=0$(或 $M=-H/4\pi$)。(超导体中 $B$ 和 $H$ 的区别见脚注 30。)(b)第二类超导体的磁化曲线。在 $H_{c1}$ 以下,行为与第一类情况相同。在 $H_{c1}$ 和 $H_{c2}$ 之间,$M$ 平滑下降至零,$B$ 平滑上升至 $H$。
A. A. Abrikosov 提出,并随后由实验(图 34.5)证实:在混合态中,磁场以细磁通丝的形式部分穿透样品。每根磁通丝内部磁场很强,材料不处于超导态。在磁通丝的芯部以外,材料仍然保持超导,磁场按照伦敦方程(见下文)确定的方式衰减。每根磁通丝周围都有屏蔽电流形成的涡旋。20
图 34.5 在垂直于表面的 80 gauss 磁场下,穿出 $\mathrm{Pb}_{.98}\mathrm{In}_{.02}$ 超导箔表面的涡旋线(vortex lines)三角排列。(由 J. Silcox 和 G. Dolan 提供。)细小铁磁颗粒的凝聚显现出涡旋。相邻涡旋间距约为半微米。
在远低于转变温度时,第一类超导体的典型临界磁场约为 $10^2$ gauss。不过,在所谓“硬”第二类超导体中,上临界磁场可高达 $10^5$ gauss,因此第二类材料在强磁场磁体设计中具有相当重要的实际意义。
元素超导体的低温临界磁场见表 34.2。
比热
在低温下,正常金属的比热具有形式 $AT+BT^3$,其中线性项来自电子激发,三次项来自晶格振动。在超导临界温度以下,这一行为发生显著变化。当温度降至 $T_c$ 以下时(零磁场中),比热跃升至更高值,随后缓慢下降,最终远低于正常金属的预期值(图 34.6)。通过施加磁场使金属进入正常态,可以比较临界温度以下超导态和正常态的比热。21这种分析表明,在超导态中,比热的线性电子贡献被一项在极低温下更快趋于零的贡献取代,其主导低温行为具有形式 $\exp(-\Delta/k_BT)$。这是激发能级与基态之间相隔能量 $2\Delta$ 的系统所特有的热行为。22理论(见式(34.19))和实验(见表 34.3)均表明能隙 $\Delta$ 为 $k_BT_c$ 的量级。
图 34.6 正常态和超导态铝的低温比热。通过施加较弱的(300-gauss)磁场,在 $T_c$ 以下获得正常相;该磁场破坏超导有序,但除此之外对比热的影响可忽略。铝的德拜温度相当高,因此在整个温度区间内,比热由电子贡献主导(正常态曲线十分接近线性这一事实可以说明这一点)。$T_c$ 处的不连续性与理论预言(34.22)$[c_s-c_n]/c_n=1.43$ 很好地符合。在远低于 $T_c$ 时,$c_s$ 降至远低于 $c_n$,暗示存在能隙。(N. E. Phillips,Phys. Rev. 114,676(1959)。)
能隙的其他表现
正常隧穿(Normal Tunneling)
使超导体与正常金属紧密接触,以致两者仅由薄绝缘层隔开,可以使两者的传导电子相互达到热平衡。23电子可通过量子力学隧穿越过该层。在热平衡下,已有足够多的电子从一种金属进入另一种金属,使两者中电子的化学势相等。24当两种金属都处于正常态时,施加电势差会使一种金属的化学势相对于另一种升高,更多电子便隧穿绝缘层。实验发现,正常金属结中的这种“隧穿电流”遵循欧姆定律。然而,当其中一种金属是温度远低于临界温度的超导体时,只有电势 $V$ 达到阈值 $eV=\Delta$ 才观测到电流(见图 34.7)。$\Delta$ 的大小与低温比热测量推得的值很好地符合,证实超导体单电子能级密度中存在能隙的图像。随着温度升至接近 $T_c$,阈值电压下降25,表明能隙本身也随温度升高而减小。
频率相关的电磁行为(Frequency Dependent Electromagnetic Behavior)
金属对电磁辐射的响应(如薄膜透射或块体样品反射)由频率相关的电导率决定。这又取决于给定频率下传导电子可以利用的能量吸收机制。由于超导态的电子激发谱以能隙 $\Delta$ 为特征,预期频率远小于 $\Delta/\hbar$ 时,交流(AC)电导率显著不同于正常态形式;频率远大于 $\Delta/\hbar$ 时,超导态与正常态的电导率则基本相同。除了非常接近临界温度的区域(见第 744 页)外,$\Delta/\hbar$ 通常处于微波与红外频率之间。在光学频率下,观察到超导态的交流(AC)行为与正常态无法区分。偏离正常态行为的现象首先在红外出现;只有到微波频率时,充分显示能隙所特有的缺乏电子吸收的交流行为才完全显现。
图 34.7 (a)电子通过两个正常金属之间的薄绝缘势垒隧穿时的电流–电压关系。电流和电压较小时,关系为线性。(b)电子通过超导体与正常金属之间的薄绝缘势垒隧穿时的电流–电压关系。该关系强烈依赖于温度。在 $T=0$ 时,存在清晰阈值;在较高温度下,超导体内电子跨越能隙的热激发使阈值变得模糊。
声衰减(Acoustic Attenuation)
声波在金属中传播时,离子位移产生的微观电场可以向费米能级附近电子传递能量,从而带走声波能量。26在远低于 $T_c$ 时,超导体中的衰减率显著低于正常金属,正如对于满足 $\hbar\omega\lt 2\Delta$ 的声波所预期的那样。
伦敦方程
F. London 和 H. London 首次定量考察了如下基本事实:超导态金属不允许其内部存在磁场。27他们的分析从 Gorter 和 Casimir 的双流体模型(two-fluid model)出发。28这里采用该模型的唯一关键假设是:在温度 $T\lt T_c$ 的超导体中,只有占传导电子总数 $n_s(T)/n$ 的电子能够参与超电流(supercurrent)。量 $n_s(T)$ 称为超导电子密度。当 $T$ 降至远低于 $T_c$ 时,它趋于全部电子密度 $n$;当 $T$ 升至 $T_c$ 时,它降为零。假定其余电子构成密度为 $n-n_s$ 的“正常流体”,不能在没有通常耗散的情况下携带电流。假定正常电流与超电流并行流动;由于后者完全没有电阻,它将承担任何微小瞬态电场所感应的全部电流,正常电子则基本不响应。因此,下述讨论忽略正常电子。
设超导体内部瞬时出现电场。超导电子将无耗散地自由加速,使其平均速度 $\boldsymbol{v}_s$ 满足29
$$ m\frac{d\boldsymbol{v}_s}{dt}=-e\boldsymbol{E}. \tag{34.1} $$
由于这些电子携带的电流密度为 $\boldsymbol{j}=-e\boldsymbol{v}_sn_s$,式(34.1)可写为
$$ \frac d{dt}\boldsymbol{j}=\frac{n_se^2}m\boldsymbol{E}. \tag{34.2} $$
注意,对式(34.2)作傅里叶变换,得到德鲁德模型中密度为 $n_s$ 的电子气在弛豫时间 $\tau$ 趋于无穷大时的通常交流电导率,即式(1.29):
$$ \begin{gathered}\boldsymbol{j}(\omega)=\sigma(\omega)\boldsymbol{E}(\omega),\\\sigma(\omega)=i\frac{n_se^2}{m\omega}.\end{gathered} \tag{34.3} $$
将式(34.2)代入法拉第感应定律,
$$ \nabla\times\boldsymbol{E}=-\frac1c\frac{\partial\boldsymbol{B}}{\partial t}, \tag{34.4} $$
得到电流密度与磁场之间的下列关系:
$$ \frac\partial{\partial t}\left(\nabla\times\boldsymbol{j}+\frac{n_se^2}{mc}\boldsymbol{B}\right)=0. \tag{34.5} $$
这一关系与麦克斯韦方程30
$$ \nabla\times\boldsymbol{B}=\frac{4\pi}c\boldsymbol{j}, \tag{34.6} $$
共同决定完全导体内部可以存在的磁场和电流密度。
特别注意,任意静态磁场 $\boldsymbol{B}$ 都通过式(34.6)确定静态电流密度 $\boldsymbol{j}$。由于任意与时间无关的 $\boldsymbol{B}$ 和 $\boldsymbol{j}$ 都显然是式(34.5)的解,这两个方程与任意静态磁场相容。这与超导体的观测行为不符,因为超导体内部不允许存在磁场。F. London 和 H. London 发现,将式(34.5)的全部解限制为满足下式的解,就可以得到超导体的这种特征行为:31
$$ \nabla\times\boldsymbol{j}=-\frac{n_se^2}{mc}\boldsymbol{B}, \tag{34.7} $$
这称为伦敦方程。式(34.5)描述任何无耗散导电介质,要求 $\nabla\times\boldsymbol{j}+(n_se^2/mc)\boldsymbol{B}$ 与时间无关;限制更强的伦敦方程专门描述超导体,并将其与单纯的“完全导体”区分,还要求这一与时间无关的值为零。用限制更强的伦敦方程替代式(34.5),是因为伦敦方程直接导出迈斯纳效应。32式(34.6)和(34.7)意味着
$$ \begin{gathered}\nabla^2\boldsymbol{B}=\frac{4\pi n_se^2}{mc^2}\boldsymbol{B},\\\nabla^2\boldsymbol{j}=\frac{4\pi n_se^2}{mc^2}\boldsymbol{j}.\end{gathered} \tag{34.8} $$
这些方程进而预言:超导体中的电流和磁场只能存在于表面厚度为 $\Lambda$ 的层内;$\Lambda$ 称为伦敦穿透深度,其值为33
$$ \Lambda=\left(\frac{mc^2}{4\pi n_se^2}\right)^{1/2}=41.9\left(\frac{r_s}{a_0}\right)^{3/2}\left(\frac n{n_s}\right)^{1/2}\ \text{Å}. \tag{34.9} $$
因此,伦敦方程不仅意味着迈斯纳效应,也给出了屏蔽外加磁场的表面电流的具体图像。这些电流存在于厚度为 $10^2$–$10^3$ Å 的表层内(在远低于 $T_c$ 时;在临界温度附近,$n_s$ 趋于零,厚度可能大得多)。在同一表层内,磁场连续降至零。厚度等于或小于穿透深度 $\Lambda$ 的超导薄膜中,磁场并非完全穿透,这一事实证实了这些预言。
微观理论:定性特征
超导的微观理论由巴丁(Bardeen)、库珀(Cooper)和施里弗(Schrieffer)于 1957 年提出。34在这种广泛概述中,无法建立充分描述其理论所需的形式体系,只能定性描述其基本物理原理和主要理论预言。
超导理论首先要求费米面附近的电子之间存在净吸引相互作用。尽管直接静电相互作用是排斥的,离子运动仍可能对库仑相互作用产生“过度屏蔽”(overscreen),从而导致净吸引。35此前已描述过这一可能性在第二十六章中。在一个简化模型中发现,允许离子响应电子运动而运动,导致波矢为 $\boldsymbol{k}$ 和 $\boldsymbol{k}'$ 的电子之间具有形式如下的净相互作用:36
$$ v^{\mathrm{eff}}_{\boldsymbol{k},\boldsymbol{k}\prime}(\boldsymbol{k},\boldsymbol{k}')=\frac{4\pi e^2}{q^2+k_0^2}\cdot\frac{\omega^2}{\omega^2-\omega_{\boldsymbol{q}}^2}, \tag{34.10} $$
校对注:式(34.10)第二个下标按本段定义补撇号,与两个电子的波矢一致。
其中 $\hbar\omega$ 是电子能量之差,$k_0$ 是托马斯–费米波矢(17.50),$\boldsymbol{q}$ 是电子波矢之差,$\omega_{\boldsymbol{q}}$ 是波矢为 $\boldsymbol{q}$ 的声子的频率。
因此,离子运动产生的屏蔽可以使能量足够接近的电子之间具有净吸引相互作用(大致是能量差小于 $\hbar\omega_D$,后者衡量典型声子能量)。这种吸引37是超导理论的基础。
既然能量相差 $O(\hbar\omega_D)$ 的电子可以受到净吸引,就可能形成束缚电子对。38这一可能性似乎值得怀疑,因为在三维中,两个粒子间的相互作用必须达到某个最小强度才能形成束缚态,而这种相当有限的有效吸引不大可能满足这一条件。然而,库珀39指出,其余 $N-2$ 个电子通过泡利不相容原理对相互作用电子对施加影响,使这一看似不合理的可能性变得相当可能。
库珀考虑了两个电子的问题:它们之间的吸引相互作用弱到在孤立情况下远不足以束缚它们。然而,他证明,在由额外电子构成的费米球存在时40,不相容原理根本改变了双电子问题,使得无论吸引多弱都存在束缚态。库珀的计算不仅表明束缚电子对的净吸引不必具有最小强度,还说明超导转变温度为何能远低于固体的其他所有特征温度。这来自其解的形式:吸引较弱时,束缚能远小于吸引势能。
库珀的论证适用于额外电子具有正常费米分布时的单个电子对。巴丁、库珀和施里弗的理论迈出了关键的进一步步骤,构造了所有电子都形成束缚对的基态。这显著推广了库珀模型,因为每个电子现在具有双重作用:通过不相容原理,对允许波矢施加必要限制,使其他电子对即使吸引很弱也能束缚;同时,该电子本身也参与一个束缚对。
电子基态波函数的 BCS 近似可描述如下:将 $N$ 个传导电子分成 $N/2$ 对41,每对用束缚态波函数 $\phi(\boldsymbol{r} s,\boldsymbol{r}'s')$ 描述,其中 $\boldsymbol{r}$ 为电子位置,$s$ 为自旋量子数。然后考虑由 $N/2$ 个相同的双电子波函数直接相乘得到的 $N$ 电子波函数:
$$ \Psi(\boldsymbol{r}_1s_1,\ldots,\boldsymbol{r}_Ns_N)=\phi(\boldsymbol{r}_1s_1,\boldsymbol{r}_2s_2)\cdots\phi(\boldsymbol{r}_{N-1}s_{N-1},\boldsymbol{r}_Ns_N). \tag{34.11} $$
这描述了所有电子成对束缚于相同双电子态的状态。然而,它缺少泡利原理所要求的对称性。为构造交换任意两个电子的空间和自旋坐标时变号的状态,必须对状态(34.11)作反对称化。这给出 BCS 基态:42
$\displaystyle \Psi_{\mathrm{BCS}}=$
$\displaystyle \Psi.$(34.12)
[校对注:式中及脚注 42 的反对称化算符保留扫描字形,其字母与字体无法可靠确认。[?]]
尽管状态(34.12)中出现的所有电子对波函数 $\phi$ 都相同,它仍满足泡利原理,这或许令人惊讶。实际上,如果用 $N$ 个相同单电子能级构造类似式(34.11)的乘积态,随后反对称化就会使它消失。反对称性的基本要求意味着:当状态是单电子能级的反对称化乘积时,没有单电子能级可以被双重占据。然而,反对称性要求并不对双电子能级的反对称化乘积态中的双电子能级占据施加相应限制。43
可以证明,若将状态(34.12)作为变分估计基态能量的试探态,那么对于任何吸引相互作用,无论多弱,$\phi$ 的最优选择都必须比斯莱特行列式(Slater determinants)的最佳选择(即最佳独立电子试探函数)给出更低能量。
在 BCS 理论中,取电子对波函数 $\phi$ 为单重态;44即电子对中的两个电子自旋相反,波函数的轨道部分 $\phi(\boldsymbol{r},\boldsymbol{r}')$ 为对称函数。如果取电子对态具有平移不变性(忽略晶格周期势可能造成的复杂情况),使 $\phi(\boldsymbol{r},\boldsymbol{r}')$ 具有形式 $\chi(\boldsymbol{r}-\boldsymbol{r}')$,则可写为
$$ \chi(\boldsymbol{r}-\boldsymbol{r}\prime)=\frac1V\sum_{\boldsymbol{k}}\chi_{\boldsymbol{k}}e^{i\boldsymbol{k}\cdot\boldsymbol{r}}e^{-i\boldsymbol{k}\cdot\boldsymbol{r}\prime}. \tag{34.13} $$
因此,$\chi$ 可视为单电子能级乘积的叠加,其中每一项都将波矢大小相等、方向相反的电子配成对。45
对 $\Psi_{\mathrm{BCS}}$ 作变分计算的一个结果是:电子对波函数46的空间范围 $\xi_0$ 远大于电子间距 $r_s$。可以如下粗略估计 $\xi_0$:电子对波函数 $\phi(\boldsymbol{r})$ 想必是能量与 $\mathcal E_F$ 相差 $O(\Delta)$ 以内的单电子能级的叠加,因为隧穿实验表明,在这一能量范围以外,单电子能级密度与正常金属中的形式相比改变很小。因此,构成电子对态的单电子能级的动量展宽由下列条件确定:
$$ \Delta=\delta\mathcal E=\delta\left(\frac{p^2}{2m}\right)=\left(\frac{p_F}m\right)\delta p\approx v_F\,\delta p. \tag{34.14} $$
于是,$\phi(\boldsymbol{r})$ 的空间范围为下列量级:
$$ \xi_0\sim\frac\hbar{\delta p}\sim\frac{\hbar v_F}\Delta\sim\frac1{k_F}\frac{\mathcal E_F}\Delta. \tag{34.15} $$
由于 $\mathcal E_F$ 通常为 $\Delta$ 的 $10^3$–$10^4$ 倍,$k_F$ 为 $10^8$ $\mathrm{cm}^{-1}$ 的量级,$\xi_0$ 通常为 $10^3$ Å。
因此,在任一给定电子对占据的区域内,会存在许多(数百万或更多)电子对的中心。这是超导态极为关键的特征:电子对不能被视为独立粒子,而是以非常复杂的方式在空间中相互交织,这对该态的稳定性至关重要。
上述描述概括了超导体电子基态的基本特征。为描述超导体的激发态、热性质或输运性质,必须采用更精细的形式体系。这里不展开这些内容,只强调基本物理图像仍是配对电子系统。在非平衡过程中,电子对态可以更复杂。在非零温度下,一部分电子对因热运动而解离,超导电子密度 $n_s$ 由仍保持配对的比例决定。此外,由于配对具有复杂的自洽性质,在非零温度下,部分电子对的热解离会使仍束缚的电子对的特征性质(如电子对函数的范围)依赖于温度。当 $T$ 升高并越过 $T_c$ 时,所有电子对都解离,基态连续恢复为独立电子近似的正常基态。
初等微观理论的定量预言
在最简单形式下,BCS 理论对描述传导电子的基本哈密顿量作了两项严重的过度简化:
- 用自由电子近似处理传导电子;忽略能带结构效应。
- 将费米能附近电子之间相当复杂的净吸引相互作用47(34.10)进一步简化为有效相互作用 $V$。电子波矢为 $\boldsymbol{k}_1$ 和 $\boldsymbol{k}_2$ 的双电子态与波矢为 $\boldsymbol{k}_3$ 和 $\boldsymbol{k}_4$ 的另一态之间的 $V$ 矩阵元,在体积 $\Omega$ 中取为
$$ \begin{aligned}\langle\boldsymbol{k}_1\boldsymbol{k}_2|V|\boldsymbol{k}_3\boldsymbol{k}_4\rangle&=-V_0/\Omega,\ \text{当 }\boldsymbol{k}_1+\boldsymbol{k}_2=\boldsymbol{k}_3+\boldsymbol{k}_4,\\ &\qquad|\mathcal E(\boldsymbol{k}_i)-\mathcal E_F|\lt \hbar\omega,\ i=1,\ldots,4,\\ &=0,\ \text{其他情况。}\end{aligned} \tag{34.16} $$
任何具有平移不变性的势都要求波矢满足上述限制;相互作用(34.16)的重要方面是:只要四个自由电子能量都处于与费米能相差 $\hbar\omega$(通常取为 $\hbar\omega_D$ 的量级)的范围内,就存在吸引。
式(34.16)严重地过度简化了实际净相互作用,任何依赖其详细特征的结果都应谨慎看待。幸运的是,理论预言了若干不含两个唯象参数 $V_0$ 和 $\hbar\omega$ 的关系。一大类超导体相当好地遵循这些关系,但存在某些显著例外(如铅和汞)。即使这些称为“强耦合超导体”(strong coupling superconductors)的例外,也已被有说服力地纳入更一般的 BCS 理论框架,只要舍弃近似相互作用(34.16)所包含的简化,以及对声子效应的某些其他过于简单的表示。48
由模型哈密顿量(34.16),BCS 理论推导出以下主要平衡性质预言:
临界温度
在零磁场下,超导有序在下列临界温度开始形成:
$$ k_BT_c=1.13\hbar\omega e^{-1/N_0V_0}, \tag{34.17} $$
其中 $N_0$ 是正常金属中单一自旋电子群的电子能级密度49,$\omega$ 和 $V_0$ 为模型哈密顿量(34.16)的参数。由于存在指数依赖,有效耦合 $V_0$ 无法被确定得足够精确,以允许由式(34.17)非常准确地计算临界温度。不过,同样的指数依赖也解释了临界温度为何很低(通常比德拜温度低一至三个数量级):虽然 $\hbar\omega$ 为 $k_B\Theta_D$ 的量级,但对 $N_0V_0$ 的强烈依赖使得 $N_0V_0$ 在 0.1 至 0.5 范围内,即 $V_0n$ 在 $0.1\mathcal E_F$ 至 $0.5\mathcal E_F$ 范围内50,便可得到观测到的临界温度范围。还应注意,无论耦合 $V_0$ 多弱,理论都预言存在转变,尽管转变温度(34.17)可能低到无法观测。
能隙
对零温能隙,理论预言了与式(34.17)类似的公式:
$$ \Delta(0)=2\hbar\omega e^{-1/N_0V_0}. \tag{34.18} $$
式(34.18)与式(34.17)之比给出不依赖唯象参数的基本公式:
$$ \frac{\Delta(0)}{k_BT_c}=1.76. \tag{34.19} $$
对许多超导体,这一结果似乎在约 10% 的误差范围内成立(表 34.3)。不符合的超导体(如铅和汞,其偏差更接近 30%)往往也系统性地偏离简单理论的其他预言;采用强耦合理论的更精细分析,可以使其更接近理论预言。
初等理论还预言,在临界温度附近(零场中),能隙按照下列普适规律趋于零:51
$$ \frac{\Delta(T)}{\Delta(0)}=1.74\left(1-\frac T{T_c}\right)^{1/2},\qquad T\approx T_c. \tag{34.20} $$
表 34.3 $2\Delta(0)/k_BT_c$ 的实测值a
| 元素 | $2\Delta(0)/k_BT_c$ |
|---|---|
| Al | 3.4 |
| Cd | 3.2 |
| Hg ($\alpha$) | 4.6 |
| In | 3.6 |
| Nb | 3.8 |
| Pb | 4.3 |
| Sn | 3.5 |
| Ta | 3.6 |
| Tl | 3.6 |
| V | 3.4 |
| Zn | 3.2 |
a $\Delta(0)$ 取自隧穿实验。注意,BCS 理论给出的这一比值为 3.53。多数所列值具有 $\pm0.1$ 的不确定度。
来源:R. Meservey 和 B. B. Schwartz,Superconductivity,R. D. Parks 编,Dekker,New York,1969。
临界磁场
初等 BCS 理论对 $H_c(T)$ 的预言,常以相对于下列经验规律的偏差表示:52
$$ \frac{H_c(T)}{H_c(0)}\approx1-\left(\frac T{T_c}\right)^2. \tag{34.21} $$
图 34.8 显示若干超导体的量 $[H_c(T)/H_c(0)]-[1-(T/T_c)^2]$,并给出 BCS 预言。在所有情况下,偏差都很小,但应注意,强耦合超导体铅和汞比其他超导体偏离更大。
比热
在临界温度处(零磁场中),初等 BCS 理论预言比热存在不连续性,也可写成不依赖模型哈密顿量(34.16)中参数的形式:53
$$ \left.\frac{c_s-c_n}{c_n}\right|_{T_c}=1.43. \tag{34.22} $$
图 34.8 相对于粗略经验关系 $H_c(T)/H_c(0)\approx1-[T/T_c]^2$ 的偏差,包括若干金属的实测结果和简单 BCS 理论的预言。“强耦合”超导体铅和汞对简单 BCS 预言的偏离更显著。(J. C. Swihart 等,Phys. Rev. Lett. 14,106(1965)。)
[图内标签译文:Lead:铅;Mercury:汞;Indium:铟;Tin:锡;Aluminium:铝;BCS weak-coupling theory:BCS 弱耦合理论。]
除强耦合超导体外,这一预言与实验同样在约 10% 的误差范围内符合(表 34.4)。
低温电子比热也可写成与参数无关的形式:
$$ \frac{c_s}{\gamma T_c}=1.34\left(\frac{\Delta(0)}{k_BT}\right)^{3/2}e^{-\Delta(0)/(k_BT)}, \tag{34.23} $$
校对注:式(34.23)两处补上 $k_B$,沿用全章的能量与开尔文温度约定。
其中 $\gamma$ 是金属正常态比热中线性项的系数(式(2.80))。注意其指数下降,下降的尺度由能隙 $\Delta(0)$ 决定。
微观理论与迈斯纳效应
存在磁场时,无论金属处于正常态还是超导态,其平衡态中都会有抗磁电流流动,不过超导体中的电流要大得多。在自由电子模型中,电流由下列形式的方程确定,精确到场的一阶:54
$$ \nabla\times\boldsymbol{j}(\boldsymbol{r})=-\int d\boldsymbol{r}^{\prime}\,K(\boldsymbol{r}-\boldsymbol{r}^{\prime})\boldsymbol{B}(\boldsymbol{r}^{\prime}). \tag{34.24} $$
如果核 $K(\boldsymbol{r})$ 满足
$$ \int d\boldsymbol{r}\,K(\boldsymbol{r})=K_0\ne0, \tag{34.25} $$
那么,当磁场在 $K(\boldsymbol{r})$ 的作用范围内变化缓慢时,式(34.24)化为
$$ \nabla\times\boldsymbol{j}(\boldsymbol{r})=-K_0\boldsymbol{B}(\boldsymbol{r}), \tag{34.26} $$
表 34.4 比值 $[(c_s-c_n)/c_n]_{T_c}$ 的实测值a
| 元素 | $\left[\dfrac{c_s-c_n}{c_n}\right]_{T_c}$ |
|---|---|
| Al | 1.4 |
| Cd | 1.4 |
| Ga | 1.4 |
| Hg | 2.4 |
| In | 1.7 |
| La (HCP) | 1.5 |
| Nb | 1.9 |
| Pb | 2.7 |
| Sn | 1.6 |
| Ta | 1.6 |
| Tl | 1.5 |
| V | 1.5 |
| Zn | 1.3 |
a 简单 BCS 理论的预言为 $[(c_s-c_n)/c_n]_{T_c}=1.43$。
来源:R. Meservey 和 B. B. Schwartz,Superconductivity,R. D. Parks 编,Dekker,New York,1969。
这就是伦敦方程(34.7),其中 $n_s$ 为
$$ n_s=\frac{mc}{e^2}K_0. \tag{34.27} $$
由于伦敦方程意味着迈斯纳效应,因此在正常金属中,常数 $K_0$ 必须为零。要证明 BCS 理论意味着迈斯纳效应,就对外加场使用微扰理论计算核 $K(\boldsymbol{r})$,并明确验证 $K_0\ne0$。
实际证明 $K_0\ne0$ 是 BCS 理论的一项相当复杂的应用。不过,在伦敦方程刚提出时,伦敦兄弟便给出了更直观的解释。借助 V. L. Ginzburg 和 L. D. Landau 的唯象理论55,这一解释可以更具说服力。该理论虽比 BCS 理论早七年提出,却可很自然地用微观理论的一些基本概念描述。
金兹堡–朗道理论(Ginzburg–Landau theory)
金兹堡和朗道提出,超导态可由复“序参量”(order parameter)$\psi(\boldsymbol{r})$ 表征;它在 $T_c$ 以上为零,而其模在 $T_c$ 以下衡量位置 $\boldsymbol{r}$ 处超导有序的程度。56从 BCS 理论的观点看,序参量可视为描述库珀对质心位置的单粒子波函数。由于所有库珀对都处于同一个双电子态,一个函数就足够了。由于序参量不涉及对中两个电子的相对坐标,用 $\psi(\boldsymbol{r})$ 描述超导体仅适用于在电子对尺寸这一尺度57上变化缓慢的现象。
在超导体基态中,每一对都处于不依赖质心坐标的平移不变态;也就是说,序参量为常数。当有电流流动或出现外加场时,序参量便呈现出有趣的结构。金兹堡–朗道理论的一项基本假设是:在以矢势(vector potential)$\boldsymbol{A}(\boldsymbol{r})$ 给出的磁场中,以序参量 $\psi(\boldsymbol{r})$ 表征的超导体所承载的电流,由通常的量子力学电流公式给出,对应于波函数为 $\psi(\boldsymbol{r})$、电荷为 $-2e$、质量为 $2m$ 的粒子(即库珀对本身),也就是
$$ \boldsymbol{j}=-\frac e{2m}\left[\psi^*\left\{\left(\frac\hbar i\nabla+\frac{2e}c\boldsymbol{A}\right)\psi\right\}+\left\{\left(\frac\hbar i\nabla+\frac{2e}c\boldsymbol{A}\right)\psi\right\}^*\psi\right]. \tag{34.28} $$
若进一步假设序参量 $\psi=|\psi|e^{i\phi}$ 的显著空间变化体现在相位 $\phi$,而不是模 $|\psi|$ 上,则由式(34.28)得到伦敦方程(34.7)。由于序参量的模衡量超导有序的程度,这一假设将所考虑的扰动限制为:库珀对密度相对于其均匀热平衡值没有显著改变的扰动。在电子对可以流动、但不能积聚或被破坏的现象中,应满足这一条件。58
在这一假设下,电流关系式(34.28)简化为
$$ \boldsymbol{j}=-\left[\frac{2e^2}{mc}\boldsymbol{A}+\frac{e\hbar}m\nabla\phi\right]|\psi|^2. \tag{34.29} $$
由于任何梯度的旋度都为零,且 $|\psi|^2$ 基本上为常数,只要将超流体密度 $n_s$ 与 $2|\psi|^2$ 对应,便立即得到伦敦方程(34.7)。鉴于 $\psi$ 被解释为表征电荷为 $2e$ 的粒子的波函数,这一对应是合理的。
磁通量子化(flux quantization)
式(34.29)包含的结论比伦敦方程更为引人注目。考虑一个环形超导体(图 34.9)。若沿着位于超导材料深处、围绕环中孔洞的路径对式(34.29)
图 34.9 超导材料环,图中示出绕孔洞且位于超导体内部深处的路径。
积分,那么,由于显著电流只能在超导体表面附近流动,得到
$$ 0=\oint\boldsymbol{j}\cdot d\boldsymbol{\ell}=\oint\left(\frac{2e^2}{mc}\boldsymbol{A}+\frac{e\hbar}m\nabla\phi\right)\cdot d\boldsymbol{\ell}. \tag{34.30} $$
斯托克斯定理(Stokes theorem)给出
$$ \int\boldsymbol{A}\cdot d\boldsymbol{\ell}=\int\nabla\times\boldsymbol{A}\cdot d\boldsymbol{S}=\int\boldsymbol{B}\cdot d\boldsymbol{S}=\Phi, \tag{34.31} $$
其中 $\Phi$ 为环所围的磁通量。59此外,由于序参量是单值的,绕环一周时,其相位必须改变 $2\pi$ 乘以整数 $n$:
$$ \oint\nabla\phi\cdot d\boldsymbol{\ell}=\Delta\phi=2\pi n. \tag{34.32} $$
结合这些结果,得出环所围的磁通量必须量子化:
$$ |\Phi|=\frac{nhc}{2e}=n\Phi_0. \tag{34.33} $$
量 $\Phi_0=hc/2e=2.0679\times10^{-7}$ gauss-cm$^2$ 称为磁通子(fluxoid)或磁通量子(flux quantum)。磁通量子化已经被观测到,是以复序参量描述超导体有效性的最有力证据之一。60
微观理论与持续电流
遗憾的是,超导体得名的这一性质,是最难从微观理论中推导出来的性质之一。从某种意义上说,迈斯纳效应意味着完全导电性,因为在平衡态中,为了屏蔽宏观磁场,宏观电流必须无耗散地流动。事实上,直接从微观理论推导持续电流,与推导迈斯纳效应并无很大差别。计算电场诱导的电流至线性阶,并证明 AC 电导率中有一部分具有式(34.3)的形式,该形式适用于无耗散的电子气。为此,只需证明61
$$ \lim_{\omega\to0}\omega\operatorname{Im}\sigma(\omega)\ne0. \tag{34.34} $$
与式(34.3)比较,便可由这个非零常数的值确定超导电子密度 $n_s$ 的值。
证明式(34.34)成立,比证明迈斯纳效应的条件(34.25)更复杂,因为必须计入散射效应:若没有散射,任何金属都满足式(34.34),但即使没有散射,对正常金属抗磁性的计算也不会得到迈斯纳效应。不过,这一计算已经完成62,结果表明,由低频电导率推得的 $n_s$ 值与由迈斯纳效应计算推得的值一致。遗憾的是,这一计算相当形式化,无法直观解释一个引人注目的事实:一旦超导金属中建立起电流,所有熟悉的散射机制都不能有效地使其衰减。以下推理至少提示了这样一种直观解释:63
假设用电场在金属中建立电流,然后撤去电场,考察电流如何衰减。在正常金属中,电流可通过逐个电子的过程衰减;也就是说,单个电子与杂质、声子、缺陷等发生一系列碰撞,散射过程便可减小电子系统的总动量。平均而言,每次碰撞都使动量分布返回其平衡形式,在该形式下总电流为零。在超导体中建立电流时,所有库珀对一起运动:描述每一对的那个双电子态具有非零的质心动量。64可能会认为,这样的电流可通过单对碰撞而衰减,类似于正常金属中的单电子碰撞:各个电子对的质心动量通过碰撞降回零。然而,这种设想没有计及电子对之间微妙的相互依赖。65配对态的稳定性必不可少地依赖于这样一个事实:其他所有电子对都存在,并由相同的电子对波函数描述。因此,不能单独改变各个电子对的波函数而不彻底破坏配对态;彻底破坏配对态需要付出巨大的自由能代价。
使超电流衰减且自由能代价最小的转变,一般依赖于样品的几何形状,但通常要求破坏样品某个宏观部分的配对。这类过程是可能的,不过其自由能代价一般如此之大,以至于在任何实际的时间尺度上,超电流的寿命都是无限长的。66
超电流隧穿(supercurrent tunneling);约瑟夫森效应(Josephson effects)
我们已描述了单个电子从超导金属穿过薄绝缘势垒隧穿进入正常金属的过程(第 735 页),并说明了隧穿电流测量如何提供超导体单电子能级密度的信息。当两种金属都是超导体时,也可以测量隧穿电流;假设两种金属的单电子能级密度均具有 BCS 理论预言的形式,就能很好地拟合实验结果(图 34.10)。1962 年,Josephson67 预言,除了这种单电子的“正常隧穿”外,隧穿电流中还应有由配对电子承载的另一分量:只要势垒不太厚,电子对就能从一个超导体穿过结进入另一个,而不发生解离。
这一认识的直接结果之一是:在没有任何外加电场时,电子对的超电流应流过结(DC 约瑟夫森效应)。由于两个超导体之间只有弱耦合(即配对电子必须穿过非超导材料的间隙),流过结的典型隧穿电流将远小于单个样品的典型临界电流。
图 34.10 两个超导体(锡与铅)之间的正常隧穿电流。实线为 BCS 理论的预言。(S. Shapiro 等,IBM J. Res. Develop. 6,34(1962)。)
[图内标签译文:Current:电流;Voltage (millivolts):电压(毫伏);Calculated curve:计算曲线;Experimental points fitted at the point marked by the arrow:在箭头标示点处拟合的实验数据点。]
Josephson 假设结两侧的超导有序可以用单一序参量 $\psi(\boldsymbol{r})$ 描述,由此预言了其他一系列效应。他证明,隧穿电流由序参量跨越结时的相位变化决定。此外,利用规范不变性(gauge invariance)将序参量的相位与外加矢势的值联系起来,他能够证明,隧穿电流对结中存在的任何磁场都敏感地依赖。具体而言,存在磁场时,隧穿电流应具有下列形式:
$$ I=I_0\frac{\sin\pi\Phi/\Phi_0}{\pi\Phi/\Phi_0}, \tag{34.35} $$
其中 $\Phi$ 是结中的总磁通量,$\Phi_0$ 是磁通量子 $hc/2e$;$I_0$ 依赖于温度和结的结构,但不依赖于磁场。随后观测到了这些效应(图 34.11),从而进一步有力地证实了用序参量描述超导态在根本上的有效性,也证实了 Josephson 对理论富于想象力的应用。68
类似的考虑使 Josephson 进一步预言:若在这样的结两端施加 DC 电势,诱导的超电流将发生振荡(AC 约瑟夫森效应),其角频率为
$$ \omega_J=\frac{2eV}\hbar. \tag{34.36} $$
这一引人注目的结果——DC 电场应诱导交变电流——不仅被观测到,还成为高精度测量电压以及基本常数 $e/h$ 精确值的技术基础。69
图 34.11 Sn-SnO-Sn 结中,约瑟夫森隧穿电流随磁场的变化。(R. C. Jaklevic,引自 James E. Mercereau,Superconductivity,第 1 卷,R. D. Parks 编,Dekker,New York,1969,第 393 页。)
[图内标签译文:Josephson current:约瑟夫森电流;Magnetic field (gauss):磁场(高斯)。]
本书以这段粗略而又令人意犹未尽的超导概述结束,是恰当的。过去二十年间,为解释超导现象而成功发展起来的丰富且极富原创性的微观理论与唯象理论,显示出当代固体理论根基的稳健,以及其未来的前景。尽管超导理论所依赖的概念新颖,而且有时复杂得令人望而却步,仍不应忘记:它建立在广泛的基础之上,几乎涵盖了我们在前面各章考察过的固体理论的所有重要领域。在其他任何主题中,固体物理的两个基本分支——电子动力学与离子晶格振动——都没有如此紧密地融合,并产生如此引人注目的结果。
习题
1. 超导态的热力学
均匀磁场中超导体的平衡态,由温度 $T$ 和磁场大小 $H$ 决定。(假设压强 $P$ 固定,且超导体为平行于磁场的长圆柱,因此退磁效应不重要。)热力学恒等式用吉布斯自由能(Gibbs free energy)$G$ 表示更为方便:
$$ dG=-S\,dT-\mathfrak M\,dH \tag{34.37} $$
其中 $S$ 为熵,$\mathfrak M$ 为总磁化($\mathfrak M=MV$,其中 $M$ 为磁化密度)。$H$-$T$ 平面上,超导态与正常态之间的相界由临界磁场曲线 $H_c(T)$ 给出(图 34.3)。
(a) 根据 $G$ 在相界两侧连续这一事实,推导
$$ \frac{dH_c(T)}{dT}=\frac{S_n-S_s}{\mathfrak M_s-\mathfrak M_n}, \tag{34.38} $$
(其中下标 $s$ 和 $n$ 分别表示超导相和正常相中的值)。
(b) 利用超导态表现完全抗磁性($B=0$)、而正常态抗磁性可忽略($M\approx0$)这一事实,由式(34.38)证明跨越相界时的熵不连续量为
$$ S_n-S_s=-\frac V{4\pi}H_c\frac{dH_c}{dT}, \tag{34.39} $$
从而,在磁场中发生转变时,潜热为
$$ Q=-TV\frac{H_c}{4\pi}\frac{dH_c}{dT}. \tag{34.40} $$
(c) 证明,在零场下(即在临界点)发生转变时,存在下列比热不连续量:
$$ (c_p)_n-(c_p)_s=-\frac T{4\pi}\left(\frac{dH_c}{dT}\right)^2. \tag{34.41} $$
2. 超导平板的伦敦方程
考虑一个无限大超导平板,其边界为位于 $y=\pm d$ 处、垂直于 $y$ 轴的两个平行平面。沿 $z$ 轴施加磁场强度为 $H_0$ 的均匀磁场。
(a) 以 $\boldsymbol{B}$ 平行分量在表面连续为边界条件,由伦敦方程(34.7)和麦克斯韦方程(34.6)推导出超导体内
$$ \boldsymbol{B}=B(y)\boldsymbol{\hat z},\qquad B(y)=H_0\frac{\cosh(y/\Lambda)}{\cosh(d/\Lambda)}. \tag{34.42} $$
(b) 证明,在平衡态中流动的抗磁电流密度为
$$ \boldsymbol{j}=j(y)\boldsymbol{\hat x},\qquad j(y)=\frac c{4\pi\Lambda}H_0\frac{\sinh(y/\Lambda)}{\cosh(d/\Lambda)}. $$
(c) 平板内一点的磁化密度为 $\boldsymbol{M}(y)=(\boldsymbol{B}(y)-\boldsymbol{H}_0)/4\pi$。证明,平均磁化密度(在平板厚度上取平均)为
$$ \overline M=-\frac{H_0}{4\pi}\left(1-\frac\Lambda d\tanh\frac d\Lambda\right), \tag{34.43} $$
并给出平板很厚($d\gg\Lambda$)和很薄($d\ll\Lambda$)时磁化率的极限形式。
3. 圆柱导线中的临界电流
半径为 $r$ cm 的圆柱形超导导线中流过 $I$ 安培的电流。证明,当电流在紧邻导线外侧产生的磁场为 $H_c$(以 gauss 为单位)时,有
$$ I=5rH_c. \tag{34.44} $$
4. 库珀问题(Cooper problem)
考虑处于单重态的一对电子,其对称空间波函数为
$$ \phi(\boldsymbol{r}-\boldsymbol{r}^{\prime})=\int\frac{d\boldsymbol{k}}{(2\pi)^3}\chi(\boldsymbol{k})e^{i\boldsymbol{k}\cdot(\boldsymbol{r}-\boldsymbol{r}^{\prime})}. \tag{34.45} $$
在动量表象中,薛定谔方程具有下列形式:
$$ \left(E-2\frac{\hbar^2k^2}{2m}\right)\chi(\boldsymbol{k})=\int\frac{d\boldsymbol{k}^{\prime}}{(2\pi)^3}V(\boldsymbol{k},\boldsymbol{k}^{\prime})\chi(\boldsymbol{k}^{\prime}). \tag{34.46} $$
假设两个电子在简并自由电子气中相互作用;电子气的存在只通过不相容原理起作用:两个电子均不能占据 $k\lt k_F$ 的电子能级,由此得到约束:
$$ \chi(\boldsymbol{k})=0,\qquad k\lt k_F. \tag{34.47} $$
取电子对的相互作用为简单的吸引形式(参见式(34.16)):
$$ \begin{aligned}V(\boldsymbol{k}_1,\boldsymbol{k}_2)&\equiv-V,\qquad\mathcal E_F\le\frac{\hbar^2k_i^2}{2m}\le\mathcal E_F+\hbar\omega,\quad i=1,2;\\&=0,\qquad\text{其他情况,}\end{aligned} \tag{34.48} $$
并寻找与约束(34.47)相容的薛定谔方程(34.46)的束缚态解。由于我们只考虑在没有吸引作用时能量高于 $2\mathcal E_F$ 的单电子能级,束缚态的能量 $E$ 将低于 $2\mathcal E_F$,其结合能为
$$ \Delta=2\mathcal E_F-E. \tag{34.49} $$
(a) 证明,若满足下式,则存在能量为 $E$ 的束缚态:
$$ 1=V\int_{\mathcal E_F}^{\mathcal E_F+\hbar\omega}\frac{N(\mathcal E)\,d\mathcal E}{2\mathcal E-E}, \tag{34.50} $$
其中 $N(\mathcal E)$ 为给定自旋的单电子能级密度。
(b) 证明,只要 $N(\mathcal E_F)\ne0$,对于任意弱的 $V$,式(34.50)都存在 $E\lt 2\mathcal E_F$ 的解。(注意不相容原理所起的关键作用:若下限截断不是 $\mathcal E_F$,而是 0,那么由于 $N(0)=0$,对于任意弱耦合就不会存在解。)
(c) 假设在 $\mathcal E_F\lt \mathcal E\lt \mathcal E_F+\hbar\omega$ 范围内,$N(\mathcal E)$ 与 $N(\mathcal E_F)$ 的差别可忽略,证明结合能为
$$ \Delta=2\hbar\omega\frac{e^{-2/N(\mathcal E_F)V}}{1-e^{-2/N(\mathcal E_F)V}}, \tag{34.51} $$
或者,在弱耦合极限下为
$$ \Delta=2\hbar\omega e^{-2/N(\mathcal E_F)V}. \tag{34.52} $$
例如施加高压,或将样品制成非常薄的薄膜。铋提供了一个超导性以意想不到的方式得到增强的突出例子:非晶铋成为超导体的温度高于晶态铋,这在独立电子近似中完全无法理解。
见 B. W. Roberts,Progr. Cryog. 4,161(1964)。
纪录似乎为 $2\tfrac12$ 年;S. C. Collins,由 E. A. Lynton,Superconductivity,Wiley,New York,1969 引述。
在多种特殊条件下,超导也可以在没有能隙的情况下出现。例如,引入适当浓度的磁性杂质可以产生无能隙超导(gapless superconductivity)。综述见 K. Maki,载于 Superconductivity,R. D. Parks 编,Dekker,New York,1969。在超导语境中,“能隙”一词始终指量 $\Delta$。
这一现象在电子隧穿实验中观察得最为直接,下面将与能隙的其他表现一起介绍这些实验。
关于唯象理论的两部基础参考文献是 F. London,Superfluids,vol. 1,Wiley,New York,1954,以及 Dover,New York,1954;还有 D. Shoenberg,Superconductivity,Cambridge,1962。E. A. Lynton,Superconductivity,Methuen,London,1969 给出了非常简短的概述。微观理论见 J. R. Schrieffer,Superconductivity,W. A. Benjamin,New York,1964;以及 A. A. Abrikosov、L. P. Gorkov 和 I. E. Dzyaloshinski,Methods of Quantum Field Theory in Statistical Physics,Prentice-Hall,Englewood Cliffs,N.J.,1963 的最后一章。G. Rickayzen,Theory of Superconductivity,Interscience,New York,1965 对本领域的理论方面作了详细综述;P. de Gennes,Superconductivity of Metals and Alloys,W. A. Benjamin,Menlo Park,Calif.,1966 的综述稍简略一些。由本领域许多主要专家撰写、涵盖理论和实验各方面的综述著作为 Superconductivity,R. D. Parks 编,Dekker,New York,1969。
临界温度是没有外加磁场时发生转变的温度。存在磁场时(见下文),转变发生于更低温度,且转变性质从二阶变为一阶;即非零磁场下存在潜热。
截至目前,寻找超导性所达到的最低温度。
例如,$\mathcal E_F\sim10$ eV,$\hbar\omega_D\sim0.1$ eV。
假定不存在磁性杂质;见第 687 页。
当安培(Ampère)最初提出可以用单个分子内流动的电流来理解磁性时,有人反对,理由是尚无已知电流能无耗散地流动。安培坚持其观点,并得到量子理论的证实:量子理论允许存在净电流流动的分子定态(见第三十一章)。处于超导态的固体表现得如同一个巨大的分子。超导体中无耗散电流的存在,是量子力学的一种显著宏观表现。
见习题 3。
见第 253 页。
这一性质用于制造热开关。
声子传热的效能想必没有减弱,但声子对热导率的贡献通常不如传导电子重要。
正常金属仅具有微弱的顺磁性或抗磁性(没有磁有序金属是超导体),外加磁场可以穿透它。
正常金属特有的微小抗磁和顺磁效应除外。
习题 1 探讨这一行为的一些定量热力学后果。
不要与中间态(intermediate state)混淆。第一类超导体的形状比平行于磁场的圆柱更复杂时,可能形成中间态:宏观超导区域与正常区域交错排列,使磁场能的降低超过正常区域的自由能代价。
“涡旋”(vortex)一词既常用于指磁通丝本身,也用于指各磁通丝附近的电流结构。可以证明,每个涡旋包围的磁通恰好等于磁通量子(magnetic flux quantum)$hc/2e$(见脚注 60)。
磁场的存在不会显著影响正常态比热。
见第 727 页第 3 点。
例如两个样品表面的薄氧化层。
见第 360 页。
由于存在隧穿所需能量较少的热激发电子,阈值也会变得模糊。
见第 275–277 页。
F. London 和 H. London,Proc. Roy. Soc.(London),A149,71(1935),Physica 2,341(1935);F. London,Superfluids vol. 1,Wiley,New York,1954,以及 Dover,New York,1954。
双流体模型也用于描述超流氦-4;F. London 的 Superfluids 第 1、2 卷均有介绍,Ibid.。
本章始终忽略能带结构效应,用自由电子动力学描述电子。
假定随时间变化的速率足够慢,可以忽略位移电流。还取式(34.6)中的磁场为 $\boldsymbol{B}$ 而非 $\boldsymbol{H}$,因为 $\boldsymbol{j}$ 表示超导体内流动的平均微观电流。只有将 $\boldsymbol{j}$ 用满足 $\nabla\times\boldsymbol{M}=\boldsymbol{j}/c$ 的有效磁化密度表示,并按通常方式定义 $\boldsymbol{H}=\boldsymbol{B}-4\pi\boldsymbol{M}$ 时,才会出现磁场 $\boldsymbol{H}$。此时式(34.6)将被方程 $\nabla\times\boldsymbol{H}=0$ 替代。在这些 $\boldsymbol{H}$ 和 $\boldsymbol{M}$ 的定义下,这是完全等价的表述。
这是局域关系(local relation),即点 $\boldsymbol{r}$ 处的电流与同一点的磁场相关。A. B. Pippard 指出,更一般地,$\boldsymbol{r}$ 处的电流应由该点邻域内的磁场决定,关系具有形式
$$ \nabla\times\boldsymbol{j}(\boldsymbol{r})=-\int d\boldsymbol{r}'\,K(\boldsymbol{r}-\boldsymbol{r}')\boldsymbol{B}(\boldsymbol{r}'), $$
其中核 $K(\boldsymbol{r})$ 只有在 $r$ 小于长度 $\xi_0$ 时才显著。距离 $\xi_0$ 是表征超导体的若干基本长度之一;遗憾的是,这些长度都被不加区分地称为“相干长度”(coherence length)。在远低于临界温度的纯材料中,这些相干长度都相同,但在 $T_c$ 附近,或杂质平均自由程很短的材料中,“相干长度”可能随语境而不同。为避开这种混杂,仅讨论低温纯材料情况下它的意义,此时所有相干长度一致。在这种情况下,超导体属于第一类还是第二类的判据,是相干长度相对于伦敦穿透深度(London penetration depth)$\Lambda$(式(34.9))较大(第一类)还是较小(第二类)。
下面将看到,微观电子有序的某些特征也提示了伦敦方程。
例如,考虑占据半空间 $x\gt 0$ 的半无限超导体。式(34.8)意味着物理解呈指数衰减:
$$ B(x)=B(0)e^{-x/\Lambda}. $$
其他几何形状见习题 2。
J. Bardeen、L. N. Cooper 和 J. R. Schrieffer,Phys. Rev. 108,1175(1957)。这一理论通常称为 BCS 理论。
同位素效应(isotope effect)提供了离子运动参与建立超导性的直接证据:给定金属元素的不同同位素具有不同临界温度,通常(但并非总是)与离子质量的平方根成反比。只要存在对离子质量的依赖,就表明离子在转变中不能仅起静态作用,而必须动态参与。
见第 518–519 页。H. Fröhlich 最先强调,这种吸引可能存在,并可能是超导性的来源。
任何其他能使费米面附近电子之间产生净吸引相互作用的机制,也会在足够低温下导致超导态。不过,金属中尚无因其他机制产生超导性的实例得到有说服力的确证。
更一般地,可以考察 $n$ 个电子结合的可能性,但弱相互作用和泡利不相容原理使 $n=2$ 的情况最有希望。
L. N. Cooper,Phys. Rev. 104,1189(1956)。
假定额外电子的简并费米分布除禁止两个电子占据波矢小于 $k_F$ 的任何能级外,不起其他作用。因此,库珀计算基本上是双电子计算,只是分析被限制为由排除了所有波矢小于 $k_F$ 的平面波的单电子能级所构成的状态。见习题 4。
大系统极限下,未配对的那个电子(若 $N$ 为奇数)无关紧要。
反对称化算符
只是将自变量的所有可能置换产生的另外 $N!-1$ 个函数加到其作用的函数上,并根据置换由偶数次或奇数次成对交换构成,分别赋予权重 $+1$ 或 $-1$。
因此,一对费米子可以在统计上表现得像玻色子。实际上,如果每对的束缚能如此强,以致电子对尺寸远小于粒子间距 $r_s$,那么基态将由 $N/2$ 个玻色子构成,全部凝聚于同一个双电子能级。然而,下面将看到,库珀对(Cooper pair)的尺寸远大于 $r_s$,将库珀对看作独立玻色子可能极具误导性。
如果电子对态为三重态(自旋 1),将意味着出现某些未被观测到的特征磁性质。不过,液态氦-3 中已观测到三重态配对;它是与金属电子气有许多相似之处的简并费米液体。例如见 Nobel Symposium 24, Collective Properties of Physical Systems,B. Lundqvist 和 S. Lundqvist 编,Academic Press,New York,1973,第 84–120 页。
基态的这一方面常被强调,并被表述为自旋和波矢相反的电子束缚成对。这种表述的准确程度,不高于也不低于如下表述:任意两个相同粒子的平移不变束缚态都将它们配为动量大小相等、方向相反的对。即该表述正确地使人注意到束缚对的总动量为零,却误导性地使人忽略该态是这种对的叠加,因此在相对位置坐标上局域(与单个平面波乘积不同)。
在远低于 $T_c$ 的纯超导体中,这与脚注 31 描述的相干长度相同,因此用同一符号表示。
不能忘记,即使式(34.10)也只是对声子诱导电子间详细动态相互作用的相对粗略表示。在所谓强耦合超导体中(见下文),甚至式(34.10)也不充分。
强耦合超导理论处理完整的电子–声子系统,而不在一开始就试图消去声子、代以形式(34.16)甚至(34.10)的有效相互作用。因此,电子间净相互作用更复杂,不再是瞬时的,而是推迟的(retarded)。此外,与费米能级相差 $\hbar\omega_D$ 以内的电子能级,因电子–声子散射而产生的寿命可能很短,以致用明确的单电子能级构成电子对的图像也需要修改。
量 $N_0$ 就是 $g(\mathcal E_F)/2$。超导文献广泛采用这种能级密度记号。
量 $N_0$ 为 $n/\mathcal E_F$ 的量级。见式(2.65)。
式(34.20)是平均场理论的典型结果(参见平均场理论预言自发磁化按 $(T_c-T)^{1/2}$ 趋于零,第三十三章习题 6)。已知在足够接近临界温度的铁磁体中,平均场理论不正确。想必在超导体中足够接近 $T_c$ 时也会失效,但已有论证指出,平均场理论失效的区域极小(通常 $(T_c-T)/T_c\approx10^{-8}$)。超导体提供了相变的一个少见例子:即使非常接近临界点,平均场理论仍能很好地描述它。
BCS 预言同样可以写成与参数无关的形式。在低温下为 $H_c(T)/H_c(0)\approx1-1.06(T/T_c)^2$,而在 $T_c$ 附近为 $H_c(T)/H_c(0)=1.74[1-(T/T_c)]$。
$T_c$ 处的比热不连续性也是平均场理论的典型结果。想必非常接近 $T_c$ 时,比热可能发散。
这与脚注 31 中提到的核 $K$ 相同。
V. L. Ginzburg 和 L. D. Landau,Zh. Eksp. Teor. Fiz. 20,1064(1950)。
有时,将其与海森堡铁磁体类比会有所帮助;在后者中,序参量可视为局域自旋 $\boldsymbol{s}(\boldsymbol{r})$ 的平均值。在 $T_c$ 以上,$\boldsymbol{s}(\boldsymbol{r})$ 为零;在其以下,它给出自发磁化的局域值。在基态中,$\boldsymbol{s}(\boldsymbol{r})$ 与 $\boldsymbol{r}$ 无关(相应地,在不承载电流的均匀超导体中,$\psi(\boldsymbol{r})$ 为常数)。不过,可以考虑铁磁体更复杂的构型,例如外加场使棒两端的磁化指向不同方向。依赖于位置的 $\boldsymbol{s}(\boldsymbol{r})$ 在研究磁畴结构的某些方面时也很有用。同样,依赖于位置的 $\psi(\boldsymbol{r})$ 用于研究超导体承载电流的构型。
远低于 $T_c$ 时,这就是第 742 页描述的长度 $\xi_0$。
更一般地,当超导有序的程度确有显著空间变化时,必须将第二个金兹堡–朗道方程与式(34.28)联立,以确定 $\psi$ 和电流。第二个方程将序参量的空间变化率与矢势联系起来,形式上与单粒子薛定谔方程相似(这种相似在一定程度上会造成误导)。例如,描述第二类超导体中的涡旋时,必须使用完整的金兹堡–朗道方程组,因为在涡旋芯部,序参量的模迅速降至零,从而形成磁通量显著的区域。
由于磁场不能穿透超导材料,只要路径位于材料深处,所围磁通量就不依赖于路径的选择。
B. S. Deaver 和 W. M. Fairbank,Phys. Rev. Lett. 7,43(1961);R. Doll 和 M. Näbauer,Phys. Rev. Lett. 7,51(1961)。金兹堡–朗道理论还预言、且实验已经证实,第二类超导体中的每个涡旋包含一个磁通量子。
式(34.34)的结构与迈斯纳效应的条件(34.25)相似。核 $K$ 在整个空间上的积分,等于其空间傅里叶变换的 $k=0$ 极限。在这两种情况下,都必须证明某个电磁响应函数在适当的长波长或低频极限下不趋于零。
例如参见 A. A. Abrikosov、L. P. Gorkov 和 I. E. Dzyaloshinski,Methods of Quantum Field Theory in Statistical Physics,Prentice-Hall,Englewood Cliffs,N.J.,1963,第 334–341 页。
关于这一点,已有各种“直观”论证,其中许多并不成立。例如,有一种论证(基于朗道早先解释 $^4$He 超流动性的论证)声称可由单电子激发谱中存在能隙推导出持续电流。但这只解释了为何单电子激发不能使电流衰减,仍留下了电流通过逐对电子的过程衰减的可能性。这里指出的论证有多种表述形式,与“波函数的刚性”(rigidity of the wave function)、“非对角长程有序”(off-diagonal long-range order)或“长程相位相干”(long-range phase coherence)等概念相关。
超导体中不存在热电效应(见第 730 页),证实了将所有配对电子视为占据同一个量子态能够很好地描述超电流。若超电流类似于构成正常金属电流的电子无序流动,则会伴随热流(珀耳帖效应)。
回顾第 742 页,在给定电子对的半径范围内,会有数百万个其他电子对的中心。
因此,原则上,承载超电流的态只是亚稳态。对于适当的几何形状(即在一个或多个空间维度上非常小的样品),破坏超电流所需的涨落未必具有极低的发生概率,因此可以观测到“持续电流”的衰减。V. Ambegaokar 和 J. S. Langer 在 Phys. Rev. 164,498(1967)中给出了这类过程一个很有吸引力的微观图像。
B. D. Josephson,Phys. Lett. 1,251(1962)。另见 Josephson 和 Mercereau 在 Superconductivity,R. D. Parks 编,Dekker,New York,1969 中的文章。
磁通量子的值极小,因此这一效应作为一种测量磁场强度的高灵敏度方法,具有实际意义。
W. H. Parker 等,Phys. Rev. 177,639(1969)。