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

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

问题陈述
拓扑不变量对于分类量子多体基态至关重要。虽然贝里曲率（一种二形式）及其相关的陈数成功地表征了零维基态族（例如布洛赫态在布里渊区中），但近期的理论发展表明，这一框架可以推广到空间扩展系统。在该框架下，$d$ 维空间中的短程纠缠态族由更高阶的微分形式（例如一维系统的三形式）来表征。具体而言，对于由三维空间参数化的一维系统，相关的量是不变量 Dixmier–Douady–Kapustin–Spodyneiko (DDKS) 数 $\nu$，它源自于一种高阶贝里曲率。

目前存在一个显著差距，即如何对通用的相互作用或自由费米子系统进行显式且具有流形量化特性的高阶不变量计算。虽然对于特定的自由费米子模型（如四维狄拉克模型）存在解析表达式，但通过格点规范理论或实空间磁场响应来计算四维中第二陈数 ($ch_2$) 的数值方法通常缺乏“流形量化”特性，这意味着它们仅在连续极限下收敛到整数值，或者需要极大规模的网格细化，而不是通过构造本身返回整数。此外，在三维空间中，DDKS 数的物理诠释，特别是其与轨道磁电耦合及 Chern-Simons 轴子角 ($\theta_{CS}$) 之间的关系，仍是一个开放性问题。

方法论
作者通过将四维（4D）陈绝缘体的格点模型重写为由三维布里渊区 ($T^3$) 参数化的平移不变无限一维（1D）链族来解决这些问题。

1. 模型重构： 从四维狄拉克哈密顿量出发，作者对一个动量分量 ($k_4$) 进行部分傅里叶变换，得到依赖于参数 $\lambda = (\lambda_1, \lambda_2, \lambda_3)$ 的一维哈密顿量族 $\{H_\lambda\}$。这些哈密顿量描述了每个原胞含有两个轨道的一维链。
2. 数值框架 (iMPS)： 这些一维链的基态使用无限矩阵乘积态 (iMPS) 进行近似。作者利用了 Ref. [13] 中提出的算法，该算法推广了 Fukui-Hatsugai-Suzuki 方法以计算高阶贝里曲率。
3. 离散高阶贝里曲率：
   • 将参数空间 $T^3$ 离散化并三角剖分为四面体。
   • 对于每个顶点，通过 iDMRG 获取 iMPS 基态。
   • 计算相邻 iMPS 张量之间的转移函数，使用的是混合转移矩阵的主特征向量。
   • 定义了一个对于三角形面 $\Delta_2$ 的规范不变相位 $\phi(\Delta_2)$，并通过 Schmidt 谱进行加权以处理变化的键维数。
   • 针对四面体 $\Delta_3$ 的高阶贝里曲率 $F(\Delta_3)$ 被定义为其边界上这些相位的定向和。
   • 通过对整个 $T^3$ 的三角剖分进行求和，得到 DDKS 数 $\nu$，从而产生一个显式量化的整数。
4. 对比分析： 作者使用标准的格点规范理论方法（Ref. [34]）计算了同一模型的第二陈数 $ch_2$，并将结果进行对比。
5. 磁电耦合： 本研究调查了 DDKS 数与控制三维轨道磁电耦合的 Chern-Simons 轴子角 $\theta_{CS}$ 之间的关系。他们比较了在泵浦循环过程中，由 $\theta_{CS}$（源自解析第二 Chern 形式）的变化与由“高阶贝里相位” $\theta_{HBP}$（源自 iMPS 构建）的变化之间的关系。

核心贡献与结果

• 第二陈数 $ch_2$ 的显式量化计算： 作者证明，通过基于 iMPS 的高阶贝里曲率算法计算得到的 DDKS 数 $\nu$，生成的四维陈绝缘体拓扑相图与第二陈数 $ch_2$ 的相图完全一致。这种一致性在不同的质量机制 ($m$) 下均成立，包括标准格点规范理论计算 $ch_2$ 时会受到显著有限网格误差影响的区域。iMPS 方法即使在粗糙网格（如 $4^3$）上也能产生精确的整数，而格点规范方法则需要更精细的网格（如 $10^4$ 个点）才能收敛。
• 通过 Künneth 公式建立的理论联系： 本文利用 Künneth 公式建立了 4D 第二陈数与 3D DDKS 数之间的严格联系。通过将 4D 环面 $T^4$ 分解为 $S^1 \times T^3$，作者表明 $T^4$ 上的第二 Chern 类对应于 $H^1(S^1, \mathbb{Z})$ 的生成元与由 $\nu$ 测量的 $H^3(T^3, \mathbb{Z})$ 类之积。因此，对于这种特定的分解，$ch_2 = \nu$。
• 四维狄拉克模型中的磁电耦合： 对于四维狄拉克模型，作者发现沿泵浦循环的变化量 $\Delta\theta_{CS}$ 被 $\Delta\theta_{HBP}$ 的变化精确地重现。这表明，在布里渊区切片的第一 Chern 数消失的情况下，高阶贝里相位与磁电响应之间存在直接对应关系。
• Hopf 绝缘体中的反例： 作者在 Hopf 绝缘体（一种具有非平凡 Hopf 不变量的两能带模型）上测试了这种对应关系。在这种情况下，虽然 Chern-Simons 角由 Hopf 不变量量化，但计算出的高阶贝里相位 $\theta_{HBP}$ 既不是量化的，也不等于 $\theta_{CS}$。作者将这种差异归因于 Hopf 不变量的“微妙”性质——它严格定义仅适用于动量空间中的两能带模型，而这一条件在热力学极限的实空间 iMPS 描述中可能无法保持。

意义与主张
本文声称提供了首次利用高阶贝里曲率和 iMPS 对 4D 格点模型进行显式、显式量化计算第二陈数的数值方法。这验证了高阶贝里曲率可以作为空间扩展态族的高阶拓扑不变量探测器的理论框架。

作者将他们的结果置于量子几何与拓扑属性的背景下，指出虽然高阶贝里相位与磁电耦合的对应关系在 4D 狄拉克模型中成立，但它并非普适的，Hopf 绝缘体的情况便证明了这一点。他们强调，其方法特别适用于单粒子布洛赫描述失效的相互作用系统，尽管目前主要侧重于自由费米子基准测试。这项工作凸显了合成维度和准周期驱动在实验平台中实现并探测这些高维拓扑现象的潜力。

论文最后谦逊地总结道，关于高阶贝里曲率与三维物理响应之间的关系（特别是与 $\theta_{CS}$ 联系的普适性）仍是一个开放性问题。他们指出，高阶贝里相位可能捕捉到不同于标准 Chern-Simons 响应的几何信息，值得进一步研究其在非线性输运和基态物理中的作用。