镭单硫化物 RaX（X=O,S,Se）和 RaO± 离子的电子性质：基于先进量子化学方法的深度解析

来源论文: https://arxiv.org/abs/2603.24590v1 生成时间: Mar 26, 2026 06:17

执行摘要

这项研究是对镭单硫化物（RaX，其中 X = O, S, Se）及其离子 RaO± 的电子性质进行的全面理论调查。作者团队巧妙地结合了全相对论和部分相对论的量子化学方法，旨在揭示这些重原子分子的微观结构和宏观可观测量。研究的核心在于运用精确双分量哈密顿量结合耦合簇（CCSD(T)+X2C）方法计算偶极矩和极化率，并采用内部收缩多参考组态相互作用（MRCI+Q+ECP+SO）方法计算势能曲线。核心发现包括中性 RaX 分子具有极大的永久偶极矩和显著的偶极极化率，但其最低电子态之间的弗兰克-康登因子（FCFs）高度非对角，这直接影响了它们作为激光冷却候选分子的适宜性。这些特性均被归因于中性物种的二价化学键合性质。此外，研究还对 RaO± 离子进行了分析，展示了离子物种极小的激发能。这项工作不仅提供了对这些奇异分子化学键合的新见解，也为未来在精确物理实验中利用这些分子提供了重要的理论基础，尽管它们在激光冷却方面存在挑战。

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

1.1 核心科学问题

镭单硫化物 RaX (X = O, S, Se) 及其离子 RaO± 的研究，站在了当前物理与化学交叉领域的前沿。这项工作聚焦于解决一系列核心科学问题，这些问题不仅对于理解重元素化学键合的本质至关重要，也对探索基础物理学中的新现象具有深远意义。

首先，重原子分子作为探针是当前物理学研究的热点。含有重放射性原子的双原子分子，如 RaX，被认为是探测与宇称（P）和时间反演（T）对称性破缺相关的新物理效应的独特体系 [1-3]。这类分子中，强相对论效应使得内部电场极度增强，从而放大可能由超出标准模型的物理引起的基本对称性破缺信号（例如电子电偶极矩）。因此，对这些分子的电子结构进行高精度理论计算，是指导和解释未来实验的关键。

其次，镭单硫化物相比镭单卤化物研究不足。尽管镭单卤化物（如 RaF, RaCl, RaBr, RaI）已经受到了广泛关注，并被认为是激光冷却的理想候选者 [4, 5, 10-14]，但镭单硫化物（RaO, RaS, RaSe）的研究相对较少。这主要是因为镭作为重元素，其复杂的电子结构和强烈的相对论效应给理论计算带来了巨大挑战。现有的文献仅限于 RaO 基态中 P, T 奇数参数的计算 [22-24]，而对整个系列及其激发态的全面电子结构和性质的深入研究尚属空白。

第三，理解复杂化学键合的本质。镭是碱土金属（AEM）中最重的元素之一，其电正性极高。硫族元素（X=O, S, Se）的电负性也很高。这种组合通常导致高度极性的化学键合，其中 AEM 的两个价电子会转移到硫族元素。这种特殊的键合机制，尤其是在重元素存在下，会受到强烈的相对论效应影响。因此，深入探究 RaX 分子的势能曲线、偶极矩、极化率和分散系数，对于揭示这些奇异分子的化学键合特性，以及它们如何与宏观可观测性质（如激光冷却的适宜性）相关联，具有重要意义。此外，对 RaO± 离子的研究，则进一步扩展了对带电重原子分子体系中离子键合和电子结构的理解。

1.2 理论基础

该研究的理论基础建立在两个核心支柱之上：量子化学和相对论理论。对于含有重原子（如镭）的体系，这两者都不可或缺。

量子化学提供了描述分子电子结构和性质的数学框架。它基于薛定谔方程，通过求解多电子波函数来获得能量和其他可观测性质。然而，精确求解多电子体系的薛定谔方程是不可能的，因此需要采用近似方法。本研究采用的量子化学方法包括：

1. 多参考组态相互作用（MRCI）：这是一种处理强电子关联效应的先进方法，尤其适用于描述基态和激发态具有多参考特征的体系。当一个体系的电子波函数不能被单个哈特里-福克行列式很好地描述时（例如，某些激发态或键断裂过程），MRCI 能够通过包含多个参考组态函数来捕捉这种多参考特性。+Q 修正（Davidson correction）被用于近似地校正 MRCI 计算中由于截断组态空间而引起的完整组态相互作用效应，从而提高能量的准确性。

2. 耦合簇理论（Coupled-Cluster, CC）：特别是 CCSD(T) 方法，被认为是计算单参考体系电子关联效应的“金标准”。它通过指数算符来描述电子波函数，该算符包含单激发（S）、双激发（D）和通过微扰理论处理的近似三重激发（T）。CCSD(T) 在处理动态关联方面表现出色，并能提供高精度的能量和性质。

3. 方程运动耦合簇（Equation-of-Motion Coupled-Cluster, EOM-CC）：EOM-CC 是 CC 理论的扩展，专门用于计算激发态能量和性质。EOM-CCSD 特别适用于描述主要由单激发主导的激发态，能够提供与 CCSD(T) 相当的精度。

4. 时间无关耦合簇极化传播子（TI-CCSD）：该方法用于计算分子的动态极化率，特别是在虚频率下的极化率，这是计算分散系数的关键。

相对论理论对于含有重原子（如 Ra，原子序数 Z=88）的体系至关重要。随着原子序数的增加，内层电子的速度接近光速，导致其质量和能量发生显著变化。这些效应如果不加以考虑，将导致化学键长、解离能和激发能等性质的严重不准确。本研究中应用的相对论方法包括：

1. 精确双分量哈密顿量（Exact Two-Component Hamiltonian, X2C）：X2C 是一种先进的全电子相对论方法，它将四分量狄拉克哈密顿量有效转化为双分量形式，从而在计算效率和精度之间取得平衡。X2C 能够自洽地处理相对论效应，包括质量-速度效应、达尔文项和自旋-轨道耦合，并将其整合到量子化学计算中。

2. 小核赝势（Small-Core Pseudopotentials, ECP）：为了降低计算成本，尤其是在处理极重原子时，赝势方法被广泛采用。ECP 用一个有效势来替代内层电子和原子核的综合作用，从而只显式处理价电子。本研究使用的 ECP78MDF（针对 Ra）和 ECP10MDF（针对 Se）不仅简化了计算，还在其中隐含了部分相对论效应。结合 ECP 后，可以将 Pauli-Breit 算符的对角化引入以处理剩余的自旋-轨道耦合效应。

3. 自旋-轨道耦合（Spin-Orbit Coupling, SO）：SO 耦合是相对论效应的一部分，描述了电子自旋与轨道角动量之间的相互作用。对于重原子，SO 耦合会导致简并态的分裂，并改变势能曲线的形状。本研究通过对 Pauli-Breit 算符的对角化，以微扰方式将 SO 耦合纳入考虑。

1.3 技术难点

这项研究面临多重技术挑战，这些挑战主要源于镭原子本身的性质及其在分子体系中的复杂行为。

1. 重元素的计算成本和精度平衡：镭（Ra）是元素周期表中非常重的元素（Z=88），这意味着其内层电子运动速度极快，必须采用相对论量子化学方法。全电子四分量相对论计算虽然精度最高，但计算成本极高，通常无法应用于中等大小的分子体系。因此，如何在保证合理精度的前提下，使计算在可接受的时间内完成，是最大的挑战之一。本研究通过结合 X2C 全电子方法和 ECP 赝势方法来解决这一问题，但这种混合方法本身也带来了如何准确评估其精度和局限性的问题。

2. 电子关联和相对论效应的协同处理：对于 RaX 这样的极性分子，电子关联（特别是多参考关联）和相对论效应都非常重要且相互影响。如何在一个统一的框架内准确描述这两种效应，是理论方法的关键挑战。MRCI+Q 能够处理多参考关联，但其结合 ECP+SO 是一种部分相对论处理。而 CCSD(T)+X2C 提供了高精度的单参考相对论关联，但可能不适用于强多参考特征的体系。选择合适的组合来最大化精度同时控制成本是关键。

3. 激发态和势能曲线的准确描述：准确计算基态和多个激发态的势能曲线是理解分子光谱学和化学反应性的基础。这对于镭单硫化物尤其困难，因为其电子结构复杂，激发能级可能彼此接近，导致态与态之间存在强烈的混合。此外，势能曲线在长程区域的渐进行为对于确定解离能和分子稳定性很重要，但本研究为了节省计算成本，并未完全解析渐进行为。

4. 弗兰克-康登因子（FCFs）的准确性：FCFs 是评估分子是否适合激光冷却的关键指标。对于镭单硫化物，由于其独特的二价键合特征，电子激发可能导致键长发生显著变化，进而导致 FCFs 高度非对角。准确捕捉这些键长变化，并计算可靠的 FCFs 矩阵，是评估其激光冷却潜力的技术挑战。

5. 高精度偶极矩和极化率的计算：这些分子被预测具有极大的偶极矩，这对于外部电场中的分子操纵和探测对称性破缺效应至关重要。然而，高精度地计算这些性质，特别是动态极化率和分散系数，需要考虑复杂的电子关联和相对论效应，并依赖于精确的波函数。

6. 基组的选择和收敛：对于重原子，选择合适的原子轨道基组是一个细致的工作。基组必须足够大以描述价电子的极化和扩散效应，同时不能过大以致计算无法进行。本研究使用了各种增强型相关一致性基组（aug-cc-pVnZ, aug-cc-pwCVnZ, dyall.avnZ），并通过多套基组和软件进行比较，以评估结果的可靠性。

7. 自旋-轨道耦合的复杂性：重原子中的 SO 耦合效应非常显著，它会引起势能曲线的重排和分裂，从 A-S 态到 Ω 态的转变是复杂的。准确地处理 SO 耦合，尤其是在多参考背景下，需要精细的理论和计算策略。

1.4 方法细节

该研究采用了多种先进的量子化学方法组合，以应对上述技术挑战，旨在在计算可行性与结果精度之间取得最佳平衡。主要方法及其具体应用如下：

1.4.1 势能曲线 (PES) 和光谱常数计算

• 方法名称：内部收缩多参考组态相互作用结合 Davidson 修正 (MRCI+Q)，并使用小核赝势 (ECP) 和自旋-轨道 (SO) 耦合的微扰处理 (MRCI+Q+ECP+SO)。
• 核心思想：
  • MRCI 适用于描述基态和激发态具有多参考特征的体系，它通过包含多个参考组态来捕捉静态电子关联。+Q 修正用于近似校正截断组态空间带来的动态关联误差。
  • ECP78MDF 赝势（针对镭）和 ECP10MDF 赝势（针对硒）用于将重原子内层电子效应纳入考虑，大幅降低计算成本。这些赝势已隐含了大部分标量相对论效应。
  • SO 耦合通过对 Pauli-Breit 算符的对角化进行微扰处理，以描述重原子中自旋和轨道角动量之间的相互作用，从而得到 Ω 态的势能曲线和光谱常数。
• 基组选择：
  • 镭 (Ra)：ECP78MDF 小核赝势，并辅以 aug-cc-pVQZ-PP 基组来描述剩余的 10 个价电子。PP 表示与赝势兼容。
  • 氧 (O)：Dunning 的 aug-cc-pwCV5Z 基组，包含 s, p, d, f, g 函数，用于描述所有 8 个电子。
  • 硫 (S)：aug-cc-pwCV5Z 基组，用于描述所有 16 个电子。
  • 硒 (Se)：ECP10MDF 小核赝势，并辅以 aug-cc-pV5Z-PP 基组来描述剩余的 24 个电子。
• 活性空间定义：
  • RaO：最初活性空间由氧的 2s2p 和镭的 6s6p7s 轨道组成。经过验证，最终决定包含氧的 1s 轨道和镭的 6s6p 轨道作为闭壳轨道，其余作为活性空间 ((6,3,3,1) for neutral RaO)。对于 RaO+，活性空间更紧凑 ((8,3,3,1))。对于 RaO-，活性空间为 (6,2,2,0)。
  • RaS：活性空间包含硫的 3s3p 和镭的 7s 轨道，并考虑了镭的 7p 和 6d 轨道对较低电子态的贡献。最终活性空间为 (6,3,3,1)，闭壳轨道为 (5,2,2,0)。
  • RaSe：活性空间包含硒的 3s3p3d4s4p 轨道和镭的 7s 轨道，并考虑了镭的 7p 轨道。活性空间为 (12,6,6,2)。
• 计算细节：
  • 使用 MOLPRO 软件包进行计算。由于 MOLPRO 不支持无限点群 C∞v，所有计算均在 C2v 点群下进行。分子轨道 (MOs) 和电子态的对称性均在 C2v 下分类。
  • 势能曲线：在感兴趣的核间距范围内密集取点。对基态和激发态的平衡键长 Re、谐振频率 ωe、旋转常数 Be 和电子激发能 Te 进行计算。
  • 振动能量和波函数：使用离散变量表示 (DVR) 方法计算。
  • 弗兰克-康登因子 (FCFs)：基于振动波函数，计算最低五个振动能级间的 FCFs 矩阵。

1.4.2 偶极矩和静态极化率计算

• 方法名称：
  1. 精确双分量哈密顿量结合耦合簇 CCSD(T) 方法 (CCSD(T)+X2C)。
  2. 小核赝势结合耦合簇 CCSD(T) 方法 (CCSD(T)+ECP)。
• 核心思想：
  • CCSD(T) 用于描述高精度电子动态关联。
  • X2C 是一种先进的全电子相对论方法，能自洽地处理相对论效应，并在 DIRAC 软件包中实现。
  • ECP 赝势在 MOLPRO 中用于处理镭的相对论效应。
  • 有限场 (FF) 方法：通过在分子上施加微扰电场 E，然后计算分子能量 U(E) 对电场的导数来获得偶极矩和极化率。一阶导数给出偶极矩，二阶导数给出极化率。
• 基组选择：
  • MOLPRO (CCSD(T)+X2C)：镭使用 aug-cc-pVQZ-X2C 全电子基组；氧和硫使用 aug-cc-pVQZ 基组；硒使用 aug-cc-pV5Z-PP 基组。
  • DIRAC (CCSD(T)+X2C)：所有原子均使用 Dyall 的 relativistic augmented valence n-zeta (dyall.avnZ) 基组 (n=2,3)。明确只包含价-价和亚价-价电子关联。
• 偶极矩和极化率定义：
  • 偶极矩 dz。
  • 平行极化率 α|| = azz。
  • 垂直极化率 α⊥ = 1/2 (axx + ayy)。
  • 平均极化率 ᾱ = (α|| + 2α⊥)/3。

1.4.3 动态极化率和分散系数计算

• 方法名称：时间无关耦合簇极化传播子结合单双激发近似 (TI-CCSD)，并结合 ECP。
• 核心思想：
  • TI-CCSD 用于计算在虚频率 iω 下的动态极化率 α(iω)。
  • 通过 Casimir-Polder 积分 (方程 4) 将两个分子在虚频率下的动态极化率联系起来，计算分子间相互作用的分散系数 Cdisp。
• 基组和赝势：与静态极化率的 CCSD(T)+ECP 方法类似。

1.4.4 激发能基准计算

• 方法名称：方程运动耦合簇结合单双激发近似 (EOM-CCSD)，并结合 ECP。
• 核心思想：用于计算中性 RaX 分子的 X¹Σ+ → a³Π 和 X¹Σ+ → A¹Π 跃迁的激发能。EOM-CCSD 在描述主要由单激发主导的激发态时表现良好。
• 计算目的：评估 MRCI+Q+ECP+SO 方法计算激发能的准确性。

通过对这些方法的综合运用，研究人员能够从不同角度和不同精度级别审视 RaX 体系的复杂性，为揭示其独特的电子性质提供了坚实的理论支撑。

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

这项研究不仅对镭单硫化物进行了详尽的理论计算，还通过与已有的基准体系（如 RaF, SrF, BaF, BaO）进行比较，评估了所用方法的准确性和可靠性。计算所得数据涵盖了势能曲线、光谱常数、偶极矩、极化率和分散系数，这些数据共同描绘了这些奇异分子的全面电子性质。

2.1 中性 RaX 分子 (RaO, RaS, RaSe)

2.1.1 势能曲线和电子态 (无自旋-轨道耦合)

• 整体相似性与原子通道差异 (图 1)：

  • RaO, RaS, RaSe 三种分子的 A-S 势能曲线呈现出显著的相似性，主要区别在于原子通道之间的能量差。这表明其基本的化学键合机制在这三种分子中是相似的。
  • RaO (图 1a)：基态 X¹Σ+ 的电子组态主要为 1σ²2σ²1π⁴。分子轨道（MOs）主要由氧的 2s 和 2p 原子轨道贡献，这表明镭的 7s 电子显著转移到氧的 2p 空穴，形成二价离子键合。最低的激发态 a³Π 和 A¹Π 位于基态上方约 10000 cm⁻¹ 处，共享相似的电子组态，而 B¹Σ+ 态则在 A¹Π 态上方 10234.01 cm⁻¹ 处。镭的 7s 原子轨道主导了 3σ 分子轨道。
  • RaS (图 1b)：基态 X¹Σ+ 电子组态与 RaO 相似 (1σ²2σ²1π⁴)，但硫的电负性较低，导致镭原子吸引电子密度的能力更强。因此，1σ 轨道具有显著的镭 6pz 原子轨道贡献。a³Π 位于基态上方 7846 cm⁻¹，A¹Π 位于约 6000 cm⁻¹，B¹Σ+ 则在 12753.58 cm⁻¹ 处。
  • RaSe (图 1c)：基态 X¹Σ+ 组态与 RaS 相似。a³Π 和 A¹Π 态与基态更接近，分别为 5979.822 cm⁻¹ 和 6764.435 cm⁻¹。B¹Σ+ 态能量更高，位于 24039.00 cm⁻¹ 处，这反映了其跃迁机制（2σ → 2π+）与 RaO/RaS 不同。

• 分子轨道分析 (图 2)：

  • 图 2 直观展示了 RaO (a,b,c)、RaS (d,e,f) 和 RaSe (g,h,i) 的价电子分子轨道。所有情况下，绿色原子为镭，红色原子为硫族元素。这些轨道展示了 X¹Σ+ 基态的 1σ, 2σ, 1π 轨道。RaO 中，轨道主要集中在氧原子上，显示出高度离子化的特征。而在 RaS 和 RaSe 中，镭的 6pz 轨道对 1σ 轨道有显著贡献，表明键合中共享性质的增加。

• 光谱常数 (表 V)：

  • 表格列出了 RaO, RaS, RaSe 最低电子态的平衡键长 Re、谐振频率 ωe、旋转常数 Be 和电子激发能 Te。这些结果由 MRCI+Q+ECP 计算得出，并与 CCSD(T)+ECP 及 EOM-CCSD+ECP 的基准数据进行了比较。
  • 对于 RaO，基态 X¹Σ+ 的 Re 为 2.101 Å，ωe 为 593.304 cm⁻¹。与 CCSD(T)+ECP 相比，键长误差估计约为 1.15%，能量分离偏差约为 8.12% (对于 A¹Π 态)。
  • RaS 和 RaSe 的键长和能量与 RaO 具有相似趋势，但绝对值不同。这些比较对于评估 MRCI+Q+ECP 方法的可靠性至关重要。

2.1.2 跃迁偶极矩和弗兰克-康登因子

• 跃迁偶极矩 (TDM) (图 3)：

  • 图 3 显示了 RaO 和 RaS 在势能最小值附近，X¹Σ+ → B¹Σ+ 和 X¹Σ+ → A¹Π 跃迁的 TDM 绝对值。Σ → Σ 跃迁的 TDM 大于 Σ → Π 跃迁，且在平衡位置附近具有显著数值。

• 弗兰克-康登因子 (FCFs) (图 4)：

  • 图 4 展示了 RaO 和 RaS 在最低五个振动能级间的 FCFs 矩阵。显著的特点是这些矩阵呈现出高度非对角线结构。
  • 解释与激光冷却：这种非对角线结构与 Ra 单卤化物（如 RaF, RaCl）中观察到的高度对角线 FCFs 形成鲜明对比 [4, 5]。Ra 单卤化物中，氟原子的高电子亲和力使得 AEM 的一个 ns² 电子转移到氟的 2p⁵ 轨道中，剩余的价 s 电子仍强烈定域在金属原子周围，因此激发时键长变化不大，导致 FCFs 对角。而 RaX 分子中的二价键合意味着氧的 2s²2p⁴ 组态吸引了镭 7s² 轨道中的两个电子。一个电子强烈定域在氧原子周围，而另一个则在两个原子之间共享，形成了极性共价键。因此，电子激发时键长会显著变化，导致振动波函数重叠差，FCFs 非对角。这意味着 RaX 分子不适合直接激光冷却，因为非对角线 FCFs 会导致激发态衰变到基态的多个振动能级，需要更多的激光频率来循环衰变，效率低下。

2.1.3 自旋-轨道 (SO) 耦合效应

• Ω 态势能曲线 (图 8)：

  • 图 8 展示了 RaO, RaS, RaSe 最低七个 Ω 电子态的势能曲线。这些曲线在不同分子间非常相似，表明 SO 耦合对整体形状的影响是普遍性的。

• 光谱常数 (表 VII)：

  • 表格列出了 Ω 态的光谱常数。与 A-S 态相比，相对论效应并未显著改变这些值，尽管在 Π 态中观察到一些收缩，这可能归因于氧 2p₁/₂ 轨道的收缩。能量移动通常在数百 cm⁻¹ 范围内，且在多数情况下，降低了 Te 值。
  • 这种现象可能是硫族元素的 s 和 p 价轨道相对论稳定化的结果。

2.2 离子 RaO± 分子

2.2.1 RaO+ (阳离子)

• 势能曲线 (图 5)：显示了 RaO+ 的最低六个 A-S 电子态。基态 X²Σ+ 具有离子键合特征，其价电子组态为 1σ²2σ²1π⁴。第一个激发态 A²Π 位于基态上方仅 134.245 cm⁻¹，非常接近。A²Π 态的电子组态为 1σ²1π³2σ²，其激发主要涉及 2πχ → 2σ 跃迁。
• 光谱常数 (表 VI)：列出了 RaO+ 的光谱常数。基态和第一个激发态之间的平衡键长差异约为 0.12 Å，这意味着振动波函数重叠差。
• 单占据分子轨道 (图 7)：展示了 RaF, RaO+, RaO⁻ 的单占据分子轨道。RaO+ 的键合主要为离子键，氧原子占据 Ra⁺ 的一个 7s 和一个 6p 电子。

2.2.2 RaO⁻ (阴离子)

• 势能曲线 (图 6)：显示了 RaO⁻ 的基态 X²Σ+ 势能曲线。原子通道为 Ra(¹S)+O⁻(²P°)。虽然是开壳结构，但平衡时仍有两个电子参与键合，导致其势阱深度不如中性 RaX 和 RaO⁺。
• 光谱常数 (表 VI)：列出了 RaO⁻ 的光谱常数。

2.3 极化率，偶极矩和分散系数

2.3.1 静态偶极矩 (d₀) 和极化率 (α||, α⊥, ᾱ) (表 VIII)

• 结果比较：表格展示了使用 CCSD(T)+ECP, CCSD(T)+X2C (pVQZ-X2C), CCSD(T)+X2C (dyall.av3z) 方法计算的静态偶极矩和极化率。
• 方法一致性：三种计算方法的结果表现出良好的一致性。X2C 计算之间的差异主要源于基组选择（全核-价关联 vs. 仅价-价和亚价-价关联）以及所包含的关联效应。
• 极大的偶极矩：中性 RaX 分子具有惊人的大偶极矩，特别是 RaS 和 RaSe，其偶极矩超过 10 D。这远远大于异核碱金属二聚体（通常 < 5 D）或镭单卤化物（如 RaF, RaCl）[31, 32]。大的偶极矩表明这些分子在外部电场中具有很强的响应性，使其成为操控分子或探测电子电偶极矩的潜在候选者。RaO 的偶极矩略小，但仍远超大多数双原子分子。
• 偶极束缚态：RaO 的大偶极矩以及 RaS 和 RaSe 类似，预示着可能存在激发态的偶极束缚态 (DBS) [33, 34]，这对于未来研究是一个有趣的发现。
• 显著的极化率各向异性：平行 (α||) 和垂直 (α⊥) 极化率之间存在显著差异，这源于低能量 σ→σ 和 σ→π 跃迁的巨大贡献。这意味着分子对沿不同方向的电场响应不同。
• 平均极化率：RaX 分子的平均极化率 (ᾱ) 小于异核碱金属二聚体。例如，LiNa 的平均极化率为 237.7 a.u.，RbCs 为 638.6 a.u. [35]，而 RaX 分子的 ᾱ 值更小。

2.3.2 分散系数 (Cdisp) (表 IX)

• 计算结果：表格列出了 RaO, RaS, RaSe 之间相互作用的分散系数 Cdisp。RaO-RaO 的 Cdisp 为 1394.853 a.u.，RaS-RaS 为 3378.980 a.u.，RaSe-RaSe 为 3985.470 a.u.，跨物种相互作用也已计算。
• 与碱金属二聚体的比较：与异核碱金属二聚体相比，RaX 分子的分散系数更小，这与它们相对较小的平均极化率一致。分散系数的大小反映了分子之间长程色散力的强度，这对于理解冷分子碰撞和凝聚态物理学很重要。

2.4 性能数据和方法评估

• MRCI+Q+ECP+SO 方法的准确性：
  • 在附录中，研究人员将该方法应用于 RaF 基准分子，并与已发表的全相对论计算和实验数据进行比较 (表 X)。对于 RaF，计算所得的 Re, ωe, Be, Te 与参考值吻合良好。例如，X²Σ+1/2 → ²Π1/2 跃迁的能量偏差为 8.44%，优于此前的一些理论计算。
  • 对于 SrF, BaF, RaF 的偶极矩计算，CCSD(T)+ECP 方法与实验值和全相对论计算结果吻合极好 (表 XI)，相对误差不超过 1%。
• EOM-CCSD+ECP 方法的性能：
  • 该方法应用于 BaO 的激发能计算，并与实验值和准相对论赝势计算进行比较 (表 XII)。结果显示 EOM-CCSD+ECP 方法与实验跃迁能吻合良好，尽管存在轻微的高估趋势。与一些非相对论赝势计算相比，本研究通过 ECP46MDF 处理钡原子的内层电子相对论效应，使能量差异减小了 2000 cm⁻¹，显著提高了精度。
• 方法局限性：
  • MRCI+Q+ECP+SO 方法仍是部分相对论处理，可能存在一定程度的精度损失。例如，对于 RaO，基态键长比 Ref [23] 报告的值大 0.06 Å。单重态-单重态激发能与 EOM-CCSD+ECP 结果相比，偏差达到 17.65%。
  • 势能曲线的渐进行为未完全解析，意味着在解离能方面可能存在不确定性。但整体而言，所采用的方法在精度和计算成本之间取得了很好的平衡，为镭单硫化物体系提供了可靠的理论数据。

总的来说，这项研究通过严谨的方法学和详尽的数据呈现，成功地揭示了镭单硫化物及其离子的电子性质，并为未来在精确物理实验和基础物理研究中的应用奠定了基础。

3. 代码实现细节，复现指南，所用的软件包及开源 repo link

这项研究涉及多项复杂的量子化学计算，需要利用专业的软件包和精心构建的计算流程。虽然论文没有提供具体的输入文件或代码，但我们可以根据其描述重构其实现细节和大致复现指南。理解这些实现细节对于任何希望在该领域进行类似研究的科研人员都至关重要。

3.1 所用软件包及功能

本研究主要依赖于两个在量子化学界广受认可的软件包：MOLPRO 和 DIRAC。

3.1.1 MOLPRO

MOLPRO 是一个功能强大的从头算（ab initio）量子化学程序包，特别擅长于多参考方法和耦合簇方法。本研究中 MOLPRO 的关键应用包括：

• 势能曲线 (PES) 和光谱常数：
  • MRCI+Q：用于计算势能曲线和光谱常数，特别是对基态和激发态具有强电子关联效应的体系。MOLPRO 的 MRCI 模块是其核心优势之一，能够执行内部收缩 MRCI 计算，并支持 Davidson 修正（+Q）。
  • 自旋-轨道耦合 (SO)：通过对 Pauli-Breit 算符的对角化，以微扰方式将 SO 耦合整合到 MRCI 计算中。MOLPRO 允许用户在 MRCI 计算后进行 SO 耦合矩阵的构建和对角化。
  • 赝势 (ECP)：MOLPRO 支持使用各种赝势，如 ECP78MDF 和 ECP10MDF，这对于处理重原子（如 Ra, Se）至关重要，因为它们能有效地模拟内层电子和核的综合效应，并包含标量相对论效应。
• 偶极矩和静态极化率：
  • CCSD(T)：用于高精度计算偶极矩和静态极化率。MOLPRO 的 CCSD(T) 模块是其“金标准”功能之一，能够处理电子动态关联。在 ECP 框架下进行 CCSD(T) 计算。
  • 有限场 (Finite-Field, FF) 方法：MOLPRO 允许通过施加外部电场并计算能量对电场的导数来计算电学性质。通过改变电场强度并拟合能量随电场的变化，可以获得偶极矩和极化率。
• 动态极化率和分散系数：
  • TI-CCSD (Time-Independent Coupled-Cluster Polarization Propagator in Singles and Doubles approximation)：MOLPRO 提供了计算极化传播子的功能，这对于获得动态极化率（在虚频率下）至关重要。TI-CCSD 模块能够有效地处理动态关联对极化率的贡献。
• 激发能：
  • EOM-CCSD：用于计算激发态能量。MOLPRO 的 EOM-CCSD 模块能够提供高精度的激发能，特别适用于单激发主导的激发态。

3.1.2 DIRAC 23.0

DIRAC 是一个专门用于全电子（all-electron）相对论量子化学计算的程序包，以其严谨的相对论处理而闻名。本研究中 DIRAC 的应用主要集中在：

• 偶极矩和静态极化率：
  • CCSD(T)+X2C：DIRAC 能够实现精确双分量哈密顿量（X2C）与 CCSD(T) 理论的结合。X2C 方法是 DIRAC 的核心优势之一，它将四分量狄拉克哈密顿量有效地约化为双分量形式，从而在计算效率和精度之间取得平衡，并自洽地包含所有相对论效应，包括自旋-轨道耦合。
  • 全电子相对论处理：与 MOLPRO 中使用的 ECP 不同，DIRAC 在 X2C 框架下通常进行全电子计算，这意味着所有电子（包括内层电子）都以相对论方式显式处理，从而获得更高的精度。
  • 有限场方法：DIRAC 也支持通过有限场方法计算偶极矩和极化率，与 MOLPRO 类似，通过能量对电场的导数来获取。

3.2 基组选择

基组是量子化学计算中描述原子轨道的基础。针对重原子和特定性质，基组的选择至关重要。

• 镭 (Ra)：
  • ECP78MDF 赝势：在 MOLPRO 中，对于镭原子，其 78 个内层电子用 ECP78MDF 赝势代替，并结合 aug-cc-pVQZ-PP 基组描述剩余的 10 个价电子。PP 表示该基组是针对赝势进行优化的。这种选择在保证计算可行的同时，包含了主要的相对论效应。
  • aug-cc-pVQZ-X2C：在 MOLPRO 中进行全电子 X2C 计算时，镭使用专门设计的 aug-cc-pVQZ-X2C 基组。该基组是球谐表示的 (38s,35p,25d,18f,2g)/[11s,10p,8d,5f,2g])，适用于 X2C 哈密顿量。
  • dyall.avnZ (n=2,3)：在 DIRAC 中进行全电子相对论计算时，镭使用 Dyall 的 relativistic augmented valence n-zeta 基组。这些基组是专门为相对论计算设计的，能够更好地描述重原子的电子结构。
• 硫族元素 (O, S, Se)：
  • 氧 (O)：aug-cc-pwCV5Z 基组，用于描述其 s, p, d, f, g 函数。pwCV 表示极化加权核价基组，能够更好地描述核-价关联。
  • 硫 (S)：aug-cc-pwCV5Z 基组，用于描述 spdf 函数。
  • 硒 (Se)：ECP10MDF 赝势，并结合 aug-cc-pV5Z-PP 基组。这与镭的赝势策略类似，用于降低计算成本。
  • dyall.avnZ (n=2,3)：在 DIRAC 中，氧、硫、硒同样使用 Dyall 的相对论基组。
• 氟 (F) (用于基准计算)：aug-cc-pwCV5Z 基组。

基组选择的考量：

• aug 前缀：表示基组中添加了弥散函数，这对于描述负离子、极性分子以及分子间相互作用（如分散力）至关重要。
• pVnZ / pwCVnZ：表示极化价基组（n 表示“zeta”级别，如 QZ 为四zeta，5Z 为五zeta），pwCV 则进一步优化了对核-价关联的描述。
• 赝势 (ECP) vs. 全电子 (X2C)：ECP 是一种折衷方案，它大大降低了计算成本，但可能牺牲部分精度。X2C 是一种更严格的全电子相对论方法，提供了更高的精度，但计算成本也更高。研究中同时使用这两种方法进行交叉验证，以评估 ECP 方法的可靠性。

3.3 活性空间定义

对于多参考方法（如 MRCI），活性空间的定义是决定计算结果准确性和计算成本的关键。

• 基本原则：活性空间通常包含那些对体系的电子结构和化学键合至关重要的价轨道，以及一些可能参与激发或键断裂过程的低能量未占据轨道。闭壳轨道则包含能量较低、对键合贡献较小的内层轨道。
• RaO (中性)：总的基态轨道空间为 (6,2,2,0)。活性空间定义为 (6,3,3,1)，闭壳轨道为 (3,1,1,0)（氧的 1s 和镭的 6s6p）。为了描述镭的激发态，额外增加了 7p 和 6d 轨道。对氧的激发态 O(¹D)，则保持 2s²2p⁴ 组态，无需额外轨道。
• RaO+ (阳离子)：为了描述 Ra+ 的 6s²6p⁶7s¹ 组态和 Ra+ 的 (2D) 激发态，活性空间调整为 (8,3,3,1)。闭壳轨道为镭的 6sp 和氧的 1s ((3,1,1,0))。
• RaO- (阴离子)：只考虑基态通道，镭的 10 个赝核和价电子，以及氧负离子 O⁻(²P°) 的 9 个电子。活性空间为 (6,2,2,0)。闭壳轨道为镭的 6sp 和氧的 1s。
• RaS：活性空间包含硫的 3s3p 和镭的 7s 电子，并加入了镭的 7p 轨道。初步 MCSCF 计算表明镭的 6d 轨道对较低电子态也有贡献。最终活性空间为 (6,3,3,1)，闭壳轨道为 (5,2,2,0)。
• RaSe：活性空间包含了硒的 3s²3p⁶3d¹⁰4s²4p⁴ 组态的 24 个电子以及镭的 10 个电子和激发原子态。总轨道空间为 (12,6,6,2)。

3.4 计算流程与复现指南 (概念性)

虽然无法提供具体的输入文件，但可以勾勒出大致的计算流程和复现步骤：

1. 准备阶段：

   • 选择合适的计算资源：这些计算需要高性能计算集群，特别是对于 MOLPRO 的 MRCI 和 DIRAC 的全电子相对论计算。
   • 安装 MOLPRO 和 DIRAC 软件包，确保所有必要的模块和库均已正确配置。
   • 熟悉 MOLPRO 和 DIRAC 的输入文件格式和关键词。

2. 分子几何优化与势能曲线计算：

   • 定义分子结构：对于双原子分子 RaX，需要定义核间距 R 作为唯一几何参数。
   • 扫描核间距：在感兴趣的核间距范围内（例如 2.0 Å 到 4.0 Å）选择一系列 R 值进行计算，以构建势能曲线。点的密度应在势阱附近更高。
   • 基组和赝势设置：在 MOLPRO 输入文件中指定镭、氧、硫、硒的基组和赝势（如 ECP78MDF, aug-cc-pVQZ-PP 等）。
   • 自洽场 (SCF) 和多组态自洽场 (MCSCF)：首先进行 Hartree-Fock 或 DFT SCF 计算，然后进行 MCSCF 以获得优化的分子轨道和参考波函数。在此阶段定义活性空间。
   • MRCI+Q 计算：基于 MCSCF 轨道，进行 MRCI+Q 计算，以获取基态和激发态的能量。
   • 自旋-轨道耦合：在 MRCI 能量计算后，利用 MOLPRO 的 SO 模块（例如 RS2 或 SO-CI 功能）对 Pauli-Breit 算符进行对角化，得到 Ω 态的能量。选择合适的 A-S 态作为 SO 耦合的参考态。
   • 提取数据：从输出文件中提取每个核间距对应的能量，绘制势能曲线。

3. 光谱常数和振动性质计算：

   • 拟合势能曲线：将势能曲线数据点拟合到合适的解析函数，或直接使用离散点数据。
   • DVR 方法：利用 MOLPRO 的 DVR 功能（或独立的振动程序）计算振动能级和波函数。
   • 光谱常数：从振动能级和拟合曲线中提取平衡键长 Re、谐振频率 ωe、旋转常数 Be 和激发能 Te。
   • 弗兰克-康登因子：使用 DVR 得到的振动波函数，计算不同电子态之间振动能级的重叠积分，得到 FCFs 矩阵。

4. 偶极矩和极化率计算：

   • 选择方法：可选择 MOLPRO 中的 CCSD(T)+ECP 或 DIRAC 中的 CCSD(T)+X2C。
   • 有限场设置：
     • 在输入文件中指定外部电场强度（例如，沿 z 轴施加 +E, -E 和 0 场，E 的值通常较小，如 0.001 a.u.）。
     • 确保在 CCSD(T) 计算中包含 Field 关键词或相应设置。
   • 能量计算：在每个电场强度下执行 CCSD(T) 或 CCSD(T)+X2C 能量计算。
   • 数值微分：根据论文中的公式 (5)，利用不同电场下的能量值进行数值微分，计算偶极矩和极化率分量。

5. 动态极化率和分散系数计算：

   • TI-CCSD 计算：在 MOLPRO 中使用 TI-CCSD 模块计算在虚频率 iω 下的动态极化率 α(iω)。需要指定虚频率网格点。
   • Casimir-Polder 积分：利用论文中的公式 (4)，对 α(iω) 进行数值积分，计算分散系数 Cdisp。这可能需要独立的脚本来实现 Gauss-Legendre 积分。

6. 结果分析和比较：

   • 将所有计算结果汇总，与现有的实验数据和更高级别的理论计算进行比较，评估方法的准确性。
   • 特别注意不同方法（ECP vs. X2C）和基组（pVQZ-X2C vs. dyall.avnZ）之间的差异及其来源。

3.5 开源 Repo 链接

需要明确的是，MOLPRO 是一款商业软件包，通常需要许可才能使用，但其在学术界广泛流行。DIRAC 是一款免费的学术用途软件，通常需要注册才能获取和使用。

• MOLPRO 官方网站：https://www.molpro.net/
• DIRAC 官方网站：https://www.diracprogram.org/

鉴于论文没有提供用于生成结果的具体输入文件或脚本，无法直接指向一个包含复现代码的开源仓库。对于希望复现或扩展此项工作的研究人员，需要自行编写 MOLPRO 和 DIRAC 的输入文件，并根据论文中详细描述的方法学参数（如基组、赝势、活性空间等）进行配置。理解这些软件包的官方文档和教程将是成功复现的关键。

3.6 计算性能考量

重原子相对论量子化学计算对计算资源的需求非常高。主要因素包括：

• 原子序数 (Z)：镭的 Z=88 意味着需要大量的电子和原子轨道来描述，导致计算的规模迅速增大。
• 基组大小：aug-cc-pVnZ 和 dyall.avnZ 基组包含大量的基函数，增加了计算的内存和 CPU 时间需求。
• 电子关联处理：MRCI 和 CCSD(T) 等方法本身就是计算密集型的。MRCI 的计算成本随着活性空间和组态函数数量的增加呈指数增长。CCSD(T) 虽然是“金标准”，但对内存和 CPU 的需求也很大。
• 相对论处理：X2C 等全电子相对论方法虽然比四分量方法高效，但仍然比非相对论或标量相对论 ECP 方法昂贵。
• 自旋-轨道耦合：对 Pauli-Breit 算符的对角化需要额外的计算步骤。
• 几何扫描：为了构建势能曲线，需要对多个核间距点进行独立计算，增加了总计算时间。
• 有限场方法：对偶极矩和极化率的计算需要进行多次不同电场下的能量计算。

因此，这类计算通常在拥有大量 CPU 核、高内存和高速互联网络（如 InfiniBand）的并行计算集群上进行。优化计算参数、选择合适的活性空间和基组截断策略对于控制计算成本至关重要。

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

4.1 关键引用文献解析

该研究建立在深厚的理论和实验基础之上，其引用的文献不仅提供了背景知识，也支撑了方法学的选择和结果的解释。以下是一些关键引用文献及其重要性：

1. [1-3] M. S. Safronova et al. (2017), ACME Collaboration (2018), G. Arrowsmith-Kron et al. (2024): 这些文献确立了重原子和分子在探测 P、T 奇数对称性破缺效应（如电子电偶极矩）方面的关键作用。它们强调了镭等重元素中相对论增强效应的重要性，为 RaX 体系的研究提供了核心动机，即将其作为新物理的潜在探针。

2. [6] R. F. Garcia Ruiz et al. (2020): 这篇 Nature 论文报道了短寿命放射性分子光谱学的实验进展，为 RaX 等分子在实验室中的研究打开了大门。它强调了高精度从头算电子结构计算对于指导实验实现的重要性，是连接理论与实验的关键桥梁。

3. [4, 5, 10-13] T. A. Isaev et al. (2010), T. Isaev et al. (2020), A. D. Kudashov et al. (2014), A. Zaitsevskii et al. (2022), Y. Osika et al. (2021, 2022): 这些研究聚焦于镭单卤化物（如 RaF, RaCl, RaBr, RaI）的电子结构和性质。它们为 RaX 提供了直接的比较对象，特别是关于激光冷却适宜性（FCFs 对角性）。这些文献展示了全相对论或准相对论方法在这些体系中的应用。

4. [17-21] H. Khalil et al. (2013), M. M. F. de Moraes et al. (2018), A. Ivana et al. (1999), M. Khatib et al. (2017, 2018): 这些文献关注更广泛的碱土金属单硫化物（AEM-MX）体系。它们建立了这些分子中离子和二价化学键合的一般特征，为理解 RaX 的键合模式提供了理论背景，并指出了 Ra-MX 体系中强相对论效应带来的独特挑战。

5. [22-24] V. V. Flambaum et al. (2007, 2013, 2020): 这些工作专门讨论了 RaO 基态中的 P、T 奇数参数计算，为本研究的RaO 分析提供了直接的基准和灵感。它们强调了单参考相对论耦合簇方法在该领域的应用。

6. [7-9] K. G. Dyall and K. Faegri (2007), W. Liu (2020), A. V. Oleynichenko et al. (2020): 这些书籍和综述文章提供了相对论量子化学的全面背景，包括理论框架和计算方法（如狄拉克方程、赝势、自旋-轨道耦合、X2C 哈密顿量）。它们是本研究方法学选择的理论支柱。

7. [25, 26] T. Fleig (2012), L. Cheng (2024): 这些文献讨论了全相对论四分量哈密顿量和精确双分量 (X2C) 方法，特别是 Fock-Space Relativistic Coupled-Cluster (FS-RCC) 方法在处理重元素体系时的计算需求。它们解释了为何本研究采用计算成本更低的 ECP 方法作为替代方案。

8. [27] H. Stoll et al. (2002): 这篇论文详细介绍了赝势在降低重原子计算成本方面的应用。它解释了 ECP 的优点和局限性，为本研究选择 ECP 方法提供了理论依据。

9. [37-39] W. Ernst et al. (1985, 1986), S. Haldar et al. (2021): 这些文献提供了 SrF, BaF, RaF 等分子偶极矩的实验和理论基准数据。本研究通过与这些数据进行比较，验证了其 CCSD(T)+ECP 偶极矩计算的准确性。

10. [40-42] R. W. Field et al. (1975, 1980, 1974): 这些论文提供了 BaO 跃迁能的实验和理论数据。本研究通过 EOM-CCSD+ECP 对 BaO 进行基准计算，评估了其激发能计算的准确性。

这些文献共同构成了本研究的知识基础，使其能够利用最先进的量子化学工具，在重元素分子的研究领域取得进展。

4.2 对这项工作局限性的评论

尽管这项研究在探索镭单硫化物及其离子方面取得了重要进展，并采用了先进的量子化学方法，但像所有科学工作一样，它也存在一些局限性。客观地认识这些局限性有助于未来工作的改进和更深入的理解。

1. 部分相对论处理的固有局限性：

   • MRCI+Q+ECP+SO 的精度妥协：论文明确指出，MRCI+Q+ECP+SO 方法是计算可行性与精度之间的权衡结果，虽然 ECP 包含了大部分标量相对论效应，但 SO 耦合仅通过 Pauli-Breit 算符的微扰对角化来处理，而非在一个完全相对论的变分框架中。这意味着在 MRCI 计算中，自旋-轨道效应没有完全自洽地包含在波函数优化过程中，可能导致一些微妙的耦合效应被忽略。论文在附录中也承认，与全相对论基准结果相比，键长、激发能等仍存在一定的偏差（例如，RaO 键长误差约 1.15%，激发能偏差达 8.12% 甚至 17.65%）。结论部分也提到：“基于全四分量相对论计算的进一步改进，仍可能提高势能曲线和跃迁能量的准确性。”
   • X2C 方法的近似性：尽管 X2C 是一个强大的全电子相对论方法，但它仍是对四分量狄拉克方程的近似，尽管是非常好的近似。一些非常细微的相对论效应可能仍未被完全捕捉，尤其是在涉及重元素最内层电子的超高精度计算中。

2. 基组收敛的挑战：

   • 对于镭这样的大原子，即使使用了 aug-cc-pVnZ 和 dyall.avnZ 等高质量的基组，实现完全的基组收敛仍然是一个巨大的计算挑战。基组的截断可能会对极化率、偶极矩和分散系数等对基组大小和弥散函数敏感的性质产生影响。论文在 Table VIII 中展示了 CCSD(T)+X2C 在不同 dyall.avnZ 基组和 aug-cc-pVQZ-X2C 基组下的结果差异，这表明基组的选择确实会影响最终数值，突显了基组收敛的复杂性。

3. 电子关联处理的局限：

   • MRCI+Q 的动态关联近似：MRCI+Q 在处理静态关联方面表现出色，但其动态关联主要通过 Davidson 修正（+Q）进行近似校正。虽然 +Q 是一个有用的经验修正，但它不是变分的，在某些情况下可能不够精确，尤其是在要求极高精度的激发态能量计算中。论文也指出，EOM-CCSD+ECP 方法“倾向于高估最低单重态-单重态激发能”。
   • CCSD(T) 的单参考假设：CCSD(T) 被认为是单参考体系的“金标准”，但如果体系的基态或某些激发态具有显著的多参考特征，CCSD(T) 的准确性会下降。尽管论文使用了 MRCI 来处理多参考态，但 CCSD(T) 在某些特定性质的计算中仍然依赖于单参考框架，可能无法完全捕捉所有多参考效应。

4. 势能曲线渐进行为的未解析：

   • 论文明确指出，势能曲线的计算集中在“势阱周围区域，而未解析渐进行为”（第 2 页）。这意味着对于长程相互作用和准确的解离能，这项研究可能无法提供最高的精度。在某些应用中，例如冷分子碰撞动力学研究，准确的长程势是至关重要的。

5. 离子通道的简化：

   • 在中性 RaX 分子研究中，论文基于镭和硫族元素的高电离势，忽略了离子通道的贡献。这在初步计算中是一个合理的简化，但更全面的分析可能需要考虑这些通道，尤其是在涉及键断裂的物理过程中。

6. 部分激发态的限制性考察：

   • 对于 RaS 和 RaSe，未考虑与 Ra(³D) 相关的激发通道。对于 RaO，也未将 O(¹D) 组态包含在所有渐近线中。五重态被认为“离解性且能量更高”，因此未予分析。这种选择限制了对更高能量激发态或特定解离途径的理解。

7. 缺乏直接的实验数据验证：

   • 对于镭单硫化物本身，目前缺乏实验数据，这意味着本研究的理论结果主要通过与镭单卤化物和钡单氧化物等已充分研究的基准体系进行比较来验证。虽然这些基准计算提供了信心，但直接的实验验证对于这些新颖的分子仍然是最终的“金标准”。

8. C2v 点群近似：

   • 在 MOLPRO 中，为了方便计算，将线性分子（C∞v 点群）在 C2v 点群下进行处理。虽然这是一种常见的实际近似，但可能引入微小的对称性破缺，并对结果产生细微影响。

总而言之，这些局限性反映了在量子化学计算中精度、计算成本和实际约束之间的固有矛盾。然而，这些局限性并未削弱该研究的整体价值，它为镭单硫化物这一复杂而重要的分子体系提供了宝贵的初步深入见解，并为未来的更精确和更全面的研究奠定了基础。

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

5.1 广阔的研究背景与意义

这项对镭单硫化物 RaX 及其离子 RaO± 的理论研究，其重要性远超简单的分子性质计算，它触及了物理学和化学的多个前沿领域，并为未来的实验和理论探索打开了新的大门。

5.1.1 探测新物理：P、T 对称性破缺

• 增强效应的分子：重极性分子，特别是含有镭（Z=88）这种极重原子的分子，被认为是探测超越标准模型的物理现象（如电子电偶极矩（eEDM）、核磁四极矩矩、宇称（P）和时间反演（T）对称性破缺）的理想体系。在这些分子中，强烈的相对论效应导致内部有效电场极度增强（通常是原子内部电场的 10^4 到 10^5 倍），从而极大地放大任何潜在的 eEDM 信号。这使得重分子成为探索基本粒子物理学的重要探针。
• RaX 的独特性：镭单硫化物具有大的原子核电荷和极性化学键，这使得它们在理论上被预测具有显著的 eEDM 增强因子。本研究对这些分子偶极矩的准确计算，特别是 RaS 和 RaSe 超过 10 D 的极高偶极矩，进一步证实了其作为 eEDM 实验候选者的潜力。高偶极矩意味着分子可以通过外部电场进行有效地极化和操控，这对于 eEDM 实验中的自旋极化和去相干时间控制至关重要。因此，尽管论文没有直接计算 P、T 奇数参数，但其电子结构和电学性质的深入分析，为此类精确物理实验奠定了坚实的基础。

5.1.2 激光冷却的挑战与机遇

• 激光冷却标准：实现分子激光冷却的关键在于分子需要具有高度对角线的弗兰克-康登因子（FCFs），这意味着激发态主要衰变回基态的相同振动能级。此外，分子应具有短的激发态寿命和简单的能级结构，以便有效地循环激发和发射光子。
• RaX 的激光冷却困境：本研究发现 RaX 分子（RaO, RaS, RaSe）的 FCFs 矩阵高度非对角线（图 4），这与此前在镭单卤化物（如 RaF, RaCl）中观察到的对角线 FCFs 形成鲜明对比 [4, 5]。这种非对角性意味着 RaX 分子不适合直接的激光冷却，因为一个激发态会衰变到基态的多个振动能级，导致激光冷却所需的激光束数量过多，效率低下。
• 化学键合差异的根源：论文深入解释了这种差异的化学键合根源。镭单卤化物（如 RaF）的键合是单价的，其中镭的一个 s 电子转移到氟的 p 轨道，剩余的镭 s 电子仍强烈定域在镭原子附近，导致激发时键长变化不大，从而产生对角线 FCFs。而 RaX 分子则表现出二价键合特征，镭的两个 s 电子转移到硫族元素，形成极性共价键。这种二价键合意味着电子激发会引起键长的显著变化，导致 FCFs 的非对角性。这种深入的化学键合分析对于理解分子性质和指导实验具有重要意义。
• 其他操控手段：尽管不能直接激光冷却，RaX 分子极大的偶极矩可能使其适合通过外部电场进行操控和减速，例如利用 Stark 减速器，这为实现冷分子样本提供了另一种途径。

5.1.3 化学键合与相对论效应的新见解

• 二价键合的本质：该研究明确了镭单硫化物中的二价键合特征，这与碱土金属单硫化物的一般规律一致，但由于镭是重元素，其相对论效应在此类键合中扮演了关键角色。镭（碱土金属）倾向于失去两个电子形成 Ra²⁺ 离子，而硫族元素（O, S, Se）倾向于接收两个电子形成 X²⁻ 离子。因此，RaX 键合具有很强的离子性，但由于共享电子的存在，也具有共价性，形成了极性共价键。
• 分子轨道可视化：图 2 的分子轨道分析直观地展示了这种键合特征，例如 RaO 中电子密度高度定域在氧原子上，而 RaS 和 RaSe 中镭的 6pz 轨道对键合的贡献增加，表明了共价特征的增强。这有助于理解键合强度、键长和电学性质的差异。
• 相对论效应的深远影响：相对论效应不仅影响键长和激发能，还导致了自旋-轨道耦合效应，使 A-S 态分裂为 Ω 态。论文详细分析了这些影响，并指出 ECP 方法在平衡键长和振动频率方面已能很好地捕捉大部分相对论修正，尽管在激发能方面仍有改进空间。这些见解对于理解重元素化学、开发更准确的理论模型至关重要。

5.2 未来展望与实验方向

这项理论工作为未来的实验和理论研究指明了几个关键方向：

5.2.1 实验合成与表征

• 合成方案设想：论文提出了一个富有前景的 RaX 分子合成方案：通过激光烧蚀固体镭靶，并在含有硫族元素供体气体（如 O₂）的环境中进行。这一思路借鉴了 CaH 分子合成的成功经验（Ca + H₂ → CaH + H）[36]。尽管 Ra + O₂ → RaO + O 是放热反应，但通过对热镭原子进行振动激发，可能有助于 RaO 的形成。随后，可以通过氦缓冲气冷却这些分子，将其制备成冷分子样本。
• 光谱学验证：一旦 RaX 分子在实验中合成，高精度光谱学测量将是关键。理论预测的势能曲线、激发能、偶极矩和光谱常数将是识别和表征这些分子的重要依据。与理论预测值的对比将进一步验证量子化学方法的准确性。

5.2.2 深入的理论研究

• 全四分量相对论计算：虽然 ECP 和 X2C 方法取得了良好的平衡，但未来的研究可以探索更精确的全四分量相对论方法（如 DIRAC 或 ReSpect 软件包中的狄拉克-哈特里-福克或狄拉克-耦合簇），以消除赝势和 X2C 框架的近似性，为所有性质提供最终的基准精度。
• P、T 奇数参数计算：鉴于 RaX 分子在探测新物理方面的潜力，下一步关键的工作是直接计算这些体系中的 P、T 奇数参数，如电子电偶极矩有效场 Eeff 和核自旋-电子标量-赝标量耦合常数 W_s。这些计算将为实验提供更具体的指导。
• 偶极束缚态的详细研究：RaX 分子极大的偶极矩预示着偶极束缚态（DBS）的存在。未来的理论工作可以深入探究这些 DBS 的能量、波函数和寿命，这对于理解分子与电场的相互作用以及在外部场中操控分子具有重要意义。
• 精确解离能和长程相互作用：为了更全面地理解 RaX 分子的行为，未来的研究应扩展势能曲线的计算范围，精确解析其长程渐进行为，从而获得准确的解离能和长程相互作用势。这对于理解分子碰撞和超冷化学至关重要。
• 更广泛的激发态分析：除了最低的几个激发态，未来的研究可以探索更多的激发态，特别是在更高能量区域，以全面理解其电子结构和潜在的光谱特征，以及与不同解离极限的关联。
• 超精细结构和斯塔克效应：计算超精细相互作用常数（由核磁矩和电四极矩引起）和 Stark 效应（偶极矩与电场相互作用）将提供更详细的能级结构信息，对于精确光谱测量和外部场操控具有直接指导意义。

5.2.3 冷分子实验应用

• 外部场操控：利用 RaX 分子巨大的偶极矩，设计和实施基于 Stark 效应的外部电场操控实验，例如 Stark 减速，以制备冷分子束。这对于高精度测量和研究分子间的相互作用至关重要。
• 精密测量实验：将冷 RaX 分子应用于 eEDM 等精密测量实验，利用其对 P、T 奇数效应的增强敏感性，寻找超出标准模型的物理信号。这代表了将先进量子化学与基础物理学研究相结合的终极目标。

5.3 技术复杂性与方法论创新

该研究的另一个亮点在于其方法论的复杂性和创新性，这体现在对多种先进量子化学技术的整合和应用上：

• 多层次方法学：通过结合 MRCI+Q+ECP+SO（处理势能曲线和光谱学，侧重多参考和标量相对论与微扰 SO）与 CCSD(T)+X2C（处理高精度偶极矩和极化率，侧重全电子相对论和动态关联），以及 EOM-CCSD 和 TI-CCSD，研究团队展示了如何针对不同性质选择最优的计算策略。这种多层次的方法学确保了计算的效率和精度。
• 基准化与验证：论文在附录中详细介绍了其方法学在 RaF, SrF, BaF, BaO 等基准体系上的性能，通过与实验数据和更高级别的理论计算进行比较，严格验证了所用方法的可靠性。这种严谨的基准化过程是高质量量子化学研究的标志。
• 重原子处理：对镭这样极重元素的准确处理本身就是一项重大技术挑战。ECP 和 X2C 方法的结合应用，以及专门设计的相对论基组（如 aug-cc-pVQZ-X2C 和 dyall.avnZ），体现了在处理强相对论效应方面的专业知识。

这项工作不仅在学术上对镭单硫化物这一新兴分子家族提供了全面且深入的理解，更从方法论和应用层面为重元素分子在基础物理和化学前沿领域的未来研究，树立了一个典范。