{"id":"f2556dde-86f4-4879-b8ad-8a303d57daa8","arxiv_id":"2504.19631","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Complex organic molecules in the center of NGC 253 emit from two gas components and most likely trace large-scale shocks, with nitrogen-bearing species depleted toward GMC 7.","lead":"Astronomers mapped five complex organic molecules in the starburst galaxy NGC 253 and found their emission concentrated in the inner 100 parsecs, split between cold gas spread over giant clouds and warmer gas in compact star-forming clumps. The results suggest these molecules, possible building blocks for life, are released by shock waves much as in the center of our own galaxy.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Two-component spatial decomposition rests on rotation-diagram curvature that the paper itself allows could be non-LTE (Sec. 4.2.1); until a non-LTE test is done for CH3CHO/CH3NH2, the GMC-scale shock conclusion and Fig. 6 correlations are not secure.","rationale":"I agree with the reader that the weakest assumption is the attribution of the two rotational-diagram components to GMC and pSSC scales. The strongest single concern is that this attribution is degenerate with non-LTE and optical-depth effects; the paper's own statements in Sec. 4.2.1 and its CH2NH LVG results make this concrete rather than speculative. This does not overturn the paper's cautious qualitative scenario: the moment maps do show emission becoming more compact with increasing Eu, and the CH2NH LVG analysis independently finds a compact source size and a low-Eu line pair that may trace a second component. However, it does mean the quantitative column densities and the chemical correlations in Fig. 6 should not be used as secure constraints until a non-LTE test for the other species is performed. Because the reader's verdict is already CONDITIONAL and identifies essentially this same concern, I recommend UNCHANGED; the non-LTE test described above is the concrete step that would settle whether the two-component spatial interpretation is required. The secondary fragility from having only four correlation points and excluding GMC 5 reinforces, but does not replace, this main concern.","tokens_in":66452,"tokens_out":7943,"duration_ms":86635,"concrete_test":"Run a non-LTE radiative transfer calculation (e.g., RADEX or grelvg) for CH3CHO and CH3NH2 as a single component using the physical conditions derived for CH2NH (Tkin ~ 30–150 K, n(H2) ~ 10^5–10^6 cm^-3, source size ~0.2″, line width ~50 km/s) and column densities near the Table C.1 values, adopting the best available collisional rate coefficients (EMAA/CDMS, or scaled from CH3OH/CH3CN if necessary). Generate synthetic integrated intensities, run them through the same rotation-diagram pipeline, and compare the two-component fit to the single-component non-LTE fit with an information criterion. If the single-component model reproduces the observed high-Eu curvature within the 15% calibration and spectral rms uncertainties, the Sec. 4.2 spatial two-component attribution and the Fig. 6 scale-dependent correlations are not uniquely supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing premise is that the two components fitted in the rotation diagrams correspond to two physical spatial scales: a cold, extended GMC-scale component and a warmer, compact pSSC-scale component (Sec. 4.2). This mapping is insecure for the species that drive the analysis. In Sec. 4.2.1 the authors state that 'the presence of two components is only clear for NH2CHO whilst it is not clear for CH3CHO and CH3NH2,' and for CH3NH2 they allow that 'the deviation in the RD could be due to non-LTE effects.' Nonetheless, CH3NH2 and CH3CHO are fitted with two components, assigned to GMC vs pSSC scales in Table C.1, and those scale-separated column densities are used in the Fig. 6 correlations and in the GMC-scale shock discussion. The only species for which a non-LTE analysis was possible, CH2NH, shows optical depths up to 3.4 and Trot below Tkin (Sec. 4.3), demonstrating that this environment can produce RD curvature from non-LTE/opacity rather than from two distinct spatial components. Because no collisional rates were used for the other species, the two-component-to-two-scale decomposition is degenerate with a single-component non-LTE/opacity model. The qualitative picture may survive, but the quantitative column densities, abundance ratios, and correlation trends built on them are not uniquely determined.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a resolved study of five interstellar complex organic molecules (CH3CHO, C2H5OH, NH2CHO, CH2NH, and CH3NH2) toward the central molecular zone of NGC 253, using 1.6 arcsec (about 27 pc) ALCHEMI ALMA data. The authors image the emission, fit Gaussian line parameters in four regions (GMC 7, GMC 6, pSSC 5, pSSC 2), and derive rotational temperatures and column densities from rotation diagrams under LTE and optically thin assumptions, with beam-filling factors corresponding to GMC scales (1.6 arcsec) and pSSC scales (0.12 arcsec). For CH2NH alone, they additionally perform a non-LTE LVG analysis. The central claim is that the iCOM emission can be reproduced by two gas components: a cold extended component at GMC scales and a warmer compact component at pSSC scales, with the GMC-scale component plausibly tracing large-scale shocks. The paper also reports column-density correlation trends and chemical differences between regions, including an O- versus N-bearing dichotomy toward GMC 7.","tokens_in":66748,"tokens_out":4679,"duration_ms":47064,"significance":"If the central interpretation is correct, this is one of the first resolved extragalactic studies showing that iCOM emission in a starburst CMZ can be decomposed into distinct physical scales and that the chemistry resembles the shock-dominated, sub-thermally excited iCOM chemistry of the Galactic Center. The paper is transparent and carefully hedged in several places: it uses public ALCHEMI data, provides extensive line tables and spectra in appendices, and includes a genuine non-LTE LVG analysis for CH2NH. The observational mapping of five iCOMs across individual GMCs in an external galaxy is a valuable contribution. The main weakness is that the quantitative backbone, especially the two-component-to-two-scale decomposition and the resulting column-density correlations, rests on rotation-diagram curvature that the paper itself acknowledges could be caused by non-LTE or opacity effects for the species driving the analysis.","major_comments":[{"comment":"The load-bearing assumption that the two fitted rotation-diagram components correspond to two distinct physical spatial scales (cold extended GMC-scale gas and warm compact pSSC-scale gas) is not established for CH3CHO and CH3NH2. The paper states in Sec. 4.2.1 that \"the presence of two components is only clear for NH2CHO whilst it is not clear for CH3CHO and CH3NH2,\" and for CH3NH2 that \"the deviation in the RD could be due to non-LTE effects.\" Despite this, both species are fitted with two components, assigned to GMC and pSSC scales in Table C.1, and those scale-separated column densities are used in Fig. 6 and in the shock discussion of Sec. 5.2. Since no non-LTE or opacity treatment is applied to CH3CHO, CH3NH2, NH2CHO, or C2H5OH, the rotation-diagram curvature is degenerate with a single-component non-LTE/opacity model. The paper's own CH2NH LVG analysis (Sec. 4.3) demonstrates that this environment produces optical depths up to 3.4 and rotational temperatures below the kinetic temperature. The two-component-to-two-scale mapping, and every quantitative quantity built on it, is therefore not uniquely constrained; this needs either a non-LTE analysis for the other species or a substantial reframing of the quantitative results as model-dependent.","section":"Sec. 4.2, Sec. 4.2.1, Table C.1"},{"comment":"The chemical correlation analysis and the quoted abundance ratios (e.g., [C2H5OH/CH3CHO] ~ 3 and [CH3NH2/NH2CHO] ~ 3) are computed from column densities whose systematic uncertainties are not propagated. The Ntot values in Table C.1 depend on the assumed source size, beam-filling factor, component assignment, and the assumed 50/50 split of blended E/A and overlapped lines. These choices are correlated across species and regions, so the Pearson coefficients and p-values in Fig. 6 overstate the statistical support for chemical links. The trends may be real, but as presented they do not provide quantitative constraints on formation pathways. The paper should propagate the systematic uncertainties or explicitly present the correlations as illustrative trends rather than measured abundance ratios.","section":"Sec. 5.1, Fig. 6"},{"comment":"For CH2NH, the paper performs an LVG analysis that yields a column density of about 1.5-2 x 10^16 cm^-2 and a compact emission size of 0.15-0.6 arcsec, but then proceeds to use the rotation-diagram column densities of CH2NH in the discussion and in Fig. 6, rather than the LVG-based values. The paper explicitly says the population-diagram correction was only meant to evaluate non-LTE effects, but the RD values differ substantially from the LVG values (e.g., for GMC 6, the RD at pSSC scale gives 4.8 x 10^16 cm^-2 versus the LVG value of 1.5 x 10^16 cm^-2). Because the correlation between CH2NH and CH3NH2 in panels (e) and (j) of Fig. 6 is used as evidence for a chemical link or shared physical process, the choice of CH2NH column density is load-bearing. The authors should either use the LVG-based column densities for CH2NH in the correlation analysis or clearly state that those correlations are based on uncorrected LTE values.","section":"Sec. 4.3, Sec. 5.1.3, Fig. 6"}],"minor_comments":[{"comment":"The caption labels the top panels as \"CH3OH\" but the text in Sec. 4.2.1 and the plotted data refer to CH3CHO; this is likely a typo and should be corrected.","section":"Fig. 3 caption"},{"comment":"The phrase \"we did use them in the analysis\" should read \"we did not use them in the analysis,\" since the text says these blended lines were excluded.","section":"Figs. B.1 and B.10 captions"},{"comment":"The warm component of CH3NH2 toward pSSC 2 has a rotational temperature of 65.8 +/- 37.8 K; this very large uncertainty should be noted in the text when this value is discussed, as it weakens the comparison with other regions.","section":"Table C.1, pSSC 2 CH3NH2"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things before you read it. First, this is the first resolved (~27 pc) per-region excitation analysis of five iCOMs in an extragalactic starburst CMZ, and it includes a non-LTE LVG fit for CH2NH that is new to extragalactic work. Second, its central quantitative claim — a clean split into cold GMC-scale and warm pSSC-scale components — is shakier than the abstract suggests. The authors admit in Sec. 4.2.1 that the two-component decomposition is only clearly supported for NH2CHO, and that the CH3NH2 curvature could be non-LTE. Their own CH2NH LVG finds optical depths up to 3.4 and sub-thermal excitation, which is exactly the sort of effect that can mimic a two-component rotation diagram.\n\nWhat it does well: it is transparent, methodical, and honest. The line lists, tables, extraction details, and caveats are all there. The authors flag the LTE and optically-thin assumptions, the four-region sample, and the shock/heating degeneracy. The CH2NH LVG, while limited by the available collisional rates, is a genuine step beyond rotation diagrams. The emission maps and the FWHM/Vpeak checks are careful.\n\nWhere it is soft: the load-bearing two-component-to-two-scale mapping is not independently tested for CH3CHO, CH3NH2, or C2H5OH. Since the RD curvature can come from opacity or non-LTE excitation, the inferred GMC/pSSC column densities and the Fig. 6 correlation trends built on them are not unique. With only four spatial points, those correlations are suggestive at best. This doesn't sink the qualitative conclusions — sub-thermal excitation, an N-bearing deficiency in GMC 7, and compact CH2NH emission — but it does mean the quantitative results should be read as conditional.\n\nThe paper deserves a serious referee, and I would send it out. The main revision I'd want is a re-analysis or a much more cautious framing of the two-component interpretation, plus a sensitivity test for non-LTE effects on the species without collisional rates. The observational work is valuable and the data are public; a solid revision would make it a reference point for extragalactic iCOM studies.","headline":"First resolved iCOM excitation study in an extragalactic starburst CMZ; the two-component GMC/pSSC split is plausible but not secured against non-LTE, so treat the quantitative column densities as conditional.","tokens_in":67439,"tokens_out":4273,"would_cite":true,"duration_ms":45250,"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":"Five interstellar organic molecules in NGC 253 trace two gas components, one cold and cloud-wide, one warm and compact.","keywords":["NGC 253","central molecular zone","interstellar complex organic molecules","astrochemistry","starburst galaxies","sub-thermal excitation","giant molecular clouds","shock chemistry"],"falsifier":"Resolve the same transitions at sub-arcsecond angular resolution across the four regions: if the high-energy emission is spatially extended rather than concentrated at the pSSC positions, the two-scale decomposition fails; separately, measuring CH3NH2, CH3CHO, and NH2CHO with a non-LTE analysis that includes collisional rates would show whether their rotation-diagram two-component fits are artifacts of optical depth and sub-thermal excitation.","tokens_in":66151,"feed_emoji":"🧪","tokens_out":6917,"duration_ms":69522,"temperature":0.7,"pith_summary":"This paper uses high-resolution millimetre observations of the central molecular zone of the nearby starburst galaxy NGC 253 to map five complex organic molecules—acetaldehyde, ethanol, formamide, methanimine, and methylamine—and asks where their emission comes from and what it traces. The authors find that the emission is concentrated in the inner part of the zone and cannot be explained by a single gas component: a cold, extended component at giant-molecular-cloud scales and a warmer, compact component at the scale of a few parsecs are both needed to reproduce the observed line intensities. If the two-component picture holds, it would mean that resolved extragalactic observations can recover the same shock-dominated, sub-thermally excited organic chemistry seen in the Galactic Center, and that individual molecular clouds in a starburst can differ chemically. The paper also reports chemical differences between the clouds, including a deficit of nitrogen-bearing species toward one cloud and different excitation for two species toward two other regions.","feed_headline":"Two gas components drive organic molecules in NGC 253","feed_subtitle":"Cold extended gas points to large-scale shocks; warmer compact gas marks star-forming clusters.","key_machinery":"The argument is carried by rotational-diagram analysis of the detected transitions, in which line intensities are used to derive rotational temperatures and column densities under LTE and optically thin assumptions, with a beam-filling factor chosen from the assumed emission size. For CH2NH, a non-LTE large-velocity-gradient calculation with published collisional rates replaces the LTE assumption, and it returns a compact emission size of about 2.5–10 pc together with optical depths up to 3.4, showing that the simple assumptions fail for at least one species. Column-density correlation plots between chemically suspected pairs—CH3CHO/C2H5OH, CH3NH2/NH2CHO, and CH3NH2/CH2NH—provide the basis for the formation-path discussion.","core_discovery":"Across the four surveyed regions of the NGC 253 central molecular zone, the emission of CH3CHO, C2H5OH, NH2CHO, CH2NH, and CH3NH2 is shown to be dominated by two distinct gas components. The lower-energy transitions trace cold gas ($T_{\\rm rot}\\lesssim 20$ K) at the scale of giant molecular clouds (about 27 pc), while higher-energy transitions require a warmer component ($T_{\\rm rot}\\sim 30$–80 K) concentrated on scales of a few parsecs, matching the sizes of proto-super-star clusters embedded in the clouds. The authors argue that at cloud scales the widespread, sub-thermally excited emission most plausibly traces large-scale shocks within the clouds, analogous to the Galactic Center, while at the smaller scales both shock and heating from ongoing star formation can act. They also find that column-density correlations support more than one formation path for the organic molecules, and that the four regions are not chemically identical: GMC 7 is depleted in nitrogen-bearing species relative to oxygen-bearing ones, and CH3CHO and NH2CHO show opposite excitation behaviour toward two of the pSSC positions.","pith_inferences":["If the shock interpretation is right, the same sub-thermally excited, shock-dominated iCOM chemistry should be observable toward other starburst nuclear zones with comparable beam sizes, and the N-bearing to O-bearing species ratio may serve as an approximate shock-age indicator.","The paper's two-scale attribution could be tested directly with higher angular resolution: if the warm high-energy component remains spatially extended at sub-arcsecond scales, the assignment of that component to few-parsec proto-clusters would need revision.","The column-density correlations used here could be turned into quantitative formation-path tests by running chemical models that simultaneously predict the ratios [C2H5OH/CH3CHO], [CH3NH2/NH2CHO], and [CH3NH2/CH2NH] under shock and non-shock conditions.","The failure of LTE and optically thin assumptions for CH2NH raises the possibility that the same failure affects the other species; measuring their collisional rates would allow the same non-LTE treatment and would check whether the two-component fits are excitation artifacts."],"forward_implications":["The iCOM emission toward the central molecular zone of NGC 253 is concentrated in the inner ~100 pc around the kinematic centre, between GMC 7 and pSSC 2.","All derived rotational temperatures are below 100 K, mostly 10–40 K, so the molecules are sub-thermally excited rather than tracing hot-core-like gas.","For CH2NH, the non-LTE analysis gives a compact emission size of 0.15–0.6 arcseconds (2.5–10 pc) and gas densities of 10^5–10^6 cm^-3, supporting the picture of emission from shocked, dense gas rather than from hot cores.","At GMC scales, the most favourable explanation is large-scale shocks within the clouds, which would make the organic chemistry of NGC 253 analogous to that of the Galactic Center.","The detected chemical differences between GMCs imply that individual clouds in a starburst nucleus can differ in nitrogen chemistry and excitation, and chemical models will need to reproduce those differences."],"supporting_citations":[{"why":"Provides the survey observations, data cubes, and line identifications on which all measurements rest.","marker":"Martín et al. 2021"},{"why":"Supplies the pSSC positions and sizes (about 0.12 arcseconds) used to set the compact component's beam-filling factor.","marker":"Leroy et al. 2018"},{"why":"Identifies the GMCs, their line widths, and the 35–40 K dust temperatures at cloud scales.","marker":"Leroy et al. 2015"},{"why":"Gives the rotational-diagram method used to derive temperatures and column densities under LTE and optically thin assumptions.","marker":"Goldsmith & Langer 1999"},{"why":"Provides the large-velocity-gradient code used for the non-LTE analysis of CH2NH.","marker":"Ceccarelli et al. 2003"},{"why":"Supplies the CH2NH–H2 collisional rate coefficients that make the non-LTE analysis possible.","marker":"Xue et al. 2024"},{"why":"Provides the HNCO shock tracers and shock timescales used to compare with the iCOM emission.","marker":"Huang et al. 2023"},{"why":"Establishes the Galactic Center template of widespread sub-thermally excited iCOM emission attributed to large-scale shocks.","marker":"Requena-Torres et al. 2006"}],"fun_headline_variants":["Cold and warm gas forge organics in NGC 253","Two gas components cook organic molecules in NGC 253","Shocks and star formation drive NGC 253's organics","NGC 253's organics: two gas, two origins"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the cold low-energy and warm high-energy parts of each molecule's rotation diagram come from two physically separate spatial components—one at the ~27-pc cloud scale and one at the few-parsec pSSC scale—under LTE and optically thin emission; the paper's own non-LTE analysis of CH2NH finds optical depths up to 3.4 and non-LTE excitation, which shows these assumptions do not all hold for at least one species.","fun_headline_variants_meta":{"raw":{"variants":["Cold and warm gas forge organics in NGC 253","Two gas components cook organic molecules in NGC 253","Shocks and star formation drive NGC 253's organics","NGC 253's organics: two gas, two origins"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001246,"raw_usage":{"total_tokens":5224,"prompt_tokens":1171,"completion_tokens":4053,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":787,"completion_tokens_details":{"reasoning_tokens":3985}},"tokens_in":787,"tokens_out":4053,"duration_ms":30022,"temperature":1.0,"reasoning_tokens":3985,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:48:32.924574+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Resolve the same transitions at sub-arcsecond angular resolution across the four regions: if the high-energy emission is spatially extended rather than concentrated at the pSSC positions, the two-scale decomposition fails; separately, measuring CH3NH2, CH3CHO, and NH2CHO with a non-LTE analysis that includes collisional rates would show whether their rotation-diagram two-component fits are artifacts of optical depth and sub-thermal excitation.","supporting_citations":[{"cited_title":"2024, The Astrophysical Journal, 967, 164","cited_arxiv_id":null,"evidence_quote":"Supplies the CH2NH–H2 collisional rate coefficients that make the non-LTE analysis possible."}],"review_version":1}