全同粒子、二次量子化、Bose-Einstein凝聚、超导与超流的量子多体理论
多体量子系统（Many-Body Quantum System）是指由大量微观粒子（电子、原子、分子等）组成的、必须用量子力学描述的物理系统。与单体系统不同，多体系统中粒子之间的相互作用以及全同粒子的不可区分性，导致了一系列涌现现象（emergent phenomena），如超导、超流、Bose-Einstein凝聚、量子磁性、量子自旋液体等。
多体量子理论是凝聚态物理、原子分子物理、核物理和量子化学的共同理论基础。从日常生活中的磁铁到超导量子计算机，从低温液氦的超流现象到中子星内部的超流态，从高温超导电缆到光学晶格中的量子模拟器，多体量子系统无处不在。
多体问题之所以困难，在于其自由度随粒子数指数增长。即使只有几十个电子，完整的量子力学波函数也无法用经典计算机对角化。因此，多体物理的核心任务是发展有效的近似方法和概念框架，以理解大量粒子共同行为中涌现的规律。
量子力学中最深刻的原理之一：全同粒子（identical particles）在本质上是不可区分的。在经典力学中，即使两个粒子完全相同（具有相同的质量和电荷），我们仍然可以通过追踪它们的运动轨迹来区分它们。但在量子力学中，由于不确定性原理，粒子的轨迹不再是定义良好的物理量。
具体来说，考虑一个包含两个全同粒子的系统（例如两个电子），其波函数为 ψ(r1,s1,r2,s2)，其中 ri 和 si 分别表示第 i 个粒子的空间坐标和自旋量子数。交换两个粒子后，概率分布必须保持不变：
∣ψ(1,2)∣2=∣ψ(2,1)∣2
这意味着交换操作的效果只能是一个相位因子：
ψ(2,1)=eiθψ(1,2)
在三维空间中，两次交换后回到原状态，因此 e2iθ=1，所以 eiθ=±1。因此：
ψ(q1,q2)=±ψ(q2,q1)
按照交换对称性，自然界中所有基本粒子可分为两大类：
| 粒子类型 |
交换对称性 |
自旋统计 |
示例 |
统计分布 |
| 玻色子 (Bosons) |
对称：(+) 交换波函数不变 |
整数自旋 (0,1,2,…) |
光子、胶子、声子、4He原子、W/Z玻色子、介子 |
Bose-Einstein分布 |
| 费米子 (Fermions) |
反对称：(−) 交换波函数变号 |
半整数自旋 (1/2,3/2,…) |
电子、质子、中子、中微子、夸克、3He原子 |
Fermi-Dirac分布 |
泡利不相容原理：两个全同费米子不能占据完全相同的量子态。如果两个费米子的自旋不同（一上一下），它们可以占据相同的空间轨道，但如果自旋相同，则必须占据不同的轨道。泡利原理是原子壳层结构、化学键形成、固体能带论等一系列物理现象的根源。
玻色子的"抱团"特性：由于玻色子倾向于占据相同的量子态，当温度降低到临界值以下时，大量玻色子会"凝聚"到最低能量状态，形成Bose-Einstein凝聚。
对于处于热平衡的粒子系统，量子统计分布由温度和化学势决定：
Bose-Einstein分布（玻色子）：
⟨ni⟩=e(ϵi−μ)/(kBT)−11
Fermi-Dirac分布（费米子）：
⟨ni⟩=e(ϵi−μ)/(kBT)+11
经典极限：当占据数很低时（e(ϵi−μ)/(kBT)≫1），两者都趋近于经典的Maxwell-Boltzmann分布：
⟨ni⟩≈e−(ϵi−μ)/(kBT)
数值对比：假设 ϵk=0.1 eV，μ=0，kBT=0.025 eV（室温），则：
| 分布类型 |
公式 |
占据数 ⟨n⟩ |
| Maxwell-Boltzmann |
e−ϵ/(kBT) |
e−4≈0.0183 |
| Bose-Einstein |
[eϵ/(kBT)−1]−1 |
[e4−1]−1≈0.0186 |
| Fermi-Dirac |
[eϵ/(kBT)+1]−1 |
[e4+1]−1≈0.0179 |
在高温低密度下三种分布几乎一致，但在低温高密度下差异巨大。例如当 ϵ=0，μ→0− 时，Bose-Einstein分布会出现发散——这正是BEC相变的信号。而Fermi-Dirac分布的指数上为正号，这意味着费米子的占据数永远不会超过1。
以氦同位素为例，直观展示自旋统计对宏观物理性质的深刻影响：
- 4He：原子核含2个质子和2个中子，总自旋为0 → 玻色子
- 3He：原子核含2个质子和1个中子，总自旋为1/2 → 费米子
虽然化学性质几乎相同（都是惰性气体），但它们的低温行为截然不同：
| 物理性质 |
4He (玻色子) |
3He (费米子) |
| 大气自然丰度 |
99.99986% |
0.00014% |
| 沸点 (1 atm) |
4.22 K |
3.19 K |
| 超流转变温度 |
2.17 K (λ 点) |
约 0.93 mK (需Pauli配对) |
| 超流机制 |
BEC型直接凝聚 |
类似BCS理论的配对凝聚 |
| 0 K时的基态 |
全部凝聚在基态 |
形成Fermi海，费米能级 EF |
4He在常压下降到2.17 K以下时，会自发进入超流态——液体可以无摩擦地流动，甚至沿容器壁爬升并自动流出容器（著名的"爬壁效应"）。而3He必须冷却到约1 mK（比4He低2000倍）才能通过类似BCS的配对机制进入超流态。这一巨大差异完全源自两种原子的量子统计性质不同。
进一步有趣的是：4He的λ相变（比热在Tλ处呈尖锐的λ形）是统计物理中第二类相变的经典范例，其临界指数（如 α≈−0.013）与三维XY模型完全吻合。
描述多体系统时，通常的坐标空间波函数 Ψ(r1,r2,…,rN) 面临两个严重问题：
- 粒子数可变：许多重要物理过程涉及粒子数变化（光吸收/发射、电子-空穴对产生、声子产生/湮灭），但坐标波函数形式要求固定的 N。
- 对称性约束操作复杂：对 N 个全同粒子，波函数必须满足对称/反对称条件，这在 N 增大时操作复杂度爆炸式增长。例如，Slater行列式的展开项数为 N!，N=10 时已有约 3.6×106 项。
二次量子化（Second Quantization） 以场算符取代多粒子波函数，将粒子数变化和对称性约束自动纳入形式体系，是现代多体物理的"母语"。
引入产生算符 a^i† 和 湮灭算符 a^i，其中下标 i 标记单粒子量子态（包括动量、自旋、能级等）。
[a^i,a^j†]=δij,[a^i,a^j]=[a^i†,a^j†]=0
{c^i,c^j†}=δij,{c^i,c^j}={c^i†,c^j†}=0
其中 [A^,B^]=A^B^−B^A^，{A^,B^}=A^B^+B^A^。
重要区别：玻色子的对易关系允许多个粒子占据同一态（a^i†∣n⟩=n+1∣n+1⟩），而费米子的反对易关系导致 Pauli原理（c^i†c^i†=0，不能两次产生同一费米子）。
单粒子算符（动能、外势能等）：
O^1=i,j∑⟨i∣o^∣j⟩a^i†a^j
其中矩阵元 ⟨i∣o^∣j⟩=∫ϕi∗(x)o^(x)ϕj(x)dx。
两粒子算符（库仑相互作用等）：
O^2=21i,j,k,l∑⟨ij∣v^∣kl⟩a^i†a^j†a^la^k
其中 ⟨ij∣v^∣kl⟩=∬ϕi∗(x1)ϕj∗(x2)v(x1−x2)ϕk(x2)ϕl(x1)dx1dx2。
注意算符的排列顺序：a^i†a^j†a^la^k 确保了费米子的反对称性（交换产生算符会引入负号）。
引入连续坐标空间中的场算符：
ψ^(x)=i∑ϕi(x)a^i,ψ^†(x)=i∑ϕi∗(x)a^i†
场算符满足：
[ψ^(x),ψ^†(x′)]∓=δ(x−x′)
其中 − 号对应玻色子，+ 号对应费米子。
用场算符重新表达哈密顿量：
H^=∫dxψ^†(x)[−2mℏ2∇2+V(x)]ψ^(x)+21∬dxdx′ψ^†(x)ψ^†(x′)v(x−x′)ψ^(x′)ψ^(x)
考虑金属中自由电子的哈密顿量，这是理解固体电子结构的起点：
H^=H^0+H^int
动能部分：
H^0=k,σ∑ϵkc^k,σ†c^k,σ
其中 ϵk=ℏ2k2/(2m)，σ=↑,↓。
库仑相互作用部分：
H^int=2V1k,k′,q=0∑σ,σ′∑q24πe2c^k+q,σ†c^k′−q,σ′†c^k′,σ′c^k,σ
符号含义：
| 符号 |
含义 |
典型值 |
| ϵk=ℏ2k2/(2m) |
电子动能 |
EF∼1-10 eV |
| V |
系统体积 |
宏观量 |
| 4πe2/q2 |
库仑作用的Fourier变换 |
在 q→0 时发散 |
| q=0 |
排除了均匀背景的自相互作用 |
Hartree项的抵消 |
关键观察：q=0 项对应均匀的电荷密度，被正离子背景抵消，所以求和排除 q=0。这就是著名的电子气模型。
平均场理论（Mean-Field Theory, MFT）是多体物理中最基本且最广泛使用的近似方法。其核心思想：将多体问题简化为单体在有效场中的运动，其中有效场由所有其他粒子的平均效应构成。
这种方法在各种相互作用系统中都能使用。虽然平均场忽略了涨落，但因为它能捕捉对称性破缺和相变的基本特征，通常作为理解复杂系统的第一步。
对于相互作用的费米子系统，Hartree-Fock方法将多体波函数近似为单个Slater行列式。在二次量子化语言中，这相当于在相互作用项中做"收缩"（Wick定理）：
c^i†c^j†c^lc^k≈⟨c^i†c^k⟩c^j†c^l+⟨c^j†c^l⟩c^i†c^k−⟨c^i†c^l⟩c^j†c^k−⟨c^j†c^k⟩c^i†c^l
前两项是Hartree项（直接库仑作用），后两项是Fock项（交换作用，带有负号）。
有效单粒子哈密顿量：
H^HF=k,σ∑(ϵk+ΣH(k)+ΣF(k))c^kσ†c^kσ
其中 Hartree自能 ΣH 和 Fock自能 ΣF 需要自洽求解。
以Hubbard模型（描述格点上的相互作用电子）为例，说明自洽循环：
H^=−t⟨i,j⟩,σ∑c^iσ†c^jσ+Ui∑n^i↑n^i↓
步骤1：给定初始猜测。例如假设 ⟨ni↑⟩=0.5，⟨ni↓⟩=0.5。
步骤2：计算有效势。在平均场近似下，Uni↑ni↓≈U(⟨ni↑⟩ni↓+ni↑⟨ni↓⟩−⟨ni↑⟩⟨ni↓⟩)，所以自旋向上和向下的有效势不同：
Veff,σ(i)=U⟨ni,−σ⟩
步骤3：解有效哈密顿量，得到新的 ⟨niσ⟩。
步骤4：迭代至收敛。通常要求相邻两次迭代的密度矩阵之差小于 10−6。
数值示例：一维4格点链，U/t=4，半填充：
| 迭代次数 |
⟨n1↑⟩ |
⟨n1↓⟩ |
⟨n2↑⟩ |
⟨n2↓⟩ |
| 0 (初猜) |
0.50 |
0.50 |
0.50 |
0.50 |
| 1 |
0.62 |
0.38 |
0.62 |
0.38 |
| 2 |
0.68 |
0.32 |
0.68 |
0.32 |
| 3 |
0.71 |
0.29 |
0.71 |
0.29 |
| 4 |
0.72 |
0.28 |
0.72 |
0.28 |
| 5 (收敛) |
0.73 |
0.27 |
0.73 |
0.27 |
系统出现了反铁磁序——自旋向上和向下在同一个格点上的占据数不再相等，且相邻格点的磁矩方向相反。这是一种自发对称性破缺：哈密顿量本身是自旋旋转对称的，但基态自发选择了特定方向。
| 优点 |
局限性 |
| 计算简单，数学上可解析处理 |
完全忽略量子涨落和热涨落 |
| 适用弱到中等相互作用强度 |
在临界点（Tc）附近严重不准确 |
| 能捕捉对称性破缺 |
无法描述的强关联现象（如Mott绝缘体） |
| 是许多高级方法的出发点 |
低估有效质量和能隙值 |
| 对长程系统较准确（涨落被平均掉） |
在低维系统中失效（涨落增强） |
实例：Ising模型的平均场临界温度
Ising模型的平均场近似给出：
TcMF=zJ
其中 z 是配位数，J 是交换耦合。与精确解的对比：
| 维度 |
z |
TcMF/J |
精确值 Tc/J |
误差 |
| 1D |
2 |
2.0 |
0 |
∞（不存在相变） |
| 2D（正方晶格） |
4 |
4.0 |
2.269 |
76% |
| 3D（立方晶格） |
6 |
6.0 |
4.511 |
33% |
| 无穷维（团簇） |
∞ |
∞ |
外推 Tc∝z |
趋于准确 |
Mermin-Wagner定理指出：在二维或更低维的短程相互作用系统中，连续对称性在有限温度下不会自发破缺。这意味着平均场在低维不准确不仅是一个定量问题，而是一个定性错误。
Bose-Einstein凝聚（BEC）是指大量玻色子占据同一单粒子基态的宏观量子现象。1924年，Satyendra Nath Bose首先推导了光子的统计分布，随后爱因斯坦将这一统计推广到物质粒子，并预言了相变。
关键条件：在绝对温度 T 下，德布罗意热波长为：
λdB=2πmkBTh
当 λdB 与粒子间的平均距离 n−1/3 相当时，量子统计效应变得显著。标准是：
nλdB3⪆2.612
这个无量纲参数 nλdB3 称为相空间密度。
对于自由玻色气体，临界温度由相空间密度条件确定：
Tc=mkB2πℏ2(ζ(3/2)n)2/3
其中 ζ(3/2)≈2.612。
具体数值计算：以87Rb原子为例
| 物理量 |
符号 |
数值 |
| 原子质量 |
m |
1.443×10−25 kg |
| 典型数密度 |
n |
1020 m−3 |
| 热波长 |
λdB=h/2πmkBT |
在 Tc 时为约 0.29 μm |
| 平均间距 |
n−1/3 |
约 0.21 μm |
| Tc |
公式计算 |
约 276 nK |
这就是为什么BEC只能在纳开尔文（10−9 K）量级下实现——原子质量远比电子大，因此临界温度极低。对于光子（质量为零），临界温度可高得多。
在BEC中，引入宏观波函数（即序参量）：
Ψ(r)=n0(r)eiθ(r)
其中 n0(r)=∣Ψ(r)∣2 是凝聚体的局域密度，θ(r) 是相位。宏观波函数的出现代表了U(1)规范对称性的自发破缺。
凝聚体的动力学由Gross-Pitaevskii方程（GPE）描述：
iℏ∂t∂Ψ(r,t)=[−2mℏ2∇2+Vext(r)+g∣Ψ(r,t)∣2]Ψ(r,t)
其中 g=4πℏ2a/m 是有效相互作用强度，a 是s波散射长度。Vext 是外势（如磁阱或光阱）。
Thomas-Fermi近似：当相互作用占主导时，可以忽略动能项，得到凝聚体的密度分布：
n0(r)=gμ−Vext(r)
其中 μ 是化学势。对于谐振子势 Vext=21m(ωx2x2+ωy2y2+ωz2z2)，密度分布呈现倒抛物面形状：
n0(r)=n0(0)(1−Rx2x2−Ry2y2−Rz2z2)
其中 Ri=2μ/(mωi2) 是Thomas-Fermi半径。
1995年，JILA的Cornell和Wieman团队首次在87Rb中实现了BEC。次年，MIT的Ketterle在23Na中也实现了BEC。三人在2001年获得诺贝尔物理学奖。
| 实验参数 |
87Rb (JILA 1995) |
23Na (MIT 1995) |
| 初始原子数 |
106 |
2×107 |
| 冷却技术 |
激光冷却 + 蒸发冷却 |
激光冷却 + 蒸发冷却 |
| 最终原子数 |
约 2000 |
约 106 |
| 临界温度 |
约 170 nK |
约 2 μK |
| 凝聚寿命 |
约 15秒 |
约 20秒 |
| 囚禁阱 |
磁阱 (Ioffe-Pritchard) |
磁阱 (Cloverleaf) |
| 观测方法 |
吸收成像 |
吸收成像 |
实验信号特征：在温度从 T>Tc 降低到 T<Tc 的过程中，动量分布发生急剧变化：
- T>Tc：热气体的动量分布是各向同性的高斯型
- T≈Tc：中心出现双峰结构——热峰+凝聚峰
- T≪Tc：几乎全是凝聚峰，呈各向异性（反映阱的几何形状）
干涉实验：1997年，Ketterle小组将两个独立的BEC合并，观察到了清晰的干涉条纹。条纹间距约为15 μm，对比度超过80%，直接证明了BEC的宏观相位相干性。
| 系统 |
粒子 |
Tc (近似) |
特点 |
| 稀释原子气体 |
87Rb, 23Na, 7Li |
100 nK-2 μK |
可调相互作用（Feshbach共振），高度可控 |
| 液氦-4 |
4He |
2.17 K |
强相互作用BEC，只约10%凝聚 |
| 激子BEC |
电子-空穴对 |
1-10 K |
半导体中实现 |
| 光BEC |
光子 |
室温 |
在染料腔中实现，2010年 |
| 磁振子BEC |
自旋波量子 |
10-100 K |
固态系统中的BEC |
| 库珀对BEC |
电子对（超导） |
1-100 K |
弱耦合BCS极限 |
值得特别指出：BCS超导态和BEC之间存在连续演化的联系。从弱耦合BCS到强耦合BEC的"BCS-BEC跨越"是当前量子多体物理的重要研究方向，已在超冷原子Feshbach共振实验中得到验证。
超流（Superfluidity）是指液体在低温下表现出零黏度的奇特流动现象。1937年，Kapitza和Allen几乎同时独立发现了液4He的超流特性（两人后来都因此获得诺贝尔奖）。
当4He冷却到2.17 K以下时，一系列惊人的现象出现：
- 零黏度流动：通过比头发丝还细的毛细管（直径约0.5 μm）时毫无阻力。如果将一个装满超流4He的杯子放在开口处，液体会在数秒内全部流空——即使玻璃杯没有任何缝隙！
- 喷泉效应：用光照加热超流4He的一个局部区域，液体会从细管中向上喷出，形成高度可达数厘米的喷泉。
- 爬壁现象：超流4He沿容器壁形成约30 nm厚的薄膜（Rollin膜），以约20 cm/s的速度"爬"过器壁，自动流出容器。
- 热导率巨高：超流4He的热导率是室温铜的数千倍，1 cm3的超流4He可通过≈ 1 kW的热流而温差极小。
为了理解这些看似矛盾的现象，Tisza（1938）和Landau（1941）提出了二流体模型：超流氦由两种互穿但互不作用的成分组成。
| 成分 |
密度 |
熵 |
黏度 |
流动特性 |
| 超流分量 |
ρs |
零 |
零 |
无旋、无摩擦流动 |
| 正常分量 |
ρn |
有 |
有 |
常规黏滞流动 |
总密度 ρ=ρs+ρn。两种成分的比例随温度变化：
| 温度 (K) |
ρs/ρ |
ρn/ρ |
物理状态 |
| 0 |
1.0 |
0 |
全部超流 |
| 1.0 |
0.82 |
0.18 |
以超流为主 |
| 1.5 |
0.55 |
0.45 |
接近各半 |
| 2.0 |
0.25 |
0.75 |
以正常为主 |
| 2.17 (Tλ) |
0 |
1.0 |
正常流体 |
热力学关系：超流速度 vs 与相位梯度相关：
vs=mℏ∇θ
正常速度 vn 则服从经典的Navier-Stokes方程。
**第二声波：二流体模型中存在两种声波模式：
- 第一声波：通常的密度波（压力振荡），正常和超流分量同向运动
- 第二声波：温度波（熵振荡），正常和超流分量反向运动，总密度不变
第二声波的速度在 T/Tλ≈0.7 时达到约 20 m/s，是超流氦独有的标志性特征。
超流最深刻的理论预言来自Onsager（1949）和Feynman（1955）：超流的环流是量子化的。
由于序参量 Ψ=ρseiθ 必须是单值函数，绕闭合回路一周相位变化只能是 2π 的整数倍：
∮vs⋅dl=mh×n
其中 n=0,±1,±2,… 是涡旋量子数，h/m≈9.97×10−4 cm2/s 是4He的环流量子。
单涡旋的结构：围绕涡旋芯（正常态，直径约1 Å），超流速度随距离变化：
vs(r)=mrnℏ
涡旋芯的能量为：
Ev=πρs(mnℏ)2ln(ξR)
其中 R 是系统尺寸，ξ 是相干长度（约1 Å）。对于半径为 1 mm 的容器：
En=1≈4×104kB（约 5.5×10−18 J）
当BEC被旋转时，会形成规则的涡旋阵列。旋转频率 Ω 与涡旋面密度的关系为：
nv=πℏmΩ
| 旋转频率 Ω (2 π rad/s) |
理论涡旋密度 (cm−2) |
典型实验中的涡旋数 |
| 10 |
1.4×107 |
0-3 (单个涡旋出现) |
| 30 |
4.1×107 |
8-15 (开始形成晶格) |
| 50 |
6.9×107 |
25-45 (三角形晶格) |
| 80 |
1.1×108 |
80-120 (完美的三角形晶格) |
涡旋晶格的稳定构型是三角形晶格（与第二类超导体中的Abrikosov磁通涡旋晶格完全类似）。这是通过最小化涡旋之间的相互作用能得到的：每个涡旋在超流中产生速度场，其他涡旋在这个速度场中感受到Lorentz力，最终排列在最低能量构型上。
超导（Superconductivity）于1911年由H.K. Onnes在汞（Hg）中发现：当汞冷却到4.2 K以下时，电阻突然完全消失，测量精度达到 10−15 Ω·cm。Onnes因此获得1913年诺贝尔物理学奖。
超导的四个基本特征：
- 零电阻：直流电阻精确为零。如果在一个超导环中通过电流，电流可以持续数年不衰减（实验观测到超过2年）。
- 迈斯纳效应（Meissner effect, 1933）：超导体完全排斥内部磁通线。这是热力学平衡态行为，不仅仅是零电阻的结果。
- 临界温度 Tc：高于此温度，超导被破坏，材料回到正常态。
- 临界磁场 Hc（或 Hc1,Hc2）：外磁场超过临界值后超导消失。
1957年，Bardeen、Cooper和Schrieffer发表了里程碑式的BCS理论，首次从微观角度解释了常规超导。这是凝聚态物理史上最成功的理论之一，三人共同获得1972年诺贝尔物理学奖。
BCS理论的物理图像：电子通过晶格振动（声子）媒介产生有效吸引。
具体机制：一个电子穿过晶格时，吸引正离子，在电子经过的路径上造成局部的正电荷过剩（晶格畸变）。这个正电荷过剩区域吸引了另一个电子。在净效果上，两个电子之间产生了一个有效吸引力。这个吸引力虽然很弱（约0.01 eV量级），但在费米面附近可以压倒屏蔽后的库仑排斥，使得两个电子形成束缚态——库珀对。
Cooper定理：无论在费米面上方的吸引多弱，总存在一个两电子束缚态（束缚能 Δ）。
BCS的能隙方程是自洽方程：
Δk=−k′∑Vkk′2Ek′Δk′tanh(2kBTEk′)
其中 Ek=(ϵk−μ)2+Δk2 是准粒子激发谱，Δk 是超导能隙。
简化形式：考虑简单接触吸引势（在费米面 ∣ϵk−μ∣<ℏωD 范围内为常数 V，其他为0），能隙与温度无关：
Δ=ℏωDsinh(N(0)V1)1
在弱耦合极限 N(0)V≪1 下：
Δ≈2ℏωDe−1/(N(0)V)
临界温度公式：
kBTc≈1.14ℏωDe−1/(N(0)V)
以铝（Al）为例：
| 参数 |
符号 |
数值 |
| 态密度 |
N(0) |
1.2×1022 eV−1cm−3 |
| 有效吸引势 |
V |
约 0.15 eV·Ω |
| 耦合常数 |
N(0)V |
0.175 |
| Debye温度 |
ΘD |
428 K |
| Debye频率 |
ωD |
5.6×1013 rad/s |
代入公式：
Tc≈1.45ΘDe−1/N(0)V≈1.45428e−5.71≈295×0.0033≈0.97 K
实验值 Tc=1.2 K，理论与实验吻合。同样方法对其他元素：
| 金属 |
N(0)V (BCS拟合) |
ΘD (K) |
Tc 理论 (K) |
Tc 实验 (K) |
| Al |
0.175 |
428 |
0.97 |
1.20 |
| Sn |
0.245 |
195 |
4.8 |
3.72 |
| Pb |
0.395 |
105 |
7.9 |
7.19 |
| Nb |
0.325 |
277 |
12.3 |
9.25 |
BCS理论的一个直接推论：临界温度与原子的同位素质量有关：
Tc∝M−α,α≈0.5
因为 Debye频率 ωD∝M−1/2，而 Tc∝ωD。
实验验证（Hg同位素）：
| 同位素 |
质量 M |
Tc (K) |
计算 α |
| 198Hg |
198 |
4.185 |
— |
| 199Hg |
199 |
4.173 |
0.49 |
| 200Hg |
200 |
4.161 |
0.50 |
| 202Hg |
202 |
4.139 |
0.50 |
| 204Hg |
204 |
4.119 |
0.51 |
α 的观测值与BCS理论的0.50非常接近，有力地证明了超导的声子机制。
根据对磁场的响应，超导体分为两类：
| 类型 |
磁响应 |
临界场 |
典型材料 |
相变特征 |
| 第一类 |
完全迈斯纳（完全排磁），直到 Hc 时突变到正常态 |
一个临界场 Hc |
Pb, Hg, Sn, Al |
一级相变 |
| 第二类 |
H<Hc1 完全排磁，Hc1<H<Hc2 混合态 |
两个临界场 Hc1,Hc2 |
Nb, NbTi, Nb3Sn, YBCO |
二级相变 |
混合态中的磁通涡旋：第二类超导体在 Hc1<H<Hc2 区间，磁通线以量子化磁通管形式穿透超导体。每个磁通量子：
Φ0=2eh≈2.07×10−15 Wb
类似超流中的涡旋，磁通涡旋也排列成三角形晶格（Abrikosov晶格）。
| 类别 |
配对机制 |
配对对称性 |
Tc 范围 |
例子 |
| s波（常规） |
电子-声子 |
s波（各向同性） |
< 30 K |
几乎所有元素超导体，NbTi |
| d波（高温） |
磁机制（争议） |
d波（有节点） |
< 165 K |
铜氧化物（YBCO, BSCCO） |
| p波 |
铁磁涨落 |
p波（自旋三重态） |
< 3 K |
3He, Sr2RuO4 |
| s±波 |
自旋涨落 |
符号变化的s波 |
< 56 K |
铁基超导体 |
1986年，Bednorz和Müller在IBM苏黎世发现La-Ba-Cu-O陶瓷在约35 K出现超导，开创了高温超导时代。1987年，吴茂昆和朱经武发现YBa2Cu3O7（Tc=92 K），首次突破液氮沸点（77 K）。
高温超导体的关键性质：
- 层状结构：导电发生在CuO2平面内
- d波配对对称性：能隙在 kx=±ky 方向有节点
- 存在赝能隙相（pseudogap phase）：在 Tc 以上就出现能隙
- 相图复杂：从反铁磁Mott绝缘体到d波超导的完整相图
| 年代 |
材料 |
Tc (K) |
突破 |
| 1911 |
Hg |
4.2 |
首次发现超导 |
| 1954 |
Nb3Sn |
18 |
A15结构首次实用 |
| 1973 |
Nb3Ge |
23.2 |
常规BCS超导体的记录 |
| 1986 |
La2−xBaxCuO4 |
35 |
首次发现高温超导 |
| 1987 |
YBa2Cu3O7 |
92 |
超过液氮温度（77 K） |
| 1988 |
Bi2Sr2Ca2Cu3O10 |
110 |
BSCCO多铜氧面层 |
| 1993 |
HgBa2Ca2Cu3O8 |
134 |
常压最高记录 |
| 2008 |
SmFeAsO1−xFx |
55 |
铁基超导体 |
| 2015 |
H3S (高压155 GPa) |
203 |
高压氢化物超导 |
| 2019 |
LaH10 (高压170 GPa) |
250 |
近室温超导 |
高温超导机理的待解之谜：虽然BCS理论完美解释了常规超导，但高温超导的配对机制仍有争议。主要观点包括自旋涨落介导配对、电荷涨落、激子机制等。这是过去近40年凝聚态物理中最重要的未解问题之一。
面对多体问题的指数复杂性，物理学家开发了多种数值计算方法：
| 方法 |
英文缩写 |
基本思想 |
适用场景 |
计算复杂度 |
| Hartree-Fock |
HF |
单Slater行列式近似 |
弱关联，绝缘体 |
O(N3) |
| 密度泛函理论 |
DFT |
电子密度作为基本变量 |
基态性质，固体能带 |
O(N3) |
| 精确对角化 |
ED |
完全对角化有限系统 |
小系统参考解 |
O(eN) |
| 配置相互作用 |
CI |
多组态展开 |
分子精确计算 |
O(N!) |
| 耦合簇方法 |
CC |
指数形式关联 |
量子化学 |
O(N6) |
| 动力学平均场 |
DMFT |
局域自能近似 |
强关联，Mott绝缘体 |
O(N4) |
| 量子蒙特卡洛 |
QMC |
随机采样 |
精确统计 |
取决于符号问题 |
| 密度矩阵重正化群 |
DMRG |
低纠缠态近似 |
一维系统 |
O(N) |
| 张量网络方法 |
TN |
量子态的张量表示 |
有限纠缠系统 |
O(Nχ3) |
用不同方法计算H2分子的结合能，以完全CI（Full CI/FCI）结果为基准：
| 方法 |
结合能 (eV) |
与FCI的绝对误差 |
与实验的误差 |
计算时间 |
| Hartree-Fock (HF) |
3.64 |
1.11 eV |
1.11 eV |
< 0.1秒 |
| DFT-LDA |
4.91 |
-0.16 eV |
-0.16 eV |
0.1秒 |
| DFT-B3LYP |
4.86 |
-0.11 eV |
-0.11 eV |
0.2秒 |
| MP2 |
4.64 |
0.11 eV |
0.11 eV |
0.5秒 |
| CCSD |
4.73 |
0.02 eV |
0.02 eV |
2秒 |
| CCSD(T) |
4.76 |
-0.01 eV |
-0.01 eV |
5秒 |
| 完全CI (FCI) |
4.75 |
0.00 eV |
0.00 eV |
60秒 |
| DMRG (χ=1024) |
4.75 |
0.00 eV |
0.00 eV |
5秒 |
| 实验值 |
4.75 |
— |
— |
— |
量子蒙特卡洛（Quantum Monte Carlo）家族包含多种方法。以变分蒙特卡洛（VMC）为例：
步骤分解：
-
构造试探波函数 ΨT(R;α)，其中 R=(r1,r2,…,rN) 是构型，α 是变分参数。常用形式是Slater-Jastrow波函数：
ΨT=D↑D↓exp[i<j∑u(rij)]
其中 Dσ 是自旋 σ 的Slater行列式，指数部分描述短程关联。
-
计算局域能量：
EL(R)=ΨT(R)H^ΨT(R)
-
Metropolis采样：按照概率分布 ∣ΨT(R)∣2/∫∣ΨT(R′)∣2dR′ 随机行走，每一步提议新的构型 R′→R′，以概率接受：
Paccept=min(1,∣ΨT(R)∣2∣ΨT(R′)∣2)
-
计算能量期望值：
E(α)=⟨EL⟩∣ΨT∣2=M1k=1∑MEL(Rk)
-
优化变分参数 α 使 E(α) 最小化。
示例：对均匀电子气（rs=5，密度参数）：
| 方法 |
关联能 (eV/电子) |
统计误差 |
| RPA |
-0.105 |
— |
| VMC (Slater-Jastrow) |
-0.182 |
$\pm$0.005 |
| DMC (扩散蒙特卡洛) |
-0.193 |
$\pm$0.003 |
| 精确结果 |
-0.194 |
$\pm$0.001 |
费米子符号问题（Fermion Sign Problem）：对费米子系统，扩散蒙特卡洛（DMC）在投影过程中会遭遇符号问题——波函数节点的近似导致指数增长的方差。这是QMC方法面临的核心理论挑战。
超冷原子在光晶格中模拟凝聚态物理的哈密顿量，是当前多体物理的前沿方向。用激光驻波形成周期势阱，原子在其中运动，通过调节激光强度和磁场，可以连续控制模型参数。
Bose-Hubbard模型是光晶格玻色子的标准模型：
H^=−t⟨i,j⟩∑b^i†b^j+2Ui∑n^i(n^i−1)−μi∑n^i
U/t 比决定了相：
| U/t 范围 |
相 |
ni 的涨落 |
相位相干性 |
| U/t<2.2 (1D) |
超流相 (SF) |
大（Δn∼⟨n⟩） |
全局相干 |
| U/t>3.3 (1D) |
Mott绝缘体 (MI) |
小（Δn≈0） |
无序 |
| 中间区域 |
临界 |
发散（量子临界） |
— |
2002年，Greiner等在慕尼黑首次通过改变激光强度，观察到 87Rb从超流到Mott绝缘体的量子相变。这是量子模拟的标志性实验。
| 应用领域 |
相关多体系统 |
应用详情 |
当前状态 |
| 量子计算 |
超导量子比特（Josephson结阵列） |
IBM、Google已实现100+量子比特处理器 |
商业化 |
| 量子传感 |
BEC干涉仪 |
可测重力梯度 10−9 g/Hz |
实验室原型 |
| 电力传输 |
高温超导电缆 |
载流密度比铜高100倍以上，零损耗 |
商业示范 |
| 医疗诊断 |
超导MRI磁体 |
3T、7T全身MRI使用NbTi超导线 |
大规模商用 |
| 粒子物理 |
超导加速器腔 |
LHC使用Nb腔实现16 MV/m加速梯度 |
运行中 |
| 磁悬浮 |
高温超导块材 |
实用化钉扎磁悬浮，载重可达100 kg/m2 |
商业示范 |
| 凝聚态研究 |
量子自旋液体 |
候选材料（κ-(BEDT-TTF)2Cu2(CN)3） |
活跃研究 |
量子多体物理仍在快速发展，几个前沿方向：
- 高温超导机制：铜氧化物和铁基超导的配对机制仍是最大的理论挑战
- 量子自旋液体：非传统磁性序的寻找和验证
- 拓扑量子物态：Majorana费米子、拓扑绝缘体/超导体
- 非平衡量子多体动力学：淬火、热化、多体局域化
- 量子模拟走向实用：超越经典计算机的量子模拟器
- 室温超导的探索：高压氢化物提供了新路径
多体量子系统理论构成了现代凝聚态物理的基石，其关键概念总结如下：
- 全同粒子与量子统计是微观物质行为的基本框架，玻色子和费米子由交换对称性区分，决定了截然不同的宏观表现
- 二次量子化提供了描述粒子数可变系统的强大数学语言，场算符框架是多体理论的"母语"
- 平均场理论是处理弱关联系统的主要工具，但在强关联区域存在根本局限
- BEC、超流和超导是多体量子系统中最具代表性的涌现现象，展现了宏观量子相干性
- 从BCS到高温超导展示了从已理解到待解问题的广阔研究空间
- 数值计算方法（HF、DFT、DMFT、QMC、DMRG）的发展使研究复杂真实材料成为可能
- 超冷原子量子模拟开启了多体物理研究的新范式
多体量子物理仍在快速演变，高温超导机理、量子自旋液体、拓扑量子物态等前沿挑战着当前的理论框架，同时驱动着量子计算、量子传感等革命性技术的发展。
- P. Coleman, Introduction to Many-Body Physics, Cambridge University Press (2015)
- A. Altland & B. Simons, Condensed Matter Field Theory, 2nd ed., Cambridge University Press (2010)
- J. W. Negele & H. Orland, Quantum Many-Particle Systems, Addison-Wesley (1988)
- G. D. Mahan, Many-Particle Physics, 3rd ed., Springer (2000)
- A. J. Leggett, Quantum Liquids: Bose Condensation and Cooper Pairing in Condensed-Matter Systems, Oxford University Press (2006)
- M. H. Anderson et al., "Observation of Bose-Einstein condensation in a dilute atomic vapor", Science 269, 198 (1995)
- M. Greiner et al., "Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms", Nature 415, 39 (2002)
- J. Bardeen, L. N. Cooper & J. R. Schrieffer, "Theory of Superconductivity", Phys. Rev. 108, 1175 (1957)
- I. Bloch, J. Dalibard & W. Zwerger, "Many-body physics with ultracold gases", Rev. Mod. Phys. 80, 885 (2008)
- P. W. Anderson, "More is different", Science 177, 393 (1972)