Researcher profile

O. Cucciati

O. Cucciati contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

15 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. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

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

preprint2026arXiv

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

preprint2023arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2021arXiv

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

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

preprint2020arXiv

Discovery of molecular gas fueling galaxy growth in a protocluster at z=1.7

Based on ALMA Band 3 observations of the CO(2-1) line transition, we report the discovery of three new gas-rich (M_H2 ~ 1.5-4.8 x 10^10 M_sun, SFRs in the range ~5-100 M_sun/yr) galaxies in an overdense region at z=1.7, that already contains eight spectroscopically confirmed members. This leads to a total of 11 confirmed overdensity members, within a projected distance of ~ 1.15 Mpc and in a redshift range of Dz = 0.012. Under simple assumptions, we estimate that the system has a total mass of >= 3-6 x 10^13 M_sun, and show that it will likely evolve into a >~ 10^14 M_sun cluster at z = 0. The overdensity includes a powerful Compton-thick Fanaroff-Riley type II (FRII) radio-galaxy, around which we discovered a large molecular gas reservoir (M_H2 ~ 2 x 10^11 M_sun). We fitted the FRII resolved CO emission with a 2-D Gaussian model with major (minor) axis of ~ 27 (~ 17) kpc, that is a factor of ~3 larger than the optical rest-frame emission. Under the assumption of a simple edge-on disk morphology, we find that the galaxy interstellar medium produces a column density towards the nucleus of ~ 5.5 x 10^23 cm^-2. Such a dense ISM may then contribute significantly to the total nuclear obscuration measured in the X-rays (N_(H,X) ~ 1.5 x 10^24 cm^-2) in addition to a small, pc-scale absorber around the central engine. The velocity map of this source unveils a rotational motion of the gas that is perpendicular to the radio-jets. The FRII is located at the center of the projected spatial distribution of the structure members, and its velocity offset from the peak of the redshift distribution is well within the structure&#39;s velocity dispersion. All this, coupled with the large amount of gas around the FRII, its stellar mass of ~ 3 x 10^11 M_sun, SFR of ~ 200-600 M_sun/yr, and powerful radio-to-X-ray emission, suggests that this source is the likely progenitor of the future brightest cluster galaxy.

preprint2020arXiv

MOONRISE: The Main MOONS GTO Extragalactic Survey

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

preprint2020arXiv

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

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.

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.