Search

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\c{c}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\c{c}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$\sigma$ 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., Comment: 27 pages, 21 figures, and 6 tables; accepted for publication in Astronomy and Astrophysics (A&A); language edits implemented; typos corrected

Detecting the 21-cm signal from the Epoch of Reionisation (EoR) is challenging due to the strong astrophysical foregrounds, ionospheric effects, radio frequency interference and instrumental effects. Understanding and calibrating these effects are crucial for the detection. In this work, we introduce a newly developed direction-dependent (DD) calibration algorithm DDECAL and compare its performance with an existing algorithm, SAGECAL, in the context of the LOFAR-EoR 21-cm power spectrum experiment. In our data set, the North Celestial Pole (NCP) and its flanking fields were observed simultaneously. We analyse the NCP and one of its flanking fields. The NCP field is calibrated by the standard pipeline, using SAGECAL with an extensive sky model and 122 directions, and the flanking field is calibrated by DDECAL and SAGECAL with a simpler sky model and 22 directions. Additionally, two strategies are used for subtracting Cassiopeia A and Cygnus A. The results show that DDECAL performs better at subtracting sources in the primary beam region due to the application of a beam model, while SAGECAL performs better at subtracting Cassiopeia A and Cygnus A. This indicates that including a beam model during DD calibration significantly improves the performance. The benefit is obvious in the primary beam region. We also compare the 21-cm power spectra on two different fields. The results show that the flanking field produces better upper limits compared to the NCP in this particular observation. Despite the minor differences between DDECAL and SAGECAL due to the beam application, we find that the two algorithms yield comparable 21-cm power spectra on the LOFAR-EoR data after foreground removal. Hence, the current LOFAR-EoR 21-cm power spectrum limits are not likely to depend on the DD calibration method., Comment: 28 pages, 14 figures, accepted for publication in A&A; a few equations in section 3 were corrected in this version

The detection of the 21 cm signal of neutral hydrogen from the Epoch of Reionization (EoR) is challenging due to bright foreground sources, radio frequency interference (RFI), the ionosphere, and instrumental effects. Even after correcting for these effects in the calibration step and applying foreground removal techniques, the remaining residuals in the observed 21 cm power spectra are still above the thermal noise, which is referred to as the "excess variance." We study potential causes of this excess variance based on 13 nights of data obtained with the Low-Frequency Array (LOFAR). We focused on the impact of gain errors, the sky model, and ionospheric effects on the excess variance by correlating the relevant parameters such as the gain variance over time or frequency, local sidereal time (LST), diffractive scale, and phase structure-function slope with the level of excess variance. Our analysis shows that excess variance has an LST dependence, which is related to the power from the sky. And the simulated Stokes I power spectra from bright sources and the excess variance show a similar progression over LST with the minimum power appearing at LST bin 6h to 9h. This LST dependence is also present in sky images of the residual Stokes I of the observations. In very-wide sky images, we demonstrate that the extra power comes exactly from the direction of bright and distant sources Cassiopeia A and Cygnus A with the array beam patterns. These results suggest that the level of excess variance in the 21 cm signal power spectra is related to sky effects and, hence, it depends on LST. In particular, very bright and distant sources such as Cassiopeia A and Cygnus A can dominate the effect. This is in line with earlier studies and offers a path forward toward a solution since the correlation between the sky-related effects and the excess variance is non-negligible., Comment: 29 pages, 14 figures, accepted for publication in A&A

Aims: Contamination from bright diffuse Galactic thermal and non-thermal radio emission poses crucial challenges in experiments aiming to measure the 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization. If not included in calibration, this diffuse emission can severely impact the analysis and signal extraction in 21-cm experiments. We examine large-scale diffuse Galactic emission at 122~MHz, around the North Celestial Pole, using the AARTFAAC-HBA system. Methods: In this pilot project, we present the first-ever wide-field image produced with a single sub-band of the data recorded with the AARTFAAC-HBA. We demonstrate two methods: multiscale CLEAN and shapelet decomposition, to model the diffuse emission revealed in the image. We use angular power spectrum metrics to quantify different components of the emission and compare the performance of the two diffuse structure modelling approaches. Results: We observe that the point sources dominate the angular power spectrum ($\ell(\ell+1)C_{\ell}/2\pi\equiv \Delta^2(\ell)$) of the emission in the field on scales $\ell\gtrsim 60$ ($\lesssim 3$~degree). The angular power spectrum after subtraction of compact sources is flat within $20\lesssim \ell \lesssim200$ range, suggesting that the residual power is dominated by the diffuse emission on scales $\ell\lesssim200$. The residual diffuse emission has a brightness temperature variance of $\Delta^2_{\ell=180} = (145.64 \pm 13.61)~{\rm K}^2$ at 122~MHz on angular scales of 1~degree, and is consistent with a power-law following $C_{\ell}\propto \ell^{-2.0}$ in $20\lesssim \ell \lesssim200$ range. We also find that, in the current setup, the multiscale CLEAN is suitable to model the compact and diffuse structures on a wide range of angular scales, whereas the shapelet decomposition method better models the large scales, which are of the order of a few degrees and wider., Comment: 11 pages, 9 figures, accepted (A&A)

We investigate systematic effects in direction dependent gain calibration in the context of the Low-Frequency Array (LOFAR) 21-cm Epoch of Reionization (EoR) experiment. The LOFAR EoR Key Science Project aims to detect the 21-cm signal of neutral hydrogen on interferometric baselines of $50-250 \lambda$. We show that suppression of faint signals can effectively be avoided by calibrating these short baselines using only the longer baselines. However, this approach causes an excess variance on the short baselines due to small gain errors induced by overfitting during calibration. We apply a regularised expectation-maximisation algorithm with consensus optimisation (sagecal-co) to real data with simulated signals to show that overfitting can be largely mitigated by penalising spectrally non-smooth gain solutions during calibration. This reduces the excess power with about a factor 4 in the simulations. Our results agree with earlier theoretical analysis of this bias-variance trade off and support the gain-calibration approach to the LOFAR 21-cm signal data., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

We compare various foreground removal techniques that are being utilised to remove bright foregrounds in various experiments aiming to detect the redshifted 21cm signal of neutral hydrogen from the Epoch of Reionization. In this work, we test the performance of removal techniques (FastICA, GMCA, and GPR) on 10 nights of LOFAR data and investigate the possibility of recovering the latest upper limit on the 21cm signal. Interestingly, we find that GMCA and FastICA reproduce the most recent 2$\sigma$ upper limit of $\Delta^2_{21} <$ (73)$^2$ mK$^2$ at $k=0.075~ h \mathrm{cMpc}^{-1}$, which resulted from the application of GPR. We also find that FastICA and GMCA begin to deviate from the noise-limit at \textit{k}-scales larger than $\sim 0.1 ~h \mathrm{cMpc}^{-1}$. We then replicate the data via simulations to see the source of FastICA and GMCA's limitations, by testing them against various instrumental effects. We find that no single instrumental effect, such as primary beam effects or mode-mixing, can explain the poorer recovery by FastICA and GMCA at larger \textit{k}-scales. We then test scale-independence of FastICA and GMCA, and find that lower \textit{k}-scales can be modelled by a smaller number of independent components. For larger scales ($k \gtrsim 0.1~h \mathrm{cMpc}^{-1}$), more independent components are needed to fit the foregrounds. We conclude that, the current usage of GPR by the LOFAR collaboration is the appropriate removal technique. It is both robust and less prone to overfitting, with future improvements to GPR's fitting optimisation to yield deeper limits., Comment: 15 pages, 11 figures, accepted by MNRAS on 02 November 2020

The ARCADE2 and LWA1 experiments have claimed an excess over the Cosmic Microwave Background (CMB) at low radio frequencies. If the cosmological high-redshift contribution to this radio background is between 0.1% and 22% of the CMB at 1.42 GHz, it could explain the tentative EDGES Low-Band detection of the anomalously deep absorption in the 21-cm signal of neutral hydrogen. We use the upper limit on the 21-cm signal from the Epoch of Reionization ($z=9.1$) based on 141 hours of observations with LOFAR to evaluate the contribution of the high redshift Universe to the detected radio background. Marginalizing over astrophysical properties of star-forming halos, we find (at 95% C.L.) that the cosmological radio background can be at most 9.6% of the CMB at 1.42 GHz. This limit rules out strong contribution of the high-redshift Universe to the ARCADE2 and LWA1 measurements. Even though LOFAR places limit on the extra radio background, excess of $0.1-9.6$% over the CMB (at 1.42 GHz) is still allowed and could explain the EDGES Low-Band detection. We also constrain the thermal and ionization state of the gas at $z = 9.1$ and put limits on the properties of the first star-forming objects. We find that, in agreement with the limits from EDGES High-Band data, LOFAR data constrain scenarios with inefficient X-ray sources and cases where the Universe was ionized by stars in massive halos only., Comment: Published in Monthly Notices of the Royal Astronomical Society (MNRAS). Available at https://doi.org/10.1093/mnras/staa2422

We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY and a Bayesian inference framework to constrain the parameters which describe the physical state of the IGM. We find that, if the gas heating remains negligible, an IGM with ionized fraction $\gtrsim 0.13$ and a distribution of the ionized regions with a characteristic size $\gtrsim 8 ~h^{-1}$ comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) $\gtrsim 16 ~h^{-1}$ Mpc is ruled out. For an IGM with a uniform spin temperature $T_{\rm S} \gtrsim 3$ K, no constraints on the ionized component can be computed. If the large-scale fluctuations of the signal are driven by spin temperature fluctuations, an IGM with a volume fraction $\lesssim 0.34$ of heated regions with a temperature larger than CMB, average gas temperature 7-160 K and a distribution of the heated regions with characteristic size 3.5-70 $h^{-1}$ Mpc and FWHM of $\lesssim 110$ $h^{-1}$ Mpc is ruled out. These constraints are within the 95 per cent credible intervals. With more stringent future upper limits from LOFAR at multiple redshifts, the constraints will become tighter and will exclude an increasingly large region of the parameter space., Comment: 22 pages, 15 Figures, 5 tables, Accepted for publication in MNRAS

A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due to spectral structure in the gain solutions has markedly reduced but some excess power still remains with a spectral correlation distinct from thermal noise. This excess has a spectral coherence scale of $0.25 - 0.45$\,MHz and is partially correlated between nights, especially in the foreground wedge region. The correlation is stronger between nights covering similar local sidereal times. A best 2-$\sigma$ upper limit of $\Delta^2_{21} < (73)^2\,\mathrm{mK^2}$ at $k = 0.075\,\mathrm{h\,cMpc^{-1}}$ is found, an improvement by a factor $\approx 8$ in power compared to the previously reported upper limit. The remaining excess power could be due to residual foreground emission from sources or diffuse emission far away from the phase centre, polarization leakage, chromatic calibration errors, ionosphere, or low-level radio-frequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power., Comment: 27 pages, 14 figues, accepted in MNRAS (updated with reference to accompanying paper)