Source author record

E. Merlin

E. Merlin 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

51works
7topics
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

51 published item(s)

preprint2026arXiv

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

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

preprint2026arXiv

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

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

preprint2026arXiv

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

We model intrinsic alignments (IA) in Euclid'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

A Lyman continuum analysis for $\sim 100$ galaxies at $z_{\text{spec}} \sim 3$ in the Abell 2744 Cluster Field

Identifying Lyman continuum (LyC) leakers at intermediate redshifts is crucial for understanding the properties of cosmic reionizers, as the opacity of the intergalactic medium (IGM) prevents direct detection of LyC emission from sources during the Epoch of Reionization (EoR). In this study, we confirm two new LyC candidate leakers at $z \sim 3$ in the Abell 2744 cluster field, with absolute escape fractions ($f_{\text{esc}}$) of $0.90^{+0.07}_{-0.86}$ and $0.60^{+0.37}_{-0.56}$, respectively. The LyC emission was detected using HST/WFC3/F275W and F336W imaging. These two candidate leakers appear faint ($ M_{\text{UV}} = -18.1 \pm 0.1 \text{ and } -17.81 \pm 0.11$), exhibit blue UV continuum slopes ($β= -2.42 \pm 0.05 \text{ and } -1.78 \pm 0.19$), have low masses ($M_\star \sim 10^{7.73} \pm 0.1 \text{ and } 10^{7.07} \pm 0.05 M_\odot$) and show \lya\ equivalent widths of $90 \pm 3$ Å and $28 \pm 12$ Å, respectively. The discovery of these two LyC candidate leakers was achieved in a catalog of 91 spectroscopically confirmed sources using JWST and/or MUSE public spectra. We also analyze properties that have been proposed as indirect indicators of LyC emission, like \lya, O32 ratio, and $M_\star$: we create subsample of galaxies selected according to such properties, stack the LyC observations of these subsample and assess the limits in escape fractions in the stacks. By analysing the individual candidates and the stacks, in the context of the currently limited sample of known LyC leakers at $z \sim 3$, we aim to enhance our understanding of LyC escape mechanisms and improve our predictions of the LyC $f_{\text{esc}}$ during the EoR.

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.

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

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)

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

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]

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

Multiple AGN activity during the BCG assembly of XDCPJ0044.0-2033 at z~1.6

Undisturbed galaxy clusters are characterized by a massive and large elliptical galaxy at their center, i.e. the Brightest Cluster Galaxy (BCG). How these central galaxies form is still debated. According to most models, a typical epoch for their assembly is z~1-2. We have performed a detailed multi-wavelength analysis of the core of XDCPJ0044.0-2033 (XDCP0044), one of the most massive and densest galaxy clusters currently known at redshift z~1.6, whose central galaxy population shows high star formation compared to lower-z clusters and an X-ray AGN located close to its center. SINFONI J-, H- and KMOS YJ-, H- bands spectroscopic data have been analyzed, together with deep archival HST photometric data in F105W, F140W, and F160W bands, Chandra X-ray, radio JVLA data at 1-2 GHz, and ALMA band-6 observations. In the central region of the cluster (~70x70 kpc^2), 2 systems of interacting galaxies have been identified and studied (Complex A and B), with a total of 7 confirmed cluster members. These galaxies show perturbed morphologies and 3 of them show signs of AGN activity. In particular, 2 type-1 AGN with typical broad lines have been found at the center of each complex (both of them X-ray obscured and highly accreting; Eddington ratio ~0.5), while a type-2 AGN has been discovered in Complex A. The AGN at the center of Complex B is also detected in X-ray while the other 2 are spatially related to radio emission. The 3 AGN provide one of the closest AGN triple at z>1 revealed so far with a minimum(maximum) projected distance of ~10(40) kpc. The observation of high star formation, merger signatures and nuclear activity in the core of XDCP0044 suggests that all these processes are key ingredients in shaping the nascent BCG. According to our data, XDCP0044 could form a typical massive galaxy of 10^12 Msun, hosting a Black Hole of 2x10^8-10^9 Msun, in a time scale of the order of ~2.5 Gyrs.

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

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

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]

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

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

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

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

preprint2012arXiv

Complexity reduction of astrochemical networks

We present a new computational scheme aimed at reducing the complexity of the chemical networks in astrophysical models, one which is shown to markedly improve their computational efficiency. It contains a flux-reduction scheme that permits to deal with both large and small systems. This procedure is shown to yield a large speed-up of the corresponding numerical codes and provides good accord with the full network results. We analyse and discuss two examples involving chemistry networks of the interstellar medium and show that the results from the present reduction technique reproduce very well the results from fuller calculations.

preprint2012arXiv

SPIDER VII - Revealing the Stellar Population Content of Massive Early-type Galaxies out to 8Re

Radial trends of stellar populations in galaxies provide a valuable tool to understand the mechanisms of galaxy growth. In this paper, we present the first comprehensive analysis of optical-optical and optical-NIR colours, as a function of galaxy mass, out to the halo region (8Re) of early-type galaxies (ETGs). We select a sample of 674 massive ETGs (M*>3x10^10MSun) from the SDSS-based SPIDER survey. By comparing with a large range of population synthesis models, we derive robust constraints on the radial trends in age and metallicity. Metallicity is unambiguously found to decrease outwards, with a measurable steepening of the slope in the outer regions (Re<R<8Re). The gradients in stellar age are found to be more sensitive to the models used, but in general, the outer regions of ETGs feature older populations compared to the cores. This trend is strongest for the most massive galaxies in our sample (M*>10^11MSun). Furthermore, when segregating with respect to large scale environment, the age gradient is more significant in ETGs residing in higher density regions. These results shed light on the processes leading from the formation of the central core to the growth of the stellar envelope of massive galaxies. The fact that the populations in the outer regions are older and more metal-poor than in the core suggests a process whereby the envelope of massive galaxies is made up of accreted small satellites (i.e. minor mergers) whose stars were born during the first stages of galaxy formation.

preprint2011arXiv

Formation and Evolution of the Dust in Galaxies. I. The Condensation Efficiencies

The growing interest in the high-z universe, where strongly obscured objects are present, has determined an effort to improve the simulations of dust formation and evolution in galaxies. Three main basic ingredients enter the problem influencing the total dust budget and the mixture of the dust grains: the types and amounts of dust injected by AGB stars and SNae and the accretion/destruction processes of dust in the ISM. They govern the relative abundances of the gas and dust in the ISM. In this study, we focus on star-dust and present a database of condensation efficiencies for the refractory elements C, O, Mg, Si, S, Ca and Fe in AGB stars and SNae that can be easily applied to the gaseous ejecta, in order to determine the amount and kind of refractory elements locally embedded into dust and injected into the ISM. The best theoretical recipes available nowadays in literature to estimate the amount of dust produced by SNae and AGB stars have been discussed and for SNae compared to the observations to get clues on the problem. The condensation efficiencies have been analyzed in the context of a classical chemical model of dust formation and evolution in the Solar Neighbourhood and Galactic Disk. Tables of coefficients are presented for (i) AGB stars at varying the metallicity and (ii) SNae at varying the density n_H of the ISM where the SNa explosions took place. In particular, we show how the controversial CNT approximation widely adopted to form dust in SNae, still gives good results and agrees with some clues coming from the observations. A new generation of dust formation models in SNae is however required to solve some contradictions recently emerged. A simple database of condensation efficiencies is set up to be used in chemical models including the effect of dust and meant to simulate galaxies going from primordial proto-galaxies to those currently seen in the local universe.

preprint2011arXiv

Formation and Evolution of the Dust in Galaxies. II. The Solar Neighbourhood

Over the past decade a new generation of chemical models have included the dust in the treatment of the ISM. This major accomplishment has been spurred by the growing amounts of data on the highly obscured high-z Universe and the intriguing local properties of the Solar Neighbourhood (SoNE). We present here a new model able to simulate the formation and evolution of dust in the ISM. The model follows the evolution of 16 elemental species, with particular attention to those that are simultaneously present in form of gas and dust, e.g. C, N, O, Mg, Si, S, Ca and Fe. In this study we focus on the SoNe and the MW Disk as a whole which are considered as laboratories to test the physical ingredients governing the dust evolution. Infall of primordial gas, birth and death of stars, radial flows of matter between contiguous shells, presence of a central bar, star-dust emission by SNae and AGB stars, dust destruction and accretion are taken into account. The model reproduces the local depletion of the elements in the gas, and simultaneously satisfies other constraints obtained from the observations. The evolution of the element abundances in the gas and dust has been well reproduced for plausible choices of the parameters. The Mg/Si ratio, in particular, drives the formation of silicates. We show that for most of the evolution of the MW, the main process for dust enrichment is the accretion in the cold regions of the ISM. SNae dominate in the early phases of the evolution. We have examined the main factors controlling the temporal window in which SNae govern the dust budget both in low and high star forming environments. The role played by AGB stars is discussed. We find that IMFs with regular slope in the range of massive stars better reproduce the observed depletions. The results obtained for the SoNe lead us to safely extend the model.

preprint2011arXiv

Formation and Evolution of the Dust in Galaxies. III. The Disk of the Milky Way

Models of chemical evolution of galaxies including the dust are nowadays required to decipher the high-z universe. In a series of three papers we have tackled the problem and set a modern chemical evolution model. In the first paper (Piovan et al., 2011a) we revised the condensation coefficients for the elements that typically are present in the dust. In the second paper (Piovan et al., 2011b) we have implemented the dust into the Padova chemical model and tested it against the observational data for the Solar Neighbourhood. In this paper we extend it to the whole Disk of the Milky Way (MW). The Disk is used as a laboratory to analyze the spatial and temporal behaviour of (i) several dust grain families with the aid of which we can describe the ISM, (ii) the abundances in the gas, dust, and total ISM of the elements present in the dust and (iii) the depletion of the same elements. The temporal evolution of the dust and gas across the Disk is calculated under the effect of radial flows and a central Bar. The gradients of the abundances of C, N, O, Mg, Si, S, Ca and Fe in gas and dust across the Disk are derived as a function of time. The theoretical gradients nicely reproduce those derived from Cepheids, OB stars, Red Giants and HII regions. This provides the backbone for the companion processes of dust formation and evolution across the Disk. We examine in detail the contributions to dust by AGB stars, SNae and grain accretion in the ISM at different galacto-centric distances. Furthermore, we examine the variation of the ratio silicates/carbonaceous grains with time and position in the Disk. Finally, some hints about the depletion of the elements in regions of high and low SFR (inner and outer Disk) are presented. The results obtained make it possible to extend the model to other astrophysical situations or different theoretical models like the chemo-dynamical N-Body simulations.

preprint2011arXiv

MaNN: Multiple Artificial Neural Networks for modelling the Interstellar Medium

Modelling the complex physics of the Interstellar Medium (ISM) in the context of large-scale numerical simulations is a challenging task. A number of methods have been proposed to embed a description of the ISM into different codes. We propose a new way to achieve this task: Artificial Neural Networks (ANNs). The ANN has been trained on a pre-compiled model database, and its predictions have been compared to the expected theoretical ones, finding good agreement both in static and in dynamical tests run using the Padova Tree-SPH code \textsc{EvoL}. A neural network can reproduce the details of the interstellar gas evolution, requiring limited computational resources. We suggest that such an algorithm can replace a real-time calculation of mass elements chemical evolution in hydrodynamical codes.

preprint2011arXiv

ROBO: a Model and a Code for the Study of the Interstellar Medium

We present ROBO, a model and its companion code for the study of the interstellar medium (ISM). The aim is to provide an accurate description of the physical evolution of the ISM and to set the ground for an ancillary tool to be inserted in NBody-Tree-SPH (NB-TSPH) simulations of large scale structures in cosmological context or of the formation and evolution of individual galaxies. The ISM model consists of gas and dust. The gas chemical composition is regulated by a network of reactions that includes a large number of species (hydrogen and deuterium based molecules, helium, and metals). New reaction rates for the charge transfer in $\mathrm H^+$ and $\mathrm H_2$ collisions are presented. The dust contains the standard mixture of carbonaceous grains (graphite grains and PAHs) and silicates of which the model follows the formation and destruction by several processes. The model takes into account an accurate treatment of the cooling process, based on several physical mechanisms, and cooling functions recently reported in the literature. The model is applied to a wide range of the input parameters and the results for important quantities describing the physical state of the gas and dust are presented. The results are organized in a database suited to the artificial neural networks (ANNs). Once trained, the ANNs yield the same results obtained by ROBO, with great accuracy. We plan to develop ANNs suitably tailored for applications to NB-TSPH simulations of cosmological structures and/or galaxies.