A&A 388, 741-757 (2002)
DOI: 10.1051/0004-6361:20020548
The phase-diagram of cosmological baryons
P. Valageas1 - R. Schaeffer1 - J. Silk2
1 - Service de Physique Théorique, CEN Saclay, 91191 Gif-sur-Yvette, France
2 - Astrophysics, Denys Wilkinson Building, Keble Road, Oxford OX1 3RH, UK
Received 28 November 2001 / Accepted 5 April 2002
Abstract
We investigate the behaviour of cosmological baryons at low redshifts
after reionization, through analytic means. In particular, we study the
density-temperature phase-diagram that describes the history of the
gas. We show how the location of the matter in this
diagram
expresses the various constraints implied by usual hierarchical
scenarios. This yields robust model-independent results that agree with
numerical simulations. The IGM is seen to be formed via two phases: a
"cool'' photo-ionized component and a "warm'' component governed by
shock-heating. We also briefly describe how the remainder of the matter is
distributed over galaxies, groups and clusters. We recover the fraction of
matter and the spatial clustering computed by numerical simulations. We also
check that the soft X-ray background due to the "warm'' IGM component is
consistent with observations. We find in the present universe a baryon
fraction of 7% in hot gas, 24% in the warm IGM, 38% in the cool IGM, 9%
within star-like objects and, as a still un-observed component, 22%
of dark baryons associated with collapsed structures,
with a relative uncertainty no larger than 30% on these numbers.
Key words: cosmology: theory - large-scale
structure of Universe - galaxies: intergalactic medium
1 Introduction
As is well known, the mass of baryons embedded within stars or galactic disks
in the current universe is quite small since it yields a baryonic parameter
(e.g., Fukugita et al. 1998) while standard
nucleosynthesis calculations give
(e.g.,
Tytler et al. 2000). Therefore, most of the baryonic matter should lie in the
intergalactic medium. This agrees rather well with the fact that at higher
redshift one observes a large amount of mass in the Lyman-
forest
which consists of moderate density fluctuations ionized by the background UV
flux emitted by distant galaxies. Thus, one gets
at
(e.g., Fukugita et al. 1998). However, as noticed in
Cen & Ostriker (1999) the mass within the Lyman-
forest decreases with time
and at z=0 summing over all observed contributions one obtains
which falls short of the required total baryonic mass. Hence at z=0
a large part of the baryons must lie in a new intergalactic component
beyond the usual Lyman-
forest clouds. As argued in Cen & Ostriker (1999) and
Dave et al. (2001) this could be part of a "warm'' phase of the intergalactic
medium (IGM), with temperatures in the range
105 < T <107 K.
The latter conclusion was reached from numerical simulations. In this
article, we reconsider this problem in order to derive the properties of the
IGM by analytic means. In particular, we wish to investigate whether one can
understand this behaviour in a quantitative manner from robust,
model-independent, arguments.
First, in Sect. 2 we study the
phase-diagram of cosmological baryons. While the Lyman-
forest is described by a well-defined Equation of State the "warm'' IGM
component shows a broad scatter (e.g., Dave et al. 1999; Dave et al. 2001) since its
temperature depends through shock-heating on the neighbouring gravitational
potential which is a stochastic field. Nevertheless, we show that it is
constrained to lie in a well-defined domain in the
plane, and
determine its average location in this plane, which may be considered as the
"Equation of State'' of the "warm'' IGM. We also give the location in this
diagram of galaxies, groups and clusters.
Next, in Sect. 3 we use our results to compute the
redshift evolution of the fraction of matter enclosed within the different
phases. Then, in Sect. 4 we estimate the
two-point correlation function and the clumping factor of the "warm''
IGM. Finally, in Sect. 5 we check that the X-ray background emitted by the "warm'' component agrees with observations.
2 The phase-diagram of cosmological baryons
In this section we investigate the evolution of the physical properties of the baryons at moderate redshifts
,
as gravitational clustering builds up. To do so, we study the distribution of matter in a density-temperature plane and determine the
phase-diagram. Indeed, as we describe below baryons are located within specific regions in this
plane which reflect the various processes (cooling, radiative heating, shock-heating, gravitational clustering, ...) which take place in the universe. We display our results in Fig. 1. Note that we use the overdensity
rather than the density
in the figures, where
is the mean density of the universe.
In this article, all our numerical results are obtained with the following cosmological parameters. We consider a low-density flat universe with
and
.
The baryonic density parameter is
and the Hubble constant is H0=65 km s-1/Mpc. The
CDM power-spectrum of the linear density fluctuations is normalized by
.
We consider two processes for the possible heating of the gas. First, the gas temperature may be governed by a local shock-heating due to gravitational processes. Thus, as large-scale structures build up we expect tidal effects and virialization processes to heat the gas within and in the neighbourhood of collapsed halos and filaments up to the virial temperature of the underlying gravitational potential well. The latter scales as
,
which defines the characteristic scale R, hence we write:
 |
(1) |
where
is the mean molecular weight of the gas and
the proton mass.
Second, the gas temperature may be set by an external heating source, like a UV background radiation. Then, the temperature is no longer given by Eq. (1) since it is not fixed by gravitational energy. However, the pressure of the gas erases the baryonic density fluctuations over a scale R given by:
 |
(2) |
where
is the sound velocity and
.
This is simply the distance over which sound waves can travel in a Hubble time and build pressure equilibrium.
Therefore, the gas sees the dark matter density field smoothed over this scale R.
There is a characteristic density
at which shock heating starts to play a role, that is where Eq. (1) is to be used in place of Eq. (2). It is expected to correspond to a density contrast of a few units at least. The detailed discussion of the relation between the two regimes around
may be found in Sect. 2.4 below.
2.2 Exclusion constraints
First, we note that the distribution of matter in the
plane is strongly constrained by the existence of some exclusion regions where no particles can be found. The advantage of such conditions is that they are very robust in the sense they do not depend on the detailed history of the baryons. They apply to any hierarchical scenario of structure formation and provide the first guideline for the properties of baryons. We discuss these various constraints in the sections below. For sake of clarity, we turn within the diagrams shown in Fig. 1 in counter-clockwise order.
2.2.1 Cooling constraint
The dashed-line on the far-right with an inclined "V-shape'' is the familiar cooling constraint (see also Rees & Ostriker 1977; Silk 1977). It expresses the fact that high-density ionized hot gas (T > 104 K) cools very rapidly. Thus, baryons cannot remain for long in regions located to the right of this dashed-line. More precisely, we define the "effective'' cooling time
set by atomic physics as:
 |
(3) |
where
is the specific energy and
is the baryon number density. The cooling time
describes collisional excitation, collisional ionization, recombination, molecular hydrogen cooling, bremsstrahlung and Compton cooling or heating (e.g., Anninos et al. 1997). The heating time
is the radiative heating time due to a UV background radiation of the form
.
As usual, we normalize the UV flux by its value J21 at the HI ionization threshold (912 Å) in units of 10-21 erg s-1 Hz-1 cm-2 sr-1. The values we choose for J21 are given in Table 1 as a function of redshift. They are consistent with observations (e.g., Giallongo et al. 1996; Cooke et al. 1997; Vogel et al. 1995) and are the same as the ones used for modeling the evolution of Lyman-
clouds Valageas et al. (1999). Then, the cooling curve shown in Fig. 1 is given by the condition:
 |
(4) |
Note that in Eq. (4) we neglected shock-heating due to gravitational processes (e.g., virialization) and possible adiabatic cooling due to the expansion of the fluid elements. However, these processes are governed by the time
which measures the expansion rate of the universe and the time-scale over which new non-linear structures form.
Therefore, they do not modify the form of Eq. (4): they may only change the ratio between both sides of Eq. (4) (which we set equal to one) by a factor of order unity. Moreover, we have
as radiative heating does play a strong role here. Hence this constraint (4) is actually quite robust. Then, the regions located to the right of the cooling curve have a cooling time which is less than the Hubble time. As a consequence, no particle can remain in this part of the
plane for a Hubble time. This process is actually at the origin of galaxy formation as this gas cools and falls to the center of the gravitational well to form a disk and stars (see Sect. 2.7 and Valageas & Schaeffer 1999). This means that the matter which entered this region at some time should accumulate on the lower branch of the cooling curve (at
K) under the form of high-density cool gas which evolves on a Hubble time (but some of this gas also turns into stars). Note that this accumulation of baryons along the lower branch is clearly seen in the numerical simulations displayed in Fig. 11 in Dave et al. (1999) and in Fig. 9 in Springel & Hernquist (2001).
Table 1:
Redshift evolution of the background UV flux J21(z) used
in this article.
z=0 |
1 |
2 |
3 |
4 |
5 |
J21 = 0.05 |
0.5 |
0.8 |
0.4 |
0.2 |
0.1 |
2.2.2 High-density fluctuations
A second constraint is set by the properties of the dark matter density field itself. Indeed, at a given scale R beyond some threshold
large densities have a negligible probability and only involve a very small amount of matter. In the linear regime for instance, at large scales R, the probability distribution function (pdf)
of the density contrast
over a spherical cell of radius R is a Gaussian of variance
.
Here we note as usual
the rms linear density fluctuation and we assume Gaussian initial density fluctuations. In the non-linear regime, for highly non-Gaussian probabilities, the analogue of this upper threshold can also be calculated. It is given in Appendix A where we detail our model for the pdf
,
where
is the overdensity. The high-density cutoff defined by Eq. (A.8) obeys the asymptotic behaviour (A.10) in the quasi-linear and highly non-linear regimes. Note that this asymptotic behaviour actually is model-independent. Indeed, in the limit
we recover the Gaussian cutoff implied by the initial conditions (where
). On the other hand, in the limit
we recover the fact that most of the matter is enclosed within high-density halos of density contrast
which occupy a small fraction of the volume (
)
of the universe.
Therefore, the constraint given by the high-density cutoff of the pdf
in the
plane is quite robust and it applies to all hierarchical scenarios of structure formation.
We assume that the baryon density
scales as the dark matter density through
 |
(5) |
This is valid even in the cooling region discussed above in Sect. 2.2.1 if we consider also the cold baryons. We then can obtain an upper cutoff for the baryon overdensity
at any given scale R. It will turn out to be convenient to express this cutoff
in terms of a density-temperature relation, in order to draw its consequences for the phase-diagram
of the IGM. To do so we simply need the relation
.
Substituting Eq. (1) into the relation
derived in Eq. (A.8) we obtain a curve
.
The subscript "loc'' in the overdensity cutoff
refers to the fact that this is a "local'' heating process. It is due to the gravitational interaction with neighbouring structures. This high-density cutoff
corresponds to the dot-dashed curve shown in the diagrams in Fig. 1 which runs from
up to
and which crosses the cooling curve.
Substituting Eq. (2) into the relation
obtained from Eq. (A.8) we get a curve
.
This corresponds in Fig. 1 to the short branch which runs upward from the curve
at
.
Indeed, we note that for
the gas is located close to non-linear structures so that local shock-heating must be taken into account. However, for large densities we have
since Eq. (2) yields
.
In this case external heating plays no role. This is why we only plot the curve
up to the characteristic density contrast where
.
At higher densities the threshold
becomes irrelevant. In fact, we see in Fig. 1 that the curve
due to external heating plays no role since at moderate densities
it is repelled to quasi-linear scales (see Sect. 2.2.3).
Then, the region to the upper-right of the curves
and
in the
plane shown in Fig. 1 corresponds to rare high-density fluctuations which are located in the tail of the pdf
.
Therefore, there should be very few particles beyond these lines. Thus, this defines a second exclusion region.
Note that at z=0 this constraint yields an upper bound
K
keV. Of course, there exist some halos with a larger temperature: massive X-ray clusters. However, these are rare objects which only contain a small fraction of the baryonic matter content of the universe (typically
)
and they indeed correspond to the high-mass tail of the mass function. This is obviously consistent with the description of the baryonic matter which is worked out in this article, see Sect. 2.7.
2.2.3 Transition to non-linearity
A characteristic scale which enters the problem we investigate here is set by the transition to the non-linear regime. At a given epoch, we define this scale R0(z) by the relation:
 |
(6) |
Thus, scales larger than R0 are still within the linear regime. Then, there can be no shock-heating due to gravitational clustering on these large scales which have not turned non-linear yet. Using Eq. (1) we obtain a characteristic temperature
:
 |
(7) |
which describes the transition to non-linear scales (which lie at
). This is shown by the straight dashed-line (
)
plotted in Fig. 1. Besides, we note that known sources of external heating (e.g., the UV background radiation) cannot heat the IGM up to such high temperatures. Therefore, no fluid element can be located in the region to the upper-left of the curve
since no physical process which is active on these large scales can heat the gas to these high temperatures. This yields a third exclusion region.
2.2.4 Low-density fluctuations
Finally, a fourth constraint on the distribution of matter is given by the low-density cutoff of the pdf
of the dark matter density field. This is the analog of the high-density cutoff discussed in Sect. 2.2.2. The low-density cutoff
is derived in Appendix A, from Eq. (A.8). It obeys the asymptotic behaviours (A.11). In the quasi-linear regime we again recover the usual Gaussian cutoff
.
In the highly non-linear regime the small overdensity
expresses the formation of extreme underdensities on small scales. Contrary to the high-density cutoff it is somewhat model-dependent through the exponent
but this has no strong effect on the
phase-diagram (we typically have
). As in Sect. 2.2.2 we need a relation T(R) in order to derive a condition of the form
.
We again consider both cases of local and external heating, described by Eqs. (1) and (2). This yields the curves
and
shown by the two steep parallel dashed-lines in Fig. 1, at
.
The curve associated with external heating is the left one (i.e. lower densities or higher temperature), as can be seen from Eqs. (1) and (2). Thus, this defines a fourth exclusion region to the left of these curves.
2.3 Equation of State of the Lyman-
forest
Thus, so far we have obtained constraints on the distribution of matter in the
plane by drawing four exclusion regions. This already gives quite useful information about the properties of the IGM which are very robust. Now, we investigate a different point, seeking the location of the gas in the
plane. This amounts to deriving an equation of State for this component.
As shown in Hui & Gnedin (1997) the low-density photo-ionized IGM exhibits such an equation of State as the gas follows a specific relation
with a rather small scatter. This was derived in Hui & Gnedin (1997) from the Zel'dovich approximation (Zel'dovich 1970) which applies up to the moderately non-linear regime (
). Here we reconsider this problem and we show that this Equation of State is rather robust with respect to the past history of the gas and applies independently of the validity of the Zel'dovich approximation. First, we assume photo-ionization equilibrium (we restrict ourselves to
after reionization) and we only take into account hydrogen. Therefore, the ionization equilibrium reads:
 |
(8) |
where
is the photo-ionization rate and
is the recombination rate. They are given by:
 |
(9) |
and:
 |
(10) |
where we defined:
 |
(11) |
The temperature
we introduced in Eq. (11) is the characteristic temperature reached by the gas through the heating due to the UV background radiation flux. It is given by:
 |
(12) |
where
eV is the hydrogen ionization threshold. Note that this temperature
does not depend on the amplitude J21 of the UV background. Moreover, it is fixed by atomic physics, independently of cosmological parameters. Next, the temperature T of a given fluid element evolves as:
 |
(13) |
where
is the Lagrangian time derivative. The heating time
is given by:
where we take
since we approximate the gas as fully ionized hydrogen. The evolution Eq. (13) is the same as Eq. (3), except that we neglect cooling (which is justified here since we consider here moderate densities and temperatures) and gravitational shock-heating which is irrelevant. On the other hand, the term
represents the pressure work, which takes into account the expansion of the fluid element. Substituting Eq. (14) into Eq. (13) we obtain:
 |
(15) |
where we define:
 |
(16) |
and:
 |
(17) |
In Eq. (15) we used the fact that hydrogen is almost entirely ionized at the low redshifts which we consider here (z < 5, after reionization), as shown by the Gunn-Peterson test. Next, it is convenient to introduce the adiabat K defined as the exponential of the specific entropy s. More precisely, we define K and s by:
 |
(18) |
Then, Eq. (15) may be written:
 |
(19) |
Here we note that
is a small number. Therefore, the specific entropy of the fluid shows a weak dependence on the evolution of its density
.
Hence we approximate the solution of Eq. (19) by:
 |
(20) |
where
is the density at the time t we consider. In terms of the temperature-density relation, this yields:
 |
(21) |
where we used Eq. (18). Thus, the expression (21) gives the first-order term for the expansion of
in terms of the small parameter
.
From Eq. (16) we note that at
we have
and
hence the temperature of the photo-ionized IGM will be of order
K.
Note that Eq. (21) is independent of the normalization of the UV flux
.
Indeed, the efficiency of radiative heating is proportional to
but the density of neutral hydrogen scales as
(at ionization equilibrium for almost fully ionized gas) so that
cancels out. Therefore, the result (21) is quite robust since it does not depend on the value of the UV flux. In particular, Eq. (21) still holds even if the UV background is inhomogeneous: we only need to assume local ionization equilibrium. This explains why there is only a very small scatter around the equation of State (21) since the actual physical conditions within these small clouds are almost independent of the actual history of each fluid element (i.e. the evolution of its density and local UV flux). This indeed agrees with the results of numerical simulations (e.g., Dave et al. 1999).
We display in Fig. 1 the equation of State (21) as the solid line which runs through
K at
.
Beyond the overdensity
defined in Eq. (22) below, we plot this line as a dashed-line, until it enters the cooling region described in Sect. 2.2.1. Indeed, as we explained in Sect. 2.2.2 at large densities the "virial temperature''
becomes larger than the temperature
due to some external energy source (here photo-ionization heating by the UV background radiation). This means that for these regions, which have already reached the non-linear regime as
(see Eq. (22)), shock-heating due to the gravitational dynamics can no longer be neglected and it actually becomes dominant. Therefore, for high densities with
the gas should no longer fall onto the curve (21). Nevertheless, since shock-heating can only increase the temperature of the gas the relation (21) now provides a lower bound to the temperature T. Hence, the region below the curve (21) is excluded in the
diagram. This holds until we enter the cooling region discussed in Sect. 2.2.1. On the low-density side, as described in Sect. 2.2.4 we are constrained by the low-density cutoff
.
Therefore, we predict that we should have two phases for the IGM. A first "cool'' phase is described by the Equation of State (21) with intermediate densities
.
It is photo-ionized gas heated up to
K by the background UV flux. This corresponds to the moderate density fluctuations which form the Lyman-
forest. A second "warm'' phase is made of higher-density regions which have already experienced some shock-heating due to the building of gravitational structures but which have not entered the cooling region yet. These particles should be located in the
plane above the curve (21) and within the constraints described in the previous sections.
![\begin{figure}
\par\includegraphics[width=8.5cm,height=6.3cm,clip]{2143.f1}\hspa...
...e*{6mm}
\includegraphics[width=8.5cm,height=6.3cm,clip]{2143.f4}
\end{figure}](/articles/aa/full/2002/24/aa2143/Timg132.gif) |
Figure 1:
The phase-diagram of the IGM from z=3 down to z=0. The straight solid line
shows the Equation of State of the "cool'' IGM (Lyman-
forest). The curved solid line
shows, as a mean trend, the "Equation of State'' of the "warm'' IGM which is shock-heated through the building of non-linear gravitational structures. The dashed curves draw exclusion regions around the allowed domain for this "warm'' IGM. In counter-clockwise order, starting from the right side, they correspond to 1) a fast-cooling region where the gas cannot remain over a Hubble time, 2) a high-density and high-temperature domain within the exponential tail of the pdf
which only contains very rare massive halos, 3) large linear scales where gravitational shock-heating has not appeared yet, 4) a low-density region within the tail of the pdf
which is associated with very rare voids and 5) a lower-bound for the temperature
set by radiative heating from the UV background. The vertical dotted line is the density threshold
of just-virialized halos. The points show the results of numerical simulations from Dave et al. (1999) (Fig. 11). |
Open with DEXTER |
2.4 Warm IGM
Thus, the previous discussion shows that a "warm'' component of the IGM should no longer follow the equation of state (21). It however remains constrained by the exclusion regions obtained in Sect. 2.2. This leaves a closed allowed region in the
plane. Since the temperature of the gas should be governed by the shock-heating due to the gravitational dynamics we can expect a rather large scatter. Indeed, the history of these particles depends on the properties of their neighbouring dark matter density field. This introduces a stochastic component which gives rise to a wide variety of possible histories so that there is no longer a unique temperature-density relation. This agrees well with Fig. 11 in Dave et al. (1999) which shows indeed a broad cloud of points in the
plane within the allowed region we obtained in the previous sections. Note also that at z=0 some of the gas which has been shock-heated to T > 104 K actually lies at low density contrasts
in the numerical simulation. This is consistent with our results shown in Fig. 1. This simply means that on small non-linear scales even regions with
have experienced shell-crossing and gravitational shock-heating (this also implies that they are not just described by a collapse dynamics).
However, the distribution of matter obtained in Dave et al. (1999) from numerical simulations does not fill entirely this allowed region (note indeed that this is not implied by our previous considerations) and it is not uniform. Hence it would be convenient to derive a curve in the
plane which would describe the "mean'' behaviour of this warm phase of the IGM. To this purpose, we first compute the gas temperature associated with a given halo of mass M which forms at redshift z. As noticed above, this does not apply to low-density regions with
but this model should provide a useful estimate for the higher density regions which actually contain most of the mass. Gravitational heating is effective in the collapsing phase of the halo which may be set at
.
One expects
to be of the order of a few units. We set it to
 |
(22) |
which is close to the value given by the spherical collapse dynamics at turn-around. We consider this not to be an accidental coincidence. Indeed, the kinetic energy of the underlying overdense large-scale structure is not thermalized during its expansion phase and it only provides an energy source for the heating of the gas within the collapse stage. Even for an expanding cloud which is not perfectly spherical, such a phase of maximum expansion is reached, where the kinetic energy in the collective motion is small. Afterwards, the cloud enters a collapse phase where shocks develop. As is usual for shocks, the latter transform collective large-scale kinetic energy into small-scale thermal energy. The collective kinetic energy per particle v2 in the collapse phase (i.e. the squared velocity of the inflow) is of order
where the constant term v02 expresses the fact that at turn-around v2=0. Then, this kinetic energy is locally transformed into heat in the shocked regions. The temperature reached there
is of the order of this kinetic energy since it is set by the inflow velocity v through the matching conditions at the interface between the inner warm shocked region and the outer collapsing region. We thus obtain a relation between the small-scale temperature (which applies to the small shocked regions) and the large-scale kinetic energy per particle. This yields (for a large-scale structure of mass M):
![\begin{displaymath}%
T_{\rm gh}\sim T_{\rm c}\left[ \left(\frac{\rho}{\rho_{\rm c}}\right)^{1/3} -1 \right] ,
\end{displaymath}](/articles/aa/full/2002/24/aa2143/img139.gif) |
(23) |
with
 |
(24) |
Of course, the characteristic temperature
depends on the mass M of the collapsing cloud. The subscript "gh'' refers to the fact that the temperature
is governed by "gravitational heating''. This provides an estimate
of the location of the warm spots in the
plane, for fluid elements parameterized by the mass M of the underlying larger-scale structure which is collapsing. Here we assumed that the density within the shocked regions scales as the overall density of the collapsing structure (note that a strong shock only increases the density by a factor 4 at most).
Now, seeking an average location of these warm spots, we estimate the typical mass M0(z) which collapses at each redshift z by:
 |
(25) |
This also yields the scale R0(z) defined in Eq. (6), with
.
Among all the
relations (23) which have been seen to depend on the halo mass M, we can then find a "mean'' locus by taking a characteristic temperature
,
Eq. (24), calculated with the mass M0.
The mass M0 and the scale R0 characterize the non-linear structures which just turned non-linear. For instance, at z=0 we have
Mpc. At virialization when the density contrast reaches the threshold
this yields a radius
Mpc as for large clusters. On the other hand, the shocked regions, at temperature
,
correspond to much smaller scales. Their size R can be estimated by using approximate hydrostatic equilibrium within the filaments along the transverse direction
.
As a consequence, the Eq. (1) still holds, up to factors of order unity, if R is taken to be the thickness of the filament. This is a natural result since the inflow and the shock actually are at the origin of the filament. Also, the thickness R of the filaments turns out, as a rule, to be much smaller than R0. For instance, at z=0 we find it is of order of a few hundred kpc, see Sect. 4 below. Note that the same processes are at work for the dark matter density field, except that shocked regions now correspond to areas where shell-crossing governs the dynamics and builds a large velocity dispersion. Therefore, we can again assume that the gas follows the dark matter density field and the gravitational potential is dominated by the dark matter. Our findings agree with the results of numerical simulations which show that the shocked regions are associated with the filaments which appear in the dark matter density field. Thus, as shown by the simulation map in Fig. 3 in Dave et al. (2001), the "warm'' IGM is an intricate network of filaments which extends over a few Mpc (the scale R0) but the thickness of the filaments (the scale R) is much smaller than this global scale.
In Sect. 4 we will see in addition that the scale R of the warm regions shows up as a small-scale cutoff for the "warm'' IGM two-point correlation function. We obtain a length of a few hundred kpc at z=0 which is much smaller than the few Mpc which characterize the underlying global structure which is collapsing. This also agrees with the numerical simulations of Dave et al. (2001).
Therefore, the "warm'' IGM component is described by the curve (23) in the
plane, together with Eq. (1) which yields the local size R of the clouds (or the thickness of the filaments). Our model relies on three major points: 1) there is a critical overdensity associated to the "warm'' IGM which reflects the turn-around of patches just going non-linear, 2) shock-heating locally transforms into heat the collective kinetic energy and 3) the scale of the shocked regions is given by local approximate hydrostatic equilibrium at the latter temperature. This yields 1) the locus in parameter space
where the "warm'' IGM appears, 2) the size of the non-linear network of shocked regions (filament network) which is the scale R0 just turning non-linear and 3) the thickness R of the filaments which appears to be much smaller than the previous scale.
We also see that at the beginning of the cloud collapse the shocked regions are predicted to have low temperatures correlated with a rather small spatial extension. As the overall collapse proceeds their temperature and their size increase and the highest temperatures are reached when the cloud is being virialized with scales of the order of the radius of the halo. This is more or less what could be expected to occur. Thus, we model the mean trend of the "warm'' IGM by Eq. (23) which appears as a curved solid line labeled "
'' in Fig. 1. It starts from the end-point of the "Cool'' IGM (at
)
and it goes towards larger densities and temperatures until it reaches the cooling region (at
)
or the density threshold
(at low z). Indeed, beyond this point we consider that we have cooled objects (galaxies) or hot virialized halos (clusters) which are not part of the "warm'' IGM.
As seen from Fig. 1, at low z gravitational heating rapidly dominates over the UV heating, for densities close to
and temperatures around
.
On the other hand, the "warm'' IGM clouds are not the same as the UV heated objects which form the "cool'' component. Note that the curve (23) used with the mass (25) only reflects the average trend. Indeed, the stochastic character of the dark matter density field leads to a broad variety of masses which are just collapsing, and hence of
trajectories. Moreover, the local properties of the shocks also provide for some additional scatter. This induces a (rather large) dispersion of the points of the "warm'' IGM in the
plane, which agrees with numerical simulations (e.g., Dave et al. 2001). However the "warm'' IGM should remain enclosed within the allowed region defined by the constraints discussed in the previous sections. Indeed, we must point out that the validity of Eq. (1) ensures that the constraints obtained in Sects. 2.2.2-2.2.4 still hold.
2.5 The different phases of the IGM
In the previous sections we have shown that one can distinguish two components in the IGM and we have determined their location in the
phase-diagram. Our results are displayed in Fig. 1.
Firstly, there is a "cool'' IGM phase (
K) which corresponds to the Lyman-
forest. These are moderate density fluctuations (
)
of photo-ionized gas. They are described by the Equation of State (21) which arises from the heating of the gas by the UV background and the cooling due to the expansion (i.e. pressure work). Thus, this component lies on a well-defined curve in the
plane. This curve is bounded towards low densities by the cutoff of the pdf
,
which expresses the fact that the dark matter density field arising from Gaussian initial conditions exhibits a finite range of densities which occur with a significant probability. Note that this lower bound indeed agrees with the points obtained from numerical simulations shown in Fig. 1. On the other hand, the high-density bound
is due to gravitational shock-heating which becomes the dominant energy source for dense regions.
Secondly, there is a "warm'' IGM phase (
K) which describes the gas heated by shocks arising from the gravitational energy of just collapsing objects. Because of the stochastic character of this energy source there is a broad scatter for this component around the "Equation of State'' we have derived. This gas is restricted to a specific allowed region in the
phase-diagram. This expresses cooling and heating constraints as well as the properties of the underlying dark matter density field. We nevertheless obtained a curve, Eq. (23), which follows the mean trend of this "warm'' phase. Its low density bound is set by the transition near
with the "cool'' IGM phase dominated by radiative heating from the UV background. The high-density bound is given by the intersection with the cooling curve (where bremsstrahlung cooling becomes dominant) or the density threshold
(beyond this point we have groups or clusters of galaxies). Note that these results are consistent with the calculations of (Nath & Silk 2001) based on the Zel'dovich approximation.
We can note that our results shown in Fig. 1 agree reasonably
well with the outcome of numerical simulations as displayed in Fig. 11 in
Dave et al. (1999) (also shown by the points in our Fig. 1) and Fig. 9
in Springel & Hernquist (2001). There is a small offset at low redshift for the
normalization of the equation of state (21) of the "cool'' IGM and
for the cooling region defined in Eq. (4). Note that the latter
could be remedied by adjusting the ratio
which we simply set
equal to unity in Eq. (4). Similarly, we could obtain a better fit
to the numerical results for the "warm'' IGM by tuning the r.h.s. in
Eq. (25) which defines the mass M0. However, our goal is not to
get the best fit to a specific numerical simulation (which would be of
little value) but to explain the physics of the IGM. Moreover,
Fig. 1 shows that the simple procedure detailed in the previous
sections already provides a good qualitative and quantitative
description which should be sufficient for most purposes. Besides, as
explained above it should be quite robust. In particular, it could be readily used with any cosmological parameters.
2.6 Redshift evolution
Finally, the four diagrams displayed in Fig. 1 show the evolution with redshift of the properties of the IGM, from z=3 down to z=0. We can see that the pattern does not evolve much qualitatively although the curves exhibit a quantitative shift with z. The characteristic temperature of the "cool'' IGM decreases slightly with time, in agreement with Eq. (21) since the density declines faster than
.
Indeed, the cooling rate due to the expansion of the universe scales as
while the recombination rate (which yields the density of neutral hydrogen involved in radiative heating) scales as
.
On the other hand, the cooling constraint defined by
in Sect. 2.2.1 does not evolve much with time. In particular, its lower branch at
K is set by the atomic physics of hydrogen ionization. Next, the characteristic temperature of the curve
associated with shock-heating due to gravitational clustering grows with time. This expresses the fact that the virial temperature associated with larger scales which turn non-linear later is higher. For instance, the velocity dispersion associated with galaxies is of order
km s-1 while for clusters it is
km s-1. Following this evolution, the mean curve
which describes the "warm'' IGM enters the cooling region at high z while at low z it first crosses the high-density threshold
.
This actually expresses the fact that at high z the collapsed halos built by gravitational clustering form galaxies since the gas undergoes a very efficient cooling, while at low z typical just-virialized halos are clusters which remain hot over a Hubble time and are still strong X-ray emitters (e.g., Valageas & Schaeffer 2000).
2.7 The IGM versus virialized halos
![\begin{figure}
\par\includegraphics[width=8.5cm,height=6.3cm,clip]{2143.f5}\hspace*{8mm}
\includegraphics[width=8.5cm,height=6.3cm,clip]{2143.f6}\end{figure}](/articles/aa/full/2002/24/aa2143/Timg161.gif) |
Figure 2:
The phase-diagram of cosmological baryons at z=3 and z=0. As in Fig. 1 the curves at
are the equations of state of the "cool'' and "warm'' phases of the IGM. The dashed curve with a vertical part at
and a branch towards higher
and T corresponds to galaxies. The two upper parts of the vertical line at
are groups (solid line at z=0) and rare clusters (dotted line). The upper right dot-dashed line corresponds to the cool cores of groups and clusters where cooling has had time to develop. |
Open with DEXTER |
As explained in the previous sections, we restrict the "mean equation of state'' of the "warm'' IGM to density contrasts below the threshold
.
Indeed, we consider that larger densities correspond to galaxies or clusters, that is virialized halos which are not part of the diffuse IGM. We display in Fig. 2 the location in the
phase-diagram of these collapsed objects and of the IGM at redshifts z=3 and z=0.
Let us first consider the right panel, obtained for z=0. The dashed line at
is the "cool'' IGM while the solid curve at
is the "warm'' IGM. The equations of state of these two components were derived in the previous sections and these curves are identical to those shown in Fig. 1. For clarity we do not plot in Fig. 2 the boundaries discussed in Fig. 1 which constrain the large scatter of the "warm'' IGM phase.
Next, the vertical line at
shows the overall density contrast associated with just-virialized objects. We divide this line into three parts.
The first is in the low temperatures
K region. These objects with
are located within the cooling region shown in Fig. 1. Hence their cooling time
is smaller than the Hubble time
,
see Eq. (4). Therefore, within these halos the gas undergoes a very efficient cooling which leads to the formation of stars. As a consequence, these objects are small galaxies which are just being formed. This part of the vertical line
which is enclosed within the cooling region is shown by the lower vertical dashed line in Fig. 2.
The remaining high temperature region may still be subdivided. The halos at
which have a higher temperature do not cool over a Hubble time since they obey
(except in their center). Hence they correspond to groups or clusters which still contain hot gas which can be observed through its X-ray emission. One may divide this part into two components: groups versus clusters. We identify groups with the low temperature halos located below the cutoff
shown in Fig. 1. Therefore, they correspond to the typical just-virialized halos which form at z=0. On the other hand, we identify clusters with the high temperature halos located above the cutoff
.
Hence they are rare massive objects which probe the high-density tail of the pdf
.
In other words, they correspond to the high-mass tail of the mass function. Thus, our distinction between groups and clusters is only based on the abundance of these objects (whether they correspond to rare or typical density fluctuations). However, we shall explain in a future paper (Valageas et al. 2002) that this subdivision also marks the boundary between objects which, depending on the depth of their gravitational potential, are affected or not by a preheating of the gas through the gravitational processes which yield the IGM. Groups (resp. clusters) are shown in Fig. 2 by the vertical solid (resp. dotted) line. Thus, these three parts of the vertical line at
describe the location of small galaxies, groups and clusters.
For completeness, we must point out that the lower part of the line at
with
K does not give all galaxies. Indeed, it is clear that there are some galaxies with a higher temperature. As discussed in Valageas & Schaeffer (1999) these objects must still satisfy a somewhat stronger cooling constraint:
 |
(26) |
where
is the redshift of formation of the object. Halos with
simply have
(i.e. they are just being virialized) while higher density contrasts correspond to older objects which collapsed at earlier times. The constraint (26) ensures that cooling was very efficient when these halos formed so that their gas was able to condense fast enough in order to be left unperturbed by the next generation of forming structures. This allows these objects to retain their individuality, stars and possibly a gaseous disk, even if they become embedded within larger structures (e.g., filaments or clusters) at later times. The condition (26) at
is shown by the dashed line in Fig. 2, where we used the approximation
(i.e.
scales as the dynamical time of the halo). Thus, we are led to consider that the properties of these galactic halos (mass and size) do not evolve much after they form and reach the virialization density contrast
.
This method yields predictions for the properties of galaxies (luminosity function, masses, metallicity, ...) which agree with observations. Therefore, galaxies are defined by the dashed line which shows a first vertical part at
and
K (small galaxies which are being formed) and a second part at higher densities and temperatures (older and more massive galaxies). We can check that we obtain at T = 106 K a radius
kpc which is indeed the typical size of present galaxies. By contrast, a halo with the same virial temperature at
would give a size
kpc (and a mass
)
which is uncomfortably large for a galaxy. This actually corresponds to a small group of galaxies.
Finally, we note that the upper part of the cooling curve (4) shown in Fig. 1 at
describes the cool cores of groups and clusters. Indeed, it corresponds to high-density regions which have had time to cool today. This yields the core radius of present groups and clusters where cooling has just had time to come into play and possibly induce cooling flows. For instance, at T = 107 K we obtain a radius
kpc, which agrees with observations, while the virial radius is
Mpc. We display this curve as the dot-dashed line in Fig. 2. Note that this assumes that baryons follow the dark matter. However, at low z this is not necessarily the case since some "preheating'' may modify the physics of the gas in low-temperature clusters (i.e. groups), see for instance Valageas & Silk (1999) and Valageas & Schaeffer (2000). This is not important for our purpose here which is mainly to describe the physics of the IGM. A specific study devoted to this problem and the entropy of the gas in the light of the results described in this paper is presented in Valageas et al. (2002).
Note that for these collapsed objects (galaxies, groups, clusters and cool cores) the density contrast
and the temperature T shown in Fig. 2 only refer to the mean density contrast and virial temperature of the halo over the relevant radius R (which may be different from the virial radius for galaxies and cool cores). At the center of the halo the gas density is larger. Moreover, for galaxies and cool cores the gas temperature can be significantly smaller since cooling is very efficient (and we should take into account feedback from supernovae). Therefore, contrary to the IGM where the location in the
phase-diagram directly gave the properties of the gas, here the
plane only shows the overall properties of the dark matter halos associated with each object. Finally, we must point out that these classes of objects are not exclusive of each other. Indeed, cooled cores are obviously embedded within groups and clusters while galaxies can be found within filaments, groups and clusters. Therefore, some of the mass associated with galaxies or cool cores is counted within the matter attached to filaments, groups or clusters (this can be handled using the methods of Valageas & Schaeffer 1997, see Valageas & Schaeffer 1999 for galaxies and Valageas & Schaeffer 2000 for clusters).
Thus, as explained above the
diagram shown in Fig. 2 describes the physics of cosmological baryons and their distribution between different phases.
The redshift evolution up to z=3 can be easily derived from Fig. 1 and for illustration we show in the left panel in Fig. 2 our result at z=3. The location of the various curves evolves as explained in Sect. 2.6 but we can also note some qualitative changes. Firstly, with our definitions we see that there are no more groups. This means that the typical objects which collapse at z=3 form galaxies as they exhibit efficient cooling. To be more precise, one can still find clusters along
,
at high temperatures where cooling is not very efficient. However, they are located within the far tail of the mass function and they correspond to extremely rare events. Secondly, we note that the curve which describes the mean "equation of state'' of the "warm'' IGM stops below
.
Indeed, as seen in Fig. 1 its high-density end-point is now given by the intersection with the cooling region. This is merely another consequence of the fact that we typically form galaxies and not groups (this feature coincides with the "disappearance'' of groups). At an even higher redshift
the "warm'' IGM component almost disappears as the equation of state of the "cool'' IGM extends up to the cooling region. However there still remains a "warm'' phase because the allowed region in the
plane has not vanished
. At these high redshifts most of the mass is within a roughly uniform "cool'' phase. We still have some rare collapsed halos which correspond to galaxies and some non-linear structures where gravitationally induced shocks heat the gas but these latter regions are severely restricted by the high efficiency of cooling processes.
3 Distribution of matter
3.1 Mass fraction within different phases
Thus, in Sect. 2 we have
discussed how the properties of cosmological baryons could be seen through a
phase-diagram. In particular, we distinguished a "cool'' and a
"warm'' IGM phase. In this section, we derive the fraction of matter
enclosed within these various components.
Let us first discuss the "cool'' IGM phase. As described in Sect. 2.3
this corresponds to the Lyman-
forest, that is moderate density
fluctuations governed by the ionization and the heating due to the UV
background. Since its temperature is not zero this gas probes the dark matter
density field over the scale R defined in Eq. (2), which describes
the length-scale over which pressure can homogenize the baryonic matter
distribution. Note that this length depends on the temperature since
.
Then, we wish to express the fraction of matter within the
"cool'' IGM in terms of the pdf
over the scales R associated
with these clouds. To do so, we first note that these scales are within the
non-linear regime as shown by the curve
in Fig. 1
which marks the transition to the linear regime. Then, as discussed in
Valageas & Schaeffer (1997), in the highly non-linear regime the pdf
shows the
scaling:
 |
(27) |
which defines the scaling function h(x). This holds for densities
larger than the low-density cutoffs
discussed in
Sect. 2.2.4. Note that if the stable-clustering
ansatz is valid this function h(x) does not depend on redshift, as long as
one remains in the highly non-linear regime. Then, the fraction of matter
enclosed within spherical cells of radius R with a density between the
thresholds
and
is:
 |
(28) |
The advantage of the scaling variable x is that the last term in Eq. (28) still provides an estimate of the fraction of matter
associated with a curve
which is not necessarily at constant radius, as discussed in Valageas & Schaeffer (1997). This prediction is compared with numerical simulations in Valageas et al. (2000) where we investigate the two cases of a constant density and a constant radius. Therefore, in order to estimate the fraction of matter associated with the Lyman-
forest we use the variable x and we write:
 |
(29) |
with
 |
(30) |
where
is the value of the pdf
at scale R and redshift z and the integral is taken along the relevant curve
.
The pdf
is obtained from Eq. (A.2) where the parameter
which enters Eq. (A.5) is computed at the relevant scale R. Thus, the expression (29) would be exact within the framework described in Appendix A.1 if the Lyman-
forest were defined by a constant scale R. The use of the variable x in Eq. (29) is only meant to handle the fact that R actually varies with
.
In particular, we must point out that our procedure (29) with the definition (30) does not assume that the stable-clustering ansatz is valid nor that we consider highly non-linear scales. Indeed, the approximation (29) merely expresses the mass fraction associated with certain clouds of size R and density
in terms of the pdf
measured at this point
.
As described in Appendix A.1 our approximation for the pdf applies both to the linear and non-linear regimes. Moreover, this formalism does not imply that these clouds are virialized objects which have reached an equilibrium state. The "extended'' function h(x;R,z) defined in Eq. (30) depends on scale and redshift. It is merely another way to write the pdf
(which goes over to the one relevant for the highly non-linear regime in this limit and is thus used by analogy). An important point of Eq. (29) is that it yields the mass fraction associated with various objects from the pdf
of the non-linear density field. This is a strong advantage since, as described in Valageas & Schaeffer (1997), one can show that in the highly non-linear regime, assuming the stable-clustering ansatz is valid, the cloud-in-cloud problem can be handled in a satisfactory way. We can expect this property to extend to the regime we consider here. Note also that Eq. (29) allows us to count objects which are not necessarily defined by a constant density contrast.
![\begin{figure}
\par\includegraphics[width=8.8cm,height=6.5cm,clip]{2143.f7}\end{figure}](/articles/aa/full/2002/24/aa2143/Timg196.gif) |
Figure 3:
The distribution of baryonic matter. The dashed line which increases with redshift is the fraction of matter within the "cool'' IGM phase, as computed from Eqs. (29) and (21). The dot-dashed line gives the mass fraction associated with the "warm'' IGM component. The solid line is the fraction of matter within collapsed objects. The filled symbols show the results of the numerical simulations from Dave et al. (2001) (panel D2 in Fig. 1) for the "cool'' IGM (squares), the "warm'' IGM (triangles) and condensed gas (circles). The empty squares show the results of the simulation presented in Dave et al. (1999) for the "cool'' IGM phase (Fig. 12). |
Open with DEXTER |
Thus, the integration of Eq. (29) along the Equation of State (21) yields the fraction of matter
within the "cool'' IGM phase. Next, we could try to use the same method for the "warm'' IGM, that is to integrate Eq. (29) along the mean "equation of state'' (23). However, this procedure is not satisfactory since the filaments and large non-linear structures also contain galaxies and groups which should be subtracted. As a consequence we rather compute the mass fraction associated with virialized halos at
.
That is, we integrate Eq. (29) along the vertical line at
shown in Fig. 2, starting from the lower point which is given by the lower intersection with the cooling curve. This yields the mass fraction
associated with virialized halos (galaxies, groups and clusters). Then, we simply define the mass fraction
within the "warm'' IGM by:
 |
(31) |
so that the sum is unity. This provides the distribution of matter within these three components: the "cool'' IGM, the "warm'' IGM and collapsed halos.
We show our results in Fig. 3. The Lyman-
forest corresponds to the dashed line which grows at larger redshift. Indeed, at higher z gravitational clustering was less advanced so that a smaller fraction of matter had been shock-heated to high temperatures or embedded within collapsed objects through the building of large scale structures. Note that the fractions of matter within collapsed halos and the "warm'' IGM are of the same order. Indeed, both components are related to the formation of non-linear gravitational structures. As noticed in Sect. 2.7 we obtain
at z > 4. This actually means that the mass fraction within the "warm'' IGM is very small and within the inaccuracy of our computation.
We can check that our results agree with the outcome of the numerical simulations described in Dave et al. (1999) and Dave et al. (2001). The mass fractions we obtain for the "warm'' IGM and virialized objects are also in reasonable agreement with simulations although it is difficult to make a detailed comparison. Indeed, our separation between different components is not exactly the same as in the simulations
and the predictions of the various simulations exhibit a significant scatter as shown by the comparison between the filled and empty squares in Fig. 3 (see also the various panels in Fig. 1 in Dave et al. 2001).
3.2 Temperature distribution of the warm phase
![\begin{figure}
\includegraphics[width=8.8cm,height=6.5cm,clip]{2143.f8}\end{figure}](/articles/aa/full/2002/24/aa2143/Timg202.gif) |
Figure 4:
The temperature distribution
of the "warm'' IGM component at redshifts z=2 and z=0 (normalized to unity). The mean temperature grows with time as larger scale structures form. |
Open with DEXTER |
We have shown in Fig. 3 the fraction of matter within the "warm'' IGM. It is also of interest to consider more precisely the temperature distribution within this "warm'' IGM. As we explained in Sect. 3.1, the mass function which describes the "warm'' IGM cannot be obtained by directly integrating Eq. (29) along the mean "equation of state'' (23) because the regions which are about to collapse also contain smaller virialized halos (galaxies, groups). Thus, we first need to subtract these substructures which yields Eq. (31) for the mass fraction. This subtraction is more difficult when we consider the detailed temperature distribution within the "warm'' IGM but we shall assume that it can be described by a simple overall renormalization. This would be exact if the fraction of matter enclosed within such virialized objects would be the same for all "warm'' IGM regions (i.e. filaments). In fact, we can expect some correlation between these substructures and the larger-scale filaments but this procedure can be seen as a zeroth-order approximation. As a consequence, following Eq. (28) we write the fraction of matter
within the range
as:
 |
(32) |
which we compute along the curve (23) which describes the mean "equation of state'' of the "warm'' IGM. In Eq. (32) and in the following we do not write explicitly the time and scale dependence of h(x;R,z). The normalization constant in Eq. (32) is set by the constraint
,
which is the total mass fraction within the "warm'' IGM derived in Sect. 3.1.
We display our result normalized to unity in Fig. 4. Of course, we recover the fact that the mean temperature increases with time as larger non-linear structures build up. This could also be seen from Fig. 1. In particular, the mean trend curve of the "warm'' IGM in Fig. 1 at z=0 is very steep because the characteristic temperature of the non-linear structures which collapse at that time is of order
K which is much larger than the temperature reached by radiative heating from the UV background
K. Therefore, gravitational shock-heating easily heats the gas up to
K. However, note that there is still a significant fraction of the matter at lower temperatures T < 106 K, in agreement with Dave et al. (2001). Although the procedure (32) is only a simple approximation, the comparison of Fig. 4 with Fig. 5 in Dave et al. (2001) shows that it captures the main trend. Of course, as explained in Sect. 2.4 the "warm'' IGM shows a broad scatter in the
plane so that the actual boundaries of the temperature distribution are not as sharp as in Fig. 4. On the other hand, note that in our description we also have some matter at T>107 K but we associate this gas with massive X-ray clusters. In any case, our result displayed in Fig. 4 appears to provide a reasonable description of the "warm'' IGM properties.
4 Warm IGM two-point correlation
Another check on our description of the "warm'' IGM component is provided by its two-point correlation function
which measures its spatial clustering. Hence in this section we evaluate
,
which is defined by:
 |
(33) |
where
is the density of the "warm'' IGM component at point
.
In Eq. (33) the averages
are over the realizations of the density fields or over space (assuming ergodicity). First, following the procedure described in Sect. 3.2 we can write the mean "warm'' IGM density as:
 |
(34) |
where M, R and V are the mass, the radius and the volume associated with the point
in the phase-diagram. The normalization constant
is defined by:
 |
(35) |
Note that we always have
since the "warm'' IGM component does not contain all the baryonic matter. Next, the second-order moment is given by:
 |
(36) |
where
is the mean "warm'' IGM density at distance
from a point embedded within the "warm'' IGM component with the parameter x1. Then, as in Valageas et al. (2001a), we write
as the sum over two contributions. First, there is a probability F(r/R1) that the second point at
belongs to the same region of thickness R1 of the "warm'' IGM. Assuming spherical regions we have:
 |
(37) |
This gives the first contribution
as
.
Second, we write the contribution
to
of other regions of the "warm'' IGM component with parameter x2 as:
![$\displaystyle %
\rho_{\rm w}^{(2)}({\vec r},x_1) = \left[ 1 - F(r/R_1) \right] ...
... \frac{{\rm d} x_2}{x_2} x_2^2 h(x_2) \left[ 1 + b(x_1) b(x_2) \xi(r) \right] .$](/articles/aa/full/2002/24/aa2143/img225.gif) |
|
|
(38) |
Here
is the two-point correlation of the underlying dark matter density field while
b(x1) b(x2) is the bias associated with the regions x1 and x2. This bias is computed from the generating function introduced in Eq. (A.2) as described in Bernardeau & Schaeffer (1992) (see also Valageas et al. 2001b for a comparison with observations). Thus, we write the second-order moment as:
![$\displaystyle %
\langle\rho_{\rm w}({\vec r}_1) \rho_{\rm w}({\vec r}_1+{\vec r...
...x_2} x_2^2 h(x_2) \left[ 1 {+} b(x_1) b(x_2) \xi(r) \right] \biggl \rbrace\cdot$](/articles/aa/full/2002/24/aa2143/img226.gif) |
|
|
(39) |
This yields the "warm'' IGM two-point correlation
through Eq. (33).
![\begin{figure}
\par\includegraphics[width=8.8cm,height=6.5cm,clip]{2143.f9}\end{figure}](/articles/aa/full/2002/24/aa2143/Timg227.gif) |
Figure 5:
The "warm'' IGM two-point correlation
at redshift z=0 (dashed line). The solid line shows the two-point correlation
of the dark matter density field. |
Open with DEXTER |
We show our result at z=0 in Fig. 5. On large scales the "warm'' IGM correlation closely follows the dark matter correlation because i) we assumed that baryonic matter follows the dark matter density field smoothed over the scale R which describes the processes at work (R is set by the location on the
plane) and ii) the "warm'' IGM component is mainly made of typical density fluctuations. This is seen by the fact that, by definition, most of its matter is within the allowed region in the
phase-diagram, below the high-density cutoff
.
This implies that the biases which enter the expression (39) are close to unity since we do not select rare high densities or voids by looking at the "warm'' IGM phase. In fact, as seen from Fig. 5 we obtain a mean bias which is slightly smaller than unity. However, this does not take into account the broad scatter in the
plane of the "warm'' IGM component.
On the other hand, at small scales
kpc the two-point correlation function
flattens and it reaches a finite limit at r=0. This expresses the fact that the non-zero temperature of the gas reached by radiative heating and shock-heating (T > 104 K) gives rise to pressure effects which homogenize the baryonic distribution over some scale R. A lower limit to this length R is set by the scale
associated with the critical point
in the phase-diagram. Indeed, as shown by the curved solid line in Fig. 1, at redshift z=0 the length scales associated with the "warm'' IGM run from
up to scales of order R0 which are turning non-linear. This sets the knee of the two-point correlation
at
kpc at z=0. As noticed in Sect. 2.4, the result displayed in Fig. 5 shows that the characteristic thickness of the filaments which build the "warm'' IGM is of order 300 kpc at z=0. We find this, at a given redshift, to be indeed much smaller than the scale which is just turning non-linear, in agreement with numerical simulations.
The value at r=0 yields the clumping factor
which we define by:
 |
(40) |
We obtain
at z=0. Finally, we can check that our results agree with the outcome of the numerical simulations described in Dave et al. (2001) (see their Fig. 7). In particular, we recover the knee at
kpc and a reasonable value for the clumping factor
.
Note however that the dispersion in the
phase-diagram may actually lead to a slightly higher value.
5 Soft X-ray background
Finally, we check in this section whether the soft X-ray background emitted by the "warm'' IGM in our model is consistent with observations. Experiments show that the extragalactic soft X-ray background flux in the 0.1-0.4 keV band is of order 20-35 keV cm-2 s-1 sr-1 keV-1 (e.g., Warwick & Roberts 1998). However, most of this flux is due to AGN so that the contribution from diffuse gas should obey the constraint
keV cm-2 s-1 sr-1 keV-1 in this frequency range (e.g., see the discussion in Wu et al. 2001). Using the Press-Schechter formalism (Press & Schechter 1974), assuming that the X-ray emitting gas is embedded within just-virialized halos with a density contrast
,
some previous studies obtained a flux which is higher than this upper bound (e.g., Pen 1999; Wu et al. 2001). This led these authors to infer that a non-gravitational heating source is needed to unbind the gas from groups and small clusters in order to reduce their contribution to the soft X-ray background. However, the diffuse "warm'' IGM component is made of filamentary structures with density contrasts lower than
.
This gives a smaller X-ray background since the X-ray emission is proportional to the squared density of the gas. Moreover, the fraction of matter within this "warm'' diffuse phase is smaller than
,
see Fig. 3. As argued in Dave et al. 2001 this reduced X-ray background could then be consistent with observational constraints. This conclusion also agrees with the outcome of numerical simulations (Croft et al. 2001; Phillips et al. 2001).
Therefore, we compute here the soft X-ray background due to the diffuse "warm'' IGM phase, in the 0.1-0.4 keV band. Note that the X-ray background due to resolved sources (AGN, clusters, cooling galaxies) was already studied in Valageas & Schaeffer (2000) from an analytic model similar to the one used here to describe the underlying dark matter density field. The X-ray flux
within this frequency band can be written:
 |
|
|
(41) |
where the temperature T is in Kelvin units. Here we used again Eq. (32) with the normalization
defined in Eq. (35). The factor
in Eq. (41) comes from the fact that the X-ray emission is proportional to the squared density
.
The integral over x in Eq. (41) follows the curve
,
from Eq. (23), in the
phase-diagram. Then, our calculation gives an X-ray background flux
keV cm-2 s-1 sr-1 in the 0.1-0.4 keV band. This yields a mean differential flux
keV cm-2 s-1 sr-1 keV-1 in this band. Thus, we can check that our result is consistent with observations. Note that the scatter of the "warm'' IGM in the
plane should lead to a slightly larger flux. Therefore, this component might be observable through its X-ray emission in the near future.
In this article we have shown that the
phase-diagram of cosmological baryons can be understood in very simple terms. It is strongly model-independent since the relevant scales are set by atomic physics (hydrogen ionization) and the basic cosmological setting (the rate of expansion and the scale which marks the transition to the non-linear regime). Building a bridge between the modeling of galaxies, groups and clusters, as well as Lyman-
absorbers, we have constructed in this paper a phase-diagram which includes also the IGM, and thus provides a phase-diagram of all baryons. In this diagram, all these objects form well-defined populations.
Firstly, we have distinguished a "cool'' IGM phase corresponding
to the Lyman-
forest. It follows a well-defined equation of state on
the
plane. We have shown that the properties of this gas do not
depend much on its previous history because of the form of the recombination
coefficient. This explains why the scatter of this component obtained in
numerical simulations is quite small. The fraction of baryons in this components, as well as its evolution with redshift is seen to be consistent with earlier modelling of the frequency of occurrence of the Lyman-
lines. We however obtain (consistently with the numerical simulations) 38%
of the baryons in the cool IGM at z=0.
This may be a little high with
a mass fraction closer to 30%
being consistent with observations.
This must be considered as being within the error bars of our analytical calculation (and the error bars of the simulations!).
Secondly, we find that a "warm'' IGM phase is formed by the gas which has
been shock-heated to larger temperatures
K as non-linear
gravitational structures appear. The dependence on the stochastic
gravitational potential entails a broad scatter in the
plane. However, we have explained that robust constraints only allow a
closed region in the
phase-diagram for this component. We have
also defined a simple curve which represents its mean behaviour and plays
the role of the "warm'' IGM Equation of State. The latter is found to be
consistent with the outcome of numerical simulations (Dave et al. 2001), but
differs from the simple power-law fit of a mean curve done there. This is
because its locus is based on physical considerations, which we have
thoroughly justified. The "warm'' IGM is seen to be due to large structures
on their way of collapsing under the action of gravity. The quite simple
picture behind this model is that the collapse of an object turning
non-linear (that is with average overdensity above 5) induces shocks which
heat patches of much smaller size, to a temperature of the order of the
locally available kinetic energy of the collapse, creating local conditions within hot spots
close to hydrostatic equilibrium. We are able to estimate this size, as well
as the temperature of the latter, and to show they are in agreement with the
simulations. Definitely, these objects are not just Lyman-
absorbers at larger densities. The hot gas which lies
within clusters and galactic halos, on the other hand, corresponds to the high-density
"continuation'' (e.g.,
)
of the "warm'' IGM curve, in a
region of phase-space where virialization insures that all the available gas
has been gravitationally heated.
Next, we find two quite different regimes for the "warm'' IGM. At high
redshift, some of the "warm'' component enters the cooling region of the
plane: this gives rise to galactic disks and stars.
Some of the baryonic
matter, also, should lie on the low temperature branch of this cooling region,
in agreement with the results of numerical simulations (Dave et al. 1999). At low
redshift, the "warm'' IGM no longer goes into the cooled region,
and it provides
the origin of the hot gas in clusters. The evolution with redshift of
the baryon phase-diagram provides a natural confirmation that at high z
collapsed halos form galaxies while at low redshift they build groups or
clusters. The warm phase occupies 24%
of the baryon fraction at z=0 in
our analytic model.
Then, we have checked that our results for the fraction of matter enclosed
within the various phases and the two-point correlation function of the
"warm'' IGM component agree with numerical simulations (e.g., Dave et al. 1999
and Dave et al. 2001). This confirms the validity of our analysis. Note however
that the latter relies on simple physical considerations and it is
independent of the findings of the former. Besides, the soft X-ray background
due to the "warm'' IGM is consistent with the upper bound set by
observations. Our prediction actually is not far from the observational
limit: this offers the prospect of measuring this X-ray emission in future
observations.
We predict that approximately 60%
of the baryons are accounted for at z=0 in the cool/warm phases of the intergalactic medium. The remaining baryons (i.e. 40%)
are in collapsed structures. With a total of
,
this amonts to
.
The observed luminous component of the baryons (stars, remnants and gas) is estimated to be about
,
as noted in the introduction, and the hot gas in clusters represents
.
Therefore, there remains an "unidentified'' or dark baryon component,
.
This is at least twice the stellar-like component but only about 1/4 of the "identified'' baryons (IGM, stellar components and very hot gas). In our model this dark baryonic matter corresponds to (cool?) gas within galactic halos and groups, which has not been observed yet.
Acknowledgements
We thank R. Dave for providing us with some of the data published in Dave et al. 1999.
Appendix A: Dark matter density field
In this appendix we describe the analytic procedure we use to model the dark matter density field. Our method provides a simple means to estimate the pdf of the density contrast from the quasi-linear up to the highly non-linear regime, parameterized by only one variable (the skewness). This should be sufficient for our purpose which is merely to obtain the main properties of the IGM. Note that by construction the variance and the skewness of the pdf we obtain below agree with the results of numerical simulations, in both the quasi-linear and the highly non-linear regimes.
A.1 The pdf
of the overdensity at scale R
First, we note
the moment of order two of the density contrast
over a spherical cell of radius R and volume V. It can be expressed in terms of the two-point correlation function
as:
 |
(A.1) |
The two-point correlation
is the Fourier transform of the non-linear power-spectrum P(k) which we compute from the linear power-spectrum
using the analytic formulae obtained by Peacock & Dodds (1996) from fits to N-body simulations. Of course, in the quasi-linear limit
we have
.
Then, it is convenient to introduce the generating function
defined by the inverse Laplace transform:
![\begin{displaymath}%
{\cal P}(\rho_{\rm R}) \equiv \int_{-i\infty}^{+i\infty}\fr...
...i}} \; {\rm e}^{[\rho_{\rm R}y - \varphi(y)]/\overline{\xi}} ,
\end{displaymath}](/articles/aa/full/2002/24/aa2143/img261.gif) |
(A.2) |
where
is the pdf of the overdensity over a cell of radius R:
 |
(A.3) |
The advantage of the generating function
is that it reaches a finite limit in the quasi-linear regime
.
Moreover, in this regime one can derive
in a rigorous manner from the equations of motion (Bernardeau 1992, 1994; Valageas 2002). On the other hand, in the highly non-linear regime
the generating function
should also reach a finite (astonishingly close to the former, but seemingly not identical) limit if the stable-clustering ansatz is valid (Balian & Schaeffer 1989). This is seen to agree reasonably well with the results of numerical simulations (e.g., Valageas et al. 2000). Therefore, hereafter we use the following simple parameterization for the generating function
:
 |
(A.4) |
with:
 |
(A.5) |
This yields a family of functions
which depends only on the parameter
.
Such a relationship arises naturally in the exact derivation of
in the quasi-linear limit (Bernardeau 1992, 1994; Valageas 2002). Besides, it also appears in the non-linear regime within the framework of the "minimal tree-model'' (Bernardeau & Schaeffer 1992; Schaeffer 1985). It can be seen that the phenomenological form (A.5) for
also yields reasonable results.
As in Bernardeau & Valageas (2000), we choose the parameter
which enters Eq. (A.5) so as to recover the skewness
predicted by HEPT (Scoccimarro & Frieman 1999) in the non-linear regime. More precisely, we use:
 |
(A.6) |
and we interpolate between both values in the range
.
Here we note n the slope of the linear-power spectrum at the scale of interest. Note that this procedure was seen in Bernardeau & Valageas (2000) to match the results of numerical simulations with respect to weak gravitational lensing effects. Therefore, it should be quite sufficient for our present purpose which is mainly to investigate the main trends of the IGM properties.
A.2 High and low density cutoffs
(R)
Now, we can estimate the high and low density cutoffs
of the pdf
as follows. To each overdensity
we associate the Laplace variable
(and
through Eq. (A.4)) defined as the saddle-point of the exponent in Eq. (A.2). This can be written:
 |
(A.7) |
where we used Eq. (A.4) to derive the last equality. Note that
and the pair
are of opposite signs. The reason for the introduction of
is that for extreme overdensities
or
,
that is in the tails of the pdf
,
the integral (A.2) is dominated by the contribution from
(e.g., Balian & Schaeffer 1989). Then, it is clear that the density cutoffs of the pdf are given by the points
where the argument of the exponential in Eq. (A.2) at the saddle-point is of order unity and negative (note that it is zero at the mean overdensity
). Indeed, for larger
the argument becomes more negative than -1 which leads to the exponential falloffs of
.
Therefore, we define the high and low density cutoffs
by the condition:
 |
(A.8) |
where
as given by Eq. (A.7). This relation can be most conveniently written in terms of the variable
which yields:
 |
(A.9) |
where we used Eq. (A.5). In the highly non-linear regime Eq. (A.9) yields
.
However, as discussed in Bernardeau & Schaeffer 1992 the implicit system (A.4) is singular at a finite value
which leads to an exponential tail for the pdf
.
Hence for
we take
and
.
Thus, we obtain for the high-density cutoff
:
 |
(A.10) |
and for the low-density cutoff
:
 |
(A.11) |
Note that in the quasi-linear regime we recover
as required for Gaussian initial conditions. Then, the density cutoffs
are close to the mean:
for
.
On the other hand, in the highly non-linear regime the cutoffs are displaced to very low densities (
:
rare voids) and high densities (
:
collapsed halos).
-
Anninos, P., Zhang, Y., Abel, T., & Norman, M. L. 1997, NewA, 2, 209
In the text
NASA ADS
-
Balian, R., & Schaeffer, R. 1989, A&A, 220, 1
In the text
NASA ADS
-
Bernardeau, F. 1992, ApJ, 392, 1
In the text
NASA ADS
-
Bernardeau, F. 1994, A&A, 291, 697
In the text
NASA ADS
-
Bernardeau, F., & Schaeffer, R. 1992, A&A, 255, 1
In the text
NASA ADS
-
Bernardeau, F., & Valageas, P. 2000, A&A, 364, 1
In the text
NASA ADS
-
Cen, R., & Ostriker, J. P. 1999, ApJ, 514, 1
In the text
NASA ADS
-
Cooke, A. J., Espey, B., & Carswell, R. F. 1997, MNRAS, 284, 552
In the text
NASA ADS
-
Croft, R. A. C., di Matteo, T., Dave, R., et al. 2001, ApJ, 557, 67
In the text
NASA ADS
-
Dave, R., Hernquist, L., Katz, N., & Weinberg, D. H. 1999, ApJ, 511, 521
In the text
NASA ADS
-
Dave, R., Cen, R., Ostriker, J. P., et al. 2001, ApJ, 552, 473
In the text
NASA ADS
-
Fukugita, M., Hogan, C. J., & Peebles, P. J. E. 1998, ApJ, 503, 518
In the text
NASA ADS
-
Giallongo, E., Cristiani, S., D'Odorico, S., Fontana, A., & Savaglio, S.
1996, ApJ, 466, 46
In the text
NASA ADS
-
Hui, L., & Gnedin, N. Y. 1997, MNRAS, 292, 27
In the text
NASA ADS
-
Nath, B. B., & Silk, J. 2001, MNRAS, submitted [astro-ph/0107394]
In the text
-
Peacock, J. A., & Dodds, S. J. 1996, MNRAS, 280, L19
In the text
NASA ADS
-
Pen, U.-L. 1999, ApJ, 510, L1
In the text
NASA ADS
-
Phillips, L. A., Ostriker, J. P., & Cen, R. 2001, ApJ, 554, L9
In the text
NASA ADS
-
Press, W. H., & Schechter, P. 1974, ApJ, 187, 425
In the text
NASA ADS
-
Rees, M. J., & Ostriker, J. P. 1977, MNRAS, 179, 541
In the text
NASA ADS
-
Scoccimarro, R., & Frieman, J. A. 1999, ApJ, 520, 35
In the text
NASA ADS
-
Schaeffer, R. 1985, A&A, 144, L1
In the text
NASA ADS
-
Silk, J. 1977, ApJ, 211, 638
In the text
NASA ADS
-
Springel, V., & Hernquist, L. 2001, MNRAS, submitted [astro-ph/0111016]
In the text
-
Tytler, D., O'Meara, J. M., Suzuki, N., & Lubin, D. 2000, Phys. Scr. T, 85, 12 In the text
-
Valageas, P. 2002, A&A, 382, 412
In the text
-
Valageas, P., & Schaeffer, R. 1997, A&A, 328, 435
In the text
NASA ADS
-
Valageas, P., & Schaeffer, R. 1999, A&A, 345, 329
In the text
NASA ADS
-
Valageas, P., & Schaeffer, R. 2000, A&A, 359, 821
In the text
NASA ADS
-
Valageas, P., & Silk, J. 1999, A&A, 350, 725
In the text
NASA ADS
-
Valageas, P., Schaeffer, R., & Silk, J. 1999, A&A, 345, 691
In the text
NASA ADS
-
Valageas, P., Lacey, C., & Schaeffer, R. 2000, MNRAS, 311, 234
In the text
NASA ADS
-
Valageas, P., Balbi, A., & Silk, J. 2001a, A&A, 367, 1
In the text
NASA ADS
-
Valageas, P., Schaeffer, R., & Silk, J. 2001b, A&A, 366, 363
In the text
NASA ADS
-
Valageas, P., Schaeffer, R., & Silk, J. 2002, in preparation
In the text
-
Vogel, S. N., Weymann, R., Rauch, M., & Hamilton, T. 1995, ApJ, 441, 162
In the text
NASA ADS
-
Warwick, R. S., & Roberts, T. P. 1998, Astron. Nachrichten, 319, 59
In the text
NASA ADS
-
Wu, K. K. S., Fabian, A. C., & Nulsen, P. E. J. 2001, MNRAS, 324, 95
In the text
NASA ADS
-
Zel'dovich, Ya. B. 1970, A&A, 5, 84
In the text
NASA ADS
Copyright ESO 2002