Hydrodynamical simulations of coupled and uncoupled quintessence models I: Halo properties and the cosmic web
Abstract
We present the results of a series of adiabatic hydrodynamical simulations of several quintessence models (both with a free and an interacting scalar field) in comparison to a standard CDM cosmology. For each we use particles in a periodic box assuming WMAP7 cosmology. In this work we focus on the properties of haloes in the cosmic web at . The web is classified into voids, sheets, filaments and knots depending on the eigenvalues of the velocity shear tensor, which are an excellent proxy for the underlying overdensity distribution. We find that the properties of objects classified according to their surrounding environment shows a substantial dependence on the underlying cosmology; for example, while shows average deviations of per cent across the different models when considering the full halo sample, comparing objects classified according to their environment, the size of the deviation can be as large as per cent.
We also find that halo spin parameters are positively correlated to the coupling, whereas halo concentrations show the opposite behaviour. Furthermore, when studying the concentrationmass relation in different environments, we find that in all cosmologies underdense regions have a larger normalization and a shallower slope. While this behaviour is found to characterize all the models, differences in the bestfit relations are enhanced in (coupled) dark energy models, thus providing a clearer prediction for this class of models.
keywords:
methods:body simulations – galaxies: haloes – cosmology: theory – dark matter1 Introduction
Over more than 15 years, since observations of highredshift Supernovae of type Ia (see Riess:1999; Perlmutter:1999) first indicated that the Universe is undergoing an accelerated expansion, a large number of cosmological probes, including cosmic microwave background (CMB) anisotropies (Wmap:2011; Sherwin:2011), weak lensing (Huterer:2010), baryon acoustic oscillations (BAO) (Beutler:2011) and large scale structure (LSS) surveys (SDSS:2009), have confirmed this startling claim and shown that the Universe is spatially flat. To explain these diverse observation, cosmology requires the presence of a fluid, called dark energy (DE), which permeates the whole Universe and exerts a negative pressure, eventually overcoming the gravitational pull that would otherwise dominate. The standard model of cosmology, referred to as CDM, provides the simplest possible explanation for DE, assuming that DE is played by a constant called which possesses a constant equation of state, such that . However, despite its simplicity and observational viability, CDM still lacks of appeal from a purely theoretical point of view, due to fine tuning and coincidence problems; the first refers to the fact that, if we assume that is the zeropoint energy of a fundamental quantum field, to be compatible with the aforementioned cosmological constraints its density requires an unnatural finetuning of several tens of orders of magnitude. The second problem arises from the difficulty in explaining in a satisfactory way the fact that matter and dark energy densities today have comparable values, although throughout most of the cosmic history their evolutions have followed completely different patterns.
It is thus natural to explore the possibility that dark energy does not take the form of a cosmological constant, , but is instead a dynamical component of the universe, whose energy density evolves with time, eventually dominating in the present epoch. In this sense, a large number of different models, such as Chaplygin gas (Kamenshchik:2001), vector dark energy (BeltranMaroto:2008; Carlesi:2012), essence (Picon:2000) and quintessence (Wetterich:1995; Caldwell:1998; Copeland:1998; Zlatev:1999) have been proposed to overcome the perceived theoretical shortcomings of the standard cosmology. In particular, quintessence (or scalar field) models are viable and likely candidates for dynamical dark energy (see Tsujikawa:2013), as they can reproduce current observational data without being plagued by the finetuning problem of CDM, since their expansion history  at least for a set of different potentials  is almost insensitive to the particular choice of the field’s initial conditions. An interesting subset of quintessence theories is represented by coupled models, where it is assumed that the scalar field has a nonnegligible interaction to the dark matter sector (Amendola:2000) and is thus expected to leave a strong imprint on structure formation.
While both classes of quintessence models have been already studied numerically by means of body simulations (see for instance Klypin:2003, DeBoni:2011 for free and Baldi:2010a, Li:2011 for coupled models), in the present work we aim at investigating and highlighting the differences arising among coupled and uncoupled scalar fields with the same potential, and compare our results to a benchmark CDM cosmology. Our aim is to disentangle the effects due to the fifth force acting on dark matter particles from those caused to the dynamical nature of dark energy, when considering the deeply nonlinear regime of the models. This way we will discern strategies to observationally distinguish between coupled and uncoupled forms of quintessence and thus to provide new tools for model selection.
Using a suitably modified version of the publicly available SPH/body code GADGET2 (Springel:2005) we undertake a series of simulations of different quintessence models with a RatraPeebles (Ratra:1988) potential and several values of the coupling parameter allowed by current observational constraints. The box size (), the number of particles () and the use of adiabatic smoothed particle hydrodynamics allow us to analyse a large amount of different properties with a good resolution and statistics. In this first of a series of papers, we will consider largescale structures (LSS) and its environment, with the physics of galaxy clusters presented in a followup paper (Paper II). In the present work, we analyze in particular the structure of the cosmic web and the correlations between the environment, dark matter haloes and gas across these different cosmological models.
The paper is organized as follows; in Section 2, we briefly recall the general features of the quintessence models considered in this work as well as of the recipes necessary for their simulation using body techniques. In Section 3, we discuss the settings of our particular simulations as well as those of the halo finder, together with the classification of the cosmic web. In Section 4 we present LSS properties of the modified frameworks, Section 5 is dedicated to the general features of the cosmic web while in Section 6 we describe the results of the correlation of haloes to their environment. A summary of the results obtained and a discussion on their implications is then presented in section Section 7.
2 Prerequisites
Here, we will briefly recall the basic properties of quintessence models and their implementation into a cosmological body algorithm. We refer the reader to the works of Wetterich:1995; Amendola:2000; Amendola:2003; Pettorino:2012; Chiba:2013 for discussions on the theoretical and observational properties of (coupled) quintessence models, and to Maccio:2004; Baldi:2010a; Li:2011d for a thorough description of the numerical approaches.
2.1 The models
In quintessence models the role of dark energy is played by a cosmological scalar field whose Lagrangian can be generally written as:
(1) 
where in principle can interact with the dark matter field through its mass term, meaning that, in general, dark matter particles will have a timevarying mass. With a suitable choice of the potential , quintessence cosmologies can account for the late time accelerated expansion of the universe both in the interacting and non interacting case. In the present work we have focused on the so called RatraPeebles (see Ratra:1988) self interaction potential:
(2) 
where is the Planck mass, while and are two constants whose values can be fixed by fitting the model to observational data (see Wang:2012; Chiba:2013).
In Eq. (1), we allowed for the scalar field to interact with matter through the mass term ; a popular choice (Pettorino:2012, see) for the function is:
(3) 
which is also the one assumed in this paper.
In the following, we have taken into account a constant interaction term , which from Eq. (3) implies an energy flow from the dark matter to the dark energy sector and thus a diminishing mass for dark matter particles.
Model 


uDE  
cDE033  
cDE066  
cDE099 
In Table 1 we list the values for , and of Eq. (2) and Eq. (3) as used in the four nonstandard cosmologies under investigation  an uncoupled Dark Energy (uDE) model and three coupled Dark Energy (cDE) ones. The latter differ only by the choice of the coupling and have been named accordingly. The particular values used in all the implementations have been selected according to the Cosmic Microwave Background (CMB) constraints discussed in Pettorino:2012, to ensure the cosmologies under investigation to be compatible with the WMAP7 dataset (wmap7). However, more recent results obtained using Planck data (see Pettorino:2013), provide even tighter constraints on the free parameters of these models, which shall be the object of subsequent investigation.
2.2 Numerical implementation
Parameter  Value 

The first simulations of interacting dark energy models were performed by Maccio:2004, who described the basic steps for implementing interacting quintessence into the ART code. In our case, we built our implementation on PGADGET2, a modified version of the publicly available code GADGET2 (Springel:2001; Springel:2005). This version has first been developed to simulate vector dark energy models (see Carlesi:2011; Carlesi:2012) and was then extended to generic dynamical dark energy as well as coupled dark energy cosmologies. The algorithm used is based on the standard TreePM solver with some modifications added to take into account the additional longrange interactions due to the coupled scalar field which effectively act as a rescaling of the gravitational constant. For the implementation of these features nonstandard models we followed closely the recipe described in Baldi:2010a, to which the reader is referred.
This approach requires that a number of quantities, namely:

the full evolution of the scalar field and its derivative ,

the variation mass of cold dark matter particles , and

the background expansion .
have to be computed in advance and then interpolated at run time. We therefore implemented background and first order Newtonian perturbation equations into the publicly available Boltzmann code CMBEASY (Doran:2005) to generate the tables containing the aforementioned quantities. The starting background densities were chosen in order to ensure the same values at for the cosmological parameters listed in Table 2; linear perturbations have been solved assuming adiabatic initial conditions.
Finally, in the case of nonstandard cosmologies it is necessary to properly generate the initial conditions of the body simulations taking into account not only the different matter power spectra but also the altered growth factors and logarithmic growth rates, respectively. These are in fact the necessary ingredients to compute the initial particles’ displacements and velocities on a uniform Cartesian grid using the first order Zel’dovich approximation (Zeldovich:1970). We implemented these changes into the publicly available NGenic^{1}^{1}1http://www.mpagarching.mpg.de/gadget/ MPI code, which is suitable for generating GADGET format initial conditions. Again, the matter power spectra, growth factors and logarithmic growth rates have been computed for the four nonstandard cosmologies using the modified CMBEASY package.
All the above changes have been carefully tested against theoretical predictions and the previous results existing in the literature to ensure the consistency and reliability of our modifications.
3 The simulations
3.1 Settings
Parameter  Value 

Our set of body simulation has been devised in order to allow us to compare and quantitatively study the peculiarities of the different models in the physics of galaxy clusters and the properties of the cosmic web. To do this, we have chosen a box of side length (comoving) where we expect to be able to analyze with adequate resolution a statistically significant () number of galaxy clusters () as well as the properties of the different cosmic environments, classified as voids, sheets, filaments and knots. The parameters chosen to set up the simulations, which are common to all the six models under investigation, are listed in Table 3.
In this series of simulations we implemented adiabatic SPH only, thus neglecting the effects of all sources of radiative effects (Monaghan:1992; Springel:2010) This way we are able to establish a clear basis for the differences induced on baryons by the different cosmologies, without the need to take into account the additional layer of complexity introduced of radiative physics, which in itself requires a substantial degree of modeling. The publicly available version of GADGET2 performs a Lagrangian sampling of the continuous fluid quantities using a set of discrete tracer particles. Gas dynamics equations are then solved using the SPH entropy conservation scheme described in Springel:2005. In our case, continuous fluid quantities are computed using a number of smoothing neighbours . Gas pressure and density are related through the relation , where under the adiabatic assumption.
In addition to the four quintessence models (whose parameters have been given in Table 1) we simulated a CDM cosmology, which we use as a benchmark to pinpoint deviations from the standard paradigm. The initial conditions for all the simulations have been generated using the same random phase realization for the Gaussian fluctuations, which enables us to consistently crosscorrelate properties enforcing the same values at present time for and (cf. Table 2), across different simulations.
3.2 Halo identification
We identified haloes in our simulation using the open source halo finder AHF^{2}^{2}2AHF stands for Amiga Halo Finder, which can be downloaded freely from http://www.popia.ft.uam.es/AHF described in Knollmann:2009; this code improves the MHF halo finder (Gill:2004) and has been widely compared to a large number of alternative halo finding methods (Knebe:2011; Onions:2012; Knebe:2013). AHF computes the density field and locates the prospective halo centres at the local overdensities. For each of these density peaks, it determines the gravitationally bound particles, retaining only peaks with at least 20 of them, which are then considered as haloes and further analyzed.
The mass is computed via the equation
(4) 
so that is defined as the total mass contained within a radius at which the halo matter overdensity reaches times the critical value . Since the critical density of the universe is a function of redshift, we must be careful when considering its definition, which reads
(5) 
as the evolution of the Hubble parameter, differs at all redshifts in the five models. In the latest version of AHF this problem is solved reading the for the cDE and uDE models in from a precomputed table, which then allows to compute the consistently in each case. For all models we assume .
3.3 Classification of the cosmic web
As we intend to correlate halo properties with the environment, it is necessary to introduce the algorithm used for the classification of the cosmic web into voids, sheets, filaments and knots. Using the term cosmic web (Bond:1996) we refer to the complex visual appearance of the large scale structure of the universe, characterized by thin linear filaments and compact knots crossing regions of very low density (Massey:2007; Kitaura:2009; Jasche:2010).
The exact mathematical formulation for describing the visual impression of the web is highly nontrivial and can be implemented using two different approaches, the geometric one and the dynamic one. The first one relies on the spatial distribution of haloes in simulations (Novikov:2006; AragonCalvo:2007) disregarding the dynamical context. The second approach starts with the classification of Hahn:2007, where they identified the type of environment using the eigenvalues of the tidal tensor (i.e. the Hessian of the gravitational potential), rather than studying the matter density distribution.
However, these particular approaches are unable to resolve the web on scales smaller than a few megaparsecs (ForeroRomero:2009). While retaining the original idea of dynamical classification, Hoffman:2012 proposed to replace the tidal tensor with the velocity shear, showing that this approach has a much finer resolution on the smaller scales while reproducing the large scale results of the other approach. Defining the velocity shear tensor as
(6) 
and diagonalizing it, we obtain the eigenvalues and . Taking the trace of we obtain
(7) 
from which we see that there is indeed a direct relationship between the eigenvalues of the velocity shear tensor and the matter overdensity. In practice the eigenvalue is related to the intensity of the inflow (outflow) of matter along the th axis in the base where is diagonal.
We therefore proceed to classify the cosmic web ordering the eigenvalues and defining the different points on the web as (Hoffman:2012; Libeskind:2012; Libeskind:2013):

voids, if

sheets, if

filaments, if

knots, if
where is a free threshold parameter (to be specified below).
The computation of the eigenvalues has been performed on a regular grid, corresponding to a cell size of . We use a triangularshaped cloud (TSC) prescription for the assignment of the particles (Hockney:1988) and then compute the overdensity and the eigenvalues of the velocity shear tensor for every grid cell. Using the AHF catalogues, we assign every halo to the nearest grid point hence providing us with a measure of environment for every object.
At this stage we still have not explicitly classified the cosmic web, as we lack a clear theoretical prescription for the value of . In our case, we have fixed to the highest value which ensures that no halo with belongs to a void in any simulation. At a first glance, this kind of constraint might seem redundant, as it would be implied in any standard definition of void as an underdense region. However, we must recall here that our definition of the cosmic web relies solely on the dynamical properties of the matter distribution (being related to the magnitude of its inflow or outflow in a given node) and may in principle overlook its net density content. It is thus necessary to enforce this principle explicitly tuning our free parameter to , which is the value which in this case satisfies the aforementioned condition and has been used in Section 5. For a more elaborate discussion of we refer the reader to Hoffman:2012; we only note that our choice is close to their proposed value.
4 Large scale clustering and general properties
Before presenting the results relative to the properties of the cosmic web and the correlation of halo properties to the environment, we will describe some aspects of large scale structure (LSS) and general halo properties in our simulations. This should give a more traditional overview of the effects of (coupled) dark energy models.
4.1 Halo mass function
The halo mass function (HMF) in coupled dark energy cosmologies has already been studied by Maccio:2004; Nusser:2004; Baldi:2010a; Li:2011; Cui:2012 so that we will only briefly comment on the topic. Our results reproduce the earlier findings of Baldi’s in the overlapping regions of mass and space, thus providing an additional proof of the correct functioning of our modified implementation.
In Fig. 1 we show the cumulative mass (left) and velocity functions (right) as well as the ratio to CDM for the four quintessence models. Singling out the region from to , and neglecting the higher mass end, where the statistics is unreliable due to the low number of objects, we can see that the largest difference in number counts amounts to for the strongest coupled models, gradually decreasing for smaller couplings. In the velocity function, this suppression reaches , thus slightly enhancing the magnitude of the effect.
We can compare our results for the HMF with those of Cui:2012, who modeled the Jenkins:2001 and Tinker:2008 mass functions for a series of similar coupled dark energy models, using FriendsofFriends (FoF) and Spherical Overdensity (SO) algorithms to build up their halo catalogues. Even though in their simulations they used different normalizations, fitting the analytical HMFs to the numerical results they were able to extend the predictions for cDE cosmology to arbitrary values. Using the same CDM normalization, then, they also found a suppression of the HMF of cDE, in perfect agreement with our results.
Although not shown here, we have also verified that these results match the analytical prediction of the Tinker mass function (Tinker:2008), provided the correct input power spectra and normalizations are used. We can safely conclude that the presence of coupled and uncoupled quintessence of the kind described here is expected to produce differences from CDM predictions up to a factor in present day’s HMF. Remarkably, this estimate is qualitatively independent of the algorithm used for the halo identification as we have seen comparing our results to the work of Cui:2012.
Mass cut  CDM  uDE  cDE033  cDE066  cDE099 

4.2 Halo properties
To study internal halo properties (such as spin parameter and concentration) we first need to define a statistically sound sample of objects, in order to reduce the impact of spurious effects on the results. This means that we need to constrain our analysis to structures which satisfy some conditions on both resolution and relaxation.
The first condition means that we have to restrict our analysis to objects with a number of particles above a given threshold, taking into account the existing tradeoff between the quality and the size of the halo sample. The second criterion needs to be applied as we want to focus on structures as close as possible to a state of dynamical equilibrium. In fact, many phenomena, such as infalling matter and major mergers, may take place, driving the structure out of equilibrium. In this case, then, the determination of quantities such as density profiles and concentrations becomes unreliable (see for instance Maccio:2007 and MunozCuartas:2011).
Following Prada:2012, we will define as relaxed only the haloes that obey to the condition
(8) 
without introducing other selection parameters; for alternative ways of identifying unrelaxed objects we refer for instance to Maccio:2007; Bett:2007; Neto:2007; Knebe:2008; Prada:2012; MunozCuartas:2011; Power:2011. For the moment, we neglect the impact of uDE and cDE on the definition of the virial ratio since this effect is of just a few percent (Abdalla:2010; Pace:2010) and is thus subleading in our case, where we are removing objects off by more than from the standard relation.
Now that we have established the rules that will shape our halo sample, we proceed to study some internal properties of dark matter haloes, namely, spin and concentration – as a function of halo mass – enforcing one additional criterion for the halo selection: the number of particles in it. When studying the spin parameter, we will restrict ourselves to haloes with , i.e. composed of at least baryon and dark matter particles, following the choices of Bett:2007, Maccio:2007, MunozCuartas:2011 and Prada:2012. In the case of halo concentrations, we applied a stricter criterion, using (or particles), due to the fact that the computation of halo concentration requires a better resolution of the central regions, as we will discuss in the dedicated subsection.
4.2.1 Spin parameter
We can study the rotational properties of haloes introducing the socalled spin parameter (e.g. Barnes:1987; Warren:1992), a dimensionless number that measures the degree of rotational support of the halo. Following Bullock:2001, we define it as
(9) 
where the quantities the total angular momentum , the total mass , the circular velocity , and the radius are all taken as defined by Eq. (4), with ; in cosmological simulations, the distribution of this parameter is found to be described as lognormal (e.g. Barnes:1987; Warren:1992; Cole:1996; Gardner:2001; Bullock:2001; Maccio:2007; Maccio:2008; MunozCuartas:2011)
(10) 
even though some authors (e.g. Bett:2007) claim that this should be slightly modified.
Due to the nonGaussian nature of this distribution, instead of the average value we plot in Fig. 2 the median value of the spin parameter as a function of halo mass. A weak negative correlation of spin to the halo mass can be observed here for haloes above , (as noted for instance by Maccio:2007 and KnebePower:2008) while the relation is positive below that threshold. However, cDE models have on average a higher value (per mass bin) compared to uDE and CDM. Albeit small, this increase in is clearly a coupling related effect, the magnitude of which is directly proportional to the value of . Given the small error bars (due to the large number of objects used in this analysis) we are confident that this is a real effect. Moreover, a similar result has been found by Hellwing:2011, that also claimed to have observed a link between fifth force and larger s.
A deeper investigation of the physical link between the coupling and increased rotational support is left to an upcoming work (Carlesi et al., in prep.) where the evolution of different parameters under different cosmologies will be analyzed. For the moment it is important to note that there appears to be some evidence of a link between the coupling strength of the fifth force and the corresponding degree of rotational support in dark matter haloes.
4.2.2 Concentration
Model  CDM  uDE  cDE033  cDE066  cDE099 

Dark matter density profiles can be described by a NavarroFrenkWhite (NFW) profile (NFW:1996), of the form
(11) 
where the , the so called scale radius, and the density are in principle two free parameters that depend on the particular halo structure. Using Eq. (11) we can define
(12) 
which is the concentration of the halo, relating the radius to the scale radius . Fitting Eq. (11) to our halo sample we observe that no substantial difference can be seen in the different simulations, that is, the NFW formula describes (on average) equally well dark matter halo profiles in CDM as in the other (coupled) dark energy models. While this is in contrast with the early findings of Maccio:2004, it is however in good agreement with the subsequent works of Baldi:2010a and Li:2011, who also found the NFW profile to be a valid description of DM haloes in interacting cosmologies. Thus, defining concentrations using Eq. (12) will not pose any problems nor introduce any systematic effect due to the fact that the NFW profile might only be valid for CDM dark matter haloes.
In Fig. 3 we now show the median concentration for objects in a certain mass bin: cDE cosmologies have a smaller concentration than CDM i.e. the larger the the smaller the ; whereas the opposite is true for the uDE model. This can partly be explained by the fact that concentrations are related to the formation time of the halo, since structures that collapsed earlier tend to have a more compact centre due to the fact that it has more time to accrete matter from the outer parts. Dynamical dark energy cosmologies generically imply larger values as a consequence of earlier structure formation, as found in works like those by Dolag:2004, Bartelmann:2006 and Grossi:2009. In fact, since the presence of early dark energy usually suppresses structure growth, in order to reproduce current observations we need to trigger an earlier start of the formation process, which on average yields a higher value for the halo concentrations. However, as explained in Baldi:2010a, smaller concentrations in cDE models are not related to the formation time of dark matter halos, but to the fact that one of the effects of coupled quintessence is to effectively act as a positive friction term. This means that dark matter particles have an increased kinetic energy, which moves the system out of virial equilibrium and causes a slight expansion, resulting in a lowering of the concentration.
In the hierarchical picture of structure formation, concentrations are usually inversely correlated to the halo mass as more massive objects form later; body simulations (Dolag:2004; MunozCuartas:2011; Prada:2012) and observations (Comerford:2007; Okabe:2010; Sereno:2011) have in fact shown that the relation between the two quantities can be written as a power law of the form
(13) 
where and can have explicit parametrizations as functions of redshift and cosmology (see Neto:2007; MunozCuartas:2011; Prada:2012). When we fit our halo sample to this relation using and as free parameters we obtain the bestfit values as shown in Table 5. Our values are qualitatively in good agreement with the ones found by, for instance, Maccio:2008 and MunozCuartas:2011 for CDM; but we do find some tension with the findings of Prada:2012. However, since they use a different algorithm for the determination of (which, according to them, leads to higher concentration values) and a different normalizations we cannot directly compare our results to theirs. On the other hand, uDE values are generally in agreement with Dolag:2004, DeBoni:2013 although in both cases there are again some discrepancies in the bestfit result, most probably due to the much different used in their simulations. For cDE we cannot directly compare our concentrationmass relation to the one obtained by Baldi:2010a since they do not provide any fit to Eq. (13).
5 Properties of the cosmic web
cell type  CDM  cDE  cDE033  cDE066  cDE099 

void  0.103  0.103  0.102  0.103  0.102 
sheet  0.343  0.343  0.344  0.344  0.343 
filament  0.437  0.438  0.443  0.442  0.443 
knot  0.116  0.115  0.109  0.111  0.118 
cell type  CDM  cDE  cDE033  cDE066  cDE099 

void  0.103  0.103  0.103  0.104  0.102 
sheet  0.349  0.348  0.348  0.347  0.346 
filament  0.449  0.450  0.450  0.449  0.453 
knot  0.097  0.097  0.097  0.098  0.098 
cell type  CDM  cDE  cDE033  cDE066  cDE099 

void  0.337  0.338  0.338  0.337  0.334 
sheet  0.460  0.456  0.460  0.460  0.461 
filament  0.185  0.184  0.185  0.185  0.186 
knot  0.017  0.017  0.017  0.017  0.018 
We now turn to the study of the cosmic web, as defined in Section 3.3, in CDM, uDE, cDE033, cDE066 and cDE099. In Fig. 4 and Fig. 5 we give a visual impression of the web classification (left) and the underlying dark matter density field (right) for a slice of thickness one cell (i.e. ) using a logarithmic colouring scheme for the density. From Figs. 4 & 5 it is evident that there is, in general, a very close correspondence between and filamentary and knotlike regions; just like between and void and sheetlike ones, so that the kinetic classification does provide in general a faithful description of the underlying density distribution – as shown in Hoffman:2012. Nonetheless, a minor number of cells do indeed violate this principle. In fact, as also noted by Hoffman:2012, in a very limited number of cases it happens that, for cells placed in the interior of a of a large dark matter halo, the velocity field will be determined by the motion of its virialized particles and not reflecting the cosmic web, respectively. On top of that, we must not forget that the freedom in the choice of the threshold , and the fixed spacing of the grid account for the fact that on scales smaller than we cannot properly resolve the complex shape of the web, which would probably require a more flexible grid implementation (Platen:2011). However, all these shortcomings do not seriously invalidate this description, as the number of such cells is generally small (for example, points defined as voids with sum up to less than of the total in all simulation, and independent of the simulation). In fact, the latter is the most important condition that we need to ensure, so that the existence of small biases disappears when considering ratios to CDM, which is at the core of the analysis we are carrying.
In Tabs. 8, 8 & 8 we show the mass and volume filling fractions as a function of cell type and cosmologies. These values are estimated simply summing all masses and volumes contained in cells belonging to the same kind of environment. What is clear by looking at these results is that the general structure of the cosmic web is almost left unchanged across models. In fact, discrepancies among different cosmologies are much less than in this regard, thus making it hard to detect deviations from CDM by simply considering the volume and the mass associated to the various kinds of environment. The same conclusion can be drawn if we look at Fig. 6, which shows the distributions of the three eigenvalues of , that appear to be identical and thus provide no leverage to distinguish the models under investigation here.
The gas distribution through the different node types seems also to be largely unaffected by the different cosmology: As we can see from Table 8, the mass fractions of gas are substantially identical throughout all the models, without any significant discrepancy. Comparing to the distribution of dark matter, we do notice a slight increase in the fraction of gas belonging to sheets and filaments paralleled by its reduction on knots, a pattern which is observable in all the models to the same extent.
We remark that our results for uDE agree with Bos:2012, who also found that quintessence cosmologies with RatraPeebles potentials do not lead to significant changes in the general properties of the cosmic web. We also emphasize that our findings relative to void regions are largely independent of the choice of . Using different threshold values we have been able to test this and see that void distributions are affected to the same degree in all the different models, confirming this particular result does not depend on our .
6 Halo properties in different environments
We now turn to the study of halo properties classified according to their environment; this kind of analysis has already been done for CDM using both geometrical (e.g. AvilaReese:2005; Maccio:2007) and dynamical (e.g. Hahn:2007; Libeskind:2012; Libeskind:2013) web classifications, finding in general a correlation between halo properties such as spin and shape to its surrounding environment.
Using the information from the halo catalogues we proceed to assign each halo to the nearest grid point and build up four different halo samples, one for each cell type. Then we repeat the analysis presented in Section 4 for the halo counts (velocity and mass), spin and concentrations. We will see that this kind of separation of haloes enhances some of the differences already seen in general among different cosmological models and is therefore of great importance when trying to constrain more effectively coupled and uncoupled scalar field cosmologies.
Voids and sheets are readily identified with underdense regions, as has also been confirmed by the analysis presented in the previous Section 5. And the fact that for these cells at most one eigenvalue of the shear tensor has a value above means that in two or more spatial directions there is a net outflow of matter, which is in turn associated with a matter density below the average. For overdense regions (i.e. filaments and knots) there is a net inflow of matter towards the center of the cell from at least two directions. Following this we partition the subsequent study into underdense regions on the one hand (using voids and sheets) and overdense regions (i.e. filaments and knots) on the other.
Underdense regions in CDM are usually associated with lower spins and slightly larger halo concentrations (Maccio:2007), raising the question whether this still holds for (coupled) dark energy cosmologies.
Environment  CDM  uDE  cDE033  cDE066  cDE099 

void  
sheet  
filament  
knot 
6.1 Halo number counts
Even though, by definition, underdense regions are less populated, nonnegligible fractions of the total halo count can be still found in voids and sheets, as shown in Table 9, ensuring that the samples used are reasonably large, and allow us to draw credible conclusions.
6.1.1 Underdense regions: voids and sheets
We notice that in underdense regions (i.e. voids and sheets, shown in the upper two plots) the trend persists that the number of objects is smaller for cDE models than for CDM, something also observed for the general halo sample. However, it is important to remark that singling out and counting the objects belonging to the underdense parts only, we end up observing larger differences among the models. This effect also appears to be much stronger in cDE than in uDE. In fact, whereas the differences in number counts of objects does not exceed , when restricting halo counts to void regions only, we can see that cDE models’ underprediction is much larger and peaks at (ignoring the higher mass ends, where only a small number of objects is found). It is also clear from Fig. 8 that while the sign of the effect is very similar in both voids and sheets, its strength is slightly reduced in the latter type of web, suggesting that there exists at least a mild dependence of this phenomenon on the specific kind of environment. Although we are not showing it here, we have also carefully checked that this result is substantially independent from the kind of chosen. In fact, repeating our computation using higher threshold values, we see that the magnitude of the effect does not change substantially. The physical mechanism behind this effect is understood and provides a consistent framework for interpreting our results. In fact, as first explained by Keselman:2010 and subsequently confirmed by Li:2011, fifth forces enhance the gravitational pull towards the overdense regions, quickly evacuating matter from underdense regions. This causes these environments to have less structures, so that in the end the number of haloes left in voids will be comparatively smaller than in the noninteracting cases, as found in our simulations.
However, we need to make an additional remark on this result before proceeding to the next section. In fact, we note that our choice of the normalization, which is taken to be the same at , plays an important role in the result just described. It is in fact known (see e.g. Baldi:2010) that using a different normalization prescription (for instance, at the redshift of the CMB for the matter density fluctuations), coupled models end up predicting (in total) more objects than CDM. Hence, at this stage we cannot completely disentangle the influence of our choice of the normalization of the initial conditions from the genuine influence of the additional interaction.
6.1.2 Overdense regions: filaments and knots
Like in the case of underdensities, we notice for overdense regions (shown in the lower two plots) that the trend of suppression which characterizes the general halo counting still holds, even though now the strength of this effect is slightly smaller across all cosmologies. This is not unexpected, since the effect seen in the HMF discussed in Section 4 has to be obtained from a combination of both underdense and overdense structures, and should therefore result in an intermediate value for cDE halo underproduction. Again, we have checked that the chosen threshold for the eigenvalues of the velocity shear tensor does not substantially affect this conclusion.
We can therefore state that there is a progression towards smoothing out the differences among different cosmologies while moving to increasingly higher density regions. This is a very important result that indicates that underdense regions should be the target of choice when searching for the effects of additional long range gravitationallike forces.
This result is in line with what has been already found for other fifthforce cosmologies (see Martino:2009; Keselman:2010; Li:2012; Winther:2012), where the environmental dependence and in particular the properties of voids were stressed as powerful tests for additional interactions and modifications of standard Newtonian gravity. It is in fact well known that void properties are extremely sensitive to cosmology (Lee:2007; Lavaux:2009; Bos:2012; Sutter:2012) and hence provide a powerful probe of alternative models. In particular, when the extra coupling in the dark sector is weak (as in the cases analyzed here) the complex evolution and phenomena that characterize the overdense regions may conceal its imprints, while void regions, whose dynamics is comparatively simpler, are expected to be more directly linked to the underlying cosmology.
Parameter  CDM  cDE  cDE033  cDE066  cDE099 

(v)  
(v)  
(s)  
(s)  
(f)  
(f)  
(k)  
(k) 
6.2 Spin and concentration
We now turn to the nondimensional spin parameter, , and dark matter halo concentrations, investigating how they will change across different environments and cosmologies. In the latter case, we will also pay particular attention to the environmentrelated changes to the relation of Eq. (13). In both cases we refer to the definitions introduced in Section 4.
6.2.1 Underdense regions: voids and sheets
Looking at the median spin parameters shown in the upper panel of Fig. 9 we can again draw the conclusion that, just like in the general case, cDE cosmologies lead to larger spins and that this increase is proportional to the coupling parameter . On the other hand, the value of for haloes in uDE cosmologies is, on average, indistinguishable from CDM. We can therefore confirm the observation that underdense region contain haloes with lower spins, just as found by Maccio:2007. However, the reduction in the median value is of the same order in all models, so that combining the information of the environment does not put tight constraints on the parameters of the model.
Concentrations, too, show a remarkable behaviour for haloes belonging to underdense regions. In Table 10 we show the results of fitting the median concentration per mass bin to a power law, i.e. Eq. (13). The first thing we observe is that the correlation between and (as measured by the powerlaw index ) is weaker than what we observed in the general case. This, combined with the fact that in the lower mass bins median concentration do not change with respect to the general case, in turn leads to observed larger values for , although the errors are also large due to the small number statistics. However, some care must be taken when considering this relation for void haloes since the fit is based upon a small mass range only and also gives more weight to lower mass objects (Prada:2012; DeBoni:2013).
6.2.2 Overdense regions: filaments and knots
In the bottom panels of Fig. 9 and Fig. 10 we plot the spin and concentrationmass relation; the bestfit values to Eq. (13) are again provided in Table 10.
In the case of spins, we find that dark matter haloes in coupled cosmologies tend to be characterized by larger values of . However, haloes located in filamentary structures show, at least in the lower mass bins, sharper differences between cDE models and CDM than what is revealed by knots. This is also due to the smaller number of low mass haloes living in knots, which visibly affects the statistics of the parameter.
Concentrations instead show two slightly different patterns in filaments and in knots. In the former environment, all cosmologies seem to be characterized by a flatter slope, which averages around and seems not to be connected to the underlying model. In the latter environment, a steeper correlation is found, with – much closer to the general case discussed in Section 4. Not only the slope but also the normalization of Eq. (13) changes when considering filamentary or knotlike environments: in the former case we find that this parameter is substantially larger than in the latter.
Our results therefore indicate that the concentrationmass relation is not only affected by the cosmological model but also by the environment the haloes under consideration live in: gets flatter while increases for decreasing densities. However, at odds with what we found for halo number counts, we find here that environment does not play a role in strengthening the magnitude of modeldependent properties of haloes. While the effect of dark energy can still be clearly seen in the higher spins and lower concentrations of dark matter haloes, these cDEinduced characteristics are not enhanced by the environment. In fact, whereas the halo content of the different regions depends on the model and reinforces the trends observed in Section 4.2, the properties of the haloes themselves, while still being correlated to the underlying cosmology, seem to be shifted by the same amount as CDM. This suggests that environmental effects, in these cases, influence to the same extent both quintessential and standard models, and do not provide a stronger modelspecific kind of prediction.
7 Conclusions
In the present work – which forms part of a series of studies of (coupled) dark energy models – we have discussed the properties of largescale structures and the cosmic web as they emerge in a series of different quintessence models, systematically comparing the results of a coupled scalar field to those obtained for a free field and the standard CDM cosmology.
We performed the following threefold analysis:

we studied halo mass function and general halo properties (mass, spin and concentrations),

we investigated the general properties of the cosmic web, using a kinetic classification algorithm,

we correlated halo properties to the environment.
First, we have studied several aspects of cDE and uDE cosmologies looking at the full halo sample. At this stage, our results proved to be in line with those of Baldi:2010a; Li:2011; Cui:2012, finding that the analytical formulae for the halo mass functions and dark matter profiles are valid also in this class of models.
Examining concentrations we found that, while uDE cosmology is characterized by haloes with higher values for , for cDE models the opposite is in general the case – in accordance with the results of Baldi:2010a and DeBoni:2013. Interestingly, in the case of spin parameters we observe a weak dependence on the coupling, since we can see that their value is mildly enhanced by larger values of , as was also noted by Hellwing:2011 in the context of other fifth force cosmological models.
The cosmic web investigated as part of this study is characterized by the eigenvalues of the velocity shear tensor, a novel method recently proposed by Hoffman:2012 and successfully applied to various simulations by Libeskind:2012. Computing the fraction of total mass and volume belonging to each type of environment in our cosmologies, we find that the structure of the cosmic web itself does not reveal any particular difference among the models. The same conclusion can be drawn when investigating the global distribution of the shear tensor eigenvalues.
This notwithstanding, the classification of the cosmic web can be extremely useful when married with the halo catalogue. Combining the two, in fact, we were able to show that many of the differences observed in some halo properties when studying a global sample of relaxed structures above a threshold mass are in fact due to objects belonging to a certain type of environment. This happens in particular in voids and sheets, where the differences among cDE and CDM are up to three times as large as they are in the general case. We have been able to verify how the magnitude of this effect is closely dependent on the coupling: while cDE cosmologies’ underproduction of haloes in these regions is largely amplified, the overproduction that characterizes the uDE model investigated here is only weakly enhanced. This means that;

one should focus on voids and sheets (underdense regions) when looking for signatures of (coupled) dark energy, and

the magnitude of the deviations from CDM allows us to place constraints on cDE cosmologies, or even detect them.
We have also seen how the standard concentrationmass relation is substantially affected when fitted for halo samples classified according to the environment they are located in. While the standard functional form of Eq. (13) still holds in underdense regions for all the models, it does so with a much steeper slope (the change is from an average of to ) and a substantial increase in the average concentration. In addition to this, we note again an amplification of the difference between the values in cDE and CDM obtained when fitting Eq. (13) in voids and sheets (up to ) with respect to the global one ( per cent).
The fact that these results are mostly visible when restricting the halo sample to underdense environments tells us the importance of the relative weights to be attached when performing global analyses. Indeed, when referring to halo properties in general, we do in fact hide a large number of peculiar features which can be seen only in a narrower subset. In the particular case of cDE we have seen how structures in voids play a major role, in the same direction of Li:2011, who also highlighted the importance of underdense regions in the context of similar cosmological models.
As a concluding remark, we would like to emphasize that a great amount of effects observed here still deserve a more indepth study. In particular, the analysis of the temporal evolution of halo parameters, will shed more light on the mechanisms that result in the previously discussed differences at and increase the observational features that can be used to constrain quintessence models. We shall turn to these in future contributions in this series.
Acknowledgements
EC is supported by the Spanish Ministerio de Economía y Competitividad (MINECO) under grant no. AYA201231101, and MultiDark Consolider project under grant CSD200900064. He further thanks Georg Robbers for providing an updated, nonpublic version of CMBEASY.
This work was undertaken as part of the Survey Simulation Pipeline (SSimPL: ssimpluniverse.tk) and GFL acknowledges support from ARC/DP 130100117
AK is supported by the Spanish Ministerio de Ciencia e Innovación (MICINN) in Spain through the Ramón y Cajal programme as well as the grants AYA 200913875C0302, CSD200900064, CAM S2009/ESP1496 (from the ASTROMADRID network) and the Ministerio de Economía y Competitividad (MINECO) through grant AYA201231101. He further thanks Emily for reflect on rye.
GY acknowledges support from MINECO under research grants AYA201231101, FPA201234694, Consolider Ingenio SyeC CSD20070050 and from Comunidad de Madrid under ASTROMADRID project (S2009/ESP1496).
The authors thankfully acknowledge the computer resources, technical expertise and assistance provided by the Red Española de Supercomputación.
We further acknowledge partial support from the European Union FP7 ITN INVISIBLES (Marie Curie Actions, PITNGA2011289442).
All the simulations used in this work were performed in the Marenostrum supercomputer at Barcelona Supercomputing Center (BSC).