Search

We present the first spectroscopic confirmation of a dust-obscured dwarf galaxy, CEERS-14821. The analysis is performed combining JWST NIRCam broad-band photometry and NIRSpec/PRISM spectroscopic data. From the detection of multiple rest-frame optical lines, we derive that CEERS-14821 is located at $z=4.883\pm0.003$. Moreover, from a secure detection of the $H_{\alpha}$ and $H_{\beta}$ we derived that the galaxy has a dust extinction ranging from Av=2.2 to Av=3.3, depending on the assumed reddening law. This value is extremely large given that we estimated a low stellar mass around log(M/Mo)=8.0-8.2. Moreover, using different metallicity tracers, we verify that the galaxy is also metal-rich, with 12+log(O/H)>8.3. This is well above the expectation from both the mass-metallicity relation and the fundamental mass-metalliticy relation. CEERS-14821 is going through a burst of star formation, there are no indications of a strong contribution from an active galactic nuclei (f(AGN)<0.5 with respect to the total dust luminosity). Based on the rest-frame optical images, this source has a size compatible with galaxies of similar stellar mass and redshift. Finally, with the current data, it seems that there are galaxies closely interacting with CEERS-14821., Comment: 14 pages, 11 figures. Submitted to A&A

The MIGHTEE survey utilises the South African MeerKAT radio telescope to observe four extragalactic deep fields, with the aim of advancing our understanding of the formation and evolution of galaxies across cosmic time. MIGHTEE's frequency coverage encompasses the $\textrm{H}\scriptstyle\mathrm{I}$ line to a redshift of z $\simeq$ 0.58, and OH megamasers to z $\simeq$ 0.9. We present the MIGHTEE-$\textrm{H}\scriptstyle\mathrm{I}$ imaging products for the COSMOS field, using a total of 94.2 h on-target and a close-packed mosaic of 15 individual pointings. The spectral imaging covers two broad, relatively interference-free regions (960-1150 and 1290-1520~MHz) within MeerKAT's L-band, with up to 26 kHz spectral resolution (5.5 km s$^{-1}$ at $z$ = 0). The median noise in the highest spectral resolution data is 74 $\mu$Jy beam$^{-1}$, corresponding to a 5$\sigma$ $\textrm{H}\scriptstyle\mathrm{I}$ mass limit of 10$^{8.5}$ M$_{\odot}$ for a 300 km s$^{-1}$ line at $z$ = 0.07. The mosaics cover $>$4 deg$^{2}$, provided at multiple angular resolution / sensitivity pairings, with an angular resolution for $\textrm{H}\scriptstyle\mathrm{I}$ at $z$ = 0 of 12$''$. We describe the spectral line processing workflow that will be the basis for future MIGHTEE-$\textrm{H}\scriptstyle\mathrm{I}$ products, and validation of, and some early results from, the spectral imaging of the COSMOS field. We find no evidence for line emission at the position of the $z$ = 0.376 \HI~line reported from the CHILES survey at a $>$94 per cent confidence level, placing a 3$\sigma$ upper limit of 8.1 $\times$ 10$^{9}$ M$_{\odot}$ on $M_{\mathrm{HI}}$ for this galaxy. A public data release accompanies this article., Comment: 22 pages, 19 figures

Euclid will cover over 14000 $deg^{2}$ with two optical and near-infrared spectro-photometric instruments, and is expected to detect around ten million active galactic nuclei (AGN). This unique data set will make a considerable impact on our understanding of galaxy evolution and AGN. In this work we identify the best colour selection criteria for AGN, based only on Euclid photometry or including ancillary photometric observations, such as the data that will be available with the Rubin legacy survey of space and time (LSST) and observations already available from Spitzer/IRAC. The analysis is performed for unobscured AGN, obscured AGN, and composite (AGN and star-forming) objects. We make use of the spectro-photometric realisations of infrared-selected targets at all-z (SPRITZ) to create mock catalogues mimicking both the Euclid Wide Survey (EWS) and the Euclid Deep Survey (EDS). Using these catalogues we estimate the best colour selection, maximising the harmonic mean (F1) of completeness and purity. The selection of unobscured AGN in both Euclid surveys is possible with Euclid photometry alone with F1=0.22-0.23, which can increase to F1=0.43-0.38 if we limit at z>0.7. Such selection is improved once the Rubin/LSST filters (a combination of the u, g, r, or z filters) are considered, reaching F1=0.84 and 0.86 for the EDS and EWS, respectively. The combination of a Euclid colour with the [3.6]-[4.5] colour, which is possible only in the EDS, results in an F1-score of 0.59, improving the results using only Euclid filters, but worse than the selection combining Euclid and LSST. The selection of composite ($f_{{\rm AGN}}$=0.05-0.65 at 8-40 $\mu m$) and obscured AGN is challenging, with F1<0.3 even when including ancillary data. This is driven by the similarities between the broad-band spectral energy distribution of these AGN and star-forming galaxies in the wavelength range 0.3-5 $\mu m$., Comment: 25 pages, 28 figures, accepted for publication on A&A

In the past, high-z AGNs were given a minor role as possible drivers of reionization, despite initial evidences in favor of their large space densities at low luminosities by Chandra and HST. Recent observations from JWST are finding relatively large numbers of faint AGNs at z>4, convincingly confirming these early results. We present a sample of z~5 AGNs (both from wide, shallow ground-based surveys and from deep, pencil-beam observations from JWST), allowing to estimate their space densities with unprecedented accuracy. The bright end (M1450<-26) of the z~5 AGN luminosity function is well constrained, with a rather steep slope. The faint end (M1450>-22) indicates a high space density, the scatter is significant and the knee (M1450~-24) is mostly undetermined. Comparisons with state-of-the-art models find reasonable agreement with the observed AGN luminosity function at z=5, while the predicted space density evolution at higher redshifts appears to be too fast with respect to observational constraints. Given the large variance at the faint end, we consider different options in fitting the luminosity functions and deriving the ionizing emissivity. Even in the most conservative scenario, the photo-ionization rate produced by z~5 AGNs is consistent with the UV background measurements. A slow evolution of the space density of faint AGNs is observed, indicating that active SMBHs are probably producing large amounts of ionizing photons at z>6, well into the epoch of reionization. This is an important indication that high-z AGNs could be the major contributors to the reionization of the Universe., Comment: Accepted for publication on ApJ, 20 pages, 9 figures, 3 tables

Euclid will collect an enormous amount of data during the mission's lifetime, observing billions of galaxies in the extragalactic sky. Along with traditional template-fitting methods, numerous machine learning algorithms have been presented for computing their photometric redshifts and physical parameters (PPs), requiring significantly less computing effort while producing equivalent performance measures. However, their performance is limited by the quality and amount of input information, to the point where the recovery of some well-established physical relationships between parameters might not be guaranteed. To forecast the reliability of Euclid photo-$z$s and PPs calculations, we produced two mock catalogs simulating Euclid photometry. We simulated the Euclid Wide Survey (EWS) and Euclid Deep Fields (EDF). We tested the performance of a template-fitting algorithm (Phosphoros) and four ML methods in recovering photo-$z$s, PPs (stellar masses and star formation rates), and the SFMS. To mimic the Euclid processing as closely as possible, the models were trained with Phosphoros-recovered labels. For the EWS, we found that the best results are achieved with a mixed labels approach, training the models with wide survey features and labels from the Phosphoros results on deeper photometry, that is, with the best possible set of labels for a given photometry. This imposes a prior, helping the models to better discern cases in degenerate regions of feature space, that is, when galaxies have similar magnitudes and colors but different redshifts and PPs, with performance metrics even better than those found with Phosphoros. We found no more than 3% performance degradation using a COSMOS-like reference sample or removing u band data, which will not be available until after data release DR1. The best results are obtained for the EDF, with appropriate recovery of photo-$z$, PPs, and the SFMS., Comment: 26 pages, 13 figures. Accepted for publication on A&A

This paper presents a search for high redshift galaxies from the Euclid Early Release Observations program "Magnifying Lens." The 1.5 deg$^2$ area covered by the twin Abell lensing cluster fields is comparable in size to the few other deep near-infrared surveys such as COSMOS, and so provides an opportunity to significantly increase known samples of rare UV-bright galaxies at $z\approx6-8$ ($M_{\rm UV}\lesssim-22$). Beyond their still uncertain role in reionisation, these UV-bright galaxies are ideal laboratories from which to study galaxy formation and constrain the bright-end of the UV luminosity function. Of the 501994 sources detected from a combined $Y_{\rm E}$, $J_{\rm E}$, and $H_{\rm E}$ NISP detection image, 168 do not have any appreciable VIS/$I_{\rm E}$ flux. These objects span a range in spectral colours, separated into two classes: 139 extremely red sources; and 29 Lyman-break galaxy candidates. Best-fit redshifts and spectral templates suggest the former is composed of both $z\gtrsim5$ dusty star-forming galaxies and $z\approx1-3$ quiescent systems. The latter is composed of more homogeneous Lyman break galaxies at $z\approx6-8$. In both cases, contamination by L- and T-type dwarfs cannot be ruled out with Euclid images alone. Additional contamination from instrumental persistence is investigated using a novel time series analysis. This work lays the foundation for future searches within the Euclid Deep Fields, where thousands more $z\gtrsim6$ Lyman break systems and extremely red sources will be identified., Comment: 22 pages, 13 figures, paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

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

Near- to mid-Infrared observations (from Spitzer and JWST) have revealed a hidden population of galaxies at redshift z=3-6, called optically-dark objects, which are believed to be massive and dusty star-formers. While optically-dark sources are widely recognized as a significant component of the stellar mass function, the history of their stellar mass assembly remains unexplored. However, they are thought to be the progenitors of the more massive early-type galaxies found in present-day groups and clusters. It is thus important to examine the possible connection between dark sources and merging events, in order to understand the environment in which they live. Here, we report our search for close companions in a sample of 19 optically-dark objects identified in the SMACS0723 JWST deep field. They were selected in the NIRCam F444W band and undetected below 2mu. We restrict our analysis to the reddest (i.e. F277W-F444W> 1.3) and brightest (F444W< 26 mag) objects. We have identified an optically-dark source showing a very close companion (<0.5"). The spatially resolved SED fitting procedure indicates that all components lying within 1.5" from the dark source are indeed at z~5.7. Tidal features (leading to a whale shaped morphology) corroborate the hypothesis that the dark source is the most massive (log(M/Msun)>10.3) and dusty (Av~3 at the core) system of an ongoing merger with a mass ratio of ~10. Similar merging systems are identified in the SERRA simulations, allowing us to reconstruct their stellar mass assembly history and predict their molecular gas properties The discovery of mergers within dark galaxies at the end of the Epoch of Reionization underscores the importance of conducting a statistical search for additional candidates in deep NIRCam fields. Such research will aid in understanding the role of merging processes during the obscured phase of stellar mass accumulation., Comment: Accepted for publication in Astronomy & Astrophysics

Stellar bars are common morphological structures in the local Universe; according to optical and NIR surveys, they are present in about two-thirds of disc galaxies. These elongated structures are also believed to play a crucial role in secular evolutionary processes, because they are able to efficiently redistribute gas, stars, and angular momentum within their hosts, although it remains unclear as to whether they enhance or suppress star formation. A useful tool to investigate this ambiguity is the main sequence (MS) relation, which tightly links stellar mass ($M_{\star}$) and star formation rate (SFR). The main goal of this work is to explore star-formation processes in barred galaxies in order to assess the relevance of bars in star-formation quenching and whether or not they affect the typical log-linear trend of the resolved MS. To this purpose, we carried out a spatially resolved analysis on subkiloparsec (subkpc) scales for a sample of six nearby barred galaxies. We collected multi-wavelength photometric data from far-ultraviolet (FUV) to far-infrared (FIR) from the DustPedia database and applied a panchromatic spectral energy distribution (SED) fitting procedure on square apertures of fixed angular size (8" $\times$ 8") using the magphys code. For each galaxy, we obtain the distributions of stellar mass and SFR surface density and relate them in the $\log \Sigma_{\star}$ - $\log \Sigma_{\rm SFR}$ plane, deriving the spatially resolved MS relation. Although significant galaxy-to-galaxy variations are in place, we infer the presence of a common anti-correlation track in correspondence with the bar-hosting region, which shows systematically lower SFRs. This central quiescent signature can be interpreted as the result of a bar-driven depletion of gas reservoirs and a consequent halting of star formation. Our findings appear to support an inside-out quenching scenario., Comment: 15 pages, 6 figures, corrected typos, published in A&A

Recent discoveries of copious amounts of dust in quiescent galaxies (QGs) at high redshifts ($z\gtrsim 1-2$) challenge the conventional view that these objects have poor interstellar medium (ISM) in proportion to their stellar mass. We use the SIMBA cosmological simulation to explore the evolution of dust and cold gas content in QGs in relation to the quenching processes affecting them. We track the changes in the ISM dust abundance across the evolutionary history of QGs identified at $0 \lesssim z \lesssim2$ in the field and cluster environments. The QGs quench via diverse pathways, both rapid and slow, and exhibit a wide range of times elapsed between the quenching event and cold gas removal (from $\sim650$ Myr to $\sim8$ Gyr). We find that quenching modes attributed to the feedback from active galactic nuclei (AGN) do not affect dust and cold gas within the same timescales. Remarkably, QGs may replenish their dust content in the quenched phase primarily due to internal processes and marginally by external factors such as minor mergers. The key mechanism for re-formation of dust is prolonged grain growth on gas-phase metals, it is effective within $\sim100$ Myr after the quenching event, and rapidly increases the dust-to-gas mass ratio in QGs above the standard values ($\delta_{\rm DGR}\gtrsim1/100$). As a result, despite heavily depleted cold gas reservoirs, roughly half of QGs maintain little evolution in their ISM dust with stellar age within the first 2 Gyr following the quenching. Overall, we predict that relatively dusty QGs ($M_{\rm dust}/M_{\star}\gtrsim10^{-3}-10^{-4}$) arise from both fast and slow quenchers, and are prevalent in systems of intermediate and low stellar masses ($9<\log(M_{\star}/M_{\odot})<10.5$). This prediction poses an immediate quest for observational synergy between e.g., James Webb Space Telescope (JWST) and the Atacama Large Millimeter Array (ALMA)., Comment: 21 pages, 10 figures, 1 table, submitted to A&A

We present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. These sources are likely heavily-reddened AGN that trace a previously-hidden phase of dust-obscured black hole growth in the early Universe. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifting bandpasses to sample the same rest-frame emission blueward and redward of the Balmer break. This approach allows us to identify LRDs over a wider redshift range and is less susceptible to contamination from galaxies with strong breaks that otherwise lack a rising red continuum. The redshift distribution of our sample increases at $z<8$ and then undergoes a rapid decline at $z\sim4.5$, which may tie the emergence, and obscuration, of these sources to the inside-out growth that galaxies experience during this epoch. We find that LRDs are 2-3 dex more numerous than bright quasars at $z\sim5-7$, but their number density is only 0.6-1 dex higher than X-ray and UV selected AGN at these redshifts. Within our sample, we have identified the first X-ray detected LRDs at $z=3.1$ and $z=4.66$. An X-ray spectral analysis confirms that these AGN are moderately obscured with $\log\,(N_{\rm H}/{\rm cm}^{2}$) of $23.3^{+0.4}_{-1.3}$ and $22.72^{+0.13}_{-0.16}$. Our analysis reveals that reddened AGN emission dominates their rest-optical light, while the rest-UV originates from their host galaxies. We also present NIRSpec follow-up spectroscopy of 17 LRDs that show broad emission lines consistent with AGN activity. The confirmed AGN fraction of our sample is $71\%$ for sources with F444W$<26.5$. In addition, we find three LRDs with narrow blue-shifted Balmer absorption features in their spectra, suggesting an outflow of high-density, low ionization gas from near the central engine of these faint, red AGN., Comment: 23 pages, 17 figures, submitted to ApJ

The main sequence (MS) of star-forming galaxies (SFGs) is the tight relation between the galaxy stellar mass and its star formation rate (SFR) and was observed up to z ~ 6. The MS relation can be used as a reference for understanding the differences among galaxies, characterised by different rates of stellar production (starbursts, SFGs, and passive galaxies), and those inside a galaxy made up of different components (bulge, disk, and halo). To investigate peculiar features found in our sample galaxies, we focus here on their star formation history (SFH). We performed a spectral energy distribution fitting procedure that accounted for the energetic balance between UV and far-IR radiation on a sample of eight nearby face-on spiral galaxies from the DustPedia sample. This approach allowed us to study the spatially resolved MS of the sample and to recover the past SFH. By exploiting the BAGPIPES code, we constrained the SFHs for each galaxy with a delayed exponentially declining model to derive their mass-weighted age (tMW). A central old region (tMW up to~7Gyr, consistent with the presence of a bulge for various systems) is followed by younger regions in which the disks are still forming stars (tMW~4Gyr). At larger distances, tMW increases mildly in general. Strikingly, in two galaxies (NGC4321 and NGC5194), we found a steep increase in tMW that reached levels similar to those of the bulge. These old stellar populations in the very galaxy outskirts are unexpected. We discuss their potential origin by considering the different gas phases of the source with the most prominent quenched ring, NGC4321, and argue for two main possibilities: 1) some environmental effect (e.g. starvation) or 2) the circumgalactic medium of sources outside of high-density clusters might have stopped to supply pristine gas to the galaxy (e.g. if its specific angular moment is too high for being accreted)., Comment: Accepted for publication on Astronomy and Astrophysics

We present the first measurement of HI mass of star-forming galaxies in different large scale structure environments from a blind survey at $z\sim 0.37$. In particular, we carry out a spectral line stacking analysis considering $2875$ spectra of colour-selected star-forming galaxies undetected in HI at $0.23 < z < 0.49$ in the COSMOS field, extracted from the MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies belonging to different subsamples depending on three different definitions of large scale structure environment: local galaxy overdensity, position inside the host dark matter halo (central, satellite, or isolated), and cosmic web type (field, filament, or knot). We first stack the full star-forming galaxy sample and find a robust HI detection yielding an average galaxy HI mass of $M_{\rm HI}=(8.12\pm 0.75)\times 10^9\, {\rm M}_\odot$ at $\sim 11.8\sigma$. Next, we investigate the different subsamples finding a negligible difference in $M_{\rm HI}$ as a function of the galaxy overdensity. We report an HI excess compared to the full sample in satellite galaxies ($M_{\rm HI}=(11.31\pm1.22)\times 10^9$, at $\sim 10.2 \sigma$) and in filaments ($M_{\rm HI}=(11.62\pm 0.90)\times 10^9$. Conversely, we report non-detections for the central and knot galaxies subsamples, which appear to be HI-deficient. We find the same qualitative results also when stacking in units of HI fraction ($f_{\rm HI}$). We conclude that the HI amount in star-forming galaxies at the studied redshifts correlates with the large scale structure environment., Comment: Accepted for publication in MNRAS. 15 figures, 3 tables

Observations at UV and optical wavelengths have revealed that galaxies at z~1-4 host star-forming regions, dubbed "clumps", which are believed to form due to the fragmentation of gravitationally unstable, gas-rich disks. However, the detection of the parent molecular clouds that give birth to such clumps is still possible only in a minority of galaxies, mostly at z~1. We investigated the [CII] and dust morphology of a z~3.4 lensed galaxy hosting four clumps detected in the UV continuum. We aimed to observe the [CII] emission of individual clumps that, unlike the UV, is not affected by dust extinction, to probe their nature and cold gas content. We conducted ALMA observations probing scales down to ~300 pc and detected three [CII] clumps. One (dubbed "NE") coincides with the brightest UV clump, while the other two ("SW" and "C") are not detected in the UV continuum. We do not detect the dust continuum. We converted the [CII] luminosity of individual clumps into molecular gas mass and found Mmol~10^8 Msun. By complementing it with the star formation rate (SFR) estimate from the UV continuum, we estimated the gas depletion time (tdep) of clumps and investigated their location in the Schmidt-Kennicutt plane. While the NE clump has a short tdep=0.16 Gyr, comparable with high-redshift starbursts, the SW and C clumps instead have longer tdep>0.65 Gyr and are likely probing the initial phases of star formation. The lack of dust continuum detection is consistent with the blue UV continuum slope estimated for this galaxy (beta~-2.5) and it indicates that dust inhomogeneities do not significantly affect the detection of UV clumps in this target. We pushed the observation of the cold gas content of individual clumps up to z~3.4 and showed that the [C II] line emission is a promising tracer of molecular clouds at high redshift, allowing the detection of clumps with a large range of depletion times., Comment: Accepted for publication in A&A

Dusty, distant, massive ($M_*\gtrsim 10^{11}\,\rm M_\odot$) galaxies are usually found to show a remarkable star-formation activity, contributing on the order of $25\%$ of the cosmic star-formation rate density at $z\approx3$--$5$, and up to $30\%$ at $z\sim7$ from ALMA observations. Nonetheless, they are elusive in classical optical surveys, and current near-infrared surveys are able to detect them only in very small sky areas. Since these objects have low space densities, deep and wide surveys are necessary to obtain statistically relevant results about them. Euclid will be potentially capable of delivering the required information, but, given the lack of spectroscopic features at these distances within its bands, it is still unclear if it will be possible to identify and characterize these objects. The goal of this work is to assess the capability of Euclid, together with ancillary optical and near-infrared data, to identify these distant, dusty and massive galaxies, based on broadband photometry. We used a gradient-boosting algorithm to predict both the redshift and spectral type of objects at high $z$. To perform such an analysis we make use of simulated photometric observations derived using the SPRITZ software. The gradient-boosting algorithm was found to be accurate in predicting both the redshift and spectral type of objects within the Euclid Deep Survey simulated catalog at $z>2$. In particular, we study the analog of HIEROs (i.e. sources with $H-[4.5]>2.25$), combining Euclid and Spitzer data at the depth of the Deep Fields. We found that the dusty population at $3\lesssim z\lesssim 7$ is well identified, with a redshift RMS and OLF of only $0.55$ and $8.5\%$ ($H_E\leq26$), respectively. Our findings suggest that with Euclid we will obtain meaningful insights into the role of massive and dusty galaxies in the cosmic star-formation rate over time., Comment: 18 pages, 13 figures, accepted in A&A

We present the star-formation-rate -- stellar-mass (SFR-M$_\ast$) relation for galaxies in the CEERS survey at $4.5\leq z\leq 12$. We model the \jwst\ and \hst\ rest-UV and rest-optical photometry of galaxies with flexible star-formation histories (SFHs) using \bagpipes. We consider SFRs averaged from the SFHs over 10~Myr (\sfrten) and 100~Myr (\sfrcen), where the photometry probes SFRs on these timescales, effectively tracing nebular emission lines in the rest-optical (on $\sim10$~Myr timescales) and the UV/optical continuum (on $\sim100$ Myr timescales). We measure the slope, normalization and intrinsic scatter of the SFR-M$_\ast$ relation, taking into account the uncertainty and the covariance of galaxy SFRs and $M_\ast$. From $z\sim 5-9$ there is larger scatter in the $\sfrten-M_\ast$ relation, with $\sigma(\log \sfrcen)=0.4$~dex, compared to the $\sfrcen-M_\ast$ relation, with $\sigma(\log \sfrten)=0.1$~dex. This scatter increases with redshift and increasing stellar mass, at least out to $z\sim 7$. These results can be explained if galaxies at higher redshift experience an increase in star-formation variability and form primarily in short, active periods, followed by a lull in star formation (i.e. ``napping'' phases). We see a significant trend in the ratio $R_\mathrm{SFR}=\log(\sfrten/\sfrcen)$ in which, on average, $R_\mathrm{SFR}$ decreases with increasing stellar mass and increasing redshift. This yields a star-formation ``duty cycle'' of $\sim40\%$ for galaxies with $\log M_\ast/M_\odot\geq 9.3$, at $z\sim5$, declining to $\sim20\%$ at $z\sim9$. Galaxies also experience longer lulls in star formation at higher redshift and at higher stellar mass, such that galaxies transition from periods of higher SFR variability at $z\gtrsim~6$ to smoother SFR evolution at $z\lesssim~4.5$., Comment: 28 pages, 12 figures, 2 Appendix figures

We present a sample of 1165 extreme emission-line galaxies (EELGs) at 45000 AA. JWST/NIRSpec spectroscopic observations of a subset (34) of the photometrically selected EELGs validate our selection method: all spectroscopically observed EELGs confirm our photometric identification of extreme emission, including some cases where the SED-derived photometric redshifts are incorrect. We find that the medium-band F410M filter in CEERS is particularly efficient at identifying EELGs, both in terms of including emission lines in the filter and in correctly identifying the continuum between Hb + [OIII] and Ha in the neighboring broad-band filters. We present examples of EELGs that could be incorrectly classified at ultra-high redshift (z>12) as a result of extreme Hb + [OIII] emission blended across the reddest photometric filters. We compare the EELGs to the broader (sub-extreme) galaxy population in the same redshift range and find that they are consistent with being the bluer, high equivalent width tail of a broader population of emission-line galaxies. The highest-EW EELGs tend to have more compact emission-line sizes than continuum sizes, suggesting that active galactic nuclei are responsible for at least some of the most extreme EELGs. Photometrically inferred emission-line ratios are consistent with ISM conditions with high ionization and moderately low metallicity, consistent with previous spectroscopic studies., Comment: 24 pages, 25 figures, submitted to ApJ

We present the discovery of the most distant OH megamaser to be observed in the main lines, using data from the MeerKAT International Giga-Hertz Tiered Extragalactic Exploration (MIGHTEE) survey. At a newly measured redshift of $z = 0.7092$, the system has strong emission in both the 1665MHz ($L \approx 2500$ L$_{\odot}$) and 1667 MHz ($L \approx 4.5\times10^4$ L$_{\odot}$) transitions, with both narrow and broad components. We interpret the broad line as a high-velocity-dispersion component of the 1667 MHz transition, with velocity $v \sim 330$km s$^{-1}$ with respect to the systemic velocity. The host galaxy has a stellar mass of $M_{\star} = 2.95 \times 10^{10}$ M$_{\odot}$ and a star-formation rate of SFR = 371 M$_{\odot}$yr$^{-1}$, placing it $\sim 1.5$dex above the main sequence for star-forming galaxies at this redshift, and can be classified as an ultra-luminous infrared galaxy. Alongside the optical imaging data, which exhibits evidence for a tidal tail, this suggests that the OH megamaser arises from a system that is currently undergoing a merger, which is stimulating star formation and providing the necessary conditions for pumping the OH molecule to saturation. The OHM is likely to be lensed, with a magnification factor of $\sim 2.5$, and perhaps more if the maser emitting region is compact and suitably offset relative to the centroid of its host galaxy's optical light. This discovery demonstrates that spectral line mapping with the new generation of radio interferometers may provide important information on the cosmic merger history of galaxies., Comment: 11 pages, accepted for publication in MNRAS

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 present a sample of 88 candidate z~8.5-14.5 galaxies selected from the completed NIRCam imaging from the Cosmic Evolution Early Release Science (CEERS) survey. These data cover ~90 arcmin^2 (10 NIRCam pointings) in six broad-band and one medium-band imaging filter. With this sample we confirm at higher confidence early JWST conclusions that bright galaxies in this epoch are more abundant than predicted by most theoretical models. We construct the rest-frame ultraviolet luminosity functions at z~9, 11 and 14, and show that the space density of bright (M_UV=-20) galaxies changes only modestly from z~14 to z~9, compared to a steeper increase from z~8 to z~4. While our candidates are photometrically selected, spectroscopic followup has now confirmed 13 of them, with only one significant interloper, implying that the fidelity of this sample is high. Successfully explaining the evidence for a flatter evolution in the number densities of UV-bright z>10 galaxies may thus require changes to the dominant physical processes regulating star formation. While our results indicate that significant variations of dust attenuation with redshift are unlikely to be the dominant factor at these high redshifts, they are consistent with predictions from models which naturally have enhanced star-formation efficiency and/or stochasticity. An evolving stellar initial mass function could also bring model predictions into better agreement with our results. Deep spectroscopic followup of a large sample of early galaxies can distinguish between these competing scenarios., Comment: Submitted to ApJL. Main paper is 33 pages, 13 figures, 5 tables. Two appendices with additional figures and tables

We have carried out HST snapshot observations at 1.1 $\mu$m of 281 candidate strongly lensed galaxies identified in the wide-area extragalactic surveys conducted with the Herschel space observatory. Our candidates comprise systems with flux densities at $500\,\mu$m$ S_{500}\geq 80$ mJy. We model and subtract the surface brightness distribution for 130 systems, where we identify a candidate for the foreground lens candidate. After combining visual inspection, archival high-resolution observations, and lens subtraction, we divide the systems into different classes according to their lensing likelihood. We confirm 65 systems to be lensed. Of these, 30 are new discoveries. We successfully perform lens modelling and source reconstruction on 23 systems, where the foreground lenses are isolated galaxies and the background sources are detected in the HST images. All the systems are successfully modelled as a singular isothermal ellipsoid. The Einstein radii of the lenses and the magnifications of the background sources are consistent with previous studies. However, the background source circularised radii (between 0.34 kpc and 1.30 kpc) are $\sim$3 times smaller than the ones measured in the sub-mm/mm for a similarly selected and partially overlapping sample. We compare our lenses with those in the SLACS survey, confirming that our lens-independent selection is more effective at picking up fainter and diffuse galaxies and group lenses. This sample represents the first step towards characterising the near-IR properties and stellar masses of the gravitationally lensed dusty star-forming galaxies., Comment: 57 pages, 18 figures, 11 tables. Accepted for publication in MNRAS

We test the relationship between UV-derived star formation rates (SFRs) and the 7.7 ${\mu}$m polycyclic aromatic hydrocarbon (PAH) luminosities from the integrated emission of galaxies at z ~ 0 - 2. We utilize multi-band photometry covering 0.2 - 160 ${\mu}$m from HST, CFHT, JWST, Spitzer, and Herschel for galaxies in the Cosmic Evolution Early Release Science (CEERS) Survey. We perform spectral energy distribution (SED) modeling of these data to measure dust-corrected far-UV (FUV) luminosities, $L_{FUV}$, and UV-derived SFRs. We then fit SED models to the JWST/MIRI 7.7 - 21 ${\mu}$m CEERS data to derive rest-frame 7.7 ${\mu}$m luminosities, $L_{770}$, using the average flux density in the rest-frame MIRI F770W bandpass. We observe a correlation between $L_{770}$ and $L_{FUV}$, where log $L_{770}$ is proportional to (1.27+/-0.04) log $L_{FUV}$. $L_{770}$ diverges from this relation for galaxies at lower metallicities, lower dust obscuration, and for galaxies dominated by evolved stellar populations. We derive a "single-wavelength" SFR calibration for $L_{770}$ which has a scatter from model estimated SFRs (${{\sigma}_{{\Delta}SFR}}$) of 0.24 dex. We derive a "multi-wavelength" calibration for the linear combination of the observed FUV luminosity (uncorrected for dust) and the rest-frame 7.7 ${\mu}$m luminosity, which has a scatter of ${{\sigma}_{{\Delta}SFR}}$ = 0.21 dex. The relatively small decrease in ${\sigma}$ suggests this is near the systematic accuracy of the total SFRs using either calibration. These results demonstrate that the rest-frame 7.7 ${\mu}$m emission constrained by JWST/MIRI is a tracer of the SFR for distant galaxies to this accuracy, provided the galaxies are dominated by star-formation with moderate-to-high levels of attenuation and metallicity., Comment: 20 pages, 11 figures, 2 tables, submitted to ApJ

Surprise Machines is a project of experimental museology that sets out to visualize the entire image collection of the Harvard Art Museums, intending to open up unexpected vistas on more than 200,000 objects usually inaccessible to visitors. Part of the exhibition Curatorial A(i)gents organized by metaLAB (at) Harvard, the project explores the limits of artificial intelligence to display a large set of images and create surprise among visitors. To achieve such a feeling of surprise, a choreographic interface was designed to connect the audience's movement with several unique views of the collection., Comment: 14 pages and 7 figures

We take advantage of the NIRCam photometric observations available as part of the Cosmic Evolution Early Release Science survey (CEERS) to identify and analyse very red sources in an effort to discover very dusty star forming galaxies. We select red galaxies as objects with a S/N>3 at 4.4 $\mu$m and a S/N<2 in all JWST and HST filters at $\lambda\leq2\mu$m, which corresponds to [F200W]-[F444W]>1.2 considering CEERS depths. This selection is ideal to identify very dusty (Av>1 mag) galaxies with stellar masses between $10^6$ to $10^{10}\, \rm M_{\odot}$ at z<5, more massive dusty galaxies at z=5-18 and galaxies at z>18 due to the Lyman absorption, independently of their dust extinction. Our sample of F200W-dropouts contains no strong candidates at z>6.5, instead it consists almost completely (~81%) of z<2 low-mass galaxies, with a median stellar mass of $10^{7.3} \rm M_{\odot}$. These galaxies show an exceptional dust extinction with median value of Av=4.9 mag, completely unexpected given their low stellar mass. The remaining galaxies, which are at z<6.5, show similar large dust extinction (Av>1), but they are generally more massive $>10^{7.5}\rm M_{\odot}$., Comment: 30 pages, 1 table, 19 figures, accepted for publication in A&A

This work focuses on the pilot run of a simulation campaign aimed at investigating the spectroscopic capabilities of the Euclid Near-Infrared Spectrometer and Photometer (NISP), in terms of continuum and emission line detection in the context of galaxy evolutionary studies. To this purpose we constructed, emulated, and analysed the spectra of 4992 star-forming galaxies at $0.3 \leq z \leq 2.5$ using the NISP pixel-level simulator. We built the spectral library starting from public multi-wavelength galaxy catalogues, with value-added information on spectral energy distribution (SED) fitting results, and from Bruzual and Charlot (2003) stellar population templates. Rest-frame optical and near-IR nebular emission lines were included using empirical and theoretical relations. We inferred the 3.5$\sigma$ NISP red grism spectroscopic detection limit of the continuum measured in the $H$ band for star-forming galaxies with a median disk half-light radius of \ang{;;0.4} at magnitude $H= 19.5\pm0.2\,$AB$\,$mag for the Euclid Wide Survey and at $H = 20.8\pm0.6\,$AB$\,$mag for the Euclid Deep Survey. We found a very good agreement with the red grism emission line detection limit requirement for the Wide and Deep surveys. We characterised the effect of the galaxy shape on the detection capability of the red grism and highlighted the degradation of the quality of the extracted spectra as the disk size increases. In particular, we found that the extracted emission line signal to noise ratio (SNR) drops by $\sim\,$45$\%$ when the disk size ranges from \ang{;;0.25} to \ang{;;1}. These trends lead to a correlation between the emission line SNR and the stellar mass of the galaxy and we demonstrate the effect in a stacking analysis unveiling emission lines otherwise too faint to detect., Comment: 24 pages, 21 figures

The evolution of the luminosity function (LF) of Active Galactic Nuclei (AGNs) at $z \gtrsim 5$ represents a key constraint to understand their contribution to the ionizing photon budget necessary to trigger the last phase transition in the Universe, i.e. the epoch of Reionization. Recent searches for bright high-z AGNs suggest that the space densities of this population at $z>4$ has to be revised upwards, and sparks new questions about their evolutionary paths. Gas accretion is the key physical mechanism to understand both the distribution of luminous sources and the growth of central Super-Massive Black Holes (SMBHs). In this work, we model the high-z AGN-LF assuming that high-z luminous AGN shine at their Eddington limit: we derive the expected evolution as a function of the ``duty-cycle'' ($f_{\rm dc}$), i.e. the fraction of life-time that a given SMBH spends accreting at the Eddington rate. Our results show that intermediate values ($f_{\rm dc} \simeq 0.1$) predict the best agreement with the ionizing background and photoionization rate, but do not provide enough ionizing photons to account for the observed evolution of the hydrogen neutral fraction. Smaller values ($f_{\rm dc} \lesssim 0.05$) are required for AGNs to be the dominant population responsible for Hydrogen reionization in the Early Universe. We then show that this low-$f_{\rm dc}$ evolution can be reconciled with the current constraints on Helium reionization, although it implies a relatively large number of inactive SMBHs at $z\gtrsim5$, in tension with SMBH growth models based on heavy seeding., Comment: 11 pages, 5 figures, MNRAS accepted

Observational studies have widely demonstrated that galaxy physical properties are strongly affected by the surrounding environment. On one side, gas inflows provide galaxies with new fuel for star formation. On the other side, the high temperatures and densities of the medium are expected to induce quenching in the star formation. Observations of large structures, in particular filaments at the cluster outskirts (r>2r$_{200}$), are currently limited to the low redshift Universe. We present a multi-band dataset for the cluster MACS J0416.1-2403 (z=0.397), observed in the context of the Galaxy Assembly as a function of Mass and Environment with VST (VST-GAME) survey. The project aims at gathering deep ($r$<24.4) and wide (20x20Mpc$^2$) observations at optical wavelengths for six massive galaxy clusters at 0.2

The origin of the Lyman-${\alpha}$ (Ly${\alpha}$) emission in galaxies is a long-standing issue: despite several processes known to originate this line (e.g. AGN, star formation, cold accretion, shock heating), it is difficult to discriminate among these phenomena based on observations. Recent studies have suggested that the comparison of the ultraviolet (UV) and optical properties of these sources could solve the riddle. For this reason, we investigate the rest-frame UV and optical properties of A2895b, a strongly lensed Ly${\alpha}$-emitter at redshift z ~ 3.7. From this study, we find that our target is a compact (r ~ 1.2 pkpc) star-forming (star formation rate ~ 11 M$_{\odot}$/yr) galaxy having a young stellar population. Interestingly, we measure a high ratio of the H${\beta}$ and the UV continuum monochromatic luminosities (L(H${\beta}$)/L(UV) ~ 100). Based on tracks of theoretical stellar models (Starburst99, BPASS), we can only partially explain this result by assuming a recent (< 10 Myr), bursty episode of star-formation and considering models characterised by binary stars, a top-heavy initial-mass function (IMF) and sub-solar metallicities (Z < 0.01 Z$_{\odot}$). These assumptions also explain the observed low (C/O) abundance of our target (~ 0.23(C/O)$_{\odot}$). By comparing the UV and optical datasets, we find that the Ly${\alpha}$ and UV continuum are more extended (x2) than the Balmer lines, and that the peak of the Ly${\alpha}$ is offset (~ 0.6 pkpc). The multi-wavelength results of our analysis suggest that the observed Ly${\alpha}$ emission originates from a recent star-formation burst, likely taking place in an off-centre clump., Comment: 18 pages (including Appendix), 5 figures. Paper accepted for publication on MNRAS. Accepted 2022 November 2. Received 2022 October 22; in original form 2022 August 3

We study the $M_{\rm HI}-M_{\star}$ relation over the last billion years using the MIGHTEE-HI sample. We first model the upper envelope of the $M_{\rm HI}-M_{\star}$ relation with a Bayesian technique applied to a total number of 249 HI-selected galaxies, without binning the datasets, while taking account of the intrinsic scatter. We fit the envelope with both linear and non-linear models, and find that the non-linear model is preferred over the linear one with a measured transition stellar mass of $\log_{10}(M_\star$/$M_{\odot})$ = $9.15\pm0.87$, beyond which the slope flattens. This finding supports the view that the lack of HI gas is ultimately responsible for the decreasing star formation rate observed in the massive main-sequence galaxies. For spirals alone, which are biased towards the massive galaxies in our sample, the slope beyond the transition mass is shallower than for the full sample, indicative of distinct gas processes ongoing for the spirals/high-mass galaxies from other types with lower stellar masses. We then create mock catalogues for the MIGHTEE-HI detections and non-detections with two main galaxy populations of late- and early-type galaxies to measure the underlying $M_{\rm HI}-M_{\star}$ relation. We find that the turnover in this relation persists whether considering the two galaxy populations as a whole or separately. We note that an underlying linear relation could mimic this turnover in the observed scaling relation, but a model with a turnover is strongly preferred. Measurements on the logarithmic average of HI masses against the stellar mass are provided as a benchmark for future studies., Comment: Revised manuscript accepted for publication in MNRAS

Mapping the average AGN luminosity across galaxy populations and over time encapsulates important clues on the interplay between supermassive black hole (SMBH) and galaxy growth. This paper presents the demography, mean power and cosmic evolution of radio AGN across star-forming galaxies (SFGs) of different stellar masses (${M_{*}}$). We exploit deep VLA-COSMOS 3 GHz data to build the rest-frame 1.4 GHz AGN luminosity functions at 0.1$\leq$$z$$\leq$4.5 hosted in SFGs. Splitting the AGN luminosity function into different ${M_{*}}$ bins reveals that, at all redshifts, radio AGN are both more frequent and more luminous in higher ${M_*}$ than in lower ${M_*}$ galaxies. The cumulative kinetic luminosity density exerted by radio AGN in SFGs peaks at $z$$\sim$2, and it is mostly driven by galaxies with 10.5$\leq$$\log$(${M_{*}}$/${M_{\odot}}$)$<$11. Averaging the cumulative radio AGN activity across all SFGs at each (${M_{*}}$,$z$) results in a "radio-AGN main sequence" that links the time-averaged radio-AGN power $\langle$$L_{1.4}^{{AGN}}$$\rangle$ and galaxy stellar mass, in the form: $\log$$\langle$[$L_{1.4}^{{AGN}}$/ W Hz$^{-1}]\rangle$ = (20.97$\pm$0.16) + (2.51$\pm$0.34)$\cdot$$\log$(1+$z$) + (1.41$\pm$0.09)$\cdot$($\log$[${M_{*}}$/${M_{\odot}}$] -10). The super-linear dependence on ${M_{*}}$, at fixed redshift, suggests enhanced radio-AGN activity in more massive SFGs, as compared to star formation. We ascribe this enhancement to both a higher radio AGN duty cycle and a brighter radio-AGN phase in more massive SFGs. A remarkably consistent ${M_{*}}$ dependence is seen for the evolving X-ray AGN population in SFGs. This similarity is interpreted as possibly driven by secular cold gas accretion fueling both radio and X-ray AGN activity in a similar fashion over the galaxy's lifetime., Comment: Accepted for publication in A&A. 20 pages + Appendices

We present a new prediction of the luminosity functions of the [CII] line at 158 $\mu$m, of the CO lines from J=0 to J=24, and of the molecular gas mass density up to z=10, using the Spectro-Photometric Realisations of Infrared-selected Targets at all-z (SPRITZ) simulation (Bisigello et al. 2021). We update the state-of-the-art phenomenological simulation SPRITZ to include both the CO ($J\leq24$) and [CII] line luminosities. This has been performed using different empirical and theoretical relations to convert the total infrared luminosity (or star formation rate) to [CII] or CO luminosity. The resulting line luminosity functions have been compared for validation with a large set of observations available in the literature. We then used the derived CO and [CII] line luminosities to estimate the molecular gas mass density and compare it with available observations. The CO and [CII] luminosity functions presented here are well in agreement with all the available observations. In particular, the best results for [CII] are obtained deriving the [CII] luminosity directly from the star formation rate, but considering a dependence of this relation on the gas metallicity. For all the CO luminosity functions, the estimates favoured by the data are derived considering different relations, depending on the ionisation mechanism dominating each galaxy, i.e. star formation or active galactic nuclei, and, moreover, deriving the $J\geq4$ CO lines directly from the [CII] luminosity. However, further data are necessary to fully discriminate between models. Finally, the best agreement with observations of the molecular gas mass density are derived by converting the [CII] luminosity to H2 mass, using a [CII]-to-H2 conversion ~130 $\rm M_{\odot}/{\rm L}_{\odot}$. All the line luminosity functions, useful for planning and interpreting future observations, are made publicly available., Comment: 18 pages, 16 figures, 8 tables, accepted for publication in A&A

A wealth of extragalactic populations completely missed at UV-optical wavelengths has been identified in the last decade, combining the deepest HST and Spitzer observations. These dark sources are thought to be very dusty and star-forming systems at 3

We present the first measurements of HI galaxy scaling relations from a blind survey at $z>0.15$. We perform spectral stacking of 9023 spectra of star-forming galaxies undetected in HI at $0.23

We present an assessment of the accuracy of common operations performed in $21$-cm spectral line stacking experiments. To this end, we generate mock interferometric data surveying the 21-cm emission at frequency $1310<\nu<1420$ MHz ($0.005

(Abridged) The Euclid mission is expected to discover thousands of z>6 galaxies in three Deep Fields, which together will cover a ~40 deg2 area. However, the limited number of Euclid bands and availability of ancillary data could make the identification of z>6 galaxies challenging. In this work, we assess the degree of contamination by intermediate-redshift galaxies (z=1-5.8) expected for z>6 galaxies within the Euclid Deep Survey. This study is based on ~176,000 real galaxies at z=1-8 in a ~0.7 deg2 area selected from the UltraVISTA ultra-deep survey, and ~96,000 mock galaxies with 25.3$\leq$H<27.0, which altogether cover the range of magnitudes to be probed in the Euclid Deep Survey. We simulate Euclid and ancillary photometry from the fiducial, 28-band photometry, and fit spectral energy distributions (SEDs) to various combinations of these simulated data. Our study demonstrates that identifying z>6 with Euclid data alone will be very effective, with a z>6 recovery of 91(88)% for bright (faint) galaxies. For the UltraVISTA-like bright sample, the percentage of z=1-5.8 contaminants amongst apparent z>6 galaxies as observed with Euclid alone is 18%, which is reduced to 4(13)% by including ultra-deep Rubin (Spitzer) photometry. Conversely, for the faint mock sample, the contamination fraction with Euclid alone is considerably higher at 39%, and minimized to 7% when including ultra-deep Rubin data. For UltraVISTA-like bright galaxies, we find that Euclid (I-Y)>2.8 and (Y-J)<1.4 colour criteria can separate contaminants from true z>6 galaxies, although these are applicable to only 54% of the contaminants, as many have unconstrained (I-Y) colours. In the most optimistic scenario, these cuts reduce the contamination fraction to 1% whilst preserving 81% of the fiducial z>6 sample. For the faint mock sample, colour cuts are infeasible., Comment: 27 pages, 14 figures, 9 tables; version accepted for publication by A&A

By compiling a comprehensive census of literature studies, we investigate the evolution of the Main Sequence (MS) of star-forming galaxies (SFGs) in the widest range of redshift ($0 < z < 6$) and stellar mass ($10^{8.5}-10^{11.5}$ $M_{\odot}$) ever probed. We convert all observations to a common calibration and find a remarkable consensus on the variation of the MS shape and normalization across cosmic time. The relation exhibits a curvature towards the high stellar masses at all redshifts. The best functional form is governed by two parameters: the evolution of the normalization and the turnover mass ($M_0(t)$), which both evolve as a power law of the Universe age. The turn-over mass determines the MS shape. It marginally evolves with time, making the MS slightly steeper towards $z\sim4-6$. At stellar masses below $M_0(t)$, SFGs have a constant specific SFR (sSFR), while above $M_0(t)$ the sSFR is suppressed. We find that the MS is dominated by central galaxies. This allows to turn $M_0(t)$ into the corresponding host halo mass. This evolves as the halo mass threshold between cold and hot accretion regimes, as predicted by the theory of accretion, where the central galaxy is fed or starved of cold gas supply, respectively. We, thus, argue that the progressive MS bending as a function of the Universe age is caused by the lower availability of cold gas in halos entering the hot accretion phase, in addition to black hole feedback. We also find qualitatively the same trend in the largest sample of star-forming galaxies provided by the IllustrisTNG simulation. Nevertheless, we still note large quantitative discrepancies with respect to observations, in particular at the high mass end. These can not be easily ascribed to biases or systematics in the observed SFRs and the derived MS., Comment: Accepted for publication by MNRAS, corrected table 2

A mild correlation exists in active galaxies between the mean black hole accretion, as traced by the mean X-ray luminosity $\left<{\rm L}_{\rm X}\right>$, and the host galaxy stellar mass M$_*$, characterised by a normalisation steadily decreasing with cosmic time and lower in more quiescent galaxies. We create comprehensive semi-empirical mock catalogues of active black holes to pin down which parameters control the shape and evolution of the $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation of X-ray detected active galaxies. We find that the normalisation of the $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation is largely independent of the fraction of active galaxies (the duty cycle), but strongly dependent on the mean Eddington ratio, when adopting a constant underlying M$_{\rm BH}-{\rm M}_*$ relation as suggested by observational studies. The data point to a decreasing mean Eddington ratio with cosmic time and with galaxy stellar mass at fixed redshift. Our data can be reproduced by black holes and galaxies evolving on similar M$_{\rm BH}-{\rm M}_*$ relations but progressively decreasing their average Eddington ratios, mean X-ray luminosities, and specific star formation rates, when moving from the starburst to the quiescent phase. Models consistent with the observed $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation and independent measurements of the mean Eddington ratios, are characterised by M$_{\rm BH}-{\rm M}_*$ relations lower than those derived from dynamically measured local black holes. Our results point to the $\left<{\rm L}_{\rm X}\right>-{\rm M}_*$ relation as a powerful diagnostic to: 1) probe black hole-galaxy scaling relations and the level of accretion onto black holes; 2) efficiently break the degeneracies between duty cycles and accretion rates in cosmological models of black holes., Comment: Accepted for publication on MNRAS. 11 pages, 5 figures

The Star Formation Rate Density (SFRD) history of the Universe is well constrained up to redshift $z \sim 2$. At earlier cosmic epochs, the picture has been largely inferred from UV-selected galaxies (e.g. Lyman-break galaxies, LBGs). However, LBGs' inferred SFRs strongly depend on the assumed dust extinction correction, which is not well-constrained at high-$z$, while observations in the radio domain are not affected by this issue. In this work we measure the SFRD from a 1.4 GHz-selected sample of $\sim$600 galaxies in the GOODS-N field up to redshift $\sim 3.5$. We take into account the contribution of Active Galactic Nuclei from the Infrared-Radio correlation. We measure the radio luminosity function, fitted with a modified Schechter function, and derive the SFRD. The cosmic SFRD shows a rise up to $z \sim 2$ and then an almost flat plateau up to $z \sim 3.5$. Our SFRD is in agreement with the ones from other FIR/radio surveys and a factor 2 higher than those from LBG samples. We also estimate that galaxies lacking a counterpart in the HST/WFC3 H-band ($H$-dark) make up $\sim 25\%$ of the $\phi$-integrated SFRD relative to the full sample at z $\sim 3.2$, and up to $58\%$ relative to LBG samples., Comment: 20 pages, 11 figures. Accepted for publication in ApJ

We use a sample of 706 galaxies, selected as [OII]$\lambda$3727 ([OII]) emitters in the Survey for High-$z$ Absorption Red and Dead Sources (SHARDS) on the CANDELS/GOODS-N field, to study the differential attenuation of the nebular emission with respect to the stellar continuum. The sample includes only galaxies with a counterpart in the infrared and $\mathrm{log}_{10}(M_{*}/\mathrm{M}_{\odot})$ $>$ 9, over the redshift interval 0.3 $\lesssim$ $z$ $\lesssim$ 1.5. Our methodology consists in the comparison of the star formation rates inferred from [OII] and H$\alpha$ emission lines with a robust quantification of the total star-forming activity (${SFR}_{\mathrm{TOT}}$) that is independently estimated based on both infrared and ultraviolet (UV) luminosities. We obtain $f$$=$$E(B-V)_{\mathrm{stellar}}$/$E(B-V)_{\mathrm{nebular}}$ $=$ 0.69$^{0.71}_{0.69}$ and 0.55$^{0.56}_{0.53}$ for [OII] and H$\alpha$, respectively. Our resulting $f$-factors display a significant positive correlation with the UV attenuation and shallower or not-significant trends with the stellar mass, the $SFR_{\mathrm{TOT}}$, the distance to the main sequence, and the redshift. Finally, our results favour an average nebular attenuation curve similar in shape to the typical dust curve of local starbursts., Comment: 23 pages, 10 figures, 3 tables; accepted for publication in MNRAS

Mid- and far-infrared (IR) photometric and spectroscopic observations are fundamental to a full understanding of the dust-obscured Universe and the evolution of both star formation and black hole accretion in galaxies. In this work, using the specifications of the SPace Infrared telescope for Cosmology and Astrophysics (SPICA) as a baseline, we investigate the capability to study the dust-obscured Universe of mid- and far-IR photometry at 34 and 70 $\mu$m and low-resolution spectroscopy at 17-36 $\mu$m using the state-of-the-art Spectro-Photometric Realisations of Infrared-selected Targets at all-z (SPRITZ) simulation. This investigation is also compared to the expected performance of the Origins Space Telescope and the Galaxy Evolution Probe. The photometric view of the Universe of a SPICA-like mission could cover not only bright objects (e.g. L$_{IR}$>10$^{12}$L$_\odot$) up to z=10, but also normal galaxies(L$_{IR}$<10$^{11}$L$_\odot$) up to z$\sim$4. At the same time, the spectroscopic observations of such mission could also allow us to estimate the redshifts and study the physical properties for thousands of star-forming galaxies and active galactic nuclei by observing the polycyclic aromatic hydrocarbons and a large set of IR nebular emission lines. In this way, a cold, 2.5-m size space telescope with spectro-photometric capability analogous to SPICA, could provide us with a complete three-dimensional (i.e. images and integrated spectra) view of the dust-obscured Universe and the physics governing galaxy evolution up to z$\sim$4., Comment: 22 pages, 17 figures, 4 tables, accepted for publication in PASA

The identification of an emission line is unambiguous when multiple spectral features are clearly visible in the same spectrum. However, in many cases, only one line is detected, making it difficult to correctly determine the redshift. We developed a freely available unsupervised machine-learning algorithm based on unbiased topology (UMLAUT) that can be used in a very wide variety of contexts, including the identification of single emission lines. To this purpose, the algorithm combines different sources of information, such as the apparent magnitude, size and color of the emitting source, and the equivalent width and wavelength of the detected line. In each specific case, the algorithm automatically identifies the most relevant ones (i.e., those able to minimize the dispersion associated with the output parameter). The outputs can be easily integrated into different algorithms, allowing us to combine supervised and unsupervised techniques and increasing the overall accuracy. We tested our software on WISP (WFC3 IR Spectroscopic Parallel) survey data. WISP represents one of the closest existing analogs to the near-IR spectroscopic surveys that are going to be performed by the future Euclid and Roman missions. These missions will investigate the large-scale structure of the universe by surveying a large portion of the extragalactic sky in near-IR slitless spectroscopy, detecting a relevant fraction of single emission lines. In our tests, UMLAUT correctly identifies real lines in 83.2% of the cases. The accuracy is slightly higher (84.4%) when combining our unsupervised approach with a supervised approach we previously developed.

Giant star-forming regions (clumps) are widespread features of galaxies at $z \approx 1-4$. Theory predicts that they can play a crucial role in galaxy evolution if they survive to stellar feedback for > 50 Myr. Numerical simulations show that clumps' survival depends on the stellar feedback recipes that are adopted. Up to date, observational constraints on both clumps' outflows strength and gas removal timescale are still uncertain. In this context, we study a line-emitting galaxy at redshift $z \simeq 3.4$ lensed by the foreground galaxy cluster Abell 2895. Four compact clumps with sizes $\lesssim$ 280 pc and representative of the low-mass end of clumps' mass distribution (stellar masses $\lesssim 2\times10^8\ {\rm M}_\odot$) dominate the galaxy morphology. The clumps are likely forming stars in a starbursting mode and have a young stellar population ($\sim$ 10 Myr). The properties of the Lyman-$\alpha$ (Ly$\alpha$) emission and nebular far-ultraviolet absorption lines indicate the presence of ejected material with global outflowing velocities of $\sim$ 200-300 km/s. Assuming that the detected outflows are the consequence of star formation feedback, we infer an average mass loading factor ($\eta$) for the clumps of $\sim$ 1.8 - 2.4 consistent with results obtained from hydro-dynamical simulations of clumpy galaxies that assume relatively strong stellar feedback. Assuming no gas inflows (semi-closed box model), the estimates of $\eta$ suggest that the timescale over which the outflows expel the molecular gas reservoir ($\simeq 7\times 10^8\ \text{M}_\odot$) of the four detected low-mass clumps is $\lesssim$ 50 Myr., Comment: 19 pages, 11 figures + Appendix. Published in MNRAS

Using the CLASH-VLT survey, we assembled an unprecedented sample of 1234 spectroscopically confirmed members in Abell~S1063, finding a dynamically complex structure at z_cl=0.3457 with a velocity dispersion \sigma_v=1380 -32 +26 km s^-1. We investigate cluster environmental and dynamical effects by analysing the projected phase-space diagram and the orbits as a function of galaxy spectral properties. We classify cluster galaxies according to the presence and strength of the [OII] emission line, the strength of the H$\delta$ absorption line, and colours. We investigate the relationship between the spectral classes of galaxies and their position in the projected phase-space diagram. We analyse separately red and blue galaxy orbits. By correlating the observed positions and velocities with the projected phase-space constructed from simulations, we constrain the accretion redshift of galaxies with different spectral types. Passive galaxies are mainly located in the virialised region, while emission-line galaxies are outside r_200, and are accreted later into the cluster. Emission-lines and post-starbursts show an asymmetric distribution in projected phase-space within r_200, with the first being prominent at Delta_v/sigma <~-1.5$, and the second at Delta_v/ sigma >~ 1.5, suggesting that backsplash galaxies lie at large positive velocities. We find that low-mass passive galaxies are accreted in the cluster before the high-mass ones. This suggests that we observe as passives only the low-mass galaxies accreted early in the cluster as blue galaxies, that had the time to quench their star formation. We also find that red galaxies move on more radial orbits than blue galaxies. This can be explained if infalling galaxies can remain blue moving on tangential orbits., Comment: 24 pages, 28 figures, 3 tables, Accepted for publication in Astronomy and Astrophysics

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

Galaxy mergers are thought to be one of the main mechanisms of the mass assembly of galaxies. Recently, many works have suggested a possible increase in the fraction of major mergers in the early Universe, reviving the debate on which processes (e.g., cold accretion, star formation, mergers) most contribute to the mass build-up of galaxies through cosmic time. To estimate the importance of major mergers in this context, we make use of the new data collected by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE), which observed the [CII] 158 $\mu$m emission line from a sample of 75 main-sequence star-forming galaxies at 4.4 < z < 5.9. We used, for the first time, the morpho-kinematic information provided by the [CII] emission to obtain the fraction of major mergers ($f_{MM}$) at z~5. By adopting different prescriptions for the merger timescales ($T_{MM}$), we converted this fraction into the merger rate per galaxy ($R_{MM}$) and per volume ($\Gamma_{MM}$). We then combined our results with those at lower redshifts from the literature, computing the cosmic evolution of the merger fraction. This is described by a rapid increase from z~0 to higher redshifts, a peak at z~3, and a slow decrease towards earlier epochs. Depending on the timescale prescription used, this fraction translates into a merger rate ranging between ~0.1 and ~4.0 Gyr$^{-1}$ at z~5. Finally, we compare the specific star formation and star-formation rate density with the analogous quantities from major mergers. Our new ALPINE data reveal the presence of a significant merging activity in the early Universe. However, whether this population of mergers can provide a relevant contribution to the galaxy mass assembly at these redshifts and through the cosmic epochs is strongly dependent on the assumption of the merger timescale., Comment: Accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics

In the BH-galaxy co-evolution framework, most of the star-formation (SF) and the black hole (BH) accretion is expected to take place in highly obscured conditions. Thus, obscured AGN are difficult to identify in optical or X-ray bands, but shine bright in the IR. Moreover, X-ray background (XRB) synthesis models predict that a large fraction of the yet-unresolved XRB is due to the most obscured (Compton thick, CT) of these AGN. In this work, we investigate the synergies between putative IR missions (using SPICA, proposed for ESA/M5 but withdrawn in October 2020, and Origins Space Telescope, OST, as `templates') and the X-ray mission Athena, which should fly in early 2030s, in detecting and characterizing AGN, with a particular focus on the most obscured ones. Using an XRB synthesis model, we estimated the number of AGN and the number of those which will be detected in the X-rays. For each AGN we associated an optical-to-FIR SED from observed AGN with both X-ray data and SED decomposition, and used these SEDs to check if the AGN will be detected by SPICA-like or OST at IR wavelengths. We expect that, with the deepest Athena and SPICA-like (or OST) surveys, we will be able to detect in the IR more than $90\,\%$ of all the AGN (down to L$_{2-10\text{keV}} \sim 10^{42}\,$erg/s and up to $z \sim 10$) predicted by XRB synthesis modeling, and we will detect at least half of them in the X-rays. Athena will be extremely powerful in detecting and discerning moderate- and high-luminosity AGN. We find that the most obscured and elusive CT-AGN will be exquisitely sampled by SPICA-like mission or OST and that Athena will allow a fine characterization of the most-luminous ones. This will provide a significant step forward in the process of placing stronger constraints on the yet-unresolved XRB and investigating the BH accretion rate evolution up to very high redshift ($z \ge 4$)., Comment: Accepted for publication in PASA

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

Studying the resolved stellar populations of the different structural components which build massive galaxies directly unveils their assembly history. We aim at characterizing the stellar population properties of a representative sample of bulges and pure spheroids in massive galaxies ($M_{\star}>10^{10}$ M$_{\odot}$) in the GOODS-N field. We take advantage of the spectral and spatial information provided by SHARDS and HST data to perform the multi-image spectro-photometrical decoupling of the galaxy light. We derive the spectral energy distribution separately for bulges and disks in the redshift range $0.14

We use the SPace Infrared telescope for Cosmology and Astrophysics (SPICA) project as a template to demonstrate how deep spectrophotometric surveys covering large cosmological volumes over extended fields (1-15 square degrees) with a mid-IR imaging spectrometer (17-36 micron) in conjunction with deep 70 micron photometry with a far-IR camera, at wavelengths which are not affected by dust extinction can answer the most crucial questions in current galaxy evolution studies. A SPICA-like mission will be able for the first time to provide an unobscured three dimensional (3-D, i.e. x, y and redshift z) view of galaxy evolution back to an age of the Universe of less than ~2 Gyrs, in the mid-IR rest-frame. This survey strategy will produce a full census of the Star formation Rate (SFR) in the Universe, using Polycyclic Aromatic Hydrocarbons (PAH) bands and fine-structure ionic lines, reaching the characteristic knee of the galaxy luminosity function, where the bulk of the population is distributed, at any redshift up to z ~3.5. Deep follow-up pointed spectroscopic observations with grating spectrometers { onboard the satellite}, across the full IR spectral range (17-210 micron), would simultaneously measure Black Hole Accretion Rate (BHAR), from high-ionization fine-structure lines, and SFR, from PAH and low- to mid-ionization lines in thousands of galaxies from solar to low metallicities, down to the knee of their luminosity functions. The analysis of the resulting atlas of IR spectra will reveal the physical processes at play in evolving galaxies across cosmic time, especially its heavily dust-embedded phase during the activity peak at the cosmic noon (z ~1-3), through IR emission lines and features that are insensitive to the dust obscuration., Comment: Paper accepted for publication on PASA on 5th March 2021, as part of the SPICA Special Issue