No Arabic abstract
Observations of chemical species can provide an insight into the physical conditions of the emitting gas but it is important to understand how their abundances and excitation vary within different heating environments. C$_2$H is a molecule typically found in PDR regions of our own Galaxy but there is evidence to suggest it also traces other regions undergoing energetic processing in extragalactic environments. As part of the ALCHEMI ALMA large program, the emission of C$_2$H in the central molecular zone of the nearby starburst galaxy NGC 253 was mapped at 1.6 (28 pc) resolution and characterized to understand its chemical origins. Spectral modelling of the N=1-0 through N=4-3 rotational transitions of C$_2$H was used to derive the C$_2$H column densities towards the dense clouds in NGC 253. Chemical modelling, including PDR, dense cloud, and shock models were then used to investigate the chemical processes and physical conditions that are producing the molecular emission. We find high C$_2$H column densities of $sim 10^{15} cm^{-3}$ detected towards the dense regions of NGC 253. We further find that these column densities cannot be reproduced by assuming that the emission arises from the PDR regions at the edge of the clouds. Instead, we find that the C$_2$H abundance remains high even in the high visual extinction interior of these clouds and that this is most likely caused by a high cosmic-ray ionization rate.
Molecular abundances are sensitive to 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 ALMA large program ALCHEMI. The distribution of the HOC$^+$(1-0) integrated intensity shows its association with superbubbles, cavities created either by supernovae or expanding HII regions. The observed HCO$^+$/HOC$^+$ abundance ratios are $sim 10-150$, and the fractional abundance of HOC$^+$ relative to H$_2$ is $sim 1.5times 10^{-11} - 6times 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 $G_0gtrsim 10^3$ if the maximum visual extinction is $gtrsim 5$, or a cosmic-ray ionization rate of $zeta gtrsim 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.
We present new radio continuum observations of NGC253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of central starburst and extended emission. The central component, corresponding to the inner 500pc of the starburst region of the galaxy, is best modelled as an internally free-free absorbed synchrotron plasma, with a turnover frequency around 230 MHz. The extended emission component of the NGC253 spectrum is best described as a synchrotron emission flattening at low radio frequencies. We find that 34% of the extended emission (outside the central starburst region) at 1 GHz becomes partially absorbed at low radio frequencies. Most of this flattening occurs in the western region of the SE halo, and may be indicative of synchrotron self-absorption of shock re-accelerated electrons or an intrinsic low-energy cut off of the electron distribution. Furthermore, we detect the large-scale synchrotron radio halo of NGC253 in our radio images. At 154 - 231 MHz the halo displays the well known X-shaped/horn-like structure, and extends out to ~8kpc in z-direction (from major axis).
We study the metallicity dependence of the H/H$_2$ and C$^+$/C/CO distributions in a self-regulated interstellar medium (ISM) across a broad range of metallicities ($0.1 < Z/Z_odot < 3$). To this end, we conduct high-resolution (particle mass of $1 {rm M_odot}$) hydrodynamical simulations coupled with a time-dependent H$_2$ chemistry network. The results are then post-processed with an accurate chemistry network to model the associated C$^+$/C/CO abundances, based on the time-dependent non-steady-state (``non-equilibrium) H$_2$ abundances. We find that the time-averaged star formation rate and the ISM structure are insensitive to metallicity. The column densities relevant for molecular shielding appear correlated with the volume densities in gravitationally unstable gas. As metallicity decreases, H$_2$ progressively deviates from steady state (``equilibrium) and shows shallow abundance profiles until they sharply truncate at the photodissociation fronts. In contrast, the CO profile is sharp and controlled by photodissociation as CO quickly reaches steady state. We construct effective one-dimensional cloud models that successfully capture the time-averaged chemical distributions in simulations. At low metallicities, the steady-state model significantly overestimates the abundance of H$_2$ in the diffuse medium. The overestimated H$_2$, however, has little impact on CO. Consequently, the mass fraction of CO-dark H$_2$ gas is significantly lower than what a fully steady-state model predicts. The mass ratios of H$_2$/C$^+$ and H$_2$/C both show a weaker dependence on $Z^{prime}$ than H$_2$/CO, which potentially indicates that C$^+$ and C could be alternative tracers for H$_2$ at low $Z^{prime}$ in terms of mass budget. Our chemistry code for post-processing is publicly available.
This work presents swarm parameters of electrons (the bulk drift velocity, the bulk longitudinal component of the diffusion tensor, and the effective ionization frequency) in C$_2$H$_n$, with $n =$ 2, 4 and 6, measured in a scanning drift tube apparatus under time-of-flight conditions over a wide range of the reduced electric field, 1 Td $leq,E/N,leq$ 1790 Td (1 Td = $10^{-21}$ Vm$^2$). The effective steady-state Townsend ionization coefficient is also derived from the experimental data. A kinetic simulation of the experimental drift cell allows estimating the uncertainties introduced in the data acquisition procedure and provides a correction factor to each of the measured swarm parameters. These parameters are compared to results of previous experimental studies, as well as to results of various kinetic swarm calculations: solutions of the electron Boltzmann equation under different approximations (multiterm and density gradient expansions) and Monte Carlo simulations. The experimental data are consistent with most of the swarm parameters obtained in earlier studies. In the case of C$_2$H$_2$, the swarm calculations show that the thermally excited vibrational population should not be neglected, in particular, in the fitting of cross sections to swarm results.
H$_2$CO ice on dust grains is an important precursor of complex organic molecules (COMs). H$_2$CO gas can be readily observed in protoplanetary disks and may be used to trace COM chemistry. However, its utility as a COM probe is currently limited by a lack of constraints on the relative contributions of two different formation pathways: on icy grain-surfaces and in the gas-phase. We use archival ALMA observations of the resolved distribution of H$_2$CO emission in the disk around the young low-mass star DM Tau to assess the relative importance of these formation routes. The observed H$_2$CO emission has a centrally peaked and radially broad brightness profile (extending out to 500 AU). We compare these observations with disk chemistry models with and without grain-surface formation reactions, and find that both gas and grain-surface chemistry are necessary to explain the spatial distribution of the emission. Gas-phase H$_2$CO production is responsible for the observed central peak, while grain-surface chemistry is required to reproduce the emission exterior to the CO snowline (where H$_2$CO mainly forms through the hydrogenation of CO ice before being non-thermally desorbed). These observations demonstrate that both gas and grain-surface pathways contribute to the observed H$_2$CO in disks, and that their relative contributions depend strongly on distance from the host star.