FastChem 4: New chemical elements and improved convergence behaviour
Pith reviewed 2026-05-20 00:15 UTC · model grok-4.3
The pith
FastChem 4 stabilizes chemical equilibrium calculations for non-solar abundances and identifies silicon monoxide as a condensate in a limited range.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the new solver methods and expanded chemical network enable robust calculations across a wide range of elemental compositions. This leads to the discovery that silicon monoxide forms as a stable condensate under carbon-rich conditions within a limited pressure-temperature window, consistent with recent observations.
What carries the argument
The multidimensional Newton-method in logarithmic element densities for the gas-phase solver, which improves convergence and eliminates the need for quad-precision arithmetic.
If this is right
- The code achieves a strong increase in computational performance and stability.
- It reproduces the classical solar-composition condensation sequence.
- Marked shifts in the condensation sequence appear under carbon-rich conditions.
- Silicon monoxide remains stable as a condensate over a limited pressure-temperature range.
Where Pith is reading between the lines
- The stability of silicon monoxide as a condensate could explain certain spectral features seen in brown dwarf atmospheres.
- Improved handling of non-solar abundances may help model the chemistry of exoplanets with unusual compositions inferred from observations.
- Future extensions could incorporate even more elements or couple the equilibrium solver to atmospheric dynamics models.
Load-bearing premise
The tabulated thermochemical data for the 44 elements and 1311 species are assumed to be accurate enough that any errors in the Gibbs energies would not change the reported condensation sequences or the stability range for silicon monoxide.
What would settle it
A mismatch between the predicted pressure-temperature range for silicon monoxide condensation and direct spectroscopic observations of carbon-rich brown dwarfs would falsify the stability result.
Figures
read the original abstract
Chemical equilibrium calculations are a key ingredient for modelling and interpreting spectroscopic observations of (exo)planets, brown dwarfs, cool stars, and protoplanetary disks. As these applications increasingly probe non-solar elemental abundances and previously underrepresented elements, equilibrium chemistry solvers must be both numerically robust and capable of handling complex chemical systems. Here we present FastChem 4, a major update to the open-source FastChem equilibrium chemistry code. We extend the gas-phase solver with a multidimensional Newton-method that mitigates the slow convergence previously encountered for strongly non-solar elemental abundances. We further reformulate the gas-phase equations in logarithmic element densities, removing the dependence on quad-precision arithmetic and allowing FastChem to be applied at low temperatures on any platform supporting double precision. The condensate solver is upgraded with adaptive Levenberg-Marquardt regularisation, a perturbed-Hessian fallback, and a combined gas-condensate Newton solver. These changes lead to a strong increase in computational performance and stability. The thermochemical data is expanded using thermochemical data from the NIST-JANAF tables and the Barin compilation, and now comprises 800 gas-phase molecules and ions and 511 condensates spanning 44 elements. We apply the updated code to a wide pressure-temperature grid for both solar and carbon-rich (C/O = 2) elemental compositions. The resulting grids reproduce the classical solar-composition condensation sequence and reveal the marked shifts that occur under carbon-rich conditions. We also find that silicon monoxide is stable as a condensate over a limited pressure-temperature range, consistent with recent JWST observations of brown dwarfs. FastChem 4 is released under the GPLv3 licence, together with a pre-compiled Python package.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents FastChem 4, a major update to the open-source FastChem chemical equilibrium code. It introduces a multidimensional Newton solver for the gas phase, reformulates equations in logarithmic element densities to eliminate quad-precision requirements, and upgrades the condensate solver with adaptive Levenberg-Marquardt regularization, perturbed-Hessian fallback, and a combined gas-condensate Newton method. The thermochemical database is expanded to 800 gas-phase species and 511 condensates spanning 44 elements using NIST-JANAF and Barin data. Application to pressure-temperature grids for solar and C/O=2 compositions reproduces the classical solar condensation sequence, shows marked shifts under carbon-rich conditions, and identifies a limited P-T domain where SiO is stable as a condensate, consistent with recent JWST observations of brown dwarfs. The code is released under GPLv3 with a pre-compiled Python package.
Significance. If the numerical improvements and reported sequences hold, this provides a more robust and broadly applicable tool for equilibrium chemistry modeling of exoplanet atmospheres, brown dwarfs, cool stars, and disks, especially for non-solar abundances and underrepresented elements. The open-source release with Python package and the explicit demonstration of SiO stability as a potential observational link are strengths that support wider use and reproducibility.
major comments (3)
- [Results, P-T grid application] In the section presenting the P-T grids and SiO stability (results on carbon-rich and solar sequences): the reported limited pressure-temperature range where SiO appears as a stable condensate depends directly on the Gibbs energies for the new SiO condensate entry from the Barin compilation. No sensitivity test to plausible variations in formation enthalpy or heat capacity is shown, which is load-bearing for the claim of consistency with JWST observations.
- [Thermochemical data expansion] In the methods description of the expanded thermochemical data: the addition of 511 new condensates and 800 gas species assumes NIST-JANAF and Barin tables are free of systematic errors across the full P-T range of interest. No cross-validation against FastChem 3 results for the overlapping species subset or against an independent equilibrium code is reported, preventing isolation of data error from the claimed solver improvements.
- [Solver improvements and performance] In the section on solver upgrades and performance: the abstract and text claim a 'strong increase in computational performance and stability,' yet no quantitative benchmarks (e.g., iteration counts, failure rates, or wall-time ratios on the same grids) comparing version 3 to version 4 are provided for the tested solar and C/O=2 cases.
minor comments (2)
- [Figures] Figure captions for the condensation sequences should explicitly state the exact elemental abundances and pressure range used to allow direct reproduction.
- [Abstract] The abstract states performance gains qualitatively; adding a sentence with specific metrics (e.g., factor of X speedup or reduction in non-convergent cases) would strengthen the summary.
Simulated Author's Rebuttal
We thank the referee for their thorough and constructive review of our manuscript. We address each of the major comments point by point below, indicating where revisions will be made to improve the paper.
read point-by-point responses
-
Referee: [Results, P-T grid application] In the section presenting the P-T grids and SiO stability (results on carbon-rich and solar sequences): the reported limited pressure-temperature range where SiO appears as a stable condensate depends directly on the Gibbs energies for the new SiO condensate entry from the Barin compilation. No sensitivity test to plausible variations in formation enthalpy or heat capacity is shown, which is load-bearing for the claim of consistency with JWST observations.
Authors: We agree that the stability range of SiO as a condensate is sensitive to the thermochemical data adopted from the Barin tables. To strengthen this aspect of the manuscript, we will add a sensitivity test in the revised version. Specifically, we will vary the standard enthalpy of formation and heat capacity data for SiO within estimated uncertainties (e.g., ±5-10 kJ/mol for enthalpy) and demonstrate how this affects the P-T domain of stability. This will provide a more robust basis for the claimed consistency with JWST observations of brown dwarfs. revision: yes
-
Referee: [Thermochemical data expansion] In the methods description of the expanded thermochemical data: the addition of 511 new condensates and 800 gas species assumes NIST-JANAF and Barin tables are free of systematic errors across the full P-T range of interest. No cross-validation against FastChem 3 results for the overlapping species subset or against an independent equilibrium code is reported, preventing isolation of data error from the claimed solver improvements.
Authors: We acknowledge the importance of validating the thermochemical database to distinguish between improvements due to the new solver and potential data differences. In the revised manuscript, we will include a cross-validation section comparing equilibrium compositions from FastChem 4 with those from FastChem 3 for the subset of species present in both versions, using identical P-T conditions and elemental abundances. We will also discuss any discrepancies and their origins. A comparison with an independent code such as TEA or GGchem could be added if space permits, but the primary validation will be against the previous FastChem version. revision: yes
-
Referee: [Solver improvements and performance] In the section on solver upgrades and performance: the abstract and text claim a 'strong increase in computational performance and stability,' yet no quantitative benchmarks (e.g., iteration counts, failure rates, or wall-time ratios on the same grids) comparing version 3 to version 4 are provided for the tested solar and C/O=2 cases.
Authors: We appreciate this observation regarding the lack of quantitative performance metrics. Although the improvements are described in the text, we will add quantitative benchmarks in the revised manuscript. This will include a table or figure presenting average iteration counts, convergence success rates, and relative wall-clock times for FastChem 3 versus FastChem 4 on the same solar and C/O=2 P-T grids. These metrics will directly support the claims of enhanced performance and stability. revision: yes
Circularity Check
No circularity in FastChem 4 code update and condensation results
full rationale
The paper describes numerical upgrades to the FastChem equilibrium solver (multidimensional Newton, logarithmic element densities, adaptive Levenberg-Marquardt) and the addition of external thermochemical data from NIST-JANAF and Barin tables for 800 gas species and 511 condensates. The reported solar and C-rich condensation sequences, including the limited P-T domain for SiO stability, are direct numerical outputs from applying the updated solver to these independent input tables. No step reduces a claimed result to a fitted parameter or self-citation by construction; the consistency note with JWST observations functions as an external cross-check rather than an input. The derivation chain is therefore self-contained against external benchmarks and data sources.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Newton's method converges to the physical solution for the chosen initial guesses and regularization schedule across the full P-T-abundance domain.
- domain assumption Gibbs free energies tabulated in NIST-JANAF and Barin compilations are accurate enough for all 44 elements and 1311 species at the pressures and temperatures considered.
Reference graph
Works this paper leans on
- [1]
-
[2]
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
work page 2009
-
[3]
Bachmann, N., Kreidberg, L., Mollière, P., Deming, D., & Tsai, S.-M. 2025, A&A, 700, A105
work page 2025
-
[4]
1995, Thermochemical Data of Pure Substances (VCH Weinheim)
Barin, I. 1995, Thermochemical Data of Pure Substances (VCH Weinheim)
work page 1995
-
[5]
Barklem, P. S. & Collet, R. 2016, A&A, 588, A96
work page 2016
-
[6]
Blecic, J., Harrington, J., & Bowman, M. O. 2016, ApJS, 225, 4
work page 2016
-
[7]
Borsato, N. W., Hoeijmakers, H. J., Prinoth, B., et al. 2023, A&A, 673, A158
work page 2023
-
[8]
Bower, D. J., Thompson, M. A., Hakim, K., Tian, M., & Sossi, P. A. 2025, ApJ, 995, 59
work page 2025
-
[9]
Burcat, A. & Ruscic, B. 2005, Third Millennium Ideal Gas and Condensed Phase Thermochemical Database for Combustion with updates from Active Ther- mochemical Tables (Argonne National Laboratory)
work page 2005
-
[10]
Burrows, A., Dulick, M., Bauschlicher, C. W., J., et al. 2005, ApJ, 624, 988
work page 2005
- [11]
-
[12]
1998, NIST-JANAF Thermochemical Tables (American Institute of Physics)
Chase, M. 1998, NIST-JANAF Thermochemical Tables (American Institute of Physics)
work page 1998
-
[13]
Dabbs, D. & Lynch, D. 1983, Metallurgical Transactions B, 14, 502 de Regt, S., Snellen, I. A. G., Allard, N. F., et al. 2025, A&A, 696, A225
work page 1983
-
[14]
Dorofeeva, O. V ., Iorish, V . S., Novikov, V . P., & Neumann, D. B. 2003, Journal of Physical and Chemical Reference Data, 32, 879
work page 2003
-
[15]
Dorofeeva, O. V ., Novikov, V . P., & Neumann, D. B. 2001, Journal of Physical and Chemical Reference Data, 30, 475
work page 2001
-
[16]
Vapor Pressure and Antoine Con- stants for Nitrogen Containing Organic Compounds
Dykyj, J., Svoboda, J., Wilhoit, R. C., Frenkel, M., & Hall, K. R. 2001, 2 Or- ganic Compounds, C0 to C84: Datasheet from Landolt-Börnstein - Group IV Physical Chemistry · V olume 20C: “Vapor Pressure and Antoine Con- stants for Nitrogen Containing Organic Compounds” in SpringerMaterials (https://doi.org/10.1007/10688591_3) (Springer-Verlag Berlin Heidelb...
-
[17]
Fisher, C. E., Hooton, M. J., Gressier, A., et al. 2026, MNRAS, 545, staf2187
work page 2026
-
[18]
Gail, H.-P. & Sedlmayr, E. 2013, Cambridge Astrophysics Series, V ol. 52, Physics and Chemistry of Circumstellar Dust Shells (Cambridge, UK: Cam- bridge University Press)
work page 2013
-
[19]
P., Wetzel, S., Pucci, A., & Tamanai, A
Gail, H. P., Wetzel, S., Pucci, A., & Tamanai, A. 2013, A&A, 555, A119
work page 2013
-
[20]
Gokcen, N. A. 1989, Bulletin of Alloy Phase Diagrams, 10, 11
work page 1989
-
[21]
Goodwin, R. D. 1985, Journal of Physical and Chemical Reference Data, 14, 849
work page 1985
-
[22]
Goos, E., Burcat, A., & Ruscic, B. 2023, Technion-IIT FTP Server
work page 2023
-
[23]
Gordon, S. & McBride, B. 1994, Computer program for calculation of complex chemical equilibrium compositions and applications, NASA RP-1311
work page 1994
-
[24]
Grasser, N., Snellen, I. A. G., de Regt, S., et al. 2025, A&A, 698, A252
work page 2025
-
[25]
1972, Geochimica et Cosmochimica Acta, 36, 597
Grossman, L. 1972, Geochimica et Cosmochimica Acta, 36, 597
work page 1972
-
[26]
Haar, L. & Gallagher, J. S. 1978, Journal of Physical and Chemical Reference Data, 7, 635
work page 1978
-
[27]
Hakim, K., Bower, D. J., Seidler, F. L., & Sossi, P. A. 2026, MNRAS, 546, stag133
work page 2026
-
[28]
Hill, E. A. 1900, Journal of the American Chemical Society, 22, 478
work page 1900
-
[29]
J., Ehrenreich, D., Kitzmann, D., et al
Hoeijmakers, H. J., Ehrenreich, D., Kitzmann, D., et al. 2019, A&A, 627, A165
work page 2019
-
[30]
1973, Selected values of the thermodynamic properties of the ele- ments
Hultgren, R. 1973, Selected values of the thermodynamic properties of the ele- ments. (Metals Park, Ohio: American Society for Metals)
work page 1973
-
[31]
Ito, Y ., Kimura, T., Ohno, K., Fujii, Y ., & Ikoma, M. 2025, ApJ, 987, 174
work page 2025
-
[32]
Kitzmann, D., Hoeijmakers, H. J., Grimm, S. L., et al. 2023, A&A, 669, A113
work page 2023
-
[33]
Kitzmann, D., Stock, J. W., & Patzer, A. B. C. 2024, MNRAS, 527, 7263
work page 2024
-
[34]
Larimer, J. W. & Bartholomay, M. 1979, Geochimica et Cosmochimica Acta, 43, 1455
work page 1979
-
[35]
Leal, A. M. M., Kulik, D. A., Kosakowski, G., & Saar, M. O. 2016, Advances in Water Resources, 96, 405
work page 2016
-
[36]
Lewis, J. S. 1969, Icarus, 10, 365
work page 1969
-
[37]
Lewis, J. S. 1972, Icarus, 16, 241
work page 1972
-
[38]
2009, CRC Handbook of Chemistry and Physics, 90th Edition (Taylor & Francis)
Lide, D. 2009, CRC Handbook of Chemistry and Physics, 90th Edition (Taylor & Francis)
work page 2009
-
[39]
1999, Journal of Physical and Chemical Reference Data, 28, 1705
Lodders, K. 1999, Journal of Physical and Chemical Reference Data, 28, 1705
work page 1999
- [40]
- [41]
-
[42]
2004, Journal of Physical and Chemical Reference Data, 33, 357
Lodders, K. 2004, Journal of Physical and Chemical Reference Data, 33, 357
work page 2004
-
[43]
Lodders, K., Bergemann, M., & Palme, H. 2025, Space Sci. Rev., 221, 23
work page 2025
-
[44]
Lodders, K. & Fegley, B. 1993, Earth and Planetary Science Letters, 117, 125
work page 1993
- [45]
-
[46]
Lodders, K. & Fegley, B., J. 1997, The Astrophysical Journal Letters, 484, L71
work page 1997
-
[47]
Lodders, K. & Fegley, B., J. 1999, in IAU Symposium, V ol. 191, Asymptotic Giant Branch Stars, ed. T. Le Bertre, A. Lèbre, & C. Waelkens (Astronomical Society of the Pacific), 279–290
work page 1999
-
[48]
McBride, B., Zehe, M., & Gordon, S. 2002, NASA Glenn coefficients for cal- culating thermodynamic properties of individual species, NASA TP-2002- 211556
work page 2002
-
[49]
McBride, B. J. & Gordon, S. 1996, Computer Program for Calculation of Com- plex Chemical Equilibrium Compositions and Applications II. Users Manual and Program Description., Tech. rep., Glenn Research Center - NASA Mollière, P., Kühnle, H., Matthews, E. C., et al. 2025, A&A, 703, A79
work page 1996
-
[50]
Moses, J. I., Allen, M., & Yung, Y . L. 1992, Icarus, 99, 318
work page 1992
-
[51]
Moses, J. I., Visscher, C., Fortney, J. J., et al. 2011, ApJ, 737, 15
work page 2011
-
[52]
Murphy, D. M. & Koop, T. 2005, Quarterly Journal of the Royal Meteorological Society, 131, 1539
work page 2005
- [53]
-
[54]
Pankratz, L. B. 1984, Thermodynamic data for mineral technology, Bulletin (United States. Bureau of Mines) ; 677 (Washington, D.C.]: U.S. Dept. of the Interior, Bureau of Mines : For sale by the Supt. of Docs., U.S. G.P.O.)
work page 1984
- [55]
-
[56]
G., de Regt, S., Gandhi, S., Grasser, N., & Snellen, I
Picos, D. G., de Regt, S., Gandhi, S., Grasser, N., & Snellen, I. A. G. 2025, A&A, 703, A65
work page 2025
- [57]
-
[58]
Rimmer, P. B. & Helling, C. 2016, ApJS, 224, 9
work page 2016
-
[59]
Rossini, F. D., Wagman, D. D., Evans, W. H., Levine, S., & Jaffe, I. 1952, Se- lected values of chemical thermodynamic properties, Tech. Rep. Circular 500, National Bureau of Standards, Washington, DC
work page 1952
-
[60]
Scott, M. G., Triaud, A. H. M. J., Barkaoui, K., et al. 2025, MNRAS, 540, 1909
work page 2025
-
[61]
Sharp, C. M. & Huebner, W. F. 1990, ApJS, 72, 417
work page 1990
-
[62]
Shenyavskaya, E. A. & Yungman, V . S. 2004, Journal of Physical and Chemical Reference Data, 33, 923
work page 2004
-
[63]
Stephenson, C. C. & Giauque, W. F. 1937, J. Chem. Phys., 5, 149
work page 1937
-
[64]
Stephenson, C. C. & Zettlemoyer, A. C. 1944, Journal of the American Chemical Society, 66, 1405
work page 1944
-
[65]
Stock, J. W., Kitzmann, D., & Patzer, A. B. C. 2022, MNRAS, 517, 4070
work page 2022
-
[66]
Stock, J. W., Kitzmann, D., Patzer, A. B. C., & Sedlmayr, E. 2018, MNRAS, 479, 865
work page 2018
-
[67]
Stull, D. R. 1947, Industrial & Engineering Chemistry, 39, 517
work page 1947
- [68]
- [69]
-
[70]
Verma, A., Goyal, J., Avarsekar, S., & Shukla, G. 2025, AJ, 170, 69
work page 2025
-
[71]
2006, The Astrophysical Journal, 648, 1181
Visscher, C., Lodders, K., & Fegley, B., J. 2006, The Astrophysical Journal, 648, 1181
work page 2006
-
[72]
Wagner, W. & Kretzschmar, H.-J. 2008, International Steam Tables: Proper- ties of Water and Steam Based on the Industrial Formulation IAPWS-IF97 (Springer-Verlag Berlin Heidelberg)
work page 2008
- [73]
-
[74]
Yaws, C. 1999, Chemical Properties Handbook: Physical, Thermodynamics, En- vironmental Transport, Safety & Health Related Properties for Organic &, Chemical engineering books (McGraw-Hill Education)
work page 1999
-
[75]
Zhang, Y ., Evans, J. R. G., & Yang, S. 2011, Journal of Chemical & Engineering Data, 56, 328 Article number, page 12 D. Kitzmann , J. W. Stock, and A.B.C. Patzer: FastChem 4: New chemical elements and improved convergence behaviour Appendix A: Overview of gas-phase and condensate species included in FastChem Table A.1.List of all gas-phase species includ...
work page 2011
-
[76]
FN, FNO, FNO2, FNO3, FNa,FNi i, FO, FOTi, OFO, FOO FP, FP –, FP+, FPS, FS, FS –, FS+, FSi, FTi,FZn i, F2, F2Fe, F2H2, F2H2Si, F2K–, F2K2, F2Mg, F2Mg+, F2N, cis-F2N2,trans-F 2N2 F2Na–, F2Na2, F2O, F2OS, F2OSi, F2OTi, F2O2, F2O2S, F2P, F2P–, F2P+, F2S, F2S–, F2S+, FS2F, SSF2 F2Si, F2Ti, F3Fe, F3HSi, F3H3, F3N, F3NO, F3OP, F3P, F3PS, F3S, F3S–, F3S+, F3Si, F...
work page 1973
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.