[論文レビュー] Chemistry in Disks. IV. Benchmarking gas-grain chemical models with surface reactions

D. Semenov, F. Hersant+8 more|SPIRE - Sciences Po Institutional REpository|Jul 14, 2010

Astrophysics and Star Formation Studies参考文献 48被引用数 118

ひとこと要約

本論文は、OSU_03_2008 率ファイルに基づく統一された反応ネットワークを用いて、寒冷分子雲、ホットコア、およびプロト惑星系円盤領域という多様な天体物理学的環境において、最先端のガス・グライン化学モデル ALCHEMIC と NAUTILUS をベンチマークした。コード実装の差異にもかかわらず、物理的パラメータ、グライン特性、表面反応処理をきめ細かく一致させた結果、全物質の時間依存的混合比について完全な一致が達成され、今後の宇宙化学的シミュレーションおよび ALMA 観測データ解釈のための公開可能な基準フレームワークを確立した。

ABSTRACT

Abridged: We detail and benchmark two sophisticated chemical models developed by the Heidelberg and Bordeaux astrochemistry groups. The main goal of this study is to elaborate on a few well-described tests for state-of-the-art astrochemical codes covering a range of physical conditions and chemical processes, in particular those aimed at constraining current and future interferometric observations of protoplanetary disks. We consider three physical models: a cold molecular cloud core, a hot core, and an outer region of a T Tauri disk. Our chemical network (for both models) is based on the original gas-phase osu_03_2008 ratefile and includes gas-grain interactions and a set of surface reactions for the H-, O-, C-, S-, and N-bearing molecules. The benchmarking is performed with the increasing complexity of the considered processes: (1) the pure gas-phase chemistry, (2) the gas-phase chemistry with accretion and desorption, and (3) the full gas-grain model with surface reactions. Using atomic initial abundances with heavily depleted metals and hydrogen in its molecular form, the chemical evolution is modeled within 10^9 years. The time-dependent abundances calculated with the two chemical models are essentially the same for all considered physical cases and for all species, including the most complex polyatomic ions and organic molecules. This result however required a lot of efforts to make all necessary details consistent through the model runs, e.g. definition of the gas particle density, density of grain surface sites, the strength and shape of the UV radiation field, etc. The reference models and the benchmark setup, along with the two chemical codes and resulting time-dependent abundances are made publicly available in the Internet: http://www.mpia.de/homes/semenov/Chemistry_benchmark/home.html

研究の動機と目的

• プロト惑星系円盤および恒星間媒体の研究に用いられる高度なガス・グライン化学モデルの標準化・再現可能性のあるベンチマークを確立すること。
• 物理的パラメータ、反応ネットワーク、数値処理の一致を図ることで、代表的な2つの宇宙化学的コード（ALCHEMIC と NAUTILUS）の乖離を解消すること。
• 今後のシミュレーション開発および検証のための、公開可能で詳細なドキュメント化された基準モデルセットを提供すること。
• ALMA や今後の高分解能干渉観測に対して、分子混合比を予測するモデルの信頼性と一貫性を保証すること。
• 初期条件、反応速度、グライン表面物理学の不一致が引き起こす長年のモデル比較の困難を解決すること。

提案手法

• OSU_03_2008 ガス相率ファイルに基づく統一化学ネットワークの実装を、H、O、C、N、S、およびS含有種のガス・グライン相互作用および表面反応を含む形に拡張。
• 5つの代表的物理モデルの使用：寒冷分子雲コア（TMC1）、ホットコア（Hot Corino）、および温度、密度、UV放射場が異なる約100 AUの円盤領域3か所。
• ガス粒子密度の定義、UV放射場スケーリング、ダスト粒子表面サイト密度、および接着係数といったモデルパラメータの段階的一致。
• 表面種の脱着および拡散エネルギー値の標準化、および電子接着とグライン帯電の一貫した取り扱い。
• 同一単位および換算係数（例：N_H と A_V の間）の採用により、コード間の数値的一致性を確保。
• 表面反応に同一の率式を採用し、均一系反応と多相反応の率の違い（例：H + H → H₂）を慎重に補正。

実験結果

リサーチクエスチョン

• RQ1 同一の物理的・化学的条件下で、独立した2つの宇宙化学的コード（ALCHEMIC と NAUTILUS）が、時間依存的分子混合比をどれほど同一に生成するか？
• RQ2 高度なガス・グライン化学モデル間で完全な一致を達成するために、どのような物理的・化学的パラメータを標準化する必要があるか？
• RQ3 表面反応、脱着機構、グライン特性が、プロト惑星系円盤環境における化学的進化にどのように影響を与えるか？
• RQ4 再現可能性と比較可能性を向上させるために、一貫性があり公開可能なベンチマークモデルを確立できるか？
• RQ5 ALMA や今後の高分解能干渉観測の信頼できる予測を保証するため、どのような実装詳細を調和させる必要があるか？

主な発見

• ALCHEMIC と NAUTILUS の2つの化学モデルは、全5つの物理的モデルにおいて、複雑な多原子イオンおよび有機分子を含む全物質の時間依存的混合比を同一に生成した。
• ガス粒子密度の定義、UV放射場スケーリング、ダスト粒子表面サイト密度といった物理的パラメータの広範な一致がなされた後、のみ完全な一致が達成された。
• 初期モデル走行での乖離は、原子質量定義（同位体 vs. 主要同位体）および接着係数処理の違いに起因することが特定された。
• 表面種の脱着および拡散エネルギー値が、モデル収束を達成するための重要なパラメータであることが同定された。
• ベンチマークプロセスにより、両コードにわずかな誤りが発見され、是正されたことで、それらの予測信頼性が向上した。
• 反応ネットワーク、物理的パラメータ、出力形式を含む最終的な調和化モデルセットは、今後のモデル開発および観測解析を支援する目的で公開された。