We investigate the feasibility of detecting 21cm absorption features in the afterglow spectra of high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic reionization, together with the instrumental characteristics of
the LOw Frequency ARray (LOFAR). We find that absorption features could be marginally (with a S/N larger than a few) detected by LOFAR at z>7 if the GRB originated from PopIII stars, while the detection would be easier if the noise were reduced by one order of magnitude, i.e. similar to what is expected for the first phase of the Square Kilometer Array (SKA1-low). On the other hand, more standard GRBs are too dim to be detected even with ten times the sensitivity of SKA1-low, and only in the most optimistic case can a S/N larger than a few be reached at z>9.
Using a combination of N-body simulations, semi-analytic models and radiative transfer calculations, we have estimated the theoretical cross power spectrum between galaxies and the 21cm emission from neutral hydrogen during the epoch of reionization.
In accordance with previous studies, we find that the 21cm emission is initially correlated with halos on large scales (> 30 Mpc), anti-correlated on intermediate (~ 5 Mpc), and uncorrelated on small (< 3 Mpc) scales. This picture quickly changes as reionization proceeds and the two fields become anti-correlated on large scales. The normalization of the cross power spectrum can be used to set constraints on the average neutral fraction in the intergalactic medium and its shape can be a tool to study the topology of reionization. When we apply a drop-out technique to select galaxies and add to the 21cm signal the noise expected from the LOFAR telescope, we find that while the normalization of the cross power spectrum remains a useful tool for probing reionization, its shape becomes too noisy to be informative. On the other hand, for a Lyalpha Emitter (LAE) survey both the normalization and the shape of the cross power spectrum are suitable probes of reionization. A closer look at a specific planned LAE observing program using Subaru Hyper-Suprime Cam reveals concerns about the strength of the 21cm signal at the planned redshifts. If the ionized fraction at z ~ 7 is lower that the one estimated here, then using the cross power spectrum may be a useful exercise given that at higher redshifts and neutral fractions it is able to distinguish between two toy models with different topologies.
We discuss the feasibility of the detection of the 21cm forest in the diffuse IGM with the radio telescope LOFAR. The optical depth to the 21cm line has been derived using simulations of reionization which include detailed radiative transfer of ioniz
ing photons. We find that the spectra from reionization models with similar total comoving hydrogen ionizing emissivity but different frequency distribution look remarkably similar. Thus, unless the reionization histories are very different from each other (e.g. a predominance of UV vs. x-ray heating) we do not expect to distinguish them by means of observations of the 21cm forest. Because the presence of a strong x-ray background would make the detection of 21cm line absorption impossible, the lack of absorption could be used as a probe of the presence/intensity of the x-ray background and the thermal history of the universe. Along a random line of sight LOFAR could detect a global suppression of the spectrum from z>12, when the IGM is still mostly neutral and cold, in contrast with the more well-defined, albeit broad, absorption features visible at lower redshift. Sharp, strong absorption features associated with rare, high density pockets of gas could be detected also at z~7 along preferential lines of sight.
We investigate the effect of helium on hydrogen reionisation using a hydrodynamical simulation combined with the cosmological radiative transfer code CRASH. The simulations are run in a 35.12/h comoving Mpc box using a variety of assumptions for the
amplitude and power-law extreme-UV (EUV) spectral index, alpha, of the ionising emissivity. We use an empirically motivated prescription for ionising sources which ensures all of the models are consistent with constraints on the Thomson scattering optical depth and the hydrogen photo-ionisation rate at z=6. The inclusion of helium slightly delays reionisation due to the small number of ionising photons which reionise neutral helium instead of hydrogen. However, helium has a significant impact on the thermal state of the IGM. Models with alpha=3 produce IGM temperatures at the mean density at z=6 which are about 20 % higher compared to models without helium photo-heating. Harder EUV indices produce even larger IGM temperature boosts. A comparison to recent observational estimates of the IGM temperature at z=5 - 6 suggests that hydrogen reionisation was primarily driven by pop-II stellar sources with a soft EUV index, alpha<3. We also find that faint, as yet undetected galaxies, characterised by a luminosity function with a steepening faint-end slope and an increasing Lyman continuum escape fraction (fesc=0.5), are required to reproduce the ionising emissivity used in our simulations at z>6. Finally, we note there is some tension between recent observational constraints which indicate the IGM is > 10% neutral by volume z=7, and estimates of the ionising emissivity at z=6 which indicate only between 1 and 3 ionising photons are emitted per hydrogen atom over a Hubble time. This tension may be alleviated by either a lower neutral fraction at z=7 or an IGM which still remains a few % neutral by volume at z=6.
In this paper we report on the improvements implemented in the cosmological radiative transfer code CRASH. In particular we present a new multi-frequency algorithm for spectra sampling which makes use of colored photon packets: we discuss the need fo
r the multi-frequency approach, describe its implementation and present the improved CRASH performance in reproducing the effects of ionizing radiation with an arbitrary spectrum. We further discuss minor changes in the code implementation which allow for more efficient performance and an increased precision.
In this paper we present CRASH_alpha, the first radiative transfer code for cosmological application that follows the parallel propagation of Ly_alpha and ionizing photons. CRASH_alpha is a version of the continuum radiative transfer code CRASH with
a new algorithm to follow the propagation of Ly_alpha photons through a gas configuration whose ionization structure is evolving. The implementation introduces the time evolution for Ly_alpha photons (a feature commonly neglected in line radiative transfer codes) and, to reduce the computational time needed to follow each scattering, adopts a statistical approach to the Ly_alpha treatment by making extensive use of pre-compiled tables. With this statistical approach we experience a drastic increase of the computational speed and, at the same time, an excellent agreement with the full Ly_alpha radiative transfer computations of the code MCLy_alpha. We find that the emerging spectra keep memory of the ionization history which generates a given ionization configuration of the gas and, to properly account for this effect, a self-consistent joint evolution of line and ionizing continuum radiation as implemented in CRASH_alpha is necessary. A comparison between the results from our code and from Ly_alpha scattering alone on a fixed HI density field shows that the extent of the difference between the emerging spectra depends on the particular configuration considered, but it can be substantial and can thus affect the physical interpretation of the problem at hand. These differences should furthermore be taken into account when computing the impact of the Ly_alpha radiation on e.g. the observability of the 21 cm line from neutral hydrogen at epochs preceeding complete reionization.
Cooling is the main process leading to the condensation of gas in the dark matter potential wells and consequently to star and structure formation. In a metal-free environment, the main available coolants are H, He, H$_2$ and HD; once the gas is enri
ched with metals, these also become important in defining the cooling properties of the gas. We discuss the implementation in Gadget-2 of molecular and metal cooling at temperatures lower that $rm10^4 K$, following the time dependent properties of the gas and pollution from stellar evolution. We have checked the validity of our scheme comparing the results of some test runs with previous calculations of cosmic abundance evolution and structure formation, finding excellent agreement. We have also investigated the relevance of molecule and metal cooling in some specific cases, finding that inclusion of HD cooling results in a higher clumping factor of the gas at high redshifts, while metal cooling at low temperatures can have a significant impact on the formation and evolution of cold objects.
This submission has been withdrawn by arXiv administrators because it is a duplicate of 0704.2182.
