Search

We use a set of three-dimensional (3D) radiative transfer simulations to study the effect that a large fraction of binary stars in galaxies during the epoch of reionization has on the physical properties of the intergalactic medium (i.e. the gas temperature and the ionization state of hydrogen and helium), on the topology of the ionized bubbles and on the 21 cm power spectra. Consistently to previous literature, we find that the inclusion of binary stars can speed up the reionization process of H  i and He  i , while He  ii reionization is still dominated by more energetic sources, especially accreting black holes. The earlier ionization attained with binary stars allows for more time for cooling and recombination so that gas fully ionized by binary stars is typically colder than gas ionized by single stars at any given redshift. With the same volume averaged ionization fraction, the inclusion of binary stars results in fewer small ionized bubbles and more large ones, with visible effects also on the large scales of the 21 cm power spectrum. [ABSTRACT FROM AUTHOR]

The presence of charged dust grains is known to have a profound impact on the physical evolution of the multiphase interstellar medium (ISM). Despite its importance, this process is still poorly explored in numerical simulations due to its complex physics and the tight dependence on the environment. Here, we introduce a novel implementation of grain charging in the cosmological radiative transfer code crash. We first benchmark the code predictions on a series of idealized dusty H  ii regions created by a single star, in order to assess the impact of grain properties on the resulting spatial distribution of charges. Secondly, we perform a realistic radiative transfer simulation of a star-forming region extracted from a dusty galaxy evolving in the Epoch of Reionization. We find that ∼13 per cent of the total dust mass gets negatively charged, mainly silicate and graphite grains of radius 10−3  |$\mu$| m. A complex spatial distribution of grain charges is also found, primarily depending on the exposure to stellar radiation and strongly varying along different lines of sight, as a result of radiative transfer effects. We finally assess the impact of grain properties (both chemical composition and size) on the resulting charge distribution. The new implementation described here will open up a wide range of possible studies investigating the physical evolution of the dusty ISM, nowadays accessible to observations of high- and low- redshift galaxies. [ABSTRACT FROM AUTHOR]

Using the latest upper limits on the 21-cm power spectrum at z ≈ 9.1 from the Low Frequency Array (LOFAR), we explore the regions of parameter space which are inconsistent with the data. We use 21cmmc , a Monte Carlo Markov chain sampler of 21cmfast which directly forward models the three dimensional (3D) cosmic 21-cm signal in a fully Bayesian framework. We use the astrophysical parametrization from 21cmfast , which includes mass-dependent star formation rates and ionizing escape fractions as well as soft-band X-ray luminosities to place limits on the properties of the high- z galaxies. Further, we connect the disfavoured regions of parameter space with existing observational constraints on the Epoch of Reionization such as ultra-violet (UV) luminosity functions, background UV photoionization rate, intergalactic medium (IGM) neutral fraction, and the electron scattering optical depth. We find that all models exceeding the 21-cm signal limits set by LOFAR at z ≈ 9.1 are excluded at ≳2 σ by other probes. Finally, we place limits on the IGM spin temperature from LOFAR, disfavouring at 95 per cent confidence spin temperatures below ∼2.6 K across an IGM neutral fraction range of |$0.15 \lesssim \bar{x}_{\rm H\, \rm {\small I}} \lesssim 0.6$|⁠. Note, these limits are only obtained from 141 h of data in a single redshift bin. With tighter upper limits, across multiple redshift bins expected in the near future from LOFAR, more viable models will be ruled out. Our approach demonstrates the potential of forward modelling tools such as 21cmmc in combining 21-cm observations with other high- z probes to constrain the astrophysics of galaxies. [ABSTRACT FROM AUTHOR]

The hyperfine transition of 3He+ at 3.5 cm has been thought as a probe of the high- z IGM, since it offers a unique insight into the evolution of the helium component of the gas, as well as potentially give an independent constraint on the 21 cm signal from neutral hydrogen. In this paper, we use radiative transfer simulations of reionization driven by sources such as stars, X-ray binaries, accreting black holes and shock heated interstellar medium, and simulations of a high- z quasar to characterize the signal and analyse its prospects of detection. We find that the peak of the signal lies in the range ∼1–50  μ K for both environments, but while around the quasar it is always in emission, in the case of cosmic reionization a brief period of absorption is expected. As the evolution of He  ii is determined by stars, we find that it is not possible to distinguish reionization histories driven by more energetic sources. On the other hand, while a bright QSO produces a signal in 21 cm that is very similar to the one from a large collection of galaxies, its signature in 3.5 cm is very peculiar and could be a powerful probe to identify the presence of the QSO. We analyse the prospects of the signal's detectability using SKA1-mid as our reference telescope. We find that the noise power spectrum dominates over the power spectrum of the signal, although a modest signal-to-noise ratio can be obtained when the wavenumber bin width and the survey volume are sufficiently large. [ABSTRACT FROM AUTHOR]

The Square Kilometre Array (SKA) 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 centred 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 z ≤ 10. 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 these redshifts. [ABSTRACT FROM AUTHOR]

We present a new dust extension to the Monte Carlo radiative transfer code crash, which enables it to simulate the propagation of ionizing radiation through mixtures of gas and dust. The new code is applied to study the impact of dust absorption on idealized galactic H  ii regions and on small-scale reionization. We find that H  ii regions are reduced in size by the presence of dust, while their inner temperature and ionization structure remain largely unaffected. In the small-scale reionization simulation, dust hardens ionization fronts and delays the overlap of ionized bubbles. This effect is found to depend only weakly on the assumed abundance of dust in underdense regions. [ABSTRACT FROM AUTHOR]

Recent results have suggested that active galactic nuclei (AGN) could provide enough photons to reionize the Universe. We assess the viability of this scenario using a semi-numerical framework for modelling reionization, to which we add a quasar contribution by constructing a Quasar Halo Occupancy Distribution (QHOD) based on Giallongo et al. observations. Assuming a constant QHOD, we find that an AGN-only model cannot simultaneously match observations of the optical depth τe, neutral fraction and ionizing emissivity. Such a model predicts τ e too low by ~2σ relative to Planck constraints, and reionizes the Universe at z ~ 5. Arbitrarily increasing the AGN emissivity to match these results yields a strong mismatch with the observed ionizing emissivity at z ~ 5. If we instead assume a redshift-independent AGN luminosity function yielding an emissivity evolution like that assumed in Madau & Haardt model, then we can match τe albeit with late reionization; however, such evolution is inconsistent with observations at z ~ 4-6 and poorly motivated physically. These results arise because AGN are more biased towards massive haloes than typical reionizing galaxies, resulting in stronger clustering and later formation times. AGN-dominated models produce larger ionizing bubbles that are reflected in ~×2 more 21 cm power on all scales. A model with equal part galaxies and AGN contribution is still (barely) consistent with observations, but could be distinguished using next-generation 21 cm experiments such as Hydrogen Epoch of Reionization Array and SKA-low. We conclude that, even with recent claims of more faint AGN than previously thought, AGN are highly unlikely to dominate the ionizing photon budget for reionization. [ABSTRACT FROM AUTHOR]

In an effort to understand the cause of the apparent depletion in the number density of radio-loud active galactic nuclei (AGNs) at z > 3, this work investigates the viability of the so-called cosmic microwave background (CMB) quenching mechanism of intrinsically jetted, high-z AGNs, whereby inverse Compton scattering of CMB photons off electrons within the extended lobes results in a substantial dimming of the lobe synchrotron emission at GHz frequencies, while simultaneously boosting their diffuse X-ray signal. We focus on five z > 3.5 radio galaxies that have sufficiently deep Chandra exposure (>50 ks) to warrant a meaningful investigation of any extended X-ray emission. For those objects with evidence for statistically significant extended X-ray lobes (4C 41.17 and 4C 03.24), we combine the Chandra measurements with literature data at lower frequencies to assemble the systems' spectral energy distributions (SEDs), and utilize state-of-the-art SED modelling - including emission from the disc, torus, jet, hotspots and lobes - to infer their physical parameters. For both radio galaxies, the magnetic energy density in the hotspots is found to exceed the energy density in CMB photons, whereas the opposite is true for the lobes. This implies that any extended synchrotron emission likely originates from the hotspots themselves, rather than the lobes. Conversely, inverse Compton scattering of CMB photons dominates the extended X-ray emission from the lobes, which are effectively 'radio-quenched'. As a result, CMB quenching is effective in these systems in spite of the fact that the observed X-ray to radio luminosity ratio does not bear the signature (1 + z)4 dependence of the CMB energy density. [ABSTRACT FROM AUTHOR]

We attempt to interpret existing data on the evolution of the UV luminosity function and UV colours, β, of galaxies at 5 ≤ z ≤ 8, to improve our understanding of their dust content and interstellar medium properties. To this aim, we post-process the results of a cosmological hydrodynamical simulation with a chemical evolution model, which includes dust formation by supernovae and intermediate-mass stars, dust destruction in supernova shocks, and grain growth by accretion of gas-phase elements in dense gas. We find that observations require a steep, Small Magellanic Cloud-like extinction curve and a clumpy dust distribution, where stellar populations younger than 15 Myr are still embedded in their dusty natal clouds. Investigating the scatter in the colour distribution and stellar mass, we find that the observed trends can be explained by the presence of two populations: younger, less massive galaxies where dust enrichment is mainly due to stellar sources, and massive, more chemically evolved ones, where efficient grain growth provides the dominant contribution to the total dust mass. Computing the IR-excess--UV colour relation, we find that all but the dustiest model galaxies follow a relation shallower than the Meurer et al. one, usually adopted to correct the observed UV luminosities of high-z galaxies for the effects of dust extinction. As a result, their total star formation rates might have been overestimated. Our study illustrates the importance to incorporate a proper treatment of dust in simulations of high-z galaxies, and that massive, dusty, UV-faint galaxies might have already appeared at z ... 7. [ABSTRACT FROM AUTHOR]