The Euclid mission designed to map the geometry of the dark Universe presents an unprecedented opportunity for advancing our understanding of the cosmos through its photometric galaxy cluster survey. Central to this endeavor is the accurate calibration of the mass and redshiftdependent halo bias HB which is the focus of this paper. Our aim is to enhance the precision of HB predictions which is crucial for deriving cosmological constraints from the clustering of galaxy clusters. Our study is based on the peakbackground split PBS model linked to the halo mass function HMF and it extends it with a parametric correction to precisely align with results from an extended set of Nbody simulations carried out with the OpenGADGET3 code. Employing simulations with fixed and paired initial conditions we meticulously analyzed the matterhalo crossspectrum and modeled its covariance using a large number of mock catalogs generated with Lagrangian perturbation theory simulations with the PINOCCHIO code. This ensures a comprehensive understanding of the uncertainties in our HB calibration. Our findings indicate that the calibrated HB model is remarkably resilient against changes in cosmological parameters including those involving massive neutrinos. The robustness and adaptability of our calibrated HB model provide an important contribution to the cosmological exploitation of the cluster surveys to be provided by the Euclid mission. This study highlights the necessity of continuously refining the calibration of cosmological tools such as the HB to match the advancing quality of observational data. As we project the impact of our calibrated model on cosmological constraints we find that given the sensitivity of the Euclid survey a miscalibration of the HB could introduce biases in cluster cosmology analysis. Our work fills this critical gap ensuring the HB calibration matches the expected precision of the Euclid survey.
['cosmological parameters', 'cosmology: theory', 'large-scale structure of Universe', 'Astrophysics - Cosmology and Nongalactic Astrophysics']
Euclid preparation L. Calibration of the halo linear bias in vCDM cosmologies
{'L. Calibration of the halo linear bias in Λ ( ν ) CDM cosmologies': "Euclid Collaboration: T. Castro ⋆ 1 , 2 , 3 , 4 , A. Fumagalli 5 , 3 , R. E. Angulo 6 , 7 , S. Bocquet 8 , S. Borgani 9 , 3 , 1 , 2 , M. Costanzi 9 , 1 , 3 , J. Dakin 10 , K. Dolag 8 , P. Monaco 9 , 1 , 2 , 3 , A. Saro 9 , 3 , 1 , 2 , 4 , E. Sefusatti 1 , 3 , 2 , N. Aghanim 11 , L. Amendola 12 , S. Andreon 13 , C. Baccigalupi 3 , 1 , 2 , 14 , M. Baldi 15 , 16 , 17 , C. Bodendorf 18 , D. Bonino 19 , E. Branchini 20 , 21 , 13 , M. Brescia 22 , 23 , 24 , A. Caillat 25 , S. Camera 26 , 27 , 19 , V. Capobianco 19 , C. Carbone 28 , J. Carretero 29 , 30 , S. Casas 31 , M. Castellano 32 , G. Castignani 16 , S. Cavuoti 23 , 24 , A. Cimatti 33 , C. Colodro-Conde 34 , G. Congedo 35 , C. J. Conselice 36 , L. Conversi 37 , 38 , Y. Copin 39 , A. Costille 25 , F. Courbin 40 , H. M. Courtois 41 , A. Da Silva 42 , 43 , H. Degaudenzi 44 , G. De Lucia 1 , A. M. Di Giorgio 45 , M. Douspis 11 , X. Dupac 38 , S. Dusini 46 , M. Farina 45 , S. Farrens 47 , S. Ferriol 39 , P. Fosalba 48 , 49 , M. Frailis 1 , E. Franceschi 16 , M. Fumana 28 , S. Galeotta 1 , B. Gillis 35 , C. Giocoli 16 , 50 , P. Gómez-Alvarez 51 , 38 , A. Grazian 52 , F. Grupp 18 , 8 , L. Guzzo 53 , 13 , S. V. H. Haugan 54 , W. Holmes 55 , F. Hormuth 56 , A. Hornstrup 57 , 58 , S. Ili'c 59 , 60 , K. Jahnke 61 , M. Jhabvala 62 , B. Joachimi 63 , E. Keihänen 64 , S. Kermiche 65 , A. Kiessling 55 , M. Kilbinger 47 , B. Kubik 39 , M. Kunz 66 , H. Kurki-Suonio 67 , 68 , P. B. Lilje 54 , V. Lindholm 67 , 68 , I. Lloro 69 , E. Maiorano 16 , O. Mansutti 1 , O. Marggraf 70 , K. Markovic 55 , M. Martinelli 32 , 71 , N. Martinet 25 , F. Marulli 72 , 16 , 17 , R. Massey 73 , S. Maurogordato 74 , E. Medinaceli 16 , M. Melchior 75 , Y. Mellier 76 , 77 , M. Meneghetti 16 , 17 , E. Merlin 32 , G. Meylan 40 , L. Moscardini 72 , 16 , 17 , E. Munari 1 , 3 , S.-M. Niemi 78 , C. Padilla 79 , S. Paltani 44 , F. Pasian 1 , K. Pedersen 80 , W. J. Percival 81 , 82 , 83 , V. Pettorino 78 , S. Pires 47 , G. Polenta 84 , M. Poncet 85 , L. A. Popa 86 , L. Pozzetti 16 , F. Raison 18 , A. Renzi 87 , 46 , G. Riccio 23 , E. Romelli 1 , M. Roncarelli 16 , R. Saglia 8 , 18 , Z. Sakr 12 , 60 , 88 , J.-C. Salvignol 78 , A. G. Sánchez 18 , D. Sapone 89 , B. Sartoris 8 , 1 , M. Schirmer 61 , A. Secroun 65 , S. Serrano 48 , 90 , 49 , C. Sirignano 87 , 46 , G. Sirri 17 , L. Stanco 46 , J. Steinwagner 18 , P. Tallada-Crespí 29 , 30 , A. N. Taylor 35 , I. Tereno 42 , 91 , R. Toledo-Moreo 92 , F. Torradeflot 30 , 29 , I. Tutusaus 60 , L. Valenziano 16 , 93 , T. Vassallo 8 , 1 , G. Verdoes Kleijn 94 , Y. Wang 95 , J. Weller 8 , 18 , A. Zacchei 1 , 3 , G. Zamorani 16 , E. Zucca 16 , A. Biviano 1 , 3 , M. Bolzonella 16 , E. Bozzo 44 , C. Burigana 96 , 93 , M. Calabrese 97 , 28 , D. Di Ferdinando 17 , J. A. Escartin Vigo 18 , F. Finelli 16 , 93 , J. Gracia-Carpio 18 , S. Matthew 35 , N. Mauri 33 , 17 , A. Pezzotta 18 , M. Pöntinen 67 , C. Porciani 70 , V. Scottez 76 , 98 , M. Tenti 17 , M. Viel 3 , 1 , 14 , 2 , 4 , M. Wiesmann 54 , Y. Akrami 99 , 100 , V. Allevato 23 , S. Anselmi 46 , 87 , 101 , M. Archidiacono 53 , 102 , F. Atrio-Barandela 103 , A. Balaguera-Antolinez 34 , 104 , M. Ballardini 105 , 16 , 106 , D. Bertacca 87 , 52 , 46 , M. Bethermin 107 , 25 , A. Blanchard 60 , L. Blot 108 , 101 , H. Böhringer 18 , 5 , 109 , S. Bruton 110 , R. Cabanac 60 , A. Calabro 32 , G. Cañas-Herrera 78 , 111 , A. Cappi 16 , 74 , F. Caro 32 , C. S. Carvalho 91 , K. C. Chambers 112 , A. R. Cooray 113 , B. De Caro 28 , S. de la Torre 25 , G. Desprez 114 , A. Díaz-Sánchez 115 , J. J. Diaz 116 , S. Di Domizio 20 , 21 , H. Dole 11 , S. Esco ffi er 65 , A. G. Ferrari 33 , 17 , P. G. Ferreira 117 , I. Ferrero 54 , A. Finoguenov 67 , A. Fontana 32 , F. Fornari 93 , L. Gabarra 117 , K. Ganga 118 , J. García-Bellido 99 , T. Gasparetto 1 , V. Gautard 119 , E. Gaztanaga 49 , 48 , 120 , F. Giacomini 17 , F. Gianotti 16 , G. Gozaliasl 121 , 67 , C. M. Gutierrez 122 , A. Hall 35 , H. Hildebrandt 123 , J. Hjorth 80 , A. Jimenez Muñoz 124 , J. J. E. Kajava 125 , 126 , V. Kansal 127 , 128 , D. Karagiannis 129 , 130 , C. C. Kirkpatrick 64 , A. M. C. Le Brun 101 , J. Le Graet 65 , L. Legrand 131 , J. Lesgourgues 31 , T. I. Liaudat 132 , A. Loureiro 133 , 134 , G. Maggio 1 , M. Magliocchetti 45 , F. Mannucci 135 , R. Maoli 136 , 32 , C. J. A. P. Martins 137 , 138 , L. Maurin 11 , R. B. Metcalf 72 , 16 , M. Miluzio 38 , 139 , A. Montoro 49 , 48 , A. Mora 140 , C. Moretti 14 , 4 , 1 , 3 , 2 , G. Morgante 16 , S. Nadathur 120 , Nicholas A. Walton 141 , L. Pagano 105 , 106 , L. Patrizii 17 , V. Popa 86 , D. Potter 10 , I. Risso 142 , P.-F. Rocci 11 , M. Sahlén 143 , E. Sarpa 14 , 4 , 2 , A. Schneider 10 , M. Sereno 16 , 17 , A. Spurio Mancini 144 , 145 , J. Stadel 10 , K. Tanidis 117 , C. Tao 65 , N. Tessore 63 , G. Testera 21 , R. Teyssier 146 , S. Toft 147 , 148 , S. Tosi 20 , 21 , A. Troja 87 , 46 , M. Tucci 44 , C. Valieri 17 , J. Valiviita 67 , 68 , D. Vergani 16 , G. Verza 149 , 150 , and P. Vielzeuf 65 arXiv:2409.01877v1  [astro-ph.CO]  3 Sep 2024  \n(A ffi liations can be found after the references)  \nSeptember 4, 2024", 'ABSTRACT': 'The Euclid mission, designed to map the geometry of the dark Universe, presents an unprecedented opportunity for advancing our understanding of the cosmos through its photometric galaxy cluster survey. Central to this endeavor is the accurate calibration of the mass- and redshift-dependent halo bias (HB), which is the focus of this paper. We enhance the precision of HB predictions, which is crucial for deriving cosmological constraints from the clustering of galaxy clusters. Our study is based on the peak-background split (PBS) model linked to the halo mass function (HMF),  \nand extends it with a parametric correction to precisely align with results from an extended set of N -body simulations carried out with the OpenGADGET3 code. Employing simulations with fixed and paired initial conditions, we meticulously analyze the matter-halo cross-spectrum and model its covariance using a large number of mock catalogs generated with Lagrangian Perturbation Theory simulations with the PINOCCHIO code. This ensures a comprehensive understanding of the uncertainties in our HB calibration. Our findings indicate that the calibrated HB model is remarkably resilient against changes in cosmological parameters including those involving massive neutrinos. The robustness and adaptability of our calibrated HB model provide an important contribution to the cosmological exploitation of the cluster surveys to be provided by the Euclid mission. This study highlights the necessity of continuously refining the calibration of cosmological tools like the HB to match the advancing quality of observational data. As we project the impact of our calibrated model on cosmological constraints, we find that, given the sensitivity of the Euclid survey, a miscalibration of the HB could introduce biases in cluster cosmology analyses. Our work fills this critical gap, ensuring the HB calibration matches the expected precision of the Euclid survey.  \nKey words. Galaxy Clusters, Cosmology, Large-Scale Structure', '1. Introduction': "The structure formation process in the Universe is hierarchical, with smaller structures collapsing and merging to form larger ones. Galaxy clusters, the most massive virialized objects in the Universe, lie at the apex of this hierarchy. They serve as valuable cosmological probes, o ff ering insights into the growth of density perturbations and the geometry of the Universe (see, for instance, Allen et al. 2011; Kravtsov & Borgani 2012, for reviews).  \nThe cosmological exploitation of cluster surveys is primarily based on number count analysis. This method involves comparing the observed number of clusters in a survey, as a function of redshift and of a given observable quantity, to the theoretical prediction of the halo mass function (HMF) within a cosmological model, thus enabling the derivation of constraints on cosmological parameters. Numerous studies have been conducted in this area (e.g. Borgani et al. 2001; Holder et al. 2001; Rozo et al. 2010; Hasselfield et al. 2013; Planck Collaboration XX: Ade et al. 2014; Bocquet et al. 2015; Mantz et al. 2015; Planck Collaboration XXIV: Ade et al. 2016; Bocquet et al. 2019; Abbott et al. 2020; Costanzi et al. 2021; Lesci et al. 2022a).  \nComplementing cluster count analysis is the cluster clustering statistics, which examines their spatial distribution in the Universe (Mana et al. 2013; Castro et al. 2016; Baxter et al. 2016; To et al. 2021; Lesci et al. 2022b; Sunayama et al. 2023; Romanello et al. 2024; Fumagalli et al. 2024). The halo bias (HB) is a fundamental concept in this analysis since it reflects the ratio between the number overdensity of cluster and of the matter distribution. This relationship is expected to bring cosmological information through the mass- and redshift-dependence of the HB and to be linear on large scales, i.e., ≳ 100 Mpc, to guarantee the scale-independence of the HB.  \nThe Euclid mission (Laureijs et al. 2011; Euclid Collaboration: Scaramella et al. 2022; Euclid Collaboration: Mellier et al. 2024) is projected to provide significant advancements to cluster cosmology. Sartoris et al. (2016) have forecasted that the combined cluster count and clustering analysis by the Euclid mission will provide constraints on the amplitude of the matter power spectrum and the mass density parameters independent and competitive with other cosmological probes, underlining the potential of galaxy clusters as cosmological probes for ongoing and future missions.  \nAt the heart of cluster cosmology are the theoretical models for the HMF and HB. Simplified models based on linear perturbation theory and spherical collapse have provided invaluable insights into the potential of cluster counts and clustering as cosmological probes (see, e.g., Press & Schechter 1974; Bond et al. 1991). However, given the complexity and strongly non-linear nature of cluster formation dynamics, a refinement of these mod-  \nels to the precision level required by available and forthcoming surveys has to rely on cosmological simulations as the primary method to capture such complexity. Several studies have been dedicated to calibrating semi-analytical models for the HMF and HB, aiming to align these models' predictions with the results from extensive sets of simulations (see, for instance, Sheth & Tormen 1999; Sheth et al. 2001; Jenkins et al. 2001; Warren et al. 2006; Tinker et al. 2008, 2010; Bhattacharya et al. 2011; Watson et al. 2013; Despali et al. 2016; Comparat et al. 2017; Euclid Collaboration: Castro et al. 2023). These simulations not only accurately describe the gravitational interactions that predominantly drive structure formation but also attempt to account for the e ff ects of baryonic matter.  \nThe influence of baryons, albeit a minor component in the Universe's overall composition, plays a significant role in the formation of structures, particularly in the context of these simulations (Cui et al. 2014; Velliscig et al. 2014; Bocquet et al. 2016; Castro et al. 2020). Given the sensitivity of baryon evolution to the inclusion and modelling of astrophysical processes occurring at scales much smaller than those resolved in simulations, the modeling of baryonic feedback in hydrodynamical simulations remains a subject of active debate. At the scale of galaxy clusters, for instance, baryonic feedback is known to reorganize the mass density profile of halos without disrupting the structures, thereby altering the mass enclosed within a given radius compared to predictions from collisionless N -body simulations. Owing to the substantially greater computational demands of hydrodynamical simulations, it has become standard practice to derive the HMF from gravitational N -body simulations, with subsequent post-processing to account for baryonic e ff ects (see, e.g. Schneider & Teyssier 2015; Aricò et al. 2021). In this paper, we concentrate on the initial step of calibrating the HB using collisionless simulations. This approach is intended as a foundational phase, with the baryonic e ff ects being integrated later, employing the methodology akin to that used for the HMF (see, Castro et al. 2020; Euclid Collaboration: Castro et al. 2024). This strategy underscores our commitment to systematically exploring the cosmological parameter space, acknowledging the importance of baryonic e ff ects while methodically building toward their inclusion in our analysis.  \nSystematic errors in the calibration of the HMF and HB can significantly impact the final cosmological constraints. Studies like those by Salvati et al. (2020), Artis et al. (2021), and Euclid Collaboration: Castro et al. (2023) have highlighted how inaccuracies in theoretical models can propagate biases into cosmological parameter inferences. In response to these challenges, Euclid Collaboration: Castro et al. (2023) presented a new, rigorously studied calibration of the HMF based on a suitably designed set of N -body simulations, o ff ering the required accuracy to analyze Euclid cluster count data.  \nSemi-analytical modeling typically starts with a simplified physical model, such as the peak-background split (PBS, Mo & White 1996), which is then extended and refined by adding more degrees of freedom. These additional degrees of freedom are subsequently fitted to simulations. Conceptually, PBS links the HB to the HMF by decomposing the density field into high- and low-frequency modes. The high-frequency modes that cross the collapsing barrier describe the collapse of structures. In contrast, the low-frequency modes modulate the density field fluctuations, thereby enhancing the number of peaks that cross the collapse threshold, therefore linking the clustering of collapsed objects with the local density field. Despite its qualitative consistency with simulations, PBS must enhance its quantitative precision, especially in the context of the Euclid mission's requirements.  \nIn this paper, we address the challenge of enhancing the accuracy of HB predictions to the level required to fully exploit the cosmological potential of the 2-point clustering statistics from the Euclid photometric cluster survey. Our approach involves calibrating a semi-analytical model to quantify the discrepancies between PBS predictions and simulation results. This calibration aims to refine HB predictions, improving the reliability of cosmological parameter estimation derived from cluster counts and clustering.  \nThis paper is organized as follows: we revisit the theoretical aspects used in this paper in Sect. 2. In Sect. 3, we describe the methodology used in our analysis. We present the HB model and its calibration in Sect. 4, along with an assessment of our model's impact in a forecast Euclid cluster cosmology analysis. Final remarks are made in Sect. 5. The implementation of our model is publicly available in /github https: // github.com / TiagoBsCastro / CCToolkit and presented in Sect. 6.", '2.1. The halo mass function': "The di ff erential HMF is given by  \nd n d M d M = ρ m M ν f ( ν ) d ln ν , (1)  \nwhere n is the comoving number density of halos with mass in the range [ M , M + d M ], ν is the peak height, the function ν f ( ν ) is known as the multiplicity function, and ρ m is the comoving cosmic mean matter density  \nρ m = 3 H 2 0 Ω m , 0 8 π G , (2)  \nwhere H 0 and Ω m , 0 are the current value of the Hubble parameter and the matter density parameter, and G is the gravitational constant. The peak height is defined as ν = δ c /σ ( M , z ), where δ c is the critical density for spherical collapse (Peebles 2020) and σ 2 ( M , z ) is the filtered mass variance at redshift z , so that it measures how rare a halo is. The mass variance is expressed in terms of the linear matter power spectrum P m( k , z ) as  \nσ 2 ( M , z ) = 1 2 π 2 Z ∞ 0 d k k 2 P m( k , z ) W 2 ( k R ) , (3)  \nwhere R ( M ) = GLYPH<2> 3 M / (4 πρ m) GLYPH<3> 1 / 3 is the Lagrangian radius of a sphere containing the mass M , and W ( k R ) is the Fourier transform of a top-hat filter of radius R .  \nThe multiplicity function is considered universal if its cosmological dependence is solely through the peak height. However, numerous studies based on N -body simulations have challenged this assumption. These analyses reveal that, while the initial approximation of HMF universality is generally valid, systematic deviations from this universality become evident at late times in the universe's evolution. This deviation has been demonstrated in various independent investigations, each indicating a nuanced understanding of the HMF's behavior (e.g., Crocce et al. 2010; Courtin et al. 2011; Watson et al. 2013; Diemer 2020; Ondaro-Mallea et al. 2021; Euclid Collaboration: Castro et al. 2023).  \nThe non-universality of the HMF is a ff ected by both the halo definition and the residual dependence of δ c on cosmology. Various studies have shown this dependence, including Watson et al. (2013), Despali et al. (2016), Diemer (2020), and OndaroMallea et al. (2021) for the dependence on the halo definition, and Courtin et al. (2011) for the cosmology dependence of δ c. In our study, we define halos as spherical overdensities (SO) with an average enclosed mean density equal to ∆ vir( z ) times the background density, where ∆ vir( z ) is the non-linear density contrast of virialized structures as predicted by spherical collapse (Eke et al. 1996; Bryan & Norman 1998). The multiplicity function for halo masses computed at the virial radius has been shown to preserve universality better than other commonly assumed definition of halo radii (Despali et al. 2016; Diemer 2020; OndaroMallea et al. 2021). As for δ c, we use the fitting formula introduced by Kitayama & Suto (1996) that ignores the e ff ect of massive neutrinos; however, for the adopted values for total neutrino masses in this work, the fitting formula is still percent level accurate (LoVerde 2014).  \nIn this paper, we use the HMF presented in Euclid Collaboration: Castro et al. (2023)  \nν f ( ν ) = A ( p , q ) r 2 a ν 2 π e -a ν 2 / 2 1 + 1 ( a ν 2 ) p ! ( ν √ a ) q -1 , (4)  \nwhere the parameters { a , p , q } depend on background evolution and power spectrum shape as  \na = aR Ω az m ( z ) , (5)  \np = p 1 + p 2 d ln σ d ln R + 0 . 5 ! , (6)  \nq = qR Ω qz m ( z ) , (7)  \nwhere Ω m is the fractional density of matter in the Universe as a function of redshift, encompassing both baryonic and dark matter contributions  \naR = a 1 + a 2 d ln σ d ln R + 0 . 6125 ! 2 , (8)  \nqR = q 1 + q 2 d ln σ d ln R + 0 . 5 ! . (9)  \nLastly, the normalization parameter A is not a free parameter but a function of the other parameters  \nA ( p , q ) = ( 2 -1 / 2 -p + q / 2 √ π GLYPH<20> 2 p Γ GLYPH<18> q 2 GLYPH<19> + Γ GLYPH<18> -p + q 2 GLYPH<19>GLYPH<21> ) -1 , (10)  \nwhere Γ denotes the Gamma function. The adopted values for the HMF parameters are presented in Table 4 of Euclid Collaboration: Castro et al. (2023) and depend on the halo-finder used. In this work, we mostly use the ROCKSTAR calibration. The SUBFIND calibration is also used in Sect. 4.1.2 to assess the impact of the halo-finder in our model.", '3. Methodology': 'The overdensity of halos of mass M at the position r at redshift z  \nδ h( r , M , z ) = n ( r , M , z ) / ¯ n ( M , z ) -1 , (11)  \nis expressed in terms of the corresponding local number halo number density, n ( r , M , z ), and of the cosmic mean number density of such halos, ¯ n ( M , z ). In linear theory, it is related to the matter density contrast δ m( r , z ) as  \nδ h( r , M , z ) = b ( M , z ) δ m( r , z ) + ϵ ( r , M , z ) , (12)  \nwhere b ( M , z ) is the linear halo bias and ϵ is a stochastic term that in the following we assume to be associated with shot-noise.  \nIt follows from Eq. (12) that the halo-halo, P h, and halomatter power spectrum, P hm, are written as a function of the linear matter power spectrum, P m, for su ffi ciently large scales, as  \nP h( k , M , z ) = b 2 ( M , z ) P m( k , z ) + P SN , (13)  \nP hm( k , M , z ) = b ( M , z ) P m( k , z ) , (14)  \nwhere P SN represents the shot-noise component. Under the assumption that halos o ff er a discrete Poisson sample of the underlying continuous matter density field, P SN denotes a shot-noise component commonly assumed to be equivalent to the Poisson term, P SN = 1 / ¯ n , where ¯ n is the mean number density of tracers. On the other hand, halos are known to correspond with highdensity peaks of the underlying matter distribution. Therefore, they are expected not to provide a purely Poissonian sampling of this continuous density field. In fact, Casas-Miranda et al. (2002) and Hamaus et al. (2010) showed that positive and negative corrections to Poisson shot-noise are expected for low- and highmass halos. In this paper, we parameterize P SN as  \nP SN = 1 -α ¯ n , (15)  \nwhere α , a fitting parameter that we calibrate through simulations, controls the deviation from the assumption of Poisson noise.  \nAssuming the universality of the HMF, Mo & White (1996) derived the HB b ( M , z ) directly from the HMF by following the PBS framework. The PBS prediction for the HB as a function of the peak height ν reads  \nb PBS( ν ) = 1 -1 δ c d ln ν f ( ν ) d ln ν . (16)  \nAlthough the PBS provides an estimate of the bias that presents the correct qualitative behavior, Tinker et al. (2010) claims a relatively poor performance of the PBS in reproducing results from N -body simulations, with an accuracy of about 10-20 percent.  \nGiven the correct qualitative behavior of the PBS prescription, in this paper, we aim to improve the prediction of the bias by calibrating a model for the bias as a function of b PBS - i.e., we assume Eq. (16) to be valid, with the HMF (4) - and model its di ff erence with respect to the simulations.', '3.1.1. N -body simulations': 'Table 1 presents the adopted values for the matter density parameter and baryonic density parameter at redshift 0 ( Ω m , 0 and Ω b , 0), the dimensionless Hubble parameter h , the spectral index of the primordial power spectrum n s, and the amplitude of matter density fluctuations on scales of 8 h -1 Mpc σ 8 for the N -body simulations used in this work. We extended the set of PICCOLO simulations introduced and used by Euclid Collaboration: Castro et al. (2023) to calibrate the HMF model. We maintain the same technical configurations as the original PICCOLO simulations and refer to the above-mentioned HMF paper for further details while summarizing the main aspects. The set comprises 69 cosmological boxes, each with a comoving size of 2000 h -1 Mpc, and 4 × 1280 3 dark matter particles. The simulations were conducted using OpenGADGET3 , with initial conditions generated by monofonIC (Michaux et al. 2021), based on third-order Lagrangian Perturbation Theory (3LPT) at a redshift of z = 24. The adopted gravitational softening is equivalent to one-fortieth of the mean inter-particle distance.  \nThe original PICCOLO set of simulations included nine di ff erent choices for cosmological parameters, which were randomly chosen from the 95 percent confidence level hyper-volume of the joint SPT and DES cluster abundance constraints (Costanzi et al. 2021). Those represent the cosmologies C 0 to C 8 in Table 1. To further stress our modeling and guarantee its robustness, we also add the cosmologies C 9 and C 10, which sample the ( Ω m , 0 , σ 8) plane in the direction orthogonal to the degeneracy direction of the constraints from Costanzi et al. (2021), and in significant tension with such constraints. For each cosmology, two white-noise realizations were created to generate initial conditions. For each noise realization, a pair was generated by fixing the amplitudes of the Fourier modes of the density fluctuation field and pairing the phases (Angulo & Pontzen 2016), except for the reference C 0 model, which had 20 realizations.  \nWefurther add three simulations with Einstein-de Sitter cosmology (EdS, i.e., Ω m = 1 and ΩΛ = 0) with power-law initial matter power spectrum P m( k ) ∝ k n s with n s ∈ {-1 . 5 , -2 . 0 , -2 . 5 } . Those simulations, which have the same box size and the same number of particles as the PICCOLO set, are only instrumental for the modeling and will not be used for the calibration as they are far from the regime used to calibrate the model of Euclid Collaboration: Castro et al. (2023).  \nLastly, we carried out three pairs of simulations with massive neutrinos, again using the same box size as the PICCOLO set. Each pair has a total neutrino mass P m ν ∈ { 0 . 15 , 0 . 30 , 0 . 60 } eV. The simulation set-up for the neutrino simulations is the same used for the OpenGADGET3 simulations extensively validated in Adamek et al. (2023). The baseline cosmological parameters is the C 0 set, with Ω ν, 0 subtracted from Ω m , 0. The simulations have the same primordial amplitude A s, resulting in a lower σ 8 for increasing neutrino mass (see Table 1).  \nThe initial conditions (ICs) for the neutrino simulations were generated using the FastDF (Elbers 2022) implementation to monofonIC (Michaux et al. 2021). 1 The forked repository with the FastDF integration can be found here. 2 We employed the same number of neutrino particles as the number of grid resolution elements used for the cold dark matter particles. The total  \nTable 1. The cosmological parameters of the PICCOLO set of simulations.  \nNotes. The parameters of the C 0 to C 8 models have been uniformly drawn from the 95 per cent confidence level hyper-volume of the cluster abundance constraints presented in Costanzi et al. (2021). The parameters for the C 9 and C 10 were specifically selected to sample the parameter space in the direction orthogonal to the direction of degeneracy in the ( Ω m , 0 , σ 8) plane of the constraints of Costanzi et al. (2021), and in significant tension with them.  \nneutrino mass specified is attributed to a single massive neutrino species.', '3.1.2. Approximate methods: PINOCCHIO': 'In this paper, we also analysed 200 (100 pairs with each pair having fixed amplitudes and paired phases) halo catalogs simulated with the approximate LPT-based method implemented in the PINOCCHIO code (Monaco et al. 2002, 2013; Munari et al. 2017). All these simulations have been carried out under the assumption of the C 0 cosmological parameters. The rationale for this extra set of simulated catalogs is to model the impact of fixing and pairing the Fourier mode amplitudes in the ICs on the cluster clustering.', '3.2. Halo finders': 'Euclid Collaboration: Castro et al. (2023) showed that the halofinder adopted for the analysis of the N -body simulations can significantly alter the HMF. To understand if the halo-finder also impacts the HB, we selected two algorithms to extract halo catalogs: ROCKSTAR (Behroozi et al. 2013a) 3 and SUBFIND (Springel et al. 2001; Dolag et al. 2009; Springel et al. 2021). Although all these algorithms rely on the SO method to define halo boundaries, they di ff er in the method used to identify the center from which the spheres are grown and the criteria to classify between structures and sub-structures.  \nROCKSTAR divides the simulation volume into 3D friend-offriends (FOF; see, for instance, Davis et al. 1985) groups and runs a recursive 6D FOF algorithm on each group to create a hierarchy of FOF subgroups. Halo centers are determined by averaging the positions of particles in the innermost subgroup. To im-  \nprove consistency, we apply the CONSISTENT algorithm, which dynamically tracks halo progenitors, to the extracted ROCKSTAR catalogs as demonstrated in Behroozi et al. (2013b). SUBFIND also determines halo centers using a parallel implementation of the 3D FOF algorithm but directly assigns it to the particle with the lowest gravitational potential.  \nAmong the halo-finders studied by Euclid Collaboration: Castro et al. (2023), ROCKSTAR and SUBFIND are good representatives of the heterogeneity of possible results as they are close to the extremes, with SUBFIND suppressing the abundance of objects more massive than 10 13 M ⊙ h -1 by roughly 10 percent. See Knebe et al. (2011) for a more detailed comparison between the halo-finding algorithms.', '3.3. Measuring the halo bias': 'To measure the HB, we bin the halo distribution in log 10 ( M vir / M ⊙ h -1 ) with equispaced width intervals of 0 . 1 dex at each redshift. If the number of halos inside a bin was less than 10 000, we merged it with its neighbor to avoid having bins where the power-spectrum measurements were primarily dominated by shot noise. We measured the cross-spectrum P hm between the halos in each mass bin and the matter distribution traced by particles. The bias is then obtained as the ratio between this cross-spectrum and the matter power spectrum P m( k ). The matter density field was computed at the initial conditions and rescaled according to the linear growth factor for the simulations without massive neutrinos. For the simulations with massive neutrinos, the matter density field was built from the particle data from the same snapshot from which the halo catalog was extracted to consider the scale-dependent growth factor. We used the PYLIANS 4 Python libraries to construct the density field and compute power spectra on a 1024 3 piecewise cubic spline mesh grid. PYLIANS averages the power spectra in k -space shells with the width given by the fundamental mode of the box, k f ≡ 2 π/ L . Following Castro et al. (2020), we only considered modes with k values smaller than 0 . 05 h Mpc -1 to measure the HB, to ensure the validity of the linear approximation. The maximum k used for the measurements corresponds to the 16 th harmonic of the box and is much smaller than the Nyquist frequency of the grid used to compute the power spectrum.  \nTo calculate the bias for each mass bin, we used the ratio of the halo-matter cross-spectrum and the matter power spectrum  \nb sim i , j = P hm( kj ) i / P m( kj ) , (17)  \nwhere i and j are the mass bin and the k -shell indexes, respectively.', '3.4. Halo bias calibration': "We used a Bayesian approach with uninformative uniform priors on all parameters to fit our model for the linear HB parameters to simulation results. The best fits were obtained using the Dual Annealing method to find the posterior maximum as implemented by Virtanen et al. (2020), and covariance between parameters was estimated using PyMC (Salvatier et al. 2016). 5 The No-U-Turn sampler (NUTS) (Ho ff man et al. 2014) was automatically assigned internally by PyMC to sample the likelihood.  \nWe assume a Gaussian likelihood for the bias, since the power spectrum estimation for a Gaussian field realization follows a χ 2 -distribution when averaged over a shell. Since the  \nnumber of modes Nk inside the shell increases rapidly with k , this distribution approaches a Gaussian distribution by the central limit theorem. However, the number of modes is small for the first bins, and deviation from the Gaussian approximation is evident. To avoid this issue, we re-bin the measurement of the first 3 k -bins by merging them, ensuring that the bin with the fewest modes still contains 117 modes to recover the Gaussian approximation's validity.  \nNote that di ff erently than Castro et al. (2020), we used ICs with fixed amplitudes. Therefore, the distribution of the simulated bias is not approximately a ratio of two Gaussian distributions but approximately Gaussian itself since the denominator of Eq. (17) is not a random variable. The variance of the shellaverage halo-matter cross-spectrum on simulations with random Gaussian initial conditions is given by  \nσ P hm P m ! 2 = 1 Nk b 2 + 1 -α ¯ n P m ! , (18)  \nwhere we assume that the shot-noise contribution to the halo power spectrum follows Eq. (15) and a linear HB b . However, Zhang et al. (2022) showed that the predictions for random Gaussian initial conditions overestimate the variance observed for biased tracers on simulations with fixed amplitudes. Therefore, we modify Eq. (18) as follows  \nσ P hm P m ! 2 ≡ σ 2 b = 1 Nk β b 2 + 1 -α ¯ n P m ! + b 2 σ 2 sys . (19)  \nHere, β and σ sys are parameters we marginalize over, which control the variance suppression and the relative error due to the limited accuracy of the HMF used as the backbone for the PBS prescription. We note that, based on the halo model (see, for instance Cooray & Sheth 2002), one should expect a value for β close to zero as in the limit where all halos are considered, the shot-noise term on the right-hand side of Eq. (19) tends to zero, so that one should recover the matter power spectrum that by construction has zero variance. Furthermore, the HMF presented in Euclid Collaboration: Castro et al. (2023) was shown to have percent-level accuracy; thus, σ sys is expected to assume similar values during the calibration, presuming the PBS framework is valid. However, it is crucial to note that should σ sys significantly deviate from zero, such an occurrence could indicate a potential violation or limitation within the PBS framework, underscoring the necessity for careful interpretation of these parameters.  \nWe obtain the total log-likelihood by summing up all mass bins, modes, redshifts, and simulations. Following Euclid Collaboration: Castro et al. (2023), we use the redshifts z ∈ { 2 . 00 , 1 . 25 , 0 . 90 , 0 . 52 , 0 . 29 , 0 . 14 , 0 . 0 } , translating to a timespacing of about 1 . 7 Gyr. This spacing is larger than the characteristic dynamical time of galaxy clusters and e ff ectively suppresses the correlation between the results of di ff erent snapshots. Similarly, we assume that the correlations between di ff erent mass bins and modes are negligible. This assumption is justified by linear theory, which posits that di ff erent modes evolve independently in the linear regime, thus minimizing their mutual influence. In our analysis, we have three fitting parameters: α , β , and σ sys, in addition to the parameters of the bias model that are subject to calibration (to be discussed in Sect. 4.1.3). This approach allows for a comprehensive calibration of the bias model, taking into account the shot-noise correction α , the suppression of variance β , and the systematic uncertainties σ sys inherent in our method.", "3.5. Forecasting Euclid's cluster counts and cluster clustering observations": "To understand the impact of the HB calibration on cosmological constraints, it is important to realistically forecast the cosmological information to be extracted from the Euclid photometric cluster survey. For this purpose, we first quantify the impact of the HB on cluster counts and cluster clustering analyses. More precisely, the HB enters the modeling of the cluster counts covariance, which we compute analytically following the model of Hu & Kravtsov (2003), as validated in Fumagalli et al. (2021). Regarding the cluster clustering, the HB enters both in the computation of the mean value (power spectrum or two-point correlation function), and in the associated covariance matrix; in this work, we test the e ff ect on the real-space two-point correlation function and its covariance, following the model presented and validated in Euclid Collaboration: Fumagalli et al. (2024).  \nAfter assessing the impact on the two statistics, we forecast how the accuracy of the HB calibration propagates on the cosmological constraints obtained by cluster counts and cluster clustering experiments. We generate synthetic cluster abundance data as described in Section 2.5 of Euclid Collaboration: Castro et al. (2023), assuming the HMF calibrated and the HB of this work as benchmarks. Through a likelihood analyses, we constrain the cosmological parameters Ω m , 0 and σ 8, and the mass-observable relation (MOR) parameters A λ, B λ, C λ, D λ (see Section 2.5 of Euclid Collaboration: Castro et al. 2024), assuming flat priors for all the parameters. The MOR parameters describe the optical richness λ distribution as a function of the halo mass according to  \n⟨ ln λ | M vir , z ⟩ = ln A λ + B λ ln M vir 3 × 10 14 h -1 M ⊙ ! + C λ ln H ( z ) H ( z = 0 . 6) ! , (20)  \nwhere H ( z ) denotes the Hubble parameter at redshift z . The range for richness λ is considered to be between 20 and 2000, with the variation in logarithmic richness for a given virial mass and redshift expressed as a log-normal scatter  \nσ 2 ln λ | M vir , z = D 2 λ . (21)  \nThe reference values for the parameters are A λ = 37 . 8, B λ = 1 . 16, C λ = 0 . 91, and D λ = 0 . 15. These values have been obtained refitting the parameters presented in Saro et al. (2015) to the virial mass definition, with the assumption that halo profiles follow an NFW profile with a mass-concentration relationship as outlined in Diemer & Joyce (2019). The parameter values adopted in this work are consistent with the model presented by Castignani & Benoist (2016) to assign cluster membership probabilities to galaxies from photometric surveys.  \nWe perform the analysis on the synthetic catalogs, comparing them with the predictions made by using our bias and the one of Tinker et al. (2010), and compare the resulting posteriors with two estimators: we quantify the broadening / tightening of the posterior's amplitude by computing the di ff erence of the figure of merit ( ∆ FOM; see Huterer & Turner 2001; Albrecht et al. 2006) in the Ω m , 0 σ 8 plane, and the shift of the posterior's position by computing the posterior agreement (Bocquet et al. 2019), which determines whether the di ff erence between two posterior distributions is consistent with zero.  \nFig. 1. The constraints on the parameters in Eq. (19) fitted to the unbiased standard deviation of the 200 PINOCCHIO mocks with C 0 cosmology. We fitted α and β considering k ≤ 0 . 05 h Mpc -1 assuming a Gaussian likelihood with error bar estimated from the measurements for 0 . 05 ≤ k / h Mpc -1 ≤ 0 . 2 assuming it to be constant in k and equivalent to the unbiased standard deviation.  \n<!-- image -->", '4.1.1. Biased tracers statistics on fixed and paired simulations': 'Before modeling and calibrating the HB to the simulations, we investigate the impact of the variance suppression technique on the halo matter cross-spectrum.  \nIn Fig.1, we present the constraints on the parameters in Eq. (19) fitted to the unbiased standard deviation of the 200 PINOCCHIO mocks with C 0 cosmology. We fixed σ sys to zero, as this exercise does not involve modeling errors on the bias. The parameters α and β were fitted using only modes with k ≤ 0 . 05 h Mpc -1 , under the assumption of a Gaussian likelihood. We estimated the error bars from the measurements within the range 0 . 05 ≤ k / h Mpc -1 ≤ 0 . 2, treating it as constant and equivalent to the unbiased standard deviation. This approach prevents the estimation of the expectation of the mean and the error from using the same data. Figure 1 reveals a positive correlation between the parameters, with both assuming positive values. As expected, the posterior for the β parameter peaks close to zero, validating the e ff ectiveness of the fixed amplitudes technique in reducing the variance of biased tracer statistics. Notably, the small value of α found in our analysis suggests that the shot noise is well modeled as Poissonian to within 1-3 percent.  \nIn Fig. 2, we present the relative di ff erence between Eq. (19) computed with the best-fit values of α and β , and the unbiased standard deviation of the measurements, for di ff erent mass bins and redshift values. We present the results for three values z ∈ 0 . 0 , 1 . 0 , 2 . 0 spanning the redshift range of interest, while for the masses, we present the first, the intermediate, and the last occupied bin for each redshift. We observe that the residuals of the  \nfit always oscillate around zero, with no statistically significant mass, redshift, or k dependence.  \nIn Fig. 3, we present the Pearson correlation coe ffi cient ρ between the measurements of the halo-matter cross-power spectrum on a given simulation and its paired realization. The correlation coe ffi cient between two random variables X and Y is defined as  \nρ ( X , Y ) = ⟨ ( X - ⟨ X ⟩ )( Y - ⟨ Y ⟩ ) ⟩ σ X σ Y , (22)  \nwhere σ X and σ Y are the standard deviation of the random variables X and Y . The cross-power spectrum was measured for k ≤ 0 . 05 h Mpc -1 , for di ff erent ranges of halo masses (as reported in each panel) and redshift values (di ff erent columns). We also present the correlation coe ffi cient between simulations that assumes uncorrelated white-noise realizations for comparison. We note that paired simulations do not show a statistically significant di ff erence in their correlation with respect to simulations that assume uncorrelated noise realizations. The same conclusion is obtained by running a p-value test on all mass and redshift bins. This result justifies the assumption that di ff erent simulations are, in fact, independent, even if they have the same amplitudes and paired phases. This conclusion aligns with the claims of Villaescusa-Navarro et al. (2018) that variance suppression techniques are unlikely to a ff ect the halo abundance distribution. As the shot-noise term in Eq. (19) dominates over the other terms for our sample selection, one could anticipate the independence of the bias results from fixed and paired simulations due to the independence of abundance fluctuations and modes paring.  \nLastly, we use the PINOCCHIO catalogs to assess the impact of neglecting the correlation between di ff erent mass bins and Fourier modes on the calibration likelihood presented in Sect. 3.4. We measured the bias on the PINOCCHIO catalogs by applying the same mass and modes binning we used in our calibration and measured the correlation matrix between di ff erent simulations. In line with the results presented in Euclid Collaboration: Fumagalli et al. (2024) for the two-point correlation function, we explicitly verified in our analysis that the o ff -diagonal terms are sub-dominant and of the order of 10 percent, validating our calibration likelihood.', '4.1.2. Impact of the halo finder on the peak-background split performance': "In Fig. 4, the impact of the halo finder on the PBS is examined through the bias ratio measured in halo catalogs generated from the same simulation using either ROCKSTAR or SUBFIND , alongside the corresponding PBS prescription. For the ROCKSTAR catalogs, the standard error of the mean is further illustrated using an additional 19 realizations, with the assumed cosmology being C 0. The analysis spans di ff erent redshifts within the z ∈ [0 , 2] range. It is observed that the PBS tends to underestimate the simulation-derived bias by approximately 10% at higher redshifts, though it shows improved accuracy at z = 0. Given the minimal impact of the choice of halo finder on PBS's overall accuracy, subsequent results will focus exclusively on the ROCKSTAR halo catalogs. Despite PBS's tendency to underpredict the HB compared to simulations, it is noteworthy that the deviation remains consistently within 5-15 percent across all explored values of ν and redshifts. Our future e ff orts will aim to refine the PBS model by addressing these discrepancies, with the objective of achieving a simulation-calibrated HB model that is precise to within a few percent.  \nFig. 2. The relative di ff erence between Eq. (19) best-fit and the unbiased standard deviation of the PINOCCHIO measurements. Di ff erent columns are for di ff erent redshifts, and the corresponding mass bins are shown in each panel. The vertical dotted line demarcates two distinct sets of measurements: those to the left of the line were utilized as data points in the parameter fitting process, while the scatter of the points to the right was analyzed to estimate the variance. The grey regions highlight areas within a 5% deviation from the expected values.  \n<!-- image -->", '4.1.3. Modelling the halo bias': "In Fig. 5, we present the mean of the ratio of the measured bias with respect to the PBS prescription for di ff erent simulations. In the left panel, we use all the C 0 runs and show the ratio of the bias as a function of ν/ (1 + z ) for redshifts 0, 0 . 29, and 1 . 25. The factor (1 + z ) is only used to scale the results from di ff erent redshifts of the C 0 model to the same range. In the middle panel, we show the mean ratio as a function of ν for the three EdS cosmologies. Lastly, in the right panel, we present the mean ratio as a function of the background evolution Ω m( z ) for the C 0, C 9, and C 10 cosmologies.  \nThe left panel of Fig. 5 illustrates that the performance of PBS is influenced by the cosmological background evolution, encapsulated by Ω m( z ), yet appears una ff ected by variations in ν . This contrasts with the observations in the central panel, where PBS performance varies with both ν and changes in the power spectrum's shape. This sensitivity to ν is attributed to the limitations of the HMF calibration by Euclid Collaboration: Castro et al. (2023) when applied to an EdS cosmology that is far from its calibration regime, introducing artificial dependency.  \nHowever, it is important to note that, although these extrapolations to EdS scenarios are beyond the initial calibration range  \nof the HMF model, the model's accuracy is not disproportionately a ff ected across di ff erent values of n s. Indeed, in Euclid Collaboration: Castro et al. (2023), it has been demonstrated that the HMFmodel retains a consistent level of precision across various EdS cosmologies characterized by scale-free linear power spectra. Therefore, the dependence on the shape of the power spectrum is more likely related to the varying degree of accuracy of the PBS bias model as the shape of the power spectrum changes.  \nWe interpret the dependence of the PBS performance as a function of the shape of the power spectrum as follows. The extrapolation of the results on the central panel of Fig. 5 indicates that the PBS performance on EdS cosmologies with a steeper power spectrum ( n s < -2 . 5) degrades with decreasing n s. Within the PBS framework, the mass variance σ ( R L) smoothed on the scale of the Lagrangian patch R L is assumed to be dominated by the contribution of scales R LSS ≫ R L. For a power-law power spectrum, it is  \nd ln σ d ln R = -( n + 3) 2 . (23)  \nTherefore, the ratio σ ( R L) /σ ( R LSS) tends to unity as n s tends to -3. On the one hand, this explains why n s = -2 . 5 presents better performance than the other cases as one of the PBS assumptions  \nFig. 3. The correlation coe ffi cient ρ between the measurements of the halo-matter power spectrum, P hm( k ) in di ff erent simulations for k ≤ 0 . 05 h Mpc -1 , as a function of halo mass bin and redshift. The red histograms show the distribution of the correlation coe ffi cient between a simulation and its paired realization, while blue histograms are for the correlation coe ffi cient between simulations with uncorrelated white-noise realizations.  \n<!-- image -->  \nis better suited. On the other hand, for n s ≡ -3, perturbations on all scales reach the collapse at the same time, and it is no longer possible to distinguish between peaks and a large-scale modulation of a background perturbation, thus breaking PBS's fundamental assumption.  \nThe right panel of Fig. 5 shows that the residuals with respect to the PBS prediction increase linearly with the value of the density parameter Ω m( z ). While the slope of this linear dependence is similar for the three cosmologies, the normalization is a decreasing function of the clustering amplitude S 8. In fact, C 9 is the simulation with the lowest S 8 = GLYPH<16> p Ω m , 0 / 0 . 3 σ 8 GLYPH<17> = 0 . 438 while C 10 has the highest clustering amplitude S 8 = 1 . 07. The C 1 to C 8 simulations are not shown in this panel for better readability, but they cluster around C 0 as they have similar S 8 values.  \nThe better performance of the PBS in cosmologies with more clustering suggests that the di ff erence between this model and the simulation results is related to the connection between Lagrangian patches in the initial density field and the collapsed structures identified by the halo-finder. Collapsed structures stand out more clearly from the non-linear density field in more evolved and clustered cosmologies. For less clustered models, large halos are still forming and overlapping due to ongoing  \nmergers. This makes it more challenging to identify and link them to their corresponding Lagrangian patches clearly. Not surprisingly, in the EdS cosmologies, the PBS best performance is for n s = -2 . 5, where the evolution of the power spectrum is the steepest.  \nThe above line of reasoning suggests then that an SO algorithm could not be accurate in providing a one-to-one mapping between Lagrangian patches, destined to form virialized halos according to spherical collapse, and for which the PBS method predicts the bias, and halos identified in the non-linearly evolved density field. In this vein, since both ROCKSTAR and SUBFIND are based on an SO algorithm, it is not surprising that they predict similar deviations from PBS (see Fig. 4).  \nOn the other hand, we expect that collapsed structures have had more time to relax in cosmological models characterized by a higher value of S 8. As a consequence, they are more likely to be spherical. Again, this is in line with the better performance of the PBS on evolved cosmologies, as shown in the right panel of Fig. 5. Although suggestive, this interpretation of the deviations of PBS predictions would require a dedicated analysis to track their origin in detail, which goes beyond the scope of the analysis presented here.  \nFig. 4. The relative di ff erence between the bias measured in halo catalogs and the bias predicted by the PBS model of Eq. (16), at di ff erent redshifts in the range z ∈ [0 , 2]. Results refer to simulations carried out for the C 0 cosmology. In each panel, blue and red lines refer to the results obtained for the halo catalogs base on the application of ROCKSTAR and SUBFIND , respectively. For the ROCKSTAR catalogs, we also show the standard error of the mean using other 19 realizations of the same cosmology.  \n<!-- image -->  \nFrom the results shown in Fig. 5, it emerges that deviations from the PBS should depend on cosmic evolution, parameterized by Ω m( z ), on the slope of the linear power spectrum, and on the clustering amplitude S 8. 6 To capture such dependencies, we adopt the following description of the correction to the PBS prediction for the linear HB  \nb b PBS B f Ω m( z ) , d ln σ d ln R , S 8 ! = A 0 f 0 ( Ω m( z )) f 1 d ln σ d ln R ! f 2( S 8) . (24)  \nIn the above expression, we assume the following functional forms for the three above dependencies  \nf 0( x ) = 1 + a 1 x , (25)  \nf 1( x ) = 1 + b 1 x + b 2 x 2 , (26)  \nf 2( x ) = 1 + c 1 x . (27)  \nThe parameters A 0, a 1, b 1, b 2, and c 1 are calibrated in the next section through a detailed comparison with simulations. A balance between simplicity and empirical accuracy drives the parameterization chosen for these contributions. Specifically, we opted for a linear relationship for Ω m( z ) and S 8 while modeling the shape of the power spectrum using a quadratic function. In order to assess the potential parameter redundancy in using an extra parameter for the shape of the power spectrum, we performed Watanabe-Akaike Information Criterion (WAIC) (Watanabe 2010) and Pareto Smoothed Importance Sampling Leave-One-Out Cross-Validation (PSIS-LOO) (Vehtari et al. 2017) analyses comparing the model with b 2 free to vary with respect to b 2 fixed to 0. Both analyses confirmed that the fewer degrees of freedom adopted by the surrogate model do not compensate for the decrease in the model's predictability.  \nLastly, we report that we do not observe any significant correlation between the model prediction residual and other cosmological parameters assumed in the simulations. Thus, we conclude that the three fi components in Eq. (24) used in our analysis are su ffi cient to achieve our goal.", '4.2. Calibration of the halo bias correction': "In Fig. 6, we present the marginalized two-dimensional and unidimensional constraints on the model parameters presented in Eqs. (24) to (27). The best fit and 95 percent limits are reported in Table 2. We calibrate our model using the sub-set of 60 PICCOLO simulations covering the cosmological parameters from C 0 to C 10.  \nIn Fig. 7, we present the ratio between the mean of the observations on the 20 C 0 simulations with respect to our model predictions. Di ff erent rows correspond to di ff erent redshifts, while each panel corresponds to a di ff erent mass bin. The presented mass bins were selected as before to span from the least massive to the most massive occupied bin in that redshift. The shaded region in red corresponds to the error on the mean, assuming that each measurement follows Eq. (19). The shaded regions in grey correspond to two and four percent regions. As can be seen, our model's prediction presents a performance below two percent for  \nFig. 5. The mean of the ratio of the bias measured in simulations with respect to the PBS prescription. Left: the ratio of the bias as a function of ν/ (1 + z ) for di ff erent z , labeled by Ω m( z ), for all the C 0 runs. Center: the mean ratio as a function of ν for the three EdS cosmologies of pure power-law shapes of the linear power spectrum at z = 0. We report the value of the three spectral indexes in the inset. Right: the mean ratio as a function of the background evolution Ω m( z ) for the C 0, C 9, and C 10 cosmologies with varying S 8.  \n<!-- image -->  \nTable 2. The best fit and 95 percent limits on the model parameters presented in Eqs. (24) to (27).  \nNotes. We calibrate our model using the sub-set of PICCOLO simulations C 0C 10.  \ndi ff erent masses and redshift regimes when not primarily dominated by the sample variance, as it happens, for instance, for k ≲ 10 -2 h Mpc -1 . This accuracy holds over the k range up to which the onset of non-linearity occurs and the approximation of scale-independent bias breaks down, i.e., k ≳ 0 . 05 h Mpc -1 (marked by a vertical line).  \nWe note that at large k values, non-linearity e ff ects cause the bias measured from the simulations to take a scale dependence in exceeding the model prediction. As expected, this effect is smaller at higher redshift, consistent with the e ff ect of non-linearity shifting to larger k values.  \nSimilarly to Fig. 7, in Fig. 8, we present the ratio between the mean of the observations on the C 9 and C 10 simulations with respect to our model predictions. C 9 and C 10 are the simulations with lowest and highest S 8, respectively. Even for such extreme scenarios, our model performs well, thus confirming that our linear bias model, with the previously described calibration, can reproduce results from simulations with an accuracy of a few percent for Λ CDMcosmologies.", '4.3. Cosmologies with massive neutrinos': "We present our model's performance when considering simulations with massive neutrinos in Fig. 9. This allows us to assess the performance of our HB calibration for this minimal extension of Λ CDM. In this case, the simulation's bias has been computed by comparing it to the linear power spectrum of the corresponding cosmological model, which includes only the contributions from cold dark matter and baryons. For consistency, the same choice of considering only CDM and baryon contribution is also made for the computation of the HMF entering in our model for the HB (see, for instance, Castorina et al. 2014; Costanzi et al. 2013). From Fig. 9, it is evident that our model also precisely de-  \ncribes the bias in Λ ( ν )CDM models, despite the fact that such models have not been used during the HB calibration.  \nWenote that, unlike for the pure Λ CDMmodels, in this case, the measured bias underpredicts the model bias at large k . We also note that the dependence of this e ff ect on redshift, if any, goes in the direction of being larger at higher z . Also, there is some hint for it to be slightly smaller for smaller values of m ν and, therefore, of Ω ν . These e ff ects align with the expectation that such deviations are not dominated by non-linear evolution but rather by the e ff ect of neutrino-free streaming (Castorina et al. 2014).", '4.4. Comparison with previous models': 'In Fig. 10, we compare our model prediction with other models in the literature: Cole & Kaiser (1989), Sheth et al. (2001), Tinker et al. (2010), and Comparat et al. (2017). We present both our benchmark model as well as the PBS predictions based on the HMFmodel of Euclid Collaboration: Castro et al. (2023), which we use as a baseline of our model. Di ff erent columns correspond to di ff erent redshifts. The relative di ff erence to our benchmark model is presented in the panels in the second row. The predictions of the external models have been computed using the COLOSSUS toolkit (Diemer 2018). 7  \nTo ensure a fair comparison, we adopt the Planck-like C 0 cosmology, where the compared models have their peak performance among the PICCOLO cosmologies. All compared models have degraded performance as we move from this benchmark cosmology while we have shown the robustness of our model in Figs. 7, 8, and 9. That is due to the models assuming either the universality of the bias relation or a redshift-only dependence, while our method explicitly models the cosmology de-  \nFig. 6. Marginalized 68 and 95 confidence level contours on the model parameters presented in Eqs. (24) to (27). We calibrate our model using the sub-set of PICCOLO simulations C 0C 10. See Table 2 for the best fit and confidence levels.  \n<!-- image -->  \npendence of this relation. As such, comparisons between these models should be interpreted cautiously, as the underlying cosmology influences the exact figures.  \nAs for the comparison with the PBS prediction based on the HMF calibration by Euclid Collaboration: Castro et al. (2023), the results shown here confirm those shown in Fig.5: the PBSbased predictions underestimate our simulations-based calibration by about 5-10%, almost independently of ν . The models of Cole & Kaiser (1989) and Sheth et al. (2001) over- and under-estimate the bias by ∼ 10%, respectively. Their relatively poor performance is not surprising. The prediction by Cole & Kaiser (1989) corresponds to the PBS prediction when using the  \nHMFfrom Press & Schechter (1974). As the Press & Schechter (1974) HMF only qualitatively explains the abundance of halos, it is expected that the bias from the PBS will not perform much better. On the other hand, the Sheth et al. (2001) model was calibrated on simulations. However, such simulations had a resolution that allowed those authors to cover a dynamic range significantly narrower than that accessible to our simulations. The prediction by Cole & Kaiser (1989) di ff ers from ours by an amount almost independent of ν and redshift. On the other hand, the HB by Sheth et al. (2001) di ff ers from ours in a ν - and z -dependent way.  \nFig. 7. The ratio between the mean of the observations on the 20 C 0 simulations with respect to our model predictions. Di ff erent rows correspond to di ff erent redshifts, while each panel corresponds to a di ff erent mass bin. The shaded region in red corresponds to the error on the mean, assuming that each measurement follows Eq. (19). The shaded regions in grey correspond to two and four percent regions.  \n<!-- image -->  \nNotably, the model of Tinker et al. (2010) is only superior to the abovementioned models for low peak-height. The di ff erences with respect to our model grow with redshift and peak-height. This could be due to the heterogeneity of the simulations used to calibrate the model of Tinker et al. (2010) and the possible limitation of the model itself to capture the cosmological dependence of the HB. In this paper, we calibrate to a set of simulations that have been run with the same code and setup. On the other hand, Tinker et al. (2010) based their analysis on a collection of simulations carried out with di ff erent codes and configurations in terms of resolution and box sizes. Also, their model assumes a redshift dependency for the evolution, while from Fig. 6 we note that a parametrization of this evolution through Ω m( z ) is more universal.  \nLastly, the model by Comparat et al. (2017) shows a good agreement with our model. The most significant di ff erences are at low-redshift where the model of Comparat et al. (2017) predicts a bias that is smaller than ours by about 5 percent. This di ff erence reduces to a sub-percent at high redshift. The agreement is unsurprising as the model of Comparat et al. (2017) was also calibrated on ROCKSTAR catalogs.', '4.5. Impact on cluster cosmology analysis': 'In this section, we forecast the impact of the HB calibration on cosmological analyses of cluster counts and cluster clustering from Euclid . We present the results for the bias model of Tinker et al. (2010) and the model calibrated in this study. The rationale for assuming Tinker et al. (2010) is to use a widely used model in cluster cosmology representative of the di ff erence in the bias  \nmodels presented in Fig. 10. Nonetheless, we do not expect the results to change significantly if we assumed the model of Comparat et al. (2017) that presents a better concordance with our model at high redshift but worse at low redshift. Therefore, the overall impact on cosmological constraints will compensate partially as the clustering cosmological signal for Euclid peaks at lower redshifts.  \nAs described in Sect. 3.5, we assess the e ff ect of the HB calibration in a more realistic scenario, performing a likelihood analysis of both the individual analysis of number counts and cluster clustering statistics and the combination of these probes. In all scenarios, the observable-mass relation (Eq. 20) is calibrated by combining the probes with weak lensing (WL) mass estimates, assuming a constant error of 1%. The mass calibration is the primary source of systematic uncertainty in cluster cosmology studies, and a 1% calibration represents the goal for stage IV surveys. Therefore, the chosen setting o ff ers a forecast of the utmost cosmological bias resulting from inaccurate modeling of the HB. Lastly, we assume three independent Gaussian likelihoods (Fumagalli et al. 2024) for number counts, clustering, and WLmasses.  \nIn the left panel of Fig. 11, we start by presenting the percentage residuals of the number counts covariance. We show the full covariance matrix for the number counts analysis, with the mass dependence within each redshift bin. Notably, the impact of the HB model is minimal at low redshift but becomes significant, reaching up to 20%, at higher redshifts, especially in the high-mass bins. However, the impact of a di ff erent bias calibration is mitigated by the shot-noise contribution when the latter becomes dominant along the diagonal, as the HB only plays a  \nFig. 8. Similar to Fig. 7 but for the C 9 and C 10 cosmological parameters. Among the PICCOLO set, C 9 and C 10 correspond to the cosmologies with the lowest and the highest S 8, respectively.  \n<!-- image -->  \nFig. 9. Similar to Fig. 7 but for simulations with massive neutrinos. For better plot readability, we only show the uncertainties (red shaded regions) for the simulation with total neutrino masses equal to 0 . 15 eV.  \n<!-- image -->  \nFig. 10. Comparison between the HB predicted by our model with predictions from other models presented in the literature: Cole & Kaiser (1989), Sheth et al. (2001), Tinker et al. (2010), and Comparat et al. (2017). We present both our benchmark model as well as the PBS predictions based on the HMF model of Euclid Collaboration: Castro et al. (2023) used as a baseline of our model. Di ff erent columns correspond to di ff erent redshifts. The relative di ff erence with respect to our benchmark model is presented in the panels in the second row. We adopt a composite scale for the residual plot to show the dynamic range of di ff erences between the models: the scale is linear for values between [ -10 , 10] % and symmetric log outside. For reference, we show the zero line in black. The predictions of the models from the literature have been computed using the COLOSSUS toolkit (Diemer 2018).  \n<!-- image -->  \nrole in the computation of sample variance. To quantify the impact of such a discrepancy on parameter posteriors, we perform the likelihood analysis for a number counts experiment, as described in Sect. 3.5. From the comparison of the two posteriors, we obtain ∆ FOM = -0 . 67 and a perfect agreement between the positions of the two contours, meaning that the impact of the HB calibration is below other systematics.  \nAs for the analysis of cluster counts, in the central and right panels of Fig. 11, we present the percent residuals for the clustering covariance as a function of the radial separation in both a low-redshift (central panel) and high redshift (right panel) bin. Similar to our findings for the number counts, the most significant impact is observed on the o ff -diagonal elements, particularly at high redshift. However, in contrast to number counts covariance, the inclusion of shot-noise in cluster clustering contributes to the o ff -diagonal elements, helping to mitigate the effect of di ff erent HB calibrations across all scales. This results in a di ff erence that never exceeds 10%. Importantly, in the case of cluster clustering, the HB also a ff ects the expected signal - either the two-point correlation function or power spectrum - leading to a di ff erence independent of the radial separation but increasing with redshift, reaching a 10% level in the high-redshift interval.  \nThe cosmological forecasts from the clustering experiment show that the minimal variation in covariance terms produces a negligible di ff erence in the posterior amplitude, equal to ∆ FOM = 1 . 03. However, comparing the two posteriors reveals an agreement at only 0 . 68 σ . This implies that the di ff erences in the two-point correlation function translate into a sizeable shift in the cosmological constraints. Notably, the di ff erence in the posteriors induced by the di ff erent calibrations of the HB alone surpasses the 0 . 25 σ threshold commonly employed in other studies (Euclid Collaboration: Deshpande et al. 2023) to flag systematic errors that, if exceeded, could accumulate and lead to a collectively significant di ff erence.  \nThe combined analysis (cluster counts + cluster clustering) posteriors are shown in Fig. 12. We notice that assuming the HB calibration by Tinker et al. (2010) still causes a shift in the posteriors with respect to the HB calibration presented in this paper. This aligns with the forecast results for the cluster clustering analysis presented above. Although the di ff erence is reduced to 0 . 39 σ , the combination with number counts and weak lensing masses cannot compensate for the impact of the HB calibration that a ff ects mostly the cluster clustering.', '5. Conclusions': "This paper presents a calibrated semi-analytical model for the HB in view of the joint cosmological exploitation of number counts and clustering of galaxy clusters from the Euclid survey. Our approach begins with the PBS model, based on the HMF of Euclid Collaboration: Castro et al. (2023), and extends this approach by introducing a novel parametric correction. This correction is designed to align the PBS prediction with the results from an extended and homogenous set of N -body simulations that we carried out for both vanilla Λ CDM models, varying cosmological parameters, and for Λ ( ν )CDMmodels by varying the sum of the neutrino masses. The simulations employed fixed and paired initial conditions (see, Angulo & Pontzen 2016), providing a robust reduced variance framework for our analysis. We measured the HB by examining the ratio of the matterhalo cross-spectrum. Additionally, we modeled the covariance of these measurements using 200 mock catalogs of the Euclid cluster survey, based on the approximate LPT-based PINOCCHIO code. This ensures a thorough understanding of the uncertainties involved in our calibration of the HB. The key findings and implications of our study are summarized in the following points.  \n- -The use of fixed and paired initial conditions for the simulations analysed in our study proved highly advantageous  \nFig. 11. Percentage residuals of cluster counts ( left panel) and cluster clustering covariance matrices ( central and right panels), computed with the bias from Tinker et al. (2010) in comparison to the one calculated using the bias calibrated in this study (Eq. 24). We show the full covariance matrix for number counts in mass and redshift bins. For the two-point correlation function of galaxy clusters, we show two blocks of the full covariance (low and high redshift bins) as a function of the radial separation.  \n<!-- image -->  \nFig. 12. Parameter posteriors at 68 and 95 percent of confidence level obtained by analyzing number counts, weak lensing masses, and cluster clustering computed with the HB calibrated in this work (cyan contours) and the bias from Tinker et al. (2010) (orange contours). The error associated with weak lensing mass is set at 1% of the mass.  \n<!-- image -->  \nfor estimating the bias of tracers. By parametrizing the covariance of the bias measurements with two parameters one controlling the shot-noise contribution and the other for suppression due to fixing, respectively α and β in Eq. (19) - we observed significant e ff ectiveness in the variance suppression term. This was demonstrated in the constraints on the terms describing the variance in the halo-matter power spectrum, P hm( k ), shown in Fig. 1. Furthermore, our analysis of the measurements of P hm( k ) between paired simulations, as illustrated in Fig. 3, revealed no significant correlation between them. This finding underscores the e ffi cacy of the fixed and paired simulation approach in providing reliable estimates of the bias factor characterizing the distribu-  \non of tracers (i.e., halos), which is free from the influences of inherent correlations that could a ff ect the results.  \n- -The impact of the choice of the halo finder used in the analysis of the N -body simulations on the performance of the PBS is shown in Fig. 4. While comparing the ROCKSTAR and SUBFIND halo finders, it was observed that the PBS generally underestimates the measured bias from simulations, particularly at higher redshifts. However, the impact of the halo finder choice on the PBS prescription's performance is almost negligible, thus reinforcing the robustness of our approach in assessing the PBS performance across di ff erent redshift ranges.\n- -Our modeling of the HB, as illustrated in Fig. 5, reveals significant insights into the cosmological dependency of the performance of the PBS. The background cosmological evolution influences the PBS performance more than the peak height parameter ν . This was particularly evident in di ff erent cosmologies, where PBS's e ff ectiveness varied with the shape of the power spectrum and the degree of clustering evolution, as described by the S 8 parameter. Notably, in more clustered cosmologies, PBS improved its performance. This suggests a possible link between the ease of identifying collapsed structures in cosmologies with more evolved clustering and their corresponding Lagrangian patches. This result led us to develop a refined model for the PBS correction, expressed in Eq. (24), which incorporates terms depending on Ω m( z ), the local slope of the power spectrum, and S 8.\n- -The calibration of our model parameters, with best-fit presented in Fig. 6 and Table 2, demonstrated its robust performance across a range of cosmological conditions. Figures 7 and 8 illustrate our model prediction performance on the reference C 0 simulations and on the C 9 and C 10 simulations. The accuracy of our model is particularly noteworthy, always remaining below a two percent deviation for di ff erent masses and redshift regimes, with the possible exception of unrealistic cases largely influenced by sample variance. Quite remarkably, this level of precision is maintained even in extreme scenarios represented by the C 9 and C 10 simulations, which have the lowest and highest S 8 values, respectively.\n- -The robustness of our HB calibration is further demonstrated in scenarios involving massive neutrinos, as showcased in Fig. 9. Despite not incorporating massive neutrino simulations during the calibration phase, our model accurately predicts the HB in these cosmologies. Neutrinos are treated ac-  \ncording to the model presented by Castorina et al. (2014, see also Costanzi et al. 2013) and the measurements of the bias with respect to the matter power spectrum of cold dark matter and baryons, as it was done for the HMF in simulations with massive neutrinos in Euclid Collaboration: Castro et al. (2023). The ability of our model to adapt and perform reliably in such scenarios without the need for recalibration highlights its robustness and versatility.  \n- -As for the comparison with HB models already introduced in the literature (see Figure 10), the models by Cole & Kaiser (1989) and Sheth et al. (2001) show significant deviations from our results, likely due to their calibration on simulations covering narrower dynamic ranges and variety of cosmological models. The HB model by Tinker et al. (2010) shows increasing discrepancies with our model at higher redshifts and peak heights. Such di ff erences could be attributed to its calibration on a heterogeneous set of simulations and an inadequate accounting for the cosmological dependence of the HB. In contrast, the model by Comparat et al. (2017) aligns more closely with our findings, particularly at higher redshifts. This agreement is expected as their model was also calibrated using ROCKSTAR catalogs.\n- -As for the impact of changing the calibration of the HB on the derived cosmological posteriors, we showed in Fig. 11 the di ff erences of the covariance matrices for a Euclid cluster count and cluster clustering analysis using both our calibration and the one provided by Tinker et al. (2010). While the impact on number counts covariance is minimal at low redshifts, it becomes substantial, up to 20%, at higher redshifts. However, the presence of shot-noise in the analysis helps mitigating this e ff ect. In cluster clustering, we observed that the HB calibration could lead to di ff erences in the two-point correlation function, particularly at high redshifts. This difference can potentially bias cosmological constraints beyond the 0 . 25 σ threshold commonly used to flag significant systematic errors (Adamek et al. 2023). Moreover, the combined analysis of number counts, cluster clustering, and weak lensing masses demonstrates that even with these additional data, the calibration of HB can not be entirely compensated. This highlights the importance of precise HB calibration in cluster cosmology, especially for a survey like the one being provided by Euclid that will reach an unprecedented sensitivity and level of statistics.  \nIn summary, the analysis presented in this paper has systematically calibrated and tested the HB for a range of cosmological scenarios, demonstrating its critical impact on the precision of cosmological analyses based on galaxy clusters for the Euclid mission. The resilience of our HB model against variations of cosmological models, including the presence of massive neutrinos and di ff erent degrees of clustering amplitude, highlights its robustness and adaptability. Importantly, our model is robust against the halo finder definition, inheriting its dependence through the HMF only. This is a remarkable feature as the correspondence between halos in N -body simulations and real clusters in surveys remains a complex issue, with uncertainties in halo identification and characterization potentially influencing the extraction of cosmological parameters. Future research should focus on understanding and quantifying these uncertainties, especially concerning observational challenges such as projection e ff ects and the mass-observable relation. As we move forward, extending this precision to departures from the standard Λ ( ν )CDM framework will be crucial in fully harnessing the capabilities of next-generation cosmological surveys.", '6. Data availability': 'In Castro & Fumagalli (2024), we implement the model presented in this paper, together with the models for the HMF presented in Euclid Collaboration: Castro et al. (2023) and for the impact of baryonic feedback on cluster masses presented in Euclid Collaboration: Castro et al. (2024). The source code can be accessed in /github https: // github.com / TiagoBsCastro / CCToolkit.  \nAcknowledgements. It is a pleasure to thank Valerio Marra for constructive comments during the production of this work, Fabio Pitari and Caterina Caravita for support with the CINECA environment, Peter Berhoozi for the support with ROCKSTAR , Oliver Hahn for the support with monofonIC , and Luca Di Mascolo for the support with PyMC . TC is supported by the Agenzia Spaziale Italiana (ASI) under - Euclid-FASE D Attivita\' scientifica per la missione - Accordo attuativo ASI-INAF n. 2018-23-HH.0. SB, TC, PM, and AS are supported by the PRIN 2022 PNRR project "Space-based cosmology with Euclid: the role of High-Performance Computing" (code no. P202259YAF), by the Italian Research Center on High-Performance Computing Big Data and Quantum Computing (ICSC), project funded by European Union - NextGenerationEU - and National Recovery and Resilience Plan (NRRP) - Mission 4 Component 2, within the activities of Spoke 3, Astrophysics and Cosmos Observations, and by the INFN INDARK PD51 grant. TC and AS are also supported by the FARE MIUR grant \'ClustersXEuclid\' R165SBKTMA. AS is also supported by the ERC \'ClustersXCosmo\' grant agreement 716762. MC and TC are supported by the PRIN 2022 project EMC2 - Euclid Mission Cluster Cosmology: unlock the full cosmological utility of the Euclid photometric cluster catalog (code no. J53D23001620006). KD acknowledges support by the DFG (EXC-2094 - 390783311) as well as support through the COMPLEX project from the European Research Council (ERC) under the European Union\'s Horizon 2020 research and innovation program grant agreement ERC-2019-AdG 882679. We acknowledge the computing centre of CINECA and INAF, under the coordination of the \'Accordo Quadro (MoU) per lo svolgimento di attività congiunta di ricerca Nuove frontiere in Astrofisica: HPC e Data Exploration di nuova generazione\', for the availability of computing resources and support. We acknowledge the use of the HOTCAT computing infrastructure of the Astronomical Observatory of Trieste - National Institute for Astrophysics (INAF, Italy) (see Bertocco et al. 2020; Ta ff oni et al. 2020). 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 Agenzia Spaziale Italiana, the Austrian Forschungsförderungsgesellschaft funded through BMK, the Belgian Science Policy, the Canadian Euclid Consortium, the Deutsches Zentrum für Luft- und Raumfahrt, the DTU Space and the Niels Bohr Institute in Denmark, the French Centre National d\'Etudes Spatiales, the Fundação para a Ciência e a Tecnologia, the Hungarian Academy of Sciences, the Ministerio de Ciencia, Innovación y Universidades, the National Aeronautics and Space Administration, the National Astronomical Observatory of Japan, the Netherlandse Onderzoekschool Voor Astronomie, the Norwegian Space Agency, the Research Council of Finland, the Romanian Space Agency, the State Secretariat for Education, Research, and Innovation (SERI) at the Swiss Space O ffi ce (SSO), and the United Kingdom Space Agency. A complete and detailed list is available on the Euclid web site ( www.euclid-ec.org ).', 'References': 'Abbott, T. M. C., Aguena, M., Alarcon, A., et al. 2020, PRD, 102, 023509 Adamek, J., Angulo, R. E., Arnold, C., et al. 2023, JCAP, 06, 035 Albrecht, A., Bernstein, G., Cahn, R., et al. 2006, arXiv:0609591 Allen, S. W., Evrard, A. E., & Mantz, A. B. 2011, ARA&A, 49, 409 Angulo, R. E. & Pontzen, A. 2016, MNRAS, 462, L1 Aricò, G., Angulo, R. E., Contreras, S., et al. 2021, MNRAS, 506, 4070 Artis, E., Melin, J.-B., Bartlett, J. G., & Murray, C. 2021, A&A, 649, A47 Baxter, E. J., Rozo, E., Jain, B., Ryko ff , E., & Wechsler, R. H. 2016, MNRAS, 463, 205  \n- Behroozi, P. S., Wechsler, R. H., & Wu, H.-Y. 2013a, ApJ, 762, 109\n- Behroozi, P. S., Wechsler, R. H., Wu, H.-Y., et al. 2013b, ApJ, 763, 18\n- Bertocco, S., Goz, D., Tornatore, L., et al. 2020, in Astronomical Society of the Pacific Conference Series, Vol. 527, 303\n- Bhattacharya, S., Heitmann, K., White, M., et al. 2011, ApJ, 732, 122\n- Bocquet, S., Dietrich, J. P., Schrabback, T., et al. 2019, ApJ, 878, 55\n- Bocquet, S., Saro, A., Dolag, K., & Mohr, J. J. 2016, MNRAS, 456, 2361\n- Bocquet, S., Saro, A., Mohr, J. J., et al. 2015, ApJ, 799, 214\n- Bond, J. R., Cole, S., Efstathiou, G., & Kaiser, N. 1991, ApJ, 379, 440\n- Borgani, S., Rosati, P., Tozzi, P., et al. 2001, ApJ, 561, 13\n- Bryan, G. L. & Norman, M. L. 1998, ApJ, 495, 80\n- Casas-Miranda, R., Mo, H. J., Sheth, R. K., & Boerner, G. 2002, MNRAS, 333, 730  \nCastignani, G. & Benoist, C. 2016, A&A, 595, A111  \n- Sheth, R. K. & Tormen, G. 1999, MNRAS, 308, 119\n- Castorina, E., Sefusatti, E., Sheth, R. K., Villaescusa-Navarro, F., & Viel, M. 2014, JCAP, 02, 049\n- Castro, T., Borgani, S., Dolag, K., et al. 2020, MNRAS, 500, 2316\n- Castro, T. & Fumagalli, A. 2024, CCToolkit: A Python Package for Cluster Cosmology Calculations  \nCastro, T., Marra, V., & Quartin, M. 2016, MNRAS, 463, 1666  \nCole, S. & Kaiser, N. 1989, MNRAS, 237, 1127  \n- Comparat, J., Prada, F., Yepes, G., & Klypin, A. 2017, MNRAS, 469, 4157, [Erratum: MNRAS 483, 2561-2562 (2019)]  \nCooray, A. & Sheth, R. K. 2002, Phys. Rept., 372, 1  \nCostanzi, M., Saro, A., Bocquet, S., et al. 2021, PRD, 103, 043522  \n- Costanzi, M., Villaescusa-Navarro, F., Viel, M., et al. 2013, JCAP, 12, 012\n- Courtin, J., Rasera, Y., Alimi, J. M., et al. 2011, MNRAS, 410, 1911\n- Crocce, M., Fosalba, P., Castander, F. J., & Gaztanaga, E. 2010, MNRAS, 403, 1353  \nCui, W., Borgani, S., & Murante, G. 2014, MNRAS, 441, 1769  \nDavis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 1985, ApJ, 292, 371  \nDespali, G., Giocoli, C., Angulo, R. E., et al. 2016, MNRAS, 456, 2486  \n- Diemer, B. 2018, ApJ Suppl., 239, 35\n- Diemer, B. 2020, ApJ, 903, 87\n- Diemer, B. & Joyce, M. 2019, ApJ, 871, 168\n- Dolag, K., Borgani, S., Murante, G., & Springel, V. 2009, MNRAS, 399, 497  \nEke, V. R., Cole, S., & Frenk, C. S. 1996, MNRAS, 282, 263  \n- Elbers, W. 2022, JCAP, 11, 058\n- Euclid Collaboration: Castro, T., Borgani, S., Costanzi, M., et al. 2024, A&A, 685, A109\n- Euclid Collaboration: Castro, T., Fumagalli, A., Angulo, R. E., et al. 2023, A&A,  \n671, A100  \n- Euclid Collaboration: Deshpande, A. C., Kitching, T., Hall, A., et al. 2023, arXiv: 2302.04507\n- Euclid Collaboration: Fumagalli, A., Saro, A., Borgani, S., et al. 2024, A&A, 683, A253\n- Euclid Collaboration: Mellier, Y., Abdurro\'uf, Acevedo Barroso, J. A., et al. 2024, arXiv:2405.13491\n- Euclid Collaboration: Scaramella, R., Amiaux, J., Mellier, Y., et al. 2022, A&A, 662, A112\n- Fumagalli, A., Costanzi, M., Saro, A., Castro, T., & Borgani, S. 2024, A&A, 682, A148  \nFumagalli, A., Saro, A., Borgani, S., et al. 2021, A&A, 652, A21  \n- Hamaus, N., Seljak, U., Desjacques, V., Smith, R. E., & Baldauf, T. 2010, PRD, D82, 043515  \nHasselfield, M., Hilton, M., Marriage, T. A., et al. 2013, JCAP, 07, 008  \nHo  \nff  \nman, M. D., Gelman, A., et al. 2014, J. Mach. Learn. Res., 15, 1593  \nHolder, G., Haiman, Z., & Mohr, J. 2001, ApJ Lett., 560, L111  \nHu, W. & Kravtsov, A. V. 2003, ApJ, 584, 702  \nHuterer, D. & Turner, M. S. 2001, PRD, 64, 123527  \nJenkins, A., Frenk, C. S., White, S. D. M., et al. 2001, MNRAS, 321, 372 Kitayama, T. & Suto, Y. 1996, ApJ, 469, 480  \nKnebe, A., Knollmann, S. R., Muldrew, S. I., et al. 2011, MNRAS, 415, 2293 Kravtsov, A. V. & Borgani, S. 2012, ARA&A, 50, 353  \n- Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, arXiv:1110.3193  \nLesci, G. F., Marulli, F., Moscardini, L., et al. 2022a, A&A, 659, A88  \nLesci, G. F., Nanni, L., Marulli, F., et al. 2022b, A&A, 665, A100  \nLoVerde, M. 2014, Phys. Rev. D, 90, 083518  \nMana, A., Giannantonio, T., Weller, J., et al. 2013, MNRAS, 434, 684  \nMantz, A. B., von der Linden, A., Allen, S. W., et al. 2015, MNRAS, 446, 2205  \nMichaux, M., Hahn, O., Rampf, C., & Angulo, R. E. 2021, MNRAS, 500, 663  \n- Mo, H. J. & White, S. D. M. 1996, MNRAS, 282, 347\n- Monaco, P., Sefusatti, E., Borgani, S., et al. 2013, MNRAS, 433, 2389\n- Monaco, P., Theuns, T., & Ta ff oni, G. 2002, MNRAS, 331, 587\n- Munari, E., Monaco, P., Sefusatti, E., et al. 2017, MNRAS, 465, 4658\n- Ondaro-Mallea, L., Angulo, R. E., Zennaro, M., Contreras, S., & Aricò, G. 2021, MNRAS, 509, 6077  \nPeebles, P. J. E. 2020, The large-scale structure of the universe, Vol. 98 (Princeton university press)  \nPlanck Collaboration XX: Ade, P. A. R., Aghanim, N., Armitage-Caplan, C., et al. 2014, A&A, 571, A20  \n- Planck Collaboration XXIV: Ade, P. A. R., Aghanim, N., Arnaud, M., et al. 2016,  \nA&A, 594, A24  \nPress, W. H. & Schechter, P. 1974, ApJ, 187, 425  \nRomanello, M., Marulli, F., Moscardini, L., et al. 2024, A&A, 682, A72  \nRozo, E., Wechsler, R. H., Ryko ff , E. S., et al. 2010, ApJ, 708, 645  \nSalvati, L., Douspis, M., & Aghanim, N. 2020, A&A, 643, A20  \nSalvatier, J., Wiecki, T. V., & Fonnesbeck, C. 2016, Probabilistic programming  \nin python using pymc3. PeerJ Computer Science, 2, e55  \nSaro, A., Bocquet, S., Rozo, E., et al. 2015, MNRAS, 454, 2305  \nSartoris, B., Biviano, A., Fedeli, C., et al. 2016, MNRAS, 459, 1764  \nSchneider, A. & Teyssier, R. 2015, JCAP, 12, 049  \nSheth, R. K., Mo, H. J., & Tormen, G. 2001, MNRAS, 323, 1  \nArticle number, page 18 of 20  \n- Springel, V., Pakmor, R., Zier, O., & Reinecke, M. 2021, MNRAS, 506, 2871 Springel, V., White, S. D. M., Tormen, G., & Kau ff mann, G. 2001, MNRAS, 328, 726\n- Sunayama, T., Miyatake, H., Sugiyama, S., More, S., et al. 2023, arXiv: 2309.13025\n- Ta ff oni, G., Becciani, U., Garilli, B., et al. 2020, in Astronomical Society of the Pacific Conference Series, Vol. 527, 307  \nTinker, J. L., Kravtsov, A. V., Klypin, A., et al. 2008, ApJ, 688, 709 Tinker, J. L., Robertson, B. E., Kravtsov, A. V., et al. 2010, ApJ, 724, 878 To, C., Krause, E., Rozo, E., et al. 2021, PRL, 126, 141301 Vehtari, A., Gelman, A., & Gabry, J. 2017, Statistics and Computing, 27, 1413 Velliscig, M., van Daalen, M. P., Schaye, J., et al. 2014, MNRAS, 442, 2641 Villaescusa-Navarro, F., Naess, S., Genel, S., et al. 2018, ApJ, 867, 137 Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261 Warren, M. S., Abazajian, K., Holz, D. E., & Teodoro, L. 2006, ApJ, 646, 881 Watanabe, S. 2010, Journal of Machine Learning Research, 11, 3571 Watson, W. A., Iliev, I. T., D\'Aloisio, A., et al. 2013, MNRAS, 433, 1230 Zhang, T., Chuang, C.-H., Wechsler, R. H., Alam, S., et al. 2022, MNRAS, 518,  \n- 3737\n- 1 INAF-Osservatorio Astronomico di Trieste, Via G. B. Tiepolo 11, 34143 Trieste, Italy\n- 2 INFN, Sezione di Trieste, Via Valerio 2, 34127 Trieste TS, Italy\n- 3 IFPU, Institute for Fundamental Physics of the Universe, via Beirut 2, 34151 Trieste, Italy\n- 4 ICSC - Centro Nazionale di Ricerca in High Performance Computing, Big Data e Quantum Computing, Via Magnanelli 2, Bologna, Italy\n- 5 Ludwig-Maximilians-University, Schellingstrasse 4, 80799 Munich, Germany\n- 6 Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal, 4, 20018, Donostia-San Sebastián, Guipuzkoa, Spain\n- 7 IKERBASQUE, Basque Foundation for Science, 48013, Bilbao, Spain\n- 8 Universitäts-Sternwarte München, Fakultät für Physik, LudwigMaximilians-Universität München, Scheinerstrasse 1, 81679 München, Germany\n- 9 Dipartimento di Fisica - Sezione di Astronomia, Università di Trieste, Via Tiepolo 11, 34131 Trieste, Italy\n- 10 Department of Astrophysics, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland\n- 11 Université Paris-Saclay, CNRS, Institut d\'astrophysique spatiale, 91405, Orsay, France\n- 12 Institut für Theoretische Physik, University of Heidelberg, Philosophenweg 16, 69120 Heidelberg, Germany\n- 13 INAF-Osservatorio Astronomico di Brera, Via Brera 28, 20122 Milano, Italy\n- 14 SISSA, International School for Advanced Studies, Via Bonomea 265, 34136 Trieste TS, Italy\n- 15 Dipartimento di Fisica e Astronomia, Università di Bologna, Via Gobetti 93 / 2, 40129 Bologna, Italy\n- 16 INAF-Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, Via Piero Gobetti 93 / 3, 40129 Bologna, Italy\n- 17 INFN-Sezione di Bologna, Viale Berti Pichat 6 / 2, 40127 Bologna, Italy\n- 18 Max Planck Institute for Extraterrestrial Physics, Giessenbachstr. 1, 85748 Garching, Germany\n- 19 INAF-Osservatorio Astrofisico di Torino, Via Osservatorio 20, 10025 Pino Torinese (TO), Italy\n- 20 Dipartimento di Fisica, Università di Genova, Via Dodecaneso 33, 16146, Genova, Italy\n- 21 INFN-Sezione di Genova, Via Dodecaneso 33, 16146, Genova, Italy\n- 22 Department of Physics "E. Pancini", University Federico II, Via Cinthia 6, 80126, Napoli, Italy\n- 23 INAF-Osservatorio Astronomico di Capodimonte, Via Moiariello 16, 80131 Napoli, Italy\n- 24 INFN section of Naples, Via Cinthia 6, 80126, Napoli, Italy\n- 25 Aix-Marseille Université, CNRS, CNES, LAM, Marseille, France\n- 26 Dipartimento di Fisica, Università degli Studi di Torino, Via P. Giuria 1, 10125 Torino, Italy  \n- 27 INFN-Sezione di Torino, Via P. Giuria 1, 10125 Torino, Italy\n- 28 INAF-IASF Milano, Via Alfonso Corti 12, 20133 Milano, Italy\n- 29 Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Avenida Complutense 40, 28040 Madrid, Spain\n- 30 Port d\'Informació Científica, Campus UAB, C. Albareda s / n, 08193 Bellaterra (Barcelona), Spain\n- 31 Institute for Theoretical Particle Physics and Cosmology (TTK), RWTH Aachen University, 52056 Aachen, Germany\n- 32 INAF-Osservatorio Astronomico di Roma, Via Frascati 33, 00078 Monteporzio Catone, Italy\n- 33 Dipartimento di Fisica e Astronomia "Augusto Righi" - Alma Mater Studiorum Università di Bologna, Viale Berti Pichat 6 / 2, 40127 Bologna, Italy\n- 34 Instituto de Astrofísica de Canarias, Calle Vía Láctea s / n, 38204, San Cristóbal de La Laguna, Tenerife, Spain\n- 35 Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK\n- 36 Jodrell Bank Centre for Astrophysics, Department of Physics and Astronomy, University of Manchester, Oxford Road, Manchester M13 9PL, UK\n- 37 European Space Agency / ESRIN, Largo Galileo Galilei 1, 00044 Frascati, Roma, Italy\n- 38 ESAC / ESA, Camino Bajo del Castillo, s / n., Urb. Villafranca del Castillo, 28692 Villanueva de la Cañada, Madrid, Spain\n- 39 Université Claude Bernard Lyon 1, CNRS / IN2P3, IP2I Lyon, UMR 5822, Villeurbanne, F-69100, France\n- 40 Institute of Physics, Laboratory of Astrophysics, Ecole Polytechnique Fédérale de Lausanne (EPFL), Observatoire de Sauverny, 1290 Versoix, Switzerland\n- 41 UCB Lyon 1, CNRS / IN2P3, IUF, IP2I Lyon, 4 rue Enrico Fermi, 69622 Villeurbanne, France\n- 42 Departamento de Física, Faculdade de Ciências, Universidade de Lisboa, Edifício C8, Campo Grande, PT1749-016 Lisboa, Portugal\n- 43 Instituto de Astrofísica e Ciências do Espaço, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal\n- 44 Department of Astronomy, University of Geneva, ch. d\'Ecogia 16, 1290 Versoix, Switzerland\n- 45 INAF-Istituto di Astrofisica e Planetologia Spaziali, via del Fosso del Cavaliere, 100, 00100 Roma, Italy\n- 46 INFN-Padova, Via Marzolo 8, 35131 Padova, Italy\n- 47 Université Paris-Saclay, Université Paris Cité, CEA, CNRS, AIM, 91191, Gif-sur-Yvette, France\n- 48 Institut d\'Estudis Espacials de Catalunya (IEEC), Edifici RDIT, Campus UPC, 08860 Castelldefels, Barcelona, Spain\n- 49 Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Can Magrans, s / n, 08193 Barcelona, Spain\n- 50 Istituto Nazionale di Fisica Nucleare, Sezione di Bologna, Via Irnerio 46, 40126 Bologna, Italy\n- 51 FRACTAL S.L.N.E., calle Tulipán 2, Portal 13 1A, 28231, Las Rozas de Madrid, Spain\n- 52 INAF-Osservatorio Astronomico di Padova, Via dell\'Osservatorio 5, 35122 Padova, Italy\n- 53 Dipartimento di Fisica "Aldo Pontremoli", Università degli Studi di Milano, Via Celoria 16, 20133 Milano, Italy\n- 54 Institute of Theoretical Astrophysics, University of Oslo, P.O. Box 1029 Blindern, 0315 Oslo, Norway\n- 55 Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA, 91109, USA\n- 56 Felix Hormuth Engineering, Goethestr. 17, 69181 Leimen, Germany\n- 57 Technical University of Denmark, Elektrovej 327, 2800 Kgs. Lyngby, Denmark\n- 58 Cosmic Dawn Center (DAWN), Denmark\n- 59 Université Paris-Saclay, CNRS / IN2P3, IJCLab, 91405 Orsay, France\n- 60 Institut de Recherche en Astrophysique et Planétologie (IRAP), Université de Toulouse, CNRS, UPS, CNES, 14 Av. Edouard Belin, 31400 Toulouse, France\n- 61 Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany\n- 62 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA  \n63  \nDepartment of Physics and Astronomy, University College Lon-  \ndon, Gower Street, London WC1E 6BT, UK  \n- 64 Department of Physics and Helsinki Institute of Physics, Gustaf Hällströmin katu 2, 00014 University of Helsinki, Finland\n- 65 Aix-Marseille Université, CNRS / IN2P3, CPPM, Marseille, France\n- 66 Université de Genève, Département de Physique Théorique and Centre for Astroparticle Physics, 24 quai Ernest-Ansermet, CH1211 Genève 4, Switzerland\n- 67 Department of Physics, P.O. Box 64, 00014 University of Helsinki, Finland\n- 68 Helsinki Institute of Physics, Gustaf Hällströmin katu 2, University of Helsinki, Helsinki, Finland\n- 69 NOVA optical infrared instrumentation group at ASTRON, Oude Hoogeveensedijk 4, 7991PD, Dwingeloo, The Netherlands\n- 70 Universität Bonn, Argelander-Institut für Astronomie, Auf dem Hügel 71, 53121 Bonn, Germany\n- 71 INFN-Sezione di Roma, Piazzale Aldo Moro, 2 - c / o Dipartimento di Fisica, Edificio G. Marconi, 00185 Roma, Italy\n- 72 Dipartimento di Fisica e Astronomia "Augusto Righi" - Alma Mater Studiorum Università di Bologna, via Piero Gobetti 93 / 2, 40129 Bologna, Italy\n- 73 Department of Physics, Institute for Computational Cosmology, Durham University, South Road, DH1 3LE, UK\n- 74 Université Côte d\'Azur, Observatoire de la Côte d\'Azur, CNRS, Laboratoire Lagrange, Bd de l\'Observatoire, CS 34229, 06304 Nice cedex 4, France\n- 75 University of Applied Sciences and Arts of Northwestern Switzerland, School of Engineering, 5210 Windisch, Switzerland\n- 76 Institut d\'Astrophysique de Paris, 98bis Boulevard Arago, 75014, Paris, France\n- 77 Institut d\'Astrophysique de Paris, UMR 7095, CNRS, and Sorbonne Université, 98 bis boulevard Arago, 75014 Paris, France\n- 78 European Space Agency / ESTEC, Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands\n- 79 Institut de Física d\'Altes Energies (IFAE), The Barcelona Institute of Science and Technology, Campus UAB, 08193 Bellaterra (Barcelona), Spain\n- 80 DARK, Niels Bohr Institute, University of Copenhagen, Jagtvej 155, 2200 Copenhagen, Denmark\n- 81 Waterloo Centre for Astrophysics, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada\n- 82 Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada\n- 83 Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada\n- 84 Space Science Data Center, Italian Space Agency, via del Politecnico snc, 00133 Roma, Italy\n- 85 Centre National d\'Etudes Spatiales - Centre spatial de Toulouse, 18 avenue Edouard Belin, 31401 Toulouse Cedex 9, France\n- 86 Institute of Space Science, Str. Atomistilor, nr. 409 M˘agurele, Ilfov, 077125, Romania\n- 87 Dipartimento di Fisica e Astronomia "G. Galilei", Università di Padova, Via Marzolo 8, 35131 Padova, Italy\n- 88 Université St Joseph; Faculty of Sciences, Beirut, Lebanon\n- 89 Departamento de Física, FCFM, Universidad de Chile, Blanco Encalada 2008, Santiago, Chile\n- 90 Satlantis, University Science Park, Sede Bld 48940, Leioa-Bilbao, Spain\n- 91 Instituto de Astrofísica e Ciências do Espaço, Faculdade de Ciências, Universidade de Lisboa, Tapada da Ajuda, 1349-018 Lisboa, Portugal\n- 92 Universidad Politécnica de Cartagena, Departamento de Electrónica y Tecnología de Computadoras, Plaza del Hospital 1, 30202 Cartagena, Spain\n- 93 INFN-Bologna, Via Irnerio 46, 40126 Bologna, Italy\n- 94 Kapteyn Astronomical Institute, University of Groningen, PO Box 800, 9700 AV Groningen, The Netherlands  \n- 95 Infrared Processing and Analysis Center, California Institute of Technology, Pasadena, CA 91125, USA\n- 96 INAF, Istituto di Radioastronomia, Via Piero Gobetti 101, 40129 Bologna, Italy\n- 97 Astronomical Observatory of the Autonomous Region of the Aosta Valley (OAVdA), Loc. Lignan 39, I-11020, Nus (Aosta Valley), Italy\n- 98 Junia, EPA department, 41 Bd Vauban, 59800 Lille, France\n- 99 Instituto de Física Teórica UAM-CSIC, Campus de Cantoblanco, 28049 Madrid, Spain\n- 100 CERCA / ISO, Department of Physics, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, OH 44106, USA\n- 101 Laboratoire Univers et Théorie, Observatoire de Paris, Université PSL, Université Paris Cité, CNRS, 92190 Meudon, France\n- 102 INFN-Sezione di Milano, Via Celoria 16, 20133 Milano, Italy\n- 103 Departamento de Física Fundamental. Universidad de Salamanca. Plaza de la Merced s / n. 37008 Salamanca, Spain\n- 104 Departamento de Astrofísica, Universidad de La Laguna, 38206, La Laguna, Tenerife, Spain\n- 105 Dipartimento di Fisica e Scienze della Terra, Università degli Studi di Ferrara, Via Giuseppe Saragat 1, 44122 Ferrara, Italy\n- 106 Istituto Nazionale di Fisica Nucleare, Sezione di Ferrara, Via Giuseppe Saragat 1, 44122 Ferrara, Italy\n- 107 Université de Strasbourg, CNRS, Observatoire astronomique de Strasbourg, UMR 7550, 67000 Strasbourg, France\n- 108 Center for Data-Driven Discovery, Kavli IPMU (WPI), UTIAS, The University of Tokyo, Kashiwa, Chiba 277-8583, Japan\n- 109 Max-Planck-Institut für Physik, Boltzmannstr. 8, 85748 Garching, Germany\n- 110 Minnesota Institute for Astrophysics, University of Minnesota, 116 Church St SE, Minneapolis, MN 55455, USA\n- 111 Institute Lorentz, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, The Netherlands\n- 112 Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822, USA\n- 113 Department of Physics & Astronomy, University of California Irvine, Irvine CA 92697, USA\n- 114 Department of Astronomy & Physics and Institute for Computational Astrophysics, Saint Mary\'s University, 923 Robie Street, Halifax, Nova Scotia, B3H 3C3, Canada\n- 115 Departamento Física Aplicada, Universidad Politécnica de Carta-\n- gena, Campus Muralla del Mar, 30202 Cartagena, Murcia, Spain\n- 116 Instituto de Astrofísica de Canarias (IAC); Departamento de Astrofísica, Universidad de La Laguna (ULL), 38200, La Laguna, Tenerife, Spain\n- 117 Department of Physics, Oxford University, Keble Road, Oxford OX1 3RH, UK\n- 118 Université Paris Cité, CNRS, Astroparticule et Cosmologie, 75013 Paris, France\n- 119 CEA Saclay, DFR / IRFU, Service d\'Astrophysique, Bat. 709, 91191 Gif-sur-Yvette, France\n- 120 Institute of Cosmology and Gravitation, University of Portsmouth, Portsmouth PO1 3FX, UK\n- 121 Department of Computer Science, Aalto University, PO Box 15400, Espoo, FI-00 076, Finland\n- 122 Instituto de Astrofísica de Canarias, c / Via Lactea s / n, La Laguna E-38200, Spain. Departamento de Astrofísica de la Universidad de La Laguna, Avda. Francisco Sanchez, La Laguna, E-38200, Spain\n- 123 Ruhr University Bochum, Faculty of Physics and Astronomy, Astronomical Institute (AIRUB), German Centre for Cosmological Lensing (GCCL), 44780 Bochum, Germany\n- 124 Univ. Grenoble Alpes, CNRS, Grenoble INP, LPSC-IN2P3, 53, Avenue des Martyrs, 38000, Grenoble, France\n- 125 Department of Physics and Astronomy, Vesilinnantie 5, 20014 University of Turku, Finland\n- 126 Serco for European Space Agency (ESA), Camino bajo del Castillo, s / n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain\n- 127 ARC Centre of Excellence for Dark Matter Particle Physics, Melbourne, Australia  \n- 128 Centre for Astrophysics & Supercomputing, Swinburne University of Technology, Hawthorn, Victoria 3122, Australia\n- 129 School of Physics and Astronomy, Queen Mary University of London, Mile End Road, London E1 4NS, UK\n- 130 Department of Physics and Astronomy, University of the Western Cape, Bellville, Cape Town, 7535, South Africa\n- 131 ICTP South American Institute for Fundamental Research, Instituto de Física Teórica, Universidade Estadual Paulista, São Paulo, Brazil\n- 132 IRFU, CEA, Université Paris-Saclay 91191 Gif-sur-Yvette Cedex, France\n- 133 Oskar Klein Centre for Cosmoparticle Physics, Department of Physics, Stockholm University, Stockholm, SE-106 91, Sweden\n- 134 Astrophysics Group, Blackett Laboratory, Imperial College London, London SW7 2AZ, UK\n- 135 INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125, Firenze, Italy\n- 136 Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 2, 00185 Roma, Italy\n- 137 Centro de Astrofísica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal\n- 138 Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, PT4150-762 Porto, Portugal\n- 139 HE Space for European Space Agency (ESA), Camino bajo del Castillo, s / n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain\n- 140 Aurora Technology for European Space Agency (ESA), Camino bajo del Castillo, s / n, Urbanizacion Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain\n- 141 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK\n- 142 Dipartimento di Fisica, Università degli studi di Genova, and INFN-Sezione di Genova, via Dodecaneso 33, 16146, Genova, Italy\n- 143 Theoretical astrophysics, Department of Physics and Astronomy, Uppsala University, Box 515, 751 20 Uppsala, Sweden\n- 144 Department of Physics, Royal Holloway, University of London, TW20 0EX, UK\n- 145 Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey RH5 6NT, UK\n- 146 Department of Astrophysical Sciences, Peyton Hall, Princeton University, Princeton, NJ 08544, USA\n- 147 Cosmic Dawn Center (DAWN)\n- 148 Niels Bohr Institute, University of Copenhagen, Jagtvej 128, 2200 Copenhagen, Denmark\n- 149 Center for Cosmology and Particle Physics, Department of Physics, New York University, New York, NY 10003, USA\n- 150 Center for Computational Astrophysics, Flatiron Institute, 162 5th Avenue, 10010, New York, NY, USA'}