Reionization and cosmic dawn astrophysics from the Square Kilometre Array:impact of observing strategies

Interferometry of the cosmic 21-cm signal is set to revolutionize our understanding of the epoch of reionization (EoR) and the cosmic dawn (CD). The culmination of ongoing efforts will be the upcoming Square Kilometre Array (SKA), which will provide tomography of the 21-cm signal from the first billion years of our Universe. Using a galaxy formation model informed by high-z luminosity functions, here we forecast the accuracy with which the first phase of SKA-low (SKA1-low) can constrain the properties of the unseen galaxies driving the astrophysics of the EoR and CD. We consider three observing strategies: (i) deep (1000 h on a single field); (ii) medium-deep (100 h on 10 independent fields); and (iii) shallow (10 h on 100 independent fields). Using the 21-cm power spectrum as a summary statistic, and conservatively only using the 21-cm signal above the foreground wedge, we predict that all three observing strategies should recover astrophysical parameters to a fractional precision of 3c0.1-10 per cent. The reionization history is recovered to an uncertainty of \u394z 7e 0.1 (1\u3c3 ) for the bulk of its duration. The medium-deep strategy, balancing thermal noise against cosmic variance, results in the tightest constraints, slightly outperforming the deep strategy. The shallow observational strategy performs the worst, with up to an 3c10-60 per cent increase in the recovered uncertainty. We note, however, that non-Gaussian summary statistics, tomography, as well as unbiased foreground removal would likely favour the deep strategy

The signature of the first stars in atomic hydrogen at redshift 20

Dark and baryonic matter moved at different velocities in the early Universe, which strongly suppressed star formation in some regions. This was estimated to imprint a large-scale fluctuation signal of about 2 mK in the 21-cm spectral line of atomic hydrogen associated with stars at a redshift of 20, although this estimate ignored the critical contribution of gas heating due to X-rays and major enhancements of the suppression. A large velocity difference reduces the abundance of halos and requires the first stars to form in halos of about a million solar masses, substantially greater than previously expected. Here we report a simulation of the distribution of the first stars at z=20 (cosmic age of ~180 Myr), incorporating all these ingredients within a 400 Mpc box. We find that the 21-cm signature of these stars is an enhanced (10 mK) fluctuation signal on the 100-Mpc scale, characterized by a flat power spectrum with prominent baryon acoustic oscillations. The required sensitivity to see this signal is achievable with an integration time of a thousand hours with an instrument like the Murchison Wide-field Array or the Low Frequency Array but designed to operate in the range of 50-100 MHz.Comment: 27 pages, 5 figures, close (but not exact) match to accepted version. Basic results unchanged from first submitted version, but justification strengthened, title and abstract modified, and substantial Supplementary Material added. Originally first submitted for publication on Oct. 12, 201

Observing the First Stars and Black Holes

The high sensitivity of JWST will open a new window on the end of the cosmological dark ages. Small stellar clusters, with a stellar mass of several 10^6 M_sun, and low-mass black holes (BHs), with a mass of several 10^5 M_sun should be directly detectable out to redshift z=10, and individual supernovae (SNe) and gamma ray burst (GRB) afterglows are bright enough to be visible beyond this redshift. Dense primordial gas, in the process of collapsing from large scales to form protogalaxies, may also be possible to image through diffuse recombination line emission, possibly even before stars or BHs are formed. In this article, I discuss the key physical processes that are expected to have determined the sizes of the first star-clusters and black holes, and the prospect of studying these objects by direct detections with JWST and with other instruments. The direct light emitted by the very first stellar clusters and intermediate-mass black holes at z>10 will likely fall below JWST's detection threshold. However, JWST could reveal a decline at the faint-end of the high-redshift luminosity function, and thereby shed light on radiative and other feedback effects that operate at these early epochs. JWST will also have the sensitivity to detect individual SNe from beyond z=10. In a dedicated survey lasting for several weeks, thousands of SNe could be detected at z>6, with a redshift distribution extending to the formation of the very first stars at z>15. Using these SNe as tracers may be the only method to map out the earliest stages of the cosmic star-formation history. Finally, we point out that studying the earliest objects at high redshift will also offer a new window on the primordial power spectrum, on 100 times smaller scales than probed by current large-scale structure data.Comment: Invited contribution to "Astrophysics in the Next Decade: JWST and Concurrent Facilities", Astrophysics & Space Science Library, Eds. H. Thronson, A. Tielens, M. Stiavelli, Springer: Dordrecht (2008

Measuring the density fields around bright quasars at z ∼6 with XQR-30 Spectra

Measuring the density of the intergalactic medium using quasar sight lines in the epoch of reionization is challenging due to the saturation of Lyα absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Lyα absorption is not saturated. In this study, we use 10 high-resolution (R â 3 10,000) z ∼6 quasar spectra from the extended XQR-30 sample to measure the density field in the quasar proximity zones. We find a variety of environments within 3 pMpc distance from the quasars. We compare the observed density cumulative distribution function (CDF) with models from the Cosmic Reionization on Computers simulation and find a good agreement between 1.5 and 3 pMpc from the quasar. This region is far away from the quasar hosts and hence approaching the mean density of the universe, which allows us to use the CDF to set constraints on the cosmological parameter σ 8 = 0.6 ± 0.3. The uncertainty is mainly due to the limited number of high-quality quasar sight lines currently available. Utilizing the more than 200 known quasars at z â 3 6, this method will allow us to tighten the constraint on σ 8 to the percent level in the future. In the region closer to the quasar within 1.5 pMpc, we find that the density is higher than predicted in the simulation by 1.23 ± 0.17, suggesting that the typical host dark matter halo mass of a bright quasar (M 1450 <-26.5) at z ∼6 is log10(Mh/M⊠)=12.5-0.7+0.4

Fast Large-Scale Reionization Simulations

We present an efficient method to generate large simulations of the Epoch of Reionization (EoR) without the need for a full 3-dimensional radiative transfer code. Large dark-matter-only simulations are post-processed to produce maps of the redshifted 21cm emission from neutral hydrogen. Dark matter haloes are embedded with sources of radiation whose properties are either based on semi-analytical prescriptions or derived from hydrodynamical simulations. These sources could either be stars or power-law sources with varying spectral indices. Assuming spherical symmetry, ionized bubbles are created around these sources, whose radial ionized fraction and temperature profiles are derived from a catalogue of 1-D radiative transfer experiments. In case of overlap of these spheres, photons are conserved by redistributing them around the connected ionized regions corresponding to the spheres. The efficiency with which these maps are created allows us to span the large parameter space typically encountered in reionization simulations. We compare our results with other, more accurate, 3-D radiative transfer simulations and find excellent agreement for the redshifts and the spatial scales of interest to upcoming 21cm experiments. We generate a contiguous observational cube spanning redshift 6 to 12 and use these simulations to study the differences in the reionization histories between stars and quasars. Finally, the signal is convolved with the LOFAR beam response and its effects are analyzed and quantified. Statistics performed on this mock data set shed light on possible observational strategies for LOFAR.Comment: 18 pages, 21 figures, submitted to MNRAS For high-resolution images follow "http://www.astro.rug.nl/~thomas/eormap.pdf

GRB 050904 at redshift 6.3: observations of the oldest cosmic explosion after the Big Bang

We present optical and near-infrared observations of the afterglow of the gamma-ray burst GRB 050904. We derive a photometric redshift z = 6.3, estimated from the presence of the Lyman break falling between the I and J filters. This is by far the most distant GRB known to date. Its isotropic-equivalent energy is 3.4x10^53 erg in the rest-frame 110-1100 keV energy band. Despite the high redshift, both the prompt and the afterglow emission are not peculiar with respect to other GRBs. We find a break in the J-band light curve at t_b = 2.6 +- 1.0 d (observer frame). If we assume this is the jet break, we derive a beaming-corrected energy E_gamma = (4-12)x10^51 erg. This limit shows that GRB 050904 is consistent with the Amati and Ghirlanda relations. This detection is consistent with the expected number of GRBs at z > 6 and shows that GRBs are a powerful tool to study the star formation history up to very high redshift.Comment: 3 figures, 5 pages, accepted for publication in A&A Letters. One figure added, minor modifications. Full author list in the pape

Imaging neutral hydrogen on large-scales during the Epoch of Reionization with LOFAR

The first generation of redshifted 21 cm detection experiments, carried out with arrays like LOFAR, MWA and GMRT, will have a very low signal-to-noise ratio per resolution element (\sim 0.2). In addition, whereas the variance of the cosmological signal decreases on scales larger than the typical size of ionization bubbles, the variance of the formidable galactic foregrounds increases, making it hard to disentangle the two on such large scales. The poor sensitivity on small scales on the one hand, and the foregrounds effect on large scales on the other hand, make direct imaging of the Epoch of Reionization of the Universe very difficult, and detection of the signal therefore is expected to be statistical.Despite these hurdles, in this paper we argue that for many reionization scenarios low resolution images could be obtained from the expected data. This is because at the later stages of the process one still finds very large pockets of neutral regions in the IGM, reflecting the clustering of the large-scale structure, which stays strong up to scales of \sim 120 comoving Mpc/h (\sim 1 degree). The coherence of the emission on those scales allows us to reach sufficient S/N (\sim 3) so as to obtain reionization 21 cm images. Such images will be extremely valuable for answering many cosmological questions but above all they will be a very powerful tool to test our control of the systematics in the data. The existence of this typical scale (\sim 120 comoving Mpc/h) also argues for designing future EoR experiments, e.g., with SKA, with a field of view of at least 4 degree.Comment: Replaced with final version (minor changes), 9 figures, 11 pages, accepted for publication in MNRA

Early star-forming galaxies and the reionization of the Universe

Star forming galaxies represent a valuable tracer of cosmic history. Recent observational progress with Hubble Space Telescope has led to the discovery and study of the earliest-known galaxies corresponding to a period when the Universe was only ~800 million years old. Intense ultraviolet radiation from these early galaxies probably induced a major event in cosmic history: the reionization of intergalactic hydrogen. New techniques are being developed to understand the properties of these most distant galaxies and determine their influence on the evolution of the universe.Comment: Review article appearing in Nature. This posting reflects a submitted version of the review formatted by the authors, in accordance with Nature publication policies. For the official, published version of the review, please see http://www.nature.com/nature/archive/index.htm

The bright end of the z ~ 7 UV Luminosity Function from a wide and deep HAWK-I survey

(Abridged) We present here the second half of an ESO Large Programme, which exploits the unique combination of area and sensitivity provided in the near-IR by the camera Hawk-I at the VLT. We have obtained - 30 observing hours with Hawk-I in the Y-band of two high galactic latitude fields. We combined the Y-band data with deep J and K Hawk-I observations, and with FORS1/FORS2 U, B, V, R, I, and Z observations to select z-drop galaxies having Z - Y > 1, no optical detection and flat Y - J and Y - K colour terms. We detect 8 high-quality candidates in the magnitude range Y = 25.5 - 26.5 that we add to the z-drop candidates selected in two Hawk-I pointings over the GOODS-South field. We use this full sample of 15 objects found in -161 arcmin^2 of our survey to constrain the average physical properties and the evolution of the number density of z ~ 7 LBGs. A stacking analysis yields a best-fit SED with photometric redshift z= 6.85 +0.20 -0.15 and an E(B-V)=0.05 +0.15 -0.05. We compute a binned estimate of the z ~ 7 LF and explore the effects of photometric scatter and model uncertainties on the statistical constraints. After accounting for the expected incompleteness through MonteCarlo simulations, we strengthen our previous finding that a Schechter luminosity function constant from z=6 to z=7 is ruled out at a >99% confidence level, even including the effects of cosmic variance. For galaxies brighter than M_1500= -19.0, we derive a luminosity density rho_UV = 1.5^{+2.1}{-0.8} x 10^25 erg/s/Hz/Mpc^3, implying a decrease by a factor 3.5 from z=6 to z=6.8. We find that, under standard assumptions, the emission rate of ionizing photons coming from UV bright galaxies is lower by at least a factor of two than the value required for reionization. Finally, we exploit deep Hawk-I J and K band observations to derive an upper limit on the number density of M1500<~ -22.0 LBGs at z-8 (Y-dropouts).Comment: 12 pages, 8 figures. Accepted for publication in Astronomy & Astrophysic