We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift ...evolution of (H2). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at z > 1. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from ∼2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to ∼90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at z = 1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H2 conversion factor, we constrain the evolution of the cosmic molecular gas density (H2): this cosmic gas density peaks at z ∼ 1.5 and drops by a factor of to the value measured locally. The observed evolution in (H2), therefore, closely matches the evolution of the cosmic star formation rate density SFR. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of (H2) and SFR leaves only little room for a significant evolution of the average star formation efficiency in galaxies since z ∼ 3 (85% of cosmic history).
We analyze the interstellar medium properties of a sample of 16 bright CO line emitting galaxies identified in the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Large Program. ...This CO−selected galaxy sample is complemented by two additional CO line emitters in the UDF that are identified based on their Multi-Unit Spectroscopic Explorer (MUSE) optical spectroscopic redshifts. The ASPECS CO−selected galaxies cover a larger range of star formation rates (SFRs) and stellar masses compared to literature CO emitting galaxies at z > 1 for which scaling relations have been established previously. Most of ASPECS CO-selected galaxies follow these established relations in terms of gas depletion timescales and gas fractions as a function of redshift, as well as the SFR-stellar mass relation ("galaxy main sequence"). However, we find that ∼30% of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with ∼12% (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low SFRs, and yet show substantial molecular gas reservoirs, yielding long gas depletion timescales. Capitalizing on the well-defined cosmic volume probed by our observations, we measure the contribution of galaxies above, below, and on the galaxy main sequence to the total cosmic molecular gas density at different lookback times. We conclude that main-sequence galaxies are the largest contributors to the molecular gas density at any redshift probed by our observations (z ∼ 1−3). The respective contribution by starburst galaxies above the main sequence decreases from z ∼ 2.5 to z ∼ 1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequence.
The Atacama Large Millimeter/submillimeter Array (ALMA) SPECtroscopic Survey in the Hubble Ultra Deep Field (HUDF) is an ALMA large program that obtained a frequency scan in the 3 mm band to detect ...emission lines from the molecular gas in distant galaxies. Here we present our search strategy for emission lines and continuum sources in the HUDF. We compare several line search algorithms used in the literature, and critically account for the line widths of the emission line candidates when assessing significance. We identify 16 emission lines at high fidelity in our search. Comparing these sources to multiwavelength data we find that all sources have optical/infrared counterparts. Our search also recovers candidates of lower significance that can be used statistically to derive, e.g., the CO luminosity function. We apply the same detection algorithm to obtain a sample of six 3 mm continuum sources. All of these are also detected in the 1.2 mm continuum with optical/near-infrared counterparts. We use the continuum sources to compute 3 mm number counts in the sub-millijansky regime, and find them to be higher by an order of magnitude than expected for synchrotron-dominated sources. However, the number counts are consistent with those derived at shorter wavelengths (0.85-1.3 mm) once extrapolating to 3 mm with a dust emissivity index of β = 1.5, dust temperature of 35 K, and an average redshift of z = 2.5. These results represent the best constraints to date on the faint end of the 3 mm number counts.
We present a characterization of the physical properties of a sample of 35 securely detected, dusty galaxies in the deep ALMA 1.2 mm image obtained as part of the ALMA Spectroscopic Survey in the ...Hubble Ultra Deep Field (ASPECS) Large Program. This sample is complemented by 26 additional sources identified via an optical/infrared source positional prior. Using their well-characterized spectral energy distributions, we derive median stellar masses and star formation rates (SFR) of and 30 M☉ yr−1, respectively, and interquartile ranges of (2.4-11.7) × 1010 M☉ and 20-50 M☉ yr−1. We derive a median spectroscopic redshift of 1.8 with an interquartile range 1.1-2.6, significantly lower than submillimeter galaxies detected in shallower, wide-field surveys. We find that 59% 13%, 6% 4%, and 34% 9% of our sources are within, above, and below 0.4 dex from the SFR-stellar-mass relation or main sequence (MS), respectively. The ASPECS galaxies closely follow the SFR-molecular gas mass relation and other previously established scaling relations, confirming a factor of five increase of the gas-to-stellar-mass ratio from z = 0.5 to 2.5 and a mild evolution of the gas depletion timescales with a typical value of 0.7 Gyr at z = 1-3. ASPECS galaxies located significantly below the MS, a poorly exploited parameter space, have low gas-to-stellar-mass ratios of ∼0.1-0.2 and long depletion timescales >1 Gyr. Galaxies along the MS dominate the cosmic density of molecular gas at all redshifts. Systems above the MS have an increasing contribution to the total gas reservoirs from z < 1 to z = 2.5, while the opposite is found for galaxies below the MS.
We discuss the nature and physical properties of gas-mass selected galaxies in the ALMA spectroscopic survey (ASPECS) of the Hubble Ultra Deep Field (HUDF). We capitalize on the deep optical ...integral-field spectroscopy from the Multi Unit Spectroscopic Explorer (MUSE) HUDF Survey and multiwavelength data to uniquely associate all 16 line emitters, detected in the ALMA data without preselection, with rotational transitions of carbon monoxide (CO). We identify 10 as CO(2-1) at 1 < z < 2, 5 as CO(3-2) at 2 < z < 3, and 1 as CO(4-3) at z = 3.6. Using the MUSE data as a prior, we identify two additional CO(2-1) emitters, increasing the total sample size to 18. We infer metallicities consistent with (super-)solar for the CO-detected galaxies at z ≤ 1.5, motivating our choice of a Galactic conversion factor between CO luminosity and molecular gas mass for these galaxies. Using deep Chandra imaging of the HUDF, we determine an X-ray AGN fraction of 20% and 60% among the CO emitters at z ∼ 1.4 and z ∼ 2.6, respectively. Being a CO-flux-limited survey, ASPECS-LP detects molecular gas in galaxies on, above, and below the main sequence (MS) at z ∼ 1.4. For stellar masses ≥1010 (1010.5) , we detect about 40% (50%) of all galaxies in the HUDF at 1 < z < 2 (2 < z < 3). The combination of ALMA and MUSE integral-field spectroscopy thus enables an unprecedented view of MS galaxies during the peak of galaxy formation.
We investigate the CO excitation and interstellar medium (ISM) conditions in a cold gas mass-selected sample of 22 star-forming galaxies at z = 0.46-3.60, observed as part of the ALMA Spectroscopic ...Survey in the Hubble Ultra Deep Field (ASPECS). Combined with Very Large Array follow-up observations, we detect a total of 34 CO transitions with J = 1 up to 8 (and an additional 21 upper limits, up to J = 10) and 6 and transitions (and 12 upper limits). The CO(2-1) and CO(3-2)-selected galaxies, at and 2.5, respectively, exhibit a range in excitation in their mid-J = 4, 5 and high-J = 7, 8 lines, on average lower than ( -brighter) BzK-color- and submillimeter-selected galaxies at similar redshifts. The former implies that a warm ISM component is not necessarily prevalent in gas mass-selected galaxies at . We use stacking and Large Velocity Gradient models to measure and predict the average CO ladders at z < 2 and z ≥ 2, finding and , respectively. From the models, we infer that the galaxies at z ≥ 2 have intrinsically higher excitation than those at z < 2. This fits a picture in which the global excitation is driven by an increase in the star formation rate surface density of galaxies with redshift. We derive a neutral atomic carbon abundance of , comparable to the Milky Way and main-sequence galaxies at similar redshifts, and fairly high densities (≥104 cm−3), consistent with the low-J CO excitation. Our results imply a decrease in the cosmic molecular gas mass density at z ≥ 2 compared to previous ASPECS measurements.
We present a CO and atomic fine-structure line-luminosity function analysis using the ALMA Spectroscopic Survey (ASPECS) in the Hubble Ultra Deep Field. ASPECS consists of two spatially overlapping ...mosaics that cover the entire ALMA 3 mm and 1.2 mm bands. We combine the results of a line-candidate search of the 1.2 mm data cube with those previously obtained from the 3 mm cube. Our analysis shows that ∼80% of the line flux observed at 3 mm arises from CO(2-1) or CO(3-2) emitters at z = 1-3 ("cosmic noon"). At 1.2 mm, more than half of the line flux arises from intermediate-J CO transitions (Jup = 3-6); ∼12% from neutral carbon lines; and <1% from singly ionized carbon, C ii. This implies that future C ii intensity mapping surveys in the epoch of reionization will need to account for a highly significant CO foreground. The CO luminosity functions probed at 1.2 mm show a decrease in the number density at a given line luminosity (in units of L′) at increasing Jup and redshift. Comparisons between the CO luminosity functions for different CO transitions at a fixed redshift reveal subthermal conditions on average in galaxies up to z ∼ 4. In addition, the comparison of the CO luminosity functions for the same transition at different redshifts reveals that the evolution is not driven by excitation. The cosmic density of molecular gas in galaxies, H2, shows a redshift evolution with an increase from high redshift up to z ∼ 1.5 followed by a factor ∼6 drop down to the present day. This is in qualitative agreement with the evolution of the cosmic star formation rate density, suggesting that the molecular gas depletion time is approximately constant with redshift, after averaging over the star-forming galaxy population.
Abstract We present detailed spectral energy distribution (SED) modeling of 14 local ultraluminous infrared galaxies (ULIRGs) with outstanding photometric data from the literature covering the ...ultraviolet–infrared (FIR) and radio bands (∼50 MHz to ∼30 GHz). We employ the CIGALE SED fitting code to model the ultraviolet–FIR–radio SED. For the radio-only SED modeling, we use the UltraNest package, leveraging its nested sampling algorithm. Combining the results from our previous study on 11 luminous infrared galaxies (LIRGs), we discuss the global astrophysical properties of a sample of 25 starburst galaxies ( z < 0.5). Their radio spectra are frequently characterized by bends and turnovers, with no indication of ULIRGs exhibiting more complicated SEDs than LIRGs despite showing more signs of interactions. Including radio measurements in the CIGALE modeling constrained the dust luminosity and star formation rate (SFR) estimates by more than 1 order of magnitude better than previously reported for starburst galaxies. We show that total and nonthermal radio luminosity at 1.4 and 4.8 GHz frequencies can be good estimators of recent SFRs for all LIRGs and those ULIRGS with an insignificant influence of active galactic nuclei. A weaker but still significant correlation is observed between radio SFRs at 1.4 GHz and old (averaged over 100 Myr) SFRs based on SED modeling, indicative of multiple episodes of starburst activity during their lifetime. The thermal radio luminosity at 4.8 GHz is a better tracer of recent star formation than the thermal luminosity at 1.4 GHz. Statistically, our modeled nonthermal radio spectral indices do not significantly correlate with redshift, stellar mass, SFR, specific SFR, and dust mass.
The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) provides new constraints for galaxy formation models on the molecular gas properties of galaxies. We compare results from ASPECS ...to predictions from two cosmological galaxy formation models: the IllustrisTNG hydrodynamical simulations and the Santa Cruz semianalytic model (SC SAM). We explore several recipes to model the H2 content of galaxies, finding them to be consistent with one another, and take into account the sensitivity limits and survey area of ASPECS. For a canonical CO-to-H2 conversion factor of CO = 3.6 M /(K km s−1 pc2) the results of our work include: (1) the H2 mass of z > 1 galaxies predicted by the models as a function of their stellar mass is a factor of 2-3 lower than observed; (2) the models do not reproduce the number of H2-rich ( ) galaxies observed by ASPECS; (3) the H2 cosmic density evolution predicted by IllustrisTNG (the SC SAM) is in tension (in tension but with less disagreement than IllustrisTNG) with the observed cosmic density, even after accounting for the ASPECS selection function and field-to-field variance effects. The tension between models and observations at z > 1 can be alleviated by adopting a CO-to-H2 conversion factor in the range CO = 2.0-0.8 M /(K km s−1 pc2). Additional work on constraining the CO-to-H2 conversion factor and CO excitation conditions of galaxies through observations and theory will be necessary to more robustly test the success of galaxy formation models.
We present the results from the 1.2 mm continuum image obtained as part of the Atacama Large Millimeter/submillimeter Array Spectroscopic Survey in the Hubble Ultra Deep Field. The 1.2 mm continuum ...image has a size of 2.9 (4.2) arcmin2 within a primary beam response of 50% (10%) and an rms value of . We detect 35 sources at high significance (Fidelity ≥0.5); 32 have well-characterized near-infrared Hubble Space Telescope counterparts. We estimate the 1.2 mm number counts to flux levels of in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below S ∼ 0.1 mJy. In a second approach, we constrain the number counts using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jackknifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting ∼93% (∼100% if we include the lower fidelity detections) of the total continuum dust emission associated with galaxies in the Hubble Ultra Deep Field. The ancillary data allow us to study the dependence of the 1.2 mm number counts on redshift (z = 0−4), galaxy dust mass ( ), stellar mass ( ), and star formation rate ( ). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift.