Papers › Toward More Reliable Analytic Thermochemical-equilibrium Abundances
Toward More Reliable Analytic Thermochemical-equilibrium Abundances
Patricio Cubillos, Jasmina Blecic, Ian Dobbs-Dixon
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
Heng & Tsai (2016) developed an analytic framework to obtain thermochemical-equilibrium abundances for H₂O, CO, CO₂, CH₄, C₂H₂, C₂H₄, HCN, NH₃, and N₂ for a system with known temperature, pressure, and elemental abundances (hydrogen, carbon, nitrogen, and oxygen). However, the implementation of their approach can become numerically unstable under certain circumstances, leading to inaccurate solutions (e.g., C/O ≥1 atmospheres at low pressures). Building up on their approach, we identified the conditions that prompt inaccurate solutions, and developed a new framework to avoid them, providing a reliable implementation for arbitrary values of temperature (200 to ∼2000 K), pressure (10⁻⁸ to 10³ bar), and CNO abundances (10⁻³ to ∼10²× solar elemental abundances), for hydrogen-dominated atmospheres. The accuracy our analytic framework is better than 10% for the more abundant species that have mixing fractions larger than 10⁻¹⁰, whereas the accuracy is better than 50% for the less abundant species. Additionally, we added the equilibrium-abundance calculation of atomic and molecular hydrogen into the system, and explored the physical limitations of this approach. Efficient and reliable tools, such as this one, are highly valuable for atmospheric Bayesian studies, which need to evaluate a large number of models. We implemented our analytic framework into the \textsc{rate} Python open-source package, available at https://github.com/pcubillos/rate .
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections