来源论文: https://arxiv.org/abs/2603.02011v1 生成时间: Mar 03, 2026 02:12

0. 执行摘要

在凝聚态物理与量子化学的交汇点，准晶（Quasicrystal）与超级莫尔（Super-moiré）系统的激子动力学模拟一直被视为“计算天花板”。由于激子希尔伯特空间的大小随格点数 $N$ 呈平方级增长（$N^2$），传统的 Bethe-Salpeter 方程（BSE）求解器在格点数达到万级时便会因内存与计算量爆炸而失效。然而，超级莫尔系统的物理特性往往需要在数百万甚至数十亿个格点的尺度上才能完全显现。

近期，Anouar Moustaj 等人在其论文《Tensor-network methodology for super-moiré excitons beyond one billion sites》中提出了一种革命性的张量网络（TN）框架。该方法结合了实空间 BSE 形式、Quantics 张量网络（QTN）编码以及高阶 Chebyshev 核多项式算法。核心突破在于：通过交错排序（Interleaved Ordering）策略，将原本非局域的电子-空穴相互作用转化为张量网络中的局域算符，从而在保持低键维度的前提下，实现了对维度高达 $10^{18}$ 的有效激子 Hamiltonian 的直接处理。这使得在超过 10 亿个格点的 1D 和 2D 系统中观测到激子能带分裂与空间受限现象成为可能，为模拟大规模非周期性量子物质开辟了新路径。

1. 核心科学问题，理论基础与技术细节

1.1 核心科学问题：尺度爆炸与非周期性灾难

激子（Exciton）作为半导体中受库仑力束缚的电子-空穴对，其物理性质由 Bethe-Salpeter Equation (BSE) 决定。在标准的周期性晶体中，动量空间形式的 BSE 能够高效求解。然而，在以下两种系统中，动量守恒失效：

  1. 超级莫尔系统：由多个扭转角度或晶格失配层叠加而成，其特征长度尺度可比单一莫尔系统大几个数量级。
  2. 准晶系统：具有长程有序但缺乏平移对称性。

在这些系统中，必须使用实空间 BSE。如果系统有 $N$ 个格点，激子基组 $|l, m\rangle$（表示空穴在 $l$ 位，电子在 $m$ 位）的大小为 $N^2$。当 $N = 10^9$ 时，激子 Hamiltonian 矩阵大小为 $10^{18} \times 10^{18}$。即使是存储该矩阵的对角元，目前的超级计算机也无法做到。

1.2 理论基础：实空间激子 Hamiltonian

作者从一个两带的 Wannier 紧束缚模型出发，定义总 Hamiltonian 为：

$$\hat{H} = \hat{H}_c + \hat{H}_v + \hat{H}_U$$

通过投影到激子子空间，得到有效激子 Hamiltonian $\hat{H}_X$：

$$\hat{H}_X = \hat{T}_c \otimes \mathbb{I} - \mathbb{I} \otimes \hat{T}_v - \hat{U}$$

其中 $\hat{T}_c$ 和 $\hat{T}_v$ 分别包含电子和空穴的跳跃（Hopping）项及单粒子势能。$\hat{U}$ 是电子-空穴相互作用。关键难点在于 $\hat{U}$ 项，它耦合了电子和空穴的坐标，通常在实空间表现为 $U_{klmn} = U_k \delta_{kl}\delta_{mn}\delta_{ln}$，即局部相互作用仅当电子和空穴处于同一格点时才存在。

1.3 技术细节：Quantics 张量网络编码

该方法的核心是 Quantics Tensor Network (QTN) 编码。它将坐标 $x \in [0, 2^L-1]$ 映射为 $L$ 个二进制位，从而将一个大维度的向量转化为一个具有 $L$ 个物理索引的张量。对于 $10^9$ 个格点，$L \approx 30$。这意味着系统的大小仅以 $\log N$ 的方式增长。

关键创新：交错排序（Interleaved Ordering） 在传统的张量网络表示中，如果先排列所有电子位再排列所有空穴位，电子-空穴相互作用算符（MPO）的键维度（Bond Dimension）会随系统规模呈指数级增长。作者提出将电子和空穴的二进制位交错排列：$(s_{e,1}, s_{h,1}, s_{e,2}, s_{h,2}, \dots, s_{e,L}, s_{h,L})$。在这种表示下：

  • 局部相互作用算符 $\hat{U}$ 变为张量网络上的局域算符。
  • 算符的键维度 $\chi$ 保持为常数（$\chi=2$），极大地降低了计算复杂度。

1.4 技术难点：高分辨率光谱解析

计算激子局域态密度（LDOS） $\rho(\mathbf{r}, E) = \langle \mathbf{r} | \delta(E - \hat{H}_X) | \mathbf{r} \rangle$ 需要高效的谱方法。作者采用了 核多项式方法 (KPM)，利用 Chebyshev 多项式展开 $\delta$ 函数。然而，对于 2D 系统，Chebyshev 矩的增长会导致张量压缩变得困难。为此，他们引入了 高阶 Delta-Chebyshev (HODC) 核，通过有理逼近技术，在较低的展开阶数下获得了极高的能量分辨率，并有效地抑制了 Gibbs 振荡。

2. 关键 Benchmark 体系与计算数据

2.1 1D 不可通约超级莫尔势能系统

系统规模：$N = 2^{30}$ 个格点（约 10.7 亿）。 模型参数：势能函数包含两个不可通约的余弦项，分别模拟原子尺度调制和大尺度超级莫尔调制。参数设为 $U=5t, V_0=3t$。 计算结果：

  • 激子能带分裂：在 LDOS 图谱中清晰地观察到了由小尺度调制引起的“微带（Miniband）”分裂。
  • 空间局域化：激子 LDOS 完美遵循了势能阱的轮廓。即使在十亿格点的尺度上，原子级的空间分辨率依然得以保持。
  • 性能对比：传统对角化方法无法处理超过 $10^5$ 维度的矩阵，而该方法处理了 $10^{18}$ 维度的有效空间。

2.2 2D 准晶超级莫尔系统

系统规模：$N_x = N_y = 2^{15}$，总格点数 $N \approx 1.07 \times 10^9$。 物理背景：模拟具有八重对称性的扭转 GeX/SnX 莫尔超晶格。势能函数由四组波矢的叠加组成。 计算数据：

  • 激子受限模式：在能量 $E_X = -4t$ 处，激子的 LDOS 展现出明显的八重对称性。光谱权重集中在势能极小值点，形成了激子量子点阵列。
  • 多尺度分辨：计算不仅捕捉到了跨越数万个格点的宏观莫尔图案，还通过缩放（Zoom-in）展示了原子格点层次上的激子局域化细节。
  • 键维度表现：虽然 2D 系统中的算符键维度增长比 1D 快，但通过 HODC 核，在仅使用约 40 个 Chebyshev 矩的情况下，成功解析了关键光谱特征。

3. 代码实现细节与复现指南

3.1 核心软件包

该研究主要基于 Julia 语言生态系统，利用了张量网络计算的高性能库：

  • ITensor.jl: 用于处理张量操作、矩阵乘积态（MPS）和矩阵乘积算符（MPO）的核心库。
  • TensorCrossInterpolation.jl: 用于实现 Quantics Tensor Cross Interpolation (QTCI)，这是构建复杂空间势能 MPO 的关键。它能够自动找到函数的低阶张量表示。

3.2 关键算法实现流程

  1. Hamiltonian 构建：
    • 使用 QTCI 算法将实空间势能 $V(\mathbf{r})$ 转化为 MPO 形式。
    • 构建电子项 $T_c$ 和空穴项 $T_v$ 的 MPO。
    • 按照“交错排序”规则重组 MPO 核心张量。
  2. Chebyshev 迭代：
    • 初始化激子态 $|\nu_0\rangle = |\mathbf{r}, \mathbf{r}\rangle$（表示电子和空穴处于同一位置）。
    • 执行迭代：$|\nu_1\rangle = \hat{H}_X |\nu_0\rangle$，$|\nu_n\rangle = 2\hat{H}_X |\nu_{n-1}\rangle - |\nu_{n-2}\rangle$。
    • 注意：每一步迭代后需进行张量截断（Truncation），控制奇异值舍弃精度以平衡内存与准确性。
  3. HODC 核处理：
    • 根据论文附录中的公式 (4)-(9)，通过求解 Vandermonde 线性方程组获取权重 $w_l$ 和极点 $z_l$。
    • 应用离散余弦变换（DCT）获取修正后的 Chebyshev 系数 $\nu_k$。

3.3 开源仓库链接

作者已将复现代码开源（基于 Julia）： https://github.com/mousanouar/Tensor-Network-Approach-to-Moire-Excitons

4. 关键引用文献与局限性评论

4.1 关键引用文献

  1. White (1992) [32]: 密度矩阵重整化群（DMRG）的奠基之作，张量网络方法的物理根源。
  2. Oseledets (2011) [45]: 提出了 Tensor Train 分解，即 QTCI 所依赖的数学框架。
  3. Weiße et al. (2006) [67]: 核多项式方法（KPM）的标准指南。
  4. Yi et al. (2025) [69]: 本文使用的高阶 Delta-Chebyshev (HODC) 核的原始理论来源。

4.2 局限性评论

尽管该工作在格点规模上取得了前所未有的突破，但仍存在以下局限：

  1. 2D 压缩极限：在 2D 系统中，随着 Chebyshev 迭代阶数的增加，纠缠熵迅速增长，导致 MPS/MPO 的键维度爆炸。目前仅能稳定计算约 40 个矩，这限制了在极窄带系统中的能量分辨率。
  2. 非正定性风险：HODC 核虽然收敛快，但由于其数学特性，可能会产生微小的负密度伪影（Negative Artifacts），在解释物理数据时需格外小心，可能需要手动“清理”数据。
  3. 相互作用范围：当前演示主要聚焦于局部库仑作用（On-site U）。如果引入长程库仑相互作用，MPO 的键维度将显著增加，这对算法的压缩效率提出了更高挑战。
  4. 算符局域化假设：该方法的成功高度依赖于电子-空穴相互作用在交错表示下的局域性，对于某些高度非局域的激子效应（如强极化激元），其优势可能减弱。

5. 补充内容：从 Wannier 到实空间 BSE 的推导解析

为了帮助量子化学背景的读者理解，我们补充一下论文附录中关于有效激子 Hamiltonian 的推导逻辑。

5.1 激子子空间的投影

定义多体基态为填满的价带（Fermi sea）$|\Omega\rangle = \prod_{l=1}^N c_{v,l}^\dagger |0\rangle$。激子态定义为从价带 $l$ 位激发一个电子到导带 $m$ 位：$|l, m\rangle = c_{v,l} c_{c,m}^\dagger |Omega\rangle$。

计算动力学项 $H_T$ 的矩阵元：

$$H_{lk, l'k'}^T = \langle l, k | \hat{H}_T | l', k' \rangle$$

利用费米子对易关系 $ \{c_i, c_j^\dagger\} = \delta_{ij} $，经过繁琐的对易子展开，可以证明：

$$H_{lk, l'k'}^T = t_{kk'}^c \delta_{ll'} - t_{l'l}^v \delta_{kk'} + \left( \sum_i t_{ii}^v \right) \delta_{kk'} \delta_{ll'}$$

这里的第一项是电子跳跃，第二项是空穴跳跃（注意索引交换和负号，反映了空穴的电荷特性），第三项是抵消项，通常作为参考能量零点减去。

5.2 相互作用项的演化

对于局部相互作用 $\hat{H}_U = \sum_i U_i c_{c,i}^\dagger c_{c,i} c_{v,i} c_{v,i}^\dagger$（注意空穴表示为 $c_v$ 的缺失），其矩阵元为：

$$H_{lk, l'k'}^U = U_k \delta_{kk'} \delta_{ll'} - U_k \delta_{kk'} \delta_{ll'} \delta_{kl}$$

这清晰地显示了激子是如何被“驱动”形成的：第二项仅在电子和空穴处于同一格点（$k=l$）时生效，且符号与势能项相反，表现为吸引力。正是这一项在超级莫尔势阱中捕获了电子-空穴对。

5.3 未来展望

这种张量网络方法不仅限于激子。它提供了一个通用的、对数级缩放的框架，可以直接应用于：

  • 双极化子 (Bipolarons)：两个电荷与格点振动的耦合。
  • 双激子 (Biexcitons)：莫尔系统中更复杂的多体关联。
  • 超导配对动力学：在非周期性实空间系统中的 Cooper 对传播。

随着张量网络算法（如树状张量网络 Tree TN）的进一步优化，我们有望在笔记本电脑上完成原本需要超级计算机才能处理的凝聚态多体模拟。