Recent simulations and observations of massive galaxy cluster evolution predict that the majority of stellar mass build up happens within cluster members by $z=2$, before cluster virialization. Protoclusters rich with dusty, star-forming galaxies (DSFGs) at $z>3$ are the favored candidate progenitors for these massive galaxy clusters at $z\sim0$. We present here the first study analyzing stellar emission along with cold dust and gas continuum emission in a spectroscopically confirmed $z=4.002$ protocluster core rich with DSFGs, the Distant Red Core (DRC). We combine new HST and Spitzer data with existing Gemini, Herschel, and ALMA observations to derive individual galaxy-level properties, and compare them to coeval field and other protocluster galaxies. All of the protocluster members are massive ($>10^{10}$ M$_\odot$), but not significantly more so than their coeval field counterparts. Within uncertainty, all are nearly indistinguishable from galaxies on the star-forming vs. stellar mass main-sequence relationship, and on the star formation efficiency plane. Assuming no future major influx of fresh gas, we estimate that these gaseous DSFGs will deplete their gas reservoirs in $\sim300$ Myr, becoming the massive quiescent ellipticals dominating cluster cores by $z\sim3$. Using various methodologies, we derive a total $z=4$ halo mass of $\sim10^{14}$ M$_\odot$, and estimate that the DRC will evolve to become an ultra-massive cluster core of mass $\gtrsim10^{15}$ M$_\odot$ by $z=0$.

We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1 data set). Using standard methodologies on the optical data only we make a significant (10 sigma) detection of the weak lensing signal (a shear power spectrum) due to the massive supercluster of galaxies in the targeted region. For the radio data we develop a new method to measure the shapes of galaxies from the interferometric data, and we construct a simulation pipeline to validate this method. We then apply this analysis to our radio observations, treating the e-MERLIN and VLA data independently. We achieve source densities of 0.5 per square arcmin in the VLA data and 0.06 per square arcmin in the e-MERLIN data, numbers which prove too small to allow a detection of a weak lensing signal in either the radio data alone or in cross-correlation with the optical data. Finally, we show preliminary results from a visibility-plane combination of the data from e-MERLIN and VLA which will be used for the forthcoming full SuperCLASS data release. This approach to data combination is expected to enhance both the number density of weak lensing sources available and the fidelity with which their shapes can be measured.

The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers.

We present optical and near-infrared imaging covering a $\sim$1.53 deg$^2$ region in the Super-Cluster Assisted Shear Survey (SuperCLASS) field, which aims to make the first robust weak lensing measurement at radio wavelengths. We derive photometric redshifts for $\approx$176,000 sources down to $i^\prime_{\rm AB}\sim24$ and present photometric redshifts for 1.4 GHz $e$-MERLIN and VLA detected radio sources found in the central 0.26 deg$^{2}$. We compile an initial catalog of 149 radio sources brighter than S$_{1.4}>75$ $\mu$Jy and find their photometric redshifts span $0<z_{\rm phot}<4$ with radio luminosities between $10^{21}-10^{25}$ W Hz$^{-1}$, with medians of $\langle z \rangle =0.55$ and $\langle L_{1.4}\rangle =1.9\times10^{23}$ W Hz$^{-1}$ respectively. We find 95% of the \uJy\ radio source sample (141/149) have SEDs best fit by star-forming templates while 5% (8/149) are better fit by AGN. Spectral indices are calculated for sources with radio observations from VLA and GMRT at 325 MHz, with an average spectral slope of $\alpha=0.59\pm0.04$. Using the full photometric redshift catalog we construct a density map at the redshift of the known galaxy clusters, $z=0.20\pm0.08$. Four of the five clusters are prominently detected at $>7 \sigma$ in the density map and we confirm the photometric redshifts are consistent with previously measured spectra from a few galaxies at the cluster centers.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

The main sequence of galaxies, a correlation between the star formation rates and stellar masses of galaxies, has been observed out to $z\sim4$. Galaxies within the scatter of the correlation are typically interpreted to be secularly evolving while galaxies with star formation rates elevated above the main sequence are interpreted to be undergoing interactions or to be Toomre-unstable disks with starbursting clumps. In this paper we investigate the recent merger histories of three dusty star forming galaxies, identified by their bright submillimeter emission at $z\sim1.5$. We analyze rest-frame optical and UV imaging, rest-frame optical emission line kinematics using slit spectra obtained with MOSFIRE on Keck I, and calculate Gini and M$_{20}$ statistics for each galaxy and conclude two are merger-driven while the third is an isolated disk galaxy. The disk galaxy lies $\sim$4$\times$ above the main sequence, one merger lies within the scatter of the main sequence, and one merger lies $\sim$4$\times$ below the main sequence. This hints that the location of a galaxy with respect to the main sequence may not be a useful discriminator of the recent star formation history of high-M$_{\star}$ galaxies at $z\sim1$.

We combine observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) to assess the redshift and to study the star formation conditions in AzTEC2: one of the brightest sub-millimeter galaxies (SMGs) in the COSMOS field ($S_{\rm 1.1mm}=10.5\pm1.4$mJy). Our high-resolution observations confirm that AzTEC2 splits into two components (namely AzTEC2-A and AzTEC2-B) for which we detect [C$\,$II] and $^{12}$CO(5$\to$4) line emission, implying a redshift of $4.626\pm0.001$ ($4.633\pm0.001$) for AzTEC2-A (AzTEC2-B) and ruling out previous associations with a galaxy at $z\sim1$. We use the $^{12}$CO(5$\to$4) line emission and adopt typical SMG-like gas excitation conditions to estimate the molecular gas mass, which is $M_{\rm gas}(\alpha_{\rm CO}/2.5)=2.1\pm0.4 \times10^{11}{\rm M}_\odot$ for AzTEC2-A, and a factor four lower for AzTEC2-B. With the infrared-derived star formation rate of AzTEC2-A ($1920\pm100 \,M_\odot{\rm \, yr}^{-1}$) and AzTEC2-B ($710\pm 35\,M_\odot{\rm \,yr}^{-1}$), they both will consume their current gas reservoir within $(30-200)\,$Myr. We find evidence of a rotation-dominated [C$\,$II] disk in AzTEC2-A, with a de-projected rotational velocity of $v_{\rm rot}(i=39^\circ)=660\pm130{\rm \,km\,s}^{-1}$, velocity dispersion $\lesssim100{\rm \,km\,s}^{-1}$, and dynamical mass of $M_{\rm dyn}(i=39^\circ)=2.6^{+1.2}_{-0.9}\times10^{11}\,M_\odot$. We propose that an elevated gas accretion rate from the cosmic web might be the main driver of the intense levels of star formation in AzTEC2-A, which might be further enhanced by gravitational torques induced by its minor companion (AzTEC2-B). These results strengthen the picture whereby the population of single-dish selected SMGs is rather heterogeneous, including a population of pairs of massive, highly-active galaxies in a pre-coalescence phase.

We present the highest resolution --- 15 pc (0.03'') --- ALMA $^{12}$CO(2-1) line emission and 1.3mm continuum maps, tracers of the molecular gas and dust, respectively, in the nearby merging galaxy system NGC 6240, that hosts two supermassive black holes growing simultaneously. These observations provide an excellent spatial match to existing Hubble optical and near-infrared observations of this system. A significant molecular gas mass, $\sim$9$\times$10$^9$M$_\odot$, is located in between the two nuclei, forming a clumpy stream kinematically dominated by turbulence, rather than a smooth rotating disk as previously assumed from lower resolution data. Evidence for rotation is seen in the gas surrounding the southern nucleus, but not in the northern one. Dynamical shells can be seen, likely associated with nuclear supernovae remnants. We further detect the presence of significant high velocity outflows, some of them reaching velocities $>$500 km/s, affecting a significant fraction, $\sim$11\% of the molecular gas in the nuclear region. Inside the spheres of influence of the northern and southern supermassive black holes we find molecular masses of 7.4$\times$10$^8$M$_\odot$ and 3.3$\times$10$^9$M$_\odot$, respectively. We are thus directly imaging the reservoir of gas that can accrete onto each supermassive black hole. These new ALMA maps highlight the critical need for high resolution observations of molecular gas in order to understand the feeding of supermassive black holes and its connection to galaxy evolution in the context of a major galaxy merger.

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 present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$\mu$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $\tau_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr).

We present ALMA Band 6 (nu=233GHz, lambda=1.3mm) continuum observations towards 68 'normal' star-forming galaxies within two Coma-like progenitor structures at z=2.10 and 2.47, from which ISM masses are derived, providing the largest census of molecular gas mass in overdense environments at these redshifts. Our sample comprises galaxies with a stellar mass range of 1x10^9M_sun - 4x10^11M_sun with a mean M_*~6x10^10M_sun. Combining these measurements with multiwavelength observations and SED modeling, we characterize the gas mass fraction and the star formation efficiency, and infer the impact of the environment on galaxies' evolution. Most of our detected galaxies (~70%) have star formation efficiencies and gas fractions similar to those found for coeval field galaxies and in agreement with the field scaling relations. However, we do find that the proto-clusters contain an increased fraction of massive, gas-poor galaxies, with low gas fractions (f_gas<6-10%) and red rest-frame ultraviolet/optical colors typical of post-starburst and passive galaxies. The relatively high abundance of passive galaxies suggests an accelerated evolution of massive galaxies in proto-cluster environments. The large fraction of quenched galaxies in these overdense structures also implies that environmental quenching takes place during the early phases of cluster assembly, even before virialization. From our data, we derive a quenching efficiency of E_q~0.45 and an upper limit on the quenching timescale of T_q<1Gyr.

We examine the role of environment on the in situ star formation (SF) hosted by the progenitors of the most massive galaxies in the present-day universe, the brightest cluster galaxies (BCGs), from $z \sim 3$ to present in the COSMOS field. Progenitors are selected from the COSMOS field using a stellar mass cut motivated by the evolving cumulative comoving number density of progenitors within the Illustris simulation, as well as the Millennium-II simulation and a constant comoving number density method for comparison. We characterize each progenitor using far-ultraviolet--far-infrared observations taken from the COSMOS field and fitting stellar, dust, and active galactic nucleus components to their spectral energy distributions. Additionally, we compare the SF rates of our progenitor sample to the local density maps of the COSMOS field to identify the effects of environment. We find that BCG progenitors evolve in three stages, starting with an in situ SF dominated phase ($z > 2.25$). This is followed by a phase until $z \sim 1.25$ where mass growth is driven by in situ SF and stellar mass deposited by mergers (both gas rich and poor) on the same order of magnitude independent of local environment. Finally, at low redshift dry mergers are the dominant stellar mass generation process. We also identify this final transition period as the time when progenitors quench, exhibiting quiescent NUV\emph{rJ} colors.

We test the use of long-wavelength dust continuum emission as a molecular gas tracer at high redshift, via a unique sample of 12, z~2 galaxies with observations of both the dust continuum and CO(1-0) line emission (obtained with the Atacama Large Millimeter Array and Karl G. Jansky Very Large Array, respectively). Our work is motivated by recent, high redshift studies that measure molecular gas masses (\ensuremath{\rm{M}_{\rm{mol}}}) via a calibration of the rest-frame $850\mu$m luminosity ($L_\mathrm{850\mu m,rest}$) against the CO(1-0)-derived \ensuremath{\rm{M}_{\rm{mol}}}\ of star-forming galaxies. We hereby test whether this method is valid for the types of high-redshift, star-forming galaxies to which it has been applied. We recover a clear correlation between the rest-frame $850\mu$m luminosity, inferred from the single-band, long-wavelength flux, and the CO(1-0) line luminosity, consistent with the samples used to perform the $850\mu$m calibration. The molecular gas masses, derived from $L_\mathrm{850\mu m,rest}$, agree to within a factor of two with those derived from CO(1-0). We show that this factor of two uncertainty can arise from the values of the dust emissivity index and temperature that need to be assumed in order to extrapolate from the observed frequency to the rest-frame at 850$\mathrm{\mu m}$. The extrapolation to 850$\mathrm{\mu m}$ therefore has a smaller effect on the accuracy of \Mmol\ derived via single-band dust-continuum observations than the assumed CO(1-0)-to-\ensuremath{\rm{M}_{\rm{mol}}}\ conversion factor. We therefore conclude that single-band observations of long-wavelength dust emission can be used to reliably constrain the molecular gas masses of massive, star-forming galaxies at $z\gtrsim2$.

We present a suite of 34 high-resolution cosmological zoom-in simulations consisting of thousands of halos up to M_halo~10^12 M_sun (M_star~10^10.5 M_sun) at z>=5 from the Feedback in Realistic Environments project. We post-process our simulations with a three-dimensional Monte Carlo dust radiative transfer code to study dust extinction, dust emission, and dust temperature within these simulated z>=5 galaxies. Our sample forms a tight correlation between infrared excess (IRX=F_IR/F_UV) and ultraviolet (UV)-continuum slope (beta_UV), despite the patchy, clumpy dust geometry shown in our simulations. We find that the IRX-beta_UV relation is mainly determined by the shape of the extinction curve and is independent of its normalization (set by the dust-to-gas ratio). The bolometric IR luminosity (L_IR) correlates with the intrinsic UV luminosity and the star formation rate (SFR) averaged over the past 10 Myr. We predict that at a given L_IR, the peak wavelength of the dust spectral energy distributions for z>=5 galaxies is smaller by a factor of 2 (due to higher dust temperatures on average) than at z=0. The higher dust temperatures are driven by higher specific SFRs and SFR surface densities with increasing redshift. We derive the galaxy UV luminosity functions (LFs) at z=5-10 from our simulations and confirm that a heavy attenuation is required to reproduce the observed bright-end UVLFs. We also predict the IRLFs and UV luminosity densities at z=5-10. We discuss the implications of our results on current and future observations probing dust attenuation and emission in z>=5 galaxies.

Dust offers a unique probe of the interstellar medium (ISM) across multiple size, density, and temperature scales. Dust is detected in outflows of evolved stars, star-forming molecular clouds, planet-forming disks, and even in galaxies at the dawn of the Universe. These grains also have a profound effect on various astrophysical phenomena from thermal balance and extinction in galaxies to the building blocks for planets, and changes in dust grain properties will affect all of these phenomena. A full understanding of dust in all of its forms and stages requires a multi-disciplinary investigation of the dust life cycle. Such an investigation can be achieved with a statistical study of dust properties across stellar evolution, star and planet formation, and redshift. Current and future instrumentation will enable this investigation through fast and sensitive observations in dust continuum, polarization, and spectroscopy from near-infrared to millimeter wavelengths.

We present deep continuum observations at a wavelength of 2mm centered on the COSMOS field using the Goddard IRAM Superconducting Millimeter Observer (GISMO) at the IRAM 30m-telescope. These data constitute the widest deep 2mm survey to-date, reaching a uniform $\sigma\sim0.23$mJy beam$^{-1}$ sensitivity over $\sim250$arcmin$^2$ at $\sim\,24''$ resolution. We detect four sources at high significance (S/N$\,\geq\,$4.4) with an expected number of false detection of 0.09 sources, and five sources at $4.4>\,$S/N$\,\geq\,3.7$ with an expected number of false detection of 1.65 sources. Combined with deep GISMO observations in GOODS-N, we constrain the 2mm number counts over one decade in flux density. These measurements agree with most galaxy evolution models tested here, except those with large population of dusty star-forming galaxies at $z>7$. Five GISMO sources have counterparts in (sub-)millimeter catalogs available in COSMOS. Their redshifts suggest that all but one lie above $z\sim3$. These four high-redshift ($z>3$) galaxies have $\tilde{z}=3.9$, SFRs $\sim\,$400 - 1200M$_{\odot}\,$yr$^{-1}$ and $M_{\rm dust}\sim10^{9.5}\,$M$_{\odot}$. They provide a relatively complete selection ($\sim66\%$) of the most luminous ($L_{\rm IR}>10^{12.6}\,$L$_{\odot}$) and highest redshift ($z>3$) galaxies detected within our survey area by AzTEC at 1.1mm. We thus conclude that 2mm surveys favor the selection of massive, vigorously star-forming, high-redshift galaxies. This is corroborated by GISMO-C4, a source with a low false detection probability ($\sim\,6.2\%$), for which the absence of a (sub-)millimeter counterpart supports a high redshift origin ($z\gtrsim3$).

There is wide consensus that galaxy outflows are one of the most important processes determining the evolution of galaxies through cosmic time, for example playing a key role in shaping the galaxy mass function. Our understanding of outflows and their drivers, however, is in its infancy --- this is particularly true for the cold (neutral atomic and molecular) phases of outflows, which present observational and modeling challenges. Here we outline several key open questions, briefly discussing the requirements of the observations necessary to make progress, and the relevance of several existing and planned facilities. It is clear that galaxy outflows, and particularly cold outflows, will remain a topic of active research for the next decade and beyond.

We present far-IR photometry and infrared spectrum of the z=3.9114 quasar/starburst composite system APM 08279+5255 obtained using the Stratospheric Observatory for Infrared Astronomy (SOFIA)/HAWC+ and the Spitzer Space Telescope Infrared Spectrograph (IRS). We decompose the IR-to-radio spectral energy distribution (SED), sampled in 51 bands, using (i) a model comprised of two-temperature modified blackbodies (MBB) and radio power-laws and (ii) a semi-analytic model, which also accounts for emission from a clumpy torus. The latter is more realistic but requires a well-sampled SED, as possible here. In the former model, we find temperatures of T_warm = 296^17_15 K and T_cold = 110^3_3 K for the warm and cold dust components, respectively. This model suggests that the cold dust component dominates the FIR energy budget (66%) but contributes only 17% to the total IR luminosity. Based on the torus models, we infer an inclination angle of i=15^8_8 degree and the presence of silicate emission, in accordance with the Type-1 active galactic nucleus nature of APM 08279+5255. Accounting for the torus' contribution to the FIR luminosity, we find a lensing-corrected star formation rate of SFR=3075x(4/mu_L) Msun yr^-1. We find that the central quasar contributes 30% to the FIR luminosity but dominates the total IR luminosity (93%). The 30% correction is in contrast to the 90% reported in previous work. In addition, the IR luminosity inferred from the torus model is a factor of two higher. These differences highlight the importance of adopting physically motivated models to properly account for IR emission in high-z quasars, which is now possible with SOFIA/HAWC+.

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.

We present the science case for mapping several thousand galaxy (proto)clusters at z=1-10 with a large aperture single dish sub-mm facility, producing a high-redshift counterpart to local large surveys of rich clusters like the well-studied Abell catalogue. Principal goals of a large survey of distant clusters are the evolution of galaxy clusters over cosmic time and the impact of environment on the evolution and formation of galaxies. To make a big leap forward in this emerging research field, the community would benefit from a large-format, wide-band, direct-detection spectrometer (e.g., based on MKID technology), covering a wide field of ~1 square degree and a frequency coverage from 70 to 700 GHz.

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.

The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure.

Galaxy evolution is driven by many complex interrelated processes as galaxies accrete gas, form new stars, grow their stellar masses and central black holes, and subsequently quench. The processes that drive these transformations is poorly understood, but it is clear that the local environment on multiple scales plays a significant role. Today's massive clusters are dominated by spheroidal galaxies with low levels of star formation while those in the field are mostly still actively forming their stars. In order to understand the physical processes that drive both the mass build up in galaxies and the quenching of star formation, we need to investigate galaxies and their surrounding gas within and around the precursors of today's massive galaxy clusters -- protoclusters at z>2. The transition period before protoclusters began to quench and become the massive clusters we observe today is a crucial time to investigate their properties and the mechanisms driving their evolution. However, until now, progress characterizing the galaxies within protoclusters has been slow, due the difficulty of obtaining highly complete spectroscopic observations of faint galaxies at z>2 over large areas of the sky. The next decade will see a transformational shift in our understanding of protoclusters as deep spectroscopy over wide fields of view will be possible in conjunction with high resolution deep imaging in the optical and near-infrared.

Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time.

The epoch of reionization (6 < z < 10) marks the period in our universe when the first large galaxies grew to fruition, and began to affect the universe around them. Massive stars, and potentially accreting supermassive black holes, filled the universe with ionizing radiation, burning off the haze of neutral gas that had filled the intergalactic medium (IGM) since recombination (z~1000). The evolution of this process constrains key properties of these earliest luminous sources, thus observationally constraining reionization is a key science goal for the next decade. The measurement of Lyman-alpha emission from photometrically-identified galaxies is a highly constraining probe of reionization, as a neutral IGM will resonantly scatter these photons, reducing detectability. While significant work has been done with 8-10m telescopes, these observations require extremely large telescopes (ELTs); the flux limits available from today's 10m class telescopes are sufficient for only the brightest known galaxies (m < 26). Ultra-deep surveys with the Giant Magellan Telescope (GMT) and Thirty Meter Telescope (TMT) will be capable of detecting Lyman-alpha emission from galaxies 2-3 magnitudes fainter than today's deepest surveys. Wide-field fiber spectroscopy on the GMT combined with narrow-field AO-assisted slit spectroscopy on the TMT will be able to probe the expected size of ionized bubbles throughout the epoch of reionization, following up degree scale deep imaging surveys with the Wide Field Infrared Space Telescope. These data will provide the first resolved Lyman-alpha-based maps of the ionized intergalactic medium throughout the epoch of reionization, constraining models of both the temporal and spatial evolution of this phase change.

Observations of the baryon to dark matter fraction in galaxies through cosmic time are a fundamental test for galaxy formation models. Recent observational studies have suggested that some disk galaxies at $z>1$ host declining rotation curves, in contrast with observations of low redshift disk galaxies where stellar or HI rotation curves flatten at large radii. We present an observational counterexample, a galaxy named DSFG850.95 at $z=1.555$ (4.1 Gyr after the big bang) that hosts a flat rotation curve between radii of $\sim$6--14 kpc (1.2--2.8 disk scale lengths) and has a dark matter fraction of $0.44\pm0.08$ at the H-band half light radius, similar to the Milky Way. We create position-velocity and position-dispersion diagrams using Keck/MOSFIRE spectroscopic observations of H$\alpha$ and [NII] emission features, which reveal a flat rotation velocity of $V_{\rm flat}=285\pm12$ km/s and an ionized gas velocity dispersion of $\sigma_{0}=48\pm4$ km/s. This galaxy has a rotation-dominated velocity field with $V_{\rm flat}/\sigma_{0}\sim6$. Ground-based H-band imaging reveals a disk with S\'ersic index of $1.29\pm0.03$, an edge-on inclination angle of $87\pm2^{\circ}$, and an H-band half light radius of $8.4\pm0.1$ kpc. Our results point to DSFG850.95 being a massive, rotationally-supported disk galaxy with a high dark-matter-to-baryon fraction in the outer galaxy, similar to disk galaxies at low redshift.

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 next generation Very Large Array (ngVLA) will revolutionize our understanding of the distant Universe via the detection of cold molecular gas in the first galaxies. Its impact on studies of galaxy characterization via detailed gas dynamics will provide crucial insight on dominant physical drivers for star-formation in high redshift galaxies, including the exchange of gas from scales of the circumgalactic medium down to resolved clouds on mass scales of $\sim10^{5}\,M_\odot$. In this study, we employ a series of high-resolution, cosmological, hydrodynamic zoom simulations from the MUFASA simulation suite and a CASA simulator to generate mock ngVLA observations of a $z\sim4.5$ gas rich star-forming galaxy. Using the DESPOTIC radiative transfer code that encompasses simultaneous thermal, chemical, and statistical equilibrium in calculating the molecular and atomic level transitions of CO from ALMA for comparison. We find that observations of CO(1-0) are especially important for tracing the systemic redshift of the galaxy and the total mass of the well-shielded molecular gas reservoir, while even CO(2-1) can predominantly trace denser gas regions distinct from CO(1-0). The factor of 100 times improvement in mapping speed for the ngVLA beyond the Jansky VLA and the proposed ALMA Band 1 will make these detailed, high-resolution imaging and kinematic studies of CO(1-0) routine at $z\sim2-5$.

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 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 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 science case and associated science requirements for a next-generation Very Large Array (ngVLA) are described, highlighting the five key science goals developed out of a community-driven vision of the highest scientific priorities in the next decade. Building on the superb cm observing conditions and existing infrastructure of the VLA site in the U.S. Southwest, the ngVLA is envisaged to be an interferometric array with more than 10 times the sensitivity and spatial resolution of the current VLA and ALMA, operating at frequencies spanning $\sim1.2 - 116$\,GHz with extended baselines reaching across North America. The ngVLA will be optimized for observations at wavelengths between the exquisite performance of ALMA at submm wavelengths, and the future SKA-1 at decimeter to meter wavelengths, thus lending itself to be highly complementary with these facilities. The ngVLA will be the only facility in the world that can tackle a broad range of outstanding scientific questions in modern astronomy by simultaneously delivering the capability to: (1) unveil the formation of Solar System analogues; (2) probe the initial conditions for planetary systems and life with astrochemistry; (3) characterize the assembly, structure, and evolution of galaxies from the first billion years to the present; (4) use pulsars in the Galactic center as fundamental tests of gravity; and (5) understand the formation and evolution of stellar and supermassive blackholes in the era of multi-messenger astronomy.

Bright high-redshift quasars ($z>6$), hosting supermassive black holes ($\rm M_{BH} > 10^8 M_{\odot}$), are expected to reside in massive host galaxies embedded within some of the earliest and most massive galaxy overdensities. We analyze 1.2\,mm ALMA dust continuum maps of 35 bright quasars at $6<z<7$ and search the primary beam for excess dust continuum emission from sources with L$_{\rm IR} \gtrsim{ 10^{12}}$\,L$_{\odot}$ as evidence for early protoclusters. We compare the detection rates of continuum sources at $\geq5\sigma$ significance in the fields surrounding the quasars (A$_{\rm eff}$ = 4.3\,arcmin$^2$) with millimeter number counts in blank field surveys. We discover 15 millimeter sources in the fields excluding the quasars themselves, corresponding to an overdensity $\delta_{\rm gal} \equiv (N_{\rm gal} - N_{\rm exp})/N_{\rm exp} = -0.07\pm0.56$, consistent with no detected overdensity of dusty galaxies within 140 physical kpc of the quasars. However, the apparent lack of continuum overdensity does not negate the hypothesis that quasars live in overdense environments, as evidenced by strong [CII] overdensities found on the same scales to similarly-selected quasars. The small field of view of ALMA could miss a true overdensity if it exists on scales larger than 1\,cMpc, if the quasar is not centered in the overdensity, or if quasar feedback plays a role close to the quasar, but it is most likely that the large line of sight volume probed by a continuum survey will wash out a true overdensity signal. We discuss the necessary factors in determining the bias with which dusty star-forming galaxies trace true dark matter overdensities in order to improve upon overdensity searches in dust continuum.

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.

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.

Interferometric observations have demonstrated that a significant fraction of single-dish submillimetre (submm) sources are blends of multiple submm galaxies (SMGs), but the nature of this multiplicity, i.e. whether the galaxies are physically associated or chance projections, has not been determined. We performed spectroscopy of 11 SMGs in six multi-component submm sources, obtaining spectroscopic redshifts for nine of them. For an additional two component SMGs, we detected continuum emission but no obvious features. We supplement our observed sources with four sources from the literature. This sample allows us to statistically constrain the physical nature of single-dish submm source multiplicity for the first time. In three [3/7, or 43 (-33/+39) per cent at 95% confidence] of the single-dish sources for which the nature of the blending is unambiguous, the components for which spectroscopic redshifts are available are physically associated, whereas 4/7 [57 (-39/+33) per cent] have at least one unassociated component. When components whose spectra exhibit continuum but no features and for which the photometric redshift is significantly different from the spectroscopic redshift of the other component are also considered, 6/9 [67 (-37/+26) per cent] of the single-dish sources are comprised of at least one unassociated component. The nature of the multiplicity of one single-dish source is ambiguous. We conclude that physically associated systems and chance projections both contribute to the multi-component single-dish submm source population. This result contradicts the conventional wisdom that bright submm sources are solely a result of merger-induced starbursts, as blending of unassociated galaxies is also important.

We present optical and near-IR Integral Field Unit (IFU) and ALMA band 6 observations of the nearby dual Active Galactic Nuclei (AGN) Mrk 463. At a distance of 210 Mpc, and a nuclear separation of $\sim$4 kpc, Mrk 463 is an excellent laboratory to study the gas dynamics, star formation processes and supermassive black hole (SMBH) accretion in a late-stage gas-rich major galaxy merger. The IFU observations reveal a complex morphology, including tidal tails, star-forming clumps, and emission line regions. The optical data, which map the full extent of the merger, show evidence for a biconical outflow and material outflowing at $>$600 km s$^{-1}$, both associated with the Mrk 463E nucleus, together with large scale gradients likely related to the ongoing galaxy merger. We further find an emission line region $\sim$11 kpc south of Mrk 463E that is consistent with being photoionized by an AGN. Compared to the current AGN luminosity, the energy budget of the cloud implies a luminosity drop in Mrk 463E by a factor 3-20 over the last 40,000 years. The ALMA observations of $^{12}$CO(2-1) and adjacent 1mm continuum reveal the presence of $\sim$10$^{9}$M$_\odot$ in molecular gas in the system. The molecular gas shows velocity gradients of $\sim$800 km/s and $\sim$400 km/s around the Mrk 463E and 463W nuclei, respectively. We conclude that in this system the infall of $\sim$100s $M_\odot$/yr of molecular gas is in rough balance with the removal of ionized gas by a biconical outflow being fueled by a relatively small, $<$0.01% of accretion onto each SMBH.

We present a deep survey of the SuperCLASS super-cluster - a region of sky known to contain five Abell clusters at redshift $z\sim0.2$ - performed using the Arcminute Microkelvin Imager (AMI) Large Array (LA) at 15.5$~$GHz. Our survey covers an area of approximately 0.9 square degrees. We achieve a nominal sensitivity of $32.0~\mu$Jy beam$^{-1}$ toward the field centre, finding 80 sources above a $5\sigma$ threshold. We derive the radio colour-colour distribution for sources common to three surveys that cover the field and identify three sources with strongly curved spectra - a high-frequency-peaked source and two GHz-peaked-spectrum sources. The differential source count (i) agrees well with previous deep radio source count, (ii) exhibits no evidence of an emerging population of star-forming galaxies, down to a limit of 0.24$~$mJy, and (iii) disagrees with some models of the 15$~$GHz source population. However, our source count is in agreement with recent work that provides an analytical correction to the source count from the SKADS Simulated Sky, supporting the suggestion that this discrepancy is caused by an abundance of flat-spectrum galaxy cores as-yet not included in source population models.

This document describes some of the fundamental astrophysical problems that require observing capabilities at millimeter- and centimeter wavelengths well beyond those of existing, or already planned, telescopes. The results summarized in this report follow a solicitation from the National Radio Astronomy Observatory to develop key science cases for a future U. S.-led radio telescope, the "next generation Very Large Array" (ngVLA). The ngVLA will have roughly 10 times the collecting area of the Jansky VLA, operate at frequencies from 1 GHz to 116 GHz with up to 20 GHz of bandwidth, possess a compact core for high surface-brightness sensitivity, and extended baselines of at least hundreds of kilometers and ultimately across the continent to provide high-resolution imaging. The ngVLA builds on the scientific and technical legacy of the Jansky VLA and ALMA, and will be designed to provide the next leap forward in our understanding of planets, galaxies, and black holes.

The total star formation budget of galaxies consists of the sum of the unobscured star formation, as observed in the rest-frame ultraviolet (UV), together with the obscured component that is absorbed and re-radiated by dust grains in the infrared. We explore how the fraction of obscured star formation depends on stellar mass for mass-complete samples of galaxies at $0<z<2.5$. We combine GALEX and WISE photometry for SDSS-selected galaxies with the 3D-HST treasury program and Spitzer/MIPS 24$\mu$m photometry in the well-studied 5 extragalactic CANDELS fields. We find a strong dependence of the fraction of obscured star formation ($f_{\mathrm{obscured}}$=SFR$_{\mathrm{IR}}$/SFR$_{\mathrm{UV+IR}}$) on stellar mass, with remarkably little evolution in this fraction with redshift out to $z$=2.5. 50\% of star formation is obscured for galaxies with log(M/M$_{\odot}$)=9.4; although unobscured star formation dominates the budget at lower masses, there exists a tail of low mass extremely obscured star-forming galaxies at $z>1$. For log(M/M$_{\odot}$)$>$10.5, $>$90\% of star formation is obscured at all redshifts. We also show that at fixed total SFR, $f_{\mathrm{obscured}}$ is lower at higher redshift. At fixed mass, high-redshift galaxies are observed to have more compact sizes and much higher star formation rates, gas fractions and hence surface densities (implying higher dust obscuration), yet we observe no redshift evolution in $f_{\mathrm{obscured}}$ with stellar mass. This poses a challenge to theoretical models to reproduce, where the observed compact sizes at high redshift seem in tension with lower dust obscuration.

Recent literature suggests that there are two modes through which galaxies grow their stellar mass - a normal mode characterized by quasi-steady star formation, and a highly efficient starburst mode possibly triggered by stochastic events such as galaxy mergers. While these differences are established for extreme cases, the population of galaxies in-between these two regimes is poorly studied and it is not clear where the transition between these two modes of star formation occurs. We utilize ALMA observations of the CO J=3-2 line luminosity in a sample of 20 infrared luminous galaxies that lie in the intermediate range between normal and starburst galaxies at z ~ 0.25-0.6 in the COSMOS field to examine the gas content and star formation efficiency of these galaxies. We compare these quantities to the galaxies' deviation from the well-studied "main sequence" correlation between star formation rate and stellar mass (MS) and find that at log($SFR/SFR_{MS}$) < 0.6, a galaxy's distance to the main sequence is mostly driven by increased gas content, and not a more efficient star formation process.

We carried out targeted ALMA observations of 129 fields in the COSMOS region at 1.25 mm, detecting 152 galaxies at S/N$\geq$5 with an average continuum RMS of 150 $\mu$Jy. These fields represent a S/N-limited sample of AzTEC / ASTE sources with 1.1 mm S/N$\geq$4 over an area of 0.72 square degrees. Given ALMA's fine resolution and the exceptional spectroscopic and multiwavelength photometric data available in COSMOS, this survey allows us unprecedented power in identifying submillimeter galaxy counterparts and determining their redshifts through spectroscopic or photometric means. In addition to 30 sources with prior spectroscopic redshifts, we identified redshifts for 113 galaxies through photometric methods and an additional nine sources with lower limits, which allowed a statistically robust determination of the redshift distribution. We have resolved 33 AzTEC sources into multi-component systems and our redshifts suggest that nine are likely to be physically associated. Our overall redshift distribution peaks at $z\sim$2.0 with a high redshift tail skewing the median redshift to $\tilde{z}$=2.48$\pm$0.05. We find that brighter millimeter sources are preferentially found at higher redshifts. Our faintest sources, with S$_{1.25 \rm mm}$<1.25 mJy, have a median redshift of $\tilde{z}$=2.18$\pm$0.09, while the brightest sources, S$_{1.25 \rm mm}$>1.8 mJy, have a median redshift of $\tilde{z}$=3.08$\pm$0.17. After accounting for spectral energy distribution shape and selection effects these results are consistent with several previous submillimeter galaxy surveys, and moreover, support the conclusion that the submillimeter galaxy redshift distribution is sensitive to survey depth.

We examine the rest-frame ultra-violet (UV) properties of 10 [CII]$\lambda158\,{\rm \mu m}$$-$detected galaxies at $z\sim5.5$ in COSMOS using new HST/WFC3 near-infrared imaging. Together with pre-existing $158\,{\rm \mu m}-$continuum and [CII] line measurements by ALMA, we study their dust attenuation properties on the IRX-$\beta$ diagram, which connects the total dust emission ($\propto$ IRX=log($L_{FIR}/L_{1600}$)) to the line-of-sight dust column ($\propto\beta$). We find systematically bluer UV continuum spectral slopes ($\beta$) compared to previous low-resolution ground-based measurements, which relieves some of the tension between models of dust attenuation and observations at high redshifts. While most of the galaxies are consistent with local starburst or Small Magellanic cloud like dust properties, we find galaxies with low IRX values and a large range in $\beta$ that cannot be explained by models of a uniform dust distribution well mixed with stars. A stacking analysis of Keck/DEIMOS optical spectra indicates that these galaxies are metal-poor with young stellar populations which could significantly alter their spatial dust distribution.

We determine the physical properties of a sample of SMGs in the COSMOS field that were pre-selected at the observed wavelength of $\lambda_{\rm obs}=1.1$ mm, and followed up at $\lambda_{\rm obs}=1.3$ mm with ALMA. We used MAGPHYS to fit the panchromatic (ultraviolet to radio) SEDs of 124 of the target SMGs, 19.4% of which are spectroscopically confirmed. The SED analysis was complemented by estimating the gas masses of the SMGs by using the $\lambda_{\rm obs}=1.3$ mm emission as a tracer of the molecular gas. The sample median and 16th-84th percentile ranges of the stellar masses, SFRs, dust temperatures, and dust and gas masses were derived to be $\log(M_{\star}/{\rm M}_{\odot})=11.09^{+0.41}_{-0.53}$, ${\rm SFR}=402^{+661}_{-233}$ ${\rm M}_{\odot}~{\rm yr}^{-1}$, $T_{\rm dust}=39.7^{+9.7}_{-7.4}$ K, $\log(M_{\rm dust}/{\rm M}_{\odot})=9.01^{+0.20}_{-0.31}$, and $\log(M_{\rm gas}/{\rm M}_{\odot})=11.34^{+0.20}_{-0.23}$, respectively. The median gas-to-dust ratio and gas fraction were found to be $120^{+73}_{-30}$ and $0.62^{+0.27}_{-0.23}$, respectively. We found that 57.3% of our SMGs populate the main sequence (MS) of star-forming galaxies, while 41.9% of the sources lie above the MS by a factor of >3 (one source lies below the MS). The largest 3 GHz radio sizes are found among the MS sources. Those SMGs that appear irregular in the rest-frame UV are predominantly starbursts, while the MS SMGs are mostly disk-like. The larger radio-emitting sizes of the MS SMGs compared to starbursts is a likely indication of their more widespread, less intense star formation. The irregular UV morphologies of the starburst SMGs are likely to echo their merger nature. Our results suggest that the transition from high-$z$ SMGs to local ellipticals via compact, quiescent galaxies (cQGs) at $z \sim 2$ might not be universal, and the latter population might also descend from the so-called blue nuggets.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 870um observations of 29 bright Herschel sources near high-redshift QSOs. The observations confirm that 20 of the Herschel sources are submillimeter-bright galaxies (SMGs) and identify 16 new SMG-QSO pairs that are useful to studies of the circumgalactic medium (CGM) of SMGs. Eight out of the 20 SMGs are blends of multiple 870um sources. The angular separations for six of the Herschel-QSO pairs are less than 10", comparable to the sizes of the Herschel beam and the ALMA primary beam. We find that four of these six "pairs" are actually QSOs hosted by SMGs. No additional submillimeter companions are detected around these QSOs and the rest-frame ultraviolet spectra of the QSOs show no evidence of significant reddening. Black hole accretion and star formation contribute almost equally in bolometric luminosity in these galaxies. The SMGs hosting QSOs show similar source sizes, dust surface densities, and SFR surface densities as other SMGs in the sample. We find that the black holes are growing $\sim$3$\times$ faster than the galaxies when compared to the present-day black-hole-galaxy mass ratio, suggesting a QSO duty cycle of $\lesssim$30% in SMGs at z ~ 3. The remaining two Herschel-detected QSOs are undetected at 870um but each has an SMG "companion" only 9" and 12" away (71 and 95 kpc at z = 3). They could be either merging or projected pairs. If the former, they would represent a rare class of "wet-dry" mergers. If the latter, the QSOs would, for the first time, probe the CGM of SMGs at impact parameters below 100 kpc.

We present the results of combined deep Keck/NIRC2, HST/WFC3 near-infrared and Herschel far infrared observations of an extremely star forming dusty lensed galaxy identified from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS J133542.9+300401). The galaxy is gravitationally lensed by a massive WISE identified galaxy cluster at $z\sim1$. The lensed galaxy is spectroscopically confirmed at $z=2.685$ from detection of $\rm {CO (1 \rightarrow 0)}$ by GBT and from detection of $\rm {CO (3 \rightarrow 2)}$ obtained with CARMA. We use the combined spectroscopic and imaging observations to construct a detailed lens model of the background dusty star-forming galaxy (DSFG) which allows us to study the source plane properties of the target. The best-fit lens model provide magnification of $\mu_{\rm star}=2.10\pm0.11$ and $\mu_{\rm dust}=2.02\pm0.06$ for the stellar and dust components respectively. Multi-band data yields a magnification corrected star formation rate of $1900(\pm200)\,M_{\odot}{\rm yr^{-1}}$ and stellar mass of $6.8_{-2.7}^{+0.9}\times10^{11}\,M_{\odot}$ consistent with a main sequence of star formation at $z\sim2.6$. The CO observations yield a molecular gas mass of $8.3(\pm1.0)\times10^{10}\,M_{\odot}$, similar to the most massive star-forming galaxies, which together with the high star-formation efficiency are responsible for the intense observed star formation rates. The lensed DSFG has a very short gas depletion time scale of $\sim40$ Myr. The high stellar mass and small gas fractions observed indicate that the lensed DSFG likely has already formed most of its stellar mass and could be a progenitor of the most massive elliptical galaxies found in the local Universe.

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.

ALMA observations of the long wavelength dust continuum are used to estimate the interstellar medium (ISM) masses in a sample of 708 galaxies at z = 0.3 to 4.5 in the COSMOS field. The galaxy sample has known far-infrared luminosities and, hence, star formation rates (SFRs), and stellar masses (M$_{\rm *}$) from the optical-infrared spectrum fitting. The galaxies sample SFRs from the main sequence (MS) to 50 times above the MS. The derived ISM masses are used to determine the dependence of gas mass on redshift, M$_{\rm *}$, and specific SFR (sSFR) relative to the MS. The ISM masses increase approximately 0.63 power of the rate of increase in SFRs with redshift and the 0.32 power of the sSFR/sSFR$_MS$. The SF efficiencies also increase as the 0.36 power of the SFR redshift evolutionary and the 0.7 power of the elevation above the MS; thus the increased activities at early epochs are driven by both increased ISM masses and SF efficiency. Using the derived ISM mass function we estimate the accretion rates of gas required to maintain continuity of the MS evolution ($>100$ \msun yr$^{-1}$ at z $>$ 2.5). Simple power-law dependences are similarly derived for the gas accretion rates. We argue that the overall evolution of galaxies is driven by the rates of gas accretion. The cosmic evolution of total ISM mass is estimated and linked to the evolution of SF and AGN activity at early epochs.

We present near-infrared and optical spectroscopic observations of a sample of 450$\mu$m and 850$\mu$m-selected dusty star-forming galaxies (DSFGs) identified in a 400 arcmin$^2$ area in the COSMOS field. Thirty-one sources of the 102 targets were spectroscopically confirmed at $0.2<z<4$, identified primarily in the near-infrared with Keck MOSFIRE and some in the optical with Keck LRIS and DEIMOS. The low rate of confirmation is attributable both to high rest-frame optical obscuration in our targets and limited sensitivity to certain redshift ranges. The high-quality photometric redshifts available in the COSMOS field allow us to test the robustness of photometric redshifts for DSFGs. We find a subset (11/31$\approx35$%) of DSFGs with inaccurate ($\Delta z/(1+z)>0.2$) or non-existent photometric redshifts; these have very distinct spectral energy distributions from the remaining DSFGs, suggesting a decoupling of highly obscured and unobscured components. We present a composite rest-frame 4300--7300\AA\ spectrum for DSFGs, and find evidence of 200$\pm$30 km s$^{-1}$ gas outflows. Nebular line emission for a sub-sample of our detections indicate that hard ionizing radiation fields are ubiquitous in high-z DSFGs, even more so than typical mass or UV-selected high-z galaxies. We also confirm the extreme level of dust obscuration in DSFGs, measuring very high Balmer decrements, and very high ratios of IR to UV and IR to H$\alpha$ luminosities. This work demonstrates the need to broaden the use of wide bandwidth technology in the millimeter to the spectroscopic confirmations of large samples of high-z DSFGs, as the difficulty in confirming such sources at optical/near-infrared wavelengths is exceedingly challenging given their obscuration.

We determine the radio size distribution of a large sample of 152 SMGs in COSMOS that were detected with ALMA at 1.3 mm. For this purpose, we used the observations taken by the VLA-COSMOS 3 GHz Large Project. One hundred and fifteen of the 152 target SMGs were found to have a 3 GHz counterpart. The median value of the major axis FWHM at 3 GHz is derived to be $4.6\pm0.4$ kpc. The radio sizes show no evolutionary trend with redshift, or difference between different galaxy morphologies. We also derived the spectral indices between 1.4 and 3 GHz, and 3 GHz brightness temperatures for the sources, and the median values were found to be $\alpha=-0.67$ and $T_{\rm B}=12.6\pm2$ K. Three of the target SMGs, which are also detected with the VLBA, show clearly higher brightness temperatures than the typical values. Although the observed radio emission appears to be predominantly powered by star formation and supernova activity, our results provide a strong indication of the presence of an AGN in the VLBA and X-ray-detected SMG AzTEC/C61. The median radio-emitting size we have derived is 1.5-3 times larger than the typical FIR dust-emitting sizes of SMGs, but similar to that of the SMGs' molecular gas component traced through mid-$J$ line emission of CO. The physical conditions of SMGs probably render the diffusion of cosmic-ray electrons inefficient, and hence an unlikely process to lead to the observed extended radio sizes. Instead, our results point towards a scenario where SMGs are driven by galaxy interactions and mergers. Besides triggering vigorous starbursts, galaxy collisions can also pull out the magnetised fluids from the interacting disks, and give rise to a taffy-like synchrotron-emitting bridge. This provides an explanation for the spatially extended radio emission of SMGs, and can also cause a deviation from the well-known IR-radio correlation.

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.

A well-instrumented Green Bank Telescope (GBT) operating at high frequency represents a unique scientific resource for the US community. As a filled-aperture, 100m-diameter telescope, the GBT is ideally suited to fast mapping of extended, low surface brightness emission with excellent instantaneous frequency coverage. This capability makes the GBT a key facility for a range of cutting edge science described in this document, only possible at these frequencies. We note that the ability to perform the necessary observations is unique and highly complementary to the capabilities offered by interferometers, and should be preserved. We argue that rather than divesting from this exceptional resource, it makes sense for the US community to invest moderately to maintain GBT operations and to instrument it in an optimal manner, enabling it to become an extraordinary complement to existing and future radio interferometers. Adequately instrumented, the GBT would be a pillar for 20-115 GHz science in the US and the world.

We present the first results from an on-going survey to characterize the circumgalactic medium (CGM) of the massive high-redshift galaxies detected as submillimeter galaxies (SMGs). We constructed a parent sample of 163 SMG-QSO pairs with separations less than $\sim$36 arcsec by cross-matching far-infrared-selected galaxies from Herschel with spectroscopically confirmed QSOs. The Herschel sources were selected to match the properties of SMGs. We determined the sub-arcsecond positions of six Herschel sources with the Very Large Array and obtained secure redshift identification for three of those with near-infrared spectroscopy. The QSO sightlines probe transverse proper distances of 112, 157, and 198 kpc at foreground redshifts of 2.043, 2.515, and 2.184, respectively, which are comparable to the virial radius of the $\sim10^{13}$ Msun halos expected to host SMGs. High-quality absorption-line spectroscopy of the QSOs reveals systematically strong HI Lyman-alpha absorption around all three SMGs, with rest-frame equivalent widths of $\sim2-3$ \AA. However, none of the three absorbers exhibits compelling evidence for optically thick HI gas or metal absorption, in contrast to the dominance of strong neutral absorbers in the CGM of luminous $z \sim 2$ QSOs. The low covering factor of optically thick HI gas around SMGs tentatively indicates that SMGs may not have as prominent cool gas reservoirs in their halos as the co-eval QSOs and that they may inhabit less massive halos than previously thought.

We present the results of 325 MHz GMRT observations of a super-cluster field, known to contain five Abell clusters at redshift $z \sim 0.2$. We achieve a nominal sensitivity of $34\,\mu$Jy beam$^{-1}$ toward the phase centre. We compile a catalogue of 3257 sources with flux densities in the range $183\,\mu\rm{Jy}\,-\,1.5\,\rm{Jy}$ within the entire $\sim 6.5$ square degree field of view. Subsequently, we use available survey data at other frequencies to derive the spectral index distribution for a sub-sample of these sources, recovering two distinct populations -- a dominant population which exhibit spectral index trends typical of steep-spectrum synchrotron emission, and a smaller population of sources with typically flat or rising spectra. We identify a number of sources with ultra-steep spectra or rising spectra for further analysis, finding two candidate high-redshift radio galaxies and three gigahertz-peaked-spectrum radio sources. Finally, we derive the Euclidean-normalised differential source counts using the catalogue compiled in this work, for sources with flux densities in excess of $223 \, \mu$Jy. Our differential source counts are consistent with both previous observations at this frequency and models of the low-frequency source population. These represent the deepest source counts yet derived at 325 MHz. Our source counts exhibit the well-known flattening at mJy flux densities, consistent with an emerging population of star-forming galaxies; we also find marginal evidence of a downturn at flux densities below $308 \, \mu$Jy, a feature so far only seen at 1.4 GHz.

We present near-infrared integral-field spectroscopic observations targeting H$\alpha$ in eight sub-millimeter galaxies (SMGs) at $z$=1.3-2.5 using VLT/SINFONI, obtaining significant detections for six of them. The star formation rates derived from the H$\alpha$ emission are $\sim$100 M$_\odot$yr$^{-1}$, which account for only $\sim$ 20-30\% of the infrared-derived values, thus suggesting that these systems are very dusty. Two of these systems present [NII]/H$\alpha$ ratios indicative of the presence of an Active Galactic Nucleus (AGN). We mapped the spatial distribution and kinematics of the star forming regions in these galaxies on kpc-scales. In general, the H$\alpha$ morphologies tend to be highly irregular and/or clumpy, showing spatial extents of $\sim$3-11~kpc. We find evidence for significant spatial offsets, of $\sim$0.1-0.4$"$ or 1.2-3.4 kpc, between the H$\alpha$ and the continuum emission in three of the sources. Performing a kinemetry analysis we conclude that the majority of the sample is not consistent with disk-like rotation-dominated kinematics. Instead, they tend to show irregular and/or clumpy and turbulent velocity and velocity dispersion fields. This can be interpreted as evidence for scenario in which these extreme star formation episodes are triggered by galaxy-galaxy interactions and major mergers. In contrast to recent results for SMGs, these sources appear to follow the same relations between gas and star forming rate densities as less luminous and/or normal star forming galaxies.

The most prodigious starburst galaxies are absent in massive galaxy clusters today, but their connection with large scale environments is less clear at $z\gtrsim2$. We present a search of large scale structure around a galaxy cluster core at $z=2.095$ using a set of spectroscopically confirmed galaxies. We find that both color-selected star-forming galaxies (SFGs) and dusty star-forming galaxies (DSFGs) show significant overdensities around the $z=2.095$ cluster. A total of 8 DSFGs (including 3 X-ray luminous active galactic nuclei, AGNs) and 34 SFGs are found within a 10 arcmin radius (corresponds to $\sim$15 cMpc at $z\sim2.1$) from the cluster center and within a redshift range of $\Delta z=0.02$, which leads to galaxy overdensities of $\delta_{\rm DSFG}\sim12.3$ and $\delta_{\rm SFG}\sim2.8$. The cluster core and the extended DSFG- and SFG-rich structure together demonstrate an active cluster formation phase, in which the cluster is accreting a significant amount of material from large scale structure while the more mature core may begin to virialize. Our finding of this DSFG-rich structure, along with a number of other protoclusters with excess DSFGs and AGNs found to date, suggest that the overdensities of these rare sources indeed trace significant mass overdensities. However, it remains puzzling how these intense star formers are triggered concurrently. Although an increased probability of galaxy interactions and/or enhanced gas supply can trigger the excess of DSFGs, our stacking analysis based on 850 $\mu$m images and morphological analysis based on rest-frame optical imaging do not show such enhancements of merger fraction and gas content in this structure.

ALMA Cycle 2 observations of the long wavelength dust emission in 145 star-forming galaxies are used to probe the evolution of star-forming ISM. We also develop the physical basis and empirical calibration (with 72 low-z and z ~ 2 galaxies) for using the dust continuum as a quantitative probe of interstellar medium (ISM) masses. The galaxies with highest star formation rates (SFRs) at <z> = 2.2 and 4.4 have gas masses up to 100 times that of the Milky Way and gas mass fractions reaching 50 to 80%, i.e. gas masses 1 - 4 times their stellar masses. We find a single high-z star formation law: SFR = 35 M_ mol^0.89 x (1+z)_{z=2}^0.95 x (sSFR)_{MS}^0.23 \msun yr^-1 -- an approximately linear dependence on the ISM mass and an increased star formation efficiency per unit gas mass at higher redshift. Galaxies above the Main Sequence (MS) have larger gas masses but are converting their ISM into stars on a timescale only slightly shorter than those on the MS -- thus these 'starbursts' are largely the result of having greatly increased gas masses rather than and increased efficiency for converting gas to stars. At z $> 1$, the entire population of star-forming galaxies has $\sim$ 2 - 5 times shorter gas depletion times than low-z galaxies. These shorter depletion times indicate a different mode of star formation in the early universe -- most likely dynamically driven by compressive, high-dispersion gas motions -- a natural consequence of the high gas accretion rates.

The Universe's largest galaxy clusters likely built the majority of their massive $>10^{11} M_\odot$ galaxies in simultaneous, short-lived bursts of activity well before virialization. This conclusion is reached from emerging datasets on $z>2$ proto-clusters and the characteristics of their member galaxies, in particular, rare starbursts and ultraluminous active galactic nuclei (AGN). The most challenging observational hurdle in identifying such structures is their very large volumes, $\sim10^{4}$ comoving Mpc$^3$ at $z>2$, subtending areas $\sim$half a degree on the sky. Thus the contrast afforded by an overabundance of very rare galaxies in comparison to the background can more easily distinguish overdense structures from the surrounding, normal density field. Five $2<z<3$ proto-clusters from the literature are discussed in detail and are found to contain up to 12 dusty starbursts or luminous AGN galaxies each, a phenomenon that is unlikely to occur by chance even in overdense environments. Measurements of gas depletion times suggest that they are indeed short-lived on $\sim$100 Myr timescales, and the probability of finding a structure containing more than 8 such systems is $\sim$0.2\%, unless their `triggering' is correlated on very large spatial scales, $\sim$10 Mpc across. The volume density of starburst-rich proto-clusters is found to be comparable to all $>10^{15} M_\odot$ galaxy clusters in the nearby Universe, a factor of five larger than expected in some simulations. Some tension yet exists between measurements of their volume density of starburst-rich proto-clusters and the expectation that they are generated via short-lived episodes. However, improved observations of proto-clusters over large regions of sky will certainly shed more light on the assembly of galaxy clusters, and whether or not they build their galaxies through episodic bursts as suggested here. [abridged]

ALMA Cycle 2 observations of the long wavelength dust emission in 180 star-forming (SF) galaxies are used to investigate the evolution of ISM masses at z = 1 to 6.4. The ISM masses exhibit strong increases from z = 0 to $\rm <z>$ = 1.15 and further to $\rm <z>$ = 2.2 and 4.8, particularly amongst galaxies above the SF galaxy main sequence (MS). The galaxies with highest SFRs at $\rm <z>$ = 2.2 and 4.8 have gas masses 100 times that of the Milky Way and gas mass fractions reaching 50 to 80\%, i.e. gas masses 1 - 4$\times$ their stellar masses. For the full sample of galaxies, we find a single, very simple SF law: $\rm SFR \propto M_{\rm ISM}^{0.9}$, i.e. a `linear' dependence on the ISM mass -- on and above the MS. Thus, the galaxies above the MS are converting their larger ISM masses into stars on a timescale similar to those on the MS. At z $> 1$, the entire population of star-forming galaxies has $\sim$5 - 10$\times$ shorter gas depletion times ($\sim0.2$ Gyr) than galaxies at low redshift. These {\bf shorter depletion times are due to a different, dominant mode of SF in the early universe} -- dynamically driven by compressive, high dispersion gas motions and/or galaxy interactions. The dispersive gas motions are a natural consequence of the extraordinarily high gas accretion rates which must occur to maintain the prodigious SF.

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.

We describe the search for Lyman-break galaxies (LBGs) near the sub-millimeter bright starburst galaxy HFLS3 at $z$$=$6.34 and a study on the environment of this massive galaxy during the end of reionization.We performed two independent selections of LBGs on images obtained with the \textit{Gran Telescopio Canarias} (GTC) and the \textit{Hubble Space Telescope} (HST) by combining non-detections in bands blueward of the Lyman-break and color selection. A total of 10 objects fulfilling the LBG selection criteria at $z$$>$5.5 were selected over the 4.54 and 55.5 arcmin$^2$ covered by our HST and GTC images, respectively. The photometric redshift, UV luminosity, and the star-formation rate of these sources were estimated with models of their spectral energy distribution. These $z$$\sim$6 candidates have physical properties and number densities in agreement with previous results. The UV luminosity function at $z$$\sim$6 and a Voronoi tessellation analysis of this field shows no strong evidence for an overdensity of relatively bright objects (m$_{F105W}$$<$25.9) associated with \textit{HFLS3}. However, the over-density parameter deduced from this field and the surface density of objects can not excluded definitively the LBG over-density hypothesis. Moreover we identified three faint objects at less than three arcseconds from \textit{HFLS3} with color consistent with those expected for $z$$\sim$6 galaxies. Deeper data are needed to confirm their redshifts and to study their association with \textit{HFLS3} and the galaxy merger that may be responsible for the massive starburst.

Numerical simulations of cosmological structure formation show that the Universe's most massive clusters, and the galaxies living in those clusters, assemble rapidly at early times (2.5 < z < 4). While more than twenty proto-clusters have been observed at z > 2 based on associations of 5-40 galaxies around rare sources, the observational evidence for rapid cluster formation is weak. Here we report observations of an asymmetric, filamentary structure at z = 2.47 containing seven starbursting, submillimeter-luminous galaxies and five additional AGN within a comoving volume of 15000 Mpc$^{3}$. As the expected lifetime of both the luminous AGN and starburst phase of a galaxy is ~100 Myr, we conclude that these sources were likely triggered in rapid succession by environmental factors, or, alternatively, the duration of these cosmologically rare phenomena is much longer than prior direct measurements suggest. The stellar mass already built up in the structure is $\sim10^{12}M_{\odot}$ and we estimate that the cluster mass will exceed that of the Coma supercluster at $z \sim 0$. The filamentary structure is in line with hierarchical growth simulations which predict that the peak of cluster activity occurs rapidly at z > 2.

Many of the transformative processes in the Universe have taken place in regions obscured by dust, and are best studied with far-IR spectroscopy. We present the Cryogenic-Aperture Large Infrared-Submillimeter Telescope Observatory (CALISTO), a 5-meter class, space-borne telescope actively cooled to 4 K, emphasizing moderate-resolution spectroscopy in the crucial 35 to 600 micron band. CALISTO will enable NASA and the world to study the rise of heavy elements in the Universe's first billion years, chart star formation and black hole growth in dust-obscured galaxies through cosmic time, and conduct a census of forming planetary systems in our region of the Galaxy. CALISTO will capitalize on rapid progress in both format and sensitivity of far-IR detectors. Arrays with a total count of a few 100,000 detector pixels will form the heart of a suite of imaging spectrometers in which each detector reaches the photon background limit. This document contains a large overview paper on CALISTO, as well as six 2-3 page scientific white papers, all prepared in response to NASA's Cosmic Origins Program Analysis Group (COPAG's) request for input on future mission concepts. The Far-IR Science Interest Group will meet from 3-5 June 2015 with the intention of reaching consensus on the architecture for the Far-IR Surveyor mission. This white paper describes one of the architectures to be considered by the community. One or more companion papers will describe alternative architectures.

We present the rest-frame optical spectral energy distribution and stellar masses of six Herschel- selected gravitationally lensed dusty, star-forming galaxies (DSFGs) at 1 < z < 3. These galaxies were first identified with Herschel/SPIRE imaging data from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). The targets were observed with Spitzer/IRAC at 3.6 and 4.5um. Due to the spatial resolution of the IRAC observations at the level of 2 arcseconds, the lensing features of a background DSFG in the near-infrared are blended with the flux from the foreground lensing galaxy in the IRAC imaging data. We make use of higher resolution Hubble/WFC3 or Keck/NIRC2 Adaptive Optics imaging data to fit light profiles of the foreground lensing galaxy (or galaxies) as a way to model the foreground components, in order to successfully disentangle the foreground lens and background source flux densities in the IRAC images. The flux density measurements at 3.6 and 4.5um, once combined with Hubble/WFC3 and Keck/NIRC2 data, provide important constraints on the rest-frame optical spectral energy distribution of the Herschel-selected lensed DSFGs. We model the combined UV- to millimeter-wavelength SEDs to establish the stellar mass, dust mass, star-formation rate, visual extinction, and other parameters for each of these Herschel-selected DSFGs. These systems have inferred stellar masses in the range 8 x 10^10 to 4 x 10^11 Msun and star-formation rates of around 100 Msun yr-1. This puts these lensed sub-millimeter systems well above the SFR-M* relation observed for normal star-forming galaxies at similar redshifts. The high values of SFR inferred for these systems are consistent with a major merger-driven scenario for star formation.

We present Hubble Space Telescope (HST) WFC3 imaging and grism spectroscopy observations of the {\it Herschel}-selected gravitationally-lensed starburst galaxy HATLASJ1429-0028. The lensing system consists of an edge-on foreground disk galaxy at $z=0.218$ with a nearly complete Einstein ring of the infrared luminous galaxy at $z=1.027$. The WFC3 spectroscopy with G102 and G141 grisms, covering the wavelength range of 0.8 to 1.7 $\mu$m, resulted in detections of H$\alpha$+[NII], H$\beta$, [SII], and [OIII] for the background galaxy from which we measure line fluxes and ratios. The Balmer line ratio H$\alpha$/H$\beta$ of $7.5 \pm 4.4$, when corrected for [NII], results in an extinction for the starburst galaxy of $E(B-V)=0.8 \pm 0.5$. The $H\alpha$ based star-formation rate, when corrected for extinction, is $60 \pm 50$ M$_{\odot}$ yr$^{-1}$, lower than the instantaneous star-formation rate of 390 $\pm$ 90 M$_{\odot}$ yr$^{-1}$ from the total IR luminosity. We also compare the nebular line ratios of HATLASJ1429-0028 with other star-forming and sub-mm bright galaxies. The nebular line ratios are consistent with an intrinsic ultra-luminous infrared galaxy with no evidence for excitation by an active galactic nucleus (AGN). We estimate the metallicity, $12 + log(O/H)$, of HATLASJ1429-0028 to be 8.49 $\pm$ 0.16. Such a low value is below the average relations for stellar mass vs. metallicity of galaxies at $z \sim 1$ for a galaxy with stellar mass of $\sim 2 \times 10^{11}$ M$_{\odot}$. The combination of high stellar mass, lack of AGN indicators, low metallicity, and the high star-formation rate of HATLASJ1429-0028 suggest that this galaxy is currently undergoing a rapid formation.

Evolution in the measured rest frame ultraviolet spectral slope and ultraviolet to optical flux ratios indicate a rapid evolution in the dust obscuration of galaxies during the first 3 billion years of cosmic time (z>4). This evolution implies a change in the average interstellar medium properties, but the measurements are systematically uncertain due to untested assumptions, and the inability to measure heavily obscured regions of the galaxies. Previous attempts to directly measure the interstellar medium in normal galaxies at these redshifts have failed for a number of reasons with one notable exception. Here we report measurements of the [CII] gas and dust emission in 9 typical (~1-4L*) star-forming galaxies ~1 billon years after the big bang (z~5-6). We find these galaxies have >12x less thermal emission compared with similar systems ~2 billion years later, and enhanced [CII] emission relative to the far-infrared continuum, confirming a strong evolution in the interstellar medium properties in the early universe. The gas is distributed over scales of 1-8 kpc, and shows diverse dynamics within the sample. These results are consistent with early galaxies having significantly less dust than typical galaxies seen at z<3 and being comparable to local low-metallicity systems.

The classification of galaxy mergers and isolated disks is key for understanding the relative importance of galaxy interactions and secular evolution during the assembly of galaxies. The kinematic properties of galaxies as traced by emission lines have been used to suggest the existence of a significant population of high-z star-forming galaxies consistent with isolated rotating disks. However, recent studies have cautioned that post-coalescence mergers may also display disk-like kinematics. To further investigate the robustness of merger/disk classifications based on kinematic properties, we carry out a systematic classification of 24 local (U)LIRGs spanning a range of galaxy morphologies: from isolated spiral galaxies, ongoing interacting systems, to fully merged remnants. We artificially redshift the WiFeS observations of these local (U)LIRGs to z=1.5 to make a realistic comparison with observations at high-z, and also to ensure that all galaxies have the same spatial sampling of ~900 pc. Using both kinemetry-based and visual classifications, we find that the reliability of kinematic classification shows a strong trend with the interaction stage of galaxies. Mergers with two nuclei and tidal tails have the most distinct kinematic properties compared to isolated disks, whereas a significant population of the interacting disks and merger remnants are indistinguishable from isolated disks. The high fraction of late-stage mergers showing disk-like kinematics reflects the complexity of the dynamics during galaxy interactions. However, the exact fractions of misidentified disks and mergers depend on the definition of kinematic asymmetries and the classification threshold when using kinemetry-based classifications. Our results suggest that additional indicators such as morphologies traced by stars or molecular gas are required to further constrain the merger/disk classifications at high-z.

The relationship between galaxy star formation rates (SFR) and stellar masses ($M_\ast$) is re-examined using a mass-selected sample of $\sim$62,000 star-forming galaxies at $z \le 1.3$ in the COSMOS 2-deg$^2$ field. Using new far-infrared photometry from $Herschel$-PACS and SPIRE and $Spitzer$-MIPS 24 $\mu$m, along with derived infrared luminosities from the NRK method based on galaxies' locations in the restframe color-color diagram $(NUV - r)$ vs. $(r - K)$, we are able to more accurately determine total SFRs for our complete sample. At all redshifts, the relationship between median $SFR$ and $M_\ast$ follows a power-law at low stellar masses, and flattens to nearly constant SFR at high stellar masses. We describe a new parameterization that provides the best fit to the main sequence and characterizes the low mass power-law slope, turnover mass, and overall scaling. The turnover in the main sequence occurs at a characteristic mass of about $M_{0} \sim 10^{10} M_{\odot}$ at all redshifts. The low mass power-law slope ranges from 0.9-1.3 and the overall scaling rises in SFR as a function of $(1+z)^{4.12 \pm 0.10}$. A broken power-law fit below and above the turnover mass gives relationships of $SFR \propto M_{*}^{0.88 \pm 0.06}$ below the turnover mass and $SFR \propto M_{*}^{0.27 \pm 0.04}$ above the turnover mass. Galaxies more massive than $M_\ast \gtrsim 10^{10}\ M_{\rm \odot}$ have on average, a much lower specific star formation rate (sSFR) than would be expected by simply extrapolating the traditional linear fit to the main sequence found for less massive galaxies.

We present stringent constraints on the average mid-, far-infrared and radio emissions of $\sim$14200 quiescent galaxies (QGs), identified out to $z=3$ in the COSMOS field via their rest-frame NUV$-$r and r$-$J colors, and with stellar masses $M_{\star}=10^{9.8-12.2} \,M_{\odot} $. Stacking in deep Spitzer (MIPS $24\,\mu$m), Herschel (PACS and SPIRE), and VLA (1.4 GHz) maps reveals extremely low dust-obscured star formation rates for QGs (SFR $<0.1-3\,M_{\odot}$yr$^{-1}$ at $z \leqslant 2$ and $<6-18\,M_{\odot}$yr$^{-1}$ at $z > 2$), consistent with the low unobscured SFRs ($<0.01-1.2\,M_{\odot}$yr$^{-1}$) inferred from modeling their ultraviolet-to-near-infrared photometry. The average SFRs of QGs are $>10\times$ below those of star-forming galaxies (SFGs) within the $M_{\star}$- and $z$-ranges considered. The stacked 1.4 GHz signals (S/N $> 5$) are, if attributed solely to star formation, in excess of the total (obscured plus unobscured) SFR limits, suggestive of a widespread presence of low-luminosity active galactic nuclei (AGN) among QGs. Our results reaffirm the existence of a significant population QGs out to $z = 3$, thus corroborating the need for powerful quenching mechanism(s) to terminate star formation in galaxies at earlier epochs.

We present Keck-Adaptive Optics and Hubble Space Telescope high resolution near-infrared (IR) imaging for 500 um-bright candidate lensing systems identified by the Herschel Multi-tiered Extra-galactic Survey (HerMES) and Herschel Astrophysical Terahertz Survey (H-ATLAS). Out of 87 candidates with near-IR imaging, 15 (~17%) display clear near-IR lensing morphologies. We present near-IR lens models to reconstruct and recover basic rest-frame optical morphological properties of the background galaxies from 12 new systems. Sources with the largest near-IR magnification factors also tend to be the most compact, consistent with the size bias predicted from simulations and pre- vious lensing models for sub-millimeter galaxies. For four new sources that also have high-resolution sub-mm maps, we test for differential lensing between the stellar and dust components and find that the 880 um magnification factor (u_880) is ~1.5 times higher than the near-IR magnification factor (u_NIR), on average. We also find that the stellar emission is ~2 times more extended in size than dust. The rest-frame optical properties of our sample of Herschel-selected lensed SMGs are consistent with those of unlensed SMGs, which suggests that the two populations are similar.

Galaxies' rest-frame ultraviolet (UV) properties are often used to directly infer the degree to which dust obscuration affects the measurement of star formation rates. While much recent work has focused on calibrating dust attenuation in galaxies selected at rest-frame ultraviolet wavelengths, locally and at high-$z$, here we investigate attenuation in dusty, star-forming galaxies (DSFGs) selected at far-infrared wavelengths. By combining multiwavelength coverage across 0.15--500\,$\mu$m in the COSMOS field, in particular making use of {\it Herschel} imaging, and a rich dataset on local galaxies, we find a empirical variation in the relationship between rest-frame UV slope ($\beta$) and ratio of infrared-to-ultraviolet emission ($L_{\rm IR}/L_{\rm UV}\equiv\,IRX$) as a function of infrared luminosity, or total star formation rate, SFR. Both locally and at high-$z$, galaxies above SFR$\gt$50\,M$_\odot$\,yr$^{-1}$ deviate from the nominal $IRX-\beta$ relation towards bluer colors by a factor proportional to their increasing IR luminosity. We also estimate contamination rates of DSFGs on high-$z$ dropout searches of $\ll1$\%\ at $z\lt4-10$, providing independent verification that contamination from very dusty foreground galaxies is low in LBG searches. Overall, our results are consistent with the physical interpretation that DSFGs, e.g. galaxies with $>50$\,M$_\odot$\,yr$^{-1}$, are dominated at all epochs by short-lived, extreme burst events, producing many young O and B stars that are primarily, yet not entirely, enshrouded in thick dust cocoons. The blue rest-frame UV slopes of DSFGs are inconsistent with the suggestion that most DSFGs at $z\sim2$ exhibit steady-state star formation in secular disks.

Ultraluminous and luminous infrared galaxies (ULIRGs and LIRGs) are the most extreme star-forming galaxies in the universe, and dominate the total star formation rate density at z>1. In the local universe (z<0.3), the majority of ULIRGs and a significant portion of LIRGs are triggered by interactions between gas-rich spiral galaxies, yet it is unclear if this is still the case at high-z. To investigate the relative importance of galaxy interactions in infrared luminous galaxies, we carry out a comparison of optical morphological properties between local (U)LIRGs and (U)LIRGs at z=0.5-1.5 based on the same sample selection, morphology classification scheme, and optical morphology at similar rest-frame wavelengths. In addition, we quantify the systematics in comparing local and high-z datasets by constructing a redshifted dataset from local (U)LIRGs, in which its data quality mimics the high-z dataset. Based on the Gini-M20 classification scheme, we find that the fraction of interacting systems decreases by ~8% from local to z<~1, and it is consistent with the reduction between local and redshifted datasets (6(+14-6)%). Based on visual classifications, the merger fraction of local ULIRGs is found to be ~20% lower compared to published results, and the reduction due to redshifiting is 15(+10-8)%. Consequently, the differences of merger fractions between local and z<~1 (U)LIRGs is only ~17%. These results demonstrate that there is no strong evolution in the fraction of (U)LIRGs classified as mergers at least out to z~1. At z>1, the morphology types of ~30% of (U)LIRGs can not be determined due to their faintness in the F814W-band, and thus the merger fraction measured at z>1 suffers from large uncertainties.

We discuss the restframe UV emission from the starbursting galaxy HFLS3 at z=6.34, discovered in Herschel/SPIRE data due to its red color in the submm wavelengths from 250-500 um. The apparent inst. SFR of HFLS3 inferred from the total FIR luminosity measured with over 15 photometric data points between 100 to 1000 um is 2900 Msun/yr. Keck/NIRC2 Ks band adaptive optics imaging data showed two potential NIR counterparts near HFLS3. Previously, the northern galaxy was taken to be in the foreground at z=2.1 while the southern galaxy was assumed to HFLS3's NIR counterpart. New HST/WFC3 and ACS imaging data show both optically bright galaxies are in the foreground at z<6. A new lensing model based on HST data and mm-wave continuum emission yields a magnification of 2.2+/-0.3. The lack of multiple imaging constrains the lensing magnification to be lower than either 2.7 or 3.5 at the 95% confidence level for the two scenarios, which attribute one or two components to HFLS3 in the source plane. Correcting for gravitational lensing, the inst. SFR is 1320 Msun/yr with the 95% confidence lower limit around 830 Msun/yr. Using models for the restframe UV to FIR SED, the ave. SFR over the last 100 Myr is around 660 Msun/yr. The dust and stellar masses of HFLS3 from the same SED models are 3x10^8 Msun and ~5x10^10 Msun, respectively, with large systematic uncertainties on assumptions related to the SED model. With HST/WFC3 images we also find diffuse NIR emission about 0.5" (~3 kpc) SW of HFLS3 that remains undetected in the ACS data. The emission has a photometric redshift consistent with either z~6 or a dusty galaxy template at z~2. If at the same redshift as HFLS3 the detected diffuse emission could be part of the complex merger system that could be triggering the starburst. Alternatively, it could be part of the foreground structure at z~2.1 that is responsible for lensing of HFLS3.

Far-infrared and submillimeter wavelength surveys have now established the important role of dusty, star-forming galaxies (DSFGs) in the assembly of stellar mass and the evolution of massive galaxies in the Universe. The brightest of these galaxies have infrared luminosities in excess of 10$^{13}$ L$_{\odot}$ with implied star-formation rates of thousands of solar masses per year. They represent the most intense starbursts in the Universe, yet many are completely optically obscured. Their easy detection at submm wavelengths is due to dust heated by ultraviolet radiation of newly forming stars. When summed up, all of the dusty, star-forming galaxies in the Universe produce an infrared radiation field that has an equal energy density as the direct starlight emission from all galaxies visible at ultraviolet and optical wavelengths. The bulk of this infrared extragalactic background light emanates from galaxies as diverse as gas-rich disks to mergers of intense starbursting galaxies. Major advances in far-infrared instrumentation in recent years, both space-based and ground-based, has led to the detection of nearly a million DSFGs, yet our understanding of the underlying astrophysics that govern the start and end of the dusty starburst phase is still in nascent stage. This review is aimed at summarizing the current status of DSFG studies, focusing especially on the detailed characterization of the best-understood subset (submillimeter galaxies, who were summarized in the last review of this field over a decade ago, Blain et al., 2002), but also the selection and characterization of more recently discovered DSFG populations. We review DSFG population statistics, their physical properties including dust, gas and stellar contents, their environments, and current theoretical models related to the formation and evolution of these galaxies.

We use the James Clerk Maxwell Telescope's SCUBA-2 camera to image a 400 arcmin^2 area surrounding the GOODS-N field. The 850 micron rms noise ranges from a value of 0.49 mJy in the central region to 3.5 mJy at the outside edge. From these data, we construct an 850 micron source catalog to 2 mJy containing 49 sources detected above the 4-sigma level. We use an ultradeep (11.5 uJy at 5-sigma) 1.4 GHz image obtained with the Karl G. Jansky Very Large Array together with observations made with the Submillimeter Array to identify counterparts to the submillimeter galaxies. For most cases of multiple radio counterparts, we can identify the correct counterpart from new and existing Submillimeter Array data. We have spectroscopic redshifts for 62% of the radio sources in the 9 arcmin radius highest sensitivity region (556/894) and 67% of the radio sources in the GOODS-N region (367/543). We supplement these with a modest number of additional photometric redshifts in the GOODS-N region (30). We measure millimetric redshifts from the radio to submillimeter flux ratios for the unidentified submillimeter sample, assuming an Arp 220 spectral energy distribution. We find a radio flux dependent K-z relation for the radio sources, which we use to estimate redshifts for the remaining radio sources. We determine the star formation rates (SFRs) of the submillimeter sources based on their radio powers and their submillimeter and find that they agree well. The radio data are deep enough to detect star-forming galaxies with SFRs >2000 solar masses per year to z~6. We find galaxies with SFRs up to ~6,000 solar masses per year over the redshift range z=1.5-6, but we see evidence for a turn-down in the SFR distribution function above 2000 solar masses per year.

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

We investigate the claim that all dwarf spheroidal galaxies (dSphs) reside within halos that share a common, universal mass profile as has been derived for dSphs of the Galaxy. By folding in kinematic information for 25 Andromeda dSphs, more than doubling the previous sample size, we find that a singular mass profile can not be found to fit all the observations well. Further, the best-fit dark matter density profile measured for solely the Milky Way dSphs is marginally discrepant (at just beyond the 1 sigma level) with that of the Andromeda dSphs, where a profile with lower maximum circular velocity, and hence mass, is preferred. The agreement is significantly better when three extreme Andromeda outliers, And XIX, XXI and XXV, all of which have large half-light radii (>600pc) and low velocity dispersions (sigma_v < 5km/s) are omitted from the sample. We argue that the unusual properties of these outliers are likely caused by tidal interactions with the host galaxy.

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.

We combine Herschel PACS and SPIRE maps of the full 2 deg^2 COSMOS field with existing multi-wavelength data to obtain template and model-independent optical-to-far-infrared spectral energy distributions (SEDs) for 4,218 Herschel-selected sources with log(L_IR/L_sun) = 9.4-13.6 and z = 0.02-3.54. Median SEDs are created by binning the optical to far-infrared (FIR) bands available in COSMOS as a function of infrared luminosity. Herschel probes rest-frame wavelengths where the bulk of the infrared radiation is emitted, allowing us to more accurately determine fundamental dust properties of our sample of infrared luminous galaxies. We find that the SED peak wavelength (lambda_peak) decreases and the dust mass (M_dust) increases with increasing total infrared luminosity (L_IR). In the lowest infrared luminosity galaxies (log(L_IR/L_sun) = 10.0-11.5), we see evidence of Polycyclic Aromatic Hydrocarbons (PAH) features (lambda=7-9 um), while in the highest infrared luminosity galaxies (L_IR > 10^12 L_sun) we see an increasing contribution of hot dust and/or power-law emission, consistent with the presence of heating from an active galactic nucleus (AGN). We study the relationship between stellar mass and star formation rate of our sample of infrared luminous galaxies and find no evidence that Herschel-selected galaxies follow the SFR/M_* "main sequence" as previously determined from studies of optically selected, star-forming galaxies. Finally, we compare the mid-infrared (MIR) to FIR properties of our infrared luminous galaxies using the previously defined diagnostic, IR8 = L_IR / L_8, and find that galaxies with L_IR > 10^11.3 L_sun tend to systematically lie above (~3-5) the IR8 "infrared main sequence", suggesting either suppressed PAH emission or an increasing contribution from AGN heating.

Galaxy interactions/mergers have been shown to dominate the population of IR luminous galaxies (log(LIR)>11.6Lsun) in the local Universe (z<0.25). Recent studies based on the relation between galaxies' star formation rates and stellar mass (the SFR-M relation or the galaxy main sequence (MS)) have suggested that galaxy interaction/mergers may only become significant when galaxies fall well above the galaxy MS. Since the typical SFR at given M increases with redshift, the existence of galaxy MS implies that massive, IR-luminous galaxies at high-z may not necessarily be driven by galaxy interactions. We examine the role of galaxy interactions in the SFR-M relation by carrying out a morphological analysis of 2084 Herschel-selected galaxies at 0.2 < z < 1.5 in the COSMOS field. Herschel-PACS and -SPIRE observations covering the full 2-deg^2 COSMOS field provide one of the largest far-IR selected samples of high-redshift galaxies with well-determined redshifts to date, with sufficient sensitivity at z ~ 1, to sample objects lying on and above the galaxy MS. Using a detailed visual classification scheme, we show that the fraction of "disk galaxies" decreases and the fraction of "irregular" galaxies increases systematically with increasing LIR out to z ~ 1.5 and z ~ 1.0, respectively. At log(LIR) > 11.5 Lsun, >50% of the objects show evident features of strongly interacting/merger systems, where this percentage is similar to the studies of local IR-luminous galaxies. The fraction of interacting/merger systems also systematically increases with the deviation from the SFR-M relation, supporting the view that galaxies fall above the MS are more dominated by mergers than the MS galaxies. Meanwhile, we find that ~18% of massive IR-luminous MS galaxies are classified as interacting systems, where this population may not evolve through the evolutionary track predicted by a simple gas exhaustion model.

We use the SCUBA-2 submillimeter camera mounted on the JCMT to obtain extremely deep number counts at 450 and 850um. We combine data on two cluster lensing fields, A1689 and A370, and three blank fields, CDF-N, CDF-S, and COSMOS, to measure the counts over a wide flux range at each wavelength. We use statistical fits to broken power law representations to determine the number counts. This allows us to probe to the deepest possible level in the data. At both wavelengths our results agree well with the literature in the flux range over which they have been measured, with the exception of the 850um counts in CDF-S, where we do not observe the counts deficit found by previous single-dish observations. At 450um, we detect significant counts down to ~1mJy, an unprecedented depth at this wavelength. By integrating the number counts above this flux limit, we measure 113.9^{+49.7}_{-28.4} Jydeg^{-2} of the 450um extragalactic background light (EBL). The majority of this contribution is from sources with S_450um between 1-10mJy, and these sources are likely to be the ones that are analogous to the local luminous infrared galaxies (LIRGs). At 850um, we measure 37.3^{+21.1}_{-12.9} Jydeg^{-2} of the EBL. Because of the large systematic uncertainties on the COBE measurements, the percentage of the EBL we resolve could range from 48%-153% (44%-178%) at 450 (850)um. Based on high-resolution SMA observations of around half of the 4sigma 850um sample in CDF-N, we find that 12.5^{+12.1}_{-6.8}% of the sources are blends of multiple fainter sources. This is a low multiple fraction, and we find no significant difference between our original SCUBA-2 850um counts and the multiplicity corrected counts.

We present deep 450um and 850um observations of a large, uniformly covered 394arcmin^2 area in the COSMOS field obtained with the SCUBA-2 instrument on the James Clerk Maxwell Telescope (JCMT). We achieve root-mean-square noise values of 4.13mJy at 450um and 0.80mJy at 850um. The differential and cumulative number counts are presented and compared to similar previous works. Individual point sources are identified at >3.6sigma significance, a threshold corresponding to a 3-5% sample contamination rate. We identify 78 sources at 450um and 99 at 850um, with flux densities S450=13-37mJy and S850=2-16mJy. Only 62-76% of 450um sources are 850um detected and 61-81% of 850um sources are 450um detected. The positional uncertainties at 450um are small (1-2.5") and therefore allow a precise identification of multiwavelength counterparts without reliance on detection at 24um or radio wavelengths; we find that only 44% of 450um-selected galaxies and 60% of 850um-sources have 24um or radio counterparts. 450um-selected galaxies peak at <z>=1.95+-0.19 and 850um=selected galaxies peak at <z>=2.16+-0.11. The two samples occupy similar parameter space in redshift and luminosity, while their median SED peak wavelengths differ by ~10-50um (translating to deltaTdust =8-12K, where 450um-selected galaxies are warmer). The similarities of the 450um and 850um populations, yet lack of direct overlap between them, suggests that submillimeter surveys conducted at any single far-infrared wavelength will be significantly incomplete (~>30%) at censusing infrared-luminous star formation at high-z.

Dust-obscured galaxies (DOGs) are a UV-faint, IR-bright galaxy population that reside at z~2 and are believed to be in a phase of dusty star-forming and AGN activity. We present far-IR observations of a complete sample of DOGs in the 2 deg2 of COSMOS. The 3077 DOGs have <z>=1.9+/-0.3 and are selected from 24um and r+ observations using a color cut of r+ - [24]>=7.5 (AB mag) and S24>=100uJy. Based on the near-IR SEDs, 47% are star-formation dominated and 10% are AGN-dominated. We use SPIRE far-IR photometry from HerMES to calculate the IR luminosity and characteristic dust temperature for the 1572 (51%) DOGs that are detected at 250um (>=3{\sigma}). For the remaining 1505 (49%) that are undetected, we perform a median stacking analysis to probe fainter luminosities. Detected and undetected DOGs have average IR luminosities of (2.8+/-0.4) x 1012 LSun and (0.77+-0.08)x10^12LSun, and dust temperatures of 34+/-7 K and 37+/-3 K, respectively. The IR luminosity function of DOGs with S24>=100uJy is calculated, using far-IR observations and stacking. DOGs contribute 10-30% to the total star formation rate density of the Universe at z=1.5-2.5, dominated by 250um detected and bump DOGs. For comparison, DOGs contribute 30% to the star-formation rate density of the Universe for all galaxies at z=1.5-2.5 with S24>=100uJy. DOGs have a large scatter about the star-formation main sequence and their specific star-formation rates show that the observed phase of star-formation could be responsible for their observed stellar mass at z~2.

We present the results from a large 850 micron survey of the sigma Orionis cluster using the SCUBA-2 camera on the James Clerk Maxwell Telescope. The 0.5-degree diameter circular region we surveyed contains 297 young stellar objects with an age estimated at about 3Myr. We detect 9 of these objects, 8 of which have infrared excesses from an inner disc. We also serendipitously detect 3 non-stellar sources at > 5sigma that are likely background submillimetre galaxies. The 9 detected stars have inferred disc masses ranging from 5 to about 17MJup, assuming similar dust properties as Taurus discs and an ISM gas-to-dust ratio of 100. There is a net positive signal toward the positions of the individually undetected infrared excess sources indicating a mean disc mass of 0.5 MJup . Stacking the emission toward those stars without infrared excesses constrains their mean disc mass to less than 0.3MJup, or an equivalent Earth mass in dust. The submillimetre luminosity distribution is significantly different from that in the younger Taurus region, indicating disc mass evolution as star forming regions age and the infrared excess fraction decreases. Submillimeter Array observations reveal CO emission toward 4 sources demonstrating that some, but probably not much, molecular gas remains in these relatively evolved discs. These observations provide new constraints on the dust and gas mass of protoplanetary discs during the giant planet building phase and provide a reference level for future studies of disc evolution.

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

Stellar archeology shows that massive elliptical galaxies today formed rapidly about ten billion years ago with star formation rates above several hundreds solar masses per year (M_sun/yr). Their progenitors are likely the sub-millimeter-bright galaxies (SMGs) at redshifts (z) greater than 2. While SMGs' mean molecular gas mass of 5x10^10 M_sun can explain the formation of typical elliptical galaxies, it is inadequate to form ellipticals that already have stellar masses above 2x10^11 M_sun at z ~ 2. Here we report multi-wavelength high-resolution observations of a rare merger of two massive SMGs at z = 2.3. The system is currently forming stars at a tremendous rate of 2,000 M_sun/yr. With a star formation efficiency an order-of-magnitude greater than that of normal galaxies, it will quench the star formation by exhausting the gas reservoir in only ~200 million years. At a projected separation of 19 kiloparsecs, the two massive starbursts are about to merge and form a passive elliptical galaxy with a stellar mass of ~4x10^11 M_sun. Our observations show that gas-rich major galaxy mergers, concurrent with intense star formation, can form the most massive elliptical galaxies by z ~ 1.5.

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.

We present a homogeneous kinematic analysis of red giant branch stars within 18 of the 28 Andromeda dwarf spheroidal (dSph) galaxies, obtained using the Keck I LRIS and Keck II DEIMOS spectrographs. Based on their g-i colors (taken with the CFHT MegaCam imager), physical positions on the sky, and radial velocities, we assign probabilities of dSph membership to each observed star. Using this information, the velocity dispersions, central masses and central densities of the dark matter halos are calculated for these objects, and compared with the properties of the Milky Way dSph population. We also measure the average metallicity ([Fe/H]) from the co-added spectra of member stars for each M31 dSph and find that they are consistent with the trend of decreasing [Fe/H] with luminosity observed in the Milky Way population. We find that three of our studied M31 dSphs appear as significant outliers in terms of their central velocity dispersion, And XIX, XXI and XXV, all of which have large half-light radii (>700 pc) and low velocity dispersions (sigma_v<5 km/s). In addition, And XXV has a mass-to-light ratio within its half-light radius of just [M/L]_{half}=10.3^{+7.0}_{-6.7}, making it consistent with a simple stellar system with no appreciable dark matter component within its 1 sigma uncertainties. We suggest that the structure of the dark matter halos of these outliers have been significantly altered by tides.

We present the first analysis of the all-sky Swift BAT ultra hard X-ray (14-195 keV) data for a targeted list of objects. We find the BAT data can be studied at 3x fainter limits than in previous blind detection catalogs based on prior knowledge of source positions and using smaller energy ranges for source detection. We determine the AGN fraction in 134 nearby (z<0.05) luminous infrared galaxies (LIRGS) from the GOALS sample. We find that LIRGs have a higher detection frequency than galaxies matched in stellar mass and redshift at 14-195 keV and 24-35 keV. In agreement with work at other wavelengths, the AGN detection fraction increases strongly at high IR luminosity with half of high luminosity LIRGs (50%, 6/12, log L_IR/L_sun>11.8) detected. The BAT AGN classification shows 97% (37/38) agreement with Chandra and XMM AGN classification using hardness ratios or detection of a iron K-alpha line. This confirms our statistical analysis and supports the use of the Swift BAT all-sky survey to study fainter populations of any category of sources in the ultra hard X-ray band. BAT AGN in LIRGs tend to show higher column densities with 40\pm9% showing 14-195 keV/2-10 keV hardness flux ratios suggestive of high or Compton-thick column densities (log N_H>24 cm^-2), compared to only 12\pm5% of non-LIRG BAT AGN. We also find that using specific energy ranges of the BAT detector can yield additional sources over total band detections with 24% (5/21) of detections in LIRGs at 24-35 keV not detected at 14-195 keV.

We present the results of SCUBA2 observations at 450 micron and 850 micron of the field lensed by the massive cluster A370. With a total survey area > 100 arcmin2 and 1 sigma sensitivities of 3.92 and 0.82 mJy/beam at 450 and 850 micron respectively, we find a secure sample of 20 sources at 450 micron and 26 sources at 850 micron with a signal-to-noise ratio > 4. Using the latest lensing model of A370 and Monte Carlo simulations, we derive the number counts at both wavelengths. The 450 micron number counts probe a factor of four deeper than the counts recently obtained from the Herschel Space Telescope at similar wavelengths, and we estimate that ~47-61% of the 450 micron extragalactic background light (EBL) resolved into individual sources with 450 micron fluxes greater than 4.5 mJy. The faint 450 micron sources in the 4 sigma sample have positional accuracies of 3 arcseconds, while brighter sources (signal-to-noise > 6 sigma) are good to 1.4 arcseconds. Using the deep radio map (1 sigma ~ 6 uJy) we find that the percentage of submillimeter sources having secure radio counterparts is 85% for 450 micron sources with intrinsic fluxes > 6 mJy and 67% for 850 micron sources with intrinsic fluxes > 4 mJy. We also find that 67% of the > 4 sigma 450 micron sources are detected at 850 micron, while the recovery rate at 450 micron of > 4 sigma 850 micron sources is 54%. Combined with the source redshifts estimated using millimetric flux ratios, the recovered rate is consistent with the scenario where both 450 micron and 20 cm emission preferentially select lower redshift dusty sources, while 850 micron emission traces a higher fraction of dusty sources at higher redshifts. [Abridge]

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

We present Keck spectroscopic observations and redshifts for a sample of 767 Herschel-SPIRE selected galaxies (HSGs) at 250, 350, and 500um, taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey (HerMES) peaks at z=0.85, with 731 sources at z<2 and a tail of sources out to z~5. We measure more significant disagreement between photometric and spectroscopic redshifts (<delta_z>/(1+z)>=0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. We estimate that the vast majority (72-83%) of z<2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1mm. We estimate the luminosity function and implied star-formation rate density contribution of HSGs at z<1.6 and find overall agreement with work based on 24um extrapolations of the LIRG, ULIRG and total infrared contributions. This work significantly increased the number of spectroscopically confirmed infrared-luminous galaxies at z>>0 and demonstrates the growing importance of dusty starbursts for galaxy evolution studies and the build-up of stellar mass throughout cosmic time. [abridged]

We present spectroscopic observations for a sample of 36 Herschel-SPIRE 250-500um selected galaxies (HSGs) at 2<z<5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Redshifts are confirmed as part of a large redshift survey of Herschel-SPIRE-selected sources covering ~0.93deg^2 in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts at 24um or 1.4GHz using a cross-identification likelihood matching method. Individual source luminosities range from log(L_IR/Lsun)=12.5-13.6 (corresponding to star formation rates 500-9000Msun/yr, assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well characterized across large areas of sky in contrast to most z>2 HyLIRGs in the literature which are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at (7+-2)x10^(-3)Msun/yr h^3Mpc^(-3) at z~2.5, which is >10% of the estimated total star formation rate density (SFRD) of the Universe from optical surveys. The contribution at z~4 has a lower limit of 3x10^(-3)Msun/yr h^3 Mpc^(-3), ~>20% of the estimated total SFRD. This highlights the importance of extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, even at the earliest epochs.

Luminous and ultraluminous infrared galaxies ((U)LIRGs) are the most extreme star forming galaxies in the universe. The local (U)LIRGs provide a unique opportunity to study their multi-wavelength properties in detail for comparison to their more numerous counterparts at high redshifts. We present common large aperture photometry at radio through X-ray wavelengths, and spectral energy distributions (SEDs) for a sample of 53 nearby LIRGs and 11 ULIRGs spanning log (LIR/Lsun) = 11.14-12.57 from the flux-limited Great Observatories All-sky LIRG Survey (GOALS). The SEDs for all objects are similar in that they show a broad, thermal stellar peak and a dominant FIR thermal dust peak, where nuLnu(60um) / nuLnu(V) increases from ~2-30 with increasing LIR. When normalized at IRAS-60um, the largest range in the luminosity ratio, R(lambda)=log[nuLnu(lambda)/nuLnu(60um)] observed over the full sample is seen in the Hard X-rays (HX=2-10 keV). A small range is found in the Radio (1.4GHz), where the mean ratio is largest. Total infrared luminosities, LIR(8-1000um), dust temperatures, and dust masses were computed from fitting thermal dust emission modified blackbodies to the mid-infrared (MIR) through submillimeter SEDs. The new results reflect an overall ~0.02 dex lower luminosity than the original IRAS values. Total stellar masses were computed by fitting stellar population synthesis models to the observed near-infrared (NIR) through ultraviolet (UV) SEDs. Mean stellar masses are found to be log(M/Msun) = 10.79+/-0.40. Star formation rates have been determined from the infrared (SFR_IR~45Msun/yr) and from the monochromatic UV luminosities (SFR_UV~1.3Msun/yr), respectively. Multiwavelength AGN indicators have be used to select putative AGN: about 60% of the ULIRGs would have been classified as an AGN by at least one of the selection criteria.

We investigate the properties (e.g. star formation rate, dust attentuation, stellar mass and metallicity) of a sample of infrared luminous galaxies at z \sim 1 via near-IR spectroscopy with Subaru-FMOS. Our sample consists of Herschel SPIRE and Spitzer MIPS selected sources in the COSMOS field with photometric redshifts in the range 0.7 < z-phot < 1.8, which have been targeted in 2 pointings (0.5 sq. deg.) with FMOS. We find a modest success rate for emission line detections, with candidate H{\alpha} emission lines detected for 57 of 168 SPIRE sources (34 per cent). By stacking the near-IR spectra we directly measure the mean Balmer decrement for the H{\alpha} and H{\beta} lines, finding a value of <E(B-V)> = 0.51\pm0.27 for <LIR> = 10^12 Lsol sources at <z> = 1.36. By comparing star formation rates estimated from the IR and from the dust uncorrected H{\alpha} line we find a strong relationship between dust attenuation and star formation rate. This relation is broadly consistent with that previously seen in star-forming galaxies at z ~ 0.1. Finally, we investigate the metallicity via the N2 ratio, finding that z ~ 1 IR-selected sources are indistinguishable from the local mass-metallicity relation. We also find a strong correlation between dust attentuation and metallicity, with the most metal-rich IR-sources experiencing the largest levels of dust attenuation.

Spectral Energy Distribution (SED) fitting in the far-infrared (FIR) is greatly limited by a dearth of data and an excess of free parameters - from galaxies' dust composition, temperature, mass, orientation, opacity, to heating from AGN. This paper presents a simple FIR SED fitting technique joining a modified, single dust temperature greybody, representing the reprocessed starburst emission in the whole galaxy, to a mid-infrared powerlaw, which approximates hot-dust emission from AGN heating or clumpy, hot starbursting regions. This FIR SED can be used to measure infrared luminosities, dust temperatures and dust masses for both local and high-z galaxies with 3 to 10+ FIR photometric measurements. This fitting method is compared to infrared template SEDs in the literature using photometric data on 65 local luminous and ultraluminous infrared galaxies, (U)LIRGs. Despite relying only on 2-4 free parameters, the coupled greybody/powerlaw SED fitting described here produces better fits to photometric measurements than best-fit literature template SEDs (with residuals a factor of ~2 lower). A mean emissivity index of beta=1.60+-0.38 and mid-infrared powerlaw slope of alpha=2.0+-0.5 is measured; the former agrees with the widely presumed emissivity index of beta=1.5 and the latter is indicative of an optically-thin dust medium with a shallow radial density profile, ~r^-0.5. Adopting characteristic dust temperature as the inverse wavelength where the SED peaks, dust temperatures ~25-45K are measured for local (U)LIRGs, ~5-15K colder than previous estimates using only simple greybodies. This comparative study highlights the impact of SED fitting assumptions on the measurement of physical properties such as infrared luminosity (and thereby infrared-based star formation rate), dust temperature and dust mass, for both local and high-redshift galaxies. [abridged]

We present two-dimensional, integral field spectroscopy covering the rest-frame wavelengths of strong optical emission lines in nine sub-mm-luminous galaxies (SMGs) at 2.0<z<2.7. The GEMINI-NIFS and VLT-SINFONI imaging spectroscopy allows the mapping of the gas morphologies and dynamics within the sources, and we measure an average Halpha velocity dispersion of sigma=220+-80km/s and an average half light radius of r=3.7+-0.8kpc. The average dynamical measure, V_obs/2sigma=0.9+-0.1 for the SMGs, is higher than in more quiescent star-forming galaxies at the same redshift, highlighting a difference in the dynamics of the two populations. The SMGs' far-infrared SFRs, measured using Herschel-SPIRE far-infrared photometry, are on average 370+-90Mo/yr which is ~2 times higher than the extinction corrected SFRs of the more quiescent star-forming galaxies. Six of the SMGs in our sample show strong evidence for kinematically distinct multiple components with average velocity offsets of 200+-100km/s and average projected spatial offsets of 8+-2kpc, which we attribute to systems in the early stages of major mergers. Indeed all SMGs are classified as mergers from a kinemetry analysis of the velocity and dispersion field asymmetry. We bring together our sample with seven other SMGs with IFU observations to describe the ionized gas morphologies and kinematics in a sample of 16 SMGs. By comparing the velocity and spatial offsets of the SMG Halpha components with sub-halo offsets in the Millennium simulation database we infer an average halo mass for SMGs of 13<log(M[h^-1Mo])<14. Finally we explore the relationship between the velocity dispersion and star formation intensity within the SMGs, finding the gas motions are consistent with the Kennicutt-Schmidt law and a range of extinction corrections, although might also be driven by the tidal torques from merging or even the star formation itself.

We present the results from a survey for 12CO emission in 40 luminous sub-millimetre galaxies (SMGs), with 850um fluxes of S850 = 4 - 20 mJy, conducted with the Plateau de Bure Interferometer. We detect 12CO emission in 32 SMGs at z~1.2 - 4.1, including 16 SMGs not previously published. Using multiple 12CO line (J_up =2 - 7) observations, we derive a median spectral line energy distribution for luminous SMGs and use this to estimate a mean gas mass of (5.3 +/- 1.0) \times 10^10 Msun. We report the discovery of a fundamental relationship between 12CO FWHM and 12CO line luminosity in high-redshift starbursts, which we interpret as a natural consequence of the baryon-dominated dynamics within the regions probed by our observations. We use far-infrared luminosities to assess the star-formation efficiency in our SMGs, finding a steepening of the L'CO-LFIR relation as a function of increasing 12CO J_up transition. We derive dynamical masses and molecular gas masses, and use these to determine the redshift evolution of the gas content of SMGs, finding that they do not appear to be significantly more gas rich than less vigorously star-forming galaxies at high redshifts. Finally, we collate X-ray observations, and study the interdependence of gas and dynamical properties of SMGs with their AGN activity and supermassive black hole masses (MBH), finding that SMGs lie significantly below the local M_BH-sigma relation. We conclude that SMGs represent a class of massive, gas-rich ultraluminous galaxies with somewhat heterogeneous properties, ranging from starbursting disc-like systems with L~10^12 L_sun, to the most highly star-forming mergers in the Universe.

During our Herschel Lensing Survey (HLS) of massive galaxy clusters, we have discovered an exceptionally bright source behind the z=0.22 cluster Abell 773, which appears to be a strongly lensed submillimeter galaxy (SMG) at z=5.2429. This source is unusual compared to most other lensed sources discovered by Herschel so far, because of its higher submm flux (\sim 200mJy at 500\micron) and its high redshift. The dominant lens is a foreground z=0.63 galaxy, not the cluster itself. The source has a far-infrared (FIR) luminosity of L_FIR= 1.1 10^{14}/\mu Lo, where \mu is the magnification factor, likely \sim 11. We report here the redshift identification through CO lines with the IRAM-30m, and the analysis of the gas excitation, based on CO(7-6), CO(6-5), CO(5-4) detected at IRAM and the CO(2-1) at the EVLA. All lines decompose into a wide and strong red component, and a narrower and weaker blue component, 540\kms apart. Assuming the ultraluminous galaxy (ULIRG) CO-to-H2 conversion ratio, the H2 mass is 5.8 10^{11}/\mu Mo, of which one third is in a cool component. From the CI line we derive a CI/H2 number abundance of 6 10^{-5} similar to that in other ULIRGs. The H2O line is strong only in the red velocity component, with an intensity ratio I(H_2O)/I(CO) \sim 0.5, suggesting a strong local FIR radiation field, possibly from an active nucleus (AGN) component. We detect the [NII]205\mics line for the first time at high-z. It shows comparable blue and red components, with a strikingly broad blue one, suggesting strong ionized gas flows.

[abridged] We present interferometric CO observations of twelve z~2 submillimetre-faint, star-forming radio galaxies (SFRGs) which are thought to be ultraluminous infrared galaxies (ULIRGs) possibly dominated by warmer dust (T_dust ~> 40 K) than submillimetre galaxies (SMGs) of similar luminosities. Four other CO-observed SFRGs are included from the literature, and all observations are taken at the Plateau de Bure Interferometer (PdBI) in the compact configuration. Ten of the sixteen SFRGs observed in CO (63%) are detected at >4sigma with a mean inferred molecular gas mass of ~2*10^10 M_sun. SFRGs trend slightly above the local ULIRG L_FIR-L'_CO relation. Since SFRGs are about two times fainter in radio luminosity but exhibit similar CO luminosities to SMGs, this suggests SFRGs are slightly more efficient star formers than SMGs at the same redshifts. SFRGs also have a narrow mean CO line width, 320+-80km/s. SFRGs bridge the gap between properties of very luminous >5*10^12 L_sun SMGs and those of local ULIRGs and are consistent with intermediate stage major mergers. We suspect that more moderate-luminosity SMGs, not yet surveyed in CO, would show similar molecular gas properties to SFRGs. The AGN fraction of SFRGs is consistent with SMGs and is estimated to be 0.3+-0.1, suggesting that SFRGs are observed near the peak phase of star formation activity and not in a later, post-SMG enhanced AGN phase. This CO survey of SFRGs serves as a pilot project for the much more extensive survey of Herschel and SCUBA-2 selected sources which only partially overlap with SMGs. Better constraints on CO properties of a diverse high-z ULIRG population are needed from ALMA to determine the evolutionary origin of extreme starbursts, and what role ULIRGs serve in catalyzing the formation of massive stellar systems in the early Universe.

We present near-infrared spectroscopic observations from VLT ISAAC of thirteen 250\mu m-luminous galaxies in the CDF-S, seven of which have confirmed redshifts which average to <z > = 2.0 \pm 0.4. Another two sources of the 13 have tentative z > 1 identifications. Eight of the nine redshifts were identified with H{\alpha} detection in H- and K-bands, three of which are confirmed redshifts from previous spectroscopic surveys. We use their near-IR spectra to measure H{\alpha} line widths and luminosities, which average to 415 \pm 20 km/s and 3 \times 10^35 W (implying SFR(H{\alpha})~200 M_\odot /yr), both similar to the H{\alpha} properties of SMGs. Just like SMGs, 250 \mu m-luminous galaxies have large H{\alpha} to far-infrared (FIR) extinction factors such that the H{\alpha} SFRs underestimate the FIR SFRs by ~8-80 times. Far-infrared photometric points from observed 24\mu m through 870\mu m are used to constrain the spectral energy distributions (SEDs) even though uncertainty caused by FIR confusion in the BLAST bands is significant. The population has a mean dust temperature of Td = 52 \pm 6 K, emissivity {\beta} = 1.73 \pm 0.13, and FIR luminosity LFIR = 3 \times 10^13 L_\odot. Although selection at 250\mu m allows for the detection of much hotter dust dominated HyLIRGs than SMG selection (at 850\mu m), we do not find any >60 K 'hot-dust' HyLIRGs. We have shown that near-infrared spectroscopy combined with good photometric redshifts is an efficient way to spectroscopically identify and characterise these rare, extreme systems, hundreds of which are being discovered by the newest generation of IR observatories including the Herschel Space Observatory.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We present high spatial resolution (0.4", ~3.5 kpc) PdBI interferometric data on three ultra-luminous infrared galaxies (ULIRGs) at z~2: two submillimetre galaxies and one submillimetre faint star forming radio galaxy. The three galaxies have been ro- bustly detected in CO rotational transitions, either 12CO(J=4-3) or 12CO(J=3-2), allowing their sizes and gas masses to be accurately constrained. These are the highest spatial resolution observations observed to date (by a factor of ~2) for intermediate-excitation CO emission in z~2 ULIRGs. The galaxies appear extended over several resolution elements, having a mean radius of 3.7 kpc. High-resolution (0.3") combined MERLIN-VLA observations of their radio continua allow an analysis of the star formation behaviour of these galaxies, on comparable spatial scales to that of the CO observations. This 'matched beam' approach sheds light on the spatial distribution of both molecular gas and star formation, and we can therefore calculate accurate star formation rates and gas surface densities: this allows us to place the three systems in the context of a Kennicutt-Schmidt (KS)-style star formation law. We find a difference in size between the CO and radio emission regions, and as such we suggest that using the spatial extent of the CO emission region to estimate the surface density of star formation may lead to error. This size difference also causes the star formation efficiencies within systems to vary by up to a factor of 5. We also find, with our new accurate sizes, that SMGs lie significantly above the KS relation, indicating that stars are formed more efficiently in these extreme systems than in other high-z star forming galaxies.

Using the IRAM Plateau-de-Bure Interferometer (PdBI) we have searched for the upper fine structure line of neutral carbon (CI(2->1), nu_rest = 809 GHz) and 12CO(J=7->6) (nu_rest=806GHz) towards the submillimetre galaxies (SMGs) GN20 (SMM J123711.9+622212, z = 4.055) and GN20.2 (SMM J123708.8+622202, z = 4.051). The far-infrared (FIR) continuum is detected at 8sigma significance in GN20, with a flux density of S_1.8mm = 1.9+-0.2mJy, while no continuum is detected in GN20.2. Both sources are statistically undetected in both CI and CO(7-6) lines; we derive line luminosity limits for both CI and CO of L'<~2x10^10K km s^-1 pc^2. Assuming carbon excitation temperatures of T_ex = 30K (the galaxies' measured dust temperatures), we infer CI mass limits of M_CI < 5.4x10^6 Msun (GN20) and M_CI < 6.8x10^6 Msun (GN20.2). The derived CI abundance limits are <1.8x10^-5 for GN20 and <3.8x10^-5 for GN20.2 implying that the systems have Milky Way level neutral carbon enrichment (X[CI]/X[H2]) or lower, similar to high-redshift carbon-detected systems (at 5x10^-5) but about 50 times less than the neutral carbon enrichment of local starburst galaxies. Observations of GN20 and GN20.2 in high-resolution MERLIN+VLA radio maps of GOODS-N are used to further constrain the sizes and locations of active regions. We conclude that the physical gas properties of young rapidly evolving systems like GN20 and GN20.2 are likely significantly different than starburst/ULIRG environments in the local Universe yet similar to z~2 SMGs. Unless gravitationally amplified examples can be found, observations of galaxies like GN20 will require the order of magnitude increase in sensitivity of the Atacama Large Millimetre Array (ALMA) to constrain their CI and high-J CO content, despite the fact that they are the brightest systems at z~4.

We identify eight z>1 radio sources undetected at 850um but robustly detected at 70um, confirming that they represent ultraluminous infrared galaxies (ULIRGs) with hotter dust temperatures (T_dust=52+-10 K) than Submillimetre Galaxies (SMGs) at similar luminosities and redshifts. These galaxies share many properties with SMGs: ultra-violet (UV) spectra consistent with starbursts, high stellar masses and radio luminosities. We can attribute their radio emission to star formation since high-resolution MERLIN radio maps show extended emission regions (with half light radii of 2-3kpc), which are unlikely to be generated by AGN activity. These observations provide the first direct confirmation of hot, dusty ULIRGs which are missed by current submillimetre surveys. They have significant implications for future observations from the Herschel Space Observatory and SCUBA2, which will select high redshift luminous galaxies with less selection biases.

We present confusion limited submillimetre (submm) observations with the SCUBA camera on the JCMT of the z=2.83 Lyman-break galaxy (LBG), Westphal-MM8, reaching an 850um sensitivity even greater than that achieved in the SCUBA map of the Hubble Deep Field region. The detection of Westphal-MM8 (S_{850um} = 1.98+-0.48mJy), along with the literature submm detections of lensed LBGs, suggest that the LBG population may contribute significantly to the source counts of submm selected galaxies in the 1-2mJy regime. Additionally, submm-luminous LBGs are a viable progenitor population for the recently discovered evolved galaxies at z~2-2.5. These observations represent an important baseline for SCUBA2 observations which will regularly map large regions of the sky to this depth.

One of the six extended X-ray sources found in the Chandra DeepField North is centred on HDF130, which has recently been shown to be a massive galaxy at z=1.99 with a compact radio nucleus. The X-ray source has a roughly double-lobed structure with each lobe about 41 arcsec long, or 345 kpc at the redshift of HDF130. We have analyzed the 2 Ms X-ray image and spectrum of the source and find that it is well fit by a power-law continuum of photon index 2.65 and has a 2--10 keV luminosity of 5.4x10^{43}ergps (if at z=1.99). Any further extended emission within a radius of 60 arcsec has a luminosity less than half this value, which is contrary to what is expected from a cluster of galaxies. The source is best explained as an inverse Compton ghost of a giant radio source, which is no longer being powered, and for which Compton losses have downgraded the energetic electrons, \gamma> 10^4, required for high-frequency radio emission. The lower energy electrons, \gamma~1000, produce X-rays by inverse Compton scattering on the Cosmic Microwave Background. Depending on the magnetic field strength, some low frequency radio emission may remain. Further inverse Compton ghosts may exist in the Chandra deep fields.

We present high spatial resolution MERLIN 1.4GHz radio observations of two high redshift (z~2) sources, RGJ123623 (HDF147) and RGJ123617 (HDF130), selected as the brightest radio sources from a sample of submillimetre-faint radio galaxies. They have starburst classifications from their rest-frame UV spectra. However, their radio morphologies are remarkably compact (<80mas and <65mas respectively), demanding that the radio luminosity be dominated by Active Galactic Nuclei (AGN) rather than starbursts. Near-IR imaging (HST NICMOS F160W) shows large scale sizes (R_(1/2)~0.75", diameters ~12kpc) and SED fitting to photometric points (optical through the mid-IR) reveals massive (~5x10^(11) M_sun), old (a few Gyr) stellar populations. Both sources have low flux densities at observed 24um and are undetected in observed 70um and 850um, suggesting a low mass of interstellar dust. They are also formally undetected in the ultra-deep 2Ms Chandra data, suggesting that any AGN activity is likely intrinsically weak. We suggest both galaxies have evolved stellar populations, low star formation rates, and low accretion rates onto massive black holes (10^(8.6) M_sun) whose radio luminosity is weakly beamed (by factors of a few). A cluster-like environment has been identified near HDF130 by an over-density of galaxies at z=1.99, reinforcing the claim that clusters lead to more rapid evolution in galaxy populations. These observations suggest that high-resolution radio (MERLIN) can be a superb diagnostic tool of AGN in the diverse galaxy populations at z~2.

This submission has been withdrawn.

This is the first measurement and detection of coherence in the intergalactic medium (IGM) at substantially high redshift (z~3.8) and on large physical scales (~2.5 h^-1 Mpc). We perform the measurement by presenting new observations from Keck LRIS of the high redshift quasar pair PC 1643+4631A, B and their Ly-alpha absorber coincidences. This experiment extends multiple sightline quasar absorber studies to higher redshift, higher opacity, larger transverse separation, and into a regime where coherence across the IGM becomes weak and difficult to detect. We fit 222 discrete Ly-alpha absorbers to sightline A and 211 to sightline B. Relative to a Monte Carlo pairing test (using symmetric, nearest neighbor matching) the data exhibit a 4sigma excess of pairs at low velocity splitting (<150 km/s), thus detecting coherence on transverse scales of ~2.5 h^-1 Mpc. We use spectra extracted from an SPH simulation to analyze symmetric pair matching, transmission distributions as a function of redshift and compute zero-lag cross-correlations to compare with the quasar pair data. The simulations agree with the data with the same strength (~4sigma) at similarly low velocity splitting above random chance pairings. In cross-correlation tests, the simulations agree when the mean flux (as a function of redshift) is assumed to follow the prescription given by Kirkman et al. (2005). While the detection of flux correlation (measured through coincident absorbers and cross-correlation amplitude) is only marginally significant, the agreement between data and simulations is encouraging for future work in which even better quality data will provide the best insight into the overarching structure of the IGM and its understanding as shown by SPH simulations.

We outline a strategy to select faint (i<24.5) type 1 AGN candidates down to the Seyfert/QSO boundary for spectroscopic targeting in the COSMOS field, picking candidates by their nonstellar colors in broadband ground-based photometry and morphological properties extracted from HST-ACS. AGN optical color selection has not been applied to such faint magnitudes in such a large continuous part of the sky. Hot stars are known to be the dominant contaminant for bright AGN candidate selection at z<2, but we anticipate the highest color contamination at all redshifts to be from faint starburst and compact galaxies. Morphological selection via the Gini Coefficient separates most potential AGN from these faint blue galaxies. Recent models of the quasar luminosity function are used to estimate quasar surface densities, and studies of stellar populations in the COSMOS field infer stellar contamination. We use 292 spectroscopically confirmed type 1 AGN and quasar templates to predict AGN colors with redshift, and contrast those predictions with the colors of known contaminating populations. The motivation of this study and subsequent spectroscopic follow-up is to populate and refine the faint end of the QLF where the population of type 1 AGN is presently not well known. The anticipated AGN observations will add to the ~300 already known AGN in the COSMOS field, making COSMOS a densely packed field of quasars to be used to understand supermassive black holes and probe the structure of the intergalactic medium in the intervening volume.