来源论文: https://arxiv.org/abs/1911.06300 生成时间: Mar 01, 2026 10:40

0. 执行摘要

在现代超冷原子与离子混合实验的版图中，锶（Strontium）因其丰富的能级结构和优异的光晶格特性而备受关注。然而，作为冷原子物理研究中的核心体系之一，锶二聚体阳离子（Sr₂⁺）的电子结构在很长一段时间内缺乏高精度的理论描述。波兰华沙大学的 Michal Śmiałkowski、Tatiana Korona 和 Michal Tomza 在 2019 年发表的这项研究（arXiv:1911.06300v1），填补了这一空白。

本研究利用高精度的从头算（ab initio）量子化学方法，包括多参考配置相互作用（MRCISD）和耦合簇（CCSD(T)/CCSDT）方法，系统地计算了 Sr₂⁺ 分子离子的基态及 41 个激发表态的势能曲线（PECs）。研究揭示了该体系虽然仅有三个价电子，但由于高能态密度和显著的多参考特性，其电子结构的精确描述极具挑战性。核心成果包括：

  1. 确定了 42 个电子态的势能曲线及相关光谱常数。
  2. 发现并在物理上解释了 $1^2\Sigma_g^+$ 态的双阱结构。
  3. 计算了各电子态间的电偶极跃迁矩（EDTMs）。
  4. 分析了 Sr₂ 分子光电离产生 Sr₂⁺ 的 Franck-Condon 因子，为冷分子离子实验制备提供了直接指导。

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

1.1 核心科学问题

Sr₂⁺ 体系的科学重要性主要体现在冷混合实验中。当激光冷却的 Sr⁺ 离子被置于超冷的 Sr 原子云中，或者通过 Sr₂ 分子的光电离，会产生 Sr₂⁺ 分子离子。理解其电子结构是控制受控化学反应、高精度光谱测量以及量子模拟的前提。然而，碱土金属分子离子由于其电子态密度极大，且在短程区域存在强烈的多参考效应（Multireference nature），传统的单参考方法（如常规 CCSD）往往在描述激发态或特定平衡位置时失效。

1.2 理论基础：Born-Oppenheimer 与 相对论效应

研究在 Born-Oppenheimer 近似框架下展开。对于锶（原子序数 38）这种中等质量的元素，标量相对论效应不可忽略。作者采用了 ECP28MDF 小芯（Small-core）相对论能量一致伪势（Pseudopotential），将 28 个内层电子用有效核心势替代，从而显著降低了计算复杂性，同时捕捉了核心电子的相对论效应。剩下的 19 个电子（Sr 和 Sr⁺ 的价电子及次价电子层 $4s^2 4p^6 5s^2/5s^1$）被显式地进行相关计算。

1.3 技术难点：多参考特性的挑战

Sr₂⁺ 虽然名义上只有三个价电子（相对于两个 Sr 原子原本的四个价电子），但其激发态涉及到 $5p$、$4d$ 轨道的大量参与。这些轨道能量接近，导致波函数表现出极强的构型混合。作者指出，在 $1^2\Sigma_g^+$ 态中，短程区域的稳定性源于特定的电子激发构型，这种构型在 CCSD 方法中被严重低估，必须使用多参考方法才能正确还原双阱（Double-well）结构。

1.4 方法细节：计算流水线

  1. 轨道优化：使用多构型自洽场（MCSCF/CASSCF）方法。活性空间包含两个原子的 $5s, 5p, 5d, 6s$ 轨道构成的全分子轨道，确保了激发态计算的平衡性。
  2. 多参考相关：采用 MRCISD（多参考单双激发配置相互作用）。为了修正尺寸一致性（Size-consistency）误差，应用了 Davidson 修正（MRCISD+Q）。
  3. 单参考基准：对于基态及最低激发态，使用了 RCCSD(T) 和 RCCSDT（显式包含三激发项）进行标定，验证多参考方法的准确性。
  4. 基组设置：采用了极大的 [14s11p6d5f4g] 高斯基组，并额外添加了 [3s3p2d2f1g] 键函数（Bond functions），这在描述极化率和长程色散力方面至关重要。

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

2.1 原子基准测试 (Table I)

为了验证方法的可靠性，作者首先计算了单原子 Sr 和 Sr⁺ 的能级：

  • 电离能 (IP)：CCSD(T) 计算值为 45820 cm⁻¹，与实验值 45932 cm⁻¹ 误差小于 0.25%。
  • 激发能：对于 Sr(¹S → ³P) 和 Sr⁺(²S → ²D) 的预测，CCSD(T) 的误差控制在 0.5% 到 5% 之间。MRCISD 误差略大（约 3.8%），这反映了该方法在处理原子激发时相较于耦合簇的局限性，但在分子态计算中，其多参考优势弥补了这一点。

2.2 分子基态 $X^2\Sigma_u^+$ (Table III)

  • 平衡键长 $R_e$：MRCISD 计算值为 7.99 bohr，与实验推导值的一致性极高。
  • 阱深 $D_e$：MRCISD+Q 给出的解离能为 8801 cm⁻¹，完美契合了实验值 8800(130) cm⁻¹。这一结果证明了计算方案在处理该体系时的极高精度。

2.3 异常的 $1^2\Sigma_g^+$ 态分析

这是本文最引人注目的发现之一。该态呈现典型的双阱特征：

  • 短程极小值：位于 $R \approx 9.2$ bohr，深度约 678 cm⁻¹。这由多参考构型交互作用稳定。
  • 长程极小值：位于 $R \approx 17.0$ bohr，深度约 184 cm⁻¹。主要受极化感应相互作用控制。
  • 物理机制：作者通过分析激发振幅发现，短程阱的形成主要源于从最高占据轨道 $\sigma_u^*$ 到最低未占据轨道 $\sigma_g$ 的电子激发。CCSD 方法由于低估了这类激发，无法预测出短程阱，导致在早期文献中该态常被误认为是完全排斥的。

2.4 势能曲线图 (Fig 2-6)

研究详细绘制了按解离限分类的 42 个电子态 PECs。这些曲线展示了极其复杂的**避免交叉（Avoided crossings）**现象，特别是在高能区域，这预示着强烈的非绝热耦合效应，对冷碰撞动力学有重要影响。

3.1 核心软件包：MOLPRO

本项工作的所有电子结构计算均使用 MOLPRO 程序包完成（版本 2012.1）。MOLPRO 是处理高度相关体系和执行高精度多参考计算的工业标准工具。

3.2 复现关键参数

  • 积分处理：采用 IntOpt 控制，并显式包含 ECP 积分。
  • CASSCF 步骤：在 $D_{2h}$ 点群对称性下执行。必须为每个不可约表示（$A_g, B_{3u}, B_{2u}, B_{1g}$ 等）设置正确的轨道占据数，以防止轨道塌陷。
  • MRCISD 步骤：在轨道优化后，使用 MRCI 模块。注意需要开启 closed 轨道（即不参与激发的内层占据轨道）和 occ 轨道（参与激发的核心轨道）的明确定义。
  • 基组输入：由于采用了自定义的 [14s11p6d5f4g] 基组，复现者需要参考作者在参考文献 [97, 119] 中提到的基组定义文件。键函数（Bond functions）通常放置在原子间距的中点。

3.3 脚本化与自动化建议

由于涉及 42 个电子态和大量 $R$ 点（通常从 5 bohr 到 50 bohr），复现此项工作建议：

  1. 编写 Shell 或 Python 脚本自动化生成多点 MOLPRO 输入文件。
  2. 采用 Restart 功能。先从长程 $R$ 点计算，利用收敛的波函数作为短程计算的初猜，以应对激发态收敛困难的问题。
  3. 态平均（State-averaged）CASSCF：在同一个对称性下计算多个激发态时，必须使用态平均权重，否则 PEC 曲线会出现不连续点。

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

4.1 关键引用

  • [111] Dugourd et al. (1992): 提供了早期的 Sr₂⁺ 实验解离能数据，是本文最主要的实验参考。
  • [128] ECP28MDF (Stuttgart/Cologne): 本文相对论效应处理的基石。
  • [124] Boys-Bernardi Counterpoise Correction: 用于消除基组重叠误差（BSSE）的标准方法，确保了长程势能的准确性。

4.2 局限性分析与评论

尽管本研究达到了极高的精度，但从现代量子动力学角度看，仍存在以下局限：

  1. 标量相对论与自旋轨道耦合（SOC）：研究主要考虑标量相对论。虽然碱土金属离子的 SOC 相对较小，但在处理精细结构分裂和特定的非绝热跃迁时，缺乏 SOC 可能会导致 PEC 交叉处的描述不够完美。
  2. 尺寸不一致性问题：即使有 Davidson 修正，MRCI 依然存在尺寸不一致性的固有缺陷。在解离限附近，PECs 的平滑度可能受到微小数值噪声的影响，这在 Table II 的平均激发能误差中有所体现。
  3. 高能态密度：在第七个解离限以上，态密度变得极其恐怖，MRCISD 的收敛性变差。这意味着更高能量的冷碰撞过程（如高能光致解离）超出了本文的预测范围。

5. 其他补充：冷原子应用前景

5.1 Franck-Condon 因子 (Fig 8)

本文计算了从 Sr₂($X^1\Sigma_g^+$) 基态到 Sr₂⁺($X^2\Sigma_u^+$) 和 $1^2\Sigma_g^+$ 的 Franck-Condon 因子。结果显示：

  • 基态到基态的跃迁具有强烈的对角分布（对角线上值最大），意味着通过光电离可以高效产生振动基态的 Sr₂⁺。
  • 到 $1^2\Sigma_g^+$ 态的跃迁由于双阱结构，FC 因子分布更广，这为实验制备高振动激发态或特定阱中的分子离子提供了能量地图。

5.2 跃迁偶极矩 (Fig 7)

EDTM 的计算揭示了某些激发态之间存在极强的跃迁强度（高达 9 Debye）。这些数据对于设计激光冷却循环（如果存在闭合回路）或执行受激拉曼断热通路（STIRAP）制备超冷离子至关重要。

5.3 总结

Tomza 团队的这项工作不仅是重元素碱土金属体系从头算的杰作，更是连接计算化学与实验室冷原子物理的桥梁。对于从事冷分子研究的团队来说，这份详尽的“能级说明书”是不可多得的参考资料。