Early Reionization Research Articles

Accelerated Structure Formation: The Early Emergence of Massive Galaxies and Clusters of Galaxies

Abstract Galaxies in the early Universe appear to have grown too big too fast, assembling into massive, monolithic objects more rapidly than anticipated in the hierarchical Lambda cold dark matter (ΛCDM) structure formation paradigm. The available photometric data are consistent with there being a population of massive galaxies that form early (z ≳ 10) and quench rapidly over a short (≲1 Gyr) timescale, consistent with the traditional picture for the evolution of giant elliptical galaxies. Similarly, kinematic observations as a function of redshift show that massive spirals and their scaling relations were in place at early times. Explaining the early emergence of massive galaxies requires either an extremely efficient conversion of baryons into stars at z > 10 or a more rapid assembly of baryons than anticipated in ΛCDM. The latter possibility was explicitly predicted in advance by modified Newtonian dynamics (MOND). We discuss some further predictions of MOND, such as the early emergence of clusters of galaxies and early reionization.

An Estimate of the Impact of Reionization on Supermassive Black Hole Growth

The supermassive black holes (SMBHs) that power active galactic nuclei found at z ≥ 6 were formed during the Epoch of Reionization. Because reionization is an inhomogeneous process, the physical properties of SMBH host-galaxy environments will vary spatially during reionization. We construct a semi-analytic model to estimate the impact of reionization on SMBH growth. Using a series of merger trees, reionization models, and black hole (BH) growth models, we find that early reionization can reduce a SMBH’s mass by up to [50, 70, 90]% within dark matter halos of mass [1012, 1011, 1010] M ⊙ by z = 6. Our findings also suggest that the redshift range in which BH growth is impacted by reionization strongly depends on whether the Eddington accretion rate can be exceeded. If so, we find that BH masses are significantly suppressed principally during the early phases of reionization (z ≳ 10), while they are more readily suppressed across the full redshift range if super-Eddington growth is not allowed. We find that the global average impact of reionization may be to reduce the masses of BHs residing in ≲1011 M ⊙ halos by a factor of ≳2. The census of SMBHs being uncovered by the JWST may offer a means to test the basic prediction that more massive BHs reside in cosmological volumes that are reionized at later times.

SPICE: the connection between cosmic reionization and stellar feedback in the first galaxies

ABSTRACT We present SPICE, a new suite of radiation-hydrodynamic, cosmological simulations targeting the epoch of reionization. The goal of these simulations is to systematically probe a variety of stellar feedback models, including ‘bursty’ and ‘smooth’ forms of supernova energy injection, as well as poorly explored physical scenarios such as hypernova explosions and radiation pressure on dust. We show that even subtle differences in the behaviour of supernova feedback drive profound differences in reionization histories, with burstier forms of feedback causing earlier reionization. However, we also find that some global galaxy properties, such as the dust-attenuated luminosity functions and star formation main sequence, remain degenerate between models. In particular, we show that stellar feedback and its strength determine the morphological mix of galaxies emerging by $z \, = \, 5$ and that the reionization history is inextricably connected to intrinsic properties such as galaxy kinematics and morphology. While star-forming, massive discs are prevalent if supernova feedback is ‘smooth’, ‘bursty’ feedback preferentially generates dispersion-dominated systems. Different modes of feedback produce different strengths of outflows, altering the interstellar/circumgalactic medium in different ways, and in turn strongly affecting the escape of Lyman continuum (LyC) photons. We establish a correlation between galaxy morphology and LyC escape fraction, revealing that dispersion-dominated systems have escape fractions 10–50 times higher than their rotation-dominated counterparts at all redshifts. At the same intrinsic luminosity, dispersion-dominated systems should thus preferentially generate large H ii regions as compared to their rotation-dominated counterparts. Since dispersion-dominated systems are more prevalent if stellar feedback is more explosive, reionization occurs earlier in our simulation with burstier feedback. We argue that statistical samples of post-reionization galaxy morphologies (using both stellar and gaseous components) probed with telescopes such as JWST, ALMA, and MUSE can constrain stellar feedback at z > 5 and models of cosmic reionization.

Constraints on the Epoch of Reionization with Roman Space Telescope and the Void Probability Function of Lyα Emitters

We use large simulations of Lyα emitters with different fractions of ionized intergalactic medium to quantify the clustering of Lyα emitters as measured by the void probability function (VPF), and how it evolves under different ionization scenarios. We quantify how well we might be able to distinguish between these scenarios with a deep spectroscopic survey using the future Nancy Grace Roman Space Telescope. Because Roman will be able to carry out blind spectroscopic surveys of Lyα emitters continuously between 7 < z < 12 to sensitivities of at least 10−17 erg s−1 over a wide field of view, it can measure the epoch of reionization as well as the pace of ionization of the intergalactic medium (IGM). We compare deep Roman surveys covering roughly 1, 4, and 16 deg2, and quantify what constraints on reionization the VPF may find for these surveys. A survey of 1 deg2 would distinguish between very late reionization and early reionization to 3σ near z = 7.7 with the VPF. The VPF of a 4 deg2 survey can distinguish between slow versus fast, as well as early versus late, reionization at >3σ–4σ at several redshifts between 7 < z < 9. However, a survey of 13–16 deg2 would allow the VPF to give several robust constraints (>5–8σ) across the epoch of reionization, and would yield a detailed history of the reionization of the IGM and its effect on Lyα emitter clustering.

Gravitational Wave Backgrounds from Coalescing Black Hole Binaries at Cosmic Dawn: An Upper Bound

Abstract The successive discoveries of binary merger events by Advanced LIGO-Virgo have been revealing the statistical properties of binary black hole (BBH) populations. A stochastic gravitational wave background (GWB) is a useful tool to probe the cosmological evolution of those compact mergers. In this paper, we study the upper bound on a GWB produced by BBH mergers, whose stellar progenitors dominate the reionization process at the cosmic dawn. Since early reionization by those progenitors yields a high optical depth of the universe inconsistent with the Planck measurements, the cumulative mass density is limited to ρ ⋆ ≲ 107 M ⊙ Mpc−3. Even with this upper bound, the amplitude of a GWB owing to the high-z BBH mergers is expected to be as high as Ω gw ≃ 1.48 − 1.27 + 1.80 × 10 − 9 at f ≃ 25 Hz, while their merger rate at the present day is consistent or lower than the observed GW event rate. This level of GWB is detectable at the design sensitivity of Advanced LIGO-Virgo and would indicate a major contribution of the high-z BBH population to the local GW events. The spectral index is expected to be substantially flatter than the canonical value of ≃2/3 generically produced by lower-redshift and less-massive BBHs. Moreover, if their mass function is more top heavy than in the local universe, the GWB spectrum is even more skewed toward lower frequencies, which would allow us to extract information on the mass function of merging BBHs at high redshifts.

Stacking redshifted 21 cm images of H ii regions around high-redshift galaxies as a probe of early reionization

ABSTRACT A number of current and future experiments aim to detect the reionization of neutral hydrogen by the first stars and galaxies in the Universe via the redshifted 21 cm line. Using the bluetides simulation, we investigate the measurement of an average ionized region towards the beginning of reionization by stacking redshifted 21 cm images around optically identified bright galaxies using mock observations. We find that with an SKA 1000 h observation, assuming perfect foreground subtraction, a 5σ detection of a stacked H ii region can be made with 30 images around some of the brightest galaxies in bluetides (brighter than MUV &amp;lt; −22.75) at z = 9 (corresponding to a neutral fraction of 90.1 per cent in our model). We present simulated relationships between the UV magnitude of galaxies, the sizes of the ionized regions they reside in, and the shape of the stacked profiles. These mock observations can also distinguish between scenarios where the intergalactic medium is in net emission or absorption of 21 cm photons. Once 21 cm foreground contamination is included, we find that even with up to 200 images around these rare, bright galaxies, only a tentative &amp;gt;1σ detection will be possible. However, partial foreground subtraction substantially improves signal to noise. For example, we predict that reducing the area of Fourier space dominated by foregrounds by 50 (80) per cent will allow &amp;gt;3σ (&amp;gt;5σ) detections of ionized regions at z = 9.

Constraining the primordial magnetic field with dwarf galaxy simulations

Using a set of cosmological hydro-dynamical simulations, we constrained the properties of primordial magnetic fields by studying their impact on the formation and evolution of dwarf galaxies. We performed a large set of simulations (8 dark matter only and 72 chemo-hydrodynamical) including primordial magnetic fields through the extra density fluctuations they induce at small length scales (k ≥ 10 h Mpc−1) in the matter power spectrum. Our sample of dwarfs includes nine systems selected out of the initial (3.4 Mpc h−1)3 parent box, resimulated from z = 200 to z = 0 using a zoom-in technique and including the physics of baryons. We explored a wide variety of primordial magnetic fields with strength Bλ ranging from 0.05 to 0.50 nG and magnetic energy spectrum slopes nB from −2.9 to −2.1. Strong magnetic fields characterized by a high amplitude (Bλ = 0.50, 0.20 nG with nB = −2.9) or by a steep initial power spectrum slope (nB = −2.1, −2.4, with Bλ = 0.05 nG) induce perturbations on mass scales from 107 to 109 M⊙. In this context emerging galaxies see their star formation rates strongly boosted. They become more luminous and metal rich than their counterparts without primordial magnetic fields. Such strong fields are ruled out by their inability to reproduce the observed scaling relations of dwarf galaxies. They predict that dwarf galaxies are at the origin of an unrealistically early reionization of the Universe and that they also overproduce luminous satellites in the Local Group. Weaker magnetic fields impacting the primordial density field at corresponding masses ≲106 M⊙, produce a large number of mini dark matter halos orbiting the dwarfs, however out of reach for current lensing observations. This study allows us, for the first time, to constrain the properties of primordial magnetic fields based on realistic cosmological simulations of dwarf galaxies.

Fully Coupled Simulation of Cosmic Reionization. III. Stochastic Early Reionization by the Smallest Galaxies

Abstract Previously we identified a new class of early galaxy that we estimate contributes up to 30% of the ionizing photons responsible for reionization. These are low-mass halos in the range = 106.5–108 that have been chemically enriched by supernova ejecta from prior Population III star formation. Despite their low star formation rates, these metal cooling halos (MCs) are significant sources of ionizing radiation, especially at the onset of reionization, owing to their high number density and ionizing escape fractions. Here we present a fully coupled radiation hydrodynamic simulation of reionization that includes these MCs, as well the more massive hydrogen atomic line cooling halos. Our method is novel: we perform halo finding in line with the radiation hydrodynamical simulation and assign escaping ionizing fluxes to halos using a probability distribution function (pdf) measured from the galaxy-resolving Renaissance Simulations. The pdf captures the mass dependence of the ionizing escape fraction, as well as the probability that a halo is actively forming stars. With MCs, reionization starts earlier than if only halos of 108 and above are included; however, the redshift when reionization completes is only marginally affected, as this is driven by more massive galaxies. Because star formation is intermittent in MCs, the earliest phase of reionization exhibits a stochastic nature, with small H ii regions forming and recombining. Only later, once halos of mass ∼109 and above begin to dominate the ionizing emissivity, does reionization proceed smoothly in the usual manner deduced from previous studies. This occurs at z ≈ 10 in our simulation.

Cosmic microwave background constraints in light of priors over reionization histories

Non-parametric reconstruction or marginalization over the history of reionization using cosmic microwave background data necessarily assumes a prior over possible histories. We show that different but reasonable choices of priors can shift current and future constraints on the reionization optical depth, τ, or correlated parameters such as the neutrino mass sum, Σmν, at the level of 0.3–0.4 σ, meaning that this analysis is somewhat prior dependent. We point out some prior-related problems with the commonly used principal component reconstruction, concluding that the significance of some recent hints of early reionization in Planck 2015 data has been overestimated. We also present the first non-parametric reconstruction applied to newer Planck intermediate (2016) data and find that the hints of early reionization disappear entirely in this more precise dataset. These results limit possible explanations of the EDGES 21cm signal which would have also significantly reionized the universe at z &gt; 15. Our findings about the dependence on priors motivate the pursuit of improved data or searches for physical reionization models which can reduce the prior volume. The discussion here of priors is of general applicability to other non-parametric reconstructions, for example of the primordial power spectrum, of the recombination history, or of the expansion rate.

Was there an early reionization component in our universe?

A deep understanding of the epoch of reionization is still missing in our knowledge of the universe. While future probes will allow us to test the precise evolution of the free electron fraction from redshifts between z≃ 6 and 0z≃ 2, at present one could ask what kind of reionization processes are allowed by present cosmic microwave background temperature and polarization measurements. An early contribution to reionization could imply a departure from the standard picture where star formation determines the reionization onset. By considering a broad class of possible reionization parameterizations, we find that current data do not require an early reionization component in our universe and that only one marginal class of models, based on a particular realization of reionization, may point to that. In addition, the frequentist Akaike information criterion (AIC) provides strong evidence against alternative reionization histories, favoring the most simple reionization scenario, which describes reionization by means of only one (constant) reionization optical depth τ.

Constraints from microlensing experiments on clustered primordial black holes

Abstract It has recently been proposed that massive primordial black holes (PBH) could constitute all of the dark matter, providing a novel scenario of structure formation, with early reionization and a rapid growth of the massive black holes at the center of galaxies and dark matter halos. The scenario arises from broad peaks in the primordial power spectrum that give both a spatially clustered and an extended mass distribution of PBH. The constraints from the observed microlensing events on the extended mass function have already been addressed. Here we study the impact of spatial clustering on the microlensing constraints. We find that the bounds can be relaxed significantly for relatively broad mass distributions if the number of primordial black holes within each cluster is typically above one hundred. On the other hand, even if they arise from individual black holes within the cluster, the bounds from CMB anisotropies are less stringent due to the enhanced black hole velocity in such dense clusters. This way, the window between a few and ten solar masses has opened up for PBH to comprise the totality of the dark matter.

Effect of the early reionization on the cosmic microwave background and cosmological parameter estimates

The early reionization (ERE) is supposed to be a physical process which happens after recombination, but before the instantaneous reionization caused by the first generation of stars. We investigate the effect of the ERE on the temperature and polarization power spectra of cosmic microwave background (CMB), and adopt principal components analysis (PCA) to model-independently reconstruct the ionization history during the ERE. In addition, we also discuss how the ERE affects the cosmological parameter estimates, and find that the ERE does not impose any significant influences on the tensor-to-scalar ratio r and the neutrino mass at the sensitivities of current experiments. The better CMB polarization data can be used to give a tighter constraint on the ERE and might be important for more precisely constraining cosmological parameters in the future.

Reionization during the dark ages from a cosmic axion background

Recently it has been pointed out that a cosmic background of relativistic axion-like particles (ALPs) would be produced by the primordial decays of heavy fields in the post-inflation epoch, contributing to the extra-radiation content in the Universe today. Primordial magnetic fields would trigger conversions of these ALPs into sub-MeV photons during the dark ages. This photon flux would produce an early reionization of the Universe, leaving a significant imprint on the total optical depth to recombination τ. Using the current measurement of τ and the limit on the extra-radiation content Δ Neff by the Planck experiment we put a strong bound on the ALP-photon conversions. Namely we obtain upper limits on the product of the photon-ALP coupling constant gaγ times the magnetic field strength B down to gaγ B ≳ 6 × 10−18 GeV−1 nG for ultralight ALPs.

Formation of massive Population III galaxies through photoionization feedback: a possible explanation for CR 7

Abstract We explore the formation of massive high-redshift Population III (Pop III) galaxies through photoionization feedback. We consider dark matter haloes formed from progenitors that have undergone no star formation as a result of early reionization and photoevaporation caused by a nearby galaxy. Once such a halo reaches ≈109 M⊙, corresponding to the Jeans mass of the photoheated intergalactic medium at z ≈ 7, pristine gas is able to collapse into the halo, potentially producing a massive Pop III starburst. We suggest that this scenario may explain the recent observation of strong He ii 1640 Å line emission in CR 7, which is consistent with ∼107 M⊙ of young Pop III stars. Such a large mass of Pop III stars is unlikely without the photoionization feedback scenario, because star formation is expected to inject metals into haloes above the atomic cooling threshold (∼108 M⊙ at z ≈ 7). We use merger trees to analytically estimate the abundance of observable Pop III galaxies formed through this channel, and find a number density of ≈10−7 Mpc−3 at z = 6.6 (the redshift of CR 7). This is approximately a factor of 10 lower than the density of Ly α emitters as bright as CR 7.

Lyman–Werner UV escape fractions from primordial haloes

Population III stars can regulate star formation in the primordial Universe in several ways. They can ionize nearby halos, and even if their ionizing photons are trapped by their own halos, their Lyman-Werner (LW) photons can still escape and destroy H$_2$ in other halos, preventing them from cooling and forming stars. LW escape fractions are thus a key parameter in cosmological simulations of early reionization and star formation but have not yet been parametrized for realistic halos by halo or stellar mass. To do so, we perform radiation hydrodynamical simulations of LW UV escape from 9--120 M$_{\odot}$ Pop III stars in $10^5$ to $10^7$ M$_{\odot}$ halos with ZEUS-MP. We find that photons in the LW lines (i.e. those responsible for destroying H$_{2}$ in nearby systems) have escape fractions ranging from 0% to 85%. No LW photons escape the most massive halo in our sample, even from the most massive star. Escape fractions for photons elsewhere in the 11.18--13.6~eV energy range, which can be redshifted into the LW lines at cosmological distances, are generally much higher, being above 60% for all but the least massive stars in the most massive halos. We find that shielding of H$_2$ by neutral hydrogen, which has been neglected in most studies to date, produces escape fractions that are up to a factor of three smaller than those predicted by H$_2$ self-shielding alone.

ON MODELING AND MEASURING THE TEMPERATURE OF THEz∼ 5 INTERGALACTIC MEDIUM

The temperature of the low-density intergalactic medium (IGM) at high redshift is sensitive to the timing and nature of hydrogen and HeII reionization, and can be measured from Lyman-alpha forest absorption spectra. Since the memory of intergalactic gas to heating during reionization gradually fades, measurements as close as possible to reionization are desirable. In addition, measuring the IGM temperature at sufficiently high redshifts should help to isolate the effects of hydrogen reionization since HeII reionization starts later, at lower redshift. Motivated by this, we model the IGM temperature at z>5 using semi-numeric models of patchy reionization. We construct mock Lyman-alpha forest spectra from these models and consider their observable implications. We find that the small-scale structure in the Lyman-alpha forest is sensitive to the temperature of the IGM even at redshifts where the average absorption in the forest is as high as 90%. We forecast the accuracy at which the z~5 IGM temperature can be measured using existing samples of high resolution quasar spectra, and find that interesting constraints are possible. For example, an early reionization model in which reionization ends at z~10 should be distinguishable -- at high statistical significance -- from a lower redshift model where reionization completes at z~6. We discuss improvements to our modeling that may be required to robustly interpret future measurements.

Radiative feedback from high-mass X-ray binaries on the formation of the first galaxies and early reionization

Recent work suggests that the first generation of stars, the so-called Population III (Pop III), could have formed primarily in binaries or as members of small multiple systems. Here we investigate the impact of X-ray feedback from High-Mass X-ray Binaries (HMXBs) left behind in stellar binary systems after the primary forms a black hole (BH), accreting gas at a high rate from the companion, a process that is thought to be favored at the low metallicities characteristic of high-redshift gas. Thanks to their large mean free path, X-rays are capable of preionizing and preheating the gas in the intergalactic medium (IGM) and in haloes long before the reionization of the Universe is complete, and thus could have strongly affected the formation of subsequent generations of stars as well as reionization. We have carried out zoomed hydrodynamical cosmological simulations of minihaloes, accounting for the formation of Pop III stars and their collapse into BHs and HMXBs, and the associated radiation-hydrodynamic feedback from UV and X-ray photons. We find no strong net feedback from HMXBs on the simulated star formation history. On the other hand, the preheating of the IGM by HMXBs leads to a strong suppression of small-scale structures and significantly lowers the recombination rate in the IGM, thus yielding a net positive feedback on reionization. We further show that X-ray feedback from HMXBs can augment the ionizing feedback from the Pop III progenitor stars to suppress gas accretion onto the first BHs, limiting their growth into supermassive BHs. Finally, we show that X-ray ionization by HMXBs leaves distinct signatures in the properties of the high-redshift hydrogen that may be probed in upcoming observations of the redshifted 21cm spin-flip line.

ON THE POSSIBLE ORIGIN OF THE LARGE SCALE COSMIC MAGNETIC FIELD

The possibility that the large scale cosmic magnetic field is directly generated at microgauss, equipartition levels during the reionization epoch by collisionless shocks that are forced to satisfy a downstream shear flow boundary condition is investigated through the development of two models—the accretion of an ionized plasma onto a weakly ionized cool galactic disk and onto a cool filament of the cosmic web. The dynamical structure and the physical parameters of the models are synthesized from recent cosmological simulations of the early reionization era after the formation of the first stars. The collisionless shock stands upstream of the disk and filament, and its dissipation is determined by ion inertial length Weibel turbulence. The downstream shear boundary condition is determined by the rotational neutral gas flow in the disk and the inward accretion flow along the filament. The shocked plasma is accelerated to the downstream shear flow velocity by the Weibel turbulence, and the relative shearing motion between the electrons and ions produces a strong, ion inertial scale current sheet that generates an equipartition strength, large scale downstream magnetic field, ~10–6 G for the disk and ~6 × 10–8 G for the filament. By assumption, hydrodynamic turbulence transports the shear-shock generated magnetic flux throughout the disk and filament volume.

Searching for the cosmic dawn

The Hubble Space Telescope, teaming up with a 'cosmic lens', has revealed a highly magnified galaxy thought to date back to 500 million years after the Big Bang. The find provides a glimpse of the first stages of galaxy formation. See Letter p.406 Young galaxies at a cosmic age of less than 500 million years remain largely unexplored because they are at or beyond the sensitivity limits of current large telescopes. This paper reports the use of strong gravitational lensing from a massive cluster of galaxies to observe a galaxy from the early Universe, at a redshift of z ≈ 9.6, equivalent to a cosmic age of approximately 490 million years. The authors suggest that because faint galaxies seem to be abundant at such a young cosmic age they are probably the dominant source for the early re-ionization of the intergalactic medium.

Constraints on superconducting cosmic strings from early reionization

Electromagnetic radiation from superconducting cosmic string loops can reionize neutral hydrogen in the universe at very early epochs, and affect the cosmic microwave background (CMB) temperature and polarization correlation functions at large angular scales. We constrain the string tension and current using WMAP7 data, and compare with earlier constraints that employed CMB spectral distortions. Over a wide range of string tensions, the current on the string has to be less than 10^7 GeV.