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.

Despite much recent theoretical and observational progress in our knowledge of the early universe, many fundamental questions remain only partially answered. Here, we review the latest achievements and persisting problems in the understanding of first cosmic structure formation.

Radiative transfer (RT) simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come within reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non‐trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on five test problems: (0) basic physics; (1) isothermal H ii region expansion; (2) H ii region expansion with evolving temperature; (3) I‐front trapping and shadowing by a dense clump and (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analysed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multifrequency treatment approach, resulting in different thicknesses of the ionized‐neutral transition regions and the temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To further this aim all test inputs and outputs are made publicly available in digital form.

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.