Search

Precision calibration poses challenges to experiments probing the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization (z~30-6). In both interferometric and global signal experiments, systematic calibration is the leading source of error. Though many aspects of calibration have been studied, the overlap between the two types of instruments has received less attention. We investigate the sky based calibration of total power measurements with a HERA dish and an EDGES style antenna to understand the role of auto-correlations in the calibration of an interferometer and the role of sky in calibrating a total power instrument. Using simulations we study various scenarios such as time variable gain, incomplete sky calibration model, and primary beam model. We find that temporal gain drifts, sky model incompleteness, and beam inaccuracies cause biases in the receiver gain amplitude and the receiver temperature estimates. In some cases, these biases mix spectral structure between beam and sky resulting in spectrally variable gain errors. Applying the calibration method to the HERA and EDGES data, we find good agreement with calibration via the more standard methods. Although instrumental gains are consistent with beam and sky errors similar in scale to those simulated, the receiver temperatures show significant deviations from expected values. While we show that it is possible to partially mitigate biases due to model inaccuracies by incorporating a time-dependent gain model in calibration, the resulting errors on calibration products are larger and more correlated. Completely addressing these biases will require more accurate sky and primary beam models., Comment: 16 pages, 13 figures, 1 table; accepted for publication in MNRAS main journal

In recent years, simple models of galaxy formation have been shown to provide reasonably good matches to available data on high-redshift luminosity functions. However, these prescriptions are primarily phenomenological, with only crude connections to the physics of galaxy evolution. Here we introduce a set of galaxy models that are based on a simple physical framework but incorporate more sophisticated models of feedback, star formation, and other processes. We apply these models to the high-redshift regime, showing that most of the generic predictions of the simplest models remain valid. In particular, the stellar mass--halo mass relation depends almost entirely on the physics of feedback (and is thus independent of the details of small-scale star formation) and the specific star formation rate is a simple multiple of the cosmological accretion rate. We also show that, in contrast, the galaxy's gas mass is sensitive to the physics of star formation, although the inclusion of feedback-driven star formation laws significantly changes the naive expectations. While these models are far from detailed enough to describe every aspect of galaxy formation, they inform our understanding of galaxy formation by illustrating several generic aspects of that process, and they provide a physically-grounded basis for extrapolating predictions to faint galaxies and high redshifts currently out of reach of observations. If observations show violations from these simple trends, they would indicate new physics occurring inside the earliest generations of galaxies., 20 pages, 13 figures, in press at MNRAS

Radio Frequency Interference (RFI) is one of the systematic challenges preventing 21cm interferometric instruments from detecting the Epoch of Reionization. To mitigate the effects of RFI on data analysis pipelines, numerous inpaint techniques have been developed to restore RFI corrupted data. We examine the qualitative and quantitative errors introduced into the visibilities and power spectrum due to inpainting. We perform our analysis on simulated data as well as real data from the Hydrogen Epoch of Reionization Array (HERA) Phase 1 upper limits. We also introduce a convolutional neural network that capable of inpainting RFI corrupted data in interferometric instruments. We train our network on simulated data and show that our network is capable at inpainting real data without requiring to be retrained. We find that techniques that incorporate high wavenumbers in delay space in their modeling are best suited for inpainting over narrowband RFI. We also show that with our fiducial parameters Discrete Prolate Spheroidal Sequences (DPSS) and CLEAN provide the best performance for intermittent ``narrowband'' RFI while Gaussian Progress Regression (GPR) and Least Squares Spectral Analysis (LSSA) provide the best performance for larger RFI gaps. However we caution that these qualitative conclusions are sensitive to the chosen hyperparameters of each inpainting technique. We find these results to be consistent in both simulated and real visibilities. We show that all inpainting techniques reliably reproduce foreground dominated modes in the power spectrum. Since the inpainting techniques should not be capable of reproducing noise realizations, we find that the largest errors occur in the noise dominated delay modes. We show that in the future, as the noise level of the data comes down, CLEAN and DPSS are most capable of reproducing the fine frequency structure in the visibilities of HERA data., Comment: 21 pages, 13 figures

We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for all metrics used, and present tailored visualizations that aid us in clearly identifying new and existing systematics. We implement these techniques using data from 105 antennas in the Hydrogen Epoch of Reionization Array (HERA) as a case study. Finally, we provide a detailed algorithm for implementing these metrics as flagging tools on real data sets., Comment: 31 pages, 17 figures

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz., 15 pages, 15 figures, 1 table, Accepted to MNRAS

FARSIDE (Farside Array for Radio Science Investigations of the Dark ages and Exoplanets) is a Probe-class concept to place a low radio frequency interferometric array on the farside of the Moon. A NASA-funded design study, focused on the instrument, a deployment rover, the lander and base station, delivered an architecture broadly consistent with the requirements for a Probe mission. This notional architecture consists of 128 dual polarization antennas deployed across a 10 km area by a rover, and tethered to a base station for central processing, power and data transmission to the Lunar Gateway. FARSIDE would provide the capability to image the entire sky each minute in 1400 channels spanning frequencies from 100 kHz to 40 MHz, extending down two orders of magnitude below bands accessible to ground-based radio astronomy. The lunar farside can simultaneously provide isolation from terrestrial radio frequency interference, auroral kilometric radiation, and plasma noise from the solar wind. This would enable near-continuous monitoring of the nearest stellar systems in the search for the radio signatures of coronal mass ejections and energetic particle events, and would also detect the magnetospheres for the nearest candidate habitable exoplanets. Simultaneously, FARSIDE would be used to characterize similar activity in our own solar system, from the Sun to the outer planets, including the hypothetical Planet Nine. Through precision calibration via an orbiting beacon, and exquisite foreground characterization, FARSIDE would also measure the Dark Ages global 21-cm signal at redshifts z=50-100. The unique observational window offered by FARSIDE would enable an abundance of additional science ranging from sounding of the lunar subsurface to characterization of the interstellar medium in the solar system neighborhood.

Precise measurements of the 21 cm power spectrum are crucial for understanding the physical processes of hydrogen reionization. Currently, this probe is being pursued by low-frequency radio interferometer arrays. As these experiments come closer to making a first detection of the signal, error estimation will play an increasingly important role in setting robust measurements. Using the delay power spectrum approach, we have produced a critical examination of different ways that one can estimate error bars on the power spectrum. We do this through a synthesis of analytic work, simulations of toy models, and tests on small amounts of real data. We find that, although computed independently, the different error bar methodologies are in good agreement with each other in the noise-dominated regime of the power spectrum. For our preferred methodology, the predicted probability distribution function is consistent with the empirical noise power distributions from both simulated and real data. This diagnosis is mainly in support of the forthcoming HERA upper limit, and also is expected to be more generally applicable., Comment: 35 Pages, 9 Figures, 4 Tables. Replaced with accepted ApJ version; some clarifying text added in response to referee comments with no changes to results

As current- and next-generation astronomical instruments come online, they will generate an unprecedented deluge of data. Analyzing these data in real time presents unique conceptual and computational challenges, and their long-term storage and archiving is scientifically essential for generating reliable, reproducible results. We present here the real-time processing (RTP) system for the Hydrogen Epoch of Reionization Array (HERA), a radio interferometer endeavoring to provide the first detection of the highly redshifted 21 cm signal from Cosmic Dawn and the Epoch of Reionization by an interferometer. The RTP system consists of analysis routines run on raw data shortly after they are acquired, such as calibration and detection of radio-frequency interference (RFI) events. RTP works closely with the Librarian, the HERA data storage and transfer manager which automatically ingests data and transfers copies to other clusters for post-processing analysis. Both the RTP system and the Librarian are public and open source software, which allows for them to be modified for use in other scientific collaborations. When fully constructed, HERA is projected to generate over 50 terabytes (TB) of data each night, and the RTP system enables the successful scientific analysis of these data., Comment: 15 pages, 3 figures, published in Astronomy and Computing

The cosmic near-infrared background (NIRB) offers a powerful integral probe of radiative processes at different cosmic epochs, including the pre-reionization era when metal-free, Population III (Pop III) stars first formed. While the radiation from metal-enriched, Population II (Pop II) stars likely dominates the contribution to the observed NIRB from the reionization era, Pop III stars -- if formed efficiently -- might leave characteristic imprints on the NIRB thanks to their strong Ly$��$ emission. Using a physically-motivated model of first star formation, we provide an analysis of the NIRB mean spectrum and anisotropy contributed by stellar populations at $z>5$. We find that in circumstances where massive Pop III stars persistently form in molecular cooling haloes at a rate of a few times $10^{-3}\,M_\odot \ \mathrm{yr}^{-1}$, before being suppressed towards the epoch of reionization (EoR) by the accumulated Lyman-Werner background, a unique spectral signature shows up redward of $1\,��$m in the observed NIRB spectrum sourced by galaxies at $z>5$. While the detailed shape and amplitude of the spectral signature depend on various factors including the star formation histories, IMF, LyC escape fraction and so forth, the most interesting scenarios with efficient Pop III star formation are within the reach of forthcoming facilities such as the Spectro-Photometer for the History of the Universe, Epoch of Reionization and Ices Explorer (SPHEREx). As a result, new constraints on the abundance and formation history of Pop III stars at high redshifts will be available through precise measurements of the NIRB in the next few years., 19 pages, 10 figures, accepted for publication in MNRAS

The reionization epoch concludes when ionizing photons reach every corner of the Universe. Reionization has generally been assumed to be limited primarily by the rate at which galaxies produce ionizing photons, but the recent measurement of a surprisingly short ionizing photon mean free path of $0.75^{+0.65}_{-0.45}$ proper Mpc at $z = 6$ by Becker et al. (2021) suggests that absorption by residual neutral hydrogen in the otherwise ionized intergalactic medium may play a much larger role than previously expected. Here we show that consistency between this short mean free path and the coeval dark pixel fraction in the Ly$\alpha$ forest requires a cumulative output of $6.1^{+11}_{-2.4}$ ionizing photons per baryon by reionization's end. This represents a dramatic increase in the ionizing photon budget over previous estimates, greatly exacerbating the tension with measurements of the ionizing output from galaxies at later times. Translating this constraint into the instantaneous ionizing production from galaxies in our model, we find $\log_{10}f_{\rm esc}\xi_{\rm ion}/\text{(erg/Hz)}^{-1} =25.02_{-0.21}^{+0.45}$ at $z\sim6$. Even with optimistic assumptions about the ionizing production efficiency of early stellar populations, and assuming the galaxy luminosity function extends to extremely faint sources ($M_{\text{UV}}\leq-11$), complete reionization requires the escape fraction of ionizing photons to exceed $20\%$ across the galaxy population. This is far larger than observed in any galaxy population at lower redshifts, requiring rapid evolution in galaxy properties after the first billion years of cosmic time. This tension cannot be completely relieved within existing observational constraints on the hydrogen neutral fraction and mean free path., Comment: 8 pages, 5 figures, submitted to ApJL

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from very bright but spectrally-smooth astrophysical foregrounds. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as "redundant-baseline calibration" resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions--both in data and in simulations--and present strategies for mitigating that structure., Comment: 24 Pages, 19 Figures. Updated to match the accepted MNRAS version

Author(s): Mebane, Richard H; Mirocha, Jordan; Furlanetto, Steven R | Abstract: ABSTRACT We investigate the effects of Population III (Pop III) stars and their remnants on the cosmological 21-cm global signal. By combining a semi-analytic model of Pop III star formation with a global 21-cm simulation code, we investigate how X-ray and radio emission from accreting Pop III black holes may affect both the timing and depth of the 21-cm absorption feature that follows the initial onset of star formation during the Cosmic Dawn. We compare our results to the findings of the EDGES experiment, which has reported the first detection of a cosmic 21-cm signal. In general, we find that our fiducial Pop III models, which have peak star formation rate densities of ∼10−4 M⊙ yr−1 Mpc−3 between z ∼ 10 and z ∼ 15, are able to match the timing of the EDGES signal quite well, in contrast to models that ignore Pop III stars. To match the unexpectedly large depth of the EDGES signal without recourse to exotic physics, we vary the parameters of emission from accreting black holes (formed as Pop III remnants) including the intrinsic strength of X-ray and radio emission as well as the local column density of neutral gas. We find that models with strong radio emission and relatively weak X-ray emission can self-consistently match the EDGES signal, though this solution requires fine-tuning. We are only able to produce signals with sharp features similar to the EDGES signal if the Pop III IMF is peaked narrowly around $140 \, \mathrm{M}_\odot$.

The experimental efforts to detect the redshifted 21 cm signal from the Epoch of Reionization (EoR) are limited predominantly by the chromatic instrumental systematic effect. The delay spectrum methodology for 21 cm power spectrum measurements brought new attention to the critical impact of an antenna's chromaticity on the viability of making this measurement. This methodology established a straightforward relationship between time-domain response of an instrument and the power spectrum modes accessible to a 21 cm EoR experiment. We examine the performance of a prototype of the Hydrogen Epoch of Reionization Array (HERA) array element that is currently observing in Karoo desert, South Africa. We present a mathematical framework to derive the beam integrated frequency response of a HERA prototype element in reception from the return loss measurements between 100-200 MHz and determined the extent of additional foreground contamination in the delay space. The measurement reveals excess spectral structures in comparison to the simulation studies of the HERA element. Combined with the HERA data analysis pipeline that incorporates inverse covariance weighting in optimal quadratic estimation of power spectrum, we find that in spite of its departure from the simulated response, HERA prototype element satisfies the necessary criteria posed by the foreground attenuation limits and potentially can measure the power spectrum at spatial modes as low as $k_{\parallel} > 0.1h$~Mpc$^{-1}$. The work highlights a straightforward method for directly measuring an instrument response and assessing its impact on 21 cm EoR power spectrum measurements for future experiments that will use reflector-type antenna., 13 pages, 10 figures, Submitted to APJ

In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark matter simulations, subhalo abundance matching and empirical models, and includes galaxy positions, morphologies, and spectral energy distributions (SEDs). The resulting synthetic catalog extends to redshifts $z \sim 12 $, and galaxy masses $\log_{10}(M/M_{\odot}) = 5 $ covering an area of $1\, {\rm deg}^2$ on the sky. We use DREaM to explore the science returns of a $1\, {\rm deg}^2$ \emph{Roman} UDF, and to provide a resource for optimizing ultra-deep survey designs. We find that a \emph{Roman} UDF to $\sim 30\,m_{\rm AB}$ will potentially detect more than $10^6$ $M_{\rm UV}7$, offering an unparalleled dataset for constraining galaxy properties during the EoR. Our synthetic catalogs and simulated images are made publicly available to provide the community with a tool to prepare for upcoming data., Comment: 39 pages, 27 figures. Submitted to ApJ. Comments welcome

The observed large-scale scatter in Lyman $\alpha$ opacity of the intergalactic medium at $z, Comment: Accepted to ApJ. 17 pages, 13 figures

Cosmic variance is the intrinsic scatter in the number density of galaxies due to fluctuations in the large-scale dark matter density field. In this work, we present a simple analytic model of cosmic variance in the high redshift Universe ($z\sim5-15$). We assume that galaxies grow according to the evolution of the halo mass function, which we allow to vary with large-scale environment. Our model produces a reasonable match to the observed ultraviolet luminosity functions in this era by regulating star formation through stellar feedback and assuming that the UV luminosity function is dominated by recent star formation. We find that cosmic variance in the UVLF is dominated by the variance in the underlying dark matter halo population, and not by differences in halo accretion or the specifics of our stellar feedback model. We also find that cosmic variance dominates over Poisson noise for future high-$z$ surveys except for the brightest sources or at very high redshifts ($z \gtrsim 12$). We provide a linear approximation of cosmic variance for a variety of redshifts, magnitudes, and survey areas through the public Python package galcv. Finally, we introduce a new method for incorporating priors on cosmic variance into estimates of the galaxy luminosity function and demonstrate that it significantly improves constraints on that important observable., Comment: 16 pages, 15 figures, accepted by MNRAS

Radio Frequency Interference (RFI) is an ever-present limiting factor among radio telescopes even in the most remote observing locations. When looking to retain the maximum amount of sensitivity and reduce contamination for Epoch of Reionization studies, the identification and removal of RFI is especially important. In addition to improved RFI identification, we must also take into account computational efficiency of the RFI-Identification algorithm as radio interferometer arrays such as the Hydrogen Epoch of Reionization Array grow larger in number of receivers. To address this, we present a Deep Fully Convolutional Neural Network (DFCN) that is comprehensive in its use of interferometric data, where both amplitude and phase information are used jointly for identifying RFI. We train the network using simulated HERA visibilities containing mock RFI, yielding a known "ground truth" dataset for evaluating the accuracy of various RFI algorithms. Evaluation of the DFCN model is performed on observations from the 67 dish build-out, HERA-67, and achieves a data throughput of 1.6$\times 10^{5}$ HERA time-ordered 1024 channeled visibilities per hour per GPU. We determine that relative to an amplitude only network including visibility phase adds important adjacent time-frequency context which increases discrimination between RFI and Non-RFI. The inclusion of phase when predicting achieves a Recall of 0.81, Precision of 0.58, and $F_{2}$ score of 0.75 as applied to our HERA-67 observations., Comment: 11 pages, 7 figures

We present deep Keck/MOSFIRE $H$-band spectroscopic observations covering the [CIII],CIII]$\lambda\lambda1907,1909$ doublet for three $z\sim8$ galaxy candidates in the AEGIS field. Along with non-detections in two galaxies, we obtain one of the highest-redshift detections to-date of [CIII]$\lambda 1907$ for the galaxy AEGIS-33376, finding $z_{\rm spec}=7.945\pm0.001$. We measure a [CIII]$\lambda$1907 flux of $2.24\pm0.71\times10^{-18} \mbox{ erg}\mbox{ s}^{-1} \mbox{ cm}^{-2}$, corresponding to a rest-frame equivalent width of $20.3\pm6.5 \unicode{x212B}$ for the single line. Given the not very constraining upper limit for CIII]$\lambda 1909$ based on strong sky-line contamination, we assume a [CIII]$\lambda$1907/CIII]$\lambda 1909$ doublet ratio of 1.5 and infer a total [CIII],CIII]$\lambda\lambda1907,1909$ equivalent width of $33.7\pm 10.8 \unicode{x212B}$. We repeat the same reductions and analysis on multiple subsets of our raw data divided on the basis of time and observing conditions, verifying that the [CIII]$\lambda 1907$ emission is present for AEGIS-33376 throughout our observations. We also confirm that the significance of the [CIII]$\lambda 1907$ detection in different subsets of our data tracks that of brighter emission features detected on the same multi-slit mask. These multiple tests suggest that the observed emission line is real and associated with the $z\sim 8$ target. The strong observed [CIII],CIII]$\lambda\lambda1907,1909$ in AEGIS-33376 likely indicates ISM conditions of low metallicity, high ionization parameter, and a hard ionizing spectrum, although AGN contributions are possible. This single detection represents a sizable increase in the current sample [CIII],CIII]$\lambda\lambda1907,1909$ detections at $z>7$, while $\textit{JWST}$ will provide the first statistical samples of such measurements at these redshifts., Comment: 6 pages, 3 figures, 1 table; Accepted to ApJL

A key physical quantity during reionization is the size of HII regions. Previous studies found a characteristic bubble size which increases rapidly during reionization, with apparent agreement between simulations and analytic excursion set theory. Using four different methods, we critically examine this claim. In particular, we introduce the use of the watershed algorithm -- widely used for void finding in galaxy surveys -- which we show to be an unbiased method with the lowest dispersion and best performance on Monte-Carlo realizations of a known bubble size PDF. We find that a friends-of-friends algorithm declares most of the ionized volume to be occupied by a network of volume-filling regions connected by narrow tunnels. For methods tuned to detect the volume-filling regions, previous apparent agreement between simulations and theory is spurious, and due to a failure to correctly account for the window function of measurement schemes. The discrepancy is already obvious from visual inspection. Instead, HII regions in simulations are significantly larger (by factors of 10-1000 in volume) than analytic predictions. The size PDF is narrower, and evolves more slowly with time, than predicted. It becomes more sharply peaked as reionization progresses. These effects are likely caused by bubble mergers, which are inadequately modeled by analytic theory. Our results have important consequences for high-redshift 21cm observations, the mean free path of ionizing photons, and the visibility of Ly-alpha emitters, and point to a fundamental failure in our understanding of the characteristic scales of the reionization process., replaced to match published version in MNRAS, conclusions unchanged

The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz., 22 pages, 22 figures, accepted for publication in ApJ

The Experiment to Detect the Global Epoch of Reionization Signature (EDGES) recently reported a strong 21-cm absorption signal relative to the cosmic microwave background at $z \sim 18$. While its anomalous amplitude may indicate new physics, in this work we focus on the timing of the signal, as it alone provides an important constraint on galaxy formation models. Whereas rest-frame ultraviolet luminosity functions (UVLFs) over a broad range of redshifts are well fit by simple models in which galaxy star formation histories track the assembly of dark matter halos, we find that these same models, with reasonable assumptions about X-ray production in star-forming galaxies, cannot generate a narrow absorption trough at $z \sim 18$. If verified, the EDGES signal therefore requires the fundamental inputs of galaxy formation models to evolve rapidly at $z \gtrsim 10$. Unless extremely faint sources residing in halos below the atomic cooling threshold are responsible for the EDGES signal, star formation in $\sim 10^8$-$10^{10} \ M_{\odot}$ halos must be more efficient than expected, implying that the faint-end of the UVLF at $M_{\mathrm{UV}} \lesssim -12$ must steepen at the highest redshifts. This steepening provides a concrete test for future galaxy surveys with the James Webb Space Telescope and ongoing efforts in lensed fields, and is required regardless of whether the amplitude of the EDGES signal is due to new cooling channels or a strong radio background in the early Universe. However, the radio background solution requires that galaxies at $z > 15$ emit 1-2 GHz photons with an efficiency $\sim 10^3$ times greater than local star-forming galaxies, posing a challenge for models of low-frequency photon production in the early Universe., 13 pages, 5 figures, submitted to MNRAS

The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Applied to the Hydrogen Epoch of Reionization Array (HERA), this paper combines detailed electromagnetic and electrical simulations in order to model the chromatic effects of the instrument, and quantify its frequency and time responses. In particular, the effects of the analogue receiver, transmission cables, and mutual coupling are included. These simulations are able to accurately predict the intensity of the reflections occurring in the 150-m cable which links the antenna to the back-end. They also show that electromagnetic waves can propagate from one dish to another one through large sections of the array due to mutual coupling. The simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h\;\rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity., Comment: 13 pages, 17 figures - Submitted to MNRAS - 2nd revision

We present a semi-analytic model of star formation in the early universe, beginning with the first metal-free stars. By employing a completely feedback-limited star formation prescription, stars form at maximum efficiency until the self-consistently calculated feedback processes halt formation. We account for a number of feedback processes including a meta-galactic Lyman-Werner background, supernovae, photoionization, and chemical feedback. Halos are evolved combining mass accretion rates found through abundance matching with our feedback-limited star formation prescription, allowing for a variety of Population III (Pop III) initial mass functions (IMFs). We find that, for a number of models, massive Pop III star formation can continue on until at least $z \sim 20$ and potentially past $z \sim 6$ at rates of around $10^{-4}$ to $10^{-5}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$, assuming these stars form in isolation. At this point Lyman-Werner feedback pushes the minimum halo mass for star formation above the atomic cooling threshold, cutting off the formation of massive Pop III stars. We find that, in most models, Pop II and Pop III star formation co-exist over cosmological time-scales, with the total star formation rate density and resulting radiation background strongly dominated by the former before Pop III star formation finally ends. These halos form at most $\sim 10^3$ M$_\odot$ of massive Pop III stars during this phase and typically have absolute magnitudes in the range of $M_\text{AB} = -5 $ to $ -10$. We also briefly discuss how future observations from telescopes such as JWST or WFIRST and 21-cm experiments may be able to constrain unknown parameters in our model such as the IMF, star formation prescription, or the physics of massive Pop III stars., Comment: 16 pages, 13 figures, submitted to MNRAS

Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted HI emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to $k_\parallel$) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the "wedge"), and the spectrally structured emission from HI in the EoR (the "EoR window"). Faraday rotation can inject spectral structure into otherwise smooth polarized foreground emission, which through instrument effects or miscalibration could possibly pollute the EoR window. Using data from the HERA 19-element commissioning array, we investigate the polarization response of this new instrument in the power spectrum domain. We perform a simple image-based calibration based on the unpolarized diffuse emission of the Global Sky Model, and show that it achieves qualitative redundancy between the nominally-redundant baselines of the array and reasonable amplitude accuracy. We construct power spectra of all fully polarized coherencies in all pseudo-Stokes parameters. We compare to simulations based on an unpolarized diffuse sky model and detailed electromagnetic simulations of the dish and feed, confirming that in Stokes I, the calibration does not add significant spectral structure beyond the expected level. Further, this calibration is stable over the 8 days of observations considered. Excess power is seen in the power spectra of the linear polarization Stokes parameters which is not easily attributable to leakage via the primary beam, and results from some combination of residual calibration errors and actual polarized emission. Stokes V is found to be highly discrepant from the expectation of zero power, strongly pointing to the need for more accurate polarized calibration., Comment: 21 pages, 12 figures, submitted to ApJ

The observed scatter in intergalactic Lyman-$\alpha$ opacity at $z \lesssim 6$ requires large-scale fluctuations in the neutral fraction of the intergalactic medium (IGM) after the expected end of reionization. Post-reionization models that explain this scatter invoke fluctuations in either the ionizing ultraviolet background (UVB) or IGM temperature. These models make very different predictions, however, for the relationship between Lyman-$\alpha$ opacity and local density. Here we test these models using Lyman-$\alpha$ emitting galaxies (LAEs) to trace the density field surrounding the longest and most opaque known Lyman-$\alpha$ trough at $z < 6$. Using deep Subaru Hyper Suprime-Cam narrow-band imaging, we find a highly significant deficit of $z \simeq 5.7$ LAEs within 20 Mpc/$h$ of the trough. The results are consistent with a model in which the scatter in Lyman-$\alpha$ opacity near $z \sim 6$ is driven by large-scale UVB fluctuations, and disfavor a scenario in which the scatter is primarily driven by variations in IGM temperature. UVB fluctuations at this epoch present a boundary condition for reionization models, and may help shed light on the nature of the ionizing sources., Comment: 13 pages, 9 figure, submitted to ApJ, comments welcome

The redshifted 21-cm monopole is expected to be a powerful probe of the epoch of the first stars and galaxies ($10, Comment: 13 pages, 8 figures. Updates to this paper include: a new Section 2 (summary of observational strategy), a revised description of DARE instrument in Section 5, and an extended version of Section 6 (extracting the 21-cm spectrum) including 2 new figures (Figs. 5 and 6)

Large variations in the effective optical depth of the He II Ly$\alpha$ forest have been observed at $z\gtrsim2.7$, but the physical nature of these variations is uncertain: either the Universe is still undergoing the process of He II reionization, or the Universe is highly ionized but the He II-ionizing background fluctuates significantly on large scales. In an effort to build upon our understanding of the latter scenario, we present a novel model for the evolution of ionizing background fluctuations. Previous models have assumed the mean free path of ionizing photons to be spatially uniform, ignoring the dependence of that scale on the local ionization state of the intergalactic medium (IGM). This assumption is reasonable when the mean free path is large compared to the average distance between the primary sources of He II-ionizing photons, $\gtrsim L_\star$ quasars. However, when this is no longer the case, the background fluctuations become more severe, and an accurate description of the average propagation of ionizing photons through the IGM requires additionally accounting for the fluctuations in opacity. We demonstrate the importance of this effect by constructing 3D semi-analytic models of the helium ionizing background from $z=2.5$-$3.5$ that explicitly include a spatially varying mean free path of ionizing photons. The resulting distribution of effective optical depths at large scales in the He II Ly$\alpha$ forest is very similar to the latest observations with HST/COS at $2.5 \lesssim z \lesssim 3.5$., Comment: 10 pages, 7 figures

The spin temperature of neutral hydrogen, which determines the optical depth and brightness of the 21 cm line, is determined by the competition between radiative and collisional processes. Here we examine the role of proton-hydrogen collisions in setting the spin temperature. We use recent fully quantum mechanical calculations of the relevant cross sections, which allow us to present accurate results over the entire physically relevant temperature range 1-10,000 K. For kinetic temperatures T_K>100 K, the proton-hydrogen rate coefficient exceeds that for hydrogen-hydrogen collisions by about a factor of two. However, at low temperatures (T_K < 5 K) H-p collisions become several thousand times more efficient than H-H and even more important than H-e^- collisions., submitted to MNRAS, 5 pages, 2 figures, typos corrected

We test the ability of equilibrium galactic disk and one-zone interstellar medium models to describe the physical and emission properties of quasar hosts, submillimeter galaxies, and Lyman-alpha emitters at z>~6. The size, line widths, star formation rates, black hole accretion rates, gas masses and temperatures, and the relationships between these properties are all well-described by our model, and we provide approximate fitting formulae for comparison with future observations. However, comparing our carbon line predictions to observations reveals differences between the ISM at low and high redshifts. Our underestimate of the [CII] line emission indicates either higher star formation efficiencies in high-redshift molecular clouds or less depletion of metals into dust at fixed metallicity. Further, our over-prediction of the CO(6-5)/CO(1-0) ratio suggests that molecular clouds in real high-redshift galaxies have a lower turbulent Mach number and more subthermal CO(6-5) emission than expected owing either to sizes smaller than the local Jeans mass or to a pressure support mechanism other than turbulence., Accepted in MNRAS; 19 pages; 10 figures; 4 tables

Interpretation of He II Ly{\alpha} absorption spectra after the epoch of He II reionization requires knowledge of the He II ionizing background. While past work has modelled the evolution of the average background, the standard cosmological radiative transfer technique assumes a uniform radiation field despite the discrete nature of the (rare) bright quasars that dominate the background. We implement a cosmological radiative transfer model that includes the most recent constraints on the ionizing spectra and luminosity function of quasars and the distribution of IGM absorbers. We also estimate, for the first time, the effects of fluctuations on the evolving continuum opacity in two ways: by incorporating the complete distribution of ionizing background amplitudes into the standard approach, and by explicitly treating the quasars as discrete -- but isolated -- sources. Our model results in a He II ionization rate that evolves steeply with redshift, increasing by a factor ~2 from z=3.0 to z=2.5. This causes rapid evolution in the mean He II Ly{\alpha} optical depth -- as recently observed -- without appealing to the reionization of He II. The observed behaviour could instead result from rapid evolution in the mean free path of ionizing photons as the helium in higher H I column density absorbers becomes fully ionized., Comment: 14 pages, 9 figures. Accepted by MNRAS; significantly modified from previous version

Photonic Doppler Velocimetry is an interferometric technique for measuring the beat frequency of a moving surface, from which the calculated velocity profile of the surface can be used to describe the physical changes the material undergoes after high-impact shock. Such a technique may also be used to characterize the performance of small detonators and determine the time at which the surface began moving. In this work, we develop a semi-automated technique for extracting the time of initial movement from a normalized lineout of the power spectrogram near the offset frequency of each probe. We characterize the response bias of this method and compare with the time of initial movement obtained by hand calculation of the raw voltage data. Results are shown on data from shock experiments such as gas gun setups and explosives-driven flyer plates.

We investigate the effects of Population III stars on the sky-averaged 21-cm background radiation, which traces the collective emission from all sources of ultraviolet and X-ray photons before reionization is complete. While UV photons from PopIII stars can in principle shift the onset of radiative coupling of the 21-cm transition -- and potentially reionization -- to early times, we find that the remnants of PopIII stars are likely to have a more discernible impact on the 21-cm signal than PopIII stars themselves. The X-rays from such sources preferentially heat the IGM at early times, which elongates the epoch of reheating and results in a more gradual transition from an absorption signal to emission. This gradual heating gives rise to broad, asymmetric wings in the absorption signal, which stand in contrast to the relatively sharp, symmetric signals that arise in models treating PopII sources only. A stronger signature of PopIII, in which the position of the absorption minimum becomes inconsistent with PopII-only models, requires extreme star-forming events that may not be physically plausible, lending further credence to predictions of relatively high frequency absorption troughs, $\nu_{\min} \sim 100$ MHz. As a result, though the trough location alone may not be enough to indicate the presence of PopIII, the asymmetric wings should arise even if only a few PopIII stars form in each halo before the transition to PopII star formation occurs, provided that the PopIII IMF is sufficiently top-heavy and at least some PopIII stars form in binaries., Comment: 17 pages, 8 figures, submitted to MNRAS

We develop a radiation pressure-balanced model for the interstellar medium of high-redshift galaxies that describes many facets of galaxy formation at z>~6, including star formation rates and distributions and gas accretion onto central black holes. We first show that the vertical gravitational force in the disk of such a model is dominated by the disk self-gravity supported by the radiation pressure of ionizing starlight on gas. Constraining our model to reproduce the UV luminosity function of Lyman-break galaxies (LBGs), we limit the available parameter-space to wind mass-loading factors 1--4 times the canonical value for momentum-driven winds. We then focus our study by exploring the effects of different angular momentum transport mechanisms in the galactic disk and find that accretion driven by gravitational torques, such as from linear spiral waves or non-linear orbit crossings, can build up black hole masses by z=6 consistent with the canonical M-sigma relation with a duty cycle of unity, while accretion mediated by a local viscosity such as in an alpha-disk results in negligible BH accretion. Both gravitational torque models produce X-ray emission from active galactic nuclei (AGN) in high-redshift LBGs in excess of the estimated contribution from high-mass X-ray binaries. Using a recent analysis of deep Chandra observations by Cowie et al., we can already begin to rule out the most extreme regions of our parameter-space: the inflow velocity of gas through the disk must either be less than one percent of the disk circular velocity or the X-ray luminosity of the AGN must be substantially obscured. Moderately deeper future observations or larger sample sizes will be able to probe the more reasonable range of angular momentum transport models and obscuring geometries.

We use a new method to model fluctuations of the Lyman-Werner (LW) and Lyman-alpha radiation backgrounds at high redshift. At these early epochs the backgrounds are symptoms of a universe newly lit with its first stars. LW photons (11.5-13.6 eV) are of particular interest because they dissociate molecular hydrogen, the primary coolant in the first minihalos. By using a variation of the halo model, we efficiently generate power spectra for any choice of radiation background. We find that the LW power spectrum typically traces the matter power spectrum at large scales but turns over at the scale corresponding to the effective `horizon' of LW photons (~100 comoving Mpc), unless the sources are extremely rare. The series of horizons that characterize the Lyman-alpha flux profile shape the fluctuations of that background in a similar fashion, though those imprints are washed out once one considers fluctuations in the brightness temperature of the 21-cm signal. The Lyman-alpha background strongly affects the redshifted 21-cm signal at just about the time the LW background begins to dissociate molecular hydrogen, so measuring that background's properties will reveal important information about the transition from early Population III stars to more normal stars. Around this time we find that fluctuations in the LW background are weak; the fractional standard deviation is less than ~0.5 on scales > 10 cMpc, only rising to be of order unity on scales < 1 cMpc. This should not lead to substantial spatial fluctuations in molecular hydrogen content, except at the earliest times. Even then, most halos form far from other sources, so the transition from star formation in low-mass to high-mass halos is rather homogeneous across the universe.

We introduce a powerful semi-numeric modeling tool, 21cmFAST, designed to efficiently simulate the cosmological 21-cm signal. Our code generates 3D realizations of evolved density, ionization, peculiar velocity, and spin temperature fields, which it then combines to compute the 21-cm brightness temperature. Although the physical processes are treated with approximate methods, we compare our results to a state-of-the-art large-scale hydrodynamic simulation, and find good agreement on scales pertinent to the upcoming observations (>~ 1 Mpc). The power spectra from 21cmFAST agree with those generated from the numerical simulation to within 10s of percent, down to the Nyquist frequency. We show results from a 1 Gpc simulation which tracks the cosmic 21-cm signal down from z=250, highlighting the various interesting epochs. Depending on the desired resolution, 21cmFAST can compute a redshift realization on a single processor in just a few minutes. Our code is fast, efficient, customizable and publicly available, making it a useful tool for 21-cm parameter studies.

Lyman-limit absorption systems can play many important roles during and after cosmological reionization. Unfortunately, due to the prohibitively large dynamic range required, it is impossible to self-consistently include these systems in cosmological simulations. Using fast and versatile seminumeric simulations, we systematically explore the spatial distribution of absorption systems during and following reionization. We self-calibrate the resulting number of absorbers to the mean free path (mfp) of the ionizing ultraviolet background (UVB), and present results at a given mfp and neutral hydrogen fraction. We use a simple optical depth criterion to identify the locations of absorbers. Our approach is fairly robust to uncertainties such as missing subgrid structure. Unlike at lower redshifts where the UVB is relatively uniform, at higher redshifts the fluctuations in the UVB and the H II morphology of reionization can drive the large-scale distribution of absorption systems. Specifically, we find that absorbers are highly correlated with the density field on small scales, and then become anti-correlated with the UVB on large scales. After reionization, the large-scale power spectrum of the absorbers traces the UVB power spectrum, which can be predicted with a simple analytic extension of the halo model. During reionization, absorbers tend to preferentially lie inside overdensities (i.e. filaments) of the recently ionized intergalactic medium. Absorbers may also dominate the small-scale (k 1M pc −1 ) 21-cm power during and after reionization. Conversely, they smooth the contrast on moderate scales. Once the H II regions grow to surpass the mfp, the absorbers add to the large-scale 21-cm power. Our results should prove useful in interpreting future observations of the reionization epoch.

We examine large-scale fluctuations in the He?II Ly? forest transmission during and after He?II reionization. We use a simple Monte Carlo model to distribute quasars throughout a large volume and compute the resulting radiation field along one-dimensional skewers. In agreement with previous studies, we find that the rarity of these sources induces order unity fluctuations in the mean optical depth after reionization, even when averaged over large segments (~10-100Mpc across). We compare our models to existing data along five He?II Ly? forest lines of sight spanning z ~ 2-3.2. The large cosmic variance contained in our model plausibly explains many of the observed fluctuations at z 2.5. But our models cannot accommodate the large fluctuations toward high optical depths on ~30Mpc scales observed at z ~ 2.7-2.9, and the measured optical depths (?eff 4) at z>2.9 are difficult to explain with a smoothly evolving mean radiation field. In order to better understand these data, we construct a toy model of He?II reionization, in which we assume that regions with the smallest radiation fields in a post-reionization universe (or farthest from strong ionizing sources) are completely dark during reionization. The observed fluctuations fit much more comfortably into this model, and we therefore argue that, according to present data, He?II reionization does not complete until z 2.9.

We use an analytic model to study how inhomogeneous hydrogen reionization affects the temperature distribution of the intergalactic medium (IGM). During this process, the residual energy of each ionizing photon is deposited in the IGM as heat, increasing its temperature to 20,000-30,000 K; subsequent expansion of the Universe then cools the gas. Because reionization most likely proceeds from high to low densities, underdense voids are ionized last, have less time to cool, and are (on average) warmer than mean-density gas immediately after reionization is complete (an "inverted" equation of state). From this initial configuration, the low-density gas cools quickly and eventually returns to a more normal equation of state. The rapidly evolving temperature introduces systematic uncertainties in measurements of the ionizing background at z~6. For example, late reionization implies rapid cooling, so that the ionizing background would have to evolve even more rapidly at z ~5-6 than typically claimed. This degeneracy is difficult to disentangle, because the Lyman-alpha forest probes only a narrow range in densities (over which the gas is nearly isothermal). However, higher Lyman-series transitions probe wider density ranges, sampling different effective temperatures, and offer a new way to measure the IGM equation of state that should work where nearly saturated absorption precludes other methods. This will help to separate evolution in temperature from that in the ionizing background. While more detailed study with hydrodynamic simulations is needed, we show that such measurements could potentially distinguish early and late reionization using only a handful of lines of sight.

We use a combination of analytic and numerical arguments to consider the impact of quasar photo-heating during HeII reionisation on the thermal evolution of the intergalactic medium (IGM). We demonstrate that rapid (\Delta z< 0.1-0.2), strong (\Delta T > 10^4 K) photo-heating is difficult to achieve across the entire IGM unless quasar spectra are significantly harder than implied by current observational constraints. Although filtering of intrinsic quasar radiation through dense regions in the IGM does increase the mean excess energy per HeII photo-ionisation, it also weakens the radiation intensity and lowers the photo-ionisation rate, preventing rapid heating over time intervals shorter than the local photo-ionisation timescale. Moreover, the hard photons responsible for the strongest heating are more likely to deposit their energy inside dense clumps. The abundance of such clumps is, however, uncertain and model-dependent, leading to a fairly large uncertainty in the photo-heating rates. Nevertheless, although some of the IGM may be exposed to a hardened and weakened ionising background for long periods, most of the IGM must instead be reionised by the more abundant, softer photons and with accordingly modest heating rates (\Delta T < 10^4 K). The repeated ionisation of fossil quasar HeIII regions does not increase the net heating because the recombination times in these regions typically exceed the IGM cooling times and the average time lag between successive rounds of quasar activity. Detailed line-of-sight radiative transfer simulations confirm these expectations and predict a rich thermal structure in the IGM during HeII reionisation. [Abridged], Comment: 20 pages, 6 figures, accepted by MNRAS

Weak gravitational lensing provides a potentially powerful method for the detection of clusters. In addition to cluster candidates, a large number of objects with possibly no optical or X-ray component have been detected in shear-selected samples. We develop an analytic model to investigate the claim of Weinberg & Kamionkowski (2002) that unvirialised protoclusters account for a significant number of these so-called "dark" lenses. In our model, a protocluster consists of a small virialised region surrounded by in-falling matter. We find that, in order for a protocluster to simultaneously escape X-ray detection and create a detectable weak lensing signal, it must have a small virial mass (~10^{13} \Msun) and large total mass (~ 10^{15} \Msun), with a relatively flat density profile outside of the virial radius. Such objects would be characterized by rising tangential shear profiles well beyond the virial radius. We use a semi-analytic approach based on the excursion set formalism to estimate the abundance of lensing protoclusters with a low probability of X-ray detection. We find that they are extremely rare, accounting for less than 0.4 per cent of the total lenses in a survey with background galaxy density n = 30 arcmin^{-2} and an intrinsic ellipticity dispersion of 0.3. We conclude that lensing protoclusters with undetectable X-Ray luminosities are too rare to account for a significant number of dark lenses., 18 pages, 10 figures, version accepted by MNRAS (minor changes in response to referee)

The ionizing background of cosmic hydrogen is an important probe of the sources and absorbers of ionizing radiation in the post-reionization universe. Previous studies show that the ionization rate should be very sensitive to changes in the source population: as the emissivity rises, absorbers shrink in size, increasing the ionizing mean free path and, hence, the ionizing background. By contrast, observations of the ionizing background find a very flat evolution from z~2-5, before falling precipitously at z~6. We resolve this puzzling discrepancy by pointing out that, at z~2-5, optically thick absorbers are associated with the same collapsed halos that host ionizing sources. Thus, an increasing abundance of galaxies is compensated for by a corresponding increase in the absorber population, which moderates the instability in the ionizing background. However, by z~5-6, gas outside of halos dominates the absorption, the coupling between sources and absorbers is lost, and the ionizing background evolves rapidly. Our halo based model reproduces observations of the ionizing background, its flatness and sudden decline, as well as the redshift evolution of the ionizing mean free path. Our work suggests that, through much of their history, both star formation and photoelectric opacity in the universe track halo growth., 16 pages, 7 figures, submitted to MNRAS, revised methods leading to same but more convincing conclusions

One of the most exciting probes of the early phases of structure formation and reionization is the “spin-flip” line of neutral hydrogen, with a rest wavelength of 21 cm. This chapter introduces the physics of this transition and the astrophysical parameters upon which it depends, including discussions of the radiation fields that permeate the intergalactic medium that fix the brightness of this transition. We describe the critical points in the evolution of the 21-cm background and focus on the sky-averaged brightness and the power spectrum as representative measurements. Finally, we include a discussion of observations and the challenges they face in the near future.

Recent observations have shown that the scatter in opacities among coeval segments of the Lyman-alpha forest increases rapidly at z > 5. In this paper, we assess whether the large scatter can be explained by fluctuations in the ionizing background in the post-reionization intergalactic medium. We find that matching the observed scatter at z ~ 5.5 requires a short spatially averaged mean free path of < 15 comoving Mpc/h, a factor of > 3 shorter than direct measurements at z ~ 5.2. We argue that such rapid evolution in the mean free path is difficult to reconcile with our measurements of the global H I photoionization rate, which stay approximately constant over the interval z ~ 4.8 - 5.5. However, we also show that measurements of the mean free path at z > 5 are likely biased towards higher values by the quasar proximity effect. This bias can reconcile the short values of the mean free path that are required to explain the large scatter in opacities. We discuss the implications of this scenario for cosmological reionization. Finally, we investigate whether other statistics applied to the z > 5 Lyman-alpha forest can shed light on the origin of the scatter. Compared to a model with a uniform ionizing background, models that successfully account for the scatter lead to enhanced power in the line-of-sight flux power spectrum on scales k < 0.1 h/Mpc. We find tentative evidence for this enhancement in observations of the high-redshift Lyman-alpha forest., Matches version published by MNRAS with clarifications and expanded discussion

Motivated by recent progress in studies of the high-$z$ Universe, we build a new model for the global 21-cm signal that is explicitly calibrated to measurements of the galaxy luminosity function (LF) and further tuned to match the Thomson scattering optical depth of the cosmic microwave background, $\tau_e$. Assuming that the $z \lesssim 8$ galaxy population can be smoothly extrapolated to higher redshifts, the recent decline in best-fit values of $\tau_e$ and the inefficient heating induced by X-ray binaries (HMXBs; the presumptive sources of the X-ray background at high-$z$) imply that the entirety of cosmic reionization and reheating occurs at redshifts $z \lesssim 12$. In contrast to past global 21-cm models, whose $z \sim 20$ ($\nu \sim 70$ MHz) absorption features and strong $\sim 25$ mK emission features were driven largely by the assumption of efficient early star-formation and X-ray heating, our new fiducial model peaks in absorption at $\nu \sim 110$ MHz at a depth of $\sim -160$ mK and has a negligible emission component. As a result, a strong emission signal would provide convincing evidence that HMXBs are not the only drivers of cosmic reheating. Shallow absorption troughs should accompany strong heating scenarios, but could also be caused by a low escape fraction of Lyman-Werner photons. Generating signals with troughs at $\nu \lesssim 95$ MHz requires a floor in the star-formation efficiency in halos below $\sim 10^{9} M_{\odot}$, which is equivalent to steepening the faint-end of the galaxy LF. These findings demonstrate that the global 21-cm signal is a powerful complement to current and future galaxy surveys and efforts to better understand the interstellar medium in high-$z$ galaxies., Comment: 17 pages, 9 figures, in press

The Hydrogen Epoch of Reionization Array (HERA) is a staged experiment to measure 21 cm emission from the primordial intergalactic medium (IGM) throughout cosmic reionization ($z=6-12$), and to explore earlier epochs of our Cosmic Dawn ($z\sim30$). During these epochs, early stars and black holes heated and ionized the IGM, introducing fluctuations in 21 cm emission. HERA is designed to characterize the evolution of the 21 cm power spectrum to constrain the timing and morphology of reionization, the properties of the first galaxies, the evolution of large-scale structure, and the early sources of heating. The full HERA instrument will be a 350-element interferometer in South Africa consisting of 14-m parabolic dishes observing from 50 to 250 MHz. Currently, 19 dishes have been deployed on site and the next 18 are under construction. HERA has been designated as an SKA Precursor instrument. In this paper, we summarize HERA's scientific context and provide forecasts for its key science results. After reviewing the current state of the art in foreground mitigation, we use the delay-spectrum technique to motivate high-level performance requirements for the HERA instrument. Next, we present the HERA instrument design, along with the subsystem specifications that ensure that HERA meets its performance requirements. Finally, we summarize the schedule and status of the project. We conclude by suggesting that, given the realities of foreground contamination, current-generation 21 cm instruments are approaching their sensitivity limits. HERA is designed to bring both the sensitivity and the precision to deliver its primary science on the basis of proven foreground filtering techniques, while developing new subtraction techniques to unlock new capabilities. The result will be a major step toward realizing the widely recognized scientific potential of 21 cm cosmology., Comment: 26 pages, 24 figures, 2 tables

We examine the effect of inhomogeneous reionization on the galaxy power spectrum and the consequences for probing dark energy. To model feedback during reionization, we apply an ansatz setting the galaxy overdensity proportional to the underlying ionization field. Thus, inhomogeneous reionization may leave an imprint in the galaxy power spectrum. We evolve this imprint to low redshift and use the Fisher-matrix formalism to assess the effect on parameter estimation. We show that a combination of low- (z=0.3) and high- (z=3) redshift galaxy surveys can constrain the size of cosmological HII regions during reionization. This imprint can also cause confusion when using baryon oscillations or other features of the galaxy power spectrum to probe the dark energy. We show that when bubbles are large, and hence detectable, our ability to constrain w can be degraded by up to 50%. When bubbles are small, the imprint has little or no effect on measuring dark-energy parameters., Comment: 9 pages, 8 figures, version accepted by MNRAS

We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the intergalactic medium, to detail the growth of HII regions around high redshift galaxies. Our simulation tracks 1024^3 dark matter particles, in a box of co-moving side length 65.6 Mpc/h. This large volume allows us to accurately characterize the size distribution of HII regions throughout most of the reionization process. At the same time, our simulation resolves many of the small galaxies likely responsible for reionization. It confirms a picture anticipated by analytic models: HII regions grow collectively around highly-clustered sources, and have a well-defined characteristic size, which evolves from a sub-Mpc scale at the beginning of reionization to R > 10 Mpc towards the end. We present a detailed statistical description of our results, and compare them with a numerical scheme based on the analytic model by Furlanetto, Zaldarriaga, and Hernquist. We find that the analytic calculation reproduces the size distribution of HII regions and the 21 cm power spectrum of the radiative transfer simulation remarkably well. The ionization field from the simulation, however, has more small scale structure than the analytic calculation, owing to Poisson scatter in the simulated abundance of galaxies on small scales. We propose and validate a simple scheme to incorporate this scatter into our calculations. Our results suggest that analytic calculations are sufficiently accurate to aid in predicting and interpreting the results of future 21 cm surveys. In particular, our fast numerical scheme is useful for forecasting constraints from future 21 cm surveys, and in constructing mock surveys to test data analysis procedures., 15 pages, 11 figures, submitted to ApJ