Source author record

A. G. de Bruyn

A. G. de Bruyn appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

57works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

57 published item(s)

preprint2016arXiv

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

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.

preprint2016arXiv

Probing Ionospheric Structures using the LOFAR radio telescope

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.

preprint2016arXiv

Systematic biases in low frequency radio interferometric data due to calibration: the LOFAR EoR case

The redshifted 21 cm line of neutral hydrogen is a promising probe of the Epoch of Reionization (EoR). However, its detection requires a thorough understanding and control of the systematic errors. We study two systematic biases observed in the LOFAR EoR residual data after calibration and subtraction of bright discrete foreground sources. The first effect is a suppression in the diffuse foregrounds, which could potentially mean a suppression of the 21 cm signal. The second effect is an excess of noise beyond the thermal noise. The excess noise shows fluctuations on small frequency scales, and hence it can not be easily removed by foreground removal or avoidance methods. Our analysis suggests that sidelobes of residual sources due to the chromatic point spread function and ionospheric scintillation can not be the dominant causes of the excess noise. Rather, both the suppression of diffuse foregrounds and the excess noise can occur due to calibration with an incomplete sky model containing predominantly bright discrete sources. We show that calibrating only on bright sources can cause suppression of other signals and introduce an excess noise in the data. The levels of the suppression and excess noise depend on the relative flux of sources which are not included in the model with respect to the flux of modeled sources. We discuss possible solutions such as using only long baselines to calibrate the interferometric gain solutions as well as simultaneous multi-frequency calibration along with their benefits and shortcomings.

preprint2015arXiv

Galactic interstellar filaments as probed by LOFAR and Planck

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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}$.

preprint2015arXiv

Realistic Simulations of the Galactic Polarized Foreground: Consequences for 21-cm Reionization Detection Experiments

Experiments designed to measure the redshifted 21~cm line from the Epoch of Reionization (EoR) are challenged by strong astrophysical foreground contamination, ionospheric distortions, complex instrumental response and other different types of noise (e.g. radio frequency interference). The astrophysical foregrounds are dominated by diffuse synchrotron emission from our Galaxy. Here we present a simulation of the Galactic emission used as a foreground module for the LOFAR- EoR key science project end-to-end simulations. The simulation produces total and polarized intensity over $10^\circ \times 10^\circ$ maps of the Galactic synchrotron and free-free emission, including all observed characteristics of the emission: spatial fluctuations of amplitude and spectral index of the synchrotron emission, together with Faraday rotation effects. The importance of these simulations arise from the fact that the Galactic polarized emission could behave in a manner similar to the EoR signal along the frequency direction. As a consequence, an improper instrumental calibration will give rise to leakages of the polarized to the total signal and mask the desired EoR signal. In this paper we address this for the first time through realistic simulations.

preprint2015arXiv

Simulating the 21cm forest detectable with LOFAR and SKA in the spectra of high-z GRBs

We investigate the feasibility of detecting 21cm absorption features in the afterglow spectra of high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic reionization, together with the instrumental characteristics of the LOw Frequency ARray (LOFAR). We find that absorption features could be marginally (with a S/N larger than a few) detected by LOFAR at z>7 if the GRB originated from PopIII stars, while the detection would be easier if the noise were reduced by one order of magnitude, i.e. similar to what is expected for the first phase of the Square Kilometer Array (SKA1-low). On the other hand, more standard GRBs are too dim to be detected even with ten times the sensitivity of SKA1-low, and only in the most optimistic case can a S/N larger than a few be reached at z>9.

preprint2015arXiv

The LOFAR Multifrequency Snapshot Sky Survey (MSSS) I. Survey description and first results

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

preprint2015arXiv

Wide-field LOFAR imaging of the field around the double-double radio galaxy B1834+620: A fresh view on a restarted AGN and doubeltjes

The existence of double-double radio galaxies (DDRGs) is evidence for recurrent jet activity in AGN, as expected from standard accretion models. A detailed study of these rare sources provides new perspectives for investigating the AGN duty cycle, AGN-galaxy feedback, and accretion mechanisms. Large catalogues of radio sources provide statistical information about the evolution of the radio-loud AGN population out to high redshifts. Using wide-field imaging with the LOFAR telescope, we study both a well-known DDRG as well as a large number of radio sources in the field of view. We present a high resolution image of the DDRG B1834+620 obtained at 144 MHz using LOFAR commissioning data. Our image covers about 100 square degrees and contains over 1000 sources. The four components of the DDRG B1834+620 have been resolved for the first time at 144 MHz. Inner lobes were found to point towards the direction of the outer lobes, unlike standard FR~II sources. Polarized emission was detected in the northern outer lobe. The high spatial resolution allows the identification of a large number of small double-lobed radio sources; roughly 10% of all sources in the field are doubles with a separation smaller than 1 arcmin. The spectral fit of the four components is consistent with a scenario in which the outer lobes are still active or the jets recently switched off, while emission of the inner lobes is the result of a mix-up of new and old jet activity. From the presence of the newly extended features in the inner lobes of the DDRG, we can infer that the mechanism responsible for their formation is the bow shock that is driven by the newly launched jet. We find that the density of the small doubles exceeds the density of FR-II sources with similar properties at 1.4 GHz, but this difference becomes smaller for low flux densities.

preprint2014arXiv

A Comprehensive Radio View of the Extremely Bright Gamma-Ray Burst 130427A

GRB 130427A was extremely bright as a result of occurring at low redshift whilst the energetics were more typical of high-redshift gamma-ray bursts (GRBs). We collected well-sampled light curves at 1.4 and 4.8 GHz of GRB 130427A with the Westerbork Synthesis Radio Telescope (WSRT); and we obtained its most accurate position with the European Very Long Baseline Interferometry Network (EVN). Our flux density measurements are combined with all the data available at radio, optical and X-ray frequencies to perform broadband modeling in the framework of a reverse-forward shock model and a two-component jet model, and we discuss the implications and limitations of both models. The low density inferred from the modeling implies that the GRB 130427A progenitor is either a very low-metallicity Wolf-Rayet star, or a rapidly rotating, low-metallicity O star. We also find that the fraction of the energy in electrons is evolving over time, and that the fraction of electrons participating in a relativistic power-law energy distribution is less than 15%. We observed intraday variability during the earliest WSRT observations, and the source sizes inferred from our modeling are consistent with this variability being due to interstellar scintillation effects. Finally, we present and discuss our limits on the linear and circular polarization, which are among the deepest limits of GRB radio polarization to date.

preprint2014arXiv

Calibration artefacts in radio interferometry. I. Ghost sources in WSRT data

This work investigates a particular class of artefacts, or ghost sources, in radio interferometric images. Earlier observations with (and simulations of) the Westerbork Synthesis Radio Telescope (WSRT) suggested that these were due to calibration with incomplete sky models. A theoretical framework is derived that validates this suggestion, and provides predictions of ghost formation in a two-source scenario. The predictions are found to accurately match the result of simulations, and qualitatively reproduce the ghosts previously seen in observational data. The theory also provides explanations for many previously puzzling features of these artefacts (regular geometry, PSF-like sidelobes, seeming independence on model flux), and shows that the observed phenomenon of flux suppression affecting unmodelled sources is due to the same mechanism. We demonstrate that this ghost formation mechanism is a fundamental feature of calibration, and exhibits a particularly strong and localized signature due to array redundancy. To some extent this mechanism will affect all observations (including those with non-redundant arrays), though in most cases the ghosts remain hidden below the noise or masked by other instrumental artefacts. The implications of such errors on future deep observations are discussed.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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%.

preprint2014arXiv

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

We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<ν<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$).

preprint2014arXiv

The intra-hour variable quasar J1819$+$3845: 13-year evolution, jet polarization structure and interstellar scattering screen properties

We examine the long-term evolution of the intra-hour variable quasar, J1819+3845, whose variations have been attributed to interstellar scintillation by extremely local turbulent plasma, located only 1-3pc from Earth. The variations in this source ceased some time between June 2006 and February 2007. The evolution of the source spectrum and the long-term lightcurve, and the persistent compactness of the source VLBI structure indicates that the cessation of rapid variability was associated with the passage of the scattering material out of the line of sight to the quasar. We present an analysis of the linear polarization variations and their relation to total intensity variations. The proper motion of polarized features in the quasar jet is found to be subluminal. Systematic time delays between Stokes I, Q and U, in combination with the structure of the source obtained from 8.4GHz VLBI data, confirm the estimate of the screen distance: 1-2pc, making the screen one of the nearest objects to the Solar System. We determine the physical properties of this scattering material. The electron density in the scattering region is extremely high with respect to the warm ionized ISM, with an estimated density of $n_e \sim 97 \, l_0^{1/3} {ΔL}_{100}^{-1/2}$cm$^{-3}$, where $l_0$ is the outer scale of the turbulence in AU and $ΔL = 100 ΔL_{100}$ AU is the depth of the scattering region. If this plasma is in pressure balance with the local magnetic field, one expects a ~2 rad/m^2 rotation measure change associated with the passage of this material past the quasar. We examine the rotation measures of sources and the diffuse polarized emission in the surrounding region. We place a limit of 10 rad/m^2 on the RM change. The variability of sources near J1819+3845 is used to deduce that the screen must therefore be either very small (~100 AU) or patchy.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

LEDDB: LOFAR Epoch of Reionization Diagnostic Database

One of the key science projects of the Low-Frequency Array (LOFAR) is the detection of the cosmological signal coming from the Epoch of Reionization (EoR). Here we present the LOFAR EoR Diagnostic Database (LEDDB) that is used in the storage, management, processing and analysis of the LOFAR EoR observations. It stores referencing information of the observations and diagnostic parameters extracted from their calibration. This stored data is used to ease the pipeline processing, monitor the performance of the telescope and visualize the diagnostic parameters which facilitates the analysis of the several contamination effects on the signals. It is implemented with PostgreSQL and accessed through the psycopg2 python module. We have developed a very flexible query engine, which is used by a web user interface to access the database, and a very extensive set of tools for the visualization of the diagnostic parameters through all their multiple dimensions.

preprint2013arXiv

LOFAR imaging capabilities and system sensitivity

This document describes the general astronomical capabilities of the LOw Frequency ARray (LOFAR). The frequency range covered by LOFAR is split into two bands denoted as low band (LB, 10 - 80 MHz) and high band (HB, 120 - 240 MHz). LOFAR stations are spread over a 100 km sized region in the northern part of the Netherlands. In addition to the Dutch stations there will be European stations providing baselines between 200 km and 1000 km. Most of the results in this document, however, are limited to the Dutch array. In section 2 the LBA and HBA station configurations are reviewed. Section 3 reviews the imaging capabilities as station Full Width Half Maximum, station Field of View, and array resolution. In section 4 the system sensitivity is considered and thermal noise levels in LOFAR images are given. Finally, in section 5 some cautionary notes are collected. We advise you to read these notes carefully before using the numbers that are presented in this report.

preprint2013arXiv

LOFAR insights into the epoch of reionization from the cross power spectrum of 21cm emission and galaxies

Using a combination of N-body simulations, semi-analytic models and radiative transfer calculations, we have estimated the theoretical cross power spectrum between galaxies and the 21cm emission from neutral hydrogen during the epoch of reionization. In accordance with previous studies, we find that the 21cm emission is initially correlated with halos on large scales (> 30 Mpc), anti-correlated on intermediate (~ 5 Mpc), and uncorrelated on small (< 3 Mpc) scales. This picture quickly changes as reionization proceeds and the two fields become anti-correlated on large scales. The normalization of the cross power spectrum can be used to set constraints on the average neutral fraction in the intergalactic medium and its shape can be a tool to study the topology of reionization. When we apply a drop-out technique to select galaxies and add to the 21cm signal the noise expected from the LOFAR telescope, we find that while the normalization of the cross power spectrum remains a useful tool for probing reionization, its shape becomes too noisy to be informative. On the other hand, for a Lyalpha Emitter (LAE) survey both the normalization and the shape of the cross power spectrum are suitable probes of reionization. A closer look at a specific planned LAE observing program using Subaru Hyper-Suprime Cam reveals concerns about the strength of the 21cm signal at the planned redshifts. If the ionized fraction at z ~ 7 is lower that the one estimated here, then using the cross power spectrum may be a useful exercise given that at higher redshifts and neutral fractions it is able to distinguish between two toy models with different topologies.

preprint2013arXiv

LOFAR: The LOw-Frequency ARray

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

preprint2013arXiv

Probing reionization with LOFAR using 21-cm redshift space distortions

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

preprint2013arXiv

The brightness and spatial distributions of terrestrial radio sources

Faint undetected sources of radio-frequency interference (RFI) might become visible in long radio observations when they are consistently present over time. Thereby, they might obstruct the detection of the weak astronomical signals of interest. This issue is especially important for Epoch of Reionisation (EoR) projects that try to detect the faint redshifted HI signals from the time of the earliest structures in the Universe. We explore the RFI situation at 30-163 MHz by studying brightness histograms of visibility data observed with LOFAR, similar to radio-source-count analyses that are used in cosmology. An empirical RFI distribution model is derived that allows the simulation of RFI in radio observations. The brightness histograms show an RFI distribution that follows a power-law distribution with an estimated exponent around -1.5. With several assumptions, this can be explained with a uniform distribution of terrestrial radio sources whose radiation follows existing propagation models. Extrapolation of the power law implies that the current LOFAR EoR observations should be severely RFI limited if the strength of RFI sources remains strong after time integration. This is in contrast with actual observations, which almost reach the thermal noise and are thought not to be limited by RFI. Therefore, we conclude that it is unlikely that there are undetected RFI sources that will become visible in long observations. Consequently, there is no indication that RFI will prevent an EoR detection with LOFAR.

preprint2013arXiv

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

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

preprint2013arXiv

The spin temperature of high-redshift damped Lyman-$α$ systems

We report results from a programme aimed at investigating the temperature of neutral gas in high-redshift damped Lyman-$α$ absorbers (DLAs). This involved (1) HI 21cm absorption studies of a large DLA sample, (2) VLBI studies to measure the low-frequency quasar core fractions, and (3) optical/ultraviolet spectroscopy to determine DLA metallicities and velocity widths. Including literature data, our sample consists of 37 DLAs with estimates of the spin temperature $T_s$ and the covering factor. We find a strong $4σ$) difference between the $T_s$ distributions in high-z (z>2.4) and low-z (z<2.4) DLA samples. The high-z sample contains more systems with high $T_s$ values, $\gtrsim 1000$ K. The $T_s$ distributions in DLAs and the Galaxy are also clearly (~$6σ$) different, with more high-$T_s$ sightlines in DLAs than in the Milky Way. The high $T_s$ values in the high-z DLAs of our sample arise due to low fractions of the cold neutral medium. For 29 DLAs with metallicity [Z/H] estimates, we confirm the presence of an anti-correlation between $T_s$ and [Z/H], at $3.5σ$ significance via a non-parametric Kendall-tau test. This result was obtained with the assumption that the DLA covering factor is equal to the core fraction. Monte Carlo simulations show that the significance of the result is only marginally decreased if the covering factor and the core fraction are uncorrelated, or if there is a random error in the inferred covering factor. We also find evidence for redshift evolution in DLA $T_s$ values even for the z>1 sub-sample. Since z>1 DLAs have angular diameter distances comparable to or larger than those of the background quasars, they have similar efficiency in covering the quasars. Low covering factors in high-z DLAs thus cannot account for the observed redshift evolution in spin temperatures. (Abstract abridged.)

preprint2012arXiv

Foreground Removal using FastICA: A Showcase of LOFAR-EoR

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.

preprint2012arXiv

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

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.

preprint2012arXiv

Post-correlation filtering techniques for off-axis source and RFI removal

Techniques to improve the data quality of interferometric radio observations are considered. Fundaments of fringe frequencies in the uv-plane are discussed and filters are used to attenuate radio-frequency interference (RFI) and off-axis sources. Several new applications of filters are introduced and tested. A low-pass filter in time and frequency direction on single baseline data is successfully used to lower the noise in the area of interest and to remove sidelobes coming from unmodelled off-axis sources and RFI. Related side effects of data integration, averaging and gridding are analysed, and shown to be able to cause ghosts and an increase in noise, especially when using long baselines or interferometric elements that have a large field of view. A novel projected fringe low-pass filter is shown to be potentially useful for first order source separation. Initial tests show that the filters can be several factors faster compared to common source separation techniques such as peeling and a variant of peeling that is currently being tested on LOFAR observations called "demixed peeling". Further testing is required to support the performance of the filters.

preprint2012arXiv

Prospects for detecting the 21cm forest from the diffuse intergalactic medium with LOFAR

We discuss the feasibility of the detection of the 21cm forest in the diffuse IGM with the radio telescope LOFAR. The optical depth to the 21cm line has been derived using simulations of reionization which include detailed radiative transfer of ionizing photons. We find that the spectra from reionization models with similar total comoving hydrogen ionizing emissivity but different frequency distribution look remarkably similar. Thus, unless the reionization histories are very different from each other (e.g. a predominance of UV vs. x-ray heating) we do not expect to distinguish them by means of observations of the 21cm forest. Because the presence of a strong x-ray background would make the detection of 21cm line absorption impossible, the lack of absorption could be used as a probe of the presence/intensity of the x-ray background and the thermal history of the universe. Along a random line of sight LOFAR could detect a global suppression of the spectrum from z>12, when the IGM is still mostly neutral and cold, in contrast with the more well-defined, albeit broad, absorption features visible at lower redshift. Sharp, strong absorption features associated with rare, high density pockets of gas could be detected also at z~7 along preferential lines of sight.

preprint2012arXiv

The LOFAR radio environment

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.

preprint2011arXiv

Early radio and X-ray observations of the youngest nearby type Ia supernova PTF11kly (SN 2011fe)

On August 24 (UT) the Palomar Transient Factory (PTF) discovered PTF11kly (SN 2011fe), the youngest and most nearby type Ia supernova (SN Ia) in decades. We followed this event up in the radio (centimeter and millimeter bands) and X-ray bands, starting about a day after the estimated explosion time. We present our analysis of the radio and X-ray observations, yielding the tightest constraints yet placed on the pre-explosion mass-loss rate from the progenitor system of this supernova. We find a robust limit of dM/dt<10^-8 (w/100 km/s) [M_solar/yr] from sensitive X-ray non-detections, as well as a similar limit from radio data, which depends, however, on assumptions about microphysical parameters. We discuss our results in the context of single-degenerate models for SNe Ia and find that our observations modestly disfavor symbiotic progenitor models involving a red giant donor, but cannot constrain systems accreting from main-sequence or sub-giant stars, including the popular supersoft channel. In view of the proximity of PTF11kly and the sensitivity of our prompt observations we would have to wait for a long time (decade or longer) in order to more meaningfully probe the circumstellar matter of Ia supernovae.

preprint2011arXiv

MeqTrees and direction-dependent effects

Direction-dependent effects (DDEs) represent a major challenge both for calibration of new radio telescopes, and for advancing the state of the art of existing ones. Various approaches to the problem are currently in development. This paper describes several WSRT observations where a pointing error -- commonly thought to be a major contributor to DDEs at higher frequencies -- was deliberately introduced. These observations were reduced using the MeqTrees package, using two different approaches to correct for DDEs: differential gain solutions, and a direct solution for pointing errors. The results of this are analysed and compared.

preprint2011arXiv

Observing pulsars and fast transients with LOFAR

Low frequency radio waves, while challenging to observe, are a rich source of information about pulsars. The LOw Frequency ARray (LOFAR) is a new radio interferometer operating in the lowest 4 octaves of the ionospheric "radio window": 10-240MHz, that will greatly facilitate observing pulsars at low radio frequencies. Through the huge collecting area, long baselines, and flexible digital hardware, it is expected that LOFAR will revolutionize radio astronomy at the lowest frequencies visible from Earth. LOFAR is a next-generation radio telescope and a pathfinder to the Square Kilometre Array (SKA), in that it incorporates advanced multi-beaming techniques between thousands of individual elements. We discuss the motivation for low-frequency pulsar observations in general and the potential of LOFAR in addressing these science goals. We present LOFAR as it is designed to perform high-time-resolution observations of pulsars and other fast transients, and outline the various relevant observing modes and data reduction pipelines that are already or will soon be implemented to facilitate these observations. A number of results obtained from commissioning observations are presented to demonstrate the exciting potential of the telescope. This paper outlines the case for low frequency pulsar observations and is also intended to serve as a reference for upcoming pulsar/fast transient science papers with LOFAR.

preprint2011arXiv

Radio Interferometric Calibration Using The SAGE Algorithm

The aim of the new generation of radio synthesis arrays such as LOFAR and SKA is to achieve much higher sensitivity, resolution and frequency coverage than what is available now, especially at low frequencies. To accomplish this goal, the accuracy of the calibration techniques used is of considerable importance. Moreover, since these telescopes produce huge amounts of data, speed of convergence of calibration is a major bottleneck. The errors in calibration are due to system noise (sky and instrumental) as well as the estimation errors introduced by the calibration technique itself, which we call solver noise. We define solver noise as the distance between the optimal solution (the true value of the unknowns, uncorrupted by the system noise) and the solution obtained by calibration. We present the Space Alternating Generalized Expectation Maximization (SAGE) calibration technique, which is a modification of the Expectation Maximization algorithm, and compare its performance with the traditional Least Squares calibration based on the level of solver noise introduced by each technique. For this purpose, we develop statistical methods that use the calibrated solutions to estimate the level of solver noise. The SAGE calibration algorithm yields very promising results both in terms of accuracy and speed of convergence. The comparison approaches we adopt introduce a new framework for assessing the performance of different calibration schemes.

preprint2010arXiv

A LOFAR RFI detection pipeline and its first results

Radio astronomy is entering a new era with new and future radio observatories such as the Low Frequency Array and the Square Kilometer Array. We describe in detail an automated flagging pipeline and evaluate its performance. With only a fraction of the computational cost of correlation and its use of the previously introduced SumThreshold method, it is found to be both fast and unrivalled in its high accuracy. The LOFAR radio environment is analysed with the help of this pipeline. The high time and spectral resolution of LOFAR have resulted in an observatory where only a few percent of the data is lost due to RFI.

preprint2010arXiv

Constraints on the flux of Ultra-High Energy neutrinos from WSRT observations

Ultra-high energy (UHE) neutrinos and cosmic rays initiate particle cascades underneath the Moon's surface. These cascades have a negative charge excess and radiate Cherenkov radio emission in a process known as the Askaryan effect. The optimal frequency window for observation of these pulses with radio telescopes on the Earth is around 150 MHz. By observing the Moon with the Westerbork Synthesis Radio Telescope array we are able to set a new limit on the UHE neutrino flux. The PuMa II backend is used to monitor the Moon in 4 frequency bands between 113 and 175 MHz with a sampling frequency of 40 MHz. The narrowband radio interference is digitally filtered out and the dispersive effect of the Earth's ionosphere is compensated for. A trigger system is implemented to search for short pulses. By inserting simulated pulses in the raw data, the detection efficiency for pulses of various strength is calculated. With 47.6 hours of observation time, we are able to set a limit on the UHE neutrino flux. This new limit is an order of magnitude lower than existing limits. In the near future, the digital radio array LOFAR will be used to achieve an even lower limit.

preprint2010arXiv

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

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.

preprint2010arXiv

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

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]

preprint2010arXiv

Post-correlation radio frequency interference classification methods

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

preprint2010arXiv

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

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.

preprint2009arXiv

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

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.

preprint2009arXiv

Improved flux limits for neutrinos with energies above 10$^{22}$ eV from observations with the Westerbork Synthesis Radio Telescope

Particle cascades initiated by ultra-high energy (UHE) neutrinos in the lunar regolith will emit an electromagnetic pulse with a time duration of the order of nano seconds through a process known as the Askaryan effect. It has been shown that in an observing window around 150 MHz there is a maximum chance for detecting this radiation with radio telescopes commonly used in astronomy. In 50 hours of observation time with the Westerbork Synthesis Radio Telescope array we have set a new limit on the flux of neutrinos, summed over all flavors, with energies in excess of $4\times10^{22}$ eV.

preprint2009arXiv

Radio spectral study of the cluster of galaxies Abell 2255

Spectral index studies of halos, relics, and radio galaxies provide useful information on their origin and connection with merger processes. We present WSRT multi-wavelength observations of the galaxy cluster Abell 2255 at 25 cm, 85 cm, and 2 m. The spectral index images allowed us to study the integrated spectrum of halo and relic and to investigate the physical properties of the Beaver head-tail radio galaxy belonging to the cluster. In the radio halo, the spectral index is steeper at the center and flatter at the locations of the radio filaments, clearly detected at 25 cm. In the relics, the spectral index flattens, moving away from the cluster center. For the Beaver radio galaxy, the spectrum severely steepens from the head towards the end of the tail, because of the energy losses suffered by the relativistic particles. In the 2 m map, which is the first high-sensitivity image presented in the literature at such a long wavelength, a new Mpc-size emission region is detected between the known radio halo and the NW relic. Not detecting this feature in the more sensitive 85 cm observations implies that it must have a very steep spectrum (alpha <= -2.6). The observational properties of the radio halo suggest that either we are looking at a superposition of different structures (filaments in the foreground plus real halo in the background) seen in projection across the cluster center or that the halo is intrinsically peculiar. The newly detected extended region to the NW of the halo could be considered as an asymmetric extension of the halo itself. However, since radio halos are known in the literature as structures showing a regular morphology, the new feature could represent the first example of steep Mpc-size diffuse structures (MDS), detected around clusters at very low frequencies.

preprint2008arXiv

WSRT Faraday tomography of the Galactic ISM at λ~ 0.86 m

We investigate the properties of the Galactic ISM by applying Faraday tomography to a radio polarization data set in the direction of the Galactic anti-centre. We address the problem of missing large-scale structure in our data, and show that this does not play an important role for the results we present. The main peak of the Faraday depth spectra in our data set is not measurably resolved for about 8% of the lines of sight. An unresolved peak indicates a separation between the regions with Faraday rotation and synchrotron emission. However, cosmic rays pervade the ISM, and synchrotron emission would therefore also be produced where there is Faraday rotation. We suggest that the orientation of the magnetic field can separate the two effects. By modelling the thermal electron contribution to the Faraday depth, we map the strength of the magnetic field component along the line of sight. Polarized point sources in our data set have rotation measures that are comparable to the Faraday depths of the diffuse emission in our data. Our Faraday depth maps show narrow canals of low polarized intensity. We conclude that depolarization over the telescope beam produces at least some of these canals. Finally, we investigate the properties of one conspicuous region in this data set and argue that it is created by a decrease in line-of-sight depolarization compared to its surroundings.

preprint2007arXiv

WSRT Faraday tomography of the Galactic ISM at λ\sim 0.86 m

We investigate the distribution and properties of Faraday rotating and synchrotron emitting regions in the Galactic ISM in the direction of the Galactic anti-centre. We apply Faraday tomography to a radio polarization dataset that we obtained with the WSRT. We developed a new method to calculate a linear fit to periodic data, which we use to determine rotation measures from our polarization angle data. From simulations of a Faraday screen + noise we could determine how compatible the data are with Faraday screens. An unexpectedly large fraction of 14% of the lines-of-sight in our dataset show an unresolved main component in the Faraday depth spectrum. For lines-of-sight with a single unresolved component we demonstrate that a Faraday screen in front of a synchrotron emitting region that contains a turbulent magnetic field component can explain the data.

preprint2006arXiv

Detecting UHE Cosmics & Neutrinos off the Moon; an Optimal Radio Window

We show that at wavelengths comparable to the length of the shower produced by an Ultra-High Energy cosmic ray or neutrino, radio signals are an extremely efficient way to detect these particles. Through an example it is shown that this new approach offers, for the first time, the realistic possibility of measuring UHE neutrino fluxes below the Waxman-Bahcall limit. It is shown that in only one month of observing with the upcoming LOFAR radio telescope, cosmic-ray events can be measured beyond the GZK-limit, at a sensitivity level of two orders of magnitude below the extrapolated values.

preprint2006arXiv

Optimal Radio Window for the Detection of Ultra-High-Energy Cosmic Rays and Neutrinos off the Moon

When high-energy cosmic rays impinge on a dense dielectric medium, radio waves are produced through the Askaryan effect. We show that at wavelengths comparable to the length of the shower produced by an Ultra-High Energy cosmic ray or neutrino, radio signals are an extremely efficient way to detect these particles. Through an example it is shown that this new approach offers, for the first time, the realistic possibility of measuring UHE neutrino fluxes below the Waxman-Bahcall limit. It is shown that in only one month of observing with the upcoming LOFAR radio telescope, cosmic-ray events can be measured beyond the GZK-limit, at a sensitivity level of two orders of magnitude below the extrapolated values.

preprint2003arXiv

Characteristics of the structure in the Galactic polarized radio background at 350 MHz

Angular power spectra and structure functions of the Stokes parameters Q and U and polarized intensity P are derived from three sets of radio polarimetric observations. Two of the observed fields have been studied at multiple frequencies, allowing determination of power spectra and structure functions of rotation measure RM as well. The third field extends over a large part of the northern sky, so that the variation of the power spectra over Galactic latitude and longitude can be studied. The power spectra of Q and U are steeper than those of P, probably because a foreground Faraday screen creates extra structure in Q and U, but not in P. The extra structure in Q and U occurs on large scales, and therefore causes a steeper spectrum. The derived slope of the power spectrum of P is the multipole spectral index, and is consistent with earlier estimates. The multipole spectral index decreases with Galactic latitude (i.e. the spectrum becomes flatter), but is consistent with a constant value over Galactic longitude. Power spectra of the rotation measure RM show a spectral index of about 1, while the structure function of RM is approximately flat. The structure function is flatter than earlier estimates from polarized extragalactic sources, which could be due to the fact that extragalactic source RM probes the complete line of sight through the Galaxy, whereas as a result of depolarization diffuse radio polarization only probes the nearby ISM.

preprint2003arXiv

Galactic magnetic fields, from radio polarimetry of the WIM

Multi-frequency radio polarimetry of the diffuse Galactic synchrotron background gives new viewpoints on the Galactic magnetic field. Rotation measure maps reveal magnetic structures on arcminute to degree scales, such as a ring in polarization that we interpret as a magnetic tunnel. A complication using this technique is depolarization across the beam and along the line of sight. The influence of beam depolarization has been estimated using numerical models of the magneto-ionic ISM, through which polarized radiation propagates. The models show that depolarization canals similar to those observed can be caused by beam depolarization, and that the one-dimensional gradients in RM needed to produce these canals are ubiquitous in the medium.

preprint2002arXiv

Interstellar scintillation as the origin of rapid radio variability in the quasar J1819+3845

Quasars shine brightly due to the liberation of gravitational energy as matter falls onto a supermassive black hole in the centre of a galaxy. Variations in the radiation received from active galactic nuclei (AGN) are studied at all wavelengths, revealing the tiny dimensions of the region and the processes of fuelling the black hole. Some AGN are variable at optical and shorter wavelengths, and display radio outbursts over years and decades. These AGN often also show faster variations at radio wavelengths (intraday variability, IDV) which have been the subject of much debate. The simplest explanation, supported by a correlation in some sources between the optical (intrinsic) and faster radio variations, is that the rapid radio variations are intrinsic. However, this explanation implies physically difficult brightness temperatures, suggesting that the variations may be due to scattering of the incident radiation in the interstellar medium of our Galaxy. Here we present results which show unambiguously that the variations in one extreme case are due to interstellar scintillation. We also measure the transverse velocity of the scattering material, revealing a surprising high velocity plasma close to the Solar System.

preprint2001arXiv

Parsec-scale structure in the warm ISM from polarized galactic radio background observations

We present multi-frequency polarization observations of the diffuse radio synchrotron background modulated by Faraday rotation, in two directions of positive latitude. No extended total intensity I is observed, which implies that total intensity has no structure on scales smaller than approximately a degree. Polarized intensity and polarization angle, however, show abundant small-scale structure on scales from arcminutes to degrees. Rotation Measure (RM) maps show coherent structure over many synthesized beams, but also abrupt large changes over one beam. RM's from polarized extragalactic point sources are correlated over the field in each of the two fields, indicating a galactic component to the RM, but show no correlation with the RM map of the diffuse radiation. The upper limit in structure in I puts constraints on the random and regular components of the magnetic field in the galactic interstellar medium and halo. The emission is partly depolarized so that the observed polarization mostly originates from a nearby part of the medium. This explains the lack of correlation between RM from diffuse emission and from extragalactic point sources as the latter is built up over the entire path length through the medium.

preprint1999arXiv

HI and OH absorption at z=0.89

We report on WSRT observations of HI and OH absorption at z=0.885 towards the radio lens PKS 1830-21, mm wave transitions of several molecular species have already been observed at this redshift. At mm wavelengths the source structure is dominated by two extremely compact components, the northeast (NE) and southwest (SW) components. At lower frequencies the continuum emission is much more extended and there is also a broad Einstein ring connecting the NE and SW components. This means that the HI and OH spectra sample a much larger region of the absorber than the mm wave spectra. The HI spectrum that we obtain is asymmetric, with a peak at -147 km/s with respect to the main molecular line redshift of z=0.88582. Weak mm wave molecular absorption has also been detected towards the NE component at this same velocity. The HI absorption, however, covers a total velocity width of 300 km/sec, i.e. including velocities well to the red of molecular features suggesting that it is spatially widespread. In OH we detect both the 1667 and the 1665 MHz transitions. The OH spectrum has a velocity width comparable to that of the HI spectrum, suggesting that it too is widespread in the absorber. The lack of a prominent HI peak in the spectrum at the velocity corresponding to the SW component, suggests that the galaxy responsible for the absorption at z=0.885 has a central molecular disk many kpc in size, and that HI is deficient in this central region. Our observations are sensitive to the large scale kinematics of the absorber, and to first order the implied dynamical mass is consistent with the lens models of Nair et. al. (1993).

preprint1999arXiv

Measuring Cosmological Parameters with the JVAS and CLASS Gravitational Lens Surveys

The JVAS (Jodrell Bank-VLA Astrometric Survey) and CLASS (Cosmic Lens All-Sky Survey) are well-defined surveys containing about ten thousand flat-spectrum radio sources. For many reasons, flat-spectrum radio sources are particularly well-suited as a population from which one can obtain unbiased samples of gravitational lenses. These are by far the largest gravitational (macro)lens surveys, and particular attention was paid to constructing a cleanly-defined sample for the survey itself and for the underlying luminosity function. Here we present the constraints on cosmological parameters, particularly the cosmological constant, derived from JVAS and combine them with constraints from optical gravitational lens surveys, `direct' measurements of $Ω_{0}$, $H_{0}$ and the age of the universe, and constraints derived from CMB anisotropies, before putting this final result into the context of the latest results from other, independent cosmological tests.

preprint1998arXiv

The 1.4 GHz light curve of GRB 970508

We report on Westerbork 1.4 GHz radio observations of the radio counterpart to $γ$-ray burst GRB~970508, between 0.80 and 138 days after this event. The 1.4 GHz light curve shows a transition from optically thick to thin emission between 39 and 54 days after the event. We derive the slope $p$ of the spectrum of injected electrons ($dN/dγ_{e}\proptoγ_{e}^{-p}$) in two independent ways which yield values very close to $p=2.2$. This is in agreement with a relativistic dynamically near-adiabatic blast wave model whose emission is dominated by synchrotron radiation and in which a significant fraction of the electrons cool fast.

preprint1996arXiv

WN1626+5153: A Giant Radio Galaxy from the WENSS survey

An FRII radio source with an angular extension of 18.4 arcminutes has been discovered in the Westerbork Northern Sky Survey that is being carried out at 325 MHz. This source is identified with Mrk 1498 at $z=0.056$ indicating that its projected extension is 1.6 \h50\ Mpc. The optical spectrum of Mrk 1498 has the typical emission lines of a narrow line radio galaxy; H$α$ is the only permitted line that has a faint broad component. The discovery shows that the good sensitivity of WENSS at low frequencies to low surface brightness features at scales of a few minutes, is uniquely suited to provide well defined and relatively large samples of giant radio galaxies. The future definition of such samples will be important for our understanding of giant radio galaxies.

preprint1995arXiv

A radio Search for high redshift HI absorption

Ground based optical observations have yielded considerable information on the statistics of damped-lyman alpha systems. In particular these systems are known to be the dominant repository of the observed neutral gas at high redshift. However, particularly at high redshift, there is the possibility that optical observations could be biased due to the exclusion of damped-lyman alpha systems that contain moderate to significant amounts of dust. Independent observational constraints on the neutral hydrogen content at high redshifts and the amount of dust in high redshift systems can be obtained from a radio search against the bright lobes of distant radio galaxies (which is less affected by the presence of dust in foreground damped-lyman alpha systems). We describe here a pilot radio survey along the line of sight to a small sample of high redshift radio galaxies, and also present some preliminary results. The survey uses a novel observing mode at the WSRT which enables one to make sensitive searches of a large redshift interval in a modest amount of telescope time.