No Arabic abstract
We describe the CO Luminosity Density at High-z (COLDz) survey, the first spectral line deep field targeting CO(1-0) emission from galaxies at $z=1.95-2.85$ and CO(2-1) at $z=4.91-6.70$. The main goal of COLDz is to constrain the cosmic density of molecular gas at the peak epoch of cosmic star formation. By targeting both a wide ($sim$51 arcmin$^2$) and a deep area ($sim$9 arcmin$^2$), the survey is designed to robustly constrain the bright end and the characteristic luminosity of the CO(1-0) luminosity function. An extensive analysis of the reliability of our line candidates, and new techniques provide detailed completeness and statistical corrections as necessary to determine the best constraints to date on the CO luminosity function. Our blind search for CO(1-0) uniformly selects starbursts and massive Main Sequence galaxies based on their cold molecular gas masses. Our search also detects CO(2-1) line emission from optically dark, dusty star-forming galaxies at $z>5$. We find a range of spatial sizes for the CO-traced gas reservoirs up to $sim40$ kpc, suggesting that spatially extended cold molecular gas reservoirs may be common in massive, gas-rich galaxies at $zsim2$. Through CO line stacking, we constrain the gas mass fraction in previously known typical star-forming galaxies at $z=2$-3. The stacked CO detection suggests lower molecular gas mass fractions than expected for massive Main Sequence galaxies by a factor of $sim3-6$. We find total CO line brightness at $sim34,$GHz of $0.45pm0.2,mu$K, which constrains future line intensity mapping and CMB experiments.
The broad spectral bandwidth at mm and cm-wavelengths provided by the recent upgrades to the Karl G. Jansky Very Large Array (VLA) has made it possible to conduct unbiased searches for molecular CO line emission at redshifts, z > 1.31. We present the discovery of a gas-rich, star-forming galaxy at z = 2.48, through the detection of CO(1-0) line emission in the COLDz survey, through a sensitive, Ka-band (31 to 39 GHz) VLA survey of a 6.5 square arcminute region of the COSMOS field. We argue that the broad line (FWHM ~570 +/- 80 km/s) is most likely to be CO(1-0) at z=2.48, as the integrated emission is spatially coincident with an infrared-detected galaxy with a photometric redshift estimate of z = 3.2 +/- 0.4. The CO(1-0) line luminosity is L_CO = (2.2 +/- 0.3) x 10^{10} K km/s pc^2, suggesting a cold molecular gas mass of M_gas ~ (2 - 8)x10^{10}M_solar depending on the assumed value of the molecular gas mass to CO luminosity ratio alpha_CO. The estimated infrared luminosity from the (rest-frame) far-infrared spectral energy distribution (SED) is L_IR = 2.5x10^{12} L_solar and the star-formation rate is ~250 M_solar/yr, with the SED shape indicating substantial dust obscuration of the stellar light. The infrared to CO line luminosity ratio is ~114+/-19 L_solar/(K km/s pc^2), similar to galaxies with similar SFRs selected at UV/optical to radio wavelengths. This discovery confirms the potential for molecular emission line surveys as a route to study populations of gas-rich galaxies in the future.
We report the first detailed measurement of the shape of the CO luminosity function at high redshift, based on $>$320 hr of the NSFs Karl G. Jansky Very Large Array (VLA) observations over an area of $sim$60 arcmin$^2$ taken as part of the CO Luminosity Density at High Redshift (COLDz) survey. COLDz blindly selects galaxies based on their cold gas content through CO($J$=1$to$0) emission at $z$$sim$2-3 and CO($J$=2$to$1) at $z$$sim$5-7 down to a CO luminosity limit of log($L_{rm CO}$/K km s$^{-1}$ pc$^2$)$simeq$9.5. We find that the characteristic luminosity and bright end of the CO luminosity function are substantially higher than predicted by semi-analytical models, but consistent with empirical estimates based on the infrared luminosity function at $z$$sim$2. We also present the currently most reliable measurement of the cosmic density of cold gas in galaxies at early epochs, i.e., the cold gas history of the universe, as determined over a large cosmic volume of $sim$375,000 Mpc$^3$. Our measurements are in agreement with an increase of the cold gas density from $z$$sim$0 to $z$$sim$2-3, followed by a possible decline towards $z$$sim$5-7. These findings are consistent with recent surveys based on higher-$J$ CO line measurements, upon which COLDz improves in terms of statistical uncertainties by probing $sim$50-100 times larger areas and in the reliability of total gas mass estimates by probing the low-$J$ CO lines accessible to the VLA. Our results thus appear to suggest that the cosmic star-formation rate density follows an increased cold molecular gas content in galaxies towards its peak about 10 billion years ago, and that its decline towards the earliest epochs is likely related to a lower overall amount of cold molecular gas (as traced by CO) bound in galaxies towards the first billion years after the Big Bang.
We present Karl G. Jansky Very Large Array (VLA) observations of 44 GHz continuum and CO J=2-1 line emission in BR1202-0725 at z=4.7 (a starburst galaxy and quasar pair) and BRI1335-0417 at z=4.4 (also hosting a quasar). With the full 8 GHz bandwidth capabilities of the upgraded VLA, we study the (rest-frame) 250 GHz thermal dust continuum emission for the first time along with the cold molecular gas traced by the Low-J CO line emission. The measured CO J=2-1 line luminosities of BR1202-0725 are L(CO) = (8.7+/-0.8)x10^10 K km/s pc^2 and L(CO) = (6.0+/-0.5)x10^10 K km/s pc^2 for the submm galaxy (SMG) and quasar, which are equal to previous measurements of the CO J=5-4 line luminosities implying thermalized line emission and we estimate a combined cold molecular gas mass of ~9x10^10 Msun. In BRI1335-0417 we measure L(CO) = (7.3+/-0.6)x10^10 K km/s pc^2. We detect continuum emission in the SMG BR1202-0725 North (S(44GHz) = 51+/-6 microJy), while the quasar is detected with S(44GHz) = 24+/-6 microJy and in BRI1335-0417 we measure S(44GHz) = 40+/-7 microJy. Combining our continuum observations with previous data at (rest-frame) far-infrared and cm-wavelengths, we fit three component models in order to estimate the star-formation rates. This spectral energy distribution fitting suggests that the dominant contribution to the observed 44~GHz continuum is thermal dust emission, while either thermal free-free or synchrotron emission contributes less than 30%.
The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($approx$2.5), sensitivity (a 1$sigma$ goal of 70 $mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an on the fly interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.
We present a population of 20 radio-luminous supernovae (SNe) with emission reaching $L_{ u}{sim}10^{26}-10^{29}rm{erg s^{-1} Hz^{-1}}$ in the first epoch of the Very Large Array Sky Survey (VLASS) at $2-4$ GHz. Our sample includes one long Gamma-Ray Burst, SN 2017iuk/GRB171205A, and 19 core-collapse SNe detected at $approx (1-60)$ years after explosion. No thermonuclear explosion shows evidence for bright radio emission, and hydrogen-poor progenitors dominate the sub-sample of core-collapse events with spectroscopic classification at the time of explosion (73%). We interpret these findings into the context of the expected radio emission from the forward shock interaction with the circumstellar medium (CSM). We conclude that these observations require a departure from the single wind-like density profile (i.e., $rho_{rm{CSM}}propto r^{-2}$) that is expected around massive stars and/or a departure from a spherical Newtonian shock. Viable alternatives include the shock interaction with a detached, dense shell of CSM formed by a large effective progenitor mass-loss rate $dot M sim (10^{-4}-10^{-1})$ M$_{odot}$ yr$^{-1}$ (for an assumed wind velocity of $1000,rm{km,s^{-1}}$); emission from an off-axis relativistic jet entering our line of sight; or the emergence of emission from a newly-born pulsar-wind nebula. The relativistic SN,2012ap that is detected 5.7 and 8.5 years after explosion with $L_{ u}{sim}10^{28}$ erg s$^{-1}$ Hz$^{-1}$ might constitute the first detections of an off-axis jet+cocoon system in a massive star. Future multi-wavelength observations will distinguish among these scenarios. Our VLASS source catalogs, which were used to perform the VLASS cross matching, are publicly available at https://doi.org/10.5281/zenodo.4895112.