We compare the extent of the dust, molecular gas and stars in three star-forming galaxies, at $z= 1.4, 1.6$ and $2.7$, selected from the Hubble Ultra Deep Field based on their bright CO and dust-continuum emission as well as their large rest-frame optical sizes. The galaxies have high stellar masses, $\mathrm{M}_*>10^{11}\mathrm{M}_\odot$, and reside on, or slightly below, the main sequence of star-forming galaxies at their respective redshifts. We probe the dust and molecular gas using subarcsecond Atacama Large Millimeter/submillimeter Array observations of the 1.3 mm continuum and CO line emission, respectively, and probe the stellar distribution using \emph{Hubble Space Telescope} observations at 1.6 \textmu m. We find that for all three galaxies the CO emission appears $\gtrsim 30\%$ more compact than the stellar emission. For the $z= 1.4$ and $2.7$ galaxies, the dust emission is also more compact, by $\gtrsim 50\%$, than the stellar emission, whereas for the $z=1.6$ galaxy, the dust and stellar emission have similar spatial extents. This similar spatial extent is consistent with observations of local disk galaxies. However, most high redshift observations show more compact dust emission, likely due to the ubiquity of central starbursts at high redshift and the limited sensitivity of many of these observations. Using the CO emission line, we also investigate the kinematics of the cold interstellar medium in the galaxies, and find that all three have kinematics consistent with a rotation-dominated disk.

Using the NSF's Karl G. Jansky Very Large Array (VLA), we report six detections of CO(J=1-0) emission and one upper limit in z=2-3 galaxies originally detected in higher-J CO emission in the Atacama Large submillimeter/Millimeter Array (ALMA) Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). From the CO(J=1-0) line strengths, we measure total cold molecular gas masses of M_gas = 2.4-11.6 x 10^10 (alpha_CO/3.6) Msun. We also measure a median CO(J=3-2) to CO(J=1-0) line brightness temperature ratio of r_31 = 0.84 +/- 0.26, and a CO(J=7-6) to CO(J=1-0) ratio range of r_71 <0.05 to 0.17. These results suggest that CO(J=3-2) selected galaxies may have a higher CO line excitation on average than CO(J=1-0) selected galaxies, based on the limited, currently available samples from the ASPECS and VLA CO Luminosity Density at High Redshift (COLDz) surveys. This implies that previous estimates of the cosmic density of cold gas in galaxies based on CO(J=3-2) measurements should be revised down by a factor of ~=2 on average based on assumptions regarding CO excitation alone. This correction further improves the agreement between the best currently existing constraints on the cold gas density evolution across cosmic history from line scan surveys, and the implied characteristic gas depletion times.

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 a sub-kpc resolved study of the interstellar medium properties in SDP.81, a z=3.042 strongly gravitationally lensed dusty star-forming galaxy, based on high-resolution, multi-band ALMA observations of the FIR continuum, CO ladder and the [CII] line. Using a visibility-plane lens modelling code, we achieve a median source-plane resolution of ~200 pc. We use photon-dominated region (PDR) models to infer the physical conditions - far-UV field strength, density, and PDR surface temperature - of the star-forming gas on 200-pc scales, finding a FUV field strength of ~10^3-10^4 G0, gas density of ~10^5 cm^-3 and cloud surface temperatures up to 1500 K, similar to those in the Orion Trapezium region. The [CII] emission is significantly more extended than that FIR continuum: ~50 per cent of [CII] emission arises outside the FIR-bright region. The resolved [CII]/FIR ratio varies by almost 2 dex across the source, down to ~2x10^-4 in the star-forming clumps. The observed [CII]/FIR deficit trend is consistent with thermal saturation of the C+ fine-structure level occupancy at high gas temperatures. We make the source-plane reconstructions of all emission lines and continuum data publicly available.

The Atacama Large Millimetre/submillimetre Array (ALMA) is currently in the process of transforming our view of star-forming galaxies in the distant ($z\gtrsim1$) universe. Before ALMA, most of what we knew about dust-obscured star formation in distant galaxies was limited to the brightest submillimetre sources -- the so-called submillimetre galaxies (SMGs) -- and even the information on those sources was sparse, with resolved (i.e., sub-galactic) observations of the obscured star formation and gas reservoirs typically restricted to the most extreme and/or strongly lensed sources. Starting with the beginning of early science operations in 2011, the last nine years of ALMA observations have ushered in a new era for studies of high-redshift star formation. With its long baselines, ALMA has allowed observations of distant dust-obscured star formation with angular resolutions comparable to -- or even far surpassing -- the best current optical telescopes. With its bandwidth and frequency coverage, it has provided an unprecedented look at the associated molecular and atomic gas in these distant galaxies through targeted follow-up and serendipitous detections/blind line scans. Finally, with its leap in sensitivity compared to previous (sub-)millimetre arrays, it has enabled the detection of these powerful dust/gas tracers much further down the luminosity function through both statistical studies of color/mass-selected galaxy populations and dedicated deep fields. We review the main advances ALMA has helped bring about in our understanding of the dust and gas properties of high-redshift ($z\gtrsim1$) star-forming galaxies during these first nine years of its science operations, and we highlight the interesting questions that may be answered by ALMA in the years to come.

We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K ($\beta=1.9$) derived under the common assumption of an optically thin far-infrared dust emission, but fully consistent with $T_{\rm d}=52\pm5$ K of a general opacity model where the optical depth ($\tau$) reaches unity at a wavelength of $\lambda_0=170\pm23$ $\mu$m. Moreover, the general opacity solution returns a factor of $\sim 2\times$ lower dust mass and, hence, a lower molecular gas mass for a fixed gas-to-dust ratio, than with the optically thin dust model. The derived properties of GN20 thus provide an appealing solution to the puzzling discovery of starbursts appearing colder than main-sequence galaxies above $z>2.5$, in addition to a lower dust-to-stellar mass ratio that approaches the physical value predicted for starburst galaxies.

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 report a ground-based detection of the [OI] 63-um line in a z=6.027 gravitationally lensed dusty star-forming galaxy (DSFG) G09.83808 using the APEX SEPIA 660 receiver, the first unambiguous detection of the [OI]63 line beyond redshift 3, and the first obtained from the ground. The [OI]63 line is robustly detected at 22$\pm$5 Jy km s$^{-1}$, corresponding to an intrinsic (de-lensed) luminosity of $(5.4\pm1.3)\times10^{9}$ L$_\odot$. With the [OI]63/[CII] luminosity ratio of 4, the [OI]63 line is the main coolant of the neutral gas in this galaxy, in agreement with model predictions. The high [OI]63 luminosity compensates for the pronounced [CII] deficit ([CII]/FIR$\simeq4\times10^{-4}$). Using photon-dominated region models, we derive a source-averaged gas density $n=10^{4.0}$ cm$^{-3}$, and far-UV field strength $G=10^4 G_0$, comparable to the z=2-4 DSFG population. If G09.83808 represents a typical high-redshift DSFG, the [OI]63 line from z=6 non-lensed DSFGs should be routinely detectable in ALMA Band 9 observations with $\sim$15 min on-source, opening a new window to study the properties of the earliest DSFGs.

We analyse the physical properties of a large, homogeneously selected sample of ALMA-located sub-mm galaxies (SMGs) detected in the SCUBA-2 CLS 850-$\mu$m map of the UKIDSS/UDS field. This survey, AS2UDS, identified 707 SMGs across the ~1 sq.deg. field, including ~17 per cent which are undetected in the optical/near-infrared to $K$>~25.7 mag. We interpret the UV-to-radio data using a physically motivated model, MAGPHYS and determine a median photometric redshift of z=2.61+-0.08, with a 68th percentile range of z=1.8-3.4 and just ~6 per cent at z>4. The redshift distribution is well fit by a model combining evolution of the gas fraction in halos with the growth of halo mass past a threshold of ~4x10$^{12}$M$_\odot$, thus SMGs may represent the highly efficient collapse of gas-rich massive halos. Our survey provides a sample of the most massive, dusty galaxies at z>~1, with median dust and stellar masses of $M_d$=(6.8+-0.3)x10$^{8}$M$_\odot$ (thus, gas masses of ~10$^{11}$M$_\odot$) and $M_\ast=$(1.26+-0.05)x10$^{11}$M$_\odot$. These galaxies have gas fractions of $f_{gas}=$0.41+-0.02 with depletion timescales of ~150Myr. The gas mass function evolution at high masses is consistent with constraints at lower masses from blind CO-surveys, with an increase to z~2-3 and then a decline at higher redshifts. The space density and masses of SMGs suggests that almost all galaxies with $M_\ast$>~2x10$^{11}$M$_\odot$ have passed through an SMG-like phase. We find no evolution in dust temperature at a constant far-infrared luminosity across z~1.5-4. We show that SMGs appear to behave as simple homologous systems in the far-infrared, having properties consistent with a centrally illuminated starburst. Our study provides strong support for an evolutionary link between the active, gas-rich SMG population at z>1 and the formation of massive, bulge-dominated galaxies across the history of the Universe.

We analyze the interstellar medium properties of a sample of sixteen 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 a couple of additional CO line emitters in the UDF that are identified based on their MUSE optical spectroscopic redshifts. The ASPECS CO$-$selected galaxies cover a larger range of star-formation rates 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 star-formation rate-stellar mass relation (`galaxy main sequence'). However, we find that $\sim30\%$ of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with $\sim12\%$ (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low star-formation rates, 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 contributor to the molecular gas density at any redshift probed by our observations (z$\sim$1$-$3). The respective contribution by starburst galaxies above the main sequence decreases from z$\sim$2.5 to z$\sim$1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequence.

We combine high-resolution ALMA and HST/CANDELS observations of 20 submillimeter galaxies (SMGs) predominantly from the AS2UDS survey at z~2 with bright rest-frame optical counterparts (Ks < 22.9) to investigate the resolved structural properties of their dust and stellar components. We derive two-dimensional stellar-mass distributions that are inferred from spatial mass-to-light ratio (M/L) corrections based on rest-frame optical colors. Due to the high central column densities of dust in our SMGs, our mass distributions likely represent a lower limit to the true central mass density. The centroid positions between the inferred stellar-mass and the dust distributions agree within 1.1 kpc, indicating an overall good spatial agreement between the two components. The majority of our sources exhibit compact dust configurations relative to the stellar component (with a median ratio of effective radii Re,dust/Re,Mstar = 0.6). This ratio does not change with specific star-formation rate (sSFR) over the factor of 30 spanned by our targets, sampling the locus of "normal" main sequence galaxies up to the starburst regime, log(sSFR/sSFRMS) > 0.5. Our results imply that massive SMGs are experiencing centrally enhanced star formation unlike typical spiral galaxies in the local Universe. The sizes and stellar densities of our SMGs are in agreement with those of the passive population at z=1.5, consistent with these systems being the descendants of z~2 SMGs.

We present sub-kpc-scale mapping of the 870 $\mu$m ALMA continuum emission in six luminous ($L_{\rm IR}~\sim~5~\times10^{12}$ L$_{\odot}$) submillimeter galaxies (SMGs) from the ALESS survey of the Extended Chandra Deep Field South. Our high-fidelity 0.07$''$-resolution imaging ($\sim$500 pc) reveals robust evidence for structures with deconvolved sizes of $\lesssim$0.5-1 kpc embedded within (dominant) exponential dust disks. The large-scale morphologies of the structures within some of the galaxies show clear curvature and/or clump-like structures bracketing elongated nuclear emission, suggestive of bars, star-forming rings, and spiral arms. In this interpretation, the ratio of the `ring' and `bar' radii (1.9$\pm$0.3) agrees with that measured for such features in local galaxies. These potential spiral/ring/bar structures would be consistent with the idea of tidal disturbances, with their detailed properties implying flat inner rotation curves and Toomre-unstable disks (Q<1). The inferred one-dimensional velocity dispersions ($\sigma_{\rm r}\lesssim$ 70-160 km s$^{-1}$) are marginally consistent with the limits implied if the sizes of the largest structures are comparable to the Jeans length. We create maps of the star formation rate density ($\Sigma_{\rm SFR}$) on $\sim$500 pc scales and show that the SMGs are able to sustain a given (galaxy-averaged) $\Sigma_{\rm SFR}$ over much larger physical scales than local (ultra-)luminous infrared galaxies. However, on 500 pc scales, they do not exceed the Eddington limit set by radiation pressure on dust. If confirmed by kinematics, the potential presence of non-axisymmetric structures would provide a means for net angular momentum loss and efficient star formation, helping to explain the very high star formation rates measured in SMGs.

The ALMA SPECtroscopic Survey in the {\it Hubble} Ultra Deep Field 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. We here 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 sixteen emission lines at high fidelity in our search. Comparing these sources to multi-wavelength data we find that all sources have optical/infrared counterparts. Our search also recovers candidates that have 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-mJy 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 $\beta=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 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 $\rho$(H$_2$). 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 $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$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-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm 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 $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history).

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 MUSE HUDF Survey and multi-wavelength data to uniquely associate all 16 line-emitters, detected in the ALMA data without preselection, with rotational transitions of carbon monoxide (CO). We identify ten as CO(2-1) at $1 < z < 2$, five as CO(3-2) at $2 < z < 3$ and one 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 \le 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 \sim 1.4$ and $z \sim 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 \sim 1.4$. For stellar masses $\ge 10^{10} (10^{10.5})$ M$_{\odot}$, 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 on MS galaxies during the peak of galaxy formation.

One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it is converted to stars. Over the last decade, important progress has been made in understanding the relative importance of these two factors along with the role of feedback, and the first measurements of the volume density of cold gas out to redshift 4, (the "cold gas history of the Universe") has been obtained. To match the precision of measurements of the star formation and black-hole accretion histories over the coming decades, a two orders of magnitude improvement in molecular line survey speeds is required compared to what is possible with current facilities. Possible pathways towards such large gains include significant upgrades to current facilities like ALMA by 2030 (and beyond), and eventually the construction of a new generation of radio-to-millimeter wavelength facilities, such as the next generation Very Large Array (ngVLA) concept.

Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions.

We present the catalogue and properties of sources in AS2UDS, an 870-$\mu$m continuum survey with the Atacama Large Millimetre/sub-millimetre Array (ALMA) of 716 single-dish sub-millimetre sources detected in the UKIDSS/UDS field by the SCUBA-2 Cosmology Legacy Survey. In our sensitive ALMA follow-up observations we detect 708 sub-millimetre galaxies (SMGs) at $>$\,4.3$\sigma$ significance across the $\sim$\,1-degree diameter field. We combine our precise ALMA positions with the extensive multi-wavelength coverage in the UDS field to fit the spectral energy distributions of our SMGs to derive a median redshift of $z_{\rm phot}=$\,2.61$\pm$0.09. This large sample reveals a statistically significant trend of increasing sub-millimetre flux with redshift suggestive of galaxy downsizing. 101 ALMA maps do not show a $>$\,4.3$\sigma$ SMG, but we demonstrate from stacking {\it Herschel} SPIRE observations at these positions, that the vast majority of these blank maps correspond to real single-dish sub-millimetre sources. We further show that these blank maps contain an excess of galaxies at $z_{\rm phot}=$\,1.5--4 compared to random fields, similar to the redshift range of the ALMA-detected SMGs. In addition, we combine X-ray and mid-infrared active galaxy nuclei activity (AGN) indicators to yield a likely range for the AGN fraction of 8--28\,\% in our sample. Finally, we compare the redshifts of this population of high-redshift, strongly star-forming galaxies with the inferred formation redshifts of massive, passive galaxies being found out to $z\sim$\,2, finding reasonable agreement -- in support of an evolutionary connection between these two classes of massive galaxy.

Recent observations of lensed sources have shown that the faintest ($M_{\mathrm{UV}} \approx -15\,\mathrm{mag}$) galaxies observed at z=6-8 appear to be extremely compact. Some of them have inferred sizes of less than 40 pc for stellar masses between $10^6$ and $10^7\,\mathrm{M}_{\odot}$, comparable to individual super star clusters or star cluster complexes at low redshift. High-redshift, low-mass galaxies are expected to show a clumpy, irregular morphology and if star clusters form in each of these well-separated clumps, the observed galaxy size would be much larger than the size of an individual star forming region. As supernova explosions impact the galaxy with a minimum delay time that exceeds the time required to form a massive star cluster, other processes are required to explain the absence of additional massive star forming regions. In this work we investigate whether the radiation of a young massive star cluster can suppress the formation of other detectable clusters within the same galaxy already before supernova feedback can affect the galaxy. We find that in low-mass ($M_{200} \lesssim 10^{10}\,\mathrm{M}_{\odot}$) haloes, the radiation from a compact star forming region with an initial mass of $10^{7}\,\mathrm{M}_{\odot}$ can keep gas clumps with Jeans masses larger than $\approx 10^{7}\,\mathrm{M}_{\odot}$ warm and ionized throughout the galaxy. In this picture, the small intrinsic sizes measured in the faintest $z=6-8$ galaxies are a natural consequence of the strong radiation field that stabilises massive gas clumps. A prediction of this mechanism is that the escape fraction for ionizing radiation is high for the extremely compact, high-z sources.

We present 0.15-arcsec (1 kpc) resolution ALMA observations of the [CII] 157.74 um line and rest-frame 160-um continuum emission in two z~3 dusty, star-forming galaxies - ALESS 49.1 and ALESS 57.1, combined with resolved CO(3-2) observations. In both sources, the [CII] surface brightness distribution is dominated by a compact core $\leq$1 kpc in radius, a factor of 2-3 smaller than the extent of the CO(3-2) emission. In ALESS 49.1, we find an additional extended (8-kpc radius), low surface-brightness [CII] component. Based on an analysis of mock ALMA observations, the [CII] and 160-um continuum surface brightness distributions are inconsistent with a single-Gaussian surface brightness distribution with the same size as the CO(3-2) emission. The [CII] rotation curves flatten at $\simeq$2 kpc radius, suggesting the kinematics of the central regions are dominated by a baryonic disc. Both galaxies exhibit a strong [CII]/FIR deficit on 1-kpc scales, with FIR-surface-brightness to [CII]/FIR slope steeper than in local star-forming galaxies. A comparison of the [CII]/CO(3-2) observations with PDR models suggests a strong FUV radiation field ($G_0\sim10^4$) and high gas density ($n\mathrm{(H)}\sim10^4-10^5$ cm$^{-3}$) in the central regions of ALESS 49.1 and 57.1. The most direct interpretation of the pronounced [CII]/FIR deficit is a thermal saturation of the C+ fine-structure levels at temperatures $\geq$500 K, driven by the strong FUV field.

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 carry out a blind search of 3mm continuum sources using the ALMA Science Archive to derive the first galaxy number counts at this wavelength. The analyzed data are drawn from observations towards three extragalactic legacy fields: COSMOS, CDF-S, and the UDS comprising more than 130 individual ALMA Band 3 pointings and an effective survey area of ~200 sq. arcmin with a continuum sensitivity that allows for the direct detection of unlensed Dusty Star-Forming Galaxies (DSFGs) dust emission beyond the epoch of reionization. We present a catalog of 16 sources detected at >5sigma with flux densities S_3mm = 60 - 600 uJy from which number counts are derived. These number counts are then used to place constraints on the volume density of DSFGs with an empirical backward evolution model. Our measured 3mm number counts indicate that the contribution of DSFGs to the cosmic star formation rate density at z >~ 4 is non-negligible. This is contrary to the generally adopted assumption of a sharply decreasing contribution of obscured galaxies at z > 4 as inferred by optical and near-infrared surveys. This work demonstrates the power of ALMA 3mm observations which can reach outstanding continuum sensitivities during typical spectral line science programs. Further constraints on 3mm-selected galaxies will be essential to refine models of galaxy formation and evolution as well as models of early Universe dust production mechanisms.

The goal of this science case is to address the use of a ngVLA as a CO redshift machine for dust-obscured high-redshift galaxies which lack of clear counterparts at other wavelengths. Thanks to its unprecedentedly large simultaneous bandwidth and sensitivity, the ngVLA will be able to detect low--J CO transitions at virtually any $z>1$. In particular, at $z>4.76$ two CO transitions will be covered in a single frequency setting, thus ensuring unambiguous line identification. The ngVLA capabilities fill in a redshift range where other approaches (e.g., photometric redshifts, search for optical/radio counterparts, etc) typically fail due to the combination of intrinsically faint emission and increasing luminosity distance. This will allow us to explore the formation of massive galaxies in the early cosmic times.

The goal of this science case is to accurately pin down the molecular gas content of high redshift galaxies. By targeting the CO ground transition, we circumvent uncertainties related to CO excitation. The ngVLA can observe the CO(1-0) line at virtually any $z>1.5$, thus exposing the evolution of gaseous reservoirs from the earliest epochs down to the peak of the cosmic history of star formation. The order-of-magnitude improvement in the number of CO detections with respect to state-of-the-art observational campaigns will provide a unique insight on the evolution of galaxies through cosmic time.

The goal of this science case is to study physical conditions of the interstellar medium (ISM) in distant galaxies. In particular, its densest component is associated with the inner cores of clouds -- this is where star formation takes place. Carbon monoxide is usually used to trace molecular gas emission; however, its transitions are practically opaque, thus preventing astronomers from piercing through the clouds, into the deepest layers that are most intimately connected with the formation of stars. Other dense gas tracers are required, although they are typically too faint and/or at too low frequencies to be effectively observed in high redshift galaxies. The ngVLA will offer for the first time the sensitivity at radio frequencies that is needed to target [CI]$_{1-0}$ (at $z>5$), as well as the ground transitions of dense gas tracers of the ISM such as HCN, HNC, HCO+ (at various redshifts $z>1$), beyond the tip of the iceberg of the hyper-luminous sources that could be studied up to now. These new tools will critically contribute to our understanding of the intimate interplay between gas clouds and star formation in different environments and cosmic epochs.

The interstellar medium feeds both the formation of stars and the growth of black holes, making it a key ingredient in the evolution of galaxies. With the advent of the Atacama Large Millimeter/ submillimeter Array (ALMA), we can now probe the interstellar medium within high-redshift galaxies in increasingly exquisite detail. Our recent ALMA observations map the molecular gas and dust continuum emission in sub-millimetre-selected galaxies on 1-5 kpc scales, revealing significant differences in how the gas, dust continuum, and existing stellar emission are distributed within the galaxies. This study demonstrates the power of ALMA to shed new light on the structure and kinematics of the interstellar medium in the early Universe, suggesting that the interpretation of such observations is more complex than typically assumed.

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.

Deep, pencil-beam surveys from ALMA at 1.1-1.3mm have uncovered an apparent absence of high-redshift dusty galaxies, with existing redshift distributions peaking around $z\sim1.5-2.5$. This has led to a perceived dearth of dusty systems at $z>4$, and the conclusion, according to some models, that the early Universe was relatively dust-poor. In this paper, we extend the backward evolution galaxy model described by Casey et al. (2018) to the ALMA regime (in depth and area) and determine that the measured number counts and redshift distributions from ALMA deep field surveys are fully consistent with constraints of the infrared luminosity function (IRLF) at $z<2.5$ determined by single-dish submillimeter and millimeter surveys conducted on much larger angular scales ($\sim1-10$deg$^{2}$). We find that measured 1.1-1.3mm number counts are most constraining for the measurement of the faint-end slope of the IRLF at $z<2.5$ instead of the prevalence of dusty galaxies at $z>4$. Recent studies have suggested that UV-selected galaxies at $z>4$ may be particularly dust-poor, but we find their millimeter-wave emission cannot rule out consistency with the Calzetti dust attenuation law even by assuming relatively typical, cold-dust ($T_{\rm dust}\approx30\,$K) SEDs. Our models suggest that the design of ALMA deep fields requires substantial revision to constrain the prevalence of $z>4$ early Universe obscured starbursts. The most promising avenue for detection and characterization of such early dusty galaxies will come from future ALMA 2mm blank field surveys covering a few hundred arcmin$^{2}$ and the combination of existing and future dual-purpose 3mm datasets.

We present ALMA observations of the mid-J 12CO emission from six single-dish selected 870-micron sources in the Extended Chandra Deep Field-South (ECDFS) and UKIDSS Ultra-Deep Survey (UDS) fields. These six single-dish submillimetre sources were selected based on previous ALMA continuum observations, which showed that each comprised a blend of emission from two or more individual submillimetre galaxies (SMGs), separated on 5--10 arcsec scales. The six single-dish submillimetre sources targeted correspond to a total of 14 individual SMGs, of which seven have previously-measured robust optical/near-infrared spectroscopic redshifts, which were used to tune our ALMA observations. We detect CO(3-2) or CO(4-3) at z=2.3--3.7 in seven of the 14 SMGs, and in addition serendipitously detect line emission from three gas-rich companion galaxies, as well as identify four new 3.3-mm selected continuum sources in the six fields. Joint analysis of our CO spectroscopy and existing data suggests that 64 \pm 18% of the SMGs in blended submillimetre sources are unlikely to be physically associated. However, three of the SMG fields (50%) contain new, serendipitously-detected CO-emitting (but submillimetre-faint) sources at similar redshifts to the 870-micron selected SMGs we targeted. These data suggest that the SMGs inhabit overdense regions, but that these are not sufficiently overdense on ~100 kpc scales to influence the source blending given the short lifetimes of SMGs. We find that 21 \pm 12% of SMGs have spatially-distinct and kinematically-close companion galaxies (~8--150 kpc and <~300 km/s), which may have enhanced their star-formation via gravitational interactions.

Though half of cosmic starlight is absorbed by dust and reradiated at long wavelengths (3$\mu$m-3mm), constraints on the infrared through millimeter galaxy luminosity function (the `IRLF') are poor in comparison to the rest-frame ultraviolet and optical galaxy luminosity function, particularly at z>2.5. Here we present a backward evolution model for interpreting number counts, redshift distributions, and cross-band flux density correlations in the infrared and submillimeter sky, from 70$\mu$m-2mm, using a model for the IRLF out to the epoch of reionization. Mock submillimeter maps are generated by injecting sources according to the prescribed IRLF and flux densities drawn from model spectral energy distributions that mirror the distribution of SEDs observed in $0<z<5$ dusty star-forming galaxies (DSFGs). We explore two extreme hypothetical case-studies: a dust-poor early Universe model, where DSFGs contribute negligibly ($<$10%) to the integrated star-formation rate density at $z>4$, and an alternate dust-rich early Universe model, where DSFGs dominate $\sim$90% of $z>4$ star-formation. We find that current submm/mm datasets do not clearly rule out either of these extreme models. We suggest that future surveys at 2mm will be crucial to measuring the IRLF beyond $z\sim4$. The model framework developed in this paper serves as a unique tool for the interpretation of multiwavelength IR/submm extragalactic datasets and will enable more refined constraints on the IRLF than can be made from direct measurements of individual galaxies' integrated dust emission.

We use ALMA observations of four sub-millimetre galaxies (SMGs) at $z\sim2-3$ to investigate the spatially resolved properties of the inter-stellar medium (ISM) at scales of 1--5 kpc (0.1--0.6$''$). The velocity fields of our sources, traced by the $^{12}$CO($J$=3-2) emission, are consistent with disk rotation to first order, implying average dynamical masses of $\sim$3$\times10^{11}$M$_{\odot}$ within two half-light radii. Through a Bayesian approach we investigate the uncertainties inherent to dynamically constraining total gas masses. We explore the covariance between the stellar mass-to-light ratio and CO-to-H$_{2}$ conversion factor, $\alpha_{\rm CO}$, finding values of $\alpha_{\rm CO}=1.1^{+0.8}_{-0.7}$ for dark matter fractions of 15 \%. We show that the resolved spatial distribution of the gas and dust continuum can be uncorrelated to the stellar emission, challenging energy balance assumptions in global SED fitting. Through a stacking analysis of the resolved radial profiles of the CO(3-2), stellar and dust continuum emission in SMG samples, we find that the cool molecular gas emission in these sources (radii $\sim$5--14 kpc) is clearly more extended than the rest-frame $\sim$250 $\mu$m dust continuum by a factor $>2$. We propose that assuming a constant dust-to-gas ratio, this apparent difference in sizes can be explained by temperature and optical-depth gradients alone. Our results suggest that caution must be exercised when extrapolating morphological properties of dust continuum observations to conclusions about the molecular gas phase of the ISM.

We present deep high resolution (0.03", 200pc) ALMA Band 7 observations covering the dust continuum and [CII] $\lambda157.7\mu$m emission in four $z\sim4.4-4.8$ sub-millimeter galaxies (SMGs) selected from the ALESS and AS2UDS surveys. The data show that the rest-frame 160$\mu$m (observed 345 GHz) dust emission is consistent with smooth morphologies on kpc scales for three of the sources. One source, UDS47.0, displays apparent substructure but this is also consistent with a smooth morphology, as indicated by simulations showing that smooth exponential disks can appear clumpy when observed at high angular resolution (0.03") and depth of these observations ($\sigma_{345\text{GHz}} \sim27-47\mu$Jy beam$^{-1}$). The four SMGs are bright [CII] emitters, and we extract [CII] spectra from the high resolution data, and recover $\sim20-100$% of the [CII] flux and $\sim40-80$% of the dust continuum emission, compared to the previous lower resolution observations. When tapered to 0.2" resolution our maps recover $\sim80-100$% of the continuum emission, indicating that $\sim60$% of the emission is resolved out on $\sim200$pc scales. We find that the [CII] emission in high-redshift galaxies is more spatially extended than the rest-frame 160$\mu$m dust continuum by a factor of $1.6\pm0.4$. By considering the $L_{\text{[CII]}}$/$L_{\text{FIR}}$ ratio as a function of the star-formation rate surface density ($\Sigma_{\text{SFR}}$) we revisit the [CII] deficit, and suggest that the decline in the $L_{\text{[CII]}}$/$L_{\text{FIR}}$ ratio as a function of $\Sigma_{\text{SFR}}$ is consistent with local processes. We also explore the physical drivers that may be responsible for these trends and can give rise to the properties found in the densest regions of SMGs.

The earliest galaxies are expected to emerge in the first billion years of the Universe during the Epoch of Reionization. However, both the spectroscopic confirmation of photometrically-selected galaxies at this epoch and the characterization of their early dynamical state has been hindered by the lack of bright, accessible lines to probe the velocity structure of their interstellar medium. We present the first ALMA spectroscopic confirmation of such sources at z > 6 using the far-infrared [C II]{\lambda}157.74{\mu}m emission line, and, for the first time, measurement of the velocity structure, for two galaxies at z = 6.8540+/-0.0003 and z = 6.8076+/-0.0002. Remarkably, the [C II] line luminosity from these galaxies is higher than previously found in `normal' star-forming galaxies at z > 6.5. This suggests that we are sampling a part of the galaxy population different from the galaxies found through detection of the Ly{\alpha} line. The luminous and extended [C II] detections reveal clear velocity gradients that, if interpreted as rotation, would suggest these galaxies have turbulent, yet rotation-dominated disks, with similar stellar-to-dynamical mass fractions as observed for H{\alpha} emitting galaxies 2 Gyr later at cosmic noon. Our novel approach for confirming galaxies during Reionization paves the way for larger studies of distant galaxies with spectroscopic redshifts from ALMA. Particularly important, this opens up opportunities for high angular-resolution [C II] dynamics in galaxies less than one billion years after the Big Bang.

We present an extremely deep CO(1-0) observation of a confirmed $z=1.62$ galaxy cluster. We detect two spectroscopically confirmed cluster members in CO(1-0) with $S/N>5$. Both galaxies have log(${\cal M_{\star}}$/\msol)$>11$ and are gas rich, with ${\cal M}_{\rm mol}$/(${\cal M_{\star}}+{\cal M}_{\rm mol}$)$\sim 0.17-0.45$. One of these galaxies lies on the star formation rate (SFR)-${\cal M_{\star}}$ sequence while the other lies an order of magnitude below. We compare the cluster galaxies to other SFR-selected galaxies with CO measurements and find that they have CO luminosities consistent with expectations given their infrared luminosities. We also find that they have comparable gas fractions and star formation efficiencies (SFE) to what is expected from published field galaxy scaling relations. The galaxies are compact in their stellar light distribution, at the extreme end for all high redshift star-forming galaxies. However, their SFE is consistent with other field galaxies at comparable compactness. This is similar to two other sources selected in a blind CO survey of the HDF-N. Despite living in a highly quenched proto-cluster core, the molecular gas properties of these two galaxies, one of which may be in the processes of quenching, appear entirely consistent with field scaling relations between the molecular gas content, stellar mass, star formation rate, and redshift. We speculate that these cluster galaxies cannot have any further substantive gas accretion if they are to become members of the dominant passive population in $z<1$ clusters.

We present detailed studies of a $z=2.12$ submillimeter galaxy, ALESS67.1, using sub-arcsecond resolution ALMA, AO-aided VLT/SINFONI, and HST/CANDELS data to investigate the kinematics and spatial distributions of dust emission (870 $\mu$m continuum), $^{12}$CO($J$=3-2), strong optical emission lines, and visible stars. Dynamical modelling of the optical emission lines suggests that ALESS67.1 is not a pure rotating disk but a merger, consistent with the apparent tidal features revealed in the HST imaging. Our sub-arcsecond resolution dataset allow us to measure half-light radii for all the tracers, and we find a factor of 4-6 smaller sizes in dust continuum compared to all the other tracers, including $^{12}$CO, and UV and H$\alpha$ emission is significantly offset from the dust continuum. The spatial mismatch between UV continuum and the cold dust and gas reservoir supports the explanation that geometrical effects are responsible for the offset of dusty galaxy on the IRX-$\beta$ diagram. Using a dynamical method we derive an $\alpha_{\rm CO}=1.8\pm1.0$, consistent with other SMGs that also have resolved CO and dust measurements. Assuming a single $\alpha_{\rm CO}$ value we also derive resolved gas and star-formation rate surface densities, and find that the core region of the galaxy ($\lesssim5$ kpc) follows the trend of mergers on the Schmidt-Kennicutt relationship, whereas the outskirts ($\gtrsim5$ kpc) lie on the locus of normal star-forming galaxies, suggesting different star-formation efficiencies within one galaxy. Our results caution against using single size or morphology for different tracers of the star-formation activity and gas content of galaxies, and therefore argue the need to use spatially-resolved, multi-wavelength observations to interpret the properties of SMGs, and perhaps even for $z>1$ galaxies in general.

We present single-dish CO(1-0) and CO(2-1) observations for 14 low-redshift quasi-stellar objects (QSOs). In combination with optical integral field spectroscopy we study how the cold gas content relates to the star formation rate (SFR) and black hole accretion rate. CO(1-0) is detected in 8 of 14 targets and CO(2-1) is detected in 7 out of 11 cases. The majority of disc-dominated QSOs reveal gas fractions and depletion times well matching normal star forming systems. Two gas-rich major mergers show clear starburst signatures with higher than average gas fractions and shorter depletion times. Bulge-dominated QSO hosts are mainly undetected in CO(1-0) which corresponds, on average, to lower gas fractions than in disc-dominated counterparts. Their SFRs however imply shorter than average depletion times and higher star formation efficiencies. Negative QSO feedback through removal of cold gas seems to play a negligible role in our sample. We find a trend between black hole accretion rate and total molecular gas content for disc-dominated QSOs when combined with literature samples. We interpret this as an upper envelope for nuclear activity which is well represented by a scaling relation between the total and circum-nuclear gas reservoir accessible for accretion. Bulge-dominated QSOs significantly differ from that scaling relation and appear uncorrelated with the total molecular gas content. This could be explained either by a more compact gas reservoir, blow out of the gas envelope through outflows, or a different ISM phase composition.

We present spectroscopic redshifts of S(870)>2mJy submillimetre galaxies (SMGs) which have been identified from the ALMA follow-up observations of 870um detected sources in the Extended Chandra Deep Field South (the ALMA-LESS survey). We derive spectroscopic redshifts for 52 SMGs, with a median of z=2.4+/-0.1. However, the distribution features a high redshift tail, with ~25% of the SMGs at z>3. Spectral diagnostics suggest that the SMGs are young starbursts, and the velocity offsets between the nebular emission and UV ISM absorption lines suggest that many are driving winds, with velocity offsets up to 2000km/s. Using the spectroscopic redshifts and the extensive UV-to-radio photometry in this field, we produce optimised spectral energy distributions (SEDs) using Magphys, and use the SEDs to infer a median stellar mass of M*=(6+/-1)x10^{10}Msol for our SMGs with spectroscopic redshifts. By combining these stellar masses with the star-formation rates (measured from the far-infrared SEDs), we show that SMGs (on average) lie a factor ~5 above the main-sequence at z~2. We provide this library of 52 template fits with robust and well-sampled SEDs available as a resource for future studies of SMGs, and also release the spectroscopic catalog of ~2000 (mostly infrared-selected) galaxies targeted as part of the spectroscopic campaign.

We present ALMA observations of the [CII] fine structure line and the underlying far-infrared (FIR) dust continuum emission in J1120+0641, the most distant quasar currently known (z=7.1). We also present observations targeting the CO(2-1), CO(7-6) and [CI] 369 micron lines in the same source obtained at the VLA and PdBI. We find a [CII] line flux of F_[CII]=1.11+/-0.10 Jy km/s and a continuum flux density of S_227GHz=0.53+/-0.04 mJy/beam, consistent with previous unresolved measurements. No other source is detected in continuum or [CII] emission in the field covered by ALMA (~25"). At the resolution of our ALMA observations (0.23", or 1.2 kpc, a factor ~70 smaller beam area compared to previous measurements), we find that the majority of the emission is very compact: a high fraction (~80%) of the total line and continuum flux is associated with a region 1-1.5 kpc in diameter. The remaining ~20% of the emission is distributed over a larger area with radius <4 kpc. The [CII] emission does not exhibit ordered motion on kpc-scales: applying the virial theorem yields an upper limit on the dynamical mass of the host galaxy of (4.3+/-0.9)x10^10 M_sun, only ~20x higher than the central black hole. The other targeted lines (CO(2-1), CO(7-6) and [CI]) are not detected, but the limits of the line ratios with respect to the [CII] emission imply that the heating in the quasar host is dominated by star formation, and not by the accreting black hole. The star-formation rate implied by the FIR continuum is 105-340 M_sun/yr, with a resulting star-formation rate surface density of ~100-350 M_sun/yr/kpc^2, well below the value for Eddington-accretion-limited star formation.

We present the first results from our on-going Australia Telescope Compact Array survey of CO(1-0) in ALMA-identified submillimetre galaxies in the Extended Chandra Deep Field South. Strong detections of CO(1-0) emission from two submillimetre galaxies, ALESS 122.1 (z = 2.0232) and ALESS 67.1 (z = 2.1230), were obtained. We estimate gas masses of M_gas ~ 1.3 \times 10^{11} M_odot and M_gas ~ 1.0 \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively, adopting alpha_CO = 1.0. Dynamical mass estimates from the kinematics of the CO(1-0) line yields M_dyn (sin i)^2 = 2.1 +- 1.1 \times 10^{11} M_odot and (3.2 +- 0.9) \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively. This is consistent with the total baryonic mass estimates of these two systems. We examine star formation efficiency using the L_FIR versus L'_CO(1-0) relation for samples of local ULIRGs and LIRGs, and more distant star-forming galaxies, with CO(1-0) detections. We find some evidence of a shallower slope for ULIRGs and SMGs compared to less luminous systems, but a larger sample is required for definite conclusions. We determine gas-to-dust ratios of 170 +- 30 and 140 +- 30 for ALESS 122.1 and ALESS 67.1, respectively, showing ALESS 122.1 has an unusually large gas reservoir. By combining the 38.1 GHz continuum detection of ALESS 122.1 with 1.4 and 5.5 GHz data, we estimate that the free-free contribution to radio emission at 38.1 GHz is 34 +- 17 microJy, yielding a star formation rate (1400 +- 700 M_\odot yr^{-1}) consistent with that from the infrared luminosity.

We present a detailed study of a molecular outflow feature in the nearby starburst galaxy NGC 253 using ALMA. We find that this feature is clearly associated with the edge of NGC 253's prominent ionized outflow, has a projected length of ~300 pc, with a width of ~50 pc and a velocity dispersion of ~40 km s^-1, consistent with an ejection from the disk about 1 Myr ago. The kinematics of the molecular gas in this feature can be interpreted (albeit not uniquely) as accelerating at a rate of 1 km s^-1 pc^-1. In this scenario, the gas is approaching escape velocity at the last measured point. Strikingly, bright tracers of dense molecular gas (HCN, CN, HCO+, CS) are also detected in the molecular outflow: We measure an HCN(1-0)/CO(1-0) line ratio of ~1/10 in the outflow, similar to that in the central starburst region of NGC 253 and other starburst galaxies. By contrast, the HCN/CO line ratio in the NGC 253 disk is significantly lower (~1/30), similar to other nearby galaxy disks. This strongly suggests that the streamer gas originates from the starburst, and that its physical state does not change significantly over timescales of ~1 Myr during its entrainment in the outflow. Simple calculations indicate that radiation pressure is not the main mechanism for driving the outflow. The presence of such dense material in molecular outflows needs to be accounted for in simulations of galactic outflows.

We present high-resolution (0.16$"$) 870um Atacama Large Millimeter/submillimeter Array (ALMA) imaging of 16 luminous (L_IR ~ 4 x 10^12 L_sun) submillimeter galaxies (SMGs) from the ALESS survey of the Extended Chandra Deep Field South. This dust imaging traces the dust-obscured star formation in these z~2.5 galaxies on ~1.3 kpc scales. The emission has a median effective radius of $R_e=0.24" \pm 0.02"$, corresponding to a typical physical size of $R_{e}=1.8\pm$0.2 kpc. We derive a median S\'ersic index of $n=0.9\pm0.2$, implying that the dust emission is remarkably disk-like at the current resolution and sensitivity. We use different weighting schemes with the visibilities to search for clumps on 0.12$"$ (~1.0 kpc) scales, but we find no significant evidence for clumping in the majority of cases. Indeed, we demonstrate using simulations that the observed morphologies are generally consistent with smooth exponential disks, suggesting that caution should be exercised when identifying candidate clumps in even moderate S/N interferometric data. We compare our maps to comparable-resolution HST H$_{160}$-band images, finding that the stellar morphologies appear significantly more extended and disturbed, and suggesting that major mergers may be responsible for driving the formation of the compact dust disks we observe. The stark contrast between the obscured and unobscured morphologies may also have implications for SED fitting routines that assume the dust is co-located with the optical/near-IR continuum emission. Finally, we discuss the potential of the current bursts of star formation to transform the observed galaxy sizes and light profiles, showing that the $z\sim0$ descendants of these SMGs are expected to have stellar masses, effective radii, and gas surface densities consistent with the most compact massive (M_* ~ 1-2 x 10^11 M_sun) early-type galaxies observed locally.

We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

We present an analysis of a deep (1$\sigma$=13 $\mu$Jy) cosmological 1.2-mm continuum map based on ASPECS, the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field. In the 1 arcmin$^2$ covered by ASPECS we detect nine sources at $>3.5\sigma$ significance at 1.2-mm. Our ALMA--selected sample has a median redshift of $z=1.6\pm0.4$, with only one galaxy detected at z$>$2 within the survey area. This value is significantly lower than that found in millimeter samples selected at a higher flux density cut-off and similar frequencies. Most galaxies have specific star formation rates similar to that of main sequence galaxies at the same epoch, and we find median values of stellar mass and star formation rates of $4.0\times10^{10}\ M_\odot$ and $\sim40~M_\odot$ yr$^{-1}$, respectively. Using the dust emission as a tracer for the ISM mass, we derive depletion times that are typically longer than 300 Myr, and we find molecular gas fractions ranging from $\sim$0.1 to 1.0. As noted by previous studies, these values are lower than using CO--based ISM estimates by a factor $\sim$2. The 1\,mm number counts (corrected for fidelity and completeness) are in agreement with previous studies that were typically restricted to brighter sources. With our individual detections only, we recover $55\pm4\%$ of the extragalactic background light (EBL) at 1.2 mm measured by the Planck satellite, and we recover $80\pm7\%$ of this EBL if we include the bright end of the number counts and additional detections from stacking. The stacked contribution is dominated by galaxies at $z\sim1-2$, with stellar masses of (1-3)$\times$10$^{10}$ M$_\odot$. For the first time, we are able to characterize the population of galaxies that dominate the EBL at 1.2 mm.

We measure the sizes of redshift ~2 star-forming galaxies by stacking data from the Atacama Large Millimeter/submillimeter Array (ALMA). We use a uv-stacking algorithm in combination with model fitting in the uv-domain and show that this allows for robust measures of the sizes of marginally resolved sources. The analysis is primarily based on the 344 GHz ALMA continuum observations centred on 88 sub-millimeter galaxies in the LABOCA ECDFS Submillimeter Survey (ALESS). We study several samples of galaxies at z~2 with $M_* \sim{} 5\times{}10^{10} $ M$_\odot$ , selected using near-infrared photometry (distant red galaxies, extremely red objects, sBzK-galaxies, and galaxies selected on photometric redshift). We find that the typical sizes of these galaxies are ~0.6 arcsec which corresponds to ~5 kpc at z~2, this agrees well with the median sizes measured in the near-infrared z-band (~0.6 arcsec). We find errors on our size estimates of ~0.1-0.2 arcsec, which agree well with the expected errors for model fitting at the given signal-to-noise ratio. With the uv-coverage of our observations (18-160 m), the size and flux density measurements are sensitive to scales out to 2 . We compare this to a simulated ALMA Cycle 3 dataset with intermediate length baseline coverage, and we find that, using only these baselines, the measured stacked flux density would be an order of magnitude fainter. This highlights the importance of short baselines to recover the full flux density of high-redshift galaxies.

We present a search for [CII] line and dust continuum emission from optical dropout galaxies at $z>6$ using ASPECS, our ALMA Spectroscopic Survey in the Hubble Ultra-Deep Field (UDF). Our observations, which cover the frequency range $212-272$ GHz, encompass approximately the range $6<z<8$ for [CII] line emission and reach a limiting luminosity of L$_{\rm [CII]}\sim$(1.6-2.5)$\times$10$^{8}$ L$_{\odot}$. We identify fourteen [CII] line emitting candidates in this redshift range with significances $>$4.5 $\sigma$, two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60\% of our candidates are expected to be spurious. For one of our blindly selected [CII] line candidates, we tentatively detect the CO(6-5) line in our parallel 3-mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all [CII] candidates results in a tentative detection with $S_{1.2mm}=14\pm5\mu$Jy. This implies a dust-obscured star formation rate (SFR) of $(3\pm1)$ M$_\odot$ yr$^{-1}$. We find that the two highest--SFR objects have candidate [CII] lines with luminosities that are consistent with the low-redshift $L_{\rm [CII]}$ vs. SFR relation. The other candidates have significantly higher [CII] luminosities than expected from their UV--based SFR. At the current sensitivity it is unclear whether the majority of these sources are intrinsically bright [CII] emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for [CII] emitters at $6<z<8$ that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.

We present direct estimates of the mean sky brightness temperature in observing bands around 99GHz and 242GHz due to line emission from distant galaxies. These values are calculated from the summed line emission observed in a blind, deep survey for specrtal line emission from high redshift galaxies using ALMA (the 'ASPECS' survey). In the 99 GHz band, the mean brightness will be dominated by rotational transitions of CO from intermediate and high redshift galaxies. In the 242GHz band, the emission could be a combination of higher order CO lines, and possibly [CII] 158$\mu$m line emission from very high redshift galaxies ($z \sim 6$ to 7). The mean line surface brightness is a quantity that is relevant to measurements of spectral distortions of the cosmic microwave background, and as a potential tool for studying large-scale structures in the early Universe using intensity mapping. While the cosmic volume and the number of detections are admittedly small, this pilot survey provides a direct measure of the mean line surface brightness, independent of conversion factors, excitation, or other galaxy formation model assumptions. The mean surface brightness in the 99GHZ band is: $T_B = 0.94\pm 0.09$ $\mu$K. In the 242GHz band, the mean brightness is: $T_B = 0.55\pm 0.033$ $\mu$K. These should be interpreted as lower limits on the average sky signal, since we only include lines detected individually in the blind survey, while in a low resolution intensity mapping experiment, there will also be the summed contribution from lower luminosity galaxies that cannot be detected individually in the current blind survey.

We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%.

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $\mu$Jy beam$^{-1}$ and 12.7 $\mu$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

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.

We summarize the design, capabilities, and some of the priority science goals of a next generation Very Large Array (ngVLA). The ngVLA is an interferometric array with 10x larger effective collecting area and 10x higher spatial resolution than the current VLA and the Atacama Large Millimeter Array (ALMA), optimized for operation in the wavelength range 0.3cm to 3cm. The ngVLA opens a new window on the Universe through ultra-sensitive imaging of thermal line and continuum emission down to milliarcecond resolution, as well as unprecedented broad band continuum polarimetric imaging of non-thermal processes. The continuum resolution will reach 9mas at 1cm, with a brightness temperature sensitivity of 6K in 1 hour. For spectral lines, the array at 1" resolution will reach 0.3K surface brightness sensitivity at 1cm and 10 km/s spectral resolution in 1 hour. These capabilities are the only means with which to answer a broad range of critical scientific questions in modern astronomy, including direct imaging of planet formation in the terrestrial-zone, studies of dust-obscured star formation and the cosmic baryon cycle down to pc-scales out to the Virgo cluster, making a cosmic census of the molecular gas which fuels star formation back to first light and cosmic reionization, and novel techniques for exploring temporal phenomena from milliseconds to years. The ngVLA is optimized for observations at wavelengths between the superb performance of ALMA at submm wavelengths, and the future SKA1 at few centimeter and longer wavelengths. This memo introduces the project. The science capabilities are outlined in a parallel series of white papers. We emphasize that this initial set of science goals are simply a starting point for the project. We invite comment on these programs, as well as new ideas, through our public forum link on the ngVLA web page https://science.nrao.edu/futures/ngvla

The Next-Generation Very Large Array (ngVLA) will be critical for understanding how galaxies are built and evolve at the earliest epochs. The sensitivity and frequency coverage will allow for the detection of cold gas and dust in `normal' distant galaxies, including the low-J transitions of molecular gas tracers such as CO, HNC, and HCO+; synchrotron and free-free continuum emission; and even the exciting possibility of thermal dust emission at the highest (z~7) redshifts. In particular, by enabling the total molecular gas reservoirs to be traced to unprecedented sensitivities across a huge range of epochs simultaneously -- something no other radio or submillimeter facility will be capable of -- the detection of the crucial low-J transitions of CO in a diverse body of galaxies will be the cornerstone of ngVLA's contribution to high-redshift galaxy evolution science. The ultra-wide bandwidths will allow a complete sampling of radio SEDs, as well as the detection of emission lines necessary for spectroscopic confirmation of elusive dusty starbursts. The ngVLA will also deliver unique contributions to our understanding of cosmic magnetism and to science accessible through microwave polarimetry. Finally, the superb angular resolution will move the field beyond detection experiments and allow detailed studies of the morphology and dynamics of these systems, including dynamical modeling of disks/mergers, determining the properties of outflows, measuring black hole masses from gas disks, and resolving multiple AGN nuclei. We explore the contribution of a ngVLA to these areas and more, as well as synergies with current and upcoming facilities including ALMA, SKA, large single-dish submillimeter observatories, GMT/TMT, and JWST.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

[abridged] The ALESS survey has followed-up a sample of 122 sub-millimeter sources in the Extended Chandra Deep Field South at 870um with ALMA, allowing to pinpoint the positions of sub-millimeter galaxies (SMGs) to 0.3'' and to find their precise counterparts at different wavelengths. This enabled the first compilation of the multi-wavelength spectral energy distributions (SEDs) of a statistically reliable survey of SMGs. In this paper, we present a new calibration of the MAGPHYS modelling code that is optimized to fit these UV-to-radio SEDs of z>1 star-forming galaxies using an energy balance technique to connect the emission from stellar populations, dust attenuation and dust emission in a physically consistent way. We derive statistically and physically robust estimates of the photometric redshifts and physical parameters for the ALESS SMGs. We find that they have a median stellar mass $M_\ast=(8.9\pm0.1)\times10^{10} M_\odot$, SFR$=280\pm70 M_\odot$/yr, overall V-band dust attenuation $A_V=1.9\pm0.2$ mag, dust mass $M_\rm{dust}=(5.6\pm1.0)\times10^8 M_\odot$, and average dust temperature Tdust~40 K. The average intrinsic SED of the ALESS SMGs resembles that of local ULIRGs in the IR range, but the stellar emission of our average SMG is brighter and bluer, indicating lower dust attenuation, possibly because they are more extended. We explore how the average SEDs vary with different parameters, and we provide a new set of SMG templates. To put the ALESS SMGs into context, we compare their stellar masses and SFRs with those of less actively star-forming galaxies at the same redshifts. At z~2, about half of the SMGs lie above the star-forming main sequence, while half are at the high-mass end of the sequence. At higher redshifts (z~3.5), the SMGs tend to have higher SFR and Mstar, but the fraction of SMGs that lie significantly above the main sequence decreases to less than a third.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum images of the asteroid 3 Juno obtained with an angular resolution of 0.042 arcseconds (60 km at 1.97 AU). The data were obtained over a single 4.4 hr interval, which covers 60% of the 7.2 hr rotation period, approximately centered on local transit. A sequence of ten consecutive images reveals continuous changes in the asteroid's profile and apparent shape, in good agreement with the sky projection of the three-dimensional model of the Database of Asteroid Models from Inversion Techniques. We measure a geometric mean diameter of 259pm4 km, in good agreement with past estimates from a variety of techniques and wavelengths. Due to the viewing angle and inclination of the rotational pole, the southern hemisphere dominates all of the images. The median peak brightness temperature is 215pm13 K, while the median over the whole surface is 197pm15 K. With the unprecedented resolution of ALMA, we find that the brightness temperature varies across the surface with higher values correlated to the subsolar point and afternoon areas, and lower values beyond the evening terminator. The dominance of the subsolar point is accentuated in the final four images, suggesting a reduction in the thermal inertia of the regolith at the corresponding longitudes, which are possibly correlated to the location of the putative large impact crater. These results demonstrate ALMA's potential to resolve thermal emission from the surface of main belt asteroids, and to measure accurately their position, geometric shape, rotational period, and soil characteristics.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($\alpha$), which ranges from $\alpha\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

Radio-loud Active Galactic Nuclei at z~2-4 are typically located in dense environments and their host galaxies are among the most massive systems at those redshifts, providing key insights for galaxy evolution. Finding radio-loud quasars at the highest accessible redshifts (z~6) is important to study their properties and environments at even earlier cosmic time. They would also serve as background sources for radio surveys intended to study the intergalactic medium beyond the epoch of reionization in HI 21 cm absorption. Currently, only five radio-loud ($R=f_{\nu,5{\rm GHz}}/f_{\nu,4400\AA}>10$) quasars are known at z~6. In this paper we search for 5.5 < z < 7.2 quasars by cross-matching the optical Pan-STARRS1 and radio FIRST surveys. The radio information allows identification of quasars missed by typical color-based selections. While we find no good 6.4 < z <7.2 quasar candidates at the sensitivities of these surveys, we discover two new radio-loud quasars at z~6. Furthermore, we identify two additional z~6 radio-loud quasars which were not previously known to be radio-loud, nearly doubling the current z~6 sample. We show the importance of having infrared photometry for z>5.5 quasars to robustly classify them as radio-quiet or radio-loud. Based on this, we reclassify the quasar J0203+0012 (z=5.72), previously considered radio-loud, to be radio-quiet. Using the available data in the literature, we constrain the radio-loud fraction of quasars at z~6, using the Kaplan--Meier estimator, to be $8.1^{+5.0}_{-3.2}\%$. This result is consistent with there being no evolution of the radio-loud fraction with redshift, in contrast to what has been suggested by some studies at lower redshifts.

We investigate the CO excitation of normal star forming disk galaxies at z=1.5 using IRAM PdBI observations of the CO[2-1], CO[3-2] and CO[5-4] transitions for 4 galaxies, including VLA observations of CO[1-0] for 3 of them, with the aim of constraining the average state of H2 gas. Exploiting prior knowledge of the velocity range, spatial extent and size of the CO emission we measure reliable line fluxes with S/N>4-7 for individual transitions. While the average CO Spectral Line Energy Distribution (SLED) has a sub-thermal excitation similar to the Milky Way (MW) up to CO[3-2], we show that the average CO[5-4] emission is four times stronger than assuming MW excitation. This demonstrates the presence of an additional component of more excited, denser and possibly warmer molecular gas. The ratio of CO[5-4] to lower-J CO emission is lower than in local (U)LIRGs and high-redshift SMGs, correlating closely with the average intensity of the radiation field <U> and with the star formation surface density, but not with the SF efficiency (SFE). The CO[5-4] luminosity correlates linearly with LIR over 4 orders of magnitudes, with z=1.5 BzK galaxies following the same trend as local spirals and (U)LIRGs and high redshift star bursting SMGs. The CO[5-4] luminosity is thus empirically related to the dense gas, and might be a more convenient way to probe it than standard high--density tracers that are much fainter than CO. We see excitation variations among our sample galaxies, that can be linked to their evolutionary state and clumpiness in optical rest frame images. In one galaxy we see spatially resolved excitation variations, where the more highly excited part of the galaxy corresponds to the location of massive SF clumps. This supports to models that suggest that giant clumps are the main source of the high excitation CO emission in high redshift disk-like galaxies.

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.

We present high-resolution observations of the 880 $\mu$m (rest-frame FIR) continuum emission in the z$=$4.05 submillimeter galaxy GN20 from the IRAM Plateau de Bure Interferometer (PdBI). These data resolve the obscured star formation in this unlensed galaxy on scales of 0.3$^{\prime\prime}$$\times$0.2$^{\prime\prime}$ ($\sim$2.1$\times$1.3 kpc). The observations reveal a bright (16$\pm$1 mJy) dusty starburst centered on the cold molecular gas reservoir and showing a bar-like extension along the major axis. The striking anti-correlation with the HST/WFC3 imaging suggests that the copious dust surrounding the starburst heavily obscures the rest-frame UV/optical emission. A comparison with 1.2 mm PdBI continuum data reveals no evidence for variations in the dust properties across the source within the uncertainties, consistent with extended star formation, and the peak star formation rate surface density (119$\pm$8 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$) implies that the star formation in GN20 remains sub-Eddington on scales down to 3 kpc$^2$. We find that the star formation efficiency is highest in the central regions of GN20, leading to a resolved star formation law with a power law slope of $\Sigma_{\rm SFR}$ $\sim$ $\Sigma_{\rm H_2}^{\rm 2.1\pm1.0}$, and that GN20 lies above the sequence of normal star-forming disks, implying that the dispersion in the star formation law is not due solely to morphology or choice of conversion factor. These data extend previous evidence for a fixed star formation efficiency per free-fall time to include the star-forming medium on $\sim$kpc-scales in a galaxy 12 Gyr ago.

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 use ALMA to derive the mass, length, and time scales associated with the nuclear starburst in NGC 253. This region forms ~2 M_sun/yr of stars and resembles other starbursts in scaling relations, with star formation consuming the gas reservoir 10 times faster than in galaxy disks. We present observations of CO, the high effective density transitions HCN(1-0), HCO+(1-0), CS(2-1), and their isotopologues. We identify ten clouds that appear as peaks in line emission and enhancements in the HCN-to-CO ratio. These clouds are massive (~10^7 M_sun) structures with sizes (~30 pc) similar to GMCs in other systems. Compared to disk galaxy GMCs, they show high line widths (~20-40 km/s) given their size, with implied Mach numbers ~90. The clouds also show high surface (~6,000 M_sun/pc^2) and volume densities (n_H2~2,000 cm^-3). Given these, self-gravity can explain the line widths. This short free fall time (~0.7 Myr) helps explain the more efficient star formation in NGC 253. We also consider the starburst region as a whole. The geometry is confused by the high inclination, but simple models support a non-axisymmetric, bar-like geometry with a compact, clumpy region of high gas density embedded in an extended CO distribution. Even for the whole region, the surface density still exceeds that of a disk galaxy GMC. The orbital time (~10 Myr), disk free fall time (<~ 3 Myr), and disk crossing time (<~ 3 Myr) are each much shorter than in a normal spiral galaxy disk. Some but not all aspects of the structure correspond to predictions from assuming vertical dynamical equilibrium or a marginally stable rotating disk. Finally, the CO-to-H2 conversion factor implied by our cloud calculations is approximately Galactic, contrasting with results showing a low value for the whole starburst region. The contrast provides resolved support for the idea of mixed molecular ISM phases in starburst galaxies.

We present ALLSMOG, the APEX Low-redshift Legacy Survey for MOlecular Gas. ALLSMOG is a survey designed to observe the CO(2-1) emission line with the APEX telescope, in a sample of local galaxies (0.01 < z < 0.03), with stellar masses in the range 8.5 < log(M*/Msun) < 10. This paper is a data release and initial analysis of the first two semesters of observations, consisting of 42 galaxies observed in CO(2-1). By combining these new CO(2-1) emission line data with archival HI data and SDSS optical spectroscopy, we compile a sample of low-mass galaxies with well defined molecular gas masses, atomic gas masses, and gas-phase metallicities. We explore scaling relations of gas fraction and gas consumption timescale, and test the extent to which our findings are dependent on a varying CO/H2 conversion factor. We find an increase in the H2/HI mass ratio with stellar mass which closely matches semi-analytic predictions. We find a mean molecular gas fraction for ALLSMOG galaxies of MH2/M* = (0.09 - 0.13), which decreases with stellar mass. We measure a mean molecular gas consumption timescale for ALLSMOG galaxies of 0.4 - 0.7 Gyr. We also confirm the non-universality of the molecular gas consumption timescale, which varies (with stellar mass) from ~100 Myr to ~2 Gyr. Importantly, we find that the trends in the H2/HI mass ratio, gas fraction, and the non-universal molecular gas consumption timescale are all robust to a range of recent metallicity-dependent CO/H2 conversion factors.

We measure the dust and gas content of the three sub-millimeter galaxies (SMGs) in the GN20 proto-cluster at $z=4.05$ using new IRAM Plateau de Bure interferometer (PdBI) CO(4-3) and 1.2-3.3 mm continuum observations. All these three SMGs are heavily dust obscured, with UV-based star formation rate (SFR) estimates significantly smaller than the ones derived from the bolometric infrared (IR), consistent with the spatial offsets revealed by HST and CO imaging. Based also on evaluations of the specific SFR, CO-to-H$_2$ conversion factor and gas depletion timescale, we classify all the three galaxies as starbursts (SBs), although with a lower confidence for GN20.2b that might be a later stage merging event. We place our measurements in the context of the evolutionary properties of main sequence (MS) and SB galaxies. ULIRGs have 3-5 times larger $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios than $z=0$ MS galaxies, but by $z\sim2$ the difference appears to be blurred, probably due to differential metallicity evolution. SB galaxies appear to slowly evolve in their $L'_{\rm CO}/M_{\rm dust}$ and $M_{\rm dust}/M_\star$ ratios all the way to $z>6$ (consistent with rapid enrichment of SB events), while MS galaxies rapidly increase in $M_{\rm dust}/M_\star$ from $z=0$ to 2 (due to gas fraction increase, compensated by a decrease of metallicities). While no IR/submm continuum detection is available for indisputably normal massive galaxies at $z>2.5$, we show that if metallicity indeed decrease rapidly for these systems at $z>3$ as claimed in the literature, we should expect a strong decrease of their $M_{\rm dust}/M_\star$, consistent with recent PdBI and ALMA upper limits. We conclude that the $M_{\rm dust}/M_\star$ ratio could be a powerful tool for distinguishing starbursts from normal galaxies at $z>4$.

We present ALMA observations of the [CII] line and far-infrared (FIR) continuum of a normally star-forming galaxy in the reionization epoch, the z=6.96 Ly-alpha emitter (LAE) IOK-1. Probing to sensitivities of sigma_line = 240 micro-Jy/beam (40 km/s channel) and sigma_cont = 21 micro-Jy/beam, we found the galaxy undetected in both [CII] and continuum. Comparison of UV - FIR spectral energy distribution (SED) of IOK-1, including our ALMA limit, with those of several types of local galaxies (including the effects of the cosmic microwave background, CMB, on the FIR continuum) suggests that IOK-1 is similar to local dwarf/irregular galaxies in SED shape rather than highly dusty/obscured galaxies. Moreover, our 3 sigma FIR continuum limit, corrected for CMB effects, implies intrinsic dust mass M_dust < 6.4 x 10^7 M_sun, FIR luminosity L_FIR < 3.7 x 10^{10} L_sun (42.5 - 122.5 micron), total IR luminosity L_IR < 5.7 x 10^{10} L_sun (8 - 1000 micron) and dust-obscured star formation rate (SFR) < 10 M_sun/yr, if we assume that IOK-1 has a dust temperature and emissivity index typical of local dwarf galaxies. This SFR is 2.4 times lower than one estimated from the UV continuum, suggesting that < 29% of the star formation is obscured by dust. Meanwhile, our 3 sigma [CII] flux limit translates into [CII] luminosity, L_[CII] < 3.4 x 10^7 L_sun. Locations of IOK-1 and previously observed LAEs on the L_[CII] vs. SFR and L_[CII]/L_FIR vs. L_FIR diagrams imply that LAEs in the reionization epoch have significantly lower gas and dust enrichment than AGN-powered systems and starbursts at similar/lower redshifts, as well as local star-forming galaxies.

We present detections of the CO(J= 1-0) emission line in a sample of four massive star-forming galaxies at z~1.5-2.2 obtained with the Karl G. Jansky Very Large Array (VLA). Combining these observations with previous CO(2-1) and CO(3-2) detections of these galaxies, we study the excitation properties of the molecular gas in our sample sources. We find an average line brightness temperature ratios of R_{21}=0.70+\-0.16 and R_{31}=0.50+\-0.29, based on measurements for three and two galaxies, respectively. These results provide additional support to previous indications of sub-thermal gas excitation for the CO(3-2) line with a typically assumed line ratio R_{31}~0.5. For one of our targets, BzK-21000, we present spatially resolved CO line maps. At the resolution of 0.18 arcsec (1.5 kpc), most of the emission is resolved out except for some clumpy structure. From this, we attempt to identify molecular gas clumps in the data cube, finding 4 possible candidates. We estimate that <40 % of the molecular gas is confined to giant clumps (~1.5 kpc in size), and thus most of the gas could be distributed in small fainter clouds or in fairly diffuse extended regions of lower brightness temperatures than our sensitivity limit.

We present and release co-added images of the Sloan Digital Sky Survey (SDSS) Stripe 82. Stripe 82 covers an area of 300 deg^2 on the Celestial Equator, and has been repeatedly scanned 70-90 times in the ugriz bands by the SDSS imaging survey. By making use of all available data in the SDSS archive, our co-added images are optimized for depth. Input single-epoch frames were properly processed and weighted based on seeing, sky transparency, and background noise before co-addition. The resultant products are co-added science images and their associated weight images that record relative weights at individual pixels. The depths of the co-adds, measured as the 5 sigma detection limits of the aperture (3.2 arcsec diameter) magnitudes for point sources, are roughly 23.9, 25.1, 24.6, 24.1, and 22.8 AB magnitudes in the five bands, respectively. They are 1.9-2.2 mag deeper than the best SDSS single-epoch data. The co-added images have good image quality, with an average point-spread function FWHM of ~1 arcsec in the r, i, and z bands. We also release object catalogs that were made with SExtractor. These co-added products have many potential uses for studies of galaxies, quasars, and Galactic structure. We further present and release near-IR J-band images that cover ~90 deg^2 of Stripe 82. These images were obtained using the NEWFIRM camera on the NOAO 4-m Mayall telescope, and have a depth of about 20.0--20.5 Vega magnitudes (also 5 sigma detection limits for point sources).

We present a study of the radio properties of 870$\mu$m-selected submillimetre galaxies (SMGs), observed at high resolution with ALMA in the Extended Chandra Deep Field South. From our initial sample of 76 ALMA SMGs, we detect 52 SMGs at $>3\sigma$ significance in VLA 1400MHz imaging, of which 35 are also detected at $>3\sigma$ in new 610MHz GMRT imaging. Within this sample of radio-detected SMGs, we measure a median radio spectral index $\alpha_{610}^{1400} = -0.79 \pm 0.06$, (with inter-quartile range $\alpha=[-1.16,-0.56]$) and investigate the far-infrared/radio correlation via the parameter $q_{\rm IR}$, the logarithmic ratio of the rest-frame 8-1000$\mu$m flux and monochromatic radio flux. Our median $q_{\rm IR} = 2.56 \pm 0.05$ (inter-quartile range $q_{\rm IR}=[2.42,2.78]$) is higher than that typically seen in single-dish 870$\mu$m-selected sources ($q_{\rm IR} \sim 2.4$), which may reflect the fact that our ALMA-based study is not biased to radio-bright counterparts, as previous samples were. Finally, we search for evidence that $q_{\rm IR}$ and $\alpha$ evolve with age in a co-dependent manner, as predicted by starburst models: the data populate the predicted region of parameter space, with the stellar mass tending to increase along tracks of $q_{\rm IR}$ versus $\alpha$ in the direction expected, providing the first observational evidence in support of these models.

We present the first photometric redshift distribution for a large unbiased sample of 870um selected submillimeter galaxies (SMGs) with robust identifications based on observations with the Atacama Large Millimeter Array (ALMA). In our analysis we consider 96 SMGs in the Extended Chandra Deep Field South, 77 of which have 4-19 band, optical-near-infrared, photometry. We model the Spectral Energy Distributions (SEDs) for these 77 SMGs, deriving a median photometric redshift of z=2.3+/-0.1. The remaining 19 SMGs have insufficient optical or near-infrared photometry to derive photometric redshifts, but a stacking analysis of IRAC and Herschel observations confirms they are not spurious. Assuming these sources have an absolute H-band magnitude distribution comparable to that of a complete sample of z~1-2 SMGs, we demonstrate that the undetected SMGs lie at higher redshifts, raising the median redshift for SMGs to z=2.5+/-0.2. More critically we show that the proportion of galaxies undergoing an SMG phase at z>3 is 35+/-5% of the total population. We derive a median stellar mass for SMGs of Mstar=(8+/-1)x10^10Mo, but caution that there are significant systematic uncertainties in our stellar mass estimate, up to x5 for individual sources. We compare our sample of SMGs to a volume-limited, morphologically classified sample of ellipticals in the local Universe. Assuming the star formation activity in SMGs has a timescale of ~100Myr we show that their descendants at z~0 would have a space density and M_H distribution which are in good agreement with those of local ellipticals. In addition the inferred mass-weighted ages of the local ellipticals broadly agree with the look-back times of the SMG events. Taken together, these results are consistent with a simple model that identifies SMGs as events that form most of the stars seen in the majority of luminous elliptical galaxies at the present day.

We exploit ALMA 870um (345GHz) observations of submillimetre sources in the Extended Chandra Deep Field South to investigate the far-infrared properties of high-redshift submillimetre galaxies (SMGs). Using the precisely located 870um ALMA positions of 99 SMGs, together with 24um and radio imaging of this field, we deblend the Herschel/SPIRE imaging of this region to extract their far-infrared fluxes and colours. The median photometric redshifts for ALMA LESS (ALESS) SMGs which are detected in at least two SPIRE bands increases with wavelength of the peak in their SEDs, with z=2.3+/-0.2, 2.5+/-0.3 and 3.5+/-0.5 for the 250, 350 and 500-um peakers respectively. We find that 34 ALESS SMGs do not have a >3-sigma counterpart at 250, 350 or 500-um. These galaxies have a median photometric redshift of z=3.3+/-0.5, which is higher than the full ALESS SMG sample; z=2.5+/-0.2. Using the photometric redshifts together with the 250-870um photometry, we estimate the far-infrared luminosities and characteristic dust temperature of each SMG. The median infrared luminosity of the S_870um>2mJy SMGs is L_IR=(3.0+/-0.3)x10^{12}Lo(SFR=300+/-30Mo/yr). At a fixed luminosity, the characteristic dust temperature of these high-redshift SMGs is 2-3K lower than comparably luminous galaxies at z=0, reflecting the more extended star formation occurring in these systems. By extrapolating the 870um number counts to S_ 870um=1mJy, we show that the contribution of S_870um>1mJy SMGs to the cosmic star formation budget is 20% of the total over the redshift range z~1-4. We derive a median dust mass for these SMGs of M_d=(3.6+/-0.3)x10^8Mo and by adopting an appropriate gas-to-dust ratio, we estimate an average molecular mass of M_H2=(4.2+/-0.4)x10^{10}Mo. Finally, we use our estimates of the H2 masses to show that SMGs with S_870um>1mJy contain ~10% of the z~2 volume-averaged H2 mass density at this epoch.

We report the detection of 6.2um polycyclic aromatic hydrocarbon (PAH) and rest-frame 4-7um continuum emission in the z=4.055 submillimeter galaxy GN20, using the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope. This represents the first detection of PAH emission at z>4. The strength of the PAH emission feature is consistent with a very high star formation rate of ~1600Msun/yr. We find that this intense starburst powers at least ~1/3 of the faint underlying 6um continuum emission, with an additional, significant (and perhaps dominant) contribution due to a power-law-like hot dust source, which we interpret to likely be a faint, dust-obscured active galactic nucleus (AGN). The inferred 6um AGN continuum luminosity is consistent with a sensitive upper limit on the hard X-ray emission as measured by the Chandra X-Ray Observatory if the previously undetected AGN is Compton-thick. This is in agreement with the finding at optical/infrared wavelengths that the galaxy and its nucleus are heavily dust-obscured. Despite the strong power-law component enhancing the mid-infrared continuum emission, the intense starburst associated with the photon-dominated regions that give rise to the PAH emission appears to dominate the total energy output in the infrared. GN20 is one of the most luminous starburst galaxies known at any redshift, embedded in a rich protocluster of star-forming galaxies. This investigation provides an improved understanding of the energy sources that power such exceptional systems, which represent the extreme end of massive galaxy formation at early cosmic times.

We study the sub-mm properties of color-selected galaxies via a stacking analysis applied for the first time to interferometric data at sub-mm wavelengths. We base our study on 344 GHz ALMA continuum observations of ~20''-wide fields centered on 86 sub-mm sources detected in the LABOCA Extended Chandra Deep Field South Sub-mm Survey (LESS). We select various classes of galaxies (K-selected, star-forming sBzK galaxies, extremely red objects and distant red galaxies) according to their optical/NIR fluxes. We find clear, >10-sigma detections in the stacked images of all these galaxy classes. We include in our stacking analysis Herschel/SPIRE data to constrain the dust SED of these galaxies. We find that their dust emission is well described by a modified black body with Tdust~30 K and beta=1.6 and IR luminosities of (5-11)x10^{11} Lsun, or implied star formation rates of 75-140 Msun/yr. We compare our results with those of previous studies based on single-dish observations at 870 micron and find that our flux densities are a factor 2-3 higher than previous estimates. The discrepancy is observed also after removing sources individually detected in ALESS maps. We report a similar discrepancy by repeating our analysis on 1.4\,GHz observations of the whole ECDFS. Hence we find tentative evidence that galaxies that are associated in projected and redshift space with sub-mm bright sources are brighter than the average population. Finally, we put our findings in the context of the cosmic star formation rate density as a function of redshift.

The large gas and dust reservoirs of submm galaxies (SMGs) could potentially provide ample fuel to trigger an Active Galactic Nucleus (AGN), but previous studies of the AGN fraction in SMGs have been controversial largely due to the inhomogeneity and limited angular resolution of the available submillimeter surveys. Here we set improved constraints on the AGN fraction and X-ray properties of the SMGs with ALMA and Chandra observations in the Extended Chandra Deep Field-South (E-CDF-S). This study is the first among similar works to have unambiguously identified the X-ray counterparts of SMGs; this is accomplished using the fully submm-identified, statistically reliable SMG catalog with 99 SMGs from the ALMA LABOCA E-CDF-S Submillimeter Survey (ALESS). We found 10 X-ray sources associated with SMGs (median redshift z = 2.3), of which 8 were identified as AGNs using several techniques that enable cross-checking. The other 2 X-ray detected SMGs have levels of X-ray emission that can be plausibly explained by their star-formation activity. 6 of the 8 SMG-AGNs are moderately/highly absorbed, with N_H > 10e23 cm-2. An analysis of the AGN fraction, taking into account the spatial variation of X-ray sensitivity, yields an AGN fraction of 17+16-6% for AGNs with rest-frame 0.5-8 keV absorption-corrected luminosity >7.8x10e42 erg s-1; we provide estimated AGN fractions as a function of X-ray flux and luminosity. ALMA's high angular resolution also enables direct X-ray stacking at the precise positions of SMGs for the first time, and we found 4 potential SMG-AGNs in our stacking sample.

We present deep IRAM Plateau de Bure Interferometer (PdBI) observations, searching for CO-emission toward two massive, non-lensed Lyman break galaxies (LBGs) at z=3.216 and 4.058. With one low significance CO detection (3.5 sigma) and one sensitive upper limit, we find that the CO lines are >~ 3-4 times weaker than expected based on the relation between IR and CO luminosities followed by similarly, massive galaxies at z=0-2.5. This is consistent with a scenario in which these galaxies have low metallicity, causing an increased CO-to-H_2 conversion factor, i.e., weaker CO-emission for a given molecular (H_2) mass. The required metallicities at z>3 are lower than predicted by the fundamental metallicity relation (FMR) at these redshifts, consistent with independent evidence. Unless our galaxies are atypical in this respect, detecting molecular gas in normal galaxies at z>3 may thus remain challenging even with ALMA.

The under-abundance of very massive galaxies in the universe is frequently attributed to the effect of galactic winds. Although ionized galactic winds are readily observable most of the expelled mass is likely in cooler atomic and molecular phases. Expanding molecular shells observed in starburst systems such as NGC 253 and M 82 may facilitate the entrainment of molecular gas in the wind. While shell properties are well constrained, determining the amount of outflowing gas emerging from such shells and the connection between this gas and the ionized wind requires spatial resolution <100 pc coupled with sensitivity to a wide range of spatial scales, hitherto not available. Here we report observations of NGC 253, a nearby starburst galaxy (D~3.4 Mpc) known to possess a wind, which trace the cool molecular wind at 50 pc resolution. At this resolution the extraplanar molecular gas closely tracks the H{\alpha} filaments, and it appears connected to molecular expanding shells located in the starburst region. These observations allow us to directly measure the molecular outflow rate to be > 3 Msun/yr and likely ~9 Msun/yr. This implies a ratio of mass-outflow rate to star formation rate of at least {\eta}~1-3, establishing the importance of the starburst-driven wind in limiting the star formation activity and the final stellar content.

We present a study of the formation of clustered, massive galaxies at large look-back times via spectroscopic imaging of CO in the unique GN20 proto-cluster at z = 4.05. Existing observations show that this is a dense concentration of gas-rich, very active star forming galaxies, including multiple bright submillimeter galaxies (SMGs). Using deep, high-resolution VLA CO(2-1) observations, we image the molecular gas with a resolution of ~1 kpc just 1.6 Gyr after the Big Bang. The SMGs GN20.2a and GN20.2b have deconvolved sizes of ~5 kpc X 3 kpc and ~8 kpc X 5 kpc (Gaussian FWHM) in CO(2-1), respectively, and we measure gas surface densities up to ~12,700/1,700X(sin i) (\alpha_CO/0.8) M_sun/pc^2 for GN20.2a/GN20.2b in the highest-resolution maps. Dynamical mass estimates allow us to constrain the CO-to-H_2 conversion factor to \alpha_CO = 1.7+/-0.8 M_sun (K km s^{-1} pc^2)^-1 for GN20.2a and \alpha_CO = 1.1+/-^{1.5}_{1.1} M_sun (K km/s pc^2)^-1 for GN20.2b. We measure significant offsets (0.5"-1") between the CO and optical emission, indicating either dust obscuration on scales of tens of kpc or that the emission originates from distinct galaxies. CO spectral line energy distributions imply physical conditions comparable to other SMGs and reveal further evidence that GN20.2a and GN20.2b are in different merging stages. We carry out a targeted search for CO emission from the 14 known B-band Lyman break galaxies (LBGs) in the field, tentatively detecting CO in a previously-undetected LBG and placing 3sigma upper limits on the CO luminosities of those that may lie within our bandpass. A blind search for emission line sources down to a 5sigma limiting CO luminosity of L'_CO(2-1) = 8 X 10^9 K km/s pc^2 and covering \Delta z = 0.0273 (~20 comoving Mpc) produces no other strong contenders associated with the proto-cluster.

We report on a blind survey for extragalactic radio variability that was carried out by comparing two epochs of data from the FIRST survey with a third epoch from a new 1.4 GHz survey of SDSS Stripe 82. The three epochs are spaced seven years apart and have an overlapping area of 60 deg^2. We uncover 89 variable sources down to the millijansky level, 75 of which are newly-identified, and we find no evidence for transient phenomena. This new sample of variable sources allows us to infer an upper limit to the mean characteristic timescale of AGN radio variability of 14 years. We find that only 1% of extragalactic sources have fractional variability f_var >3, while 44% of Galactic sources vary by this much. The variable sample contains a larger fraction of quasars than a comparable non-variable control sample, though the majority of the variable sources appear to be extended galaxies in the optical. This implies that either quasars are not the dominant contributor to the variability of the sample, or that the deep optical data allow us to detect the host galaxies of some low-z quasars. We use the new, higher resolution data to report on the morphology of the variable sources. Finally, we show that the fraction of sources that are variable remains constant or increases at low flux densities. This may imply that next generation radio surveys with telescopes like the Australian Square Kilometer Array Pathfinder and MeerKAT will see a constant or even increasing fraction of variable sources down into the submillijansky regime.

We present an Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 survey of 126 submillimeter sources from the LABOCA ECDFS Submillimeter Survey (LESS). Our 870 micron survey with ALMA (ALESS) has produced maps ~3X deeper and with a beam area ~200X smaller than the original LESS observations, doubling the current number of interferometrically-observed submillimeter sources. The high resolution of these maps allows us to resolve sources that were previously blended and accurately identify the origin of the submillimeter emission. We discuss the creation of the ALESS submillimeter galaxy (SMG) catalog, including the main sample of 99 SMGs and a supplementary sample of 32 SMGs. We find that at least 35% (possibly up to 50%) of the detected LABOCA sources have been resolved into multiple SMGs, and that the average number of SMGs per LESS source increases with LESS flux density. Using the (now precisely known) SMG positions, we empirically test the theoretical expectation for the uncertainty in the single-dish source positions. We also compare our catalog to the previously predicted radio/mid-infrared counterparts, finding that 45% of the ALESS SMGs were missed by this method. Our ~1.6" resolution allows us to measure a size of ~9 kpc X 5 kpc for the rest-frame ~300 um emission region in one resolved SMG, implying a star formation rate surface density of 80 M_sol yr^-1 kpc^-2, and we constrain the emission regions in the remaining SMGs to be <10 kpc. As the first statistically reliable survey of SMGs, this will provide the basis for an unbiased multiwavelength study of SMG properties.

We present CO(2-1) observations of the submillimetre galaxy ALESS65.1 performed with the Australia Telescope Compact Array at 42.3 GHz. A previous ALMA study of submillimetre galaxies in the Extended Chandra Deep Field South detected [CII] 157.74 micron emission from this galaxy at a redshift of z = 4.44. No CO(2-1) emission was detected but we derive a firm upper limit to the cold gas mass in ALESS65.1 of M_gas < 1.7 x 10^10 M_odot. The estimated gas depletion timescale is <50 Myr, which is similar to other high redshift SMGs, and consistent with z > 4 SMGs being the likely progenitors of massive red-and-dead galaxies at z > 2. The ratio of the [CII], CO and far-infrared luminosities implies a strong far-ultraviolet field of G_0 > 10^3, as seen in Galactic star forming regions or local ULIRGs. The observed L_[CII]/L_FIR = 2.3 x 10^{-3} is high compared to local ULIRGs and, combined with L_[CII]/L_CO > 2700, it is consistent with ALESS65.1 either having an extended (several kpc) [CII] emitting region or lower than solar metallicity.

We report the first counts of faint submillimetre galaxies (SMG) in the 870-um band derived from arcsecond resolution observations with the Atacama Large Millimeter Array (ALMA). We have used ALMA to map a sample of 122 870-um-selected submillimetre sources drawn from the (0.5x0.5)deg^2 LABOCA Extended Chandra Deep Field South Submillimetre Survey (LESS). These ALMA maps have an average depth of sigma(870um)~0.4mJy, some ~3x deeper than the original LABOCA survey and critically the angular resolution is more than an order of magnitude higher, FWHM of ~1.5" compared to ~19" for the LABOCA discovery map. This combination of sensitivity and resolution allows us to precisely pin-point the SMGs contributing to the submillimetre sources from the LABOCA map, free from the effects of confusion. We show that our ALMA-derived SMG counts broadly agree with the submillimetre source counts from previous, lower-resolution single-dish surveys, demonstrating that the bulk of the submillimetre sources are not caused by blending of unresolved SMGs. The difficulty which well-constrained theoretical models have in reproducing the high-surface densities of SMGs, thus remains. However, our observations do show that all of the very brightest sources in the LESS sample, S(870um)>12mJy, comprise emission from multiple, fainter SMGs, each with 870-um fluxes of <9mJy. This implies a natural limit to the star-formation rate in SMGs of <10^3 M_Sun/yr, which in turn suggests that the space densities of z>1 galaxies with gas masses in excess of ~5x10^10 M_Sun is <10^-5 Mpc^-3. We also discuss the influence of this blending on the identification and characterisation of the SMG counterparts to these bright submillimetre sources and suggest that it may be responsible for previous claims that they lie at higher redshifts than fainter SMGs.

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

We present Karl G. Jansky Very Large Array (VLA) observations of the CO(2-1) emission in the z=4.05 submillimeter galaxy (SMG) GN20. These high-resolution data allow us to image the molecular gas at 1.3 kpc resolution just 1.6 Gyr after the Big Bang. The data reveal a clumpy, extended gas reservoir, 14 +/- 4 kpc in diameter, in unprecedented detail. A dynamical analysis shows that the data are consistent with a rotating disk of total dynamical mass 5.4 +/- 2.4 X 10^11 M_sun. We use this dynamical mass estimate to constrain the CO-to-H_2 mass conversion factor (alpha_CO), finding alpha_CO=1.1 +/- 0.6 M_sun (K km s^-1 pc^2)^-1. We identify five distinct molecular gas clumps in the disk of GN20 with masses a few percent of the total gas mass, brightness temperatures of 16-31K, and surface densities of >3,200-4,500 X (alpha_CO/0.8) M_sun pc^-2. Virial mass estimates indicate they could be self-gravitating, and we constrain their CO-to-H_2 mass conversion factor to be <0.2-0.7 M_sun (K km s^-1 pc^2)^-1. A multiwavelength comparison demonstrates that the molecular gas is concentrated in a region of the galaxy that is heavily obscured in the rest-frame UV/optical. We investigate the spatially-resolved gas excitation and find that the CO(6-5)/CO(2-1) ratio is constant with radius, consistent with star formation occuring over a large portion of the disk. We discuss the implications of our results in the context of different fueling scenarios for SMGs.

We present ALMA 870-um (345GHz) observations of two sub-millimetre galaxies (SMGs) drawn from an ALMA study of the 126 sub-millimeter sources from the LABOCA Extended Chandra Deep Field South Survey (LESS). The ALMA data identify the counterparts to these previously unidentified sub-millimeter sources and serendipitously detect bright emission lines in their spectra which we show are most likely to be [C II]157.74um emission yielding redshifts of z=4.42 and z=4.44. This blind detection rate within the 7.5-GHz bandpass of ALMA is consistent with the previously derived photometric redshift distribution of SMGs and suggests a modest, but not dominant (<25%), tail of 870-um selected SMGs at z>4. We find that the ratio of L_CII/L_FIR in these SMGs is much higher than seen for similarly far-infrared-luminous galaxies at z~0, which is attributed to the more extended gas reservoirs in these high-redshift ULIRGs. Indeed, in one system we show that the [C II] emission shows hints of extended emission on >3kpc scales. Finally, we use the volume probed by our ALMA survey to show that the bright end of the [C II] luminosity function evolves strongly between z=0 and z~4.4, reflecting the increased ISM cooling in galaxies as a result of their higher star-formation rates. These observations demonstrate that even with short integrations, ALMA is able to detect the dominant fine structure cooling lines from high-redshift ULIRGs, measure their energetics and trace their evolution with redshift.

We present Herschel-PACS spectroscopy of the [OI]63um far-infrared cooling line from a sample of six unlensed and spectroscopically-confirmed 870um-selected submillimetre (submm) galaxies (SMGs) at 1.1<z<1.6 from the LABOCA Extended Chandra Deep Field South (ECDFS) Submm Survey (LESS). This is the first survey of [OI]63um, one of the main photodissociation region (PDR) cooling lines, in SMGs. New high-resolution ALMA interferometric 870um continuum imaging confirms that these six Herschel-targeted SMG counterparts are bona fide sources of submm emission. We detect [OI]63um in two SMGs with a SNR >3, tentatively detect [OI]63um in one SMG, and constrain the line flux for the non-detections. We also exploit the combination of submm continuum photometry from 250-870um and our new PACS continuum measurements to constrain the far-infrared (FIR) luminosity, L_FIR, in these SMGs to < 30%. We find that SMGs do not show a deficit in their [OI]63um-to-far-infrared continuum luminosity ratios (with ratios ranging from ~0.5-1.5%), similar to what was seen previously for the [CII]158um-to-FIR ratios in SMGs. These observed ratios are about an order of magnitude higher than what is seen typically for local ultra luminous infrared galaxies (ULIRGs), which adds to the growing body of evidence that SMGs are not simply `scaled up' versions of local ULIRGs. Rather, the PDR line-to-L_FIR ratios suggest that the star formation modes of SMGs are likely more akin to that of local normal (lower-luminosity) star-forming galaxies, with the bulk of the star formation occurring in extended regions, galaxy-scale (~kpc) in size. These observations represent the first step towards a census of the major PDR cooling lines in typical SMGs that will be attainable with ALMA, enabling detailed modelling to probe the global properties of the star formation and the evolutionary status of SMGs.

The Hubble Deep Field (HDF) is a region in the sky that provides one of the deepest multi-wavelength views of the distant universe and has led to the detection of thousands of galaxies seen throughout cosmic time. An early map of the HDF at a wavelength of 850 microns that is sensitive to dust emission powered by star formation revealed the brightest source in the field, dubbed HDF850.1. For more than a decade, this source remained elusive and, despite significant efforts, no counterpart at shorter wavelengths, and thus no redshift, size or mass, could be identified. Here we report, using a millimeter wave molecular line scan, an unambiguous redshift determination for HDF850.1 of z=5.183. This places HDF850.1 in a galaxy overdensity at z~5.2 in the HDF, corresponding to a cosmic age of only 1.1 Gyr after the Big Bang. This redshift is significantly higher than earlier estimates and higher than most of the >100 sub-millimeter bright galaxies identified to date. The source has a star formation rate of 850 M_sun/yr and is spatially resolved on scales of 5 kpc, with an implied dynamical mass of ~1.3x10^11 M_sun, a significant fraction of which is present in the form of molecular gas. Despite our accurate redshift and position, a counterpart arising from starlight remains elusive.

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.

We present a study of the formation of clustered, massive galaxies at large look--back times, via high resolution spectroscopic imaging of CO in the unique GN20 proto-cluster. The data reveal evidence for rich structure and gas dynamics in unprecedented detail, allowing us to image the molecular gas with a resolution of only 1 kpc just 1.5 Gyr after the Big Bang. These state-of-the-art data give new insight into the detailed physical processes involved in early massive galaxy formation, and they provide a first glimpse of the morphological studies that will become feasible on a regular basis with ALMA.

We explore the gas-to-dust mass ratio (G/D) and the CO luminosity-to-Mgas conversion factor (a_co) of two well studied galaxies in the GOODS-N field, that are expected to have different star forming modes, the starburst GN20 at z=4.05 and the normal star-forming galaxy BzK-21000 at z=1.52. Detailed sampling is available for their Rayleigh-Jeans emission via ground based mm interferometry (1.1-6.6mm) along with Herschel, PACS and SPIRE data that probe the peak of their infrared emission. Using the physically motivated Draine & Li (2007) models, as well as a modified black body function, we measure the dust mass (Md) of the sources and find 2.0^{+0.7}_{-0.6} x 10^{9} Msun for GN20 and 8.6^{+0.6}_{-0.9} x 10^{8} Msun for BzK-21000. The addition of mm data reduces the uncertainties of the derived Md by a factor of ~2, allowing the use of the local G\D vs metallicity relation to place constraints on the a_co values of the two sources. For GN20 we derive a conversion factor of a_co < 1.0 Msun pc^{-2}(K km s^{-1})^{-1}, consistent with that of local ULIRGs, while for BzK-21000 we find a considerably higher value, a_co ~4.0 Msun pc^{-2}(K km s^{-1})^{-1}, in agreement with an independent kinematic derivation reported previously. The implied star formation efficiency is ~25 Lsun/Msun for BzK-21000, a factor of ~5-10 lower than that of GN20. The findings for these two sources support the existence of different disk-like and starburst star-formation modes in distant galaxies, although a larger sample is required to draw statistically robust results

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 the detection of spatially resolved CO(1-0) emission in the z~3.4 submillimeter galaxies (SMGs) SMM J09431+4700 and SMM J13120+4242, using the Expanded Very Large Array (EVLA). SMM J09431+4700 is resolved into the two previously reported millimeter sources H6 and H7, separated by ~30kpc in projection. We derive CO(1-0) line luminosities of L'(CO 1-0) = (2.49+/-0.86) and (5.82+/-1.22) x 10^10 K km/s pc^2 for H6 and H7, and L'(CO 1-0) = (23.4+/-4.1) x 10^10 K km/s pc^2 for SMM J13120+4242. These are ~1.5-4.5x higher than what is expected from simple excitation modeling of higher-J CO lines, suggesting the presence of copious amounts of low-excitation gas. This is supported by the finding that the CO(1-0) line in SMM J13120+4242, the system with lowest CO excitation, appears to have a broader profile and more extended spatial structure than seen in higher-J CO lines (which is less prominently seen in SMM J09431+4700). Based on L'(CO 1-0) and excitation modeling, we find M_gas = 2.0-4.3 and 4.7-12.7 x 10^10 Msun for H6 and H7, and M_gas = 18.7-69.4 x 10^10 Msun for SMM J13120+4242. The observed CO(1-0) properties are consistent with the picture that SMM J09431+4700 represents an early-stage, gas-rich major merger, and that SMM J13120+4242 represents such a system in an advanced stage. This study thus highlights the importance of spatially and dynamically resolved CO(1-0) observations of SMGs to further understand the gas physics that drive star formation in these distant galaxies, which becomes possible only now that the EVLA rises to its full capabilities.

We present observations of the molecular gas in the GN20 proto-cluster of galaxies at $z =4.05$ using the Expanded Very Large Array (EVLA). This group of galaxies is the ideal laboratory for studying the formation of massive galaxies via luminous, gas-rich starbursts within 1.6 Gyr of the Big Bang. We detect three galaxies in the proto-cluster in CO 2-1 emission, with gas masses (H$_2$) between $10^{10}$ and $10^{11} \times (\alpha/0.8)$ M$_\odot$. The emission from the brightest source, GN20, is resolved with a size $\sim 2"$, and has a clear north-south velocity gradient, possibly indicating ordered rotation. The gas mass in GN20 is comparable to the stellar mass ($1.3\times 10^{11} \times (\alpha/0.8)$ M$_\odot$ and $2.3\times 10^{11}$ M$_\odot$, respectively), and the sum of gas plus stellar mass is comparable to the dynamical mass of the system ($\sim 3.4\times 10^{11} [sin(i)/sin(45^o)]^{-2}$ M$_\odot$), within a 5kpc radius. There is also evidence for a tidal tail extending another $2"$ north of the galaxy with a narrow velocity dispersion. GN20 may be a massive, gas rich disk that is gravitationally disturbed, but not completely disrupted. There is one Lyman-break galaxy (BD29079) in the GN20 proto-cluster with an optical spectroscopic redshift within our search volume, and we set a 3$\sigma$ limit to the molecular gas mass of this galaxy of $1.1\times 10^{10} \times (\alpha/0.8)$ M$_\odot$.

We present the discovery of only the second radio-selected, z ~ 6 quasar. We identified SDSS J222843.54+011032.2 (z=5.95) by matching the optical detections of the deep Sloan Digital Sky Survey (SDSS) Stripe 82 with their radio counterparts in the Stripe82 VLA Survey. We also matched the Canadian-France-Hawaiian Telescope Legacy Survey Wide (CFHTLS Wide) with the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey but have yet to find any z ~ 6 quasars in this survey area. The discovered quasar is optically-faint, z = 22.3 and M_{1450} ~ -24.5, but radio-bright, with a flux density of f$_{1.4GHz, peak}$ = 0.31mJy and a radio-loudness of R ~ 1100 (where R = f_{5GHz}/f_{2500}). The i-z color of the discovered quasar places it outside the color selection criteria for existing optical surveys. We conclude by discussing the need for deeper wide-area radio surveys in the context of high-redshift quasars.

We report the detection of spatially extended CO 1-0 and 5-4 emission in the z=2.49 submillimeter galaxy (SMG) J123707+6214, using the Expanded Very Large Array and the Plateau de Bure Interferometer. The large molecular gas reservoir is spatially resolved into two CO(1-0) components (north-east and south-west; previously identified in CO 3-2 emission) with gas masses of 4.3 and 3.5 x 10^10 (alpha_CO/0.8) Msun. We thus find that the optically invisible north-east component slightly dominates the gas mass in this system. The total molecular gas mass derived from the CO(1-0) observations is ~2.5 times larger than estimated from CO(3-2). The two components are at approximately the same redshift, but separated by ~20 kpc in projection. The morphology is consistent with that of an early-stage merger. The total amount of molecular gas is sufficient to maintain the intense 500 Msun/yr starburst in this system for at least ~160 Myr. We derive line brightness temperature ratios of r_31=0.39+/-0.09 and 0.37+/-0.10, and r_51=0.26+/-0.07 and 0.25+/-0.08 in the two components, respectively, suggesting that the J>=3 lines are substantially subthermally excited. This also suggests comparable conditions for star formation in both components. Given the similar gas masses of both components, this is consistent with the comparable starburst strengths observed in the radio continuum emission. Our findings are consistent with other recent studies that find evidence for lower CO excitation in SMGs than in high-z quasar host galaxies with comparable gas masses. This may provide supporting evidence that both populations correspond to different evolutionary stages in the formation of massive galaxies.

We present a high-resolution radio survey of the Sloan Digital Sky Survey (SDSS) Southern Equatorial Stripe, a.k.a. Stripe 82. This 1.4 GHz survey was conducted with the Very Large Array (VLA) primarily in the A-configuration, with supplemental B-configuration data to increase sensitivity to extended structure. The survey has an angular resolution of 1.8" and achieves a median rms noise of 52 microJy/bm over 92 deg^2. This is the deepest 1.4 GHz survey to achieve this large of an area, filling a gap in the phase space between small, deep and large, shallow surveys. It also serves as a pilot project for a larger high-resolution survey with the Expanded Very Large Array (EVLA). We discuss the technical design of the survey and details of the observations, and we outline our method for data reduction. We present a catalog of 17,969 isolated radio components, for an overall source density of ~195 sources/deg^2. The astrometric accuracy of the data is excellent, with an internal check utilizing multiply-observed sources yielding an rms scatter of 0.19" in both right ascension and declination. A comparison to the SDSS DR7 Quasar Catalog further confirms that the astrometry is well tied to the optical reference frame, with mean offsets of 0.02" +/- 0.01" in right ascension, and 0.01" +/- 0.02" in declination. A check of our photometry reveals a small, negative CLEAN-like bias on the level of 35 microJy. We report on the catalog completeness, finding that 97% of FIRST-detected quasars are recovered in the new Stripe 82 radio catalog, while faint, extended sources are more likely to be resolved out by the resolution bias. We conclude with a discussion of the optical counterparts to the catalog sources, including 76 newly-detected radio quasars. The full catalog as well as a search page and cutout server are available online at http://third.ucllnl.org/cgi-bin/stripe82cutout.

Using a high resolution radio image, we successfully resolve the two fold image components B and C of the quasar lens system SDSS J1029+2623. The flux anomalies associated with these two components in the optical regime persist, albeit less strongly, in our radio observations, suggesting that the cluster must be modeled by something more than a single central potential. We argue that placing substructure close to one of the components can account for a flux anomaly with negligible changes in the component positions. Our best fit model has a substructure mass of ~10^8 solar masses up to the mass-sheet degeneracy, located roughly 0.1 arcsecs West and 0.1 arcsecs North of component B. We demonstrate that a positional offset between the centers of the source components can explain the differences between the optical and radio flux ratios.

We detect and study the properties of faint radio AGN in Luminous Red Galaxies (LRGs). The LRG sample comprises 760,000 objects from a catalog of LRG photometric redshifts constructed from the Sloan Digital Sky Survey (SDSS) imaging data, and 65,000 LRGs from the SDSS spectroscopic sample. These galaxies have typical 1.4 GHz flux densities in the 10s-100s of microJy, with the contribution from a low-luminosity AGN dominating any contribution from star formation. To probe the radio properties of such faint objects, we employ a stacking technique whereby FIRST survey image cutouts at each optical LRG position are sorted by the parameter of interest and median-combined within bins. We find that median radio luminosity scales with optical luminosity (L_opt) as L_1.4 GHz ~ L_opt^(beta), where beta appears to decrease from beta ~ 1 at z = 0.4 to beta ~ 0 at z = 0.7, a result which could be indicative of AGN cosmic downsizing. We also find that the overall LRG population, which is dominated by low-luminosity AGN, experiences significant cosmic evolution between z = 0.2 and z = 0.7. This implies a considerable increase in total AGN heating for these massive ellipticals with redshift. By matching against the FIRST catalog, we investigate the incidence and properties of LRGs associated with double-lobed (FR I/II) radio galaxies. (Abridged)

We present results on a survey to find extremely dust-reddened Type-1 Quasars. Combining the FIRST radio survey, the 2MASS Infrared Survey and the Sloan Digital Sky Survey, we have selected a candidate list of 122 potential red quasars. With more than 80% spectroscopically identified objects, well over 50% are classified as dust-reddened Type 1 quasars, whose reddenings (E(B-V)) range from approximately 0.1 to 1.5 magnitudes. They lie well off the color selection windows usually used to detect quasars and many fall within the stellar locus, which would have made it impossible to find these objects with traditional color selection techniques. The reddenings found are much more consistent with obscuration happening in the host galaxy rather than stemming from the dust torus. We find an unusually high fraction of Broad Absorption Line (BAL) quasars at high redshift, all but one of them belonging to the Low Ionization BAL (LoBAL) class and many also showing absorption the metastable FeII line (FeLoBAL). The discovery of further examples of dust-reddened LoBAL quasars provides more support for the hypothesis that BAL quasars (at least LoBAL quasars) represent an early stage in the lifetime of the quasar. The fact that we see such a high fraction of BALs could indicate that the quasar is in a young phase in which quasar feedback from the BAL winds is suppressing star formation in the host galaxy.

We investigate the radio emission of ~185,000 quiescent (optically unclassifiable) galaxies selected from the Sloan Digital Sky Survey (SDSS). By median-stacking FIRST cutouts centered on the optically-selected sources, we are able to reach flux densities down to the 10s of microJy. The quiescent galaxy sample is composed of two subgroups inhabiting vastly different regimes: those targeted for the SDSS MAIN Galaxy Sample (~55%), and those targeted for the Luminous Red Galaxy (LRG) sample (~45%). To investigate the star-formation rates (SFRs) of these quiescent galaxies, we calibrate a radio-SFR conversion using a third sample of star-forming galaxies. Comparing this SFR-indicator with indicators in the optical and UV, we derive conflicting SFR estimates for the MAIN sample quiescent galaxies. These radio-derived SFRs intersect those calculated using the 4000-Angstrom break (D4000) around an SFR of 1 Msun/yr and agree to within a factor of 3 over the range of SFRs. However, we find that the radio-derived SFRs are too high relative to the SFRs estimated for similar populations of galaxies using analysis of UV emission, implying either contamination of the radio by Active Galactic Nuclei (AGN) or incomplete dust modeling. If AGN activity is dominant in these galaxies, then a relation between AGN radio luminosity and galaxy mass is required to explain the observed trends. For the LRGs, on the other hand, we find the radio luminosity to be independent of SFR as derived from D4000, indicating an AGN component dominates their radio emission. AGN-based radio emission often implies the existence of radio jets, providing evidence of a mechanism for low-level feedback in these quiescent LRGs. (Abridged)

We investigate faint radio emission from low- to high-luminosity Active Galactic Nuclei (AGN) selected from the Sloan Digital Sky Survey (SDSS). Their radio properties are inferred by co-adding large ensembles of radio image cut-outs from the FIRST survey, as almost all of the sources are individually undetected. We correlate the median radio flux densities against a range of other sample properties, including median values for redshift, [OIII] luminosity, emission line ratios, and the strength of the 4000A break. We detect a strong trend for sources that are actively undergoing star-formation to have excess radio emission beyond the ~10^28 ergs/s/Hz level found for sources without any discernible star-formation. Furthermore, this additional radio emission correlates well with the strength of the 4000A break in the optical spectrum, and may be used to assess the age of the star-forming component. We examine two subsamples, one containing the systems with emission line ratios most like star-forming systems, and one with the sources that have characteristic AGN ratios. This division also separates the mechanism responsible for the radio emission (star-formation vs. AGN). For both cases we find a strong, almost identical, correlation between [OIII] and radio luminosity, with the AGN sample extending toward lower, and the star-formation sample toward higher luminosities. A clearer separation between the two subsamples is seen as function of the central velocity dispersion of the host galaxy. For systems with similar redshifts and velocity dispersions, the star-formation subsample is brighter than the AGN in the radio by an order of magnitude. This underlines the notion that the radio emission in star-forming systems can dominate the emission associated with the AGN.

The unexplained sunward acceleration $a_\mathrm{P}$ of the Pioneer 10 (P10) and the Pioneer 11 (P11) spacecraft remains a mystery. A scalar potential model (SPM) that derived from considerations of galaxy clusters, of redshift, and of H{\scriptsize{I}} rotation curves of spiral galaxies is applied to the Pioneer Anomaly. Matter is posited to warp the scalar potential $\rho$ field. The changing $\rho$ field along the light path causes the Pioneer Anomaly. The SPM is consistent with the general value of $a_\mathrm{P}$; with the annual periodicity; with the differing $a_\mathrm{P}$ between the spacecraft; with the discrepancy between {\textit{Sigma}} and CHASMP programs at P10 (I) and their closer agreement at P10 (III); with the slowly declining $a_\mathrm{P}$; with the low value of $a_\mathrm{P}$ immediately before the P11's Saturn encounter; with the high uncertainty in the value of $a_\mathrm{P}$ obtained during and after the P11's Saturn encounter; and with the cosmological connection suggested by $a_\mathrm{P} \approx cH_\mathrm{o}$. The effect of the $\rho$ field warp appears as the curvature of space proposed by general relativity (GR). The Hubble Law and $a_\mathrm{P} \approx cH_\mathrm{o}$ are manifestations of the Newtonian spherical property. Therefore, gravitational attraction, the equivalence principle, and the planet ephemeris remain as described by GR. GR corresponds to the SPM in the limit in which the Sources and Sinks may be replaced by a flat and static $\rho$ field such as between cluster cells and on the Solar System scale at a relatively large distance from a Source or Sink.

A scalar potential model (SPM) was developed from considerations of galaxy clusters and of redshift. The SPM is applied to HI rotation curves (RCs) and RC asymmetry of spiral galaxies. The resulting model adds the force of a scalar potential of the host galaxy and of neighboring galaxies to the Newtonian rotation velocity equation. The RC is partitioned radially into regions. The form of the equation for each parameter of each region is the same with differing proportionality constants. Integer values of each equation are determined empirically for each galaxy. Among the sample galaxies, the global properties of galaxies of B band luminosity, of position, and of orientation determine the RC and RC asymmetry. The Source of the scalar field acts as a monopole at distances of a few kpc from the centre of spiral galaxies. The scalar potential field causes Newtonian mechanics to considerably underestimate the mass in galaxies, which is the ``missing mass problem''. The SPM is consistent with RC and RC asymmetry observations of the sample spiral galaxies.

The galaxy central mass $M_\mathrm{c}$ and central velocity dispersion $\sigma_\mathrm{c}$ have been found to correlate with large scale galaxy parameters for samples of galaxies with a limited range of characteristics. A scalar potential model (SPM) that derived from considerations of galaxy clusters, of redshift, of discrete redshift, of H{\scriptsize{I}} rotation curves (RCs) of spiral galaxies and of RC asymmetry is applied to central region parameters. The $\sigma_\mathrm{c}$ and $ M_\mathrm{c}$ are found to correlate to the host galaxy's and neighboring galaxy's B band luminosity. The sample included galaxies with rising, flat and declining RCs; galaxies with a wide range of characteristics; and galaxies excluded from samples of other studies of $\sigma_\mathrm{c}$ relationships. The equations have the same form as the SPM equations for the parameters of the H{\scriptsize{I}} RCs. Because the SPM is consistent with $M_\mathrm{c}$ and $\sigma_\mathrm{c}$ observations of the sample galaxies, the Sources and Sinks act as monopoles at the center of the galaxies around them. This suggests the outward scalar potential force of a Source holds the $M_\mathrm{c}$ from collapse into a supermassive black hole.

We model the complex distribution of atomic hydrogen (HI) in the interstellar medium (ISM) assuming that it is driven entirely by supernovae (SN). We develop and assess two different models. In the first approach, the simulated volume is randomly populated with non-overlapping voids of a range of sizes. This may relate to a snapshot distribution of supernova-remnant voids, although somewhat artificially constrained by the non-overlap criterion. In the second approach, a simplified time evolution (considering momentum conservation as the only governing constraint during interactions) is followed as SN populate the space with the associated input mass and energy. We describe these simulations and present our results in the form of images of the mass and velocity distributions and the associated power spectra. The latter are compared with trends indicated by available observations. In both approaches, we find remarkable correspondence with the observed statistical description of well-studied components of the ISM, wherein the spatial spectra have been found to show significant deviations from the Kolmogorov spectrum. One of the key indications from this study, regardless of whether or not the SN-induced turbulence is the dominant process in the ISM, is that the apparent non-Kolmogorov spectral characteristics (of HI and/or electron column density across thick or thin screens) needed to explain related observations may not at all be in conflict with the underlying turbulence (i.e. the velocity structure) being of Kolmogorov nature. We briefly discuss the limitations of our simulations and the various implications of our results.

A derivation of a theoretical, time average, cosmic microwave background (CMB), Planckian temperature V of the universe remains a challenge. A scalar potential model (SPM) that resulted from considerations of galaxy cells is applied to deriving a value for V. The heat equation is solved for a cell with the boundary conditions of SPM Source and Sink characteristics, with simplified cell characteristics, and with zero initial temperature. The universe is a collection of cells. The CMB radiation is black body radiation with the cells acting as radiators and absorbers. Conventional thermodynamics is applied to calculate V = 2.718 K. The temperature and matter content of cells are finely controlled by a feedback mechanism. Because time is required for matter to flow from Sources to Sinks, the radiation temperature of cells cycles about V after an initial growth phase. If the universe is like an ideal gas in free expansion and is not in thermal equilibrium, then the pressure and volume follow the measured CMB temperature vm = 2.725 K. Therefore, increasing vm >V equates to an expansion pressure on matter and expanding volume.

On the galactic scale the universe is inhomogeneous and redshift $z$ is occasionally less than zero. A scalar potential model (SPM) that links the galaxy scale $z$ to the cosmological scale $z$ of the Hubble Law is postulated. Several differences among galaxy types suggest that spiral galaxies are Sources and that early type, lenticular, and irregular galaxies are Sinks of a scalar potential field. The morphology-radius and the intragalactic medium cluster observations support the movement of matter from Source galaxies to Sink galaxies. A cell structure of galaxy groups and clusters is proposed to resolve a paradox concerning the scalar potential like the Olber's paradox concerning light. For the sample galaxies, the ratio of the luminosity of Source galaxies to the luminosity of Sink galaxies approaches $2.7 \pm 0.1$. An equation is derived from sample data, which is anisotropic and inhomogeneous, relating $z$ of and the distance $D$ to galaxies. The calculated $z$ has a correlation coefficient of 0.88 with the measured $z$ for a sample of 32 spiral galaxies with $D$ calculated using Cepheid variable stars. The equation is consistent with $z<0$ observations of close galaxies. At low cosmological distances, the equation reduces to $z \approx \exp(KD) \, -1 \approx KD$, where $K$ is a constant, positive value. The equation predicts $z$ from galaxies over 18 Gpc distant approaches a constant value on the order of 500. The SPM of $z$ provides a physical basis for the $z$ of particle photons. Further, the SPM qualitatively suggests the discrete variations in $z$, which was reported by W. G. Tifft, 1997, Astrophy. J. 485, 465 and confirmed by others, are consistent with the SPM.

The changing universe model (CUM) describes galaxy parameter relationships (SESAPS '03, session EB 2). The CUM must be successfully applied to cosmological scale observations to be considered a cosmological model. A major component of current cosmological models is the Hubble constant $H_\mathrm{o}$. An equation is derived using the CUM model relating redshift $z$ and the distance $D$ to galaxies and is applied to a sample of 32 spiral galaxies with $D$ calculated using Cepheid variable stars. The equation predicts a galaxy may have $z<0$ in special circumstances. Three elliptical galaxies with peculiar characteristics are discovered to be CUM Sinks. The Sinks give a physical explanation of the ``Virgocentric infall'' and ``Great Attractor'' observations without a large, unobserved mass. At low cosmological distances, the equation reduces to $z \approx \exp(KD) -1 \approx KD$, where $K$ is a constant, positive value. The equation predicts $z$ from galaxies over 23 Gpc distant approaches a constant value on the order of 1000. The CUM gives a physical basis for the Doppler shift of particle photons.

Although the current galaxy models yield calculations consistent with much of the data, many irregularities exist, exceptions have been found to the current models, the $\Lambda$CDM model apparently fails on galaxy scales, dark matter remains elusive, the phenomena at the center of galaxies are only beginning to be addressed, and many observations and empirical relationships are unexplained. A changing universe model (CUM) is proposed that posits the stuff of our universe is continually erupting into our universe from sources at the center of galaxies. In this first test of the CUM, a single equation that describes the rotation velocity curve of spiral galaxies is derived. The equation is also used to correlate the measured mass of the theorized, central supermassive black hole with other galaxy parameters. The equation adds one parameter for each mature galaxy and one parameter for each particle (atom nucleus and smaller) species to Newtonian dynamics. Rising, flat, and declining rotation velocity curves are explained without unknown dark matter. Ten previously unknown parametric relationships are discovered. Also, the CUM integrates several heretofore unrelated observations. The CUM model is in an early development stage and has already predicted new relationships among galaxy parameters.

A correlation between (1) the difference between the Tully-Fisher calculated distance and Cepheid calculated distance for a target galaxy and (2) the magnitude and distance of galaxies close to the target galaxy is described. The result is based on a sample of 31 galaxies that have published Cepheid distances with a wide range of characteristics and distances from 2.02 Mpc to 49.7 Mpc. The energy impinging on a target galaxy from neighboring galaxies is related to the residual between Tully-Fisher and Cepheid distance calculations of the target galaxy. These relations have four different zones dependent on the value of the impinging energy. The correlation coefficients are 0.77, 0.87, 0.81, and 0.98, respectively. This relationship is of interest not only for its ability to reduce the Tully-Fisher to Cepheid distance residual but also because it suggests neighboring galaxies have a strong influence on the apparent luminosity of galaxies and it suggests the Tully-Fisher relationship assumption that the intrinsic mass to intrinsic luminosity ratio is constant among galaxies is valid for z < 0.006.

The rotation velocity asymmetry v observed in spiral galaxy HI rotation curves linearly correlates with the effective potential force from the 10 closest neighboring galaxies normalized for the test particle mass and the gravitational constant. The magnitude of the potential force from a close galaxy is proportional to the luminosity of the close galaxy and inversely proportional to the square of the distance from the close galaxy to the target galaxy. The correlation coefficient is 0.99 and F test is 0.99. Also, the slope of the rotation curve in the disk region of a galaxy from rising to flat to declining is qualitatively correlated with increasing asymmetry and, hence, to the net force from other galaxies. The result is based on a sample of nine spiral galaxies with published Cepheid distances and rotation curves and with a wide range of characteristics. These relationships are interesting not only for their predictive power but also because (1) they suggest a galaxy's dynamics and the shape of its rotation curve are related to the potential energy and distance of neighboring galaxies, (2) they suggest a rising rotation curve in the disk region is intrinsic, (3) they suggest the Tully-Fisher relationship's assumption that the mass to intrinsic luminosity ratio is constant among galaxies is valid, and (4) they are inconsistent with MOND, dark matter, and the Linear Potential Model.