Search

Photometric redshifts are widely used in studies of dusty star-forming galaxies (DSFGs), but catastrophic photo-$z$ failure can undermine all redshift-dependent results. Here we report the spectroscopic redshift confirmation of COSBO-7, a strongly lensed DSFG in the COSMOS-PRIMER field. Recently, a photometric redshift solution of $z\gtrsim7.0$ was reported for COSBO-7 based on ten bands of {\it James Webb} Space Telescope (JWST) NIRCam and MIRI imaging data. This $z$ value was favored by four independent spectral energy distribution (SED) fitting codes, and the result provided an appealing candidate for the most distant massive DSFG known to date. This photo-$z$ solution was also supported by a single line detection in Atacama Large Millimeter Array (ALMA) Band 3 consistent with CO(7-6) at $z=7.46$. However, our new ALMA observations robustly detect two lines in Band 6 identified as CO(7-6) and [CI](2-1) at $z_{\rm spec}=2.625$, and thus the Band 3 line as CO(3-2). These three robust line detections decidedly place COSBO-7 at $z=2.625$, refuting the photo-$z$ solution. We derive physical parameters by fitting near-infrared(NIR)-to-millimeter(mm) photometry and lens modeling, revealing that COSBO-7 is a main sequence galaxy. We examine possible reasons for this photo-$z$ failure and attribute it to (1) the likely underestimation of photometric uncertainties at 0.9\,$\mu$m and 1.15 \,$\mu$m; and (2) the lack of photometry at wavelengths beyond 20\,$\mu$m. Notably, we recover a bona fide $z_{\rm phot}\sim 2.3$ by including the existing MIPS $24\,\mu$m photometry, demonstrating the critical importance of mid-infrared (MIR) data in bolstering photo-$z$ measurements. This work highlights a common challenge in modeling the SEDs of DSFGs, and provides a cautionary tale regarding the reliability of photometric redshifts as well as pseudo-spectroscopic redshifts based on single line detection., Comment: Published in A&A

The NOEMA formIng Cluster survEy (NICE) is a large program targeting 69 massive galaxy group candidates at $z>2$ in six deep fields. We report spectroscopic confirmation of eight groups at $1.65\leq z\leq3.61$ in COSMOS. Homogeneously selected as significant overdensities of red IRAC sources with red Herschel colors, four groups are confirmed by CO and [CI] with NOEMA 3mm observations, three are confirmed with ALMA, and one is confirmed by H$\alpha$ from Subaru/FMOS. We constructed the integrated FIR SEDs for the eight groups, obtaining total IR SFR $=260-1300~{\rm M_\odot}$~yr$^{-1}$. We adopted six methods to estimate the dark matter masses, including stellar mass to halo mass relations, overdensity with galaxy bias, and NFW profile fitting to radial stellar mass density. We found the radial stellar mass density are consistent with a NFW profile, supporting that they are collapsed structures hosted by a single dark matter halo. The best halo mass estimates are $\log(M_{\rm h}/{\rm M_\odot})=12.8-13.7$ with uncertainty of 0.3 dex. From halo mass estimates, we derive baryonic accretion rate ${\rm BAR}=(1-8)\times10^{3}\,{\rm M_{\odot}/yr}$ for this sample. We find a quasi-linear correlation between the integrated SFR/BAR and the theoretical halo mass limit for cold streams, $M_{\rm stream}/M_{\rm h}$, with ${\rm SFR/BAR}=10^{-0.46\pm0.22}\left({M_{\rm stream}/M_{\rm h}}\right)^{0.71\pm0.16}$ with a scatter of $0.40\,{\rm dex}$. Further, we compare halo masses and stellar masses with simulations, and find all structures are consistent with being progenitors of $M_{\rm h}(z=0)>10^{14}\,{\rm M_{\odot}}$ galaxy clusters, and the most massive central galaxies have stellar masses consistent with brightest cluster galaxies (BCGs) progenitors in the TNG300 simulation. The results strongly suggest these structures are forming massive galaxy clusters via baryonic and dark matter accretion., Comment: 44 pages (27pp appendix), 32 figures, 18 tables, accepted for publication in A&A

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

Strong gravitational magnification by massive galaxy clusters enable us to detect faint background sources, resolve their detailed internal structures, and in the most extreme cases identify and study individual stars in distant galaxies. Highly magnified individual stars allow for a wide range of applications, including studies of stellar populations in distant galaxies and constraining small-scale dark matter structures. However, these applications have been hampered by the small number of events observed, as typically one or a few stars are identified from each distant galaxy. Here, we report the discovery of 46 significant microlensed stars in a single strongly-lensed high-redshift galaxy behind the Abell 370 cluster at redshift of 0.725 when the Universe was half of its current age (dubbed the ``Dragon arc''), based on two observations separated by one year with the James Webb Space Telescope ({\it JWST}). These events are mostly found near the expected lensing critical curves, suggesting that these are magnified individual stars that appear as transients from intracluster stellar microlenses. Through multi-wavelength photometry and colors, we constrain stellar types and find that many of them are consistent with red giants/supergiants magnified by factors of thousands. This finding reveals an unprecedented high occurrence of microlensing events in the Dragon arc, and proves that {\it JWST}'s time-domain observations open up the possibility of conducting statistical studies of high-redshift stars and subgalactic scale perturbations in the lensing dark matter field., Comment: 15 pages, 4 figures, 1 table submitted to Nature Astronomy

Sub/millimeter galaxies are a key population for the study of galaxy evolution because the majority of star formation at high redshifts occurred in galaxies deeply embedded in dust. To search for this population, we have performed an extensive survey with ALMA, called the ALMA Lensing Cluster Survey (ALCS). This survey covers 133 arcmin^2 area and securely detects 180 sources at z~0.5-6 with a flux limit of ~0.2 mJy at 1.2 mm (Fujimoto et al. 2023). Here we report the results of multi-wavelength spectral energy distribution (SED) analysis of the whole ALCS sample, utilizing the observed-frame UV to millimeter photometry. We find that the majority of the ALCS sources lie on the star-forming main sequence, with a smaller fraction showing intense starburst activities. The ALCS sample contains high infrared-excess sources IRX=log(Ldust/LUV)>1), including two extremely dust-obscured galaxies (IRX>5). We also confirm that the ALCS sample probes a broader range in lower dust mass than conventional SMG samples in the same redshift range. We identify six heavily obscured AGN candidates that are not detected in the archival Chandra data in addition to the three X-ray AGNs reported by Uematsu et al. (2023). The inferred AGN luminosity density shows a possible excess at z=2-3 compared with that determined from X-ray surveys below 10 keV., Comment: 27 pages, 20 figures

Our knowledge of relations between supermassive black holes and their host galaxies at $z\gtrsim1$ is still limited, even though being actively sought out to $z\sim6$. Here, we use the high resolution and sensitivity of JWST to measure the host galaxy properties for 107 X-ray-selected type-I AGNs at $0.68

Understanding the gas content in galaxies, its consumption and replenishment, remains pivotal in our comprehension of the evolution of the Universe. Numerous studies have addressed this, utilizing various observational tools and analytical methods. These include examining low-transition $^{12}$CO millimeter rotational lines and exploring the far-infrared and the (sub-)millimeter emission of galaxies. With the capabilities of present-day facilities, much of this research has been centered on relatively bright galaxies. We aim at exploring the gas reservoirs of a more general type of galaxy population at $1.0\leq z\leq 3.0$. We stack ALMA 1.1 mm data to measure the gas content of a mass-complete sample down to $\sim10^{8.6}$ M$_\odot$ at $z=1$ ($\sim10^{9.2}$ M$_\odot$ at $z=3$), extracted from the HST/CANDELS sample in GOODS-S. The sample is composed of 5,530 on average blue ($\sim0.12$ mag, $\sim0.81$ mag), star-forming main sequence objects ($\Delta$MS$\sim-0.03$). We report measurements at $10^{10-11}$ M$_\odot$ and upper limits for the gas fractions at $10^{8-10}$ M$_\odot$. At $10^{10-11}$ M$_\odot$, our f$_{\mathrm{gas}}$, ranging from 0.32 to 0.48, agree well with other studies based on mass-complete samples down to $10^{10}$ M$_\odot$, and are lower than expected according to other works more biased to individual detections. At $10^{9-10}$ M$_\odot$, we obtain 3$\sigma$ upper limits for f$_{\mathrm{gas}}$ ranging from 0.69 to 0.77. These upper limits are on the level of the extrapolations of scaling relations based on mass-complete samples down to $10^{10}$ M$_\odot$. As such, it suggests that the gas content of low-mass galaxies is at most what is extrapolated from literature scaling relations. The comparison of our results with previous works reflects how the inclusion of bluer, less obscured, and more MS-like objects progressively pushes the gas level to lower values., Comment: 18 pages, 8 figures, 4 tables. Accepted for publication in A&A

We present observations and analysis of the starburst, PACS-819, at z=1.45 ($M_*=10^{10.7}$ M$_{ \odot}$), using high-resolution ($0^{\prime \prime}.1$; 0.8 kpc) ALMA and multi-wavelength JWST images from the COSMOS-Web program. Dissimilar to HST/ACS images in the rest-frame UV, the redder NIRCam and MIRI images reveal a smooth central mass concentration and spiral-like features, atypical for such an intense starburst. Through dynamical modeling of the CO J=5--4 emission with ALMA, PACS-819 is rotation-dominated thus has a disk-like nature. However, kinematic anomalies in CO and asymmetric features in the bluer JWST bands (e.g., F150W) support a more disturbed nature likely due to interactions. The JWST imaging further enables us to map the distribution of stellar mass and dust attenuation, thus clarifying the relationships between different structural components, not discernable in the previous HST images. The CO J = 5 -- 4 and FIR dust continuum emission are co-spatial with a heavily-obscured starbursting core (<1 kpc) which is partially surrounded by much less obscured star-forming structures including a prominent arc, possibly a tidally-distorted dwarf galaxy, and a clump, either a sign of an ongoing violent disk instability or a recently accreted low-mass satellite. With spatially-resolved maps, we find a high molecular gas fraction in the central area reaching $\sim3$ ($M_{\text{gas}}$/$M_*$) and short depletion times ($M_{\text{gas}}/SFR\sim$ 120 Myrs) across the entire system. These observations provide insights into the complex nature of starbursts in the distant universe and underscore the wealth of complementary information from high-resolution observations with both ALMA and JWST., Comment: 21 pages(including appendix), 13 figures, accepted by ApJ

We report the discovery of a large-scale structure at z=3.44 revealed by JWST data in the Extended Groth Strip (EGS) field. This structure, called the Cosmic Vine, consists of 20 galaxies with spectroscopic redshifts at 3.4370\%$). Comparisons with simulations suggest that the Cosmic Vine would form a cluster with halo mass $M_{\rm halo}>10^{14}M_\odot$ at z=0, and the two massive galaxies are likely forming the brightest cluster galaxies (BCGs). The results unambiguously reveal that massive quiescent galaxies can form in growing large-scale structures at z>3, thus disfavoring the environmental quenching mechanisms that require a virialized cluster core. Instead, as suggested by the interacting and bulge-dominated morphologies, the two galaxies are likely quenched by merger-triggered starburst or active galactic nucleus (AGN) feedback before falling into a cluster core. Moreover, we found that the observed specific star formation rates of massive quiescent galaxies in z>3 dense environments are one to two orders of magnitude lower than that of the BCGs in the TNG300 simulation. This discrepancy potentially poses a challenge to the models of massive cluster galaxy formation. Future studies comparing a large sample with dedicated cluster simulations are required to solve the problem., Comment: A&A Letters in press, open access at https://doi.org/10.1051/0004-6361/202348540

The study of distant galaxy groups and clusters at the peak epoch of star formation is limited by the lack of a statistically and homogeneously selected and spectroscopically confirmed sample. Recent discoveries of concentrated starburst activities in cluster cores have opened a new window to hunt for these structures based on their integrated IR luminosities. Hereby we carry out the large NOEMA (NOrthern Extended Millimeter Array) program targeting a statistical sample of infrared-luminous sources associated with overdensities of massive galaxies at z>2, the Noema formIng Cluster survEy (NICE). We present the first result from the ongoing NICE survey, a compact group at z=3.95 in the Lockman Hole field (LH-SBC3), confirmed via four massive (M_star>10^10.5M_sun) galaxies detected in CO(4-3) and [CI](1-0) lines. The four CO-detected members of LH-SBC3 are distributed over a 180 kpc physical scale, and the entire structure has an estimated halo mass of ~10^13Msun and total star formation rate (SFR) of ~4000Msun/yr. In addition, the most massive galaxy hosts a radio-loud AGN with L_1.4GHz, rest = 3.0*10^25W/Hz. The discovery of LH-SBC3 demonstrates the feasibility of our method to efficiently identify high-z compact groups or forming cluster cores. The existence of these starbursting cluster cores up to z~4 provides critical insights into the mass assembly history of the central massive galaxies in clusters., Comment: 10 pages, 8 figures, published by A&A

We report galaxy MACS0416-Y3 behind the lensing cluster MACSJ0416.1--2403 as a tentative rotating disk at $z=8.34$ detected through its [OIII]$\lambda5007$ emission in JWST NIRCam wide-field slitless spectroscopic observations. The discovery is based on our new grism dynamical modeling methodology for JWST NIRCam slitless spectroscopy, using the data from ``Median-band Astrophysics with the Grism of NIRCam in Frontier Fields'' (MAGNIF), a JWST Cycle-2 program. The [OIII]$\lambda5007$ emission line morphology in grism data shows velocity offsets compared to the F480M direct imaging, suggestive of rotation. Assuming a geometrically thin disk model, we constrain the rotation velocity of $v_{\rm rot}=58^{+53}_{-35}$ km s$^{-1}$ via forward modeling of the two-dimensional (2D) spectrum. We obtain the kinematic ratio of $v_{\rm rot}/\sigma_v=1.6^{+1.9}_{-0.9}$, where $\sigma_v$ is the velocity dispersion, in line with a quasi-stable thin disk. The resulting dynamical mass is estimated to be $\log(M_{\rm dyn}/M_{\odot})=8.4^{+0.5}_{-0.7}$. If the rotation confirmed, our discovery suggests that rotating gaseous disks may have already existed within 600 million years after Big Bang., Comment: 15 pages, 6 figures. Comments welcome

We present the results of a systematic study of the rest-frame optical morphology of quiescent galaxies at $z \geq 3$ using the Near-Infrared Camera (NIRCam) onboard $JWST$. Based on a sample selected by $UVJ$ color or $NUVUVJ$ color, we focus on 26 quiescent galaxies with $9.8<\log{(M_\star/M_\odot)}<11.4$ at $2.810.3$ implies that our size - stellar mass relations are below those at lower redshifts, with the amplitude of $\sim0.6\, {\rm kpc}$ at $M_\star = 5\times 10^{10}\, M_\odot$. This value agrees with the extrapolation from the size evolution of quiescent galaxies at $z<3$ in the literature, implying that the size of quiescent galaxies increases monotonically from $z\sim3-5$. Our sample is mainly composed of galaxies with bulge-like structures according to their median S\'ersic index and axis ratio of $n\sim3-4$ and $q\sim0.6-0.8$, respectively. On the other hand, there is a trend of increasing fraction of galaxies with low S\'ersic index, suggesting $3

We present NOEMA CO (2-1) line and ALMA 870 $\mu$m continuum observations of a main-sequence galaxy at $z=1.38$. The galaxy was initially selected as a "gas-giant", based on the gas mass derived from sub-mm continuum (log$(M_{\rm gas}/M_{\odot})=11.20\pm0.20$), however the gas mass derived from CO (2-1) luminosity brings down the gas mass to a value consistent with typical star-forming galaxies at that redshift (log$(M_{\rm gas}/M_{\odot})=10.84\pm0.03$). Despite that the dust-to-stellar mass ratio remains elevated above the scaling relations by a factor of 5. We explore the potential physical picture and consider an underestimated stellar mass and optically thick dust as possible causes. Based on the updated gas-to-stellar mass ratio we rule out the former, and while the latter can contribute to the dust mass overestimate it is still not sufficient to explain the observed physical picture. Instead, possible explanations include enhanced HI reservoirs, CO-dark H$_2$ gas, an unusually high metallicity, or the presence of an optically dark, dusty contaminant. Using the ALMA data at 870 $\mu$m coupled with $HST$/ACS imaging, we find extended morphology in dust continuum and clumpy star-formation in rest-frame UV in this galaxy, and a tentative $\sim 10$ kpc dusty arm is found bridging the galaxy center and a clump in F814W image. The galaxy shows levels of dust obscuration similar to the so-called $HST$-dark galaxies at higher redshifts, and would fall into the optically faint/dark $JWST$ color-color selection at $z>2$. It is therefore possible that our object could serve as low-$z$ analog of the $HST$-dark populations. This galaxy serves as a caveat to the gas masses based on the continuum alone, with a larger sample required to unveil the full picture., Comment: 9 pages, 5 figures, 1 table. Submitted to A&A Letters

We present new CO ($J=5-4$ and $7-6$) and [CI] ($^3P_2\,-\, ^3P_1$ and $^3P_1\,-\, ^3P_0$) emission line observations of the star-forming galaxy D49 at the massive end of the Main Sequence at $z=3$. We incorporate previous CO ($J=3-2$) and optical-to-millimetre continuum observations to fit its spectral energy distribution (SED). Our results hint at high-$J$ CO luminosities exceeding the expected location on the empirical correlations with the infrared luminosity. [CI] emission fully consistent with the literature trends is found. We do not retrieve any signatures of a bright active galactic nucleus that could boost the $J=5-4,\,7-6$ lines in either the infrared or X-ray bands, but warm photon-dominated regions, shocks or turbulence could in principle do so. We suggest that mechanical heating could be a favourable mechanism able to enhance the gas emission at fixed infrared luminosity in D49 and other main-sequence star-forming galaxies at high redshift, but further investigation is necessary to confirm this explanation. We derive molecular gas masses from dust, CO, and [CI] that all agree within the uncertainties. Given its large star formation rate (SFR) $\sim 500~M_\odot~{\rm yr}^{-1}$ and stellar mass $>10^{11.5}~M_\odot$, the short depletion time scale of $<0.3$ Gyr might indicate that D49 is experiencing its last growth spurt and will soon transit to quiescence., Comment: 7 pages, 4 figures, accepted for publication in A&A Letter

We present a search for extremely red, dust-obscured, $z>7$ galaxies with $\textit{JWST}$/NIRCam+MIRI imaging over the first 20 arcmin$^2$ of publicly-available Cycle 1 data from the COSMOS-Web, CEERS, and PRIMER surveys. Based on their red color in F277W$-$F444W ($\sim 2.5$ mag) and detection in MIRI/F770W ($\sim 25$ mag), we identify two galaxies$\unicode{x2014}$COS-z8M1 and CEERS-z7M1$\unicode{x2014}$which have best-fit photometric redshifts of $z=8.5^{+0.3}_{-0.4}$ and $z=7.6^{+0.1}_{-0.1}$, respectively. We perform SED fitting with a variety of codes (including BAGPIPES, PROSPECTOR, BEAGLE, and CIGALE), and find a $>95\%$ probability that these indeed lie at $z>7$. Both sources are compact ($R_{\rm eff} \lesssim 200$ pc), highly obscured ($A_V \sim 1.5$$\unicode{x2013}$$2.5$), and, at our best-fit redshift estimates, likely have strong [OIII]+H$\beta$ emission contributing to their $4.4\,\mu$m photometry. We estimate stellar masses of $\sim 10^{10}~M_\odot$ for both sources; by virtue of detection in MIRI at $7.7\,\mu$m, these measurements are robust to the inclusion of bright emission lines, for example, from an AGN. We identify a marginal (2.9$\sigma$) ALMA detection at 2 mm within $0.5''$ of COS-z8M1, which if real, would suggest a remarkably high IR luminosity of $\sim 10^{12} L_\odot$. These two galaxies, if confirmed at $z\sim 8$, would be extreme in their stellar and dust masses, and may be representative of a substantial population of modestly dust-obscured galaxies at cosmic dawn., Comment: 23 pages, 8 figures, submitted to ApJL, comments welcome

In recent years, observations have uncovered a population of massive galaxies that are invisible or very faint in deep optical/near-infrared (near-IR) surveys but brighter at longer wavelengths. However, the nature of these optically dark or faint galaxies (OFGs; one of several names given to these objects) is highly uncertain. In this work, we investigate the drivers of dust attenuation in the JWST era. In particular, we study the role of stellar mass, size, and orientation in obscuring star-forming galaxies (SFGs) at $3 < z < 7.5$, focusing on the question of why OFGs and similar galaxies are so faint at optical/near-IR wavelengths. We find that stellar mass is the primary proxy for dust attenuation, among the properties studied. Effective radius and axis ratio do not show a clear link with dust attenuation, with the effect of orientation being close to random. However, there is a subset of highly dust attenuated ($A_V > 1$, typically) SFGs, of which OFGs are a specific case. For this subset, we find that the key distinctive feature is their compact size (for massive systems with $\log (M_{*}/M_{\odot}) > 10$); OFGs exhibit a 30% smaller effective radius than the average SFG at the same stellar mass and redshift. On the contrary, OFGs do not exhibit a preference for low axis ratios (i.e., edge-on disks). The results in this work show that stellar mass is the primary proxy for dust attenuation and compact stellar light profiles behind the thick dust columns obscuring typical massive SFGs., Comment: Accepted by A&A. 13 pages, 9 figures

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

We report the discovery of a triply imaged active galactic nucleus (AGN), lensed by the galaxy cluster MACS J0035.4-2015 ($z_{\mathrm{d}}=0.352$). The object is detected in Hubble Space Telescope imaging taken for the RELICS program. It appears to have a quasi-stellar nucleus consistent with a point-source, with a de-magnified radius of $r_e\lesssim100$ pc. The object is spectroscopically confirmed to be an AGN at $z_{\mathrm{spec}}=2.063\pm0.005$ showing broad rest-frame UV emission lines, and is detected in both X-ray observations with Chandra and in ALCS ALMA band 6 (1.2 mm) imaging. It has a relatively faint rest-frame UV luminosity for a quasar-like object, $M_{\mathrm{UV},1450}=-19.7\pm0.2$. The object adds to just a few quasars or other X-ray sources known to be multiply lensed by a galaxy cluster. Some diffuse emission from the host galaxy is faintly seen around the nucleus and there is a faint object nearby sharing the same multiple-imaging symmetry and geometric redshift, possibly an interacting galaxy or a star-forming knot in the host. We present an accompanying lens model, calculate the magnifications and time delays, and infer physical properties for the source. We find the rest-frame UV continuum and emission lines to be dominated by the AGN, and the optical emission to be dominated by the host galaxy of modest stellar mass $M_{\star}\simeq10^{9.2} \mathrm{M}_{\odot}$. We also observe some variation in the AGN emission with time, which may suggest that the AGN used to be more active. This object adds a low-redshift counterpart to several relatively faint AGN recently uncovered at high redshifts with HST and JWST., Comment: Published in MNRAS. Corrected miss-assignment of affiliations in the author list

We present the results of a systematic search for candidate quiescent galaxies in the distant Universe in eleven $JWST$ fields with publicly available observations collected during the first three months of operations and covering an effective sky area of $\sim145$ arcmin$^2$. We homogeneously reduce the new $JWST$ data and combine them with existing observations from the $Hubble\,Space\,Telescope$. We select a robust sample of $\sim80$ candidate quiescent and quenching galaxies at $3 < z < 5$ using two methods: (1) based on their rest-frame $UVJ$ colors, and (2) a novel quantitative approach based on Gaussian Mixture Modeling of the $NUV-U$, $U-V$, and $V-J$ rest-frame color space, which is more sensitive to recently quenched objects. We measure comoving number densities of massive ($M_\star\geq 10^{10.6} M_\odot$) quiescent galaxies consistent with previous estimates relying on ground-based observations, after homogenizing the results in the literature with our mass and redshift intervals. However, we find significant field-to-field variations of the number densities up to a factor of $2-3$, highlighting the effect of cosmic variance and suggesting the presence of overdensities of red quiescent galaxies at $z>3$, as it could be expected for highly clustered massive systems. Importantly, $JWST$ enables the robust identification of quenching/quiescent galaxy candidates at lower masses and higher redshifts than before, challenging standard formation scenarios. All data products, including the literature compilation, are made publicly available., Comment: 15 pages, 6 Figures + Appendix. Accepted for publication in ApJ on Feb, 9. Data release: - Reduced HST+JWST mosaics + photometric catalogs and Eazy-py modeling: https://erda.ku.dk/archives/7166d013c1ca1371aac3c57b9e73190d/published-archive.html - Supplementary material and tables: https://zenodo.org/record/7614908#.Y-4ZruzMLmE - MAST: https://doi.org/10.17909/g3nt-a370 - See also Gould et al. 2023

We report the multi-wavelength properties of millimeter galaxies hosting X-ray detected active galactic nuclei (AGNs) from the ALMA Lensing Cluster Survey (ALCS). ALCS is an extensive survey of well-studied lensing clusters with ALMA, covering an area of 133 arcmin$^2$ over 33 clusters with a 1.2 mm flux-density limit of ${\sim}$60 $\mathrm{\mu Jy}$ ($1\sigma$). Utilizing the archival data of Chandra, we identify three AGNs at $z=$1.06, 2.09, and 2.84 among the 180 millimeter sources securely detected in the ALCS (of which 155 are inside the coverage of Chandra). The X-ray spectral analysis shows that two AGNs are not significantly absorbed ($\log N_{\mathrm{H}}/\mathrm{cm}^{-2} < 23$), while the other shows signs of moderate absorption ($\log N_{\mathrm{H}}/\mathrm{cm}^{-2}\sim 23.5$). We also perform spectral energy distribution (SED) modelling of X-ray to millimeter photometry. We find that our X-ray AGN sample shows both high mass accretion rates (intrinsic 0.5--8 keV X-ray luminosities of ${\sim}10^{\text{44--45}}\,\mathrm{erg\ s^{-1}}$) and star-formation rates (${\gtrsim}100\,M_{\odot}\,\mathrm{yr}^{-1}$). This demonstrates that a wide-area survey with ALMA and Chandra can selectively detect intense growth of both galaxies and supermassive black holes (SMBHs) in the high-redshift universe., Comment: 14 pages, 3 figures, 2 tables

Using the novel $JWST$/NIRCam observations in the Abell 2744 field, we present a first spatially resolved overview of an $HST$-dark galaxy, spectroscopically confirmed at $z=2.58$ with magnification $\mu\approx1.9$. While being largely invisible at $\sim$1 $\mu$m with NIRCam, except for sparse clumpy sub-structures, the object is well-detected and resolved in the long-wavelength bands with a spiral shape clearly visible in F277W. By combining ancillary ALMA and $Herschel$ data, we infer that this object is an edge-on dusty spiral with an intrinsic stellar mass log$(M_*/M_\odot)\sim11.3$ and a dust-obscured SFR$\sim 300~M_\odot$~yr$^{-1}$. A massive quiescent galaxy (log$(M_*/M_\odot)\sim10.8$) with tidal features lies 2\farcs{0} away ($r$$\sim$9 kpc), at a consistent redshift as inferred by $JWST$ photometry, indicating a potential major merger. The dusty spiral lies on the main-sequence of star formation, and shows high dust attenuation in the optical ($34$ would be only detectable in F356W and F444W in UNCOVER-like survey, and become totally $JWST$-dark at $z\sim6$. This suggests that detecting highly attenuated galaxies in the Epoch of Reionization might be a challenging task for $JWST$., Comment: 15 pages, 5 figures, 1 table. Accepted to ApJL

We report the discovery of a compact group of galaxies, CGG-z5, at z~5.2 in the EGS field covered by the JWST/CEERS survey. CGG-z5 was selected as the highest overdensity of galaxies at z>2 in recent JWST public surveys and it consists of six candidate members lying within a projected area of $1.5"\times3"$ (10$\times$20~kpc$^2$). All group members are HST/F435W and HST/F606W dropouts while securely detected in the JWST/NIRCam bands, yielding a narrow range of robust photometric redshifts $5.011$) at z~1. This implies that CGG-z5 could be a "proto-massive galaxy" captured during a short-lived phase of massive galaxy formation., Comment: Published in A&A Letter

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

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

We report a galaxy group candidate HPC1001 at $z\approx3.7$ in the COSMOS field. This structure was selected as a high galaxy overdensity at $z>3$ in the COSMOS2020 catalog. It contains ten candidate members, of which eight are assembled in a $10''\times10''$ area with the highest sky density among known protoclusters and groups at $z>3$. Four out of ten sources were also detected at 1.2$~$mm with Atacama Large Millimeter Array continuum observations. Photometric redshifts, measured by four independent methods, fall within a narrow range of $3.5

We present a set of multi-wavelength mosaics and photometric catalogs in the ALMA lensing cluster survey (ALCS) fields. The catalogs were built by reprocessing of archival data from the CHArGE compilation, taken by the $\textit{Hubble Space Telescope}$ ($\textit{HST}$) in the RELICS, CLASH and Hubble Frontier Fields. Additionally we have reconstructed the $\textit{Spitzer}$ IRAC 3.6 and 4.5 $\mu$m mosaics, by utilising all the available archival IRSA/SHA exposures. To alleviate the effect of blending in such a crowded region, we have modelled the $\textit{Spitzer}$ photometry by convolving the $\textit{HST}$ detection image with the $\textit{Spitzer}$ PSF using the novel $\texttt{golfir}$ software. The final catalogs contain 218,000 sources, covering a combined area of 690 arcmin$^2$. These catalogs will serve as an important tool in aiding the search of the sub-mm galaxies in future ALMA surveys, as well as follow ups of the $\textit{HST}$ dark - IRAC sources. Coupled with the available $\textit{HST}$ photometry the addition of the 3.6 and 4.5 $\mu$m bands will allow us to place a better constraint on photometric redshifts and stellar masses of these objects, thus giving us an opportunity to identify high-redshift candidates for spectroscopic follow ups and answer the important questions regarding the epoch of reionization and formation of first galaxies., Comment: 41 pages, 16 figures, 8 tables. Accepted to ApJS. Mosaics and photometric catalogs can be accessed online https://github.com/dawn-cph/alcs-clusters

Using NOEMA and ALMA 3mm line scans, we measure spectroscopic redshifts of six new dusty galaxies at 3.56, based on template IR spectrum energy distribution (SED) from known submillimeter galaxies at z=4--6. Dust SED analyses explain the photo-z overestimate from seemingly cold dust temperatures (Td) and steep Rayleigh-Jeans (RJ) slopes, providing additional examples of cold dusty galaxies impacted by the Cosmic Microwave Background (CMB). We therefore study the general properties of the enlarged sample of 10 ``cold" dusty galaxies over 3.5~4, together with the severe impact of the CMB on their RJ observables, paving the way for the diagnostics of optically thick dust in the early universe. Conventional gas mass estimates based on RJ dust continuum luminosities implicitly assume an optically thin case, overestimating gas masses by a factor of 2--3 on average in compact dusty star-forming galaxies., Comment: Published version in A&A

We present a detailed study of a galaxy merger taking place at $z=1.89$ in the GOODS-S field. Here we analyze Keck/MOSFIRE spectroscopic observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey along with multi-wavelength photometry assembled by the 3D-HST survey. The combined dataset is modeled to infer the past star-formation histories (SFHs) of both merging galaxies. They are found to be massive, with log$_{10}(M_{\ast}/M_{\odot}) > 11$, with a close mass ratio satisfying the typical major-merger definition. Additionally, in the context of delayed-$\tau$ models, GOODS-S 43114 and GOODS-S 43683 have similar SFHs and low star-formation rates (log$_{10}$(SFR(SED)/$M_{\odot}/\rm{yr}^{-1}$) $<$ 1.0) compared to their past averages. The best-fit model SEDs show elevated H$\delta_{\rm{A}}$ values for both galaxies, indicating that their stellar spectra are dominated by A-type stars, and that star formation peaked $\sim0.5-1$ Gyr ago and has recently declined. Additionally, based on SED fitting both merging galaxies turned on and shut off star formation within a few hundred Myr of each other, suggesting that their bursts of star formation may be linked. Combining the SFHs and H$\delta_{\rm{A}}$ results with recent galaxy merger simulations, we infer that these galaxies have recently completed their first pericentric passage and are moving apart. Finally, the relatively low second velocity moment of GOODS-S 43114 given its stellar mass, suggests a disk-like structure. However, including the geometry of the galaxy in the modeling does not completely resolve the discrepancy between the dynamical and stellar masses. Future work is needed to resolve this inconsistency in mass., Comment: 12 pages, 7 figures, 1 table, Accepted for publication in MNRAS

We investigate the fine-structure [C${\rm \scriptsize II}$] line at $158\,\mu$m as a molecular gas tracer by analyzing the relationship between molecular gas mass ($M_{\rm mol}$) and [C${\rm \scriptsize II}$] line luminosity ($L_{\rm [CII]}$) in 11,125 $z\simeq 6$ star-forming, main sequence galaxies from the SIMBA simulations, with line emission modeled by S\'IGAME. Though most ($\sim 50-100\,\%$) of the gas mass in our simulations is ionized, the bulk ($> 50\,\%$) of the [C${\rm \scriptsize II}$] emission comes from the molecular phase. We find a sub-linear (slope $0.78\pm 0.01$) $\log L_{\rm [CII]}-\log M_{\rm mol}$ relation, in contrast with the linear relation derived from observational samples of more massive, metal-rich galaxies at $z \lesssim 6$. We derive a median [C${\rm \scriptsize II}$]-to-$M_{\rm mol}$ conversion factor of $\alpha_{\rm [CII]} \simeq 18\,{\rm M_{\rm \odot}/L_{\rm \odot}}$. This is lower than the average value of $\simeq 30\,{\rm M_{\rm \odot}/L_{\rm \odot}}$ derived from observations, which we attribute to lower gas-phase metallicities in our simulations. Thus, a lower, luminosity-dependent, conversion factor must be applied when inferring molecular gas masses from [C${\rm \scriptsize II}$] observations of low-mass galaxies. For our simulations, [C${\rm \scriptsize II}$] is a better tracer of the molecular gas than CO $J=1-0$, especially at the lowest metallicities, where much of the gas is 'CO-dark'. We find that $L_{\rm [CII]}$ is more tightly correlated with $M_{\rm mol}$ than with star-formation rate (${\rm SFR}$), and both the $\log L_{\rm [CII]}-\log M_{\rm mol}$ and $\log L_{\rm [CII]}-\log {\rm SFR}$ relations arise from the Kennicutt-Schmidt relation. Our findings suggest that $L_{\rm [CII]}$ is a promising tracer of the molecular gas at the earliest cosmic epochs., Comment: 13 pages, 9 figures. Accepted for publication in ApJ

Observations of cold molecular gas reservoirs are critical for understanding the shutdown of star formation in massive galaxies. While dust continuum is an efficient and affordable tracer, this method relies upon the assumption of a "normal" molecular-gas to dust mass ratio, $\delta_{\mathrm{GDR}}$, typically of order one hundred. Recent null detections of quiescent galaxies in deep dust continuum observations support a picture where the cold gas and dust has been rapidly depleted or expelled. In this work, we present another viable explanation: a significant fraction of galaxies with low star formation per unit stellar mass are predicted to have extreme $\delta_{\mathrm{GDR}}$ ratios. We show that simulated massive quiescent galaxies at $0 < z < 3$ in the \textsc{simba} cosmological simulations have $\delta_{\mathrm{GDR}}$ values that extend $>$4 orders of magnitude. The dust in most simulated quiescent galaxies is destroyed significantly more rapidly than the molecular gas depletes, and cannot be replenished. The transition from star-forming to quiescent halts dust formation via star formation processes, with dust subsequently destroyed by supernova shocks and thermal sputtering of dust grains embedded in hot plasma. After this point, the dust growth rate in the models is not sufficient to overcome the loss of $>$3 orders of magnitude in dust mass to return to normal values of $\delta_{\mathrm{GDR}}$ despite having high metallicity. Our results indicate that it is not straight forward to use a single observational indicator to robustly pre-select exotic versus normal ratios. These simulations make strong predictions that can be tested with millimeter facilities., Comment: 9 pages, 5 figures, version accepted for publication in the Astrophysical Journal Letters

The 2mm Mapping Obscuration to Reionization with ALMA (MORA) Survey was designed to detect high redshift ($z\gtrsim4$), massive, dusty star-forming galaxies (DSFGs). Here we present two, likely high redshift sources, identified in the survey whose physical characteristics are consistent with a class of optical/near-infrared (OIR) invisible DSFGs found elsewhere in the literature. We first perform a rigorous analysis of all available photometric data to fit spectral energy distributions and estimate redshifts before deriving physical properties based on our findings. Our results suggest the two galaxies, called MORA-5 and MORA-9, represent two extremes of the "OIR-dark" class of DSFGs. MORA-5 ($z_{\rm phot}=4.3^{+1.5}_{-1.3}$) is a significantly more active starburst with a star-formation rate of 830$^{+340}_{-190}$M$_\odot$yr$^{-1}$ compared to MORA-9 ($z_{\rm phot}=4.3^{+1.3}_{-1.0}$) whose star-formation rate is a modest 200$^{+250}_{-60}$M$_\odot$yr$^{-1}$. Based on the stellar masses (M$_{\star}\approx10^{10-11}$M$_\odot$), space density ($n\sim(5\pm2)\times10^{-6}$Mpc$^{-3}$, which incorporates two other spectroscopically confirmed OIR-dark DSFGs in the MORA sample at $z=4.6$ and $z=5.9$), and gas depletion timescales ($<1$Gyr) of these sources, we find evidence supporting the theory that OIR-dark DSFGs are the progenitors of recently discovered $3

We present the characteristics of 2mm-selected sources from the largest Atacama Large Millimeter and submillimeter Array (ALMA) blank-field contiguous survey conducted to-date, the Mapping Obscuration to Reionization with ALMA (MORA) survey covering 184arcmin$^2$ at 2mm. Twelve of the thirteen detections above 5$\sigma$ are attributed to emission from galaxies, eleven of which are dominated by cold dust emission. These sources have a median redshift of $\langle z_{\rm 2mm}\rangle=3.6^{+0.4}_{-0.3}$ primarily based on optical/near-infrared (OIR) photometric redshifts with some spectroscopic redshifts, with 77$\pm$11% of sources at $z>3$ and 38$\pm$12% of sources at $z>4$. This implies that 2mm selection is an efficient method for identifying the highest redshift dusty star-forming galaxies (DSFGs). Lower redshift DSFGs ($z<3$) are far more numerous than those at $z>3$ yet likely to drop out at 2mm. MORA shows that DSFGs with star-formation rates in excess of 300M$_\odot$ yr$^{-1}$ and relative rarity of $\sim$10$^{-5}$ Mpc$^{-3}$ contribute $\sim$30% to the integrated star-formation rate density between $32$. Analysis of MORA sources' spectral energy distributions hint at steeper empirically-measured dust emissivity indices than typical literature studies, with $\langle\beta\rangle=2.2^{+0.5}_{-0.4}$. The MORA survey represents an important step in taking census of obscured star-formation in the Universe's first few billion years, but larger area 2mm surveys are needed to more fully characterize this rare population and push to the detection of the Universe's first dusty galaxies., Comment: 42 pages, 14 figures, accepted for publication in ApJ

Star formation in half of massive galaxies was quenched by the time the Universe was three billion years old. Very low amounts of molecular gas appear responsible for this, at least in some cases, though morphological gas stabilization, shock heating, or activity associated with accretion onto a central supermassive black hole is invoked in other cases. Recent studies of quenching by gas depletion have been based upon upper limits that are insufficiently sensitive to determine this robustly, or stacked emission with its problems of averaging. Here we report 1.3mm observations of dust emission from six strongly lensed galaxies where star formation has been quenched, with magnifications of up to a factor of 30. Four of the six galaxies are undetected in dust emission, with an estimated upper limit on the dust mass of 0.0001 times the stellar mass, and by proxy (assuming a Milky Way molecular gas-to-dust ratio) 0.01 times the stellar mass in molecular gas. This is two orders of magnitude less molecular gas per unit stellar mass than seen in star forming galaxies at similar redshifts. It remains difficult to extrapolate from these small samples, but these observations establish that gas depletion is responsible for a cessation of star formation in some fraction of high-redshift galaxies., Comment: 17 pages, 3 figures. Authors' version. Published online by Nature on September 22, 2021

In the framework of a systematic ALMA study of IR-selected main-sequence and starburst galaxies at z~1-1.7 at typical ~1" resolution, we report on the effects of mid-IR- and X-ray-detected active galactic nuclei (AGN) on the reservoirs and excitation of molecular gas in a sample of 55 objects. We find detectable nuclear activity in ~30% of the sample. The presence of dusty tori influences the IR SED of galaxies, as highlighted by the strong correlation among the AGN contribution to the total IR luminosity budget (fAGN = LIR,AGN/LIR), its hard X-ray emission, and the Rayleigh-Jeans to mid-IR (S1.2mm/S24um) observed color, with consequences on the empirical SFR estimates. Nevertheless, we find only marginal effects of AGN on the CO (J=2,4,5,7) or neutral carbon ([CI](1-0), [CI](2-1)) line luminosities and on the derived molecular gas excitation as gauged by line ratios and the full SLEDs. The [CI] and CO emission up to J=5,7 thus primarily traces the properties of the host in typical IR luminous galaxies. However, we highlight the existence of a large variety of line luminosities and ratios despite the homogeneous selection. In particular, we find a sparse group of AGN-dominated sources with the highest LIR,AGN/LIR,SFR ratios, >3, that are more luminous in CO(5-4) than what is predicted by the L'CO(5-4)-LIR,SFR relation, which might be the result of the nuclear activity. For the general population, our findings translate into AGN having minimal effects on quantities such as gas and dust fractions and SFEs. If anything, we find hints of a marginal tendency of AGN hosts to be compact at far-IR wavelengths and to display 1.8x larger dust optical depths. In general, this is consistent with a marginal impact of the nuclear activity on the gas reservoirs and star formation in average star-forming AGN hosts with LIR>5e11 Lsun, typically underrepresented in surveys of quasars and SMGs., Comment: 16 pages, 10 figures + Appendix, accepted in A&A on August 28th, 2021. The data compilation will be available on Vizier. (Minor) updates for the AGN population with respect to the data release of Valentino et al. 2020, A&A, 641, 155. For an early access, please contact the corresponding author. Abstract slightly modified to adjust to arXiv's requirements

We present a Spitzer/IRAC survey of H-faint ($H_{160} \gtrsim 26.4$, $<5\sigma$) sources in 101 lensing cluster fields. Across a CANDELS/Wide-like survey area of $\sim$648 arcmin$^2$ (effectively $\sim$221 arcmin$^2$ in the source plane), we have securely discovered 53 sources in the IRAC Channel-2 band (CH2, 4.5 $\mathrm{\mu m}$; median CH2$=22.46\pm0.11$ AB mag) that lack robust HST/WFC3-IR F160W counterparts. The most remarkable source in our sample, namely ES-009 in the field of Abell 2813, is the brightest H-faint galaxy at 4.5 $\mathrm{\mu m}$ known so far ($\mathrm{CH2}=20.48\pm0.03$ AB mag). We show that the H-faint sources in our sample are massive (median $M_\mathrm{star} = 10^{10.3\pm 0.3}$ $M_\odot$), star-forming (median star formation rate $=100_{-40}^{+60}$ $M_\odot$yr$^{-1}$) and dust-obscured ($A_V=2.6\pm0.3$) galaxies around a median photometric redshift of $z=3.9\pm0.4$. The stellar continua of 14 H-faint galaxies can be resolved in the CH2 band, suggesting a median circularized effective radius ($R_\mathrm{e,circ}$; lensing corrected) of $1.9\pm0.2$ kpc and $<1.5$ kpc for the resolved and whole samples, respectively. This is consistent with the sizes of massive unobscured galaxies at $z\sim4$, indicating that H-faint galaxies represent the dusty tail of the distribution of a wider galaxy population. Comparing with the ALMA dust continuum sizes of similar galaxies reported previously, we conclude that the heavy dust obscuration in H-faint galaxies is related to the compactness of both stellar and dust continua ($R_\mathrm{e,circ}\sim 1$ kpc). These H-faint galaxies make up $16_{-7}^{+13}$% of the galaxies in the stellar mass range of $10^{10}-10^{11.2}$ $M_\odot$ at $z=3\sim5$, contributing to $8_{-4}^{+8}$% of the cosmic star formation rate density in this epoch and likely tracing the early phase of massive galaxy formation., Comment: 26 pages, 10 figures, 3 tables. Accepted for publication in ApJ

We characterise the basic far-IR (FIR) properties and the gas mass fraction of massive ( ~ 11.0) quiescent galaxies (QGs) and explore how these evolve from z = 2.0 to the present day. We use robust, multi-wavelength (mid- to far-IR and sub-millimetre to radio) stacking ensembles of homogeneously selected and mass complete samples of log(M*/Msun) > 10.8 QGs. We find that the dust to stellar mass ratio (Md/M*) rises steeply as a function of redshift up to z~1.0 and then remains flat at least out to z = 2.0. Using Md as a proxy of gas mass (Mgas), we find a similar trend for the evolution of the gas mass fraction (fgas) with z > 1.0 QGs having fgas ~ 7.0% (for solar metallicity). This fgas is 3 - 10 times lower than that of normal star forming galaxies (SFGs) at their corresponding redshift but ~3 and ~10 times larger compared to that of z = 0.5 and local QGs. Furthermore, the inferred gas depletion time scales are comparable to that of local SFGs and systematically longer than that of main sequence galaxies at their corresponding redshifts. Our analysis also reveals that the average dust temperature (Td) of massive QGs remains roughly constant (< Td > = 21.0 \pm 2.0K) at least out to z ~ 2.0 and is substantially colder (~ 10K) compared to that of z > 0 SFGs. This motivated us to construct and release a redshift-invariant template IR SED, that we use to make predictions for ALMA observations and to explore systematic effects in the Mgas estimates of massive, high-z QGs. Finally, we discuss how a simple model that considers progenitor-bias can effectively reproduce the observed evolution of Md/M* and fgas. Our results indicate universal initial interstellar medium conditions for quenched galaxies and a large degree of uniformity in their internal processes across cosmic time., Comment: Accepted for publication in A&A

We present an analysis of the chemical abundance properties of $\approx$650 star-forming galaxies at $z \approx0.6-1.8$. Using integral-field observations from the $K$-band Multi-Object Spectrograph (KMOS), we quantify the [NII]/H$\alpha$ emission-line ratio, a proxy for the gas-phase Oxygen abundance within the interstellar medium. We define the stellar mass-metallicity relation at $z \approx0.6-1.0$ and $z \approx1.2-1.8$ and analyse the correlation between the scatter in the relation and fundamental galaxy properties (e.g. H$\alpha$ star-formation rate, H$\alpha$ specific star-formation rate, rotation dominance, stellar continuum half-light radius and Hubble-type morphology). We find that for a given stellar mass, more highly star-forming, larger and irregular galaxies have lower gas-phase metallicities, which may be attributable to their lower surface mass densities and the higher gas fractions of irregular systems. We measure the radial dependence of gas-phase metallicity in the galaxies, establishing a median, beam smearing-corrected, metallicity gradient of $ \Delta Z/ \Delta R=0.002 \pm0.004$ dex kpc$^{-1}$, indicating on average there is no significant dependence on radius. The metallicity gradient of a galaxy is independent of its rest-frame optical morphology, whilst correlating with its stellar mass and specific star-formation rate, in agreement with an inside-out model of galaxy evolution, as well as its rotation dominance. We quantify the evolution of metallicity gradients, comparing the distribution of $\Delta Z/ \Delta R$ in our sample with numerical simulations and observations at $z \approx0-3$. Galaxies in our sample exhibit flatter metallicity gradients than local star-forming galaxies, in agreement with numerical models in which stellar feedback plays a crucial role redistributing metals., Comment: Accepted for publication in MNRAS. Comments welcome. 20 pages, 12 figures. Revised version with updated introduction and FMR evolution following comments

We present a method for recovering the intrinsic (extinction-corrected) luminosity of the 11.2 micron PAH band in galaxy spectra. Using 105 high S/N Spitzer/IRS spectra of star-forming galaxies, we show that the equivalent width ratio of the 12.7 and 11.2 micron PAH bands is independent on the optical depth, with small dispersion of ~5% indicative of a nearly constant intrinsic flux ratio R_int = (f_12.7/f_11.2)_int = 0.377 +/- 0.020. Conversely, the observed flux ratio, R_obs = (f_12.7/f_11.2)_obs strongly correlates with the silicate strength (S_sil) confirming that differences in R_obs reflect variation in the optical depth. The relation between R_obs and S_sil reproduces predictions for the Galactic Centre extinction law but disagrees with other laws. We calibrate the total extinction affecting the 11.2 micron PAH from R_obs, which we apply to another sample of 215 galaxies with accurate measurements of the total infrared luminosity (L_IR) to investigate the impact of extinction on L_11.2/L_IR. Correlation between L_11.2/L_IR and R_obs independently on L_IR suggests that increased extinction explains the well known decrease in the average L_11.2/L_IR at high L_IR. The extinction-corrected L_11.2 is proportional to L_IR in the range L_IR/L_sun = 10^9--10^13. These results consolidate L_11.2 as a robust tracer of star formation in galaxies., Comment: 14 pages, 9 figures, 2 tables. Accepted for publication in MNRAS

We thoroughly explore the properties of (sub)-millimeter (mm) selected galaxies (SMGs) in the Shark semi-analytic model of galaxy formation. Compared to observations, the predicted number counts at wavelengths (lambda) 0.6-2mm and redshift distributions at 0.1-2mm, agree well. At the bright end (>1mJy), Shark galaxies are a mix of mergers and disk instabilities. These galaxies display a stacked FUV-to-FIR spectrum that agrees well with observations. We predict that current optical/NIR surveys are deep enough to detect bright (>1mJy) lambda=0.85-2mm-selected galaxies at z<5, but too shallow to detect counterparts at higher redshift. A James Webb Space Telescope 10,000s survey should detect all counterparts for galaxies with $S_{\rm 0.85mm}>0.01$mJy. We predict SMG's disks contribute significantly (negligibly) to the rest-frame UV (IR). We investigate the 01mJy lambda=0.85-2mm-selected galaxies finding their: (i) stellar masses are $>10^{10.2}M_{\odot}$, with the 2mm ones tracing the most massive galaxies ($>10^{11}M_{\odot}$); (ii) specific star formation rates (SFR) are mildly (~3-10x) above the main sequence (MS); (iii) host halo masses are $\gtrsim 10^{12.3}M_{\odot}$, with 2mm galaxies tracing the most massive halos (proto-clusters); (iv) SMGs have lower dust masses ($\approx 10^{8}M_{\odot}$), higher dust temperatures ($\approx 40-45$K) and higher rest-frame V-band attenuation (>1.5) than MS galaxies; (v) sizes decrease with redshift, from 4kpc at z=1 to <1kpc at z=4; (vi) the Carbon Monoxide line spectra of $S_{\rm 0.85mm}>1$mJy sources peak at 4->3. Finally, we study the contribution of SMGs to the molecular gas and cosmic SFR density at 01mJy sources make a negligible contribution at z>3 and z>5, respectively, suggesting current observations have unveiled the majority of the star formation at 0

We present the detection of multiple CO line transitions with ALMA in a few tens of infrared-selected galaxies on and above the main sequence at z=1.1-1.7. We reliably detected the emission of CO(5-4), CO(2-1), and CO(7-6)+[CI](2-1) in 50, 33, and 13 galaxies, respectively, and we complemented this information with available CO(4-3) and [CI](1-0) fluxes for part of the sample, and modeling of the optical-to-mm SEDs. We retrieve a quasi-linear relation between LIR and CO(5-4) or CO(7-6) for main-sequence galaxies and starbursts, corroborating the hypothesis that these transitions can be used as SFR tracers. We find the CO excitation to steadily increase as a function of the star formation efficiency, the mean intensity of the radiation field warming the dust, the surface density of SFR, and, less distinctly, with the distance from the main sequence. This adds to the tentative evidence for higher excitation of the CO+[CI] SLED of starbursts relative to that for main-sequence objects, where the dust opacities play a minor role in shaping the high-J CO transitions in our sample. However, the distinction between the average SLED of upper main-sequence and starburst galaxies is blurred, driven by a wide variety of intrinsic shapes. LVG radiative transfer modeling demonstrates the existence of a highly excited component that elevates the CO SLED of high-redshift main-sequence and starbursting galaxies above the typical values observed in the disk of the Milky Way. This excited component is dense and it encloses ~50% of the total molecular gas mass in main-sequence objects. We interpret the observed trends involving the CO excitation as mainly driven by a combination of large SFRs and compact sizes, naturally connected with enhanced dense molecular gas fractions and higher dust and gas temperatures, due to increasing UV radiation fields, cosmic ray rates, and dust/gas coupling. [Abridged], Comment: 18 pages, 9 figures + Appendix, accepted in A&A. The data compilation will be released as a .fits table on the journal's website. For an early access, please contact the corresponding author. The plots of spectra and far-infrared SEDs will be available on the journal's website and temporarily at this link: https://bit.ly/spectra_and_seds_valentino_2020c

We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K ($\beta=1.9$) derived under the common assumption of an optically thin far-infrared dust emission, but fully consistent with $T_{\rm d}=52\pm5$ K of a general opacity model where the optical depth ($\tau$) reaches unity at a wavelength of $\lambda_0=170\pm23$ $\mu$m. Moreover, the general opacity solution returns a factor of $\sim 2\times$ lower dust mass and, hence, a lower molecular gas mass for a fixed gas-to-dust ratio, than with the optically thin dust model. The derived properties of GN20 thus provide an appealing solution to the puzzling discovery of starbursts appearing colder than main-sequence galaxies above $z>2.5$, in addition to a lower dust-to-stellar mass ratio that approaches the physical value predicted for starburst galaxies., Comment: Comments: 6 pages, 4 figures, accepted for publication in Astronomy and Astrophysics Letters on February 4th 2020

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters., Comment: Accepted ApJS; MAST archive will be live concurrent with publication

We report ALMA observations of the neutral atomic carbon transitions [CI] and multiple CO lines in a sample of $\sim30$ main sequence galaxies at $z\sim1$, including novel information on [CI](2-1) and CO(7-6) for 7 of such normal objects. We complement our observations with a collection of $>200$ galaxies with coverage of similar transitions, spanning the $z=0-4$ redshift interval and a variety of ambient conditions from local to high-redshift starbursts. We find systematic variations in the [CI]/IR and [CI]/high-$J$ ($J=7$) CO luminosity ratios among the various samples. We interpret these differences as increased dense molecular gas fractions and star formation efficiencies in the strongest high-redshift starbursts with respect to normal main sequence galaxies. We further report constant $L'_{\rm [CI]2-1}$/$L'_{\rm [CI]1-0}$ ratios across the galaxy populations and redshifts, suggesting that gas temperatures $T_{\rm exc}$ traced by [CI] do not strongly vary. We find only a mild correlation with $T_{\rm dust}$ and that, generally, $T_{\rm exc} \lesssim T_{\rm dust}$. We fit the line ratios with classical PDR models, retrieving consistently larger densities and intensities of the UV radiation fields in submm galaxies than in main sequence and local objects. However, these simple models fall short in representing the complexity of a multiphase interstellar medium and should be treated with caution. Finally, we compare our observations with the Santa Cruz semi-analytical model of galaxy evolution, recently extended to simulate submm emission. While we confirm the success in reproducing the CO lines, we find systematically larger [CI] luminosities at fixed IR luminosity than predicted theoretically. This highlights the necessity of improving our understanding of the mechanisms regulating the [CI] emission on galactic scales. We release our data compilation to the community., Comment: 20 pages + Appendix, accepted for publication in the Astrophysical Journal on Dec. 20th, 2019. The data compilation is available in the online version of the article

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years., Comment: 30 pages, 10 figures, accepted in ApJ

We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$\mu$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $\tau_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr)., Comment: 24 pages, 8 figures, accepted for publication in ApJ

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

We present the first stellar velocity dispersion measurement of a massive quenching galaxy at z=4.01. The galaxy is first identified as a massive z>~4 galaxy with suppressed star formation from photometric redshifts based on deep multi-band data in the UKIDSS Ultra Deep Survey field. A follow-up spectroscopic observation with MOSFIRE on Keck revealed strong multiple absorption features, which are identified as Balmer absorption lines, giving a secure redshift of z=4.01. Thanks to the high S/N of the spectrum, we are able to estimate the stellar velocity dispersion, sigma=268+/-59 km/s. This velocity dispersion is consistent with that of massive galaxies today, implying no significant evolution in stellar velocity dispersion over the last 12 Gyr. Based on an upper limit on its physical size from deep optical images (r_eff<1.3 kpc), we find that its dynamical mass is consistent with the stellar mass inferred from photometry. Furthermore, the galaxy is located on the mass fundamental plane extrapolated from lower redshift galaxies. Combining all these results, we find that the velocity dispersion does not significantly evolve with redshift, although the size and mass of massive quenched galaxies do. This suggests that the mass in the core of massive galaxies does not evolve significantly, while most of the mass growth occurs in the outskirts of the galaxies, which also increases the size. This picture is consistent with a two-phase formation scenario in which mass and size growth is due to accretion in the outskirts of galaxies via mergers., Comment: Published in the Astrophysical Journal letters. Fixed an error in dynamical mass

We report two secure ($z=3.775, 4.012$) and one tentative ($z\approx3.767$) spectroscopic confirmations of massive and quiescent galaxies through $K$-band observations with Keck/MOSFIRE and VLT/X-Shooter. The stellar continuum emission, the absence of strong nebular emission lines and the lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of $\mathrm{log}(M_{\star}/M_\odot)\sim11$ and ongoing star formation rates placing these galaxies $\gtrsim 1-2$ dex below the main sequence at their redshifts. The adopted parametrization of the star formation history suggests that these sources experienced a strong ($\langle \rm SFR \rangle \sim 1200-3500\,M_\odot\,\mathrm{yr}^{-1}$) and short ($\sim 50$ Myr) burst of star formation, peaking $\sim 150-500$ Myr before the time of observation, all properties reminiscent of the characteristics of sub-millimeter galaxies (SMGs) at $z>4$. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest sub-mm surveys detecting not only the most extreme starbursts, but also more normal galaxies. We support these findings by further exploring the Illustris-TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at $z\sim3-4$ and SMGs at $z\sim4-5$, with number densities and properties in agreement with the observations at $z\sim3$, but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past and, similarly, not all bright SMGs quench by $z\sim3$, both fractions depending on the threshold assumed to define the SMGs themselves., Comment: 20 pages (+appendix), 11 figures, accepted for publication in the Astrophysical Journal (ApJ). This work comes with a companion letter by Tanaka, Valentino, Toft et al. 2019, ApJL, 885, L34