PRL 121, 086803 · Yao & Wang · 2018

非厄米拓扑为什么必须离开普通布里渊区

姚顺宇与汪忠这篇 PRL 的范式突破在于:非厄米系统的体态本身会“贴边”,因此不能再用普通 Bloch 波和普通 winding number 去预言边缘态。真正的拓扑不变量必须定义在广义布里渊区上。

非厄米趋肤效应 广义布里渊区 Cβ non-Bloch winding number non-Bloch bulk-boundary correspondence
Hermitian chain:体态延展 Non-Hermitian chain:所有体态也向一侧堆积 普通 Bloch 因子 e^{ik} 被 β = r e^{ik} 替代,r ≠ 1。

论文身份

论文题为 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。

核心对象非厄米 SSH 链。它保留手性对称性,但跃迁振幅带有非互易/增益损耗效应。
核心发现开边界下,所有体态都可指数局域到边界,形成非厄米趋肤效应。
核心方案把拓扑不变量从普通 Bloch 圆移动到广义布里渊区 Cβ 上定义。

1. 这篇论文要解决什么难题

在普通厄米拓扑材料里,有一条核心原则叫体边对应:只看无限大体材料的 Bloch Hamiltonian,算出拓扑不变量,就能预言有限样品边界上有没有稳健边缘态。这个原则是拓扑绝缘体、SSH 链、量子霍尔等理论的基石。

但在非厄米系统里,数值上出现了令人困惑的现象:周期边界条件和开边界条件下的谱完全不一样;用普通 Bloch Hamiltonian 算出来的拓扑数,不能正确预言开链的零能边缘态。这看起来像是体边对应彻底失效。

问题不是“体边对应没了”,而是我们以前选错了“体”。普通 Bloch 体态假设 |e^{ik}|=1,但非厄米开链的体态并不是这种波。

姚顺宇和汪忠的答案是肯定的:存在广义体边对应,只是它不是 Bloch 体边对应,而是 non-Bloch bulk-boundary correspondence。

2. 非厄米 SSH 模型:最小战场

论文研究的是一个非厄米版本的 SSH 链。SSH 链本来是最经典的一维拓扑模型:一个晶胞里有 A/B 两个子格,胞内跃迁和胞间跃迁竞争,决定是否有零能边缘态。

非厄米改造后,Bloch Hamiltonian 写成:

H(k) = dx σx + (dy + i γ/2) σy

其中 γ 表示非厄米强度。直观上,它让向左和向右、A 到 B 和 B 到 A 的跃迁不再像厄米系统那样互为共轭。系统仍有手性对称性,所以能谱成对出现:如果有能量 E,也有 -E。

Non-Hermitian SSH chain:A/B 子格,胞内跃迁变成不对称 A B unit cell A B A B A B t1 + γ/2 t1 - γ/2 t2

3. 普通 Bloch 理论错在哪里

如果你只看周期边界条件下的 H(k),能隙关闭点会出现在:

t1 = ± t2 ± γ/2

普通 Bloch 拓扑数只能在这些点发生跳变。问题是,论文数值和解析解都显示,开链真正的零模相变点不是这些点,而是:

t1 = ± √( t22 + (γ/2)2 )

这就是论文的震撼点:周期边界下看起来很自然的拓扑相变点,不能解释开边界零模。比如论文使用 t2=1, γ=4/3 时,真正相变点约为 1.20,而普通 Bloch 间隙关闭点是 t2 ± γ/2

通俗比喻:你用“环形跑道”的地图去预测“有墙边界的走廊”里人会聚到哪里。厄米系统里这通常没问题;非厄米系统里,墙本身改变了所有体态的形状,所以地图要重画。

4. 非厄米趋肤效应:不是边缘态,是所有体态都贴边

论文解析开链时发现:当 γ ≠ 0,开边界下的体态不再是延展的 Bloch 波,而是整体指数局域到某一端。这就是后来被广泛称为 non-Hermitian skin effect 的现象。

在 t3=0 的简化情况下,可以通过相似变换把非厄米 SSH 变成普通 SSH,但代价是波函数多了一个指数因子。普通 Bloch 因子 e^{ik} 被替换为:

β = r eik,    r = √ | (t1 - γ/2) / (t1 + γ/2) |

厄米系统里 r=1,所以 β 在单位圆上;非厄米系统里 r 一般不等于 1,波函数会多出指数衰减或指数增长因子。开链能接受的是衰减方向,所以体态堆到边界。

普通 Bloch:|e^{ik}| = 1 非布洛赫:|β| = r e^{ik} β 灰虚线:旧单位圆

5. 广义布里渊区:Cβ 才是真正的体

论文用更一般的方法重新求解开链。令波函数在每个晶胞上按 β^n 变化,把两个 β 解记为 β1(E)β2(E)。在长链极限下,体能谱必须满足:

1(E)| = |β2(E)|

这个条件非常关键。它不是额外假设,而是由开边界条件和长链极限推出的。如果两个 β 的模不同,某一项会在长链中指数压倒另一项,无法形成真正的体谱。

在 t3=0 时,条件给出一个半径为 r 的圆;在 t3 非零时, 不再是圆,而是一条更一般的闭合曲线。这条曲线就是 generalized Brillouin zone。

一句话:普通布里渊区是假定体态是纯相位波;广义布里渊区承认体态带指数衰减。非厄米开链的拓扑必须在后者上计算。

6. 非布洛赫 winding number

有了 β,就把 Bloch Hamiltonian 改写成 non-Bloch Hamiltonian。做法非常直接:

eik → β,    e-ik → β-1

对 t3=0 模型,论文写成:

H(β) = (t1 - γ/2 + β t2) σ- + (t1 + γ/2 + β-1 t2) σ+

非厄米矩阵有左右本征矢,所以论文用 biorthogonal 的方式定义 Q 矩阵。手性对称性让 Q 是非对角形式,里面有一个复函数 q。最终的拓扑不变量是:

W = ( i / 2π ) ∮ q-1 dq

这就是 non-Bloch winding number。它和普通 winding number 形式很像,但积分路径不再是单位圆,而是广义布里渊区

项目普通 Bloch 拓扑非布洛赫拓扑
波函数因子e^{ik}β = r e^{ik} 或更一般的 Cβ 曲线
积分路径单位圆广义布里渊区 Cβ
适用边界周期边界体谱开边界体谱
是否能预言零模非厄米情形一般失败论文中正确预言稳健零模数

7. 零模到底怎么被预言

论文得到的结论是:2W 计数左右两端稳健零能模式的总数。比如在 t2=1, γ=4/3 的例子里:

当 |t1| < √(t22 + (γ/2)2) 时,有两个零模;否则没有。

还有一个很细但很重要的点:普通 Bloch gap closing 点并不一定是“零模出现/消失”的点,而可能只是零模从左端迁移到右端的点。论文指出,在某些参数区间,两个零模都在左端;中间区间一个在左端一个在右端;另一侧两个都在右端。总数不变,只是位置迁移。

-√(t2²+(γ/2)²) +√(t2²+(γ/2)²) Bloch gap closing Bloch gap closing 非布洛赫拓扑区:零模总数 = 2 红线是真正零模产生/消失的开边界相变;蓝虚线可对应零模左右迁移。

8. 为什么以前会被“正常化条件”误导

论文专门讨论了一个容易犯的错误。找左端零模时,人们通常要求波函数可归一化,即某个 |β| < 1。把临界条件写成 |β| = 1,会得到与 Bloch gap closing 一致的结果,看起来很合理。

但这正是陷阱:它偷偷假设零模在相变时并入的是普通 Bloch 体态,也就是 |β|=1 的体态。非厄米开链的体态实际在 |β|=r 上。所以真正的并入条件不是 |β|=1,而是:

zero| = r

这一步改正,直接把相变点从普通 Bloch 点移到 ±√(t2²+(γ/2)²)

9. t3 非零:广义布里渊区不一定是圆

论文主文后半部分加入 t3,使模型更一般。此时 β 方程变成四次方程,存在四个根。体谱由某一对根的模相等条件决定,即 |βi(E)|=|βj(E)|

这时 不再是半径固定的圆,不同能量对应不同 |β|,但沿这条 Cβ 曲线计算 non-Bloch winding number,仍然正确预言零模数。这个部分非常重要,因为它说明论文不是只发现了一个可解特例,而是给出了可推广的程序。

推广程序:写出开链 bulk 方程 → 求 β 根 → 用模相等条件找 Cβ → 用 H(β) 构造 Q → 沿 Cβ 算 winding number → 预言开边界零模。

10. 这篇论文的范式意义

这篇论文的贡献可以分成三层。

  1. 物理现象层:指出非厄米开链中体态会整体局域到边界,也就是非厄米趋肤效应。
  2. 数学对象层:用 β 和 Cβ 替代普通动量 k 与单位圆,建立广义布里渊区。
  3. 拓扑原则层:提出 non-Bloch winding number,并建立非布洛赫体边对应。

如果说普通拓扑能带理论的基础是“体态是 Bloch 波”,那么这篇论文告诉我们:非厄米系统里,这个基础要换成“体态是 β 波”。这就是它被视为范式级突破的原因。

一张心智地图

开边界谱不同于周期谱 所有体态指数贴边 普通 Bloch拓扑数失效 β = r e^{ik}Cβ 广义布里渊区 H(β) 与 Q(β)非布洛赫 winding 正确预言零模数 核心:不是边缘态特殊,而是“体态”已经变了。

来源与说明

本文基于 arXiv 版本和 PRL 公开页反复阅读后整理,所有示意图为重新绘制的解释图,并非论文原图复刻。公式为方便阅读做了排版简化。