-
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.
-
The LoReLi database: 21 cm signal inference with 3D radiative hydrodynamics simulations
Authors:
Romain Meriot,
Benoit Semelin
Abstract:
The Square Kilometer array is expected to measure the 21cm signal from the Epoch of Reionization (EoR) in the coming decade, and its pathfinders may provide a statistical detection even earlier. The currently reported upper limits provide tentative constraints on the astrophysical parameters of the models of the EoR.
In order to interpret such data with 3D radiative hydrodynamics simulations usi…
▽ More
The Square Kilometer array is expected to measure the 21cm signal from the Epoch of Reionization (EoR) in the coming decade, and its pathfinders may provide a statistical detection even earlier. The currently reported upper limits provide tentative constraints on the astrophysical parameters of the models of the EoR.
In order to interpret such data with 3D radiative hydrodynamics simulations using Bayesian inference, we present the latest developments of the \textsc{Licorice} code. Relying on an implementation of the halo conditional mass function to account for unresolved star formation, this code now allows accurate simulations of the EoR at $256^3$ resolution. We use this version of \textsc{Licorice} to produce the first iteration of \textsc{LoReLi}, a public dataset now containing hundreds of 21cm signals computed from radiative hydrodynamics simulations. We train a neural network on \textsc{LoReLi} to provide a fast emulator of the \textsc{Licorice} power spectra, \textsc{LorEMU}, which has $\sim 5\%$ rms error relative to the simulated signals. \textsc{LorEMU} is used in a Markov Chain Monte Carlo framework to perform Bayesian inference, first on a mock observation composed of a simulated signal and thermal noise corresponding to 100h observations with the SKA. We then apply our inference pipeline to the latest measurements from the HERA interferometer. We report constraints on the X-ray emissivity, and confirm that cold reionization scenarios are unlikely to accurately represent our Universe.
△ Less
Submitted 1 December, 2023; v1 submitted 4 October, 2023;
originally announced October 2023.
-
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.
-
The Cosmic Mach Number as an Environment Measure for the Underlying Dark Matter Density Field
Authors:
Romain Meriot,
Sadegh Khochfar,
Jose Onorbe,
Britton Smith
Abstract:
Using cosmological dark matter only simulations of a $(1.6$ Gpc$/h)^3$ volume from the Legacy simulation project, we calculate Cosmic Mach Numbers (CMN) and perform a theoretical investigation of their relation with halo properties and features of the density field to gauge their use as an measure of the environment.
CMNs calculated on individual spheres show correlations with both the overdensi…
▽ More
Using cosmological dark matter only simulations of a $(1.6$ Gpc$/h)^3$ volume from the Legacy simulation project, we calculate Cosmic Mach Numbers (CMN) and perform a theoretical investigation of their relation with halo properties and features of the density field to gauge their use as an measure of the environment.
CMNs calculated on individual spheres show correlations with both the overdensity in a region and the density gradient in the direction of the bulk flow around that region. To reduce the scatter around the median of these correlations, we introduce a new measure, the rank ordered Cosmic Mach number ($\hat{\mathcal{M}}_g$), which shows a tight correlations with the overdensity $δ=\frac{ρ-\barρ}{\barρ}$. Measures of the large scale density gradient as well as other average properties of the halo population in a region show tight correlations with $\hat{\mathcal{M}}_g$ as well. Our results in this first empirical study suggest that $\hat{\mathcal{M}}_g$ is an excellent proxy for the underlying density field and hence environment that can circumvent reliance on number density counts in a region. For scales between $10$ and $100 Mpc$/h, Mach numbers calculated using dark matter halos $(> 10^{12}$ M$_{\odot})$ that would typically host massive galaxies are consistent with theoretical predictions of the linear matter power spectrum at a level of $10\%$ due to non-linear effects of gravity. At redshifts $z\geq 3$, these deviations disappear. We also quantify errors due to missing large scale modes in simulations. Simulations of box size $\leq 1 $ Gpc/$h$ typically predict CMNs 10-30\% too small on scales of$\sim 100$ Mpc$/h$.
△ Less
Submitted 13 February, 2022; v1 submitted 8 February, 2022;
originally announced February 2022.