FastChem: A computer program for efficient complex chemical equilibrium calculations in the neutral/ionized gas phase with applications to stellar and planetary atmospheres
Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 kernel pith:IZUSV4ZIrecord.jsonopen to challenge →
read the original abstract
For the calculation of complex neutral/ionized gas phase chemical equilibria, we present a semi-analytical versatile and efficient computer program, called FastChem. The applied method is based on the solution of a system of coupled nonlinear (and linear) algebraic equations, namely the law of mass action and the element conservation equations including charge balance, in many variables. Specifically, the system of equations is decomposed into a set of coupled nonlinear equations in one variable each, which are solved analytically whenever feasible to reduce computation time. Notably, the electron density is determined by using the method of Nelder and Mead at low temperatures. The program is written in object-oriented C++ which makes it easy to couple the code with other programs, although a stand-alone version is provided. FastChem can be used in parallel or sequentially and is available under the GNU General Public License version 3 at https://github.com/exoclime/FastChem together with several sample applications. The code has been successfully validated against previous studies and its convergence behavior has been tested even for extreme physical parameter ranges down to 100 K and up to 1000 bar. FastChem converges stable and robust in even most demanding chemical situations, which posed sometimes extreme challenges for previous algorithms.
This paper has not been read by Pith yet.
Forward citations
Cited by 7 Pith papers
-
The 35-Myr old infant planet TOI-837 b has a mildly misaligned orbit
TOI-837 b has a true obliquity of 25.9+7.5-6.3 deg, the first planet younger than 100 Myr with accessible ψ incompatible with an aligned orbit, favoring primordial disc torque followed by disc-driven migration.
-
Native-resolution retrievals of VHS 1256-1257 b spanning the JWST/NIRSpec wavelength range: Chemical composition of a partially cloudy atmosphere
Native-resolution retrievals on 0.97-5.27 micron JWST spectra of VHS 1256 b find a ~79% cloud deck, solar-like metallicity and C/O, and depleted 18O in a partially cloudy disequilibrium atmosphere.
-
Unraveling the Mystery of the Peculiar and Young Hot Jupiter CoRoT-2b II: Phase Resolved Emission Spectroscopy with VLT/CRIRES+ and Gemini-S/IGRINS
Phase-resolved high-resolution spectroscopy of CoRoT-2b measures sub-synchronous rotation at 2.6-sigma significance, consistent with its western hotspot offset.
-
The CRIMSON survey I: super-stellar SiO in the directly imaged companion TWA 5 B from high-resolution M-band spectroscopy
High-resolution M-band spectroscopy detects super-stellar SiO in TWA 5 B, implying no significant magnesium-silicate clouds and formation consistent with core accretion beyond the CO snowline or gravitational instabil...
-
A parameterised approach to disequilibrium retrievals in the JWST era: Application to NIRCam observations of HD 189733b
A quench-pressure parameterisation with two free parameters recovers unbiased C/O and [M/H] from synthetic JWST spectra with vertical mixing, and tentatively detects quenching and photochemical H2S depletion in HD 189...
-
The GAPS programme at TNG ?. TOI-1533: a compact system hosting a super-Neptune-mass pair with disparate radii
TOI-1533 hosts an inner sub-Neptune (P=3.63 d, R=3.15 R⊕) and outer super-Neptune-mass hot giant (P=8.06 d, R>7.5 R⊕, M≈40 M⊕, ρ<0.48 g cm⁻³) both transiting an active K-dwarf.
-
The ESO SupJup Survey XI. Atmospheric properties of six isolated M- and L-type dwarfs with CRIRES+
Retrievals on six isolated brown dwarfs yield near-solar C/O (0.51-0.63), metallicities, and 12C/13C ratios (91-155) supporting molecular cloud fragmentation origin.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.