Search

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 study the interstellar medium (ISM) properties as a function of the molecular gas size for 77 infrared-selected galaxies at $z \sim 1.3$. Molecular gas sizes are measured on ALMA images that combine CO(2-1), CO(5-4) and underlying continuum observations, and include CO(4-3), CO(7-6)+[CI]($^3 P_2-^3P_1$), [CI]($^3 P_1-^3P_0$) observations for a subset of the sample. The $\gtrsim 46 \%$ of our galaxies have a compact molecular gas reservoir, and lie below the optical disks mass-size relation. Compact galaxies on and above the main sequence have higher CO excitation and star formation efficiency than galaxies with extended molecular gas reservoirs, as traced by CO(5-4)/CO(2-1) and CO(2-1)/$L_{\rm IR, SF}$ ratios. Average CO+[CI] spectral line energy distributions indicate higher excitation in compacts relative to extended sources. Using CO(2-1) and dust masses as molecular gas mass tracers, and conversion factors tailored to their ISM conditions, we measure lower gas fractions in compact main-sequence galaxies compared to extended sources. We suggest that the sub-millimetre compactness, defined as the ratio between the molecular gas and the stellar size, is an unavoidable information to be used with the main sequence offset to describe the ISM properties of galaxies, at least above $M_{\star} \geqslant 10^{10.6}$ M$_{\odot}$, where our observations fully probe the main sequence scatter. Our results are consistent with mergers driving the gas in the nuclear regions, enhancing the CO excitation and star formation efficiency. Compact main-sequence galaxies are consistent with being an early post-starburst population following a merger-driven starburst episode, stressing the important role of mergers in the evolution of massive galaxies., Accepted on MNRAS on October 11th 2021. Added Figure 10, 11, 12 following referee's comments. Abstract slightly modified to fit to arXiv's requirements

We conduct a systematic search for protocluster candidates at $z \geq 6$ in the COSMOS field using the recently released COSMOS2020 source catalog. We select galaxies using a number of selection criteria to obtain a sample of galaxies that have a high probability of being inside a given redshift bin. We then apply overdensity analysis to the bins using two density estimators, a Weighted Adaptive Kernel Estimator and a Weighted Voronoi Tessellation Estimator. We have found 15 significant ($>4\sigma$) candidate galaxy overdensities across the redshift range $6\le z\le7.7$. The majority of the galaxies appear to be on the galaxy main sequence at their respective epochs. We use multiple stellar-mass-to-halo-mass conversion methods to obtain a range of dark matter halo mass estimates for the overdensities in the range of $\sim10^{11-13}\,M_{\rm \odot}$, at the respective redshifts of the overdensities. The number and the masses of the halos associated with our protocluster candidates are consistent with what is expected from the area of a COSMOS-like survey in a standard $\Lambda$CDM cosmology. Through comparison with simulation, we expect that all the overdensities at $z\simeq6$ will evolve into a Virgo-/Coma-like clusters at present (i.e., with masses $\sim 10^{14}-10^{15}\,M_{\rm \odot}$). Compared to other overdensities identified at $z \geq 6$ via narrow-band selection techniques, the overdensities presented appear to have $\sim10\times$ higher stellar masses and star-formation rates. We compare the evolution in the total star-formation rate and stellar mass content of the protocluster candidates across the redshift range $6\le z\le7.7$ and find agreement with the total average star-formation rate from simulations., Comment: 52 pages, 32 figues, 18 tables, main text is 30 pages, appendix is 22 pages, to be published in ApJ

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

One of the most prominent features of galaxy clusters is the presence of a dominant population of massive ellipticals in their cores. Stellar archaeology suggests that these gigantic beasts assembled most of their stars in the early Universe via starbursts. However, the role of dense environments and their detailed physical mechanisms in triggering starburst activities remain unknown. Here we report spatially resolved Atacama Large Millimeter/submillimeter Array (ALMA) observations of the CO $J= 3-2$ emission line, with a resolution of about 2.5 kiloparsecs, toward a forming galaxy cluster core with starburst galaxies at $z=2.51$. In contrast to starburst galaxies in the field often associated with galaxy mergers or highly turbulent gaseous disks, our observations show that the two starbursts in the cluster exhibit dynamically cold (rotation-dominated) gas-rich disks. Their gas disks have extremely low velocity dispersion ($\sigma_{\mathrm{0}} \sim 20-30$ km s$^{-1}$), which is three times lower than their field counterparts at similar redshifts. The high gas fraction and suppressed velocity dispersion yield gravitationally unstable gas disks, which enables highly efficient star formation. The suppressed velocity dispersion, likely induced by the accretion of corotating and coplanar cold gas, might serve as an essential avenue to trigger starbursts in massive halos at high redshifts., Comment: 19 pages, 14 figures, 3 tables. Accepted for publication in A&A

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: A&A Letter in press

The MAMMOTH-1 nebula at $z=2.317$ is an enormous Ly$\alpha$ nebula (ELAN) extending to a $\sim$440 kpc scale at the center of the extreme galaxy overdensity BOSS 1441. In this paper, we present observations of the $\rm CO(3-2)$ and 250 GHz dust-continuum emission from the MAMMOTH-1 using the IRAM NOrthern Extended Millimeter Array. Our observations show that $\rm CO(3-2)$ emission in this ELAN has not extended widespread emission into the circum- and inter-galactic media. We also find a remarkable concentration of six massive galaxies in $\rm CO(3-2)$ emission in the central $\sim$100 kpc region of the ELAN. Their velocity dispersions suggest a total halo mass of $M_{200c} \sim 10^{13.1} M_{\odot}$, marking a possible protocluster core associated with the ELAN. The peak position of the $\rm CO(3-2)$ line emission from the obscured AGN is consistent with the location of the intensity peak of MAMMOTH-1 in the rest-frame UV band. Its luminosity line ratio between the $\rm CO(3-2)$ and $\rm CO(1-0)$ $r_{3,1}$ is 0.61$\pm$0.17. The other five galaxies have $\rm CO(3-2)$ luminosities in the range of (2.1-7.1)$\times 10^9$ K $\rm km\,s^{-1}$ pc$^2$, with the star-formation rates derived from the 250GHz continuum of ($, Comment: 14 pages, 5 figures, 2 tables, accepted for publication in ApJ

We study the radio spectral properties of 2,094 star-forming galaxies (SFGs) by combining our early science data from the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey with VLA, GMRT radio data, and rich ancillary data in the COSMOS field. These SFGs are selected at VLA 3GHz, and their flux densities from MeerKAT 1.3GHz and GMRT 325MHz imaging data are extracted using the "super-deblending" technique. The median radio spectral index is $\alpha_{\rm 1.3GHz}^{\rm 3GHz}=-0.80\pm0.01$ without significant variation across the rest-frame frequencies ~1.3-10GHz, indicating radio spectra dominated by synchrotron radiation. On average, the radio spectrum at observer-frame 1.3-3GHz slightly steepens with increasing stellar mass with a linear fitted slope of $\beta=-0.08\pm0.01$, which could be explained by age-related synchrotron losses. Due to the sensitivity of GMRT 325MHz data, we apply a further flux density cut at 3GHz ($S_{\rm 3GHz}\ge50\,\mu$Jy) and obtain a sample of 166 SFGs with measured flux densities at 325MHz, 1.3GHz, and 3GHz. On average, the radio spectrum of SFGs flattens at low frequency with the median spectral indices of $\alpha^{\rm 1.3GHz}_{\rm 325MHz}=-0.59^{+0.02}_{-0.03}$ and $\alpha^{\rm 3.0GHz}_{\rm 1.3GHz}=-0.74^{+0.01}_{-0.02}$. At low frequency, our stacking analyses show that the radio spectrum also slightly steepens with increasing stellar mass. By comparing the far-infrared-radio correlations of SFGs based on different radio spectral indices, we find that adopting $\alpha_{\rm 1.3GHz}^{\rm 3GHz}$ for $k$-corrections will significantly underestimate the infrared-to-radio luminosity ratio ($q_{\rm IR}$) for >17% of the SFGs with measured flux density at the three radio frequencies in our sample, because their radio spectra are significantly flatter at low frequency (0.33-1.3GHz)., Comment: 15 pages, 6 figures, 1 table; accepted for publication in MNRAS on 2021 Aug 4

Tracing the emergence of the massive quiescent galaxy (QG) population requires the build-up of reliable quenched samples. We present Hubble Space Telescope WFC3/G141 grism spectra of 10 quiescent galaxy candidates selected at $2.5, Comment: Accepted for publication in A&A. 28 pages, 16 figures + additional plots in appendix A

Early-type galaxies (ETGs) are of crucial importance to trace back the galaxy mass assembly across cosmic time, yet their formation and quenching remain remarkably elusive. The discoveries of massive, dead galaxies at ever-growing redshifts provided compelling evidence to push their formation up to redshift > 4–5 when the Universe was barely 1 Gyr old. In this talk I will present our results on the ages of a new sample of ETGs at z ∼ 3, built by exploiting HST WFC3/G141 rest-frame optical/near-UV grism spectroscopy to study the nature of 10 passive galaxy candidates at 2.5 < z < 3.5 in COSMOS.This work is part of a PhD project aimed at quantifying the parent space density of distant genuine passive galaxies. I will also discuss the importance of multi-wavelength data in clarifying the degree of contamination by dusty star-forming galaxies affecting the color selection.

We have discovered a 300kpc-wide giant Lya nebula centered on the massive galaxy group RO-1001 at z=2.91 in the COSMOS field. Keck Cosmic Web Imager observations reveal three cold gas filaments converging into the center of the potential well of its ~4x10^13Msun dark matter halo, hosting 1200Msun/yr of star formation as probed by ALMA and NOEMA observations. The nebula morphological and kinematics properties and the prevalence of blueshifted components in the Lya spectra are consistent with a scenario of gas accretion. The upper limits on AGN activity and overall energetics favor gravity as the primary Lya powering source and infall as the main source of gas flows to the system. Although interpretational difficulties remain, with outflows and likely also photoionization with ensuing recombination still playing a role, this finding provides arguably an ideal environment to quantitatively test models of cold gas accretion and galaxy feeding inside an actively star-forming massive halo at high redshift., Comment: 20 pages, 12 figures, 3 tables. Accepted for publication in A&A

(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

Feedback-driven winds from star formation or active galactic nuclei might be a relevant channel for the abrupt quenching star formation in massive galaxies. However, both observations and simulations support the idea that these processes are non-conflictingly co-evolving and self-regulating. Furthermore, evidence of disruptive events that are capable of fast quenching is rare, and constraints on their statistical prevalence are lacking. Here we present a massive starburst galaxy at z=1.4 which is ejecting $46 \pm 13$\% of its molecular gas mass at a startling rate of $\gtrsim 10,000$ M$_{\odot}{\rm yr}^{-1}$. A broad component that is red-shifted from the galaxy emission is detected in four (low- and high-J) CO and [CI] transitions and in the ionized phase, which ensures a robust estimate of the expelled gas mass. The implied statistics suggest that similar events are potentially a major star-formation quenching channel. However, our observations provide compelling evidence that this is not a feedback-driven wind, but rather material from a merger that has been probably tidally ejected. This finding challenges some literature studies in which the role of feedback-driven winds might be overstated., Published in Nature Astronomy as an Article on 11 January 2021. The published version is available at https://www.nature.com/articles/s41550-020-01268-x (DOI 10.1038/s41550-020-01268-x)

We analyse the far-infrared properties of $\sim$ 5,000 star-forming galaxies at $z, 39 pages, 21 figures, 3 tables. Accepted to ApJ July 10, 2021

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., 21 pages + Appendices, 9 figures, 6 tables. Accepted for publication in A&A

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

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

Based on optical/near-IR Magellan FIRE spectra of 25 starburst galaxies at 0.5 < z < 0.9, Calabr�� et al.(2018) showed that their attenuation properties can be explained by a single-parameter sequence of total obscurations ranging from A(V)=2 to A(V)=30 towards the starburst core centers in a mixed stars and dust configuration. We investigate here the origin of this sequence for the same sample. We show that total attenuations anti-correlate with the starburst sizes in radio (3 GHz) with a significance larger than 5sigma and a scatter of 0.26 dex. More obscured and compact starbursts also show enhanced N2 (=[NII]/Halpha) ratios and larger line velocity widths that we attribute to an increasing shock contribution toward later merger phases, driven by deeper gravitational potential wells at the coalescence. Additionally, the attenuation is also linked to the equivalent width (EW) of hydrogen recombination lines, which is sensitive to the luminosity weighted age of the relatively unobscured stellar populations. Overall, the correlations among A(V), radio size, line width, N2 and EW of Balmer/Paschen lines converge towards suggesting an evolutionary sequence of merger stages: all of these quantities are likely to be good time-tracers of the merger phenomenon, and their large spanned range appears to be characteristic of the different merger phases. Half of our sample at higher obscurations have radio sizes approximately 3 times smaller than early type galaxies at the same redshift, suggesting that, in analogy with local Ultraluminous Infrared galaxies (ULIRGs), these cores cannot be directly forming elliptical galaxies. Finally, we detect mid-IR AGN torus for half of our sample and additional X-ray emission for 6 starbursts; intriguingly, the latter have systematically more compact sizes, suggestive of emerging AGNs towards later merger stages, possibly precursors of a later QSO phase., 30 pages, 21 figures and 2 tables. Accepted for publication to A&A. The key result is in Fig. 3, while a Cartoon of the merger sequence is shown in Fig. 9

Context. The upcoming new generation of optical spectrographs on four-meter-class telescopes, with their huge multiplexing capabilities, excellent spectral resolution, and unprecedented wavelength coverage, will provide invaluable information for reconstructing the history of star formation in individual galaxies up to redshifts of about 0.7. Aims. We aim at defining simple but robust and meaningful physical parameters that can be used to trace the coexistence of widely diverse stellar components: younger stellar populations superimposed on the bulk of older ones. Methods. We produced spectra of galaxies closely mimicking data from the forthcoming Stellar Populations at intermediate redshifts Survey (StePS), a survey that uses the WEAVE spectrograph on the William Herschel Telescope. First, we assessed our ability to reliably measure both ultraviolet and optical spectral indices in galaxies of different spectral types for typically expected signal-to-noise ratios. We then analyzed such mock spectra with a Bayesian approach, deriving the probability density function of r- and u-band light-weighted ages as well as of their difference. Results. We find that the ultraviolet indices significantly narrow the uncertainties in estimating the r- and u-band light-weighted ages and their difference in individual galaxies. These diagnostics, robustly retrievable for large galaxy samples even when observed at moderate signal-to-noise ratios, allow us to identify secondary episodes of star formation up to an age of ∼0.1 Gyr for stellar populations older than ∼1.5 Gyr, pushing up to an age of ∼1 Gyr for stellar populations older than ∼5 Gyr. Conclusions. The difference between r-band and u-band light-weighted ages is shown to be a powerful diagnostic to characterize and constrain extended star-formation histories and the presence of young stellar populations on top of older ones. This parameter can be used to explore the interplay between different galaxy star-formation histories and physical parameters such as galaxy mass, size, morphology, and environment.© ESO 2019., L. C. wish to acknowledge financial support from Premiale 2015 MITiC, program 1.05.06.10 and Comunidad de Madrid under Atraccion de Talento grant 2018-T2/TIC-11612. A. I., S. Z., M. L., A. G., A. M., and C. T. acknowledge the financial support from the INAF PRIN-SKA 2017 program 1.05.01.88.04.ESKAPE-HI. I. L. acknowledges financial support from INAF-WEAVE funds, program 1.05.03.04.05. R. G. B. acknowledges support from the Spanish Ministry of Economy and Competitiveness through grant 205 AYA2016-77846-P and the State Agency for Research of the Spanish MCIU through the >Center of Excellence Severo Ochoa> award to 2017 the Instituto de Astrofisica de Andalucia (SEV-2017-0709). A. F. M. has received financial support through the Postdoctoral Junior Leader Fellowship Programme from >La Caixa> Banking Foundation (LCF/BQ/LI18/11630007). B. V. acknowledges the financial support from INAF Main Stream 2018 (P. I.: B. Vulcani).

We infer and compare the specific X-ray luminosity distributions for a sample of massive (i.e. $\log_{10} (M*/M\odot) > 10.5$) galaxies split according to their far-infrared-derived star-forming properties (i.e., starburst and non-starburst) and redshift. We model each distribution as a power-law with an upper and lower turnover, and adopt a maximum likelihood method to include information from non-detections in the form of upper limits. When we use our inferred distributions to calculate the ratios of high to low sLx AGN (corresponding to above and below $0.1\lambda_{\text{Edd}}$, respectively) we find that starbursts have significantly higher proportions of high sLx AGN compared to their non-starburst counterparts. These findings help explain the increase in average X-ray luminosity in bins of increasing SFR reported by previous studies., Comment: 14 pages, 10 figures, accepted for publication in MNRAS

While the formation of stellar clumps in distant galaxies is usually attributed to gravitational violent disk instabilities, we show here that major mergers also represent a competitive mechanism to form bright clumps. Using ~0.1'' resolution ACS F814W images in the entire COSMOS field, we measure the fraction of clumpy emission in 109 main sequence (MS) and 79 Herschel-detected starbursts (off-MS) galaxies at 0.5 < z < 0.9, representative of normal versus merger induced star-forming activity, respectively. We additionally identify merger samples from visual inspection and from Gini-M20 morphological parameters. Regardless of the merger criteria adopted, the clumpiness distribution of merging systems is different from that of normal isolated disks at > 99.5 % confidence level, with the former reaching higher clumpiness values, up to 20 % of the total galaxy emission. We confirm the merger induced clumpiness enhancement with novel hydrodynamical simulations of colliding galaxies with gas fractions typical of z ~ 0.7. Multi-wavelength images of three starbursts in the CANDELS field support the young nature of clumps, which are likely merger products rather than older pre-existing structures. Finally, for a subset of 19 starbursts with existing near-IR rest frame spectroscopy, we find that the clumpiness is mildly anti-correlated with the merger phase, decreasing towards final coalescence. Our result can explain recent ALMA detections of clumps in hyperluminous high-z starbursts, while normal objects are smooth. This work raises a question on the role of mergers on the origin of clumps in high redshift galaxies in general., 24 pages, 19 figures and 1 table ; Accepted for publication in Astronomy & Astrophysics

We present an enhanced version of the multiwavelength spectral modeling code MAGPHYS that allows the estimation of galaxy photometric redshift and physical properties (e.g., stellar mass, star formation rate, dust attenuation) simultaneously, together with robust characterization of their uncertainties. The self-consistent modeling over ultraviolet to radio wavelengths in MAGPHYS+photo-z is unique compared to standard photometric redshift codes. The broader wavelength consideration is particularly useful for breaking certain degeneracies in color vs. redshift for dusty galaxies with limited observer-frame ultraviolet and optical data (or upper limits). We demonstrate the success of the code in estimating redshifts and physical properties for over 4,000 infrared-detected galaxies at 0.4, Comment: 27 pages, 14 figures, 3 tables. Accepted for publication in ApJ. The MAGPHYS+photo-z code can be downloaded at http://www.iap.fr/magphys/download.html

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 $��_{\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$., 7 pages, 5 figures. Accepted for publication in ApJ Letters

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., Accepted for publication in A&A

We study the Carbon Monoxide (CO) excitation, mean molecular gas density and interstellar radiation field (ISRF) intensity in a comprehensive sample of 76 galaxies from local to high redshift (z~0-6), selected based on detections of their CO transitions J=2-1 and 5-4 and their optical/infrared/(sub-)millimeter spectral energy distributions (SEDs). We confirm the existence of a tight correlation between CO excitation as traced by the CO(5-4)/(2-1) line ratio (R52), and the mean ISRF intensity U as derived from infrared SED fitting using dust SED templates. By modeling the molecular gas density probability distribution function (PDF) in galaxies and predicting CO line ratios with large velocity gradient radiative transfer calculations, we present a framework linking global CO line ratios to the mean molecular hydrogen gas density nH2 and kinetic temperature Tkin. Mapping in this way observed R52 ratios to nH2 and Tkin probability distributions, we obtain positive U-nH2 and U-Tkin correlations, which imply a scenario in which the ISRF in galaxies is mainly regulated by Tkin and (non-linearly) by nH2. A small fraction of starburst galaxies showing enhanced nH2 could be due to merger-driven compaction. Our work demonstrates that ISRF and CO excitation are tightly coupled, and that density-PDF modeling is a promising tool for probing detailed ISM properties inside galaxies., Comment: Accepted for publication in ApJ

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., 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 an analysis of the molecular gas properties, based on CO(2 - 1) emission, of twelve starburst galaxies at z~1.6 selected by having a boost (>~4x) in their star formation rate (SFR) above the average star-forming galaxy at an equivalent stellar mass. ALMA observations are acquired of six additional galaxies than previously reported through our effort. As a result of the larger statistical sample, we significantly detect, for the first time at high-z, a systematically lower L'_CO/L_IR ratio in galaxies lying above the star-forming `main sequence' (MS). Based on an estimate of alpha_CO (i.e., the ratio of molecular gas mass to L'_CO(1-0)), we convert the observational quantities (e.g., L'_CO/L_IR) to physical units (M_gas/SFR) that represent the gas depletion time or its inverse, the star formation efficiency. We interpret the results as indicative of the star formation efficiency increasing in a continuous fashion from the MS to the starburst regime, whereas the gas fractions remain comparable to those of MS galaxies. Although, the balance between an increase in star-formation efficiency or gas fraction depends on the adopted value of alpha_CO as discussed., 16 pages, 8 figures

We explore the distribution of RMS=SFR/SFR_MS (where SFR_MS is the star formation rate of "Main Sequence" star-forming galaxies) for AGN hosts at z=1. We split our sample into two bins of X-ray luminosity divided at Lx=2x10^43erg s-1 to investigate whether the RMS distribution changes as a function of AGN power. Our main results suggest that, when the RMS distribution of AGN hosts is modelled as a log-normal distribution (i.e. the same shape as that of MS galaxies), galaxies hosting more powerful X-ray AGNs (i.e. Lx>2x10^43erg s-1) display a narrower RMS distribution that is shifted to higher values compared to their lower Lx counterparts. In addition, we find that more powerful X-ray AGNs have SFRs that are more consistent with that of MS galaxies compared to lower Lx AGNs. Despite this, the mean SFRs (as opposed to RMS) measured from these distributions are consistent with the previously observed flat relationship between SFR and Lx. Our results suggest that the typical star-forming properties of AGN hosts change with Lx , and that more powerful AGNs typically reside in more MS-like star-forming galaxies compared to lower Lx AGNs., 6 pages, 4 figures, Accepted for publication in MNRAS Letter

We present a new technique to measure multi-wavelength "Super-deblended" photometry from highly confused images, which we apply to Herschel and ground-based far-infrared (FIR) and (sub-)millimeter (mm) data in the northern field of the Great Observatories Origins Deep Survey (GOODS). There are two key novelties. First, starting with a large database of deep Spitzer 24{\mu}m and VLA 20cm detections that are used to define prior positions for fitting the FIR/submm data, we perform an active selection of useful priors independently at each frequency band, moving from less to more confused bands. Exploiting knowledge of redshift and all available photometry, we identify hopelessly faint priors that we remove from the fitting pool. This approach significantly reduces blending degeneracies and allows reliable photometry to be obtained for galaxies in FIR+mm bands. Second, we obtain well-behaved, nearly Gaussian flux density uncertainties, individually tailored to all fitted priors in each band. This is done by exploiting extensive simulations that allow us to calibrate the conversion of formal fitting uncertainties to realistic uncertainties depending on quantities directly measurable. We achieve deeper detection limits with high fidelity measurements and uncertainties at FIR+mm bands. As an illustration of the utility of these measurements, we identify 70 galaxies with z>3 and reliable FIR+mm detections. We present new constraints on the cosmic star formation rate density at 33 dusty galaxies that are missed by optical-to-near-infrared color selection. Photometric measurements for 3306 priors, including over 1000 FIR+mm detections are released publicly with our catalog., Comment: Published in The Astrophysical Journal, Volume 853, Number 2. The entire machine-readable catalog is available from ApJ at http://iopscience.iop.org/0004-637X/853/2/172/suppdata/apjaaa600t2_mrt.txt . We also provide the same catalog in FITS format at https://goo.gl/aiJVs9

Early type galaxies (ETG) contain most of the stars present in the local Universe and, above a stellar mass of ~5e10 Msun, vastly outnumber spiral galaxies like the Milky Way. These massive spheroidal galaxies have, in the present day, very little gas or dust, and their stellar populations have been evolving passively for over 10 billion years. The physical mechanisms that led to the termination of star formation in these galaxies and depletion of their interstellar medium remain largely conjectural. In particular, there are currently no direct measurements of the amount of residual gas that might be still present in newly quiescent spheroids at high redshift. Here we show that quiescent ETGs at z~1.8, close to their epoch of quenching, contained 2-3 orders of magnitude more dust at fixed stellar mass than local ETGs. This implies the presence of substantial amounts of gas (5-10%), which was however consumed less efficiently than in more active galaxies, probably due to their spheroidal morphology, and consistently with our simulations. This lower star formation efficiency, and an extended hot gas halo possibly maintained by persistent feedback from an active galactic nucleus (AGN), combine to keep ETGs mostly passive throughout cosmic time., Comment: 16 pages, 4 main figures main, 7 supplementary figures, 2 supplementary tables; accepted for publication in Nature Astronomy

We present a "super-deblended" far-infrared to (sub)millimeter photometric catalog in the Cosmic Evolution Survey (COSMOS), prepared with the method recently developed by Liu et al. 2018, with key adaptations. We obtain point spread function (PSF) fitting photometry at fixed prior positions including 88,008 galaxies detected in either VLA 1.4~GHz, 3~GHz and/or MIPS 24~$��$m images. By adding a specifically carved mass-selected sample (with an evolving stellar mass limit), a highly complete prior sample of 194,428 galaxies is achieved for deblending FIR/(sub)mm images. We performed ``active' removal of non relevant priors at FIR/(sub)mm bands using spectral energy distribution (SED) fitting and redshift information. In order to cope with the shallower COSMOS data we subtract from the maps the flux of faint non-fitted priors and explicitly account for the uncertainty of this step. The resulting photometry (including data from Spitzer, Herschel, SCUBA2, AzTEC, MAMBO and NSF's Karl G. Jansky Very Large Array at 3~GHz and 1.4~GHz) displays well behaved quasi-Gaussian uncertainties, calibrated from Monte Carlo simulations and tailored to observables (crowding, residual maps). Comparison to ALMA photometry for hundreds of sources provide a remarkable validation of the technique. We detect 11,220 galaxies over the 100--1200~$��$m range, extending to $z_{\rm phot}\sim7$. We conservatively selected a sample of 85 $z>4$ high redshift candidates, significantly detected in the FIR/(sub)mm, often with secure radio and/or Spitzer/IRAC counterparts. This provides a chance to investigate the first generation of vigorous starburst galaxies (SFRs$\sim1000\mathrm{M}_\odot$~yr$^{-1}$). The photometric and value added catalogs are publicly released., 49 pages (+appendix), 48 figures, 3 tables, accepted for publication in ApJ

We obtained optical/near-IR rest-frame Magellan FIRE spectra (including Pa$\beta$ and Pa$\gamma$) of 25 starburst galaxies at 0.5, Comment: ApJ Letters in press; the key result is in Figure 4 (left)

We show that the most distant X-ray detected cluster known to date, ClJ1001 at z=2.506, hosts a strong overdensity of radio sources. Six of them are individually detected (within 10") in deep 0.75" resolution VLA 3GHz imaging, with S(3GHz)>8uJy. Of the six, AGN likely affects the radio emission in two galaxies while star formation is the dominant source powering the remaining four. We searched for cluster candidates over the full COSMOS 2-square degree field using radio-detected 3GHz sources and looking for peaks in Sigma5 density maps. ClJ1001 is the strongest overdensity by far with >10sigma, with a simple z_phot>1.5 preselection. A cruder photometric rejection of z2.5. Samples of hundreds such high-redshift clusters could potentially constrain cosmological parameters and test cluster and galaxy formation models., ApJ Letters in press

We present rest-frame optical spectra from the FMOS-COSMOS survey of twelve $z \sim 1.6$ \textit{Herschel} starburst galaxies, with Star Formation Rate (SFR) elevated by $\times$8, on average, above the star-forming Main Sequence (MS). Comparing the H$\alpha$ to IR luminosity ratio and the Balmer Decrement we find that the optically-thin regions of the sources contain on average only $\sim 10$ percent of the total SFR whereas $\sim90$ percent comes from an extremely obscured component which is revealed only by far-IR observations and is optically-thick even in H$\alpha$. We measure the [NII]$_{6583}$/H$\alpha$ ratio, suggesting that the less obscured regions have a metal content similar to that of the MS population at the same stellar masses and redshifts. However, our objects appear to be metal-rich outliers from the metallicity-SFR anticorrelation observed at fixed stellar mass for the MS population. The [SII]$_{6732}$/[SII]$_{6717}$ ratio from the average spectrum indicates an electron density $n_{\rm e} \sim 1,100\ \mathrm{cm}^{-3}$, larger than what estimated for MS galaxies but only at the 1.5$\sigma$ level. Our results provide supporting evidence that high-$z$ MS outliers are the analogous of local ULIRGs, and are consistent with a major merger origin for the starburst event., Comment: 6 pages, 4 figures, Accepted for publication in ApJ Letters

We have analyzed the distributions in the color-magnitude diagram (CMD) of a large sample of face-on galaxies to minimize the effect of dust extinctions on galaxy color. About 300 thousand galaxies with $log(a/b) < $ 0.2 and redshift $z < 0.2$ are selected from the SDSS DR7 catalog. Two methods are employed to investigate the distributions of galaxies in the CMD including 1-D Gaussian fitting to the distributions in individual magnitude bins and 2-D Gaussian mixture model (GMM) fitting to galaxies as a whole. We find that in the 1-D fitting only two Gaussians are not enough to fit galaxies with the excess present between the blue cloud and the red sequence. The fitting to this excess defines the centre of the green-valley in the local universe to be $(u-r)_{0.1} = -0.121M_{r,0.1}-0.061$. The fraction of blue cloud and red sequence galaxies turns over around $M_{r,0.1} \sim -20.1$ mag, corresponding to stellar mass of $3\times10^{10}M_\odot$. For the 2-D GMM fitting, a total of four Gaussians are required, one for the blue cloud, one for the red sequence and the additional two for the green valley. The fact that two Gaussians are needed to describe the distributions of galaxies in the green valley is consistent with some models that argue for two different evolutionary paths from the blue cloud to the red sequence.

We present a sample of 34 spectroscopically confirmed BzK-selected ~1e11 Msun quiescent galaxies (pBzK) in the COSMOS field. The targets were initially observed with VIMOS on the VLT to facilitate the calibration of the photometric redshifts of massive galaxies at z >~ 1.5. Here we describe the reduction and analysis of the data, and the spectrophotometric properties of these pBzK galaxies. In particular, using a spatially resolved median 2D spectrum, we find that the fraction of stellar populations with ages, 12 pages, 10 figures, accepted for publication in Astronomy & Astrophysics