[論文レビュー] Genesis of Polyatomic Molecules in Dark Clouds: CO$_2$ Formation on Cold Amorphous Solid Water

Meenu Upadhyay, Marco Pezzella+1 more|arXiv (Cornell University)|Apr 23, 2021

Astrophysics and Star Formation Studies被引用数 1

ひとこと要約

本研究では、高水準のポテンシャルエネルギー表面（PES）を用いた反応性分子動力学シミュレーションを用い、宇宙空間条件下におけるアモルファス・ソリッド・ワーエル（ASW）上でのO(1D) + CO(1Σ+)の再結合によるCO2生成を調査した。反応はサブナノ秒スケールでエネルギー障壁なしに進行し、CO2は内部エネルギーの迅速な再分配によってASW表面に安定化される。一方、高精度のRKHS PESは、COO中間体の形成や原子交換・脱吸着の詳細なダイナミクスを明らかにした。

ABSTRACT

Understanding the formation of molecules under conditions relevant to interstellar chemistry is fundamental to characterize the chemical evolution of the universe. Using reactive molecular dynamics simulations with model-based or high-quality potential energy surfaces provides a means to specifically and quantitatively probe individual reaction channels at a molecular level. The formation of CO<sub>2</sub> from collision of CO­(<sup>1</sup>Σ) and O­(<sup>1</sup>D) is characterized on amorphous solid water (ASW) under conditions typical in cold molecular clouds. Recombination takes place on the subnanosecond time scale and internal energy redistribution leads to stabilization of the product with CO<sub>2</sub> remaining adsorbed on the ASW on extended time scales. Using a high-level, reproducing kernel-based potential energy surface for CO<sub>2</sub>, formation into and stabilization of CO<sub>2</sub> and COO are observed.

研究の動機と目的

• 寒冷な宇宙空間条件下におけるアモルファス・ソリッド・ワーエル（ASW）上でのO(1D) + CO(1Σ+)の再結合によるCO2生成の分子レベルのメカニズムを調査すること。
• ポテンシャルエネルギー表面（PES）が、CO2の生成やCOOのような可能性のある中間体を正確に記述する上で果たす役割を評価すること。
• CO2生成、安定化、および原子交換や脱吸着といった競合経路の効率を定量化すること。
• 表面構造および温度がASW上でのO(1D)およびCOの反応性およびダイナミクスに与える影響を評価すること。
• 暗黒分子雲の化学的進化に関連する定量的で原子レベルのCO2生成の記述を提供すること。

提案手法

• 計算効率の高いモース-モース-ハーモニック（MMH）PESと、CCSD(T)-F12のab initioデータの再現型ヒルバート空間（RKHS）表現による高精度PESを用いた反応性分子動力学（RMD）シミュレーションを実施した。
• MMH-PESを用いて定性的なダイナミクスと初期スクリーニングを行い、RKHS-PESを用いてCO2およびCOO生成を含む反応経路の高精度な分析を実施した。
• 寒冷分子雲に相当する温度（約15 K）におけるASW上でのO(1D)およびCOの軌道をシミュレートし、原子間距離、角度、エネルギー再分配を追跡した。
• 再結合中にエネルギー保存則と系の温度変化をモニタリングし、シミュレーションの正確性を検証した。
• 初期O-CO間隔（R）およびOCO角（θ）を変化させた1000本以上の軌道を分析し、CO2生成、原子交換、脱吸着、複合体形成なしの結果に分類した。
• 径方向分布関数および時系列データを用いて、再結合および安定化過程における構造的・動的進化を分析した。

実験結果

リサーチクエスチョン

• RQ1 宇宙空間条件下におけるアモルファス・ソリッド・ワーエル上でのO(1D) + CO(1Σ+)再結合によるCO2生成のダイナミクスおよび時間スケールは何か？
• RQ2 異なるポテンシャルエネルギー表面（MMH対RKHS）が予測する反応経路および結果にどのように影響を与えるか？
• RQ3 CO2生成物がASW表面に安定化される過程で内部エネルギー再分配が果たす役割は何か？
• RQ4 COO中間体の形成や原子交換といった代替反応経路は観察されるか？その頻度はどの程度か？
• RQ5 初期条件（Rおよびθ）がCO2生成確率と比較して脱吸着や原子交換の確率にどの程度影響を与えるか？

主な発見

• CO2生成は、O(1D)とCO(1Σ+)の間でエネルギー障壁なしの再結合経路により、サブナノ秒スケール（約8 ps以内）で進行する。
• 再結合後、CO2はASW表面に長時間吸着し続け、内部エネルギー再分配が安定化をもたらす。
• 高精度のRKHSポテンシャルエネルギー表面は、CO2に加えCOO中間体の形成も明らかにした。シミュレーションではCOOの再吸着も観察された。
• 原子交換は初期条件に応じて1～17％の軌道で発生し、0～100％の系でCO、O、または両方の脱吸着が観察された。
• 初期OCO角が180°でR = 3.66 Åの条件下では、53％の軌道でCO2が生成され、292％の軌道で衝突複合体が形成されず、135％の軌道で脱吸着が観察された（複数の結果カテゴリが重複するため合計が100％を超える）。
• RKHS-PESではエネルギー保存が0.0026 eVのRMS誤差（R² = 1.0）の範囲内に維持され、シミュレーションの高精度が確認された。