{"id":"89daed99-b736-404e-8def-00709aac0351","arxiv_id":"2502.10271","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"JWST spectroscopy of Arp 220's western nucleus reveals hot-core chemistry and shocks in dense molecular gas, with no direct evidence for an embedded AGN.","lead":"Using JWST spectrographs, the authors separated the two nuclei of Arp 220 and read the molecular absorption signature of its western nucleus. The data reveal a hot-core-like, shocked chemical environment around a compact hot dust core, with no clear sign of a hidden black hole.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative column densities and the 330 K HCN component depend on pure-absorption LTE fits with a scalar background fraction; the paper's own evidence for line emission and the unresolved 7/14 µm factor-2 discrepancy leave Table 1 systematics open, though the order-of-magnitude result is…","rationale":"The reader's weakest-assumption analysis correctly identifies the pure-absorption LTE assumption with a scalar background fraction as the key fragile point. I agree that this is the most load-bearing concern because the reported column densities, the warm HCN component, and the inferred foreground-dilution geometry all depend on it. The paper is honest about the limitation, and the authors even argue that line emission would make the column densities lower limits, which protects the order-of-magnitude claim against one direction of systematic error. However, the same non-LTE effects could alter the cold/warm HCN decomposition and the dilution factor used to infer the radiation field, and the unresolved factor-1.6–2.0 discrepancy between the 7 µm and 14 µm bands (Table 1) shows that the LTE framework is not quantitatively closed. Because the qualitative conclusions—hot-core-like chemistry, shocks, and no clear AGN signature in the absorbing gas—are robust to these systematics, the verdict should remain CONDITIONAL rather than being strengthened or reversed. The proposed non-LTE multi-band test would directly settle whether the quantitative claims survive a more physical treatment.","tokens_in":30477,"tokens_out":14640,"duration_ms":156010,"concrete_test":"Run a non-LTE rovibrational radiative transfer model for the HCN ν2 and 2ν2 bands and the C2H2 ν5 and ν4+ν5 bands simultaneously, with a hot (≳500 K) background continuum, absorbing/emitting molecular gas, and foreground dust, letting f_bg vary with wavelength. Require the model to reproduce the observed 7 and 14 µm band depths and column-density ratio, the Q-branch shapes, and the f_bg values at 14 and 21 µm. If the 330 K HCN component persists with N(HCN) ≳ 10^17 cm^-2 and a single dust geometry reproduces the f_bg trend, the central claims stand; if not, the LTE/scalar-f_bg framework is inadequate and the quantitative conclusions in Table 1 would need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims rest on the LTE, pure-absorption model with a single scalar background fraction f_bg per band (Section 3, Section 4.1.2). This assumption is load-bearing because f_bg = 0.28 for HCN/C2H2 at 14 µm and the derived column densities are what produce the order-of-magnitude increase over Spitzer (Section 4.3). The paper itself provides evidence that the assumption is imperfect: the P/R-branch asymmetries in H2O, CO, and the tentatively identified HCN 3ν2−ν2 band indicate line emission partially filling the absorption lines (Sections 3.2, 3.3, 3.4.2). If line emission is significant, the fitted f_bg values are lower than the true values and the column densities are lower limits, so the order-of-magnitude claim is not invalidated by emission alone. However, the same effect can change the decomposition into cold and warm HCN components and the inferred dilution factor from Eq. 1, which underpin the claim of a radiatively excited 330 K component tracing a hot Td > 500 K inner core. A second internal tension is that the 7 µm HCN and C2H2 bands yield column densities a factor 1.6–2.0 higher than the 14 µm bands (Table 1), and the quantitative toy model of Section 4.1.1 is not shown to reproduce the observed f_bg values at both 14 and 21 µm. Thus the absolute column densities and the warm-HCN interpretation are plausible but not fully secured until these systematics are quantified with a non-LTE, multi-band model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents 3–28 µm JWST MIRI/MRS and NIRSpec/IFU spectra of the western nucleus of Arp 220, the first spatially resolved mid-infrared spectroscopic view of this nucleus. The authors identify and model a large set of molecular rovibrational absorption bands, using LTE spectral models with MCMC fitting and a single scalar background fraction f_bg per band. They derive column densities and rotational temperatures for 10 species. The main claims are: (i) HCN, C2H2 and CO2 column densities an order of magnitude higher than previous Spitzer estimates because the Q-branches are optically thick and f_bg≈0.28; (ii) a warm 330 K HCN component radiatively excited by a hot (Td>500 K) embedded core; and (iii) a chemical footprint resembling Galactic hot cores with additional shock tracers, and no evidence for AGN-driven X-ray chemistry. The paper is careful in stating its assumptions and explicitly discusses the possibility of line emission filling in the absorption bands.","tokens_in":30882,"tokens_out":12927,"duration_ms":111038,"significance":"If correct, these results provide the deepest mid-infrared view of the Arp 220 western nucleus and establish rovibrational absorption spectroscopy as a powerful probe of compact obscured nuclei. The order-of-magnitude increase over Spitzer would revise our picture of the molecular content of ULIRG nuclei. The C2H2 ortho-to-para ratio of 1.71±0.05 is the first extragalactic measurement and offers a new diagnostic of nuclear chemistry. The absence of X-ray-driven chemical signatures in gas close to the putative buried AGN is an important, albeit not decisive, contribution to the AGN/starburst debate. The paper's strengths include the careful use of MCMC fitting, the explicit treatment of the degeneracy between f_bg and column density, and the transparent discussion of systematic uncertainties, including the self-identified limitations of the pure-absorption LTE approach.","major_comments":[{"comment":"The pure-absorption LTE assumption with a single scalar background fraction is load-bearing for the central quantitative results. The paper itself provides evidence that this assumption is imperfect: P/R-branch asymmetries in H2O (Fig. 4) and CO (Fig. 5), and the tentative HCN 3ν2−ν2 hot band (Fig. 7) indicate line emission partially filling the absorption. In Section 4.1.2 the authors show that such emission can fill in the 14 µm ν2 band while leaving the 7 µm band nearly unchanged. Under these conditions the fitted f_bg values are lower limits, and the decomposition into 50 K and 330 K HCN components in Table 1 may not be unique. Because the warm 330 K component and the dilution factors derived from Eq. (1) in Section 4.3 underlie the main geometric conclusion, the paper would be substantially strengthened by a non-LTE or two-layer model, or by an explicit sensitivity test showing how the derived parameters change when emission filling is included.","section":"Section 3; Section 4.1.2; Table 1; Eq. (1)"},{"comment":"The factor 1.6–2.0 discrepancy between the 7 µm and 14 µm column densities of HCN and C2H2 is acknowledged but not resolved. The toy model of Section 4.1.1 (Figs. 11–13) is invoked to explain the difference qualitatively, but it is not fitted to the data and does not demonstrate that a single set of physical parameters can reproduce the observed f_bg at 7, 14 and 21 µm simultaneously. Since the absolute column densities in Table 1 are central to the order-of-magnitude claim relative to Spitzer and to the abundance ratios in Section 4.4, the systematic uncertainty of roughly a factor of two should be propagated into the quoted column densities and the derived ratios.","section":"Section 4.1.1; Table 1; Sections 3.1 and 3.4.2"},{"comment":"The derivation of the dilution factor f_D = 0.2–0.5 assumes that the 330 K rotational temperature of HCN is set entirely by radiative pumping from a blackbody with T_bg = 600–1500 K. The paper argues that collisional excitation is unlikely because of the high critical densities of the J≥10 levels, but it does not quantify the competition between collisional and radiative pumping, including the effect of photon trapping which it mentions. A simple LVG or two-level estimate would show whether the derived f_D and the resulting statement that the absorbing gas lies 'at or near the edge of the optically thick background source' are robust. Without this, the geometric interpretation is plausible but not fully secured.","section":"Section 4.3; Eq. (1)"},{"comment":"The H2O column density is listed in Table 1 and used in the chemical comparison of Section 4.4 (e.g., the [H2O]/[CO] ratio of ~1), despite the fact that the LTE model fails to reproduce the R-branch lines (Fig. 4) and the authors state that far-infrared radiative excitation is important for H2O. The quoted N(H2O) = (2.4–7.0)×10^18 cm^-2 should therefore be regarded as a rough estimate with large systematic uncertainty. Since H2O is one of the two most abundant molecules in the sample, this uncertainty propagates into the hot-core/shock chemistry conclusions. At minimum, the paper should either restrict the H2O-based claims to a lower limit or present a non-LTE excitation model for H2O.","section":"Section 3.2; Table 1; Section 4.4"}],"minor_comments":[{"comment":"The statement in Section 3.1 that 'we do not detect any clear hot band absorption' may be read as contradicting the tentative 3ν2−ν2 hot band detection in Section 3.4.2; clarifying that the former refers specifically to the 14 µm ν2 fundamental region would avoid confusion.","section":"Section 3.1; Section 3.4.2"},{"comment":"For CO, the text states that the cold and warm components are not uniquely determined by the fits because they vary strongly with small changes to the assumed kinematics, but the table lists only narrow ranges for their column densities and temperatures. It would be helpful to indicate in the table that these ranges do not include the kinematic degeneracy.","section":"Table 1; Section 3.3"},{"comment":"The abundance ratios in Section 4.4 compare CO, whose column density is measured from the 0.15 arcsec NIRSpec aperture (Section 2.2), with species measured from the 0.435 arcsec MIRI aperture. Since the CO absorption is known to be aperture-dependent, the aperture mismatch should be stated when presenting ratios involving CO.","section":"Section 4.4; Section 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a strong observational paper well matched to A&A. The authors are transparent about the limitations of their LTE pure-absorption modeling, and the main concerns raised in this report are about the quantitative impact of those limitations on the central claims. I see no novelty or attribution issues."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this is a genuinely New dataset and the analysis is careful. JWST separates the two Arp 220 nuclei across 3-28 um, and the molecular absorption bands are fitted with LTE models plus MCMC, with explicit discussion of degeneracies. The headline results—HCN, C2H2, CO2 columns an order of magnitude above Spitzer, a 330 K HCN component, first extragalactic C2H2 ortho-to-para ratio of 1.71±0.05, and a tentative HCN 3nu2-nu2 hot band—are new and will be cited.\n\nThe paper does well in several respects. The authors are transparent about their assumptions: pure absorption, LTE, scalar background fraction f_bg. They flag P/R-branch asymmetries in H2O, CO, and the hot band as evidence for line emission partially filling the bands. They acknowledge the factor 1.6-2.0 discrepancy between 7 and 14 um column densities and offer a plausible explanation involving foreground dust and line emission, but they do not close it with a quantitative model. The toy model in Section 4.1.1 is illustrative; it shows a mechanism, not a fit. No circularity: the independent submillimeter column densities are used only for context.\n\nThe soft spots are real but not fatal. The absolute column densities and the warm HCN decomposition depend on the f_bg = 0.28 pure-absorption assumption; if emission fills in the bands, the columns are lower limits and the warm component's inferred dilution factor changes. The authors say this themselves but do not quantify it. The 7-14 um discrepancy sits at the same place: acknowledged but unresolved. The HCO+ detection is marginal, and the CH4 band sits on the 7.7 um PAH complex, so those numbers carry large systematics. None of this undermines the qualitative picture of hot-core chemistry with shocks and no obvious AGN signature; that conclusion is robust to reasonable changes in the model.\n\nWho is this for? Anyone working on ULIRGs, compact obscured nuclei, or mid-IR molecular absorption as a probe of nuclear ISM. It deserves a serious referee and will get one. My advice: accept with major revisions asking for a sensitivity analysis around the emission-filling and f_bg degeneracy, plus a clearer statement of which columns are robust lower limits. The paper is honest about its limitations, which makes the remaining uncertainties acceptable.","headline":"Genuinely new JWST data and careful fits make this a valuable paper on Arp 220's western nucleus, but the absolute column densities and the warm HCN component rest on pure-absorption LTE assumptions that the authors flag but do not fully quantify.","tokens_in":31487,"tokens_out":2224,"would_cite":true,"duration_ms":21956,"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 spectra show the western nucleus of Arp 220 is a shocked hot core with no AGN signature.","keywords":["Arp 220","ultraluminous infrared galaxies","JWST","mid-infrared spectroscopy","molecular absorption bands","hot cores","starburst galaxy","AGN signatures"],"falsifier":"A resolved, higher-signal observation of the 14 µm HCN ν2 band that shows its Q-branch in emission rather than absorption—or a non-LTE radiative-transfer calculation demonstrating that the P/R-branch asymmetry requires substantial emission filling—would falsify the pure-absorption column densities and the low background fractions.","tokens_in":30320,"feed_emoji":"🔭","tokens_out":7444,"duration_ms":65612,"temperature":0.7,"pith_summary":"This paper analyzes the full 3–28 µm JWST spectra of the western nucleus of Arp 220, the nearest ultraluminous infrared galaxy, and extracts column densities and rotational temperatures for fourteen molecular species from their rovibrational absorption bands. It argues that the absorbing gas is not a cold foreground screen but sits in the inner regions of the 100 pc starburst disk, directly around a hot 20 pc core with dust temperature $T_\\mathrm{d} \\gtrsim 500$ K. The inferred HCN, C$_2$H$_2$, and CO$_2$ column densities are an order of magnitude higher than previous Spitzer estimates, and a warm 330 K HCN component is interpreted as radiatively excited by that hot inner core. The chemical footprint matches Galactic hot cores with an added shock signature, and the paper finds no evidence for X-ray-driven chemistry or extreme excitation that would betray a buried AGN. If correct, the result reshapes how the dense nuclear gas in obscured mergers is probed and what it can say about the nature of the central engine.","feed_headline":"Arp 220's western nucleus shows hot-core chemistry, not an AGN","feed_subtitle":"JWST finds gas columns ten times higher and a 330 K HCN layer around a hot core.","key_machinery":"The load-bearing machinery is the LTE spectral fitting of rovibrational bands with a free scalar background fraction $f_\\mathrm{bg}$, which parametrizes the fraction of observed continuum that passes through the absorbing gas and is constrained when Q-branches saturate. The critical identity is the dilution-factor relation $f_D = (e^{E_{ul}/k_B T_\\mathrm{bg}}-1)/(e^{E_{ul}/k_B T_\\mathrm{rot}}-1)$, which converts the measured 330 K HCN rotational temperature into a geometric dilution factor for a 600–1500 K background, placing the gas at the edge of the opaque core. The supporting toy model—a hot 1000 K background blackbody behind a cool 100 K foreground dust shell—explains the wavelength-dependent background fractions and the partial filling of the 14 µm bands by line emission while leaving the 7 µm bands nearly intact. This combination of saturated-feature fitting, a dilution-factor identity, and a foreground-dilution toy model carries the central argument that the absorbing gas is deep inside the nuclear disk rather than in a cold foreground screen.","core_discovery":"Using MIRI/MRS and NIRSpec/IFU spectra that for the first time separate the two nuclei of Arp 220, the paper detects rovibrational absorption bands of C$_2$H$_2$, HCN, HNC, CO$_2$, H$_2$O, CO, CH$_4$, C$_2$H, NO, N$_2$H$^+$, and other species along the line of sight to the western nucleus. The key quantitative claim is that the optically thick Q-branches of C$_2$H$_2$, HCN, and HNC force a background fraction $f_\\mathrm{bg}$ well below unity (about 0.28 at 14 µm and 0.06 for HNC at 21.6 µm), so the true column densities are much larger than the absorption depths naively suggest: HCN, C$_2$H$_2$, and CO$_2$ columns come out an order of magnitude above the earlier Spitzer values. A second, warm HCN component with $T_\\mathrm{rot}=330$ K is required to fit the 14 µm Q-branch shape; the paper attributes it to radiative excitation by the hot inner nucleus seen with a dilution factor $f_D\\approx0.2$–$0.5$. The derived chemistry is hot-core-like with a low C$_2$H$_2$ ortho-to-para ratio of $1.71\\pm0.05$ and an HCN/HNC ratio above unity, and the paper concludes there is no sign of an AGN in either the chemistry or the excitation.","pith_inferences":["A testable extension is to compare the vibrational temperatures implied by the 7 and 14 µm bands: if the column-density discrepancy is indeed due to emission filling, non-LTE models should predict a measurable 14 µm emission component whose strength correlates with the 330 K rotational temperature.","The same foreground-dilution toy model could be applied to other compact obscured nuclei observed by JWST, predicting a correlation between band wavelength and derived background fraction; a uniform pattern would strengthen the hot-background interpretation.","If the hot core is purely starburst-powered, the inferred H$_2$ column density of $10^{23}$–$10^{24}$ cm$^{-2}$ in the inner disk constrains the star formation rate surface density there and could be compared to the supernova rate traced by the radio VLBI sources.","The tentative HCN 3ν2−ν2 hot-band detection, if confirmed, opens a direct probe of the vibrational temperature and would allow a purely radiative-excitation model of the HCN level populations to be tested against the data."],"forward_implications":["The order-of-magnitude higher column densities mean the nuclear gas mass and abundances derived from mid-infrared absorption in ULIRGs need revision: C$_2$H$_2$ and HCN reach abundances near $10^{-6}$ relative to H$_2$.","The 330 K HCN component is a new diagnostic of the local radiation field and can be used to measure the brightness and dilution of an embedded hot core in other compact obscured nuclei.","The non-detection of the fast molecular outflow in absorption, despite its column density being above the detection limit, implies a covering-factor dilution that must be accounted for in future outflow studies using pencil-beam absorption.","The C$_2$H$_2$ ortho-to-para ratio of $1.71\\pm0.05$, the first extragalactic measurement, provides a new probe of the past thermal history of the gas, analogous to hot-core measurements.","If no AGN is present, the entire energy budget of the western nucleus can be powered by star formation in a compact, highly obscured starburst disk, sharpening the debate on what powers ULIRG nuclei."],"supporting_citations":[{"why":"Provides the earlier Spitzer-based column densities for C2H2, HCN, and CO2 that this paper supersedes by an order of magnitude.","marker":"Lahuis et al. 2007"},{"why":"Identifies the opaque 20 pc core with Td ≳ 500 K that serves as the background source in the proposed geometry.","marker":"Sakamoto et al. 2017"},{"why":"Reports the fast collimated molecular outflow whose non-detection in the MIR absorption is explained by covering-factor dilution.","marker":"Barcos-Muñoz et al. 2018"},{"why":"Constrains the compact emission size and column densities of vibrationally excited HCN, used to model the tentative 3ν2−ν2 hot band.","marker":"Sakamoto et al. 2021b"},{"why":"Provides the Herschel/SPIRE HCN absorption measurements and a 330 K rotational temperature that the warm HCN component is compared to.","marker":"Rangwala et al. 2011"},{"why":"Supplies the hot-core C2H2 ortho-to-para ratios against which the measured 1.71±0.05 value is compared.","marker":"Barr et al. 2020"},{"why":"Demonstrates C-type shock destruction of CO2, supporting the paper's explanation for the low CO2 abundance in the nucleus.","marker":"Charnley & Kaufman 2000"}],"fun_headline_variants":["JWST reveals hot-core chemistry in Arp 220's western nucleus","Arp 220 west: no AGN, just hot-core chemistry","JWST separates nuclei, finds 10x gas columns and 330 K HCN","Arp 220's hidden core: shocked hot cores, no AGN (JWST)","JWST: western nucleus of Arp 220 is hot-core, not AGN"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes every band is a pure absorption line under LTE with a single scalar background fraction $f_\\mathrm{bg}$; if line emission fills in the bands, as the P/R-branch asymmetries and the tentative HCN hot band suggest, then the derived column densities are lower limits and the low derived background fractions change.","fun_headline_variants_meta":{"raw":{"variants":["JWST reveals hot-core chemistry in Arp 220's western nucleus","Arp 220 west: no AGN, just hot-core chemistry","JWST separates nuclei, finds 10x gas columns and 330 K HCN","Arp 220's hidden core: shocked hot cores, no AGN (JWST)","JWST: western nucleus of Arp 220 is hot-core, not AGN"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000838,"raw_usage":{"total_tokens":3819,"prompt_tokens":1277,"completion_tokens":2542,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":893,"completion_tokens_details":{"reasoning_tokens":2436}},"tokens_in":893,"tokens_out":2542,"duration_ms":17418,"temperature":1.0,"reasoning_tokens":2436,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T18:41:23.190739+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A resolved, higher-signal observation of the 14 µm HCN ν2 band that shows its Q-branch in emission rather than absorption—or a non-LTE radiative-transfer calculation demonstrating that the P/R-branch asymmetry requires substantial emission filling—would falsify the pure-absorption column densities and the low background fractions.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the earlier Spitzer-based column densities for C2H2, HCN, and CO2 that this paper supersedes by an order of magnitude."},{"cited_title":"2017, ApJ, 849, 14","cited_arxiv_id":null,"evidence_quote":"Identifies the opaque 20 pc core with Td ≳ 500 K that serves as the background source in the proposed geometry."},{"cited_title":"R., Glenn, J., et al","cited_arxiv_id":null,"evidence_quote":"Provides the Herschel/SPIRE HCN absorption measurements and a 330 K rotational temperature that the warm HCN component is compared to."},{"cited_title":"G., Boogert, A., DeWitt, C","cited_arxiv_id":null,"evidence_quote":"Supplies the hot-core C2H2 ortho-to-para ratios against which the measured 1.71±0.05 value is compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Demonstrates C-type shock destruction of CO2, supporting the paper's explanation for the low CO2 abundance in the nucleus."}],"review_version":1}