Enabling the next generation of cm-wavelength studies of high-redshift molecular gas with the SKA

The Square Kilometre Array will be a revolutionary instrument for the study of gas in the distant Universe. SKA1 will have sufficient sensitivity to detect and image atomic 21 cm HI in individual galaxies at significant cosmological distances, complementing ongoing ALMA imaging of redshifted high-J CO line emission and far-infrared interstellar medium lines such as [CII] 157.7 um. At frequencies below ~50 GHz, observations of redshifted emission from low-J transitions of CO, HCN, HCO+, HNC, H2O and CS provide insight into the kinematics and mass budget of the cold, dense star-forming gas in galaxies. In advance of ALMA band 1 deployment (35 to 52 GHz), the most sensitive facility for high-redshift studies of molecular gas operating below 50~GHz is the Karl G. Jansky Very Large Array (VLA). Here, we present an overview of the role that the SKA could play in molecular emission line studies during SKA1 and SKA2, with an emphasis on studies of the dense gas tracers directly probing regions of active star-formation.

COLDz: Karl G. Jansky Very Large Array discovery of a gas-rich galaxy in COSMOS

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.

Karl G. Jansky Very Large Array observations of cold dust and molecular gas in starbursting quasar host galaxies at z~4.5

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%.

[CII] line emission in massive star-forming galaxies at z=4.7

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the [CII] 157.7micron fine structure line and thermal dust continuum emission from a pair of gas-rich galaxies at z=4.7, BR1202-0725. This system consists of a luminous quasar host galaxy and a bright submm galaxy (SMG), while a fainter star-forming galaxy is also spatially coincident within a 4 (25 kpc) region. All three galaxies are detected in the submm continuum, indicating FIR luminosities in excess of 10^13 Lsun for the two most luminous objects. The SMG and the quasar host galaxy are both detected in [CII] line emission with luminosities, L([CII]) = (10.0 +/- 1.5)x10^9 Lsun and L([CII]) = (6.5+/-1.0)x10^9 Lsun, respectively. We estimate a luminosity ratio, L([CII])/L(FIR) = (8.3+/-1.2)x10^-4 for the starburst SMG to the North, and L([CII])/L(FIR) = (2.5+/-0.4)x10^-4 for the quasar host galaxy, in agreement with previous high-redshift studies that suggest lower [CII]-to-FIR luminosity ratios in quasars than in starburst galaxies. The third fainter object with a flux density, S(340GHz) = 1.9+/-0.3 mJy, is coincident with a Ly-Alpha emitter and is detected in HST ACS F775W and F814W images but has no clear counterpart in the H-band. Even if this third companion does not lie at a similar redshift to BR1202-0725, the quasar and the SMG represent an overdensity of massive, infrared luminous star-forming galaxies within 1.3 Gyr of the Big Bang.

An EVLA search for water megamaser emission in the submm galaxy SMM J16359+6612 at z=2.5

Using the Expanded Very Large Array, we have conducted a search for 22.2 GHz H2O megamaser emission in the strongly lensed submm galaxy, SMM J16359+6612 at z=2.517. This object is lensed into three components, and after a correction for magnification is applied to its submm-wavelength flux density, it is typical of the bulk of the high-redshift, submm galaxy population responsible for the 850 um extragalactic background (S(850um)~1 mJy). We do not detect any H2O megamaser emission, but the lensing allows us to place an interesting constraint on the luminosity of any megamasers present, L(H2O) < 5305 solar luminosities for an assumed linewidth of 80 km/s. Because the far-infrared luminosity in submm galaxies is mainly powered by star formation, and very luminous H2O megamasers are more commonly associated with quasar activity, it could be that blind searches for H2O megamasers will not be an effective means of determining redshifts for less luminous members of the submm galaxy population.

Thermal Emission from Warm Dust in the Most Distant Quasars

We report new continuum observations of fourteen z~6 quasars at 250 GHz and fourteen quasars at 1.4 GHz. We summarize all recent millimeter and radio observations of the sample of the thirty-three quasars known with 5.71<=z<=6.43, and present a study of the rest frame far-infrared (FIR) properties of this sample. These quasars were observed with the Max Plank Millimeter Bolometer Array (MAMBO) at 250 GHz with mJy sensitivity, and 30% of them were detected. We also recover the average 250 GHz flux density of the MAMBO undetected sources at 4 sigma, by stacking the on-source measurements. The derived mean radio-to-UV spectral energy distributions (SEDs) of the full sample and the 250 GHz non-detections show no significant difference from that of lower-redshift optical quasars. Obvious FIR excesses are seen in the individual SEDs of the strong 250 GHz detections, with FIR-to-radio emission ratios consistent with that of typical star forming galaxies. Most 250 GHz-detected sources follow the L_{FIR}--L_{bol} relationship derived from a sample of local IR luminous quasars (L_{IR}>10^{12}L_{odot}), while the average L_{FIR}/L_{bol} ratio of the non-detections is consistent with that of the optically-selected PG quasars. The MAMBO detections also tend to have weaker Lyalpha emission than the non-detected sources. We discuss possible FIR dust heating sources, and critically assess the possibility of active star formation in the host galaxies of the z~6 quasars. The average star formation rate of the MAMBO non-detections is likely to be less than a few hundred M_{odot} yr^{-1}, but in the strong detections, the host galaxy star formation is probably at a rate of gtrsim10^{3} M_{odot} yr^{-1}, which dominates the FIR dust heating.