来源论文: https://arxiv.org/abs/2607.08178v1 生成时间: Jul 10, 2026 00:23

0. 执行摘要

传统的密度矩阵嵌入理论 (DMET) 在处理局部激发态时，由于其嵌入空间通常基于基态Slater行列式构建，因此对基态和激发态的描述存在固有的不平衡性，往往高估激发能。本文作者Zhe-Bin Guan和Hong Jiang*（北京大学）提出了一种创新的态平均密度矩阵嵌入理论 (SA-DMET)，通过将DMET的低级参考点从单态Hartree-Fock扩展到态平均完全活性空间自洽场 (SA-CASSCF) 波函数，从而在构建浴轨道时显式地纳入了激发态信息。此外，研究还探索了组态平均Hartree-Fock (CAHF) 作为DMET起始点的性能，发现其在特定条件下等效于SA-CASSCF，并具备更高的计算效率。为进一步提升精度，本文还测试了非正交原子轨道DMET (AO-DMET) 的性能，该方法通过直接使用原子轨道定义杂质，避免了正交化引入的配体污染，更真实地保留了活性轨道的金属中心特征。研究通过对3d/4f单离子磁体、Ce(III)发光配合物的4f-5d跃迁以及多种过渡金属配合物的自旋态能量学计算，系统性地评估了SA-DMET结合CASSCF和NEVPT2方法的性能。结果表明，SA-DMET在所有测试体系中均显著优于传统的单态DMET，而与AO-DMET和CAHF的结合则能进一步将平均绝对误差降至极低水平，验证了这一新范式在强关联金属配合物局部激发态研究中的强大潜力与广阔应用前景。

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

1.1 核心科学问题

在量子化学领域，准确描述过渡金属和镧系配合物等强关联体系的电子激发态和磁性性质一直是一项重大挑战。这些体系通常涉及部分填充的d或f壳层，其电子结构表现出复杂的静态关联和动态关联效应。现有的从头算方法，如多组态自洽场 (MCSCF) 和多参考微扰理论 (MRPT)，虽然能够提供高精度结果，但其计算成本随体系规模呈指数增长，对于含有庞大配体基团的复杂配合物而言，往往难以承受。量子嵌入方法应运而生，通过将整个体系划分为一个小而关键的“杂质”区域和一个较大的“环境”区域，并对杂质区域应用高精度方法，同时以低级方法处理环境，从而实现精度与效率的平衡。

然而，现有的密度矩阵嵌入理论 (DMET) 在处理激发态，特别是局部激发态时，面临一个核心问题：其嵌入子空间通常基于低级基态波函数（如Hartree-Fock Slater行列式）构建，导致所生成的浴轨道固有的偏向于基态。这种偏置使得DMET对激发态的描述不平衡，常常高估激发能，甚至在某些情况下给出错误的基态自旋多重度预测，从而限制了其在自旋态能量学、磁各向异性以及光物理性质等研究中的应用。因此，如何构建一个能够平衡描述基态和激发态的嵌入子空间，成为DMET方法发展的一个关键瓶颈。

1.2 理论基础与DMET基本原理

DMET是由Chan及其合作者提出的一种强耦合量子嵌入理论，其核心思想在于利用一个低级、全体系计算的单电子约化密度矩阵 (1-RDM) 来定义一个高精度、小规模的“嵌入杂质哈密顿量”。该理论的数学严谨性在于通过施密特分解 (Schmidt decomposition) 的概念，将杂质区域和环境区域的轨道进行缠绕，从而提取出与杂质相互作用的“浴轨道” (bath orbitals)。

DMET的步骤如下：

  1. 体系划分：将整个体系的局域正交轨道 (LOs) 划分为杂质区域 (I) 和环境区域 (E)。
  2. 低级参考计算：对整个体系进行一个低级计算（如Hartree-Fock或ROHF），得到全体系的基态波函数 $\Phi_0$ 或1-RDM。
  3. 浴轨道构建：根据低级计算的1-RDM，在环境区域E中构建与杂质I缠绕的浴轨道。这些浴轨道与杂质轨道共同构成了“嵌入杂质空间” ($I_{emb}$)。
  4. 嵌入杂质哈密顿量：将全体系哈密顿量投影到$I_{emb}$张成的Fock空间，构建一个有效杂质哈密顿量 ($H_{emb}$)，其中包含单粒子项 ($h$) 和双电子相互作用项。单粒子项 $h$ (Eq. 2) 包括动能、核势以及核轨道a的库仑和交换算符 $J_a, K_a$。
  5. 高精度求解：使用高水平量子化学方法（如FCI, CASSCF, DMRG, NEVPT2等）求解$H_{emb}$，获得杂质区域的高精度电子结构信息。

DMET的一个重要特性是，当高水平求解器是变分量子化学方法时（如FCI和CASSCF），通过增大杂质区域的尺寸，DMET的精度可以系统性地提高。这使得DMET在平衡计算成本和精度方面具有灵活性。

1.3 技术难点与传统DMET的局限性

传统DMET的主要技术难点在于浴轨道的构建高度依赖于低级参考波函数。当使用基态Hartree-Fock或ROHF作为参考时，所产生的浴轨道会优先优化基态的描述，从而导致对激发态的描述不足。这在以下方面表现出局限性：

  • 激发能的过高估计：基态偏置使得激发态在嵌入空间中的能量被高估。
  • 自旋态多重度预测不准确：对于涉及多个自旋态的体系，传统DMET可能无法正确预测基态的自旋多重度，也无法准确描述不同自旋态之间的能量差。
  • 起始点依赖性强：ROHF参考波函数的自旋多重度选择对计算结果影响巨大，缺乏普适性。
  • 对镧系元素体系的收敛性问题：ROHF在镧系配合物中进行自洽场 (SCF) 收敛可能非常困难，需要复杂的策略。

1.4 方法细节：SA-DMET, CAHF和AO-DMET

为解决上述局限性，本文提出了多项创新：

1.4.1 态平均密度矩阵嵌入理论 (SA-DMET)

SA-DMET是本文的核心创新点。它通过将DMET的低级参考点从单态Hartree-Fock扩展到态平均完全活性空间自洽场 (SA-CASSCF) 波函数。SA-CASSCF是一种多参考方法，它同时优化多个目标态的能量加权平均值，从而在波函数中自然地包含了激发态的信息。

在SA-CASSCF中，要最小化的能量泛函是目标态能量的加权平均值 (Eq. 3)： $E^{SA} = \sum_I w_I \langle \Psi_I | \hat{H} | \Psi_I \rangle$ 其中，$w_I$ 是目标态的权重，$\sum w_I = 1$，而 $|\Psi_I \rangle$ 是活性空间中由所有可能组态 ($\Phi_\Gamma$) 展开的正交多组态波函数 (Eq. 4)： $|\Psi_I \rangle = \sum_\Gamma C_{\Gamma I} |\Phi_\Gamma \rangle$

在SA-DMET中，浴轨道构建方式与HF-DMET相同，但其环境区域的单电子约化密度矩阵 (1-RDM) 是基于态平均的。通过对这个态平均1-RDM的环境块进行对角化 (Eq. 5, 6)，并利用其本征解来形成环境中局域正交轨道的线性组合 (Eq. 7)，从而得到三类环境轨道（浴轨道、核心轨道、虚拟环境轨道）。

SA-CASSCF的引入使得浴轨道的构建包含了激发态的信息，从而实现了对基态和激发态的更平衡描述。由于典型的SA-CASSCF计算只涉及有限数量的活性轨道（中心金属离子和少数配体轨道的d/f轨道），嵌入子空间的大小仍然远小于整个系统，因此SA-DMET的计算效率与原始单态DMET相当。

值得注意的是，CASSCF-based DMET保留了HF-based DMET的一个重要特性：当使用与全电子CASSCF相同活性空间选择的CASSCF来求解嵌入杂质哈密顿量时，它可以精确恢复全电子CASSCF解。这在图1中有示意性证明，说明了SA-DMET在包含激发态信息的同时，仍能保持理论上的严格性。

1.4.2 组态平均Hartree-Fock (CAHF) 理论

SA-DMET的起始点是SA-CASSCF，这对于大型体系而言计算量仍然较大。本文进一步提出了使用组态平均Hartree-Fock (CAHF) 作为DMET的起始点。CAHF在特定情况下是SA-CASSCF的极限情况，即当活性空间中的所有自旋态都以相等权重平均时。在多参考计算中，特别是对于过渡金属和镧系配合物，一种常用的态平均选择是平均所有可能的dn或fn自旋态，这与从头算配体场理论 (AILFT) 中的做法一致。

CAHF的能量泛函定义如Eq. 8所示，它通过为活性轨道分配分数占据数来实现对所有相关自旋态的平均。Roothaan在早期工作中已经提出了这种为保持原子和某些分子中电子密度旋转对称性的分式占据方案。CAHF的优点在于其能量泛函相对于轨道旋转更平滑，因此更容易收敛，尤其是在处理镧系配合物时，其SCF收敛性远优于ROHF。这使得CAHF成为SA-DMET构建浴轨道的一个更高效、更鲁棒的选择。

1.4.3 非正交原子轨道密度矩阵嵌入理论 (AO-DMET)

作为传统LO-DMET的推广，Ai等人最近开发了一种非正交原子轨道DMET (AO-DMET)。与LO-DMET使用局域正交轨道定义杂质不同，AO-DMET直接使用非正交原子轨道来划分杂质和环境。这种方法的优势在于：

  • 避免配体污染：LO-DMET中使用的Löwdin正交化过程可能会将配体字符引入到被认为是金属中心的杂质轨道中，从而影响计算精度。AO-DMET通过直接使用原子轨道，更真实地保留了活性轨道的金属中心特征。
  • 更好的杂质空间描述：通过奇异值分解 (SVD) 分析表明，AO-DMET构建的嵌入子空间与全电子活性空间更加接近，尤其是在描述那些通过LO-DMET被低估的轨道时 (图3)。

AO-DMET的引入为DMET的精度提升提供了另一条途径。它与SA-DMET分别从不同的角度（杂质定义与激发态信息纳入）优化嵌入理论，但都服务于同一个物理要求：即嵌入子空间应尽可能完整地涵盖全电子活性空间。将这两种策略结合起来，有望在适度增加浴轨道数量的情况下，获得接近全电子计算的精度。

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

为了全面评估SA-DMET、CAHF以及AO-DMET的性能，本文选择了三类典型的单中心强关联体系进行基准测试，涵盖了磁各向异性、光致发光和自旋态能量学等关键性质。所有计算均在PySCF软件包的本地扩展版本中进行，并考虑了自旋自由的eXact-2-component (SFX2C) 哈密顿量来处理标量相对论效应，以及通过自旋轨道平均场 (SOMF) 近似处理Breit-Pauli哈密顿量中的自旋轨道耦合。

2.1 3d和4f单离子磁体 (SIMs) 的磁各向异性

SIMs是近年来在信息存储、量子传感和量子计算领域备受关注的材料。其关键特征是磁各向异性，源于自旋轨道耦合 (SOC) 和晶体场分裂的相互作用。本研究主要关注零场分裂 (ZFS) 参数 D 的预测。

  • 3d-SIMs：选择了13个过渡金属配合物（8个Co(II)配合物1-8，5个Fe(II), Ni(II), Mn(III), V(III), Mn(II)配合物9-13）。活性空间主要由中心金属离子的3d轨道构成。杂质区域选择为过渡金属原子。

    • LO-DMET性能 (表1)：
      • LO-DMET@ROHF：相对于全电子结果，平均绝对误差 (MAE) 为4.13 cm⁻¹，平均绝对相对误差 (MARE) 为9.1%。特别是在配合物4, 5, 7中存在较大误差，表明基态ROHF参考对激发态的描述不足。
      • LO-DMET@CAHF：MAE显著降低至1.48 cm⁻¹，MARE降至3.3%。结果明显优于LO-DMET@ROHF，尤其是对4, 5等误差较大的体系改进显著。这归因于CAHF的态平均策略能够更平衡地描述所有自旋自由激发态。
    • AO-DMET性能 (表1)：
      • AO-DMET@ROHF：MAE为1.95 cm⁻¹，MARE为3.6%。相较于LO-DMET@ROHF有明显提升，说明原子轨道定义杂质的优势。
      • AO-DMET@CAHF：MAE进一步降至0.93 cm⁻¹，MARE为2.1%。这是所有DMET变体中表现最好的，将DMET@CAHF与AO-DMET结合，实现了最佳精度。

  • 4f-SIMs (表2)：以镝茂阳离子 $[(Cp^{iPr5})Dy(Cp^*)]^+$ (DyCp) 为例，其展现出高的磁化反转有效势垒和磁阻塞温度。活性空间为七个Dy-4f轨道 (9e, 7o)。

    • LO-DMET性能 (表2)：
      • LO-DMET@ROHF：MAE为12.8 cm⁻¹，与全电子结果基本一致，表明在某些情况下ROHF参考也能提供合理结果。
      • LO-DMET@CAHF：MAE为11.6 cm⁻¹，精度与LO-DMET@ROHF相当。
    • AO-DMET性能 (表2)：
      • AO-DMET@ROHF：MAE显著降低至0.9 cm⁻¹。
      • AO-DMET@CAHF：MAE为5.0 cm⁻¹。
    • 结论：对于镧系SIMs，CAHF作为起始点在精度上与ROHF相比没有显著优势，但其SCF收敛性更鲁棒，具有实际操作优势。

2.2 Ce(III) 配合物的4f-5d跃迁

以Ce-BpMe配合物为例，研究了其4f-5d激发态能量，这类跃迁与光致发光性能密切相关。活性空间包括Ce的4f和5d轨道 (1e, 12o)。杂质区域选择为中心Ce原子和6个配位N原子。

  • DMET性能 (表3)：
    • LO-DMET@ROHF：MAE为0.56 eV。
    • LO-DMET@CAHF(f)：仅对4f轨道进行组态平均，MAE为0.46 eV。
    • LO-DMET@CAHF(fd)：对4f和5d轨道都进行组态平均，等效于态平均CASSCF，MAE显著降低至0.22 eV。
    • AO-DMET@ROHF：MAE为0.19 eV，优于所有LO-DMET@ROHF和LO-DMET@CAHF(f)结果。
    • AO-DMET@CAHF(f)：MAE为0.07 eV。
    • AO-DMET@CAHF(fd)：MAE仅为0.03 eV，实现了极高的精度。这突出表明了结合态平均嵌入和AO-DMET在描述局部激发态方面的巨大优势。

2.3 过渡金属配合物的自旋态能量学

自旋态能量学对于许多生物无机化学体系至关重要，也是理论计算的挑战。研究选择了三例过渡金属配合物：Fe(III)L2(OH)、[V(III)(ddpd)2]³⁺和[Co(III)(N3S3)]³⁺。

  • Fe(III)L2(OH) (表4)：计算了低自旋 (LS) S=1/2 和中间自旋 (IS) S=3/2 态相对于高自旋 (HS) S=5/2 态的相对能量。

    • LO-DMET@ROHF：结果对ROHF参考的自旋态选择高度依赖。使用ROHF(1/2)或ROHF(3/2)作为起始点时，LS和IS态的相对能量被显著低估，甚至错误地预测IS态为基态。ROHF(5/2)的误差约为0.1 eV。
    • LO-DMET@SA-CAS/CAHF：与全电子结果的误差在0.05 eV以内，展现出显著改进，提供了准确的自旋态能量。
    • AO-DMET：起始点依赖性显著降低。AO-DMET@ROHF(S=1/2)和ROHF(S=3/2)也能给出定性正确结果，所有起始点的误差均在0.1 eV以内。
    • 轨道重叠分析 (表5)：通过全电子活性空间与DMET嵌入子空间之间的RMS奇异值重叠偏差1-σF，进一步证实AO-DMET能够构建更接近全电子活性空间的嵌入子空间，从而降低了参考依赖性。

  • [V(III)(ddpd)2]³⁺ (表6)：计算了最低三个单重态和三重态相对于三重态基态T0的相对能量。

    • LO-DMET@RHF/ROHF(1)：未能给出令人满意的结果。DMET@RHF显著低估S1能量，高估T2。DMET@ROHF(1)高估了T1, T2, S3能量。
    • LO-DMET@SA-CAS/CAHF：能够准确估算所有低能自旋态的相对能量。
    • AO-DMET：所有不同起始点的结果都与全电子结果吻合良好，因为CASSCF活性空间（V(III)的3d和2p轨道）已通过原子轨道定义包含在杂质中，显著降低了起始点的影响。

  • [Co(III)(N3S3)]³⁺ (表7)：该配合物以单重态 (S=0) 为基态，三个三重态 (S=1) 由于空间对称性而接近三重简并。

    • LO-DMET@RHF：单重态-三重态能隙被严重高估了2.7 eV。
    • LO-DMET@ROHF(1)：虽然有所改进，但破坏了三重态的简并性。
    • LO-DMET@SA-CAS/CAHF：结果与全电子计算非常接近，并保留了三重态的简并性，误差在0.15 eV以内。
    • 结论：对于自旋态能量学，DMET@SA-CAS和DMET@CAHF仍是最佳选择，它们避免了选择起始点的模糊性。

综上所述，SA-DMET，特别是与CAHF和AO-DMET结合时，在所有测试体系中都表现出显著的精度提升。SA-DMET通过引入激发态信息解决了传统DMET的基态偏置问题，而AO-DMET则通过更真实的杂质定义避免了正交化引入的误差，两者协同作用，共同提升了DMET的普适性和准确性。

3.1 软件包与基础框架

本研究中所有的量子化学计算都通过一个本地扩展版本的 PySCF 软件包 (参考文献 73) 来执行。PySCF是一个功能强大的开源量子化学程序包，提供了从头算方法（如Hartree-Fock, DFT, CASSCF）到多参考方法（如NEVPT2）的广泛功能。其模块化的设计使得研究人员可以方便地进行二次开发和方法扩展，为SA-DMET的实现提供了坚实的基础。

在处理相对论效应方面，研究采用了自旋自由的eXact-2-component (SFX2C) 哈密顿量 (参考文献 74-76) 来纳入标量相对论效应。对于自旋轨道耦合 (SOC)，则通过自旋轨道平均场 (SOMF) 近似 (参考文献 77, 78) 来处理Breit-Pauli哈密顿量，这是一种在多参考波函数框架下处理SOC的有效方法。

3.2 DMET的实现细节

本工作中的DMET实现采用了其单次单杂质 (one-shot single-impurity) 变体，与作者之前的研究 (参考文献 20-22, 38) 类似。核心实现包括：

  1. 体系轨道局域化：将分子轨道（或原子轨道，对于AO-DMET）局域化，以便清晰地定义杂质区域和环境区域。
  2. 低级参考波函数计算：使用PySCF计算全体系的低级波函数。对于SA-DMET，这涉及SA-CASSCF计算；对于CAHF-DMET，则进行CAHF计算；对于单态DMET，则使用RHF/ROHF。
    • SA-CASSCF/CAHF的自洽场 (SCF) 迭代：采用了二阶自洽场 (SOSCF) 算法 (参考文献 70) 来优化轨道，以确保SCF的快速和稳定收敛。为了生成一个良好的初始猜测轨道，研究借鉴了原子价活性空间 (AVAS) 方法 (参考文献 71)，这有助于在复杂体系中识别出合适的活性空间轨道并加速收敛。
  3. 态平均1-RDM的构建：根据低级参考波函数（无论是单态HF/ROHF还是SA-CASSCF/CAHF）计算全体系的1-RDM。对于SA-DMET和CAHF-DMET，这个1-RDM是态平均的。
  4. 浴轨道构建：根据杂质区域的定义和全体系1-RDM的环境块，通过对环境块密度矩阵的对角化（如Eq. 5, 6, 7所示），提取出浴轨道。浴轨道选择的标准参数 $\varepsilon_{bath}$ 通常取一个小的正值，本研究中默认为$10^{-13}$。
  5. 嵌入杂质哈密顿量的构建：将全体系哈密顿量投影到由杂质轨道和浴轨道构成的嵌入杂质空间，形成 $H_{emb}$ (Eq. 1)。
  6. 高精度求解器：利用PySCF中的CASSCF和NEVPT2模块，对嵌入杂质哈密顿量进行高精度求解。CASSCF用于处理静态关联，而NEVPT2则进一步纳入了动态关联效应。

3.3 AO-DMET的特定实现

AO-DMET的实现与LO-DMET的主要区别在于杂质的定义方式。在AO-DMET中，杂质直接通过非正交原子轨道来定义，而非局域正交轨道。这意味着在构建杂质区域时，避免了Löwdin正交化等可能引入配体字符的步骤，从而更直接地保留了金属中心的物理性质。其余的浴轨道构建和高精度求解步骤与LO-DMET类似。

3.4 复现指南与开源链接

为了确保研究的透明度和可复现性，作者明确指出：

用于生成和分析结果的脚本和代码，以及本研究中使用的基础数据，均可在以下GitHub仓库中公开获取：

https://github.com/ccme-tmc/SA-DMET-data

对于希望复现或进一步开发相关方法的科研人员，可以按照以下步骤进行：

  1. 下载PySCF：首先，需要安装PySCF量子化学软件包及其依赖项。PySCF的官方文档提供了详细的安装指南。
  2. 获取本地扩展代码：从上述GitHub仓库下载本研究的PySCF本地扩展代码。这些扩展可能包括对DMET、SA-CASSCF和CAHF的特定实现，以及对SFX2C和SOMF的接口。
  3. 数据准备：仓库中应包含所有基准测试体系的分子结构（例如，XYZ格式）和基组信息。确保基组文件在PySCF可访问的路径中。
  4. 配置计算：根据论文中2.4节“Computation details”和Supporting Information中S4节“Additional computational details”描述的计算参数（如活性空间选择、CAHF参数、DMET浴轨道选择标准等），配置各个体系的输入文件。
  5. 运行计算：使用提供的脚本运行DMET、SA-CASSCF、CAHF和NEVPT2计算。脚本通常会管理PySCF模块的调用、结果的解析和后处理。
  6. 结果分析：将计算结果与论文中提供的全电子基准数据进行比较，验证方法的精度和性能。GitHub仓库中可能也包含用于数据分析和绘图的脚本。

通过遵循这些指南和利用开源仓库，其他研究人员可以独立验证本文的发现，并在此基础上开展进一步的计算研究和方法开发。

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

4.1 关键引用文献

本研究建立在量子化学和量子嵌入理论的深厚基础上，并引用了大量经典和前沿文献。以下是其中一些对于理解本工作至关重要的引用：

  1. DMET的提出与发展 (参考文献 9, 17, 18, 50)：
    • Knizia, G.; Chan, G. K.-L. (2012, 2013): 首次提出密度矩阵嵌入理论，奠定了DMET的理论框架和基本思想。这是DMET领域的开山之作。
    • Wouters, S.; Jiménez-Hoyos, C. A.; Sun, Q.; Chan, G. K.-L. (2016): 提供了DMET在量子化学中实际应用的详细指南，涵盖了方法学细节和实现考量。
  2. 量子嵌入理论综述 (参考文献 1-4)：
    • Sun, Q.; Chan, G. K.-L. (2016) 和 Jones, L. O. et al. (2020)：对量子嵌入理论进行了全面的综述，帮助读者理解DMET在更广泛的嵌入方法家族中的位置和重要性。
  3. 多参考波函数理论与CASSCF (参考文献 55, 56, 57)：
    • Andersson, K. et al. (1992) 和 Angeli, C. et al. (2001)：NEVPT2（二阶微扰理论）和 n-electron valence states perturbation theory (NEVPT) 是处理动态关联的关键多参考方法，本文的高级求解器基于此。
    • Werner, H.-J.; Meyer, W. (1981): 提出了二次收敛的MCSCF方法，是SA-CASSCF的基础，对于理解其能量泛函的最小化至关重要。
  4. CAHF理论 (参考文献 66, 67, 69)：
    • Roothaan, C. C. J. (1960): 早期关于开壳层电子系统自洽场理论的奠基性工作，引入了分数占据的概念，是CAHF理论的理论源头。
    • Van Den Heuvel, W. et al. (2016): 描述了镧系(III)配合物中低能晶体场水平的组态平均4f轨道从头算计算，与本工作中使用CAHF的想法密切相关。
    • Zerner, M. C. (1989): 讨论了组态平均Hartree-Fock程序。
  5. AO-DMET的提出 (参考文献 38)：
    • Ai, Y.; Li, Z.-W.; Jiang, H. (2025): 提出了非正交Slater行列式分解的量子嵌入方法，即AO-DMET，是本工作采用AO-DMET方法的基础。这项工作展示了AO-DMET在镧系发光配合物激发能计算中的潜力。
  6. 过渡金属/镧系配合物的应用背景 (参考文献 41-45, 60, 64)：
    • Khedkar, A.; Roemelt, M. (2021) 和 Feldt, M.; Phung, Q. M. (2022): 综述了现代多参考方法在过渡金属化学中的应用。
    • Chibotaru, L. F. (2023): 关于分子纳米磁体计算建模的综述，提供了SIMs磁各向异性研究的背景。
    • Rado’n, M. et al. (2024): 介绍了用于基准测试自旋态能量学的量子化学方法。

4.2 本项工作的局限性评论

尽管SA-DMET结合AO-DMET和CAHF在描述强关联体系的局部激发态方面取得了显著进展，但仍存在一些局限性，为未来的研究指明了方向：

  1. 几何结构弛豫的缺失：本研究主要关注固定分子几何结构下的垂直激发能。然而，许多自旋态能量学和激发态过程会伴随着显著的几何结构变化（结构弛豫）。忽略结构弛豫可能导致对实际激发能和热力学性质的不准确预测。未来的工作需要将结构优化和弛豫状态的激发能计算纳入DMET框架。
  2. 活性空间选择的挑战：虽然AVAS方法有助于活性空间的选择，但在更复杂的镧系配合物中，可能涉及多壳层或多中心活性空间。例如，对于多个f电子或d电子相互作用的体系，准确选择一个既能捕捉所有关键关联，又不会导致CASSCF计算复杂度过高的活性空间，仍然是一个技术挑战。SA-DMET的性能在很大程度上依赖于活性空间选择的合理性。
  3. CAHF的普适性限制：CAHF作为SA-CASSCF的等效情况，仅适用于活性空间中所有自旋态都以相等权重平均的情况。对于需要非等权重或特定自旋态平均的场景，SA-CASSCF仍然是更通用的选择。虽然CAHF计算效率更高且收敛性更好，但其适用范围不如SA-CASSCF广泛。未来的工作可以探索如何扩展CAHF以处理更灵活的态平均方案。
  4. 浴轨道数量的适度增加：论文指出，AO-DMET和SA-DMET的结合会“适度增加浴轨道的数量”。虽然这种增加在计算成本可控范围内，但对于非常庞大的体系或需要极高精度的计算，浴轨道数量的增加仍可能对高水平求解器的计算效率产生影响。优化浴轨道选择策略以在保证精度的同时最小化嵌入空间大小是持续的研究方向。
  5. 强关联电子体系的维度挑战：即使是DMET，当嵌入杂质空间（杂质轨道+浴轨道）变得非常大时，高水平求解器（如CASSCF或NEVPT2）的计算成本仍然会变得很高。对于一些极端强关联的体系，这可能仍然构成挑战。
  6. 理论普适性与实际应用的平衡：虽然SA-DMET提供了更平衡的描述，但其基石SA-CASSCF本身就是一项计算成本较高的任务，尤其是对于拥有大型活性空间或众多态需要平均的体系。因此，在理论的普适性和实际的计算可承受性之间寻找最佳平衡点，是所有高性能量子化学方法的共同挑战。

5. 其他你认为必要的补充

5.1 DMET在强关联体系研究中的重要性

密度矩阵嵌入理论 (DMET) 作为一种前沿的量子嵌入方法，在处理强关联体系，特别是过渡金属和镧系配合物的电子结构方面，展现出巨大的潜力。这些体系因其独特的磁性、催化和光物理性质，在材料科学、生物无机化学和量子信息等领域具有广泛应用。然而，其复杂的电子结构（如强烈的静态和动态关联效应）使得传统的单参考方法（如Hartree-Fock或DFT）难以准确描述，而高精度的多参考方法计算成本又过于高昂。DMET通过“分而治之”的策略，使得研究人员能够在保持计算可承受性的前提下，对体系的关键局部区域进行高精度处理，极大地拓宽了强关联体系研究的范围。

5.2 SA-DMET和AO-DMET的协同效应

本研究最显著的贡献之一是强调了SA-DMET和AO-DMET的协同作用。两者虽然在嵌入理论中解决了不同的问题，但共同指向了同一个目标：构建一个能够最完整、最真实地代表全电子体系活性空间的嵌入子空间。

  • SA-DMET：解决了传统DMET在处理激发态时基态偏置的问题。通过引入态平均的概念，它确保了在构建浴轨道时，激发态的信息能够被公平地纳入考量，从而实现对基态和激发态的平衡描述。这对于准确计算激发能、自旋态能量差和磁各向异性等性质至关重要，因为这些性质本身就涉及多个电子态的相互作用。
  • AO-DMET：解决了传统LO-DMET在定义杂质时可能引入的轨道污染问题。通过直接使用非正交原子轨道作为杂质，它更忠实地保留了金属中心活性轨道的固有特征，避免了Löwdin正交化等操作可能导致的配体字符混入。这种更“物理”的杂质定义方式，使得嵌入子空间能够更准确地反映局部强关联区域的真实电子结构。

当SA-DMET与AO-DMET结合时，它们的效果是互补而非简单的叠加。SA-DMET确保了嵌入子空间在能量维度上能够同时捕捉基态和激发态的特征，而AO-DMET则确保了嵌入子空间在空间维度上能够更准确地表示中心金属离子的活性轨道。这种双重优化策略，使得最终的计算结果在精度上能够大幅逼近全电子计算，且计算成本仅适度增加。例如，在Ce(III)配合物的4f-5d跃迁计算中，AO-DMET@CAHF(fd)实现了惊人的0.03 eV平均绝对误差，充分展示了这种协同作用的强大效力。

5.3 CAHF作为DMET起始点的实用优势

CAHF在此项工作中被证明是SA-DMET的一个高效且鲁棒的替代方案，尤其是在处理特定类型的态平均问题时。其主要优势体现在：

  • 计算效率高：CAHF相对于完全的SA-CASSCF计算，通常具有更低的计算复杂度，因为它将多组态的关联以平均占据数的Hartree-Fock形式来处理。
  • SCF收敛性好：CAHF的能量泛函相对于轨道旋转更平滑，使得其自洽场 (SCF) 迭代更容易收敛。这对于处理镧系元素等收敛性差的体系尤为重要，甚至比ROHF具有显著优势。ROHF在镧系体系中经常面临收敛困难，需要复杂的初始化和步进策略。
  • 精度可接受：在许多情况下，特别是当活性空间的所有自旋态都被等权重平均时，CAHF能够提供与SA-CASSCF相当的精度。这使得它成为在追求高效率和稳定性时的优选起始点。

5.4 未来发展方向

本研究为DMET在激发态和强关联体系中的应用开辟了新途径，同时也指出了几个重要的未来研究方向：

  1. 结构弛豫的纳入：将结构弛豫（即考虑分子几何结构在激发态下的变化）纳入SA-DMET框架，是提升其在实际应用中预测能力的下一步。这将允许研究人员计算绝热激发能和研究光物理过程中的势能面。
  2. 高通量筛选与自旋-声子耦合：对于单离子磁体，SA-DMET@CAHF+NEVPT2的低计算成本优势，使其有潜力应用于高通量筛选具有特定磁性性质的分子，并研究自旋-声子耦合等重要动力学过程，这对于设计新型分子磁体至关重要。
  3. 复杂镧系体系的扩展：本研究中的Ce(III)配合物是f¹体系的简单案例。将该方法扩展到更复杂的镧系配合物，涉及多f电子（如Dy(III)的f⁹）以及可能存在的多中心活性空间，将是更具挑战性的任务，需要进一步的方法学发展。
  4. 与其他高水平方法的结合：DMET框架的灵活性允许其与各种高水平求解器结合。未来可以探索SA-DMET与更先进的多参考方法（如全CI、张量网络态方法）或基于机器学习的量子态表示方法 (如NNQS) 的结合，以进一步提升精度和处理能力。

5.5 结论与展望

本工作通过引入SA-DMET、AO-DMET和CAHF，成功克服了传统DMET在描述强关联体系局部激发态时的主要局限性。这一新范式不仅显著提高了计算精度，还为处理复杂体系提供了更高的效率和鲁棒性。SA-DMET的提出，标志着量子嵌入方法在平衡基态和激发态描述方面迈出了关键一步，为分子磁体、发光材料和催化剂等领域的研究提供了强大的理论工具。随着未来在结构弛豫、多中心体系和新求解器结合等方面的进一步发展，DMET有望在解决量子化学中最具挑战性的问题方面发挥越来越重要的作用。