# Soft X-ray Absorption by High-Redshift Intergalactic Helium

## 1 Introduction

The Ly$`\alpha `$ absorption by intergalactic He ii has now been observed in four quasars at redshifts $`2.4<z<3.2`$ (Jakobsen et al. 1994; Davidsen et al. 1996; Reimers et al. 1997; Anderson et al. 1999; Heap et al. 1999). While it is known that hydrogen (together with He i) was ionized at redshift $`z5`$, given the presence of transmitted flux to the blue of the hydrogen Ly$`\alpha `$ wavelength in sources up to that redshift, the double ionization of helium can take place at a lower redshift. The reason is simple: even though there are only $`8`$ helium atoms for every 100 hydrogen atoms, He ii recombines at a rate $`5.5`$ times faster than hydrogen. Consequently, if a large number of recombinations for each helium ion take place, then the He ii reionization will be delayed relative to the hydrogen reionization as long as the number of photons emitted to the intergalactic medium (hereafter IGM) above $`54.4`$ eV is less than 0.44 times the number of photons emitted above $`13.6`$ eV. This condition is satisfied by the known sources of ionizing radiation (quasars and galaxies). In practice, the mean number of recombinations of He ii may not be very large: for the uniform IGM, the He ii recombination rate is $`4`$ times the Hubble rate at $`z=4`$, and then a somewhat softer emission spectrum is required for a delayed He ii reionization. However, the recombination rate is enhanced by the clumping factor of the ionized gas, which increases as progressively denser gas is reionized, as described recently in Miralda-Escudé, Haehnelt, & Rees (1999; hereafter MHR).

The present state of the observations of intergalactic He ii can be summarized as follows: at redshift $`z2.8`$, there is a He ii “Ly$`\alpha `$ forest”, with a flux decrement that is consistent with the hydrogen Ly$`\alpha `$ forest flux decrement, and a background spectrum produced by quasars, plus a possible contribution from galaxies to increase the ratio $`J_{\mathrm{H}\mathrm{i}}/J_{\mathrm{He}\mathrm{ii}}`$ of intensities at the H i and He ii ionization edges (Davidsen et al. 1996, Miralda-Escudé et al. 1996, Bi & Davidsen 1997, Croft et al. 1997). At $`z2.8`$, the He ii Ly$`\alpha `$ spectra appear to be divided into two different types of regions: those where transmitted flux is still observed (with a similar $`J_{\mathrm{H}\mathrm{i}}/J_{\mathrm{He}\mathrm{ii}}`$ implied for the ionizing background), and those where the transmitted flux is undetectable or, at least, much smaller (see Heap et al. 1999 and references therein). The regions with greater transmitted flux occupy a small fraction of the spectrum at $`z3`$, and they have typical widths of $`10^3\mathrm{km}\mathrm{s}^1`$. Some of them can be attributed to the proximity effect of the source being observed (e.g., Hogan et al. 1997).

These observations are not yet conclusive in clarifying the ionization history of He ii in the IGM. The reason is that the optical depth of a homogeneous medium with the mean baryon density of the universe, where all the helium is He ii, is very large:

$$\tau _{0,\mathrm{He}\mathrm{ii}}=1.7\times 10^3\frac{\mathrm{\Omega }_bhY}{0.007}\frac{H_0(1+z)^{3/2}}{H(z)}\left(\frac{1+z}{4}\right)^{3/2}.$$
(1)   
Therefore, if only a fraction $`10^3`$ of the average IGM helium density is in the form of He ii, the flux transmission can already be reduced to very low levels. Even when the line of sight crosses a void with density $`\rho =0.1\overline{\rho }`$ (about the lowest densities in the photoionized IGM at $`z=3`$ according to numerical simulations of Cold Dark Matter models; see MHR and references therein), the void will still be optically thick to He ii Ly$`\alpha `$ photons if its He ii fraction is greater than $`0.01`$. The implication is that even a very stringent upper limit to the He ii Ly$`\alpha `$ flux transmission does not imply that most of the He ii had not yet been reionized at the observed redshift.

Although the observations quoted earlier could be simply interpreted as revealing an IGM with a patchy ionization of He ii, with He iii regions surrounding luminous quasars and a pure He ii IGM filling the volume between the He iii regions, the correct picture is probably much more complicated. In the regions where no transmitted flux is detected, the helium might also be doubly ionized over most of the volume, although with a much lower ionizing background intensity than in the regions near luminous quasars. The higher He ii fraction implied in the regions with low intensity, although still much smaller than unity, could be enough to completely absorb the Ly$`\alpha `$ photons. The presence of this low intensity, but widespread, background above the He ii edge is quite plausible, as discussed in MHR, if there is a modest contribution to the emissivity from sources of luminosity much lower than quasars.

A more direct observational probe to the mean fraction of He ii in the IGM at different redshifts would be highly valuable to test models for the reionization history, which is in turn important for the history of galaxy formation owing to the strong heating of the IGM gas produced by the He ii reionization (e.g., Efstathiou 1992; Miralda-Escudé & Rees 1994). It is shown in this paper that this new probe may be found in the soft X-ray continuum absorption spectra of high-redshift quasars.

The distributed He ii in the IGM should cause, in addition to the He ii forest absorption in Ly$`\alpha `$ and the other Lyman series lines, the continuum absorption (due to photoionization) below a rest-frame wavelength $`228`$ Å. This continuum absorption should be very large near the edge, but the flux should recover at shorter wavelengths. We shall see that when most of the helium in the IGM is He ii, this recovery of the flux should occur at $`0.5\mathrm{keV}`$, or an observed frequency of $`0.1\mathrm{keV}`$ at $`z=4`$.

## 2 X-Ray Absorption by Intergalactic Helium

Let $`F(z)`$ be the mean fraction of helium in the IGM in the form of He ii as a function of redshift. For the moment, we consider a homogeneous IGM (with a uniform fraction $`F(z)`$), but because we shall consider continuum absorption only, the results will be valid also for a clumpy IGM, as long as the clumps have a large covering factor. The continuum absorption spectrum by He ii on a source at redshift $`z_s`$ is given by

$$\tau (\nu )=_0^{z_s}𝑑z\frac{dl}{dz}n_{\mathrm{He},0}(1+z)^3F(z)\sigma _{\mathrm{He}}[\nu (1+z)],$$
(2)   
where $`dl`$ is the proper length element along the line of sight, $`n_{\mathrm{He},0}`$ is the mean primordial helium density extrapolated to the present epoch, $`\sigma _{\mathrm{He}}(\nu )`$ is the photoionization cross section of He ii, and $`\nu `$ is the observed frequency. If He ii is highly ionized at low redshift, most of the contribution to $`\tau (\nu )`$ will be from high redshift, where we can use $`dl/dzcH_0^1\mathrm{\Omega }_0^{1/2}(1+z)^{5/2}`$. We also define $`\mathrm{\Sigma }_{\mathrm{He}}(\nu )\sigma _{\mathrm{He}}(\nu )(\nu /\nu _{\mathrm{He}})^3`$, where $`\nu _{\mathrm{He}}`$ is the frequency at the ionization edge of He ii (the function $`\mathrm{\Sigma }_{\mathrm{He}}`$ varies only slowly with frequency and is plotted below in Fig. 1). Equation (1) can be rewritten as

$$\tau (\nu )=\frac{cn_{\mathrm{He},0}}{H_0\mathrm{\Omega }_0^{1/2}}\sigma _{\mathrm{He}}[\nu (1+z_s)]\left(1+z_s\right)^3_0^{z_s}\frac{dzF(z)}{(1+z)^{5/2}}\frac{\mathrm{\Sigma }_{\mathrm{He}}[\nu (1+z)]}{\mathrm{\Sigma }_{\mathrm{He}}[\nu (1+z_s)]},$$
(3)   
or, using $`\mathrm{\Omega }_0=0.3`$, $`\mathrm{\Omega }_bh^2=0.019`$, $`Y=0.25`$, $`H_0=65\mathrm{km}\mathrm{s}^1\mathrm{Mpc}^1`$:

$$\tau (\nu )=550\left(\frac{\nu _{\mathrm{He}}}{\nu }\right)^3\frac{\mathrm{\Sigma }_{\mathrm{He}}[\nu (1+z_s)]}{\mathrm{\Sigma }_{\mathrm{He}}(\nu _{\mathrm{He}})}_0^{z_s}\frac{dzF(z)}{(1+z)^{5/2}}\frac{\mathrm{\Sigma }_{\mathrm{He}}[\nu (1+z)]}{\mathrm{\Sigma }_{\mathrm{He}}[\nu (1+z_s)]}.$$
(4)   
As an example, we consider the simple case of a sudden and complete reionization of He ii at $`z=z_i`$: $`F=1`$ at $`z>z_i`$ and $`F=0`$ at $`z<z_i`$. The integral in the above equation is then well approximated as $`(z_sz_i)/(1+z_s)^{5/2}`$ (for small $`z_sz_i`$). For $`z_i=3`$ and $`z_s=4`$, $`\tau (\nu )10(\nu _{\mathrm{He}}/\nu )^3`$, so the optical depth reaches unity at an observed frequency $`\nu 0.12`$ keV.

## 3 Comparison to Galactic absorption

Any extragalactic source will also be absorbed in soft X-rays by Galactic hydrogen and helium, with a hydrogen column density $`N_{\mathrm{H}\mathrm{i}}`$ that can be determined from the H i 21 cm emission. The absorption optical depth is

$$\tau _G(\nu )=N_{\mathrm{H}\mathrm{i}}\left[\sigma _{\mathrm{H}\mathrm{i}}(\nu )+0.084\sigma _{\mathrm{He}\mathrm{i}}(\nu )\right]N_{\mathrm{H}\mathrm{i}}\sigma _{\mathrm{H}\mathrm{i}}(\nu )R_{\mathrm{He}}(\nu ),$$
(5)   
where the ratio of helium to hydrogen atoms is $`(Y/3.97)/(1Y)=0.084`$, $`\sigma _{\mathrm{H}\mathrm{i}}`$ and $`\sigma _{\mathrm{He}\mathrm{i}}`$ are the photoionization cross sections of H i and He i, and we have defined the quantity $`R_{\mathrm{He}}(\nu )`$. At the frequencies that will be of interest to us, $`h\nu 0.2`$ keV, this quantity is $`R_{\mathrm{He}}(\nu )3`$ (this is shown in detail below, in Fig. 1).

Dividing the optical depth due to the high-redshift IGM from (4) by the Galactic optical depth in (5), we obtain the ratio to be:

$$\frac{\tau (\nu )}{\tau _G(\nu )}=\frac{56}{N_{\mathrm{H}\mathrm{i},20}R_{\mathrm{He}}(\nu )}\frac{\mathrm{\Sigma }_{\mathrm{He}}[\nu (1+z_s)]}{\mathrm{\Sigma }_{\mathrm{He}}(4\nu )}_0^{z_s}\frac{dzF(z)}{(1+z)^{5/2}}\frac{\mathrm{\Sigma }_{\mathrm{He}}[\nu (1+z)]}{\mathrm{\Sigma }_{\mathrm{He}}[\nu (1+z_s)]},$$
(6)   
where $`N_{\mathrm{H}\mathrm{i},20}=N_{\mathrm{H}\mathrm{i}}/(10^{20}\mathrm{cm}^2)`$. This fiducial H i column density is around the lowest value that is usually reached at high Galactic latitude (e.g., Laor et al. 1997). Using the same example as in §2, where the integral is given by $`(z_sz_i)/(1+z_s)^{5/2}0.018`$, we find $`\tau (\nu )/\tau _G(\nu )1/[N_{\mathrm{H}\mathrm{i},20}R_{\mathrm{He}}(\nu )]1/(3N_{\mathrm{H}\mathrm{i},20})`$. Thus, we conclude that the absorption by high-redshift intergalactic helium should be smaller than the Galactic absorption only by a factor $`3`$ when the Galactic column density has the lowest value normally found in sources, $`N_{\mathrm{H}\mathrm{i}}10^{20}\mathrm{cm}^2`$.

Laor et al. (1997) measured the Galactic column densities in a sample of 23 low-redshift quasars from the absorption in soft X-rays. Comparing the column densities they derived with those measured from the 21 cm emission, they show that there are no significant differences within observational error. In particular, the two quasars with the smallest error (10%) in the column density determined from the soft X-ray spectrum by Laor et al. also agree with the 21 cm column density (see their Fig. 2). Provided that the Galactic absorption can be subtracted from the measured X-ray absorption on a high-redshift quasar, given the column density measured from 21cm emission, it should be possible to detect the continuum absorption of intergalactic helium in the soft X-rays.

A natural question to ask here is if the shape of the absorption can be used to distinguish Galactic absorption from the high-redshift He ii absorption. Unfortunately, the shape is almost identical in the two cases, as we shall now see. In Figure 1, the solid line is the hydrogen cross section $`\sigma _{\mathrm{H}\mathrm{i}}(\nu )`$, multiplied by $`(\nu /\nu _{\mathrm{H}\mathrm{i}})^3`$, and the dashed line is $`[\sigma _{\mathrm{H}\mathrm{i}}(\nu )+n_{\mathrm{He}}/n_\mathrm{H}\sigma _{\mathrm{He}\mathrm{i}}(\nu )](\nu /\nu _{\mathrm{H}\mathrm{i}})^3`$ \[so the ratio of the dashed line to the solid line is $`R_{\mathrm{He}}(\nu )`$, defined in eq. (5)\]. We have used the exact analytical expression for $`\sigma _{\mathrm{H}\mathrm{i}}(\nu )`$ (e.g., Spitzer 1978), and the fit of Verner et al. (1996) for $`\sigma _{\mathrm{He}\mathrm{i}}(\nu )`$. The dashed line is therefore the expected shape of the Galactic absorption, except for the fact that some of the ionized gas in the Galaxy may not have the helium doubly ionized, and there could then be some additional He ii absorption at $`z=0`$. This is shown by the dotted line, equal to $`\sigma _{\mathrm{He}}(\nu )(n_{\mathrm{He}}/n_\mathrm{H})(\nu /\nu _{\mathrm{H}\mathrm{i}})^3`$, and we see that it has almost the same slope as the dashed line at energies $`0.15`$ keV.

The intergalactic He ii absorption should have the same shape as the solid line, at redshift $`z=3`$ (because the He ii cross section is identical to the H i cross section shifted to a frequency 4 times higher). Notice from equation (4) that the frequency dependence of $`\tau (\nu )`$ is basically the same as the redshifted cross section $`\sigma _{\mathrm{He}}(\nu )`$. Thus, we see from Fig. 1 that the high redshift absorption has a slightly steeper cross section than the Galactic one. In the range from $`0.1`$ keV to $`0.2`$ keV, the Galactic cross section is $`\sigma \nu ^{3.0}`$, and the intergalactic helium cross section is $`\sigma _{\mathrm{He}}\nu ^{3.15}`$. In practice, it will be extremely difficult to reach the high signal-to-noise required to distinguish between these two slopes; moreover, the intrinsic emission spectrum of the source introduces an additional uncertainty. Therefore, the high-redshift helium absorption can probably be detected only as an excess of absorption over that expected from the Galactic column density derived from the 21 cm emission.

## 4 A Reionization Model

In any realistic model, the reionization of all the He ii in the universe will take place over a substantial length of time, of the order of a Hubble time. In the reionization model proposed in MHR, the low-density gas is ionized earlier and reionization advances outside-in. The stage of reionization depends at every redshift on the gas overdensity $`\mathrm{\Delta }_{\mathrm{He}}`$ up to which the helium is mostly ionized to He iii. Using results of numerical simulations for the density distribution and the photon mean free path in the IGM, the model in MHR predicts the He ii Ly$`\alpha `$ flux decrement and the fraction of gas at overdensities greater than $`\mathrm{\Delta }_{\mathrm{He}}`$, assumed to be still in the form of He ii. The observed flux decrement of $`0.7`$ at $`z=2.6`$ (Davidsen et al. 1996) requires $`\mathrm{\Delta }_{\mathrm{He}}70`$, and the fraction of He ii $`F0.2`$ (see Figs. 2 and 9 in MHR). At the slightly higher redshift $`z3.1`$, a flux decrement as high as $`0.99`$ (see Heap et al. 1999) requires $`\mathrm{\Delta }_{\mathrm{He}}12`$, and the fraction of gas at greater overdensities increases only to $`F0.3`$, showing how a small increase in He ii ionization can result in a dramatic decline of the transmitted flux due to the effect that the most underdense voids rapidly become optically thick to the He ii Ly$`\alpha `$ photons.

As an example of a reasonable model matching the above conditions, where the He ii reionization takes place over an extended period of time, we consider the simple case $`F=[(1+z)/7]^{5/2}`$, where the first He ii sources would turn on at $`z=6`$. For this case, equation (6) yields (neglecting the ratios of the function $`\mathrm{\Sigma }_{\mathrm{He}}`$) $`\tau (\nu )/\tau _G(\nu )0.43z_s/N_{\mathrm{H}\mathrm{i},20}/R_{\mathrm{He}}(\nu )`$ (for $`z_s<6`$).

It must be born in mind here that whatever contribution arises from low redshifts to the integral in equation (6) should originate mostly in dense systems, where the helium is self-shielded and remains in the form of He ii (these are the observed Lyman limit and damped absorption systems). When the mean free path between these absorption systems is not much smaller than the Hubble length (i.e., when the absorbers do not have a large covering factor), the He ii column density on a particular line of sight will no longer be equal to the mean value, but will have large fluctuations: it will most often be smaller than the mean value, but it will exceed this mean value when a rare, strong absorption system is intersected. The presence of these rare absorption systems, which should dominate the total He ii content of the universe at low redshift only, can in principle be determined independently for every line of sight from observations of the H i absorption spectrum. In practice, though, the Lyman limit absorption from higher redshift systems, and the He ii Ly$`\alpha `$ absorption itself, will prevent the determination of the H i column density of any low-redshift absorber in a high-redshift source. The absence of strong metal lines can still be a good indication that a strong H i absorber is not present.

In order for a low-redshift absorber to cause soft X-ray absorption comparable to that of the high-redshift IGM, it must have a column density $`N_{\mathrm{H}\mathrm{i}}10^{19.5}(1+z)^3\mathrm{cm}^2`$. Only $``$ 20% of the lines of sight have an absorber with this strength (e.g., Storrie-Lombardi, Irwin, & McMahon 1996). Nevertheless, the fluctuation in the low-redshift absorption will be an additional source of uncertainty (especially if the low-redshift H i Ly$`\alpha `$ absorption cannot be measured), implying that the excess soft X-ray absorption will need to be detected in several sources before one can be sure that the effect of the helium in the low-density IGM at high redshift has been detected.

The fluctuations in the He iicolumn density due to high-redshift absorbers of high column density are less important, because these absorbers do not contain a large fraction of the baryons. For example, at $`z=3`$ we need an absorber with $`N_{\mathrm{H}\mathrm{i}}10^{21}\mathrm{cm}^2`$ to produce an absorption similar to the IGM. These absorbers are rare, and they can also be readily identified in the H i Ly$`\alpha `$ spectrum when their redshift is not much lower than the source redshift.

## 5 Conclusions

Soft X-ray absorption can be a new powerful tool to measure the fraction of helium in the form of He ii as a function of redshift. This direct determination of the He ii fraction provides a straightforward test of any model of reionization based on the observed emitting sources (quasars and galaxies). Although the observations of the He ii Ly$`\alpha `$ absorption spectra provide a much greater wealth of information, their interpretation is more complicated due to the complexities introduced by the highly inhomogeneous IGM and the fluctuating intensity of the ionizing background.

The main challenge in detecting this effect will be to find enough high-redshift, X-ray bright sources, and to accurately measure the 21 cm H i column density to subtract the Galactic contribution to the soft X-ray absorption. A systematic uncertainty that will be faced is the possible existence of some Galactic He ii along the line of sight, which can also produce excess soft X-ray absorption not accounted for by the Galactic H i column density. In addition, extragalactic damped Ly$`\alpha `$ absorbers at low redshifts can also make a significant contribution, which may be difficult to correct for in high-redshift sources. These systematic uncertainties can be put under control once the effect is measured in many sources over a wide range of redshift.

I am grateful to Andy Fabian for encouraging discussions on the possibilities for observing the effect proposed here.