eaa.iop.org DOI: 10.1888/0333750888/2140

Lyman Alpha Forest Michael Rauch

From Encyclopedia of Astronomy & Astrophysics P. Murdin

© IOP Publishing Ltd 2006

ISBN: 0333750888

Institute of Physics Publishing Bristol and Philadelphia

Downloaded on Thu Mar 02 23:49:43 GMT 2006 [131.215.103.76]

Terms and Conditions Lyman Alpha Forest E NCYCLOPEDIA OF A STRONOMY AND A STROPHYSICS

Lyman Alpha Forest Lyα forest systems in this restricted definition, but we have theoretical reasons to believe that there is a genuine The Lyman alpha forest is an absorption phenomenon seen dichotomy between intergalactic gas (represented by the in the spectra of high redshift QSOs and galaxies (figure Lyα forest), even if partly polluted by metals, and the 1). It is the only direct observational evidence we have of invariably metal-enriched higher column density systems, the existence and properties of the general INTERGALACTIC thought to be related to galaxies. MEDIUM, and, as we have reason to believe, of most of the Lyα forest absorption systems have now been baryonic matter contents of the universe. observed from REDSHIFT zero (with UV satellites) up to the On its way to us the light of a bright, distant highest redshifts at which background light sources (QSOs QSO passes through intervening intergalactic gas and and galaxies) can still be found (currently z ∼ 5–6). through gas clouds associated with foreground galaxies. Absorption by the gas modifies the spectra of the The Gunn–Peterson effect: where does the background objects and imprints a record of the gas absorption come from? α clouds’ physical and chemical states on the observed Ly forest absorption in a QSO spectrum was predicted background QSO and galaxy spectra. The whole and first detected by Gunn and Peterson (1965). The arrangement is reminiscent of a giant cosmic slide basic idea is as follows: going back in time an increasing projector, where a QSO plays the role of the light bulb, fraction of the total baryonic mass of the universe must and the intervening gas clouds are the slides, changing the be in the form of gas. The absorption cross-section of α colors of the light source by absorbing parts of the (white) the Ly line of neutral is large enough that spectrum. even if only a small fraction of the total mass of the α The name ‘Lyman α forest’ refers to the appearance universe were in the form of H I the redshifted Ly of the optical QSO spectra, which show a forest of lines should completely absorb a part of the spectrum hundreds of sharp absorption lines, mostly from the of any background light source. The absorption should neutral hydrogen (H I) Lyman α line, superimposed on the essentially assume the shape of an absorption trough α more smoothly varying QSO continuum (figure 2). Almost in a QSO spectrum, extending blueward from the Ly all the lines in the Lyman α forest correspond to the same emission of the QSO. This particular absorption pattern is atomic transition (which at 1215.67 Å is in the referred to in the literature as the ‘Gunn–Peterson effect’. region). The phenomenon was first observed Gunn and Peterson did detect such a trough but the light in the optical waveband (∼ 4000–9000Å) implying that the of the QSO was not completely absorbed and there was gas clouds causing the absorption are highly redshifted by some residual light left in the spectral region in question. the Hubble expansion. The absorption systems appear The relative weakness of the absorption could mean two spread out into a ‘forest’ of lines because each line is things: (a) there is little hydrogen left in intergalactic space α redshifted by a different amount in proportion to the and by the time the Ly forest is observed most of the absorbing cloud’s distance from us. The stronger ones matter has already condensed into galaxies. Or (b) most among the absorption systems do show further spectral of the hydrogen is not in neutral form, where it can produce α signatures in addition to their Lyα line: higher-order Ly absorption, but is fully ionized. Lyman series lines begin to be detectable for absorption QSO surveys (see QUASISTELLAR OBJECTS: SURVEYS) later systems with Lyα close to saturation. Clouds with H I showed that the second possibility is more important: column densities larger than N ∼ 1017 cm−2 start showing the combined ionizing radiation output from all known a discontinuity due to continuous absorption at a rest QSOs at high redshift amounts to a UV radiation field frame wavelength 912 Å, beyond the limit of the Lyman probably strong enough to keep most of the baryonic series. These ‘Lyman limit systems’ occupy a column matter in the universe highly ionized—i.e., if the baryons density regime where a gas cloud starts shielding itself are predominantly in the form of a more or less against ionizing radiation from the outside. Clouds with homogeneously distributed gas. even higher column densities (N>1019 cm−2) exhibit This conclusion is based on the assumption that the damping wings caused by the internal finite lifetime the gas is in approximate photoionization equilibrium of the Lyα transition. The gas in these ‘damped Lyα’ with the cosmic UV background field, i.e., the rate of recombinations of with protons to form neutral systems (see LYMAN ALPHA ABSORPTION: THE DAMPED SYSTEMS)is hydrogen balances the rate of ionizations from the ground almost completely self-shielded and mostly neutral. Most . − state of H I, absorption systems with N>1014 5 cm 2 also show metal n n α(T ) = n . (1) absorption lines (triply ionized carbon and silicon, and e p HI some other common elements and ionization stages). For Here is the rate of photoionizations per neutral hydrogen that reason the higher column density systems are usually , caused by the hydrogen-ionizing portion of the UV n n n referred to as ‘metal’ or ‘heavy element’ systems. Here background field. The quantities e, HI and p refer to the we are concerned only with the low column density gas, number densities of electrons, neutral hydrogen , i.e., those absorption systems where the Lyα line is not and protons (= ionized hydrogen), respectively. saturated, which we will refer to as the Lyα forest proper. Later observations of QSOs with higher spectral There definitely is an overlap between metal systems and resolution (< a few hundred km s−1) showed that

Copyright © Nature Publishing Group 2001 Brunel Road, Houndmills, Basingstoke, Hampshire, RG21 6XS, UK Registered No. 785998 and Institute of Physics Publishing 2001 Dirac House, Temple Back, Bristol, BS1 6BE, UK 1 Lyman Alpha Forest E NCYCLOPEDIA OF A STRONOMY AND A STROPHYSICS

Figure 1. High resolution spectrum of the Lyα forest part of a redshift 3.63 QSO, taken with the HIRES spectrograph on the Keck 10 m telescope in Hawaii. The plot shows the flux of the QSO in arbitrary units versus the observed wavelength in units of Angstroms. The noise level can be judged from the longer wavelength wing of the broad, intrinsic Lyα emission line of the QSO (near 5650 Å). All of the ragged features are high redshift absorption lines. Most of the lines between the QSO’s Lyα and Lyβ emission lines (the humps at 5650 and 4750 Å) are due to Lyα absorption by intervening gas. The actual rate of incidence of absorbers decreases towards shorter (lower redshifts). Nevertheless, the line density increases to the blue because higher order absorption lines from the Lyman series appear and overlap randomly with Lyα lines of systems at lower redshift.

QSO’s flux I ∝ e−τ depends exponentially on the Lyα optical depth τ, which itself depends almost quadratically on the gas density (or the density ne; for a highly ionized gas at constant temperature). Thus small density fluctuations produce enhanced fluctuations in the optical depth. The Lyα forest absorption is observed in velocity space, and a convergent velocity field (e.g., a collapsing gas cloud) could also produce absorption ‘lines’. Caustics in velocity space may form if several gas volume elements are moving at the same velocity relative to the observer. Indeed, if Lyα clouds are produced by gravitational collapse, both overdense regions and infall should contribute to an absorption line. Spectroscopy of the Lyα line is an incredibly sensitive method to detect baryonic matter at any redshift. The photoionization cross section of neutral hydrogen is so large that an extremely tenuous gas at or below the Figure 2. Detailed section of the previous spectrum. The image mean density of the universe can be detected easily shows a number of absorption lines all corresponding to the in absorption. The method of choice for studying the α neutral hydrogen (H I) Ly 1215.67 Å transition. Lines close to Lyα forest is optical high resolution spectroscopy, with a saturation (= zero flux in the line center) have neutral hydrogen . − λ/λ > α column density typically around N ∼ 1014 2 cm 2, spectral resolution 30 000 sufficient to resolve Ly corresponding to a gas density enhanced by roughly an order of lines thermally broadened by the photoionization heating magnitude with respect to the mean density of the universe. The from the UV background. With 8 m class telescopes, a mean redshift of the stretch shown is z = 3.248. The spectral α − spectrum of a QSO suitable for further analysis of the Ly region extends over 3480 km s 1. For a flat = 1 universe this −1 forest absorption can be obtained within a few hours of corresponds to a spatial extent of approximately 9.6 h Mpc observing time. along the line-of-sight. Basic observational properties of Lyα absorbers what appeared as a smooth absorption ‘trough’ to Early models of Lyα absorption systems envisaged the earlier observers is in fact a ragged ‘forest’ of hundreds absorption lines as arising from discrete ‘clouds’ of gas of individual absorption lines. In other words, the in intergalactic space. The clouds would be analoguous to distribution of neutral hydrogen in the universe is galaxies, but the gas densities appeared too tenuous and inhomogeneous on scales down to the width of a typical too highly ionized to give rise to star formation. Moreover, Lyα line (see below). The degree of clumpiness appears the average clustering of the lines in velocity space was magnified by the absorption pattern in the QSO spectra, too weak for the clouds to be related directly to galaxies. because the residual (= unabsorbed) portion of the Thus the gas giving rise to Lyα absorption systems came

Copyright © Nature Publishing Group 2001 Brunel Road, Houndmills, Basingstoke, Hampshire, RG21 6XS, UK Registered No. 785998 and Institute of Physics Publishing 2001 Dirac House, Temple Back, Bristol, BS1 6BE, UK 2 Lyman Alpha Forest E NCYCLOPEDIA OF A STRONOMY AND A STROPHYSICS to be referred to as ‘Lyα clouds’ or ‘intergalactic clouds’, astrophysics (gravity, galactic outflows, nearby sources of to distinguish it from gas associated with galaxies. It it ionizing radiation) to affect the clustering of the lines in is worth keeping in mind that even the term ‘cloud’ does velocity space. already imply a prejudice as to the spatial distribution of (a) The evolution of the number of absorption lines N with the gas. time or redshift z (and with a column density above a Before the advent of cosmogonies dominated by certain threshold) traditionally has been approximated by DARK MATTER (which could have supplied the gravitational a power-law in (1+z), attraction necessary to hold together the gas), it was dN γ thought more likely that the intergalactic clouds would ∝ (1+z) . (3) be confined by the pressure of an even hotter and more dz tenuous intercloud medium. In this model, a cloud For a universe where Lyα clouds are non-evolving, 0 < would be homogeneous, static or slowly expanding and, γ<1/2, depending on the cosmological model. With at least for some of its lifetime, in thermal photoionization the observational baseline now extending from the local equilibrium with the UV background. universe out to redshift five it has become clear that γ is Such an object has a simple observational signature in not constant with time. In the local universe (z ∼ 0) up to absorption. With only the intrinsic atomic line width of the redshift ∼ 1.5, γ is consistent with the above no-evolution Lyα transition and purely thermal motions contributing values. However, beyond redshift z ∼ 1.5 there is a sharp to the absorption line broadening, an individual line increase in the number of absorption lines, with γ rising profile consists of a convolution of a Lorentzian resonance up to a value γ ∼ 5 by redshift z ∼ 4. curve and a Maxwell–Boltzmann velocity distribution. There is a general trend for the number of absorption The result is a so-called Voigt profile. Each such profile lines to increase with redshift, but in any individual QSO is fully specified by only 3 parameters: its position spectrum there is a relative lack of absorption systems close in velocity space (redshift), its column density (i.e. the to the redshift of the QSO. This so-called ‘proximity effect’ number of neutral hydrogen atoms per unit area), and its has been ascribed to the ionizing radiation produced by the line width (traditionally expressed in terms of the Doppler √ QSO itself, which reduces the neutral hydrogen fraction for parameter b = 2σ, where σ is the standard deviation α gas close to the QSO. The effect can be used to measure the of the Maxwellian velocity distribution). The Ly forest strength of the ionizing background radiation: consider a as a whole can then be characterized in terms of the point at a distance from the QSO, where the number of Lyα distribution functions of these three quantities. systems has declined by half as compared to the average line density in the Lyα forest. At this point the intensity of Column densities ionizing radiation from the QSO must equal the intensity Let us define dN /dN as the number of absorption lines, of the general UV background. Knowing the luminosity of dN , per unit redshift with an H I column density between the QSO we can compute the flux at that point, and thus, N and N +dN. This function tells us how likely it is for the intensity of the ionizing background. This method, our line-of-sight to a background QSO to intersect a cloud N though fraught with many uncertainties, has yielded the with a given H I column density . Observationally, for first measurement of the ionizing UV background. H I column densities N spanning the amazing range from 1012 to 1022 cm−2 the distribution was found to be well- (b) The small-scale distribution of absorption lines in parametrized by a single power law in column density, velocity space along the line-of-sight, often referred to as ‘clustering’, has yielded only limited information α dN − . about the nature of the Ly forest systems. If the ∝ N 1 5. (2) dN clustering is measured with the same methods used for galaxy surveys, namely by applying the two-point- To give an idea of the normalization, there are hundreds correlation function to discrete absorption lines, there is of unsaturated lines (N<1014.5 cm−2) per unit redshift a significant signal only for the highest column density but, typically, less than one damped Lyα system (N> (N>1015 cm−2) systems. These stronger absorption lines 1021 cm−2). For comparison, a QSO spectrum at redshift are invariably accompanied by absorption from heavier 3 covers z ∼ 0.6oftheLyα forest (between the QSO elements (carbon, silicon, oxygen, iron in various stages redshift and the redshift corresponding to the onset of the of ionization), so the high column density clouds must Lyβ series). in some way be more closely associated with galaxies as the production sites of the metals. The observed stronger The distribution of absorption systems in velocity space clustering of this gas may arise when gas clouds move in We may consider two limiting cases. (a) The distribution a galactic gravitational potential well, or they may just of absorbers over very large (Gigaparsec) scales; this is reflect stellar ejecta in the of those is equivalent to a distribution in time or redshift, as the galaxies. However, most of the Lyα forest consists of universe changes considerably during the time it takes a weaker lines, which do not cluster along the line-of-sight light ray to traverse a Gigaparsec. (b) The distribution appreciably. Gas densities inferred for the weaker lines on small (Megaparsec) scales; here we expect local are likely to be much lower than for galaxies. If gravity is

Copyright © Nature Publishing Group 2001 Brunel Road, Houndmills, Basingstoke, Hampshire, RG21 6XS, UK Registered No. 785998 and Institute of Physics Publishing 2001 Dirac House, Temple Back, Bristol, BS1 6BE, UK 3 Lyman Alpha Forest E NCYCLOPEDIA OF A STRONOMY AND A STROPHYSICS the main structure forming agent it is not very surprising Observations of close lines-of-sight to groups of QSOs that lower density regions are clustered more weakly. or to multiple images of gravitationally lensed QSOs can Searches for voids in the Lyα forest similar to those restore to some degree the missing second dimension seen in the spatial distribution of galaxies have been and provide a measurement of the cloud sizes, at least equally unsuccessful, implying again that most Lyα as projected in the plane of the sky. Imagine that we absorption systems are more homogeneously distributed had two lines of sight with a known separation. For in space than galaxies. Lyα ‘clouds’ with a given size, some of the clouds will intersect both lines-of-sight, others only one of them. In a Absorption line widths statistical sense, the typical size of Lyα forest absorbers can The width of an absorption line is a measure of the be determined from the numbers of ‘hits’ and ‘misses’, i.e., total velocity distribution in the gas. Both microscopic from knowing how often an absorption system seen in one (thermal motion) and macroscopic processes (turbulence, of the lines-of-sight also appears in the second one. Such bulk flows, broadening of an extended object by the estimates have yielded astonishingly large transverse sizes Hubble expansion) may contribute to the line profile. on the order of 0.1–1 Mpc proper separation. These sizes Without an a priori theory about line formation the total enable us to compute rough estimates of the ionization width b of the line profile can only be used as an upper fraction and density of the gas. limit to the width due to any individual process. As an Assuming that the neutral fraction x of hydrogen is example, the Doppler parameter of a thermally broadened determined by photoionization equilibrium (cf equation absorption line with an additional Gaussian contribution (1)), we have describing turbulence in the gas is given exactly by     T −0.35 I −0.5  x = . × −6 3 9 10 − 2kT 3 × 104 10 21 b = + b2 (4)     m turb N 0.5 D −0.5 × , (7) 1014 100 kpc where T is the temperature of the gas, k is Boltzmann’s m b constant, is the mass of the hydrogen atom, and turb is where T is the gas temperature in K, I is the intensity of the turbulent contribution to the Doppler parameter. the ionizing UV background in units of ergs Hz−1 sr−1 s−1 ∼ At redshift 3 observations show the Doppler cm−2, N the H I column density in cm−2, and D the b parameters to be distributed according to a Gaussian thickness of the cloud (or the path length of our line of b −1 b with a mean between 25 and 30 km s , and a cutoff cut sight through it). −1 below about 19 km s , i.e., Likewise, the total number density of the gas (the   number of protons per cm3) is given by N (b − b)2 d ∝ b ≥ b exp for cut (5)     db 2σ 2 T 0.35 I 0.5 b n = x−1ND−1 = . × −5 8 3 10 − = 0 for b

Size, density, and ionization state of Lyα clouds The intergalactic medium as a cosmic fluid Absorption line studies suffer from the limitation that in So far we have used only simple astrophysical arguments general only one-dimensional information along the line- to infer the basic physical properties of the Lyα absorbers, of-sight is available. From the absorption line itself we without explicit reference to a model of structure measure the temperature and column density of the gas, formation. The decomposition into Voigt profiles had and from observations of the QSO luminosity function originally been justified by the assumption that Lyα we can compute an estimate of the ionizing background absorption systems are discrete, pressure-confined gas radiation, but to obtain the density and the ionized fraction clouds in hydrostatic, thermal and ionization equilibrium, of the gas an additional independent estimate of the size floating in intergalactic space. This picture provided of the absorbing gas cloud is needed. an analytically tractable model which made a host of

Copyright © Nature Publishing Group 2001 Brunel Road, Houndmills, Basingstoke, Hampshire, RG21 6XS, UK Registered No. 785998 and Institute of Physics Publishing 2001 Dirac House, Temple Back, Bristol, BS1 6BE, UK 4 Lyman Alpha Forest E NCYCLOPEDIA OF A STRONOMY AND A STROPHYSICS observationally testable predictions, several of which unfortunately contradicted the observational evidence. For example, the large range in column densities observed, and the evolution of the number of absorption systems with redshift are difficult to reproduce in that model. Moreover, observations of the spectral shape of the cosmic microwave background radiation with the COBE satellite have provided arguments against the existence of a hot intercloud medium necessary to provide the confinement pressure for the Lyα clouds. The pressure-confinement model finally fell from grace because of its lack of a convincing theory of the formation of the cloudlets, and the rise (in popularity) of structure formation cosmogonies dominated by the presence of large amounts of weakly interacting dark matter. The gravitational collapse of dark matter would have trapped large amounts of baryonic gas as well. The most popular, cold dark matter (CDM) structure formation scenario predicts a large abundance of collapsed CDM halos with individual masses too small to form stars and turn into galaxies. Warm photoionized intergalactic gas sinks into these ‘mini’-halos or accretes onto dark matter Figure 3. filaments and sheets. The thermal gas pressure prevents Cosmological hydro-simulation: projection of the baryonic density distribution in a simulated box at z = 3 (same the gas from further collapse; relatively stable gaseous spatial extent as for figure 2) showing the ‘cosmic web’ configurations are formed. These structures are visible filamentary structure of the Lyα forest (image courtesy Michael only in absorption since there are no stars which could Norman, University of Illinois). produce any light. A breakthrough in cosmology occured in the early 1990s when hydrodynamic cosmological simulations settles in a dense, central cooling layer, presumably to became available. As far as the comparison with data form stars in some of the denser regions. At the lowest is concerned the hydrodynamics simulations represent column densities gas remains unshocked and just bounces an advance over pure dark matter (so-called N-body) back because of the hydrostatic pressure. The gas is simulations, as they attempt to directly predict observable partly confined by dark matter gravity and partly by ram- astrophysical quantities. The hydro-codes include, in pressure. Higher column density clouds arise in more a simplified way, gas dynamics, elementary radiative filamentary structures, with column density contours of − processes, heating and cooling, and some schematic stellar N ∼ 1014 cm 2 extending continuously and at relatively feedback, all of which are essential to an understanding of constant thickness (∼ 40–100 kpc proper) over Mpc the formation of structure for ordinary (baryonic) matter distances. With increasing column density the absorber like galaxies, stars and intergalactic gas clouds. The geometry becomes rounder; column density contours at − ultimate hope of the hydro-simulations is to reproduce N  1016cm 2 invariably are spherical. Such absorbers realistic galaxies, but at the current level of detail and more closely correspond to the aforementioned minihalos; spatial resolution possible the intergalactic medium with there the enclosed gas column is high enough to make the its simple physics is perhaps the most promising target for absorption system appear as a Lyman limit or damped Lyα quantitative modelling. system. Looking at the higher column density, optically thick gas on scales of several Mpc one gets a somewhat The cosmic web different impression of chains of mini- or larger halos, If the underlying cosmological picture (a universe lining up like pearls on a string, quite similar to the dominated by cold dark matter) is correct, then the hydro- structure seen in N-body simulations of the dark matter simulations are telling us that the spatial distribution and distribution. To produce as much absorption as observed, physical state of the Lyα forest gas is more complex than a large fraction of all baryons (80–90%) is required to reside previously thought. The gas is arranged in filaments in the low column density Lyα forest, mostly in the column and sheets, in what has been called the ‘cosmic web’, density range 1014

Copyright © Nature Publishing Group 2001 Brunel Road, Houndmills, Basingstoke, Hampshire, RG21 6XS, UK Registered No. 785998 and Institute of Physics Publishing 2001 Dirac House, Temple Back, Bristol, BS1 6BE, UK 5 Lyman Alpha Forest E NCYCLOPEDIA OF A STRONOMY AND A STROPHYSICS with redshift is largely a consequence of the expansion of This equation relates the optical depth to the mean the gas with the Hubble flow. The statistical distribution baryonic density (in gas) in units of the critical density, −τ of the flux level I = e in the Lyα forest is extremely well b, the Hubble constant at redshift z, H(z), the average matched by these models, as can be seen from a compar- gas temperature T , the proper baryon density ρ, the ison of the observed distribution with one from a simu- photoionization rate , and the gradient of the local  v / r lated CDM universe (figure 4). The spatial correlations peculiar velocity d pec d along the line-of-sight. A as evident from the flux power-spectrum are equally well further convolution with a Voigt profile is necessary to reproduced. The most complex piece of observational in- include the proper thermal velocity broadening. The formation, the distribution of Doppler parameters is quali- exponent α (α = 2 for an isothermal gas) takes account tatively understood, but as it depends on a variety of initial of the fact that in denser regions of the universe the conditions (epoch of , baryonic density of the gas is typically warmer because it is more effectively universe, cosmological model) and possible stellar energy heated by photoionization, but α also depends on the feedback, several of these parameters need to be tuned reionization history of the gas and the amount of adiabatic carefully to get acceptable quantitative agreement. expansion/compression. Values of α ≈ 1.6–1.8 are given Perhaps most importantly, the successes of the hydro- in the literature. α models show that the majority of the Ly absorbers are To turn this relation into a complete description of consistent with being low density condensations formed the observed Lyα forest, cosmology has to predict the by gravitational collapse of the intergalactic medium. In cosmic density and velocity fields, the fraction of the α this picture the gas seen as Ly forest is the original closure density in the form of gas, the equation of state reservoir of matter from which galaxies are condensing. of the intergalactic medium, and the ionizing radiation Conversely, the success in reproducing the observa- field. Measuring cosmological parameters then boils α tions of the Ly forest provides perhaps the best obser- down to creating Lyα forest spectra according to a given vational evidence for hierarchical structure formation we cosmological prescription, and iterating with varying have to date. input parameters until good agreement between observed and predicted properties is obtained. At the time of writing Cosmology with the Lyα forest Most of the intergalactic gas observed at high redshift these techniques are just beginning to be explored, but it has experienced only mild gravitational collapse. A has become obvious that there are at least three main areas typical region of space has undergone little processing to benefit from such studies. other than heating by photoionization and compression, From equation (9) the optical depth is directly ( h2 )2/ in competition with adiabatic cooling by the Hubble proportional to the ratio b 50 . In other words, a expansion. It turns out that the weaker, unsaturated higher density produces stronger absorption, but a higher absorption lines (N<1014 cm−2) are still on or near photoionization rate reduces the neutral fraction of the gas the linear regime of gravitational collapse. Even for and decreases the absorption. With a suitable hypothesis lines close to saturation the overdensities with respect for or an independent measurement of the photoionization to the mean density of the universe are less than about rate (e.g. estimating the density of ionizing photons from a factor δ ∼ 10–15. On spatial scales on the order of QSO surveys) the baryon density b of the universe may one Mpc and larger where the thermal pressure of the be constrained. gas is not important the intergalactic medium traces the The ‘equation of state’ of the intergalactic medium, underlying mass distribution much more closely than the as the statistical distribution of the gas volume elements T − ρ stellar light of galaxies observed in emission. This fact in temperature-density ( ) space has become known, and the simple astrophysics involved should make the has a more complex and subtle influence on the optical intergalactic medium an ideal cosmological laboratory. depth distribution: a change in the slope dT/dρ changes Recently there have been various attempts at tapping the value of the temperature of a given volume element the cosmological potential of the Lyα forest. The link and thus the recombination coefficient α(T ), the neutral between the observable appearance of the Lyα forest fraction and column density, and the absorption line and the various cosmological input parameters can be width. It appears that the distribution of the Lyα forest described approximately by the Gunn–Peterson relation lines in Doppler parameter–column density (b − N) space for the H I optical depth τ, generalized to include an is a distorted map of the density–temperature relation, and inhomogenous density and velocity field. As long as the can be used to constrain the latter. T(ρ), in turn, contains gas is highly ionized and in photoionization equilibrium information on the epoch of reionization and reheating (not necessarily thermal equilibrium), and the gas is and on the sources of ionizing radiation. unshocked, the optical depth for Lyα absorption at redshift Finally, it is clear that the fluctuation properties z is proportional to (amplitude, spatial correlation) of the density ρ and   peculiar velocity v fields propagate through to the ( h2 )2 H(0) pec τ(z) ∝ b 50 T −0.7 optical depth (cf equation (9)). At least in the regime h H(z) α    50  where the Ly lines are strong enough to be detected but α v −1 ρ d pec not yet too strong to be saturated, the spectra provide a × (1+z)6 1+ (9) ρ H(z)dr record of the initial conditions of gravitational structure

Copyright © Nature Publishing Group 2001 Brunel Road, Houndmills, Basingstoke, Hampshire, RG21 6XS, UK Registered No. 785998 and Institute of Physics Publishing 2001 Dirac House, Temple Back, Bristol, BS1 6BE, UK 6 Lyman Alpha Forest E NCYCLOPEDIA OF A STRONOMY AND A STROPHYSICS

Figure 4. Observed and simulated statistical distributions of the absorbed flux level in the Lyα forest, for three different mean redshifts. The solid lines show the distribution for a sample of QSOs observed with the Keck telescope, the dotted lines are from a simulation of a CDM universe (performed by Renyue Cen et al at Princeton University). The good agreement in the shape of the distributions lends observational support to the Lyα forest being a by-product of hierarchical structure formation.

formation. In principle the Lyα forest spectra can be AGN ABSORPTION LINES). These include metal absorption inverted to provide a measurement of the power spectrum systems in general, and the correspondence between shape and amplitude of the initial fluctuations. absorption systems and various galactic or interstellar environments, as seen in local galaxies. Damped Lyα The helium Lyα forest systems, the absorbers most relevant to high redshift We conclude with a glance at one of the most intriguing galaxy formation, are treated in the separate article LYMAN α new topics of research, the Ly forest absorption due ALPHA ABSORPTION: THE DAMPED SYSTEMS. to singly ionized helium (He II). The He II forest has only recently become acessible as the transition is in Bibliography the far ultraviolet (304 Å) and even if redshifted out to Reviews of the subject at a level accessible to students, and z ∼ 3 still needs to be observed with a UV spectrograph detailed references to the literature are given in from space. In addition, it is very difficult to find a QSO with its UV continuum intact, as the Lyman limit Weymann R J, Carswell R F and Smith M G 1981 continuum absorption from higher redshift systems tends Absorption lines in the spectra of quasistellar objects to obliterate the far UV in most QSOs. Nevertheless, there Ann. Rev. Astron. Astrophys. 19 41 are potentially big scientific rewards to be gained: (1) since Rauch M 1998 The Lyman Alpha Forest in the spectra of He II and H I have different ionization potentials, looking QSOs Ann. Rev. Astron. Astrophys. 36 267 at the ratio of the column densities from their respective Lyα forests one can constrain the shape of the UV spectrum The scientific literature on the subject begins with the first ionizing the two species. (2) Helium and hydrogen atomic measurement of the Gunn–Peterson effect, published as masses differ sufficiently that a comparison of the Doppler Gunn J E and Peterson B A 1965 On the density of neutral parameters of both can measure the amount of non- hydrogen in intergalactic space Astrophys. J. 142 1633 thermal line broadening. (3) The first observations of He II forests are showing that He II ionization appears A number of conference proceedings give very useful to be patchy and not fully developed by redshift 3–4. This overviews at a somewhat more technical level than the observation could help us to understand how and when reviews; they are indispensible as guides to the many the universe was first ionized. topics related to QSO absorption lines which could not There are many other topics related to the rapidly be treated here: growing field of QSO absorption lines for which we refer the interested reader to the literature given below (see also Blades J C, Turnshek D A and Norman C (eds) 1988 QUASISTELLAR OBJECTS: INTERVENING ABSORPTION LINES; INTRINSIC QSO Absorption Lines: Probing the Universe (Proc. QSO

Copyright © Nature Publishing Group 2001 Brunel Road, Houndmills, Basingstoke, Hampshire, RG21 6XS, UK Registered No. 785998 and Institute of Physics Publishing 2001 Dirac House, Temple Back, Bristol, BS1 6BE, UK 7 Lyman Alpha Forest E NCYCLOPEDIA OF A STRONOMY AND A STROPHYSICS

Absorption Line Meeting, Baltimore 1987) (Cambridge: Cambridge University Press) Meylan G (ed) 1995 QSO Absorption Lines (Proc. ESO Workshop, November 21–24, 1994) (Berlin: Springer) Petitjean P and Charlot S 1997 Structure and Evolution of the Intergalactic Medium from QSO Absorption Lines (Proc. 13th IAP Astrophysics Colloquium) (Paris: Editions Fronti`eres)

Recent work on cosmological applications of the Lyα forest is discussed in articles by Weinberg et al, Hui, Haehnelt and Nusser in

Banday A J, Sheth R K and Da Costa L N (eds) 1999 Evolution of Large Scale Structure from Recombination to Garching (Proc. MPA-Garching Cosmology Conference) (Enschede: PrintPartners Ipskamp)

Michael Rauch

Copyright © Nature Publishing Group 2001 Brunel Road, Houndmills, Basingstoke, Hampshire, RG21 6XS, UK Registered No. 785998 and Institute of Physics Publishing 2001 Dirac House, Temple Back, Bristol, BS1 6BE, UK 8