ABSTRACT
We measure escape fractions, fesc, of ionizing radiation from galaxies in the sphinx suite of cosmological radiation-hydrodynamical simulations of reionization, resolving haloes with $M_{\rm ...vir}\gtrsim 7.5 \times 10^7 \ {\rm {M}_{\odot }}$ with a minimum cell width of ≈10 pc. Our new and largest 20 co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with halo masses up to a few times 1011 M⊙. The simulated galaxies agree well with observational constraints of the ultraviolet (UV) luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over time-scales of a few Myr, due to its regulation by supernova and radiation feedback, and at any given time a tiny fraction of star-forming galaxies emits a large fraction of the ionizing radiation escaping into the intergalactic medium. Statistically, fesc peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies (M* ≈ 107 M⊙, M1500 ≈ −17, Z ≲ 5 × 10−3 Z⊙), dropping strongly for lower and higher masses, brighter and dimmer galaxies, and more metal-rich galaxies. The escape fraction correlates positively with both the short-term and long-term specific star formation rate. According to sphinx, galaxies too dim to be yet observed, with ${M_{1500}}\gtrsim -17$, provide about 55 per cent of the photons contributing to reionization. The global averaged fesc naturally decreases with decreasing redshift, as predicted by UV background models and low-redshift observations. This evolution is driven by decreasing specific star formation rates over cosmic time.
ABSTRACT
A long (110 cMpc h−1) and deep absorption trough in the Ly α forest has been observed extending down to redshift 5.5 in the spectrum of ULAS J0148+0600. Although no Ly α transmission is ...detected, Ly β spikes are present which has led to claims that the gas along this trough must be ionized. Using high-resolution cosmological radiative transfer simulations in large volumes, we show that in a scenario where reionization ends late (z ∼ 5.2), our simulations can reproduce troughs as long as observed. In this model, we find that the troughs are caused by islands of neutral hydrogen. Small ionized holes within the neutral islands allow for the transmission of Ly β. We have also modelled the Ly α emitter population around the simulated troughs, and show that there is a deficit of Ly α emitters close to the trough as is observed.
Abstract
Fast radio bursts offer the opportunity to place new constraints on the mass and density profile of hot and ionized gas in galactic haloes. We test here the X-ray emission and dispersion ...measure predicted by different gas profiles for the halo of the Milky Way. We examine a range of models, including entropy stability conditions and external pressure continuity. We find that incorporating constraints from X-ray observations leads to favouring dispersion measures on the lower end of the range given by these models. We show that the dispersion measure of the Milky Way halo could be less than 10 cm−3 pc in the most extreme model we consider, which is based on constraints from O vii absorption lines. However, the models allowed by the soft X-ray constraints span more than an order of magnitude in dispersion measures. Additional information on the distribution of gas in the Milky Way halo could be obtained from the signature of a dipole in the dispersion measure of fast radio bursts across the sky, but this will be a small effect for most cases.
Abstract
The thermal state of the post-reionization IGM is sensitive to the timing of reionization and the nature of the ionizing sources. We have modelled here the thermal state of the IGM in ...cosmological radiative transfer simulations of a realistic, extended, spatially inhomogeneous hydrogen reionization process, carefully calibrated with Ly α forest data. We compare these with cosmological simulations run using a spatially homogeneous ionizing background. The simulations with a realistic growth of ionized regions and a realistic spread in reionization redshifts show, as expected, significant spatial fluctuations in the temperature–density relation (TDR) of the post-reionization IGM. The most recently ionized regions are hottest and exhibit a flatter TDR. In simulations consistent with the average TDR inferred from Ly α forest data, these spatial fluctuations have a moderate but noticeable effect on the statistical properties of the Ly α opacity of the IGM at z ∼ 4–6. This should be taken into account in accurate measurements of the thermal properties of the IGM and the free-streaming of dark matter from Ly α forest data in this redshift range. The spatial variations of the TDR predicted by our simulations are, however, smaller by about a factor of 2 than would be necessary to explain the observed large spatial opacity fluctuations on large (≥50 h−1 comoving Mpc) scales atz ≳ 5.5.
Abstract High-redshift quasi-stellar object (QSO) spectra show large spatial fluctuations in the Ly α opacity of the intergalactic medium on surprisingly large scales at $z$ ≳ 5.5. We present a ...radiative transfer simulation of cosmic reionization driven by galaxies that reproduces this large scatter and the rapid evolution of the Ly α opacity distribution at 5 < $z$ < 6. The simulation also reproduces the low Thomson scattering optical depth reported by the latest cosmic microwave background (CMB) measurement and is consistent with the observed short near-zones and strong red damping wings in the highest redshift QSOs. It also matches the rapid disappearance of observed Ly α emission by galaxies at $z$ ≳ 6. Reionization is complete at $z$ = 5.3 in our model, and 50 per cent of the volume of the Universe is ionized at $z$ = 7. Agreement with the Ly α forest data in such a late reionization model requires a rapid evolution of the ionizing emissivity of galaxies that peaks at $z$ ∼ 6.8. The late end of reionization results in a large scatter in the photoionization rate and the neutral hydrogen fraction at redshifts as low as $z$ ≲ 5.5 with large residual neutral ‘islands’ that can produce very long Gunn–Peterson troughs resembling those seen in the data.
ABSTRACT
Proximity zones of high-redshift quasars are unique probes of supermassive black hole formation, but simultaneously explaining proximity zone sizes and black hole masses has proved to be ...challenging. We study the robustness of some of the assumptions that are usually made to infer quasar lifetimes from proximity zone sizes. We show that small proximity zones can be readily explained by quasars that vary in brightness with a short duty cycle of fduty ∼ 0.1 and short bright periods of ton ∼ 104 yr, even for long lifetimes. We further show that reconciling this with black hole mass estimates requires the black hole to continue to grow and accrete during its obscured phase. The consequent obscured fractions of ≳ 0.7 or higher are consistent with low-redshift measurements and models of black hole accretion. Such short duty cycles and long obscured phases are also consistent with observations of large proximity zones, thus providing a simple, unified model for proximity zones of all sizes. The large dynamic range of our simulation, and its calibration to the Lyman-α forest, allows us to investigate the influence of the large-scale topology of reionization and the quasar’s host halo mass on proximity zones. We find that incomplete reionization can impede the growth of proximity zones and make them smaller up to 30 per cent, but the quasar host halo mass only affects proximity zones weakly and indirectly. Our work suggests that high-redshift proximity zones can be an effective tool to study quasar variability and black hole growth.
ABSTRACT
Our understanding of the intergalactic medium at redshifts z = 5–6 has improved considerably in the last few years due to the discovery of quasars with z > 6 that enable Lyman-α forest ...studies at these redshifts. A realization from this has been that hydrogen reionization could end much later than previously thought, so that large ‘islands’ of cold, neutral hydrogen could exist in the IGM at redshifts z = 5–6. By using radiative transfer simulations of the IGM, we consider the implications of the presence of these neutral hydrogen islands for the 21-cm power spectrum signal and its potential detection by experiments such as hera, ska, lofar, and mwa. In contrast with previous models of the 21-cm signal, we find that thanks to the late end of reionization the 21-cm power in our simulation continues to be as high as $\Delta ^2_{21}=10~\mathrm{mK}^2$ at k ∼ 0.1 h cMpc−1 at z = 5–6. This value of the power spectrum is several orders of magnitude higher than that in conventional models considered in the literature for these redshifts. Such high values of the 21-cm power spectrum should be detectable by hera and ska1-low in ∼1000 h, assuming optimistic foreground subtraction. This redshift range is also attractive due to relatively low sky temperature and potentially greater abundance of multiwavelength data.
Abstract
A cosmological zoom-in simulation that develops into a Milky Way-like halo begins at redshift 7. The initial dark matter distribution is seeded with dense star clusters of median mass 5 × 10
...5
M
⊙
, placed in the largest subhalos present, which have a median peak circular velocity of 25 km s
−1
. Three simulations are initialized using the same dark matter distribution with the star clusters starting on approximately circular orbits having initial median radii 6.8, 0.14 kpc, and, at the exact center of the subhalos. The simulations are evolved to the current epoch at which time the median galactic orbital radii of the three sets of clusters are 30, 5, and 16 kpc, with the clusters losing about 2%, 50%, and 15% of their mass, respectively. Clusters beginning at small orbital radii have so much tidal forcing that they are often not in equilibrium. Clusters that start at larger subhalo radii have a velocity dispersion that declines smoothly to ≃20% of the central value at ≃20 half-mass radii. The clusters that begin in the subhalo centers can show a rise in velocity dispersion beyond 3–5 half-mass radii. That is, the clusters that form without local dark matter always have stellar-mass-dominated kinematics at all radii, whereas about 25% of the clusters that begin in subhalo centers have remnant local dark matter.
QSO near zones are an important probe of the ionization state of the intergalactic medium (IGM) at z ∼ 6–7, at the end of reionization. We present here high-resolution cosmological 3D radiative ...transfer simulations of QSO environments for a wide range of host halo masses, 1010–12.5 M⊙. Our simulated near zones reproduce both the overall decrease of observed near-zone sizes at 6 < z < 7 and their scatter. The observable near-zone properties in our simulations depend only very weakly on the mass of the host halo. The size of the H ii region expanding into the IGM is generally limited by (super-)Lyman Limit systems loosely associated with (low-mass) dark matter haloes. This leads to a strong dependence of near-zone size on direction and drives the large observed scatter. In the simulation centred on our most massive host halo, many sightlines show strong red damping wings even for initial volume averaged neutral hydrogen fractions as low as ∼10−3. For QSO lifetimes long enough to allow growth of the central supermassive black hole while optically bright, we can reproduce the observed near zone of ULAS J1120+0641 only with an IGM that is initially neutral. Our results suggest that larger samples of z > 7 QSOs will provide important constraints on the evolution of the neutral hydrogen fraction and thus on how late reionization ends.
ABSTRACT
We compare a sample of five high-resolution, high S/N Ly α forest spectra of bright 6 < z < ∼6.5 QSOs aimed at spectrally resolving the last remaining transmission spikes at z > 5 with ...those obtained from mock absorption spectra from the Sherwoodand Sherwood–Relics simulation suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile-fitting procedure for the inverted transmitted flux, 1 − F, similar to the widely used Voigt profile fitting of the transmitted flux F at lower redshifts, to characterize the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominent in low temperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution (≤ 8 km s−1) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of T0 = (11 000 ± 1600, 10 500 ± 2100, 12 000 ± 2200) K at z = (5.4, 5.6, 5.8). This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to z ∼ 5.3, the width distribution of the observed transmission spikes is consistent with the range of T0 caused by spatial fluctuations in the temperature–density relation.