行星冰体中碳酸的氢-碳双超离子输运通道

《SCIENCE ADVANCES》：Hydrogen-carbon doubly superionic conduits of carbonic acids in planetary ices

　　

近期天文观测显示，二氧化碳（CO2）在行星天体上广泛存在，并常与水（H2O）共存。由CO2-H2O相互作用产生的晶态碳酸（H2CO3）已在高压下被预测和合成，但其在行星内部条件下的动态行为仍知之甚少。本研究调查了0至500吉帕（GPa）压力范围内CO2-H2O行星冰的稳定性。研究人员的模拟表明，H2CO3和正碳酸（orthocarbonic acid, H4CO4）成为主要的碳-氧-氢（C─O─H）化合物，并随着压力增加从分子晶体演化为三维（3D）固体。特别地，这两种碳酸均转变为氢（H）扩散的超离子（superionic）相，随后转变为由相互连接的氧多面体空隙实现的碳-氢（C─H）双超离子通道，其中离子间的关联运动增强了离子电导率。此外，H2CO3中的氢表现出强烈的各向异性行为，这可能对冰巨行星的非轴对称磁场有贡献。这些发现进一步表明，碳酸冰可以维持氢-碳传输，可能增强巨行星内部对流挥发分的循环。

本研究《Hydrogen-carbon doubly superionic conduits of carbonic acids in planetary ices》发表在《SCIENCE ADVANCES》期刊，旨在探究行星内部高压高温极端条件下，由水和二氧化碳相互作用形成的碳酸化合物的结构演化、动态特性及其在行星物理化学过程中的作用。

研究背景与动机

氢、碳、氧是太阳系中最丰富且化学性质多样的元素，其组合形成H₂O、CO₂、CH₄等关键分子物种，通常被归类为行星冰或轻元素冰。这些冰不仅构成冰巨行星“填充冰”（filled ices）的主要成分，也是地球地幔和过渡带的重要挥发分。在行星内部高压高温条件下，这些小分子可相互作用形成复杂物相，例如水和水在超高压下可转变为超离子态和金属态。近年观测证实CO₂和H₂O在行星表面和内部的存在，并可能在行星动态特性中起关键作用。碳酸（H₂CO₃）作为H₂O和CO₂最简单的结合物，其高压下的晶体结构已有研究，但其在行星内部高温高压条件下的动态行为，特别是轻元素如碳和氢的输运性质，仍很大程度上未被探索，限制了对它们在行星内部对流过程中作用的理解。因此，阐明轻元素冰在极端条件下的结构、动态和输运特性，对于模拟行星（特别是冰巨行星）的形成、演化和磁场行为至关重要。

研究方法概述

为开展此项研究，研究人员主要运用了以下几项关键技术方法：首先，采用基于密度泛函理论（DFT）的第一性原理计算，通过维也纳从头算模拟软件包（VASP）实现，并使用optB86b泛函处理氢键或范德华相互作用。其次，利用粒子群优化技术，通过CALYPSO代码在CO₂-H₂O体系中进行了全面的晶体结构搜索，压力范围覆盖25、50、140和300 GPa。第三，使用从头算分子动力学（AIMD）模拟研究了稳定轻元素冰的动态性质，包括在宽泛温压条件下的离子扩散行为和相变。第四，通过电子定域函数（ELF）、集成晶体轨道汉密尔顿布居（ICOHP）、量子理论中的原子（QTAIM）拓扑分析等方法，对化学键演变、电子结构和扩散机制进行了深入分析。第五，基于格林-久保（Green-Kubo）公式和能斯特-爱因斯坦（Nernst-Einstein）关系，计算了离子的扩散系数和离子电导率，并分解了离子间的关联运动对电导率的贡献。

研究结果与发现

RESULTS（结果）

• •
  相图与结构演化：研究人员在CO₂-H₂O体系中进行了全面的晶体结构搜索。结果显示，在广阔的压强范围内，只有H₂CO₃在形成能上是稳定的。H₂CO₃经历了一系列结构转变：在低压下结晶为具有顺-顺（cis-cis）构象分子的Pnma相；在约36 GPa以上，聚合的Cmc2₁相出现并保持稳定；在233 GPa时，Cmc2₁相内部发生同构转变，形成三维（3D）结构。此外，研究还揭示了两个先前未被认识的亚稳态同质异形体（空间群P2₁和Fdd2），在32-35 GPa范围内与稳定相能量相近。在压力高于320 GPa时，H₂CO₃可进一步与H₂O反应形成正碳酸I4₁/a-H₄CO₄。H₂CO₃的晶体结构表现出独特的维度交叉行为：从分子晶体到一维链，再到二维（2D）平面，最终在高压下形成完全三维的聚合网络。这一过程伴随着氢键的对称化和化学键的演变。
• •
  多种超离子态：AIMD模拟揭示了H₂CO₃在高温下的一系列不同超离子态。在100 GPa条件下，随着温度从500 K升至4000 K，体系从固态（所有原子在平衡位置振荡）先后转变为氢超离子态（氢扩散，碳氧固定）和氢-碳双超离子态（氢和碳均扩散，氧骨架保持）。在更高温度下，体系最终熔化。研究特别识别出两种氢超离子态：一种是准二维超离子态，氢离子仅限制在特定的bc平面内扩散，形成受限的二维通道；另一种是完全三维氢超离子态，氢离子在整个晶格中扩散，但在中等温度下仍表现出显著的各向异性行为。径向分布函数（RDF）进一步区分了这些相变。
• •
  相图与扩散机制：研究人员将压强范围扩展到30至500 GPa，构建了H₂CO₃的相图。根据氢、碳、氧原子的熔融行为及其维度相关的扩散特征，相图可分为五个区域：固态、二维氢超离子态（SI-I）、三维氢超离子态（SI-II）、氢-碳双超离子态（SI-III）和熔融态。随着压力增加，超离子态的起始温度逐渐升高，而熔化温度则急剧上升，导致高压下超离子态的温度窗口变宽。H₄CO₄在高温下也经历了多种超离子转变，但其扩散不表现出强各向异性。对原子空间概率分布的分析阐明了超离子性出现的原子机制。在二维氢超离子态中，氢密度被限制在层间区域，形成明确的二维扩散平面，氧原子形成的四面体和八面体空隙促进了面内氢扩散。在三维氢超离子态中，热能为氢克服链间势垒提供了足够能量，O₆八面体形成的空隙提供了额外的面外扩散通道。在双超离子态中，碳不仅在其原始CO₄四面体位点附近涨落，还占据氧骨架内相邻的O₄四面体和O₆八面体空隙，这些相互连接的多面体空隙为氢和碳的扩散创造了通道。
• •
  离子电导率：H₂CO₃中的双超离子通道对其离子电导率有贡献。氢的扩散系数遵循阿伦尼乌斯（Arrhenius）关系，约为~10^?4cm²/s，与其他超离子导体报道的值接近。碳的扩散系数远低于氢，但与铁-碳（Fe-C）合金中的值相当。基于格林-久保公式计算的离子电导率（σ_GK）在~10至~150 S/cm之间，与超离子H₂O和NH₃相的电导率相当。通过分解电导率贡献发现，除氢-碳（H-C）对的交叉关联项为较小的负值（表明碳氢离子的反关联运动降低了电导率）外，其他所有交叉关联贡献均为正，意味着这些离子物种间的关联运动增强了电荷输运。对于氢离子，其自身不同离子间的关联贡献（σ^distinct）为正，进一步表明氢离子间的关联对电荷输运有正面贡献。这些关联扩散效应使得由格林-久保公式得到的离子电导率高于简化的能斯特-爱因斯坦关系计算值。

讨论与结论

在讨论部分，研究人员将H₂CO₃的相图与H₂O的相图进行了比较。虽然H₂CO₃的密度比冰大，但其超离子转变的起始温度通常低于被充分研究的H₂O相。这种差异可以从化学键角度解释：在所研究的压力范围内，H₂CO₃中H─O键的ICOHP值持续低于H₂O，表明H₂CO₃中的H─O键更弱，因此需要更少的热能来激活氢扩散，从而降低了超离子转变温度。

研究指出了H₂CO₃和H₄CO₄在行星内部碳循环中的潜在作用。先前提出的“钻石雨”假说预测，在冰巨行星深部，甲烷等还原碳会分解为更氧化的钻石和氢，导致碳沉降至核幔边界（CMB）。研究人员通过第一性原理计算检验了化学反应活性，提出在高于240 GPa的压力下，通过反应4Fe + C + 3H₂O = H₂CO₃+ 4FeH和4Fe + C + 4H₂O = H₄CO₄+ 4FeH（Fe来自核心，H₂O来自地幔）形成H₂CO₃和H₄CO₄在能量上是有利的。在此情景下，降至冰巨行星核幔边界的钻石可能通过形成氧化的SI-III H₂CO₃-H₄CO₄混合物而被活化。这表明冰巨行星内部碳的化学状态丰富，并可能通过C─H超离子相在核幔边界发生活跃的碳交换。

此外，H₂CO₃和H₄CO₄超离子相的存在也可能对冰行星磁场的生成有贡献，因为高于~10 S/cm的电导率通常被认为足以维持行星磁场。特别是，由于H₂CO₃固有的层状结构，即使在三维扩散区域，离子扩散也具有高度各向异性，导致平行或垂直于平面的离子电导率值不相等。这种方向性电导率可能诱发双分量磁场，从而可能对在天王星和海王星观测到的非偶极和非轴对称磁场形态有所贡献。

研究结论翻译：

本研究表明，H₂CO₃和H₄CO₄是CO₂-H₂O体系中热力学稳定的化合物，在高压下经历从分子晶体到三维固体的结构演化。这些碳酸在高温下转变为氢扩散超离子相，并进一步演化为氢-碳双超离子通道，其中相互连接的氧多面体空隙实现了离子的关联输运，显著增强了离子电导率。H₂CO₃中的氢扩散表现出强烈的各向异性。这些发现表明，碳酸冰能够在巨行星内部维持氢和碳的传输，可能通过促进挥发分的对流循环，并对冰巨行星的非轴对称磁场产生潜在影响，从而深化了对行星内部动态过程和磁场起源的理解。