The Euclid mission -with its spectroscopic galaxy survey covering a sky area over 15 000 deg 2 in the redshift range 0.9 < z < 1.8 -will provide a sample of tens of thousands of cosmic voids. This paper thoroughly explores for the first time the constraining power of the void size function on the properties of dark energy (DE) from a survey mock catalogue, the official Euclid Flagship simulation. We identify voids in the Flagship light-cone, which closely matches the features of the upcoming Euclid spectroscopic data set. We model the void size function considering a state-of-the art methodology: we rely on the volume conserving (Vdn) model, a modification of the popular Sheth & van de Weygaert model for void number counts, extended by means of a linear function of the large-scale galaxy bias. We find an excellent agreement between model predictions and measured mock void number counts. We compute updated forecasts for the Euclid mission on DE from the void size function and provide reliable void number estimates to serve as a basis for further forecasts of cosmological applications using voids. We analyse two different cosmological models for DE: the first described by a constant DE equation of state parameter, w, and the second by a dynamic equation of state with coefficients w 0 and w a . We forecast 1σ errors on w lower than the 10%, and we estimate an expected figure of merit (FoM) for the dynamical DE scenario FoM w 0 ,wa = 17 when considering only the neutrino mass as additional free parameter of the model. The analysis is based on conservative assumptions to ensure full robustness, and is a pathfinder for future enhancements of the technique. Our results showcase the impressive constraining power of the void size function from the Euclid spectroscopic sample, both as a stand-alone probe, and to be combined with other Euclid cosmological probes.

We have investigated the cosmological constraints that can be expected from measurement of the cross-correlation of galaxies with cosmic voids identified in the Euclid spectroscopic survey, which will include spectroscopic information for tens of millions of galaxies over 15 000 deg2 of the sky in the redshift range 0.9 ≤ z < 1.8. We have done this using simulated measurements obtained from the Flagship mock catalogue, the official Euclid mock that closely matches the expected properties of the spectroscopic dataset. To mitigate anisotropic selection-bias effects, we have used a velocity field reconstruction method to remove large-scale redshift-space distortions from the galaxy field before void-finding. This allowed us to accurately model contributions to the observed anisotropy of the cross-correlation function arising from galaxy velocities around voids as well as from the Alcock–Paczynski effect, and we studied the dependence of constraints on the efficiency of reconstruction. We find that Euclid voids will be able to constrain the ratio of the transverse comoving distance DM and Hubble distance DH to a relative precision of about 0.3%, and the growth rate fσ8 to a precision of between 5% and 8% in each of the four redshift bins covering the full redshift range. In the standard cosmological model, this translates to a statistical uncertainty ΔΩm = ±0.0028 on the matter density parameter from voids, which is better than what can be achieved from either Euclid galaxy clustering and weak lensing individually. We also find that voids alone can measure the dark energy equation of state to a 6% precision.