Panchromatic View of the Frigid Jovian Exoplanet COCONUTS-2 b
Pith reviewed 2026-05-10 18:10 UTC · model grok-4.3
The pith
Panchromatic 1-15 micron observations with ATMO2020++ models constrain COCONUTS-2 b to 496 K, 1.03 Jupiter radii, and 7.3 Jupiter masses.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The MIRI-LRS data (5.45-11 microns, accounting for 41 percent of the bolometric flux) complete the spectral energy distribution when merged with Gemini/FLAMINGOS-2, JWST/NIRSpec, WISE, and Spitzer observations. Fitting the ATMO2020++ atmospheric model grid to this combined dataset, after accounting for correlated noise via Gaussian processes, produces T_eff = 496^{+5}_{-3} K, log(g) = 4.30^{+0.04}_{-0.02} dex, [M/H] = -0.02^{+0.03}_{-0.02} dex, and R = 1.03^{+0.01}_{-0.02} R_jup. Combined with the system age of 414 ± 23 Myr, evolutionary cooling models then predict a mass of 7.3 ± 0.3 M_jup and a precise luminosity of log(L/L_⊙) = -6.166 ± 0.002 dex.
What carries the argument
Self-consistent atmospheric model grids, especially ATMO2020++, fitted to the full 1-15 micron spectral energy distribution to extract temperature, gravity, metallicity, radius, luminosity, and mass.
If this is right
- The MIRI contribution enables a bolometric luminosity precise to 0.002 dex, tightening the mass estimate from evolutionary models.
- All five tested model grids reproduce the main molecular absorption features but share a systematic flux deficit in the Y and N bands.
- The derived parameters remain consistent with earlier studies once correlated noise is modeled with Gaussian processes.
- The presence of strong NH3 absorption confirms ammonia as a detectable constituent in atmospheres near 500 K.
Where Pith is reading between the lines
- Improving the underpredicted Y- and N-band fluxes in the models would allow tighter constraints on additional opacity sources or vertical mixing.
- The same multi-instrument SED approach could be applied to other cold directly imaged companions to build a comparative sample of their atmospheric properties.
- A revised system age would scale the inferred mass linearly through the cooling tracks, altering formation-scenario interpretations.
- The clear NH3 detection at this temperature provides a calibration point for ammonia chemistry and cloud formation in future model grids.
Load-bearing premise
The self-consistent atmospheric models correctly capture the temperature structure, molecular chemistry, and radiative transfer throughout this cold atmosphere even though they systematically underpredict the observed flux in the Y and N bands.
What would settle it
An independent dynamical mass measurement for COCONUTS-2 b that lies more than 1 M_jup away from 7.3 M_jup at the 3-sigma level would falsify the combination of the derived atmospheric parameters and the cooling-model mass inference.
Figures
read the original abstract
We use a high signal-to-noise MIRI-LRS spectrum (5.45 - 11 $\mu$m, R$_\lambda$ $\sim100$) of COCONUTS-2~b revealing prominent molecular features of H$_2$O, CH$_4$ and NH$_3$. This dataset is combined with spectra from Gemini/FLAMINGOS-2 and JWST/NIRSpec (G395H), as well as photometry from WISE and Spitzer, resulting in almost continuous wavelength coverage from 1 to 15 $\mu$m. We analyze the data using five grids of self-consistent atmospheric models, spanning a wide range of T$_\text{eff}$, log(g), and [M/H]. We also investigate the use of Gaussian Processes to account for correlated noise either caused by the spectrograph or by systematic departures of models in the inversion framework. All models manage to fit the overall combined observations but predict fainter flux in Y- and N-bands. Classical model comparison suggests that the ATMO2020++ synthetic specra (with and without PH$_3$) are statistically preferred. Fitting for the correlated noise of the three spectroscopic instruments, ATMO2020++ models yields constraints consistent with previous studies and evolutionary models predictions: T$_\text{eff}$ $=496^{+5}_{-3}$ K, log(g) $=4.30^{+0.04}_{-0.02}$ dex, [M/H] $=-0.02^{+0.03}_{-0.02}$ dex, and R $=1.03^{+0.01}_{-0.02}$ R$_\text{jup}$. The extended wavelength coverage provided by MIRI (accounting for 41% of the bolometric flux) completes the SED, yielding a precise luminosity estimation of log(L/L$_{\odot}$) $=-6.166\pm0.002$ dex. Combined with a previous estimate of the system age ($414\pm23$ Myr), cooling models predict a mass of M $=7.3\pm0.3$ M$_\text{jup}$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a panchromatic SED analysis of the directly imaged exoplanet COCONUTS-2 b, combining a new high-S/N JWST/MIRI-LRS spectrum (5.45–11 μm) with NIRSpec G395H, Gemini/FLAMINGOS-2 spectra, and WISE/Spitzer photometry to achieve near-continuous 1–15 μm coverage. Five self-consistent atmospheric model grids are fitted to the data, with Gaussian Processes used to model correlated noise or systematic model–data residuals; the preferred ATMO2020++ grid yields T_eff = 496^{+5}_{-3} K, log g = 4.30^{+0.04}_{-0.02}, [M/H] = −0.02^{+0.03}_{-0.02}, R = 1.03^{+0.01}_{-0.02} R_Jup. Bolometric luminosity is obtained by direct integration of the observed SED (log L/L_⊙ = −6.166 ± 0.002), and mass (7.3 ± 0.3 M_Jup) follows from cooling tracks given the system age.
Significance. If the atmospheric-parameter posteriors are robust, the work supplies one of the most precise characterizations of a cold (~500 K) Jovian exoplanet, with MIRI supplying 41 % of the bolometric flux and completing the SED. Strengths include the use of multiple independent model grids, explicit GP treatment of correlated noise, and direct (model-independent) luminosity integration; these elements allow a clean comparison with evolutionary predictions and prior studies.
major comments (2)
- [Abstract and §3] Abstract and §3 (atmospheric modeling): All five grids, including the statistically preferred ATMO2020++ models, systematically underpredict flux in the Y- and N-bands. Because the Y-band lies on the Wien tail and the N-band overlaps the MIRI range that supplies 41 % of the bolometric flux, any residual mismatch not fully absorbed by the GP kernel can shift the joint posterior on T_eff, log g, [M/H], and R. The manuscript does not report quantitative tests (e.g., fits with those bands masked or with varied GP length-scale priors) to demonstrate that the quoted uncertainties (T_eff = 496^{+5}_{-3} K etc.) remain unbiased.
- [§4] §4 (results and fitting procedure): The exact likelihood, priors on the free parameters (T_eff, log g, [M/H], R), GP kernel forms for each of the three spectroscopic instruments, and full error budget (including hyperparameter marginalization) are not specified in sufficient detail. Without these, it is impossible to assess whether the reported parameter uncertainties fully propagate the known model–data discrepancies in the Y- and N-bands.
minor comments (2)
- [Abstract] The abstract states that MIRI accounts for 41 % of the bolometric flux but does not explicitly list the wavelength intervals used for the integration; a short table or sentence in §2 would improve clarity.
- Notation for surface gravity is given as both “log(g)” and “log g” in different places; adopt a single convention throughout.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed comments. We address each major point below and describe the revisions that will be incorporated to improve clarity and robustness.
read point-by-point responses
-
Referee: [Abstract and §3] All five grids, including the statistically preferred ATMO2020++ models, systematically underpredict flux in the Y- and N-bands. Because the Y-band lies on the Wien tail and the N-band overlaps the MIRI range that supplies 41 % of the bolometric flux, any residual mismatch not fully absorbed by the GP kernel can shift the joint posterior on T_eff, log g, [M/H], and R. The manuscript does not report quantitative tests (e.g., fits with those bands masked or with varied GP length-scale priors) to demonstrate that the quoted uncertainties (T_eff = 496^{+5}_{-3} K etc.) remain unbiased.
Authors: We explicitly note in the manuscript that all model grids underpredict flux in the Y- and N-bands. The Gaussian Process terms fitted to each spectroscopic dataset are designed to capture correlated residuals, including systematic model-data mismatches in specific wavelength regions. The derived parameters remain consistent with independent literature values and with evolutionary-model predictions at the system age, providing supporting evidence that the GP marginalization has mitigated bias. To directly address the referee's concern, we will add quantitative robustness tests in the revised manuscript: (i) refits with the Y- and N-band photometry masked, and (ii) refits with varied GP length-scale priors. These additional results will be presented to confirm that the reported posterior uncertainties are not materially affected. revision: yes
-
Referee: [§4] The exact likelihood, priors on the free parameters (T_eff, log g, [M/H], R), GP kernel forms for each of the three spectroscopic instruments, and full error budget (including hyperparameter marginalization) are not specified in sufficient detail. Without these, it is impossible to assess whether the reported parameter uncertainties fully propagate the known model–data discrepancies in the Y- and N-bands.
Authors: We agree that §4 would benefit from a more complete technical description of the fitting procedure. In the revised manuscript we will expand this section to specify: the exact likelihood function employed, the prior distributions placed on T_eff, log g, [M/H], and R, the functional form of the GP kernel (including length-scale and amplitude hyperparameters) adopted for the NIRSpec, FLAMINGOS-2, and MIRI spectra separately, and the method used to marginalize or optimize the GP hyperparameters. This added detail will make the full error budget transparent and allow readers to evaluate how model-data discrepancies are propagated into the final parameter uncertainties. revision: yes
Circularity Check
No circularity: luminosity integrated directly, parameters fitted to data, mass from independent cooling tracks
full rationale
The derivation begins with direct integration of the observed 1–15 μm SED (MIRI contributing 41% of bolometric flux) to obtain log(L/L⊙) = −6.166 ± 0.002. Atmospheric parameters are obtained by χ² minimization of five external self-consistent model grids (including ATMO2020++) against the combined spectra and photometry, with Gaussian Processes added to absorb correlated residuals in Y- and N-bands where models underpredict flux. Radius, T_eff, log g and [M/H] are therefore outputs of the fit, not inputs. Mass follows from standard evolutionary cooling models evaluated at the independently measured system age (414 ± 23 Myr). No equation defines a quantity in terms of itself, no fitted parameter is relabeled as a prediction, and no load-bearing premise rests on a self-citation chain or imported uniqueness theorem. The analysis explicitly acknowledges model–data mismatches and cross-checks results across grids, keeping the chain self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (5)
- effective temperature =
496 K
- surface gravity =
4.30 dex
- metallicity =
-0.02 dex
- radius =
1.03 Rjup
- luminosity
axioms (2)
- domain assumption Self-consistent atmospheric models (ATMO2020++ and others) accurately represent the temperature, chemistry, and radiative transfer in this cold exoplanet atmosphere
- domain assumption Evolutionary cooling models correctly map luminosity and system age to planetary mass for this object
Reference graph
Works this paper leans on
-
[1]
Ackerman , A. S. & Marley , M. S. 2001, , 556, 872
work page 2001
-
[2]
1974, IEEE Transactions on Automatic Control, 19, 716
Akaike, H. 1974, IEEE Transactions on Automatic Control, 19, 716
work page 1974
-
[3]
Alejandro Merchan , S., Faherty , J. K., Su \'a rez , G., et al. 2025, , 989, 80
work page 2025
- [4]
-
[5]
Alibert , Y., Mordasini , C., Benz , W., & Winisdoerffer , C. 2005, , 434, 343
work page 2005
-
[6]
2012, Philosophical Transactions of the Royal Society of London Series A, 370, 2765
Allard , F., Homeier , D., & Freytag , B. 2012, Philosophical Transactions of the Royal Society of London Series A, 370, 2765
work page 2012
-
[7]
Allers , K. N. & Liu , M. C. 2013, , 772, 79
work page 2013
-
[8]
2011, American Journal of Mathematical and Management Sciences, 31, 13
Ando, T. 2011, American Journal of Mathematical and Management Sciences, 31, 13
work page 2011
- [9]
-
[10]
Bailer-Jones , C. A. L., Rybizki , J., Fouesneau , M., Demleitner , M., & Andrae , R. 2021, , 161, 147
work page 2021
-
[11]
Batalha , N. E., Marley , M. S., Lewis , N. K., & Fortney , J. J. 2019, , 878, 70
work page 2019
-
[12]
A., Mukherjee , S., Cushing , M
Beiler , S. A., Mukherjee , S., Cushing , M. C., et al. 2024, , 973, 60
work page 2024
- [13]
-
[14]
2014, in Protostars and Planets VI, ed
Benz , W., Ida , S., Alibert , Y., Lin , D., & Mordasini , C. 2014, in Protostars and Planets VI, ed. H. Beuther , R. S. Klessen , C. P. Dullemond , & T. Henning , 691--713
work page 2014
-
[15]
2024, The Journal of Open Source Software, 9, 7028
Blain , D., Molli \`e re , P., & Nasedkin , E. 2024, The Journal of Open Source Software, 9, 7028
work page 2024
-
[16]
Boss , A. P. 1997, Science, 276, 1836
work page 1997
-
[17]
2016, Statistics and Computing, 26, 383
Buchner , J. 2016, Statistics and Computing, 26, 383
work page 2016
-
[18]
Cameron , A. G. W. 1978, Moon and Planets, 18, 5
work page 1978
-
[19]
Czekala , I., Andrews , S. M., Mandel , K. S., Hogg , D. W., & Green , G. M. 2015, , 812, 128
work page 2015
-
[20]
de Regt , S., Gandhi , S., Snellen , I. A. G., et al. 2024, , 688, A116
work page 2024
-
[21]
de Regt , S., Snellen , I. A. G., Allard , N. F., et al. 2025, , 696, A225
work page 2025
- [22]
-
[23]
Dupuy , T. J. & Liu , M. C. 2011, , 733, 122
work page 2011
-
[24]
K., Burningham , B., Gagn \'e , J., et al
Faherty , J. K., Burningham , B., Gagn \'e , J., et al. 2024, , 628, 511
work page 2024
- [25]
- [26]
-
[27]
Gaia Collaboration , Brown , A. G. A., Vallenari , A., et al. 2021, , 649, A1
work page 2021
-
[28]
Gaidos , E., Mann , A. W., L \'e pine , S., et al. 2014, , 443, 2561
work page 2014
-
[29]
P., Aigrain , S., Roberts , S., et al
Gibson , N. P., Aigrain , S., Roberts , S., et al. 2012, , 419, 2683
work page 2012
- [30]
- [31]
- [32]
- [33]
-
[34]
JWST TECERST , Ahrer , E.-M., Alderson , L., et al. 2023, , 614, 649
work page 2023
-
[35]
Kiman , R., Beichman , C. A., Ruiz Diaz , A., et al. 2026, , 171, 60
work page 2026
-
[36]
2025, arXiv e-prints, arXiv:2506.05392, doi: 10.48550/arXiv.2506.05392
Kipping , D. & Benneke , B. 2025, arXiv e-prints, arXiv:2506.05392
-
[37]
Kirkpatrick , J. D., Cushing , M. C., Gelino , C. R., et al. 2011, , 197, 19
work page 2011
-
[38]
Kirkpatrick , J. D., Martin , E. C., Smart , R. L., et al. 2019, , 240, 19
work page 2019
-
[39]
Kuiper , G. P. 1951, Proceedings of the National Academy of Science, 37, 1
work page 1951
- [40]
-
[41]
Landman , R., Stolker , T., Snellen , I. A. G., et al. 2024, , 682, A48
work page 2024
-
[42]
K., Tremblin , P., Phillips , M
Leggett , S. K., Tremblin , P., Phillips , M. W., et al. 2021, , 918, 11
work page 2021
-
[43]
Line , M. R., Brogi , M., Bean , J. L., et al. 2021, , 598, 580
work page 2021
- [44]
-
[45]
S., Saumon , D., Visscher , C., et al
Marley , M. S., Saumon , D., Visscher , C., et al. 2021, , 920, 85
work page 2021
-
[46]
Marocco , F., Kirkpatrick , J. D., Schneider , A. C., et al. 2024, , 967, 147
work page 2024
-
[47]
Meisner , A. M., Leggett , S. K., Logsdon , S. E., et al. 2023, , 166, 57
work page 2023
- [48]
- [49]
-
[50]
Molli \`e re , P., Wardenier , J. P., van Boekel , R., et al. 2019, , 627, A67
work page 2019
-
[51]
V., Mukherjee , S., Marley , M
Morley , C. V., Mukherjee , S., Marley , M. S., et al. 2024, , 975, 59
work page 2024
-
[52]
Morley , C. V., Skemer , A. J., Allers , K. N., et al. 2018, , 858, 97
work page 2018
-
[53]
Mukherjee , S., Batalha , N. E., Fortney , J. J., & Marley , M. S. 2023, , 942, 71
work page 2023
-
[54]
Mukherjee , S., Fortney , J. J., Batalha , N. E., et al. 2022, , 938, 107
work page 2022
-
[55]
Mukherjee , S., Fortney , J. J., Morley , C. V., et al. 2024, , 963, 73
work page 2024
- [56]
- [57]
- [58]
-
[59]
W., Tremblin , P., Baraffe , I., et al
Phillips , M. W., Tremblin , P., Baraffe , I., et al. 2020, , 637, A38
work page 2020
-
[60]
Piette , A. A. A. & Madhusudhan , N. 2020, , 497, 5136
work page 2020
-
[61]
B., Hubickyj , O., Bodenheimer , P., et al
Pollack , J. B., Hubickyj , O., Bodenheimer , P., et al. 1996, , 124, 62
work page 1996
- [62]
- [63]
- [64]
-
[65]
Ruffio , J.-B., Xuan , J. W., Chachan , Y., et al. 2026, Nature Astronomy [ [arXiv] 2601.08227 ]
-
[66]
Schneider , A. C., Shkolnik , E. L., Allers , K. N., et al. 2019, , 157, 234
work page 2019
-
[67]
Schneider , A. D. & Bitsch , B. 2021, , 654, A71
work page 2021
-
[68]
Sellke, T., Bayarri, M. J., & Berger, J. O. 2001, The American Statistician, 55, 62, publisher: [American Statistical Association, Taylor & Francis, Ltd.]
work page 2001
-
[69]
2004, in American Institute of Physics Conference Series, Vol
Skilling , J. 2004, in American Institute of Physics Conference Series, Vol. 735, Bayesian Inference and Maximum Entropy Methods in Science and Engineering: 24th International Workshop on Bayesian Inference and Maximum Entropy Methods in Science and Engineering, ed. R. Fischer , R. Preuss , & U. V. Toussaint (AIP), 395--405
work page 2004
-
[70]
Spiegel , D. S. & Burrows , A. 2012, , 745, 174
work page 2012
-
[71]
2002, Journal of the Royal Statistical Society
Spiegelhalter, D., Best, N., Carlin, B., & \ Van Der Linde\ , A. 2002, Journal of the Royal Statistical Society. Series B: Statistical Methodology, 64, 583
work page 2002
-
[72]
Su \'a rez , G., Vos , J. M., Metchev , S., Faherty , J. K., & Cruz , K. 2023, , 954, L6
work page 2023
- [73]
- [74]
- [75]
-
[76]
Torres , C. A. O., Quast , G. R., da Silva , L., et al. 2006, , 460, 695
work page 2006
- [77]
- [78]
- [79]
-
[80]
Vorobyov , E. I. 2013, , 552, A129
work page 2013
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.