{"id":"93853534-178a-4861-843b-1d4ff960ded7","arxiv_id":"2506.01771","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"JWST observations of WASP-121b reveal SiO and a super-stellar C/O ratio on the dayside, plus methane on the nightside, indicating a volatile-rich envelope shaped by vertical mixing.","lead":"Astronomers report detections of water, carbon monoxide, silicon monoxide, and methane in the atmosphere of the ultrahot giant planet WASP-121b from a single JWST observation. The planet's carbon, oxygen, and silicon abundances are higher than its host star's, pointing to formation through both drifting pebbles and rocky planetesimals.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the stellar O non-LTE normalization shifts the quoted enrichment factors but not the super-stellar direction.","rationale":"The paper's central claims are the dayside detections of H2O, CO, and SiO, the nightside CH4 detection, and super-stellar C/H, O/H, Si/H, and C/O. The detection claims are unusually well supported: five independent retrieval codes, two independent reductions, explicit with/without-molecule Bayesian comparisons, and consistency of derived elemental ratios across free-chemistry and equilibrium retrievals. The weakest link is the stellar O abundance normalization, exactly as the reader identified. The non-LTE correction is large, but the direction of the super-stellar C/O and O/H results survives even in the LTE limit. To invert the C/O conclusion, the non-LTE correction would need to be wrong by roughly another 0.33 dex in the opposite sense, which is not suggested by any available evidence. Thus the concern is not load-bearing. The only honest stress-test outcome is that the paper's central argument holds; the recommended verdict is unchanged from ACCEPT. The proposed concrete test is still worth running because it would place the enrichment factors on firmer footing, but it would not be expected to change the qualitative conclusions.","tokens_in":39530,"tokens_out":9383,"duration_ms":109081,"concrete_test":"Recompute the WASP-121 stellar C and O abundances with an independent non-LTE grid and with additional O I lines, then re-derive (C/O)/(C/O)* and (O/H)/(O/H)* using both LTE and non-LTE values. If the ratios remain above unity in all plausible cases, the central super-stellar C/O and O/H claims are insensitive to the stellar O normalization.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The only candidate concern that could move the headline enrichment factors is the host-star oxygen normalization. The non-LTE O I correction is large (A(O) = 9.238 LTE to 8.855 non-LTE), and it changes the reported (O/H)/(O/H)* from roughly 5 to 12.2 and (C/O)/(C/O)* from roughly 4.5 to 1.96. However, the central claim of super-stellar C/O and O/H does not depend on this correction: using the paper's LTE stellar C/O of 0.203, the planetary C/O of 0.92 is still super-stellar by a factor of about 4.5, and using LTE stellar O still gives (O/H)/(O/H)* > 1. The C/H and Si/H ratios are almost unaffected by the correction. The molecular detections are supported by multiple independent retrieval codes and two independent data reductions; the ATMO equilibrium retrieval constrains elemental abundances directly, and the derived C/O, C/Si, and O/Si ratios are consistent across the free-chemistry retrievals. I therefore find no load-bearing internal weakness: the stellar non-LTE O normalization is a real systematic, but it changes the magnitude of the enrichment factors, not the qualitative super-stellar claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents JWST NIRSpec G395H phase-curve observations of WASP-121b. From the dayside disc-integrated spectrum the authors report detections of H2O (5.5–13.5σ), CO (10.8–12.8σ), and SiO (5.7–6.2σ), and from the nightside spectrum a CH4 detection (3.1–5.1σ). Retrievals with five independent codes, including an equilibrium-chemistry retrieval (ATMO), yield super-stellar atmospheric C/H, O/H, and Si/H ratios relative to the host star, with C/O = 0.92, corresponding to a C/O ratio 1.96 times the adopted non-LTE stellar value. The authors interpret the super-stellar C/O and volatile/refractory ratios as evidence for pebble-driven gas enrichment and post-formation accretion of rocky material, and propose vertical mixing to explain the nightside CH4 abundance.","tokens_in":39808,"tokens_out":7949,"duration_ms":80546,"significance":"The result is significant because it demonstrates, from a single JWST observation, simultaneous constraints on a refractory species (SiO) and major volatiles (H2O, CO, CH4), opening a new window on giant-planet formation. The analysis is unusually thorough: five independent retrieval codes (ATMO, NEMESIS, CHIMERA, HyDRA, PETRA) give consistent elemental ratios; two independent data reductions (FIREFly and Eureka!) agree; detection significances are quantified via Bayesian evidence; and line-list differences are explicitly discussed. Data and code are public. The central caveat is the large non-LTE correction to the stellar oxygen abundance, which changes the magnitude of the quoted enrichment factors (e.g., (O/H)/(O/H)* from roughly 5 to 12 and (C/O)/(C/O)* from roughly 4.5 to 1.96) but does not reverse the qualitative super-stellar conclusion; the manuscript discloses this sensitivity in the Methods.","major_comments":[],"minor_comments":[{"comment":"The quoted enrichment factors, such as (O/H)/(O/H)star = 12.19 and (C/O)/(C/O)star = 1.96, assume non-LTE stellar oxygen abundances; because the O I non-LTE correction is large, please state explicitly in the main text that adopting the LTE stellar C/O of 0.203 instead of the non-LTE value of 0.470 would change (C/O)/(C/O)star from about 2 to about 4.5, and similarly for (O/H)/(O/H)star, while still leaving the ratios super-stellar.","section":"Results (first paragraph of 'Results')"},{"comment":"The estimate of about 21 Earth masses of accreted rocky material depends directly on the assumed core mass (15 MEarth), the envelope hydrogen mass fraction (75%), and the assumption that the primordial envelope had a stellar Si/H ratio; the text gives one alternative (23.7 MEarth) for zero initial silicon, but it would be helpful to state how the result scales with a plausible range of core masses from the literature.","section":"Methods — Accretion of rocky material"},{"comment":"The nightside CH4 detection significance ranges from 3.1σ to 5.1σ across the retrieval codes, partly because of differing CH4 line lists; the paper should state explicitly that this spread is a systematic uncertainty from the opacity treatment, not just a statistical fluctuation, so that readers do not over-rely on the favored 4.7σ value.","section":"Methods — Atmospheric retrieval analyses"},{"comment":"The inferred eddy diffusion coefficients Kzz of roughly 10^8–10^11 cm2/s are compared only with broad literature constraints; if possible, include a comparison with the Kzz profile expected from the retrieved nightside thermal structure or a simple mixing-length estimate, which would strengthen the physical plausibility argument.","section":"Discussion (vertical mixing paragraph)"}],"recommendation":"minor_revision","confidential_remarks":"This is a strong empirical paper with robust multi-code, multi-reduction validation. The only substantive systematic is the non-LTE stellar oxygen normalization; the authors already report both LTE and non-LTE stellar values, and the qualitative super-stellar C/O conclusion survives both. I recommend minor revision to make this sensitivity more prominent in the main text. No concerns about novelty, citation practices, or scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a strong paper and a genuine advance. It gives the first simultaneous refractory (SiO) and volatile (H2O, CO, CH4) abundance constraints from a single JWST observation of an exoplanet, and the infrared SiO detection is new and convincing. The super-stellar C/O for WASP-121b was already reported, but this paper adds Si/H and the nightside CH4 detection, and the combination is what matters for formation studies. The analysis is careful: five independent retrieval codes agree on the key ratios (C/O, C/Si, O/Si), two independent reductions give spectra consistent to better than 2 sigma, and the detection significances are backed by Bayesian evidence with a sensible treatment of line-list differences. The nightside CH4 detection is softer (3.1-5.1 sigma) and partly line-list dependent, but it is supported by multiple codes and by the consistency check with the dayside C/O: at C/O ~ 0.92, equilibrium chemistry predicts a CH4-rich nightside, so the detection is not isolated. I do not find the circularity concern compelling; the paper uses the dayside-derived C/O as a plausibility check, not as a retrieval input. The softest spot is the host-star oxygen abundance: the non-LTE correction is large (A(O) drops from 9.238 to 8.855), and the quoted enrichment factors ((O/H)/(O/H)_star = 12.19 and (C/O)/(C/O)_star = 1.96) scale with it. But using their LTE stellar C/O of 0.203, the planetary C/O of 0.92 is still super-stellar by a factor of about 4.5, so the direction is robust. The rocky-material estimate (15-33 Earth masses) rests on a simple one-zone calculation with assumptions about core mass, mixing, and Earth-like composition; it is clearly labeled as interpretive and should not be read as a precision measurement. The formation narrative (pebbles plus planetesimals, vertical mixing) is speculative, but the paper says so and cites competing scenarios. This paper is for anyone working on hot Jupiter atmospheres, retrieval methodology, or giant planet formation. It deserves a serious referee. My recommendation: accept after minor revisions, mainly asking the authors to state more explicitly how much the enrichment factors depend on the stellar non-LTE corrections and to report LTE-normalized values in an appendix.","headline":"Strong paper: first simultaneous refractory+volatile abundances from one JWST phase curve, with a credible new SiO detection; the stellar O non-LTE caveat shifts enrichment magnitudes but not the super-stellar conclusion.","tokens_in":40424,"tokens_out":1756,"would_cite":true,"duration_ms":17633,"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":"JWST finds SiO and a C/O ratio 1.96 times stellar on WASP-121b.","keywords":["exoplanet atmospheres","ultrahot Jupiter","WASP-121b","JWST NIRSpec phase curve","silicon monoxide","carbon-to-oxygen ratio","refractory elements","vertical mixing"],"falsifier":"Recompute the WASP-121 stellar oxygen abundance with an independent non-LTE or 3D model grid: if the resulting stellar C/O rises to near the planetary 0.92, the central claim that WASP-121b has a super-stellar C/O loses its quantitative support. A second check is to observe the nightside at higher spectral resolution or longer wavelengths to confirm that the 3.3 micron feature is CH4 rather than a P-T profile artifact.","tokens_in":39365,"feed_emoji":"🪐","tokens_out":6367,"duration_ms":61108,"temperature":0.7,"pith_summary":"The paper claims that a single 37.8-hour JWST/NIRSpec phase curve of the ultrahot giant WASP-121b carries simultaneous detections of volatile molecules (H$_2$O and CO) and a refractory molecule (SiO) in the dayside atmosphere, plus CH$_4$ on the nightside. From those detections it derives super-stellar C/H, O/H, Si/H, and a C/O ratio 1.96 times the host star's value, arguing the planet built its envelope from volatile-enriched gas and later rocky additions. If correct, this is the first abundance set that jointly constrains refractory and volatile elements from one exoplanet observation, and it ties the presence of nightside CH$_4$ and the survival of dayside SiO to strong vertical mixing rather than horizontal quenching.","feed_headline":"SiO and a 1.96x stellar C/O seen in one JWST phase curve of WASP-121b","feed_subtitle":"Single 37.8-hour orbit yields H2O, CO, SiO on the dayside and CH4 on the nightside, with super-stellar C, O, Si abundances.","key_machinery":"The load-bearing object is the phase-resolved emission spectrum built from 349 wavelength channels of the JWST/NIRSpec G395H phase curve, fitted with a spherical-harmonic brightness map and a quadratic limb-darkening law. The retrieval machinery is a set of five independent atmospheric retrieval codes (ATMO, NEMESIS, CHIMERA, HyDRA, PETRA) that convert those spectra into abundances; the central identifications are the 4.0$-$4.3 $\\mu$m SiO band on the dayside and the 3.3 $\\mu$m CH$_4$ features on the nightside. The vertical-mixing mechanism is quantified with a C$-$H$-$N$-$O$-$S non-equilibrium chemical network, whose eddy diffusion coefficient $K_{zz}$ is matched to the observed nightside CH$_4$ abundance, and the stellar normalization uses non-LTE abundance corrections for the host star.","core_discovery":"Using the G395H grism of JWST's NIRSpec instrument to watch WASP-121b for 37.8 hours, the authors extract dayside and nightside emission spectra from the phase curve and fit them with five independent retrieval codes. They report decisive detections of dayside H$_2$O (5.5$-$13.5$\\sigma$), CO (10.8$-$12.8$\\sigma$), and SiO (5.7$-$6.2$\\sigma$), and nightside CH$_4$ (3.1$-$5.1$\\sigma$). Relative to updated stellar abundances, the dayside ratios are $(\\mathrm{C/H})/(\\mathrm{C/H})_\\star=23.96$, $(\\mathrm{O/H})/(\\mathrm{O/H})_\\star=12.19$, $(\\mathrm{Si/H})/(\\mathrm{Si/H})_\\star=9.89$, and C/O $=0.92$, which is 1.96 times the stellar C/O of 0.47; all four ratios are super-stellar at $>99.99\\%$ probability. The same data show a CH$_4$-rich nightside that cannot be explained by equilibrium chemistry alone, and the authors argue that vertical mixing ($K_{zz}\\sim10^9$-$10^{11}\\,\\mathrm{cm^2\\,s^{-1}}$) transports CH$_4$-rich gas from deep layers to the nightside photosphere while also keeping silicon from being cold-trapped, preserving the dayside SiO signal.","pith_inferences":["Even if the non-LTE correction to the stellar oxygen abundance is wrong, the LTE stellar C/O of 0.203 still leaves the planetary C/O of 0.92 super-stellar, so the qualitative formation story survives; only the magnitude of the enrichment changes.","The paper reports a wavelength-dependent phase offset whose flat-line fit gives $\\chi^2_\\nu=1.74$, hinting at longitudinal or spectral variations in the brightness map that are not interpreted; a dedicated $\\Delta\\phi(\\lambda)$ mapping could test the vertical-mixing picture.","The nightside CH$_4$ abundance predicted by vertical mixing could be verified by observing WASP-121b at longer wavelengths where CH$_4$ bands are stronger, or by comparing the inferred $K_{zz}$ with values from other tracers such as CO/CO$_2$ ratios.","If SiO is a strong shortwave absorber, the retrieved Si abundance should affect the altitude of the dayside thermal inversion; coupling that abundance into self-consistent GCMs would test whether the inversion pressure matches the data."],"forward_implications":["If the detections hold, WASP-121b becomes the first exoplanet whose atmosphere has been measured for both refractory (Si) and volatile (C, O) elements in a single observation, tying SiO to the dayside thermal inversion.","The super-stellar C/H and C/O support formation by gas enriched with inward-drifting, CH$_4$-rich pebbles evaporating beyond the H$_2$O ice line, while the super-stellar Si/H requires an additional roughly 21 Earth masses of rocky planetesimals.","The nightside CH$_4$ detection implies strong vertical mixing ($K_{zz}\\sim10^9$-$10^{11}\\,\\mathrm{cm^2\\,s^{-1}}$) rather than horizontal quenching, and predicts that the same mixing prevents refractory cold-trapping, explaining why SiO is seen on the dayside.","Future JWST observations should find SiO bands in other ultrahot giants with silicon- and oxygen-enriched atmospheres, complementing UV searches that suffer from dust extinction and low UV flux.","The measured C/O of 0.92 places the nightside opacity under CH$_4$ control, explaining the non-detection of nightside H$_2$O and CO."],"supporting_citations":[{"why":"Provides the JWST NIRSpec phase-curve observation and the FIREFly reduction approach on which this study's data products and white-light fit are built.","marker":"[5]"},{"why":"Supplies the GCM spectra and thermal-dissociation parameterization used to interpret muted H2O and SiO dayside features and to compare nightside P-T profiles.","marker":"[11]"},{"why":"Previous near-UV SiO candidate in WASP-121b and WASP-178b that this NIRSpec detection verifies unambiguously.","marker":"[13]"},{"why":"Pebble drift/evaporation model invoked to explain the super-stellar C/H and C/O as volatile-rich gas accretion beyond the H2O ice line.","marker":"[22]"},{"why":"Pebble accretion simulations predicting super-stellar C/H, O/H, C/O with high (C+O)/Si, used to argue that refractory solids are needed to explain Si/H.","marker":"[24]"},{"why":"C-H-N-O-S non-equilibrium chemical network (VULCAN) used to model vertical mixing and reproduce the nightside CH4 abundance.","marker":"[35]"},{"why":"Latest ExoMol SiO line list used by HyDRA, giving the preferred 6.2-sigma SiO detection significance.","marker":"[82]"},{"why":"HITEMP CH4 line list used by CHIMERA, the preferred nightside CH4 detection.","marker":"[78]"},{"why":"Non-LTE abundance correction grids applied to stellar C, O, Al, Ca, Si, K, Fe; the oxygen correction sets the stellar C/O used for the super-stellar ratios.","marker":"[159-162]"}],"fun_headline_variants":["JWST finds SiO and super-stellar C/O on WASP-121b","One JWST orbit reveals SiO and 1.96x stellar C/O on WASP-121b","Super-stellar C/O with SiO in WASP-121b's JWST phase curve","SiO detection and super-stellar C/O ratio in WASP-121b atmosphere","WASP-121b's air: JWST sees SiO and C/O 1.96x stellar"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The enrichment factors rest on the non-LTE correction to the host star's oxygen abundance; if that correction is wrong, the sizes of the super-stellar O/H and C/O ratios change, though their direction survives even with LTE stellar abundances.","fun_headline_variants_meta":{"raw":{"variants":["JWST finds SiO and super-stellar C/O on WASP-121b","One JWST orbit reveals SiO and 1.96x stellar C/O on WASP-121b","Super-stellar C/O with SiO in WASP-121b's JWST phase curve","SiO detection and super-stellar C/O ratio in WASP-121b atmosphere","WASP-121b's air: JWST sees SiO and C/O 1.96x stellar"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001036,"raw_usage":{"total_tokens":4455,"prompt_tokens":1132,"completion_tokens":3323,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":748,"completion_tokens_details":{"reasoning_tokens":3205}},"tokens_in":748,"tokens_out":3323,"duration_ms":22300,"temperature":1.0,"reasoning_tokens":3205,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:34:13.831313+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the WASP-121 stellar oxygen abundance with an independent non-LTE or 3D model grid: if the resulting stellar C/O rises to near the planetary 0.92, the central claim that WASP-121b has a super-stellar C/O loses its quantitative support. A second check is to observe the nightside at higher spectral resolution or longer wavelengths to confirm that the 3.3 micron feature is CH4 rather than a P-T profile artifact.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"HITEMP CH4 line list used by CHIMERA, the preferred nightside CH4 detection."}],"review_version":1}