Search

The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts., 18 pages, 7 figures, v.2: Accepted for publications in MNRAS

This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources., 6 pages, 3 figures

An alternative to both the tomography technique and the power spectrum approach is to search for the 21cm forest, that is the 21cm absorption features against high-z radio loud sources caused by the intervening cold neutral intergalactic medium (IGM) and collapsed structures. Although the existence of high-z radio loud sources has not been confirmed yet, SKA-low would be the instrument of choice to find such sources as they are expected to have spectra steeper than their lower-z counterparts. Since the strongest absorption features arise from small scale structures (few tens of physical kpc, or even lower), the 21cm forest can probe the HI density power spectrum on small scales not amenable to measurements by any other means. Also, it can be a unique probe of the heating process and the thermal history of the early universe, as the signal is strongly dependent on the IGM temperature. Here we show what SKA1-low could do in terms of detecting the 21cm forest in the redshift range z = 7.5-15., Accepted for publication in the SKA Science Book 'Advancing Astrophysics with the Square Kilometre Array', to appear in 2015; 10 pages, 5 figures; the manuscript is based on Ciardi et al., 2013, MNRAS, 428, 1755

We interpret recent ALMA observations of z > 6 normal star forming galaxies by means of a semi-numerical method, which couples the output of a cosmological hydrodynamical simulation with a chemical evolution model which accounts for the contribution to dust enrichment from supernovae, asymptotic giant branch stars and grain growth in the interstellar medium. We find that while stellar sources dominate the dust mass of small galaxies, the higher level of metal enrichment experienced by galaxies with Mstar > 10^9 Msun allows efficient grain growth, which provides the dominant contribution to the dust mass. Even assuming maximally efficient supernova dust production, the observed dust mass of the z = 7.5 galaxy A1689-zD1 requires very efficient grain growth. This, in turn, implies that in this galaxy the average density of the cold and dense gas, where grain growth occurs, is comparable to that inferred from observations of QSO host galaxies at similar redshifts. Although plausible, the upper limits on the dust continuum emission of galaxies at 6.5 < z < 7.5 show that these conditions must not apply to the bulk of the high redshift galaxy population, 5 pages, 2 figures, accepted for publication as a Letter to MNRAS

We perform a systematic study of how the inhomogeneities in the Inter-Galactic Medium (IGM) affect the observability of Lyman Alpha Emitters (LAEs) around the Epoch of Reionization. We focus on the IGM close to the galaxies as the detailed ionization distribution and velocity fields of this region could significantly influence the scattering of Ly-alpha photons off neutral H atoms as they traverse the IGM after escaping from the galaxy. We simulate the surface brightness (SB) maps and spectra of more than 100 LAEs at z=7.7 as seen by an observer at z=0. To achieve this, we extract the source properties of galaxies and their surrounding IGM from cosmological simulations of box sizes 5-30 Mpc/h and follow the coupled radiative transfer of ionizing and Ly-alpha radiation through the IGM using CRASH-alpha. We find that the simulated SB profiles are extended and their detailed structure is affected by inhomogeneities in the IGM, especially at high neutral fractions. The detectability of LAEs and the fraction of the flux observed depend heavily on the shape of the SB profile and the SB threshold (SB_th) of the observational campaign. Only ultradeep observations (e.g. SB_th ~ 10^-23 ergs/s/cm^2/arcsec^2) would be able to obtain the true underlying mass-luminosity relation and luminosity functions of LAEs. The details of our results depend on whether Ly-alpha photons are significantly shifted in the galaxy to longer wavelengths, the mean ionization fraction in the IGM and the clustering of ionizing sources. These effects can lead to an easier escape of Ly-alpha photons with less scattering in the IGM and a concentrated SB profile similar to the one of a point source. Finally, we show that the SB profiles are steeper at high ionization fraction for the same LAE sample which can potentially be observed from the stacked profile of a large number of LAEs., Comment: 22 pages, 23 figures, 2 tables, Accepted by MNRAS. Minor changes

We propose to use high-redshift long $\gamma$-ray bursts (GRBs) as cosmological tools to constrain the amount of primordial non-Gaussianity in the density field. By using numerical, N-body, hydrodynamic, chemistry simulations of different cosmological volumes with various Gaussian and non-Gaussian models, we self-consistently relate the cosmic star formation rate density to the corresponding GRB rate. Assuming that GRBs are fair tracers of cosmic star formation, we find that positive local non-Gaussianities, described in terms of the non-linear parameter, \fnl, might boost significantly the GRB rate at high redshift, $z \gg 6$. Deviations with respect to the Gaussian case account for a few orders of magnitude if \fnl$\sim 1000$, one order of magnitude for \fnl$\sim 100$, and a factor of $\sim 2$ for \fnl$\sim 50$. These differences are found only at large redshift, while at later times the rates tend to converge. Furthermore, a comparison between our predictions and the observed GRB data at $z > 6$ allows to exclude large negative \fnl, consistently with previous works. Future detections of any long GRB at extremely high redshift ($z\sim 15-20$) could favor non-Gaussian scenarios with positive \fnl. More stringent constraints require much larger high-$z$ GRB complete samples, currently not available in literature. By distinguishing the contributions to the GRB rate from the metal-poor population III regime, and the metal-enriched population II-I regime, we conclude that the latter is a more solid tracer of the underlying matter distribution, while the former is strongly dominated by feedback mechanisms from the first, massive, short-lived stars, rather than by possible non-Gaussian fluctuations. This holds quite independently of the assumed population III initial mass function., Comment: 12 pages; MNRAS in press. Chi ha paura muore tutti i giorni

We present the new parallel version (pCRASH2) of the cosmological radiative transfer code CRASH2 for distributed memory supercomputing facilities. The code is based on a static domain decomposition strategy inspired by geometric dilution of photons in the optical thin case that ensures a favourable performance speed-up with increasing number of computational cores. Linear speed-up is ensured as long as the number of radiation sources is equal to the number of computational cores or larger. The propagation of rays is segmented and rays are only propagated through one sub-domain per time step to guarantee an optimal balance between communication and computation. We have extensively checked pCRASH2 with a standardised set of test cases to validate the parallelisation scheme. The parallel version of CRASH2 can easily handle the propagation of radiation from a large number of sources and is ready for the extension of the ionisation network to species other than hydrogen and helium., Comment: Accepted by Monthly Notices of the Royal Astronomical Society. 18 pages, 20 figures

We discuss expectations for the absorption of high-energy gamma-rays by gamma-gamma pair production with intergalactic radiation fields (IRFs) at very high redshifts (z~5-20), and the prospects thereof for probing the cosmic reionization era. For the evolving IRF, a semi-analytical model incorporating both Population II and Population III stars is employed, which is consistent with a wide variety of existing high-z observations including QSO spectral measurements, WMAP Thomson depth constraints, near-IR source count limits, etc. We find that the UV IRF below the Lyman edge energy with intensities in the range of a few times 10^{-19} erg cm^{-2} s^{-1} Hz^{-1} sr^{-1} can cause appreciable attenuation above ~12 GeV at z~5, down to ~6-8 GeV at z>~8-10. This may be observable in the spectra of blazars or gamma-ray bursts by the Fermi Gamma-ray Space Telescope or next generation facilities such as the Cherenkov Telescope Array, Advanced Gamma-ray Imaging System or 5@5, providing invaluable insight into early star formation and cosmic reionization., MNRAS in press with minor revisions, 5 pages, 5 figures. Numerical data of the model results will be available at http://www-tap.scphys.kyoto-u.ac.jp/~inoue/hizabs/

The morphology of the total gas and the neutral hydrogen (H I) distributions during the cosmological epoch of reionization can be quantified with Minkowski Functionals (MFs) of isodensity surfaces. We compute the MFs from the output of a high-resolution numerical simulation which includes explicit treatment of transfer of UV ionizing radiation. "Galaxies" identified in the simulation using semi-analytic models of galaxy formation are assumed to be the sole sources of UV photons. The MFs of the total gas distribution are well described by the analytic expressions derived for lognormal random fields. The statistical properties of the diffuse H I depend on the gas distribution and on the way ionized regions propagate in the inter-galactic medium (IGM). The deviations of the MFs of the H I distribution from those of a lognormal random field are, therefore, caused by reionization. We use the MFs to discriminate between the various stages of reionization in the simulation. We suggest a simple model of reionization which reproduces the MFs derived from this simulation. Using random realizations of lognormal density fields, we also assess the ability of MFs to distinguish between different reionization scenarios. Our results are relevant to the analysis of future cosmological twenty-one centimeter maps., Comment: 12 pages, 8 figures; Minor content changes. To be published in MNRAS

We calculate the expected Mid-IR molecular hydrogen line emission from the first objects in the universe. Due to their low masses, the stellar feedback from massive stars is able to blowaway their gas content and collect it into a cooling shell where H_2 rapidly forms and IR roto-vibrational (as for example the restframe 2.12 micron) lines carry away a large fraction (up to 10%) of the explosion energy. The fluxes from these sources are in the range 10^{-21}-10^{-17} erg/s/cm^2. The highest number counts are expected in the 20 micron band, where about 10^5 sources/deg^2 are predicted at the limiting flux of 3 \times 10^{-18} erg/s/cm^2. Among the planned observational facilities, we find that the best detection perspectives are offered by NGST, which should be able to reveal about 200 first objects in one hour observation time at its limiting flux in the above band. Therefore, Mid-IR instruments appear to represent perfect tools to trace star formation and stellar feedback in the high (z>5) redshift universe., 6 pages, 4 figures. Revised version, accepted for publication by MNRAS

A number of experiments are set to measure the 21-cm signal of neutral hydrogen from the Epoch of Reionization (EoR). The common denominator of these experiments are the large data sets produced, contaminated by various instrumental effects, ionospheric distortions, RFI and strong Galactic and extragalactic foregrounds. In this paper, the first in a series, we present the Data Model that will be the basis of the signal analysis for the LOFAR (Low Frequency Array) EoR Key Science Project (LOFAR EoR KSP). Using this data model we simulate realistic visibility data sets over a wide frequency band, taking properly into account all currently known instrumental corruptions (e.g. direction-dependent gains, complex gains, polarization effects, noise, etc). We then apply primary calibration errors to the data in a statistical sense, assuming that the calibration errors are random Gaussian variates at a level consistent with our current knowledge based on observations with the LOFAR Core Station 1. Our aim is to demonstrate how the systematics of an interferometric measurement affect the quality of the calibrated data, how errors correlate and propagate, and in the long run how this can lead to new calibration strategies. We present results of these simulations and the inversion process and extraction procedure. We also discuss some general properties of the coherency matrix and Jones formalism that might prove useful in solving the calibration problem of aperture synthesis arrays. We conclude that even in the presence of realistic noise and instrumental errors, the statistical signature of the EoR signal can be detected by LOFAR with relatively small errors. A detailed study of the statistical properties of our data model and more complex instrumental models will be considered in the future., 20 pages, 9 figures

The visibility of LyA emitting galaxies during the Epoch of Reionization is controlled by both diffuse HI patches in large-scale bubble morphology and small-scale absorbers. To investigate the impact on LyA photons, we apply a novel combination of analytic and numerical calculations to three scenarios: (i) the `bubble' model, where only diffuse HI outside ionized bubbles is present; (ii) the `web' model, where HI exists only in overdense self-shielded gas; and (iii) the more realistic 'web-bubble' model, which contains both. Our analysis confirms that there is a degeneracy between the ionization structure of the intergalactic medium (IGM) and the HI fraction inferred from LyA surveys, as the three models suppress LyA flux equally with very different HI fractions. We argue that a joint analysis of the LyA luminosity function and the rest-frame equivalent width distribution/LyA fraction can break this degeneracy and provide constraints on the reionization history and its topology. We further show that constraints can improve if we consider the full shape of the M_UV-dependent redshift evolution of the LyA fraction of Lyman break galaxies. Contrary to conventional wisdom, we find that (i) a drop of LyA fraction larger for UV-faint than for UV-bright galaxies can be reproduced with web and web-bubble models and therefore does not provide exclusive evidence of patchy reionization, and (ii) the IGM-transmission PDF is unimodal for bubble models and bimodal in web models. We further highlight the importance of galaxy-absorber cross-correlation. Comparing our models to observations, the neutral fraction at z~7 is likely to be of order of tens of per cent when interpreted with bubble or web-bubble models. Alternatively, we obtain a conservative lower limit ~1% in the web models, if we allow for a drop in the photoionization rate by a factor of ~100 from the post-reionized universe. [abridged], Comment: 23 pages, 20 figures, submitted to MNRAS, Abstract abridged for arXiv submission

Estimating the intergalactic medium ionization level of a region needs proper treatment of the reionization process for a large representative volume of the universe. The clumping factor, a parameter which accounts for the effect of recombinations in unresolved, small-scale structures, aids in achieving the required accuracy for the reionization history even in simulations with low spatial resolution. In this paper, we study for the first time the redshift evolution of clumping factors of different ionized species of H and He in a small but very high resolution simulation of the reionization process. We investigate the dependence of the value and redshift evolution of clumping factors on their definition, the ionization level of the gas, the grid resolution, box size and mean dimensionless density of the simulations., Comment: 12 pages, 10 figures, 1 table. Accepted by MNRAS

Radiative transfer simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come to the reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray-tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non-trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on 5 test problems: (0) basic physics, (1) isothermal H II region expansion and (2) H II region expansion with evolving temperature, (3) I-front trapping and shadowing by a dense clump, (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analyzed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multi-frequency treatment approach, resulting in different thicknesses of the ionized-neutral transition regions and different temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To this aim all test inputs and outputs are made publicly available in digital form., Comment: 32 pages, 39 figures (all color), comments welcome