Neutrinos and Galaxy Custering

**Chapter 12**

**Abstract**

Cosmologies

severely affect the constraints.

galaxy clustering

**1. Introduction**

**195**

*Jorge Enrique García-Farieta*

*and Rigoberto Ángel Casas Miranda*

Massive Neutrinos and Galaxy

Cosmic neutrinos have been playing a key role in cosmology since the discovery of their mass. They can affect cosmological observables and have several implications being the only hot dark matter candidates that we currently know to exist. The combination of massive neutrinos and an adequate theory of gravity provide a perfect scenario to address questions on the dark sector that have remained unanswered for years. In particular, in the era of precision cosmology, galaxy clustering and redshift-space distortions afford one of the most powerful tools to characterise the spatial distribution of cosmic tracers and to extract robust constraints on neutrino masses. In this chapter, we study how massive neutrinos affect the galaxy clustering and investigate whether the cosmological effects of massive neutrinos might be degenerate with *f*(*R*) gravity cosmologies, which would

**Keywords:** massive neutrinos, gravity theories, structure formation,

the formation and evolution of cosmic structures.

From first principles, it is well known that a theory of gravity is needed to describe the spatial properties and dynamics of the large-scale structures (LSS) of the universe. The observational data collected for several decades provide strong support to the concordance model Lambda cold dark matter (ΛCDM), which yields a consistent description of the main properties of the LSS [1–4]. However, since cosmological observations have entered in an unprecedented precision era, one of the current aims is to test some of the most fundamental assumptions of the concordance model of the universe. In this sense, the ΛCDM model assumes (i) the general theory of relativity (GR) as the theory describing gravitational interactions at large scales, (ii) the standard model of particles and (iii) the cosmological principle. Moreover, in this framework, the universe is currently dominated by dark energy (DE) in the form of a cosmological constant, responsible for the late-time cosmic acceleration [5–7] and by a cold dark matter (CDM) component that drives

Recently, several shortcomings have been found in the ΛCDM scenario, like a possible tension in the parameter constraint of *H*<sup>0</sup> and *σ*<sup>8</sup> when different probes are

Clustering in *f*(*R*) Gravity

## **Chapter 12**

## Massive Neutrinos and Galaxy Clustering in *f*(*R*) Gravity Cosmologies

*Jorge Enrique García-Farieta and Rigoberto Ángel Casas Miranda*

## **Abstract**

Cosmic neutrinos have been playing a key role in cosmology since the discovery of their mass. They can affect cosmological observables and have several implications being the only hot dark matter candidates that we currently know to exist. The combination of massive neutrinos and an adequate theory of gravity provide a perfect scenario to address questions on the dark sector that have remained unanswered for years. In particular, in the era of precision cosmology, galaxy clustering and redshift-space distortions afford one of the most powerful tools to characterise the spatial distribution of cosmic tracers and to extract robust constraints on neutrino masses. In this chapter, we study how massive neutrinos affect the galaxy clustering and investigate whether the cosmological effects of massive neutrinos might be degenerate with *f*(*R*) gravity cosmologies, which would severely affect the constraints.

**Keywords:** massive neutrinos, gravity theories, structure formation, galaxy clustering

## **1. Introduction**

From first principles, it is well known that a theory of gravity is needed to describe the spatial properties and dynamics of the large-scale structures (LSS) of the universe. The observational data collected for several decades provide strong support to the concordance model Lambda cold dark matter (ΛCDM), which yields a consistent description of the main properties of the LSS [1–4]. However, since cosmological observations have entered in an unprecedented precision era, one of the current aims is to test some of the most fundamental assumptions of the concordance model of the universe. In this sense, the ΛCDM model assumes (i) the general theory of relativity (GR) as the theory describing gravitational interactions at large scales, (ii) the standard model of particles and (iii) the cosmological principle. Moreover, in this framework, the universe is currently dominated by dark energy (DE) in the form of a cosmological constant, responsible for the late-time cosmic acceleration [5–7] and by a cold dark matter (CDM) component that drives the formation and evolution of cosmic structures.

Recently, several shortcomings have been found in the ΛCDM scenario, like a possible tension in the parameter constraint of *H*<sup>0</sup> and *σ*<sup>8</sup> when different probes are used [4, 8–10]. This has motivated the interest on theoretical models beyond GR. Among them, the models based on *f*(*R*) gravity are the favourite ones because of their generality and rich phenomenology [11, 12]. Moreover, modified gravity (MG) models represent one of the most viable alternatives to explain cosmic acceleration [13] that require satisfying simultaneously solar system constraints and to be consistent with the measured accelerated cosmic expansion and large-scale constraints [14–17]. An extra motivation to study MG models is given by the fact that massive neutrinos, the only (hot) dark matter candidates we actually know to exist, can affect these observables and have several cosmological implications [18]. However, the degeneracy between some MG models and the total neutrino mass [19–21] give rise to a limitation of many standard cosmological statistics [22, 23]. In this context, the clustering analysis and redshift-space galaxy clustering have been proven to be a powerful cosmological probe to discriminate among MG scenarios with massive neutrinos, as will be discussed in the following sections.

Regarding the dynamics of background universe in the standard framework, it is well described by the Friedmann-Lemaître-Robertson-Walker (FLRW) metric, whose line element in natural units *c* ¼ 1 is given by

$$ds^2 = -dt^2 + a^2(t) \left[ \frac{dr^2}{1 - kr^2} + r^2d\theta^2 + r^2 \sin^2\theta d\phi^2 \right].\tag{1}$$

where *a*(*t*) is the so-called scale factor, (*r*, *θ*, *ϕ*) are dimensionless sphericalpolar coordinates and *k* defines the geometry of the universe under consideration to be flat (*k* ¼ 0), open (*k*<0) or closed (*k*>0). The equations of motion that describe the time evolution of *a*(*t*) and the dynamic growth of the universe are called *Friedmann equations* and are given by

$$
\left(\frac{\dot{a}}{a}\right)^2 = \frac{8\pi G}{3}\rho - \frac{k}{a^2} + \frac{\Lambda}{3},\tag{2}
$$

assuming a constant density energy so that the corresponding EoS is *p* ¼ �*ρ*. The case of a flat universe (*k* ¼ 0) is interesting because of its agreement with many observational results [4]; it also implies a special value of the matter density in the universe that allows to introduce naturally a critical density *ρcr* in terms of the Hubble parameter. The critical density is also useful to define the dimensionless density parameter Ω*<sup>i</sup>* for the various species *i*, related to radiation Ω*r*, matter Ω*m*, cosmological constant ΩΛ and curvature Ω*k*. It is easy to verify that the sum of these parameters is equal to unity, as it can be expected from the Friedmann equations; in

*(a) The current composition of the universe derived from Planck data [4], (b) evolution of the energy density*

fact, Ω*<sup>m</sup>* þ Ω*<sup>r</sup>* þ ΩΛ þ Ω*<sup>k</sup>* ¼ 1 is known as the *cosmic sum rule*. Therefore, a Friedmann universe can be described by the cosmological parameters

given by

**Figure 1.**

*H*2

universe, see Refs. [25–27].

**2.1 Modified gravity models**

**197**

ð Þ¼ *<sup>a</sup> <sup>H</sup>*<sup>2</sup>

*parameter ρ normalized to ρcr*,0 *as a function of the scale factor* a*.*

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies*

*DOI: http://dx.doi.org/10.5772/intechopen.92205*

**2. Modified gravity and massive neutrinos**

*H*0, Ω*m*, Ω*<sup>r</sup>* ð Þ , Ω*k*, ΩΛ such that the expansion rate as a function of the scale factor is

This equation is usually written as a dimensionless function defined by *E a*ð Þ¼ *H a*ð Þ*=H*0. The last constraints on cosmological parameters obtained by the Planck satellite show that *w* ¼ �1*:*03 � 0*:*03, consistent with a cosmological constant, and Ω*<sup>m</sup>* ¼ 0*:*3153 � 0*:*0073, ΩΛ ¼ 0*:*692 � 0*:*012, Ω*<sup>k</sup>* ¼ 0*:*000 � 0*:*005, where Ω*<sup>m</sup>* contains the density of baryons ð Þ Ω*<sup>b</sup>* and cold dark matter ð Þ Ω*CDM* . Additionally, in the last few years, it became usual to include the energy density of neutrinos Ω*ν*; they contribute to the radiation density at early times but behave as matter after the nonrelativistic transition at late times [24], so that for a flat universe, the total energy density is given by *ρ* ¼ *ργ* þ *ρ*CDM þ *ρ*<sup>b</sup> þ *ρν* þ *ρ*Λ. **Figure 1a** shows the percentages, derived from [4] data, in which each species contributes to the total content of the universe. **Figure 1b** shows the evolution of the density parameter *ρ* (in units of *ρcr*,0) as a function of the scale factor *a*. For further details concerning the background

One of the most interesting modifications of GR is that which modifies the

Einstein-Hilbert action by introducing a scalar function, *f*(*R*), as follows:

<sup>0</sup> <sup>Ω</sup>*ra*�<sup>4</sup> <sup>þ</sup> <sup>Ω</sup>*ma*�<sup>3</sup> <sup>þ</sup> <sup>Ω</sup>*ka*�<sup>2</sup> <sup>þ</sup> ΩΛ

*:* (6)

$$2\frac{\ddot{a}}{a} + \left(\frac{\dot{a}}{a}\right)^2 = -8\pi G p - \frac{k}{a^2} + \Lambda,\tag{3}$$

that can be re-expressed as

$$
\dot{\rho} = -3\frac{\dot{a}}{a}(\rho + p),
\tag{4}
$$

$$\frac{\ddot{a}}{a} = -\frac{4\pi G}{3}(\rho + 3p) + \frac{\Lambda}{3} \tag{5}$$

after eliminating the curvature term. These equations lead to the definition of the Hubble parameter as *<sup>H</sup>* � *<sup>a</sup>*\_ *<sup>a</sup>*, which drives the expansion rate of the universe and usually is represented in terms of the dimensionless factor *h*, defined by the expression H0 = 100 h km s�<sup>1</sup> Mp c�<sup>1</sup> at the present epoch. In order to have a full description of the background universe, it is necessary an equation of state (EoS) of the cosmic fluid, considering that it has three principal components: baryonic matter, dark matter and radiation (see **Figure 1a**). A first approximation consists in assuming a linear relationship between *ρ* and *p*; thus, the equation of state can be written as *p* ¼ *wρ*, where *w* is a parameter that in principle can be timedependent *w* ¼ *w t*ð Þ, but the simplest approach is to consider it as a constant. Under this assumption, Eq. (4) is easy to solve, resulting in *<sup>ρ</sup>*ð Þ¼ *<sup>a</sup> <sup>a</sup>*�3 1ð Þ <sup>þ</sup>*<sup>w</sup>* . The case with *w* ¼ �1 corresponds to the so-called cosmological constant; it is obtained *Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies DOI: http://dx.doi.org/10.5772/intechopen.92205*

**Figure 1.**

used [4, 8–10]. This has motivated the interest on theoretical models beyond GR. Among them, the models based on *f*(*R*) gravity are the favourite ones because of their generality and rich phenomenology [11, 12]. Moreover, modified gravity (MG) models represent one of the most viable alternatives to explain cosmic acceleration [13] that require satisfying simultaneously solar system constraints and to be consistent with the measured accelerated cosmic expansion and large-scale constraints [14–17]. An extra motivation to study MG models is given by the fact that massive neutrinos, the only (hot) dark matter candidates we actually know to exist, can affect these observables and have several cosmological implications [18]. However, the degeneracy between some MG models and the total neutrino mass [19–21] give rise to a limitation of many standard cosmological statistics [22, 23]. In this context, the clustering analysis and redshift-space galaxy clustering have been proven to be a powerful cosmological probe to discriminate among MG scenarios

with massive neutrinos, as will be discussed in the following sections.

ð Þ*t*

*a*\_ *a* � �<sup>2</sup>

*a*€

2 *a*€ *a* þ *a*\_ *a* � �<sup>2</sup>

whose line element in natural units *c* ¼ 1 is given by

*ds*<sup>2</sup> ¼ �*dt*<sup>2</sup> <sup>þ</sup> *<sup>a</sup>*<sup>2</sup>

*Friedmann equations* and are given by

*Progress in Fine Particle Plasmas*

that can be re-expressed as

the Hubble parameter as *<sup>H</sup>* � *<sup>a</sup>*\_

**196**

Regarding the dynamics of background universe in the standard framework, it is

2 *<sup>d</sup>θ*<sup>2</sup> <sup>þ</sup> *<sup>r</sup>*

" #

<sup>2</sup> sin <sup>2</sup>

*θdϕ*<sup>2</sup>

<sup>3</sup> , (2)

*<sup>a</sup>*<sup>2</sup> <sup>þ</sup> <sup>Λ</sup>, (3)

<sup>3</sup> (5)

*<sup>a</sup>* ð Þ *<sup>ρ</sup>* <sup>þ</sup> *<sup>p</sup>* , (4)

*:* (1)

well described by the Friedmann-Lemaître-Robertson-Walker (FLRW) metric,

*dr*<sup>2</sup> <sup>1</sup> � *kr*<sup>2</sup> <sup>þ</sup> *<sup>r</sup>*

<sup>¼</sup> <sup>8</sup>*π<sup>G</sup>*

*ρ*\_ ¼ �3

*<sup>a</sup>* ¼ � <sup>4</sup>*π<sup>G</sup>*

*a*\_

<sup>3</sup> *<sup>ρ</sup>* � *<sup>k</sup>*

¼ �8*πGp* � *<sup>k</sup>*

<sup>3</sup> ð Þþ *<sup>ρ</sup>* <sup>þ</sup> <sup>3</sup>*<sup>p</sup>*

after eliminating the curvature term. These equations lead to the definition of

and usually is represented in terms of the dimensionless factor *h*, defined by the expression H0 = 100 h km s�<sup>1</sup> Mp c�<sup>1</sup> at the present epoch. In order to have a full description of the background universe, it is necessary an equation of state (EoS) of the cosmic fluid, considering that it has three principal components: baryonic matter, dark matter and radiation (see **Figure 1a**). A first approximation consists in assuming a linear relationship between *ρ* and *p*; thus, the equation of state can be written as *p* ¼ *wρ*, where *w* is a parameter that in principle can be timedependent *w* ¼ *w t*ð Þ, but the simplest approach is to consider it as a constant. Under this assumption, Eq. (4) is easy to solve, resulting in *<sup>ρ</sup>*ð Þ¼ *<sup>a</sup> <sup>a</sup>*�3 1ð Þ <sup>þ</sup>*<sup>w</sup>* . The case with *w* ¼ �1 corresponds to the so-called cosmological constant; it is obtained

Λ

*<sup>a</sup>*, which drives the expansion rate of the universe

*a*<sup>2</sup> þ Λ

where *a*(*t*) is the so-called scale factor, (*r*, *θ*, *ϕ*) are dimensionless sphericalpolar coordinates and *k* defines the geometry of the universe under consideration to be flat (*k* ¼ 0), open (*k*<0) or closed (*k*>0). The equations of motion that describe the time evolution of *a*(*t*) and the dynamic growth of the universe are called

*(a) The current composition of the universe derived from Planck data [4], (b) evolution of the energy density parameter ρ normalized to ρcr*,0 *as a function of the scale factor* a*.*

assuming a constant density energy so that the corresponding EoS is *p* ¼ �*ρ*. The case of a flat universe (*k* ¼ 0) is interesting because of its agreement with many observational results [4]; it also implies a special value of the matter density in the universe that allows to introduce naturally a critical density *ρcr* in terms of the Hubble parameter. The critical density is also useful to define the dimensionless density parameter Ω*<sup>i</sup>* for the various species *i*, related to radiation Ω*r*, matter Ω*m*, cosmological constant ΩΛ and curvature Ω*k*. It is easy to verify that the sum of these parameters is equal to unity, as it can be expected from the Friedmann equations; in fact, Ω*<sup>m</sup>* þ Ω*<sup>r</sup>* þ ΩΛ þ Ω*<sup>k</sup>* ¼ 1 is known as the *cosmic sum rule*. Therefore, a Friedmann universe can be described by the cosmological parameters *H*0, Ω*m*, Ω*<sup>r</sup>* ð Þ , Ω*k*, ΩΛ such that the expansion rate as a function of the scale factor is given by

$$H^2(a) = H\_0^2 \left[\Omega\_r a^{-4} + \Omega\_m a^{-3} + \Omega\_k a^{-2} + \Omega\_\Lambda\right]. \tag{6}$$

This equation is usually written as a dimensionless function defined by *E a*ð Þ¼ *H a*ð Þ*=H*0. The last constraints on cosmological parameters obtained by the Planck satellite show that *w* ¼ �1*:*03 � 0*:*03, consistent with a cosmological constant, and Ω*<sup>m</sup>* ¼ 0*:*3153 � 0*:*0073, ΩΛ ¼ 0*:*692 � 0*:*012, Ω*<sup>k</sup>* ¼ 0*:*000 � 0*:*005, where Ω*<sup>m</sup>* contains the density of baryons ð Þ Ω*<sup>b</sup>* and cold dark matter ð Þ Ω*CDM* . Additionally, in the last few years, it became usual to include the energy density of neutrinos Ω*ν*; they contribute to the radiation density at early times but behave as matter after the nonrelativistic transition at late times [24], so that for a flat universe, the total energy density is given by *ρ* ¼ *ργ* þ *ρ*CDM þ *ρ*<sup>b</sup> þ *ρν* þ *ρ*Λ. **Figure 1a** shows the percentages, derived from [4] data, in which each species contributes to the total content of the universe. **Figure 1b** shows the evolution of the density parameter *ρ* (in units of *ρcr*,0) as a function of the scale factor *a*. For further details concerning the background universe, see Refs. [25–27].

## **2. Modified gravity and massive neutrinos**

### **2.1 Modified gravity models**

One of the most interesting modifications of GR is that which modifies the Einstein-Hilbert action by introducing a scalar function, *f*(*R*), as follows:

$$\mathbf{S} = \int \mathbf{d}^4 \mathbf{x} \sqrt{-\mathbf{g}} \left( \frac{\mathbf{R} + f(\mathbf{R})}{16\pi G} + \mathcal{L}\_m \right),\tag{7}$$

specific combinations of the *f*(*R*) function and of the total neutrino mass *m<sup>ν</sup>* � Σ*m<sup>ν</sup>*,*i*, standard statistics would not distinguish them from the standard ΛCDM expectations (see [21, 22, 36]). In addition, since the degeneracy is mostly driven by the non-linear behaviour of both the MG and the massive neutrinos effects on the LSS, the linear tools are not suitable to properly disentangle the combined parameter space [36].

Motivated by the apparent violation of energy, momentum and spin in *β*-decay processes, Pauli proposed the existence of neutrinos in 1930 to keep the conservation laws safe. Eventually, 26 years after they have been theoretically postulated, the neutrinos were detected for the first time by Cowan [37]. Neutrinos are classified in three 'flavours' in the standard model of particles; they were considered to be massless for some time until the discovery of the neutrino oscillation phenomena, i.e. related to the change of flavour [38]. Since then, it is known that at least two of the three neutrino families are massive, in contrast to the particle standard model assumption; however, measuring the absolute masses of the neutrinos is not easy, which makes this a very active field of research today, both for cosmology and particle physics. In a cosmological context, the neutrinos leave detectable imprints on observations that can then be used to constrain their properties; in particular, the presence of massive neutrinos impacts the background evolution of the universe and the growth of structures [39]. In the early universe, massive neutrinos are relativistic and indistinguishable from the massless ones, behaving like photons, meaning that their energy density drops like ∝*a*�4. In this stage, neutrinos are in thermal equilibrium, and their momentum follows the standard Fermi-Dirac

*n p*ð Þd*<sup>p</sup>* <sup>¼</sup> <sup>4</sup>*πg<sup>ν</sup>*

ð Þ <sup>2</sup>*π*ℏ*<sup>c</sup>* <sup>3</sup>

where *n p*ð Þ is the number of cosmic neutrinos with momentum between *p* and *p* þ *dp*, *g<sup>ν</sup>* is the number of neutrino spin states, *Tν*ð Þ*z* is the neutrino temperature at redshift *z* and *kB* is the Boltzmann constant. In principle, in the momentum distribution function, the chemical potential should be also included; however, it has been shown to be negligible for cosmological neutrinos [40]. The temperature of the cosmic neutrino background and the one from CMB are related by *Tν*ð Þ¼ *z* ¼ 0

*T<sup>γ</sup>* ð Þ *z* ¼ 0 [41], such that the temperature of the neutrino background at

After neutrinos with mass *m<sup>ν</sup>* decouple from the rest of the plasma at redshift *znr*,

the number density per flavour is fixed by the temperature, so that the universe is currently filled by a relic neutrino background, uniformly distributed, with a density of 113 part/cm<sup>3</sup> per species and an average temperature of 1.95 K. As

1eV

, (14)

<sup>1</sup> <sup>þ</sup> *znr*ð Þ *<sup>m</sup><sup>ν</sup>* <sup>≃</sup> <sup>1890</sup> *<sup>m</sup><sup>ν</sup>*

certain redshift *z* is given by *Tν*ð Þffi *z* 1*:*95 1ð Þ þ *z* K. Then, when the average momentum of neutrinos drops below a certain mass, they become non-relativistic, and their energy density drops like ∝*a*�3, behaving like baryons and cold dark matter. **Figure 2** shows the evolution of the massive neutrino density, normalised to the today's critical density, as a function of the scale factor from its early stage to the

*p*2d*p*

, (13)

*e p kBT<sup>ν</sup>* þ 1

**2.2 Massive neutrinos and the large-scale structure**

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies*

*DOI: http://dx.doi.org/10.5772/intechopen.92205*

distribution

4 <sup>11</sup> <sup>1</sup>*=*<sup>3</sup>

**199**

late universe.

as shown by Eq. (14)

where *R* is the Ricci scalar, *G* is the Newton's gravitational constant, *g* is the determinant of the metric tensor *gμν* and L*<sup>m</sup>* is the Lagrangian density of all matter fields. For a classification of *f*(*R*) theories of gravity and the assumptions needed to arrive to the various versions of *f*(*R*) gravity and GR, see, e.g., [28]. Thus, for a general *f*(*R*) model, one can consider a spatially flat FLRW universe with metric *ds*<sup>2</sup> ¼ �*dt*<sup>2</sup> <sup>þ</sup> *<sup>a</sup>*2ð Þ*<sup>t</sup> <sup>d</sup>***x**2, so that varying the Einstein-Hilbert action with respect to *<sup>g</sup>μν*, one can get a general form of the modified Einstein field equations. Consequently, the corresponding modified Friedmann equations are given by

$$2\pi \mathcal{H}^2(\mathbb{1} + f\_R) = 8\pi \mathcal{G}(\rho\_m + \rho\_{rad}) + \frac{1}{2} \left(f\_R R - f\right) - 3H\dot{f}\_R \tag{8}$$

$$-2\left(\mathbf{1} + f\_R\right)\dot{H} = 8\pi\mathbf{G}\left(\rho\_m + \frac{4}{3}\rho\_{rad}\right) + \ddot{f}\_R - H\dot{f}\_R\,\,,\tag{9}$$

where *<sup>R</sup>* <sup>¼</sup> 6 2*H*<sup>2</sup> <sup>þ</sup> *<sup>H</sup>*\_ � � and the over-dot denotes a derivative with respect to the cosmic time *t*. In general, the background evolution of a viable *f*(*R*) is not simple as it has been shown by [29–31]. However it is possible to get an approximation in a way that is analogous to the DE models, by neglecting the higher derivative and the non-linear terms. By defining the growth rate as *f* � *d* ln *δ=d* ln *a*, the equation that describes the growth of matter perturbations in terms of the density contrast *δ* in a *f* (*R*) model is approximated by [29, 32]

$$-\mathfrak{A}\mathfrak{\Omega}\_{m}(a)(\mathfrak{1}-\mathfrak{\Omega}\_{m}(a))\frac{df}{d\mathfrak{\Omega}\_{m}(a)} + f^{2} + \left(2 - \frac{3}{2}\mathfrak{\Omega}\_{m}(a)\right)f = \frac{3}{2}\frac{G\_{\text{eff}}}{G}\mathfrak{\Omega}\_{m}(a). \tag{10}$$

where *Geff* is the effective gravitational constant that can be written as [33]. A plausible *f*(*R*) function able to satisfy the solar system constraints, to mimic the ΛCDM model at high-redshift regime where it is well tested by the CMB and, at the same time, to accelerate the expansion of the universe at low redshift but without a cosmological constant [34], suggests that

$$\lim\_{R \to \infty} f(R) = \text{const}, \quad \lim\_{R \to 0} f(R) = 0,\tag{11}$$

which can be satisfied by a broken power law function such that

$$f(R) = -m^2 \frac{c\_1 \left(\frac{R}{m^2}\right)^n}{c\_2 \left(\frac{R}{m^2}\right)^n + 1},\tag{12}$$

where the mass scale *<sup>m</sup>* is defined as *<sup>m</sup>*<sup>2</sup> � *<sup>H</sup>*<sup>2</sup> <sup>0</sup>Ω<sup>M</sup> and *c*1, *c*<sup>2</sup> and *n* are nonnegative free parameters of the model [34]. For this *f*(*R*) model, the background expansion history is consistent with the ΛCDM case by choosing *c*1*=c*<sup>2</sup> ¼ 6ΩΛ*=*ΩM, where ΩΛ and Ω<sup>M</sup> are the dimensionless density parameters for vacuum and matter, respectively.

Nowadays it is generally accepted that some MG theories such as the [34] *f*(*R*) are strongly degenerated in a wide range of their observables with the effects of massive neutrinos; see, e.g., [19–21, 35]. This represents a serious challenge constraining cosmological models from current and future galaxy surveys requiring robust and reliable methods to disentangle both phenomena. Furthermore, for some specific combinations of the *f*(*R*) function and of the total neutrino mass *m<sup>ν</sup>* � Σ*m<sup>ν</sup>*,*i*, standard statistics would not distinguish them from the standard ΛCDM expectations (see [21, 22, 36]). In addition, since the degeneracy is mostly driven by the non-linear behaviour of both the MG and the massive neutrinos effects on the LSS, the linear tools are not suitable to properly disentangle the combined parameter space [36].

### **2.2 Massive neutrinos and the large-scale structure**

*S* ¼ ð d4 *x* ffiffiffiffiffiffi

*Progress in Fine Particle Plasmas*

corresponding modified Friedmann equations are given by

*d*Ω*m*ð Þ *a*

lim *R*!∞

where the mass scale *<sup>m</sup>* is defined as *<sup>m</sup>*<sup>2</sup> � *<sup>H</sup>*<sup>2</sup>

þ *f*

� � <sup>¼</sup> <sup>8</sup>*π<sup>G</sup> <sup>ρ</sup><sup>m</sup>* <sup>þ</sup> *<sup>ρ</sup>rad* ð Þþ <sup>1</sup>

� �*H*\_ <sup>¼</sup> <sup>8</sup>*π<sup>G</sup> <sup>ρ</sup><sup>m</sup>* <sup>þ</sup>

<sup>3</sup>*H*<sup>2</sup> <sup>1</sup> <sup>þ</sup> *<sup>f</sup> <sup>R</sup>*

(*R*) model is approximated by [29, 32]

cosmological constant [34], suggests that

respectively.

**198**

�3Ω*m*ð Þ *<sup>a</sup>* ð Þ <sup>1</sup> � <sup>Ω</sup>*m*ð Þ *<sup>a</sup> df*

�2 1 þ *f <sup>R</sup>*

�*<sup>g</sup>* <sup>p</sup> *<sup>R</sup>* <sup>þ</sup> *f R*ð Þ

where *R* is the Ricci scalar, *G* is the Newton's gravitational constant, *g* is the determinant of the metric tensor *gμν* and L*<sup>m</sup>* is the Lagrangian density of all matter fields. For a classification of *f*(*R*) theories of gravity and the assumptions needed to arrive to the various versions of *f*(*R*) gravity and GR, see, e.g., [28]. Thus, for a general *f*(*R*) model, one can consider a spatially flat FLRW universe with metric *ds*<sup>2</sup> ¼ �*dt*<sup>2</sup> <sup>þ</sup> *<sup>a</sup>*2ð Þ*<sup>t</sup> <sup>d</sup>***x**2, so that varying the Einstein-Hilbert action with respect to *<sup>g</sup>μν*, one can get a general form of the modified Einstein field equations. Consequently, the

16*πG* þ L*<sup>m</sup>* � �

4 <sup>3</sup> *<sup>ρ</sup>rad* � �

where *<sup>R</sup>* <sup>¼</sup> 6 2*H*<sup>2</sup> <sup>þ</sup> *<sup>H</sup>*\_ � � and the over-dot denotes a derivative with respect to the cosmic time *t*. In general, the background evolution of a viable *f*(*R*) is not simple as it has been shown by [29–31]. However it is possible to get an approximation in a way that is analogous to the DE models, by neglecting the higher derivative and the non-linear terms. By defining the growth rate as *f* � *d* ln *δ=d* ln *a*, the equation that describes the growth of matter perturbations in terms of the density contrast *δ* in a *f*

<sup>2</sup> <sup>þ</sup> <sup>2</sup> � <sup>3</sup>

where *Geff* is the effective gravitational constant that can be written as [33]. A plausible *f*(*R*) function able to satisfy the solar system constraints, to mimic the ΛCDM model at high-redshift regime where it is well tested by the CMB and, at the same time, to accelerate the expansion of the universe at low redshift but without a

*f R*ð Þ¼ const, lim

which can be satisfied by a broken power law function such that

*f R*ð Þ¼�*m*<sup>2</sup> *<sup>c</sup>*<sup>1</sup>

*c*2 *R m*<sup>2</sup> � �*<sup>n</sup>* <sup>þ</sup> <sup>1</sup>

Nowadays it is generally accepted that some MG theories such as the [34] *f*(*R*) are strongly degenerated in a wide range of their observables with the effects of massive neutrinos; see, e.g., [19–21, 35]. This represents a serious challenge

constraining cosmological models from current and future galaxy surveys requiring robust and reliable methods to disentangle both phenomena. Furthermore, for some

negative free parameters of the model [34]. For this *f*(*R*) model, the background expansion history is consistent with the ΛCDM case by choosing *c*1*=c*<sup>2</sup> ¼ 6ΩΛ*=*ΩM, where ΩΛ and Ω<sup>M</sup> are the dimensionless density parameters for vacuum and matter,

2 Ω*m*ð Þ *a* � �

*R*!0

*R m*<sup>2</sup> � �*<sup>n</sup>*

<sup>2</sup> *<sup>f</sup> <sup>R</sup><sup>R</sup>* � *<sup>f</sup>* � � � <sup>3</sup>*H*\_

<sup>þ</sup> €*<sup>f</sup> <sup>R</sup>* � *<sup>H</sup>*\_

*<sup>f</sup>* <sup>¼</sup> <sup>3</sup> 2 *G*eff

, (7)

*f <sup>R</sup>* (8)

*f <sup>R</sup>* , (9)

*<sup>G</sup>* <sup>Ω</sup>*m*ð Þ *<sup>a</sup> :* (10)

*f R*ð Þ¼ 0, (11)

, (12)

<sup>0</sup>Ω<sup>M</sup> and *c*1, *c*<sup>2</sup> and *n* are non-

Motivated by the apparent violation of energy, momentum and spin in *β*-decay processes, Pauli proposed the existence of neutrinos in 1930 to keep the conservation laws safe. Eventually, 26 years after they have been theoretically postulated, the neutrinos were detected for the first time by Cowan [37]. Neutrinos are classified in three 'flavours' in the standard model of particles; they were considered to be massless for some time until the discovery of the neutrino oscillation phenomena, i.e. related to the change of flavour [38]. Since then, it is known that at least two of the three neutrino families are massive, in contrast to the particle standard model assumption; however, measuring the absolute masses of the neutrinos is not easy, which makes this a very active field of research today, both for cosmology and particle physics. In a cosmological context, the neutrinos leave detectable imprints on observations that can then be used to constrain their properties; in particular, the presence of massive neutrinos impacts the background evolution of the universe and the growth of structures [39]. In the early universe, massive neutrinos are relativistic and indistinguishable from the massless ones, behaving like photons, meaning that their energy density drops like ∝*a*�4. In this stage, neutrinos are in thermal equilibrium, and their momentum follows the standard Fermi-Dirac distribution

$$n(p)\mathrm{d}p = \frac{4\pi\mathrm{g}\_{\nu}}{\left(2\pi\hbar c\right)^{3}}\frac{p^{2}\mathrm{d}p}{\mathrm{e}^{\mathrm{k}\_{B}\mathrm{T}}+\mathbf{1}},\tag{13}$$

where *n p*ð Þ is the number of cosmic neutrinos with momentum between *p* and *p* þ *dp*, *g<sup>ν</sup>* is the number of neutrino spin states, *Tν*ð Þ*z* is the neutrino temperature at redshift *z* and *kB* is the Boltzmann constant. In principle, in the momentum distribution function, the chemical potential should be also included; however, it has been shown to be negligible for cosmological neutrinos [40]. The temperature of the cosmic neutrino background and the one from CMB are related by *Tν*ð Þ¼ *z* ¼ 0 4 <sup>11</sup> <sup>1</sup>*=*<sup>3</sup> *T<sup>γ</sup>* ð Þ *z* ¼ 0 [41], such that the temperature of the neutrino background at certain redshift *z* is given by *Tν*ð Þffi *z* 1*:*95 1ð Þ þ *z* K. Then, when the average momentum of neutrinos drops below a certain mass, they become non-relativistic, and their energy density drops like ∝*a*�3, behaving like baryons and cold dark matter. **Figure 2** shows the evolution of the massive neutrino density, normalised to the today's critical density, as a function of the scale factor from its early stage to the late universe.

After neutrinos with mass *m<sup>ν</sup>* decouple from the rest of the plasma at redshift *znr*, as shown by Eq. (14)

$$1 + z\_{nr}(m\_{\nu}) \simeq 1890 \left(\frac{m\_{\nu}}{1 \text{eV}}\right),\tag{14}$$

the number density per flavour is fixed by the temperature, so that the universe is currently filled by a relic neutrino background, uniformly distributed, with a density of 113 part/cm<sup>3</sup> per species and an average temperature of 1.95 K. As

In this context, the formation and evolution of cosmic structures can be understood as a dynamical system of many particles, which trace the underlying mass distribution in a certain cosmological model. The N-body simulations, methods and algorithms have progressed continuously, achieving a high resolution to resolve finer structures with millions of particles, reducing the gap between theory and observations. For a detailed description on fundamentals of cosmological simulations, see,

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies*

*DOI: http://dx.doi.org/10.5772/intechopen.92205*

Since the formation and evolution of cosmic structures is based on the growth of

small fluctuations in the density field, it is expected that the amplitude of these initial perturbations have the correct value at late times to match the observed clustering today. An analytical development based on perturbation theory makes possible to follow the growth of structures to a certain extent using the linear approximation, being valid as long as *δρ* ≪ *ρ*. Nevertheless, these calculations are limited and cannot be extrapolated to explain completely the observational data; they break down on a scale where the density contrast *δ* � 1. Moreover, beyond the linear regime, the observed structures have a density contrast in a wide range from cosmic voids with *<sup>δ</sup>* � �1 to *<sup>δ</sup>* � <sup>10</sup><sup>6</sup> and larger. It makes necessary a more elaborated description of the perturbations in the non-linear regime, which can be achieved using higher-order perturbation theory or numerical simulations [48]. In this section we show a complementary analysis to the one performed in [52] in order to investigate the clustering in the context of modified gravity with massive neutrinos. We used a subset of the DUSTGRAIN-*pathfinder* runs [36], which implement the Hu-Sawicki *f*(*R*) model including massive neutrinos and whose cosmological parameters are consistent with Planck 2015 constraints [53].

As CDM haloes form from collapsing regions that detach from the background density field, their abundance can be related to the volume fraction of a Gaussian density smoothed on a radius *R* above a critical collapse threshold *δ<sup>c</sup>* [54]. The comoving number density of the haloes is strictly related to underlining cosmological model, such

> *ρ M*

*ρ*. The simplest argument to compute analytically the multiplicity func-

*<sup>M</sup>*, resulting in

*c* 2*σ*<sup>2</sup>ð Þ *M*

*<sup>σ</sup>*ð Þ *<sup>M</sup>* exp � *<sup>δ</sup>*<sup>2</sup>

Another approach to determine *f*ð Þ *σ* is given by accurate fitting functions, like the proposed by [55], which extends phenomenologically the results of [54]. For [55] the function *f*ð Þ *σ* is expected to be universal to the changes in redshift and

where *f*ð Þ *σ* is the multiplicity function, *σ* the RMS variance of the linear density

*f*ð Þ *σ*ð Þ *M*, *z ρ* quantifies the amount of mass contained in fluctuations of typical mass

tion *f*ð Þ *σ* comes from the spherical collapse theory, following the [54] formalism, such that a perturbation is supposed to collapse when it reaches the threshold *δ<sup>c</sup>* ≃1*:*68, by assuming that the probability distribution for a perturbation on a scale

d ln *<sup>σ</sup>*ð Þ *<sup>M</sup>*, *<sup>z</sup>* �<sup>1</sup>

<sup>d</sup>*<sup>M</sup>* , (17)

� �*:* (18)

that within a mass interval ½ � *M*, *M* þ *dM* , the halo mass function is given by

field smoothed on scale *R*(*M*) and *ρ* the mean matter density. The product

ffiffiffi 2 *π* r *δ<sup>c</sup>*

<sup>d</sup>*<sup>M</sup>* <sup>¼</sup> *<sup>f</sup>*ð Þ *<sup>σ</sup>*ð Þ *<sup>M</sup>*, *<sup>z</sup>*

d*n M*ð Þ , *z*

*M* is a Gaussian function with variance *σ*<sup>2</sup>

cosmology and is parameterized as follows:

*f*ð Þ¼ *σ*ð Þ *M*

e.g., [48–51].

**3.1 Halo mass function**

*<sup>M</sup>* <sup>¼</sup> <sup>4</sup> <sup>3</sup> *<sup>π</sup>R*<sup>3</sup>

**201**

#### **Figure 2.**

*Evolution of the massive neutrino density, normalised with respect today's critical density, as a function of the scale factor* a*, for a neutrino mass* m? *= 0.1 eV. The pink solid line shows the correct solution for ρν*ð Þ *a , while the green and blue dashed lines show the limit at early and late time, respectively. The early time limit is proportional to* a*<sup>4</sup> , underlining that at early times neutrinos are relativistic and behave like radiation. The late time limit, instead, is proportional to* a*<sup>3</sup> , following the same evolution of matter. The short vertical dashed line (black) marks the non-relativistic transition for neutrinos of this given mass. The figure was computed by using the CAMB code.*

neutrinos are non-relativistic particles at late times, they contribute to the total matter density of the universe Ω*M*, so that Ω*<sup>M</sup>* ¼ ΩC*DM* þ Ω*<sup>b</sup>* þ Ω*ν*, where ΩC*DM*, Ω*<sup>b</sup>* and Ω*<sup>ν</sup>* are the dimensionless density parameters for CDM, baryons and neutrinos, respectively. The density background is affected by massive neutrinos such that a perturbation in the density field is well described by [39, 42] as follows:

$$
\delta\_m = \left(\mathbb{1} - f\_\nu\right) \delta\_{\text{CDM}} + f\_\nu \delta\_\nu, \quad \text{where } f\_\nu \equiv \frac{\Omega\_\nu}{\Omega\_{\text{CDM}} + \Omega\_\nu}, \tag{15}
$$

*δν* being the neutrino perturbations and Ω*<sup>ν</sup>* the density contribution related to P massive neutrinos that can be expressed in terms of the total neutrino mass, *m<sup>ν</sup>* � *i mν<sup>i</sup>* , as follows:

$$
\Omega\_{\nu} = \frac{\Sigma\_{l} m\_{\nu\_{i}}}{93.14 \, h^{2} \, eV}. \tag{16}
$$

Currently, several observations provide limits on the total neutrino mass under the assumption of standard GR [43]. Depending of their mass, neutrinos can affect different quantities such as the matter-radiation equality and at the same time imprint features in cosmological observables like the clustering, the matter power spectrum, the halo mass function and the redshift-space distortions [18, 44–47].

## **3. Halo mass function and clustering analysis**

Considering the impact of MG and massive neutrinos in the clustering, a powerful cosmological test to discriminate among these scenarios is provided by the redshift-space distortions (RSD), that is, the shift in the position of the tracers due to their peculiar motions. For this purpose, cosmological simulations have become a powerful tool for testing theoretical predictions and to lead observational projects.

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies DOI: http://dx.doi.org/10.5772/intechopen.92205*

In this context, the formation and evolution of cosmic structures can be understood as a dynamical system of many particles, which trace the underlying mass distribution in a certain cosmological model. The N-body simulations, methods and algorithms have progressed continuously, achieving a high resolution to resolve finer structures with millions of particles, reducing the gap between theory and observations. For a detailed description on fundamentals of cosmological simulations, see, e.g., [48–51].

Since the formation and evolution of cosmic structures is based on the growth of small fluctuations in the density field, it is expected that the amplitude of these initial perturbations have the correct value at late times to match the observed clustering today. An analytical development based on perturbation theory makes possible to follow the growth of structures to a certain extent using the linear approximation, being valid as long as *δρ* ≪ *ρ*. Nevertheless, these calculations are limited and cannot be extrapolated to explain completely the observational data; they break down on a scale where the density contrast *δ* � 1. Moreover, beyond the linear regime, the observed structures have a density contrast in a wide range from cosmic voids with *<sup>δ</sup>* � �1 to *<sup>δ</sup>* � <sup>10</sup><sup>6</sup> and larger. It makes necessary a more elaborated description of the perturbations in the non-linear regime, which can be achieved using higher-order perturbation theory or numerical simulations [48].

In this section we show a complementary analysis to the one performed in [52] in order to investigate the clustering in the context of modified gravity with massive neutrinos. We used a subset of the DUSTGRAIN-*pathfinder* runs [36], which implement the Hu-Sawicki *f*(*R*) model including massive neutrinos and whose cosmological parameters are consistent with Planck 2015 constraints [53].

### **3.1 Halo mass function**

neutrinos are non-relativistic particles at late times, they contribute to the total matter density of the universe Ω*M*, so that Ω*<sup>M</sup>* ¼ ΩC*DM* þ Ω*<sup>b</sup>* þ Ω*ν*, where ΩC*DM*, Ω*<sup>b</sup>* and Ω*<sup>ν</sup>* are the dimensionless density parameters for CDM, baryons and neutrinos, respectively. The density background is affected by massive neutrinos such that a

*Evolution of the massive neutrino density, normalised with respect today's critical density, as a function of the scale factor* a*, for a neutrino mass* m? *= 0.1 eV. The pink solid line shows the correct solution for ρν*ð Þ *a , while the green and blue dashed lines show the limit at early and late time, respectively. The early time limit is*

*(black) marks the non-relativistic transition for neutrinos of this given mass. The figure was computed by using*

*, underlining that at early times neutrinos are relativistic and behave like radiation. The late*

*, following the same evolution of matter. The short vertical dashed line*

� �*δCDM* <sup>þ</sup> *<sup>f</sup> <sup>ν</sup>δν*, where *<sup>f</sup> <sup>ν</sup>* � <sup>Ω</sup>*<sup>ν</sup>*

*δν* being the neutrino perturbations and Ω*<sup>ν</sup>* the density contribution related to

Currently, several observations provide limits on the total neutrino mass under the assumption of standard GR [43]. Depending of their mass, neutrinos can affect different quantities such as the matter-radiation equality and at the same time imprint features in cosmological observables like the clustering, the matter power spectrum, the halo mass function and the redshift-space distortions [18, 44–47].

Considering the impact of MG and massive neutrinos in the clustering, a powerful cosmological test to discriminate among these scenarios is provided by the redshift-space distortions (RSD), that is, the shift in the position of the tracers due to their peculiar motions. For this purpose, cosmological simulations have become a powerful tool for testing theoretical predictions and to lead observational projects.

massive neutrinos that can be expressed in terms of the total neutrino mass, *m<sup>ν</sup>* �

<sup>Ω</sup>*<sup>ν</sup>* <sup>¼</sup> <sup>Σ</sup>*im<sup>ν</sup><sup>i</sup>*

Ω*CDM* þ Ω*<sup>ν</sup>*

<sup>93</sup>*:*<sup>14</sup> *<sup>h</sup>*<sup>2</sup> *eV :* (16)

, (15)

perturbation in the density field is well described by [39, 42] as follows:

*δ<sup>m</sup>* ¼ 1 � *f <sup>ν</sup>*

**3. Halo mass function and clustering analysis**

P

**200**

*i mν<sup>i</sup>*

**Figure 2.**

*proportional to* a*<sup>4</sup>*

*the CAMB code.*

*time limit, instead, is proportional to* a*<sup>3</sup>*

*Progress in Fine Particle Plasmas*

, as follows:

As CDM haloes form from collapsing regions that detach from the background density field, their abundance can be related to the volume fraction of a Gaussian density smoothed on a radius *R* above a critical collapse threshold *δ<sup>c</sup>* [54]. The comoving number density of the haloes is strictly related to underlining cosmological model, such that within a mass interval ½ � *M*, *M* þ *dM* , the halo mass function is given by

$$\frac{\mathrm{d}n(M,z)}{\mathrm{d}M} = f(\sigma(M,z)) \frac{\overline{\rho}}{M} \frac{\mathrm{d} \ln \sigma(M,z)}{\mathrm{d}M},\tag{17}$$

where *f*ð Þ *σ* is the multiplicity function, *σ* the RMS variance of the linear density field smoothed on scale *R*(*M*) and *ρ* the mean matter density. The product *f*ð Þ *σ*ð Þ *M*, *z ρ* quantifies the amount of mass contained in fluctuations of typical mass *<sup>M</sup>* <sup>¼</sup> <sup>4</sup> <sup>3</sup> *<sup>π</sup>R*<sup>3</sup> *ρ*. The simplest argument to compute analytically the multiplicity function *f*ð Þ *σ* comes from the spherical collapse theory, following the [54] formalism, such that a perturbation is supposed to collapse when it reaches the threshold *δ<sup>c</sup>* ≃1*:*68, by assuming that the probability distribution for a perturbation on a scale *M* is a Gaussian function with variance *σ*<sup>2</sup> *<sup>M</sup>*, resulting in

$$f(\sigma(\mathcal{M})) = \sqrt{\frac{2}{\pi}} \frac{\delta\_c}{\sigma(\mathcal{M})} \exp\left(-\frac{\delta\_c^2}{2\sigma^2(\mathcal{M})}\right). \tag{18}$$

Another approach to determine *f*ð Þ *σ* is given by accurate fitting functions, like the proposed by [55], which extends phenomenologically the results of [54]. For [55] the function *f*ð Þ *σ* is expected to be universal to the changes in redshift and cosmology and is parameterized as follows:

$$f(\sigma) = A \left[ \left( \frac{\sigma}{b} \right)^{-d} + \sigma^{-\epsilon} \right] \exp \left( -d/\sigma^2 \right), \tag{19}$$

**Figure 4** compares the halo mass functions of the different DUSTGRAIN-*pathfinder* simulations, computed at *z* ¼ 0 (left column), *z* ¼ 1 (central column) and *z* ¼ 1*:*6 (right column). The lower panels show the percentage difference with respect to the ΛCDM model. It is possible to see that the effect of *f*(*R*) and massive neutrinos on the dynamical evolution of the matter density field results in different halo formation epochs and different number density of collapsed systems. In particular, the *fR*4 model (blue) is the most deviated model from the standard scenario, whereas the *fR*6, *fR*6\_0*:*06*eV* and *fR*6\_0*:*1*eV* models mimic the ΛCDM behaviour

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies*

*DOI: http://dx.doi.org/10.5772/intechopen.92205*

To quantify the halo clustering, we used the two-point correlation function (2PCF) that can be defined as the joint probability of finding a pair of objects at certain spatial separation given a tracer distribution. To measure the full 2PCF

*DD r*ð Þ� , *μ* 2*DR r*ð Þþ , *μ RR r*ð Þ , *μ*

*ξ*ð Þ� *s*, *μ ξ*0ð Þ*s L*0ð Þþ *μ ξ*2ð Þ*s L*2ð Þþ *μ ξ*4ð Þ*s L*4ð Þ *μ* , (23)

where *μ* is the cosine of the angle between the line of sight and the comoving halo pair separation, *r*, and *DD*, *RR* and *DR* represent the normalised number of data-data, random-random and data-random pairs, respectively. This estimator is almost unbiased with minimum variance; this is the reason why it is preferred over the other estimators regarding the clustering measurements. Since the possible deviations from GR are more evident on small scales, we consider an intermediate non-linear range from 1 *h*�<sup>1</sup> Mpc to 50 *h*�<sup>1</sup> Mpc and random samples ten times larger than the halo ones. Then, in order to examine how significant is the RSD correction, it is convenient to expand the 2D 2PCF in the orthonormal basis of the

*RR r*ð Þ , *<sup>μ</sup>* , (22)

*dμξ*ð Þ *r*, *μ Ll*ð Þ *μ :* (24)

denoted as *ξ*ð Þ *r*, *μ* , we used the Landy-Szalay estimator given by [56]:

where each coefficient corresponds to the *l*th multipole moment:

2*l* þ 1 2

ðþ<sup>1</sup> �1

*The mass function of CDM haloes for all the models of the DUSTGRAIN-pathfinder at three different redshifts: z* ¼ 0 *(left column), z* ¼ 1 *(central column) and z* ¼ 1*:*6 *(right column). The lower panels show the percentage difference with respect to the ΛCDM model. As in Figure 3 the dashed line represents the theoretical prediction*

*ξl*ð Þ¼ *r*

^*ξ*ð Þ¼ *<sup>r</sup>*, *<sup>μ</sup>*

Legendre polynomials *Ll*ð Þ *μ* [57] such that

**Figure 4.**

*by [55].*

**203**

over a wide range of masses.

**3.2 Clustering analysis**

where *A* is an amplitude of the mass function and *a*, *b*, *c* and *d* are free parameters that depend on halo definition. The variance *σ*<sup>2</sup> is usually given by

$$\sigma^2(R(M), z) = \int \frac{P(k, z)}{2\pi^2} W^2(kR(M)) k^2 \mathrm{d}k,\tag{20}$$

where *P k*ð Þ is the linear matter power spectrum as a function of the wave number *k* and *W* is the Fourier transform of the real-space top-hat window function of radius *R*. A fundamental feature of the mass function is that it decreases monotonically with increasing masses; furthermore, its dependency on cosmology is encoded in the variance *σ*<sup>2</sup> , as shown by the integrand of Eq. (20). From the point of view of N-body simulations, an approach to compute the mass function is given straightforward from Eq. (17), by counting the number of haloes *Nh* above a certain mass threshold *Mmin* in a comoving volume *V* such as

$$N\_h = A \int\_{x\_{\rm min}}^{x\_{\rm max}} \int\_{M\_{\rm min}}^{\infty} \frac{\mathrm{d}V}{\mathrm{d}z} n(M, z) \mathrm{d}M \mathrm{d}z,\tag{21}$$

where *A* is the area, *zmin* and *zmax* are the redshift boundaries and *dV=dz* is the comoving volume element.

**Figure 3** shows the mass function of CDM haloes measured for all models of the DUSTGRAIN-*pathfinder* runs at six different redshifts *z* ¼ 0, 0*:*5, 1, 1*:*4, 1*:*6, 2. Each panel contains the mass function, per each model as labelled, to track its evolution in redshift. As reference, the black dashed line represents the theoretical expectation by [55] for a flat ΛCDM model. As expected, massive haloes are less abundant with respect to smaller ones in a fixed comoving volume. The mass function decreases with redshift, since at earlier times the density field is smoother than at late times. The plot is logarithmic, meaning that the number density of large mass haloes falls off by several orders of magnitude over the range of redshifts shown. The *f*(*R*) models both with and without neutrinos reproduce in very well agreement this pattern, but only at really high masses, significant differences appear.

#### **Figure 3.**

*The mass function of CDM haloes per each model of the DUSTGRAIN-pathfinder project at six different redshifts z* ¼ 0, 0*:*5, 1, 1*:*4, 1*:*6, 2 *as labelled. Each panel corresponds to one model as labelled, and the dashed line represents the theoretical expectation by [55], assuming a flat ΛCDM model.*

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies DOI: http://dx.doi.org/10.5772/intechopen.92205*

**Figure 4** compares the halo mass functions of the different DUSTGRAIN-*pathfinder* simulations, computed at *z* ¼ 0 (left column), *z* ¼ 1 (central column) and *z* ¼ 1*:*6 (right column). The lower panels show the percentage difference with respect to the ΛCDM model. It is possible to see that the effect of *f*(*R*) and massive neutrinos on the dynamical evolution of the matter density field results in different halo formation epochs and different number density of collapsed systems. In particular, the *fR*4 model (blue) is the most deviated model from the standard scenario, whereas the *fR*6, *fR*6\_0*:*06*eV* and *fR*6\_0*:*1*eV* models mimic the ΛCDM behaviour over a wide range of masses.

### **3.2 Clustering analysis**

*<sup>f</sup>*ð Þ¼ *<sup>σ</sup> <sup>A</sup> <sup>σ</sup>*

ð Þ¼ *R M*ð Þ, *z*

mass threshold *Mmin* in a comoving volume *V* such as

*Nh* ¼ *A*

ð*zmax zmin*

ð<sup>∞</sup> *Mmin*

this pattern, but only at really high masses, significant differences appear.

*The mass function of CDM haloes per each model of the DUSTGRAIN-pathfinder project at six different redshifts z* ¼ 0, 0*:*5, 1, 1*:*4, 1*:*6, 2 *as labelled. Each panel corresponds to one model as labelled, and the dashed*

*line represents the theoretical expectation by [55], assuming a flat ΛCDM model.*

*σ*2

encoded in the variance *σ*<sup>2</sup>

*Progress in Fine Particle Plasmas*

comoving volume element.

**Figure 3.**

**202**

*b* � ��*<sup>a</sup>* <sup>þ</sup> *<sup>σ</sup>*�*<sup>c</sup>* h i

eters that depend on halo definition. The variance *σ*<sup>2</sup> is usually given by

where *A* is an amplitude of the mass function and *a*, *b*, *c* and *d* are free param-

ð Þ *kR M*ð Þ *<sup>k</sup>*<sup>2</sup>

, as shown by the integrand of Eq. (20). From the point of

<sup>ð</sup> *P k*ð Þ , *<sup>z</sup>* <sup>2</sup>*π*<sup>2</sup> *<sup>W</sup>*<sup>2</sup>

where *P k*ð Þ is the linear matter power spectrum as a function of the wave number *k* and *W* is the Fourier transform of the real-space top-hat window function of radius *R*. A fundamental feature of the mass function is that it decreases monotonically with increasing masses; furthermore, its dependency on cosmology is

view of N-body simulations, an approach to compute the mass function is given straightforward from Eq. (17), by counting the number of haloes *Nh* above a certain

> d*V* d*z*

where *A* is the area, *zmin* and *zmax* are the redshift boundaries and *dV=dz* is the

**Figure 3** shows the mass function of CDM haloes measured for all models of the DUSTGRAIN-*pathfinder* runs at six different redshifts *z* ¼ 0, 0*:*5, 1, 1*:*4, 1*:*6, 2. Each panel contains the mass function, per each model as labelled, to track its evolution in redshift. As reference, the black dashed line represents the theoretical expectation by [55] for a flat ΛCDM model. As expected, massive haloes are less abundant with respect to smaller ones in a fixed comoving volume. The mass function decreases with redshift, since at earlier times the density field is smoother than at late times. The plot is logarithmic, meaning that the number density of large mass haloes falls off by several orders of magnitude over the range of redshifts shown. The *f*(*R*) models both with and without neutrinos reproduce in very well agreement

exp �*d=σ*<sup>2</sup> � �, (19)

*n M*ð Þ , *z* d*M*d*z*, (21)

d*k*, (20)

To quantify the halo clustering, we used the two-point correlation function (2PCF) that can be defined as the joint probability of finding a pair of objects at certain spatial separation given a tracer distribution. To measure the full 2PCF denoted as *ξ*ð Þ *r*, *μ* , we used the Landy-Szalay estimator given by [56]:

$$\hat{\xi}(r,\mu) = \frac{DD(r,\mu) - 2DR(r,\mu) + RR(r,\mu)}{RR(r,\mu)},\tag{22}$$

where *μ* is the cosine of the angle between the line of sight and the comoving halo pair separation, *r*, and *DD*, *RR* and *DR* represent the normalised number of data-data, random-random and data-random pairs, respectively. This estimator is almost unbiased with minimum variance; this is the reason why it is preferred over the other estimators regarding the clustering measurements. Since the possible deviations from GR are more evident on small scales, we consider an intermediate non-linear range from 1 *h*�<sup>1</sup> Mpc to 50 *h*�<sup>1</sup> Mpc and random samples ten times larger than the halo ones. Then, in order to examine how significant is the RSD correction, it is convenient to expand the 2D 2PCF in the orthonormal basis of the Legendre polynomials *Ll*ð Þ *μ* [57] such that

$$\xi(\mathfrak{s},\mu) \equiv \xi\_0(\mathfrak{s})L\_0(\mu) + \xi\_2(\mathfrak{s})L\_2(\mu) + \xi\_4(\mathfrak{s})L\_4(\mu),\tag{23}$$

where each coefficient corresponds to the *l*th multipole moment:

$$\xi\_l(r) = \frac{2l+1}{2} \int\_{-1}^{+1} d\mu \,\xi(r,\mu) L\_l(\mu). \tag{24}$$

#### **Figure 4.**

*The mass function of CDM haloes for all the models of the DUSTGRAIN-pathfinder at three different redshifts: z* ¼ 0 *(left column), z* ¼ 1 *(central column) and z* ¼ 1*:*6 *(right column). The lower panels show the percentage difference with respect to the ΛCDM model. As in Figure 3 the dashed line represents the theoretical prediction by [55].*

**Figure 5** shows the RSD effects on the iso-correlation curves of the 2D two-point correlation function (2PCF) in the plane (*r*⊥,*r*∥), where *r*<sup>⊥</sup> and *r*<sup>∥</sup> coordinates are, respectively, the perpendicular and parallel components along the line of sight of the observer. The 2PCF has been computed for the ΛCDM catalogue of the DUSTGRAIN-*pathfinder* simulations, in real space (*left panel*), and for the corresponding sample in redshift space (*right panel*). The contours are drawn at the iso-correlation levels *ξ s*⊥, *s*<sup>∥</sup> <sup>¼</sup> <sup>0</sup>*:*3, 0*:*5, 1*:*0, 1*:*4, 2*:*2, 3*:*6, 7*:*2, 21*:*6. In real space the correlation function is undistorted, describing circular curves in this plane. In redshift space, the effect caused by the RSD is clearly visible on small scales, where the 2PCF is stretched in the direction of *r*<sup>∥</sup> (Fingers-of-God effect), and in the infall effect on large scales, the contours are squashed along the perpendicular direction (Kaiser effect).

The symmetry of the 2PCF in real space means that the full clustering signal is encoded in the monopole moment *ξ*0ð Þ*r* , while the rest of the multipole moments are statistically equal to zero. **Figure 6** shows the monopole moment, *ξ*0ð Þ*r* , of the CDM haloes for all models considered in the DUSTGRAIN-*pathfinder* project, at three different redshifts *z* ¼ 0, 1, 2. Subpanels show the percentage difference between MG models [*f*(*R*) with and without massive neutrinos] and the ΛCDM model. The monopole moments of the 2PCF of the *fR*4 and *fR*4\_0*:*3*eV* models are the ones that deviate most from ΛCDM at low redshift. This behaviour is also present in the models *fR*5, *fR*5\_0*:*15*eV* and *fR*5\_0*:*1*eV*, but it is less significant. In general, it is observed that massive neutrinos increase the clustering signal for all models, especially at high redshift, the *fR*6, *fR*6\_0*:*06*eV* and *fR*6\_0*:*1*eV* models being degenerated with respect to ΛCDM. The quadrupole and hexadecapole moments are consistent with zero at 1*σ* error bar.

Another important feature to take into account in the clustering analysis is related to the bias that is introduced when the ΛCDM model is wrongly assumed to predict the DM clustering of a *f*(*R*) universe with massive neutrinos [58]. This

> effective halo bias, h i *b* , allows to characterise the relation between the halo clustering and the underlying mass distribution. By using the theoretical effective bias proposed by [59] and its corresponding mass function, it is possible to disentangle the degeneracy with respect to *σ*8. The level of agreement with the measurements obtained from the *f*(*R*) and *f R*ð Þþ *m<sup>ν</sup>* scenarios give us a better understating of the discrepancy with the ΛCDM ones. To compute the theoretical mass function, we consider the linear CDM + baryon power spectrum as have been stated in the CDM prescription [60] and replacing *ρ<sup>m</sup>* with *ρ*CDM [42]. These quantities can be obtained with CAMB or another Boltzmann code; the linear power-spectrum for CDM,

> *The real-space 2PCF r*<sup>2</sup>*ξ*<sup>0</sup> *of CDM haloes for all the models of the DUSTGRAIN-pathfinder project at three different redshifts: z* ¼ 0 *(left column), z* ¼ 1 *(central column) and z* ¼ 2 *(right column). From top to bottom, the panels show the fR*4*, fR*5 *and fR*6 *models, respectively, compared to the results of the ΛCDM model. The error bars, shown only for the ΛCDM model for clarity reasons, are the diagonal values of the bootstrap covariance matrices used for the statistical analysis. Percentage differences between* f*(*R*), f R*ð Þþ *m<sup>ν</sup> and ΛCDM predictions are in the subpanels, while the shaded regions represent the deviation at 1σ confidence level.*

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies*

*DOI: http://dx.doi.org/10.5772/intechopen.92205*

*lin*ð Þ*<sup>k</sup>* , is expressed in terms of *<sup>P</sup>*CDMþ*<sup>b</sup>*

*T*CDMð Þ*k* and *Tb*ð Þ*k* being the transfer functions. It implies, as shown by [61], that the effect of neutrinos on the cluster abundance is well captured by rescaling the

<sup>ð</sup> *<sup>P</sup>*CDMð Þ *<sup>k</sup>*, *<sup>z</sup>*

with the CDM power spectrum obtained by rescaling the total matter power spectrum with the corresponding transfer functions, *TCDM* and *Tb*, weighted by the

<sup>2</sup>*π*<sup>2</sup> *<sup>W</sup>*<sup>2</sup>

ΩCDM*T*CDMð Þþ *k*, *z* Ωb*T*bð Þ *k*, *z T*mð Þ *k*, *z* ð Þ ΩCDM þ Ω<sup>b</sup> � �<sup>2</sup>

ð Þ *kR <sup>k</sup>*<sup>2</sup>

*lin* ð Þ*<sup>k</sup>* and *<sup>P</sup><sup>m</sup>*

*lin*ð Þ*k* , with

*:* (26)

d*k* (25)

*P*CDM

**205**

**Figure 6.**

*lin* ð Þ¼ *<sup>k</sup> <sup>T</sup>*<sup>2</sup>

CDM*=T*<sup>2</sup>

smoothed density field such that

density of each species so that

*mP<sup>m</sup>*

*<sup>σ</sup>*<sup>2</sup> ! *<sup>σ</sup>*<sup>2</sup>

*P*CDMð Þ¼ *k*, *z P*mð Þ *k*, *z*

CDMð Þ¼ *z*

#### **Figure 5.**

*Contours of the 2D two-point correlation function (2PCF) in the plane (r*⊥,*r*∥*), the r*<sup>⊥</sup> *and r*<sup>∥</sup> *coordinates are, respectively, the perpendicular and parallel components along the line of sight of the observer. The 2PCF has been computed from an N-body simulation, in real space (left panel), and for the corresponding sample in redshift space (right panel). The contours are draw at the iso-correlation levels ξ s*⊥, *s*<sup>∥</sup> <sup>¼</sup> <sup>0</sup>*:*3, 0*:*5, 1*:*0, 1*:*4, 2*:*2, 3*:*6, 7*:*2, 21*:*<sup>6</sup> *as indicated by the colour bar. Both the effect of redshift space distortions on small scales (Finger-of-God effect) and the infall effect at large scales (Kaiser effect) are clearly visible in the left panel. On small scales the 2PCF is stretched in the direction of r*∥*, and on large scales the contours are squashed along the perpendicular direction.*

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies DOI: http://dx.doi.org/10.5772/intechopen.92205*

#### **Figure 6.**

**Figure 5** shows the RSD effects on the iso-correlation curves of the 2D two-point correlation function (2PCF) in the plane (*r*⊥,*r*∥), where *r*<sup>⊥</sup> and *r*<sup>∥</sup> coordinates are, respectively, the perpendicular and parallel components along the line of sight of the observer. The 2PCF has been computed for the ΛCDM catalogue of the DUSTGRAIN-*pathfinder* simulations, in real space (*left panel*), and for the

corresponding sample in redshift space (*right panel*). The contours are drawn at the

The symmetry of the 2PCF in real space means that the full clustering signal is encoded in the monopole moment *ξ*0ð Þ*r* , while the rest of the multipole moments are statistically equal to zero. **Figure 6** shows the monopole moment, *ξ*0ð Þ*r* , of the CDM haloes for all models considered in the DUSTGRAIN-*pathfinder* project, at three different redshifts *z* ¼ 0, 1, 2. Subpanels show the percentage difference between MG models [*f*(*R*) with and without massive neutrinos] and the ΛCDM model. The monopole moments of the 2PCF of the *fR*4 and *fR*4\_0*:*3*eV* models are the ones that deviate most from ΛCDM at low redshift. This behaviour is also present in the models *fR*5, *fR*5\_0*:*15*eV* and *fR*5\_0*:*1*eV*, but it is less significant. In general, it is observed that massive neutrinos increase the clustering signal for all models, especially at high redshift, the *fR*6, *fR*6\_0*:*06*eV* and *fR*6\_0*:*1*eV* models being degenerated with respect to ΛCDM. The quadrupole and hexadecapole

Another important feature to take into account in the clustering analysis is related to the bias that is introduced when the ΛCDM model is wrongly assumed to predict the DM clustering of a *f*(*R*) universe with massive neutrinos [58]. This

*Contours of the 2D two-point correlation function (2PCF) in the plane (r*⊥,*r*∥*), the r*<sup>⊥</sup> *and r*<sup>∥</sup> *coordinates are, respectively, the perpendicular and parallel components along the line of sight of the observer. The 2PCF has been computed from an N-body simulation, in real space (left panel), and for the corresponding sample in redshift space*

*as indicated by the colour bar. Both the effect of redshift space distortions on small scales (Finger-of-God effect) and the infall effect at large scales (Kaiser effect) are clearly visible in the left panel. On small scales the 2PCF is stretched in the direction of r*∥*, and on large scales the contours are squashed along the perpendicular*

<sup>¼</sup> <sup>0</sup>*:*3, 0*:*5, 1*:*0, 1*:*4, 2*:*2, 3*:*6, 7*:*2, 21*:*<sup>6</sup>

correlation function is undistorted, describing circular curves in this plane. In redshift space, the effect caused by the RSD is clearly visible on small scales, where the 2PCF is stretched in the direction of *r*<sup>∥</sup> (Fingers-of-God effect), and in the infall effect on large scales, the contours are squashed along the perpendicular direction

<sup>¼</sup> <sup>0</sup>*:*3, 0*:*5, 1*:*0, 1*:*4, 2*:*2, 3*:*6, 7*:*2, 21*:*6. In real space the

iso-correlation levels *ξ s*⊥, *s*<sup>∥</sup>

*Progress in Fine Particle Plasmas*

moments are consistent with zero at 1*σ* error bar.

*(right panel). The contours are draw at the iso-correlation levels ξ s*⊥, *s*<sup>∥</sup>

(Kaiser effect).

**Figure 5.**

*direction.*

**204**

*The real-space 2PCF r*<sup>2</sup>*ξ*<sup>0</sup> *of CDM haloes for all the models of the DUSTGRAIN-pathfinder project at three different redshifts: z* ¼ 0 *(left column), z* ¼ 1 *(central column) and z* ¼ 2 *(right column). From top to bottom, the panels show the fR*4*, fR*5 *and fR*6 *models, respectively, compared to the results of the ΛCDM model. The error bars, shown only for the ΛCDM model for clarity reasons, are the diagonal values of the bootstrap covariance matrices used for the statistical analysis. Percentage differences between* f*(*R*), f R*ð Þþ *m<sup>ν</sup> and ΛCDM predictions are in the subpanels, while the shaded regions represent the deviation at 1σ confidence level.*

effective halo bias, h i *b* , allows to characterise the relation between the halo clustering and the underlying mass distribution. By using the theoretical effective bias proposed by [59] and its corresponding mass function, it is possible to disentangle the degeneracy with respect to *σ*8. The level of agreement with the measurements obtained from the *f*(*R*) and *f R*ð Þþ *m<sup>ν</sup>* scenarios give us a better understating of the discrepancy with the ΛCDM ones. To compute the theoretical mass function, we consider the linear CDM + baryon power spectrum as have been stated in the CDM prescription [60] and replacing *ρ<sup>m</sup>* with *ρ*CDM [42]. These quantities can be obtained with CAMB or another Boltzmann code; the linear power-spectrum for CDM, *P*CDM *lin* ð Þ¼ *<sup>k</sup> <sup>T</sup>*<sup>2</sup> CDM*=T*<sup>2</sup> *mP<sup>m</sup> lin*ð Þ*<sup>k</sup>* , is expressed in terms of *<sup>P</sup>*CDMþ*<sup>b</sup> lin* ð Þ*<sup>k</sup>* and *<sup>P</sup><sup>m</sup> lin*ð Þ*k* , with *T*CDMð Þ*k* and *Tb*ð Þ*k* being the transfer functions. It implies, as shown by [61], that the effect of neutrinos on the cluster abundance is well captured by rescaling the smoothed density field such that

$$
\sigma^2 \to \sigma\_{\rm CDM}^2(z) = \int \frac{P\_{\rm CDM}(k, z)}{2\pi^2} W^2(kR) k^2 \,\mathrm{d}k \tag{25}
$$

with the CDM power spectrum obtained by rescaling the total matter power spectrum with the corresponding transfer functions, *TCDM* and *Tb*, weighted by the density of each species so that

$$P\_{\rm CDM}(k,z) = P\_{\rm m}(k,z) \left( \frac{\Omega\_{\rm CDM} T\_{\rm CDM}(k,z) + \Omega\_{\rm b} T\_{\rm b}(k,z)}{T\_{\rm m}(k,z)(\Omega\_{\rm CDM} + \Omega\_{\rm b})} \right)^2. \tag{26}$$

#### **Figure 7.**

*The coloured solid lines represent the apparent effective halo bias, b*h i*, as a function of redshift, averaged in the range 10 h*�<sup>1</sup> *Mpc* <*r* <*50 h*�<sup>1</sup> *Mpc. Black lines show the theoretical ΛCDM effective bias predicted by [59] (dashed), normalised to the σ*<sup>8</sup> *values of each DUSTGRAIN-pathfinder simulation, while the cyan-shaded areas show a 10% error. The upper set of panels shows the results, considering the total power spectrum, while the lower set of panels shows the results when the CDM+baryon power spectrum is used instead.*

rate *f*ð Þ Ω*<sup>m</sup>* and the linear bias. All numerical tasks were performed with the

*panels), at three different redshifts z* ¼ 0*:*5, 1*:*0 *and 1.6 from left to right.*

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies*

*DOI: http://dx.doi.org/10.5772/intechopen.92205*

The Kaiser formula is a good description of the RSD only at very large scales, where non-linear effects can be neglected, but it does not describe accurately the non-linear regime. Thus, with the aim of extracting information from the RSD signal at non-linear regime and considering the increasing precision of recent and upcoming surveys, many more approaches have been proposed. There is a vast literature that shows the efforts to model the RSD beyond the linear Kaiser model [63–66], some of them making use of a phenomenological description of the velocity field and others, instead, taking into account higher orders in perturbation theory since, in principle, there is no reason to stop at linear order. Other

*Covariance matrices from the analysis of the redshift-space monopole and quadrupole moments of CDM haloes with Bootstrap errors. The models correspond to ΛCDM (upper panels),* fR*4 (central panels) and* fR*5 (bottom*

approaches do a combination of both frameworks. A simple alternative to model the redshift-space 2PCF at small scales consists of extending the Kaiser formula, by adding a phenomenological damping factor that plays the role of a pairwise velocity distribution. It can account for both linear and non-linear dynamics. Therefore, to construct the likelihood, we consider this model sometimes called *dispersion model* [67], which introduces a damping function to describe the distortions in the clustering at small scales (Fingers-of-God). This model is enough accurate to quantify the relative differences between *f*(*R*) models with massive neutrinos and ΛCDM. For the Bayesian analysis, the dispersion model is fully described by three parameters, *f σ*8, *bσ*<sup>8</sup> and Σ*S*, that we constrain by minimising numerically the

**Figure 8** shows the normalised covariance matrices *Ci*,*<sup>j</sup>=* ffiffiffiffiffiffiffiffiffiffiffiffiffiffi

redshift-space monopole and quadrupole moments of CDM haloes with bootstrap errors resampling at three different redshifts *z* ¼ 0*:*5, 1*:*0 and 1.6. As it can be appreciated, the covariance matrices represent how the scatter propagates into the

<sup>1</sup> Freely available at the public GitHub repository https://github.com/federicomarulli/CosmoBolognaLib.

*Ci*,*iC <sup>j</sup>*,*<sup>j</sup>* � � p of the

CosmoBolognaLib<sup>1</sup> [62].

**Figure 8.**

negative log-likelihood.

**207**

The impact on the bias when the CDM prescription is not considered can be appreciated in **Figure 7**. In most cases this correction is small with the exception of the *fR*4\_0*:*3 eV model. A detailed discussion on these results can be found in [52].

#### **3.3 Modelling the redshift-space distortions**

In a realistic case, spectroscopic surveys observe a combination of density and velocity fields in redshift space. The observed redshift is a combination of cosmological effects plus an additional term caused by the peculiar motions along the line of sight of the observer. This combination makes the redshift-space catalogues appear distorted with respect to the real-space ones, and they can be reproduced from N-body simulations since the positions and velocities are known. Currently, the modelling of the redshift-space distortions provide a powerful tool to test the gravity theory by exploring the spatial statistics encoded in the 2PCF, which is anisotropic due to the dynamic distortions.

We consider the first two even multipoles of the 2PCF, *ξ*<sup>0</sup> and *ξ*2, taking into account that odd multipole moments vanish by symmetry. The analysis consisted of modelling (i) the individual multipole moments and (ii) both multipole moments simultaneously. To derive cosmological constraints from the clustering signal and quantify the effects of *f*(*R*) gravity and massive neutrinos on RSD, we performed a Bayesian analysis to set constraints on the linear growth *Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies DOI: http://dx.doi.org/10.5772/intechopen.92205*

#### **Figure 8.**

The impact on the bias when the CDM prescription is not considered can be appreciated in **Figure 7**. In most cases this correction is small with the exception of the *fR*4\_0*:*3 eV model. A detailed discussion on these results can be found

*the lower set of panels shows the results when the CDM+baryon power spectrum is used instead.*

*The coloured solid lines represent the apparent effective halo bias, b*h i*, as a function of redshift, averaged in the range 10 h*�<sup>1</sup> *Mpc* <*r* <*50 h*�<sup>1</sup> *Mpc. Black lines show the theoretical ΛCDM effective bias predicted by [59] (dashed), normalised to the σ*<sup>8</sup> *values of each DUSTGRAIN-pathfinder simulation, while the cyan-shaded areas show a 10% error. The upper set of panels shows the results, considering the total power spectrum, while*

In a realistic case, spectroscopic surveys observe a combination of density and velocity fields in redshift space. The observed redshift is a combination of cosmological effects plus an additional term caused by the peculiar motions along the line of sight of the observer. This combination makes the redshift-space catalogues appear distorted with respect to the real-space ones, and they can be reproduced from N-body simulations since the positions and velocities are known. Currently, the modelling of the redshift-space distortions provide a powerful tool to test the gravity theory by exploring the spatial statistics encoded in the 2PCF, which is

We consider the first two even multipoles of the 2PCF, *ξ*<sup>0</sup> and *ξ*2, taking into account that odd multipole moments vanish by symmetry. The analysis consisted of modelling (i) the individual multipole moments and (ii) both multipole moments simultaneously. To derive cosmological constraints from the clustering signal and quantify the effects of *f*(*R*) gravity and massive neutrinos on RSD, we performed a Bayesian analysis to set constraints on the linear growth

in [52].

**206**

**Figure 7.**

*Progress in Fine Particle Plasmas*

**3.3 Modelling the redshift-space distortions**

anisotropic due to the dynamic distortions.

*Covariance matrices from the analysis of the redshift-space monopole and quadrupole moments of CDM haloes with Bootstrap errors. The models correspond to ΛCDM (upper panels),* fR*4 (central panels) and* fR*5 (bottom panels), at three different redshifts z* ¼ 0*:*5, 1*:*0 *and 1.6 from left to right.*

rate *f*ð Þ Ω*<sup>m</sup>* and the linear bias. All numerical tasks were performed with the CosmoBolognaLib<sup>1</sup> [62].

The Kaiser formula is a good description of the RSD only at very large scales, where non-linear effects can be neglected, but it does not describe accurately the non-linear regime. Thus, with the aim of extracting information from the RSD signal at non-linear regime and considering the increasing precision of recent and upcoming surveys, many more approaches have been proposed. There is a vast literature that shows the efforts to model the RSD beyond the linear Kaiser model [63–66], some of them making use of a phenomenological description of the velocity field and others, instead, taking into account higher orders in perturbation theory since, in principle, there is no reason to stop at linear order. Other approaches do a combination of both frameworks. A simple alternative to model the redshift-space 2PCF at small scales consists of extending the Kaiser formula, by adding a phenomenological damping factor that plays the role of a pairwise velocity distribution. It can account for both linear and non-linear dynamics. Therefore, to construct the likelihood, we consider this model sometimes called *dispersion model* [67], which introduces a damping function to describe the distortions in the clustering at small scales (Fingers-of-God). This model is enough accurate to quantify the relative differences between *f*(*R*) models with massive neutrinos and ΛCDM. For the Bayesian analysis, the dispersion model is fully described by three parameters, *f σ*8, *bσ*<sup>8</sup> and Σ*S*, that we constrain by minimising numerically the negative log-likelihood.

**Figure 8** shows the normalised covariance matrices *Ci*,*<sup>j</sup>=* ffiffiffiffiffiffiffiffiffiffiffiffiffiffi *Ci*,*iC <sup>j</sup>*,*<sup>j</sup>* � � p of the redshift-space monopole and quadrupole moments of CDM haloes with bootstrap errors resampling at three different redshifts *z* ¼ 0*:*5, 1*:*0 and 1.6. As it can be appreciated, the covariance matrices represent how the scatter propagates into the

<sup>1</sup> Freely available at the public GitHub repository https://github.com/federicomarulli/CosmoBolognaLib.

likelihood and on the final posterior probabilities of the parameters. Then, to assess the posterior distributions of the three model parameters, we perform a MCMC analysis. The fitting analysis is limited to the scale range 10≤ *r* Mpc h�<sup>1</sup> ≤ 50, assuming flat priors in the ranges 0 ≤*f σ*<sup>8</sup> ≤2, 0≤*bσ*<sup>8</sup> ≤ 3 and 0 ≤Σ*<sup>S</sup>* ≤2. **Figure 9** shows the *f σ*8-*bσ*<sup>8</sup> posterior constraints obtained from the MCMC analysis of *ξ*0, *ξ*<sup>2</sup> and *ξ*<sup>0</sup> þ *ξ*<sup>2</sup> for each mock catalogue of the DUSTGRAIN-*pathfinder.* The figure shows the constraints for all models considered in this work at *z* ¼ 1*:*0 using the monopole (orange), quadrupole (green) and monopole plus quadrupole (blue), while the intersection regions correspond to the joint analysis of the multipoles. Thus, the joint analysis of the redshift-space monopole and quadrupole is able to break the degeneracy between *f σ*<sup>8</sup> and *bσ*<sup>8</sup> even in the presence of massive neutrinos. In **Figure 10**, we show the theoretical behaviour of the linear distortion parameter and the growth factor as a function of the redshift for each family of models, assuming a flat ΛCDM model.

From the *f σ*8-*bσ*<sup>8</sup> posterior contours, at 1 � 2*σ* confidence levels, the results suggest that the clustering information encoded in the two first non-null multipole moments of the 2PCF can discriminate the alternative MG models considered in this work at *z*≳1. At low redshifts the *f*(*R*) models studied are statistically indistinguishable from ΛCDM, and further studies are required to break this degeneracy.

**4. Conclusions**

*a flat ΛCDM model.*

**Figure 10.**

models.

**209**

In this chapter we have introduced the theoretical framework of modern cosmology in present massive neutrinos. We emphasize on the structure formation and on the statistical description of the density field as well as the measurements of galaxy clustering and discuss the redshift-space distortions and the differences between clustering in real and redshift space, considering that in the recent years, the spatial distribution of matter on cosmological scales has become one of the most efficient probes to investigate the properties of the universe, such as test gravity theories on large scales, to explore the dark sector and the origin of the accelerated expansion of the universe as well as a probe to constrain alternative cosmological

*Theoretical expectation of the linear distortion parameter β (upper panel) and the growth rate f σ*<sup>8</sup> *(lower panel) as a function of the redshift for each family of mocks from the DUSTGRAIN-pathfinder runs, assuming*

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies*

*DOI: http://dx.doi.org/10.5772/intechopen.92205*

In the context of models based on modified gravity and massive neutrino cosmologies, we investigated the spatial properties of the large-scale structure by exploiting the DUSTGRAIN-*pathfinder* simulations that follow, simultaneously, the effects of *f*(*R*) gravity and massive neutrinos. These are two of the most interesting scenarios that have been recently explored to account for possible observational deviations from the standard ΛCDM model. In particular, we studied whether redshift-space distortions in the 2PCF multipole moments can be effective, breaking the cosmic degeneracy between these two effects. We analysed the redshiftspace distortions in the clustering of dark matter haloes at different redshifts, focusing on the monopole and quadrupole moments of the two-point correlation function, both in real and redshift space. The deviations with respect to ΛCDM model have been quantified in terms of the linear growth rate parameter. We found that multipole moments of the 2PCF from redshift-space distortions provide a useful probe to discriminate between ΛCDM and modified gravity models, especially at high redshifts (*z*≳1), even in the presence of massive neutrinos. The linear

#### **Figure 9.**

*Contours at* 1 � 2*σ confidence level of the f σ*<sup>8</sup> � *bσ*<sup>8</sup> *posterior distributions, obtained from the MCMC analysis in redshift space for 2PCF multipoles of CDM haloes. The contours correspond to z* ¼ 1*:*0 *with the monopole in orange and the quadrupole in green.*

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies DOI: http://dx.doi.org/10.5772/intechopen.92205*

**Figure 10.**

likelihood and on the final posterior probabilities of the parameters. Then, to assess the posterior distributions of the three model parameters, we perform a MCMC analysis. The fitting analysis is limited to the scale range 10≤ *r* Mpc h�<sup>1</sup> ≤ 50, assuming flat priors in the ranges 0 ≤*f σ*<sup>8</sup> ≤2, 0≤*bσ*<sup>8</sup> ≤ 3 and 0 ≤Σ*<sup>S</sup>* ≤2. **Figure 9** shows the *f σ*8-*bσ*<sup>8</sup> posterior constraints obtained from the MCMC analysis of *ξ*0, *ξ*<sup>2</sup> and *ξ*<sup>0</sup> þ *ξ*<sup>2</sup> for each mock catalogue of the DUSTGRAIN-*pathfinder.* The figure shows the constraints for all models considered in this work at *z* ¼ 1*:*0 using the monopole (orange), quadrupole (green) and monopole plus quadrupole (blue), while the intersection regions correspond to the joint analysis of the multipoles. Thus, the joint analysis of the redshift-space monopole and quadrupole is able to break the degeneracy between *f σ*<sup>8</sup> and *bσ*<sup>8</sup> even in the presence of massive neutrinos. In **Figure 10**, we show the theoretical behaviour of the linear distortion parameter and the growth factor as a function of the redshift for each family of

From the *f σ*8-*bσ*<sup>8</sup> posterior contours, at 1 � 2*σ* confidence levels, the results suggest that the clustering information encoded in the two first non-null multipole moments of the 2PCF can discriminate the alternative MG models considered in this work at *z*≳1. At low redshifts the *f*(*R*) models studied are statistically indistinguishable from ΛCDM, and further studies are required to break this degeneracy.

*Contours at* 1 � 2*σ confidence level of the f σ*<sup>8</sup> � *bσ*<sup>8</sup> *posterior distributions, obtained from the MCMC analysis in redshift space for 2PCF multipoles of CDM haloes. The contours correspond to z* ¼ 1*:*0 *with the monopole in*

models, assuming a flat ΛCDM model.

*Progress in Fine Particle Plasmas*

**Figure 9.**

**208**

*orange and the quadrupole in green.*

*Theoretical expectation of the linear distortion parameter β (upper panel) and the growth rate f σ*<sup>8</sup> *(lower panel) as a function of the redshift for each family of mocks from the DUSTGRAIN-pathfinder runs, assuming a flat ΛCDM model.*

## **4. Conclusions**

In this chapter we have introduced the theoretical framework of modern cosmology in present massive neutrinos. We emphasize on the structure formation and on the statistical description of the density field as well as the measurements of galaxy clustering and discuss the redshift-space distortions and the differences between clustering in real and redshift space, considering that in the recent years, the spatial distribution of matter on cosmological scales has become one of the most efficient probes to investigate the properties of the universe, such as test gravity theories on large scales, to explore the dark sector and the origin of the accelerated expansion of the universe as well as a probe to constrain alternative cosmological models.

In the context of models based on modified gravity and massive neutrino cosmologies, we investigated the spatial properties of the large-scale structure by exploiting the DUSTGRAIN-*pathfinder* simulations that follow, simultaneously, the effects of *f*(*R*) gravity and massive neutrinos. These are two of the most interesting scenarios that have been recently explored to account for possible observational deviations from the standard ΛCDM model. In particular, we studied whether redshift-space distortions in the 2PCF multipole moments can be effective, breaking the cosmic degeneracy between these two effects. We analysed the redshiftspace distortions in the clustering of dark matter haloes at different redshifts, focusing on the monopole and quadrupole moments of the two-point correlation function, both in real and redshift space. The deviations with respect to ΛCDM model have been quantified in terms of the linear growth rate parameter. We found that multipole moments of the 2PCF from redshift-space distortions provide a useful probe to discriminate between ΛCDM and modified gravity models, especially at high redshifts (*z*≳1), even in the presence of massive neutrinos. The linear

growth rate constraints that we obtain from all the analysed *f R*ð Þþ *m<sup>ν</sup>* mock catalogues are statistically distinguishable from ΛCDM predictions at high redshifts.

**References**

[1] Tonry JL et al. Cosmological results

*DOI: http://dx.doi.org/10.5772/intechopen.92205*

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies*

[11] Sotiriou TP, Faraoni V. *f*(*R*) theories of gravity. Reviews of Modern Physics.

[12] De Felice A, Tsujikawa S. *f*(*R*) theories. Living Reviews in Relativity.

[13] Joyce A, Lombriser L, Schmidt F. Dark energy versus modified gravity. Annual Review of Nuclear and Particle

[14] Uzan J-P. Testing general relativity: from local to cosmological scales. Philosophical Transactions of the Royal Society of London A: Mathematical Physical and Engineering Sciences.

Science. 2016;**66**(1):95-122

2011;**369**(1957):5042-5057

[15] Will CM. The confrontation between general relativity and experiment. Living Reviews in Relativity. 2014;**17**(1):4

Extragalactic Redshift Survey

[17] Collett TE et al. A precise

[16] Pezzotta A et al. The VIMOS Public

(VIPERS). The growth of structure at 0.5 < z < 1.2 from redshift-space distortions in the clustering of the PDR-2 final sample. Astronomy & Astrophysics. 2017;**604**:A33

extragalactic test of general relativity. Science. 2018;**360**(6395):1342-1346

[18] Lesgourgues J, Pastor S. Massive neutrinos and cosmology. Physics Reports. 2006;**429**(6):307-379

[19] Motohashi H, Starobinsky AA, Yokoyama J. Cosmology based on *f*(*r*) gravity admits 1 eV sterile neutrinos. Physical Review Letters. 2013;**110**:

[20] He J-H. Weighing neutrinos in *f*(*r*) gravity. Physical Review D. 2013;**88**:

121302

103523

2010;**82**(1):451-497

2010;**13**(1):3

Astrophysical Journal. 2003;**594**(1):1-24

[2] Hamana T, Sakurai J, Koike M, Miller L. Cosmological constraints from subaru weak lensing cluster counts. Publications of the Astronomical Society

[3] Abbott TMC et al. Dark energy survey year 1 results: Cosmological constraints from galaxy clustering and weak lensing. Physical Review D. 2018;

[4] Planck Collaboration; Aghanim N

et al. Planck 2018 results. VI. Cosmological parameters. arXiv eprints, arXiv:1807.06209. July 2018

[5] Riess AG et al. Observational evidence from supernovae for an accelerating universe and a cosmological constant. The Astronomical Journal.

[6] Schmidt BP et al. The high-Z supernova search: Measuring cosmic deceleration and global curvature of the universe using type IA supernovae. Astrophysical Journal. 1998;**507**:46-63

*ω* and *λ* from 42 high-redshift

[7] Perlmutter S et al. Measurements of

supernovae. The Astrophysical Journal.

[8] Planck Collaboration, Ade PAR, et al. Planck 2015 results—xxiv. Cosmology from Sunyaev-Zeldovich cluster counts. Astronomy & Astrophysics. 2016;**594**:A24

[9] Riess AG et al. A 2.4% determination of the local value of the hubble constant. The Astrophysical Journal. 2016;**826**(1):56

[10] Luis Bernal J, Verde L, Riess AG. The trouble with h0. Journal of Cosmology and Astroparticle Physics.

1998;**116**(3):1009

1999;**517**(2):565

2016;**2016**(10):019

**211**

from high-zSupernovae. The

of Japan. 2015;**67**(3):34

**98**:043526

## **Acknowledgements**

We thank Lauro Moscardini, Federico Marulli and Alfonso Veropalumbo for the continuous development of the CosmoBolognaLib and their suggestions during this project. We also thank Carlo Giocoli and Marco Baldi for the crucial work performing the DUSTGRAIN-*pathfinder* runs.

## **Conflict of interest**

The authors declare no conflict of interest.

## **Author details**

Jorge Enrique García-Farieta\* and Rigoberto Ángel Casas Miranda Universidad Nacional de Colombia, Bogotá, Colombia

\*Address all correspondence to: joegarciafa@unal.edu.co

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies DOI: http://dx.doi.org/10.5772/intechopen.92205*

## **References**

growth rate constraints that we obtain from all the analysed *f R*ð Þþ *m<sup>ν</sup>* mock catalogues are statistically distinguishable from ΛCDM predictions at high redshifts.

We thank Lauro Moscardini, Federico Marulli and Alfonso Veropalumbo for the continuous development of the CosmoBolognaLib and their suggestions during this

project. We also thank Carlo Giocoli and Marco Baldi for the crucial work

Jorge Enrique García-Farieta\* and Rigoberto Ángel Casas Miranda

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

Universidad Nacional de Colombia, Bogotá, Colombia

provided the original work is properly cited.

\*Address all correspondence to: joegarciafa@unal.edu.co

**Acknowledgements**

*Progress in Fine Particle Plasmas*

**Conflict of interest**

**Author details**

**210**

performing the DUSTGRAIN-*pathfinder* runs.

The authors declare no conflict of interest.

[1] Tonry JL et al. Cosmological results from high-zSupernovae. The Astrophysical Journal. 2003;**594**(1):1-24

[2] Hamana T, Sakurai J, Koike M, Miller L. Cosmological constraints from subaru weak lensing cluster counts. Publications of the Astronomical Society of Japan. 2015;**67**(3):34

[3] Abbott TMC et al. Dark energy survey year 1 results: Cosmological constraints from galaxy clustering and weak lensing. Physical Review D. 2018; **98**:043526

[4] Planck Collaboration; Aghanim N et al. Planck 2018 results. VI. Cosmological parameters. arXiv eprints, arXiv:1807.06209. July 2018

[5] Riess AG et al. Observational evidence from supernovae for an accelerating universe and a cosmological constant. The Astronomical Journal. 1998;**116**(3):1009

[6] Schmidt BP et al. The high-Z supernova search: Measuring cosmic deceleration and global curvature of the universe using type IA supernovae. Astrophysical Journal. 1998;**507**:46-63

[7] Perlmutter S et al. Measurements of *ω* and *λ* from 42 high-redshift supernovae. The Astrophysical Journal. 1999;**517**(2):565

[8] Planck Collaboration, Ade PAR, et al. Planck 2015 results—xxiv. Cosmology from Sunyaev-Zeldovich cluster counts. Astronomy & Astrophysics. 2016;**594**:A24

[9] Riess AG et al. A 2.4% determination of the local value of the hubble constant. The Astrophysical Journal. 2016;**826**(1):56

[10] Luis Bernal J, Verde L, Riess AG. The trouble with h0. Journal of Cosmology and Astroparticle Physics. 2016;**2016**(10):019

[11] Sotiriou TP, Faraoni V. *f*(*R*) theories of gravity. Reviews of Modern Physics. 2010;**82**(1):451-497

[12] De Felice A, Tsujikawa S. *f*(*R*) theories. Living Reviews in Relativity. 2010;**13**(1):3

[13] Joyce A, Lombriser L, Schmidt F. Dark energy versus modified gravity. Annual Review of Nuclear and Particle Science. 2016;**66**(1):95-122

[14] Uzan J-P. Testing general relativity: from local to cosmological scales. Philosophical Transactions of the Royal Society of London A: Mathematical Physical and Engineering Sciences. 2011;**369**(1957):5042-5057

[15] Will CM. The confrontation between general relativity and experiment. Living Reviews in Relativity. 2014;**17**(1):4

[16] Pezzotta A et al. The VIMOS Public Extragalactic Redshift Survey (VIPERS). The growth of structure at 0.5 < z < 1.2 from redshift-space distortions in the clustering of the PDR-2 final sample. Astronomy & Astrophysics. 2017;**604**:A33

[17] Collett TE et al. A precise extragalactic test of general relativity. Science. 2018;**360**(6395):1342-1346

[18] Lesgourgues J, Pastor S. Massive neutrinos and cosmology. Physics Reports. 2006;**429**(6):307-379

[19] Motohashi H, Starobinsky AA, Yokoyama J. Cosmology based on *f*(*r*) gravity admits 1 eV sterile neutrinos. Physical Review Letters. 2013;**110**: 121302

[20] He J-H. Weighing neutrinos in *f*(*r*) gravity. Physical Review D. 2013;**88**: 103523

[21] Baldi M et al. Cosmic degeneracies— I. Joint n-body simulations of modified gravity and massive neutrinos. Monthly Notices of the Royal Astronomical Society. 2014;**440**(1):75-88

[22] Peel A, Pettorino V, Giocoli C, Starck J-L, Baldi M. Breaking degeneracies in modified gravity with higher (than 2nd) order weak-lensing statistics. Astronomy & Astrophysics. 2018;**619**:A38

[23] Hagstotz S, Costanzi M, Baldi M, Weller J. Joint halo mass function for modified gravity and massive neutrinos. I: Simulations and cosmological forecasts. arXiv e-prints, arXiv: 1806.07400. June 2018

[24] Lesgourgues J, Pastor S. Neutrino mass from cosmology. Advances in High Energy Physics. 2012;**2012**

[25] Peebles PJE. Principles of Physical Cosmology. Princeton, New Jersey: Princeton University Press; 1993

[26] Dodelson S. Modern Cosmology. San Diego, California: Academic Press; 2003

[27] Mukhanov V. Physical Foundations of Cosmology. Cambridge, UK: Cambridge University Press; 2005

[28] Sotiriou TP. *f*(*R*) gravity and scalar tensor theory. Classical and Quantum Gravity. 2006;**23**(17):5117-5128

[29] Boisseau B, Esposito-Farèse G, Polarski D, Starobinsky AA. Reconstruction of a scalar-tensor theory of gravity in an accelerating universe. Physical Review Letters. 2000;**85**(11): 2236-2239

[30] Nojiri S, Odintsov SD, Sáez-Gómez D. Cosmological reconstruction of realistic modified F(R) gravities. Physics Letters B. 2009;**681**(1):74-80

[31] Cognola G, Elizalde E, Odintsov SD, Tretyakov P, Zerbini S. Initial and final

de Sitter universes from modified *f*(*R*) gravity. Physical Review D. 2009;**79**(4): 044001

degeneracy improved by flavor oscillations. Nuclear Physics B. June

[41] Weinberg S. Cosmology. Oxford: Oxford University Press; 2008

*DOI: http://dx.doi.org/10.5772/intechopen.92205*

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies*

Modern Cosmology, Volume 252 of Astronomical Society of the Pacific Conference Series. 2001. p. 355

[49] Knebe A. How to simulate the universe in a computer. Publications of the Astronomical Society of Australia.

[50] Moscardini L, Dolag K, Colpi M, Gorini V, Moschella U. Cosmology with Numerical Simulations. Dordrecht: Springer Netherlands; 2011. pp. 217-237

[51] Dolag K, Borgani S, Schindler S, Diaferio A, Bykov AM. Simulation techniques for cosmological simulations. Space Science Reviews. 2008;**134**(1):

[52] Garca-Farieta JE et al. Clustering and redshift-space distortions in modified gravity models with massive neutrinos. Monthly Notices of the Royal

[53] Planck Collaboration, Ade PAR, et al. Planck 2015 results. xiii.

Cosmological parameters. Astronomy &

[54] Press WH, Schechter P. Formation of galaxies and clusters of galaxies by selfsimilar gravitational condensation. The Astrophysical Journal. 1974;**187**:425-438

[55] Tinker J, Kravtsov AV, Klypin A, Abazajian K, Warren M, Yepes G, et al. Toward a halo mass function for precision cosmology: The limits of universality. The Astrophysical Journal.

[56] Landy SD, Szalay AS. Bias and variance of angular correlation functions. The Astrophysical Journal.

[57] Hamilton AJS. Measuring Omega and the real correlation function from the redshift correlation function. The Astrophysical Journal Letters. 1992;**385**:

Astronomical Society. 2019:07

Astrophysics. 2016;**594**:A13

2008;**688**(2):709-728

1993;**412**:64-71

L5-L8

2005;**22**(3):184189

229-268

[42] Castorina E, Sefusatti E, Sheth RK,

Cosmology with massive neutrinos. II: On the universality of the halo mass function and bias. Journal of Cosmology

Villaescusa-Navarro F, Viel M.

and Astroparticle Physics. 2014;

[43] Seljak U, Slosar A, McDonald P. Cosmological parameters from combining the Lyman-*α* forest with CMB, galaxy clustering and SN constraints. Journal of Cosmology and Astroparticle Physics. 2006;**10**:014

[44] Brandbyge J, Hannestad S, Troels H, Wong YYY. Neutrinos in non-linear structure formation—The effect on halo properties. Journal of Cosmology and Astroparticle Physics. 2010;**2010**(09):

[45] Saito S, Takada M, Taruya A. Impact of massive neutrinos on the nonlinear matter power spectrum. Physical Review Letters. 2008;**100**:

Physics. 2013;**2013**(03):019

distribution as probed by the intergalactic medium. Journal of Cosmology and Astroparticle Physics.

[48] Yepes G. The universe in a

simulations in cosmology. In:

computer: The importance of numerical

Martnez VJ, Trimble V, Pons-Bordera MJ, editors. Historical Development of

2010;**2010**(06):015

**213**

[46] Villaescusa-Navarro F, Bird S, Peña-Garay C, Viel M. Non-linear evolution of the cosmic neutrino background. Journal of Cosmology and Astroparticle

[47] Viel M, Haehnelt MG, Springel V. The effect of neutrinos on the matter

2002;**632**:363-382

**2014**(02):049

014

191301

[32] Narikawa T, Yamamoto K. Characterizing the linear growth rate of cosmological density perturbations in an *f*(*R*) model. Physical Review D. 2010; **85**(11):2236-2239

[33] Tsujikawa S. Matter density perturbations and effective gravitational constant in modified gravity models of dark energy. Physical Review D. 2007; **76**(2):023514

[34] Hu W, Sawicki I. Models of *f*(*r*) cosmic acceleration that evade solar system tests. Physical Review D. 2007; **76**:064004

[35] Wright BS, Winther HA, Koyama K. Cola with massive neutrinos. Journal of Cosmology and Astroparticle Physics. 2017;**2017**(10):054

[36] Giocoli C, Baldi M, Moscardini L. Weak lensing light-cones in modified gravity simulations with and without massive neutrinos. Monthly Notices of the Royal Astronomical Society. 2018; **481**(2):2813-2828

[37] Cowan CL Jr, Reines F, Harrison FB, Kruse HW, McGuire AD. Detection of the free neutrino: A confirmation. Science. 1956;**124**(3212):103-104

[38] Cleveland BT, Daily T, Davis R Jr, Distel JR, Lande K, Lee CK, et al. Measurement of the solar electron neutrino flux with the homestake chlorine detector. Astrophysics Journal. 1998;**496**:505-526

[39] Lesgourgues J, Mangano G, Miele G, Pastor S. Neutrino Cosmology. Cambridge, UK: Cambridge University Press; 2013

[40] Dolgov AD, Hansen SH, Pastor S, Petcov ST, Raffelt GG, Semikoz DV. Cosmological bounds on neutrino

*Massive Neutrinos and Galaxy Clustering in* f*(*R*) Gravity Cosmologies DOI: http://dx.doi.org/10.5772/intechopen.92205*

degeneracy improved by flavor oscillations. Nuclear Physics B. June 2002;**632**:363-382

[21] Baldi M et al. Cosmic degeneracies— I. Joint n-body simulations of modified gravity and massive neutrinos. Monthly Notices of the Royal Astronomical Society. 2014;**440**(1):75-88

*Progress in Fine Particle Plasmas*

de Sitter universes from modified *f*(*R*) gravity. Physical Review D. 2009;**79**(4):

Characterizing the linear growth rate of cosmological density perturbations in an *f*(*R*) model. Physical Review D. 2010;

perturbations and effective gravitational constant in modified gravity models of dark energy. Physical Review D. 2007;

[34] Hu W, Sawicki I. Models of *f*(*r*) cosmic acceleration that evade solar system tests. Physical Review D. 2007;

[35] Wright BS, Winther HA, Koyama K. Cola with massive neutrinos. Journal of Cosmology and Astroparticle Physics.

[36] Giocoli C, Baldi M, Moscardini L. Weak lensing light-cones in modified gravity simulations with and without massive neutrinos. Monthly Notices of the Royal Astronomical Society. 2018;

[37] Cowan CL Jr, Reines F, Harrison FB, Kruse HW, McGuire AD. Detection of the free neutrino: A confirmation. Science. 1956;**124**(3212):103-104

[38] Cleveland BT, Daily T, Davis R Jr, Distel JR, Lande K, Lee CK, et al. Measurement of the solar electron neutrino flux with the homestake chlorine detector. Astrophysics Journal.

[39] Lesgourgues J, Mangano G, Miele G,

Cambridge, UK: Cambridge University

[40] Dolgov AD, Hansen SH, Pastor S, Petcov ST, Raffelt GG, Semikoz DV. Cosmological bounds on neutrino

Pastor S. Neutrino Cosmology.

[32] Narikawa T, Yamamoto K.

[33] Tsujikawa S. Matter density

044001

**85**(11):2236-2239

**76**(2):023514

**76**:064004

2017;**2017**(10):054

**481**(2):2813-2828

1998;**496**:505-526

Press; 2013

[22] Peel A, Pettorino V, Giocoli C, Starck J-L, Baldi M. Breaking

2018;**619**:A38

degeneracies in modified gravity with higher (than 2nd) order weak-lensing statistics. Astronomy & Astrophysics.

[23] Hagstotz S, Costanzi M, Baldi M, Weller J. Joint halo mass function for modified gravity and massive neutrinos.

[24] Lesgourgues J, Pastor S. Neutrino mass from cosmology. Advances in High

[25] Peebles PJE. Principles of Physical Cosmology. Princeton, New Jersey: Princeton University Press; 1993

[26] Dodelson S. Modern Cosmology. San Diego, California: Academic Press; 2003

[27] Mukhanov V. Physical Foundations

[28] Sotiriou TP. *f*(*R*) gravity and scalar tensor theory. Classical and Quantum Gravity. 2006;**23**(17):5117-5128

Reconstruction of a scalar-tensor theory of gravity in an accelerating universe. Physical Review Letters. 2000;**85**(11):

[31] Cognola G, Elizalde E, Odintsov SD, Tretyakov P, Zerbini S. Initial and final

[29] Boisseau B, Esposito-Farèse G, Polarski D, Starobinsky AA.

[30] Nojiri S, Odintsov SD, Sáez-Gómez D. Cosmological reconstruction of realistic modified F(R) gravities. Physics Letters B. 2009;**681**(1):74-80

2236-2239

**212**

of Cosmology. Cambridge, UK: Cambridge University Press; 2005

I: Simulations and cosmological forecasts. arXiv e-prints, arXiv:

1806.07400. June 2018

Energy Physics. 2012;**2012**

[41] Weinberg S. Cosmology. Oxford: Oxford University Press; 2008

[42] Castorina E, Sefusatti E, Sheth RK, Villaescusa-Navarro F, Viel M. Cosmology with massive neutrinos. II: On the universality of the halo mass function and bias. Journal of Cosmology and Astroparticle Physics. 2014; **2014**(02):049

[43] Seljak U, Slosar A, McDonald P. Cosmological parameters from combining the Lyman-*α* forest with CMB, galaxy clustering and SN constraints. Journal of Cosmology and Astroparticle Physics. 2006;**10**:014

[44] Brandbyge J, Hannestad S, Troels H, Wong YYY. Neutrinos in non-linear structure formation—The effect on halo properties. Journal of Cosmology and Astroparticle Physics. 2010;**2010**(09): 014

[45] Saito S, Takada M, Taruya A. Impact of massive neutrinos on the nonlinear matter power spectrum. Physical Review Letters. 2008;**100**: 191301

[46] Villaescusa-Navarro F, Bird S, Peña-Garay C, Viel M. Non-linear evolution of the cosmic neutrino background. Journal of Cosmology and Astroparticle Physics. 2013;**2013**(03):019

[47] Viel M, Haehnelt MG, Springel V. The effect of neutrinos on the matter distribution as probed by the intergalactic medium. Journal of Cosmology and Astroparticle Physics. 2010;**2010**(06):015

[48] Yepes G. The universe in a computer: The importance of numerical simulations in cosmology. In: Martnez VJ, Trimble V, Pons-Bordera MJ, editors. Historical Development of

Modern Cosmology, Volume 252 of Astronomical Society of the Pacific Conference Series. 2001. p. 355

[49] Knebe A. How to simulate the universe in a computer. Publications of the Astronomical Society of Australia. 2005;**22**(3):184189

[50] Moscardini L, Dolag K, Colpi M, Gorini V, Moschella U. Cosmology with Numerical Simulations. Dordrecht: Springer Netherlands; 2011. pp. 217-237

[51] Dolag K, Borgani S, Schindler S, Diaferio A, Bykov AM. Simulation techniques for cosmological simulations. Space Science Reviews. 2008;**134**(1): 229-268

[52] Garca-Farieta JE et al. Clustering and redshift-space distortions in modified gravity models with massive neutrinos. Monthly Notices of the Royal Astronomical Society. 2019:07

[53] Planck Collaboration, Ade PAR, et al. Planck 2015 results. xiii. Cosmological parameters. Astronomy & Astrophysics. 2016;**594**:A13

[54] Press WH, Schechter P. Formation of galaxies and clusters of galaxies by selfsimilar gravitational condensation. The Astrophysical Journal. 1974;**187**:425-438

[55] Tinker J, Kravtsov AV, Klypin A, Abazajian K, Warren M, Yepes G, et al. Toward a halo mass function for precision cosmology: The limits of universality. The Astrophysical Journal. 2008;**688**(2):709-728

[56] Landy SD, Szalay AS. Bias and variance of angular correlation functions. The Astrophysical Journal. 1993;**412**:64-71

[57] Hamilton AJS. Measuring Omega and the real correlation function from the redshift correlation function. The Astrophysical Journal Letters. 1992;**385**: L5-L8

[58] Marulli F, Baldi M, Moscardini L. Clustering and redshift-space distortions in interacting dark energy cosmologies. Monthly Notices of the Royal Astronomical Society. 2012; **420**(3):2377-2386

[59] Tinker JL, Robertson BE, Kravtsov AV, Klypin A, Warren MS, Yepes G, et al. The large-scale bias of dark matter halos: Numerical calibration and model tests. The Astrophysical Journal. 2010;**724**:878-886

[60] Villaescusa-Navarro F et al. Cosmology with massive neutrinos. I: Towards a realistic modeling of the relation between matter, haloes and galaxies. Journal of Cosmology and Astroparticle Physics. 2014;**3**:011

[61] Costanzi M et al. Cosmology with massive neutrinos. III: The halo mass function and an application to galaxy clusters. Journal of Cosmology and Astroparticle Physics. 2013;**2013**(12): 012

[62] Marulli F, Veropalumbo A, Moresco M. CosmoBolognaLib: C++ libraries for cosmological calculations. Astronomy and Computing. 2016;**14**: 35-42

[63] Scoccimarro R, Couchman HMP, Frieman JA. The bispectrum as a signature of gravitational instability in redshift space. The Astrophysics Journal. 1999;**517**(2):531

[64] Scoccimarro R. Redshift-space distortions, pairwise velocities, and nonlinearities. Physical Review D. 2004; **70**(8):083007

[65] Taruya A, Nishimichi T, Saito S. Baryon acoustic oscillations in 2D: Modeling redshift-space power spectrum from perturbation theory. Physical Review D. 2010;**82**(6):063522

[66] Reid BA, White M. Towards an accurate model of the redshift-space clustering of haloes in the quasi-linear regime. Monthly Notices of the Royal Astronomical Society. 2011;**417**(3): 1913-1927

[67] Peacock JA, Dodds SJ. Reconstructing the linear power spectrum of cosmological mass fluctuations. Monthly Notices of the Royal Astronomical Society. 1994;**267**: 1020

[58] Marulli F, Baldi M, Moscardini L.

clustering of haloes in the quasi-linear regime. Monthly Notices of the Royal Astronomical Society. 2011;**417**(3):

[67] Peacock JA, Dodds SJ. Reconstructing the linear power spectrum of cosmological mass fluctuations. Monthly Notices of the Royal Astronomical Society. 1994;**267**:

1913-1927

1020

distortions in interacting dark energy cosmologies. Monthly Notices of the Royal Astronomical Society. 2012;

Kravtsov AV, Klypin A, Warren MS, Yepes G, et al. The large-scale bias of dark matter halos: Numerical calibration and model tests. The Astrophysical

Clustering and redshift-space

*Progress in Fine Particle Plasmas*

[59] Tinker JL, Robertson BE,

Journal. 2010;**724**:878-886

[60] Villaescusa-Navarro F et al. Cosmology with massive neutrinos. I: Towards a realistic modeling of the relation between matter, haloes and galaxies. Journal of Cosmology and Astroparticle Physics. 2014;**3**:011

[61] Costanzi M et al. Cosmology with massive neutrinos. III: The halo mass function and an application to galaxy clusters. Journal of Cosmology and Astroparticle Physics. 2013;**2013**(12):

[62] Marulli F, Veropalumbo A, Moresco M. CosmoBolognaLib: C++ libraries for cosmological calculations. Astronomy and Computing. 2016;**14**:

[63] Scoccimarro R, Couchman HMP, Frieman JA. The bispectrum as a signature of gravitational instability in redshift space. The Astrophysics

[64] Scoccimarro R. Redshift-space distortions, pairwise velocities, and nonlinearities. Physical Review D. 2004;

[65] Taruya A, Nishimichi T, Saito S. Baryon acoustic oscillations in 2D: Modeling redshift-space power spectrum from perturbation theory. Physical Review D. 2010;**82**(6):063522

[66] Reid BA, White M. Towards an accurate model of the redshift-space

Journal. 1999;**517**(2):531

**70**(8):083007

**214**

**420**(3):2377-2386

012

35-42

## *Edited by Tetsu Mieno, Yasuaki Hayashi and Kun Xue*

In the field of plasma physics, plasmas (including charged fine particles) have been actively studied for more than 40 years, and special features of wave phenomena, selforganizations of the particles, potential formations, fluid-like motions of the particles, generations of fine particles in the plasmas, etc. have been investigated. Here, these plasmas are called "fine particle plasmas", which are also called "dusty plasmas" and "complex plasmas". This book intends to provide the reader with the recent progress of studies of fine particle plasmas from the viewpoints of wide and interdisciplinary directions, such as self-organized fine particles, Coulomb crystal formation, behaviors of fine particles, their stability, and syntheses of nano-sized particles in reactive plasmas. Further, the phenomena of dense grain particles and the effects of massive neutrinos in galaxy clustering are included.

Published in London, UK © 2020 IntechOpen © ppart / iStock

Progress in Fine Particle Plasmas

Progress in

Fine Particle Plasmas

*Edited by Tetsu Mieno,* 

*Yasuaki Hayashi and Kun Xue*