Showing 1–50 of 102 results for author: Fialkov, A

Gravitational waves from mergers of Population III binary black holes: roles played by two evolution channels

Authors: Boyuan Liu, Tilman Hartwig, Nina S. Sartorio, Irina Dvorkin, Guglielmo Costa, Filippo Santoliquido, Anastasia Fialkov, Ralf S. Klessen, Volker Bromm

Abstract: The gravitational wave (GW) signal from binary black hole (BBH) mergers is a promising probe of Population III (Pop III) stars. To fully unleash the power of the GW probe, one important step is to understand the relative importance and features of different BBH evolution channels. We model two channels, isolated binary stellar evolution (IBSE) and nuclear star cluster-dynamical hardening (NSC-DH),… ▽ More The gravitational wave (GW) signal from binary black hole (BBH) mergers is a promising probe of Population III (Pop III) stars. To fully unleash the power of the GW probe, one important step is to understand the relative importance and features of different BBH evolution channels. We model two channels, isolated binary stellar evolution (IBSE) and nuclear star cluster-dynamical hardening (NSC-DH), in one theoretical framework based on the semi-analytical code A-SLOTH, under various assumptions on Pop III initial mass function (IMF), initial binary statistics and high-$z$ nuclear star clusters (NSCs). The NSC-DH channel contributes $\sim 8-95\%$ of Pop III BBH mergers across cosmic history, with higher contributions achieved by initially wider binary stars, more top-heavy IMFs, and more abundant high-$z$ NSCs. The dimensionless stochastic GW background (SGWB) produced by Pop III BBH mergers has peak values $Ω^{\rm peak}_{\rm GW}\sim 10^{-11}-8\times 10^{-11}$ around observer-frame frequencies $ν\sim 10-100\ \rm Hz$. The Pop III contribution can be a non-negligible ($\sim 2-32\%$) component in the total SGWB at $ν\lesssim 10\ \rm Hz$. The estimated detection rates of Pop III BBH mergers by the Einstein Telescope are $\sim 6-230\ \rm yr^{-1}$ and $\sim 30-1230\ \rm yr^{-1}$ for the NSC-DH and IBSE channels, respectively. Pop III BBH mergers in NSCs are more massive than those from IBSE, so they dominate the Pop III SGWB below $20$ Hz in most cases. Besides, the detection rate of Pop III BBH mergers involving at least one intermediate-mass BH above $100\ \rm M_\odot$ by the Einstein Telescope is $\sim 0.5-200\ \rm yr^{-1}$ in NSCs but remains below $0.1\ \rm yr^{-1}$ for IBSE. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 28 + 5 pages, 25 figures, submitted to MNRAS

Understanding spectral artefacts in SKA-LOW 21-cm cosmology experiments: the impact of cable reflections

Authors: Oscar Sage David O'Hara, Fred Dulwich, Eloy de Lera Acedo, Jiten Dhandha, Thomas Gessey-Jones, Dominic Anstey, Anastasia Fialkov

Abstract: The Cosmic Dawn marks the first star formations and preceded the Epoch-of-Reionization, when the Universe underwent a fundamental transformation propelled by the radiation from these first stars and galaxies. Interferometric 21-cm experiments aim to probe redshifted neutral hydrogen signals from these periods, constraining the conditions of the early Universe. The SKA-LOW instrument of the Square… ▽ More The Cosmic Dawn marks the first star formations and preceded the Epoch-of-Reionization, when the Universe underwent a fundamental transformation propelled by the radiation from these first stars and galaxies. Interferometric 21-cm experiments aim to probe redshifted neutral hydrogen signals from these periods, constraining the conditions of the early Universe. The SKA-LOW instrument of the Square Kilometre Array telescope is envisaged to be the largest and most sensitive radio telescope at m and cm wavelengths. The latest Aperture Array Verification Systems feature 7m coaxial transmission lines connecting the Low Noise Amplifiers to optical transmitters at the front of the analogue-receiving chain. An impedance mismatch between these components results in a partially reflected electromagnetic signal, which introduces chromatic aberrations in the instrument bandpass. This causes power from the foreground signals to appear at higher delays, potentially contaminating the EoR window, a region at which the 21-cm signal should be detectable. We present an end-to-end simulation pipeline for SKA-LOW using a composite sky model combining radio foregrounds from The GLEAM Survey, Haslam $408$MHz, and a $1.5$cGpc 21-cm brightness temperature cube generated with the 21cmSPACE simulator. Iterating a parametric approach, we derive a model for the scattering parameters of a coaxial transmission line in terms of its specifications and bulk material properties. Assuming identical cables of length $\leq 15.0$m with impedance mismatch $\leq 10Ω$ confines the reflection to k-modes below the EoR window. However, we demonstrate that even a $0.1$% length tolerance introduces contamination with an RMSE of $\sim 10$% across all accessible k-modes. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: 16 pages, 11 figures

Constraining the clustering and 21-cm signature of radio galaxies at cosmic dawn

Authors: Sudipta Sikder, Rennan Barkana, Anastasia Fialkov

Abstract: The efficiency of radio emission is an important unknown parameter of early galaxies at cosmic dawn, as models with high efficiency have been shown to modify the cosmological 21-cm signal substantially, deepening the absorption trough and boosting the 21-cm power spectrum. Such models have been previously directly constrained by the overall extragalactic radio background as observed by ARCADE-2 an… ▽ More The efficiency of radio emission is an important unknown parameter of early galaxies at cosmic dawn, as models with high efficiency have been shown to modify the cosmological 21-cm signal substantially, deepening the absorption trough and boosting the 21-cm power spectrum. Such models have been previously directly constrained by the overall extragalactic radio background as observed by ARCADE-2 and LWA-1. In this work, we constrain the clustering of high redshift radio sources by utilizing the observed upper limits on arcminute-scale anisotropy from the VLA at 4.9~GHz and ATCA at 8.7~GHz. Using a semi-numerical simulation of a plausible astrophysical model for illustration, we show that the clustering constraints on the radio efficiency are much stronger than those from the overall background intensity, by a factor that varies from 12 at redshift 7 to 30 at redshift 22. As a result, the predicted maximum depth of the global 21-cm signal is lowered by a factor of 5 (to 1700~mK), and the maximum 21-cm power spectrum peak at cosmic dawn is lowered by a factor of 24 (to $2\times 10^5$~mK$^2$). We conclude that the observed clustering is the strongest current direct constraint on such models, but strong early radio emission from galaxies remains viable for producing a strongly enhanced 21-cm signal from cosmic dawn. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Comments: 8 pages, 3 figures

On the Constraints on Superconducting Cosmic Strings from 21-cm Cosmology

Authors: T. Gessey-Jones, S. Pochinda, H. T. J. Bevins, A. Fialkov, W. J. Handley, E. de Lera Acedo, S. Singh, R. Barkana

Abstract: Constraints on the potential properties of superconducting cosmic strings provide an indirect probe of physics beyond the standard model at energies inaccessible to terrestrial particle colliders. In this study, we perform the first joint Bayesian analysis to extract constraints on superconducting cosmic strings from current 21-cm signal measurements while accounting rigorously for the uncertainti… ▽ More Constraints on the potential properties of superconducting cosmic strings provide an indirect probe of physics beyond the standard model at energies inaccessible to terrestrial particle colliders. In this study, we perform the first joint Bayesian analysis to extract constraints on superconducting cosmic strings from current 21-cm signal measurements while accounting rigorously for the uncertainties in foregrounds and high redshift astrophysics. We include the latest publicly available 21-cm power spectrum upper limits from HERA, 21-cm global signal data from SARAS 3, and the synergistic probe of the unresolved X-ray background in our final analysis. This paper thus constitutes the first attempt to use 21-cm power spectrum data to probe cosmic strings. In contrast to previous works, we find no strong constraints can be placed on superconducting cosmic strings from current 21-cm measurements. This is because of uncertainties in the X-ray emission efficiency of the first galaxies, with X-ray emissivities greater than $3 \times 10^{40}$erg s$^{-1}$ M$_{\odot}^{-1}$ yr able to mask the presence of cosmic strings in the 21-cm signal. We conclude by discussing the prospects for future constraints from definitive 21-cm signal measurements and argue that the recently proposed soft photon heating should be cause for optimism due to its potential to break degeneracies that would have otherwise made the signatures of cosmic strings difficult to distinguish from those of astrophysical origin. △ Less

Submitted 15 February, 2024; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: 18 pages, 10 figures. Companion paper to arXiv:2312.08095. Accepted for publication in MNRAS, updated to accepted version

Constraining the properties of Population III galaxies with multi-wavelength observations

Authors: S. Pochinda, T. Gessey-Jones, H. T. J. Bevins, A. Fialkov, S. Heimersheim, I. Abril-Cabezas, E. de Lera Acedo, S. Singh, S. Sikder, R. Barkana

Abstract: The early Universe, spanning 400,000 to 400 million years after the Big Bang ($z\approx1100-11$), has been left largely unexplored as the light from luminous objects is too faint to be observed directly. While new experiments are pushing the redshift limit of direct observations, measurements in the low-frequency radio band promise to probe early star and black hole formation via observations of t… ▽ More The early Universe, spanning 400,000 to 400 million years after the Big Bang ($z\approx1100-11$), has been left largely unexplored as the light from luminous objects is too faint to be observed directly. While new experiments are pushing the redshift limit of direct observations, measurements in the low-frequency radio band promise to probe early star and black hole formation via observations of the hydrogen 21-cm line. In this work we explore synergies between 21-cm data from the HERA and SARAS 3 experiments and observations of the unresolved radio and X-ray backgrounds using multi-wavelength Bayesian analysis. We use the combined data set to constrain properties of Population II and Population III stars as well as early X-ray and radio sources. The joint fit reveals a 68 percentile disfavouring of Population III star formation efficiencies $\gtrsim5.7\%$. We also show how the 21-cm and the X-ray background data synergistically constrain opposite ends of the X-ray efficiency prior distribution to produce a peak in the 1D posterior of the X-ray luminosity per star formation rate. We find (at 68\% confidence) that early galaxies were likely 0.3 to 318 times as X-ray efficient as present-day starburst galaxies. We also show that the functional posteriors from our joint fit rule out global 21-cm signals deeper than $\lesssim-203\ \mathrm{mK}$ and power spectrum amplitudes at $k=0.34\ h\mathrm{Mpc^{-1}}$ greater than $Δ_{21}^2 \gtrsim 946\ \mathrm{mK}^2$ with $3σ$ confidence. △ Less

Submitted 1 May, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

FlexKnot as a Generalised Model of the Sky-averaged 21-cm Signal at z ~ 6-30 in the Presence of Systematics

Authors: Emma Shen, Dominic Anstey, Eloy de Lera Acedo, Anastasia Fialkov

Abstract: Global 21-cm experiments are built to study the evolution of the Universe between the cosmic dawn and the epoch of reionisation. FlexKnot is a function parameterised by freely moving knots stringed together by splines. Adopting the FlexKnot function as the signal model has the potential to separate the global 21-cm signal from the foregrounds and systematics while being capable of recovering the c… ▽ More Global 21-cm experiments are built to study the evolution of the Universe between the cosmic dawn and the epoch of reionisation. FlexKnot is a function parameterised by freely moving knots stringed together by splines. Adopting the FlexKnot function as the signal model has the potential to separate the global 21-cm signal from the foregrounds and systematics while being capable of recovering the crucial features given by theoretical predictions. In this paper, we implement the FlexKnot method by integrating twice over a function of freely moving knots interpolated linearly. The function is also constrained at the lower frequencies corresponding to the dark ages by theoretical values. The FlexKnot model is tested in the framework of the realistic data analysis pipeline of the REACH global signal experiment using simulated antenna temperature data. We demonstrate that the FlexKnot model performs better than existing signal models, e.g. the Gaussian signal model, at reconstructing the shape of the true signals present in the simulated REACH data, especially for injected signals with complex structures. The capabilities of the FlexKnot signal model is also tested by introducing various systematics and simulated global signals of different types. These tests show that four to five knots are sufficient to recover the general shape of most realistic injected signals, with or without sinusoidal systematics. We show that true signals whose absorption trough of amplitude between 120 to 450 mK can be well recovered with systematics up to about 50 mK. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 12 pages, 11 figures

Cosmic mysteries and the hydrogen 21-cm line: bridging the gap with lunar observations

Authors: Anastasia Fialkov, Thomas Gessey-Jones, Jiten Dhandha

Abstract: The hydrogen 21-cm signal is predicted to be the richest probe of the young Universe including eras known as the cosmic Dark Ages, the Cosmic Dawn when the first star and black hole formed, and the Epoch of Reionization. This signal holds the key to deciphering processes that take place at the early stages of cosmic history. In this opinion piece, we discuss the potential scientific merit of lunar… ▽ More The hydrogen 21-cm signal is predicted to be the richest probe of the young Universe including eras known as the cosmic Dark Ages, the Cosmic Dawn when the first star and black hole formed, and the Epoch of Reionization. This signal holds the key to deciphering processes that take place at the early stages of cosmic history. In this opinion piece, we discuss the potential scientific merit of lunar observations of the 21-cm signal and their advantages over more affordable terrestrial efforts. The moon is a prime location for radio cosmology which will enable precision observations of the low-frequency radio sky. The uniqueness of such observations is that they will provide an unparalleled opportunity to test cosmology and the nature of dark matter using the Dark Ages 21-cm signal. No less enticing is the opportunity to obtain a much clearer picture of Cosmic Dawn than what is achievable from the ground, which will allow us to probe properties of the first stars and black holes. △ Less

Submitted 17 January, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 15 pages, 3 figures, accepted for publication in Philosophical Transactions A

First upper limits on the 21 cm signal power spectrum from cosmic dawn from one night of observations with NenuFAR

Authors: S. Munshi, F. G. Mertens, L. V. E. Koopmans, A. R. Offringa, B. Semelin, D. Aubert, R. Barkana, A. Bracco, S. A. Brackenhoff, B. Cecconi, E. Ceccotti, S. Corbel, A. Fialkov, B. K. Gehlot, R. Ghara, J. N. Girard, J. M. Grießmeier, C. Höfer, I. Hothi, R. Mériot, M. Mevius, P. Ocvirk, A. K. Shaw, G. Theureau, S. Yatawatta , et al. (2 additional authors not shown)

Abstract: The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nançay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift of $z$ = 2… ▽ More The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nançay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift of $z$ = 20.3. NenuFAR is a new low-frequency radio interferometer, operating in the 10-85 MHz frequency range, currently under construction at the Nançay Radio Observatory in France. It is a phased array instrument with a very dense uv coverage at short baselines, making it one of the most sensitive instruments for 21 cm cosmology analyses at these frequencies. Our analysis adopts the foreground subtraction approach, in which sky sources are modeled and subtracted through calibration and residual foregrounds are subsequently removed using Gaussian process regression. The final power spectra are constructed from the gridded residual data cubes in the uv plane. Signal injection tests are performed at each step of the analysis pipeline, the relevant pipeline settings are optimized to ensure minimal signal loss, and any signal suppression is accounted for through a bias correction on our final upper limits. We obtain a best 2$σ$ upper limit of $2.4\times 10^7$ $\text{mK}^{2}$ at $z$ = 20.3 and $k$ = 0.041 $h\,\text{cMpc}^{-1}$. We see a strong excess power in the data, making our upper limits two orders of magnitude higher than the thermal noise limit. We investigate the origin and nature of this excess power and discuss further improvements to the analysis pipeline that can potentially mitigate it and consequently allow us to reach thermal noise sensitivity when multiple nights of observations are processed in the future. △ Less

Submitted 30 April, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 27 pages, 21 figures, and 6 tables; accepted for publication in Astronomy and Astrophysics (A&A); language edits implemented; typos corrected

Journal ref: A&A 681, A62 (2024)

Constraining exotic dark matter models with the dark ages 21-cm signal

Authors: Rajesh Mondal, Rennan Barkana, Anastasia Fialkov

Abstract: The dark ages 21-cm signal is a powerful tool for precision cosmology and probing new physics. We study two non-standard models: an excess radio background (ERB) model (possibly generated by dark matter decay) and the millicharged dark matter (mDM) model. These models were inspired by the possible EDGES detection of a strong global 21-cm absorption during cosmic dawn, but more generally they provi… ▽ More The dark ages 21-cm signal is a powerful tool for precision cosmology and probing new physics. We study two non-standard models: an excess radio background (ERB) model (possibly generated by dark matter decay) and the millicharged dark matter (mDM) model. These models were inspired by the possible EDGES detection of a strong global 21-cm absorption during cosmic dawn, but more generally they provide a way to anticipate the potential discovery space. During the dark ages the 21-cm global signal in the ERB model reaches a saturated form for an amplitude $A_{\rm r}=0.4$, where $A_{\rm r}$ is the radio background intensity at cosmic dawn relative to the cosmic microwave background. This amplitude is one-fifth of the minimum required to explain the EDGES signal, and corresponds to just 0.1% of the observed extragalactic background; it would give a signal that can be detected at 5.9$σ$ significance (compared to $4.1\,σ$ for the standard signal) and can be distinguished from the standard (no ERB) signal at $8.5\,σ$, all with a 1,000 hr global signal measurement. The 21-cm power spectrum has potentially more information, but far greater resources would be required for comparable constraints. For the mDM model, over a range of viable parameters, the global signal detection significance would be $4.7-7.2\,σ$, and it could be distinguished from the standard at $2.2-9.3\,σ$. With an array of global signal antennas achieving an effective 100,000 hr integration, the significance would be $10\,\times$ better. Our analysis helps motivate the development of lunar and space-based dark ages experiments. △ Less

Submitted 8 November, 2023; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: Published in Monthly Notices of the Royal Astronomical Society

Journal ref: MNRAS, Volume 527, Issue 1, January 2024, Pages 1461-1471

Modeling Post-Reionization HI Distributions in Fuzzy Dark Matter Cosmologies Using Conditional Normalizing Flows

Authors: Tibor Dome, Rumail Azhar, Anastasia Fialkov

Abstract: Upcoming 21 cm intensity mapping experiments like the Square Kilometer Array (SKA) hold significant potential to constrain the properties of dark matter. In this work, we model neutral hydrogen (HI) distributions using high-resolution hydrodynamical $N$-body simulations of both cold dark matter (CDM) and fuzzy dark matter (FDM) cosmologies in the post-reionization redshift range of $z=3.42-4.94$.… ▽ More Upcoming 21 cm intensity mapping experiments like the Square Kilometer Array (SKA) hold significant potential to constrain the properties of dark matter. In this work, we model neutral hydrogen (HI) distributions using high-resolution hydrodynamical $N$-body simulations of both cold dark matter (CDM) and fuzzy dark matter (FDM) cosmologies in the post-reionization redshift range of $z=3.42-4.94$. We show that the HI abundance decreases in FDM-like cosmologies. Extreme FDM models with $m\sim 10^{-22}$ eV are at odds with a range of measurements. Due to the increased halo bias, the HI bias increases, paralleled by the damped Lyman-$α$ (DLA) bias which we infer from the cross-section of DLAs. The distribution of the latter in extreme FDM models has a high median at the low-mass end, which can be traced to the high column density of cosmic filaments. FDM models exhibit a very similar abundance of DLAs compared to CDM while sub-DLAs are already less abundant. We study the prospects of detecting the brightest HI peaks with SKA1-Low at $z=4.94$, indicating moderate signal-to-noise ratios (SNR) at angular resolution $θ_A = 2^{\prime}$ with a rapidly declining SNR for lower values of $θ_{A}$. After training the conditional normalizing flow network HIGlow on 2D HI maps, we interpolate its latent space of axion masses to predict the peak flux for a new, synthetic FDM cosmology, finding good agreement with expectations. This work thus underscores the potential of normalizing flows in capturing complex, non-linear structures within HI maps, offering a versatile tool for conditional sample generation and prediction tasks. △ Less

Submitted 30 December, 2023; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 18 pages, 10 figures, 4 tables

Journal ref: MNRAS, 527, 10397-10415 (2023)

FlexKnot and Gaussian Process for 21 cm global signal analysis and foreground separation

Authors: Stefan Heimersheim, Leiv Rønneberg, Henry Linton, Filippo Pagani, Anastasia Fialkov

Abstract: The cosmological 21 cm signal is one of the most promising avenues to study the Epoch of Reionization. One class of experiments aiming to detect this signal is global signal experiments measuring the sky-averaged 21 cm brightness temperature as a function of frequency. A crucial step in the interpretation and analysis of such measurements is separating foreground contributions from the remainder o… ▽ More The cosmological 21 cm signal is one of the most promising avenues to study the Epoch of Reionization. One class of experiments aiming to detect this signal is global signal experiments measuring the sky-averaged 21 cm brightness temperature as a function of frequency. A crucial step in the interpretation and analysis of such measurements is separating foreground contributions from the remainder of the signal, requiring accurate models for both components. Current models for the signal (non-foreground) component, which may contain cosmological and systematic contributions, are incomplete and unable to capture the full signal. We propose two new methods for extracting this component from the data: Firstly, we employ a foreground-orthogonal Gaussian Process to extract the part of the signal that cannot be explained by the foregrounds. Secondly, we use a FlexKnot parameterization to model the full signal component in a free-form manner, not assuming any particular shape or functional form. This method uses Bayesian model selection to find the simplest signal that can explain the data. We test our methods on both, synthetic data and publicly available EDGES low-band data. We find that the Gaussian Process can clearly capture the foreground-orthogonal signal component of both data sets. The FlexKnot method correctly recovers the full shape of the input signal used in the synthetic data and yields a multi-modal distribution of different signal shapes that can explain the EDGES observations. △ Less

Submitted 17 December, 2023; v1 submitted 9 October, 2023; originally announced October 2023.

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

Population synthesis of Be X-ray binaries: metallicity dependence of total X-ray outputs

Authors: Boyuan Liu, Nina S. Sartorio, Robert G. Izzard, Anastasia Fialkov

Abstract: X-ray binaries (XRBs) are thought to regulate cosmic thermal and ionization histories during the Epoch of Reionization and Cosmic Dawn ($z\sim 5-30$).Theoretical predictions of the X-ray emission from XRBs are important for modelling such early cosmic evolution. Nevertheless, the contribution from Be-XRBs, powered by accretion of compact objects from decretion disks around rapidly rotating O/B sta… ▽ More X-ray binaries (XRBs) are thought to regulate cosmic thermal and ionization histories during the Epoch of Reionization and Cosmic Dawn ($z\sim 5-30$).Theoretical predictions of the X-ray emission from XRBs are important for modelling such early cosmic evolution. Nevertheless, the contribution from Be-XRBs, powered by accretion of compact objects from decretion disks around rapidly rotating O/B stars, has not been investigated systematically. Be-XRBs are the largest class of high-mass XRBs (HMXBs) identified in local observations and are expected to play even more important roles in metal-poor environments at high redshifts. In light of this, we build a physically motivated model for Be-XRBs based on recent hydrodynamic simulations and observations of decretion disks. Our model is able to reproduce the observed population of Be-XRBs in the Small Magellanic Cloud with appropriate initial conditions and binary stellar evolution parameters. We derive the X-ray output from Be-XRBs as a function of metallicity in the (absolute) metallicity range $Z\in [10^{-4},0.03]$ with a large suite of binary population synthesis (BPS) simulations. The simulated Be-XRBs can explain a non-negligible fraction ($\gtrsim 30\%$) of the total X-ray output from HMXBs observed in nearby galaxies for $Z\sim 0.0003-0.02$. The X-ray luminosity per unit star formation rate from Be-XRBs in our fiducial model increases by a factor of $\sim 8$ from $Z=0.02$ to $Z=0.0003$, which is similar to the trend seen in observations of all types of HMXBs. We conclude that Be-XRBs are potentially important X-ray sources that deserve greater attention in BPS of XRBs. △ Less

Submitted 9 November, 2023; v1 submitted 11 August, 2023; originally announced August 2023.

Comments: 26 pages, 16+7 figures, accepted for publication in MNRAS, added discussions on supernova kicks, stellar winds and treatments of super-Eddington accretion and angular momentum loss during mass transfer

Receiver design for the REACH global 21-cm signal experiment

Authors: Nima Razavi-Ghods, Ian L. V. Roque, Steven H. Carey, John A. Ely, Will Handley, Alessio Magro, Riccardo Chiello, Tian Huang, P. Alexander, D. Anstey, G. Bernardi, H. T. J. Bevins, J. Cavillot, W. Croukamp, J. Cumner, E. de Lera Acedo, D. I. L. de Villiers, A. Fialkov, T. Gessey-Jones, Q. Gueuning, A. T. Josaitis, G. Kulkarni, S. A. K. Leeney, R. Maiolino, P. D. Meerburg , et al. (13 additional authors not shown)

Abstract: We detail the the REACH radiometric system designed to enable measurements of the 21-cm neutral hydrogen line. Included is the radiometer architecture and end-to-end system simulations as well as a discussion of the challenges intrinsic to highly-calibratable system development. Following this, we share laboratory results based on the calculation of noise wave parameters utilising an over-constrai… ▽ More We detail the the REACH radiometric system designed to enable measurements of the 21-cm neutral hydrogen line. Included is the radiometer architecture and end-to-end system simulations as well as a discussion of the challenges intrinsic to highly-calibratable system development. Following this, we share laboratory results based on the calculation of noise wave parameters utilising an over-constrained least squares approach demonstrating a calibration RMSE of 80 mK for five hours of integration on a custom-made source with comparable impedance to that of the antenna used in the field. This paper therefore documents the state of the calibrator and data analysis in December 2022 in Cambridge before shipping to South Africa. △ Less

Submitted 14 July, 2023; v1 submitted 30 June, 2023; originally announced July 2023.

Comments: 30 pages, 19 figures

Mini-quenching of $z=4-8$ galaxies by bursty star formation

Authors: Tibor Dome, Sandro Tacchella, Anastasia Fialkov, Daniel Ceverino, Avishai Dekel, Omri Ginzburg, Sharon Lapiner, Tobias J. Looser

Abstract: The recent reported discovery of a low-mass $z=5.2$ and an intermediate-mass $z=7.3$ quenched galaxy with JWST/NIRSpec is the first evidence of halted star formation above $z\approx 5$. Here we show how bursty star formation at $z=4-8$ gives rise to temporarily quenched, or mini-quenched galaxies in the mass range $M_{\star} = 10^7-10^9 \ M_{\odot}$ using four models of galaxy formation: the perio… ▽ More The recent reported discovery of a low-mass $z=5.2$ and an intermediate-mass $z=7.3$ quenched galaxy with JWST/NIRSpec is the first evidence of halted star formation above $z\approx 5$. Here we show how bursty star formation at $z=4-8$ gives rise to temporarily quenched, or mini-quenched galaxies in the mass range $M_{\star} = 10^7-10^9 \ M_{\odot}$ using four models of galaxy formation: the periodic box simulation IllustrisTNG, the zoom-in simulations VELA and FirstLight and an empirical halo model. The main causes for mini-quenching are stellar feedback, lack of gas accretion onto galaxies and galaxy-galaxy interactions. The abundance of (mini-)quenched galaxies agrees across the models: the population first appears below $z\approx 8$, after which their proportion increases with cosmic time, from $\sim 0.5-1.0$% at $z=7$ to $\sim 2-4$% at $z=4$, corresponding to comoving number densities of $\sim 10^{-5}$ Mpc$^{-3}$ and $\sim 10^{-3}$ Mpc$^{-3}$, respectively. These numbers are consistent with star formation rate duty cycles inferred for VELA and FirstLight galaxies. Their star formation histories (SFHs) suggest that mini-quenching at $z=4-8$ is short-lived with a duration of $\sim 20-40$ Myr, which is close to the free-fall timescale of the inner halo. However, mock spectral energy distributions of mini-quenched galaxies in IllustrisTNG and VELA do not match JADES-GS-z7-01-QU photometry, unless their SFHs are artificially altered to be more bursty on timescales of $\sim 40$ Myr. Studying mini-quenched galaxies might aid in calibrating sub-grid models governing galaxy formation, as these may not generate sufficient burstiness at high redshift to explain the SFH inferred for JADES-GS-z7-01-QU. △ Less

Submitted 21 November, 2023; v1 submitted 11 May, 2023; originally announced May 2023.

Comments: 13 pages, 5 figures, 5 tables

Journal ref: MNRAS, 527, 2139-2151 (2023)

Signatures of Cosmic Ray Heating in 21-cm Observables

Authors: T. Gessey-Jones, A. Fialkov, E. de Lera Acedo, W. J. Handley, R. Barkana

Abstract: Cosmic rays generated by supernovae carry away a significant portion of the lifetime energy emission of their parent star, making them a plausible mechanism for heating the early universe intergalactic medium (IGM). Following a review of the existing literature on cosmic ray heating, we develop a flexible model of this heating mechanism for use in 3D semi-numerical 21-cm signal simulations and con… ▽ More Cosmic rays generated by supernovae carry away a significant portion of the lifetime energy emission of their parent star, making them a plausible mechanism for heating the early universe intergalactic medium (IGM). Following a review of the existing literature on cosmic ray heating, we develop a flexible model of this heating mechanism for use in 3D semi-numerical 21-cm signal simulations and conduct the first investigations of the signatures it imprints on the 21-cm power spectrum and tomographic maps. We find that cosmic ray heating of the IGM is short-ranged, leading to heating clustered around star-forming sites, and a sharp contrast between heated regions of 21-cm emission and unheated regions of absorption. This contrast results in greater small-scale power for cosmic ray heated scenarios compared to what is found for X-ray heating, thus suggesting a way to test the nature of IGM heating with future 21-cm observations. Finally, we find an unexpectedly rich thermal history in models where cosmic rays can only escape efficiently from low-mass halos, such as in scenarios where these energetic particles originate from population III star supernovae remnants. The interplay of heating and the Lyman-Werner feedback in these models can produce a local peak in the IGM kinetic temperature and, for a limited parameter range, a flattened absorption trough in the global 21-cm signal. △ Less

Submitted 3 October, 2023; v1 submitted 14 April, 2023; originally announced April 2023.

Comments: 23 pages, 10 figures. Accepted for publication in MNRAS, updated to accepted version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 526, Issue 3, December 2023, Pages 4262-4284

Population III X-ray Binaries and their Impact on the Early Universe

Authors: Nina S. Sartorio, A. Fialkov, T. Hartwig, G. M. Mirouh, R. G. Izzard, M. Magg, R. S. Klessen, S. C. O. Glover, L. Chen, Y. Tarumi, D. D. Hendriks

Abstract: The first population of X-ray binaries (XRBs) is expected to affect the thermal and ionization states of the gas in the early Universe. Although these X-ray sources are predicted to have important implications for high-redshift observable signals, such as the hydrogen 21-cm signal from cosmic dawn and the cosmic X-ray background, their properties are poorly explored, leaving theoretical models lar… ▽ More The first population of X-ray binaries (XRBs) is expected to affect the thermal and ionization states of the gas in the early Universe. Although these X-ray sources are predicted to have important implications for high-redshift observable signals, such as the hydrogen 21-cm signal from cosmic dawn and the cosmic X-ray background, their properties are poorly explored, leaving theoretical models largely uninformed. In this paper we model a population of X-ray binaries arising from zero metallicity stars. We explore how their properties depend on the adopted initial mass function (IMF) of primordial stars, finding a strong effect on their number and X-ray production efficiency. We also present scaling relations between XRBs and their X-ray emission with the local star formation rate, which can be used in sub-grid models in numerical simulations to improve the X-ray feedback prescriptions. Specifically, we find that the uniformity and strength of the X-ray feedback in the intergalactic medium is strongly dependant on the IMF. Bottom-heavy IMFs result in a smoother distribution of XRBs, but have a luminosity orders of magnitude lower than more top-heavy IMFs. Top-heavy IMFs lead to more spatially uneven, albeit strong, X-ray emission. An intermediate IMF has a strong X-ray feedback while sustaining an even emission across the intergalactic medium. These differences in X-ray feedback could be probed in the future with measurements of the cosmic dawn 21-cm line of neutral hydrogen, which offers us a new way of constraining population III IMF. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: Accepted for publication in MNRAS, 17 pages, 9 figures

Accurate modelling of the Lyman-$α$ coupling for the 21-cm signal, observability with NenuFAR and SKA

Authors: Benoit Semelin, Romain Mériot, Florent Mertens, Léon V. E. Koopmans, Dominique Aubert, Rennan Barkana, Anastasia Fialkov, Satyapan Munshi, Pierre Ocvirk

Abstract: The measurement of the $21$ cm signal from the Cosmic Dawn is a major goal for several existing and upcoming radio interferometers such as NenuFAR and the SKA. During this era before the beginning of the Epoch of Reionization, the signal is more difficult to observe due to brighter foregrounds but reveals additional information on the underlying astrophysical processes encoded in the spatial fluct… ▽ More The measurement of the $21$ cm signal from the Cosmic Dawn is a major goal for several existing and upcoming radio interferometers such as NenuFAR and the SKA. During this era before the beginning of the Epoch of Reionization, the signal is more difficult to observe due to brighter foregrounds but reveals additional information on the underlying astrophysical processes encoded in the spatial fluctuations of the spin temperature of hydrogen. To interpret future measurements, controlling the level of accuracy of the Lyman-$α$ flux modelling is mandatory. In this work, we evaluate the impact of various approximations that exist in the main fast modelling approach compared to the results of a costly full radiative transfer simulation. The fast SPINTER code, presented in this work, computes the Lyman-$α$ flux including the effect of wing scatterings for an inhomogeneous emissivity field, but assuming an otherwise homogeneous expanding universe. The LICORICE code computes the full radiative transfer in the Lyman-$α$ line without any substantial approximation. We find that the difference between homogeneous and inhomogeneous gas density and temperature is very small for the computed flux. On the contrary, neglecting the effect of gas velocities produces a significant change in the computed flux. We identify the causes (mainly Doppler shifts due to velocity gradients) and quantify the magnitude of the effect in both an idealised setup and a realistic cosmological situation. We find that the amplitude of the effect, up to a factor of $\sim 2$ on the $21$ cm signal power spectrum on some scales (depending on both other model parameters and the redshift), can be easily discriminated with an SKA-like survey and already be approached, particularly for exotic signals, by the ongoing NenuFAR Cosmic Dawn Key Science Program. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 14 pages, accepted for publication in A&A

Journal ref: A&A 672, A162 (2023)

Cosmological Structure Formation and Soliton Phase Transition in Fuzzy Dark Matter with Axion Self-Interactions

Authors: Philip Mocz, Anastasia Fialkov, Mark Vogelsberger, Michael Boylan-Kolchin, Pierre-Henri Chavanis, Mustafa A. Amin, Sownak Bose, Tibor Dome, Lars Hernquist, Lachlan Lancaster, Matthew Notis, Connor Painter, Victor H. Robles, Jesus Zavala

Abstract: We investigate cosmological structure formation in Fuzzy Dark Matter (FDM) with an attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum 'pressure' and alter structure formation. We f… ▽ More We investigate cosmological structure formation in Fuzzy Dark Matter (FDM) with an attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum 'pressure' and alter structure formation. We find in our simulations that the SI can enhance small-scale structure formation, and soliton cores above a critical mass undergo a phase transition, transforming from dilute to dense solitons. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: 10 pages, 3 figures, submitted to mnras

Cosmic Web Dissection in Fuzzy Dark Matter Cosmologies

Authors: Tibor Dome, Anastasia Fialkov, Nina Sartorio, Philip Mocz

Abstract: On large cosmological scales, anisotropic gravitational collapse is manifest in the dark cosmic web. Its statistical properties are little known for alternative dark matter models such as fuzzy dark matter (FDM). In this work, we assess for the first time the relative importance of cosmic nodes, filaments, walls and voids in a cosmology with primordial small-scale suppression of power. We post-pro… ▽ More On large cosmological scales, anisotropic gravitational collapse is manifest in the dark cosmic web. Its statistical properties are little known for alternative dark matter models such as fuzzy dark matter (FDM). In this work, we assess for the first time the relative importance of cosmic nodes, filaments, walls and voids in a cosmology with primordial small-scale suppression of power. We post-process $N$-body simulations of FDM-like cosmologies with varying axion mass $m$ at redshifts $z\sim 1.0-5.6$ using the NEXUS+ Multiscale Morphology Filter technique at smoothing scale $Δx = 0.04 \ h^{-1}$Mpc. The formation of wall and void halos is more suppressed than naively expected from the half-mode mass $M_{1/2}$. Also, we quantify the mass and volume filling fraction of cosmic environments and find that 2D cosmic sheets host a larger share of the matter content of the Universe as $m$ is reduced, with an $\sim 8-12$\% increase for the $m=7 \times 10^{-22}$ eV model compared to CDM. We show that in FDM-like cosmologies, filaments, walls and voids are cleaner and more pronounced structures than in CDM, revealed by a strong mid-range peak in the conditioned overdensity PDFs $P(δ)$. At high redshift, low-density regions are more suppressed than high-density regions. Furthermore, skewness estimates $S_3$ of the total overdensity PDF in FDM-like cosmologies are consistently higher than in CDM, especially at high redshift $z\sim 5.6$ where the $m=10^{-22}$ eV model differs from CDM by $\sim 6 σ$. Accordingly, we advocate for the usage of $P(δ)$ as a testbed for constraining FDM and other alternative dark matter models. △ Less

Submitted 16 August, 2023; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: 16 pages, 7 figures

Journal ref: MNRAS, 525, 348-363 (2023)

Strong 21-cm fluctuations and anisotropy due to the line-of-sight effect of radio galaxies at cosmic dawn

Authors: Sudipta Sikder, Rennan Barkana, Anastasia Fialkov, Itamar Reis

Abstract: The reported detection of the global 21-cm signal by the EDGES collaboration is significantly stronger than standard astrophysical predictions. One possible explanation is an early radio excess above the cosmic microwave background. Such a radio background could have been produced by high redshift galaxies, if they were especially efficient in producing low-frequency synchrotron radiation. We have… ▽ More The reported detection of the global 21-cm signal by the EDGES collaboration is significantly stronger than standard astrophysical predictions. One possible explanation is an early radio excess above the cosmic microwave background. Such a radio background could have been produced by high redshift galaxies, if they were especially efficient in producing low-frequency synchrotron radiation. We have previously studied the effects of such an inhomogeneous radio background on the 21-cm signal; however, we made a simplifying assumption of isotropy of the background seen by each hydrogen cloud. Here we perform a complete calculation that accounts for the fact that the 21-cm absorption occurs along the line of sight, and is therefore sensitive to radio sources lying behind each absorbing cloud. We find that the complete calculation strongly enhances the 21-cm power spectrum during cosmic dawn, by up to two orders of magnitude; on the other hand, the effect on the global 21-cm signal is only at the $5\%$ level. In addition to making the high-redshift 21-cm fluctuations potentially more easily observable, the line of sight radio effect induces a new anisotropy in the 21-cm power spectrum. While these effects are particularly large for the case of an extremely-enhanced radio efficiency, they make it more feasible to detect even a moderately-enhanced radio efficiency in early galaxies. This is especially relevant since the EDGES signal has been contested by the SARAS experiment. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: Submitted to MNRAS

Joint analysis constraints on the physics of the first galaxies with low frequency radio astronomy data

Authors: Harry T. J. Bevins, Stefan Heimersheim, Irene Abril-Cabezas, Anastasia Fialkov, Eloy de Lera Acedo, William Handley, Saurabh Singh, Rennan Barkana

Abstract: The first billion years of cosmic history remains largely unobserved. We demonstrate, using a novel machine learning technique, how combining upper limits on the spatial fluctuations in the 21-cm signal with observations of the sky-averaged 21-cm signal from neutral hydrogen can improve our understanding of this epoch. By jointly analysing data from SARAS3 (redshift $z\approx15-25$) and limits fro… ▽ More The first billion years of cosmic history remains largely unobserved. We demonstrate, using a novel machine learning technique, how combining upper limits on the spatial fluctuations in the 21-cm signal with observations of the sky-averaged 21-cm signal from neutral hydrogen can improve our understanding of this epoch. By jointly analysing data from SARAS3 (redshift $z\approx15-25$) and limits from HERA ($z\approx8$ and $10$), we show that such a synergetic analysis provides tighter constraints on the astrophysics of galaxies 200 million years after the Big Bang than can be achieved with the individual data sets. Although our constraints are weak, this is the first time data from a sky-averaged 21-cm experiment and power spectrum experiment have been analysed together. In synergy, the two experiments leave only $64.9^{+0.3}_{-0.1}$% of the explored broad theoretical parameter space to be consistent with the joint data set, in comparison to $92.3^{+0.3}_{-0.1}$% for SARAS3 and $79.0^{+0.5}_{-0.2}$% for HERA alone. We use the joint analysis to constrain star formation efficiency, minimum halo mass for star formation, X-ray luminosity of early emitters and the radio luminosity of early galaxies. The joint analysis disfavours at 68% confidence a combination of galaxies with X-ray emission that is $\lesssim 33$ and radio emission that is $\gtrsim 32$ times as efficient as present day galaxies. We disfavour at $95$% confidence scenarios in which power spectra are $\geq126$ mK$^{2}$ at $z=25$ and the sky-averaged signals are $\leq-277$ mK. △ Less

Submitted 18 December, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: Accepted for MNRAS

Anticipating a New-Physics Signal in Upcoming 21-cm Power Spectrum Observations

Authors: Rennan Barkana, Anastasia Fialkov, Hongwan Liu, Nadav Joseph Outmezguine

Abstract: Dark matter-baryon interactions can cool the baryonic fluid, which has been shown to modify the cosmological 21-cm global signal. We show that in a two-component dark sector with an interacting millicharged component, dark matter-baryon scattering can produce a 21-cm power spectrum signal with acoustic oscillations. The signal can be up to three orders of magnitude larger than expected in $Λ$CDM c… ▽ More Dark matter-baryon interactions can cool the baryonic fluid, which has been shown to modify the cosmological 21-cm global signal. We show that in a two-component dark sector with an interacting millicharged component, dark matter-baryon scattering can produce a 21-cm power spectrum signal with acoustic oscillations. The signal can be up to three orders of magnitude larger than expected in $Λ$CDM cosmology, given realistic astrophysical models. This model provides a new-physics target for near-future experiments such as HERA or NenuFAR, which can potentially discover or strongly constrain the dark matter explanation of the putative EDGES anomaly. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Comments: 5 pages + appendices, 2 figures

Astrophysical Constraints from the SARAS3 non-detection of the Cosmic Dawn Sky-Averaged 21-cm Signal

Authors: H. T. J. Bevins, A. Fialkov, E. de Lera Acedo, W. J. Handley, S. Singh, R. Subrahmanyan, R. Barkana

Abstract: Observations of the redshifted 21-cm line of atomic hydrogen have resulted in several upper limits on the 21-cm power spectrum and a tentative detection of the sky-averaged signal at $z\sim17$. Made with the EDGES Low-Band antenna, this claim was recently disputed by the SARAS3 experiment, which reported a non-detection and is the only available upper limit strong enough to constrain cosmic dawn a… ▽ More Observations of the redshifted 21-cm line of atomic hydrogen have resulted in several upper limits on the 21-cm power spectrum and a tentative detection of the sky-averaged signal at $z\sim17$. Made with the EDGES Low-Band antenna, this claim was recently disputed by the SARAS3 experiment, which reported a non-detection and is the only available upper limit strong enough to constrain cosmic dawn astrophysics. We use these data to constrain a population of radio-luminous galaxies $\sim 200$ million years after the Big Bang ($z\approx 20$). We find, using Bayesian data analysis, that the data disfavours (at 68% confidence) radio-luminous galaxies in dark matter halos with masses of $4.4\times10^{5}$ M$_\odot \lesssim M \lesssim 1.1\times10^{7}$M$_\odot$ (where $M_\odot$ is the mass of the Sun) at $z = 20$ and galaxies in which $>5$% of the gas is converted into stars. The data disfavour galaxies with radio luminosity per star formation rate of $L_\mathrm{r}/\mathrm{SFR} \gtrsim 1.549 \times 10^{25}$ W Hz$^{-1}$M$_\odot^{-1}$ yr at 150 MHz, a thousand times brighter than today, and, separately, a synchrotron radio background in excess of the CMB by $\gtrsim 6%$ at 1.42 GHz. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: Published in Nature Astronomy

The Second Radio Synchrotron Background Workshop: Conference Summary and Report

Authors: J. Singal, N. Fornengo, M. Regis, G. Bernardi, D. Bordenave, E. Branchini, N. Cappelluti, A. Caputo, I. P. Carucci, J. Chluba, A. Cuoco, C. DiLullo, A. Fialkov, C. Hale, S. E. Harper, S. Heston, G. Holder, A. Kogut, M. G. H. Krause, J. P. Leahy, S. Mittal, R. A. Monsalve, G. Piccirilli, E. Pinetti, S. Recchia , et al. (2 additional authors not shown)

Abstract: We summarize the second radio synchrotron background workshop, which took place June 15-17, 2022 in Barolo, Italy. This meeting was convened because available measurements of the diffuse radio zero level continue to suggest that it is several times higher than can be attributed to known Galactic and extragalactic sources and processes, rendering it the least well understood electromagnetic backgro… ▽ More We summarize the second radio synchrotron background workshop, which took place June 15-17, 2022 in Barolo, Italy. This meeting was convened because available measurements of the diffuse radio zero level continue to suggest that it is several times higher than can be attributed to known Galactic and extragalactic sources and processes, rendering it the least well understood electromagnetic background at present and a major outstanding question in astrophysics. The workshop agreed on the next priorities for investigations of this phenomenon, which include searching for evidence of the Radio Sunyaev-Zel'dovich effect, carrying out cross-correlation analyses of radio emission with other tracers, and supporting the completion of the 310 MHz absolutely calibrated sky map project. △ Less

Submitted 1 March, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: 34 pages, 19 figures, 3 tables, accepted to PASP, updated to in press version

Journal ref: 2003, PASP, 135, 136001

Controlling ball progression in soccer

Authors: Catherine Pfaff, Emily Hunter, Haozhi Hong, Daniel Forestell, Ari Fialkov, Zoey Drassinower, Timothy Chan

Abstract: In this paper, we examine how soccer players can use their spatial relationships to control parts of the field and safely move play up the field via chains of ``safe configurations,'' i.e. configurations of players on a team ensuring the possessor of the ball has a collection of open passing options all connected by open passing lanes. An underlying philosophy behind our work is that it is most di… ▽ More In this paper, we examine how soccer players can use their spatial relationships to control parts of the field and safely move play up the field via chains of ``safe configurations,'' i.e. configurations of players on a team ensuring the possessor of the ball has a collection of open passing options all connected by open passing lanes. An underlying philosophy behind our work is that it is most difficult to disrupt an attacking team's progression forward (with the ball) when this attacking team has multiple ``good'' options of how to proceed at each moment in time. We provide some evidence of this. Our main construction is a directed weighted graph where the nodes encode the configurations of players, the directed edges encode transformations between these configurations, and the weights encode the relative frequencies of the transformations. We conclude with a few applications and proposed further investigations. We believe that our work can serve as a launching platform for significant further investigation into how teams can ``safely progress'' the ball up the field, strategy development, and sophisticated decision making metrics. For coaches and players, we aim to streamline the process of moving safely up the field. In particular, we aim to construct a framework for creating new successful patterns of movement and aim that our framework allows for the development of new strategy. At the same time, our work could help identify configurations on the field which often result in turnovers or that allow for a lot of strategic flexibility (also impeding defensive containment of the attacking team). *Because the contributions of the female authors to this paper were by no means less than those of the male authors, we have chosen to reverse convention and to list the authors in reverse alphabetical order to ensure that the female authors were not listed only after the male authors. △ Less

Submitted 17 November, 2022; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: correction to bug in code (leading to new figures and example) and to a few references

The REACH radiometer for detecting the 21-cm hydrogen signal from redshift 7.5 to 28

Authors: E. de Lera Acedo, D. I. L. de Villiers, N. Razavi-Ghods, W. Handley, A. Fialkov, A. Magro, D. Anstey, H. T. J. Bevins, R. Chiello, J. Cumner, A. T. Josaitis, I. L. V. Roque, P. H. Sims, K. H. Scheutwinkel, P. Alexander, G. Bernardi, S. Carey, J. Cavillot, W. Croukamp, J. A. Ely, T. Gessey-Jones, Q. Gueuning, R. Hills, G. Kulkarni, R. Maiolino , et al. (9 additional authors not shown)

Abstract: Observations of the 21-cm line from primordial hydrogen promise to be one of the best tools to study the early epochs of the Universe: the Dark Ages, the Cosmic Dawn, and the subsequent Epoch of Reionization. In 2018, the EDGES experiment caught the attention of the cosmology community with a potential detection of an absorption feature in the sky-averaged radio spectrum centred at 78 MHz. The fea… ▽ More Observations of the 21-cm line from primordial hydrogen promise to be one of the best tools to study the early epochs of the Universe: the Dark Ages, the Cosmic Dawn, and the subsequent Epoch of Reionization. In 2018, the EDGES experiment caught the attention of the cosmology community with a potential detection of an absorption feature in the sky-averaged radio spectrum centred at 78 MHz. The feature is deeper than expected, and, if confirmed, would call for new physics. However, different groups have re-analyzed the EDGES data and questioned the reliability of the signal. The Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) is a sky-averaged 21-cm experiment aiming at improving the current observations by tackling the issues faced by current instruments related to residual systematic signals in the data. The novel experimental approach focuses on detecting and jointly explaining these systematics together with the foregrounds and the cosmological signal using Bayesian statistics. To achieve this, REACH features simultaneous observations with two different antennas, an ultra wideband system (redshift range 7.5 to 28), and a receiver calibrator based on in-field measurements. Simulated observations forecast percent-level constraints on astrophysical parameters, potentially opening up a new window to the infant Universe. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: Publication date: 21 July 2022

Journal ref: Nature Astronomy 6, p 984 998, 2022

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

On the Cosmic Web Elongation in Fuzzy Dark Matter Cosmologies: Effects on Density Profiles, Shapes and Alignments of Halos

Authors: Tibor Dome, Anastasia Fialkov, Philip Mocz, Björn Malte Schäfer, Michael Boylan-Kolchin, Mark Vogelsberger

Abstract: The fuzzy dark matter (FDM) scenario has received increased attention in recent years due to the small-scale challenges of the vanilla Lambda cold dark matter ($Λ$CDM) cosmological model and the lack of any experimental evidence for any candidate particle. In this study, we use cosmological $N$-body simulations to investigate high-redshift dark matter halos and their responsiveness to an FDM-like… ▽ More The fuzzy dark matter (FDM) scenario has received increased attention in recent years due to the small-scale challenges of the vanilla Lambda cold dark matter ($Λ$CDM) cosmological model and the lack of any experimental evidence for any candidate particle. In this study, we use cosmological $N$-body simulations to investigate high-redshift dark matter halos and their responsiveness to an FDM-like power spectrum cutoff on small scales in the primordial density perturbations. We study halo density profiles, shapes and alignments in FDM-like cosmologies (the latter two for the first time) by providing fits and quantifying departures from $Λ$CDM as a function of the particle mass $m$. Compared to $Λ$CDM, the concentrations of FDM-like halos are lower, peaking at an $m$-dependent halo mass and thus breaking the approximate universality of density profiles in $Λ$CDM. The intermediate-to-major and minor-to-major shape parameter profiles are monotonically increasing with ellipsoidal radius in $N$-body simulations of $Λ$CDM. In FDM-like cosmologies, the monotonicity is broken, halos are more elongated around the virial radius than their $Λ$CDM counterparts and less elongated closer to the center. Finally, intrinsic alignment correlations, stemming from the deformation of initially spherically collapsing halos in an ambient gravitational tidal field, become stronger with decreasing $m$. At $z\sim 4$, we find a $6.4 σ$-significance in the fractional differences between the isotropised linear alignment magnitudes $D_{\text{iso}}$ in the $m=10^{-22}$ eV model and $Λ$CDM. Such FDM-like imprints on the internal properties of virialised halos are expected to be strikingly visible in the high-$z$ Universe. △ Less

Submitted 19 January, 2023; v1 submitted 7 August, 2022; originally announced August 2022.

Comments: 20 pages, 10 figures, 3 tables

Journal ref: MNRAS, 519, 4183-4202 (2022)

Relative baryon-dark matter velocities in cosmological zoom simulations

Authors: Luke Conaboy, Ilian T. Iliev, Anastasia Fialkov, Keri L. Dixon, David Sullivan

Abstract: Supersonic relative motion between baryons and dark matter due to the decoupling of baryons from the primordial plasma after recombination affects the growth of the first small-scale structures. Large box sizes (greater than a few hundred Mpc) are required to sample the full range of scales pertinent to the relative velocity, while the effect of the relative velocity is strongest on small scales (… ▽ More Supersonic relative motion between baryons and dark matter due to the decoupling of baryons from the primordial plasma after recombination affects the growth of the first small-scale structures. Large box sizes (greater than a few hundred Mpc) are required to sample the full range of scales pertinent to the relative velocity, while the effect of the relative velocity is strongest on small scales (less than a few hundred kpc). This separation of scales naturally lends itself to the use of `zoom' simulations, and here we present our methodology to self-consistently incorporate the relative velocity in zoom simulations, including its cumulative effect from recombination through to the start time of the simulation. We apply our methodology to a large-scale cosmological zoom simulation, finding that the inclusion of relative velocities suppresses the halo baryon fraction by $46$--$23$ per cent between $z=13.6$ and $11.2$, in qualitative agreement with previous works. In addition, we find that including the relative velocity delays the formation of star particles by $\sim 20 {~\rm Myr}$ Myr on average (of the order of the lifetime of a $\sim 9~{\rm M}_\odot$ Population III star) and suppresses the final stellar mass by as much as $79$ per cent at $z=11.2$. △ Less

Submitted 14 September, 2023; v1 submitted 23 July, 2022; originally announced July 2022.

Comments: 14 pages, 12 figures. Accepted for publication in MNRAS

Marginal Bayesian Statistics Using Masked Autoregressive Flows and Kernel Density Estimators with Examples in Cosmology

Authors: Harry Bevins, Will Handley, Pablo Lemos, Peter Sims, Eloy de Lera Acedo, Anastasia Fialkov

Abstract: Cosmological experiments often employ Bayesian workflows to derive constraints on cosmological and astrophysical parameters from their data. It has been shown that these constraints can be combined across different probes such as Planck and the Dark Energy Survey and that this can be a valuable exercise to improve our understanding of the universe and quantify tension between multiple experiments.… ▽ More Cosmological experiments often employ Bayesian workflows to derive constraints on cosmological and astrophysical parameters from their data. It has been shown that these constraints can be combined across different probes such as Planck and the Dark Energy Survey and that this can be a valuable exercise to improve our understanding of the universe and quantify tension between multiple experiments. However, these experiments are typically plagued by differing systematics, instrumental effects and contaminating signals, which we collectively refer to as `nuisance' components, that have to be modelled alongside target signals of interest. This leads to high dimensional parameter spaces, especially when combining data sets, with > 20 dimensions of which only around 5 correspond to key physical quantities. We present a means by which to combine constraints from different data sets in a computationally efficient manner by generating rapid, reusable and reliable marginal probability density estimators, giving us access to nuisance-free likelihoods. This is possible through the unique combination of nested sampling, which gives us access to Bayesian evidences, and the marginal Bayesian statistics code MARGARINE. Our method is lossless in the signal parameters, resulting in the same posterior distributions as would be found from a full nested sampling run over all nuisance parameters, and typically quicker than evaluating full likelihoods. We demonstrate our approach by applying it to the combination of posteriors from the Dark Energy Survey and Planck. △ Less

Submitted 24 November, 2022; v1 submitted 23 July, 2022; originally announced July 2022.

Comments: Published in Phys. Sci. Forum 2022, 5(1), 1; https://doi.org/10.3390/psf2022005001

Antenna beam characterisation for the global 21cm experiment LEDA and its impact on signal model parameter reconstruction

Authors: M. Spinelli, G. Kyriakou, G. Bernardi, P. Bolli, L. J. Greenhill, A. Fialkov, H. Garsden

Abstract: Cosmic Dawn, the onset of star formation in the early universe, can in principle be studied via the 21cm transition of neutral hydrogen, for which a sky-averaged absorption signal, redshifted to MHz frequencies, is predicted to be {\it O}(10-100)\,mK. Detection requires separation of the 21cm signal from bright chromatic foreground emission due to Galactic structure, and the characterisation of ho… ▽ More Cosmic Dawn, the onset of star formation in the early universe, can in principle be studied via the 21cm transition of neutral hydrogen, for which a sky-averaged absorption signal, redshifted to MHz frequencies, is predicted to be {\it O}(10-100)\,mK. Detection requires separation of the 21cm signal from bright chromatic foreground emission due to Galactic structure, and the characterisation of how it couples to instrumental response. In this work, we present characterisation of antenna gain patterns for the Large-aperture Experiment to detect the Dark Ages (LEDA) via simulations, assessing the effects of the antenna ground-plane geometries used, and measured soil properties. We then investigate the impact of beam pattern uncertainties on the reconstruction of a Gaussian absorption feature. Assuming the pattern is known and correcting for the chromaticity of the instrument, the foregrounds can be modelled with a log-polynomial, and the 21cm signal identified with high accuracy. However, uncertainties on the soil properties lead to \textperthousand\ changes in the chromaticity that can bias the signal recovery. The bias can be up to a factor of two in amplitude and up to few \% in the frequency location. These effects do not appear to be mitigated by larger ground planes, conversely gain patterns with larger ground planes exhibit more complex frequency structure, significantly compromising the parameter reconstruction. Our results, consistent with findings from other antenna design studies, emphasise the importance of chromatic response and suggest caution in assuming log-polynomial foreground models in global signal experiments. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS

Marginal Post Processing of Bayesian Inference Products with Normalizing Flows and Kernel Density Estimators

Authors: Harry T. J. Bevins, William J. Handley, Pablo Lemos, Peter H. Sims, Eloy de Lera Acedo, Anastasia Fialkov, Justin Alsing

Abstract: Bayesian analysis has become an indispensable tool across many different cosmological fields including the study of gravitational waves, the Cosmic Microwave Background and the 21-cm signal from the Cosmic Dawn among other phenomena. The method provides a way to fit complex models to data describing key cosmological and astrophysical signals and a whole host of contaminating signals and instrument… ▽ More Bayesian analysis has become an indispensable tool across many different cosmological fields including the study of gravitational waves, the Cosmic Microwave Background and the 21-cm signal from the Cosmic Dawn among other phenomena. The method provides a way to fit complex models to data describing key cosmological and astrophysical signals and a whole host of contaminating signals and instrumental effects modelled with `nuisance parameters'. In this paper, we summarise a method that uses Masked Autoregressive Flows and Kernel Density Estimators to learn marginal posterior densities corresponding to core science parameters. We find that the marginal or 'nuisance-free' posteriors and the associated likelihoods have an abundance of applications including; the calculation of previously intractable marginal Kullback-Leibler divergences and marginal Bayesian Model Dimensionalities, likelihood emulation and prior emulation. We demonstrate each application using toy examples, examples from the field of 21-cm cosmology and samples from the Dark Energy Survey. We discuss how marginal summary statistics like the Kullback-Leibler divergences and Bayesian Model Dimensionalities can be used to examine the constraining power of different experiments and how we can perform efficient joint analysis by taking advantage of marginal prior and likelihood emulators. We package our multipurpose code up in the pip-installable code margarine for use in the wider scientific community. △ Less

Submitted 18 December, 2023; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: Accepted for MNRAS

Bayesian Data Analysis for Sky-averaged 21-cm Experiments in the Presence of Ionospheric Effects

Authors: Emma Shen, Dominic Anstey, Eloy de Lera Acedo, Anastasia Fialkov

Abstract: The ionosphere introduces chromatic distortions on low frequency radio waves, and thus poses a hurdle for 21-cm cosmology. In this paper we introduce time-varying chromatic ionospheric effects on simulated antenna temperature data of a global 21-cm data analysis pipeline, and try to detect the injected global signal. We demonstrate that given turbulent ionospheric conditions, more than 5\% error i… ▽ More The ionosphere introduces chromatic distortions on low frequency radio waves, and thus poses a hurdle for 21-cm cosmology. In this paper we introduce time-varying chromatic ionospheric effects on simulated antenna temperature data of a global 21-cm data analysis pipeline, and try to detect the injected global signal. We demonstrate that given turbulent ionospheric conditions, more than 5\% error in our knowledge of the ionospheric parameters could lead to comparatively low evidence and high root-mean-square error (RMSE), suggesting a false or null detection. When using a constant antenna beam for cases that include data at different times, the significance of the detection lowers as the number of time samples increases. It is also shown that for observations that include data at different times, readjusting beam configurations according to the time-varying ionospheric conditions should greatly improve the significance of a detection, yielding higher evidences and lower RMSE, and that it is a necessary procedure for a successful detection when the ionospheric conditions are not ideal. △ Less

Submitted 22 April, 2022; originally announced April 2022.

Comments: 8 pages, 7 figures, 7 tables

Impact of the Primordial Stellar Initial Mass Function on the 21-cm Signal

Authors: T. Gessey-Jones, N. S. Sartorio, A. Fialkov, G. M. Mirouh, M. Magg, R. G. Izzard, E. de Lera Acedo, W. J. Handley, R. Barkana

Abstract: Properties of the first generation of stars (Pop III), such as their initial mass function (IMF), are poorly constrained by observations and have yet to converge between simulations. The cosmological 21-cm signal of neutral hydrogen is predicted to be sensitive to Lyman-band photons produced by these stars, thus providing a unique way to probe the first stellar population. In this paper, we invest… ▽ More Properties of the first generation of stars (Pop III), such as their initial mass function (IMF), are poorly constrained by observations and have yet to converge between simulations. The cosmological 21-cm signal of neutral hydrogen is predicted to be sensitive to Lyman-band photons produced by these stars, thus providing a unique way to probe the first stellar population. In this paper, we investigate the impacts of the Pop III IMF on the cosmic dawn 21-cm signal via the Wouthuysen-Field effect, Lyman-Werner feedback, Ly-alpha heating, and CMB heating. We calculate the emission spectra of star-forming halos for different IMFs by integrating over individual metal-free stellar spectra, computed from a set of stellar evolution histories and stellar atmospheres, and taking into account variability of the spectra with stellar age. Through this study, we therefore relax two common assumptions: that the zero-age main sequence emission rate of a Pop III star is representative of its lifetime mean emission rate, and that Pop III emission can be treated as instantaneous. Exploring a bottom-heavy, a top-heavy, and intermediate IMFs, we show that variations in the 21-cm signal are driven by stars lighter than 20 solar masses. For the explored models we find maximum relative differences of 59% in the cosmic dawn global 21-cm signal, and 131% between power spectra. Although this impact is modest, precise modelling of the first stars and their evolution is necessary for accurate prediction and interpretation of the 21-cm signal. △ Less

Submitted 21 July, 2022; v1 submitted 4 February, 2022; originally announced February 2022.

Comments: 20 pages, 16 figures. Accepted for publication in MNRAS, updated to accepted version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 516, Issue 1, October 2022, Pages 841-860

A Comprehensive Bayesian re-analysis of the SARAS2 data from the Epoch of Reionization

Authors: H. T. J. Bevins, E. de Lera Acedo, A. Fialkov, W. J. Handley, S. Singh, R. Subrahmanyan, R. Barkana

Abstract: We present a Bayesian re-analysis of the sky-averaged 21-cm experimental data from SARAS2 using nested sampling implemented with polychord, spectrally smooth foreground modelling implemented with maxsmooth, detailed systematic modelling and rapid signal emulation with globalemu. Our analysis differs from previous analysis of the SARAS2 data through the use of a full Bayesian framework and separate… ▽ More We present a Bayesian re-analysis of the sky-averaged 21-cm experimental data from SARAS2 using nested sampling implemented with polychord, spectrally smooth foreground modelling implemented with maxsmooth, detailed systematic modelling and rapid signal emulation with globalemu. Our analysis differs from previous analysis of the SARAS2 data through the use of a full Bayesian framework and separate modelling of the foreground and non-smooth systematics. We use the most up-to-date signal models including Lyman-$α$ and CMB heating parameterised by astrophysical parameters such as star formation efficiency, X-ray heating efficiency, minimal virial circular velocity of star forming galaxies, CMB optical depth and the low energy cutoff of the X-ray spectral energy distribution. We consider models with an excess radio background above the CMB produced via radio emission from early galaxies and parameterised by a radio production efficiency. A non-smooth systematic is identified and modelled as both a frequency damped sinusoid introduced by the electronics and separately from the sky. The latter is modulated by the total efficiency of the antenna and marginally favoured by the data. We consider three different models for the noise in the data. The SARAS2 constraints on individual astrophysical parameters are extremely weak however we identify classes of disfavoured signals. We weakly disfavour standard astrophysical models with high Lyman-$α$ fluxes and weak heating and more confidently disfavour exotic models with high Lyman-$α$ fluxes, low X-ray efficiencies and high radio production efficiencies in early galaxies. △ Less

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

Comments: 19 Pages, 13 Figures. Accepted by MNRAS

Emulation of the Cosmic Dawn 21-cm Power Spectrum and Classification of Excess Radio Models Using an Artificial Neural Network

Authors: Sudipta Sikder, Rennan Barkana, Itamar Reis, Anastasia Fialkov

Abstract: The cosmic 21-cm line of hydrogen is expected to be measured in detail by the next generation of radio telescopes. The enormous dataset from future 21-cm surveys will revolutionize our understanding of early cosmic times. We present a machine learning approach based on an Artificial Neural Network that uses emulation in order to uncover the astrophysics in the epoch of reionization and cosmic dawn… ▽ More The cosmic 21-cm line of hydrogen is expected to be measured in detail by the next generation of radio telescopes. The enormous dataset from future 21-cm surveys will revolutionize our understanding of early cosmic times. We present a machine learning approach based on an Artificial Neural Network that uses emulation in order to uncover the astrophysics in the epoch of reionization and cosmic dawn. Using a seven-parameter astrophysical model that covers a very wide range of possible 21-cm signals, over the redshift range 6 to 30 and wavenumber range $0.05$ to $1 \ \rm{Mpc}^{-1}$ we emulate the 21-cm power spectrum with a typical accuracy of $10 - 20\%$. As a realistic example, we train an emulator using the power spectrum with an optimistic noise model of the Square Kilometre Array (SKA). Fitting to mock SKA data results in a typical measurement accuracy of $2.8\%$ in the optical depth to the cosmic microwave background, $34\%$ in the star-formation efficiency of galactic halos, and a factor of 9.6 in the X-ray efficiency of galactic halos. Also, with our modeling we reconstruct the true 21-cm power spectrum from the mock SKA data with a typical accuracy of $15 - 30\%$. In addition to standard astrophysical models, we consider two exotic possibilities of strong excess radio backgrounds at high redshifts. We use a neural network to identify the type of radio background present in the 21-cm power spectrum, with an accuracy of $87\%$ for mock SKA data. △ Less

Submitted 10 January, 2024; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: Revision to match MNRAS published version

Effect of the cosmological transition to metal-enriched star-formation on the hydrogen 21-cm signal

Authors: Mattis Magg, Itamar Reis, Anastasia Fialkov, Rennan Barkana, Ralf S. Klessen, Simon C. O. Glover, Li-Hsin Chen, Tilman Hartwig, Anna T. P. Schauer

Abstract: Mapping Cosmic Dawn with 21-cm tomography offers an exciting new window into the era of primordial star formation. However, self-consistent implementation of both the process of star formation and the related 21-cm signal is challenging, due to the multi-scale nature of the problem. In this study, we develop a flexible semi-analytical model to follow the formation of the first stars and the proces… ▽ More Mapping Cosmic Dawn with 21-cm tomography offers an exciting new window into the era of primordial star formation. However, self-consistent implementation of both the process of star formation and the related 21-cm signal is challenging, due to the multi-scale nature of the problem. In this study, we develop a flexible semi-analytical model to follow the formation of the first stars and the process of gradual transition from primordial to metal-enriched star formation. For this transition we use different in scenarios with varying time-delays (or recovery times) between the first supernovae and the formation of the second generation of stars. We use recovery times between 10 and 100\,Myr and find that these delays have a strong impact on the redshift at which the transition to metal-enriched star formation occurs. We then explore the effect of this transition on the 21-cm signal and find that the recovery time has a distinctive imprint in the signal. Together with an improved understanding of how this time-delay relates to the properties of Population~III stars, future 21-cm observations can give independent constraints on the earliest epoch of star formation. △ Less

Submitted 29 October, 2021; originally announced October 2021.

Comments: 17 pages, 11 Figures, submitted to MNRAS, comments very welcome

EuCAPT White Paper: Opportunities and Challenges for Theoretical Astroparticle Physics in the Next Decade

Authors: R. Alves Batista, M. A. Amin, G. Barenboim, N. Bartolo, D. Baumann, A. Bauswein, E. Bellini, D. Benisty, G. Bertone, P. Blasi, C. G. Böhmer, Ž. Bošnjak, T. Bringmann, C. Burrage, M. Bustamante, J. Calderón Bustillo, C. T. Byrnes, F. Calore, R. Catena, D. G. Cerdeño, S. S. Cerri, M. Chianese, K. Clough, A. Cole, P. Coloma , et al. (112 additional authors not shown)

Abstract: Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, und… ▽ More Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, under the coordination of the European Consortium for Astroparticle Theory (EuCAPT). Addressed to the whole astroparticle physics community, it explores upcoming theoretical opportunities and challenges for our field of research, with particular emphasis on the possible synergies among different subfields, and the prospects for solving the most fundamental open questions with multi-messenger observations. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: White paper of the European Consortium for Astroparticle Theory (EuCAPT). 135 authors, 400 endorsers, 133 pages, 1382 references

Radio antenna design for sky-averaged 21 cm cosmology experiments: the REACH case

Authors: J. Cumner, E. De Lera Acedo, D. I. L. de Villiers, D. Anstey, C. I. Kolitsidas, B. Gurdon, N. Fagnoni, P. Alexander, G. Bernardi, H. T. J. Bevins, S. Carey, J. Cavillot, R. Chiello, C. Craeye, W. Croukamp, J. A. Ely, A. Fialkov, T. Gessey-Jones, Q. Gueuning, W. Handley, R. Hills, A. T. Josaitis, G. Kulkarni, A. Magro, R. Maiolino , et al. (13 additional authors not shown)

Abstract: Following the reported detection of an absorption profile associated with the 21~cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21~cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experim… ▽ More Following the reported detection of an absorption profile associated with the 21~cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21~cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experiment will seek to understand and compensate for systematic errors present using detailed modelling and characterization of the instrumentation. There is detailed the quantitative figures of merit and numerical modelling used to assist the design process of the REACH dipole antenna (one of the 2 antenna designs for REACH Phase I). This design process produced a 2.5:1 frequency bandwidth dipole. The aim of this design was to balance spectral smoothness and low impedance reflections with the ability to describe and understand the antenna response to the sky signal to inform the critically important calibration during observation and data analysis. △ Less

Submitted 12 January, 2023; v1 submitted 21 September, 2021; originally announced September 2021.

Comments: 32 pages, 30 figures, to be submitted to the Journal of Astronomical Instrumentation

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

What it Takes to Measure Reionization with Fast Radio Bursts

Authors: Stefan Heimersheim, Nina S. Sartorio, Anastasia Fialkov, Duncan R. Lorimer

Abstract: Fast Radio Bursts (FRBs) are extra-galactic radio transients which exhibit a distance-dependent dispersion of their signal, and thus can be used as cosmological probes. In this article we, for the first time, apply a model-independent approach to measure reionization from synthetic FRB data assuming these signals are detected beyond redshift 5. This method allows us to constrain the full shape of… ▽ More Fast Radio Bursts (FRBs) are extra-galactic radio transients which exhibit a distance-dependent dispersion of their signal, and thus can be used as cosmological probes. In this article we, for the first time, apply a model-independent approach to measure reionization from synthetic FRB data assuming these signals are detected beyond redshift 5. This method allows us to constrain the full shape of the reionization history as well as the CMB optical depth $τ$ while avoiding the problems of commonly used model-based techniques. 100 localized FRBs, originating from redshifts 5-15, could constrain (at 68% confidence level) the CMB optical depth to within 11%, and the midpoint of reionization to 4%, surpassing current state-of-the-art CMB bounds and quasar limits. Owing to the higher numbers of expected FRBs at lower redshifts, the $τ$ constraints are asymmetric (+14%, -7%) providing a much stronger lower limit. Finally, we show that the independent constraints on reionization from FRBs will improve limits on other cosmological parameters such as the amplitude of the power spectrum of primordial fluctuations. △ Less

Submitted 18 May, 2022; v1 submitted 29 July, 2021; originally announced July 2021.

Comments: 13 pages, 5 figures. Accepted by ApJ. Code and data available at https://zenodo.org/record/6542596

21cmVAE: A Very Accurate Emulator of the 21-cm Global Signal

Authors: Christian H. Bye, Stephen K. N. Portillo, Anastasia Fialkov

Abstract: Considerable observational efforts are being dedicated to measuring the sky-averaged (global) 21-cm signal of neutral hydrogen from Cosmic Dawn and the Epoch of Reionization. Deriving observational constraints on the astrophysics of this era requires modeling tools that can quickly and accurately generate theoretical signals across the wide astrophysical parameter space. For this purpose artificia… ▽ More Considerable observational efforts are being dedicated to measuring the sky-averaged (global) 21-cm signal of neutral hydrogen from Cosmic Dawn and the Epoch of Reionization. Deriving observational constraints on the astrophysics of this era requires modeling tools that can quickly and accurately generate theoretical signals across the wide astrophysical parameter space. For this purpose artificial neural networks were used to create the only two existing global signal emulators, 21cmGEM and globalemu. In this paper we introduce 21cmVAE, a neural network-based global signal emulator, trained on the same dataset of ~30,000 global signals as the other two emulators, but with a more direct prediction algorithm that prioritizes accuracy and simplicity. Using neural networks, we compute derivatives of the signals with respect to the astrophysical parameters and establish the most important astrophysical processes that drive the global 21-cm signal at different epochs. 21cmVAE has a relative rms error of only 0.34 - equivalently 0.54 mK - on average, which is a significant improvement compared to the existing emulators, and a run time of 0.04 seconds per parameter set. The emulator, the code, and the processed datasets are publicly available at https://github.com/christianhbye/21cmVAE and through https://zenodo.org/record/5904939. △ Less

Submitted 16 April, 2022; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: 11 pages, 7 figures, 2 tables (including appendices). In press (ApJ). Code available at https://github.com/christianhbye/21cmVAE

Shot noise and scatter in the star formation efficiency as a source of 21-cm fluctuations

Authors: Itamar Reis, Rennan Barkana, Anastasia Fialkov

Abstract: The 21-cm signal from cosmic dawn and the epoch of reionization (EoR) probes the characteristics of the high redshift galaxy population. Many of the astrophysical properties of galaxies at high redshifts are currently unconstrained due to the lack of observations. This creates a vast space of possible astrophysical scenarios where the 21-cm signal needs to be modeled in order to plan for, and even… ▽ More The 21-cm signal from cosmic dawn and the epoch of reionization (EoR) probes the characteristics of the high redshift galaxy population. Many of the astrophysical properties of galaxies at high redshifts are currently unconstrained due to the lack of observations. This creates a vast space of possible astrophysical scenarios where the 21-cm signal needs to be modeled in order to plan for, and eventually fit, future observations. This is done with fast numerical methods which make simplifying approximations for the underlying physical processes. In this work we quantify the effect of Poisson fluctuations and scatter in the star formation efficiency; while Poisson fluctuations are included in some works and not in others, scatter in the star formation efficiency is usually neglected, and all galaxies of a given mass are assumed to have the same properties. We show that both features can have a significant effect on the 21-cm power spectrum, most importantly in scenarios where the signal is dominated by massive galaxies. Scatter in the star formation efficiency does not simply enhance the effect of Poisson fluctuations; for example we show that the power spectrum shape at cosmic dawn has a feature corresponding to the width of the galaxy brightness distribution. We also discuss some of the consequences for 21-cm imaging, and the signature of reduced correlation between the density and radiation fields. △ Less

Submitted 24 June, 2021; originally announced June 2021.

Comments: Submitted to MNRAS on June 7th

GLOBALEMU: A novel and robust approach for emulating the sky-averaged 21-cm signal from the cosmic dawn and epoch of reionisation

Authors: H. T. J. Bevins, W. J. Handley, A. Fialkov, E. de Lera Acedo, K. Javid

Abstract: Emulation of the Global (sky-averaged) 21-cm signal with neural networks has been shown to be an essential tool for physical signal modelling. In this paper we present globalemu, a Global 21-cm signal emulator that uses redshift as a character defining variable alongside a set of astrophysical parameters to estimate the signal brightness temperature. Combined with physically-motivated data pre-pro… ▽ More Emulation of the Global (sky-averaged) 21-cm signal with neural networks has been shown to be an essential tool for physical signal modelling. In this paper we present globalemu, a Global 21-cm signal emulator that uses redshift as a character defining variable alongside a set of astrophysical parameters to estimate the signal brightness temperature. Combined with physically-motivated data pre-processing this makes for a reliable and fast emulator that is relatively insensitive to the network design. globalemu can emulate a high resolution signal in 1.3 ms in comparison to 133 ms, a factor of 102 improvement, when using the existing public state of the art 21cmGEM. We illustrate, with the standard astrophysical models used to train 21cmGEM, that globalemu is almost twice as accurate and for a test set of $\approx1,700$ signals we achieve a mean RMSE of 2.52 mK across the band z=7-28 ($\approx$10 per cent the expected noise of the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH)). The models are parameterised by the star formation efficiency, $f_*$, minimum virial circular velocity, $V_c$, X-ray efficiency, $f_X$, CMB optical depth, $τ$, the slope and low energy cut off of the X-ray spectral energy density, $α$ and $ν_\mathrm{min}$, and the mean free path of ionizing photons, $R_\mathrm{mfp}$. globalemu provides a flexible framework for easily emulating updated simulations of the Global signal and in addition the neutral fraction history. The emulator is pip installable and available at: \url{https://github.com/htjb/globalemu}. globalemu will be used extensively by the REACH collaboration. △ Less

Submitted 15 October, 2021; v1 submitted 9 April, 2021; originally announced April 2021.

Comments: 14 pages, 14 figures

A 21-cm power spectrum at 48 MHz, using the Owens Valley Long Wavelength Array

Authors: Hugh Garsden, Lincoln Greenhill, Gianni Bernardi, Anastasia Fialkov, Daniel C. Price, Daniel Mitchell, Jayce Dowell, Marta Spinelli, Frank K. Schinzel

Abstract: The Large-aperture Experiment to detect the Dark Age (LEDA) was designed to measure the 21-cm signal from neutral hydrogen at Cosmic Dawn, $z \approx $15-30. Using observations made with the $\approx $ 200 m diameter core of the Owens Valley Long Wavelength Array (OVRO-LWA), we present a 2-D cylindrical spatial power spectrum for data at 43.1-53.5 MHz ($z_{\rm median}\approx 28$) incoherently inte… ▽ More The Large-aperture Experiment to detect the Dark Age (LEDA) was designed to measure the 21-cm signal from neutral hydrogen at Cosmic Dawn, $z \approx $15-30. Using observations made with the $\approx $ 200 m diameter core of the Owens Valley Long Wavelength Array (OVRO-LWA), we present a 2-D cylindrical spatial power spectrum for data at 43.1-53.5 MHz ($z_{\rm median}\approx 28$) incoherently integrated for 4 hours, and an analysis of the array sensitivity. Power from foregrounds is localized to a "wedge" within $k_\perp, k_\parallel$ space. After calibration of visibilities using 5 bright compact sources including VirA, we measure $Δ^2(k) \approx 2 \times 10^{12}\ \mathrm{mK}^2$ outside the foreground wedge, where an uncontaminated cosmological signal would lie, in principle. The measured $Δ^2(k)$ is an upper limit that reflects a combination of thermal instrumental and sky noise, and unmodelled systematics that scatter power from the wedge, as will be discussed. By differencing calibrated visibilities for close pairs of frequency channels, we suppress foreground sky structure and systematics, extract thermal noise, and use a mix of coherent and incoherent integration to simulate a noise-dominated power spectrum for a 3000 h observation and $z = $16-37. For suitable calibration quality, the resulting noise level, $Δ^2(k) \approx 100$ mK$^2$ (k = 0.3 Mpc$^{-1}$), would be sufficient to detect peaks in the 21-cm spatial power spectrum due to early Ly-$α$ and X-ray sources, as predicted for a range of theoretical model parameters. △ Less

Submitted 22 June, 2021; v1 submitted 18 February, 2021; originally announced February 2021.

Comments: 16 pages, 11 figures. Accepted for MNRAS; replaced with accepted version

The subtlety of Ly-a photons: changing the expected range of the 21-cm signal

Authors: Itamar Reis, Anastasia Fialkov, Rennan Barkana

Abstract: We present the evolution of the 21-cm signal from cosmic dawn and the epoch of reionization (EoR) in an upgraded model including three subtle effects of Ly-a radiation: Ly-a heating, CMB heating (mediated by Ly-a photons), and multiple scattering of Ly-a photons. Taking these effects into account we explore a wide range of astrophysical models and quantify the impact of these processes on the glob… ▽ More We present the evolution of the 21-cm signal from cosmic dawn and the epoch of reionization (EoR) in an upgraded model including three subtle effects of Ly-a radiation: Ly-a heating, CMB heating (mediated by Ly-a photons), and multiple scattering of Ly-a photons. Taking these effects into account we explore a wide range of astrophysical models and quantify the impact of these processes on the global 21-cm signal and its power spectrum at observable scales and redshifts. We find that, in agreement with the literature, Ly-a and CMB heating raise the gas temperature by up to $\mathcal{O}(100)$ degrees in models with weak X-ray heating and, thus, suppress the predicted 21-cm signals. Varying the astrophysical parameters over broad ranges, we find that in the upgraded model the absorption trough of the global signal reaches a lowest floor of $-165$ mK at redshifts $z\approx 15-19$. This is in contrast with the predictions for a pure adiabatically cooling Universe, for which the deepest possible absorption is a monotonically decreasing function of cosmic time and is $-178$ mK at $z = 19$ and $-216$ mK at $z=15$, dropping to even lower values at lower redshifts (e.g. $-264$ mK at $z = 10$). With the Ly-a and CMB heating included we also observe a strong suppression in the low-redshift power spectra, with the maximum possible power (evaluated over the ensemble of models) attenuated by a factor of $6.6$ at $z=9$ and $k = 0.1$ Mpc$^{-1}$. Finally, we find that at high redshifts corresponding to cosmic dawn, the heating terms have a subdominant effect while multiple scattering of Ly-a photons is important, leading to an amplification of the power spectrum by a factor of $\sim 2-5$. △ Less

Submitted 12 January, 2021; v1 submitted 5 January, 2021; originally announced January 2021.

Quantifying Ionospheric Effects on Global 21-cm Observations

Authors: Emma Shen, Dominic Anstey, Eloy de Lera Acedo, Anastasia Fialkov, Will Handley

Abstract: We modelled the two major layer of Earth's ionosphere, the F-layer and the D-layer, by a simplified spatial model with temporal variance to study the chromatic ionospheric effects on global 21-cm observations. From the analyses, we found that the magnitude of the ionospheric disruptions due to ionospheric refraction and absorption can be greater than the expected global 21-cm signal, and the varia… ▽ More We modelled the two major layer of Earth's ionosphere, the F-layer and the D-layer, by a simplified spatial model with temporal variance to study the chromatic ionospheric effects on global 21-cm observations. From the analyses, we found that the magnitude of the ionospheric disruptions due to ionospheric refraction and absorption can be greater than the expected global 21-cm signal, and the variation of its magnitude can differ, depending on the ionospheric conditions. Within the parameter space adopted in the model, the shape of the global 21-cm signal is distorted after propagating through the ionosphere, while its amplitude is weakened. It is observed that the ionospheric effects do not cancel out over time, and thus should be accounted for in the foreground calibration at each timestep to account for the chromaticity introduced by the ionosphere. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: 11 pages, 13 figures, 4 tables

Spectral index of the Galactic foreground emission in the 50-87 MHz range

Authors: Marta Spinelli, Gianni Bernardi, Hugh Garsden, Lincoln J. Greehill, Anastasia Fialkov, Jayce Dowell, Daniel C. Price

Abstract: Total-power radiometry with individual meter-wave antennas is a potentially effective way to study the Cosmic Dawn ($z\sim20$) through measurement of sky brightness arising from the $21$~cm transition of neutral hydrogen, provided this can be disentangled from much stronger Galactic and extra-galactic foregrounds. In the process, measured spectra of integrated sky brightness temperature can be use… ▽ More Total-power radiometry with individual meter-wave antennas is a potentially effective way to study the Cosmic Dawn ($z\sim20$) through measurement of sky brightness arising from the $21$~cm transition of neutral hydrogen, provided this can be disentangled from much stronger Galactic and extra-galactic foregrounds. In the process, measured spectra of integrated sky brightness temperature can be used to quantify the foreground emission properties. In this work, we analyze a subset of data from the Large-aperture Experiment to Detect the Dark Age (LEDA) in the range $50-87$~MHz and constrain the foreground spectral index $β$ in the northern sky visible from mid-latitudes. We focus on two zenith-directed LEDA radiometers and study how estimates of $β$ vary with local sidereal time (LST). We correct for the effect of gain pattern chromaticity and compare estimated absolute temperatures with simulations. We develop a reference dataset consisting of 14 days of optimal condition observations. Using this dataset we estimate, for one radiometer, that $β$ varies from $-2.55$ at LST~$<6$~h to a steeper $-2.58$ at LST~$\sim13$~h, consistently with sky models and previous southern sky measurements. In the LST~$=13-24$~h range, however, we find that $β$ fluctuates between $-2.55$ and $-2.61$ (data scatter $\sim0.01$). We observe a similar $β$ vs. LST trend for the second radiometer, although with slightly smaller $|β|$, in the $-2.46<β<-2.43$ range, over $24$~h of LST (data scatter $\sim0.02$). Combining all data gathered during the extended campaign between mid-2018 to mid-2019, and focusing on the LST~$=9-12.5$~h range, we infer good instrument stability and find $-2.56<β<-2.50$ with $0.09<Δβ<0.12$. △ Less

Submitted 24 April, 2021; v1 submitted 8 November, 2020; originally announced November 2020.

Comments: revised version, new figures

Mapping discrete galaxies at cosmic dawn with 21-centimeter observations

Authors: Itamar Reis, Rennan Barkana, Anastasia Fialkov

Abstract: At cosmic dawn, the 21-centimeter signal from intergalactic hydrogen was driven by Lyman-$α$ photons from some of the earliest stars, producing a spatial pattern that reflected the distribution of galaxies at that time. Due to the large foreground, it is thought that around redshift 20 it is only observationally feasible to detect 21-cm fluctuations statistically, yielding a limited, indirect prob… ▽ More At cosmic dawn, the 21-centimeter signal from intergalactic hydrogen was driven by Lyman-$α$ photons from some of the earliest stars, producing a spatial pattern that reflected the distribution of galaxies at that time. Due to the large foreground, it is thought that around redshift 20 it is only observationally feasible to detect 21-cm fluctuations statistically, yielding a limited, indirect probe of early galaxies. Here we show that 21-cm images at cosmic dawn should actually be dominated by large (tens of comoving megaparsecs), high contrast bubbles surrounding individual galaxies. We demonstrate this using a substantially upgraded semi-numerical simulation code that realistically captures the formation and 21-cm effects of the small galaxies expected during this era. Small number statistics associated with the rarity of early galaxies, combined with the multiple scattering of photons in the blue wing of the Lyman-$α$ line, create the large bubbles and also enhance the 21-cm power spectrum by a factor of 2--7 and add to it a feature that measures the typical brightness of galaxies. These various signatures of discrete early galaxies are potentially detectable with planned experiments such as the Square Kilometer Array or the Hydrogen Epoch of Reionization Array, even if the early stars formed in dark matter halos with masses as low as $10^8\, M_\odot$, ten thousand times smaller than the Milky Way halo. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Comments: First submitted for publication on Apr. 5th, 2020

High-redshift radio galaxies: a potential new source of 21-cm fluctuations

Authors: Itamar Reis, Anastasia Fialkov, Rennan Barkana

Abstract: Radio sources are expected to have formed at high redshifts, producing an excess radiation background above the cosmic microwave background (CMB) at low frequencies. Their effect on the redshifted 21-cm signal of neutral hydrogen is usually neglected, as it is assumed that the associated background is small. Recently, an excess radio background above the level of the CMB has been proposed as one o… ▽ More Radio sources are expected to have formed at high redshifts, producing an excess radiation background above the cosmic microwave background (CMB) at low frequencies. Their effect on the redshifted 21-cm signal of neutral hydrogen is usually neglected, as it is assumed that the associated background is small. Recently, an excess radio background above the level of the CMB has been proposed as one of the possible explanations for the unusually strong 21-cm signal from redshift $z\sim 17$ reported by the EDGES collaboration. As a result, the implications of a smooth and extremely strong excess radio background on both the sky-averaged (global) 21-cm signal and its fluctuations have been considered. Here we take into account the inhomogeneity of the radio background created by a population of high-redshift galaxies, and show that it adds a new type of 21-cm fluctuations to the well-known contributions of density, velocity, Ly-$α$ coupling, heating and reionization. We find that a population of high-redshift galaxies even with a moderately-enhanced radio efficiency (unrelated to the EDGES result) can have a significant effect on the 21-cm power spectrum and global signal in models with weak X-ray heating. For models that can explain the EDGES data, we conduct a large parameter survey to explore their signatures. We show that in such models the 21-cm power spectrum at $z\sim 17$ is enhanced by up to two orders of magnitude compared to the CMB-only standard case, and the shape and time evolution of the power spectrum is significantly modified by the radio fluctuations. These fluctuations are within reach of upcoming radio interferometers. We also find that these models can be significantly constrained by current and future observations of radio sources. △ Less

Submitted 10 August, 2020; originally announced August 2020.

Comments: Comments are welcome