Sonora Elf Owl: Disequilibrium Chemistry Model

• Sonora Elf Owl is a cloud-free, 1D disequilibrium-chemistry model that simulates substellar atmospheres by incorporating free parameters like metallicity, C/O ratio, and eddy diffusion (K_zz).
• It employs self-consistent radiative–convective equilibrium via PICASO 3.0, extending previous grids by coupling quench chemistry and secondary CO2 corrections over a broad temperature and gravity range.
• The model spans effective temperatures of 275–2400 K and has been validated through spectral fittings with JWST and SpeX, while noting limitations in reproducing cloudy L-dwarf spectra.

Sonora Elf Owl is the disequilibrium-chemistry branch of the Sonora substellar atmosphere model family, designed for brown dwarfs and self-luminous giant exoplanets. It comprises self-consistent, cloud-free, one-dimensional radiative–convective equilibrium atmospheres built with PICASO 3.0, with free metallicity, carbon-to-oxygen ratio, and eddy diffusion coefficient $K_{\rm zz}$, and was developed for spectra in regimes where vertical mixing drives departures from local chemical equilibrium. Within the Sonora series, it extends earlier cloud-free grids by coupling composition variation and quench chemistry over a broad parameter space, and version 2 further corrected the treatment of disequilibrium CO$_2$ and removed PH$_3$ opacity from the released spectra (Mukherjee et al., 2024, Wogan et al., 6 May 2025).

1. Position within the Sonora model family

Sonora Elf Owl is the fourth major release in the Sonora Substellar Atmosphere Models series. Sonora Bobcat provides cloud-free rainout equilibrium atmospheres spanning $200$–$2400$ K with metallicity and C/O variations but no disequilibrium chemistry. Sonora Cholla adds cloud-free disequilibrium chemistry at solar composition for approximately $500$–$1300$ K. Sonora Elf Owl extends disequilibrium chemistry self-consistently over $275$–$2400$ K, multiple gravities, seven orders of magnitude in $K_{\rm zz}$, sub- to super-solar metallicities, and C/O from oxygen-rich to carbon-rich. Opacity mixing is performed “on-the-fly” for all gases, which improves self-consistency relative to previous grids.

Its intended regime is the cloud-free low-temperature atmosphere domain, especially mid- and late-T dwarfs and cold Y dwarfs, where condensates are largely absent from the photosphere. Later applications explicitly describe it as a cloud-free grid optimized especially for T–Y dwarfs. The same cloud-free assumption also defines its principal failure mode: it does not reproduce dusty L-dwarf spectral energy distributions or the L/T transition, where condensate opacity and cloud clearing modify the continuum shape and band depths. In those regimes, studies comparing Sonora Elf Owl with cloudy alternatives report that cloud-free models compensate by pushing metallicity and C/O toward grid edges or by returning physically problematic surface gravities, while residuals remain large (Mukherjee et al., 2024, Lueber et al., 26 May 2025, Tu et al., 2 Oct 2025).

2. Physical framework and disequilibrium chemistry

The grid is formulated as a $_2$0D, plane-parallel atmosphere with $_2$1 pressure layers under hydrostatic equilibrium,

$_2$2

PICASO iterates the temperature–pressure profile until radiative–convective equilibrium is achieved. Convective energy transport is handled by mixing-length theory inside the equilibrium solver, while disequilibrium abundances are imposed through vertical mixing parameterized by a constant-with-pressure eddy diffusion coefficient $_2$3 for each model.

Tracer transport is represented by an eddy-diffusion flux,

$_2$4

where $_2$5 is the total number density and $_2$6 the mixing ratio of species $_2$7. The vertical mixing timescale is

$_2$8

and quenching occurs where

$_2$9

Below the quench pressure, equilibrium is maintained; above it, abundances are frozen in.

The disequilibrium network includes CH$_3$0, CO, CO$_3$1, NH$_3$2, N$_3$3, HCN, PH$_3$4, and H$_3$5O in the original release. The governing net reactions include

$_3$6

$_3$7

$_3$8

$_3$9

$200$0

$200$1

and, for phosphorus,

$200$2

Chemical timescales are taken from Zahnle & Marley (2014), with Visscher & Fegley formulations for phosphorus chemistry. Equilibrium abundances below quench levels are obtained from precomputed tables across [M/H] and C/O, using solar elemental abundances from Lodders (2009), and C/O is varied by holding C+O fixed and changing their ratio before metallicity scaling so that C/O and [M/H] are not conflated (Mukherjee et al., 2024).

3. Grid architecture and parameterization

The full Sonora Elf Owl grid spans effective temperature $200$3–$200$4 K, surface gravity $200$5–$200$6 in cgs, $200$7, metallicities [M/H] $200$8, and C/O $200$9–$2400$0. The temperature sampling is $2400$1 K from $2400$2–$2400$3 K, $2400$4 K from $2400$5–$2400$6 K, and $2400$7 K from $2400$8–$2400$9 K. The atmosphere files provide $500$0 structures, equilibrium and quenched abundance profiles, and emission spectra from $500$1–$500$2 at $500$3, for a total of $500$4 models.

Published fitting studies use subsets or instrument-matched interpolants rather than the raw grid alone. In a comparative study of low-temperature atmosphere fitting methodologies, the working Sonora Elf Owl subset covered $500$5–$500$6 K, $500$7–$500$8, [M/H] $500$9 to $1300$0, C/O $1300$1–$1300$2 relative to solar, and $1300$3–$1300$4, with $1300$5 models after down-selection. In the JWST study of $1300$6 cold brown dwarfs, forward modeling used linear interpolation across eight fitted parameters: $1300$7, $1300$8, [M/H], C/O, $1300$9, two MIRI wavelength-correction parameters, and a global flux-scaling parameter $275$0. In the later JWST/NIRSpec distant-brown-dwarf survey, Sonora Elf Owl version 2 was interpolated and convolved to PRISM/CLEAR resolution over $275$1–$275$2, with flat priors on the physical parameters and a logarithmic prior on the scale factor (Mukherjee et al., 2024, Tu et al., 2024, Tu et al., 2 Oct 2025).

4. Secondary CO$275$3 quenching and the version 2 correction

Version 1 applied quench chemistry species-by-species: for each species $275$4, a quench pressure $275$5 was obtained from $275$6, the atmosphere was assumed to remain in chemical equilibrium for $275$7, and the volume mixing ratio for $275$8 was held constant at the equilibrium value evaluated at $275$9. This procedure worked well for species such as CO, CH$2400$0, and NH$2400$1, but it systematically underpredicted CO$2400$2 because the observable CO$2400$3 abundance is controlled by equilibrium with CO and H$2400$4O, while CO commonly quenches deeper than CO$2400$5.

The controlling reaction is

$2400$6

with mass-action relation

$2400$7

and

$2400$8

For the CO$2400$9CH$K_{\rm zz}$0 system, the characteristic timescale used to locate the CO quench level is

$K_{\rm zz}$1

with $K_{\rm zz}$2 in seconds, $K_{\rm zz}$3 in bars, and $K_{\rm zz}$4 in Kelvin. The version 2 correction therefore adopts “secondary quenching”: first determine the CO and H$K_{\rm zz}$5O quench level, freeze their abundances there, then compute CO$K_{\rm zz}$6 in equilibrium with those already-quenched abundances until CO$K_{\rm zz}$7 itself quenches, after which CO$K_{\rm zz}$8 is frozen.

Operationally, version 2 implements four steps. First, compute $K_{\rm zz}$9 from $_2$00 and $_2$01 and locate quench pressures for CO/H$_2$02O and CO$_2$03. Second, at $_2$04, set $_2$05 and $_2$06 above that level equal to their equilibrium values at $_2$07. Third, for pressures between $_2$08 and $_2$09, set

$_2$10

Fourth, freeze CO$_2$11 at its own quench level:

$_2$12

The quantitative effect is substantial for the coldest objects. The correction increases CO$_2$13 by over one order of magnitude for objects with $_2$14 K and by up to two orders of magnitude in some cases. In an illustrative $_2$15 K brown dwarf model with $_2$16, metallicity $_2$17 solar, solar C/O, and vertically constant $_2$18, a full-kinetics Photochem calculation shows CO quenching near $_2$19 bar and CO$_2$20 quenching near $_2$21 bar, with CO$_2$22 freezing at $_2$23. In that case, version 1 underpredicted CO$_2$24 above $_2$25 bar by roughly two orders of magnitude. Version 2 also removes PH$_2$26 opacity because observations indicated that version 1 consistently produced too much PH$_2$27 absorption. The updated abundances and spectra were posted as new versions of the original Zenodo records, but the radiative–convective climate simulations were not re-run, so the published version 2 pressure–temperature profiles are not strictly self-consistent with the corrected composition (Wogan et al., 6 May 2025).

5. Spectral behavior, parameter sensitivities, and degeneracies

The influence of $_2$28 on the thermal structure and emergent spectrum depends strongly on metallicity and $_2$29. At solar metallicity, vigorous mixing with $_2$30 cools the $_2$31 profile by approximately $_2$32 K relative to chemical equilibrium. At [M/H] $_2$33, the same mixing can cool the profile by approximately $_2$34–$_2$35 K. At [M/H] $_2$36, the effect on $_2$37 is minimal because the atmosphere is more transparent and disequilibrium versus equilibrium optical depths differ by less than an order of magnitude across most wavelengths. The temperature range of strongest sensitivity also shifts with composition: metal-poor atmospheres are most $_2$38-sensitive at approximately $_2$39–$_2$40 K, solar-composition atmospheres at approximately $_2$41–$_2$42 K, and super-solar atmospheres at approximately $_2$43–$_2$44 K.

Key wavelength regions encode different combinations of physics. The $_2$45–$_2$46 region is dominated by CH$_2$47; increasing $_2$48 weakens CH$_2$49 bands by lowering CH$_2$50 aloft, higher metallicity can also weaken them, and high C/O deepens them. The $_2$51–$_2$52 region contains CO$_2$53, which grows strongly with metallicity and has weaker dependence on $_2$54 and C/O; this is the principal band for breaking the $_2$55–[M/H] degeneracy. The $_2$56–$_2$57 region contains the CO fundamental and is degenerate because both higher $_2$58 and higher metallicity enhance CO absorption. The $_2$59–$_2$60 region is dominated by H$_2$61O and alkali line wings and is more sensitive to metallicity than to $_2$62. The $_2$63–$_2$64 region probes NH$_2$65 and becomes important in combined near- and mid-infrared analyses.

Synthetic fitting tests define a practical observing strategy. M band alone at $_2$66 generally cannot uniquely constrain $_2$67, [M/H], C/O, or $_2$68. Coverage from $_2$69–$_2$70 can constrain $_2$71, $_2$72, [M/H], and C/O, but $_2$73 remains loose. Coverage from $_2$74–$_2$75 constrains all parameters, including gravity, and JWST NIRSpec Prism with broad $_2$76–$_2$77 coverage yields the tightest posteriors in the published tests. In the original nine-object T-dwarf application, inferred $_2$78 values were consistently low, approximately $_2$79–$_2$80, and were interpreted as evidence that the observations probe quenching in a deep detached radiative zone rather than vigorous convective transport (Mukherjee et al., 2024).

6. Inference workflows and observational performance

Sonora Elf Owl has been used in both direct grid fitting and machine-learning-assisted inference. In benchmark SpeX prism analyses at $_2$81–$_2$82 and $_2$83, one study compared a Markov Chain Monte Carlo algorithm that interpolates across spectral fluxes with a Random Forest Retrieval trained on the grid. The MCMC method uses a $_2$84 statistic with an optimal multiplicative scale factor $_2$85, linear interpolation in log flux over logarithmic parameter axes, a Metropolis–Hastings acceptance scheme, and posterior sampling under

$_2$86

The Random Forest Retrieval uses absolute-calibrated flux densities as features and predicts $_2$87, $_2$88, [M/H], C/O, $_2$89, and a flux-scaling parameter $_2$90, with ensemble averaging

$_2$91

That study found the MCMC approach to yield higher fit quality and more precise parameters, while the Random Forest Retrieval was orders of magnitude faster after training. The recommended workflow was therefore rapid multi-grid screening with the random forest followed by MCMC refinement on the chosen grid. In that benchmark sample, Sonora Elf Owl was generally optimal for mid- and late-T dwarfs and not preferred for L dwarfs or many early-T dwarfs, because the absence of cloud opacity is decisive in the latter regimes (Lueber et al., 26 May 2025).

JWST applications confirm both the utility and the systematic structure of the grid. In a sample of $_2$92 T and Y dwarfs fitted jointly to NIRSpec CLEAR/PRISM and MIRI LRS spectra, Sonora Elf Owl and ATMO2020++ returned relatively consistent effective temperatures, typically differing by tens of Kelvin, and both yielded metallicities close to solar values for nearby objects. The largest disagreement was in surface gravity: Sonora Elf Owl values were typically about $_2$93 dex lower than ATMO2020++, with $_2$94 of $_2$95 objects at the lower Elf Owl grid boundary $_2$96 and $_2$97 objects having $_2$98. Radii inferred from the scaling parameter and parallax were typically $_2$99–$_3$00, and the corresponding masses and ages derived from evolutionary tracks indicated that most objects were below $_3$01 and younger than $_3$02 Gyr, with Y dwarfs at $_3$03–$_3$04 and $_3$05–$_3$06 Gyr. The same study found that NIRSpec-only fits usually recovered $_3$07 and $_3$08 close to the combined-spectrum values, whereas MIRI-only fits gave more uncertain and often lower gravities unless $_3$09 was fixed (Tu et al., 2024).

A later JWST/NIRSpec PRISM/CLEAR survey of $_3$10 public spectra used Sonora Elf Owl version 2 alongside LOWZ and SAND. In that pipeline, the spectra were truncated at the blue end where $_3$11, cross-calibrated to NIRCam photometry to mitigate MOS mis-centering, convolved to the wavelength-dependent PRISM resolution, and analyzed with UltraNest nested sampling using a jitter term in the likelihood. Sonora Elf Owl version 2 performed strongly for cool T and Y dwarfs, but cloud-free grids failed to reproduce L/T transition spectra, whereas SAND provided a better match in those cases. For the metallicity-gradient analysis in that survey, Sonora-based results showed no evident trend of [M/H] with Galactic height $_3$12 after excluding cluster sources, L/T transition objects, and objects with $_3$13 K (Tu et al., 2 Oct 2025).

7. Limitations, systematic biases, and future development

The main limitations of Sonora Elf Owl follow directly from its simplifying assumptions. It is a cloud-free, one-dimensional grid with constant-with-pressure $_3$14 in each model. Real atmospheres can have pressure-dependent mixing that differs by orders of magnitude between radiative and convective zones, and clouds can alter both the thermal profile and quenched abundances by warming deep adiabats. These omissions are central to the grid’s poor performance for L dwarfs, the L/T transition, and some early-T objects.

Several systematic biases recur across applications. Surface gravities inferred with Sonora Elf Owl are often lower than those from ATMO2020++ or from evolutionary expectations derived from benchmark primaries. In the JWST $_3$15-object study, this discrepancy was attributed to differences in convective treatment and opacity databases, particularly the ATMO2020++ use of an effective adiabatic index $_3$16, which reddens colors and lowers J-band peaks in a way partly degenerate with gravity. In benchmark SpeX fitting, late-T dwarf gravities from Sonora Elf Owl were again systematically low relative to expectations, a behavior the authors associated with pressure-dependent chemistry or collision-induced H$_3$17 absorption modeling. For very cold objects, metallicities derived from NIRSpec-only Sonora fits can differ from those obtained with combined NIRSpec+MIRI coverage, and both C/O and $_3$18 are especially sensitive to interpolation error and spectral mismatch.

The uncertainties quoted in fits are also often optimistic if only observed flux errors are propagated. Published interpolation tests show that most parameters behave well under linear interpolation, but the $_3$19–$_3$20 region is more nonlinear because CO$_3$21, CO, CH$_3$22, and H$_3$23O overlap there, and biases are strongest for $_3$24 and C/O. The large NIRSpec survey therefore added half-grid spacing as an uncertainty term and used a free jitter parameter; including that jitter increased posterior uncertainties by approximately $_3$25–$_3$26 depending on the parameter.

Version 2 resolves the specific CO$_3$27 underprediction in the original release and removes PH$_3$28 opacity, but it does not yet provide fully self-consistent recomputed pressure–temperature profiles. The correction note states that a fully self-consistent grid incorporating corrected CO$_3$29 will be released in future work. Other suggested developments across the literature include more finely sampled grids, improved interpolation or emulation, pressure-sensitive opacity revisions, methane-chemistry updates, expansion or validation of metallicity coverage for NIRSpec-only analyses, and incorporation of condensate physics across the L/T transition (Wogan et al., 6 May 2025, Lueber et al., 26 May 2025, Tu et al., 2024).