Source author record

M. Castellano

M. Castellano appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

123works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

123 published item(s)

preprint2026arXiv

A Pristine Star-Forming Complex at z=4.19

We report the discovery of a faint (M_1700 ~ -12.2), oxygen-deficient strongly-lensed ionizing source -- dubbed LAP2 (Lensed And Pristine 2) -- at a spectroscopic redshift of z=4.19. LAP2 appears to be isolated and lies very close to the caustic produced by the lensing galaxy cluster Abell 2744. It was observed with the James Webb Space Telescope (JWST) NIRSpec MSA in prism mode as part of the UNCOVER program. The NIRSpec spectrum reveals prominent Lya emission (7.1 sigma), clear Ha emission (6.2 sigma), tentative Hb emission (2.8 sigma), and no detectable [OIII]4959,5007 (~ 7 times fainter than Ha). The inferred [OIII] 2 sigma upper limit corresponds to an R3 index <0.85 (assuming the Ha/Hb ~ 2.86 case~B recombination ratio), which, under high-ionization conditions, implies a metallicity of Z < 6 X 10^(-3) Z_sun. The combination of faint ultraviolet luminosity, large rest-frame Ha equivalent width (~ 650 A), and extremely compact size (< 10 pc) suggests that LAP2 is being caught in an early, pristine formation phase consistent with an instantaneous-burst scenario, with an estimated stellar mass of at most a few 10^4 Msun. Deep VLT/MUSE observations further reveal copious Lya emission forming an arclet that straddles the critical line. LAP2 thus joins the rare class of extremely metal-poor star-forming complexes that JWST has started to unveil at redshifts 3 - 7, and provides a rare glimpse into a still very poorly explored faint-luminosity regime.

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's Flagship simulation to investigate its impact on Euclid'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'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'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'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

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') 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's Data Release 1 (DR1).

preprint2025arXiv

Pushing JWST to the extremes: search and scrutiny of bright galaxy candidates at z$\simeq$15-30

We designed customized Lyman-break color selection techniques to identify galaxy candidates in the redshift ranges $15 \leq z \leq 20$ and $20 \leq z \leq 28$. The selection was performed on the ASTRODEEP-JWST multi-band catalogs of the CEERS, Abell-2744, JADES, NGDEEP, and PRIMER survey fields, covering a total area of $\sim0.2$ sq. deg. We identify five candidates at $15 \leq z \leq 20$, while no objects are found based on the $z\gtrsim20$ color selection criteria. Despite exhibiting a $>$1.5 mag break, all the objects display multimodal redshift probability distributions across different SED-fitting codes and methodologies. The alternative solutions correspond to poorly understood populations of low-mass quiescent or dusty galaxies at z$\sim$3-7. This conclusion is supported by the analysis of five F200W-dropout objects that we find to be interlopers on the basis of NIRSpec PRISM spectra: four dusty star-forming galaxies at z$\sim$2.2-6.6, and a passive galaxy at z=4.91 with log$(M_{\rm star}/{\rm M}_{\odot}) \lesssim$ 9. We measured the UV luminosity function under different assumptions on the contamination level within our sample. We find that if even a fraction of the candidates is indeed at $z\gtrsim15$, the resulting UV LF points to a very mild evolution compared to estimates at $z<15$, implying a significant tension with existing theoretical models. In particular, confirming our bright ($M_{\text{UV}}<-21$) candidates would require substantial revisions to the theoretical framework. In turn, if all these candidates will be confirmed to be interlopers, we conclude that future surveys may need ten times wider areas to select $M_{\text{UV}}\lesssim-20$ galaxies at $z>15$. Observations in the F150W and F200W filters at depths comparable to those in the NIRCam LW bands are also required to mitigate contamination from rare red objects at z$\lesssim$8.

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

Early results from GLASS-JWST. X: Rest-frame UV-optical properties of galaxies at 7 < z < 9

We present the first James Webb Space Telescope/NIRCam-led determination of $7<z<9$ galaxy properties based on broadband imaging from 0.8 to 5~$\mathrm{μm}$ as part of the GLASS-JWST Early Release Science program. This is the deepest dataset acquired at these wavelengths to date, with an angular resolution $\lesssim0.14$ arcsec. We robustly identify 13 galaxies with $S/N\gtrsim8$ in F444W from 8 arcmin$^2$ of data at $m_{AB}\leq 28$ from a combination of dropout and photometric redshift selection. From simulated data modeling, we estimate the dropout sample purity to be $\gtrsim90\%$. We find that the number density of these F444W-selected sources is broadly consistent with expectations from the UV luminosity function determined from Hubble Space Telescope data. We characterize galaxy physical properties using a Bayesian Spectral Energy Distribution fitting method, finding median stellar mass $10^{8.5}M_\odot$ and age 140 Myr, indicating they started ionizing their surroundings at redshift $z>9.5$. Their star formation main sequence is consistent with predictions from simulations. Lastly, we introduce an analytical framework to constrain main-sequence evolution at $z>7$ based on galaxy ages and basic assumptions, through which we find results consistent with expectations from cosmological simulations. While this work only gives a glimpse of the properties of typical galaxies that are thought to drive the reionization of the universe, it clearly shows the potential of JWST to unveil unprecedented details on galaxy formation in the first billion years.

preprint2022arXiv

Early results from GLASS-JWST. XI: Stellar masses and mass-to-light ratio of z>7 galaxies

We exploit James Webb Space Telescope (JWST) NIRCam observations from the GLASS-JWST-Early Release Science program to investigate galaxy stellar masses at z>7. We first show that JWST observations reduce the uncertainties on the stellar mass by a factor of at least 5-10, when compared with the highest quality data sets available to date. We then study the UV mass-to-light ratio, finding that galaxies exhibit a two orders of magnitude range of M/L_UV values for a given luminosity, indicative of a broad variety of physical conditions and star formation histories. As a consequence, previous estimates of the cosmic star stellar mass density - based on an average correlation between UV luminosity and stellar mass - can be biased by as much as a factor of ~6. Our first exploration demonstrates that JWST represents a new era in our understanding of stellar masses at z>7, and therefore of the growth of galaxies prior to cosmic reionization.

preprint2022arXiv

Early Results From GLASS-JWST. XII: The Morphology of Galaxies at the Epoch of Reionization

Star-forming galaxies can exhibit strong morphological differences between the rest-frame far-UV and optical, reflecting inhomogeneities in star-formation and dust attenuation. We exploit deep, high resolution NIRCAM 7-band observations to take a first look at the morphology of galaxies in the epoch of reionization ($z>7$), and its variation in the rest-frame wavelength range between Lyman $α$ and 6000-4000Å, at $z=7-12$. We find no dramatic variations in morphology with wavelength -- of the kind that would have overturned anything we have learned from the Hubble Space Telescope. No significant trends between morphology and wavelengths are detected using standard quantitative morphology statistics. We detect signatures of mergers/interactions in 4/19 galaxies. Our results are consistent with a scenario in which Lyman Break galaxies -- observed when the universe is only 400-800 Myrs old - are growing via a combination of rapid galaxy-scale star formation supplemented by accretion of star forming clumps and interactions.

preprint2022arXiv

Early results from GLASS-JWST. XIII. A faint, distant, and cold brown dwarf

We present the serendipitous discovery of a late T-type brown dwarf candidate in JWST NIRCam observations of the Early Release Science Abell 2744 parallel field. The discovery was enabled by the sensitivity of JWST at 4~$μ$m wavelengths and the panchromatic 0.9--4.5~$μ$m coverage of the spectral energy distribution. The unresolved point source has magnitudes F115W = 27.95$\pm$0.15 and F444W = 25.84$\pm$0.01 (AB), and its F115W$-$F444W and F356W$-$F444W colors match those expected for other, known T dwarfs. We can exclude it as a reddened background star, high redshift quasar, or very high redshift galaxy. Comparison with stellar atmospheric models indicates a temperature of $T_{eff}$ $\approx$ 600~K and surface gravity $\log{g}$ $\approx$ 5, implying a mass of 0.03~M$_{\odot}$ and age of 5~Gyr. We estimate the distance of this candidate to be 570--720~pc in a direction perpendicular to the Galactic plane, making it a likely thick disk or halo brown dwarf. These observations underscore the power of JWST to probe the very low-mass end of the substellar mass function in the Galactic thick disk and halo.

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

Exploring the physical properties of lensed star-forming clumps at $2\lesssim z \lesssim6$

We study the physical properties (size, stellar mass, luminosity, star formation rate) and scaling relations for a sample of 166 star-forming clumps with redshift $z \sim 2-6.2$. They are magnified by the Hubble Frontier Field galaxy cluster MACS~J0416 and have robust lensing magnification ($2\lesssim μ\lesssim 82$) computed by using our high-precision lens model, based on 182 multiple images. Our sample extends by $\sim 3$ times the number of spectroscopically-confirmed lensed clumps at $z \gtrsim 2$. We identify clumps in ultraviolet continuum images and find that, whenever the effective spatial resolution (enhanced by gravitational lensing) increases, they fragment into smaller entities, likely reflecting the hierarchically-organized nature of star formation. Kpc-scale clumps, most commonly observed in field, are not found in our sample. The physical properties of our sample extend the parameter space typically probed by $z \gtrsim 1$ field observations and simulations, by populating the low mass (M$_\star \lesssim 10^7$ M$_\odot$), low star formation rate (SFR $\lesssim 0.5$ M$_\odot$ yr$^{-1}$), and small size (R$_\mathrm{eff} \lesssim 100$ pc) regime. The new domain probed by our study approaches the regime of compact stellar complexes and star clusters. In the mass-size plane, our sample spans the region between galaxies and globular clusters, with a few clumps in the region populated by young star clusters and globular-clusters. For the bulk of our sample, we measure star-formation rates which are higher than those observed locally in compact stellar systems, indicating different conditions for star formation at high redshift than in the local Universe.

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' 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 environmental dependence of the stellar and gas-phase mass-metallicity relation at 2 < z < 4

In the local universe, galaxies in clusters show different properties compared to more isolated systems. Understanding how this difference originates and whether it is already in place at high redshift is still a matter of debate. Thanks to uniquely deep optical spectra from the VANDELS survey, we investigate environmental effects on the stellar mass-metallicity relation (MZR) for a sample of ~1000 star-forming galaxies at 2<z<4. We complement our dataset with MOSFIRE follow-up of 21 galaxies to study the environmental dependence of the gas-phase MZR. Robust stellar and gas metallicities are derived, respectively, from well-calibrated photospheric absorptions features at 1501 and 1719 Åin the stacked spectra, and from optical emission lines ([OII]3726-3729, [OIII]5007, and Hbeta) in individual systems. We characterize the environment through multiple criteria by using the local galaxy density maps previously derived in VANDELS. We find that environmental effects are weak at these redshifts, and more important around the densest overdensity structures, where galaxies have a lower stellar Z (by 0.2 dex) and a lower gas-phase Z (by 0.1 dex) compared to the field, with a significance of 1 and 2 sigma, respectively. Crucially, this offset cannot be explained by a selection effect due to a higher SFR, a fainter UV continuum, or different dust attenuations and stellar ages. Despite the still low S/N of our results, we propose a combination of increased mergers and high-speed encounters, more efficient AGN feedback in dense cores, and cold gas inflows as viable mechanisms diluting the metal content of overdense galaxies or expelling their metals to the IGM. Finally, some tensions remain with semi-analytic models and hydrodynamical simulations, which predict no significant offset as a function of host halo mass, suggesting that an explicit implementation of environmental processes is needed.

preprint2022arXiv

The GLASS James Webb Space Telescope Early Release Science Program. I. Survey Design and Release Plans

The GLASS James Webb Space Telescope Early Release Science (hereafter GLASS-JWST-ERS) Program will obtain and make publicly available the deepest extragalactic data of the ERS campaign. It is primarily designed to address two key science questions, namely, "what sources ionized the universe and when?" and "how do baryons cycle through galaxies?", while also enabling a broad variety of first look scientific investigations. In primary mode, it will obtain NIRISS and NIRSpec spectroscopy of galaxies lensed by the foreground Hubble Frontier Field cluster, Abell 2744. In parallel, it will use NIRCam to observe two fields that are offset from the cluster center, where lensing magnification is negligible, and which can thus be effectively considered blank fields. In order to prepare the community for access to this unprecedented data, we describe the scientific rationale, the survey design (including target selection and observational setups), and present pre-commissioning estimates of the expected sensitivity. In addition, we describe the planned public releases of high-level data products, for use by the wider astronomical community.

preprint2022arXiv

The ionizing properties of two bright Ly$α$ emitters in the BDF reionized bubble at z=7

We investigate the ionizing properties of the pair of bright Ly$α$ emitting galaxies BDF521 and BDF2195 at z=7.012 in order to constrain their contribution to the formation of the BDF "reionized bubble". We obtain constraints on UV emission lines (CIV$λ1548$ doublet, HeII$λ1640$, OIII]$λ1660$ doublet, and CIII]$λ1909$ doublet) from deep VLT-XSHOOTER observations and compare them to those available for other high-redshift objects, and to models with mixed stellar and AGN emission. We use this spectroscopic information together with the photometry available in the field to constrain the physical properties of the two objects using the spectro-photometric fitting code BEAGLE. We do not detect any significant emission at the expected position of UV lines, with 3$σ$ upper limits of EW$\lesssim$2-7 AA rest-frame. We find that the two objects have lower CIII] emission than expected on the basis of the correlation between the Ly$α$ and CIII] EWs. The EW limits on CIV and HeII emission exclude pure AGN templates at $\sim2-3σ$ significance, and only models with a $\lesssim$40% AGN contribution are compatible with the observations. The two objects are found to be relatively young ($\sim$20-30 Myrs) and metal-poor ($\lesssim 0.3 Z_{\odot}$) with stellar masses of a few $10^9M_{\odot}$. Their production rate of hydrogen ionizing photons per intrinsic UV luminosity is log($ξ_{ion}^*$/Hz erg$^{-1}$)=25.02-25.26, consistent with values typically found in high-redshift galaxies, but more than twice lower than values measured in $z>$7 galaxies with strong CIII] and/or optical line emission ($\simeq$25.6-25.7). The two BDF emitters have no evidence of higher than average ionizing capabilities and are not capable of reionizing their surroundings by their own means under realistic assumptions on the escape fraction of ionizing photons. (abridged)

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.

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

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'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: Forecast constraints on the cosmic distance duality relation with complementary external probes

In metric theories of gravity with photon number conservation, the luminosity and angular diameter distances are related via the Etherington relation, also known as the distance-duality relation (DDR). A violation of this relation would rule out the standard cosmological paradigm and point at the presence of new physics. We quantify the ability of Euclid, in combination with contemporary surveys, to improve the current constraints on deviations from the DDR in the redshift range $0<z<1.6$. We start by an analysis of the latest available data, improving previously reported constraints by a factor of 2.5. We then present a detailed analysis of simulated Euclid and external data products, using both standard parametric methods (relying on phenomenological descriptions of possible DDR violations) and a machine learning reconstruction using Genetic Algorithms. We find that for parametric methods Euclid can (in combination with external probes) improve current constraints by approximately a factor of six, while for non-parametric methods Euclid can improve current constraints by a factor of three. Our results highlight the importance of surveys like Euclid in accurately testing the pillars of the current cosmological paradigm and constraining physics beyond the standard cosmological model.

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

High star cluster formation efficiency in the strongly lensed Sunburst Lyman-continuum galaxy at z=2.37

We investigate the strongly lensed (μx10-100) Lyman continuum (LyC) galaxy, dubbed Sunburst, at z=2.37, taking advantage of a new accurate model of the lens. A characterization of the intrinsic (delensed) properties of the galaxy yields a size of ~3 sq.kpc, a luminosity Muv=-20.3,and a stellar mass M~10^9 Msun;16% of the ultraviolet light is located in a 3 Myr old gravitationally-bound young massive star cluster (YMC) with an effective radius of Re~8 pc and a dynamical mass of ~10^7 Msun (similar to the stellar mass), from which LyC radiation is detected (λ< 912A). The inferred outflowing gas velocity (>300 km/s) exceeds the escape velocity of the star cluster. The resulting escape fraction of the ionizing radiation emerging from the Sunburst galaxy is >6-12%, whilst it is >46-93% if inferred from the YMC. 12 additional likely star clusters with 3<Re<20 pc are identified in the galaxy from which we derive a cluster formation efficiency Γ>~30%, which is consistent with the high Γderived in local galaxies experiencing extreme gas physical conditions. The presence of the YMC influences the morphology (nucleation), photometry (photometric jumps) and spectroscopic output (nebular emission) of the entire galaxy. The de-lensed LyC and UV (1600A) magnitudes of the YMC are ~30.6 and ~26.9, whilst the galaxy has m1600~24.8. A relatively large rest-frame equivalent width of EWrest(Hb+[OIII]4959-5007)~450A emerges from the galaxy with the YMC contributing to ~30%. If O-type stars are mainly forged in star clusters, then such engines were the key ionizing agents during reionization and the increasing occurrence of high EW lines (Hb+[OIII]) observed at z>6.5 might be an indirect signature of a high Γat reionization.Future facilities (like VLT/MAVIS or ELT), will probe bound clusters on moderately magnified (μ<5-10) galaxies across cosmic epochs up to reionization[ABRIDGED]

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

Candidate Population III stellar complex at z=6.629 in the MUSE Deep Lensed Field

We discovered a strongly lensed (μ>40) Lya emission at z=6.629 (S/N~18) in the MUSE Deep Lensed Field (MDLF) targeting the Hubble Frontier Field galaxy cluster MACS~J0416. Dedicated lensing simulations imply that the Lya emitting region necessarily crosses the caustic. The arc-like shape of the Lya extends 3 arcsec on the observed plane and is the result of two merged multiple images, each one with a de-lensed Lya luminosity L<~2.8 x 10^(40) erg/s arising from a confined region (< 150 pc effective radius). A spatially unresolved HST counterpart is barely detected at S/N~2 after stacking the near-infrared bands, corresponding to an observed(intrinsic) magnitude m_(1500)>~30.8(>~35.0). The inferred rest-frame Lya equivalent width is EWo > 1120 A if the IGM transmission is T(IGM)<0.5. The low luminosities and the extremely large Lya EWo match the case of a Population~III star complex made of several dozens stars (~ 10^4 Msun) which irradiate a HII region crossing the caustic. While the Lya and stellar continuum are among the faintest ever observed at this redshift, the continuum and the Lya emissions could be affected by differential magnification, possibly biasing the EWo estimate. The aforementioned tentative HST detection tend to favor a large EWo, making such a faint Pop~III candidate a key target for the James Webb Space Telescope and Extremely Large Telescopes.

preprint2020arXiv

Constraints on Dynamical Dark Energy Models from the Abundance of Massive Galaxies at High Redshifts

We compare the maximal abundance of massive systems predicted in different dynamical dark energy (DDE) models at high redshifts z = 4-7 with the measured abundance of the most massive galaxies observed to be already in place at such redshifts. The aim is to derive constraints for the evolution of the dark energy equation of state parameter w which are complementary to existing probes. We adopt the standard parametrization for the DDE evolution in terms of the local value w_0 and of the look-back time derivative w_a of the equation of state. We derive constraints on combinations (w_0, w_a) in the different DDE models by using three different, independent probes: (i) the observed stellar mass function of massive objects at z = 6 derived from the CANDELS survey; (ii) the estimated volume density of massive halos derived from the observation of massive, star-forming galaxies detected in the submillimeter range at z = 4; (iii) The rareness of he most massive system (estimated gas mass exceeding 3 10^11 M_sun) observed to be in place at z = 7, a far-infrared-luminous object recently detected in the South Pole Telescope (SPT) survey. Finally, we show that the combination of our results from the three above probes excludes a sizable fraction of the DDE parameter space w_a > -3/4 - (w_0 + 3/2) presently allowed (or even favored) by existing probes.

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."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: Reconstruction of Weak Lensing mass maps for non-Gaussianity studies

Weak lensing, which is the deflection of light by matter along the line of sight, has proven to be an efficient method for constraining models of structure formation and reveal the nature of dark energy. So far, most weak-lensing studies have focused on the shear field that can be measured directly from the ellipticity of background galaxies. However, within the context of forthcoming full-sky weak-lensing surveys such as Euclid, convergence maps (mass maps) offer an important advantage over shear fields in terms of cosmological exploitation. While it carry the same information, the lensing signal is more compressed in the convergence maps than in the shear field. This simplifies otherwise computationally expensive analyses, for instance, non-Gaussianity studies. However, the inversion of the non-local shear field requires accurate control of systematic effects caused by holes in the data field, field borders, shape noise, and the fact that the shear is not a direct observable (reduced shear). We present the two mass-inversion methods that are included in the official Euclid data-processing pipeline: the standard Kaiser & Squires method (KS), and a new mass-inversion method (KS+) that aims to reduce the information loss during the mass inversion. This new method is based on the KS method and includes corrections for mass-mapping systematic effects. The results of the KS+ method are compared to the original implementation of the KS method in its simplest form, using the Euclid Flagship mock galaxy catalogue. In particular, we estimate the quality of the reconstruction by comparing the two-point correlation functions and third- and fourth- order moments obtained from shear and convergence maps, and we analyse each systematic effect independently and simultaneously...

preprint2020arXiv

Euclid: The reduced shear approximation and magnification bias for Stage IV cosmic shear experiments

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

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

The MUSE Deep Lensed Field on the Hubble Frontier Field MACS~J0416

Context: A census of faint and tiny star forming complexes at high redshift is key to improving our understanding of reionizing sources, galaxy growth and the formation of globular clusters. Aims: We present the MUSE Deep Lensed Field (MDLF) program. Methods: We describe Deep MUSE observations of 17.1 hours integration on a single pointing over the Hubble Frontier Field galaxy cluster MACS~J0416. Results: We confirm spectroscopic redshifts for all 136 multiple images of 48 source galaxies at 0.9<z<6.2. Within those galaxies, we securely identify 182 multiple images of 66 galaxy components that we use to constrain our lens model. We identify 116 clumps belonging to background high-z galaxies; the majority of them are multiple images and span magnitude, size and redshift intervals of [-18,-10], [~400-3] parsec and 1<z<6.6, respectively, with the most magnified ones probing possible single gravitationally bound star clusters. The depth of the MDLF combined with lensing magnification lead us to reach a detection limit for unresolved emission lines of a few 10$^{-20}$ erg/s/cm2, after correction for lensing magnification. Ultraviolet high-ionization metal lines (and HeII1640) are detected with S/N>10 for individual objects down to de-lensed magnitude 28-30 suggesting that they are common in such faint sources. Conclusions:Deep MUSE observations, in combination with existing HST imaging, allowed us to:(1) confirm redshifts for extremely faint high-z sources;(2) peer into their internal clumps (down to 100-200 pc scale);(3) in some cases break down such clumps into bound star clusters (<20 pc scale);(4) double the number of constraints for the lens model,reaching an unprecedented set of 182 bona-fide multiple images and confirming up to 213 galaxy cluster members. These results demonstrate the power that JWST and future ELTs will have when combined to study gravitational telescopes.[abridged]

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

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

Euclid preparation: VI. Verifying the Performance of Cosmic Shear Experiments

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

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

preprint2016arXiv

A Stringent Limit on the Warm Dark Matter Particle Masses from the Abundance of z=6 Galaxies in the Hubble Frontier Fields

We show that the recently measured UV luminosity functions of ultra-faint lensed galaxies at z= 6 in the Hubble Frontier Fields provide an unprecedented probe for the mass m_X of the Warm Dark Matter candidates independent of baryonic physics. Comparing the measured abundance of the faintest galaxies with the maximum number density of dark matter halos in WDM cosmologies sets a robust limit m_X> 2.9 keV for the mass of thermal relic WDM particles at a 1-sigma confidence level, m_X> 2.4 keV at 2-sigma, and m_X> 2.1 keV at 3-sigma. These constitute the tightest constraints on WDM particle mass derived to date independently of the baryonic physics involved in galaxy formation. We discuss the impact of our results on the production mechanism of sterile neutrinos. In particular, if sterile neutrinos are responsible for the 3.5 keV line reported in observations of X-ray clusters, our results firmly rule out the Dodelson-Widrow production mechanism, and yield m_{sterile}> 6.1 keV for sterile neutrinos produced via the Shi-Fuller mechanism.

preprint2016arXiv

Chandra counterparts of CANDELS GOODS-S sources

Improving the capabilities of detecting faint X-ray sources is fundamental to increase the statistics on faint high-z AGN and star-forming galaxies. We performed a simultaneous Maximum Likelihood PSF fit in the [0.5-2] keV and [2-7] keV energy bands of the 4 Ms{\em Chandra} Deep Field South (CDFS) data at the position of the 34930 CANDELS H-band selected galaxies. For each detected source we provide X-ray photometry and optical counterpart validation. We validated this technique by means of a raytracing simulation. We detected a total of 698 X-ray point-sources with a likelihood $\mathcal{L}$$>$4.98 (i.e. $>$2.7$σ$). We show that the prior knowledge of a deep sample of Optical-NIR galaxies leads to a significant increase of the detection of faint (i.e. $\sim$10$^{-17}$ cgs in the [0.5-2] keV band) sources with respect to "blind" X-ray detections. By including previous X-ray catalogs, this work increases the total number of X-ray sources detected in the 4 Ms CDFS, CANDELS area to 793, which represents the largest sample of extremely faint X-ray sources assembled to date. Our results suggest that a large fraction of the optical counterparts of our X-ray sources determined by likelihood ratio actually coincides with the priors used for the source detection. Most of the new detected sources are likely star-forming galaxies or faint absorbed AGN. We identified a few sources sources with putative photometric redshift z$>$4. Despite the low number statistics and the uncertainties on the photo-z, this sample significantly increases the number of X--ray selected candidate high-z AGN.

preprint2016arXiv

Characterizing elusive, faint dusty star-forming galaxies: a lensed, optically undetected ALMA galaxy at z~3.3

We present the serendipitous ALMA detection of a faint submillimeter galaxy (SMG) lensed by a foreground z~1 galaxy. By optimizing the source detection to deblend the system, we accurately build the full spectral energy distribution of the distant galaxy from the I814 band to radio wavelengths. It is extremely red, with a I-K colour larger than 2.5. We estimate a photometric redshift of 3.28 and determine the physical parameters. The distant galaxy turns out to be magnified by the foreground lens by a factor of ~1.5, which implies an intrinsic Ks-band magnitude of ~24.5, a submillimeter flux at 870um of ~2.5 mJy and a SFR of ~150-300Msun/yr, depending on the adopted tracer. These values place our source towards the faint end of the distribution of observed SMGs, and in particular among the still few faint SMGs with a fully characterized spectral energy distribution, which allows us not only to accurately estimate its redshift but also to measure its stellar mass and other physical properties. The galaxy studied in this work is a representative of the population of faint SMGs, of which only few objects are known to date, that are undetected in optical and therefore are not typically accounted for when measuring the cosmic star formation history (SFH). This faint galaxy population thus likely represents an important and missing piece in our understanding of the cosmic SFH. Its observation and characterization is of major importance to achieve a solid picture of galaxy evolution.

preprint2016arXiv

Constraining the Warm Dark Matter Particle Mass through Ultra-Deep UV Luminosity Functions at z=2

We compute the mass function of galactic dark matter halos for different values of the Warm Dark Matter (WDM) particle mass m_X and compare it with the abundance of ultra-faint galaxies derived from the deepest UV luminosity function available so far at redshift z~2. The magnitude limit M_UV=-13 reached by such observations allows us to probe the WDM mass functions down to scales close to or smaller than the half-mass mode mass scale ~10^9 M_sun. This allowed for an efficient discrimination among predictions for different m_X which turn out to be independent of the star formation efficiency adopted to associate the observed UV luminosities of galaxies to the corresponding dark matter masses. Adopting a conservative approach to take into account the existing theoretical uncertainties in the galaxy halo mass function, we derive a robust limit m_X>1.8 keV for the mass of thermal relic WDM particles when comparing with the measured abundance of the faintest galaxies, while m_X>1.5 keV is obtained when we compare with the Schechter fit to the observed luminosity function. The corresponding lower limit for sterile neutrinos depends on the modeling of the production mechanism; for instance m_sterile > 4 keV holds for the Shi-Fuller mechanism. We discuss the impact of observational uncertainties on the above bound on m_X. As a baseline for comparison with forthcoming observations from the HST Frontier Field, we provide predictions for the abundance of faint galaxies with M_UV=-13 for different values of m_X and of the star formation efficiency, valid up to z~4.

preprint2016arXiv

Constraints on photoionization feedback from number counts of ultra-faint high-redshift galaxies in the Frontier Fields

We exploit a sample of ultra-faint high-redshift galaxies (demagnified HST $H_{160}$ magnitude $>30$) in the Frontier Fields clusters A2744 and M0416 to constrain a theoretical model for the UV luminosity function (LF) in the presence of photoionization feedback. The objects have been selected on the basis of accurate photometric redshifts computed from multi-band photometry including 7 HST bands and deep $K_s$ and IRAC observations. Magnification is computed on an object-by-object basis from all available lensing models of the two clusters. We take into account source detection completeness as a function of luminosity and size, magnification effects and systematics in the lens modeling of the clusters under investigation. We find that our sample of high-$z$ galaxies constrain the cut-off halo circular velocity below which star-formation is suppressed by photo-ionization feedback to $v_c^{\rm cut} < 50$ km s$^{-1}$. This circular velocity corresponds to a halo mass of $\approx5.6\times10^9~M_\odot$ and $\approx2.3\times10^9~M_\odot$ at $z=5$ and 10 respectively: higher mass halos can thus sustain continuous star formation activity without being quenched by external ionizing flux. More stringent constraints are prevented by the uncertainty in the modeling of the cluster lens, as embodied by systematic differences among the lens models available.

preprint2016arXiv

First observational support for overlapping reionized bubbles generated by a galaxy overdensity

We present the analysis of deep HST multi-band imaging of the BDF field specifically designed to identify faint companions around two of the few Ly-alpha emitting galaxies spectroscopically confirmed at z~7 (Vanzella et al. 2011). Although separated by only 4.4 proper Mpc these galaxies cannot generate HII regions large enough to explain visibility of their Ly-alpha line, thus requiring a population of fainter ionizing sources in their vicinity. We use deep HST and VLT-Hawk-I data to select z~7 Lyman break galaxies around the emitters. We select 6 new robust z~7 LBGs at Y~26.5-27.5 whose average spectral energy distribution is consistent with the objects being at the redshift of the close-by Ly-alpha emitters. The resulting number density of z~7 LBGs in the BDF field is a factor ~3-4 higher than expected in random pointings of the same size. We compare these findings with cosmological hydrodynamic plus radiative transfer simulations of a universe with a half neutral IGM: we find that indeed Ly-alpha emitter pairs are only found in completely ionized regions characterized by significant LBG overdensities. Our findings match the theoretical prediction that the first ionization fronts are generated within significant galaxy overdensities and support a scenario where faint, "normal" star-forming galaxies are responsible for reionization.

preprint2016arXiv

High-resolution spectroscopy of a young, low-metallicity optically-thin L=0.02L* star-forming galaxy at z=3.12

We present VLT/X-Shooter and MUSE spectroscopy of an faint F814W=28.60+/-0.33 (Muv=-17.0), low mass (~<10^7 Msun) and compact (Reff=62pc) freshly star-forming galaxy at z=3.1169 magnified (16x) by the Hubble Frontier Fields galaxy cluster Abell S1063. Gravitational lensing allows for a significant jump toward low-luminosity regimes, in moderately high resolution spectroscopy (R=lambda/dlambda ~ 3000-7400). We measured CIV1548,1550, HeII1640, OIII]1661,1666, CIII]1907,1909, Hbeta, [OIII]4959,5007, emission lines with FWHM< 50 km/s and (de-lensed) fluxes spanning the interval 1.0x10^-19 - 2.0x10^-18 erg/s/cm2 at S/N=4-30. The double peaked Lya emission with Delta_v(red-blue) = 280(+/-7)km/s and de-lensed fluxes 2.4_(blue)|8.5_(red)x10^-18 erg/s/cm2 (S/N=38_(blue)|110_(red)) indicate a low column density of neutral hydrogen gas consistent with a highly ionized interstellar medium as also inferred from the large [OIII]5007/[OII]3727>10 ratio. We detect CIV1548,1550 resonant doublet in emission, each component with FWHM ~< 45 km/s, and redshifted by +51(+/-10)km/s relative to the systemic redshift. We interpret this as nebular emission tracing an expanding optically-thin interstellar medium. Both CIV1548,1550 and HeII1640 suggest the presence of hot and massive stars (with a possible faint AGN). The ultraviolet slope is remarkably blue, beta =-2.95 +/- 0.20 (F_lambda=lambda^beta), consistent with a dust-free and young ~<20 Myr galaxy. Line ratios suggest an oxygen abundance 12+log(O/H)<7.8. We are witnessing an early episode of star-formation in which a relatively low NHI and negligible dust attenuation might favor a leakage of ionizing radiation. This galaxy currently represents a unique low-luminosity reference object for future studies of the reionization epoch with JWST.

preprint2016arXiv

Hubble imaging of the ionizing radiation from a star-forming galaxy at z=3.2 with fesc>50%

Star-forming galaxies are considered to be the leading candidate sources that dominate the cosmic reionization at z>7, and the search for analogs at moderate redshift showing Lyman continuum (LyC) leakage is currently a active line of research. We have observed a star-forming galaxy at z=3.2 with Hubble/WFC3 in the F336W filter, corresponding to the 730-890A rest-frame, and detect LyC emission. This galaxy is very compact and also has large Oxygen ratio [OIII]5007/[OII]3727 (>=10). No nuclear activity is revealed from optical/near-infrared spectroscopy and deep multi-band photometry (including the 6Ms X-ray, Chandra). The measured escape fraction of ionizing radiation spans the range 50-100\%, depending on the IGM attenuation. The LyC emission is detected at S/N=10 with m(F336W)=27.57+/-0.11 and it is spatially unresolved, with effective radius R_e<200pc. Predictions from photoionization and radiative transfer models are in line with the properties reported here, indicating that stellar winds and supernova explosions in a nucleated star-forming region can blow cavities generating density-bounded conditions compatible with optically thin media. Irrespective to the nature of the ionizing radiation, spectral signatures of these sources over the entire electromagnetic spectrum are of central importance for their identification during the epoch of reionization, when the LyC is unobservable. Intriguingly, the Spitzer/IRAC photometric signature of intense rest-frame optical emissions ([OIII]+Hbeta) observed recently at z~7.5-8.5 is similar to what is observed in this galaxy. Only the James Webb Space Telescope will measure optical line ratios at z>7 allowing a direct comparison with lower redshift LyC emitters, as reported here.

preprint2016arXiv

Limits on the LyC signal from z~3 sources with secure redshift and HST coverage in the E-CDFS field

Aim: We aim to measure the LyC signal from a sample of sources in the Chandra deep field south. We collect star-forming galaxies (SFGs) and active galactic nuclei (AGN) with accurate spectroscopic redshifts, for which Hubble Space Telescope (HST) coverage and multi-wavelength photometry are available. Method: We selected a sample of about 200 sources at z~3. Taking advantage of HST resolution, we applied a careful cleaning procedure and rejected sources showing nearby clumps with different colours, which could be lower-z interlopers. Our clean sample consisted of 86 SFGs (including 19 narrow-band selected Lya emitters) and 8 AGN (including 6 detected in X-rays). We measured the LyC flux from aperture photometry in four narrow-band filters covering wavelengths below a 912 A rest frame (3.11<z<3.53). We estimated the ratio between ionizing (LyC flux) and 1400 A non-ionizing emissions for AGN and galaxies. Results: By running population synthesis models, we assume an average intrinsic L(1400 A)/L(900 A) ratio of 5 as the representative value for our sample. With this value and an average treatment of the lines of sight of the inter-galactic medium, we estimate the LyC escape fraction relative to the intrinsic value (fesc_rel(LyC)). We do not directly detect ionizing radiation from any individual SFG, but we are able to set a 1(2)sigma upper limit of fesc_rel(LyC)<12(24)%. This result is consistent with other non-detections published in the literature. No meaningful limits can be calculated for the sub-sample of Lya emitters. We obtain one significant direct detection for an AGN at z=3.46, with fesc_rel(LyC) = (72+/-18)%. Conclusions: Our upper limit on fescrel(LyC) implies that the SFGs studied here do not present either the physical properties or the geometric conditions suitable for efficient LyC-photon escape.

preprint2016arXiv

T-PHOT version 2.0: improved algorithms for background subtraction, local convolution, kernel registration, and new options

We present the new release v2.0 of T-PHOT, a publicly available software package developed to perform PSF-matched, prior-based, multiwavelength deconfusion photometry of extragalactic fields. New features included in the code are presented and discussed: background estimation, fitting using position dependent kernels, flux prioring, diagnostical statistics on the residual image, exclusion of selected sources from the model and residual images, individual registration of fitted objects. These new options improve on the performance of the code, allowing for more accurate results and providing useful aids for diagnostics.

preprint2016arXiv

The ASTRODEEP Frontier Fields Catalogues: I - Multiwavelength photometry of Abell-2744 and MACS-J0416

We present multiwavelength photometric catalogues (HST, Spitzer and Hawk-I K band) for the first two of the Frontier Fields, Abell2744 and MACSJ0416 (plus their parallel fields). To detect faint sources even in the central regions of the clusters, we develop a robust and repeatable procedure that uses the public codes Galapagos and Galfit to model and remove most of the light contribution from both the brightest cluster members as well as the ICL. We perform the detection on the HST H160 processed image to obtain a pure H-selected sample. We also add a sample of sources which are undetected in the H160 image but appear on a stacked infrared image. Photometry in the other HST bands is obtained using SExtractor, performed again on residual images after the Galfit procedure for foreground light removal. Photometry on the Hawk-I and IRAC bands has been obtained using our PSF-matching deconfusion code T-PHOT. A similar procedure, but without the need for the foreground light removal, is adopted for the Parallel fields. The procedure allows for the detection and the photometric measurements of ~2500 sources per field. We deliver and release complete photometric H-detected catalogues, with the addition of a complementary sample of infrared-detected sources. All objects have multiwavelength coverage including B to H HST bands, plus K band from Hawk-I, and 3.6 - 4.5 μm from Spitzer. Full and detailed treatment of photometric errors is included. We perform basic sanity checks on the reliability of our results. The multiwavelength catalogues are publicly available and are ready to be used for scientific purposes. Our procedures allows for the detection of outshined objects near the bright galaxies, which, coupled with the magnification effect of the clusters, can reveal extremely faint high redshift sources. Full analysis on photometric redshifts is presented in a companion Paper II. [abridged]

preprint2016arXiv

The ASTRODEEP Frontier Fields Catalogues: II - Photometric redshifts and rest-frame properties in Abell-2744 and MACS-J0416

We present the first public release of photometric redshifts, galaxy rest-frame properties and associated magnification values in the cluster and parallel pointings of the first two Frontier Fields, Abell-2744 and MACS-J0416. We exploit a multi-wavelength catalogue ranging from HST to ground-based K and Spitzer IRAC which is specifically designed to enable detection and measurement of accurate fluxes in crowded cluster regions. The multi-band information is used to derive photometric redshifts and physical properties of sources detected either in the H-band image alone or from a stack of four WFC3 bands. To minimize systematics median photometric redshifts are assembled from six different approaches to photo-z estimates. Their reliability is assessed through a comparison with available spectroscopic samples. State of the art lensing models are used to derive magnification values on an object-by-object basis by taking into account sources positions and redshifts. We show that photometric redshifts reach a remarkable ~3-5% accuracy. After accounting for magnification the H band number counts are found in agreement at bright magnitudes with number counts from the CANDELS fields, while extending the presently available samples to galaxies intrinsically as faint as H160~32-33 thanks to strong gravitational lensing. The Frontier Fields allow to probe the galaxy stellar mass distribution at 0.5-1.5 dex lower masses, depending on magnification, with respect to extragalactic wide fields, including sources at Mstar~ 10^7-10^8 Msun at z>5. Similarly, they allow the detection of objects with intrinsic SFRs>1dex lower than in the CANDELS fields reaching 0.1-1 Msun/yr at z~6-10. [abridged]

preprint2016arXiv

The VIMOS Ultra Deep Survey: Ly$α$ Emission and Stellar Populations of Star-Forming Galaxies at 2<z<2.5

The aim of this paper is to investigate spectral and photometric properties of 854 faint ($i_{AB}$<~25 mag) star-forming galaxies (SFGs) at 2<z<2.5 using the VIMOS Ultra-Deep Survey (VUDS) spectroscopic data and deep multi-wavelength photometric data in three extensively studied extragalactic fields (ECDFS, VVDS, COSMOS). These SFGs were targeted for spectroscopy based on their photometric redshifts. The VUDS spectra are used to measure the UV spectral slopes ($β$) as well as Ly$α$ equivalent widths (EW). On average, the spectroscopically measured $β$ (-1.36$\pm$0.02), is comparable to the photometrically measured $β$ (-1.32$\pm$0.02), and has smaller measurement uncertainties. The positive correlation of $β$ with the Spectral Energy Distribution (SED)-based measurement of dust extinction, E$_{\rm s}$(B-V), emphasizes the importance of $β$ as an alternative dust indicator at high redshifts. To make a proper comparison, we divide these SFGs into three subgroups based on their rest-frame Ly$α$ EW: SFGs with no Ly$α$ emission (SFG$_{\rm N}$; EW$\le$0Å), SFGs with Ly$α$ emission (SFG$_{\rm L}$; EW$>$0Å), and Ly$α$ emitters (LAEs; EW$\ge$20Å). The fraction of LAEs at these redshifts is $\sim$10%, which is consistent with previous observations. We compared best-fit SED-estimated stellar parameters of the SFG$_{\rm N}$, SFG$_{\rm L}$ and LAE samples. For the luminosities probed here ($\sim$L$^*$), we find that galaxies with and without Ly$α$ in emission have small but significant differences in their SED-based properties. We find that LAEs have less dust, and lower star-formation rates (SFR) compared to non-LAEs. We also find that LAEs are less massive compared to non-LAEs, though the difference is smaller and less significant compared to the SFR and E$_{\rm s}$(B-V). [abridged]

preprint2016arXiv

Tracing the reionization epoch with ALMA: [CII] emission in z~7 galaxies

We present new results on [CII]158$μ$ m emission from four galaxies in the reionization epoch. These galaxies were previously confirmed to be at redshifts between 6.6 and 7.15 from the presence of the Ly$α$ emission line in their spectra. The Ly$α$ emission line is redshifted by 100-200 km/s compared to the systemic redshift given by the [CII] line. These velocity offsets are smaller than what is observed in z~3 Lyman break galaxies with similar UV luminosities and emission line properties. Smaller velocity shifts reduce the visibility of Ly$α$ and hence somewhat alleviate the need for a very neutral IGM at z~7 to explain the drop in the fraction of Ly$α$ emitters observed at this epoch. The galaxies show [CII] emission with L[CII]=0.6-1.6 x10$^8 L_\odot$: these luminosities place them consistently below the SFR-L[CII] relation observed for low redshift star forming and metal poor galaxies and also below z =5.5 Lyman break galaxies with similar star formation rates. We argue that previous undetections of [CII] in z~7 galaxies with similar or smaller star formation rates are due to selection effects: previous targets were mostly strong Ly$α$ emitters and therefore probably metal poor systems, while our galaxies are more representative of the general high redshift star forming population .

preprint2015arXiv

ALMA constraints on the faint millimetre source number counts and their contribution to the cosmic infrared background

We have analysed 18 ALMA continuum maps in Bands 6 and 7, with rms down to 7.8$μ$Jy, to derive differential number counts down to 60$μ$Jy and 100$μ$Jy at $λ=$1.3 mm and $λ=$1.1 mm, respectively. The area covered by the combined fields is $\rm 9.5\times10^{-4}deg^2$ at 1.1mm and $\rm 6.6\times10^{-4}deg^{2}$ at 1.3mm. We improved the source extraction method by requiring that the dimension of the detected sources be consistent with the beam size. This method enabled us to remove spurious detections that have plagued the purity of the catalogues in previous studies. We detected 50 faint sources with S/N$>$3.5 down to 60$μ$Jy, hence improving the statistics by a factor of four relative to previous studies. The inferred differential number counts are $\rm dN/d(Log_{10}S)=1\times10^5~deg^2$ at a 1.1 mm flux $S_{λ= 1.1~mm} = 130~μ$Jy, and $\rm dN/d(Log_{10}S)=1.1\times10^5~deg^2$ at a 1.3 mm flux $\rm S_{λ= 1.3~mm} = 60~μ$Jy. At the faintest flux limits, i.e. 30$μ$Jy and 40$μ$Jy, we obtain upper limits on the differential number counts of $\rm dN/d(Log_{10}S) < 7\times10^5~deg^2$ and $\rm dN/d(Log_{10}S)<3\times10^5~deg^2$, respectively. Our results provide a new lower limit to CIB intensity of 17.2${\rm Jy\ deg^{-2}}$ at 1.1mm and of 12.9${\rm Jy\ deg^{-2}}$ at 1.3mm. Moreover, the flattening of the integrated number counts at faint fluxes strongly suggests that we are probably close to the CIB intensity. Our data imply that galaxies with SFR$<40~M_{\odot}/yr$ certainly contribute less than 50% to the CIB while more than 50% of the CIB must be produced by galaxies with $\rm SFR>40~M_{\odot}/yr$. The differential number counts are in nice agreement with recent semi-analytical models of galaxy formation even as low as our faint fluxes. Consequently, this supports the galaxy evolutionary scenarios and assumptions made in these models.

preprint2015arXiv

An extreme [OIII] emitter at $z=3.2$: a low metallicity Lyman continuum source

[Abridged] We investigate the physical properties of a Lyman continuum emitter candidate at $z=3.212$ with photometric coverage from $U$ to MIPS 24$μ$m band and VIMOS/VLT and MOSFIRE/Keck spectroscopy. Investigation of the UV spectrum confirms a direct spectroscopic detection of the Lyman continuum emission with $S/N>5$. Non-zero Ly$α$ flux at the systemic redshift and high Lyman-$α$ escape fraction suggest a low HI column density. The weak C and Si low-ionization absorption lines are also consistent with a low covering fraction along the line of sight. The [OIII]$λ\lambda4959,5007+\mathrm{H}β$ equivalent width is one of the largest reported for a galaxy at $z>3$ ($\mathrm{EW}([\mathrm{OIII}]λ\lambda4959,5007+\mathrm{H}β) \simeq 1600Å$, rest-frame) and the NIR spectrum shows that this is mainly due to an extremely strong [OIII] emission. The large observed [OIII]/[OII] ratio ($>10$) and high ionization parameter are consistent with prediction from photoionization models in case of a density-bounded nebula scenario. Furthermore, the $\mathrm{EW}([\mathrm{OIII}]λ\lambda4959,5007+\mathrm{H}β)$ is comparable to recent measurements reported at $z\sim7-9$, in the reionization epoch. We also investigate the possibility of an AGN contribution to explain the ionizing emission but most of the AGN identification diagnostics suggest that stellar emission dominates instead. This source is currently the first high-$z$ example of a Lyman continuum emitter exhibiting indirect and direct evidences of a Lyman continuum leakage and having physical properties consistent with theoretical expectation from Lyman continuum emission from a density-bounded nebula.

preprint2015arXiv

Discovering extremely compact and metal-poor, star-forming dwarf galaxies out to z ~ 0.9 in the VIMOS Ultra-Deep Survey

We report the discovery of 31 low-luminosity (-14.5 > M_{AB}(B) > -18.8), extreme emission line galaxies (EELGs) at 0.2 < z < 0.9 identified by their unusually high rest-frame equivalent widths (100 < EW[OIII] < 1700 A) as part of the VIMOS Ultra Deep Survey (VUDS). VIMOS optical spectra of unprecedented sensitivity ($I_{AB}$ ~ 25 mag) along with multiwavelength photometry and HST imaging are used to investigate spectrophotometric properties of this unique sample and explore, for the first time, the very low stellar mass end (M* < 10^8 M$_{\odot}$) of the luminosity-metallicity (LZR) and mass-metallicity (MZR) relations at z < 1. Characterized by their extreme compactness (R50 < 1 kpc), low stellar mass and enhanced specific star formation rates (SFR/M* ~ 10^{-9} - 10^{-7} yr^{-1}), the VUDS EELGs are blue dwarf galaxies likely experiencing the first stages of a vigorous galaxy-wide starburst. Using T_e-sensitive direct and strong-line methods, we find that VUDS EELGs are low-metallicity (7.5 < 12+log(O/H) < 8.3) galaxies with high ionization conditions, including at least three EELGs showing HeII 4686A emission and four EELGs of extremely metal-poor (<10% solar) galaxies. The LZR and MZR followed by EELGs show relatively large scatter, being broadly consistent with the extrapolation toward low luminosity and mass from previous studies at similar redshift. However, we find evidences that galaxies with younger and more vigorous star formation -- as characterized by their larger EWs, ionization and sSFR -- tend to be more metal-poor at a given stellar mass.

preprint2015arXiv

Electromagnetic follow-up of gravitational wave candidates: perspectives in INAF

The electromagnetic (EM) emission associated with a gravitational wave (GW) signal is one of the main goal of future astronomy. Merger of neutron stars and/or black holes and core-collapse of massive stars are expected to cause rapid transient electromagnetic signals. The EM follow-up of GW signals will have to deal with large position uncertainties. The gravitational sky localization is expected to be tens to hundreds of square degrees. Wide-field cameras and rapid follow-up observations will be crucial to characterize the EM candidates for the first EM counterpart identification. We present some of the activities that we are currently carrying on to optimize the response of the INAF network of facilities to expected GW triggers. The INAF network will represent an efficient operational framework capable of fast reaction on large error box triggers and direct identification and characterization of the candidates.

preprint2015arXiv

Faint AGNs at z>4 in the CANDELS GOODS-S field: looking for contributors to the reionization of the Universe

In order to derive the AGN contribution to the cosmological ionizing emissivity we have selected faint AGN candidates at $z>4$ in the CANDELS GOODS-South field which is one of the deepest fields with extensive multiwavelength coverage from Chandra, HST, Spitzer and various groundbased telescopes. We have adopted a relatively novel criterion. As a first step high redshift galaxies are selected in the NIR $H$ band down to very faint levels ($H\leq27$) using reliable photometric redshifts. This corresponds at $z>4$ to a selection criterion based on the galaxy rest-frame UV flux. AGN candidates are then picked up from this parent sample if they show X-ray fluxes above a threshold of $F_X\sim 1.5\times 10^{-17}$ cgs (0.5-2 keV). We have found 22 AGN candidates at $z>4$ and we have derived the first estimate of the UV luminosity function in the redshift interval $4<z<6.5$ and absolute magnitude interval $-22.5\lesssim M_{1450} \lesssim -18.5$ typical of local Seyfert galaxies. The faint end of the derived luminosity function is about two/four magnitudes fainter at $z\sim 4-6$ than that derived from previous UV surveys. We have then estimated ionizing emissivities and hydrogen photoionization rates in the same redshift interval under reasonable assumptions and after discussion of possible caveats, the most important being the large uncertainties involved in the estimate of photometric redshift for sources with featureless, almost power-law SEDs and/or low average escape fraction of ionizing photons from the AGN host galaxies. We argue that, under reasonable evaluations of possible biases, the probed AGN population can produce at $z=4-6.5$ photoionization rates consistent with that required to keep highly ionized the intergalactic medium observed in the Lyman-$α$ forest of high redshift QSO spectra, providing an important contribution to the cosmic reionization.

preprint2015arXiv

Identification of z~>2 Herschel 500 micron sources using color-deconfusion

We present a new method to search for candidate z~>2 Herschel 500μm sources in the GOODS-North field, using a S500μm/S24μm "color deconfusion" technique. Potential high-z sources are selected against low-redshift ones from their large 500μm to 24μm flux density ratios. By effectively reducing the contribution from low-redshift populations to the observed 500μm emission, we are able to identify counterparts to high-z 500μm sources whose 24μm fluxes are relatively faint. The recovery of known z~4 starbursts confirms the efficiency of this approach in selecting high-z Herschel sources. The resulting sample consists of 34 dusty star-forming galaxies at z~>2. The inferred infrared luminosities are in the range 1.5x10^12-1.8x10^13 Lsun, corresponding to dust-obscured star formation rates (SFRs) of ~260-3100 Msun/yr for a Salpeter IMF. Comparison with previous SCUBA 850μm-selected galaxy samples shows that our method is more efficient at selecting high-z dusty galaxies with a median redshift of z=3.07+/-0.83 and 10 of the sources at z~>4. We find that at a fixed luminosity, the dust temperature is ~5K cooler than that expected from the Td-LIR relation at z<1, though different temperature selection effects should be taken into account. The radio-detected subsample (excluding three strong AGN) follows the far-infrared/radio correlation at lower redshifts, and no evolution with redshift is observed out to z~5, suggesting that the far-infrared emission is star formation dominated. The contribution of the high-z Herschel 500μm sources to the cosmic SFR density is comparable to that of SMG populations at z~2.5 and at least 40% of the extinction-corrected UV samples at z~4 (abridged).

preprint2015arXiv

Infrared color selection of massive galaxies at z > 3

We introduce a new color-selection technique to identify high-redshift, massive galaxies that are systematically missed by Lyman-break selection. The new selection is based on the H_{160} and IRAC 4.5um bands, specifically H - [4.5] > 2.25 mag. These galaxies, dubbed "HIEROs", include two major populations that can be separated with an additional J - H color. The populations are massive and dusty star-forming galaxies at z > 3 (JH-blue) and extremely dusty galaxies at z < 3 (JH-red). The 350 arcmin^2 of the GOODS-N and GOODS-S fields with the deepest HST/WFC3 and IRAC data contain 285 HIEROs down to [4.5] < 24 mag. We focus here primarily on JH-blue (z > 3) HIEROs, which have a median photometric redshift z ~4.4 and stellar massM_{*}~10^{10.6} Msun, and are much fainter in the rest-frame UV than similarly massive Lyman-break galaxies (LBGs). Their star formation rates (SFRs) reaches ~240 Msun yr^{-1} leading to a specific SFR, sSFR ~4.2 Gyr^{-1}, suggesting that the sSFRs for massive galaxies continue to grow at z > 2 but at a lower growth rate than from z=0 to z=2. With a median half-light radius of 2 kpc, including ~20% as compact as quiescent galaxies at similar redshifts, JH-blue HIEROs represent perfect star-forming progenitors of the most massive (M_{*} > 10^{11.2} Msun) compact quiescent galaxies at z ~ 3 and have the right number density. HIEROs make up ~60% of all galaxies with M_{*} > 10^{10.5} Msun identified at z > 3 from their photometric redshifts. This is five times more than LBGs with nearly no overlap between the two populations. While HIEROs make up 15-25% of the total SFR density at z ~ 4-5, they completely dominate the SFR density taking place in M_{*} >10^{10.5} Msun galaxies, and are therefore crucial to understanding the very early phase of massive galaxy formation.

preprint2015arXiv

New observations of z~7 galaxies: evidence for a patchy reionization

We present new results from our search for z~7 galaxies from deep spectroscopic observations of candidate z-dropouts in the CANDELS fields. Despite the extremely low flux limits achieved by our sensitive observations, only 2 galaxies have robust redshift identifications, one from its Lyalpha emission line at z=6.65, the other from its Lyman-break, i.e. the continuum discontinuity at the Lyalpha wavelength consistent with a redshift 6.42, but with no emission line. In addition, for 23 galaxies we present deep limits in the Lyalpha EW derived from the non detections in ultra-deep observations. Using this new data as well as previous samples, we assemble a total of 68 candidate z~7 galaxies with deep spectroscopic observations, of which 12 have a line detection. With this much enlarged sample we can place solid constraints on the declining fraction of Ly$α$ emission in z~7 Lyman break galaxies compared to z~6, both for bright and faint galaxies. Applying a simple analytical model, we show that the present data favor a patchy reionization process rather than a smooth one.

preprint2015arXiv

Peering through the holes: the far UV color of star-forming galaxies at z~3-4 and the escaping fraction of ionizing radiation

We aim to investigate the effect of the escaping ionizing radiation on the color selection of high redshift galaxies and identify candidate Lyman continuum (LyC) emitters. The intergalactic medium prescription of Inoue et al.(2014) and galaxy synthesis models of Bruzual&Charlot (2003) have been used to properly treat the ultraviolet stellar emission, the stochasticity of the intergalactic transmission and mean free path in the ionizing regime. Color tracks are computed by turning on/off the escape fraction of ionizing radiation. At variance with recent studies, a careful treatment of IGM transmission leads to no significant effects on the high-redshift broad-band color selection. The decreasing mean free path of ionizing photons with increasing redshift further diminishes the contribution of the LyC to broad-band colors. We also demonstrate that prominent LyC sources can be selected under suitable conditions by calculating the probability of a null escaping ionizing radiation. The method is applied to a sample of galaxies extracted from the GOODS-S field. A known LyC source at z=3.795 is successfully recovered as a LyC emitter candidate and another convincing candidate at z=3.212 is reported. A detailed analysis of the two sources (including their variability and morphology) suggests a possible mixture of stellar and non-stellar (AGN) contribution in the ultraviolet. Conclusions: Classical broad-band color selection of 2.5<z<4.5 galaxies does not prevent the inclusion of LyC emitters in the selected samples. Large fesc in relatively bright galaxies (L>0.1L*) could be favored by the presence of a faint AGN not easily detected at any wavelength. A hybrid stellar and non-stellar (AGN) ionizing emission could coexist in these systems and explain the tensions found among the UV excess and the stellar population synthesis models reported in literature.

preprint2015arXiv

Stellar masses from the CANDELS survey: the GOODS-South and UDS fields

We present the public release of the stellar mass catalogs for the GOODS-S and UDS fields obtained using some of the deepest near-IR images available, achieved as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) project. We combine the effort from ten different teams, who computed the stellar masses using the same photometry and the same redshifts. Each team adopted their preferred fitting code, assumptions, priors, and parameter grid. The combination of results using the same underlying stellar isochrones reduces the systematics associated with the fitting code and other choices. Thanks to the availability of different estimates, we can test the effect of some specific parameters and assumptions on the stellar mass estimate. The choice of the stellar isochrone library turns out to have the largest effect on the galaxy stellar mass estimates, resulting in the largest distributions around the median value (with a semi interquartile range larger than 0.1 dex). On the other hand, for most galaxies, the stellar mass estimates are relatively insensitive to the different parameterizations of the star formation history. The inclusion of nebular emission in the model spectra does not have a significant impact for the majority of galaxies (less than a factor of 2 for ~80% of the sample). Nevertheless, the stellar mass for the subsample of young galaxies (age < 100 Myr), especially in particular redshift ranges (e.g., 2.2 < z < 2.4, 3.2 < z < 3.6, and 5.5 < z < 6.5), can be seriously overestimated (by up to a factor of 10 for < 20 Myr sources) if nebular contribution is ignored.

preprint2015arXiv

T-PHOT: A new code for PSF-matched, prior-based, multiwavelength extragalactic deconfusion photometry

We present T-PHOT, a publicly available software aimed at extracting accurate photometry from low-resolution images of deep extragalactic fields, where the blending of sources can be a serious problem for the accurate and unbiased measurement of fluxes and colours. T-PHOT has been developed within the ASTRODEEP project and it can be considered as the next generation to TFIT, providing significant improvements above it and other similar codes. T-PHOT gathers data from a high-resolution image of a region of the sky, and uses it to obtain priors for the photometric analysis of a lower resolution image of the same field. It can handle different types of datasets as input priors: i) a list of objects that will be used to obtain cutouts from the real high-resolution image; ii) a set of analytical models; iii) a list of unresolved, point-like sources, useful e.g. for far-infrared wavelength domains. We show that T-PHOT yields accurate estimations of fluxes within the intrinsic uncertainties of the method, when systematic errors are taken into account (which can be done thanks to a flagging code given in the output). T-PHOT is many times faster than similar codes like TFIT and CONVPHOT (up to hundreds, depending on the problem and the method adopted), whilst at the same time being more robust and more versatile. This makes it an optimal choice for the analysis of large datasets. In addition we show how the use of different settings and methods significantly enhances the performance. Given its versatility and robustness, T-PHOT can be considered the preferred choice for combined photometric analysis of current and forthcoming extragalactic optical to far-infrared imaging surveys. [abridged]

preprint2015arXiv

The assembly of "normal" galaxies at z=7 probed by ALMA

We report new deep ALMA observations aimed at investigating the [CII]158um line and continuum emission in three spectroscopically confirmed Lyman Break Galaxies at 6.8<z<7.1, i.e. well within the re-ionization epoch. With Star Formation Rates of SFR ~ 5-15 Msun/yr these systems are much more representative of the high-z galaxy population than other systems targeted in the past by millimeter observations. For the galaxy with the deepest observation we detect [CII] emission at redshift z=7.107, fully consistent with the Lyalpha redshift, but spatially offset by 0.7" (4 kpc) from the optical emission. At the location of the optical emission, tracing both the Lyalpha line and the far-UV continuum, no [CII] emission is detected in any of the three galaxies, with 3sigma upper limits significantly lower than the [CII] emission observed in lower reshift galaxies. These results suggest that molecular clouds in the central parts of primordial galaxies are rapidly disrupted by stellar feedback. As a result, [CII] emission mostly arises from more external accreting/satellite clumps of neutral gas. These findings are in agreement with recent models of galaxy formation. Thermal far-infrared continuum is not detected in any of the three galaxies. However, the upper limits on the infrared-to-UV emission ratio do not exceed those derived in metal- and dust-poor galaxies.

preprint2015arXiv

The Grism Lens-Amplified Survey from Space (GLASS). IV. Mass reconstruction of the lensing cluster Abell 2744 from frontier field imaging and GLASS spectroscopy

We present a strong and weak lensing reconstruction of the massive cluster Abell 2744, the first cluster for which deep Hubble Frontier Field (HFF) images and spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS) are available. By performing a targeted search for emission lines in multiply imaged sources using the GLASS spectra, we obtain 5 high-confidence spectroscopic redshifts and 2 tentative ones. We confirm 1 strongly lensed system by detecting the same emission lines in all 3 multiple images. We also search for additional line emitters blindly and use the full GLASS spectroscopic catalog to test reliability of photometric redshifts for faint line emitters. We see a reasonable agreement between our photometric and spectroscopic redshift measurements, when including nebular emission in photometric redshift estimations. We introduce a stringent procedure to identify only secure multiple image sets based on colors, morphology, and spectroscopy. By combining 7 multiple image systems with secure spectroscopic redshifts (at 5 distinct redshift planes) with 18 multiple image systems with secure photometric redshifts, we reconstruct the gravitational potential of the cluster pixellated on an adaptive grid, using a total of 72 images. The resulting mass map is compared with a stellar mass map obtained from the deep Spitzer Frontier Fields data to study the relative distribution of stars and dark matter in the cluster. We find that the stellar to total mass ratio varies substantially across the cluster field, suggesting that stars do not trace exactly the total mass in this interacting system. The maps of convergence, shear, and magnification are made available in the standard HFF format.

preprint2015arXiv

The Lyman Continuum escape fraction of galaxies at z=3.3 in the VUDS-LBC/COSMOS field

The Lyman continuum (LyC) flux escaping from high-z galaxies into the IGM is a fundamental quantity to understand the physical processes involved in the reionization epoch. We have investigated a sample of star-forming galaxies at z~3.3 in order to search for possible detections of LyC photons escaping from galaxy halos. UV deep imaging in the COSMOS field obtained with the prime focus camera LBC at the LBT telescope was used together with a catalog of spectroscopic redshifts obtained by the VIMOS Ultra Deep Survey (VUDS) to build a sample of 45 galaxies at z~3.3 with L>0.5L*. We obtained deep LBC images of galaxies with spectroscopic redshifts in the interval 3.27<z<3.40 both in the R and deep U bands. A sub-sample of 10 galaxies apparently shows escape fractions>28% but a detailed analysis of their properties reveals that, with the exception of two marginal detections (S/N~2) in the U band, all the other 8 galaxies are most likely contaminated by the UV flux of low-z interlopers located close to the high-z targets. The average escape fraction derived from the stacking of the cleaned sample was constrained to fesc_rel<2%. The implied HI photo-ionization rate is a factor two lower than that needed to keep the IGM ionized at z~3, as observed in the Lyman forest of high-z QSO spectra or by the proximity effect. These results support a scenario where high redshift, relatively bright (L>0.5L*) star-forming galaxies alone are unable to sustain the level of ionization observed in the cosmic IGM at z~3. Star-forming galaxies at higher redshift and at fainter luminosities (L<<L*) can be the major contributors to the reionization of the Universe only if their physical properties are subject to rapid changes from z~3 to z~6-10. Alternatively, ionizing sources could be discovered looking for fainter sources among the AGN population at high-z.

preprint2014arXiv

A 52 hours VLT/FORS2 spectrum of a bright z~7 HUDF galaxy: no Ly-alpha emission

We aim to determine the redshift of GDS-1408, the most solid z~7 galaxy candidate lying in the Hubble Ultra Deep Field. We have used all the VLT spectra of GDS-1408 collected by us and two other groups with FORS2 at VLT in the last five years, for a total integration time of 52hr. The combined spectrum is the deepest ever obtained of a galaxy in the Reionization epoch. We do not detect any emission line or continuum over the whole wavelength range, up to 10100A. Based on an accurate set of simulations, we are able to put a stringent upper limit of f(Lya)<3x10^(-18) erg/s/cm2 at 3-9 sigma in the explored wavelength range, corresponding to a rest-frame equivalent width EW<9A. Combining this limit with the SED modelling we refine the redshift to be z=6.82+/- 0.1 (1-sigma). The same SED fitting indicates that GDS-1408 is relatively extinct (A1600~1) with a dust corrected star formation rate of ~ 20 Msol/yr. The comparison between the un-attenuated equivalent width predicted by the case-B recombination theory and the observed upper limit, provides a limit on the effective Lya escape fraction of f_(esc)^(eff)(Lya)<8%. Even though we cannot rule out a major contribution of the inter/circum galactic medium in damping the line, a plausible interpretation is that GDS-1408 is moderately evolved and contains sufficient gas and dust to attenuate the Lya emission, before it reaches the intergalactic medium. The redshift confirmation of even the best z~7 candidates is very hard to achieve (unless the Lya or unusually strong rest-UV nebular emission lines are present) with the current generation of 8-10m class telescopes. We show that both JWST and E-ELT will be necessary to make decisive progresses. Currently, the increased redshift accuracy obtained with this kind of analysis makes ALMA an interesting option for the redshift confirmation.

preprint2014arXiv

A mass threshold in the number density of passive galaxies at z$\sim$2

The process that quenched star formation in galaxies at intermediate and high redshift is still the subject of considerable debate. One way to investigate this puzzling issue is to study the number density of quiescent galaxies at z~2, and its dependence on mass. Here we present the results of a new study based on very deep Ks-band imaging (with the HAWK-I instrument on the VLT) of two HST CANDELS fields (the UKIDSS Ultra-deep survey (UDS) field and GOODS-South). The new HAWK-I data (taken as part of the HUGS VLT Large Program) reach detection limits of Ks>26 (AB mag). We select a sample of passively-evolving galaxies in the redshift range 1.4<z<2.5. Thanks to the depth and large area coverage of our imaging, we have been able to extend the selection of quiescent galaxies a magnitude fainter than previous analyses. Through extensive simulations we demonstrate, for the first time, that the observed turn-over in the number of quiescent galaxies at K>22 is real. This has enabled us to establish unambiguously that the number counts of quiescent galaxies at z~2 flatten and slightly decline at magnitudes fainter than Ks~22(AB mag.). We show that this trend corresponds to a stellar mass threshold $M_*10^{10.8}\,{\rm M_{\odot}}$ below which the mechanism that halts the star formation in high-redshift galaxies seems to be inefficient. Finally we compare the observed pBzK number counts with those of quiescent galaxies extracted from four different semi-analytic models. We find that none of the models provides a statistically acceptable description of the number density of quiescent galaxies at these redshifts. We conclude that the mass function of quiescent galaxies as a function of redshift continues to present a key and demanding challenge for proposed models of galaxy formation and evolution.

preprint2014arXiv

Characterizing faint galaxies in the reionization epoch: LBT confirms two L<0.2L* sources at z=6.4 behind the CLASH/Frontier Fields cluster MACS0717.5+3745

We report the LBT/MODS1 spectroscopic confirmation of two images of faint Lyman alpha emitters at $z=6.4$ behind the Frontier Fields galaxy cluster MACSJ0717.5+3745. A wide range of lens models suggests that the two images are highly magnified, with a strong lower limit of mu>5. These are the faintest z>6 candidates spectroscopically confirmed to date. These may be also multiple images of the same z=6.4 source as supported by their similar intrinsic properties, but the lens models are inconclusive regarding this interpretation. To be cautious, we derive the physical properties of each image individually. Thanks to the high magnification, the observed near-infrared (restframe ultraviolet) part of the spectral energy distributions and Ly-alpha lines are well detected with S/N(m_1500)>~10 and S/N(Ly-alpha)~10-15. Adopting mu>5, the absolute magnitudes, M_1500, and Ly-alpha fluxes, are fainter than -18.7 and 2.8x10^(-18)erg/s/cm2, respectively. We find a very steep ultraviolet spectral slope beta=-3.0+/-0.5 (F_lambda=lambda^(beta)), implying that these are very young, dust-free and low metallicity objects, made of standard stellar populations or even extremely metal poor stars (age<~30Myr, E(B-V)=0 and metallicity 0.0-0.2 Z/Zsolar). The objects are compact (< 1 kpc^(2)), and with a stellar mass M* < 10^(8) M_solar. The very steep beta, the presence of the Ly-alpha line and the intrinsic FWHM (<300 km/s) of these newborn objects do not exclude a possible leakage of ionizing radiation. We discuss the possibility that such faint galaxies may resemble those responsible for cosmic reionization.

preprint2014arXiv

Constraints on the star-formation rate of z~3 LBGs with measured metallicity in the CANDELS GOODS-South field

We analyse 14 LBGs at z~2.8-3.8 constituting the only sample where both a spectroscopic measurement of their metallicity and deep IR observations (CANDELS+HUGS survey) are available. Fixing the metallicity of population synthesis models to the observed values, we determine best-fit physical parameters under different assumptions about the star-formation history and also consider the effect of nebular emission. For comparison we determine the UV slope of the objects, and use it to estimate their SFR_UV99 by correcting the UV luminosity following Meurer et al. (1999). A comparison between SFR obtained through SED-fitting (SFR_fit) and the SFR_UV99 shows that the latter are underestimated by a factor 2-10, regardless of the assumed SFH. Other SFR indicators (radio, far-IR, X-ray, recombination lines) coherently indicate SFRs a factor of 2-4 larger than SFR_UV99 and in closer agreement with SFR_fit. This discrepancy is due to the solar metallicity implied by the usual beta-A1600 conversion factor. We propose a refined relation, appropriate for sub-solar metallicity LBGs: A1600 = 5.32+1.99beta. This relation reconciles the dust-corrected UV with the SED-fitting and the other SFR indicators. We show that the fact that z~3 galaxies have sub-solar metallicity implies an upward revision by a factor of ~1.5-2 of the global SFRD, depending on the assumptions about the age of the stellar populations. We find very young best-fit ages (10-500 Myrs) for all our objects. From a careful examination of the uncertainties in the fit and the amplitude of the Balmer break we conclude that there is little evidence of the presence of old stellar population in at least half of the LBGs in our sample, suggesting that these objects are probably caught during a huge star-formation burst, rather than being the result of a smooth evolution.

preprint2014arXiv

Deep R-band counts of z~3 Lyman break galaxy candidates with the LBT

Aims. We present a deep multiwavelength imaging survey (UGR) in 3 different fields, Q0933, Q1623, and COSMOS, for a total area of ~1500arcmin^2. The data were obtained with the Large Binocular Camera on the Large Binocular Telescope. Methods. To select our Lyman break galaxy (LBG) candidates, we adopted the well established and widely used color-selection criterion (U-G vs. G-R). One of the main advantages of our survey is that it has a wider dynamic color range for U-dropout selection than in previous studies. This allows us to fully exploit the depth of our R-band images, obtaining a robust sample with few interlopers. In addition, for 2 of our fields we have spectroscopic redshift information that is needed to better estimate the completeness of our sample and interloper fraction. Results. Our limiting magnitudes reach 27.0(AB) in the R band (5σ) and 28.6(AB) in the U band (1σ). This dataset was used to derive LBG candidates at z~3. We obtained a catalog with a total of 12264 sources down to the 50% completeness magnitude limit in the R band for each field. We find a surface density of ~3 LBG candidates arcmin^2 down to R=25.5, where completeness is >=95% for all 3 fields. This number is higher than the original studies, but consistent with more recent samples.

preprint2014arXiv

Discovery of a rich proto-cluster at z=2.9 and associated diffuse cold gas in the VIMOS Ultra-Deep Survey (VUDS)

[Abridged] We characterise a massive proto-cluster at z=2.895 that we found in the COSMOS field using the spectroscopic sample of the VIMOS Ultra-Deep Survey (VUDS). This is one of the rare structures at z~3 not identified around AGNs or radio galaxies, so it is an ideal laboratory to study galaxy formation in dense environments. The structure comprises 12 galaxies with secure spectroscopic redshift in an area of 7'x8', in a z bin of Dz=0.016. The measured galaxy number overdensity is delta_g=12+/-2. This overdensity has total mass of M~8.1x10^(14)M_sun in a volume of 13x15x17 Mpc^3. Simulations indicate that such an overdensity at z~2.9 is a proto-cluster that will collapse in a cluster of total mass M~2.5x10^(15)M_sun at z=0. We compare the properties of the galaxies within the overdensity with a control sample at the same z but outside the overdensity. We did not find any statistically significant difference between the properties (stellar mass, SFR, sSFR, NUV-r, r-K) of the galaxies inside and outside the overdensity. The stacked spectrum of galaxies in the overdensity background shows a significant absorption feature at the wavelength of Lya redshifted at z=2.895 (lambda=4736 A), with a rest frame EW = 4+/- 1.4 A. Stacking only background galaxies without intervening sources at z~2.9 along their line of sight, we find that this absorption feature has a rest frame EW of 10.8+/-3.7 A, with a detection S/N of ~4. These EW values imply a high column density (N(HI)~3-20x10^(19)cm^(-2)), consistent with a scenario where such absorption is due to intervening cold gas streams, falling into the halo potential wells of the proto-cluster galaxies. However, we cannot exclude the hypothesis that this absorption is due to the diffuse gas within the overdensity.

preprint2014arXiv

Evidence of very low metallicity and high ionization state in a strongly lensed, star-forming dwarf galaxy at z=3.417

We investigate the gas-phase metallicity and Lyman Continuum (LyC) escape fraction of a strongly gravitationally lensed, extreme emission-line galaxy at z=3.417, J1000+0221S, recently discovered by the CANDELS team. We derive ionization and metallicity sensitive emission-line ratios from H+K band LBT/LUCI medium resolution spectroscopy. J1000+0221S shows high ionization conditions, as evidenced by its enhanced [OIII]/[OII] and [OIII]/Hbeta ratios. Consistently, strong-line methods based on the available line ratios suggest that J1000+0221S is an extremely metal-poor galaxy, with a metallicity of 12+log(O/H) < 7.44 (< 5% solar), placing it among the most metal-poor star-forming galaxies at z > 3 discovered so far. In combination with its low stellar mass (2x10^8 Msun) and high star formation rate (5 Msun/yr), the metallicity of J1000+0221S is consistent with the extrapolation to low masses of the mass-metallicity relation traced by Lyman-break galaxies at z > 3, but it is 0.55 dex lower than predicted by the fundamental metallicity relation at z < 2.5. These observations suggest the picture of a rapidly growing galaxy, possibly fed by the massive accretion of pristine gas. Additionally, deep LBT/LBC in the UGR bands are used to derive a limit to the LyC escape fraction, thus allowing us to explore for the first time the regime of sub-L* galaxies at z > 3. We find a 1sigma upper limit to the escape fraction of 23%, which adds a new observational constraint to recent theoretical models predicting that sub-L* galaxies at high-z have high escape fractions and thus are the responsible for the reioization of the Universe.

preprint2014arXiv

Galaxy population properties of the massive X-ray luminous galaxy cluster XDCP J0044.0-2033 at z=1.58: red-sequence formation, massive galaxy assembly, and central star formation activity

We investigate various galaxy population properties of the massive X-ray luminous galaxy cluster XDCP J0044.0-2033 at z=1.58, which constitutes the most extreme matter density peak at this redshift currently known. We analyze deep VLT/HAWK-I NIR data in the J- and Ks-bands, complemented by Subaru imaging in i and V, Spitzer observations at 4.5 micron, and new spectroscopic observations with VLT/FORS2. We detect a cluster-associated excess population of about 90 galaxies, which follows a centrally peaked, compact NFW galaxy surface density profile with a concentration of c200~10. Based on the Spitzer 4.5 micron imaging data, we measure a stellar mass fraction of fstar,500=(3.3+-1.4)% consistent with local values. The total J- and Ks-band galaxy luminosity functions of the core region yield characteristic magnitudes J* and Ks* consistent with expectations from simple z_f=3 burst models. However, a detailed look at the morphologies and color distributions of the spectroscopically confirmed members reveals that the most massive galaxies are undergoing a very active mass assembly epoch through merging processes. Consequently, the bright end of the cluster red-sequence is not in place, while at intermediate magnitudes [Ks*,Ks*+1.6] a red-locus population is present, which is then sharply truncated at magnitudes fainter than Ks*+1.6. The dominant cluster core population comprises post-quenched galaxies transitioning towards the red-sequence at intermediate magnitudes, while additionally a significant blue cloud population of faint star-forming galaxies is present even in the densest central regions. Our observations lend support to the scenario in which the dominant effect of the dense z~1.6 cluster environment is an accelerated mass assembly timescale through merging activity that is responsible for driving core galaxies across the mass quenching threshold of log(Mstar/Msun)~10.4.

preprint2014arXiv

Technical Design Report EuroGammaS proposal for the ELI-NP Gamma beam System

The machine described in this document is an advanced Source of up to 20 MeV Gamma Rays based on Compton back-scattering, i.e. collision of an intense high power laser beam and a high brightness electron beam with maximum kinetic energy of about 720 MeV. Fully equipped with collimation and characterization systems, in order to generate, form and fully measure the physical characteristics of the produced Gamma Ray beam. The quality, i.e. phase space density, of the two colliding beams will be such that the emitted Gamma ray beam is characterized by energy tunability, spectral density, bandwidth, polarization, divergence and brilliance compatible with the requested performances of the ELI-NP user facility, to be built in Romania as the Nuclear Physics oriented Pillar of the European Extreme Light Infrastructure. This document illustrates the Technical Design finally produced by the EuroGammaS Collaboration, after a thorough investigation of the machine expected performances within the constraints imposed by the ELI-NP tender for the Gamma Beam System (ELI-NP-GBS), in terms of available budget, deadlines for machine completion and performance achievement, compatibility with lay-out and characteristics of the planned civil engineering.

preprint2014arXiv

The evolution of clustering length, large-scale bias and host halo mass at 2<z<5 in the VIMOS Ultra Deep Survey (VUDS)

We investigate the evolution of galaxy clustering for galaxies in the redshift range 2.0<$z$<5.0 using the VIMOS Ultra Deep Survey (VUDS). We present the projected (real-space) two-point correlation function $w_p(r_p)$ measured by using 3022 galaxies with robust spectroscopic redshifts in two independent fields (COSMOS and VVDS-02h) covering in total 0.8 deg$^2$. We quantify how the scale dependent clustering amplitude $r_0$ changes with redshift making use of mock samples to evaluate and correct the survey selection function. Using a power-law model $ξ(r) = (r/r_0)^{-γ}$ we find that the correlation function for the general population is best fit by a model with a clustering length $r_0$=3.95$^{+0.48}_{-0.54}$ h$^{-1}$Mpc and slope $γ$=1.8$^{+0.02}_{-0.06}$ at $z$~2.5, $r_0$=4.35$\pm$0.60 h$^{-1}$Mpc and $γ$=1.6$^{+0.12}_{-0.13}$ at $z$~3.5. We use these clustering parameters to derive the large-scale linear galaxy bias $b_L^{PL}$, between galaxies and dark matter. We find $b_L^{PL}$ = 2.68$\pm$0.22 at redshift $z$~3 (assuming $σ_8$ = 0.8), significantly higher than found at intermediate and low redshifts. We fit an HOD model to the data and we obtain that the average halo mass at redshift $z$~3 is $M_h$=10$^{11.75\pm0.23}$ h$^{-1}$M$_{\odot}$. From this fit we confirm that the large-scale linear galaxy bias is relatively high at $b_L^{HOD}$ = 2.82$\pm$0.27. Comparing these measurements with similar measurements at lower redshifts we infer that the star-forming population of galaxies at $z$~3 should evolve into the massive and bright ($M_r$<-21.5) galaxy population which typically occupy haloes of mass $\langle M_h\rangle$ = 10$^{13.9}$ h$^{-1}$ $M_{\odot}$ at redshift $z$=0.

preprint2014arXiv

The evolving SFR-M_star relation and sSFR since z~5 from the VUDS spectroscopic survey

We study the evolution of the star formation rate (SFR) - stellar mass (M_star) relation and specific star formation rate (sSFR) of star forming galaxies (SFGs) since a redshift z~5.5 using 2435 (4531) galaxies with highly reliable (reliable) spectroscopic redshifts in the VIMOS Ultra-Deep Survey (VUDS). It is the first time that these relations can be followed over such a large redshift range from a single homogeneously selected sample of galaxies with spectroscopic redshifts. The log(SFR) - log(M_star) relation for SFGs remains roughly linear all the way up to z=5 but the SFR steadily increases at fixed mass with increasing redshift. We find that for stellar masses M_star>3.2 x 10^9 M_sun the SFR increases by a factor ~13 between z=0.4 and z=2.3. We extend this relation up to z=5, finding an additional increase in SFR by a factor 1.7 from z=2.3 to z=4.8 for masses M_star > 10^10 M_sun. We observe a turn-off in the SFR-M_star relation at the highest mass end up to a redshift z~3.5. We interpret this turn-off as the signature of a strong on-going quenching mechanism and rapid mass growth. The sSFR increases strongly up to z~2 but it grows much less rapidly in 2<z<5. We find that the shape of the sSFR evolution is not well reproduced by cold gas accretion-driven models or the latest hydrodynamical models. Below z~2 these models have a flatter evolution (1+z)^{Phi} with Phi=2-2.25 compared to the data which evolves more rapidly with Phi=2.8+-0.2. Above z~2, the reverse is happening with the data evolving more slowly with Phi=1.2+-0.1. The observed sSFR evolution over a large redshift range 0<z<5 and our finding of a non linear main sequence at high mass both indicate that the evolution of SFR and M_star is not solely driven by gas accretion. The results presented in this paper emphasize the need to invoke a more complex mix of physical processes {abridge}

preprint2014arXiv

The galaxy stellar mass function at 3.5<z<7.5 in the CANDELS/UDS, GOODS-South, and HUDF fields

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.5<z<7.5. We utilised the HST WFC3/IR NIR imaging from CANDELS and HUDF09, reaching H~27-28.5 over a total area of 369 arcmin2, in combination with associated deep HST ACS optical data, deep Spitzer IRAC imaging from the SEDS programme, and deep Y and K-band VLT Hawk-I images from the HUGS programme, to select a galaxy sample with high-quality photometric redshifts. These have been calibrated with more than 150 spectroscopic redshifts in the range 3.5<z<7.5, resulting in an overall precision of sigma_z/(1+z)~0.037. We have determined the low-mass end of the high-z GSMF with unprecedented precision, reaching down to masses as low as M*~10^9 Msun at z=4 and ~6x10^9 Msun at z=7. We find that the GSMF at 3.5<z<7.5 depends only slightly on the recipes adopted to measure the stellar masses, namely the photo-z, the SFHs, the nebular contribution or the presence of AGN on the parent sample. The low-mass end of the GSMF is steeper than has been found at lower redshifts, but appears to be unchanged over the redshift range probed here. Our results are very different from previous GSMF estimates based on converting UV galaxy luminosity functions into mass functions via tight M/L relations. Integrating our evolving GSMF over mass, we find that the growth of stellar mass density is barely consistent with the time-integral of the SFR density over cosmic time at z>4. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts.

preprint2014arXiv

The Hawk-I UDS and GOODS Survey (HUGS): Survey design and deep K-band number counts

We present the results of a new, ultra-deep, near-infrared imaging survey executed with the Hawk-I imager at the ESO VLT, of which we make all the data public. This survey, named HUGS (Hawk-I UDS and GOODS Survey), provides deep, high-quality imaging in the K and Y bands over the CANDELS UDS and GOODS-South fields. We describe here the survey strategy, the data reduction process, and the data quality. HUGS delivers the deepest and highest quality K-band images ever collected over areas of cosmological interest, and ideally complements the CANDELS data set in terms of image quality and depth. The seeing is exceptional and homogeneous, confined to the range 0.38"-0.43". In the deepest region of the GOODS-S field, (which includes most of the HUDF) the K-band exposure time exceeds 80 hours of integration, yielding a 1-sigma magnitude limit of ~28.0 mag/sqarcsec. In the UDS field the survey matches the shallower depth of the CANDELS images reaching a 1-sigma limit per sq.arcsec of ~27.3mag in the K band and ~28.3mag in the Y-band, We show that the HUGS observations are well matched to the depth of the CANDELS WFC3/IR data, since the majority of even the faintest galaxies detected in the CANDELS H-band images are also detected in HUGS. We present the K-band galaxy number counts produced by combining the HUGS data from the two fields. We show that the slope of the number counts depends sensitively on the assumed distribution of galaxy sizes, with potential impact on the estimated extra-galactic background light (abridged).

preprint2014arXiv

The stellar mass - halo mass relation from galaxy clustering in VUDS: a high star formation efficiency at z~3

The relation between the galaxy stellar mass M_star and the dark matter halo mass M_h gives important information on the efficiency in forming stars and assembling stellar mass in galaxies. We present the stellar mass to halo mass ratio (SMHR) measurements at redshifts 2<z<5, obtained from the VIMOS Ultra Deep Survey. We use halo occupation distribution (HOD) modelling of clustering measurements on ~3000 galaxies with spectroscopic redshifts to derive the dark matter halo mass M_h, and SED fitting over a large set of multi-wavelength data to derive the stellar mass M_star and compute the SMHR=M_star/M_h. We find that the SMHR ranges from 1% to 2.5% for galaxies with M_star=1.3x10^9 M_sun to M_star=7.4x10^9 M_sun in DM halos with M_h=1.3x10^{11} M_sun} to M_h=3x10^{11} M_sun. We derive the integrated star formation efficiency (ISFE) of these galaxies and find that the star formation efficiency is a moderate 6-9% for lower mass galaxies while it is relatively high at 16% for galaxies with the median stellar mass of the sample ~7x10^9 M_sun. The lower ISFE at lower masses may indicate that some efficient means of suppressing star formation is at work (like SNe feedback), while the high ISFE for the average galaxy at z~3 is indicating that these galaxies are efficiently building-up their stellar mass at a key epoch in the mass assembly process. We further infer that the average mass galaxy at z~3 will start experiencing star formation quenching within a few hundred millions years.

preprint2014arXiv

The VIMOS Ultra-Deep Survey (VUDS): fast increase of the fraction of strong Lyman alpha emitters from z=2 to z=6

Aims. The aim of this work is to constrain the evolution of the fraction of Lya emitters among UV selected star forming galaxies at 2<z<6, and to measure the stellar escape fraction of Lya photons over the same redshift range. Methods. We exploit the ultradeep spectroscopic observations collected by the VIMOS Ultra Deep Survey (VUDS) to build an unique, complete and unbiased sample of 4000 spectroscopically confirmed star forming galaxies at 2<z<6. Our galaxy sample UV luminosities brighter than M* at 2<z<6, and luminosities down to one magnitude fainter than M* at 2<z<3.5. Results. We find that 80% of the star forming galaxies in our sample have EW0(Lya)<10A, and correspondingly fesc(Lya)<1%. By comparing these results with literature, we conclude that the bulk of the Lya luminosity at 2<z<6 comes from galaxies that are fainter in the UV than those we sample in this work. The strong Lya emitters constitute, at each redshift, the tail of the distribution of the galaxies with extreme EW0(Lya) and fesc(Lya) . This tail of large EW0 and fesc(Lya) becomes more important as the redshift increases, and causes the fraction of Lya with EW0> 25A to increase from 5% at z=2 to 30% at z=6, with the increase being relatively stronger beyond z=4. We observe no difference, for the narrow range of UV luminosities explored in this work, between the fraction of strong Lya emitters among galaxies fainter or brighter than M*, although the fraction for the FUV faint galaxies evolves faster, at 2<z<3.5, than for the bright ones. We do observe an anticorrelation between E(B-V) and fesc(Lya): generally galaxies with high fesc(Lya) have also small amounts of dust (and viceversa). However, when the dust content is low (E(B-V)<0.05) we observe a very broad range of fesc(Lya), ranging from 10^-3 to 1. This implies that the dust alone is not the only regulator of the amount of escaping Lya photons.

preprint2014arXiv

The VIMOS Ultra-Deep Survey: ~10,000 galaxies with spectroscopic redshifts to study galaxy assembly at early epochs 2<z<~6

We present the VIMOS Ultra Deep Survey (VUDS), a spectroscopic redshift survey of ~10.000 very faint galaxies to study the major phase of galaxy assembly 2<z<~6. The survey covers 1 deg^2 in 3 separate fields: COSMOS, ECDFS and VVDS-02h, with targets selection based on an inclusive combination of photometric redshifts and color properties. Spectra covering 3650<lambda<9350 A are obtained with VIMOS on the ESO-VLT with integration times of 14h. Here we present the survey strategy, the target selection, the data processing, as well as the redshift measurement process, emphasizing the specific methods adapted to this high redshift range. The spectra quality and redshift reliability are discussed, and we derive a completeness in redshift measurement of 91%, or 74% for the most reliable measurements, down to i_AB=25, and measurements are performed all the way down to i_AB=27. The redshift distribution of the main sample peaks at z=3-4 and extends over a large redshift range mainly in 2 < z < 6. At 3<z<5, the galaxies cover a large range of luminosities -23< M_U < -20.5, stellar mass 10^9 M_sun< M_star < 10^{11} M_sun, and star formation rates 1 M_sun/yr< SFR < 10^3 M_sun/yr. We discuss the spectral properties of galaxies using individual as well as stacked spectra. The comparison between spectroscopic and photometric redshifts as well as color selection demonstrate the effectiveness of our selection scheme. With ~6000 galaxies with reliable spectroscopic redshifts in 2<z<6 expected when complete, this survey is the largest at these redshifts and offers the opportunity for unprecedented studies of the star-forming galaxy population and its distribution in large scale structures during the major phase of galaxy assembly.

preprint2014arXiv

VIMOS Ultra-Deep Survey (VUDS): Witnessing the Assembly of a Massive Cluster at z~3.3

Using new spectroscopic observations obtained as part of the VIMOS Ultra-Deep Survey (VUDS), we perform a systematic search for overdense environments in the early universe ($z>2$) and report here on the discovery of Cl J0227-0421, a massive protocluster at $z=3.29$. This protocluster is characterized by both the large overdensity of spectroscopically confirmed members, $δ_{gal}=10.5\pm2.8$, and a significant overdensity in photometric redshift members. The halo mass of this protocluster is estimated, by a variety of methods, to be roughly $3\times10^{14}$ $\mathcal{M}_{\odot}$ at $z\sim3.3$, which, evolved to $z=0$ results in a halo mass rivaling or exceeding that of the Coma cluster. The properties of 19 spectroscopically confirmed member galaxies are compared with a large sample of VUDS/VVDS galaxies in lower density field environments at similar redshifts. We find tentative evidence for an excess of redder, brighter, and more massive galaxies within the confines of the protocluster relative to the field population, which suggests that we may be observing the beginning of environmentally-induced quenching. The properties of these galaxies are investigated, including a discussion of the brightest protocluster galaxy which appears to be undergoing vigorous coeval nuclear and starburst activity. The remaining member galaxies appear to have characteristics which are largely similar to the field population. Though we find weaker evidence of the suppression of the median star formation rates amongst and differences in stacked spectra of member galaxies with respect to the field, we defer any conclusions of these trends to future work with the ensemble of protostructures that are found in the full VUDS sample.

preprint2013arXiv

Diffuse optical intracluster light as a measure of stellar tidal stripping: the cluster CL0024+17 at $z\sim $0.4 observed at LBT

We have evaluated the diffuse intracluster light (ICL) in the central core of the galaxy cluster CL0024+17 at $z\sim 0.4$ observed with the prime focus camera (LBC) at LBT. The measure required an accurate removal of the galaxies light within $\sim 200$ kpc from the center. The residual background intensity has then been integrated in circular apertures to derive the average ICL intensity profile. The latter shows an approximate exponential decline as expected from theoretical cold dark matter models. The radial profile of the ICL over the galaxies intensity ratio (ICL fraction) is increasing with decreasing radius but near the cluster center it starts to bend and then decreases where the overlap of the halos of the brightest cluster galaxies becomes dominant. Theoretical expectations in a simplified CDM scenario show that the ICL fraction profile can be estimated from the stripped over galaxy stellar mass ratio in the cluster. It is possible to show that the latter quantity is almost independent of the properties of the individual host galaxies but mainly depends on the average cluster properties. The predicted ICL fraction profile is thus very sensitive to the assumed CDM profile, total mass and concentration parameter of the cluster. Adopting values very similar to those derived from the most recent lensing analysis in CL0024+17 we find a good agreement with the observed ICL fraction profile. The galaxy counts in the cluster core have then been compared with that derived from composite cluster samples in larger volumes, up to the clusters virial radius. The galaxy counts in the CL0024+17 core appear flatter and the amount of bending respect to the average cluster galaxy counts imply a loss of total emissivity in broad agreement with the measured ICL fraction.

preprint2013arXiv

IRIDE White Book, An Interdisciplinary Research Infrastructure based on Dual Electron linacs&lasers

This report describes the scientific aims and potentials as well as the preliminary technical design of IRIDE, an innovative tool for multi-disciplinary investigations in a wide field of scientific, technological and industrial applications. IRIDE will be a high intensity 'particle factory', based on a combination of a high duty cycle radio-frequency superconducting electron linac and of high energy lasers. Conceived to provide unique research possibilities for particle physics, for condensed matter physics, chemistry and material science, for structural biology and industrial applications, IRIDE will open completely new research possibilities and advance our knowledge in many branches of science and technology. IRIDE will contribute to open new avenues of discoveries and to address most important riddles: What does matter consist of? What is the structure of proteins that have a fundamental role in life processes? What can we learn from protein structure to improve the treatment of diseases and to design more efficient drugs? But also how does an electronic chip behave under the effect of radiations? How can the heat flow in a large heat exchanger be optimized? The scientific potential of IRIDE is far reaching and justifies the construction of such a large facility in Italy in synergy with the national research institutes and companies and in the framework of the European and international research. It will impact also on R&D work for ILC, FEL, and will be complementarity to other large scale accelerator projects. IRIDE is also intended to be realized in subsequent stages of development depending on the assigned priorities.

preprint2013arXiv

Serendipitous discovery of a massive cD galaxy at z=1.096: Implications for the early formation and late evolution of cD galaxies

We have made a serendipitous discovery of a massive cD galaxy at z=1.096 in a candidate rich cluster in the HUDF area of GOODS-South. This brightest cluster galaxy is the most distant cD galaxy confirmed to date. Ultra-deep HST/WFC3 images reveal an extended envelope starting from ~10 kpc and reaching ~70 kpc in radius along the semi-major axis. The spectral energy distributions indicate that both its inner component and outer envelope are composed of an old, passively-evolving stellar population. The cD galaxy lies on the same mass-size relation as the bulk of quiescent galaxies at similar redshifts. The cD galaxy has a higher stellar mass surface density but a similar velocity dispersion to those of more-massive, nearby cDs. If the cD galaxy is one of the progenitors of today's more massive cDs, its size and stellar mass have had to increase on average by factors of $3.4\pm1.1$ and $3.3\pm1.3$ over the past ~8 Gyrs, respectively. Such increases in size and stellar mass without being accompanied by significant increases in velocity dispersion are consistent with evolutionary scenarios driven by both major and minor dry mergers. If such cD envelopes originate from dry mergers, our discovery of even one example proves that some BCGs entered the dry merger phase at epochs earlier than z=1. Our data match theoretical models which predict that the continuance of dry mergers at z<1 can result in structures similar to those of massive cD galaxies seen today. Moreover, our discovery is a surprise given that the extreme depth of the HUDF is essential to reveal such an extended cD envelope at z>1 and, yet, the HUDF covers only a minuscule region of sky. Adding that cDs are rare, Our serendipitous discovery hints that such cDs may be more common than expected. [Abridged]

preprint2013arXiv

The evolution of the AGN content in groups up to z~1

Determining the AGN content in structures of different mass/velocity dispersion and comparing them to higher mass/lower redshift analogs is important to understand how the AGN formation process is related to environmental properties. We use our well-tested cluster finding algorithm to identify structures in the GOODS North and South fields, exploiting the available spectroscopic redshifts and accurate photometric redshifts. We identify 9 structures in GOODS-south (presented in a previous paper) and 8 new structures in GOODS-north. We only consider structures where at least 2/3 of the members brighter than M_R=-20 have a spectroscopic redshift. For those group members that coincide with X-ray sources in the 4 and 2 Msec Chandra source catalogs respectively, we determine if the X-ray emission originates from AGN activity or it is related to the galaxies' star-formation activity. We find that the fraction of AGN with Log L_H > 42 erg/s in galaxies with M_R < -20 is on average 6.3+-1.3%, much higher than in lower redshift groups of similar mass and more than double the fraction found in massive clusters at a similarly high redshift. We then explore the spatial distribution of AGN in the structures and find that they preferentially populate the outer regions. The colors of AGN host galaxies in structures tend to be confined to the green valley, thus avoiding the blue cloud and, partially, also the red-sequence, contrary to what happens in the field. We finally compare our results to the predictions of two sets of semi analytic models to investigate the evolution of AGN and evaluate potential triggering and fueling mechanisms. The outcome of this comparison attests the importance of galaxy encounters, not necessarily leading to mergers, as an efficient AGN triggering mechanism. (abridged)

preprint2012arXiv

AGNs as main contributors to the UV ionizing emissivity at high redshifts: predictions from a Lambda-CDM model with linked AGN/galaxy evolution

We have evaluated the contribution of the AGN population to the ionization history of the Universe based on a semi-analytic model of galaxy formation and evolution in the CDM cosmological scenario. The model connects the growth of black holes and of the ensuing AGN activity to galaxy interactions. In the model we have included a self consistent physical description of the escape of ionizing UV photons; this is based on the blast-wave model for the AGN feedback we developed in a previous paper to explain the distribution of hydrogen column densities in AGNs of various redshifts and luminosities, due to absorption by the host galaxy gas. The model predicts UV luminosity functions for AGNs which are in good agreement with those derived from the observations especially at low and intermediate redshifts (z=3). At higher redshifts (z>5) the model tends to overestimate the data at faint luminosities. Critical biases both in the data and in the model are discussed to explain such apparent discrepancies. The predicted hydrogen photoionization rate as a function of redshift is found to be consistent with that derived from the observations. All that suggests to reconsider the role of the AGNs as the main driver of the ionization history of the Universe.

preprint2012arXiv

The blue UV slopes of z~4 Lyman break galaxies: implications for the corrected star formation rate density

We plan to analyse dust extinction in Lyman Break Galaxies (LBGs) by introducing a new and more reliable approach to their selection and to the characterization of their distribution of UV slopes beta, using deep IR images from HST. We exploit deep WFC3 IR observations of the ERS and HUDF fields over GOODS-South, combined with HST-ACS optical data, to select z~4 LBGs through a new (B-V) vs. (V-H) colour diagram. The UV slope of the selected galaxies is robustly determined by a linear fit over their observed I, Z, Y, J magnitudes, coherently with the original definition of beta. The same fit is used to determine their rest-frame UV magnitudes M1600 through a simple interpolation. We estimate the effect of observational uncertainties with detailed simulations that we also exploit, under a parametric maximum-likelihood approach, to constrain the probability density function of UV slopes PDF(beta) as a function of rest-frame magnitude. We find 142 and 25 robust LBGs in the ERS and HUDF fields respectively, limiting our sample to S/N(H)>10 objects. Our newly defined criteria improve the selection of z~4 LBGs and allow us to exclude red interlopers at lower redshift, especially z~3-3.5 objects. We find that z~4 LBGs are characterized by blue UV slopes, suggesting a low dust extinction: all L<L* galaxies have an average UV slope <beta>\simeq-2.1, while brighter objects only are slightly redder (<beta> -1.9). We find an intrinsic dispersion ~ 0.3 for PDF(beta) at all magnitudes. The SFRD at z~4 corrected according to these estimates turns out to be lower than previously found: log(SFRD)\simeq-1.09 M_sun/yr/Mpc^3. Finally, we discuss how the UV slope of z~4 galaxies changes as a function of the dust-corrected UV magnitude (i.e. SFR) showing that most galaxies with a high SFR (> 80 M_sun/yr) are highly extincted objects. [Abridged]

preprint2012arXiv

The evolving slope of the stellar mass function at 0.6 <= z < 4.5 from deep WFC3 data

We used Early Release Science (ERS) observations taken with the Wide Field Camera 3 (WFC3) in the GOODS-S field to study the galaxy stellar mass function (GSMF) at 0.6<=z<4.5. Deep WFC3 near-IR data (for Y as faint as 27.3, J and H as faint as 27.4 AB mag at 5 sigma), as well as deep Ks (as faint as 25.5 at 5 sigma) Hawk-I band data, provide an exquisite data set with which determine in an unprecedented way the low-mass end of the GSMF, allowing an accurate probe of masses as low as M~7.6 10^9 Msun at z~3. Although the area used is relatively small (~33 arcmin^2), we found generally good agreement with previous studies on the entire mass range. Our results show that the slope of the faint-end increases with redshift, from alpha=-1.44+/-0.03 at z~0.8 to alpha=-1.86+/-0.16 at z~3, although indications exist that it does not steepen further between z~3 and z~4. This result is insensitive to any uncertainty in the M* parameter. The steepness of the GSMF faint-end solves the well-known disagreement between the stellar mass density (SMD) and the integrated star formation history at z>2. However, we confirm the that there appears to be an excess of integrated star formation with respect to the SMD at z<2, by a factor of ~2-3. Our comparison of the observations with theoretical predictions shows that the models forecast a greater abundance of low mass galaxies, at least up to z~3, as well as a dearth of massive galaxies at z~4 with respect to the data, and that the predicted SMD is generally overestimated at z<~2.

preprint2012arXiv

The size-luminosity relation at z=7 in CANDELS and its implication on reionization

The exploration of the relation between galaxy sizes and other physical parameters has provided important clues for understanding galaxy formation. We use the CANDELS Deep+Wide surveys in the GOODS-South, UDS and EGS fields, complemented by data from the HUDF09 program, to address the relation between size and luminosity at z\sim7. We select 153 z-band drop-out galaxies in six different fields characterized by a wide combination of depth and areal coverage, ideally suited to sample without biases the observed size-magnitude plane. Detailed simulations allow us to derive the completeness as a function of size and magnitude and to quantify measurements errors/biases. We find that the half light radius distribution function of z\sim7 galaxies fainter than J=26.6 is peaked at <0.1 arcsec (or equivalently 0.5 kpc proper), while at brighter magnitudes high-z galaxies are typically larger than ~0.15 arcsec. We also find a well defined size-luminosity relation, Rh\simL^1/2. We compute the Luminosity Function in the HUDF and P12HUDF fields, finding large spatial variation on the number density of faint galaxies. Adopting the size distribution and the size-luminosity relation found for faint galaxies at z=7, we derive a mean slope of -1.7\pm0.1 for the luminosity function of LBGs at this redshift. Using this LF, we find that the amount of ionizing photons cannot keep the Universe re-ionized if the IGM is clumpy (C_HII>3) and the Lyman continuum escape fraction of high-z LBGs is relatively low (f_esc<0.3). With future CANDELS data, we can put severe limits to the role of galaxies in the reionization of the Universe.

preprint2011arXiv

A Compton-thick AGN at z~5 in the 4 Ms Chandra Deep Field South

We report the discovery of a Compton-thick Active Galactic Nucleus (AGN) at z=4.76 in the 4 Ms Chandra Deep Field South. This object was selected as a V-band dropout in HST/ACS images and previously recognized as an AGN from optical spectroscopy. The 4 Ms Chandra observations show a significant (~4.2sigma) X-ray detection at the V-band dropout position. The X-ray source displays a hardness ratio of HR=0.23+-0.24, which, for a source at z~5, is highly suggestive of Compton-thick absorption. The source X-ray spectrum is seen above the background level in the energy range of ~0.9-4 keV, i.e., in the rest-frame energy range of ~5-23 keV. When fixing the photon index to Gamma=1.8, the measured column density is N_H=1.4^{+0.9}_{-0.5} x 10^{24} cm^{-2}, which is Compton-thick. To our knowledge, this is the most distant heavily obscured AGN, confirmed by X-ray spectral analysis, discovered so far. The intrinsic (de-absorbed), rest-frame luminosity in the 2-10 keV band is ~2.5 x 10^{44} erg s^{-1}, which places this object among type-2 quasars. The Spectral Energy Distribution shows that massive star formation is associated with obscured black hole accretion. This system may have then been caught during a major co-eval episode of black hole and stellar mass assembly at early times. The measure of the number density of heavily obscured AGN at high redshifts will be crucial to reconstruct the BH/galaxy evolution history from the beginning.

preprint2011arXiv

A critical analysis of the UV Luminosity Function at redshift~7 from deep WFC3 data

The study of the Luminosity Function (LF) of Lyman Break Galaxies (LBGs) at z=7 is important for ascertaining their role in the reionization of the Universe. We perform a detailed and critical analysis of the statistical and systematic errors in the z~7 LF determination: we have assembled a large sample of candidate LBGs at z~7 from different surveys, spanning a large variety of areas and depths. In particular, we have combined data from the deep (J<27.4) and ultradeep (J<29.2) surveys recently acquired with the new WFC3 NIR camera on HST, over the GOODS-ERS and the HUDF fields, with ground based surveys in wide and shallow areas from VLT and Subaru. We have used public ACS images in the z-band to select z-dropout galaxies, and other public data both in the blue (BVI) and in the red bands to reject possible low-redshift interlopers. We have compared our results with extensive simulations to quantify the observational effects of our selection criteria as well as the effects of photometric scatter, color selections or the morphology of the candidates. We have found that the number density of faint LBGs at z~7 is only marginally sensitive to the color selection adopted, but it is strongly dependent from the assumption made on the half light distributions of the simulated galaxies, used to correct the observed sample for incompleteness. The slope of the faint end of the LBGs LF has thus a rather large uncertainty, due to the unknown distribution of physical sizes of the z~7 LBGs. We conclude that galaxies at z~7 are unable to reionize the Universe unless there is a significant evolution in the clumpiness of the IGM or in the escape fraction of ionising photons or, alternatively, there is a large population of z~7 LBGs with large physical dimensions but still not detected by the present observations.

preprint2011arXiv

A low escape fraction of ionizing photons of L>L* Lyman break galaxies at z=3.3

We present an upper limit for the relative escape fraction (f_{esc}^{rel}) of ionizing radiation at z~3.3 using a sample of 11 Lyman Break Galaxies (LBGs) with deep imaging in the U band obtained with the Large Binocular Camera, mounted on the prime focus of the Large Binocular Telescope. We selected 11 LBGs with secure redshift in the range 3.27<z<3.35, from 3 independent fields. We stacked the images of our sources in the R and U band, which correspond to an effective rest-frame wavelength of 1500Åand 900Årespectively, obtaining a limit in the U band image of >=30.7(AB)mag at 1 sigma. We derive a 1 sigma upper limit of f_{esc}^{rel}~5%, which is one of the lowest values found in the literature so far at z~3.3. Assuming that the upper limit for the escape fraction that we derived from our sample holds for all galaxies at this redshift, the hydrogen ionization rate that we obtain (Gamma_{-12}<0.3 s^{-1}) is not enough to keep the IGM ionized and a substantial contribution to the UV background by faint AGNs is required. Since our sample is clearly still limited in size, larger z~3 LBG samples, at similar or even greater depths are necessary to confirm these results on a more firm statistical basis.

preprint2011arXiv

An X-ray underluminous cluster of galaxies in the 4Ms CDFS observations

[Abridged] We use the large public spectroscopic database available in the GOODS-South field to estimate the dynamical mass and the virialization status of cluster ClG 0332-2747 at z=0.734. Cluster members selected from their photometric redshift are used with spectroscopic ones to analyse the galaxy population of the cluster. In the newly released Chandra 4Ms observations we detect a faint extended X-ray emission associated to the cluster. Finally, we compare the optical and X-ray properties of ClG 0332-2747 with the predictions of a well tested semianalytic model. We estimate the velocity dispersion and the virial mass considering all 44 spectroscopic members, or 20 red-sequence members only. We obtain sigma_v=634 +/- 105 Km/s, M_200=3.07 ^{+1.57}_{-1.16}~10^{14} M_sun in the former case, and slightly lower values in the latter case. The cluster appears to have reached the virial equilibrium: it shows a perfectly Gaussian velocity distribution and no evidence for substructures. ClG 0332-2747 contains a high fraction of bright red galaxies, and is dominated by a very massive (1.1 x 10^{12} M_sun) old brightest cluster galaxy (BCG), suggesting that it formed at an early epoch. We detect a faint extended X-ray source centered on the BCG, with a total X-ray luminosity of L_X ~ 2 x 10^{42} erg s^-1 (0.1-2.4 keV). This L_X is lower by a factor of ~10-20 than expected according to the M-L_X relation. We provide a possible explanation of this discrepancy as due to the effects of AGN feedback on the ICM: the semianalytic model reproduces the M-L_X relation measured from "X-ray bright" clusters, and it predicts a high scatter at low masses due to heating and expulsion of the cluster gas. Interestingly, in the model clusters with an evolved galaxy population like ClG 0332-2747 present the largest scatter in X-ray luminosity.

preprint2011arXiv

Faint high-redshift AGN in the Chandra Deep Field South: the evolution of the AGN luminosity function and black hole demography

We present detection and analysis of faint X-ray sources in the Chandra deep field south (CDFS) using the 4 Msec Chandra observation and adopting a new detection algorithm, based on a targeted search at the position of known high-z galaxies. This optimized technique results in the identification of 54 z>3 AGNs, 29 of which are new detections. Applying stringent completeness criteria, we derive AGN luminosity functions in the redshift bins 3-4, 4-5 and >5.8 and for 42.75<log L(2-10 keV)<44.5. We join this data with the luminous AGN luminosity functions from optical surveys and find that the evolution of the high-z, wide luminosity range luminosity function can be best modeled by pure luminosity evolution with L* decreasing from 6x10^44 ergs/s at z=3 to L*=2x10^44 ergs/s at z=6. We compare the high-z luminosity function with the prediction of theoretical models using galaxy interactions as AGN triggering mechanism. We find that these models are broadly able to reproduce the high-z AGN luminosity functions. A better agreement is found assuming a minimum Dark Matter halo mass for black hole formation and growth. We compare our AGN luminosity functions with galaxy mass functions to derive high-z AGN duty cycle using observed Eddington ratio distributions to derive black hole masses. We find that the duty cycle increases with galaxy stellar mass and redshift by a factor 10-30 from z=0.25 to z=4-5. We also report on the detection of a large fraction of highly obscured, Compton thick AGN at z>3 (18+17-10%). Their optical counterparts are not strongly reddened and we thus conclude that the size of the X-ray absorber is likely smaller than the dust sublimation radius. We finally report the discovery of a highly star-forming galaxy at z=3.47. If confirmed, this would be one of the farthest objects in which stellar sources are detected in X-rays.

preprint2011arXiv

Spectroscopic confirmation of two Lyman break galaxies at redshift beyond 7

We report the spectroscopic confirmation of two Lyman break galaxies at redshift > 7. The galaxies were observed as part of an ultra-deep spectroscopic campaign with FORS2 at the ESO/VLT for the confirmation of z~7 ``z--band dropout'' candidates selected from our VLT/Hawk-I imaging survey. Both galaxies show a prominent emission line at 9735A and 9858A respectively: the lines have fluxes of ~ 1.6-1.2 x 10^(-17) erg/s/cm2 and exhibit a sharp decline on the blue side and a tail on the red side. The asymmetry is quantitatively comparable to the observed asymmetry in z ~ 6 Ly-alpha lines, where absorption by neutral hydrogen in the IGM truncates the blue side of the emission line profile. We carefully evaluate the possibility that the galaxies are instead at lower redshift and we are observing either [OII], [OIII] or H-alpha emission: however from the spectroscopic and the photometric data we conclude that there are no other plausible identifications, except for Ly-alpha at redshift >7, implying that these are two of the most robust redshift determination for galaxies in the reionization epoch. Based on their redshifts and broad--band photometry, we derive limits on the star formation rate and on the ultraviolet spectral slopes of the two galaxies. We argue that these two galaxies alone are unlikely to have ionized the IGM in their surroundings.

preprint2011arXiv

Spectroscopic confirmation of z~7 LBGs: probing the earliest galaxies and the epoch of reionization

We present the final results from our ultra-deep spectroscopic campaign with FORS2 at the ESO/VLT for the confirmation of z~7 "z--band dropout" candidates selected from our VLT/Hawk-I imaging survey over three independent fields. In particular we report on two newly discovered galaxies at redshift ~6.7 in the NTT deep field: both galaxies show a Ly-alpha emission line with rest-frame EWs of the order 15-20 A and luminosities of 2-4 X 10^{42} erg/s. We also present the results of ultra-deep observations of a sample of i-dropout galaxies, from which we set a solid upper limit on the fraction of interlopers. Out of the 20 z-dropouts observed we confirm 5 galaxies at 6.7 < z < 7.1. This is systematically below the expectations drawn on the basis of lower redshift observations: in particular there is a significant lack of objects with intermediate Ly-alpha EWs (between 20 and 55 A). We conclude that the trend for the fraction of Ly-alpha emission in LBGs that is constantly increasing from z~3 to z~6 is most probably reversed from z~6 to z~7. Explaining the observed rapid change in the LAE fraction among the drop-out population with reionization requires a fast evolution of the neutral fraction of hydrogen in the Universe. Assuming that the Universe is completely ionized at z=6 and adopting the semi-analytical models of Dijkstra et al. (2011), we find that our data require a change of the neutral hydrogen fraction of the order Delta chi_{HI} ~ 0.6 in a time Delta z ~ 1, provided that the escape fraction does not increase dramatically over the same redshift interval.

preprint2010arXiv

Physical and morphological properties of z~3 LBGs: dependence on Lyalpha line emission

We investigate the physical and morphological properties of LBGs at z ~2.5 to ~3.5, to determine if and how they depend on the nature and strength of the Lyalpha emission. We selected U-dropout galaxies from the z-detected GOODS MUSIC catalog, by adapting the classical Lyman Break criteria on the GOODS filter set. We kept only those galaxies with spectroscopic confirmation, mainly from VIMOS and FORS public observations. Using the full multi-wavelength 14-bands photometry, we determined the physical properties of the galaxies, through a standard spectral energy distribution fitting with the updated Charlot & Bruzual (2009) templates. We also added other relevant observations, i.e. the 24mu m observations from Spitzer/MIPS and the 2 MSec Chandra X-ray observations. Finally, using non parametric diagnostics (Gini, Concentration, Asymmetry, M_20 and ellipticity), we characterized the rest-frame UV morphology of the galaxies. We then analyzed how these physical and morphological properties correlate with the presence of the Lyalpha line in the optical spectra. We find that, unlike at higher redshift, the dependence of physical properties on the Lyalpha line is milder: galaxies without Lyalpha in emission tend to be more massive and dustier than the rest of the sample, but all other parameters, ages, SFRs, X-ray emission as well as UV morphology do not depend strongly on the presence of the line emission. A simple scenario where all LBGs have intrinsically high Lyalpha emission, but where dust and neutral hydrogen content (which shape the final appearance of the Lyalpha) depend on the mass of the galaxies, is able to reproduce the majority of the observed properties at z~3. Some modification might be needed to account for the observed evolution of these properties with cosmic epoch, which is also discussed.

preprint2010arXiv

Proposal for taking data with the KLOE-2 detector at the DA$Φ$NE collider upgraded in energy

This document reviews the physics program of the KLOE-2 detector at DA$Φ$NE upgraded in energy and provides a simple solution to run the collider above the $ϕ$-peak (up to 2, possibly 2.5 GeV). It is shown how a precise measurement of the multihadronic cross section in the energy region up to 2 (possibly 2.5) GeV would have a major impact on the tests of the Standard Model through a precise determination of the anomalous magnetic moment of the muon and the effective fine-structure constant at the $M_Z$ scale. With a luminosity of about $10^{32}$cm$^{-2}$s$^{-1}$, DA$Φ$NE upgraded in energy can perform a scan in the region from 1 to 2.5 GeV in one year by collecting an integrated luminosity of 20 pb$^{-1}$ (corresponding to a few days of data taking) for single point, assuming an energy step of 25 MeV. A few years of data taking in this region would provide important tests of QCD and effective theories by $γγ$ physics with open thresholds for pseudo-scalar (like the $η'$), scalar ($f_0,f'_0$, etc...) and axial-vector ($a_1$, etc...) mesons; vector-mesons spectroscopy and baryon form factors; tests of CVC and searches for exotics. In the final part of the document a technical solution for the energy upgrade of DA$Φ$NE is proposed.

preprint2010arXiv

The bright end of the z ~ 7 UV Luminosity Function from a wide and deep HAWK-I survey

(Abridged) We present here the second half of an ESO Large Programme, which exploits the unique combination of area and sensitivity provided in the near-IR by the camera Hawk-I at the VLT. We have obtained - 30 observing hours with Hawk-I in the Y-band of two high galactic latitude fields. We combined the Y-band data with deep J and K Hawk-I observations, and with FORS1/FORS2 U, B, V, R, I, and Z observations to select z-drop galaxies having Z - Y > 1, no optical detection and flat Y - J and Y - K colour terms. We detect 8 high-quality candidates in the magnitude range Y = 25.5 - 26.5 that we add to the z-drop candidates selected in two Hawk-I pointings over the GOODS-South field. We use this full sample of 15 objects found in -161 arcmin^2 of our survey to constrain the average physical properties and the evolution of the number density of z ~ 7 LBGs. A stacking analysis yields a best-fit SED with photometric redshift z= 6.85 +0.20 -0.15 and an E(B-V)=0.05 +0.15 -0.05. We compute a binned estimate of the z ~ 7 LF and explore the effects of photometric scatter and model uncertainties on the statistical constraints. After accounting for the expected incompleteness through MonteCarlo simulations, we strengthen our previous finding that a Schechter luminosity function constant from z=6 to z=7 is ruled out at a >99% confidence level, even including the effects of cosmic variance. For galaxies brighter than M_1500= -19.0, we derive a luminosity density rho_UV = 1.5^{+2.1}{-0.8} x 10^25 erg/s/Hz/Mpc^3, implying a decrease by a factor 3.5 from z=6 to z=6.8. We find that, under standard assumptions, the emission rate of ionizing photons coming from UV bright galaxies is lower by at least a factor of two than the value required for reionization. Finally, we exploit deep Hawk-I J and K band observations to derive an upper limit on the number density of M1500<~ -22.0 LBGs at z-8 (Y-dropouts).

preprint2010arXiv

The lack of intense Lyman~alpha in ultradeep spectra of z=7 candidates in GOODS-S: imprint of reionization?

We present ultradeep optical spectroscopy obtained with FORS2 on VLT of seven Lyman-break galaxy (LBG) candidates at z>6.5 selected in the GOODS-S field from Hawk-I/VLT and WFC3/HST imaging. For one galaxy we detect a low significance emission line (S/N< 7), located at 9691.5 +/- 0.5A and with flux 3.4 x 10^(-18)erg/cm^2/s. If identified as Lyman alpha, it places the LBG at redshift z=6.972+/- 0.002, with a rest-frame equivalent width EW}=13A. Using Monte Carlo simulations and conservative EW distribution functions at 2<z<6, we estimate that the probability of observing no galaxies in our data with S/N>10 is ~ 2%, and that of observing only one galaxy out of seven with S/N=5 is ~4%, but these can be as small as ~1E-3, depending on the details of the EW distribution. We conclude that either a significant fraction of the candidates is not at high redshift or that some physical mechanism quenches the Lyman alpha emission emerging from the galaxies at z>6.5, abruptly reversing the trend of the increasing fraction of strong emitters with increasing redshift observed up to z~ 6.5. We discuss the possibility that an increasingly neutral intergalactic medium is responsible for such quenching.

preprint2009arXiv

A Comprehensive Study of Large Scale Structures in the GOODS-SOUTH Field up to z \sim 2.5

The aim of this paper is to identify and study the properties and galactic content of groups and clusters in the GOODS-South field up to z\sim2.5, and to analyse the physical properties of galaxies as a continuous function of environmental density up to high redshift. We use the deep (z850\sim26), multi-wavelength GOODS-MUSIC catalogue, which has a 15% of spectroscopic redshifts and accurate photometric redshifts for the remaining fraction. On these data, we apply a (2+1)D algorithm, previously developed by our group, that provides an adaptive estimate of the 3D density field. We support our analysis with simulations to evaluate the purity and the completeness of the cluster catalogue produced by our algorithm. We find several high density peaks embedded in larger structures in the redshift range 0.4-2.5. From the analysis of their physical properties (mass profile, M200, \sigmav, LX, U-B vs. B diagram), we derive that most of them are groups of galaxies, while two are poor clusters with masses of few times 10^14 Mo. For these two clusters we find, from the Chandra 2Ms data, an X-ray emission significantly lower than expected from their optical properties, suggesting that the two clusters are either not virialised or gas poor. We also analyse the dependance on environment of galaxy colours, luminosities, stellar masses, ages and star formations. We find that galaxies in high density regions are, on average, more luminous and massive than field galaxies up to z\sim 2. The fraction of red galaxies increases with luminosity and with density up to z\sim 1.2. At higher z this dependance on density disappears. The variation of galaxy properties as a function of redshift and density suggests that a significant change occurs at z\sim 1.5-2. (abridged)

preprint2009arXiv

Evidence for a fast evolution of the UV luminosity function beyond redshift 6 from a deep HAWK-I survey of the GOODS-S field

We perform a deep search for galaxies in the redshift range 6.5<z<7.5, to measure the evolution of the number density of luminous galaxies in this redshift range and derive useful constraints on the evolution of their Luminosity Function. We present here the first results of an ESO Large Program, that exploits the unique combination of area and sensitivity provided in the near-IR by the camera Hawk-I at the VLT. We have obtained two Hawk-I pointings on the GOODS South field for a total of 32 observing hours, covering ~90 arcmin2. The images reach Y=26.7 mags for the two fields. We have used public ACS images in the z band to select z-dropout galaxies with the colour criteria Z-Y>1, Y-J<1.5 and Y-K<2. The other public data in the UBVRIJHK bands are used to reject possible low redshift interlopers. The output has been compared with extensive Monte Carlo simulations to quantify the observational effects of our selection criteria as well as the effects of photometric errors. We detect 7 high quality candidates in the magnitude range Y=25.5-26.7. This interval samples the critical range for M* at z>6 (M_1500 ~- 19.5 to -21.5). After accounting for the expected incompleteness, we rule out at a 99% confidence level a Luminosity Function constant from z=6 to z=7, even including the effects of cosmic variance. For galaxies brighter than M_1500=-19.0 we derive a luminosity density rho_UV=1.5^{+2.0}_{-0.9} 10^25 erg/s/Hz/Mpc3, implying a decrease by a factor 3.5 from z=6 to z~6.8. On the basis of our findings, we make predictions for the surface densities expected in future surveys surveys, based on ULTRA-VISTA, HST-WFC3 or JWST-NIRCam, evaluating the best observational strategy to maximise their impact.

preprint2009arXiv

Star formation and mass assembly in high redshift galaxies

We study the star formation and the mass assembly process of 0.3<=z<2.5 galaxies using their IR emission from MIPS 24um band. We used an updated version of the GOODS-MUSIC catalog, extended by the addition of mid-IR fluxes. We compared two different estimators of the Star Formation Rate: the total infrared emission derived from 24um, estimated using both synthetic and empirical IR templates, and the multiwavelength fit to the full galaxy SED. For both estimates, we computed the SFR Density and the Specific SFR. The two SFR tracers are roughly consistent, given the uncertainties involved. However, they show a systematic trend, IR-based estimates exceeding the fit-based ones as the SFR increases. We show that: a) at z>0.3, the SFR is well correlated with stellar mass, and this relationship seems to steepen with redshift (using IR-based SFRs); b) the contribution to the global SFRD by massive galaxies increases with redshift up to ~2.5, more rapidly than for galaxies of lower mass, but appears to flatten at higher z; c) despite this increase, the most important contributors to the SFRD at any z are galaxies of about, or immediately lower than, the characteristic stellar mass; d) at z~2, massive galaxies are actively star-forming, with a median SFR 300 Msun/yr. During this epoch, they assemble a substantial part of their final stellar mass; e) the SSFR shows a clear bimodal distribution. The analysis of the SFRD and the SSFR seems to support the downsizing scenario, according to which high mass galaxies have formed their stars earlier and faster than their low mass counterparts. A comparison with theoretical models indicates that they follow the global increase in the SSFR with redshift and predict the existence of quiescent galaxies even at z>1.5, but they systematically underpredict the average SSFR.

preprint2009arXiv

The Fraction of Quiescent Massive Galaxies in the Early Universe

Aims: The aim of this work is to collect a complete, mass--selected sample of galaxies with very low specific star formation rate, for a comparison with the prediction of recent theoretical models. Method: We use the 24/K flux ratio, complemented by the SED fitting to the full 0.35-8.0 mum spectral distribution, to select quiescent galaxies from z~0.4 to z~4 in the GOODS--MUSIC sample. Our observational selection can be translated into thresholds on the specific star formation rate SFR/M_*, that can be used to compare with the theoretical predictions. Results: We find that, in the framework of the well known global decline of the quiescent fraction with redshift, a non-negligible fraction ~15-20% of massive galaxies with very low specific star formation rate exists up to z~4, including a tail of "Red&Dead" galaxies with SFR/M_*<10^{-11}/yr. Recent theoretical models vary to a large extent in the prediction of the fraction of galaxies with very low specific star formation rates, and are unable to provide a global match to our data.

preprint2008arXiv

Large Scale Structures at High Redshift in the GOODS Field

We present a catalogue of overdensities in the GOODS-South field. We find various high density peaks that are embedded in structures diffused on the entire field, up to z ~ 2.5. The slope of their colour-magnitude relation does not show significative evolution with z. We find evidence that galaxies forming these structures are more massive than galaxies located in low density regions. We also analyse the variation of galaxy properties with the associated environmental density and we find that the segregation of red galaxies with density is stronger at low redshift and at high luminosities while it gets much weaker for increasing z.