format: text

text:
# Detecting X-ray filaments in the low redshift Universe with XEUS and Constellation-X

## 1 Introduction

The census of baryons in the low redshift Universe shows that a large fraction of them has not been detected yet. While almost 80% of the baryons reside in the Ly$`\alpha `$ forest at $`z2`$ (Rauch 1998), estimates at lower redshift imply that the observed baryon density is significantly lower (Fukugita et al. 1998; Nicastro et al. 2002b) and does not match the nucleosynthesis constraints (Burles & Tytler 1998). Hydro-dynamical simulations suggest that a large fraction of baryons at $`z0`$ could be in a warm-hot intergalactic medium (WHIM), with a temperature between $`10^5`$ and $`10^7`$ K and with moderate overdensities $`\delta 10100`$ (Ostriker & Cen 1996; Cen & Ostriker 1999; Davé et al. 2001; Cen et al. 2001). Detection of this component is particularly difficult and requires very sensitive UV and X-ray satellites.

Recently, various attempts to model analytically the WHIM component have been made. Perna & Loeb (1998) used the Press & Schechter theory (Press & Schechter 1974) and an isothermal model for the gas component to address the detectability in absorption of the hot and highly ionized intergalactic medium (IGM), which resides in the outskirts of clusters and groups of galaxies. Valageas, Schaeffer & Silk (2002) showed that the temperature-density relation found in hydro-dynamical simulations can be reproduced using straightforward physical arguments. In this simple picture the IGM consists of two phases: the cool photo-ionized component which gives rise to the Ly$`\alpha `$ forest and the WHIM, dominated by shock heating.

Hydro-dynamical simulations of structure formation in the context of cold dark matter models have been the major contribution to the theoretical understanding of the WHIM. Hellsten et al. (1998), referred to the network of filaments and sheet-like structures, possibly seen in the spectra of background bright sources in the X-ray band, as the ‘X-ray forest’, in analogy with the Ly$`\alpha `$ forest. Their main result is that OVII and OVIII absorptions by the intervening IGM can be seen by future X-ray missions. Cen et al. (2001) used large box-size hydro-dynamical simulations to detect OVI absorption by the WHIM in the UV band and find results in reasonably good agreement with observations. Fang et al. (2002a) performed hydro-dynamical simulations to study the metal distribution in the IGM and, assuming collisionally ionized gas, compared their results with semi-analytical models based on the Press & Schechter formalism. Kravtsov et al. 2002 addressed WHIM detectability using hydro-dynamical simulations to model the gas distribution in the Local Supercluster region and imposing constraints from the MARK III catalogue of galaxy peculiar velocity.

More recently, Chen et al. (2002) have made an extensive study of the detectability of OVII and OVIII absorbers using the $`z=0`$ output of a hydro-dynamical simulation of a $`\mathrm{\Lambda }`$CDM Universe, and assuming a uniform metallicity of 0.1 solar for the IGM. Their results confirm that the detectability of these absorbers is challenging. It has been shown that the number of intervening OVI absorbers per unit redshift is quite high at low redshift (Tripp et al. 2000). Thus, a very promising technique to detect OVII and OVIII systems is to search for corresponding UV OVI absorptions in the spectra: this allows to lower the S/N ratio needed for OVII and/or OVIII detectability. Nicastro et al. (2002a) presented the X-ray 5$`\sigma `$ level detection by Chandra of a resonant absorption from WHIM along the line of sight towards the blazar PKS 2155-304. Moreover, Fang et al. (2002b) reported a detection of an OVIII absorption line along the sight-line towards PKS 2155-304 with Chandra. They constrained the gas that gives rise to this line to have an overdensity in the range 30-350 and a temperature of 4-5 $`\times 10^6`$ K, in agreement with hydro-simulations. Also emission by the WHIM could be important and several models have been proposed (e.g. Scharf et al. 2000; Phillips et al. 2001; Zappacosta et al. 2002). Here, we will focus on absorption: the bulk of the missing baryons should reside in a medium which is too rarefied to be detectable through its X-ray emission or Sunyaev-Zel’dovich effect (Ostriker & Cen 1996).

These very recent observations are very encouraging and show that we have just started to detect the first ‘trees’ of the X-ray forest. Future X-ray missions such as XEUS and Constellation-X will be much more accurate a probe of the current cosmological models and will allow a detailed study of the WHIM properties.

In this paper we will present a semi-analytical technique to simulate X-ray absorption spectra. There are two advantages in using a semi-analytical technique instead of the more accurate hydro-dynamical simulations. First, semi-analytical models are not computationally expensive and are not limited by resolution and box-size effects. Second, semi-analytical techniques can be very efficient when coupled with detailed hydro-dynamical simulations when the latter can be used to calibrate the former and allow one to explore very quickly the parameter space. For example, it will be straightforward to study the effects of the metallicity of the IGM, different UV and X-ray backgrounds, scatter in the temperature-density and in the metallicity-density relation, peculiar velocities etc. For the estimates we want to make here, semi-analytical techniques, which use the inputs of hydro-dynamical simulations, will show to be successful in reproducing potentially observable quantities such as the column densities of OVII and OVIII absorbers.

In particular, we will focus on the simulations of absorption spectra of AGN pairs and we will try to determine the characteristic size of the absorbers by identifying absorption features in the simulated spectra. The estimate of the extent of the absorbers, using the information contained in the transverse direction, could be a very powerful probe of the filamentary network of WHIM. An extensive analysis of the correlations of Ly$`\alpha `$ clouds at $`z>2`$ using hydro-dynamical simulations has been made in Miralda-Escudé et al. (1996) and Cen & Simcoe (1997). McGill (1990), Crotts et al. (1994), Bechtold et al. (1994), Charlton et al. (1997) performed similar analysis for observed Ly$`\alpha `$ clouds and showed that the coherence length in the transverse direction of the clouds could be significantly larger than the size of a single cloud. This argument has been clarified recently using simple physical arguments (Schaye 2001) and hydro-dynamical simulations (e.g. Theuns et al. 1998), which indeed show that the clouds are embedded in a network of filaments. The idea is very simple: while a single line-of-sight (LOS) allows to probe only the size of a single absorber, pairs of LOS are sensitive to the random orientation of the filaments and diffuse blobs and can better sample the WHIM. Another advantage of using the information contained in a sample of pairs is that this method minimizes the contamination of the WHIM signal we have to detect with the intrinsic lines of the background source. This is due to the fact that the positions of the intrinsic absorption or emission features of the source along the LOS are random. Thus, their effect can be statistically modelled and removed as it has been done recently to recover the linear dark-matter power spectrum from simulated Ly$`\alpha `$ absorptions (Viel et al. 2002a).

The outline of the paper is as follows. In Section 2 we describe the ingredients needed to simulate the X-ray AGN spectra focussing both on the inputs from hydro-dynamical simulations (subsection 2.1) and on the analytical modelling (subsection 2.2). In Section 3 we test how well the semi-analytical model matches the properties of simulated lines extracted directly from hydro-simulations. This comparison is done in terms of the predicted number counts of absorbers per unit redshift range. In Section 4 we discuss a possible way to recover the characteristic size of these absorbers by using the so called ‘hits-and-misses’ statistics and discuss detectability of these absorption features by satellites of the next generation such as XEUS and Constellation-X. Section 5 contains our main conclusions.

## 2 Semi-analytical model of the X-Ray forest

This Section presents the semi-analytical method used to simulate X-ray spectra of AGN. In Subsection 2.1 we will present the hydro-dynamical simulation used, focussing on the basic ingredients which will constitute the inputs for the semi-analytical model: i) the probability distribution function (pdf) of the gas; ii) the gas metallicity - density relation; iii) the gas temperature - density relation; iv) the X-ray background. In Subsection 2.2 we will describe in detail the technique used, with the relative mathematical implementation (see also Appendix A1), to produce the simulated spectra of AGN pairs.

### 2.1 Model parameters from hydro-dynamical simulations

The hydro-dynamical simulation used here is a $`\mathrm{\Lambda }`$CDM model with the following parameters: $`H_0=100h`$ km s-1Mpc-1 with $`h=0.67`$, $`\mathrm{\Omega }_{0m}=0.30`$, $`\mathrm{\Omega }_{0b}=0.035`$, $`\mathrm{\Lambda }=0.70`$, $`\sigma _8=0.90`$, and the spectral index of the primordial power spectrum $`n=1`$. The box-size of the simulation is $`25`$ comoving $`h^1`$ Mpc on a uniform mesh with $`768^3`$ cells. The comoving cell size is 32.6 $`h^1`$ kpc and the mass of each dark matter particle is $`2\times 10^7`$ M (further details can be found in Cen et al. 2001). The simulation includes galaxy and star formation, energy feedback from supernova explosions, ionization radiation from massive stars and metal recycling due to SNe/galactic winds. Metals are ejected into the local gas cells where stellar particles are located using a yield $`Y=0.02`$ and are followed as a separate variable adopting the standard solar composition (some slices through the simulation box at various redshift outputs can be found at http://astro.Princeton.EDU/$``$cen/PROJECTS/p2/p2.html).

To properly simulate X-ray absorptions we have to know the actual X-ray and UV background of the simulation at $`z=0`$, which is computed self-consistently given the sinks and the sources in the simulation box. We have checked that there are very small differences between the hydro-simulation background and the estimate from Shull et al. (1999): $`I_{UV}=I_{UV}^0(E/13.6\mathrm{eV})^{1.8}`$, with $`I_{UV}^0=2.3\times 10^{23}`$ erg cm-2 Hz-1 sr-1 s-1 for the UV background (which includes a contribution from AGNs and starburst galaxies), and $`I_X=I_X^0(E/E_X)^{1.29}\mathrm{exp}(E/E_X)`$, with $`I_X^0=1.75\times 10^{26}`$ erg cm-2 Hz-1 sr-1 s-1 and $`E_X=40`$ keV (Boldt 1987; Fabian & Barcons 1992; Hellsten et al. 1998; Chen et al. 2002), for the X-ray background. The difference between the amplitude of the UV and X-ray background of the hydro-simulations and the estimates above is less than a factor of 2, in the range 0.003 $``$$`<`$ E (keV) $``$$`<`$ 4. In Figure 1 we plot the X-ray and UV background of the hydro-simulations at $`z=0`$ (dashed line), the estimates of Shull et al. (1999) (long dashed line) and Fabian & Barcons 1992 (dot-dashed line), and the sum of these two (continuous line).

Among the two backgrounds, the X-ray one plays a major role in determining the statistics of the simulated spectra. We will assume that the amplitude of the ionizing backgrounds scales like $`(1+z)^3`$, in the redshift range $`0<z<1`$. This assumption will not influence significantly our results: given the fact that we will generate mock-spectra at $`z<0.3`$, this will determine a scaling of the amplitude by a factor of less than 2.

Along many randomly chosen sight-lines parallel to one of the axes of the simulation box, we compute the gas density and peculiar velocity, the gas temperature, and the metal density. This allows us to obtain a sample of 1000 simulated absorption spectra.

From the $`z=0`$ output we extract the pdf of the gas density, which is shown in Figure 2, this will be implemented in the semi-analytical (hereafter SA) model. Starting from the linear density field, generated with the right correlation properties predicted by linear theory (Subsection 2.2), we assign each pixel a new non-linear density in such a way that the resulting non-linear pdf is that extracted from hydro-dynamical simulations.

The gas density-temperature relation is shown in panel $`(a)`$ of Figure 3, where we report the mean value in each density bin and the scatter. This result is in reasonable agreement with the recent results obtained by Yoshida et al. (2002 - their Figure 2) and also by Davé et al. (1999), who studied the gas cooling processes in the framework of Smoothed Particle Hydro-dynamics (SPH) simulations. We note that the scatter in temperature can be very high for values of $`\mathrm{log}(1+\delta )>1`$, spanning more than one order of magnitude. We report the gas density-metallicity relation in panel $`(b)`$. The trend here is as expected: low density regions are the most metal poor, while at large densities the metallicity can be solar. In this plot the scatter is extremely large for the intermediate density regime, i.e. $`1<\mathrm{log}(1+\delta )<0.5`$, while it is significantly smaller for larger overdensities. This plot has to be compared with similar plots shown in Cen & Ostriker (1999b), but we find that here the scatter is somewhat larger than in previous simulations.

These curves, with the relative scatter, assumed to be drawn from a Gaussian distribution, will be implemented into the SA model.

### 2.2 Simulation of absorption features in X-ray spectra of AGN pairs

We start here from the model introduced by Bi and collaborators (Bi et al. 1992, 1995; Bi 1993; Bi & Davidsen 1997 - hereafter, BD97) for generating a Ly$`\alpha `$ absorption spectrum along a line-of-sight. This simple model predicts many properties of the absorption lines, including the column density distribution and the distribution of line widths ($`b`$ parameters), which can be directly compared with observations.

The BD97 model is based on the assumption that the low-column density Ly$`\alpha `$ forest is produced by smooth fluctuations in the intergalactic medium which arise as a result of gravitational growth of perturbations. Linear density perturbations in the intergalactic medium $`\delta _0^{\mathrm{IGM}}(𝐱,z)`$ are related to the underlying linear dark-matter (DM) perturbations by a convolution, which models the effects of gas pressure. In Fourier space one has: $`\delta _0^{\mathrm{IGM}}(𝐤,z)=W_{\mathrm{IGM}}(k,z)D_+(z)\delta _0^{\mathrm{DM}}(𝐤)`$, where $`D_+(z)`$ is the growing mode of density perturbations (normalised so that $`D_+(0)=1`$) and $`\delta _0^{\mathrm{DM}}(𝐤)`$ is the Fourier transformed DM linear over-density at $`z=0`$. A commonly adopted low-pass filter is $`W_{\mathrm{IGM}}(k,z)=(1+k^2/k_J^2)^1`$, which depends on the comoving Jeans wave-length $`\lambda _J(z)2\pi H_0^1\left[\frac{2\gamma k_BT_m(z)}{3\mu m_p\mathrm{\Omega }_{0m}(1+z)}\right]^{1/2}`$, where $`k_B`$ is Boltzmann’s constant, $`T_m`$ the gas temperature, $`\mu `$ the molecular weight and $`\gamma `$ the ratio of specific heats. $`\mathrm{\Omega }_{0m}`$ is the present-day matter density parameter. BD97 adopt a simple lognormal transformation to obtain the IGM density in the mildly non-linear regime from the linear one.

Several approximations are involved in this modelling. A first point which can be important is that the Jeans length depends on density and temperature, and therefore Jeans smoothing should be an adaptive smoothing of the density field, while the Fourier implementation cannot take into account this local aspect (Viel et al. 2002b). In addition, a filtering of the dark-matter density field at the Jeans scale may be not very accurate: different filtering functions result in very different gas density fields and the real ‘filtering scale’ which reproduces better the gas distribution, in the linear regime, depends on the thermal history of the IGM (Hui & Gnedin 1998; Matarrese & Mohayaee 2002; Gnedin et al. 2002).

The modelling of the gas component using the lognormal transformation and a smoothing at the Jeans length would determine a poor description of the gas distribution at $`z0`$, where the density fluctuations are in the non-linear regime. Even the tight relation between density and temperature for the gas responsible of Ly$`\alpha `$ absorption is no longer valid for the warm-hot intergalactic medium. In this case, hydro-simulations suggest that the scatter in the $`T\delta `$ relation, for the WHIM, can be very large (Davé et al. 1999; Yoshida et al. 2002).

At least three inputs, taken from hydro-dynamical simulations, can be used to improve the model:  i) the probability distribution function (pdf) of the gas density;  ii) the gas density-temperature relation;  iii) the gas density-metallicity relation. Viel et al. (2002b) showed that by using directly the pdf of the gas density extracted from the hydro-dynamical simulations, the pdf of the flux predicted by the semi-analytical model is in significant better agreement with that in the hydro-simulations. The other two inputs will allow to take into account the physical state (temperature and metallicity) of the structures responsible for the absorptions.

Here we will follow closely the prescriptions used in Viel et al. (2002a, 2002b). The modelling is based on the generation in Fourier space of 1D density and velocity fields, with the right correlation properties predicted by the theory.

Starting from a linear power spectrum $`P(k)`$ we generate 1D LOS linear random fields for the IGM density and peculiar velocity, with the low-pass filter window function $`W_{IGM}`$. If we draw a LOS in the $`x_{}`$ direction at a fixed coordinate $`x_{}`$ the Fourier component at redshift $`z`$ of the linear density contrast and of the linear peculiar velocity field of the intergalactic medium (IGM) along the LOS are:

$`\delta _0^{\mathrm{IGM}}(k_{},z|𝐱_{})`$ $`=`$ $`D_+(z){\displaystyle \frac{d^2k_{}}{(2\pi )^2}e^{i𝐤_{}𝐱_{}}W_{\mathrm{IGM}}(\sqrt{k_{}^2+k_{}^2},z)\delta _0^{\mathrm{DM}}(k_{},𝐤_{})},`$ (1)
$`v_{}^{\mathrm{IGM}}(k_{},z|𝐱_{})`$ $`=`$ $`ik_{}E_+(z){\displaystyle \frac{d^2k_{}}{(2\pi )^2}e^{i𝐤_{}𝐱_{}}\frac{1}{k_{}^2+k_{}^2}W_{\mathrm{IGM}}(\sqrt{k_{}^2+k_{}^2},z)\delta _0^{\mathrm{DM}}(k_{},𝐤_{})},`$ (2)
with $`E_+(z)=H(z)f(\mathrm{\Omega }_m,\mathrm{\Omega }_\mathrm{\Lambda })D_+(z)/(1+z)`$. Here $`f(\mathrm{\Omega }_m,\mathrm{\Omega }_\mathrm{\Lambda })d\mathrm{ln}D_+(z)/d\mathrm{ln}(1+z)`$), and $`\mathrm{\Omega }_m`$ and $`\mathrm{\Omega }_\mathrm{\Lambda }`$ are the matter and vacuum-energy contribution to the cosmic density (e.g. Lahav et al. 1991), respectively; $`H(z)`$ is the Hubble parameter at redshift $`z`$: $`H(z)=H_0\sqrt{\mathrm{\Omega }_{0m}(1+z)^3+\mathrm{\Omega }_{0r}(1+z)^2+\mathrm{\Omega }_{0\mathrm{\Lambda }}}`$, with $`\mathrm{\Omega }_{0r}=1\mathrm{\Omega }_{0m}\mathrm{\Omega }_{0\mathrm{\Lambda }}`$. In our case we use the $`\mathrm{\Lambda }`$CDM parameters of the previous subsection.

We then simulate the corresponding density and peculiar velocity at a distance $`r_{}`$ from the first LOS. The correlation properties between these 4 fields are computed using an algebraic implementation, as described in Appendix A1.

The linear density field is mapped into the non-linear one using directly the pdf of the gas distribution obtained from hydro-simulations. This has been done with a rank-ordering technique between the linear density field and the non-linear one, which intrinsically consists in a monotonic and deterministic mapping of each generated field. To account for the scatter, we assign each pixel a new temperature and metallicity value by adding to the mean value (continuous lines of Figure 3) a random number taken from a Gaussian distribution with the same standard deviation as the hydro-simulations.

Peculiar velocities are approximated by those predicted by linear theory, as the peculiar velocity field is known to keep linear even on scales where the density contrast gets mildly non-linear. This simplification should not be critical in simulating the spectra, as the effect of peculiar velocities results in a shift of the line and in an alteration of the profile (BD97), which will have small impact on the recovered column density of the absorber.

After having simulated the gas distribution, the photoionization code CLOUDY (Ferland et al. 1998) is used to compute the ionization states of metals. In this analysis, we have considered both the collisional ionization and the photo-ionization. The ionizing background consists of the sum of a UV background, arising from AGNs and galaxies, and a X-ray background, probably produced by a population of moderately X-ray luminous AGN at $`z2`$ (Fabian & Barcons 1992), while for the soft X-ray part ($`<1`$ keV) a significant contribution comes from the WHIM itself (Cen et al. 1995; Phillips et al. 2001).

Among the heavy elements that can produce absorption lines in the X-rays Oxygen is the most abundant one and produces the strongest lines. A more complete treatment should also include other elements such as C, N, Ne, Fe and Si but, in our simplified scheme, we will concentrate only on Oxygen lines. For a full list of resonant absorption lines we refer to Verner et al. (1996). We simulate absorption spectra for the two strongest transitions of the ions OVII ($`E=0.57`$ keV, $`\lambda =21.6\AA `$, $`f=0.7`$) and OVIII ($`E=0.65`$ keV, $`\lambda =18.97\AA `$, $`f=0.42`$), with $`f`$ the oscillator strength.

The density of each ion is obtained through: $`n_I(x)=n_H(x)X_IY_Z_{}Z/Z_{}`$, where $`X_I`$ is the ionization fraction of the ion as determined by CLOUDY and depends on gas temperature, gas density and ionizing background, $`n_H`$ is the density of hydrogen atoms, Z is the metallicity of the element and $`Y_Z_{}`$ is the solar abundance of the element. In Figure 4 we plot the ionization fraction for OVII and OVIII as a function of the temperature for two different values of the IGM density corresponding to an overdensity $`\delta _{IGM}=6`$ (left panel) and $`\delta _{IGM}=120`$ (right panel), corresponding to densities of $`1.4\times 10^6`$ cm-3 and $`3\times 10^5`$ cm-3. To check if the metals considered here are in photoionization equilibrium we compute the recombination time $`t_{rec}=(n_{IGM}\alpha _{rec})^1`$, with $`\alpha _{rec}`$ the recombination rate. For OVII and OVIII, in the temperature ranges considered here, one has that $`\alpha _{rec}2\times 10^{11}`$ cm3 s-1 (Mazzotta et al. 1998), where both the dielectronic and radiative recombinations are taken into account. This means that the Hubble time is larger than the recombination time roughly for densities $`>\text{ }10^7`$ cm-3. Thus, for the WHIM, the approximation of photoionization equilibrium should be reasonable.

Given the density of a given ion along the LOS, the optical depth in redshift-space at velocity $`u`$ (in km s-1) is

$$\tau _I(u)=\frac{\sigma _{0,I}c}{H(z)}_{\mathrm{}}^{\mathrm{}}𝑑yn_\mathrm{I}(y)𝒱[uyv_{}^\mathrm{I}(y),b(y)]$$
(3)
where $`\sigma _{0,I}`$ is the cross-section for the resonant absorption and depends on $`\lambda _I`$ and $`f_I`$, $`y`$ is the real-space coordinate (in km s-1), $`𝒱`$ is the standard Voigt profile normalised in real-space, $`b=(2k_BT/m_Ic^2)^{1/2}`$ is the thermal width and we assume that $`v_I=v_{IGM}`$. Velocity $`v`$ and redshift $`z`$ are related through $`d\lambda /\lambda =dv/c`$, where $`\lambda =\lambda _I(1+z)`$. For the low column-density systems considered here, the Voigt profile is well approximated by a Gaussian: $`𝒱=(\sqrt{\pi }b)^1\mathrm{exp}[(uyv_{}^\mathrm{I}(y))^2/b^2]`$. The X-ray optical depth $`\tau `$ will be the sum of the source continuum, which we assume we can determine, and that of eq. (3). Finally, the transmitted flux is simply $`=\mathrm{exp}(\tau )`$.

We generate 400 spectra with our semi-analytical technique at a median redshift of $`z0.05`$ and $`12000`$ km/s long. Both for the semi-analytical model and for the hydro-simulations spectra have been produced with a resolution of $`3`$ km s-1. The goal is to reproduce number and characteristic size of these absorbers with the simulation of AGN spectra by a simple semi-analytical recipe, and address their detectability. Figure 5 shows the OVII distribution at $`z=0`$ extracted from the hydro-simulation. We can see that the OVII connects high density regions and displays a network of filaments and diffuse blobs. From this Figure we can argue that strong absorptions arise in virialized halos, while the weaker ones, possibly detected by future X-ray missions, sample filaments connecting high density regions. The OVIII distribution is quite similar to the OVII one (Chen et al. 2002).

## 3 Testing semi-analytical models with hydro-dynamical simulations

In this Section we show how well the semi-analytical model works, after having calibrated it with the inputs obtained from hydro-dynamical simulations: the pdf of the gas, the gas density-temperature relation and the gas density-metallicity relation (with the relative scatter) and the X-ray background. We compare the column density distribution function for OVII and OVIII predicted by the semi-analytical model with those extracted directly from the hydro-dynamical simulations. To this aim, the SA method has been run cutting off the large scale power for scales larger than the box-size of the hydro simulation (25$`h^1`$ Mpc). This allows a more direct comparison between the hydro-simulation and the semi-analytical model. We checked that the statistics considered here are only marginally influenced by this assumption, by running another set of SA spectra without the cut in the power-spectrum. We conclude that the box-size of the hydro-simulations is large enough to study reliably these absorbers. Column densities have been computed by integrating directly the simulated optical depth and not using any fitting routine available, both for the hydro-simulations and for the SA method. Usually, the number of absorbers found is small and they can be very easily identified in our mock-spectra. However, we checked that the estimated column densities are generally in good agreement with the more accurate results obtained with fitting routines.

In Figure 6 we show the cumulative number of absorbers per unit redshift range as a function of their column density (bottom axis), i.e. the number of systems with column density larger than a given value, and as a function of the corresponding equivalent width (top axis), for OVII (left panel) and OVIII (right panel): the continuous lines represent the results obtained from hydro-dynamical simulations, while the points are from the SA method. The absorption equivalent width, in velocity units, is related to the column density by the simple relation:

$$W=\frac{\pi e^2}{(m_ec)}\lambda fN$$
(4)
where $`N`$ is the column density, $`f`$ the oscillator strength and $`\lambda `$ the wavelength of the transition (Sarazin 1989, Chen et al. 2002). This relation formally holds in the limit of optically thin absorber but is also a good approximation for column densities smaller than $`10^{15}\mathrm{cm}^2`$.

To quantify the role of the scatter in our simulations, we decided to run another set of SA spectra and switch-off any form of scatter, both in the temperature and metallicity. Results without scatter are shown as empty circles, while those that include the scatter are represented by filled circles. From Figure 6 we see that, when the scatter is considered, the agreement between hydro simulations and SA is good. A first conclusion is that the inputs we used for the SA method, taken from hydro-simulations, are enough to have a good description of the number of absorbers per unit redshift range for the column density range $`10^{12.5}`$cm$`{}_{}{}^{2}<N_{OVII,OVIII}<10^{15.5}`$cm-2, corresponding to $`0.1\mathrm{km}/\mathrm{s}<\text{ }\mathrm{W}<\text{ }35`$ km/s for $`O_{VII}`$.

The values found are in agreement with the column densities distributions of Chen et al. (2002), which are reported in Figure 6 and are represented as a dashed line. In the case in which they take into account the metallicity-density relation, parameterized by $`\mathrm{log}Z/Z_{}=1.66+0.36\mathrm{log}\delta `$ and consider a scatter drawn randomly from a lognormal distribution with $`\mathrm{log}Z/Z_{}=1`$ and $`\sigma _{\mathrm{log}Z}=0.4`$, the results are quite similar. We note a slightly steeper column density distribution function in our analysis compared to the one in Chen et al. (2002), which results in a smaller number of absorbers at intermediate column densities in the range $`10^{13}10^{15}`$ cm-2. This is probably determined by the fact that our gas-density metallicity relation predicts lower metallicities than the one they used.

Before making any consideration about the detectability of these systems, we make some comments about the effect of the scatter. If we do not take into account the scatter, we cannot reproduce all the systems at column densities larger than $`10^{14}`$ cm-2 for both the ions. For smaller column densities the effect of the scatter is smaller, but it seems that also in this range the scatter produces an increase in the number of systems.

We checked the different role played by the metallicity and temperature scatter by running another set of SA spectra switching off only one of them. It turns out that the scatter that influences most the column density of the absorbers is the one in metallicity. This is probably due to the fact that the column density of the absorbers is directly proportional to the gas metallicity (Section 2.2). The dependence on the temperature is in a sense less direct, as the latter determines the ionization fraction, computed with CLOUDY including both the collisional ionization and the photo-ionization contributions (see Section 2.2), which depends on the particular X-ray background chosen in a non-trivial way (see for example Figure 4). Another reason could be that the scatter in metallicity is slightly larger than the temperature one for the regions with $`5<\delta <10`$, a low density WHIM, which is responsible for a significant fraction of the absorptions.

We expect $`>\text{ }30`$ OVII systems with column densities larger than $`10^{13}\mathrm{cm}^2`$ corresponding to equivalent widths larger than 2 km/s, per unit redshift range at $`z0`$. While the number of systems with column densities larger than $`10^{14.5}\mathrm{cm}^2`$, corresponding to equivalent widths larger than $`10`$ km/s, drops to 5. The number of OVIII systems is similar to that of OVII ones showing that both ions trace structures at the same density, but the corresponding OVIII equivalent width is almost a factor 2 lower.

## 4 Results

To better understand the physical state of the gas which determines the absorptions we plot in Figure 7 and Figure 8 two AGN pairs, for two different separations of 0.2 $`h^1`$ Mpc and 1.2 comoving $`h^1`$ Mpc (continuous line and dashed line represent the two members of the pair) simulated with the semi-analytical technique. Assuming that the absorbers are placed at a median redshift of $`z0.05`$, these distances correspond to an angular separation of $`3.5`$ arcmin and $`42`$ arcmin, respectively, for a $`\mathrm{\Lambda }`$CDM Universe. If the absorbers are placed at a median redshift of $`z0.1`$ the angular separation is $`2`$ arcmin and $`30`$ arcmin, respectively.

In Figure 7 and 8 it is possible to see that the absorptions arise from gas at a temperature larger than $`10^5`$ K which has a density larger than $`10^6`$ cm-3, i.e. overdensity $`\delta >6`$. One can see that, while for the smaller separation (Figure 7) the two AGNs show almost identical features both in density and temperature, which determine very similar spectra for OVII and OVIII, for the largest separation of Figure 8 the differences are larger. This is due to the fact that now the correlations are weaker, resulting in more different density and temperature profiles for the two AGN spectra. Also the spectra now show some coincident absorptions, for example at 5600 km s-1, and some anticoincidences, such as the strong one at $`10000`$ km s-1, while at smaller separations there are only coincident features (Figure 7).

To quantify the amount of correlation between the spectra of the AGNs of the simulated pairs we use the ’hits-and-misses’ statistics described by McGill (1990), Crotts et al. (1994), Bechtold et al. (1994), Charlton et al. (1997). A coincidence is defined as the case in which an absorption line is present in both the spectra within a given velocity difference $`\mathrm{\Delta }v`$ and above some signal-to-noise ratio. An anticoincidence is defined when a line is present in one spectrum but not in the other one. If there are two lines within $`\mathrm{\Delta }v`$, a pretty rare event given the incidence of OVII and OVIII absorbers, we count only one coincidence and no anticoincidence, as in Fang et al. (1996). We also generate a set of spectra uncorrelated in the transverse direction and compute the same statistics to take into account the level of random coincidences. Despite the arbitrariness of the definition, which involves the choice of a velocity binning for the spectra, it has been shown that it is still useful to define a characteristic size for the observed Ly$`\alpha `$ absorbers (Charlton et al. 1997; Fang et al. 1996). The inferred ’characteristic size’ of the Ly$`\alpha `$ clouds is in agreement with the results of other hydro-dynamical simulations (Theuns et al. 1998) and with semi-analytical models of the Ly$`\alpha `$ forest (Viel et al. 2002a, Viel et al. 2002b).

Results are shown in Figure 9 where we plot the coefficient $`f_{co}`$, the fraction of coincidences, which is defined as the ratio between the number of ’hits’ and the total number of ’hits’ and ’misses’: the larger the value of this coefficient the larger the correlation between the two spectra. These plots have been obtained by simulating absorption spectra of 20 AGN pairs at a median redshift of $`z0.05`$ for 12 different separations, spanning the range $`0.12`$ comoving $`h^1`$ Mpc, that is angular separations from 3 to 70 arcmin. OVII and OVIII absorbers are represented by empty and filled cirlces, respectively. This coefficient has been computed only for absorbers with column densities $`>10^{13}\mathrm{cm}^2`$, which will be detectable by XEUS and Constellation-X. We have decided to show this result selecting all the absorbers detectable by future X-ray missions, i.e. all the absorbers with column densities larger than $`10^{13}\mathrm{cm}^2`$. This assumption is not crucial: the values assumed by the coefficient $`f_{co}`$ are pretty robust and change significantly only if we include the absorbers with lower column densities. In fact, for these absorbers the level of random coincidences is higher and determines higher values of $`f_{co}`$. The values assumed by $`f_{co}`$ are mainly determined by the absorbers with column densities in the range $`10^{13.5}10^{15}`$ cm-2. For larger column densities $`f_{co}`$ becomes smaller because the number of random coincidences found is smaller.

In the left panel a bin in velocity of 100 km s-1 has been chosen for the hits-and-misses statistics, while in the right one the bin is of 200 km s-1. The dashed line is the level of random coincidences calculated from the ’hits-and-misses’ analysis of fully uncorrelated AGN spectra. For both absorbers the coefficient peaks within $`0.2`$ comoving $`h^1`$ Mpc and then drops. In both panels the value at which the coefficient becomes comparable to the random coincidence level is at $`1.2`$ comoving $`h^1`$ Mpc. Of course the number of random coincidences depends on the binning chosen, being smaller for smaller $`\mathrm{\Delta }v`$. However, we stress that the number of random coincidences is very small, given the small number of systems expected per unit redshift path.

Assuming that the noise-level for $`f_{co}`$ is of the order of 0.05, as we can infer from the right panel of Figure 9, a 3$`\sigma `$ detection of this filamentary structure will require at least 20 pairs with separations $`<\text{ }1`$ comoving $`h^1`$ Mpc corresponding to the range $`<\text{ }20`$ arcmin if the absorbers are placed at $`z<\text{ }0.1`$. Another conclusion is that a sample of at least 20 AGN pairs with this separation is necessary to discriminate between random coincidences and ’real’ signal from this filamentary structure.

### 4.1 XEUS and Constellation-X detectability of OVII absorbers

In this Section we quantify the detectability of OVII absorbers, following the treatment and formalism of Sarazin (1989). In particular, we compute detection thresholds of the filamentary network for single AGN spectra and for AGN pairs. This means that in the latter case we will need 20 pairs of sources with separation $`<\text{ }20`$ arcmin, in a given portion of the sky, to probe the characteristic size of the absorbers by analysing coincident absorption features (see Section 4). We decided to focus only on OVII absorbers. As stressed above OVII and OVIII absorbers probe the same structures, but, for a given column density the corresponding OVII absorber equivalent width is a factor 2 higher than the OVIII one. This means that OVII absorbers can probably be detected in slightly lower density regions than OVIII ones, so they can better trace the bulk of the filamentary structure.

In the following discussion all the fluxes are expressed in the ROSAT All-Sky Survey (RASS) 0.1-2.4 keV energy band.

In order to detect absorption lines, with a specified signal-to-noise level $`S/N`$, the detector must collect, within the instrument energy resolution $`\mathrm{\Delta }E`$, a number of continuum photons given by $`N_{photons}>\text{ }\left(\frac{S}{N}\right)^2\left(\frac{W_{abs}}{\mathrm{\Delta }E}\right)^2`$, where $`W_{abs}`$ is the absorber equivalent width. The corresponding flux at the line energy E is:

$$F>\text{ }1.3\times 10^{12}\left(\frac{E}{1\mathrm{keV}}\right)^1\left(\frac{S/N}{3}\right)^2\left(\frac{\mathrm{\Delta }E}{10\mathrm{eV}}\right)\left(\frac{W_{abs}}{100\mathrm{km}/\mathrm{s}}\right)^2\left(\frac{AQ}{10000\mathrm{cm}^2}\right)^1\left(\frac{t}{100\mathrm{ksec}}\right)^1\mathrm{erg}\mathrm{cm}^2\mathrm{s}^1\mathrm{keV}^1,$$
(5)
where $`AQ`$ is the effective collecting area ($`Q`$ is the quantum efficiency of the detector) and $`t`$ is the integration time.

We estimate the source flux needed to detect OVII absorption lines ($`E=0.57`$ keV, in the observer rest frame) using different X-ray telescopes and exposure times of 500 ksec. Although a conservative line detection requires a $`S/N=5`$, in all our calculations below we choose $`S/N=3`$.

We find that the detection of these absorbers is quite difficult. The resolution and collecting areas at $`0.5`$ keV of the present generation X-ray satellites are: $`R=E/\mathrm{\Delta }E=500`$, $`AQ=25`$ cm2 and $`R=500`$ and $`AQ=70`$ cm2 for Chandra and XMM-Newton, respectively. If we substitute these values in eq. (5) we find that column densities of $`10^{14}`$ cm-2, corresponding to equivalent widths $`W<\text{ }10`$ km s-1, cannot in practice be detected by these satellites: they either require exceptionally bright background sources with fluxes $`10^9\mathrm{erg}\mathrm{cm}^2\mathrm{s}^1`$ (calculated assuming a source spectrum F $`E^{0.7}`$) or unrealistic exposure time of the order of the Chandra life-time. On the other hand, if we consider more realistic, although still quite bright and rare, X-ray sources with $`F10^{11}\mathrm{erg}\mathrm{cm}^2\mathrm{s}^1`$ we find that the only observable absorption features detectable by Chandra and XMM-Newton are the ones with $`W>200`$ km/s (see also Mathur et al. 2002; Fang et al. 2002b).

Conversely, the low density network of filaments can be detectable by next generation satellites, with much higher collecting area, such as Constellation-X 111http://constellation.gsfc.nasa.gov and especially XEUS 222http://astro.estec.esa.nl/SA-general/Projects/XEUS. For these satellites the effective collecting area and resolution are: $`R=1500`$, $`AQ=3000`$ cm2 for Constellation-X and $`R=800`$ and $`AQ=40000`$ cm2 for XEUS. We mention here that the value $`R=1500`$ for constellation-X is that expected from the off-plane option for the Reflection Grating Spectrometer (W. Cash 2002, http://constellation.gsfc.nasa.gov) and that for XEUS the ultimate satellite configuration XEUS-2 will allow to reach an effective area of $`AQ=2\times 10^5`$ cm2.

Assuming that sources with a flux $`10^{11}\mathrm{erg}\mathrm{cm}^2\mathrm{s}^1`$ are observed, we find that the minimum detectable equivalent width is $`6`$ km s-1 and $`3`$ km s-1, for Constellation-X and XEUS, respectively. With XEUS-2 the minimum detectable equivalent width will be $`1.5`$ km s-1. These equivalent widths will probe the range of absorbers of column densities of $`>\text{ }10^{14.6}`$ cm-2 for Constellation-X, $`>\text{ }10^{13.8}`$ cm-2 for XEUS and $`>\text{ }10^{13.5}`$ cm-2 for XEUS-2.

We now estimate the number of AGN with fluxes sufficient to allow detection of absorbers of given $`W_{abs}`$. We use the fit to the log N - log S relation $`𝒩(>S)5.1\times 10^{16}(S/\mathrm{erg}\mathrm{cm}^2`$ s$`{}_{}{}^{1})^{1.5}\mathrm{sr}^1`$, originally calculated in the Einstein band and here renormalized to the ROSAT band (see Maccacaro et al. 1982; Hasinger et al. 1993). Formally this fit is valid for fluxes in the range $`10^{13}10^{11}\mathrm{erg}\mathrm{cm}^2\mathrm{s}^1`$ and extrapolation to lower fluxes can be dangerous as we know, for example, that the ROSAT Log N - Log S function in the Lockman Hole has a shallower slope than the one used here at lower fluxes (Hasinger et al. 1998). Nevertheless, as we are interested in very bright sources, the Maccacaro et al. fit should be a good approximation. Assuming that the typical quasar spectrum is given by $`FE^{0.7}`$ from eq. (5) and the adopted Log N - Log S relation we find:

$$𝒩(>F)100\left(\frac{E}{1\mathrm{keV}}\right)^{0.45}\left(\frac{S/N}{3}\right)^3\left(\frac{\mathrm{\Delta }E}{10\mathrm{eV}}\right)^{1.5}\left(\frac{W_{abs}}{100\mathrm{km}/\mathrm{s}}\right)^3\left(\frac{AQ}{10000\mathrm{cm}^2}\right)^{1.5}\left(\frac{t}{100\mathrm{ksec}}\right)^{1.5}\mathrm{sr}^1.$$
(6)
In Figure 10 we plot the number of sources per steradian with enough flux to allow the detection of absorbers of a given equivalent width. The solid line refers to XEUS, the dotted one to XEUS-2 and the dashed one on the right refers to Constellation-X. In the same Figure we also report the corresponding OVII column densities.

From Figure 10 one can see that XEUS can detect roughly 10 sources per steradian with equivalent width $`2.5`$ km/s, corresponding to column densities of $`10^{13.8}`$ cm-2 for OVII. For these sources, at least a number of 30 features can be seen in absorption per unit redshift range if the current models for the WHIM are correct (Figure 8). XEUS-2 will find a factor of 10 more sources per steradian than XEUS, which will allow to probe the same WHIM structures. Constellation-X will not be able to probe the same very low density range. For example, in this case 10 single sources per unit solid angle, bright enough to probe absorbers with an equivalent width of $`10`$ km s-1, are predicted to be in the sky. The number of absorption features with this equivalent width is approximately 1 per unit redshift range.

We estimate now the number of AGN pairs with flux above a given value and whose relative separation is between $`\theta _{min}`$ and $`\theta _{max}`$. To this aim we consider a survey area $`\mathrm{\Sigma }`$ and we use the X-ray sources angular correlation function $`w(\theta )=(\theta /\theta _0)^{0.8}`$, with $`\theta _03`$ arcmin (Akylas et al. 2000).

If we indicate with $`𝒩_\mathrm{\Sigma }=_\mathrm{\Sigma }𝒩(>F)d\mathrm{\Omega }`$ the total number of sources in the survey with flux above a given value, then the number of distinct pairs whose members have a relative separation between $`\theta _{min}`$ and $`\theta _{max}`$ is:

$$𝒩_{pairs}(\theta _{min},\theta _{max})=\pi 𝒩𝒩_\mathrm{\Sigma }_{\theta _{min}}^{\theta _{max}}(1+w(\theta ))\mathrm{sin}\theta d\theta \frac{\pi 𝒩𝒩_\mathrm{\Sigma }}{2}[\theta _{max}^2+1.7\theta _0^{0.8}\theta _{max}^{1.2}],$$
(7)
where the last expression is valid for small separations and for $`\theta _{min}=0`$.

We choose the limits $`\theta _{min}=0`$ and $`\theta _{max}=20`$ arcmin, in computing the number of pairs. For these angular separations, as we have found in Section 4, we expect to find significant coincidence features in the spectra of AGN pairs.

The vertical bars in Figure 10 represent the minimum absorber equivalent width which can be probed by using AGN pairs, instead of single AGN spectra. Detection of the filamentary structure with future missions can be achieved by finding coincident absorptions in AGN pairs, as we showed in Section 4.

The first threshold on the left in Figure 10 corresponds to the equivalent width which can be probed from a sample of 20 pairs, with the right separations, found in the whole sky, while the threshold to the right corresponds to a field of 4 deg2. This means that, if XEUS or Constellation-X will observe a random field of 4 deg2 in the sky, they will find 20 AGN pairs with the wanted separation with fluxes high enough to probe the WHIM structures at the corresponding equivalent width. If this is the case, the expected number of OVII absorption features seen by XEUS per unit redshift range is $`2.5`$. Thus, if these objects are placed at $`z0.5`$ we predict that a sample of 20 spectra of AGN pairs will possibly show a total of $`50`$ absorption OVII features, while Constellation-X will detect a factor of 3 less absorption features. The more powerful XEUS-2 will detect $`100`$ absorptions in the same conditions. An observation of the whole sky will possibly allow these satellites to detect WHIM structures at smaller column densities.

To conclude we verify if the ROSAT all-sky survey bright source catalogue (Voges et al. 1999) may contain the positions of a sufficient number of pairs to probe the WHIM structure. This catalogue contains a sample of 18,811 sources down to a limiting ROSAT PSPC count- rate of 0.05 cts s-1 in the 0.1 - 2.4 keV band, which corresponds to a flux limit $`5.2\times 10^{13}`$ erg cm-2 s-1 (assuming a photon spectrum which can be modeled by a power law with slope -1 and no absorption, 1 cts/s $`1\times 10^{11}`$ erg cm-2 s-1). We find that the number of pairs of sources whose separation is smaller than 20 arcmin is $`100`$ for fluxes larger than $`3.3\times 10^{12}`$ erg cm-2 s-1, and is $`300`$ for fluxes larger than $`2.2\times 10^{12}`$ erg cm2 s-1. If we consider that approximately 16 % of the X-ray sources have unique counterpart of extragalactic origin (Voges et al. 1999), we deduce that the RASS catalogue should contain the positions of $`48`$ pairs of extra-galactic sources with fluxes $`>\text{ }2\times 10^{12}`$ erg cm -2 s-1. These bright sources will allow XEUS to probe absorbers of $`W7`$ km s-1 (see eq. (5)).

## 5 Discussion and conclusions

In this paper we have presented a semi-analytical technique to simulate X-ray absorption spectra of distant AGNs. The Warm Hot Intergalactic Medium (WHIM) can be detected both in absorption and in emission. Recently, a number of missions devoted to the study of the WHIM through the observation of emission lines (see e.g. the missions SPIDR 333http://www.bu.edu/spidr/ or IMBOSS 444http://www.ias.rm.cnr.it/IMXS/iasimxs.html) have been proposed. Unlike the line emission process, which is proportional to the density square, the line absorption has the advantage of being proportional to the density of the intervening absorber. Furthermore, absorption measurements are not subjected to background contamination. Absorption measurements, then, happen to be more sensitive than emission ones to the lower density regions of the WHIM. Having a larger volume filling factor, these low regions constitute the bulk of the filamentary structure. Thus, absorption measurements, besides representing an important tool complementary to emission ones, are probably more promising to study the large scale structures traced by the WHIM.

The method proposed here uses as input four results obtained from the $`z=0`$ output of a large box-size $`\mathrm{\Lambda }`$CDM hydro-dynamical simulations: the probability distribution function of the gas, the gas temperature-density relation, the gas metallicity-density relation. These quantities are sufficient to give a good description of the number of OVII and OVIII absorbers expected in the low redshift Universe. Simulations of these two ions are very important: they are the strongest trace of the WHIM. If we sum up all the baryons at $`z2`$ and compare this value with the baryons seen at $`z0`$, we discover that a significant fraction of the baryons is missing (Cen & Ostriker 1999). This fraction is thought to be in the WHIM. Thus, the detectability of this network of OVII and OVIII filaments is of fundamental importance.

We have extended the simulation of a line-of-sight (LOS) to the simulation of pairs of LOS with the mathematical implementation reported in Appendix A1 (see also Viel et al. 2002a). This has been done in order to give an estimate of the characteristic size of the absorbers by using the information contained in the transverse direction. While single LOS are sensitive to the size of single absorbers, pairs of LOS should be more reliable a probe of the network of filaments traced by the metals as they are more sensitive to the random orientation of each filament. We have performed the so called ‘hits-and-misses’ statistics (McGill 1990; Charlton et al. 1997) on a sample of simulated pairs to obtain an estimate of this size.

The main result of the first part of the paper can be summarized as follows: i) semi-analytical models successfully reproduce the number of absorbers per unit redshift path of hydro-dynamical simulations and we expect $`>\text{ }30`$ absorbers per unit redshift with an OVII column density $`>10^{13.5}`$ cm-2; ii) scatter in the metallicity-gas density relation plays a significant role in increasing the number of absorbers with column densities $`>\text{ }10^{14}`$ cm-2; iii) the characteristic size of these absorbers at $`z0`$ obtained by analysing coincident absorptions in simulated spectra of AGN pairs, is $`1`$.

In the second part we have discussed the detectability of the absorbers by present X-ray missions, such as Chandra and XMM-Newton and future ones, such as XEUS and Constellation-X. Up to now only few ‘trees’ of the ‘X-ray forest’ have been detected by Chandra (Fang et al. 2002b, Nicastro et al. 2002a). Most likely, these absorbers sample the high column density tail of the distribution shown in Figure 6. A detection of the lower column density absorbers, whose volume filling factor is larger, can be achieved only by future X-ray missions whose collecting area and sensitivity will be higher.

We have found that at least a background of 20 pairs of bright extragalactic AGN-like sources at a median redshift of $`z0.5`$ and with relative separations lower than 20 arcmin is needed to probe this filamentary structures with ‘hits-and-misses’ statistics. We have made a theoretical estimate to address the detectability of these sources by XEUS and Constellation-X taking into account the bright X-ray sources angular correlation function (Akylas et al. 2000). The main conclusions are: i) by observing a random field of 4 deg2 XEUS and Constellation-X will be able to probe WHIM absorbers, with ‘hits-and-misses’ statistics, which will correspond to OVII column densities $`10^{14.8}`$ cm-2 and $`10^{15.5}`$ cm-2, respectively; ii) in this field XEUS and XEUS-2 will possibly find $`50`$ and $`80`$ absorption features, respectively, if we assume that the sources are placed at a median redshift of $`z0.5`$; iii) XEUS will detect at least 10 sources per steradian bright enough to probe structures with an absorber equivalent width of 2.5 km/s, corresponding to an OVII column density of $`10^{13.8}`$ cm-2; iv) Constellation-X will detect few thousand sources per steradian bright enough to probe structure with an equivalent width of $`100`$ km/s. We have checked in the ROSAT all-sky survey bright source catalogue (Voges et al. 1999) to determine if these sources had already been detected. We have found that the required number of extragalactic sources with the necessary angular separation is present in the catalogue. Re-observation of these sources will allow future X-ray missions to detect baryons in the low redshift Universe and to unveil the filamentary structure where they preferentially reside.

## Acknowledgements.

We thank M. Haehnelt, V. Mainieri, G. Matt, M. Mendez, G. Perola, C. Porciani for useful discussions and technical help. We thank the anonymous referee for helpful comments and X. Chen for kindly providing his data. MV acknowledges partial financial support from an EARA Marie Curie Fellowship under contract HPMT-CT-2000-00132. This work has been partially supported by the European Community Research and Training Network ‘The Physics of the Intergalactic Medium’ and by grant ASC97-40300.

## Appendix A Correlation procedure to generate Multiple lines of sight density and velocity fields

We briefly report here the correlation procedure of Viel et al. (2002a).

Let us start by obtaining auto and cross-spectra for the 1D random fields which are defined along single or multiple LOS. In the Gaussian case these quantities fully determine the statistical properties of the generated fields. If $`\psi (𝐱)`$ is a 3D random field with Fourier transform $`\psi (𝐤)`$ and 3D power spectrum $`P(|k|)`$, one can define the LOS random field $`\psi _{}(x_{},𝐱_{})`$ as the 1D Fourier transform

$$\psi _{}(k_{}|𝐱_{})\frac{d^2k_{}}{(2\pi )^2}e^{i𝐤_{}𝐱_{}}\psi (k_{},𝐤_{}).$$
(8)
The cross-spectrum $`\pi (|k_{}||r_{})`$ for this LOS random field along parallel LOS, separated by a transverse distance $`r_{}`$, is defined by

$$\psi _{}(k_{}|𝐱_{})\psi _{}(k_{}^{}|𝐱_{}+𝐫_{})=2\pi \delta _D(k_{}+k_{}^{})\pi (|k_{}||r_{}),$$
(9)
where $`\delta _\mathrm{D}`$ is the Dirac delta function and $`\pi (k|r_{})`$ can be related to the 3D power spectrum as follows

$$\pi (k|r_{})=\frac{d^2k_{}}{(2\pi )^2}e^{i𝐤_{}𝐫_{}}P(\sqrt{k_{}^2+k^2}).$$
(10)
Integrating over angles and shifting the integration variable yields

$$\pi (k|r_{})=\frac{1}{2\pi }_k^{\mathrm{}}𝑑qqJ_0(r_{}\sqrt{q^2k^2})P(q),$$
(11)
where $`J_0`$ is the Bessel function of order $`0`$.

In the limit of vanishing distance between the two LOS $`J_01`$ and the above formula reduces to the standard relation for the LOS (1D) auto-spectrum in terms of the 3D power spectrum (Lumsden et al. 1989)

$$p(k)\pi (k|r_{}=0)=\frac{1}{2\pi }_k^{\mathrm{}}𝑑qqP(q).$$
(12)
Our next problem is how to generate the two random fields $`\delta ^{\mathrm{IGM}}`$ and $`v^{\mathrm{IGM}}`$ in 1D Fourier space ($`\mathrm{}<k_{}<\mathrm{}`$). These random fields have non-vanishing cross-correlations but unlike their 3D Fourier space counterparts they cannot be related by simple algebraic transformations.

We gan generally write any M-dimensional Gaussian random vector $`𝐕`$ with correlation matrix $`𝐂`$ and components $`c_{ij}=V_iV_j`$, as a linear combination of another M-dimensional Gaussian random vector $`𝐗`$ with diagonal correlation matrix, which we can take as the identity $`𝐈`$ without any loss of generality. The transformation involves the M $`\times `$ M matrix $`𝐀`$, with components $`\alpha _{ij}`$, as follows: $`𝐕=\mathrm{𝐀𝐗}`$. One gets $`𝐂=\mathrm{𝐀𝐀}^T`$, i.e. $`c_{ij}=_k\alpha _{ik}\alpha _{jk}`$. There is a slight complication because $`𝐕`$ is a random vector $`\mathrm{𝑓𝑖𝑒𝑙𝑑}`$ defined in 1D Fourier space. We can, however, extend the above formalism to vector fields, assuming that $`𝐗`$ is a Gaussian vector field with white-noise power spectrum,

$$X_i(k_{})X_j(k_{}^{})=2\pi \delta _{ij}\delta _D(k_{}+k_{}^{})$$
(13)
($`\delta _{ij}`$ is the Kronecker symbol). We then have $`V_i=_{j=1}^2\alpha _{ij}X_j`$ where 3 of the 4 $`\alpha _{ij}`$ components are determined by the conditions $`_{k=1}^2\alpha _{ik}\alpha _{jk}=p_{ij}`$. The remaining freedom (due to the symmetry of the original correlation matrix) can be used to simplify the calculations.

It is straightforward to extend our formalism to simulate the IGM properties along parallel LOS. Let $`𝐕(k_{})`$ and $`𝐖(k_{})`$ be two 1D Gaussian random vector fields obtained as in Section 2, each with the same set of coefficients $`\alpha _{ij}`$ but starting from two independent white-noise vector fields $`𝐗`$ and $`𝐘`$ (i.e. such that $`X_iY_j=0`$). Then both $`𝐕`$ and $`𝐖`$ have the correct LOS auto-spectra by construction while their mutual cross-spectra vanish: $`V_i(k_{})W_j(k_{}^{})=0`$.

Let us further define a new vector $`𝐕^{}(k_{}|r_{})`$ with components $`V_i^{}=_{k=1}^2\left(\beta _{ik}V_k+\gamma _{ik}W_k\right)`$, such that its auto and cross-spectra components are given by,

$`V_i(k_{})V_j(k_{}^{})`$ $`=`$ $`V_i^{}(k_{}|r_{})V_j^{}(k_{}^{}|r_{})=2\pi \delta _D(k_{}+k_{}^{})p_{ij}(|k_{}|),`$
$`V_i(k_{})V_j^{}(k_{}^{}|r_{})`$ $`=`$ $`2\pi \delta _D(k_{}+k_{}^{})\pi _{ij}(|k_{}||r_{}).`$ (14)
The vectors $`𝐕`$ and $`𝐕^{}`$ will then represent our physical IGM linear fields on parallel LOS at a distance $`r_{}`$. They will be statistically indistinguishable from those obtained by drawing two parallel LOS separated by $`r_{}`$ in a 3D realization of the linear IGM density and velocity fields.

The transformation coefficients are determined by the equations $`_{k,\mathrm{}=1}^2\left(\beta _{ik}\beta _j\mathrm{}+\gamma _{ik}\gamma _j\mathrm{}\right)p_k\mathrm{}=p_{ij}`$ and $`_{k=1}^2\beta _{ik}p_{kj}=\pi _{ij}`$. A complete set of coefficients can be found in Viel et al. (2002a).