[논문 리뷰] FastChem 2 : an improved computer program to determine the gas-phase chemical equilibrium composition for arbitrary element distributions

Joachim W. Stock, Daniel Kitzmann외 1명|arXiv (Cornell University)|2022. 06. 16.

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의 온도 범위에서 복잡한 열화학적 평형 계산의 계산 효율성을 향상시키기 위해.
• HNC, FeH, TiH와 같은 외계행성 및 별 대기에서 중요한 분자들을 포함해 총 28개 원소와 523종의 종을 추가하여 열화학 데이터베이스를 강화하기 위해.
• FastChem 1과의 후행 호환성을 유지하면서도 사용자 제공 데이터로 종을 필터링하거나 확장할 수 있도록 성능 튜닝 기능을 제공하기 위해.
• 수소가 없는 대기나 헬륨 농도가 낮은 상황을 포함한 다양한 대기 조성에서 신뢰성 있고 정확한 성능을 입증하기 위해.

제안 방법

• 원소 보존 조건 하에서 자유 에너지 최소화 문제를 재구성하여 임의의 원소 농도 분포를 처리할 수 있도록 반분석 알고리즘을 재설계하였다.
• 각 종에 대해 연관된 비선형 방정식 시스템으로 분해되는 변수 기반 반복적 해법을 구현하였으며, 수정된 뉴턴-라프슨 방법을 사용하였다.
• 종의 기준 상태를 위한 새로운 파arameterization을 도입하여 활성 계수를 조정함으로써 비수소 기반 시스템에 대한 민감도를 향상시켰다.
• 동적 종 필터링 기능을 도입하여 사용자가 성능 향상을 위해 종의 수를 줄이거나, 열화학적 데이터가 확보된 경우 새로운 종을 추가할 수 있도록 하였다.
• C++ 코어의 성능을 그대로 유지하면서도 과학 워크플로우에 원활하게 통합할 수 있도록 Python 인터페이스(pyfastchem)를 개발하였다.
• 523종의 종을 포함한 28개 원소로 구성된 열화학 데이터베이스를 업데이트하였으며, 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 1에서는 계산 시간이 3.2~5.2배 감소하지만, FastChem 2에서는 32~67% 감소하여 알고리즘 효율성이 향상됨을 시사한다.
• 저수소, 저헬륨 조건에서는 FastChem 2가 동일한 계산에 대해 FastChem 1 대비 3% 미만의 런타임을 소요하여 극적인 효율성 향상을 입증한다.
• FastChem 1과의 비교를 통해 검증되었으며, 질소, 탄소, 산소 기반 대기 모델링에 성공적으로 적용되어 다양한 조성에서 정확성과 견고성을 확인하였다.