{"id":"1c2d3a80-f36e-4c80-8e6d-93da5f817f91","arxiv_id":"2412.17028","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Gaseous and solid methane are detected simultaneously toward IRAS 23385+6053, with solid methane best fit by CH4:CO2 = 1:5 and CH4:H2O = 1:10 ices.","lead":"This paper re-analyzes JWST spectra of the star-forming region IRAS 23385+6053 and finds both solid and gaseous methane in the 7.7 micron band. It proposes a new ice composition: methane embedded mostly in carbon dioxide rather than water, plus a tentative nitrous oxide feature.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Gas-phase CH4 detection rests on a 2.3σ single-pixel feature; no integrated detection significance is reported, so the simultaneous gas+ice decomposition is not yet secure.","rationale":"The reader's weakest_assumption pinpoints the same issue: the gas detection is anchored to a 2.3σ feature and an arbitrary S/N=3 area threshold. My stress-test confirms that this is the weakest link because the claim of simultaneity is unique; without the gas detection, the paper would simply be a solid-methane decomposition of the 7.7 µm band, and the quoted gas temperature and column density would be meaningless. The concern is not that the authors are wrong, but that the evidence as presented does not meet the standard for a secure detection. The suggested check (integrated S/N and false-alarm rate) can be done with the public MAST data, so the issue is addressable. I find no additional internal inconsistency, and the authors' honesty about the N2O region is a positive sign. Therefore the appropriate verdict remains CONDITIONAL; no change from the reader's verdict is needed.","tokens_in":14290,"tokens_out":5517,"duration_ms":51056,"concrete_test":"Re-analyze the MAST level-3 MRS cube using the paper's extraction and background subtraction, then measure the S/N of the 1306 cm^-1 line in the extracted spectrum by comparing the integrated line flux to the local noise RMS. Also compute the false-alarm probability via 10^4 Monte Carlo noise realizations of the same spectral window. Report the integrated S/N and false-alarm probability. If the integrated S/N is ≥3 and false-alarm probability <1%, the detection stands; otherwise the gas-phase parameters and the simultaneous decomposition are not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing concern is the statistical support for the gaseous methane detection. In Section 5 the authors state that single pixels show S/N ≤ 2.3 for the strongest feature at 1306 cm^-1, and they set an S/N threshold of 3 to define the lower emitting-area limit. The paper never reports the S/N of the 1306 cm^-1 line after aperture integration, nor a detection significance for the full multi-line fit. Equation (1) uses an optically thin approximation and yields only the product S·N, so the derived N = 0.78e17 cm^-2 and R = 2940 au depend on an area definition that is not statistically justified. If the 1306 cm^-1 peak is a noise fluctuation or a residual of the LOESS continuum subtraction, the gas model in Figure 4 is spurious, and the subsequent subtraction of that model from the optical depth changes the ice decomposition and the CH4:CO2 versus CH4:H2O balance. Thus the paper's central claim of a simultaneous gas and solid methane measurement rests on an unquantified, marginal feature.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes the 7.7 μm region of the JWST MIRI/MRS spectrum of the high-mass star-forming region IRAS 23385+6053. The authors interpret the 7.58–7.8 μm band as a superposition of broad solid methane absorption and sharp gaseous methane emission. Using an LTE slab model with an external line list (TheoReTS) they derive gas parameters R = 2940 au, T = 103(+13/-11) K, N = 0.78(+6.18/-0.64) × 10^17 cm^-2, and using new laboratory ice spectra (ISEAge) plus literature data they derive solid methane column densities for CH4:CO2 = 1:5 at 27.4 K and CH4:H2O = 1:10 at 8.4 K, with a total ice column of 3.99 × 10^17 cm^-2. The paper also tentatively assigns a residual to N2O in various ice environments.","tokens_in":14539,"tokens_out":3574,"duration_ms":34711,"significance":"If the gas-phase detection and the gas–ice decomposition are secure, this would be the first simultaneous quantification of gaseous and solid methane toward IRAS 23385+6053, with implications for methane partitioning and for methane ice environments (favoring CO2-rich rather than water-rich mantles). The paper has clear strengths: it uses publicly available JWST data, combines an external spectroscopic line list with new laboratory measurements, and makes an explicit effort to treat residuals statistically with Anderson–Darling tests. However, the central gas detection is statistically fragile, and the current presentation does not provide an integrated detection significance or a robustness analysis of the continuum subtraction, so the central claim is not yet fully supported.","major_comments":[{"comment":"The gaseous methane detection rests on single-pixel S/N values ≤ 2.3, and the emitting area is defined by adopting an S/N threshold of 3, but no integrated detection significance is reported for the 1306 cm^-1 feature or for the full multi-line fit. Because the subsequent gas–ice decomposition depends on the reality of this emission, please compute and report the S/N after aperture integration over the chosen region, and the detection significance from the full spectral fit. If the integrated feature is not significant at the adopted threshold, the gas detection and the derived R, T, and N values are not secure.","section":"Section 5"},{"comment":"Equation (1) constrains only the product S·N, and the paper separates S from N by thresholding the single-pixel S/N map. The uncertainty in S is not propagated into the reported N = 0.78(+6.18/-0.64) × 10^17 cm^-2, whose upper error bar is almost an order of magnitude larger than the central value. Please quantify the uncertainty on S, propagate it properly into N and T, and show confidence contours for the parameters (e.g., Δχ² maps in the S–N plane).","section":"Equation (1) and Section 6"},{"comment":"The ice continuum is estimated by LOESS smoothing with a 20-point window after masking 'prominent gaseous features and heavy outliers of unknown origin'. This procedure can in principle introduce or suppress narrow residuals at the gas line positions. Please test the robustness of the 1306 cm^-1 feature and the derived gas parameters against the choice of smoothing window and mask selection, or alternatively fit the gas, ice, and continuum simultaneously. Without such a test, the sharp 1306 cm^-1 feature may be a continuum artifact rather than methane emission.","section":"Section 5, continuum subtraction"},{"comment":"The conclusion that solid methane is predominantly in a CO2-rich environment (CH4:CO2 = 1:5) rather than a water-rich environment rests on a restricted library of laboratory mixtures: only CH4:CO2 = 1:5 and CH4:H2O = 1:10, deposited at 6.7 K and warmed at 0.5 K/min, plus literature spectra. Please justify the choice of these specific ratios and temperatures, or test additional compositions (e.g., CH4:CO2 at other ratios, CH4:NH3, or different warm-up rates) to show that the CO2 dominance is not an artifact of the limited library.","section":"Section 6 and Table 1"}],"minor_comments":[{"comment":"The band strengths for CH4:CO2 and CH4:H2O are listed as 'This work', but the derivation of these band strengths is not described in the main text; please add a brief description or a reference to the follow-up paper.","section":"Table 1"},{"comment":"The red circle and red cross in Figure 1 are not defined in the caption; please explain that they indicate the methane-emitting area and the gas emission peak, respectively.","section":"Figure 1 caption"},{"comment":"The phrase 'for the first time gaseous and solid methane were analyzed simultaneously' should be qualified relative to Rocha et al. (2024), who already modeled the ice component toward this source; please clarify that the new contribution is the simultaneous treatment and the gas-phase detection.","section":"Abstract and Introduction"},{"comment":"The text uses the symbol S both for the emitting area in Equation (1) and for signal-to-noise ratio in the S/N discussion; please use a distinct symbol (e.g., A) for the area to avoid ambiguity.","section":"Section 5"},{"comment":"The phrase 'Plank’s law' should be corrected to 'Planck’s law'.","section":"Section 5"},{"comment":"The Anderson–Darling filtering uses α = 12% and 10% for the best fit and confidence regions, respectively; please state whether these thresholds are fixed a priori or chosen to obtain the reported fit, and how the results change with modest variations of α.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The paper is potentially suitable for publication in an astrochemistry/astrophysics journal after the statistical support for the gas detection is strengthened and the robustness of the gas–ice decomposition is demonstrated. The single-pixel S/N of 2.3 for the key gas feature is a genuine concern that the authors must address quantitatively. The ice-fitting library limitation is also important given the paper's central claim of a CO2-rich methane ice environment. I recommend major revision rather than rejection because the overall approach is reasonable and the new laboratory data are valuable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. This is the first simultaneous gas and solid methane analysis for IRAS 23385+6053, and it brings new ISEAge lab spectra of CH4:CO2 and CH4:H2O ices to the problem. The central ice result—that the 7.7 µm band is mostly CH4 in a CO2-rich environment, not the usual water-rich assumption—is plausible and argued with new laboratory data. But the gas-phase detection driving the decomposition is statistically fragile: the strongest single-pixel feature has S/N 2.3, and the emitting area is defined by adopting S/N > 3. That threshold directly sets R = 2940 au, and since the LTE slab model only constrains the product S·N, the reported column density N = 0.78e17 cm^-2 is not independent of the area choice. The paper never reports an integrated detection significance for the full multi-line feature.\n\nCredit where due: the continuum subtraction is handled carefully (LOESS with masked emission lines), the analysis uses external line lists (TheoReTS) and external laboratory data, and the authors are transparent about the low S/N. The Anderson-Darling test on residuals is a good safeguard against overfitting ice mixtures. The total solid CH4 column (3.99e17 cm^-2) matches Rocha et al. within errors, so the ice part is consistent with prior work. The N2O assignment is explicitly tentative, and that caution is appropriate given the missing 4.5 µm feature.\n\nThe soft spots, in order of severity. First, the gas detection needs a proper significance after aperture integration, and the area/column degeneracy should be explored with a S/N map rather than one threshold. Second, the new laboratory spectra are not yet public, though the experimental description is detailed enough that the measurements seem reproducible. Third, the ice fit is limited to two methane mixtures plus LIDA components, so the CO2-rich conclusion is only as strong as the library; the doublet shape does however argue against water-rich ices alone.\n\nThe intended reader is anyone working on interstellar methane ices or JWST ice inventories. The paper deserves a serious referee. I would send it out, asking for a quantitative gas detection significance and a sensitivity analysis of the area assumption, but I see no reason to desk reject.","headline":"A careful and honest lab-plus-JWST analysis worth refereeing, but the gas-phase CH4 column rests on an unquantified 2.3σ detection and a threshold-defined emitting area.","tokens_in":15125,"tokens_out":2416,"would_cite":true,"duration_ms":22950,"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":"The 7.7 µm band toward IRAS 23385+6053 is a superposition of solid methane absorption and warm gaseous methane emission, fit simultaneously with open JWST MIRI/MRS data.","keywords":["methane ice","gaseous methane","IRAS 23385+6053","JWST MIRI/MRS","interstellar ices","infrared spectroscopy","protostar","astrochemistry"],"falsifier":"Take a JWST MIRI/MRS spectrum of the 3.3 µm ν3 methane band toward the same 1.2 arcsec aperture: the gas-phase interpretation predicts an emission feature at the same temperature and column density, while a solid-state origin for the 1306 $cm^{-1}$ feature would show a different or absent 3.3 µm feature.","tokens_in":14102,"feed_emoji":"🌌","tokens_out":9026,"duration_ms":76311,"temperature":0.7,"pith_summary":"This paper reinterprets the 7.7 µm band toward the high-mass star-forming region IRAS 23385+6053, a feature usually attributed to methane ice alone. Using public JWST MIRI/MRS spectra, the authors argue the band is actually a wide solid-state absorption feature of methane ice overlapped by sharp emission lines from warm gaseous methane. They report the first simultaneous quantification of both phases in this source: gas at temperature 103 K, column density 0.78 × $10^{17}$ $cm^{-2}$, and emitting radius 2940 au, plus solid methane with total column density 3.99 × $10^{17}$ $cm^{-2}$. The ice is best reproduced by laboratory mixtures of CH4:CO2 = 1:5 at 27.4 K and CH4:H2O = 1:10 at 8.4 K, implying that methane ice in this object is predominantly carbon-dioxide-rich rather than water-rich. A residual near 1283–1297 $cm^{-1}$ is tentatively assigned to nitrous oxide in various ice environments.","feed_headline":"JWST data split methane's 7.7 micron band into gas and ice","feed_subtitle":"First simultaneous gas-and-ice measurement finds 103 K gas plus CO2-rich methane ice toward a high-mass protostar.","key_machinery":"The band decomposition rests on a two-step fitting procedure. First, an LTE slab emission model built from the TheoReTS methane line list and partition function is fit to the flux spectrum after masking the ice continuum; the emitting area is fixed from the flux map using a S/N = 3 aperture, and temperature and column density come from χ² minimization with confidence regions from Δχ² maps. Then the solid-state part is fit in optical depth space with a linear combination of laboratory spectra of CH4-bearing ices made on the ISEAge setup — CH4:CO2 = 1:5 and CH4:H2O = 1:10 deposited at 6.7 K and warmed at 0.5 K per minute — together with OCN⁻, CH3CH2OH:H2O, and SO2:CH3OH reference spectra. The laboratory CH4:CO2 mixture uniquely reproduces the doublet width and both wings, and that match carries the CO2-rich conclusion.","core_discovery":"The central claim is that the resolved doublet structure of the 7.58–7.8 µm methane ν4 deformation band in IRAS 23385+6053 is not a single ice feature. The sharp component at 1306 $cm^{-1}$ coincides with a ν4 gaseous methane transition and is modeled with an LTE slab emission model, while the broad wings are solid methane in a CO2-dominated matrix with a smaller water-rich component. The best-fit gas parameters are R = 2940 au, T = 103(+13/-11) K, and N = 0.78(+6.18/-0.64) × $10^{17}$ $cm^{-2}$; the solid methane column densities are 2.97 × $10^{17}$ $cm^{-2}$ in CH4:CO2 = 1:5 at 27.4 K and 1.02 × $10^{17}$ $cm^{-2}$ in CH4:H2O = 1:10 at 8.4 K, for a total solid column of 3.99 × $10^{17}$ $cm^{-2}$. The authors conclude that methane ice is mostly surrounded by CO2 rather than H2O, contrasting with the water-methane mixtures commonly used in previous literature.","pith_inferences":["If the decomposition holds, the same doublet morphology in other protostellar MIRI/MRS spectra could serve as a general diagnostic for separating warm methane gas from methane ice, and the CO2-dominated environment would point to methane formation or processing in CO2-rich mantles rather than water-rich ones.","The laboratory result that a N2O:CO2 = 1:20 mixture shifts the N2O feature by 14.4 cm^-1 provides a diagnostic for identifying N2O in CO2-rich ices; future JWST observations of the stronger 4.5 µm N2O band would turn the tentative 7.7 µm detection into a secure one.","Because the gas-phase column density has very asymmetric error bars (0.78 +6.18/-0.64 × 10^17 cm^-2), the gas-to-ice methane ratio is not yet tightly constrained; a higher-S/N observation of the same source would make that ratio a useful probe of methane desorption and warm carbon-chain chemistry."],"forward_implications":["With the gas-phase contribution removed, the residual 7.7 µm band is a single broad ice absorption that can be described by CH4:CO2 = 1:5 at 27.4 K plus CH4:H2O = 1:10 at 8.4 K, so solid methane toward this source is predominantly in a CO2-rich rather than water-rich environment.","The total solid CH4 column density of 3.99 × 10^17 cm^-2 gives Nice(CH4)/Nice(H2O) ≈ 2.44%, consistent with the earlier estimate of 2.72%, meaning the new decomposition preserves the overall methane abundance while changing its phase split.","Gaseous methane is spatially offset from the continuum peak and is present only near the most massive source, implying the emitting gas sits at the edge of the accretion disk or in the envelope rather than in the cold core.","If this doublet structure appears in other JWST spectra, the same two-step procedure offers a way to separate gas and ice methane columns without assuming the whole ν4 feature is solid-state.","The unassigned 1283–1297 cm^-1 residual is tentatively attributed to N2O in various ice matrices, making IRAS 23385+6053 a new candidate for interstellar nitrous oxide."],"supporting_citations":[{"why":"Supplies the calibrated JWST MIRI/MRS spectrum, the adopted global continuum, silicate subtraction, and the prior 6.8–8.6 µm COM analysis that this paper reinterprets.","marker":"Rocha, W. R. M. et al. 2024"},{"why":"Establishes the first secure detection of interstellar gaseous methane and possible solid methane, anchoring the ν4 band identification.","marker":"Lacy et al. 1991"},{"why":"Documents the doublet structure of methane ice in a CO2 matrix on warming, the precedent for the CH4:CO2 assignment.","marker":"Boogert et al. 1997b"},{"why":"Provides band strengths and amorphous solid methane spectra needed to convert fitted optical depths into column densities.","marker":"Gerakines & Hudson 2015"},{"why":"Describes the ISEAge laboratory setup used to grow and measure the CH4:CO2 and CH4:H2O ice mixtures that compose the solid-state fit.","marker":"Ozhiganov et al. 2024"},{"why":"Supplies the TheoReTS methane line list and partition function used to model the gaseous emission spectrum at 103 K.","marker":"Rey et al. 2013"},{"why":"Gives the Δχ² method used to construct the reported 1σ confidence intervals on gas and ice parameters.","marker":"Avni 1976"},{"why":"Represents the Spitzer c2d water-methane ice paradigm and the CH4–CO2 correlation that this work contrasts with a CO2-dominated fit.","marker":"Oberg et al. 2008"},{"why":"Sets the adopted 4.9 kpc distance and astrometry that enter the emitting-area and column-density estimates.","marker":"Molinari et al. 1998, 2002"}],"fun_headline_variants":["JWST untangles methane gas and ice in one infrared band","First joint gas-ice methane detection in a massive protostar","CO2-rich methane ice unveiled by JWST around a protostar","JWST sees methane both frozen and gaseous in a protostar"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole gas-ice decomposition hangs on the sharp feature near 1306 $cm^{-1}$ being real gaseous methane emission; if it is instead a solid-state band, a noise fluctuation, or a residual continuum artifact, the gas parameters and the phase split collapse.","fun_headline_variants_meta":{"raw":{"variants":["JWST untangles methane gas and ice in one infrared band","First joint gas-ice methane detection in a massive protostar","CO2-rich methane ice unveiled by JWST around a protostar","JWST sees methane both frozen and gaseous in a protostar"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00067,"raw_usage":{"total_tokens":3243,"prompt_tokens":1323,"completion_tokens":1920,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":939,"completion_tokens_details":{"reasoning_tokens":1845}},"tokens_in":939,"tokens_out":1920,"duration_ms":14176,"temperature":1.0,"reasoning_tokens":1845,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:50:30.636056+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a JWST MIRI/MRS spectrum of the 3.3 µm ν3 methane band toward the same 1.2 arcsec aperture: the gas-phase interpretation predicts an emission feature at the same temperature and column density, while a solid-state origin for the 1306 $cm^{-1}$ feature would show a different or absent 3.3 µm feature.","supporting_citations":[{"cited_title":"A., & Hudson, R","cited_arxiv_id":null,"evidence_quote":"Provides band strengths and amorphous solid methane spectra needed to convert fitted optical depths into column densities."},{"cited_title":"V., & Tyuterev, V","cited_arxiv_id":null,"evidence_quote":"Supplies the TheoReTS methane line list and partition function used to model the gaseous emission spectrum at 103 K."},{"cited_title":"1998, , 505, L39, 10.1086/311591","cited_arxiv_id":null,"evidence_quote":"Sets the adopted 4.9 kpc distance and astrometry that enter the emitting-area and column-density estimates."}],"review_version":1}