We present interferometric observations of the CN(1–0) line emission in Mrk 231 and combine them with previous observations of CO and other H
2
gas tracers to study the physical properties of the ...massive molecular outflow. We find a strong boost of the CN/CO(1–0) line luminosity ratio in the outflow of Mrk 231, which is unprecedented compared to any other known Galactic or extragalactic astronomical source. For the dense gas phase in the outflow traced by the HCN and CN emissions, we infer
X
CN
≡ CN/H
2
>
X
HCN
by at least a factor of three, with H
2
gas densities of
n
H
2
∼ 10
5−6
cm
−3
. In addition, we resolve for the first time narrow spectral features in the HCN(1–0) and HCO
+
(1–0) high-velocity line wings tracing the dense phase of the outflow. The velocity dispersions of these spectral features,
σ
v
∼ 7−20 km s
−1
, are consistent with those of massive extragalactic giant molecular clouds detected in nearby starburst nuclei. The H
2
gas masses inferred from the HCN data are quite high,
M
mol
∼ 0.3−5 × 10
8
M
⊙
. Our results suggest that massive complexes of denser molecular gas survive embedded into the more diffuse H
2
phase of the outflow, and that the chemistry of these outflowing dense clouds is strongly affected by UV radiation.
Context. The relationship between outflow launching and the formation of accretion disks around young stellar objects is still not entirely understood, which is why spectrally and spatially resolved ...observations are needed. Recently, the Atacama Large Millimetre/sub-millimetre Array (ALMA) carried out long-baseline observations towards a handful of young sources, revealing connections between outflows and the inner regions of disks. Aims. Here we aim to determine the small-scale kinematical and morphological properties of the outflow from the isolated protostar B335 for which no Keplerian disk has, so far, been observed on scales down to 10 au. Methods. We used ALMA in its longest-baseline configuration to observe emission from CO isotopologues, SiO, SO2, and CH3OH. The proximity of B335 provides a resolution of ~3 au (0.03′′). We also combined our long-baseline data with archival observations to produce a high-fidelity image covering scales up to 700 au (7′′). Results. 12CO has an X-shaped morphology with arms ~50 au in width that we associate with the walls of an outflow cavity, similar to what is observed on larger scales. Long-baseline continuum emission is confined to <7 au from the protostar, while short-baseline continuum emission follows the 12CO outflow and cavity walls. Methanol is detected within ~30 au of the protostar. SiO is also detected in the vicinity of the protostar, but extended along the outflow. Conclusions. The 12CO outflow does not show any clear signs of rotation at distances ≳30 au from the protostar. SiO traces the protostellar jet on small scales, but without obvious rotation. CH3OH and SO2 trace a region <16 au in diameter, centred on the continuum peak, which is clearly rotating. Using episodic, high-velocity, 12CO features, we estimate the launching radius of the outflow to be <0.1 au and dynamical timescales of the order of a few years.
ABSTRACT
We are just starting to understand the physical processes driving the dramatic change in cosmic star formation rate between z ∼ 2 and the present day. A quantity directly linked to star ...formation is the molecular gas density, which should be measured through independent methods to explore variations due to cosmic variance and systematic uncertainties. We use intervening CO absorption lines in the spectra of mm-bright background sources to provide a census of the molecular gas mass density of the Universe. The data used in this work are taken from ALMACAL, a wide and deep survey utilizing the ALMA calibrator archive. While we report multiple Galactic absorption lines and one intrinsic absorber, no extragalactic intervening molecular absorbers are detected. However, due to the large redshift path surveyed (Δz = 182), we provide constraints on the molecular column density distribution function beyond z ∼ 0. In addition, we probe column densities of N(H2) > 1016 atoms cm−2, 5 orders of magnitude lower than in previous studies. We use the cosmological hydrodynamical simulation IllustrisTNG to show that our upper limits of $\rho ({\rm H}_2)\lesssim 10^{8.3}\, \text{M}_{\odot }\, \text{Mpc}^{-3}$ at 0 < z ≤ 1.7 already provide new constraints on current theoretical predictions of the cold molecular phase of the gas. These results are in agreement with recent CO emission-line surveys and are complementary to those studies. The combined constraints indicate that the present decrease of the cosmic star formation rate history is consistent with an increasing depletion of molecular gas in galaxies compared to z ∼ 2.
Context. Argonium has recently been detected as a ubiquitous molecule in our Galaxy. Model calculations indicate that its abundance peaks at molecular fractions in the range of 10-4 to 10-3 and that ...the observed column densities require high values of the cosmic ray ionization rate. Therefore, this molecular cation may serve as an excellent tracer of the very diffuse interstellar medium (ISM), as well as an indicator of the cosmic ray ionization rate. Aims. We attempted to detect ArH+ in extragalactic sources to evaluate its diagnostic power as a tracer of the almost purely atomic ISM in distant galaxies. Methods. We obtained ALMA observations of a foreground galaxy at z = 0.89 in the direction of the lensed blazar PKS 1830−211. Results. Two isotopologs of argonium, 36ArH+ and 38ArH+, were detected in absorption along two different lines of sight toward PKS 1830−211, known as the SW and NE images of the background blazar. The argonium absorption is clearly enhanced on the more diffuse line of sight (NE) compared to other molecular species. The isotopic ratio 36Ar/38Ar is 3.46 ± 0.16 toward the SW image, i.e., significantly lower than the solar value of 5.5. Conclusions. Our results demonstrate the suitability of argonium as a tracer of the almost purely atomic, diffuse ISM in high-redshift sources. The evolution of the isotopic ratio with redshift may help to constrain nucleosynthetic scenarios in the early Universe.
Abstract
We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple ...HCN and HNC rotational transitions observed using the Atacama Large Millimeter/submillimeter Array via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested sampling techniques applied to chemical and radiative transfer models, which are constrained using our HCN and HNC measurements. We find that the volume density
n
H
2
and cosmic-ray ionization rate (CRIR)
ζ
are enhanced by about an order of magnitude in the galaxy’s central regions as compared to those further from the nucleus. In NGC 253's central giant molecular clouds (GMCs), where observed HCN/HNC abundance ratios are the lowest,
n
∼ 10
5.5
cm
−3
and
ζ
∼ 10
−12
s
−1
(greater than 10
4
times the average Galactic rate). We find a positive correlation in the association of both density and CRIR with the number of star formation-related heating sources (supernova remnants, H
ii
regions, and super hot cores) located in each GMC, as well as a correlation between CRIRs and supernova rates. Additionally, we see an anticorrelation between the HCN/HNC ratio and CRIR, indicating that this ratio will be lower in regions where
ζ
is higher. Though previous studies suggested HCN and HNC may reveal strong mechanical heating processes in NGC 253's CMZ, we find cosmic-ray heating dominates the heating budget, and mechanical heating does not play a significant role in the HCN and HNC chemistry.
The cosmic-ray ionization rate (CRIR) is a key parameter in understanding the physical and chemical processes in the interstellar medium. Cosmic rays are a significant source of energy in star ...formation regions, impacting the physical and chemical processes that drive the formation of stars. Previous studies of the circum-molecular zone of the starburst galaxy NGC 253 have found evidence for a high CRIR value: 103–106 times the average CRIR within the Milky Way. This is a broad constraint, and one goal of this study is to determine this value with much higher precision. We exploit ALMA observations toward the central molecular zone of NGC 253 to measure the CRIR. We first demonstrate that the abundance ratio of H3O+ and SO is strongly sensitive to the CRIR. We then combine chemical and radiative transfer models with nested sampling to infer the gas properties and CRIR of several star-forming regions in NGC 253 from emission from their transitions. We find that each of the four regions modeled has a CRIR in the range (1–80) × 10−14 s−1 and that this result adequately fits the abundances of other species that are believed to be sensitive to cosmic rays, including C2H, HCO+, HOC+, and CO. From shock and photon-dominated/X-ray dominated region models, we further find that neither UV-/X-ray-driven nor shock-dominated chemistry is a viable single alternative as none of these processes can adequately fit the abundances of all of these species.
Abstract
The centers of starburst galaxies may be characterized by a specific gas and ice chemistry due to their gas dynamics and the presence of various ice desorption mechanisms. This may result in ...a peculiar observable composition. We analyse the abundances of CO
2
, a reliable tracer of ice chemistry, from data collected as part of the Atacama Large Millimeter/submillimeter Array large program ALCHEMI, a wide-frequency spectral scan toward the starburst galaxy NGC 253 with an angular resolution of 1.″6. We constrain the CO
2
abundances in the gas phase using its protonated form HOCO
+
. The distribution of HOCO
+
is similar to that of methanol, which suggests that HOCO
+
is indeed produced from the protonation of CO
2
sublimated from ice. The HOCO
+
fractional abundances are found to be (1–2) × 10
−9
at the outer part of the central molecular zone (CMZ), while they are lower (∼10
−10
) near the kinematic center. This peak fractional abundance at the outer CMZ is comparable to that in the Milky Way CMZ, and orders of magnitude higher than that in Galactic disk, star-forming regions. From the range of HOCO
+
/CO
2
ratios suggested from chemical models, the gas-phase CO
2
fractional abundance is estimated to be (1–20) × 10
−7
at the outer CMZ, and orders of magnitude lower near the center. We estimate the CO
2
ice fractional abundances at the outer CMZ to be (2–5) × 10
−6
from the literature. A comparison between the ice and gas CO
2
abundances suggests an efficient sublimation mechanism. This sublimation is attributed to large-scale shocks at the orbital intersections of the bar and CMZ.
Abstract
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 H
2
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
G
0
≳ 10
3
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.
Understanding the nature of high-redshift dusty galaxies requires a comprehensive view of their interstellar medium (ISM) and molecular complexity. However, the molecular ISM at high redshifts is ...commonly studied using only a few species beyond
12
C
16
O, limiting our understanding. In this paper, we present the results of deep 3 mm spectral line surveys using the NOrthern Extended Millimeter Array (NOEMA) targeting two strongly lensed dusty galaxies observed when the Universe was less than 1.8 Gyr old: APM 08279+5255, a quasar at redshift
z
= 3.911, and NCv1.143 (
H
-ATLAS J125632.7+233625), a
z
= 3.565 starburst galaxy. The spectral line surveys cover rest-frame frequencies from about 330 to 550 GHz for both galaxies. We report the detection of 38 and 25 emission lines in APM 08279+5255 and NCv1.143, respectively. These lines originate from 17 species, namely CO,
13
CO, C
18
O, CN, CCH, HCN, HCO
+
, HNC, CS, C
34
S, H
2
O, H
3
O
+
, NO, N
2
H
+
, CH, c-C
3
H
2
, and the vibrationally excited HCN and neutral carbon. The spectra reveal the chemical richness and the complexity of the physical properties of the ISM. By comparing the spectra of the two sources and combining the analysis of the molecular gas excitation, we find that the physical properties and the chemical imprints of the ISM are different: the molecular gas is more excited in APM 08279+5255, which exhibits higher molecular gas temperatures and densities compared to NCv1.143; the molecular abundances in APM 08279+5255 are akin to the values of local active galactic nuclei (AGN), showing boosted relative abundances of the dense gas tracers that might be related to high-temperature chemistry and/or the X-ray-dominated regions, while NCv1.143 more closely resembles local starburst galaxies. The most significant differences between the two sources are found in H
2
O: the 448 GHz ortho-H
2
O(4
23
− 3
30
) line is significantly brighter in APM 08279+5255, which is likely linked to the intense far-infrared radiation from the dust powered by AGN. Our astrochemical model suggests that, at such high column densities, far-ultraviolet radiation is less important in regulating the ISM, while cosmic rays (and/or X-rays and shocks) are the key players in shaping the molecular abundances and the initial conditions of star formation. Both our observed CO isotopologs line ratios and the derived extreme ISM conditions (high gas temperatures, densities, and cosmic-ray ionization rates) suggest the presence of a top-heavy stellar initial mass function. From the ∼330–550 GHz continuum, we also find evidence of nonthermal millimeter flux excess in APM 08279+5255 that might be related to the central supermassive black hole. Such deep spectral line surveys open a new window into the physics and chemistry of the ISM and the radiation field of galaxies in the early Universe.
Abstract
We present a spatially resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the Atacama Large Millimeter/submillimeter Array ...Comprehensive High-resolution Extragalactic Molecular Inventory, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way’s Galactic center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3–7 transitions from nine molecular species to delineate the position–position–velocity distributions of column density (
N
H
2
), volume density (
n
H
2
), and temperature (
T
kin
) at 27 pc resolution. Two distinct components are detected: a low-density component with
(
n
H
2
,
T
kin
)
∼
(
10
3.3
cm
−
3
,
85
K
)
and a high-density component with
(
n
H
2
,
T
kin
)
∼
(
10
4.4
cm
−
3
,
110
K
)
, separated at
n
H
2
∼
10
3.8
cm
−
3
. NGC 253 has ∼10 times the high-density gas mass and ∼3 times the dense-gas mass fraction of the GC. These properties are consistent with their HCN/CO ratio but cannot alone explain the factor of ∼30 difference in their star formation efficiencies (SFEs), contradicting the dense-gas mass to star formation rate scaling law. The
n
H
2
histogram toward NGC 253 exhibits a shallow declining slope up to
n
H
2
∼
10
6
cm
−
3
, while that of the GC steeply drops in
n
H
2
≳
10
4.5
cm
−
3
and vanishes at 10
5
cm
−3
. Their dense-gas mass fraction ratio becomes consistent with their SFEs when the threshold
n
H
2
for the dense gas is taken at ∼10
4.2−4.6
cm
−3
. The rich abundance of gas above this density range in the NGC 253 CMZ, or its scarcity in the GC, is likely to be the critical difference characterizing the contrasting star formation in the centers of the two galaxies.