Starburst Energy Feedback Seen through HCO+/HOC+ Emission in NGC 253 from ALCHEMI

Molecular abundances are sensitive to the UV photon flux and cosmic-ray ionization rate. In starburst environments, the effects of high-energy photons and particles are expected to be stronger. We examine these astrochemical signatures through multiple transitions of HCO+ and its metastable isomer HOC+ in the center of the starburst galaxy NGC 253 using data from the Atacama Large Millimeter/submillimeter Array large program ALMA Comprehensive High-resolution Extragalactic Molecular inventory. The distribution of the HOC+(1−0) integrated intensity shows its association with "superbubbles," cavities created either by supernovae or expanding H ii regions. The observed HCO+/HOC+ abundance ratios are ∼10–150, and the fractional abundance of HOC+ relative to H2 is ∼1.5 × 10−11–6 × 10−10, which implies that the HOC+ abundance in the center of NGC 253 is significantly higher than in quiescent spiral arm dark clouds in the Galaxy and the Galactic center clouds. Comparison with chemical models implies either an interstellar radiation field of G0 ≳ 103 if the maximum visual extinction is ≳5, or a cosmic-ray ionization rate of ζ ≳ 10−14 s−1 (3–4 orders of magnitude higher than that within clouds in the Galactic spiral arms) to reproduce the observed results. From the difference in formation routes of HOC+, we propose that a low-excitation line of HOC+ traces cosmic-ray dominated regions, while high-excitation lines trace photodissociation regions. Our results suggest that the interstellar medium in the center of NGC 253 is significantly affected by energy input from UV photons and cosmic rays, sources of energy feedback

Reconstructing the shock history in the CMZ of NGC 253 with ALCHEMI

Context: HNCO and SiO are well-known shock tracers and have been observed in nearby galaxies, including the nearby (D = 3.5 Mpc) starburst galaxy NGC 253. The simultaneous detection of these two species in regions where the star-formation rate is high may be used to study the shock history of the gas. // Aims: We perform a multi-line molecular study of NGC 253 using the shock tracers SiO and HNCO and aim to characterize its gas properties. We also explore the possibility of reconstructing the shock history in the central molecular zone (CMZ) of the galaxy. // Methods: Six SiO transitions and eleven HNCO transitions were imaged at high resolution 1.″6 (28 pc) with the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the ALCHEMI Large Programme. Both non local thermaldynamic equilibrium (non-LTE) radiative transfer analysis and chemical modeling were performed in order to characterize the gas properties and investigate the chemical origin of the emission. // Results: The nonLTE radiative transfer analysis coupled with Bayesian inference shows clear evidence that the gas traced by SiO has different densities and temperatures than that traced by HNCO, with an indication that shocks are needed to produce both species. Chemical modeling further confirms such a scenario and suggests that fast and slow shocks are responsible for SiO and HNCO production, respectively, in most GMCs. We are also able to infer the physical characteristics of the shocks traced by SiO and HNCO for each GMC. // Conclusions: Radiative transfer and chemical analysis of the SiO and HNCO in the CMZ of NGC 253 reveal a complex picture whereby most of the GMCs are subjected to shocks. We speculate on the possible shock scenarios responsible for the observed emission and provide potential history and timescales for each shock scenario. Observations of higher spatial resolution for these two species are required in order to quantitatively differentiate between the possible scenarios