Cosmological reionization around the first stars: Monte Carlo radiative transfer

We study the evolution of ionization fronts around the first protogalaxies by using high-resolution numerical cosmological (Λ+ cold dark matter, CDM, model) simulations and Monte Carlo radiative transfer methods. We present the numerical scheme in detail and show the results of test runs from which we conclude that the scheme is both fast and accurate. As an example of interesting cosmological application, we study the reionization produced by a stellar source of total mass M =2×10⁸ M_⊙ turning on at z ≈12, located at a node of the cosmic web. The study includes a spectral energy distribution of a zero-metallicity stellar population, and two initial mass functions (IMFs; Salpeter/Larson). The expansion of the ionization front (I-front) is followed as it breaks out from the galaxy and is channelled by the filaments into the voids, assuming (in a 2D representation) a characteristic butterfly shape. The ionization evolution is very well tracked by our scheme, as realized by the correct treatment of the channelling and shadowing effects resulting from overdensities. We confirm previous claims that both the shape of the IMF and the ionizing power metallicity dependence are important to correctly determine the reionization of the Universe.     

Zero-metallicity stellar sources and the reionization epoch

We reconsider the problem of the cosmological reionization owing to stellar sources. Using a method similar to that developed by Haiman & Loeb, we investigate the effect of changing the stellar models and the stellar spectra adopted for deriving the ionizing photon production rate. In particular, we study the consequences of adopting zero-metallicity stars, which is the natural choice for the first stellar populations. We construct young isochrones representative of Population III stars from existing sets of evolutionary models (by Forieri and Cassisi & Castellani) and calculate a suitable library of zero-metallicity model atmospheres. The number of ionizing photons emitted by such a zero-metal population is about 40 per cent higher than that produced by standard metal-poor isochrones. We find that adopting suitable zero-metallicity models modifies the reionization epoch. However the latter is still largely affected by current uncertainties in other important physical processes such as the efficiency of the star formation and the fraction of escaping UV photons.     

The reionization of He II and the temperature evolution of the intergalactic medium

A number of observations suggest that He II in the intergalactic medium(IGM) was fully ionized at z ～ 3, probably by quasi-stellar objects(QSOs). Here we construct a simple model of a QSO to study the reionization of He II and the corresponding thermal evolution of the IGM. We assume that QSOs are triggered by major mergers of dark matter halos, and the luminosity evolution of individual QSOs is described by an initial accretion stage with a constant Eddington ratio and then a powerlaw decay driven by long term disk evolution or fueling. Once a QSO is triggered, it immediately ionizes its surrounding area as an ionized bubble. The resulting changes in size and volume of the bubble are determined by the luminosity evolution of the central QSO. With the emergence of more and more bubbles, they eventually overlap each other and finally permeate the whole universe. During the He II reionization,the IGM temperature increases due to the photoheating by the ionization processes.Applying the bubble model and considering various heating and cooling mechanisms,we trace the thermal evolution of the IGM and obtain the average IGM temperature as a function of redshift, which is very consistent with observations. The increase in IGM temperature due to the reionization of He II may be determined more accurately in the future, which may put robust constraints on the QSO model and the physics of He II reionization.     

Scintillation as a probe of the intergalactic medium

Most of the baryons in the low-redshift Universe reside in a warm/hot component which is difficult to detect with standard absorption/emission-line techniques. We propose to use quasar refractive scintillation as a useful, complementary probe for such ionized, intergalactic gas. In particular, an application to the case of the intracluster medium is presented. We show that clusters located at z ≈0.02 should produce a source rms intensity fluctuation at 50–100 GHz of several tens of per cent and on time-scales ranging from days to months, depending on the projected location of the source on the foreground cluster. However, in order to produce such a signal, the source needs to be very compact. This effect, if observed, can be used as an independent test of the baryonic mass fraction in clusters.     

Metal enrichment of the intergalactic medium

I demonstrate by means of high-resolution cosmological simulations, which include modelling of a two-phase interstellar medium, that the dominant mechanism for transporting heavy elements from protogalaxies into the intergalactic medium (IGM) is the merger mechanism as discovered by Gnedin & Ostriker. Direct ejection of the interstellar gas by supernovae plays only a minor role in transporting metals into the IGM: for a realistic cosmological scenario only a small fraction of all metals in the IGM is delivered by the supernova-driven winds, while most of the metals in the IGM are transported by the merger mechanism. As a result, the metallicity distribution in the IGM is highly inhomogeneous, in agreement with studies of the QSO metal absorption systems, and the predicted metallicity distribution of Lyman alpha absorbers as a function of their column density is in excellent agreement with the observational data.     

Source mergers and bubble growth during reionization

The recently introduced models of reionization bubbles based on extended Press–Schechter theory (by Furlanetto, Hernquist & Zaldarriaga) are generalized to include mergers of ionization sources. Sources with a recent major merger are taken to have enhanced photon production due to star formation, and accretion on to a central black hole if a black hole is present. This produces a scatter in the number of ionized photons corresponding to a halo of a given mass and a change in photon production over time for any given halo mass. By extending previous methods, photon production histories, bubble distributions and ionization histories are computed for several different parameter and recombination assumptions. The resulting distributions interpolate between previously calculated limiting cases.     

Radiation from early black holes – I. Effects on the neutral intergalactic medium

In the pre-reionization Universe, the regions of the intergalactic medium (IGM) which are far from luminous sources are the last to undergo reionization. Until then, they should be scarcely affected by stellar radiation; instead, the X-ray emission from an early black hole (BH) population can have much larger influence. We investigate the effects of such emission, looking at a number of BH model populations (differing for the cosmological density evolution of BHs, the BH properties, and the spectral energy distribution of the BH emission). We find that BH radiation can easily heat the IGM to  10³–10⁴ K  , while achieving partial ionization. The most interesting consequence of this heating is that BHs are expected to induce a 21-cm signal (  δ T _b∼ 20–30 mK  at   z ≲ 12  ) which should be observable with forthcoming experiments (e.g. LOFAR). We also find that at   z ≲ 10  BH emission strongly increases the critical mass separating star-forming and non-star-forming haloes.     

Probing the neutral fraction of the IGM with GRBs during the epoch of reionization

We show that near-infrared observations of the red side of the Lyα line from a single gamma-ray burst (GRB) afterglow cannot be used to constrain the global neutral fraction of the intergalactic medium (IGM),     , at the GRB's redshift to better than     . Some GRB sightlines will encounter more neutral hydrogen than others at fixed     owing to the patchiness of reionization. GRBs during the epoch of reionization will often bear no discernible signature of a neutral IGM in their afterglow spectra. We discuss the constraints on     from the   z = 6.3  burst, GRB050904, and quantify the probability of detecting a neutral IGM using future spectroscopic observations of high-redshift, near-infrared GRB afterglows. Assuming an observation with signal-to-noise ratio similar to the Subaru FOCAS spectrum of GRB050904 and that the column density distribution of damped Lyα absorbers is the same as measured at lower redshifts, a GRB from an epoch when     can be used to detect a partly neutral IGM at 97 per cent confidence level ≈10 per cent of the time (and, for an observation with three times the sensitivity, ≈30 per cent of the time).