We present new C i(1-0) and 12CO(4-3) Plateau de Bure Interferometer observations of five sub-millimeter galaxies (SMGs) and combine these with all available C i(1-0) literature detections in SMGs to ...probe the gas distribution within a sample of 14 systems. We explore the C i(1-0) properties of the SMG population, particularly investigating the ratio of the C i(1-0) luminosity to various 12CO transition and far-infrared luminosities. We find that the SMGs with new observations extend the spread of L
Ci(1-0)/L
FIR to much higher values than found before, with our complete sample providing a good representation of the diverse z > 2 SMG population. We compare the line ratios to the outputs of photodissociation region (PDR) models to constrain the physical conditions in the interstellar medium (ISM) of the SMGs, finding an average density of 〈log (n/cm−3)〉 = 4.3 ± 0.2 and an average radiation field (in terms of the local field value, G
0) of 〈log (G
0)〉 = 3.9 ± 0.4. Overall, we find the SMGs are most comparable to local ultraluminous infrared galaxies (ULIRGs) in G
0 and n; however, a significant tail of 5 of the 14 SMGs are likely best compared to less compact, local starburst galaxies, providing new evidence that many SMGs have extended star formation distributions and are therefore not simply scaled up versions of local ULIRGs. We derive the ISM properties of a sample of quasars also finding that they have higher densities and radiation fields on average than the SMGs, consistent with the more extreme local ULIRGs, and reinforcing their interpretation as transition objects. We explore the limitations of using simple PDR models to understand C i, which may be concomitant with the bulk H2 mass rather than PDR distributed. We therefore also assess C i as a tracer of H2, finding that for our sample SMGs, the H2 masses derived from C i are often consistent with those determined from low excitation 12CO. We conclude that C i observations provide a useful tool to probe the bulk gas and gas processes occurring within merging SMGs, however more detailed, resolved observations are required to fully exploit C i as a diagnostic.
ABSTRACT
Cosmic dust is an essential component shaping both the evolution of galaxies and their observational signatures. How quickly dust builds up in the early Universe remains an open question ...that requires deep observations at (sub-)millimetre wavelengths to resolve. Here, we use Atacama Large Millimeter Array observations of 45 galaxies from the Reionization Era Bright Emission Line Survey (REBELS) and its pilot programs, designed to target C ii and dust emission in UV-selected galaxies at z ∼ 7, to investigate the dust content of high-redshift galaxies through a stacking analysis. We find that the typical fraction of obscured star formation fobs = SFRIR/SFRUV+IR depends on stellar mass, similar to what is observed at lower redshift, and ranges from fobs ≈ 0.3 − 0.6 for galaxies with log10(M⋆/M⊙) = 9.4–10.4. We further adopt the z ∼ 7 stellar mass function from the literature to extract the obscured cosmic star formation rate density (SFRD) from the REBELS survey. Our results suggest only a modest decrease in the SFRD between 3 ≲ z ≲ 7, with dust-obscured star formation still contributing ${\sim}30{{\ \rm per\ cent}}$ at z ∼ 7. While we extensively discuss potential caveats, our analysis highlights the continued importance of dust-obscured star formation even well into the epoch of reionization.
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
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.
Abstract
Gas-rich minor mergers contribute significantly to the gas reservoir of early-type galaxies (ETGs) at low redshift, yet the star formation efficiency (SFE; the star formation rate divided by ...the molecular gas mass) appears to be strongly suppressed following some of these events, in contrast to the more well-known merger-driven starbursts. We present observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of six ETGs, which have each recently undergone a gas-rich minor merger, as evidenced by their disturbed stellar morphologies. These galaxies were selected because they exhibit extremely low SFEs. We use the resolving power of ALMA to study the morphology and kinematics of the molecular gas. The majority of our galaxies exhibit spatial and kinematical irregularities, such as detached gas clouds, warps, and other asymmetries. These asymmetries support the interpretation that the suppression of the SFE is caused by dynamical effects stabilizing the gas against gravitational collapse. Through kinematic modelling we derive high velocity dispersions and Toomre Q stability parameters for the gas, but caution that such measurements in edge-on galaxies suffer from degeneracies. We estimate merger ages to be about 100 Myr based on the observed disturbances in the gas distribution. Furthermore, we determine that these galaxies lie, on average, two orders of magnitude below the Kennicutt–Schmidt relation for star-forming galaxies as well as below the relation for relaxed ETGs. We discuss potential dynamical processes responsible for this strong suppression of star formation surface density at fixed molecular gas surface density.
Abstract
We describe a
Herschel Space Observatory
194–671
μ
m spectroscopic survey of a sample of 121 local luminous infrared galaxies and report the fluxes of the CO
J
to
J
–1 rotational transitions ...for
, the N
ii
205
μ
m line, the C
i
lines at 609 and 370
μ
m, as well as additional and usually fainter lines. The CO spectral line energy distributions (SLEDs) presented here are consistent with our earlier work, which was based on a smaller sample, that calls for two distinct molecular gas components in general: (i) a cold component, which emits CO lines primarily at
J
≲ 4 and likely represents the same gas phase traced by CO (1−0), and (ii) a warm component, which dominates over the mid-
J
regime (4 <
J
≲ 10) and is intimately related to current star formation. We present evidence that the CO line emission associated with an active galactic nucleus is significant only at
J
> 10. The flux ratios of the two C
i
lines imply modest excitation temperatures of 15–30 K; the C
i
370
μ
m line scales more linearly in flux with CO (4−3) than with CO (7−6). These findings suggest that the C
i
emission is predominantly associated with the gas component defined in (i) above. Our analysis of the stacked spectra in different far-infrared (FIR) color bins reveals an evolution of the SLED of the rotational transitions of
vapor as a function of the FIR color in a direction consistent with infrared photon pumping.
Abstract
Feedback and outflows in galaxies that are associated with a quasar phase are expected to be pivotal in quenching the most massive galaxies. However, observations targeting the molecular ...outflow phase, which dominates both the mass and momentum and removes the immediate fuel for star formation, are limited in high-
z
QSO hosts. Massive quiescent galaxies found at
z
∼ 4 are predicted to have quenched star formation already by
z
∼ 5 and undergone their most intense growth at
z
> 6. Here, we present two Atacama Large Millimeter/submillimeter Array (ALMA) detections of molecular outflows, traced by blueshifted absorption of the OH 119
μ
m doublet, from a sample of three
z
> 6 infrared luminous QSO hosts: J2310+1855 and P183+05. OH 119
μ
m is also detected in emission from P183+05, and tentatively in the third source: P036+03. Using similar assumptions as for high-
z
dusty star-forming galaxy outflows, we find that our QSOs drive molecular outflows with comparable mass outflow rates, which are comparably energetic except for J2310+1855's significantly lower outflow energy flux. We do not find evidence, nor require additional input from the central active galactic nucleus (AGN) to drive the molecular outflow in J2310+1855, but we cannot rule out an AGN contribution in P183+05 if a significant AGN contribution to
L
FIR
is assumed and/or if the outflow covering fraction is high (≥53%), which evidence from the literature suggests is unlikely in these sources. Differences observed in the blueshifted absorption spectral properties may instead be caused by the QSO hosts’ more compact dust continuums, limiting observations to lower altitude and more central regions of the outflow.
Abstract
We report the detection of a massive neutral gas outflow in the
z
= 2.09 gravitationally lensed dusty star-forming galaxy HATLAS J085358.9+015537 (G09v1.40), seen in absorption with the OH
+
...(1
1
−1
0
) transition using spatially resolved (0.″5 × 0.″4) Atacama Large Millimeter/submillimeter Array (ALMA) observations. The blueshifted OH
+
line is observed simultaneously with the CO(9−8) emission line and underlying dust continuum. These data are complemented by high-angular-resolution (0.″17 × 0.″13) ALMA observations of CH
+
(1−0) and underlying dust continuum, and Keck 2.2
μ
m imaging tracing the stellar emission. The neutral outflow, dust, dense molecular gas, and stars all show spatial offsets from each other. The total atomic gas mass of the observed outflow is 6.7 × 10
9
M
⊙
, >25% as massive as the gas mass of the galaxy. We find that a conical outflow geometry best describes the OH
+
kinematics and morphology and derive deprojected outflow properties as functions of possible inclination (0.°38–64°). The neutral gas mass outflow rate is between 83 and 25,400
M
⊙
yr
−1
, exceeding the star formation rate (788 ± 300
M
⊙
yr
−1
) if the inclination is >3.°6 (mass-loading factor = 0.3–4.7). Kinetic energy and momentum fluxes span (4.4–290) × 10
9
L
⊙
and (0.1–3.7) × 10
37
dyne, respectively (energy-loading factor = 0.013–16), indicating that the feedback mechanisms required to drive the outflow depend on the inclination assumed. We derive a gas depletion time between 29 and 1 Myr, but find that the neutral outflow is likely to remain bound to the galaxy unless the inclination is small and may be reaccreted if additional feedback processes do not occur.