J Phys Astr, Volume: 8( 3) DOI: 10.37532/2320-6756.2020.8(3).194
Study of Cosmological Bias of Galaxies in Quasi Linear Regime
- *Correspondence author :
- Sujata Mohanty
Research Scholar, CMR University, Bangalore-43, India
E-Mail: [email protected]
Received: May 13, 2020; Accepted: June 18, 2020; Published: June 27, 2020
Citation: Sujata M, Harish R. Study of Cosmological Bias of Galaxies in Quasi Linear Regime. J Phys Astron. 2020;8(3):194.
Power spectrum of density variation of matter dominated era is calculated to prove that structure formation in the universe starts in the quasi-linear regime. Dimensionless power spectra for 1-halo and 2-halo for baryonic matter are calculated and compared with the 1-halo and 2-halo power spectra for dark matter. The higher values obtained for baryonic matter indicates over dense regions of baryonic matter leading to collapse to form large scale structures. The power spectrum values lead to cosmological bias values in unperturbed and perturbed universe. The obtained bias values being approaching unity indicate formation of large scale structure of galaxies starts in the quasi-linear regimes
Cosmological bias; Large scale structure; Quasi linear regime; Power spectrum; Dark matter; Baryonic matter
The concept of galaxies being biased tracers with respect to dark matter was introduced by Kaiser Nick . Bias, which depends on the length scale, refers to the ratio of the baryonic matter density to the underlying dark matter density. Dark matter is invisible, as it does not emit radiation and its presence is inferred by its gravitational effect. In order to study bias, knowledge of the inhomogeneity in density is required. Models such as the spherical collapse model, halo model, local bias model and linear and non-linear perturbation models are required to obtain density contrast of baryonic matter and dark matter. On large scale, the universe is homogeneous and isotropic and hence does not lead to bias. But bias is prominent in small scale. Gravity equally affect both baryonic matter and dark matter. Baryonic matter is acted upon by some other forces like radiation pressure, magnetic field etc. While dark matter only affected by gravity and behaves like a potential well inside the halo into which the baryonic matter is going on accumulating to form small structures first further leading to formation of large scale structures. So the distribution of galaxies and the dark matter within a given halo can be used to estimate the statistical properties of the dark matter density in the quasi-linear regime, where the density fluctuation (δ) approaches unity. These galaxies are the biased tracers of the dark matter in the halo. Under these conditions the matter is in virial equilibrium, where the desired physical properties of the halo can be estimated . The spherical collapse model was first studied by Gunn and Gott . This is the simplest non-trivial model for the way an object like galaxy or cluster of galaxies breaks away from the general expansion. In the model, the universe is spherically symmetric about one spot and the matter is an ideal fluid with zero pressure. At the initial epoch everything is expanding smoothly as the universe is expanding though gravity acts there. In this paper, dimensionless power spectra for 1-halo and 2-halo for baryonic matter are calculated and compared with the 1-halo and 2-halo power spectra for dark matter. The higher values obtained for baryonic matter indicate over dense regions of baryonic matter leading to form large scale structures. The power spectrum values lead to cosmological bias values in unperturbed and perturbed universe. The obtained bias value being approaching unity indicates the start of formation of large scale structures of galaxies in the quasi-linear regime.
Analysis of halo power spectrum
Fry and Gatzanaga  assume that the density contrast in halo distribution, dh; can be expressedas a non-linear function of the local density contrast of dark matter (DM) ,dm:
In this section we evaluate the bias values (b1,b2,b3) under certain assumptions, which follow and show using these values that in the quasi-linear regime, the clustering of galaxies is more due to many forces like gravity, magnetic field, radiation pressure etc., than clustering of dark matter which is only affected by gravity.
We assume an initial spectrum of the density fluctuation field as
In accordance with , where the authors have considered this spectrum as a possible shape of the power spectrum on cluster-like scales in Cold Dark Matter (CDM) models. In terms of variance s2 (R) on scale R, the dimension less power spectrum as per 
Where, is the critical density for an object to collapse.
We assume that the density around a virialized halo is given by 
Where is the core radius of the halo profile, is the halo concentration parameter,Δnl is the non-linear dimensionless power spectrum. For which the normalized Fourier transform is
Where ? is a numerical coefficient independent of scale. The dimensionless 1-halo and 2-halo power spectrum for massive and less concentrated halo, for , are given by
B(k) is known as bi-spectrum which is given by
Where is the critical density for an object to collapse. Dimensionless power spectrum for 1-halo term and 2-halo term for more massive and more concentrated halo have been derived analytically  which are as follows:
Equation-4 and 5 explains the variation of concentration parameter with halo mass (m). Which further contributes to the total power. At small scale most of the power comes from 2-halo terms. At large k 1-halo term dominates the power.
Where is called as kappa and is given by
Kappa is a distance in the unit of scale radius. The non-linear dimensionless power spectrum Assuming halos to be virialized at half the turn around time R and m are the initial size and mass of the halo. R*and m* are the characteristic size and mass of the halo respectively [5-7]. In Einstein de-Sitter Cosmology, and for initial power spectrum at the turn around epoch, since
We have Where defines the characteristic mas of the halo. Now straightforwardly we have
Where is a numerical coefficient independent of scale. The value of or more massive and more concentrated halo. From equation 6,Kappa .The variance is given by 
The value of A can be obtained as
from equation 3 can evaluated by using as using equation-4 and using equation-5. We have
The density contrast in halo matter power spectrum is then
Analysis of dark matter power spectrum
For an Einstein de-sitter universe the solution to the density field
can be expressed in terms of the Fourier transform of the initial fluctuation 
The window function is and the delta function For R=8h-1 MPC,k=1h MPC-1, The first order density spectrum for matter is The dimensionless power spectrum for matter is
The unit of power spectrum is (h−1MPC)3.
The second order density field of matter is where
Now 2nd order density power spectrum for matter is and hence dimensionless power spectrum for matter is
Calculation of bias
It is the square root of ratio of halo matter power spectrum to the dark matter power spectrum.
Using equation (1), b0can evaluated as b0=0.809 (28)
We calculated bias in small scales as most of the power comes from both the halo terms (TABLE 1 and FIG.1).
Table 1: [Δ_h]2 represents dimensionless power spectrum of Halo profile and [Δ_m]2 represents dimensionless matter power Spectrum, b represents bias values.
Figure 1: It represents the graph of power spectrum versus bias. It compares the value of bias of power spectrum of halo profile and that of dark matter.
Result and Discussion
The values of bias has been studied analytically and found asb0=0.809,b1=9.639, b2=12.916. The value of b0clearly indicates that clustering of baryonic matter and clustering of dark matter are almost equal. That means this is the quasi-linear regime where large scale structure formation starts and being transformed from linear to nonlinear regime which is indicated by the increased value of bias valuesb1and b2 From the graph one can observe that with increase of bias, the halo matter power spectrum increases. But the dark matter power spectrum decreases with increase of bias. That means the clustering of galaxies are more under the action of different forces including gravity than the clustering of dark matter under the force of gravity only to form the large scale structures .
We have calculated the bias valuesb0, b1,b2 for baryonic matter and dark matter using the 1-halo and 2-halo power spectra. From our study, we find that the dimensionless power spectrum for baryonic matter increases at a faster rate than that for dark matter. In equation (1), in Taylor expansion of density perturbation the cosmological bias values increases. The lowest value b0 approaching 1, indicates the quasi- linear regime where the baryonic matter (under the action of gravity as well as radiation force and magnetic field etc) accumulates into the potential well of dark matter (affected by gravity only) to form small structures first. Which further leads to large scale structures we see today. The largest values of bias indicate the non-linearity domain where the formation of large scale structure formation is more complex and can’t studied by analytic method. It can be obtained by simulation only.
- Nick K. On the spatial correlations of Abel clusters. Astrophysics. 1984;284:12.
- Press WH, Schechter P. Formation of galaxies and clusters of galaxies by selfsimilar gravitational condensation. APJ. 1974;187:425.
- James EG, Richard GJ. On the in-fall of matter into clusters of galaxies and some effects on their evolution. APJ. 1972;176:1-19.
- Asantha C, Ravi S. Halo models of large scale structures. Astrophysic. 2002;508:1-129.
- Bernardeau F. The effects of smoothing on the statistical properties of large scale cosmic fields. A and A.1994;291:697B.
- arxiv:1402.7073v2 [astro-ph.co], 2014.
- Fry JN, Gatzanaga E. Biasing and hierarchical structures in large-scale structures. APJ. 1993;413:47.
- Manera M, Gaztanaga E. The local bias model the large scale halo distribution. RAS. 2011;415:383-98.
- Julio FN, Carlos SF, Simon DMW. The structure of cold dark matter halos, APJ. 1996;462:563.
- Peacock JA. Cambridge University Press. 2012.
- Roman S, Joshua F. Loop corrections in non linear cosmological perturbation Theory. Astrophysical Journal. 1995;473:620.
- Sheath RK, Torment G. MNRAS, 199:323.