[论文解读] Inversion symmetry of singular values and a new orbital ordering method in tensor train approximations for quantum chemistry

Mi-Song Dupuy, Gero Friesecke|arXiv (Cornell University)|Feb 18, 2020

Tensor decomposition and applications参考文献 28被引用 6

一句话总结

本文提出了一种用于量子化学张量轨道（TT）近似的新轨道排序方法，该方法利用了非相互作用哈密顿量基态奇异值中意外出现的反演对称性。通过优化由此对称性导出的隐藏不变量，该方法将奇异值尾部衰减降低了4至5个数量级，显著优于现有的Fiedler排序等方法。

ABSTRACT

The tensor train approximation of electronic wave functions lies at the core of the QC-DMRG (Quantum Chemistry Density Matrix Renormalization Group) method, a recent state-of-the-art method for numerically solving the $N$-electron Schr\\"odinger equation. It is well known that the accuracy of TT approximations is governed by the tail of the associated singular values, which in turn strongly depends on the ordering of the one-body basis. Here we find that the singular values $s_1\\ge s_2\\ge ... \\ge s_d$ of tensors representing ground states of noninteracting Hamiltonians possess a surprising inversion symmetry, $s_1s_d=s_2s_{d-1}=s_3s_{d-2}=...$, thus reducing the tail behaviour to a single hidden invariant, which moreover depends explicitly on the ordering of the basis. For correlated wavefunctions, we find that the tail is upper bounded by a suitable superposition of the invariants. Optimizing the invariants or their superposition thus provides a new ordering scheme for QC-DMRG. Numerical tests on simple examples, i.e. linear combinations of a few Slater determinants, show that the new scheme reduces the tail of the singular values by several orders of magnitudes over existing methods, including the widely used Fiedler order.

研究动机与目标

• 解决张量轨道（TT）近似在量子化学中基底排序的关键挑战，该挑战决定了QC-DMRG方法的精度。
• 研究非相互作用哈密顿量中奇异值在TT近似下的行为，特别是其尾部衰减如何依赖于轨道排序。
• 开发一种理论基础坚实的优化轨道排序方案，以最小化关联波函数的奇异值尾部衰减。
• 提供一种实用且经过数值验证的方法，其在奇异值衰减和TT近似质量方面优于现有方法（如Fiedler排序）。

提出的方法

• 识别非相互作用基态张量表示中奇异值的新型反演对称性：$ s_1 s_d = s_2 s_{d-1} = \text{const} $，其中 $ s_1 \to s_d $ 为按降序排列的奇异值。
• 定义一个隐藏不变量 $ I = s_1 s_d $，其可完全表征给定轨道排序下奇异值尾部的行为。
• 提出“最优前因子排序”——一种通过最小化该不变量以抑制奇异值衰减的轨道重排策略。
• 通过利用构成斯莱特行列式成分的不变量的叠加，将该方法扩展至关联波函数，以对尾部衰减进行上界估计。
• 在少数斯莱特行列式的线性组合上实现并测试该方法，以模拟典型的量子化学态。
• 通过数值实验，比较该新方法与Fiedler排序及其他标准方案在奇异值衰减方面的表现。

实验结果

研究问题

• RQ1 在非相互作用系统中，张量轨道近似下奇异值尾部衰减对单体轨道排序的依赖关系如何？
• RQ2 非相互作用基态的奇异值谱是否表现出可被利用以改善近似的隐藏对称性？
• RQ3 能否利用来自斯莱特行列式成分的不变量，对关联波函数中奇异值的尾部行为进行上界估计？
• RQ4 基于最小化单一不变量的轨道排序方案，是否能显著优于现有方法（如Fiedler排序）的TT近似质量？

主要发现

• 非相互作用基态（斯莱特行列式）的张量轨道表示的奇异值表现出精确的反演对称性：$ s_1 s_d = s_2 s_{d-1} = \text{const} $，其明确依赖于轨道排序。
• 隐藏不变量 $ I = s_1 s_d $ 完全决定了奇异值的尾部行为，从而将尾部控制的复杂性简化为单一标量量。
• 通过“最优前因子排序”优化该不变量，与标准方法相比，奇异值尾部衰减降低了4至5个数量级。
• 该新方法优于Fiedler排序，后者仅在奇异值衰减方面实现约1个数量级的改善。
• 对于关联波函数，尾部衰减由各斯莱特行列式成分的不变量之和上界控制，从而支持推广的排序方案。
• 在少数斯莱特行列式的线性组合上的数值测试结果表明，该新方法在所有测试案例中均一致实现了更优的奇异值衰减。