Access the full text.

Sign up today, get DeepDyve free for 14 days.

Astrophysics
, Volume 2020 (1607) – Jul 25, 2016

/lp/arxiv-cornell-university/evolution-of-density-and-velocity-profiles-of-matter-in-large-voids-UrIQGTXC9d

- ISSN
- 2227-1481
- eISSN
- ARCH-3330
- DOI
- 10.17721/2227-1481.6.28-33
- Publisher site
- See Article on Publisher Site

We analyse the evolution of cosmological perturbations which leads to the formation of large voids in the distribution of galaxies. We assume that perturbations are spherical and all components of the Universe — radiation, matter and dark energy - are continuous media with ideal ﬂuid energy-momentum tensors, which interact only gravitationally. Equations of the evolution of perturbations in the comoving to cosmological background reference frame for every component are obtained from equations of conservation and Einstein’s ones and are integrated by modiﬁed Euler method. Initial conditions are set at the early stage of evolution in the radiation-dominated epoch, when the scale of perturbation is mush larger than the particle horizon. Results show how the proﬁles of density and velocity of matter in spherical voids with diﬀerent overdensity shells are formed. Key words: cosmology: dark energy, large-scale structure of Universe introduction Large voids in the spatial distribution of galaxies are elements of large scale structure, study of which can give important information about the hidden components of the Universe - dark matter and dark energy. Usually it is assumed, that dark energy is unperturbed in the voids or, at least, impact of its density perturbations on the peculiar motion and spatial distribution of galaxies is negligibly small. In this paper we investigate the evolution of perturbations of density and velocity of matter together with dark energy ones from the early stage, when the scale of initial perturbation is much larger than particle horizon, up to current epoch. We analyse the inﬂuence of dynamical dark energy on such evolution and its dependence on initial conditions. We point attention to the evolution of density and velocity proﬁles of matter during void formation. For this we have developed the program for integrating the system of equation, obtained for description of evolution of spherical perturbation in 3- component medium - radiation, matter and dark energy, - from the equations of relativistic hydrodynamics and gravitation [3]. The component “matter” consists of dark matter (25% of total density) and typical baryonic matter (5%) the dynamics of which is well described by the dust-like medium approach in the large scales. model of spherical void and initial conditions We assume that voids in spatial distribution of galaxies are formed as the result of the evolution of cosmo- logical density perturbations with a negative initial amplitude. It is believed that such perturbations are the result of quantum ﬂuctuations of space-time metric in the inﬂationary epoch. They are randomly distributed in amplitude with normal distribution and are symmetrical by sign of density perturbation from the average in diﬀerent regions of space. We consider only scalar mode of perturbations, in which perturbations of density δ (t, r) and velocity v (t, r) in every component N are correlated because of survival of the growing solution N N only at the stage when the scale of perturbation was larger than particle horizon. Positive perturbations lead to the formation of galaxies and galaxy clusters and negative ones - to the formation of voids. Formation of structures with positive perturbation is well described by Press-Schechter formalism, theory of Gaussian peaks and halo theory of structure formation and their modern modiﬁcations based on the numerical N-body simulations. Although the evolution of voids in the distribution of galaxies is much simpler than evolution of galaxy clusters, since it is described by the quasilinear theory, there is no complete theory of voids formation. Here we analyse the development of negative cosmological density perturbations, which form the voids. The tsizh@astro.franko.lviv.ua arXiv:1607.07404v3 [astro-ph.CO] 30 Sep 2016 Advances in Astronomy and Space Physics mathematical base of their description is the system of 7 diﬀerential equations in partial derivatives for 7 unknown functions of 2 independent variables δ (a, r), δ (a, r), δ (a, r), v (a, r), v (a, r), v (a, r), ν(a, r), de m r de m r which were obtained in [3] (equations (17)-(22)). Here Ω-s denote the mean densities of the components in the unit of the critical one at the current epoch, w ≡ p /ρ is the equation of state parameter of dark de de energy, c is the eﬀective speed of sound of dark energy in its proper frame, H(a) ≡ d ln a/dt is the Hubble parameter, which deﬁnes the rate of the expansion of the Universe and is known function of time for given −4 −3 −3(1+w) cosmology and the model of dark energy (H(a) = H Ω a + Ω a + Ω a ) and H is its today 0 r m de 0 value (Hubble constant). The independent variables are scale factor a and radial comoving coordinate r, which deﬁne the interval in Friedman-Robertson-Walker 4-space: 2 ν(t,r) 2 2 −ν(t,r) 2 2 2 2 2 ds = e dt − a (t)e [dr + r (dθ + sin θdϕ )]. (1) It is assumed that geometry of 3-space of the Universe (unperturbed cosmological background) is Euclidean. The metric function ν(t, r) at the late stages, when the scale of perturbation is much smaller than the particle horizon, is the doubled gravitational potential in the Newtonian approximation of eq. (17) in the paper [3]. The density and 3-velocity perturbations δ and v are deﬁned in coordinates, which are comoving to the N N unperturbed cosmological background (see paragraph 2.2 in [3]). Thus, the velocity perturbation coincide with deﬁnition of peculiar velocity of galaxies (see, for example, [4]). To solve the system of equations (17)-(22) from [3] the initial conditions must be set. Let us relate the initial amplitude of given perturbation with mean-square one given by power spectrum of cosmological perturbations. For this we deﬁne the initial conditions in the early Universe, when ρ ≫ ρ ≫ ρ , and r m de physical size of the perturbation aλ ≫ ct. In that time the perturbations are linear (δ, v, ν ≪ 1), so without loss of generality the solution can be presented in the form of separated variables: ν(a, r) = ν˜(a)f(r), δ (a, r) = δ (a)f(r), v (a, r) = v˜ (a)f (r), (2) N N N N where f(0) = 1 and f (r) ∝ r near the center r = 0. Ordinary diﬀerential equations for amplitudes ν˜(a), δ (a), v˜ (a) are obtained from general system of equations (17)-(22) from [3] by their expansion in Taylor N N series near the center. The analytical solutions of equations for the amplitudes for the radiation-dominated epoch (matter and dark energy can be treated as test components) in the “superhorizon” asymptotic give the simple relation for them: 4 4 C init init init init init init init ˜ ˜ ˜ δ = δ = δ = −ν˜ = C, v˜ = v˜ = v˜ = , (3) r m de r m de 3 3(1 + w) 4a H(a ) init init where C is integration constant, which is deﬁned by initial conditions. We set the value of C in the units of mean-square amplitude of perturbations, which is implied from modern observations. The Planck + HST + WiggleZ + SNLS3 data (see [6] and references therein) tell that amplitude A and spectral index n −1 n of power spectrum of initial perturbations of curvature P (k) = A (k/0.05) are the following [6]: s R s −9 −1 A = 2.224 · 10 , n = 0.963. Since for perturbations with ak ≫ ct the power spectrum perturbations of s s curvature P ≡< ν · ν > is constant in time in the matter- and radiation-dominated epochs, in the range of −5 scales 0.01 ≤ k ≤ 0.1 the initial amplitude which correspondent to mean-square one is: σ ≡ A ≈ 4.7·10 . k s −4 −6 Hereafter we put in our computations C = −1 · 10 ≈ 2σ at a = 10 . init numerical integration For numerical integration of the system of equations (17)-(22) from [3] with initial conditions (3) we have created a computer code npdes.f, which implements the modiﬁed Euler method taking into account the derivatives from the forthcoming step and improving the results by iterations. This scheme of integration is the most resistant to the numerical spurious oscillations, is the fast and precise enough. For example, the Hamming method of prediction and correction of 4-order of precision with 5 iterations at each step need 3 times more processor time for the same precision of ﬁnal result. The step of integration was posed as variable: da = a/N , where number N was picked up so that the numerical precision of the result of integration at a a a = 1 was not worse than 0.1%. In all calculations presented here we took N = 3 · 10 . The numerical derivatives with respect to r in the grid with constant step dr = R /N , where R is m r m radius of spatial region of integration, were evaluated with help of 3-rd order polynomial by method of Savitzky-Golay convolution [5]: y = [3(y − y )/4 − (y − y )/12]/dr. The method was tested by i+1 i−1 i+2 i−2 2 Advances in Astronomy and Space Physics comparing the derivatives of analytical functions of the initial proﬁles of density and velocity perturbations. The value of step dr was estimated so that the diﬀerence between numerical and analytical derivatives do −5 not exceed ∼ 10 of their values. To take into account the Silk damping eﬀect for radiation we have added into equations of evolution of 2 2 δ and v the terms δ k /H/a and v k /H/a accordingly, where the scale of damping k was computed r r r D r D D by formula (10) from [1]. If the values of eﬀective speed of sound in dark energy is c > 0.01c, then the spurious oscillations with growing amplitude appear in this component. Their cause consist in no perfect scheme of integration by time, the numeric derivatives on spatial coordinates and accumulation of numerical errors. To remove them we used the Savitzky-Golay convolution ﬁlter [5] with parameters n = 12, n = 12, m = 6, by which the space- l r dependences of derivatives δ and v˙ were smoothing at each step of integration by a. Such smoothing de de practically does not inﬂuence on the ﬁnal result of integration, which is conﬁrmed by comparison of the results with smoothing and without it for case of the dark energy model with c = 0, for which spurious oscillations do not appear. The maximum diﬀerence is less than 4% for density perturbation and 1% for velocity perturbation of dark energy in the region of maximum amplitude of velocity perturbation. The input parameters of the program are: the Hubble parameter H , the density parameters of all components Ω , Ω , Ω = 1 − Ω − Ω , the equation of state parameter of dark energy w, the speed of r m r de de sound of dark energy c , the initial amplitude of perturbation C, the parameters of proﬁle f(r) of initial perturbation, the parameter of step N in a, the size of integration region R and number of steps of the a m spatial grid N . The computer code npdes.f has been tested by comparison of the results of the integration by code with 1) known analytical solutions for density and velocity perturbations in conformal-Newtonian frame for radiation- and matter-dominated Universes [8], 2) results of integration of linear perturbation by CAMB code [2] and 3) results of integration by dedmhalo.f code [3], developed on the basis of dverk.f for perturbation in the central region of the spherical perturbation. In all cases deviations did not exceed a few tenths of a percent, which means, that precision of the integration is better then 1 %, and hence is high enough for our studies. formation of voids in the cosmological models with dark energy From our previous studies and studies of other authors we know that the values of density parameter and the equation of state parameter of dark energy are well constrained by current observational data, while the value of eﬀective speed of sound of dark energy is not constrained (see, for example, [6] and citation therein). That is why in this work we analyse the formation of voids in the cosmological models with dark energy with Ω = 0.7, w = −0.9 and diﬀerent values of c ∈ [0, 1]. Other cosmological parameters in computations are de s −5 ﬁxed too: Ω = 4.17 · 10 , Ω = 0.3 − Ω , H = 70 km/s·Mpc. r m r 0 2 −βr In this work we study the formation of the spherical voids with initial proﬁle f(r) = (1−αr )e , where α gives the size of the void r = 1/ α and β deﬁnes the initial amplitude of shell overdensity around the −1 −1−β/α void: δ = −αβ Ce . For comparison of the results of this paper with the results of accompanying one [7], let us set α = (k/π) and β = 3α/4. This is a proto-void, which is surrounded with overdensity shell with δ ≈ δ(r = 0)/8. For comparison we will also analyse the evolution of the void with shells with smaller amplitudes of overdensity in 2 and 4 times. −1 In ﬁg. 1 we show the formation of the spherical void with r = 31.4 Mpc (k = 0.1 Mpc ) in the matter and dark energy with c = 0: δ (a , r) and v (a , r) for a = a , ..., a = 1. Black lines denote m,de i m,de i i init 30 the initial proﬁles of density and velocity perturbations of both components, red lines denote the ﬁnal ones. The ﬁgure on the right depicts the evolution of absolute values of amplitudes of perturbations in the central point of spherical void (matter - thick solid line, dark energy - thick dashed line) and in the overdensity shell (matter - thin solid line, dark energy - thin dashed line). Velocity perturbation (central panel) are given for the ﬁrst maximum (thick lines) and ﬁrst minimum (thin lines). Dotted line denotes the Silk damping for the radiation component. One can see, that in this dark energy model the perturbations of matter and dark energy grow monotonically after entering the horizon: the black lines are internal, the red lines are external. We also note, that the amplitude of the density perturbation of dark energy is approximately 40 times smaller than the matter one. The values of velocity perturbations of matter and dark energy in this model of dark energy are the same throughout the evolution of the void. They increase monotonically from a to a ≈ 0.56. It is easy to see that the latter value corresponds to the moment of change from decelerated init http://camb.info http://www.cs.toronto.edu/NA/dverk.f.gz 3 Advances in Astronomy and Space Physics Figure 1: Void formation in dark matter (left column) and dark energy with c = 0 (central column). On the right - evolution of absolute values of amplitudes of density (top panel) and velocity (bottom panel) perturbations; solid lines in the top and middle panels - dark matter, dashed lines - dark energy, point lines - radiation; thick black solid and dashed lines - for central point, green ones - for overdensity shell. Figure 2: Void formation in dark matter (left column) and dark energy with c = 0.1 (central column). On the right - evolution of absolute values of amplitudes of density (top panel), velocity (bottom panel) perturbations in the central point (thick lines) and overdensity shell (thin lines). 4 Advances in Astronomy and Space Physics −2 Figure 3: Evolution of proﬁles of matter density and velocity perturbations with initial parameters α = (31.4) −2 Mpc and β = 3α/4 (left), β = 3α/2 (central), β = 3α (right). Thick black line is initial proﬁle, thick red one is ﬁnal proﬁle. expansion of the Universe to the accelerated one. The evolution of the absolute values of density and velocity perturbations of matter and dark energy in the overdensity shell is similar to those in the center. Similar results of modelling of the void formation in the matter and dark energy with c = 0.1 are shown in ﬁg. 2. “The picture” of the evolution of the matter density and velocity perturbations has not changed, while for dark energy it has changed drastically. The ﬁnal proﬁles of dark energy perturbations (red lines) are lying on the zero line now. The right ﬁgure explains such behaviour of dark energy during the void formation: the velocity perturbation after the entering into horizon decrease quickly, and density perturbation slightly changes during all stages and in the current epoch doesn’t diﬀer practically from the background value: −5 δ (1, 0) ≈ −2· 10 . The matter density perturbation in the central part of this void at the current epoch is de δ (1, 0) ≈ −0.7. We see also that the evolution of the absolute values of density and velocity perturbations of dark energy in the overdensity shell slightly diﬀer from the evolution of ones in the center of the void. The perturbation of dark energy with larger values of eﬀective speed of sound after entering the particle horizon is smoothed out even faster. Therefore, the ratio of densities of dark energy and matter in the center of the void is ρ (1, 0) 1 + δ (1, 0) Ω de de de = , ρ (1, 0) 1 + δ (1, 0) Ω m m m and in the case of evolution with considered initial condition this ratio is 3 time larger than on cosmological background. This points to the importance of studying of the voids for establishing the nature of dark energy. Study of the evolution of spatial proﬁles of matter density and velocity perturbations is important for understanding of the formation of voids. They can be obtained by normalization of every curve in the left columns in ﬁg. 1-2 by its amplitude. The result is given in ﬁgs. 3 for initial proﬁles with three diﬀerent values of parameter β. They show, that in our model the perturbation with initial density proﬁle −4 2 −βr −2 δ (r) = −1·10 [1−(r/r ) ]e with r = 31.4 Mpc and β = 3r /4 (left column) leads to the formation init v v of void with radius in comoving coordinates ≈ 38 Mpc with amplitude of density perturbation in the center δ (1, 0) ≈ −0.68 and the overdensity shell around it (δ > 0) with thickness ≈ 36 Mpc and amplitude m m −2 of density perturbation δ ≈ 0.33. In the case of β = 3r /2 (central column) the radius of the central void is ≈ 31 Mpc, the amplitude of density perturbation in the center is δ (1, 0) ≈ −0.69, the shell of overdensity has thickness ≈ 30 Mpc and the amplitude of density perturbation δ ≈ 0.16. In the case of −2 β = 3r (right column) the void has the following parameters: the radius ≈ 25 Mpc, the amplitude of 5 Advances in Astronomy and Space Physics density perturbation in the center δ (1, 0) ≈ −0.73, the maximum of overdensity in the shell δ ≈ 0.09 is at m e distance ≈ 31 Mpc from the center, external bound of shell, where sign of perturbation changes from ”+“ to ”-“ is absent. Important is the dependence of proﬁles of peculiar velocity of matter in the void and around it on the model parameters and initial perturbation. From ﬁg. 1-3 one can see, that the ﬁrst positive peak of peculiar velocity (from the center) is at the edge of the void and the second negative (velocity towards −2 the center) is at the edge of the overdensity shell. The values of velocities in the case of β = 3r /4 are −2 v (30 Mpc) ≈ 200 km/s, v (65 Mpc) ≈ −70 km/s, in the case of β = 3r /2: v (25 Mpc) ≈ 170 m−v m−e m−v −2 km/s, v (60 Mpc) ≈ −66 km/s, and in the case of β = 3r : v (20 Mpc) ≈ 150 km/s and the second m−e m−v negative peak is absent. Note, that ﬁnal value of the amplitude of the perturbation in the shell is the greater, the greater is its initial value (smaller value of β), for the same value of initial amplitude in the center. One can see also, that overdensity shell appears in the process of evolution of void even if its amplitude was very small in the initial proﬁle (ﬁgure on the right), or absence at all (α = 0, Gaussian initial proﬁle). The evolution of matter density and velocity proﬁles points that for interpretation of the observational data on the distribution of void galaxies in the phase space the non-linear theory should be used (see also table in [7]). conclusion The large voids in the spatial distribution of galaxies are formed from the negative cosmological density perturbations of matter. The amplitude of the density perturbation in the central part of the void at the current epoch is deﬁned by the depth of dip of Gaussian ﬁeld of the initial matter density perturbations, the parameters of the cosmological model and parameters of the initial proﬁle. For example, in the cosmological model with quintessential dark energy the initial negative density perturbation with proﬁles similar to the Gaussian ones lead to the formation of the voids with the overdensity shells. In such voids with r ≈ 30 Mpc and δ (1, 0) ≈ −0.7 the maximal values of the peculiar velocity of galaxies are ∼ 150 − 200 km/s (movement from the center in the comoving coordinates) and are reached near the boundary. In the shells such velocity is directed to the center, however its value does not exceed ∼ 70 km/s. With increasing the parameter of initial proﬁle β for the same r the amplitudes of the density and velocity perturbations in the shells decrease. The density and velocity perturbations of the dark energy evolve similarly to the perturbations of matter at the stage when their scales are much larger than the particle horizon. After they enter the particle horizon their evolution depends on the value of the eﬀective speed of sound c . If c = 0, then similarity is conserved s s with the diﬀerence that the amplitude of density perturbation of dark energy is smaller in factor 1 + w. At the later epoch, when the dark energy density dominates, this diﬀerence increased yet in ≈ 4− 5 times more. If 0 < c ≤ 1, then the amplitude of velocity perturbation of dark energy after entering the horizon decreases rapidly, the amplitude of the density perturbation doesn’t increase or even decreases too. Therefore, in the voids the density of quintessential dark energy is approximately the same as in cosmological background. The ratio of the densities of dark energy and matter is in 1/(1 + δ ) larger than in the cosmological background. The more hollow void is the larger this ratio is. That is why the large voids are important elements of large-scale structure of the Universe for testing models of dark energy and gravity modiﬁcations. References [1] Hu W. & Sugiyama N. Anisotropies in the cosmic microwave background: an analytic approach. 1995, Astrophys. J. 444, [2] Lewis A., Challinor A. & Lasenby A. Eﬃcient Computation of Cosmic Microwave Background Anisotropies in Closed Friedmann-Robertson-Walker Models. 2000, Astrophys. J. 538, 473; http://camb.info. [3] Novosyadlyj B., Tsizh M. & Kulinich Yu. Dynamics of minimally coupled dark energy in spherical halos of dark matter. 2016, Gen. Relat. Grav. 48, 3. [4] Peebles P. J. E. The large scale structure of the Universe. 1980, Princeton University Press, Princeton, 422 p., ISBN:9780691082400 [5] Savitzky A. & Golay M.J.E. Smoothing and Diﬀerentiation of Data by Simpliﬁed Least Squares Procedures. 1964, Analytical Chemistry, 36(8), 1627-39. [6] Sergijenko O. & Novosyadlyj B. Sound speed of scalar ﬁeld dark energy: weak eﬀects and large uncertainties. 2015, Phys. Rev. D 91, 083007. [7] Tsizh M. & Novosyadlyj B. Formation of the spherical voids in models of Universe with dark energy. 2016, Visnyk Kyivsogo Universytety. Astronomia, 53, 32. 6 Advances in Astronomy and Space Physics [8] Novosyadlyj B. Formation of the large scale structure of the Universe: theory and observations. 2007, J.Phys.Stud. 11, 226-257.

Astrophysics – arXiv (Cornell University)

**Published: ** Jul 25, 2016

Loading...

You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!

Read and print from thousands of top scholarly journals.

System error. Please try again!

Already have an account? Log in

Bookmark this article. You can see your Bookmarks on your DeepDyve Library.

To save an article, **log in** first, or **sign up** for a DeepDyve account if you don’t already have one.

Copy and paste the desired citation format or use the link below to download a file formatted for EndNote

Access the full text.

Sign up today, get DeepDyve free for 14 days.

All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.