论文身份
论文题为 Edge States and Topological Invariants of Non-Hermitian Systems,作者是 Shunyu Yao 与 Zhong Wang。arXiv 版本 2018 年 3 月 5 日提交,PRL 版本发表于 Physical Review Letters 121, 086803 (2018),并被标注为 Editors' Suggestion。
1. 这篇论文要解决什么难题
在普通厄米拓扑材料里,有一条核心原则叫体边对应:只看无限大体材料的 Bloch Hamiltonian,算出拓扑不变量,就能预言有限样品边界上有没有稳健边缘态。这个原则是拓扑绝缘体、SSH 链、量子霍尔等理论的基石。
但在非厄米系统里,数值上出现了令人困惑的现象:周期边界条件和开边界条件下的谱完全不一样;用普通 Bloch Hamiltonian 算出来的拓扑数,不能正确预言开链的零能边缘态。这看起来像是体边对应彻底失效。
问题不是“体边对应没了”,而是我们以前选错了“体”。普通 Bloch 体态假设 |e^{ik}|=1,但非厄米开链的体态并不是这种波。
姚顺宇和汪忠的答案是肯定的:存在广义体边对应,只是它不是 Bloch 体边对应,而是 non-Bloch bulk-boundary correspondence。
2. 非厄米 SSH 模型:最小战场
论文研究的是一个非厄米版本的 SSH 链。SSH 链本来是最经典的一维拓扑模型:一个晶胞里有 A/B 两个子格,胞内跃迁和胞间跃迁竞争,决定是否有零能边缘态。
非厄米改造后,Bloch Hamiltonian 写成:
其中 γ 表示非厄米强度。直观上,它让向左和向右、A 到 B 和 B 到 A 的跃迁不再像厄米系统那样互为共轭。系统仍有手性对称性,所以能谱成对出现:如果有能量 E,也有 -E。
3. 普通 Bloch 理论错在哪里
如果你只看周期边界条件下的 H(k),能隙关闭点会出现在:
普通 Bloch 拓扑数只能在这些点发生跳变。问题是,论文数值和解析解都显示,开链真正的零模相变点不是这些点,而是:
这就是论文的震撼点:周期边界下看起来很自然的拓扑相变点,不能解释开边界零模。比如论文使用 t2=1, γ=4/3 时,真正相变点约为 1.20,而普通 Bloch 间隙关闭点是 t2 ± γ/2。
4. 非厄米趋肤效应:不是边缘态,是所有体态都贴边
论文解析开链时发现:当 γ ≠ 0,开边界下的体态不再是延展的 Bloch 波,而是整体指数局域到某一端。这就是后来被广泛称为 non-Hermitian skin effect 的现象。
在 t3=0 的简化情况下,可以通过相似变换把非厄米 SSH 变成普通 SSH,但代价是波函数多了一个指数因子。普通 Bloch 因子 e^{ik} 被替换为:
厄米系统里 r=1,所以 β 在单位圆上;非厄米系统里 r 一般不等于 1,波函数会多出指数衰减或指数增长因子。开链能接受的是衰减方向,所以体态堆到边界。
5. 广义布里渊区:Cβ 才是真正的体
论文用更一般的方法重新求解开链。令波函数在每个晶胞上按 β^n 变化,把两个 β 解记为 β1(E) 和 β2(E)。在长链极限下,体能谱必须满足:
这个条件非常关键。它不是额外假设,而是由开边界条件和长链极限推出的。如果两个 β 的模不同,某一项会在长链中指数压倒另一项,无法形成真正的体谱。
在 t3=0 时,条件给出一个半径为 r 的圆;在 t3 非零时,Cβ 不再是圆,而是一条更一般的闭合曲线。这条曲线就是 generalized Brillouin zone。
6. 非布洛赫 winding number
有了 β,就把 Bloch Hamiltonian 改写成 non-Bloch Hamiltonian。做法非常直接:
对 t3=0 模型,论文写成:
非厄米矩阵有左右本征矢,所以论文用 biorthogonal 的方式定义 Q 矩阵。手性对称性让 Q 是非对角形式,里面有一个复函数 q。最终的拓扑不变量是:
这就是 non-Bloch winding number。它和普通 winding number 形式很像,但积分路径不再是单位圆,而是广义布里渊区 Cβ。
| 项目 | 普通 Bloch 拓扑 | 非布洛赫拓扑 |
|---|---|---|
| 波函数因子 | e^{ik} | β = r e^{ik} 或更一般的 Cβ 曲线 |
| 积分路径 | 单位圆 | 广义布里渊区 Cβ |
| 适用边界 | 周期边界体谱 | 开边界体谱 |
| 是否能预言零模 | 非厄米情形一般失败 | 论文中正确预言稳健零模数 |
7. 零模到底怎么被预言
论文得到的结论是:2W 计数左右两端稳健零能模式的总数。比如在 t2=1, γ=4/3 的例子里:
还有一个很细但很重要的点:普通 Bloch gap closing 点并不一定是“零模出现/消失”的点,而可能只是零模从左端迁移到右端的点。论文指出,在某些参数区间,两个零模都在左端;中间区间一个在左端一个在右端;另一侧两个都在右端。总数不变,只是位置迁移。
8. 为什么以前会被“正常化条件”误导
论文专门讨论了一个容易犯的错误。找左端零模时,人们通常要求波函数可归一化,即某个 |β| < 1。把临界条件写成 |β| = 1,会得到与 Bloch gap closing 一致的结果,看起来很合理。
但这正是陷阱:它偷偷假设零模在相变时并入的是普通 Bloch 体态,也就是 |β|=1 的体态。非厄米开链的体态实际在 |β|=r 上。所以真正的并入条件不是 |β|=1,而是:
这一步改正,直接把相变点从普通 Bloch 点移到 ±√(t2²+(γ/2)²)。
9. t3 非零:广义布里渊区不一定是圆
论文主文后半部分加入 t3,使模型更一般。此时 β 方程变成四次方程,存在四个根。体谱由某一对根的模相等条件决定,即 |βi(E)|=|βj(E)|。
这时 Cβ 不再是半径固定的圆,不同能量对应不同 |β|,但沿这条 Cβ 曲线计算 non-Bloch winding number,仍然正确预言零模数。这个部分非常重要,因为它说明论文不是只发现了一个可解特例,而是给出了可推广的程序。
10. 这篇论文的范式意义
这篇论文的贡献可以分成三层。
- 物理现象层:指出非厄米开链中体态会整体局域到边界,也就是非厄米趋肤效应。
- 数学对象层:用 β 和 Cβ 替代普通动量 k 与单位圆,建立广义布里渊区。
- 拓扑原则层:提出 non-Bloch winding number,并建立非布洛赫体边对应。
如果说普通拓扑能带理论的基础是“体态是 Bloch 波”,那么这篇论文告诉我们:非厄米系统里,这个基础要换成“体态是 β 波”。这就是它被视为范式级突破的原因。
一张心智地图
来源与说明
本文基于 arXiv 版本和 PRL 公开页反复阅读后整理,所有示意图为重新绘制的解释图,并非论文原图复刻。公式为方便阅读做了排版简化。