Search

We obtain CO(1-0) molecular gas measurements with ATCA on a sample of 43 spectroscopically confirmed H$\alpha$ emitters in the Spiderweb protocluster at $z=2.16$ and investigate the relation between their star formation and cold gas reservoirs as a function of environment. We achieve a CO(1-0) detection rate of $\sim23\pm12\%$ with 10 dual CO(1-0) and H$\alpha$ detections at $10<\log M_{*}/M_\odot<11.5$. In addition, we obtain upper limits for the remaining sources. In terms of total gas fractions ($F_{gas}$), our sample is divided into two different regimes with a steep transition at $\log M_{*}/M_\odot\approx10.5$. Galaxies below that threshold have gas fractions that in some cases are close to unity, indicating that their gas reservoir has been replenished by inflows from the cosmic web. However, objects at $\log M_{*}/M_\odot>10.5$ display significantly lower gas fractions and are dominated by AGN (12 out of 20). Stacking results yield $F_{gas}\approx0.55$ for massive emitters excluding AGN, and $F_{gas}\approx0.35$ when examining only AGN candidates. Furthermore, depletion times show that most H$\alpha$ emitters may become passive by $1

Massive, intensely star-forming galaxies at high redshift require a supply of molecular gas from their gas reservoirs, replenished by infall from the surrounding circumgalactic medium, to sustain their immense star-formation rates. However, our knowledge of the extent and morphology of their cold-gas reservoirs is still in its infancy. We present the results of stacking 80 hours of JVLA observations of CO(1--0) emission -- which traces the cold molecular gas -- in 19 $z=2.0-4.5$ dusty, star-forming galaxies from the AS2VLA survey. The visibility-plane stack reveals extended emission with a half-light radius of $3.8\pm0.5$~kpc, 2--3$\times$ more extended than the dust-obscured star formation and $1.4\pm0.2\times$ more extended than the stellar emission. Similarly, stacking the [CI](1--0) observations for a subsample of our galaxies yields sizes consistent with CO(1--0). The CO(1--0) size is comparable to the [CII] halos detected around high-redshift star-forming galaxies.The bulk (up to 80\%) of molecular gas resides outside the star-forming region; only a small part of their molecular gas reservoir directly contributes to their current star formation. Photon-dissociation region modelling indicates that the extended CO(1--0) emission arises from clumpy, dense clouds rather than smooth, diffuse gas., Comment: Submitted to A&A. 9 pages, 5 figures

We report the results of the deep and wide Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm mapping of the Spiderweb protocluster at $z=2.16$. The observations were divided into six contiguous fields covering a survey area of 19.3\,arcmin$^2$. With $\sim$13h on-source time, the final maps in the six fields reach the 1$\sigma$ rms noise in a range of $40.3-57.1 \mu$Jy at a spatial resolution of $0.5-0.9$ arcsec. By using different source extraction codes and careful visual inspection, we detect 47 ALMA sources at a significance higher than 4$\sigma$. We construct the differential and cumulative number counts down to $\sim0.2$ mJy after the correction for purity and completeness obtained from Monte Carlo simulations. The ALMA 1.2 mm number counts of dusty star-forming galaxies (DSFGs) in the Spiderweb protocluster are overall two times that of general fields, some fields/regions showing even higher overdensities (more than a factor of 3). This is consistent with the results from previous studies over a larger scale using single-dish instruments. Comparison of the spatial distributions between different populations indicates that our ALMA sources are likely drawn from the same distribution as CO(1-0) emitters from the COALAS large program, but distinct from that of H$\alpha$ emitters. The cosmic SFR density of the ALMA sources is consistent with previous results (e.g. LABOCA 870 $\mu$m observations) after accounting for the difference in volume. We show that molecular gas masses estimates from dust measurements are not consistent with the ones derived from CO(1-0) and thus have to be taken with caution. The multiplicity fraction of single-dish DSFGs is higher than that of the field. Moreover, two extreme concentrations of ALMA sources are found on the outskirts of the Spiderweb protocluster, with an excess of more than 12 times that of general fields., Comment: 22 pages, 13 figures, 6 tables. Accepted for publication in A&A

High-redshift protoclusters consisting of dusty starbursts are thought to play an important role in galaxy evolution. Their dusty nature makes them bright in the FIR/submm but difficult to find in optical/NIR surveys. Radio observations are an excellent way to study these dusty starbursts, as dust is transparent in the radio and there is a tight correlation between the FIR and radio emission of a galaxy. Here, we present MeerKAT 1.28 GHz radio imaging of 3 Herschel candidate protoclusters, with a synthesised beam size of ~$7.5'' \times 6.6''$ and a central thermal noise down to $4.35~\mu$Jy/beam. Our source counts are consistent with other radio counts with no evidence of overdensities. Around 95% of the Herschel sources have 1.28 GHz IDs. Using the Herschel 250 micron primary beam size as the searching radius, we find 54.2% Herschel sources have multiple 1.28 GHz IDs. Our average FIR-radio correlation coefficient $q_{250\mu m}$ is $2.33\pm 0.26$. Adding $q_{250\mu m}$ as a new constraint, the probability of finding chance-aligned sources is reduced by a factor of ~6, but with the risk of discarding true identifications of radio-loud/quiet sources. With accurate MeerKAT positions, we cross-match our Herschel sources to optical/NIR data followed by photometric redshift estimations. By removing $z<1$ sources, the density contrasts of two of the candidate protoclusters increase, suggestive of them being real protoclusters at $z>1$. There is also potentially a $0.9

The Cosmic Dawn Survey (DAWN survey) provides multiwavelength (UV/optical to mid-IR) data across the combined 59 deg$^{2}$ of the Euclid Deep and Auxiliary fields (EDFs and EAFs). Here, the first public data release (DR1) from the DAWN survey is presented. DR1 catalogues are made available for a subset of the full DAWN survey that consists of two Euclid Deep fields: Euclid Deep Field North (EDF-N) and Euclid Deep Field Fornax (EDF-F). The DAWN survey DR1 catalogues do not include $Euclid$ data as they are not yet public for these fields. Nonetheless, each field has been covered by the ongoing Hawaii Twenty Square Degree Survey (H20), which includes imaging from CFHT MegaCam in the new $u$ filter and from Subaru Hyper Suprime-Cam (HSC) in the $griz$ filters. Each field is further covered by $Spitzer$/IRAC 3.6-4.5$\mu$m imaging spanning 10 deg$^{2}$ and reaching $\sim$25 mag AB (5$\sigma$). All present H20 imaging and all publicly available imaging from the aforementioned facilities are combined with the deep $Spitzer$/IRAC data to create source catalogues spanning a total area of 16.87 deg$^{2}$ in EDF-N and 2.85 deg$^{2}$ in EDF-F for this first release. Photometry is measured using The Farmer, a well-validated model-based photometry code. Photometric redshifts and stellar masses are computed using two independent codes for modeling spectral energy distributions: EAZY and LePhare. Photometric redshifts show good agreement with spectroscopic redshifts ($\sigma_{\rm NMAD} \sim 0.5, \eta < 8\%$ at $i < 25$). Number counts, photometric redshifts, and stellar masses are further validated in comparison to the COSMOS2020 catalogue. The DAWN survey DR1 catalogues are designed to be of immediate use in these two EDFs and will be continuously updated. Future data releases will provide catalogues of all EDFs and EAFs and include $Euclid$ data.

Galaxy proto-clusters are receiving an increased interest since most of the processes shaping the structure of clusters of galaxies and their galaxy population are happening at early stages of their formation. The Euclid Survey will provide a unique opportunity to discover a large number of proto-clusters over a large fraction of the sky (14 500 square degrees). In this paper, we explore the expected observational properties of proto-clusters in the Euclid Wide Survey by means of theoretical models and simulations. We provide an overview of the predicted proto-cluster extent, galaxy density profiles, mass-richness relations, abundance, and sky-filling as a function of redshift. Useful analytical approximations for the functions of these properties are provided. The focus is on the redshift range z= 1.5 to 4. We discuss in particular the density contrast with which proto-clusters can be observed against the background in the galaxy distribution if photometric galaxy redshifts are used as supplied by the ESA Euclid mission together with the ground-based photometric surveys. We show that the obtainable detection significance is sufficient to find large numbers of interesting proto-cluster candidates. For quantitative studies, additional spectroscopic follow-up is required to confirm the proto-clusters and establish their richness., Comment: Submitted to Astronomy and Astrophysics, 24 pages, 28 figures

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance., Comment: Accepted for publication in the A&A special issue`Euclid on Sky'

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

Active galactic nuclei (AGNs) are known to exhibit optical/UV variability and most of them can be well modeled by the damped random walks. Physical processes that are not related to the accretion disk, such as tidal disruption events (TDE) or moving foreground dusty clouds, can cause flare-like and eclipse-like features in the optical light curve. Both long-term and high-cadence monitoring are needed to identify such features. By combining the Sloan Digital Sky Survey (SDSS), Panoramic Survey Telescope, and Rapid Response System (Pan-STARRS) with the Zwicky Transient Facility (ZTF) survey, we are able to identify a rare sample (11) out of the SDSS quasar catalog (around 83, 000). These quasars exhibit more or less constant brightness but show rapid optical variation in the ZTF DR2 epochs. To investigate the possible origins of these flare/eclipse-like variabilities, we propose the second epoch spectroscopic observations with the Gran Telescopio CANARIAS (GTC). We find that the change in accretion rate plays a significant role in these quasar variabilities. Among them, we identify two Changing-Look Active Galactic Nuclei (CL-AGN) candidates: SDSS J1427+2930 and SDSS J1420+3757. The luminosity change of the former may be caused by the enhanced SMBH accretion or the tidal disruption event, while the latter is more related to the change in the accretion rate., Comment: 11 pages, 13 figures, 2 tables, accepted for publication in MNRAS

We present a new method for modelling the kinematics of galaxies from interferometric observations by performing the optimization of the kinematic model parameters directly in visibility-space instead of the conventional approach of fitting velocity fields produced with the CLEAN algorithm in real-space. We demonstrate our method on ALMA observations of $^{12}$CO (2$-$1), (3$-$2) or (4$-$3) emission lines from an initial sample of 30 massive 850$\mu$m-selected dusty star-forming galaxies with far-infrared luminosities $\gtrsim$$\,10^{12}\,$L$_{\odot}$ in the redshift range $z \sim\,$1.2$-$4.7. Using the results from our modelling analysis for the 12 sources with the highest signal-to-noise emission lines and disk-like kinematics, we conclude the following: (i) Our sample prefers a CO-to-$H_2$ conversion factor, of $\alpha_{\rm CO} = 0.92 \pm 0.36$; (ii) These far-infrared luminous galaxies follow a similar Tully$-$Fisher relation between the circularized velocity, $V_{\rm circ}$, and baryonic mass, $M_{\rm b}$, as more typical star-forming samples at high redshift, but extend this relation to much higher masses $-$ showing that these are some of the most massive disk-like galaxies in the Universe; (iii) Finally, we demonstrate support for an evolutionary link between massive high-redshift dusty star-forming galaxies and the formation of local early-type galaxies using the both the distributions of the baryonic and kinematic masses of these two populations on the $M_{\rm b}\,-\,\sigma$ plane and their relative space densities.

We present the detailed analysis of the thermal, diffuse emission of the proto-intracluster medium (ICM) detected in the halo of the Spiderweb Galaxy at z=2.16, within a radius of $\sim$ 150 kpc. We combined deep X-ray data from Chandra and millimeter observations of the Sunyaev-Zeldovich (SZ) effect obtained by ALMA. Thanks to independent measurements of the pressure profile from ALMA SZ observation and the electron density profile from the available X-ray data, we derived, for the first time, the temperature profile in the ICM of a z>2 protocluster. It reveals the presence of a strong cool core (comparable to the local ones) that may host a significant mass deposition flow, consistent with measured local star formation values. We also find mild evidence of an asymmetry in the X-ray surface brightness distribution, which may be tentatively associated with a cavity carved into the proto-ICM by the radio jets or, alternatively, may be due to the young dynamical status of the halo. The cooling time of baryons in the core of the Spiderweb Protocluster is estimated to be $\sim$ 0.1 Gyr, implying that the baryon cycle in the first stages of the protocluster formation is characterised by a high-duty cycle and a very active environment. In the case of the Spiderweb protocluster, we are witnessing the presence of a strongly peaked core that is possibily hosting a cooling flow with a mass deposition rate up to 250-1000 $M_{\odot}$/yr, responsible for feeding both the central supermassive black hole and the high star formation rate observed in the Spiderweb Galaxy. This phase is expected to be rapidly followed by active galactic nucleus feedback events, whose onset may have already left an imprint in the radio and X-ray appearance of the Spiderweb protocluster, eventually driving the ICM into a self-regulated, long-term evolution in less than one Gyr.

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

There is controversy whether dusty starbursts selected at submillimeter wavelengths can trace galaxy overdensities. We perform the first systematic search for protoclusters around a homogeneously selected sample of 12 spectroscopically confirmed submillimeter galaxies (SMGs) at $z\sim1.2-5.3$ in the GOODS-N field. We applied the Poisson Probability Method (PPM) to search for Mpc scale overdensities around these SMGs using three photometric redshift catalogs. We detect galaxy overdensities for 11 out of the 12 SMGs ($92\%\pm8$\%), distributed over eight protoclusters. We confirm three previously discovered protoclusters, and we detect five new ones around the SMGs SMMJ123634 ($z=1.225$), ID.19 ($z=2.047$), SMMJ123607 ($z=2.487$), SMMJ123606 ($z=2.505$), and GN10 ($z=5.303$). A wavelet-based analysis shows that the SMGs live in protocluster cores with a complex morphology (compact, filamentary, or clumpy) and an average size of $\sim(0.4-1)$Mpc. By comparing the PPM results obtained using independently the three redshift catalogs, we possibly witness a transitioning phase at $z\gtrsim4$ for the galaxy populations. While $z\lesssim4$ protoclusters appear to be populated by dusty galaxies, those at highest redshifts $z\sim5$ are detected as overdensities of Lyman$\alpha$ emitters or Lyman break galaxies. We also find a good correlation between the molecular (H$_2$) gas mass of the SMG and the overdensity significance. To explain the overall phenomenology, we suggest that galaxy interactions in dense environments likely triggered the starburst and gas-rich phase of the SMGs. Altogether, we support the scenario that SMGs are excellent tracers of distant protoclusters. Those presented in this work are excellent targets for the {\it James Webb Space Telescope.} Surveys with forthcoming facilities (e.g., {\it Euclid}, LSST) can be tuned to detect even larger samples of distant protoclusters., Comment: 22 pages, Astronomy & Astrophysics in press. HTML files of Figure 3 are available at https://obswww.unige.ch/~castigna/Figure3_html/

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 850um imaging of the XMM-LSS field observed for 170 hours as part of the James Clerk Maxwell Telescope SCUBA-2 Large eXtragalactic Survey (S2LXS). S2LXS XMM-LSS maps an area of 9 square degrees, reaching a moderate depth of 1-sigma ~ 4 mJy/beam. This is the largest contiguous area of extragalactic sky mapped by JCMT at 850um to date. The wide area of the S2LXS XMM-LSS survey allows us to probe the ultra-bright (S_850um > 15 mJy), yet rare submillimetre population. We present the S2LXS XMM-LSS catalogue, which comprises 40 sources detected at >5-sigma significance, with deboosted flux densities in the range of 7 mJy to 48 mJy. We robustly measure the bright-end of the 850um number counts at flux densities >7 mJy, reducing the Poisson errors compared to existing measurements. The S2LXS XMM-LSS observed number counts show the characteristic upturn at bright fluxes, expected to be motivated by local sources of submillimetre emission and high-redshift strongly lensed galaxies. We find that the observed 850um number counts are best reproduced by model predictions that include either strong lensing or source blending from a 15 arcsec beam, indicating that both may make an important contribution to the observed over-abundance of bright single dish 850um selected sources. We make the S2LXS XMM-LSS 850um map and >5-sigma catalogue presented here publicly available., Comment: Accepted for publication in MNRAS

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

We present the X-ray imaging and spectral analysis of the diffuse emission around the Spiderweb galaxy at z=2.16 and of its nuclear emission, based on a deep (700 ks) Chandra observation. We characterize the nuclear emission and computed the contamination in the surrounding regions due to the wings of the instrument PSF. Then, we quantified the extended emission within 12". We find that the Spiderweb galaxy hosts a mildly absorbed quasar, with modest yet significant variability on a timescale of ~1 year. We find that the emission in the jet regions is well described by a power law with Gamma~2-2.5, and it is consistent with IC upscattering of the CMB photons by the relativistic electrons. We also find a roughly symmetric, diffuse emission within a radius of ~100 kpc. This emission is consistent with thermal bremsstrahlung from a hot ICM with a temperature of kT=2.0_{-0.4}^{+0.7} keV, and a metallicity of Z<1.6Z_sun. The average electron density within 100 kpc is n_e=(1.51+-0.24+-0.14)E-2 cm^{-3}, corresponding to an upper limit for the total ICM mass of <=(1.76+-0.30+-0.17)E+12 M_sun (where error bars are 1 sigma statistical and systematic, respectively). If we apply hydrostatic equilibrium to the ICM, we measure a total gravitational mass M(<100 kpc)=(1.5^{+0.5}_{-0.3})E+13 M_sun and, extrapolating at larger radii, we estimate a total mass M_{500}=(3.2^{+1.1}_{-0.6})E+13 M_sun within a radius of r_{500}=(220+-30) kpc. We conclude that the Spiderweb protocluster shows significant diffuse emission within a radius of 12 arcsec, whose major contribution is provided by IC scattering associated with the radio jets. Outside the jet regions, we also identified thermal emission within a radius of ~100 kpc, revealing the presence of hot, diffuse baryons that may represent the embryonic virialized halo of the forming cluster., Comment: 19 pages, 16 figures, accepted in Astronomy & Astrophysics

While unobscured and radio-quiet active galactic nuclei are regularly being found at redshifts $z > 6$, their obscured and radio-loud counterparts remain elusive. We build upon our successful pilot study, presenting a new sample of low-frequency-selected candidate high-redshift radio galaxies (HzRGs) over a sky area twenty times larger. We have refined our selection technique, in which we select sources with curved radio spectra between 72-231 MHz from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. In combination with the requirements that our GLEAM-selected HzRG candidates have compact radio morphologies and be undetected in near-infrared $K_{\rm s}$-band imaging from the Visible and Infrared Survey Telescope for Astronomy Kilo-degree Infrared Galaxy (VIKING) survey, we find 51 new candidate HzRGs over a sky area of approximately 1200 deg$^2$. Our sample also includes two sources from the pilot study: the second-most distant radio galaxy currently known, at $z=5.55$, with another source potentially at $z \sim 8$. We present our refined selection technique and analyse the properties of the sample. We model the broadband radio spectra between 74 MHz and 9 GHz by supplementing the GLEAM data with both publicly available data and new observations from the Australia Telescope Compact Array at 5.5 and 9 GHz. In addition, deep $K_{\rm s}$-band imaging from the High-Acuity Widefield $K$-band Imager (HAWK-I) on the Very Large Telescope and from the Southern Herschel Astrophysical Terahertz Large Area Survey Regions $K_{\rm s}$-band Survey (SHARKS) is presented for five sources. We discuss the prospects of finding very distant radio galaxies in our sample, potentially within the epoch of reionisation at $z \gtrsim 6.5$., Comment: 49 pages, 3 figures (one of which is a multi-page figure with 102 separate panels), 9 tables, accepted for publication in PASA

(Abridged) We present an analysis of the 700 ks Chandra ACIS-S observation of the field around the Spiderweb Galaxy at z=2.156, focusing on the nuclear activity in the associated large-scale environment. We identify unresolved X-ray sources down to flux limits of 1.3X10^{-16} and 3.9X10^{-16} erg/s/cm^2 in the soft and hard band, respectively. We search for counterparts in the optical, NIR and submm bands to identify X-ray sources belonging to the protocluster. We detect 107 X-ray unresolved sources within 5 arcmin (corresponding to 2.5 Mpc) of J1140-2629, among which 13 have optical counterparts with spectroscopic redshift 2.11~1.84+-0.04. The best-fit intrinsic absorption for 5 protocluster X-ray members is N_H>10^{23} cm^{-2}, while other 6 have upper limits of the order of fewX10^{22} cm^{-2}. Two sources can only be fitted with very flat \Gamma<=1, and are therefore considered Compton-thick candidates. Their 0.5-10 keV rest frame luminosities are larger than 2X10^{43} erg/s, significantly greater than X-ray luminosities expected from star formation activity. The X-ray luminosity function of AGN in the volume associated to the Spiderweb protocluster in the range 10^{43}10.5, corresponding to an enhancement of 6.0^{+9.0}_{-3.0} with respect to the COSMOS field at comparable redshifts and stellar mass range. We conclude that the galaxy population in the Spiderweb Protocluster is characterized by enhanced X-ray nuclear activity triggered by environmental effects on Mpc scales., Comment: Astronomy & Astrophysics. Minor changes in this version. Version accepted in A&A

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

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

We report an ultra-bright lensed submillimeter galaxy at $z_{spec}=1.4796$, identified as a result of a full-sky cross-correlation of the AllWISE and Planck compact source catalogs aimed to search for bright submillimeter galaxies at $z \sim 1.5-2.8$. APEX/LABOCA observations of the candidate galaxy reveal a source with flux (S$_{870 \mu m}= 54\pm 8$ mJy). The position of the APEX source coincides with the position of the AllWISE mid-IR source, and with the Einstein ring GAL-CLUS-022058s, observed with the HST. Archival VLT/FORS observations reveal the redshift of this Einstein ring, $z_{spec}=1.4796$, and detection of the CO(5-4) line at $z_{spec} = 1.4802$ with APEX/nFLASH230 confirms the redshift of the submillimeter emission. The lensed source appears to be gravitationally magnified by a massive foreground galaxy cluster lens at $z = 0.36$. We use Lenstool to model the gravitational lensing, which is near to a "fold arc" configuration for an elliptical mass distribution of the central halo, where four images of the lensed galaxy are seen; the mean magnification is $\mu_{\rm L} =18\pm 4$. We have determined an intrinsic rest-frame infrared luminosity of $L_{IR} \approx 10^{12} L_\odot $ and a likely star formation rate of $\sim 70-170$ $M_\odot\ yr^{-1}$. The molecular gas mass is $M_{mol} \sim 2.6 \times 10^{10} M_\odot$ and the gas fraction is $f = 0.34\pm 0.07$. We also obtain a stellar mass log$(M_\ast/M_\odot) = 10.7 \pm 0.1$ and a specific star formation rate log$(sSFR/Gyr^{-1})=0.15 \pm 0.03$. This galaxy lies on the so-called main sequence of star-forming galaxies at this redshift.

We present Atacama Large Millimetre Array (ALMA) 2mm continuum observations of a complete and unbiased sample of 99 870micron-selected sub-millimeter galaxies (SMGs) in the Extended Chandra Deep Field South (ALESS). Our observations of each SMG reach average sensitivities of 53 microJy/beam. We measure the flux densities for 70 sources, for which we obtain a typical 870micron-to-2mm flux ratio of 14 +/- 5. We do not find a redshift dependence of this flux ratio, which would be expected if the dust emission properties of our SMGs were the same at all redshifts. By combining our ALMA measurements with existing Herschel/SPIRE observations, we construct a (biased) subset of 27 galaxies for which the cool dust emission is sufficiently well sampled to obtain precise constraints on their dust properties using simple isothermal models. Thanks to our new 2mm observations, the dust emissivity index is well-constrained and robust against different dust opacity assumptions. The median dust emissivity index of our SMGs is $\beta\simeq1.9\pm0.4$, consistent with the emissivity index of dust in the Milky Way and other local and high-redshift galaxies, as well as classical dust grain model predictions. We also find a negative correlation between the dust temperature and $\beta$, similar to low-redshift observational and theoretical studies. Our results indicate that $\beta\simeq2$ in high-redshift dusty star-forming galaxies, implying little evolution in dust grain properties between our SMGs and local dusty galaxy samples, and suggesting these high-mass and high-metallicity galaxies have dust reservoirs driven by grain growth in their ISM., Comment: 31 pages, 20 figures, accepted for publication in the Astrophysical Journal

We present LMT/AzTEC 1.1mm observations of $\sim100$ luminous high-redshift dusty star-forming galaxy candidates from the $\sim600\,$sq.deg $Herschel$-ATLAS survey, selected on the basis of their SPIRE red far-infrared colours and with $S_{500\mu\rm m}=35-80$ mJy. With an effective $\theta_{\rm FWHM}\approx9.5\,$ arcsec angular resolution, our observations reveal that at least 9 per cent of the targets break into multiple systems with SNR $\geq 4$ members. The fraction of multiple systems increases to $\sim23\,$ per cent (or more) if some non-detected targets are considered multiples, as suggested by the data. Combining the new AzTEC and deblended $Herschel$ photometry we derive photometric redshifts, IR luminosities, and star formation rates. While the median redshifts of the multiple and single systems are similar $(z_{\rm med}\approx3.6)$, the redshift distribution of the latter is skewed towards higher redshifts. Of the AzTEC sources $\sim85\,$ per cent lie at $z_{\rm phot}>3$ while $\sim33\,$ per cent are at $z_{\rm phot}>4$. This corresponds to a lower limit on the space density of ultra-red sources at $44$ galaxies., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 23 pages, 18 figures

We report a detailed CO(1-0) survey of a galaxy protocluster field at $z=2.16$, based on 475 hours of observations with the Australia Telescope Compact Array. We constructed a large mosaic of 13 individual pointings, covering an area of 21 arcmin$^2$ and $\pm6500$ km/s range in velocity. We obtain a robust sample of 46 CO(1-0) detections spanning $z=2.09-2.22$, constituting the largest sample of molecular gas measurements in protoclusters to date. The CO emitters show an overdensity at $z=2.12-2.21$, suggesting a galaxy super-protocluster or a protocluster connected to large-scale filaments with ~120 cMpc size. We find that 90% CO emitters have distances $>0'.5-4'$ to the center galaxy, indicating that small area surveys would miss the majority of gas reservoirs in similar structures. Half of the CO emitters have velocities larger than escape velocities, which appears gravitationally unbound to the cluster core. These unbound sources are barely found within the $R_{200}$ radius around the center, which is consistent with a picture in which the cluster core is collapsed while outer regions are still in formation. Compared to other protoclusters, this structure contains relatively more CO emitters with relatively narrow line width and high luminosity, indicating galaxy mergers. We use these CO emitters to place the first constraint on the CO luminosity function and molecular gas density in an overdense environment. The amplitude of the CO luminosity function is 1.6$\pm$0.5 orders of magnitudes higher than observed for field galaxy samples at $z\sim2$, and one order of magnitude higher than predictions for galaxy protoclusters from semi-analytical SHARK models. We derive a high molecular gas density of $0.6-1.3\times10^{9}$ $M_\odot$ cMpc$^{-3}$ for this structure, consistent with predictions for cold gas density of massive structures from hydro-dynamical DIANOGA simulations., Comment: 21 pages + Appendices, 9 figures, 6 tables. Accepted for publication in A&A

We present a conceptual study of a large format imaging spectrograph for the Large Submillimeter Telescope (LST) and the Atacama Large Aperture Submillimeter Telescope (AtLAST). Recent observations of high-redshift galaxies indicate the onset of earliest star formation just a few 100 million years after the Big Bang (i.e., z = 12--15), and LST/AtLAST will provide a unique pathway to uncover spectroscopically-identified first forming galaxies in the pre-reionization era, once it will be equipped with a large format imaging spectrograph. We propose a 3-band (200, 255, and 350 GHz), medium resolution (R = 2,000) imaging spectrograph with 1.5 M detectors in total based on the KATANA concept (Karatsu et al. 2019), which exploits technologies of the integrated superconducting spectrometer (ISS) and a large-format imaging array. A 1-deg2 drilling survey (3,500 hr) will capture a large number of [O III] 88 um (and [C II] 158 um) emitters at z = 8--9, and constrain [O III] luminosity functions at z > 12., Comment: Presented at SPIE Astronomical Telescopes + Instrumentation 2020. Full published paper and video available at http://dx.doi.org/10.1117/12.2561238 See also https://www.lstobservatory.org for the LST project

We report the results of a pilot spectroscopic program of a region at z=5.2 in the GOODS-N field containing an overdensity of galaxies around the well-known submillimeter galaxy HDF850.1. We have selected candidate cluster members from the optical 25 medium-band photometric catalog of the project SHARDS (Survey for High-z Absorption Red and Dead Sources). 17 rest-frame UV selected galaxies (LAEs and LBGs) with 5.1510^15 M_sun)., Comment: Accepted by MNRAS. 20 pages, 11 figures and 3 tables

(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

Dusty star-forming galaxies at high redshift (1 < z < 3) represent the most intense star-forming regions in the Universe. Key aspects to these processes are the gas heating and cooling mechanisms. Although it is well known that these galaxies are gas-rich, little is known about the gas excitation conditions. Here we examine these processes in a sample of 24 strongly lensed star-forming galaxies identified by the \textit{Planck} satellite (LPs) at z ~ 1.1 - 3.5. We analyze 162 CO rotational transitions (ranging from Jupper = 1 - 12) and 37 atomic carbon fine-structure lines ([CI]) in order to characterize the physical conditions of the gas in sample of LPs. We simultaneously fit the CO and [CI] lines, and the dust continuum emission, using two different non-LTE, radiative transfer models. The first model represents a two component gas density, while the second assumes a turbulence driven log-normal gas density distribution. These LPs are among the most gas-rich, infrared (IR) luminous galaxies ever observed ($\mu_{\rm L}$L$_{\rm IR(8-1000\mu m) } \sim 10^{13-14.6} $\Lsun; $< \mu_{\rm L}$M$_{\rm ISM}> = 2.7 \pm 1.2 \times 10^{12}$ \Msun, with $\mu_{\rm L} \sim 10-30$ the average lens magnification factor). Our results suggest that the turbulent ISM present in the LPs can be well-characterized by a high turbulent velocity dispersion ($<\Delta V_{\rm turb}> \sim 100 $ \kms) and gas kinetic temperature to dust temperature ratios $ \sim 2.5$, sustained on scales larger than a few kpc. We speculate that the average surface density of the molecular gas mass and IR luminosity $\Sigma_{\rm M_{\rm ISM}}$ $\sim 10^{3 - 4}$ \Msun pc$^{-2}$ and $\Sigma_{\rm L_{\rm IR}}$ $\sim 10^{11 - 12}$ \Lsun kpc$^{-2}$, arise from both stellar mechanical feedback and a steady momentum injection from the accretion of intergalactic gas., Comment: 54 pages, 18 figures. Accepted for publication in The Astrophysical Journal (as of Oct. 30, 2020). Please feel free to view the supplementary figures here (which can later be found online in the ApJ after the full publication procedure): https://drive.google.com/drive/folders/1CN3rqlUDcNi5JSDH2GhFk_SLSlo10han?usp=sharing

We analyse the physical properties of 121 SNR $\geq$ 5 sub-millimetre galaxies (SMGs) from the STUDIES 450-$\mu$m survey. We model their UV-to-radio spectral energy distributions using MAGPHYS+photo-$z$ and compare the results to similar modelling of 850-$\mu$m-selected SMG sample from AS2UDS, to understand the fundamental physical differences between the two populations at the observed depths. The redshift distribution of the 450-$\mu$m sample has a median of $z$ = 1.85 $\pm$ 0.12 and can be described by strong evolution of the far-infrared luminosity function. The fainter 450-$\mu$m sample has $\sim$14 times higher space density than the brighter 850-$\mu$m sample at $z$ $\lesssim$2, and a comparable space density at $z$ = 2-3, before rapidly declining, suggesting LIRGs are the main obscured population at $z$ $\sim$ 1-2, while ULIRGs dominate at higher redshifts. We construct rest-frame $\sim$ 180-$\mu$m-selected and dust-mass-matched samples at $z$ = 1-2 and $z$ = 3-4 from the 450-$\mu$m and 850-$\mu$m samples, respectively, to probe the evolution of a uniform sample of galaxies spanning the cosmic noon era. Using far-infrared luminosity, dust masses and an optically-thick dust model, we suggest that higher-redshift sources have higher dust densities due to inferred dust continuum sizes which are roughly half of those for the lower-redshift population at a given dust mass, leading to higher dust attenuation. We track the evolution in the cosmic dust mass density and suggest that the dust content of galaxies is governed by a combination of both the variation of gas content and dust destruction timescale., Comment: 20 pages, 7 figures, MNRAS accepted

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

Understanding how galaxies cease to form stars represents an outstanding challenge for galaxy evolution theories. This process of "star formation quenching" has been related to various causes, including Active Galactic Nuclei (AGN) activity, the influence of large-scale dynamics, and the environment in which galaxies live. In this paper, we present the first results from a follow-up of CALIFA survey galaxies with observations of molecular gas obtained with the APEX telescope. Together with EDGE survey CARMA observations, we collect $^{12}$CO observations that cover approximately one effective radius in 472 CALIFA galaxies. We observe that the deficit of galaxy star formation with respect to the star formation main sequence (SFMS) increases with the absence of molecular gas and with a reduced efficiency of conversion of molecular gas into stars, in line with results of other integrated studies. However, by dividing the sample into galaxies dominated by star formation and galaxies quenched in their centres (as indicated by the average value of the H$\alpha$ equivalent width), we find that this deficit increases sharply once a certain level of gas consumption is reached, indicating that different mechanisms drive separation from the SFMS in star-forming and quenched galaxies. Our results indicate that differences in the amount of molecular gas at a fixed stellar mass are the primary driver for the dispersion in the SFMS, and the most likely explanation for the start of star-formation quenching. However, once a galaxy is quenched, changes in star formation efficiency drive how much a retired galaxy separates in star formation rate from star-forming ones of similar masses. In other words, once a paucity of molecular gas has significantly reduced star formation, changes in the star formation efficiency are what drives a galaxy deeper into the red cloud, retiring it., Comment: 12 pages, 4 figures, 1 table

We present a survey of the molecular gas in 61 submillimetre galaxies (SMGs) selected from 870$\mu$m continuum surveys of the COSMOS, UDS and ECDFS fields, using the Atacama Large Millimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA). 46 $^{12}$CO ($J=$2-5) emission lines are detected in 45 of the targets at $z=$1.2-4.8, with redshifts indicating that those which are submillimetre bright and undetected/faint in the optical/near-infrared typically lie at higher redshifts, with a gradient of $\Delta z/\Delta S_{870}=$0.11$\pm$0.04mJy$^{-1}$. We also supplement our data with literature sources to construct a statistical CO spectral line energy distribution and find the $^{12}$CO line luminosities in SMGs peak at $J_{\rm up}\sim$6, consistent with the Cosmic Eyelash, among similar studies. Our SMGs lie mostly on or just above the main sequence, displaying a decrease in their gas depletion timescales $t_{\rm dep} = M_{\rm gas}/{\rm SFR}$ with redshift in the range $z\sim$1-5 and a median of 200$\pm$50Myr at $z\sim$2.8. This coincides with an increase in molecular gas fraction $\mu_{\rm gas} = M_{\rm gas}/M_\ast$ across the same redshift range. Finally we demonstrate that the $M_{\rm baryon}$-$\sigma$ distribution of our SMGs is consistent with that followed by early-type galaxies in the Coma cluster, providing strong support to the suggestion that SMGs are progenitors of massive local spheroidal galaxies. On the basis of this we suggest that the SMG populations above and below an 870-$\mu$m flux limit of $S_{870}\sim$5mJy may correspond to the division between slow- and fast-rotators seen in local early-type galaxies., Comment: Submitted to MNRAS

We present new gas kinematic observations with the OSIRIS instrument at the GTC for galaxies in the Cl1604 cluster system at z=0.9. These observations together with a collection of other cluster samples at different epochs analyzed by our group are used to study the evolution of the Tully-Fisher, velocity-size and stellar mass-angular momentum relations in dense environments over cosmic time. We use 2D and 3D spectroscopy to analyze the kinematics of our cluster galaxies and extract their maximum rotation velocities (Vmax). Our methods are consistently applied to all our cluster samples which make them ideal for an evolutionary comparison. Up to redshift one, our cluster samples show evolutionary trends compatible with previous observational results in the field and in accordance with semianalytical models and hydrodynamical simulations concerning the Tully-Fisher and velocity-size relations. However, we find a factor 3 drop in disk sizes and an average B-band luminosity enhancement of 2 mag by z=1.5. We discuss the role that different cluster-specific interactions may play in producing this observational result. In addition, we find that our intermediate-to-high redshift cluster galaxies follow parallel sequences with respect to the local specific angular momentum-stellar mass relation, although displaying lower angular momentum values in comparison with field samples at similar redshifts. This can be understood by the stronger interacting nature of dense environments with respect to the field., Comment: 26 pages, 13 figures. Accepted for publication in A&A

Using the EMIR instrument on the IRAM 30m telescope, we conducted a spectroscopic redshift search of seven z$_{\rm phot}$ $\sim$ 4 sub-millimetre bright galaxies selected from the Herschel Bright Sources (HerBS) sample with fluxes at 500 $\mu$m greater than 80 mJy. For four sources, we obtained spectroscopic redshifts between 3.4 < z < 4.1 through the detection of multiple CO-spectral lines with J $\leq$ 3. Later, we detected low-J transitions for two of these sources with the GBT including the CO(1-0) transition. For the remaining three sources, more data are needed to determine the spectroscopic redshift unambiguously. The measured CO luminosities and line widths suggest that all these sources are gravitationally lensed. These observations demonstrate that the 2 mm window is indispensable to confirm robust spectroscopic redshifts for z < 4 sources. Finally, we present an efficient graphical method to correctly identify spectroscopic redshifts., Comment: 20 Pages, 8 Figures + 7 Appendix Figures + 1 animated figure. Resubmitted in order to reflect the version as printed in MNRAS

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

Using data from ALMA and near-infrared (NIR) integral field spectrographs including both SINFONI and KMOS on the VLT, we investigate the two-dimensional distributions of H$\alpha$ and rest-frame far-infrared (FIR) continuum in six submillimeter galaxies at $z\sim2$. At a similar spatial resolution ($\sim$0.5" FWHM; $\sim$4.5 kpc at $z=2$), we find that the half-light radius of H$\alpha$ is significantly larger than that of the FIR continuum in half of the sample, and on average H$\alpha$ is a median factor of $2.0\pm0.4$ larger. Having explored various ways to correct for the attenuation, we find that the attenuation-corrected H$\alpha$-based SFRs are systematically lower than the IR-based SFRs by at least a median factor of $3\pm1$, which cannot be explained by the difference in half-light radius alone. In addition, we find that in 40% of cases the total $V$-band attenuation ($A_V$) derived from energy balance modeling of the full ultraviolet(UV)-to-FIR spectral energy distributions (SEDs) is significantly higher than that derived from SED modeling using only the UV-to-NIR part of the SEDs, and the discrepancy appears to increase with increasing total infrared luminosity. Finally, considering all our findings along with the studies in the literature, we postulate that the dust distributions in SMGs, and possibly also in less IR luminous $z\sim2$ massive star-forming galaxies, can be decomposed into three main components; the diffuse dust heated by older stellar populations, the more obscured and extended young star-forming HII regions, and the heavily obscured central regions that have a low filling factor but dominate the infrared luminosity in which the majority of attenuation cannot be probed via UV-to-NIR emissions., Comment: 15 pages, 8 figures, A&A in press

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

We present ALMA Band 6 (nu=233GHz, lambda=1.3mm) continuum observations towards 68 'normal' star-forming galaxies within two Coma-like progenitor structures at z=2.10 and 2.47, from which ISM masses are derived, providing the largest census of molecular gas mass in overdense environments at these redshifts. Our sample comprises galaxies with a stellar mass range of 1x10^9M_sun - 4x10^11M_sun with a mean M_*~6x10^10M_sun. Combining these measurements with multiwavelength observations and SED modeling, we characterize the gas mass fraction and the star formation efficiency, and infer the impact of the environment on galaxies' evolution. Most of our detected galaxies (~70%) have star formation efficiencies and gas fractions similar to those found for coeval field galaxies and in agreement with the field scaling relations. However, we do find that the proto-clusters contain an increased fraction of massive, gas-poor galaxies, with low gas fractions (f_gas<6-10%) and red rest-frame ultraviolet/optical colors typical of post-starburst and passive galaxies. The relatively high abundance of passive galaxies suggests an accelerated evolution of massive galaxies in proto-cluster environments. The large fraction of quenched galaxies in these overdense structures also implies that environmental quenching takes place during the early phases of cluster assembly, even before virialization. From our data, we derive a quenching efficiency of E_q~0.45 and an upper limit on the quenching timescale of T_q<1Gyr., Comment: Accepted for publication in The Astrophysical Journal

Dusty star-forming galaxies (DSFGs) detected at $z > 4$ provide important examples of the first generations of massive galaxies. However, few examples with spectroscopic confirmation are currently known, with Hershel struggling to detect significant numbers of $z > 6$ DSFGs. NGP6_D1 is a bright 850 $\mu m$ source (12.3 $\pm$ 2.5 mJy) with no counterparts at shorter wavelengths (a SPIRE dropout). Interferometric observations confirm it is a single source, with no evidence for any optical or NIR emission, or nearby likely foreground lensing sources. No $>3\sigma$ detected lines are seen in both LMT RSR and IRAM 30m EMIR spectra of NGP6_D1 across 32 $GHz$ of bandwidth despite reaching detection limits of $\sim 1 mJy/500 km~s^{-1}$, so the redshift remains unknown. Template fitting suggests that NGP6_D1 is most likely between $z = 5.8$ and 8.3. SED analysis finds that NGP6_D1 is a ULIRG, with a dust mass $\sim 10^8$ - $10^9$ $M_{\odot}$ and a SFR of $\sim$ 500 $M_{\odot}~yr^{-1}$. We place upper limits on the gas mass of NGP6_D1 of $M_{H2}$ $ < (1.1~\pm~3.5) \times 10^{11}$ $M_{\odot}$, consistent with a gas-to-dust ratio of $\sim$ 100 - 1000. We discuss the nature of NGP6_D1 in the context of the broader submm population, and find that comparable SPIRE dropouts account for $\sim$ 20% of all SCUBA-2 detected sources, but with a similar flux density distribution to the general population., Comment: Accepoted for publication in MNRAS

We present SCUBA-2 850-$\mu$m observations of 13 candidate starbursting protoclusters selected using Planck and Herschel data. The cumulative number counts of the 850-$\mu$m sources in 9/13 of these candidate protoclusters show significant overdensities compared to the field, with the probability $<$10$^{-2}$ assuming the sources are randomly distributed in the sky. Using the 250-, 350-, 500- and 850-$\mu$m flux densities, we estimate the photometric redshifts of individual SCUBA-2 sources by fitting spectral energy distribution (SED) templates with an MCMC method. The photometric redshift distribution, peaking at $22$ protoclusters and the peak of the cosmic star-formation rate density (SFRD). We find that the 850-$\mu$m sources in our candidate protoclusters have infrared luminosities of $L_{\mathrm{IR}}\gtrsim$10$^{12}L_{\odot}$ and star-formation rates of SFR=(500-1,500)$M_{\odot}$yr$^{-1}$. By comparing with results in the literature considering only Herschel photometry, we conclude that our 13 candidate protoclusters can be categorised into four groups: six of them being high-redshift starbursting protoclusters, one being a lower-redshift cluster/protocluster, three being protoclusters that contain lensed DSFG(s) or are rich in 850-$\mu$m sources, and three regions without significant Herschel or SCUBA-2 source overdensities. The total SFRs of the candidate protoclusters are found to be comparable or higher than those of known protoclusters, suggesting our sample contains some of the most extreme protocluster population. We infer that cross-matching Planck and Herschel data is a robust method for selecting candidate protoclusters with overdensities of 850-$\mu$m sources., Comment: 30 pages, 8 figures, 5 tables, accepted for publication in MNRAS

We present spectra and near-infrared images of a sample of faint radio sources initially selected as promising high-redshift radio galaxy (HzRG) candidates. We have determined redshifts for a total of 13 radio galaxies with redshifts ranging from $0.52\le z \le 5.72$. Our sample probes radio luminosities that are almost an order of magnitude fainter than previous large samples at the highest redshifts. We use near-infrared photometry for a subsample of these galaxies to calculate stellar masses using simple stellar population models, and find stellar masses to be in the range $10^{10.8} - 10^{11.7} M_\odot$. We then compare our faint radio galaxies with brighter radio galaxies at $z\ge2$ from the literature. We find that fainter radio galaxies have lower Ly$\alpha$ luminosities and narrower full width at half maxima compared to the bright ones, implying photoionisation by weaker AGN. We also rule out the presence of strong shocks in faint HzRGs. The stellar masses determined for faint HzRGs are lower than what has been observed for brighter ones. We find that faint HzRG population in the redshift range $2-4$ forms a bridge between star-forming and narrow-line AGN, whereas the ones at $z>4$ are dominated by star-formation, and may be building up their stellar mass through cold accretion of gas. Finally, we show that the overall redshift evolution of radio sizes at $z>2$ is fully compatible with increased inverse Compton scattering losses at high redshifts., Comment: 26 pages (including 2 appendices), 17 figures, accepted for publication in MNRAS

ALMA 870$\mu$m continuum imaging has uncovered a population of blends of multiple dusty star-forming galaxies (DSFGs) in sources originally detected with the Herschel Space Observatory. However, their pairwise separations are much smaller that what is found by ALMA follow-up of other single-dish surveys or expected from theoretical simulations. Using ALMA and VLA, we have targeted three of these systems to confirm whether the multiple 870$\mu$m continuum sources lie at the same redshift, successfully detecting $^{12}$CO($J = 3$-2) and $^{12}$CO($J = 1$-0) lines and being able to confirm that in the three cases all the multiple DSFGs are likely physically associated within the same structure. Therefore, we report the discovery of two new gas-rich dusty protocluster cores (HELAISS02, $z = 2.171 \pm 0.004$; HXMM20, $z = 2.602 \pm 0.002$). The third target is located in the well known COSMOS overdensity at $z = 2.51$ (named CL J1001+0220 in the literature), for which we do not find any new secure CO(1-0) detection, although some of its members show only tentative detections and require further confirmation. From the gas, dust, and stellar properties of the two new protocluster cores, we find very large molecular gas fractions yet low stellar masses, pushing the sources above the main sequence, while not enhancing their star formation efficiency. We suggest that the sources might be newly formed galaxies migrating to the main sequence. The properties of the three systems compared to each other and to field galaxies may suggest a different evolutionary stage between systems., Comment: Accepted by ApJ. 24 pages, 9 figures

ALMA measurements for 93 $Herschel$-selected galaxies at $1.1 \leqslant z \leqslant 1.7$ in COSMOS reveal a sizable ($>29$\%) population with compact star formation (SF) sizes, lying on average $> \times 3.6$ below the optical stellar mass ($M_{\star}$)-size relation of disks. This sample widely spans the star-forming Main Sequence (MS), having $10^{8} \leqslant M_{\star} \leqslant 10^{11.5} \ M_{\odot}$ and $20 \leqslant SFR \leqslant 680 \ M_{\odot} \rm yr^{-1}$. The 32 size measurements and 61 upper limits are measured on ALMA images that combine observations of CO(5-4), CO(4-3), CO(2-1) and $\lambda_{\rm obs} \sim 1.1-1.3 \ \rm mm$ continuum, all tracing the star-forming molecular gas. These compact galaxies have instead normally extended $K_{band}$ sizes, suggesting strong specific $SFR$ gradients. Compact galaxies comprise the $50\pm18 \%$ of MS galaxies at $M_{\star} > 10^{11} M_{\odot}$. This is not expected in standard bi-modal scenarios where MS galaxies are mostly steadily-growing extended disks. We suggest that compact MS objects are early post-starburst galaxies in which the merger-driven boost of SF has subsided. They retain their compact SF size until either further gas accretion restores pre-merger galaxy-wide SF, or until becoming quenched. The fraction of merger-affected SF inside the MS seems thus larger than anticipated and might reach $\sim 50$\% at the highest $M_{\star}$. The presence of large galaxies above the MS demonstrates an overall poor correlation between galaxy SF size and specific $SFR$., Comment: 7 pages, 5 figures. Accepted for publication in ApJ Letters

High-redshift ($z$) proto-clusters will serve as testing grounds to probe the gas supply furnishing the emerging metals, stars, and large-scale structures we see at the current epoch. This work focuses on the major role large radio/millimeter (mm) single dish facilities will have in constraining the bulk, cold (T $= 10^{1-4}$K) molecular and atomic gas content. To highlight the need for large radio/mm single dishes, we calculate how the high-sensitivity of the Green Bank Telescope's (GBT) unblocked 100m aperture provides vital interferometric short-spacing coverage to support higher-resolution ngVLA observations of the cold neutral gas at the largest scales. These combined observations are optimal for revealing low-surface brightness emission, and thus aid in the total baryonic mass estimates across cosmic time., Comment: Submitted for the March 11, 2019 Astro2020 Decadal Survey deadline