Showing 1–50 of 102 results for author: Bull, P

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

Constraining Post-Newtonian Parameters with the Cosmic Microwave Background

Authors: Daniel B. Thomas, Theodore Anton, Timothy Clifton, Philip Bull

Abstract: The Parameterised Post-Newtonian (PPN) approach is the default framework for performing precision tests of gravity in nearby astrophysical systems. In recent works we have extended this approach for cosmological applications, and in this paper we use observations of the anisotropies in the Cosmic Microwave Background to constrain the time variation of the PPN parameters $α$ and $γ$ between last sc… ▽ More The Parameterised Post-Newtonian (PPN) approach is the default framework for performing precision tests of gravity in nearby astrophysical systems. In recent works we have extended this approach for cosmological applications, and in this paper we use observations of the anisotropies in the Cosmic Microwave Background to constrain the time variation of the PPN parameters $α$ and $γ$ between last scattering and the present day. We find their time-averages over cosmological history should be within $\sim 20\%$ of their values in GR, with $\barα=0.89^{+0.08}_{-0.09}$ and $\barγ=0.90^{+0.07}_{-0.08}$ at the $68\%$ confidence level. We also constrain the time derivatives of these parameters, and find that their present-day values should be within a factor of two of the best Solar System constraints. Many of these results have no counter-part from Solar System observations, and are entirely new constraints on the gravitational interaction. In all cases, we find that the data strongly prefer $\barα\simeq \barγ$, meaning that observers would typically find local gravitational physics to be compatible with GR, despite considerable variation of $α$ and $γ$ being allowed over cosmic history. This study lays the groundwork for future precision tests of gravity that combine observations made over all cosmological and astrophysical scales of length and time. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: Comments welcome

Radio Frequency Interference from Radio Navigation Satellite Systems: simulations and comparison to MeerKAT single-dish data

Authors: Brandon Engelbrecht, Mario G. Santos, José Fonseca, Yichao Li, Jingying Wang, Melis O. Irfan, Stuart E. Harper, Keith Grainge, Philip Bull, Isabella P. Carucci, Steven Cunnington, Alkistis Pourtsidou, Marta Spinelli, Laura Wolz

Abstract: Radio Frequency Interference (RFI) is emitted from various sources, terrestrial or orbital, and create a nuisance for ground-based 21cm experiments. In particular, single-dish 21cm intensity mapping experiments will be highly susceptible to contamination from these sources due to its wide primary beam and sensitivity. This work aims to simulate the contamination effects emitted from orbital source… ▽ More Radio Frequency Interference (RFI) is emitted from various sources, terrestrial or orbital, and create a nuisance for ground-based 21cm experiments. In particular, single-dish 21cm intensity mapping experiments will be highly susceptible to contamination from these sources due to its wide primary beam and sensitivity. This work aims to simulate the contamination effects emitted from orbital sources in the Radio Navigational Satellite System within the 1100-1350 MHz frequency. This simulation can be split into two parts: (I) satellite positioning, emission power, and beam response on the telescope and (II) fitting of the satellite signal to data in order to improve the original model. We use previously observed single dish MeerKAT L-band data which needs to be specially calibrated to include data contaminated by satellite-based RFI. We find that due to non-linearity effects, it becomes non-trivial to fit the satellite power. However, when masking regions where this non-linearity is problematic, we can recreate the satellite contamination with high accuracy around its peak frequencies. The simulation can predict satellite movements and signal for past and future observations, which can help in RFI avoidance and testing novel cleaning methods. The predicted signal from simulations sits below the noise in the target cosmology window for the L-band (970 - 1015 MHz) making it difficult to confirm any out-of-band emission from satellites. However, a power spectrum analysis shows that such signal can still contaminate the 21cm power spectrum at these frequencies. In our simulations, this contamination overwhelms the auto-power spectrum but still allows for a clean detection of the signal in cross-correlations with mild foreground cleaning. Whether such contamination does exist one will require further characterization of the satellite signals far away from their peak frequencies. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: 20 pages, 27 figures, 4 tables, Comments welcome

Modes of the Dark Ages 21cm field accessible to a lunar radio interferometer

Authors: Philip Bull, Caroline Guandalin, Chris Addis

Abstract: At redshifts beyond $z \gtrsim 30$, the 21cm line from neutral hydrogen is expected to be essentially the only viable probe of the 3D matter distribution. The lunar far-side is an extremely appealing site for future radio arrays that target this signal, as it is protected from terrestrial radio frequency interference, and has no ionosphere to attenuate and absorb radio emission at low frequencies… ▽ More At redshifts beyond $z \gtrsim 30$, the 21cm line from neutral hydrogen is expected to be essentially the only viable probe of the 3D matter distribution. The lunar far-side is an extremely appealing site for future radio arrays that target this signal, as it is protected from terrestrial radio frequency interference, and has no ionosphere to attenuate and absorb radio emission at low frequencies (tens of MHz and below). We forecast the sensitivity of low-frequency lunar radio arrays to the bispectrum of the 21cm brightness temperature field, which can in turn be used to probe primordial non-Gaussianity generated by particular early universe models. We account for the loss of particular regions of Fourier space due to instrumental limitations and systematic effects, and predict the sensitivity of different representative array designs to local-type non-Gaussianity in the bispectrum, parametrised by $f_{\rm NL}$. Under the most optimistic assumption of sample variance-limited observations, we find that $σ(f_{\rm NL}) \lesssim 0.01$ could be achieved for several broad redshift bins at $z \gtrsim 30$ if foregrounds can be removed effectively. These values degrade to between $σ(f_{\rm NL}) \sim 0.03$ and $0.7$ for $z=30$ to $z=170$ respectively when a large foreground wedge region is excluded. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 14 pages, 5 figures, published as part of Phil. Trans. R. Soc. A discussion meeting issue "Astronomy from the Moon: the next decades (part 2)"

Journal ref: Phil.Trans.R.Soc.A 382, 20230072 (2024)

High-dimensional inference of radio interferometer beam patterns I: Parametric model of the HERA beams

Authors: Michael J. Wilensky, Jacob Burba, Philip Bull, Hugh Garsden, Katrine A. Glasscock, Nicolas Fagnoni, Eloy de Lera Acedo, David R. DeBoer, Nima Razavi-Ghods

Abstract: Accurate modelling of the primary beam is an important but difficult task in radio astronomy. For high dynamic range problems such as 21cm intensity mapping, small modelling errors in the sidelobes and spectral structure of the beams can translate into significant systematic errors. Realistic beams exhibit complex spatial and spectral structure, presenting a major challenge for beam measurement an… ▽ More Accurate modelling of the primary beam is an important but difficult task in radio astronomy. For high dynamic range problems such as 21cm intensity mapping, small modelling errors in the sidelobes and spectral structure of the beams can translate into significant systematic errors. Realistic beams exhibit complex spatial and spectral structure, presenting a major challenge for beam measurement and calibration methods. In this paper series, we present a Bayesian framework to infer per-element beam patterns from the interferometric visibilities for large arrays with complex beam structure, assuming a particular (but potentially uncertain) sky model and calibration solution. In this first paper, we develop a compact basis for the beam so that the Bayesian computation is tractable with high-dimensional sampling methods. We use the Hydrogen Epoch of Reionization Array (HERA) as an example, verifying that the basis is capable of describing its single-element E-field beam (i.e. without considering array effects like mutual coupling) with a relatively small number of coefficients. We find that 32 coefficients per feed, incident polarization, and frequency, are sufficient to give percent-level and $\sim$10\% errors in the mainlobe and sidelobes respectively for the current HERA Vivaldi feeds, improving to $\sim 0.1\%$ and $\sim 1\%$ for 128 coefficients. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Disentangling the anisotropic radio sky: Fisher forecasts for 21cm arrays

Authors: Zheng Zhang, Philip Bull, Katrine A. Glasscock

Abstract: The existence of a radio synchrotron background (RSB) excess is implied by a number of measurements, including excess emission seen by the ARCADE~2 and LWA experiments. Highly sensitive wideband radio arrays, of the kind used to measure the cosmic 21cm signal, provide a promising way to further constrain the RSB excess through its anisotropy, providing additional insight into its origin. We presen… ▽ More The existence of a radio synchrotron background (RSB) excess is implied by a number of measurements, including excess emission seen by the ARCADE~2 and LWA experiments. Highly sensitive wideband radio arrays, of the kind used to measure the cosmic 21cm signal, provide a promising way to further constrain the RSB excess through its anisotropy, providing additional insight into its origin. We present a framework for evaluating the potential of 21cm arrays to disentangle different components of the diffuse radio sky based on the combination of their frequency spectrum and angular power spectrum (APS). The formalism is designed to calculate uncertainties due to the intrinsic cosmic variance alone or together with instrumental noise. In particular, we predict the potential for measuring the anisotropy of a broad generalised class of excess radio background models using the low-frequency HERA array as an example. We find that a HERA-like array can distinguish an RSB excess from other sky components based on its angular clustering and spectral dependence, even if these are quite similar to one or more of the other components -- but only in the case that the RSB excess is relatively bright. △ Less

Submitted 1 May, 2024; v1 submitted 20 March, 2024; originally announced March 2024.

Comments: 10 pages, 5 figures

Sensitivity of Bayesian 21 cm power spectrum estimation to foreground model errors

Authors: Jacob Burba, Philip Bull, Michael J. Wilensky, Fraser Kennedy, Hugh Garsden, Katrine A. Glasscock

Abstract: Power spectrum estimators are an important tool in efforts to detect the 21 cm brightness temperature fluctuations from neutral hydrogen at early times. An initial detection will likely be statistical in nature, meaning that it will not be possible to make a coherent map of the brightness temperature fluctuations; instead, only their variance will be measured against a background of noise and resi… ▽ More Power spectrum estimators are an important tool in efforts to detect the 21 cm brightness temperature fluctuations from neutral hydrogen at early times. An initial detection will likely be statistical in nature, meaning that it will not be possible to make a coherent map of the brightness temperature fluctuations; instead, only their variance will be measured against a background of noise and residual systematic effects. Optimal Quadratic Estimator (OQE)-based methods often apply an inverse covariance weighting to the data. However, inaccurate covariance modelling can lead to reduced sensitivity and, in some cases, severe signal loss. We recently proposed a Bayesian method to jointly estimate the 21 cm fluctuations, their power spectrum, and foreground emission. Instead of requiring a fixed a priori estimate of the covariance, we estimate the covariance as part of the inference. Choices of parametrization, particularly of the foregrounds, are subject to model errors and could lead to biases and other ill effects if not properly controlled. In this paper, we investigate the effects of inaccurate foreground models on 21 cm power spectrum recovery. Using simulated visibilities, we find that, even in the most extreme scenarios tested, our approach is capable of recovering 21 cm delay power spectrum estimates consistent with a known input signal for delays $\gtrsim300$ ns ($\sim$88\% of the available Fourier modes). This is true even when using foreground models derived from modified foreground catalogs containing spatial and spectral perturbations at the quoted level of uncertainty on our foreground catalogs. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 13 pages, 9 figures

Statistical estimation of full-sky radio maps from 21cm array visibility data using Gaussian Constrained Realisations

Authors: Katrine A. Glasscock, Philip Bull, Jacob Burba, Hugh Garsden, Michael J. Wilensky

Abstract: An important application of next-generation wide-field radio interferometers is making high dynamic range maps of radio emission. Traditional deconvolution methods like CLEAN can give poor recovery of diffuse structure, prompting the development of wide-field alternatives like Direct Optimal Mapping and $m$-mode analysis. In this paper, we propose an alternative Bayesian method to infer the coeffi… ▽ More An important application of next-generation wide-field radio interferometers is making high dynamic range maps of radio emission. Traditional deconvolution methods like CLEAN can give poor recovery of diffuse structure, prompting the development of wide-field alternatives like Direct Optimal Mapping and $m$-mode analysis. In this paper, we propose an alternative Bayesian method to infer the coefficients of a full-sky spherical harmonic basis for a drift-scan telescope with potentially thousands of baselines. The can precisely encode the uncertainties and correlations between the parameters used to build the recovered image. We use Gaussian Constrained Realisations (GCR) to efficiently draw samples of the spherical harmonic coefficients, despite the very large parameter space and extensive sky-regions of missing data. Each GCR solution provides a complete, statistically-consistent gap-free realisation of a full-sky map conditioned on the available data, even when the interferometer's field of view is small. Many realisations can be generated and used for further analysis and robust propagation of statistical uncertainties. In this paper, we present the mathematical formalism of the spherical harmonic GCR-method for radio interferometers. We focus on the recovery of diffuse emission as a use case, along with validation of the method against simulations with a known diffuse emission component. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 16 pages, 12 figures, submitted to RASTI

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

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

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

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

The foreground transfer function for HI intensity mapping signal reconstruction: MeerKLASS and precision cosmology applications

Authors: Steven Cunnington, Laura Wolz, Philip Bull, Isabella P. Carucci, Keith Grainge, Melis O. Irfan, Yichao Li, Alkistis Pourtsidou, Mario G. Santos, Marta Spinelli, Jingying Wang

Abstract: Blind cleaning methods are currently the preferred strategy for handling foreground contamination in single-dish HI intensity mapping surveys. Despite the increasing sophistication of blind techniques, some signal loss will be inevitable across all scales. Constructing a corrective transfer function using mock signal injection into the contaminated data has been a practice relied on for HI intensi… ▽ More Blind cleaning methods are currently the preferred strategy for handling foreground contamination in single-dish HI intensity mapping surveys. Despite the increasing sophistication of blind techniques, some signal loss will be inevitable across all scales. Constructing a corrective transfer function using mock signal injection into the contaminated data has been a practice relied on for HI intensity mapping experiments. However, assessing whether this approach is viable for future intensity mapping surveys where precision cosmology is the aim, remains unexplored. In this work, using simulations, we validate for the first time the use of a foreground transfer function to reconstruct power spectra of foreground-cleaned low-redshift intensity maps and look to expose any limitations. We reveal that even when aggressive foreground cleaning is required, which causes ${>}\,50\%$ negative bias on the largest scales, the power spectrum can be reconstructed using a transfer function to within sub-percent accuracy. We specifically outline the recipe for constructing an unbiased transfer function, highlighting the pitfalls if one deviates from this recipe, and also correctly identify how a transfer function should be applied in an auto-correlation power spectrum. We validate a method that utilises the transfer function variance for error estimation in foreground-cleaned power spectra. Finally, we demonstrate how incorrect fiducial parameter assumptions (up to ${\pm}100\%$ bias) in the generation of mocks, used in the construction of the transfer function, do not significantly bias signal reconstruction or parameter inference (inducing ${<}\,5\%$ bias in recovered values). △ Less

Submitted 23 May, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 25 pages, 20 figures. See Figure 4 for the main demonstration of the transfer function's performance for reconstructing signal loss from foreground cleaning. Accepted for publication in MNRAS

Mitigating the effect of 1/f noise on the detection of the HI intensity mapping power spectrum from single-dish measurements

Authors: Melis O. Irfan, Yichao Li, Mario G. Santos, Philip Bull, Junhua Gu, Steven Cunnington, Keith Grainge, Jingying Wang

Abstract: We present and compare several methods to mitigate time-correlated (1/f) noise within the HI intensity mapping component of the MeerKAT Large Area Synoptic Survey (MeerKLASS). By simulating scan strategies, the HI signal, foreground emissions, white and correlated noise, we assess the ability of various data processing pipelines to recover the power spectrum of HI brightness temperature fluctuatio… ▽ More We present and compare several methods to mitigate time-correlated (1/f) noise within the HI intensity mapping component of the MeerKAT Large Area Synoptic Survey (MeerKLASS). By simulating scan strategies, the HI signal, foreground emissions, white and correlated noise, we assess the ability of various data processing pipelines to recover the power spectrum of HI brightness temperature fluctuations. We use MeerKAT pilot data to assess the level of 1/f noise expected for the MeerKLASS survey and use these measurements to create realistic levels of time-correlated noise for our simulations. We find the time-correlated noise component within the pilot data to be between 10 and 20 times higher than the white noise level at the scale of k = 0.04 Mpc^-1. Having determined that the MeerKAT 1/f noise is partially correlated across all the frequency channels, we employ Singular Value Decomposition (SVD) as a technique to remove both the 1/f noise and Galactic foregrounds but find that over-cleaning results in the removal of HI power at large (angular and radial) scales; a power loss of 40 per cent is seen for a 3-mode SVD clean at the scale of k = 0.04 Mpc^-1. We compare the impact of map-making using weighting by the full noise covariance (i.e. including a 1/f component), as opposed to just a simple unweighted binning, finding that including the time-correlated noise information reduces the excess power added by 1/f noise by up to 30 per cent. △ Less

Submitted 9 November, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: 14 pages, 12 figures, accepted by MNRAS. This replacement features updated simulations which now include foreground emissions, telescope slew time and a receiver and elevation-dependant temperature contribution to each receiver which changes over frequency

Statistical recovery of 21cm visibilities and their power spectra with Gaussian constrained realisations and Gibbs sampling

Authors: Fraser Kennedy, Philip Bull, Michael Wilensky, Jacob Burba, Samir Choudhuri

Abstract: Radio interferometers designed to probe the 21cm signal from Cosmic Dawn and the Epoch of Reionisation must contend with systematic effects that make it difficult to achieve sufficient dynamic range to separate the 21cm signal from foreground emission and other effects. For instance, the instrument's chromatic response modulates the otherwise spectrally smooth foregrounds, making them difficult to… ▽ More Radio interferometers designed to probe the 21cm signal from Cosmic Dawn and the Epoch of Reionisation must contend with systematic effects that make it difficult to achieve sufficient dynamic range to separate the 21cm signal from foreground emission and other effects. For instance, the instrument's chromatic response modulates the otherwise spectrally smooth foregrounds, making them difficult to model, while a significant fraction of the data must be excised due to the presence of radio frequency interference (RFI), leaving gaps in the data. Errors in modelling the (modulated and gappy) foregrounds can easily generate spurious contamination of what should otherwise be 21cm signal-dominated modes. Various approaches have been developed to mitigate these issues by (e.g.) using non-parametric reconstruction of the foregrounds, in-painting the gaps, and weighting the data to reduce the level of contamination. We present a Bayesian statistical method that combines these approaches, using the coupled techniques of Gaussian constrained realisations (GCR) and Gibbs sampling. This provides a way of drawing samples from the joint posterior distribution of the 21cm signal modes and their power spectrum in the presence of gappy data and an uncertain foreground model in a computationally scalable manner. The data are weighted by an inverse covariance matrix that is estimated as part of the inference, along with a foreground model that can then be marginalised over. We demonstrate the application of this technique on a simulated HERA-like delay spectrum analysis, comparing three different approaches for accounting for the foreground components. △ Less

Submitted 17 March, 2023; v1 submitted 9 November, 2022; originally announced November 2022.

Comments: Updated to match accepted version. Note minor changes to power spectrum plots due to a small bugfix

Bayesian jackknife tests with a small number of subsets: Application to HERA 21cm power spectrum upper limits

Authors: Michael J. Wilensky, Fraser Kennedy, Philip Bull, Joshua S. Dillon, The HERA Collaboration

Abstract: We present a Bayesian jackknife test for assessing the probability that a data set contains biased subsets, and, if so, which of the subsets are likely to be biased. The test can be used to assess the presence and likely source of statistical tension between different measurements of the same quantities in an automated manner. Under certain broadly applicable assumptions, the test is analytically… ▽ More We present a Bayesian jackknife test for assessing the probability that a data set contains biased subsets, and, if so, which of the subsets are likely to be biased. The test can be used to assess the presence and likely source of statistical tension between different measurements of the same quantities in an automated manner. Under certain broadly applicable assumptions, the test is analytically tractable. We also provide an open source code, CHIBORG, that performs both analytic and numerical computations of the test on general Gaussian-distributed data. After exploring the information theoretical aspects of the test and its performance with an array of simulations, we apply it to data from the Hydrogen Epoch of Reionization Array (HERA) to assess whether different sub-seasons of observing can justifiably be combined to produce a deeper 21cm power spectrum upper limit. We find that, with a handful of exceptions, the HERA data in question are statistically consistent and this decision is justified. We conclude by pointing out the wide applicability of this test, including to CMB experiments and the $H_0$ tension. △ Less

Submitted 31 October, 2022; originally announced October 2022.

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

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

Characterizing line-of-sight variability of polarized dust emission with future CMB experiments

Authors: Lisa McBride, Philip Bull, Brandon S. Hensley

Abstract: Galactic dust emission is often accounted for in cosmic microwave background (CMB) analyses by fitting a simple two-parameter modified blackbody (MBB) model in each pixel, which nominally accounts for the temperature and opacity of the dust. While this may be a good approximation for individual dust clouds, typically a number of such clouds are found along each line of sight and within each angula… ▽ More Galactic dust emission is often accounted for in cosmic microwave background (CMB) analyses by fitting a simple two-parameter modified blackbody (MBB) model in each pixel, which nominally accounts for the temperature and opacity of the dust. While this may be a good approximation for individual dust clouds, typically a number of such clouds are found along each line of sight and within each angular pixel, resulting in a superposition of their spectra. In this paper, we study the effects of this superposition on pixel-based foreground fitting strategies by modelling the spectral energy distribution (SED) in each pixel as the integral of individual MBB spectra over various physically-motivated statistical distributions of dust cloud properties. We show that fitting these SEDs with the simple two-parameter MBB model generally results in unbiased estimates of the CMB Stokes Q and U amplitudes in each pixel, unless there are significant changes in both the dust SED and polarization angle along the line of sight, in which case significant ($ > 10σ$) biases are observed in an illustrative model. We also find that the best-fit values of the dust temperature, $T_d$, and spectral index, $β_d$, are significantly biased away from the mean/median of the corresponding statistical distributions when the distributions are broad, suggesting that MBB model fits can give an unrepresentative picture of the physical properties of the dust at microwave wavelengths if not interpreted carefully. Using a Fisher matrix analysis, we also determine the experimental sensitivity required to recover the parameters of the $T_d$ and $β_d$ distributions themselves by fitting a probabilistic MBB model, finding that only the parameters of broad distributions can be measured by SED fitting on a single line of sight. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Comments: 15 pages, 10 figures; see Table 3 for main results

HI intensity mapping with MeerKAT: power spectrum detection in cross-correlation with WiggleZ galaxies

Authors: Steven Cunnington, Yichao Li, Mario G. Santos, Jingying Wang, Isabella P. Carucci, Melis O. Irfan, Alkistis Pourtsidou, Marta Spinelli, Laura Wolz, Paula S. Soares, Chris Blake, Philip Bull, Brandon Engelbrecht, José Fonseca, Keith Grainge, Yin-Zhe Ma

Abstract: We present a detection of correlated clustering between MeerKAT radio intensity maps and galaxies from the WiggleZ Dark Energy Survey. We find a $7.7σ$ detection of the cross-correlation power spectrum, the amplitude of which is proportional to the product of the HI density fraction ($Ω_{\rm HI}$), HI bias ($b_{\rm HI}$) and the cross-correlation coefficient ($r$). We therefore obtain the constrai… ▽ More We present a detection of correlated clustering between MeerKAT radio intensity maps and galaxies from the WiggleZ Dark Energy Survey. We find a $7.7σ$ detection of the cross-correlation power spectrum, the amplitude of which is proportional to the product of the HI density fraction ($Ω_{\rm HI}$), HI bias ($b_{\rm HI}$) and the cross-correlation coefficient ($r$). We therefore obtain the constraint $Ω_{\rm HI} b_{\rm HI} r\,{=}\,[0.86\,{\pm}\,0.10\,({\rm stat})\,{\pm}\,0.12\,({\rm sys})]\,{\times}\,10^{-3}$, at an effective scale of $k_{\rm eff}\,{\sim}\,0.13\,h\,{\rm Mpc}^{-1}$. The intensity maps were obtained from a pilot survey with the MeerKAT telescope, a 64-dish pathfinder array to the SKA Observatory (SKAO). The data were collected from 10.5 hours of observations using MeerKAT's L-band receivers over six nights covering the 11hr field of WiggleZ, in the frequency range $1015-973\,{\rm MHz}$ (0.400$\,{<}\,z\,{<}\,$0.459 in redshift). This detection is the first practical demonstration of the multi-dish auto-correlation intensity mapping technique for cosmology. This marks an important milestone in the roadmap for the cosmology science case with the full SKAO. △ Less

Submitted 25 October, 2022; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 11 pages, 7 figures. See Fig.5 for the main result, showing the cross-power spectrum between MeerKAT and WiggleZ, revealing the 7.7σdetection. v2 is the version accepted for publication by MNRAS. The main changes from v1 were in introducing the radio intensity mapping data

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

Advancing an Interdisciplinary Science of Conversation: Insights from a Large Multimodal Corpus of Human Speech

Authors: Andrew Reece, Gus Cooney, Peter Bull, Christine Chung, Bryn Dawson, Casey Fitzpatrick, Tamara Glazer, Dean Knox, Alex Liebscher, Sebastian Marin

Abstract: People spend a substantial portion of their lives engaged in conversation, and yet our scientific understanding of conversation is still in its infancy. In this report we advance an interdisciplinary science of conversation, with findings from a large, novel, multimodal corpus of 1,656 recorded conversations in spoken English. This 7+ million word, 850 hour corpus totals over 1TB of audio, video,… ▽ More People spend a substantial portion of their lives engaged in conversation, and yet our scientific understanding of conversation is still in its infancy. In this report we advance an interdisciplinary science of conversation, with findings from a large, novel, multimodal corpus of 1,656 recorded conversations in spoken English. This 7+ million word, 850 hour corpus totals over 1TB of audio, video, and transcripts, with moment-to-moment measures of vocal, facial, and semantic expression, along with an extensive survey of speaker post conversation reflections. We leverage the considerable scope of the corpus to (1) extend key findings from the literature, such as the cooperativeness of human turn-taking; (2) define novel algorithmic procedures for the segmentation of speech into conversational turns; (3) apply machine learning insights across various textual, auditory, and visual features to analyze what makes conversations succeed or fail; and (4) explore how conversations are related to well-being across the lifespan. We also report (5) a comprehensive mixed-method report, based on quantitative analysis and qualitative review of each recording, that showcases how individuals from diverse backgrounds alter their communication patterns and find ways to connect. We conclude with a discussion of how this large-scale public dataset may offer new directions for future research, especially across disciplinary boundaries, as scholars from a variety of fields appear increasingly interested in the study of conversation. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: 116 pages; Cross-posted to PsyArxiv; To request data access, go to https://betterup-data-requests.herokuapp.com/

Measurements of the diffuse Galactic synchrotron spectral index and curvature from MeerKLASS pilot data

Authors: Melis O. Irfan, Philip Bull, Mario G. Santos, Jingying Wang, Keith Grainge, Yichao Li, Isabella P. Carucci, Marta Spinelli, Steven Cunnington

Abstract: 21cm intensity mapping experiments are bringing an influx of high spectral resolution observational data in the $\sim100$ MHz $- 1$ GHz regime. We use pilot $971-1075$ MHz data from MeerKAT in single-dish mode, recently used to test the calibration and data reduction scheme of the upcoming MeerKLASS survey, to probe the spectral index of diffuse synchrotron emission below 1 GHz within… ▽ More 21cm intensity mapping experiments are bringing an influx of high spectral resolution observational data in the $\sim100$ MHz $- 1$ GHz regime. We use pilot $971-1075$ MHz data from MeerKAT in single-dish mode, recently used to test the calibration and data reduction scheme of the upcoming MeerKLASS survey, to probe the spectral index of diffuse synchrotron emission below 1 GHz within $145^{\circ} < α< 180^{\circ}$, $-1^{\circ} < δ< 8^{\circ}$. Through comparisons with data from the OVRO Long Wavelength Array and the Maipu and MU surveys, we find an average spectral index of $-2.75 < β< -2.71$ between 45 and 1055 MHz. By fitting for spectral curvature with a spectral index of the form $β+ c \, {\rm{ln}}(ν/ 73~{\rm MHz})$, we measure $β= -2.55 \pm 0.13$ and $c = -0.12 \pm 0.05$ within our target field. Our results are in good agreement (within $1σ$) with existing measurements from experiments such as ARCADE2 and EDGES. These results show the calibration accuracy of current data and demonstrate that MeerKLASS will also be capable of achieving a secondary science goal of probing the interstellar medium. △ Less

Submitted 11 January, 2022; v1 submitted 16 November, 2021; originally announced November 2021.

Comments: 17 pages, 13 figures, accepted for publication in MNRAS. Updated to match published paper (additional references and acknowledgements)

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

HERA Phase I Limits on the Cosmic 21-cm Signal: Constraints on Astrophysics and Cosmology During the Epoch of Reionization

Authors: The HERA Collaboration, Zara Abdurashidova, James E. Aguirre, Paul Alexander, Zaki Ali, Yanga Balfour, Rennan Barkana, Adam Beardsley, Gianni Bernardi, Tashalee Billings, Judd Bowman, Richard Bradley, Phillip Bull, Jacob Burba, Steven Carey, Christopher Carilli, Carina Cheng, David DeBoer, Matthew Dexter, Eloy de Lera Acedo, Joshua Dillon, John Ely, Aaron Ewall-Wice, Nicolas Fagnoni, Anastasia Fialkov , et al. (59 additional authors not shown)

Abstract: Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated… ▽ More Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated above the adiabatic cooling threshold by z~8, independent of uncertainties about the IGM ionization state and the nature of the radio background. Combining HERA limits with galaxy and EoR observations constrains the spin temperature of the z~8 neutral IGM to 27 K < T_S < 630 K (2.3 K < T_S < 640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the CMB dominates the z~8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones (with soft band X-ray luminosities per star formation rate constrained to L_X/SFR = { 10^40.2, 10^41.9 } erg/s/(M_sun/yr) at 68% confidence), consistent with expectations of X-ray binaries in low-metallicity environments. The z~10 limits require even earlier heating if dark-matter interactions (e.g., through millicharges) cool down the hydrogen gas. Using a model in which an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-ray luminosities of L_{r,ν}/SFR > 3.9 x 10^24 W/Hz/(M_sun/yr) and L_X/SFR<10^40 erg/s/(M_sun/yr). The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent EDGES detection. The analysis framework described here provides a foundation for the interpretation of future HERA results. △ Less

Submitted 20 December, 2022; v1 submitted 16 August, 2021; originally announced August 2021.

Comments: 40 pages, 19 figures, accepted to ApJ

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/

Spin-based removal of instrumental systematics in 21cm intensity mapping surveys

Authors: Nialh McCallum, Daniel B. Thomas, Philip Bull, Michael L. Brown

Abstract: Upcoming cosmological intensity mapping surveys will open new windows on the Universe, but they must first overcome a number of significant systematic effects, including polarization leakage. We present a formalism that uses scan strategy information to model the effect of different instrumental systematics on the recovered cosmological intensity signal for `single-dish' (autocorrelation) surveys.… ▽ More Upcoming cosmological intensity mapping surveys will open new windows on the Universe, but they must first overcome a number of significant systematic effects, including polarization leakage. We present a formalism that uses scan strategy information to model the effect of different instrumental systematics on the recovered cosmological intensity signal for `single-dish' (autocorrelation) surveys. This modelling classifies different systematics according to their spin symmetry, making it particularly relevant for dealing with polarization leakage. We show how to use this formalism to calculate the expected contamination from different systematics as a function of the scanning strategy. Most importantly, we show how systematics can be disentangled from the intensity signal based on their spin properties via map-making. We illustrate this, using a set of toy models, for some simple instrumental systematics, demonstrating the ability to significantly reduce the contamination to the observed intensity signal. Crucially, unlike existing foreground removal techniques, this approach works for signals that are non-smooth in frequency, e.g. polarized foregrounds. These map-making approaches are simple to apply and represent an orthogonal and complementary approach to existing techniques for removing systematics from upcoming 21cm intensity mapping surveys. △ Less

Submitted 14 December, 2021; v1 submitted 16 July, 2021; originally announced July 2021.

Comments: 19 pages, 14 Figures, 2 Tables, published in MNRAS

Cleaning foregrounds from single-dish 21cm intensity maps with Kernel Principal Component Analysis

Authors: Melis O. Irfan, Philip Bull

Abstract: The high dynamic range between contaminating foreground emission and the fluctuating 21cm brightness temperature field is one of the most problematic characteristics of 21cm intensity mapping data. While these components would ordinarily have distinctive frequency spectra, making it relatively easy to separate them, instrumental effects and calibration errors further complicate matters by modulati… ▽ More The high dynamic range between contaminating foreground emission and the fluctuating 21cm brightness temperature field is one of the most problematic characteristics of 21cm intensity mapping data. While these components would ordinarily have distinctive frequency spectra, making it relatively easy to separate them, instrumental effects and calibration errors further complicate matters by modulating and mixing them together. A popular class of foreground cleaning method are unsupervised techniques related to Principal Component Analysis (PCA), which exploit the different shapes and amplitudes of each component's contribution to the covariance of the data in order to segregate the signals. These methods have been shown to be effective at removing foregrounds, while also unavoidably filtering out some of the 21cm signal too. In this paper we examine, for the first time in the context of 21cm intensity mapping, a generalised method called Kernel PCA, which instead operates on the covariance of non-linear transformations of the data. This allows more flexible functional bases to be constructed, in principle allowing a cleaner separation between foregrounds and the 21cm signal to be found. We show that Kernel PCA is effective when applied to simulated single-dish (autocorrelation) 21cm data under a variety of assumptions about foregrounds models, instrumental effects etc. It presents a different set of behaviours to PCA, e.g. in terms of sensitivity to the data resolution and smoothing scale, outperforming it on intermediate to large scales in most scenarios. △ Less

Submitted 12 October, 2021; v1 submitted 5 July, 2021; originally announced July 2021.

Comments: 14 pages + 3 new pages, 11 + 5 new figures. Updated to match version accepted by MNRAS

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

Statistical recovery of the BAO scale from multipoles of the beam-convolved 21cm correlation function

Authors: Fraser Kennedy, Philip Bull

Abstract: Despite being designed as an interferometer, the MeerKAT radio array (an SKA pathfinder) can also be used in autocorrelation (`single-dish') mode, where each dish scans the sky independently. Operating in this mode allows extremely high survey speeds to be achieved, albeit at significantly lower angular resolution. We investigate the recovery of the baryon acoustic oscillation (BAO) scale from mul… ▽ More Despite being designed as an interferometer, the MeerKAT radio array (an SKA pathfinder) can also be used in autocorrelation (`single-dish') mode, where each dish scans the sky independently. Operating in this mode allows extremely high survey speeds to be achieved, albeit at significantly lower angular resolution. We investigate the recovery of the baryon acoustic oscillation (BAO) scale from multipoles of the redshift-space correlation function as measured by a low angular resolution 21cm intensity mapping survey of this kind. Our approach is to construct an analytic model of the multipoles of the correlation function and their covariance matrix that includes foreground contamination and beam resolution effects, which we then use to generate an ensemble of mock data vectors from which we attempt to recover the BAO scale. In line with previous studies, we find that recovery of the transverse BAO scale $α_{\perp}$ is hampered by the strong smoothing effect of the instrumental beam with increasing redshift, while the radial scale $α_\parallel$ is much more robust. The multipole formalism naturally incorporates transverse information when it is available however, and so there is no need to perform a radial-only analysis. In particular, the quadrupole of the correlation function preserves a distinctive BAO `bump' feature even for large smoothing scales. We also investigate the robustness of BAO scale recovery to beam model accuracy, severity of the foreground removal cuts, and accuracy of the covariance matrix model, finding in all cases that the radial BAO scale can be recovered in an accurate, unbiased manner. △ Less

Submitted 23 June, 2021; v1 submitted 15 March, 2021; originally announced March 2021.

Comments: Updated to MNRAS accepted version. 20 pages, 14 figures. For the busy reader: see Figs. 4, 5, and 8

Patterns of primary beam non-redundancy in close-packed 21 cm array observations

Authors: Samir Choudhuri, Philip Bull, Hugh Garsden

Abstract: Radio interferometer arrays such as HERA consist of many close-packed dishes arranged in a regular pattern, giving rise to a large number of `redundant' baselines with the same length and orientation. Since identical baselines should see an identical sky signal, this provides a way of finding a relative gain/bandpass calibration without needing an explicit sky model. In reality, there are many rea… ▽ More Radio interferometer arrays such as HERA consist of many close-packed dishes arranged in a regular pattern, giving rise to a large number of `redundant' baselines with the same length and orientation. Since identical baselines should see an identical sky signal, this provides a way of finding a relative gain/bandpass calibration without needing an explicit sky model. In reality, there are many reasons why baselines will not be exactly identical, giving rise to a host of effects that spoil the redundancy of the array and induce spurious structure in the calibration solutions if not accounted for. In this paper, we seek to build an understanding of how differences in the primary beam response between antennas affect redundantly-calibrated interferometric visibilities and their resulting frequency (delay-space) power spectra. We use simulations to study several generic types of primary beam variation, including differences in the width of the main lobe, the angular and frequency structure of the sidelobes, and the beam ellipticity and orientation. For all of these types, we find that additional temporal structure is induced in the gain solutions, particularly when bright point sources pass through the beam. In comparison, only a low level of additional spectral structure is induced. The temporal structure modulates the cosmological 21cm power spectrum, but only at the level of a few percent in our simulations. We also investigate the possibility of signal loss due to decoherence effects when non-redundant visibilities are averaged together, finding that the decoherence is worst when bright point sources pass through the beam, and that its magnitude varies significantly between baseline groups and types of primary beam variation. Redundant calibration absorbs some of the decoherence effect however, reducing its impact compared to if the visibilities were perfectly calibrated. △ Less

Submitted 22 June, 2021; v1 submitted 7 January, 2021; originally announced January 2021.

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

HI intensity mapping with MeerKAT: Calibration pipeline for multi-dish autocorrelation observations

Authors: Jingying Wang, Mario G. Santos, Philip Bull, Keith Grainge, Steven Cunnington, Jose Fonseca, Melis O. Irfan, Yichao Li, Alkistis Pourtsidou, Paula S. Soares, Marta Spinelli, Gianni Bernardi, Brandon Engelbrecht

Abstract: While most purpose-built 21cm intensity mapping experiments are close-packed interferometer arrays, general-purpose dish arrays should also be capable of measuring the cosmological 21cm signal. This can be achieved most efficiently if the array is used as a collection of scanning autocorrelation dishes rather than as an interferometer. As a first step towards demonstrating the feasibility of this… ▽ More While most purpose-built 21cm intensity mapping experiments are close-packed interferometer arrays, general-purpose dish arrays should also be capable of measuring the cosmological 21cm signal. This can be achieved most efficiently if the array is used as a collection of scanning autocorrelation dishes rather than as an interferometer. As a first step towards demonstrating the feasibility of this observing strategy, we show that we are able to successfully calibrate dual-polarisation autocorrelation data from 64 MeerKAT dishes in the L-band (856-1712 MHz, 4096 channels), with 10.5 hours of data retained from six nights of observing. We describe our calibration pipeline, which is based on multi-level RFI flagging, periodic noise diode injection to stabilise gain drifts and an absolute calibration based on a multi-component sky model. We show that it is sufficiently accurate to recover maps of diffuse celestial emission and point sources over a 10 deg x 30 deg patch of the sky overlapping with the WiggleZ 11hr field. The reconstructed maps have a good level of consistency between per-dish maps and external datasets, with the estimated thermal noise limited to 1.4 x the theoretical noise level (~ 2 mK). The residual maps have rms amplitudes below 0.1 K, corresponding to <1% of the model temperature. The reconstructed Galactic HI intensity map shows excellent agreement with the Effelsberg-Bonn HI Survey, and the flux of the radio galaxy 4C+03.18 is recovered to within 3.6%, which demonstrates that the autocorrelation can be successfully calibrated to give the zero-spacing flux and potentially help in the imaging of MeerKAT interferometric data. Our results provide a positive indication towards the feasibility of using MeerKAT and the future SKA to measure the HI intensity mapping signal and probe cosmology on degree scales and above. △ Less

Submitted 13 May, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 24 pages, 34 figures, accepted for publication in MNRAS

Observing relativistic features in large-scale structure surveys -- II: Doppler magnification in an ensemble of relativistic simulations

Authors: Louis Coates, Julian Adamek, Philip Bull, Caroline Guandalin, Chris Clarkson

Abstract: The standard cosmological model is inherently relativistic, and yet a wide range of cosmological observations can be predicted accurately from essentially Newtonian theory. This is not the case on `ultra-large' distance scales, around the cosmic horizon size, however, where relativistic effects can no longer be neglected. In this paper, we present a novel suite of 53 fully relativistic simulations… ▽ More The standard cosmological model is inherently relativistic, and yet a wide range of cosmological observations can be predicted accurately from essentially Newtonian theory. This is not the case on `ultra-large' distance scales, around the cosmic horizon size, however, where relativistic effects can no longer be neglected. In this paper, we present a novel suite of 53 fully relativistic simulations generated using the gevolution code, each covering the full sky out to $z \approx$ 0.85, and approximately 1930 square degrees out to $z \approx$ 3.55. These include a relativistic treatment of massive neutrinos, as well as the gravitational potential that can be used to exactly calculate observables on the past light cone. The simulations are divided into two sets, the first being a set of 39 simulations of the same fiducial cosmology (based on the Euclid Flagship 2 cosmology) with different realisations of the initial conditions, and the second which fixes the initial conditions, but varies each of seven cosmological parameters in turn. Taken together, these simulations allow us to perform statistical studies and calculate derivatives of any relativistic observable with respect to cosmological parameters. As an example application, we compute the cross-correlation between the Doppler magnification term in the convergence, $κ_v$, and the CDM+baryon density contrast, $δ_{\rm cb}$, which arises only in a (special) relativistic treatment. We are able to accurately recover this term as predicted by relativistic perturbation theory, and study its sample variance and derivatives with respect to cosmological parameters. △ Less

Submitted 2 June, 2021; v1 submitted 25 November, 2020; originally announced November 2020.

Comments: 9 pages, 6 figures, corresponds to published version

Journal ref: MNRAS,504(3),2021,3534-3543

Observing relativistic features in large-scale structure surveys -- I: Multipoles of the power spectrum

Authors: Caroline Guandalin, Julian Adamek, Philip Bull, Chris Clarkson, L. Raul Abramo, Louis Coates

Abstract: Planned efforts to probe the largest observable distance scales in future cosmological surveys are motivated by a desire to detect relic correlations left over from inflation, and the possibility of constraining novel gravitational phenomena beyond General Relativity (GR). On such large scales, the usual Newtonian approaches to modelling summary statistics like the power spectrum and bispectrum ar… ▽ More Planned efforts to probe the largest observable distance scales in future cosmological surveys are motivated by a desire to detect relic correlations left over from inflation, and the possibility of constraining novel gravitational phenomena beyond General Relativity (GR). On such large scales, the usual Newtonian approaches to modelling summary statistics like the power spectrum and bispectrum are insufficient, and we must consider a fully relativistic and gauge-independent treatment of observables such as galaxy number counts in order to avoid subtle biases, e.g. in the determination of the $f_{\rm NL}$ parameter. In this work, we present an initial application of an analysis pipeline capable of accurately modelling and recovering relativistic spectra and correlation functions. As a proof of concept, we focus on the non-zero dipole of the redshift-space power spectrum that arises in the cross-correlation of different mass bins of dark matter halos, using strictly gauge-independent observable quantities evaluated on the past light cone of a fully relativistic N-body simulation in a redshift bin $1.7 \le z \le 2.9$. We pay particular attention to the correct estimation of power spectrum multipoles, comparing different methods of accounting for complications such as the survey geometry (window function) and evolution/bias effects on the past light cone, and discuss how our results compare with previous attempts at extracting novel GR signatures from relativistic simulations. △ Less

Submitted 8 January, 2021; v1 submitted 4 September, 2020; originally announced September 2020.

Comments: Minor corrections in the text of version 2

Journal ref: MNRAS, 501(2), 2021, 2547-2561

Searching for dark energy in the matter-dominated era

Authors: Philip Bull, Martin White, Anže Slosar

Abstract: Most efforts to detect signatures of dynamical dark energy are focused on late times, $z \lesssim 2$, where the dark energy component begins to dominate the cosmic energy density. Many theoretical models involving dynamical dark energy exhibit a 'freezing' equation of state however, where $w \to -1$ at late times, with a transition to a 'tracking' behaviour at earlier times (with $w \gg -1$ at suf… ▽ More Most efforts to detect signatures of dynamical dark energy are focused on late times, $z \lesssim 2$, where the dark energy component begins to dominate the cosmic energy density. Many theoretical models involving dynamical dark energy exhibit a 'freezing' equation of state however, where $w \to -1$ at late times, with a transition to a 'tracking' behaviour at earlier times (with $w \gg -1$ at sufficiently high redshift). In this paper, we study whether large-scale structure surveys in the post-reionisation matter-dominated regime, $2 \lesssim z \lesssim 6$, are sensitive to this behaviour, on the basis that the dark energy component should remain detectable (despite being strongly subdominant) in this redshift range given sufficiently precise observations. Using phenomenological models inspired by parameter space studies of Horndeski (generalised scalar-tensor) theories, we show how existing CMB and large-scale structure measurements constrain the dark energy equation of state in the matter-dominated era, and examine how forthcoming galaxy surveys and 21cm intensity mapping instruments can improve constraints in this regime. We also find that the combination of existing CMB and LSS constraints with DESI will already come close to offering the best possible constraints on $H_0$ using BAO/galaxy power spectrum measurements, and that either a spectroscopic follow-up of the LSST galaxy sample (e.g. along the lines of MegaMapper or SpecTel) or a Stage 2/PUMA-like intensity mapping survey, both at $z \gtrsim 2$, would offer better constraints on the class of dark energy models considered here than a comparable cosmic variance-limited galaxy survey at $z \lesssim 1.5$. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: 13 pages, 10 figures; comments welcome

DAYENU: A Simple Filter of Smooth Foregrounds for Intensity Mapping Power Spectra

Authors: Aaron Ewall-Wice, Nicholas Kern, Joshua S. Dillon, Adrian Liu, Aaron Parsons, Saurabh Singh, Adam Lanman, Paul La Plante, Nicolas Fagnoni, Eloy de Lera Acedo, David R. DeBoer, Chuneeta Nunhokee, Philip Bull, Tzu-Ching Chang, T. Joseph Lazio, James Aguirre, Sean Weinberg

Abstract: We introduce DAYENU, a linear, spectral filter for HI intensity mapping that achieves the desirable foreground mitigation and error minimization properties of inverse co-variance weighting with minimal modeling of the underlying data. Beyond 21 cm power-spectrum estimation, our filter is suitable for any analysis where high dynamic-range removal of spectrally smooth foregrounds in irregularly (or… ▽ More We introduce DAYENU, a linear, spectral filter for HI intensity mapping that achieves the desirable foreground mitigation and error minimization properties of inverse co-variance weighting with minimal modeling of the underlying data. Beyond 21 cm power-spectrum estimation, our filter is suitable for any analysis where high dynamic-range removal of spectrally smooth foregrounds in irregularly (or regularly) sampled data is required, something required by many other intensity mapping techniques. Our filtering matrix is diagonalized by Discrete Prolate Spheroidal Sequences which are an optimal basis to model band-limited foregrounds in 21 cm intensity mapping experiments in the sense that they maximally concentrate power within a finite region of Fourier space. We show that DAYENU enables the access of large-scale line-of-sight modes that are inaccessible to tapered DFT estimators. Since these modes have the largest SNRs, DAYENU significantly increases the sensitivity of 21 cm analyses over tapered Fourier transforms. Slight modifications allow us to use DAYENU as a linear replacement for iterative delay CLEANing (DAYENUREST). We refer readers to the Code section at the end of this paper for links to examples and code. △ Less

Submitted 25 October, 2020; v1 submitted 23 April, 2020; originally announced April 2020.

Comments: Accepted to MNRAS, 19 pages, 17 figures. Replacement matches accepted version

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

Mitigating Internal Instrument Coupling II: A Method Demonstration with the Hydrogen Epoch of Reionization Array

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

Abstract: We present a study of internal reflection and cross coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. T… ▽ More We present a study of internal reflection and cross coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA's first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted EoR line-of-sight modes in the range 0.2 < k_parallel < 0.5\ h^-1 Mpc. After systematic removal, we find we can recover these modes in the power spectrum down to the integrated noise-floor of a nightly observation, achieving a dynamic range in the EoR window of 10^-6 in power (mK^2 units) with respect to the bright galactic foreground signal. In the absence of other systematics and assuming the systematic suppression demonstrated here continues to lower noise levels, our results suggest that fully-integrated HERA Phase I may have the capacity to set competitive upper limits on the 21 cm power spectrum. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field. △ Less

Submitted 29 October, 2019; v1 submitted 25 September, 2019; originally announced September 2019.

Comments: Submitted to ApJ

Understanding the HERA Phase I receiver system with simulations and its impact on the detectability of the EoR delay power spectrum

Authors: Nicolas Fagnoni, Eloy de Lera Acedo, David R. DeBoer, 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, Phil Bull, Jacob Burba, Chris L. Carilli, Carina Cheng, Matt Dexter, Joshua S. Dillon, Aaron Ewall-Wice, Randall Fritz, Steve R. Furlanetto, Kingsley Gale-Sides, Brian Glendenning, Deepthi Gorthi , et al. (45 additional authors not shown)

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

Submitted 25 August, 2020; v1 submitted 6 August, 2019; originally announced August 2019.

Comments: 13 pages, 17 figures - Submitted to MNRAS - 2nd revision

Research and Development for HI Intensity Mapping

Authors: Zeeshan Ahmed, David Alonso, Mustafa A. Amin, Réza Ansari, Evan J. Arena, Kevin Bandura, Adam Beardsley, Philip Bull, Emanuele Castorina, Tzu-Ching Chang, Romeel Davé, Joshua S. Dillon, Alexander van Engelen, Aaron Ewall-Wice, Simone Ferraro, Simon Foreman, Josef Frisch, Daniel Green, Gilbert Holder, Daniel Jacobs, Dionysios Karagiannis, Alexander A. Kaurov, Lloyd Knox, Emily Kuhn, Adrian Liu , et al. (29 additional authors not shown)

Abstract: Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and an… ▽ More Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift. △ Less

Submitted 29 July, 2019; originally announced July 2019.

Comments: 10 pages + references, 2 figures, 1 table; APC white paper submitted to the Astro2020 decadal survey. arXiv admin note: text overlap with arXiv:1810.09572

A Roadmap for Astrophysics and Cosmology with High-Redshift 21 cm Intensity Mapping

Authors: The Hydrogen Epoch of Reionization Array, Collaboration, James E. Aguirre, Adam P. Beardsley, Gianni Bernardi, Judd D. Bowman, Philip Bull, Chris L. Carilli, Wei-Ming Dai, David R. DeBoer, Joshua S. Dillon, Aaron Ewall-Wice, Steve R. Furlanetto, Bharat K. Gehlot, Deepthi Gorthi, Bradley Greig, Bryna J. Hazelton, Jacqueline N. Hewitt, Daniel C. Jacobs, Nicholas S. Kern, Piyanat Kittiwisit, Matthew Kolopanis, Paul La Plante, Adrian Liu, Yin-Zhe Ma , et al. (7 additional authors not shown)

Abstract: In this white paper, we lay out a US roadmap for high-redshift 21 cm cosmology (30 < z < 6) in the 2020s. Beginning with the currently-funded HERA and MWA Phase II projects and advancing through the decade with a coordinated program of small-scale instrumentation, software, and analysis projects targeting technology development, this roadmap incorporates our current best understanding of the syste… ▽ More In this white paper, we lay out a US roadmap for high-redshift 21 cm cosmology (30 < z < 6) in the 2020s. Beginning with the currently-funded HERA and MWA Phase II projects and advancing through the decade with a coordinated program of small-scale instrumentation, software, and analysis projects targeting technology development, this roadmap incorporates our current best understanding of the systematics confronting 21 cm cosmology into a plan for overcoming them, enabling next-generation, mid-scale 21 cm arrays to be proposed late in the decade. Submitted for consideration by the Astro2020 Decadal Survey Program Panel for Radio, Millimeter, and Submillimeter Observations from the Ground as a Medium-Sized Project. △ Less

Submitted 15 July, 2019; originally announced July 2019.

Comments: 10 pages (plus a cover page and references), 6 figures. Submitted as a APC White Paper for Astro2020

Modified Gravity and Dark Energy models Beyond $w(z)$CDM Testable by LSST

Authors: Mustapha Ishak, Tessa Baker, Philip Bull, Eske M. Pedersen, Jonathan Blazek, Pedro G. Ferreira, C. Danielle Leonard, Weikang Lin, Eric Linder, Kris Pardo, Georgios Valogiannis

Abstract: One of the main science goals of the Large Synoptic Survey Telescope (LSST) is to uncover the nature of cosmic acceleration. In the base analysis, possible deviations from the Lambda-Cold-Dark-Matter ($Λ$CDM) background evolution will be probed by fitting a $w(z)$CDM model, which allows for a redshift-dependent dark energy equation of state with $w(z)$, within general relativity (GR). A rich array… ▽ More One of the main science goals of the Large Synoptic Survey Telescope (LSST) is to uncover the nature of cosmic acceleration. In the base analysis, possible deviations from the Lambda-Cold-Dark-Matter ($Λ$CDM) background evolution will be probed by fitting a $w(z)$CDM model, which allows for a redshift-dependent dark energy equation of state with $w(z)$, within general relativity (GR). A rich array of other phenomena can arise due to deviations from the standard $Λ$CDM+GR model though, including modifications to the growth rate of structure and lensing, and novel screening effects on non-linear scales. Concrete physical models are needed to provide consistent predictions for these (potentially small) effects, to give us the best chance of detecting them and separating them from astrophysical systematics. A complex plethora of possible models has been constructed over the past few decades, with none emerging as a particular favorite. This document prioritizes a subset of these models along with rationales for further study and inclusion into the LSST Dark Energy Science Collaboration (DESC) data analysis pipelines, based on their observational viability, theoretical plausibility, and level of theoretical development. We provide references and theoretical expressions to aid the integration of these models into DESC software and simulations, and give justifications for why other models were not prioritized. While DESC efforts are free to pursue other models, we provide here guidelines on which theories appear to have higher priority for collaboration efforts due to their perceived promise and greater instructional value. △ Less

Submitted 6 September, 2019; v1 submitted 23 May, 2019; originally announced May 2019.

Comments: 61 pages. Some acknowledgments and references added. This is version-1.1 of an internal collaboration document of LSST-DESC that is being made public and is not planned for submission to a journal

Inflation and Dark Energy from spectroscopy at $z > 2$

Authors: Simone Ferraro, Michael J. Wilson, Muntazir Abidi, David Alonso, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Arturo Avelino, Carlo Baccigalupi, Kevin Bandura, Nicholas Battaglia, Chetan Bavdhankar, José Luis Bernal, Florian Beutler, Matteo Biagetti, Guillermo A. Blanc, Jonathan Blazek, Adam S. Bolton, Julian Borrill, Brenda Frye, Elizabeth Buckley-Geer, Philip Bull, Cliff Burgess, Christian T. Byrnes, Zheng Cai , et al. (118 additional authors not shown)

Abstract: The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at… ▽ More The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $σ(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey