来源论文: https://arxiv.org/abs/2605.11794v1 生成时间: May 13, 2026 13:31

莫尔分数陈绝缘体中的激子-旋子模式：量子几何与光学探测

0. 执行摘要

本研究深入探讨了莫尔分数陈绝缘体（FCIs）中的集体激发行为，首次揭示了一种独特的激子-旋子模式。该模式是由系统内禀的能带内磁旋子激发与能带间激发的强耦合（杂化）所形成，其强度受量子几何性质严格调控，并产生可被光学实验直接探测到的双峰光谱特征，为理解和探测拓扑量子流体提供了新途径。

1. 核心科学问题，理论基础，技术难点，方法细节

1.1 核心科学问题

分数陈绝缘体（FCIs）作为在零磁场下实现分数量子霍尔（FQH）效应的量子物质新范式，其基态性质已被广泛研究。然而，与传统FQH体系中相对成熟的理论描述相比，FCIs中的集体激发，特别是其动力学响应，仍然是一个有待深入探索的领域。主要问题在于：

  1. 激发的本质： FCIs中是否存在类似于FQH体系的标志性低能中性激发——磁旋子（magneto-roton）？如果存在，其性质有何不同？
  2. 能带间效应的引入： 与简单的朗道能级不同，莫尔超晶格赋予了FCIs复杂的能带结构，包括近邻和远程能带。这些能带间激发（interband excitations）如何与传统的能带内激发（intraband excitations）相互作用？这种能带间耦合对FCIs的激发谱会产生怎样的独特影响？
  3. 可探测性： FCIs中的激发模式是否具有可探测的光学信号？尤其对于通常光学惰性的磁旋子，是否存在机制使其变得光学活跃，从而提供实验探针？

本研究的核心目标正是揭示莫尔FCIs中能带内与能带间激发的相互作用机制，识别由此产生的独特激发模式，并预测其可观测的光学特征，以此为莫尔FCIs的实验探测提供理论指导。

1.2 理论基础：两能级模型与量子几何

为了理解FCIs中能带内和能带间激发的相互作用，研究团队首先提出了一个简洁而富有洞察力的两能级模型。该模型将基态 $|GS angle$ 的低能单粒子-空穴激发分为两类：

  1. 能带内激发 (Intraband Sector)： 表示为 $|\Psi_{ ext{intra}}(\mathbf{q}) angle = C^\dagger_{v,\mathbf{k}+\mathbf{q}} C_{v,\mathbf{k}} |GS angle$ 的叠加态。其中 $C^\dagger_{v,\mathbf{k}+\mathbf{q}}$ 创建一个位于价带中动量为 $\mathbf{k}+\mathbf{q}$ 的粒子，$C_{v,\mathbf{k}}$ 湮灭一个位于价带中动量为 $\mathbf{k}$ 的空穴。这些集体叠加态构成了磁旋子模式 [19, 22]，在纯能带内情况下通常是光学惰性的。
  2. 能带间激发 (Interband Sector)： 表示为 $|\Psi_{ ext{inter}}(\mathbf{q}) angle = C^\dagger_{c,\mathbf{k}'+\mathbf{q}} C_{v,\mathbf{k}'} |GS angle$ 的叠加态。其中 $C^\dagger_{c,\mathbf{k}'+\mathbf{q}}$ 创建一个位于导带中动量为 $\mathbf{k}'+\mathbf{q}$ 的粒子，$C_{v,\mathbf{k}'}$ 湮灭一个位于价带中动量为 $\mathbf{k}'$ 的空穴。这些模式描述了从莫尔价带到导带的光学跃迁 [37]。

在非相互作用极限下，这两类激发是正交的。然而，莫尔体系中的库仑相互作用会明确地将它们耦合起来。值得注意的是，在扭曲二硫化钼（tMoTe$_2$）等体系中，库仑相互作用强度尺度 $U_c = e^2 / (4\pi\varepsilon\varepsilon_0 a_M)$ （约 50 meV）与能带间隙 $E_{ ext{bg}}$ （约 15 meV）相当（参见补充材料第一节 [41]），这种强耦合足以驱动显著的杂化，从而彻底改变传统的激发结构。

为了捕捉这种相互作用，研究团队构建了一个简化版的有效两能级哈密顿量（公式1）：

$$H_{ ext{eff}}(\mathbf{q}) = egin{pmatrix} E_{ ext{intra}}(\mathbf{q}) & \Delta(\mathbf{q}) \ \Delta^*(\mathbf{q}) & E_{ ext{inter}}(\mathbf{q}) \end{pmatrix}$$

其中 $E_{ ext{intra}}(\mathbf{q})$ 和 $E_{ ext{inter}}(\mathbf{q})$ 是各自解耦模式的色散关系。核心物理量是离对角耦合项 $\Delta(\mathbf{q}) = \langle \Psi_{ ext{inter}}(\mathbf{q}) | H | \Psi_{ ext{intra}}(\mathbf{q}) angle$，它负责混合这两支裸激发模式。该耦合项主要来源于库仑相互作用中导致能带混合的部分。

量子几何的关键作用：

研究发现，该耦合项的平均场估计（公式2）与量子几何性质密切相关：

$$\Delta_{kk'}(\mathbf{q}) \sim \sum_{\mathbf{G}} V(\mathbf{q}+\mathbf{G}) I_{cv}(\mathbf{k}',\mathbf{q}+\mathbf{G}) I_{vv}(\mathbf{k},\mathbf{q}+\mathbf{G})$$

其中 $I_{nm}(\mathbf{k},\mathbf{q}) = \langle u_{n,\mathbf{k}+\mathbf{q}} | u_{m,\mathbf{k}} angle$ 是形式因子（form factor）。更重要的是，在小 $\mathbf{q}$ 极限下，耦合项的平方近似满足（公式3）：

$$\langle |\Delta_{kk'}|^2 angle_0 \propto ext{tr} g_{\mathbf{k}} \ge |\Omega_\mathbf{k}|$$

这表明角度平均的耦合强度与量子度规（quantum metric $g_{\mu u}$）和贝里曲率（Berry curvature $\Omega_k$）有关。对于拓扑平带系统，这保证了 $\Delta(\mathbf{q})$ 通常是非零的，从而在有效哈密顿量中产生鲁棒的能级排斥。这种排斥导致较低的杂化能支被推向更低的能量，形成激子-旋子模式。因此，激子-旋子模式的形成和性质与系统的量子几何属性内在相关。

1.3 技术难点与方法细节

技术难点：

  1. 多体效应的精确捕捉： 莫尔FCIs是一个强关联多体系统，其激发谱不仅涉及单粒子-空穴对，还可能涉及高阶关联和多粒子-空穴激发。传统的平均场理论不足以捕捉这些复杂的效应。
  2. 能带混合与晶格效应： 与连续FQH系统不同，莫尔FCIs具有周期性晶格势，这导致伽利略不变性破缺，使得能带混合效应在动量空间中更为复杂，并允许在 $\mathbf{q}=0$ 处能带内和能带间激发的耦合。
  3. 大尺度计算： 采用精确对角化（ED）方法计算多体激发谱，特别是对于包含能带内和能带间通道的杂化激发，需要处理庞大的希尔伯特空间和复杂的矩阵元计算。
  4. 识别与区分： 在复杂的激发谱中，精确区分纯能带内激发、纯能带间激发以及它们杂化后的模式，并量化其各自的贡献。

方法细节：ED-变分Bethe-Salpeter方程（BSE）

为了克服上述技术挑战并精确捕捉强耦合体系中的多体效应，研究团队采用了基于精确对角化的变分Bethe-Salpeter方程（BSE）方法。该方法步骤如下：

  1. 基态准备： 首先，通过对扭曲MoTe$_2$的两能带相互作用哈密顿量进行精确对角化，获得精确的FCIs基态 $|GS angle$。这确保了后续激发计算的物理基础是可靠的。

  2. 激发算符定义： 对于动量 $\mathbf{q}$，构建两类粒子-空穴激发算符：

    • 能带内算符：$O_{\mathbf{k}, ext{intra}}(\mathbf{q}) = C^\dagger_{v,\mathbf{k}+\mathbf{q}} C_{v,\mathbf{k}}$
    • 能带间算符：$O_{\mathbf{k}, ext{inter}}(\mathbf{q}) = C^\dagger_{c,\mathbf{k}+\mathbf{q}} C_{v,\mathbf{k}}$ 其中 $C^\dagger, C$ 分别是电子的产生和湮灭算符，下标 $v,c$ 分别表示价带和导带。

  3. 变分激发态： 假设集体激发态 $|\Psi_\lambda(\mathbf{q}) angle$ 可以表示为这些粒子-空穴对的相干叠加：$|\Psi_\lambda(\mathbf{q}) angle = \sum_{\mathbf{k},\alpha} \psi_{\mathbf{k},\alpha}(\mathbf{q}) O_{\mathbf{k},\alpha}(\mathbf{q}) |GS angle$。这里的 $\alpha \in \{ ext{intra, inter}\}$ 标签区分能带内和能带间激发。

  4. Bethe-Salpeter方程求解： 通过最小化该变分激发态的能量，可以得到如下的BSE（公式4）：

    $$\sum_{\mathbf{k}',eta} [H^{\alphaeta}_{\mathbf{kk}'}(\mathbf{q}) - \varepsilon_\lambda(\mathbf{q}) S^{\alphaeta}_{\mathbf{kk}'}(\mathbf{q})] \psi_{\mathbf{k}',eta}(\mathbf{q}) = 0$$

    这是一个广义的本征值问题，其中：

    • $\varepsilon_\lambda(\mathbf{q})$ 是第 $\lambda$ 模式的激发能量。
    • $H^{\alphaeta}_{\mathbf{kk}'}(\mathbf{q}) = \langle GS | O^\dagger_{\mathbf{k},\alpha}(\mathbf{q}) (H - E_{GS}) O_{\mathbf{k}',eta}(\mathbf{q}) | GS angle$ 是微观多体哈密顿量 $H$ 的投影核（kernel），$E_{GS}$ 是基态能量。
    • $S^{\alphaeta}_{\mathbf{kk}'}(\mathbf{q}) = \langle GS | O^\dagger_{\mathbf{k},\alpha}(\mathbf{q}) O_{\mathbf{k}',eta}(\mathbf{q}) | GS angle$ 是重叠矩阵。

    $H$ 和 $S$ 矩阵的计算都是基于用于精确对角化的相同两能带微观模型进行的，这确保了模型的一致性，并能捕捉能带混合和晶格效应。

  5. 变分空间的限制： 为了分离和研究杂化起源，研究人员分别在纯能带内空间（裸磁旋子 $|\Psi_{ ext{intra}} angle$）、纯能带间空间（裸能带间跃迁 $|\Psi_{ ext{inter}} angle$）以及完整的两能带变分空间（杂化激子-旋子模式 $|\chi angle$）中进行计算。

微观模型细节（来自补充材料）：

  • 单谷连续哈密顿量： 使用扭曲MoTe$_2$的单谷连续哈密顿量 $H_0$（公式6），其中包含动能项、莫尔层内势 $\Delta_{t,b}(\mathbf{r})$ 和层间隧穿项 $\Delta_ au(\mathbf{r})$（公式7）。这些参数（如 $V=20.8 ext{ meV}, \psi=107.7^\circ, w=-23.8 ext{ meV}$）均来源于第一性原理计算 [9] 的拟合结果。
  • 相互作用哈密顿量： 完整的相互作用哈密顿量 $H$（公式8）在 $H_0$ 的基础上，加入了双栅极屏蔽库仑相互作用项 $V(\mathbf{q})$。库仑相互作用强度 $U_c \sim 50 ext{ meV}$，莫尔带隙 $E_{ ext{bg}} \sim 15 ext{ meV}$。

与FQH体系的比较： 为了突出莫尔FCIs中杂化现象的独特之处，研究团队还进行了周期性边界条件下FQH体系的ED-BSE计算，使用最低朗道能级（LLL）和第一激发朗道能级（1LL）。通过对比，揭示了FCIs中伽利略不变性破缺对激发谱的关键影响。

2. 关键 benchmark 体系，计算所得数据，性能数据

2.1 关键 Benchmark 体系

本研究选用的主要benchmark体系是扭曲二硫化钼（twisted MoTe$_2$），它是一种在实验上已被证实可以形成分数陈绝缘体的莫尔超晶格材料 [1,2]。研究参数设置如下：

  • 扭曲角 (Twist Angle)： $ heta = 3.5^\circ$。这是一个典型的可以实现FCI相的扭曲角。
  • 相互作用强度 (Interaction Strength)： 通过介电常数 $\varepsilon$ 来调节。论文关注两个主要区域：
    • 弱能带混合区域： $\varepsilon = 15$。此时库仑相互作用相对较弱，能带间隙相对较大，能带混合效应不显著。
    • 强能带混合区域： $\varepsilon = 5$。此时库仑相互作用相对较强，能带间隙相对较小，能带混合效应显著。
  • 计算集群： 为了进行精确对角化，采用的是一个有限动量空间集群，最大系统尺寸为 $N=18$（例如 $3 imes 6$）。

2.2 计算所得数据与性能分析

激发色散关系（图2a,b）：

  • 弱能带混合区域 ($ heta=3.5^\circ, \varepsilon=15$，图2a)：

    • 长波长能带间激发能隙 $E_{ ext{inter}} = \varepsilon_{ ext{inter}}(0) \approx 15.04 ext{ meV}$，仅略高于非相互作用能带间隙 $E_{ ext{bg}} \approx 12.32 ext{ meV}$。这表明在弱相互作用下，能带间跃迁与FQH体系的Kohn模式行为相似。
    • 能带内分支与能带间分支能量间隔较大，能带内激发保留了磁旋子的特征。
    • 杂化程度较弱。最低杂化模式 $|\chi angle$ 与能带内分支几乎重合，仅含有少量能带间成分，其贡献为 $|\langle \Psi_{ ext{inter}} | \chi angle|^2 \approx 4.56\%$（图2g），表明其光学活性非常微弱。

  • 强能带混合区域 ($ heta=3.5^\circ, \varepsilon=5$，图2b)：

    • 能带间和能带内分支的能量接近，并在 $q \approx 0.02 ext{ Å}^{-1}$ 处发生能量交叉。这意味着在长波长极限下，它们的耦合变得非常强。
    • 杂化导致显著的能级排斥，混合模式的能量被向下推移了 $\Delta E_{ ext{inter}} - E_{\chi} \approx 1.86 ext{ meV}$，与两能级图像中的能级排斥预期一致。
    • 结果是，杂化模式（激子-旋子）的能量低于能带间跃迁的能隙，并在有限动量 $q_{ ext{min}} \approx 0.06 ext{ Å}^{-1}$ 处发展出旋子最小值。这一模式同时具有激子和旋子特征，因此被命名为激子-旋子模式。

多体激发谱（图2e,f）：

  • 弱混合区域 ($\varepsilon=15$，图2e)： 激发谱呈现了FCIs预期的三重基态简并，与FQH系统相似。GMP模式和纯能带内模式的能量显著提高。最低杂化模式主要由能带内成分构成。
  • 强混合区域 ($\varepsilon=5$，图2f)： 激发谱与弱混合情况明显不同。GMP模式和纯能带内模式的能量被大幅提升。杂化分支在旋子最小值附近仍是最低的中性激发。这表明能带间成分对于理解FCIs的低能激发谱至关重要。

与能隙相关的激发能与重叠度演化（图2c,d）：

  • FCIs（图2c）： 莫尔FCIs中，随着能带间隙 $E_{ ext{bg}}$ 的调谐，能带间激发 $E_{ ext{inter}}$ 和能带内激发 $E_{ ext{intra}}$ 的能量以及杂化模式与纯能带间模式的重叠度 $|\langle \Psi_{ ext{inter}} | \chi angle|^2$ 呈现平滑连续的演化。在 $q=0$ 处，能带间激发能 $E_{ ext{inter}}$ 几乎线性地跟踪能带间隙 $E_{ ext{bg}}$，而能带内激发能 $E_{ ext{intra}}$ 则对 $E_{ ext{bg}}$ 不敏感。这表明莫尔FCIs中即使在 $q=0$ 处也存在鲁棒的能带间-能带内杂化。
  • FQH系统（图2d）： 为了进行比较，研究团队还计算了FQH系统的激发谱。FQH系统在 $q=0$ 处展示了简单的能级交叉，而不是平滑的杂化。这是由于伽利略不变性（Kohn定理 [30]）的存在，导致 $q=0$ 处的质心运动与相对运动激发严格解耦。最低杂化态 $|\chi angle$ 在 $q=0$ 处与两个裸分支中的一个重合，而不是它们的混合态，导致 $|\langle \Psi_{ ext{inter}} | \chi angle|^2$ 出现骤然切换。

结论： FCIs中长波长处（$q=0$）的平滑杂化是一种真正的晶格效应，它打破了伽利略不变性，从而允许能带间和能带内激发的耦合。

莫尔参数对杂化模式光学活性的调控（图2g）：

  • 该图展示了 $q=0$ 处长波长杂化模式的光学活性，通过 $|\langle \Psi_{ ext{inter}} | \chi angle|^2$ 来衡量。随着相互作用强度 $\varepsilon$ 的减小（即相对能带间隙，库仑相互作用增强），光学活性显著增加。在 $\varepsilon=15$ 时仅为约 $4.56\%$，而在 $\varepsilon=5$ 时则高达约 $65.55\%$。这表明激子-旋子模式的光学活性可以通过莫尔带隙和相互作用强度进行有效调控。

光学电导率 Re$\sigma(\omega)$ 与双峰光谱特征（图3a,b,c）：

  • 通过Kubo公式（公式5）计算了FCIs的纵向光学电导率 Re$\sigma(\omega)$，以预测其光学响应。
  • 弱混合区域 ($\varepsilon \ge 10$，图3a)： 光谱主要表现为一个主导峰（Peak 2），对应于能带间跃迁的连续谱。最低的杂化模式（Peak 1）几乎是光学惰性的，与图2g中其较小的能带间权重一致。
  • 强混合区域 ($\varepsilon \le 6$，图3b)： 进入强混合区域后，激发谱发生了剧烈重构。强杂化导致两个分支之间发生大规模的能级排斥。较低的能支（Peak 1）被向下推移，并从连续谱中分离出来，形成了能带间隙内离散且稳定的束缚激子态。
  • 伴随这一形成过程的是光谱权重的显著转移：Peak 1 迅速变得明亮（在 $\varepsilon=5$ 时，其光学活性增至约 $65.55\%$）。
  • 演化过程（图3c）： 展示了从 $\varepsilon=7$ 到 $\varepsilon=3$ 逐渐增强的相互作用下，Re$\sigma(\omega)$ 如何从单峰结构演变为独特的双峰结构。Peak 1（激子-旋子模式）变得越来越明显，并向低能量移动，而Peak 2（能带间连续谱）的位置保持相对稳定。这个双峰结构是激子-旋子模式的关键实验指纹。

实验可探测性：

  • 在 $ heta=3.5^\circ, \varepsilon=5$ 的条件下，两个杂化能支的能量间隔 $\Delta \approx 1.9 ext{ meV}$，位于 $\hbar\omega \sim 15-20 ext{ meV}$ 范围内。这正好处于太赫兹（THz）到远红外（3-5 THz）光谱窗，目前的时域THz光谱和傅里叶变换红外吸收技术均能达到亚毫电子伏特（sub-meV）的分辨率。
  • 这种 $1.9 ext{ meV}$ 的分裂也远超FCIs工作温度（$T \le 2 ext{ K}$）下的热展宽，确保了双峰在光谱上是可分辨的。

有限尺寸效应（补充材料图4）：

  • 研究团队评估了ED计算中的有限尺寸效应，通过对比不同系统尺寸（$N=12, 15, 18$）下激发的色散关系。结果显示，在两种相互作用区域（强混合和弱混合）中，不同系统尺寸的数据点很好地落在共同曲线上，表明磁旋子最小值的位置、能带间模式的能量尺度以及长波长行为对系统尺寸不敏感。这证实了文章中提出的激子-旋子现象不是有限尺寸集群的产物。

全面光学电导率与能带内贡献（补充材料图5）：

  • 补充材料图5对比了完整的多体光学电导率 Re$\sigma_{ ext{full}}(\omega)$ 与仅包含能带内贡献的 Re$\sigma_{ ext{intra}}(\omega)$。结果显示，Re$\sigma_{ ext{intra}}(\omega)$ 在整个频率窗口内都远小于 Re$\sigma_{ ext{full}}(\omega)$。这支持了主要的结论，即集体模式的光学权重主要来源于能带间贝里联络（Berry connection）通道，能带内激发本身是光学惰性的。即使在周期势折叠下，能带内理论 [21,52] 也仅能产生微弱的光学活性，再次强调了能带间混合效应在强混合区域中的主导地位。

2.3 性能数据

本研究的所有计算均在天津国家超级计算中心的天河新一代超级计算机上进行。精确对角化和BSE计算对于多体系统是计算密集型的。尽管论文中未详细说明具体的计算耗时、CPU/GPU核时、内存使用量等性能指标，但使用国家级超算资源表明了这些计算的复杂性和所需的计算规模。ED计算的系统尺寸为$N=18$，这对于精确对角化来说已是相当大的规模，但为了确保结果的可靠性，ED-BSE方法通常需要高效的并行算法和优化实现。

3.1 核心计算方法与实现细节

本研究的核心计算方法是**精确对角化（Exact Diagonalization, ED）与变分Bethe-Salpeter方程（Bethe-Salpeter Equation, BSE）**的结合。这两种方法在强关联量子多体物理中是标准但计算要求极高的技术。

  1. 基态精确对角化：

    • 哈密顿量构建： 首先，需要根据莫尔扭曲MoTe$_2$的连续哈密顿量 $H_0$（公式6，7）和库仑相互作用项构建完整的两能带相互作用哈密顿量 $H$（公式8）。这涉及将实空间莫尔势和隧穿项傅里叶变换到动量空间，并投影到 Bloch 波函数基上。
    • 多体基态求解： 在一个有限的动量空间集群（例如 $N=18$ 个电子的 $3 imes 6$ 集群）中，构建多体哈密顿量的矩阵表示。对于分数陈绝缘体，这通常涉及一个特定填充因子下的自旋极化电子体系。随后，通过稀疏矩阵对角化算法（例如 Lanczos 算法），找到哈密顿量的最低本征态，即FCI基态 $|GS angle$ 及其能量 $E_{GS}$。这一步是后续激发计算的基础。

  2. 激发谱的Bethe-Salpeter方程求解：

    • 粒子-空穴激发算符： 基于能带的 Bloch 波函数，定义能带内（$O_{\mathbf{k}, ext{intra}}(\mathbf{q})$）和能带间（$O_{\mathbf{k}, ext{inter}}(\mathbf{q})$）粒子-空穴激发算符。这些算符通过电子在价带和导带之间的跃迁来构建。
    • 矩阵元计算： BSE 的核心在于构建哈密顿量投影核 $H^{\alphaeta}_{\mathbf{kk}'}(\mathbf{q})$ 和重叠矩阵 $S^{\alphaeta}_{\mathbf{kk}'}(\mathbf{q})$（公式4）。这些矩阵元涉及复杂的四费米子算符在基态 $|GS angle$ 上的期望值，这需要精确计算 Bloch 波函数之间的形式因子 $I_{nm}(\mathbf{k},\mathbf{q})$ 和库仑相互作用势 $V(\mathbf{q})$。特别是，$H^{\alphaeta}_{\mathbf{kk}'}(\mathbf{q})$ 包含了 $(H-E_{GS})$ 项，这意味着它考虑了基态能量相对于激发态能量的差异以及多体相互作用的贡献。
    • 广义本征值问题： 将 $H$ 和 $S$ 矩阵构建完成后，对于每个动量 $\mathbf{q}$，求解广义本征值问题，得到激发能量 $\varepsilon_\lambda(\mathbf{q})$ 和对应的变分激发态波函数 $\psi_\lambda(\mathbf{q})$。这些波函数描述了集体激发的组成，包括能带内和能带间成分的混合比例。

  3. 光学电导率计算：

    • Kubo 公式： 利用 Kubo 公式（公式5）计算纵向光学电导率 Re$\sigma(\omega)$。这需要计算速度算符 $\hat{v}$ 在基态 $|GS angle$ 和每个激发态 $|\Psi_n angle$ 之间的矩阵元 $\langle \Psi_n | \hat{v} | \Psi_0 angle$。
    • 速度算符投影： 速度算符 $\hat{v}$ 会被投影到低能子空间，该子空间由顶价带和最低导带构成。在 Bloch 基中，速度算符可以分解为两部分：能带内群速度项 $\hat{v}_{ ext{intra}}$ 和能带间贝里联络项 $\hat{v}_{ ext{inter}}$（公式15）。其中，能带间贝里联络项与能带间贝里联络 $A_{cv}(\mathbf{k})$ 成正比，是产生光学活性的关键。
    • 光谱展宽： 为了模拟实验测量，Dirac delta 函数通常用 Lorentzian 峰形替代，引入一个展宽参数 $\gamma$（例如 $\gamma = 0.5 ext{ meV}$）。

3.2 复现指南

要复现本研究结果，需要遵循以下主要步骤：

  1. 设置莫尔哈密顿量：

    • 莫尔势： 根据论文补充材料中的公式6和7，实现扭曲MoTe$_2$的单谷连续哈密顿量 $H_0$。这包括计算莫尔势 $\Delta_{t,b}(\mathbf{r})$ 和层间隧穿项 $\Delta_ au(\mathbf{r})$。所有参数（$V, \psi, w, m^*, a_0$）均已在补充材料中给出，可直接使用。
    • 布里渊区离散化： 在莫尔布里渊区内离散化动量空间，构建单粒子哈密顿量矩阵，并通过对角化获得 Bloch 波函数 $u_{n,\mathbf{k}}$ 和单粒子能带色散 $E_n(\mathbf{k})$。

  2. 计算相互作用矩阵元：

    • 库仑势： 实现双栅极屏蔽库仑相互作用势 $V(\mathbf{q})$。论文中给出了其表达式 $V(\mathbf{q}) = e^2 anh(q|d) / (A \pi \varepsilon \varepsilon_0 q)$。
    • 形式因子： 计算所有必要的能带间和能带内形式因子 $I_{nm}(\mathbf{k},\mathbf{q}) = \langle u_{n,\mathbf{k}+\mathbf{q}} | u_{m,\mathbf{k}} angle$。这些因子是莫尔哈密顿量中相互作用项投影到 Bloch 基上的关键。

  3. 构建多体哈密顿量并进行精确对角化：

    • 多体基： 选择合适的粒子数 $N$ 和集群尺寸（例如 $N=18$），构建多体 Fock 空间基态。对于FCIs，需要考虑填充因子 $ u=1/3$。
    • 哈密顿量矩阵： 将单粒子能带哈密顿量和相互作用项（包含能带混合）转换到多体基中，构建多体哈密顿量 $H$ 的稀疏矩阵表示。
    • 基态求解： 对 $H$ 进行精确对角化，获取最低本征态作为FCIs基态 $|GS angle$。

  4. Bethe-Salpeter方程计算：

    • 算符矩阵元： 对于每一个动量 $\mathbf{q}$，计算粒子-空穴激发算符 $O_{\mathbf{k},\alpha}(\mathbf{q})$ 在基态上的投影核 $H^{\alphaeta}_{\mathbf{kk}'}(\mathbf{q})$ 和重叠矩阵 $S^{\alphaeta}_{\mathbf{kk}'}(\mathbf{q})$。这是计算最密集的部分，需要精细处理多体矩阵元。
    • 求解 BSE： 求解广义本征值问题 $[H(\mathbf{q}) - \varepsilon(\mathbf{q})S(\mathbf{q})] \psi(\mathbf{q}) = 0$，获取激发能量 $\varepsilon_\lambda(\mathbf{q})$ 和相应的激发波函数 $\psi_\lambda(\mathbf{q})$。

  5. 光学电导率计算：

    • 速度算符： 计算速度算符 $\hat{v}$ 在 Bloch 基下的矩阵元，特别关注能带间贝里联络项。
    • Kubo 公式： 利用 Kubo 公式（公式5），结合基态和激发态的能量及速度算符矩阵元，计算 Re$\sigma(\omega)$。应用 Lorentzian 展宽来模拟实验峰形。

本研究论文未提供其计算代码的开源仓库链接。这在计算凝聚态物理领域较为常见，因为这些代码通常是高度专业化且基于内部开发框架的。然而，要实现上述计算，通常会依赖以下类型的软件包和库：

  • 数值计算核心库：

    • NumPy / SciPy (Python)： 提供高效的数组操作、线性代数、稀疏矩阵处理和傅里叶变换功能。
    • LAPACK / BLAS (Fortran/C)： 底层高性能线性代数库，通常用于矩阵对角化。
    • ARPACK / SLEPc (C/Fortran)： 用于大规模稀疏矩阵的本征值问题求解（例如 Lanczos 或 Arnoldi 算法）。

  • 量子多体物理库：

    • QuSpin (Python)： 一个用于精确对角化和量子动力学模拟的开源库，特别适用于自旋链和费米子系统。
    • OpenFermion (Python)： 尽管主要用于量子计算中的费米子哈密顿量，但其哈密顿量构建工具也可能对多体哈密顿量的构造有帮助。
    • 定制ED代码： 许多研究团队会基于 Fortran, C++ 或 Python 自行开发高度优化的ED代码，以适应其特定的哈密顿量结构和计算需求。本研究很可能使用了此类定制的内部代码。

  • 数据分析与可视化：

    • Matplotlib (Python)： 用于生成论文中所示的各种色散关系图和光谱图。

  • 高性能计算环境：

    • 本研究明确指出计算在“天津国家超级计算中心的天河新一代超级计算机”上进行。这意味着代码需要支持并行计算（例如 MPI 或 OpenMP），以利用大规模计算资源。

由于没有提供具体的代码链接，任何复现工作都需要研究人员投入大量精力从头开始实现ED-BSE框架，或者基于现有开源库进行大量定制和扩展。这要求对量子多体理论、数值方法和高性能计算有深刻的理解。

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献及其相关性

本研究建立在凝聚态物理学，特别是分数量子霍尔效应（FQH）和莫尔超晶格领域的深厚基础上。以下是论文中一些关键引用文献及其对本工作的支撑作用：

  1. [1, 2] J. Cai et al. (2023), H. Park et al. (2023)： 这些是关于在扭曲MoTe$_2$中观察到分数量子反常霍尔效应的里程碑式实验论文。它们证实了扭曲MoTe$_2$是实现FCI的优异平台，为本研究的理论预测提供了重要的实验背景和材料基础。
  2. [19] X. Shen et al. (2026)： 这是作者团队的早期工作，专注于莫尔分数陈绝缘体中的磁旋子。本研究在此基础上，通过引入能带间激发，极大地扩展了对FCI集体激发的理解，揭示了磁旋子与能带间激发的杂化。
  3. [22] S. M. Girvin, A. H. MacDonald, and P. M. Platzman (1986)： 这篇经典论文建立了FQH效应中磁旋子的理论框架（GMP理论）。本研究将FCI中的激发与GMP理论进行对比，突出了莫尔晶格效应（特别是伽利略不变性破缺）导致的独特现象。
  4. [24] F. D. M. Haldane (2011)： 这篇论文提出了FQH效应的几何描述，强调了量子几何在理解拓扑相中的重要性。本研究进一步揭示了量子几何如何直接控制能带内和能带间激发的耦合，从而影响集体激发谱。
  5. [30] W. Kohn (1961)： Kohn定理是连续系统中质心运动与相对运动解耦的基石。本研究通过与FQH系统的对比（图2d），明确展示了莫尔势如何打破伽利略不变性，使得FCIs中能带内和能带间激发在 $q=0$ 处发生杂化，这是与传统FQH系统最显著的区别之一。
  6. [37] W.-X. Qiu and F. Wu (2024)： 该文探讨了陈绝缘体中手性激子光学响应如何揭示量子几何。这与本研究中激子-旋子模式的光学活性预测紧密相关，为实验探测提供了理论视角。
  7. [42] J. Wang et al. (2021)： 该研究提供了在平带中几何和相互作用的朗道能级精确描述，为本工作处理莫尔平带中的相互作用和量子几何提供了理论工具。

4.2 本工作局限性评论

尽管本研究对莫尔FCIs的集体激发提供了深刻的见解和重要的理论预测，但仍存在一些局限性：

  1. 变分空间的限制与高阶关联效应：

    • 本研究主要关注单粒子-空穴激发组成的变分空间。论文本身在第三页脚注中指出，一些接近 $\Gamma$ 点的低能激发模式（例如四极模式 [45-48] 或手性引力子模式 [65]）未被当前的BSE变分空间完全捕获。虽然补充材料中图5的结果表明这些模式的光学活性可能较弱，但一个更完整的激发谱理论仍需考虑这些高阶或多粒子-空穴激发。
    • 两能级模型的物理解释虽然直观，但它简化了多能带系统中的复杂多体相互作用。论文中也提到了两能级模型与完整ED计算在交叉参数（如 $\varepsilon$）估计上的差异，这表明更复杂的关联效应在精确定量上可能仍有影响。

  2. 有限尺寸效应与热力学极限：

    • 精确对角化计算是在有限的动量空间集群上进行的（例如 $N=18$ 个电子）。尽管补充材料图4展示了主要定性特征对系统尺寸不敏感，但有限尺寸计算始终无法完全替代热力学极限。某些长程关联效应或临界现象在小系统中可能难以准确捕捉。

  3. 材料特异性与普适性：

    • 本研究主要集中在扭曲MoTe$_2$体系。虽然理论框架具有通用性，但激子-旋子模式的出现及其具体性质（如杂化强度、光学活性）与MoTe$_2$的特定材料参数（能带隙、相互作用强度、形式因子）紧密相关。将结论推广到其他莫尔FCI材料需要进一步的研究，因为不同材料可能具有不同的能带结构和量子几何特征。

  4. 缺乏开源代码与可复现性：

    • 论文中未提供其计算代码的开源链接或详细的实现手册。这使得其他研究团队直接复现本工作面临挑战，需要投入大量时间和精力独立开发或重新实现复杂的ED-BSE框架。这在一定程度上限制了研究的透明度和可复现性。

  5. 理论预测的实验验证：

    • 本研究提出了明确的实验指纹——双峰光学光谱特征，并建议了可行的实验技术（THz光谱）。然而，这仍然是一个理论预测，最终的验证仍需依赖未来的实验观测。实验中的非理想因素，如缺陷、无序或应变，可能会影响实际观测结果与理论预测的符合程度。

这些局限性并不削弱本研究的重大贡献，而是为未来的研究指明了方向，鼓励更全面、更普适、更精确的理论探索和实验验证。

5. 其他必要补充

5.1 广泛影响与科学意义

本研究在理解莫尔分数陈绝缘体（FCIs）的集体激发方面迈出了重要一步，其科学意义深远，具有广泛影响：

  1. 开启FCI动力学研究的新范式： 莫尔FCIs在零磁场下实现分数量子霍尔物理，为探索新颖的拓扑量子物态提供了独特平台。此前研究主要集中于基态性质，而本工作则将焦点转向了其动力学响应。通过揭示激子-旋子模式，本研究为FCI的动力学研究奠定了基础，拓宽了我们对拓扑量子物质复杂行为的理解。
  2. 超越传统FQH的物理： 与具有简单朗道能级的传统FQH系统不同，莫尔FCIs因其周期性晶格势而打破了伽利略不变性，并拥有复杂的能带间结构。这种差异使得FCIs能够支持在连续FQH系统中不存在的集体激发模式，例如本研究发现的激子-旋子模式，极大地丰富了拓扑量子物质的物理内涵。
  3. 量子几何的深层作用： 本研究明确建立了集体激发杂化强度与量子几何性质（如量子度规和贝里曲率）之间的内在联系。这强调了本征能带的几何结构，而不仅仅是能量，在调控多体集体激发中的关键作用。这一发现为“能带工程”提供了新的设计维度，即通过精细调整莫尔结构或外部条件来塑造量子几何，进而控制新奇的量子现象。
  4. 直接的光谱学探针： 预测到的独特双峰光学光谱特征为实验探测FCIs中的集体激发提供了一个具体且直接的方法。这一发现使得研究人员可以从基态的间接测量转向对激发态的直接观测，极大地推动了莫尔拓扑量子材料的实验研究进展。
  5. 促进多学科交叉： 本工作连接了凝聚态物理的多个前沿领域，包括拓扑物理、强关联物理、莫尔材料科学以及量子光学，有望激发更多跨学科的合作与创新。

5.2 未来工作与开放问题

本研究虽然取得了重要突破，但也为未来的探索留下了广阔空间：

  1. 更全面的激发谱： 除了激子-旋子模式，FCIs中可能还存在其他低能多体激发，如四极模式、手性引力子模式 [45-48, 65]、莫尔激子 [63, 66-68] 等。对这些模式的详细研究及其光学特征的识别，将有助于构建一个更完整的FCIs激发图景。
  2. 探索其他FCI类型： 将当前框架应用于其他类型的分数陈绝缘体，例如Jain序列FCIs [9, 18]、具有更高陈数的FCIs [69, 70] 或非阿贝尔FCIs [71-74]，可能会揭示更丰富、更复杂的激子-旋子物理，甚至可能涉及新型的拓扑物态。
  3. 任意子-激子复合体（Anyon-Trion Complexes）： 论文提到了磁旋子与其他激子类型耦合形成“任意子-激子”复合体的可能性 [28, 60-63]。这是一个引人注目的方向，可以将集体激发与准粒子相互作用联系起来，深入理解分数统计量的表现。
  4. 实验验证与材料调控： 最紧迫的下一步是实验上观测到预测的双峰光学信号。此外，通过莫尔带隙工程（例如调整扭曲角 [55-58]、施加静水压力 [54]、改变介电屏蔽 [53]）来系统调控激子-旋子模式，并对其动力学行为进行更深入的实验研究。
  5. 无序和应变的影响： 真实的莫尔体系不可避免地存在无序和应变。研究这些因素对激子-旋子模式及其光学信号的影响，对于实验观测和实际应用具有重要意义。

5.3 量子几何与莫尔势的详细阐述

量子几何的深层机制：

本研究的核心创新之一在于明确指出能带内与能带间激发的杂化强度 $\Delta(\mathbf{q})$ 与量子度规 $g_{\mu u}$ 和贝里联络 $A_{cv}$ 密切相关。在 Bloch 理论中，量子度规 $g_{\mu u}$ 描述了 Bloch 波函数在动量空间中的几何“形状”和距离，而贝里联络 $A_{cv} = i\langle u_{c,\mathbf{k}} | abla_{\mathbf{k}} u_{v,\mathbf{k}} angle$ 则表征了不同能带间 Bloch 波函数在动量空间中的连通性。形式因子 $I_{nm}(\mathbf{k},\mathbf{q})$ 本质上是 Bloch 波函数在不同动量下的重叠，其小 $\mathbf{q}$ 展开直接包含了贝里联络和量子度规的信息。因此，杂化强度 $\Delta(\mathbf{q})$ 并非仅仅取决于能带的能量色散，更是由这些能带的内在几何属性所决定。这意味着研究人员可以通过调控莫尔超晶格的结构参数（如扭曲角、应变）来改变能带的量子几何，从而精确控制集体激发的耦合和性质。

莫尔势对伽利略不变性的破缺：

在连续介质的FQH体系中，Kohn定理是一个基本原理，它指出在 $q=0$ 极限下，电子的质心运动与内部相对运动是完全解耦的。这意味着在 $q=0$ 处，能带内激发（如磁旋子）和能带间激发（如激子）之间不会发生耦合。FQH体系的激发谱因此表现为能级交叉（图2d）。然而，莫尔超晶格为FCIs引入了一个周期性的晶格势。这个势打破了体系的连续平移对称性，进而破缺了伽利略不变性。伽利略不变性的破缺是本研究中激子-旋子模式出现的关键。它使得能带内和能带间激发即使在 $q=0$ 处也能够通过库仑相互作用发生耦合，从而导致了莫尔FCIs中平滑的杂化过程和激子-旋子模式的形成。这种晶格效应是莫尔FCIs独有的，使其集体激发行为与传统的FQH体系截然不同。

5.4 技术影响

本研究的理论发现不仅具有深刻的科学价值，还蕴含着潜在的技术应用：

  1. 可调谐光学器件： 激子-旋子模式的光学活性可以被莫尔带隙和相互作用强度高效调控。这种在太赫兹/远红外波段具有可调谐光学响应的特性，为开发新型光电器件（如调制器、探测器或新型量子传感器）提供了可能性。想象一下，基于“量子几何开关”来控制特定频率光的吸收或发射。
  2. 新型量子材料平台： 莫尔FCIs作为一种新兴的量子材料平台，结合了拓扑和强关联效应。本研究提供的集体激发图景，将有助于设计和优化这些材料，以实现更复杂的量子功能或信息处理。对激子-旋子模式的理解可能促进构建基于集体激发的量子器件。
  3. 复杂量子物质的诊断工具： 预测到的双峰光谱特征为实验物理学家提供了一个强大的诊断工具，用于探测和表征莫尔FCIs中复杂的多体动力学。通过分析这些光谱特征，可以更深入地理解拓扑量子物质的微观机制，甚至可能识别出此前未知的激发类型。

5.5 个人见解与展望

本研究在莫尔分数陈绝缘体（FCIs）领域取得了显著进展，它不仅仅是对特定材料的计算研究，更提供了一个理解复杂量子多体系统集体激发的一般性框架。其最令人兴奋之处在于：

  • 揭示隐藏的模式： 通过将能带内磁旋子和能带间激发的相互作用统一处理，研究团队成功揭示了一种全新的激子-旋子模式。这种模式在传统理论中是被忽略的，它的发现拓展了我们对拓扑量子物态中激发谱的认知。
  • 量子几何的直观联系： 将杂化强度与量子几何性质直接联系起来，为理论物理学提供了一个优雅的工具，也为材料科学家提供了新的设计思路。这意味着我们可以通过精细的莫尔结构调控，来“雕刻”能带的几何形状，进而定制量子材料的宏观功能。
  • 实验可达性： 双峰光谱特征的清晰预测，以及现有实验技术（太赫兹光谱）的适配性，使得本研究具有很高的实验验证潜力。这有望在不久的将来，将理论预测转化为激动人心的实验发现，从而进一步推动莫尔量子材料领域的发展。

总而言之，这项工作不仅深化了我们对FCIs基本物理的理解，也为未来在拓扑量子材料中实现功能化量子器件提供了理论基础和实验指引。它再次证明了，在多体强关联和晶格效应并存的复杂系统中，仍有丰富的、出人意料的物理等待我们去发现和利用。