ABSTRACT ALMA Cycle 2 observations of long-wavelength dust emission in 145 star-forming galaxies are used to probe the evolution of the star-forming interstellar medium (ISM). We also develop a ...physical basis and empirical calibration (with 72 low-z and z ∼ 2 galaxies) for using the dust continuum as a quantitative probe of ISM masses. The galaxies with the highest star formation rates (SFRs) at = 2.2 and 4.4 have gas masses up to 100 times that of the Milky Way and gas mass fractions reaching 50%-80%, i.e., gas masses 1-4× their stellar masses. We find a single high-z star formation law: yr−1-an approximately linear dependence on the ISM mass and an increased star formation efficiency per unit gas mass at higher redshift. Galaxies above the main sequence (MS) have larger gas masses but are converting their ISM into stars on a timescale only slightly shorter than those on the MS; thus, these "starbursts" are largely the result of having greatly increased gas masses rather than an increased efficiency of converting gas to stars. At z > 1, the entire population of star-forming galaxies has ∼2-5 times shorter gas depletion times than low-z galaxies. These shorter depletion times indicate a different mode of star formation in the early universe-most likely dynamically driven by compressive, high-dispersion gas motions-a natural consequence of the high gas accretion rates.
ALMA observations of the long wavelength dust continuum are used to estimate the interstellar medium (ISM) masses in a sample of 708 galaxies at z = 0.3 to 4.5 in the COSMOS field. The galaxy sample ...has known far-infrared luminosities and, hence, star formation rates (SFRs) and stellar masses ( ) from the optical-infrared spectrum fitting. The galaxies sample SFRs from the main sequence (MS) to 50 times above the MS. The derived ISM masses are used to determine the dependence of gas mass on redshift, , and specific SFR (sSFR) relative to the MS. The ISM masses increase approximately with the 0.63 power of the rate of increase in SFRs with redshift and the 0.32 power of the sSFR/sSFRMS. The SF efficiencies also increase as the 0.36 power of the SFR redshift evolution and the 0.7 power of the elevation above the MS; thus the increased activities at early epochs are driven by both increased ISM masses and SF efficiency. Using the derived ISM mass function, we estimate the accretion rates of gas required to maintain continuity of the MS evolution ( yr−1 at z > 2.5). Simple power-law dependencies are similarly derived for the gas accretion rates. We argue that the overall evolution of galaxies is driven by the rates of gas accretion. The cosmic evolution of total ISM mass is estimated and linked to the evolution of SF and active galactic nucleus activity at early epochs.
Abstract
We present a catalog of 10,718 objects in the COSMOS field, observed through multi-slit spectroscopy with the Deep Imaging Multi-Object Spectrograph (DEIMOS) on the Keck II telescope in the ...wavelength range ∼5500–9800 Å. The catalog contains 6617 objects with high-quality spectra (two or more spectral features), and 1798 objects with a single spectroscopic feature confirmed by the photometric redshift. For 2024 typically faint objects, we could not obtain reliable redshifts. The objects have been selected from a variety of input catalogs based on multi-wavelength observations in the field, and thus have a diverse selection function, which enables the study of the diversity in the galaxy population. The magnitude distribution of our objects is peaked at
I
AB
∼ 23 and
K
AB
∼ 21, with a secondary peak at
K
AB
∼ 24. We sample a broad redshift distribution in the range 0 <
z
< 6, with one peak at
z
∼ 1, and another one around
z
∼ 4. We have identified 13 redshift spikes at
z
> 0.65 with chance probabilities < 4 × 10
−4
, some of which are clearly related to protocluster structures of sizes >10 Mpc. An object-to-object comparison with a multitude of other spectroscopic samples in the same field shows that our DEIMOS sample is among the best in terms of fraction of spectroscopic failures and relative redshift accuracy. We have determined the fraction of spectroscopic blends to about 0.8% in our sample. This is likely a lower limit and at any rate well below the most pessimistic expectations. Interestingly, we find evidence for strong lensing of Ly
α
background emitters within the slits of 12 of our target galaxies, increasing their apparent density by about a factor of 4.
Using a compilation of 25 studies from the literature, we investigate the evolution of the star-forming galaxy (SFG) main sequence (MS) in stellar mass and star formation rate (SFR) out to z ~ 6. ...After converting all observations to a common set of calibrations, we find a remarkable consensus among MS observations (~0.1 dex 1sigma interpublication scatter). By fitting for time evolution of the MS in bins of constant mass, we deconvolve the observed scatter about the MS within each observed redshift bin. After accounting for observed scatter between different SFR indicators, we find the width of the MS distribution is ~0.2 dex and remains constant over cosmic time. Our best fits indicate the slope of the MS is likely time-dependent, with our best-fit log SFR(Mlow *,t) = (0.84 + or - 0.02 - 0.026 + or - 0.003 x t) log Mlow *-(6.51 + or - 0.24 - 0.11 + or - 0.03 x t), where t is the age of the universe in Gyr. We use our fits to create empirical evolutionary tracks in order to constrain MS galaxy star formation histories (SFHs), finding that (1) the most accurate representations of MS SFHs are given by delayed-Tau models, (2) the decline in fractional stellar mass growth for a "typical" MS galaxy today is approximately linear for most of its lifetime, and (3) scatter about the MS can be generated by galaxies evolving along identical evolutionary tracks assuming an initial 1Sigma spread in formation times of ~1.4 Gyr.
ABSTRACT We study the relationship between stellar mass, star formation rate (SFR), ionization state, and gas-phase metallicity for a sample of 41 normal star-forming galaxies at 3 z 3.7. The ...gas-phase oxygen abundance, ionization parameter, and electron density of ionized gas are derived from rest-frame optical strong emission lines measured on near-infrared spectra obtained with Keck/Multi-Object Spectrograph for Infra-Red Exploration. We remove the effect of these strong emission lines in the broadband fluxes to compute stellar masses via spectral energy distribution fitting, while the SFR is derived from the dust-corrected ultraviolet luminosity. The ionization parameter is weakly correlated with the specific SFR, but otherwise the ionization parameter and electron density do not correlate with other global galaxy properties such as stellar mass, SFR, and metallicity. The mass-metallicity relation (MZR) at z 3.3 shows lower metallicity by 0.7 dex than that at z = 0 at the same stellar mass. Our sample shows an offset by 0.3 dex from the locally defined mass-metallicity-SFR relation, indicating that simply extrapolating such a relation to higher redshift may predict an incorrect evolution of MZR. Furthermore, within the uncertainties we find no SFR-metallicity correlation, suggesting a less important role of SFR in controlling the metallicity at high redshift. We finally investigate the redshift evolution of the MZR by using the model by Lilly et al., finding that the observed evolution from z = 0 to z 3.3 can be accounted for by the model assuming a weak redshift evolution of the star formation efficiency.
We examine the behaviour of the infrared-radio correlation (IRRC) over the range 0 <z ≲ 6 using new, highly sensitive 3 GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared ...data from the Herschel Space Observatory in the 2 deg2 COSMOS field. We distinguish between objects where emission is believed to arise solely from star-formation, and those where an active galactic nucleus (AGN) is thought to be present. We account for non-detections in the radio or in the infrared using a doubly-censored survival analysis. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4 GHz radio luminosity ratio (qTIR), decreases with increasing redshift: qTIR(z) = (2.88 ± 0.03)(1 + z)− 0.19 ± 0.01. This is consistent with several previous results from the literature. Moderate-to-high radiative luminosity AGN do not follow the same qTIR(z) trend as star-forming galaxies, having a lower normalisation and steeper decrease with redshift. We cannot rule out the possibility that unidentified AGN contributions only to the radio regime may be steepening the observed qTIR(z) trend of the star-forming galaxy population. We demonstrate that the choice of the average radio spectral index directly affects the normalisation, as well as the derived trend with redshift of the IRRC. An increasing fractional contribution to the observed 3 GHz flux by free-free emission of star-forming galaxies may also affect the derived evolution. However, we find that the standard (M82-based) assumption of the typical radio spectral energy distribution (SED) for star-forming galaxies is inconsistent with our results. This suggests a more complex shape of the typical radio SED for star-forming galaxies, and that imperfect K corrections in the radio may govern the derived trend of decreasing qTIR with increasing redshift. A more detailed understanding of the radio spectrum is therefore required for robust K corrections in the radio and to fully understand the cosmic evolution of the IRRC. Lastly, we present a redshift-dependent relation between rest-frame 1.4 GHz radio luminosity and star formation rate taking the derived redshift trend into account.
We study the composition of the faint radio population selected from the Karl G. Jansky Very Large Array Cosmic Evolution Survey (VLA-COSMOS) 3 GHz Large Project, which is a radio continuum survey ...performed at 10 cm wavelength. The survey covers a 2.6 square degree area with a mean rms of ~ 2.3 μJy/beam, cataloging 10 830 sources above 5σ, and enclosing the full 2 square degree COSMOS field. By combining these radio data with optical, near-infrared (UltraVISTA), and mid-infrared (Spitzer/IRAC) data, as well as X-ray data (Chandra), we find counterparts to radio sources for ~93% of the total radio sample reaching out to z ≲ 6; these sources are found in the unmasked areas of the COSMOS field, i.e., those not affected by saturated or bright sources in the optical to near-infrared (NIR) bands. We further classify the sources as star-forming galaxies or AGN based on various criteria, such as X-ray luminosity; observed mid-infrared color; UV–far-infrared spectral energy distribution; rest-frame, near-UV optical color that is corrected for dust extinction; and radio excess relative to that expected from the star formation rate of the hosts. We separate the AGN into subsamples dominated by low-to-moderate and moderate-to-high radiative luminosity AGN, i.e., candidates for high-redshift analogs to local low- and high-excitation emission line AGN, respectively. We study the fractional contributions of these subpopulations down to radio flux levels of ~11 μJy at 3 GHz (or ~20 μJy at 1.4 GHz assuming a spectral index of –0.7). We find that the dominant fraction at 1.4 GHz flux densities above ~200 μJy is constituted of low-to-moderate radiative luminosity AGN. Below densities of ~100 μJy the fraction of star-forming galaxies increases to ~ 60%, followed by the moderate-to-high radiative luminosity AGN (~ 20%) and low-to-moderate radiative luminosity AGN (~ 20%). Based on this observational evidence, we extrapolate the fractions down to sensitivities of the Square Kilometer Array (SKA). Our estimates suggest that at the faint flux limits to be reached by the (Wide, Deep, and UltraDeep) SKA1 surveys, a selection based only on radio flux limits can provide a simple tool to efficiently identify samples highly (>75%) dominated by star-forming galaxies.
We study the efficiency of galactic feedback in the early Universe by stacking the C II 158
μ
m emission in a large sample of normal star-forming galaxies at 4 <
z
< 6 from the ALMA Large Program ...to INvestigate C II at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked C II profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked C II spectrum, with velocities of |
v
|≲500 km s
−1
. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFR
med
= 25
M
⊙
yr
−1
), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at
z
> 4. From the stacking analysis of the datacubes, we find that the combined C II core emission (|
v
|< 200 km s
−1
) of the higher-SFR galaxies is extended on physical sizes of ∼30 kpc (diameter scale), well beyond the analogous C II core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.
We use the large COSMOS sample of galaxies to study in an internally self-consistent way the change in the number densities of quenched early-type galaxies (Q-ETGs) of a given size over the redshift ...interval 0.2 < z < 1 in order to study the claimed size evolution of these galaxies. In a stellar mass bin at 10 super(10.5) < M sub(galaxy) < 10 super(11) M sub(middot in circle), we see no change in the number density of compact Q-ETGs over this redshift range, while in a higher mass bin at >10 super(11) M sub(middot in circle), where we would expect merging to be more significant, we find a small decrease, by ~30%. In both mass bins, the increase of the median sizes of Q-ETGs with time is primarily caused by the addition to the size function of larger and more diffuse Q-ETGs. At all masses, compact Q-ETGs become systematically redder toward later epochs, with a (U - V) color difference which is consistent with a passive evolution of their stellar populations, indicating that they are a stable population that does not appreciably evolve in size. We find furthermore, at all epochs, that the larger Q-ETGs (at least in the lower mass bin) have average rest-frame colors that are systematically bluer than those of the more compact Q-ETGs, suggesting that the former are indeed younger than the latter. The idea that new, large, Q-ETGs are responsible for the observed growth in the median size of the population at a given mass is also supported by analysis of the sizes and number of the star-forming galaxies that are expected to be the progenitors of the new Q-ETGs over the same period. In the low mass bin, the new Q-ETGs appear to have similar to 30% smaller half-light radii than their star-forming progenitors. This is likely due to the fading of their disks after they cease star formation. Comparison with higher redshifts shows that the median size of newly quenched galaxies roughly scales, at constant mass, as (1 + z) super(-1). We conclude that the dominant cause of the size evolution seen in the Q-ETG population is that the average sizes and thus stellar densities of individual Q-ETGs roughly scale with the average density of the universe at the time when they were quenched, and that subsequent size changes in individual objects, through merging or other processes, are of secondary importance, especially at masses below 10 super(11) M sub(middot in circle).
The rest-frame ultraviolet properties of galaxies during the first three billion years of cosmic time (redshift z > 4) indicate a rapid evolution in the dust obscuration of such galaxies. This ...evolution implies a change in the average properties of the interstellar medium, but the measurements are systematically uncertain owing to untested assumptions and the inability to detect heavily obscured regions of the galaxies. Previous attempts to measure the interstellar medium directly in normal galaxies at these redshifts have failed for a number of reasons, with two notable exceptions. Here we report measurements of the forbidden C ii emission (that is, C II) from gas, and the far-infrared emission from dust, in nine typical star-forming galaxies about one billion years after the Big Bang (z ≈ 5-6). We find that these galaxies have thermal emission that is less than 1/12 that of similar systems about two billion years later, and enhanced C II emission relative to the far-infrared continuum, confirming a strong evolution in the properties of the interstellar medium in the early Universe. The gas is distributed over scales of one to eight kiloparsecs, and shows diverse dynamics within the sample. These results are consistent with early galaxies having significantly less dust than typical galaxies seen at z < 3 and being comparable in dust content to local low-metallicity systems.