Search

After over three decades of unsuccessful attempts, we report the first detection of molecular gas emission in Malin 1, the largest spiral galaxy observed to date, and one of the most iconic giant low surface brightness galaxies. Using ALMA, we detect significant $^{12}$CO(J=1-0) emission in the galaxy's central region and tentatively identify CO emission across three regions on the disc. These observations allow for a better estimate of the H$_2$ mass and molecular gas mass surface density, both of which are remarkably low given the galaxy's scale. By integrating data on its HI mass, we derive a very low molecular-to-atomic gas mass ratio. Overall, our results highlight the minimal presence of molecular gas in Malin 1, contrasting sharply with its extensive, homogeneous atomic gas reservoir. For the first time, we position Malin 1 on the Kennicutt-Schmidt (K-S) diagram, where it falls below the main sequence for normal spirals, consistent with previous upper limits but now with more accurate figures. These findings are crucial for constraining our understanding of star formation processes in environments characterized by extremely low molecular gas densities and for refining models of galaxy formation, thereby improving predictions concerning the formation, evolution, and distribution of these giant, elusive galaxies., Comment: Accepted for publication in The Astrophysical Journal Letters. Published

Using a combination of HST, JWST, and ALMA data, we perform spatially resolved spectral energy distributions (SED) fitting of fourteen 4

We investigate the spatial extent and structure of the [CII] line emission in a sample of 34 galaxies at $z=4-6$ from the ALMA-CRISTAL Survey. By modeling the [CII] line emission in the interferometric visibility, we derive the effective radius of [CII] line emission assuming an exponential profile. The [CII] line radius ranges from 0.5 to 3.5 kpc with an average value of $\langle R_{e,[CII]}\rangle=1.90$ kpc. We compare the [CII] sizes with the sizes of rest-frame UV and FIR continua, which were measured from the HST F160W images and ALMA Band-7 continuum images, respectively. We confirm that the [CII] line emission is more spatially extended than the continuum emission, with average size ratios of $\langle R_{e,[CII]}/R_{e,UV}\rangle=2.90$ and $\langle R_{e,[CII]}/R_{e,FIR}\rangle=1.54$, although about half of the FIR-detected sample show comparable spatial extent between [CII] line and FIR continuum emission. The residual visibility of the best-fit model do not show statistical evidence of flux excess, indicating that the [CII] line emission in star-forming galaxies can be characterized by an extended exponential profile. Overall, our results suggest that the spatial extent of [CII] line emission can primarily be explained by PDRs associated with star formation activity, while the contribution from diffuse neutral medium (atomic gas) and the effects of past merger events may further expand the [CII] line distributions, causing their variations. Finally, we report the negative correlation between $\Sigma_{[CII]}$ and EW$_{Ly\alpha}$, and possible negative correlation between $R_{e,[CII]}/R_{e,UV}$ and EW$_{Ly\alpha}$, which may be in line with the scenario that atomic gas largely contributes to the extended [CII] line emission. Future 3-D analysis of Ly$\alpha$ and H$\alpha$ lines will shed light on the association of the extended [CII] line emission with atomic gas and outflows., Comment: Resubmitted to A&A (23 pages, 14 figures, 5 tables) with some updates on the analyses

We present results from Atacama Large Millimeter/submillimeter Array (ALMA) spectral line-scan observations at 3-mm and 2-mm bands of three near-infrared-dark (NIR-dark) galaxies behind two massive lensing clusters MACS J0417.5-1154 and RXC J0032.1+1808. Each of these three sources is a faint (de-lensed $S_{\text{1.2 mm}}$ $<$ 1 mJy) triply lensed system originally discovered in the ALMA Lensing Cluster Survey. We have successfully detected CO and [C I] emission lines and confirmed that their spectroscopic redshifts are $z=3.652$, 2.391, and 2.985. By utilizing a rich multi-wavelength data set, we find that the NIR-dark galaxies are located on the star formation main sequence in the intrinsic stellar mass range of log ($M_*$/$M_\odot$) = 9.8 - 10.4, which is about one order of magnitude lower than that of typical submillimeter galaxies (SMGs). These NIR-dark galaxies show a variety in gas depletion times and spatial extent of dust emission. One of the three is a normal star-forming galaxy with gas depletion time consistent with a scaling relation, and its infrared surface brightness is an order of magnitude smaller than that of typical SMGs. Since this galaxy has an elongated axis ratio of $\sim 0.17$, we argue that normal star-forming galaxies in an edge-on configuration can be heavily dust-obscured. This implies that existing deep WFC3/F160W surveys may miss a fraction of typical star-forming main-sequence galaxies due to their edge-on orientation., Comment: 23 pages, 10 figures, Submitted to ApJ

We present the morphological parameters and global properties of dust-obscured star formation in typical star-forming galaxies at z=4-6. Among 26 galaxies composed of 20 galaxies observed by the Cycle-8 ALMA Large Program, CRISTAL, and six galaxies from archival data, we have individually detected rest-frame 158$\mu$m dust continuum emission from 19 galaxies, nine of which are reported for the first time. The derived far-infrared luminosities are in the range $\log_{10} L_{\rm IR}\,[L_{\odot}]=$10.9-12.4, an order of magnitude lower than previously detected massive dusty star-forming galaxies (DSFGs). The average relationship between the fraction of dust-obscured star formation ($f_{\rm obs}$) and the stellar mass is consistent with previous results at z=4-6 in a mass range of $\log_{10} M_{\ast}\,[M_{\odot}]\sim$9.5-11.0 and show potential evolution from z=6-9. The individual $f_{\rm obs}$ exhibits a significant diversity, and it shows a correlation with the spatial offset between the dust and the UV continuum, suggesting the inhomogeneous dust reddening may cause the source-to-source scatter in $f_{\rm obs}$. The effective radii of the dust emission are on average $\sim$1.5 kpc and are $\sim2$ times more extended than the rest-frame UV. The infrared surface densities of these galaxies ($\Sigma_{\rm IR}\sim2.0\times10^{10}\,L_{\odot}\,{\rm kpc}^{-2}$) are one order of magnitude lower than those of DSFGs that host compact central starbursts. On the basis of the comparable contribution of dust-obscured and dust-unobscured star formation along with their similar spatial extent, we suggest that typical star-forming galaxies at z=4-6 form stars throughout the entirety of their disks.

We develop new tools for continuum and spectral stacking of ALMA data, and apply these to the ALMA Lensing Cluster Survey (ALCS). We derive average dust masses, gas masses and star formation rates (SFR) from the stacked observed 260~GHz continuum of 3402 individually undetected star-forming galaxies, of which 1450 are cluster galaxies and 1952 field galaxies, over three redshift and stellar mass bins (over $z = 0$-1.6 and log $M_{*} [M_{\odot}] = 8$-11.7), and derive the average molecular gas content by stacking the emission line spectra in a SFR-selected subsample. The average SFRs and specific SFRs of both cluster and field galaxies are lower than those expected for Main Sequence (MS) star-forming galaxies, and only galaxies with stellar mass of log $M_{*} [M_{\odot}] = 9.35$-10.6 show dust and gas fractions comparable to those in the MS. The ALMA-traced average `highly obscured' SFRs are typically lower than the SFRs observed from optical to near-IR spectral analysis. Cluster and field galaxies show similar trends in their contents of dust and gas, even when field galaxies were brighter in the stacked maps. From spectral stacking we find a potential CO ($J=4\to3$) line emission (SNR $\sim4$) when stacking cluster and field galaxies with the highest SFRs., Comment: 13 pages, 8 figures, 1 table (+4 pages, 4 figures in Appendix). Accepted for publication in MNRAS

Although radar and communications signal classification are usually treated separately, they share similar characteristics, and methods applied in one domain can be potentially applied in the other. We propose a simple and unified scheme for the classification of radar and communications signals using Long Short-Term Memory (LSTM) neural networks. This proposal provides an improvement of the state of the art on radar signals where LSTM models are starting to be applied within schemes of higher complexity. To date, there is no standard public dataset for radar signals. Therefore, we propose DeepRadar2022, a radar dataset used in our systematic evaluations that is available publicly and will facilitate a standard comparison between methods., Comment: This paper was submitted to the Radar Conference 2023

We present a statistical study of 180 dust continuum sources identified in 33 massive cluster fields by the ALMA Lensing Cluster Survey (ALCS) over a total of 133 arcmin$^{2}$ area, homogeneously observed at 1.2 mm. ALCS enables us to detect extremely faint mm sources by lensing magnification, including near-infrared (NIR) dark objects showing no counterparts in existing {\it Hubble Space Telescope} and {\it Spitzer} images. The dust continuum sources belong to a blind sample ($N=141$) with S/N $\gtrsim$ 5.0 (a purity of $>$ 0.99) or a secondary sample ($N=39$) with S/N= $4.0-5.0$ screened by priors. With the blind sample, we securely derive 1.2-mm number counts down to $\sim7$ $\mu$Jy, and find that the total integrated 1.2mm flux is 20.7$^{+8.5}_{-6.5}$ Jy deg$^{-2}$, resolving $\simeq$ 80 % of the cosmic infrared background light. The resolved fraction varies by a factor of $0.6-1.1$ due to the completeness correction depending on the spatial size of the mm emission. We also derive infrared (IR) luminosity functions (LFs) at $z=0.6-7.5$ with the $1/V_{\rm max}$ method, finding the redshift evolution of IR LFs characterized by positive luminosity and negative density evolution. The total (=UV+IR) cosmic star-formation rate density (SFRD) at $z>4$ is estimated to be $161^{+25}_{-21}$ % of the established measurements, which were almost exclusively based on optical$-$NIR surveys. Although our general understanding of the cosmic SFRD is unlikely to change beyond a factor of 2, these results add to the weight of evidence for an additional ($\approx 60$ %) SFRD component contributed by the faint-mm population, including NIR dark objects., Comment: ApJS in press, 34 pages, 8 Tables, and 18 Figures (+38 pages, 5 Tables, and 6 Figures in Appendix)

We present ALMA deep spectroscopy for a lensed galaxy at $z_{\rm spec}=8.496$ with $\log(M_{\rm star}/M_{\odot})\sim7.8$ whose optical nebular lines and stellar continuum are detected by JWST/NIRSpec and NIRCam Early Release Observations in SMACS0723. Our ALMA spectrum shows [OIII]88$\mu$m and [CII]158$\mu$m line detections at $4.0\sigma$ and $4.5\sigma$, respectively. The redshift and position of the [OIII] line coincide with those of the JWST source, while the [CII] line is blue-shifted by 90 km s$^{-1}$ with a spatial offset of $0.''5$ ($\approx0.5$ kpc in source plane) from the JWST source. The NIRCam F444W image, including [OIII]$\lambda$5007 and H$\beta$ line emission, spatially extends beyond the stellar components by a factor of $>8$. This indicates that the $z=8.5$ galaxy has already experienced strong outflows whose oxygen and carbon produce the extended [OIII]$\lambda$5007 and the offset [CII] emission, which would promote ionizing photon escape and facilitate reionization. With careful slit-loss corrections and removals of emission spatially outside the galaxy, we evaluate the [OIII]88$\mu$m/$\lambda$5007 line ratio, and derive the electron density $n_{\rm e}$ by photoionization modeling to be $220^{+170}_{-100}$ cm$^{-3}$, which is comparable with those of $z\sim2-3$ galaxies. We estimate an [OIII]88$\mu$m/[CII]158$\mu$m line ratio in the galaxy of $>4$, as high as those of known $z\sim6-9$ galaxies. This high [OIII]88$\mu$m/[CII]158$\mu$m line ratio is generally explained by the high $n_{\rm e}$ as well as the low metallicity ($Z_{\rm gas}/Z_{\odot}=0.04^{+0.02}_{-0.02}$), high ionization parameter ($\log U > -2.27$), and low carbon-to-oxygen abundance ratio ($\log$(C/O) $=[-0.52:-0.24]$) obtained from the JWST/NIRSpec data; further [CII] follow-up observations will constrain the covering fraction of photodissociation regions., Comment: 27 pages, 11 figures, 5 tables. ApJ in press

Spiral arms serve crucial purposes in star formation and galaxy evolution. In this paper, we report the identification of A2744-DSG-$z3$, a dusty, multiarm spiral galaxy at $z=3.059$ using the James Webb Space Telescope (JWST) NIRISS imaging and grism spectroscopy. A2744-DSG-$z3$ was discovered as a gravitationally lensed sub-millimeter galaxy with ALMA. This is the most distant stellar spiral structure seen thus far, consistent with cosmological simulations which suggest $z\approx3$ as the epoch when spirals emerge. Thanks to the gravitational lensing and excellent spatial resolution of JWST, the spiral arms are resolved with a spatial resolution of $\approx290$\,pc. Based on SED fitting, the spiral galaxy has a de-lensed star formation rate of $85\pm30 \ M_{\odot}$ yr$^{-1}$, and a stellar mass of $\approx10^{10.6}\ M_{\odot}$, indicating that A2744-DSG-$z3$ is a main-sequence galaxy. After fitting the spiral arms, we find a stellar effective radius ($R_{e, \rm{star}}$) of $5.0\pm1.5$ kpc. Combing with ALMA measurements, we find that the effective radii ratio between dust and stars is $\approx0.4$, similar to {those} of massive SFGs at $z\sim2$, indicating a compact dusty core in A2744-DSG-$z3$. Moreover, this galaxy appears to be living in a group environment: including A2744-DSG-$z3$, at least three galaxies at $z=3.05 - 3.06$ {are} spectroscopically confirmed by JWST/NIRISS and ALMA, residing within a lensing-corrected projected scale of $\approx 70$ kpc. This, along with the asymmetric brightness profile, further suggests that the spiral arms may be triggered by minor merger events at $z\gtrsim3$., Comment: Main text 9 pages, four figures in the main text, and three figures in the appendix; accepted for publication by ApJL

Spatially extended halos of H I Ly$\alpha$ emission are now ubiquitously found around high-redshift star-forming galaxies. But our understanding of the nature and powering mechanisms of these halos is still hampered by the complex radiative transfer effects of the Ly$\alpha$ line and limited angular resolution. In this paper, we present resolved Multi Unit Spectroscopic Explorer (MUSE) observations of SGAS J122651.3+215220, a strongly-lensed pair of $L^{*}$ galaxies at $z=2.92$ embedded in a Ly$\alpha$ halo of $L_{Ly\alpha}=(6.2\pm1.3)\times10^{42}$ erg s$^{-1}$. Globally, the system shows a line profile that is markedly asymmetric and redshifted, but its width and peak shift vary significantly across the halo. By fitting the spatially binned Ly$\alpha$ spectra with a collection of radiative transfer galactic wind models, we infer a mean outflow expansion velocity of $\approx 211$ km s$^{-1}$, with higher values preferentially found on both sides of the system's major axis. The velocity of the outflow is validated with the blueshift of low-ionization metal absorption lines in the spectra of the central galaxies. We also identify a faint ($M_{1500} \approx -16.7$) companion detected in both Ly$\alpha$ and the continuum, whose properties are in agreement with a predicted population of satellite galaxies that contribute to the extended Ly$\alpha$ emission. Finally, we briefly discuss the impact of the interaction between the central galaxies on the properties of the halo and the possibility of in situ fluorescent Ly$\alpha$ production., Comment: Published in ApJ. 26 pages, 14 figures

We present an ALMA-Herschel joint analysis of sources detected by the ALMA Lensing Cluster Survey (ALCS) at 1.15 mm. Herschel/PACS and SPIRE data at 100-500 $\mu$m are deblended for 180 ALMA sources in 33 lensing cluster fields that are either detected securely (141 sources; in our main sample) or tentatively at S/N$\geq$4 with cross-matched HST/Spitzer counterparts, down to a delensed 1.15-mm flux density of $\sim0.02$ mJy. We performed far-infrared spectral energy distribution modeling and derived the physical properties of dusty star formation for 125 sources (109 independently) that are detected at $>2\sigma$ in at least one Herschel band. 27 secure ALCS sources are not detected in any Herschel bands, including 17 optical/near-IR-dark sources that likely reside at $z=4.2\pm1.2$. The 16-50-84 percentiles of the redshift distribution are 1.15-2.08-3.59 for ALCS sources in the main sample, suggesting an increasing fraction of $z\simeq1-2$ galaxies among fainter millimeter sources ($f_{1150}\sim 0.1$ mJy). With a median lensing magnification factor of $\mu = 2.6_{-0.8}^{+2.6}$, ALCS sources in the main sample exhibit a median intrinsic star-formation rate of $94_{-54}^{+84}\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$, lower than that of conventional submillimeter galaxies at similar redshifts by a factor of $\sim$3. Our study suggests weak or no redshift evolution of dust temperature with $L_\mathrm{IR}<10^{12}\,\mathrm{L}_\odot$ galaxies within our sample at $z \simeq 0 - 2$. At $L_\mathrm{IR}>10^{12}\,\mathrm{L}_\odot$, the dust temperatures show no evolution across $z \simeq 1 -4$ while being lower than those in the local Universe. For the highest-redshift source in our sample ($z=6.07$), we can rule out an extreme dust temperature ($>$80 K) that was reported for MACS0416 Y1 at $z=8.31$., Comment: 44 pages, 23 figures, 4 tables, accepted for publication in ApJ

We report spectroscopic confirmations of 15 Lyman-alpha galaxies at $z\sim7$, implying a spectroscopic confirmation rate of $\sim$80% on candidates selected from LAGER (Lyman-Alpha Galaxies in the Epoch of Reionization), which is the largest (24 deg$^2$) survey aimed at finding Lyman-alpha emitters (LAEs) at $z\sim7$ using deep narrow-band imaging from DECam at CTIO. LAEs at high-redshifts are sensitive probes of cosmic reionization and narrow-band imaging is a robust and effective method for selecting a large number of LAEs. In this work, we present results from the spectroscopic follow-up of LAE candidates in two LAGER fields, COSMOS and WIDE-12, using observations from Keck/LRIS. We report the successful detection of Ly$\alpha$ emission in 15 candidates (11 in COSMOS and 4 in WIDE-12 fields). Three of these in COSMOS have matching confirmations from a previous LAGER spectroscopic follow-up and are part of the overdense region, LAGER-$z7$OD1. Additionally, two candidates that were not detected in the LRIS observations have prior spectroscopic confirmations from Magellan. Including these, we obtain a spectroscopic confirmation success rate of $\sim$$80$% for LAGER LAE candidates. Apart from Ly$\alpha$, we do not detect any other UV nebular lines in our LRIS spectra; however, we estimate a 2$\sigma$ upper limit for the ratio of NV/Ly$\alpha$, $f_{NV}/f_{Ly\alpha} \lesssim 0.27$, which implies that ionizing emission from these sources is mostly dominated by star formation. Including confirmations from this work, a total of 33 LAE sources from LAGER are now spectroscopically confirmed. LAGER has more than doubled the sample of spectroscopically confirmed LAE sources at $z\sim7$., Comment: 8 pages, 4 figures; Submitted to ApJ

We present an Atacama Large Millimeter/submillimeter Array (ALMA) survey of CO(4--3) line emitting galaxies in 17 quasar fields at $z\sim4$ aimed to perform the first systematic search of dusty galaxies in high$-z$ quasar environments. Our blind search of galaxies around the quasars results in 5 CO emitters with S/N$\geq5.6$ within a projected radius of $R\lesssim1.5\,h^{-1}\,$cMpc and a velocity range of $\rm\Delta\,v=\pm1000\,$km\,s$^{-1}$ around the quasar. In blank fields, we expect to detect only 0.28 CO emitters within the same volume, implying a total overdensity of $17.6^{+11.9}_{-7.6}$ in our fields, and indicating that quasars trace massive structures in the early universe. We quantify this overdensity by measuring the small-scale clustering of CO emitters around quasars, resulting in a cross-correlation length of $r_{\rm 0,QG}=8.37^{+2.42}_{-2.04}\,h^{-1}\,$cMpc, assuming a fixed slope $\gamma=1.8$. This contradicts the reported mild overdensities (x1.4) of Lyman alpha emitters (LAEs) in the same fields at scales of $R\lesssim7\,h^{-1}\,$cMpc which is well described by a cross-correlation length 3 times lower than that measured for CO emitters. We discuss some possibilities to explain this discrepancy, including low star formation efficiency, and excess of dust in galaxies around quasars. Finally, we constrain, for the first time, the clustering of CO emitters at $z\sim4$, finding an auto-correlation length of $r_{\rm 0,CO}=3.14\pm1.71\,h^{-1}\,$cMpc (with $\gamma=1.8$). Our work together with the previous study of LAEs around quasars traces simultaneously the clustering properties of both optical and dusty galaxy populations in quasars fields, stressing the importance of multi-wavelength studies, and highlighting important questions about galaxy properties in high$-z$ dense environments., Comment: 17 pages, 8 figures, submitted to the 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

We present a new measurement of the Ly$\alpha$ luminosity function at redshift $z=6.9$, finding moderate evolution from $z=5.7$ that is consistent with a fully or largely ionized $z\sim7$ intergalactic medium. Our result is based on four fields of the LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) project. Our survey volume of $6.1\times10^{6}$ Mpc$^{3}$ is double that of the next largest $z\sim 7$ survey. We combine two new LAGER fields (WIDE12 and GAMA15A) with two previously reported LAGER fields (COSMOS and CDFS). In the new fields, we identify $N=95$ new $z=6.9$ Ly$\alpha$ emitters (LAEs); characterize our survey's completeness and reliability; and compute Ly$\alpha$ luminosity functions. The best-fit Schechter luminosity function parameters for all four LAGER fields are in good general agreement. Two fields (COSMOS and WIDE12) show evidence for a bright-end excess above the Schechter function fit. We find that the Ly$\alpha$ luminosity density declines at the same rate as the UV continuum LF from $z=5.7$ to $z=6.9$. This is consistent with an intergalactic medium that was fully ionized as early as redshift $z\sim 7$, or with a volume-averaged neutral hydrogen fraction of $x_{HI} < 0.33$ at $1\sigma$., Comment: 19 pages, 11 figures; accepted by ApJ

Hot, dust-obscured galaxies (Hot DOGs) are a population of hyper-luminous obscured quasars identified by WISE. We present ALMA observations of the [CII] fine-structure line and underlying dust continuum emission in a sample of seven of the most extremely luminous (EL; L$_{\rm bol}$ $\ge$ 10$^{14}$ L$_\odot$) Hot DOGs, at redshifts z ~ 3.0-4.6. The [CII] line is robustly detected in four objects, tentatively in one, and likely red-shifted out of the spectral window in the remaining two based on additional data. On average, [CII] is red-shifted by ~ 780 km/s from rest-frame ultraviolet emission lines. EL Hot DOGs exhibit consistently very high ionized gas surface densities, with $\Sigma_{\rm [CII]}$ ~ 1-2 x 10$^{9}$ L$_\odot$ kpc$^{-2}$; as high as the most extreme cases seen in other high-redshift quasars. As a population, EL Hot DOG hosts seem to be roughly centered on the main-sequence of star forming galaxies, but the uncertainties are substantial and individual sources can fall above and below. The average, intrinsic [CII] and dust continuum sizes (FWHMs) are ~ 2.1 kpc and ~ 1.6 kpc, respectively, with a very narrow range of line-to-continuum size ratios, 1.61 $\pm$ 0.10, suggesting they could be linearly proportional. The [CII] velocity fields of EL Hot DOGs are diverse: from barely rotating structures, to resolved hosts with ordered, circular motions, to complex, disturbed systems that are likely the result of ongoing mergers. In contrast, all sources display large line-velocity dispersions, FWHM $\gtrsim$ 500 km/s, which on average are larger than optically and IR-selected quasars at similar or higher redshifts. We argue that one possible hypothesis for the lack of a common velocity structure, the systematically large dispersion of the ionized gas, and the presence of nearby companion galaxies may be that, rather than a single event, the EL Hot DOG phase could be recurrent., Comment: 31 pages, 13 figures. Accepted for publication in A&A

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

Protoclusters, the progenitors of the most massive structures in the Universe, have been identified at redshifts of up to 6.6. Besides exploring early structure formation, searching for protoclusters at even higher redshifts is particularly useful to probe the reionization. Here we report the discovery of the protocluster LAGER-z7OD1 at a redshift of 6.93, when the Universe was only 770 million years old and could be experiencing rapid evolution of the neutral hydrogen fraction in the intergalactic medium. The protocluster is identified by an overdensity of 6 times the average galaxy density, and with 21 narrowband selected Lyman-$\alpha$ galaxies, among which 16 have been spectroscopically confirmed. At redshifts similar to or above this record, smaller protogroups with fewer members have been reported. LAGER-z7OD1 shows an elongated shape and consists of two subprotoclusters, which would have merged into one massive cluster with a present-day mass of $3.7 \times 10^{15}$ solar masses. The total volume of the ionized bubbles generated by its member galaxies is found to be comparable to the volume of the protocluster itself, indicating that we are witnessing the merging of the individual bubbles and that the intergalactic medium within the protocluster is almost fully ionized. LAGER-z7OD1 thus provides a unique natural laboratory to investigate the reionization process., Comment: 9 pages, 5 figures, and 1 table. Accepted by Nature Astronomy

We present bright [CII] 158 $\mu$m line detections from a strongly magnified and multiply-imaged ($\mu\sim20-160$) sub-$L^{*}$ ($M_{\rm UV}$ = $-19.75^{+0.55}_{-0.44}$) Lyman-break galaxy (LBG) at $z=6.0719\pm0.0004$ from the ALMA Lensing Cluster Survey (ALCS). Emission lines are identified at 268.7 GHz at $\geq$ 8$\sigma$ exactly at positions of two multiple images of the LBG behind the massive galaxy cluster RXCJ0600$-$2007. Our lens models, updated with the latest spectroscopy from VLT/MUSE, indicate that a sub region of the LBG crosses the caustic and is lensed into a long ($\sim6''$) arc with a local magnification of $\mu\sim 160$, for which the [CII] line is also significantly detected. The source-plane reconstruction resolves the interstellar medium (ISM) structure, showing that the [CII] line is co-spatial with the rest-frame UV continuum at the scale of $\sim$300 pc. The [CII] line properties suggest that the LBG is a rotation-dominated system whose velocity gradient explains a slight difference of redshifts between the whole LBG and its sub region. The star formation rate (SFR)-$L_{\rm [CII]}$ relations from the sub to the whole regions of the LBG are consistent with those of local galaxies. We evaluate the lower limit of the faint-end of the [CII] luminosity function at $z=6$, and find that it is consistent with predictions from semi analytical models and from the local SFR-$L_{\rm [CII]}$ relation with a SFR function at $z=6$. These results imply that the local SFR-$L_{\rm [CII]}$ relation is universal for a wide range of scales including the spatially resolved ISM, the whole region of galaxy, and the cosmic scale, even in the epoch of reionization., Comment: 27 pages, 10 figures, 6 tables, ApJ in press

We combine the recent determination of the evolution of the cosmic density of molecular gas (H_2) using deep, volumetric surveys, with previous estimates of the cosmic density of stellar mass, star formation rate and atomic gas (HI), to constrain the evolution of baryons associated with galaxies averaged over cosmic time and space. The cosmic HI and H_2 densities are roughly equal at z~1.5. The H_2 density then decreases by a factor 6^{+3}_{-2} to today's value, whereas the HI density stays approximately constant. The stellar mass density is increasing continuously with time and surpasses that of the total gas density (HI and H_2) at redshift z~1.5. The growth in stellar mass cannot be accounted for by the decrease in cosmic H_2 density, necessitating significant accretion of additional gas onto galaxies. With the new H_2 constraints, we postulate and put observational constraints on a two step gas accretion process: (i) a net infall of ionized gas from the intergalactic/circumgalactic medium to refuel the extended HI reservoirs, and (ii) a net inflow of HI and subsequent conversion to H_2 in the galaxy centers. Both the infall and inflow rate densities have decreased by almost an order of magnitude since z~2. Assuming that the current trends continue, the cosmic molecular gas density will further decrease by about a factor of two over the next 5 Gyr, the stellar mass will increase by approximately 10%, and cosmic star formation activity will decline steadily toward zero, as the gas infall and accretion shut down., Comment: Accepted for publication in The Astrophysical Journal

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 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)

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

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

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., Comment: 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 used the Plateau De Bure Interferometer to observe multiple CO and neutral carbon transitions in a z=2.2 main sequence disk galaxy, BX610. Our observation of CO(7-6), CO(4-3), and both far-infrared(FIR) [CI] lines complements previous observations of H$\alpha$ and low-J CO, and reveals a galaxy that is vigorously forming stars with UV fields (Log($G$ G$_0^{-1}) \lesssim3.25);$ although less vigorously than local ultra-luminous infrared galaxies or most starbursting submillimeter galaxies in the early universe. Our observations allow new independent estimates of the cold gas mass which indicate $M_\textrm{gas}\sim2\times10^{11}$M$_\odot$, and suggest a modestly larger $\alpha_{\textrm{CO}}$ value of $\sim$8.2. The corresponding gas depletion timescale is $\sim$1.5 Gyr. In addition to gas of modest density (Log($n$ cm$^3)\lesssim3$ ) heated by star formation, BX610 shows evidence for a significant second gas component responsible for the strong high-J CO emission. This second component might either be a high-density molecular gas component heated by star formation in a typical photodissociation region, or could be molecular gas excited by low-velocity C shocks. The CO(7-6)-to-FIR luminosity ratio we observe is significantly higher than typical star-forming galaxies and suggests that CO(7-6) is not a reliable star-formation tracer in this galaxy., Comment: 9 pages, 6 figures, 2 tables, accepted for publication in Astronomy & Astrophysics. Data from the spectra in Fig. 1 are available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/

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

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., Comment: Re-submitted to ApJ after addressing the first round of comments by the referee, other comments welcome

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

We analyze the >4-sigma sources in the most sensitive 100 arcmin^2 area (rms <0.56 mJy) of a SCUBA-2 850 micron survey of the GOODS-S and present the 75 band 7 ALMA sources (>4.5-sigma) obtained from high-resolution interferometric follow-up observations. The SCUBA-2---and hence ALMA---samples should be complete to 2.25 mJy. Of the 53 SCUBA-2 sources in this complete sample, only five have no ALMA detections, while 13% (68% confidence range 7-19%) have multiple ALMA counterparts. Color-based high-redshift dusty galaxy selection techniques find at most 55% of the total ALMA sample. In addition to using literature spectroscopic and optical/NIR photometric redshifts, we estimate FIR photometric redshifts based on an Arp 220 template. We identify seven z>4 candidates. We see the expected decline with redshift of the 4.5 micron and 24 micron to 850 micron flux ratios, confirming these as good diagnostics of z>4 candidates. We visually classify 52 ALMA sources, finding 44% (68% confidence range 35-53%) to be apparent mergers. We calculate rest-frame 2-8 keV and 8-28 keV luminosities using the 7 Ms Chandra X-ray image. Nearly all of the ALMA sources detected at 0.5-2 keV are consistent with a known X-ray luminosity to 850 micron flux relation for star-forming galaxies, while most of those detected at 2-7 keV are moderate luminosity AGNs that lie just above the 2-7 keV detection threshold. The latter largely have substantial obscurations of log N_H = 23-24 cm^-2, but two of the high-redshift candidates may even be Compton thick., Comment: 44 pages, 35 figures. Submitted to ApJ

Binary neutron-star mergers (BNSMs) are among the most readily detectable gravitational-wave (GW) sources with LIGO. They are also thought to produce short $\gamma$-ray bursts (SGRBs), and kilonovae that are powered by r-process nuclei. Detecting these phenomena simultaneously would provide an unprecedented view of the physics during and after the merger of two compact objects. Such a Rosetta Stone event was detected by LIGO/Virgo on 17 August 2017 at a distance of $\sim 44$ Mpc. We monitored the position of the BNSM with ALMA at 338.5 GHz and GMRT at 1.4 GHz, from 1.4 to 44 days after the merger. Our observations rule out any afterglow more luminous than $3\times 10^{26}~{\rm erg\,s}^{-1}\,{\rm Hz}^{-1}$ in these bands, probing $>$2--4 dex fainter than previous SGRB limits. We match these limits, in conjunction with public data announcing the appearance of X-ray and radio emission in the weeks after the GW event, to templates of off-axis afterglows. Our broadband modeling suggests that GW170817 was accompanied by a SGRB and that the GRB jet, powered by $E_{\rm AG,\,iso}\sim10^{50}$~erg, had a half-opening angle of $\sim20^\circ$, and was misaligned by $\sim41^\circ$ from our line of sight. The data are also consistent with a more collimated jet: $E_{\rm AG,\,iso}\sim10^{51}$~erg, $\theta_{1/2,\,\rm jet}\sim5^\circ$, $\theta_{\rm obs}\sim17^\circ$. This is the most conclusive detection of an off-axis GRB afterglow and the first associated with a BNSM-GW event to date. Assuming a uniform top-hat jet, we use the viewing angle estimates to infer the initial bulk Lorentz factor and true energy release of the burst., Comment: Accepted for publication in the Astrophysical Journal Letter. 12 pages, 4 figures

We present Atacama Large Millimeter/submillimeter Array observations of CO lines and dust continuum emission of the source RCSGA 032727--132609, a young $z=1.7$ low-metallicity starburst galaxy. The CO(3-2) and CO(6-5) lines, and continuum at rest-frame $450\,\mu m$ are detected and show a resolved structure in the image plane. We use the corresponding lensing model to obtain a source plane reconstruction of the detected emissions revealing intrinsic flux density of $S_{450\,\mu m}=23.5_{-8.1}^{+26.8}$ $\mu$Jy and intrinsic CO luminosities $L'_{\rm CO(3-2)}=2.90_{-0.23}^{+0.21}\times10^{8}$ ${\rm K\,km\,s^{-1}\,pc^{2}}$ and $L'_{\rm CO(6-5)}=8.0_{-1.3}^{+1.4}\times10^{7}$ ${\rm K\,km\,s^{-1}\,pc^{2}}$. We used the resolved properties in the source plane to obtain molecular gas and star-formation rate surface densities of $\Sigma_{\rm H2}=16.2_{-3.5}^{+5.8}\,{\rm M}_{\odot}\,{\rm pc}^{-2}$ and $\Sigma_{\rm SFR}=0.54_{-0.27}^{+0.89}\,{\rm M}_{\odot}\,{\rm yr}^{-1}\,{\rm kpc}^{-2}$ respectively. The intrinsic properties of RCSGA 032727--132609 show an enhanced star-formation activity compared to local spiral galaxies with similar molecular gas densities, supporting the ongoing merger-starburst phase scenario. RCSGA 032727--132609 also appears to be a low--density starburst galaxy similar to local blue compact dwarf galaxies, which have been suggested as local analogs to high-redshift low-metallicity starburst systems. Finally, the CO excitation level in the galaxy is consistent with having the peak at ${\rm J}\sim5$, with a higher excitation concentrated in the star-forming clumps., Comment: 11 pages, 4 figures, accepted for publication in ApJL

We have investigated the structure of Co$_2$CO$_3$(OH)$_2$ using Density Functional Theory (DFT) simulations as well as Powder X-Ray Diffraction (PXRD) measurements on samples synthesized via deposition from aqueous solution. We consider two possible monoclinic phases, with closely related but symmetrically different crystal structures, based on those of the minerals malachite and rosasite, as well as an orthorhombic phase that can be seen as a common parent structure for the two monoclinic phases, and a triclinic phase with the structure of the mineral kolwezite. Our DFT simulations predict that the rosasite-like and the malachite-like phases are two different local minima of the potential energy landscape for Co$_2$CO$_3$(OH)$_2$, and are practically degenerate in energy, while the orthorhombic and triclinic structures are unstable and experience barrierless transformations to the malachite phase upon relaxation. The best fit to the PXRD data is obtained using a rosasite model (monoclinic with space group P1121/n and cell parameters a = 3.1408(4) {\AA}, b = 12.2914(17) {\AA}, c = 9.3311(16) {\AA}, $\gamma$ = 82.299(16)$^{\circ}$). However, some features of the PXRD pattern are still not well accounted for by this refinement and the residual parameters are relatively poor. We discuss the relationship between the rosasite and malachite phases of Co$_2$CO$_3$(OH)$_2$ and show that they can be seen as polytypes. Based on the similar calculated stability of these two polytypes, we speculate that some level of stacking disorder could account for the poor fit of our PXRD data. The possibility that Co$_2$CO$_3$(OH)$_2$ could crystallize, under different growth conditions, as either rosasite or malachite, or even as a stacking-disordered phase intermediate between the two, requires further investigation., Comment: To appear in Acta Crystallographica Section B

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$)., Comment: Accepted for publication in the Astrophysical Journal

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

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., Comment: Accepted for publication in the Astrophysical Journal

We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%., Comment: 33 pages, 19 figures, 14 tables, accepted for publication in ApJ, updated to match version in press, paper VI of the ASPECS (ALMA Spectroscopic Survey in the Hubble Ultra Deep Field) series

We present a survey of $z\sim3$ Ly$\alpha$ emitters (LAEs) within the fields of the VLT LBG Redshift Survey. The data encompasses 5 independent survey fields co-spatial with spectroscopic LBG data and covering a larger total area than previously analysed for LAE number counts and clustering. This affords an improved analysis over previous work by minimising the effects of cosmic variance and allowing the cross-clustering analysis of LAEs and LBGs. Our photometric sample consists of $\approx600$ LAE candidates, over an area of 1.07~deg$^2$, with equivalent widths of $\gtrsim65$~\AA\ and a flux limit of $\approx2\times10^{-17}$~erg~cm$^{-2}$~s$^{-1}$. From spectroscopic follow-up, we measured a success rate of $78\pm18\%$. We find the $R$-band continuum luminosity function to be $\sim10\times$ lower than the luminosity function of LBGs at this redshift, consistent with previous studies. Exploiting the large area of the survey, we estimate the LAE auto-correlation function and find a clustering length of $r_0=2.86\pm0.33~h^{-1}$~Mpc, low compared to the $z\sim3$ LBG population, but somewhat higher than previous LAE measurements. This corresponds to a median halo mass of $M_{\rm DM}=10^{11.0\pm0.3}~h^{-1}~$M$_{\odot}$. We present an analysis of clustering length versus continuum magnitude and find that the measurements for LAEs and LBGs are consistent at faint magnitudes. Our combined dataset of LAEs and LBGs allows us to measure, for the first time, the LBG-LAE cross-correlation, finding a clustering length of $r_0=3.29\pm0.57~h^{-1}$~Mpc and a LAE halo mass of $10^{11.1\pm0.4}~h^{-1}$~M$_{\odot}$. Overall, we conclude that LAEs inhabit primarily low mass halos, but form a relatively small proportion of the galaxy population found in such halos., Comment: 20 pages, 15 figures, published in MNRAS