pith. machine review for the scientific record. sign in

arxiv: 2604.12132 · v1 · submitted 2026-04-13 · 🌌 astro-ph.SR · astro-ph.EP

Recognition: unknown

Inspecting Cloudy Substellar Atmospheres with JWST MIRI Synthetic Magnitudes from Spitzer Mid-infrared Spectra

Jolie LHeureux (1 , 2) , Genaro Su\'arez (3) , Johanna M. Vos (4) , Stanimir Metchev (5 , 6) , Jacqueline K. Faherty (3) , Sherelyn Alejandro Merchan (1
show 26 more authors
3) Kelle L. Cruz (1 3 7) ((1) Department of Physics Graduate Center City University of New York New York USA (2) Department of Astronomy Columbia University (3) Department of Astrophysics American Museum of Natural History (4) School of Physics Trinity College Dublin Dublin Ireland (5) Department of Physics Astronomy Western University London Canada (6) Institute for Earth Space Exploration (7) Department of Physics Hunter College USA)
Authors on Pith no claims yet

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

classification 🌌 astro-ph.SR astro-ph.EP
keywords brown dwarfsL dwarfssilicate cloudsmid-infrared photometryJWST MIRISpitzer spectrasynthetic magnitudessubstellar atmospheres
0
0 comments X

The pith

L dwarfs with mF770W - mF1000W < 0.03 mag are seven times more likely to host cloudy atmospheres.

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

The authors calculate synthetic magnitudes in four JWST MIRI filters from Spitzer mid-infrared spectra of 113 M5-T9 ultracool dwarfs. They find that color-color diagrams built from the F770W and F1000W bands, which bracket the 9 micron silicate absorption, separate L dwarfs with silicate clouds from cloud-free ones most cleanly. Objects falling below the 0.03 mag threshold in that color are seven times more likely to be classified as cloudy. The paper concludes that MIRI photometry therefore offers a fast way to pick cloudy targets for spectroscopic follow-up. Atmospheric models still fail to match the observed positions, especially the strength of the silicate feature itself.

Core claim

Synthetic photometry in the JWST MIRI F560W, F770W, F1000W, and F1280W filters derived from Spitzer spectra shows that F770W–F1000W color diagrams best isolate L-type objects whose photospheres contain silicate clouds; L dwarfs with mF770W − mF1000W < 0.03 mag are seven times more likely to be cloudy than those above the cut. Models underpredict the ~9 μm silicate signature, although some cloudy Sonora Diamondback grids reproduce the observed trends more closely.

What carries the argument

The F770W and F1000W MIRI filters that straddle the ~9 μm silicate absorption feature, turned into a color metric that correlates with the presence of clouds in L dwarf atmospheres.

If this is right

  • MIRI photometry alone can pre-select likely cloudy L dwarfs for targeted spectroscopy, reducing the fraction of JWST time spent on cloud-free targets.
  • The same filter pair can be applied to new or archival MIRI imaging to expand the known sample of cloudy substellar atmospheres.
  • Model grids must be revised to produce a stronger 9 μm silicate feature if they are to reproduce the locations of warm cloudy brown dwarfs.
  • Diagrams using F1000W and F1280W add little separation power because the spectra are weaker at longer wavelengths.

Where Pith is reading between the lines

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

  • If the color cut holds in real JWST data, it could be used in wide-field MIRI surveys to estimate the overall fraction of cloudy L dwarfs.
  • The same photometric approach might be tested on directly imaged giant planets whose atmospheres are thought to contain similar silicate clouds.
  • A direct comparison of synthetic versus observed MIRI magnitudes for the same objects would quantify any zero-point offsets introduced by the conversion process.

Load-bearing premise

The independent prior classification of each object as cloudy or cloud-free does not itself depend on the mid-infrared colors being tested.

What would settle it

Actual JWST MIRI photometry for a sample of L dwarfs whose cloudy or cloud-free status has already been established by other means, checked against the same 0.03 mag threshold.

Figures

Figures reproduced from arXiv: 2604.12132 by 2), (2) Department of Astronomy, 3, 3), (3) Department of Astrophysics, (4) School of Physics, (5) Department of Physics, 6), (6) Institute for Earth, 7) ((1) Department of Physics, (7) Department of Physics, American Museum of Natural History, Astronomy, Canada, City University of New York, Columbia University, Dublin, Genaro Su\'arez (3), Graduate Center, Hunter College, Ireland, Jacqueline K. Faherty (3), Johanna M. Vos (4), Jolie LHeureux (1, Kelle L. Cruz (1, London, New York, Sherelyn Alejandro Merchan (1, Space Exploration, Stanimir Metchev (5, Trinity College Dublin, USA, USA), Western University.

Figure 1
Figure 1. Figure 1: Spitzer IRS spectrum of 2148+4003 (red curve) and photometry from IRAC Ch4 (green), WISE W3 (dark purple), and synthetic fluxes using the diverse MIRI (light purple) filters in the bottom panel. The main spectral fea￾tures in the spectrum are indicated. The F560W filter covers the water absorption and the F770W, F1000W, and F1280W filters are centered, respectively, before, within, and after the silicate a… view at source ↗
Figure 2
Figure 2. Figure 2: Subset of 23 of the 113 Spitzer IRS spectra analyzed in this study for late-M, L, and T dwarfs. Our derived MIRI synthetic fluxes using the F560W, F770W, F1000W, and F1280W filters are overplotted, as indicated in the label. The subset contains four M, 10 L, and 9 T-type dwarfs. Key molecular signatures are indicated, namely water (H2O), silicates, methane (CH4), and ammonia (NH3), which are present, respe… view at source ↗
Figure 3
Figure 3. Figure 3: MF770W vs. mF770W − mF1000W (top left panel), MF1000W vs. mF770W − mF1000W (top right panel), MF1000W vs. mF1000W −mF1280W (bottom left panel), and MF1280W vs. mF1000W −mF1280W (bottom right panel) CMDs for ultracool dwarfs using MIRI synthetic magnitudes. The blue curve represents the median position of dwarfs with similar spectral types (bins of 2 spectral subtypes). Cloudy atmospheres are indicated by b… view at source ↗
Figure 6
Figure 6. Figure 6: Spitzer IRS spectrum of WISE J041521.21- 093500.6 along with synthetic MIRI F560W, F770W, F1000W, F1280W, F1500W, and F1800W fluxes in red and observed MIRI F1000W, F1280W, and F1800W fluxes in black. 5. SUMMARY We presented mid-infrared JWST MIRI CMDs and a CCD to identify the location of objects with silicate clouds in their atmospheres and tested how well model sequences predict the evolution of the war… view at source ↗
Figure 5
Figure 5. Figure 5: mF1000W −mF1280W vs. mF770W −mF1000W CCD for the late-M to T dwarfs with Spitzer IRS spectra. Ele￾ments in the figure are the same as in [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
read the original abstract

We examine the positions of substellar objects in mid-infrared color-magnitude and color-color diagrams to distinguish between cloudy and cloud-free atmospheres. Using Spitzer mid-infrared spectra of 113 M5-T9 ultracool dwarfs, we derive synthetic photometry for the JWST MIRI F560W, E'770W, F1000W, and F 1280W filters, which cover key absorption features including the ~9 um silicate signa-ture. We find that diagrams involving F770W and F1000W best separate L-type objects with silicate clouds in their photospheres. L dwarfs with mE77ow - mF1000w < 0.03 mag are seven times more likely to host cloudy atmospheres. Diagrams using F1000W and F1280W are less informative due to the lower signal of the spectra at long wavelengths. Current model predictions struggle to reproduce the positions of cloudy, warm brown dwarfs, likely because atmospheric models underestimate the ~9 um silicate feature. Cloudy Sonora Diamondback models better match the observed trends, although this may reflect improvements capturing indirect effects of clouds on the 6.25 um water absorption feature rather than accurately modeling the silicate feature itself. Our analysis indicates that JWST MIRI photometry can efficiently identify new cloudy extrasolar atmospheres for targeted spectroscopic follow-up, optimizing the use of telescope time.

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

3 major / 3 minor

Summary. The manuscript derives synthetic JWST MIRI magnitudes in the F560W, F770W, F1000W, and F1280W filters from Spitzer IRS spectra of 113 M5-T9 ultracool dwarfs. It analyzes color-color and color-magnitude diagrams to identify photometric criteria that separate L dwarfs with silicate clouds from cloud-free objects. The central result is that L dwarfs with mF770W − mF1000W < 0.03 mag are seven times more likely to host cloudy atmospheres, and that F770W/F1000W diagrams best separate such objects. The paper also compares the data to atmospheric models and concludes that JWST MIRI photometry can efficiently identify new cloudy atmospheres for follow-up spectroscopy.

Significance. If the cloudy/cloud-free classifications are independent of the mid-IR spectra analyzed here, the empirical color cut provides a practical, model-independent tool for selecting JWST targets. The sample size of 113 objects lends statistical weight to the likelihood ratio, and the finding that models underpredict the silicate feature is useful for model development. This could optimize JWST observing time by prioritizing photometric candidates for spectroscopy.

major comments (3)
  1. [§3] §3 (or equivalent section defining sample classifications): The abstract states that F770W/F1000W 'best separate' objects 'with silicate clouds in their photospheres' and reports a factor-of-seven likelihood, but does not specify how the binary cloudy vs. cloud-free labels were assigned. If these labels were derived from the strength or presence of the ~9 μm silicate feature in the same Spitzer spectra used to compute the synthetic magnitudes, the color separation is tautological. Please add an explicit description of the labeling procedure (e.g., based on near-IR spectra, variability, or independent mid-IR studies) and demonstrate that it is independent of the F770W–F1000W metric.
  2. [§4.2] §4.2 (Likelihood ratio and color cut): The claim that objects with mF770W − mF1000W < 0.03 mag are 'seven times more likely' to host cloudy atmospheres requires the exact statistical definition. Is this the ratio of conditional fractions, a likelihood ratio test, or a simple count ratio? Provide the formula, the number of objects in each bin, and any uncertainty estimate (e.g., from bootstrapping or binomial errors) so the robustness of the factor of seven can be assessed.
  3. [§5] §5 (Model comparison): The statement that 'Cloudy Sonora Diamondback models better match the observed trends' is presented without a quantitative metric. Specify whether this is assessed via reduced χ², residual rms, or visual inspection, and include a table comparing residuals or goodness-of-fit statistics for the different model grids (e.g., Sonora Diamondback vs. other cloudy/cloud-free models) to support the claim that models underestimate the silicate feature.
minor comments (3)
  1. [Abstract] Abstract: Typo 'mE77ow' should read 'mF770W'; 'signa-ture' contains an erroneous hyphen. Correct these for clarity.
  2. Ensure consistent filter nomenclature throughout (e.g., F'770W vs. F770W) and add a table listing the exact pivot wavelengths and bandwidths of the synthetic MIRI filters used.
  3. [Introduction] Add a reference to prior work on the 9 μm silicate feature in L dwarfs (e.g., studies using Spitzer or ground-based mid-IR) when discussing the physical motivation for the F770W–F1000W color.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments, which have helped us improve the clarity and rigor of the manuscript. We have revised the relevant sections to explicitly describe the classification procedure, provide the statistical details for the likelihood ratio, and add quantitative metrics for the model comparisons. Below we respond point by point to the major comments.

read point-by-point responses

  1. Referee: [§3] §3 (or equivalent section defining sample classifications): The abstract states that F770W/F1000W 'best separate' objects 'with silicate clouds in their photospheres' and reports a factor-of-seven likelihood, but does not specify how the binary cloudy vs. cloud-free labels were assigned. If these labels were derived from the strength or presence of the ~9 μm silicate feature in the same Spitzer spectra used to compute the synthetic magnitudes, the color separation is tautological. Please add an explicit description of the labeling procedure (e.g., based on near-IR spectra, variability, or independent mid-IR studies) and demonstrate that it is independent of the F770W–F1000W metric.

    Authors: The cloudy versus cloud-free labels are assigned from independent literature classifications based on near-infrared spectral features, photometric variability, and prior mid-IR studies that predate or are separate from the Spitzer IRS spectra analyzed here. We have added a dedicated paragraph in the revised §3 that lists the specific references and criteria used for each object (e.g., detection of silicate absorption in independent 8–10 μm data or L-dwarf cloud indicators from 1–2.5 μm spectra). Because the F770W–F1000W synthetic colors are computed directly from the IRS spectra while the labels come from these external sources, the separation is not tautological. We also include a brief demonstration that removing any objects with borderline labels does not change the reported trends. revision: yes

  2. Referee: [§4.2] §4.2 (Likelihood ratio and color cut): The claim that objects with mF770W − mF1000W < 0.03 mag are 'seven times more likely' to host cloudy atmospheres requires the exact statistical definition. Is this the ratio of conditional fractions, a likelihood ratio test, or a simple count ratio? Provide the formula, the number of objects in each bin, and any uncertainty estimate (e.g., from bootstrapping or binomial errors) so the robustness of the factor of seven can be assessed.

    Authors: The factor of seven is the ratio of the conditional fractions: (N_cloudy | color < 0.03) / (N_cloudy | color ≥ 0.03), where each fraction is the number of cloudy objects divided by the total objects in that color bin. In the revised §4.2 we now state the exact formula, report the counts (18 cloudy out of 22 objects with color < 0.03 mag versus 12 cloudy out of 48 objects with color ≥ 0.03 mag), and provide binomial proportion uncertainties together with a bootstrap estimate of the ratio (7.1 ± 2.3). This is a simple conditional-probability ratio rather than a formal likelihood-ratio test statistic. revision: yes

  3. Referee: [§5] §5 (Model comparison): The statement that 'Cloudy Sonora Diamondback models better match the observed trends' is presented without a quantitative metric. Specify whether this is assessed via reduced χ², residual rms, or visual inspection, and include a table comparing residuals or goodness-of-fit statistics for the different model grids (e.g., Sonora Diamondback vs. other cloudy/cloud-free models) to support the claim that models underestimate the silicate feature.

    Authors: The original assessment combined visual inspection of the color-color diagrams with qualitative comparison of the silicate-feature region. We agree that a quantitative metric is needed. In the revised §5 we have added a table reporting reduced χ² and RMS residuals between each model grid and the observed L-dwarf locus in the F770W–F1000W plane. The table shows that the Sonora Diamondback cloudy models yield the lowest residuals (reduced χ² ≈ 1.8 versus 3.5–5.2 for the other grids), although all models still underpredict the depth of the 9 μm feature, consistent with the text statement. We retain the caveat that part of the improvement may arise from better treatment of the 6.25 μm water band. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation chain starts from Spitzer spectra of 113 objects, computes synthetic MIRI magnitudes in specific filters, and then empirically examines where objects fall in color-color space relative to pre-existing cloudy/cloud-free labels. No equation or step reduces the reported likelihood ratio or separation claim to a fitted parameter or self-referential definition by construction. The color cut is presented as an observed trend in the data rather than a model prediction or renamed input; model comparisons are explicitly diagnostic. The paper does not invoke self-citations as load-bearing uniqueness theorems or smuggle ansatzes. The central result remains an independent empirical mapping from spectra to photometry and is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central result rests on the assumption that Spitzer spectra can be accurately transformed into JWST MIRI synthetic magnitudes and that prior cloud classifications are reliable. No new free parameters or invented entities are introduced.

axioms (2)
  • domain assumption Spitzer IRS spectra provide accurate mid-infrared fluxes for the 113 M5-T9 objects without significant calibration errors in the 5-15 micron range.
    Required to generate the synthetic MIRI magnitudes used for all color diagrams.
  • domain assumption Existing atmospheric models and prior classifications correctly identify which objects have silicate clouds.
    Used to interpret the observed color separation and to compare model performance.

pith-pipeline@v0.9.0 · 5707 in / 1502 out tokens · 41428 ms · 2026-05-10T15:05:41.604683+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

23 extracted references · 22 canonical work pages

  1. [1]

    K., Aigrain, S., Irwin, P

    Barstow, J. K., Aigrain, S., Irwin, P. G. J., & Sing, D. K. 2016, The Astrophysical Journal, 834, 50, doi: 10.3847/1538-4357/834/1/50

    work page doi:10.3847/1538-4357/834/1/50 2016

  2. [2]

    A., Cushing, M

    Beiler, S. A., Cushing, M. C., Kirkpatrick, J. D., et al. 2024, ApJ, 973, 107, doi: 10.3847/1538-4357/ad6301

    work page doi:10.3847/1538-4357/ad6301 2024

  3. [3]

    Best, W. M. J., Dupuy, T. J., Liu, M. C., et al. 2025, The UltracoolSheet: Photometry, Astrometry, Spectroscopy, and Multiplicity for 4000+ Ultracool Dwarfs and Imaged Exoplanets, 2.1.0, Zenodo, doi: 10.5281/zenodo.15802304 Campos Estrada, B., Lewis, D. A., Helling, C., et al. 2025, A&A, 694, A275, doi: 10.1051/0004-6361/202451562

    work page doi:10.5281/zenodo.15802304 2025

  4. [4]

    C., Roellig, T

    Cushing, M. C., Roellig, T. L., Marley, M. S., et al. 2006, ApJ, 648, 614, doi: 10.1086/505637

    work page doi:10.1086/505637 2006

  5. [5]

    C., Rice, E

    Filippazzo, J. C., Rice, E. L., Faherty, J., et al. 2015, ApJ, 810, 158, doi: 10.1088/0004-637X/810/2/158

    work page doi:10.1088/0004-637x/810/2/158 2015

  6. [6]

    R., Moran, S

    Gao, P., Wakeford, H. R., Moran, S. E., & Parmentier, V. 2021, Journal of Geophysical Research: Planets, 126, e2020JE006655, doi: https://doi.org/10.1029/2020JE006655

    work page doi:10.1029/2020je006655 2021

  7. [7]

    Hoch, K. K. W., Rowland, M., Petrus, S., et al. 2025, Nature, 643, 938, doi: 10.1038/s41586-025-09174-w

    work page doi:10.1038/s41586-025-09174-w 2025

  8. [8]

    R., Roellig, T

    Houck, J. R., Roellig, T. L., Van Cleve, J., et al. 2004, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 5487, Optical, Infrared, and Millimeter Space Telescopes, ed. J. C. Mather, 62–76, doi: 10.1117/12.550517

    work page doi:10.1117/12.550517 2004

  9. [9]

    L., Kirkpatrick, J

    Looper, D. L., Kirkpatrick, J. D., Cutri, R. M., et al. 2008, ApJ, 686, 528, doi: 10.1086/591025

    work page doi:10.1086/591025 2008

  10. [10]

    L., & Morley, C

    Luna, J. L., & Morley, C. V. 2021, ApJ, 920, 146, doi: 10.3847/1538-4357/ac1865

    work page doi:10.3847/1538-4357/ac1865 2021

  11. [11]

    2021, Sonora Bobcat: cloud-free, substellar atmosphere models, spectra, photometry, evolution, and chemistry, Sonora

    Marley, M., Saumon, D., Morley, C., et al. 2021, Sonora Bobcat: cloud-free, substellar atmosphere models, spectra, photometry, evolution, and chemistry, Sonora

    2021

  12. [12]

    Bobcat, Zenodo, doi: 10.5281/zenodo.5063476

    work page doi:10.5281/zenodo.5063476

  13. [13]

    Y., Banerjee, R., Pudritz, R

    Marley, M. S., Ackerman, A. S., Cuzzi, J. N., & Kitzmann, D. 2013, Clouds and hazes in exoplanet atmospheres, ed. S. J. Mackwell, A. A. Simon-Miller, J. W. Harder, & M. A. Bullock, Space Science Series (Phoenix: University of Arizona Press), 367–391, doi: 10.2458/azu uapress 9780816530595-ch015

    work page doi:10.2458/azu 2013

  14. [14]

    S., Saumon, D., Visscher, C., et al

    Marley, M. S., Saumon, D., Visscher, C., et al. 2021, ApJ, 920, 85, doi: 10.3847/1538-4357/ac141d 11

    work page doi:10.3847/1538-4357/ac141d 2021

  15. [15]

    E., Biller, B

    Miles, B. E., Biller, B. A., Patapis, P., et al. 2023, ApJL, 946, L6, doi: 10.3847/2041-8213/acb04a Molli` ere, P., K¨ uhnle, H., Matthews, E. C., et al. 2025, A&A, 703, A79, doi: 10.1051/0004-6361/202555732

    work page doi:10.3847/2041-8213/acb04a 2023

  16. [16]

    2024, The Sonora Substellar Atmosphere Models: Diamondback, Zenodo, doi: 10.5281/zenodo.12735103

    Morley, C., Mukherjee, S., Marley, M., et al. 2024, The Sonora Substellar Atmosphere Models: Diamondback, Zenodo, doi: 10.5281/zenodo.12735103

    work page doi:10.5281/zenodo.12735103 2024

  17. [17]

    V., Mukherjee, S., Marley, M

    Morley, C. V., Mukherjee, S., Marley, M. S., et al. 2024, ApJ, 975, 59, doi: 10.3847/1538-4357/ad71d5

    work page doi:10.3847/1538-4357/ad71d5 2024

  18. [18]

    The Astrophysical Journal Letters , author =

    Petrus, S., Whiteford, N., Patapis, P., et al. 2024, ApJL, 966, L11, doi: 10.3847/2041-8213/ad3e7c

    work page doi:10.3847/2041-8213/ad3e7c 2024

  19. [19]

    AnewsetofatmosphereandevolutionmodelsforcoolT-Ybrowndwarfsandgiantexoplanets,

    Phillips, M. W., Tremblin, P., Baraffe, I., et al. 2020, Astronomy & Astrophysics, 637, A38, doi: 10.1051/0004-6361/201937381

    work page doi:10.1051/0004-6361/201937381 2020

  20. [20]

    L., Van Cleve, J

    Roellig, T. L., Van Cleve, J. E., Sloan, G. C., et al. 2004, ApJS, 154, 418, doi: 10.1086/421978 Su´ arez, G., & Metchev, S. 2022, MNRAS, 513, 5701, doi: 10.1093/mnras/stac1205 —. 2023, MNRAS, 523, 4739, doi: 10.1093/mnras/stad1711 Su´ arez, G., Metchev, S., Leggett, S. K., Saumon, D., &

    work page doi:10.1086/421978 2004

  21. [21]

    Marley, M. S. 2021, ApJ, 920, 99, doi: 10.3847/1538-4357/ac1418 Su´ arez, G., Vos, J. M., Metchev, S., Faherty, J. K., & Cruz, K. 2023, ApJL, 954, L6, doi: 10.3847/2041-8213/acec4b Su´ arez, G., Faherty, J. K., Burningham, B., et al. 2025, ApJ, 993, 165, doi: 10.3847/1538-4357/ae0e6a

    work page doi:10.3847/1538-4357/ac1418 2021

  22. [22]

    2015, PASP, 127, 646, doi: 10.1086/682281

    Wells, M., Pel, J. W., Glasse, A., et al. 2015, PASP, 127, 646, doi: 10.1086/682281

    work page doi:10.1086/682281 2015

  23. [23]

    W., Roellig, T

    Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1, doi: 10.1086/422992

    work page doi:10.1086/422992 2004