[論文レビュー] FastChem 2 : an improved computer program to determine the gas-phase chemical equilibrium composition for arbitrary element distributions

Joachim W. Stock, Daniel Kitzmann+1 more|arXiv (Cornell University)|Jun 16, 2022

Atmospheric Ozone and Climate被引用数 2

ひとこと要約

FastChem 2 は、水素に支配される系に限定されない、任意の元素組成に対応できる、オープンソースの C++ および Python 互換コードであり、気相化学平衡組成を高速に計算する。FastChem 1 よりも最大 50 倍高速であり、HNC や FeH、TiH といった新しい分子を含む 523 種の物質をサポートしており、系外惑星や冷却星の大気における窒素、炭素、酸素に支配される大気の正確なモデル化を可能にする。

ABSTRACT

The computation of complex neutral/ionised chemical equilibrium compositions is invaluable to obtain scientific insights of, for example, the atmospheres of extrasolar planets and cool stars. We present FastChem 2, a new version of the established semi-analytical thermochemical equilibrium code FastChem. Whereas the original version is limited to atmospheres containing a significant amount of hydrogen, FastChem 2 is also applicable to chemical mixtures dominated by any other species, such as CO$_2$ or N$_2$. The new C++ code and an optional Python module are publicly available under the GPLv3 license at https://github.com/exoclime/FastChem. The program is backward compatible so that the previous version can be easily substituted. We updated the thermochemical database by adding HNC, FeH, TiH, Ca$^-$, and some organic molecules. In total 523 species are now in the thermochemical database including 28 chemical elements. The user can reduce the total number of species to, for example, increase the computation performance or can add further species if the thermochemical data are available. The program is validated against its previous version and extensively tested over an extended pressure-temperature grid with pressures ranging from $10^{-13}\,\mathrm{bar}$ up to $10^3\,\mathrm{bar}$ and temperatures between $100\,\mathrm{K}$ and $6000\,\mathrm{K}$. FastChem 2 is successfully applied to a number of different scenarios including nitrogen, carbon, and oxygen-dominated atmospheres as well as test cases without hydrogen and helium. Averaged over the extended pressure-temperature grid FastChem 2 is up to 50 times faster than the previous version and is also applicable to situations not treatable with version 1.

研究の動機と目的

• FastChem コードの適用範囲を水素に支配される系に限定されないよう拡張し、CO2 や N2 に支配される大気を含める。
• 広範な圧力範囲（10−13 から 103 bar）および温度範囲（100–6000 K）における複雑な熱力学的平衡計算の計算効率を向上させる。
• 28 種の元素と 523 種の物質を追加することで、系外惑星および恒星の大気において重要な HNC、FeH、TiH などの主要分子を含む熱力学的データベースを強化する。
• FastChem 1 との後方互換性を保ちつつ、物質のフィルタリングやユーザーが提供するデータによる拡張により、パフォーマンスチューニングを可能にする。
• 水素を含まない、ヘリウムが乏しい状況を含む多様な大気組成に対して、新しいコードの妥当性と正確性を検証する。

提案手法

• 元素保存則を満たすガウス自由エネルギー最小化問題を再定式化することで、任意の元素組成に対応できるように、半解析的アルゴリズムを再設計した。
• 各物質ごとに結合された非線形方程式系に分解する変数ベースの反復解法を実装し、修正ニュートン・ラプソン法を用いた。
• 物質の基準状態の新しいパrameterization を導入し、活性係数の調整により、水素に支配されない系に対しても柔軟に取り扱えるようにした。
• 動的物質フィルタリングを実装：パフォーマンス向上のため物質数を減らすか、熱力学的データが入手可能な場合に新しい物質を追加可能。
• Python インターフェース（pyfastchem）を開発し、科学的ワークフローへのシームレスな統合を可能にした。C++ コアのパフォーマンスを完全に活用。
• 28 要素にわたる 523 種の物質を含む熱力学的データベースを更新。FeH、TiH、Ca−、HNC などを含み、検証済みの熱力学的データを提供。

実験結果

リサーチクエスチョン

• RQ1 FastChem 2 は、CO2 や N2 や O2 に支配されるような水素由来でない物質に支配される大気に対し、化学平衡組成を正確に計算できるか？
• RQ2 FastChem 1 と比較して、FastChem 2 のパフォーマンスは、広範な圧力・温度範囲（10−13 bar から 103 bar、100 K から 6000 K）でどのようにスケーリングするか？
• RQ3 イオン種を除外することで、異なる元素組成の範囲で計算時間と正確性にどのような影響を与えるか？
• RQ4 特に低水素またはヘリウムが乏しい環境において、FastChem 2 は FastChem 1 よりもどれほど速く動作するか？
• RQ5 FastChem 2 は、水素やヘリウムが極めて少ない状況を含め、任意の元素組成に耐えうるか、信頼性高く大気をモデル化できるか？

主な発見

• FastChem 2 は、全圧力・温度グリッドにおいて、FastChem 1 より最大 50 倍高速であり、低水素環境においては性能向上がさらに顕著に現れる。
• 太陽の光球層における水素の豊度では、FastChem 2 は FastChem 1 が要する時間の 10–20% のみで実行される。
• 水素の豊度が太陽のヘリウムの豊度未満の場合、FastChem 1 とは逆に、FastChem 2 では反復収束ダイナミクスの変化により、計算時間の挙動が逆転する。
• イオン化学を除外することで、FastChem 1 では計算時間が 3.2～5.2 倍短縮されるが、FastChem 2 では 32～67% の短縮が達成され、アルゴリズムの効率性が向上していることが示された。
• 低水素・低ヘリウム環境では、FastChem 2 は同じ計算に対して FastChem 1 の 3% 未満の実行時間で完了し、顕著な効率向上が確認された。
• FastChem 1 との比較と、窒素、炭素、酸素に支配される大気への応用を通じて、コードの正確性と耐障害性が検証され、多様な組成にわたって妥当であることが確認された。