{"id":"e8eb7f56-f70e-4173-9904-9eb00c26d423","arxiv_id":"2509.02120","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Airless, thin residual, and hazy-inversion atmospheres all match TRAPPIST-1 b's eclipse and phase curve data, so a single 15 µm eclipse depth is not enough to conclude a planet is airless.","lead":"TRAPPIST-1 b's bright JWST eclipse can be explained by an airless surface or by several other atmospheres, including thin CO2-poor ones and hazy CO2-rich ones. A 15 µm phase curve rules out many of these, but a single photometric point cannot decide if a rocky exoplanet has an atmosphere.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified","rationale":"The paper's central claim is that a single 15 µm eclipse depth is degenerate: several atmospheric scenarios—airless, thin N2-CO2 residual atmospheres, and CO2 atmospheres with high-altitude hazes—can reproduce the high dayside flux. The load-bearing condition is that at least one non-airless scenario matches the eclipse data and cannot be dismissed a priori. This is satisfied by the thin N2-CO2 grid (Table 1), which includes models with 1 ppm CO2 at 0.1 and 1 bar that match both eclipses and the phase curve. These models are not fine-tuned; they span a broad region of the grid and are physically plausible as residual atmospheres. I considered whether the CO2 collapse prediction (Sect. 3.1.6) is load-bearing, as the reader suggested. It is not: collapse is used to explain the evolutionary origin of the residual atmospheres and to mark some unstable grid points (e.g., N2+100 ppm CO2 at 0.01 bar), but those unstable points are already consistent with the data and, if anything, their instability converts them into the same thin residual atmospheres. The pure-CO2 and CO2-rich thin atmospheres that the collapse might 'rule out' are already excluded by the eclipse depths themselves, independent of any stability argument (Table 1). Thus, even a large error in the collapse threshold would not remove the degeneracy; it would only add more viable atmospheric scenarios. I also weighed the phase-curve selection (Sect. 2.4) and the relative χ2 without model uncertainties. While these are genuine limitations, they affect which scenarios survive the phase-curve cut and how 'most likely' is ranked, not the core non-uniqueness of a single photometric point. The paper explicitly acknowledges these limitations (Sect. 4). I therefore find no load-bearing concern that would change the reader's conditional verdict.","tokens_in":27225,"tokens_out":16184,"duration_ms":179914,"concrete_test":"Re-run the 0.1-bar pure-CO2 collapse simulation (Sect. 3.1.6, Fig. 7) in an independent GCM (e.g., ROCKE-3D) and in the Generic PCM at 2× horizontal resolution; check whether total collapse to a few Pa over ~10 yr persists. If it does not, the evolutionary narrative for residual atmospheres changes, but the single-point degeneracy conclusion remains because the thin N2-CO2 fits in Table 1 are independent of the collapse.","verdict_should_be":"UNCHANGED","load_bearing_attack":"I find no load-bearing flaw in the central claim that a single 15 µm eclipse depth is degenerate. The degeneracy is established by Table 1: thin N2-CO2 atmospheres (e.g., N2+1 ppm CO2 at 0.1–1 bar) match both eclipse depths and the phase curve with reduced χ2 ~1.6, as does the airless model. This result does not depend on the CO2 collapse prediction (Sect. 3.1.6). The collapse is invoked to explain how a thicker CO2 atmosphere might evolve into such a residual atmosphere, not to remove data-compatible models; pure-CO2 0.1 bar and N2+1% CO2 0.1 bar are already excluded by the eclipse photometry itself (Table 1: χ2_eclipse 17 and 15). If collapse were overestimated, additional thin CO2-bearing models (e.g., N2+100 ppm CO2 0.01 bar, currently flagged as collapsing) would remain viable, which would strengthen—not weaken—the caution against single-point conclusions. Acknowledged limitations (idealized hazes, one selected phase-curve analysis, relative χ2 without model uncertainties) affect the 'most likely airless' ranking but not the qualitative non-uniqueness. The central claim is robust.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript uses 1D radiative-convective (exo_k) and 3D Generic PCM simulations to explore which atmospheres on TRAPPIST-1 b can match the JWST secondary-eclipse depths at 12.8 and 15 µm (Greene et al. 2023; Ducrot et al. 2024) and the 15 µm thermal phase curve (Gillon et al. 2025). The authors identify four families that can yield a bright 15 µm dayside: thin residual N2-CO2 atmospheres, thick transparent atmospheres, thick greenhouse-efficient atmospheres with an opacity window near 12.8/15 µm, and atmospheres with a strong thermal inversion (CH4+CO2, idealized hazes, or dust/tholins). After computing synthetic eclipse spectra and phase curves and comparing them with a relative chi-square metric (Sect. 2.4), they find that the airless model and several atmospheric models remain compatible with all the available data. The paper concludes that the observations are consistent with an airless planet, with a thin CO2-poor residual atmosphere as the main alternative, and that a single 15 µm eclipse depth cannot be used to infer airlessness.","tokens_in":27511,"tokens_out":13142,"duration_ms":145895,"significance":"The central conclusion is important: it provides a quantitative demonstration that JWST emission photometry of TRAPPIST-1 b, including the phase curve, is degenerate between airless, thin-residual, and fine-tuned hazy/inverted atmospheres. This is a valuable warning for the many upcoming JWST programs targeting rocky exoplanets in emission. The paper's strengths are the broad N2-CO2 pressure/mixing-ratio grid (Table 1, Fig. 5), explicit 3D treatment of heat redistribution in the upper atmosphere, GCM simulations of CO2 cold-trap collapse (Fig. 7, Appendix A.2), an honest discussion of model limitations (Sect. 4), and concrete, falsifiable predictions for future 10 and 18 µm observations (Fig. 14) and transit spectra (Appendix D). The non-uniqueness conclusion does not depend on the collapse prediction: the thin N2+1 ppm CO2 and N2+0.1 ppm CO2 cases already match the eclipse and phase-curve points in Table 1. However, the quantitative ranking of the remaining scenarios (\"most likely airless\") is not supported by the relative chi-square metric used.","major_comments":[{"comment":"The paper states that the chi-square is only a relative metric and \"cannot provide an absolute quantification of the mismatch\" (Sect. 2.4), yet Table 1 labels all entries as reduced chi-square, including values below 1 (e.g., CO2+Hazes Model 1, chi2_joint = 0.4), and Sect. 3.3.3 calls Model 1 \"the best agreement\". Because the two haze parameters of Model 1 are tuned to the same eclipse/phase-curve data and no degrees-of-freedom or model-covariance treatment is given, these numbers cannot support the \"most likely airless\" ranking in the abstract and conclusion. The qualitative degeneracy remains, but the ranking needs to be either replaced by a proper model comparison or removed and presented as a prior-based or qualitative preference.","section":"§2.4, Eq. (2), Table 1"},{"comment":"The haze-inversion models that pass all data are the least physical: Model 1 mimics a single-scattering albedo of 0.5 by arbitrarily reducing the stellar flux by a factor 2, and Model 2 uses a haze factor two orders of magnitude below the value previously advocated (Sect. 3.3.2). Appendix C.1 shows that realistic aerosol optical properties (tholins, H2SO4, martian dust) do not produce a CO2 emission band and that the simplified-haze case is fine-tuned. Since the abstract lists \"CO2 atmospheres rich in hazes\" alongside the other families without these caveats, the paper should explicitly distinguish proof-of-concept/idealized models from physically self-consistent models in the abstract and conclusion. The caveats in Sect. 4.1 are not carried through to the paper's headline statements.","section":"§3.3.2, §4.1, Appendix C.1"},{"comment":"Models flagged with asterisks in Table 1 are time-dependent collapsing states, but they are assigned chi-square values on the same footing as stable equilibrium models. For example, N2+100 ppm CO2 at 0.01 bar has chi2_joint = 1.7 and is listed as collapsing; its observables are snapshots during an ongoing collapse, not a stationary state. The conclusion that collapse \"allowed us to rule out some atmospheres\" (Sect. 5) is also in tension with the statement that the data-compatible thin CO2-poor atmospheres are the end products of partial collapse (Sect. 3.1.6). The text should clarify whether collapse is used as a dynamical selection criterion or as a formation mechanism. This clarification matters, but it would not overturn the degeneracy claim: if the collapse were overestimated, additional thin CO2-bearing models would remain viable, strengthening the paper's main warning.","section":"§3.1.6, Table 1, Appendix A.2"}],"minor_comments":[{"comment":"The unit for the planetary radius is given as \"1.116 (M⊕)\"; it should be \"1.116 (R⊕)\".","section":"Table 2"},{"comment":"The N2+CO2 k-tables are from \"Chaverot et al. in prep.\" and the phase-curve data are from Gillon et al. (2025), \"under review\". For reproducibility, provide public access to the opacity tables and phase-curve data products (or an updated reference) before publication.","section":"Table 3 / §2.1.3"},{"comment":"The phrase \"the most likely scenario\" in the abstract is not derivable from the relative chi-square comparison alone; if retained, it should be explicitly qualified as a prior-based or qualitative assessment.","section":"§2.4"},{"comment":"The discussion of tidal heat flux (Table 4) is clear, but the text should state whether the quoted internal flux range (4e-2 to 500 W/m2) is used in any of the GCM simulations or only as a post hoc stability argument. It appears to be the latter.","section":"§3.1.6 / §3.1.7"},{"comment":"The non-convergence of the deep atmosphere in the 3D dust simulation is honestly disclosed; consider adding a convergence criterion or a statement that the upper-atmosphere emission layers are converged to within the quoted precision.","section":"Appendix C.2"}],"recommendation":"major_revision","confidential_remarks":"The central non-uniqueness result is robust and should be published, so I do not see a basis for rejection. The major revision is about framing and quantitative support for the \"most likely airless\" statement, the status of the idealized haze models, and the treatment of collapsing atmospheres in Table 1. For the editor: acceptance should be contingent on the Gillon et al. (2025) phase-curve paper (or its data products) being available, since it is a primary input to the comparison."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper is the strongest case so far that TRAPPIST-1 b's high 15-µm eclipse depth does not by itself imply airless. The authors run a systematic 1D/3D grid (Generic PCM + exo_k) and find four families that match the two eclipse points: thin residual N2-CO2, thick transparent, thick greenhouse with a window, and thermal-inversion with hazes. The 15-µm phase curve then kills off some, leaving bare rock, thin residual, and the fine-tuned hazy inversion. That conclusion is robust to the main caveats because the non-uniqueness is shown directly in Table 1 by forward models, not by a retrieval.\n\nWhat's genuinely new are the 3D processes. The CO2 collapse calculation (0.1 bar pure CO2 collapses in ~10 yr, no secondary equilibrium) is a real finding, and the sensitivity to substellar circulation and internal heat is handled honestly. The result that thick 'transparent' N2 atmospheres still emit on the nightside via N2-N2 CIA and therefore fail the phase curve is a nice 3D-only constraint. The paper also shows the 1D full-redistribution assumption is wrong for thermal inversions because low-pressure dayside has a short radiative timescale. These effects transfer to other close-in rocky targets.\n\nSoft spots: the chi-squared is explicitly relative, not absolute, so '2σ' in Table 1 isn't calibrated for model-data mismatch. Only two eclipse points anchor the dayside spectrum. The haze cases are fine-tuned (f_haze, single-scattering albedo); realistic aerosols (tholins, dust) do not reproduce the CO2 emission band, and the authors admit this. The phase curve comparison uses one extraction (analysis #1-MG) of a combined b+c phase curve; they justify it by comparing three robust parameters, but there is pipeline dependence. No data or code is released. None of these break the central degeneracy claim, but the metrics should be clearly labeled as relative.\n\nOne caveat on the collapse: the prediction that thin pure CO2 collapses is model-sensitive (nightside cold trap, turbulent mixing), but the paper's own Table 1 shows those atmospheres are already excluded by eclipse photometry, so even if collapse is overestimated, the degeneracy warning is not affected. That shows the authors chose their load-bearing arguments well.\n\nRecommendation: worth serious peer review. Send it. A careful referee should check the selection of phase-curve extraction and the haze tuning. I'd cite it in any paper discussing single-point eclipse limits.","headline":"The paper's central warning—don't infer an airless rocky planet from a single 15-µm eclipse depth—survives scrutiny, and the 3D modeling alone is worth the read.","tokens_in":28086,"tokens_out":2128,"would_cite":true,"duration_ms":24285,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"TRAPPIST-1 b is most likely airless, but 3D climate modeling shows several atmospheres — thin residual gas, hazy thermal inversions — produce the same bright 15-micron dayside, so a single eclipse point cannot settle whether an atmosphere e","keywords":["TRAPPIST-1 b","rocky exoplanet atmospheres","secondary eclipse","thermal phase curve","3D climate modeling","atmospheric collapse","thermal inversion","JWST mid-infrared photometry"],"falsifier":"Measure the 10 µm secondary eclipse of TRAPPIST-1 b. The bare-rock and hazy CO2-emission models diverge most at this wavelength, and the paper's sensitivity calculation says roughly two JWST visits would separate them; eight visits at 18 µm would also work. A 10 µm depth matching the hazy model would overturn the airless-most-likely conclusion; one matching bare rock would confirm it.","tokens_in":27139,"feed_emoji":"🪐","tokens_out":13760,"duration_ms":120497,"temperature":0.7,"pith_summary":"TRAPPIST-1 b, the innermost rocky planet of the TRAPPIST-1 system, shows a dayside so bright in JWST mid-infrared observations that the simplest reading is a bare, airless rock. The paper asks whether that reading is unique, and shows it is not. Using 1D and 3D climate simulations compared against two eclipse depths and a full 15-micron phase curve, the authors find that an airless planet is the most likely configuration, with a thin CO2-poor residual atmosphere as the leading alternative — but also that thin pure-CO2 atmospheres would freeze out onto the nightside within about ten years, and that haze-shrouded atmospheres with a hot upper layer (a thermal inversion) can mimic the bare-rock signal. The broader payoff is a warning: concluding a planet is airless from a single photometric point is hazardous, and phase-curve or multi-wavelength observations are needed to back up such verdicts.","feed_headline":"One infrared point can't prove TRAPPIST-1 b is airless","feed_subtitle":"3D climate models find four atmospheres matching the same bright dayside; the 15-micron phase curve narrows them down.","key_machinery":"The central object is the Generic Planetary Climate Model, a 3D global climate model with correlated-k radiative transfer and an active CO2 condensation scheme, post-processed by the Pytmosph3R tool into synthetic eclipse spectra and 15 µm phase curves. It carries the argument because heat redistribution is emergent, not assumed: the model self-consistently produces the nightside cold trap that collapses thin pure-CO2 atmospheres, and the inefficient high-altitude redistribution that lets hazy atmospheres keep a hot dayside and cold nightside. The 1D exo_k model is the fast screening stage that maps the composition-pressure space before the 3D stage.","core_discovery":"The paper's central claim: all JWST emission data on TRAPPIST-1 b — 12.8 and 15 µm eclipse depths plus a 15 µm phase curve — are consistent with an airless planet, the most likely scenario; a thin CO2-poor residual atmosphere is the leading alternative. Four families of atmospheres reproduce the bright dayside: thin bare-rock-like ones, thick transparent ones, reduced greenhouse atmospheres with opacity windows at the observed wavelengths, and thermal-inversion atmospheres where CO2 emits. 3D modeling adds effects 1D cannot see: pure-CO2 atmospheres at 0.1 bar or less collapse into the nightside cold trap within a decade, and the high layers probed at 15 µm redistribute heat poorly even in h","pith_inferences":["The collapse threshold is the natural target for model intercomparison: if another GCM's nightside cold trap is warmer because of different turbulent-mixing or surface parameterizations, thin CO2 atmospheres re-enter the viable set, and the 'most likely airless' conclusion weakens.","The paper's single-photometric-point warning extends to the wider JWST rocky-planet eclipse surveys: targets flagged as airless from one eclipse may need phase-curve or multi-wavelength follow-up before being used as clean reference bodies for stellar-contamination correction.","Only idealized hazes with tuned scattering properties match the data; realistic tholin, sulfate, and dust optics do not. Laboratory or microphysical constraints on high-altitude haze optical properties around cool M dwarfs could settle whether the hazy scenario is physically plausible.","A joint fit of planets b and c in the observed phase curve — noted by the authors as the cleaner but out-of-scope approach — could sharpen the extraction of b's nightside emission and further separate the surviving scenarios."],"forward_implications":["Thin pure-CO2 atmospheres (about 0.1 bar or less) are not viable on TRAPPIST-1 b: they freeze out onto the nightside within roughly a decade, so any surviving CO2-dominated atmosphere must be thick or stabilized by other heating.","The 15 µm phase curve discriminates where eclipse depth cannot: thick, heat-redistributing atmospheres (reduced N2-CH4-C2H4 and thick transparent cases) are ruled out by their flat phase curves.","After all available emission data, the viable set shrinks to an airless planet, a thin CO2-poor residual atmosphere, and a fine-tuned hazy thermal-inversion atmosphere.","The same eclipse-then-phase-curve methodology and the same degeneracy of a single high 15 µm eclipse depth apply to other temperate rocky exoplanets observed in emission with JWST.","Two JWST visits at 10 µm could, by the paper's sensitivity estimates, separate the bare-rock scenario from the hazy CO2-emission scenario."],"supporting_citations":[{"why":"the 15 µm MIRI secondary-eclipse measurement that established the bright dayside this study must reproduce","marker":"Greene et al. 2023"},{"why":"the 12.8 µm eclipse point and the original hazy-CO2 hypothesis, re-tested here in 3D","marker":"Ducrot et al. 2024"},{"why":"the 15 µm phase curve of TRAPPIST-1 b and c from which dayside flux, nightside flux, and peak offset are extracted","marker":"Gillon et al. 2025"},{"why":"the prior 1D grid of N2-CO2 atmospheres whose exclusion limits this study redraws and extends","marker":"Ih et al. 2023"},{"why":"the cold-trap mechanism predicting CO2 collapse on tidally locked planets, confirmed here with 3D simulations","marker":"Wordsworth et al. 2015"},{"why":"the CO2 condensation scheme built into the Generic PCM used to simulate collapse","marker":"Forget et al. 2013"},{"why":"exo_k, the 1D radiative-convective model used for the fast wide-range composition exploration","marker":"Leconte 2021"},{"why":"Pytmosph3R, the post-processor that turns GCM outputs into synthetic eclipse spectra and phase curves","marker":"Falco et al. 2022"},{"why":"the adopted masses, radii, and orbital parameters for TRAPPIST-1 b","marker":"Agol et al. 2021"}],"fun_headline_variants":["Airless is best, but four atmospheres fit TRAPPIST-1 b's dayside","One bright spot doesn't rule out atmospheres on TRAPPIST-1 b","TRAPPIST-1 b: thin CO2-poor atmosphere still in the running","3D models show TRAPPIST-1 b's phase curve can't pick a winner","TRAPPIST-1 b's eclipse data: many atmospheres survive, airless leads"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The argument depends on the model's prediction that pure CO2 atmospheres of about 0.1 bar or less freeze out completely onto the cold nightside within a decade; if real turbulence, internal heat, or orbital motion keeps the nightside warmer, thin CO2 atmospheres would remain compatible with the observations.","fun_headline_variants_meta":{"raw":{"variants":["Airless is best, but four atmospheres fit TRAPPIST-1 b's dayside","One bright spot doesn't rule out atmospheres on TRAPPIST-1 b","TRAPPIST-1 b: thin CO2-poor atmosphere still in the running","3D models show TRAPPIST-1 b's phase curve can't pick a winner","TRAPPIST-1 b's eclipse data: many atmospheres survive, airless leads"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000365,"raw_usage":{"total_tokens":1885,"prompt_tokens":915,"completion_tokens":970,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":659,"completion_tokens_details":{"reasoning_tokens":856}},"tokens_in":659,"tokens_out":970,"duration_ms":9155,"temperature":1.0,"reasoning_tokens":856,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T11:50:50.268413+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the 10 µm secondary eclipse of TRAPPIST-1 b. The bare-rock and hazy CO2-emission models diverge most at this wavelength, and the paper's sensitivity calculation says roughly two JWST visits would separate them; eight visits at 18 µm would also work. A 10 µm depth matching the hazy model would overturn the airless-most-likely conclusion; one matching bare rock would confirm it.","supporting_citations":[{"cited_title":"P., Bell , T","cited_arxiv_id":null,"evidence_quote":"the 15 µm MIRI secondary-eclipse measurement that established the bright dayside this study must reproduce"},{"cited_title":"2024, Nat","cited_arxiv_id":null,"evidence_quote":"the 12.8 µm eclipse point and the original hazy-CO2 hypothesis, re-tested here in 3D"},{"cited_title":"J., et al","cited_arxiv_id":null,"evidence_quote":"the 15 µm phase curve of TRAPPIST-1 b and c from which dayside flux, nightside flux, and peak offset are extracted"},{"cited_title":"D., Kerber , L., Pierrehumbert , R","cited_arxiv_id":null,"evidence_quote":"the cold-trap mechanism predicting CO2 collapse on tidally locked planets, confirmed here with 3D simulations"},{"cited_title":"2013, , 222, 81","cited_arxiv_id":null,"evidence_quote":"the CO2 condensation scheme built into the Generic PCM used to simulate collapse"},{"cited_title":"2021, , 645, A20","cited_arxiv_id":null,"evidence_quote":"exo_k, the 1D radiative-convective model used for the fast wide-range composition exploration"},{"cited_title":"2022, , 658, A41","cited_arxiv_id":null,"evidence_quote":"Pytmosph3R, the post-processor that turns GCM outputs into synthetic eclipse spectra and phase curves"},{"cited_title":"L., et al","cited_arxiv_id":null,"evidence_quote":"the adopted masses, radii, and orbital parameters for TRAPPIST-1 b"}],"review_version":1}