A New Global Chemical Equilibrium Code: Refractory Element Signatures in Super-Earths and Sub-Neptunes

技術的概要：新しい全球化学平衡コード：スーパーアースとサブネプチューンにおける難揮発性元素のシグネチャ

問題提起
スーパーアースおよびサブネプチューンの大気は、長寿命のマグマ海洋フェーズ中に溶融した内部との化学的交換によって著しく変化する。ジェームズ・ウェッブ宇宙望遠鏡（JWST）などの施設からの大気観測を解釈するには、大気化学を惑星内部の組成と自己整合的に結合させるモデルが必要である。Schlichting & Young (2022) によって導入された全球化学平衡（GCE）枠組みは、金属 - 珪酸塩 - ガス系をモデル化する物理的に整合的な手法を提供したが、その元の実装はマルコフ連鎖モンテカルロ（MCMC）手法に依存していた。このアプローチは計算コストが高く、系外惑星の多様な集団を探索するために必要な大規模なパラメータ研究の実施を制限していた。

手法
著者らは、前バージョンの計算上のボトルネックに対処する GCE 枠組みの更新された実装を提示する。

• 数値ソルバ： 本研究では、MCMC ソルバを勾配ベースのオプティマイザ（ADAM）に置き換えた。数値的安定性を確保するため、ソルバはモル分率（$x_i$）に対して対数空間で動作する。モル分率は多数の桁にわたる可能性がある。
• 最適化ワークフロー： このコードは、Python の記号数学パッケージ SymPy を利用して支配方程式の微分を自動的に計算する。これらの微分は共通部分式除去（CSE）を用いて簡略化され、実行速度を最大化するために自動的に C++ コードに変換される。
• 熱力学的枠組み： このモデルは、3 相系（ガス、珪酸塩溶融物、金属）の平衡条件と質量保存の制約を表す連立非線形方程式系を解く。この系には 26 種の化学種と 19 の独立反応が含まれており、金属相への炭素の分配を含むように以前の定式化を拡張している。
• 水素の溶解度： 本研究は、珪酸塩溶融物中の水素溶解度に関する 2 つの代替処方箋を検討する。固定分配係数（$k_D$）と固定平衡定数（$k_{eq}$）である。$k_{eq}$処方箋は、圧力に伴う無制限の水素溶解を許容するものであり、最近の熱力学的制約が示唆する高い溶解度レベルを反映するため、主要な結果にはこれが採用された。
• 合成集団： この枠組みは、質量が 1 から 10 $M_\oplus$ の惑星の合成集団に適用される。本研究では、恒星のばらつきに基づいた範囲内で、取り込まれた水素質量分率（$w_{H_2}$）、水質量分率（$w_{H_2O}$）、およびバルクの難揮発性元素比（Mg/Si および Fe/Si）を変化させた。シミュレーションは、大気 - マグマ海洋界面（AMOI）温度 2000 K および 3000 K で実行された。

主な貢献

• 計算効率： 新しいソルバは、元の MCMC 実装と比較して化学ネットワークの解くための計算コストを 2 桁以上削減する。この効率性により、以前は実行不可能だった大規模なパラメータ研究が可能になった。
• オープンソースコード： 更新された GCE コードはオープンソースソフトウェアとしてリリースされ、モジュール設計となっており、CPU 実行向けに最適化されている。これにより、コアソルバを変更することなく、代替の反応セットや化学種の組み込みが可能である。
• 体系的な分析： 本論文は、バルクの難揮発性元素組成（Mg/Si、Fe/Si）および熱的状態が大気特性にどのように影響するか、特に顕著な水保有量を持つ惑星と持たない惑星について、包括的な分析を提供する。

結果

• 大気質量分率： 大気質量分率は主に AMOI 温度と惑星の水予算によって制御される。高い AMOI 温度は、揮発性の内部への溶解を促進し、大気質量分率を低下させる。$T_{AMOI} = 3000$ K の惑星では、取り込まれた水素質量に関わらず、水素の大部分が珪酸塩相に溶解するため、大気質量分率は約 1% 未満に低下する。取り込まれた水素質量分率は、水素の大部分が内部に分配されるため、最終的な大気質量を決定する上でわずかな役割しか果たさない。
• 大気金属質量分率： この指標は惑星の熱的状態（より高い温度はより高い金属量をもたらす）および水予算と強く相関する。水を蓄積した惑星は、乾燥した惑星に比べて 4 倍から 6 倍高い大気金属質量分率を示す。
• 難揮発性元素の影響： バルクの Mg/Si および Fe/Si 比が大気組成に及ぼす影響は、水の存在に強く依存する。
  • 乾燥惑星（$w_{H_2O} = 0$）： 難揮発性元素比は大気組成に弱い影響しか及ぼさない。炭素は Mg/Si または Fe/Si 比に関わらず、主に金属相に分配される。
  • 水に富む惑星（$w_{H_2O} > 0$）： 難揮発性元素比は炭素の分配を強く規定する。Mg/Si および Fe/Si 比が 1 未満の場合、炭素はガス相（CH$_4$、CO、CO$_2$として）に留まり、高い大気金属質量分率をもたらす。これらの比が 1 を超えると、炭素は金属相に分配され、気体炭素種が最大 5 桁にわたって劇的に減少し、それに伴って大気金属質量分率および C/O 比が低下する。
• C/O 比： 水に富む惑星の場合、大気 C/O 比は難揮発性元素比および水素溶解度に敏感である。採用された $k_{eq}$ 処方箋の下では、取り込まれた水素質量を増加させても、水素が内部に溶解するため、大気の還元力は著しく増加しない。その結果、水に富む惑星の C/O 比は一般に乾燥惑星よりも高いが、具体的な値は Mg/Si および Fe/Si 比に依存する。

重要性
本論文は、サブネプチューンの大気組成が単に揮発性物質の保有量によって決定されるのではなく、その岩石物質（難揮発性元素）のバルク組成および熱的進化に敏感に依存することを示している。著者らは、大気観測を解釈するには、大気と下部のマグマ海洋との間の地球化学的相互作用を考慮する必要があると主張する。

具体的には、本研究は将来の観測調査が、広範な Mg/Si および Fe/Si 比を持つ恒星を周回するサブネプチューンを対象とすべきであることを示唆している。宿主星の難揮発性元素比と惑星の大気 C/O 比との相関は、水氷線に対する形成場所の頑健な指標となり得る。本論文は、難揮発性組成が内部と大気間の炭素や酸素などの揮発性物質の分配に直接影響を与えるため、恒星の存在量制約がこれらの惑星の大気シグネチャを解釈する上で不可欠であると結論づけている。