IF50.5｜亚海王星与超级地球的核幔混溶性

硅酸盐与氢在核幔界面完全混溶，颠覆传统行星结构模型。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

亚海王星和超级地球是银河系中最常见的行星类型，半径介于地球和海王星之间，但太阳系中没有对应物。通常认为超级地球具有岩石成分，而亚海王星由岩石核心和富氢包层组成。在核幔界面，压力超过数GPa、温度超过数千K，远超硅酸盐熔点，核心形成岩浆海，包层为氢流体。此前实验表明氢在硅酸盐熔体中有一定溶解度，但条件远不如行星内部极端。理论模型多假设氢和硅酸盐不混溶或溶解度有限，这直接影响对行星结构、演化和大气成分的预测。

科学问题

卡在哪里： 此前对氢-硅酸盐体系的研究要么基于低压实验外推，要么假设理想混合或惰性气体溶解度模型，无法描述高压高温下的真实化学相互作用。核幔界面条件下氢与硅酸盐是否发生大规模反应、能否完全混溶，以及混溶的临界条件，均缺乏第一性原理层面的定量认识。

本文要回答： 利用第一性原理分子动力学模拟，确定氢流体与MgSiO3熔体在核幔界面条件下的相平衡，揭示二者混溶的临界温度、压力依赖性和微观化学机制。

技术路线

作者主要用了：第一性原理分子动力学、相共存法、密度泛函理论、径向分布函数、2PT-MF。

作者采用相共存法，构建由纯氢和纯MgSiO3流体组成的初始双层体系，在DFT精度下进行分子动力学模拟，观察两相自发反应和平衡。通过一维密度剖面和三维粗粒化密度场定量提取两相组成，绘制相图并确定临界温度。同时进行均相模拟，计算混合焓、熵和体积，拟合非对称规则溶液模型以解释混溶行为。此外，通过物种分析和径向分布函数揭示氢与硅酸盐之间的氧化还原反应产物。

图 1　Fig. 1。图1展示相共存模拟的初始和平衡快照，以及一维密度剖面。初始时氢和硅酸盐完全分离，平衡后形成两个组成不同的相，界面清晰。密度剖面拟合双曲正切函数，给出两相组成和界面宽度。该图直观证明相分离存在，并为后续定量分析奠定基础。

核心亮点

亮点 1｜氢与硅酸盐完全混溶

在2–10 GPa、2500–4000 K范围内，相共存模拟显示氢流体与MgSiO3熔体形成两个平衡相：富氢相和贫氢相。随温度升高，两相组成逐渐接近，在临界温度TC以上合并为单一均相流体。TC随压力升高而降低，从2 GPa时的3500 K降至10 GPa时的2600 K。这一结果与以往基于亨利定律或理想混合的溶解度估计截然不同，表明在典型亚海王星核幔界面条件下，氢和硅酸盐可以任意比例混合。

图 2　Fig. 2。图2a展示两相组成随温度的变化，红色、灰色和蓝色符号分别对应2、5和10 GPa。随温度升高，两相组成相互靠近，在临界温度合并。曲线为规则溶液模型拟合结果，阴影为不确定性。图2b显示5 GPa、4000 K下混合焓和熵随组成的变化，焓为正、熵显著大于理想混合熵，解释混溶的热力学起源。

亮点 2｜混溶源于强非理想混合

均相模拟显示，混合焓为正（不利于混合），但超额混合熵远大于理想混合熵，导致溶线曲线顶部平坦，即两相组成在冷却时迅速偏离。混合体积为负，使得TC随压力升高而降低。非对称规则溶液模型拟合得到参数A=75 kJ/mol、B=424 kJ/mol、τ=4400 K、π=-35 GPa，成功再现相平衡和临界曲线。这种强非理想性源于氢与硅酸盐之间的化学反应，而非简单的物理溶解。

图 3　Fig. 3。图3为压力-温度相图，黑点表示两相共存条件，白点表示均相条件，黑线为临界曲线TC(P)，阴影为不确定性。虚线为MgSiO3液相线。图中标出典型亚海王星核幔边界条件，显示其位于临界曲线之上，表明完全混溶。该图将模拟结果与行星内部条件直接关联。

亮点 3｜氢与硅酸盐发生氧化还原反应

物种分析显示，在混合过程中H2减少，生成OH、H2O、SiO和SiHn等物种。反应可概括为SiO2 + H2 → SiO + H2O 和 SiO2 + 3H2 → SiH4 + 2OH，即硅被还原、氢被氧化。径向分布函数证实O-H和Si-H键的存在。在富氢相中，硅烷（SiH4）比低配位SiHn更丰富。这些反应产物可能被JWST观测到，为解释亚海王星大气中的水富集和潜在硅物种提供内生来源。

图 4　Fig. 4。图4展示5 GPa、4000 K下物种演化及均相流体的径向分布函数。插图显示H2比例下降，OH、H2O、SiO和SiHn比例上升。主图显示O-H和Si-H径向分布函数在短距离处有明显峰，证实化学键形成。该图揭示混溶的微观机制是氢与硅酸盐之间的氧化还原反应。

亮点 4｜行星演化模型需考虑核幔混溶

典型亚海王星核幔边界条件（如4.5 GPa、3600 K）位于临界曲线之上，处于完全混溶区。这意味着传统模型中假设核幔化学不反应是不合理的。在吸积阶段、沸腾丢失阶段和冷却阶段，氢在核与幔之间的分配会随温度压力变化而改变，影响包层质量、行星半径和大气成分。对于质量小于4 M⊕、年龄大于10 Gyr或包层质量分数小于0.02的亚海王星，可能不满足混溶条件，但大多数已观测行星处于混溶区。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 第一性原理分子动力学、相共存法、密度泛函理论、径向分布函数、2PT-MF
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 研究系外行星内部结构、大气演化和高温高压材料性质的行星科学家和计算材料学者。

带走一句

行星内部模拟必须考虑氢与硅酸盐的化学反应和完全混溶，否则会严重高估包层质量、低估核中氢含量，影响对行星演化和大气观测的解释。