Superatomic hydrogen: achieving effective aggregation of hydrogen atoms at pressures lower than that of metallic hydrogen

技术摘要：超原子氢

问题陈述
金属氢以其电子离域和分子键断裂的特性，在理论上被公认为受控核聚变和高温超导性的潜在介质。然而，合成体相金属氢需要超过 500 GPa 的极高压力，这一条件在实验上仍难以实现，并构成了实际应用的重大障碍。虽然类金属氢化物（如以镧为中心的结构、$\text{CaH}_6$、$\text{CeH}_9$）可以在较低压力（数十 GPa）下合成，但它们包含金属原子，因此不适用于核聚变应用，因为核聚变需要纯元素氢。核心挑战在于寻找一种路径，能够在显著降低压力的条件下，实现纯元素氢的有效聚集和类金属特性。

研究方法
作者提出了一种基于氢笼结构（$\text{H}_{13}$）的“超原子氢”模型，该结构由一个中心氢原子和由 12 个氢原子组成的壳层（$\text{H}_{12}$）组成，模拟了稳定的碱金属簇（$\text{M}_{13}$，其中 $\text{M} = \text{Li, Na, K}$）的二十面体（$I_h$）对称性。

本研究采用高精度从头算（ab initio）计算来研究各向同性压缩下 $\text{H}_{13}$ 系统的结构演化和电子特性：

1. 结构压缩： 中心原子与 $\text{H}_{12}$ 笼之间的径向距离（$r$）从 4.0 Å 系统地减小至 0.9 Å。
2. 电子结构分析：
   • 多态完全活性空间二阶微扰理论 (MS-CASPT2)： 用于扫描势能曲线并确定不同压缩范围内的基态自旋多重度。
   • 耦合簇单双激发及微扰三激发理论 (CCSD(T))： 用于对最低能量自旋态进行单点能计算，以获得可靠的分子轨道（MO）能级。
   • 径向分布函数 (RDF) 和电子密度： 用于追踪电子的空间分布以及从定域化到离域化的转变。
   • 定域轨道定位器 (LOL)： 用于进一步表征电子定域化特征。
3. 压力估算： 外部压力（$P$）通过公式 $P = |dE/dV|$ 进行估算，其中 $dE$ 是总能量变化，$dV$ 是范德华体积变化（通过 Bader 方法计算）。
4. 对比模型： 研究还分析了椭球形 $\text{H}_{13}$、球形 $\text{H}_{12}$ 和椭球形 $\text{H}_{12}$ 系统，以评估对称性破缺和缺失中心原子所带来的影响。

关键结果

• 超原子态出现： $\text{H}_{13}$ 系统在径向距离约为 1.9 Å 时转变为超原子态。在此状态下，电子构型遵循超原子壳层填充规则（例如，在 1.3 Å 处为 $1S^2 1P^6 1D^5$，在 0 9 Å 处为 $1S^2 1P^6 2S^1 1D^4$）。
• 电子离域： 原子级分析表明，在压缩的超原子态中，中心氢原子将其电子捐赠给周围的笼层。因此，所有 13 个电子都离域到超原子分子轨道（SAMOs）上，表现出类似于金属氢的特性。
• 压力降低： 实现该超原子态所需的压力显著低于体相金属氢。具体而言，$\text{H}_{13}$ 系统在 1.9 Å 半径处达到超原子态所需的压力计算为 6.2 GPa。即使在 1.3 Å 半径处，压力也仅为 82.7 GPa。这些数值大约比通常引用的 ~500 GPa 金属氢阈值低两个数量级。
• 对变形的鲁棒性： 对低对称性模型（椭球形 $\text{H}_{13}$ 和 $\text{H}_{12}$）的分析表明，尽管具体的压力阈值有所不同，但这些系统也在远低于 500 GPa 的压力下实现了超原子态。例如，椭球形 $\text{H}_{13}$ 系统在 1.7 Å 处于 2.1 GPa 时即可达到超原子态。

意义与主张
本文声称提供了一个理论框架，用于在远低于传统金属氢所需压力的情况下实现元素氢的有效聚集。通过利用超原子的概念（即由原子簇模拟单个原子的电子壳层结构），作者证明了氢笼（$\text{H}_{13}$）可以在 ~6.2 GPa 下表现出类金属的电子离域特性。

作者认为，这一发现“原则上为超原子氢在核聚变等领域的应用开辟了前景”，为目前被认为必需的极端条件提供了一种潜在的替代方案。该研究表明，构建基于超原子态的氢笼可能是降低合成金属氢难度的可行策略，从而推动能源和材料科学的应用。研究结论指出，尽管构建此类空间限制系统仍是一个科学挑战，但超原子视角为未来的研究提供了一条充满希望的途径。