{"id":"e0cf1c86-f892-48e6-9001-8dbd106a8bd1","arxiv_id":"2608.06454","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A new cloud-free grid of atmospheric and evolutionary models for ultra-cool substellar objects spanning 50 to 2400 K, with updated disequilibrium chemistry and public data.","lead":"This paper presents Flame Skimmer, a new public suite of computer models for the atmospheres and cooling histories of brown dwarfs and giant planets, reaching colder temperatures, lower gravities, and higher metallicities than earlier Sonora grids. The models are built to interpret JWST observations of the coldest substellar objects and to guide searches for true Solar System analogs.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The [M/H]=+2.0 DBMM/HBMM values inherit an acknowledged pure-H/He convective adiabat; a metallicity-dependent lapse rate could shift Table 3 enough that the headline burning limits should be re-derived.","rationale":"The reader's weakest-assumption analysis identified the same load-bearing premise: the radiative-convective equilibrium solver's pure-H/He adiabatic lapse rate applied at metallicities up to 100x solar. The manuscript itself flags this limitation in Section 4.1, so the concern is grounded in the paper's own text rather than an external consensus dispute. The central claims that the grid is valid, publicly available, and that high metallicity slows cooling are independently supported by careful benchmarking against Sonora Bobcat and Elf Owl, by the self-consistent second-quench CO2 treatment, and by the qualitative behavior of the evolutionary tracks. Those supports do not, however, retire the concern: the benchmark comparisons in Appendix A cover only [M/H] = -0.5 to +0.5 or +1.0, and therefore do not validate the thermal structure at [M/H] = +2.0 where the adiabat assumption matters most. The high-metallicity boundary conditions directly feed the evolution models that produce the headline DBMM and HBMM values, so an error in the convective lapse rate would propagate into Table 3. A concrete numerical test, recomputing the high-metallicity grid with a metallicity-dependent adiabat and regenerating the evolution tracks, would settle whether the effect is large enough to change the reported burning limits. Pending that test, the conditional verdict is appropriate: the grid should be released and used, but the extreme-metallicity burning-limit numbers should be treated as provisional.","tokens_in":37717,"tokens_out":4768,"duration_ms":46230,"concrete_test":"Recompute the equilibrium [M/H]=+1.5 and +2.0, C/O=0.458 atmospheric structures at Teff = 1200-2400 K and log(g) = 3.5-5.5 using a metallicity-dependent convective adiabat (for example, evaluating the adiabatic gradient from the same water-plus-H/He EOS used in the interior, with the appropriate mean molecular weight, instead of the pure solar H/He lapse rate), then rerun the same evolutionary code and regenerate Table 3. If the DBMM shifts by more than about 0.5 M_J or the HBMM by more than about 3 M_J at [M/H]=+2.0, the headline burning-limit values should be reported with a caveat or re-derived; if they shift by less, the current values are robust to this assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's own Section 4.1 states that PICASO's RCE solver uses an adiabatic lapse rate computed assuming a pure H/He mixture with solar abundances, even for metal-rich atmospheres, and that at the highest metallicities deviations from a non-solar H/He adiabat may alter the thermal structure and locations of convective zones. This is not a cosmetic detail: the RCE pressure-temperature profile provides the entropy-at-10-bar versus flux boundary condition for the Chachan et al. evolution models, and the DBMM (5.39 M_J) and HBMM (45.03 M_J) in Table 3 are evaluated from those tracks. At [M/H]=+2.0 the heavy-element mass fraction is large enough that both the mean molecular weight and the adiabatic gradient differ materially from the solar H/He case, and the paper reports convergence problems in exactly this regime (early L dwarfs, [M/H]=+1.5 and +2.0) with extended convective zones. Since the high-metallicity burning-limit numbers are among the most striking quantitative claims, this assumption is load-bearing. A secondary concern is the SPHINX parametric extrapolation from metallicity limits of about +0.5 to +2.0 for Teff > 2400 K initial conditions, but that chiefly affects young massive tracks rather than the 10-Gyr burning limits.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents Sonora Flame Skimmer, a new grid of cloud-free 1D atmospheric and evolutionary models for substellar objects, spanning Teff = 50–2400 K, log(g) = 2.0–5.5, [M/H] = −1.0 to +2.0, C/O = 0.229–1.14, and Kzz = 10^2–10^9 cm^2/s for disequilibrium models. The atmospheric models are computed with PICASO, including updated opacities, metallicity-dependent chemical timescales, volatile condensation in disequilibrium models, and a self-consistent treatment of the CO2 'second quench'. The evolutionary models use Chachan et al. (2025) structure equations with the 10-bar entropy matching to the atmospheric grid, updated H/He and water EOS, and core prescriptions for low masses. Headline results include slower cooling at high metallicity, negligible Kzz impact on evolution, and metallicity-dependent deuterium- and hydrogen-burning minimum masses that drop to 5.39 and 45.03 MJ at [M/H] = +2.0 (Table 3). The full grid is publicly available on Zenodo.","tokens_in":37957,"tokens_out":3499,"duration_ms":31734,"significance":"If the grid and its trends hold, this is a substantial community resource for JWST and Roman interpretation of cold brown dwarfs and directly imaged temperate planets. The explicit public data release, self-consistent treatment of volatile condensation in disequilibrium chemistry, and careful benchmarking against Sonora Bobcat and Elf Owl are clear strengths. The most striking quantitative claims are the metallicity-dependent DBMM and HBMM values, which have physically plausible mechanisms (opacity-driven cooling delay and altered interior entropy) and are qualitatively consistent with earlier work. However, the paper itself acknowledges a load-bearing simplification in the radiative-convective solver at the highest metallicities, and the quantitative impact of that simplification on the headline burning-limit values is not assessed. The SPHINX-based boundary extension also introduces unquantified extrapolation uncertainty for young massive tracks, though this is less central to the 10-Gyr burning limits.","major_comments":[{"comment":"The paper states that the PICASO RCE solver computes the convective adiabatic lapse rate assuming a pure H/He mixture with solar abundances even for models up to [M/H] = +2.0, and that deviations from a non-solar H/He adiabat may alter thermal structure and convective-zone locations. Since the 10-bar entropy boundary from these P-T profiles directly feeds the evolutionary tracks, the Table 3 values for DBMM (5.39 MJ) and HBMM (45.03 MJ) at [M/H] = +2.0 inherit this uncertainty. This is not a cosmetic issue: at 100x solar metallicity the mean molecular weight and adiabatic gradient differ materially from the solar H/He value, and the paper reports convergence difficulties precisely in this regime. I request a quantitative sensitivity test: either recompute the [M/H] = +1.5 and +2.0 P-T profiles (or a representative subset) with a metal-dependent adiabat, or provide an order-of-magnitude estimate of the resulting change in the 10-bar entropy and in the derived DBMM/HBMM. Without such a bound, the extreme-metallicity burning-limit values should be presented as provisional, and the central claim of 'strikingly low' HBMM needs qualification.","section":"Section 4.1"},{"comment":"The parametric fit to SPHINX models is calibrated over log(g) = 3–5.5 and [M/H] = −0.5 to +0.5 but is then used to extend the boundary condition to [M/H] = +2.0 and to low gravities. The paper notes that this chiefly affects young massive tracks, but the HBMM is initially approached by massive tracks that pass through Teff > 2400 K at early ages. The sensitivity of the derived HBMM at high metallicity to reasonable variations in this extrapolation (e.g., clipping instead of extrapolating) should be quantified, or the affected mass/age region should be excluded from the headline claims. As written, the uncertainty in this boundary condition is acknowledged but not propagated into the reported burning-limit values.","section":"Appendix B"},{"comment":"The DBMM comparison to Sonora Diamondback at [M/H] = +0.5 is described as 'increases by 0.12 MJ relative to Sonora Diamondback', yet the surrounding text also reports a 0.05 MJ offset at solar metallicity. Please clarify whether these offsets are computed at the same DBMM definition (50% vs 90% deuterium depletion) and the same age criterion, and whether the Diamondback comparison uses the same core-mass and EOS assumptions. A mismatch in these definitions could make the differences appear smaller or larger than physically meaningful.","section":"Section 3.3, Table 3"}],"minor_comments":[{"comment":"The phrase '10x sub-solar to 100x super-solar' is ambiguous; use '0.1x solar to 100x solar' or '10x sub-solar' is unclear.","section":"Abstract"},{"comment":"The sentence 'deviations from a non-solar H/He adiabatic lapse rate' is confusing; likely intended is 'deviations from a solar H/He adiabat' or 'deviations arising from non-solar H/He composition'. Please reword.","section":"Section 4.1"},{"comment":"The text says 'we adopt the same constant Kzz values of 10^2, 10^4, 10^7, 10^8, and 10^9', but the list is missing 10^6; Table 1 lists log(Kzz) = 2, 4, 7, 8, 9. This is consistent with Elf Owl but should be stated exactly once without the incomplete enumeration.","section":"Section 2.1.1"},{"comment":"The panel labels 'Flame Skimmer (Quench + Rainout)' and 'Elf Owl (Quench + Constant)' are clear, but the dashed lines in the left panel are condensation curves while the right panel uses dashed lines for Elf Owl; please ensure the legend distinguishes the two uses.","section":"Figure 3"},{"comment":"The text says 'the equilibrium model having a deeper convective zone compared to the disequilibrium model' for the 1300 K case; Figure 4 appears to show this for the middle panel, but the left panel is not labeled with a marker. Please specify which panel is being discussed.","section":"Section 3.1"},{"comment":"In the appendix text, the phrase 'reproduces the Sonora Elf Owl P–T profiles' is followed by 'differences emerge primarily at the lowest temperatures' - please make the degree of agreement explicit in the caption of Figure 21 as well, since the shaded regions suggest small but nonzero deviations at several temperatures.","section":"Appendix A"},{"comment":"The manuscript cites several 'submitted' and 'in prep' works (e.g., Miles et al., Kothari et al., Strampelli et al.) in the introduction and Section 4.3. While acceptable for context, the central claims should not rely on unpublished results; the DBMM/HBMM comparisons use published grids, which is good.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is clearly within scope for ApJ and represents a useful extension of the Sonora model family. The central grid and qualitative trends are defensible, and the public data release is a strong point. The main risk is that the most attention-grabbing numbers (DBMM/HBMM at [M/H]=+2.0) rest on an acknowledged simplification in the convective lapse rate, and the paper does not quantify the sensitivity. I would not reject over this if the authors can add a sensitivity estimate or soften the headline wording. The heavy self-citation pattern is noticeable but does not appear to distort the scientific content. No concerns about novelty or overlap with other groups' unpublished work."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"James, here's my read on 2608.06454. The paper delivers what it says: a cloud-free atmospheric and evolutionary grid from 50 to 2400 K, log(g) 2-5.5, metallicities -1 to +2, C/O up to 2.5x solar, with equilibrium and disequilibrium chemistry, and all of it public on Zenodo. That alone makes it a standard resource for JWST work on Y dwarfs and temperate giant planets. The genuinely new pieces are the self-consistent second-quench treatment for CO2, volatile condensation in the disequilibrium models, and the first Sonora evolution tracks with disequilibrium chemistry. The benchmarks against Bobcat and Elf Owl are honest and give me confidence that the framework changes don't break continuity with the earlier grids.\n\nThe soft spots are real but narrower than the headline hype. Section 4.1 states plainly that the RCE solver's convective lapse rate assumes a pure H/He mixture at solar abundances even for [M/H]=+2.0. The paper flags this itself, but then Table 3 quotes DBMM=5.39 and HBMM=45.03 MJ without an error bar. At 100x solar, the metal mass fraction is high enough that both the mean molecular weight and the adiabatic gradient should differ non-negligibly. The paper also reports convergence trouble in exactly that regime (early L, [M/H]=+1.5/+2.0). I would not call the central argument broken; the qualitative result that metallicity lowers the burning limits is almost certainly robust. But the specific numbers in Table 3 should be treated as provisional until someone runs cases with a metallicity-dependent adiabat or at least quantifies the shift. The SPHINX six-parameter fit for Teff>2400 is a secondary concern; it affects young massive tracks, not the 10-Gyr burning limits that carry the interesting claim. The FeH opacity fudge factor is empirical but used consistently with prior Sonora work and doesn't contaminate the evolution claims.\n\nThe citation pattern is heavy on the Sonora/PICASO family, but that's the lineage, not self-promotion; the benchmarks to external grids are there. Nobody is fitting the DBMM/HBMM to a target. This is a serious, reproducible piece of work. I'd send it to referees and push for one additional figure or table: how much the high-metallicity adiabat choice shifts the burning limits. Then accept.","headline":"A solid, community-ready grid that fills a real gap for JWST-era cold substellar work; the high-metallicity burning-limit numbers should carry a quantified caveat about the solar-H/He convective adiabat before being quoted as headline results.","tokens_in":38612,"tokens_out":2183,"would_cite":true,"duration_ms":18776,"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":"Metal-rich atmospheres lower the mass needed for deuterium burning to about 5.4 Jupiter masses, the paper argues.","keywords":["brown dwarfs","Y dwarfs","substellar atmosphere models","disequilibrium chemistry","deuterium burning minimum mass","metallicity","evolutionary models","JWST"],"falsifier":"Recompute the [M/H] = +2.0 grid with a convective lapse rate that includes the heavy-element composition (a non-solar adiabat), and check whether the deuterium-burning minimum mass remains close to 5.39 $M_{\\rm J}$; if it moves by more than about one Jupiter mass, the headline number is not robust. A complementary test is to measure CO$_2$ abundances in a metal-rich T/Y dwarf and see whether the second-quench prediction of one to two orders of magnitude enhancement holds.","tokens_in":37468,"feed_emoji":"🔥","tokens_out":5262,"duration_ms":40112,"temperature":0.7,"pith_summary":"The paper presents Sonora Flame Skimmer, a new grid of cloud-free atmospheric and evolutionary models for the coldest substellar objects, spanning effective temperatures from 50 to 2400 K, surface gravities down to $\\log(g)=2$, metallicities from 0.1x to 100x solar, and C/O ratios up to 2.5 times solar. The grid is the first in the Sonora family to pair disequilibrium chemistry with self-consistent evolutionary tracks, and it updates the treatment of volatile condensation and CO$_2$ quenching. Its central result is that atmospheric metallicity strongly shapes cooling: metal-rich objects cool more slowly, and at 100x solar metallicity the minimum masses for deuterium and hydrogen burning drop to 5.39 and 45.03 $M_{\\rm J}$. These models give observers a tool for interpreting JWST spectra of Y dwarfs and temperate giant planets, including the emerging population of directly imaged Neptune- and Saturn-like objects.","feed_headline":"Metal-rich worlds can burn deuterium at 5.4 Jupiter masses","feed_subtitle":"New Sonora models span 50-2400 K and show metallicity rewrites substellar burning limits.","key_machinery":"The load-bearing machinery is the coupled atmosphere-interior calculation: the PICASO solver iterates to radiative-convective equilibrium with quench chemistry, then the evolution code matches the atmospheric entropy at 10 bar to a homogeneous adiabatic interior, so the atmosphere's opacity and chemistry set the cooling rate. Within the atmosphere, the new elements are the metallicity-dependent quench timescales from Zahnle and Marley, the CO$_2$ second-quench treatment, the removal of PH$_3$, and the rainout condensation of H$_2$O, NH$_3$, and CH$_4$ above their saturation pressures.","core_discovery":"The central claim is that the Sonora Flame Skimmer grid, computed with the PICASO radiative-convective equilibrium code and matched at 10 bar to interior evolution models, is a valid and public description of substellar atmospheres from 2400 K down to 50 K. The paper argues that high-metallicity atmospheres slow cooling by raising opacity, that vertical mixing strength ($K_{\\rm zz}$) barely affects evolutionary tracks, and that the extreme 100x solar case lowers the deuterium-burning minimum mass to 5.39 $M_{\\rm J}$ and the hydrogen-burning minimum mass to 45.03 $M_{\\rm J}$, roughly half the sub-solar value. It also introduces a physically motivated 'second quench' for CO$_2$, letting CO quench deep and then convert to CO$_2$ before CO$_2$ itself freezes out, raising CO$_2$ abundances by one to two orders of magnitude relative to earlier grids.","pith_inferences":["If high metallicity truly lowers the deuterium-burning minimum mass this much, surveys for old, metal-rich 'planetary-mass' objects may find deuterium-depleted brown dwarfs with masses around 5-6 $M_{\\rm J}$, blurring the usual planet/brown-dwarf division.","The paper's convective lapse rate assumes a solar H/He mixture even at 100x metallicity; testing with a metal-dependent adiabat could shift the quoted 5.39 $M_{\\rm J}$ value, so the number should be read as a model-dependent estimate rather than a fixed threshold.","The same grid could be used to search for metallicity signatures in JWST color-magnitude diagrams: the paper predicts that high-metallicity objects are fainter in F444W and bluer in 4-5 micron colors due to CO$_2$, which offers a testable selection criterion.","Future cloudy versions of these cold models will likely slow cooling even further, meaning the cloud-free tracks presented here may bracket the fastest-cooling end of the substellar population."],"forward_implications":["High-metallicity substellar objects (up to 100x solar) cool slower and stay brighter and hotter at a given age than solar-metallicity objects of the same mass.","The deuterium-burning minimum mass drops from 12.21 $M_{\\rm J}$ at solar metallicity to 5.39 $M_{\\rm J}$ at 100x solar, and the hydrogen-burning minimum mass drops from 77.92 to 45.03 $M_{\\rm J}$, so composition alone can move objects across conventional brown-dwarf mass boundaries.","Disequilibrium chemistry shifts spectra significantly (especially via enhanced CO$_2$ at 4-5 microns) but has only a minor effect on radius and cooling tracks, so evolutionary ages are largely insensitive to $K_{\\rm zz}$.","The condensation treatment removes the spuriously strong H$_2$O, NH$_3$, and CH$_4$ absorption features that earlier disequilibrium grids predicted for $T_{\\rm eff}$ below roughly 300 K.","The grid extends low enough in temperature and gravity (50 K, $\\log(g)=2$) to serve as a baseline for interpreting JWST observations of directly imaged temperate planets and Y dwarfs."],"supporting_citations":[{"why":"Supplies the PICASO atmospheric model and its radiative-convective equilibrium solver that generates the grid.","marker":"Mukherjee et al. (2023)"},{"why":"Provides the metallicity-dependent chemical timescales and quench prescriptions for CO, CH4, H2O, and HCN used in the disequilibrium models.","marker":"Zahnle & Marley (2014)"},{"why":"Sonora Bobcat, the baseline grid that Flame Skimmer reproduces and extends to colder temperatures, lower gravities, and higher metallicities.","marker":"Marley et al. (2021)"},{"why":"Sonora Elf Owl v2, which introduced the CO2 second-quench correction that Flame Skimmer makes self-consistent within the atmospheric structure.","marker":"Wogan et al. (2025)"},{"why":"Provides the evolution model framework, including core-mass choices and the entropy matching at 10 bar used for the tracks.","marker":"Chachan et al. (2025)"},{"why":"Supplies the SPHINX models used to extend the atmospheric boundary condition above 2400 K for early, massive objects.","marker":"Davis et al. (2025)"},{"why":"Supplies the updated H/He equation of state used in the interior evolution models, replacing the older SCvH EOS.","marker":"Chabrier & Debras (2021)"},{"why":"Sonora Elf Owl, the earlier disequilibrium grid whose chemistry treatment Flame Skimmer refines with volatile condensation and updated quenching.","marker":"Mukherjee et al. (2024)"}],"fun_headline_variants":["High metallicity drops deuterium-burning mass to 5.4 Jupiter masses","Sonora Flame Skimmer models span 50-2400 K with metal-rich cooling","Metallicity, not mixing, controls substellar cooling in new models","New grid shows CO2 quench boosts abundances at cold extremes","Coldest substellar models: cloud-free down to 50 K"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The weakest spot is the assumption that a convective lapse rate computed for a pure solar-abundance H/He mixture stays accurate for atmospheres with up to 100 times solar metallicity, since the paper itself notes that a non-solar adiabat could shift the thermal structure and convective-zone locations.","fun_headline_variants_meta":{"raw":{"variants":["High metallicity drops deuterium-burning mass to 5.4 Jupiter masses","Sonora Flame Skimmer models span 50-2400 K with metal-rich cooling","Metallicity, not mixing, controls substellar cooling in new models","New grid shows CO2 quench boosts abundances at cold extremes","Coldest substellar models: cloud-free down to 50 K"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000762,"raw_usage":{"total_tokens":3449,"prompt_tokens":1077,"completion_tokens":2372,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":693,"completion_tokens_details":{"reasoning_tokens":2274}},"tokens_in":693,"tokens_out":2372,"duration_ms":13859,"temperature":1.0,"reasoning_tokens":2274,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:33:48.527686+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the [M/H] = +2.0 grid with a convective lapse rate that includes the heavy-element composition (a non-solar adiabat), and check whether the deuterium-burning minimum mass remains close to 5.39 $M_{\\rm J}$; if it moves by more than about one Jupiter mass, the headline number is not robust. A complementary test is to measure CO$_2$ abundances in a metal-rich T/Y dwarf and see whether the second-quench prediction of one to two orders of magnitude enhancement holds.","supporting_citations":[],"review_version":2}