Showing 1–50 of 112 results for author: Furlanetto, S R .

The AURORA Survey: A New Era of Emission-line Diagrams with JWST/NIRSpec

Authors: Alice E. Shapley, Ryan L. Sanders, Michael W. Topping, Naveen A. Reddy, Danielle A. Berg, Rychard J. Bouwens, Gabriel Brammer, Adam C. Carnall, Fergus Cullen, Romeel Davé, James S. Dunlop, Richard S. Ellis, N. M. Förster Schreiber, Steven R . Furlanetto, Karl Glazebrook, Garth D. Illingworth, Tucker Jones, Mariska Kriek, Derek J. McLeod, Ross J. McLure, Desika Narayanan, Pascal Oesch, Anthony J. Pahl, Max Pettini, Daniel Schaerer , et al. (6 additional authors not shown)

Abstract: We present results on the emission-line properties of z=1.4-7.5 star-forming galaxies in the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) Cycle 1 JWST/NIRSpec program. Based on its depth, continuous wavelength coverage from 1--5 microns, and medium spectral resolution (R~1000), AURORA includes detections of a large suite of nebular emission lines spanning a broad… ▽ More We present results on the emission-line properties of z=1.4-7.5 star-forming galaxies in the Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (AURORA) Cycle 1 JWST/NIRSpec program. Based on its depth, continuous wavelength coverage from 1--5 microns, and medium spectral resolution (R~1000), AURORA includes detections of a large suite of nebular emission lines spanning a broad range in rest wavelength. We investigate the locations of AURORA galaxies in multiple different emission-line diagrams, including traditional "BPT" diagrams of [OIII]/Hbeta vs. [NII]/Halpha, [SII]/Halpha, and [OI]/Halpha, and the "ionization-metallicity" diagram of [OIII]/[OII] (O32) vs. ([OIII]+[OII])/Hbeta (R23). We also consider a bluer rest-frame "ionization-metallicity" diagram introduced recently to characterize z>10 galaxies: [NeIII]/[OII] vs. ([NeIII]+[OII])/Hdelta; as well as longer-wavelength diagnostic diagrams extending into the rest-frame near-IR: [OIII]/Hbeta vs. [SIII]/[SII] (S32); and HeI/Pagamma and [SIII]/Pagamma vs. [FeII]/Pabeta. With a significant boost in signal-to-noise and large, representative samples of individual galaxy detections, the AURORA emission-line diagrams presented here definitively confirm a physical picture in which chemically-young, alpha-enhanced, massive stars photoionize the ISM in distant galaxies with a harder ionizing spectrum at fixed nebular metallicity than in their z~0 counterparts. We also uncover previously unseen evolution prior to z~2 in the [OIII]/Hbeta vs. [NII]/Halpha diagram, which motivates deep NIRSpec observations at even higher redshift. Finally, we present the first statistical sample of rest-frame near-IR emission-line diagnostics in star-forming galaxies at high redshift. In order to truly interpret rest-frame near-IR line ratios including [FeII], we must obtain better constraints on dust depletion in the high-redshift ISM. △ Less

Submitted 2 July, 2024; v1 submitted 28 June, 2024; originally announced July 2024.

Comments: 17 pages, 8 figures, submitted to ApJ

The Predicament of Absorption-dominated Reionization II: Observational Estimate of the Clumping Factor at the End of Reionization

Authors: Frederick B. Davies, Sarah E. I. Bosman, Steven R. Furlanetto

Abstract: The history of reionization reflects the cumulative injection of ionizing photons by sources and the absorption of ionizing photons by sinks. The latter process is traditionally described in terms of a "clumping factor" which encodes the average quadratic increase in the recombination rate of dense gas within the cosmic web. The recent measurement of a short mean free path of ionizing photons from… ▽ More The history of reionization reflects the cumulative injection of ionizing photons by sources and the absorption of ionizing photons by sinks. The latter process is traditionally described in terms of a "clumping factor" which encodes the average quadratic increase in the recombination rate of dense gas within the cosmic web. The recent measurement of a short mean free path of ionizing photons from stacked quasar spectra at $z\simeq6$ has placed the importance of sinks under increased scrutiny, but its connection to the recombination rate is not immediately obvious. Here we present analytic arguments to connect the clumping factor to the mean free path by invoking ionization equilibrium within the ionized phase of the intergalactic medium at the end of (and after) reionization. We find that the latest mean free path and hydrogen photoionization rate measurements at $z=5-6$ imply a global clumping factor $C\approx12$, much higher than previous determinations from radiation-hydrodynamic simulations of the reionization process. Similar values of $C$ are also derived when applying the same procedure to observations at $2<z<5$. Compared to the traditional assumption of $C=3$, high-redshift galaxies must produce roughly twice as many ionizing photons ($\approx3$ photons per baryon) to reionize the universe by $z\sim6$. This additional requirement on the ionizing photon budget may help to reconcile the reionization history with JWST observations that suggest a far greater output of ionizing photons by the most distant galaxy populations. △ Less

Submitted 11 July, 2024; v1 submitted 26 June, 2024; originally announced June 2024.

Comments: 11 pages, 5+1 figures. Typos fixed, discussion improved, appendix added

Report number: NORDITA 2024-021

Investigating Mutual Coupling in the Hydrogen Epoch of Reionization Array and Mitigating its Effects on the 21-cm Power Spectrum

Authors: E. Rath, R. Pascua, A. T. Josaitis, A. Ewall-Wice, N. Fagnoni, E. de Lera Acedo, Z. E. Martinot, Z. Abdurashidova, T. Adams, J. E. Aguirre, R. Baartman, A. P. Beardsley, L. M. Berkhout, G. Bernardi, T. S. Billings, J. D. Bowman, P. Bull, J. Burba, R. Byrne, S. Carey, K. -F. Chen, S. Choudhuri, T. Cox, D. R. DeBoer, M. Dexter , et al. (56 additional authors not shown)

Abstract: Interferometric experiments designed to detect the highly redshifted 21-cm signal from neutral hydrogen are producing increasingly stringent constraints on the 21-cm power spectrum, but some k-modes remain systematics-dominated. Mutual coupling is a major systematic that must be overcome in order to detect the 21-cm signal, and simulations that reproduce effects seen in the data can guide strategi… ▽ More Interferometric experiments designed to detect the highly redshifted 21-cm signal from neutral hydrogen are producing increasingly stringent constraints on the 21-cm power spectrum, but some k-modes remain systematics-dominated. Mutual coupling is a major systematic that must be overcome in order to detect the 21-cm signal, and simulations that reproduce effects seen in the data can guide strategies for mitigating mutual coupling. In this paper, we analyse 12 nights of data from the Hydrogen Epoch of Reionization Array and compare the data against simulations that include a computationally efficient and physically motivated semi-analytic treatment of mutual coupling. We find that simulated coupling features qualitatively agree with coupling features in the data; however, coupling features in the data are brighter than the simulated features, indicating the presence of additional coupling mechanisms not captured by our model. We explore the use of fringe-rate filters as mutual coupling mitigation tools and use our simulations to investigate the effects of mutual coupling on a simulated cosmological 21-cm power spectrum in a "worst case" scenario where the foregrounds are particularly bright. We find that mutual coupling contaminates a large portion of the "EoR Window", and the contamination is several orders-of-magnitude larger than our simulated cosmic signal across a wide range of cosmological Fourier modes. While our fiducial fringe-rate filtering strategy reduces mutual coupling by roughly a factor of 100 in power, a non-negligible amount of coupling cannot be excised with fringe-rate filters, so more sophisticated mitigation strategies are required. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 19 pages, 12 figures, submitted to MNRAS

A hidden population of active galactic nuclei can explain the overabundance of luminous $z>10$ objects observed by JWST

Authors: Sahil Hegde, Michael M. Wyatt, Steven R. Furlanetto

Abstract: The first wave of observations with JWST has revealed a striking overabundance of luminous galaxies at early times ($z>10$) compared to models of galaxies calibrated to pre-JWST data. Early observations have also uncovered a large population of supermassive black holes (SMBHs) at $z>6$. Because many of the high-$z$ objects appear extended, the contribution of active galactic nuclei (AGNs) to the t… ▽ More The first wave of observations with JWST has revealed a striking overabundance of luminous galaxies at early times ($z>10$) compared to models of galaxies calibrated to pre-JWST data. Early observations have also uncovered a large population of supermassive black holes (SMBHs) at $z>6$. Because many of the high-$z$ objects appear extended, the contribution of active galactic nuclei (AGNs) to the total luminosity has been assumed to be negligible. In this work, we use a semi-empirical model for assigning AGNs to galaxies to show that active galaxies can boost the stellar luminosity function (LF) enough to solve the overabundance problem while simultaneously remaining consistent with the observed morphologies of high-$z$ sources. We construct a model for the composite AGN+galaxy LF by connecting dark matter halo masses to galaxy and SMBH masses and luminosities, accounting for dispersion in the mapping between host galaxy and SMBH mass and luminosity. By calibrating the model parameters -- which characterize the $M_\bullet-M_\star$ relation -- to a compilation of $z>10$ JWST UVLF data, we show that AGN emission can account for the excess luminosity under a variety of scenarios, including one where 10\% of galaxies host BHs of comparable luminosities to their stellar components. Using a sample of simulated objects and real observations, we demonstrate that such low-luminosity AGNs can be `hidden' in their host galaxies and be missed in common morphological analyses. We find that for this explanation to be viable, our model requires a population of BHs that are overmassive ($M_\bullet/M_\star\sim10^{-2}$) with respect to their host galaxies compared to the local relation and are more consistent with the observed relation at $z=4-8$. We explore the implications of this model for BH seed properties and comment on observational diagnostics necessary to further investigate this explanation. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 34 pages, 12 figures; submitted to JCAP

Reionization after JWST: a photon budget crisis?

Authors: Julian B. Muñoz, Jordan Mirocha, John Chisholm, Steven R. Furlanetto, Charlotte Mason

Abstract: New James Webb Space Telescope (JWST) observations are revealing the first galaxies to be prolific producers of ionizing photons, which we argue gives rise to a tension between different probes of reionization. For hydrogen reionization to proceed there must be enough ionizing photons for all of the hydrogen atoms, including their recombinations. Over the last two decades a consensus has emerged w… ▽ More New James Webb Space Telescope (JWST) observations are revealing the first galaxies to be prolific producers of ionizing photons, which we argue gives rise to a tension between different probes of reionization. For hydrogen reionization to proceed there must be enough ionizing photons for all of the hydrogen atoms, including their recombinations. Over the last two decades a consensus has emerged where star-forming galaxies are able to generate these photons, given reasonable values for their number densities, ionizing efficiencies $ξ_{\rm ion}$ (per unit UV luminosity), and escape fractions $f_{\rm esc}$. However, new JWST observations infer high values of $ξ_{\rm ion}$ during reionization and an enhanced abundance of earlier ($z\gtrsim 9$) galaxies, dramatically increasing the number of ionizing photons produced at high $z$. Simultaneously, recent low-$z$ studies predict significant escape fractions for faint reionization-era galaxies. Put together, we show that the galaxies we have directly observed ($M_{\rm UV} < -15$) not only can drive reionization, but would end it too early. That is, our current galaxy observations, taken at face value, imply an excess of ionizing photons and thus a process of reionization in tension with the cosmic microwave background (CMB) and Lyman-$α$ forest. Considering galaxies down to $M_{\rm UV}\approx -11$, below current observational limits, only worsens this tension, requiring on average $f_{\rm esc}\approx 3\%$, far lower than expected for early galaxies from post-reionization studies. We discuss possible avenues to resolve this photon budget crisis, including missing astrophysical or observational selection effects, as well as enhanced recombinations. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 6+3 pages, 3+2 figures, comments welcome

A demonstration of the effect of fringe-rate filtering in the Hydrogen Epoch of Reionization Array delay power spectrum pipeline

Authors: Hugh Garsden, Philip Bull, Mike Wilensky, Zuhra Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Lindsay M. Berkhout, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Steven Carey, Chris L. Carilli, Kai-Feng Chen, Carina Cheng, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter , et al. (72 additional authors not shown)

Abstract: Radio interferometers targeting the 21cm brightness temperature fluctuations at high redshift are subject to systematic effects that operate over a range of different timescales. These can be isolated by designing appropriate Fourier filters that operate in fringe-rate (FR) space, the Fourier pair of local sidereal time (LST). Applications of FR filtering include separating effects that are correl… ▽ More Radio interferometers targeting the 21cm brightness temperature fluctuations at high redshift are subject to systematic effects that operate over a range of different timescales. These can be isolated by designing appropriate Fourier filters that operate in fringe-rate (FR) space, the Fourier pair of local sidereal time (LST). Applications of FR filtering include separating effects that are correlated with the rotating sky vs. those relative to the ground, down-weighting emission in the primary beam sidelobes, and suppressing noise. FR filtering causes the noise contributions to the visibility data to become correlated in time however, making interpretation of subsequent averaging and error estimation steps more subtle. In this paper, we describe fringe rate filters that are implemented using discrete prolate spheroidal sequences, and designed for two different purposes -- beam sidelobe/horizon suppression (the `mainlobe' filter), and ground-locked systematics removal (the `notch' filter). We apply these to simulated data, and study how their properties affect visibilities and power spectra generated from the simulations. Included is an introduction to fringe-rate filtering and a demonstration of fringe-rate filters applied to simple situations to aid understanding. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 21 pages, 18 figures, submitted to Monthly Notices of the Royal Astronomical Society

Dust in High-Redshift Galaxies: Reconciling UV Attenuation and IR Emission

Authors: Jianhan Zhao, Steven R. Furlanetto

Abstract: Dust is a key component of galaxies, but its properties during the earliest eras of structure formation remain elusive. Here we present a simple semi-analytic model of the dust distribution in galaxies at $z \gtrsim 5$. We calibrate the free parameters of this model to estimates of the UV attenuation (using the IRX-$β$ relation between infrared emission and the UV spectral slope) and to ALMA measu… ▽ More Dust is a key component of galaxies, but its properties during the earliest eras of structure formation remain elusive. Here we present a simple semi-analytic model of the dust distribution in galaxies at $z \gtrsim 5$. We calibrate the free parameters of this model to estimates of the UV attenuation (using the IRX-$β$ relation between infrared emission and the UV spectral slope) and to ALMA measurements of dust emission. We find that the observed dust emission requires that most of the dust expected in these galaxies is retained (assuming a similar yield to lower-redshift sources), but if the dust is spherically distributed, the modest attenuation requires that it be significantly more extended than the stars. Interestingly, the retention fraction is larger for less massive galaxies in our model. However, the required radius is a significant fraction of the host's virial radius and is larger than the estimated extent of dust emission from stacked high-$z$ galaxies. These can be reconciled if the dust is distributed anisotropically, with typical covering fractions of $\sim 0.2$--0.7 in bright galaxies and $\lesssim 0.1$ in fainter ones. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: 32 pages, 15 figures, submitted to JCAP

Hydrogen Epoch of Reionization Array (HERA) Phase II Deployment and Commissioning

Authors: Lindsay M. Berkhout, Daniel C. Jacobs, Zuhra Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steven Carey, Chris L. Carilli, Kai-Feng Chen, Carina Cheng, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon , et al. (71 additional authors not shown)

Abstract: This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system an… ▽ More This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system and discuss progress on commissioning and future upgrades. As HERA is a designated Square Kilometer Array (SKA) pathfinder instrument, we also show a number of "case studies" that investigate systematics seen while commissioning the phase II system, which may be of use in the design and operation of future arrays. Common pathologies are likely to manifest in similar ways across instruments, and many of these sources of contamination can be mitigated once the source is identified. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Journal ref: PASP 2024 136 045002

Effects of Bursty Star Formation on [C II] Line Intensity Mapping of High-redshift Galaxies

Authors: Lun-Jun Liu, Guochao Sun, Tzu-Ching Chang, Steven R. Furlanetto, Charles M. Bradford

Abstract: Bursty star formation -- a key prediction for high-redshift galaxies from cosmological simulations explicitly resolving stellar feedback in the interstellar medium -- has recently been observed to prevail among galaxies at redshift $z \gtrsim 6$. Line intensity mapping (LIM) of the 158 $μ$m [C II] line as a star formation rate indicator offers unique opportunities to tomographically constrain cosm… ▽ More Bursty star formation -- a key prediction for high-redshift galaxies from cosmological simulations explicitly resolving stellar feedback in the interstellar medium -- has recently been observed to prevail among galaxies at redshift $z \gtrsim 6$. Line intensity mapping (LIM) of the 158 $μ$m [C II] line as a star formation rate indicator offers unique opportunities to tomographically constrain cosmic star formation at high redshift, as an alternative to observations of individually detected galaxies. To understand effects of bursty star formation on [C II] LIM, which remain unexplored in previous studies, we present an analytic modeling framework for high-$z$ galaxy formation and [C II] LIM signals that accounts for bursty star formation histories induced by delayed supernova feedback. We use it to explore and characterize how bursty star formation can impact and thus complicate the interpretation of the [C II] luminosity function and power spectrum. Our simple analytic model indicates that bursty star formation is most important for low halo masses, and in the power spectrum it can create a substantial excess in the large-scale clustering term. This distortion results in a power spectrum shape which cannot be explained by invoking a mass-independent scatter. We conclude that burstiness must be accounted for when modeling and analyzing [C II] datasets from the early universe, and that in the extreme, the signature of burstiness may be detectable with first-generation experiments such as TIME, CONCERTO, and CCAT-DSS. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 12 pages, 4 figures, submitted to AAS journals; comments welcome

matvis: A matrix-based visibility simulator for fast forward modelling of many-element 21 cm arrays

Authors: Piyanat Kittiwisit, Steven G. Murray, Hugh Garsden, Philip Bull, Christopher Cain, Aaron R. Parsons, Jackson Sipple, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Lindsay M. Berkhout, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Steven Carey, Chris L. Carilli, Kai-Feng Chen, Carina Cheng , et al. (73 additional authors not shown)

Abstract: Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability… ▽ More Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability to perform high-fidelity simulations of the kinds of data that are produced by the large, many-element, radio interferometric arrays that have been purpose-built for these studies. The large scale of these arrays presents a computational challenge, as one must simulate a detailed sky and instrumental model across many hundreds of frequency channels, thousands of time samples, and tens of thousands of baselines for arrays with hundreds of antennas. In this paper, we present a fast matrix-based method for simulating radio interferometric measurements (visibilities) at the necessary scale. We achieve this through judicious use of primary beam interpolation, fast approximations for coordinate transforms, and a vectorised outer product to expand per-antenna quantities to per-baseline visibilities, coupled with standard parallelisation techniques. We validate the results of this method, implemented in the publicly-available matvis code, against a high-precision reference simulator, and explore its computational scaling on a variety of problems. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 25 pages, 20 figures, submitted to RAS Techniques and Instruments, matvis is publicly available at https://github.com/HERA-Team/matvis

Bayesian estimation of cross-coupling and reflection systematics in 21cm array visibility data

Authors: Geoff G. Murphy, Philip Bull, Mario G. Santos, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Christopher Cain, Steven Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon, Nico Eksteen , et al. (54 additional authors not shown)

Abstract: Observations with radio arrays that target the 21-cm signal originating from the early Universe suffer from a variety of systematic effects. An important class of these are reflections and spurious couplings between antennas. We apply a Hamiltonian Monte Carlo sampler to the modelling and mitigation of these systematics in simulated Hydrogen Epoch of Reionisation Array (HERA) data. This method all… ▽ More Observations with radio arrays that target the 21-cm signal originating from the early Universe suffer from a variety of systematic effects. An important class of these are reflections and spurious couplings between antennas. We apply a Hamiltonian Monte Carlo sampler to the modelling and mitigation of these systematics in simulated Hydrogen Epoch of Reionisation Array (HERA) data. This method allows us to form statistical uncertainty estimates for both our models and the recovered visibilities, which is an important ingredient in establishing robust upper limits on the Epoch of Reionisation (EoR) power spectrum. In cases where the noise is large compared to the EoR signal, this approach can constrain the systematics well enough to mitigate them down to the noise level for both systematics studied. Where the noise is smaller than the EoR, our modelling can mitigate the majority of the reflections with there being only a minor level of residual systematics, while cross-coupling sees essentially complete mitigation. Our approach performs similarly to existing filtering/fitting techniques used in the HERA pipeline, but with the added benefit of rigorously propagating uncertainties. In all cases it does not significantly attenuate the underlying signal. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 19 pages, 14 figures, submitted to MNRAS

Direct Optimal Mapping Image Power Spectrum and its Window Functions

Authors: Zhilei Xu, Honggeun Kim, Jacqueline N. Hewitt, Kai-Feng Chen, Nicholas S. Kern, Eleanor Rath, Ruby Byrne, Adélie Gorce, Robert Pascua, Zachary E. Martinot, Joshua S. Dillon, Bryna J. Hazelton, Adrian Liu, Miguel F. Morales, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman , et al. (57 additional authors not shown)

Abstract: The key to detecting neutral hydrogen during the epoch of reionization (EoR) is to separate the cosmological signal from the dominating foreground radiation. We developed direct optimal mapping (DOM) to map interferometric visibilities; it contains only linear operations, with full knowledge of point spread functions from visibilities to images. Here, we demonstrate a fast Fourier transform-based… ▽ More The key to detecting neutral hydrogen during the epoch of reionization (EoR) is to separate the cosmological signal from the dominating foreground radiation. We developed direct optimal mapping (DOM) to map interferometric visibilities; it contains only linear operations, with full knowledge of point spread functions from visibilities to images. Here, we demonstrate a fast Fourier transform-based image power spectrum and its window functions computed from the DOM images. We use noiseless simulation, based on the Hydrogen Epoch of Reionization Array Phase I configuration, to study the image power spectrum properties. The window functions show $<10^{-11}$ of the integrated power leaks from the foreground-dominated region into the EoR window; the 2D and 1D power spectra also verify the separation between the foregrounds and the EoR. △ Less

Submitted 5 July, 2024; v1 submitted 17 November, 2023; originally announced November 2023.

Comments: Published in ApJ

The relationship between IGM Lyman-alpha opacity and galaxy density near the end of reionization

Authors: Holly M. Christenson, George D. Becker, Anson D'Aloisio, Frederick B. Davies, Yongda Zhu, Elisa Boera, Fahad Nasir, Steven R. Furlanetto, Matthew A. Malkan

Abstract: Observed scatter in the Lyman-alpha opacity of quasar sightlines at $z<6$ has motivated measurements of the correlation between Ly$α$ opacity and galaxy density, as models that predict this scatter make strong and sometimes opposite predictions for how they should be related. Our previous work associated two highly opaque Ly$α$ troughs at $z\sim5.7$ with a deficit of Lyman-$α$ emitting galaxies (L… ▽ More Observed scatter in the Lyman-alpha opacity of quasar sightlines at $z<6$ has motivated measurements of the correlation between Ly$α$ opacity and galaxy density, as models that predict this scatter make strong and sometimes opposite predictions for how they should be related. Our previous work associated two highly opaque Ly$α$ troughs at $z\sim5.7$ with a deficit of Lyman-$α$ emitting galaxies (LAEs). In this work, we survey two of the most highly transmissive lines of sight at this redshift, towards the $z=6.02$ quasar SDSS J1306+0356 and the $z=6.17$ quasar PSO J359-06. We find that both fields are underdense in LAEs within 10 $h^{-1}$ Mpc of the quasar sightline, somewhat less extensive than underdensities associated with Ly$α$ troughs. We combine our observations with three additional fields from the literature, and find that while fields with extreme opacities are generally underdense, moderate opacities span a wider density range. The results at high opacities are consistent with models that invoke UV background fluctuations and/or late reionization to explain the observed scatter in IGM Ly$α$ opacities. There is tension at low opacities, however, as the models tend to associate lower IGM Ly$α$ opacities with higher densities. Although the number of fields surveyed is still small, the low-opacity results may support a scenario in which the ionizing background in low-density regions increases more rapidly than some models suggest after becoming ionized. Elevated gas temperatures from recent reionization may also be making these regions more transparent. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Comments: 20 pages, 17 figures. Accepted to ApJ

A self-consistent semi-analytic model for Population III star formation in minihalos

Authors: Sahil Hegde, Steven R. Furlanetto

Abstract: The formation of the first stars marks a watershed moment in the history of our universe. As the first luminous structures, these stars (also known as Population III, or Pop III stars) seed the first galaxies and begin the process of reionization. We construct an analytic model to self-consistently trace the formation of Pop III stars inside minihalos in the presence of the fluctuating ultraviolet… ▽ More The formation of the first stars marks a watershed moment in the history of our universe. As the first luminous structures, these stars (also known as Population III, or Pop III stars) seed the first galaxies and begin the process of reionization. We construct an analytic model to self-consistently trace the formation of Pop III stars inside minihalos in the presence of the fluctuating ultraviolet background, relic dark matter-baryon relative velocities from the early universe, and an X-ray background, which largely work to suppress cooling of gas and delay the formation of this first generation of stars. We demonstrate the utility of this framework in a semi-analytic model for early star formation that also follows the transition between Pop III and Pop II star formation inside these halos. Using our new prescription for the criteria allowing Pop III star formation, we follow a population of dark matter halos from $z=50$ through $z=6$ and examine the global star formation history, finding that each process defines its own key epoch: (i) the stream velocity dominates at the highest redshifts ($z\gtrsim30$), (ii) the UV background sets the tone at intermediate times ($30\gtrsim z\gtrsim15$), and (iii) X-rays control the end of Pop III star formation at the latest times ($z\lesssim 15$). In all of our models, Pop III stars continue to form down to $z\sim 7-10$, when their supernovae will be potentially observable with forthcoming instruments. Finally, we identify the signatures of variations in the Pop III physics in the global 21-cm spin-flip signal of atomic hydrogen. △ Less

Submitted 4 September, 2023; v1 submitted 6 April, 2023; originally announced April 2023.

Comments: 19 pages, 14 figures, published in MNRAS

Search for the Epoch of Reionisation with HERA: Upper Limits on the Closure Phase Delay Power Spectrum

Authors: Pascal M. Keller, Bojan Nikolic, Nithyanandan Thyagarajan, Chris L. Carilli, Gianni Bernardi, Ntsikelelo Charles, Landman Bester, Oleg M. Smirnov, Nicholas S. Kern, Joshua S. Dillon, Bryna J. Hazelton, Miguel F. Morales, Daniel C. Jacobs, Aaron R. Parsons, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley , et al. (58 additional authors not shown)

Abstract: Radio interferometers aiming to measure the power spectrum of the redshifted 21 cm line during the Epoch of Reionisation (EoR) need to achieve an unprecedented dynamic range to separate the weak signal from overwhelming foreground emissions. Calibration inaccuracies can compromise the sensitivity of these measurements to the effect that a detection of the EoR is precluded. An alternative to standa… ▽ More Radio interferometers aiming to measure the power spectrum of the redshifted 21 cm line during the Epoch of Reionisation (EoR) need to achieve an unprecedented dynamic range to separate the weak signal from overwhelming foreground emissions. Calibration inaccuracies can compromise the sensitivity of these measurements to the effect that a detection of the EoR is precluded. An alternative to standard analysis techniques makes use of the closure phase, which allows one to bypass antenna-based direction-independent calibration. Similarly to standard approaches, we use a delay spectrum technique to search for the EoR signal. Using 94 nights of data observed with Phase I of the Hydrogen Epoch of Reionization Array (HERA), we place approximate constraints on the 21 cm power spectrum at $z=7.7$. We find at 95% confidence that the 21 cm EoR brightness temperature is $\le$(372)$^2$ "pseudo" mK$^2$ at 1.14 "pseudo" $h$ Mpc$^{-1}$, where the "pseudo" emphasises that these limits are to be interpreted as approximations to the actual distance scales and brightness temperatures. Using a fiducial EoR model, we demonstrate the feasibility of detecting the EoR with the full array. Compared to standard methods, the closure phase processing is relatively simple, thereby providing an important independent check on results derived using visibility intensities, or related. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 16 pages, 14 figures, accepted for publication by MNRAS

Characterization Of Inpaint Residuals In Interferometric Measurements of the Epoch Of Reionization

Authors: Michael Pagano, Jing Liu, Adrian Liu, Nicholas S. Kern, Aaron Ewall-Wice, Philip Bull, Robert Pascua, Siamak Ravanbakhsh, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Steven Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer , et al. (53 additional authors not shown)

Abstract: 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 du… ▽ More 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. △ Less

Submitted 20 February, 2023; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: 21 pages, 13 figures

Lyman-alpha Emitters in Ionized Bubbles: Constraining the Environment and Ionized Fraction

Authors: A. C. Trapp, Steven R. Furlanetto, Frederick B. Davies

Abstract: Lyman-alpha emitters (LAEs) are excellent probes of the reionization process, as they must be surrounded by large ionized bubbles in order to be visible during the reionization era. Large ionized regions are thought to correspond to over-dense regions and may be protoclusters, making them interesting test-beds for early massive structures. Close associations containing several LAEs are often assum… ▽ More Lyman-alpha emitters (LAEs) are excellent probes of the reionization process, as they must be surrounded by large ionized bubbles in order to be visible during the reionization era. Large ionized regions are thought to correspond to over-dense regions and may be protoclusters, making them interesting test-beds for early massive structures. Close associations containing several LAEs are often assumed to mark over-dense, ionized bubbles. Here, we develop the first framework to quantify the ionization and density fields of high-z galaxy associations. We explore the interplay between (i) the large-scale density of a survey field, (ii) Poisson noise due to the small number density of bright sources at high redshifts (z~7), and (iii) the effects of the ionized fraction on the observation of LAEs. We use Bayesian statistics, a simple model of reionization, and a Monte-Carlo simulation to construct a more comprehensive method for calculating the large-scale density of LAE regions than previous works. We find that Poisson noise has a strong effect on the inferred density of a region and show how the ionized fraction can be inferred. We then apply our framework to the strongest association yet identified: Hu et al. (2021) found 14 LAEs in a volume of ~50,000 cMpc^3 inside the COSMOS field at z~7. We show that this is most likely a 2.5-sigma over-density inside of an ionized or nearly ionized bubble. We also show that this LAE association implies that the global ionized fraction is Q = 0.60 (+0.08,-0.09), within the context of a simple reionization model. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: 14 pages, 14 figures, submitted to MNRAS

Improved Constraints on the 21 cm EoR Power Spectrum and the X-Ray Heating of the IGM with HERA Phase I Observations

Authors: The HERA Collaboration, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Rennan Barkana, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Daniela Breitman, Philip Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon , et al. (70 additional authors not shown)

Abstract: We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $Δ^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that… ▽ More We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $Δ^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that $Δ^2 (k = 0.36$ $h$ Mpc$^{-1}) \leq 3,496$ mK$^2$ at $z = 10.4$, an improvement by a factor of 2.1 and 2.6 respectively. These limits are mostly consistent with thermal noise over a wide range of $k$ after our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration (2022b), we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early as $z = 10.4$, ruling out a broad set of so-called "cold reionization" scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result's 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars. △ Less

Submitted 19 January, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 57 pages, 37 figures. Updated to match the accepted ApJ version. Corresponding author: Joshua S. Dillon

Journal ref: 2023 ApJ 945 124

Impact of instrument and data characteristics in the interferometric reconstruction of the 21 cm power spectrum

Authors: Adélie Gorce, Samskruthi Ganjam, Adrian Liu, Steven G. Murray, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steven Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon , et al. (53 additional authors not shown)

Abstract: Combining the visibilities measured by an interferometer to form a cosmological power spectrum is a complicated process. In a delay-based analysis, the mapping between instrumental and cosmological space is not a one-to-one relation. Instead, neighbouring modes contribute to the power measured at one point, with their respective contributions encoded in the window functions. To better understand t… ▽ More Combining the visibilities measured by an interferometer to form a cosmological power spectrum is a complicated process. In a delay-based analysis, the mapping between instrumental and cosmological space is not a one-to-one relation. Instead, neighbouring modes contribute to the power measured at one point, with their respective contributions encoded in the window functions. To better understand the power measured by an interferometer, we assess the impact of instrument characteristics and analysis choices on these window functions. Focusing on the Hydrogen Epoch of Reionization Array (HERA) as a case study, we find that long-baseline observations correspond to enhanced low-k tails of the window functions, which facilitate foreground leakage, whilst an informed choice of bandwidth and frequency taper can reduce said tails. With simple test cases and realistic simulations, we show that, apart from tracing mode mixing, the window functions help accurately reconstruct the power spectrum estimator of simulated visibilities. The window functions depend strongly on the beam chromaticity, and less on its spatial structure - a Gaussian approximation, ignoring side lobes, is sufficient. Finally, we investigate the potential of asymmetric window functions, down-weighting the contribution of low-k power to avoid foreground leakage. The window functions presented here correspond to the latest HERA upper limits for the full Phase I data. They allow an accurate reconstruction of the power spectrum measured by the instrument and will be used in future analyses to confront theoretical models and data directly in cylindrical space. △ Less

Submitted 11 January, 2023; v1 submitted 7 October, 2022; originally announced October 2022.

Comments: 18 pages, 19 figures, accepted for publication in MNRAS

The extent of intergalactic metal enrichment from galactic winds during the Cosmic Dawn

Authors: Natsuko Yamaguchi, Steven R. Furlanetto, A. C. Trapp

Abstract: One of the key processes driving galaxy evolution during the Cosmic Dawn is supernova feedback. This likely helps regulate star formation inside of galaxies, but it can also drive winds that influence the large-scale intergalactic medium. Here, we present a simple semi-analytic model of supernova-driven galactic winds and explore the contributions of different phases of galaxy evolution to cosmic… ▽ More One of the key processes driving galaxy evolution during the Cosmic Dawn is supernova feedback. This likely helps regulate star formation inside of galaxies, but it can also drive winds that influence the large-scale intergalactic medium. Here, we present a simple semi-analytic model of supernova-driven galactic winds and explore the contributions of different phases of galaxy evolution to cosmic metal enrichment in the high-redshift (z > 6) Universe. We show that models calibrated to the observed galaxy luminosity function at z~6-8 have filling factors ~1% at z~6 and ~0.1% at z~12, with different star formation prescriptions providing about an order of magnitude uncertainty. Despite the small fraction of space filled by winds, these scenarios predict an upper limit to the abundance of metal-line absorbers in quasar spectra at z>5 which is comfortably above that currently observed. We also consider enrichment through winds driven by Pop III star formation in minihalos. We find that these can dominate the total filling factor at z>10 and even compete with winds from normal galaxies at z~6, at least in terms of the total enriched volume. But these regions have much lower overall metallicities, because each one is generated by a small burst of star formation. Finally, we show that Compton cooling of these supernova-driven winds at z>6 has only a small effect on the cosmic microwave background. △ Less

Submitted 18 January, 2024; v1 submitted 19 September, 2022; originally announced September 2022.

Comments: 13 pages, 8 figures, accepted by MNRAS

On the expected purity of photometric galaxy surveys targeting the Cosmic Dawn

Authors: Steven R. Furlanetto, Jordan Mirocha

Abstract: Over the last three decades, photometric galaxy selection using the Lyman-break technique has transformed our understanding of the high-z Universe, providing large samples of galaxies at 3 < z < 8 with relatively small contamination. With the advent of the James Webb Space Telescope, the Lyman-break technique has now been extended to z ~ 17. However, the purity of the resulting samples has not bee… ▽ More Over the last three decades, photometric galaxy selection using the Lyman-break technique has transformed our understanding of the high-z Universe, providing large samples of galaxies at 3 < z < 8 with relatively small contamination. With the advent of the James Webb Space Telescope, the Lyman-break technique has now been extended to z ~ 17. However, the purity of the resulting samples has not been tested. Here we use a simple model, built on the robust foundation of the dark matter halo mass function, to show that the expected level of contamination rises dramatically at z > 10, especially for luminous galaxies, placing stringent requirements on the selection process. The most luminous sources at z > 12 are likely at least ten thousand times rarer than potential contaminants, so extensive spectroscopic followup campaigns may be required to identify a small number of target sources. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: 6 pages, 2 figures, submitted to MNRAS

Balancing the efficiency and stochasticity of star formation with dust extinction in z > 10 galaxies observed by JWST

Authors: Jordan Mirocha, Steven R. Furlanetto

Abstract: Early observations with JWST indicate an over-abundance of bright galaxies at redshifts $z \gtrsim 10$ relative to Hubble-calibrated model predictions. More puzzling still is the apparent lack of evolution in the abundance of such objects between $z \sim 9$ and the highest redshifts yet probed, $z \sim 13$-$17$. In this study, we first show that, despite a poor match with JWST LFs, semi-empirical… ▽ More Early observations with JWST indicate an over-abundance of bright galaxies at redshifts $z \gtrsim 10$ relative to Hubble-calibrated model predictions. More puzzling still is the apparent lack of evolution in the abundance of such objects between $z \sim 9$ and the highest redshifts yet probed, $z \sim 13$-$17$. In this study, we first show that, despite a poor match with JWST LFs, semi-empirical models calibrated to UVLFs and colours at $4 \lesssim z \lesssim 8$ are largely consistent with constraints on the properties of individual JWST galaxies, including their stellar masses, ages, and rest-ultraviolet spectral slopes. We then show that order-of-magnitude scatter in the star formation rate of galaxies (at fixed halo mass) can indeed boost the abundance of bright galaxies, provided that star formation is more efficient than expected in low-mass halos. However, this solution to the abundance problem introduces tension elsewhere: because it relies on the up-scattering of low-mass halos into bright magnitude bins, one expects typical ages, masses, and spectral slopes to be much lower than constraints from galaxies observed thus far. This tension can be alleviated by non-negligible reddening, suggesting that -- if the first batch of photometrically-selected candidates are confirmed -- star formation \textit{and} dust production could be more efficient than expected in galaxies at $z \gtrsim 10$. △ Less

Submitted 19 December, 2022; v1 submitted 26 August, 2022; originally announced August 2022.

Comments: 12 pages, 5 figures, MNRAS accepted

A joint measurement of galaxy luminosity functions and large-scale field densities during the Epoch of Reionization

Authors: A. C. Trapp, Steven R. Furlanetto

Abstract: One of the most exciting advances of the current generation of telescopes has been the detection of galaxies during the epoch of reionization, using deep fields that have pushed these instruments to their limits. It is essential to optimize our analyses of these fields in order to extract as much information as possible from them. In particular, standard methods of measuring the galaxy luminosity… ▽ More One of the most exciting advances of the current generation of telescopes has been the detection of galaxies during the epoch of reionization, using deep fields that have pushed these instruments to their limits. It is essential to optimize our analyses of these fields in order to extract as much information as possible from them. In particular, standard methods of measuring the galaxy luminosity function discard information on large-scale dark matter density fluctuations, even though this large-scale structure drives galaxy formation and reionization during the Cosmic Dawn. Measuring these densities would provide a bedrock observable, connecting galaxy surveys to theoretical models of the reionization process and structure formation. Here, we use existing Hubble deep field data to simultaneously fit the universal luminosity function and measure large-scale densities for each Hubble deep field at $z =$ 6--8 by directly incorporating priors on the large-scale density field and galaxy bias. Our fit of the universal luminosity function is consistent with previous methods but differs in the details. For the first time, we measure the underlying densities of the survey fields, including the most over/under-dense Hubble fields. We show that the distribution of densities is consistent with current predictions for cosmic variance. This analysis on just 17 fields is a small sample of what will be possible with the James Webb Space Telescope, which will measure hundreds of fields at comparable (or better) depths and at higher redshifts. △ Less

Submitted 11 August, 2022; originally announced August 2022.

LIMFAST. II. Line Intensity Mapping as a Probe of High-Redshift Galaxy Formation

Authors: Guochao Sun, Lluís Mas-Ribas, Tzu-Ching Chang, Steven R. Furlanetto, Richard H. Mebane, Michael O. Gonzalez, Jasmine Parsons, A. C. Trapp

Abstract: The epoch of reionization (EoR) offers a unique window into the dawn of galaxy formation, through which high-redshift galaxies can be studied by observations of both themselves and their impact on the intergalactic medium. Line intensity mapping (LIM) promises to explore cosmic reionization and its driving sources by measuring intensity fluctuations of emission lines tracing the cosmic gas in vary… ▽ More The epoch of reionization (EoR) offers a unique window into the dawn of galaxy formation, through which high-redshift galaxies can be studied by observations of both themselves and their impact on the intergalactic medium. Line intensity mapping (LIM) promises to explore cosmic reionization and its driving sources by measuring intensity fluctuations of emission lines tracing the cosmic gas in varying phases. Using LIMFAST, a novel semi-numerical tool designed to self-consistently simulate LIM signals of multiple EoR probes, we investigate how building blocks of galaxy formation and evolution theory, such as feedback-regulated star formation and chemical enrichment, might be studied with multi-tracer LIM during the EoR. On galaxy scales, we show that the star formation law and the feedback associated with star formation can be indicated by both the shape and redshift evolution of LIM power spectra. For a baseline model of metal production that traces star formation, we find that lines highly sensitive to metallicity are generally better probes of galaxy formation models. On larger scales, we demonstrate that inferring ionized bubble sizes from cross-correlations between tracers of ionized and neutral gas requires a detailed understanding of the astrophysics that shape the line luminosity--halo mass relation. Despite various modeling and observational challenges, wide-area, multi-tracer LIM surveys will provide important high-redshift tests for the fundamentals of galaxy formation theory, especially the interplay between star formation and feedback by accessing statistically the entire low-mass population of galaxies as ideal laboratories, complementary to upcoming surveys of individual sources by new-generation telescopes. △ Less

Submitted 23 March, 2023; v1 submitted 28 June, 2022; originally announced June 2022.

Comments: 24 pages, 15 figures; to appear in ApJ

Direct Optimal Mapping for 21cm Cosmology: A Demonstration with the Hydrogen Epoch of Reionization Array

Authors: Zhilei Xu, Jacqueline N. Hewitt, Kai-Feng Chen, Honggeun Kim, Joshua S. Dillon, Nicholas S. Kern, Miguel F. Morales, Bryna J. Hazelton, Ruby Byrne, Nicolas Fagnoni, Eloy de Lera Acedo, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba , et al. (56 additional authors not shown)

Abstract: Motivated by the desire for wide-field images with well-defined statistical properties for 21cm cosmology, we implement an optimal mapping pipeline that computes a maximum likelihood estimator for the sky using the interferometric measurement equation. We demonstrate this direct optimal mapping with data from the Hydrogen Epoch of Reionization (HERA) Phase I observations. After validating the pipe… ▽ More Motivated by the desire for wide-field images with well-defined statistical properties for 21cm cosmology, we implement an optimal mapping pipeline that computes a maximum likelihood estimator for the sky using the interferometric measurement equation. We demonstrate this direct optimal mapping with data from the Hydrogen Epoch of Reionization (HERA) Phase I observations. After validating the pipeline with simulated data, we develop a maximum likelihood figure-of-merit for comparing four sky models at 166MHz with a bandwidth of 100kHz. The HERA data agree with the GLEAM catalogs to <10%. After subtracting the GLEAM point sources, the HERA data discriminate between the different continuum sky models, providing most support for the model of Byrne et al. 2021. We report the computation cost for mapping the HERA Phase I data and project the computation for the HERA 320-antenna data; both are feasible with a modern server. The algorithm is broadly applicable to other interferometers and is valid for wide-field and non-coplanar arrays. △ Less

Submitted 26 October, 2022; v1 submitted 12 April, 2022; originally announced April 2022.

Comments: 16 pages, 10 figures, 2 tables, published on ApJ

A galaxy-free phenomenological model for the 21-cm power spectrum during reionization

Authors: Jordan Mirocha, Julian B. Muñoz, Steven R. Furlanetto, Adrian Liu, Andrei Mesinger

Abstract: Upper limits from the current generation of interferometers targeting the 21-cm signal from high redshifts have recently begun to rule out physically realistic, though still extreme, models of the Epoch of Reionization (EoR). While inferring the detailed properties of the first galaxies is one of the most important motivations for measuring the high-$z$ 21-cm signal, they can also provide useful c… ▽ More Upper limits from the current generation of interferometers targeting the 21-cm signal from high redshifts have recently begun to rule out physically realistic, though still extreme, models of the Epoch of Reionization (EoR). While inferring the detailed properties of the first galaxies is one of the most important motivations for measuring the high-$z$ 21-cm signal, they can also provide useful constraints on the properties of the intergalactic medium (IGM). Motivated by this, we build a simple, phenomenological model for the 21-cm power spectrum that works directly in terms of IGM properties, which bypasses the computationally expensive 3-D semi-numerical modeling generally employed in inference pipelines and avoids explicit assumptions about galaxy properties. The key simplifying assumptions are that (i) the ionization field is binary, and composed of spherical bubbles with an abundance described well by a parametric bubble size distribution, and (ii) that the spin temperature of the ``bulk'' IGM outside bubbles is uniform. Despite the simplicity of the model, the mean ionized fraction and spin temperature of the IGM recovered from mock 21-cm power spectra generated with \textsc{21cmfast} are generally in good agreement with the true input values. This suggests that it is possible to obtain comparable constraints on the IGM using models with very different assumptions, parameters, and priors. Our approach will thus be complementary to semi-numerical models as upper limits continue to improve in the coming years. △ Less

Submitted 27 May, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 23 pages, 14 figures, accepted for publication in MNRAS, code publicly available at https://github.com/mirochaj/micro21cm

Deep Realistic Extragalactic Model (DREaM) Galaxy Catalogs: Predictions for a Roman Ultra-Deep Field

Authors: Nicole E. Drakos, Bruno Villasenor, Brant E. Robertson, Ryan Hausen, Mark E. Dickinson, Henry C. Ferguson, Steven R. Furlanetto, Jenny E. Greene, Piero Madau, Alice E. Shapley, Daniel P. Stark, Risa H. Wechsler

Abstract: 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… ▽ More 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}<-17$ galaxies, with more than $10^4$ at redshifts $z>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. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: 39 pages, 27 figures. Submitted to ApJ. Comments welcome

A Framework for Simultaneously Measuring Field Densities and the High-z Luminosity Function

Authors: A. C. Trapp, Steven R. Furlanetto, Jinghong Yang

Abstract: Cosmic variance from large-scale structure will be a major source of uncertainty for galaxy surveys at z > 6, but that same structure will also provide an opportunity to identify and study dense environments in the early Universe. Using a robust model for galaxy clustering, we directly incorporate large-scale densities into an inference framework that simultaneously measures the high-z (z > 6) UV… ▽ More Cosmic variance from large-scale structure will be a major source of uncertainty for galaxy surveys at z > 6, but that same structure will also provide an opportunity to identify and study dense environments in the early Universe. Using a robust model for galaxy clustering, we directly incorporate large-scale densities into an inference framework that simultaneously measures the high-z (z > 6) UV luminosity function and the average matter density of each distinct volume in a survey. Through this framework, we forecast the performance of several major upcoming James Webb Space Telescope (JWST) galaxy surveys. We find that they can constrain field matter densities down to the theoretical limit imposed by Poisson noise and unambiguously identify over-dense (and under-dense) regions on transverse scales of tens of comoving Mpc. We also predict JWST will measure the luminosity function with a precision at z = 12 comparable to existing Hubble Space Telescope's constraints at z = 8 (and even better for the faint-end slope). We also find that wide-field surveys are especially important in distinguishing luminosity function models. △ Less

Submitted 11 October, 2021; originally announced October 2021.

Constraints on the End of Reionization from the Density Fields Surrounding Two Highly Opaque Quasar Sightlines

Authors: Holly M. Christenson, George D. Becker, Steven R. Furlanetto, Frederick B. Davies, Matthew A. Malkan, Yongda Zhu, Elisa Boera, Adam Trapp

Abstract: The observed large-scale scatter in Lyman $α$ opacity of the intergalactic medium at $z<6$ implies large fluctuations in the neutral hydrogen fraction that are unexpected long after reionization has ended. A number of models have emerged to explain these fluctuations that make testable predictions for the relationship between Ly$α$ opacity and density. We present selections of $z=5.7$ Ly$α$-emitti… ▽ More The observed large-scale scatter in Lyman $α$ opacity of the intergalactic medium at $z<6$ implies large fluctuations in the neutral hydrogen fraction that are unexpected long after reionization has ended. A number of models have emerged to explain these fluctuations that make testable predictions for the relationship between Ly$α$ opacity and density. We present selections of $z=5.7$ Ly$α$-emitting galaxies (LAEs) in the fields surrounding two highly opaque quasar sightlines with long Ly$α$ troughs. The fields lie towards the $z=6.0$ quasar ULAS J0148+0600, for which we re-analyze previously published results using improved photometric selection, and towards the $z=6.15$ quasar SDSS J1250+3130, for which results are presented here for the first time. In both fields, we report a deficit of LAEs within 20 $h^{-1}$ Mpc of the quasar. The association of highly opaque sightlines with galaxy underdensities in these two fields is consistent with models in which the scatter in Ly$α$ opacity is driven by large-scale fluctuations in the ionizing UV background, or by an ultra-late reionization that has not yet concluded at $z=5.7$. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: Accepted to ApJ. 17 pages, 13 figures

Automated Detection of Antenna Malfunctions in Large-N Interferometers: A Case Study with the Hydrogen Epoch of Reionization Array

Authors: Dara Storer, Joshua S. Dillon, Daniel C. Jacobs, Miguel F. Morales, Bryna J. Hazelton, Aaron Ewall-Wice, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steven Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Scott Dynes , et al. (53 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 4 May, 2022; v1 submitted 26 September, 2021; originally announced September 2021.

Comments: 31 pages, 17 figures

Journal ref: Radio Science, vol. 57, no. 1, 2022

Bursty star formation during the Cosmic Dawn driven by delayed stellar feedback

Authors: Steven R. Furlanetto, Jordan Mirocha

Abstract: In recent years, several analytic models have demonstrated that simple assumptions about halo growth and feedback-regulated star formation can match the (limited) existing observational data on galaxies at z>6. By extending such models, we demonstrate that imposing a time delay on stellar feedback (as inevitably occurs in the case of supernova explosions) induces burstiness in small galaxies. Alth… ▽ More In recent years, several analytic models have demonstrated that simple assumptions about halo growth and feedback-regulated star formation can match the (limited) existing observational data on galaxies at z>6. By extending such models, we demonstrate that imposing a time delay on stellar feedback (as inevitably occurs in the case of supernova explosions) induces burstiness in small galaxies. Although supernova progenitors have short lifetimes (~5-30 Myr), the delay exceeds the dynamical time of galaxies at such high redshifts. As a result, star formation proceeds unimpeded by feedback for several cycles and "overshoots" the expectations of feedback-regulated star formation models. We show that such overshoot is expected even in atomic cooling halos, with masses up to ~10^10.5 Msun at z>6. However, these burst cycles damp out quickly in massive galaxies, because large haloes are more resistant to feedback so retain a continuous gas supply. Bursts in small galaxies - largely beyond the reach of existing observations - induce a scatter in the luminosity of these haloes (of ~1 mag) and increase the time-averaged star formation efficiency by up to an order of magnitude. This kind of burstiness can have substantial effects on the earliest phases of star formation and reionization. △ Less

Submitted 9 September, 2021; originally announced September 2021.

Comments: 15 pages, 10 figures, submitted to MNRAS

First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum

Authors: The HERA Collaboration, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, Taylor Dibblee-Barkman, Joshua S. Dillon, John Ely, Aaron Ewall-Wice, Nicolas Fagnoni, Randall Fritz , et al. (52 additional authors not shown)

Abstract: We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground… ▽ More We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of $Δ^2_{21} \le (30.76)^2\ {\rm mK}^2$ at $k=0.192\ h\ {\rm Mpc}^{-1}$ at $z=7.9$, and also $Δ^2_{21} \le (95.74)^2\ {\rm mK}^2$ at $k=0.256\ h\ {\rm Mpc}^{-1}$ at $z=10.4$. At $z=7.9$, these limits are the most sensitive to-date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier $k_\parallel$ modes, at high $k_\parallel$ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable sub-reflections, or residual instrumental cross-coupling, and warrant further study. This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms of loss or bias associated with the pipeline. Overall, these results are a promising first step in the development of a tuned, instrument-specific analysis pipeline for HERA, particularly as Phase II construction is completed en route to reaching the full sensitivity of the experiment. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: Accepted to ApJ. https://reionization.org/science/public-data-release-1/

Revealing the formation histories of the first stars with the cosmic near-infrared background

Authors: Guochao Sun, Jordan Mirocha, Richard H. Mebane, Steven R. Furlanetto

Abstract: 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 form… ▽ More 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. △ Less

Submitted 15 September, 2021; v1 submitted 20 July, 2021; originally announced July 2021.

Comments: 19 pages, 10 figures, accepted for publication in MNRAS

The Detection of Ionized Carbon Emission at z~8

Authors: Michael W. Topping, Alice E. Shapley, Daniel P. Stark, Ryan Endsley, Brant Robertson, Jenny E. Greene, Steven R. Furlanetto, Mengtao Tang

Abstract: We present deep Keck/MOSFIRE $H$-band spectroscopic observations covering the [CIII],CIII]$λ\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]$λ1907$ for the galaxy AEGIS-33376, finding $z_{\rm spec}=7.945\pm0.001$. We measure a [CIII]$λ$1907 flux of… ▽ More We present deep Keck/MOSFIRE $H$-band spectroscopic observations covering the [CIII],CIII]$λ\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]$λ1907$ for the galaxy AEGIS-33376, finding $z_{\rm spec}=7.945\pm0.001$. We measure a [CIII]$λ$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]$λ1909$ based on strong sky-line contamination, we assume a [CIII]$λ$1907/CIII]$λ1909$ doublet ratio of 1.5 and infer a total [CIII],CIII]$λ\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]$λ1907$ emission is present for AEGIS-33376 throughout our observations. We also confirm that the significance of the [CIII]$λ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]$λ\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]$λ\lambda1907,1909$ detections at $z>7$, while $\textit{JWST}$ will provide the first statistical samples of such measurements at these redshifts. △ Less

Submitted 4 August, 2021; v1 submitted 13 July, 2021; originally announced July 2021.

Comments: 6 pages, 3 figures, 1 table; Accepted to ApJL

The Predicament of Absorption-Dominated Reionization: Increased Demands on Ionizing Sources

Authors: Frederick B. Davies, Sarah E. I. Bosman, Steven R. Furlanetto, George D. Becker, Anson D'Aloisio

Abstract: 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 residua… ▽ More 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$α$ 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}ξ_{\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. △ Less

Submitted 21 May, 2021; originally announced May 2021.

Comments: 8 pages, 5 figures, submitted to ApJL

Effects of model incompleteness on the drift-scan calibration of radio telescopes

Authors: Bharat K. Gehlot, Daniel C. Jacobs, Judd D. Bowman, Nivedita Mahesh, Steven G. Murray, Matthew Kolopanis, Adam P. Beardsley, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Gianni Bernardi, Tashalee S. Billings, Richard F. Bradley, Phil Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, Joshua S. Dillon, John Ely , et al. (54 additional authors not shown)

Abstract: 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 at… ▽ More 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. △ Less

Submitted 15 July, 2021; v1 submitted 25 April, 2021; originally announced April 2021.

Comments: 16 pages, 13 figures, 1 table; accepted for publication in MNRAS main journal

Validation of the HERA Phase I Epoch of Reionization 21 cm Power Spectrum Software Pipeline

Authors: James E. Aguirre, Steven G. Murray, Robert Pascua, Zachary E. Martinot, Jacob Burba, Joshua S. Dillon, Daniel C. Jacobs, Nicholas S. Kern, Piyanat Kittiwisit, Matthew Kolopanis, Adam Lanman, Adrian Liu, Lily Whitler, Zara Abdurashidova, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Steve Carey, Chris L. Carilli , et al. (51 additional authors not shown)

Abstract: We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the… ▽ More We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: 32 pages, 20 figures. Submitted to the Astrophysical Journal

A Real Time Processing System for Big Data in Astronomy: Applications to HERA

Authors: Paul La Plante, Peter K. G. Williams, Matthew Kolopanis, Joshua S. Dillon, Adam P. Beardsley, Nicholas S. Kern, Michael Wilensky, Zaki S. Ali, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Yanga Balfour, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Phil Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, John Ely , et al. (50 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 30 September, 2021; v1 submitted 8 April, 2021; originally announced April 2021.

Comments: 15 pages, 3 figures, published in Astronomy and Computing

Journal ref: Astronomy and Computing 36 (2021) 100489

Methods of Error Estimation for Delay Power Spectra in $21\,\textrm{cm}$ Cosmology

Authors: Jianrong Tan, Adrian Liu, Nicholas S. Kern, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steven Carey, Christopher L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, Joshua S. Dillon, John Ely, Aaron Ewall-Wice, Nicolas Fagnoni , et al. (49 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 25 May, 2021; v1 submitted 17 March, 2021; originally announced March 2021.

Comments: 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

Improved treatments of the ionizing photon mean free path in semi-numerical simulations of reionization

Authors: Frederick B. Davies, Steven R. Furlanetto

Abstract: Efficient and accurate simulations of the reionization epoch are crucial to exploring the vast uncharted parameter space that will soon be constrained by measurements of the 21 cm power spectrum. One of these parameters, $R_{\rm max}$, is meant to characterize the absorption of photons by residual neutral gas inside of ionized regions, but has historically been implemented in a very simplistic fas… ▽ More Efficient and accurate simulations of the reionization epoch are crucial to exploring the vast uncharted parameter space that will soon be constrained by measurements of the 21 cm power spectrum. One of these parameters, $R_{\rm max}$, is meant to characterize the absorption of photons by residual neutral gas inside of ionized regions, but has historically been implemented in a very simplistic fashion acting only as a maximum filtering scale. We leverage the correspondence between excursion set methods and the integrated flux from ionizing sources to define two physically-motivated prescriptions of the mean free path of ionizing photons that smoothly attenuate the contribution from distant sources. Implementation of these methods in semi-numerical reionization codes requires only modest additional computational effort due to the fact that spatial filtering is still performed on scales larger than the characteristic absorption distance. We find that our smoothly-defined mean free path prescriptions more effectively suppress large-scale structures in the ionization field in semi-numerical reionization simulations compared to the standard $R_{\rm max}$ approach, and the magnitude of the mean free path modulates the power spectrum in a much smoother manner. We show that this suppression of large-scale power is significant enough to be relevant for upcoming 21 cm power spectrum observations. Finally, we show that in our model the mean free path plays a larger role in regulating the reionization history than in models using $R_{\rm max}$. △ Less

Submitted 7 April, 2022; v1 submitted 17 March, 2021; originally announced March 2021.

Comments: 13 pages, 12+2 figures, accepted for publication in MNRAS

Quasi-equilibrium models of high-redshift disc galaxy evolution

Authors: Steven R. Furlanetto

Abstract: 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 sophisti… ▽ More 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. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: 20 pages, 13 figures, in press at MNRAS

A Flexible Analytic Model of Cosmic Variance in the First Billion Years

Authors: A. C. Trapp, Steven R. Furlanetto

Abstract: 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… ▽ More 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. △ Less

Submitted 20 October, 2020; v1 submitted 10 September, 2020; originally announced September 2020.

Comments: 16 pages, 15 figures, accepted by MNRAS

Measuring HERA's primary beam in-situ: methodology and first results

Authors: Chuneeta D. Nunhokee, Aaron R. Parsons, Nicholas S. Kern, Bojan Nikolic, Jonathan C. Pober, Gianni Bernardi, Chris L. Carilli, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de~Lera~Acedo, Joshua S. Dillon, Aaron Ewall-Wice, Nicolas Fagnoni, Randall Fritz , et al. (42 additional authors not shown)

Abstract: 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 partic… ▽ More 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. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 22 pages, 22 figures, accepted for publication in ApJ

Detection of Cosmic Structures using the Bispectrum Phase. II. First Results from Application to Cosmic Reionization Using the Hydrogen Epoch of Reionization Array

Authors: Nithyanandan Thyagarajan, Chris L. Carilli, Bojan Nikolic, James Kent, Andrei Mesinger, Nicholas S. Kern, Gianni Bernardi, Siyanda Matika, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Steve Carey, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, Joshua S. Dillon, John Ely , et al. (47 additional authors not shown)

Abstract: Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirem… ▽ More Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1σ$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $κ_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $κ_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues. △ Less

Submitted 2 July, 2020; v1 submitted 20 May, 2020; originally announced May 2020.

Comments: 22 pages, 12 figures (including sub-figures). Published in PhRvD. Abstract may be slightly abridged compared to the actual manuscript due to length limitations on arXiv

Journal ref: Phys. Rev. D 102, 022002 (2020)

Foreground modelling via Gaussian process regression: an application to HERA data

Authors: Abhik Ghosh, Florent Mertens, Gianni Bernardi, Mário G. Santos, Nicholas S. Kern, Christopher L. Carilli, Trienko L. Grobler, Léon V. E. Koopmans, Daniel C. Jacobs, Adrian Liu, Aaron R. Parsons, Miguel F. Morales, James E. Aguirre, Joshua S. Dillon, Bryna J. Hazelton, Oleg M. Smirnov, Bharat K. Gehlot, Siyanda Matika, Paul Alexander, Zaki S. Ali, Adam P. Beardsley, Roshan K. Benefo, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley , et al. (48 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 12 May, 2020; v1 submitted 13 April, 2020; originally announced April 2020.

Comments: 15 pages, 15 figures, 1 table, Accepted to MNRAS

Redundant-Baseline Calibration of the Hydrogen Epoch of Reionization Array

Authors: Joshua S. Dillon, Max Lee, Zaki S. Ali, Aaron R. Parsons, Naomi Orosz, Chuneeta Devi Nunhokee, Paul La Plante, Adam P. Beardsley, Nicholas S. Kern, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Yanga Balfour, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Phil Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo , et al. (54 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 3 November, 2020; v1 submitted 18 March, 2020; originally announced March 2020.

Comments: 24 Pages, 19 Figures. Updated to match the accepted MNRAS version

Absolute Calibration Strategies for the Hydrogen Epoch of Reionization Array and Their Impact on the 21 cm Power Spectrum

Authors: Nicholas S. Kern, Joshua S. Dillon, Aaron R. Parsons, Christopher L. Carilli, Gianni Bernardi, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki S. Ali, Yanga Balfour, Adam P. Beardsley, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steven Carey, Carina Cheng, David R. DeBoer, Matt Dexter, Eloy de Lera Acedo, John Ely, Aaron Ewall-Wice, Nicolas Fagnoni, Randall Fritz , et al. (47 additional authors not shown)

Abstract: We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the… ▽ More We discuss absolute calibration strategies for Phase I of the Hydrogen Epoch of Reionization Array (HERA), which aims to measure the cosmological 21 cm signal from the Epoch of Reionization (EoR). HERA is a drift-scan array with a 10 degree wide field of view, meaning bright, well-characterized point source transits are scarce. This, combined with HERA's redundant sampling of the uv plane and the modest angular resolution of the Phase I instrument, make traditional sky-based and self-calibration techniques difficult to implement with high dynamic range. Nonetheless, in this work we demonstrate calibration for HERA using point source catalogues and electromagnetic simulations of its primary beam. We show that unmodeled diffuse flux and instrumental contaminants can corrupt the gain solutions, and present a gain smoothing approach for mitigating their impact on the 21 cm power spectrum. We also demonstrate a hybrid sky and redundant calibration scheme and compare it to pure sky-based calibration, showing only a marginal improvement to the gain solutions at intermediate delay scales. Our work suggests that the HERA Phase I system can be well-calibrated for a foreground-avoidance power spectrum estimator by applying direction-independent gains with a small set of degrees of freedom across the frequency and time axes. △ Less

Submitted 4 January, 2020; v1 submitted 28 October, 2019; originally announced October 2019.

Comments: Accepted to ApJ

The Effects of Population III X-ray and Radio Backgrounds on the Cosmological 21-cm Signal

Authors: Richard H. Mebane, Jordan Mirocha, Steven R. Furlanetto

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… ▽ More 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 $\sim 10^{-4}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ between $z \sim 10$ and $z \sim 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 \, M_\odot$. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 11 pages, 8 figures, submitted to MNRAS

The Fundamentals of the 21-cm Line

Authors: Steven R. Furlanetto

Abstract: We review some of the fundamental physics necessary for computing the highly-redshifted spin-flip background. We first discuss the radiative transfer of the 21-cm line and define the crucial quantities of interest. We then review the processes that set the spin temperature of the transition, with a particular focus on Wouthuysen-Field coupling, which is likely to be the most important process duri… ▽ More We review some of the fundamental physics necessary for computing the highly-redshifted spin-flip background. We first discuss the radiative transfer of the 21-cm line and define the crucial quantities of interest. We then review the processes that set the spin temperature of the transition, with a particular focus on Wouthuysen-Field coupling, which is likely to be the most important process during and after the Cosmic Dawn. Finally, we discuss processes that heat the intergalactic medium during the Cosmic Dawn, including the scattering of Lyman-alpha, cosmic microwave background, and X-ray photons. △ Less

Submitted 26 September, 2019; originally announced September 2019.

Comments: To appear as a book chapter in "The Cosmic 21-cm Revolution: Charting the first billion years of our Universe," Ed Andrei Mesinger (Bristol: IOP Publishing Ltd) AAS-IOP ebooks http://www.iopscience.org/books/aas. arXiv admin note: text overlap with arXiv:astro-ph/0608032

Physical Cosmology From the 21-cm Line

Authors: Steven R. Furlanetto

Abstract: We describe how the high-redshift 21-cm background can be used to improve both our understanding of the fundamental cosmological parameters of our Universe and exotic processes originating in the dark sector. The 21-cm background emerging during the cosmological Dark Ages, the era between hydrogen recombination and the formation of the first luminous sources (likely at z ~ 30), is difficult to mea… ▽ More We describe how the high-redshift 21-cm background can be used to improve both our understanding of the fundamental cosmological parameters of our Universe and exotic processes originating in the dark sector. The 21-cm background emerging during the cosmological Dark Ages, the era between hydrogen recombination and the formation of the first luminous sources (likely at z ~ 30), is difficult to measure but provides several powerful advantages for these purposes: in addition to the lack of astrophysical contamination, it will allow probes of very small scale structure over a very large volume. Additionally, the 21-cm background is sensitive to the thermal state of the intergalactic hydrogen and therefore probes any exotic processes (including, e.g., dark matter scattering or decay and primordial black holes) during that era. After astrophysical sources have formed, cosmological information can be separated from astrophysical effects on the 21-cm background through methods such as redshift space distortions, joint modeling, and by searching for indirect effects on the astrophysical sources themselves. △ Less

Submitted 26 September, 2019; originally announced September 2019.

Comments: To appear as a book chapter in "The Cosmic 21-cm Revolution: Charting the first billion years of our Universe," Ed Andrei Mesinger (Bristol: IOP Publishing Ltd) AAS-IOP ebooks http://www.iopscience.org/books/aas