Majorana 费米子与拓扑超导 · 学习笔记
基于 Kitaev 2001 年论文 “Unpaired Majorana fermions in quantum wires” 的系列问答整理
目录
- U(1) → Z₂ 对称性破缺
- 为何需要超导系统
- Kitaev 链:拓扑相与平庸相
- 体能隙与准粒子局域化
- BCS 与 Ginzburg-Landau 中的 U(1) 破缺
- 一维 SOC 纳米线
- 拓扑超导中的量子涡旋
- Majorana 零能模与基态简并
- Majorana 编织与非阿贝尔统计
- 拓扑量子比特:最小配置
- 动量守恒的困境与破局
- U(1) → Z₂ 对称性破缺
核心机制
| 操作 | 数学形式 | 物理含义 |
|---|
| U(1) 变换 | aj↦eiϕaj | 电荷(粒子数)守恒 |
| 超导配对项 | ajaj+1→e2iϕajaj+1 | 破坏任意 ϕ 的 U(1) 对称 |
| Z₂ 残留 | aj↦−aj(ϕ=π) | 费米子宇称守恒 |
Kitaev 链哈密顿量: H=∑j[−t(aj†aj+1+h.c.)−μaj†aj+Δ(ajaj+1+h.c.)]
结论: 超导序参量将电荷守恒定律降级为费米子宇称守恒定律。守恒量变为: P=(−1)Ntot,Ntot=∑jaj†aj
- 为何需要超导系统
关键论证:U(1) 禁止孤立 Majorana
复费米子可分解:a=21(γ1+iγ2)
在 U(1) 对称系统中,a→eiϕa 导致 Majorana 算符混合: (γ1′ γ2′)=(cosϕsinϕ−sinϕcosϕ)(γ1γ2)
Z₂ 对称性(超导)下: 变换矩阵为 −I(对角),γ1 和 γ2 独立变换,不相互混合。
逻辑链条:
- 拓扑量子计算需要非局域自由度(unpaired Majorana)
- 要局域化,对称性变换不能混合 Majorana 对
- U(1) 会混合 → 必须破缺 U(1)
- 超导将 U(1) 降级为 Z₂ → Majorana 可独立存在
- ∴ 必须在超导系统中寻找 Majorana
- Kitaev 链:拓扑相与平庸相
两个极端情形对比
| 平庸相 (a) | 拓扑相 (b) |
|---|
| 参数 | Δ=t=0, μ<0 | Δ=t>0, μ=0 |
| 哈密顿量 | Ha=−i2μ∑jc2j−1c2j | Hb=it∑jc2jc2j+1 |
| 配对方式 | 格点内(Intra-site) | 格点间(Inter-site) |
| 边界态 | ✗ 无 | ✓ 未配对 Majorana |
为何存在边界态
物理格点不可分割: 一个电子格点由 c2j−1,c2j 紧紧绑定。 在拓扑相中,配对跨越了原子边界,链端点处的 c1 和 c2L 必然”落单”。
有限尺寸效应
端点 Majorana 的波函数以相干长度指数衰减: ϕL(x)∼e−x/ξ,ϕR(x)∼e−(L−x)/ξ
有效耦合(隧穿振幅): Heff≈iδE,b′b′′,δE∝e−L/ξ
∴ 链越长,拓扑保护越强。
- 体能隙与准粒子局域化
等价关系
体能隙⟺零能波函数指数衰减⟺隧穿振幅指数抑制
BdG 谱: E(k)=(2tcosk+μ)2+4Δ2sin2k
衰减长度(相干长度): l0=ΔgapvF
- 能隙越大 → l0 越短 → 端点 Majorana 越难重叠 → 保护越强
- 能隙关闭 Δ→0 → l0→∞ → 拓扑保护失效
- BCS 与 Ginzburg-Landau 中的 U(1) 破缺
BCS:平均场”强行”破缺 U(1)
HMF∼∑k(Δ∗c−k↓ck↑+Δck↑†c−k↓†)
- 包含 c†c†(产生2个电子)和 cc(湮灭2个电子)
- [HMF,N^]=0 → 粒子数不守恒
- BCS 基态 ∣ΨBCS⟩∼∏k(uk+vkck†c−k†)∣0⟩ 是不同粒子数态的线性叠加
GL 方程:保留 U(1) 对称性
F=α∣Ψ∣2+β∣Ψ∣4+2m∗1∣(−iℏ∇−2eA)Ψ∣2
物理定律(自由能泛函)是对称的,对称性破缺发生在系统选择基态时。
序参量如何”选择”相位
- 不稳定性: T<Tc 后,原点变为局域极大值(“帽顶”)
- 随机涨落: 热/量子涨落给系统一个微小随机推力
- 正反馈: 梯度力迅速将系统推向某个随机相位 θ
- 相位相干: 梯度能 Fgrad∝∣∇θ∣2 将相位锁定至全局一致
可观测量: 单个超导体的绝对相位不可观测;约瑟夫森结中的相位差 φ=θA−θB 是物理量,决定超流电流 I=Icsinφ。
- 一维 SOC 纳米线
哈密顿量
H=2m∗px2+αRpxσy+Vzσz+超导项
| 项 | 作用 |
|---|
| Vzσz(塞曼场,沿 z 轴) | 打破时间反演对称性,产生涡旋 |
| αRpxσy(Rashba SOC) | 自旋-动量锁定,有效场在 xy 平面 |
关键几何关系
Beff(k)∝k×z^
- 沿 x 方向运动 → 感受到沿 y 方向的有效磁场(σy)
- 沿 y 方向运动 → 感受到沿 −x 方向的有效磁场(−σx)
- B∗Z⊥B∗eff → 竞争作用打开拓扑能隙
无磁场时:偶数费米面(时间反演对称)
色散关系配方: E±(k)=2m∗ℏ2(k±kSO)2−ESO
Kramers 定理保证费米点数量总是偶数 → 无法实现拓扑相。
- 拓扑超导中的量子涡旋
3个涡旋的情形:3+1模式
在有限二维拓扑超导岛(p+ip 超导)中引入 3 个涡旋:
- 3个 Majorana 局域在涡旋核心(γ1,γ2,γ3)
- 1个 Majorana 沿外边界离域分布(γedge)
原因: 希尔伯特空间要求 Majorana 总数为偶数;体-边对应要求外边界承载手性边缘模。
量子信息编码
| 费米子 | 构成 | 编码 |
|---|
| f12=(γ1+iγ2)/2 | 涡旋1+2 | 逻辑量子比特 |
| f3e=(γ3+iγedge)/2 | 涡旋3+边界 | 辅助自由度 |
- Majorana 零能模与基态简并
为什么 b′,b′′ 不出现在哈密顿量中 = 零能
在 sweet spot(t=Δ,μ=0): H=2t∑j=1L−1(dj†dj−21)+Zero Mode0⋅a~†a~
- ∣ψ0⟩(空态)和 ∣ψ1⟩(占据态)能量完全相同
- 费米子宇称守恒 [H,P]=0 → 两态被对称性”隔离”,不会自发跃迁
宇称本征值
−ib′b′′∣ψ0⟩=+∣ψ0⟩(偶宇称) −ib′b′′∣ψ1⟩=−∣ψ1⟩(奇宇称)
- Majorana 编织与非阿贝尔统计
Ivanov 规则(涡旋交换变换)
γ1→γ2,γ2→−γ1
负号来源: γ2 在运动路径中穿过涡旋1的拓扑”割线”(Branch Cut),获得相位 eiπ=−1。
宇称守恒验证:
Pafter=−i(γ2)(−γ1)=iγ2γ1=−iγ1γ2=Pbefore✓
编织酉算符
U12=exp(4πγ1γ2)=21(1+γ1γ2)
推导(Clifford 代数 + Taylor 展开):
由于 (γnγm)2=−1(类比虚数单位 i),类比欧拉公式: exp(θγnγm)=cosθ+γnγmsinθ
代入 θ=π/4 即得上式。
作用在量子态上
U12∣0⟩=eiπ/4∣0⟩,U12∣1⟩=e−iπ/4∣1⟩
对叠加态 α∣0⟩+β∣1⟩:
∣ψfinal⟩=eiπ/4(α∣0⟩−iβ∣1⟩)
相对相位差为 π/2 → 等价于绕 Z 轴旋转 90°(S 门)。
矩阵表示
U12(对角,同组交换): U12=diag(e−iπ/4, eiπ/4, e−iπ/4, eiπ/4)
U23(跨组交换,产生量子纠缠): U23=211i00i100001i00i1
非相邻交换的合成(“华容道”)
U13≡exp(4πγ1γ3)=U12U23†U12†
物理意义: 即使芯片只支持最近邻操作,也可通过编译等效任意长程量子门:
- U12:把 A 从位置1搬到位置2
- U23†:把 A 从位置2搬到位置3,C 被推到中间
- U12†:把 C 从位置2归还到位置1
- 拓扑量子比特:最小配置
为什么需要 4 个 Majorana(2 个准费米子)
核心限制:费米子宇称超选律
| 系统 | 希尔伯特空间 | 问题 |
|---|
| 2个 Majorana(1个准费米子) | ${ | 0\rangle, |
| 4个 Majorana(2个准费米子) | ${ | 00\rangle, |
逻辑量子比特定义
∣0⟩L≡∣00⟩(偶宇称) ∣1⟩L≡∣11⟩(偶宇称)
拓扑保护: 任何局域噪声(只作用于单端 Majorana)只能将系统踢入奇宇称子空间,被能隙保护阻止,无法区分 ∣0⟩L 和 ∣1⟩L。
- 动量守恒的困境与破局
困境:动量失配
在一维量子线中加塞曼场(无 SOC):
kF↑=ℏ2m∗(μ+VZ)=kF↓=ℏ2m∗(μ−VZ)
s-波超导体只能提供 P=0 的库珀对,而纳米线配对需要 Δk=kF↑−kF↓=0,违背动量守恒。
破局:Rashba SOC
SOC 使色散关系发生水平偏移,形成自旋-动量锁定。在强 SOC + 塞曼场下:
- 费米点仍在 ±kF,但自旋是上/下的混合态
- s-波超导体可以找到动量互反、含相反自旋分量的准粒子 → 完美实现 P=0 配对
一句话:SOC 是修复磁场造成动量失配的”补丁”,是实验上实现拓扑超导的关键。
为何不用连续一维模型(Kitaev 的选择)
Nielsen-Ninomiya 定理(费米子加倍定理): 在任意严格的一维晶格模型中,右行模式与左行模式数量必须严格相等。
- 1D 系统边界是 0D 的点,无法承载”手性”(定向传播)
- 要实现手性边缘态,必须依赖高维体态作为支撑(2D 量子霍尔效应)
∴ Kitaev 选择格点紧束缚模型,利用晶格 Umklapp 过程自然处理动量失配。
关键公式速查
| 物理量 | 公式 |
|---|
| 相干长度 | ξ=ℏvF/Δgap |
| 有限尺寸能量劈裂 | δE∝e−L/ξ |
| 编织算符 | U12=21(1+γ1γ2) |
| Majorana 积分关系 | γ1γ2=i(1−2n^1) |
| 费米子宇称算符 | P^=1−2N^ |
| BdG 谱 | E(k)=(2tcosk+μ)2+4Δ2sin2k |
整理自 Gemini 3 Pro 对 Kitaev 2001 论文的系列解读