On the implementation of the spherical collapse model for dark energy models
Abstract
In this work we review the theory of the spherical collapse model and critically analyse the aspects of the numerical implementation of its fundamental equations. By extending a recent work by [1], we show how different aspects, such as the initial integration time, the definition of constant infinity and the criterion for the extrapolation method (how close the inverse of the overdensity has to be to zero at the collapse time) can lead to an erroneous estimation (a few per mill error which translates to a few percent in the mass function) of the key quantity in the spherical collapse model: the linear critical overdensity δ{sub c}, which plays a crucial role for the mass function of halos. We provide a better recipe to adopt in designing a code suitable to a generic smooth dark energy model and we compare our numerical results with analytic predictions for the EdS and the ΛCDM models. We further discuss the evolution of δ{sub c} for selected classes of dark energy models as a general test of the robustness of our implementation. We finally outline which modifications need to be taken into account to extend the code to more general classes of models, suchmore »
 Authors:
 Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, The University of Manchester, Manchester, M13 9PL (United Kingdom)
 Zentrum für Astronomie der Universität Heidelberg, Institut für theoretische Astrophysik, Philosophenweg 12, D69120, Heidelberg (Germany)
 Publication Date:
 OSTI Identifier:
 22667628
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Cosmology and Astroparticle Physics; Journal Volume: 2017; Journal Issue: 10; Other Information: Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; COMPARATIVE EVALUATIONS; DESIGN; EQUATIONS; EXTRAPOLATION; FORECASTING; GRAVITATIONAL COLLAPSE; MASS; MODIFICATIONS; NONLUMINOUS MATTER; SPHERICAL CONFIGURATION
Citation Formats
Pace, Francesco, Meyer, Sven, and Bartelmann, Matthias, Email: francesco.pace@manchester.ac.uk, Email: sven.meyer@uniheidelberg.de, Email: bartelmann@uniheidelberg.de. On the implementation of the spherical collapse model for dark energy models. United States: N. p., 2017.
Web. doi:10.1088/14757516/2017/10/040.
Pace, Francesco, Meyer, Sven, & Bartelmann, Matthias, Email: francesco.pace@manchester.ac.uk, Email: sven.meyer@uniheidelberg.de, Email: bartelmann@uniheidelberg.de. On the implementation of the spherical collapse model for dark energy models. United States. doi:10.1088/14757516/2017/10/040.
Pace, Francesco, Meyer, Sven, and Bartelmann, Matthias, Email: francesco.pace@manchester.ac.uk, Email: sven.meyer@uniheidelberg.de, Email: bartelmann@uniheidelberg.de. 2017.
"On the implementation of the spherical collapse model for dark energy models". United States.
doi:10.1088/14757516/2017/10/040.
@article{osti_22667628,
title = {On the implementation of the spherical collapse model for dark energy models},
author = {Pace, Francesco and Meyer, Sven and Bartelmann, Matthias, Email: francesco.pace@manchester.ac.uk, Email: sven.meyer@uniheidelberg.de, Email: bartelmann@uniheidelberg.de},
abstractNote = {In this work we review the theory of the spherical collapse model and critically analyse the aspects of the numerical implementation of its fundamental equations. By extending a recent work by [1], we show how different aspects, such as the initial integration time, the definition of constant infinity and the criterion for the extrapolation method (how close the inverse of the overdensity has to be to zero at the collapse time) can lead to an erroneous estimation (a few per mill error which translates to a few percent in the mass function) of the key quantity in the spherical collapse model: the linear critical overdensity δ{sub c}, which plays a crucial role for the mass function of halos. We provide a better recipe to adopt in designing a code suitable to a generic smooth dark energy model and we compare our numerical results with analytic predictions for the EdS and the ΛCDM models. We further discuss the evolution of δ{sub c} for selected classes of dark energy models as a general test of the robustness of our implementation. We finally outline which modifications need to be taken into account to extend the code to more general classes of models, such as clustering dark energy models and nonminimally coupled models.},
doi = {10.1088/14757516/2017/10/040},
journal = {Journal of Cosmology and Astroparticle Physics},
number = 10,
volume = 2017,
place = {United States},
year = 2017,
month =
}

We generalize the small scale dynamics of the Universe by taking into account models with an equation of state which evolves with time, and provide a complete formulation of the cluster virialization attempting to address the nonlinear regime of structure formation. In the context of the current dark energy models, we find that galaxy clusters appear to form at z{approx}12, in agreement with previous studies. Also, we investigate thoroughly the evolution of spherical matter perturbations, as the latter decouple from the background expansion and start to 'turn around' and finally collapse. Within this framework, we find that the concentration parametermore »

Spherical collapse model with nonclustering dark energy
We investigate a spherical overdensity model for the nonclustering dark energy (DE) with the constant equation of state, ω{sub de} in a flat universe. In this case, the exact solution for the evolution of the scale factor is obtained for general ω{sub de}. We also obtain the exact (when ω{sub de} = −1/3) and the approximate (when ω{sub de}≠−1/3) solutions for the ratio of the overdensity radius to its value at the turnaround epoch (y≡R/R{sub ta}) for general cosmological parameters. Also the exact and approximate solutions of the overdensity at the turnaround epoch ζ = (ρ{sub cluster}/ρ{sub m}){sub z} {submore » 
Clarifying spherical collapse in coupled dark energy cosmologies
The spherical collapse model is often used to follow the evolution of overdensities into the nonlinear regime. We describe the correct approach to be used in coupled dark energy cosmologies, where a fifth force, different from gravity and mediated by the dark energy scalar field, influences the collapse. We reformulate the spherical collapse description by deriving it directly from the set of nonlinear hydrodynamical NavierStokes equations. By comparing with the corresponding relativistic equations, we show how the fifth force should be taken into account within the spherical collapse picture and clarify the problems arising when an inhomogeneous scalar field ismore » 
Spherical collapse of dark energy with an arbitrary sound speed
We consider a generic type of dark energy fluid, characterised by a constant equation of state parameter w and sound speed c{sub s}, and investigate the impact of dark energy clustering on cosmic structure formation using the spherical collapse model. Along the way, we also discuss in detail the evolution of dark energy perturbations in the linear regime. We find that the introduction of a finite sound speed into the picture necessarily induces a scaledependence in the dark energy clustering, which in turn affects the dynamics of the spherical collapse in a scaledependent way. As with other, more conventional fluids,more » 
Simple implementation of general dark energy models
We present a formalism for the numerical implementation of general theories of dark energy, combining the computational simplicity of the equation of state for perturbations approach with the generality of the effective field theory approach. An effective fluid description is employed, based on a general action describing singlescalar field models. The formalism is developed from first principles, and constructed keeping the goal of a simple implementation into CAMB in mind. Benefits of this approach include its straightforward implementation, the generality of the underlying theory, the fact that the evolved variables are physical quantities, and that modelindependent phenomenological descriptions may bemore »