Search

Context. The Fred Young Submillimeter Telescope (FYST) line intensity mapping (LIM) survey will measure the power spectrum (PS) of the singly ionized carbon 158 $\rm \mu$m fine-structure line, [CII], to trace the appearance of the first galaxies that emerged during and right after the epoch of reionization (EoR, $6

We use the Atacama Large sub/Millimetre Array (ALMA) to efficiently observe spectral lines across Bands 3, 4, 5, 6, 7, and 8 at high-resolution (0.5" - 0.1") for 16 bright southern Herschel sources at $1.5 < z < 4.2$. With only six and a half hours of observations, we reveal 66 spectral lines in 17 galaxies. These observations detect emission from CO (3-2) to CO(18-17), as well as atomic ([CI](1-0), (2-1), [OI] 145 $\mu$m and [NII] 205 $\mu$m) lines. Additional molecular lines are seen in emission (${\rm H_2O}$ and ${\rm H_2O^+}$) and absorption (OH$^+$ and CH$^+$). The morphologies based on dust continuum ranges from extended sources to strong lensed galaxies with magnifications between 2 and 30. CO line transitions indicate a diverse set of excitation conditions with a fraction of the sources ($\sim 35$%) showcasing dense, warm gas. The resolved gas to star-formation surface densities vary strongly per source, and suggest that the observed diversity of dusty star-forming galaxies could be a combination of lensed, compact dusty starbursts and extended, potentially-merging galaxies. The predicted gas depletion timescales are consistent with 100 Myr to 1 Gyr, but require efficient fueling from the extended gas reservoirs onto the more central starbursts, in line with the Doppler-shifted absorption lines that indicate inflowing gas for two out of six sources. This pilot paper explores a successful new method of observing spectral lines in large samples of galaxies, supports future studies of larger samples, and finds that the efficiency of this new observational method will be further improved with the planned ALMA Wideband Sensitivity Upgrade., Comment: Accepted for publication in MNRAS. 42 pages, including 20 figures + spectra. Comments and discussion are warmly welcomed

Herschel surveys have found large numbers of sources with red far-IR colours, and spectral energy distributions (SEDs) rising from 250 to 500$\mu$m: 500 risers. The nature and role of these sources is not fully understood. We here present Submillimeter Array (SMA) interferometric imaging at 200 GHz of a complete sample of five 500 risers with F500 $>$ 44 mJy selected within a 4.5 square degree region of the XMMLSS field. These observations can resolve the separate components of multiple sources and allow cross identification at other wavelengths using the extensive optical-to-IR data in this field. Of our five targets we find that two are likely gravitationally lensed, two are multiple sources, and one an isolated single source. Photometric redshifts, using optical-to-IR data and far-IR/submm data, suggest they lie at redshifts $z \sim 2.5 - 3.5$. Star formation rates and stellar masses estimated from the SEDs show that the majority of our sources lie on the star-formation rate-stellar mass `main sequence', though with outliers both above and below this relation. Of particular interest is our most multiple source which consists of three submm emitters and one submm-undetected optical companion within a 7 arcsecond region, all with photometric redshifts $\sim$ 3. One of the submm emitters in this group lies above the `main sequence', while the optical companion lies well below the relation, and has an estimated stellar mass of 3.3$\pm 1.3 \times$10$^{11}$ M$_{\odot}$. We suggest this object is a forming brightest cluster galaxy (BCG) in the process of accreting actively star forming companions., Comment: Accepted for publication in MNRAS

We present new JWST/NIRSpec IFU observations of the J1000+0234 system at $z=4.54$, the dense core of a galaxy protocluster hosting a massive, dusty star forming galaxy (DSFG) with a low luminosity radio counterpart. The new data reveals two extended, high equivalent width (EW$_0 > 1000$ {\AA}) nebulae at each side of the DSFG disk along its minor axis (namely O3-N and O3-S). On one hand, O3-N's spectrum shows a prominent FWHM $\sim1300$ km s$^{-1}$ broad and blueshifted component, suggesting an outflow origin. On the other hand, O3-S stretches over parsec and has a velocity gradient that spans $800$ km s$^{-1}$ but no evidence of a broad component. Both sources, however, seem to be powered at least partially by an active galactic nucleus (AGN), so we classify them as extended emission-line regions (EELRs). The strongest evidence comes from the detection of the high-ionization [Ne V] $\lambda3427$ line toward O3-N, which paired with the non-detection of hard X-rays implies an obscuring column density above the Compton-thick regime. In O3-S, the [Ne V] line is not detected, but we measure a He II well above the expectation for star formation. We interpret this as O3-S being externally irradiated by the AGN, akin to the famous Hanny's Voorwerp object in the local Universe. In addition, more classical line ratio diagnostics (e.g. [O III]/H$\beta$ vs [N II]/H$\alpha$) put the DSFG itself in the AGN region of the diagrams, and hence the most probable host of the AGN. These results showcase the ability of JWST of unveiling highly obscured AGN at high redshifts., Comment: 5 pages, 4 figures plus 5 appendices (incl. 3 extra figures and one table). Submitted to A&A on July 17th 2024

We report on deep SCUBA-2 observations at 850$\mu$m and NOEMA spectroscopic measurements at 2 mm of the environment surrounding the luminous, massive ($M_{*} \approx 2 \times 10^{11}$ M$_{\odot}$) Herschel-selected source HerBS-70. This source was revealed by previous NOEMA observations to be a binary system of dusty star-forming galaxies at $z= 2.3$, with the East component (HerBS-70E) hosting an Active Galactic Nucleus (AGN). The SCUBA-2 observations detected, in addition to the binary system, twenty-one sources at $> 3.5 \sigma$ over an area of $\sim 25$ square comoving Mpc with a sensitivity of $\sigma_{850} = 0.75$ mJy. The surface density of continuum sources around HerBS-70 is three times higher than for field galaxies. The NOEMA spectroscopic measurements confirm the protocluster membership of three of the nine brightest sources through their CO(4 - 3) line emission, yielding a volume density 36 times higher than for field galaxies. All five confirmed sub-mm galaxies in the HerBS-70 system have relatively short gas depletion times ($80 - 500$ Myr), indicating the onset of quenching for this protocluster core due to the depletion of gas. The dark matter halo mass of the HerBS-70 system is estimated around $5 \times{} 10^{13}$ M$_{\odot}$, with a projected current-day mass of $10^{15}$ M$_{\odot}$, similar to the local Virgo and Coma clusters. These observations support the claim that DSFGs, in particular the ones with observed multiplicity, can trace cosmic overdensities., Comment: 19 pages, 13 figures, accepted for publication in MNRAS

We have implemented a novel method to create simulated [CII] emission line intensity mapping (LIM) data cubes using COSMOS 2020 galaxy catalogue data. It allows us to provide solid lower limits for previous simulation-based model predictions and the expected signal strength of upcoming surveys. We applied [CII]158$\mu$m luminosity models to COSMOS 2020 to create LIM cubes covering a 1.2$\times$1.2 deg$^2$ sky area. These models were derived using galaxy bulk property data from the ALPINE-ALMA survey over the redshift range $4.41 when using FYST parameters. However, our high-redshift results remain inconclusive because of the poor completeness of COSMOS 2020 at $z$>6.3. These predictions will be improved on the basis of future JWST data., Comment: 28 pages, 22 figures, accepted by A&A

The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions., Comment: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

We present results of dust continuum and [CII]$\,158\,{\rm \mu m}$ emission line observations of a remarkably UV-luminous ($M_{\rm UV}=-21.6$) galaxy at $z=10.603$: GN-z11. Using the Northern Extended Millimeter Array (NOEMA), observations have been carried out over multiple observing cycles. We achieved a high sensitivity resulting in a $\lambda_{\rm rest}=160\,{\rm \mu m}$ continuum $1\,\sigma$ depth of $13.0\,\rm{\mu Jy/beam}$ and a [CII] emission line $1\,\sigma$ sensitivity of $31\,\rm{mJy/beam\,km/s}$ using $50\,\rm{km/s}$ binning with a $\sim 2\,{\rm arcsec}$ synthesized beam. Neither dust continuum nor [CII]$\,158\,{\rm \mu m}$ line emission are detected at the expected frequency of $\nu_{\rm [CII]} = 163.791\,\rm{GHz}$ and the sky location of GN-z11. The upper limits show that GN-z11 is neither luminous in $L_{\rm IR}$ nor $L_{\rm [CII]}$, with a dust mass $3\,\sigma$ limit of ${\rm log}(M_{\rm dust}/{\rm M_{\odot}}) < 6.5-6.9$ and with a [CII] based molecular gas mass $3\,\sigma$ limit of ${\rm log}(M_{\rm mol,[CII]}/{\rm M_{\odot}}) < 9.3$. Together with radiative transfer calculations, we also investigated the possible cause of the dust poor nature of the GN-z11 showed by the blue color in the UV continuum of GN-z11 ($\beta_{\rm UV}=-2.4$), and found that $\gtrsim3\times$ deeper observations are crucial to study dust production at very high-redshift. Nevertheless, our observations show the crucial role of deep mm/submm observations of very high redshift galaxies to constrain multiple phases in the interstellar medium., Comment: submitted to MNRAS, 7 pages, 4 figures, 1 table

The z-GAL survey observed 137 bright Herschel-selected targets with the IRAM NOrthern Extended Millimeter Array, with the aim to measure their redshift and study their properties. Several of them have been resolved into multiple sources. Consequently, robust spectroscopic redshifts have been measured for 165 individual galaxies in the range 0.8

(Abridged) We present the dust properties of 125 bright Herschel galaxies selected from the z-GAL survey. The large instantaneous bandwidth of NOEMA provides an exquisite sampling of the underlying dust continuum emission at 2 and 3 mm in the observed frame, with flux densities in at least four side bands for each source. Together with the available Herschel 250, 350, and 500 micron and SCUBA-2 850 micron flux densities, the spectral energy distribution of each source can be analyzed from the far-infrared to the millimeter, with a fine sampling of the Rayleigh-Jeans tail. This wealth of data provides a solid basis to derive robust dust properties, in particular the dust emissivity index, beta, and the dust temperature, T(dust). In order to demonstrate our ability to constrain the dust properties, we used a flux-generated mock catalog and analyzed the results under the assumption of an optically thin and optically thick modified black body emission. For the z-GAL sources, we report a range of dust emissivities with beta ~ 1.5 - 3 estimated up to high precision with relative uncertainties that vary in the range 7% - 15%, and an average of 2.2 +/- 0.3. We find dust temperatures varying from 20 to 50 K with an average of T(dust) ~ 30 K for the optically thin case and ~38 K in the optically thick case. For all the sources, we estimate the dust masses and apparent infrared luminosities (based on the optically thin approach). An inverse correlation is found between T(dust) and beta, which is similar to what is seen in the local Universe. Finally, we report an increasing trend in the dust temperature as a function of redshift at a rate of 6.5 +/- 0.5 K/z for this 500 micron-selected sample. Based on this study, future prospects are outlined to further explore the evolution of dust temperature across cosmic time., Comment: Accepted for publication on A&A; 32 pages; 26 figures

(Abridged) Using the IRAM NOEMA interferometer, we measures the redshifts of 126 bright galaxies detected in the Herschel H-ATLAS, HeLMS, and HerS surveys. We report reliable spectroscopic redshifts for a total of 124 of the Herschel-selected galaxies. The redshifts are estimated from scans of the 3 and 2-mm bands (and, in one case, the 1-mm band) and are based on the detection of at least two emission lines. Together with the Pilot Programme (Neri et al. 2020), including spectroscopic redshifts of 11 sources, our survey has derived precise redshifts for 135 bright Herschel-selected galaxies, making it the largest sample of high-z galaxies with robust redshifts to date. Most emission lines detected are from 12CO (mainly from J=2-1 to 5-4), with some sources seen in [CI] and H2O emission lines. The spectroscopic redshifts are in the range 0.8

We present Atacama Large Millimeter/Submillimeter Array (ALMA) [CII] and $\sim158$ $\rm\mu m$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; $L_{\rm IR}=1.5^{+0.8}_{-0.5}\times10^{12}$ L$_\odot$) at $z=7.31$, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a significant stellar mass of $M_{*}=8^{+4}_{-2}\times10^{9}$ M$_\odot$. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR $=199^{+101}_{-63}$ M$_\odot$ yr$^{-1}$. This is about four times the SFR estimated from an extrapolated main-sequence. We also infer a [CII]-based molecular gas mass of $M_{\rm H_2}=5.1^{+5.1}_{-2.6}\times10^{10}$ $M_\odot$, implying a molecular gas depletion time of $ t_{\rm depl, H_2}=0.3^{+0.3}_{-0.2}$ Gyr. We observe a [CII] velocity gradient consistent with disc rotation, but given the current resolution we cannot rule out a more complex velocity structure such as a merger. The spectrum exhibits excess [CII] emission at large positive velocities ($\sim500$ km s$^{-1}$), which we interpret as either a merging companion or an outflow. In the outflow scenario, we derive a lower limit of the mass outflow rate of 200 M$_\odot$ yr$^{-1}$, which is consistent with expectations for a star formation-driven outflow. Given its large stellar mass, SFR and molecular gas reservoir $\sim700$ Myr after the Big Bang, we explore the future evolution of REBELS-25. Considering a simple, conservative model assuming an exponentially declining star formation history, constant star formation efficiency, and no additional gas inflow, we find that REBELS-25 has the potential to evolve into a galaxy consistent with the properties of high-mass quiescent galaxies recently observed at $z\sim4$., Comment: Accepted for publication in MNRAS. 21 pages, 8 figures

We analyse the molecular and atomic emission lines of 71 bright Herschel-selected galaxies between redshifts 1.4 to 4.6 detected by the Atacama Large Millimetre/submillimetre Array. These lines include a total of 156 CO, [C I], and H2O emission lines. For 46 galaxies, we detect two transitions of CO lines, and for these galaxies we find gas properties similar to those of other dusty star-forming galaxy (DSFG) samples. A comparison to photo-dissociation models suggests that most of Herschel-selected galaxies have similar interstellar medium conditions as local infrared-luminous galaxies and high-redshift DSFGs, although with denser gas and more intense far-ultraviolet radiation fields than normal star-forming galaxies. The line luminosities agree with the luminosity scaling relations across five orders of magnitude, although the star-formation and gas surface density distributions (i.e., Schmidt-Kennicutt relation) suggest a different star-formation phase in our galaxies (and other DSFGs) compared to local and low-redshift gas-rich, normal star-forming systems. The gas-to-dust ratios of these galaxies are similar to Milky Way values, with no apparent redshift evolution. Four of 46 sources appear to have CO line ratios in excess of the expected maximum (thermalized) profile, suggesting a rare phase in the evolution of DSFGs. Finally, we create a deep stacked spectrum over a wide rest-frame frequency (220-890 GHz) that reveals faint transitions from HCN and CH, in line with previous stacking experiments., Comment: 30 pages, 17 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal. Comments are warmly welcomed

We present the results of a survey of CO(1-0) emission in 14 infrared luminous dusty star forming galaxies (DSFGs) at 2 < z < 4 with the NSF's Karl G. Jansky Very Large Array. All sources are detected in CO(1-0), with an ~1arcsec angular resolution. Seven sources show extended and complex structure. We measure CO luminosities of $({\mu})L'_{CO(1-0)}=0.4-2.9x10^{11}$ K km s$^{-1}$ pc$^2$, and molecular gas masses of (${\mu}$)M$_{H2}$ = 1.3 - 8.6 x 10$^{11}$ Mo, where ({\mu}) is the magnification factor. The derived molecular gas depletion times of t$_{\rm dep}$ = 40 - 460 Myr, cover the expected range of both normal star forming galaxies and starbursts. Comparing to the higher-J CO transitions previously observed for the same sources, we find CO temperature brightness ratios of r$_{32/10}$ = 0.4 - 1.4, r$_{43/10}$ = 0.4 - 1.7, and r$_{54/10}$ = 0.3 - 1.3. We find a wide range of CO spectral line energy distributions (SLEDs), in agreement with other high-z DSFGs, with the exception of three sources that are most comparable to the Cloverleaf and APM08279+5255. Based on radiative transfer modelling of the CO SLEDs we determine densities of n$_{H2}$ = 0.3 - 8.5 x 10$^3$ cm$^{-3}$ and temperatures of T$_K$ = 100 - 200 K. Lastly, four sources are detected in the continuum, three have radio emission consistent with their infrared derived star formation rates, while HerBS-70E requires an additional synchrotron radiation component from an active galactic nucleus. Overall, we find that even though the sample is similarly luminous in the infrared, by tracing the CO(1-0) emission a diversity of galaxy and excitation properties are revealed, demonstrating the importance of CO(1-0) observations in combination to higher-J transitions., Comment: 25pages , 13 figures, Accepted for publication at ApJ

We present 101 and 151 GHz ALMA continuum images for 85 fields selected from Herschel observations that have 500 micron flux densities >80 mJy and 250-500 micron colours consistent with z > 2, most of which are expected to be gravitationally lensed or hyperluminous infrared galaxies. Approximately half of the Herschel 500 micron sources were resolved into multiple ALMA sources, but 11 of the 15 brightest 500 micron Herschel sources correspond to individual ALMA sources. For the 37 fields containing either a single source with a spectroscopic redshift or two sources with the same spectroscopic redshift, we examined the colour temperatures and dust emissivity indices. The colour temperatures only vary weakly with redshift and are statistically consistent with no redshift-dependent temperature variations, which generally corresponds to results from other samples selected in far-infrared, submillimetre, or millimetre bands but not to results from samples selected in optical or near-infrared bands. The dust emissivity indices, with very few exceptions, are largely consistent with a value of 2. We also compared spectroscopic redshifts to photometric redshifts based on spectral energy distribution templates designed for infrared-bright high-redshift galaxies. While the templates systematically underestimate the redshifts by ~15%, the inclusion of ALMA data decreases the scatter in the predicted redshifts by a factor of ~2, illustrating the potential usefulness of these millimetre data for estimating photometric redshifts., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Water vapor (H$_{2}$O) is one of the brightest molecular emitters after carbon monoxide (CO) in galaxies with high infrared (IR) luminosity, and allows us to investigate the warm dense phase of the interstellar medium (ISM) where star formation occurs. However, due to the complexity of its radiative spectrum, H$_{2}$O is not frequently exploited as an ISM tracer in distant galaxies. Therefore, H$_{2}$O studies of the warm and dense gas at high-$z$ remains largely unexplored. In this work we present observations conducted with the Northern Extended Millimeter Array (NOEMA) toward three $z>6$ IR-bright quasars J2310+1855, J1148+5251, and J0439+1634 targeted in their multiple para-/ortho-H$_{2}$O transitions ($3_{12}-3_{03}$, $1_{11}-0_{00}$, $2_{20}-2_{11}$, and $4_{22}-4_{13}$), as well as their far-IR (FIR) dust continuum. By combining our data with previous measurements from the literature we estimate dust masses and temperatures, continuum optical depths, IR luminosities, and the star-formation rates from the FIR continuum. We model the H$_{2}$O lines using the MOLPOP-CEP radiative transfer code and find that water vapor lines in our quasar host galaxies are primarily excited in warm dense (gas kinetic temperature and density of $T_{\rm kin} = 50\,{\rm K}$, $n_{\rm H_{2}}\sim 10^{4.5}-10^{5}\,{\rm cm^{-3}}$) molecular medium with water vapor column density of $N_{\rm H_{2}O}\sim 2\times10^{17}-3\times10^{18}\,{\rm cm^{-3}}$. High-$J$ H$_{2}$O lines are mainly radiatively pumped by the intense optically-thin far-IR radiation field associated with a warm dust component with temperatures of $T_{\rm dust}\sim 80-190\,{\rm K}$ that account for $<5-10\%$ of the total dust mass. Our results are in agreement with expectations based on the H$_{2}$O spectral line energy distribution of local and high-$z$ ultra-luminous IR galaxies and AGN. [abridged], Comment: Accepted for publication in A&A: 18 pages, 9 figures

Dusty star-forming galaxies (DSFGs) at redshift z$\geq$1 are among the most vigorously star-forming galaxies in the Universe. However, their dense ($\geq$10$^5$ cm$^{-3}$ ) gas phase - typically traced by HCN(1-0) - remains almost entirely unexplored: only two DSFGs have been detected in HCN(1-0) to date. We present results of a JVLA survey of the J=1-0 transition of HCN, HCO+, and HNC(1-0) in six strongly lensed DSFGs at z = 2.5 - 3.3, effectively doubling the number of DSFGs with deep observations of these lines. We detect HCN(1-0) emission in one source (J1202+5354, 4.4$\sigma$), with a tentative HCO+ (1-0) detection in another (J1609+6045, 3.3$\sigma$). Spectral stacking yields strict upper limits on the HCN/FIR ($\leq$3.6$\times$10$^{-4}$) and HCN/CO(1-0) ratios ($\leq$0.045). The inferred HCN/FIR ratios (a proxy for the star-formation efficiency) are consistent with those in z$\sim$0 FIR-luminous starbursts. However, the HCN/CO ratios - a proxy for the dense-gas fraction - are a factor of a few lower than suggested by the two previous DSFG detections. Our results imply that most DSFGs have low dense-gas fractions. A comparison with Krumholz & Thompson (2007) models of star-forming galaxies indicates that the bulk of gas in DSFGs is at lower densities ($\approx$10$^2$ cm$^{-3}$ ), similar to "normal" star-forming galaxies, rather than ultraluminous starbursts., Comment: Submitted to A&A 28th April 2022, accepted 19th July 2022

We make use of ultra-deep 3 GHz Karl G. Jansky Very Large Array observations of the COSMOS field from the multi-band COSMOS-XS survey to infer radio luminosity functions (LFs) of star-forming galaxies (SFGs). Using $\sim$1300 SFGs with redshifts out to $z\sim4.6$, and fixing the faint and bright end shape of the radio LF to the local values, we find a strong redshift trend that can be fitted by pure luminosity evolution with the luminosity parameter given by $\alpha_L \propto (3.40 \pm 0.11) - (0.48 \pm 0.06)z$. We then combine the ultra-deep COSMOS-XS data-set with the shallower VLA-COSMOS $\mathrm{3\,GHz}$ large project data-set over the wider COSMOS field in order to fit for joint density+luminosity evolution, finding evidence for significant density evolution. By comparing the radio LFs to the observed far-infrared (FIR) and ultraviolet (UV) LFs, we find evidence of a significant underestimation of the UV LF by $21.6\%\, \pm \, 14.3 \, \%$ at high redshift ($3.3\,<\,z\,<\,4.6$, integrated down to $0.03\,L^{\star}_{z=3}$). We derive the cosmic star formation rate density (SFRD) by integrating the fitted radio LFs and find that the SFRD rises up to $z\,\sim\,1.8$ and then declines more rapidly than previous radio-based estimates. A direct comparison between the radio SFRD and a recent UV-based SFRD, where we integrate both LFs down to a consistent limit ($0.038\,L^{\star}_{z=3}$), reveals that the discrepancy between the radio and UV LFs translates to a significant ($\sim$1 dex) discrepancy in the derived SFRD at $z>3$, even assuming the latest dust corrections and without accounting for optically dark sources., Comment: Submitted to ApJ; 27 pages, 13 figures, 4 tables

The luminosity functions at z < 4 - 5 suggest that most galaxies have a relatively low stellar mass (logM_star = 10) and a low dust attenuation (A_FUV = 1.0). The physical properties of these objects are quite homogeneous. We used an approach where we combined their rest-frame far-infrared and submillimeter emissions and utilized the universe and the redshift as a spectrograph to increase the amount of information in a collective way. From a subsample of 27 ALMA-detected galaxies at z > 4.5, we built an infrared spectral energy distribution composite template. It was used to fit, with CIGALE, the 105 galaxies (detections and upper limits) in the sample from the FUV to the FIR. The derived physical parameters provide information to decipher the nature of the dust cycle and of the stellar populations in these galaxies. The derived IR composite template is consistent with the galaxies in the studied sample. A delayed star formation history with tau_main = 500 Myrs is slightly favored by the statistical analysis as compared to a delayed with a final burst or a continuous star formation history. The position of the sample in the star formation rate (SFR)- M_star diagram is consistent with previous papers. The redshift evolution of the log M_star versus A_FUV relation is in agreement with evolution in the redshift of this relation. This evolution is necessary to explain the cosmic evolution of the average dust attenuation of galaxies. Evolution is also observed in the L_dust/ L_FUV (IRX) versus UV slope beta_FUV diagram: younger galaxies have bluer beta_FUV. We modeled the shift of galaxies in the IRX versus the beta_FUV diagram with the mass-weighted age as a free parameter, and we provide an equation to make predictions., Comment: 40 pages, 13 figures

We present SCUBA-2 and SMA follow-up observations of four candidate high redshift Dusty Star-Forming Galaxies, selected as sources with rising SEDs in the 250, 350 and 500$\mu$m Herschel SPIRE bands. Previous SMA observations showed no counterparts to these sources, but in our deeper sub-mm observations we detect counterparts to all four 500$\mu$m risers, with three resolving into multiple systems. For these three multiple systems, the SMA 345GHz ($\approx 870\mu$m) observations recover $123 \pm 73\%$, $60 \pm 15\%$ and $19 \pm 4\%$ respectively of the integrated 850$\mu$m flux density from SCUBA-2, indicating that there may be additional sources below our SMA detection limit making up a dense, protocluster core. The fourth 500$\mu$m riser was observed at a lower frequency and so we cannot make a similar comparison. We estimate photometric redshifts based on FIR/sub-mm colours, finding that 3/4 likely lie at $z \geq 2$. This fits with the interpretation that the 500$\mu$m riser selection criterion selects both intrinsically red, individual galaxies at $z > 4$, and multiple systems at more moderate redshifts, artificially reddened by the effects of blending. We use the SCUBA-2 850$\mu$m maps to investigate the environments of these 500$\mu$m risers. By constructing cumulative number counts and estimating photometric redshifts for surrounding SCUBA-2 detections, we find that one of our 500$\mu$m risers could plausibly reside in a $z \geq 2$ protocluster. We infer that bright 500$\mu$m risers with faint 850$\mu$m flux densities are typically multiple systems at $z \geq 2$ that may reside in overdensities of bright sub-mm galaxies., Comment: 20 pages, 13 figures, 3 tables. Submitted to MNRAS

We include a fully coupled treatment of metal and dust enrichment into the Delphi semi-analytic model of galaxy formation to explain the dust content of 13 Lyman Break Galaxies (LBGs) detected by the Atacama Large millimetre Array (ALMA) REBELS Large Program at $z\simeq 7$. We find that the galaxy dust mass, $M_d$, is regulated by the combination of SNII dust production, astration, shock destruction, and ejection in outflows; grain growth (with a standard timescale $\tau_0= 30$ Myr) plays a negligible role. The model predicts a dust-to-stellar mass ratio of $\sim 0.07-0.1\%$ and a UV-to-total star formation rate relation such that $log (\psi_{\rm UV}) = -0.05 ~[log (\psi)]^{2} + 0.86 ~log(\psi) -0.05$ (implying that 55-80\% of the star formation is obscured) for REBELS galaxies with stellar mass $M_* = 10^{9-10} M_\odot$. This relation reconciles the intrinsic UV luminosity of LBGs with their observed luminosity function at $z=7$. However, 2 out of the 13 systems show dust-to-stellar mass ratios ($\sim 0.94-1.1\%$) that are up to $18\times$ larger than expected from the fiducial relation. Due to the physical coupling between dust and metal enrichment, even decreasing $\tau_0$ to very low values (0.3 Myr) only increases the dust-to-stellar mass ratio by a factor $ \sim 2$. Given that grain growth is not a viable explanation for such high observed ratios of the dust-to-stellar mass, we propose alternative solutions., Comment: Accepted to MNRAS

We analyse FIR dust continuum measurements for 14 galaxies ($z\approx 7$) in the ALMA REBELS LP to derive their physical properties. Our model uses three input data: (a) the UV spectral slope, $\beta$, (b) the observed UV continuum flux at $1500$A, $F_{\rm UV}$, (c) the observed continuum flux at $\approx 158\mu$m, $F_{158}$, and considers Milky Way (MW) and SMC extinction curves, along with different dust geometries. We find that REBELS galaxies have (28-90.5)% of their star formation obscured; the total (UV+IR) star formation rates are in the range $31.5 < {\rm SFR}/ (M_\odot {\rm yr}^{-1}) < 129.5$. The sample-averaged dust mass and temperature are $(1.3\pm 1.1)\times 10^7 M_\odot$ and $52 \pm 11$ K, respectively. In some galaxies dust is abundant (REBELS-14, $M'_d \approx 3.4 \times 10^7 M_\odot$), or hot (REBELS-18, $T'_d \approx 67$ K). The dust distribution is compact ($<0.3$ kpc for 70% of the galaxies). The dust yield per supernova is $0.1 \le y_d/M_\odot \le 3.3$, with 70% of the galaxies requiring $y_d < 0.25 M_\odot$. Three galaxies (REBELS-12, 14, 39) require $y_d > 1 M_\odot$. With the SFR predicted by the model and a MW extinction curve, REBELS galaxies detected in [CII] nicely follow the local $L_{\rm CII}-$SFR relation, and are approximately located on the Kennicutt-Schmidt relation. The sample-averaged gas depletion time is of $0.11\, y_P^{-2}$ Gyr, where $y_P$ is the ratio of the gas-to-stellar distribution radius. For some systems a solution simultaneously matching the observed ($\beta, F_{\rm UV}, F_{158}$) values cannot be found. This occurs when the index $I_m = (F_{158}/F_{\rm UV})/(\beta-\beta_{\rm int})$, where $\beta_{\rm int}$ is the intrinsic UV slope, exceeds $I_m^*\approx 1120$ for a MW curve. For these objects we argue that the FIR and UV emitting regions are not co-spatial, questioning the use of the IRX-$\beta$ relation., Comment: 22 pages, 7 figures, accepted by MNRAS. Comments welcome

We study the members of the protocluster around AzTEC3 submillimeter galaxy at z=5.3. We analyzed the data from the MUSE instrument in an area of 1.4x1.4 arcmin^2 around AzTEC3 and derived information on the Lya line in emission. We compared the Lya profile of various regions of the environment with the zELDA radiative transfer model, revealing the neutral gas distribution and kinematics. We identified 10 Lya emitting sources, including 2 regions with extended emission: one embedding AzTEC3 and LBG3, a star-forming galaxy located 12 kpc north of the SMG and another toward LBG-1, a star-forming galaxy located 90 kpc to the southeast. The sources appear distributed in an elongated configuration of about 70'' in extent. The number of sources confirms the overdensity around AzTEC3. For the AzTEC3+LBG3 system, the Lya emission appears redshifted and more spatially extended than the [CII] line emission. Similarly, the Lya line spectrum is broader in velocity than [CII] for LBG1. In the former spectrum, the Lya emission is elongated to the north of LBG3 and to the south of AzTEC3, where a faint Lya emitting galaxy is also located. The elongated structures could resemble tidal features due to the interaction of the two galaxies with AzTEC3. Also, we find a bridge of gas, revealed by the Lya emission between AzTEC3 and LBG3. The Lya emission toward LBG1 embeds its three components. The HI kinematics support the idea of a merger of the three components. Given the availability of CO and [CII] observations from previous campaigns, and our Lya information, we find evidence of starburst-driven phenomena and interactions around AzTEC-3. The stellar mass of the galaxies of the overdensity and the Lya luminosity of the HI nebula associated with AzTEC-3 imply a dark matter halo of 10^12 Msun at z=5.3 that could evolve into a cluster of 2x10^14 Msun at z=0., Comment: Submitted to A&A, Accepted on Feb 4th, 2022

ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [CII] 158$\mu$m line and underlying dust continuum measurements, we derive $T_ d$ in the continuum and [CII] detected $z\approx 7$ galaxies in the ALMA Large Project REBELS sample. We find $39\ \mathrm{K} < T_d < 58\ \mathrm{K}$, and dust masses in the narrow range $M_d = (0.9-3.6)\times 10^7 M_{\odot}$. These results allow us to extend for the first time the reported $T_d(z)$ relation into the Epoch of Reionization. We produce a new physical model that explains the increasing $T_ d(z)$ trend with the decrease of gas depletion time, $t_{dep}=M_g/\mathrm{SFR}$, induced by the higher cosmological accretion rate at early times; this hypothesis yields $T_d \propto (1+z)^{0.4}$. The model also explains the observed $T_d$ scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, $T_d$ only depends on the gas column density (metallicity), $T_d \propto N_H^{1/6}$ ($T_d \propto Z^{-1/6}$). REBELS galaxies are on average relatively transparent, with effective gas column densities around $N_H \simeq (0.03-1)\times 10^{21} \mathrm{cm}^{-2}$. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated $T_d \gg 60$ K, are significantly obscured, low-metallicity systems. In fact $T_d$ is higher in metal-poor systems due to their smaller dust content, which for fixed $L_{ IR}$ results in warmer temperatures., Comment: Accepted for publication in MNRAS

We present spectroscopic measurements for 71 galaxies associated with 62 of the brightest high-redshift submillimeter sources from the Southern fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), while targeting 85 sources which resolved into 142. We have obtained robust redshift measurements for all sources using the 12-m Array and an efficient tuning of ALMA to optimise its use as a redshift hunter, with 73 per cent of the sources having a robust redshift identification. Nine of these redshift identifications also rely on observations from the Atacama Compact Array. The spectroscopic redshifts span a range $1.41

We have carried out the first spatially-resolved investigation of the multi-phase interstellar medium (ISM) at high redshift, using the z=4.24 strongly-lensed sub-millimetre galaxy H-ATLASJ142413.9+022303 (ID141). We present high-resolution (down to ~350 pc) ALMA observations in dust continuum emission and in the CO(7-6), H_2O (2_{1,1} - 2_{0,2}), CI(1-0) and CI(2-1) lines, the latter two allowing us to spatially resolve the cool phase of the ISM for the first time. Our modelling of the kinematics reveals that the system appears to be dominated by a rotationally-supported gas disk with evidence of a nearby perturber. We find that the CI(1-0) line has a very different distribution to the other lines, showing the existence of a reservoir of cool gas that might have been missed in studies of other galaxies. We have estimated the mass of the ISM using four different tracers, always obtaining an estimate in the range (3.2-3.8) x 10^{11} M_sol, significantly higher than our dynamical mass estimate of (0.8-1.3) x 10^{11} M_sol. We suggest that this conflict and other similar conflicts reported in the literature is because the gas-to-tracer ratios are ~4 times lower than the Galactic values used to calibrate the ISM in high-redshift galaxies. We demonstrate that this could result from a top-heavy initial mass function and strong chemical evolution. Using a variety of quantitative indicators, we show that, extreme though it is at z=4.24, ID141 will likely join the population of quiescent galaxies that appears in the Universe at z~3., Comment: This article has been accepted for publication in MNRAS, published by Oxford University Press on behalf of the Royal Astronomical Society. 24 pages

The Cosmic Evolution Survey (COSMOS) has become a cornerstone of extragalactic astronomy. Since the last public catalog in 2015, a wealth of new imaging and spectroscopic data has been collected in the COSMOS field. This paper describes the collection, processing, and analysis of this new imaging data to produce a new reference photometric redshift catalog. Source detection and multi-wavelength photometry is performed for 1.7 million sources across the $2\,\mathrm{deg}^{2}$ of the COSMOS field, $\sim$966,000 of which are measured with all available broad-band data using both traditional aperture photometric methods and a new profile-fitting photometric extraction tool, The Farmer, which we have developed. A detailed comparison of the two resulting photometric catalogs is presented. Photometric redshifts are computed for all sources in each catalog utilizing two independent photometric redshift codes. Finally, a comparison is made between the performance of the photometric methodologies and of the redshift codes to demonstrate an exceptional degree of self-consistency in the resulting photometric redshifts. The $i<21$ sources have sub-percent photometric redshift accuracy and even the faintest sources at $25

Over the past decades, rest-frame ultraviolet (UV) observations have provided large samples of UV luminous galaxies at redshift (z) greater than 6, during the so-called epoch of reionization. While a few of these UV identified galaxies revealed significant dust reservoirs, very heavily dust-obscured sources at these early times have remained elusive. They are limited to a rare population of extreme starburst galaxies, and companions of rare quasars. These studies conclude that the contribution of dust-obscured galaxies to the cosmic star formation rate density at $z>6$ is sub-dominant. Recent ALMA and Spitzer observations have identified a more abundant, less extreme population of obscured galaxies at $z=3-6$. However, this population has not been confirmed in the reionization epoch so far. Here, we report the discovery of two dust-obscured star forming galaxies at $z=6.6813\pm0.0005$ and $z=7.3521\pm0.0005$. These objects are not detected in existing rest-frame UV data, and were only discovered through their far-infrared [CII] lines and dust continuum emission as companions to typical UV-luminous galaxies at the same redshift. The two galaxies exhibit lower infrared luminosities and star-formation rates than extreme starbursts, in line with typical star-forming galaxies at $z\sim7$. This population of heavily dust-obscured galaxies appears to contribute 10-25 per cent to the $z>6$ cosmic star formation rate density., Comment: 16 pages, 7 figures, 4 tables, published in Nature

We present multi-band observations of an extremely dusty star-forming lensed galaxy (HERS1) at $z=2.553$. High-resolution maps of \textit{HST}/WFC3, SMA, and ALMA show a partial Einstein-ring with a radius of $\sim$3$^{\prime\prime}$. The deeper HST observations also show the presence of a lensing arc feature associated with a second lens source, identified to be at the same redshift as the bright arc based on a detection of the [NII] 205$\mu$m emission line with ALMA. A detailed model of the lensing system is constructed using the high-resolution HST/WFC3 image, which allows us to study the source plane properties and connect rest-frame optical emission with properties of the galaxy as seen in sub-millimeter and millimeter wavelengths. Corrected for lensing magnification, the spectral energy distribution fitting results yield an intrinsic star formation rate of about $1000\pm260$ ${\rm M_{\odot}}$yr$^{-1}$, a stellar mass ${\rm M_*}=4.3^{+2.2}_{-1.0}\times10^{11} {\rm M_{\odot}}$, and a dust temperature ${\rm T}_{\rm d}=35^{+2}_{-1}$ K. The intrinsic CO emission line ($J_{\rm up}=3,4,5,6,7,9$) flux densities and CO spectral line energy distribution are derived based on the velocity-dependent magnification factors. We apply a radiative transfer model using the large velocity gradient method with two excitation components to study the gas properties. The low-excitation component has a gas density $n_{\rm H_2}=10^{3.1\pm0.6}$ cm$^{-3}$ and kinetic temperature ${\rm T}_{\rm k}=19^{+7}_{-5}$ K and a high-excitation component has $n_{\rm H_2}=10^{2.8\pm0.3}$ cm$^{-3}$ and ${\rm T}_{\rm k}=550^{+260}_{-220}$ K. Additionally, HERS1 has a gas fraction of about $0.4\pm0.2$ and is expected to last 250 Myr. These properties offer a detailed view of a typical sub-millimeter galaxy during the peak epoch of star-formation activity., Comment: 16 pages, 14 figures

We present a detailed overview of the science goals and predictions for the Prime-Cam direct detection camera/spectrometer being constructed by the CCAT-prime collaboration for dedicated use on the Fred Young Submillimeter Telescope (FYST). The FYST is a wide-field, 6-m aperture submillimeter telescope being built (first light in mid-2024) by an international consortium of institutions led by Cornell University and sited at more than 5600 meters on Cerro Chajnantor in northern Chile. Prime-Cam is one of two instruments planned for FYST and will provide unprecedented spectroscopic and broadband measurement capabilities to address important astrophysical questions ranging from Big Bang cosmology through reionization and the formation of the first galaxies to star formation within our own Milky Way galaxy. Prime-Cam on the FYST will have a mapping speed that is over ten times greater than existing and near-term facilities for high-redshift science and broadband polarimetric imaging at frequencies above 300 GHz. We describe details of the science program enabled by this system and our preliminary survey strategies., Comment: 61 pages, 16 figures. Resubmitted to ApJSS July 11, 2022

The Reionization Era Bright Emission Line Survey (REBELS) is a cycle-7 ALMA Large Program (LP) that is identifying and performing a first characterization of many of the most luminous star-forming galaxies known in the z>6.5 universe. REBELS is providing this probe by systematically scanning 40 of the brightest UV-selected galaxies identified over a 7-deg**2 area for bright 158-micron [CII] and 88-micron [OIII] lines and dust-continuum emission. Selection of the 40 REBELS targets was done by combining our own and other photometric selections, each of which is subject to extensive vetting using three completely independent sets of photometry and template-fitting codes. Building on the observational strategy deployed in two pilot programs, we are increasing the number of massive interstellar medium (ISM) reservoirs known at z>6.5 by ~4-5x to >30. In this manuscript, we motivate the observational strategy deployed in the REBELS program and present initial results. Based on the 60.6 hours of ALMA observations taken in the first year of the program (November 2019 to January 2020), 18 highly significant >~7sigma [CII] lines have already been discovered, the bulk of which (13/18) also show >~3.3 sigma dust-continuum emission. These newly discovered lines more than triple the number of bright ISM-cooling lines known in the z>6.5 universe, such that the number of ALMA-derived redshifts at z>6.5 already rival Lya redshift discoveries. An analysis of the completeness of our search results vs. star formation rate (SFR) suggests an ~79% efficiency in scanning for [CII] when the SFR(UV+IR) is in excess of 28 M_sol/yr. These new LP results further demonstrate ALMA's efficiency as a "redshift machine", particularly in the epoch of reionization., Comment: 30 pages, 4 tables, 18 figures, accepted for publication in ApJ

The dust content of normal galaxies and the dust mass density (DMD) at high-z (z>4) are unconstrained given the source confusion and the sensitivity limitations of previous observations. The ALMA Large Program to INvestigate [CII] at Early Times (ALPINE), which targeted 118 UV-selected star-forming galaxies at 4.4

We present a multi-line survey of the interstellar medium (ISM) in two $z>6$ quasar (QSO) host galaxies, PJ231-20 ($z=6.59$) and PJ308-21 ($z=6.23$), and their two companion galaxies. Observations were carried out using the Atacama Large (sub-)Millimeter Array (ALMA). We targeted eleven transitions including atomic fine structure lines (FSLs) and molecular lines: [NII]$_{\rm 205\mu m}$, [CI]$_{\rm 369\mu m}$, CO ($J_{\rm up} = 7, 10, 15, 16$), H$_2$O $3_{12}-2_{21}$, $3_{21}-3_{12}$, $3_{03}-2_{12}$, and the OH$_{\rm 163\mu m}$ doublet. The underlying far-infrared (FIR) continuum samples the Rayleigh-Jeans tail of the respective dust emission. By combining this information with our earlier ALMA [CII]$_{\rm 158\mu m}$ observations, we explore the effects of star formation and black hole feedback on the galaxies' ISM using the CLOUDY radiative transfer models. We estimate dust masses, spectral indexes, IR luminosities, and star-formation rates from the FIR continuum. The analysis of the FSLs indicates that the [CII]$_{\rm 158\mu m}$ and [CI]$_{\rm 369\mu m}$ emission arises predominantly from the neutral medium in photodissociation regions (PDRs). We find that line deficits are in agreement with those of local luminous infrared galaxies. The CO spectral line energy distributions (SLEDs), reveal significant high-$J$ CO excitation in both quasar hosts. Our CO SLED modeling of the quasar PJ231-20 shows that PDRs dominate the molecular mass and CO luminosities for $J_{\rm up}\le 7$, while the $J_{\rm up}\ge10$ CO emission is likely driven by X-ray dissociation regions produced by the active galactic nucleus (AGN) at the very center of the quasar host [abridged]., Comment: 26 pages, 17 figures, accepted for publication in A&A

The high-frequency radio sky has historically remained largely unexplored due to the typical faintness of sources in this regime, and the modest survey speed compared to observations at lower frequencies. However, high-frequency radio surveys present an invaluable tracer of high-redshift star-formation, as they directly target the faint radio free-free emission. We present deep continuum observations at 34 GHz in the COSMOS and GOODS-North fields from the Karl G. Jansky Very Large Array (VLA), as part of the COLDz survey. The deep COSMOS mosaic spans $\sim10\text{arcmin}^2$ down to $\sigma=1.3\mu\text{Jy beam}^{-1}$, while the wider GOODS-N observations cover $\sim50\text{arcmin}^2$ to $\sigma=5.3\mu\text{Jy beam}^{-1}$. We present the deepest 34 GHz radio number counts to date, with five and thirteen continuum detections in COSMOS and GOODS-N, respectively. Nine galaxies show 34 GHz continuum emission that is originating from star-formation, although for two sources this is likely due to thermal emission from dust. Utilizing deep ancillary radio data at 1.4, 3, 5 and 10 GHz, we decompose the spectra of the remaining seven star-forming galaxies into their synchrotron and thermal free-free components, finding typical thermal fractions and synchrotron spectral indices comparable to those observed in local star-forming galaxies. Using calibrations from the literature, we determine free-free star-formation rates (SFRs), and show that these are in agreement with SFRs from spectral energy distribution fitting and the far-infrared/radio correlation. Our observations place strong direct constraints on the high-frequency radio emission in typical galaxies at high-redshift, and provide some of the first insight in what is set to become a key area of study with future radio facilities as the Square Kilometer Array Phase 1 and next-generation VLA., Comment: Submitted to ApJ; 22 pages + appendices; 13 figures in main text

(Abridged) Exploiting the sensitivity and broad band width of NOEMA, we have studied the molecular gas and dust in the galaxy HerBS-89a, at z=2.95. High angular resolution images reveal a partial 1.0" diameter Einstein ring in the dust continuum emission and the molecular emission lines of 12CO(9-8) and H2O(2_02-1_11). We report the detection of the three fundamental transitions of the molecular ion OH+, seen in absorption; the molecular ion CH+(1-0) seen in absorption (and tentatively in emission); two transitions of amidogen (NH2), seen in emission; and HCN(11-10) and/or NH(1_2-0_1) seen in absorption. The NOEMA data are complemented with VLA data tracing the 12CO(1-0) emission line, which provides a measurement of the total mass of molecular gas and an anchor for a CO excitation analysis. In addition, we present HST imaging that reveals the foreground lensing galaxy in the near-infrared. Together with data from the GTC, we derive a photometric redshift of z(phot)~0.9 for the foreground lensing galaxy. Modelling the lensing of HerBS-89a, we reconstruct the dust continuum and molecular emission lines (magnified by a factor ~4-5) in the source plane. The 12CO(9-8) and H2O emission lines have comparable spatial and kinematic distributions; the source-plane reconstructions do not clearly distinguish between a one-component and a two-component scenario, but the latter accounts for the observed broad line widths. HerBS-89a is a powerful star forming galaxy with a dust-to-gas ratio delta(GDR)~80, a SFR = 614 +/- 59 Msun/yr and a depletion timescale tau(depl) = (3.4 +/- 1.0) 1e8 years. The OH+ and CH+ absorption lines, all have their main velocity component red-shifted by \Delta(V)~100 km/s relative to the global CO reservoir. We argue that these absorption lines trace a rare example of gas inflow towards the center of the galaxy., Comment: 33 pages, Accepted for publication on A&A

We present ultra-deep, matched-resolution Karl G. Jansky Very Large Array (VLA) observations at 10 and $3$ GHz in the COSMOS field: the COSMOS-XS survey. The final 10 and $3$ GHz images cover $\sim16\rm{arcmin}^{2}$ and $\sim180\rm{arcmin}^{2}$ and reach median rms values of $0.41\mu\rm{Jy\,beam}^{-1}$ and $0.53\mu\rm{Jy\,beam}^{-1}$, respectively. Both images have an angular resolution of $\sim 2.0''$. To fully account for the spectral shape and resolution variations across the broad bands, we image all data with a multi-scale, multi-frequency synthesis algorithm. We present source catalogs for the 10 and $3$ GHz image with 91 and 1498 sources, respectively, above a peak brightness threshold of $5\sigma$. We present source counts with completeness corrections included that are computed via Monte Carlo simulations. Our corrected radio counts at $3$ GHz with direct detections down to $\sim2.8\mu$Jy are consistent within the uncertainties with other results at 3 and 1.4 GHz, but extend to fainter flux densities than previous direct detections. The ultra-faint $3$ GHz number counts are found to exceed the counts predicted by the semi-empirical radio sky simulations developed in the framework of the SKA Simulated Skies project, consistent with previous P(D) analyses. Our measured source counts suggest a steeper luminosity function evolution for these faint star-forming sources. The semi-empirical Tiered Radio Extragalactic Continuum Simulation (T-RECS) predicts this steeper evolution and is in better agreement with our results. The $10$ GHz radio number counts also agree with the counts predicted by the T-RECS simulation within the expected variations from cosmic variance. In summary, the multi-band, matched-resolution COSMOS-XS survey in the well-studied COSMOS field provides a high-resolution view of the ultra-faint radio sky that can help guide next generation radio facilities., Comment: Resubmitted to ApJ after incorporating suggestions from the referee; 20 pages, 17 figures. This is a companion paper to Algera et al., on arXiv today

Ultra-deep radio surveys are an invaluable probe of dust-obscured star formation, but require a clear understanding of the relative contribution from radio AGN to be used to their fullest potential. We study the composition of the $\mu$Jy radio population detected in the Karl G. Jansky Very Large Array COSMOS-XS survey based on a sample of 1540 sources detected at 3 GHz over an area of $\sim350\text{arcmin}^2$. This ultra-deep survey consists of a single pointing in the well-studied COSMOS field at both 3 and 10 GHz and reaches RMS-sensitivities of $0.53$ and $0.41\mu$Jy beam$^{-1}$, respectively. We find multi-wavelength counterparts for $97\%$ of radio sources, based on a combination of near-UV/optical to sub-mm data, and through a stacking analysis at optical/near-infrared wavelengths we further show that the sources lacking such counterparts are likely to be high-redshift in nature (typical $z\sim4-5$). Utilizing the multi-wavelength data over COSMOS, we identify AGN through a variety of diagnostics and find these to make up $23.2\pm1.3\%$ of our sample, with the remainder constituting uncontaminated star-forming galaxies. However, more than half of the AGN exhibit radio emission consistent with originating from star-formation, with only $8.8\pm0.8\%$ of radio sources showing a clear excess in radio luminosity. At flux densities of $\sim30\mu$Jy at 3 GHz, the fraction of star-formation powered sources reaches $\sim90\%$, and this fraction is consistent with unity at even lower flux densities. Overall, our findings imply that ultra-deep radio surveys such as COSMOS-XS constitute a highly effective means of obtaining clean samples of star-formation powered radio sources., Comment: Accepted for publication in ApJ; 23 pages + appendices; 15 figures in main text. This is a companion paper to Van der Vlugt et al., on arXiv today

We present Herschel-PACS spectroscopy of four main-sequence star-forming galaxies at z~1.5. We detect [OI]63micron line emission in BzK-21000 at z=1.5213, and measure a line luminosity, L([OI]63micron) = (3.9+/-0.7)x1.E+9 Lsun. Our PDR modelling of the interstellar medium in BzK-21000 suggests a UV radiation field strength, G~320 G0, and gas density, n~1800 cm-3, consistent with previous LVG modelling of the molecular CO line excitation. The other three targets in our sample are individually undetected in these data, and we perform a spectral stacking analysis which yields a detection of their average emission and an [OI]63micron line luminosity, L([OI]63micron) =(1.1+/-0.2)x1E+9 Lsun. We find that the implied luminosity ratio, L([OI]63micron)/L(IR), of the undetected BzK-selected star-forming galaxies broadly agrees with that of low-redshift star-forming galaxies, while BzK-21000 has a similar ratio to that of a dusty star-forming galaxy at z~6. The high [OI]63micron line luminosities observed in BzK-21000 and the $z \sim 1 -3$ dusty and submm luminous star-forming galaxies may be associated with extended reservoirs of low density, cool neutral gas., Comment: accepted for publication in MNRAS

We present MUSE observations of the gas-rich major-merger BR1202-0725 at z~4.7, which constitutes one of the most overdense fields known in the early Universe. We utilise these data in conjunction with existing ALMA observations to compare and contrast the spatially resolved ionised- and cool-gas content of this system which hosts a quasar (QSO), a sub-millimeter galaxy (SMG), the two known optical companions ("LAE1", "LAE2"), and an additional companion discovered in this work "LAE3" just 5 arcsec to the North of the QSO. We find that QSO BR1202-0725 exhibits a large Ly$\alpha$ halo, covering $\approx55$ pkpc on-sky at surface brightness levels of SB$\geq$1E-17 erg/s/cm$^2$/arcsec$^2$. In contrast, the SMG, of similar far-infrared luminosity and star formation rate (SFR), does not exhibit such a Ly$\alpha$ halo. The QSO's halo exhibits high velocity widths ($\sim1000$ km/s) but the gas motion is to some extent kinematically coupled with the previously observed [CII] bridge between the QSO and the SMG. We note that the object known in the literature as LAE2 shows no local peak of Ly$\alpha$ emission, rather, its profile is more consistent with being part of the QSO's extended Ly$\alpha$ halo. The properties of LAE3 are typical of high-redshift LAEs; we measure F$_{\rm{Ly\alpha}}$(LAE3) = $0.24\pm$0.03E-16 erg/s/cm$^2$, corresponding to SFR$_{\rm{Ly\alpha}}\approx\ $5.0$\pm$0.5 M${_\odot}$/yr. The velocity width is $\Delta v$(LAE3) $\approx 400$ km/s, and equivalent width EW$_0$(Ly$\alpha_{\,5\sigma}^{\,lim})\geq 34.05$ $\\A$, consistent with star formation being the primary driver of Ly$\alpha$ emission. We also note a coherent absorption feature at $\sim -400$km/s in spectra from at least three objects; the QSO, LAE1 and "LAE2" which could imply the presence of an expanding neutral gas shell with an extent of at least $24$ pkpc., Comment: 23 pages, accepted for publication in ApJ

Star formation rate (SFR) measurements at z>4 have relied mostly on rest-frame far-ultraviolet (FUV) observations. The corrections for dust attenuation based on IRX-$\beta$ relation are highly uncertain and are still debated in the literature. Hence, rest-frame far-infrared (FIR) observations are necessary to constrain the dust-obscured component of the SFR. In this paper, we exploit the rest-frame FIR continuum observations collected by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) to directly constrain the obscured SFR in galaxies at 4.4

We present SMA observations at resolutions from 0.35 to 3 arcseconds of a sample of 34 candidate high redshift dusty star forming galaxies (DSFGs). These sources were selected from the HerMES Herschel survey catalogues to have SEDs rising from 250 to 350 to 500$\mu$m, a population termed 500-risers. We detect counterparts to 24 of these sources, with four having two counterparts. We conclude that the remaining ten sources that lack detected counterparts are likely to have three or more associated sources which blend together to produce the observed Herschel source. We examine the role of lensing, which is predicted to dominate the brightest (F500 $>$ 60mJy) half of our sample. We find that while lensing plays a role, at least 35% of the bright sources are likely to be multiple sources rather than the result of lensing. At fainter fluxes we find a blending rate comparable to, or greater than, the predicted 40%. We determine far-IR luminosities and star formation rates for the non-multiple sources in our sample and conclude that, in the absence of strong lensing, our 500-risers are very luminous systems with L$_{FIR} > 10^{13}$L$_{\odot}$ and star formation rates $> 1000$M$_{\odot}$/yr., Comment: Accepted for publication in MNRAS

We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA band 7, as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE) in COSMOS and ECDFS. These sources have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z=6. We have looked for counterparts in all the available multi-wavelength and photometric redshift catalogues, and in deeper near- and mid-IR source lists and maps, to identify optically dark sources with no matches in the public catalogues. Our ALMA blind survey allows us to push further the study of the nature and evolution of dusty galaxies at high-z, identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint-end of the infrared LF, showing little evolution from z=2.5 to z=6, and a flat slope up to the highest redshifts. The SFRD obtained by integrating the luminosity function remains almost constant between z=2 and 6, and significantly higher than the optical/UV derivations, showing an important contribution of dusty galaxies and obscured star formation up to high-z. About 14 per cent of the ALPINE serendipitous continuum sources are optically+near-IR dark (six show a counterpart only in the mid-IR and no HST or near-IR identification, while two are detected as [CII] emitters at z=5). The six HST and near-IR dark galaxies with mid-IR counterpart contribute for about 17 per cent of the total SFRD at z=5 and dominate the high-mass end of the stellar mass function at z>3., Comment: 26 pages, 18 figures. Accepted for publication in A&A

We present the [CII]158$\mu$m line luminosity functions (LFs) at $z\sim4-6$ using the ALMA observations of 118 sources, which are selected to have UV luminosity $M_{1500A}<-20.2$ and optical spectroscopic redshifts in COSMOS and ECDF-S. Of the 118 targets, 75 have significant [CII] detections and 43 are upper limits. This is by far the largest sample of [CII] detections which allows us to set constraints to the volume density of [CII] emitters at $z\sim4-6$. But because this is a UV-selected sample, we are missing [CII]-bright but UV-faint sources making our constraints strict lower limits. Our derived LFs are statistically consistent with the $z\sim0$ [CII] LF at $10^{8.25} - 10^{9.75}L_\odot$. We compare our results with the upper limits of the [CII] LF derived from serendipitous sources in the ALPINE maps (Loiacono et al. 2020). We also infer the [CII] LFs based on published far-IR and CO LFs at $z\sim4-6$. Combining our robust lower limits with these additional estimates, we set further constraints to the true number density of [CII] emitters at $z\sim 4 - 6$. These additional LF estimates are largely above our LF at $L_{[CII]}>10^9L_{\odot}$, suggesting that UV-faint but [CII]-bright sources likely make a significant contributions to the [CII] emitter volume density. When we include all the LF estimates, we find that available model predictions underestimate the number densities of [CII] emitters at $z\sim4-6$. Finally, we set a constraint on the molecular gas mass density at $z\sim4-6$, with $\rho_{mol} \sim (2-7)\times10^7M_\odot$\,Mpc$^{-3}$. This is broadly consistent with previous studies., Comment: 17 pages, 9 Figures, Update to match with the published version. Accepted for the publication in ApJ

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

We measure the 850-$\mu$m source densities of 46 candidate protoclusters selected from the Planck High-z catalogue (PHz) and the Planck Catalogue of Compact Sources (PCCS) that were followed up with Herschel-SPIRE and SCUBA-2. This paper aims to search for overdensities of 850-$\mu$m sources in order to select the fields that are most likely to be genuine protoclusters. Of the 46 candidate protoclusters, 25 have significant overdensities ($>$5 times the field counts), 11 have intermediate overdensities (3--5 times the field counts) and 10 have no overdensity ($<$3 times the field counts) of 850-$\mu$m sources. We find that the enhanced number densities are unlikely to be the result of sample variance. Compared with the number counts of another sample selected from Planck's compact source catalogues, this [PHz+PCCS]-selected sample has a higher fraction of candidate protoclusters with significant overdensities, though both samples show overdensities of 850-$\mu$m sources above intermediate level. Based on the estimated star-formation rate densities (SFRDs), we suggest that both samples can efficiently select protoclusters with starbursting galaxies near the redshift at which the global field SFRD peaks ($2 < z < 3$). Based on the confirmation of overdensities found here, future follow-up observations on other PHz targets may greatly increase the number of genuine DSFG-rich clusters/protoclusters., Comment: 8 pages, 1 figure, 1 table, accepted for publication in MNRAS

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

The ALPINE-ALMA large program targets the [CII] 158 $\mu$m line and the far-infrared continuum in 118 spectroscopically confirmed star-forming galaxies between z=4.4 and z=5.9. It represents the first large [CII] statistical sample built in this redshift range. We present details of the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates, we measured the conversion factor from the ALMA 158 $\mu$m rest-frame dust continuum luminosity to the total infrared luminosity (L$_{\rm IR}$) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median L$_{\rm IR}$ of 4.4$\times 10^{11}$ L$_\odot$. We also detected 57 additional continuum sources in our ALMA pointings. They are at lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5$\pm$0.2. We measured the 850 $\mu$m number counts between 0.35 and 3.5 mJy, improving the current interferometric constraints in this flux density range. We found a slope break in the number counts around 3 mJy with a shallower slope below this value. More than 40 % of the cosmic infrared background is emitted by sources brighter than 0.35 mJy. Finally, we detected the [CII] line in 75 of our targets. Their median [CII] luminosity is 4.8$\times$10$^8$ L$_\odot$ and their median full width at half maximum is 252 km/s. After measuring the mean obscured SFR in various [CII] luminosity bins by stacking ALPINE continuum data, we find a good agreement between our data and the local and predicted SFR-L$_{\rm [CII]}$ relations of De Looze et al. (2014) and Lagache et al. (2018)., Comment: 44 pages, 29 figures, 14 tables, accepted by A&A, ALPINE products released at https://cesam.lam.fr/a2c2s/, luminosities in table B2 corrected

The Lya line in the UV and the [CII] line in the FIR are widely used tools to identify galaxies and to obtain insights into ISM properties in the early Universe. By combining data obtained with ALMA in band 7 at ~ 320 GHz as part of the ALMA Large Program to INvestigate [CII] at Early Times (ALPINE) with spectroscopic data from DEIMOS at Keck, VIMOS and FORS2 at the VLT, we assembled a unique sample of 53 main-sequence star-forming galaxies at 4.4 < z < 6 in which we detect both the Lya line and the [CII]. We used [CII], observed with ALMA, as a tracer of the systemic velocity of the galaxies, and we find that 90% of the selected objects have Lya-[CII] velocity offsets in the range 0 < Dv_Lya-[CII] < 400 km/s, in line with the few measurements available so far in the early Universe, and significantly smaller than those observed at lower z. We observe ISM-[CII] offsets in the range -500 < Dv_ISM-[CII] < 0 km/s, in line with values at all redshifts. We find significant anticorrelations between Dv_Lya-[CII] and the Lya rest-frame equivalent width EW0(Lya) (or equivalently, the Lya escape fraction f_esc(Lya)). According to available models for the radiative transfer of Lya photons, the escape of Lya photons would be favored in galaxies with high outflow velocities, in agreement with our observations. The uniform shell model would also predict that the Lya escape in galaxies with slow outflows (0 < v_out < 300 km/s) is mainly determined by the neutral hydrogen column density (NHI), while the alternative model by Steidel+10 would favor a combination of NHI and covering fraction as driver of the Lya escape. We suggest that the observed increase in Lya escape that is observed between z~2 and z~6 is not due to a higher incidence of fast outflows at high redshift, but rather to a decrease in average NHI along the line of sight, or alternatively, a decrease in HI covering fraction. [abridged], Comment: 22 pages, 6 figures, accepted for publication in A&A. This version replaces the previous one

We report the serendipitous discovery of a bright galaxy (Gal-A) observed as part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE). While this galaxy is detected both in line and continuum emission in ALMA Band 7, it is completely dark in UV/optical filters and only presents a marginal detection in the UltraVISTA Ks band. We discuss the nature of the observed ALMA line, i.e. whether the emission comes from [CII] at z~4.6, or from high-J CO transitions at z~2.2. In the first case we find a [CII]-to-FIR luminosity ratio of log(L_[CII]/L_FIR)=-2.5, consistent with the average value for local star-forming galaxies (SFGs); in the second case, instead, the source would lie outside of the empirical relations between L_CO and L_FIR found in the literature. At both redshifts, we derive the star-formation rate (SFR) from the ALMA continuum, and the stellar mass (M*) by using stellar population synthesis models as input for LePHARE spectral energy distribution (SED) fitting. Exploiting our results, we believe that Gal-A is a Main-Sequence (MS), dusty SFG at z=4.6 (i.e. [CII] emitter) with log(SFR/[M/yr])~1.4 and log(M*/M)~9.7. This work underlines the crucial role of the ALPINE survey in making a census of this class of objects, in order to unveil their contribution to the global star-formation rate density (SFRD) of the Universe at the end of the Reionisation epoch., Comment: 13 pages, 6 figures. Submitted to MNRAS

Using the IRAM NOrthern Extended Millimeter Array (NOEMA), we conducted a program to measure redshifts for 13 bright galaxies detected in the Herschel Astrophysical Large Area Survey (H-ATLAS) with $S_{500{\mu}\rm m}\ge$80 mJy. We report reliable spectroscopic redshifts for 12 individual sources, which are derived from scans of the 3 and 2 mm bands, covering up to 31 GHz in each band, and are based on the detection of at least two emission lines. The spectroscopic redshifts are in the range $2.08 10^{13}\,L_\odot$ galaxies. We also present a reanalysis of the spectral energy distributions including the continuum flux densities measured at 3 and 2 mm to derive the overall properties of the sources. Future prospects based on these efficient measurements of redshifts of high-z galaxies using NOEMA are outlined, including a comprehensive survey of all the brightest Herschel galaxies., Comment: 18 pages, 21 figures