-
The 21-cm signal during the end stages of reionization
Authors:
Sambit K. Giri,
Michele Bianco,
Timothée Schaeffer,
Ilian T. Iliev,
Garrelt Mellema,
Aurel Schneider
Abstract:
During the epoch of reionization (EoR), the 21-cm signal allows direct observation of the neutral hydrogen (\hi{}) in the intergalactic medium (IGM). In the post-reionization era, this signal instead probes \hi{} in galaxies, which traces the dark matter density distribution. With new numerical simulations, we investigated the end stages of reionization to elucidate the transition of our Universe…
▽ More
During the epoch of reionization (EoR), the 21-cm signal allows direct observation of the neutral hydrogen (\hi{}) in the intergalactic medium (IGM). In the post-reionization era, this signal instead probes \hi{} in galaxies, which traces the dark matter density distribution. With new numerical simulations, we investigated the end stages of reionization to elucidate the transition of our Universe into the post-reionization era. Our models are consistent with the latest high-redshift measurements, including ultraviolet (UV) luminosity functions. Notably, these models consistently reproduced the evolution of the UV photon background, which is constrained from Lyman-$α$ absorption spectra. We studied the dependence of this background on the nature of photon sinks in the IGM, requiring mean free path of UV photons to be $\sim$10 comoving-megaparsecs (cMpc) during the EoR that increases gradually with time during late stages ($z\lesssim 6$). Our models revealed that the reionization of the IGM transitioned from an \textit{inside-out} to an \textit{outside-in} process when the Universe is less than 0.01 per cent neutral. During this epoch, the 21-cm signal also shifted from probing predominantly the \hi{} in the IGM to that in galaxies. Furthermore, we identified a statistically significant number of large neutral islands (with sizes up to 40 cMpc) persisting until very late stages ($5 \lesssim z \lesssim 6$) that can imprint features in Lyman-$α$ absorption spectra and also produce a knee-like feature in the 21-cm power spectrum.
△ Less
Submitted 7 March, 2024;
originally announced March 2024.
-
pyC$^2$Ray: A flexible and GPU-accelerated Radiative Transfer Framework for Simulating the Cosmic Epoch of Reionization
Authors:
Patrick Hirling,
Michele Bianco,
Sambit K. Giri,
Ilian T. Iliev,
Garrelt Mellema,
Jean-Paul Kneib
Abstract:
Detailed modelling of the evolution of neutral hydrogen in the intergalactic medium during the Epoch of Reionization, $5 \leq z \leq 20$, is critical in interpreting the cosmological signals from current and upcoming 21-cm experiments such as Low-Frequency Array (LOFAR) and the Square Kilometre Array (SKA). Numerical radiative transfer codes offer the most physically motivated approach for simulat…
▽ More
Detailed modelling of the evolution of neutral hydrogen in the intergalactic medium during the Epoch of Reionization, $5 \leq z \leq 20$, is critical in interpreting the cosmological signals from current and upcoming 21-cm experiments such as Low-Frequency Array (LOFAR) and the Square Kilometre Array (SKA). Numerical radiative transfer codes offer the most physically motivated approach for simulating the reionization process. However, they are computationally expensive as they must encompass enormous cosmological volumes while accurately capturing astrophysical processes occurring at small scales ($\lesssim\rm Mpc$). Here, we present pyC$^2$Ray, an updated version of the massively parallel ray-tracing and chemistry code, C$^2$Ray, which has been extensively employed in reionization simulations. The most time-consuming part of the code is calculating the hydrogen column density along the path of the ionizing photons. Here, we present the Accelerated Short-characteristics Octhaedral RAytracing (ASORA) method, a ray-tracing algorithm specifically designed to run on graphical processing units (GPUs). We include a modern Python interface, allowing easy and customized use of the code without compromising computational efficiency. We test pyC$^2$Ray on a series of standard ray-tracing tests and a complete cosmological simulation with volume size $(349\,\rm Mpc)^3$, mesh size of $250^3$ and approximately $10^6$ sources. Compared to the original code, pyC$^2$Ray achieves the same results with negligible fractional differences, $\sim 10^{-5}$, and a speedup factor of two orders of magnitude. Benchmark analysis shows that ASORA takes a few nanoseconds per source per voxel and scales linearly for an increasing number of sources and voxels within the ray-tracing radii.
△ Less
Submitted 2 November, 2023;
originally announced November 2023.
-
The stability of deep learning for 21cm foreground removal across various sky models and frequency-dependent systematics
Authors:
T. Chen,
M. Bianco,
E. Tolley,
M. Spinelli,
D. Forero-Sanchez,
J. P. Kneib
Abstract:
Deep learning (DL) has recently been proposed as a novel approach for 21cm foreground removal. Before applying DL to real observations, it is essential to assess its consistency with established methods, its performance across various simulation models and its robustness against instrumental systematics. This study develops a commonly used U-Net and evaluates its performance for post-reionisation…
▽ More
Deep learning (DL) has recently been proposed as a novel approach for 21cm foreground removal. Before applying DL to real observations, it is essential to assess its consistency with established methods, its performance across various simulation models and its robustness against instrumental systematics. This study develops a commonly used U-Net and evaluates its performance for post-reionisation foreground removal across three distinct sky simulation models based on pure Gaussian realisations, the Lagrangian perturbation theory, and the Planck sky model. Stable outcomes across the models are achieved provided that training and testing data align with the same model. On average, the residual foreground in the U-Net reconstructed data is $\sim$10% of the signal across angular scales at the considered redshift range. Comparable results are found with traditional approaches. However, blindly using a network trained on one model for data from another model yields inaccurate reconstructions, emphasising the need for consistent training data. The study then introduces frequency-dependent Gaussian beams and gain drifts to the test data. The network struggles to denoise data affected by "unexpected" systematics without prior information. However, after re-training consistently with systematics-contaminated data, the network effectively restores its reconstruction accuracy. This highlights the importance of incorporating prior systematics knowledge during training for successful denoising. Our work provides critical guidelines for using DL for 21cm foreground removal, tailored to specific data attributes. Notably, it is the first time that DL has been applied to the Planck sky model being most realistic foregrounds at present.
△ Less
Submitted 1 November, 2023;
originally announced November 2023.
-
BIPP: An efficient HPC implementation of the Bluebild algorithm for radio astronomy
Authors:
Emma Tolley,
Simon Frasch,
Etienne Orliac,
Shreyam Krishna,
Michele Bianco,
Sepand Kashani,
Paul Hurley,
Matthieu Simeoni,
Jean-Paul Kneib
Abstract:
The Bluebild algorithm is a new technique for image synthesis in radio astronomy which forms a least-squares estimate of the sky intensity function using the theory of sampling and interpolation operators. We present an HPC implementation of the Bluebild algorithm for radio-interferometric imaging: Bluebild Imaging++ (BIPP). BIPP is a spherical imager that leverages functional PCA to decompose the…
▽ More
The Bluebild algorithm is a new technique for image synthesis in radio astronomy which forms a least-squares estimate of the sky intensity function using the theory of sampling and interpolation operators. We present an HPC implementation of the Bluebild algorithm for radio-interferometric imaging: Bluebild Imaging++ (BIPP). BIPP is a spherical imager that leverages functional PCA to decompose the sky into distinct energy levels. The library features interfaces to C++, C and Python and is designed with seamless GPU acceleration in mind. We evaluate the accuracy and performance of BIPP on simulated observations of the upcoming Square Kilometer Array Observatory and real data from the Low-Frequency Array (LOFAR) telescope. We find that BIPP offers accurate wide-field imaging with no need for a w-term approximation and has comparable execution time with respect to the interferometric imaging libraries CASA and WSClean. Futhermore, due to the energy level decomposition, images produced with BIPP can reveal information about faint and diffuse structures before any cleaning iterations. The source code of BIPP is publicly released.
△ Less
Submitted 13 October, 2023;
originally announced October 2023.
-
Deep learning approach for identification of HII regions during reionization in 21-cm observations -- II. foreground contamination
Authors:
Michele Bianco,
Sambit. K. Giri,
David Prelogović,
Tianyue Chen,
Florent G. Mertens,
Emma Tolley,
Andrei Mesinger,
Jean-Paul Kneib
Abstract:
The upcoming Square Kilometre Array Observatory (SKAO) will produce images of neutral hydrogen distribution during the epoch of reionization by observing the corresponding 21-cm signal. However, the 21-cm signal will be subject to instrumental limitations such as noise and galactic foreground contamination which pose a challenge for accurate detection. In this study, we present the SegU-Net v2 fra…
▽ More
The upcoming Square Kilometre Array Observatory (SKAO) will produce images of neutral hydrogen distribution during the epoch of reionization by observing the corresponding 21-cm signal. However, the 21-cm signal will be subject to instrumental limitations such as noise and galactic foreground contamination which pose a challenge for accurate detection. In this study, we present the SegU-Net v2 framework, an enhanced version of our convolutional neural network, built to identify neutral and ionized regions in the 21-cm signal contaminated with foreground emission. We trained our neural network on 21-cm image data processed by a foreground removal method based on Principal Component Analysis achieving an average classification accuracy of 71 per cent between redshift $z=7$ to $11$. We tested SegU-Net v2 against various foreground removal methods, including Gaussian Process Regression, Polynomial Fitting, and Foreground-Wedge Removal. Results show comparable performance, highlighting SegU-Net v2's independence on these pre-processing methods. Statistical analysis shows that a perfect classification score with $AUC=95\%$ is possible for $8<z<10$. While the network prediction lacks the ability to correctly identify ionized regions at higher redshift and differentiate well the few remaining neutral regions at lower redshift due to low contrast between 21-cm signal, noise and foreground residual in images. Moreover, as the photon sources driving reionization are expected to be located inside ionised regions, we show that SegU-Net v2 can be used to correctly identify and measure the volume of isolated bubbles with $V_{\rm ion}>(10\, {\rm cMpc})^3$ at $z>9$, for follow-up studies with infrared/optical telescopes to detect these sources.
△ Less
Submitted 28 February, 2024; v1 submitted 5 April, 2023;
originally announced April 2023.
-
PINION: Physics-informed neural network for accelerating radiative transfer simulations for cosmic reionization
Authors:
Damien Korber,
Michele Bianco,
Emma Tolley,
Jean-Paul Kneib
Abstract:
With the advent of the Square Kilometre Array Observatory (SKAO), scientists will be able to directly observe the Epoch of Reionization by mapping the distribution of neutral hydrogen at different redshifts. While physically motivated results can be simulated with radiative transfer codes, these simulations are computationally expensive and can not readily produce the required scale and resolution…
▽ More
With the advent of the Square Kilometre Array Observatory (SKAO), scientists will be able to directly observe the Epoch of Reionization by mapping the distribution of neutral hydrogen at different redshifts. While physically motivated results can be simulated with radiative transfer codes, these simulations are computationally expensive and can not readily produce the required scale and resolution simultaneously. Here we introduce the Physics-Informed neural Network for reIONization (PINION), which can accurately and swiftly predict the complete 4-D hydrogen fraction evolution from the smoothed gas and mass density fields from pre-computed N-body simulation. We trained PINION on the C$^2$-Ray simulation outputs and a physics constraint on the reionization chemistry equation is enforced. With only five redshift snapshots and a propagation mask as a simplistic approximation of the ionizing photon mean free path, PINION can accurately predict the entire reionization history between $z=6$ and $12$. We evaluate the accuracy of our predictions by analysing the dimensionless power spectra and morphology statistics estimations against C$^2$-Ray results. We show that while the network's predictions are in good agreement with simulation to redshift $z>7$, the network's accuracy suffers for $z<7$ primarily due to the oversimplified propagation mask. We motivate how PINION performance can be drastically improved and potentially generalized to large-scale simulations.
△ Less
Submitted 9 March, 2023; v1 submitted 29 August, 2022;
originally announced August 2022.
-
Antenna characterization for the HIRAX experiment
Authors:
Emily R. Kuhn,
Benjamin R. B. Saliwanchik,
Kevin Bandura,
Michele Bianco,
H. Cynthia Chiang,
Devin Crichton,
Meiling Deng,
Sindhu Gaddam,
Kit Gerodias,
Austin Gumba,
Maile Harris,
Kavilan Moodley,
V. Mugundhan,
Laura Newburgh,
Jeffrey Peterson,
Elizabeth Pieters,
Anna R. Polish,
Alexandre Refregier,
Ajith Sampath,
Mario G. Santos,
Onkabetse Sengate,
Jonathan Sievers,
Ema Smith,
Will Tyndall,
Anthony Walters
, et al. (2 additional authors not shown)
Abstract:
The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) aims to improve constraints on the dark energy equation of state through measurements of large-scale structure at high redshift ($0.8<z<2.5$), while serving as a state-of-the-art fast radio burst detector. Bright galactic foregrounds contaminate the 400--800~MHz HIRAX frequency band, so meeting the science goals will require precise…
▽ More
The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) aims to improve constraints on the dark energy equation of state through measurements of large-scale structure at high redshift ($0.8<z<2.5$), while serving as a state-of-the-art fast radio burst detector. Bright galactic foregrounds contaminate the 400--800~MHz HIRAX frequency band, so meeting the science goals will require precise instrument characterization. In this paper we describe characterization of the HIRAX antenna, focusing on measurements of the antenna beam and antenna noise temperature.
Beam measurements of the current HIRAX antenna design were performed in an anechoic chamber and compared to simulations. We report measurement techniques and results, which find a broad and symmetric antenna beam for $ν<$650MHz, and elevated cross-polarization levels and beam asymmetries for $ν>$700MHz. Noise temperature measurements of the HIRAX feeds were performed in a custom apparatus built at Yale. In this system, identical loads, one cryogenic and the other at room temperature, are used to take a differential (Y-factor) measurement from which the noise of the system is inferred. Several measurement sets have been conducted using the system, involving CHIME feeds as well as four of the HIRAX active feeds. These measurements give the first noise temperature measurements of the HIRAX feed, revealing a $\sim$60K noise temperature (relative to 30K target) with 40K peak- to-peak frequency-dependent features, and provide the first demonstration of feed repeatability. Both findings inform current and future feed designs.
△ Less
Submitted 25 July, 2022;
originally announced July 2022.
-
Crucial Factors for Lyα Transmission in the Reionizing Intergalactic Medium: Infall Motion, HII Bubble Size, and Self-shielded Systems
Authors:
Hyunbae Park,
Intae Jung,
Hyunmi Song,
Pierre Ocvirk,
Paul R. Shapiro,
Taha Dawoodbhoy,
Ilian T. Iliev,
Kyungjin Ahn,
Michele Bianco,
Hyo Jeong Kim
Abstract:
Using the CoDa II simulation, we study the Ly$α$ transmissivity of the intergalactic medium (IGM) during reionization. At $z>6$, a typical galaxy without an active galactic nucleus fails to form a proximity zone around itself due to the overdensity of the surrounding IGM. The gravitational infall motion in the IGM makes the resonance absorption extend to the red side of Ly$α$, suppressing the tran…
▽ More
Using the CoDa II simulation, we study the Ly$α$ transmissivity of the intergalactic medium (IGM) during reionization. At $z>6$, a typical galaxy without an active galactic nucleus fails to form a proximity zone around itself due to the overdensity of the surrounding IGM. The gravitational infall motion in the IGM makes the resonance absorption extend to the red side of Ly$α$, suppressing the transmission up to roughly the circular velocity of the galaxy. In some sight lines, an optically thin blob generated by a supernova in a neighboring galaxy results in a peak feature, which can be mistaken for a blue peak. Redward of the resonance absorption, the damping-wing opacity correlates with the global IGM neutral fraction and the UV magnitude of the source galaxy. Brighter galaxies tend to suffer lower opacity because they tend to reside in larger HII regions, and the surrounding IGM transmits redder photons, which are less susceptible to attenuation, owing to stronger infall velocity. The HII regions are highly nonspherical, causing both sight-line-to-sight-line and galaxy-to-galaxy variation in opacity. Also, self-shielded systems within HII regions strongly attenuate the emission for certain sight lines. All these factors add to the transmissivity variation, requiring a large sample size to constrain the average transmission. The variation is largest for fainter galaxies at higher redshift. The 68\% range of the transmissivity is similar to or greater than the median for galaxies with $M_{\rm UV}\ge-21$ at $z\ge7$, implying that more than a hundred galaxies would be needed to measure the transmission to 10\% accuracy.
△ Less
Submitted 6 December, 2021; v1 submitted 22 May, 2021;
originally announced May 2021.
-
Deep learning approach for identification of HII regions during reionization in 21-cm observations
Authors:
Michele Bianco,
Sambit. K. Giri,
Ilian T. Iliev,
Garrelt Mellema
Abstract:
The upcoming Square Kilometre Array (SKA-Low) will map the distribution of neutral hydrogen during reionization, and produce a tremendous amount of 3D tomographic data. These images cubes will be subject to instrumental limitations, such as noise and limited resolution. Here we present SegU-Net, a stable and reliable method for identification of neutral and ionized regions in these images. SegU-Ne…
▽ More
The upcoming Square Kilometre Array (SKA-Low) will map the distribution of neutral hydrogen during reionization, and produce a tremendous amount of 3D tomographic data. These images cubes will be subject to instrumental limitations, such as noise and limited resolution. Here we present SegU-Net, a stable and reliable method for identification of neutral and ionized regions in these images. SegU-Net is a U-Net architecture based convolutional neural network (CNN) for image segmentation. It is capable of segmenting our image data into meaningful features (ionized and neutral regions) with greater accuracy compared to previous methods. We can estimate the true ionization history from our mock observation of SKA with an observation time of 1000 h with more than 87 per cent accuracy. We also show that SegU-Net can be used to recover various topological summary statistics, such as size distributions and Betti numbers, with a relative difference of only a few per cent. These summary statistics characterise the non-Gaussian nature of the reionization process.
△ Less
Submitted 24 June, 2021; v1 submitted 12 February, 2021;
originally announced February 2021.
-
The impact of inhomogeneous subgrid clumping on cosmic reionization II: modelling stochasticity
Authors:
Michele Bianco,
Ilian T. Iliev,
Kyungjin Ahn,
Sambit K. Giri,
Yi Mao,
Hyunbae Park,
Paul R. Shapiro
Abstract:
Small-scale density fluctuations can significantly affect reionization but are typically modelled quite crudely. Unresolved fluctuations in numerical simulations and analytical calculations are included using a gas clumping factor, typically assumed to be independent of the local environment. In Paper I, we presented an improved, local density-dependent model for the sub-grid gas clumping. Here we…
▽ More
Small-scale density fluctuations can significantly affect reionization but are typically modelled quite crudely. Unresolved fluctuations in numerical simulations and analytical calculations are included using a gas clumping factor, typically assumed to be independent of the local environment. In Paper I, we presented an improved, local density-dependent model for the sub-grid gas clumping. Here we extend this using an empirical stochastic model based on the results from high-resolution numerical simulations which fully resolve all relevant fluctuations. Our model reproduces well both the mean density-clumping relation and its scatter. We applied our stochastic model, along with the mean clumping one and the Paper I deterministic model, to create a large-volume realisation of the clumping field, and used these in radiative transfer simulations of cosmic reionization. Our results show that the simplistic mean clumping model delays reionization compared to local density-dependent models, despite producing fewer recombinations overall. This is due to the very different spatial distribution of clumping, resulting in much higher photoionization rates in the latter cases. The mean clumping model produces smaller HII regions throughout most of the reionization, but those percolate faster at late times. It also causes a significant delay in the 21-cm fluctuations peak and yields lower non-Gaussianity and many fewer bright pixels in the PDF distribution. The stochastic density-dependent model shows relatively minor differences from the deterministic one, mostly concentrated around overlap, where it significantly suppresses the 21-cm fluctuations, and at the bright tail of the 21-cm PDFs, where it produces noticeably more bright pixels.
△ Less
Submitted 5 May, 2021; v1 submitted 5 January, 2021;
originally announced January 2021.
-
The impact of inhomogeneous subgrid clumping on cosmic reionization
Authors:
Yi Mao,
Jun Koda,
Paul R. Shapiro,
Ilian T. Iliev,
Garrelt Mellema,
Hyunbae Park,
Kyungjin Ahn,
Michele Bianco
Abstract:
Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corr…
▽ More
Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation.
△ Less
Submitted 12 November, 2019; v1 submitted 6 June, 2019;
originally announced June 2019.
-
Supergravity, $α$-attractors and primordial non-Gaussianity
Authors:
Nicola Bartolo,
Domenico Matteo Bianco,
Raul Jimenez,
Sabino Matarrese,
Licia Verde
Abstract:
We compute in detail how deviations from Einstein gravity at the inflation energy scale could appear as non-Gaussian features in the sky. To illustrate this we use multi-field $α-$attractor models in the framework of supergravity to realise inflation. We find no obvious obstacle for having choices of model's parameters that generate non-Gaussian features of the equilateral and local type at the…
▽ More
We compute in detail how deviations from Einstein gravity at the inflation energy scale could appear as non-Gaussian features in the sky. To illustrate this we use multi-field $α-$attractor models in the framework of supergravity to realise inflation. We find no obvious obstacle for having choices of model's parameters that generate non-Gaussian features of the equilateral and local type at the $\cal O$(1) level in the $f_{\rm NL}$ non-Gaussianity parameter, thus being potentially detectable in future cosmological surveys. This non-Gaussianity has its origin in either the non-canonical kinetic term (which, in turn, is an immediate consequence of assuming an hyperbolic geometry of the moduli space), the interactions of the fields in the potential or the $α-$parameter, or a combination of these three. This opens up the exciting possibility of constraining the law of gravity at energy scales close to the Planck one.
△ Less
Submitted 5 October, 2018; v1 submitted 11 May, 2018;
originally announced May 2018.
