# Euclid: Forecasts from the void-lensing cross-correlation

M. Bonici, C. Carbone, S. Davini, P. Vielzeuf, L. Paganin, V. Cardone, N. Hamaus, A. Pisani, A.J. Hawken, A. Kovacs, et al.

► To cite this version:

M. Bonici, C. Carbone, S. Davini, P. Vielzeuf, L. Paganin, et al. Euclid: Forecasts from the void-lensing cross-correlation. Astronomy & Astrophysics - A&A, 2023, 670, pp.A47. 10.1051/0004-6361/20224445. hal-03737493

HAL Id: hal-03737493

https://hal.science/hal-03737493v1

Submitted on 10 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Euclid: Forecasts from the void-lensing cross-correlation*

M. Bonici<sup>1,2,●</sup>, C. Carbone<sup>2</sup>, S. Davini<sup>3</sup>, P. Vielzeuf<sup>4</sup>, L. Paganin<sup>1,3</sup>, V. Cardone<sup>5,6</sup>, N. Hamaus<sup>7</sup>, A. Pisani<sup>8</sup>,
A. J. Hawken<sup>4</sup>, A. Kovacs<sup>9,10</sup>, S. Nadathur<sup>11</sup>, S. Contarini<sup>12,13,14</sup>, G. Verza<sup>15,16</sup>, I. Tutusaus<sup>17,18,19,20</sup>,
F. Marulli<sup>12,13,14</sup>, L. Moscardini<sup>12,13,14</sup>, M. Aubert<sup>21</sup>, C. Giocoli<sup>22,23</sup>, A. Pourtsidou<sup>24</sup>, S. Camera<sup>25,26,27</sup>,
S. Escoffier<sup>4</sup>, A. Caminata<sup>3</sup>, S. Di Domizio<sup>1</sup>, M. Martinelli<sup>28</sup>, M. Pallavicini<sup>1,3</sup>, V. Pettorino<sup>29,30</sup>, Z. Sakr<sup>20,31</sup>,
D. Sapone<sup>32</sup>, G. Testera<sup>3</sup>, S. Tosi<sup>1</sup>, V. Yankelevich<sup>33</sup>, A. Amara<sup>34</sup>, N. Auricchio<sup>13</sup>, M. Baldi<sup>35,13,14</sup>, D. Bonino<sup>27</sup>,
E. Branchini<sup>36,37</sup>, M. Brescia<sup>38</sup>, J. Brinchmann<sup>39</sup>, V. Capobianco<sup>27</sup>, J. Carretero<sup>40,41</sup>, M. Castellano<sup>5</sup>,
S. Cavuoti<sup>38,42,43</sup>, R. Cledassou<sup>44,45</sup>, G. Congedo<sup>24</sup>, L. Conversi<sup>46,47</sup>, Y. Copin<sup>48</sup>, L. Corcione<sup>27</sup>, F. Courbin<sup>49</sup>,
M. Cropper<sup>50</sup>, A. Da Silva<sup>51,52</sup>, H. Degaudenzi<sup>53</sup>, M. Douspis<sup>54</sup>, F. Dubath<sup>53</sup>, C. A. J. Duncan<sup>55</sup>, X. Dupac<sup>47</sup>,
S. Dusini<sup>16</sup>, A. Ealet<sup>48</sup>, S. Farrens<sup>29</sup>, S. Ferriol<sup>48</sup>, P. Fosalba<sup>19,18</sup>, M. Frailis<sup>56</sup>, E. Franceschi<sup>13</sup>, M. Fumana<sup>2</sup>,
P. Gómez-Alvarez<sup>57,47</sup>, B. Garilli<sup>2</sup>, B. Gillis<sup>24</sup>, A. Grazian<sup>58</sup>, F. Grupp<sup>59,7</sup>, L. Guzzo<sup>60,61,62</sup>, S. V. H. Haugan<sup>63</sup>,
W. Holmes<sup>64</sup>, F. Hormuth<sup>65</sup>, A. Hornstrup<sup>66</sup>, K. Jahnke<sup>67</sup>, M. Kümmel<sup>7</sup>, S. Kermiche<sup>4</sup>, A. Kiessling<sup>64</sup>,
M. Kilbinger<sup>30</sup>, M. Kunz<sup>17</sup>, H. Kurki-Suonio<sup>68</sup>, R. Laureijs<sup>69</sup>, S. Ligori<sup>27</sup>, P. B. Lilje<sup>63</sup>, I. Lloro<sup>70</sup>, E. Maiorano<sup>13</sup>,
O. Mansutti<sup>56</sup>, O. Marggraf<sup>71</sup>, K. Markovic<sup>64</sup>, R. Massey<sup>72</sup>, E. Medinaceli<sup>13</sup>, M. Melchior<sup>73</sup>, M. Meneghetti<sup>13,74</sup>,
G. Meylan<sup>49</sup>, M. Moresco<sup>12,13</sup>, E. Munari<sup>56</sup>, S. M. Niemi<sup>69</sup>, C. Padilla<sup>40</sup>, S. Paltani<sup>53</sup>, F. Pasian<sup>56</sup>, K. Pedersen<sup>75</sup>,
W. J. Percival<sup>76,77,78</sup>, S. Pires<sup>29</sup>, G. Polenta<sup>79</sup>, M. Poncet<sup>44</sup>, L. Popa<sup>80</sup>, F. Raison<sup>59</sup>, R. Rebolo<sup>10,81</sup>, A. Renzi<sup>16,15</sup>,
J. Rhodes<sup>64</sup>, E. Rossetti<sup>12</sup>, R. Saglia<sup>59,7</sup>, B. Sartoris<sup>56,82</sup>, M. Scodeggio<sup>2</sup>, A. Secroun<sup>4</sup>, G. Seidel<sup>67</sup>, C. Sirignano<sup>15,16</sup>,
G. Sirri<sup>14</sup>, L. Stanco<sup>16</sup>, J.-L. Starck<sup>29</sup>, C. Surace<sup>83</sup>, P. Tallada-Crespi<sup>84,41</sup>, D. Tavagnacco<sup>56</sup>, A. N. Taylor<sup>24</sup>,
I. Tereno<sup>51,85</sup>, R. Toledo-Moreo<sup>86</sup>, F. Torradeflot<sup>84,41</sup>, E. A. Valentijn<sup>87</sup>, L. Valenziano<sup>13,14</sup>, Y. Wang<sup>88</sup>, J. Weller<sup>59,7</sup>,
G. Zamorani<sup>13</sup>, J. Zoubian<sup>4</sup>, and S. Andreon<sup>61</sup>

*(Affiliations can be found after the references)*

Received 7 July 2022 / Accepted 29 October 2022

## ABSTRACT

The *Euclid* space telescope will survey a large dataset of cosmic voids traced by dense samples of galaxies. In this work we estimate its expected performance when exploiting angular photometric void clustering, galaxy weak lensing, and their cross-correlation. To this aim, we implemented a Fisher matrix approach tailored for voids from the *Euclid* photometric dataset and we present the first forecasts on cosmological parameters that include the void-lensing correlation. We examined two different probe settings, pessimistic and optimistic, both for void clustering and galaxy lensing. We carried out forecast analyses in four model cosmologies, accounting for a varying total neutrino mass, $M_ν$, and a dynamical dark energy (DE) equation of state, $w(z)$, described by the popular Chevallier-Polarski-Linder parametrization. We find that void clustering constraints on $h$ and $\Omega_b$ are competitive with galaxy lensing alone, while errors on $n_s$ decrease thanks to the orthogonality of the two probes in the 2D-projected parameter space. We also note that, as a whole, with respect to assuming the two probes as independent, the inclusion of the void-lensing cross-correlation signal improves parameter constraints by 10–15%, and enhances the joint void clustering and galaxy lensing figure of merit (FoM) by 10% and 25%, in the pessimistic and optimistic scenarios, respectively. Finally, when further combining with the spectroscopic galaxy clustering, assumed as an independent probe, we find that, in the most competitive case, the FoM increases by a factor of 4 with respect to the combination of weak lensing and spectroscopic galaxy clustering taken as independent probes. The forecasts presented in this work show that photometric void clustering and its cross-correlation with galaxy lensing deserve to be exploited in the data analysis of the *Euclid* galaxy survey and promise to improve its constraining power, especially on $h$, $\Omega_b$, the neutrino mass, and the DE evolution.

**Key words.** gravitational lensing: weak – cosmological parameters – large-scale structure of Universe

## 1. Introduction

The late-time cosmic acceleration of the Universe on large scales is an established observational fact (Riess et al. 1998; Perlmutter et al. 1998). The cosmological standard model, the so-called ΛCDM model, explains the acceleration as the effect of a cosmological constant Λ, postulating the existence of a dark energy (DE) component with negative pressure, whose physical nature remains poorly understood. Minimal extensions to such a model include a variation of DE with time, parametrized by an equation of state following Chevallier-Polarski-Linder

(CPL): $w(z) = w_0 + z/(1+z)w_a$ (Chevallier & Polarski 2001; Linder 2002). Such pressure impacts not only the background evolution of the Universe, but also the growth of large-scale structure (LSS): the greater the acceleration, the slower the structures can grow. While different DE models have been proposed (Ratra & Peebles 1988; Steinhardt et al. 1999; Armendariz-Picon et al. 2000; Copeland et al. 2006), the physical nature of DE remains unknown. Precise measurements of $w(z)$ promise to shed light on its features.

Massive relic neutrinos can also affect the growth of structure and the expansion history of the Universe. Due to their free streaming with hot thermal velocities

<footnote>* This paper is published on behalf of the Euclid Consortium.</footnote>
(Lesgourgues & Pastor 2006, 2012, 2014; Kiakotou et al. 2008), they alter the epoch of matter-radiation equality and suppress the growth of structure at mildly nonlinear and small scales. The Universe can thus be used as a laboratory to constrain the neutrino mass scale.

Upcoming galaxy surveys are designed to probe the DE equation of state $w$ and the neutrino mass scale, with direct observations and measurements of the growth of structure and of distance scales. Cosmological probes used in these investigations include galaxy clustering (GC) and weak lensing (WL).

At the GC level, measurements of the galaxy correlation function from LSS surveys are used to constrain cosmological parameters, exploiting the sensitivity of the galaxy density fluctuations to the underlying dark matter (DM) density field. At the WL level, images of large ensembles of galaxies provide the so-called cosmic shear: tiny distortions in the shapes of galaxies due to the gravitational potential produced by intervening density perturbations, crossed by light propagating from the source to the observer. Images of galaxy shapes, complemented by an estimation of their redshifts, allow us to measure the structure growth and improve the inference on cosmological parameters.

Recently, aside from traditional GC, under-dense regions known as cosmic voids have been used to extract cosmological information (Lavaux & Wandelt 2010; Sutter et al. 2014; Hamaus et al. 2016, 2017, 2020; Nadathur et al. 2019; Pisani et al. 2019). Their sizes range from about ten to a few hundreds of Mpc (Sheth & Van De Weygaert 2004). Galaxy redshift surveys allow us to build cosmic void catalogs (Sutter et al. 2012; Micheletti et al. 2014; Clampitt & Jain 2015; Sánchez et al. 2017; Mao et al. 2017; Hawken et al. 2020; Pollina et al. 2019; Hamaus et al. 2020; Aubert et al. 2020; Nadathur et al. 2020). As for galaxies and galaxy clusters, it is possible to study the correlation function of cosmic voids and their possible cross-correlation with other cosmological probes (Granett et al. 2008, 2015; Ilić et al. 2013; Hamaus et al. 2014a). In particular, it is possible to measure the cross-correlation between WL and under-dense regions, and exploit this signal to infer cosmological information, and break possible parameter degeneracy (Krause et al. 2013; Melchior et al. 2014; Clampitt & Jain 2015; Sánchez et al. 2017; Fang et al. 2019).

The use of cosmic voids as a cosmological probe presents several advantages. Firstly, voids are less affected than DM halos by shell-crossing and virialization effects, which makes their dynamical evolution more amenable to theoretical models and easier to describe (Sheth & Van De Weygaert 2004; Hamaus et al. 2014b; Stopyra et al. 2020).

Secondly, recent works in the literature (e.g., Lavaux & Wandelt 2012; Pisani et al. 2015, 2019; Massara et al. 2015; Kreisch et al. 2019, 2021; Schuster et al. 2019; Verza et al. 2019; Bos et al. 2012; Lee & Park 2009) have shown that void formation and evolution are sensitive to the DE equation of state, and to the total neutrino mass, in a qualitatively and quantitatively different way with respect to over-dense structures, because of the different scales involved and their different nature.

For instance, the DE contribution within voids can be more dominant than in the Universe, on average, allowing such a probe to have a large sensitivity to its amount and evolution. In this work we consider DE as an effective quintessence “fluid” that does not cluster significantly and is nearly homogeneously distributed in space. Due to their hot thermal velocities, massive neutrinos have a free-streaming length that, depending on their mass and redshift, can range between one hundred to a few tens of $h^{-1}$ Mpc (Lesgourgues & Pastor 2006), thus matching the typical void sizes (Kreisch et al. 2019).

<table><caption>**Table 1.** Cosmological parameters varied in the forecast analysis, together with their values in the fiducial cosmology assumed in this work.</caption><thead><tr><th>Ω<sub>b</sub></th><th>Ω<sub>m</sub></th><th>w<sub>0</sub></th><th>w<sub>a</sub></th><th>h</th><th>M<sub>ν</sub> [eV]</th><th>n<sub>s</sub></th><th>σ<sub>8</sub></th></tr></thead><tbody><tr><td>0.05</td><td>0.32</td><td>-1</td><td>0</td><td>0.67</td><td>0.06</td><td>0.96</td><td>0.816</td></tr></tbody></table>

This paper presents the first forecasts on cosmological parameters based on the combination of WL, void-lensing cross-correlation, and angular void clustering, where voids are found in the galaxy photometric catalog of the Euclid survey, an ESA medium-class mission, currently scheduled for launch in 2023. The paper belongs to a series of companion papers investigating the scientific return that can be expected from voids from the Euclid mission (Hamaus et al. 2021; Contarini et al. 2022).

It is expected that the Euclid photometric galaxy catalog will be characterized by a surface density of $n_g = 30 \text{ arcmin}^{-2}$ (Laureijs et al. 2011), and the imaging Euclid catalog will contain the shapes of about 1.5 billion galaxies, observed in the visible range with the VIS instrument (Cropper et al. 2016). The redshifts of such galaxies will be measured in photometric mode, using the Near Infrared Spectroscopic Photometric (NISP) instrument (Costille et al. 2019), complemented by ground-based observations in different bands.

The parameter forecasts presented here closely follow the recipe of the Euclid Inter Science Task-Force for Forecasts (IST:F) group (Euclid Collaboration 2020, hereafter EC20) in the case of WL and photometric galaxy clustering (GC<sub>ph</sub>): the covariance between cosmological parameters is evaluated with the Fisher matrix approach, using angular power spectra computed within the Limber approximation (Limber 1953), and exploiting a tomography technique. Forecasts are given for spatially flat ΛCDM, vΛCDM, $w_0w_a$CDM, and $vw_0w_a$CDM model cosmologies, adopting the following set of cosmological parameters: the reduced Hubble constant $h$ defined via $H_0 = h \text{ 100 km s}^{-1} \text{ Mpc}^{-1}$, the baryon density parameter $\Omega_b$ at present time, the total matter density parameter $\Omega_m$ at present time, the sum of the three active neutrino masses $M_ν = \sum_i m_{ν_i}$, the massive neutrino density parameter $\Omega_ν = M_ν (\text{eV})/(93.14 h^2)$ at present time, the cold DM density parameter $\Omega_c = \Omega_m - \Omega_b - \Omega_ν$ at present time, where $\Omega_m = \Omega_b + \Omega_c + \Omega_ν$, the DE density parameter $\Omega_{DE} = 1 - \Omega_m$ at present time, the parameters of the DE equation of state, $w_0$ and $w_a$, the rms of the matter linear density fluctuations, $\sigma_8$, inside a radius of $8 h^{-1}$ Mpc, and the scalar spectral index $n_s$. The cosmological parameters varied in the analysis, together with their fiducial values, are summarized in Table 1.

The article is organized as follows. Section 2 illustrates the theoretical modeling of void clustering and bias. Section 3 describes the angular power spectra exploited as observables for parameter forecasts. Section 4 details the evaluation of the Fisher matrix. Section 5 reports the results and interpretation of our analysis. Finally, Sect. 6 presents our concluding remarks.

## 2. The clustering of cosmic voids

We evaluated the angular power spectrum, $C_\ell^{vv}$, of cosmic voids from the void auto-power spectrum, $P_{vv}(k, z)$, while the angular cross-power spectra, $C_\ell^{\gamma\gamma}$, between voids and WL were evaluated from the void-matter cross-spectrum $P_{vm}(k, z)$. The latter are shown, together with the nonlinear matter power spectrum $P_{mm}(k, z)$, in the top left panel of Fig. 1. $P_{vv}(k, z)$ was obtained from the matter power spectrum, $P_{mm}(k, z)$, assuming
the Poisson shot noise (Hamaus et al. 2014a; Chan et al. 2014, 2019; Jamieson & Loverde 2019) and the “effective void bias” $b_v^{\text{eff}}(z)$ defined in Eq. (19). To this aim we adopted the following relations:

$$ P_{vv}(k, z) = [1 - S_N(k)] b_v^{\text{eff}}(z)^2 P_{\text{mm}}(k, z) + 1/\bar{n}^v(z) \\ \equiv b_v^{\text{eff}}(z)^2 \hat{P}_{\text{mm}}(k, z) + 1/\bar{n}^v(z), \quad (1) $$

where the void bias, $b_v^{\text{eff}}(z)$, is assumed to be scale-independent<fnref n="1" />, $\bar{n}^v(z)$ is the void number density, and $\hat{P}_{\text{mm}}(k, z)$ is the nonlinear matter power spectrum when an additional wavenumber filter, $S_N(k)$, is applied to cut small scales. The low-$k$ pass filter, $S_N(k)$, is necessary since in this work we focus on void clustering at large scales, as, in this case, the behavior of the void power spectrum, $P_{vv}(k, z)$, is well described by linear theory via a simple multiplicative factor for the void bias. This can be understood in the context of the halo model formalism (Cooray & Sheth 2002), which has an analog formulation for voids (Voivodic et al. 2020). Such an assumption is inaccurate at small scales (excluded from the analyses in the present work), which are sensitive to the void density profile (Hamaus et al. 2014a). While the inclusion of small scales could improve the results presented in this work, this poses a serious modeling challenge, as to date in the literature there are no accurate models of the cosmological dependence of the void density profile. Furthermore, we point out that the high shot noise would likely reduce the contribution of the small scales. Therefore, we decided for a more conservative approach, considering only scales independent of the particular void profile, which we can confidently include in our analysis.

Analogously, the void-matter cross-spectrum $P_{\text{vm}}(k, z)$ is defined as

$$ P_{\text{vm}}(k, z) = [1 - S_N(k)] b_v^{\text{eff}}(z) P_{\text{mm}}(k, z) \equiv b_v^{\text{eff}}(z) \hat{P}_{\text{mm}}(k, z). \quad (2) $$

To avoid numerical instabilities, we prefer the low-$k$ pass filter rather than a sharp cut in $k$. $S_N(k)$ is the so-called smoothstep function, defined as

$$ S_N(k) = \begin{cases} 0 & \text{if } \frac{k}{k_{\text{ref}}} \le a \\ \left(\frac{k}{k_{\text{ref}}}\right)^{N+1} \sum_{n=0}^{N} \left(\frac{2N+1}{N-n}\right) \left(\frac{N+n}{n}\right) \left(-\frac{k}{k_{\text{ref}}}\right)^n & \text{if } a \le \frac{k}{k_{\text{ref}}} \le b \\ 1 & \text{if } \frac{k}{k_{\text{ref}}} \ge b \end{cases} \quad (3) $$

where $N$ measures the degree of smoothness of the function itself: the first discontinuous derivative of $S_N(k)$ is the $(N+1)\text{th}$ derivative. The order used in $S_N$ for producing $P_{vv}(k, z)$ and $P_{\text{vm}}(k, z)$ is $N=3$. The smoothstep parameters are set to: $a=1$, $b=1.8$, and $k_{\text{ref}} = 0.25 \, \text{h Mpc}^{-1}$; the filter suppresses $P_{vv}(k, z)$ and $P_{\text{vm}}(k, z)$ for $k > k_{\text{ref}}$. The cutoff scale $k_{\text{ref}}$ was determined combining the mean radius of voids $\bar{r}_v$ in the catalog with the void exclusion principle, from which we know that the cutoff scale is given by $k \approx 2\pi/\bar{r}_v$ (Hamaus et al. 2014a). We chose this particular filter because it allowed us to better control its effect, since this function has some more parameters. Furthermore, we checked that the exact shape of the filter did not alter sensibly our results, as it was used to remove numerical artifacts.

Summarizing, we define the probe-dependent power spectrum $\hat{P}_{\text{AB}}$ as

$$ \hat{P}_{\text{AB}}(k, z) = \begin{cases} P_{\text{mm}}(k, z) & \text{if } A = B = \gamma \\ [1 - S_N(k)] P_{\text{mm}}(k, z) & \text{else,} \end{cases} \quad (4) $$

where $\gamma$ stands for cosmic shear.

As we describe below, to obtain the effective void bias, $b_v^{\text{eff}}(z)$, we used the peak-background split (PBS) formalism (Sheth & Van De Weygaert 2004), and weighted the void bias, $b_v(z)$, over the void size function (Sheth & Van De Weygaert 2004; Jennings et al. 2013), that is to say, the comoving number density of voids per radius interval. The void size function was predicted from the “excursion set formalism” in the form of the Sheth & van de Weygart (SvdW) model (Sheth & Van De Weygaert 2004), according to which voids are treated as isolated objects, whose evolution is, therefore, not affected by the environment and is assumed to be spherically symmetric.

Under such assumptions, the prediction of the void abundance is completely described by two physical parameters. The former is the density threshold for structure collapse, $\delta_c$, linearly extrapolated to the present time as $\delta_c(z) = \delta_c(0)/G(z)$, where $G(z)$ is the so-called linear growth factor normalized at $z=0$, and $\delta_c(0)$ slightly depends on cosmology, here assumed to be $\delta_c(0) = 1.686$. The latter is the density threshold for the void formation $\delta_v$, linearly extrapolated to the present time as $\delta_v(z) = \delta_v(0)/G(z)$, with $\delta_v(0) = -0.9$ (see, e.g., Chan et al. 2014; Ronconi & Marulli 2017; Contarini et al. 2019; Verza et al. 2019, recently showing that various values for the threshold $\delta_v(0)$ can be used to model voids more reliably). We chose a value of $\delta_v(0)$ by comparing our effective bias prediction against the measurement from the *Euclid* Flagship simulation (more on this in Sect. 3).

In the void abundance characterization, both $\delta_v$ and $\delta_c$ have an important role. The void formation threshold, $\delta_v$, is the value an under-dense region needs to overcome to turn into a void. The density threshold, $\delta_c$, is necessary to account for the “void-in-cloud problem”, where an under-dense region inside an over-dense one will not become a void, since the latter will collapse.

The SvdW model predicts the void abundance to be a function of both $\delta_c$ and $\delta_v$, via the so-called void-in-cloud parameter, $\mathcal{D}$, and the dimensionless parameter, $x$, defined as

$$ \mathcal{D} = \frac{|\delta_v|}{\delta_c + |\delta_v|}, \quad x = \frac{\mathcal{D}}{|\delta_v|} \sigma, \quad (5) $$

where $\sigma^2$ is the variance of the filtered linear density field on a scale $R$

$$ \sigma^2(R) \equiv S(R) = \int \frac{dk}{k} \frac{k^3 P_{\text{lin}}(k)}{2\pi^2} |W(kR)|^2. \quad (6) $$

Here $P_{\text{lin}}(k)$ is the linear matter power spectrum at $z=0$, and $W(kR)$ is the top-hat filter in Fourier space:

$$ W(kR) = \begin{cases} 1 & kR \le 1 \\ 0 & kR > 1. \end{cases} \quad (7) $$

The abundance of voids with mass $M$ is predicted by the SvdW model as

$$ \frac{dn}{d \ln M} = \frac{\rho_m}{M} f_{\ln \sigma}(\sigma) \frac{d \ln \sigma^{-1}}{d \ln M}, \quad (8) $$

where $\rho_m$ is the mean background density of matter and the multiplicity function, $f_{\ln \sigma}$, represents the fraction of voids in a unit

<footnote n="1">¹ In the presence of massive neutrinos, the void bias could become scale dependent already at the linear level (Schuster et al. 2019), as is also the case for DM halos and galaxies. However, given the small values of the total neutrino mass considered in this work, we assumed the void bias to be scale independent, following the same approach adopted for galaxies in EC20.</footnote>
range of ln M, and it is shown to be (Sheth & Van De Weygaert 2004)

$$f_{\ln \sigma}(\sigma) = 2 \sum_{j=1}^{\infty} \exp\left(-\frac{(j\pi x)^2}{2}\right) j\pi x^2 \sin(j\pi D). \quad (9)$$

Dealing with observations, it is more convenient to work in terms of void volumes rather than void masses. In linear theory, the relation between the void volume and its mass is

$$\frac{\rho_m}{M(r_L)} = \frac{1}{V(r_L)}, \quad (10)$$

where $V(r_L) \equiv 4\pi r_L^3/3$ is the volume of a sphere with radius $r_L$ computed in the linear regime.

Even if the SvdW approach accurately predicts the void mass abundance, the predicted volume fraction of the Universe inside voids is larger than one. In order to fix this issue, the SvdW approach has been extended to the V dn size function from Jennings et al. (2013), who noticed that such a volume fraction issue originates from the incorrect matching between linear and corresponding nonlinear quantities. Since this approach correctly predicts the volume fraction of the Universe inside voids, it is known as the volume-conserving void size function.

Jennings et al. (2013) corrected this match by requiring the void volume fraction of the Universe in the nonlinear regime to have the same value as in the linear one:

$$V(r_v)dn = V(r_L) dn_L|_{r_L=r_L(r_v)}, \quad (11)$$

where the subscript L denotes quantities computed in the linear regime. The V dn void size function then reads (Jennings et al. 2013)

$$\frac{dn}{d \ln r_v} = \frac{f_{\ln \sigma}(\sigma)}{V(r_v)} \frac{d \ln \sigma^{-1}}{d \ln r_v} \frac{d \ln r_L}{d \ln r_v} \Big|_{r_L=r_L(r_v)}, \quad (12)$$

where $r_v$ indicates the actual void radius. The conversion from $r_L$ to $r_v$ is given by

$$\frac{r_v}{r_L} = \left( \frac{\rho_m}{\rho_v} \right)^{1/3}, \quad (13)$$

where $\rho_v$ is the average density within the void.

However, the expression of the multiplicity function, $f_{\ln \sigma}$, in Eq. (9) is difficult to evaluate in practice. Therefore, Jennings et al. (2013) suggested the following approximation that works at 0.2% accuracy:

$$f_{\ln \sigma}(\sigma) \approx \begin{cases} \sqrt{\frac{2}{\pi} \frac{|\delta_v|}{\sigma}} \exp\left(-\frac{\delta_v^2}{2\sigma^2}\right) & x \le 0.276 \\ 2 \sum_{j=1}^4 \exp\left(-\frac{(j\pi x)^2}{2}\right) j\pi x^2 \sin(j\pi D) & x > 0.276. \end{cases} \quad (14)$$

Finally, the V dn size function is given by

$$\frac{dn}{d \ln r_v} = \frac{f_{\ln \sigma}(\sigma)}{V(r_v)} \frac{d \ln \sigma^{-1}}{d \ln r_v} \Big|_{r_L=r_v/1.7}, \quad (15)$$

and, after some rearrangements, its final expression is

$$\frac{dn}{d r_v} = \frac{f_{\ln \sigma}(\sigma)}{r_v V(r_v)} \left[ -\frac{r_L}{\sigma} \frac{d\sigma}{d r_L} \right]_{r_L=r_v/1.7}, \quad (16)$$

where $r_v = 1.7 r_L$ is predicted by considering a spherical expansion.

As mentioned previously, the void bias, $b_v$, at a fixed void radius, $r_v$, is computed using the PBS formalism (Sheth & Van De Weygaert 2004; Chan et al. 2014):

$$b_v(z, r_v) = 1 + \frac{v^2 - 1}{\delta_v} + \frac{\delta_v D}{4\delta_c^2 v^2}, \quad (17)$$

where the z dependence is encapsulated in $\delta_c$ and $\delta_v$, and the $r_v$ dependence in $v$:

$$v \equiv \frac{|\delta_v|}{\sigma(r_v)}. \quad (18)$$

We defined the effective void bias by weighting Eq. (17) with the void size function V dn in Eq. (16):

$$b_v^{\text{eff}}(z) = \frac{\int_{r_{v,\min}}^{r_{v,\max}} \frac{dr_v}{dr_v} b_v(z, r_v)}{\int_{r_{v,\min}}^{r_{v,\max}} \frac{dr_v}{dr_v} \frac{dn}{dr_v}}, \quad (19)$$

exploiting in this way all void radii $r \ge r_{v,\min} h^{-1}$ Mpc. In this work, we assumed $r_{v,\min} = 25 h^{-1}$ Mpc, and we provide forecasts treating the evolution of the void bias in two main scenarios. A “pessimistic void bias” scenario, where the redshift evolution of the void bias is supposed to be known from linear theory, but its absolute normalization is unknown. In this case, $b_v^{\text{eff}}$ is evolved with the growth factor of the fiducial cosmology. Its fiducial value today is assumed to be $b_v^{\text{eff}}(z = 0) = -11.9$, according to Eq. (19), and it is marginalized over as a nuisance parameter. This void bias scenario will be combined with the pessimistic WL scenario as described in Sect. 3. An “optimistic void bias” scenario, where we exploit the cosmology dependence of both the void size function and void bias. In particular, such dependence is encapsulated in the growth factor evolution of the density thresholds $\delta_c$ and $\delta_v$, which appear both in the void bias definition and in the void size function. This configuration is considered to be optimistic, since we assume that the value of $b_v$, entering $b_v^{\text{eff}}$, is known from linear theory, and therefore its dependence on the cosmological parameters is exploitable. This void bias scenario will be combined with the optimistic WL scenario as described in Sect. 3. To summarize, we present forecasts in two scenarios: a pessimistic one, where we treat $b_v^{\text{eff}}(z = 0)$ as an unknown quantity to be marginalized over, fixing its evolution to the fiducial cosmology; and an optimistic one, where we assume we are able to model the cosmological dependence of $b_v^{\text{eff}}(z)$ and use it as a further source of information on cosmological parameters. Providing forecasts in both cases is, in our opinion, the fairest approach, since we present parameter constraints spanning from the worst to the best case scenario.

### 3. Observables: Angular power spectra

This section summarizes the observables used for cosmological parameter forecasts based on WL, voids, and their cross-correlation. Here we describe in particular the angular power spectra $C(\ell)$, which are functions of the multipole number $\ell$, and are defined as the spherical harmonic transform of a two-point correlation function. The $a_{\ell m}$, the coefficients of the spherical harmonics decomposition of a field $\mathcal{F}$ on the sphere (e.g., the cosmic shear or void density fields), are defined on the sky as

$$a_{\ell m} = \int d\Omega \mathcal{F}(\theta, \phi) Y_{\ell m}(\theta, \phi). \quad (20)$$

The $C(\ell)$ represent the two-point correlation function of the $a_{\ell m}$ and they are diagonal in $\ell$ and independent of $m$ because of statistical isotropy.
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Figure containing the power spectra used in this work. **Top left:** nonlinear matter auto-power spectrum $P_{\text{mm}}$ (solid red line), the void auto-power spectrum $P_{\text{vv}}$ without the shot noise (dashed blue line), computed assuming the pessimistic void bias scenario with $b_{\text{eff}} \approx -11.9$, and the absolute value of the void-matter cross-power spectrum $P_{\text{vm}}$ (dot-dashed black line), all computed at redshift $z = 0.001$. **Top right:** tomographic void angular auto-power spectra $C_{ij}^{\text{vv}}(\ell)$ for the diagonal ($i = j$) tomographic bins, without shot noise. **Bottom left:** tomographic lensing auto-power spectra $C_{ij}^{\gamma}(\ell)$ for the diagonal ($i = j$) tomographic bins, without shape noise. **Bottom right:** tomographic void-lensing angular cross-spectra $C_{ij}^{\gamma\gamma}(\ell)$ for the diagonal ($i = j$) tomographic bins. All the spectra are theoretically evaluated in the reference cosmology reported in Table 1 and the bin endpoints are given in Eq. (24).</figcaption></figure>

In the forecasts presented in this work, we used three kinds of $C(\ell)$: the void-void auto-correlation $C^{\text{vv}}(\ell)$, the lensing-lensing auto-correlation $C^{\gamma\gamma}(\ell)$, and the void-lensing cross-correlation $C^{\gamma}(\ell)$. We computed the $C(\ell)$ “tomographically” (Hu 1999) in a set of redshift bins, that is to say they were evaluated over two redshift bins $i$ and $j$, and were generically denoted as $C_{ij}^{\text{AB}}(\ell)$. For example, $C_{ij}^{\gamma\gamma}(\ell)$ is the spherical harmonic transform of the correlation function between the lensing signal in the $i$th bin and the void signal in the $j$th bin.

We evaluated $C^{\text{vv}}(\ell)$ and $C^{\gamma\gamma}(\ell)$ for $\ell \in [10, 1500]$, while for $C^{\gamma\gamma}(\ell)$, two scenarios were considered: a “pessimistic” WL scenario, with $\ell \in [10, 1500]$, and an “optimistic” WL scenario, with $\ell \in [10, 5000]$.

The Limber approximation (Limber 1953; LoVerde & Afshordi 2008) allowed us to write a simple integral expression for the tomographic $C(\ell)$, where the power spectrum $\hat{P}_{\text{AB}}(k, z)$, defined in Eq. (4), enters the integral as

$$C_{ij}^{\text{AB}}(\ell) \approx \frac{c}{H_0} \int_{z_{\text{min}}}^{z_{\text{max}}} dz \frac{W_i^A(z) W_j^B(z)}{E(z) r^2(z)} \hat{P}_{\text{AB}}\left[\frac{\ell + 1/2}{r(z)}, z\right], \quad (21)$$

where $W_i^{A/B}(z)$ are suitable “weight” functions for the different probes A and B defined in Eqs. (27) and (30), $E(z)$ is the

dimensionless Hubble parameter,

$$\begin{aligned} E(z) &= \frac{H(z)}{H_0} \\ &= \sqrt{\Omega_m(1+z)^3 + (1-\Omega_m)(1+z)^{3(1+w_0+w_a)}e^{-3w_a z/(1+z)}}, \end{aligned} \quad (22)$$

and $r(z)$ is the comoving distance,

$$r(z) = \frac{c}{H_0} \int_{0}^{z} \frac{dz'}{E(z')}. \quad (23)$$

Throughout the forecast, the integral expression of $C_{ij}^{\text{AB}}(\ell)$ was computed numerically with the Simpson method (Brew & Atkinson 1979). The signals $C_{ij}^{\text{vv}}(\ell)$, $C_{ij}^{\gamma\gamma}(\ell)$, and $C_{ij}^{\gamma\gamma}(\ell)$, evaluated using the reference cosmology in Table 1, are shown in Fig. 1.

In the *Euclid* photometric survey, the minimum and the maximum redshifts are $z_{\text{min}} = 0.001$ and $z_{\text{max}} = 2.5$, respectively. The $C(\ell)$ and the distributions of galaxies and voids are sampled in ten redshift bins $i$, with the following endpoints:

$$\{z_e\} = \{0.001, 0.418, 0.560, 0.678, 0.789, 0.900, 1.019, 1.155, 1.324, 1.576, 2.50\}. \quad (24)$$

These bins were chosen to be equi-populated in the photometric galaxy catalog. Although, as explained in Sect. 3, the mock
void catalog was split in ten “equi-spaced” redshift bins, we preferred to convert the redshift void distribution to the same binning adopted for galaxies.

In order to simplify the notation of the integrand in Eq. (21), we defined the “kernel” functions as

$$K_{ij}^{AB}(z) = \frac{c}{H_0} \frac{W_i^A(z)W_j^B(z)}{E(z)r^2(z)}, \qquad (25)$$

with *i*, *j* being the tomographic indices, and A, B the considered probes.

The tomographic void-void angular power spectra $C_{ij}^{vv}(l)$ were computed as

$$C_{ij}^{vv}(l) = \int_{z_{\min}}^{z_{\max}} dz K_{ij}^{vv}(z) \hat{P}_{\text{mm}} \left[ \frac{l + 1/2}{r(z)}, z \right], \qquad (26)$$

where $\hat{P}_{\text{mm}}(k, z)$ is defined in Eq. (1) and we omitted the shot-noise additive term $1/\bar{n}_i^v$ defined in Eq. (40).

The void weight function is defined as

$$W_i^v(z) = \frac{H(z)}{c} n_i^v(z) b_{v,i}^{\text{eff}}. \qquad (27)$$

The void effective bias, $b_{\text{veff},i}$, entering in the void weight function, is given by $b_{\text{veff}}(z)$, as defined in Eq. (19), evaluated in the center of the ith tomographic bin. The un-normalized projected void density distribution in redshift was measured directly from the *Euclid* Flagship mock photometric galaxy catalog<fnref n="2" />, as explained in Sect. 3.

The tomographic void-lensing angular power spectra, $C_{ij}^{\gamma\gamma}(l)$, and the lensing-lensing angular power spectra, $C_{ij}^{\gamma\gamma}(l)$, were computed, respectively, as

$$C_{ij}^{\gamma\gamma}(l) = \int_{z_{\min}}^{z_{\max}} dz K_{ij}^{\gamma\gamma}(z) \hat{P}_{\text{mm}} \left[ \frac{l + 1/2}{r(z)}, z \right] \qquad (28)$$

and

$$C_{ij}^{\gamma\gamma}(l) = \int_{z_{\min}}^{z_{\max}} dz K_{ij}^{\gamma\gamma}(z) P_{\text{mm}} \left[ \frac{l + 1/2}{r(z)}, z \right], \qquad (29)$$

where in $C_{ij}^{\gamma\gamma}(l)$ we omitted the shape-noise term given by Eq. (40). Following EC20, in the kernels of Eqs. (28) and (29), the lensing weight function reads

$$W_i^\gamma(z) = \frac{3}{2} \left(\frac{H_0}{c}\right)^2 \Omega_m (1+z)r(z) \tilde{W}_i^\gamma(z), \qquad (30)$$

with

$$\tilde{W}_i^\gamma(z) = \int_z^{z_{\max}} dz' n_i^g(z') \left[ 1 - \frac{r(z)}{r(z')} \right]. \qquad (31)$$

Here $\tilde{W}_i^\gamma(z)$ is the lensing efficiency, defined as an integral expression of the observed galaxy distribution $n_i^g(z)$.

The density distribution of the “observed” galaxies in the ith tomographic bin, $n_i^g(z)$, was computed as a convolution of the galaxy distribution, $n^g(z)$, and the photometric instrument response, as described below. According to EC20, the galaxy redshift distribution, $n^g(z)$, adopted in this work is

$$n^g(z) \propto \left(\frac{z}{z_0}\right)^2 \exp\left[-\left(\frac{z}{z_0}\right)^{3/2}\right], \qquad (32)$$

<table><caption>**Table 2.** Parameters used in the photometric redshift distribution $p_{\text{ph}}(z_{\text{p}}|z)$ of Eq. (33).</caption><thead><tr><th>c<sub>b</sub></th><th>z<sub>b</sub></th><th>σ<sub>b</sub></th><th>c<sub>o</sub></th><th>z<sub>o</sub></th><th>σ<sub>o</sub></th><th>f<sub>out</sub></th></tr></thead><tbody><tr><td>1.0</td><td>0.0</td><td>0.05</td><td>1.0</td><td>0.1</td><td>0.05</td><td>0.1</td></tr></tbody></table>

where $z_0 = 0.9/\sqrt{2}$ (Laureijs et al. 2011). Redshift measurements are affected by experimental and reconstruction effects. We modeled the probability $p_{\text{ph}}(z_{\text{p}}|z)$ that a true galaxy redshift $z$ will be reconstructed with a photometric redshift $z_{\text{p}}$ as

$$p_{\text{ph}}(z_{\text{p}}|z) = \frac{1 - f_{\text{out}}}{\sqrt{2\pi}\sigma_{\text{b}}(1+z)} \exp\left\{-\frac{1}{2}\left[\frac{z - c_{\text{b}}z_{\text{p}} - z_{\text{b}}}{\sigma_{\text{b}}(1+z)}\right]^2\right\} \\ + \frac{f_{\text{out}}}{\sqrt{2\pi}\sigma_{\text{o}}(1+z)} \exp\left\{-\frac{1}{2}\left[\frac{z - c_{\text{o}}z_{\text{p}} - z_{\text{o}}}{\sigma_{\text{o}}(1+z)}\right]^2\right\}. \qquad (33)$$

The modeling and the parameters entering $p_{\text{ph}}(z_{\text{p}}|z)$ are reported in Table 2, and were taken from EC20. The convolution of $n^g(z)$ and $p_{\text{ph}}(z_{\text{p}}|z)$ estimates the galaxy redshift distribution as measured by NISP in the photometric mode, together with ground-based observation. The density distribution, $n_i^g(z)$, of the “observed” galaxies in the ith tomographic bin is obtained convolving $p_{\text{ph}}(z_{\text{p}}|z)$ with $n^g(z)$ in the ith bin

$$n_i^g(z) = \frac{\int_{z_i^-}^{z_i^+} dz_p n^g(z) p_{\text{ph}}(z_p|z)}{\int_{z_{\min}}^{z_{\max}} dz \int_{z_i^-}^{z_i^+} dz_p n^g(z) p_{\text{ph}}(z_p|z)}, \qquad (34)$$

where $z_i^+$ and $z_i^-$ are the edges of the ith bin as defined in Eq. (24). The density distributions, $n_i^g(z)$, in the different tomographic bins are shown in the left panel of Fig. 2.

The matter power spectrum, $P_{\text{mm}}(k, z)$, was computed numerically using as Boltzmann solver the “Code for Anisotropies in the Microwave Background”<fnref n="3" /> (CAMB; Lewis et al. 1999)<fnref n="4" />. In order to obtain the cosmological parameter forecasts presented in this work, we computed the matter power spectra, $P_{\text{mm}}(k, z)$, for each set of input cosmological parameters, which were varied with respect to the fiducial cosmology in Table 1 to also compute the derivatives entering the Fisher matrix approach (Fisher 1935). CAMB evaluates $P_{\text{mm}}(k, z)$ over a $k-z$ grid. According to EC20, the range of redshifts chosen in the present analysis is $z \in [0.001, 2.5]$, the range of wave numbers is $k \in [10^{-5}, 400]$ Mpc$^{-1}$, and at least 300 equispaced points in redshift are used. CAMB performs a logarithmic $k$-binning. In the forecast, at least 50 steps per decade are used. The impact of different $k$ and $z$ binnings on the final results is reported in Sect. 5.6.

Since both $n_i^g(z)$ and $n_i^v(z)$ are normalized quantities, the variation of the cosmological parameters in the Fisher forecast computation does not (or negligibly) affect these quantities. In fact, the main effect on number counts from changing the model cosmology is due to the associated volume change, which cancels out for normalized quantities. Instead, the shot-noise matrix in Eq. (40) is un-normalized and computed at the fiducial cosmology assumed in the Fisher forecast computation. The cosmological parameter variation is needed for the computations of the

<footnote n="2">² Euclid Collaboration (in prep.).</footnote>

<footnote n="3">³ https://camb.info/</footnote>

<footnote n="4">⁴ The CAMB version used in this work is v1.0.4. The different CAMB options were set to include neutrino effects and also nonlinear evolution via the Takahashi halofit recipe (Takahashi et al. 2012; Bird et al. 2011).</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Galaxy and void distributions used in this work. Left: dimensionless normalized projected galaxy distributions $n_i^g(z)$ over the tomographic bins i. Right: dimensionless normalized projected void distribution in equi-spaced redshift bins. The corresponding un-normalized projected galaxy and void distributions are reported in Table 3 in units of arcmin$^{-2}$.</figcaption></figure>

$C(\ell)$ derivatives, and enters the matter power spectra, $P_{\text{mm}}(k, z)$, the growth factor, $G(z)$, the Hubble parameter, $H(z)$, and the comoving distance, $r(z)$, where the latter two terms change the weight functions $W_i^A$.

In order to use realistic estimates of void bias and void distribution from the *Euclid* photometric galaxy sample, in the present analysis we employed the Flagship mock galaxy catalog (Euclid Collaboration, in prep.). This is based on an N-body simulation of 12 600$^3$ DM particles in a periodic box of 3780 h$^{-1}$ Mpc on a side (Potter et al. 2017), with a flat $\Lambda$CDM cosmology very similar to that reported in Table 1, namely $\Omega_m = 0.319$, $\Omega_b = 0.049$, $\Omega_\Lambda = 0.681$, $\sigma_8 = 0.83$, and $n_s = 0.96$, $h = 0.67$, as obtained by Planck (Planck Collaboration XIII 2016). Dark matter halos were identified with the ROCKSTAR halo finder (Behroozi et al. 2013), and populated with central and satellite galaxies using a halo occupation distribution (HOD) framework to reproduce the relevant observables for the *Euclid* main cosmological probes. The HOD algorithm (Carretero et al. 2015; Crocce et al. 2015) was calibrated exploiting several observational constraints, such as the local luminosity function for the faintest galaxies (Blanton et al. 2003, 2005) and GC statistics as a function of luminosity and color (Zehavi et al. 2011). The resulting Flagship galaxy mock lightcone spanned one octant of the sky and simulated both spectroscopic and photometric *Euclid* galaxy samples. In this paper we consider the latter, which extends up to redshift $z = 2.3$ and in which a Gaussian photometric redshift error of $\Delta z = 0.05(1 + z)$ was applied to each galaxy.

To identify cosmic voids in this catalog, we applied the 2D void finder of Sánchez et al. (2017), Vielzeuf et al. (2019) to the photometric galaxy sample, computing the number of voids in ten “equi-spaced” redshift bins, as shown in the right panel of Fig. 2, where the ratio between the void number in each bin and the bin width represents the void projected density at the bin center. We find that the 2D void population traced by photometric galaxies in the Flagship catalog extends from $r_{v,\text{min}} \sim 25$ h$^{-1}$ Mpc up to $r_{v,\text{max}} \sim 300$ h$^{-1}$ Mpc, and its projected spatial density, $n^v(z)$, is obtained by interpolating the void density in each bin center.

Moreover, to verify that the bias modeling, Eq. (19), used in our forecasts was representative of *Euclid* observations, we measured the void bias redshift evolution in the obtained Flagship photometric void catalog. To this aim, we followed the methodology presented in Hamaus et al. (2014a). We used

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Effective void bias (green error bars) and galaxy bias (red error bars) as measured in the Flagship simulation versus the theoretical effective void bias (blue line) of Eq. (19). Void bias measurements have a bigger error bar, compared to galaxy bias measurements, since there are fewer voids than galaxies.</figcaption></figure>

the open-source code nbodykit (Hand et al. 2018) to compute the void auto power-spectrum in eight redshift bins ($z = [0.3, 0.5, 0.75, 1.0, 1.25, 1.5, 1.75, 2.]$) with a bin size of $\Delta z = 0.2$. Then, from Eq. (1), we inferred the void bias as

$$b_v^{\text{eff}}(z) = \sqrt{\frac{P_{vv}(k, z) - 1/\bar{n}^v}{P_{\text{mm}}(k, z)}}, \quad (35)$$

where the void power spectrum was measured in the Flagship photometric void catalog. Then, we averaged the measured bias in the range $0.05 < k(h^{-1} \text{Mpc}) < 0.1$, that it to say, in the regime in which the void bias is constant as a function of the scale (see Hamaus et al. 2014a).

The evolution of the void bias as a function of redshift in the Flagship photometric void catalog (green error bars) is shown in Fig. 3. The measurement errors were estimated using a jackknife resampling. The red stars in Fig. 3 show the galaxy bias measured from the Flagship mock catalog using Eq. (35). Finally the solid blue line shows the theoretical bias modeling, Eq. (19), obtained via the PBS formalism and adopted for our forecast computation. We note that the measured negative void bias shows a redshift evolution close enough to the theoretical one. Therefore, we consider the PBS prescription of the void bias
evolution acceptable for the purposes of the forecast analysis presented in this work. We accounted for any theoretical uncertainty on the void threshold and size function by marginalizing over the normalization of the effective bias in the pessimistic case. We verified that the agreement between the theoretical void size function and the one measured from the photometric void sample in the Flagship simulation was good enough for the forecasting purposes of this work, and we postpone to a future work a detailed calibration against simulated data, when more realistic mocks, accounting for detailed Euclid survey specification and systematics, will be available.

## 4. The void-lensing Fisher matrix

The Fisher matrix formalism (Fisher 1935) is used to predict uncertainties of cosmological parameter measurements from different probes. This section summarizes the approach and provides details on its application to the cases considered in this work.

The Fisher matrix is defined as the expectation value of the second derivatives of the logarithm of the likelihood L:

$$F_{\alpha\beta} = - \left\langle \frac{\partial^2 \ln L}{\partial \alpha \partial \beta} \right\rangle, \qquad (36)$$

where $\alpha$ and $\beta$ are the parameters of interest. The expected “error covariance matrix” is the inverse of the Fisher matrix:

$$C_{\alpha\beta} = (F^{-1})_{\alpha\beta}. \qquad (37)$$

The diagonal elements of the error covariance matrix are the squares of the marginalized 1-$\sigma$ errors on the parameters:

$$\sigma_{\alpha} = \sqrt{C_{\alpha\alpha}}. \qquad (38)$$

As we deal with angular power spectra, it is convenient to introduce the matrix $\Sigma_{ij}^{AB}$ associated with a given $C(\ell)$:

$$\Sigma_{ij}^{AB}(\ell) = \sqrt{\frac{2}{(2\ell + 1)\Delta\ell f_{\text{sky}}}} [C_{ij}^{AB}(\ell) + N_{ij}^{AB}(\ell)], \qquad (39)$$

where $\Delta\ell$ is the multipole bin width, and $f_{\text{sky}}$ the sky fraction covered by the survey. The “shot-noise matrix”, $N_{ij}^{AB}(\ell)$, depends on the particular probe combination. For the $C(\ell)$ of the void clustering, WL, and void-lensing, assuming Poisson statistics, they are respectively

$$N_{ij}^{\gamma\gamma} = \frac{1}{\bar{n}_i^\gamma} \delta_{ij}, \quad N_{ij}^{\gamma\gamma} = \frac{\sigma_\epsilon^2}{\bar{n}_i^\gamma} \delta_{ij}, \quad N_{ij}^{\gamma\gamma} = N_{ij}^{\gamma\gamma} = 0, \qquad (40)$$

where $\delta_{ij}$ is the Kronecker delta, $\sigma_\epsilon^2$ is the galaxy shape noise, and $\bar{n}_i^\gamma$ and $\bar{n}_i^\gamma$ are the un-normalized average galaxy and void surface densities in the $i$th tomographic bin computed in the fiducial cosmology, respectively. The survey specifications used to compute the covariances are shown in Table 3.

In the case of a single probe (A=B), the covariance matrix of the $a_{\ell m}$ is simply given by Eq. (39). When two or more probes are combined together, one needs to construct the full covariance matrix, $\Sigma^{\text{XC}}$, composed of matrix blocks defined in Eq. (39). For the probes $\gamma$ and $\nu$ considered in this work, $\Sigma^{\text{XC}}$ is given by

$$\Sigma^{\text{XC}}(\ell) = \begin{pmatrix} \Sigma^{\gamma\gamma}(\ell) & \Sigma^{\gamma\nu}(\ell) \\ \Sigma^{\nu\gamma}(\ell) & \Sigma^{\nu\nu}(\ell) \end{pmatrix}. \qquad (41)$$

The covariance and the $C(\ell)$ matrix share the same structure, and therefore the block matrix reads

$$\mathbf{C}^{\text{XC}}(\ell) = \begin{pmatrix} \mathbf{C}^{\gamma\gamma}(\ell) & \mathbf{C}^{\gamma\nu}(\ell) \\ \mathbf{C}^{\nu\gamma}(\ell) & \mathbf{C}^{\nu\nu}(\ell) \end{pmatrix}. \qquad (42)$$

Assuming that the $a_{\ell m}$, defined in Eq. (20), follow a multivariate Gaussian distribution, an analytical expression for the Fisher matrix elements is given by<fnref n="5" />

$$F_{\alpha\beta} = \sum_{\ell=\ell_{\min}}^{\ell_{\max}} \text{Tr} \left\{ [\Sigma(\ell)]^{-1} \frac{\partial \mathbf{C}(\ell)}{\partial \alpha} [\Sigma(\ell)]^{-1} \frac{\partial \mathbf{C}(\ell)}{\partial \beta} \right\}. \qquad (43)$$

The formula above applies both to the single- and two-probe correlation case, provided that the $C(\ell)$ and the covariance matrices are chosen accordingly. The Fisher matrix computation involves the derivatives of the $C(\ell)$ with respect to cosmological and nuisance parameters. In order to ensure reliable results, numerical derivatives were computed with<fnref n="6" /> a numerical approach based on the SteM fitting procedure (Camera et al. 2017), which is based on a iterative linear regression, and a semi-analytical approach, in which both analytical and numerical derivatives are used.

## 5. Results: Correlations, errors, and the inclusion of the void bias evolution

In this section we present the results obtained implementing the analysis described in previous sections. We discuss the forecasts of cosmological parameters, expressed as Fisher matrix marginalized contours, obtained combining the expected Euclid WL and angular void clustering.

We report parameter forecasts in two different scenarios, namely the pessimistic and optimistic setups, described in Sect. 2 for the void bias, and in Sect. 3 for the WL, respectively. They are also summarized in Table 4. The 1-$\sigma$ errors, together with the DE figure of merit<fnref n="7" /> (FoM), in both scenarios, and different cosmological models, are reported in Table 5. The corresponding contour plots are shown in Figs. 4–8. Focusing on the DE equation of state and on the total neutrino mass, the contour plots for $w_0$, $w_a$, $M_\nu$, evaluated respectively for the pessimistic and optimistic scenarios, are also shown in the top panels of Fig. 9.

### 5.1. The standard ΛCDM cosmology

The contour plots for the ΛCDM cosmology are reported in Fig. 4, the pessimistic case in the top panel and the optimistic case in the bottom panel. In both the pessimistic and optimistic scenarios, we can observe that WL has a larger constraining power than photometric void clustering, except for $h$ and $\Omega_b$ in the optimistic case, where for $h$ void clustering provides better constraints than WL by a factor of ~3, while for $\Omega_b$ the constraints from the two single probes look comparable. This is due to the form of the WL kernel and the integration along

<footnote n="5">⁵ In Appendix A we show that this expression is valid both in the so-called “field” and “estimator” perspectives.</footnote>

<footnote n="6">⁶ In an early stage of the analysis, numeric derivatives were also performed with a five-point stencil method. However, this method has shown numeric instabilities, and was therefore discarded.</footnote>

<footnote n="7">⁷ The FoM is defined as the inverse of the square root of the determinant of the covariance matrix $C_{w_0 w_a}$ and is inversely proportional to the area of the contour ellipse $w_0 - w_a$ in the marginalized parameter plane (EC20). The covariance matrix is the inverse of the Fisher matrix.</footnote>
<table><caption>**Table 3.** Survey specifications entering Eqs. (39) and (40).</caption><thead><tr><th>f<sub>sky</sub></th><th>σ<sub>ε</sub></th><th>n&#773;<sub>i</sub><sup>g</sup></th><th>n&#773;<sub>1</sub><sup>v</sup> × 10<sup>4</sup></th><th>n&#773;<sub>2</sub><sup>v</sup> × 10<sup>4</sup></th><th>n&#773;<sub>3</sub><sup>v</sup> × 10<sup>4</sup></th><th>n&#773;<sub>4</sub><sup>v</sup> × 10<sup>4</sup></th><th>n&#773;<sub>5</sub><sup>v</sup> × 10<sup>4</sup></th><th>n&#773;<sub>6</sub><sup>v</sup> × 10<sup>4</sup></th><th>n&#773;<sub>7</sub><sup>v</sup> × 10<sup>4</sup></th><th>n&#773;<sub>8</sub><sup>v</sup> × 10<sup>3</sup></th><th>n&#773;<sub>9</sub><sup>v</sup> × 10<sup>3</sup></th><th>n&#773;<sub>10</sub><sup>v</sup> × 10<sup>3</sup></th></tr></thead><tbody><tr><td>0.3636</td><td>0.3</td><td>3</td><td>3.03</td><td>3.52</td><td>3.68</td><td>4.09</td><td>4.7</td><td>5.71</td><td>7.39</td><td>1.05</td><td>1.82</td><td>9.36</td></tr></tbody></table>

Notes. Here $\bar{n}_i^g$ and $\bar{n}_i^v$ are in units of arcmin$^{-2}$.

<table><caption>**Table 4.** Probe configurations in the pessimistic and optimistic scenarios for the combination V + WL, both considered for two values of $r_{v,min}$.</caption><thead><tr><th>Configuration</th><th>V + WL – Pessimistic</th><th>V + WL – Optimistic</th></tr></thead><tbody><tr><td>ℓ-range (WL)</td><td>10–1500</td><td>10–5000</td></tr><tr><td>b<sub>v</sub> modeling</td><td>b<sub>v</sub><sup>eff</sup>(z = 0) nuisance</td><td>b<sub>v</sub><sup>eff</sup>(z) exploited</td></tr><tr><td>r<sub>v,min</sub></td><td colspan="2">25 h<sup>-1</sup> Mpc</td></tr><tr><td>ℓ-range (V)</td><td colspan="2">10–1500</td></tr><tr><td># of evaluated k per decade (CAMB)</td><td colspan="2">60</td></tr><tr><td># of evaluated z (CAMB)</td><td colspan="2">450</td></tr><tr><td>Differentiation method</td><td colspan="2">SteM</td></tr></tbody></table>

the line of sight, and to the presence of baryon acoustic oscillations (BAO) in the void angular auto power spectrum, $C^{vv}(\ell)$, and the void-lensing cross spectrum, $C^{\gamma\gamma}(\ell)$, (see Fig. 1), which are instead completely washed out in the WL angular power spectrum, $C^{\gamma\gamma}(\ell)$. We also remark that, in general, both in the pessimistic and optimistic scenarios, the parameters for which the constraints improve most, when combining void clustering and WL, are not only $h$ and $\Omega_b$, but also $n_s$; similar conclusions (without WL) were found by Kreisch et al. (2021).

In the case of $n_s$, this happens since the WL and void clustering ellipses in the $n_s-\Omega_b$ plane happen to be orthogonal. We verified that, even if $n_s$ and $\Omega_b$ produce an increment in the same direction for $C^{\gamma\gamma}(\ell)$ and $C^{vv}(\ell)$, which would imply that the two parameters are negatively correlated, this orthogonality comes from projecting the large parameter space onto the $n_s-\Omega_b$ 2D space; in other words, it is due to the Fisher matrix inversion.

Finally, when moving from the pessimistic to the optimistic scenario, all marginalized errors decrease: in particular the error on $n_s$ decreases by a factor of ~2, the constraints on $\sigma_8$ and $\Omega_m$ get ~40% tighter, and the uncertainty on $h$ decreases by a factor of ~3. The increase in the constraining power on $\sigma_8$ is somehow expected: the void bias affects the overall amplitude of the $C(\ell)$, so it is degenerate with $\sigma_8$. On the one hand, in the pessimistic scenario, both the value of the bias parameters and $\sigma_8$ were assumed to be measured from data, and this increases the uncertainties on $\sigma_8$. On the other hand, in the optimistic scenario, we assumed the cosmological dependence of the effective void bias to be known and we exploited it in the forecasts of the cosmological parameters.

## 5.2. The massive neutrino cosmology

Here we present the results in the case where the total neutrino mass, $M_ν$, is not kept fixed to its fiducial value, but it is a quantity whose value has to be determined from the data. Strictly speaking, in the case of the neutrino mass, the Fisher matrix approach goes beyond its range of validity, as the associated likelihood is non-Gaussian being truncated at $M_ν = 0$. However, as shown in Fig. 1 of Brinckmann & Lesgourgues (2019), the posterior obtained from the MCMC is in good agreement with the one obtained from the truncated Fisher matrix. Furthermore, we

notice that, in the case of WL, the neutrino constraints presented in this work do not represent the final Euclid-WL constraints, as intrinsic alignments are not taken into account in this analysis. Moreover, a varying $M_ν$ is not considered when the spectroscopic galaxy clustering, $GC_{sp}$, is included in Sect. 5.4. The contour plots for this cosmology are reported in Fig. 5, with the pessimistic case on the left and the optimistic one on the right. On the one hand, adding $M_ν$ as a free parameter mainly impacts the measurement of $\sigma_8$, weakening its constraints by a factor of ~2 with respect to the baseline ΛCDM case. This is expected since the effects of $\sigma_8$ and massive neutrinos on the matter power spectrum are similar: the former regulates the normalization of $P_{mm}$, and the latter suppresses $P_{mm}$ in a scale-dependent way, due to neutrinos free streaming (Lesgourgues & Pastor 2006). On the other hand, such a suppression allows void clustering to impact positively on the neutrino mass measurements: the error on $M_ν$ decreases with respect to the WL case alone, by ~5% and ~15% in the pessimistic and optimistic scenarios, respectively. In addition to considering a fiducial total mass of $M_ν = 0.06$ eV in the normal hierarchy scenario, we also computed parameter forecasts assuming a neutrino mass degenerate scenario with fiducial total mass $M_ν = 0.15$ eV$^8$, from the best-fit value in Pellejero-Ibanez et al. (2017). With this choice, the marginalized errors on $M_ν$ decrease by ~5%. This is expected: the higher the neutrino mass, the more the matter power spectrum is sensitive to its effects (Lesgourgues & Pastor 2006).

## 5.3. The dynamical dark energy cosmology

In this cosmological scenario, we analyzed the $w_0w_a$CDM model, considering a time-dependent DE equation of state. The contour plots for this cosmology are reported in Fig. 6, both in the pessimistic and optimistic scenarios. Again, adding $w_0$ and $w_a$ mainly affects $h$, $\sigma_8$, and $\Omega_m$, increasing their errors. The impact on $h$ and $\Omega_m$ is explained considering that $w_0$ and $w_a$ enter the Hubble parameter $H(z)$ (22), where the dependence from the DE equation of state has been exploited. Moreover,

<footnote>$^8$ The value $M_ν = 0.15$ eV is compatible with either an inverted or a normal hierarchy (Jimenez et al. 2010); we chose a normal hierarchy in this case, as implemented in EC20.</footnote>
<table><caption>**Table 5.** Marginalized 1-σ errors in different cosmological scenarios for the galaxy weak lensing (WL) and angular void clustering (V) probes, together with their combinations in the case they are assumed to be independent, WL + V, and when their cross-correlation is included.</caption><thead><tr><th>Probe</th><th><i>h</i></th><th>Ω<sub>m</sub></th><th>Ω<sub>b</sub></th><th>σ<sub>8</sub></th><th>n<sub>s</sub></th><th><i>M</i><sub>v</sub> [eV]</th><th><i>w</i><sub>0</sub></th><th><i>w</i><sub>a</sub></th><th>FoM</th></tr></thead><tbody><tr><td colspan="10" style="text-align:center; font-weight:bold;">ΛCDM</td></tr><tr><td rowspan="2">WL</td><td>0.141</td><td>0.00494</td><td>0.0244</td><td>0.00708</td><td>0.0327</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>0.136</td><td>0.00339</td><td>0.0236</td><td>0.00528</td><td>0.0298</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td rowspan="2">V</td><td>0.215</td><td>0.252</td><td>0.0625</td><td>0.735</td><td>0.581</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>0.0398</td><td>0.0198</td><td>0.0145</td><td>0.110</td><td>0.103</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td rowspan="2">WL + V</td><td>0.107</td><td>0.00483</td><td>0.0187</td><td>0.00655</td><td>0.0257</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>0.00735</td><td>0.00304</td><td>0.00389</td><td>0.00363</td><td>0.00432</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td rowspan="2">WL + V + XC</td><td>0.0216</td><td>0.00466</td><td>0.00530</td><td>0.00575</td><td>0.0104</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>0.00604</td><td>0.00285</td><td>0.00321</td><td>0.00338</td><td>0.00409</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td colspan="10" style="text-align:center; font-weight:bold;">νΛCDM</td></tr><tr><td rowspan="2">WL</td><td>0.144</td><td>0.00526</td><td>0.0265</td><td>0.0129</td><td>0.0348</td><td>0.310</td><td>—</td><td>—</td><td>—</td></tr><tr><td>0.138</td><td>0.00342</td><td>0.0248</td><td>0.00873</td><td>0.0304</td><td>0.229</td><td>—</td><td>—</td><td>—</td></tr><tr><td rowspan="2">V</td><td>0.216</td><td>0.254</td><td>0.0667</td><td>0.812</td><td>0.707</td><td>2.46</td><td>—</td><td>—</td><td>—</td></tr><tr><td>0.0402</td><td>0.0205</td><td>0.0154</td><td>0.113</td><td>0.131</td><td>0.593</td><td>—</td><td>—</td><td>—</td></tr><tr><td rowspan="2">WL + V</td><td>0.109</td><td>0.00519</td><td>0.0203</td><td>0.0127</td><td>0.0272</td><td>0.302</td><td>—</td><td>—</td><td>—</td></tr><tr><td>0.00741</td><td>0.00309</td><td>0.00579</td><td>0.00752</td><td>0.00461</td><td>0.202</td><td>—</td><td>—</td><td>—</td></tr><tr><td rowspan="2">WL + V + XC</td><td>0.0216</td><td>0.00509</td><td>0.00718</td><td>0.0124</td><td>0.0118</td><td>0.292</td><td>—</td><td>—</td><td>—</td></tr><tr><td>0.00605</td><td>0.00289</td><td>0.00505</td><td>0.00705</td><td>0.00443</td><td>0.193</td><td>—</td><td>—</td><td>—</td></tr><tr><td colspan="10" style="text-align:center; font-weight:bold;">w<sub>0</sub>w<sub>a</sub>CDM</td></tr><tr><td rowspan="2">WL</td><td>0.141</td><td>0.0130</td><td>0.0246</td><td>0.0152</td><td>0.0348</td><td>—</td><td>0.147</td><td>0.552</td><td>26.9</td></tr><tr><td>0.138</td><td>0.00949</td><td>0.0237</td><td>0.0104</td><td>0.0299</td><td>—</td><td>0.121</td><td>0.426</td><td>53.9</td></tr><tr><td rowspan="2">V</td><td>0.316</td><td>0.394</td><td>0.0793</td><td>0.804</td><td>0.674</td><td>—</td><td>2.30</td><td>6.59</td><td>0.253</td></tr><tr><td>0.0636</td><td>0.0934</td><td>0.0181</td><td>0.113</td><td>0.166</td><td>—</td><td>0.390</td><td>1.55</td><td>5.11</td></tr><tr><td rowspan="2">WL + V</td><td>0.108</td><td>0.0128</td><td>0.0191</td><td>0.0148</td><td>0.0279</td><td>—</td><td>0.146</td><td>0.543</td><td>27.7</td></tr><tr><td>0.0131</td><td>0.00744</td><td>0.00407</td><td>0.00748</td><td>0.00461</td><td>—</td><td>0.0912</td><td>0.296</td><td>84.7</td></tr><tr><td rowspan="2">WL + V + XC</td><td>0.0244</td><td>0.0124</td><td>0.00638</td><td>0.0139</td><td>0.0143</td><td>—</td><td>0.144</td><td>0.526</td><td>29.3</td></tr><tr><td>0.0105</td><td>0.00649</td><td>0.00348</td><td>0.00650</td><td>0.00425</td><td>—</td><td>0.0791</td><td>0.251</td><td>106</td></tr><tr><td colspan="10" style="text-align:center; font-weight:bold;">νw<sub>0</sub>w<sub>a</sub>CDM</td></tr><tr><td rowspan="2">WL</td><td>0.148</td><td>0.0148</td><td>0.0270</td><td>0.0273</td><td>0.0430</td><td>0.413</td><td>0.149</td><td>0.630</td><td>20.2</td></tr><tr><td>0.139</td><td>0.00950</td><td>0.0250</td><td>0.0120</td><td>0.0305</td><td>0.231</td><td>0.122</td><td>0.429</td><td>53.6</td></tr><tr><td rowspan="2">V</td><td>0.351</td><td>0.482</td><td>0.0801</td><td>1.26</td><td>0.710</td><td>4.21</td><td>2.43</td><td>7.98</td><td>0.148</td></tr><tr><td>0.0676</td><td>0.0987</td><td>0.0182</td><td>0.119</td><td>0.169</td><td>0.655</td><td>0.409</td><td>1.67</td><td>4.63</td></tr><tr><td rowspan="2">WL + V</td><td>0.110</td><td>0.0146</td><td>0.0204</td><td>0.0267</td><td>0.0341</td><td>0.397</td><td>0.148</td><td>0.614</td><td>20.9</td></tr><tr><td>0.0131</td><td>0.00747</td><td>0.00584</td><td>0.0102</td><td>0.00492</td><td>0.203</td><td>0.0912</td><td>0.296</td><td>84.4</td></tr><tr><td rowspan="2">WL + V + XC</td><td>0.0244</td><td>0.0142</td><td>0.00721</td><td>0.0259</td><td>0.0200</td><td>0.380</td><td>0.145</td><td>0.591</td><td>22.6</td></tr><tr><td>0.0105</td><td>0.00652</td><td>0.00514</td><td>0.00929</td><td>0.00463</td><td>0.182</td><td>0.0791</td><td>0.251</td><td>105</td></tr></tbody></table>

**Notes.** In each probe block, the first row shows the errors for the pessimistic scenario, while the second row corresponds to the optimistic one.

$w_0$ and $w_a$ also enter the linear growth factor (Linder & Jenkins 2003), and therefore impact $\Omega_m$ and the normalization $\sigma_8$ again. When moving from the pessimistic to the optimistic scenario, the constraints on $w_0$ and $w_a$ improve, and consequently the FoM is enhanced by a factor of ~3. This can be explained since, in the optimistic scenario, we exploited both the cosmology dependence of the growth factor in the void bias evolution, $b_v$, as well as the cosmology dependence of the void size function used to compute the effective void bias in Eq. (19). This is confirmed by previous works (Pisani et al. 2015; Verza et al. 2019), although these focused on upcoming spectroscopic data and found that the void size function is sensitive to the DE equation of state. We also stress the increase in the constraining power when the two probes are combined together: comparing from Table 5 the constraints produced by WL alone against the ones obtained from its combination with the void clustering and void-lensing cross-correlation, we can observe that the FoM is enhanced by ~10% in the pessimistic scenario, and by a factor of ~2 in the optimistic one.

Finally, in order to measure the impact on parameter forecasts of the void-lensing cross-correlation signal, we also evaluated the marginalized errors when WL and void clustering are treated as independent probes, that is, when the Fisher matrices of the single probes are directly summed up. Then, we compared the FoM against the case when the void-lensing cross-correlation is included in the analysis. We found that, in the latter case, the DE FoM is enhanced by 5% in the pessimistic scenario, and by 20% in the optimistic one, as reported in Table 5.

## 5.4. Including the spectroscopic galaxy clustering

In this section we combine the spectroscopic galaxy clustering, $GC_{sp}$, to the WL, void clustering, and void-lensing cross-correlation probes. To this purpose, we considered $GC_{sp}$ as a probe independent of the others, and made use of the $GC_{sp}$ Fisher matrix as provided by EC20 for the observed anisotropic galaxy power spectrum, adding it to the ones computed in this work. Here our goal was to evaluate if and by how much the void
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Fisher matrix marginalized contours for the (baseline) ΛCDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios.</figcaption></figure>

clustering and the void-lensing correlation can still improve the *Euclid* performance even when its two primary probes are both accounted for. The corresponding results are reported in Fig. 7 and summarized in Table 6.

In the pessimistic scenario, constraints given by $GC_{sp} + V$ are weaker than the ones given by $GC_{sp} + WL$; for instance, the FoM has a value of 15 in the former case, and of 112 in the latter one. We find that the constraints given by $GC_{sp} + WL + V$ are only slightly improved with respect to $GC_{sp} + WL$. Nevertheless, the constraints given by $GC_{sp} + WL + V + XC$ are tighter than in the $GC_{sp} + WL$ case. In particular, the FoM increases by ~5% when considering $GC_{sp} + WL + V + XC$ against $GC_{sp} + WL$.

In the optimistic scenario, constraints given by $GC_{sp} + V$ are comparable to the ones given by $GC_{sp} + WL$. Moreover, the FoM increases from 192 for $GC_{sp} + WL$, up to 702 for $GC_{sp} + WL + V$, in which case constraints on the other cosmological parameters become tighter: by ~10% for $h$, by ~30% for $n_s$, by ~50–55% for $\Omega_b$ and $\sigma_8$, and by ~70% for $\Omega_m$.

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Fisher matrix marginalized contours for the $\nu\Lambda$CDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios.</figcaption></figure>

Finally, when also including the void-lensing cross-correlation, while in the pessimistic case there is at most a 10% improvement in the FoM, in the optimistic setup we find that the most improved parameters are $\Omega_m$ and $\sigma_8$, with their constraints improved by a factor of 2 and 3, respectively; the total FoM increases by ~10%, that is, from 702 for $GC_{sp} + WL + V$, up to 791 for $GC_{sp} + WL + V + XC$.

### 5.5. Combining massive neutrinos with dynamical dark energy

In this section we consider $M_\nu$, $w_0$, and $w_a$ as free parameters. The contour plots for this model cosmology are reported in Fig. 8, both for the pessimistic and optimistic cases. As expected, in this scenario the constraint on $\sigma_8$ gets even weaker with respect to the baseline ΛCDM case, due to the combined variation of $M_\nu$, $w_0$, and $w_a$, entering the matter power spectrum, the Hubble parameter, and the growth factor. Comparing
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Fisher matrix marginalized contours for the $w_0w_a$ CDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios.</figcaption></figure>

the pessimistic and optimistic scenarios, we observe improvements on the constraints in the latter case, but the orientation of the ellipses remains unchanged: this is reassuring as it means that the results do not change qualitatively when the two different setups are considered. In addition, we also performed a forecast for this model with a fiducial neutrino mass of $M_ν = 0.33$ eV with a degenerate mass spectrum, using the best-fit value in Pellejero-Ibanez et al. (2017): with this choice, the marginalized errors on $M_ν$ decrease by ~10% with respect to the fiducial case with $M_ν = 0.06$ eV. Again, this can be explained since, by increasing the value of $M_ν$, we increase the impact of neutrinos on the matter power spectrum, which suffers a larger suppression with respect to the fiducial neutrino mass case, becoming more sensitive to the presence of neutrinos.

From Table 5, it is possible to observe that, when we consider the case including the void-lensing cross-correlation signal compared to the one where the two probes are joined independently,

the error on $M_ν$ gets reduced by ~5% and ~10%, and the FoM is enhanced by 10% and 25%, in the pessimistic and optimistic scenarios, respectively.

As a final consideration, we comment on the mutual impact of $M_ν$ and DE. Including $w_0$ and $w_a$ increases the error on $M_ν$ and vice versa. The increase in the $w_0$ and $w_a$ errors, due to the variation of $M_ν$, is at most ~15%, while the error on $M_ν$ increases at most by 70% due to the variation of $w_0$ and $w_a$. However, especially for $w_0$ and $w_a$, the constraints are less impacted when the void-lensing cross-correlation is included in the analysis, compared to the constraints given by the single probes alone, as shown in the top panels of Fig. 9.

## 5.6. Systematics checks

The analysis presented in this work relies on various choices for the computation of the Fisher matrices, such as the number of bins in multipoles $ℓ$, redshift $z$, and wave-number $k$, as well as on differentiation methods. In this section we present the stability of the final results against different computation choices.

We performed stability tests where the multipole range was divided into 80, 100, and 120 logarithmically equispaced bins, where the $C(ℓ)$ were evaluated in the linear center of the bin. The impact of this $ℓ$-binning was negligible: the marginalized 1-$\sigma$ errors stayed unchanged up to the third digit. We also performed a forecast with a reduced multipole range, from $ℓ = 20$ to $ℓ = 1350$. Both for the pessimistic and optimistic void bias scenarios, the impact of the multipole range reduction on the FoM and $M_ν$ 1-$\sigma$ uncertainty was about 10%.

The impact of the $k$-binning was negligible too. CAMB evaluates the matter power spectrum $P_{\text{mm}}(k, z)$ with a logarithmic $k$-binning grid. For the reported forecasts, $P_{\text{mm}}(k, z)$ was evaluated with 60 $k$ values per decade. Stability tests were performed also using 50 and 90 $k$ values per decade; the marginalized 1-$\sigma$ errors remained unchanged up to the third digit.

The impact of the $z$-binning was minor. CAMB evaluates the matter power spectrum $P_{\text{mm}}(k, z)$ with a linear $z$-binning grid. For the reported forecast, $P_{\text{mm}}(k, z)$ was evaluated for 450 $z$ values. Stability tests were also performed using 300 and 600 $z$-bins: the differences in the marginalized 1-$\sigma$ errors were always smaller than 3%.

The impact of the differentiation method was also negligible. In the reported forecast, the SteM derivative method was used. A forecast was also performed with the semi-analytical derivative method; the marginalized 1-$\sigma$ errors stayed unchanged up to the fourth digit.

We also performed a forecast using equi-populated redshift bins for voids, rather than the equi-populated redshift bins for galaxies. This modification led to slightly (about 10–15%) worse constraints. This can be understood looking at the void distribution in Fig. 2, which increases with the redshift $z$. Thus the requirement of equi-populated redshift bins translates into broader (and so sparser) bins at low $z$, and narrower (and so denser) bins at high $z$, with respect to the equi-populated galaxy configuration. Since the effective void bias decreases in absolute value with $z$, this leads to more bins with a lower void bias, resulting in a decrease in the $C(ℓ)$ in those bins, and hence to weaker constraints, with, as expected, a larger impact in the optimistic case than in the pessimistic case, where we marginalize over the void bias. Overall we conclude that our results are robust against various computation choices.
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Full contour plots when $GC_{sp}$ is combined with WL + V + XC for the $w_0w_a$ CDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios.</figcaption></figure>

**6. Conclusion and outlook**

In this work we present the first forecast on cosmological parameter inference obtained combining the void-lensing cross-correlation with the WL and void angular two-point correlations,

as will be measured from the *Euclid* photometric galaxy catalogs. In order to reach this goal we considered the following approach:

- we measured the projected void density distribution in redshift, $n^v(z)$, directly from the *Euclid* Flagship photometric
<figure><img src="image_9.png" /><figcaption>Fig. 8. Fisher matrix marginalized contours for the $\nu w_0 w_a$ CDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Forecast contours for the $vw_0w_a$ CDM model, in the pessimistic (**left**) and in the optimistic (**right**) scenarios.</figcaption></figure>

<table><caption>**Table 6.** Marginalized $1-\sigma$ errors in different cosmological scenarios for WL + V + GC\_sp, when they are assumed to be independent, together with their combinations when the void-lensing cross-correlation is included.</caption><thead><tr><th>Probe</th><th><i>h</i></th><th>Ω<sub>m</sub></th><th>Ω<sub>b</sub></th><th>σ<sub>8</sub></th><th><i>n</i><sub>s</sub></th><th><i>M</i><sub>v</sub> [eV]</th><th><i>w</i><sub>0</sub></th><th><i>w</i><sub>a</sub></th><th>FoM</th></tr></thead><tbody><tr><td colspan="10" style="text-align:center;"><i>w</i><sub>0</sub><i>w</i><sub>a</sub> CDM</td></tr><tr><td rowspan="2">WL + V + GC<sub>sp</sub></td><td>0.00435</td><td>0.00931</td><td>0.00301</td><td>0.00947</td><td>0.00832</td><td>—</td><td>0.0923</td><td>0.31</td><td>113</td></tr><tr><td>0.00133</td><td>0.00395</td><td>0.00128</td><td>0.0065</td><td>0.0036</td><td>—</td><td>0.0256</td><td>0.196</td><td>702</td></tr><tr><td rowspan="2">WL + V + XC + GC<sub>sp</sub></td><td>0.00420</td><td>0.00856</td><td>0.00295</td><td>0.0092</td><td>0.00811</td><td>—</td><td>0.0883</td><td>0.295</td><td>124</td></tr><tr><td>0.0013</td><td>0.0018</td><td>0.00088</td><td>0.00205</td><td>0.00332</td><td>—</td><td>0.0228</td><td>0.113</td><td>791</td></tr></tbody></table>

**Notes.** In each probe block, the first row shows the errors for the pessimistic scenario, while the second row corresponds to the optimistic one.

mock galaxy catalog, using the 2D void finder of Sánchez et al. (2017);

- we evaluated the $C(l)$ numerical derivatives, which enter the Fisher matrix expression, using the SteM technique (Camera et al. 2017). We tested the stability of results against different differentiation techniques, and found them to be robust;

- we used the void bias as obtained from the PBS and excursion set formalisms (Sheth & Van De Weygaert 2004), and considered the volume-conserving $Vdn$ void size function model (Jennings et al. 2013).

We present our results for two different choices of the effective void bias $b_v^{\text{eff}}(z)$: an “optimistic scenario”, where we exploited the cosmology dependence of both the void size function and void bias, and a “pessimistic scenario”, where we assumed the void bias evolution to be given by the growth factor in the reference cosmology, but its absolute normalization was supposed to be unknown, so that $b_v^{\text{eff}}(z = 0)$ was marginalized over as a nuisance parameter.

We present parameter forecasts for different cosmological models: starting from a flat $\Lambda$CDM model, we first separately added as free parameters the total neutrino mass and the CPL parametrization of the DE equation of state (Chevallier & Polarski 2001; Linder 2002), then we add both in combination. Our main findings are presented in Table 5 and Figs. 4–8, and can be summarized as follows:

- $\Lambda$CDM cosmology: the WL angular correlation function is able to constrain all the cosmological parameters better than

the photometric void angular spectrum; however, void constraints on $h$ and $\Omega_b$ are competitive in the optimistic scenario. For the Hubble constant, $h$, this is due to the form of the weak lensing kernel and the integration along the line of sight. For $\Omega_b$ this is due to the presence of the BAO features in the void angular spectrum, which on the other hand are washed out in the lensing angular spectrum (see, e.g., Fig. 1). In general, these, together with $n_s$, are the parameters for which the constraints improve most when void clustering is combined with WL.

- $\Lambda$CDM cosmology: adding to the analysis the total neutrino mass, $M_v$, as a free parameter mainly impacts the constraints on $\sigma_8$, as expected, since $\sigma_8$ and $M_v$ both affect the amplitude of the matter power spectrum (even if in the $M_v$-case this is a scale-dependent effect). The void clustering, together with the void-lensing cross-correlation, improves the determination of the neutrino mass scale with respect to WL alone: the constraint on $M_v$ becomes stronger by ~5% and ~15%, respectively, in the pessimistic and optimistic scenarios.

- $w_0w_a$CDM cosmology: when the CPL parametrization of the DE equation of state is considered, we find that the most affected parameters are $h$, $\Omega_m$, and $\sigma_8$. The expression of the Hubble parameter Eq. (22) explains the effect on $h$ and $\Omega_m$, while the impact of $w_0$ and $w_a$ on the linear growth factor (Linder & Jenkins 2003) accounts for the effect on $\sigma_8$ and, again, $\Omega_m$. Also, in this scenario void clustering helps to tighten the constraints on the cosmological parameters. In particular, the DE FoM increases by ~10% in the pessimistic scenario and by a factor of ~2 in the optimistic scenario.
- *vw*₀*w*ₐCDM cosmology: in this model we evaluated, in particular, the mutual impact among *w*₀, *w*ₐ, and *M*ᵥ. On the one hand, we find that, when adding *M*ᵥ, the marginalized errors of *w*₀ and *w*ₐ increase at most by ~15%. On the other hand, when adding *w*₀, *w*ₐ as free parameters, the impact on the *M*ᵥ constraints is higher, as its error increases even by 70%. However, especially for *w*₀ and *w*ₐ, the constraints are less impacted when the void-lensing cross-correlation is included, compared to the constraints provided by the two probes combined independently.

- *FoM and the void-lensing cross-signal:* in order to evaluate the improvement in the constraining power provided by adding the void-lensing cross-correlation signal to void clustering and weak lensing, we also considered the case where these two probes were assumed to be independent. When including the void-lensing correlation, the error on *M*ᵥ was reduced by 5% and 10%, in the pessimistic and optimistic scenarios, respectively. The FoM was enhanced by 10% and 25%, in the pessimistic and optimistic scenarios, respectively.

- *Adding* GC<sub>sp</sub>: finally, we combined the galaxy lensing, void clustering, and void-lensing cross-correlation probes with the spectroscopic galaxy clustering, GC<sub>sp</sub>. In this case, we exploited the Fisher matrices as obtained by EC20 in the so-called “pessimistic” and “optimistic” settings, and considered the flat *w*₀*w*ₐCDM scenario, since GC<sub>sp</sub> forecasts are computed keeping *M*ᵥ = 0.06 eV fixed. When combining WL + V + XC with GC<sub>sp</sub>, both assumed in the pessimistic configuration, we find that the FoM increases from 23 for WL + V + XC alone up to 117 for WL + V + XC + GC<sub>sp</sub>, and from 105 up to 791 in the optimistic case. In particular in the latter case, results are promising and competitive with other kinds of probe combination, as, for example, when galaxy lensing is combined with photometric galaxy clustering, GC<sub>ph</sub>, and their cross-correlation (Tutusaus et al. 2020). As intrinsic alignments are not included in this analysis, we cannot make a proper comparison with the results reported in Table 14 in EC20, where the authors show a total FoM of 377 and 1257 in the pessimistic and optimistic settings, respectively. The results are, nevertheless, extremely encouraging: we show that including void clustering and the void-lensing cross-correlation in the total likelihood analysis will allow us to considerably improve Euclid's performance.

The forecasts presented in this work show that photometric void clustering and its cross-correlation with WL deserve to be exploited in the data analysis of the *Euclid* galaxy survey, as they could be able to improve constraints on several cosmological parameters, such as, in particular, the total neutrino mass (Kreisch et al. 2021) and the DE equation of state. In this respect, it is worth noting that, for a full comparison with the *Euclid* performance from primary probes, in our analysis we should have included also the information from photometric galaxy clustering. However, we decided not to add this probe to the analysis in order to avoid double counting the information from the photometric sample. A way to include GC<sub>ph</sub>, without incurring in this issue, is offered by a full-field inference approach, which allows us to self-consistently analyze GC and voids in the DM distribution (Leclercq et al. 2015). It would be important to also include GC<sub>ph</sub> since, as shown in this work, the combination of void-clustering, void-lensing, and GC<sub>sp</sub> does not allow us to reach the same *Euclid* performance as the combination of *Euclid* primary probes (EC20). However, this does not imply that the information coming from the photometric void sample should be neglected, as voids and galaxies are tracers with a very different

bias, and hence probe different regimes of LSS (Wang & Zhao 2020), therefore providing therefore complementary information about our Universe. Moreover, in this paper we only consider only some observational and astrophysical uncertainties, neglecting baryonic physics and intrinsic alignment of galaxies.

While the aim of this work was to show that the inclusion of the void-void auto-correlation and void-lensing cross-correlation will improve the *Euclid* survey performance, and we performed some checks against the Flagship simulation to ensure the reliability of our theoretical modeling, additional work needs to be carried out to prepare the pipelines for the analyses of forthcoming real data. First and foremost, the theoretical model should be further improved. As already mentioned, we did not include nonlinear scales, since this would require a modeling of the void density profile and its dependence on cosmological parameters. Second, we need to rely on more realistic mock data, which will account for more effects: *Euclid* survey specification and systematics and possibly misSpecifications induced by the HOD, which could induce some errors in the identified void catalogs. Finally, the framework presented here can also be easily extended to investigate non-flat scenarios, more specific DE models, theories of modified gravity, and primordial non-Gaussianities (Chan et al. 2019).

*Acknowledgements.* M.B. acknowledges financial support from the ASI agreement n. I/023/12/0 "Euclid attività relativa alla fase B2/C". N.H. is supported by the Excellence Cluster ORIGINS, which is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy - EXC-2094 - 390783311. A.P. is supported by NASA ROSES grant 12-EUCLID12-0004, and NASA grant 15-WFIRST15-0008 to the *Nancy Grace Roman Space Telescope Science Investigation Team* "Cosmology with the High Latitude Survey". A.K. has been supported by a Juan de la Cierva Incaperación fellowship with project number IJC2018-037730-I, and funding for this project was also available in part through SEV-2015-0548 and AYA2017-89891-P. S.N. acknowledges support from an STFC Ernest Rutherford Fellowship, grant reference ST/T005009/1. The Euclid Consortium acknowledges the European Space Agency and a number of agencies and institutes that have supported the development of *Euclid*, in particular the Academy of Finland, the Agenzia Spaziale Italiana, the Belgian Science Policy, the Canadian Euclid Consortium, the French Centre National d'Etudes Spatiales, the Deutsches Zentrum für Luft- und Raumfahrt, the Danish Space Research Institute, the Fundação para a Ciência e a Tecnologia, the Ministerio de Economia y Competitividad, the National Aeronautics and Space Administration, the National Astronomical Observatory of Japan, the Netherlandse Onderzoeksg school Voor Astronomie, the Norwegian Space Agency, the Romanian Space Agency, the State Secretariat for Education, Research and Innovation (SERI) at the Swiss Space Office (SSO), and the United Kingdom Space Agency. A complete and detailed list is available on the *Euclid* website (http://www.euclid-ec.org).

## References

- Armendariz-Picon, C., Mukhanov, V., & Steinhardt, P. J. 2000, Phys. Rev. Lett., 85, 4438

- Aubert, M., Cousinou, M.-C., Escoffier, S., et al. 2020, MNRAS, 513, 186

- Aheroozi, P. S., Wechsler, R. H., Wu, H.-Y., et al. 2013, ApJ, 762, 109

- Bird, S., Viel, M., & Haehnelt, M. G. 2011, MNRAS, 420, 2551

- Blanton, M. R., Hogg, D. W., Bahcall, N. A., et al. 2003, ApJ, 592, 819

- Blanton, M. R., Lupton, R. H., Schlegel, D. J., et al. 2005, ApJ, 631, 208

- Bos, E. G., van de Weygaert, R., Dolag, K., & Pettorino, V. 2012, MNRAS, 426, 440

- Brew, E., & Atkinson, K. E. 1979, The Mathematical Gazette, 63

- Brinckmann, T., & Lesgourgues, J. 2019, Phys. Dark Univ., 24, 100260

- Camera, S., Harrison, I., Bonaldi, A., & Brown, M. L. 2017, MNRAS, 464, 4747

- Carretero, J., Castander, F. J., Gaztañaga, E., Crocce, M., & Fosalba, P. 2015, MNRAS, 447, 646

- Chan, K. C., Hamaus, N., & Desjacques, V. 2014, Phys. Rev. D, 90, 103521

- Chan, K. C., Hamaus, N., & Biagetti, M. 2019, Phys. Rev. D, 99, 121304

- Chevallier, M., & Polarski, D. 2001, Int. J. Mod. Phys., D, 10, 213

- Clampitt, J., & Jain, B. 2015, MNRAS, 454, 3357

- Contarini, S., Ronconi, T., Marulli, F., et al. 2019, MNRAS, 488, 3526

- Contarini, S., Verza, G., Pisani, A., et al. 2022, A&A, 667, A162

- Cooray, A., & Sheth, R. 2002, Phys. Rep., 372, 1
Copeland, E. J., Sami, M., & Tsujikawa, S. 2006, Int. J. Mod. Phys. D, 15, 1753

Costille, A., Caillat, A., Rossin, C., et al. 2019, SPIE Conf. Ser., 11180, 1118013

Crocce, M., Castander, F. J., Gaztañaga, E., Fosalba, P., & Carretero, J. 2015, MNRAS, 453, 1513

Cropper, M., Pottinger, S., Niemi, S., et al. 2016, SPIE Conf. Ser., 9904, 99040Q

Euclid Collaboration (Blanchard, A., et al.) 2020, A&A, 642, A191

Fang, Y., Hamaus, N., Jain, B., et al. 2019, MNRAS, 490, 3573

Fisher, R. A. 1935, J. R. Stat. Soc., 98, 39

Granett, B. R., Neyrinck, M. C., & Szapudi, I. 2008, ApJ, 683, L99

Granett, B. R., Kovács, A., & Hawken, A. J. 2015, MNRAS, 454, 2804

Hamaus, N., Wandelt, B. D., Sutter, P. M., Lavaux, G., & Warren, M. S. 2014a, Phys. Rev. Lett., 112, 041304

Hamaus, N., Sutter, P. M., & Wandelt, B. D. 2014b, Phys. Rev. Lett., 112, 251302

Hamaus, N., Pisani, A., Sutter, P. M., et al. 2016, Phys. Rev. Lett., 117, 091302

Hamaus, N., Cousinou, M.-C., Pisani, A., et al. 2017, JCAP, 07, 014

Hamaus, N., Pisani, A., Choi, J.-A., et al. 2020, JCAP, 12, 023

Hamaus, N., Aubert, M., Pisani, A., et al. 2021, A&A, 658, A20

Hand, N., Feng, Y., Beutler, F., et al. 2018, AJ, 156, 160

Hawken, A. J., Aubert, M., Pisani, A., et al. 2020, JCAP, 06, 012

Hu, W. 1999, ApJ, 522, L21

Ilić, S., Langer, M., & Douspis, M. 2013, A&A, 556, A51

Jamieson, D., & Loverde, M. 2019, Phys. Rev. D, 100, 123528

Jennings, E., Li, Y., & Hu, W. 2013, MNRAS, 434, 2167

Jimenez, R., Kitching, T., Pena-Garay, C., & Verde, L. 2010, JCAP, 05, 035

Kiakotou, A., Elgarøy, Ø., & Lahav, O. 2008, Phys. Rev. D, 77, 063005

Krause, E., Chang, T.-C., Doré, O., & Umetsu, K. 2013, ApJ, 762, L20

Kreisch, C. D., Pisani, A., Carbone, C., et al. 2019, MNRAS, 488, 4413

Kreisch, C. D., Pisani, A., Villaescusa-Navarro, F., et al. 2021, ApJ, 935, 100

Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, ArXiv e-prints [arXiv:1110.3193]

Lavaux, G., & Wandelt, B. D. 2010, MNRAS, 403, 1392

Lavaux, G., & Wandelt, B. D. 2012, ApJ, 754, 109

Leclercq, F., Jasche, J., Sutter, P. M., Hamaus, N., & Wandelt, B. 2015, JCAP, 03, 047

Lee, J., & Park, D. 2009, ApJ, 696, L10

Lesgourgues, J., & Pastor, S. 2006, Phys. Rep., 429, 307

Lesgourgues, J., & Pastor, S. 2012, Adv. High Energy Phys., 2012, 608515

Lesgourgues, J., & Pastor, S. 2014, New J. Phys., 16, 065002

Lewis, A., Challinor, A., & Lasenby, A. 1999, ApJ, 538, 473

Limber, D. N. 1953, ApJ, 117, 134

Linder, E. V. 2002, Phys. Rev. Lett., 90, 4

Linder, E. V., & Jenkins, A. 2003, MNRAS, 346, 573

LoVerde, M., & Afshordi, N. 2008, Phys. Rev. D, 78, 123506

Mao, Q., Berlind, A. A., Scherrer, R. J., et al. 2017, ApJ, 835, 161

Massara, E., Villaescusa-Navarro, F., Viel, M., & Sutter, P. M. 2015, JCAP, 11, 018

Melchior, P., Sutter, P. M., Sheldon, E. S., Krause, E., & Wandelt, B. D. 2014, MNRAS, 440, 2922

Micheletti, D., Iovino, A., Hawken, A. J., et al. 2014, A&A, 570, A106

Nadathur, S., Carter, P. M., Percival, W. J., Winther, H. A., & Bautista, J. E. 2019, Phys. Rev. D, 100, 023504

Nadathur, S., Woodfinden, A., Percival, W. J., et al. 2020, MNRAS, 499, 4140

Pellejero-Ibanez, M., Chuang, C.-H., Rubiño-Martín, J. A., et al. 2017, MNRAS, 468, 4116

Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1998, ApJ, 517, 565

Pisani, A., Sutter, P. M., Hamaus, N., et al. 2015, Phys. Rev. D, 92, 083531

Pisani, A., Massara, E., Spergel, D. N., et al. 2019, BAAS, 51, 40

Planck Collaboration XIII. 2016, A&A, 594, A13

Pollina, G., Hamaus, N., Paech, K., et al. 2019, MNRAS, 487, 2836

Potter, D., Stadel, J., & Teyssier, R. 2017, Comput. Astrophys. Cosmol., 4, 2

Ratra, B., & Peebles, P. J. 1988, Phys. Rev. D, 37, 3406

Riess, A. G., Filippenko, A. V., Challis, P., et al. 1998, AJ, 116, 1009

Ronconi, T., & Marulli, F. 2017, A&A, 607, A24

Sánchez, C., Clampitt, J., Kovacs, A., et al. 2017, MNRAS, 465, 746

Schuster, N., Hamaus, N., Pisani, A., et al. 2019, JCAP, 12, 055

Sheth, R. K., & Van De Weygaert, R. 2004, MNRAS, 350, 517

Steeb, W.-H. H. 1997, Matrix Calculus and the Kronecker Product with Applications and C++ Programs (USA: World Scientific Publishing Co., Inc.)

Steinhardt, P. J., Wang, L., & Zlatev, I. 1999, Phys. Rev. D, 59, 123504

Stopyra, S., Peiris, H. V., & Pontzen, A. 2020, MNRAS, 500, 4173

Sutter, P. M., Lavaux, G., Wandelt, B. D., & Weinberg, D. H. 2012, ApJ, 761, 44

Sutter, P. M., Pisani, A., Wandelt, B. D., & Weinberg, D. H. 2014, MNRAS, 443, 2983

Takahashi, R., Sato, M., Nishimichi, T., Taruya, A., & Oguri, M. 2012, ApJ, 761, 152

Tutusaus, I., Martinelli, M., Cardone, V. F., et al. 2020, A&A, 643, A70

Verza, G., Pisani, A., Carbone, C., Hamaus, N., & Guzzo, L. 2019, JCAP, 12, 040

Vielzeuf, P., Kovács, A., Demirbozan, U., et al. 2019, MNRAS, 500, 464

Voivodic, R., Rubira, H., & Lima, M. 2020, JCAP, 10, 033

Wang, Y., & Zhao, G.-B. 2020, RAA, 20, 158

Zehavi, I., Zheng, Z., Weinberg, D. H., et al. 2011, ApJ, 736, 59

¹ Dipartimento di Fisica, Università degli Studi di Genova, and INFN-Sezione di Genova, Via Dodecaneso 33, 16146 Genova, Italy

² INAF-IASF Milano, Via Alfonso Corti 12, 20133 Milano, Italy e-mail: marco.bonici@inaf.it

³ INFN-Sezione di Genova, Via Dodecaneso 33, 16146 Genova, Italy

⁴ Aix-Marseille Univ., CNRS/IN2P3, CPPM, Marseille, France

⁵ INAF-Osservatorio Astronomico di Roma, Via Frascati 33, 00078 Monteporzio Catone, Italy

⁶ INFN-Sezione di Roma, Piazzale Aldo Moro, 2 – c/o Dipartimento di Fisica, Edificio G. Marconi, 00185 Roma, Italy

⁷ Universitäts-Sternwarte München, Fakultät für Physik, Ludwig-Maximilians-Universität München, Scheinerstrasse 1, 81679 München, Germany

⁸ Department of Astrophysical Sciences, Peyton Hall, Princeton University, Princeton, NJ 08544, USA

⁹ Instituto de Astrofísica de Canarias (IAC); Departamento de Astrofísica, Universidad de La Laguna (ULL), 38200 La Laguna, Tenerife, Spain

¹⁰ Instituto de Astrofísica de Canarias, Calle Vía Láctea s/n, 38204 San Cristóbal de La Laguna, Tenerife, Spain

¹¹ Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK

¹² Dipartimento di Fisica e Astronomia “Augusto Righi” – Alma Mater Studiorium Università di Bologna, Via Piero Gobetti 93/2, 40129 Bologna, Italy

¹³ INAF-Osservatorio di Astrofísica e Scienza dello Spazio di Bologna, Via Piero Gobetti 93/3, 40129 Bologna, Italy

¹⁴ INFN-Sezione di Bologna, Viale Berti Pichat 6/2, 40127 Bologna, Italy

¹⁵ Dipartimento di Fisica e Astronomia “G. Galilei”, Università di Padova, Via Marzolo 8, 35131 Padova, Italy

¹⁶ INFN-Padova, Via Marzolo 8, 35131 Padova, Italy

¹⁷ Université de Genève, Département de Physique Théorique and Centre for Astroparticle Physics, 24 Quai Ernest-Ansermet, 1211 Genève 4, Switzerland

¹⁸ Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Can Magrans, s/n, 08193 Barcelona, Spain

¹⁹ Institut d'Estudis Espacials de Catalunya (IEEC), Carrer Gran Capitá 2-4, 08034 Barcelona, Spain

²⁰ Institut de Recherche en Astrophysique et Planétologie (IRAP), Université de Toulouse, CNRS, UPS, CNES, 14 Av. Edouard Belin, 31400 Toulouse, France

²¹ University of Lyon, UCB Lyon 1, CNRS/IN2P3, IUF, IP2I Lyon, France

²² Istituto Nazionale di Astrofisica (INAF) – Osservatorio di Astrofisica e Scienza dello Spazio (OAS), Via Gobetti 93/3, 40127 Bologna, Italy

²³ Istituto Nazionale di Fisica Nucleare, Sezione di Bologna, Via Irne-rio 46, 40126 Bologna, Italy

²⁴ Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK

²⁵ Dipartimento di Fisica, Università degli Studi di Torino, Via P. Giuria 1, 10125 Torino, Italy

²⁶ INFN-Sezione di Torino, Via P. Giuria 1, 10125 Torino, Italy

²⁷ INAF-Osservatorio Astrofisico di Torino, Via Osservatorio 20, 10025 Pino Torinese, TO, Italy

²⁸ Instituto de Física Teórica UAM-CSIC, Campus de Cantoblanco, 28049 Madrid, Spain

²⁹ AIM, CEA, CNRS, Université Paris-Saclay, Université de Paris, 91191 Gif-sur-Yvette, France
30 Université Paris-Saclay, Université Paris Cité, CEA, CNRS, Astro-
physique, Instrumentation et Modélisation Paris-Saclay, 91191 Gif-
sur-Yvette, France

31 Université St Joseph; Faculty of Sciences, Beirut, Lebanon

32 Departamento de Física, FCFM, Universidad de Chile, Blanco Encalada 2008, Santiago, Chile

33 Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool L3 5RF, UK

34 Institute of Cosmology and Gravitation, University of Portsmouth, Portsmouth PO1 3FX, UK

35 Dipartimento di Fisica e Astronomia, Università di Bologna, Via Gobetti 93/2, 40129 Bologna, Italy

36 Dipartimento di Fisica, Università di Genova, Via Dodecaneso 33, 16146 Genova, Italy

37 INFN-Sezione di Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy

38 INAF-Osservatorio Astronomico di Capodimonte, Via Moiariello 16, 80131 Napoli, Italy

39 Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, 4150-762 Porto, Portugal

40 Institut de Física d'Altes Energies (IFAE), The Barcelona Insti-
tute of Science and Technology, Campus UAB, 08193 Bellaterra,
Barcelona, Spain

41 Port d'Informació Científica, Campus UAB, C. Albareda s/n, 08193 Bellaterra, Barcelona, Spain

42 INFN Section of Naples, Via Cinthia 6, 80126 Napoli, Italy

43 Department of Physics “E. Pancini”, University Federico II, Via Cinthia 6, 80126 Napoli, Italy

44 Centre National d'Etudes Spatiales, 18 Av. Edouard Belin, Toulouse, France

45 Institut National de Physique Nucléaire et de Physique des Particules, 3 Rue Michel-Ange, 75794 Paris Cédex 16, France

46 European Space Agency/ESRIN, Largo Galileo Galilei 1, 00044 Frascati, Roma, Italy

47 ESAC/ESA, Camino Bajo del Castillo, s/n, Urb. Villafranca del Castillo, 28692 Villanueva de la Cañada, Madrid, Spain

48 Univ. Lyon, Univ. Claude Bernard Lyon 1, CNRS/IN2P3, IP2I Lyon, UMR 5822, 69622 Villeurbanne, France

49 Institute of Physics, Laboratory of Astrophysics, Ecole Polytechnique Fédérale de Lausanne (EPFL), Observatoire de Sauverny, 1290 Versoix, Switzerland

50 Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey RH5 6NT, UK

51 Departamento de Física, Faculdade de Ciências, Universidade de Lisboa, Edifício C8, Campo Grande, 1749-016 Lisboa, Portugal

52 Instituto de Astrofísica e Ciências do Espaço, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal

53 Department of Astronomy, University of Geneva, Ch. dÉcogia 16, 1290 Versoix, Switzerland

54 Université Paris-Saclay, CNRS, Institut d'Astrophysique Spatiale, 91405 Orsay, France

55 Department of Physics, Oxford University, Keble Road, Oxford OX1 3RH, UK

56 INAF-Osservatorio Astronomico di Trieste, Via G. B. Tiepolo 11, 34143 Trieste, Italy

57 FRACTAL S.L.N.E., Calle Tulipán 2, Portal 13 1A, 28231 Las Rozas de Madrid, Spain

58 INAF-Osservatorio Astronomico di Padova, Via dell'Osservatorio 5, 35122 Padova, Italy

59 Max Planck Institute for Extraterrestrial Physics, Giessenbachstr. 1, 85748 Garching, Germany

60 Dipartimento di Fisica “Aldo Pontremoli”, Università degli Studi di Milano, Via Celoria 16, 20133 Milano, Italy

61 INAF-Osservatorio Astronomico di Brera, Via Brera 28, 20122 Milano, Italy

62 INFN-Sezione di Milano, Via Celoria 16, 20133 Milano, Italy

63 Institute of Theoretical Astrophysics, University of Oslo, PO Box 1029, Blindern 0315, Oslo, Norway

64 Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109, USA

65 von Hoerner & Sulger GmbH, SchloßPlatz 8, 68723 Schwetzingen, Germany

66 Technical University of Denmark, Elektrovej 327, 2800 Kgs. Lyngby, Denmark

67 Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany

68 Department of Physics and Helsinki Institute of Physics, University of Helsinki, Gustaf Hällströmin katu 2, 00014 Helsinki, Finland

69 European Space Agency/ESTEC, Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands

70 NOVA Optical Infrared Instrumentation Group at ASTRON, Oude Hoogeveensedijk 4, 7991 PD Dwingeloo, The Netherlands

71 Argelander-Institut für Astronomie, Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany

72 Department of Physics, Institute for Computational Cosmology, Durham University, South Road, Durham DH1 3LE, UK

73 University of Applied Sciences and Arts of Northwestern Switzerland, School of Engineering, 5210 Windisch, Switzerland

74 INFN-Bologna, Via Irnerio 46, 40126 Bologna, Italy

75 Department of Physics and Astronomy, University of Aarhus, Ny Munkegade 120, 8000 Aarhus C, Denmark

76 Centre for Astrophysics, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada

77 Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada

78 Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada

79 Space Science Data Center, Italian Space Agency, Via del Politecnico snc, 00133 Roma, Italy

80 Institute of Space Science, Bucharest 077125, Romania

81 Departamento de Astrofísica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain

82 IFPU, Institute for Fundamental Physics of the Universe, Via Beirut 2, 34151 Trieste, Italy

83 Aix-Marseille Univ., CNRS, CNES, LAM, Marseille, France

84 Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Avenida Complutense 40, 28040 Madrid, Spain

85 Instituto de Astrofísica e Ciências do Espaço, Faculdade de Ciências, Universidade de Lisboa, Tapada da Ajuda, 1349-018 Lisboa, Portugal

86 Universidad Politécnica de Cartagena, Departamento de Electrónica y Tecnología de Computadoras, 30202 Cartagena, Spain

87 Kapteyn Astronomical Institute, University of Groningen, PO Box 800, 9700 AV Groningen, The Netherlands

88 Infrared Processing and Analysis Center, California Institute of Technology, Pasadena, CA 91125, USA
## Appendix A: Fisher matrix expression proof

This paper presents a forecast for cosmological parameter measurements using the Fisher matrix technique. This appendix is devoted to the derivation of the Fisher matrix expression considering as observables either the $a_{\ell m}$ (field perspective) or the $C(\ell)$ (estimator perspective).

### A.1. The field perspective

We consider the spherical harmonics' expansion coefficients $a_{\ell m}^{Ai}$, the coefficients of the spherical harmonics decomposition of the 2D field A in the i-th tomographic bin, which we assume to follow a multivariate gaussian distribution with zero mean and covariance matrix S

$$L(a_{\ell}|\theta) = \frac{\exp\left(-\frac{1}{2}a_{\ell}^{T}S^{-1}a_{\ell}\right)}{\sqrt{(2\pi)^n|S|}}, \qquad (A.1)$$

where $\theta$ is the vector of the parameters and $a_{\ell}$ is a vector collecting the $a_{\ell m}^{Ai}$. The logarithm of this probability distribution function is

$$\ln L(a_{\ell}|\theta) = -\frac{1}{2} \left[ n \ln(2\pi) + \ln|S| + a_{\ell}^{T} S^{-1} a_{\ell} \right] \\ = -\frac{1}{2} \left[ n \ln(2\pi) + \ln|S| + \text{Tr}(S^{-1}A) \right], \qquad (A.2)$$

with $A \equiv a_{\ell} a_{\ell}^{T}$. The computation of the second derivative, after taking the expectation value, yields

$$F_{\alpha\beta}(\ell) = -\langle \ln L(a_{\ell}|\theta)_{,\alpha\beta} \rangle = \frac{1}{2} \text{Tr} (S^{-1}S_{,\alpha}S^{-1}S_{,\beta}), \qquad (A.3)$$

where $S = \langle A \rangle$ and we used the following matrix identities:

$$(\ln|S|)_{,\alpha} = \text{Tr}(S^{-1}S_{,\alpha}), \qquad (S^{-1})_{,\alpha} = -S^{-1}S_{,\alpha}S^{-1}. \qquad (A.4)$$

Now, we specify how the $a_{\ell m}^{Ai}$ were collected into the vector $a_{\ell}$. First of all, we recall that

$$\langle a_{\ell m}^{Ai} a_{\ell' m'}^{Bj} \rangle = \sum_{ij}^{AB} (\ell) \delta_{\ell \ell'} \delta_{mm'}, \qquad \sum_{ij}^{AB} (\ell) = C_{ij}^{AB}(\ell) + N_{ij}^{AB}(\ell), \qquad (A.5)$$

where $N_{ij}^{AB}$ is the Poisson shot noise for probes combination AB in tomographic bins $ij$. Working at fixed multipole $\ell$, the vector $a_{\ell}$ has a multi-index $I = (m, i, A)$. Each of these indices runs in a different range, and the range of the tomographic index $i$ may depend on the probe A considered. To summarize:

- $m$ varies from $-\ell$ to $\ell$

- $i$ varies from 1 to $N_A$

- A varies from 1 to $N$,

with $N$ being the number of probes and $N_A$ the number of tomographic bins for probe A. We chose to order the array $a_{\ell}$ varying the three indices with a significance increasing from left to right. With this choice the matrix S has the following block form:

$$S = \Sigma(\ell) \otimes \mathbb{I}_{2\ell+1}, \qquad \Sigma(\ell) = \begin{pmatrix} \Sigma^{11}(\ell) & \Sigma^{12}(\ell) & \cdots & \Sigma^{1N}(\ell) \\ \Sigma^{21}(\ell) & \Sigma^{22}(\ell) & \cdots & \Sigma^{2N}(\ell) \\ \vdots & \vdots & \ddots & \vdots \\ \Sigma^{N1}(\ell) & \Sigma^{N2}(\ell) & \cdots & \Sigma^{NN}(\ell) \end{pmatrix}. \qquad (A.6)$$

Here $\Sigma^{AB}(\ell)$ is the tomographic covariance matrix, that is to say, $[\Sigma^{AB}(\ell)]_{ij} = \Sigma_{ij}^{AB}(\ell)$, and $\otimes$ is the Kronecker product.

Equation (A.6) can be understood as follows: the matrix S is made up of diagonal blocks of size $\ell$, each of them proportional to the identity $\mathbb{I}_{2\ell+1}$ with a different factor $\Sigma_{ij}^{AB}(\ell)$. Here the identity matrix $\mathbb{I}_{2\ell+1}$ is exactly the Kronecker delta $\delta_{mm'}$, which appears in Eq. (A.5). In order to compute the trace in Eq. (A.3) we make use of some properties of the Kronecker product (Steeb 1997):

$$(A \otimes B)(C \otimes D) = AC \otimes BD, \qquad (A.7)$$

$$(A \otimes B)^{-1} = A^{-1} \otimes B^{-1}, \qquad (A.8)$$

$$\text{Tr}(A \otimes B) = \text{Tr}(A) \text{Tr}(B). \qquad (A.9)$$

Using these properties, we have

$$\begin{aligned} \text{Tr} [S^{-1}S_{,\alpha}S^{-1}S_{,\beta}] &= \text{Tr} [\Sigma(\ell)^{-1}\Sigma(\ell)_{,\alpha}\Sigma(\ell)^{-1}\Sigma(\ell)_{,\beta}] \otimes \mathbb{I}_{2\ell+1}] \\ &= \text{Tr} [\Sigma(\ell)^{-1}\Sigma(\ell)_{,\alpha}\Sigma(\ell)^{-1}\Sigma(\ell)_{,\beta}] \text{Tr}(\mathbb{I}_{2\ell+1}) \\ &= (2\ell + 1) \text{Tr} [\Sigma(\ell)^{-1}\Sigma(\ell)_{,\alpha}\Sigma(\ell)^{-1}\Sigma(\ell)_{,\beta}] \\ &= (2\ell + 1) \text{Tr} [\Sigma(\ell)^{-1}C(\ell)_{,\alpha}\Sigma(\ell)^{-1}C(\ell)_{,\beta}], \end{aligned}$$

where in the last line we used the fact that the shot noise is independent of the cosmological parameters, and the $C(\ell)$'s block matrix $C(\ell)$ has the same form of $\Sigma(\ell)$ of Eq. (A.6), with entries $C_{ij}^{AB}(\ell)$ instead of $\Sigma(\ell)_{ij}^{AB}$. Then from Eq. (A.3), we can finally write

$$F_{\alpha\beta}(\ell) = \frac{2\ell + 1}{2} \text{Tr} \left( \Sigma(\ell)^{-1} C(\ell)_{,\alpha} \Sigma(\ell)^{-1} C(\ell)_{,\beta} \right).$$

Now, we sum over the multipoles $\ell$, considering them as independent

$$F_{\alpha\beta} = \sum_{\ell} \frac{2\ell + 1}{2} \text{Tr} [\Sigma(\ell)^{-1}C(\ell)_{,\alpha}\Sigma(\ell)^{-1}C(\ell)_{,\beta}]. \qquad (A.10)$$

In order to obtain the expression Eq. (43), it sufficient to redefine the covariance as

$$\Sigma(\ell) \rightarrow \sqrt{\frac{2}{(2\ell + 1)\Delta\ell f_{\text{sky}}}} \Sigma(\ell). \qquad (A.11)$$

In this way one accounts for a possible unequal spacing between multipole bins, weighting each term of the sum with the bin width, and with the sky fraction $f_{\text{sky}}$ covered by the survey.

### A.2. The estimator perspective

In the previous section, we evaluated the expression of the Fisher matrix for a multivariate normal distribution; then, we specialized this expression for the field perspective, when the observables are the $a_{\ell m}$'s. In this section we derive the Fisher matrix expression when the observables are the estimator $\hat{C}(\ell)$, defined as

$$\hat{C}_{ij}^{AB}(\ell) = \frac{1}{2\ell + 1} \sum_{m=-\ell}^{\ell} a_{\ell m}^{Ai} a_{\ell m}^{Bj}. \qquad (A.12)$$

One may argue that, since the $a_{\ell m}$ are normally distributed, the $\hat{C}(\ell)$ are also normally distributed. However, this ansatz is wrong; in the proceeding equation, we obtain the $\hat{C}(\ell)$ posterior distribution function and then their Fisher matrix. Since coefficients with different $\ell$ are uncorrelated, under the Gaussian assumption, we work at a fixed $\ell$. The $\hat{C}(\ell)$ likelihood is

$$L(\hat{C}|\theta) = \int d\mathbf{a}_{\ell} L(\mathbf{a}_{\ell}|\theta) \prod_{A,B} \prod_{i,j} \delta_D \left( \hat{C}_{ij}^{AB}(\ell) - \sum_{m=-\ell}^{\ell} \frac{a_{\ell m}^{Ai} a_{\ell m}^{Bj}}{2\ell + 1} \right). \qquad (A.13)$$
Using the Dirac delta function $\delta_D$, $L(a_t|\theta)$ gets out from the integral

$$L(\hat{C}|\theta) = \frac{f(\hat{C})}{\sqrt{(2\pi)^n|S|}} \exp\left(-\frac{1}{2}\text{Tr}(S^{-1}A)\right), \qquad (A.14)$$

where $f(\hat{C})$ is the result of the integral in the $a_{\ell m}$ space<fnref n="9" />. Taking the logarithm of the previous equation we obtain

$$\ln L(\hat{C}|\theta) = \ln\left(-\frac{f(\hat{C})}{\sqrt{(2\pi)^n}}\right) - \frac{1}{2}\ln(|S|) - \frac{1}{2}\text{Tr}(S^{-1}A). \qquad (A.15)$$

Proceeding as in the previous section and using $\langle A \rangle = S$ the expectation value of the second derivative becomes

$$\langle \ln L(\hat{C}|\theta, \alpha\beta) \rangle = -\frac{1}{2} \text{Tr}(S^{-1}S_{,\beta}S^{-1}S_{,\alpha}). \qquad (A.16)$$

Performing the same manipulations of the previous section, we finally obtain Eq. (A.10), showing that the Fisher matrix in the field and estimator perspective is the same.

<footnote n="9">⁹ Since it does not carry information on the parameters $\theta$, we do not derive it explicitly. If computing the integral, one would have found the Wishart probability distribution function.</footnote>