Reionization Of Cosmic Hydrogen Research Articles

A Strong He ii λ1640 Emitter with an Extremely Blue UV Spectral Slope at z = 8.16: Presence of Population III Stars?

Cosmic hydrogen reionization and cosmic production of the first metals are major phase transitions of the Universe occurring during the first billion years after the Big Bang; however, these are still underexplored observationally. Using the JWST/NIRSpec prism spectroscopy, we report the discovery of a sub-L * galaxy at z spec = 8.1623 ± 0.0007, dubbed RX J2129–z8He II, via the detection of a series of strong rest-frame UV/optical nebular emission lines and the clear Lyman break. RX J2129–z8He II shows a pronounced UV continuum with an extremely steep (i.e., blue) spectral slope of β=−2.53−0.07+0.06 , the steepest among all spectroscopically confirmed galaxies at z spec ≳ 7, in support of its very hard ionizing spectrum that could lead to a significant leakage of its ionizing flux. Therefore, RX J2129–z8He II is representative of the key galaxy population driving the cosmic reionization. More importantly, we detect a strong He II λ1640 emission line in its spectrum, one of the highest redshifts at which such a line is robustly detected. Its high rest-frame equivalent width (EW = 21 ± 4 Å) and extreme flux ratios with respect to UV metal and Balmer lines raise the possibility that part of RX J2129–z8He II’s stellar population could be Pop III (Pop III)-like. Through careful photoionization modeling, we show that the physically calibrated phenomenological models of the ionizing spectra of Pop III stars with strong mass loss can successfully reproduce the emission line flux ratios observed in RX J2129–z8He II. Assuming the Eddington limit, the total mass of the Pop III stars within this system is estimated to be 7.8 ± 1.4 × 105 M ⊙. To date, this galaxy presents the most compelling case in the early Universe where trace Pop III stars might coexist with metal-enriched populations.

The thesan project: connecting ionized bubble sizes to their local environments during the Epoch of Reionization

ABSTRACT An important characteristic of cosmic hydrogen reionization is the growth of ionized gas bubbles surrounding early luminous objects. Ionized bubble sizes are beginning to be probed using Lyman α emission from high-redshift galaxies, and will also be probed by upcoming 21 cm maps. We present results from a study of bubble sizes using the state-of-the-art thesan radiation-hydrodynamics simulation suite, which self-consistently models radiation transport and realistic galaxy formation. We employ the mean free path method and track the evolution of the effective ionized bubble size at each point (Reff) throughout the Epoch of Reionization. We show that there is a slow growth period for regions ionized early, but a rapid ‘flash ionization’ process for regions ionized later as they immediately enter a large, pre-existing bubble. We also find that bright sources are preferentially in larger bubbles, and find consistency with recent observational constraints at z ≳ 9, but tension with idealized Lyman α damping-wing models at z ≈ 7. We find that high-overdensity regions have larger characteristic bubble sizes, but the correlation decreases as reionization progresses, likely due to runaway formation of large percolated bubbles. Finally, we compare the redshift at which a region transitions from neutral to ionized (zreion) with the time it takes to reach a given bubble size and conclude that zreion is a reasonable local probe of small-scale bubble size statistics ($R_\text{eff} \lesssim 1\, \rm {cMpc}$). However, for larger bubbles, the correspondence between zreion and size statistics weakens due to the time delay between the onset of reionization and the expansion of large bubbles, particularly at high redshifts.

SILVERRUSH. XIII. A Catalog of 20,567 Lyα Emitters at z = 2−7 Identified in the Full-depth Data of the Subaru/HSC-SSP and CHORUS Surveys

We present 20,567 Lyα emitters (LAEs) at z = 2.2 − 7.3 that are photometrically identified by the SILVERRUSH program in a large survey area up to 25 deg2 with deep images of five broadband filters (grizy) and seven narrowband filters targeting Lyα lines at z = 2.2, 3.3, 4.9, 5.7, 6.6, 7.0, and 7.3 taken by the Hyper Suprime-Cam Subaru Strategic Program and the Cosmic HydrOgen Reionization Unveiled with Subaru survey. We select secure >5σ sources showing narrowband color excesses via Lyα break screening, taking into account the spatial inhomogeneity of limiting magnitudes. After removing spurious sources by careful masking and visual inspection of coadded and multiepoch images obtained over the 7 yr of the surveys, we construct LAE samples consisting of 6995, 4641, 726, 6124, 2058, 18, and 5 LAEs at z = 2.2, 3.3, 4.9, 5.7, 6.6, 7.0, and 7.3, respectively, although the z = 7.3 candidates are tentative. Our LAE catalogs contain 289 spectroscopically confirmed LAEs at the expected redshifts from previous work. We demonstrate that the number counts of our LAEs are consistent with previous studies with similar LAE selection criteria. The LAE catalogs will be made public on our project webpage with detailed descriptions of the content and ancillary information about the masks and limiting magnitudes.

Thermal and reionization history within a large-volume semi-analytic galaxy formation simulation

ABSTRACTWe predict the 21-cm global signal and power spectra during the Epoch of Reionization using the meraxes semi-analytic galaxy formation and reionization model, updated to include X-ray heating and thermal evolution of the intergalactic medium. Studying the formation and evolution of galaxies together with the reionization of cosmic hydrogen using semi-analytic models (such as M eraxes) requires N-body simulations within large volumes and high-mass resolutions. For this, we use a simulation of side-length 210 h−1 Mpc with 43203 particles resolving dark matter haloes to masses of $5\times 10^8 \rm{ }h^{-1}\, \mathrm{M_\odot }$. To reach the mass resolution of atomically cooled galaxies, thought to be the dominant population contributing to reionization, at z = 20 of $\sim 2\times 10^7 \text{ }h^{-1}\, \mathrm{M_\odot }$, we augment this simulation using the darkforest Monte Carlo merger tree algorithm (achieving an effective particle count of ∼1012). Using this augmented simulation, we explore the impact of mass resolution on the predicted reionization history as well as the impact of X-ray heating on the 21-cm global signal and the 21-cm power spectra. We also explore the cosmic variance of 21-cm statistics within 703 h−3 Mpc3 sub-volumes. We find that the midpoint of reionization varies by Δz ∼ 0.8 and that the cosmic variance on the power spectrum is underestimated by a factor of 2–4 at k ∼ 0.1–0.4 Mpc−1 due to the non-Gaussian nature of the 21-cm signal. To our knowledge, this work represents the first model of both reionization and galaxy formation which resolves low-mass atomically cooled galaxies while simultaneously sampling sufficiently large scales necessary for exploring the effects of X-rays in the early Universe.

The Sherwood–Relics simulations: overview and impact of patchy reionization and pressure smoothing on the intergalactic medium

ABSTRACT We present the Sherwood–Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for treating radiative transfer effects. This scheme follows the spatial variations in the ionizing radiation field, as well as the associated fluctuations in IGM temperature and pressure smoothing. It is computationally much cheaper than full radiation hydrodynamics simulations, and circumvents the difficult task of calibrating a galaxy formation model to observational constraints on cosmic reionization. Using this hybrid technique, we study the spatial fluctuations in IGM properties that are seeded by patchy cosmic reionization. We investigate the relevant physical processes and assess their impact on the z > 4 Lyman-α forest. Our main findings are: (i) consistent with previous studies patchy reionization causes large-scale temperature fluctuations that persist well after the end of reionization, (ii) these increase the Lyman-α forest flux power spectrum on large scales, and (iii) result in a spatially varying pressure smoothing that correlates well with the local reionization redshift. (iv) Structures evaporated or puffed up by photoheating cause notable features in the Lyman-α forest, such as flat-bottom or double-dip absorption profiles.

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) – XV. Constraining the cosmic reionization at 5.5 < z < 7

ABSTRACT Revealing the cosmic hydrogen reionization history is one of the main goals of the modern cosmology. z > 5 quasars (QSOs) have been used as back-lights to investigate the evolution of the intervening intergalactic medium (IGM) during the cosmic reionization since their first discovery. However, due to the small population of luminous QSOs (∼130 QSOs known to date), a tight constraint on the reionization history has not yet been placed. In this work, we aim to tighten the constraint using the 93 QSOs (5.5 < z < 7.1) recently discovered in the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQS) project. This is the largest QSO sample used to constrain the epoch of reionization. We measure the mean IGM Ly α transmission and the QSO near-zone size using the UV spectra of these QSOs. The mean IGM Ly α transmission rises above zero at z ≲ 6, indicating the end of the reionization. The near-zone sizes of the SHELLQs QSOs are consistent with sizes spanned by QSOs of lifetime tq ∼ 1–100 Myr in simulations. Due to the scatter created by the low signal-to-noise spectra and large Ly α redshift uncertainty, we cannot conclude whether the redshift evolution of the near-zone size is affected by the reionization effect.

The Contribution of AGN Accretion Disks to Hydrogen Reionization

Abstract We examine the contribution of high-redshift (z > 6) active galactic nuclei (AGNs) to cosmic hydrogen reionization, by tracing the growth and ionizing output of the first generation of supermassive black holes (SMBHs). Our calculations are anchored to the observed population of z ≃ 6 quasars, and trace back the evolving spectral energy distributions (SEDs) of the accretion flows that power these early AGNs and consider a variety of growth histories, including super-Eddington accretion. Compared to a fixed-shape SED, the evolving thin disks produce ionizing radiation that is higher by up to ∼80%. Across a variety of SMBH growth scenarios, the contribution of AGNs to reionization is limited to late epochs (z < 7), and remains subdominant compared to star-forming galaxies. This conclusion holds irrespective of the (still unknown) space density of low-luminosity z = 6 AGNs, and for growth scenarios that allow super-Eddington accretion. The contribution of AGNs to reionization can extend to earlier epochs (z ≳ 8) in scenarios with relatively slow SMBH mass growth, i.e., for low accretion rates and/or high spins. We finally demonstrate that our framework can reproduce the observed quasar proximity-zone sizes, and that compact proximity zones around z = 6 quasars can be explained by the late onset of super-Eddington accretion.

SILVERRUSH X: Machine Learning-aided Selection of 9318 LAEs at z = 2.2, 3.3, 4.9, 5.7, 6.6, and 7.0 from the HSC SSP and CHORUS Survey Data

Abstract We present a new catalog of 9318 Lyα emitter (LAE) candidates at z = 2.2, 3.3, 4.9, 5.7, 6.6, and 7.0 that are photometrically selected by the SILVERRUSH program with a machine learning technique from large area (up to 25.0 deg2) imaging data with six narrowband filters taken by the Subaru Strategic Program with Hyper Suprime-Cam and a Subaru intensive program, Cosmic HydrOgen Reionization Unveiled with Subaru. We construct a convolutional neural network that distinguishes between real LAEs and contaminants with a completeness of 94% and a contamination rate of 1%, enabling us to efficiently remove contaminants from the photometrically selected LAE candidates. We confirm that our LAE catalogs include 177 LAEs that have been spectroscopically identified in our SILVERRUSH programs and previous studies, ensuring the validity of our machine learning selection. In addition, we find that the object-matching rates between our LAE catalogs and our previous results are ≃80%–100% at bright NB magnitudes of ≲24 mag. We also confirm that the surface number densities of our LAE candidates are consistent with previous results. Our LAE catalogs will be made public on our project webpage.

SCORCH. III. Analytical Models of Reionization with Varying Clumping Factors

Abstract In the Simulations and Constructions of the Reionization of Cosmic Hydrogen (SCORCH) project, we compare analytical models of the hydrogen ionization fraction with radiation–hydrodynamic simulations. We derive analytical models of the mass-weighted hydrogen ionization fraction from the local ionization balance equations as a more accurate alternative to the widely adopted model based on the volume filling factor. In particular, our model has a recombination term quadratic in the ionization fraction, which is consistent with the two-body interaction nature of recombination. Then, we use the radiation–hydrodynamic simulations to study the clumping factors needed to solve the analytical equations and provide accurate fitting functions. We find that the ionized hydrogen clumping factors from our radiative transfer simulations are significantly different than those from other simulations that use a uniform photoionization background. In addition to redshift dependence, we also see the dependence of the ionized hydrogen clumping factor on ionization fraction, and we incorporate this into our fits. We calculate the reionization histories using our analytical models and clumping factors and compare with widely adopted models, and all of our models achieve <7% difference from simulation results, while the other models have >20% deviations. The Thomson optical depths from reionization calculated from our analytical models result in <5% deviation from simulations, while the previous analytical models have a >20% difference and could result in biased conclusions of the intergalactic medium reionization.

CHORUS. I. Cosmic HydrOgen Reionization Unveiled with Subaru: Overview

Abstract To determine the dominant sources for cosmic reionization, the evolution history of the global ionizing fraction, and the topology of the ionized regions, we have conducted a deep imaging survey using four narrow-band (NB) and one intermediate-band (IB) filters on the Subaru/Hyper Suprime-Cam (HSC), called Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS). The central wavelengths and full-widths-at-half-maximum of the CHORUS filters are, respectively, 386.2 nm and 5.5 nm for NB387, 526.0 nm and 7.9 nm for NB527, 717.1 nm and 11.1 nm for NB718, 946.2 nm and 33.0 nm for IB945, and 971.2 nm and 11.2 nm for NB973. This combination, including NB921 (921.5 nm and 13.5 nm) from the Subaru Strategic Program with HSC (HSC SSP), is carefully designed, as if they were playing a chorus, to observe multiple spectral features simultaneously, such as Lyman continuum, Lyα, C iv, and He ii for $z$ = 2–7. The observing field is the same as that of the deepest footprint of the HSC SSP in the COSMOS field and its effective area is about 1.6 deg2. We present an overview of the CHORUS project, which includes descriptions of the filter design philosophy, observations and data reduction, multiband photometric catalogs, assessments of the imaging quality, measurements of the number counts, and example use cases for the data. All the imaging data, photometric catalogs, masked pixel images, data of limiting magnitudes and point spread functions, results of completeness simulations, and source number counts are publicly available through the HSC SSP database.

A 16 deg2 survey of emission-line galaxies at z < 1.6 from HSC-SSP PDR2 and CHORUS

Abstract We have conducted a comprehensive survey of emission-line galaxies at z ≲ 1.6 based on narrow-band (NB) imaging data taken with Hyper Suprime-Cam (HSC) on the Subaru telescope. In this paper, we update the catalogs of Hα, [O iii], and [O ii] emission-line galaxies using the data from the second Public Data Release (PDR2) of the Subaru Strategic Program (SSP) of the HSC and Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) survey along with the spectroscopic redshifts for 2019 emission-line galaxies selected with the PDR1 data. The wider effective coverage of NB816 and NB921, 16.3 deg2 and 16.9 deg2, respectively, are available in the Deep and UltraDeep layers of HSC-SSP from the PDR2. The CHORUS survey provides us with data with additional three NBs (NB527, NB718, and NB973) in the COSMOS field in the UltraDeep layer (1.37 deg2). The five NB datasets allow us to investigate the star-forming galaxies presenting emission-lines at 14 specific redshifts ranging from z ∼ 1.6 down to z ∼ 0.05. We revisit the distribution of large-scale structures and luminosity functions (LFs) for the emission-line galaxies with the large samples of 75377 emission-line galaxies selected. The redshift revolution of LFs shows that the star formation rate densities (SFRDs) decreases monotonically from z ∼ 1.6, which is consistent with the cosmic SFRD known to-date. Our samples of emission-line galaxies covering a sufficiently large survey volume are useful to investigate the evolution of star-forming galaxies since the cosmic noon in a wide range of environments including galaxy clusters, filaments, and voids.

Bayesian approach to constraining the properties of ionized bubbles during reionization

ABSTRACT A possible way to study the reionization of cosmic hydrogen is by observing the large ionized regions (bubbles) around bright individual sources, e.g. quasars, using the redshifted 21 cm signal. It has already been shown that matched filter-based methods are not only able to detect the weak 21 cm signal from these bubbles but also aid in constraining their properties. In this work, we extend the previous studies to develop a rigorous Bayesian framework to explore the possibility of constraining the parameters that characterize the bubbles. To check the accuracy with which we can recover the bubble parameters, we apply our method on mock observations appropriate for the upcoming SKA1-low. For a region of size ≳50 cMpc around a typical quasar at redshift 7, we find that ≈20 h of integration with SKA1-low will be able to constrain the size and location of the bubbles, as well as the difference in the neutral hydrogen fraction inside and outside the bubble, with $\lesssim 10$ per cent precision. The recovery of the parameters are more precise and the signal-to-noise ratio of the detected signal is higher when the bubble sizes are larger and their shapes are close to spherical. Our method can be useful in identifying regions in the observed field that contain large ionized regions and hence are interesting for following up with deeper integration times.

A bound on the 12C/13C ratio in near-pristine gas with ESPRESSO

ABSTRACT Using science verification observations obtained with ESPRESSO at the Very Large Telescope (VLT) in 4UT mode, we report the first bound on the carbon isotope ratio $\rm ^{12}C/^{13}C$ of a quiescent, near-pristine damped Ly α (DLA) system at z = 2.34. We infer a limit $\rm log_{10}\, ^{12}C/^{13}C \gt +0.37\, (2\sigma)$. We use the abundance pattern of this DLA, combined with a stochastic chemical enrichment model, to infer the properties of the enriching stars, finding the total gas mass of this system to be $\log _{10}(M_{\rm gas}/{\rm M_{\odot }})=6.3^{+1.4}_{-0.9}$ and the total stellar mass to be log10(M⋆/M⊙) = 4.8 ± 1.3. The current observations disfavour enrichment by metal-poor asymptotic giant branch (AGB) stars with masses $\rm \lt 2.4\, M_{\odot }$, limiting the epoch at which this DLA formed most of its enriching stars. Our modelling suggests that this DLA formed very few stars until $\gtrsim 1$ Gyr after the cosmic reionization of hydrogen and, despite its very low metallicity ($\sim 1/1000$ of solar), this DLA appears to have formed most of its stars in the past few hundred Myr. Combining the inferred star formation history with evidence that some of the most metal-poor DLAs display an elevated [C/O] ratio at redshift z ≲ 3, we suggest that very metal-poor DLAs may have been affected by reionization quenching. Finally, given the simplicity and quiescence of the absorption features associated with the DLA studied here, we use these ESPRESSO data to place a bound on the possible variability of the fine-structure constant, Δα/α = (−1.2 ± 1.1) × 10−5.

Contribution from stars stripped in binaries to cosmic reionization of hydrogen and helium

Massive stars are often found in binary systems, and it has been argued that binary products boost the ionizing radiation of stellar populations. Accurate predictions for binary products are needed to understand and quantify their contribution to cosmic reionization. We investigate the contribution of stars stripped in binaries because (1) they are, arguably, the best-understood products of binary evolution, (2) we recently produced the first radiative transfer calculations for the atmospheres of these stripped stars that predict their ionizing spectra, and (3) they are very promising sources because they boost the ionizing emission of stellar populations at late times. This allows stellar feedback to clear the surroundings such that a higher fraction of their photons can escape and ionize the intergalactic medium. Combining our detailed predictions for the ionizing spectra with a simple cosmic reionization model, we estimate that stripped stars contributed tens of percent of the photons that caused cosmic reionization of hydrogen, depending on the assumed escape fractions. More importantly, stripped stars harden the ionizing emission. We estimate that the spectral index for the ionizing part of the spectrum can increase to −1 compared to ≲ − 2 for single stars. At high redshift, stripped stars and massive single stars combined dominate the He II-ionizing emission, but we expect that active galactic nuclei drive cosmic helium reionization. Further observational consequences we expect are (1) high ionization states for the intergalactic gas surrounding stellar systems, such as C IV and Si IV, and (2) additional heating of the intergalactic medium of up to a few thousand Kelvin. Quantifying these warrants the inclusion of accurate models for stripped stars and other binary products in full cosmological simulations.

CHORUS. II. Subaru/HSC Determination of the Lyα Luminosity Function at z = 7.0: Constraints on Cosmic Reionization Model Parameter

We present the Ly$\alpha$ luminosity function (LF) derived from 34 Ly$\alpha$ emitters (LAEs) at $z=7.0$ on the sky of $3.1$ deg$^2$, the largest sample compared to those in the literature obtained at a redshift $z\gtrsim7$. The LAE sample is made by deep large-area Subaru narrowband observations conducted by the Cosmic HydrOgen Reionization Unveiled with Subaru (CHORUS) project. The $z=7.0$ Ly$\alpha$ LF of our project is consistent with those of the previous DECam and Subaru studies at the bright and faint ends, respectively, while our $z=7.0$ Ly$\alpha$ LF has uncertainties significantly smaller than those of the previous study results. Exploiting the small errors of our measurements, we investigate the shape of the faint to bright-end Ly$\alpha$ LF. We find that the $z=7.0$ Ly$\alpha$ LF shape can be explained by the steep slope of $\alpha \simeq -2.5$ suggested at $z=6.6$, and that there is no clear signature of a bright-end excess at $z\simeq 7$ claimed by the previous work, which was thought to be made by the ionized bubbles around bright LAEs whose Ly$\alpha$ photons could easily escape from the partly neutral IGM at $z \simeq 7$. We estimate the Ly$\alpha$ luminosity densities (LDs) with Ly$\alpha$ LFs at $z\simeq 6-8$ given by our and the previous studies, and compare the evolution of the UV-continuum LD estimated with dropouts. The Ly$\alpha$ LD monotonically decreases from $z\sim 6$ to $8$, and evolves stronger than the UV-continuum LD, indicative of the Ly$\alpha$ damping wing absorption of the IGM towards the heart of the reionization epoch.

Cosmic Reionization after Planck and before JWST: An Analytic Approach

Abstract The reionization of cosmic hydrogen marks a critical juncture in the history of structure formation. Here we present a new formulation of the standard reionization equation for the evolution of the volume-averaged H ii fraction that is more consistent with the accepted conceptual model of inhomogeneous intergalactic absorption. The revised equation explicitly accounts for the presence of the optically thick “Lyman-limit systems” that are known to determine the mean-free path of ionizing radiation after overlap. Integration of this equation provides a better characterization of the timing of reionization by smoothly linking the pre-overlap with the post-overlap phases of such a process. We confirm the validity of the quasi-instantaneous approximation as a predictor of reionization completion/maintenance and discuss new insights on the sources of cosmic reionization using the improved formalism. A constant emission rate into the intergalactic medium (IGM) of three Lyman continuum (LyC) photons per atom per gigayear leads to a reionization history that is consistent with a number of observational constraints on the ionization state of the z = 5–9 universe. While star-forming galaxies can dominate the reionization process if the luminosity-weighted fraction of LyC photons that escape into the IGM, , exceeds 15% (for a faint magnitude cut-off of the galaxy UV luminosity function of and a LyC photon yield per unit 1500 Å luminosity of ), simple models where the product of the two unknowns is not evolving with redshift fail to reproduce the changing neutrality of the IGM observed at these epochs.

No evidence for a significant AGN contribution to cosmic hydrogen reionization

We reinvestigate a claimed sample of 22 X-ray detected active galactic nuclei (AGN) at redshifts z > 4, which has reignited the debate as to whether young galaxies or AGN reionized the Universe. These sources lie within the GOODS-S/CANDELS field, and we examine both the robustness of the claimed X-ray detections (within the Chandra 4Ms imaging) and perform an independent analysis of the photometric redshifts of the optical/infrared counterparts. We confirm the reality of only 15 of the 22 reported X-ray detections, and moreover find that only 12 of the 22 optical/infrared counterpart galaxies actually lie robustly at z > 4. Combining these results we find convincing evidence for only 7 X-ray AGN at z > 4 in the GOODS-S field, of which only one lies at z > 5. We recalculate the evolving far-UV (1500 Angstrom) luminosity density produced by AGN at high redshift, and find that it declines rapidly from z = 4 to z = 6, in agreement with several other recent studies of the evolving AGN luminosity function. The associated rapid decline in inferred hydrogen-ionizing emissivity contributed by AGN falls an order-of-magnitude short of the level required to maintain hydrogen ionization at z ~ 6. We conclude that all available evidence continues to favour a scenario in which young galaxies reionized the Universe, with AGN making, at most, a very minor contribution to cosmic hydrogen reionization.

SCORCH. I. THE GALAXY–HALO CONNECTION IN THE FIRST BILLION YEARS

SCORCH (Simulations and Constructions of the Reionization of Cosmic Hydrogen) is a new project to study the Epoch of Reionization (EoR). In this first paper, we probe the connection between observed high-redshift galaxies and simulated dark matter halos to better understand the abundance and evolution of the primary source of ionizing radiation. High-resolution N-body simulations are run to quantify the abundance of dark matter halos as a function of mass $M$, accretion rate $\dot{M}$, and redshift $z$. A new fit for the halo mass function $dn/dM$ is $\approx 20\%$ more accurate at the high-mass end where bright galaxies are expected to reside. A novel approach is used to fit the halo accretion rate function $dn/d\dot{M}$ in terms of the halo mass function. Abundance matching against the observed galaxy luminosity function is used to estimate the luminosity-mass relation and the luminosity-accretion-rate relation. The inferred star formation efficiency is not monotonic with $M$ nor $\dot{M}$, but reaches a maximum value at a characteristic mass $\sim 2 \times 10^{11}\ M_\odot$ and a characteristic accretion rate $\sim 6 \times 10^2\ M_\odot/{\rm yr}$ at $z \approx 6$. We find a universal EoR luminosity-accretion-rate relation and construct a fiducial model for the galaxy luminosity function. The Schechter parameters evolve such that $\phi_\star$ decreases, $M_\star$ is more positive (fainter), and $\alpha$ is more negative (steeper) at higher redshifts. We forecast for the upcoming James Webb Space Telescope and show that with apparent magnitude limit $m_{\rm AB} \approx 31\ (32)$, it can observe $\gtrsim 11\ (24)$ unlensed galaxies per square degree per unit redshift at least down to $M_\star$ at $z \lesssim 13\ (14)$.

Cosmic reionization of hydrogen and helium: contribution from both mini-quasars and stars

Observations on the high-redshift galaxies at $z>6$ imply that their ionizing emissivity is unable to fully reionize the Universe at $z\sim 6$. Either a high escape fraction of ionizing photons from these galaxies or a large population of faint galaxies below the detection limit are required. However, these requirements are somewhat in tension with present observations. In this work, we explored the combined contribution of mini-quasars and stars to the reionization of cosmic hydrogen and helium. Our model is roughly consistent with: (1) the low escape fractions of ionizing photons from the observed galaxies, (2) the optical depth of Cosmic Microwave Background (CMB) measured by the WMAP-7, and (3) the redshift of the end of hydrogen and helium reionization at $z\approx 6$ and $z\approx 3$, respectively. Neither an extremely high escape fraction nor a large population of fainter galaxies is required in this scenario. In our most optimistic model, more than $\sim20\%$ of the cosmic helium is reionized by $z\sim6$, and the ionized fraction of cosmic helium rapidly climbs to more than $50\%$ by $z\sim5$. These results may imply that better measurements of helium reionization, especially at high redshifts, could be helpful in constraining the growth of intermediate-mass black holes (IMBHs) in the early Universe, which would shed some light on the puzzles concerning the formation of supermassive black holes (SMBHs).

Large-scale environment of z ∼ 5.7 C iv absorption systems –II. Spectroscopy of Lyman α emitters★

The flow of baryons to and from a galaxy, which is fundamental for galaxy formation and evolution, can be studied with galaxy-metal absorption system pairs. Our search for galaxies around CIV absorption systems at $z\sim5.7$ showed an excess of photometric Lyman-$\alpha$ emitter (LAE) candidates in the fields J1030+0524 and J1137+3549. Here we present spectroscopic follow-up of 33 LAEs in both fields. In the first field, three out of the five LAEs within 10$h^{{-}1}$ projected comoving Mpc from the CIV system are within $\pm500$ km s$^{{-}1}$ from the absorption at $z_{\text{CIV}}=5.7242\pm0.0001$. The closest candidate (LAE 103027+052419) is robustly confirmed at $212.8^{+14}_{-0.4}h^{-1}$ physical kpc from the CIV system. In the second field, the LAE sample is selected at a lower redshift ($\Delta z\sim0.04$) than the CIV absorption system as a result of the filter transmission and, thus, do not trace its environment. The observed properties of LAE 103027+052419 indicate that it is near the high-mass end of the LAE distribution, probably having a large HI column density and large-scale outflows. Therefore, our results suggest that the CIV system is likely produced by a star-forming galaxy which has been injecting metals into the intergalactic medium since $z>6$. Thus, the CIV system is either produced by LAE 103027+052419, implying that outflows can enrich larger volumes at $z>6$ than at $z\sim3.5$, or an undetected dwarf galaxy. In either case, CIV systems like this one trace the ionized intergalactic medium at the end of cosmic hydrogen reionization and may trace the sources of the ionizing flux density.