2025-11-23T05:04:16.353666

High-Winding-Number Zero-Energy Edge States in Rhombohedral-Stacked Su-Schrieffer-Heeger Multilayers

Lu, Zhou, Cheng et al.
We study the topological properties of rhombohedral-stacked N-layer Su-Schrieffer-Heeger networks with interlayer coupling. We find that these systems exhibit $2N$-fold degenerate zero-energy edge states with winding number $W=N$, providing a direct route to high-winding-number topological phases where $W$ equals the layer number. Using effective Hamiltonian theory and Zak phase calculations, we demonstrate that the winding number scales linearly with $N$ through a layer-by-layer topological amplification mechanism. We introduce the Wigner entropy as a novel detection method for these edge states, showing that topological boundary states exhibit significantly enhanced Wigner entropy compared to bulk states. Our results establish rhombohedral stacking as a systematic approach for engineering high-winding-number topological insulators with potential applications in quantum information processing.
academic

High-Winding-Number Zero-Energy Edge States in Rhombohedral-Stacked Su-Schrieffer-Heeger Multilayers

基本信息

  • 论文ID: 2511.08167
  • 标题: High-Winding-Number Zero-Energy Edge States in Rhombohedral-Stacked Su-Schrieffer-Heeger Multilayers
  • 作者: Feng Lu, Ao Zhou, Shujie Cheng, Gao Xianlong (浙江师范大学物理系)
  • 分类: cond-mat.dis-nn (凝聚态物理-无序系统与神经网络)
  • 发表时间: 2025年11月11日提交至arXiv
  • 论文链接: https://arxiv.org/abs/2511.08167v1

摘要

本文研究了具有层间耦合的菱形堆叠N层Su-Schrieffer-Heeger (SSH)网络的拓扑性质。研究发现这些系统展现出2N重简并的零能量边缘态，其缠绕数W=N，提供了实现高缠绕数拓扑相的直接途径。通过有效哈密顿量理论和Zak相位计算，证明了缠绕数通过逐层拓扑放大机制与层数N呈线性标度关系。引入Wigner熵作为检测这些边缘态的新方法，表明拓扑边界态相比体态展现出显著增强的Wigner熵。研究结果确立了菱形堆叠作为工程化高缠绕数拓扑绝缘体的系统方法，在量子信息处理中具有潜在应用价值。

研究背景与动机

1. 核心问题

本文旨在解决如何系统地构造高缠绕数(|W| ≫ 2)的一维拓扑相这一关键问题。传统SSH模型通常只能实现W = ±1的拓扑相，而高缠绕数拓扑相能够承载更丰富的零能量边缘态流形。

2. 问题重要性

高缠绕数拓扑相具有多重重要意义：

  • 高度简并边缘模式：可以承载更多的拓扑保护边界态
  • 分数电荷局域化：可能实现新奇的量子现象
  • 增强的抗无序鲁棒性：对系统扰动具有更强的抵抗能力
  • 量子信息应用：为拓扑量子计算提供潜在平台

3. 现有方法的局限性

  • 大多数SSH模型研究集中在低缠绕数相(W = 0, ±1)
  • 缺乏系统构造高缠绕数拓扑相的通用方法
  • 高缠绕数相的表征和检测手段不足
  • 对层数与拓扑不变量之间关系的理解不够深入

4. 研究动机

受石墨烯多层系统启发，其中菱形堆叠的石墨烯与Haldane模型层耦合后展现出高Chern数量子反常霍尔态(|C| = N+1)。本文探索类似的层依赖拓扑放大机制是否可以在SSH模型中实现，从而建立层数N与缠绕数W之间的直接对应关系。

核心贡献

  1. 发现W=N的线性标度关系：证明了菱形堆叠N层SSH系统的缠绕数严格等于层数N，建立了层数与拓扑不变量的一一对应关系。
  2. 揭示逐层拓扑放大机制：通过有效哈密顿量理论和Schrieffer-Wolff变换，解析证明了弱层间耦合(t'⊥ ≪ t⊥)保持拓扑性质的同时实现缠绕数的线性增长。
  3. 观察到2N重简并零能态：发现系统在拓扑非平凡区域(ω > ν)展现2N重简并的零能量边缘态，简并度随层数线性增长。
  4. 引入Wigner熵检测方法：提出使用Wigner相空间分布的熵作为识别拓扑边界态的新诊断工具，边缘态的Wigner熵显著高于体态。
  5. 建立系统工程化方案：确立菱形堆叠为构造高缠绕数拓扑绝缘体的通用策略，可推广到其他一维拓扑系统。

方法详解

任务定义

研究菱形堆叠N层SSH网络的拓扑性质，具体包括：

  • 输入：层数N、层内跃迁参数(ν, ω)、层间耦合参数(t⊥, t'⊥)
  • 输出：能谱、缠绕数W、零能态简并度、Wigner熵分布
  • 约束：弱层间耦合条件 t'⊥ ≪ t⊥，拓扑非平凡条件 ω > ν

模型架构

1. 系统哈密顿量

总哈密顿量由三部分组成：

H^=∑n=1NH^n+∑n=2N−1H^n,n+1⊥+H^H⊥\hat{H} = \sum_{n=1}^{N} \hat{H}_n + \sum_{n=2}^{N-1} \hat{H}^{\perp}_{n,n+1} + \hat{H}^{\perp}_H

单层SSH哈密顿量： H^n=∑j,nνcj,A,n†cj,B,n+ωcj,B,n†cj+1,B,n+H.c.\hat{H}_n = \sum_{j,n} \nu c^\dagger_{j,A,n}c_{j,B,n} + \omega c^\dagger_{j,B,n}c_{j+1,B,n} + \text{H.c.}

其中ν为胞内跃迁强度，ω为胞间跃迁强度。

层间耦合：

  • 标准层间耦合(n ≥ 2)： H^n,n+1⊥=t⊥∑jcj,B,n†cj,A,n+1+H.c.\hat{H}^{\perp}_{n,n+1} = t_{\perp} \sum_j c^\dagger_{j,B,n}c_{j,A,n+1} + \text{H.c.}
  • 第一层与其余层的弱耦合： H^H⊥=t⊥′∑jcj,B,0†cj,A,1+H.c.\hat{H}^{\perp}_H = t'_{\perp} \sum_j c^\dagger_{j,B,0}c_{j,A,1} + \text{H.c.}

2. 动量空间表示

通过傅里叶变换，Bloch哈密顿量为2N×2N三对角矩阵：

0 & f_k & 0 & 0 & \cdots \\ f^*_k & 0 & t'_{\perp} & 0 & \cdots \\ 0 & t'_{\perp} & 0 & f_k & \cdots \\ 0 & 0 & f^*_k & 0 & \cdots \\ \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}$$ 其中 $f_k = \nu + \omega e^{-ik}$。 #### 3. 有效哈密顿量理论 采用低能理论推导(N-1)层SSH网络的有效哈密顿量： $$H^{\text{SSH}}_{N-1}(k) = \begin{pmatrix} 0 & F_{N-1}(k) \\ F^*_{N-1}(k) & 0 \end{pmatrix}$$ 其中 $F_{N-1}(k) = (-t_{\perp})^{2-N}f_k^N$。 当加入第一层SSH时，总有效哈密顿量为： $$H^{\text{eff}}(k) = \begin{pmatrix} 0 & f_k & 0 & 0 \\ f^*_k & 0 & t'_{\perp} & 0 \\ 0 & t'_{\perp} & 0 & F_{N-1}(k) \\ 0 & 0 & F^*_{N-1}(k) & 0 \end{pmatrix}$$ #### 4. Zak相位与缠绕数计算 通过Schrieffer-Wolff变换，有效哈密顿量可标度为： $$H^{\text{eff}}_{N-1}(k) \propto \begin{pmatrix} 0 & e^{-i(N-1)\varphi(k)} \\ e^{i(N-1)\varphi(k)} & 0 \end{pmatrix}$$ 其中 $\varphi(k) = \arctan\left(-\frac{\sin k}{\nu/\omega + \cos k}\right)$。 负能带波函数为： $$|u_-(k)\rangle = \frac{1}{\sqrt{2}}\begin{pmatrix} e^{-i(N-1)\varphi} \\ -1 \end{pmatrix}$$ Zak相位计算： $$\varphi_{\text{Zak}} = i\oint_{-\pi}^{\pi} dk \langle u_-(k)|\partial_k|u_-(k)\rangle = \begin{cases} 0, & \omega < \nu \\ (N-1)\pi, & \omega > \nu \end{cases}$$ 因此缠绕数 $W_{N-1} = \varphi_{\text{Zak}}/\pi = N-1$，加上第一层贡献，总缠绕数 **W = N**。 ### 技术创新点 #### 1. 菱形堆叠几何的巧妙利用 - 通过特殊的层间耦合拓扑结构实现拓扑性质的叠加 - 弱耦合参数t'⊥的引入使第一层可视为独立贡献W=1 - 强耦合参数t⊥确保其余层形成整体贡献W=N-1 #### 2. 递归证明策略 将N层系统的证明转化为(N-1)层系统的证明，建立归纳关系： - 单层SSH：W = 1（已知） - 假设(N-1)层：W = N-1 - 证明N层：W = N #### 3. Wigner相空间方法 - **预处理**：将N层的概率幅映射到单有效层 $$\tilde{\varphi}_j = \sum_{n=1}^{N} \frac{1}{N}[\phi^2_{j,A,n} + \phi^2_{j,B,n}]$$ - **Wigner函数**： $$W(x,p) = \frac{1}{2\pi\hbar}\int_{-\infty}^{\infty} \langle x-\frac{y}{2}|\rho|x+\frac{y}{2}\rangle e^{-ipy/\hbar}dy$$ - **Wigner熵**： $$W_s = -\int\int W(x,p)\ln|W(x,p)|dxdp$$ 该方法的创新在于将拓扑边界态的实空间局域化特征转化为相空间熵的增强信号。 ## 实验设置 ### 数值参数 本文主要进行数值模拟研究，参数设置如下： **系统参数**： - 晶格长度：L = 100（图2）或L = 200（图3） - 层数：N = 1（单层对比）和N = 5（多层演示） - 层内跃迁比：ω/ν在0到3之间变化 **层间耦合参数**： - 弱耦合：t'⊥ = 0.1ν - 强耦合：t⊥ = ν - 满足弱耦合条件：t'⊥ ≪ t⊥ ### 计算方法 1. **能谱计算**：通过Schmidt正交化方法数值求解Bloch哈密顿量的本征值和本征态 2. **零能态识别**：统计能量在[-ε, ε]范围内的本征态数量（ε为小量） 3. **Wigner分布计算**：对每个本征态计算其在相空间的Wigner函数 4. **Wigner熵计算**：对Wigner分布进行数值积分得到熵值 ### 对比场景 - **拓扑非平凡相**：ω = 1.5ν（ω > ν） - **拓扑平凡相**：ω = 0.5ν（ω < ν） - **单层系统**：N = 1作为基准 - **多层系统**：N = 5展示高缠绕数效应 ## 实验结果 ### 主要结果 #### 1. 能谱结构（图2） **单层SSH模型(N=1)**： - 图2(a)显示在ω > ν时出现零能态 - 图2(b)在ω = 2ν时展现**2重简并**零能态，符合W = 1 **五层堆叠系统(N=5)**： - 图2(c)显示更丰富的能带结构，零能态在ω > ν区域清晰可见 - 图2(d)在ω = 2ν时展现**10重简并**(2N = 10)零能态，验证W = 5 **关键发现**：零能态简并度严格遵循2N规律，直接验证了缠绕数W = N的理论预测。 #### 2. Wigner相空间分析（图3） **拓扑非平凡相(ω = 1.5ν)**： - 图3(a)：第1000个本征态的Wigner分布W(x,p)在x方向呈现明显的**边界局域化**，在p方向呈高斯分布 - 图3(b)：Wigner熵Ws随能级指数的变化，零能态（插图放大）显示**显著增强的熵值** **拓扑平凡相(ω = 0.5ν)**： - 图3(c)：Wigner分布在x方向呈高斯分布，在p方向呈离散分布，表明为**体态** - 图3(d)：所有态的Wigner熵保持一致的低值，无明显峰值 **定量对比**： - 边缘态的Wigner熵比体态高出约**一个数量级** - 提供了清晰的相空间诊断信号，可区分拓扑边界态与体态 #### 3. 缠绕数验证 通过Zak相位计算解析证明： $$W = \frac{\varphi_{\text{Zak}}}{\pi} = \begin{cases} 0, & \omega < \nu \text{ (平凡相)} \\ N, & \omega > \nu \text{ (拓扑相)} \end{cases}$$ 数值结果与解析预测完全一致，确认了**W = N的线性标度关系**。 ### 实验发现 1. **拓扑放大的层依赖性**：缠绕数与层数呈严格线性关系，每增加一层贡献ΔW = 1 2. **简并度的精确预测**：零能边缘态的简并度2N可通过层数直接确定 3. **弱耦合的关键作用**：t'⊥ ≪ t⊥的条件确保第一层的独立性，是实现线性标度的必要条件 4. **Wigner熵的普适性**：作为拓扑边界态的诊断工具，Wigner熵在不同层数和参数下均展现出一致的区分能力 5. **鲁棒性**：拓扑性质在合理的参数范围内保持稳定 ## 相关工作 ### SSH模型及其扩展 1. **原始SSH模型**[6]：描述聚乙炔中的电荷密度波，是一维拓扑系统的范式模型 2. **二维SSH晶格**[16-18]：展现高阶拓扑相，具有角局域化零模 3. **SSH双层研究**：探索层间耦合拓扑态，介于1D和2D行为之间 4. **扩展SSH模型**[19-21]：研究无序、相互作用对高缠绕数相的影响 ### 多层拓扑系统 1. **石墨烯多层系统**[25-29]： - 菱形堆叠石墨烯与Haldane模型耦合实现|C| = N+1的高Chern数 - 层间跃迁打破高阶Dirac点简并 - 拓扑邻近效应打开体能隙并保持非平凡边界输运 2. **拓扑异质结构**：通过堆叠诱导的杂化实现超越单层极限的拓扑数 ### 拓扑态检测方法 1. **传统方法**：Zak相位测量[14,31-38]、声学信号[32,37]、光学信号[14,33,34,36]、电学信号[35,38] 2. **量子层析技术**：光学零差层析可直接测量Wigner函数[39-42] 3. **本文创新**：首次将Wigner熵用于拓扑边界态识别 ### 本文优势 1. **系统性**：提供了构造任意高缠绕数拓扑相的通用方案 2. **解析可解性**：通过有效理论给出缠绕数的解析表达式 3. **新检测方法**：Wigner熵提供了互补的相空间视角 4. **可推广性**：方法可扩展到其他一维拓扑模型 ## 结论与讨论 ### 主要结论 1. **W = N的严格线性关系**：菱形堆叠N层SSH系统的缠绕数精确等于层数，建立了层数与拓扑不变量的一一对应 2. **2N重简并零能态**：系统在拓扑非平凡区域(ω > ν)展现2N重简并的零能量边缘态，简并度由菱形堆叠几何保护 3. **逐层拓扑放大机制**：通过有效哈密顿量理论和Zak相位计算证明，弱层间耦合保持拓扑性质的同时实现缠绕数的线性增长 4. **Wigner熵诊断方法**：边缘态相比体态展现显著增强的Wigner熵，提供了识别拓扑边界态的鲁棒实验信号 5. **普适工程化策略**：确立菱形堆叠为构造高缠绕数拓扑绝缘体的系统方法，可推广到其他一维拓扑系统 ### 局限性 1. **弱耦合假设**：要求t'⊥ ≪ t⊥，这限制了第一层与其余层的相互作用强度，可能影响实际系统的实现 2. **无相互作用近似**：模型假设无自旋和粒子间相互作用，实际系统中相互作用可能改变拓扑性质 3. **理想边界条件**：数值模拟采用开边界条件，实际器件中的边界效应可能更复杂 4. **Wigner熵的计算成本**：需要对每个本征态计算相空间分布，计算复杂度较高 5. **实验实现挑战**：虽然SSH模型已在多种平台实现，但精确控制多层菱形堆叠结构和弱耦合参数仍具挑战性 ### 未来方向 1. **相互作用效应**：研究粒子间相互作用对高缠绕数拓扑相的影响 2. **无序鲁棒性**：系统研究无序对2N重简并零能态的保护机制 3. **动力学性质**：探索高缠绕数边缘态的输运性质和时间演化 4. **其他堆叠方式**：研究AB堆叠、Bernal堆叠等其他几何结构的拓扑性质 5. **实验验证**：在冷原子光晶格、光子晶体或超导量子电路中实现并测量 6. **量子信息应用**：利用高度简并边缘态构建拓扑量子比特和量子存储 ## 深度评价 ### 优点 #### 1. 理论严谨性 - **解析证明**：不仅提供数值结果，还通过有效哈密顿量理论和Zak相位给出严格的解析证明 - **递归论证**：将N层问题归约为(N-1)层问题的证明策略优雅且有说服力 - **自洽性检验**：2N重简并与W = N的关系通过体-边对应原理相互验证 #### 2. 方法创新性 - **Wigner熵的首次应用**：开创性地将量子层析技术用于拓扑态检测，提供了新的实验测量途径 - **相空间视角**：将实空间的边界局域化转化为相空间熵的增强，物理图像清晰 - **系统构造方案**：提供了实现任意高缠绕数的通用方法，具有重要指导意义 #### 3. 结果的普适性 - **线性标度律**：W = N的简洁关系具有普适性，类似于石墨烯多层系统的Chern数标度 - **可推广性**：方法不局限于SSH模型，可扩展到其他一维拓扑系统 - **多平台实现潜力**：SSH模型已在多种实验平台实现，本文结果具有实验可行性 #### 4. 物理图像清晰 - **直观的几何理解**：菱形堆叠的层间耦合结构与拓扑放大机制的关系清晰 - **可视化效果好**：Wigner分布图直观展示边缘态与体态的差异 - **参数依赖明确**：清楚说明各参数(ν, ω, t⊥, t'⊥)对拓扑性质的影响 ### 不足 #### 1. 实验实现的挑战 - **弱耦合条件**：t'⊥ ≪ t⊥的精确控制在实验中可能困难 - **多层制备**：菱形堆叠的精确几何排列需要高精度的实验技术 - **参数调控**：需要独立调控层内和层间跃迁参数 #### 2. 理论假设的限制 - **无相互作用近似**：忽略了粒子间相互作用，可能在强关联系统中失效 - **完美晶格假设**：未考虑无序、缺陷等实际因素 - **零温近似**：有限温度效应未被讨论 #### 3. 分析的深度 - **动力学性质缺失**：仅研究了静态能谱，未涉及边缘态的输运和动力学 - **鲁棒性分析不足**：对参数偏离、无序等扰动的鲁棒性缺乏定量分析 - **相互作用效应**：未讨论相互作用对高缠绕数相的可能影响 #### 4. Wigner熵方法的局限 - **计算复杂度**：对大系统计算Wigner函数的数值成本较高 - **预处理依赖**：将多层映射到单层的预处理可能损失部分信息 - **阈值选择**：区分边缘态与体态的Wigner熵阈值需要经验确定 ### 影响力 #### 1. 对领域的贡献 - **理论突破**：建立了多层拓扑系统中层数与拓扑不变量的定量关系 - **方法学贡献**：Wigner熵为拓扑态检测提供了新工具 - **概念推广**：将石墨烯多层系统的拓扑放大机制推广到一维系统 #### 2. 实用价值 - **量子信息应用**：高度简并边缘态可用于拓扑量子比特和量子存储 - **拓扑电路设计**：为设计具有特定拓扑性质的人工晶格提供指导 - **实验指导**：为冷原子、光子、声学等平台的实验提供理论预测 #### 3. 可复现性 - **参数明确**：所有数值模拟参数清晰给出 - **方法详细**：计算方法和理论推导步骤完整 - **解析结果**：提供可验证的解析表达式 - **潜在挑战**：Wigner函数计算需要适当的数值技巧 ### 适用场景 #### 1. 理论研究 - 研究一维拓扑系统的高缠绕数相 - 探索层数依赖的拓扑相变 - 发展多层拓扑系统的有效理论 #### 2. 实验平台 - **冷原子光晶格**：已有SSH模型实现，可扩展到多层结构 - **光子晶体**：波导阵列可实现可调的层间耦合 - **超导量子电路**：可编程性强，适合实现复杂耦合结构 - **声学系统**：机械超材料可实现菱形堆叠几何 #### 3. 应用领域 - **拓扑量子计算**：利用高度简并边缘态编码量子信息 - **拓扑光子学**：设计具有多通道边界输运的光子器件 - **量子模拟**：模拟高维拓扑系统的低能有效理论 #### 4. 推广方向 - 其他一维拓扑模型（如Kitaev链、Rice-Mele模型）的多层推广 - 二维拓扑系统的层依赖拓扑不变量研究 - 拓扑半金属的多层堆叠效应 ## 参考文献（关键文献） 1. **[6] W.-P. Su et al., PRL 42, 1698 (1979)** - SSH模型的原始论文 2. **[25] Y. Zhao et al., [Submitted] (2025)** - 石墨烯多层系统的高Chern数研究，启发本文工作 3. **[39] A. I. Lvovsky & M. G. Raymer, RMP 81, 299 (2009)** - 连续变量光学量子态层析综述 4. **[45] Y. Hatsugai, PRL 71, 3697 (1993)** - 体-边对应原理的经典文献 5. **[47] J. Zak, PRL 62, 2747 (1989)** - Zak相位的原始定义 --- **总体评价**：这是一篇高质量的凝聚态理论物理论文，系统地研究了菱形堆叠SSH多层系统的拓扑性质。论文的核心贡献在于建立了缠绕数与层数的线性关系(W=N)，并引入Wigner熵作为新的检测方法。理论推导严谨，数值结果清晰，物理图像直观。尽管在实验实现和相互作用效应方面存在一些局限，但为高缠绕数拓扑相的工程化提供了重要的理论基础和实验指导，具有较高的学术价值和应用潜力。