来源论文: https://arxiv.org/abs/2603.02011v2 生成时间: Mar 11, 2026 15:12

超越十亿格点：基于张量网络技术的超莫尔激子谱计算深度解析

0. 执行摘要

在凝聚态物理和量子化学的交叉领域，理解准晶（Quasicrystals）和超莫尔（Super-moiré）系统中的激发态物理——特别是激子（Excitons）的行为，一直是理论计算的“天花板”。由于激子体系的希尔伯特空间随格点数 $N$ 呈平方级（$N^2$）增长，传统的 Bethe-Salpeter 方程（BSE）求解器在处理超过数万个格点时便会遭遇“内存爆炸”。

由 Anouar Moustaj 和 Jose L. Lado 领衔的最新研究提出了一种革命性的张量网络（Tensor-network, TN）方法。该方法巧妙地将实空间 BSE 哈密顿量编码为矩阵乘积算符（MPO），并利用交错排序（Interleaved Ordering）策略将算符的键维（Bond Dimension）降至常数级。结合高阶 Delta-切比雪夫核（HODC），该算法成功模拟了包含超过 10 亿个格点的激子系统，计算规模相比传统方法提升了数个数量级。这一突破不仅为超莫尔量子物质的模拟建立了新的范式，也为大尺度量子化学相关领域的多体相关问题提供了高效的算符压缩思路。

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

1.1 核心科学问题：尺度鸿沟与维度灾难

在二维范德华材料（如 TMDs）中，通过旋转扭角或晶格失配形成的莫尔纹（Moiré pattern）会产生一个有效势场，从而限制电子和空穴的运动。当两个莫尔纹相互叠加形成所谓的“超莫尔”结构时，其特征长度尺度可能达到微米级。这意味着有效的紧束缚模型需要包含数百万甚至数十亿个格点。

对于单粒子谱，现有的量子多体求解器（如基于 TN 的单粒子算符编码）已经能处理大尺度体系。然而，激子物理涉及电子-空穴对的相关性。在实空间 BSE 框架下，哈密顿量 $H_X$ 的维度是 $N^2 imes N^2$。对于一个拥有 $10^9$ 个格点的系统，哈密顿量矩阵的维度达到 $10^{18}$。即使是显式存储这个矩阵的万亿分之一也是不可能的。传统的动量空间方法在面对不通约（Incommensurate）的准晶或超莫尔系统时会因平移对称性缺失而失效，因此必须回到实空间寻找对数级缩放的方案。

1.2 理论基础：实空间 Bethe-Salpeter 方程

研究从一个极简的二带 Wannier 哈密顿量出发：

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

其中 $c$ 和 $v$ 分别代表导带和价带。通过将全哈密顿量投影到单电子-空穴激发子空间（即激子基底 $|l, m\rangle = c^{\dagger}_{c,l} c_{v,m} |\Omega\rangle$），得到激子哈密顿量：

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

这里的核心难点在于最后的一项 $\hat{U}$，它耦合了电子和空穴的坐标。如果采用普通的乘积形式，其在张量网络表示下的局部性非常差，会导致计算成本随系统尺寸指数增长。

1.3 技术难点：键维度的指数爆炸

在张量网络中，我们将实空间位置编码为伪自旋索引（Pseudo-spin indices）。如果简单地将所有电子位点排列在所有空穴位点之前，那么 $\hat{U}$ 项（描述电子和空穴在同一位置的相互作用）将涉及到链条两端极远距离的算符耦合。在张量收缩理论中，这意味着中间切迹的张量键维 $\chi$ 会随着系统层数 $L$ 指数级增长 $\chi \sim 2^L$。

1.4 方法细节：三大核心技术攻关

1.4.1 交错排序 MPO 构造 (Interleaved Ordering)

这是本项目最重要的工程突破。作者没有采用传统的“先电子后空穴”排序，而是采用了交错排序：$(s_{e,1}, s_{h,1}, s_{e,2}, s_{h,2}, \dots, s_{e,L}, s_{h,L})$。在这种排序下，物理上对应的电子和空穴位点在张量链中变成了邻居。通过这种重排，局域库仑相互作用 $\hat{U}$ 变成了一个纯局域算符。其结果是，无论系统规模多大（即 $L$ 有多大），哈密顿量的 MPO 键维始终保持为常数（$\chi = 2$）。这一改进直接将原本不可计算的问题转化为了线性复杂度问题。

1.4.2 Quantics 张量交叉插值 (QTCI)

为了高效地在 MPO 中编码空间变化的势场（如莫尔势场），作者使用了 Quantics 算符表示技术。利用 QuanticsTCI.jl 库，可以将复杂的解析函数（如公式 4 和 5 中的多项余弦叠加）以对数级精度自动生成 MPO 核心张量。这使得模拟能够同时捕捉原子级的微观结构和微米级的宏观调制。

1.4.3 高阶 Delta-切比雪夫核 (HODC)

传统的切比雪夫核多项式方法（KPM）使用 Jackson 核来抑制 Gibbs 振荡，但其代价是分辨率受限于 $1/N_\mu$。在处理二维激子谱时，由于键维增长迅速，迭代次数 $N_\mu$ 受限。作者引入了 HODC 技术，通过有理函数近似 $\delta(\omega - H)$，实现了比 Jackson 核快得多的收敛速度（收敛阶数为 $\eta^m$）。这使得在仅迭代 40 次的情况下，就能在 2D 系统中清晰地分辨出激子局部态密度（LDOS）的精细特征。

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

2.1 1D 不通约超莫尔势场

体系描述：一个 1D 链条，其格点数为 $N = 2^{30}$（约 10 亿个格点）。势场由两个不同周期的余弦函数叠加而成： $V_i = V_0 [1 + 0.1 \cos(\pi x_i / 3\sqrt{5}) + 0.1 \cos(5\pi x_i / N\sqrt{3})]$

计算结果：

  • 激子能带形成：LDOS 图像显示，孤立的激子带完美地跟随了宏观莫尔势场的空间轮廓。
  • 微型能带（Minibands）：通过放大局部图像（长度 $\ell=64$），可以看到由于小尺度势场调制导致的能带分裂（Miniband splitting）。
  • 多尺度分辨能力：算法成功展示了在十亿位点背景下，原子级（$x=0$ 附近）的激子局域化行为。

2.2 2D 准晶超莫尔系统

体系描述：模拟一个 2D 方格子，格点数为 $N_x = N_y = 2^{15}$（总计约 $10^9$ 格点）。电势场模拟了旋转 45 度的双层准晶势场，具有八重旋转对称性。

数据表现：

  • 激子限域模式：在束缚激子能量 $E_X = -4t$ 处，LDOS 显示激子在势场极小值点呈现出清晰的八重对称限域图案。
  • 性能对比：传统 BSE 求解器对 $10^{18}$ 维度的矩阵完全束手无策，而该算法在普通高性能计算节点上通过控制键维实现了全实空间分辨率。
  • HODC 增益：相比 Jackson 核，HODC 在 $N_\mu = 40$ 时就消除了绝大部分模糊，清晰地揭示了激子在莫尔陷阱中的原子级分布。

3.1 核心软件包

该研究的实现高度依赖于 Julia 语言生态系统：

  1. ITensor / ITensors.jl：用于基础的张量收缩、MPO 构造和 MPS 操作。ITensor 提供的自动收缩和索引管理是处理复杂交错排序的基础。
  2. QuanticsTCI.jl：这是由 M. Ritter 等人开发的库，用于执行 Quantics Tensor Cross Interpolation。它能将实空间函数 $V(r)$ 自动化地转化为张量训练（Tensor Train）格式。
  3. KernelPolynomials.jl（或自定义 KPM 模块）：用于实现切比雪夫递归逻辑。特别需要注意递归中的算符作用：$|\nu_n\rangle = 2\hat{H}|\nu_{n-1}\rangle - |\nu_{n-2}\rangle$。

3.2 代码复现逻辑

复现该工作的关键步骤如下：

  • 步骤 1：坐标编码。将 $2^L$ 个位点映射到 $L$ 个伪自旋位点。对于 2D 体系，采用 Z-ordering（Lebesgue 曲线）来保持局部性。
  • 步骤 2：构造单粒子 MPO。分别构造 $T_e$ 和 $T_h$。使用 expansion 函数配合 $\delta$ 张量（恒等映射）将它们扩展到 $2L$ 位点的联合空间。
  • 步骤 3：交错重排。通过索引置换将位点顺序改为 $(e, h, e, h, \dots)$。这是降低键维的关键代码段。
  • 步骤 4：势场插值。调用 QTCI 算法，输入解析势场函数，生成对应的 MPO。
  • 步骤 5：KPM 迭代。选择初始状态为局域化的 $|r, r\rangle$，进行多项式迭代。如果是 2D，务必使用 HODC 核以节省内存并获得高分辨率。

3.3 开源资源

作者已将核心代码库开源：

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

4.1 关键引用文献

  • [31] Sun et al. (2025)：奠定了利用 TN 处理超十亿格点单粒子问题的基础。
  • [48, 50] Ritter et al. (2024, 2025)：QTCI 技术的核心参考文献，提供了高分辨率函数编码的数学工具。
  • [69] Weiße et al. (2006)：经典的 KPM 方法综述，理解 Jackson 核的起点。
  • [71] Yi et al. (2025)：HODC 高阶切比雪夫核的来源，是 2D 高精度模拟的关键。

4.2 工作局限性评论

尽管该工作在规模上取得了巨大突破，但从量子化学和材料模拟的严谨角度看，仍存在以下局限：

  1. 短程相互作用假设：目前的模型仅考虑了局域（on-site）库仑作用。虽然论文提到可以扩展到长程库仑，但在 TN 框架下，长程作用（$\sim 1/r$）会显著增加 MPO 的键维，可能导致在 2D 系统中出现计算瓶颈。
  2. 交换项缺失：计算中忽略了激子交换相互作用（Exchange terms）。对于某些细致的能带结构（如谷分裂），交换项至关重要。
  3. 非正定性风险：HODC 核虽然分辨率高，但不保证 LDOS 的非负性。在数据处理中需要手动将负值置零，这在能量分辨要求极高的物理分析中可能引入微小伪影。
  4. 静态格点局限：模型假设晶格是刚性的，未考虑莫尔晶格常见的机械弛豫（Atomic reconstruction）效应，这在真实莫尔材料中会显著改变势场分布。

5. 其他补充：量子化学视角下的广阔前景

5.1 对大分子激子动力学的启示

这项技术对于大尺度共轭聚合物或光合作用复合物（LHC）的模拟极具参考价值。在量子化学中，这些系统通常包含数万个原子，且激子在不同发色团之间的跳跃可以用类似的紧束缚模型描述。利用张量网络的对数级缩放，我们可以模拟超长聚合物链中的激子扩散和猝灭过程。

5.2 关联电子态的普适扩展

论文在 Discussion 部分提到，该框架不仅限于激子。实际上，任何二体相关问题——如**双极化子（Bipolarons）**的形成、**双镁振子（Bimagnons）**的相关性，都可以采用这种交错 MPO 方法进行编码。这为处理强关联电子系统中的复杂复合粒子提供了一条全新的数值路径。

5.3 张量交叉插值（TCI）的未来

QTCI 技术正在成为连接解析物理和数值计算的桥梁。在量子化学中，多维势能面（PES）的拟合一直是一个难题。如果能将 TCI 引入 PES 编码，配合张量收缩，或许能解决高维振动能级计算中的维度灾难问题。

总结

Anouar Moustaj 等人的这项工作不仅是计算规模的胜利，更是物理算法设计的胜利。通过对哈密顿量结构的深刻理解（交错排序），他们将一个看似不可逾越的内存障碍化解为优雅的张量操作。这标志着实空间激发态计算正式进入了“十亿格点”时代。