Lyman-alpha emitters gone missing: evidence for late reionization?

We combine high resolution hydrodynamical simulations with an intermediate resolution, dark matter only simulation and an analytical model for the growth of ionized regions to estimate the large scale distribution and redshift evolution of the visibility of Lyman-alpha emission in 6<=z<=8 galaxies. The inhomogeneous distribution of neutral hydrogen during the reionization process results in significant fluctuations in the Lyman-alpha transmissivity on large scales. The transmissivity depends not only on the ionized fraction of the intergalactic medium by volume and the amplitude of the local ionizing background, but is also rather sensitive to the evolution of the relative velocity shift of the Lyman-alpha emission line due to resonant scattering. We reproduce a decline in the space density of Lyman-alpha emitting galaxies as rapid as observed with a rather rapidly evolving neutral fraction between z=6-8, and a typical Lyman-alpha line velocity offset of 100 km/s redward of systemic at z=6 which decreases toward higher redshift. The new (02/2015) Planck results indicate such a recent end to reionization is no longer disfavoured by constraints from the cosmic microwave background.

Reply to [arXiv:1105.5653]: Comment on Quasinormal modes in Schwarzschild-de Sitter spacetime: A simple derivation of the level spacing of the frequencies

We explain why the analysis in our paper [Phys. Rev. D 69, 064033 (2004), arXiv:gr-qc/0311064 ] is relevant and correct.

Reionization constraints using Principal Component Analysis

Using a semi-analytical model developed by Choudhury & Ferrara (2005) we study the observational constraints on reionization via a principal component analysis (PCA). Assuming that reionization at z>6 is primarily driven by stellar sources, we decompose the unknown function N_{ion}(z), representing the number of photons in the IGM per baryon in collapsed objects, into its principal components and constrain the latter using the photoionization rate obtained from Ly-alpha forest Gunn-Peterson optical depth, the WMAP7 electron scattering optical depth and the redshift distribution of Lyman-limit systems at z sim 3.5. The main findings of our analysis are: (i) It is sufficient to model N_{ion}(z) over the redshift range 2<z<14 using 5 parameters to extract the maximum information contained within the data. (ii) All quantities related to reionization can be severely constrained for z<6 because of a large number of data points whereas constraints at z>6 are relatively loose. (iii) The weak constraints on N_{ion}(z) at z>6 do not allow to disentangle different feedback models with present data. There is a clear indication that N_{ion}(z) must increase at z>6, thus ruling out reionization by a single stellar population with non-evolving IMF, and/or star-forming efficiency, and/or photon escape fraction. The data allows for non-monotonic N_{ion}(z) which may contain sharp features around z sim 7. (iv) The PCA implies that reionization must be 99% completed between 5.8<z<10.3 (95% confidence level) and is expected to be 50% complete at z approx 9.5-12. With future data sets, like those obtained by Planck, the z>6 constraints will be significantly improved.

The impact of anisotropy from finite light travel time on detecting ionized bubbles in redshifted 21-cm maps

The detection of ionized bubbles around quasars in redshifted 21-cm maps is possibly one of the most direct future probes of reionization. We consider two models for the growth of spherical ionized bubbles to study the apparent shapes of the bubbles in redshifted 21-cm maps, taking into account the finite light travel time (FLTT) across the bubble. We find that the FLTT, whose effect is particularly pronounced for large bubbles, causes the bubbles image to continue to grow well after its actual growth is over. There are two distinct FLTT distortions in the bubbles image: (i) its apparent center is shifted along the line of sight (LOS) towards the observer from the quasar; (ii) its shape is anisotropic along the LOS. The bubble initially appears elongated along the LOS. This is reversed in the later stages of growth where the bubble appears compressed. The FLTT distortions are expected to have an impact on matched filter bubble detection where it is most convenient to use a spherical template for the filter. We find that the best matched spherical filter gives a reasonably good estimate of the size and the shift in the center of the anisotropic image. The mismatch between the spherical filter and the anisotropic image causes a 10 - 20% degradation in the SNR relative to that of a spherical bubble. We conclude that a spherical filter is adequate for bubble detection. The FLTT distortions do not effect the lower limits for bubble detection with 1000 hr of GMRT observations. The smallest spherical filter for which a detection is possible has comoving radii 24 Mpc and 33 Mpc for a 3-sigma and 5-sigma detection respectively, assuming a neutral fraction 0.6 at z sim 8.

Inside-out or Outside-in: The topology of reionization in the photon-starved regime suggested by Lyman-alpha forest data

We use a set of semi-numerical simulations based on Zeldovich approximation, friends-of-friends algorithm and excursion set formalism to generate reionization maps of high dynamic range with a range of assumptions regarding the distribution and luminosity of ionizing sources and the spatial distribution of sinks for the ionizing radiation. We find that ignoring the inhomogeneous spatial distribution of regions of high gas density where recombinations are important -- as is often done in studies of this kind -- can lead to misleading conclusions regarding the topology of reionization, especially if reionization occurs in the photon-starved regime suggested by Lya forest data. The inhomogeneous spatial distribution of recombinations significantly reduces the mean free path of ionizing photons and the typical size of coherently ionized regions. Reionization proceeds then much more as an outside-in process. Low-density regions far from ionizing sources become ionized before regions of high gas density not hosting sources of ionizing radiation. The spatial distribution of sinks of ionization radiation also significantly affects shape and amplitude the power spectrum of fluctuations of 21cm emission. The slope of the 21cm power spectrum as measured by upcoming 21cm experiments should be able to distinguish to what extent the topology of reionization proceeds outside-in or inside-out while the evolution of the amplitude of the power spectrum with increasing ionized mass fraction should be sensitive to the spatial distribution and the luminosity of ionizing sources.

On the minimum mass of reionization sources

By means of carefully calibrated semi-analytical reionization models, we estimate the minimum mass of star-forming haloes required to match the current data. Models which do not include haloes of total mass M < 10^9 M_sun fail at reproducing the Gunn-Peterson and electron scattering optical depths simultaneously, as they contribute too few (many) photons at high (low, z approx 6) redshift. Marginally acceptable solutions require haloes with M approx 5 times 10^7 M_sun at z approx 10, corresponding to virial temperatures (sim 10^4K) for which cooling can be ensured by atomic transitions. However, a much better match to the data is obtained if minihaloes (M sim 10^6 M_sun) are included in the analysis. We have critically examined the assumptions made in our model and conclude that reionization in the large-galaxies-only scenario can remain viable only if metal-free stars and/or some other exotic sources at z > 6 are included.