Search

We present the results of a James Webb Space Telescope (JWST) NIRCam and NIRSpec investigation into the young massive star cluster (YMC) population of NGC 3256, the most cluster-rich luminous infrared galaxy (LIRG) in the Great Observatories All Sky LIRG Survey. We detect 3061 compact YMC candidates with a $S/N \geq 3$ at F150W, F200W, and F335M. Based on yggdrasil stellar population models, we identify 116/3061 sources with F150W - F200W $> 0.47$ and F200W - F355M $> -1.37$ colors suggesting they are young (t $\leq 5$ Myr), dusty ($A_{V} = 5 - 15$), and massive ($M_{\odot} > 10^{5}$). This increases the sample of dust-enshrouded YMCs detected in this system by an order of magnitude relative to previous HST studies. With NIRSpec IFU pointings centered on the northern and southern nucleus, we extract the Pa$\alpha$ and 3.3$\mu$m PAH equivalent widths for 8 bright and isolated YMCs. Variations in both the F200W - F335M color and 3.3$\mu$m PAH emission with the Pa$\alpha$ line strength suggest a rapid dust clearing ($< 3 - 4$ Myr) for the emerging YMCs in the nuclei of NGC 3256. Finally, with both the age and dust emission accurately measured we use yggdrasil to derive the color excess (E(B - V)) for all 8 YMCs. We demonstrate that YMCs with strong 3.3$\mu$m PAH emission (F200W - F335M $> 0$) correspond to sources with E(B - V) $> 3$, which are typically missed in UV-optical studies. This underscores the importance of deep near-infrared imaging for finding and characterizing these very young and dust-embedded sources., Comment: 18 pages, 7 figures, 1 table; Accepted for publication in ApJL

We present new multi-frequency ALMA continuum observations of the massive [$\log_{10}(M_\star/M_\odot) = 10.3_{-0.2}^{+0.1}$], UV-luminous [$M_\mathrm{UV} = -21.7 \pm 0.2$] $z=7.31$ galaxy REBELS-25 in Bands 3, 4, 5, and 9. Combining the new observations with previously-taken data in Bands 6 and 8, we cover the dust continuum emission of the galaxy in six distinct bands -- spanning rest-frame $50-350\,\mu$m -- enabling simultaneous constraints on its dust mass ($M_\mathrm{dust}$), temperature ($T_\mathrm{dust}$) and emissivity index ($\beta_\mathrm{IR}$) via modified blackbody fitting. Given a fiducial model of optically thin emission, we infer a cold dust temperature of $T_\mathrm{dust} = 32_{-6}^{+9}\,$K and a high dust mass of $\log_{10}(M_\mathrm{dust}/M_\odot) = 8.2_{-0.4}^{+0.6}$, and moderately optically thick dust does not significantly alter these estimates. If we assume dust production is solely through supernovae (SNe), the inferred dust yield would be high, $y = 0.7_{-0.4}^{+2.3}\,M_\odot$ per SN. Consequently, we argue grain growth in the interstellar medium of REBELS-25 also contributes to its dust build-up. This is supported by the steep dust emissivity index $\beta_\mathrm{IR} = 2.5 \pm 0.4$ we measure for REBELS-25, as well as by its high stellar mass, dense interstellar medium, and metal-rich nature. Our results suggest that constraining the dust emissivity indices of high-redshift galaxies is important not only to mitigate systematic uncertainties in their dust masses and obscured star formation rates, but also to assess if dust properties evolve across cosmic time. We present an efficient observing setup to do so with ALMA, combining observations of the peak and Rayleigh-Jeans tail of the dust emission., Comment: 14 pages main text, 6 figures + appendices; re-submitted to MNRAS after a positive referee report

We present high resolution ($\sim0.14$" = 710 pc) ALMA [CII] 158$\mu$m and dust continuum follow-up observations of REBELS-25, a [CII]-luminous ($L_{\mathrm{[CII]}}=(1.7\pm0.2)\times 10^9 \mathrm{L_{\odot}}$) galaxy at redshift $z=7.3065\pm0.0001$. These high resolution, high signal-to-noise observations allow us to study the sub-kpc morphology and kinematics of this massive ($M_* = 8^{+4}_{-2} \times 10^9 \mathrm{M_{\odot}}$) star-forming (SFR$_{\mathrm{UV+IR}} = 199^{+101}_{-63} \mathrm{M_{\odot}} \mathrm{yr}^{-1}$) galaxy in the Epoch of Reionisation. By modelling the kinematics with $^{\mathrm{3D}}$BAROLO, we find it has a low velocity dispersion ($\bar{\sigma} = 33 \pm 9$ km s$^{-1}$) and a high ratio of ordered-to-random motion ($V_{\mathrm{rot, ~max}}/\bar{\sigma} = 11 ^{+8}_{-4}$), indicating that REBELS-25 is a dynamically cold disc. Additionally, we find that the [CII] distribution is well fit by a near-exponential disc model, with a S\'ersic index, $n$, of $1.3 \pm 0.2$, and we see tentative evidence of more complex non-axisymmetric structures suggestive of a bar in the [CII] and dust continuum emission. By comparing to other high spatial resolution cold gas kinematic studies, we find that dynamically cold discs seem to be more common in the high redshift Universe than expected based on prevailing galaxy formation theories, which typically predict more turbulent and dispersion-dominated galaxies in the early Universe as an outcome of merger activity, gas accretion and more intense feedback. This higher degree of rotational support seems instead to be consistent with recent cosmological simulations that have highlighted the contrast between cold and warm ionised gas tracers, particularly for massive galaxies. We therefore show that dynamically settled disc galaxies can form as early as 700 Myr after the Big Bang., Comment: Submitted to MNRAS

We present James Webb Space Telescope (JWST) Integral Field Spectrograph observations of NGC 3256, a local infrared-luminous late-stage merging system with two nuclei roughly 1$\;\rm{kpc}$ apart, both of which have evidence of cold molecular outflows. Using JWST/NIRSpec and MIRI datasets, we investigate this morphologically complex system on spatial scales of $<$100$\;\rm{pc}$, where we focus on the warm molecular H$_2$ gas surrounding the nuclei. We detect collimated outflowing warm H$_2$ gas originating from the southern nucleus, though we do not find significant outflowing H$_2$ gas surrounding the northern nucleus. We measure maximum intrinsic outflow velocities of $\sim$1,000$\;\rm{km}\;\rm{s}^{-1}$, which extend out to a distance of 0.7$\;\rm{kpc}$. Based on H$_2$ S(7)/S(1) ratios, we find a larger fraction of warmer gas near the S nucleus, which decreases with increasing distance from the nucleus, signifying the southern nucleus as a primary source of H$_2$ heating. The gas mass of the warm H$_2$ outflow component is estimated to be $M\rm{_{warm,out}}=(1.4\pm0.2)\times10^6\;\rm{M}_{\odot}$, as much as 6$\%$ of the cold H$_2$ mass estimated using ALMA CO data. The outflow time scale is about $7\times10^5\;\rm{yr}$, resulting in a mass outflow rate $\dot{M}\rm{_{warm,out}}=2.0\pm0.8\;\rm{M}_{\odot}\;\rm{yr}^{-1}$ and kinetic power $P\rm{_{warm,out}}\;\sim\;4\times10^{41}\;\rm{erg}\;\rm{s}^{-1}$. Lastly, regions within our 3.0"x3.0" NIRSpec data where the outflowing gas reside show high [Fe II]/Pa$\beta$ and H$_2$/Br$\gamma$ line ratios, indicate enhanced mechanical heating caused by the outflows. The fluxes and ratios of Polycyclic Aromatic Hydrocarbons (PAH) in these regions are not significantly different compared to those elsewhere in the disk, suggesting the outflows may not significantly alter the PAH ionization state or grain size., Comment: 19 pages, 9 figures. Accepted by ApJ

The enormous increase in mid-IR sensitivity and spatial and spectral resolution provided by the JWST spectrographs enables, for the first time, detailed extragalactic studies of molecular vibrational bands. This opens an entirely new window for the study of the molecular interstellar medium in luminous infrared galaxies (LIRGs). We present a detailed analysis of rovibrational bands of gas-phase CO, H$_2$O, C$_2$H$_2$ and HCN towards the heavily-obscured eastern nucleus of the LIRG VV 114, as observed by NIRSpec and MIRI MRS. Spectra extracted from apertures of 130 pc in radius show a clear dichotomy between the obscured AGN and two intense starburst regions. We detect the 2.3 $\mu$m CO bandheads, characteristic of cool stellar atmospheres, in the star-forming regions, but not towards the AGN. Surprisingly, at 4.7 $\mathrm{\mu}$m we find highly-excited CO ($T_\mathrm{ex} \approx 700-800$ K out to at least rotational level $J = 27$) towards the star-forming regions, but only cooler gas ($T_\mathrm{ex} \approx 200$ K) towards the AGN. We conclude that only mid-infrared pumping through the rovibrational lines can account for the equilibrium conditions found for CO and H$_2$O in the deeply-embedded starbursts. Here the CO bands probe regions with an intense local radiation field inside dusty young massive star clusters or near the most massive young stars. The lack of high-excitation molecular gas towards the AGN is attributed to geometric dilution of the intense radiation from the bright point source. An overview of the relevant excitation and radiative transfer physics is provided in an appendix., Comment: 28 pages, 17 figures. Accepted for publication in ApJ. This version includes small revisions following the referee report

ALMA observations revealed the presence of significant amounts of dust in the first Gyr of Cosmic time. However, the metal and dust buildup picture remains very uncertain due to the lack of constraints on metallicity. JWST has started to reveal the metal content of high-redshift targets, which may lead to firmer constraints on high-redshift dusty galaxies evolution. In this work, we use detailed chemical and dust evolution models to explore the evolution of galaxies within the ALMA REBELS survey, testing different metallicity scenarios that could be inferred from JWST observations. In the models, we track the buildup of stellar mass by using non-parametric SFHs for REBELS galaxies. Different scenarios for metal and dust evolution are simulated by allowing different prescriptions for gas flows and dust processes. The model outputs are compared with measured dust scaling relations, by employing metallicity-dependent calibrations for the gas mass based on the [CII]158micron line. Independently of the galaxies metal content, we found no need for extreme dust prescriptions to explain the dust masses revealed by ALMA. However, different levels of metal enrichment will lead to different dominant dust production mechanisms, with stardust production dominant over other ISM dust processes only in the metal-poor case. This points out how metallicity measurements from JWST will significantly improve our understanding of the dust buildup in high-redshift galaxies. We also show that models struggle to reproduce observables such as dust-to-gas and dust-to-stellar ratios simultaneously, possibly indicating an overestimation of the gas mass through current calibrations, especially at high metallicities., Comment: 16 pages + appendices, 9 Figures, 1 Table. Resubmitted to MNRAS after moderate revision

We present new JWST-NIRSpec IFS data for the luminous infrared galaxy NGC7469: a nearby (70.6Mpc) active galaxy with a Sy 1.5 nucleus that drives a highly ionized gas outflow and a prominent nuclear star-forming ring. Using the superb sensitivity and high spatial resolution of the JWST instrument NIRSpec-IFS, we investigate the role of the Seyfert nucleus in the excitation and dynamics of the circumnuclear gas. Our analysis focuses on the [Fe ii], H2, and hydrogen recombination lines that trace the radiation/shocked-excited molecular and ionized ISM around the AGN. We investigate the gas excitation through H2/Br{\gamma} and [Fe ii]/Pa\b{eta} emission line ratios and find that photoionization by the AGN dominates within the central 300 pc of the galaxy and together with a small region show ing signatures of shock-heated gas; these shock-heated regions are likely associated with a compact radio jet. In addition, the velocity field and velocity dispersion maps reveal complex gas kinematics. Rotation is the dominant feature, but we also identify non-circular motions consistent with gas inflows as traced by the velocity residuals and the spiral pattern in the Pa{\alpha} velocity dispersion map. The inflow is consistent with the mass outflow rate and two orders of magnitude higher than the AGN accretion rate. The compact nuclear radio jet has enough power to drive the highly ionized outflow. This scenario suggests that the inflow and outflow are in a self-regulating feeding-feedback process, with a contribution from the radio jet helping to drive the outflow., Comment: 15 pages, 4 figures, 3 tables, accepted for publication in ApJ

We present James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec) integral-field spectroscopy of the nearby luminous infrared galaxy, NGC 7469. We take advantage of the high spatial/spectral resolution and wavelength coverage of JWST /NIRSpec to study the 3.3 um neutral polycyclic aromatic hydrocarbon (PAH) grain emission on ~60 pc scales. We find a clear change in the average grain properties between the star-forming ring and the central AGN. Regions in the vicinity of the AGN, with [NeIII]/[NeII]>0.25, tend to have larger grain sizes and lower aliphatic-to-aromatic (3.4/3.3) ratios indicating that smaller grains are preferentially removed by photo-destruction in the vicinity of the AGN. We find an overall suppression of the total PAH emission relative to the ionized gas in the central 1 kpc region of the AGN in NGC 7469 compared to what has been observed with Spitzer on 3 kpc scales. However, the fractional 3.3 um to total PAH power is enhanced in the starburst ring, possibly due to a variety of physical effects on sub-kpc scales, including recurrent fluorescence of small grains or multiple photon absorption by large grains. Finally, the IFU data show that while the 3.3 um PAH-derived star formation rate (SFR) in the ring is 8% higher than that inferred from the [NeII] and [NeIII] emission lines, the integrated SFR derived from the 3.3 um feature would be underestimated by a factor of two due to the deficit of PAHs around the AGN, as might occur if a composite system like NGC 7469 were to be observed at high-redshift., Comment: 14 pages, 5 figures, 2 tables, Submitted to ApJL

Surface densities of gas, dust and stars provide a window into the physics of star-formation that, until the advent of high-resolution far-infrared/sub-millimeter observations, has been historically difficult to assess amongst dusty galaxies. To study the link between infrared (IR) surface densities and dust properties, we leverage the Atacama Large Millimetre/Submillimetre Array (ALMA) archive to measure the extent of cold dust emission in 15 $z\sim2$ IR selected galaxies selected on the basis of having available mid-IR spectroscopy from Spitzer. We use the mid-IR spectra to constrain the relative balance between dust heating from star-formation and active galactic nuclei (AGN), and to measure emission from Polycylic Aromatic Hydrocarbons (PAHs) -- small dust grains that play a key role in the photoelectric heating of gas. In general, we find that dust-obscured star-formation at high IR surface densities exhibits similar properties at low- and high-redshift, namely: local luminous IR galaxies have comparable PAH luminosity to total dust mass ratios as high-$z$ galaxies, and star-formation at $z\sim0-2$ is more efficient at high IR surface densities despite the fact that our sample of high$-z$ galaxies are closer to the main-sequence than local luminous IR galaxies. High star-formation efficiencies are coincident with a decline in the PAH/IR luminosity ratio reminiscent of the deficit observed in far-infrared fine-structure lines. Changes in the gas and dust conditions arising from high star-formation surface densities might help drive the star-formation efficiency up. This could help explain high efficiencies needed to reconcile star-formation and gas volume densities in dusty galaxies at cosmic noon., Comment: 18 pages, 10 figures, accepted to ApJ

Deep VLT/MUSE optical integral field spectroscopy has recently revealed an abundant population of ultra-faint galaxies ($M_{UV} \approx -15$; 0.01 $L_{\star}$) at $z=$2.9$-$6.7 due to their strong Lyman-$\alpha$ emission with no detectable continuum. The implied Lyman-$\alpha$ equivalent widths can be in excess of 100-200 Angstrom, challenging existing models of normal star formation and indicating extremely young ages, small stellar masses, and a very low amount of metal enrichment. We use JWST/NIRSpec's microshutter array to follow-up 45 of these galaxies (11h in G235M/F170LP and 7h in G395M/F290LP), as well as 45 lower-equivalent width Lyman-$\alpha$ emitters. Our spectroscopy covers the range 1.7$-$5.1 micron in order to target strong optical emission lines: H$\alpha$, [N II], [O III], and H$\beta$. Individual measurements as well as stacks reveal line ratios consistent with a metal poor nature (2$-$40% $Z_{\odot}$, depending on the calibration). The galaxies with the highest equivalent widths of Lyman-$\alpha$, in excess of 90 Angstrom, have lower [N II]/H$\alpha$ (1.9-$\sigma$) and [O III]/H$\beta$ (2.2-$\sigma$) ratios than those with lower equivalent widths, implying lower gas-phase metallicities at a combined significance of 2.4-$\sigma$. This implies a selection based on Lyman-$\alpha$ equivalent width is an efficient technique for identifying younger, less chemically enriched systems., Comment: 14 pages, 4 appendices; ApJ Accepted

We present new ALMA Band 8 (rest-frame $90\,\mu$m) observations of three massive ($M_\star \approx 10^{10}\,M_\odot$) galaxies at $z\approx7$ previously detected in [CII]$158\,\mu$m and underlying dust continuum emission in the Reionization Era Bright Emission Line Survey (REBELS). We detect the dust continuum emission of two of our targets in Band 8 (REBELS-25 and REBELS-38), while REBELS-12 remains undetected. Through modified blackbody fitting we determine cold dust temperatures ($T_\mathrm{dust} \approx 30 - 35\,$K) in both of the dual-band detected targets, given a fiducial model of optically thin emission with $\beta = 2.0$. Their dust temperatures are lower than most $z\sim7$ galaxies in the literature, and consequently their dust masses are higher ($M_\mathrm{dust} \approx 10^{8}\,M_\odot$). Nevertheless, these large dust masses are still consistent with predictions from models of dust production in the early Universe. In addition, we target and detect [OIII]$88\,\mu$m emission in both REBELS-12 and REBELS-25, and find $L_\mathrm{[OIII]} / L_\mathrm{[CII]}$ ratios of approximately unity, low compared to the $L_\mathrm{[OIII]} / L_\mathrm{[CII]} \gtrsim 2 - 10$ observed in the known $z\gtrsim6$ population thus far. We argue the lower line ratios are due to a comparatively weaker ionizing radiation field resulting from the less starbursty nature of our targets. This low burstiness supports the cold dust temperatures and below average $\mathrm{[OIII]}\lambda\lambda4959,5007 + \mathrm{H}\beta$ equivalent widths of REBELS-25 and REBELS-38, compared to the known high-redshift population. Overall, this provides evidence for the existence of a massive, dust-rich galaxy population at $z\approx7$ which has previously experienced vigorous star formation, but is currently forming stars in a steady, as opposed to bursty, manner., Comment: 18 pages + appendices; 8 figures in main text; Submitted to MNRAS on 27 Dec 2022

We present the results of a {\it James Webb Space Telescope} NIRCam investigation into the young massive star cluster (YMC) population in the luminous infrared galaxy VV 114. We identify 374 compact YMC candidates with a $S/N \geq 3$, 5, and 5 at F150W, F200W, and F356W respectively. A direct comparison with our {\it HST} cluster catalog reveals that $\sim 20\%$ of these sources are undetected at optical wavelengths. Based on {\it yggdrasil} stellar population models, we identify 17 YMC candidates in our {\it JWST} imaging alone with F150W-F200W and F200W-F356W colors suggesting they are all very young, dusty ($A_{V} = 5 - 15$), and massive ($10^{5.8} < M_{\odot} < 10^{6.1}$). The discovery of these `hidden' sources, many of which are found in the `overlap' region between the two nuclei, quadruples the number of $t < 3$ Myr clusters, and nearly doubles the number of $t < 6$ Myr clusters detected in VV 114. Now extending the cluster age distribution ($dN/d\tau \propto \tau^{\gamma}$) to the youngest ages, we find a slope of $\gamma = -1.30 \pm 0.39$ for $10^{6} < \tau (\mathrm{yr}) < 10^{7}$, which is consistent with the previously determined value from $10^{7} < \tau (\mathrm{yr}) < 10^{8.5}$, and confirms that VV 114 has a steep age distribution slope for all massive star clusters across the entire range of cluster ages observed. Finally, the consistency between our {\it JWST}- and {\it HST}-derived age distribution slopes indicates that the balance between cluster formation and destruction has not been significantly altered in VV 114 over the last 0.5 Gyr., Comment: 13 pages, 4 figures, 1 table, published in the Astrophysical Journal

Our current understanding of the cosmic star formation history at z>3 is primarily based on UV-selected galaxies (i.e., LBGs). Recent studies of H-dropouts have revealed that we may be missing a large proportion of star formation that is taking place in massive galaxies at z>3. In this work, we extend the H-dropout criterion to lower masses to select optically dark/faint galaxies (OFGs), in order to complete the census between LBGs and H-dropouts. Our criterion (H> 26.5 mag & [4.5] < 25 mag) combined with a de-blending technique is designed to select not only extremely dust-obscured massive galaxies but also normal star-forming galaxies. In total, we identified 27 OFGs at z_phot > 3 (z_med=4.1) in the GOODS-ALMA field, covering a wide distribution of stellar masses with log($M_{\star}$/$M_{\odot}$) = 9.4-11.1. We find that up to 75% of the OFGs with log($M_{\star}$/$M_{\odot}$) = 9.5-10.5 were neglected by previous LBGs and H-dropout selection techniques. After performing stacking analyses, the OFGs exhibit shorter gas depletion timescales, slightly lower gas fractions, and lower dust temperatures than typical star-forming galaxies. Their SFR_tot (SFR_ IR+SFR_UV) is much larger than SFR_UVcorr (corrected for dust extinction), with SFR_tot/SFR_UVcorr = $8\pm1$, suggesting the presence of hidden dust regions in the OFGs that absorb all UV photons. The average dust size measured by a circular Gaussian model fit is R_e(1.13 mm)=1.01$\pm$0.05 kpc. We find that the cosmic SFRD at z>3 contributed by massive OFGs is at least two orders of magnitude higher than the one contributed by equivalently massive LBGs. Finally, we calculate the combined contribution of OFGs and LBGs to the cosmic SFRD at z=4-5 to be 4 $\times$ 10$^{-2}$ $M_{\odot}$ yr$^{-1}$Mpc$^{-3}$, which is about 0.15 dex (43%) higher than the SFRD derived from UV-selected samples alone at the same redshift., Comment: 26 pages, 14 figures, 8 tables. Accepted for publication in A&A

Observations and simulations of interacting galaxies and mergers in the local universe have shown that interactions can significantly enhance the star formation rates (SFR) and fueling of Active Galactic Nuclei (AGN). However, at higher redshift, some simulations suggest that the level of star formation enhancement induced by interactions is lower due to the higher gas fractions and already increased SFRs in these galaxies. To test this, we measure the SFR enhancement in a total of 2351 (1327) massive ($M_*>10^{10}M_\odot$) major ($1

We present James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) integral-field spectroscopy of the nearby merging, luminous infrared galaxy, NGC 7469. This galaxy hosts a Seyfert type-1.5 nucleus, a highly ionized outflow, and a bright, circumnuclear star-forming ring, making it an ideal target to study AGN feedback in the local Universe. We take advantage of the high spatial/spectral resolution of JWST/MIRI to isolate the star-forming regions surrounding the central active nucleus and study the properties of the dust and warm molecular gas on ~100 pc scales. The starburst ring exhibits prominent Polycyclic Aromatic Hydrocarbon (PAH) emission, with grain sizes and ionization states varying by only ~30%, and a total star formation rate of $\rm 10 - 30 \ M_\odot$/yr derived from fine structure and recombination emission lines. Using pure rotational lines of H2, we detect 1.2$\times$10$^{7} \rm \ M_\odot$ of warm molecular gas at a temperature higher than 200 K in the ring. All PAH bands get significantly weaker towards the central source, where larger and possibly more ionized grains dominate the emission. However, the bulk of the dust and molecular gas in the ring appears unaffected by the ionizing radiation or the outflowing wind from the AGN. These observations highlight the power of JWST to probe the inner regions of dusty, rapidly evolving galaxies for signatures of feedback and inform models that seek to explain the co-evolution of supermassive black holes and their hosts., Comment: 13 pages, 4 figures, 2 tables, Submitted to ApJL

We present James Webb Space Telescope (JWST) imaging of NGC 7469 with the Near-Infrared Camera (NIRCam) and the Mid-InfraRed Instrument (MIRI). NGC 7469 is a nearby, $z=0.01627$, luminous infrared galaxy (LIRG) that hosts both a Seyfert Type-1.5 nucleus and a circumnuclear starburst ring with a radius of $\sim$0.5 kpc. The new near-infrared (NIR) JWST imaging reveals 66 star-forming regions, 37 of which were not detected by HST observations. Twenty-eight of the 37 sources have very red NIR colors that indicate obscurations up to A$_{\rm{v}}\sim7$ and a contribution of at least 25$\%$ from hot dust emission to the 4.4$\mu$m band. Their NIR colors are also consistent with young ($<$5 Myr) stellar populations and more than half of them are coincident with the MIR emission peaks. These younger, dusty star-forming regions account for $\sim$6$\%$ and $\sim$17$\%$ of the total 1.5$\mu$m and 4.4$\mu$m luminosity of the starburst ring, respectively. Thanks to JWST, we find a significant number of young dusty sources that were previously unseen due to dust extinction. The newly identified 28 young sources are a significant increase compared to the number of HST-detected young sources (4-5). This makes the total percentage of the young population rise from $\sim$15$\%$ to 48$\%$. These results illustrate the effectiveness of JWST in identifying and characterizing previously hidden star formation in the densest star-forming environments around AGN., Comment: 4 figures, 1 table, Accepted by ApJL

The nearby, luminous infrared galaxy (LIRG) NGC 7469 hosts a Seyfert nucleus with a circumnuclear star-forming ring and is thus the ideal local laboratory for investigating the starburst--AGN connection in detail. We present integral-field observations of the central 1.3 kpc region in NGC 7469 obtained with the JWST Mid-InfraRed Instrument. Molecular and ionized gas distributions and kinematics at a resolution of {\sim}100 pc over the 4.9 - 7.6{\mu}m region are examined to study gas dynamics influenced by the central AGN. The low-ionization [Fe II] {\lambda}5.34{\mu}m and [Ar II] {\lambda}6.99{\mu}m lines are bright on the nucleus and in the starburst ring, as opposed to H2 S(5) {\lambda}6.91{\mu}m which is strongly peaked at the center and surrounding ISM. The high-ionization [Mg V] line is resolved and shows a broad, blueshifted component associated with the outflow. It has a nearly face-on geometry that is strongly peaked on the nucleus, where it reaches a maximum velocity of -650 km/s, and extends about 400 pc to the East. Regions of enhanced velocity dispersion in H2 and [Fe II] {\sim}180 pc from the AGN that also show high L(H2)/L(PAH) and L([Fe II])/L(Pf{\alpha}) ratios to the W and N of the nucleus pinpoint regions where the ionized outflow is depositing energy, via shocks, into the dense interstellar medium between the nucleus and the starburst ring. These resolved mid-infrared observations of the nuclear gas dynamics demonstrate the power of JWST and its high-sensitivity integral-field spectroscopic capability to resolve feedback processes around supermassive black holes in the dusty cores of nearby LIRGs., Comment: 14 pages, 6 figures, 1 table; accepted for publication in ApJL

James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI) images of the luminous infrared (IR) galaxy VV 114 are presented. This redshift ~ 0.020 merger has a western component (VV 114W) rich in optical star clusters and an eastern component (VV 114E) hosting a luminous mid-IR nucleus hidden at UV and optical wavelengths by dust lanes. With MIRI, the VV 114E nucleus resolves primarily into bright NE and SW cores separated by 630 pc. This nucleus comprises 45% of the 15um light of VV 114, with the NE and SW cores having IR luminosities, L_ IR (8-1000um) ~ 8+/-0.8x10^10 L_sun and ~ 5+/-0.5x10^10 L_sun, respectively, and IR densities, Sigma_IR >~ 2+/-0.2x10^13 L_sun / kpc^2 and >~ 7+/-0.7x10^12 L_sun / kpc^2, respectively -- in the range of Sigma_IR for the Orion star-forming core and the nuclei of Arp 220. The NE core, previously speculated to have an Active Galactic Nucleus (AGN), has starburst-like mid-IR colors. In contrast, the VV 114E SW has AGN-like colors. Approximately 40 star-forming knots with L_IR ~ 0.02-5x10^10 L_sun are identified, 25% of which have no optical counterpart. Finally, diffuse emission accounts for 40-60% of the mid-IR emission. Mostly notably, filamentary Poly-cyclic Aromatic Hydrocarbon (PAH) emission stochastically excited by UV and optical photons accounts for half of the 7.7um light of VV 114. This study illustrates the ability of JWST to detect obscured compact activity and distributed PAH emission in the most extreme starburst galaxies in the local Universe., Comment: 8 pages, 4 figures, 1 table, Submitted to ApJL

We have used the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST) to obtain the first spatially resolved, mid-infrared (mid-IR) images of IIZw096, a merging luminous infrared galaxy (LIRG) at $z = 0.036$. Previous observations with the Spitzer Space Telescope suggested that the vast majority of the total IR luminosity (LIR) of the system originated from a small region outside of the two merging nuclei. New observations with JWST/MIRI now allow an accurate measurement of the location and luminosity density of the source that is responsible for the bulk of the IR emission. We estimate that 40-70% of the IR bolometric luminosity, or $3-5 \times 10^{11}\,{\rm{L_{\odot}}}$, arises from a source no larger than 175pc in radius, suggesting a luminosity density of at least $3-5 \times 10^{12} \, {\rm{L_{\odot} \, kpc^{-2}}}$. In addition, we detect 11 other star forming sources, five of which were previously unknown. The MIRI F1500W/F560W colors of most of these sources, including the source responsible for the bulk of the far-IR emission, are much redder than the nuclei of local LIRGs. These observations reveal the power of JWST to disentangle the complex regions at the hearts of merging, dusty galaxies., Comment: 3 figures, 1 table, Accepted for publication in ApJL

Cosmic dust is an essential component shaping both the evolution of galaxies and their observational signatures. How quickly dust builds up in the early Universe remains an open question that requires deep observations at (sub-)millimeter wavelengths to resolve. Here we use Atacama Large Millimeter Array observations of 45 galaxies from the Reionization Era Bright Emission Line Survey (REBELS) and its pilot programs, designed to target [CII] and dust emission in UV-selected galaxies at $z\sim7$, to investigate the dust content of high-redshift galaxies through a stacking analysis. We find that the typical fraction of obscured star formation $f_\mathrm{obs} = \mathrm{SFR}_\mathrm{IR} / \mathrm{SFR}_\mathrm{UV + IR}$ depends on stellar mass, similar to what is observed at lower redshift, and ranges from $f_\mathrm{obs} \approx 0.3 - 0.6$ for galaxies with $\log_{10}\left(M_\star / M_\odot\right) = 9.4 - 10.4$. We further adopt the $z\sim7$ stellar mass function from the literature to extract the obscured cosmic star formation rate density (SFRD) from the REBELS survey. Our results suggest only a modest decrease in the SFRD between $3\lesssim z \lesssim 7$, with dust-obscured star formation still contributing $\sim30\%$ at $z\sim7$. While we extensively discuss potential caveats, our analysis highlights the continued importance of dust-obscured star formation even well into the epoch of reionization., Comment: MNRAS Accepted (01 Nov 2022)

We present the results of a stacking analysis performed on Spitzer/Infrared Spectrograph high-resolution mid-infrared spectra of luminous infrared galaxies (LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). By binning on mid-infrared active galactic nucleus (AGN) fraction and stacking spectra, we detect bright emission lines [Ne II] and [Ne III], which trace star formation, and fainter emission lines [Ne V] and [O IV], which trace AGN activity, throughout the sample. We find the [Ne II] luminosity is fairly constant across all AGN fraction bins, while the [O IV] and [Ne V] luminosities increase by over an order of magnitude. Our measured average line ratios, [Ne V]/[Ne II] and [O IV]/[Ne II], at low AGN fraction are similar to H II galaxies while the line ratios at high AGN fraction are similar to LINERs and Seyferts. We decompose the [O IV] luminosity into star-formation and AGN components by fitting the [O IV] luminosity as a function of the [Ne II] luminosity and the mid-infrared AGN fraction. The [O IV] luminosity in LIRGs is dominated by star formation for mid-infrared AGN fractions $\lesssim0.3$. With the corrected [O IV] luminosity, we calculate black hole accretion rates ranging from $10^{-5}$ M$_{\odot}$/yr at low AGN fractions to 0.2 M$_{\odot}$/yr at the highest AGN fractions. We find that using the [O IV] luminosity, without correcting for star formation, can lead to an overestimate of the BHAR by up to a factor of 30 in starburst dominated LIRGs. Finally, we show the BHAR/SFR ratio increases by more than three orders of magnitude as a function of mid-infrared AGN fraction in LIRGs., Comment: 17 pages, 15 figures. Accepted to ApJ

We report 18 dust continuum detections ($\geq 3.3\sigma$) at $\sim88{\rm \mu m}$ and $158{\rm \mu m}$ out of 49 ultraviolet(UV)-bright galaxies ($M_{\rm UV} < -21.3$ mag) at $z>6.5$, observed by the Cycle-7 ALMA Large Program, REBELS and its pilot programs. This has more than tripled the number of dust continuum detections known at $z>6.5$. Out of these 18 detections, 12 are reported for the first time as part of REBELS. In addition, 15 of the dust continuum detected galaxies also show a [CII]$_{\rm 158{\rm \mu m}}$ emission line, providing us with accurate redshifts. We anticipate more line emission detections from six targets (including three continuum detected targets) where observations are still ongoing. The dust continuum detected sources in our sample tend to have a redder UV spectral slope than the ones without a dust continuum detection. We estimate that all of the sources have an infrared (IR) luminosity ($L_{\rm IR}$) in a range of $3-8 \times 10^{11} L_\odot$, except for one with $L_{\rm IR} = 1.5^{+0.8}_{-0.5} \times 10^{12}\,L_{\odot}$. Their fraction of obscured star formation is significant at $\gtrsim 50\%$. Some of the dust continuum detected galaxies show spatial offsets ($\sim 0.5-1.5''$) between the rest-UV and far-IR emission peaks. These separations appear to have an increasing trend against an indicator that suggests spatially decoupled phases of obscured and unobscured star formation. REBELS offers the best available statistical constraints on obscured star formation in UV-bright, massive galaxies at $z > 6.5$., Comment: 17 pages, 9 figures, 4 tables, Accepted for publication in MNRAS

We present specific star-formation rates for 40 UV-bright galaxies at $z\sim7-8$ observed as part of the Reionization Era Bright Emission Line Survey (REBELS) ALMA large program. The sSFRs are derived using improved measures of SFR and stellar masses, made possible by measurements of far-infrared (FIR) continuum emission and [CII]-based spectroscopic redshifts. For each source in the sample, we derive stellar masses from SED fitting and total SFRs from calibrations of the UV and FIR emission. The median sSFR is $18_{-5}^{+7}$ Gyr$^{-1}$, significantly larger than literature measurements lacking constraints in the FIR. The increase in sSFR reflects the larger obscured SFRs we derive from the dust continuum relative to that implied by the UV+optical SED. We suggest that such differences may reflect spatial variations in dust across these luminous galaxies, with the component dominating the FIR distinct from that dominating the UV. We demonstrate that the inferred stellar masses (and hence sSFRs) are strongly-dependent on the assumed star formation history (SFH) in reionization-era galaxies. When large sSFR galaxies are modeled with non-parametric SFHs, the derived stellar masses can increase by an order of magnitude relative to constant star formation models, owing to the presence of a significant old stellar population that is outshined by the recent burst. The [CII] line widths in the largest sSFR systems are often very broad, suggesting dynamical masses that are easily able to accommodate the dominant old stellar population suggested by non-parametric models. Regardless of these systematic uncertainties in the derived parameters, we find that the sSFR increases rapidly toward higher redshifts for massive galaxies ($9.6<\log(\rm M_*/M_{\odot})<9.8$), with a power law that goes as $(1+z)^{1.7\pm0.3}$, broadly consistent with expectations from the evolving baryon accretion rates., Comment: 18 pages, 10 figures, submitted to MNRAS

Recent work has shown that UV-luminous reionization-era galaxies often exhibit strong Lyman-alpha emission despite being situated at redshifts where the IGM is thought to be substantially neutral. It has been argued that this enhanced Ly$\alpha$ transmission reflects the presence of massive galaxies in overdense regions which power large ionized bubbles. An alternative explanation is that massive galaxies shift more of their Ly$\alpha$ profile to large velocities (relative to the systemic redshift) where the IGM damping wing absorption is reduced. Such a mass-dependent trend is seen at lower redshifts, but whether one exists at $z\sim7$ remains unclear owing to the small number of existing systemic redshift measurements in the reionization era. This is now changing with the emergence of [CII]-based redshifts from ALMA. Here we report MMT/Binospec Ly$\alpha$ spectroscopy of eight UV-bright ($\mathrm{M_{UV}}^{}\sim-22$) galaxies at $z\simeq7$ selected from the ALMA REBELS survey. We detect Ly$\alpha$ in 4 of 8 galaxies and use the [CII] systemic redshifts to investigate the Ly$\alpha$ velocity profiles. The Ly$\alpha$ lines are significantly redshifted from systemic (average velocity offset=223 km/s) and broad (FWHM$\approx$300$-$650 km/s), with two sources showing emission extending to $\approx$750 km/s. We find that the broadest Ly$\alpha$ profiles are associated with the largest [CII] line widths, suggesting a potential link between the Ly$\alpha$ FWHM and the dynamical mass. Since Ly$\alpha$ photons at high velocities transmit efficiently through the $z=7$ IGM, our data suggest that velocity profiles play a significant role in boosting the Ly$\alpha$ visibility of the most UV-luminous reionization-era galaxies., Comment: 17 pages, 7 figures, submitted to MNRAS

We present the results of a Hubble Space Telescope WFC3 near-UV and ACS/WFC optical study into the star cluster populations of 10 luminous and ultra-luminous infrared galaxies (U/LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). Through integrated broadband photometry we have derived ages, masses, and extinctions for a total of 1027 star clusters in galaxies with $d_{L} <$ 110 Mpc in order to avoid issues related to cluster blending. The measured cluster age distribution slope of $dN/d\tau \propto \tau^{-0.5 +/- 0.2}$ is steeper than what has been observed in lower-luminosity star-forming galaxies. Further, differences in the slope of the observed cluster age distribution between inner- ($dN/d\tau \propto \tau^{-1.07 +/- 0.12}$) and outer-disk ($dN/d\tau \propto \tau^{-0.37 +/- 0.09}$) star clusters provides evidence of mass-dependent cluster destruction in the central regions of LIRGs driven primarily by the combined effect of strong tidal shocks and encounters with massive GMCs. Excluding the nuclear ring surrounding the Seyfert 1 nucleus in NGC 7469, the derived cluster mass function (CMF: $dN/dM \propto M^{\alpha}$) has marginal evidence for a truncation in the power-law (PL) at $M_{t} \sim 2$x$10^{6} M_{\odot}$ for our three most cluster-rich galaxies, which are all classified as early-stage mergers. Finally, we find evidence of a flattening of the CMF slope of $dN/dM \propto M^{-1.42 \pm 0.1}$ for clusters in late-stage mergers relative to early-stage ($\alpha = -1.65 \pm 0.02$), which we attribute to an increase in the formation of massive clusters over the course of the interaction., Comment: 70 pages, 21 figures, 5 tables, accepted for publication in ApJ

We study the molecular gas content of 24 star-forming galaxies at $z=3-4$, with a median stellar mass of $10^{9.1}$ M$_{\odot}$, from the MUSE Hubble Ultra Deep Field (HUDF) Survey. Selected by their Lyman-alpha-emission and H-band magnitude, the galaxies show an average EW $\approx 20$ angstrom, below the typical selection threshold for Lyman Alpha Emitters (EW $> 25$ angstrom), and a rest-frame UV spectrum similar to Lyman Break Galaxies. We use rest-frame optical spectroscopy from KMOS and MOSFIRE, and the UV features observed with MUSE, to determine the systemic redshifts, which are offset from Lyman alpha by 346 km s$^{-1}$, with a 100 to 600 km s$^{-1}$ range. Stacking CO(4-3) and [CI](1-0) (and higher-$J$ CO lines) from the ALMA Spectroscopic Survey of the HUDF (ASPECS), we determine $3\sigma$ upper limits on the line luminosities of $4.0\times10^{8}$ K km s$^{-1}$pc$^{2}$ and $5.6\times10^{8}$ K km s$^{-1}$pc$^{2}$, respectively (for a 300 km s$^{-1}$ linewidth). Stacking the 1.2 mm and 3 mm dust continuum flux densities, we find a $3\sigma$ upper limits of 9 $\mu$Jy and $1.2$ $\mu$Jy, respectively. The inferred gas fractions, under the assumption of a 'Galactic' CO-to-H$_{2}$ conversion factor and gas-to-dust ratio, are in tension with previously determined scaling relations. This implies a substantially higher $\alpha_{\rm CO} \ge 10$ and $\delta_{\rm GDR} \ge 1200$, consistent with the sub-solar metallicity estimated for these galaxies ($12 + \log(O/H) \approx 7.8 \pm 0.2$). The low metallicity of $z \ge 3$ star-forming galaxies may thus make it very challenging to unveil their cold gas through CO or dust emission, warranting further exploration of alternative tracers, such as [CII]., Comment: 26 pages, 12 figures, 3 tables. Accepted for publication in ApJ

The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Band 6 scan (212-272 GHz) covers potential [CII] emission in galaxies at $6\leq z \leq8$ throughout a 2.9 arcmin$^2$ area. By selecting on known Lyman-$\alpha$ emitters (LAEs) and photometric dropout galaxies in the field, we perform targeted searches down to a 5$\sigma$ [CII] luminosity depth $L_{\mathrm{[CII]}}\sim2.0\times10^8$ L$_{\odot}$, corresponding roughly to star formation rates (SFRs) of $10$-$20$ M$_{\odot}$ yr$^{-1}$ when applying a locally calibrated conversion for star-forming galaxies, yielding zero detections. While the majority of galaxies in this sample are characterized by lower SFRs, the resulting upper limits on [CII] luminosity in these sources are consistent with the current literature sample of targeted ALMA observations of $z=6$-$7$ LAEs and Lyman-break galaxies (LBGs), as well as the locally calibrated relations between $L_{\mathrm{[CII]}}$ and SFR -- with the exception of a single [CII]-deficient, UV luminous LBG. We also perform a blind search for [CII]-bright galaxies that may have been missed by optical selections, resulting in an upper limit on the cumulative number density of [CII] sources with $L_{\mathrm{[CII]}}>2.0\times10^8$ L$_{\odot}$ ($5\sigma $) to be less than $1.8\times10^{-4}$ Mpc$^{-3}$ (90% confidence level). At this luminosity depth and volume coverage, we present an observed evolution of the [CII] luminosity function from $z=6$-$8$ to $z\sim0$ by comparing the ASPECS measurement to literature results at lower redshift., Comment: 16 pages, 6 figures, accepted for publication in ApJ

Star formation in galaxies is regulated by the heating and cooling in the interstellar medium. In particular, the processing of molecular gas into stars will depend strongly on the ratio of gas heating to gas cooling in the neutral gas around sites of recent star-formation. In this work, we combine mid-infrared (mid-IR) observations of Polycyclic Aromatic Hydrocarbons (PAHs), the dominant heating mechanism of gas in the interstellar medium (ISM), with [C II], [O I], and [Si II] fine-structure emission, the strongest cooling channels in dense, neutral gas. The ratio of IR cooling line emission to PAH emission measures the photoelectric efficiency, a property of the ISM which dictates how much energy carried by ultraviolet photons gets transferred into the gas. We find that star-forming, IR luminous galaxies in the Great Observatories All-Sky LIRG Survey (GOALS) with high IR surface densities have low photoelectric efficiencies. These systems also have, on average, higher ratios of radiation field strength to gas densities, and larger average dust grain size distributions. The data support a scenario in which the most compact galaxies have more young star-forming regions per unit area, which exhibit less efficient gas heating. These conditions may be more common at high-z, and may help explain the higher star-formation rates at cosmic noon. We make predictions on how this can be investigated with JWST., Comment: 17 pages, 8 figures, 2 tables, accepted for publication in ApJ

We report the discovery of a 10 comoving Mpc-scale structure traced by massive submillimeter galaxies (SMGs) at z~4.6. These galaxies are selected from an emission line search of ALMA Band 7 observations targeting 184 luminous submillimeter sources ($S_{850\mu{\rm m}}\geq$ 6.2 mJy) across 1.6 degrees$^2$ in the COSMOS field. We identify four [CII] emitting SMGs and two probable [CII] emitting SMG candidates at z=4.60-4.64 with velocity-integrated signal-to-noise ratio of SNR>8. Four of the six emitters are near-infrared blank SMGs. After excluding one SMG whose emission line is falling at the edge of the spectral window, all galaxies show clear velocity gradients along the major axes that are consistent with rotating gas disks. The estimated rotation velocities of the disks are 330-550 km s$^{-1}$ and the inferred host dark-matter halo masses are ~2-8 $\times$ 10$^{12}$M$_{\odot}$. From their estimated halo masses and [CII] luminosity function, we suggest that these galaxies have a high (50-100%) duty cycle and high (~0.1) baryon conversion efficiency (SFR relative to baryon accretion rate), and that they contribute $\simeq$2% to the total star-formation rate density at z=4.6. These SMGs are concentrated within just 0.3% of the full survey volume, suggesting they are strongly clustered. The extent of this structure and the individual halo masses suggest that these SMGs will likely evolve into members of a ~10$^{15}$M$_{\odot}$ cluster at z=0. This survey reveals synchronized dusty starburst in massive halos at z>4, which could be driven by mergers or fed by smooth gas accretion.

Galaxy interactions and mergers are thought to play an important role in the evolution of galaxies. Studies in the nearby universe show a higher AGN fraction in interacting and merging galaxies than their isolated counterparts, indicating that such interactions are important contributors to black hole growth. To investigate the evolution of this role at higher redshifts, we have compiled the largest known sample of major spectroscopic galaxy pairs (2381 with $\Delta V <5000$ km s$^{-1}$) at $0.5

We present a CO and atomic fine-structure line luminosity function analysis using the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). ASPECS consists of two spatially-overlapping mosaics that cover the entire ALMA 3mm and 1.2mm bands. We combine the results of a line candidate search of the 1.2mm data cube with those previously obtained from the 3mm cube. Our analysis shows that $\sim$80% of the line flux observed at 3mm arises from CO(2-1) or CO(3-2) emitters at $z$=1-3 (`cosmic noon'). At 1.2mm, more than half of the line flux arises from intermediate-J CO transitions ($J_{\rm up}$=3-6); $\sim12$% from neutral carbon lines; and $< 1$% from singly-ionized carbon, [CII]. This implies that future [CII] intensity mapping surveys in the epoch of reionization will need to account for a highly significant CO foreground. The CO luminosity functions probed at 1.2mm show a decrease in the number density at a given line luminosity (in units of $L'$) at increasing $J_{\rm up}$ and redshift. Comparisons between the CO luminosity functions for different CO transitions at a fixed redshift reveal sub-thermal conditions on average in galaxies up to $z\sim 4$. In addition, the comparison of the CO luminosity functions for the same transition at different redshifts reveals that the evolution is not driven by excitation. The cosmic density of molecular gas in galaxies, $\rho_{\rm H2}$, shows a redshift evolution with an increase from high redshift up to $z\sim1.5$ followed by a factor $\sim 6$ drop down to the present day. This is in qualitative agreement with the evolution of the cosmic star-formation rate density, suggesting that the molecular gas depletion time is approximately constant with redshift, after averaging over the star-forming galaxy population., Comment: 23 pages, 11 figures, 5 tables. Accepted for publication in the Astrophysical Journal

We report molecular gas mass estimates obtained from a stacking analysis of CO line emission in the ALMA Spectroscopic Survey (ASPECS) using the spectroscopic redshifts from the optical integral field spectroscopic survey by the Multi Unit Spectroscopic Explorer (MUSE) of the {\it Hubble} Ultra Deep Field (HUDF). Stacking was performed on subsets of the sample of galaxies classified by their stellar mass and position relative to the main-sequence relation (on, above, below). Among all the CO emission lines, from \cotwoone to CO(6-5), with redshifts accessible via the ASPECS Band~3 and the MUSE data, \cotwoone provides the strongest constraints on the molecular gas content. We detect \cotwoone emission in galaxies down to stellar masses of $\log{(M_*/M_\odot)}=10.0$. Below this stellar mass, we present a new constraint on the molecular gas content of $z\sim1.5$ main-sequence galaxies by stacking based on the MUSE detections. We find that the molecular gas mass of main-sequence galaxies continuously decreases with stellar mass down to $\log{(M_*/M_\odot)}\approx9.0$. Assuming a metallicity-based CO--to--$\rm H_2$ conversion factor, the molecular gas-to-stellar mass ratio from $\log{(M_*/M_\odot)}\sim9.0$ to $\sim10.0$ does not seem to decrease as fast as for $\log{(M_*/M_\odot)}>10.0$, which is in line with simulations and studies at lower redshift. The inferred molecular gas density $\rho{\rm (H_2)}=(0.49\pm0.09)\times10^8\,{\rm M_\odot\,Mpc^{-3}}$ of MUSE-selected galaxies at $z\sim1.5$ is comparable with the one derived in the HUDF with a different CO selection. Using the MUSE data we recover most of the CO emission in our deep ALMA observations through stacking, demonstrating the synergy between volumetric surveys obtained at different wavebands., Comment: 26 pages, 12 figures, 6 tables accepted for publication in ApJ

We make use of sensitive (9.3 microJy/beam RMS) 1.2mm-continuum observations from the ASPECS ALMA large program of the Hubble Ultra Deep Field (HUDF) to probe dust-enshrouded star formation from 1362 Lyman-break galaxies spanning the redshift range z=1.5-10 (to ~7-28 Msolar/yr at 4 sigma over the entire range). We find that the fraction of ALMA-detected galaxies in our z=1.5-10 samples increases steeply with stellar mass, with the detection fraction rising from 0% at 10^9 Msolar to 85(-18)(+9)% at >10^{10} Msolar. Moreover, stacking all 1253 low-mass (<10^{9.25} Msolar) galaxies over the ASPECS footprint, we find a mean continuum flux of -0.1+/-0.4 microJy/beam, implying a hard upper limit on the obscured SFR of <0.6 Msolar/yr (4 sigma) in a typical low-mass galaxy. The correlation between the infrared excess IRX of UV-selected galaxies (L(IR)/L(UV)) and the UV-continuum slope is also seen in our ASPECS data and shows consistency with a Calzetti-like relation at >10^{9.5} M_{solar} and a SMC-like relation at lower masses. Using stellar-mass and beta measurements for z~2 galaxies over CANDELS, we derive a new empirical relation between beta and stellar mass and then use this correlation to show that our IRX-beta and IRX-stellar mass relations are consistent with each other. We then use these constraints to express the infrared excess as a bivariate function of beta and stellar mass. Finally, we present updated estimates of star-formation rate density determinations at z>3, leveraging current improvements in the measured infrared excess and recent probes of ultra-luminous far-IR galaxies at z>2., Comment: 32 pages, 22 figures, 11 pages, accepted for publication in the Astrophysical Journal, key results shown in Figures 7, 8, 11, and 20, updated to match version in press

We investigate the CO excitation and interstellar medium (ISM) conditions in a cold gas mass-selected sample of 22 star-forming galaxies at $z=0.46-3.60$, observed as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). Combined with VLA follow-up observations, we detect a total of 34 CO $J \rightarrow J-1$ transitions with $J=1$ up to $8$ (and an additional 21 upper limits, up to $J=10$) and six [C I] ${^3P}_1 \rightarrow {^3P}_0$ and ${^3P}_2 \rightarrow {^3P}_1$ transitions (and 12 upper limits). The CO(2-1) and CO(3-2)-selected galaxies, at $z=1.2$ and $2.5$, respectively, exhibit a range in excitation in their mid-$J=4,5$ and high-$J=7,8$ lines, on average lower than ($L_{\rm IR}$-brighter) BzK-color- and submillimeter-selected galaxies at similar redshifts. The former implies that a warm ISM component is not necessarily prevalent in gas mass-selected galaxies at $z=1.2$. We use stacking and Large Velocity Gradient models to measure and predict the average CO ladders at $z<2$ and $z\geq2$, finding $r_{21}=0.75 \pm 0.11$ and $r_{31}=0.77 \pm 0.14$, respectively. From the models, we infer that the galaxies at $z\geq2$ have intrinsically higher excitation than those at $z<2$. This fits a picture in which the global excitation is driven by an increase in the star formation rate surface density of galaxies with redshift. We derive a neutral atomic carbon abundance of $(1.9 \pm 0.4) \times 10^{-5}$, comparable to the Milky Way and main-sequence galaxies at similar redshifts, and fairly high densities ($\geq 10^4$ cm$^{-3}$), consistent with the low-$J$ CO excitation. Our results imply a decrease in the cosmic molecular gas mass density at $z\geq2$ compared to previous ASPECS measurements., Comment: 37 pages, 15 figures, 6 tables. Accepted for publication in the Astrophysical Journal

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.

We investigate the ultraviolet (UV) spectral properties of faint Lyman-$\alpha$ emitters (LAEs) in the redshift range 2.93. We compute various averaged spectra of LAEs sub-sampled on the basis of their observational (e.g., Ly$\alpha$ strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We search for UV spectral features other than Ly$\alpha$, such as higher-ionization nebular emission lines and absorption features. We successfully observe the OIII]1666 and CIII]909 collisionally excited emission lines and the HeII1640 recombination feature, as well as the resonant CIV1550 doublet either in emission or P-Cygni. We compare the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Ly$\alpha$ equivalent widths, fainter UV magnitudes, bluer UV spectral slopes and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z>2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass and faint LAEs. We publicly release the stacked spectra as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope., Comment: 26 pages, 18 figures. Accepted for publication in A&A

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 submillimeter 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., Comment: ApJ published on Sept. 24, 2020. Contains 36 pages, 17 figures. Final version matching the proof reviewed article

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

The Lya emitter (LAE) fraction, X_LAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of X_LAE are still debated. Thanks to deep data obtained with MUSE, we can measure the evolution of X_LAE homogeneously over a wide redshift range of z~3-6 for UV-faint galaxies (down to M_1500~-17.75). This is significantly fainter than in former studies, and allows us to probe the bulk of the population of high-z star-forming galaxies. We construct a UV-complete photo-redshift sample following UV luminosity functions and measure the Lya emission with MUSE using the second data release from the MUSE HUDF Survey. We derive the redshift evolution of X_LAE for M_1500 in [-21.75;-17.75] for the first time with a equivalent width range EW(Lya)>=65 A and find low values of X_ LAE<~30% at z<~6. For M_1500 in [-20.25;-18.75] and EW(Lya)<~25 A, our X_LAE values are consistent with those in the literature within 1sigma at z<~5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of X_LAE on M_1500 for EW(Lya)>~50 A at z~3-4, in contrast with previous work. The differences in X_LAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lya emission, hence X_LAE is biased towards higher values. Our results suggest either a lower increase of X_LAE towards z~6 than previously suggested, or even a turnover of X_LAE at z~5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed X_LAE much better than models where SF is a smooth function of time., Comment: 25 pages, 17 figures, 5 tables, accepted for publication in A&A

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic Reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z~4-5 continuum-faint Lyman-alpha emitters from deep VLT/MUSE spectroscopy and directly measure their Halpha emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than -16 (~0.01 L_star), implying a median Lyman-alpha equivalent width of 249 Angstroms. By combining the Halpha measurement with the UV magnitude we determine the ionizing photon production efficiency, xi_ion, a first for such faint galaxies. The measurement of log (xi_ion [Hz/erg]) = 26.28 (+0.28; -0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower-luminosity, Lyman-alpha-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4e-4 and 0.008, with light-weighted ages less than 3 Myr., Comment: 12 Pages, 6 Figures; Accepted for publication in MNRAS

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., Comment: 32 pages, 13 figures, Accepted for publication in ApJ

In this paper we present our imaging observations on the CO(7-6) line and its underlying continuum emission of the young submillimeter galaxy LESS 073 at redshift 4.755, using the Atacama Large Millimeter/submillimeter Array (ALMA). At the achieved resolution of $\sim$$1^{\prime\prime}.2\times0^{\prime\prime}.9$ ($8\times6$~kpc$^2$), the CO(7-6) emission is largely unresolved (with a deconvolved size of $1^{\prime\prime}.1(\pm0^{\prime\prime}.5) \times 0^{\prime\prime}.9(\pm0^{\prime\prime}.8)$), and the continuum emission is totally unresolved. The CO(7-6) line emission has an integrated flux of $0.86\pm0.08$~Jy km/s, and a line width of $343\pm40$ km/s. The continuum emission has a flux density of 0.51 mJy. By fitting the observed far-infrared (FIR) spectral energy distribution of LESS 073 with a single-temperature modified blackbody function, we obtained a dust temperature $T_{\rm dust}=57.6\pm3.5$ K, 60-to-100 $\mu$m flux density ratio $f_{60}/f_{100}=0.86\pm0.08$, and total infrared luminosity $L_{\rm IR}=(5.8\pm0.9) \times 10^{12}~L_\odot$. The SED-fit-based $f_{60}/f_{100}$ is consistent with those estimated from various line ratios as advocated by our earlier work, indicating that those proposed line-ratio-based method can be used to practically derive $f_{60}/f_{100}$ for high-$z$ sources. The total molecular gas mass of LESS 073 is $(3.3\pm1.7) \times10^{10}~M_\odot$, and the inferred gas depletion time is about 43 Myr., Comment: 8 pages, 6 figures, accepted for publication in ApJ

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., Comment: 31 pages, 8 Figures, accepted for publication in ApJ

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)., Comment: 22 pages, 11 figures. Paper re-submitted to ApJ after addressing the first round of comments by the referee

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., Comment: Accepted for publication in ApJ. Only minor differences respect to previous version

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., Comment: 31 pages, 20 figures, 5 tables. Resubmitted to ApJ after addressing the first round of comments by the referee

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., Comment: Accepted for publication in ApJ. 28 Pages and 13 Figures

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., Comment: 7 pages, 2 figures, Science White paper submitted to Astro2020 Decadal Survey

We study the physical properties of the ionized gas in local disks using the sample of 38 nearby $\sim10^{8.5-11.2}$M$_\odot$ Star-Forming Main Sequence (SFMS) galaxies observed so far as part of the MUSE Atlas of Disks (MAD). Specifically, we use all strong emission lines in the MUSE wavelength range 4650-9300 \AA\ to investigate the resolved ionized gas properties on $\sim$100 pc scales. This spatial resolution enables us to disentangle HII regions from the Diffuse Ionized Gas (DIG) in the computation of gas metallicities and Star Formation Rates (SFRs) of star forming regions. The gas metallicities generally decrease with radius. The metallicity of the HII regions is on average $\sim$0.1 dex higher than that of the DIG, but the metallicity radial gradient in both components is similar. The mean metallicities within the inner galaxy cores correlate with the total stellar mass of the galaxies. On our <100 pc scales, we find two correlations previously reported at kpc scales: a spatially resolved Mass-Metallicity Relation (RMZR) and a spatially resolved SFMS (RSFMS). We find no secondary dependency of the RMZR with the SFR density. We find that both resolved relations have a local origin, as they do not depend on the total stellar mass. The observational results of this paper are consistent with the inside-out scenario for the growth of galactic disks., Comment: 21 pages + appendices. Accepted for publication in MNRAS