从一维到二维的困境

非厄米准周期修饰Lieb晶格中的精确迁移率环

作者

快速导读：本文通过将二维准周期修饰Lieb晶格精确映射到广义AAH模型，解析推导出复能量平面中精确的迁移率环边界，并揭示其随修饰参数κ的拓扑演化行为。

非厄米准周期系统中，迁移率边（ME）这一区分扩展态与局域态的临界边界，在复能量平面中推广为闭合的迁移率环（MR）。一维非厄米准周期模型中的MR已有广泛研究，但二维精确可解模型长期处于空白状态。本文的核心贡献在于：利用二维准周期修饰Lieb晶格（QDLL）的多子晶格对称性，通过代数消元将二维紧束缚哈密顿量精确映射到有效二维广义AAH模型，从而解析推导出复能量平面中精确的MR闭合表达式。

这一映射的关键在于QDLL仅在顶点位置施加复准周期势，使得非顶点子晶格的自由度可以被代数消去。消元后得到的特征函数a_κ(E)以第二类Chebyshev多项式的形式出现，其模值条件|a_κ(E)V e^{iθ}|=1直接给出MR边界。该表达式揭示了一个深刻的物理事实：MR边界对非厄米相位角θ不变，仅由准周期势强度V和修饰参数κ决定。

英文题目：Exact mobility rings in non-Hermitian quasiperiodically decorated Lieb lattices

论文出处：arXiv 每日论文精选 · arXiv:2608.16719

原始论文：PDF / 论文页面

这篇论文解决什么问题？

迁移率边是安德森局域化研究中的核心概念。在厄米系统中，ME是实能量轴上分离扩展态与局域态的临界能量；在非厄米系统中，由于本征能量变为复数，实轴上的ME推广为复能量平面中的闭合迁移率环。Li & Li在Phys. Rev. B 110, L041102 (2024)中首次提出MR作为非厄米迁移率边的指纹概念，此后一维非厄米准周期模型中的MR研究迅速发展。

然而，二维非厄米系统中精确可解的MR解析框架极为稀缺。根本困难在于：二维系统的转移矩阵计算涉及高维矩阵乘积，存在严重的有限尺寸效应和收敛瓶颈。即便在数值上，要精确刻画复平面中的闭合临界边界也极具挑战。因此，寻找一个二维模型，使其MR边界可以被解析推导，成为该领域的关键开放问题。

本文选择的平台是准周期修饰Lieb晶格（QDLL）。Lieb晶格因其平带结构和多子晶格几何而备受关注，Wang et al.在Phys. Rev. B 107, L140201 (2023)中已在二维顶点修饰Lieb晶格中发现了厄米情形下的精确迁移率边。本文将该模型推广到非厄米情形，在顶点位置施加复准周期势，从而打开通往精确MR解析解的道路。

核心创新

  • 首次在二维非厄米准周期系统中解析推导出精确的迁移率环闭合表达式
  • 通过代数消元将二维多子晶格QDLL精确映射到有效二维广义AAH模型，规避了高维转移矩阵的收敛瓶颈
  • 揭示MR几何的拓扑演化：κ=2为单一圆环，κ=3在V_c=1处经Bernoulli双纽线分裂为两个独立子环
  • 证明MR边界对非厄米相位角θ不变，仅由准周期势强度V和修饰参数κ决定

方法概览

利用二维准周期修饰Lieb晶格（QDLL）的多子晶格对称性，通过代数消元将二维紧束缚哈密顿量精确映射到有效二维广义AAH模型；基于对偶变换和Avila全局理论计算Lyapunov指数，得到MR的闭合解析表达式|a_κ(E)V e^{iθ}|=1。

  • 模型构造：QDLL是扩展Lieb晶格，每个原胞包含五个子晶格（A至E），复准周期势V e^{iθ}仅施加在顶点位置（A子晶格），其余子晶格保持干净。这一选择性修饰是后续代数消元的关键前提。
  • 代数消元映射：利用非顶点子晶格（B至E）的紧束缚方程，将其自由度用顶点子晶格A的波函数代数表达，代入A子晶格方程后得到仅含A自由度的有效二维广义AAH模型。该映射在x和y方向均成立，如图1(d)所示。
  • 特征函数a_κ(E)：消元过程中产生的有效在位能项以第二类Chebyshev多项式的形式出现，记为a_κ(E)。对于κ=2，a_2(E)=E²−1；对于κ=3，a_3(E)=E³−2E。这一多项式结构直接决定了MR的几何形态。
  • 对偶变换与Avila全局理论：对有效AAH模型施加对偶变换，利用Avila全局理论计算Lyapunov指数（LE），得到MR的闭合解析表达式|a_κ(E)V e^{iθ}|=1。该表达式将复能量E、准周期势强度V和非厄米相位角θ统一在一个简洁的模值条件中。
  • MR几何的解析刻画：κ=2时MR方程为Re(E)²+Im(E)²=1/V²，是复平面中的单一圆环；κ=3时MR方程为(Re(E)²+Im(E)²)²−2(Re(E)²−Im(E)²)=1/V²−1，是Cassini卵形线。
  • 拓扑相变条件：κ=3时，当V_c=1时Cassini卵形线退化为Bernoulli双纽线，MR从单一闭合环分裂为两个独立子环。这一拓扑相变完全由解析表达式预测，无需数值拟合。

逐图理解论文

从一维到二维的困境

从一维到二维的困境
FIG. 1. Schematic illustration of the MEs, MRs and the 2D QDLL model. (a) Conceptual phase diagram in the complex energy plane. The solid red curve denotes the MR, which forms a closed boundary perfectly separating the extended states (interior) from the localized states (exterior). The dashed red lines illustrate the conventional MEs projected onto the real energy axis. (b) and (c) Structural representa- tions of the standard 2D Lieb lattice and the extended Lieb lattice, comprising three (labeled A, B, and C) and five (la- beled A through E) sub-lattices per unit cell, respectively. The complex quasiperiodic potentials are exclusively intro- duced at the vertex sites, highlighted as red spheres. (d) Schematic of the mapping from the 2D QDLL model onto an effective AAH model. Here we show the mapping in the x direction, which is similar in the y direction.

图1(a)展示复能量平面中的概念相图，红色闭合曲线为MR，内部为扩展态，外部为局域态；图1(b)-(c)展示标准Lieb晶格与扩展Lieb晶格的结构，红色球体标记施加复准周期势的顶点位置；图1(d)示意从二维QDLL到有效AAH模型的映射过程，是理解全文方法的核心图示。

拓扑演化：从圆环到双纽线分裂

拓扑演化：从圆环到双纽线分裂
FIG. 2. Energy spectrum, mobility edges (MEs), and mobility rings (MRs) for the κ = 2 case. (a) Wavefunction fractal dimension Γ as a function of the quasiperiodic potential strength V and real eigenvalues E in the Hermitian limit (θ = 0). The red dashed lines denote the critical ME boundary |Ec| = 1/V , which separating from extended to localized states. (b), (c) Spectral distributions for V = 0.8 and V = 1.5 at θ = 0. (d) Complex energy spectrum and exact MR boundary |E| = 1/V (solid red curve) for the non-Hermitian case θ = π/2. (e), (f) Spectral phase diagrams in the complex energy plane for V = 0.8 and V = 1.5 at θ = π/2. The eigenstates within the MR are extended (Γ ≈1), while the eigenstates outside the MR are localized (Γ ≈0). Here, t = 1 is set as the unit of energy, and the system size are Lx = Ly = 55 under periodic boundary conditions.

图2展示κ=2情形的完整验证：图(a)为厄米极限下分形维数Γ随V和E的分布，红色虚线为解析ME边界；图(d)为非厄米情形下的复能量谱与解析MR边界（红色实线圆环），谱点严格落在边界内外两侧，直观验证了解析表达式的准确性。

相位角不变性

相位角不变性
FIG. 5. The MRs for different θ for κ = 3 case. (a)-(b) represent the θ = π/4 and θ = 3π/4. Although the complex spectral features shift markedly with varying θ, the analytical sub-ring boundaries defined by Eq. (12) remain invariant.

图5展示κ=3时不同非厄米相位角θ=π/4和θ=3π/4下的复能量谱：尽管谱分布随θ显著移动，解析MR边界（红色曲线）保持不变，直接验证了MR对θ的不变性这一核心理论预言。

实验如何设计？

  • 数值验证采用分形维数（FD）作为态分类指标，通过逆参与比（IPR）定义Γ=−lim ln(IPR)/ln N，Γ≈1对应扩展态，Γ≈0对应局域态。
  • 厄米极限验证：在θ=0时计算波函数分形维数Γ随V和实本征能量E的分布，验证ME边界|E_c|=1/V（κ=2）和E_c=±√(1±1/V)（κ=3）。
  • 非厄米验证：在θ=π/2时计算复能量谱与分形维数，验证MR边界与解析表达式的吻合。系统尺寸为Lx=Ly=55，采用周期性边界条件。
  • 相位角不变性检验：对κ=3情形分别计算θ=π/4和θ=3π/4的复能量谱，验证尽管谱分布随θ显著移动，MR边界保持不变。
  • 典型本征态空间分布：选取MR内外的典型复本征值，计算实空间概率分布|ψ|²，直观展示扩展态与局域态的空间特征差异。

关键结果与论文证据

  • κ=2时MR为单一圆环，方程为Re(E)²+Im(E)²=1/V²，半径由1/V决定（第4页）。
  • κ=3时MR为Cassini卵形线，方程为(Re(E)²+Im(E)²)²−2(Re(E)²−Im(E)²)=1/V²−1（第5页）。
  • κ=3时在V_c=1处MR经Bernoulli双纽线发生拓扑分裂，从单一闭合环变为两个独立子环（第5页）。
  • MR边界对非厄米相位角θ不变，仅由V和κ决定；θ的变化仅移动谱分布，不改变临界边界（第5页）。
  • 厄米极限下κ=2的ME为E_c=±1/V，κ=3的ME为E_c=±√(1±1/V)，与MR在θ→0时的投影一致（第4页）。
  • 数值分形维数Γ在MR内外呈现清晰的二值分布（Γ≈1与Γ≈0），验证了解析边界的准确性（第4-6页）。

阅读时需要注意

  • 数值验证受限于有限系统尺寸（Lx=Ly=55或34），未做热力学极限外推，边界附近的有限尺寸效应未被系统评估。
  • 仅考虑周期性边界条件，未讨论开放边界条件或非厄米趋肤效应对MR的影响。
  • 模型仅针对特定修饰参数κ=2和κ=3，未系统讨论更大κ值的普适性。
  • 文中图3(a)-(b)的图注与正文描述存在不一致：图注标注κ=2和κ=3，正文描述为θ=π/4和3π/4，需注意统一。

关联工作

  • Wang et al., Phys. Rev. B 107, L140201 (2023)：二维顶点修饰Lieb晶格中精确迁移率边与鲁棒平带，是本文的厄米前驱工作。
  • Jiang et al., New J. Phys. 26, 083020 (2024)：二维Lieb晶格中虚准周期势的精确非厄米迁移率边与平带，与本文共享Lieb晶格平台但修饰方式不同。
  • Li & Li, Phys. Rev. B 110, L041102 (2024)：提出复平面环结构作为非厄米迁移率边的指纹，奠定MR概念基础。
  • Pang et al., Phys. Rev. B 111, 214205 (2025)：平带晶格中非厄米准周期势的精确迁移率线与迁移率环，与本文同期探索类似问题。

发表评论