有限基组中密度最优与能量最优单电子约化密度矩阵（1-RDM）的失配：氢原子中弱占据自然轨道振荡导致的能量代价

《Physical Chemistry Chemical Physics》：Mismatch between density- and energy-optimal one-electron reduced density matrices in finite basis setsOpen Access

　　

原则上，精确基态单电子约化密度矩阵（one-electron reduced density matrix, 1-RDM）同时给出精确电子密度和基态能量。然而，在有限基组中，一个能够紧密复现精确密度的近似1-RDM可能给出出奇差的能量。研究人员利用氢原子研究该效应的一种成因，其中不存在泛函误差和参考密度误差。研究发现，在系综N-可表示（ensemble N-representable）1-RDM上进行密度拟合会产生具有虚假振荡的弱占据自然轨道（natural orbitals），这些轨道对总能量作出异常大的贡献。将密度优化限制在纯态（pure states）可去除这些高能贡献，并使总能量更接近能量最优值。研究人员还识别出氢原子的“密度–能量一致（density–energy consistent）”基组，其中系综密度最优解本质上是纯态，但这些基组既不可转移，也在基组增广（augmentation）下不稳定。尽管纯态限制可能无法为密度–能量失配提供通用补救，但弱占据、振荡自然轨道的能量代价可能超出单电子体系。

研究背景方面，精确基态单电子约化密度矩阵（1-RDM）在完备基组极限下同时对应精确电子密度和基态能量，因此从电子密度重建有限基1-RDM在理论化学中具有重要意义，并可用于从X射线衍射数据提取量子力学信息；密度泛函理论中的Kohn–Sham反演、有限基Kohn–Sham矩阵重建和密度v-可表示性等问题也涉及相关反问题。问题在于，数值证据表明，即便对原子和简单轻元素体系，仅以密度为目标拟合1-RDM往往不能得到能量上合理的有限基1-RDM，更令人满意的1-RDM常需偏离纯密度最优；以往重建还可能存在N-可表示性越界、强制幂等性或未确保全局密度拟合极小等混杂因素，因而密度拟合导致总能量变差的明确机制尚未被识别。研究人员以氢原子为干净体系分离出一种机制：其精确参考密度和1-RDM能量泛函均已知，且系综域上的密度与能量优化均为凸问题；研究比较系综密度最优、纯态密度最优与能量最优1-RDM，并分析基组选择的影响。结论是有限基组中密度最优与能量最优1-RDM可严重失配，密度精度本身不足以衡量重建1-RDM质量；弱占据自然轨道的虚假振荡可带来不成比例的能量代价，纯态限制在单电子H原子上可显著缓解该失配，但特殊“密度–能量一致”基组脆弱且不可转移。其意义在于揭示有限基密度拟合可利用改善密度但高能量代价的自由度，并为在密度基1-RDM重建中引入能量或动能约束提供依据。

关键技术方法（不超过250字）：研究人员以氢原子精确基态密度和精确1-RDM能量泛函为参照，在未收缩STO-mG得到的原始Gaussian型轨道（Gaussian-type orbitals, GTO）(ms)基组（m=2–6）中比较系综N-可表示域内密度最优与能量最优1-RDM；密度优化为凸二次规划，重列为半正定规划（semidefinite programming, SDP）并用CVXPY与Clarabel求解，能量优化化为单电子Hartree–Fock本征值问题；纯态密度优化为非凸问题，采用多起点序列二次规划（sequential quadratic programming, SQP）并以能量最优轨道、Γ^ens_ρ的首个自然轨道和随机向量为起点；用自然轨道分解分析能量贡献，并用线性相关乘积基组考察非唯一性，用Nelder–Mead优化Gaussian指数寻找密度–能量一致基组，并转移至H₂⁺测试。无样本队列来源。

研究结果如下。3.1 密度优化与能量优化（Density optimization versus energy optimization）：通过比较(ms)基组中拟合精确H原子密度所得的系综密度最优Γ^ens_ρ、纯态密度最优Γ^pure_ρ与能量最优Γ_E的总能量和密度误差，发现Γ^ens_ρ在各基组中复现密度更准确但能量显著更高；(2s)中相对Γ_E的能量惩罚达147 mE_h，约为H原子总能量的30%，(6s)中仍比精确值高36 mE_h；用HGBS-9大基组计算的Gaussian密度作参考密度重复计算，结果基本一致，表明核尖并非大能量误差来源；将密度拟合限制为纯态1-RDM后，Γ^pure_ρ仍略具密度拟合优势而能量与Γ_E近乎相同，说明Γ^ens_ρ的大能量惩罚几乎全部与单电子纯态N-可表示性所排除的分数占据有关。3.2 H原子密度–能量失配的起源（Origin of the density–energy mismatch for the H atom）：通过自然轨道分解，将密度写为各自然轨道密度按占据数n_i加权之和、总能量写为n_i?h?_i之和，发现弱占据自然轨道因?h?_i异常高而对总能量贡献不成比例；(6s)中第三自然轨道占据数仅0.002却单独贡献约0.026 E_h；第一自然轨道接近精确1s轨道并随基组扩大趋近，第二自然轨道相对精确2s轨道出现额外振荡且在核附近保留可观振幅，导致高能量期望值；故系综N-可表示性允许密度拟合利用分数占据，产生弱占据但空间变化快的自然轨道，从而提高密度拟合质量却带来大动能/能量代价；纯态限制因要求电子占据单一轨道而去除该机制。3.3 非唯一性对系综密度最优1-RDM的影响（Effect of nonuniqueness on the ensemble density-optimal 1-RDM）：通过采用具有精确乘积线性相关的三个s型Gaussian基组（指数0.3、1.4、2.5，f₁f₃?f₂²=0），构造复现同一拟合密度的系综N-可表示1-RDM族Γ(λ)=P+λA，?0.2880≤λ≤0.2880，发现整个族中密度最优1-RDM均含分数占据且?h?_i异常高的自然轨道；E[Γ(λ)]在许可区间最小值仅?0.1767 E_h，远高于E[Γ_E]=?0.4765 E_h，说明大能量惩罚不能仅归因于线性独立乘积基组强制密度最优1-RDM唯一。3.4 密度–能量一致基组（Density–energy consistent basis sets）：通过以Hilbert–Schmidt距离D(Γ_E,Γ^ens_ρ)为目标、用Nelder–Mead优化Gaussian指数，发现可获得2s–6s密度–能量一致基组，其中Γ^ens_ρ数值上秩一、n₁>0.9999、占据自然轨道为1s型，密度最优总能量与Γ_E几乎相同；但此类基组脆弱：向优化(3s)基组加入指数1.0的s型Gaussian后，E[Γ^ens_ρ]由?0.490514急剧变为?0.369344 E_h；将H的(3s)一致基组用于r_HH=2.0a₀的H₂⁺并采用准确HF/HGBSP1–5参考密度时，D(Γ_E,Γ^ens_ρ)=0.057，密度最优能量为?0.932697 E_h，明显高于同一基组中的能量最优值?1.067399 E_h（未加核间排斥的精确电子能量为?1.1026 E_h）；因此一致基组对增广敏感、即使转移至另一单电子体系也表现差，且需系统特异的非线性优化，不具通用实用性。

讨论部分总结：H原子结果表明密度拟合不必然产生能量准确的1-RDM，罪魁是具大动能贡献的振荡自然轨道发生分数占据；更广泛原因是电子密度只依赖1-RDM对角，而动能还依赖其非对角行为，拟合密度控制势能贡献但不直接控制动能，这为在密度基1-RDM重建中采用能量或动能约束提供理由，以防止密度拟合引入对密度贡献小但能量代价过大的自然轨道。对H原子，施加纯态N-可表示性通过要求电子占据单一轨道消除虚假自然轨道的高能贡献；但该机制不预期推广到多电子体系，因为分数自然轨道占据可由纯相关基态产生而不一定表示激发态污染；且对偶数N，Coleman的系综N-可表示条件与自然自旋轨道占据数成对简并据报足以保证纯态N-可表示性，这意味着纯态限制对单重态可容许1-RDM不增加超出Coleman条件的约束，故纯态限制不太可能是通用解决方案。

研究结论部分翻译：在完备基组中，精确基态密度和能量对应同一1-RDM。研究结果表明，在有限基组中，即使参考密度和1-RDM能量泛函均精确，这种对应关系也可严重破裂。对H原子，最佳复现密度的系综N-可表示1-RDM尽管密度拟合极准，其能量仍可能显著差于能量最优1-RDM；因此，仅凭密度精度不是重建有限基1-RDM质量的可靠指标。对H原子，密度拟合的大能量惩罚来自弱占据自然轨道的虚假振荡；限制为纯态1-RDM可阻止占据这些轨道并去除几乎全部该惩罚，同时保持良好密度拟合。特别优化的基组可使密度与能量最优近乎重合，但这种一致性在基组增广下不稳定，且甚至不能转移至另一单电子体系。纯态限制对H原子有效的原因是一电子体系特有的；该限制本身不应被视为密度–能量失配的通用解决方案。更广泛的教训是，有限基密度优化可利用改善密度但携带不成比例能量代价的自由度；多电子体系是否发生类似效应以及何种机制负责，仍是开放问题。