Resolving core-envelope degeneracies in giant planets and sub-Neptunes: constraining the equivalence in the presence of dilute gradients

技術要約：巨大惑星およびサブ・ネプチューンにおけるコア・エンベロープの縮退の解決

問題提起
巨大惑星およびサブ・ネプチューンの内部構造は、惑星科学における高度に縮退した多次元の問題である。「ファジー（曖昧）」または希薄なコア（重元素が鋭い境界ではなく勾配として分布している状態）や、大規模な揮発性マントルの存在の認識は、観測された惑星半径（$R_p$）および質量（$M_p$）の解釈を複雑にしてきた。半径の谷（radius valley）を架橋するには、厚いH/He大気を持つ惑星、超臨界水海洋を持つ惑星、あるいは混合組成を持つ惑星を区別する必要がある。現在の進化モデル（MESA、APPLE、CEPAMなどのヘニイ型コード）は計算コストが高いため、静的な構造の縮退（特に組成勾配の幅 $\sigma$ に関して）の全空間をマッピングするために必要な高密度なパラメータ探索をしばしば阻害している。重元素の具体的な空間分布は、全半径および重力モーメント（$J_n$）に大きく影響するが、これらの勾配の関数的トポロジーは観測のみでは十分に制約されていない。

手法
著者らは、内部構造の縮退を高速にフォワードモデル探索するために設計された、オープンソースの静的なマルチゾーン静水圧構造積分器である fuzzycore を提示する。進化型コードとは異なり、fuzzycore は構造プロファイリングを、時間依存のエネルギー輸送、対流混合、および相平衡のクリティカルカーブから切り離している。

主な手法構成要素は以下の通りである：

• マルチゾーン静水圧平衡： 本フレームワークは、質量保存則と静水圧平衡に関する1次元球対称方程式を解く。界面における急峻な密度勾配を解像するために、グリッド解像度を動的に適応させる。
• 状態方程式（EOS）の定式化：
  • 混合コア： 鉄と多相岩石（高圧下での $\text{MgO}$ と $\text{SiO}_2$ への分離を含む $\text{MgSiO}_3$）の ab initio テーブルを加法的な体積則（アマガットの法則）を用いて混合し、単調な密度層状構造を確保する。
  • 流体エンベロープ： エンベロープを（H, He, 重元素）の三成分混合物としてモデル化し、比体積の理想的な線形混合を用いる。重元素（$Z$）は、原始的な大気から富化した基底層への遷移における平均分子量（$\mu$）の影響と圧縮率性を捉えるために、広範な範囲の $\text{H}_2\text{O}$ EOS によってプロキシされる。
• 凝縮マントル： 揮発性に富むサブ・ネプチューンのための、明確に相分離した水マントル（Ice VII/X またはスーパーイオニック状態）をサポートし、これをエンベロープの底部と熱力学的に結合させる。
• 組成勾配： 離散的な境界の代わりに、重元素質量分率 $Z(x)$ を連続的なガウス減衰関数としてモデル化する。勾配が鋭くなる（$\sigma \to 0$）際の数値的安定性を処理するため、フレームワークは解析的なサブグリッド平均を採用し、重元素の収支を保存し、非物理的な密度のスパイクを防ぐ。
• 数値積分： エンベロープは、完全に混合された対流ゾーンの「階段状」に離散化される。3D位相空間（$P, T, S$）の KD-Tree 上で局所的な線形最小二乗回帰を用いた堅牢な断熱ステップ法により、熱力学的状態を積分する。
• 境界シューティング： 境界値問題を解くために、入れ子状の境界シューティング法（ブレント法）を用い、質量、圧力、および温度の連続性を一致させるように界面圧力を反復計算する。これは、質量ターゲットモードおよび表面重力ターゲットモードの両方をサポートする。
• 安定性ガード： 低圧相転移における数値的な冷却アーティファクトを防ぐための「等温ガード」と、極めて薄いエンベロープのための適応型ダウンサンプリングを含む。

主な貢献および結果

1. 関数形式の堅牢性： 進化由来の木星プロファイル（Vazan et al. 2018）とのベンチマークにより、エンベロープ境界の金属量および積分された重元素収支が一致している場合、惑星の半径は勾配の正確な関数的トポロジーに対して頑健であることが示された。6つの異なるパラメータファミリー（ガウス型、シグモイド型、指数型、バイリニア型、シングルステップ型、およびデジタル化された流体力学的階段型）は、半径の広がりがわずか 2.3% であった。これは、滑らかな解析的パラメータ化（ガウス減衰など）を用いて縮退をマッピングすることの妥当性を検証している。
2. サブ・ネプチューンの縮退アトラス： 本フレームワークを $8 M_\oplus$ の惑星に適用し、水質量、エンベロープ質量、および希釈幅（$\sigma$）の $15 \times 15 \times 12$ グリッドにわたる「水の世界の縮退アトラス」を生成した。
   • $\mu$-コンパクション（平均分子量による緻密化）： 重元素をコンパクトなコアから希薄なエンベロープへと再分配することで、平均分子量が増加し、大気を物理的に緻密化させるプロセスを定量化している。この「$\mu$-コンパクション」は、固定された重元素収支に対して、熱膨張を抑制し、惑星半径を最大45%減少させることができる。
   • 熱ブランキング（熱的遮蔽）： 金属に富んだエンベロープは熱のスロットルとして機能し、表面とコアの間の温度差を約 5000 K 増大させる。ただし、これは静的な構造的知見であり、冷却履歴ではない。
   • 観測限界： 固定された $(M_p, R_p)$ に対して、希釈幅 $\sigma$ は実質的に制約されない。アトラスは、同じ半径を再現する一連の連続的な静的アーキテクチャ（水質量、エンベロープ質量、および $\sigma$ の変化）を明らかにしており、$(M_p, R_p)$ だけでは内部構造を解決できないことを示している。
3. 特定のターゲットへの適用：
   • GJ 1214 b（揮発性に富む）： コンパクトな水・岩石体と、薄く高度に濃縮されたエンベロープとしてモデル化され、観測値の 1% 以内の半径に収束した。
   • Kepler-11 e（ガスに富む／スーパー・パフ）： 大規模なコアと、高度に汚染された希薄なエンベロープとしてモデル化された。本モデルは、原始的な H/He エンベロープから予想される膨張に対抗し、観測されたコンパクトな半径まで大規模なエンベロープを「バラスト（重石）」として押し下げるために、深く懸濁された重元素勾配が必要であることを示している。
4. 大気と内部の結合： 本フレームワークは、大気リトリーバル（例：HADES との結合）のバックエンドとして機能するように設計されている。JWST/Ariel の大気金属量プライアを、自己整合的な表面重力と境界条件を解くことで、深部の内部構造へとリンクさせることが可能である。

意義と範囲
本論文は、fuzzycore を、完全な暗黙的進化型コードの代替としてではなく、必要な補完的ツールとして位置づけている。時間依存の混合やエネルギー輸送の計算オーバーヘッドを排除することで、fuzzycore は、ヘニイ型ソルバーでは探索が極めて困難であった、静的な構造の縮退を高解像度で体系的にマッピングすることを可能にする。

著者らは、fuzzycore が力学的に許容される静水圧構成を特定するものであることを強調している。本フレームワークは、それらの構成が、分離の物理学と熱力学的に自己整合しているか、あるいは地質学的時間スケールにおける乱流侵食に対して流体力学的に安定しているかを決定するものではない。それらの問いは、依然として進化モデリングおよび ab initio 相平衡研究の領域である。しかし、本フレームワークは、現在のおよび将来の観測データを解釈するために不可欠な「縮退アトラス」を提供しており、重元素の空間分布が、惑星特性評価における極めて重要かつ観測的に縮退した次元であることを示している。