Researcher profile

G. Zamorani

G. Zamorani contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
60works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

60 published item(s)

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

preprint2026arXiv

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space.

preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid&#39;s Flagship simulation to investigate its impact on Euclid&#39;s weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

preprint2026arXiv

Euclid: Early Release Observations -- The extended stellar component of the IC10 dwarf galaxy

We present a detailed analysis of the old, extended stellar component of the Local Group dwarf galaxy IC 10 using deep resolved-star photometry in the VIS and NISP bands of the Euclid Early Release Observations. Leveraging Euclid&#39;s unique combination of wide field of view and high spatial resolution, we trace red giant branch (RGB) stars out to $\sim$8 kpc from the galaxy centre, reaching azimuthally-averaged surface brightness levels as faint as $μ_{HE}\sim$29 mag arcsec$^{-2}$. Our analysis reveals that IC 10&#39;s stellar distribution is significantly more extended than previously thought. After correcting for foreground extinction and subtracting contamination from Milky Way stars and background galaxies, we derive a radial stellar density profile from RGB star counts. The profile shows a marked flattening beyond $\sim$5 kpc, and is best fit by a two-component (Sersic + exponential) model, yielding a total stellar mass in old (age $>$1 Gyr) stars of $M_{\star}=(6.7-8.1)\times10^8 M_{\odot}$. The origin of the outer stellar component is unclear. It may be accreted, even possibly associated with the counter-rotating HI gas in the outer regions of IC 10, or it may represent an ancient in-situ stellar halo. We tentatively detect two symmetric stellar overdensities at the edge of our imagery. These roughly align with the direction of IC 10&#39;s orbit around M31, suggesting that they may be signatures of tidal stripping. As part of our analysis, we derive a new distance to IC 10 based on the RGB tip, finding $D=762\pm 20$ kpc and the distance modulus is $(m-M)_0=24.41\pm 0.05$.

preprint2026arXiv

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Weak lensing surveys require accurate correction for the point spread function (PSF) when measuring galaxy shapes. For a diffraction-limited PSF, as arises in space-based missions, this correction depends on each galaxy SED. In the Euclid mission, galaxy SED reconstruction, a tasks of the photometric-redshift processing function (PHZ PF), relies on broad- and medium-band ancillary photometry. The limited wavelength sampling of the Euclid VIS passband and signal-to-noise ratio may affect the reconstruction accuracy and translate into biases in the weak lensing measurements. In this study, we present the methodology, which is employed in the Euclid PHZ PF, for reconstructing galaxy SEDs at 55 wavelengths, sampling the VIS passband every 10 nm, and we assess whether it fulfils the accuracy requirements imposed on the Euclid PSF model. We employ both physics- and data-driven methods, focusing on a new approach of template-based flux correction and Gaussian processes, and we introduce an SED metric whose bias propagates into PSF quadrupole moment errors. Our findings demonstrate that Gaussian processes and template fitting meet the requirements only in specific, but complementary, redshift intervals. We therefore propose a hybrid approach, which leverages both methods. This solution proves to be effective in meeting the Euclid accuracy requirements for most of the redshift range of the survey. Finally, we investigate the impact on the SED reconstruction of a new set of 16 evenly-spaced medium-band filters for the Subaru telescope, providing quasi-spectroscopic coverage of the VIS passband. This study shows promising results, ensuring accurate SED reconstruction and meeting the mission PSF requirements. This work thus provides not only the methodological foundation of galaxy SED reconstruction in the Euclid PHZ PF, but also a roadmap for future improvements using a new medium-band survey.

preprint2025arXiv

Active galactic nuclei-heated dust revealed in &#34;little red dots&#34;

Little red dots (LRDs) are a puzzling population of extragalactic sources whose origin is highly debated. In this {work}, we performed a comprehensive stacking analysis of NIRCam, MIRI, and ALMA images of a large and homogeneously selected sample of LRDs from multiple JWST Legacy fields. We report clear evidence of hot-dust emission in the median stacked spectral energy distribution (SED) that features a rising near-infrared continuum up to rest-frame $λ_{\rm rest}$$\sim$ 3$μ$m, which is best explained by a standard dusty active galactic nucleus (AGN) structure. Although LRDs are likely to be a heterogeneous population, our findings suggest that most ($\gtrsim$50 %) LRDs show AGN-heated dust emission, regardless of whether the optical and ultraviolet (UV) continua are stellar or AGN-dominated. In either case, the best-fit dusty-AGN SED, combined with the lack of X-ray detection in the deep Chandra stacks, suggests that Compton-thick ($N_{\rm H}$$>$3$\times$10$^{24}$ cm$^{-2}$) gas obscuration is common, and likely confined within the dust sublimation radius ($R$$_{\rm sub}$$\sim$0.1 pc). Therefore, we argue that AGN-heated dust does not directly obscure either the optical-UV continuum or the broad-line region emission, in order to explain the observed blue UV slopes and prominent Balmer features. While a gas-dust displacement is in line with several models, the formation scenario (in-situ or ex-situ) of this pre-enriched hot dust remains unclear.

preprint2025arXiv

Euclid preparation. Simulating thousands of Euclid spectroscopic skies

We present two extensive sets of 3500+1000 simulations of dark matter haloes on the past light cone, and two corresponding sets of simulated (`mock&#39;) galaxy catalogues that represent the Euclid spectroscopic sample. The simulations were produced with the latest version of the PINOCCHIO code, and provide the largest, public set of simulated skies. Mock galaxy catalogues were obtained by populating haloes with galaxies using an halo occupation distribution (HOD) model extracted from the Flagship galaxy catalogue provided by Euclid Collaboration. The Geppetto set of 3500 simulated skies was obtained by tiling a 1.2 Gpc/h box to cover a light-cone whose sky footprint is a circle of 30 deg radius, for an area of 2763 deg$^2$ and a minimum halo mass of $1.5\times10^{11}$ Msun/h. The relatively small box size makes this set unfit for measuring very large scales. The EuclidLargeBox set consists of 1000 simulations of 3.38 Gpc/h, with the same mass resolution and a footprint that covers half of the sky, excluding the Milky Way zone of avoidance. From this we produced a set of 1000 EuclidLargeMocks on the 30 deg radius footprint, whose comoving volume is fully contained in the simulation box. We validated the two sets of catalogues by analysing number densities, power spectra, and 2-point correlation functions, showing that the Flagship spectroscopic catalogue is consistent with being one of the realisations of the simulated sets, although we noticed small deviations limited to the quadrupole at k>0.2 h/Mpc. We show cosmological parameter inference from these catalogues and demonstrate that using one realisation of EuclidLargeMocks in place of the Flagship mock produces the same posteriors, to within the expected shift given by sample variance. These simulated skies will be used for the galaxy clustering analysis of Euclid&#39;s Data Release 1 (DR1).

preprint2023arXiv

Euclid preparation: XXIII. Derivation of galaxy physical properties with deep machine learning using mock fluxes and H-band images

Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshifts, stellar masses, and star-formation rates (SFR) can be measured with deep learning algorithms for observed galaxies within data mimicking the Euclid and Rubin/LSST surveys. We find that Deep Learning Neural Networks and Convolutional Neutral Networks (CNN), which are dependent on the parameter space of the training sample, perform well in measuring the properties of these galaxies and have a better accuracy than methods based on spectral energy distribution fitting. CNNs allow the processing of multi-band magnitudes together with $H_{\scriptscriptstyle\rm E}$-band images. We find that the estimates of stellar masses improve with the use of an image, but those of redshift and SFR do not. Our best results are deriving i) the redshift within a normalised error of less than 0.15 for 99.9$\%$ of the galaxies with S/N>3 in the $H_{\scriptscriptstyle\rm E}$-band; ii) the stellar mass within a factor of two ($\sim0.3 \rm dex$) for 99.5$\%$ of the considered galaxies; iii) the SFR within a factor of two ($\sim0.3 \rm dex$) for $\sim$70$\%$ of the sample. We discuss the implications of our work for application to surveys as well as how measurements of these galaxy parameters can be improved with deep learning.

preprint2023arXiv

Euclid preparation. XXX. Performance assessment of the NISP Red-Grism through spectroscopic simulations for the Wide and Deep surveys

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$σ$ 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.

preprint2022arXiv

A combined VANDELS and LEGA-C study: the evolution of quiescent galaxy size, stellar mass and age from $\mathbf{\textit{z} = 0.6}$ to $\mathbf{\textit{z} = 1.3}$

We study the relationships between stellar mass, size and age within the quiescent population, using two mass-complete spectroscopic samples with $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}})>10.3$, taken from VANDELS at $1.0<z<1.3$, and LEGA-C at $0.6<z<0.8$. Using robust D$_{n}$4000 values, we demonstrate that the well-known &#39;downsizing&#39; signature is already in place by $z\simeq1.1$, with D$_{n}$4000 increasing by $\simeq0.1$ across a $\simeq$ 1 dex mass interval for both VANDELS and LEGA-C. We then proceed to investigate the evolution of the quiescent galaxy stellar mass-size relation from $z\simeq1.1$ to $z\simeq0.7$. We find the median size increases by a factor of $1.9\pm{0.1}$ at $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}})=10.5$, and see tentative evidence for flattening of the relation, finding slopes of $α=0.72\pm0.06$ and $α=$ $0.56\pm0.04$ for VANDELS and LEGA-C respectively. We finally split our sample into galaxies above and below our fitted mass-size relations, to investigate how size and D$_{n}$4000 correlate. For LEGA-C, we see a clear difference, with larger galaxies found to have smaller D$_{n}$4000 at fixed stellar mass. Due to the faintness and smaller numbers of the VANDELS sample, we cannot confirm whether a similar relation exists at $z\simeq1.1$. We consider whether differences in stellar age or metallicity are most likely to drive this size-D$_{n}$4000 relation, finding that any metallicity differences are unlikely to fully explain the observed offset, meaning smaller galaxies must be older than their larger counterparts. We find the observed evolution in size, mass and D$_{n}$4000 across the $\simeq2$ Gyr from $z\sim1.1$ to $z\sim0.7$ can be explained by a simple toy model in which VANDELS galaxies evolve passively, whilst experiencing a series of minor mergers.

preprint2022arXiv

A new estimate of the cosmic star formation density from a radio-selected sample, and the contribution of $H$-dark galaxies at $z \geq 3$

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&#39; 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 $ϕ$-integrated SFRD relative to the full sample at z $\sim 3.2$, and up to $58\%$ relative to LBG samples.

preprint2022arXiv

Euclid preparation: XIII. Forecasts for galaxy morphology with the Euclid Survey using Deep Generative Models

We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of $0.4\,\rm{deg}^2$ as it will be seen by the Euclid visible imager VIS and show that galaxy structural parameters are recovered with similar accuracy as for pure analytic Sérsic profiles. Based on these simulations, we estimate that the Euclid Wide Survey will be able to resolve the internal morphological structure of galaxies down to a surface brightness of $22.5\,\rm{mag}\,\rm{arcsec}^{-2}$, and $24.9\,\rm{mag}\,\rm{arcsec}^{-2}$ for the Euclid Deep Survey. This corresponds to approximately $250$ million galaxies at the end of the mission and a $50\,\%$ complete sample for stellar masses above $10^{10.6}\,\rm{M}_\odot$ (resp. $10^{9.6}\,\rm{M}_\odot$) at a redshift $z\sim0.5$ for the wide (resp. deep) survey. The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies.

preprint2022arXiv

Euclid Preparation: XIV. The Complete Calibration of the Color-Redshift Relation (C3R2) Survey: Data Release 3

The Complete Calibration of the Color-Redshift Relation (C3R2) survey is obtaining spectroscopic redshifts in order to map the relation between galaxy color and redshift to a depth of i ~ 24.5 (AB). The primary goal is to enable sufficiently accurate photometric redshifts for Stage IV dark energy projects, particularly Euclid and the Roman Space Telescope, which are designed to constrain cosmological parameters through weak lensing. We present 676 new high-confidence spectroscopic redshifts obtained by the C3R2 survey in the 2017B-2019B semesters using the DEIMOS, LRIS, and MOSFIRE multi-object spectrographs on the Keck telescopes. Combined with the 4454 redshifts previously published by this project, the C3R2 survey has now obtained and published 5130 high-quality galaxy spectra and redshifts. If we restrict consideration to only the 0.2 < z(phot) < 2.6 range of interest for the Euclid cosmological goals, then with the current data release C3R2 has increased the spectroscopic redshift coverage of the Euclid color space from 51% (as reported by Masters et al. 2015) to the current 91%. Once completed and combined with extensive data collected by other spectroscopic surveys, C3R2 should provide the spectroscopic calibration set needed to enable photometric redshifts to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for Roman.

preprint2022arXiv

Euclid preparation: XIX. Impact of magnification on photometric galaxy clustering

We investigate the importance of lensing magnification for estimates of galaxy clustering and its cross-correlation with shear for the photometric sample of Euclid. Using updated specifications, we study the impact of lensing magnification on the constraints and the shift in the estimation of the best fitting cosmological parameters that we expect if this effect is neglected. We follow the prescriptions of the official Euclid Fisher matrix forecast for the photometric galaxy clustering analysis and the combination of photometric clustering and cosmic shear. The slope of the luminosity function (local count slope), which regulates the amplitude of the lensing magnification, and the galaxy bias have been estimated from the Euclid Flagship simulation.We find that magnification significantly affects both the best-fit estimation of cosmological parameters and the constraints in the galaxy clustering analysis of the photometric sample. In particular, including magnification in the analysis reduces the 1$σ$ errors on $Ω_{\text{m},0}, w_{0}, w_a$ at the level of 20-35%, depending on how well we will be able to independently measure the local count slope. In addition, we find that neglecting magnification in the clustering analysis leads to shifts of up to 1.6$σ$ in the best-fit parameters. In the joint analysis of galaxy clustering, cosmic shear, and galaxy-galaxy lensing, magnification does not improve precision, but it leads to an up to 6$σ$ bias if neglected. Therefore, for all models considered in this work, magnification has to be included in the analysis of galaxy clustering and its cross-correlation with the shear signal ($3\times2$pt analysis) for an accurate parameter estimation.

preprint2022arXiv

Euclid preparation: XX. The Complete Calibration of the Color-Redshift Relation survey: LBT observations and data release

The Complete Calibration of the Color-Redshift Relation survey (C3R2) is a spectroscopic programme designed to empirically calibrate the galaxy color-redshift relation to the Euclid depth (I_E=24.5), a key ingredient for the success of Stage IV dark energy projects based on weak lensing cosmology. A spectroscopic calibration sample as representative as possible of the galaxies in the Euclid weak lensing sample is being collected, selecting galaxies from a self-organizing map (SOM) representation of the galaxy color space. Here, we present the results of a near-infrared H- and K-bands spectroscopic campaign carried out using the LUCI instruments at the LBT. For a total of 251 galaxies, we present new highly-reliable redshifts in the 1.3<= z <=1.7 and 2<= z<=2.7 ranges. The newly-determined redshifts populate 49 SOM cells which previously contained no spectroscopic measurements and almost double the occupation numbers of an additional 153 SOM cells. A final optical ground-based observational effort is needed to calibrate the missing cells in particular in the redshift range 1.7<= z<=2.7 that lack spectroscopic calibration. In the end, Euclid itself will deliver telluric-free NIR spectra that can complete the calibration.

preprint2022arXiv

Euclid preparation. XVIII. The NISP photometric system

Euclid will be the first space mission to survey most of the extragalactic sky in the 0.95-2.02 $μ$m range, to a 5$σ$ point-source median depth of 24.4 AB mag. This unique photometric data set will find wide use beyond Euclid&#39;s core science. In this paper, we present accurate computations of the Euclid Y_E, J_E and H_E passbands used by the Near-Infrared Spectrometer and Photometer (NISP), and the associated photometric system. We pay particular attention to passband variations in the field of view, accounting among others for spatially variable filter transmission, and variations of the angle of incidence on the filter substrate using optical ray tracing. The response curves&#39; cut-on and cut-off wavelengths - and their variation in the field of view - are determined with 0.8 nm accuracy, essential for the photometric redshift accuracy required by Euclid. After computing the photometric zeropoints in the AB mag system, we present linear transformations from and to common ground-based near-infrared photometric systems, for normal stars, red and brown dwarfs, and galaxies separately. A Python tool to compute accurate magnitudes for arbitrary passbands and spectral energy distributions is provided. We discuss various factors from space weathering to material outgassing that may slowly alter Euclid&#39;s spectral response. At the absolute flux scale, the Euclid in-flight calibration program connects the NISP photometric system to Hubble Space Telescope spectrophotometric white dwarf standards; at the relative flux scale, the chromatic evolution of the response is tracked at the milli-mag level. In this way, we establish an accurate photometric system that is fully controlled throughout Euclid&#39;s lifetime.

preprint2022arXiv

Euclid: Constraining ensemble photometric redshift distributions with stacked spectroscopy

The ESA Euclid mission will produce photometric galaxy samples over 15000 square degrees of the sky that will be rich for clustering and weak lensing statistics. The accuracy of the cosmological constraints derived from these measurements will depend on the knowledge of the underlying redshift distributions based on photometric redshift calibrations. A new approach is proposed to use the stacked spectra from Euclid slitless spectroscopy to augment broad-band photometric information to constrain the redshift distribution with spectral energy distribution fitting. The high spectral resolution available in the stacked spectra complements the photometry and helps to break the colour-redshift degeneracy and constrain the redshift distribution of galaxy samples. We modelled the stacked spectra as a linear mixture of spectral templates. The mixture may be inverted to infer the underlying redshift distribution using constrained regression algorithms. We demonstrate the method on simulated Vera C. Rubin Observatory and Euclid mock survey data sets based on the Euclid Flagship mock galaxy catalogue. We assess the accuracy of the reconstruction by considering the inference of the baryon acoustic scale from angular two-point correlation function measurements. We selected mock photometric galaxy samples at redshift z>1 using the self-organising map algorithm. Considering the idealised case without dust attenuation, we find that the redshift distributions of these samples can be recovered with 0.5% accuracy on the baryon acoustic scale. The estimates are not significantly degraded by the spectroscopic measurement noise due to the large sample size. However, the error degrades to 2% when the dust attenuation model is left free. We find that the colour degeneracies introduced by attenuation limit the accuracy considering the wavelength coverage of Euclid near-infrared spectroscopy.

preprint2022arXiv

Euclid: Covariance of weak lensing pseudo-$C_\ell$ estimates. Calculation, comparison to simulations, and dependence on survey geometry

An accurate covariance matrix is essential for obtaining reliable cosmological results when using a Gaussian likelihood. In this paper we study the covariance of pseudo-$C_\ell$ estimates of tomographic cosmic shear power spectra. Using two existing publicly available codes in combination, we calculate the full covariance matrix, including mode-coupling contributions arising from both partial sky coverage and non-linear structure growth. For three different sky masks, we compare the theoretical covariance matrix to that estimated from publicly available N-body weak lensing simulations, finding good agreement. We find that as a more extreme sky cut is applied, a corresponding increase in both Gaussian off-diagonal covariance and non-Gaussian super-sample covariance is observed in both theory and simulations, in accordance with expectations. Studying the different contributions to the covariance in detail, we find that the Gaussian covariance dominates along the main diagonal and the closest off-diagonals, but further away from the main diagonal the super-sample covariance is dominant. Forming mock constraints in parameters describing matter clustering and dark energy, we find that neglecting non-Gaussian contributions to the covariance can lead to underestimating the true size of confidence regions by up to 70 per cent. The dominant non-Gaussian covariance component is the super-sample covariance, but neglecting the smaller connected non-Gaussian covariance can still lead to the underestimation of uncertainties by 10--20 per cent. A real cosmological analysis will require marginalisation over many nuisance parameters, which will decrease the relative importance of all cosmological contributions to the covariance, so these values should be taken as upper limits on the importance of each component.

preprint2022arXiv

Euclid: Forecast constraints on consistency tests of the $Λ$CDM model

The standard cosmological model is based on the fundamental assumptions of a spatially homogeneous and isotropic universe on large scales. An observational detection of a violation of these assumptions at any redshift would immediately indicate the presence of new physics. We quantify the ability of the Euclid mission, together with contemporary surveys, to improve the current sensitivity of null tests of the canonical cosmological constant $Λ$ and the cold dark matter (LCDM) model in the redshift range $0<z<1.8$. We considered both currently available data and simulated Euclid and external data products based on a LCDM fiducial model, an evolving dark energy model assuming the Chevallier-Polarski-Linder (CPL) parameterization or an inhomogeneous Lema\^ıtre-Tolman-Bondi model with a cosmological constant $Λ$ (LLTB), and carried out two separate but complementary analyses: a machine learning reconstruction of the null tests based on genetic algorithms, and a theory-agnostic parametric approach based on Taylor expansion and binning of the data, in order to avoid assumptions about any particular model. We find that in combination with external probes, Euclid can improve current constraints on null tests of the LCDM by approximately a factor of three when using the machine learning approach and by a further factor of two in the case of the parametric approach. However, we also find that in certain cases, the parametric approach may be biased against or missing some features of models far from LCDM. Our analysis highlights the importance of synergies between Euclid and other surveys. These synergies are crucial for providing tighter constraints over an extended redshift range for a plethora of different consistency tests of some of the main assumptions of the current cosmological paradigm.

preprint2022arXiv

IR characteristic emission and dust properties of star-forming galaxies at 4.5 $<$ z $<$ 6.2

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

preprint2022arXiv

KiDS & Euclid: Cosmological implications of a pseudo angular power spectrum analysis of KiDS-1000 cosmic shear tomography

We present a tomographic weak lensing analysis of the Kilo Degree Survey Data Release 4 (KiDS-1000), using a new pseudo angular power spectrum estimator (pseudo-$C_{\ell}$) under development for the ESA Euclid mission. Over 21 million galaxies with shape information are divided into five tomographic redshift bins, ranging from 0.1 to 1.2 in photometric redshift. We measured pseudo-$C_{\ell}$ using eight bands in the multipole range $76<\ell<1500$ for auto- and cross-power spectra between the tomographic bins. A series of tests were carried out to check for systematic contamination from a variety of observational sources including stellar number density, variations in survey depth, and point spread function properties. While some marginal correlations with these systematic tracers were observed, there is no evidence of bias in the cosmological inference. B-mode power spectra are consistent with zero signal, with no significant residual contamination from E/B-mode leakage. We performed a Bayesian analysis of the pseudo-$C_{\ell}$ estimates by forward modelling the effects of the mask. Assuming a spatially flat $Λ$CDM cosmology, we constrained the structure growth parameter $S_8 = σ_8(Ω_{\rm m}/0.3)^{1/2} = 0.754_{-0.029}^{+0.027}$. When combining cosmic shear from KiDS-1000 with baryon acoustic oscillation and redshift space distortion data from recent Sloan Digital Sky Survey (SDSS) measurements of luminous red galaxies, as well as the Lyman-$α$ forest and its cross-correlation with quasars, we tightened these constraints to $S_8 = 0.771^{+0.006}_{-0.032}$. These results are in very good agreement with previous KiDS-1000 and SDSS analyses and confirm a $\sim 3σ$ tension with early-Universe constraints from cosmic microwave background experiments.

preprint2022arXiv

No strong dependence of Lyman continuum leakage on physical properties of star-forming galaxies at $\mathbf{3.1 \lesssim z \lesssim 3.5}$

We present Lyman continuum (LyC) radiation escape fraction $f_{\rm{esc}}$ measurements for 183 spectroscopically confirmed star-forming galaxies in the redshift range $3.11 < z < 3.53$ in the \textit{Chandra} Deep Field South. We use ground-based imaging to measure $f_{\rm{esc}}$, and use ground- and space-based photometry to derive galaxy physical properties using spectral energy distribution (SED) fitting. We additionally derive [O III]+H$β$ equivalent widths (that fall in the observed K band) by including nebular emission in the SED fitting. After removing foreground contaminants, we report the discovery of 11 new candidate LyC leakers, with absolute LyC escape fractions, $f_{\rm{esc}}$ in the range $0.14-0.85$. From non-detections, we place $1σ$ upper limits of $f_{\rm{esc}}<0.12$, where the Lyman-break selected galaxies have $f_{\rm{esc}} < 0.11$ and `blindly&#39; discovered galaxies with no prior photometric selection have $f_{\rm{esc}}<0.13$. We find a slightly higher $1σ$ limit of $f_{\rm{esc}}<0.20$ for extreme emission line galaxies with rest-frame [O III]+H$β$ equivalent widths $>300$A. For candidate LyC leakers, we find a weak negative correlation between $f_{\rm{esc}}$ and galaxy stellar masses, no correlation between $f_{\rm{esc}}$ specific star-formation rates (sSFRs) and a positive correlation between $f_{\rm{esc}}$ and EW$_0$([O III]+H$β$). The weak/no correlations between stellar mass and sSFRs may be explained by misaligned viewing angles and/or non-coincident timescales of starburst activity and periods of high $f_{\rm{esc}}$. Alternatively, escaping radiation may predominantly occur in highly localised star-forming regions, or $f_{\rm{esc}}$ measurements may be impacted by stochasticity of the intervening neutral medium, obscuring any global trends with galaxy properties. These hypotheses have important consequences for models of reionisation.

preprint2022arXiv

The ALPINE-ALMA [CII] Survey: Investigation of 10 Galaxies at $z\sim4.5$ with [OII] and [CII] Line Emission $-$ ISM Properties and [OII]$-$SFR Relation

We present $10$ main-sequence ALPINE galaxies (log($M/M_{\odot}$) = 9.2-11.1 and ${\rm SFR}=23-190\,{\rm M_{\odot}\,yr^{-1}}$) at $z\sim4.5$ with optical [OII] measurements from Keck/MOSFIRE spectroscopy and Subaru/MOIRCS narrow-band imaging. This is the largest such multi-wavelength sample at these redshifts, combining various measurements in the ultra-violet, optical, and far-infrared including [CII]$_{158{\rm μm}}$ line emission and dust continuum from ALMA and H$α$ emission from Spitzer photometry. For the first time, this unique sample allows us to analyze the relation between [OII] and total star-formation rate (SFR) and the interstellar medium (ISM) properties via [OII]/[CII] and [OII]/\halpha luminosity ratios at $z\sim4.5$. The [OII]$-$SFR relation at $z\sim4.5$ cannot be described using standard local descriptions, but is consistent with a metal-dependent relation assuming metallicities around $50\%$ solar. To explain the measured dust-corrected luminosity ratios of $L[OII]/L[CII] \sim 0.98^{+0.21}_{-0.22}$ and $L[OII]/LHa \sim -0.22^{+0.13}_{-0.15}$ for our sample, ionization parameters $\log(U)< -2$ and electron densities $\log(\rm n_e / {\rm [cm^{-3}]}) \sim 2.5-3$ are required. The former is consistent with galaxies at $z\sim2-3$, however lower than at $z>6$. The latter may be slightly higher than expected given the galaxies&#39; specific SFR. The analysis of this pilot sample suggests that typical log($ M/M_{\odot})$ > 9 galaxies at $z\sim4.5$ to have broadly similar ISM properties as their descendants at $z\sim2$ and suggest a strong evolution of ISM properties since the Epoch of Reionization at $z>6$.

preprint2022arXiv

The ALPINE-ALMA [CII] Survey: the population of [CII]-undetected galaxies and their role in the $\mathrm{L_{[CII]}}$-SFR relation

The [CII] 158$~μ$m emission line represents so far one of the most profitable tools for the investigation of the high-redshift galaxies in the early Universe. Being one of the brightest cooling lines in the rest-frame far-infrared regime of star-forming galaxies, it has been successfully exploited as a tracer of star-formation rate (SFR) in local sources. The picture is more complex at higher redshifts, where its usability in this context is still under investigation. Recent results from the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey suggest that there is no (or weak) evolution of the L$\mathrm{_{[CII]}}$-SFR relation up to $z\sim6$ but their reliability is hampered by the presence of a large population of [CII] non-detected galaxies. In this work, we characterize the population of [CII] non-detections in ALPINE. By stacking their ALMA spectra, we obtain a signal detected at $\sim5.1σ$, resulting in a [CII] luminosity of $\mathrm{log(L_\mathrm{[CII]}}/\mathrm{L_{\odot}})$ $\sim7.8$. When combining this value with those from the [CII] detections, we find a $\mathrm{L_{[CII]}}$-SFR relation with a slope $b=1.14\pm0.11$, in agreement within the uncertainties both with the linear relation found in the local Universe, and with the previous findings from ALPINE at $z\sim5$. This suggests that the [CII] line can be considered a good tracer of star formation up to the distant Universe. Finally, we show that the galaxies of our sample that most deviate from the observed L$_\mathrm{[CII]}$-SFR relation could suffer from a less precise redshift estimation, perhaps artificially reducing their [CII] luminosity. In this respect, we claim that there is no evidence in favour of a deficit of [CII] content in high-z galaxies, in contrast with earlier studies.

preprint2022arXiv

The stellar metallicities of massive quiescent galaxies at 1.0 < z < 1.3 from KMOS+VANDELS

We present a rest-frame UV-optical stacked spectrum representative of massive quiescent galaxies at $1.0<z<1.3$ with log$(M_*/\rm{M_\odot})>10.8$. The stack is constructed using VANDELS survey data, combined with new KMOS observations. We apply two independent full-spectral-fitting approaches, measuring a total metallicity, [Z/H]=$-0.13\pm0.08$ with Bagpipes, and [Z/H]=$0.04\pm0.14$ with Alf, a fall of $\sim0.2-0.3$ dex compared with the local Universe. We also measure an iron abundance, [Fe/H] =$-0.18\pm0.08$, a fall of $\sim0.15$ dex compared with the the local Universe. We measure the alpha enhancement via the magnesium abundance, obtaining [Mg/Fe]=$0.23\pm$0.12, consistent with similar-mass galaxies in the local Universe, indicating no evolution in the average alpha enhancement of log$(M_*/\rm{M_\odot})=11$ quiescent galaxies over the last $\sim8$ Gyr. This suggests the very high alpha enhancements recently reported for several bright $z\sim1-2$ quiescent galaxies are due to their extreme masses, log$(M_*/\rm{M_\odot})\gtrsim11.5$, rather than being typical of the $z\gtrsim1$ population. The metallicity evolution we observe with redshift (falling [Z/H], [Fe/H], constant [Mg/Fe]) is consistent with recent studies. We recover a mean stellar age of $2.5^{+0.6}_{-0.4}$ Gyr, corresponding to a formation redshift, $z_\rm{form}=2.4^{+0.6}_{-0.3}$. Recent studies have obtained varying average formation redshifts for $z\gtrsim1$ massive quiescent galaxies, and, as these studies report consistent metallicities, we identify different star-formation-history models as the most likely cause. Larger spectroscopic samples from upcoming ground-based instruments will provide precise constraints on ages and metallicities at $z\gtrsim1$. Combining these with precise JWST $z>2$ quiescent-galaxy stellar-mass functions will provide an independent test of formation redshifts derived from spectral fitting.

preprint2022arXiv

The VANDELS survey: a measurement of the average Lyman-continuum escape fraction of star-forming galaxies at z=3.5

We present a study designed to measure the average LyC escape fraction ($\langle f_{\rm esc}\rangle$) of star-forming galaxies at z=3.5. We assemble a sample of 148 galaxies from the VANDELS survey at $3.35\leq z_{\rm spec}\leq3.95$, selected to minimize line-of-sight contamination of their photometry. For this sample, we use ultra-deep, ground-based, $U-$band imaging and HST $V-$band imaging to robustly measure the distribution of $\mathcal{R_{\rm obs}}$ $=(L_{\rm LyC}/L_{\rm UV})_{\rm obs}$. We then model the distribution as a function of $\langle f_{\rm esc}\rangle$, carefully accounting for attenuation by dust, and the IGM (and CGM). A maximum likelihood fit to the $\mathcal{R_{\rm obs}}$ distribution returns a best-fitting value of $\langle f_{\rm esc}\rangle =0.07\pm0.02$, a result confirmed using an alternative Bayesian inference technique (both exclude $\langle f_{\rm esc}\rangle=0.0$ at $> 3σ$). By splitting our sample in two, we find evidence that $\langle f_{\rm esc}\rangle$ is positively correlated with Ly$α$ equivalent width, with high and low sub-samples returning best fits of $\langle f_{\rm esc}\rangle=0.12^{+0.06}_{-0.04}$ and $\langle f_{\rm esc} \rangle=0.02^{+0.02}_{-0.01}$, respectively. In contrast, we find evidence that $\langle f_{\rm esc}\rangle$ is anti-correlated with intrinsic UV luminosity and UV dust attenuation; with low UV luminosity and dust attenuation sub-samples returning best fits in the range $0.10 \leq \langle f_{\rm esc}\rangle \leq 0.22$. We do not find evidence for a clear correlation between $f_{\rm esc}$ and galaxy stellar mass, suggesting it is not a primary indicator of leakage. Although larger samples are needed to further explore these trends, they suggest that it is entirely plausible that the low dust and metallicity galaxies found at z > 6 will display the $\langle f_{\rm esc}\rangle\geq0.1$ required to drive reionization.

preprint2022arXiv

The VIMOS Ultra Deep Survey: The reversal of the star-formation rate $-$ density relation at $2 < z < 5$

Utilizing spectroscopic observations taken for the VIMOS Ultra-Deep Survey (VUDS), new observations from Keck/DEIMOS, and publicly available observations of large samples of star-forming galaxies, we report here on the relationship between the star formation rate (SFR) and the local environment ($δ_{gal}$) of galaxies in the early universe ($2<z<5$). Unlike what is observed at lower redshifts ($z<2$), we observe a definite, nearly monotonic increase in the average SFR with increasing galaxy overdensity over more than an order of magnitude in $δ_{gal}$. The robustness of this trend is quantified by accounting for both uncertainties in our measurements and galaxy populations that are either underrepresented or not present in our sample and find that the trend remains significant under all circumstances. This trend appears to be primarily driven by the fractional increase of galaxies in high density environments that are more massive in their stellar content and are forming stars at a higher rate than their less massive counterparts. We find that, even after stellar mass effects are accounted for, there remains a weak but significant SFR-$δ_{gal}$ trend in our sample implying that additional environmentally-related processes are helping to drive this trend. We also find clear evidence that the average SFR of galaxies in the densest environments increases with increasing redshift. These results lend themselves to a picture in which massive gas-rich galaxies coalesce into proto-cluster environments at $z\sim3$, interact with other galaxies or with a forming large-scale medium, subsequently using or losing most of their gas in the process, and begin to seed the nascent red sequence that is present in clusters at slightly lower redshifts.

preprint2021arXiv

$Euclid$ preparation: XV. Forecasting cosmological constraints for the $Euclid$ and CMB joint analysis

The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological parameters of the standard cosmological model and some of its extensions. This work expands and complements the recently published forecasts based on $Euclid$-specific probes, namely galaxy clustering, weak lensing, and their cross-correlation. With some assumptions on the specifications of current and future CMB experiments, the predicted constraints are obtained from both a standard Fisher formalism and a posterior-fitting approach based on actual CMB data. Compared to a $Euclid$-only analysis, the addition of CMB data leads to a substantial impact on constraints for all cosmological parameters of the standard $Λ$-cold-dark-matter model, with improvements reaching up to a factor of ten. For the parameters of extended models, which include a redshift-dependent dark energy equation of state, non-zero curvature, and a phenomenological modification of gravity, improvements can be of the order of two to three, reaching higher than ten in some cases. The results highlight the crucial importance for cosmological constraints of the combination and cross-correlation of $Euclid$ probes with CMB data.

preprint2021arXiv

An Optical Observational Cluster Mass Function at $z\sim1$ with the ORELSE Survey

We present a new mass function of galaxy clusters and groups using optical/near-infrared wavelength spectroscopic and photometric data from the Observations of Redshift Evolution in Large-Scale Environments (ORELSE) survey. At $z\sim1$, cluster mass function studies are rare regardless of wavelength and have never been attempted from an optical/near-infrared perspective. This work serves as a proof of concept that $z\sim1$ cluster mass functions are achievable without supplemental X-ray or Sunyaev-Zel&#39;dovich (SZ) data. Measurements of the cluster mass function provide important contraints on cosmological parameters and are complementary to other probes. With ORELSE, a new cluster finding technique based on Voronoi tessellation Monte-Carlo (VMC) mapping, and rigorous purity and completeness testing, we have obtained $\sim$240 galaxy overdensity candidates in the redshift range $0.55<z<1.37$ at a mass range of $13.6<\log(M/M_{\odot})<14.8$. This mass range is comparable to existing optical cluster mass function studies for the local universe. Our candidate numbers vary based on the choice of multiple input parameters related to detection and characterization in our cluster finding algorithm, which we incorporated into the mass function analysis through a Monte-Carlo scheme. We find cosmological constraints on the matter density of $Ω_{m} = 0.250^{+0.104}_{-0.099}$ and on the amplitude of fluctuations of $σ_{8} = 1.150^{+0.260}_{-0.163}$. While our $Ω_{m}$ value is close to concordance, our $σ_{8}$ value is $\sim2σ$ higher because of the inflated observed number densities compared to theoretical mass function models owing to how our survey targeted overdense regions. With Euclid and several other large, unbiased optical surveys on the horizon, VMC mapping will enable optical/NIR cluster cosmology at redshifts much higher than what has been possible before.

preprint2021arXiv

COSMOS2020: A panchromatic view of the Universe to $z\sim10$ from two complementary catalogs

The Cosmic Evolution Survey (COSMOS) has become a cornerstone of extragalactic astronomy. Since the last public catalog in 2015, a wealth of new imaging and spectroscopic data has been collected in the COSMOS field. This paper describes the collection, processing, and analysis of this new imaging data to produce a new reference photometric redshift catalog. Source detection and multi-wavelength photometry is performed for 1.7 million sources across the $2\,\mathrm{deg}^{2}$ of the COSMOS field, $\sim$966,000 of which are measured with all available broad-band data using both traditional aperture photometric methods and a new profile-fitting photometric extraction tool, The Farmer, which we have developed. A detailed comparison of the two resulting photometric catalogs is presented. Photometric redshifts are computed for all sources in each catalog utilizing two independent photometric redshift codes. Finally, a comparison is made between the performance of the photometric methodologies and of the redshift codes to demonstrate an exceptional degree of self-consistency in the resulting photometric redshifts. The $i<21$ sources have sub-percent photometric redshift accuracy and even the faintest sources at $25<i<27$ reach a precision of $5\,\%$. Finally, these results are discussed in the context of previous, current, and future surveys in the COSMOS field. Compared to COSMOS2015, reaches the same photometric redshift precision at almost one magnitude deeper. Both photometric catalogs and their photometric redshift solutions and physical parameters will be made available through the usual astronomical archive systems (ESO Phase 3, IPAC IRSA, and CDS).

preprint2021arXiv

Euclid preparation: I. The Euclid Wide Survey

Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$.

preprint2021arXiv

Euclid preparation: XVI. Exploring the ultra low-surface brightness Universe with Euclid/VIS

While Euclid is an ESA mission specifically designed to investigate the nature of Dark Energy and Dark Matter, the planned unprecedented combination of survey area ($\sim15\,000$ deg$^2$), spatial resolution, low sky-background, and depth also make Euclid an excellent space observatory for the study of the low surface brightness Universe. Scientific exploitation of the extended low surface brightness structures requires dedicated calibration procedures yet to be tested. We investigate the capabilities of Euclid to detect extended low surface brightness structure by identifying and quantifying sky background sources and stray-light contamination. We test the feasibility of generating sky flat-fields to reduce large-scale residual gradients in order to reveal the extended emission of galaxies observed in the Euclid Survey. We simulate a realistic set of Euclid/VIS observations, taking into account both instrumental and astronomical sources of contamination, including cosmic rays, stray-light, zodiacal light, ISM, and the CIB, while simulating the effects of the presence of background sources in the FOV. We demonstrate that a combination of calibration lamps, sky flats and self-calibration would enable recovery of emission at a limiting surface brightness magnitude of $μ=29.5^{+0.08}_{-0.27} $ mag arcsec$^{-2}$ ($3σ$, $10\times10$ arcsec$^2$) in the Wide Survey, reaching regions 2 magnitudes deeper in the Deep Surveys. Euclid/VIS has the potential to be an excellent low surface brightness observatory. Covering the gap between pixel-to-pixel calibration lamp flats and self-calibration observations for large scales, the application of sky flat-fielding will enhance the sensitivity of the VIS detector at scales of larger than 1 degree, up to the size of the FOV, enabling Euclid to detect extended surface brightness structures below $μ=31$ mag arcsec$^{-2}$ and beyond.

preprint2021arXiv

Euclid preparation: XVIII. Cosmic Dawn Survey. Spitzer observations of the Euclid deep fields and calibration fields

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer&#39;s Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s.

preprint2021arXiv

Euclid: Forecasts from redshift-space distortions and the Alcock-Paczynski test with cosmic voids

Euclid is poised to survey galaxies across a cosmological volume of unprecedented size, providing observations of more than a billion objects distributed over a third of the full sky. Approximately 20 million of these galaxies will have their spectroscopy available, allowing us to map the 3D large-scale structure of the Universe in great detail. This paper investigates prospects for the detection of cosmic voids therein and the unique benefit they provide for cosmology. In particular, we study the imprints of dynamic and geometric distortions of average void shapes and their constraining power on the growth of structure and cosmological distance ratios. To this end, we made use of the Flagship mock catalog, a state-of-the-art simulation of the data expected to be observed with Euclid. We arranged the data into four adjacent redshift bins, each of which contains about 11000 voids and estimated the stacked void-galaxy cross-correlation function in every bin. Fitting a linear-theory model to the data, we obtained constraints on $f/b$ and $D_M H$, where $f$ is the linear growth rate of density fluctuations, $b$ the galaxy bias, $D_M$ the comoving angular diameter distance, and $H$ the Hubble rate. In addition, we marginalized over two nuisance parameters included in our model to account for unknown systematic effects. With this approach, Euclid will be able to reach a relative precision of about 4% on measurements of $f/b$ and 0.5% on $D_M H$ in each redshift bin. Better modeling or calibration of the nuisance parameters may further increase this precision to 1% and 0.4%, respectively. Our results show that the exploitation of cosmic voids in Euclid will provide competitive constraints on cosmology even as a stand-alone probe. For example, the equation-of-state parameter $w$ for dark energy will be measured with a precision of about 10%, consistent with previous more approximate forecasts.

preprint2021arXiv

The Type II AGN-host galaxy connection: insights from the VVDS and VIPERS surveys

We present a study of optically-selected Type II AGN at 0.5 < z < 0.9 from the VIPERS and VVDS surveys, to investigate the connection between AGN activity and physical properties of their host galaxies. The host stellar mass is estimated through spectral energy distribution fitting with the CIGALE code, and star formation rates are derived from the [OII]$λ$3727 $Å$ line luminosity. We find that 49% of the AGN host galaxies are on or above the main sequence (MS), 40% lie in the sub-MS locus, and 11% in the quiescent locus. Using the [OIII]$λ$5007 $Å$ line luminosity as a proxy of the AGN power, we find that at fixed AGN power Type II AGN host galaxies show a bimodal behaviour: systems with host galaxy stellar mass <10$^{10}$ M$_{\odot}$, reside along the MS or in the starbursts locus (high-SF Type II AGN), while systems residing in massive host-galaxies (>10$^{10}$ M$_{\odot}$) show a lower level of star formation (low-SF Type II AGN). At all stellar masses, the offset from the MS is positively correlated with the AGN power. We interpret this correlation as evidence of co-evolution between the AGN and the host, possibly due to the availability of cold gas. In the most powerful AGN with host galaxies below the MS we find a hint, though weak, of asymmetry in the [OIII] line profile, likely due to outflowing gas, consistent with a scenario in which AGN feedback removes the available gas and halts the star formation in the most massive hosts.

preprint2021arXiv

The VANDELS survey: global properties of CIII]$λ$1908Å emitting star-forming galaxies at z$\sim$3

We study the mean properties of a large representative sample of 217 galaxies showing CIII] emission at $2<z<4$, selected from a parent sample of $\sim$750 main-sequence star-forming galaxies in the VANDELS survey. These CIII] emitters have a broad range of UV luminosities, thus allowing a detailed stacking analysis to characterize their stellar mass, star formation rate (SFR) and stellar metallicity, as a function of the UV emission line ratios, EWs, and the carbon-to-oxygen (C/O) abundance ratio. Reliable CIII] detections represent $\sim$30% of the parent sample. Extreme CIII] emitters (EW(CIII])$\gtrsim$8Å) are exceedingly rare ($\sim$3%) in VANDELS. The UV line ratios of the sample suggest no ionization source other than massive stars. Stacks with larger EW(CIII]) show larger EW(Ly$α$) and lower metallicity, but not all CIII] emitters are Ly$α$ emitters. The stellar metallicities of CIII] emitters are not significantly different from that of the parent sample, increasing from $\sim$10% to $\sim$40% solar for stellar masses $\log$(M$_{\star}$/M$_{\odot})\sim$9-10.5. The stellar mass-metallicity relation of the CIII] emitters is consistent with previous works showing strong evolution from $z=0$ to $z\sim3$. The C/O abundances of the sample range 35%-150% solar, with a noticeable increase with FUV luminosity and a smooth decrease with the CIII] EW. We discuss the CIII] emitters in the C/O-Fe/H and the C/O-O/H planes and find they follow stellar and nebular abundance trends consistent with those of Milky Way halo and thick disc stars and local HII galaxies, respectively. A qualitative agreement is also found with chemical evolution models, which suggests that CIII] emitters at $z\sim$3 are experiencing an active phase of chemical enrichment.

preprint2020arXiv

Dust and gas content of high-redshift galaxies hosting obscured AGN in the CDF-S

Obscured AGN represent a significant fraction of the entire AGN population, especially at high redshift (~70% at z=3--5). They are often characterized by the presence of large gas and dust reservoirs that are thought to sustain and possibly obscure vigorous star formation processes that make these objects shine at far-IR and sub-mm wavelengths. We exploit ALMA Cycle 4 observations of the continuum (~2.1mm) and high-J CO emission of a sample of six X-ray selected SMGs hosting an obscured AGN at z_spec>2.5 in the 7 Ms CDF-S. We measured the masses and sizes of the dust and molecular gas and we derived the gas density and column density on the basis of a uniform sphere geometry. Finally, we compared the measured column densities with those derived from the Chandra X-ray spectra. We detected both the continuum and line emission for three sources for which we measured both the flux density and size. For the undetected sources, we derived an upper limit on the flux density. We found that the detected galaxies are rich in gas and dust (molecular gas mass in the range <0.5 - 2.7 x 10^10 M_sun for α_CO=0.8 and up to ~2 x 10^11~M_sun for α_CO=6.5, and dust mass <0.9 - 4.9 x 10^8 M_sun) and compact (gas major axis 2.1-3.0 kpc, dust major axis 1.4-2.7 kpc). The column densities associated with the ISM are on the order of 10^(23-24) cm-2, which is comparable with those derived from the X-ray spectra. For the detected sources we also derived dynamical masses in the range 0.8 - 3.7 x 10^10 M_sun. We conclude that the ISM of high redshift galaxies can substantially contribute to nuclear obscuration up to the Compton-thick (>10^24 cm-2) regime. In addition, we found that all the detected sources show a velocity gradient reminding one rotating system, even though two of them show peculiar features in their morphology that can be associated with a chaotic, possibly merging, structure.

preprint2020arXiv

Euclid preparation: VIII. The Complete Calibration of the Colour-Redshift Relation survey: VLT/KMOS observations and data release

The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1.&#34;2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

preprint2020arXiv

Euclid: The selection of quiescent and star-forming galaxies using observed colours

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ&#39; diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

preprint2020arXiv

Modelling the quenching of star formation activity from the evolution of the colour-magnitude relation in VIPERS

We study the evolution of the colour-magnitude relation for galaxies in the VIMOS Public Extragalactic Redshift Survey (VIPERS) by introducing the concept of the bright edge, and use this to derive constraints on the quenching of star formation activity in galaxies over the redshift range $0.5 < z < 1.1$. The bright-edge of the colour-magnitude diagram evolves with little dependence on galaxy colour, and therefore on the amount of star formation taking place in bright galaxies. We modelled this evolution with delayed exponential star formation histories (SFHs), to better understand the time-scale of the turn-off in star formation activity. We show that using SFHs without quenching, the transition from the blue cloud to the red sequence is too slow. This indicates that a scenario purely driven by the consumption of the gas inside each galaxy does not reproduce the observed evolution of the colour-magnitude bright edge. Among the quenching scenarios explored, the one that best matches the observations assumes that galaxies stop their star formation at a randomly selected time with a uniform distribution up to $2.5$ Gyr. We argue that quenching is required over a wide range of stellar masses. Qualitatively similar evolution of the bright edge is found in the predictions of a semi-analytical galaxy formation model, but quantitatively there are marked differences with the observations. This illustrates the utility of the bright edge as a test of galaxy formation models. The evolution changes and no longer matches the observed trend if feedback from heating by active galactic nuclei is turned off.

preprint2020arXiv

MOONRISE: The Main MOONS GTO Extragalactic Survey

The MOONS instrument possesses an exceptional combination of large multiplexing, high sensitivity, broad simultaneous spectral coverage (from optical to near-infrared bands), large patrol area and high fibre density. These properties provide the unprecedented potential of enabling, for the very first time, SDSS-like surveys around Cosmic Noon (z~1-2.5), when the star formation rate in the Universe peaked. The high-quality spectra delivered by MOONS will sample the same nebular and stellar diagnostics observed in extensive surveys of local galaxies, providing an accurate and consistent description of the evolution of various physical properties of galaxies, and hence a solid test of different scenarios of galaxy formation and transformation. Most importantly, by spectroscopically identifying hundreds of thousands of galaxies at high redshift, the MOONS surveys will be capable of determining the environments in which primeval galaxies lived and will reveal how such environments affected galaxy evolution. In this article, we specifically focus on the main Guaranteed Time Observation (GTO) MOONS extragalactic survey, MOONRISE, by providing an overview of its scientific goals and observing strategy.

preprint2020arXiv

SUPER III. Broad Line Region properties of AGN at z$\sim$2

The SINFONI survey for Unveiling the Physics and Effect of Radiative feedback (SUPER) was designed to conduct a blind search for AGN-driven outflows on X-ray selected AGN at redshift z$\sim$2 with high ($\sim$2 kpc) spatial resolution, and correlate them to the properties of the host galaxy and central black hole. The main aims of this paper are: a) to derive reliable estimates for the BH mass and accretion rates for the Type-1 AGN in this survey; b) to characterize the properties of the AGN driven winds in the BLR. We analyzed rest-frame optical and UV spectra of 21 Type-1 AGN. We found that the BH masses estimated from H$α$ and H$β$ lines are in agreement. We estimate BH masses in the range Log(M$\rm_{BH}/M_{\odot}$)=8.4-10.8 and Eddington ratios $\rmλ_{Edd}$ =0.04-1.3. We confirm that the CIV line width does not correlate with the Balmer lines and the peak of the line profile is blue-shifted with respect to the [OIII]-based systemic redshift. These findings support the idea that the CIV line is tracing outflowing gas in the BLR, with velocities up to $\sim$4700 km/s. We confirm the strong dependence of the BLR wind velocity with the UV-to-Xray continuum slope, L$\rm_{Bol}$ and $\rmλ_{Edd}$. We inferred BLR mass outflow rates in the range 0.005-3 M$_{\odot}$/yr, showing a correlation with the bolometric luminosity consistent with that observed for ionized winds in the NLR and X-ray winds detected in local AGN, and kinetic power $\sim$10$^{[-7:-4]}\times$ L$\rm_{Bol}$. Finally, we found an anti-correlation between the equivalent width of the [OIII] line with respect to the CIV shift, and a positive correlation with [OIII] outflow velocity. These findings, for the first time in an unbiased sample of AGN at z$\sim$2, support a scenario where BLR winds are connected to galaxy scale detected outflows, and are capable of affecting the gas in the NLR located at kpc scale.

preprint2020arXiv

The ALPINE-ALMA [C II] Survey: [C II]158micron Emission Line Luminosity Functions at $z \sim 4-6$

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

preprint2020arXiv

The ALPINE-ALMA [CII] Survey: Multi-Wavelength Ancillary Data and Basic Physical Measurements

We present the ancillary data and basic physical measurements for the galaxies in the ALMA Large Program to Investigate C+ at Early Times (ALPINE) survey - the first large multi-wavelength survey which aims at characterizing the gas and dust properties of 118 main-sequence galaxies at redshifts 4.4 < z < 5.9 via the measurement of [CII]-emission at 158 micro-meter (64% at >3.5$σ$) and the surrounding far-infrared (FIR) continuum in conjunction with a wealth of optical and near-infrared data. We outline in detail the spectroscopic data and selection of the galaxies as well as the ground- and space-based imaging products. In addition, we provide several basic measurements including stellar masses, star formation rates (SFR), rest-frame ultra-violet (UV) luminosities, UV continuum slopes ($β$), and absorption line redshifts, as well as H-alpha emission derived from Spitzer colors. We find that the ALPINE sample is representative of the 4 < z < 6 galaxy population selected by photometric methods and only slightly biased towards bluer colors ($Δβ$ ~ 0.2). Using [CII] as tracer of the systemic redshift (confirmed for one galaxy at z=4.5 out of 118 for which we obtained optical [OII]$λ$3727A emission), we confirm red shifted Ly-alpha emission and blue shifted absorption lines similar to findings at lower redshifts. By stacking the rest-frame UV spectra in the [CII] rest-frame we find that the absorption lines in galaxies with high specific SFR are more blue shifted, which could be indicative of stronger winds and outflows.

preprint2020arXiv

The ALPINE-ALMA [CII] Survey: On the nature of an extremely obscured serendipitous galaxy

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

preprint2020arXiv

The ALPINE-ALMA [CII] Survey: Size of Individual Star-Forming Galaxies at z=4-6 and their Extended Halo Structure

We present the physical extent of [CII] 158um line-emitting gas from 46 star-forming galaxies at z=4-6 from the ALMA Large Program to INvestigate CII at Early Times (ALPINE). Using exponential profile fits, we measure the effective radius of the [CII] line (r_e,[CII]) for individual galaxies and compare them with the rest-frame ultra-violet (UV) continuum (r_e,UV) from Hubble Space Telescope images. The effective radius r_e,[CII] exceeds r_e,UV by factors of ~2-3 and the ratio of r_e,[CII]/r_e,UV increases as a function of M_star. We do not find strong evidence that [CII] line, the rest-frame UV, and FIR continuum are always displaced over ~ 1-kpc scale from each other. We identify 30% of isolated ALPINE sources as having an extended [CII] component over 10-kpc scales detected at 4.1$σ$-10.9$σ$ beyond the size of rest-frame UV and far-infrared (FIR) continuum. One object has tentative rotating features up to ~10-kpc, where the 3D model fit shows the rotating [CII]-gas disk spread over 4 times larger than the rest-frame UV-emitting region. Galaxies with the extended [CII] line structure have high star-formation rate (SFR), stellar mass (M_star), low Lya equivalent-width, and more blue-shifted (red-shifted) rest-frame UV metal absorption (Lya line), as compared to galaxies without such extended [CII] structures. Although we cannot rule out the possibility that a selection bias towards luminous objects may be responsible for such trends, the star-formation driven outflow also explains all these trends. Deeper observations are essential to test whether the extended [CII] line structures are ubiquitous to high-z star-forming galaxies.

preprint2020arXiv

The deep Chandra survey in the SDSS J1030+0524 field

We present the X-ray source catalog for the 479 ks Chandra exposure of the SDSS J1030+0524 field, that is centered on a region that shows the best evidence to date of an overdensity around a z > 6 quasar, and also includes a galaxy overdensity around a Compton-thick Fanaroff-Riley type II radio galaxy at z = 1.7. Using wavdetect for initial source detection and ACIS Extract for source photometry and significance assessment, we create preliminary catalogs of sources that are detected in the full, soft, and hard bands, respectively. We produce X-ray simulations that mirror our Chandra observation to filter our preliminary catalogs and get a completeness level of > 91% and a reliability level of 95% in each band. The catalogs in the three bands are then matched into a final main catalog of 256 unique sources. Among them, 244, 193, and 208 are detected in the full, soft, and hard bands, respectively. This makes J1030 field the fifth deepest extragalactic X-ray survey to date. The field is part of the Multiwavelength Survey by Yale-Chile (MUSYC), and is also covered by optical imaging data from the Large Binocular Camera (LBC) at the Large Binocular Telescope, near-IR imaging data from the Canada France Hawaii Telescope WIRCam, and Spitzer IRAC. Thanks to its dense multi-wavelength coverage, J1030 represents a legacy field for the study of large-scale structures around distant accreting supermassive black holes. Using a likelihood ratio analysis, we associate multi-band counterparts for 252 (98.4%) of the 256 Chandra sources, with an estimated reliability of 95%. Finally, we compute the cumulative number of sources in each X-ray band, finding that they are in general agreement with the results from the Chandra Deep Fields.

preprint2020arXiv

The evolving AGN duty cycle in galaxies since z$\sim$3 as encoded in the X-ray luminosity function

We present a new modeling of the X-ray luminosity function (XLF) of Active Galactic Nuclei (AGN) out to z$\sim$3, dissecting the contribution of main-sequence (MS) and starburst (SB) galaxies. For each galaxy population, we convolved the observed galaxy stellar mass (M$_{\star}$) function with a grid of M$_{\star}$-independent Eddington ratio ($λ_{\rm EDD}$) distributions, normalised via empirical black hole accretion rate (BHAR) to star formation rate (SFR) relations. Our simple approach yields an excellent agreement with the observed XLF since z$\sim$3. We find that the redshift evolution of the observed XLF can only be reproduced through an intrinsic flattening of the $λ_{\rm EDD}$ distribution, and with a positive shift of the break $λ^{*}$, consistent with an anti-hierarchical behavior. The AGN accretion history is predominantly made by massive (10$^{10}<$M$_{\star}<$10$^{11}$ M$_{\odot}$) MS galaxies, while SB-driven BH accretion, possibly associated with galaxy mergers, becomes dominant only in bright quasars, at $\log$(L$_{\rm X}$/erg s$^{-1}$)$>$44.36 + 1.28$\cdot$(1+z). We infer that the probability of finding highly-accreting ($λ_{\rm EDD}>$ 10%) AGN significantly increases with redshift, from 0.4% (3.0%) at z=0.5 to 6.5% (15.3%) at z=3 for MS (SB) galaxies, implying a longer AGN duty cycle in the early Universe. Our results strongly favor a M$_{\star}$-dependent ratio between BHAR and SFR, as BHAR/SFR $\propto$ M$_{\star}^{0.73[+0.22,-0.29]}$, supporting a non-linear BH buildup relative to the host. Finally, this framework opens potential questions on super-Eddington BH accretion and different $λ_{\rm EDD}$ prescriptions for understanding the cosmic BH mass assembly.

preprint2020arXiv

The properties of He II 1640 emitters at z ~ 2.5-5 from the VANDELS survey

Strong He II emission is produced by low-metallicity stellar populations. Here, we aim to identify and study a sample of He II $λ1640$-emitting galaxies at redshifts of $z \sim 2.5-5$ in the deep VANDELS spectroscopic survey.. We identified a total of 33 Bright He II emitters (S/N > 2.5) and 17 Faint emitters (S/N < 2.5) in the VANDELS survey and used the available deep multi-wavelength data to study their physical properties. After identifying seven potential AGNs in our sample and discarding them from further analysis, we divided the sample of \emph{Bright} emitters into 20 \emph{Narrow} (FWHM < 1000 km s$^{-1}$) and 6 \emph{Broad} (FWHM > 1000 km s$^{-1}$) He II emitters. We created stacks of Faint, Narrow, and Broad emitters and measured other rest-frame UV lines such as O III] and C III] in both individual galaxies and stacks. We then compared the UV line ratios with the output of stellar population-synthesis models to study the ionising properties of He II emitters. We do not see a significant difference between the stellar masses, star-formation rates, and rest-frame UV magnitudes of galaxies with He II and no He II emission. The stellar population models reproduce the observed UV line ratios from metals in a consistent manner, however they under-predict the total number of \heii ionising photons, confirming earlier studies and suggesting that additional mechanisms capable of producing He II are needed, such as X-ray binaries or stripped stars. The models favour subsolar metallicities ($\sim0.1Z_\odot$) and young stellar ages ($10^6 - 10^7$ years) for the He II emitters. However, the metallicity measured for He II emitters is comparable to that of non-He II emitters at similar redshifts. We argue that galaxies with He II emission may have undergone a recent star-formation event, or may be powered by additional sources of He II ionisation.

preprint2020arXiv

The role of galaxy mass on AGN emission: a view from the VANDELS survey

We present a comparative analysis of the properties of AGN emitting at radio and X-ray wavelengths. The study is performed on 907 X-ray AGN and 100 radio AGN selected on the CDFS and UDS fields and makes use of new and ancillary data available to the VANDELS collaboration. Our results indicate that the mass of the host galaxy is a fundamental quantity which determines the level of AGN activity at the various wavelengths. Indeed large stellar masses are found to be connected with AGN radio emission, as virtually all radio-active AGN reside within galaxies of M*>10^{10} Msun. Large stellar masses also seem to favour AGN activity in the X-ray, even though X-ray AGN present a mass distribution which is more spread out and with a non-negligible tail at M*<10^{9} Msun. Stellar mass alone is also observed to play a fundamental role in simultaneous radio and X-ray emission: the percentage of AGN active at both wavelengths increases from around 1% of all X-ray AGN residing within hosts of M*<10^{11} Msun to about 13% in more massive galaxies. In the case of radio-selected AGN, such a percentage moves from about 15% to about 45% (but up to 80% in the deepest fields). Neither cosmic epoch, nor radio luminosity, X-ray luminosity, Eddington ratio or star-formation rate of the hosts are found to be connected to an enhanced probability for joint radio+X-ray emission of AGN origin. Furthermore, only a loose relation is observed between X-ray and radio luminosity in those AGN which are simultaneously active at both frequencies.

preprint2020arXiv

The VANDELS survey: A strong correlation between Ly$α$ equivalent width and stellar metallicity at $\mathbf{3\leq z \leq 5}$

We present the results of a new study investigating the relationship between observed Ly$α$ equivalent width ($W_λ$(Ly$α$)) and the metallicity of the ionizing stellar population ($Z_{\star}$) for a sample of $768$ star-forming galaxies at $3 \leq z \leq 5$ drawn from the VANDELS survey. Dividing our sample into quartiles of rest-frame $W_λ$(Ly$α$) across the range $-58 \unicode{xC5} \lesssim$ $W_λ$(Ly$α$) $\lesssim 110 \unicode{xC5}$ we determine $Z_{\star}$ from full spectral fitting of composite far-ultraviolet (FUV) spectra and find a clear anti-correlation between $W_λ$(Ly$α$) and $Z_{\star}$. Our results indicate that $Z_{\star}$ decreases by a factor $\gtrsim 3$ between the lowest $W_λ$(Ly$α$) quartile ($\langle$$W_λ$(Ly$α$)$\rangle=-18\unicode{xC5}$) and the highest $W_λ$(Ly$α$) quartile ($\langle$$W_λ$(Ly$α$)$\rangle=24\unicode{xC5}$). Similarly, galaxies typically defined as Lyman Alpha Emitters (LAEs; $W_λ$(Ly$α$) $>20\unicode{xC5}$) are, on average, metal poor with respect to the non-LAE galaxy population ($W_λ$(Ly$α$) $\leq20\unicode{xC5}$) with $Z_{\star}$$_{\rm{non-LAE}}\gtrsim 2 \times$ $Z_{\star}$$_{\rm{LAE}}$. Finally, based on the best-fitting stellar models, we estimate that the increasing strength of the stellar ionizing spectrum towards lower $Z_{\star}$ is responsible for $\simeq 15-25\%$ of the observed variation in $W_λ$(Ly$α$) across our sample, with the remaining contribution ($\simeq 75-85\%$) being due to a decrease in the HI/dust covering fractions in low $Z_{\star}$ galaxies.

preprint2020arXiv

The VANDELS survey: Discovery of massive overdensities of galaxies at z>2

We want to investigate whether we can use Lyalpha emission to obtain information on the environment properties and whether Lyalpha emitters show different characteristics as a function of their environment. We estimated local densities in the VANDELS Chandra Deep Field-South (CDFS) and UKIDSS Ultra Deep Survey (UDS) fields, by using a three-dimensional algorithm which works in the RA-dec-redshift space. We selected a sample of 131 Lyalpha-emitting galaxies (EW(Lyalpha)>0 A), unbiased with respect to environmental density, to study their location with respect to the over- or under-dense environment. We identify 13 (proto)cluster candidates in the CDFS and nine in the UDS at 2<z<4, based on photometric and spectroscopic redshifts from VANDELS and from all the available literature. No significant difference is observed in the rest-frame U-V color between field and galaxies located within the identified overdensities. We find that VANDELS Lyalpha emitters (LAEVs) lie preferentially outside of overdense regions as the majority of the galaxies with spectroscopic redshifts from VANDELS. The LAEVs in overdense regions tend to have low Lyalpha equivalent widths and low specific SFRs, and they also tend to be more massive than the LAEVs in the field. Their stacked Lyalpha profile shows a dominant red peak and a hint of a blue peak. Our results show that LAEVs are likely to be influenced by the environment and favour a scenario with outflows of low expansion velocities and high HI column densities for galaxies in overdense regions. An outflow with low expansion velocity could be related to the way galaxies are forming stars in overdense regions; the high HI column density can be a consequence of the gravitational potential of the overdensity. Therefore, Lyalpha-emitting galaxies can provide useful insights on the environment in which they reside.

preprint2020arXiv

The VANDELS survey: the relation between UV continuum slope and stellar metallicity in star-forming galaxies at z~3

The estimate of stellar metallicities (Z*) of high-z galaxies are of paramount importance in order to understand the complexity of dust effects and the reciprocal interrelations among stellar mass, dust attenuation, stellar age, and metallicity. Benefiting from uniquely deep FUV spectra of >500 star-forming galaxies at redshifts 2<z<5 extracted from the VANDELS survey and stacked in bins of stellar mass (M*) and UV continuum slope (beta), we estimate their stellar metallicities Z* from stellar photospheric absorption features at 1501 and 1719 Angstrom, which are calibrated with Starburst99 models and are largely unaffected by stellar age, dust, IMF, nebular continuum or interstellar absorption. Comparing them to photometric based spectral slopes in the range 1250-1750 Angstrom, we find that the stellar metallicity increases by ~0.5 dex from beta ~ -2 to beta ~ -1 (1 < A(1600) < 3.2), and a dependence with beta holds at fixed UV absolute luminosity M(UV) and stellar mass up to 10^(9.65) Msun. As a result, the metallicity is a fundamental ingredient for properly rescaling dust corrections based on M(UV) and M*. Using the same absorption features, we analyze the mass-metallicity relation (MZR) and find it is consistent with the previous VANDELS estimation based on a global fit of the FUV spectra. Similarly, we do not find a significant evolution between z=2 and z=3.5. Finally, the slopes of our MZR and Z*-beta relation are in agreement with the predictions of well-studied semi-analytic models of galaxy formation (SAM), while some tensions with observations remain as to the absolute metallicity normalization. The relation between UV slope and stellar metallicity is fundamental for the exploitation of large volume surveys with next-generation telescopes and for the physical characterization of galaxies in the first billion years of our Universe.

preprint2020arXiv

The VIMOS Ultra-Deep Survey: the Ly$α$ emission line morphology at $2 < z < 6$

The Lyman-$α$ (Ly$α$) emission line has been ubiquitously used to confirm and study high redshift galaxies. We report on the line morphology as seen in the 2D spectra from the VIMOS Ultra Deep Survey in a sample of 914 Ly$α$ emitters from a parent sample of 4192 star-forming galaxies at $2<z_\mathrm{spec}\lesssim6$. The study of the spatial extent of Ly$α$ emission provides insight into the escape of Ly$α$ photons from galaxies. We classify the line emission as either non-existent, coincident, projected spatial offset, or extended with respect to the observed 2D UV continuum emission. The line emitters in our sample are classified as ~45% coincident, ~24% extended and ~11% offset emitters. For galaxies with detected UV continuum, we show that extended Ly$α$ emitters (LAEs) correspond to the highest equivalent width galaxies (with an average $W_\mathrm{Lyα}\sim-22$Å). This means that this class of objects is the most common in narrow-band selected samples, which usually select high equivalent width LAEs, $<-20$Å. Extended Ly$α$ emitters are found to be less massive, less star-forming, with lower dust content, and smaller UV continuum sizes ($r_{50}\sim0.9$kpc) of all the classes considered here. We also find that galaxies with larger UV-sizes have lower fractions of Ly$α$ emitters. By stacking the spectra per emitter class we find that the weaker Ly$α$ emitters have stronger low ionization inter-stellar medium (ISM) absorption lines. Interestingly, we find that galaxies with Ly$α$ offset emission (median separation of $1.1_{-0.8}^{+1.3}$kpc from UV continuum) show similar velocity offsets in the ISM as those with no visible emission (and different from other Ly$α$ emitting classes). This class of objects may hint at episodes of gas accretion, bright offset clumps, or on-going merging activity into the larger galaxies.

preprint2020arXiv

Timing the earliest quenching events with a robust sample of massive quiescent galaxies at 2 < z < 5

We present a sample of 151 massive ($M_* > 10^{10}\mathrm{M_\odot}$) quiescent galaxies at $2 < z < 5$, based on a sophisticated Bayesian spectral energy distribution fitting analysis of the CANDELS UDS and GOODS-South fields. Our sample includes a robust sub-sample of 61 objects for which we confidently exclude low-redshift and star-forming solutions. We identify 10 robust objects at $z>3$, of which 2 are at $z>4$. We report formation redshifts, demonstrating that the oldest objects formed at $z > 6$, however individual ages from our photometric data have significant uncertainties, typically $\sim0.5$ Gyr. We demonstrate that the UVJ colours of the quiescent population evolve with redshift at $z>3$, becoming bluer and more similar to post-starburst galaxies at lower redshift. Based upon this we construct a model for the time-evolution of quiescent galaxy UVJ colours, concluding that the oldest objects are consistent with forming the bulk of their stellar mass at $z\sim6-7$ and quenching at $z\sim5$. We report spectroscopic redshifts for two of our objects at $z=3.440$ and $3.396$, which exhibit extremely weak Ly$α$ emission in ultra-deep VANDELS spectra. We calculate star-formation rates based on these line fluxes, finding that these galaxies are consistent with our quiescent selection criteria, provided their Ly$α$ escape fractions are $>3$ and $>10$ per cent respectively. We finally report that our highest-redshift robust object exhibits a continuum break at $λ\sim7000$A in a spectrum from VUDS, consistent with our photometric redshift of $z_\mathrm{phot}=4.72^{+0.06}_{-0.04}$. If confirmed as quiescent this object would be the highest-redshift known quiescent galaxy. To obtain stronger constraints on the times of the earliest quenching events, high-SNR spectroscopy must be extended to $z\gtrsim3$ quiescent objects.

preprint2020arXiv

X-ray properties of He II 1640 emitting galaxies in VANDELS

We explore X-ray emission from a sample of 18 He II 1640 emitting star-forming galaxies at z ~ 2.3-3.6 from the VANDELS survey in the Chandra Deep Field South, to set constraints on the role of X-ray sources in powering the He II emission. We find that 4 He II emitters have tentative detections with S/N ~ 2 and have X-ray luminosities, L_X = 1.5-4.9 x 10^41 erg/s. The stacked luminosity of all 18 He II emitters is 2.6 x 10^41 erg/s, and that of a subset of 13 narrow He II emitters (FHWM(He II) < 1000 km/s) is 3.1 x 10^41 erg/s. We also measure stacked L_X for non-He II emitters through bootstrapping of matched samples, and find L_X = 2.5 x 10^41 erg/s, which is not significantly different from L_X measured for He II emitters. The L_X per star-formation rate for He II emitters (log (L_X/SFR) ~ 40.0) and non-emitters (log (L_X/SFR) ~ 39.9) are also comparable and in line with the redshift evolution and metallicity dependence predicted by models. Due to the non-significant difference between the X-ray emission from galaxies with and without He II, we conclude that X-ray binaries or weak or obscured AGNs are unlikely to be the dominant producers of He II ionising photons in VANDELS star-forming galaxies at z ~ 3. Given the comparable physical properties of both He II emitters and non-emitters reported previously, alternative He II ionising mechanisms such as localised low-metallicity stellar populations, Pop-III stars, etc. may need to be explored.

preprint2019arXiv

The ALPINE-ALMA [CII] survey: Star-formation-driven outflows and circumgalactic enrichment in the early Universe

We study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 um emission in a large sample of normal star-forming galaxies at 4 < z < 6 from the ALMA Large Program to INvestigate [C II] at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked [C II] profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked [C II] spectrum, with velocities of |v|<~ 500 km/s. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFRmed = 25 Msun/yr), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at z > 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v|< 200 km/s) of the higher-SFR galaxies is extended on physical sizes of ~ 30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.

preprint2019arXiv

The galaxy&#39;s gas content regulated by the dark matter halo mass results in a super-linear M$_{\rm BH}$-M$_{\star}$ relation

Supermassive black holes (SMBHs) are tightly correlated with their hosts but the origin of such connection remains elusive. To explore the cosmic build-up of this scaling relation, we present an empirically-motivated model that tracks galaxy and SMBH growth down to z=0. Starting from a random mass seed distribution at z=10, we assume that each galaxy evolves on the star-forming &#34;main sequence&#34; (MS) and each BH follows the recently-derived stellar mass (M$_{\star}$) dependent ratio between BH accretion rate and star formation rate, going as BHAR/SFR$\propto$M$_{\star}^{0.73[+0.22,-0.29]}$. Our simple recipe naturally describes the BH-galaxy build-up in two stages. At first, the SMBH lags behind the host that evolves along the MS. Later, as the galaxy grows in M$_{\star}$, our M$_{\star}$-dependent BHAR/SFR induces a super-linear BH growth, as M$_{\rm BH}$$\propto$M$_{\star}^{1.7}$. According to this formalism, smaller BH seeds increase their relative mass faster and earlier than bigger BH seeds, at fixed M$_{\star}$, thus setting along a gradually tighter M$_{\rm BH}$-M$_{\star}$ locus towards higher M$_{\star}$. Assuming reasonable values of the radiative efficiency $ε\sim$0.1, our empirical trend agrees with both high-redshift model predictions and intrinsic M$_{\rm BH}$-M$_{\star}$ relations of local BHs. We speculate that the observed non-linear BH-galaxy build-up is reflected in a twofold behavior with dark matter halo mass (M$_{\rm DM}$), displaying a clear turnover at M$_{\rm DM}\sim$2$\times$10$^{12}$M$_{\odot}$. While Supernovae-driven feedback suppresses BH growth in smaller halos (BHAR/SFR$\propto$M$_{\rm DM}^{1.6}$), above the M$_{\rm DM}$ threshold cold gas inflows possibly fuel both BH accretion and star formation in a similar fashion (BHAR/SFR$\propto$M$_{\rm DM}^{0.3}$).

preprint2019arXiv

The Intergalactic medium transmission towards z>4 galaxies with VANDELS and the impact of dust attenuation

Aims. Our aim is to estimate the intergalactic medium transmission towards UV-selected star-forming galaxies at redshift 4 and above and study the effect of the dust attenuation on these measurements. Methods. The ultra-violet spectrum of high redshift galaxies is a combination of their intrinsic emission and the effect of the Inter-Galactic medium (IGM) absorption along their line of sight. Using data coming from the unprecedented deep spectroscopy from the VANDELS ESO public survey carried out with the VIMOS instrument we compute both the dust extinction and the mean transmission of the IGM as well as its scatter from a set of 281 galaxies at z>3.87. Because of a degeneracy between the dust content of the galaxy and the IGM, we first estimate the stellar dust extinction parameter E(B-V) and study the result as a function of the dust prescription. Using these measurements as constraint for the spectral fit we estimate the IGM transmission Tr(Lyalpha). Both photometric and spectroscopic SED fitting are done using the SPectroscopy And photometRy fiTting tool for Astronomical aNalysis (SPARTAN) that is able to fit the spectral continuum of the galaxies as well as photometric data. Results. Using the classical Calzetti&#39;s attenuation law we find that E(B-V) goes from 0.11 at z=3.99 to 0.08 at z=5.15. These results are in very good agreement with previous measurements from the literature. We estimate the IGM transmission and find that the transmission is decreasing with increasing redshift from Tr(Lyalpha)=0.53 at z=3.99 to 0.28 at z=5.15. We also find a large standard deviation around the average transmission that is more than 0.1 at every redshift. Our results are in very good agreement with both previous measurements from AGN studies and with theoretical models.

preprint2019arXiv

The UV and Ly$α$ Luminosity Functions of galaxies and the Star Formation Rate Density at the end of HI reionization from the VIMOS Ultra-Deep Survey (VUDS)

We establish a robust statistical description of the star-forming galaxy population at the end of cosmic HI reionization ($5.0\le{}z\le6.6$) from a large sample of 52 galaxies with spectroscopically confirmed redshifts from the VIMOS UltraDeep Survey. We identify galaxies with Ly$α$ either in absorption or in emission, at variance with most spectroscopic samples in the literature where Ly$α$ emitters dominate. We find that star-forming galaxies at these redshifts are distributed along a main sequence in the stellar mass vs. SFR plane. We report a flat evolution of the sSFR(z) in 3<z<6 compared to lower redshift measurements. UV-continuum slopes vary with luminosity, with a large dispersion. We determing UV and Ly$α$ luminosity functions using V$_{max}$ method and use them to derive star formation rate densities (SFRD). We find that both UV-derived and Ly$α$-derived SFRDs are in excellent agreement after correcting Ly$α$ luminosity density for IGM absorption. Our new SFRD measurements at a mean redshift z=5.6 confirm the steep decline of the SFRD at z>2. The bright end of the Ly$α$ luminosity function has a high number density, indicating a significant star formation activity concentrated in the brightest Ly$α$ emitters (LAE) at these redshifts. LAE with EW>25Å~contribute to about 75\% of the total UV-derived SFRD. While our analysis favors a low dust content in 5.0<z<6.6, uncertainties on the dust extinction correction and associated degeneracies in spectral fitting will remain an issue to estimate the total SFRD until future survey extending spectroscopy to the NIR rest-frame spectral domain, e.g. with JWST.