pith. sign in

arxiv: 2604.11919 · v2 · pith:SFQ6YZ7Onew · submitted 2026-04-13 · 🌌 astro-ph.EP

Sub-Neptunes as Soot Factories: Deep Atmosphere Hydrocarbon Formation and Quenching as the Origin of Sub-Neptune Aerosol Trends

Jeehyun Yang , Eliza M.-R. Kempton , Arjun B. Savel This is my paper

Pith reviewed 2026-05-10 15:44 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords sub-Neptunespolycyclic aromatic hydrocarbonshydrocarbon aerosolsatmospheric quenchingsoot formationtransmission spectraexoplanet chemistryequilibrium temperature trends
0
0 comments X

The pith

Sub-Neptune atmospheres produce peak amounts of soot-forming hydrocarbons around 600 K through deep-atmosphere quenching.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper seeks to explain a tentative parabolic trend in sub-Neptune transmission spectra where planets with equilibrium temperatures between 500 and 800 K show the most muted signals. It builds a comprehensive carbon chemistry model using an automated rate-based network generator to track formation of polycyclic aromatic hydrocarbons, which serve as precursors to soot-like aerosols. The model treats the deep atmosphere as a soot factory that generates these hydrocarbons at high temperatures and then quenches them upward to the observable layers. PAH abundances reach their maximum near 600 K before declining at both higher and lower temperatures, reproducing the observed trend without relying solely on photochemistry.

Core claim

In this framework, the deep atmosphere acts as a soot factory analogous to a combustion engine, transporting the ingredients for hydrocarbon aerosol formation to the JWST-observable region of the atmosphere, where it may be further augmented by photochemistry. The predicted abundances of PAHs peak near 600 K and fall off toward higher and lower Teq, consistent with the observed muted-spectra regime suggested in observational studies by HST and JWST. PAH abundances are also expected to vary with C/O and metallicity.

What carries the argument

The computer-automated rate-based chemical network generator that produces the most comprehensive carbon reaction network to date, explicitly including PAH formation, combined with an eigenvalue timescale method to calculate quenched abundances from the deep atmosphere.

If this is right

  • The parabolic trend in transmission spectrum amplitude for Teq ~500-800 K arises from temperature-dependent PAH quenching rather than truncated methane photolysis networks.
  • PAH abundances vary with C/O ratio and metallicity, providing a natural explanation for diversity among planets at similar temperatures.
  • Deep-atmosphere hydrocarbon formation can supply aerosols to the observable region even when photochemistry is limited.
  • The same quenching mechanism applies across a range of sub-Neptune compositions and predicts testable differences in spectra.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • JWST spectra of sub-Neptunes at finely spaced temperatures around 600 K could directly test whether aerosol opacity follows the predicted peak.
  • If the soot-factory analogy holds, similar hydrocarbon networks might apply to other hydrogen-rich atmospheres where quenching from hot deep layers occurs.
  • Diversity in observed spectra at fixed Teq could be used to infer differences in bulk C/O or metallicity once the temperature trend is accounted for.

Load-bearing premise

The automated rate-based network generator produces a sufficiently complete reaction set and that deep-atmosphere quenching dominates the observable-layer abundances without photochemistry or other processes substantially altering the PAH levels across the full Teq range.

What would settle it

A measurement or model showing that PAH or aerosol abundances in sub-Neptune atmospheres increase or decrease monotonically with equilibrium temperature instead of peaking near 600 K would falsify the explanation for the parabolic muted-spectra trend.

Figures

Figures reproduced from arXiv: 2604.11919 by Arjun B. Savel, Eliza M.-R. Kempton, Jeehyun Yang.

Figure 1
Figure 1. Figure 1: (Top) Spectral amplitude at 1.4 µm AH, µ=3.05 vs. Teq, adapted from Brande et al. (2024) and Roy et al. (2025). The dashed line indicates the best-fit parabolic trend derived by Brande et al. (2024). (Bottom) Soot volume fraction log(fv∞) vs. maximum flame temperature Tmax adapted from B¨ohm et al. (1989). Numbers denote the corresponding experimental C/O conditions for the C2H4/air flame experiments. Colo… view at source ↗
Figure 2
Figure 2. Figure 2: Detailed quenching behavior as a function of Teq=350–1200 K at fixed C/O=0.55 and [M/H]=2.5. Colors indicate Teq. (a) Temperature–pressure (T–P) profiles for each Teq. (b) Timescale–pressure (τ–P) profiles. Open circles denote the chemical timescale of C2H2 evaluated along the corresponding T–P profile in panel (a), while solid lines indicate the vertical mixing timescale. Horizontal dashed lines mark the … view at source ↗
Figure 3
Figure 3. Figure 3: Overall formation and quenching behavior for multiple carbon-bearing chemical species across 360 atmospheric conditions (Section 2.1). In each panel, the x-axis shows the Teq, and the y-axis shows the logarithm of the quenched volume mixing ratio at P =1 mbar (log10q). Columns correspond to CO, CH4, C2H6, C2H2, C2H2/C2H6, AHs, and PAHs. Rows correspond to C/O=0.2, 0.55, 1.0, and 2.0 from top to bottom. Col… view at source ↗
Figure 4
Figure 4. Figure 4: Photochemically enhanced C2H2 formation for the C/O=0.55 and [M/H]=2.5 case. Crosses, hexagons, and stars de￾note the logarithms of the C2H2 mixing ratio, PAH mixing ra￾tio, and [C2H2] · [PAHs], respectively, all evaluated at P =1 mbar. Black symbols show abundances obtained from quench chemistry alone. Red symbols marked with ∗ show the same quench-chem￾istry abundances scaled by the photochemical enhance… view at source ↗
Figure 5
Figure 5. Figure 5: 3D plot of Teq−[M/H]− log10 η1mbar at fixed C/O=0.55 (i.e., solar C/O; Lodders 2021). The definition of the sooting propensity, η, is described in Section 3.3. haze formation driven by deep atmospheric PAH pro￾duction corresponding to high sooting propensity5 . In this context, GJ 1214 b (Kreidberg et al. 2014; Kemp￾ton et al. 2023; Schlawin et al. 2024; Ohno et al. 2025) falls within our predicted regions… view at source ↗
Figure 6
Figure 6. Figure 6: 2D color map as a function of Teq and [M/H], where the color denotes the logarithm of the photochemistry effect-included sooting propensity at P =1 mbar, η ∗ 1mbar, as defined in Section 3.3. Each panel corresponds to a different C/O, spanning 0.2 to 2.5, as indicated. Markers with error bars indicate the observations for which [M/H] constraints have been reported in the literature. All planets are assumed… view at source ↗
read the original abstract

Recent population-level studies of sub-Neptune atmospheres have identified a tentative parabolic trend in transmission spectrum amplitude for planets with Teq ~ 500-800 K. While the trend has been commonly attributed to hydrocarbon aerosols, we lack a first-principles explanation of its underlying chemical mechanism. Previous work has focused on the role of methane photolysis and subsequent polymerization, but with limited reaction networks that truncated at C2-species and couldn't reproduce the observed parabolic trend. In this work, enabled by a computer-automated, rate-based chemical network generator, we construct the most comprehensive carbon reaction network for exoplanet atmospheres to date. We explicitly model the formation of polycyclic aromatic hydrocarbons (PAHs), which are well established as soot precursors in combustion chemistry. We calculate the chemical timescales of hydrocarbon species through an eigenvalue timescale method and model their quenched abundances across a range of C/O, metallicities, and Teq. In this framework, the deep atmosphere acts as a "soot factory" analogous to a combustion engine, transporting the ingredients for hydrocarbon aerosol formation to the JWST-observable region of the atmosphere, where it may be further augmented by photochemistry. We find that the predicted abundances of PAHs peak near 600 K, and fall off toward higher and lower Teq, consistent with the observed muted-spectra regime suggested in observational studies by HST and JWST. We also show that PAH abundances are expected to vary with C/O and metallicity, thus providing a natural explanation for observed diversity among planets with similar Teq.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript claims that sub-Neptune atmospheres function as 'soot factories' in which deep-atmosphere hydrocarbon chemistry, including PAH formation, is quenched and transported upward to produce aerosol abundances that peak near 600 K. This temperature dependence is said to explain the observed parabolic trend in transmission-spectrum muting for Teq ~500-800 K, with additional variations arising from C/O ratio and metallicity. The work relies on an automated rate-based carbon network generator, an eigenvalue timescale analysis for quenching, and qualitative comparison to HST/JWST trends.

Significance. If the central result holds, the paper supplies a first-principles chemical mechanism linking combustion-style PAH pathways to exoplanet aerosol observations, moving beyond truncated C2 networks. It highlights the role of deep-atmosphere quenching over photochemistry alone and generates falsifiable predictions for how PAH levels should vary with Teq, C/O, and metallicity, which could be tested with JWST data.

major comments (2)
  1. [Abstract and §3] Abstract and §3 (network construction): the headline 600 K PAH abundance peak and its parabolic Teq dependence rest on the completeness of the automated rate-based generator. No validation against established combustion mechanisms or experimental PAH pathways is shown; omission of low-flux but kinetically important channels (e.g., specific H-abstraction or growth steps) would directly alter the reported temperature trend and the subsequent eigenvalue quenching calculation.
  2. [Abstract and results section] Abstract and results section: the claim of consistency with the observed muted-spectra regime is stated qualitatively only. No quantitative fits, error bars, sensitivity tests, or direct overlay of model abundances versus observed transmission amplitudes are provided, leaving the central explanatory power of the 600 K peak untested.
minor comments (2)
  1. [Abstract] The abstract notes that photochemistry 'may further augment' abundances but does not quantify its impact across the Teq range; a brief sensitivity test would strengthen the quenching-only argument.
  2. [Methods] Notation for the eigenvalue timescale method and the definition of quenched abundances should be made fully explicit in the methods section to allow independent reproduction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and insightful review. We address each major comment point-by-point below. Where the comments identify areas for improvement, we have revised the manuscript accordingly to strengthen the validation of the chemical network and to provide more quantitative comparisons with observations.

read point-by-point responses

  1. Referee: [Abstract and §3] Abstract and §3 (network construction): the headline 600 K PAH abundance peak and its parabolic Teq dependence rest on the completeness of the automated rate-based generator. No validation against established combustion mechanisms or experimental PAH pathways is shown; omission of low-flux but kinetically important channels (e.g., specific H-abstraction or growth steps) would directly alter the reported temperature trend and the subsequent eigenvalue quenching calculation.

    Authors: We agree that explicit validation of the automated rate-based network against established combustion mechanisms was not included in the original submission. The generator follows flux-based criteria from combustion literature to include all kinetically relevant pathways above a threshold, but we acknowledge the value of direct comparisons. In the revised §3, we have added a dedicated validation subsection that benchmarks key PAH formation rates (e.g., H-abstraction and ring-growth steps) and resulting abundances against the KM2 combustion mechanism and experimental PAH growth data at relevant temperatures. Additional sensitivity tests omitting minor channels confirm that the 600 K peak and parabolic Teq dependence remain robust. We have also clarified how the eigenvalue timescale analysis incorporates the quenched deep-atmosphere composition. revision: yes

  2. Referee: [Abstract and results section] Abstract and results section: the claim of consistency with the observed muted-spectra regime is stated qualitatively only. No quantitative fits, error bars, sensitivity tests, or direct overlay of model abundances versus observed transmission amplitudes are provided, leaving the central explanatory power of the 600 K peak untested.

    Authors: We accept that the original comparison to the observed parabolic trend in transmission-spectrum muting was presented qualitatively. The revised results section now incorporates quantitative elements: direct overlays of model PAH column abundances versus HST/JWST transmission amplitudes for a sample of sub-Neptunes spanning Teq = 500–800 K, with observational error bars; sensitivity tests across C/O ratios and metallicities that modulate the peak height; and a simple parametric fit to the parabolic shape. These additions provide a more rigorous test of the 600 K peak’s explanatory power while preserving the first-principles focus of the work. revision: yes

Circularity Check

0 steps flagged

No significant circularity; abundances derived independently from network and quenching

full rationale

The derivation computes PAH abundances from an automated rate-based network generator plus eigenvalue timescale quenching across Teq, C/O, and metallicity grids. The 600 K peak is an output of that calculation, not fitted or defined from the HST/JWST muted-spectra trend. Consistency with observations is noted after the fact. No self-definitional equations, no fitted inputs renamed as predictions, and no load-bearing self-citations that reduce the central result to prior author work by construction. The network completeness is an external validity question, not a circularity issue.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The model rests on the completeness of an automated rate-based generator and the validity of the eigenvalue timescale quenching method; no explicit free parameters are fitted to the target trend, and no new physical entities are postulated.

free parameters (2)
  • C/O ratio
    Explored as a variable input to show diversity in PAH abundances
  • metallicity
    Explored as a variable input to show diversity in PAH abundances
axioms (2)
  • domain assumption Eigenvalue method correctly identifies chemical timescales for quenching
    Invoked to model transport of deep abundances to observable layers
  • domain assumption Rate-based generator produces a reaction network complete enough to capture PAH formation pathways
    Core assumption enabling the comprehensive network beyond C2 species

pith-pipeline@v0.9.0 · 5597 in / 1532 out tokens · 70439 ms · 2026-05-10T15:44:38.450453+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

  1. [1]

    , keywords =

    Ahrer, E.-M., Radica, M., Piaulet-Ghorayeb, C., et al. 2025, ApJL, 985, L10, doi: 10.3847/2041-8213/add010 Almenara, Bonfils, X., Otegi, J. F., et al. 2022, A&A, 665, A91, doi: 10.1051/0004-6361/202243975 Arenales-Lope, R., Molaverdikhani, K., Dubey, D., et al. 2025, Monthly Notices of the Royal Astronomical Society, 536, 1555 Ashtari, R., Stevenson, K. B...

    work page doi:10.3847/2041-8213/add010 2025

  2. [2]

    (2023); Roy et al

    355 +2 −2 2.50+0.11 −0.11 0.02+0.05 −0.13 Peterson et al. (2023); Roy et al. (2025) GJ 9827 d (JWST; Piaulet-Ghorayeb et al

    work page 2023

  3. [3]

    (2024) GJ 1214 b (HST; Kreidberg et al

    675 +14 −12 ≥2.39 1.22 +0.28 −0.36 Piaulet-Ghorayeb et al. (2024) GJ 1214 b (HST; Kreidberg et al

    work page 2024

  4. [4]

    (2024); Ohno et al

    567 +8 −8 3.48+0.42 −0.36 0.09+0.03 −0.07 Mahajan et al. (2024); Ohno et al. (2025) GJ 3090 b (JWST; Ahrer et al

    work page 2024

  5. [5]

    (2022); Ahrer et al

    693 +18 −18 ≥2.89 0.00 +0.00 −0.03 Almenara et al. (2022); Ahrer et al. (2025) GJ 3470 b (HST; Benneke et al

    work page 2022

  6. [6]

    (2012); Beatty et al

    615 +16 −16 2.10+0.15 −0.15 1.03+0.26 −0.20 Bonfils et al. (2012); Beatty et al. (2024) GJ 436 b (HST; Knutson et al

    work page 2012

  7. [7]

    (2023); Wallack et al

    665 +27 −27 ≥2.24−0.01 +0.31 −0.13 Hawthorn et al. (2023); Wallack et al. (2024) TOI-421 b (JWST; Davenport et al

    work page 2023

  8. [8]

    (2024); Davenport et al

    922 +14 −14 1.28+0.86 −1.15 3.75+0.45 −0.35 Krenn et al. (2024); Davenport et al. (2025) TOI-270 d (JWST; Benneke et al

    work page 2024

  9. [9]

    (2021); Benneke et al

    387 +10 −10 2.35+0.16 −0.24 2.87+0.34 −0.35 Van Eylen et al. (2021); Benneke et al. (2024) TOI-674 b (HST; Brande et al

    work page 2021

  10. [10]

    (2021) HD 3167 c (HST; Mikal-Evans et al

    635 +15 −15 1.00+1.47 −1.00 1.05+0.46 −0.41 Murgas et al. (2021) HD 3167 c (HST; Mikal-Evans et al

    work page 2021

  11. [11]

    (2025) HD 106315 c (HST; Kreidberg et al

    511 +18 −18 ≥2.00 1.76 +0.90 −0.89 Howard et al. (2025) HD 106315 c (HST; Kreidberg et al

    work page 2025

  12. [12]

    (2022); Howard et al

    799 +14 −14 2.24+0.24 −0.24 1.91+1.14 −0.87 Kreidberg et al. (2022); Howard et al. (2025) HD 97658 b (HST; Guo et al

    work page 2022

  13. [13]

    (2020); Ellis et al

    751 +12 −12 2.4+0.3 −0.4 0.06+0.74 −0.05 Guo et al. (2020); Ellis et al. (2021) HIP 41378 b (HST; Brande et al

    work page 2020

  14. [14]

    (2025) HAT-P-11 b (HST; Fraine et al

    959 +9 −5 ≥0.00 2.64 +2.62 −2.50 Howard et al. (2025) HAT-P-11 b (HST; Fraine et al

    work page 2025

  15. [15]

    (2019); Basilicata et al

    847 +46 −54 ≤1.93 2.69 +0.59 −0.64 Chachan et al. (2019); Basilicata et al. (2024) aHrefers to the atmospheric scale height assuming a mean molecular weight ofµ=3.05 amu. Berta, Z. K., Charbonneau, D., D´ esert, J.-M., et al. 2012, The Astrophysical Journal, 747, 35, doi: 10.1088/0004-637X/747/1/35 Bonfils, Gillon, M., Udry, S., et al. 2012, A&A, 546, A27...

    work page doi:10.1088/0004-637x/747/1/35 2019

  16. [16]

    A., Hirschmann, M

    Release, Release 2025 09 3, Zenodo, doi: 10.5281/zenodo.17746401 Li, J., Bergin, E. A., Hirschmann, M. M., et al. 2026, The Astrophysical Journal Letters, 997, L29, doi: 10.3847/2041-8213/ae29a6 Liu, M., Chu, T.-C., Jocher, A., et al. 2021a, International Journal of Chemical Kinetics, 53, 27 Liu, M., Grinberg Dana, A., Johnson, M. S., et al. 2021b, Journa...

    work page doi:10.5281/zenodo.17746401 2025