We present Gemini-S and {\it Spitzer}-IRAC optical-through-near-IR observations in the field of the SPT2349-56 proto-cluster at $z=4.3$. We detect optical/IR counterparts for only nine of the 14 submillimetre galaxies (SMGs) previously identified by ALMA in the core of SPT2349-56. In addition, we detect four $z\sim4$ Lyman-break galaxies (LBGs) in the 30 arcsec diameter region surrounding this proto-cluster core. Three of the four LBGs are new systems, while one appears to be a counterpart of one of the nine observed SMGs. We identify a candidate brightest cluster galaxy (BCG) with a stellar mass of $(3.2^{+2.5}_{-1.4})\times10^{11}\,{\rm M}_{\odot}$. The stellar masses of the eight other SMGs place them on, above, and below the main sequence of star formation at $z\approx4.5$. The cumulative stellar mass for the SPT2349-56 core is at least $(11.5\pm2.9)\times10^{11}\,{\rm M}_{\odot}$, a sizeable fraction of the stellar mass in local BCGs, and close to the universal baryon fraction (0.16) relative to the virial mass of the core ($10^{13}\,{\rm M}_{\odot}$). As all 14 of these SMGs are destined to quickly merge, we conclude that the proto-cluster core has already developed a significant stellar mass at this early stage, comparable to $z=1$ BCGs. Importantly, we also find that the SPT2349-56 core structure would be difficult to uncover in optical surveys, with none of the ALMA sources being easily identifiable or constrained through $g,r,$ and $i$ colour-selection in deep optical surveys and only a modest overdensity of LBGs over the extended core structure. SPT2349-56 therefore represents a truly dust-obscured phase of a massive cluster core under formation.

Matching the number counts of high-$z$ sub-millimetre-selected galaxies (SMGs) has been a long standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the sub-mm emission from galaxies in the SIMBA cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850$\mathrm{\mu m}$-selected sources. We find unprecedented agreement with the integrated 850$\mathrm{\mu m}$ luminosity function, along with good agreement in the redshift distribution of bright SMGs. The excellent agreement is driven primarily by SIMBA's good match to infrared measures of the star formation rate (SFR) function between $z = 2-4$ at high SFRs. Also important is the self-consistent on-the-fly dust model in SIMBA, which predicts, on average, higher dust masses (by up to a factor of 7) compared to using a fixed dust-to-metals ratio of 0.3. We construct a lightcone to investigate the effect of far-field blending, and find minimal contribution to the shape and normalisation of the luminosity function. We provide new fits to the 850$\mathrm{\mu m}$ luminosity as a function of SFR and dust mass. Our results demonstrate that exotic solutions to the discrepancy between sub-mm counts in simulations and observations, such as a top-heavy IMF, are unnecessary, and that sub-millimetre-bright phases are a natural consequence of massive galaxy evolution.

We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, Simba, EAGLE, and Illustris-TNG, versus observations from $z\sim 0-2$. These simulations all rely on similar sub-resolution prescriptions to model cold interstellar gas which they cannot represent directly, and qualitatively reproduce the observed $z\approx 0$ HI and H$_2$ mass functions (HIMF, H2MF), CO(1-0) luminosity functions (COLF), and gas scaling relations versus stellar mass, specific star formation rate, and stellar surface density $\mu_*$, with some quantitative differences. To compare to the COLF, we apply an H$_2$-to-CO conversion factor to the simulated galaxies based on their average molecular surface density and metallicity, yielding substantial variations in $\alpha_{\rm CO}$ and significant differences between models. Using this, predicted $z=0$ COLFs agree better with data than predicted H2MFs. Out to $z\sim 2$, EAGLE's and Simba's HIMF and COLF strongly increase, while TNG's HIMF declines and COLF evolves slowly. EAGLE and Simba reproduce high $L_{\rm CO1-0}$ galaxies at $z\sim 1-2$ as observed, owing partly to a median $\alpha_{\rm CO}(z=2)\sim 1$ versus $\alpha_{\rm CO}(z=0)\sim 3$. Examining \HI, H$_2$, and CO scaling relations, their trends with $M_*$ are broadly reproduced in all models, but EAGLE yields too little HI in green valley galaxies, TNG and Simba overproduce cold gas in massive galaxies, and Simba overproduces molecular gas in small systems. Using Simba variants that exclude individual AGN feedback modules, we find that Simba's AGN jet feedback is primarily responsible by lowering cold gas contents from $z\sim 1\to0$ by suppressing cold gas in $M_*> 10^{10}{\rm M}_\odot$ galaxies, while X-ray feedback suppresses the formation of high-$\mu_*$ systems.

The South Pole Telescope (SPT) has systematically identified 81 high-redshift, strongly gravitationally lensed, dusty star-forming galaxies (DSFGs) in a 2500 square degree cosmological mm-wave survey. We present the final spectroscopic redshift survey of this flux-limited ($S_{870 \, \mathrm{\mu m}} > 25 \, \mathrm{mJy}$) sample. The redshift survey was conducted with the Atacama Large Millimeter/submillimeter Array across the 3 mm spectral window, targeting carbon monoxide line emission. The SPT sample is now spectroscopically complete, with redshifts spanning the range $1.9$$<$$z$$<$$6.9$, with a median of $z=3.9 \pm 0.2$. We present the mm through far-infrared photometry and spectral energy density fits for all sources, along with their inferred intrinsic properties. Comparing the properties of the SPT sources to the unlensed DSFG population, we demonstrate that the SPT-selected DSFGs represent the most extreme infrared-luminous galaxies, even after accounting for strong gravitational lensing. The SPT sources have a median star formation rate of $2.3(2)\times 10^3 \, \mathrm{M_\odot yr^{-1}}$ and a median dust mass of $1.4(1)\times10^9\, \mathrm{M_\odot}$. However, the inferred gas depletion timescales of the SPT sources are comparable to those of unlensed DSFGs, once redshift is taken into account. This SPT sample contains roughly half of the known spectroscopically confirmed DSFGs at $z$$>$$5$, making this the largest sample of high-redshift DSFGs to-date, and enables us to measure the "high-redshift tail" of the distribution of luminous DSFGs. Though galaxy formation models struggle to account for the SPT redshift distribution, the larger sample statistics from this complete and well-defined survey will help inform future theoretical efforts.

We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.

The primary method for inferring the stellar mass ($M_*$) of a galaxy is through spectral energy distribution (SED) modeling. However, the technique rests on assumptions such as the galaxy star formation history and dust attenuation law that can severely impact the accuracy of derived physical properties from SED modeling. Here, we examine the effect that the assumed star formation history (SFH) has on the stellar properties inferred from SED fitting by ground truthing them against mock observations of high-resolution cosmological hydrodynamic galaxy formation simulations. Classically, SFHs are modeled with simplified parameterized functional forms, but these forms are unlikely to capture the true diversity of galaxy SFHs and may impose systematic biases with under-reported uncertainties on results. We demonstrate that flexible nonparametric star formation histories outperform traditional parametric forms in capturing variations in galaxy star formation histories, and as a result, lead to significantly improved stellar masses in SED fitting. We find a decrease in the average bias of 0.4 dex with a delayed-$\tau$ model to a bias of just under 0.05 dex for the nonparametric model. Similarly, using nonparametric star formation histories in SED fitting result in increased accuracy in recovered galaxy star formation rates (SFRs) and stellar ages.

The [CII] 158 micron line is one of the strongest IR emission lines, which has been shown to trace the star-formation rate (SFR) of galaxies in the nearby Universe and up to $z \sim 2$. Whether this is also the case at higher redshift and in the early Universe remains debated. The ALPINE survey, which targeted 118 star-forming galaxies at $4.4 < z< 5.9$, provides a new opportunity to examine this question with the first statistical dataset. Using the ALPINE data and earlier measurements from the literature we examine the relation between the [CII] luminosity and the SFR over the entire redshift range from $z \sim 4-8$. ALPINE galaxies, which are both detected in [CII] and dust continuum, show a good agreement with the local L([CII])-SFR relation. Galaxies undetected in the continuum with ALMA are found to be over-luminous in [CII], when the UV SFR is used. After accounting for dust-obscured star formation, by an amount SFR(IR)$\approx$SFR(UV) on average, which results from two different stacking methods and SED fitting, the ALPINE galaxies show an L([CII])-SFR relation comparable to the local one. When [CII] non-detections are taken into account, the slope may be marginally steeper at high-z, although this is still somewhat uncertain. When compared in a homogeneous manner, the $z>6 $ [CII] measurements (detections and upper limits) do not behave very differently from the $z \sim 4-6$ data. We find a weak dependence of L([CII])/SFR on the Lyman-alpha equivalent width. Finally, we find that the ratio L([CII])/LIR $\sim (1-3) \times 10^{-3}$ for the ALPINE sources, comparable to that of "normal" galaxies at lower redshift. Our analysis, which includes the largest sample ($\sim 150$ galaxies) of [CII] measurements at $z>4$ available so far, suggests no or little evolution of the L([CII])-SFR relation over the last 13 Gyr of cosmic time.

We study specific star formation rate (sSFR) and gas profiles of star forming and green valley galaxies in the Simba cosmological hydrodynamic simulation. Star-forming galaxy half-light radii ($R_{\rm half}$) at $z=0$ and their evolution ($\propto(1+z)^{-0.78}$) agree with observations. Passive galaxy $R_{\rm half}$ agree with observations at high redshift, but by $z=0$ are too large, owing to numerical heating. We compare Simba $z=0$ sSFR radial profiles for star forming and green valley galaxies to observations. Simba shows strong central depressions in star formation rate (SFR), sSFR, and gas fraction in green valley galaxies and massive star-forming systems, qualitatively as observed, owing to black hole X-ray feedback, which pushes central gas outwards. Turning off X-ray feedback leads to centrally peaked sSFR profiles as in other simulations. In conflict with observations, Simba yields green valley galaxies with strongly dropping sSFR profiles beyond $\gtrsim R_{\rm half}$, regardless of AGN feedback. The central depression owes to lowering molecular gas content; the drop in the outskirts owes to reduced star formation efficiency. Simba's satellites have higher central sSFR and lower outskirts sSFR than centrals, in qualitative agreement with observations. At $z=2$ Simba does not show central depressions in massive star-forming galaxies, suggesting Simba's X-ray feedback should be more active at high-$z$. High resolution tests indicate central sSFR suppression is not sensitive to numerical resolution. Reproducing the central sSFR depression in $z=0$ green valley galaxies represents a unique success of Simba. The remaining discrepancies highlight the importance of SFR and gas profiles in constraining quenching mechanisms.

We present the physical extent of [CII] 158um line-emitting gas from 46 star-forming galaxies at z=4-6 from the ALMA Large Program to INvestigate CII at Early Times (ALPINE). Using exponential profile fits, we measure the effective radius of the [CII] line (r_e,[CII]) for individual galaxies and compare them with the rest-frame ultra-violet (UV) continuum (r_e,UV) from Hubble Space Telescope images. The effective radius r_e,[CII] exceeds r_e,UV by factors of ~2-3 and the ratio of r_e,[CII]/r_e,UV increases as a function of M_star. We do not find strong evidence that [CII] line, the rest-frame UV, and FIR continuum are always displaced over ~ 1-kpc scale from each other. We identify 30% of isolated ALPINE sources as having an extended [CII] component over 10-kpc scales detected at 4.1$\sigma$-10.9$\sigma$ beyond the size of rest-frame UV and far-infrared (FIR) continuum. One object has tentative rotating features up to ~10-kpc, where the 3D model fit shows the rotating [CII]-gas disk spread over 4 times larger than the rest-frame UV-emitting region. Galaxies with the extended [CII] line structure have high star-formation rate (SFR), stellar mass (M_star), low Lya equivalent-width, and more blue-shifted (red-shifted) rest-frame UV metal absorption (Lya line), as compared to galaxies without such extended [CII] structures. Although we cannot rule out the possibility that a selection bias towards luminous objects may be responsible for such trends, the star-formation driven outflow also explains all these trends. Deeper observations are essential to test whether the extended [CII] line structures are ubiquitous to high-z star-forming galaxies.

We present the first predictions for the $L_{\rm [CII]}$ - SFR relation and [CII] luminosity function (LF) in the Epoch of Reionization (EoR) based on cosmological hydrodynamics simulations using the SIMBA suite plus radiative transfer calculations via S\'IGAME. The sample consists of 11,137 galaxies covering halo mass $\log M_{\rm halo}\in$[9, 12.4] $M_\odot$, star formation rate SFR$\in$[0.01, 330] $M_\odot$ yr$^{-1}$, and metallicity $<Z_{\rm gas}>_{\rm SFR}\in$[0.1, 0.9] $Z_\odot$. The simulated $L_{\rm [CII]}$-SFR relation is consistent with the range observed, but with a spread of $\simeq$0.3 dex at the high end of SFR ($>$100 $M_\odot$ yr$^{-1}$) and $\simeq$0.6 dex at the lower end, and there is tension between our predictions and the values of $L_{\rm [CII]}$ above 10$^{8.5}$ $L_\odot$ observed in some galaxies reported in the literature. The scatter in the $L_{\rm [CII]}$-SFR relation is mostly driven by galaxy properties, such that at a given SFR, galaxies with higher molecular gas mass and metallicity have higher $L_{\rm [CII]}$. The [CII] LF predicted by SIMBA is consistent with the upper limits placed by the only existing untargeted flux-limited [CII] survey at the EoR (ASPECS) and those predicted by semi-analytic models. We compare our results with existing models and discuss differences responsible for the discrepant slopes in the $L_{\rm [CII]}$-SFR relatiion.

We present dust attenuation properties of spectroscopically confirmed star forming galaxies on the main sequence at redshift ~4.4-5.8. Our analyses are based on the far infrared continuum observations of 118 galaxies at rest-frame $158\,\rm{\mu m}$ obtained with the ALMA large program ALPINE. We study the connection between the UV spectral slope ($\beta$), stellar mass ($M_{\star}$), and infrared excess (IRX$=L_{\rm{IR}}/L_{\rm{UV}}$). Twenty-three galaxies are individually detected in the continuum at $>3.5\,\sigma$ significance. We perform a stacking analysis using both detections and non-detections to study the average dust attenuation properties at z~4.4-5.8. The individual detections and stacks show that the IRX-$\beta$ relation at z~5 is consistent with a steeper dust attenuation curve than typically found at lower redshifts (z<4). The attenuation curve is similar to or even steeper than that of the extinction curve of the Small Magellanic Cloud (SMC). This systematic change of the IRX-$\beta$ relation as a function of redshift suggests an evolution of dust attenuation properties at $z>4$. Similarly, we find that our galaxies have lower IRX values up to 1 dex on average at fixed mass compared to previously studied IRX-$M_{\star}$ relations at $z\lesssim4$, albeit with significant scatter. This implies a lower obscured fraction of star-formation than at lower redshifts. Our results suggest that dust properties of UV-selected star forming galaxies at $z\gtrsim4$ are characterised by (i) a steeper attenuation curve than at $z\lesssim4$, and (ii) a rapidly decreasing dust obscured fraction of star formation as a function of redshift. Nevertheless, even among this UV-selected sample, massive galaxies ($\log M_{\star}/M_\odot > 10$) at z~5-6 already exhibit an obscured fraction of star formation of $\sim45\%$, indicating a rapid build-up of dust during the epoch of reionization.

The molecular gas content of normal galaxies at z>4 is poorly constrained, because the commonly used molecular gas tracers become hard to detect. We use the [CII]158um luminosity, recently proposed as a molecular gas tracer, to estimate the molecular gas content in a large sample of main-sequence star-forming galaxies at z=4.4-5.9, with a median stellar mass of 10^9.7 Msun, drawn from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. The agreement between molecular gas masses derived from [CII] luminosity, dynamical mass, and rest-frame 850um luminosity, extrapolated from the rest-frame 158um continuum, supports [CII] as a reliable tracer of molecular gas in our sample. We find a continuous decline of the molecular gas depletion timescale from z=0 to z=5.9, which reaches a mean value of (4.6+/-0.8)x10^8 yr at z~5.5, only a factor of 2-3 shorter than in present-day galaxies. This suggests a mild enhancement of star formation efficiency toward high redshifts, unless the molecular gas fraction significantly increases. Our estimates show that the rise in molecular gas fraction as reported previously, flattens off above z~3.7 to achieve a mean value of 63%+/-3 over z=4.4-5.9. This redshift evolution of the gas fraction is in line with the one of the specific star formation rate. We use multi-epoch abundance matching to follow the gas fraction evolution over cosmic time of progenitors of z=0 Milky Way-like galaxies in 10^13 Msun halos and of more massive z=0 galaxies in 10^14 Msun halos. Interestingly, the former progenitors show a monotonic decrease of the gas fraction with cosmic time, while the latter show a constant gas fraction from z=5.9 to z~2 and a decrease at z<2. We discuss three possible effects, namely outflows, halt of gas supplying, and over-efficient star formation, which may jointly contribute to the gas fraction plateau of the latter massive galaxies.

High-redshift Lyman-alpha blobs (LABs) are an enigmatic class of objects that have been the subject of numerous observational and theoretical investigations. It is of particular interest to determine the dominant power sources for the copious luminosity, as direct emission from HII regions, cooling gas, and fluorescence due to the presence of active galactic nuclei (AGN) can all contribute significantly. In this paper, we present the first theoretical model to consider all of these physical processes in an attempt to develop an evolutionary model for the origin of high-z LABs. This is achieved by combining a series of high-resolution cosmological zoom-in simulations with ionization and Lyman-alpha (Lya) radiative transfer models. We find that massive galaxies display a range of Lya luminosities and spatial extents (which strongly depend on the limiting surface brightness used) over the course of their lives, though regularly exhibit luminosities and sizes consistent with observed LABs. The model LABs are typically powered from a combination of recombination in star-forming galaxies, as well as cooling emission from gas associated with accretion. When AGN are included in the model, the fluorescence caused by AGN-driven ionization can be a significant contributor to the total Lya luminosity as well. We propose that the presence of an AGN may be predicted from the Gini coefficient of the blob's surface brightness. Within our modeled mass range, there are no obvious threshold physical properties that predict appearance of LABs, and only weak correlations of the luminosity with the physical properties of the host galaxy. This is because the emergent Lya luminosity from a system is a complex function of the gas temperature, ionization state, and Lya escape fraction.

Understanding the properties and physical mechanisms that shape dust attenuation curves in galaxies is one of the fundamental questions of extragalactic astrophysics, with a great practical significance for deriving the physical properties of galaxies, such as the star formation rate and stellar mass. Attenuation curves result from a combination of dust grain properties, dust content, and the spatial arrangement of dust and different populations of stars. In this review we assess the current state of the field, paying particular attention to the importance of extinction curves as the building blocks of attenuation laws. We introduce a quantitative framework to characterize and compare extinction and attenuation curves, present a theoretical foundation for interpreting empirical results, overview an array of observational methods, and review the observational state of the field at both low and high redshift. Our main conclusions are: Attenuation curves exhibit a large range of slopes, from curves with shallow (Milky Way-like) slopes to those exceeding the slope of the SMC extinction curve. The slopes of the curves correlate strongly with the effective optical opacities, in the sense that galaxies with low dust column density (lower visual attenuation) tend to have steeper slopes, whereas the galaxies with high dust column density have shallower (grey) slopes. Galaxies appear to exhibit a diverse range of 2175A UV bump strengths, but on average have suppressed bumps compared to the average Milky Way extinction curve. Theoretical studies indicate that variations in bump strength may result from similar geometric and radiative transfer effects that drive the correlation between the slope and the dust column.

We report ALMA observations of the neutral atomic carbon transitions [CI] and multiple CO lines in a sample of $\sim30$ main sequence galaxies at $z\sim1$, including novel information on [CI](2-1) and CO(7-6) for 7 of such normal objects. We complement our observations with a collection of $>200$ galaxies with coverage of similar transitions, spanning the $z=0-4$ redshift interval and a variety of ambient conditions from local to high-redshift starbursts. We find systematic variations in the [CI]/IR and [CI]/high-$J$ ($J=7$) CO luminosity ratios among the various samples. We interpret these differences as increased dense molecular gas fractions and star formation efficiencies in the strongest high-redshift starbursts with respect to normal main sequence galaxies. We further report constant $L'_{\rm [CI]2-1}$/$L'_{\rm [CI]1-0}$ ratios across the galaxy populations and redshifts, suggesting that gas temperatures $T_{\rm exc}$ traced by [CI] do not strongly vary. We find only a mild correlation with $T_{\rm dust}$ and that, generally, $T_{\rm exc} \lesssim T_{\rm dust}$. We fit the line ratios with classical PDR models, retrieving consistently larger densities and intensities of the UV radiation fields in submm galaxies than in main sequence and local objects. However, these simple models fall short in representing the complexity of a multiphase interstellar medium and should be treated with caution. Finally, we compare our observations with the Santa Cruz semi-analytical model of galaxy evolution, recently extended to simulate submm emission. While we confirm the success in reproducing the CO lines, we find systematically larger [CI] luminosities at fixed IR luminosity than predicted theoretically. This highlights the necessity of improving our understanding of the mechanisms regulating the [CI] emission on galactic scales. We release our data compilation to the community.

We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO $(8-7)$, $(9-8)$, $\rm H_{2}O (2_{0,2}-1_{1,1})$ and $\rm OH^{+} (1_{1}-0_{1})$ and NOrthern Extended Millimeter Array (NOEMA) observations of CO $(5-4)$, $(6-5)$, $(12-11)$ and $(13-12)$ towards the $z = 6.003$ quasar SDSS J231038.88+185519.7, aiming to probe the physical conditions of the molecular gas content of this source. We present the best sampled CO spectral line energy distribution (SLED) at $z = 6.003$, and analyzed it with the radiative transfer code MOLPOP-CEP. Fitting the CO SLED to a one-component model indicates a kinetic temperature $T_{\rm kin} = 228 \ \rm K$, molecular gas density $log (n(\rm H_{2})/\rm cm^{-3}$ )=4.75, and CO column density $log(N(\rm CO)/\rm cm^{-2}) =17.5$, although a two-component model better fits the data. In either case, the CO SLED is dominated by a "warm" and "dense" component. Compared to samples of local (Ultra) Luminous Infrared Galaxies ((U)LIRGs), starburst galaxies and high redshift Submillimeter Galaxies (SMGs), J2310+1855 exhibits higher CO excitation at ($J \geq 8$), like other high redshift quasars. The high CO excitation, together with the enhanced $L_{\rm H_{2}O}/ L_{IR} $, $L_{\rm H_{2}O}/ L_{CO} $ and $L_{OH^{+}}/L_{\rm H_{2}O} $ ratios, suggests that besides the UV radiation from young massive stars, other mechanisms such as shocks, cosmic rays and X-rays might also be responsible for the heating and ionization of the molecular gas. In the nuclear region probed by the molecular emissions lines, any of these mechanisms might be present due to the powerful quasar and the starburst activity.

The Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report.

We study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 um emission in a large sample of normal star-forming galaxies at 4 < z < 6 from the ALMA Large Program to INvestigate [C II] at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked [C II] profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked [C II] spectrum, with velocities of |v|<~ 500 km/s. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFRmed = 25 Msun/yr), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at z > 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v|< 200 km/s) of the higher-SFR galaxies is extended on physical sizes of ~ 30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.

We present ALMA observations of the [CII] 158 $\mu$m fine structure line and dust continuum emission from two quasars, SDSS J104433.04-012502.2 and SDSS J012958.51-003539.7, at z=5.78. The ALMA observations at 0.2'' resolution map the dust and gas on kpc scales. The spatially resolved emission shows a similar trend of decreasing [CII]-FIR ratios with increasing FIR surface brightnesses as was found in the infrared luminous galaxies with intense star formation. We confirm the velocity gradients of [CII] emission found previously in SDSS J0129-0035. No clear evidence of order motion is detected in SDSS J1044-0125. The velocity maps and PV-diagrams also suggest turbulent gas clumps in both objects. We tentatively detect a [CII] peak offset 4.9 kpc to the East of SDSS J1044-0125. This may be associated with an infalling companion, or node of gas outflow. All these results suggest significant dynamical evolution of the ISM in the nuclear region of these young quasar-starburst systems. We fit the velocity map of the [CII] emission from SDSS J0129-0035 with a rotating disk model. The result suggests a face-on system with an inclination angle of 16 +/- 20 degree and constrains the lower limit of the host galaxy dynamical mass to be 2.6x10^10 Msun within the [CII] emitting region. It is likely that SDSS J0129-0035, as well as other young quasars with super massive black hole masses on order of 10^7 to 10^8 Msun, falls close to the black hole and host galaxy mass relation defined by local galaxies.

We report the serendipitous detection of two 3 mm continuum sources found in deep ALMA Band 3 observations to study intermediate redshift galaxies in the COSMOS field. One is near a foreground galaxy at 1.3", but is a previously unknown dust-obscured star-forming galaxy (DSFG) at probable $z_{CO}=3.329$, illustrating the risk of misidentifying shorter wavelength counterparts. The optical-to-mm spectral energy distribution (SED) favors a grey $\lambda^{-0.4}$ attenuation curve and results in significantly larger stellar mass and SFR compared to a Calzetti starburst law, suggesting caution when relating progenitors and descendants based on these quantities. The other source is missing from all previous optical/near-infrared/sub-mm/radio catalogs ("ALMA-only"), and remains undetected even in stacked ultradeep optical ($>29.6$ AB) and near-infrared ($>27.9$ AB) images. Using the ALMA position as a prior reveals faint $SNR\sim3$ measurements in stacked IRAC 3.6+4.5, ultradeep SCUBA2 850$\mu$m, and VLA 3GHz, indicating the source is real. The SED is robustly reproduced by a massive $M^*=10^{10.8}$M$_\odot$ and $M_{gas}=10^{11}$M$_\odot$, highly obscured $A_V\sim4$, star forming $SFR\sim300$ M$_{\odot}$yr$^{-1}$ galaxy at redshift $z=5.5\pm$1.1. The ultrasmall 8 arcmin$^{2}$ survey area implies a large yet uncertain contribution to the cosmic star formation rate density CSFRD(z=5) $\sim0.9\times10^{-2}$ M$_{\odot}$ yr$^{-1}$ Mpc$^{-3}$, comparable to all ultraviolet-selected galaxies combined. These results indicate the existence of a prominent population of DSFGs at $z>4$, below the typical detection limit of bright galaxies found in single-dish sub-mm surveys, but with larger space densities $\sim3 \times 10^{-5}$ Mpc$^{-3}$, higher duty cycles $50-100\%$, contributing more to the CSFRD, and potentially dominating the high-mass galaxy stellar mass function.

We present predictions for the evolution of the galaxy dust-to-gas (DGR) and dust-to-metal (DTM) ratios from z=0 to 6, using a model for the production, growth, and destruction of dust grains implemented into the \simba\ cosmological hydrodynamic galaxy formation simulation. In our model, dust forms in stellar ejecta, grows by the accretion of metals, and is destroyed by thermal sputtering and supernovae. Our simulation reproduces the observed dust mass function at z=0, but modestly under-predicts the mass function by ~x3 at z ~ 1-2. The z=0 DGR vs metallicity relationship shows a tight positive correlation for star-forming galaxies, while it is uncorrelated for quenched systems. There is little evolution in the DGR-metallicity relationship between z=0-6. We use machine learning techniques to search for the galaxy physical properties that best correlate with the DGR and DTM. We find that the DGR is primarily correlated with the gas-phase metallicity, though correlations with the depletion timescale, stellar mass and gas fraction are non-negligible. We provide a crude fitting relationship for DGR and DTM vs. the gas-phase metallicity, along with a public code package that estimates the DGR and DTM given a set of galaxy physical properties.

We use the Simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\leq$ 4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass $M_*$ well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate sSFR$>t_{H}^{-1}$ to sSFR$<0.2t_{H}^{-1}$, where $t_{H}$ is the Hubble time. At $z\approx 0-3$, mergers show $\sim\times 2-3$ higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $\leq 1\%$ of the cosmic SF budget. At low masses, the increase in SFR in mergers is mostly attributed to an increase in the $H_2$ content, but for $M_*\geq 10^{10.5} M_{\odot}$ mergers also show an elevated star formation efficiency suggesting denser gas within merging galaxies. The merger rate for star-forming galaxies shows a rapid increase with redshift $\propto (1+z)^{3.5}$, but the quenching rate evolves much more slowly, $\propto (1+z)^{0.9}$; there are insufficient mergers to explain the quenching rate at $z\leq 1.5$. Simba first quenches galaxies at $z\geq 3$, with a number density in good agreement with observations. The quenching timescales $\tau_q$ are strongly bimodal, with `slow' quenchings ($\tau_q \sim 0.1t_{H}$) dominating overall, but `fast' quenchings ($\tau_q\sim 0.01 t_H$) dominating in $M_*\sim 10^{10}-10^{10.5}M_{\odot}$ galaxies, likely induced by Simba's jet-mode black hole feedback. The delay time distribution between mergers and quenching events suggests no physical connection to either fast or slow quenching. Hence, Simba predicts that major mergers induce starbursts, but are unrelated to quenching in either fast or slow mode.

We present high-fidelity, 30 milliarcsecond (200-pc) resolution ALMA rest-frame 240 $\mu$m observations of cold dust emission in three typical main-sequence star-forming galaxies (SFGs) at $z \sim 3$ in the Hubble Ultra-Deep Field (HUDF). The cold dust is distributed within the smooth disk-like central regions of star formation $1 - 3$ kpc in diameter, despite their complex and disturbed rest-frame UV and optical morphologies. No dust substructures or clumps are seen down to $\simeq 1- 3$ $M_\odot$yr$^{-1}$ (1$\sigma$) per 200-pc beam. No dust emission is observed at the locations of UV-emitting clumps, which lie $\simeq 2-10$ kpc from the bulk of star formation. Clumpy substructures can contribute no more than $1-7$% of the total star formation in these galaxies (3$\sigma$ upper limits). The lack of star-forming substructures in our HUDF galaxies is to be contrasted with the multiple substructures characteristic of submillimeter-selected galaxies (SMGs) at the same cosmic epoch, particularly the far-IR-bright SMGs with similarly high-fidelity ALMA observations of Hodge et al. (2019). Individual star-forming substructures in these SMGs contain $\sim10-30$% of their total star formation. A substructure in these SMGs is often comparably bright in the far-infrared as (or in some cases brighter than) our typical SFGs, suggesting that these SMGs originate from a class of disruptive event involving multiple objects at the scale of our HUDF galaxies. The scale of the disruptive event found in our main-sequence SFGs, characterized by the lack of star-forming substructures at our resolution and sensitivity, could be less violent, e.g., gas-rich disk instability or minor mergers.

Water ($\rm H_{2}O$), one of the most ubiquitous molecules in the universe, has bright millimeter-wave emission lines easily observed at high-redshift with the current generation of instruments. The low excitation transition of $\rm H_{2}O$, p$-$$\rm H_{2}O$(202 $-$ 111) ($\nu_{rest}$ = 987.927 GHz) is known to trace the far-infrared (FIR) radiation field independent of the presence of active galactic nuclei (AGN) over many orders-of-magnitude in FIR luminosity (L$_{\rm FIR}$). This indicates that this transition arises mainly due to star formation. In this paper, we present spatially ($\sim$0.5 arcsec corresponding to $\sim$1 kiloparsec) and spectrally resolved ($\sim$100 kms$^{-1}$) observations of p$-$$\rm H_{2}O$(202 $-$ 111) in a sample of four strong gravitationally lensed high-redshift galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA). In addition to increasing the sample of luminous ($ > $ $10^{12}$L$_{\odot}$) galaxies observed with $\rm H_{2}O$, this paper examines the L$_{\rm H_{2}O}$/L$_{\rm FIR}$ relation on resolved scales for the first time at high-redshift. We find that L$_{\rm H_{2}O}$ is correlated with L$_{\rm FIR}$ on both global and resolved kiloparsec scales within the galaxy in starbursts and AGN with average L$_{\rm H_{2}O}$/L$_{\rm FIR}$ =$2.76^{+2.15}_{-1.21}\times10^{-5}$. We find that the scatter in the observed L$_{\rm H_{2}O}$/L$_{\rm FIR}$ relation does not obviously correlate with the effective temperature of the dust spectral energy distribution (SED) or the molecular gas surface density. This is a first step in developing p$-$$\rm H_{2}O$(202 $-$ 111) as a resolved star formation rate (SFR) calibrator.

Aims: We present and study spatially resolved imaging obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) of multiple $^{12}$CO($J=$6$-$5, 8$-$7 and 9$-$8) and two H$_2$O(2$_{02}-$1$_{11}$ and 2$_{11}-$2$_{02}$) emission lines and cold dust continuum toward the gravitationally lensed dusty star forming galaxy SPT0346-52 at z=$5.656$. Methods: Using a visibility-domain source-plane reconstruction we probe the structure and dynamics of the different components of the interstellar medium (ISM) in this galaxy down to scales of 1 kpc in the source plane. Results: Measurements of the intrinsic sizes of the different CO emission lines indicate that the higher J transitions trace more compact regions in the galaxy. Similarly, we find smaller dust continuum intrinsic sizes with decreasing wavelength, based on observations at rest-frame 130, 300 and 450$\mu$m. The source shows significant velocity structure, and clear asymmetry where an elongated structure is observed in the source plane with significant variations in their reconstructed sizes. This could be attributed to a compact merger or turbulent disk rotation. The differences in velocity structure through the different line tracers, however, hint at the former scenario in agreement with previous [CII] line imaging results. Measurements of the CO line ratios and magnifications yield significant variations as a function of velocity, suggesting that modeling of the ISM using integrated values could be misinterpreted. Modeling of the ISM in SPT0346-52 based on delensed fluxes indicate a highly dense and warm medium, qualitatively similar to that observed in high redshift quasar hosts.

We describe a next major frontier in observational studies of galaxy evolution and star formation: linking the physical conditions in the cold, star-forming interstellar medium to host galaxy and local environment. The integrated gas content of galaxies has been surveyed extensively over the last decades. The link between environment and cold gas density, turbulence, excitation, dynamical state, and chemical makeup remain far less well understood. We know that these properties do vary dramatically and theoretical work posits a strong connection between the state of the gas, its ability to form stars, and the impact of stellar feedback. A next major step in the field will be to use sensitive cm-, mm-, and submm-wave spectroscopy and high resolution spectroscopic imaging to survey the state of cold gas across the whole local galaxy population. Such observations have pushed the capabilities of the current generation of telescopes. We highlight three critical elements for progress in the next decade: (1) robust support and aggressive development of ALMA, (2) the deployment of very large heterodyne receiver arrays on single dish telescopes, and (3) development of a new interferometric array that dramatically improves on current capabilities at cm- and mm-wavelengths (~ 1-115 GHz).

We introduce the Simba simulations, the next generation of the Mufasa cosmological galaxy formation simulations run with Gizmo's meshless finite mass hydrodynamics. Simba includes updates to Mufasa's sub-resolution star formation and feedback prescriptions, and introduces black hole growth via the torque-limited accretion model of Angl\'es-Alc\'azar et al. (2017) from cold gas and Bondi accretion from hot gas, along with black hole feedback via kinetic bipolar outflows and X-ray energy. Ejection velocities are taken to be ~10^3 km/s at high Eddington ratios, increasing to ~8000 km/s at Eddington ratios below 2%, with a constant momentum input of 20L/c. Simba further includes an on-the-fly dust production, growth, and destruction model. Our Simba run with (100 Mpc/h)^3 and 1024^3 gas elements reproduces numerous observables, including galaxy stellar mass functions at z=0-6, the stellar mass--star formation rate main sequence, HI and H2 fractions, the mass-metallicity relation at z=0 and z=2, star-forming galaxy sizes, hot gas fractions in massive halos, and z=0 galaxy dust properties. However, Simba also yields an insufficiently sharp truncation of the z=0 mass function, and too-large sizes for low-mass quenched galaxies. We show that Simba's jet feedback is primarily responsible for quenching massive galaxies.

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

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

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 cycling of baryons in and out of galaxies is what ultimately drives galaxy formation and evolution. The circumgalactic medium (CGM) represents the interface between the interstellar medium and the cosmic web, hence its properties are directly shaped by the baryon cycle. Although traditionally the CGM is thought to consist of warm and hot gas, recent breakthroughs are presenting a new scenario according to which an important fraction of its mass may reside in the cold atomic and molecular phase. This would represent fuel that is readily available for star formation, with crucial implications for feeding and feedback processes in galaxies. However, such cold CGM, especially in local galaxies where its projected size on sky is expected to be of several arcminutes, cannot be imaged by ALMA due to interferometric spatial scale filtering of large-scale structures. We show that the only way to probe the multiphase CGM including its coldest component is through a large (e.g. 50-m) single dish (sub-)mm telescope.

One of the greatest challenges to theoretical models of massive galaxy formation is the regulation of star formation at early times. The relative roles of molecular gas expulsion, depletion, and stabilization are uncertain as direct observational constraints of the gas reservoirs in quenched or quenching galaxies at high redshift are scant. We present ALMA observations of CO(2-1) in a massive ($\log M_{\star}/M_{\odot}=11.2$), recently quenched galaxy at $z=1.522$. The optical spectrum of this object shows strong Balmer absorption lines, which implies that star formation ceased $\sim$0.8 Gyr ago. We do not detect CO(2-1) line emission, placing an upper limit on the molecular $\mathrm{H_2}$ gas mass of 1.1$\times10^{10}\,M_{\odot}$. The implied gas fraction is $f_{\rm{H_2}}{\equiv M_{H_2}/M_{\star}}<7\%$, $\sim10\times$ lower than typical star forming galaxies at similar stellar masses at this redshift, among the lowest gas fractions at this specific star formation rate at any epoch, and the most stringent constraint on the gas contents of a $z>1$ passive galaxy to date. Our observations show that the depletion of $\mathrm{H_2}$ from the interstellar medium of quenched objects can be both efficient and fairly complete, in contrast to recent claims of significant cold gas in recently quenched galaxies. We explore the variation in observed gas fractions in high-$z$ galaxies and show that galaxies with high stellar surface density have low $f_{\rm{H_2}}$, similar to recent correlations between specific star formation rate and stellar surface density.

We present a comprehensive suite of high-resolution (parsec-scale), idealised (non-cosmological) galaxy merger simulations (24 runs, stellar mass ratio ~2.5:1) to investigate the connection between interaction-induced star formation and the evolution of the interstellar medium (ISM) in various temperature-density regimes. We use the GIZMO code and the second version of the 'Feedback in Realistic Environments' model (FIRE-2), which captures the multi-phase structure of the ISM. Our simulations are designed to represent galaxy mergers in the local Universe. In this work, we focus on the 'galaxy-pair period' between first and second pericentric passage. We split the ISM into four regimes: hot, warm, cool and cold-dense, motivated by the hot, ionised, atomic and molecular gas phases observed in real galaxies. We find that, on average, interactions enhance the star formation rate of the pair (~30%, merger-suite sample average) and elevate their cold-dense gas content (~18%). This is accompanied by a decrease in warm gas (~11%), a negligible change in cool gas (~4% increase), and a substantial increase in hot gas (~400%). The amount of cold-dense gas with densities above 1000 cm^3 (the cold ultra-dense regime) is elevated significantly (~240%), but only accounts for 0.15% (on average) of the cold-dense gas budget.

SPT0346-52 is one of the most most luminous and intensely star-forming galaxies in the universe, with L_FIR > 10^13 L_sol and Sigma_SFR ~ 4200 M_sol yr^-1 kpc^-2. In this paper, we present ~0.15'' ALMA observations of the [CII]158micron emission line in this z=5.7 dusty star-forming galaxy. We use a pixellated lensing reconstruction code to spatially and kinematically resolve the source-plane [CII] and rest-frame 158 micron dust continuum structure at ~700 pc (~0.12'') resolution. We discuss the [CII] deficit with a pixellated study of the L_[CII]/L_FIR ratio in the source plane. We find that individual pixels within the galaxy follow the same trend found using unresolved observations of other galaxies, indicating that the deficit arises on scales <700 pc. The lensing reconstruction reveals two spatially and kinematically separated components (~1 kpc and ~500 km s^-1 apart) connected by a bridge of gas. Both components are found to be globally unstable, with Toomre Q instability parameters << 1 everywhere. We argue that SPT0346-52 is undergoing a major merger, which is likely driving the intense and compact star formation.

We use a suite of cosmological zoom galaxy formation simulations and dust radiative transfer calculations to explore the use of the monochromatic $850~\mu m$ luminosity (L$_{\rm \nu,850}$) as a molecular gas mass (M$_{\rm mol}$) estimator in galaxies between $0 < z < 9.5$ for a broad range of masses. For our fiducial simulations, where we assume the dust mass is linearly related to the metal mass, we find that empirical L$_{\rm \nu,850}$-M$_{\rm mol}$ calibrations accurately recover the molecular gas mass of our model galaxies, and that the L$_{\rm \nu,850}$-dependent calibration is preferred. We argue the major driver of scatter in the L$_{\rm \nu,850}$-M$_{\rm mol}$ relation arises from variations in the molecular gas to dust mass ratio, rather than variations in the dust temperature, in agreement with the previous study of Liang et al. Emulating a realistic measurement strategy with ALMA observing bands that are dependent on the source redshift, we find that estimating S$_{\rm \nu,850}$ from continuum emission at a different frequency contributes $10-20\%$ scatter to the L$_{\rm \nu,850}$-M$_{\rm mol}$ relation. This additional scatter arises from a combination of mismatches in assumed T$_{dust}$ and $\beta$ values, as well as the fact that the SEDs are not single-temperature blackbodies.Finally we explore the impact of a dust prescription in which the dust-to-metals ratio varies with metallicity. Though the resulting mean dust temperatures are $\sim50\%$ higher, the dust mass is significantly decreased for low-metallicity halos. As a result, the observationally calibrated L$_{\rm \nu,850}$-M$_{\rm mol}$ relation holds for massive galaxies, independent of the dust model, but below L$_{\rm \nu,850}\lesssim10^{28}$ erg s$^{-1}$ (metallicities $\log_{10}({\rm Z}/{\rm Z}_{\odot})\lesssim -0.8$) we expect galaxies may deviate from literature observational calibrations by $\gtrsim0.5$ dex.

The advent of new sub-millimeter observational facilities has stimulated the desire to model the sub-mm line emission of galaxies within cosmological galaxy formation models. This is typically done in post-processing by applying sub-resolution recipes to describe the properties of the unresolved interstellar medium. At the same time, while there is freedom in how one implements these sorts of recipes, the impact of various choices has yet to be systematically explored in simulations. In this paper, we do just that. We combine a semi-analytic model of galaxy formation with chemical equilibrium networks and numerical radiative transfer models and explore how different choices for the sub-grid modeling affect the predicted CO, [CI], and [CII] emission of galaxies. We find that a key component for a successful model includes a molecular cloud mass-size relation and scaling for the ultraviolet and cosmic ray radiation field that varies based on local ISM properties. Our most successful model adopts a Plummer profile for the radial density distribution of gas within molecular clouds. On the other hand, different assumptions for the clumping of gas within molecular clouds and changes in the slope of the molecular cloud mass distribution function hardly affect the CO, [CI], and [CII] luminosities of galaxies. At fixed star-formation rate the [CII]-SFR ratio of galaxies scales inversely with the pressure acting on molecular clouds, increasing the molecular clouds density and hence decreasing the importance of [CII] line cooling. Overall we find that it is essential that a wide range of sub-mm emission lines arising in vastly different phases of the ISM are used as model constraints in order to limit the freedom in sub-grid choices and we present the details of the most successful model variant.

The origin of the high SFR observed in high-z dusty star-forming galaxies is still unknown. Large fractions of dense molecular gas might provide part of the explanation, but there are few observational constraints on the amount of dense gas in high-z systems dominated by star formation. We present the results of our ALMA program targeting dense-gas tracers (HCN(5-4), HCO+(5-4), and HNC(5-4)) in 5 strongly lensed galaxies from the SPT SMG sample. We detected two of these lines (SNR>5) in SPT-125-47 at z=2.51 and tentatively detected all three (SNR~3) in SPT0551-50 at z=3.16. Since a significant fraction of our target lines is not detected, we developed a statistical method to derive unbiased mean properties taking into account both detections and non-detections. On average, the HCN(5-4) and HCO+(5-4) luminosities of our sources are a factor of ~1.7 fainter than expected, based on the local L'HCN(5-4)-LIR relation, but this offset corresponds to only ~2 sigma. We find that both the HCO+/HCN and HNC/HCN flux ratios are compatible with unity. The first ratio is expected for PDRs while the second is consistent with PDRs or XDRs and/or mid-IR pumping of HNC. Our sources are at the high end of the local relation between the star formation efficiency, determined using the LIR/[CI] and LIR/CO ratios, and the dense gas fraction, estimated using the HCN/[CI] and HCN/CO ratios. In SPT0125-47, we found that the velocity profiles of the lines tracing dense (HCN, HCO+) and lower-density (CO, [CI]) gas are similar. In addition to these lines, we obtained one robust and one tentative detection of 13CO(4-3) and found an average I12CO(4-3)/I13CO(4-3) flux ratio of 26.1$_{-3.5}^{+4.5}$, indicating a young but not pristine interstellar medium. We argue that the combination of large and slightly enriched gas reservoirs and high dense-gas fractions could explain the prodigious star formation in these systems.

The $\rm H_2$ mass of molecular clouds has traditionally been traced by the CO(J=1-0) rotational transition line. This said, CO is relatively easily photodissociated, and can also be destroyed by cosmic rays, thus rendering some fraction of molecular gas to be "CO-dark". We investigate the amount and physical properties of CO-dark gas in two $z \sim 0$ disc galaxies, and develop predictions for the expected intensities of promising alternative tracers ([CI 609 $\mu$m and [CII] 158 $\mu$m emission). We do this by combining cosmological zoom simulations of disc galaxies with thermal-radiative-chemical equilibrium interstellar medium (ISM) calculations to model the predicted H~\textsc{i} and $\rm H_2$ abundances and CO(J=1-0), [CI] 609 $\mu$m and [CII] 158 $\mu$m emission properties. Our model treats the ISM as a collection of radially stratified clouds whose properties are dictated by their volume and column densities, the gas-phase metallicity, and the interstellar radiation field and cosmic ray ionization rates. Our main results follow. Adopting an observationally motivated definition of CO-dark gas, i.e. $\rm H_2$ gas with $W_{\rm CO} < 0.1 $ K-km/s, we find that a significant amount ($> 50\%$) of the total $\rm H_2$ mass lies in CO-dark gas, most of which is diffuse gas, poorly shielded due to low dust column density. The CO-dark molecular gas tends to be dominated by [CII], though [CI] also serves as a bright tracer of the dark gas in many instances. At the same time, [CII] also tends to trace neutral atomic gas. As a result, when we quantify the conversion factors for the three carbon-based tracers of molecular gas, we find that [CI] suffers the least contamination from diffuse atomic gas, and is relatively insensitive to secondary parameters.

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

The goal of this science case is to 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 Next-Generation Very Large Array (ngVLA) will revolutionize our understanding of the Early Universe by tracing the coldest phase of molecular gas -the raw ingredient for star formation- in the most distant galaxies and galaxy-clusters. The km-scale core of the ngVLA will be densely packed with antennas, making it a prime instrument for imaging low-surface-brightness emission from large-scale molecular gas in the high-z circum- and inter-galactic medium (CGM/IGM). Recent studies indicate that large amounts of cold molecular gas are hiding in the 10s-100 kpc environments of distant galaxies, but that technical limitations on existing telescope arrays have prevented us from efficiently detecting these large molecular reservoirs. This may have led to a severely biased view of the molecular Universe. We present the science case for low-surface-brightness CO observations of the Early Universe, and how the core of the ngVLA will reveal the cold molecular Universe to limits and at scales not currently detectable by radio telescopes. As such, the ngVLA core will be a powerful instrument for studying the cold baryon cycle that drives the early evolution of galaxies and clusters.

The Origins Space Telescope, one of four large Mission Concept studies sponsored by NASA for review in the 2020 US Astrophysics Decadal Survey, will open unprecedented discovery space in the infrared, unveiling our cosmic origins. We briefly describe in this article the key science themes and architecture for OST. With a sensitivity gain of up to a factor of 1,000 over any previous or planned mission, OST will open unprecedented discovery space, allow us to peer through an infrared window teeming with possibility. OST will fundamentally change our understanding of our cosmic origins - from the growth of galaxies and black holes, to uncovering the trail of water, to life signs in nearby Earth-size planets, and discoveries never imagined. Built to be highly adaptable, while addressing key science across many areas of astrophysics, OST will usher in a new era of infrared astronomy.

Galaxies grow inefficiently, with only a few percent of the available gas converted into stars each free-fall time. Feedback processes, such as outflowing winds driven by radiation pressure, supernovae or supermassive black hole accretion, can act to halt star formation if they heat or expel the gas supply. We report a molecular outflow launched from a dust-rich star-forming galaxy at redshift 5.3, one billion years after the Big Bang. The outflow reaches velocities up to 800 km/s relative to the galaxy, is resolved into multiple clumps, and carries mass at a rate within a factor of two of the star formation rate. Our results show that molecular outflows can remove a large fraction of the gas available for star formation from galaxies at high redshift.

Post-starburst or "E+A" galaxies are rapidly transitioning from star-forming to quiescence. While the current star formation rate of post-starbursts is already at the level of early type galaxies, we recently discovered that many have large CO-traced molecular gas reservoirs consistent with normal star forming galaxies. These observations raise the question of why these galaxies have such low star formation rates. Here we present an ALMA search for the denser gas traced by HCN (1--0) and HCO+ (1--0) in two CO-luminous, quiescent post-starburst galaxies. Intriguingly, we fail to detect either molecule. The upper limits are consistent with the low star formation rates and with early-type galaxies. The HCN/CO luminosity ratio upper limits are low compared to star-forming and even many early type galaxies. This implied low dense gas mass fraction explains the low star formation rates relative to the CO-traced molecular gas and suggests the state of the gas in post-starburst galaxies is unusual, with some mechanism inhibiting its collapse to denser states. We conclude that post-starbursts galaxies are now quiescent because little dense gas is available, in contrast to the significant CO-traced lower density gas reservoirs that still remain.

In this paper, we provide a physical model for the origin of variations in the shapes and bump strengths of dust attenuation laws in galaxies by combining a large suite of cosmological "zoom-in" galaxy formation simulations with 3D Monte Carlo dust radiative transfer calculations. We model galaxies over 3 orders of magnitude in stellar mass, ranging from Milky Way like systems through massive galaxies at high-redshift. Critically, for these calculations we employ a constant underlying dust extinction law in all cases, and examine how the role of geometry and radiative transfer effects impact the resultant attenuation curves. Our main results follow. Despite our usage of a constant dust extinction curve, we find dramatic variations in the derived attenuation laws. The slopes of normalized attenuation laws depend primarily on the complexities of star-dust geometry. Increasing fractions of unobscured young stars flatten normalized curves, while increasing fractions of unobscured old stars steepen curves. Similar to the slopes of our model attenuation laws, we find dramatic variation in the 2175 Angstrom ultraviolet (UV) bump strength, including a subset of curves with little to no bump. These bump strengths are primarily influenced by the fraction of unobscured O and B stars in our model, with the impact of scattered light having only a secondary effect. Taken together, these results lead to a natural relationship between the attenuation curve slope and 2175 Angstrom bump strength. Finally, we apply these results to a 25 Mpc/h box cosmological hydrodynamic simulation in order to model the expected dispersion in attenuation laws at integer redshifts from z=0-6. A significant dispersion is expected at low redshifts, and decreases toward z=6. We provide tabulated results for the best fit median attenuation curve at all redshifts.

The process by which massive galaxies transition from blue, star-forming disks into red, quiescent galaxies remains one of the most poorly-understood aspects of galaxy evolution. In this investigation, we attempt to gain a better understanding of how star formation is quenched by focusing on a massive post-starburst galaxy at z = 0.747. The target has a high stellar mass and a molecular gas fraction of ~30% -- unusually high for its low star formation rate. We look for indicators of star formation suppression mechanisms in the stellar kinematics and age distribution of the galaxy obtained from spatially resolved Gemini Integral-Field spectra and in the gas kinematics obtained from ALMA. We find evidence of significant rotation in the stars, but we do not detect a stellar age gradient within 5 kpc. The molecular gas is aligned with the stellar component, and we see no evidence of strong gas outflows. Our target may represent the product of a merger-induced starburst or of morphological quenching; however, our results are not completely consistent with any of the prominent quenching models.

Massive galaxy clusters are now found as early as 3 billion years after the Big Bang, containing stars that formed at even earlier epochs. The high-redshift progenitors of these galaxy clusters, termed 'protoclusters', are identified in cosmological simulations with the highest dark matter overdensities. While their observational signatures are less well defined compared to virialized clusters with a substantial hot intra-cluster medium (ICM), protoclusters are expected to contain extremely massive galaxies that can be observed as luminous starbursts. Recent claimed detections of protoclusters hosting such starbursts do not support the kind of rapid cluster core formation expected in simulations because these structures contain only a handful of starbursting galaxies spread throughout a broad structure, with poor evidence for eventual collapse into a protocluster. Here we report that the source SPT2349-56 consists of at least 14 gas-rich galaxies all lying at z = 4.31 based on sensitive observations of carbon monoxide and ionized carbon. We demonstrate that each of these galaxies is forming stars between 50 and 1000 times faster than our own Milky Way, and all are located within a projected region only $\sim$ 130 kiloparsecs in diameter. This galaxy surface density is more than 10 times the average blank field value (integrated over all redshifts) and $>$1000 times the average field volume density. The velocity dispersion ($\sim$ 410 km s$^{-1}$) of these galaxies and enormous gas and star formation densities suggest that this system represents a galaxy cluster core at an advanced stage of formation when the Universe was only 1.4 billion years old. A comparison with other known protoclusters at high redshifts shows that SPT2349-56 is a uniquely massive and dense system that could be building one of the most massive structures in the Universe today.

Modeling of the spectral line energy distribution (SLED) of the CO molecule can reveal the physical conditions (temperature, density) of molecular gas in Galactic clouds and other galaxies. Recently, the Herschel Space Observatory and ALMA have offered, for the first time, a comprehensive view of the rotational J = 4-3 through J = 13-12 lines, which arise from a complex, diverse range of physical conditions that must be simplified to one, two, or three components when modeled. Here we investigate the recoverability of physical conditions from SLEDs produced by galaxy evolution simulations containing a large dynamical range in physical properties. These simulated SLEDs were generally fit well by one component of gas whose properties largely resemble or slightly underestimate the luminosity-weighted properties of the simulations when clumping due to non-thermal velocity dispersion is taken into account. If only modeling the first three rotational lines, the median values of the marginalized parameter distributions better represent the luminosity-weighted properties of the simulations, but the uncertainties in the fitted parameters are nearly an order of magnitude, compared to approximately 0.2 dex in the "best-case" scenario of a fully sampled SLED through J = 10-9. This study demonstrates that while common CO SLED modeling techniques cannot reveal the underlying complexities of the molecular gas, they can distinguish bulk luminosity-weighted properties that vary with star formation surface densities and galaxy evolution, if a sufficient number of lines are detected and modeled.

Non-parametric morphology measures are a powerful tool for identifying galaxy mergers at low redshifts. We employ cosmological zoom simulations using Gizmo with the Mufasa feedback scheme, post-processed using 3D dust radiative transfer into mock observations, to study whether common morphological measures Gini G, M20, concentration C, and asymmetry A are effective at identifying major galaxy mergers at z ~ 2 - 4, i.e. "Cosmic Noon". Our zoom suite covers galaxies with 10^8.6 < M_* < 10^11 M_sun at z ~ 2, and broadly reproduces key global galaxy observations. Our primary result is that these morphological measures are unable to robustly pick out galaxies currently undergoing mergers during Cosmic Noon, typically performing no better than a random guess. This improves only marginally if we consider whether galaxies have undergone a merger within the last Gyr. When also considering minor mergers, galaxies display no trend of moving towards the merger regime with increasing merger ratio. From z = 4 -> 2, galaxies move from the non-merger towards the merger regime in all statistics, but this is primarily an effect of mass: Above a given noise level, higher mass galaxies display a more complex outer morphology induced by their clustered environment. We conclude that during Cosmic Noon, these morphological statistics are of limited value in identifying galaxy mergers.

According to the current understanding of cosmic structure formation, the precursors of the most massive structures in the Universe began to form shortly after the Big Bang, in regions corresponding to the largest fluctuations in the cosmic density field. Observing these structures during their period of active growth and assembly - the first few hundred million years of the Universe - is challenging because it requires surveys that are sensitive enough to detect the distant galaxies that act as signposts for these structures and wide enough to capture the rarest objects. As a result, very few such objects have been detected so far. Here we report observations of a far-infrared-luminous object at redshift 6.900 (less than 800 Myr after the Big Bang) that was discovered in a wide-field survey. High-resolution imaging reveals this source to be a pair of extremely massive star-forming galaxies. The larger of these galaxies is forming stars at a rate of 2900 solar masses per year, contains 270 billion solar masses of gas and 2.5 billion solar masses of dust, and is more massive than any other known object at a redshift of more than 6. Its rapid star formation is probably triggered by its companion galaxy at a projected separation of just 8 kiloparsecs. This merging companion hosts 35 billion solar masses of stars and has a star-formation rate of 540 solar masses per year, but has an order of magnitude less gas and dust than its neighbor and physical conditions akin to those observed in lower-metallicity galaxies in the nearby Universe. These objects suggest the presence of a dark-matter halo with a mass of more than 400 billion solar masses, making it among the rarest dark-matter haloes that should exist in the Universe at this epoch.

Of the almost 40 star forming galaxies at z>~5 (not counting QSOs) observed in [CII] to date, nearly half are either very faint in [CII], or not detected at all, and fall well below expectations based on locally derived relations between star formation rate (SFR) and [CII] luminosity. Combining cosmological zoom simulations of galaxies with SIGAME (SImulator of GAlaxy Millimeter/submillimeter Emission) we have modeled the multi-phased interstellar medium (ISM) and its emission in [CII], [OI] and [OIII], from 30 main sequence galaxies at z~6 with star formation rates ~3-23Msun/yr, stellar masses ~(0.7-8)x10^9Msun, and metallicities ~(0.1-0.4)xZsun. The simulations are able to reproduce the aforementioned [CII]-faintness at z>5, match two of the three existing z>~5 detections of [OIII], and are furthermore roughly consistent with the [OI] and [OIII] luminosity relations with SFR observed for local starburst galaxies. We find that the [CII] emission is dominated by the diffuse ionized gas phase and molecular clouds, which on average contribute ~66% and ~27%, respectively. The molecular gas, which constitutes only ~10% of the total gas mass is thus a more efficient emitter of [CII] than the ionized gas making up ~85% of the total gas mass. A principal component analysis shows that the [CII] luminosity correlates with the star formation activity as well as average metallicity. The low metallicities of our simulations together with their low molecular gas mass fractions can account for their [CII]-faintness, and we suggest these factors may also be responsible for the [CII]-faint normal galaxies observed at these early epochs.

The physical mechanisms that quench star formation, turning blue star-forming galaxies into red quiescent galaxies, remain unclear. In this Letter, we investigate the role of gas supply in suppressing star formation by studying the molecular gas content of post-starburst galaxies. Leveraging the wide area of the SDSS, we identify a sample of massive intermediate-redshift galaxies that have just ended their primary epoch of star formation. We present ALMA CO(2-1) observations of two of these post-starburst galaxies at z~0.7 with M* ~ 2x10^11 Msun. Their molecular gas reservoirs of (6.4 +/- 0.8) x 10^9 Msun and (34.0 +/- 1.6) x 10^9 Msun are an order of magnitude larger than comparable-mass galaxies in the local universe. Our observations suggest that quenching does not require the total removal or depletion of molecular gas, as many quenching models suggest. However, further observations are required both to determine if these apparently quiescent objects host highly obscured star formation and to investigate the intrinsic variation in the molecular gas properties of post-starburst galaxies.

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum and [C II] 158 $\mu$m fine structure line emission towards a far-infrared-luminous quasar, ULAS J131911.29$+$095051.4 at $z=6.13$, and combine the new Cycle 1 data with ALMA Cycle 0 data. The combined data have an angular resolution $\sim$ $0.3$, and resolve both the dust continuum and the [C II] line emission on few kpc scales. The [C II] line emission is more irregular than the dust continuum emission which suggests different distributions between the dust and [C II]-emitting gas. The combined data confirm the [C II] velocity gradient that we previously detected in lower resolution ALMA image from Cycle 0 data alone. We apply a tilted ring model to the [C II] velocity map to obtain a rotation curve, and constrain the circular velocity to be 427 $\pm$ 55 km s$^{-1}$ at a radius of 3.2 kpc with an inclination angle of 34$^\circ$. We measure the dynamical mass within the 3.2 kpc region to be 13.4$_{-5.3}^{+7.8}$ $\times 10^{10}\,M_{\odot}$. This yields a black hole and host galaxy mass ratio of 0.020$_{-0.007}^{+0.013}$, which is about 4$_{-2}^{+3}$ times higher than the present-day $M_{\rm BH}$/$M_{\rm bulge}$ ratio. This suggests that the supermassive black hole grows the bulk of its mass before the formation of the most of stellar mass in this quasar host galaxy in the early universe.

It is not yet known if the properties of molecular gas in distant protocluster galaxies are significantly affected by their environment as galaxies are in local clusters. Through a deep, 64 hours of effective on-source integration with the Australian Telescope Compact Array (ATCA), we discovered a massive, M_mol=2.0+-0.2x10^11 M_sun, extended, ~40 kpc, CO(1-0)-emitting disk in the protocluster surrounding the radio galaxy, MRC1138-262. The galaxy, at z_CO=2.1478, is a clumpy, massive disk galaxy, M_star~5x10^11 M_sun, which lies 250 kpc in projection from MRC1138-262 and is a known H-alpha emitter, HAE229. HAE229 has a molecular gas fraction of ~30%. The CO emission has a kinematic gradient along its major axis, centered on the highest surface brightness rest-frame optical emission, consistent with HAE229 being a rotating disk. Surprisingly, a significant fraction of the CO emission lies outside of the UV/optical emission. In spite of this, HAE229 follows the same relation between star-formation rate and molecular gas mass as normal field galaxies. HAE229 is the first CO(1-0) detection of an ordinary, star-forming galaxy in a protocluster. We compare a sample of cluster members at z>0.4 that are detected in low-order CO transitions with a similar sample of sources drawn from the field. We confirm findings that the CO-luminosity and FWHM are correlated in starbursts and show that this relation is valid for normal high-z galaxies as well as those in overdensities. We do not find a clear dichotomy in the integrated Schmidt-Kennicutt relation for protocluster and field galaxies. Not finding any environmental dependence in the "star-formation efficiency" or the molecular gas content, especially for such an extended CO disk, suggests that environmentally-specific processes such as ram pressure stripping are not operating efficiently in (proto)clusters. (abridged)

We report the discovery and constrain the physical conditions of the interstellar medium of the highest-redshift millimeter-selected dusty star-forming galaxy (DSFG) to date, SPT-S J031132-5823.4 (hereafter SPT0311-58), at $z=6.900 +/- 0.002$. SPT0311-58 was discovered via its 1.4mm thermal dust continuum emission in the South Pole Telescope (SPT)-SZ survey. The spectroscopic redshift was determined through an ALMA 3mm frequency scan that detected CO(6-5), CO(7-6) and [CI](2-1), and subsequently confirmed by detections of CO(3-2) with ATCA and [CII] with APEX. We constrain the properties of the ISM in SPT0311-58 with a radiative transfer analysis of the dust continuum photometry and the CO and [CI] line emission. This allows us to determine the gas content without ad hoc assumptions about gas mass scaling factors. SPT0311-58 is extremely massive, with an intrinsic gas mass of $M_{\rm gas} = 3.3 \pm 1.9 \times10^{11}\,M_{\odot}$. Its large mass and intense star formation is very rare for a source well into the Epoch of Reionization.

We utilise a series of high-resolution cosmological zoom simulations of galaxy formation to investigate the relationship between the ultraviolet (UV) slope, beta, and the ratio of the infrared luminosity to UV luminosity (IRX) in the spectral energy distributions (SEDs) of galaxies. We employ dust radiative transfer calculations in which the SEDs of the stars in galaxies propagate through the dusty interstellar medium. Our main goals are to understand the origin of, and scatter in the IRX-beta relation; to assess the efficacy of simplified stellar population synthesis screen models in capturing the essential physics in the IRX-beta relation; and to understand systematic deviations from the canonical local IRX-beta relations in particular populations of high-redshift galaxies. Our main results follow. Galaxies that have young stellar populations with relatively cospatial UV and IR emitting regions and a Milky Way-like extinction curve fall on or near the standard Meurer relation. This behaviour is well captured by simplified screen models. Scatter in the IRX-beta relation is dominated by three major effects: (i) older stellar populations drive galaxies below the relations defined for local starbursts due to a reddening of their intrinsic UV SEDs; (ii) complex geometries in high-z heavily star forming galaxies drive galaxies toward blue UV slopes owing to optically thin UV sightlines; (iii) shallow extinction curves drive galaxies downward in the IRX-beta plane due to lowered NUV/FUV extinction ratios. We use these features of the UV slopes of galaxies to derive a fitting relation that reasonably collapses the scatter back toward the canonical local relation. Finally, we use these results to develop an understanding for the location of two particularly enigmatic populations of galaxies in the IRX-beta plane: z~2-4 dusty star forming galaxies, and z>5 star forming galaxies.

Observations of ionised carbon at 158 micron ([CII]) from luminous star-forming galaxies at z~0 show that their ratios of [CII] to far infrared (FIR) luminosity are systematically lower than those of more modestly star-forming galaxies. In this paper, we provide a theory for the origin of this so called "[CII] deficit" in galaxies. Our model treats the interstellar medium as a collection of clouds with radially-stratified chemical and thermal properties, which are dictated by the clouds' volume and surface densities, as well as the interstellar radiation and cosmic ray fields to which they are exposed. [CII] emission arises from the outer, HI dominated layers of clouds, and from regions where the hydrogen is H2 but the carbon is predominantly C+. In contrast, the most shielded regions of clouds are dominated by CO and produce little [CII] emission. This provides a natural mechanism to explain the observed [CII]-star formation relation: galaxies' star formation rates are largely driven by the surface densities of their clouds. As this rises, so does the fraction of gas in the CO-dominated phase that produces little [CII] emission. Our model further suggests that the apparent offset in the [CII]-FIR relation for high-z sources compared to those at present epoch may arise from systematically larger gas masses at early times: a galaxy with a large gas mass can sustain a high star formation rate even with relatively modest surface density, allowing copious [CII] emission to coexist with rapid star formation.

We present the results of the first, deep ALMA imaging covering the full 4.5 sq arcmin of the Hubble Ultra Deep Field (HUDF) as previously imaged with WFC3/IR on HST. Using a mosaic of 45 pointings, we have obtained a homogeneous 1.3mm image of the HUDF, achieving an rms sensitivity of 35 microJy, at a resolution of 0.7 arcsec. From an initial list of ~50 >3.5sigma peaks, a rigorous analysis confirms 16 sources with flux densities S(1.3) > 120 microJy. All of these have secure galaxy counterparts with robust redshifts (<z> = 2.15), and 12 are also detected at 6GHz in new deep JVLA imaging. Due to the wealth of supporting data in this unique field, the physical properties of the ALMA sources are well constrained, including their stellar masses (M*) and UV+FIR star-formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-z Universe; indeed at z > 2 our ALMA sample contains 7 of the 9 galaxies in the HUDF with M* > 2 x 10^10 Msun and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S(1.3) ~ 10 micro-Jy. We find strong evidence for a steep `main sequence' for star-forming galaxies at z ~ 2, with SFR \propto M* and a mean specific SFR = 2.2 /Gyr. Moreover, we find that ~85% of total star formation at z ~ 2 is enshrouded in dust, with ~65% of all star formation at this epoch occurring in high-mass galaxies (M* > 2 x 10^10 Msun), for which the average obscured:unobscured SF ratio is ~200. Finally, we combine our new ALMA results with the existing HST data to revisit the cosmic evolution of star-formation rate density; we find that this peaks at z ~ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured thereafter at z ~ 4.

We present $\simeq$0$.\!\!^{\prime\prime}4$-resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at $z = 1.3-3.0$. These galaxies are selected from sensitive, blank-field surveys of the $2' \times 2'$ Hubble Ultra-Deep Field at $\lambda = 5$ cm and 1.3 mm using the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA). They have star-formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at $z \sim 2$. Morphological classification performed on spatially-resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGN), thereby representing a diversity of $z \sim 2$ SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of $4.2 \pm 1.8$ kpc. This provides direct evidence for galaxy-wide star formation in distant, blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGN, which are found to be co-spatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies whose SFRs are $3-8$ times larger, providing a constraint on the characteristic SFR ($\sim300$ M$_{\odot}$yr$^{-1}$) above which a significant population of more compact star-forming galaxies appears to emerge.

This paper describes the beginning of the Far-Infrared Surveyor mission study for NASA's Astrophysics Decadal 2020. We describe the scope of the study, and the open process approach of the Science and Technology Definition Team. We are currently developing the science cases and provide some preliminary highlights here. We note key areas for technological innovation and improvements necessary to make a Far-Infrared Surveyor mission a reality.

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant LABs.

We present S\'IGAME simulations of the [CII]157.7$\mu$ fine structure line emission from cosmological smoothed particle hydrodynamics (SPH) simulations of seven main sequence galaxies at z=2. Using sub-grid physics prescriptions the gas in our simulations is modeled as a multi-phased interstellar medium (ISM) comprised of molecular gas residing in giant molecular clouds, an atomic gas phase associated with photo-dissociation regions (PDRs) at the cloud surfaces, and a diffuse, ionized gas phase. Adopting logotropic cloud density profiles and accounting for heating by the local FUV radiation field and cosmic rays by scaling both with local star formation rate (SFR) volume density, we calculate the [CII] emission using a photon escape probability formalism. The [CII] emission peaks in the central $\lesssim$1 kpc of our galaxies as do the SFR radial profiles, with most [CII] ($\gtrsim$70%) originating in the molecular gas phase, whereas further out ($\gtrsim$2 kpc), the atomic/PDR gas dominates ($\gtrsim$90%) the [CII] emission, no longer tracing on-going star formation. Throughout, the ionized gas contribution is negligible ($\lesssim$3%). The [CII] luminosity vs. SFR ([CII]-SFR) relationship, integrated as well as spatially resolved (on scales of 1 kpc), delineated by our simulated galaxies is in good agreement with the corresponding relations observed locally and at high redshifts. In our simulations, the molecular gas dominates the [CII] budget at SFR$\gtrsim$20 M$_{\odot}$yr$^{-1}$ (SFR density $\gtrsim$0.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$), while atomic/PDR gas takes over at lower SFRs, suggesting a picture in which [CII] predominantly traces the molecular gas in high-density/pressure regions where star formation is on-going, and otherwise reveals the atomic/PDR gas phase.

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.

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.

Submillimetre-luminous galaxies at high-redshift are the most luminous, heavily star-forming galaxies in the Universe, and are characterised by prodigious emission in the far-infrared at 850 microns (S850 > 5 mJy). They reside in halos ~ 10^13Msun, have low gas fractions compared to main sequence disks at a comparable redshift, trace complex environments, and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy which simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter halos have rising star formation histories that peak at collective rates ~ 500-1000 Msun/yr at z=2-3, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fueled in part by a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a duty cycle of nearly a gigayear, our simulations show that the submillimetre-luminous phase of high-z galaxies is a drawn out one that is associated with significant mass buildup in early Universe proto-clusters, and that many submillimetre-luminous galaxies are actually composed of numerous unresolved components (for which there is some observational evidence).

Post-starburst (or "E+A") galaxies are characterized by low H$\alpha$ emission and strong Balmer absorption, suggesting a recent starburst, but little current star formation. Although many of these galaxies show evidence of recent mergers, the mechanism for ending the starburst is not yet understood. To study the fate of the molecular gas, we search for CO (1-0) and (2-1) emission with the IRAM 30m and SMT 10m telescopes in 32 nearby ($0.01<z<0.12$) post-starburst galaxies drawn from the Sloan Digital Sky Survey. We detect CO in 17 (53%). Using CO as a tracer for molecular hydrogen, and a Galactic conversion factor, we obtain molecular gas masses of $M(H_2)=10^{8.6}$-$10^{9.8} M_\odot$ and molecular gas mass to stellar mass fractions of $\sim10^{-2}$-$10^{-0.5}$, comparable to those of star-forming galaxies. The large amounts of molecular gas rule out complete gas consumption, expulsion, or starvation as the primary mechanism that ends the starburst in these galaxies. The upper limits on $M(H_2)$ for the 15 undetected galaxies range from $10^{7.7} M_\odot$ to $10^{9.7} M_\odot$, with the median more consistent with early-type galaxies than with star-forming galaxies. Upper limits on the post-starburst star formation rates (SFRs) are lower by $\sim10\times$ than for star-forming galaxies with the same $M(H_2)$. We also compare the molecular gas surface densities ($\Sigma_{\rm H_2}$) to upper limits on the SFR surface densities ($\Sigma_{\rm SFR}$), finding a significant offset, with lower $\Sigma_{\rm SFR}$ for a given $\Sigma_{\rm H_2}$ than is typical for star-forming galaxies. This offset from the Kennicutt-Schmidt relation suggests that post-starbursts have lower star formation efficiency, a low CO-to-H$_2$ conversion factor characteristic of ULIRGs, and/or a bottom-heavy initial mass function, although uncertainties in the rate and distribution of current star formation remain.

We report ALMA Band 7 (350 GHz) imaging at 0.4 - 0.6arcsec resolution and Band 9 (696 GHz) at ~0.25arcsec resolution of the luminous IR galaxies Arp 220 and NGC 6240. The long wavelength dust continuum is used to estimate ISM masses for Arp 220 East, West and NGC 6240 of 1.9, 4.2 and 1.6x10^9 msun within radii of 69, 65 and 190 pc. The HCN emission was modeled to derive the emissivity distribution as a function of radius and the kinematics of each nuclear disk, yielding dynamical masses consistent with the masses and sizes derived from the dust emission. In Arp 220, the major dust and gas concentrations are at radii less than 50 pc in both counter-rotating nuclear disks. The thickness of the disks in Arp 220estimated from the velocity dispersion and rotation velocities are 10-20 pc and the mean gas densities are n_H2 ~10^5 cm^-3 at R < 50 pc. We develop an analytic treatment for the molecular excitation (including photon trapping), yielding volume densities for both the HCN and CS emission with n_H2 ~2x10^5 cm^-3. The agreement of the mean density from the total mass and size with that required for excitation suggests that the volume is essentially filled with dense gas, i.e. it is not cloudy or like swiss cheese.

We present < 1 kpc resolution CO imaging study of 37 optically-selected local merger remnants using new and archival interferometric maps obtained with ALMA, CARMA, SMA and PdBI. We supplement a sub-sample with single-dish measurements obtained at the NRO 45 m telescope for estimating the molecular gas mass (10^7 - 10^11 M_sun), and evaluating the missing flux of the interferometric measurements. Among the sources with robust CO detections, we find that 80 % (24/30) of the sample show kinematical signatures of rotating molecular gas disks (including nuclear rings) in their velocity fields, and the sizes of these disks vary significantly from 1.1 kpc to 9.3 kpc. The size of the molecular gas disks in 54 % of the sources is more compact than the K-band effective radius. These small gas disks may have formed from a past gas inflow that was triggered by a dynamical instability during a potential merging event. On the other hand, the rest (46 %) of the sources have gas disks which are extended relative to the stellar component, possibly forming a late-type galaxy with a central stellar bulge. Our new compilation of observational data suggests that nuclear and extended molecular gas disks are common in the final stages of mergers. This finding is consistent with recent major-merger simulations of gas rich progenitor disks. Finally, we suggest that some of the rotation-supported turbulent disks observed at high redshifts may result from galaxies that have experienced a recent major merger.

Observations of molecular gas in high-z star-forming galaxies typically rely on emission from CO lines arising from states with rotational quantum numbers J > 1. Converting these observations to an estimate of the CO J=1-0 intensity, and thus inferring H2 gas masses, requires knowledge of the CO excitation ladder, or spectral line energy distribution (SLED). The few available multi-J CO observations of galaxies show a very broad range of SLEDs, even at fixed galaxy mass and star formation rate, making the conversion to J=1-0 emission and hence molecular gas mass highly uncertain. Here, we combine numerical simulations of disk galaxies and galaxy mergers with molecular line radiative transfer calculations to develop a model for the physical parameters that drive variations in CO SLEDs in galaxies. An essential feature of our model is a fully self-consistent computation of the molecular gas temperature and excitation structure. We find that, while the shape of the SLED is ultimately determined by difficult-to-observe quantities such as the gas density, temperature, and optical depth distributions, all of these quantities are well-correlated with the galaxy's mean star formation rate surface density (Sigma_SFR), which is observable. We use this result to develop a model for the CO SLED in terms of Sigma_SFR, and show that this model quantitatively reproduces the SLEDs of galaxies over a dynamic range of ~200 in SFR surface density, at redshifts from z=0-6. This model should make it possible to significantly reduce the uncertainty in deducing molecular gas masses from observations of high-J CO emission.

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.

The CO-H2 conversion factor (Xco; otherwise known as the X-factor) is observed to be remarkably constant in the Milky Way and in the Local Group (aside from the SMC). To date, our understanding of why Xco should be so constant remains poor. Using a combination of extremely high resolution (~ 1 pc) galaxy evolution simulations and molecular line radiative transfer calculations, we suggest that Xco displays a narrow range of values in the Galaxy due to the fact that molecular clouds share very similar physical properties. In our models, this is itself a consequence of stellar feedback competing against gravitational collapse. GMCs whose lifetimes are regulated by radiative feedback show a narrow range of surface densities, temperatures and velocity dispersions with values comparable to those seen in the Milky Way. As a result, the X-factors from these clouds show reasonable correspondence with observed data from the Local Group, and a relatively narrow range. On the other hand, feedback-free clouds collapse to surface densities that are larger than those seen in the Galaxy, and hence result in X-factors that are systematically too large compared to the Milky Way's. We conclude that radiative feedback within GMCs can generate cloud properties similar to those observed in the Galaxy, and hence a roughly constant Milky Way X-factor in normal, quiescent clouds.

We examine the cosmic evolution of a stellar initial mass function (IMF) in galaxies that varies with the Jeans mass in the interstellar medium, paying particular attention to the K-band stellar mass to light ratio (M/L_K) of present-epoch massive galaxies. We calculate the typical Jeans mass using high-resolution hydrodynamic simulations coupled with a fully radiative model for the ISM, which yields a parameterisation of the IMF characteristic mass as a function of galaxy star formation rate (SFR). We then calculate the star formation histories of galaxies utilising an equilibrium galaxy growth model coupled with constraints on the star formation histories set by abundance matching models. We find that at early times, energetic coupling between dust and gas drive warm conditions in the ISM, yielding bottom-light/top- heavy IMFs associated with large ISM Jeans masses for massive star-forming galaxies. Owing to the remnants of massive stars that formed during the top-heavy phases at early times, the resultant M/L_K(sigma) in massive galaxies at the present epoch is increased relative to the non- varying IMF case. At late times, lower cosmic ray fluxes allow for cooler ISM temperatures in massive galaxies, and hence newly formed clusters will exhibit bottom-heavy IMFs, further increasing M/L_K(sigma). Our central result is hence that a given massive galaxy may go through both top-heavy and bottom-heavy IMF phases during its lifetime, though the bulk of the stars form during a top-heavy phase. Qualitatively, the variations in M/L_K(sigma) with galaxy mass are in agreement with observations, however, our model may not be able to account for bottom-heavy mass functions as indicated by stellar absorption features.

We present ALMA observations of the [C II] 158 micron fine structure line and dust continuum emission from the host galaxies of five redshift 6 quasars. We also report complementary observations of 250 GHz dust continuum and CO (6-5) line emission from the z=6.00 quasar SDSS J231038.88+185519.7. The ALMA observations were carried out in the extended array at 0.7" resolution. We have detected the line and dust continuum in all five objects. The derived [C II] line luminosities are 1.6x10^{9} to 8.8x10^{9} Lsun and the [C II]-to-FIR luminosity ratios are 3.0-5.6x10^{-4}, which is comparable to the values found in other high-redshift quasar-starburst systems and local ultra-luminous infrared galaxies. The sources are marginally resolved and the intrinsic source sizes (major axis FWHM) are constrained to be 0.3" to 0.6" (i.e., 1.7 to 3.5 kpc) for the [C II] line emission and 0.2" to 0.4" (i.e., 1.2 to 2.3 kpc) for the continuum. These measurements indicate that there is vigorous star formation over the central few kpc in the quasar host galaxies. The ALMA observations also constrain the dynamical properties of the atomic gas in the starburst nuclei. The intensity-weighted velocity maps of three sources show clear velocity gradients. Such velocity gradients are consistent with a rotating, gravitationally bound gas component, although they are not uniquely interpreted as such. Under the simplifying assumption of rotation, the implied dynamical masses within the [C II]-emitting regions are of order 10^{10} to 10^{11} Msun. Given these estimates, the mass ratios between the SMBHs and the spheroidal bulge are an order of magnitude higher than the mean value found in local spheroidal galaxies, which is in agreement with results from previous CO observations of high redshift quasars.

We show that the mass fraction of GMC gas (n>100 cm^-3) in dense (n>>10^4 cm^-3) star-forming clumps, observable in dense molecular tracers (L_HCN/L_CO(1-0)), is a sensitive probe of the strength and mechanism(s) of stellar feedback. Using high-resolution galaxy-scale simulations with pc-scale resolution and explicit models for feedback from radiation pressure, photoionization heating, stellar winds, and supernovae (SNe), we make predictions for the dense molecular gas tracers as a function of GMC and galaxy properties and the efficiency of stellar feedback. In models with weak/no feedback, much of the mass in GMCs collapses into dense sub-units, predicting L_HCN/L_CO(1-0) ratios order-of-magnitude larger than observed. By contrast, models with feedback properties taken directly from stellar evolution calculations predict dense gas tracers in good agreement with observations. Changing the strength or timing of SNe tends to move systems along, rather than off, the L_HCN-L_CO relation (because SNe heat lower-density material, not the high-density gas). Changing the strength of radiation pressure (which acts efficiently in the highest density gas), however, has a much stronger effect on L_HCN than on L_CO. We predict that the fraction of dense gas (L_HCN/L_CO(1-0)) increases with increasing GMC surface density; this drives a trend in L_HCN/L_CO(1-0) with SFR and luminosity which has tentatively been observed. Our results make specific predictions for enhancements in the dense gas tracers in unusually dense environments such as ULIRGs and galactic nuclei (including the galactic center).

We use simulations with realistic models for stellar feedback to study galaxy mergers. These high resolution (1 pc) simulations follow formation and destruction of individual GMCs and star clusters. The final starburst is dominated by in situ star formation, fueled by gas which flows inwards due to global torques. The resulting high gas density results in rapid star formation. The gas is self gravitating, and forms massive (~10^10 M_sun) GMCs and subsequent super-starclusters (masses up to 10^8 M_sun). However, in contrast to some recent simulations, the bulk of new stars which eventually form the central bulge are not born in superclusters which then sink to the center of the galaxy, because feedback efficiently disperses GMCs after they turn several percent of their mass into stars. Most of the mass that reaches the nucleus does so in the form of gas. The Kennicutt-Schmidt law emerges naturally as a consequence of feedback balancing gravitational collapse, independent of the small-scale star formation microphysics. The same mechanisms that drive this relation in isolated galaxies, in particular radiation pressure from IR photons, extend over seven decades in SFR to regulate star formation in the most extreme starbursts (densities >10^4 M_sun/pc^2). Feedback also drives super-winds with large mass loss rates; but a significant fraction of the wind material falls back onto the disks at later times, leading to higher post-starburst SFRs in the presence of stellar feedback. Strong AGN feedback is required to explain sharp cutoffs in star formation rate. We compare the predicted relic structure, mass profile, morphology, and efficiency of disk survival to simulations which do not explicitly resolve GMCs or feedback. Global galaxy properties are similar, but sub-galactic properties and star formation rates can differ significantly.

We consider the effects of different star formation criteria on galactic scales, in high-resolution simulations with explicitly resolved GMCs and stellar feedback. We compare: (1) a self-gravity criterion (based on the local virial parameter and the assumption that self-gravitating gas collapses to high density in a free-fall time), (2) a fixed density threshold, (3) a molecular-gas law, (4) a temperature threshold, (5) a Jeans-instability requirement, (6) a criteria that cooling times be shorter than dynamical times, and (7) a convergent-flow criterion. We consider these both MW-like and high-density (starburst) galaxies. With feedback present, all models produce identical integrated star formation rates (SFRs), in agreement with the Kennicutt relation. Without feedback all produce orders-of-magnitude excessive SFRs. This is totally dependent on feedback and independent of the SF law. However, the spatial and density distribution of SF depend strongly on the SF criteria. Because cooling rates are generally fast and gas is turbulent, criteria (4)-(7) are 'weak' and spread SF uniformly over the disk (above densities n~0.01-0.1 cm^-3). A molecular criterion (3) localizes to higher densities, but still a wide range; for Z Z_solar, it is similar to a density threshold at n~1 cm^-3 (well below mean densities in the MW center or starbursts). Fixed density thresholds (2) can always select the highest densities, but must be adjusted for simulation resolution and galaxy properties; the same threshold that works in a MW-like simulation will select nearly all gas in a starburst. Binding criteria (1) tend to adaptively select the largest over-densities, independent of galaxy model or resolution, and automatically predict clustered SF. We argue that this SF model is most physically-motivated and presents significant numerical advantages in large-dynamic range simulations.

We study galaxy super-winds driven in major mergers, using pc-resolution simulations with detailed models for stellar feedback that can self-consistently follow the formation/destruction of GMCs and generation of winds. The models include molecular cooling, star formation at high densities in GMCs, and gas recycling and feedback from SNe (I&II), stellar winds, and radiation pressure. We study mergers of systems from SMC-like dwarfs and Milky Way analogues to z~2 starburst disks. Multi-phase super-winds are generated in all passages, with outflow rates up to ~1000 M_sun/yr. However, the wind mass-loading efficiency (outflow rate divided by SFR) is similar to that in isolated galaxy counterparts of each merger: it depends more on global galaxy properties (mass, size, escape velocity) than on the dynamical state of the merger. Winds tend to be bi- or uni-polar, but multiple 'events' build up complex morphologies with overlapping, differently-oriented bubbles/shells at a range of radii. The winds have complex velocity and phase structure, with material at a range of speeds up to ~1000 km/s, and a mix of molecular, ionized, and hot gas that depends on galaxy properties and different feedback mechanisms. These simulations resolve a problem in some 'sub-grid' models, where simple wind prescriptions can dramatically suppress merger-induced starbursts. But despite large mass-loading factors (>~10) in the winds, the peak SFRs are comparable to those in 'no wind' simulations. Wind acceleration does not act equally, so cold dense gas can still lose angular momentum and form stars, while blowing out gas that would not have participated in the starburst in the first place. Considerable wind material is not unbound, and falls back on the disk at later times post-merger, leading to higher post-starburst SFRs in the presence of stellar feedback. This may require AGN feedback to explain galaxy quenching.

Panchromatic observations of the best candidate HyLIRG from the widest Herschel extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z=2.41 galaxies across a ~100-kpc region - a cluster of starbursting proto-ellipticals. Via sub-arcsecond interferometric imaging we have measured accurate gas and star-formation surface densities. The two brightest galaxies span ~3 kpc FWHM in submm/radio continuum and CO J=4-3, and double that in CO J=1-0. The broad CO line is due partly to the multitude of constituent galaxies and partly to large rotational velocities in two counter-rotating gas disks -- a scenario predicted to lead to the most intense starbursts, which will therefore come in pairs. The disks have M(dyn) of several x 10^11 Msun, and gas fractions of ~40%. Velocity dispersions are modest so the disks are unstable, potentially on scales commensurate with their radii: these galaxies are undergoing extreme bursts of star formation, not confined to their nuclei, at close to the Eddington limit. Their specific star-formation rates place them ~>5x above the main sequence, which supposedly comprises large gas disks like these. Their high star-formation efficiencies are difficult to reconcile with a simple volumetric star-formation law. N-body and dark matter simulations suggest this system is the progenitor of a B(inary)-type ~10^14.6-Msun cluster.

I present a model for the star formation properties of z~2 starburst galaxies. Here, I discuss models for the formation of high-z Submillimeter Galaxies, as well as the CO-H2 conversion factor for these systems. I then apply these models to literature observations. I show that when using a functional form for XCO that varies smoothly with the physical properties in galaxies, galaxies at both local and high-z lie on a unimodal Kennicutt-Schmidt star formation law, with power-law index of ~2. The inferred gas fractions of these galaxies are large (fgas ~ 0.2-0.4), though a factor ~2 lower than most literature estimates that utilize locally-calibrated CO-H2 conversion factors.

We have been carrying out a systematic survey of the star formation and ISM properties in the host galaxies of z~6 quasars. Our 250 GHz observations, together with available data from the literature, yield a sample of 14 z~6 quasars that are bright in millimeter dust continuum emission with estimated FIR luminosities of a few 10^12 to 10^13 Lsun. Most of these millimeter-detected z~6 quasars have also been detected in molecular CO line emission, indicating molecular gas masses on order of 10^10 Msun. We have searched for [C II] 158 micron fine structure line emission toward four of the millimeter bright z~6 quasars with ALMA and all of them have been detected. All these results suggest massive star formation at rates of about 600 to 2000 Msun/yr over the central few kpc region of these quasar host galaxies.

High-redshift submillimetre galaxies (SMGs) are some of the most rapidly star-forming galaxies in the Universe. Historically, galaxy formation models have had difficulty explaining the observed number counts of SMGs. We combine a semi-empirical model with 3-D hydrodynamical simulations and 3-D dust radiative transfer to predict the number counts of unlensed SMGs. Because the stellar mass functions, gas and dust masses, and sizes of our galaxies are constrained to match observations, we can isolate uncertainties related to the dynamical evolution of galaxy mergers and the dust radiative transfer. The number counts and redshift distributions predicted by our model agree well with observations. Isolated disc galaxies dominate the faint (S_1.1 < ~1 mJy, or S_850 < ~2 mJy) population. The brighter sources are a mix of merger-induced starbursts and galaxy-pair SMGs; the latter subpopulation accounts for ~30 - 50 per cent of all SMGs at all S_1.1 > ~0.5 mJy (S_850 > ~1 mJy). The mean redshifts are ~3.0 - 3.5, depending on the flux cut, and the brightest sources tend to be at higher redshifts. Because the galaxy-pair SMGs will be resolved into multiple fainter sources by ALMA, the bright ALMA counts should be as much as 2 times less than those observed using single-dish telescopes. The agreement between our model, which uses a Kroupa IMF, and observations suggests that the IMF in high-redshifts starbursts need not be top-heavy; if the IMF were top-heavy, our model would over-predict the number counts. We conclude that the difficulty some models have reproducing the observed SMG counts is likely indicative of more general problems -- such as an under-prediction of the abundance of massive galaxies or a star formation rate--stellar mass relation normalisation lower than that observed -- rather than a problem specific to the SMG population.

CO measurements of z~1-4 galaxies have found that their baryonic gas fractions are significantly higher than galaxies at z=0, with values ranging from 20-80 %. Here, we suggest that the gas fractions inferred from observations of star-forming galaxies at high-z are overestimated, owing to the adoption of locally-calibrated CO-H2 conversion factors (Xco). Evidence from both observations and numerical models suggest that Xco varies smoothly with the physical properties of galaxies, and that Xco can be parameterised simply as a function of both gas phase metallicity and observed CO surface brightness. When applying this functional form, we find fgas ~10-40 % in galaxies with M*=10^10-10^12 Msun at high-z. Moreover, the scatter in the observed fgas-M* relation is lowered by a factor of two. The lower inferred gas fractions arise physically because the interstellar media of high-z galaxies have higher velocity dispersions and gas temperatures than their local counterparts, which results in an Xco that is lower than the z=0 value for both quiescent discs and starbursts. We further compare these gas fractions to those predicted by cosmological galaxy formation models. We show that while the canonically inferred gas fractions from observations are a factor of 2-3 larger at a given stellar mass than predicted by models, our rederived Xco values for z=1-4 galaxies results in revised gas fractions that agree significantly better with the simulations.

It is well-established observationally that the characteristic angular momentum axis on small scales around AGN, traced by radio jets and the putative torus, is not well-correlated with the large-scale angular momentum axis of the host galaxy. In this paper, we show that such misalignments arise naturally in high-resolution simulations in which we follow angular momentum transport and inflows from galaxy to sub-pc scales near AGN, triggered either during galaxy mergers or by instabilities in isolated disks. Sudden misalignments can sometimes be caused by single massive clumps falling into the center slightly off-axis, but more generally, they arise even when the gas inflows are smooth and trace only global gravitational instabilities. When several nested, self-gravitating modes are present, the inner ones can precess and tumble in the potential of the outer modes. Resonant angular momentum exchange can flip or re-align the spin of an inner mode on a short timescale, even without the presence of massive clumps. We therefore do not expect that AGN and their host galaxies will be preferentially aligned, nor should the relative alignment be an indicator of the AGN fueling mechanism. We discuss implications of this conclusion for AGN feedback and BH spin evolution. The misalignments may mean that even BHs accreting from smooth large-scale disks will not be spun up to maximal rotation, and so have more modest radiative efficiencies and inefficient jet formation. Even more random orientations are possible if there is further, un-resolved clumpiness in the gas, and more ordered accretion may occur if the inflow is slower and not self-gravitating.

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.

Observations of star-forming galaxies at high-z have suggested discrepancies in the inferred star formation rates (SFRs) either between data and models, or between complementary measures of the SFR. These putative discrepancies could all be alleviated if the stellar IMF is systematically weighted toward more high-mass star formation in rapidly star-forming galaxies. Here, we explore how the IMF might vary under the central assumption that the turnover mass in the IMF, Mc, scales with the Jeans mass in giant molecular clouds (GMCs), M_J. We employ hydrodynamic and radiative transfer simulations of galaxies to predict how the typical GMC Jeans mass, and hence the IMF, varies with galaxy property. We then study the impact of such an IMF on the star formation law, the SFR-M* relation, submillimetre galaxies (SMGs), and the cosmic SFR density. Our main results are: The H2 mass-weighted Jeans mass in a galaxy scales with the SFR when the SFR is greater a few M_sun/yr. SPS modeling shows that this results in a nonlinear relation between SFR and Lbol, such that SFR Lbol^0.88. Using this model relation, the inferred SFR of local ULIRGs decreases by ~2, and that of high-z SMGs decreases by ~3-5. At z 2, this results in a lowered normalisation of the SFR-M* relation in better agreement with models, a reduced discrepancy between the observed cosmic SFR density and stellar mass density evolution, and SMG SFRs that are easier to accommodate in current hierarchical structure formation models. It further results in a Schmidt relation with slope of ~1.6 when utilising a physically motivated form for the CO-H2 conversion factor. While each of the discrepancies considered here could be alleviated without appealing to a varying IMF, the modest variation implied by assuming Mc M_J is a plausible solution that simultaneously addresses numerous thorny issues regarding the SFRs of high-z galaxies.

The most common means of converting an observed CO line intensity into a molecular gas mass requires the use of a conversion factor (Xco). While in the Milky Way this quantity does not appear to vary significantly, there is good reason to believe that Xco will depend on the larger-scale galactic environment. Utilising numerical models, we investigate how varying metallicities, gas temperatures and velocity dispersions in galaxies impact the way CO line emission traces the underlying H2 gas mass, and under what circumstances Xco may differ from the Galactic mean value. We find that, due to the combined effects of increased gas temperature and velocity dispersion, Xco is depressed below the Galactic mean in high surface density environments such as ULIRGs. In contrast, in low metallicity environments, Xco tends to be higher than in the Milky Way, due to photodissociation of CO in metal-poor clouds. At higher redshifts, gas-rich discs may have gravitationally unstable clumps which are warm (due to increased star formation) and have elevated velocity dispersions. These discs tend to have Xco values ranging between present-epoch gas-rich mergers and quiescent discs at low-z. This model shows that on average, mergers do have lower Xco values than disc galaxies, though there is significant overlap. Xco varies smoothly with the local conditions within a galaxy, and is not a function of global galaxy morphology. We combine our results to provide a general fitting formula for Xco as a function of CO line intensity and metallicity. We show that replacing the traditional approach of using one constant Xco for starbursts and another for discs with our best-fit function produces star formation laws that are continuous rather than bimodal, and that have significantly reduced scatter.

The CO-H2 conversion factor in galaxies is typically described as bimodal: one value for discs and quiescent regions, and another (lower) value for mergers and starbursts. In this proceeding, I will describe both empirical observational evidence that the conversion factor varies with physical environment, as well as a theoretical model which aims to understand the physical processes which drive these variations. I present a functional form for Xco which can be applied to observations ranging in scale from ~70 pc to galaxy-wide scales, and show the consequences of the application of this model to the Kennicutt-Schmidt star formation law.

Multi-scale simulations have made it possible to follow gas inflows onto massive black holes (BHs) from galactic scales to the accretion disk. When sufficient gas is driven towards the BH, gravitational instabilities generically form lopsided, eccentric disks that propagate inwards. The lopsided stellar disk exerts a strong torque on the gas disk, driving inflows that fuel rapid BH growth. Here, we investigate whether the same gas disk is the 'torus' invoked to explain obscured AGN. The disk is generically thick and has characteristic ~1-10 pc sizes and masses resembling those required of the torus. The scale heights and obscured fractions of the predicted torii are substantial even in the absence of strong stellar feedback providing the vertical support. Rather, they can be maintained by strong bending modes and warps excited by the inflow-generating instabilities. Other properties commonly attributed to feedback processes may be explained by dynamical effects: misalignment between torus and host galaxy, correlations between local SFR and turbulent gas velocities, and dependence of obscured fractions on AGN luminosity or SFR. We compare the predicted torus properties with observations of gas surface density profiles, kinematics, scale heights, and SFR densities in AGN nuclei, and find that they are consistent. We argue that it is not possible to reproduce these observations and the observed column density (N_H) distribution without a clumpy gas distribution, but allowing for clumping on small scales the predicted N_H distribution is in good agreement with observations from 10^20-27 cm^-2. We examine how N_H scales with galaxy and AGN properties, and find that AGN feedback may be necessary to explain some trends with luminosity and/or redshift. The torus is not merely a bystander or passive fuel source for accretion, but is itself the mechanism driving accretion.

We perform 3-D dust radiative transfer (RT) calculations on hydrodynamic simulations of isolated and merging disk galaxies in order to quantitatively study the dependence of observed-frame submillimeter (submm) flux density on galaxy properties. We find that submm flux density and star formation rate (SFR) are related in dramatically different ways for quiescently star-forming galaxies and starbursts. Because the stars formed in the merger-induced starburst do not dominate the bolometric luminosity and the rapid drop in dust mass and more compact geometry cause a sharp increase in dust temperature during the burst, starbursts are very inefficient at boosting submm flux density (e.g., a $\ga16$x boost in SFR yields a $\la 2$x boost in submm flux density). Moreover, the ratio of submm flux density to SFR differs significantly between the two modes; thus one cannot assume that the galaxies with highest submm flux density are necessarily those with the highest bolometric luminosity or SFR. These results have important consequences for the bright submillimeter-selected galaxy (SMG) population. Among them are: 1. The SMG population is heterogeneous. In addition to merger-driven starbursts, there is a subpopulation of galaxy pairs, where two disks undergoing a major merger but not yet strongly interacting are blended into one submm source because of the large ($\ga 15$", or $\sim 130$ kpc at $z = 2$) beam of single-dish submm telescopes. 2. SMGs must be very massive ($M_{\star} \ga 6 \times 10^{10} \msun$). 3. The infall phase makes the SMG duty cycle a factor of a few greater than what is expected for a merger-driven starburst. (Abridged.)

Relating the observed CO emission from giant molecular clouds (GMCs) to the underlying H2 column density is a long-standing problem in astrophysics. While the Galactic CO-H2 conversion factor (Xco) appears to be reasonably constant, observations indicate that Xco may be depressed in high-surface density starburst environments. Using a multi-scale approach, we investigate the dependence of Xco on the galactic environment in numerical simulations of disc galaxies and galaxy mergers. Xco is proportional to the GMC surface density divided by the integrated CO intensity, Wco, and Wco is related to the kinetic temperature and velocity dispersion in the cloud. In disc galaxies (except within the central ~ kpc), the galactic environment is largely unimportant in setting the physical properties of GMCs provided they are gravitationally bound. The temperatures are roughly constant at ~10 K due to the balance of CO cooling and cosmic ray heating, giving a nearly constant CO-H2 conversion factor in discs. In mergers, the velocity dispersion of the gas rises dramatically during coalescence. The gas temperature also rises as it couples well to the warm (~50 K) dust at high densities (n > 10^4 cm^-3). The rise in velocity dispersion and temperature combine to offset the rise in surface density in mergers, causing Xco to drop by a factor of ~2-10 compared to the disc simulation. This model predicts that high-resolution ALMA observations of nearby ULIRGs should show velocity dispersions of ~10-100 km/s, and brightness temperatures comparable to the dust temperatures.

Recent observations of excited CO emission lines from z~2 disc galaxies have shed light on the Kennicutt-Schmidt relation at high-z via observed SFR-CO (J=2-1) and (J=3-2) relations. Here, we describe a novel methodology for utilising these observations of high-excitation CO to derive the underlying Schmidt (SFR-rho^N) relationship. To do this requires an understanding of the potential effects of differential CO excitation with SFR. If the most heavily star-forming galaxies have a larger fraction of their gas in highly excited CO states than the lower SFR galaxies, then the observed molecular SFR-CO^alpha index, alpha, will be less than the underlying (volumetric) Schmidt index, N. Utilising a combination of SPH models of galaxy evolution and molecular line radiative transfer, we present the first calculations of CO excitation in z~2 disc galaxies with the aim of developing a mapping between various observed SFR-CO relationships and the underlying Schmidt relation. We find that even in relatively luminous z~2 discs, differential excitation does indeed exist, resulting in alpha < N for highly excited CO lines. This means that an observed (e.g.) SFR-CO (J=3-2) relation does not map linearly to SFR-H2 relation. We utilise our model results to provide a mapping from alpha to N for the range of Schmidt indices N=1-2. By comparing to recent observational surveys, we find that the observed SFR-CO (J=2-1) and SFR-CO (J=3-2) relations suggest that an underlying SFR rho^1.5 relation describes z~2 disc galaxies.

Theoretical models have had difficulty matching the observed number density of sub-millimeter galaxies (SMGs), causing some authors (e.g., Baugh et al. 2005) to suggest that SMGs provide evidence for a top-heavy initial mass function (IMF). To test this claim, we have, for the first time, combined high-resolution 3-D hydrodynamic simulations of isolated and merging massive, gas-rich galaxies, radiative transfer, and a semi-empirical merger rate model to predict the number density of SMGs. Our model can reproduce the observed SMG number density even when using a standard (Kroupa) IMF. The agreement is due to a combination of relatively long sub-mm duty cycles for mergers (a few times 10^8 years for our most massive models), which owe to our combination of high-resolution 3-D hydrodynamic simulations and dust radiative transfer; sufficient number densities of massive, gas-rich mergers; and the decrease in sub-mm counts observed by recent deep/wide surveys (e.g., Austermann et al. 2010) relative to previous surveys. Our results suggest that the observed SMG number counts do not provide evidence for a top-heavy IMF at high redshift.

We study the simulated CO emission from elliptical galaxies formed in the mergers of gas-rich disk galaxies. The cold gas not consumed in the merger-driven starburst quickly resettles into a disk-like configuration. By analyzing a variety of arbitrary merger orbits that produce a range of fast to slow-rotating remnants, we find that molecular disk formation is a fairly common consequence of gas-rich galaxy mergers. Hence, if a molecular disk is observed in an early-type merger remnant, it is likely the result of a "wet merger" rather than a "dry merger". We compare the physical properties from our simulated disks (e.g. size and mass) and find reasonably good agreement with recent observations. Finally, we discuss the detectability of these disks as an aid to future observations.

We present a physical model for the origin of z~2 Dust-Obscured Galaxies (DOGs), a class of high-redshift ULIRGs selected at 24 micron which are particularly optically faint (24/R>1000). By combining N-body/SPH simulations of high redshift galaxy evolution with 3D polychromatic dust radiative transfer models, we find that luminous DOGs (with F24 > 0.3 mJy at z~2 are well-modeled as extreme gas-rich mergers in massive (~5x10^12-10^13 Msun) halos, with elevated star formation rates (~500-1000 Msun/yr) and/or significant AGN growth (Mdot > 0.5 Msun/yr), whereas less luminous DOGs are more diverse in nature. At final coalescence, merger-driven DOGs transition from being starburst dominated to AGN dominated, evolving from a "bump" to a power-law shaped mid-IR (IRAC) spectral energy distribution (SED). After the DOG phase, the galaxy settles back to exhibiting a "bump" SED with bluer colors and lower star formation rates. While canonically power-law galaxies are associated with being AGN-dominated, we find that the power-law mid-IR SED can owe both to direct AGN contribution, as well as to a heavily dust obscured stellar bump at times that the galaxy is starburst dominated. Thus power-law galaxies can be either starburst or AGN dominated. Less luminous DOGs can be well-represented either by mergers, or by massive ($M_{\rm baryon} ~5x10^11 Msun) secularly evolving gas-rich disc galaxies (with SFR > 50 Msun/yr). By utilising similar models as those employed in the SMG formation study of Narayanan et al. (2010), we investigate the connection between DOGs and SMGs. We find that the most heavily star-forming merger driven DOGs can be selected as Submillimetre Galaxies (SMGs), while both merger-driven and secularly evolving DOGs typically satisfy the BzK selection criteria.

We present a detailed analysis of the relation between infrared luminosity and molecular line luminosity, for a variety of molecular transitions, using a sample of 34 nearby galaxies spanning a broad range of infrared luminosities (10^{10} < L_{IR} < 10^{12.5} L_sun). We show that the power-law index of the relation is sensitive to the critical density of the molecular gas tracer used, and that the dominant driver in observed molecular line ratios in galaxies is the gas density. As most nearby ultraluminous infrared galaxies (ULIRGs) exhibit strong signatures of active galactic nuclei (AGN) in their center, we revisit previous claims questioning the reliability of HCN as a probe of the dense gas responsible for star formation in the presence of AGN. We find that the enhanced HCN(1-0)/CO(1-0) luminosity ratio observed in ULIRGs can be successfully reproduced using numerical models with fixed chemical abundances and without AGN-induced chemistry effects. We extend this analysis to a total of ten molecular line ratios by combining the following transitions: CO(1-0), HCO+(1-0), HCO+(3-2), HCN(1-0), and HCN(3-2). Our results suggest that AGNs reside in systems with higher dense gas fraction, and that chemistry or other effects associated with their hard radiation field may not dominate (NGC 1068 is one exception). Galaxy merger could be the underlying cause of increased dense gas fraction and the evolutionary stage of such mergers may be another determinant of the HCN/CO luminosity ratio.

We use a novel method to predict the contribution of normal star-forming galaxies, merger-induced bursts, and obscured AGN, to IR luminosity functions (LFs) and global SFR densities. We use empirical halo occupation constraints to populate halos with galaxies and determine the distribution of normal and merging galaxies. Each system can then be associated with high-resolution hydrodynamic simulations. We predict the distribution of observed luminosities and SFRs, from different galaxy classes, as a function of redshift from z=0-6. We provide fitting functions for the predicted LFs, quantify the uncertainties, and compare with observations. At all redshifts, 'normal' galaxies dominate the LF at moderate luminosities ~L* (the 'knee'). Merger-induced bursts increasingly dominate at L>>L*; at the most extreme luminosities, AGN are important. However, all populations increase in luminosity at higher redshifts, owing to increasing gas fractions. Thus the 'transition' between normal and merger-dominated sources increases from the LIRG-ULIRG threshold at z~0 to bright Hyper-LIRG thresholds at z~2. The transition to dominance by obscured AGN evolves similarly, at factor of several higher L_IR. At all redshifts, non-merging systems dominate the total luminosity/SFR density, with merger-induced bursts constituting ~5-10% and AGN ~1-5%. Bursts contribute little to scatter in the SFR-stellar mass relation. In fact, many systems identified as 'ongoing' mergers will be forming stars in their 'normal' (non-burst) mode. Counting this as 'merger-induced' star formation leads to a stronger apparent redshift evolution in the contribution of mergers to the SFR density.

We describe a model for the formation of \zsim 2 Submillimeter Galaxies (SMGs) which simultaneously accounts for both average and bright SMGs while providing a reasonable match to their mean observed spectral energy distributions (SEDs). By coupling hydrodynamic simulations of galaxy mergers with the high resolution 3D polychromatic radiative transfer code Sunrise, we find that a mass sequence of merger models which use observational constraints as physical input naturally yield objects which exhibit black hole, bulge, and H2 gas masses similar to those observed in SMGs. The dominant drivers behind the 850 micron flux are the masses of the merging galaxies and the stellar birthcloud covering fraction. The most luminous (S850 ~ 15 mJy) sources are recovered by ~10^13 Msun 1:1 major mergers with a birthcloud covering fraction close to unity, whereas more average SMGs ~5-7 mJy) may be formed in lower mass halos ~5x10^12 Msun. These models demonstrate the need for high spatial resolution hydrodynamic and radiative transfer simulations in matching both the most luminous sources as well as the full SEDs of SMGs. While these models suggest a natural formation mechanism for SMGs, they do not attempt to match cosmological statistics of galaxy populations; future efforts along this line will help ascertain the robustness of these models.

We present a model for the CO molecular line emission from high redshift Submillimeter Galaxies (SMGs). By combining hydrodynamic simulations of gas rich galaxy mergers with the polychromatic radiative transfer code, Sunrise, and the 3D non-LTE molecular line radiative transfer code, Turtlebeach, we show that if SMGs are typically a transient phase of major mergers, their observed compact CO spatial extents, broad line widths, and high excitation conditions (CO SED) are naturally explained. In this sense, SMGs can be understood as scaled-up analogs to local ULIRGs. We utilize these models to investigate the usage of CO as an indicator of physical conditions. We find that care must be taken when applying standard techniques. The usage of CO line widths as a dynamical mass estimator from SMGs can possibly overestimate the true enclosed mass by a factor ~1.5-2. At the same time, assumptions of line ratios of unity from CO J=3-2 (and higher lying lines) to CO (J=1-0) will oftentimes lead to underestimates of the inferred gas mass. We provide tests for these models by outlining predictions for experiments which are imminently feasible with the current generation of bolometer arrays and radio-wave spectrometers.

We present a merger-driven evolutionary model for the production of luminous (LIRGs) and ultraluminous infrared galaxies (ULIRGs) with warm IR colours. Our results show that simulations of gas--rich major mergers including star formation, black hole growth, and feedback can produce warm (U)LIRGs. We also find that while the warm evolutionary phase is associated with increased AGN activity, star formation alone may be sufficient to produce warm IR colours. However, the transition can be suppressed entirely - even when there is a significant AGN contribution - when we assume a single-phase ISM, which maximizes the attenuation. Finally, our evolutionary models are consistent with the 25-to-60 micron flux density ratio vs. L_HX/L_IR relation for local LIRGs and ULIRGs, and predict the observed scatter in IR color at fixed L_HX/L_IR. Therefore, our models suggest a cautionary note in the interpretation of warm IR colours: while associated with periods of active black hole growth, they are probably produced by a complex mix of star formation and AGN activity intermediate between the cold star formation dominated phase and the birth of a bright, unobscured quasar.

When, and how, did the first galaxies and supermassive black holes (SMBH) form, and how did they reionization the Universe? First galaxy formation and cosmic reionization are among the last frontiers in studies of cosmic structure formation. We delineate the detailed astrophysical probes of early galaxy and SMBH formation afforded by observations at centimeter through submillimeter wavelengths. These observations include studies of the molecular gas (= the fuel for star formation in galaxies), atomic fine structure lines (= the dominant ISM gas coolant), thermal dust continuum emission (= an ideal star formation rate estimator), and radio continuum emission from star formation and relativistic jets. High resolution spectroscopic imaging can be used to study galaxy dynamics and star formation on sub-kpc scales. These cm and mm observations are the necessary compliment to near-IR observations, which probe the stars and ionized gas, and X-ray observations, which reveal the AGN. Together, a suite of revolutionary observatories planned for the next decade from centimeter to X-ray wavelengths will provide the requisite panchromatic view of the complex processes involved in the formation of the first generation of galaxies and SMBHs, and cosmic reionization.

We present observations made with the 10m Heinrich Hertz Sub-Millimeter Telescope of HCN (3-2) emission from a sample of 30 nearby galaxies ranging in infrared luminosity from 10^10 - 10^12.5 L_sun and HCN (3-2) luminosity from 10^6 - 10^9 K km s^-1 pc^2. We examine the correlation between the infrared luminosity and HCN (3-2) luminosity and find that the best fit linear regression has a slope (in log-log space) of 0.74+/-0.12. Including recently published data from Gracia-Carpio et al. tightens the constraints on the best-fit slope to 0.79+/-0.09. This slope below unity suggests that the HCN (3-2) molecular line luminosity is not linearly tracing the amount of dense gas. Our results are consistent with predictions from recent theoretical models that find slopes below unity when the line luminosity depends upon the average gas density with a power-law index greater than a Kennicutt-Schmidt index of 1.5.

We investigate the relationship between the star formation rate (SFR) and dense molecular gas mass in the nuclei of galaxies. To do this, we utilize the observed 850 micron luminosity as a proxy for the infrared luminosity and SFR, and correlate this with the observed CO (J=3-2) luminosity. We find tentative evidence that the LIR-CO (J=3-2) index is similar to the Kennicutt-Schmidt (KS) index (N ~ 1.5) in the central ~1.7 kpc of galaxies, and flattens to a roughly linear index when including emission from the entire galaxy. This result may imply that the volumetric Schmidt relation is the underlying driver behind the observed SFR-dense gas correlations, and provides tentative confirmation for recent numerical models. While the data exclude the possibility of a constant LIR-CO (J=3-2) index for both galaxy nuclei and global measurements at the ~80% confidence level, the considerable error bars cannot preclude alternative interpretations.

We derive a physical model for the observed relations between star formation rate (SFR) and molecular line (CO and HCN) emission in galaxies, and show how these observed relations are reflective of the underlying star formation law. We do this by combining 3D non-LTE radiative transfer calculations with hydrodynamic simulations of isolated disk galaxies and galaxy mergers. We demonstrate that the observed SFR-molecular line relations are driven by the relationship between molecular line emission and gas density, and anchored by the index of the underlying Schmidt law controlling the SFR in the galaxy. Lines with low critical densities (e.g. CO J=1-0) are typically thermalized and trace the gas density faithfully. In these cases, the SFR will be related to line luminosity with an index similar to the Schmidt law index. Lines with high critical densities greater than the mean density of most of the emitting clouds in a galaxy (e.g. CO J=3-2, HCN J=1-0) will have only a small amount of thermalized gas, and consequently a superlinear relationship between molecular line luminosity and mean gas density. This results in a SFR-line luminosity index less than the Schmidt index for high critical density tracers. One observational consequence of this is a significant redistribution of light from the small pockets of dense, thermalized gas to diffuse gas along the line of sight, and prodigious emission from subthermally excited gas. At the highest star formation rates, the SFR-Lmol slope tends to the Schmidt index, regardless of the molecular transition. The fundamental relation is the Kennicutt-Schmidt law, rather than the relation between SFR and molecular line luminosity. We use these results to make imminently testable predictions for the SFR-molecular line relations of unobserved transitions.

We assess the impact of starburst and AGN feedback-driven winds on the CO emission from galaxy mergers, and, in particular, search for signatures of these winds in the simulated CO morphologies and emission line profiles. We do so by combining a 3D non-LTE molecular line radiative transfer code with smoothed particle hydrodynamics (SPH) simulations of galaxy mergers that include prescriptions for star formation, black hole growth, a multiphase interstellar medium (ISM), and the winds associated with star formation and black hole growth. Our main results are: (1) Galactic winds can drive outflows of masses ~10^8-10^9 Msun which may be imaged via CO emission line mapping. (2) AGN feedback-driven winds are able to drive imageable CO outflows for longer periods of time than starburst-driven winds owing to the greater amount of energy imparted to the ISM by AGN feedback compared to star formation. (3) Galactic winds can control the spatial extent of the CO emission in post-merger galaxies, and may serve as a physical motivation for the sub-kiloparsec scale CO emission radii observed in local advanced mergers. (4) Secondary emission peaks at velocities greater than the circular velocity are seen in the CO emission lines in all models. In models with winds, these high velocity peaks are seen to preferentially correspond to outflowing gas entrained in winds, which is not the case in the model without winds. The high velocity peaks seen in models without winds are typically confined to velocity offsets (from the systemic) < 1.7 times the circular velocity, whereas the models with AGN feedback-driven winds can drive high velocity peaks to ~2.5 times the circular velocity.

We investigate the nature of CO emission from z~6 quasars by combining non-LTE radiative transfer calculations with merger-driven models of z~6 quasar formation that arise naturally in LCDM cosmological simulations. We consider four model quasars formed in 10^12-10^13 M_sun halos from different merging histories. Our main results follow. Owing to massive starbursts and funneling of dense gas into the nuclear regions of merging galaxies, the CO is highly excited and the flux density peaks between J=5-8. The CO morphology of z~6 quasars often exhibits multiple emission peaks which arise from H2 concentrations which have not yet fully coalesced. Quasars at z~6 display a large range of sightline dependent line widths such that the lines are narrowest when the rotating H2 gas associated with the quasar is viewed face-on (when L_B is largest), and broadest when the gas is seen edge-on (when L_B is lowest). Thus for all models selection effects exist such that quasars selected for optical luminosity are preferentially face-on which may result in detected CO line widths narrower than the median. The sightline averaged line width is reflective of the circular velocity (V_c) of the host halo, and ranges from sigma~300-650 km/s. For optically selected QSOs, 10-25% (halo-mass dependant) of sightlines have narrow line widths compatible with the sole CO detection at z>6, J1148+5251. When accounting for both the temporal evolution of CO line widths, as well as the redshift evolution of halo V_c, these models self-consistently account for the CO line widths of both z~2 sub-mm galaxies and QSO's. Finally, the dynamical mass derived from the sightline averaged line widths provides a good estimate of the total mass, and allows for a stellar bulge and SMBH consistent with the local M_BH-M_bulge relation. [abridged]

In the era of high resolution submillimeter interferometers, it will soon be possible to observe the neutral circumstellar medium directly involved in gas giant planet (GGP) formation at physical scales previously unattainable. In order to explore possible signatures of gas giant planet formation via disk instabilities, we have combined a 3D, non-local thermodynamic equilibrium (LTE) radiative transfer code with a 3D, finite differences hydrodynamical code to model molecular emission lines from the vicinity of a 1.4 M_J self-gravitating proto-GGP. Here, we explore the properties of rotational transitions of the commonly observed dense gas tracer, HCO+. Our main results are the following: 1. Very high lying HCO+ transitions (e.g. HCO+ J=7-6) can trace dense planet forming clumps around circumstellar disks. Depending on the molecular abundance, the proto-GGP may be directly imageable by the Atacama Large Millimeter Array (ALMA). 2. HCO+ emission lines are heavily self-absorbed through the proto-GGP's dense molecular core. This signature is nearly ubiquitous, and only weakly dependent on assumed HCO+ abundances. The self-absorption features are most pronounced at higher angular resolutions. Dense clumps that are not self-gravitating only show minor self-absorption features. 3. Line temperatures are highest through the proto-GGP at all assumed abundances and inclination angles. Conversely, due to self-absorption in the line, the velocity-integrated intensity may not be. High angular resolution interferometers such as the Submillimeter Array (SMA) and ALMA may be able to differentiate between competing theories of gas giant planet formation.

We study the effects of feedback from active galactic nuclei (AGN) on emission from molecular gas in galaxy mergers by combining hydrodynamic simulations which include black holes with a three-dimensional, non-local thermodynamic equilibrium (LTE) radiative transfer code. We find that molecular clouds entrained in AGN winds produce an extended CO morphology with significant off-nuclear emission, which may be detectable via contour mapping. Furthermore, kinematic signatures of these molecular outflows are visible in emission line profiles when the outflow has a large line of sight velocity. Our results can help interpret current and upcoming observations of luminous infrared galaxies, as well as provide a detailed test of subresolution prescriptions for supermassive black hole growth in galaxy-scale hydrodynamic simulations.

We investigate the effect of embedded active galactic nuclei (AGN) in galaxy mergers on the CO molecular line emission by combining non-local thermodynamic equilibrium (LTE) radiative transfer calculations with hydrodynamic simulations. We find that AGN feedback energy in gas rich galaxy mergers can contribute to large molecular outflows which may be detectable via velocity-integrated emission contour maps, as well as through kinematic features in the emission line profiles.

The role of star formation in luminous and ultraluminous infrared galaxies is a hotly debated issue: while it is clear that starbursts play a large role in powering the IR luminosity in these galaxies, the relative importance of possible enshrouded AGNs is unknown. It is therefore important to better understand the role of star forming gas in contributing to the infrared luminosity in IR-bright galaxies. The J=3 level of 12CO lies 33K above ground and has a critical density of ~1.5 X 10^4 cm^-3. The 12CO(J=3-2) line serves as an effective tracer for warm-dense molecular gas heated by active star formation. Here we report on 12CO (J=3-2) observations of 17 starburst spirals, LIRGs and ULIRGs which we obtained with the Heinrich Hertz Submillimeter Telescope on Mt. Graham, Arizona. Our main results are the following: 1. We find a nearly linear relation between the infrared luminosity and warm-dense molecular gas such that the infrared luminosity increases as the warm-dense molecular gas to the power 0.92; We interpret this to be roughly consistent with the recent results of Gao & Solomon (2004a,b). 2. We find L_IR/M_H2 ratios ranging from ~10 to ~128 L_sun/M_sun using a standard CO-H2 conversion factor of 3 X 10^20 cm^-2 (K km s^-1)^-1. If this conversion factor is ~an order of magnitude less, as suggested in a recent statistical survey (Yao et al. 2003), then 2-3 of our objects may have significant contributions to the L_IR by dust-enshrouded AGNs.

Quasar spectra have a variety of absorption lines whose origins range from energetic winds expelled from the central engines to unrelated, intergalactic clouds. We present multi-epoch, medium resolution spectra of eight quasars at z~2 that have narrow ``associated'' absorption lines (AALs, within $\pm$5000 km s^{-1} of the emission redshift). Two of these quasars were also known previously to have high-velocity mini-broad absorption lines (mini-BALs). We use these data, spanning ~17 years in the observed frame with two to four observations per object, to search for line strength variations as an identifier of absorption that occurs physically near (``intrinsic'' to) the central AGN. Our main results are the following: Two out of the eight quasars with narrow AALs exhibit variable AAL strengths. Two out of two quasars with high-velocity mini-BALs exhibit variable mini-BAL strengths. We also marginally detect variability in a high-velocity narrow absorption line (NAL) system, blueshifted \~32,900 km s^{-1}$ with respect to the emission lines. No other absorption lines in these quasars appeared to vary. The outflow velocities of the variable AALs are 3140 km s^{-1} and 1490 km s^{-1}. The two mini-BALs identify much higher velocity outflows of ~28,400 km s^{-1} and ~52,000 km s^{-1}. Our temporal sampling yields upper limits on the variation time scales from 0.28 to 6.1 years in the quasar rest frames. The corresponding minimum electron densities in the variable absorbers, based on the recombination time scale, are \~40,000 cm^{-3} to ~1900 cm^{-3}. The maximum distances of the absorbers from the continuum source, assuming photoionization with no spectral shielding, range from ~1.8 kpc to ~7 kpc.