# Using large scale structure data and a halo model to constrain Generalised Dark Matter

###### Abstract

Constraints on the properties of the cosmological dark matter have previously been obtained in a model-independent fashion using the Generalised Dark Matter (GDM) framework. Here we extend that work in several directions: We consider the inclusion of WiggleZ matter power spectrum data, and show that this improves the constraints on the two perturbative GDM parameters, and , by a factor of 3, for a conservative choice of wavenumber range. A less conservative choice can yield an improvement of up to an order of magnitude compared to previous constraints. In order to examine the robustness of this result we develop a GDM halo model to explore how non-linear structure formation could proceed in this framework, since currently GDM has only been defined perturbatively and only linear theory has been used when generating constraints. We then examine how the halo model affects the constraints obtained from the matter power spectrum data. The less-conservative wavenumber range shows a significant difference between linear and non-linear modelling, with the latter favouring GDM parameters inconsistent with CDM, underlining the importance of careful non-linear modelling when using this data. We also use this halo model to establish the robustness of previously obtained constraints, particularly those that involve weak gravitational lensing of the cosmic microwave background. Additionally, we show how the inclusion of neutrino mass as a free parameter affects previous constraints on the GDM parameters.

###### Contents

- I Introduction
- II Brief review of GDM model and previous constraints
- III GDM halo model
- IV Description of datasets and methodology
- V Constraints
- VI Conclusion
- A CDM halo model
- B Use of the matter continuity equation for redshift space distortions and BAO reconstruction

## I Introduction

The CDM cosmological model does a good job of reproducing the current cosmological observations. In this model, the standard model of particle physics is supplemented by a cosmological constant and a dark matter particle. This dark matter particle is assumed to interact purely due to the influence of gravity and to have a negligible (initial) velocity dispersion, thus the name Cold Dark Matter (CDM). In perturbative calculations this is typically modelled as a pressure-less perfect fluid. As a result, many cosmological constraints on the dark matter density are, more correctly, constraints on the density of this pressure-less perfect fluid. More generally, CDM is evolved by solving the collision-less Boltzmann equation. This is done on large scales using cosmological perturbation theory (implemented in Boltzmann codes such as class and camb) and on smaller scales using N-body simulations and other non-linear methods.

Since we are entering the era of so-called “precision cosmology,” in which many cosmological parameters have been measured with 1% accuracy or better, it is timely to consider whether such an idealised and simple dark matter model is sufficient when analysing the data. There are many physical dark matter models that do not yield precisely CDM, for example Warm Dark Matter (WDM) Dodelson and Widrow (1994); Armendariz-Picon and Neelakanta (2014); Piattella et al. (2015) or ultra light axions (one example of Fuzzy Dark Matter (FDM)) Hu et al. (2000); Hlozek et al. (2015). In addition, recent work on the Effective Field Theory of Large Scale Structure (EFTofLSS) Baumann et al. (2012) shows that even an ideal CDM candidate develops a more complicated energy momentum tensor, even on linear scales, once the non-linearities that inevitably form on small scales back-react on the large scales. This causes an effective pressure and viscosity on large scales. From a non-cosmological perspective, despite a large number of direct and indirect detection experiments for dark matter, no convincing detections have been made, and many theoretically favoured regions of parameter space have been ruled out Aprile et al. (2012, 2014); Buckley et al. (2013); Olive et al. (2014); The CRESST Collaboration et al. (2015); LUX Collaboration et al. (2015). Thus, there are strong reasons to go beyond the simplest ways of modelling dark matter.

In Kopp et al. (2016), the Generalised Dark Matter (GDM) model (first proposed in Hu (1998)) was examined in some detail, notably how it relates to particular physical models. GDM adds to the CDM energy momentum tensor a background pressure, pressure perturbation and anisotropic stress. Closure relations are then postulated to match qualitative properties of known models, like massive neutrinos, and in order to de-correlate background and perturbative properties. GDM encompasses WDM, FDM and the EFTofLSS effects as well as other physical models, so it is sufficiently versatile for examining dark matter properties in a model independent fashion. In Thomas et al. (2016), all GDM parameters were constrained using Cosmic Microwave Background (CMB) data, supported by additional data on the cosmological expansion history (see section 4.3 in this paper and references therein for comparison to earlier works constraining partial or similar parameters to those we consider here, such as Müller (2005); Calabrese et al. (2009); Xu and Chang (2013)). The results showed no evidence for any non-CDM properties of dark matter. This was expanded on in Kopp et al. (2018a), where an improved freedom was given to one of the GDM parameters; this was used to demonstrate for the first time that there is no cosmological epoch where the data would favour a nonzero equation of state, and furthermore that there is no cosmological epoch where the data is consistent with zero dark matter density, thus showing the strength of the GDM approach to testing the CDM paradigm. An independent group subsequently Kunz et al. (2016) verified some of the results in Thomas et al. (2016), as well as using some late time matter clustering data; we will comment further on this later in the paper. Further work constraining the GDM parameters is now ongoing by other groups, see e.g. Tutusaus et al. (2016).

It was noted in Thomas et al. (2016) that matter power spectrum data could not only improve the constraints on the GDM parameters, but also has the potential to break a degeneracy between two of them (see section II). The robust use of such data requires a non-linear extension of the GDM model, which is not present in the literature. It was also noted in Thomas et al. (2016) that the inclusion of a non-linear extension to perturbation theory for CDM makes a difference to the CMB lensing potential. This effect is of a similar magnitude, but opposite sign, to that of GDM with parameters saturating the constraints found in Thomas et al. (2016). In this paper we develop a halo model for GDM which allows us first to test the robustness of the results in Thomas et al. (2016), and second to use matter power spectrum data from the WiggleZ survey Parkinson et al. (2012) to improve the constraints on the GDM parameters. The paper is laid out as follows: in section II we briefly review the GDM model and previous constraints, before constructing the GDM halo model in section III. We then present the methodology for our constraints in section IV and present our resulting constraints and robustness tests in section V. We conclude in section VI.

## Ii Brief review of GDM model and previous constraints

The GDM model was first proposed as an extension to the standard CDM model in Hu (1998). Here we give brief details of the model, following Kopp et al. (2016); see both this work and Hu (1998) for further details of the model, its motivation and the different physical models that it can encompass.

The standard CDM energy momentum tensor is given by

(1) |

i.e. the fluid is specified entirely by its density, , and velocity . This is then typically divided into a background part that is homogeneous and a perturbation. The GDM parameterisation adds pressure and anisotropic stress to this, giving

(2) |

The pressure and density perturbations are divided into background quantities (denoted by an overbar) and perturbed quantities as usual, and the additional scalar perturbations of and are controlled by the equation of state (background pressure), the sound speed (pressure perturbation) and the viscosity (anisotropic stress). The equation of state relates to the background quantities in the usual way: , and the additional perturbations are governed by the closure equations

(3) | |||||

(4) |

Here, and are (a particular choice of) gauge invariant density and velocity perturbations for GDM, and and are the pressure and (scalar) anisotropic stress perturbations. The adiabatic sound speed is . Note that overdots refer to conformal time . See Kopp et al. (2016) for an in depth explanation of our notation and of this choice of closure equations. Note that recovers the pressureless perfect fluid, and therefore the standard CDM cosmological model.

Both and cause a decay in the gravitational potential power spectrum on scales below in a GDM dominated universe Kopp et al. (2016). In addition, if is sufficiently larger than , then it causes oscillations in the density perturbation below the Jeans length. Although we refer to as the sound speed, this is technically only true if (see Kopp et al. (2016)). The viscosity damps the density perturbations without causing any oscillations. As expected, the equation of state changes the expansion history of the universe for a fixed . In particular, in Kopp et al. (2016) it was shown that the main effect is to change the time of matter-radiation equality, and thus to change the relative heights of the peaks in the CMB. In addition, changes the distance to the last scattering surface.

The aforementioned phenomenology of the parameters is all manifest in Thomas et al. (2016). In this work we took simple forms of the GDM parameters, giving them a single value with no time and scale dependence. We then constrained these parameters using Planck CMB data (temperature, polarisation and the lensing potential) Planck Collaboration et al. (2015a), BAO data Beutler et al. (2011); Anderson et al. (2014) and a prior from the HST key project Riess et al. (2011). We found upper bounds on and of and respectively (at the 99.7% confidence level), in line with the degeneracy expected if was primarily constrained by the CMB. We also put constraints on and found degeneracies between and due to the effects above. In all cases we found no evidence for a non-zero value of any of these parameters. The first three rows of table 1 summarise the constraints obtained in previous work. The main conclusions of Thomas et al. (2016) were independently verified in Kunz et al. (2016). Furthermore, the assumption of a single time independent value for was relaxed in Kopp et al. (2018b), which showed that a non-vanishing dark matter background density is required at every epoch, thus showing the power of the GDM formalism for constraining extensions to CDM.

In Thomas et al. (2016) it was noted that the phenomenology of the GDM model suggests that the use of late time clustering data, such as the matter power spectrum, could significantly improve the constraints on and . In principle, if the data is precise enough to determine oscillations in the decaying region, then such data could also break the degeneracy between and . However, the use of matter power spectrum data requires going beyond the scales where perturbation theory is valid. The GDM model as described in this section is only defined perturbatively and thus must be extended in order to be applicable on smaller scales. One of the main results of this paper is the development of a halo model extension to the GDM model (see III): As well as seeking to improve the constraints on GDM using matter power spectrum data, we also to seek to understand how safe this process is, i.e. how robust it is to non-linear modelling. In Kunz et al. (2016), the authors used late time matter clustering data (weak lensing data) which probes the same underlying potential power spectrum as the matter power spectrum does. However, they did not consider a non-linear modelling of the GDM in that work; instead they used halofit Smith et al. (2003) as a sanity check, whilst noting themselves that halofit has limitations when applied outside of a CDM context. Our goal is thus not just to improve the constraints on the GDM parameters using matter power spectrum data, but also to be able to quantify how much we can trust any such results.

A further goal is to investigate the robustness of the constraints obtained in Thomas et al. (2016). More precisely, it was noted that even in CDM, using halofit makes a small difference to the lensing potential and thus to the lensed temperature and polarisation s. We thus wish to determine whether inclusion of a non-linear prescription for GDM would strengthen or weaken the constraints previously obtained. Nonlinearities typically act to increase the matter power spectrum, which is the opposite effect to that caused by increasing GDM parameters, see figure 1. Hence we expect that the constraints on and could be weakened once nonlinearities are included.

Since we are focussing on the complexities introduced by the non-linearities and additional datasets, we will work with single, constant values of the GDM parameters as in Thomas et al. (2016). In particular, the assumption of time independence means that .

## Iii GDM halo model

As stated in section II, the GDM model is only defined for linear perturbations and the homogeneous background, and thus cannot be constrained by all of the currently available cosmological data. One framework that has been used to make predictions on non-linear scales is the halo model Cooray and Sheth (2002). The halo model is a semi-analytic method for computing the matter power spectrum on non-linear scales that works from the premise that the matter is organised into halos and that averaging over all of the halos gives the mean matter density in the universe , where is the halo mass function. The two point correlation of the matter field will thus depend on the halo density profile and also the halo mass function. For the former we will use the empirical Navarro-Frenk-White (NFW) profile and for the latter excursion set theory to predict the so-called multiplicity function and relate it to the halo mass function through .

For more details we refer to appendix A, where we present the CDM halo model Seljak (2000). This also serves to introduce our notation and perspective, as these can vary between presentations of the halo model. We also present our mass function in this appendix. There is another commonly used non-linear correction for CDM: halofit Smith et al. (2003), which is an extension of the halo model that is calibrated against N-body simulations. As was already noticed in Kunz et al. (2016), it does not make sense to use halofit for GDM, as GDM is not part of the cosmologies that have been used to calibrate it. Furthermore, the numerical implementation in Boltzmann codes (e.g. class) simply crashes for values of GDM parameters where the linear power spectrum falls off too quickly. See figures 1 and 2 for the differences between halofit and the halo model as implemented in this paper for a CDM cosmology; it is known that halofit and the halo model differ for CDM Smith et al. (2003). This difference is largest, up to 15%, in the interval , visible in figure 2 comparing the red and dotted orange lines. However, for , relevant for our applications the agreement between our halo model and halofit is better than 2% and sufficient.

Examining the standard CDM halo model, we can see that there are several obstacles to simply applying the framework “as is” to GDM. Firstly, due to the large drop in power on small scales when and are non-zero, the equation defining does not necessarily have a solution. In addition, a lot of the formulas that are used are calibrated against CDM N-body simulations. We thus implement the halo model for GDM using a similar approach to the Warm and Fuzzy code Marsh (2016), which is designed for warm dark matter and axions (note that both of these can fit into the GDM parameterisation Kopp et al. (2016)). In this approach, four modifications are made relative to a CDM cosmology: A) a modified linear spectrum, B) a modified halo concentration, C) a mass dependent barrier related to critical density for spherical collapse, and D) a modified mass function to account for the non-Markovian corrections to standard Press-Schechter expression. We implement these corrections differently to the Warm and Fuzzy code; we detail our implementation here. Readers not interested in the construction and definition of our halo model can skip to section IV.

### iii.1 Modified linear spectrum

The simplest and most obvious modification is that the appropriate -GDM linear theory power spectrum is used as an input for the halo model, rather than a CDM spectrum. While this is done via fitting functions for the transfer functions in Warm and Fuzzy, we instead directly use the output from the full Boltzmann calculation from the modified class code.

### iii.2 Modified concentration

Following Marsh (2016), which itself follows Schneider et al. (2012) for WDM, we calculate the CDM value of the concentration. We then apply a correction according to

(5) |

where , and is the half-mode mass defined by

(6) | |||||

(7) |

Note that this functional form and the specific values were found to give good fits to FDM and WDM simulations. As GDM contains these two models as limiting cases we adopt this “as is” to GDM. The results are not sensitive to this guess. Furthermore, in the absence of non-perturbatively defined GDM model and cosmological GDM simulations, this is the best and most conservative choice we can make, as it is known to work in two limiting cases of the GDM model.

### iii.3 Mass dependent spherical collapse density and how it relates to the mass function

The central object in the excursion set theory is the so-called multiplicity function , determined by the first upcrossing rate of the smoothed random field through a barrier as a function of the smoothing scale . These crossings are then identified with proto-halos of size , and thus fixed mass , so that determines the mass function . In its most rudimentary form that mass function is mostly sensitive to the ratio

(8) |

where is the standard deviation, Eq. (32), of the (non-relativistic and participating in structure formation) matter perturbations smoothed on a scale , and is the spherical collapse barrier. The redshift is a time where linear perturbation theory applies to all scales of interest and thus the statistics of the density field is gaussian, but well after radiation domination (i.e. such that a spherical collapse threshold can be obtained neglecting the radiation component).

It is important to note that the -dependence in is determined by non-perturbative fluid dynamics, it is not a density field, or some linear extrapolation of it. Rather, assigns a collapse redshift to each region (assumed for simplicity to be spherically symmetric tophat profiles) of the gaussian density field filtered at scale at time ; this collapse redshift is then identified with the formation time of the halo of mass . This assignment involves some linear dynamics (at times not much later than ) when the field is still linear, but more importantly it includes the fully non-perturbative collapse that defines the collapse redshift at the time when formally the density contrast . This formal divergence is then associated with the time of halo formation, and thus approximated as instantaneous event.

#### iii.3.1 Standard CDM linear extrapolation

Typically, the mass function is not written in terms of equation (8) but instead in terms of

(9) |

where

(10) |

This “linear extrapolation” from to is possible, because in a purely CDM and dominated universe the growth is scale independent. Thus,

(11) |

Writing things in this way in CDM is convenient because it happens (see Weinberg (2008), Chapter 8.2) that the linearly extrapolated , , is a constant in an EdS universe, and is only very mildly dependent on (but still independent of ) in a CDM dominated universe. More precisely

(12) |

This is the reason why it is common in the literature to use this linearly extrapolated (denoted here by ) and the variance of a fictitious linearly extrapolated density field, even though there is no physical interpretation for such an extrapolation. In order to more clearly separate linear from non-linear physics, and to make the least amount of guessing to arrive at our GDM mass function, we avoid here using the linearly extrapolated , , and let the excursion set theory unfold at .

In order to proceed this way we rewrite the multiplicity function , defined in Eq. (31), appearing in the mass function (see appendix A) using (11). The relevant term appearing in the mass function is , see equation (36). Inspired by equations (8) and (11) we multiply numerator and denominator by , and define

(13) |

where , the barrier is now written in terms of the “strongly time-dependent” spherical collapse barrier and , and in , is replaced with , making it manifest that the excursion set theory is applied to the random field at . The Markovian part (see later) of the multiplicity function of the mass function is thus given by (35), replacing with and with ,

(14) |

At this point, this mass function is generic and is not derived in the context of any particular extension to CDM cosmology.

To consider how this mass function relates to the standard CDM case, note that if the term is ignored then this mass function reduces exactly to equation (35), just without the assumption of scale independent growth. This re-writing thus makes it clearer how scale dependent growth should manifest in the formalism. For the case of scale independent growth, the new derivative term in (14) reduces to the previous form

(15) |

thus the mass function reduces exactly to the standard form for the CDM case.

#### iii.3.2 GDM approach

We will now postulate the spherical collapse barrier for GDM by reversing the direction of definition in equation (11). While in CDM, is defined via (11), leading to (10), we now assume for GDM the validity of equation (11)

(16) |

but use it to define while fixing on the right hand side to be given by (12). In the following we will drop the superscript GDM and simply write for the GDM spherical collapse barrier. The idea behind equation (16) is that this definition implements the intuitive idea that if in GDM power is removed in a scale dependent way at then the collapse should be inhibited compared to CDM, and thus the spherical collapse barrier should be increased. It is at this point unclear how to judge whether the threshold defined in this way is correct given the absence of spherical collapse simulations within a (still to be) non-linearly defined GDM model. However, in addition to increasing the threshold whenever power is removed in GDM compared the CDM, this definition smoothly and naturally reduces to the CDM prescription in the CDM limit of GDM, and can thus be considered conservative for small GDM parameters. The mean barrier for collapse at in CDM can be approximated by (and we postulate this to hold in GDM too)

(17) |

where is a purely time dependent function. The multiplicity function (14) takes into account that the mean barrier for a randomly selected point deviates from the spherical collapse barrier because collapse is not spherically symmetric, and also that the barrier is diffusive, rather than 100% absorbing. The former is parameterized by , the latter by . We saw above that a consistent choice in the context of scale dependent growth is . A conservative choice for the GDM mass function thus is

(18) | ||||

(19) |

where is the largest smoothing scale used in the halo model code, and for all reasonably small values of GDM parameters, reduces to the growth function . We remove the scale dependence from for two reasons. Firstly, it is not clear what it would mean to include scale dependence here; it is not done in any version of the excursion set that we know of. In addition, we wish our halo model to act like a standard CDM halo model in the limit that the GDM parameters are zero. As we are working directly from a Boltzmann code and not performing the standard linear extrapolation to redshift zero, the radiation component of the universe will cause a scale dependent growth even in CDM. Thus we explicitly remove the scale dependence from this term and use a growth factor that is defined to be scale-independent even in GDM.

In order to calculate the mass function for GDM (or any other cosmology with scale-dependent growth), we need to evaluate the term. The result is

(20) |

Thus, the moving barrier term is given by

(21) |

which reduces to in the case where the initial and final have the same shape. That is often taken to be exactly correct in a CDM universe. More correctly, it depends on the value of : for example, if since there is still enough radiation “contamination” left to modify the shape of the matter transfer function, such that the ratio is not independent. The modification of the mass function is however less then for CDM and can be neglected, see figure 10. However for GDM, where in general growth is scale dependent even during matter domination, we do not expect the second term in the bracket to be close to 1.

We will assume that GDM does not change the values of and , which amounts to the assumption that collapse inhibition from asphericity and the scatter of the barrier due to environmental and stochastic processes are unchanged. In principle, these quantities could be measured in FDM and WDM simulations. We are not aware of any such measurements and so using the CDM values seems to be the most sensible approach.

### iii.4 Non-Markovian corrections and asphericity of collapse

So far we only looked at the mass and time dependent spherical part of the collapse barrier. Now we turn to the mass and time dependent contributions to the barrier due to the asphericity of the collapse as well as the non-Markovian corrections to the mass function. For GDM, the Markovian part of the mass function, (35) is replaced by (14).

The non-Markovian corrections will be implemented in a similar fashion as done in Kopp et al. (2013). There a mass dependent spherical collapse barrier was obtained due to a modification of gravity that left the background cosmology unchanged. It was shown that the non-Markovian corrections could be included through a simple relation , where includes the known and calculable non-Markovian corrections in CDM. In our case we cannot use as reference CDM because the background might be different in GDM due to . For that reason we will define another non-Markovian reference mass function , using a mass-independent spherical collapse barrier closer to . We choose that reference mass-independent spherical collapse barrier to be , i.e. the GDM spherical collapse barrier evaluated at the largest mass scale used in the halo model code. The reference mass function thus corresponds to a fictitious GDM model in which the spherical collapse barrier is mass independent. The reason for choosing is that this barrier will be similar to that of a GDM model with (since in the limit the effects of and disappear). This way we don’t need to run class twice for each model. The reference multiplicity function is given by

(22) |

where

We will also include a further correction to the mass function that has been observed to fit mass functions measured in warm dark matter simulations Schneider et al. (2012); Marsh (2016). The origin of that correction is likely to also be non-Markovian in nature, and it arises due to the absence of power below the scale . If the density field does not perform a random walk as function of it can happen that the mass function suffers a cutoff, see Paranjape et al. (2012). If the power is sharply dropping for scales , then the density field at a fixed point no longer performs a random walk for varying for . Thus we expect the mass function to be more non-Markovian for those small scales, implying a cutoff determined by mass scale related to . We follow the fit of Schneider et al. (2012), which works well for WDM.

The final expression for the multiplicity function entering (31) then is

(23) |

where , the half mode mass (6), is used instead of , see Schneider et al. (2012). To summarize: the first two factors take into account non-Markovian effects. The first one that the random walk below is highly non-Markovian, the second one the standard non-Markovian corrections for a diffusive barrier of the form const+const. The last factor is the Markovian mass function for a moving diffusive barrier (17). As part of the halo model, this provides a non-linear prescription for computing the GDM matter power spectrum and comprises one of the main results of this paper. When the GDM parameters are zero, the halo model reduces to a CDM halo model as expected.

We make one final modification to the CDM prescription in appendix A: In order to apply the compensation for the 1-halo term (see appendix A) in GDM, we need to define a scale independent growth for GDM. In the code, we do so by replacing with , where is defined to be the growth on the scale , which corresponds to the largest mass value computed by the code inside the halo model routine. This is chosen to be consistent with the definition of .

### iii.5 CDM reference model for GDM halo model

In order to implement the halo model as described above, we need a reference CDM cosmology when class is run for GDM (this is used to calculate the half mode mass; see equation 6). For this, we use the Eisenstein-Hu fitting formula, for a cosmology with the same , , , and as the GDM cosmology. The two main references for the fitting formulas are Eisenstein and Hu (1998) and Eisenstein and Hu (1997). The first of these includes the effect of baryon oscillation but not neutrinos, whereas the second takes an average over the oscillations but includes the damping effects of massive neutrinos. We implement the former of these, as used for HMCODE Mead et al. (2015) to calculate the CDM spectrum at the - and - values required in class, which are then stored in an array. This spectrum is normalised to have the same value as the class GDM spectrum at , as this scale should be above the scales that are affected by GDM.

The CDM growth function from Lahav et al. (1991) is used as part of the Eisenstein-Hu formulas,

(24) | |||||

(25) | |||||

(26) |

Here, is a parameter from the Eisenstein-Hu formulas and is included for completeness, however note that it is irrelevant once the growth factor is normalised to unity today. The reference wavenumber is set to be the largest wavenumber in the table used by class that is less than , and is set to be the smallest redshift in the table used by class that is greater than 50.

### iii.6 Comparison to WDM and FDM halo models

#### iii.6.1 Theoretical comparison

In Barkana et al. (2001) the cutoff of the mass function at small masses for WDM was achieved by an additional mass dependence of the barrier (see also Marsh (2016); Benson et al. (2013)). This mass dependence of (a steep increase for small masses) was argued to be caused by the velocity dispersion, however it is unlikely that this is the true physical mechanism that suppresses the mass function below since WDM simulations have shown that the velocity dispersion is irrelevant for the large scale structure and the mass function Schneider et al. (2013); Viel et al. (2012). This disparity was explained in Schneider et al. (2013) (p. 4 last paragraph before section 3) by splitting the effects of the velocity dispersion into two distinct time periods: the accumulated effect from times , (which manifests in the usual linear theory matter power spectrum cutoff), and the late time velocity dispersion (as should be present but turns out to be negligible in N-body simulations).

In our GDM halo model, we have allowed for the possibility of both a cut-off of the matter power spectrum due to accumulated effects with , and a steepening of the barrier due to effects related to times . Physically, , equation (16), takes into account pressure and viscous effects that hinder collapse at , whereas is the integrated effect due to . If is chosen during matter domination (as in our halo model) then only the latter (integrated) effect matters for WDM, and the cutoff in the WDM mass function must originate independent of late time velocity dispersion effects on , since it is observed in simulations without any added velocity dispersion as in Schneider et al. (2012). Thus the correct implementation of the mass function cut-off^{1}^{1}1In Schneider et al. (2012) it was argued that this mass function cutoff is due to a non-hierarchical structure formation for masses below that is in conflict with the excursion set picture. However, it might be possible to show that this effect is the same as the strongly non-Markovian random walk, responsible for a mass function cut-off in Paranjape et al. (2012), such that the cut-off can be understood within the excursion set theory. due to non-Markovian behaviour caused by the linear theory power spectrum cutoff is not via the steep increase of the barrier for small masses when is chosen during matter domination. Instead it manifests through the phenomenological prefactor that is present in equation (23), which depends on .

Furthermore, when the WDM halo model is expressed in terms of our halo model, , i.e. the spherical collapse threshold reduces to the scale-independent CDM spherical collapse threshold. This is a special (approximate) property of WDM, caused by the dependence of pressure and viscosity that we do not expect to apply to more general forms of the GDM model, which is why we allow our halo model to have both a cut-off due to integrated earlier time behaviour and a later time change to the barrier.

Note that the mass dependence of we have introduced in section III.3 for GDM, equation (16), is due to the evolution of the shape of the linear power spectrum after . There is (approximately) no shape change for WDM and thus our mass function is very similar to the one used in Schneider et al. (2012). The difference is that we allow for evolution of the shape of the power spectrum after , which is expected for constant GDM parameters, and thus we have mass dependent , which then requires using a mass function that can deal with mass dependent barriers.

#### iii.6.2 Qualitative numerical comparison

In principle, we can compare our non-linear matter power spectrum to the spectra in literature (e.g. Smith and Markovic (2011), Viel et al. (2012), Schneider et al. (2012) and Marsh (2016). However, we note that we focus our work here on constant values of the GDM parameters, whereas WDM and FDM correspond to time and scale dependent parameters. Thus the impact of the non-linearities will be different, although we can nonetheless perform a qualitative comparison of how the halo model affects the predictions of our GDM model compared to how it affects the predictions for the specific cases of WDM and FDM.

For all models, we see the same qualitative behaviour that the non-linear corrections increase the matter power spectrum and reduce the differences between CDM and the modified matter content compared to the linear theory. However, there are differences in detail. For example, consider WDM with a mass of 0.25keV, where meaningful changes to the linear spectrum compared to CDM begin to occur on scales of , and these changes begin on even smaller scales as the mass increases. For FDM, in line with Marsh (2016), the changes are on even smaller scales, similar to WDM with mass 1keV. Whereas, at , our GDM models consistent with Planck constraints differ from CDM by two orders of magnitude in linear theory. This means that the GDM linear spectrum differs from CDM on scales larger than those where the non-linear corrections matter for CDM, whereas these two scales are swapped for the WDM and FDM models studied in the literature. This is because we have time-independent GDM parameters, whereas WDM would correspond to having them decay as , which causes their effects to appear only on small scales where nonlinearities are important. We expect that, once we allow for time dependent GDM parameters, the halo model will make a bigger difference relative to the linear spectrum, due to this reduced suppression in the linear theory at late times.

## Iv Description of datasets and methodology

In this section we explain the data and methodology that were used to generate our results. We used the class code Blas et al. (2011), modified as detailed in Thomas et al. (2016); Kopp et al. (2016), to evolve the GDM perturbation equations. We have added a module implementing the halo model as described in the previous section.

Our parameter constraints were obtained using the same basic methodology as in Thomas et al. (2016), see there for further details. We used the MCMC code MontePython Audren et al. (2013) and established convergence of the chains using the Gelman-Rubin criterion Andrew Gelman (1992). We constrain a 6 parameter CDM model , where is the density of the dark matter fluid, which is CDM in the CDM case and GDM otherwise. We set uniform priors on and such that . The helium fraction was set to Planck Collaboration et al. (2015b) and we assumed adiabatic initial conditions. We used two massless and one massive neutrino with mass eV keeping the effective number of neutrinos to (thus for simplicity we refer to “neutrino mass” during the analysis, although this is equivalent to the sum of the neutrino masses for this choice of parameters). The base parameter set is augmented by 3 GDM parameters for the GDM runs, and additionally also the neutrino mass (for the single massive neutrino species) for some runs.

We perform runs for CDM that are purely linear, linear+halofit and linear+halo model (where the halo model is as documented in the previous section, which is why it is important that our GDM halo model reduces to CDM for vanishing GDM parameters). The halofit Smith et al. (2003); Takahashi et al. (2012) runs are performed using the halofit model built into class. For GDM, we will perform purely linear runs and linear+halo model runs, these runs will be referred to as “HM” in the results table.

Our primary dataset is the Planck 2015 data release Planck Collaboration
et al. (2015a) of the CMB anisotropies power spectra, consisting of the low- T/E/B likelihood and the TT/TE/EE high- likelihood with the full “not-lite” set of nuisance parameters.^{2}^{2}2For full details, see the Planck papers and wiki http://wiki.cosmos.esa.int/planckpla2015/index.php/. These likelihoods combined are referred to as Planck Power Spectra (PPS). We also use BAO^{3}^{3}3In appendix B we examine a possible subtlety with the use of BAO data, which would also be relevant if GDM were constrained using redshift space distortion data. from 6dF Galaxy Survey Beutler et al. (2011) and the Baryon Oscillation Spectroscopic Survey Sloan Digital Sky Survey Anderson et al. (2014) (collectively referred to as BAO hereafter), and the Planck CMB lensing likelihood (Lens).

The key additional dataset that we use here is the WiggleZ matter power spectrum Parkinson et al. (2012) (referred to as MPS in the results table). This includes galaxy power spectrum measurements at four redshifts, . We follow the procedure laid out in Parkinson et al. (2012) for the likelihood from this data, as implemented in MontePython, with the following exceptions. We do not use the WiggleZ giggleZ non-linear prescription, as this is not valid for GDM. Instead, we will use both halofit and our halo model for CDM runs, as detailed above. For GDM, we will perform purely linear runs and runs using our halo model. Note that the WiggleZ likelihood processes the input theory spectrum, including convolving with the window function and other transformations. In particular, an analytic marginalisation over the linear bias is performed, see Parkinson et al. (2012) for details. We also consider two different subsets of the whole WiggleZ data: a conservative cut, where we use the provided k-bands up to Mpc, and a less-conservative cut using the complete data up to Mpc, the latter of which was used by the WiggleZ collaboration for their CDM results.

## V Constraints

We divide our constraints into two groupings: those without MPS data that focus on examining the robustness of previous constraints, and those using MPS data that aim to improve the constraints on GDM. The main results from these two groupings are that the previously obtained constraints are indeed robust and that the MPS data improves the constraints on and by a factor of three. The constraints from the different runs can be found in table 1.

### v.1 Robustness of previous results

Our first result follows from looking at table 1. Here we show our previous constraints on the GDM parameters using the data combinations PPS and PPS+Lens, both with and without the inclusion of the halo model (HM). It is clear that the constraints are essentially independent of the halo-model correction. This is because the halo model implemented as described above only has a small effect in the GDM matter power spectrum relative to the purely linear theory for the scales relevant to the upper limits of the constraints (see figures 1 and 2). This is partly due to the choice of constant GDM parameters; as explained above we expect that the impact relative to the linear spectrum would be more important if we chose time dependent forms for the parameters, e.g. an time dependence, or general binned functions. The relative effect of our halo model on the linear spectrum is decreased as the and parameters are increased, as can be seen by comparing the GDM and CDM halo model curves in figure 2. This is caused by the strong effects on the linear spectrum from the phenomenology, which dominate over any changes to the matter power spectrum due to the non-linear effects. Thus, we expect that any constraints on and obtained from the CMB temperature, polarisation and lensing spectra on these scales are actually more robust to potential non-linear complications than the standard CDM parameters.

Likelihood Model | (upper bound) | (upper bound) | ||||
---|---|---|---|---|---|---|

(PPS+…) (-GDM+…) | ||||||

+ Lens | ||||||

+ Lens + BAO | ||||||

+ HM | ||||||

+ Lens + HM | ||||||

+ Lens + BAO | ||||||

+ Lens + BAO + MPS () | ||||||

+ Lens + BAO + MPS ()+ HM | ||||||

+ Lens + BAO + MPS () | ||||||

+ Lens + BAO + MPS ()+ HM |

Also in table 1 we show the constraints for the data combination PPS+Lens+BAO, when the neutrino mass is both fixed or varied (). The differences between the posteriors with and without the inclusion of the neutrino mass can be seen in figure 3. The perturbative GDM parameters and are affected little by the inclusion of the neutrino mass, indeed the constraints improve very slightly, whereas the inclusion of the neutrino mass does noticeably worsen the constraints on the equation of state . These effects are both caused by the degeneracies between the neutrino mass and the GDM parameters, which can also be seen in figure 3; we shall now explore these in more detail.

The neutrino mass correlates with and in the same way that they are correlated with each other: The neutrino velocity dispersion is Lesgourgues and Pastor (2006), which causes a reduction in the lensing potential just like and . This is not surprising since massive neutrinos can be described by a GDM fluid Blas et al. (2011). This similarity between the effects of the neutrino mass and these GDM parameters can be seen by comparing the second and third panels on the left in figure 4, showing the ratio of the spectra to CDM with and without the lensing contribution. This is also shown in the right hand plot of figure 4, showing that all of these parameters result in a substantial reduction to the CMB lensing potential spectrum. The insets of figure 3 show the 3D posterior of , and . The lower of the two plots is colour coded according to the probability density, which peaks in the CDM corner, as expected due to the lack of a detection of GDM parameters and neutrinos mass. The upper of the two insets shows the 90% confidence level contour in orange, which is shown to be well modelled by constant . This is in rough agreement with the expression for , as expected if the degeneracy is indeed due to the reduction of the lensing potential described here.

Note that the geometry of this situation is slightly non-trivial: we are dealing with a “corner” of a multi-dimensional parameter space where the parameters are all required by physics to be non-negative. As noted in Heavens and Sellentin (2018), marginalising over parameters in such situations can cause some subtle effects, and this may be the source of the slight improvement to the constraints on and when the neutrino mass is included and marginalised over.

Figure 3 also shows a clear degeneracy between and the equation of state , which is due to the ability to generate cosmologies with identical but different when . This can be seen by comparing the different panels in figure 4, which show the CMB temperature power spectrum for different sets of parameters. Comparing the increased neutrino mass cosmology (red dashed) line in the third and fourth panels on the left, it can be seen that for fixed (the angular scale of the acoustic oscillations) and the main effect of the increased neutrino mass on a CDM cosmology (aside from the lensing effect discussed above) is a reduction in the ISW effect and a tilt for higher , see the “no-lensing” panel. The reduction of ISW (both early and late) as caused by the increased abundance non-relativistic matter content (compared to radiation and cosmological constant) when the neutrino mass is increased. One cannot simply compensate this with a change to , because this would adversely affect expansion history at early times. However, when the parameter is introduced, it is possible to vary , whilst also changing (and ) to keep the expansion history approximately fixed. This allows a GDM cosmology with increased neutrino mass to have the same ISW effect and high- tilt of as a CDM cosmology with lower neutrino mass, as can be seen by comparing the red (dashed) and blue (short-dashed) lines in the third panel of figure 4. This ability of to counteract these two effects of increasing drive the degeneracies between and .

Note that the degeneracies between and the GDM parameters mean that if tighter constraints are put on the neutrino mass from other experiments, then using these results as a prior on CMB analysis could further improve the constraints on the GDM parameters.

### v.2 Use of MPS data

#### v.2.1 Conservative cut - No detection of GDM

In the lower half of table 1 we show the constraints obtained when including the WiggleZ matter power spectrum data (MPS), for two ranges of wavenumbers: and . We discuss the former (more conservative) of these first. Note that for CDM, the WiggleZ team found that the linear theory appears to give a better fit than halofit (see figure 3 in Parkinson et al. (2012)), and we obtain the same result both for halofit and the halo model, although the halo model has a smaller effect on the large scales than halofit. We will focus the discussion on the perturbative parameters and , because the inclusion of the WiggleZ data has little effect on the constraints on . This is because the effect of for the scales under consideration is primarily a small change to the amplitude, with little change to the shape of the spectrum, see Kopp et al. (2016). Thus, the effect of will be almost entirely removed by the marginalisation over the bias in the WiggleZ likelihood, although note that future matter power spectrum data could constrain due to its effect on the peak location.

The inclusion of the matter power spectrum with the conservative cut has a strong effect on the perturbative parameters and , improving the constraints by a factor of three, see table 1. This can also be seen in figure 5, where the 1D posteriors narrow considerably when the WiggleZ data is included. As this figure shows, this improvement is the only significant change to the posteriors due to the extra data. The green contours in this figure are from the linear implementation of GDM, however we note that the changes due to the extra data are approximately independent of whether the linear or halo model implementation of GDM is used; this can be seen both from the constraints in table 1 and in the inset to figure 5 by comparing the green (linear) and black dotted (halo model) contours. This is important because it implies that the constraints using the conservative cut come from physics that is well understood, and is not sensitive to detailed considerations of the non-linear regime.

Left: GDM curves with and non-zero (but note that the plot for and non-zero is essentially identical, since the degeneracies between these parameters hasn’t been broken). The (dashed) blue and green curves correspond to the linear and halo-model predictions respectively, for parameters corresponding to our previous (linear) constraints in Thomas et al. (2016) (). The (solid) orange and yellow curves correspond to the linear and halo-model predictions respectively, with GDM parameters corresponding to the improved constraints in this paper when the MPS data is used with the conservative cut (). It can be seen that the previous constraints have some tension with the WiggleZ data, and thus that its inclusion improves the constraints on the GDM parameters by requiring smaller values of the parameters to reduce the tension. For all parameter choices here, the difference between the linear and non-linear spectra is small.

Right: Theoretical spectra constructed using the best fit parameters from the MCMC runs with the less-conservative cut for the WiggleZ data. The spectra correspond to linear CDM (black dot-dashed), CDM with the halo model (orange dotted), linear GDM (blue dashed) and GDM with the halo model (green solid). The linear GDM best fit parameters are and and the GDM plus halofit best fit parameters are and . Similarly to WiggleZ Parkinson et al. (2012), we get a better fit for linear CDM than when the halo model is included. The best fit model for linear GDM has small GDM parameters and a spectrum that is very similar to the CDM spectrum, as expected since the constraints in this case are consistent with CDM. For the GDM halo model case, the best fit model has a much larger sound speed, and thus deviates more from the CDM best fit.

The improvement on the and constraints is due to the decay of the matter power spectrum for , which creates a slope to the matter power spectrum that is inconsistent with the data for larger values of and . This can be seen in the left panel of figure 6, which compares theoretical spectra with various upper limits of GDM parameters to the lowest redshift bin in the WiggleZ data. Note that in both panels of this figure the theory spectra are transformed in line with how the likelihood is computed, in order to be compared to the data. This includes convolving with the survey window function, marginalising over the linear bias and accounting for the difference in background to the fiducial cosmology. See Parkinson et al. (2012) for a full description and explanation of these processes. The left hand panel of figure 6 compares GDM spectra computed with non-zero values corresponding to upper limits of previous constraints (“old”) and the constraints when the matter power spectrum data is included (“new”). There is a tension between the slope of the theoretical spectra and the data that is reduced when the lower value (associated to the constraints from the WiggleZ data) is used. The equivalent plot for non-zero is essentially identical, since the degeneracy between these parameters has not been broken (see below). In accordance with what was noted earlier, the use of the halo model makes little difference to the theoretical spectra in figure 6, and thus to the constraints obtained with and without the halo model. The constraining power here does not come from the measured amplitude of the matter power spectrum because of the marginalisation over the linear bias. We expect that these constraints could increase even further if the matter power spectrum is measured on larger scales, particularly the turnover around the peak. This would have the additional advantage of staying inside our conservative regime that we have seen is robust to non-linear modelling. Considering how the difference in slopes shown in figure 6 continues for , we expect that the less-conservative cut for the WiggleZ data will improve the constraints further, see below.

One of the motivations for considering matter clustering data for GDM constraints is the attempt to break the degeneracy between and caused by the phenomenology. The scales we are looking at with the WiggleZ data here are insufficient to break this degeneracy, see figure 7: up to , there is little difference between and . Note that from Kopp et al. (2016) we expect the value where the oscillations start to be , which is in agreement with what we see here. Even on scales down to (i.e. to the level of our less-conservative cut; see below), the difference between the spectra generated by the two parameters is not large, although beyond the difference between the linear spectra increases substantially. Interestingly, the non-linear spectra corrections act to recreate the degeneracy between the and spectra on scales below approximately (see figure 7). If this modelling of GDM non-linearities is accurate then this means that there is only a small range of scales (around ), in which data could allow us to distinguish the effects of these two parameters.

#### v.2.2 Non-conservative cut - Possible detection of GDM

As mentioned above, we also consider a less conservative cut, with . The constraints are presented in the final two lines of table 1, where it can be seen that a significant gulf opens between the constraints obtained with and without the halo model. The constraints including the halo model improve by a factor of 3 compared to the more conservative cut, amounting to a combined improvement compared to previous constraints of an order of magnitude. However, the constraints without the halo model are another factor of 3 or so stronger still. This weakening of the constraints due to the inclusion of the halo model naively seems to match our expectations of weaker results when the non-linear effects are included, however there is a deeper story here.

The inset to figure 5 shows the 2D posterior contour plots for and for the conservative and less-conservative cuts, for both linear and halo model implementations of GDM. As discussed above, the linear and halo model versions of GDM result in very similar constraints for the conservative cut to the matter power spectrum data. The less-conservative cut for linear GDM results in a similar looking set of contours, in the sense that the contours are all essentially right-angled triangles with a similar slope on the side joining the two axes. I.e. the contours are solely upper bounds with a particular slope and the two GDM parameters both being zero is consistent with the data. The only difference is that the upper bounds have been significantly reduced.

However, the shape of the less-conservative cut for GDM with the halo model is significantly different; the lower contour is now apparent, resulting in a trapezoidal contour and a clear inconsistency of the CDM point ( and ) with the contours. There is now a clear preference for a non-zero GDM parameter. This difference between the linear and halo model GDM results is further seen by looking at the right hand panel of figure 6. Here we plot the spectra generated using the best fit parameters from the different MCMC runs. The linear GDM best fit parameters are and and the GDM plus halofit best fit parameters are and . The best fit model for linear GDM has small GDM parameters and a spectrum that is very similar to the CDM spectrum, as expected since the constraints in this case are consistent with CDM. For the GDM halo model case, the best fit model has a much larger sound speed, and thus deviates more from the CDM best fit. Thus we can see the discrepancy between the linear GDM and halo model GDM manifesting here as well. Interestingly, the GDM halo model run returns a best fit spectrum that is closer to the linear CDM spectrum than the halo model CDM spectrum is. The for the two GDM curves and the linear CDM curve are almost indistinguishable,^{4}^{4}4We should caution here that an MCMC code such as MontePython is not optimised in terms of finding the lowest possible value of the likelihood, so we expect that the best fit values we have found here are close to the absolute minimum, but not precisely the lowest values. showing further that the inclusion of the halo model into GDM can compensate for higher values of the GDM parameters.

Despite the preference for a non-zero GDM parameter shown in the GDM+HM contours in the inset of figure 5, the individual 1D posteriors (not shown here as their interpretation is dubious; see below) show no preference to be non-zero due to the degeneracy between the two parameters. This means that marginalising over the other parameter results in no “detection” of a non-zero value of either parameter. Despite the resulting large difference in upper bounds between the less-conservative results with and without the halo model, the maximal width of the contours along the line marked in the plot is similar in the two cases.

The slope of the degeneracy in these contours is well understood as the direction along which remains fixed, and it is the direction perpendicular to this (the line marked in the plot) that is most constrained by the data. In principle, it would be possible to create a new parameter that describes this direction (i.e. essentially ), and then the second coordinate (i.e. , ) of this 2D space can be marginalised over, in order to create a 1D posterior for . If we had chosen uniform priors on and , we would expect the 1D-posterior of this new -parameter to peak at non-zero values for the yellow contour in Fig.5, and to peak at zero for the blue contour. And thus this parameter would allow to compactly quantify the detection of perturbative GDM parameters. However, we note that our current choice of priors (uniform priors on and ) would make peak at non-zero value for any sensible 2D posterior, in particular also for the blue contour in Fig.5. We will explore these issues in a forthcoming paper.

We note that this difference between the linear and non-linear results for the less-conservative cut is a strong justification of the motivation behind this paper, namely that correctly modelling these non-linear scales will be crucial for using late-time clustering data to constrain the GDM parameters. At the level of the modelling we have done here, we cannot be sure of our results with either the linear or non-linear modelling. Instead, we take these results to show that there is a need to look into the non-linear modelling in substantially more detail, before a detection of non-zero GDM parameters using late time clustering data could be claimed. Furthermore, even if the less-conservative linear contours agreed with the non-linear contours, then we would be hesitant to claim a detection of non-zero GDM parameters because of the caveats related to marginalising over parameter subspaces with special geometries. Nonetheless, our results show that late time matter clustering data can strongly constrain the GDM parameters, to the level where a detection is possible.

## Vi Conclusion

We have investigated how considerations around large scale structure can affect constraints on the generalised dark matter parameters. The main results of this work are the development of the halo model for GDM presented in section III and the improved constraints on the GDM parameters presented in table 1.

In section III we argued for modifying the CDM halo model in a particular way, by backtracking the “linearly-extrapolated” critical density for collapse. This allows the mass dependence of the collapse barrier to be implemented in a natural way. This halo model reduces to a standard CDM halo model in the case of scale independent growth, and produces qualitatively similar results to halofit. Having derived the halo model for GDM, we note that the non-linear corrections are much less significant for GDM (with constant and ) than for CDM, because the strong linear decay dominates over the corrections from the halo model.

We use this halo model to test the robustness of previously obtained constraints based on Planck CMB power spectra, as seems circumspect considering the magnitude of the CDM non-linear corrections and difference between CDM and GDM spectra in figure 1. We find that the GDM constraints change little, as expected from the aforementioned phenomenology of the GDM halo model. Interestingly, the perturbative GDM parameters and are less sensitive to the non-linear corrections than the standard CDM parameters, which will be increasingly important for future CMB lensing surveys, such as Simons Observatory The Simons Observatory Collaboration (2018). We additionally checked the changes to previous GDM constraints when the neutrino mass is allowed to vary as a free parameter, primarily finding a worsening of the constraints on the equation of state (of GDM) , due to the degenerate effects on the expansion history. We also elucidate the three-way degeneracy between , and , and note that the geometry of this situation requires that marginalisation over these parameters is done carefully Heavens and Sellentin (2018).

We examined the effect of including the WiggleZ matter power spectrum data when constraining the GDM parameters, finding a factor of three improvement on the sound speed and viscosity constraints when a conservative cut in wavenumber is used. When increasing the -range that is included, these constraints improve by a further factor of three, for a total improvement from the use of matter power spectrum data of an order of magnitude. This shows the value of datasets that constrain the phenomenology of GDM. Since we analytically marginalise over the linear bias, we expect that once large scale structure measurements reach the peak of the matter power spectrum, the total constraining power will be sufficient to either constrain the GDM parameters to the point of cosmological irrelevance or yield a detection of beyond CDM physics.^{5}^{5}5Note that, being more precise, this “beyond CDM physics” could just be a more careful and precise modelling of the CDM universe, according to the EFTofLSS interpretation of GDM. These improved constraints are one of the key results of this work.

The results from extending the -range that is included show some important features for future work. The first is that there is a difference between the constraints obtained by linear and non-linear modelling, thus showing the importance of robust non-linear modelling of the GDM model for the use of late time matter clustering data. This suggests that a re-evaluation of the tight constraints obtained in Kunz et al. (2016) could be interesting. This also shows that future surveys, combined with a more detailed analysis of the non-linear completion of GDM, have the potential to show that the matter power spectrum is not consistent with the CDM. Furthermore, we have shown that GDM with the halo model yields a 2D contour that is clearly inconsistent with the CDM point (), although we advise caution in the interpretation of this due to the difference between the linear and non-linear results. Even in this case, neither nor is individually detected, due to the degeneracy between these two parameters, and we have noted that it is not straightforward to quantify a detection in such cases, because marginalisation in areas of parameter space with corners can lead to biases. Doing so requires a careful analysis of priors, see e.g. Heavens and Sellentin (2018). We plan to explore this issue further in the future.

As the parameters in the halo model were originally calibrated for WDM, it may be preferable to treat them as nuisance parameters and vary them in an MCMC analysis for GDM. However, this further increases the computational demands of the already expensive codes and is unlikely to make a difference to our results as we have made no detection of GDM. Furthermore many other aspects of the halo model, like the precise form of the spherical collapse barrier, remain currently an educated guess and require a detailed study using numerical simulations of a suitably defined non-linear GDM model. Thus we leave an investigation into these issues for future work.

Throughout this work we used the simplest parameterisation of the GDM parameters; a single value with no redshift or scale dependence. We expect the halo model to have a larger impact for time dependent GDM parameters, although we leave this investigation to future work. In particular, we expect that for , corresponding to WDM and FDM, that the halo model has a large impact. The assumption of constant GDM parameters was relaxed in Kopp et al. (2018a), where the equation of state was measured in multiple redshift bins. The importance of the WiggleZ data considered here for the constraints on and