Self-Consistent Evolution Models Show Weak Double-Diffusive Mixing in Jupiter and Saturn

技术摘要：自洽演化模型显示木星与土星中的双扩散混合作用微弱

问题陈述
巨行星的内部结构，特别是木星和土星中存在“模糊”或扩展核心的问题，仍然是行星科学领域的一个重大挑战。虽然双扩散对流已被提议作为将重元素从致密核心重新分配到上方包层中的一种机制，但其在现实演化背景下的效率仍具有不确定性。以往的研究依赖于理想化的成分结构，并且未能将双扩散传输自洽地纳入行星演化计算中。此外，局部流体动力学模拟（通常预测快速的层级合并和高效混合）与必须符合现今观测约束的全球演化模型之间存在差异。

方法论
作者在行星演化代码 APPLE 中实现了一种对双扩散传输的自洽处理。

• 初始条件： 研究采用了木星（基于 Stevenson et al. 2022）和土星（基于 Bodenheimer et al. 2025）的形成后内部模型。这些模型具有扩展的成分梯度（“模糊核心”），这是由形成过程中吸积的小行星发生部分瓦解和蒸发导致的，而非由于尖锐的核心边界。
• 传输物理学： 作者引入了由 Spruit (2013) 推导出的分层对流处方。该方法将跨越对流阶梯的传输建模为一种增强的扩散过程。有效成分扩散率 ($D_{\text{eff}}$) 和热扩散率 ($\kappa_{T,\text{eff}}$) 是基于传输独立于层厚度的渐近机制推导出的，并随稳定性比 ($R_\rho$) 和微观扩散率进行缩放。
• 演化模拟： 模型演化了 4.56 Gyr（太阳系年龄）。代码使用拉格朗日质量网格同时求解热演化和成分传输。对于不稳定于翻转对流（Ledoux 判据）的区域，采用标准的混合长度理论进行处理；而对于满足 $1 < R_\rho < \tau^{-1/2}$ 的 Ledoux 稳定区域，则利用双扩散处方进行处理。
• 敏感性分析： 为了测试结果的鲁棒性，作者运行了模拟，将有效扩散率相对于基准 Spruit (2013) 值分别提高了 100 倍和 1000 倍。

关键结果

• 混合效率低下： 模拟表明，双扩散对流在深层内部与包层之间产生的混合作用非常有限。在整个 4.56 Gyr 的冷却历史中，木星和土星中重新分配的重物质质量均不足 $\sim 1 M_\oplus$。
  • 在基准模型中，仅有约 $\approx 0.11 M_\oplus$（木星）和 $\approx 0.07 M_\oplus$（土星）被重新分配。
  • 即使扩散率提高了 1000 倍，重新分配的质量仍然较少（木星约为 $\approx 2 M_\oplus$，土星约为 $\approx 1.3 M_\oplus$）。
• 自我限制机制： 这种低效性源于一个负反馈循环。增强的热传输会改变热分层，从而降低稳定性比 ($R_\rho$)，进而抑制有效扩散率。这防止了失控混合，确保了即使是极大的传输系数也无法使原始成分梯度趋于均匀。
• 能量约束： 作者确定了一个根本性的能量约束。在许多流体动力学模拟所使用的 Boussinesq 近似中，由于层厚度远小于热尺度高度，混合在能量上是廉价的。然而，在可压缩的巨行星内部，稳定区域占据了行星半径的显著比例（$H \sim 0.25\text{--}0.5 R_p$），这使得混合成分梯度所需的能量与可用的热能预算相当。因此，驱动成分传输的浮力功受到严重限制。

意义与主张
本文声称，当在全局演化背景下进行自洽处理时，双扩散对流无法显著侵蚀行星形成过程中产生的成分梯度。观察到的木星和土星重元素分布可能需要额外的传输机制或形成路径，以产生比标准吸积模型预测的更广泛的初始混合。

作者认为，“模糊”核心可能是形成过程本身的直接产物（例如，形成早期更广泛的消融混合），或者是大型碰撞事件的结果，而非形成后双扩散混合的结果。该研究强调了将局部流体动力学结果（通常基于 Boussinesq 近似）应用于巨行星全球演化时的关键局限性，并强调深层内部的能量预算是调节混合效率的主要因素。