本页目录
固体 IV · Berry 相位与拓扑能带
对标:Berry 相位、Chern 绝缘体与二维两带模型 | 前置:晶格与倒空间、Bloch 能带、两能级系统;几何符号见Berry 计算桥 “有边缘态”是一个醒目的实验信号,却不是拓扑的定义。本讲用一个能算到底的两带模型,把 bulk gap、Berry phase、Berry curvature、Chern invariant 和边界谱放进同一套账本;你会看到边缘态怎样支持 bulk-boundary correspondence,也会看到为什么一条漂亮的边界色散不能独自证明 Chern 数。
学习层:边缘态为什么不能替代拓扑不变量?
1. 具体谜题:同一条边界色散,究竟证明了什么?
二维晶体的 Bloch 波函数不仅有能量,还有随着 \(\mathbf k\) 改变的相位结构。相位可以局部改写,却可能在整个 Brillouin zone 上留下不能连续抹掉的整数。实验上,边缘输运或扫描谱常先看到一条穿过 bulk gap 的支,于是很容易把“有 edge state”直接等同于“有拓扑”。
先把问题说尖锐:
- 改变质量参数 \(m\) 时,bulk gap 在哪里必须闭合,系统才能从 \(C_-= +1\) 变成 \(C_-=-1\)?
- 固定 \(k_y=0\) 绕一圈 \(k_x\),Berry phase 是整张二维 BZ 的 Chern 数吗?
- 若某个 \(k_y\) 切片没有边缘态,它能否推翻整个二维系统的拓扑?
本实验采用 QWZ 两带模型。你先预测,再打开曲率热图、gap、Berry phase、Chern 数和 半无限边界谱;所有数值都来自同一个 Hamiltonian。
2. 先预测:让相变先发生在你的脑中
实验默认 \(m=-1,\ k_y=0\)。这里的 \(k_x,k_y\) 都以晶格常数为 1 的无量纲倒空间坐标表示。
- 把 \(m\) 从 \(-1\) 连续调到 \(+1\),若下带 Chern 数从 \(+1\) 变为 \(-1\),中途的 bulk gap 能否始终大于零?
- 在默认切片 \(k_y=0\) 上,Berry phase 更像 \(0\) 还是 \(\pi\)?注意它是这条一维闭合回路的量。
- 默认切片满足边缘条件吗?把 \(k_y\) 改变到 \(\pi\),整个二维相的拓扑会不会因此消失?
选择预测后再揭示结果。读图顺序建议固定为:先看 \(\Delta\) 是否闭合,再对照解析 \(C_-\) 与整张曲率的实际积分,最后用边界谱检查哪些 \(k_y\) 切片承载边缘支。
3. 最小机制:两带 Hamiltonian 把几何变成向量
模型写成
这里 \(\boldsymbol{\sigma}=(\sigma_x,\sigma_y,\sigma_z)\) 是轨道或赝自旋空间中的 Pauli 矩阵;它不必等同于真实电子自旋。两条能带为
因此 bulk direct gap 是
gap 什么时候可能为零?必须同时满足 \(\sin k_x=0,\ \sin k_y=0\),所以只需检查四个高对称点:
| 点 | \((k_x,k_y)\) | \(d_z\) |
|---|---|---|
| \(\Gamma\) | \((0,0)\) | \(m+2\) |
| \(X\) | \((\pi,0)\) | \(m\) |
| \(Y\) | \((0,\pi)\) | \(m\) |
| \(M\) | \((\pi,\pi)\) | \(m-2\) |
相变边界是 \(m=-2,0,2\)。本模型甚至能精确求全局最小能隙。令 \(a=\cos k_x,b=\cos k_y\),则
它对 \(a,b\in[-1,1]\) 分别是线性的,最小值必能在正方形的某个顶点取得。因此
实验直接使用这个解析最小值,临界点显示精确的零。一般模型未必有这条捷径;不含闭隙点的粗网格可能把真正的零漏掉。
4. Berry phase:一条闭合路径的相位账本
设下带归一化本征态为 \(|u_-(\mathbf k)\rangle\)。局部 Berry connection 可写为
沿固定 \(k_y\) 的闭合 \(k_x\) 回路,Berry phase 为
数值上更稳健的离散写法是相邻本征态重叠的乘积:
负号来自这里选取的 \(\mathbf A=i\langle u|\nabla u\rangle\) 约定;换约定可能让所有相位与 Chern 数同时变号,但物理预测不变。\(\gamma\) 只回答“一条回路绕了多少相位”,并不是自动等于二维 BZ 上的整数。
每个态乘任意相位时,相邻重叠中的相位因子沿闭合链逐项抵消,所以这个乘积不依赖本征求解器选的局部相位。实验比较 256 与 512 点的结果;两者之差只提供离散化诊断,不是严格误差界。全局闭隙也不代表每条回路都经过简并:\(m=0,k_y=\pi/2\) 时,\(\mathbf d=(\sin k_x,1,\cos k_x)\),整条回路能隙为 \(2\sqrt2\),仍有
但二维下带的 \(C_-\) 此时未定义。这里相位来自单位球上围绕 \(+y\) 轴的小圆所围立体角的一半;不能用“全局没有 C”推成“局部回路没有几何”。
5. Berry curvature 与 Chern invariant:局部弯曲怎样累积成整数?
对于 \(H=\mathbf d\cdot\boldsymbol{\sigma}\),下带的 Berry curvature 在本实验约定下为
二维 Chern invariant 是整个 BZ 的积分:
当 \(\Delta_{\mathrm{bulk}}>0\) 时,\(C_-\) 只能取整数,并且在连续改变参数时不能改变,除非 gap 闭合。对当前方向、下带和曲率符号约定,四个区间是
四个 Dirac 点的手性为 \(\Gamma,M\) 取 \(+1\),\(X,Y\) 取 \(-1\);各点的质量符号给出
这是本模型的解析相图。实验另列实际中点积分,逐次加密网格到最多 \(512\times512\),同时检查相邻积分差与最小质量尺度。热图只有 \(17\times17\) 个显示色块,不能拿它当积分精度。即使两次网格结果接近,也可能一起漏掉很窄的峰,所以一致性检查仍不是数学上的严格误差界。
例如 \(m=-0.001\) 的解析 \(C_-=+1\),而 \(512\times512\) 中点积分约为 \(0.2112\):这是网格漏峰,不是分数 Chern 数。实验保留该数值并提示尚未充分分辨,不用解析标签覆盖它。在 \(m=0\) 等边界,两带不能在整个 BZ 上保持隔离,绝缘带 Chern 数本身未定义。
更隐蔽的反例:整数也可能来自未分辨的网格。 在相同 \(A=i\langle u|du\rangle\) 约定下,Fukui–Hatsugai–Suzuki(FHS)格点法先取归一化重叠
其中 \(\operatorname{Arg}\) 取主值,各方向周期回接,且重叠须非零。本仓独立 NumPy 本征矢复算,在 \(m=+0.001\)、网格点 \(k_\mu=-\pi+2\pi j/N\) 时,\(N=41\) 得到 \(C_N=+1\),\(N=81\) 才得到解析相图的 \(-1\)。这是格点整数与连续带拓扑尚未对应,不是连续模型在加密时改变拓扑。需检查窄峰的分辨与加密稳定性;只检查“输出整数”不能验收。算法来源见 Fukui–Hatsugai–Suzuki 原论文;这里的具体反例来自本模型的复算,数值保存在独立参考数据中。
6. 边界谱:bulk-boundary correspondence 的一条可检验切片
沿 \(x\) 方向取半无限链、保留 \(k_y\) 为好量子数。实验分别画左、右两种半无限边界的解;有限宽条带的两侧重叠修正在这里没有加入。有效质量与可归一化条件为
满足条件时有理想边界支,其能量为
衰减因子为 \(|\lambda|=|m_{\mathrm{eff}}|\)。所以 \(|m_{\mathrm{eff}}|\) 越小,状态越集中在边界;到 \(m_{\mathrm{eff}}=0\) 时,波函数只占据最外侧一个格点,振幅衰减长度为零。这个公式描述的是干净、平移不变、理想边界的教学模型,不是任意真实表面的完整谱。
默认 \(m=-1\)、\(k_y=0\) 时 \(m_{\mathrm{eff}}=0\),故有 \(E_{\mathrm{edge}}=0\) 的边缘切片。可是同一个质量参数在 \(k_y=\pi\) 时 \(m_{\mathrm{eff}}=-2\),该切片没有理想边界支。这不是矛盾,而是“边界支只存在于一段动量区间”的直接体现。
7. 无 JavaScript 时的静态读法:把数值和公式对账
默认值是 \(m=-1,\ k_y=0\)。静态代入得到
因此默认实验应显示 gap 为 2、下带 Chern 数约为 \(+1\)、Berry phase 约为 \(3.142\ \mathrm{rad}\),并报告当前切片存在零能边缘支。解析标签为精确整数;实际积分与相位还含网格离散误差,不能把浮点舍入误差当成总误差。默认回路因对称性可特别精确,一般切片没有这个保证。
| 参数 | bulk gap \(\Delta\) | 下带 \(C_-\) | \(k_y=0\) Berry phase | \(k_y=0\) 边缘支 |
|---|---|---|---|---|
| \(m=-1\) | \(2.0\) | \(\approx+1\) | \(\pi\) | 有,\(E=0\) |
| \(m=1\) | \(2.0\) | \(\approx-1\) | \(0\) | 无;其他 \(k_y\) 可有 |
| \(m=2.6\) | \(1.2\) | \(\approx0\) | \(0\) | 无 |
| \(m=0\) | \(0\) | 不定义 | 若回路穿过简并点则不定义 | 不可用绝缘体语言判断 |
特别留意第二行:\(m=1\) 仍在非平庸区间,但默认的 \(k_y=0\) 切片不满足 \(|1+\cos0|<1\)。所以“我没在一个动量切片看到 edge”不能推出 \(C_-=0\)。
8. 误区与失败边界
- Berry phase 不是 Chern invariant。前者依赖选定的闭合路径,后者积分整个二维 BZ;只有在额外条件下,随 \(k_y\) 的 Berry phase 变化才与 Chern 数联系起来。
- 边缘态不是拓扑的定义。边界势、有限宽度、无序和表面重构会移动或混合局部谱;Chern 数是在 bulk gap 打开的前提下定义的 bulk 量。
- gap 不能只看一个点,也不能只看粗网格。本模型的候选闭合点由 \(\sin k_x=\sin k_y=0\) 给出;数值网格必须覆盖它们,或单独做高对称点检查。
- 临界点不拥有可靠的绝缘体 Chern 数。gap 闭合时下带与上带不能在整个 BZ 上保持分离,有限采样得到的近似整数没有拓扑保护。
- 符号依赖约定。交换 \(k_x,k_y\) 的方向、改用上带或改变 Berry connection 定义,可能让报告的符号改变;先写清楚坐标取向、能带和公式。
- 两带 QWZ 模型不是所有材料的身份证。真实系统还会有自旋简并、多个轨道、相互作用、非理想边界和 disorder;拓扑账本必须在相应有效 Hamiltonian 上重新建立。
9. 迁移任务:把“边界证据”升级成完整判据
设你在样品的一条边看到穿越费米能的色散支,但 bulk 光谱还有一个未确认的小 gap。请写一份最小验证清单:先如何测 bulk gap,如何在整个 BZ 或等价的参数空间上计算 Berry curvature,再如何判断边缘支是否在同一个 gap 中连接价带与导带?最后说明为什么在 gap 尚未确认打开时,不能把边缘输运直接命名为量子反常霍尔拓扑。
展开核对:怎样排除三种看起来像拓扑的误判?
先扫描整个相关动量范围的体谱,并报告实验分辨率、费米能与温度;“没有看见闭隙”不能排除低于分辨率的小隙或漏测的接触点。再给出占据带、坐标取向和 Berry 约定,用更密网格或独立方法检验全 BZ 的积分。最后改变边界条件或条带宽度,看跨隙的净手性连接能否持续存在,并区分局部边缘峰与真正连接体带的支。
用本模型做三个对照:\(m=-1,k_y=0\) 有全局隙、\(C=1\) 和零能边界态;\(m=1,k_y=0\) 没有这个切片的边界态,仍有 \(C=-1\);\(m=0,k_y=\pi/2\) 的切片有隙和 Berry 相位,半无限边界解也存在,但全局 C 未定义。最后一种情况说明,局部几何与边界证据不能替代全局绝缘条件。
要推到量子反常霍尔输运,还需孤立占据带填满、费米能处于全局隙及合适的低温输运条件;模型中的整数不会自动排除体漏电或接触误差。
若联络与曲率的符号难以跟上,先完成两能级 Berry 计算桥:它从明确的下带本征向量推出正号曲率,并用南北规范核对闭合相位。进一步的态距离见量子几何张量。
1. Bloch 态的相位为什么会留下几何
晶体平移对称性把波函数写成 Bloch 形式
其中 \(u_{n\mathbf k}\) 具有晶格周期。对每一个 \(\mathbf k\),\(|u_{n\mathbf k}\rangle\) 只确定到一个相位:\(|u\rangle\to e^{i\chi(\mathbf k)}|u\rangle\)。因此单点的相位不是可观测量,但沿路径累积的几何相位可以通过干涉、极化、轨道磁矩或输运响应体现出来。
这个“可局部改写、整体可能有障碍”的结构,和电磁学中的规范势很像。Berry connection 在规范变换下改变为
(采用本讲的 \(\mathbf A=i\langle u|\nabla_{\mathbf k}u\rangle\) 约定),所以闭合回路的 \(\gamma\) 只在 \(2\pi\) 意义下不变。曲率
则直接消除了局部规范选择。Chern 数把曲率在闭合 BZ 上积分,得到一个不依赖平滑规范的整数。
2. 两带球面映射与 gap 边界
把 \(\hat{\mathbf d}=\mathbf d/|\mathbf d|\) 看作从动量环面到单位球面的映射。只要 \(\mathbf d\) 不为零,\(\hat{\mathbf d}\) 在整个 BZ 上连续;改变 \(m\) 只能连续变形这张图,不能改变它包裹球面的整数次数。若 \(\mathbf d=0\),归一化失效,正是 band gap closure。
本模型的三个分量各承担一个角色:\(\sin k_x\) 与 \(\sin k_y\) 在高对称点为零,决定可能的 Dirac 接触;\(m+\cos k_x+\cos k_y\) 是质量项,决定这些接触点映到球面的北极还是南极。四个点的质量符号随 \(m\) 改变,导致不同 Dirac 点的半整数曲率贡献重新组合为 \(0,\pm1\)。
可以直接算导数
它们的叉积为
把它与 \(\mathbf d\) 做点积,再除以 \(|\mathbf d|^3\),就得到曲率热图中红蓝斑块的数值。靠近小质量的 Dirac 点时分母变小,曲率集中;远离相变时曲率更平缓。曲率“集中得很高”不等于 Chern 数已经很大,必须再乘面积元并在全 BZ 累积。
3. 从 Berry phase 到 Chern 数的连续关系
对每一个固定 \(k_y\),沿 \(k_x\) 定义 \(\gamma(k_y)\)。只要这条回路上的能带不闭合,\(\gamma\) 随 \(k_y\) 连续变化,允许在 \(2\pi\) 意义下展开。由 Stokes 定理,在没有规范奇点的局部区域内
因此当 \(k_y\) 从 \(-\pi\) 走到 \(\pi\),Berry phase 的总绕行数与 Chern 数相关:
这里需要选择连续的 phase branch,并处理 BZ 边界的规范拼接;直接把每个切片都压回 \([-\pi,\pi]\) 后相减,可能把一次真实的 \(2\pi\) 绕行藏掉。这正是实验台同时画曲率、显示单切片 Berry phase、再对照解析 Chern 数与整区间实际积分 的理由。
在默认 \(m=-1\)、\(k_y=0\),一维回路包围的是一个有效质量为 \(m+1=0\) 的平面 Dirac 结构,离散重叠乘积给出 \(\pi\)。这是一条很好的预测线索,却仍需结合二维积分才能回答整个系统的拓扑类别。
4. 边界、局域化与输运
为什么恰好得到正弦边界谱?固定 \(k_y\),记 \(b=m+\cos k_y\),每个格点有两个轨道。链的方程为
取左端 \(j=1\),删去 \(\psi_0\) 项。令 \(\chi_+=(1,i)^T/\sqrt2\),直接乘矩阵可核对 \(\sigma_y\chi_+=\chi_+\)、\(T^\dagger\chi_+=0\)、\(T\chi_+=\sigma_z\chi_+\)。代入指数试探解 \(\psi_j=N\lambda^{j-1}\chi_+\),得到
因此选 \(\lambda=-b\) 就消去第二项,得到左边界 \(E_L=+\sin k_y\);右端用相反的 \(\sigma_y\) 本征态得到 \(E_R=-\sin k_y\)。几何级数可归一化当且仅当 \(|\lambda|<1\),且 \(N=\sqrt{1-|\lambda|^2}\)。这是半无限链的精确解。
振幅按 \(e^{-j/\xi}\) 衰减,故 \(\xi=-1/\ln|\lambda|\)(晶格单位);概率密度的衰减长度是它的一半。\(\lambda=0\) 时态只在端点,\(|\lambda|\to1\) 时长度发散。有限宽条带中的两端态可能杂化,不能把上述两条精确正弦曲线直接说成任意宽度条带的完整本征谱。
右图浅蓝区域是固定 \(k_y\) 的体谱投影隙,它的半宽为
只需把 \(|\mathbf d|^2=\sin^2k_y+1+b^2+2b\cos k_x\) 对 \(k_x\) 最小化便能得到。边界解落在这个切片隙中;它在某些 \(k_y\) 的能量可以超出全局能隙。全局能隙还要对所有 \(k_y\) 再取最小值,所以两种隙不能混叫。
能量 \(E=\pm\sin k_y\) 在 \(k_y=0\) 和 \(\pi\) 过零,但是否真的属于 edge branch 要先检查存在条件。对于 \(m=-1\),\(|{-1+\cos k_y}|<1\) 给出一段围绕 \(k_y=0\) 的边界区间;对于 \(m=1\),存在区间转移到靠近 \(\pi\) 的位置。质量跨过 \(m=0\) 时,边界支的动量位置和传播方向发生重排,同时 bulk 在 \(X,Y\) 点关闭并重新打开。
真实边界可以加入势 \(V_{\mathrm{edge}}\)、改变 hopping 或破坏 \(k_y\) 平移对称性。这样色散形状不必仍是正弦,局域化长度也会改变,但只要 bulk gap 和相关对称条件保持,跨 gap 的净手性不能被任意小扰动消去。相反,一对相反传播的普通边界态可以互相散射并开 gap;所以“看到一条局部峰”与“存在受 Chern 差保护的净手性支”是不同强度的证据。
规定电子电荷为 \(-e\)(\(e>0\))、\(\sigma_{xy}=j_x/E_y\),零温下孤立下带填满且费米能在能隙内时,量子反常霍尔响应的理想化关系是
在同一约定下 \(\sigma_{yx}=+C_-e^2/h\)。符号可由 \(\hbar\dot{\mathbf k}=-e\mathbf E\)、\(\dot{\mathbf r}=\nabla_kE/\hbar-\dot{\mathbf k}\times\boldsymbol\Omega\) 和 \(\mathbf j=-e\int\dot{\mathbf r}\,d^2k/(2\pi)^2\) 核对。材料实验还要处理费米能位置、体漏电、接触电阻、温度和磁畴。公式给出拓扑能带对输运的贡献,不会替实验自动解决样品是否真的处于绝缘 bulk。
5. 如何报告一次可信的拓扑计算
一份可复现报告至少应包含四层信息。第一,写出 Hamiltonian、坐标方向、下带还是上带以及 Berry curvature 的符号。第二,给出 bulk gap 的最小值与搜索方法,说明是否显式检查高对称点。第三,说明曲率积分的网格、收敛性和 gap 接近零时的数值不稳定。第四,把边界条件、边界势和有限尺寸写清楚,不要只贴一条未标注坐标的谱线。
例如 \(m=-1\) 的解析 \(C_-=+1\),实际积分随网格加密接近 \(+1\);报告应分别列出这两个结果、网格尺寸和相邻差。\(m=2.6\) 的解析标签为 \(0\),实际积分会有很小的残差,残差随算法与网格改变。\(m=0\) 则由精确能隙公式给出 \(\Delta=0\),此时不是“把积分四舍五入为零”,而是全局 C 不定义。接近临界点时,小积分残差或网格一致本身都不足以证明精度。
6. 三个数量级例子
例 1(默认相) 取 \(m=-1\)。\(\Gamma\)、\(X/Y\)、\(M\) 的质量分别为 \(1,-1,-3\),最小的 \(|\mathbf d|\) 为 1,故 \(\Delta=2\);下带积分为 \(C_-=+1\)。在 \(k_y=0\),\(m_{\mathrm{eff}}=0\),边缘能量为 0。
例 2(穿过相变) 从 \(m=-0.4\) 调到 \(m=+0.4\)。在 \(m=0\),\(X\) 与 \(Y\) 的质量同时为零;gap 先降到零,再重新打开。前后两个绝缘区的 Chern 数分别为 \(+1\) 与 \(-1\),变化量为 \(-2\),其绝对值 2 对应两个相关的 Dirac 接触点贡献。若只在 \(m=-0.4\) 和 \(+0.4\) 测量而跳过临界点,很容易错过拓扑改变必须经过 gap closure 这一机制。
例 3(拓扑与切片) 取 \(m=1\)。它在 \(0<m<2\),所以 \(C_-=-1\);但 \(k_y=0\) 时 \(m_{\mathrm{eff}}=2\),没有该切片的理想 edge。改看靠近 \(k_y=\pi\) 的区间,\(m+\cos k_y\) 才可能落入 \((-1,1)\)。这一个例子足以拆掉“一个切片无 edge,所以 bulk 平庸”的错误推理。\(\blacksquare\)
下一页:从拓扑能带的几何不变量继续走向真实材料中的对称性保护、无序与测量误差。
延伸阅读:Qi–Wu–Zhang 原始论文介绍二维拓扑响应模型;Asbóth、Oroszlány 与 Pályi 的课程讲义第 6 章讨论 QWZ 模型及边界态。跨资料比较时先对齐 Berry 连接、上下带与坐标取向的符号约定。