We report the detection of CO($J$=2$\to$1) emission from three massive dusty starburst galaxies at $z$$>$5 through molecular line scans in the NSF's Karl G. Jansky Very Large Array (VLA) CO Luminosity Density at High Redshift (COLDz) survey. Redshifts for two of the sources, HDF 850.1 ($z$=5.183) and AzTEC-3 ($z$=5.298), were previously known. We revise a previous redshift estimate for the third source GN10 ($z$=5.303), which we have independently confirmed through detections of CO $J$=1$\to$0, 5$\to$4, 6$\to$5, and [CII] 158 $\mu$m emission with the VLA and the NOrthern Extended Milllimeter Array (NOEMA). We find that two currently independently confirmed CO sources in COLDz are "optically dark", and that three of them are dust-obscured galaxies at $z$$>$5. Given our survey area of $\sim$60 arcmin$^2$, our results appear to imply a $\sim$6-55 times higher space density of such distant dusty systems within the first billion years after the Big Bang than previously thought. At least two of these $z$$>$5 galaxies show star-formation rate surface densities consistent with so-called "maximum" starbursts, but we find significant differences in CO excitation between them. This result may suggest that different fractions of the massive gas reservoirs are located in the dense, star-forming nuclear regions - consistent with the more extended sizes of the [CII] emission compared to the dust continuum and higher [CII]-to-far-infrared luminosity ratios in those galaxies with lower gas excitation. We thus find substantial variations in the conditions for star formation between $z$$>$5 dusty starbursts, which typically have dust temperatures $\sim$57%$\pm$25% warmer than starbursts at $z$=2-3 due to their enhanced star formation activity.

We present ~1" resolution (~2 kpc in the source plane) observations of the CO(1-0), CO(3-2), Halpha, and [N II] lines in the strongly-lensed z=2.26 star-forming galaxy SDSS J0901+1814. We use these observations to constrain the lensing potential of a foreground group of galaxies, and our source-plane reconstructions indicate that SDSS J0901+1814 is a nearly face-on (i~30 degrees) massive disk with r_{1/2}>~4 kpc for its molecular gas. Using our new magnification factors (mu_tot~30), we find that SDSS J0901+1814 has a star formation rate (SFR) of 268^{+63}_{-61} M_sun/yr, M_gas=(1.6^{+0.3}_{-0.2})x10^11x(alpha_CO/4.6) M_sun, and M_star=(9.5^{+3.8}_{-2.8})x10^10 M_sun, which places it on the star-forming galaxy "main sequence." We use our matched high-angular resolution gas and SFR tracers (CO and Halpha, respectively) to perform a spatially resolved (pixel-by-pixel) analysis of SDSS J0901+1814 in terms of the Schmidt-Kennicutt relation. After correcting for the large fraction of obscured star formation (SFR_Halpha/SFR_TIR=0.054^{+0.015}_{-0.014}), we find SDSS J0901+1814 is offset from "normal" star-forming galaxies to higher star formation efficiencies independent of assumptions for the CO-to-H_2 conversion factor. Our mean best-fit index for the Schmidt-Kennicutt relation for SDSS J0901+1814, evaluated with different CO lines and smoothing levels, is n=1.54+/-0.13; however, the index may be affected by gravitational lensing, and we find n=1.24+/-0.02 when analyzing the source-plane reconstructions. While the Schmidt-Kennicutt index largely appears unaffected by which of the two CO transitions we use to trace the molecular gas, the source-plane reconstructions and dynamical modeling suggest that the CO(1-0) emission is more spatially extended than the CO(3-2) emission.

In order to better understand how active galactic nuclei (AGN) effect the interstellar media of their host galaxies, we perform a meta-analysis of the CO emission for a sample of $z=0.01-4$ galaxies from the literature with existing CO detections and well-constrained AGN contributions to the infrared (67 galaxies). Using either Spitzer/IRS mid-IR spectroscopy or Spitzer+Herschel colors we determine the fraction of the infrared luminosity in each galaxy that can be attributed to heating by the AGN or stars. We calculate new average CO spectral line ratios (primarily from Carilli \& Walter 2013) to uniformly scale the higher-$J$ CO detections to the ground state and accurately determine our sample's molecular gas masses. We do not find significant differences in the gas depletion timescales/star formation efficiencies (SFEs) as a function of the mid-infrared AGN strength ($f_{\rm AGN}({\rm MIR})$ or $L_{\rm IR} ({\rm AGN})$), which indicates that the presence of an IR-bright AGN is not a sufficient sign-post of galaxy quenching. We also find that the dust-to-gas ratio is consistent for all sources, regardless of AGN emission, redshift, or $L_{\rm IR}$, indicating that dust is likely a reliable tracer of gas mass for massive dusty galaxies (albeit with a large degree of scatter). Lastly, if we classify galaxies as either AGN or star formation dominated, we do not find a robust statistically significant difference between their CO excitation.

We report the first detailed measurement of the shape of the CO luminosity function at high redshift, based on $>$320 hr of the NSF's 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 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 $z\sim2$. 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.45\pm0.2\,\mu$K, which constrains future line intensity mapping and CMB experiments.

We report the serendipitous discovery of a dusty, starbursting galaxy at $z=5.667$ (hereafter called CRLE) in close physical association with the "normal" main-sequence galaxy HZ10 at $z=5.654$. CRLE was identified by detection of [CII], [NII] and CO(2-1) line emission, making it the highest redshift, most luminous starburst in the COSMOS field. This massive, dusty galaxy appears to be forming stars at a rate of at least 1500$\,M_\odot$ yr$^{-1}$ in a compact region only $\sim3$ kpc in diameter. The dynamical and dust emission properties of CRLE suggest an ongoing merger driving the starburst, in a potentially intermediate stage relative to other known dusty galaxies at the same epoch. The ratio of [CII] to [NII] may suggest that an important ($\sim15\%$) contribution to the [CII] emission comes from a diffuse ionized gas component, which could be more extended than the dense, starbursting gas. CRLE appears to be located in a significant galaxy overdensity at the same redshift, potentially associated with a large-scale cosmic structure recently identified in a Lyman Alpha Emitter survey. This overdensity suggests that CRLE and HZ10 reside in a protocluster environment, offering the tantalizing opportunity to study the effect of a massive starburst on protocluster star formation. Our findings support the interpretation that a significant fraction of the earliest galaxy formation may occur from the inside out, within the central regions of the most massive halos, while rapidly evolving into the massive galaxy clusters observed in the local Universe.

We present the detection of four far-infrared fine-structure oxygen lines, as well as strong upper limits for the CO(2-1) and [N II] 205 um lines, in 3C 368, a well-studied radio-loud galaxy at z = 1.131. These new oxygen lines, taken in conjunction with previously observed neon and carbon fine-structure lines, suggest a powerful active galactic nucleus (AGN), accompanied by vigorous and extended star formation. A starburst dominated by O8 stars, with an age of ~6.5 Myr, provides a good fit to the fine-structure line data. This estimated age of the starburst makes it nearly concurrent with the latest episode of AGN activity, suggesting a link between the growth of the supermassive black hole and stellar population in this source. We do not detect the CO(2-1) line, down to a level twelve times lower than the expected value for star forming galaxies. This lack of CO line emission is consistent with recent star formation activity if the star-forming molecular gas has low metallicity, is highly fractionated (such that CO is photodissociated through much of the clouds), or is chemically very young (such that CO has not yet had time to form). It is also possible, though we argue unlikely, that the ensemble of fine structure lines are emitted from the region heated by the AGN.

We report interferometric measurements of [NII] 205 um fine-structure line emission from a representative sample of three galaxies at z=5-6 using the Atacama Large (sub)Millimeter Array (ALMA). These galaxies were previously detected in [CII] and far-infrared continuum emission and span almost two orders of magnitude in star formation rate (SFR). Our results show at least two different regimes of ionized inter-stellar medium properties for galaxies in the first billion years of cosmic time, separated by their L_[CII]/L_[NII] ratio. We find extremely low [NII] emission compared to [CII] (L_ [CII]/L_[NII]=68 [+200/-28]) from a "typical" L*_UV star-forming galaxy, likely directly or indirectly (by its effect on the radiation field) related to low dust abundance and low metallicity. The infrared-luminous modestly star-forming Lyman Break Galaxy (LBG) in our sample is characterized by an ionized-gas fraction (L_[CII]/L_[NII]<=20) typical of local star-forming galaxies and shows evidence for spatial variations in its ionized-gas fraction across an extended gas reservoir. The extreme SFR, warm and compact dusty starburst AzTEC-3 shows an ionized fraction higher than expected given its star-formation rate surface density (L_[CII]/L_[NII]=22+/-8) suggesting that [NII] dominantly traces a diffuse ionized medium rather than star-forming HII regions in this type of galaxy. This highest redshift sample of [NII] detections provides some of the first constraints on ionized and neutral gas modeling attempts and on the structure of the inter-stellar medium at z=5-6 in "normal" galaxies and starbursts.

We present CO(1-0) observations obtained at the Karl G. Jansky Very Large Array (VLA) for 14 z~2 galaxies with existing CO(3-2) measurements, including 11 galaxies which contain active galactic nuclei (AGN) and three submillimeter galaxies (SMGs). We combine this sample with an additional 15 z~2 galaxies from the literature that have both CO(1-0) and CO(3-2) measurements in order to evaluate differences in CO excitation between SMGs and AGN host galaxies, measure the effects of CO excitation on the derived molecular gas properties of these populations, and to look for correlations between the molecular gas excitation and other physical parameters. With our expanded sample of CO(3-2)/CO(1-0) line ratio measurements, we do not find a statistically significant difference in the mean line ratio between SMGs and AGN host galaxies as found in the literature, instead finding r_3,1=1.03+/-0.50 for AGN host galaxies and r_3,1=0.78+/-0.27 for SMGs (or r_3,1=0.90+/-0.40 for both populations combined). We also do not measure a statistically significant difference between the distributions of the line ratios for these populations at the p=0.05 level, although this result is less robust. We find no excitation dependence on the index or offset of the integrated Schmidt-Kennicutt relation for the two CO lines, and obtain indices consistent with N=1 for the various sub-populations. However, including low-z "normal" galaxies increases our best-fit Schmidt-Kennicutt index to N~1.2. While we do not reproduce correlations between the CO line width and luminosity, we do reproduce correlations between CO excitation and star formation efficiency.

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 present multiline CO observations of the complex submillimeter galaxy SMM J00266+1708. Using the Zpectrometer on the Green Bank Telescope, we provide the first precise spectroscopic measurement of its redshift (z=2.742). Based on followup CO(1-0), CO(3-2), and CO(5-4) mapping, SMM J00266+1708 appears to have two distinct components separated by ~500 km/s that are nearly coincident along our line of sight. The two components show hints of different kinematics, with the blue-shifted component dispersion-dominated and the red-shifted component showing a clear velocity gradient. CO line ratios differ slightly between the two components, indicating that the physical conditions in their molecular gas may not be alike. We tentatively infer that SMM J00266+1708 is an ongoing merger with a mass ratio of (7.8+/-4.0)/sin^2(i), with its overall size and surface brightness closely resembling that of other merging systems. We perform large velocity gradient modeling of the CO emission from both components and find that each component's properties are consistent with a single phase of molecular gas (i.e., a single temperatures and density); additional multi-phase modeling of the red-shifted component, although motivated by a CO(1-0) size larger than the CO(3-2) size, is inconclusive. SMM J00266+1708 provides evidence of early stage mergers within the submillimeter galaxy population. Continuum observations of J00266 at the ~1" resolution of our observations could not have distinguished between the two components due to their separation (0.73" +/- 0.06"), illustrating that the additional velocity information provided by spectral line studies is important for addressing the prevalence of unresolved galaxy pairs in low-resolution submillimeter surveys.

We present high-resolution CO(1-0) observations of the lensed submillimeter galaxy (SMG) SMM J14011+0252 at z=2.6. Comparison to the previously-detected CO(3-2) line gives an intensity ratio of r_3,1=0.97+/-0.16 in temperature units, larger than is typical for SMGs but within the range seen in the low-z ultraluminous infrared galaxy population. Combining our new data with previous mid-J CO observations, we perform a single-phase large velocity gradient (LVG) analysis to constrain the physical conditions of the molecular gas. Acceptable models have significant degeneracies between parameters, even when we rule out all models that produce optically thin emission, but we find that the bulk of the molecular gas has T_kin=20-60 K, n_{H_2}~10^4-10^5 cm^-3, and N_CO/Delta-v=10^{17.00+/-0.25} cm^-2 km^-1 s. For our best-fit models to self-consistently recover a typical CO-to-H_2 abundance and a plausible degree of virialization, the local velocity gradient in the molecular gas must be substantially larger than its galaxy-wide average. This conclusion is consistent with a scenario in which SMM J14011+0252 has a fairly face-on orientation and a molecular ISM composed of many unresolved clouds. Using previous H-alpha observations, we find that SMM J14011+0252 has a spatially resolved star formation rate vs. molecular gas surface density relation inconsistent with those of "normal" local star-forming galaxies, even if we adopt a local "disk-like" CO-to-H_2 conversion factor as motivated by our LVG analysis. This discrepancy supports the inference of a star formation relation for high-z starbursts distinct from the local relation that is not solely due to differing choices of gas mass conversion factor.

We present a study using the Karl G. Jansky Very Large Array of 12CO J=1-0 emission in three strongly lensed submillimetre-selected galaxies (SMMJ16359, SMMJ14009 and SMMJ02399) at z=2.5-2.9. These galaxies span L(IR) = 10^11 - 10^13 Lsun, offering an opportunity to compare the interstellar medium of LIRGs and ULIRGs at high redshift. We estimate molecular gas masses in the range (2-40) x 10^9 Msun using a method that assumes canonical underlying brightness temperature ratios for star-forming and non-star-forming gas phases and a maximal star-formation efficiency. A more simplistic method using X(CO) = 0.8 yields gas masses twice as high. The observed CO(3-2)/CO(1-0) brightness temperature ratio for SMMJ14009, r(3-2)/(1-0) = (0.95 \pm 0.12), is indicative of warm star-forming gas, possibly influenced by the central AGN. We search for 12CO(1-0) emission in the Lyman-break galaxy, A2218 #384, located at z=2.517 in the same field as SMMJ16359, and assign a 3-sigma gas mass limit of <6 x 10^8 Msun. We use rest-frame 115-GHz free-free flux densities in SMMJ14009 and SMMJ02399 - measurements tied directly to the photionisation rate of massive stars and made possible by the VLA's bandwidth - to estimate star-formation rates of 400-600 Msun/yr and to estimate the fraction of L(IR) due to the AGN.

We report the detection of CO(1-0) emission in the strongly lensed high-redshift quasars IRAS F10214+4724 (z=2.286), the Cloverleaf (z=2.558), RX J0911+0551 (z=2.796), SMM J04135+10277 (z=2.846), and MG 0751+2716 (z=3.200), using the Expanded Very Large Array and the Green Bank Telescope. We report lensing-corrected CO(1-0) line luminosities of L'(CO) = 0.34-18.4 x 10^10 K km/s pc^2 and total molecular gas masses of M(H2) = 0.27-14.7 x 10^10 Msun for the sources in our sample. Based on CO line ratios relative to previously reported observations in J>=3 rotational transitions and line excitation modeling, we find that the CO(1-0) line strengths in our targets are consistent with single, highly-excited gas components with constant brightness temperature up to mid-J levels. We thus do not find any evidence for luminous extended, low excitation, low surface brightness molecular gas components. These properties are comparable to those found in z>4 quasars with existing CO(1-0) observations. These findings stand in contrast to recent CO(1-0) observations of z~2-4 submillimeter galaxies (SMGs), which have lower CO excitation and show evidence for multiple excitation components, including some low-excitation gas. These findings are consistent with the picture that gas-rich quasars and SMGs represent different stages in the early evolution of massive galaxies.

We report detections of three z ~ 2.5 submillimeter-selected galaxies (SMGs; SMM J14011+0252, SMM J14009+0252, SMM J04431+0210) in the lowest rotational transition of the carbon monoxide molecule (CO J = 1-0) and one nondetection (SMM J04433+0210). For the three galaxies we detected, we find a line-integrated brightness temperature ratio of the J = 3-2 and 1-0 lines of 0.68 +/- 0.08; the 1-0 line is stronger than predicted by the frequent assumption of equal brightnesses in the two lines and by most single-component models. The observed ratio suggests that mass estimates for SMGs based on J = 3-2 observations and J = 1-0 column density or mass conversion factors are low by a factor of 1.5. Comparison of the 1-0 line intensities with intensities of higher-J transitions indicates that single-component models for the interstellar media in SMGs are incomplete. The small dispersion in the ratio, along with published detections of CO lines with J_upper > 3 in most of the sources, indicates that the emission is from multi-component interstellar media with physical structures common to many classes of galaxies. This result tends to rule out the lowest scaling factors between CO luminosity and molecular gas mass, and further increases molecular mass estimates calibrated against observations of galaxies in the local universe. We also describe and demonstrate a statistically sound method for finding weak lines in broadband spectra that will find application in searches for molecular lines from sources at unknown redshifts.

We take advantage of good atmospheric transparency and the availability of high quality instrumentation in the 1 um near-infrared atmospheric window to present a grid of F, G, K, and M spectral standards observed at high spectral resolution (R ~ 25,000). In addition to a spectral atlas, we present a catalog of atomic line absorption features in the 0.95-1.11 um range. The catalog includes a wide range of line excitation potentials, from 0-13 eV, arising from neutral and singly ionized species, most frequently those of Fe I and Ti I at low excitation, Cr I, Fe I, and Si I at moderate excitation, and C I, S I, and Si I having relatively high excitation. The spectra also include several prominent molecular bands from CN and FeH. For the atomic species, we analyze trends in the excitation potential, line depth, and equivalent width across the grid of spectroscopic standards to identify temperature and surface gravity diagnostics near 1 um. We identify the line ratios that appear especially useful for spectral typing as those involving Ti I and C I or S I, which are temperature sensitive in opposite directions, and Sr II, which is gravity sensitive at all spectral types. Ascii versions of all spectra are available to download with the electronic version of the journal.

We combine Lyman-break colour selection with ultradeep (> 200 ks) Chandra X-ray imaging over a survey area of ~0.35 deg^2 to select high redshift AGN. Applying careful corrections for both the optical and X-ray selection functions, the data allow us to make the most accurate determination to date of the faint end of the X-ray luminosity function (XLF) at z~3. Our methodology recovers a number density of X-ray sources at this redshift which is at least as high as previous surveys, demonstrating that it is an effective way of selecting high z AGN. Comparing to results at z=1, we find no evidence that the faint slope of the XLF flattens at high z, but we do find significant (factor ~3.6) negative evolution of the space density of low luminosity AGN. Combining with bright end data from very wide surveys we also see marginal evidence for continued positive evolution of the characteristic break luminosity L*. Our data therefore support models of luminosity-dependent density evolution between z=1 and z=3. A sharp upturn in the the XLF is seen at the very lowest luminosities (Lx < 10^42.5 erg s^-1), most likely due to the contribution of pure X-ray starburst galaxies at very faint fluxes.