We present new measurements of the cosmic cold molecular gas evolution out to redshift 6 based on systematic mining of the Atacama Large Millimeter/submillimeter Array (ALMA) public archive in the ...COSMOS deep field (A3COSMOS). Our A3COSMOS data set contains ∼700 galaxies (0.3 z 6) with high-confidence ALMA detections in the (sub)millimeter continuum and multiwavelength spectral energy distributions. Multiple gas mass calibration methods are compared, and biases in band conversions (from observed ALMA wavelength to rest-frame Rayleigh-Jeans tail continuum) have been tested. Combining our A3COSMOS sample with ∼1000 CO-observed galaxies at 0 z 4 (75% at z < 0.1), we parameterize galaxies' molecular gas depletion time ( ) and molecular gas to stellar mass ratio ( ) each as a function of the stellar mass ( ), offset from the star-forming main sequence ( ) and cosmic age (or redshift). Our proposed functional form provides a statistically better fit to current data (than functional forms in the literature) and implies a "downsizing" effect (i.e., more-massive galaxies evolve earlier than less-massive ones) and "mass quenching" (gas consumption slows down with cosmic time for massive galaxies but speeds up for low-mass ones). Adopting galaxy stellar mass functions and applying our function for gas mass calculation, we for the first time infer the cosmic cold molecular gas density evolution out to redshift 6 and find agreement with CO blind surveys as well as semianalytic modeling. These together provide a coherent picture of cold molecular gas, star formation rate, and stellar mass evolution in galaxies across cosmic time.
The existence of massive (10
solar masses) elliptical galaxies by redshift z ≈ 4 (refs 1, 2, 3; when the Universe was 1.5 billion years old) necessitates the presence of galaxies with star-formation ...rates exceeding 100 solar masses per year at z > 6 (corresponding to an age of the Universe of less than 1 billion years). Surveys have discovered hundreds of galaxies at these early cosmic epochs, but their star-formation rates are more than an order of magnitude lower. The only known galaxies with very high star-formation rates at z > 6 are, with one exception, the host galaxies of quasars, but these galaxies also host accreting supermassive (more than 10
solar masses) black holes, which probably affect the properties of the galaxies. Here we report observations of an emission line of singly ionized carbon (C ii at a wavelength of 158 micrometres) in four galaxies at z > 6 that are companions of quasars, with velocity offsets of less than 600 kilometres per second and linear offsets of less than 100 kiloparsecs. The discovery of these four galaxies was serendipitous; they are close to their companion quasars and appear bright in the far-infrared. On the basis of the C ii measurements, we estimate star-formation rates in the companions of more than 100 solar masses per year. These sources are similar to the host galaxies of the quasars in C ii brightness, linewidth and implied dynamical mass, but do not show evidence for accreting supermassive black holes. Similar systems have previously been found at lower redshift. We find such close companions in four out of the twenty-five z > 6 quasars surveyed, a fraction that needs to be accounted for in simulations. If they are representative of the bright end of the C ii luminosity function, then they can account for the population of massive elliptical galaxies at z ≈ 4 in terms of the density of cosmic space.
ABSTRACT We present high-resolution (0 16) 870 m Atacama Large Millimeter/submillimeter Array (ALMA) imaging of 16 luminous ( ) submillimeter galaxies (SMGs) from the ALESS survey of the Extended ...Chandra Deep Field South. This dust imaging traces the dust-obscured star formation in these galaxies on ∼1.3 kpc scales. The emission has a median effective radius of Re = 0 24 0 02, corresponding to a typical physical size of 1.8 0.2 kpc. We derive a median Sérsic index of n = 0.9 0.2, implying that the dust emission is remarkably disk-like at the current resolution and sensitivity. We use different weighting schemes with the visibilities to search for clumps on 0 12 (∼1.0 kpc) scales, but we find no significant evidence for clumping in the majority of cases. Indeed, we demonstrate using simulations that the observed morphologies are generally consistent with smooth exponential disks, suggesting that caution should be exercised when identifying candidate clumps in even moderate signal-to-noise ratio interferometric data. We compare our maps to comparable-resolution Hubble Space Telescope -band images, finding that the stellar morphologies appear significantly more extended and disturbed, and suggesting that major mergers may be responsible for driving the formation of the compact dust disks we observe. The stark contrast between the obscured and unobscured morphologies may also have implications for SED fitting routines that assume the dust is co-located with the optical/near-IR continuum emission. Finally, we discuss the potential of the current bursts of star formation to transform the observed galaxy sizes and light profiles, showing that the descendants of these SMGs are expected to have stellar masses, effective radii, and gas surface densities consistent with the most compact massive ( 1-2 × 1011 ) early-type galaxies observed locally.
We present deep, high-resolution (0 03, 200 pc) ALMA Band 7 observations covering the dust continuum and C ii λ157.7 m emission in four z ∼ 4.4-4.8 sub-millimeter galaxies (SMGs) selected from the ...ALESS and AS2UDS surveys. The data show that the rest-frame 160 m (observed 345 GHz) dust emission is consistent with smooth morphologies on kpc scales for three of the sources. One source, UDS 47.0, displays apparent substructure, but this is also consistent with a smooth morphology-as indicated by simulations showing that smooth exponential disks can appear clumpy when observed at the high angular resolution (0 03) and depth of these observations ( Jy beam−1). The four SMGs are bright C ii emitters. We extract C ii spectra from the high-resolution data, and recover ∼20%-100% of the C ii flux and ∼40%-80% of the dust continuum emission, compared to the previous lower-resolution observations. When tapered to 0 2 resolution, our maps recover ∼80%-100% of the continuum emission, indicating that ∼60% of the emission is resolved out on ∼200 pc scales. We find that the C ii emission in high-redshift galaxies is more spatially extended than the rest-frame 160 m dust continuum by a factor of 1.6 0.4. By considering the / ratio as a function of the star formation rate surface density ( ), we revisit the C ii deficit and suggest that the decline in the / ratio as a function of is consistent with local processes. We also explore the physical drivers that may be responsible for these trends and can give rise to the properties found in the densest regions of SMGs.
Dusty star-forming galaxies at high redshift (1 < z < 3) represent the most intense star-forming regions in the universe. Key aspects to these processes are the gas heating and cooling mechanisms, ...and although it is well known that these galaxies are gas-rich, little is known about the gas excitation conditions. Only a few detailed radiative transfer studies have been carried out owing to a lack of multiple line detections per galaxy. Here we examine these processes in a sample of 24 strongly lensed star-forming galaxies identified by the Planck satellite (LPs) at z ∼ 1.1-3.5. We analyze 162 CO rotational transitions (ranging from Jup = 1 to 12) and 37 atomic carbon fine-structure lines (C i) in order to characterize the physical conditions of the gas in the sample of LPs. We simultaneously fit the CO and C i lines and the dust continuum emission, using two different non-LTE, radiative transfer models. The first model represents a two-component gas density, while the second assumes a turbulence-driven lognormal gas density distribution. These LPs are among the most gas-rich, IR-luminous galaxies ever observed ( L L IR ( 8 − 1000 m ) ∼ 10 13 − 14.6 L ; 〈 LMISM 〉 = (2.7 1.2) × 1012 M , with L ∼ 10-30 the average lens magnification factor). Our results suggest that the turbulent interstellar medium present in the LPs can be well characterized by a high turbulent velocity dispersion ( 〈 ΔVturb 〉 ∼ 100 km s−1) and ratios of gas kinetic temperature to dust temperature 〈 Tkin/Td 〉 ∼ 2.5, sustained on scales larger than a few kiloparsecs. We speculate that the average surface density of the molecular gas mass and IR luminosity, M ISM ∼ 103-4 M pc−2 and L IR ∼ 1011-12 L kpc−2, arise from both stellar mechanical feedback and a steady momentum injection from the accretion of intergalactic gas.
ABSTRACT We present an analysis of a deep (1 = 13 Jy) cosmological 1.2 mm continuum map based on ASPECS, the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field. In the 1 arcmin2 covered by ...ASPECS we detect nine sources at significance at 1.2 mm. Our ALMA-selected sample has a median redshift of , with only one galaxy detected at z > 2 within the survey area. This value is significantly lower than that found in millimeter samples selected at a higher flux density cutoff and similar frequencies. Most galaxies have specific star formation rates (SFRs) similar to that of main-sequence galaxies at the same epoch, and we find median values of stellar mass and SFRs of and yr−1, respectively. Using the dust emission as a tracer for the interstellar medium (ISM) mass, we derive depletion times that are typically longer than 300 Myr, and we find molecular gas fractions ranging from ∼0.1 to 1.0. As noted by previous studies, these values are lower than those using CO-based ISM estimates by a factor of ∼2. The 1 mm number counts (corrected for fidelity and completeness) are in agreement with previous studies that were typically restricted to brighter sources. With our individual detections only, we recover 55% 4% of the extragalactic background light (EBL) at 1.2 mm measured by the Planck satellite, and we recover 80% 7% of this EBL if we include the bright end of the number counts and additional detections from stacking. The stacked contribution is dominated by galaxies at , with stellar masses of (1-3) × 1010 M . For the first time, we are able to characterize the population of galaxies that dominate the EBL at 1.2 mm.
We present a multiline survey of the interstellar medium (ISM) in two
z
> 6 quasar host galaxies, PJ231−20 (
z
= 6.59) and PJ308−21 (
z
= 6.23), and their two companion galaxies. Observations were ...carried out using the Atacama Large (sub-)Millimeter Array (ALMA). We targeted 11 transitions including atomic fine-structure lines (FSLs) and molecular lines: NII
205 μm
, CI
369 μm
, CO (
J
up
= 7, 10, 15, 16), H
2
O 3
12
− 2
21
, 3
21
− 3
12
, 3
03
− 2
12
, and the OH
163 μm
doublet. The underlying far-infrared (FIR) continuum samples the Rayleigh-Jeans tail of the respective dust emission. By combining this information with our earlier ALMA CII
158 μm
observations, we explored the effects of star formation and black hole feedback on the ISM of the galaxies using the CLOUDY radiative transfer models. We estimated dust masses, spectral indexes, IR luminosities, and star-formation rates from the FIR continuum. The analysis of the FSLs indicates that the CII
158 μm
and CI
369 μm
emission arises predominantly from the neutral medium in photodissociation regions (PDRs). We find that line deficits agree with those of local luminous IR galaxies. The CO spectral line energy distributions (SLEDs) reveal significant high-
J
CO excitation in both quasar hosts. Our CO SLED modeling of the quasar PJ231−20 shows that PDRs dominate the molecular mass and CO luminosities for
J
up
≤ 7, while the
J
up
≥ 10 CO emission is likely driven by X-ray dissociation regions produced by the active galactic nucleus (AGN) at the very center of the quasar host. The
J
up
> 10 lines are undetected in the other galaxies in our study. The H
2
O 3
21
− 3
12
line detection in the same quasar places this object on the
L
H
2
O
−
L
TIR
relation found for low-
z
sources, thus suggesting that this water vapor transition is predominantly excited by IR pumping. Models of the H
2
O SLED and of the H
2
O-to-OH
163 μm
ratio point to PDR contributions with high volume and column density (
n
H
∼ 0.8 × 10
5
cm
−3
,
N
H
= 10
24
cm
−2
) in an intense radiation field. Our analysis suggests a less highly excited medium in the companion galaxies. However, the current data do not allow us to definitively rule out an AGN in these sources, as suggested by previous studies of the same objects. This work demonstrates the power of multiline studies of FIR diagnostics in order to dissect the physical conditions in the first massive galaxies emerging from cosmic dawn.
Aims. To study the structure of nearby (<500 pc) dense starless and star-forming cores with the particular goal to identify and understand evolutionary trends in core properties, and to explore the ...nature of Very Low Luminosity Objects (≤0.1 $L_{\odot}$; VeLLOs). Methods. Using the MAMBO bolometer array, we create maps unusually sensitive to faint (few mJy per beam) extended (≈5$ \arcmin$) thermal dust continuum emission at 1.2 mm wavelength. Complementary information on embedded stars is obtained from Spitzer, IRAS, and 2MASS. Results. Our maps are very rich in structure, and we characterize extended emission features (“subcores”) and compact intensity peaks in our data separately to pay attention to this complexity. We derive, e.g., sizes, masses, and aspect ratios for the subcores, as well as column densities and related properties for the peaks. Combination with archival infrared data then enables the derivation of bolometric luminosities and temperatures, as well as envelope masses, for the young embedded stars. Conclusions. Starless and star-forming cores occupy the same parameter space in many core properties; a picture of dense core evolution in which any dense core begins to actively form stars once it exceeds some fixed limit in, e.g., mass, density, or both, is inconsistent with our data. A concept of necessary conditions for star formation appears to provide a better description: dense cores fulfilling certain conditions can form stars, but they do not need to, respectively have not done so yet. Comparison of various evolutionary indicators for young stellar objects in our sample (e.g., bolometric temperatures) reveals inconsistencies between some of them, possibly suggesting a revision of some of these indicators. Finally, we challenge the notion that VeLLOs form in cores not expected to actively form stars, and we present a first systematic study revealing evidence for structural differences between starless and candidate VeLLO cores.
Dust-enshrouded, starbursting, submillimeter galaxies (SMGs) at z ≥ 3 have been proposed as progenitors of z ≥ 2 compact quiescent galaxies (cQGs). To test this connection, we present a detailed ...spatially resolved study of the stars, dust, and stellar mass in a sample of six submillimeter-bright starburst galaxies at z ∼ 4.5. The stellar UV emission probed by HST is extended and irregular and shows evidence of multiple components. Informed by HST, we deblend Spitzer/IRAC data at rest-frame optical, finding that the systems are undergoing minor mergers with a typical stellar mass ratio of 1:6.5. The FIR dust continuum emission traced by ALMA locates the bulk of star formation in extremely compact regions (median re = 0.70 0.29 kpc), and it is in all cases associated with the most massive component of the mergers (median ). We compare spatially resolved UV slope (β) maps with the FIR dust continuum to study the infrared excess (IRX = LIR/LUV)-β relation. The SMGs display systematically higher IRX values than expected from the nominal trend, demonstrating that the FIR and UV emissions are spatially disconnected. Finally, we show that the SMGs fall on the mass-size plane at smaller stellar masses and sizes than the cQGs at z = 2. Taking into account the expected evolution in stellar mass and size between z = 4.5 and z = 2 due to the ongoing starburst and mergers with minor companions, this is in agreement with a direct evolutionary connection between the two populations.
ABSTRACT We present a search for C ii line and dust continuum emission from optical dropout galaxies at z > 6 using ASPECS, our Atacama Large Millimeter submillimeter Array Spectroscopic Survey in ...the Hubble Ultra-deep Field (UDF). Our observations, which cover the frequency range of 212-272 GHz, encompass approximately the range of 6 < z < 8 for C ii line emission and reach a limiting luminosity of LC ii ∼ (1.6-2.5) × 108 L . We identify 14 C ii line emitting candidates in this redshift range with significances >4.5 , two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60% of our candidates are expected to be spurious. For one of our blindly selected C ii line candidates, we tentatively detect the CO(6-5) line in our parallel 3 mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all C ii candidates results in a tentative detection with S1.2 mm = 14 5 Jy. This implies a dust-obscured star-formation rate (SFR) of (3 1) M yr−1. We find that the two highest-SFR objects have candidate C ii lines with luminosities that are consistent with the low-redshift LC ii versus SFR relation. The other candidates have significantly higher C ii luminosities than expected from their UV-based SFR. At the current sensitivity, it is unclear whether the majority of these sources are intrinsically bright C ii emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for C ii emitters at 6 < z < 8 that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.