{"id":"e82fdc23-7a50-4a88-af54-00193f5d64bd","arxiv_id":"1908.09847","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A grid of 84,672 disequilibrium chemistry models shows vertical mixing shifts quenching pressures and spectral signatures, while the Methane Valley of CH4 detectability remains intact.","lead":"This paper introduces ChemKM, a chemical kinetic model for exoplanet atmospheres, and runs it across a grid of 84,672 model atmospheres to map how vertical mixing changes molecules like methane, water, and CO. The authors identify where disequilibrium chemistry should be visible to JWST and argue that the 'Methane Valley' of detectable methane persists despite mixing.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Methane Valley claim rests on TP profiles fixed at thermochemical equilibrium; if disequilibrium chemistry shifts the thermal structure, the valley boundaries and transition C/O values move, so the survival claim needs a self-consistent test.","rationale":"The reader identified the static-TP assumption as the weakest point, and my independent reading agrees that it is the single most load-bearing assumption for the headline claim. The Methane Valley survival claim (Section 3.5) is quantitative: it specifies a temperature window, a C/O threshold that varies with beta, and an observational anchor (HD 102195b). All of these quantities are computed on petitCODE TP profiles that assume thermochemical equilibrium and are held fixed while the composition evolves (Appendix A.6). Because CH4 and CO dominate the opacity in exactly the 800–1500 K band, the quenched composition would feed back into the thermal structure, changing quench pressures and therefore the abundances and spectra on which the valley boundaries are drawn. Section 3.4 contains the authors' own admission that this feedback matters at C/O ~ 0.95 ('could make the inversion ... chemically and radiatively unstable', 'a further self-consistent disequilibirum chemistry calculation must be performed'). I weighed alternative concerns and found them secondary: the gCV threshold of 0.05 is asserted, but it supports the quench-level analysis of Section 3.3, not the Methane Valley claim; and the valley is explicitly demonstrated only for Kzz = 10^12 in Figure 8, although the weak-mixing limit should trivially preserve it. I also considered whether the claim is circular (the valley inherited from the equilibrium TP grid), but it is not: the kinetic abundances do change the transition C/O ratios (Figure 8), so the persistence of the valley is a nontrivial chemical result, just one that is conditional on the static-TP setup. The paper has genuine independent support: the HD 189733b benchmark against Venot et al. (2012) agrees except in the microbar regime, and the chemistry-to-spectra pipeline is explicit and reproducible in principle. None of this is fatal, so REJECT or UNVERDICTED would be too harsh; the evidence justifies the paper's claims only with the static-TP qualifier made explicit and a self-consistent check identified as the decisive next step. That is exactly the CONDITIONAL verdict the reader issued, so I recommend no change. The concrete test — iterating ChemKM and petitCODE to a converged TP-composition solution for valley-spanning models and recomputing the valley boundaries — is feasible with the existing codes and would settle whether the temperature window and C/O thresholds survive; a shift of the boundaries by more than ~100 K or ~0.1 in C/O would require revising the recommendation window, while a null result would strengthen the claim.","tokens_in":32364,"tokens_out":6262,"duration_ms":65022,"concrete_test":"Select a subset of grid models spanning the Methane Valley and its edges (Teff = 800–1500 K; C/O = 0.5–1.25; log(g) = 2–5; [Fe/H] = -1 to +2; M5, K5, and G5 hosts; Kzz = 10^6, 10^9, 10^12 cm^2 s^-1). For each model, iterate the coupling: run ChemKM to diffusion equilibrium on the petitCODE TP profile, feed the resulting disequilibrium mixing ratios back into petitCODE to recompute a new TP profile, and repeat until the TP profile and abundances converge or a limit cycle is identified. Then recompute the transition C/O ratios and valley boundaries from the converged spectra using the same spectral decomposition as Section 3.5.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the Methane Valley (800–1500 K, above a C/O threshold) survives vertical mixing. ChemKM runs on temperature–pressure profiles generated by petitCODE under radiative-convective and thermochemical equilibrium, and those profiles are held fixed while the composition evolves: Appendix A.6 states 'the current version of ChemKM only considers the TP structure statically.' Since CH4 and CO are major opacity sources in exactly the 800–1500 K band, a disequilibrium composition that depletes or enhances CH4 would change the opacity and hence the TP profile; that would shift the quench pressures, the transition C/O ratios in Figure 8, and the 800/1500 K valley boundaries, and would also affect the HD 102195b anchoring and the 1000–1800 K JWST recommendation. This is not a speculative edge case: Section 3.4 itself concludes, for the C/O ~ 0.95 minimum-IR-opacity region, that the feedback 'could make the inversion ... chemically and radiatively unstable' and that 'a further self-consistent disequilibirum chemistry calculation must be performed.' The paper therefore flags that the quantitative boundaries of its headline result are inherited from equilibrium TP structures, not computed self-consistently. The claim is not disproven: the Venot et al. (2012) benchmark, the Appendix B consistency checks, and the qualitative persistence of the valley at Kzz = 10^12 provide genuine support. But the strongest statement justified by the presented runs is that the valley survives under a static equilibrium-TP assumption, and the quantitative window and threshold should be treated as provisional until the feedback loop is tested.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents ChemKM, a 1D chemical kinetic model for irradiated exoplanet atmospheres, and applies it to a large grid of 84,672 models built from petitCODE cloud-free, radiative-convective equilibrium temperature-pressure profiles. The model is benchmarked against the Venot et al. (2012) HD 189733b photochemical model, with good agreement except in the microbar regime, where molecular diffusion and photolysis dominate. The authors introduce a geometric coefficient of variation (gCV) and assert that gCV_i = 0.05 marks the onset of disequilibrium for species i. Using this metric, they find that quenching pressures depend on effective temperature, surface gravity, metallicity, and C/O ratio. The central claim is that the 'Methane Valley' (roughly 800-1500 K, above a C/O threshold) still exists when vertical mixing is included, supported by the CH4 detection on HD 102195b. The paper further recommends JWST targets with Teff between 1000 and 1800 K around M dwarfs with low gravity, high metallicity, and C/O near unity, and presents Spitzer color maps suggesting that the two main color populations are largely insensitive to vertical mixing, with deviations attributed to clouds.","tokens_in":32701,"tokens_out":4827,"duration_ms":54529,"significance":"If correct, the persistence of the Methane Valley under vertical mixing is an important result: it validates the authors' earlier classification scheme from Paper I, and gives observers an actionable, falsifiable statement that CH4 detections and non-detections in the 800-1500 K range are diagnostic of clouds or other processes rather than diffusive disequilibrium. The paper's strengths are its unusually broad parameter coverage, the explicit external benchmark against Venot et al. (2012), and the verification tests in Appendix B showing that kinetic steady states converge to Gibbs free-energy equilibrium in the appropriate limits. The JWST target-selection recommendations and the Spitzer color-population analysis are directly useful for planning observations. However, the significance is tempered by the static temperature-pressure treatment, which is acknowledged in Appendix A.6, and by the asserted rather than sensitivity-tested gCV threshold. These issues affect the quantitative boundaries of the Methane Valley, the transition C/O values in Figure 8, and the detectability maps in Figures 5-7, even though the qualitative survival of the valley may be robust.","major_comments":[{"comment":"The manuscript's central claim that the Methane Valley survives vertical mixing rests on TP profiles that are computed in radiative-convective/thermochemical equilibrium and held fixed while the composition evolves. Appendix A.6 states this explicitly: 'the current version of ChemKM only considers the TP structure statically.' In the 800-1500 K valley, CH4 and CO are major opacity sources, so a disequilibrium abundance change can alter the thermal structure, which in turn shifts the quenching pressures, the transition C/O lines in Figure 8, and the detectability maps in Figures 5-7. Section 3.4 itself concedes that the feedback 'could make the inversion ... chemically and radiatively unstable' and calls for a self-consistent disequilibrium calculation. I ask the authors to demonstrate, on a representative subset spanning the valley boundaries, that iterating the disequilibrium abundances back into the radiative model leaves the valley's existence and approximate boundaries unchanged, or alternatively to state explicitly how large the uncertainty in the valley boundaries is. Without this, the strongest claim in the abstract and conclusions goes beyond what the presented runs justify.","section":"Appendix A.6; Sections 3.4 and 3.5"},{"comment":"The gCV metric is a useful diagnostic, but the threshold gCV_i = 0.05 for the onset of disequilibrium is asserted without sensitivity analysis. The paper's quantitative quenching-pressure results, the parameter dependencies reported in Section 3.3, and the 'constant profiles' caveat indicated by the dotted line in Figure 4 all depend on this calibration. Please show how the quenching-pressure maps change if the threshold is varied over, say, 0.01-0.1, or compare the gCV-based quenching levels with a timescale-based criterion for a representative subset. This test would establish whether the reported dependencies on Teff, log g, [Fe/H], and C/O are robust features of the model or artifacts of the chosen threshold.","section":"Section 3.3, Eq. (4)"},{"comment":"The abstract and conclusions state that deviations of observations from the equilibrium Spitzer color maps are 'likely due to the presence of clouds and not disequilibrium processes.' This goes beyond the model, which contains no clouds and samples only three discrete Kzz values. The paper's own text is more cautious, calling clouds 'a strong contender' and deferring cloud modeling to the next paper. The authors should either soften the public-facing conclusion to 'not explained by the diffusive disequilibrium processes considered here' or provide a quantitative argument that cloud opacity in the IRAC bandpasses dominates the expected disequilibrium shifts. As written, the cloud conclusion is an interpretation rather than a result of this paper's calculations.","section":"Section 3.5 and Figure 9"}],"minor_comments":[{"comment":"The notation in Eq. (4), where s_ln is defined as the sample standard deviation of log-transformed abundances and then written as s_ln = s ln(10), is confusing because s is already a standard deviation. Please define a single symbol and state clearly whether the base-10 or natural logarithm is used throughout.","section":"Section 3.3, Eq. (4)"},{"comment":"There is a typo in the phrase 'a further self-consistent disequilibirum chemistry calculation must be performed'; 'disequilibirum' should be 'disequilibrium'.","section":"Section 3.4, paragraph 6"},{"comment":"The caption uses fragments 'T op)' and 'Bottom)' rather than complete statements; also 'T op' should be 'Top'.","section":"Figure 5 caption"},{"comment":"The color maps are individually scaled, so the 20, 50, and 100 ppm contours are not directly comparable across panels. A common color scale or a statement that each panel is individually normalized would improve readability.","section":"Figure 7"},{"comment":"The HD 189733b benchmark comparison would benefit from a quantitative statement of the agreement level, not only the statement that abundances are 'almost identical' except in the microbar regime. A table or figure listing maximum differences in the overlapping pressure range would make the benchmark more reproducible.","section":"Section 3.1 and Table 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope and represents a large, careful parameter study. The authors have been transparent about the static TP limitation in Appendix A.6, which is to their credit, but that limitation is load-bearing for the headline Methane Valley claim. The requested self-consistent TP test or explicit uncertainty quantification is feasible within the scope of a revision, so I do not think rejection is warranted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper deserves a serious referee. The genuinely new pieces are the gCV quenching metric and an 84,672-model grid run at three Kzz values; the useful result is that the Methane Valley from Paper I still shows up after vertical mixing. The benchmark against Venot et al. (2012) for HD 189733b is a real external anchor, and the independent CH4 detection on HD 102195b falling inside the valley is the kind of external hit that separates this from a self-referential exercise. The appendix verification is unusually thorough for this literature.\n\nThe soft spots are real but not fatal. The biggest one is the one the paper itself flags in A.6: ChemKM keeps the TP profiles fixed at thermochemical equilibrium while the composition evolves. Since CH4 and CO are major opacity sources in exactly the 800-1500 K band, a disequilibrium composition that changes abundances could move the thermal structure, which would shift quench pressures and the transition C/O values in Figure 8. Section 3.4 even concedes that the C/O ~ 0.95 minimum-IR-opacity region could be chemically and radiatively unstable and calls for a self-consistent calculation. So the strongest statement the presented runs justify is that the valley survives under a static equilibrium-TP assumption; the quantitative window and threshold are provisional until the feedback loop is tested. I do not think this sinks the paper: the Venot benchmark, the Appendix B checks, and the persistence of the valley at Kzz=10^12 all point the same way.\n\nTwo smaller issues. The gCV = 0.05 threshold is asserted rather than derived; a short sensitivity test around that value would harden every figure that uses it. And the Spitzer color-map discussion leads to a \"deviations are likely clouds\" conclusion that goes beyond a cloud-free grid; the authors say clouds are the next paper, which is fair, but Section 3.5 reads a bit stronger than the models support. Minor: the grid is promised to be public but the cited link is not a data release, and for a paper this grid-heavy, shipping the grid and spectra matters.\n\nWho is this for? Anyone choosing JWST targets or building retrieval priors for hot and warm Jupiters, and modelers who want a quenching metric that is easy to compute. Yes, send it to peer review. Ask for: (1) a static-vs-self-consistent TP sensitivity test, or at minimum an explicit statement that the valley boundaries are preliminary; (2) gCV threshold sensitivity; (3) the actual grid and spectra release.","headline":"A solid, honest grid study whose Methane Valley result is real but whose quantitative edges are provisional until the temperature structure is allowed to respond to the chemistry.","tokens_in":33259,"tokens_out":2441,"would_cite":true,"duration_ms":24088,"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":"A grid of 84,672 chemical kinetic models finds that the Methane Valley, the 800-1500 K band where methane should dominate, survives vertical mixing, so methane detections still encode C/O and cloud information.","keywords":["exoplanet atmospheres","chemical disequilibrium","vertical mixing","methane valley","quenching pressure","transmission spectroscopy","JWST detectability","chemical kinetic model"],"falsifier":"A self-consistent calculation in which the temperature profile responds to disequilibrium abundances, or an observing campaign that finds no preferential methane detection among 800-1500 K planets with C/O above the threshold, would settle the claim. A practical version: compare JWST spectra of two matched samples inside and outside the valley, controlling for clouds; if methane occurrence does not rise across the valley boundary, the prediction fails.","tokens_in":32163,"feed_emoji":"🪐","tokens_out":10808,"duration_ms":102622,"temperature":0.7,"pith_summary":"The paper introduces a chemical kinetic model, ChemKM, and uses it to ask whether vertical mixing erases the spectral fingerprints of chemical disequilibrium in irradiated gas-giant exoplanets. Its central case is that the Methane Valley, the band between roughly 800 and 1500 K where, above a C/O threshold, methane dominates transmission spectra, survives the addition of mixing because the model finds the same valley in diffusion-equilibrium abundances. If that is right, methane detections and non-detections in that temperature band remain readable as diagnostics of C/O and cloud formation rather than being artifacts of transport. The same grid singles out a JWST sweet spot, planets near 1000-1800 K around M dwarfs with low surface gravity, high metallicity, and C/O near unity, where disequilibrium signals should be strongest.","feed_headline":"Methane Valley survives vertical mixing in exoplanet atmospheres","feed_subtitle":"An 84,672-model kinetic grid maps where JWST can see chemical disequilibrium and methane detections should cluster","key_machinery":"The engine is ChemKM, a 1D chemical kinetic model that solves the continuity-diffusion equation for more than 100 species and 1000 reactions, including eddy and molecular diffusion, photolysis, condensation, and optional influxes. The argument is carried by combining it with a quantitative quenching metric, the geometric coefficient of variation (gCV) of each species' abundance, with $gCV_i=0.05$ marking the onset of disequilibrium, and with self-consistent temperature-pressure profiles computed under radiative-convective and thermochemical equilibrium. The named object that organizes the results is the Methane Valley, the 800-1500 K band where methane becomes the dominant transmission-spectrum feature above a C/O threshold; the model's key comparison is diffusion-equilibrium spectra versus thermochemical-equilibrium spectra on a grid of more than 84,000 cases.","core_discovery":"The load-bearing discovery is the survival of the Methane Valley under vertical mixing. Starting from 28,224 self-consistent cloud-free atmospheric models and adding three mixing strengths ($K_{zz}=10^6$, $10^9$, and $10^{12}$ cm$^2$ s$^{-1}$), the authors compute 84,672 chemical kinetic models in diffusion equilibrium. They find that quenching pressure decreases with effective temperature but scatters widely with $[{\\rm Fe/H}]$, $\\log(g)$, and ${\\rm C/O}$, and that the transmission spectra of most models change most at five wavelength windows near 1, 3.3, 4.5, 12, and 15 $\\mu$m, with the 3.3 $\\mu$m CH$_4$ feature the single most sensitive tracer. Despite these changes, the region between 800 and 1500 K where methane is expected to dominate above a C/O threshold persists, and the first robust CH$_4$ detection on an irradiated planet falls inside it. In the Spitzer IRAC color diagrams the two main populations barely move with mixing, so off-population points are attributed mainly to clouds; only planets cooler than about 900 K with C/O below 0.25 show strong mixing-induced deviations.","pith_inferences":["If the Methane Valley is as robust as claimed, the same grid logic could be turned on emission spectra: planets in the valley should show correlated CH$_4$ emission and absorption behavior that separates cloud effects from transport effects more cleanly than transmission alone.","The gCV metric could be borrowed by retrieval codes as a cheap quench indicator, replacing the constant-quench-abundance assumption; doing so might change inferred $K_{zz}$ values in re-analyses of methane-deficient planets such as GJ 436b.","A coupled calculation letting chemistry alter the temperature profile is the natural next test; if it shifts quench pressures, the valley's edges and the JWST sweet spot could move by more than the current error bars.","The cold, very-low-C/O outliers singled out by the color maps are concrete follow-up targets for high-resolution spectroscopy of CO and CH$_4$ lines, since they are the only class where mixing is predicted to push a planet off the main color populations."],"forward_implications":["Transmission spectra of planets with effective temperature between 800 and 1500 K can be read as tests of C/O and cloud formation even when vertical mixing is strong, because mixing does not erase the Methane Valley.","JWST programs that target 1000-1800 K planets around M dwarfs with low surface gravity, high metallicity, and C/O near unity maximize the chance of seeing disequilibrium fingerprints.","The five spectral windows near 1, 3.3, 4.5, 12, and 15 $\\mu$m are the most promising places to look for vertical-mixing signatures in transmission.","Observed outliers in Spitzer IRAC two-color diagrams are more likely to be caused by clouds than by vertical mixing, except for very cold planets with C/O below 0.25, where mixing can matter.","Quenching pressure is not a single number per planet: it decreases with effective temperature but varies widely with metallicity, gravity, and C/O, so retrieval recipes that assume a constant quenched abundance are not generally valid."],"supporting_citations":[{"why":"Supplies the 28,224-model equilibrium grid and the four-class classification whose Methane Valley is tested against vertical mixing.","marker":"Molaverdikhani et al. (2019)"},{"why":"Computes the self-consistent radiative-convective temperature-pressure profiles and equilibrium abundances that feed ChemKM.","marker":"Mollière et al. (2015, 2017)"},{"why":"Provides the full kinetic network and the HD 189733b benchmark against which ChemKM is validated.","marker":"Venot et al. (2012)"},{"why":"Provides the thermal structure and Kzz profile used for the HD 189733b benchmark case.","marker":"Moses et al. (2011)"},{"why":"Reports the first robust CH4 detection on an irradiated planet, placing it inside the Methane Valley as predicted.","marker":"Guilluy et al. (2019)"},{"why":"Provides the comparison trend for quench pressure as a function of effective temperature at Kzz = 10^12 cm^2/s.","marker":"Venot et al. (2018)"},{"why":"Establishes the GJ 436b methane-deficiency case that motivates vertical mixing and the non-linear CH4 response to Kzz.","marker":"Madhusudhan & Seager (2011)"},{"why":"Computes the transmission and emission spectra used to build the detectability maps and color diagrams.","marker":"Mollière et al. (2019)"}],"fun_headline_variants":["Methane Valley survives vertical mixing in hot exoplanets","JWST: target 1000-1800K M-dwarf planets for methane","Spitzer colors: clouds, not disequilibrium, scatter points","84,672 kinetic models confirm methane valley stability","Chemical disequilibrium fingerprints persist in spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The temperature-pressure profile of every modelled atmosphere is fixed at its radiative-convective, thermochemical-equilibrium value, and the chemistry is not allowed to feed back and change the temperature structure; if disequilibrium chemistry alters the temperature structure, the predicted quenching pressures and detectability maps would shift.","fun_headline_variants_meta":{"raw":{"variants":["Methane Valley survives vertical mixing in hot exoplanets","JWST: target 1000-1800K M-dwarf planets for methane","Spitzer colors: clouds, not disequilibrium, scatter points","84,672 kinetic models confirm methane valley stability","Chemical disequilibrium fingerprints persist in spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000475,"raw_usage":{"total_tokens":2465,"prompt_tokens":1160,"completion_tokens":1305,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":776,"completion_tokens_details":{"reasoning_tokens":1221}},"tokens_in":776,"tokens_out":1305,"duration_ms":12892,"temperature":1.0,"reasoning_tokens":1221,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:59:48.933709+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A self-consistent calculation in which the temperature profile responds to disequilibrium abundances, or an observing campaign that finds no preferential methane detection among 800-1500 K planets with C/O above the threshold, would settle the claim. A practical version: compare JWST spectra of two matched samples inside and outside the valley, controlling for clouds; if methane occurrence does not rise across the valley boundary, the prediction fails.","supporting_citations":[{"cited_title":"I., Visscher, C., Fortney, J","cited_arxiv_id":null,"evidence_quote":"Provides the thermal structure and Kzz profile used for the HD 189733b benchmark case."},{"cited_title":"Exoplanet atmospheres with GIANO II. Detection of molecular absorption in the dayside spectrum of HD 102195b","cited_arxiv_id":"1904.04170","evidence_quote":"Reports the first robust CH4 detection on an irradiated planet, placing it inside the Methane Valley as predicted."},{"cited_title":"2018, Experimental Astronomy, 46, 101","cited_arxiv_id":null,"evidence_quote":"Provides the comparison trend for quench pressure as a function of effective temperature at Kzz = 10^12 cm^2/s."}],"review_version":1}