2025-11-15T17:07:12.333568

A One-Particle Density Matrix Framework for Mode-Shell Correspondence: Characterizing Topology in Amorphous Higher-Order Topological Insulators

Martínez, Jezequel, Bardarson et al.
We present a framework for characterizing higher-order topological phases directly from the one-particle density matrix, without any reference to an underlying Hamiltonian. Our approach extends the mode-shell correspondence, originally formulated for single-particle Hamiltonians, to Gaussian states subject to chiral constraints. In this correspondence, the mode index counts topological boundary modes, while the shell index quantifies the bulk topology in a region surrounding the modes, providing a bulk-boundary diagnostic. In one-dimensional topological insulators, the shell index reduces to the local chiral marker, recovering the winding number in the translation-invariant limit. We apply the mode-shell correspondence to a $C_4$-symmetric higher-order topological insulator with a chiral constraint and show that a fractional shell index implies that the higher-order phase is intrinsic. The one-particle density matrix is formulated in real space, so the mode-shell correspondence also applies to models without translation invariance. By introducing structural disorder into the $C_4$-symmetric higher-order insulator, we show that the mode-shell correspondence remains a meaningful diagnostic in amorphous structures. The mode-shell correspondence generalizes to interacting states with a gapped bulk spectrum in the one-particle density matrix, providing a practical and diverse route to characterize higher-order topology from the quantum state itself.
academic

A One-Particle Density Matrix Framework for Mode-Shell Correspondence: Characterizing Topology in Amorphous Higher-Order Topological Insulators

基本信息

  • 论文ID: 2509.03632
  • 标题: A One-Particle Density Matrix Framework for Mode-Shell Correspondence: Characterizing Topology in Amorphous Higher-Order Topological Insulators
  • 作者: Miguel F. Martínez, Lucien Jezequel, Jens H. Bardarson, Thomas Klein Kvorning, Julia D. Hannukainen
  • 分类: cond-mat.mes-hall cond-mat.dis-nn
  • 发表时间: 2024年9月（arXiv提交）
  • 论文链接: https://arxiv.org/abs/2509.03632

摘要

本文提出了一个直接从单粒子密度矩阵表征高阶拓扑相的框架，无需参考底层哈密顿量。该方法将原本为单粒子哈密顿量制定的模式-壳层对应关系扩展到受手性约束的高斯态。在此对应关系中，模式指数计算拓扑边界模式，而壳层指数量化模式周围区域的体拓扑，提供体-边界诊断。在一维拓扑绝缘体中，壳层指数简化为局域手性标记，在平移不变极限下恢复绕数。作者将模式-壳层对应关系应用于具有手性约束的C4C_4对称高阶拓扑绝缘体，并证明分数壳层指数意味着高阶相是内在的。

研究背景与动机

问题定义

拓扑量子态的表征传统上依赖于哈密顿量的动量空间分析，但这种方法在处理无序系统和非晶材料时面临根本性困难。特别是：

  1. 平移不变性的缺失：在非晶材料中，传统的布里渊区概念失效，基于动量空间的拓扑不变量无法定义
  2. 高阶拓扑相的复杂性：高阶拓扑绝缘体的角模式和铰链模式需要新的表征方法
  3. 内在与外在相的区分：需要判断拓扑相是由体拓扑保证（内在）还是依赖于特定边界终止（外在）

研究重要性

该研究解决了凝聚态物理中的重要问题：

  • 为非晶拓扑材料提供了表征工具
  • 建立了直接从量子态分析拓扑性质的方法
  • 扩展了拓扑物理的应用范围到无序系统

现有方法局限性

传统方法如嵌套Wilson循环、量子化多极矩等都依赖于平移对称性，无法应用于：

  • 非晶材料
  • 强无序系统
  • 相互作用系统的某些情况

核心贡献

  1. 理论框架创新：将模式-壳层对应关系从哈密顿量表述推广到单粒子密度矩阵表述
  2. 实空间拓扑表征：提供了不依赖平移不变性的拓扑相表征方法
  3. 内在性判据：建立了基于分数壳层指数判断高阶拓扑相内在性的准则
  4. 非晶系统应用：首次在具有C4C_4对称性的非晶高阶拓扑绝缘体中验证了该方法
  5. 相互作用系统扩展：证明了该方法可扩展到具有能隙单粒子密度矩阵的相互作用系统

方法详解

单粒子密度矩阵框架

对于多体态∣Ψ⟩|\Psi\rangle，Bogoliubov-de Gennes单粒子密度矩阵具有分块形式：

ρ=(ρ~κκ†1−ρ~∗)\rho = \begin{pmatrix} \tilde{\rho} & \kappa \\ \kappa^\dagger & 1-\tilde{\rho}^* \end{pmatrix}

其中ρ~ij=⟨Ψ∣ci†cj∣Ψ⟩\tilde{\rho}_{ij} = \langle\Psi|c_i^\dagger c_j|\Psi\rangle和κij=⟨Ψ∣ci†cj†∣Ψ⟩\kappa_{ij} = \langle\Psi|c_i^\dagger c_j^\dagger|\Psi\rangle。

模式-壳层对应关系

模式指数定义

模式指数定义为： Imode=∑λ4λ(1−λ)⟨λ∣θC∣λ⟩I^{\text{mode}} = \sum_\lambda 4\lambda(1-\lambda)\langle\lambda|\theta C|\lambda\rangle

其中：

  • λ\lambda是限制单粒子密度矩阵ρA\rho_A的本征值
  • ∣λ⟩|\lambda\rangle是对应的本征态
  • θ\theta是投影到子区域Ashell⊂AA_{\text{shell}} \subset A的算符
  • CC是手性对称算符

壳层指数推导

通过迹的循环性质和手性约束{ρA,C}=C\{\rho_A, C\} = C，模式指数可重写为： Ishell=−2Tr(CρA[ρA,θ])I^{\text{shell}} = -2\text{Tr}(C\rho_A[\rho_A, \theta])

在梯度展开下： Ishell=2∑iTr(CρA[ρA,Xi]∂xiθ)I^{\text{shell}} = 2\sum_i \text{Tr}(C\rho_A[\rho_A, X_i]\partial_{x_i}\theta)

一维情况的等价性

在一维系统中，壳层指数等价于局域手性标记： ν(x)=−2tr⟨x∣ρCXρ∣x⟩\nu(x) = -2\text{tr}\langle x|\rho CX\rho|x\rangle

在体区域中：Ishell=∓νI^{\text{shell}} = \mp\nu

高阶拓扑绝缘体中的应用

BBH模型

考虑Benalcazar-Bernevig-Hughes (BBH)模型的实空间哈密顿量： hrr′=−γ(σz+σyτy)δrr′−η2(σz+iσyτz)δr′,r+ax−η2(σyτy+iσyτx)δr′,r+ayh_{rr'} = -\gamma(\sigma_z+\sigma_y\tau_y)\delta_{rr'} - \frac{\eta}{2}(\sigma_z + i\sigma_y\tau_z)\delta_{r',r+a_x} - \frac{\eta}{2}(\sigma_y\tau_y + i\sigma_y\tau_x)\delta_{r',r+a_y}

内在性判据

对于C4C_4对称的高阶拓扑绝缘体：

  • 奇数模式指数：每个边-壳交叉点贡献半整数值，表明内在拓扑相
  • 偶数模式指数：整数值贡献，无法区分内在和外在相

实验设置

晶体系统参数

  • 正方格子线性尺寸：L=20L = 20
  • BBH模型参数：γ=0.5\gamma = 0.5, η=1\eta = 1
  • 限制区域AA：下左四分之一格子(LA=10L_A = 10)

非晶系统生成

  1. 生成小正方格子(L~=L/2\tilde{L} = L/2)
  2. 通过正态分布引入无序度ww（非晶度参数）
  3. 创建四个拷贝并施加C4C_4旋转对称
  4. 设置参数：w=0.1w = 0.1, 截止距离r0=1.3r_0 = 1.3

壳层函数设计

θ(x,y)=1−[1+e−c(x2+y2−ω)]−1\theta(x,y) = 1 - [1 + e^{-c(\sqrt{x^2+y^2}-\omega)}]^{-1}

其中ω=LA/2\omega = L_A/2确定壳层位置，c=LA/2c = L_A/2（晶体）或c=2LA/3c = 2L_A/3（非晶）控制过渡锐度。

实验结果

晶体BBH模型结果

  • 模式指数：Imode=0.998≈1I^{\text{mode}} = 0.998 \approx 1
  • 壳层指数：Ishell=0.998≈1I^{\text{shell}} = 0.998 \approx 1
  • 四个λ=1/2\lambda = 1/2本征值：对应四个手性角模式
  • 壳层指数分布：仅在边界-壳层交叉点有贡献，每点贡献1/21/2

非晶BBH模型结果

  • 模式指数：Imode=0.996≈1I^{\text{mode}} = 0.996 \approx 1
  • 壳层指数：Ishell=0.996≈1I^{\text{shell}} = 0.996 \approx 1
  • 结构无序影响：更多本征值出现在(0,1)(0,1)区间内
  • 体贡献出现：壳层指数在体区域有正负贡献，但平均为零

关键发现

  1. 量化精度：即使在非晶系统中，模式-壳层对应关系仍保持良好的量化
  2. 内在性确认：分数壳层指数(1/21/2)确认了BBH模型的内在高阶拓扑性质
  3. 鲁棒性：该方法对结构无序具有鲁棒性

相关工作

实空间拓扑标记

  • 局域手性标记：Hannukainen等人提出的一维实空间拓扑不变量
  • Bott指数：用于无序系统的实空间拓扑标记
  • 其他实空间方法：包括谱局域化、局域Chern标记等

高阶拓扑绝缘体

  • 嵌套Wilson循环：传统的高阶拓扑表征方法
  • 多极矩：基于Wannier中心的拓扑不变量
  • 对称性指示器：基于能带表示的拓扑分类

非晶拓扑材料

  • 非晶Chern绝缘体：最早的非晶拓扑材料研究
  • 非晶高阶拓扑绝缘体：近期发展的研究方向

结论与讨论

主要结论

  1. 理论突破：成功将模式-壳层对应关系推广到量子态的单粒子密度矩阵表述
  2. 实用性验证：在晶体和非晶高阶拓扑绝缘体中验证了方法的有效性
  3. 内在性判据：建立了基于分数壳层指数的内在高阶拓扑相判据
  4. 广泛适用性：方法可扩展到相互作用系统和任意维度

局限性

  1. 手性约束要求：方法仅适用于具有手性对称性的系统
  2. 能隙要求：需要单粒子密度矩阵具有能隙
  3. 壳层选择：壳层区域的选择需要仔细考虑相关长度尺度
  4. 计算复杂度：对于大系统，单粒子密度矩阵的对角化可能计算昂贵

未来方向

  1. 其他对称性类别：扩展到其他Altland-Zirnbauer对称性类别
  2. 三维系统：应用到三维高阶拓扑绝缘体
  3. 强相互作用系统：进一步探索相互作用系统中的应用
  4. 实验验证：在真实材料中验证该方法

深度评价

优点

  1. 理论创新性强：首次将模式-壳层对应关系推广到密度矩阵表述，为非晶拓扑材料研究开辟新途径
  2. 方法普适性好：不依赖平移不变性，可应用于广泛的系统类型
  3. 物理图像清晰：模式指数和壳层指数的对应关系提供直观的体-边界关系理解
  4. 数值验证充分：在晶体和非晶系统中都进行了详细的数值验证
  5. 内在性判据新颖：提出了判断高阶拓扑相内在性的新方法

不足

  1. 对称性限制：方法局限于具有手性对称性的系统，限制了应用范围
  2. 计算成本：大系统中密度矩阵对角化的计算复杂度较高
  3. 参数敏感性：壳层函数θ\theta的选择可能影响结果的准确性
  4. 实验联系有限：缺乏与实际材料测量的直接联系

影响力

  1. 学术贡献：为拓扑物理理论发展提供了重要工具
  2. 应用前景：为非晶拓扑材料的设计和表征提供理论基础
  3. 方法推广：可能启发其他实空间拓扑表征方法的发展
  4. 跨学科价值：方法可能在量子信息、材料科学等领域有应用

适用场景

  1. 非晶拓扑材料：无序度较高但保持某些对称性的系统
  2. 有限尺寸系统：边界效应重要的小尺寸拓扑系统
  3. 相互作用系统：单粒子密度矩阵有隙的强关联系统
  4. 数值研究：需要实空间拓扑表征的理论计算

参考文献

本文引用了98篇参考文献，涵盖了拓扑绝缘体、实空间拓扑标记、非晶材料、高阶拓扑相等多个重要研究领域的核心文献，为该研究提供了坚实的理论基础。