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

Galaxies forming in low-mass halos are thought to be primarily responsible for reionizing the Universe during the first billion years after the Big Bang. Yet, these halos are extremely inefficient at forming stars in the nearby Universe. In this work, we address this apparent tension, and ask whether a physically motivated model of galaxy formation that reproduces the observed abundance of faint galaxies in the nearby Universe is also consistent with available observational constraints on the reionization history. By interfacing the Santa Cruz semi-analytic model for galaxy formation with an analytic reionization model, we constructed a computationally efficient pipeline that connects `ground-level' galaxy formation physics to `top-level' cosmological-scale observables. Based on photometric properties of the galaxy populations predicted up to $z=15$, we compute the reionization history of intergalactic hydrogen. We quantify the three degenerate quantities that influence the total ionizing photon budget, including the abundance of galaxies, the intrinsic production rate of ionizing photons, and the LyC escape fraction. We explore covariances between these quantities using a Markov chain Monte Carlo method. We find that our locally calibrated model is consistent with all currently available constraints on the reionization history, under reasonable assumptions about the LyC escape fraction. We quantify the fraction of ionizing photons contributed by galaxies of different luminosities and find that the galaxies expected to be detected in JWST NIRCam wide and deep surveys are responsible for producing $\sim 40$-$80\%$ of ionizing photons throughout the EoR. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

We present a characterization of the physical properties of a sample of 35 securely-detected, dusty galaxies in the deep ALMA 1.2-mm image obtained as part of the ALMA Spectroscopic Survey in the {\it Hubble} Ultra Deep Field (ASPECS) Large Program. This sample is complemented by 26 additional sources identified via an optical/infrared source positional prior. Using their well-characterized spectral energy distributions, we derive median stellar masses and star formation rates (SFR) of $4.8\times10^{10}~M_\odot$ and 30 $M_\odot$ yr$^{-1}$, and interquartile ranges of $(2.4-11.7)\times10^{10}~M_\odot$ and $20-50~M_\odot$ yr$^{-1}$, respectively. We derive a median spectroscopic redshift of 1.8 with an interquartile range $1.1-2.6$, significantly lower than sub-millimeter galaxies detected in shallower, wide-field surveys. We find that 59\%$\pm$13\%, 6\%$\pm$4\%, and 34\%$\pm$9\% of our sources are within, above and below $\pm0.4$ dex from the SFR- stellar mass relation or main-sequence (MS), respectively. The ASPECS galaxies closely follow the SFR-molecular gas mass relation and other previously established scaling relations, confirming a factor of five increase of the gas-to-stellar mass ratio from $z=0.5$ to $z=2.5$, and a mild evolution of the gas depletion timescales with a typical value of 0.7 Gyr at $z=1-3$. ASPECS galaxies located significantly below the MS, a poorly exploited parameter space, have low gas-to-stellar mass ratios of $\sim0.1-0.2$ and long depletion timescales $>1$ Gyr. Galaxies along the MS dominate the cosmic density of molecular gas at all redshifts. Systems above the main sequence have an increasing contribution to the total gas reservoirs from $z<1$ to $z=2.5$, while the opposite is found for galaxies below the MS.

Measurements of the cosmic far-infrared background (CIB) indicate that emission from many extragalactic phenomena, including star formation and black hole accretion, in the Universe can be obscured by dust. Resolving the CIB to study the population of galaxies in which this activity takes place is a major goal of submillimetre astronomy. Here, we present interferometric 650$\mu$m submillimetre number counts. Using the Band 8 data from the ALMACAL survey, we have analysed 81 ALMA calibrator fields together covering a total area of 5.5~arcmin$^2$. The typical central rms in these fields is $\sim 100 \mu$Jy~beam$^{-1}$ with the deepest maps reaching $\sigma = 47 \mu$Jy~beam$^{-1}$ at sub-arcsec resolution. Multi-wavelength coverage from ALMACAL allows us to exclude contamination from jets associated with the calibrators. However, residual contamination by jets and lensing remain a possibility. Using a signal-to-noise threshold of $4.5\sigma$, we find 21 dusty, star-forming galaxies with 650$\mu$m flux densities of $\geq 0.7 $mJy. At the detection limit we resolve $\simeq 100$ per cent of the CIB at 650$\mu$m, a significant improvement compared to low resolution studies at similar wavelength. We have therefore identified all the sources contributing to the EBL at 650 microns and predict that the contribution from objects with flux 0.7<mJy will be small.

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.

The James Webb Space Telescope (JWST) is expected to enable transformational progress in studying galaxy populations in the very early Universe, during the Epoch of Reionization (EoR). A critical parameter for understanding the sources that reionized the Universe is the Lyman-continuum production efficiency, $\xi_{\rm ion}$, defined as the rate of production of ionizing photons divided by the intrinsic UV luminosity. In this work, we combine self-consistent star formation and chemical enrichment histories predicted by semi-analytic models of galaxy formation with stellar population synthesis (SPS) models to predict the expected dependence of $\xi_{\rm ion}$ on galaxy properties and cosmic epoch from $z=4-10$. We then explore the sensitivity of the production rate of ionizing photons, $\dot{N}_{\rm ion}$, to the choice of SPS model and the treatment of stellar feedback in our galaxy formation model. We compare our results to those of other simulations, constraints from empirical models, and observations. We find that adopting SPS models that include binary stars predict about a factor of two more ionizing radiation than models that only assume single stellar populations. We find that UV-faint, low-mass galaxies have values of $\xi_{\rm ion}$ about 0.25 dex higher than those of more massive galaxies, but find weak evolution with cosmic time, about 0.2 dex from $z\sim 12-4$ at fixed rest-UV luminosity. We provide predictions of $\dot{N}_{\rm ion}$ as a function of $M_{\rm h}$ and a number of other galaxy properties. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

Using the deepest 1.2 mm continuum map to date in the Hubble Ultra Deep Field obtained as part of the ALMA Spectroscopic Survey (ASPECS) large program, we measure the cosmic density of dust and implied gas (H$_{2}+$H I) mass in galaxies as a function of look-back time. We do so by stacking the contribution from all $H$-band selected galaxies above a given stellar mass in distinct redshift bins, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$. At all redshifts, $\rho_{\rm dust}(M_\ast>M,z)$ and $\rho_{\rm gas}(M_\ast>M,z)$ grow rapidly as $M$ decreases down to $10^{10}\,M_\odot$, but this growth slows down towards lower stellar masses. This flattening implies that at our stellar mass-completeness limits ($10^8\,M_\odot$ and $10^{8.9}\,M_\odot$ at $z\sim0.4$ and $z\sim3$), both quantities converge towards the total cosmic dust and gas mass densities in galaxies. The cosmic dust and gas mass densities increase at early cosmic time, peak around $z\sim2$, and decrease by a factor $\sim4$ and 7, compared to the density of dust and molecular gas in the local universe, respectively. The contribution of quiescent galaxies -- i.e., with little on-going star-formation -- to the cosmic dust and gas mass densities is minor ($\lesssim10\%$). The redshift evolution of the cosmic gas mass density resembles that of the star-formation rate density, as previously found by CO-based measurements. This confirms that galaxies have relatively constant star-formation efficiencies (within a factor $\sim2$) across cosmic time. Our results also imply that by $z\sim0$, a large fraction ($\sim90\%$) of dust formed in galaxies across cosmic time has been destroyed or ejected to the intergalactic medium.

We present the results from the 1.2 mm continuum image obtained as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin$^2$ within a primary beam response of 50% (10%) and a rms value of $9.3\thinspace{\rm\mu Jy\thinspace beam^{-1}}$. We detect 35 sources at high significance (Fidelity $\geq0.5$), 32 of these have well characterized near-infrared HST counterparts. We estimate the 1.2 mm number counts to flux levels of $<30\thinspace{\rm\mu Jy}$ in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below $S_\nu \sim 0.1$ mJy. In a second approach, we constrain the number counts by using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jack-knifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting $\sim$93% ($\sim$100% if we include the lower fidelity detections) of the total continuum dust emission associated to galaxies in the HUDF. The ancillary data allows us to study the dependence of the 1.2 mm number counts on redshift ($z=0-4$), galaxy dust mass (${\rm M}_{\rm dust}=10^{7}-10^{9}{\rm M}_{\odot}$), stellar mass (${\rm M}_{*}=10^{9}-10^{12}{\rm M}_{\odot}$), and star-formation rate (${\rm SFR}=1-1000\thinspace{\rm M}_{\odot}\thinspace{\rm yr^{-1}}$). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift.

We present a new semi-empirical model for the dust continuum number counts of galaxies at 1.1 millimeter and 850 \micron. Our approach couples an observationally motivated model for the stellar mass and SFR distribution of galaxies with empirical scaling relations to predict the dust continuum flux density of these galaxies. Without a need to tweak the IMF, the model reproduces the currently available observations of the 1.1 millimeter and 850 \micron number counts, including the observed flattening in the 1.1 millimeter number counts below 0.3 mJy \citep{Gonzalez2019numbercounts} and the number counts in discrete bins of different galaxy properties. Predictions of our work include : (1) the galaxies that dominate the number counts at flux densities below 1 mJy (3 mJy) at 1.1 millimeter (850 $\mu$m) have redshifts between $z=1$ and $z=2$, stellar masses of $\sim 5\times10^{10}~\rm{M}_\odot$, and dust masses of $\sim 10^{8}~\rm{M}_\odot$; (2) the flattening in the observed 1.1 millimeter number counts corresponds to the knee of the 1.1 millimeter luminosity function. A similar flattening is predicted for the number counts at 850 $\mu$m; (3) the model reproduces the redshift distribution of current 1.1 millimeter detections; (4) to efficiently detect large numbers of galaxies through their dust continuum, future surveys should scan large areas once reaching a 1.1 millimeter flux density of 0.1 mJy rather than integrating to fainter fluxes. Our modeling framework also suggests that the amount of information on galaxy physics that can be extracted from the 1.1 millimeter and 850 $\mu$m number counts is almost exhausted.

We study the relation between the UV-slope, $\beta$, and the ratio between the IR and UV luminosities of galaxies, IRX, using TNG50, the latest installment of the IllustrisTNG galaxy formation simulations. We select 7280 TNG50 star-forming galaxies with stellar mass above $10^9M_\odot$ at selected redshifts, $0 \leq z \leq 4$, and perform radiative transfer calculations with SKIRT to model the effects of ISM dust on the simulated stellar light. We adopt a Milky Way (MW) type dust and a fixed dust-to-metal ratio of 0.3 throughout and find that TNG50 star-forming main-sequence galaxies agree with the empirically-derived reference IRX-$\beta$ relations at $z \lesssim 1$. There appears to be a redshift-dependent systematic offset with respect to the reference relations, with the TNG50 IRX-$\beta$ median relation steepening towards higher redshifts. This is to first order driven by variations in intrinsic UV-slopes due to different star-formation histories of galaxies selected at different cosmic epochs. Once differences in the intrinsic UV-slope are corrected for, TNG50 galaxies exhibit a scatter of 0.3 dex in IRX at fixed $\beta$ at all studied redshifts. These galaxy-to-galaxy variations are correlated with intrinsic galaxy properties, such that at fixed UV-slope galaxies with higher star--formation rates, star--formation efficiencies, gas metallicities, and stellar masses exhibit larger IRX values. We additionally demonstrate a degeneracy between stellar age and dust type. The combination of young stellar population age with a SMC type of dust for high-redshift ($z=4$) TNG50 galaxies can result in an IRX-$\beta$ relation similar to what is found for low-redshift galaxies with MW dust. This hampers the use of the IRX-$\beta$ relation as a proxy for dust type. We provide a redshift-dependent fitting function for the IRX-$\beta$ relation with MW dust based on our numerical models.

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 present a power spectrum analysis of the ALMA Spectroscopic Survey Large Program (ASPECS LP) data from 84 to 115 GHz. These data predominantly probe small-scale fluctuations ($k=10$-$100$ h Mpc$^{-1}$) in the aggregate CO emission in galaxies at $1 \lesssim z \lesssim 4$. We place an integral constraint on CO luminosity functions (LFs) in this redshift range via a direct measurement of their second moments in the three-dimensional (3D) auto-power spectrum, finding a total CO shot noise power $P_{\textrm{CO,CO}}(k_{\textrm{CO(2-1)}}) \leq 1.9\times10^2$ $\mu$K$^2$ (Mpc h$^{-1}$)$^3$. This upper limit ($3\sigma$) is consistent with the observed ASPECS CO LFs in Decarli et al. 2019, but rules out a large space in the range of $P_{\textrm{CO,CO}}(k_{\textrm{CO(2-1)}})$ inferred from these LFs, which we attribute primarily to large uncertainties in the normalization $\Phi_*$ and knee $L_*$ of the Schechter-form CO LFs at $z > 2$. Also, through power spectrum analyses of ASPECS LP data with 415 positions from galaxies with available optical spectroscopic redshifts, we find that contributions to the observed mean CO intensity and shot noise power of MUSE galaxies are largely accounted for by ASPECS blind detections, though there are $\sim20$% contributions to the CO(2-1) mean intensity due to sources previously undetected in the blind line search. Finally, we sum the fluxes from individual blind CO detections to yield a lower limit on the mean CO surface brightness at 99 GHz of $\langle T_{\textrm{CO}} \rangle = 0.55\pm0.02$ $\mu$K, which we estimate represents $68$-$80$% of the total CO surface brightness at this frequency.

We present new measurements of the cosmic cold molecular gas evolution out to redshift 6 based on systematic mining of the ALMA public archive in the COSMOS deep field (A3COSMOS). Our A3COSMOS dataset contains ~700 galaxies (0.3 < z < 6) with high-confidence ALMA detections in the (sub-)millimeter continuum and multi-wavelength spectral energy distributions (SEDs). Multiple gas mass calibration methods are compared and biases in band conversions (from observed ALMA wavelength to rest-frame Rayleigh-Jeans(RJ)-tail continuum) have been tested. Combining our A3COSMOS sample with ~1,000 CO-observed galaxies at 0 < z < 4 (75% at z < 0.1), we parameterize galaxies' molecular gas depletion time and molecular gas to stellar mass ratio (gas fraction) each as a function of the stellar mass, offset from the star-forming main sequence (Delta MS) and cosmic age (or redshift). Our proposed functional form provides a statistically better fit to current data (than functional forms in the literature), and implies a "downsizing" effect (i.e., more-massive galaxies evolve earlier than less-massive ones) and "mass-quenching" (gas consumption slows down with cosmic time for massive galaxies but speeds up for low-mass ones). Adopting galaxy stellar mass functions and applying our function for gas mass calculation, we for the first time infer the cosmic cold molecular gas density evolution out to redshift 6 and find agreement with CO blind surveys as well as semi-analytic modeling. These together provide a coherent picture of cold molecular gas, SFR and stellar mass evolution in galaxies across cosmic time.

The long-anticipated James Webb Space Telescope (JWST) will be able to directly detect large samples of galaxies at very high redshift. Using the well-established, computationally efficient Santa Cruz semi-analytic model, with recently implemented multiphase gas partitioning and H2-based star formation recipes, we make predictions for a wide variety of galaxy properties for galaxy populations at $z$ = 4-10. In this work, we provide forecasts for the physical properties of high-redshift galaxies and links to their photometric properties. With physical parameters calibrated only to $z\sim0$ observations, our model predictions are in good agreement with current observational constraints on stellar mass and star formation rate distribution functions up to $z \sim 8$. We also provide predictions representing wide, deep, and lensed JWST survey configurations. We study the redshift evolution of key galaxy properties and the scaling relations among them. Taking advantage of our models' high computational efficiency, we study the impact of systematically varying the model parameters. All distribution functions and scaling relations presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

We are just starting to understand the physical processes driving the dramatic change in cosmic star-formation rate between $z\sim 2$ and the present day. A quantity directly linked to star formation is the molecular gas density, which should be measured through independent methods to explore variations due to cosmic variance and systematic uncertainties. We use intervening CO absorption lines in the spectra of mm-bright background sources to provide a census of the molecular gas mass density of the Universe. The data used in this work are taken from ALMACAL, a wide and deep survey utilizing the ALMA calibrator archive. While we report multiple Galactic absorption lines and one intrinsic absorber, no extragalactic intervening molecular absorbers are detected. However, thanks to the large redshift path surveyed ($\Delta z=182$), we provide constraints on the molecular column density distribution function beyond $z\sim 0$. In addition, we probe column densities of N(H$_2$) > 10$^{16}$ atoms~cm$^{-2}$, five orders of magnitude lower than in previous studies. We use the cosmological hydrodynamical simulation IllustrisTNG to show that our upper limits of $\rho ({\rm H}_2)\lesssim 10^{8.3} \text{M}_{\odot} \text{Mpc}^{-3}$ at $0 < z \leq 1.7$ already provide new constraints on current theoretical predictions of the cold molecular phase of the gas. These results are in agreement with recent CO emission-line surveys and are complementary to those studies. The combined constraints indicate that the present decrease of the cosmic star-formation rate history is consistent with an increasing depletion of molecular gas in galaxies compared to $z\sim 2$.

We analyze the interstellar medium properties of a sample of sixteen bright CO line emitting galaxies identified in the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Large Program. This CO$-$selected galaxy sample is complemented by a couple of additional CO line emitters in the UDF that are identified based on their MUSE optical spectroscopic redshifts. The ASPECS CO$-$selected galaxies cover a larger range of star-formation rates and stellar masses compared to literature CO emitting galaxies at $z>1$ for which scaling relations have been established previously. Most of ASPECS CO-selected galaxies follow these established relations in terms of gas depletion timescales and gas fractions as a function of redshift, as well as the star-formation rate-stellar mass relation (`galaxy main sequence'). However, we find that $\sim30\%$ of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with $\sim12\%$ (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low star-formation rates, and yet show substantial molecular gas reservoirs, yielding long gas depletion timescales. Capitalizing on the well-defined cosmic volume probed by our observations, we measure the contribution of galaxies above, below, and on the galaxy main sequence to the total cosmic molecular gas density at different lookback times. We conclude that main sequence galaxies are the largest contributor to the molecular gas density at any redshift probed by our observations (z$\sim$1$-$3). The respective contribution by starburst galaxies above the main sequence decreases from z$\sim$2.5 to z$\sim$1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequence.

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

We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift evolution of $\rho$(H$_2$). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at $z>1$. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at $z$=1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm SFR}$. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history).

The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) provides new constraints for galaxy formation models on the molecular gas properties of galaxies. We compare results from ASPECS to predictions from two cosmological galaxy formation models: the IllustrisTNG hydrodynamical simulations and the Santa Cruz semi-analytic model (SC SAM). We explore several recipes to model the H$_2$ content of galaxies, finding them to be consistent with one another, and take into account the sensitivity limits and survey area of ASPECS. For a canonical CO-to-H$_2$ conversion factor of $\alpha_{\rm CO} = 3.6\,\rm{M}_\odot/(\rm{K}\,\rm{km/s}\,\rm{pc}^{2})$ the results of our work include: (1) the H$_2$ mass of $z>1$ galaxies predicted by the models as a function of their stellar mass is a factor of 2-3 lower than observed; (2) the models do not reproduce the number of H$_2$-rich ($M_{\rm H2} > 3\times 10^{10}\,\rm{M}_\odot$) galaxies observed by ASPECS; (3) the H$_2$ cosmic density evolution predicted by IllustrisTNG (the SC SAM) is in tension (only just agrees) with the observed cosmic density, even after accounting for the ASPECS selection function and field-to-field variance effects. The tension between models and observations at $z>1$ can be alleviated by adopting a CO-to-H$_2$ conversion factor in the range $\alpha_{\rm CO} = 2.0 - 0.8\,\rm{M}_\odot/(\rm{K}\,\rm{km/s}\,\rm{pc}^{2})$. Additional work on constraining the CO-to-H$_2$ conversion factor and CO excitation conditions of galaxies through observations and theory will be necessary to more robustly test the success of galaxy formation models.

We discuss the nature and physical properties of gas-mass selected galaxies in the ALMA spectroscopic survey (ASPECS) of the Hubble Ultra Deep Field (HUDF). We capitalize on the deep optical integral-field spectroscopy from the MUSE HUDF Survey and multi-wavelength data to uniquely associate all 16 line-emitters, detected in the ALMA data without preselection, with rotational transitions of carbon monoxide (CO). We identify ten as CO(2-1) at $1 < z < 2$, five as CO(3-2) at $2 < z < 3$ and one as CO(4-3) at $z = 3.6$. Using the MUSE data as a prior, we identify two additional CO(2-1)-emitters, increasing the total sample size to 18. We infer metallicities consistent with (super-)solar for the CO-detected galaxies at $z \le 1.5$, motivating our choice of a Galactic conversion factor between CO luminosity and molecular gas mass for these galaxies. Using deep Chandra imaging of the HUDF, we determine an X-ray AGN fraction of 20% and 60% among the CO-emitters at $z \sim 1.4$ and $z \sim 2.6$, respectively. Being a CO-flux limited survey, ASPECS-LP detects molecular gas in galaxies on, above and below the main sequence (MS) at $z \sim 1.4$. For stellar masses $\ge 10^{10} (10^{10.5})$ M$_{\odot}$, we detect about 40% (50%) of all galaxies in the HUDF at $1 < z < 2$ ($2 < z < 3$). The combination of ALMA and MUSE integral-field spectroscopy thus enables an unprecedented view on MS galaxies during the peak of galaxy formation.

The ALMA SPECtroscopic Survey in the {\it Hubble} Ultra Deep Field is an ALMA large program that obtained a frequency scan in the 3\,mm band to detect emission lines from the molecular gas in distant galaxies. We here present our search strategy for emission lines and continuum sources in the HUDF. We compare several line search algorithms used in the literature, and critically account for the line-widths of the emission line candidates when assessing significance. We identify sixteen emission lines at high fidelity in our search. Comparing these sources to multi-wavelength data we find that all sources have optical/infrared counterparts. Our search also recovers candidates that have lower significance that can be used statistically to derive, e.g. the CO luminosity function. We apply the same detection algorithm to obtain a sample of six 3 mm continuum sources. All of these are also detected in the 1.2 mm continuum with optical/near-infrared counterparts. We use the continuum sources to compute 3 mm number counts in the sub-mJy regime, and find them to be higher by an order of magnitude than expected for synchrotron-dominated sources. However, the number counts are consistent with those derived at shorter wavelengths (0.85--1.3\,mm) once extrapolating to 3\,mm with a dust emissivity index of $\beta=1.5$, dust temperature of 35\,K and an average redshift of $z=2.5$. These results represent the best constraints to date on the faint end of the 3 mm number counts.

In anticipation of the upcoming deployment of the James Webb Space Telescope (JWST), we present high-redshift predictions by the well-established Santa Cruz semi-analytic model. We update the models by re-calibrating them after adopting cosmological parameters consistent with recent constraints from Planck. We provide predictions for rest-frame UV luminosity functions for galaxy populations over a wide range of $M_\text{UV}$ from $\sim-6$ to $\sim-24$ between $z = 4 - 10$. In addition, we present the corresponding predictions for observed-frame galaxy number counts in different redshift bins in the full set of NIRCam filters. We provide predictions of the quantitative effect on these observables of varying the physical recipes implemented in the models, such as the molecular gas depletion time (star formation efficiency) scalings or the scalings of outflow rates driven by stars and supernovae with galaxy circular velocity. Based on these results, we discuss what may be learned about the physical processes that shape galaxy formation from JWST observations of galaxy number densities at different intrinsic luminosities. All data tables for the results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

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.

Large-scale cosmological simulations of galaxy formation currently do not resolve the densities at which molecular hydrogen forms, implying that the atomic-to-molecular transition must be modeled either on the fly or in postprocessing. We present an improved postprocessing framework to estimate the abundance of atomic and molecular hydrogen and apply it to the IllustrisTNG simulations. We compare five different models for the atomic-to-molecular transition, including empirical, simulation-based, and theoretical prescriptions. Most of these models rely on the surface density of neutral hydrogen and the ultraviolet (UV) flux in the Lyman-Werner band as input parameters. Computing these quantities on the kiloparsec scales resolved by the simulations emerges as the main challenge. We show that the commonly used Jeans length approximation to the column density of a system can be biased and exhibits large cell-to-cell scatter. Instead, we propose to compute all surface quantities in face-on projections and perform the modeling in two dimensions. In general, the two methods agree on average, but their predictions diverge for individual galaxies and for models based on the observed midplane pressure of galaxies. We model the UV radiation from young stars by assuming a constant escape fraction and optically thin propagation throughout the galaxy. With these improvements, we find that the five models for the atomic-to-molecular transition roughly agree on average but that the details of the modeling matter for individual galaxies and the spatial distribution of molecular hydrogen. We emphasize that the estimated molecular fractions are approximate due to the significant systematic uncertainties.

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

Gas-rich galaxies are selected efficiently via quasar absorption lines. Recently, a new perspective on such absorption-selected systems has opened up by studying the molecular gas content of absorber host galaxies using ALMA CO emission line observations. Here, we present an analysis of multiple CO transitions ($L'_{\rm CO} \sim 10^9$ K km s$^{-1}$) in two $z \sim 0.5$ galaxies associated with one Ly$\alpha$ absorber towards J0238+1636. The CO spectral line energy distribution (CO SLED) of these galaxies appear distinct from that of typical star-forming galaxies at similar redshifts and is comparable with that of luminous infrared galaxies or AGN. Indeed, these galaxies are associated with optically identified AGN activity. We infer that the CO line ratios and the $\alpha_{\rm CO}$ conversion factor differ from the Galactic values. Our findings suggest that at least a fraction of absorption selected systems shows ISM conditions deviating from those of normal star-forming galaxies. For a robust molecular gas mass calculation, it is therefore important to construct the CO SLED. Absorption-line-selection identifies systems with widely distributed gas, which may preferentially select interacting galaxies, which in turn will have more excited CO SLEDs than isolated galaxies. Furthermore, we raise the question whether quasar absorbers preferentially trace galaxy overdensities.

A tightly correlated star formation rate-stellar mass relation of star forming galaxies, or star-forming sequence (SFS), is a key feature in galaxy property-space that is predicted by modern galaxy formation models. We present a flexible data-driven approach for identifying this SFS over a wide range of star formation rates and stellar masses using Gaussian mixture modeling (GMM). Using this method, we present a consistent comparison of the $z=0$ SFSs of central galaxies in the Illustris, EAGLE, and Mufasa hydrodynamic simulations and the Santa Cruz semi-analytic model (SC-SAM), alongside data from the Sloan Digital Sky Survey. We find, surprisingly, that the amplitude of the SFS varies by up to ${\sim} 0.7\,\mathrm{dex}$ (factor of ${\sim} 5$) among the simulations with power-law slopes range from $0.7$ to $1.2$. In addition to the SFS, our GMM method also identifies sub-components in the star formation rate-stellar mass relation corresponding to star-burst, transitioning, and quiescent sub-populations. The hydrodynamic simulations are similarly dominated by SFS and quiescent sub-populations unlike the SC-SAM, which predicts substantial fractions of transitioning and star-burst galaxies at stellar masses above and below $10^{10} M_\odot$, respectively. All of the simulations also produce an abundance of low-mass quiescent central galaxies in apparent tension with observations. These results illustrate that, even among models that well reproduce many observables of the galaxy population, the $z=0$ SFS and other sub-populations still show marked differences that can provide strong constraints on galaxy formation models.

We present ALMA Band 9 observations of the [C II]158um emission for a sample of 10 main-sequence galaxies at redshift z ~ 2, with typical stellar masses (log M*/Msun ~ 10.0 - 10.9) and star formation rates (~ 35 - 115 Msun/yr). Given the strong and well understood evolution of the interstellar medium from the present to z = 2, we investigate the behaviour of the [C II] emission and empirically identify its primary driver. We detect [C II] from six galaxies (four secure, two tentative) and estimate ensemble averages including non detections. The [C II]-to-infrared luminosity ratio (L[C II]/LIR) of our sample is similar to that of local main-sequence galaxies (~ 2 x 10^-3), and ~ 10 times higher than that of starbursts. The [C II] emission has an average spatial extent of 4 - 7 kpc, consistent with the optical size. Complementing our sample with literature data, we find that the [C II] luminosity correlates with galaxies' molecular gas mass, with a mean absolute deviation of 0.2 dex and without evident systematics: the [C II]-to-H2 conversion factor (alpha_[C II] ~ 30 Msun/Lsun) is largely independent of galaxies' depletion time, metallicity, and redshift. [C II] seems therefore a convenient tracer to estimate galaxies' molecular gas content regardless of their starburst or main-sequence nature, and extending to metal-poor galaxies at low and high redshifts. The dearth of [C II] emission reported for z > 6 - 7 galaxies might suggest either a high star formation efficiency or a small fraction of UV light from star formation reprocessed by dust.

Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of "Walking the Line", a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodology. Emission lines across the spectrum from the millimeter to the UV were discussed, with most of the focus on the interstellar medium, but also some topics on the circumgalactic medium. The most important quality of a useful model is a good synergy with observations and experiments. Challenges in simulating line emission are identified, some of which are already being worked upon, and others that must be addressed in the future for models to agree with observations. Recent advances in several areas aiming at achieving that synergy are summarized here, from micro-physical to galactic and circum-galactic scale.

We present the dust mass function (DMF) of 15,750 galaxies with redshift $z< 0.1$, drawn from the overlapping area of the GAMA and {\it H-}ATLAS surveys. The DMF is derived using the density corrected $V_{\rm max}$ method, where we estimate $V_{\rm max}$ using: (i) the normal photometric selection limit ($pV_{\rm max}$) and (ii) a bivariate brightness distribution (BBD) technique, which accounts for two selection effects. We fit the data with a Schechter function, and find $M^{*}=(4.65\pm0.18)\times 10^{7}\,h^2_{70}\, M_{\odot}$, $\alpha=(1.22\pm 0.01)$, $\phi^{*}=(6.26\pm 0.28)\times 10^{-3}\,h^3_{70}\,\rm Mpc^{-3}\,dex^{-1}$. The resulting dust mass density parameter integrated down to $10^4\,M_{\odot}$ is $\Omega_{\rm d}=(1.11 \pm0.02)\times 10^{-6}$ which implies the mass fraction of baryons in dust is $f_{m_b}=(2.40\pm0.04)\times 10^{-5}$; cosmic variance adds an extra 7-17\,per\,cent uncertainty to the quoted statistical errors. Our measurements have fewer galaxies with high dust mass than predicted by semi-analytic models. This is because the models include too much dust in high stellar mass galaxies. Conversely, our measurements find more galaxies with high dust mass than predicted by hydrodynamical cosmological simulations. This is likely to be from the long timescales for grain growth assumed in the models. We calculate DMFs split by galaxy type and find dust mass densities of $\Omega_{\rm d}=(0.88\pm0.03)\times 10^{-6}$ and $\Omega_{\rm d}=(0.060\pm0.005)\times 10^{-6}$ for late-types and early-types respectively. Comparing to the equivalent galaxy stellar mass functions (GSMF) we find that the DMF for late-types is well matched by the GMSF scaled by $(8.07\pm0.35) \times 10^{-4}$.

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.

Studying the flow of baryons into and out of galaxies is an important part of understanding the evolution of galaxies over time. We present a detailed case study of the environment around an intervening Ly $\alpha$ absorption line system at $z_{\rm abs} = 0.633$, seen towards the quasar J0423$-$0130 ($z_{\rm QSO} = 0.915$). We detect with ALMA the $^{12}$CO(2--1), $^{12}$CO(3--2) and $1.2$~mm continuum emission from a galaxy at the redshift of the Ly $\alpha$ absorber at a projected distance of $135$ kpc. From the ALMA detections, we infer ISM conditions similar to those in low redshift Luminous Infrared Galaxies. DDT MUSE integral field unit observations reveal the optical counterpart of the $^{12}$CO emission line source and three additional emission line galaxies at the absorber redshift, which together form a galaxy group. The $^{12}$CO emission line detections originate from the most massive galaxy in this group. While we cannot exclude that we miss a fainter host, we reach a dust-uncorrected star-formation rate (SFR) limit of > $0.3 \text{M}_{\odot} \text{ yr}^{-1}$ within $100$ kpc from the sightline to the background quasar. We measure the dust-corrected SFR (ranging from $3$ to $50$ M$_{\odot}$ yr$^{-1}$), the morpho-kinematics and the metallicities of the four group galaxies to understand the relation between the group and the neutral gas probed in absorption. We find that the Ly $\alpha$ absorber traces either an outflow from the most massive galaxy or intra-group gas. This case study illustrates the power of combining ALMA and MUSE to obtain a census of the cool baryons in a bounded structure at intermediate redshift.

The initial mass function (IMF) is a key ingredient in many studies of galaxy formation and evolution. Although the IMF is often assumed to be universal, there is continuing evidence that it is not universal. Spectroscopic studies that derive the IMF of the unresolved stellar populations of a galaxy often assume that this spectrum can be described by a single stellar population (SSP). To alleviate these limitations, in this paper we have developed a unique hierarchical Bayesian framework for modelling composite stellar populations (CSPs). Within this framework we use a parameterized IMF prior to regulate a direct inference of the IMF. We use this new framework to determine the number of SSPs that is required to fit a set of realistic CSP mock spectra. The CSP mock spectra that we use are based on semi-analytic models and have an IMF that varies as a function of stellar velocity dispersion of the galaxy. Our results suggest that using a single SSP biases the determination of the IMF slope to a higher value than the true slope, although the trend with stellar velocity dispersion is overall recovered. If we include more SSPs in the fit, the Bayesian evidence increases significantly and the inferred IMF slopes of our mock spectra converge, within the errors, to their true values. Most of the bias is already removed by using two SSPs instead of one. We show that we can reconstruct the variable IMF of our mock spectra for signal-to-noise ratios exceeding $\sim$75.

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

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

The use of ultraviolet (UV) emission as a tracer of galaxy star-formation rate (SFR) is hampered by dust obscuration. The empirical relationship between UV slope, $\beta$, and the ratio between far-infrared and UV luminosity, IRX, is commonly employed to account for obscured UV emission. We present a simple model that explores the physical origin of variations in the IRX - $\beta$ dust attenuation relation. A relative increase in FUV attenuation compared to NUV attenuation and an increasing stellar population age cause variations towards red UV slopes for a fixed IRX. Dust geometry effects (turbulence, dust screen with holes, mixing of stars within the dust screen, two-component dust model) cause variations towards blue UV slopes. Poor photometric sampling of the UV spectrum causes additional observational variations. We provide an analytic approximation for the IRX - $\beta$ relation invoking a subset of the explored physical processes (dust type, stellar population age, turbulence). We discuss observed variations in the IRX - $\beta$ relation for local (sub-galactic scales) and high-redshift (normal and dusty star-forming galaxies, galaxies during the epoch of reionization) galaxies in the context of the physical processes explored in our model. High spatial resolution imaging of the UV and sub-mm emission of galaxies can constrain the IRX - $\beta$ dust attenuation relation for different galaxy types at different epochs, where different processes causing variations may dominate. These constraints will allow the use of the IRX - $\beta$ relation to estimate intrinsic SFRs of galaxies, despite the lack of a universal relation.

We study the dust content of galaxies from $z=0$ to $z=9$ in semi-analytic models of galaxy formation that include new recipes to track the production and destruction of dust. We include condensation of dust in stellar ejecta, the growth of dust in the interstellar medium (ISM), the destruction of dust by supernovae and in the hot halo, and dusty winds and inflows. The rate of dust growth in the ISM depends on the metallicity and density of molecular clouds. Our fiducial model reproduces the relation between dust mass and stellar mass from $z=0$ to $z=7$, the number density of galaxies with dust masses less than $10^{8.3}\,\rm{M}_\odot$, and the cosmic density of dust at $z=0$. The model accounts for the double power-law trend between dust-to-gas (DTG) ratio and gas-phase metallicity of local galaxies and the relation between DTG ratio and stellar mass. The dominant mode of dust formation is dust growth in the ISM, except for galaxies with $M_*<10^7\,\rm{M}_\odot$, where condensation of dust in supernova ejecta dominates. The dust-to-metal ratio of galaxies depends on the gas-phase metallicity, unlike what is typically assumed in cosmological simulations. Model variants including higher condensation efficiencies, a fixed timescale for dust growth in the ISM, or no growth at all reproduce some of the observed constraints, but fail to simultaneously reproduce the shape of dust scaling relations and the dust mass of high-redshift galaxies.

We present ALMA detections of the [CI] 1-0, CO J=3-2, and CO J=4-3 emission lines, as well as the ALMA band 4 continuum for a compact star-forming galaxy (cSFG) at z=2.225, 3D-HST GS30274. As is typical for cSFGs, this galaxy has a stellar mass of $1.89 \pm 0.47\,\times 10^{11}\,\rm{M}_\odot$, with a star formation rate of $214\pm44\,\rm{M}_\odot\,\rm{yr}^{-1}$ putting it on the star-forming `main-sequence', but with an H-band effective radius of 2.5 kpc, making it much smaller than the bulk of `main-sequence' star-forming galaxies. The intensity ratio of the line detections yield an ISM density (~ 6 $\times 10^{4}\,\rm{cm}^{-3}$) and a UV-radiation field ( ~2 $\times 10^4\,\rm{G}_0$), similar to the values in local starburst and ultra-luminous infrared galaxy environments. A starburst phase is consistent with the short depletion times ($t_{\rm H2, dep} \leq 140$ Myr) we find using three different proxies for the H2 mass ([CI], CO, dust mass). This depletion time is significantly shorter than in more extended SFGs with similar stellar masses and SFRs. Moreover, the gas fraction of 3D-HST GS30274 is smaller than typically found in extended galaxies. We measure the CO and [CI] kinematics and find a FWHM line width of ~$750 \pm 41 $ km s$^{-1}$. The CO and [CI] FWHM are consistent with a previously measured H$\alpha$ FWHM for this source. The line widths are consistent with gravitational motions, suggesting we are seeing a compact molecular gas reservoir. A previous merger event, as suggested by the asymmetric light profile, may be responsible for the compact distribution of gas and has triggered a central starburst event. This event gives rise to the starburst-like ISM properties and short depletion times. The centrally located and efficient star formation is quickly building up a dense core of stars, responsible for the compact distribution of stellar light in 3D-HST GS30274.

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

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

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

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

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

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

Now that ALMA is reaching its full capabilities, observations of sub-mm emission line deep fields become feasible. We couple a semi-analytic model of galaxy formation with a radiative transfer code to make predictions for the luminosity function of CO J=1-0 up to CO J=6-5 and [CII] at redshifts z=0-6. We find that: 1) our model correctly reproduces the CO and [CII] emission of low- and high-redshift galaxies and reproduces the available constraints on the CO luminosity function at z<2.75; 2) we find that the CO and [CII] luminosity functions of galaxies increase from z = 6 to z = 4, remain relatively constant till z = 1 and rapidly decrease towards z = 0. The galaxies that are brightest in CO and [CII] are found at z~2; 3) the CO J=3-2 emission line is most favourable to study the CO luminosity and global H2 mass content of galaxies, because of its brightness and observability with currently available sub-mm and radio instruments; 4) the luminosity functions of high-J CO lines show stronger evolution than the luminosity functions of low-J CO lines; 5) our model barely reproduces the available constraints on the CO and [CII] luminosity function of galaxies at z>1.5 and the CO luminosity of individual galaxies at intermediate redshifts. We argue that this is driven by a lack of cold gas in galaxies at intermediate redshifts as predicted by cosmological simulations of galaxy formation.

We derive the total cold gas, atomic hydrogen, and molecular gas masses of approximately 24 000 galaxies covering four decades in stellar mass at redshifts 0.5 < z < 3.0, taken from the CANDELS survey. Our inferences are based on the inversion of a molecular hydrogen based star formation law, coupled with a prescription to separate atomic and molecular gas. We find that: 1) there is an increasing trend between the inferred cold gas (HI and H2), HI, and H2 mass and the stellar mass of galaxies down to stellar masses of 10^8 Msun already in place at z = 3; 2) the molecular fractions of cold gas increase with increasing stellar mass and look-back time; 3) there is hardly any evolution in the mean HI content of galaxies at fixed stellar mass; 4) the cold gas fraction and relative amount of molecular hydrogen in galaxies decrease at a relatively constant rate with time, independent of stellar mass; 5) there is a large population of low-stellar mass galaxies dominated by atomic gas. These galaxies are very gas rich, but only a minor fraction of their gas is molecular; 6) the ratio between star-formation rate (SFR) and inferred total cold gas mass (HI + H2) of galaxies (i.e., star-formation efficiency; SFE) increases with star-formation at fixed stellar masses. Due to its simplicity, the presented approach is valuable to assess the impact of selection biases on small samples of directly-observed gas masses and to extend scaling relations down to stellar mass ranges and redshifts that are currently difficult to probe with direct measurements of gas content.

We implement physically motivated recipes for partitioning cold gas into different phases (atomic, molecular, and ionized) in galaxies within semi-analytic models of galaxy formation based on cosmological merger trees. We then model the conversion of molecular gas into stars using empirical recipes motivated by recent observations. We explore the impact of these new recipes on the evolution of fundamental galaxy properties such as stellar mass, star formation rate (SFR), and gas and stellar phase metallicity. We present predictions for stellar mass functions, stellar mass vs. SFR relations, and cold gas phase and stellar mass-metallicity relations for our fiducial models, from redshift $z\sim 6$ to the present day. In addition we present predictions for the global SFR, mass assembly history, and cosmic enrichment history. We find that the predicted stellar properties of galaxies (stellar mass, SFR, metallicity) are remarkably insensitive to the details of the recipes used for partitioning gas into HI and H$_2$. We see significant sensitivity to the recipes for H$_2$ formation only in very low mass halos, which host galaxies that are not detectable with current observational facilities except very nearby. The properties of low-mass galaxies are also quite insensitive to the details of the recipe used for converting H$_2$ into stars, while the formation epoch of massive galaxies does depend on this significantly. (Abridged)

We present a semi-empirical model to infer the atomic and molecular hydrogen content of galaxies as a function of halo mass and time. Our model combines the SFR-halo mass-redshift relation (constrained by galaxy abundances) with inverted SFR-surface density relations to infer galaxy H I and H2 masses. We present gas scaling relations, gas fractions, and mass functions from z = 0 to z = 3 and the gas properties of galaxies as a function of their host halo masses. Predictions of our work include: 1) there is a ~ 0.2 dex decrease in the H I mass of galaxies as a function of their stellar mass since z = 1.5, whereas the H2 mass of galaxies decreases by > 1 dex over the same period. 2) galaxy cold gas fractions and H2 fractions decrease with increasing stellar mass and time. Galaxies with M* > 10^10 Msun are dominated by their stellar content at z < 1, whereas less-massive galaxies only reach these gas fractions at z = 0. We find the strongest evolution in relative gas content at z < 1.5. 3) the SFR to gas mass ratio decreases by an order of magnitude from z = 3 to z = 0. This is consistent with lower H2 fractions; these lower fractions in combination with smaller gas reservoirs correspond to decreased present-day galaxy SFRs. 4) an H2-based star- formation relation can simultaneously fuel the evolution of the cosmic star-formation and reproduce the observed weak evolution in the cosmic HI density. 5) galaxies residing in haloes with masses near 10^12 Msun are most efficient at obtaining large gas reservoirs and forming H2 at all redshifts. These two effects lie at the origin of the high star-formation efficiencies in haloes with the same mass.

We combine a semi-analytic model of galaxy formation, tracking atomic and molecular phases of cold gas, with a three-dimensional radiative-transfer and line tracing code to study the sub-mm emission from atomic and molecular species (CO, HCN, [CI], [CII], [OI]) in galaxies. We compare the physics that drives the formation of stars at the epoch of peak star formation (SF) in the Universe (z = 2.0) with that in local galaxies. We find that normal star-forming galaxies at high redshift have much higher CO-excitation peaks than their local counterparts and that CO cooling takes place at higher excitation levels. CO line ratios increase with redshift as a function of galaxy star-formation rate, but are well correlated with H2 surface density independent of redshift. We find an increase in the [OI]/[CII] line ratio in typical star-forming galaxies at z = 1.2 and z = 2.0 with respect to counterparts at z = 0. Our model results suggest that typical star-forming galaxies at high redshift consist of much denser and warmer star-forming clouds than their local counterparts. Galaxies belonging to the tail of the SF activity peak at z = 1.2 are already less dense and cooler than counterparts during the actual peak of SF activity (z = 2.0). We use our results to discuss how future ALMA surveys can best confront our predictions and constrain models of galaxy formation.

We report Plateau de Bure Interferometer (PdBI) 1.1 mm continuum imaging towards two extremely red H-[4.5]>4 (AB) galaxies at z>3, which we have previously discovered making use of Spitzer SEDS and Hubble Space Telescope CANDELS ultra-deep images of the Ultra Deep Survey field. One of our objects is detected on the PdBI map with a 4.3 sigma significance, corresponding to Snu(1.1mm)=(0.78 +/- 0.18) mJy. By combining this detection with the Spitzer 8 and 24 micron photometry for this source, and SCUBA2 flux density upper limits, we infer that this galaxy is a composite active galactic nucleus/star-forming system. The infrared (IR)-derived star formation rate is SFR~(200 +/- 100) Msun/yr, which implies that this galaxy is a higher-redshift analogue of the ordinary ultra-luminous infrared galaxies (ULIRGs) more commonly found at z~2-3. In the field of the other target, we find a tentative 3.1 sigma detection on the PdBI 1.1 mm map, but 3.7 arcsec away of our target position, so it likely corresponds to a different object. In spite of the lower significance, the PdBI detection is supported by a close SCUBA2 3.3 sigma detection. No counterpart is found on either the deep SEDS or CANDELS maps, so, if real, the PdBI source could be similar in nature to the sub-millimetre source GN10. We conclude that the analysis of ultra-deep near- and mid-IR images offers an efficient, alternative route to discover new sites of powerful star formation activity at high redshifts.

We study the evolution of atomic and molecular gas in galaxies in semi-analytic models of galaxy formation that include new modeling of the partitioning of cold gas in galactic discs into atomic, molecular, and ionised phases. We adopt two scenarios for the formation of molecules: one pressure-based and one metallicity-based. We find that both recipes successfully reproduce the gas fractions and gas-to-stellar mass ratios of HI and H2 in local galaxies,as well as the HI and H2 disc sizes up to z<2. We reach good agreement with the locally observed HI and H2 mass function, although both recipes slightly overpredict the low-mass end of the HI mass function. Both of our models predict that the high-mass end of the HI mass function remains nearly constant at redshifts z < 2.0. The metallicity-based recipe yields a higher cosmic density of cold gas and much lower cosmic H2 fraction over the entire redshift range probed than the pressure based recipe. These strong differences in HI mass function and cosmic density between the two recipes are driven by low mass galaxies (log (M*/Msun) < 7) residing in low mass halos (log (Mvir/Msun) < 10). Both recipes predict that galaxy gas fractions remain high from z ~ 6-3 and drop rapidly at lower redshift. The galaxy H2 fractions show a similar trend, but drop even more rapidly. We provide predictions for the CO J = 1-0 luminosity of galaxies, which will be directly comparable with observations with sub-mm and radio instruments.

We investigate the properties of damped Ly{\alpha} absorption systems (DLAs) in semi-analytic models of galaxy formation, including partitioning of cold gas in galactic discs into atomic, molecular, and ionized phases with a molecular gas-based star formation recipe. We investigate two approaches for partitioning gas into these constituents: a pressure-based and a metallicity-based recipe. We identify DLAs by passing lines of sight through our simulations to compute HI column densities. We find that models with "standard" gas radial profiles - where the average specific angular momentum of the gas disc is equal to that of the host dark matter halo - fail to reproduce the observed column density distribution of DLAs. These models also fail to reproduce the distribution of velocity widths {\Delta}v, overproducing low {\Delta}v relative to high {\Delta}v systems. Models with "extended" radial gas profiles - corresponding to gas discs with higher specific angular momentum - are able to reproduce quite well the column density distribution of absorbers over the column density range 19 < log NHI < 22.5 in the redshift range 2 < z < 3.5. The model with pressure-based gas partitioning also reproduces the observed line density of DLAs, HI gas density, and {\Delta}v distribution at z < 3 remarkably well. However all of the models investigated here underproduce DLAs and the HI gas density at z > 3. If this is the case, the flatness in the number of DLAs and HI gas density over the redshift interval 0 < z < 5 may be due to a cosmic coincidence where the majority of DLAs at z > 3 arise from intergalactic gas in filaments while those at z < 3 arise predominantly in galactic discs. We further investigate the dependence of DLA metallicity on redshift and {\Delta}v, and find reasonably good agreement with the observations, particularly when including the effects of metallicity gradients (abbrv.).

One key piece of information missing from high redshift galaxy surveys is the galaxies' cold gas contents. We present a new method to indirectly determine cold gas surface densities and integrated gas masses from galaxy star formation rates and to separate the atomic and molecular gas components. Our predicted molecular and total gas surface densities and integrated masses are in very good agreement with direct measurements quoted in the literature for low and high-z galaxies. We apply this method to predict the gas content for a sample of $\sim 57000$ galaxies in the COSMOS field at $0.5 \leq z \leq 2.0$, selected to have $I_{AB} < 24$ mag. This approach allows us to investigate in detail the redshift evolution of galaxy cold and molecular gas content versus stellar mass and to provide fitting formulae for galaxy gas fractions. We find a clear trend between galaxy gas fraction, molecular gas fraction and stellar mass with redshift, suggesting that massive galaxies consume and/or expel their gas at higher redshift than less massive objects and have lower fractions of their gas in molecular form. The characteristic stellar mass separating gas- from stellar-dominated galaxies decreases with time. This indicates that massive galaxies reach a gas-poor state earlier than less massive objects. These trends can be considered to be another manifestation of downsizing in star formation activity.

NGC 2903 is a nearby barred spiral with an active starburst in the center and Hii regions distributed along its bar. We aim to analyse the star formation properties in the bar region of NGC 2903 and study the links with the typical bar morphological features. A combination of space and ground-based data from the far-ultraviolet to the sub-millimeter spectral ranges is used to create a panchromatic view of the NGC 2903 bar. We produce two catalogues: one for the current star formation regions, as traced by the halpha compact emission, and a second one for the ultraviolet (UV) emitting knots, containing positions and luminosities. From them we have obtained ultraviolet colours, star formation rates, dust attenuation and halpha EWs, and their spatial distribution have been analysed. Stellar cluster ages have been estimated using stellar population synthesis models (Starburst99). NGC 2903 is a complex galaxy, with a very different morphology on each spectral band. The CO(J=1-0) and the 3.6 micron emission trace each other in a clear barred structure, while the halpha leads both components and it has an s-shape distribution. The UV emission is patchy and does not resemble a bar. The UV emission is also characterised by a number of regions located symmetrically with respect to the galaxy center, almost perpendicular to the bar, in a spiral shape covering the inner ~2.5 kpc. These regions do not show a significant halpha nor 24 micron emission. We have estimated ages for these regions ranging from 150 to 320 Myr, being older than the rest of the UV knots, which have ages lower than 10 Myr. The SFR calculated from the UV emission is ~0.4 M$_{\odot}$/yr, compatible with the SFR as derived from halpha calibrations (M$_{\odot}$/yr).