晶体绝缘体中的高阶 Berry 曲率、第二陈数与磁电耦合

论文信息

标题: Higher Berry curvature, second Chern numbers and magnetoelectric coupling in crystalline insulators

作者: Niclas Heinsdorf, Ken Shiozaki

发布日期: 2026-06-24

arXiv ID: 2606.26096v1

PDF 链接: 下载 PDF

3 分钟速览

• 研究问题：如何用一维量子多体态的高阶 Berry 曲率直接计算四维 Chern 绝缘体的第二 Chern 数，并探讨高阶 Berry 相位与三维轨道磁电耦合（Chern–Simons 轴子角）之间的关系。
• 核心方法：将四维晶格 Dirac 模型重构为三维参数空间上的一维平移不变链族，利用无限矩阵乘积态（iMPS）计算离散的高阶三形式 Berry 曲率，得到显式量子化的 Dixmier–Douady–Kapustin–Spodyneiko（DDKS）数。
• 关键结果：DDKS 数 $\nu$ 与第二 Chern 数 $ch_2$ 在完整相图中完全一致（$\nu=ch_2$），且高阶 Berry 相位 $\theta_{\text{HBP}}$ 在 Dirac 模型中精确再现了 Chern–Simons 轴子角的变化 $\Delta\theta_{\text{CS}}$，但在 Hopf 绝缘体中二者明显偏离。
• 主要局限：算法目前仅适用于一维平移不变系统，且 iMPS 计算所需的键维数依赖于能隙；Hopf 绝缘体中 $\theta_{\text{HBP}}$ 与 $\theta_{\text{CS}}$ 的关系尚无理论解释，高阶 Berry 曲率能否作为通用磁电响应标识仍不明确。
• 适合读者：从事拓扑物态、张量网络算法、量子输运与量子模拟研究的理论物理学者，以及对多体拓扑不变量及其数值计算感兴趣的研究生。

论文背景和研究动机

拓扑不变量为量子多体基态的分类提供了强大框架。在零维参数族中，Berry 曲率及其积分（Chern 数）是标准工具；近年来，这一视角被推广到具有空间扩展性的短程纠缠态族，引出了高阶微分形式和更高阶的拓扑不变量。Kitaev 猜想，短程纠缠可逆态的拓扑分类可以推广到任意空间维度的基态族，其中 “维度” 是量子态的物理维度而非参数空间维度。在 Haldane 模型中，Chern 数是基于动量空间作为参数空间而定义，本质上处理的是一族零维系统。Kapustin 和 Spodyneiko 提出了高阶 Berry 曲率的框架：对一维系统，Berry 曲率由二形式变为三形式；对二维系统，则变为四形式等等。相应的高阶不变量被称为 DDKS 数，可以由高阶 Berry 曲率在参数空间上积分得到。

然而，直到最近，高阶 Berry 曲率的显式计算主要局限于自由费米子体系，或可约化为二体问题的自旋模型。Shiozaki 等人提出了一种基于 iMPS 的算法，通过离散化参数空间，将高阶 Berry 曲率表达为三角剖分上相位的和，从而给出显式量子化的 DDKS 数。该算法类似于计算一阶 Chern 数的 Fukui-Hatsugai-Suzuki 方法，但处理的是一维系统基态的 gerbe 结构。

本文作者将这一算法应用于四维 Chern 绝缘体。四维 Chern 绝缘体的能带拓扑由第二 Chern 数 $ch_2$ 表征，其相图已知。作者将四维晶格哈密顿量通过部分 Fourier 变换重新解释为三维参数空间上的一维链族，利用 iMPS 求得基态并计算高阶 Berry 曲率。这种对应关系有一个数学基础：利用 Künneth 公式，$H^4(S^1\times T^3,\mathbb{Z})\cong H^3(T^3,\mathbb{Z})$，因此 DDKS 数 $\nu$ 应与 $ch_2$ 相同。文章通过数值计算验证了这一等价性，并进一步探讨了高阶 Berry 相位与轨道磁电耦合中 Chern–Simons 轴子角 $\theta_{\text{CS}}$ 的联系。

核心方法和技术细节

作者从 $(4+1)$ 维 Dirac 哈密顿量的晶格正则化出发：

$$H_{\mathbf{k}}=\psi^\dagger_{\mathbf{k}}\left[\sum_\mu \sin k_\mu \Gamma^\mu + \left(m+\sum_\mu \cos k_\mu\right)\Gamma^0\right]\psi_{\mathbf{k}},$$

其中 $\mathbf{k}\in T^4$，$m$ 为质量参数。利用 $T^4=S^1\times T^3$ 的分解，将 $k_4$ 方向的 Fourier 变换保留在实空间，得到一维链族 $H_{\boldsymbol\lambda}$，$\boldsymbol\lambda=(\lambda_1,\lambda_2,\lambda_3)\in T^3$。每个 $\boldsymbol\lambda$ 对应一个平移不变的一维双轨道链。通过 Jordan-Wigner 变换和矩阵乘积算符表示，使用 iDMRG 获得基态 iMPS，其张量 $A^i_{\boldsymbol\lambda}$ 为 $d\times D\times D$ 矩阵（$d$ 为物理维度，$D$ 为键维数）。

在三维参数空间上引入 $(N+1)^3$ 网格并进行三角剖分，每个小立方体划分为六个四面体 $\Delta_3$。对于由三个顶点构成的定向三角形 $\Delta_2$，定义相位：

$$\phi(\Delta_2)=\mathrm{Arg}\,\mathrm{Tr}[\Lambda_{\lambda_0}^{2/3} V_{\lambda_0\lambda_1} \Lambda_{\lambda_1}^{2/3} V_{\lambda_1\lambda_2} \Lambda_{\lambda_2}^{2/3} V_{\lambda_2\lambda_0}],$$

其中 $\Lambda$ 为 iMPS 的 Schmidt 矩阵，$V_{01}$ 为混合转移矩阵的主特征向量（过渡函数）。这个相位在 MPS 规范变换（投射酉群）下不变，但在过渡函数的标度变换下会改变，这正是高阶 Berry 曲率作为三形式的要求。四面体 $\Delta_3$ 的高阶 Berry 曲率则定义为其四个面的相位和：

$$F(\Delta_3)=\phi(\partial\Delta_3).$$

将所有四面体的 $F(\Delta_3)$ 相加，除以 $2\pi$，即可得到显式整数的 DDKS 数 $\nu$。这一构造从几何上类似于将三维参数空间上的 gerbe 曲率积分，数值上严格避免了规范固定问题，类似于离散格点上的 Chern 数算法。

同时，作者用另一种不显式量子化的方法计算了同一模型的第二 Chern 数 $ch_2$，以做对比。对 Bloch 哈密顿量在四维动量网格（$10^4$ 点）上对角化，用链接变量构建场强并积分。在足够大的网格上，$ch_2$ 的数值趋于整数，但有限网格下存在误差。

创新点和贡献

1. 将高阶 Berry 曲率算法首次应用于四维拓扑绝缘体：尽管该算法已在一维自旋模型中展示，但将其推广到自由费米子体系并验证其等价于第二 Chern 数具有重要意义，为未来处理相互作用一维链族提供了基准。

2. 提供了 DDKS 数与第二 Chern 数完全一致的数值证据：相图中，$\nu$ 随着 $m$ 的变化与解析公式给出的 $ch_2$ 完全一致（图 4(a)），且在较粗糙的 $4^3$ 网格上 $\nu$ 就已给出准确整数，而 $ch_2$ 的常规计算需要 $10^4$ 网格才能使误差小于 0.5。这种显式量子化特性正是因为算法在构造上消除了规范不确定性。

3. 探讨了高阶 Berry 相位与磁电耦合的关系：提出了一个自然的问题：既然第二 Chern 数控制着 Chern–Simons 有效作用量，而 DDKS 数等于 $ch_2$，那么三维切片上的高阶 Berry 相位 $\theta_{\text{HBP}}$ 是否直接等于轴子角 $\theta_{\text{CS}}$？对于 Dirac 模型，当固定 $\lambda_3$ 时，积分 $F$ 得到的 $\theta_{\text{HBP}}$ 与 $\Delta\theta_{\text{CS}}$ 高度吻合（图 5(c)），且在高对称动量点被量子化。但在 Hopf 绝缘体中，$\theta_{\text{HBP}}$ 不与 $\theta_{\text{CS}}$ 一致（图 6(c)），揭示了高阶 Berry 相位包含的信息可能比轴子角更丰富，或对 delicate topology 的敏感度不同。

4. 计算了 Hopf 绝缘体的高阶 Berry 相位：Hop 绝缘体是一个两带模型，其拓扑由 Hopf 不变量 $\chi$（等于 $\theta_{\text{CS}}/\pi$）表征。作者构造了一维链族并用 iMPS 计算 $\theta_{\text{HBP}}$，发现对于 $|m|<3$，$\theta_{\text{HBP}}$ 既不量子化也不等于 $\pi\chi$，从而首次给出了 $\theta_{\text{HBP}}$ 与 $\theta_{\text{CS}}$ 不一致的反例。这暗示 delicate 拓扑的保护可能在 iMPS 描述下被破坏，或高阶 Berry 曲率探测的是不同类型的几何信息。

实验结果分析（高阶 Berry 曲率与相图）

在 Dirac 模型中，作者选取 $m=1$ 和 $m=3$ 分别对应于 $ch_2=3$ 和 $-1$ 的拓扑相。图 3 展示了参数空间上高阶 Berry 曲率 $2\pi F$ 的分布：其绝对值在能隙最小处最大，这与预期的 Berry 曲率集中在带反转点附近相符。求和后得到的 $\nu$ 分别为 $3$ 和 $-1$，与 $ch_2$ 完全匹配。相图扫描（图 4(a)）显示 $\nu$ 完美复现 $ch_2$ 的台阶结构。

对 $m=1$，收敛性测试（图 4(b)）表明，即使只用 $4^3$ 的粗糙网格，$\nu$ 就已经得到正确的整数 3，而常规 $ch_2$ 计算在 $10^4$ 网格仍有约 0.06 的偏差。这说明 iMPS 算法的显式量子化对网格尺寸不敏感。

在 Dirac 模型中进一步检验 $\theta_{\text{HBP}}$ 与 $\Delta\theta_{\text{CS}}$ 的关系。取不同 $\lambda_3$ 切片，计算 $\theta_{\text{HBP}}$ 并与解析结果比较（图 5(c)）。在 $\lambda_3$ 远离 TRIM 时，$\Delta\theta_{\text{CS}}$ 无对称性强制量子化，而 $\theta_{\text{HBP}}$ 的数值几乎落在 $\Delta\theta_{\text{CS}}$ 曲线上，展示出二者的直接对应。在 Hopf 绝缘体中，情况截然不同（图 6(c)）：$\theta_{\text{HBP}}/2\pi$ 在 $|m|<3$ 区间呈现非量子化的平台，但远远低于 $\theta_{\text{CS}}/\pi$ 值（如 $m=0.5$ 时 $\chi=1$，因此 $\theta_{\text{CS}}/\pi=1$，而 $\theta_{\text{HBP}}/2\pi\approx 0.2$），说明 $\theta_{\text{HBP}}$ 不能直接等同于轴子角。论文指出，这可能因为 Hopf 不变量是 delicate 拓扑，仅在严格的两带模型中定义，而 iMPS 表示已隐式地超越了单粒子图像。两种方法对拓扑的 “感知” 不同，深层原因有待阐明。

局限与待解决问题

本文的数值工作证实了 iMPS 算法在自由费米子模型中有效，但该方法的应用边界和物理诠释仍存在若干未解问题。

一维限制：当前高阶 Berry 曲率构造仅适用于一维系统（iMPS）。更高空间维度需要推广到 PEPS 等张量网络，相关算法虽有初步发展（如 Sommer 等 2024 年的工作），但尚未在类似模型上验证。

Hopf 模型中的不一致：$\theta_{\text{HBP}}$ 与 $\theta_{\text{CS}}$ 在 Hopf 绝缘体中的差异需要理论解释。该模型具有 delicate topology，其拓扑性质在超越两带近似时可能无法稳定。iMPS 计算中使用的键维数是否足以保护 delicate 特征，以及高阶 Berry 相位是否本质上探测的是不同于 Chern–Simons 耦合的响应量，都是有待厘清的问题。论文并未给出解决线索。

相互作用的推广：尽管 iMPS 天然适合处理强关联系统，作者只是展望了在相互作用下的应用，如 Hubbard 型相互作用可能导致质量重整化，但是否能驱动高阶拓扑相变仍未知。他们明确指出，单纯在一维链上加入局域相互作用并不等于四维相互作用 Chern 绝缘体，因为关联仅沿实空间一维方向传播，而不是跨参数空间，对应的实验平台需要准一维系统阵列，如有机超导体或 Li₀.₉Mo₆O₁₇ 等。

磁电响应的直接关联：目前尚无已知公式将高阶 Berry 曲率直接与三维输运系数对应。虽然 $\theta_{\text{HBP}}$ 在 Dirac 模型中再现了 $\Delta\theta_{\text{CS}}$，但缺乏普适证明。高阶 Berry 相位可能提供量子几何的新诊断量，但究竟对应于哪一类物理响应仍不明确。

算法的数值效率：虽然 iMPS 方法在粗糙网格上就能给出量子化整数，但 iMPS 计算本身需要对角化 $D^2\times D^2$ 的混合转移矩阵，且 iDMRG 收敛成本随 $D^3$ 增长。对于有复杂纠缠结构的系统，所需 $D$ 可能很大。论文中 $D=64$ 已足够，但更接近相变点或存在长程关联时，计算成本将上升，限制了直接扫描整个相图的可行性。

总体而言，本文为高阶 Berry 曲率在拓扑物态中的应用提供了重要范例，但也揭示了 delicate topology 或强关联情形下需要更细致的理论框架。这些开放问题将成为后续研究的焦点。