Revisiting the Equation-of-Motion Method: A Universal Framework for Correlated Quantum Systems

技术摘要：重新审视关联量子系统的运动方程法

问题陈述
有限量子系统（尤其是原子核）中激发态的描述仍然是多体物理学中的一个核心挑战。虽然现代从头算（ab initio）方法在利用手征有效场论哈密顿量描述基态可观物理量方面取得了显著成功，但激发态（特别是涉及多粒子构型的相干叠加的集体模）的描述仍然十分困难。传统的运动方程（Equation-of-Motion, EOM）方法通常建立在独立粒子参考态（如 Hartree-Fock）之上，这限制了其一致性地将基态关联传递到激发态动力学过程的能力。此外，现有的 EOM 实现（在耦合簇 CC 和内中相似重整化群 IMSRG 框架内）往往依赖于相似变换后的哈密顿量以及截断的激发空间（例如 1p-1h 或 2p-2h），这实际上将问题简化为了配置相互作用（Configuration Interaction, CI）类型的特征值方程，可能遗失了与时间相关变分原理相关的本质动力学自由度。

方法论
本文提出了一种通用的 EOM 形式化实现，旨在直接作用于由现代多体求解器生成的关联参考态，而非无关联的平均场态。

• 理论框架： 推导过程始于时间相关（TD）变分原理（狄拉克作用量的极值条件），而非时间无关（TI）薛定谔方程。这种方法是在希尔伯特空间的子流形上处理能量函数 $H[\Psi]$。通过引入由反埃尔米特算符 $X(t)$ 生成的幺正变换，该方法推导出了能量函数的正规模。
• 算符基底： 一项关键创新是在完整的单体算符基底内构建进入 EOM 方程的对易子算符。不同于将空间限制在哈特里-福克真空上的粒子-空穴（ph）激发的传统 RPA，该公式利用了单体变换的完整李代数。这包括了粒子-粒子（$pp'$）和空穴-空穴（$hh'$）构型，这些构型对于当参考态$|\Psi\rangle$ 已具有关联性时进行一致性描述是必不可少的。
• 参考态： 该方法使用 IMSRG(2) 关联参考态进行应用。双对易子的矩阵元在这一关联态上进行求值。该方法避免了对哈密顿量进行相似变换的需求；相反，关联性被显式地编码在参考态中，并传播至激发算符。
• 数值实现： 使用全单体算符空间求解广义特征值问题。计算采用包含核-核（NN）和三体（3N）相互作用的手征核哈密顿量，并在谐振子基底中进行。该方法通过对 $^4$He 的无心壳模型（NCSM）结果进行了基准测试，并应用于氧同位素（$^{16}$O, $^{22}$O, $^{24}$O）。

主要贡献

1. 统一框架： 本工作建立了一个适用于任意关联波函数的通用 EOM 框架，使其脱离了传统的独立粒子图景。它提供了一种一致的处理方式，使得基态关联能够非微扰地转移到激发态。
2. 与 EOM-CC/IMSRG 的区别： 本文阐明了该方法与标准 EOM-CC 或 EOM-IMSRG 的结构性差异。虽然后者通常简化为在无关联参考态上对相似变换后哈密顿量的 CI 类对角化，但目前的公式直接在关联态上处理响应问题，保留了时间相关变分原理的辛结构。
3. 处理复数解： 研究识别并分析了在使用关联参考态求解 EOM 方程时出现的复能量解。论文认为，这些复特征值不仅仅是病态现象，更是诊断工具，表明参考态在全套单体变化（特别是超越 ph 区块）下并非能量泛函的稳定极值点。
4. 算符构造： 推导过程提供了一种独立于特定参考态的通用 EOM 方程算符结构，便于应用于各种多体求解器。

结果

• 基准测试： 对于 $^4$He 的同位素偶矢量偶极响应的计算显示，巨偶极共振（GDR）的位置（26 MeV）和高度具有快速收敛性，与 NCSM 基准测试吻合良好。
• 氧同位素： 在 $^{16}$O 中，该方法证明了包含基态关联可以显著减少不同手征相互作用（如 NNLOsat, N3LO Texas, 1.8/2.0 EM）所预测的 GDR 位置的弥散。在关联 EOM 框架下，峰值位置的分布从标准 RPA 的 ~7 MeV 降低到了 ~4 MeV（排除异常值后为 ~1.7 MeV）。
• 富中子系统： 对 $^{22}$O 和 $^{24}$O 的计算揭示了 GDR 向低能方向的移动与中子过剩相关，并伴随着在低能区的响应碎片化增加，这符合对富中子原子核的定性预期。
• 复谱： 包含全单体空间导致了复特征值的出现。论文观察到，复强度分布随所用相互作用的不同而变化，这表明存在与费米面以上激发态的耦合。

意义与主张
论文声称，EOM 框架为超越传统平均场范式的核谱学微观描述提供了一个灵活且强大的工具。其主要意义在于，在 ab initio 背景下，为基态关联与激发态动力学之间建立了连贯的、非微扰的联系。

作者强调，复解的出现应被视为参考态缺失关联的诊断，而非公式本身的失效。这一见解为未来的发展铺平了道路，例如轨道优化技术（针对单体旋转最小化能量）和自洽 RPA 方案。该框架被呈现为一个通用的工具，适用于任何能够评估二体算符期望值的计算方法，并具有扩展到对称性破缺参考态、投影平均场形式以及混合量子-经典算法的潜力。本工作并不声称解决了所有的碎片化问题（这需要更大的算符基底），但为系统性改进集体激发描述奠定了坚实的基础。