地核中的轻氢同位素

《SCIENCE ADVANCES》：Light hydrogen isotopes in terrestrial core

　　

地球水的起源尚未明确。氢同位素为行星水历史提供了关键示踪，但地核——地球最大的内部储库——的同位素组成此前未受约束。在此，研究人员采用第一性原理计算结合机器学习加速路径积分分子动力学（ML-accelerated PIMD），量化了成核条件下硅酸盐熔体与金属熔体之间的氢同位素分馏。研究人员发现，核形成使硅酸盐地球（silicate Earth）富集氘（deuterium），同时将同位素轻的氢（isotopically light hydrogen）集中于地核，这要求原始地球（proto-Earth）的初始氘氢比（D/H）低于现今地幔所保存的值。研究人员的模型表明，全地球（bulk-Earth）D/H比可通过两种方式解释：直接吸积以顽火辉石球粒陨石（enstatite chondrite，EC）为主的物质，或通过富氘星子与太阳星云气体（solar nebular gas）相互作用发生同位素重置。这些结果表明，地球的水储量是在吸积最早阶段确立的。

**论文解读文章**

**研究背景与问题**
水是地球演化和长期宜居性的核心要素，影响地幔动力学、岩浆作用和板块构造。然而，地球水的起源仍存在争议：晚期增生模型（如碳质球粒陨石或彗星输运）、顽火辉石球粒陨石继承模型以及太阳星云气体直接吸积模型各有支持。氢同位素（D/H比）是追溯行星水历史的关键示踪剂，但地核作为地球最大的内部储库（可能储存75–95%的氢），其氢同位素组成此前从未被约束。核心缺失的参数是成核条件下（高压高温）硅酸盐熔体与金属熔体之间的平衡氢同位素分馏因子。实验难以达到这些条件，强核量子效应使理论预测复杂化。因此，研究人员开展本研究，旨在量化该分馏因子，进而约束全地球的D/H比，并甄别地球水的可能来源。论文发表在《SCIENCE ADVANCES》。

**研究工作与结论**
研究人员通过第一性原理密度泛函理论（DFT）结合机器学习加速路径积分分子动力学（PIMD），构建了高精度机器学习势（MLPs），计算了30–90 GPa、3000 K条件下硅酸盐熔体（如Mg₃₂Si₃₂O₁₀₆H₂₀、Mg₃₀Ca₂Fe₄Si₂₄Al₃NaO₉₂H₆）与金属熔体（Fe₁₀₀H₂₈、Fe₉₃Ni₅Si₆O₈C₃H₁₂）的平衡氢同位素分馏因子（10³lnα_{silicate-metal}）。结果表明，硅酸盐熔体相对富集氘，分馏因子在20–60 GPa、3000 K下为约21–26‰，显著低于此前间接估算的约50‰。基于此，研究人员模拟了核形成过程中氢同位素的演化，发现核心形成使硅酸盐地球的δD升高12–78‰（瑞利分馏）或12–25‰（平衡分馏），反推出原始地球的δD为?42至?188‰。这一数值可通过直接吸积顽火辉石球粒陨石（δD ≈ ?104‰）来解释，或通过富氘星子（如普通球粒陨石与碳质球粒陨石混合物，δD≈+168‰）与太阳星云气体（δD低至?850‰）发生同位素交换来实现。由此，地球的水储量在吸积最早阶段即已确立。

**关键技术方法**
研究人员采用了三种主要技术：（1）第一性原理分子动力学（FPMD）模拟，基于密度泛函理论（DFT）生成硅酸盐和金属熔体的构型数据集（30–90 GPa、3000 K）；（2）机器学习势（MLPs）训练，使用DeePMD-kit构建高精度势函数（能量RMSE 3.82–7.46 meV/atom，力RMSE 0.19–0.27 eV/?），以加速PIMD模拟；（3）路径积分分子动力学（PIMD）结合直接缩放坐标估计器（ZSC），计算平衡氢同位素分馏因子。样本来源为理论模拟的熔体体系，未涉及具体样本队列。

**研究结果**

**氢键合行为（Hydrogen bonding behaviors in melts）**
通过FPMD和PIMD分析熔体结构：硅酸盐熔体中H主要以羟基（OH）形式存在，H–O键长1.10–1.16 ?，压力下键长反常伸长（至1.16 ?），同时配位数（CN）增加（1.33→1.80），导致有效键强近似压力无关；金属熔体中H与Fe键合，H–Fe键长1.78–1.88 ?，压力下键长缩短、CN增加（6.53→6.63）。其他元素（Ca、Fe、Al、O、C、S）不显著改变氢键合行为。

**核形成过程中的氢同位素分馏（Hydrogen isotope fractionation during core formation）**
利用PIMD和MLPs计算了折合配分函数比（10³lnβ）和分馏因子（10³lnα_{silicate-metal}）。金属熔体的β因子随压力线性增大，硅酸盐熔体的β因子近乎压力无关。在3000 K、20–60 GPa下，10³lnα_{silicate-metal}从约26‰降至约21‰。氧逸度影响较小（约6‰差异）。研究人员将此分馏因子应用于单阶段和多阶段核形成模型，显示地核形成使硅酸盐地球δD升高12–78‰，地核δD范围为?46至?192‰。特别地，Baffin Island橄榄石熔融包裹体的低δD（?144±24‰）与预测的地核δD一致，支持核幔交换可能改写了深部地幔储库的氢同位素组成。

**讨论与结论**
研究结果识别出两条可行的地球水源路径：直接继承顽火辉石球粒陨石为主的物质，或通过太阳星云气体（δD = ?100至?850‰）与千米级熔融星子发生同位素交换，重置初始富氘混合物。后一机制与早期太阳系盘内落积过程的氢同位素分布一致，表明星云气体-星子相互作用可能是调控类地行星挥发分和同位素架构的关键普遍过程。研究结论指出，地球的水储量在吸积最早阶段（非晚期增生）即已确立，核心形成进一步分馏氢同位素，使硅酸盐地球富集氘，而地核保留轻氢同位素，最终形成了现今地幔的氢同位素特征。