{"id":"36ff8a43-749b-4496-ab67-780e3876793a","arxiv_id":"2507.05422","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A 3D atmosphere simulation of the cold Jupiter WD 1856+534 b produces temperature, chemistry, wind, and emission-spectrum predictions across six metallicity and internal-heat scenarios.","lead":"This paper runs the first 3D general circulation model of the cold Jupiter WD 1856+534 b, a giant planet orbiting a white dwarf, and computes its chemistry, winds, and synthetic emission spectra. The results give JWST observers specific molecular fingerprints that could reveal the planet's internal heat and metal content.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The cloud-free assumption is load-bearing for the H2O-based Tint and metallicity diagnostics, because the model simultaneously predicts H2O condensation in the very regimes where H2O is the proposed indicator and never quantifies the resulting gas-phase depletion.","rationale":"The reader identifies the cloud-free atmosphere assumption as the weakest load-bearing premise, and I agree. The paper's central claim is that simulated chemistry and spectra provide observational diagnostics of Tint and metallicity, with H2O featuring in both the Tint diagnostic (H2O absorption) and the high-Tint metallicity diagnostic (H2O/CO). The model simultaneously predicts H2O cloud formation in the cold cases and the warm 100x-solar case, yet includes no condensation or cloud microphysics. The resulting gas-phase H2O depletion is asserted to be 'moderate' or 'minor' without quantification, so the magnitude of the effect on the proposed spectral indicators is unknown. Radiatively active clouds would also alter the temperature structure and circulation, affecting quench levels and all abundance ratios. This is a direct attack on the headline claim, not a peripheral numerical issue. The run-length concern (2100 days versus the 50,000-day circulation changes the authors themselves cite) is real but secondary: it would shift transport and quench pressures, whereas cloud formation removes the very molecules used for diagnosis. I keep the verdict at CONDITIONAL because the paper is a legitimate first 3D modeling application with an unusually complete data deposit, but the H2O-based diagnostics should not be used for JWST interpretation until cloudy simulations and convergence checks are performed.","tokens_in":22137,"tokens_out":5146,"duration_ms":62642,"concrete_test":"Run the Tint = 500 K, [M/H] = 2 case with a cloud microphysics scheme (e.g., CARMA) coupled to Exo-FMS, or minimally post-process the GCM output by forcing gas-phase H2O to saturation above the diagnosed condensation level with rainout. Recompute the gCMCRT emission spectrum and the H2O/CO and N2/NH3 abundance ratios. If the column-integrated gas-phase H2O above the cloud base drops by more than a factor of 2 relative to the cloud-free run, or if the ordering of H2O/CO with metallicity reverses, the Sec. 4 diagnostics are not robust to the cloud-free assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Sec. 4: CH4/H2O/CO/CO2 absorption and abundance ratios as indicators of Tint and metallicity) requires that simulated gas-phase H2O, CH4, and NH3 abundances are representative of the planet. The model, however, is cloud-free while predicting its own condensates: Sec. 3 finds temperatures below the H2O dew point in all Tint = 100 K cases and in the warm 100x-solar case, and Sec. 6 concludes that H2O is 'likely removed by cloud formation and rain out' in these regimes. Because condensation and rainout are not simulated, the gas-phase H2O column in exactly the cases used to define the H2O/CO metallicity indicator (high Tint, 100x solar) and the H2O-absorption part of the Tint diagnostic is unmodelled. The authors state in Sec. 5 that they 'did not simulate haze, clouds, and condensation processes, which can have a large impact on the abundances of certain elements and chemical species,' but then treat the depletion as a minor correction ('moderately', 'minor depletions'). No calculation supports that magnitude. Radiatively active clouds would additionally change T(p), circulation, and quench levels for CO/CO2/CH4/NH3, feeding back on every proposed abundance ratio. This is not a peripheral caveat; it is the dominant unquantified systematic in the paper's observational diagnostics.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents the first 3D general circulation model (GCM) study of the cold Jupiter WD 1856+534 b, a planet orbiting a white dwarf. The authors use Exo-FMS with correlated-k radiative transfer, mixing-length theory, and the mini-chem kinetics scheme, running a 2x2 grid in internal temperature (100 K and 500 K) and metallicity (1x, 10x, 100x solar). They post-process the GCM output with gCMCRT to produce synthetic emission spectra and brightness temperature curves. The central claim is that simulated absorption features of CH4, H2O, CO, and CO2, as well as abundance ratios such as H2O/CO and N2/NH3, can serve as observational diagnostics of the planet's internal temperature and metallicity, and that these diagnostics can be tested with current and future JWST observations. The paper also compares its T-p profiles with 1D models (ATMO2020, Sonora Bobcat) and with the retrieved temperature from Limbach et al. (2025), and it discusses implications for the atmospheric dynamics and chemistry of cold giant planets.","tokens_in":22430,"tokens_out":5935,"duration_ms":69678,"significance":"If the proposed diagnostics are robust, this paper provides a useful interpretative framework for JWST observations of WD 1856+534 b, currently the only detected cold Jupiter outside the Solar System. The work is reproducible: simulation outputs are archived on Zenodo, the codes are available on GitHub, and the model setup is described in detail. The internal trends (warmer core and higher metallicity favor CO/CO2; colder cases favor H2O/CH4) are thermochemically consistent. However, the central diagnostic claim rests on cloud-free gas-phase chemistry, while the model itself predicts H2O condensation in exactly the regimes where H2O is proposed as an indicator. The manuscript acknowledges this limitation but does not quantify its effect on the abundances or the spectra. The significance is therefore conditional on the sensitivity of the diagnostics to clouds and condensation, which is currently unassessed.","major_comments":[{"comment":"The cloud-free assumption is load-bearing for the proposed Tint and metallicity diagnostics. The model predicts temperatures below the H2O dew point in all Tint = 100 K cases and in the warm 100x solar case (Section 3), and Section 6 states that H2O is 'likely removed by cloud formation and rain out' in those regimes. Yet the proposed H2O-based indicators (H2O/CO at high Tint, H2O absorption features, and the H2O abundances in the cold cases) are computed from gas-phase chemistry that does not include condensation. Section 5 acknowledges that 'we did not simulate haze, clouds, and condensation processes,' but then expects 'minor depletions' without a quantitative calculation; the Jupiter-based analogy is not justified for the much higher H2O abundances predicted here. Radiatively active clouds would also change the temperature structure, circulation, and quench levels, feeding back on every proposed abundance ratio. I request either a sensitivity test (e.g., removing condensate H2O above the cloud base and recomputing spectra) or an explicit reframing of the diagnostics as valid only for a cloud-free atmosphere.","section":"Section 5 (with Figures 1/2 and Section 6)"},{"comment":"The proposed metallicity indicators are not stated consistently. Section 4 states that at high Tint the N2/CH4 abundance ratio is a good metallicity indicator and that 'there is a similar analogy with the abundance ratio of H2O to CO,' while Section 6 lists only H2O/CO as the high-Tint metallicity indicator and omits N2/CH4. The paper should present a single, unified list of diagnostics, state which ratios apply in which Tint and metallicity regimes, and provide quantitative thresholds (e.g., the ratio values that distinguish 1x, 10x, and 100x solar) so that the claimed indicators can be tested against future data.","section":"Sections 4 and 6"},{"comment":"The comparison with Limbach et al. (2025) is qualitative and posterior: the paper states that the cold cases are 'in line with' the retrieved Teff = 184 K, but no synthetic spectrum or model photometry is quantitatively compared with the observed MIRI data. Given that the stated motivation is to support the interpretation of these JWST observations, a quantitative comparison (e.g., synthetic MIRI photometry versus the observed flux ratios, or a chi-square assessment) would materially strengthen the claim that the modeled cold cases are compatible with the data.","section":"Section 5"}],"minor_comments":[{"comment":"The object name is rendered inconsistently as 'WD-1856b+534b' in the abstract and Section 1; the standard designation is WD 1856+534 b (or WD-1856b for short after first definition).","section":"Throughout"},{"comment":"The phrase 'the the dew point' should be 'the dew point' in both figure captions.","section":"Figure 1 and Figure 2 captions"},{"comment":"The sentence 'Rhis radiative feedback from more abundant heavier elements therefore explains...' contains a typo: 'Rhis' should be 'This'.","section":"Section 5"},{"comment":"The word 'presetned' in the caption should be 'presented'.","section":"Table 3 caption"},{"comment":"The sentence 'The decreasing order of chemical species are line with Rensen et al. (2023)' should read 'The decreasing order of chemical species is in line with Rensen et al. (2023)'.","section":"Section 6"},{"comment":"The phrase 'GGM simulation output' should be 'GCM simulation output'.","section":"Data Availability"},{"comment":"The sentence 'The colourbar indicate its coordinates' in the Figure 1 caption is grammatically incomplete; it should be 'the colourbar indicates the coordinates'.","section":"Section 3"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper delivers the first 3D GCM treatment of a cold Jupiter around a white dwarf, with a clean six-case grid (1x/10x/100x solar at Tint 100/500 K), synthetic emission spectra, and a fully archived data deposit. That alone justifies a serious referee. The modeling chain is mature (Exo-FMS + mini-chem + gCMCRT), and the paper is honest about its chemistry: it compares against FastChem 2, notes the WTG behavior, and flags the missing cloud physics. The central claims about CO/CO2 versus CH4/H2O as Tint indicators are thermochemically reasonable and internally consistent.\n\nThe main soft spot is the one the authors themselves concede: no clouds, haze, or condensation are simulated, yet the paper predicts H2O cloud formation precisely in the regimes where H2O is proposed as an indicator (all Tint=100 K cases, warm 100x solar). They hand-wave the depletion as 'moderate' or 'minor' based on Jupiter analogies, but Jupiter's condensation is not a substitute for a calculation here, especially at 100x solar where H2O is a few percent. Because radiatively active clouds would feed back on T(p), circulation, and quench levels, the abundance-ratio diagnostics are not yet calibrated. I would treat them as motivated predictions, not robust observational discriminants.\n\nSecond issue: the 2100-day runs are short for this regime. The authors cite a 50,000-day simulation that produced a different circulation, so convergence evidence is needed before taking the jet structure or overturning strengths seriously. This is a minor-to-moderate issue; the main spectral conclusions are less sensitive to it.\n\nThe comparison to Limbach et al. is appropriately qualitative, with no fitting or retrieval, so no circularity. Self-citation is not a problem here because the tools are open and the output data are on Zenodo. I found no internal contradiction; the paper is clearly written and careful.\n\nBottom line: this is a solid forward-model study for a genuinely new object class. It deserves peer review, with the cloudy-simulation follow-up as a required revision before the diagnostics are used interpretively. I would cite it as the reference for WD-1856b GCM predictions.","headline":"First 3D GCM of a white-dwarf cold Jupiter with clean grid and archived data, but cloudy simulations are required before its H2O-based diagnostics are calibrated.","tokens_in":23040,"tokens_out":1651,"would_cite":true,"duration_ms":19078,"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":"The thermal spectrum of a single cold Jupiter can reveal its internal heat and metal content.","keywords":["cold Jupiter","white dwarf planets","3D general circulation model","atmospheric chemistry","emission spectroscopy","internal temperature","metallicity diagnostics","exoplanet atmospheres"],"falsifier":"Measure WD-1856b's emission spectrum between about 2 and 12 µm with JWST and compare the CH4, H2O, CO, and CO2 band strengths against the six synthetic spectra; if the observed methane-to-carbon-monoxide ratio falls between the 100 K and 500 K predictions, or if water features are far weaker than every gas-phase model, the quenching/cloud-free chemistry mechanism is insufficient.","tokens_in":79,"feed_emoji":"🪐","tokens_out":6101,"duration_ms":83047,"temperature":0.7,"pith_summary":"The paper claims that emission spectroscopy of WD-1856b, the only known cold Jupiter outside the Solar System, can determine two hidden properties of the planet: its internal temperature and its atmospheric metallicity. Across six 3D general-circulation simulations (1x, 10x, and 100x solar composition at internal temperatures of 100 K and 500 K), the authors find that the absorption strengths of CH4, H2O, CO, and CO2 in the thermal spectrum track the internal temperature, while abundance ratios such as H2O/CO at high internal temperature and N2/NH3 at low internal temperature track the metallicity. If this is right, a single JWST emission spectrum can choose among the six modeled states and turn this one object into a benchmark for understanding how Jupiter-mass planets cool and acquire heavy elements.","feed_headline":"Cold Jupiter's spectrum may expose its heat and metal content","feed_subtitle":"3D atmosphere models show methane, water, CO, and CO2 bands reveal the planet's interior state.","key_machinery":"The argument is carried by a 3D general-circulation model with online gas-phase kinetics for thirteen species, coupled to correlated-k radiative transfer, with the resulting abundances post-processed into synthetic emission spectra by a 3D Monte Carlo radiative-transfer code. The load-bearing mechanism is chemical quenching: vertical mixing freezes disequilibrium abundances at deep pressures and transports them into the photosphere, so the observable band strengths of CH4, H2O, CO, and CO2 encode the deep temperature and metallicity rather than the local equilibrium chemistry.","core_discovery":"The central discovery is a set of concrete spectral diagnostics: in the simulated planetary flux of WD-1856b, pronounced CH4 and H2O absorption marks a cold interior (100 K), especially at high metallicity, while strong CO and CO2 absorption marks a warm interior (500 K). At high internal temperature the H2O-to-CO abundance ratio is a metallicity indicator; at low internal temperature the N2-to-NH3 ratio plays that role. The paper further shows that the cold 1x and 10x solar cases produce a Jupiter-like hierarchy of H2O, CH4, and NH3, whereas the 100x solar cold case produces much more CH4 and H2O. These predictions are made from synthetic emission spectra generated by post-processing the 3D circulation and chemistry output with a Monte Carlo radiative-transfer code.","pith_inferences":["The same band-ratio diagnostics could in principle be extended to other cold giant exoplanets and directly imaged planets, where internal temperature is often unknown; the ratios offer a way to separate internal luminosity from irradiation.","If the observed spectrum shows strong CO and CO2 while the independently estimated internal temperature is low, that would imply hidden heat sources such as tidal dissipation or a recent common-envelope episode, which the paper's Tint range does not cover.","A testable extension is to check whether the predicted water-cloud depletion produces a characteristic weakening of the 2.7 and 6.3 µm water bands relative to the cloud-free models; a JWST phase curve could test this directly.","Because the paper's diagnostics rely on gas-phase quenching, adding photochemistry from the white dwarf's ultraviolet flux might change the N2/NH3 and HCN abundances, so the cold-metallicity indicator should be checked once photochemical models are available."],"forward_implications":["A JWST emission spectrum of WD-1856b can distinguish the six modeled states by comparing the CH4/H2O versus CO/CO2 band strengths.","If the planet is cold and near solar metallicity, its observable chemistry should resemble Jupiter's, with H2O, CH4, and NH3 as the main carriers of O, C, and N.","If the planet is cold and 100x solar, CH4 and H2O should rise to roughly one percent mixing ratios, producing much stronger water and methane features.","Water clouds are expected in the cold cases and in the warm 100x solar case, so gas-phase H2O and its spectral features should be somewhat depleted in the upper atmosphere.","The abundance ratios H2O/CO (warm interior) and N2/NH3 (cold interior) serve as metallicity diagnostics that can be read from the same spectrum."],"supporting_citations":[{"why":"Supplies the discovery, system parameters, and age used to set planetary radius, rotation rate, and Tint estimates.","marker":"Vanderburg et al. (2020)"},{"why":"Supplies the measured thermal emission and effective temperature (184 K) that the cold models are compared against.","marker":"Limbach et al. (2025)"},{"why":"Supplies the JWST NIRSpec observations and the planetary mass (7.4 Jupiter masses) adopted in the simulations.","marker":"MacDonald et al. (2021)"},{"why":"Supplies the GCM configuration and mixing-length theory scheme used for the white-dwarf-brown-dwarf simulations.","marker":"Lee et al. (2024)"},{"why":"Supplies the mini-chem kinetic scheme that computes the 13-species gas-phase chemistry online.","marker":"Tsai et al. (2022)"},{"why":"Supplies the 3D Monte Carlo radiative-transfer code used to post-process the GCM output into spectra.","marker":"Lee et al. (2022b)"},{"why":"Supplies the ATMO2020 1D radiative-convective equilibrium model compared against the GCM pressure-temperature profiles.","marker":"Phillips et al. (2020)"},{"why":"Supplies the Sonora Bobcat 1D atmosphere model used for the same comparison.","marker":"Marley et al. (2021)"},{"why":"Supplies the Jupiter chemical model with enrichment factors used for the cold-planet abundance comparison.","marker":"Rensen et al. (2023)"},{"why":"Supplies the FastChem 2 equilibrium chemistry code used for comparison with the kinetic mini-chem abundances.","marker":"Stock et al. (2022)"}],"fun_headline_variants":["Methane and water mark a cold core in this exoplanet","CO and CO2 betray a warm interior in a white-dwarf planet","Cold Jupiter's spectrum reveals its internal heat and metal","Spectral signs of a planet's interior from 3D models","White-dwarf planet's air traces its core temperature"],"cache_read_input_tokens":25088,"weakest_assumption_plain":"The predicted spectral diagnostics assume a cloud-free atmosphere, even though the models themselves predict water clouds in several cases; if those clouds form and are radiatively active, the temperature, circulation, and emitted spectrum could all change, weakening the abundance-ratio indicators.","fun_headline_variants_meta":{"raw":{"variants":["Methane and water mark a cold core in this exoplanet","CO and CO2 betray a warm interior in a white-dwarf planet","Cold Jupiter's spectrum reveals its internal heat and metal","Spectral signs of a planet's interior from 3D models","White-dwarf planet's air traces its core temperature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000252,"raw_usage":{"total_tokens":1598,"prompt_tokens":1018,"completion_tokens":580,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":634,"completion_tokens_details":{"reasoning_tokens":494}},"tokens_in":634,"tokens_out":580,"duration_ms":7125,"temperature":1.0,"reasoning_tokens":494,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T19:26:43.690799+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure WD-1856b's emission spectrum between about 2 and 12 µm with JWST and compare the CH4, H2O, CO, and CO2 band strengths against the six synthetic spectra; if the observed methane-to-carbon-monoxide ratio falls between the 100 K and 500 K predictions, or if water features are far weaker than every gas-phase model, the quenching/cloud-free chemistry mechanism is insufficient.","supporting_citations":[],"review_version":1}