[論文レビュー] Mercury-Ar$χ$es: a high-performance n-body code for planet formation studies

D. Turrini, S. Fonte+5 more|arXiv (Cornell University)|Jan 23, 2026

Astro and Planetary Science被引用数 0

ひとこと要約

Mercury-Arχes は Mercury+ をベースにした高性能の OpenMP 並列 n-体コードで、惑星の形成と成長する惑星、ディスクガス相互作用、微惑星体ダイナミクスをシミュレートします。

ABSTRACT

Forming planetary systems are populated by large numbers of gravitationally interacting planetary bodies, spanning from massive giant planets to small planetesimals akin to present-day asteroids and comets. All these planetary bodies are embedded in the gaseous embrace of their native protoplanetary disks, and their interactions with the disk gas play a central role in shaping their dynamical evolution and the outcomes of planet formation. These factors make realistic planet formation simulations extremely computationally demanding, which in turn means that accurately modeling the formation of planetary systems requires the use of high-performance methods. The planet formation code Mercury-Ar$χ$es was developed to address these challenges and, since its first implementation, has been used in multiple exoplanetary and Solar System studies. Mercury-Ar$χ$es is a parallel n-body code that builds on the widely used Mercury code and is capable of modeling the growth and migration of forming planets, the interactions between planetary bodies and the disk gas, as well as the evolving impact flux of planetesimals on forming planets across the different stages of their formation process. In this work we provide the up-to-date overview of its physical modeling capabilities and the first detailed description of its high-performance implementation based on the OpenMP directive-based parallelism for shared memory environments, to harness the multi-thread and vectorization features of modern processor architectures.

研究の動機と目的

• Mercury-Arχes のモデリング機能と高性能な OpenMP 実装の自己完結描述を提供する。
• Mercury+ が n-体エンジンをどのように形成し、共有メモリアーキテクチャ向けの並列化戦略をどのように設計しているかを説明する。
• Arχes 多物理ライブラリのディスクガス相互作用、惑星成長、移動、微惑星 disk モデル化を詳述する。
• 現実的なワークロードで異なるハードウェア環境を用いた Mercury-Arχes の性能とスケーリングを示す。

提案手法

• 高性能ハイブリッドシンプリティック n-体エンジンとしての Mercury+ を、OpenMP ベクトル化実装と Kepler 損失 drift の WHFAST インタフェースとともに提示する。
• Arχes ライブラリの構成要素（arxes_common.f90、arxes_disk.f90、arxes_growth.f90、arxes_migration.f90、arxes.for）と、それらがディスクガス相互作用、惑星成長、移動に果たす役割を説明する。
• 表面密度と温度プロファイル、ガス抗力（Re および Ma に依存する C_D）、 disk 重力、微惑星体へのガス減衰効果を含むガス-ディスクモデリングを説明する。
• 物質沈着後にガス放出までの二相惑星成長モデルを、質量と半径の進化の式（式13-16）とともに詳細に説明する。
• 二相移動モデル（線形およびべき乗則域）とドリフト速度（式17-18）を成長フェーズと結合させて説明する。
• 雪線の位置に基づく慣性質量を持つダイナミクス標識子としての微惑星ディスクの扱いを述べる。

実験結果

リサーチクエスチョン

• RQ1 形成惑星と出身ディスクの結合進化（成長、移動、ディスクガス相互作用を含む）を Mercury-Arχes フレームワークはどのようにモデル化するか。
• RQ2 高解像度の惑星形成シミュレーションにおける現代アーキテクチャで、OpenMP 並列化とベクタ化によって得られる性能向上とスケーラビリティはどの程度か。
• RQ3 ディスク重力、ガス抗力、惑星摂動がこのフレームワークにおける動的進化と形成結果にどのように影響するか。

主な発見

• Mercury+ は高性能なハイブリッドシンプリティックエンジンを提供し、OpenMP 実装と正確・効率的な Keplerian drift の WHFAST インタフェースを備える。
• Mercury-Arχes はディスクガス相互作用、惑星成長、移動を含む包括的な多物理モデリングをサポートし、初期惑星形成段階の自己一致した研究を可能にする。
• ガス-ディスク相互作用は特定のアーキテクチャで計算負荷を支配しうるため、ガス関連カーネルの並列化およびベクトル化を正当化する。
• コードはシリアルと並列実行の再現性を保つための設計上の工夫（ヘッダ規約とガード付き並列領域など）を採用している。
• 微惑星ディスクは質量を持たないトレーサとして慣性質量を有するモデル化を採用し、相互重力なしで減衰と輸送計算を可能にする。