Astronomy: Starbursts near and far

Observations of intensely bright star-forming galaxies both close by and in the distant Universe at first glance seem to emphasize their similarity. But look a little closer, and differences emerge.Comment: 6 pages including 1 figur

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

Kinetic temperature of massive star-forming molecular clumps measured with formaldehyde IV. The ALMA view of N113 and N159W in the LMC

We mapped the kinetic temperature structure of two massive star-forming regions, N113 and N159W, in the Large Magellanic Cloud (LMC). We have used ~1.′′6 (~0.4 pc) resolution measurements of the para-H2CO JKaKc = 303–202, 322–221, and 321–220 transitions near 218.5 GHz to constrain RADEX non local thermodynamic equilibrium models of the physical conditions. The gas kinetic temperatures derived from the para-H2CO line ratios 322–221/303–202 and 321–220/303–202 range from 28 to 105 K in N113 and 29 to 68 K in N159W. Distributions of the dense gas traced by para-H2CO agree with those of the 1.3 mm dust and Spitzer 8.0 μm emission, but they do not significantly correlate with the Hα emission. The high kinetic temperatures (Tkin ≳ 50 K) of the dense gas traced by para-H2CO appear to be correlated with the embedded infrared sources inside the clouds and/or young stellar objects in the N113 and N159W regions. The lower temperatures (Tkin < 50 K) were measured at the outskirts of the H2CO-bearing distributions of both N113 and N159W. It seems that the kinetic temperatures of the dense gas traced by para-H2CO are weakly affected by the external sources of the Hα emission. The non thermal velocity dispersions of para-H2CO are well correlated with the gas kinetic temperatures in the N113 region, implying that the higher kinetic temperature traced by para-H2CO is related to turbulence on a ~0.4 pc scale. The dense gas heating appears to be dominated by internal star formation activity, radiation, and/or turbulence. It seems that the mechanism heating the dense gas of the star-forming regions in the LMC is consistent with that in Galactic massive star-forming regions located in the Galactic plane

Atmospheric Refractive Electromagnetic Wave Bending and Propagation Delay

In this tutorial we summarize the physics and mathematics behind refractive electromagnetic wave bending and delay. Refractive bending and delay through the Earth's atmosphere at both radio/millimetric and optical/IR wavelengths are discussed, but with most emphasis on the former, and with Atacama Large Millimeter Array (ALMA) applications in mind. As modern astronomical measurements often require sub-arcsecond position accuracy, care is required when selecting refractive bending and delay algorithms. For the spherically-uniform model atmospheres generally used for all refractive bending and delay algorithms, positional accuracies $\lesssim 1^{\prime\prime}$ are achievable when observing at zenith angles $\lesssim 75^\circ$. A number of computationally economical approximate methods for atmospheric refractive bending and delay calculation are presented, appropriate for astronomical observations under these conditions. For observations under more realistic atmospheric conditions, for zenith angles $\gtrsim 75^\circ$, or when higher positional accuracy is required, more rigorous refractive bending and delay algorithms must be employed. For accurate calculation of the refractive bending, we recommend the Auer & Standish (2000) method, using numerical integration to ray-trace through a two-layer model atmosphere, with an atmospheric model determination of the atmospheric refractivity. For the delay calculation we recommend numerical integration through a model atmosphere.Comment: 18 pages with 8 figures. Published 2015 PASP 127, p. 74-91. Found error in code which produced Figure 4 and some of the associated values listed in section 3.1.4, but did not affect any other calculations. This version includes replaced text and figure provided as erratum which will appear in PASP May 201