E. Chapman

Foreground Removal using FastICA: A Showcase of LOFAR-EoR

We introduce a new implementation of the fastica algorithm on simulated Low Frequency Array Epoch of Reionization data with the aim of accurately removing the foregrounds and extracting the 21-cm reionization signal. We find that the method successfully removes the foregrounds with an average fitting error of 0.5 per cent and that the 2D and 3D power spectra are recovered across the frequency range. We find that for scales above several point spread function scales, the 21-cm variance is successfully recovered though there is evidence of noise leakage into the reconstructed foreground components. We find that this blind independent component analysis technique provides encouraging results without the danger of prior foreground assumptions.

The scale of the problem: recovering images of reionization with Generalized Morphological Component Analysis

The accurate and precise removal of 21-cm foregrounds from Epoch of Reionization redshifted 21-cm emission data is essential if we are to gain insight into an unexplored cosmological era. We apply a non-parametric technique, Generalized Morphological Component Analysis or GMCA, to simulated LOFAR-EoR data and show that it has the ability to clean the foregrounds with high accuracy. We recover the 21-cm 1D, 2D and 3D power spectra with high accuracy across an impressive range of frequencies and scales. We show that GMCA preserves the 21-cm phase information, especially when the smallest spatial scale data is discarded. While it has been shown that LOFAR-EoR image recovery is theoretically possible using image smoothing, we add that wavelet decomposition is an efficient way of recovering 21-cm signal maps to the same or greater order of accuracy with more flexibility. By comparing the GMCA output residual maps (equal to the noise, 21-cm signal and any foreground fitting errors) with the 21-cm maps at one frequency and discarding the smaller wavelet scale information, we find a correlation coefficient of 0.689, compared to 0.588 for the equivalently smoothed image. Considering only the central 50% of the maps, these coefficients improve to 0.905 and 0.605 respectively and we conclude that wavelet decomposition is a significantly more powerful method to denoise reconstructed 21-cm maps than smoothing.

Upper Limits on the 21 cm Epoch of Reionization Power Spectrum from One Night with LOFAR

We present the first limits on the Epoch of Reionization 21 cm H I power spectra, in the redshift range z = 7.9–10.6, using the Low-Frequency Array (LOFAR) High-Band Antenna (HBA). In total, 13.0 hr of data were used from observations centered on the North Celestial Pole. After subtraction of the sky model and the noise bias, we detect a non-zero Δ^2_I = (56 ± 13 mK)^2 (1-σ) excess variance and a best 2-σ upper limit of Δ^2_(21) < (79.6 mK)^2 at k = 0.053 h cMpc^(−1) in the range z = 9.6–10.6. The excess variance decreases when optimizing the smoothness of the direction- and frequency-dependent gain calibration, and with increasing the completeness of the sky model. It is likely caused by (i) residual side-lobe noise on calibration baselines, (ii) leverage due to nonlinear effects, (iii) noise and ionosphere-induced gain errors, or a combination thereof. Further analyses of the excess variance will be discussed in forthcoming publications.

The Scale of the Problem : Recovering Images of Reionization with GMCA

The accurate and precise removal of 21-cm foregrounds from Epoch of Reionization redshifted 21-cm emission data is essential if we are to gain insight into an unexplored cosmological era. We apply a non-parametric technique, Generalized Morphological Component Analysis or GMCA, to simulated LOFAR-EoR data and show that it has the ability to clean the foregrounds with high accuracy. We recover the 21-cm 1D, 2D and 3D power spectra with high accuracy across an impressive range of frequencies and scales. We show that GMCA preserves the 21-cm phase information, especially when the smallest spatial scale data is discarded. While it has been shown that LOFAR-EoR image recovery is theoretically possible using image smoothing, we add that wavelet decomposition is an efficient way of recovering 21-cm signal maps to the same or greater order of accuracy with more flexibility. By comparing the GMCA output residual maps (equal to the noise, 21-cm signal and any foreground fitting errors) with the 21-cm maps at one frequency and discarding the smaller wavelet scale information, we find a correlation coefficient of 0.689, compared to 0.588 for the equivalently smoothed image. Considering only the central 50% of the maps, these coefficients improve to 0.905 and 0.605 respectively and we conclude that wavelet decomposition is a significantly more powerful method to denoise reconstructed 21-cm maps than smoothing.

Probing reionization with LOFAR using 21-cm redshift space distortions

One of the most promising ways to study the epoch of reionization (EoR) is through radio observations of the redshifted 21-cm line emission from neutral hydrogen. These observations are complicated ...

Constraining the epoch of reionization with the variance statistic: simulations of the LOFAR case

The Epoch of Reionization (EoR) is the epoch in which most of the neutral gas in the Universe was re-ionized by the radiation from the first stars and galaxies. Many projects are underway to detect the redshifted 21 cm signal of the neutral hydrogen from the EoR with the low frequency radio telescopes. These experiments aim for a statistical detection of the signal due to its very low signal-to-noise ratio. We study extraction of the sample variance of the signal to constrain the global properties of the EoR. We show that the LOw Frequency ARray (LOFAR) should be able to detect the signal with a significance of 3 standard deviations in 600 hours of integration. Additionally, it should be able to constrain the timing and duration of reionization with uncertainties of 0.4 and 1 redshifts, respectively, with 95 percent confidence. We also show the upper limits on the signal detection obtained from the analysis of 115 hours of the observed data with LOFAR. In the observed redshift range of 7 to 11, our current best upper limits are 25 mK at 1 MHz, 12 arcmin spectral and spatial resolutions, respectively.

Linear polarization structures in LOFAR observations of the interstellar medium in the 3C 196 field

Aims. This study aims to characterize linear polarization structures in LOFAR observations of the interstellar medium (ISM) in the 3C 196 field, one of the primary fields of the LOFAR-Epoch of Reionization key science project. Methods. We have used the high band antennas (HBA) of LOFAR to image this region and rotation measure (RM) synthesis to unravel the distribution of polarized structures in Faraday depth. Results. The brightness temperature of the detected Galactic emission is 5−15 K in polarized intensity and covers the range from –3 to +8 rad m-2 in Faraday depth. The most interesting morphological feature is a strikingly straight filament at a Faraday depth of +0.5 rad m-2 running from north to south, right through the centre of the field and parallel to the Galactic plane. There is also an interesting system of linear depolarization canals conspicuous in an image showing the peaks of Faraday spectra. We used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. For the first time, we see some common morphology in the RM cubes made at 150 and 350 MHz. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies and previous LOFAR observations. Based on our results, we determined physical parameters of the ISM and proposed a simple model that may explain the observed distribution of the intervening magneto-ionic medium. Conclusions. The mean line-of-sight magnetic field component, B∥, is determined to be 0.3 ± 0.1 μG and its spatial variation across the 3C 196 field is 0.1 μG. The filamentary structure is probably an ionized filament in the ISM, located somewhere within the Local Bubble. This filamentary structure shows an excess in thermal electron density (neB∥> 6.2 cm-3μG) compared to its surroundings.

Probing ionospheric structures using the LOFAR radio telescope

LOFAR is the LOw-Frequency Radio interferometer ARray located at midlatitude (52°53′N). Here we present results on ionospheric structures derived from 29 LOFAR nighttime observations during the winters of 2012/2013 and 2013/2014. We show that LOFAR is able to determine differential ionospheric total electron content values with an accuracy better than 0.001 total electron content unit = 1016m−2 over distances ranging between 1 and 100 km. For all observations the power law behavior of the phase structure function is confirmed over a long range of baseline lengths, between 1 and 80 km, with a slope that is, in general, larger than the 5/3 expected for pure Kolmogorov turbulence. The measured average slope is 1.89 with a one standard deviation spread of 0.1. The diffractive scale, i.e., the length scale where the phase variance is 1rad2, is shown to be an easily obtained single number that represents the ionospheric quality of a radio interferometric observation. A small diffractive scale is equivalent to high phase variability over the field of view as well as a short time coherence of the signal, which limits calibration and imaging quality. For the studied observations the diffractive scales at 150 MHz vary between 3.5 and 30 km. A diffractive scale above 5 km, pertinent to about 90% of the observations, is considered sufficient for the high dynamic range imaging needed for the LOFAR epoch of reionization project. For most nights the ionospheric irregularities were anisotropic, with the structures being aligned with the Earth magnetic field in about 60% of the observations.

Lunar occultation of the diffuse radio sky: LOFAR measurements between 35 and 80 MHz

We present radio observations of the Moon between 35 and 80 MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies 35 z > 12) and the Epoch of Reionization (12 > z > 5).

Effects of the sources of reionization on 21-cm redshift-space distortions

The observed 21 cm signal from the epoch of reionization will be distorted along the line of sight by the peculiar velocities of matter particles. These redshift-space distortions will affect the contrast in the signal and will also make it anisotropic. This anisotropy contains information about the cross-correlation between the matter density field and the neutral hydrogen field, and could thus potentially be used to extract information about the sources of reionization. In this paper, we study a collection of simulated reionization scenarios assuming different models for the sources of reionization. We show that the 21 cm anisotropy is best measured by the quadrupole moment of the power spectrum. We find that, unless the properties of the reionization sources are extreme in some way, the quadrupole moment evolves very predictably as a function of global neutral fraction. This predictability implies that redshift-space distortions are not a very sensitive tool for distinguishing between reionization sources. However, the quadrupole moment can be used as a model-independent probe for constraining the reionization history. We show that such measurements can be done to some extent by first-generation instruments such as LOFAR, while the SKA should be able to measure the reionization history using the quadrupole moment of the power spectrum to great accuracy.