Issue |
A&A
Volume 670, February 2023
|
|
---|---|---|
Article Number | A47 | |
Number of page(s) | 20 | |
Section | Cosmology (including clusters of galaxies) | |
DOI | https://doi.org/10.1051/0004-6361/202244445 | |
Published online | 02 February 2023 |
Euclid: Forecasts from the void-lensing cross-correlation⋆
1
Dipartimento di Fisica, Universitá degli Studi di Genova, and INFN-Sezione di Genova, Via Dodecaneso 33, 16146 Genova, Italy
2
INAF-IASF Milano, Via Alfonso Corti 12, 20133 Milano, Italy
e-mail: marco.bonici@inaf.it
3
INFN-Sezione di Genova, Via Dodecaneso 33, 16146 Genova, Italy
4
Aix-Marseille Univ., CNRS/IN2P3, CPPM, Marseille, France
5
INAF-Osservatorio Astronomico di Roma, Via Frascati 33, 00078 Monteporzio Catone, Italy
6
INFN-Sezione di Roma, Piazzale Aldo Moro, 2 – c/o Dipartimento di Fisica, Edificio G. Marconi, 00185 Roma, Italy
7
Universitäts-Sternwarte München, Fakultät für Physik, Ludwig-Maximilians-Universität München, Scheinerstrasse 1, 81679 München, Germany
8
Department of Astrophysical Sciences, Peyton Hall, Princeton University, Princeton, NJ 08544, USA
9
Instituto de Astrofísica de Canarias (IAC); Departamento de Astrofísica, Universidad de La Laguna (ULL), 38200 La Laguna, Tenerife, Spain
10
Instituto de Astrofísica de Canarias, Calle Vía Láctea s/n, 38204 San Cristóbal de La Laguna, Tenerife, Spain
11
Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK
12
Dipartimento di Fisica e Astronomia “Augusto Righi” – Alma Mater Studiorum Università di Bologna, Via Piero Gobetti 93/2, 40129 Bologna, Italy
13
INAF-Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, Via Piero Gobetti 93/3, 40129 Bologna, Italy
14
INFN-Sezione di Bologna, Viale Berti Pichat 6/2, 40127 Bologna, Italy
15
Dipartimento di Fisica e Astronomia “G. Galilei”, Universitá di Padova, Via Marzolo 8, 35131 Padova, Italy
16
INFN-Padova, Via Marzolo 8, 35131 Padova, Italy
17
Université de Genève, Département de Physique Théorique and Centre for Astroparticle Physics, 24 Quai Ernest-Ansermet, 1211 Genève 4, Switzerland
18
Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Can Magrans, s/n, 08193 Barcelona, Spain
19
Institut d’Estudis Espacials de Catalunya (IEEC), Carrer Gran Capitá 2-4, 08034 Barcelona, Spain
20
Institut de Recherche en Astrophysique et Planétologie (IRAP), Université de Toulouse, CNRS, UPS, CNES, 14 Av. Edouard Belin, 31400 Toulouse, France
21
University of Lyon, UCB Lyon 1, CNRS/IN2P3, IUF, IP2I Lyon, France
22
Istituto Nazionale di Astrofisica (INAF) – Osservatorio di Astrofisica e Scienza dello Spazio (OAS), Via Gobetti 93/3, 40127 Bologna, Italy
23
Istituto Nazionale di Fisica Nucleare, Sezione di Bologna, Via Irnerio 46, 40126 Bologna, Italy
24
Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK
25
Dipartimento di Fisica, Universitá degli Studi di Torino, Via P. Giuria 1, 10125 Torino, Italy
26
INFN-Sezione di Torino, Via P. Giuria 1, 10125 Torino, Italy
27
INAF-Osservatorio Astrofisico di Torino, Via Osservatorio 20, 10025 Pino Torinese, TO, Italy
28
Instituto de Física Teórica UAM-CSIC, Campus de Cantoblanco, 28049 Madrid, Spain
29
AIM, CEA, CNRS, Université Paris-Saclay, Université de Paris, 91191 Gif-sur-Yvette, France
30
Université Paris-Saclay, Université Paris Cité, CEA, CNRS, Astrophysique, 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 Institute 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
Received:
7
July
2022
Accepted:
29
October
2022
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 Ωb are competitive with galaxy lensing alone, while errors on ns 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, Ωb, the neutrino mass, and the DE evolution.
Key words: gravitational lensing: weak / cosmological parameters / large-scale structure of Universe
© The Authors 2023
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe-to-Open model. Subscribe to A&A to support open access publication.
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) = w0 + z/(1 + z) wa (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 (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).
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 ng = 30 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 (GCph): 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, νΛCDM, w0waCDM, and νw0waCDM model cosmologies, adopting the following set of cosmological parameters: the reduced Hubble constant h defined via H0 = h 100 km s−1 Mpc−1, the baryon density parameter Ωb at present time, the total matter density parameter Ωm at present time, the sum of the three active neutrino masses Mν = ∑imνi, the massive neutrino density parameter Ων = Mν (eV)/(93.14 h2) at present time, the cold DM density parameter Ωc = Ωm − Ωb − Ων at present time, where Ωm = Ωb + Ωc + Ων, the DE density parameter ΩDE = 1 − Ωm at present time, the parameters of the DE equation of state, w0 and wa, the rms of the matter linear density fluctuations, σ8, inside a radius of 8 h−1 Mpc, and the scalar spectral index ns. The cosmological parameters varied in the analysis, together with their fiducial values, are summarized in Table 1.
Cosmological parameters varied in the forecast analysis, together with their values in the fiducial cosmology assumed in this work.
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, , of cosmic voids from the void auto-power spectrum, Pvv(k, z), while the angular cross-power spectra,
, between voids and WL were evaluated from the void-matter cross-spectrum Pvm(k, z). The latter are shown, together with the nonlinear matter power spectrum Pmm(k, z), in the top left panel of Fig. 1. Pvv(k, z) was obtained from the matter power spectrum, Pmm(k, z), assuming the Poisson shot noise (Hamaus et al. 2014a; Chan et al. 2014, 2019; Jamieson & Loverde 2019) and the “effective void bias”
defined in Eq. (19). To this aim we adopted the following relations:
![]() |
Fig. 1. Figure containing the power spectra used in this work. Top left: nonlinear matter auto-power spectrum Pmm (solid red line), the void auto-power spectrum Pvv without the shot noise (dashed blue line), computed assuming the pessimistic void bias scenario with bv eff ≈ −11.9, and the absolute value of the void-matter cross-power spectrum Pvm (dot-dashed black line), all computed at redshift z = 0.001. Top right: tomographic void angular auto-power spectra |
where the void bias, , is assumed to be scale-independent1,
is the void number density, and
is the nonlinear matter power spectrum when an additional wavenumber filter, SN(k), is applied to cut small scales. The low-k pass filter, SN(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, Pvv(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 Pvm(k, z) is defined as
To avoid numerical instabilities, we prefer the low-k pass filter rather than a sharp cut in k. SN(k) is the so-called smoothstep function, defined as
where N measures the degree of smoothness of the function itself: the first discontinuous derivative of SN(k) is the (N + 1)th derivative. The order used in SN for producing Pvv(k, z) and Pvm(k, z) is N = 3. The smoothstep parameters are set to: a = 1, b = 1.8, and kref = 0.25 h Mpc−1; the filter suppresses Pvv(k, z) and Pvm(k, z) for k > kref. The cutoff scale kref was determined combining the mean radius of voids in the catalog with the void exclusion principle, from which we know that the cutoff scale is given by
(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 as
where γ stands for cosmic shear.
As we describe below, to obtain the effective void bias, , we used the peak-background split (PBS) formalism (Sheth & Van De Weygaert 2004), and weighted the void bias, bv(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, δc, linearly extrapolated to the present time as δc(z) = δc(0)/G(z), where G(z) is the so-called linear growth factor normalized at z = 0, and δc(0) slightly depends on cosmology, here assumed to be δc(0) = 1.686. The latter is the density threshold for the void formation δv, linearly extrapolated to the present time as δv(z) = δv(0)/G(z), with δ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 δv(0) can be used to model voids more reliably). We chose a value of δ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 δv and δc have an important role. The void formation threshold, δv, is the value an under-dense region needs to overcome to turn into a void. The density threshold, δ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 δc and δv, via the so-called void-in-cloud parameter, 𝒟, and the dimensionless parameter, x, defined as
where σ2 is the variance of the filtered linear density field on a scale R
Here Plin(k) is the linear matter power spectrum at z = 0, and W(kR) is the top-hat filter in Fourier space:
The abundance of voids with mass M is predicted by the SvdW model as
where ρm is the mean background density of matter and the multiplicity function, flnσ, represents the fraction of voids in a unit range of lnM, and it is shown to be (Sheth & Van De Weygaert 2004)
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
where is the volume of a sphere with radius rL 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:
where the subscript L denotes quantities computed in the linear regime. The V dn void size function then reads (Jennings et al. 2013)
where rv indicates the actual void radius. The conversion from rL to rv is given by
where ρv is the average density within the void.
However, the expression of the multiplicity function, flnσ, in Eq. (9) is difficult to evaluate in practice. Therefore, Jennings et al. (2013) suggested the following approximation that works at 0.2% accuracy:
Finally, the V dn size function is given by
and, after some rearrangements, its final expression is
where rv = 1.7 rL is predicted by considering a spherical expansion.
As mentioned previously, the void bias, bv, at a fixed void radius, rv, is computed using the PBS formalism (Sheth & Van De Weygaert 2004; Chan et al. 2014):
where the z dependence is encapsulated in δc and δv, and the rv dependence in ν:
We defined the effective void bias by weighting Eq. (17) with the void size function V dn in Eq. (16):
exploiting in this way all void radii r ≥ rv, minh−1 Mpc. In this work, we assumed rv, 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, is evolved with the growth factor of the fiducial cosmology. Its fiducial value today is assumed to be
, 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 δc and δ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 bv, entering
, 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
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
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(ℓ), which are functions of the multipole number ℓ, and are defined as the spherical harmonic transform of a two-point correlation function. The aℓm, the coefficients of the spherical harmonics decomposition of a field ℱ on the sphere (e.g., the cosmic shear or void density fields), are defined on the sky as
The C(ℓ) represent the two-point correlation function of the aℓm and they are diagonal in ℓ and independent of m because of statistical isotropy.
In the forecasts presented in this work, we used three kinds of C(ℓ): the void-void auto-correlation Cvv(ℓ), the lensing-lensing auto-correlation Cγγ(ℓ), and the void-lensing cross-correlation Cvγ(ℓ). We computed the C(ℓ) “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 . For example,
is the spherical harmonic transform of the correlation function between the lensing signal in the ith bin and the void signal in the jth bin.
We evaluated Cvv(ℓ) and Cvγ(ℓ) for ℓ ∈ [10, 1500], while for Cγγ(ℓ), two scenarios were considered: a “pessimistic” WL scenario, with ℓ ∈ [10, 1500], and an “optimistic” WL scenario, with ℓ ∈ [10, 5000].
The Limber approximation (Limber 1953; LoVerde & Afshordi 2008) allowed us to write a simple integral expression for the tomographic C(ℓ), where the power spectrum , defined in Eq. (4), enters the integral as
where are suitable “weight” functions for the different probes A and B defined in Eqs. (27) and (30), E(z) is the dimensionless Hubble parameter,
and r(z) is the comoving distance,
Throughout the forecast, the integral expression of was computed numerically with the Simpson method (Brew & Atkinson 1979). The signals
,
, and
, 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 zmin = 0.001 and zmax = 2.5, respectively. The C(ℓ) and the distributions of galaxies and voids are sampled in ten redshift bins i, with the following endpoints:
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
with i, j being the tomographic indices, and A, B the considered probes.
The tomographic void-void angular power spectra were computed as
where is defined in Eq. (1) and we omitted the shot-noise additive term
defined in Eq. (40).
The void weight function is defined as
The void effective bias, bveff, i, entering in the void weight function, is given by bveff(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 catalog2, as explained in Sect. 3.
The tomographic void-lensing angular power spectra, , and the lensing-lensing angular power spectra,
, were computed, respectively, as
and
where in 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
with
Here is the lensing efficiency, defined as an integral expression of the observed galaxy distribution
.
The density distribution of the “observed” galaxies in the ith tomographic bin, , was computed as a convolution of the galaxy distribution, ng(z), and the photometric instrument response, as described below. According to EC20, the galaxy redshift distribution, ng(z), adopted in this work is
where (Laureijs et al. 2011). Redshift measurements are affected by experimental and reconstruction effects. We modeled the probability pph(zp|z) that a true galaxy redshift z will be reconstructed with a photometric redshift zp as
The modeling and the parameters entering pph(zp|z) are reported in Table 2, and were taken from EC20. The convolution of ng(z) and pph(zp|z) estimates the galaxy redshift distribution as measured by NISP in the photometric mode, together with ground-based observation. The density distribution, , of the “observed” galaxies in the ith tomographic bin is obtained convolving pph(zp|z) with ng(z) in the ith bin
where and
are the edges of the ith bin as defined in Eq. (24). The density distributions,
, in the different tomographic bins are shown in the left panel of Fig. 2.
![]() |
Fig. 2. Galaxy and void distributions used in this work. Left: dimensionless normalized projected galaxy distributions |
The matter power spectrum, Pmm(k, z), was computed numerically using as Boltzmann solver the “Code for Anisotropies in the Microwave Background”3 (CAMB; Lewis et al. 1999)4. In order to obtain the cosmological parameter forecasts presented in this work, we computed the matter power spectra, Pmm(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 Pmm(k, z) over a k − z grid. According to EC20, the range of redshifts chosen in the present analysis is z ∈ [0.001, 2.5], the range of wave numbers is k ∈ [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 and
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 C(ℓ) derivatives, and enters the matter power spectra, Pmm(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
.
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 6003 DM particles in a periodic box of 3780 h−1 Mpc on a side (Potter et al. 2017), with a flat ΛCDM cosmology very similar to that reported in Table 1, namely Ωm = 0.319, Ωb = 0.049, ΩΛ = 0.681, σ8 = 0.83, and ns = 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 Δ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 rv, min ∼ 25 h−1 Mpc up to rv, max ∼ 300 h−1 Mpc, and its projected spatial density, nv(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 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 Δz = 0.2. Then, from Eq. (1), we inferred the void bias as
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 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.
![]() |
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. |
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:
where α and β are the parameters of interest. The expected “error covariance matrix” is the inverse of the Fisher matrix:
The diagonal elements of the error covariance matrix are the squares of the marginalized 1-σ errors on the parameters:
As we deal with angular power spectra, it is convenient to introduce the matrix associated with a given C(ℓ):
where Δℓ is the multipole bin width, and fsky the sky fraction covered by the survey. The “shot-noise matrix”, , depends on the particular probe combination. For the C(ℓ) of the void clustering, WL, and void-lensing, assuming Poisson statistics, they are respectively
where δij is the Kronecker delta, is the galaxy shape noise, and
and
are the un-normalized average galaxy and void surface densities in the ith 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ℓm is simply given by Eq. (39). When two or more probes are combined together, one needs to construct the full covariance matrix, ΣXC, composed of matrix blocks defined in Eq. (39). For the probes γ and v considered in this work, ΣXC is given by
The covariance and the C(ℓ) matrix share the same structure, and therefore the block matrix reads
Assuming that the aℓm, defined in Eq. (20), follow a multivariate Gaussian distribution, an analytical expression for the Fisher matrix elements is given by5
The formula above applies both to the single- and two-probe correlation case, provided that the C(ℓ) and the covariance matrices are chosen accordingly. The Fisher matrix computation involves the derivatives of the C(ℓ) with respect to cosmological and nuisance parameters. In order to ensure reliable results, numerical derivatives were computed with6 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-σ errors, together with the DE figure of merit7 (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 w0, wa, Mν, evaluated respectively for the pessimistic and optimistic scenarios, are also shown in the top panels of Fig. 9.
![]() |
Fig. 4. Fisher matrix marginalized contours for the (baseline) ΛCDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios. |
![]() |
Fig. 5. Fisher matrix marginalized contours for the νΛCDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios. |
![]() |
Fig. 6. Fisher matrix marginalized contours for the w0waCDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios. |
![]() |
Fig. 7. Full contour plots when GCsp is combined with WL + V + XC for the w0waCDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios. |
![]() |
Fig. 8. Fisher matrix marginalized contours for the νw0waCDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios. |
Probe configurations in the pessimistic and optimistic scenarios for the combination V + WL, both considered for two values of rv, min.
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 Ωb in the optimistic case, where for h void clustering provides better constraints than WL by a factor of ∼3, while for Ωb the constraints from the two single probes look comparable. This is due to the form of the WL kernel and the integration along the line of sight, and to the presence of baryon acoustic oscillations (BAO) in the void angular auto power spectrum, Cvv(ℓ), and the void-lensing cross spectrum, Cvγ(ℓ), (see Fig. 1), which are instead completely washed out in the WL angular power spectrum, Cγγ(ℓ). 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 Ωb, but also ns; similar conclusions (without WL) were found by Kreisch et al. (2021).
In the case of ns, this happens since the WL and void clustering ellipses in the ns − Ωb plane happen to be orthogonal. We verified that, even if ns and Ωb produce an increment in the same direction for Cγγ(ℓ) and Cvv(ℓ), which would imply that the two parameters are negatively correlated, this orthogonality comes from projecting the large parameter space onto the ns − Ω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 ns decreases by a factor of ∼2, the constraints on σ8 and Ωm get ∼40% tighter, and the uncertainty on h decreases by a factor of ∼3. The increase in the constraining power on σ8 is somehow expected: the void bias affects the overall amplitude of the C(ℓ), so it is degenerate with σ8. On the one hand, in the pessimistic scenario, both the value of the bias parameters and σ8 were assumed to be measured from data, and this increases the uncertainties on σ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, GCsp, 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 σ8, weakening its constraints by a factor of ∼2 with respect to the baseline ΛCDM case. This is expected since the effects of σ8 and massive neutrinos on the matter power spectrum are similar: the former regulates the normalization of Pmm, and the latter suppresses Pmm 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 eV8, 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 w0waCDM 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 w0 and wa mainly affects h, σ8, and Ωm, increasing their errors. The impact on h and Ωm is explained considering that w0 and wa enter the Hubble parameter H(z) (22), where the dependence from the DE equation of state has been exploited. Moreover, w0 and wa also enter the linear growth factor (Linder & Jenkins 2003), and therefore impact Ωm and the normalization σ8 again. When moving from the pessimistic to the optimistic scenario, the constraints on w0 and wa 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, bv, 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.
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.
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, GCsp, to the WL, void clustering, and void-lensing cross-correlation probes. To this purpose, we considered GCsp as a probe independent of the others, and made use of the GCsp 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 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.
Marginalized 1-σ errors in different cosmological scenarios for WL + V + GCsp, when they are assumed to be independent, together with their combinations when the void-lensing cross-correlation is included.
In the pessimistic scenario, constraints given by GCsp + V are weaker than the ones given by GCsp + 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 GCsp + WL + V are only slightly improved with respect to GCsp + WL. Nevertheless, the constraints given by GCsp + WL + V + XC are tighter than in the GCsp + WL case. In particular, the FoM increases by ∼5% when considering GCsp + WL + V + XC against GCsp + WL.
In the optimistic scenario, constraints given by GCsp + V are comparable to the ones given by GCsp + WL. Moreover, the FoM increases from 192 for GCsp + WL, up to 702 for GCsp + WL + V, in which case constraints on the other cosmological parameters become tighter: by ∼10% for h, by ∼30% for ns, by ∼50 − 55% for Ωb and σ8, and by ∼70% for Ωm.
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 Ωm and σ8, with their constraints improved by a factor of 2 and 3, respectively; the total FoM increases by ∼10%, that is, from 702 for GCsp + WL + V, up to 791 for GCsp + WL + V + XC.
5.5. Combining massive neutrinos with dynamical dark energy
In this section we consider Mν, w0, and wa 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 σ8 gets even weaker with respect to the baseline ΛCDM case, due to the combined variation of Mν, w0, and wa, entering the matter power spectrum, the Hubble parameter, and the growth factor. Comparing 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 w0 and wa increases the error on Mν and vice versa. The increase in the w0 and wa 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 w0 and wa. However, especially for w0 and wa, 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.
![]() |
Fig. 9. Forecast contours for the νw0waCDM model, in the pessimistic (left) and in the optimistic (right) scenarios. |
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-σ 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-σ uncertainty was about 10%.
The impact of the k-binning was negligible too. CAMB evaluates the matter power spectrum Pmm(k, z) with a logarithmic k-binning grid. For the reported forecasts, Pmm(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-σ errors remained unchanged up to the third digit.
The impact of the z-binning was minor. CAMB evaluates the matter power spectrum Pmm(k, z) with a linear z-binning grid. For the reported forecast, Pmm(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-σ 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-σ 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.
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, nv(z), directly from the Euclid Flagship photometric mock galaxy catalog, using the 2D void finder of Sánchez et al. (2017);
-
we evaluated the C(ℓ) 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 V dn void size function model (Jennings et al. 2013).
We present our results for two different choices of the effective void bias : 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
was marginalized over as a nuisance parameter.
We present parameter forecasts for different cosmological models: starting from a flat Λ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:
-
Λ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 Ω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 Ω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 ns, are the parameters for which the constraints improve most when void clustering is combined with WL.
-
νΛCDM cosmology: adding to the analysis the total neutrino mass, Mν, as a free parameter mainly impacts the constraints on σ8, as expected, since σ8 and Mν both affect the amplitude of the matter power spectrum (even if in the Mν-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ν becomes stronger by ∼5% and ∼15%, respectively, in the pessimistic and optimistic scenarios.
-
w0waCDM cosmology: when the CPL parametrization of the DE equation of state is considered, we find that the most affected parameters are h, Ωm, and σ8. The expression of the Hubble parameter Eq. (22) explains the effect on h and Ωm, while the impact of w0 and wa on the linear growth factor (Linder & Jenkins 2003) accounts for the effect on σ8 and, again, Ω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.
-
νw0waCDM cosmology: in this model we evaluated, in particular, the mutual impact among w0, wa, and Mν. On the one hand, we find that, when adding Mν, the marginalized errors of w0 and wa increase at most by ∼15%. On the other hand, when adding w0, wa as free parameters, the impact on the Mν constraints is higher, as its error increases even by 70%. However, especially for w0 and wa, 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 GCsp: finally, we combined the galaxy lensing, void clustering, and void-lensing cross-correlation probes with the spectroscopic galaxy clustering, GCsp. In this case, we exploited the Fisher matrices as obtained by EC20 in the so-called “pessimistic” and “optimistic” settings, and considered the flat w0waCDM scenario, since GCsp forecasts are computed keeping Mν = 0.06 eV fixed. When combining WL + V + XC with GCsp, 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 + GCsp, 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, GCph, 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 GCph, 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 GCph since, as shown in this work, the combination of void-clustering, void-lensing, and GCsp 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 mispecifications 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).
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.
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).
In Appendix A we show that this expression is valid both in the so-called “field” and “estimator” perspectives.
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.
Acknowledgments
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 Incorporació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 Onderzoekschool 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 [Google Scholar]
- Aubert, M., Cousinou, M.-C., Escoffier, S., et al. 2020, MNRAS, 513, 186 [Google Scholar]
- Behroozi, P. S., Wechsler, R. H., Wu, H.-Y., et al. 2013, ApJ, 762, 109 [NASA ADS] [CrossRef] [Google Scholar]
- Bird, S., Viel, M., & Haehnelt, M. G. 2011, MNRAS, 420, 2551 [Google Scholar]
- Blanton, M. R., Hogg, D. W., Bahcall, N. A., et al. 2003, ApJ, 592, 819 [Google Scholar]
- Blanton, M. R., Lupton, R. H., Schlegel, D. J., et al. 2005, ApJ, 631, 208 [Google Scholar]
- Bos, E. G., van de Weygaert, R., Dolag, K., & Pettorino, V. 2012, MNRAS, 426, 440 [NASA ADS] [CrossRef] [Google Scholar]
- Brew, E., & Atkinson, K. E. 1979, The Mathematical Gazette, 63 [Google Scholar]
- Brinckmann, T., & Lesgourgues, J. 2019, Phys. Dark Univ., 24, 100260 [Google Scholar]
- Camera, S., Harrison, I., Bonaldi, A., & Brown, M. L. 2017, MNRAS, 464, 4747 [NASA ADS] [CrossRef] [Google Scholar]
- Carretero, J., Castander, F. J., Gaztañaga, E., Crocce, M., & Fosalba, P. 2015, MNRAS, 447, 646 [NASA ADS] [CrossRef] [Google Scholar]
- Chan, K. C., Hamaus, N., & Desjacques, V. 2014, Phys. Rev. D, 90, 103521 [NASA ADS] [CrossRef] [Google Scholar]
- Chan, K. C., Hamaus, N., & Biagetti, M. 2019, Phys. Rev. D, 99, 121304 [CrossRef] [Google Scholar]
- Chevallier, M., & Polarski, D. 2001, Int. J. Mod. Phys., D, 10, 213 [NASA ADS] [CrossRef] [Google Scholar]
- Clampitt, J., & Jain, B. 2015, MNRAS, 454, 3357 [NASA ADS] [CrossRef] [Google Scholar]
- Contarini, S., Ronconi, T., Marulli, F., et al. 2019, MNRAS, 488, 3526 [NASA ADS] [CrossRef] [Google Scholar]
- Contarini, S., Verza, G., Pisani, A., et al. 2022, A&A, 667, A162 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cooray, A., & Sheth, R. 2002, Phys. Rep., 372, 1 [Google Scholar]
- Copeland, E. J., Sami, M., & Tsujikawa, S. 2006, Int. J. Mod. Phys. D, 15, 1753 [NASA ADS] [CrossRef] [Google Scholar]
- Costille, A., Caillat, A., Rossin, C., et al. 2019, SPIE Conf. Ser., 11180, 1118013 [NASA ADS] [Google Scholar]
- Crocce, M., Castander, F. J., Gaztañaga, E., Fosalba, P., & Carretero, J. 2015, MNRAS, 453, 1513 [NASA ADS] [CrossRef] [Google Scholar]
- Cropper, M., Pottinger, S., Niemi, S., et al. 2016, SPIE Conf. Ser., 9904, 99040Q [NASA ADS] [Google Scholar]
- Euclid Collaboration (Blanchard, A., et al.) 2020, A&A, 642, A191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fang, Y., Hamaus, N., Jain, B., et al. 2019, MNRAS, 490, 3573 [NASA ADS] [CrossRef] [Google Scholar]
- Fisher, R. A. 1935, J. R. Stat. Soc., 98, 39 [CrossRef] [Google Scholar]
- Granett, B. R., Neyrinck, M. C., & Szapudi, I. 2008, ApJ, 683, L99 [NASA ADS] [CrossRef] [Google Scholar]
- Granett, B. R., Kovács, A., & Hawken, A. J. 2015, MNRAS, 454, 2804 [NASA ADS] [CrossRef] [Google Scholar]
- Hamaus, N., Wandelt, B. D., Sutter, P. M., Lavaux, G., & Warren, M. S. 2014a, Phys. Rev. Lett., 112, 041304 [NASA ADS] [CrossRef] [Google Scholar]
- Hamaus, N., Sutter, P. M., & Wandelt, B. D. 2014b, Phys. Rev. Lett., 112, 251302 [NASA ADS] [CrossRef] [Google Scholar]
- Hamaus, N., Pisani, A., Sutter, P. M., et al. 2016, Phys. Rev. Lett., 117, 091302 [NASA ADS] [CrossRef] [Google Scholar]
- Hamaus, N., Cousinou, M.-C., Pisani, A., et al. 2017, JCAP, 07, 014 [CrossRef] [Google Scholar]
- Hamaus, N., Pisani, A., Choi, J.-A., et al. 2020, JCAP, 12, 023 [CrossRef] [Google Scholar]
- Hamaus, N., Aubert, M., Pisani, A., et al. 2021, A&A, 658, A20 [Google Scholar]
- Hand, N., Feng, Y., Beutler, F., et al. 2018, AJ, 156, 160 [Google Scholar]
- Hawken, A. J., Aubert, M., Pisani, A., et al. 2020, JCAP, 06, 012 [CrossRef] [Google Scholar]
- Hu, W. 1999, ApJ, 522, L21 [Google Scholar]
- Ilić, S., Langer, M., & Douspis, M. 2013, A&A, 556, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jamieson, D., & Loverde, M. 2019, Phys. Rev. D, 100, 123528 [Google Scholar]
- Jennings, E., Li, Y., & Hu, W. 2013, MNRAS, 434, 2167 [CrossRef] [Google Scholar]
- Jimenez, R., Kitching, T., Pena-Garay, C., & Verde, L. 2010, JCAP, 05, 035 [CrossRef] [Google Scholar]
- Kiakotou, A., Elgarøy, Ø., & Lahav, O. 2008, Phys. Rev. D, 77, 063005 [Google Scholar]
- Krause, E., Chang, T.-C., Doré, O., & Umetsu, K. 2013, ApJ, 762, L20 [NASA ADS] [CrossRef] [Google Scholar]
- Kreisch, C. D., Pisani, A., Carbone, C., et al. 2019, MNRAS, 488, 4413 [NASA ADS] [CrossRef] [Google Scholar]
- Kreisch, C. D., Pisani, A., Villaescusa-Navarro, F., et al. 2021, ApJ, 935, 100 [Google Scholar]
- Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, ArXiv e-prints [arXiv:1110.3193] [Google Scholar]
- Lavaux, G., & Wandelt, B. D. 2010, MNRAS, 403, 1392 [NASA ADS] [CrossRef] [Google Scholar]
- Lavaux, G., & Wandelt, B. D. 2012, ApJ, 754, 109 [NASA ADS] [CrossRef] [Google Scholar]
- Leclercq, F., Jasche, J., Sutter, P. M., Hamaus, N., & Wandelt, B. 2015, JCAP, 03, 047 [CrossRef] [Google Scholar]
- Lee, J., & Park, D. 2009, ApJ, 696, L10 [NASA ADS] [CrossRef] [Google Scholar]
- Lesgourgues, J., & Pastor, S. 2006, Phys. Rep., 429, 307 [Google Scholar]
- Lesgourgues, J., & Pastor, S. 2012, Adv. High Energy Phys., 2012, 608515 [CrossRef] [Google Scholar]
- Lesgourgues, J., & Pastor, S. 2014, New J. Phys., 16, 065002 [NASA ADS] [CrossRef] [Google Scholar]
- Lewis, A., Challinor, A., & Lasenby, A. 1999, ApJ, 538, 473 [Google Scholar]
- Limber, D. N. 1953, ApJ, 117, 134 [NASA ADS] [CrossRef] [Google Scholar]
- Linder, E. V. 2002, Phys. Rev. Lett., 90, 4 [Google Scholar]
- Linder, E. V., & Jenkins, A. 2003, MNRAS, 346, 573 [Google Scholar]
- LoVerde, M., & Afshordi, N. 2008, Phys. Rev. D, 78, 123506 [NASA ADS] [CrossRef] [Google Scholar]
- Mao, Q., Berlind, A. A., Scherrer, R. J., et al. 2017, ApJ, 835, 161 [NASA ADS] [CrossRef] [Google Scholar]
- Massara, E., Villaescusa-Navarro, F., Viel, M., & Sutter, P. M. 2015, JCAP, 11, 018 [Google Scholar]
- Melchior, P., Sutter, P. M., Sheldon, E. S., Krause, E., & Wandelt, B. D. 2014, MNRAS, 440, 2922 [NASA ADS] [CrossRef] [Google Scholar]
- Micheletti, D., Iovino, A., Hawken, A. J., et al. 2014, A&A, 570, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nadathur, S., Carter, P. M., Percival, W. J., Winther, H. A., & Bautista, J. E. 2019, Phys. Rev. D, 100, 023504 [CrossRef] [Google Scholar]
- Nadathur, S., Woodfinden, A., Percival, W. J., et al. 2020, MNRAS, 499, 4140 [NASA ADS] [CrossRef] [Google Scholar]
- Pellejero-Ibanez, M., Chuang, C.-H., Rubiño-Martín, J. A., et al. 2017, MNRAS, 468, 4116 [NASA ADS] [CrossRef] [Google Scholar]
- Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1998, ApJ, 517, 565 [Google Scholar]
- Pisani, A., Sutter, P. M., Hamaus, N., et al. 2015, Phys. Rev. D, 92, 083531 [NASA ADS] [CrossRef] [Google Scholar]
- Pisani, A., Massara, E., Spergel, D. N., et al. 2019, BAAS, 51, 40 [Google Scholar]
- Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pollina, G., Hamaus, N., Paech, K., et al. 2019, MNRAS, 487, 2836 [NASA ADS] [CrossRef] [Google Scholar]
- Potter, D., Stadel, J., & Teyssier, R. 2017, Comput. Astrophys. Cosmol., 4, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Ratra, B., & Peebles, P. J. 1988, Phys. Rev. D, 37, 3406 [Google Scholar]
- Riess, A. G., Filippenko, A. V., Challis, P., et al. 1998, AJ, 116, 1009 [Google Scholar]
- Ronconi, T., & Marulli, F. 2017, A&A, 607, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sánchez, C., Clampitt, J., Kovacs, A., et al. 2017, MNRAS, 465, 746 [CrossRef] [Google Scholar]
- Schuster, N., Hamaus, N., Pisani, A., et al. 2019, JCAP, 12, 055 [CrossRef] [Google Scholar]
- Sheth, R. K., & Van De Weygaert, R. 2004, MNRAS, 350, 517 [NASA ADS] [CrossRef] [Google Scholar]
- Steeb, W.-H. H. 1997, Matrix Calculus and the Kronecker Product with Applications and C++ Programs (USA: World Scientific Publishing Co., Inc.) [CrossRef] [Google Scholar]
- Steinhardt, P. J., Wang, L., & Zlatev, I. 1999, Phys. Rev. D, 59, 123504 [Google Scholar]
- Stopyra, S., Peiris, H. V., & Pontzen, A. 2020, MNRAS, 500, 4173 [NASA ADS] [CrossRef] [Google Scholar]
- Sutter, P. M., Lavaux, G., Wandelt, B. D., & Weinberg, D. H. 2012, ApJ, 761, 44 [NASA ADS] [CrossRef] [Google Scholar]
- Sutter, P. M., Pisani, A., Wandelt, B. D., & Weinberg, D. H. 2014, MNRAS, 443, 2983 [NASA ADS] [CrossRef] [Google Scholar]
- Takahashi, R., Sato, M., Nishimichi, T., Taruya, A., & Oguri, M. 2012, ApJ, 761, 152 [Google Scholar]
- Tutusaus, I., Martinelli, M., Cardone, V. F., et al. 2020, A&A, 643, A70 [EDP Sciences] [Google Scholar]
- Verza, G., Pisani, A., Carbone, C., Hamaus, N., & Guzzo, L. 2019, JCAP, 12, 040 [CrossRef] [Google Scholar]
- Vielzeuf, P., Kovács, A., Demirbozan, U., et al. 2019, MNRAS, 500, 464 [Google Scholar]
- Voivodic, R., Rubira, H., & Lima, M. 2020, JCAP, 10, 033 [CrossRef] [Google Scholar]
- Wang, Y., & Zhao, G.-B. 2020, RAA, 20, 158 [NASA ADS] [Google Scholar]
- Zehavi, I., Zheng, Z., Weinberg, D. H., et al. 2011, ApJ, 736, 59 [NASA ADS] [CrossRef] [Google Scholar]
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ℓm (field perspective) or the C(ℓ) (estimator perspective).
A.1. The field perspective
We consider the spherical harmonics’ expansion coefficients , 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
where θ is the vector of the parameters and aℓ is a vector collecting the . The logarithm of this probability distribution function is
with . The computation of the second derivative, after taking the expectation value, yields
where S = ⟨A⟩ and we used the following matrix identities:
Now, we specify how the were collected into the vector aℓ. First of all, we recall that
where is the Poisson shot noise for probes combination AB in tomographic bins ij. Working at fixed multipole ℓ, the vector aℓ 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 −ℓ to ℓ
-
i varies from 1 to NA
-
A varies from 1 to 𝒩,
with 𝒩 being the number of probes and NA the number of tomographic bins for probe A. We chose to order the array aℓ varying the three indices with a significance increasing from left to right. With this choice the matrix S has the following block form:
Here ΣAB(ℓ) is the tomographic covariance matrix, that is to say, , and ⊗ is the Kronecker product. Equation (A.6) can be understood as follows: the matrix S is made up of diagonal blocks of size ℓ, each of them proportional to the identity 𝟙2ℓ+1 with a different factor
. Here the identity matrix 𝟙2ℓ+1 is exactly the Kronecker 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):
Using these properties, we have
where in the last line we used the fact that the shot noise is independent of the cosmological parameters, and the C(ℓ)’s block matrix C(ℓ) has the same form of Σ(ℓ) of Eq. (A.6), with entries instead of
. Then from Eq. (A.3), we can finally write
Now, we sum over the multipoles ℓ, considering them as independent
In order to obtain the expression Eq. (43), it sufficient to redefine the covariance as
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 fsky 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ℓm’s. In this section we derive the Fisher matrix expression when the observables are the estimator , defined as
One may argue that, since the aℓm are normally distributed, the are also normally distributed. However, this ansatz is wrong; in the proceeding equation, we obtain the
posterior distribution function and then their Fisher matrix. Since coefficients with different ℓ are uncorrelated, under the Gaussian assumption, we work at a fixed ℓ. The
likelihood is
Using the Dirac delta function δD, L(aℓ|θ) gets out from the integral
where is the result of the integral in the aℓm space9. Taking the logarithm of the previous equation we obtain
Proceeding as in the previous section and using ⟨A⟩=S the expectation value of the second derivative becomes
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.
All Tables
Cosmological parameters varied in the forecast analysis, together with their values in the fiducial cosmology assumed in this work.
Probe configurations in the pessimistic and optimistic scenarios for the combination V + WL, both considered for two values of rv, min.
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.
Marginalized 1-σ errors in different cosmological scenarios for WL + V + GCsp, when they are assumed to be independent, together with their combinations when the void-lensing cross-correlation is included.
All Figures
![]() |
Fig. 1. Figure containing the power spectra used in this work. Top left: nonlinear matter auto-power spectrum Pmm (solid red line), the void auto-power spectrum Pvv without the shot noise (dashed blue line), computed assuming the pessimistic void bias scenario with bv eff ≈ −11.9, and the absolute value of the void-matter cross-power spectrum Pvm (dot-dashed black line), all computed at redshift z = 0.001. Top right: tomographic void angular auto-power spectra |
In the text |
![]() |
Fig. 2. Galaxy and void distributions used in this work. Left: dimensionless normalized projected galaxy distributions |
In the text |
![]() |
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. |
In the text |
![]() |
Fig. 4. Fisher matrix marginalized contours for the (baseline) ΛCDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios. |
In the text |
![]() |
Fig. 5. Fisher matrix marginalized contours for the νΛCDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios. |
In the text |
![]() |
Fig. 6. Fisher matrix marginalized contours for the w0waCDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios. |
In the text |
![]() |
Fig. 7. Full contour plots when GCsp is combined with WL + V + XC for the w0waCDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios. |
In the text |
![]() |
Fig. 8. Fisher matrix marginalized contours for the νw0waCDM model, in the pessimistic (top) and in the optimistic (bottom) scenarios. |
In the text |
![]() |
Fig. 9. Forecast contours for the νw0waCDM model, in the pessimistic (left) and in the optimistic (right) scenarios. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.