{"id":"5a9f48a8-d379-47be-8307-051d01eaddfa","arxiv_id":"2509.04034","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A full-atmosphere Jupiter model with sulfur chemistry predicts a 33 ppb HCN peak in the stratosphere, testable by JUICE.","lead":"A new 1D model of Jupiter's full atmosphere predicts a thin stratospheric layer of hydrogen cyanide (HCN) reaching 33 parts per billion, a claim the JUICE mission could test. The model couples deep thermochemistry with upper-atmosphere photochemistry and includes sulfur chemistry for the first time across this pressure range.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper reports two inconsistent N2 mixing ratios (abstract: 490 ppm; Section 5: 11 ppm), and the HCN peak is produced from N2 photolysis (Eq. 36), so the 33 ppb prediction is not well-defined until this factor-45 discrepancy is resolved.","rationale":"The reader's CONDITIONAL verdict is right, but the condition should prioritize the N2 inconsistency over (or alongside) the aerosol caveat. The abstract quotes a quenched N2 mixing ratio of 490 ppm; Section 5 says 11 ppm. Both are presented as model outputs, and the difference is a factor of 45. The HCN peak—the paper's headline prediction—is generated by Eq. 36, 2CH4 + N2 -> 2HCN + 3H2, initiated by N2 photolysis, so the HCN formation rate is directly proportional to the N2 abundance in the 1e-6–6.76e-8 bar region. At those levels CH4 is still plentiful (drop-off at 1e-8 bar), so the linear scaling should hold to first order. That means the headline number 33 ppb is not even internally self-consistent: if N2 is 11 ppm, the peak should be ~0.7 ppb; if 490 ppm, it should be ~1.5 ppm. Neither value makes the prediction robust. The absence of aerosol chemistry is a real caveat but is explicitly acknowledged and is qualitatively the same direction as the observed depletion; the N2 discrepancy is a quantitative internal contradiction with no way for the reader to know which value the model actually produced. The paper also has related numerical typos (e.g., '320 mbar' likely for '320 ppm', NH3 740 ppm exceeding the total N inventory), which further argues for an editorial pass. A simple sensitivity run would resolve the concern; if the 33 ppb peak collapses or shifts by orders of magnitude with the corrected N2, the central claim must be revised. Until then, CONDITIONAL is the appropriate verdict, with the N2 check as a required condition.","tokens_in":27808,"tokens_out":9532,"duration_ms":83180,"concrete_test":"Run the ARGO model with the same setup and extract the N2 mixing ratio profile at 1e-6 bar; verify whether it is 490 ppm (as in the abstract) or 11 ppm (as in Section 5). Then recompute the HCN profile with the N2 mixing ratio fixed at the alternative value (or with a bracketing range 11–490 ppm) while keeping all other inputs identical; if the HCN peak at 2.94e-7 bar changes by more than a factor of 5, the 33 ppb claim is not robust and should be re-stated with an uncertainty estimate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract states a quenched N2 mixing ratio of 490 ppm up to 10^-6 bar, while Section 5 states \"a relatively constant value of 11 ppm from the bottom of the atmosphere up to 10^-6 bar.\" These differ by a factor of ~45. This discrepancy directly controls the paper's central claim: the newly discovered HCN layer (33 ppb peak at 2.94e-7 bar) is formed via the net reaction 2CH4 + N2 -> 2HCN + 3H2 (Eq. 36), which begins with N2 photolysis. In the 1e-6–6.76e-8 bar region, CH4 is still abundant, so HCN production is effectively first-order in N2; the HCN mixing ratio should scale roughly linearly with f_N2. If the actual N2 profile is 11 ppm, the predicted HCN peak would be ~0.7 ppb, not 33 ppb; if it is 490 ppm, the peak may exceed current observational upper limits (Weisstein & Serabyn 1996; Davis et al. 1997). The paper provides no sensitivity analysis over N2 abundance and no model output to resolve which value is correct. This is more load-bearing than the acknowledged missing aerosol chemistry, because it is an internal inconsistency in the exact quantity that drives the proposed HCN formation mechanism.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper couples a thermochemical-kinetic 1D model (ARGO with the STAND network) over 1.1e3 bar to 7.4e-11 bar to describe Jupiter's deep atmosphere, troposphere, and stratosphere in a single code. New updates include NH4SH formation/destruction, Antoine equations for NH4SH and H2S condensation, and two updated reaction rates. The model is benchmarked against several regional models (Rimmer & Helling 2016; Hu 2021; Tsai et al. 2021; Visscher et al. 2010; Moses et al. 2010; Rensen et al. 2023) and against observations. It reproduces the main hydrocarbon profiles, CO quenching (4.2 ppb at ~318 bar), the NH3-NH4SH cloud near 0.1-1 bar, and a relatively constant N2 mixing ratio (11 ppm in §5) up to 1e-6 bar. The headline result is a newly identified HCN formation region at 1e-6 to 6.76e-8 bar driven by N2 photolysis (Eq. 36), with a predicted peak of 33 ppb at 2.94e-7 bar, presented as testable by JUICE/MAJIS. The abstract additionally claims a quenched N2 mixing ratio of 490 ppm, which is inconsistent with §5 and with the N inventory in Table 1.","tokens_in":28235,"tokens_out":8916,"duration_ms":76779,"significance":"If correct, the paper would provide the first full-pressure-range coupled kinetics model for Jupiter's sulfur and nitrogen chemistry and a falsifiable stratospheric HCN prediction. Strengths include use of a published, reusable code framework; comparisons to multiple state-of-the-art codes and observations; explicit incorporation of sulfur chemistry; and a clearly articulated reaction mechanism for HCN (Eq. 36). The HCN peak is a genuine model output rather than a fit to target observables, which lowers circularity concerns. However, the N2 abundance inconsistency directly affects the magnitude of this prediction, and the NH4SH thermodynamics rest on unvalidated extrapolations. These points prevent acceptance in the current form.","major_comments":[{"comment":"The abstract states a quenched N2 mixing ratio of 490 ppm up to 10^-6 bar; §5 reports 'a relatively constant value of 11 ppm from the bottom of the atmosphere up to 10^-6 bar' — a factor-45 discrepancy. Table 1 gives N mole fraction 4.68×10^-4, so 490 ppm N2 would require 9.8×10^-4 N atoms, exceeding the total N inventory. This is load-bearing because the new HCN region (Eq. 36) begins with N2 photolysis; in the 1e-6–6.76e-8 bar window HCN production is approximately first-order in f_N2. With 11 ppm the 33 ppb peak would become ~0.7 ppb. No sensitivity analysis over f_N2 is given. Please reconcile the values, verify N conservation, and provide f_N2 sensitivity for the HCN peak.","section":"Abstract vs. §5; Eq. (36)"},{"comment":"The NH4SH equilibrium constant from Lewis (1969), Eq. (11), is explicitly valid only for 180–300 K, yet it is extrapolated to the full atmospheric temperature range, and the resulting NASA polynomials are used up to 6500 K. Similarly, the NH4SH Antoine parameters (Table 4, valid 222.1–306.4 K) are extrapolated to 110 K. The claimed NH3–NH4SH cloud layer between 0.1 and 1 bar (§4, Fig. 8) depends on these saturation relations, and the paper presents sulfur chemistry as one of its main novelties. No uncertainty, error budget, or alternative treatment is offered. Because the cloud-layer claim is a stated validation, this extrapolation needs a targeted sensitivity test (e.g., varying the equilibrium constant within its stated validity or comparing with a current NH4SH thermodynamic dataset).","section":"§2.4.1, Eq. (11); Table 4"},{"comment":"The 33 ppb HCN prediction is the central new result. The paper states 'Our model does not include aerosol chemistry and would therefore not catch this effect,' referring to Cavalié et al. (2023b), who propose heterogeneous/aerosol removal as the cause of depleted HCN. Since the new production region lies in the stratosphere (1e-6 to 6.76e-8 bar), where haze is present, the predicted peak could be an upper limit. The paper should quantify the heterogeneous-loss timescale or perform a simple sensitivity test with an adopted sticking coefficient. Without that, the 'testable by observations' claim is not yet supported. In addition, the model predicts 33 ppb while the earlier uniform-mixing upper limits of Weisstein & Serabyn (1996) (0.3 ppb) and Davis et al. (1997) (0.16 ppb) are cited; a column-density check against those limits would clarify whether the prediction is already observable.","section":"§6, Fig. 14"}],"minor_comments":[{"comment":"'down to 320 mbar at 1.55 bar' — the unit appears to be a typo for ppm (mixing ratio vs pressure).","section":"§4"},{"comment":"The units of k_f are not clearly defined; including NHTOT in Eq. (23) suggests a termolecular rate coefficient, but c1 is quoted in cm3 s-1. Define bimolecular vs termolecular coefficients to avoid dimensional confusion.","section":"§2.4.1, Eqs. (22)–(23)"},{"comment":"The 5% margin is described as applied to mixing ratios directly and then said not to correspond to a 5% change in metallicity; the interpretation is confusing.","section":"§2.2, Table 1"}],"recommendation":"major_revision","confidential_remarks":"The abstract/body discrepancy for N2 is the kind of issue that should be caught in a careful proofread but is also scientifically load-bearing because the HCN peak scales with f_N2. The self-citation to ARGO/STAND is appropriate because the code is published; no novelty concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, here's my take. The paper is a serious modeling effort: one 1D kinetics code spanning 1100 bar down to 7e-11 bar, with updated NH4SH sulfur chemistry and a consistent STAND network. The authors benchmark against VULCAN, Hu, Moses, and observations, and the full-atmosphere span is genuinely useful for giant-planet and exoplanet work. The best thing in it is the specific, falsifiable prediction of a stratospheric HCN layer forming near 3e-7 bar via N2 + CH4 radical chemistry, which JUICE could actually test. They are also transparent about limitations: no aerosol chemistry, extrapolated NH4SH thermodynamics, no heterogeneous loss for HCN.\n\nBut there is a load-bearing internal inconsistency that the paper does not resolve. The abstract reports a quenched N2 mixing ratio of 490 ppm up to 1e-6 bar; Section 5 reports a constant 11 ppm from the bottom of the atmosphere up to 1e-6 bar. That is a factor of roughly 45. Since the HCN peak is manufactured from N2 photolysis (their Eq. 36), the 33 ppb peak scales essentially linearly with N2. If the actual N2 is 11 ppm, the HCN peak is around 0.7 ppb, not 33 ppb; if 490 ppm, it may exceed existing upper limits (Weisstein and Serabyn; Davis et al.). The paper offers no sensitivity analysis on N2, so the central prediction is not well-defined. This is more serious than the acknowledged missing aerosol chemistry, and it needs to be fixed before the HCN claim can be trusted.\n\nA few other soft spots in proportion: the NH4SH thermodynamics extrapolate the Lewis (1969) equilibrium constant from 180-300 K across the full temperature range, and the Antoine parameters are also extrapolated; that is a real uncertainty but not necessarily fatal since the condensation level is roughly right. There is no code archive, only data available upon request, which makes replication difficult for a model this complex. The comparison to observations is mostly qualitative, which is fine for an initial tool paper but should be stated more clearly.\n\nWho is the audience? People working on giant-planet photochemistry and JUICE planning. The paper deserves a serious referee: the model is useful, the issues are fixable, and the HCN prediction is worth resolving. I would send it to peer review but flag the N2 discrepancy as a must-address.","headline":"Useful full-atmosphere Jupiter chemistry model with a testable HCN prediction, but the abstract's N2 value contradicts the body and undermines the HCN peak.","tokens_in":28662,"tokens_out":3179,"would_cite":false,"duration_ms":26097,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A 1D kinetics model spanning Jupiter's full atmosphere predicts a new stratospheric HCN layer peaking at 33 ppb at 2.94e-7 bar.","keywords":["Jupiter","hydrogen cyanide","photochemistry","thermochemical kinetics","chemical reaction networks","planetary atmospheres","NH4SH clouds","stratosphere"],"falsifier":"A limb or occultation retrieval resolving the 1e-6 to 6.76e-8 bar window that returns an HCN upper limit below a few ppb at ~3e-7 bar (well under the predicted 33 ppb peak) would falsify the gas-phase layer; detecting a localized rise toward tens of ppb at 2.94e-7 bar would confirm it.","tokens_in":27741,"feed_emoji":"🪐","tokens_out":15051,"duration_ms":129723,"temperature":0.7,"pith_summary":"This paper argues that Jupiter's chemistry is one connected system, not separate deep and upper atmospheres: a 1D kinetics code spanning 1.1e3 bar to 7.4e-11 bar can reproduce the main observed chemical features while also making new predictions. Its headline claim is a new stratospheric hydrogen cyanide layer: between 1e-6 and 6.76e-8 bar, photolysis of CH4 and N2 feeds radical reactions with net outcome 2CH4 + N2 -> 2HCN + 3H2, producing 33 parts per billion of HCN at a pressure of 2.94e-7 bar. The same model places a mixed NH3-NH4SH cloud between 0.1 and 1 bar, carries a nitrogen reservoir high into the atmosphere, and anchors sulfur chemistry with updated NH4SH thermodynamics. A reader should care because the HCN peak is a specific, observable target that can be checked with existing ALMA capabilities and upcoming JUICE observations, and because modeling the full atmospheric extent yields predictions that separate regional models miss.","feed_headline":"Jupiter model finds new HCN peak: 33 ppb in stratosphere","feed_subtitle":"Full-atmosphere Jupiter model predicts HCN forms via radicals at 3e-7 bar; JUICE/ALMA can test 33 ppb peak.","key_machinery":"ARGO, a 1D photochemical-thermochemical kinetics solver, running the STAND gas-phase reaction network with the new NH4SH chemistry. At each level ARGO integrates coupled continuity equations dn_i/dt = P_i - L_i - dPhi_i/dz, with eddy and molecular diffusion supplying vertical transport; reverse reaction rates are set by Gibbs free energies, which lets the deep atmosphere relax to thermochemical equilibrium. HCN's stratospheric peak is carried by a radical chain: CH4 photolysis gives CH/CH3, N2 photolysis gives N, then N + CH3 -> H2CN and H2CN -> HCN, with net 2CH4 + N2 -> 2HCN + 3H2. The NH4SH cloud layer is anchored by fitted NASA-polynomial thermodynamics derived from the NH3+H2S equilibri","core_discovery":"The paper's central object is a 1D photochemical-thermochemical kinetics model extending from 1.1e3 bar to 7.4e-11 bar, with sulfur chemistry included over this full span for the first time. After validation against regional models and spacecraft observations, it reproduces Jupiter's main chemistry, including a mixed NH3-NH4SH cloud layer between 0.1 and 1 bar. Its new claim is six vertically separated HCN regimes; in the fifth, between 1e-6 and 6.76e-8 bar, photolysis of CH4 and N2 feeds radicals that make HCN by net 2CH4 + N2 -> 2HCN + 3H2, peaking at 33 ppb at 2.94e-7 bar. The authors present this peak as observationally testable and additionally report CO frozen at 4.2 ppb near 318 bar.","pith_inferences":["The authors do not say so, but a steady background HCN source implies post-Shoemaker-Levy 9 HCN maps need a modeled baseline before cometary material can be cleanly separated; the predicted layer could serve as that baseline.","If the radical mechanism is real, any gas-rich world with overlapping CH4 and N2 photolysis regions should produce a similar HCN layer; locating where the two photolysis cutoffs cross would give a fast, transferable prediction.","Because aerosol chemistry is absent, a non-detection of the 33 ppb peak would more likely mean HCN uptake on haze than failure of the gas-phase chain; fitting an aerosol uptake coefficient to the observed profile would separate the two effects.","The peak's width and altitude encode stratospheric vertical mixing, so high-resolution spectra of the layer could constrain Kzz independently of temperature retrievals even if the absolute abundance remains uncertain."],"forward_implications":["A limb or occultation observation aimed at the 1e-6 to 6.76e-8 bar window can confirm or reject the 33 ppb HCN peak without an entry probe.","The model's gas-phase destruction transition near 2 mbar offers a purely chemical explanation for the sharp HCN drop-off seen in recent ALMA data, independent of the proposed aerosol mechanism.","NH4SH condensation between 0.1 and 1 bar sequesters most sulfur below the cloud, so deep sulfur abundance must be inferred from H2S below that layer, not from upper-atmosphere sulfur.","A CO abundance frozen at 4.2 ppb near 318 bar provides an observable linkage from Jupiter's deep water and oxygen supply to the CO seen higher in the atmosphere.","Having one code from 1,100 bar to 7.4e-11 bar means future retrievals can be compared with a single self-consistent chemical profile, rather than patched regional models."],"supporting_citations":[{"why":"Basis of the ARGO code: supplies the 1D photochemistry-kinetics-diffusion solver, condensation/evaporation routines, and the actinic flux prescription.","marker":"Rimmer and Helling (2016)"},{"why":"Earlier STAND network whose gas-phase reaction set this paper extends with updated rates and NH4SH chemistry.","marker":"Rimmer et al. (2021)"},{"why":"Source of the two updated reaction rate coefficients and a comparative temperate-giant model for HCN and hydrocarbon photochemistry.","marker":"Hu (2021)"},{"why":"Supplies the stratospheric pressure-temperature and eddy diffusion profile and a benchmark for hydrocarbon chemistry.","marker":"Moses et al. (2005)"},{"why":"Supplies the deep-atmosphere temperature profile, Kzz value, and deep chemistry baseline for CO and oxygen species.","marker":"Visscher et al. (2010)"},{"why":"Supplies the elemental abundances used as bottom-boundary initial conditions, derived from Galileo and Juno measurements.","marker":"Rensen et al. (2023)"},{"why":"ALMA retrievals of stratospheric HCN that the model is compared against and that motivate the missing aerosol-chemistry caveat.","marker":"Cavalié et al. (2023b)"},{"why":"Provides the NH4SH equilibrium constant used to fit NASA polynomial coefficients for NH4SH thermodynamics.","marker":"Lewis (1969)"},{"why":"Full-extent disequilibrium model used as a benchmark showing agreement for CH4 and C2 species across the coupled atmosphere.","marker":"Tsai et al. (2021)"}],"fun_headline_variants":["Jupiter model finds HCN peak at 33 ppb","New Jupiter model predicts HCN spike 33 ppb","Jupiter model: HCN radical peak 33 ppb, testable","HCN peak 33 ppb predicted for Jupiter's stratosphere"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The predicted 33 ppb HCN peak assumes gas-phase radical chemistry alone sets HCN between 1e-6 and 6.76e-8 bar; the paper itself notes the model excludes aerosol chemistry, while observations have proposed that aerosols remove HCN.","fun_headline_variants_meta":{"raw":{"variants":["Jupiter model finds HCN peak at 33 ppb","New Jupiter model predicts HCN spike 33 ppb","Jupiter model: HCN radical peak 33 ppb, testable","HCN peak 33 ppb predicted for Jupiter's stratosphere"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001154,"raw_usage":{"total_tokens":4686,"prompt_tokens":876,"completion_tokens":3810,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":620,"completion_tokens_details":{"reasoning_tokens":3736}},"tokens_in":620,"tokens_out":3810,"duration_ms":21663,"temperature":1.0,"reasoning_tokens":3736,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T10:26:04.774868+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A limb or occultation retrieval resolving the 1e-6 to 6.76e-8 bar window that returns an HCN upper limit below a few ppb at ~3e-7 bar (well under the predicted 33 ppb peak) would falsify the gas-phase layer; detecting a localized rise toward tens of ppb at 2.94e-7 bar would confirm it.","supporting_citations":[],"review_version":1}