四分之一磁通 Harper-Hofstadter 模型的精确解析能谱、本征态与量子几何

Exact analytical spectrum, eigenstates, and quantum geometry of the quarter-flux Harper-Hofstadter model

arXiv: 2608.26099v1

论文信息

标题: Exact analytical spectrum, eigenstates, and quantum geometry of the quarter-flux Harper-Hofstadter model

作者: Isaac Tesfaye, André Eckardt

发布日期: 2026-08-26

arXiv ID: 2608.26099v1

PDF 链接: 下载 PDF

3 分钟速览

  • 研究问题：该文要给出四分之一磁通 Harper-Hofstadter 模型的精确解析能谱、本征态和全部四带的量子几何，并评估最低陈带的分数陈绝缘体稳定性。
  • 核心方法：利用对称 2×22\times2 磁元胞中幸存的子晶格对称性，将 4×44\times4 布洛赫哈密顿量化为反块对角形式，再对剩余 2×22\times2 块做奇异值分解；同时推导任意子晶格对称系统的量子几何张量分解公式。
  • 关键结果：得到闭式能谱（式 (5)）和闭式贝里曲率（式 (8)）；最低带陈数 C1=−1C_1=-1，其 trace defect 积分只比拓扑界 2π∣C1∣2\pi|C_1| 高约 8.3%，在积分意义上接近理想陈带。
  • 主要局限：解析可解性依赖 q=4q=4 且使用对称磁元胞；对 Harper-Hofstadter 族，子晶格对称性只在 qq 可被 4 整除时于动量空间存活，q=8q=8 已需处理 4×44\times4 块，无法保证仍然得到闭式解。
  • 适合读者：研究拓扑能带、量子几何、分数陈绝缘体、超冷原子或光子量子模拟的凝聚态物理与量子模拟方向研究者。

论文背景和研究动机

量子几何张量（QGT）已经成为原子物理和凝聚态物理中描述布洛赫能带拓扑响应的统一语言。QGT 的虚部给出贝里曲率，实部给出量子度量，两者共同影响量子霍尔电导、平带超流权重、轨道磁性、非线性光学响应以及瓦尼尔函数局域化等物理量。

Harper-Hofstadter 模型描述均匀磁场下二维正方晶格上的带电粒子，是拓扑能带绝缘体的教科书模型，并在超冷原子、光子系统和超导电路中实现了量子模拟。尽管理论上对于有理磁通 α=p/q\alpha=p/q 可化为 q×qq\times q 布洛赫哈密顿量，但在本文之前，除 q=3q=3 的三带情形外，q>3q>3 情形一直没有闭合形式的能谱和量子几何解，通常只能依赖数值或微扰方法。该文针对 q=4q=4 的四分之一磁通情形，给出了全部四带能谱、本征态和量子几何的精确解析表达式。

核心方法和技术细节

该文的关键观察是：在所有四分之一磁通磁元胞选择中，只有对称的 2×22\times2 磁元胞能在动量空间中保留子晶格对称性。这里需要区分手征对称与子晶格对称：手征对称只要求幺正对合算符与哈密顿量反对易；子晶格对称则更强，要求该算符在格点基下对角，从而把晶格分成奇偶两个子晶格。真实空间中的棋盘算符 S^=(−1)m+n\hat S=(-1)^{m+n} 同时具有这两种性质，但只有对称磁元胞的磁平移算符与 S^\hat S 对易，因此 S^\hat S 才能作为与 k\bm{k} 无关的 grading 进入动量空间。

在合适排序后，布洛赫哈密顿量变为

HS(k)=(0DkDk†0),H_S(\bm{k})=\begin{pmatrix}0&D_{\bm{k}}\\D_{\bm{k}}^\dagger&0\end{pmatrix},

其中 DkD_{\bm{k}} 是 2×22\times2 的块。由于 HS2H_S^2 块对角化，能带由 DkD_{\bm{k}} 的奇异值的平方根给出。定义 Mk=Dk†Dk=h0σ0+h⋅σM_{\bm{k}}=D_{\bm{k}}^\dagger D_{\bm{k}}=h_0\sigma_0+\bm{h}\cdot\bm{\sigma}，其本征值为 h0±∣h∣h_0\pm|\bm{h}|。由此得到四个能带

E±,j(k)/J=±4+ηj26+cos⁡(4kxa)+cos⁡(4kya),E_{\pm,j}(\bm{k})/J =\pm\sqrt{4+\eta_j\sqrt{2}\sqrt{6+\cos(4k_x a)+\cos(4k_y a)}},

其中 ηj=(−1)j+1\eta_j=(-1)^{j+1}。这对应论文式 (5)。最低带与最高带陈数均为 −1-1，中间两条带在若干高对称点接触，形成总陈数 C2+3=2C_{2+3}=2 的 super-band。

本征态由 DkD_{\bm{k}} 的左右奇异向量构造：

∣Ψϵ,j⟩=12(∣uj⟩ϵ∣vj⟩),|\Psi_{\epsilon,j}\rangle =\frac{1}{\sqrt{2}} \begin{pmatrix}|u_j\rangle\\\epsilon|v_j\rangle\end{pmatrix},

其中 ∣vj⟩|v_j\rangle 是 MkM_{\bm{k}} 的本征态，∣uj⟩=Dk∣vj⟩/E+,j|u_j\rangle=D_{\bm{k}}|v_j\rangle/E_{+,j}。该奇异值分解关系固定了 ∣uj⟩|u_j\rangle 与 ∣vj⟩|v_j\rangle 的相对相位，这是得到正确本征态的关键。

在量子几何方面，该文导出了一个普适结论：对任意具有子晶格对称性的系统，阿贝尔 QGT 可以分解为两个子晶格扇区的贡献。体现在式 (6) 和式 (7) 中：

Qμνj=12(Qμν(uj)+Qμν(vj))+14ΔAμjΔAνj,Q^j_{\mu\nu} =\frac12\big(Q^{(u_j)}_{\mu\nu}+Q^{(v_j)}_{\mu\nu}\big) +\frac14\Delta A^j_\mu\Delta A^j_\nu,

其中 Q(uj)Q^{(u_j)}、Q(vj)Q^{(v_j)} 是两个扇区的 QGT，ΔAμj\Delta A^j_\mu 是贝里联络差。进一步拆分后，贝里曲率只是两个扇区的平均，而联络差只进入量子度量。因此手征对称伙伴带具有完全相同的贝里曲率和陈数。

创新点和贡献

该工作的主要贡献可以概括为三点。

第一，识别出对称 2×22\times2 磁元胞是保留子晶格 symmetry 的唯一选择，从而把四带问题约化为一个 2×22\times2 块，首次给出四分之一磁通 Harper-Hofstadter 模型的闭式能谱和本征态。论文还强调，三种磁元胞描述的是同一个 Landau 规范下的同一哈密顿量，不是不同规范选择。

第二，给出任意子晶格对称系统的阿贝尔与非阿贝尔 QGT 分解公式。这一普适公式不限于 Harper-Hofstadter 模型，只要系统具有子晶格 grading，就可以把完整量子几何分成物理子晶格扇区贡献。论文还为简并中间带导出了非阿贝尔 QGT 的对应分解，并验证了 tr Fxy=−2Ωxy1{\rm tr}\,F_{xy}=-2\Omega_{xy}^1，与所有带的总贝里曲率守恒一致。

第三，利用闭式 QGT 定量评估分数陈绝缘体稳定性。最低带的贝里曲率为

Ωxy1(k)=−2a26−cos⁡(4kxa)−cos⁡(4kya)[6+cos⁡(4kxa)+cos⁡(4kya)]3/2,\Omega^1_{xy}(\bm{k}) =-\sqrt{2}a^2 \frac{6-\cos(4k_x a)-\cos(4k_y a)} {[6+\cos(4k_x a)+\cos(4k_y a)]^{3/2}},

该式即论文式 (8)。最低带在整个布里渊区都保持负曲率，因此陈数明确为 −1-1。

解析结果分析

论文利用闭式贝里曲率和量子度量评估最低带的近理想性。二维情形下，QGT 的正半定性给出局域不等式

tr g(k)≥2det⁡g(k)≥∣Ωxy(k)∣.{\rm tr}\,g(\bm{k})\ge 2\sqrt{\det g(\bm{k})}\ge|\Omega_{xy}(\bm{k})|.

最低带的 trace defect

D1(k)=tr g1(k)−∣Ωxy1(k)∣D^1(\bm{k})={\rm tr}\,g^1(\bm{k})-|\Omega^1_{xy}(\bm{k})|

在布里渊区内点状非负，并在高对称点消失（图 2）。其积分给出全局 trace condition

THH=∫BZD1(k) d2k.\mathcal{T}_{\mathrm{HH}} =\int_{\mathrm{BZ}}D^1(\bm{k})\,d^2k.

论文计算得到 THH\mathcal{T}_{\mathrm{HH}} 仅比拓扑界 2π∣C1∣2\pi|C_1| 高约 8.3%，因此最低带在积分意义上是 “近理想” 陈带。论文还给出量子权重

K≡Kxx+Kyy=∣C1∣+THH/2π≈1.083,K\equiv K_{xx}+K_{yy} =|C_1|+\mathcal{T}_{\mathrm{HH}}/2\pi \approx 1.083,

超过拓扑下界 K≥∣C1∣=1K\ge|C_1|=1，这与静态结构因子的小 q\bm{q} 行为相关。

此外，从精确能谱得到最低带的能隙与带宽分别为

W/J=22−4+22,Δ/J=22−2,W/J=2\sqrt{2}-\sqrt{4+2\sqrt{2}}, \qquad \Delta/J=2\sqrt{2-\sqrt{2}},

因此精确的平坦度比值为 f=Δ/W≈7.1097f=\Delta/W\approx 7.1097。论文比较指出，在没有精细调节的情况下，这比 Haldane 模型在最优精细调节点的平坦度 f≈6.0136f\approx 6.0136 更平。不过论文也强调，该带并非几何平坦，其贝里曲率和量子度量在布里渊区内仍会波动，理想性只在积分意义下接近，而非逐点饱和。

局限与待解决问题

该方法的核心局限来自对称磁元胞条件的限制。在 Harper-Hofstadter 族中，子晶格对称性只在 qq 可被 4 整除时才在动量空间幸存。对于 q=4q=4，剩余块是 2×22\times2，可以解析对角化；但对于 q=8q=8，剩余块已经变成 4×44\times4，无法保证继续得到闭式解。因此本文的方法不能直接推广到一般有理磁通。

更一般地，虽然 QGT 分解公式适用于任意子晶格对称系统，但能否得到闭式结果仍取决于两个子晶格扇区是否各自可解。论文建议可将该方法用于 BBH 四极绝缘体、方八边形晶格、dice 晶格和 Lieb 晶格等模型；对于具有多数子晶格的系统，零模和平坦简并扇区需要非阿贝尔 QGT 来处理。对于更大的块矩阵，论文认为可以结合本征投影方法，但尚未在本工作中展开。

另一个值得注意的限制是，手征对称只保证半满时能谱关于零能对称，此时零能处无能隙，没有缠绕数拓扑不变量。若要讨论最低带的陈数，需要引入化学势或孤立占据最低带，这同时也会破坏手征对称，使模型进入 class A。因此最低带的拓扑性质与手征对称本身并不是同一框架下的直接结果，而是通过能隙隔离得到的陈数。论文未讨论相互作用和有限温度下这些解析量子几何量的稳定性，也未给出更大 qq 情形的近似解析方案。