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.

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

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.