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

论文信息

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

作者: Niclas Heinsdorf, Ken Shiozaki

发布日期: 2026-06-24

PDF 链接: 下载 PDF

引言

在凝聚态物理中，拓扑不变量为理解量子多体基态的分类提供了强大的工具。经典的例子是通过贝里曲率（Berry curvature）来刻画参数化家族中有限维度基态的拓扑性质。近年来，这一视角从零维系统延伸到了空间扩展的短程纠缠态，引入了更高阶的几何结构和微分形式。Kitaev 曾猜想，具有无穷空间延展的短程纠缠可逆态的分类可以推广到任意维度，这催生了对高维拓扑不变量的系统性探索。

Kapustin 和 Spodyneiko 提出了一个高阶贝里曲率框架：对于一维空间系统，通常的二维贝里曲率被三维形式所取代；对二维系统则是四维形式，以此类推。对这些高阶曲率进行积分，得到的整数不变量被称为 Dixmier–Douady–Kapustin–Spodyneiko (DDKS) 数。然而，早期的数值计算局限于自由费米子或可简化为两格点问题的自旋模型。最近，基于无限矩阵乘积态（iMPS）的算法被提出，使得一般平移不变一维哈密顿量的高阶贝里曲率计算成为可能。本文正是在这一背景下，首次用 iMPS 计算了一族四维陈绝缘体降维后的一维链的高阶贝里曲率，展示了其二阶陈数与 DDKS 数的精确对应，并深入探讨了高阶贝里相与三维轨道磁电耦合之间的关系。

高阶贝里曲率与 DDKS 数

高阶贝里曲率的核心思想源于参数空间中量子态的几何结构。对于一族一维的系统，其基态波函数在参数空间上形成一个希尔伯特空间上的纤维丛，但由于规范冗余，普通的贝里联络和曲率并不能直接推广。Kapustin 和 Spodyneiko 发现，对于一维系统，相关的几何对象是一个三形式曲率，其积分给出一个三维上同调类，即 DDKS 数。从物理上看，一阶陈数可视为参数循环中泵浦的电荷数目，而 DDKS 数则度量了参数空间中泵浦的 “陈数” 数目，起到了更高阶拓扑泵浦的作用。

计算这种高阶曲率的关键在于正确地处理规范自由度。在矩阵乘积态（MPS）的框架下，基态由一系列张量描述，其规范的局部酉变换以及全局的相位不确定性都会影响相位的定义。作者采用了一种类似于零维 Fukui-Hatsugai-Suzuki 算法的策略：将三维参数空间三角剖分为四面体，在每个四面体上定义由过渡函数和施密特权重组合而成的离散高阶贝里曲率。通过取四面体边界上四个三角形贡献的代数和，过渡函数的标度不确定性恰好抵消，从而得到了规范不变的曲率。对整个参数流形求和后，结果自然是一个整数，即 DDKS 数 $\nu$ ，对应 $H^3(T^3,\mathbb{Z})$ 中的一个类。

从四维陈绝缘体到一维参数化链

四维陈绝缘体的标准晶格哈密顿量为

H_k = \psi^\dagger_k \left[ \sum_\mu \sin k_\mu \Gamma^\mu + \left(m + \sum_\mu \cos k_\mu\right) \Gamma^0 \right] \psi_k,

其中 $k \in T^4$ ， $m$ 是质量参数。其拓扑由二阶陈数 $ch_2$ 刻画，解析表达式表明 $ch_2$ 在 $m$ 变化时取 $-3,-1,0,1,3$ 等整数值。为了将高阶贝里曲率的方法应用于此模型，作者利用布里渊区的乘积结构 $T^4 = S^1 \times T^3$ ，将第四维动量 $k_4$ 作傅里叶变换，得到一个平移不变的一维链哈密顿量，其参数由剩余三维动量 $\lambda = (k_1,k_2,k_3)$ 给出。这一族一维链保持间隙，除非 $m$ 达到临界值 $0,\pm2,\pm4$ 。在非临界区域，基态可以用 iMPS 近似表示，且由于短程纠缠的平凡性，可取为内射（injective）形式。这为后续在 $T^3$ 参数空间中计算高阶贝里曲率提供了理想场景。

数值方法与结果

在数值实现中，作者在三维参数空间中建立均匀网格，并在每个格点上通过 iDMRG 得到基态 iMPS。为了适应不同参数点间纠缠结构的差异，允许键维度 $D$ 自适应变化，文中使用了典型的 $D=64$ 。对于每两个相邻参数点，计算混合转移矩阵的主导本征向量 $V_{01}$ （即过渡函数），并结合施密特权重 $\Lambda$ 构造每个三角形面上的相位 $\phi(\Delta_2)$ 。然后按照三角剖分，将六个三角形围绕立方体的六个面组合为四面体的边界，求和得到每个小立方体内部的高阶贝里曲率 $F(\Delta_3)$ 。最终对全空间求和得到 DDKS 数 $\nu$ 。

结果显示， $\nu$ 随 $m$ 的变化与解析的二阶陈数完全一致，复制出 $ch_2 = -3,-1,0,1,3$ 的相图。特别地，在 $m=1$ 时， $\nu=3$ ，且只需 $4^3$ 网格就能得到精确整数，而传统基于四维布里渊区离散化的二阶陈数计算则需要 $10^4$ 以上网格才能达到误差小于 $0.5$ 。如图 3 所示，高阶贝里曲率集中在能隙较小的区域，这正是四维狄拉克模型在动量空间中 Berry 曲率热点的自然投影。这一结果证明，高阶贝里曲率算法本身就是量化的，无需外推，且能处理具有非平凡一阶陈数消失的情况，克服了传统方法中需要全局光滑规范的限制。

高阶贝里曲率与磁电耦合

高阶贝里曲率的物理意义不止于数学上的拓扑分类。论文进一步探讨了其与三维系统中轨道磁电耦合的联系。在三维绝缘体中，轨道磁电耦合的轴子项 $\theta_{CS}$ 由陈-西蒙斯三形式的积分给出，其二阶陈形式的关系为 $d\theta_{CS} \propto \text{Tr}[F\wedge F]$ 。若将参数空间的一个方向视为绝热泵浦参数，则在泵浦周期内 $\theta_{CS}$ 的变化等于二阶陈数的积分。类似地，DDKS 数可通过定义在每个固定参数二维切片上 “高阶贝里相” $\theta_{\text{HBP}}$ 的泵浦总和给出。因此，自然期望 $\theta_{\text{HBP}}$ 与 $\theta_{CS}$ 存在直接对应。

在四维狄拉克模型中，计算表明 $\Delta\theta_{CS}$ 和 $\Delta\theta_{\text{HBP}}$ 随泵浦参数 $\lambda_3$ 的变化完全吻合，即使在破缺时间反演对称性的非量子化区域也是如此。这给出了高阶贝里相与磁电轴子响应的一个直接实现。然而，这一关系并非普遍。当作者将同样的方法应用于 Hopf 绝缘体时，情况出现分歧。Hopf 绝缘体的陈-西蒙斯角度由 Hopf 不变量 $\chi$ 决定： $\theta_{CS} = \pi\chi \mod 2\pi$ ，且 $\chi$ 在几个区间内取 $0,1,-2$ 。但用 iMPS 计算出的 $\theta_{\text{HBP}}$ 既不等于 $\theta_{CS}$ ，也不量子化，甚至在远离相变时也显著偏离。这种差异很可能源于 Hopf 不变量仅对双带模型定义良好，而在 iMPS 描述中，一维链的热力学极限可能并未保留动量空间中精细的纤维丛结构。这揭示出高阶贝里相可能捕捉到不同于陈-西蒙斯响应的几何信息，值得进一步探索。

实际应用与未来展望

本文的工作为高阶拓扑不变量在量子多体系统中的数值研究开辟了新途径。iMPS 框架天然适用于处理强关联系统，而高阶贝里曲率的离散量化算法可以无缝对接。现实中，准一维材料（如有机超导体 Li $_{0.9}$ Mo $_6$ O $_{17}$ 等）以及量子模拟平台（冷原子、光子学中的合成维度）为检验 DDKS 数的物理效应提供了可能。例如，通过对一维链阵列施加周期性调制，可以实现等效的高维泵浦，从而测量高阶贝里相的特征。此外，高阶不变量还可用于分类含有多次简并的拓扑半金属的节点结构，最近已有研究将其应用于三次费米子等体系。

从方法学角度看，尽管此次工作集中在自由费米子模型以验证一致性，但 iMPS 方法真正发光发热的领域应当是具有相互作用的拓扑态。局部相互作用（如 Hubbard $U$ ）可以通过重整化质量项驱动拓扑相变，并可能重新分布参数空间中的高阶贝里曲率。未来的一个方向是寻找具有非平庸 DDKS 数的强关联晶格模型，并探索其在量子退火、量子存储等方案中的应用。此外，高阶贝里相与非线性输运、量子几何张量等当前热门领域的交叉也是一个令人兴奋的课题。

总结

本文通过将四维陈绝缘体转化为三维参数空间上的一维 iMPS 家族，首次成功地用高阶贝里曲率算法计算了 DDKS 数，并展示了其与二阶陈数的完美对应。该项工作不仅提供了一种无需外推、天然量化的二阶陈数计算方法，更重要的是建立了一维量子多体态的高阶拓扑与四维自由费米子拓扑之间的桥梁。对磁电耦合的考察表明，高阶贝里相与陈-西蒙斯角在特定模型中等价，但在更复杂的拓扑（如 Hopf 不变量）下分道扬镳，凸显了高阶几何结构的丰富性。随着对量子几何和拓扑泵浦理解的深入，高阶贝里曲率有望成为强关联拓扑物态研究的一项重要工具。