⚛️ quantum physics

Site Basis Excitation Ansatz for Matrix Product States

本文提出了一种基于单点基的激发变分法（SBEA），通过利用无限矩阵乘积态描述基态并结合非正交基组对角化，高效且高精度地计算了一维量子格点系统的元激发谱，同时揭示了保持基组非正交性对收敛性的关键作用及构建沃尼耶激发的方法。

原作者： Steven R. White

发布于 2026-03-03

📖 1 分钟阅读🧠 深度阅读

原作者： Steven R. White

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

这篇论文介绍了一种名为**“站点基激发 Ansatz" (SBEA)** 的新方法，用来研究一维量子系统（比如原子排成一条线的材料）中的“粒子”是如何运动的。

为了让你轻松理解，我们可以把这篇论文的核心思想想象成**“如何用最少的积木，拼出最完美的波浪”**。

1. 背景：我们在找什么？

想象你有一根很长的弹簧（这就是一维量子系统）。如果你拨动它一下，会产生一个波（这就叫“激发”或“磁振子”）。

科学家想知道： 这个波在不同频率下（动量 $k$ ）的能量是多少？这就叫“能谱”。
以前的困难： 以前计算这个，就像是要为每一个可能的频率单独去解一道极其复杂的数学题。如果频率有 500 种，你就得算 500 次，非常耗时且容易出错。

2. 核心创新：SBEA 是什么？

Steven R. White 提出的新方法（SBEA）就像是一个**“万能积木盒”**。

旧方法（EA）： 就像是你每次想造一个波浪，都要现场重新设计一块特殊的积木（激发张量 $B$ ），而且每次都要重新计算。
新方法（SBEA）： 作者发现，其实不需要每次都重新设计。他先在一个“小房间”里（单个格点）做了一次计算，找出最关键的几种基础积木（比如 7 种）。
- 这就好比：你不需要为每一首曲子都发明新的音符，你只需要准备好“Do, Re, Mi, Fa, Sol, La, Si"这 7 个基本音符。
- 一旦有了这 7 个基础积木，任何频率的波浪（激发），都可以通过混合这 7 个积木来完美表示。

比喻：
想象你要画出一幅画里所有的云彩形状。

旧方法： 每画一朵云，都要重新从泥土里捏出一块新的云团。
SBEA 方法： 你先捏好 7 种不同形状的“标准云团”（基础积木）。以后不管要画什么云，你只需要把这 7 种云团按不同比例混合一下就行了。这大大节省了时间！

3. 关键技巧：为什么要“不整齐”？

在数学计算中，通常要求所有的“积木”必须是正交的（就像互相垂直的坐标轴，互不干扰）。以前的方法都强制要求积木必须“正交”。

但作者发现了一个反直觉的秘诀：不要强迫积木“正交”！

比喻： 想象你在拼乐高。以前大家觉得，每块积木必须严丝合缝、互不重叠才能拼得稳。但作者发现，如果允许积木之间稍微有点重叠（非正交），反而能用更少的积木（比如从 100 块减少到 7 块）就拼出完美的形状。
虽然积木重叠会让计算稍微复杂一点点（就像处理重叠的阴影），但换来的是巨大的效率提升和更简单的操作。作者说，如果强行把积木摆正（施加“左规范”条件），反而会让计算变得极其困难，甚至算不出来。

4. 另一个亮点：制造“魏纳尔”积木

论文还提到了一个很酷的概念：魏纳尔激发（Wannier Excitations）。

比喻： 在物理学里，我们通常用“波”（像海浪一样扩散）来描述粒子。但魏纳尔函数就像是把波“打包”成一个个局域的小包（就像把海浪打包成一个个独立的水球）。
作者展示了如何把这些“水球”打包好。这些“水球”非常集中，不会到处乱跑。
有什么用？ 这就像把复杂的交响乐简化成了几个简单的乐器独奏。以后如果要研究两个粒子怎么打架（多粒子激发），就可以直接用这些打包好的“水球”来模拟，就像搭积木一样简单。

5. 总结：这方法好在哪？

快：以前算 500 个频率要算 500 次，现在只需要算一次“基础积木”，剩下的 500 次只是简单的“混合”计算（几秒钟搞定）。
准：在著名的海森堡自旋链模型中，它算出的结果几乎和“上帝视角”的精确解一样准。
简单： 它不需要复杂的无限系统初始化技巧，甚至可以用普通的有限系统计算来起步，门槛更低，更容易被其他科学家使用。

一句话总结：
这篇论文发明了一种**“用少量基础积木混合出任意量子波浪”的聪明办法，并且发现“允许积木稍微重叠”**反而能让计算变得更快、更准、更简单。这就像是用一套简单的调色盘，就能画出世间万变的色彩。

这是一份关于 Steven R. White 论文《Site Basis Excitation Ansatz for Matrix Product States》（矩阵乘积态的格点基激发 Ansatz）的详细技术总结。

1. 研究背景与问题 (Problem)

核心挑战：计算一维量子晶格系统的基态激发能谱（Elementary excitation spectra）对于理论理解和实验对比至关重要。然而，相比于静态性质，获取完整的色散关系（dispersion relation）在数值上更具挑战性。
现有方法的局限：
- 传统 DMRG：在有限系统中只能获得有限数量的激发态，难以构建完整的色散关系。
- 激发态 Ansatz (EA)：基于无限 MPS (iMPS) 的单模近似（Single-mode approximation）虽然高效，但存在以下技术难点：
  1. 计算成本高：传统的 EA 方法需要对每一个动量 $k$ 分别进行完整的对角化计算。
  2. 收敛性问题：变分均匀矩阵乘积态 (VUMPS) 方法在初始化或收敛性上可能遇到困难。
  3. 规范选择 (Gauge Choice)：以往工作通常强制施加左正交规范（left-orthonormal gauge），这严重阻碍了基组的收敛效率。
  4. 复杂性：原始公式涉及复杂的图示，且对近邻相互作用以外的情况处理繁琐。

2. 方法论 (Methodology)

作者提出了一种名为格点基激发 Ansatz (Site Basis Excitation Ansatz, SBEA) 的新方法，主要包含以下三个核心步骤：

A. 基于有限系统 DMRG 的无限 MPS 构建

思路：为了简化无限 MPS (iMPS) 的获取，作者提出了一种基于有限系统 DMRG 的简单替代方案，而非直接使用 iDMRG 或 VUMPS。
步骤：
1. 在长开链（ $N \gg \xi$ ， $\xi$ 为关联长度）上运行有限系统 DMRG。
2. 在链中心插入一个新格点，将正交中心置于该处。
3. 利用 Lanczos 对角化优化该新格点的张量，使其能量最小化。
4. 由于新张量连接左右正交半链，它自然形成了 $\dots \lambda \Gamma \lambda \Gamma \dots$ 形式的均匀 iMPS。
5. 通过标准的 Orus-Vidal 过程将其转换为规范形式（Canonical form），即 $\dots AAA \dots$ 形式。
优势：该方法极其简单，利用了有限系统 DMRG 的成熟性和可靠性，且计算成本几乎等同于一次有限系统 DMRG 计算。

B. 格点基激发 Ansatz (SBEA) 的核心构建

基本思想：不再针对每个 $k$ 单独优化激发张量 $B(k)$ ，而是构建一个小的格点基 (Site Basis) $\{B_\alpha\}$ 来同时表示所有动量的激发。
基组生成：
1. 在单个格点上，对有效哈密顿量进行 Lanczos 对角化。
2. 选取能量最低的 $N_\alpha$ 个本征态作为基组 $\{B_\alpha\}$ 。
3. 关键创新：不施加左正交规范条件。作者发现，保留基组的非正交性（Non-orthogonality）对于基组截断至关重要。如果施加左规范，低能态会消失，导致需要极大的基组才能收敛。
4. 预截断 (Pre-truncation)：在 Lanczos 对角化之前，利用 Schmidt 值对张量 $B$ 的环境进行截断（将键维 $\chi$ 截断为较小的 $\chi_{tr}$ ）。这大大降低了计算成本，同时保留了构建单磁子色散所需的关键物理态。

C. 非正交能带理论 (Nonorthogonal Band Theory)

求解过程：
1. 计算重叠核 (Overlap kernel) $O_{\alpha\alpha'}(j'-j)$ 和哈密顿量核 $H_{\alpha\alpha'}(j-j')$ 。由于系统有能隙，这些核随距离指数衰减，只需计算有限范围。
2. 对实空间核进行傅里叶变换，得到动量空间的矩阵 $\tilde{O}(k)$ 和 $\tilde{H}(k)$ 。
3. 对于任意动量 $k$ ，求解广义本征值问题：
  $\tilde{H}(k)c = E(k)\tilde{O}(k)c$
效率：一旦计算出 $O$ 和 $H$ 矩阵，求解任意 $k$ 点的能谱仅需对角化一个极小的矩阵（尺寸 $N_\alpha \times N_\alpha$ ），类似于非正交能带理论的计算，效率极高。

D. 沃尼尔激发 (Wannier Excitations)

作者展示了如何构建类似于能带理论中沃尼尔函数 (Wannier functions) 的沃尼尔激发。
通过投影算符将动量本征态 $|k\rangle$ 变换回实空间，构造出定域化 (Localized) 且正交的激发态 $|W_j\rangle$ 。
这些定域态可以精确重构所有动量的单磁子模式，为构建多磁子有效哈密顿量提供了基础。

3. 关键结果 (Results)

测试系统：自旋 $S=1$ 的海森堡反铁磁链（Haldane 链）。
色散关系精度：
- 使用极小的基组（例如 $N_\alpha = 7$ ，对应预截断 $\chi_{tr}=8$ ），SBEA 就能极其精确地复现整个单磁子色散曲线。
- 计算得到的 Haldane 能隙 $\Delta/J = 0.4107$ ，与精确值 $0.410479\dots$ 高度吻合。
- 随着 $N_\alpha$ 的增加，色散曲线在更高能量处依然保持收敛。
规范选择的影响：
- 实验证明，不施加左规范是方法成功的关键。施加左规范会导致低能本征态能量虚高，使得基组截断失效，收敛极慢。
- 非正交基组虽然引入了重叠矩阵，但允许使用极小的 $N_\alpha$ 即可达到高精度。
沃尼尔激发：成功构造了定域化的沃尼尔激发态，其自旋分布范围在几个关联长度内，证明了该方法的物理自洽性。

4. 主要贡献 (Key Contributions)

简化无限 MPS 获取：提出了一种基于有限系统 DMRG 和单格点 Lanczos 的简单方法来构建 iMPS，避免了 VUMPS 的复杂性和收敛困难。
SBEA 算法：引入格点基激发 Ansatz，将每个动量的独立优化问题转化为一次性的基组构建和后续的小矩阵对角化，极大地提高了计算效率。
非正交基组的必要性：明确指出并验证了不施加左正交规范对于激发态计算的重要性。保留非正交性允许对基组进行强力截断，这是该方法高效的核心。
沃尼尔激发的构建：在 MPS 框架下成功构建了定域化的沃尼尔激发，为研究多粒子激发和构建有效模型提供了新工具。
普适性：该方法框架易于推广到更大的原胞（Unit cells）和更复杂的系统（如二维系统）。

5. 意义与影响 (Significance)

计算效率：SBEA 将激发态计算的成本从“每个 $k$ 点一次完整对角化”降低为“一次基组构建 + 任意 $k$ 点的小矩阵对角化”，使得扫描整个布里渊区变得极其快速（在普通笔记本电脑上计算 500 个 $k$ 点仅需数秒）。
易用性：该方法基于成熟的有限系统 DMRG，降低了技术门槛，使得非专家也能轻松应用。
理论深度：通过引入沃尼尔激发的概念，将凝聚态物理中的能带理论思想成功移植到张量网络激发态计算中，为理解多体激发提供了新的视角。
未来应用：该方法特别适用于具有有限关联长度的有能隙系统，有望成为研究二维量子系统激发谱的有力工具，填补了当前计算能力的空白。

总结来说，这篇论文通过巧妙的数学构造（非正交基组、格点基截断）和算法简化（基于有限 DMRG 的 iMPS 构建），显著提升了 MPS 计算激发能谱的效率和精度，是一篇在数值量子多体物理领域具有重要技术价值的论文。