{"id":"73b5539b-d12f-4a65-9c74-dee72b235837","arxiv_id":"1908.06834","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"LHS 3844b's thermal phase curve rules out thick atmospheres above 10 bar and is best fit by a dark, bare rock surface.","lead":"A new Spitzer space telescope phase curve of the small rocky exoplanet LHS 3844b shows its dayside is very hot and its nightside very cold, with no sign of heat redistribution by an atmosphere. This is some of the strongest evidence yet that hot rocky planets around small stars can lose their atmospheres and be left as bare rock.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 10-bar atmosphere exclusion rests on Eq. 2, an unpublished analytic heat-redistribution scaling validated only against a semi-grey GCM; this model dependence is the paper's main load-bearing weakness.","rationale":"The paper presents a high-quality Spitzer phase curve with robust observational evidence for a large day-night temperature contrast and no substantial heat redistribution. The two independent analyses, the Allan deviation test, and the PLD cross-check support the reliability of the measured eclipse depth and amplitude. The bare-rock interpretation is compelling and does not depend heavily on the contested analytic scaling. However, the specific claim that thick atmospheres above 10 bar are ruled out is not purely observational; it relies on Eq. 2, an unpublished scaling that has only been validated against a semi-grey GCM with a simplified opacity treatment. If this scaling is biased, the pressure exclusion could shift, and the abstract currently omits the composition qualifier that appears in the main text. The reader's weakest_assumption identified this same scaling as the leading concern, along with the limited composition grid and the unknown mass. I focus on the scaling because it is the linchpin of the quantitative pressure bound, while the energy-balance model provides only a looser, separately parameterized constraint. This concern does not overturn the central inference that LHS 3844b most likely has no substantial atmosphere, but it does justify the reader's CONDITIONAL verdict: the pressure limit should be verified through independent derivation or full-physics GCM tests before being stated without qualification. The verdict remains CONDITIONAL because the core observational result is secure and the requested check is a verification rather than a rejection.","tokens_in":11668,"tokens_out":5516,"duration_ms":60449,"concrete_test":"Independently re-derive or obtain the derivation of Eq. 2 from the companion paper, then run the FMS GCM used in Extended Data Figure 6 with the actual O2/CO2/N2 opacity tables (or a validated band model) for surface pressures of 1, 10, and 30 bar. Compare the predicted Spitzer 4.5 micron eclipse depths to the scaling predictions shown in Figure 3. If the GCM eclipses differ from the scaling by more than the 40 ppm measurement uncertainty at 10 bar, the 3-sigma pressure limit is not robust and the abstract's unqualified 'above 10 bar' claim should be weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim that thick atmospheres above 10 bar are excluded depends on the analytic heat-redistribution scaling in Methods Eq. 2, which maps surface pressure and longwave optical depth to a dayside brightness temperature. The derivation is not given; it is cited to a submitted companion paper (Koll et al., submitted). Validation in Extended Data Figure 6 is against a semi-grey GCM using a linear tau_LW = ps/1bar parameterization, not the composition-dependent radiative transfer used to produce Figure 3. If Eq. 2 mispredicts day-night heat transport for thick, non-grey atmospheres (for example, due to spectral windows, high-altitude absorbers, or circulation regime changes), the predicted eclipse depth at a given pressure could shift, changing the 3-sigma exclusion pressure. The abstract further states 'above 10 bar' without the qualifier 'for all compositions we consider' that appears in the main text, overstating the reach of a model-dependent result. The independent energy-balance model gives a complementary photospheric pressure constraint, but it relies on assumed wind speed and surface gravity and does not directly reproduce the Figure 3 pressure exclusion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a 100-hour Spitzer IRAC 4.5 µm phase curve of the ultra-short-period terrestrial exoplanet LHS 3844b. The authors find a symmetric, large-amplitude phase variation with a secondary eclipse depth of 380±40 ppm and a peak-to-trough amplitude of 350±40 ppm, corresponding to a dayside brightness temperature of 1040±40 K and a nightside brightness temperature consistent with zero within 1σ. The phase-curve shape and amplitude are fitted by a zero-albedo bare-rock model, with a 2σ upper limit on Bond albedo of 0.2. Atmospheric modeling with an analytic heat-redistribution scaling and 1D radiative transfer is used to argue that thick atmospheres above 10 bar are excluded at 3σ for the O2/CO2/N2 compositions considered, and that thinner atmospheres are unstable to escape and stellar-wind erosion, leading to the conclusion that LHS 3844b is most likely a bare rock.","tokens_in":11893,"tokens_out":8410,"duration_ms":81734,"significance":"If the atmospheric exclusion holds, this is a landmark result: it provides the first thermal phase curve of a small, terrestrial exoplanet around an M dwarf and strong observational support for the theoretical prediction that hot rocky planets around small stars do not retain thick atmospheres. The photometric analysis is unusually careful: two independent reductions (spline mapping and pixel-level decorrelation) agree, red noise is characterized with an Allan deviation plot and a full-data test, and the MCMC treatment is described in enough detail to be reproduced. The paper ships reproducible code and publicly available data. The main caveat is that the 10-bar pressure bound rests on an unpublished analytic scaling; nevertheless, the observational core of the paper is sound.","major_comments":[{"comment":"The 3σ exclusion of surface pressures above 10 bar is computed with the analytic scaling in Eq. (2), but the derivation is not given in the manuscript; it is deferred to a submitted companion paper (Koll et al., submitted). The validation in Extended Data Figure 6 uses a semi-grey GCM with τLW = ps/1 bar, so it does not test the composition-dependent, non-grey radiative transfer used to produce Figure 3. Moreover, Eq. (2) contains both τLW and an explicit (ps/1 bar)^{2/3} factor; since τLW as defined in Eq. (3) already depends on surface pressure for a fixed composition, the separate role of the pressure factor is unexplained and the scaling appears to double-count pressure. Please provide the derivation or an extended validation for thick, non-grey atmospheres, and state the regime of applicability of the scaling.","section":"Methods, 'Model for atmospheric heat redistribution'; Eq. (2)"},{"comment":"The abstract states that thick atmospheres 'above 10 bar' are ruled out, but the main text qualifies this as 'for all compositions we consider' (O2/CO2 and N2/CO2 mixtures). The quoted bound is therefore a composition-conditional statement, and an atmosphere with a different composition or a high-altitude absorber could in principle evade it. The abstract should carry the same qualifier, or the authors should demonstrate that the 10-bar bound is robust across a broader class of compositions.","section":"Abstract and main text, Figure 3"},{"comment":"The claim that thin atmospheres are unstable to erosion is based on scaling ion escape rates calculated for Proxima Centauri b, assuming a constant stellar wind flux over the planet's lifetime and an unknown surface gravity (the planet mass is not measured). These assumptions make the '0.7–7 bar' erosion estimate a lower limit with large systematic uncertainty. Since this argument is part of the final 'most likely a bare rock' conclusion, the authors should either soften the claim or add a sensitivity analysis over the assumed wind parameters and planet mass.","section":"Methods, 'Atmospheric escape due to stellar wind'"}],"minor_comments":[{"comment":"The phrase 'nightside temperature consistent with zero kelvin' should be phrased as a 1σ upper limit (0–710 K) rather than a physical zero, especially because the temperature prior forces non-negative temperatures.","section":"Abstract and Results"},{"comment":"The statement that 'the nightside temperature was below zero for the best fit sinusoid' is unphysical; please say the best-fit sinusoid would imply a formally negative brightness temperature, which motivated the spherical-harmonics model with a non-negativity prior.","section":"Methods, sinusoid inversion"},{"comment":"The surface-composition comparison is based on a single 4.5 µm flux measurement; the discrimination among basaltic, feldspathic, granitoid, and ultramafic surfaces is therefore model-dependent and should be presented as tentative rather than definitive.","section":"Figure 2"}],"recommendation":"major_revision","confidential_remarks":"Editor: The observational result is strong and the paper will be of broad interest. The main concern is the unpublished scaling in Eq. (2); given that the companion paper is 'submitted' elsewhere, I recommend asking the authors to include the derivation in the Methods or to supply the companion manuscript for review. The abstract's unqualified 'above 10 bar' claim should also be corrected. These issues are fixable and do not undermine the photometric detection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this is the first thermal phase curve of a small (1.3 R_Earth) rocky exoplanet, and it is the cleanest direct evidence yet that hot terrestrial planets around M dwarfs can lose substantial atmospheres. The observation is high quality and the analysis is unusually careful.\n\nThe measurement itself is solid. Two independent reductions—the spline mapping and PLD—give consistent eclipse depths and amplitudes. An Allan deviation plot shows the residuals bin down with photon noise, and a full-data red-noise test gives consistent results. The observed phase variation is symmetric and large: dayside brightness temperature 1040±40 K, nightside consistent with zero. That alone implies very inefficient heat redistribution, and the bare-rock model with Bond albedo below 0.2 fits naturally. The surface-composition interpretation (basaltic) is a nice extra but secondary.\n\nNow the soft spots, in proportion. The quantitative claim that thick atmospheres above 10 bar are ruled out rests on Eq. 2, an analytic scaling whose derivation is not in the paper—it is cited to a submitted companion paper. Validation is against a semi-grey GCM with tau_LW = ps/1bar, not the composition-dependent radiative transfer used for Figure 3. So the exact pressure exclusion is model-dependent. The abstract also says “above 10 bar” without the main text’s qualifier “for all compositions we consider”—a minor overstatement. The composition grid is limited to O2/CO2/N2; a high-altitude absorber or different mean molecular weight could shift the limits. The stellar-wind erosion estimate is scaled from Proxima Centauri b and assumes an Earth-like mass, which is unknown. These are real caveats, but they don’t undermine the central inference: the large day-night contrast and lack of heat redistribution strongly disfavor any substantial atmosphere. The independent energy-balance model gives a photospheric pressure limit below 0.06 bar, consistent, though with its own assumptions.\n\nThis paper is for exoplanet atmosphere and habitability readers. It deserves serious peer review—I would send it out. Revisions should ask for the scaling derivation or explicit multi-band GCM validation, and a qualifying sentence in the abstract. The data and code are public, the reductions are reproducible, and the paper is honest about most of its limitations. I’d cite it and bring it to reading group.","headline":"Solid, careful phase-curve observation that makes a strong case LHS 3844b is a bare rock; the model-dependent 10-bar pressure limit needs a bit of qualification but the headline conclusion survives.","tokens_in":12513,"tokens_out":1716,"would_cite":true,"duration_ms":18155,"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":"Spitzer 4.5 µm phase curve of LHS 3844b is symmetric and large, ruling out thick atmospheres and pointing to a bare basaltic rock.","keywords":["thermal phase curve","secondary eclipse","terrestrial exoplanet","LHS 3844b","atmospheric heat redistribution","bare rock","atmospheric escape","Spitzer IRAC"],"falsifier":"A JWST secondary-eclipse spectrum of LHS 3844b from 5 to 12 µm is the direct test: a bare basaltic rock should appear nearly featureless or show mineral emission, while any leftover CO$_2$, H$_2$O, or high-altitude haze would imprint absorption features or change the phase-curve amplitude with wavelength. A radial-velocity mass measurement would independently test the assumed surface gravity of $16\\ \\mathrm{m\\,s^{-2}}$ on which the stellar-wind erosion argument rests.","tokens_in":11452,"feed_emoji":"🪐","tokens_out":10192,"duration_ms":96673,"temperature":0.7,"pith_summary":"This paper seeks to determine whether the small, hot terrestrial exoplanet LHS 3844b has retained an atmosphere. Using 100 hours of Spitzer 4.5 µm photometry, the authors measure a thermal phase curve with a large, symmetric amplitude: the dayside reaches $1040\\pm40$ K, the nightside is consistent with zero brightness, and the peak of emission sits at the substellar point. They argue that such a curve rules out thick atmospheres, with surface pressures above about 10 bar excluded at 3$\\sigma$ for the oxygen, carbon dioxide, and nitrogen compositions modeled, and that thinner atmospheres would be eroded by the stellar wind over the planet's lifetime. If the interpretation holds, LHS 3844b most likely has no substantial atmosphere and is a bare basaltic rock, supporting the prediction that hot terrestrial planets around small stars lose their atmospheres.","feed_headline":"LHS 3844b's phase curve rules out a thick atmosphere","feed_subtitle":"Symmetric 4.5-micron phase variation favors a dark basaltic surface with no air.","key_machinery":"The load-bearing mechanism is an analytic day-night heat-redistribution scaling (Eq. 2) that predicts the dayside brightness temperature, hence eclipse depth, from surface pressure, longwave optical thickness, and equilibrium temperature; the paper validates it against a semi-grey general circulation model. On the data side, the phase curve is fit with a first-degree spherical-harmonics temperature map, which naturally produces a symmetric curve peaked at the substellar point. A 1D radiative-transfer model converts each model atmosphere's composition and pressure into a 4.5 µm eclipse depth, and an independent energy-balance model with Bond albedo, radiative-to-advective timescale ratio, and greenhouse factor gives a consistent upper limit on photospheric pressure. The stellar-wind erosion estimate and the atmospheric-evolution model together supply the argument that thin atmospheres are unstable rather than merely unobserved.","core_discovery":"On the paper's own terms, the central discovery is that the 4.5 µm thermal phase curve of LHS 3844b is large and symmetric: secondary-eclipse depth $380\\pm40$ ppm, peak-to-trough phase amplitude $350\\pm40$ ppm, and peak brightness at longitude $-6\\pm6^\\circ$, implying a dayside brightness temperature of $1040\\pm40$ K and a nightside brightness temperature consistent with zero (0–710 K at 1$\\sigma$). These numbers match a synchronously rotating, low-albedo bare rock (Bond albedo below 0.2 at 2$\\sigma$) with a basaltic surface. For the oxygen, carbon dioxide, and nitrogen atmospheres considered, the same data exclude surface pressures above about 10 bar at 3$\\sigma$, exclude CO$_2$-dominated atmospheres down to roughly Mars-like 0.006 bar, and place the best-fitting pressures below 0.1 bar. The paper further argues, from atmospheric-evolution and stellar-wind-erosion models, that a thin atmosphere would not have survived 5 Gyr: the wind alone could remove 0.7–7 bars, and only initial water inventories above 240 Earth oceans could have left a thick atmosphere. The conclusion is that LHS 3844b most likely has no substantial atmosphere today.","pith_inferences":["Beyond the paper: the same eclipse-depth-versus-pressure scaling could be applied to other ultra-short-period rocky planets already observed by Spitzer, turning a single-object result into a survey that maps where the bare-rock regime ends.","Beyond the paper: the zero nightside temperature is an inference from a truncated spherical-harmonics fit; a thin, high-altitude haze that is optically thick at 4.5 µm but transports little heat could evade the pressure limits, and a longer-wavelength phase curve would test that loophole.","Beyond the paper: if the planet is bare rock, the dayside emission spectrum should show surface mineral features rather than molecular absorption, so a future JWST secondary-eclipse spectrum could confirm or refute the basaltic surface interpretation."],"forward_implications":["The measured phase curve excludes surface pressures above about 10 bar at 3$\\sigma$ for the oxygen, carbon dioxide, and nitrogen atmospheres modeled, and CO$_2$-dominated atmospheres are excluded down to pressures as low as Mars's 0.006 bar.","A nightside brightness temperature consistent with zero, together with the small radiative-to-advective timescale ratio ($\\tau_{\\rm rad}/\\tau_{\\rm adv} < 0.3$ at 2$\\sigma$), implies the planet's photosphere lies at pressures below about 0.06 bar if any gas is present.","Atmospheric evolution modeling indicates that any water-derived atmosphere would be lost over 5 Gyr; the stellar wind alone could strip 0.7–7 bars, so a present-day atmosphere would require continuous replenishment.","The surface emission spectrum is most consistent with a basaltic composition, comparable to lunar maria or Mercury, implying widespread extrusive volcanism if the bare-rock interpretation is correct.","The result strengthens the theoretical expectation that hot terrestrial planets around small stars do not keep substantial atmospheres and motivates phase-curve observations of cooler planets, where atmospheres may survive."],"supporting_citations":[{"why":"It supplies the analytic heat-redistribution scaling that converts surface pressure into a predicted dayside brightness temperature and eclipse depth.","marker":"[10]"},{"why":"It is the prior super-Earth phase curve whose offset peak motivated the search for atmospheric circulation and provides the comparison case.","marker":"[11]"},{"why":"It documents the Spitzer IRAC instrument and Channel 2 bandpass used to obtain the 4.5 µm photometry.","marker":"[12]"},{"why":"It is the code used to generate light curves from the spherical-harmonics temperature map fitted to the phase curve.","marker":"[15]"},{"why":"It is the TESS discovery paper that provides the orbital parameters, planet radius, and stellar properties used as fixed inputs.","marker":"[18]"},{"why":"It supplies the 1D radiative-transfer models used to compute Spitzer 4.5 µm eclipse depths for each model atmosphere and surface composition.","marker":"[21]"},{"why":"It provides the atmospheric-evolution model used to argue that any water-derived atmosphere would be lost over the planet's lifetime.","marker":"[24]"},{"why":"It supplies the semi-analytical energy-balance model used as an independent fit to the phase curve, yielding the radiative-to-advective timescale limit.","marker":"[25]"},{"why":"It provides the stellar-wind ion escape rates for Proxima Centauri b that are scaled to LHS 3844b to estimate 0.7–7 bars of wind erosion.","marker":"[43]"},{"why":"It is the atmospheric-collapse study used to argue that thin CO$_2$ atmospheres are stable against collapse at LHS 3844b's high insolation.","marker":"[5]"}],"fun_headline_variants":["LHS 3844b: airless rock with a 1040 K dayside","No thick atmosphere on LHS 3844b, phase curve shows bare rock","LHS 3844b's phase curve says no atmosphere: just bare rock","Airless exoplanet LHS 3844b: hot dayside, cold nightside","LHS 3844b: a scorching rock with no sky, phase curve confirms"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the analytic heat-redistribution scaling and the oxygen, carbon dioxide, and nitrogen atmospheric compositions used in the radiative-transfer models correctly predict how eclipse depth depends on surface pressure, and that the planet's unknown mass corresponds to a surface gravity near 16 m s$^{-2}$; if any of these fail, the quoted pressure limits do not hold.","fun_headline_variants_meta":{"raw":{"variants":["LHS 3844b: airless rock with a 1040 K dayside","No thick atmosphere on LHS 3844b, phase curve shows bare rock","LHS 3844b's phase curve says no atmosphere: just bare rock","Airless exoplanet LHS 3844b: hot dayside, cold nightside","LHS 3844b: a scorching rock with no sky, phase curve confirms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001175,"raw_usage":{"total_tokens":4921,"prompt_tokens":1076,"completion_tokens":3845,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":692,"completion_tokens_details":{"reasoning_tokens":3732}},"tokens_in":692,"tokens_out":3845,"duration_ms":29081,"temperature":1.0,"reasoning_tokens":3732,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:33:03.591152+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A JWST secondary-eclipse spectrum of LHS 3844b from 5 to 12 µm is the direct test: a bare basaltic rock should appear nearly featureless or show mineral emission, while any leftover CO$_2$, H$_2$O, or high-altitude haze would imprint absorption features or change the phase-curve amplitude with wavelength. A radial-velocity mass measurement would independently test the assumed surface gravity of $16\\ \\mathrm{m\\,s^{-2}}$ on which the stellar-wind erosion argument rests.","supporting_citations":[],"review_version":1}