2025-11-16T18:52:11.890181

Towards Balanced Description of Ground and Excited States with Transcorrelated F12 Methods

Masteran, Gaudel, Valeev
By correlating only the 1-particle states occupied in the reference determinant the conventional design for the single-reference R12/F12 explicitly-correlated methods biases them towards the ground state description thereby making treatment of response properties of the ground state, and energies and other properties of excited states less robust. While the use of multireference methods and/or extensions of the standard SP projected geminals can achieve a more balanced description of ground and excited states, here we show that the same goals can be achieved by extending the action of F12 correlators to the occupied and valence unoccupied 1-particle states only. This design choice reflects the strong dependence of the optimal correlation lengthscale of the F12 ansatz on the orbital energies/structure, and helps to avoid the unphysical raising of the ground state energy if the F12 geminals are used to correlate pairs of all 1-particle states. The improved F12 geminal design is incorporated into the unitary transcorrelation framework to produce a unitary 2-body Hamiltonian that incorporates the short-range dynamical correlation physics for ground and low-energy excited states in a balanced manner. This explicitly-correlated effective Hamiltonian reduces the basis set requirement on the correlation-consistent basis cardinal number by 1 or more over the uncorrelated counterpart for the description of the ground state coupled-cluster singles and doubles (CCSD) energies, the vertical excitation energies and harmonic vibrational frequencies of equation-of-motion CCSD low-energy excited states.
academic

Towards Balanced Description of Ground and Excited States with Transcorrelated F12 Methods

基本信息

  • 论文ID: 2508.20233
  • 标题: Towards Balanced Description of Ground and Excited States with Transcorrelated F12 Methods
  • 作者: Conner Masteran, Bimal Gaudel, Edward F. Valeev (Virginia Tech)
  • 分类: physics.chem-ph
  • 发表时间: 2025年8月29日
  • 论文链接: https://arxiv.org/abs/2508.20233

摘要

传统的单参考R12/F12显式相关方法仅关联参考行列式中占据的单粒子态，这使得它们偏向于基态描述，从而使基态响应性质的处理以及激发态的能量和其他性质变得不够稳健。本文证明了通过将F12相关子的作用扩展到占据轨道和价层未占据单粒子态，可以实现基态和激发态的平衡描述。这种设计反映了F12 ansatz的最优相关长度尺度对轨道能量/结构的强烈依赖性，并有助于避免在F12 geminals用于关联所有单粒子态对时基态能量的非物理升高。改进的F12 geminal设计被整合到幺正跨相关框架中，产生了一个幺正二体哈密顿量，以平衡的方式包含基态和低能激发态的短程动力学相关物理。

研究背景与动机

问题背景

  1. 显式相关方法的重要性: F12方法通过在波函数中显式包含电子间距离来大幅加速基组收敛，有效建模电子-电子尖点附近波函数的已知解析行为。
  2. 传统F12方法的局限性:
    • 绝大多数F12方法的发展和应用集中在基态性质上
    • 传统的II ansatz（仅关联占据轨道对）对激发态描述存在偏差
    • 导致垂直激发能被人为高估，如Fliegl等人首次观察到的
  3. 现有解决方案的不足:
    • 多参考方法虽然可以平等处理多个电子态，但计算成本高
    • 扩展SP ansatz需要额外参数优化，可能导致数值不稳定

研究动机

本文的核心假设是：单参数指数相关因子主要针对价层轨道调优，不适合描述包含非价层轨道的轨道对。这一观点得到了核相关研究的支持，Werner等人发现核-核、核-价和价-价轨道对的最优相关因子指数存在显著差异。

核心贡献

  1. 提出了价层约束geminal方法: 将F12相关子限制在价层虚轨道子集，避免了传统方法中基态能量的非物理升高
  2. 开发了三种新的相关子ansatz:
    • II ansatz: 传统方法，仅关联占据轨道对
    • IP ansatz: 关联占据-虚轨道对的混合方法
    • PP ansatz: 关联所有活性轨道对的完全方法
  3. 实现了平衡的基态-激发态描述: VVO PP相关子在保持基态精度的同时显著改善激发态描述
  4. 证明了基组收敛性的显著改善: 相关能和激发能的基组基数需求减少1个或更多

方法详解

任务定义

本研究旨在开发能够平衡描述分子基态和激发态的F12跨相关方法，输入为分子几何结构，输出为高精度的基态和激发态能量及性质。

核心方法架构

1. 幺正F12跨相关框架

跨相关哈密顿量通过相似变换生成：

Ĥ̄ = e^(-Â)Ĥe^Â

其中Â是包含粒子间距离依赖性的算符。

2. 三种相关子ansatz

II Ansatz (传统方法):

Â_II = 1/2(G^(α₁α₂)_(i₁i₂) Ê^(i₁i₂)_(α₁α₂) - h.c.)

IP Ansatz (混合方法):

Â_IP = (G^(α₁α₂)_(i₁p'₁) Ê^(i₁p'₁)_(α₁α₂) - h.c.) - Â_II

PP Ansatz (完全方法):

Â_PP = 1/2(G^(α₁α₂)_(p'₁p'₂) Ê^(p'₁p'₂)_(α₁α₂) - h.c.)

3. 价层虚轨道(VVO)构建

通过最小基组投影定义价层虚轨道：

  1. 对最小基组AO进行Löwdin对称正交化
  2. 计算与未占据轨道的重叠矩阵
  3. 通过SVD分解获得价层虚轨道

技术创新点

  1. 价层约束策略: 首次将F12 geminals限制在价层轨道子空间，基于轨道能量尺度的物理考虑
  2. 自动化实现: 使用SeQuant符号张量代数包自动生成和优化张量表达式
  3. 状态普适性设计: PP ansatz提供了更好的低能激发态描述，同时保持基态精度

实验设置

测试分子体系

选择了5个小分子进行详细研究：BH、CH₂O、N₂、H₂O、CH₂，涵盖不同类型的激发态特征。

基组选择

  • aug-cc-pVXZ (aXZ)
  • doubly-augmented daug-cc-pVXZ (daXZ)
  • cc-pVXZ-F12 (XZ-F12)
  • 使用推荐的geminal指数值

计算方法

  • 基态：完全重新优化Hartree-Fock轨道 + CCSD
  • 激发态：冻芯EOM-CCSD方法
  • 几何优化和振动频率通过拟合势能面获得

实验结果

主要结果

1. 垂直激发能改善

以BH分子为例（表1）：

  • II相关子使激发能高估约0.1 eV
  • VVO PP相关子在aTZ水平上将激发能误差从0.01 eV降至0.004 eV
  • 基态相关能恢复接近II ansatz水平

2. 不同分子体系的表现

  • CH₂O: VVO PP相关子显著改善Rydberg态激发能，DZ水平误差从0.2-0.3 eV降至更小值
  • N₂: II相关子高估激发能达0.5 eV，VVO IP/PP相关子显著减小误差
  • H₂O: 对于从下方收敛的激发能，F12跨相关显著改善基组收敛性

3. 激发态性质改善

键长和振动频率的基组收敛性平均提升相当于1个基数：

  • 平衡键长：F12方法与传统EOM-CCSD(F12)结果几乎一致
  • 振动频率：VVO PP相关子的{平均, 最大}绝对误差为{2.1, 5} cm⁻¹

消融实验

通过FVO vs VVO对比证明了价层约束的重要性：

  • FVO变体普遍产生较小的基态相关能
  • VVO变体在基态和激发态之间实现更好平衡

相关工作

本文建立在以下重要工作基础上：

  1. F12理论基础: Ten-no的SP ansatz和标准投影geminal
  2. 激发态F12方法: Köhn的XSP ansatz，Bokhan和Ten-no的EOM-CC-F12
  3. 跨相关方法: Yanai-Shiozaki的幺正跨相关框架
  4. 价层轨道理论: Lu等人和Schmidt等人的价层虚轨道定义

结论与讨论

主要结论

  1. VVO PP相关子最优: 在基态相关能恢复和激发态描述之间实现最佳平衡
  2. 基组收敛显著改善: 相关能基组收敛性提升相当于2个基数
  3. 方法的普适性: 适用于不同类型的激发态（价层态和Rydberg态）

局限性

  1. 单粒子基组不完整性: 某些激发态仍存在显著的1粒子基组误差
  2. geminal参数优化: 当前使用固定参数，可能不是所有体系的最优选择
  3. 体系规模限制: 目前仅在小分子体系上验证

未来方向

  1. 1粒子基组不完整性处理: 开发状态普适的单粒子修正方法
  2. 相对论效应扩展: 将方法扩展到包含特殊相对论效应
  3. 响应性质评估: 评估方法对分子响应性质的表现

深度评价

优点

  1. 理论创新性强: 首次系统性地将F12方法扩展到平衡描述基态和激发态
  2. 物理直觉清晰: 基于轨道能量尺度的价层约束策略具有清晰的物理动机
  3. 实验充分: 涵盖多种分子体系和激发态类型，结果具有说服力
  4. 实现完整: 提供了完整的自动化实现框架

不足

  1. 测试体系有限: 仅在小分子上测试，大分子的表现未知
  2. 理论分析不足: 缺乏对为什么VVO约束有效的深入理论分析
  3. 计算成本分析: 未详细讨论相对于传统方法的计算开销

影响力

  1. 方法学贡献: 为激发态F12方法提供了新的设计范式
  2. 实用价值: 可直接用于改善激发态计算精度
  3. 可复现性: 提供了详细的技术实现细节

适用场景

  1. 光化学研究: 需要高精度激发态能量和性质的研究
  2. 光谱学计算: 垂直激发能和振动频率的精确预测
  3. 方法学开发: 作为进一步开发平衡基态-激发态方法的基础

参考文献

本文引用了91篇重要文献，涵盖了F12理论、跨相关方法、激发态计算等多个领域的经典工作，为研究提供了坚实的理论基础。

总结: 这是一篇在量子化学方法学领域具有重要贡献的论文，通过巧妙的价层约束策略解决了F12方法在激发态描述中的长期问题，为平衡描述基态和激发态提供了新的解决方案。方法具有清晰的物理动机和良好的数值表现，在激发态量子化学计算领域具有重要的应用前景。