超流体与玻色-爱因斯坦凝聚:接近绝对零度时出现的宏观量子现象的物理学
在现代物理学中,接近绝对零度(0 K = -273.15 ℃)的极低温世界,是一个极大地颠覆了我们日常直觉的奇迹现象宝库。其中,“超流体(Superfluidity)”和“玻色-爱因斯坦凝聚(Bose-Einstein Condensation, BEC)”作为微观量子力学性质在宏观尺度上被直接观测到的“宏观量子现象”的代表,持续吸引着众多物理学家。
粘度完全变为零、自发爬上容器壁的液氦。受到光照时,液体从细管中喷射而出的热机械效应(喷泉效应)。还有,数万亿个原子落入完全相同的量子态,作为单一巨大物质波来运作的BEC。本文将从历史背景到高级理论框架,以极其详细且学术的视角深入探讨这些现象是如何被发现并在理论上被阐明的。
第一章:向绝对零度挑战与氦液化的历史
1.1 低温物理学的黎明与永久气体的液化
19世纪后半叶,物理学的一个重要主题是“是否所有气体都可以通过冷却和加压被液化”。通过迈克尔·法拉第等人的努力,许多气体被液化了,但氧、氮、氢和氦等被称为“永久气体”,其液化被认为极其困难。
然而,随着热力学的发展和利用焦耳-汤姆孙效应的冷却技术的进步,1877年氧气被液化,1898年詹姆斯·杜瓦实现了氢的液化(沸点约 20 K)。剩下的最后一种未被征服的气体,那就是氦。
1.2 氦的液化与卡末林·昂内斯的荣耀
在荷兰莱顿大学建立低温研究所的海克·卡末林·昂内斯(Heike Kamerlingh Onnes),进行了周密的准备并构建了巨大的液化装置。1908年7月10日,他使用液氢进行预冷,并进一步反复进行焦耳-汤姆孙膨胀,最终成功实现了氦的液化。液氦的沸点为 4.2 K,人类打开了距离绝对零度仅差4度的极低温世界的大门。由于这项成就,卡末林·昂内斯于1913年获得了诺贝尔物理学奖。
液氦的产生直接导致了后来1911年超导现象(水银的电阻在 4.2 K 时变为零的现象)的发现,但直到液氦本身的异常性质被揭示,又花了数十年的时间。
1.3 极化(Lambda)转变与超流体的发现
在1920年代到1930年代,威廉·亨德里克·基索姆等人发现,当液氦进一步冷却时,在约 2.17 K 处比热呈现对数发散。由于该比热随温度变化的图表类似于希腊字母“$\lambda$(Lambda)”,这种相变被命名为“$\lambda$(Lambda)转变”。高于 2.17 K($T_\lambda$)的相被称为“氦I”,低于此温度的相被称为“氦II”。
氦II表现出与普通液体完全不同的异常行为。最引人注目的发现是在1938年相继报道的。莫斯科的彼得·卡皮察(Pyotr Kapitza)以及剑桥大学的约翰·艾伦(John F. Allen)和唐·迈斯纳(Don Misener)独立发现,氦II可以无阻力(零粘度)地流过极细的毛细管。
卡皮察通过与超导(Superconductivity)的类比,将这种异常的流动性命名为“超流体(Superfluidity)”。在卡皮察进行的实验中,他将玻璃圆盘面对面放置,制造出微小的缝隙(几微米),并测量了液氦流过时的粘度。在氦I中观察到了作为普通粘性流体的阻力,但在低于 $T_\lambda$ 的瞬间,流出速度急剧增加,确认其表观粘性系数至少急剧下降到普通的 $10^4$ 分之一以下。
更奇妙的是,使用旋转粘度计进行的测量(在液体中旋转圆柱的实验)中,得到了矛盾的结果,氦II表现得似乎仍然具有有限的粘度。这种“在细管中粘度为零,但对旋转体施加粘性阻力”的看似不可思议的现象,后来被“双流体模型”完美解决。
第二章:玻色-爱因斯坦凝聚(BEC)的理论框架
从微观角度理解液氦超流体的关键在于量子统计力学。氦4原子由2个质子、2个中子和2个电子组成,其自旋总和为整数(0)。也就是说,氦4是“玻色子(Boson)”。当玻色子群体被冷却到绝对零度时会发生什么,我们将推导出其理论框架。
2.1 玻色统计与费米统计
在量子力学中,同类粒子在原理上是不可区分的。根据波函数对粒子交换的对称性,粒子大致分为两类。
- 玻色子(Boson):自旋为整数的粒子。对于粒子的交换,波函数是对称的(符号不改变)。不遵循泡利不相容原理,多个粒子可以同时处于同一个量子态。
- 费米子(Fermion):自旋为半整数的粒子。对于粒子的交换,波函数是反对称的(符号改变)。遵循泡利不相容原理,同一个量子态只能容纳一个粒子。
这里,为了使粒子数为非负,对于所有状态 $i$ 必须有 $\epsilon_i - \mu > 0$。如果设基态的能量为 $\epsilon_0 = 0$,则化学势必须始终满足 $\mu \leq 0$。
2.2 热德布罗意波长与宏观波函数的出现
$$ \lambda_{dB} = \sqrt{\frac{2\pi \hbar^2}{m k_B T}} $$在高温下 $\lambda_{dB}$ 非常短,远小于粒子间的平均距离 $d = (V/N)^{1/3}$。此时,粒子像经典的台球一样运动,用麦克斯韦-玻尔兹曼统计可以很好地近似。 然而,当温度 $T$ 降低时,$\lambda_{dB}$ 逐渐变长。然后,在某个临界温度 $T_c$ 处,波长 $\lambda_{dB}$ 变得与粒子间距离 $d$ 相当($\lambda_{dB} \sim d$)。此时,各个粒子的波函数在空间上开始重叠,量子力学干涉效应在宏观尺度上显现。
当温度进一步降低时,大量的玻色子会像雪崩一样落入最低能量状态(基态)。在这种状态下,无数个粒子将表现为具有完全相同相位的单一巨大物质波,即“宏观波函数(Macroscopic Wavefunction) $\Psi(\mathbf{r}) = \sqrt{n_0(\mathbf{r})} e^{i\theta(\mathbf{r})}$”。这就是玻色-爱因斯坦凝聚(BEC)。
2.3 凝聚转变温度 $T_c$ 的数学推导
$$ N = \sum_i \frac{1}{e^{\beta(\epsilon_i - \mu)} - 1} $$$$ D(\epsilon) = \frac{V}{(2\pi)^2} \left( \frac{2m}{\hbar^2} \right)^{3/2} \sqrt{\epsilon} $$$$ N_{ex} = \int_0^\infty D(\epsilon) \frac{1}{e^{\beta(\epsilon - \mu)} - 1} d\epsilon $$$$ N_{max} = \frac{V}{(2\pi)^2} \left( \frac{2m}{\hbar^2} \right)^{3/2} \int_0^\infty \frac{\sqrt{\epsilon}}{e^{\beta \epsilon} - 1} d\epsilon $$$$ \int_0^\infty \frac{\sqrt{\epsilon}}{e^{\beta \epsilon} - 1} d\epsilon = (k_B T)^{3/2} \int_0^\infty \frac{x^{1/2}}{e^x - 1} dx = (k_B T)^{3/2} \Gamma(3/2) \zeta(3/2) $$$$ N_{max} = V \left( \frac{m k_B T}{2\pi \hbar^2} \right)^{3/2} \zeta(3/2) = \frac{V}{\lambda_{dB}^3} \zeta(3/2) $$如果总粒子数 $N$ 超过 $N_{max}$ 会发生什么呢。多余的粒子 $N_0 = N - N_{max}$ 只能凝聚到不包含在积分($\epsilon > 0$)中的能量 $\epsilon = 0$ 的基态中。这就是 BEC 的数学机制。
凝聚开始的临界温度 $T_c$ 定义为使得 $N = N_{max}(T_c)$ 成立的温度。解出 $T_c$: $$ T_c = \frac{2\pi \hbar^2}{m k_B} \left( \frac{n}{\zeta(3/2)} \right)^{2/3} $$ (其中 $n = N/V$ 是数密度) 该方程给出了无相互作用的理想气体中严格的 BEC 转变温度。代入液氦的密度,得到 $T_c \approx 3.1 K$,这接近实际的 $\lambda$转变温度 2.17 K。这种一致性强烈暗示了氦II的超流性在本质上是与 BEC 相关的现象。然而,因为液氦是粒子间相互作用非常强的“强关联量子液体”,所以与理想气体模型存在偏差。
第三章:蒂萨和朗道的双流体模型
为了解释超流氦表现出的“在细管中无粘度流动,但在旋转体上显示粘性阻力”的奇妙性质,1938年拉斯洛·蒂萨(László Tisza)提出了突破性的唯象理论。之后,在1941年列夫·朗道(Lev Landau)基于量子力学的微观基础将其完善为更严谨的理论,这就是“双流体模型(Two-Fluid Model)”。
3.1 正常流体成分与超流体成分的共存
$$ \rho = \rho_n + \rho_s $$$$ \mathbf{j} = \rho_n \mathbf{v}_n + \rho_s \mathbf{v}_s $$这里:
- 正常流体成分(Normal fluid component, $\rho_n$):具有有限粘度、携带熵的成分。对应于处于热激发态的粒子群。
- 超流体成分(Superfluid component, $\rho_s$):粘度完全为零、不携带熵(与绝对零度状态相同)的成分。对应于处于宏观量子态(BEC状态)的粒子群。
在绝对零度($T = 0$)时,全部为超流体成分($\rho_s = \rho, \rho_n = 0$),但随着温度升高,正常流体成分的比例增加,在 $T_\lambda$ 处超流体成分消失($\rho_s = 0, \rho_n = \rho$)。
利用该模型,先前提到的粘度矛盾得到了完美的解决。在穿过极细管的实验中,有粘度的正常流体成分因与管壁的摩擦而被阻挡,但由于只有无粘度的超流体成分穿过管子,因此表观粘度变为零。另一方面,在体液中旋转圆柱的实验中,由于圆柱与正常流体成分产生摩擦,因此观察到有限的粘度。
3.2 热机械效应(喷泉效应)与力卡效应(Mechanocaloric Effect)
超流体成分“不携带熵”的性质引起了热与物质流动直接耦合的惊人现象。其代表就是“热机械效应(Thermomechanical effect)”或“喷泉效应(Fountain effect)”。
$$ \Delta P = \rho S \Delta T $$($S$ 是单位质量的熵)
反之,如果通过超漏将氦II压出,由于流出的仅仅是不携带熵的超流体成分,留下的液体的温度就会上升。这被称为力卡效应(Mechanocaloric effect)。
3.3 第二声(Second Sound)的物理学
在普通流体中,“声波”是密度的疏密波(压力波)。然而,在氦II中,由于存在两个独立的速度场 $\mathbf{v}_n$ 和 $\mathbf{v}_s$,因此可以传播两种类型的波。
- 第一声(First Sound):正常流体成分和超流体成分以同相位振动的波。相当于普通的密度波(压力波)。
- 第二声(Second Sound):正常流体成分和超流体成分以反相位(相互反向)振动的波。在保持密度不变的情况下,正常流体成分(熵的集合)和超流体成分的比例在空间上振荡,因此这实际上作为“温度的波(热波)”来传播。
在普通物质中,热作为扩散现象传播,但在氦II中,热作为无衰减的“波”以声速(第二声速)传播。这一现象也由朗道预言,并在后来被实验证实,从而决定性地证明了双流体模型的正确性。
第四章:朗道的超流体临界速度与元激发光谱
由于玻色-爱因斯坦凝聚理论是以理想气体为前提的,它无法解释氦原子之间强烈的相互作用。在没有相互作用的理想玻色气体中,实际上不会出现超流性(很容易因杂质散射而被激发)。列夫·朗道将液氦公式化为“处于绝对零度基态的背景液体”和叠加其上的“元激发(Elementary Excitations)气体”。
4.1 声子与旋子:能量-动量色散曲线
朗道为了解释氦II的热学性质,假定了元激发的能量 $\epsilon$ 和动量 $p$ 的关系(色散关系 $\epsilon(p)$)如下。
- $$ \epsilon(p) = c p $$
($c$ 为第一声速)
- $$ \epsilon(p) = \Delta + \frac{(p - p_0)^2}{2\mu} $$
其中,$\Delta$ 是旋子能隙(Roton gap energy),$p_0$ 是取得极小值时的动量,$\mu$ 是旋子的有效质量。
这条奇妙的色散曲线,后来通过中子非弹性散射实验被极其精确地实测出来,完美地证明了朗道天才的直觉。
4.2 通过临界速度 $v_c$ 证明无摩擦散射不可能
朗道最大的功绩在于利用这个色散关系,在力学上证明了“为什么氦II能够在没有与壁面摩擦的情况下流动”这一超流体的根本机制(朗道临界速度定理)。
$$ \Delta E = \mathbf{v} \cdot \mathbf{p} = \epsilon(p) $$$$ v \geq \frac{\epsilon(p)}{p} $$$$ v_c = \min \left[ \frac{\epsilon(p)}{p} \right] $$如果液氦是理想玻色气体,它将具有自由粒子的色散关系 $\epsilon(p) = p^2 / 2m$,因此 $v_c = \min [p/2m] = 0$,无论移动得多么慢,都会发生激发,不会成为超流体。 然而,在真实的液氦中,由于相互作用,色散曲线从线性($\epsilon = c p$)开始,所以在原点附近的斜率为 $c$(声速,约 240 m/s)。从整个色散曲线的形状求取最小值,指向旋子极小值的切线斜率给出了 $v_c$,其值约为 60 m/s。 (※实际实验中观测到的临界速度要小得多,大约只有几 cm/s,这将在后来由理查德·费曼等人作为“量子涡旋”的产生机制而被阐明。)
这样,系统的能量色散关系中存在“能隙”或“有限斜率”,是保证超流性(零摩擦流动)的绝对条件。
第五章:量子涡旋(Quantized Vortex)与拓扑缺陷
虽然朗道的临界速度理论解释了超流体的微观起源,但实验中观测到的临界速度远低于理论值的原因仍是一个谜。拉斯·昂萨格(Lars Onsager)和理查德·费曼(Richard Feynman)的“量子涡旋”理论解开了这个谜团,并将其与超流体宏观波函数的相位(拓扑学)联系起来。
5.1 宏观波函数的相位与环量的量子化
$$ \mathbf{v}_s = \frac{\hbar}{m} \nabla \theta $$$$ \kappa = \oint_C \mathbf{v}_s \cdot d\mathbf{r} = \frac{\hbar}{m} \oint_C \nabla \theta \cdot d\mathbf{r} $$$$ \Delta\theta = 2\pi n \quad (n = 0, \pm 1, \pm 2, \dots) $$$$ \kappa = \frac{\hbar}{m} (2\pi n) = n \frac{h}{m} $$也就是说,超流体中的环量不能取连续值,而是以 $h/m$(普朗克常数 $h$ 除以氦原子质量 $m$ 的值)为单位完全“量子化”了。这就是“环量的量子化”。
5.2 旋转液氦中的涡旋晶格与拓扑缺陷
当旋转装有普通粘性流体的容器时,流体会像刚体一样整体旋转,并描绘出一个抛物面。但是,超流氦(的超流体成分)由于具有 $\nabla \times \mathbf{v}_s = (\hbar/m) \nabla \times \nabla \theta = 0$,即无旋(无涡)的性质,即使慢慢旋转容器,里面的物质也不会一起旋转。 然而,当逐渐提高旋转速度时,在超过某个临界角速度的瞬间,超流体中会侵入无数条“量子涡旋(Quantized Vortex)”线(涡丝)。
在量子涡旋的中心(核心)处,由于相位成为奇点而无法定义,因此超流密度 $n_0$ 变为零(即变为正常流体或真空)。流体以量子化的环量 $\kappa = h/m$ (通常 $n=1$ 的最小状态在能量上最稳定)绕核心旋转。 如果进一步提高旋转速度,量子涡旋的数量增加,它们相互排斥,形成规则的三角晶格(类似于阿布里科索夫晶格的涡旋晶格)。从宏观上看,这无数个量子涡旋的集合体,仿佛整个流体都在进行刚体旋转。
朗道临界速度的不一致也由这个量子涡旋解释了。在流经毛细管壁时,即使流速非常小,壁上的微小凹凸也会产生量子涡旋。由于产生这种量子涡旋所需的能量远低于旋子激发,因此实验中的临界速度极大地小于朗道的理论值。量子涡旋作为破坏超流性的“拓扑缺陷”发挥作用。
第六章:激光冷却实现气体BEC及向现代物理学的扩展
由于液氦中粒子间相互作用极强,即使在绝对零度,凝聚率也只停留在约 10% 左右,距离纯粹理想气体的 BEC 相去甚远。为了实现真正意义上的“稀薄气体中的 BEC”,需要将原子间相互作用可忽略不计的稀薄气体冷却到接近绝对零度。
6.1 碱金属原子气体的 BEC 达成
在20世纪80年代到90年代期间,开发出了激光冷却(多普勒冷却、西西弗斯冷却)和磁捕集,以及蒸发冷却等创新技术,使得将原子气体冷却到纳开尔文($10^{-9}$ K)区域成为可能。 然后,在1995年,科罗拉多大学(JILA)的埃里克·康奈尔(Eric Cornell)和卡尔·威曼(Carl Wieman)使用铷($^{87}$Rb)原子气体,同年,麻省理工学院的沃尔夫冈·克特勒(Wolfgang Ketterle)使用钠($^{23}$Na)原子气体,在人类历史上首次成功直接观测到了稀薄碱金属原子气体中的玻色-爱因斯坦凝聚。他们凭借这一成就于2001年获得了诺贝尔物理学奖。
克特勒团队进行的实验中,通过让两个独立的 BEC 在空间上干涉,观测到了清晰的“干涉条纹”。这是证明宏观尺度的物质块作为单一物质波在运作的决定性证据。
6.2 光晶格中的超流体-莫特绝缘体转变
现代的 BEC 研究已经超越了单纯的现象确认,向作为“量子模拟器”的应用发展。当把 BEC 引入由相对射的激光束干涉产生的周期性势场(光晶格,Optical Lattice)中时,原子的行为就像固体晶格中的电子一样。 通过调整激光的强度,可以自由控制原子间的相互作用和跃迁(向相邻格点的隧道效应)的比例。当相互作用很弱时,原子在整个系统中处于相位一致的状态下自由运动,形成“超流体状态”;但如果增强相互作用,每个格点上将局域化一个原子,宏观相位相干性被完全破坏,发生向“莫特绝缘体(Mott Insulator)状态”的量子相变。 这种超流体-莫特绝缘体转变(由 Greiner 等人在2002年观测到)为在理想的人工系统中精确验证强关联电子系统(如高温超导体等)的哈伯德模型开辟了道路。
6.3 与超导(库珀对的 BEC)的统一理解
最后,让我们提一下超流体与超导之间的深刻联系。超导现象可以解释为金属中的传导电子(费米子)通过声子(晶格振动)媒介的引力,两两结对(库珀对),并且这些对像玻色子一样运作从而发生 BEC 的现象(BCS理论)。 近年来,使用费米原子气体,通过磁场(费什巴赫共振)控制原子间的相互作用,不断有研究在进行将分子的 BEC(BEC 极限)连续转变为库珀对的 BCS 状态(BCS 极限)的“BCS-BEC 交叉(crossover)”研究。
超流体和 BEC,是微观量子力学支配的法则扩展到日常宏观尺度上的终极物理现象。零粘度的流动、环量的量子化、物质波的干涉。这些全都是自然界在绝对零度这种极限状态下向我们展示的,宇宙中最根本和纯粹的姿态。
