在动力学阻挫体系中可调谐实现平带与例外点：以非厄米 Creutz 梯子为例

论文信息

标题: Tunably realizing flat-bands and exceptional points in kinetically frustrated systems: An example on the non-Hermitian Creutz ladder

作者: Debashish Dutta, Sayan Choudhury

发布日期: 2025-12-23

arXiv ID: 2512.20614v1

PDF 链接: 下载 PDF

3 分钟速览

研究问题：在非厄米系统中，动力学阻挫诱导的平带如何与例外点相互作用？周期和开放边界条件下能谱结构有何本质区别？
核心方法：通过局域基变换将非厄米 Creutz 梯子映射为两个解耦的非厄米 SSH 链，分别在周期边界条件 (PBC) 和开放边界条件 (OBC) 下进行严格对角化。
关键结果：发现参数空间中由例外线分隔出完全实谱、完全虚谱和复谱三种区域，并展示平带可以表现为厄米简并点（diabolical point）或非厄米例外点（exceptional flat band）两种形式。
主要局限：仅考虑 $\Delta t = \Delta \gamma = 0$ 的平衡链内跃迁情形，且未讨论相互作用、无序效应和实验实现的具体方案。
适合读者：从事非厄米拓扑物理、平带系统、凝聚态理论和量子动力学研究的学者和研究生。

论文背景和研究动机

非厄米哈密顿量已成为描述开放、耗散和非互易系统的重要框架 [bender]。与保证实能谱和正交本征态的厄米算符不同，非厄米模型展现出复能谱、实-复能谱转变和例外点（exceptional point）等独特现象 [EP1, EP2, EP5]。例外点处，本征值和本征态同时简并，哈密顿量变为缺陷矩阵，这在厄米系统中没有对应物。

与此同时，平带系统因几何阻挫导致的动能猝灭和宏观能谱简并而备受关注 [fb1, fb2, fb3]。Creutz 梯子正是这类平带系统的范例：在精细调节参数 $t_1 = t_2 = \pm t_0$ 时，晶格几何引发的破坏性干涉完全抑制了动能色散，产生两个严格平坦的能带，波包呈现仅占据相邻面片的紧凑局域化（compact localization）[creutz1, creutz2]。

将非厄米性引入 Creutz 梯子提出了一系列基本问题：非互易跃迁如何重组平带简并？平带是保持为厄米简并点还是获得例外点特征？边界敏感性如何体现在阻挫几何中？这些问题不仅具有概念上的重要性，也为在光子、力学和电路平台上实现非厄米平带物理提供了理论指引。本文的作者正是基于这些动机，构建了一个具有一般非互易跃迁的严格可解非厄米 Creutz 梯子模型。

核心方法和技术细节

模型构建与 PBC 能谱

作者从标准 Creutz 梯子出发，将每个互易跃迁过程配以相应的非互易分量。系统的实空间哈密顿量形如公式 (1)，包含沿上、下链的跃迁项和链间耦合项。每条腿上的跃迁包含互易部分（ $t_1, t_2$ ）和非互易部分（ $\gamma_1, \gamma_2$ ），链间耦合也分解为 $t_0$ 和 $\gamma_0$ 。

在周期边界条件下，通过傅里叶变换到动量空间，哈密顿量分块对角化为 $2 \times 2$ 矩阵 $h(k)$ （公式 2）。引入重整化参数 $\bar{t} = (t_1+t_2)/2$ , $\bar{\gamma} = (\gamma_1+\gamma_2)/2$ , $\Delta t = (t_1-t_2)/2$ , $\Delta \gamma = (\gamma_1-\gamma_2)/2$ 后，能谱的解析表达式如公式 (3) 所示。作者聚焦于平衡链内跃迁 ( $\Delta t = \Delta \gamma = 0$ ) 的情形，并固定 $\bar{t}=1$ 作为能量标度。

PBC 能谱的一个重要特征是存在一条精细调节线 $t_0 \gamma_0 = \bar{t} \bar{\gamma}$ ，沿线能谱完全为实数（当 $\bar{t} > \gamma_0$ 时）。偏离这条线后，能谱变为复数，但这一实-复转变并非通过例外点，而是通过厄米简并点（diabolical point）实现（见图 4（a,b））。

精确映射到非厄米 SSH 链

论文的核心技术突破在于通过局域基变换（公式 4）将梯子从 $(a, b)$ 基矢转换到 $(w, \bar{w})$ 基矢。在 $\Delta t = \Delta \gamma = 0$ 的条件下，系统严格解耦为两个独立的非厄米 SSH 链 $H_1$ 和 $H_2$ （见图 2），每个链的哈密顿量如公式 (5) 所示，其中重定义参数 $g = \bar{t}+t_0, g' = \bar{t}-t_0, f = \bar{\gamma}+\gamma_0, f' = \bar{\gamma}-\gamma_0$ 。

这一映射的重要意义在于：非厄米 SSH 链在开放边界条件下可以通过非幺正的相似性变换（imaginary gauge transformation, IGT）转化为厄米 SSH 链 [hn1igt1, hn2igt2, igt3]。相似性矩阵 $S_1$ 的结构如公式 (6) 所示，其元素由比率 $\sqrt{(f'+g')/(f'-g')}$ 和 $\sqrt{(f+g)/(f-g)}$ 递推确定。

经过 IGT 后，OBC 能谱在热力学极限下的解析形式为：

E^{\text{OBC}}_{\pm}(q) = \pm \sqrt{u^2 + v^2 + 2uv \cos q}

其中 $u = \sqrt{g^2 - f^2}$ （胞间跃迁）和 $v = \sqrt{g'^2 - f'^2}$ （胞内跃迁）。尽管此色散关系形式上与厄米 SSH 模型相同，但 $u$ 和 $v$ 可以为复数，因而能量本征值可以是实数、虚数或复数。

创新点和贡献

参数空间的精确相图

本文最重要的创新在于借助 IGT 方法得到了 OBC 下的完整参数空间相图。通过引入谱密度函数 $M$ （公式 8），作者将 $(t_0, \bar{\gamma})$ 平面划分为三个具有本质不同能谱特征的区域（见图 4（d-f））：

实能谱区（蓝色）： $u$ 和 $v$ 同为实数
虚能谱区（红色）： $u$ 和 $v$ 同为虚数
复能谱或混合区（彩色渐变）：其他情形

这三个区域由四条例外线（exceptional lines）分隔。实线橙色线对应 $u=0$ （即 $g = \pm f$ ），虚线橙色线对应 $v=0$ （即 $g' = \pm f'$ ）。沿这些例外线，相似性矩阵变为奇异，OBC 链展现出缺陷动力学。在一般点处，哈密顿量约化为 $D_2(0) \oplus J_{L/2-1}(\pm\lambda) \oplus J_{L/2}(\pm\lambda)$ 的形式，其中 $D$ 和 $J$ 分别表示对角矩阵和 Jordan 块矩阵。

三种平板的精细区分

论文揭示并精细区分了非厄米 Creutz 梯子中三种不同的平带形式：

厄米简并点平带：在 $\eta = \pm 1$ （即 $t_0/\bar{t} = \pm 1$ ）处，PBC 能谱的两个带塌缩为无色散的平带，但由于哈密顿量仍可对角化，这些平带对应的简并是 diabolical point 型而非例外点型。
例外平带（exceptional flat band）：沿 $\bar{t} = \pm \gamma_0$ 和 $t_0 = \bar{\gamma}$ 的线，所有本征态在 $E=0$ 处简并，哈密顿量成为由二阶 Jordan 块组成的缺陷矩阵 $\mathcal{H}_{\mathrm{OBC}} \cong [\oplus J_2(0)]^L$ 。这种完全的 EP2 简并构成例外平带 [efb1]，其调谐约束比 Hermitian 情形更为宽松。
三重交汇点：在 $(\pm\gamma_0, \pm\bar{t})$ 处，实线例外线与虚线例外线相交，两个 NH-SSH 链各自退化为 $L$ 阶 Jordan 块 $J_L(0)$ ，整个 OBC 能谱塌缩到复能面的原点。这些点成为实谱、虚谱和复谱三种区域的三重交汇点（图 4（d,e）中的黄色菱形标记）。

非厄米趋肤效应的严格刻画

通过 IGT 矩阵的结构，作者解析推导出本征态的逆局域长度：

\xi^{-1} = \frac{1}{2} \left[ \log(f+g) + \log(f'+g') - \log(f-g) - \log(f'-g') \right]

这一表达式直接给出了波包指数局域化的方向和强度。图 5 展示了平均方向性逆参与比（dIPR）在参数空间的分布，粉色和绿色区域分别对应右端和左端局域的本征态。

体-边对应（BBC）恢复的条件是 IGT 矩阵的有界性，这等价于生长因子为幺模：

t_0^2 + \gamma_0^2 = \bar{t}^2 + \bar{\gamma}^2 \quad \text{和} \quad t_0 \gamma_0 = \bar{t} \bar{\gamma}

这两条线（图 5 中的虚线黑线）的交点正好是 diabolical point 平带点，对应紧凑局域化态。

波包动力学的独特特征

论文通过改进的方向性逆参与比（mIPR, 公式 13）系统研究了波包在 OBC 下的时间演化（图 6）。关键发现包括：

在 diabolical point 平带处，波包呈现振荡型紧凑动力学，周期性地在相邻胞间扩展和收缩 [图 6（iii）]
在例外平带处，紧凑动力学发生质变：波包虽保持空间局域但丧失振荡恢复行为 [图 6（vii,viii）]
沿例外线，波包动力学直接反映静态的 NHSE 极性 [图 6（iv,v）]
在 BBC 恢复点，波包对称扩展，类似厄米行为 [图 6（vi）]

这些结果不仅验证了静态能谱分析得出的结论，更重要的是揭示了平带的起源（正常简并 vs. 例外简并）决定了其动力学特征，这为实验区分不同平带类型提供了清晰的测量方案。

局限与待解决问题

本文的模型设定在几个重要方面存在理想化假设：

平衡跃迁的限制：全文主动取 $\Delta t = \Delta \gamma = 0$ ，这意味着上、下链的跃迁振幅完全对称。论文仅在映射到 NH-SSH 链时简要提及此约束（见公式推导），未深入讨论当 $\Delta t \neq 0$ 或 $\Delta \gamma \neq 0$ 时链间耦合如何改变能谱结构和例外点特征。
无相互作用系统的局限：Creutz 梯子在厄米极限下是被广泛研究的强关联平带平台，可产生超导、铁磁、量子霍尔态等 [fb-sc1, fb-sc2, fb-mag1]。但本文完全未讨论非厄米情形下的相互作用效应，这使得理论距离实际关联电子系统较远。
无序和无序效应的缺失：论文分析了纯理论模型的能谱结构，未涉及无序对外部平带、例外点和 NHSE 的影响。在实际的光子晶格或电路实现中，加工缺陷不可避免。
实验实现路径的模糊性：论文只在结论部分笼统提到可能的光子、力学和电路平台实现，未给出任何具体的参数映射方案或实验可观测量的讨论。例如，在光子波导阵列中如何实现非互易跃迁 $(\gamma_0, \gamma_1, \gamma_2)$ 的具体参数配置，论文未说明。这使得模型与实验的联系仍停留在概念层面。
动力学分析的有限范围：波包动力学研究（第 V 节）仅讨论了单个初态在特定参数点的演化行为，未进行系统性扫描或讨论长时稳态的普适性质（如弛豫速率与例外点阶数的关系）。

这些局限同时也指明了有价值的后续研究方向：建立包含相互作用和非平衡效应的理论框架、研究无序对例外平带稳定性的影响、以及设计具体的实验实现方案和测量协议。