Showing 1–40 of 40 results for author: Brentjens, M A

Unintended electromagnetic radiation from Starlink satellites detected with LOFAR between 110 and 188 MHz

Authors: F. Di Vruno, B. Winkel, C. G. Bassa, G. I. G. Józsa, M. A. Brentjens, A. Jessner, S. Garrington

Abstract: We report on observations of 68 satellites belonging to the SpaceX Starlink constellation with the LOFAR radio telescope. Radiation associated with Starlink satellites was detected at observing frequencies between 110 and 188 MHz. A combination of broad-band features, covering the entire observed bandwidth, as well as narrow-band (bandwidth < 12.2 kHz) emission at frequencies of 125, 135, 143.05,… ▽ More We report on observations of 68 satellites belonging to the SpaceX Starlink constellation with the LOFAR radio telescope. Radiation associated with Starlink satellites was detected at observing frequencies between 110 and 188 MHz. A combination of broad-band features, covering the entire observed bandwidth, as well as narrow-band (bandwidth < 12.2 kHz) emission at frequencies of 125, 135, 143.05, 150, and 175 MHz, was observed. The presence and properties of both the narrow- and broad-band features vary between satellites at different orbital altitudes. While the narrow-band detections at 143.05 MHz can be attributed to reflections of radar signals from the French GRAVES Space Surveillance Radar, the signal properties of the broad- and narrow-band features at the other frequencies suggest that this radiation is intrinsic to the Starlink satellites and it is seen for 47 out of the 68 Starlink satellites that were observed. We observed spectral power flux densities vary from 0.1 to 10 Jy for broad-band radiation, to 10 to 500 Jy for some of the narrow-band radiation, equivalent to electric field strengths of up to 49 dB[uV/m] (as measured at a 10 m distance from the satellites, with a measurement bandwidth of 120 kHz). In addition, we present equivalent power flux density simulations of the full Starlink phase 1 constellation, as well as other satellite constellations, for one frequency band allocated to radio astronomy by the International Telecommunication Union (ITU). With these, we calculate the maximum radiation level that each satellite constellation would need to have to comply with regulatory limits for intended emissions in that band. However, these limits do not apply if the radiation is unintended, that is to say if it does not originate from intentionally radiated signals for radio communication or other purposes. (shortened) △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 20 pages, 14 figures, in press

The Scintillating Tail of Comet C/2020 F3 (Neowise)

Authors: R. A. Fallows, B. Forte, M. Mevius, M. A. Brentjens, C. G. Bassa, M. M. Bisi, A. Offringa, G. Shaifullah, C. Tiburzi, H. Vedantham, P. Zucca

Abstract: Context. The occultation of a radio source by the plasma tail of a comet can be used to probe structure and dynamics in the tail. Such occultations are rare, and the occurrence of scintillation, due to small-scale density variations in the tail, remains somewhat controversial. Aims. A detailed observation taken with the Low-Frequency Array (LOFAR) of a serendipitous occultation of the compact radi… ▽ More Context. The occultation of a radio source by the plasma tail of a comet can be used to probe structure and dynamics in the tail. Such occultations are rare, and the occurrence of scintillation, due to small-scale density variations in the tail, remains somewhat controversial. Aims. A detailed observation taken with the Low-Frequency Array (LOFAR) of a serendipitous occultation of the compact radio source 3C196 by the plasma tail of comet C/2020 F3 (Neowise) is presented. 3C196 tracked almost perpendicularly behind the tail, providing a unique profile cut only a short distance downstream from the cometary nucleus itself. Methods. Interplanetary scintillation (IPS) is observed as the rapid variation of the intensity received of a compact radio source due to density variations in the solar wind. IPS in the signal received from 3C196 was observed for five hours, covering the full transit behind the plasma tail of comet C/2020 F3 (Neowise) on 16 July 2020, and allowing an assessment of the solar wind in which the comet and its tail are embedded. Results. The results reveal a sudden and strong enhancement in scintillation which is unequivocally attributable to the plasma tail. The strongest scintillation is associated with the tail boundaries, weaker scintillation is seen within the tail, and previously-unreported periodic variations in scintillation are noted, possibly associated with individual filaments of plasma. Furthermore, contributions from the solar wind and comet tail are separated to measure a sharp decrease in the velocity of material within the tail, suggesting a steep velocity shear resulting in strong turbulence along the tail boundary △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: Accepted for publication in Astronomy and Astrophysics, 8 pages, 9 figures

Journal ref: A&A 667, A57 (2022)

A numerical study of 21-cm signal suppression and noise increase in direction-dependent calibration of LOFAR data

Authors: M. Mevius, F. Mertens, L. V. E. Koopmans, A. R. Offringa, S. Yatawatta, M. A. Brentjens, E. Chapman, B. Ciardi, H. Gan, B. K. Gehlot, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, G. Mellema, V. N. Pandey, S. Zaroubi

Abstract: 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 λ$. We show that suppression of faint signals can effectively be avoided by calibrating these short base… ▽ More 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 λ$. 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. △ Less

Submitted 3 November, 2021; originally announced November 2021.

Comments: This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

Journal ref: MNRAS 2021

LOFAR imaging of the solar corona during the 2015 March 20 solar eclipse

Authors: A. M. Ryan, P. T. Gallagher, E. P. Carley, M. A. Brentjens, P. C. Murphy, C. Vocks, D. E. Morosan, H. Reid, J. Magdalenic, F. Breitling, P. Zucca, R. Fallows, G. Mann, A. Kerdraon, R. Halfwerk

Abstract: The solar corona is a highly-structured plasma which can reach temperatures of more than ~2 MK. At low frequencies (decimetric and metric wavelengths), scattering and refraction of electromagnetic waves are thought to considerably increase the imaged radio source sizes (up to a few arcminutes). However, exactly how source size relates to scattering due to turbulence is still subject to investigati… ▽ More The solar corona is a highly-structured plasma which can reach temperatures of more than ~2 MK. At low frequencies (decimetric and metric wavelengths), scattering and refraction of electromagnetic waves are thought to considerably increase the imaged radio source sizes (up to a few arcminutes). However, exactly how source size relates to scattering due to turbulence is still subject to investigation. The theoretical predictions relating source broadening to propagation effects have not been fully confirmed by observations due to the rarity of high spatial resolution observations of the solar corona at low frequencies. Here, the LOw Frequency ARray (LOFAR) was used to observe the solar corona at 120-180 MHz using baselines of up to ~3.5 km (corresponding to a resolution of ~1-2') during the partial solar eclipse of 2015 March 20. A lunar de-occultation technique was used to achieve higher spatial resolution (~0.6') than that attainable via standard interferometric imaging (~2.4'). This provides a means of studying the contribution of scattering to apparent source size broadening. It was found that the de-occultation technique reveals a more structured quiet corona that is not resolved from standard imaging, implying scattering may be overestimated in this region when using standard imaging techniques. However, an active region source was measured to be ~4' using both de-occultation and standard imaging. This may be explained by the increased scattering of radio waves by turbulent density fluctuations in active regions, which is more severe than in the quiet Sun. △ Less

Submitted 10 February, 2021; originally announced February 2021.

Comments: 8 pages, 6 figures, 1 table

Journal ref: A&A 648, A43 (2021)

The AARTFAAC Cosmic Explorer: observations of the 21-cm power spectrum in the EDGES absorption trough

Authors: B. K. Gehlot, F. G. Mertens, L. V. E. Koopmans, A. R. Offringa, A. Shulevski, M. Mevius, M. A. Brentjens, M. Kuiack, V. N. Pandey, A. Rowlinson, A. M. Sardarabadi, H. K. Vedantham, R. A. M. J. Wijers, S. Yatawatta, S. Zaroubi

Abstract: The 21-cm absorption feature reported by the EDGES collaboration is several times stronger than that predicted by traditional astrophysical models. If genuine, a deeper absorption may lead to stronger fluctuations on the 21-cm signal on degree scales (up to 1~Kelvin in rms), allowing these fluctuations to be detectable in nearly 50~times shorter integration times compared to previous predictions.… ▽ More The 21-cm absorption feature reported by the EDGES collaboration is several times stronger than that predicted by traditional astrophysical models. If genuine, a deeper absorption may lead to stronger fluctuations on the 21-cm signal on degree scales (up to 1~Kelvin in rms), allowing these fluctuations to be detectable in nearly 50~times shorter integration times compared to previous predictions. We commenced the "AARTFAAC Cosmic Explorer" (ACE) program, that employs the AARTFAAC wide-field imager, to measure or set limits on the power spectrum of the 21-cm fluctuations in the redshift range $z = 17.9-18.6$ ($Δν= 72.36-75.09$~MHz) corresponding to the deep part of the EDGES absorption feature. Here, we present first results from two LST bins: 23.5-23.75h and 23.5-23.75h, each with 2~h of data, recorded in `semi drift-scan' mode. We demonstrate the application of the new ACE data-processing pipeline (adapted from the LOFAR-EoR pipeline) on the AARTFAAC data. We observe that noise estimates from the channel and time-differenced Stokes~$V$ visibilities agree with each other. After 2~h of integration and subtraction of bright foregrounds, we obtain $2σ$ upper limits on the 21-cm power spectrum of $Δ_{21}^2 < (8139~\textrm{mK})^2$ and $Δ_{21}^2 < (8549~\textrm{mK})^2$ at $k = 0.144~h\,\textrm{cMpc}^{-1}$ for the two LST bins. Incoherently averaging the noise bias-corrected power spectra for the two LST bins yields an upper limit of $Δ_{21}^2 < (7388~\textrm{mK})^2$ at $k = 0.144~h\,\textrm{cMpc}^{-1}$. These are the deepest upper limits thus far at these redshifts. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Tied-array holography with LOFAR

Authors: P. Salas, M. A. Brentjens, D. D. Bordenave, J. B. R. Oonk, H. J. A. Röttgering

Abstract: A radio interferometer uses time delays to maximize its response to radiation coming from a particular direction. These time delays compensate for differences in the time of arrival of the wavefront at the different elements of the interferometer, and for delays in the instrument's signal chain. If the radio interferometer is operated as a phased array (tied array), the time delays cannot be accou… ▽ More A radio interferometer uses time delays to maximize its response to radiation coming from a particular direction. These time delays compensate for differences in the time of arrival of the wavefront at the different elements of the interferometer, and for delays in the instrument's signal chain. If the radio interferometer is operated as a phased array (tied array), the time delays cannot be accounted for after an observation, so they must be determined in advance. Our aim is to characterize the time delays between the stations in the core of the LOw Frequency ARray (LOFAR). We used radio holography to determine the time delays for the core stations of LOFAR (innermost 3.5 km). Using the multibeaming capability of LOFAR we map the voltage beam faster than with a raster scan, while simultaneously calibrating the observed beam continuously. For short radio holographic observations (60 s and 600 s) of 3C196, 3C147, and 3C48 we are able to derive time delays with errors of less than one nanosecond. After applying the derived time delays to the beamformer, the beam shows residuals of less than $20\%$ with respect to the theoretical beam shape. Tied-array holography could be a way towards semi-real-time beam calibration for the Square Kilometer Array. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Comments: Accepted for publication in A&A. 8 pages, 10 figures

Journal ref: A&A 635, A207 (2020)

Improved upper limits on the 21-cm signal power spectrum of neutral hydrogen at $\boldsymbol{z \approx 9.1}$ from LOFAR

Authors: F. G. Mertens, M. Mevius, L. V. E Koopmans, A. R. Offringa, G. Mellema, S. Zaroubi, M. A. Brentjens, H. Gan, B. K. Gehlot, V. N. Pandey, A. M. Sardarabadi, H. K. Vedantham, S. Yatawatta, K. M. B. Asad, B. Ciardi, E. Chapman, S. Gazagnes, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, V. Jelić, R. Kooistra, R. Mondal, J. Schaye , et al. (1 additional authors not shown)

Abstract: 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… ▽ More 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-$σ$ upper limit of $Δ^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. △ Less

Submitted 19 February, 2020; v1 submitted 17 February, 2020; originally announced February 2020.

Comments: 27 pages, 14 figues, accepted in MNRAS (updated with reference to accompanying paper)

The first power spectrum limit on the 21-cm signal of neutral hydrogen during the Cosmic Dawn at z=20-25 from LOFAR

Authors: B. K. Gehlot, F. G. Mertens, L. V. E. Koopmans, M. A. Brentjens, S. Zaroubi, B. Ciardi, A. Ghosh, M. Hatef, I. T. Iliev, V. Jelić, R. Kooistra, F. Krause, G. Mellema, M. Mevius, M. Mitra, A. R. Offringa, V. N. Pandey, A. M. Sardarabadi, J. Schaye, M. B. Silva, H. K. Vedantham, S. Yatawatta

Abstract: Observations of the redshifted 21-cm hyperfine line of neutral hydrogen from early phases of the Universe such as Cosmic Dawn and the Epoch of Reionization promise to open a new window onto the early formation of stars and galaxies. We present the first upper limits on the power spectrum of redshifted 21-cm brightness temperature fluctuations in the redshift range $z = 19.8 - 25.2$ ($54-68$ MHz fr… ▽ More Observations of the redshifted 21-cm hyperfine line of neutral hydrogen from early phases of the Universe such as Cosmic Dawn and the Epoch of Reionization promise to open a new window onto the early formation of stars and galaxies. We present the first upper limits on the power spectrum of redshifted 21-cm brightness temperature fluctuations in the redshift range $z = 19.8 - 25.2$ ($54-68$ MHz frequency range) using 14 hours of data obtained with the LOFAR-Low Band Antenna (LBA) array. We also demonstrate the application of a multiple pointing calibration technique to calibrate the LOFAR-LBA dual-pointing observations centred on the North Celestial Pole and the radio galaxy 3C220.3. We observe an unexplained excess of $\sim 30-50\%$ in Stokes $I$ noise compared to Stokes $V$ for the two observed fields, which decorrelates on $\gtrsim 12$ seconds and might have a physical origin. We show that enforcing smoothness of gain errors along frequency direction during calibration reduces the additional variance in Stokes $I$ compared Stokes $V$ introduced by the calibration on sub-band level. After subtraction of smooth foregrounds, we achieve a $2σ$ upper limit on the 21-cm power spectrum of $Δ_{21}^2 < (14561\,\text{mK})^2$ at $k\sim 0.038\,h\,\text{cMpc}^{-1}$ and $Δ_{21}^2 < (14886\,\text{mK})^2$ at $k\sim 0.038 \,h\,\text{cMpc}^{-1}$ for the 3C220 and NCP fields respectively and both upper limits are consistent with each other. The upper limits for the two fields are still dominated by systematics on most $k$ modes. △ Less

Submitted 20 July, 2019; v1 submitted 18 September, 2018; originally announced September 2018.

Comments: 19 pages, 16 figures, accepted for publication in MNRAS

Cobalt: A GPU-based correlator and beamformer for LOFAR

Authors: P. Chris Broekema, J. Jan David Mol, Ronald Nijboer, Alexander S. van Amesfoort, Michiel A. Brentjens, G. Marcel Loose, Wouter F. A. Klijn, John W. Romein

Abstract: For low-frequency radio astronomy, software correlation and beamforming on general purpose hardware is a viable alternative to custom designed hardware. LOFAR, a new-generation radio telescope centered in the Netherlands with international stations in Germany, France, Ireland, Latvia, Poland, Sweden and the UK, has successfully used software real-time processors based on IBM Blue Gene technology s… ▽ More For low-frequency radio astronomy, software correlation and beamforming on general purpose hardware is a viable alternative to custom designed hardware. LOFAR, a new-generation radio telescope centered in the Netherlands with international stations in Germany, France, Ireland, Latvia, Poland, Sweden and the UK, has successfully used software real-time processors based on IBM Blue Gene technology since 2004. Since then, developments in technology have allowed us to build a system based on commercial off-the-shelf components that combines the same capabilities with lower operational cost. In this paper we describe the design and implementation of a GPU-based correlator and beamformer with the same capabilities as the Blue Gene based systems. We focus on the design approach taken, and show the challenges faced in selecting an appropriate system. The design, implementation and verification of the software system shows the value of a modern test-driven development approach. Operational experience, based on three years of operations, demonstrates that a general purpose system is a good alternative to the previous supercomputer-based system or custom-designed hardware. △ Less

Submitted 30 April, 2018; v1 submitted 15 January, 2018; originally announced January 2018.

Journal ref: Astronomy and Computing, Volume 23, April 2018, Pages 180-192

Wide-field LOFAR-LBA power-spectra analyses: Impact of calibration, polarization leakage and ionosphere

Authors: B. K. Gehlot, L. V. E. Koopmans, A. G. de Bruyn, S. Zaroubi, M. A. Brentjens, K. M. B. Asad, M. Hatef, V. Jelic, M. Mevius, A. R. Offringa, V. N. Pandey, S. Yatawatta

Abstract: Contamination due to foregrounds (Galactic and Extra-galactic), calibration errors and ionospheric effects pose major challenges in detection of the cosmic 21 cm signal in various Epoch of Reionization (EoR) experiments. We present the results of a pilot study of a field centered on 3C196 using LOFAR Low Band (56-70 MHz) observations, where we quantify various wide field and calibration effects su… ▽ More Contamination due to foregrounds (Galactic and Extra-galactic), calibration errors and ionospheric effects pose major challenges in detection of the cosmic 21 cm signal in various Epoch of Reionization (EoR) experiments. We present the results of a pilot study of a field centered on 3C196 using LOFAR Low Band (56-70 MHz) observations, where we quantify various wide field and calibration effects such as gain errors, polarized foregrounds, and ionospheric effects. We observe a `pitchfork' structure in the 2D power spectrum of the polarized intensity in delay-baseline space, which leaks into the modes beyond the instrumental horizon (EoR/CD window). We show that this structure largely arises due to strong instrumental polarization leakage ($\sim30\%$) towards {Cas\,A} ($\sim21$ kJy at 81 MHz, brightest source in northern sky), which is far away from primary field of view. We measure an extremely small ionospheric diffractive scale ($r_{\text{diff}} \approx 430$ m at 60 MHz) towards {Cas\,A} resembling pure Kolmogorov turbulence compared to $r_{\text{diff}} \sim3 - 20$ km towards zenith at 150 MHz for typical ionospheric conditions. This is one of the smallest diffractive scales ever measured at these frequencies. Our work provides insights in understanding the nature of aforementioned effects and mitigating them in future Cosmic Dawn observations (e.g. with SKA-low and HERA) in the same frequency window. △ Less

Submitted 26 April, 2018; v1 submitted 22 September, 2017; originally announced September 2017.

Comments: 20 pages, 11 figures, accepted for publication in MNRAS

Upper limits on the 21-cm Epoch of Reionization power spectrum from one night with LOFAR

Authors: A. H. Patil, S. Yatawatta, L. V. E. Koopmans, A. G. de Bruyn, M. A. Brentjens, S. Zaroubi, K. M. B. Asad, M. Hatef, V. Jelic, M. Mevius, A. R. Offringa, V. N. Pandey, H. Vedantham, F. B. Abdalla, W. N. Brouw, E. Chapman, B. Ciardi, B. K. Gehlot, A. Ghosh, G. Harker, I. T. Iliev, K. Kakiichi, S. Majumdar, M. B. Silva, G. Mellema , et al. (3 additional authors not shown)

Abstract: We present the first limits on the Epoch of Reionization (EoR) 21-cm HI power spectra, in the redshift range $z=7.9-10.6$, using the Low-Frequency Array (LOFAR) High-Band Antenna (HBA). In total 13\,h of data were used from observations centred on the North Celestial Pole (NCP). After subtraction of the sky model and the noise bias, we detect a non-zero $Δ^2_{\rm I} = (56 \pm 13 {\rm mK})^2$ (1-… ▽ More We present the first limits on the Epoch of Reionization (EoR) 21-cm HI power spectra, in the redshift range $z=7.9-10.6$, using the Low-Frequency Array (LOFAR) High-Band Antenna (HBA). In total 13\,h of data were used from observations centred on the North Celestial Pole (NCP). After subtraction of the sky model and the noise bias, we detect a non-zero $Δ^2_{\rm I} = (56 \pm 13 {\rm mK})^2$ (1-$σ$) excess variance and a best 2-$σ$ upper limit of $Δ^2_{\rm 21} < (79.6 {\rm 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 non-linear effects, (iii) noise and ionosphere-induced gain errors, or a combination thereof. Further analyses of the excess variance will be discussed in forthcoming publications. △ Less

Submitted 28 February, 2017; originally announced February 2017.

Comments: 17 pages, 8 figures, accepted by ApJ

Probing Ionospheric Structures using the LOFAR radio telescope

Authors: M. Mevius, S. van der Tol, V. N. Pandey, H. K. Vedantham, M. A. Brentjens, A. G. de Bruyn, F. B. Abdalla, K. M. B. Asad, J. D. Bregman, W. N. Brouw, S. Bus, E. Chapman, B. Ciardi, E. R. Fernandez, A. Ghosh, G. Harker, I. T. Iliev, V. Jelić, S. Kazemi, L. V. E. Koopmans, J. E. Noordam, A. R. Offringa, A. H. Patil, R. J. van Weeren, S. Wijnholds , et al. (2 additional authors not shown)

Abstract: LOFAR is the LOw Frequency Radio interferometer ARray located at mid-latitude ($52^{\circ} 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 TEC values with an accuracy better than 1 mTECU over distances ranging between 1 and 100 km.… ▽ More LOFAR is the LOw Frequency Radio interferometer ARray located at mid-latitude ($52^{\circ} 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 TEC values with an accuracy better than 1 mTECU 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 $1\, \mathrm{rad^2}$, 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. △ Less

Submitted 15 June, 2016; originally announced June 2016.

Comments: accepted for publication in Radio Science

Journal ref: Radio Sci. 51 (2016) 927-941

Polarization leakage in epoch of reionization windows: II. Primary beam model and direction dependent calibration

Authors: K. M. B. Asad, L. V. E. Koopmans, V. Jelić, A. Ghosh, F. B. Abdalla, M. A. Brentjens, A. G. de Bruyn, B. Ciardi, B. K. Gehlot, I. T. Iliev, M. Mevius, V. N. Pandey, S. Yatawatta, S. Zaroubi

Abstract: Leakage of diffuse polarized emission into Stokes I caused by the polarized primary beam of the instrument might mimic the spectral structure of the 21-cm signal coming from the epoch of reionization (EoR) making their separation difficult. Therefore, understanding polarimetric performance of the antenna is crucial for a successful detection of the EoR signal. Here, we have calculated the accuracy… ▽ More Leakage of diffuse polarized emission into Stokes I caused by the polarized primary beam of the instrument might mimic the spectral structure of the 21-cm signal coming from the epoch of reionization (EoR) making their separation difficult. Therefore, understanding polarimetric performance of the antenna is crucial for a successful detection of the EoR signal. Here, we have calculated the accuracy of the nominal model beam of LOFAR in predicting the leakage from Stokes I to Q, U by comparing them with the corresponding leakage of compact sources actually observed in the 3C295 field. We have found that the model beam has errors of less than or equal to 10% on the predicted levels of leakage of ~1% within the field of view, i. e. if the leakage is taken out perfectly using this model the leakage will reduce to $10^{-3}$ of the Stokes I flux. If similar levels of accuracy can be obtained in removing leakage from Stokes Q, U to I, we can say, based on the results of our previous paper, that the removal of this leakage using this beam model would ensure that the leakage is well below the expected EoR signal in almost the whole instrumental k-space of the cylindrical power spectrum. We have also shown here that direction dependent calibration can remove instrumentally polarized compact sources, given an unpolarized sky model, very close to the local noise level. △ Less

Submitted 27 July, 2016; v1 submitted 15 April, 2016; originally announced April 2016.

Comments: 13 pages, 15 figures, accepted in MNRAS on 26 July 2016

Galactic interstellar filaments as probed by LOFAR and Planck

Authors: S. Zaroubi, V. Jelić, A. G. de Bruyn, F. Boulanger, A. Bracco, R. Kooistra, M. I. R. Alves, M. A. Brentjens, K. Ferrière, T. Ghosh, L. V. E. Koopmans, F. Levrier, M. -A. Miville-Deschênes, L. Montier, V. N. Pandey, J. D. Soler

Abstract: Recent Low Frequency Array (LOFAR) observations at 115-175 MHz of a field at medium Galactic latitudes (centered at the bright quasar 3C196) have shown striking filamentary structures in polarization that extend over more than 4 degrees across the sky. In addition, the Planck satellite has released full sky maps of the dust emission in polarization at 353GHz. The LOFAR data resolve Faraday structu… ▽ More Recent Low Frequency Array (LOFAR) observations at 115-175 MHz of a field at medium Galactic latitudes (centered at the bright quasar 3C196) have shown striking filamentary structures in polarization that extend over more than 4 degrees across the sky. In addition, the Planck satellite has released full sky maps of the dust emission in polarization at 353GHz. The LOFAR data resolve Faraday structures along the line of sight, whereas the Planck dust polarization maps probe the orientation of the sky projected magnetic field component. Hence, no apparent correlation between the two is expected. Here we report a surprising, yet clear, correlation between the filamentary structures, detected with LOFAR, and the magnetic field orientation, probed by the Planck satellite. This finding points to a common, yet unclear, physical origin of the two measurements in this specific area in the sky. A number of follow-up multi- frequency studies are proposed to shed light on this unexpected finding. △ Less

Submitted 26 August, 2015; originally announced August 2015.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS Letters

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

Authors: V. Jelić, A. G. de Bruyn, V. N. Pandey, M. Mevius, M. Haverkorn, M. A. Brentjens, L. V. E. Koopmans, S. Zaroubi, F. B. Abdalla, K. M. B. Asad, S. Bus, E. Chapman, B. Ciardi, E. R. Fernandez, A. Ghosh, G. Harker, I. T. Iliev, H. Jensen, S. Kazemi, G. Mellema, A. R. Offringa, A. H. Patil, H. K. Vedantham, S. Yatawatta

Abstract: This study aims to characterize linear polarization structures in LOFAR observations of the interstellar medium (ISM) in the 3C196 field, one of the primary fields of the LOFAR-Epoch of Reionization key science project. 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. Th… ▽ More This study aims to characterize linear polarization structures in LOFAR observations of the interstellar medium (ISM) in the 3C196 field, one of the primary fields of the LOFAR-Epoch of Reionization key science project. 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. The brightness temperature of the detected Galactic emission is $5-15~{\rm K}$ in polarized intensity and covers the range from -3 to +8 ${\rm rad~m^{-2}}$ in Faraday depth. The most interesting morphological feature is a strikingly straight filament at a Faraday depth of $+0.5~{\rm 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~{\rm 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. The mean line-of-sight magnetic field component, $B_\parallel$, is determined to be $0.3\pm0.1~{\rm μG}$ and its spatial variation across the 3C196 field is $0.1~{\rm μ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 ($n_e B_\parallel>6.2~{\rm cm^{-3}μG}$) compared to its surroundings. △ Less

Submitted 21 September, 2015; v1 submitted 26 August, 2015; originally announced August 2015.

Comments: 16 pages, 10 figures, aceppted for publication in A&A

Journal ref: A&A 583, A137 (2015)

Polarization leakage in Epoch of Reionization windows: I. LOFAR observations of the 3C196 field

Authors: K. M. B. Asad, L. V. E. Koopmans, V. Jelić, V. N. Pandey, A. Ghosh, F. B. Abdalla, G. Bernardi, M. A. Brentjens, A. G. de Bruyn, S. Bus, B. Ciardi, E. Chapman, S. Daiboo, E. R. Fernandez, G. Harker, I. T. Iliev, H. Jensen, O. Martinez-Rubi, G. Mellema, M. Mevius, A. R. Offringa, A. H. Patil, J. Schaye, R. M. Thomas, S. van der Tol , et al. (3 additional authors not shown)

Abstract: Detection of the 21-cm signal coming from the epoch of reionization (EoR) is challenging especially because, even after removing the foregrounds, the residual Stokes $I$ maps contain leakage from polarized emission that can mimic the signal. Here, we discuss the instrumental polarization of LOFAR and present realistic simulations of the leakages between Stokes parameters. From the LOFAR observatio… ▽ More Detection of the 21-cm signal coming from the epoch of reionization (EoR) is challenging especially because, even after removing the foregrounds, the residual Stokes $I$ maps contain leakage from polarized emission that can mimic the signal. Here, we discuss the instrumental polarization of LOFAR and present realistic simulations of the leakages between Stokes parameters. From the LOFAR observations of polarized emission in the 3C196 field, we have quantified the level of polarization leakage caused by the nominal model beam of LOFAR, and compared it with the EoR signal using power spectrum analysis. We found that at 134--166 MHz, within the central 4$^\circ$ of the field the $(Q,U)\rightarrow I$ leakage power is lower than the EoR signal at $k<0.3$ Mpc$^{-1}$. The leakage was found to be localized around a Faraday depth of 0, and the rms of the leakage as a fraction of the rms of the polarized emission was shown to vary between 0.2-0.3\%, both of which could be utilized in the removal of leakage. Moreover, we could define an `EoR window' in terms of the polarization leakage in the cylindrical power spectrum above the PSF-induced wedge and below $k_\parallel\sim 0.5$ Mpc$^{-1}$, and the window extended up to $k_\parallel\sim 1$ Mpc$^{-1}$ at all $k_\perp$ when 70\% of the leakage had been removed. These LOFAR results show that even a modest polarimetric calibration over a field of view of $\lesssim 4^\circ$ in the future arrays like SKA will ensure that the polarization leakage remains well below the expected EoR signal at the scales of 0.02-1 Mpc$^{-1}$. △ Less

Submitted 13 May, 2015; v1 submitted 5 March, 2015; originally announced March 2015.

Comments: 19 pages, 13 figures. Accepted for publication in MNRAS

Initial LOFAR observations of Epoch of Reionization windows: II. Diffuse polarized emission in the ELAIS-N1 field

Authors: V. Jelic, A. G. de Bruyn, M. Mevius, F. B. Abdalla, K. M. B. Asad, G. Bernardi, M. A. Brentjens, S. Bus, E. Chapman, B. Ciardi, S. Daiboo, E. R. Fernandez, A. Ghosh, G. Harker, H. Jensen, S. Kazemi, L. V. E. Koopmans, P. Labropoulos, O. Martinez-Rubi, G. Mellema, A. R. Offringa, V. N. Pandey, A. H. Patil, R. M. Thomas, H. K. Vedantham , et al. (84 additional authors not shown)

Abstract: This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes.… ▽ More This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%. △ Less

Submitted 8 July, 2014; originally announced July 2014.

Comments: 12 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 568, A101 (2014)

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

Authors: Ajinkya H. Patil, Saleem Zaroubi, Emma Chapman, Vibor Jelić, Geraint Harker, Filipe B. Abdalla, Khan M. B. Asad, Gianni Bernardi, Michiel A. Brentjens, A. G. de Bruyn, Sander Bus, Benedetta Ciardi, Soobash Daiboo, Elizabeth R. Fernandez, Abhik Ghosh, Hannes Jensen, Sanaz Kazemi, Léon V. E. Koopmans, Panagiotis Labropoulos, Maaijke Mevius, Oscar Martinez, Garrelt Mellema, Andre. R. Offringa, Vishhambhar N. Pandey, Joop Schaye , et al. (5 additional authors not shown)

Abstract: Several experiments are underway to detect the cosmic redshifted 21-cm signal from neutral hydrogen from the Epoch of Reionization (EoR). Due to their very low signal-to-noise ratio, these observations aim for a statistical detection of the signal by measuring its power spectrum. We investigate the extraction of the variance of the signal as a first step towards detecting and constraining the glob… ▽ More Several experiments are underway to detect the cosmic redshifted 21-cm signal from neutral hydrogen from the Epoch of Reionization (EoR). Due to their very low signal-to-noise ratio, these observations aim for a statistical detection of the signal by measuring its power spectrum. We investigate the extraction of the variance of the signal as a first step towards detecting and constraining the global history of the EoR. Signal variance is the integral of the signal's power spectrum, and it is expected to be measured with a high significance. We demonstrate this through results from a simulation and parameter estimation pipeline developed for the Low Frequency Array (LOFAR)-EoR experiment. We show that LOFAR should be able to detect the EoR in 600 hours of integration using the variance statistic. Additionally, the redshift ($z_r$) and duration ($Δz$) of reionization can be constrained assuming a parametrization. We use an EoR simulation of $z_r = 7.68$ and $Δz = 0.43$ to test the pipeline. We are able to detect the simulated signal with a significance of 4 standard deviations and extract the EoR parameters as $z_r = 7.72^{+0.37}_{-0.18}$ and $Δz = 0.53^{+0.12}_{-0.23}$ in 600 hours, assuming that systematic errors can be adequately controlled. We further show that the significance of detection and constraints on EoR parameters can be improved by measuring the cross-variance of the signal by cross-correlating consecutive redshift bins. △ Less

Submitted 16 June, 2014; v1 submitted 16 January, 2014; originally announced January 2014.

Comments: 13 pages, 14 figures, Accepted for publication in MNRAS

Journal ref: MNRAS 443 1113 (2014)

Chromatic effects in the 21 cm global signal from the cosmic dawn

Authors: H. K. Vedantham, L. V. E. Koopmans, A. G. de Bruyn, S. J. Wijnholds, B. Ciardi, M. A. Brentjens

Abstract: The redshifted 21 cm brightness distribution from neutral hydrogen is a promising probe into the cosmic dark ages, cosmic dawn, and re-ionization. LOFAR's Low Band Antennas (LBA) may be used in the frequency range 45 MHz to 85 MHz (30>z>16) to measure the sky averaged redshifted 21 cm brightness temperature as a function of frequency, or equivalently, cosmic redshift. These low frequencies are aff… ▽ More The redshifted 21 cm brightness distribution from neutral hydrogen is a promising probe into the cosmic dark ages, cosmic dawn, and re-ionization. LOFAR's Low Band Antennas (LBA) may be used in the frequency range 45 MHz to 85 MHz (30>z>16) to measure the sky averaged redshifted 21 cm brightness temperature as a function of frequency, or equivalently, cosmic redshift. These low frequencies are affected by strong Galactic foreground emission that is observed through frequency dependent ionospheric and antenna beam distortions which lead to chromatic mixing of spatial structure into spectral structure. Using simple models, we show that (i) the additional antenna temperature due to ionospheric refraction and absorption are at a \sim 1% level--- 2 to 3 orders of magnitude higher than the expected 21 cm signal, and have an approximate ν^{-2} dependence, (ii) ionospheric refraction leads to a knee-like modulation on the sky spectrum at ν\approx 4\times plasma frequency. Using more realistic simulations, we show that in the measured sky spectrum, more than 50% of the 21 cm signal variance can be lost to confusion from foregrounds and chromatic effects. We conclude that foregrounds and chromatic mixing may not be subtracted as generic functions of frequency as previously thought, but must rather be carefully modeled using additional priors and interferometric measurements. △ Less

Submitted 10 June, 2013; originally announced June 2013.

Comments: 16 pages, 12 figures

Probing reionization with LOFAR using 21-cm redshift space distortions

Authors: Hannes Jensen, Kanan K. Datta, Garrelt Mellema, Emma Chapman, Filipe B. Abdalla, Ilian T. Iliev, Yi Mao, Mario G. Santos, Paul R. Shapiro, Saleem Zaroubi, G. Bernardi, M. A. Brentjens, A. G. de Bruyn, B. Ciardi, G. J. A. Harker, V. Jelić, S. Kazemi, L. V. E. Koopmans, P. Labropoulos, O. Martinez, A. R. Offringa, V. N. Pandey, J. Schaye, R. M. Thomas, V. Veligatla , et al. (2 additional authors not shown)

Abstract: 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 by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, bu… ▽ More 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 by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, but they also encode information about cosmology and astrophysics. We study the effects of redshift space distortions on the power spectrum of 21-cm radiation from the EoR using large scale $N$-body and radiative transfer simulations. We quantify the anisotropy introduced in the 21-cm power spectrum by redshift space distortions and show how it evolves as reionization progresses and how it relates to the underlying physics. We go on to study the effects of redshift space distortions on LOFAR observations, taking instrument noise and foreground subtraction into account. We find that LOFAR should be able to directly observe the power spectrum anisotropy due to redshift space distortions at spatial scales around $k \sim 0.1$ Mpc$^{-1}$ after $\gtrsim$ 1000 hours of integration time. At larger scales, sample errors become a limiting factor, while at smaller scales detector noise and foregrounds make the extraction of the signal problematic. Finally, we show how the astrophysical information contained in the evolution of the anisotropy of the 21-cm power spectrum can be extracted from LOFAR observations, and how it can be used to distinguish between different reionization scenarios. △ Less

Submitted 19 July, 2013; v1 submitted 22 March, 2013; originally announced March 2013.

Comments: Accepted for publication in MNRAS

Initial deep LOFAR observations of Epoch of Reionization windows: I. The North Celestial Pole

Authors: S. Yatawatta, A. G. de Bruyn, M. A. Brentjens, P. Labropoulos, V. N. Pandey, S. Kazemi, S. Zaroubi, L. V. E. Koopmans, A. R. Offringa, V. Jelic, O. Martinez Rubi, V. Veligatla, S. J. Wijnholds, W. N. Brouw, G. Bernardi, B. Ciardi, S. Daiboo, G. Harker, G. Mellema, J. Schaye, R. Thomas, H. Vedantham, E. Chapman, F. B. Abdalla, A. Alexov , et al. (64 additional authors not shown)

Abstract: The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospec… ▽ More The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospective observing windows for the LOFAR EoR project will be centered at the North Celestial Pole (NCP). We present results from observations of the NCP window using the LOFAR highband antenna (HBA) array in the frequency range 115 MHz to 163 MHz. The data were obtained in April 2011 during the commissioning phase of LOFAR. We used baselines up to about 30 km. With about 3 nights, of 6 hours each, effective integration we have achieved a noise level of about 100 microJy/PSF in the NCP window. Close to the NCP, the noise level increases to about 180 microJy/PSF, mainly due to additional contamination from unsubtracted nearby sources. We estimate that in our best night, we have reached a noise level only a factor of 1.4 above the thermal limit set by the noise from our Galaxy and the receivers. Our continuum images are several times deeper than have been achieved previously using the WSRT and GMRT arrays. We derive an analytical explanation for the excess noise that we believe to be mainly due to sources at large angular separation from the NCP. △ Less

Submitted 11 January, 2013; v1 submitted 8 January, 2013; originally announced January 2013.

Comments: 19 pages, 24 figures. Draft version with low resolution images. Accepted on 08/01/2013 Astronomy & Astrophysics. Abstract abridged. Version with high resolution images: http://www.astro.rug.nl/~yatawatta/ncp_eor.pdf

Journal ref: Astronomy & Astrophysics, Volume 550, id.A136, 17 pp., 2013

The LOFAR radio environment

Authors: A. R. Offringa, A. G. de Bruyn, S. Zaroubi, G. van Diepen, O. Martinez-Ruby, P. Labropoulos, M. A. Brentjens, B. Ciardi, S. Daiboo, G. Harker, V. Jelic, S. Kazemi, L. V. E. Koopmans, G. Mellema, V. N. Pandey, R. F. Pizzo, J. Schaye, H. Vedantham, V. Veligatla, S. J. Wijnholds, S. Yatawatta, P. Zarka, A. Alexov, J. Anderson, A. Asgekar , et al. (71 additional authors not shown)

Abstract: Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency… ▽ More Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills. △ Less

Submitted 1 October, 2012; originally announced October 2012.

Comments: Accepted for publication in A&A. 15 figures, 16 pages

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

Authors: Emma Chapman, Filipe B. Abdalla, J. Bobin, J. -L. Starck, Geraint Harker, Vibor Jelic, Panagiotis Labropoulos, Saleem Zaroubi, Michiel A. Brentjens, A. G. de Bruyn, L. V. E. Koopmans

Abstract: 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… ▽ More 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. △ Less

Submitted 10 January, 2013; v1 submitted 21 September, 2012; originally announced September 2012.

Comments: 13 pages, 12 figures, accepted by MNRAS

Journal ref: MNRAS, 2013, Vol. 429, Pg. 165-176

Imaging neutral hydrogen on large-scales during the Epoch of Reionization with LOFAR

Authors: S. Zaroubi, A. G. de Bruyn, G. Harker, R. M. Thomas, P. Labropolous, V. Jelic, L. V. E. Koopmans, M. A. Brentjens, G. Bernardi, B. Ciardi, S. Daiboo, S. Kazemi, O. Martinez-Rubi, G. Mellema, A. R. Offringa, V. N. Pandey, J. Schaye, V. Veligatla, H. Vedantham, S. Yatawatta

Abstract: The first generation of redshifted 21 cm detection experiments, carried out with arrays like LOFAR, MWA and GMRT, will have a very low signal-to-noise ratio per resolution element (\sim 0.2). In addition, whereas the variance of the cosmological signal decreases on scales larger than the typical size of ionization bubbles, the variance of the formidable galactic foregrounds increases, making it ha… ▽ More The first generation of redshifted 21 cm detection experiments, carried out with arrays like LOFAR, MWA and GMRT, will have a very low signal-to-noise ratio per resolution element (\sim 0.2). In addition, whereas the variance of the cosmological signal decreases on scales larger than the typical size of ionization bubbles, the variance of the formidable galactic foregrounds increases, making it hard to disentangle the two on such large scales. The poor sensitivity on small scales on the one hand, and the foregrounds effect on large scales on the other hand, make direct imaging of the Epoch of Reionization of the Universe very difficult, and detection of the signal therefore is expected to be statistical.Despite these hurdles, in this paper we argue that for many reionization scenarios low resolution images could be obtained from the expected data. This is because at the later stages of the process one still finds very large pockets of neutral regions in the IGM, reflecting the clustering of the large-scale structure, which stays strong up to scales of \sim 120 comoving Mpc/h (\sim 1 degree). The coherence of the emission on those scales allows us to reach sufficient S/N (\sim 3) so as to obtain reionization 21 cm images. Such images will be extremely valuable for answering many cosmological questions but above all they will be a very powerful tool to test our control of the systematics in the data. The existence of this typical scale (\sim 120 comoving Mpc/h) also argues for designing future EoR experiments, e.g., with SKA, with a field of view of at least 4 degree. △ Less

Submitted 12 June, 2012; v1 submitted 15 May, 2012; originally announced May 2012.

Comments: Replaced with final version (minor changes), 9 figures, 11 pages, accepted for publication in MNRAS

Foreground Removal using FastICA: A Showcase of LOFAR-EoR

Authors: Emma Chapman, Filipe B. Abdalla, Geraint Harker, Vibor Jelić, Panagiotis Labropoulos, Saleem Zaroubi, Michiel A. Brentjens, A. G. de Bruyn, L. V. E. Koopmans

Abstract: We introduce a new implementation of the FastICA algorithm on simulated LOFAR EoR 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… ▽ More We introduce a new implementation of the FastICA algorithm on simulated LOFAR EoR 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 PSF 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. △ Less

Submitted 24 May, 2012; v1 submitted 10 January, 2012; originally announced January 2012.

Comments: 16 pages, 29 figures; Accepted by MNRAS; Updated to match online version. Added section on sensitivity of FastICA method

Journal ref: MNRAS Volume 423, Issue 3, pages 2518-2532, July 2012

A Constraint on the Organization of the Galactic Center Magnetic Field Using Faraday Rotation

Authors: C. J. Law, M. A. Brentjens, G. Novak

Abstract: We present new 6 and 20 cm Very Large Array (VLA) observations of polarized continuum emission of roughly 0.5 square degrees of the Galactic center (GC) region. The 6 cm observations detect diffuse linearly-polarized emission throughout the region with a brightness of roughly 1 mJy per 15"x10" beam. The Faraday rotation measure (RM) toward this polarized emission has structure on degree size scale… ▽ More We present new 6 and 20 cm Very Large Array (VLA) observations of polarized continuum emission of roughly 0.5 square degrees of the Galactic center (GC) region. The 6 cm observations detect diffuse linearly-polarized emission throughout the region with a brightness of roughly 1 mJy per 15"x10" beam. The Faraday rotation measure (RM) toward this polarized emission has structure on degree size scales and ranges from roughly +330 rad/m2 east of the dynamical center (Sgr A) to -880 rad/m2 west of the dynamical center. This RM structure is also seen toward several nonthermal radio filaments, which implies that they have a similar magnetic field orientation and constrains models for their origin. Modeling shows that the RM and its change with Galactic longitude are best explained by the high electron density and strong magnetic field of the GC region. Considering the emissivity of the GC plasma shows that while the absolute RM values are indirect measures of the GC magnetic field, the RM longitude structure directly traces the magnetic field in the central kiloparsec of the Galaxy. Combining this result with previous work reveals a larger RM structure covering the central ~2 degrees of the Galaxy. This RM structure is similar to that proposed by Novak and coworkers, but is shifted roughly 50 pc west of the dynamical center of the Galaxy. If this RM structure originates in the GC region, it shows that the GC magnetic field is organized on ~300 pc size scales. The pattern is consistent with a predominantly poloidal field geometry, pointing from south to north, that is perturbed by the motion of gas in the Galactic disk. △ Less

Submitted 10 February, 2011; originally announced February 2011.

Comments: Accepted to ApJ. emulateapj style, 14 pages, 15 figures

Wide field polarimetry around the Perseus cluster at 350 MHz

Authors: M. A. Brentjens

Abstract: This paper investigates the fascinating diffuse polarization structures at 350 MHz that have previously been tentatively attributed to the Perseus cluster and, more specifically, tries to find out whether the structures are located at (or near) the Perseus cluster, or in the Milky Way. A wide field, eight point Westerbork Synthesis Radio Telescope mosaic of the area around the Perseus cluster was… ▽ More This paper investigates the fascinating diffuse polarization structures at 350 MHz that have previously been tentatively attributed to the Perseus cluster and, more specifically, tries to find out whether the structures are located at (or near) the Perseus cluster, or in the Milky Way. A wide field, eight point Westerbork Synthesis Radio Telescope mosaic of the area around the Perseus cluster was observed in full polarization. The frequency range was 324 to 378 MHz and the resolution of the polarization maps was 2'x3'. The maps were processed using Faraday rotation measure synthesis to counter bandwidth depolarization. The RM-cube covers Faraday depths of -384 to +381 rad m^{-2} in steps of 3 rad m^{-2}. There is emission all over the field at Faraday depths between -50 and +100 rad m^{-2}. All previously observed structures were detected. However, no compelling evidence was found supporting association of those structures with either the Perseus cluster or large scale structure formation gas flows in the Perseus-Pisces super cluster. On the contrary, one of the structures is clearly associated with a Galactic depolarization canal at 1.41 GHz. Another large structure in polarized intensity, as well as Faraday depth at a Faraday depth of +30 rad m^{-2}, coincides with a dark object in WHAM H-alpha maps at a kinematic distance of 0.5 \pm 0.5 kpc. All diffuse polarized emission at 350 MHz towards the Perseus cluster is most likely located within 1 kpc from the Sun. The layers that emit the polarized radiation are less than 40 pc/B|| thick. △ Less

Submitted 3 November, 2010; originally announced November 2010.

Comments: 16 pages, accepted for publication in A&A

Deep multi-frequency rotation measure tomography of the galaxy cluster A2255

Authors: R. F. Pizzo, A. G. de Bruyn, G. Bernardi, M. A. Brentjens

Abstract: We aim to unveil their 3-dimensional geometry of Abell 2255 through WSRT observations at 18, 21, 25, 85, and 200 cm. The polarization images of the cluster were processed through rotation measure (RM) synthesis, producing three final RM cubes. The radio galaxies and the filaments at the edges of the halo are detected in the high-frequency RM cube, obtained by combining the data at 18, 21, and 25 c… ▽ More We aim to unveil their 3-dimensional geometry of Abell 2255 through WSRT observations at 18, 21, 25, 85, and 200 cm. The polarization images of the cluster were processed through rotation measure (RM) synthesis, producing three final RM cubes. The radio galaxies and the filaments at the edges of the halo are detected in the high-frequency RM cube, obtained by combining the data at 18, 21, and 25 cm. Their Faraday spectra show different levels of complexity. The radio galaxies lying near by the cluster center have Faraday spectra with multiple peaks, while those at large distances show only one peak, as do the filaments. Similar RM distributions are observed for the external radio galaxies and for the filaments, with much lower average RM values and RM variance than those found in previous works for the central radio galaxies. The 85 cm RM cube is dominated by the Galactic foreground emission, but it also shows features associated with the cluster. At 2 m, no polarized emission from A2255 nor our Galaxy is detected. The radial trend observed in the RM distributions of the radio galaxies and in the complexity of their Faraday spectra favors the interpretation that the external Faraday screen for all the sources in A2255 is the ICM. Its differential contribution depends on the amount of medium that the radio signal crosses along the line of sight. The filaments should therefore be located at the periphery of the cluster, and their apparent central location comes from projection effects. Their high fractional polarization and morphology suggest that they are relics rather than part of a genuine radio halo. Their inferred large distance from the cluster center and their geometry could argue for an association with large-scale structure (LSS) shocks. △ Less

Submitted 20 August, 2010; originally announced August 2010.

Comments: 24 pages, 19 figures

Power spectrum extraction for redshifted 21-cm epoch of reionization experiments: the LOFAR case

Authors: Geraint Harker, Saleem Zaroubi, Gianni Bernardi, Michiel A. Brentjens, A. G. de Bruyn, Benedetta Ciardi, Vibor Jelic, Leon V. E. Koopmans, Panagiotis Labropoulos, Garrelt Mellema, Andre Offringa, V. N. Pandey, Andreas H. Pawlik, Joop Schaye, Rajat M. Thomas, Sarod Yatawatta

Abstract: One of the aims of the Low Frequency Array (LOFAR) Epoch of Reionization (EoR) project is to measure the power spectrum of variations in the intensity of redshifted 21-cm radiation from the EoR. The sensitivity with which this power spectrum can be estimated depends on the level of thermal noise and sample variance, and also on the systematic errors arising from the extraction process, in partic… ▽ More One of the aims of the Low Frequency Array (LOFAR) Epoch of Reionization (EoR) project is to measure the power spectrum of variations in the intensity of redshifted 21-cm radiation from the EoR. The sensitivity with which this power spectrum can be estimated depends on the level of thermal noise and sample variance, and also on the systematic errors arising from the extraction process, in particular from the subtraction of foreground contamination. We model the extraction process using realistic simulations of the cosmological signal, the foregrounds and noise, and so estimate the sensitivity of the LOFAR EoR experiment to the redshifted 21-cm power spectrum. Detection of emission from the EoR should be possible within 360 hours of observation with a single station beam. Integrating for longer, and synthesizing multiple station beams within the primary (tile) beam, then enables us to extract progressively more accurate estimates of the power at a greater range of scales and redshifts. We discuss different observational strategies which compromise between depth of observation, sky coverage and frequency coverage. A plan in which lower frequencies receive a larger fraction of the time appears to be promising. We also study the nature of the bias which foreground fitting errors induce on the inferred power spectrum, and discuss how to reduce and correct for this bias. The angular and line-of-sight power spectra have different merits in this respect, and we suggest considering them separately in the analysis of LOFAR data. △ Less

Submitted 4 March, 2010; originally announced March 2010.

Comments: 13 pages, 11 figures, 1 table; accepted for publication in MNRAS

Journal ref: Mon. Not. R. Astron. Soc. 405 (2010) 2492-2504

Foregrounds for observations of the cosmological 21 cm line: II. Westerbork observations of the fields around 3C196 and the North Celestial Pole

Authors: G. Bernardi, A. G. de Bruyn, G. Harker, M. A. Brentjens, B. Ciardi, V. Jelić, L. V. E. Koopmans, P. Labropoulos, A. Offringa, V. N. Pandey, J. Schaye, R. M. Thomas, S. Yatawatta, S. Zaroubi

Abstract: In the coming years a new insight into galaxy formation and the thermal history of the Universe is expected to come from the detection of the highly redshifted cosmological 21 cm line. The cosmological 21 cm line signal is buried under Galactic and extragalactic foregrounds which are likely to be a few orders of magnitude brighter. Strategies and techniques for effective subtraction of these foreg… ▽ More In the coming years a new insight into galaxy formation and the thermal history of the Universe is expected to come from the detection of the highly redshifted cosmological 21 cm line. The cosmological 21 cm line signal is buried under Galactic and extragalactic foregrounds which are likely to be a few orders of magnitude brighter. Strategies and techniques for effective subtraction of these foreground sources require a detailed knowledge of their structure in both intensity and polarization on the relevant angular scales of 1-30 arcmin. We present results from observations conducted with the Westerbork telescope in the 140-160 MHz range with 2 arcmin resolution in two fields located at intermediate Galactic latitude, centred around the bright quasar 3C196 and the North Celestial Pole. They were observed with the purpose of characterizing the foreground properties in sky areas where actual observations of the cosmological 21 cm line could be carried out. The polarization data were analysed through the rotation measure synthesis technique. We have computed total intensity and polarization angular power spectra. Total intensity maps were carefully calibrated, reaching a high dynamic range, 150000:1 in the case of the 3C196 field. [abridged] △ Less

Submitted 12 August, 2010; v1 submitted 22 February, 2010; originally announced February 2010.

Comments: 20 pages, 22 figures, accepted for publication in A&A. A version with full resolution figures is available at http://www.astro.rug.nl/~bernardi/NCP_3C196/bernardi.pdf

Foregrounds for observations of the cosmological 21 cm line: I. First Westerbork measurements of Galactic emission at 150 MHz in a low latitude field

Authors: G. Bernardi, A. G. de Bruyn, M. A. Brentjens, B. Ciardi, G. Harker, V. Jelic, L. V. E. Koopmans, P. Labropoulos, A. Offringa, V. N. Pandey, J. Schaye, R. M. Thomas, S. Yatawatta, S. Zaroubi

Abstract: 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. 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 la… ▽ More 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. 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. Those fluctuations have an $rms$ of 14 K. The total intensity power spectrum shows a power--law behaviour down to $\ell \sim 900$ with slope $β^I_\ell = -2.2 \pm 0.3$. The detection of diffuse emission at smaller angular scales is limited by residual point sources. We measured an $rms$ confusion noise of $\sim$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 shows a power--law behaviour down to $\ell \sim 2700$ with slope $β^P_\ell = -1.65 \pm 0.15$. The $rms$ of polarization fluctuations is 7.2 K on 4 arcmin scales. By extrapolating the measured spectrum of total intensity emission, we find a contamination on the cosmological signal of $δT= \sqrt{\ell (\ell+1) C^I_\ell / 2π} \sim 5.7$ K on 5 arcmin scales and a corresponding $rms$ value of $\sim$18.3 K at the same angular scale. The level of the polarization power spectrum is $δT \sim 3.3$ K on 5 arcmin scales. Given its exceptionally bright polarized signal, the Fan region is likely to represent an upper limit on the sky brightness at moderate and high Galactic latitude. △ Less

Submitted 17 July, 2009; v1 submitted 2 April, 2009; originally announced April 2009.

Comments: Minor corrections made to match the final version printed on A&A. A version with high resolution figures is available at http://www.astro.rug.nl/~bernardi/FAN/fan.pdf

Non-parametric foreground subtraction for 21cm epoch of reionization experiments

Authors: Geraint Harker, Saleem Zaroubi, Gianni Bernardi, Michiel A. Brentjens, A. G. de Bruyn, Benedetta Ciardi, Vibor Jelic, Leon V. E. Koopmans, Panagiotis Labropoulos, Garrelt Mellema, Andre Offringa, V. N. Pandey, Joop Schaye, Rajat M. Thomas, Sarod Yatawatta

Abstract: An obstacle to the detection of redshifted 21cm emission from the epoch of reionization (EoR) is the presence of foregrounds which exceed the cosmological signal in intensity by orders of magnitude. We argue that in principle it would be better to fit the foregrounds non-parametrically - allowing the data to determine their shape - rather than selecting some functional form in advance and then f… ▽ More An obstacle to the detection of redshifted 21cm emission from the epoch of reionization (EoR) is the presence of foregrounds which exceed the cosmological signal in intensity by orders of magnitude. We argue that in principle it would be better to fit the foregrounds non-parametrically - allowing the data to determine their shape - rather than selecting some functional form in advance and then fitting its parameters. Non-parametric fits often suffer from other problems, however. We discuss these before suggesting a non-parametric method, Wp smoothing, which seems to avoid some of them. After outlining the principles of Wp smoothing we describe an algorithm used to implement it. We then apply Wp smoothing to a synthetic data cube for the LOFAR EoR experiment. The performance of Wp smoothing, measured by the extent to which it is able to recover the variance of the cosmological signal and to which it avoids leakage of power from the foregrounds, is compared to that of a parametric fit, and to another non-parametric method (smoothing splines). We find that Wp smoothing is superior to smoothing splines for our application, and is competitive with parametric methods even though in the latter case we may choose the functional form of the fit with advance knowledge of the simulated foregrounds. Finally, we discuss how the quality of the fit is affected by the frequency resolution and range, by the characteristics of the cosmological signal and by edge effects. △ Less

Submitted 15 June, 2009; v1 submitted 16 March, 2009; originally announced March 2009.

Comments: 15 pages, 12 figures; lengthened and two figures added, to match version accepted by MNRAS

Journal ref: Mon. Not. Roy. Astron. Soc. 397 (2009) 1138-1152

Detection and extraction of signals from the epoch of reionization using higher order one-point statistics

Authors: Geraint J. A. Harker, Saleem Zaroubi, Rajat M. Thomas, Vibor Jelic, Panagiotis Labropoulos, Garrelt Mellema, Ilian T. Iliev, Gianni Bernardi, Michiel A. Brentjens, A. G. de Bruyn, Benedetta Ciardi, Leon V. E. Koopmans, V. N. Pandey, Andreas H. Pawlik, Joop Schaye, Sarod Yatawatta

Abstract: Detecting redshifted 21cm emission from neutral hydrogen in the early Universe promises to give direct constraints on the epoch of reionization (EoR). It will, though, be very challenging to extract the cosmological signal (CS) from foregrounds and noise which are orders of magnitude larger. Fortunately, the signal has some characteristics which differentiate it from the foregrounds and noise, a… ▽ More Detecting redshifted 21cm emission from neutral hydrogen in the early Universe promises to give direct constraints on the epoch of reionization (EoR). It will, though, be very challenging to extract the cosmological signal (CS) from foregrounds and noise which are orders of magnitude larger. Fortunately, the signal has some characteristics which differentiate it from the foregrounds and noise, and we suggest that using the correct statistics may tease out signatures of reionization. We generate mock datacubes simulating the output of the Low Frequency Array (LOFAR) EoR experiment. These cubes combine realistic models for Galactic and extragalactic foregrounds and the noise with three different simulations of the CS. We fit out the foregrounds, which are smooth in the frequency direction, to produce residual images in each frequency band. We denoise these images and study the skewness of the one-point distribution in the images as a function of frequency. We find that, under sufficiently optimistic assumptions, we can recover the main features of the redshift evolution of the skewness in the 21cm signal. We argue that some of these features - such as a dip at the onset of reionization, followed by a rise towards its later stages - may be generic, and give us a promising route to a statistical detection of reionization. △ Less

Submitted 2 March, 2009; v1 submitted 14 September, 2008; originally announced September 2008.

Comments: 11 pages, 10 figures. One figure added. Other small changes to match accepted version

Journal ref: Mon.Not.Roy.Astron.Soc. 393 (2009) 1449-1458

Fast Large-Scale Reionization Simulations

Authors: Rajat M. Thomas, Saleem Zaroubi, Benedetta Ciardi, Andreas H. Pawlik, Panagiotis Labropoulos, Vibor Jelic, Gianni Bernardi, Michiel A. Brentjens, A. G. de Bruyn, Geraint J. A. Harker, Leon V. E. Koopmans, Garrelt Mellema, V. N. Pandey, Joop Schaye, Sarod Yatawatta

Abstract: We present an efficient method to generate large simulations of the Epoch of Reionization (EoR) without the need for a full 3-dimensional radiative transfer code. Large dark-matter-only simulations are post-processed to produce maps of the redshifted 21cm emission from neutral hydrogen. Dark matter haloes are embedded with sources of radiation whose properties are either based on semi-analytical… ▽ More We present an efficient method to generate large simulations of the Epoch of Reionization (EoR) without the need for a full 3-dimensional radiative transfer code. Large dark-matter-only simulations are post-processed to produce maps of the redshifted 21cm emission from neutral hydrogen. Dark matter haloes are embedded with sources of radiation whose properties are either based on semi-analytical prescriptions or derived from hydrodynamical simulations. These sources could either be stars or power-law sources with varying spectral indices. Assuming spherical symmetry, ionized bubbles are created around these sources, whose radial ionized fraction and temperature profiles are derived from a catalogue of 1-D radiative transfer experiments. In case of overlap of these spheres, photons are conserved by redistributing them around the connected ionized regions corresponding to the spheres. The efficiency with which these maps are created allows us to span the large parameter space typically encountered in reionization simulations. We compare our results with other, more accurate, 3-D radiative transfer simulations and find excellent agreement for the redshifts and the spatial scales of interest to upcoming 21cm experiments. We generate a contiguous observational cube spanning redshift 6 to 12 and use these simulations to study the differences in the reionization histories between stars and quasars. Finally, the signal is convolved with the LOFAR beam response and its effects are analyzed and quantified. Statistics performed on this mock data set shed light on possible observational strategies for LOFAR. △ Less

Submitted 8 September, 2008; originally announced September 2008.

Comments: 18 pages, 21 figures, submitted to MNRAS For high-resolution images follow "http://www.astro.rug.nl/~thomas/eormap.pdf"

Deep Westerbork observations of Abell 2256 at 350 MHz

Authors: M. A. Brentjens

Abstract: Deep polarimetric Westerbork observations of the galaxy cluster Abell 2256 are presented, covering a frequency range of 325-377 MHz. The central halo source has a diameter of the order of 1.2 Mpc (18'), which is somewhat larger than at 1.4 GHz. With alpha=-1.61+/-0.04, the halo spectrum between 1.4 GHz and 22.25 MHz is less steep than previously thought. The centre of the ultra steep spectrum so… ▽ More Deep polarimetric Westerbork observations of the galaxy cluster Abell 2256 are presented, covering a frequency range of 325-377 MHz. The central halo source has a diameter of the order of 1.2 Mpc (18'), which is somewhat larger than at 1.4 GHz. With alpha=-1.61+/-0.04, the halo spectrum between 1.4 GHz and 22.25 MHz is less steep than previously thought. The centre of the ultra steep spectrum source in the eastern part of the cluster exhibits a spectral break near 400 MHz. It is estimated to be at least 51 million years old, but possibly older than 125 million years. A final measurement requires observations in the 10-150 MHz range. It remains uncertain whether the source is a radio tail of Fabricant galaxy 122, situated in the northeastern tip of the source. Faraday rotation measure synthesis revealed no polarized flux at all in the cluster. The polarization fraction of the brightest parts of the relic area is less than 1%. The RM-synthesis nevertheless revealed 9 polarized sources in the field enabling an accurate measurement of the Galactic Faraday rotation (-33+/-2 rad m^{-2} in front of the relic). Based on its depolarization on longer wavelengths, the line-of-sight magnetic field in relic filament G is estimated to be between 0.02 and 2 uG. A value of 0.2 uG appears most reasonable given the currently available data. △ Less

Submitted 28 July, 2008; originally announced July 2008.

Comments: 16 pages, 16 figures, accepted by A&A

Foreground simulations for the LOFAR - Epoch of Reionization Experiment

Authors: Vibor Jelic, Saleem Zaroubi, Panagiotis Labropoulos, Rajat M. Thomas, Gianni Bernardi, Michiel A. Brentjens, Ger de Bruyn, Benedetta Ciardi, Geraint Harker, Leon V. E. Koopmans, Vishambhar N. Pandey, Joop Schaye, Sarod Yatawatta

Abstract: 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… ▽ More 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 datacubes and discuss their possible statistical effects on the data. The foreground maps are produced assuming 5 deg x 5 deg 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 300 hours 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 idealised 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. △ Less

Submitted 1 July, 2008; v1 submitted 8 April, 2008; originally announced April 2008.

Comments: 18 figures, 3 tables, accepted to be published in MNRAS

Journal ref: 2008MNRAS.389.1319J

A very brief description of LOFAR - the Low Frequency Array

Authors: H. Falcke, M. P. van Haarlem, A. G. de Bruyn, R. Braun, H. J. A. Röttgering, B. Stappers, W. H. W. M. Boland, H. R. Butcher, E. J. de Geus, L. Koopmans, R. Fender, J. Kuijpers, G. K. Miley, R. T. Schilizzi, C. Vogt, R. A. M. J. Wijers, M. Wise, W. N. Brouw, J. P. Hamaker, J. E. Noordam, T. Oosterloo, L. Bähren, M. A. Brentjens, S. J. Wijnholds, J. D. Bregman , et al. (6 additional authors not shown)

Abstract: LOFAR (Low Frequency Array) is an innovative radio telescope optimized for the frequency range 30-240 MHz. The telescope is realized as a phased aperture array without any moving parts. Digital beam forming allows the telescope to point to any part of the sky within a second. Transient buffering makes retrospective imaging of explosive short-term events possible. The scientific focus of LOFAR wi… ▽ More LOFAR (Low Frequency Array) is an innovative radio telescope optimized for the frequency range 30-240 MHz. The telescope is realized as a phased aperture array without any moving parts. Digital beam forming allows the telescope to point to any part of the sky within a second. Transient buffering makes retrospective imaging of explosive short-term events possible. The scientific focus of LOFAR will initially be on four key science projects (KSPs): 1) detection of the formation of the very first stars and galaxies in the universe during the so-called epoch of reionization by measuring the power spectrum of the neutral hydrogen 21-cm line (Shaver et al. 1999) on the ~5' scale; 2) low-frequency surveys of the sky with of order $10^8$ expected new sources; 3) all-sky monitoring and detection of transient radio sources such as gamma-ray bursts, x-ray binaries, and exo-planets (Farrell et al. 2004); and 4) radio detection of ultra-high energy cosmic rays and neutrinos (Falcke & Gorham 2003) allowing for the first time access to particles beyond 10^21 eV (Scholten et al. 2006). Apart from the KSPs open access for smaller projects is also planned. Here we give a brief description of the telescope. △ Less

Submitted 21 October, 2006; originally announced October 2006.

Comments: 2 pages, IAU GA 2006, Highlights of Astronomy, Volume 14, K.A. van der Hucht, ed

LOFAR - Opening up a new window on the Universe

Authors: H. J. A. Rottgering, R. Braun, P. D. Barthel, M. P. van Haarlem, G. K. Miley, R. Morganti, I. Snellen, H. Falcke, A. G. de Bruyn, R. B. Stappers, W. H. W. M. Boland, H. R. Butcher, E. J. de Geus, L. Koopmans, R. Fender, J. Kuijpers, R. T. Schilizzi, C. Vogt, R. A. M. J. Wijers, M. Wise, W. N. Brouw, J. P. Hamaker, J. E. Noordam, T. Oosterloo, L. Bahren , et al. (9 additional authors not shown)

Abstract: LOFAR, the Low Frequency Array, is a next-generation radio telescope that is being built in Northern Europe and expected to be fully operational at the end of this decade. It will operate at frequencies from 15 to 240 MHz (corresponding to wavelengths of 20 to 1.2 m). Its superb sensitivity, high angular resolution, large field of view and flexible spectroscopic capabilities will represent a dra… ▽ More LOFAR, the Low Frequency Array, is a next-generation radio telescope that is being built in Northern Europe and expected to be fully operational at the end of this decade. It will operate at frequencies from 15 to 240 MHz (corresponding to wavelengths of 20 to 1.2 m). Its superb sensitivity, high angular resolution, large field of view and flexible spectroscopic capabilities will represent a dramatic improvement over previous facilities at these wavelengths. As such, LOFAR will carry out a broad range of fundamental astrophysical studies. The design of LOFAR has been driven by four fundamental astrophysical applications: (i) The Epoch of Reionisation, (ii) Extragalactic Surveys and their exploitation to study the formation and evolution of clusters, galaxies and black holes, (iii) Transient Sources and their association with high energy objects such as gamma ray bursts, and (iv) Cosmic Ray showers and their exploitation to study the origin of ultra-high energy cosmic rays. In this conference the foreseen LOFAR work on the epoch of reionisation has been covered by de Bruyn and on cosmic ray showers by Falcke. During this contribution we will first present the LOFAR project with an emphasis on the challenges faced when carrying out sensitive imaging at low radio frequencies. Subsequently, we will discuss LOFAR's capabilities to survey the low-frequency radio sky. Main aims for the planned surveys are studies of z>6 radio galaxies, diffuse emission associated with distant clusters and starbursting galaxies at z>2. △ Less

Submitted 29 January, 2007; v1 submitted 19 October, 2006; originally announced October 2006.

Comments: 11 pages, 5 figures, to appear in the proceedings of the conference "Cosmology, galaxy formation and astroparticle physics on the pathway to the SKA", Oxford, April 10-12 2006

Diffuse polarized emission associated with the Perseus cluster

Authors: A. G. de Bruyn, M. A. Brentjens

Abstract: We report on full-polarization radio observations of the Perseus cluster (Abell 426) using the Westerbork Synthesis Radio Telescope (WSRT) at wavelengths from 81-95 cm. We have employed a novel technique, Rotation Measure synthesis (Brentjens and de Bruyn, 2005) to unravel the polarization properties of the emission across the full field of view and detect polarized emission over a wide range of… ▽ More We report on full-polarization radio observations of the Perseus cluster (Abell 426) using the Westerbork Synthesis Radio Telescope (WSRT) at wavelengths from 81-95 cm. We have employed a novel technique, Rotation Measure synthesis (Brentjens and de Bruyn, 2005) to unravel the polarization properties of the emission across the full field of view and detect polarized emission over a wide range of RM from about 0 to 90 rad m^-2. The low RM emission is associated with our Galaxy, while the high RM emission is associated with the Perseus cluster. The latter reaches typical surface brightness levels of 0.5-1 mJy per beam and must be rather highly polarized. Most of the peripheral polarized emission appears too bright, by about 1-2 orders of magnitude, to be explainable as Thomson scattered emission of the central radio source off the thermal electrons in the cluster. The bulk of the emission associated with the Perseus cluster is probably related to buoyant bubbles of relativistic plasma, probably relics from still active or now dormant AGN within the cluster. A lenticular shaped structure measuring 0.5-1 Mpc is strikingly similar to the structures predicted by Ensslin et al. (1998). At the western edge of the cluster, we detect very long, linear structures that may be related to shocks caused by infall of gas into the Perseus cluster. △ Less

Submitted 2 November, 2005; v1 submitted 14 July, 2005; originally announced July 2005.

Comments: 18 pages, 17 figures, accepted by A&A, corrected small typo, added reference

Journal ref: Astron.Astrophys.441:931-947, 2005

Faraday Rotation Measure Synthesis

Authors: M. A. Brentjens, A. G. de Bruyn

Abstract: We extend the rotation measure work of Burn (1966) to the cases of limited sampling of lambda squared space and non-constant emission spectra. We introduce the rotation measure transfer function (RMTF), which is an excellent predictor of n-pi ambiguity problems with the lambda squared coverage. Rotation measure synthesis can be implemented very efficiently on modern computers. Because the analys… ▽ More We extend the rotation measure work of Burn (1966) to the cases of limited sampling of lambda squared space and non-constant emission spectra. We introduce the rotation measure transfer function (RMTF), which is an excellent predictor of n-pi ambiguity problems with the lambda squared coverage. Rotation measure synthesis can be implemented very efficiently on modern computers. Because the analysis is easily applied to wide fields, one can conduct very fast RM surveys of weak spatially extended sources. Difficult situations, for example multiple sources along the line of sight, are easily detected and transparently handled. Under certain conditions, it is even possible to recover the emission as a function of Faraday depth within a single cloud of ionized gas. Rotation measure synthesis has already been successful in discovering widespread, weak, polarized emission associated with the Perseus cluster (De Bruyn and Brentjens, 2005). In simple, high signal to noise situations it is as good as traditional linear fits to polarization angle versus lambda squared plots. However, when the situation is more complex or very weak polarized emission at high rotation measures is expected, it is the only viable option. △ Less

Submitted 2 November, 2005; v1 submitted 14 July, 2005; originally announced July 2005.

Comments: 17 pages, 14 figures, accepted by A&A, added references, corrected typos

Journal ref: Astron.Astrophys.441:1217-1228, 2005