LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

Future high-redshift 21-cm experiments will suffer from a high degree of contamination, due both to astrophysical foregrounds and to non-astrophysical and instrumental effects. In order to reliably extract the cosmological signal from the observed data, it is essential to understand very well all data components and their influence on the extracted signal. Here we present simulated astrophysical foregrounds data cubes and discuss their possible statistical effects on the data. The foreground maps are produced assuming 5 • × 5 • windows that match those expected to be observed by the LOFAR epoch of reionization (EoR) key science project. We show that with the expected LOFAR-EoR sky and receiver noise levels, which amount to ≈52 mK at 150 MHz after 400 h of total observing time, a simple polynomial fit allows a statistical reconstruction of the signal. We also show that the polynomial fitting will work for maps with realistic yet idealized instrument response, i.e. a response that includes only a uniform uv coverage as a function of frequency and ignores many other uncertainties. Polarized Galactic synchrotron maps that include internal polarization and a number of Faraday screens along the line of sight are also simulated. The importance of these stems from the fact that the LOFAR instrument, in common with all current interferometric EoR experiments, has an instrumentally polarized response. 5 We assume a Lambda cold dark matter ( CDM) universe with b = 0.04, m = 0.26, = 0.738 and H 0 = 70.8 k ms −1 Mpc −1 .

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 56 13 mK D < ( ) 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 Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window both on the time of the formation of the first stars and accreting black holes and the subsequent period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an Executive SummaryThe Square Kilometre Array (SKA) will have a low frequency component (AA-low/SKA-low 1 ) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe (see SKA Memo 125). It is during this phase that the first building blocks of the galaxies that we see around us today, including our own Milky Way, were formed. It is a crucial period for understanding the history of the Universe and one for which we have currently very little observational data.We divide the period into two different phases based on the physical processes which affect the Intergalactic Medium. The first period, which we call the Cosmic Dawn, saw the formation of the first stars and accreting black holes, which changed the quantum state of the still neutral Intergalactic Medium. The second period, known as the Epoch of Reionization, is the one during which large areas between the galaxies were photo-ionized by the radiation produced in galaxies and which ended when the Intergalactic Medium had become completely ionized.Observations of the redshifted 21-cm line with SKA will provide a new and unique window on the entire period of Cosmic Dawn and Reionization. The signal is sensitive to the emergence of the first stellar populations, radiation from growing massive black holes and the formation of larger groups of galaxies and bright quasars. At the same time it maps the distribution of most of the baryonic matter in the Universe. The study of the redshifted 21cm line will teach us fundamental new things about the earliest phases of structure formation and cosmology. It even has the potential to lead to the discovery of new physical phenomena. Here we present an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.The redshifted 21cm signal will be analyzed with different techniques, which each come with their own requirements for the SKA: (i) Tomography, (ii) power-spectra and higher-order statistics, (iii) hydrogen absorption, (iv) global/total-intensity signal. Whereas all precursors/pathfinders aim to study the signal statistically through its power spectrum, SKA will be able to image the neutral hydrogen distribution directly and its focus will therefore be more on tomograph...

A new upper limit on the 21-cm signal power spectrum at a redshift of z ≈ 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-σ upper limit of ∆ 2 21 < (73) 2 mK 2 at k = 0.075 h cMpc −1 is found, an improvement by a factor ≈ 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.

The formalism of Wiener ltering is developed here for the purpose of reconstructing the large scale structure of the universe from noisy, sparse and incomplete data. The method is based on a linear minimum variance solution, given data and an assumed prior model which speci es the covariance matrix of the eld to be reconstructed. While earlier applications of the Wiener ler have focused on estimation, namely suppressing the noise in the measured quantities, we extend the method here to perform both prediction and dynamical reconstruction. The Wiener lter is used to predict the values of unmeasured quantities, such as the density eld in un-sampled regions of space, or to deconvolve blurred data. The method is developed, within the context of linear gravitational instability theory, to perform dynamical reconstruction of one eld which is dynamically related to some other observed eld. This is the case, for example, in the reconstruction of the real space galaxy distribution from its redshift distribution or the prediction of the radial velocity eld from the observed density eld.When the eld to be reconstructed is a Gaussian random eld, such as the primordial perturbation eld predicted by the canonical model of cosmology, the Wiener lter can 1 be pushed to its fullest potential. In such a case the Wiener estimator coincides with the Bayesian estimator designed to maximize the posterior probability. The Wiener lter can be also derived by assuming a quadratic regularization function, in analogy with thè Maximum Entropy' method. The mean eld obtained by the minimal variance solution can be supplemented with constrained realizations of the Gaussian eld to create random realizations of the residual from the mean. Subject Headings: cosmology: large scale structure of universe; analytical{ methods: data analysis

Context. The cosmological 21 cm line promises to be a formidable tool for cosmology, allowing the investigation of the end of the so-called dark ages, when the first galaxies formed. Aims. Astrophysical foregrounds are expected to be about three orders of magnitude greater than the cosmological signal and therefore represent a serious contamination of the cosmological 21 cm line. Detailed knowledge of both their intensity and polarization structure on the relevant angular scale of 1-30 arcmin will be essential for extracting the cosmological signal from the data. Methods. We present the first results from a series of observations conducted with the Westerbork telescope in the 140-160 MHz range with a 2 arcmin resolution aimed at characterizing the properties of the foregrounds for epoch of reionization experiments. The polarization data were analysed through the rotation measure synthesis technique. We computed total intensity and polarization angular power spectra. Results. For the first time we have detected fluctuations in the Galactic diffuse emission on scales greater than 13 arcmin at 150 MHz, in the low Galactic latitude area known as Fan region, centred at α = 3 h 10 m , δ = 65 • 30 . Those fluctuations have an rms of 14 K. The total intensity power spectrum shows a power-law behaviour down to ∼900 with slope β I = −2.2 ± 0.3. The detection of diffuse emission at smaller angular scales is limited by residual point sources. We measured an rms confusion noise of ∼3 mJy beam −1 . Diffuse polarized emission was also detected for the first time at this frequency. The polarized signal shows complex structure both spatially and along the line of sight. The polarization power spectrum is not affected by residual point sources and is only limited by the thermal noise. It shows a power-law behaviour down to ∼2700 with slope β P = −1.65±0.15. The rms of polarization fluctuations is 7.2 K on 4 arcmin scales. Conclusions. The measured total intensity fluctuations are used to estimate the foreground contamination on the cosmological signal.By extrapolating the spectrum of total intensity emission, we find a contamination of δT = ( + 1)C I /2π ∼ 5.7 K on 5 arcmin scales and a corresponding rms value of ∼18.3 K at the same angular scale. The level of the polarization power spectrum is δT ∼ 3.3 K on 5 arcmin scales. However, the Fan region cannot be taken as representative of other sky regions, given its exceptionally bright polarized signal, but is likely to represent an upper limit on the sky brightness at moderate and high Galactic latitude.

We describe and compare several post-correlation radio frequency interference classification methods. As data sizes of observations grow with new and improved telescopes, the need for completely automated, robust methods for radio frequency interference mitigation is pressing. We investigated several classification methods and find that, for the data sets we used, the most accurate among them is the SumThreshold method. This is a new method formed from a combination of existing techniques, including a new way of thresholding. This iterative method estimates the astronomical signal by carrying out a surface fit in the time-frequency plane. With a theoretical accuracy of 95% recognition and an approximately 0.1% false probability rate in simple simulated cases, the method is in practice as good as the human eye in finding RFI. In addition it is fast, robust, does not need a data model before it can be executed and works in almost all configurations with its default parameters. The method has been compared using simulated data with several other mitigation techniques, including one based upon the singular value decomposition of the time-frequency matrix, and has shown better results than the rest.Comment: 14 pages, 12 figures (11 in colour). The software that was used in the article can be downloaded from http://www.astro.rug.nl/rfi-software