Cosmological evolution of matter with interacting dust fluids

We split the total matter fluid into a bound (halo) component and an unbound (free particles) fluid component that is accreted by the halos. We adopt a different framework that treats the structure formation problem as a gravitational interaction between these virialised cold dark matter halos and the unbound inter-halo cold dark matter (and cold baryon) particles. This interaction involves in general an exchange of energy and momentum during the accretion process. We then explore the evolution of the average matter density and of large-scale structure formation, using a simplified phenomenological model that is based on results from extended Press-Schechter and N-body simulations. At high redshifts most matter is in diffuse form and is not part of the halos. As particles are accreted by the virialised halos, the particle number density decreases and that of the bound matter increases. We also present a general analysis of the background and linear perturbations for the interacting fluids, showing in detail the energy and momentum exchange terms.


Introduction
In the standard ΛCDM cosmological model, at late times when radiation may be neglected, the observed large-scale structure in the universe is understood to be based on cold dark matter (CDM) halos which act as gravitational sinks for baryonic matter, leading to the formation of stars, galaxies and higher-level clustering of luminous matter.Here we consider a simplified matter model in which the CDM is an elementary particle (e.g. the neutralino) and the baryons are cold.We neglect the complexities associated with non-cold baryonic matter by confining our analysis to linear scales.In this scenario, the first collapsed bound (virialised) objects are CDM halos of mass determined by the mass of the CDM particle -e.g., approximately Earth mass for a neutralino particle [1].Subsequently, halos grow by accretion of the unbound interhalo particles.Moving forward in time, an increasing fraction of the total matter is captured in bound structures, with a spectrum of masses [1][2][3][4][5][6].
Analytical models of structure formation are based on a Friedmann-Lemâitre-Robertson-Walker (FLRW) background with perturbations that grow under gravitational instability.The implicit assumption involved here is that the evolving mixture of particles of different (and changing) masses may be treated as a single pressure-free 'dust' model, all of whose particles are considered as test particles -i.e.effectively massless particles.Standard perturbation theory does not specify the masses of the dust particles and their changing mass spectrum.The agreement between the results of the perturbative analysis and those of N-body simulations on large enough scales suggests that this assumption is reasonable.Nevertheless, it is interesting to ask what the nature of the evolving spectrum of masses is and whether and how it affects cosmic evolution and structure formation.This question can be answered by N-body simulations [1][2][3][4][5][6] but an analytical approach is also of interest.In this paper we adopt a simplified analytical approach to tackle this question (see e.g.[7][8][9] for more sophisticated approaches).
Studies of interacting fluids (see e.g [10] for the theoretical foundations and [11,12] for recent work) have so far mainly focused on the "dark sector" interaction between dark matter and dark energy and have provided useful insights into the coincidence problem and structure formation, within general relativity. 1 A major difficulty confronting these models is the absence of a physically well-motivated non-gravitational interaction.The interaction that we propose here is not intended in any sense as non-gravitational energy-momentum exchange -on the contrary, it is the purely gravitational interaction of accretion.
In other words, we treat the structure formation problem as a gravitational interaction between virialised CDM halos and the unbound inter-halo CDM (and cold baryon) particles.This interaction is simply an exchange of energy and momentum during the accretion process.However, given the nonlinear complexities of accretion in structure formation, we use a simple phenomenological model that is based on an extended Press-Schechter analysis, tested with N-body simulations.To this end, we use a simple fit to the results from [1].These results indicate that by redshift zero, ∼ 80−95% of matter is within halos of all possible masses -with ∼ 60−70% of matter in halos of mass ≳ 3×10 9 M ⊙ (see Fig. 3 in [1]).The percentages depend on the details assumed in their models.These details are not important for our purposes, since our focus is on a simple analytical model that reflects the key qualitative features.

The 2-fluid model
The two fluids are the halo fluid A = h and the free (unbound) particle fluid A = f , each with a dust equation of state, w A = 0, and vanishing speed of sound c sA = 0.The total matter density is which satisfies the background and perturbation equations of the standard ΛCDM model.For example, in the background where the prime denotes a conformal time derivative and H = a ′ /a is the conformal Hubble rate.
We assume a simplified model for the density transfer from free particles to halos: In other words, the fraction of mass in halos is given by a function of the total density.
The general background equations for two interacting fluids (see e.g. the extensive review in [10]) reduce in our case to Here QA are the rates of energy density transfer: Clearly Qh is positive since the halo fluid is accreting the free-particle fluid.We can also rewrite these equations in terms of the dimensionless density parameters:

Linear perturbations
Linear perturbation of (3) gives From this we obtain It follows from these equations that the power spectra of the halo fluid and the freeparticle fluid are where P m is the matter power spectrum.Equations ( 9)-( 15) allow us to evade the solving of the perturbed energymomentum conservation equations.Nevertheless, it is instructive to see how our equations correspond to the general perturbation equations for coupled fluids.In order to make this comparison, we choose the Newtonian gauge.The perturbed line element at late times is then where Φ is the gravitational potential.The four velocity of fluid A is where υ A is the peculiar velocity potential, which vanishes in the background.The general perturbation equations for two interacting fluids are based on [10] where Q µ A is the energy-momentum transfer four-vector, which satisfies A Q µ A = 0 from conservation of total energy-momentum.In general the four-vector Here F µ A is the momentum density transfer rate relative to the total four-velocity u µ = a −1 (1 − Φ, ∂ i υ) and f A is the intrinsic momentum transfer potential.It follows that In the case of 2 dust fluids, it is natural to choose the total 4-velocity as that of matter: υ = υ m .Then the equations reduce to The total matter density contrast evolves according to the equation Our simplified model does not include the small effect of the gravitational potential on the density contrast.Since |Φ| ≪ |δ A | on the scales of interest, we can neglect the gravitational potential in ( 22) and ( 24), so that Since our model is confined to scales where linear perturbation theory is valid, we can ignore the velocity bias between the halos and the unbound particles.This bias arises on small nonlinear scales because the halos represent peaks in the matter distribution, with a threshold for formation [7,9,15].Therefore it is reasonable to assume Consequently, there is no intrinsic momentum transfer in the matter frame: With these additional simplifications, the perturbed conservation equations reduce to By comparing ( 11)-( 13) with ( 29)-(30), we find that Here we used δQ f = −δQ h , and the linear growth rate is defined as

A simple model
A simple power-law fit to the results of [1] (their Fig. 3, ellipsoidal collapse model) is where the constants are F 0 = 0.9 and α = 0.07.Note that in the background, F = F 0 (1 + z) −3αz , where z is the cosmological redshift.The background equations ( 9

and (10) now become
The numerical integration of these equations leads to the background energy density parameters for this model, which are displayed in Figure 1.The plot shows how the matter density of the bound structures grows while that of the unclustered matter decreases over cosmic time.
The key features of the plot are as follows.The high-redshift growth and decay in the matter densities happens until Ω h -Ω f equality is reached at z ∼ 13 (a ∼ 0.07), when Ω h = Ω f ∼ 0.5, which is when the clustered matter existed in the same quantity as the unclustered matter.Past this time, Ω h continues to grow while Ω f continues with its decline.The halo fluid reaches its highest density (Ω h ∼ 0.8) when z ∼ 5.The total matter energy density, Ω m , remains nearly constant, up to this point, whereafter it starts to decay due to the effect of dark energy, before reaching its current value of Ω m0 .
Plots of the interaction function (34) are displayed in Figure 2 where the behaviour of the mass halo fraction with respect to the redshift and the scale factor over cosmic time can be observed.The fraction of all matter predicted to be in halos at z = 0, i.e. to be part of clustered matter, is close to hundred per cent.However, at high redshifts or early times a much higher fraction of cold dark matter and cold baryon particles is in diffuse form and is not part of the halos.

Discussion
We presented a new phenomenological model to describe the evolution of structure formation by treating the matter as composed of two pressure-free components with a purely gravitational interaction: a free-particle fluid (cold dark matter and cold baryons) that is accreted onto the halos, which are the 'clustered' particles of the halo fluid.We derived the background and perturbation equations for the case of two interacting fluids and solved them for the model in which the halo mass fraction F, i.e. the ratio between the bound and unbound matter, is a power law, approximating the results from the extended Press-Schechter formalism and N-body simulations [1].
The clustered matter grows due to accretion of cold particles, while the unclustered matter density decays over cosmic time.On the other hand, the total matter density, Ω m , follows the standard evolution, remaining nearly constant for a large part of the cosmic evolution, and decreasing at late times, due to the effect of dark energy.
The halo mass fraction, F, grows over cosmic time.This means that the fraction of all matter predicted to be in diffuse form is much larger at early time while that predicted to be in gravitationally bound structures is much larger at late times, close to hundred per cent at the current epoch.At low values of a, i.e. high redshift, a large fraction of cold particles is free, up to Ω h -Ω f equality at z ∼ 13, beyond which a much higher fraction of cold particles has become part of halos.These results show that in the framework of hierarchical structure formation, gravitationally bound structures can form not only from the assembly of smaller halos, but directly from free diffuse cold particles (see also [1]).
Our simplistic model is of course no substitute for nonlinear perturbations and simulations, but it does allow for an alternative perspective on some of the general properties of large-scale structure formation.For example, Figure 1 also reveals the existence of a maximum in halo density at z ∼ 5. Since galaxies are hosted in halos this feature could be tested against high-redshift galaxy observations, such as those currently underway with JWST (see e.g.[16]).

Fig. 1
Fig.1The background density parameters of the bound (halo) fluid h and the unbound (freeparticle) fluid f , shown together with the total matter density Ωm, for the model (34)-(36).