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.

The radio source J1819+3845 underwent a period of extreme interstellar scintillation between circa 1999 and 2007. The plasma structure responsible for this scintillation was determined to be just $1$-$3\,$pc from the solar system and to posses a density of $n_e\sim 10^2\,$cm$^{-3}$ that is three orders of magnitude higher than the ambient interstellar density (de Bruyn & Macquart 2015). Here we present radio-polarimetric images of the field towards J1819+3845 at wavelengths of 0.2, 0.92 and 2$\,$m. We detect an elliptical plasma globule of approximate size $1^\circ \times \gtrsim 2^\circ$ (major-axis position angle of $\approx -40^\circ$), via its Faraday-rotation imprint ($\approx 15\,$rad$\,$m$^{-2}$) on the diffuse Galactic synchrotron emission. The extreme scintillation of J1819+3845 was most likely caused at the turbulent boundary of the globule (J1819+3845 is currently occulted by the globule). The origin and precise nature of the globule remain unknown. Our observations are the first time plasma structures that likely cause extreme scintillation have been directly imaged.

Leakage of polarized Galactic diffuse emission into total intensity can potentially mimic the 21-cm signal coming from the epoch of reionization (EoR), as both of them might have fluctuating spectral structure. Although we are sensitive to the EoR signal only in small fields of view, chromatic sidelobes from further away can contaminate the inner region. Here, we explore the effects of leakage into the 'EoR window' of the cylindrically averaged power spectra (PS) within wide fields of view using both observation and simulation of the 3C196 and NCP fields, two observing fields of the LOFAR-EoR project. We present the polarization PS of two one-night observations of the two fields and find that the NCP field has higher fluctuations along frequency, and consequently exhibits more power at high-$k_\parallel$ that could potentially leak to Stokes $I$. Subsequently, we simulate LOFAR observations of Galactic diffuse polarized emission based on a model to assess what fraction of polarized power leaks into Stokes $I$ because of the primary beam. We find that the rms fractional leakage over the instrumental $k$-space is $0.35\%$ in the 3C196 field and $0.27\%$ in the NCP field, and it does not change significantly within the diameters of $15^\circ$, $9^\circ$ and $4^\circ$. Based on the observed PS and simulated fractional leakage, we show that a similar level of leakage into Stokes $I$ is expected in the 3C196 and NCP fields, and the leakage can be considered to be a bias in the PS.

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 $\Delta^2_{\rm I} = (56 \pm 13 {\rm mK})^2$ (1-$\sigma$) excess variance and a best 2-$\sigma$ upper limit of $\Delta^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.

The new generation of low-frequency radio telescopes, such as the Low Frequency Array (LOFAR: a Square Kilometre Array-low pathfinder), provides advancements in our capability of probing Galactic magnetism through low-frequency polarimetry. Maps of diffuse polarized radio emission and Faraday rotation can be used to infer properties of, and trace structure in, the magnetic fields in the ISM. However, to date very little of the sky has been probed at high angular and Faraday depth resolution. We observed a 5x5 degree region centred on the nearby galaxy IC342 using LOFAR in the frequency range 115-178 MHz at 4 arcmin resolution and performed Faraday tomography to detect foreground Galactic polarized synchrotron emission separated by Faraday depth (different amounts of Faraday rotation). Our Faraday depth cube shows rich polarized structure, with up to 30 K of polarized emission at 150 MHz. We detect two overlapping diffuse polarized features that are clearly separated in Faraday depth. Faraday-thick structures at such low frequencies would be too strongly depolarized to explain the observations and are therefore rejected. Only Faraday thin structures will not be strongly depolarized; producing such structures requires localized variations in the ratio of synchrotron emissivity to Faraday depth per unit distance, which can arise from several physical phenomena, such as a transition between regions of ionized and neutral gas. We conclude that the observed polarized emission is Faraday thin, and propose that the emission originates from two neutral clouds in the local ISM. We have modeled the Faraday rotation for this line of sight and estimated that the line of sight component of magnetic field of the local ISM for this direction varies between -0.86 and +0.12 uG. We propose that this may be a useful method for mapping magnetic fields within the local ISM.

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.

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.

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.

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.

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 \mu G}$ and its spatial variation across the 3C196 field is $0.1~{\rm \mu 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}\mu G}$) compared to its surroundings.

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.

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.

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.

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.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

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

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.

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} {\Delta L}_{100}^{-1/2}$cm$^{-3}$, where $l_0$ is the outer scale of the turbulence in AU and $\Delta L = 100 \Delta 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.

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<\nu<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$).

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

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 ($\Delta z$) of reionization can be constrained assuming a parametrization. We use an EoR simulation of $z_r = 7.68$ and $\Delta 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 $\Delta 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.

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.

We report results from a programme aimed at investigating the temperature of neutral gas in high-redshift damped Lyman-$\alpha$ 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\sigma$) 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\sigma$) 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\sigma$ 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.)

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.

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.

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.

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 \nu^{-2} dependence, (ii) ionospheric refraction leads to a knee-like modulation on the sky spectrum at \nu\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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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]

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.

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.

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.

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.

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.

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.

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 $\beta^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 $\beta^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 $\delta T= \sqrt{\ell (\ell+1) C^I_\ell / 2\pi} \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 $\delta 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.

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.

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.

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.

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. First results are presented of an analysis of 20 hours of observation data for NuMoon project using the Westerbork Synthesis Radio Telescope to search for short radio pulses from the Moon. A limit on the neutrino flux is set that is a factor four better than the current one (based on FORTE).

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.

Low-frequency radio observations of galaxy clusters are the key to detecting the diffuse extended emission associated with them. The presence and properties of such radio sources in galaxy clusters reveal the existence of magnetic fields on a large scale and allow theories to be tested concerning both the origin of relativistic particles in the ICM and their propagation. A deep radio bservation of the A2255 galaxy cluster was carried out at 85 cm with the WSRT. The good UV-coverage and sensitivity achieved by these observations allowed us to image the complex structure of the low-brightness, extended cluster sources (radio halo and relic). These sources show a larger extent than what has been imaged so far at this frequency, with two new structures located SW and NW of the cluster center and at projected distances of 2 Mpc from it. The physical properties of the newly detected structures, together with the active dynamical state of the cluster, support a connection with large-scale structure (LSS) formation shocks.

We present Very Large Array images of a "Double-Double Radio Galaxy", a class of objects in which two pairs of lobes are aligned either side of the nucleus. In this object, B0925+420, we discover a third pair of lobes, close to the core and again in alignment with the other lobes. This first-known "Triple-Double" object strongly increases the likelihood that these lobes represent mutiple episodes of jet activity, as opposed to knots in an underlying jet. We model the lobes in terms of their dynamical evolution. We find that the inner pair of lobes is consistent with the outer pair having been displaced buoyantly by the ambient medium. The middle pair of lobes is more problematic - to the extent where an alternative model interpreting the middle and inner "lobes" as additional bow shocks within the outer lobes may be more appropriate - and we discuss the implications of this on our understanding of the density of the ambient medium.

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.

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.

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.

In this paper we use the diffuse Galactic synchrotron emission to study the properties of the Galactic magnetic field. We combined data from the low-frequency WENSS survey with single-dish observations carried out by Brouw and Spoelstra to (partially) fill in the central gap in the (u,v) plane that is missing from the interferometer data. The small bandwidth of the WENSS data meant that we could not determine rotation measures (RM) directly. Instead we interpreted the polarization angle gradients that we derived as gradients in RM, and to do this we found a new way to efficiently and reliably fit linear gradients to periodic data. RM are available for the single-dish observations, and we found that we have to rescale these single-dish RM to match the RM gradients that we derive from the combined WENSS/single-dish dataset with the single-dish only RM gradients. We tentatively show that the difference in beamsize between the datasets can be responsible for this. We interpret the scaled-up RM we find in terms of a simple toy model, and by combining our results with those by Haverkorn et al. (2004) we reconstruct the full 3D magnetic field vector along a number of lines-of-sight. For these lines-of-sight we derive the properties of the magnetic field component perpendicular to the Galactic plane.

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.

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.

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.

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.

We report the discovery of a new gravitational lens system from the CLASS survey, CLASS B0631+519. VLA, MERLIN and VLBA observations show a doubly-imaged radio core, a doubly-imaged lobe and a second lobe that is probably quadruply-imaged. The maximum image separation is 1.16 arcseconds. The VLBA resolves the most magnified image of the flat-spectrum radio core into a number of sub-components spread across approximately 20 milli-arcseconds. Optical and near-infrared imaging with the ACS and NICMOS cameras on the HST show that there are two galaxies along the line of sight to the lensed source, as was previously discovered by optical spectroscopy. The nearer galaxy at z=0.0896 is a small blue irregular, while the more distant galaxy at z=0.6196 is an elliptical type and appears to contribute most of the lensing effect. The host galaxy of the lensed source is visible in NICMOS imaging as a set of arcs that form an almost complete Einstein ring. Mass modelling using non-parametric techniques can reproduce the ring and indicates that the irregular galaxy has a (localised) effect on the flux density distribution in the Einstein ring at the 5-10% level.

We present observations of CLASS B2108+213, the widest separation gravitational lens system discovered by the Cosmic Lens All-Sky Survey. Radio imaging using the VLA at 8.46 GHz and MERLIN at 5 GHz shows two compact components separated by 4.56 arcsec with a faint third component in between which we believe is emission from a lensing galaxy. 5-GHz VLBA observations reveal milliarcsecond-scale structure in the two lensed images that is consistent with gravitational lensing. Optical emission from the two lensed images and two lensing galaxies within the Einstein radius is detected in Hubble Space Telescope imaging. Furthermore, an optical gravitational arc, associated with the strongest lensed component, has been detected. Surrounding the system are a number of faint galaxies which may help explain the wide image separation. A plausible mass distribution model for CLASS B2108+213 is also presented.

Radio waves are imprinted with propagation effects from ionized media through which they pass. Owing to electron density fluctuations, compact sources (pulsars, masers, and compact extragalactic sources) can display a wide variety of scattering effects. These scattering effects, particularly interstellar scintillation, can be exploited to provide *superresolution*, with achievable angular resolutions (<~ 1 microarcsecond) far in excess of what can be obtained by very long baseline interferometry on terrestrial baselines. Scattering effects also provide a powerful sub-AU probe of the microphysics of the interstellar medium, potentially to spatial scales smaller than 100 km, as well as a tracer of the Galactic distribution of energy input into the interstellar medium through a variety of integrated measures. Coupled with future gamma-ray observations, SKA observations also may provide a means of detecting fainter compact gamma-ray sources. Though it is not yet clear that propagation effects due to the intergalactic medium are significant, the SKA will either detect or place stringent constraints on intergalactic scattering.

The polarized component of the diffuse radio synchrotron emission of our Galaxy shows structure, which is apparently unrelated to the structure in total intensity, on many scales. The structure in the polarized emission can be due to several processes or mechanisms. Some of those are related to the observational setup, such as beam depolarization -- the vector combination and (partial) cancellation of polarization vectors within a synthesized beam --, or the insensitivity of a synthesis telescope to structure on large scales, also known as the 'missing short spacings problem'. Other causes for structure in the polarization maps are intrinsic to the radiative transfer of the emission in the warm ISM, which induces Faraday rotation and depolarization. We use data obtained with the Westerbork Synthesis Radio Telescope at 5 frequencies near 350 MHz to estimate the importance of the various mechanisms in producing structure in the linearly polarized emission. In the two regions studied here, which are both at positive latitudes in the second Galactic quadrant, the effect of 'missing short spacings' is not important. The properties of the narrow depolarization 'canals' that are observed in abundance lead us to conclude that they are mostly due to beam depolarization, and that they separate regions with different rotation measures. As beam depolarization only creates structure on the scale of the synthesized beam, most of the structure on larger scales must be due to depth depolarization. We do not discuss that aspect of the observations here, but in a companion paper we derive information about the properties of the ISM from the structure of the polarized emission.

We describe a first attempt to derive properties of the regular and turbulent Galactic magnetic field from multi-frequency polarimetric observations of the diffuse Galactic synchrotron background. A single-cell-size model of the thin Galactic disk is constructed which includes random and regular magnetic fields and thermal and relativistic electrons. The disk is irradiated from behind with a uniform partially polarized background. Radiation from the background and from the thin disk is Faraday rotated and depolarized while propagating through the medium. The model parameters are estimated from a comparison with 350 MHz observations in two regions at intermediate latitudes done with the Westerbork Synthesis Radio Telescope. We obtain good consistency between the estimates for the random and regular magnetic field strengths and typical scales of structure in the two regions. The regular magnetic field strength found is a few microGauss, and the ratio of random to regular magnetic field strength is 0.7 +/- 0.5, for a typical scale of the random component of 15 +/- 10 pc. Furthermore, the regular magnetic field is directed almost perpendicular to the line of sight. This modeling is a potentially powerful method to estimate the structure of the Galactic magnetic field, especially when more polarimetric observations of the diffuse synchrotron background at intermediate latitudes become available.

We present a kinematic model of the neutral-hydrogen in the spiral galaxy of the lens system PKS1830-211, based on a Multi Element Radio-Linked Interferometer (MERLIN) 1.4-GHz radio map and the integrated and redshifted 21-cm hydrogen absorption-line profile as measured with the Westerbork Synthesis Radio Telescope (WSRT). Degeneracies in the models do not allow a unique determination of the kinematic center and forthcoming deeper Hubble Space Telescope observations with the Advanced Camera for Surveys (ACS) are required to break this degeneracy. Even so, we measure the inclination of the hydrogen disk: i=(17-32) deg., indicating a close-to face-on spiral galaxy. The optical depth increases with radius over the extent of the Einstein ring, suggesting HI depletion towards the lens center. The latter could be due to star formation or conversion of HI in to molecular hydrogen because of a higher metalicity/dust content in the galaxy center. The neutral hydrogen optical depth gives N_HI= 2x10^21 cm^-2 at r=5.0 h_{70}^-1 kpc in the disk (T_s=100K), comparable to local spiral galaxies. Our study shows that planned new radio telescopes (i.e. ALMA, LOFAR and SKA) -- which will discover large numbers of similar lens systems -- are powerful new tools to probe the internal structure, kinematics and evolution of spiral galaxies in detail to z>=1, as well as their neutral-hydrogen (and molecular) content, thereby complementing studies of HI emission from spirals at z<=1 and damped Ly-alpha systems to z>>1.

We present deep GMRT OH and WSRT HI absorption spectra of the $z = 0.6846$ gravitational lens toward B0218+357. Both the 1665 MHz and 1667 MHz OH lines are clearly detected for the first time while a new wide absorption component was detected (at low significance) in the HI spectrum; the OH spectra yielded an OH column density of $N_{OH} = 2.3 \times 10^{15}$ cm^{-2}. The ratio of 1667 and 1665 MHz equivalent widths is $\sim 1.8$ while the redshift of peak OH absorption is $z = 0.68468 \pm 0.000008$ for both lines; this redshift agrees with that obtained from the HI line. The velocity spread (between nulls) of the HI absorption is $\sim 140$ km/s, while that of both OH lines is $\sim 100$ km/s; the HI and OH spectra are broadly similar in that they each have two principal narrow components and a wide absorption trough. We argue that the wide absorption is likely to arise from source components in the Einstein ring and derive a rotation velocity of $\sim 150$ km/s at a distance of 1.5 kpc from the centre of the $z \sim 0.6846$ galaxy.

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.

LOFAR, the Low Frequency Array, is a large radio telescope consisting about 100 soccer field sized antenna stations spread over a region of 400 km in diameter. It will operate in the frequency range from ~10 to 240 MHz, with a resolution at 240 MHz of better than an arcsecond. Its superb sensitivity will allow for a broad range of astrophysical studies. In this contribution we first discuss four major areas of astrophysical research in which LOFAR will undoubtedly make important contributions: reionisation, distant galaxies and AGNs, transient radio sources and cosmic rays. Subsequently, we will discuss the technical concept of the instrument and the status of the LOFAR project

We present flux-ratio curves of the fold and cusp (i.e. close multiple) images of six JVAS/CLASS gravitational lens systems. The data were obtained over a period of 8.5 months in 2001 with the Multi-Element Radio-Linked Interferometer Network (MERLIN) at 5-GHz with 50 mas resolution, as part of a MERLIN Key-Project. Even though the time delays between the fold and cusp images are small (<~1 day) compared to the time-scale of intrinsic source variability, all six lens systems show evidence that suggests the presence of extrinsic variability. In particular, the cusp images of B2045+265 -- regarded as the strongest case of the violation of the cusp relation (i.e. the sum of the magnifications of the three cusp images add to zero) -- show extrinsic variations in their flux-ratios up to ~40 percent peak-to-peak on time scales of several months. Its low Galactic latitude of b=-10 degree and a line-of-sight toward the Cygnus superbubble region suggest that Galactic scintillation is the most likely cause. The cusp images of B1422+231 at b=+69 degree do not show strong extrinsic variability. Galactic scintillation can therefore cause significant scatter in the cusp and fold relations of some radio lens systems (up to 10 percent rms), even though these relations remain violated when averaged over a <~1 year time baseline.

We study a conspicuous ring-like structure with a radius of about 1.4 degrees which was observed with the Westerbork Synthesis Radio Telescope (WSRT) at 5 frequencies around 350 MHz. This ring is very prominent in Stokes Q and U, and less so in polarized intensity P. No corresponding structure is visible in total intensity Stokes I, which indicates that the ring is created by Faraday rotation and depolarization processes. The polarization angle changes regularly from the center of the ring outwards to a radius of more than 1.7 degrees. Thus, the structure in polarization angle is not ring-like but resembles a disk, and it is larger than the ring in P. The rotation measure (RM) decreases almost continuously over the disk, from RM = 0 rad/m2 at the edge, to -8 rad/m2 in the center, while outside the ring the RM is slightly positive. This radial variation of RM yields stringent constraints on the nature of the ring-like structure, because it rules out any spherically symmetrical magnetic field configuration, such as might be expected from supernova remnants or wind-blown bubbles. We discuss several possible connections between the ring and known objects in the ISM, and conclude that the ring is a predominantly magnetic funnel-like structure. This description can explain both the field reversal from outside to inside the ring, and the increase in magnitude of RM towards the center of the ring. The ring-structure in P is most likely caused by a lack of depolarization due to a very uniform RM distribution at that radius. Beyond the ring, the RM gradient increases, depolarizing the polarized emission, so that P decreases. In the SW corner of the field a pattern of narrow filaments of low polarization, aligned with Galactic longitude, is observed, indicative of beam depolarization due to abrupt changes in RM.

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.

We present two years of monitoring observations of the extremely variable quasar J1819+3845. We observe large yearly changes in the timescale of the variations (from ~ 1 hour to ~ 10 hours at 5GHz). This annual effect can only be explained if the variations are caused by a propagation effect, and thus affected by the Earth's relative speed through the projected intensity pattern. To account for this effect, the scattering plasma must have a transverse velocity with respect to the local standard of rest. The velocity calculated from these observations is in good agreement with that obtained from a two telescope delay experiment (Dennett-Thorpe & de Bruyn 2001). We also show that either the source itself is elongated, or that the scattering plasma is anisotropic, with an axial ratio of >6:1. As the source is extended on scales relevant to the scattering phenomenon, it seems plausible that the anisotropy is due to the source itself, but this remains to be investigated. From the scintillation characteristics we find that the scattering material is a very strong, thin scatterer within ~ten parsecs. We determine a source size at 5GHz of 100 to 900microarcsecs, and associated brightness temperatures of 10^{10} to 10^{12}K.

We report the final results of the search for gravitationally lensed flat-spectrum radio sources found in the combination of CLASS (Cosmic Lens All Sky Survey) and JVAS (Jodrell-Bank VLA Astrometric Survey). VLA observations of 16,503 sources have been made, resulting in the largest sample of arcsec-scale lens systems available. Contained within the 16,503 sources is a complete sample of 11,685 sources having two-point spectral indices between 1.4 and 5 GHz flatter than -0.5 and 5 GHz flux densities $\geq$30 mJy. A subset of 8,958 sources form a well-defined statistical sample suitable for analysis of the lens statistics. We describe the systematic process by which 149 candidate lensed sources were picked from the statistical sample on the basis of possessing multiple compact components in the 0.2 arcsec-resolution VLA maps. Candidates were followed up with 0.05 arcsec resolution MERLIN and 0.003 arcsec VLBA observations at 5 GHz and rejected as lens systems if they failed well-defined surface brightness and/or morphological tests. Maps for all the candidates can be found on the World Wide Web at http://www.jb.man.ac.uk/research/gravlens/index.html We summarize the properties of each of the 22 gravitational lens systems in JVAS/CLASS. Twelve are double-image systems, nine are four-image systems and one is a six-image system. Thirteen constitute a statistically well-defined sample giving a point-source lensing rate of 1:690$\pm$190. The interpretation of the results in terms of the properties of the lensing galaxy population and cosmological parameters will be published elsewhere. (Abridged)

The Cosmic Lens All-Sky Survey (CLASS) is an international collaborative program which has obtained high-resolution radio images of over 10000 flat-spectrum radio sources in order to create the largest and best studied statistical sample of radio-loud gravitationally lensed systems. With this survey, combined with detailed studies of the lenses found therein, constraints can be placed on the expansion rate, matter density, and dark energy (e.g. cosmological constant, quintessence) content of the Universe that are complementary to and independent of those obtained through other methods. CLASS is aimed at identifying lenses where multiple images are formed from compact flat-spectrum radio sources, which should be easily identifiable in the radio maps. Because CLASS is radio-based, dust obscuration in lensing galaxies is not a factor, and the relative insensitivity of the instrument to environmental conditions leads to nearly uniform sensitivity and resolution over the entire survey. In four observing seasons from 1994-1999, CLASS has observed 13783 radio sources with the VLA at 8.4 GHz at 0.2 arcsecond resolution. When combined with the JVAS survey, the CLASS sample contains over 16,000 images. A complete sample of 11685 flat-spectrum sources was observed, selected from GB6 catalogue at 4.85 GHz and the NVSS at 1.4 GHz. So far, CLASS has found 16 new gravitational lens systems, and the JVAS/CLASS survey contains a total of 22 lenses. (Abridged)

We report the discovery of a new gravitational lens system from the CLASS survey. Radio observations with the VLA, the WSRT and MERLIN show that the radio source B0850+054 is comprised of two compact components with identical spectra, a separation of 0.7 arcsec and a flux density ratio of 6:1. VLBA observations at 5 GHz reveal structures that are consistent with the gravitational lens hypothesis. The brighter of the two images is resolved into a linear string of at least six sub-components whilst the weaker image is radially stretched towards the lens galaxy. UKIRT K-band imaging detects an 18.7 mag extended object, but the resolution of the observations is not sufficient to resolve the lensed images and the lens galaxy. Mass modelling has not been possible with the present data and the acquisition of high-resolution optical data is a priority for this system.

A new two-image gravitational lens system has been discovered as a result of the Cosmic Lens All-Sky Survey (CLASS). Radio observations with the VLA, MERLIN and the VLBA at increasingly higher resolutions all show two components with a flux density ratio of ~7:1 and a separation of 1.34". Both components are compact and have the same spectral index. Followup observations made with the VLA at 8.4 GHz show evidence of a feature to the south-east of the brighter component and a corresponding extension of the weaker component to the north-west. Optical observations with the WHT show ~1.7" extended emission aligned in approximately the same direction as the separation between the radio components with an R-band magnitude of 21.8 +/- 0.4.

We derive constraints on cosmological parameters and the properties of the lensing galaxies from gravitational lens statistics based on the final Cosmic Lens All Sky Survey (CLASS) data. For a flat universe with a classical cosmological constant, we find that the present matter fraction of the critical density is $\Omega_{\rm m}=0.31^{+0.27}_{-0.14}$ (68%) $^{+0.12}_{-0.10}$ (systematic). For a flat universe with a constant equation of state for dark energy $w = p_x({pressure})/\rho_x({energy density})$, we find $w < -0.55^{+0.18}_{-0.11}$ (68%).

The 352 MHz Westerbork In the Southern Hemisphere (WISH) survey is the southern extension of the WENSS, covering 1.60 sr between -9 < DEC < -26 to a limiting flux density of ~18 mJy (5sigma). Due to the very low elevation of the observations, the survey has a much lower resolution in declination than in right ascension (54" x 54"cosec(DEC)). A correlation with the 1.4 GHz NVSS shows that the positional accuracy is less constrained in declination than in right ascension, but there is no significant systematic error. We present a source list containing 73570 sources. We correlate this WISH catalogue with the NVSS to construct a sample of faint Ultra Steep Spectrum (USS) sources, which is accessible for follow-up studies with large optical telescopes in the southern hemisphere. This sample is aimed at increasing the number of known high redshift radio galaxies to allow detailed follow-up studies of these massive galaxies and their environments in the early Universe.

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.

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.

We present multi-frequency polarization observations of the diffuse radio synchrotron background modulated by Faraday rotation. No total intensity is observed, indicating that total intensity does not vary on scales below approximately a degree. However, polarized intensity and polarization angle show abundant small-scale structure due to Faraday rotation in the warm ionized disk on small scales. The distribution of Rotation Measures enables us to estimate structure in magnetic field, weighted with electron density, in the warm ionized disk.

We have used the Westerbork Northern Sky Survey (WENSS) to define a sample of 47 low redshift (z < 0.4) giant radio galaxies. This is the largest sample yet of such radio sources originating from a single survey. We present radio maps of the newly discovered giants and optical images and spectra of their host galaxies. We present some of their properties and discuss the selection effects. We investigate the distribution of the sources in the radio power -- linear size (P-D) diagram, and how these parameters relate to the redshifts of the sources in our sample. We find a strong drop in the number of sources with a linear size above 2 Mpc. We argue that this is not a result of selection effects, but that it indicates either a strong luminosity evolution of radio sources of such a size, or that a large fraction of radio sources `switch off' before they are able to grow to a size of 2 Mpc.

We report the discovery of a new two-image gravitational lens system from the Cosmic Lens All-Sky Survey, CLASS B2319+051. Radio imaging with the Very Large Array (VLA) and Multi-Element Radio-Linked Interferometer Network (MERLIN) shows two compact components with a flux density ratio of 5:1, separated by 1.36 arcsec. Observations with the Very Long Baseline Array (VLBA) resolve each of the radio components into a pair of parity-reversed subcomponents. Hubble Space Telescope (HST) observations with the Near-Infrared Camera and Multi-Object Spectrometer (NICMOS) show a bright elliptical galaxy (G1) coincident with the radio position, and a second irregular galaxy (G2) 3.4 arcsec to the northwest. Previous spectroscopic studies have indicated that these galaxies are at different redshifts: z(G1) = 0.624, z(G2) = 0.588. Infrared counterparts to the lensed radio components are not detected in the NICMOS image, and the source redshift has not yet been determined. Preliminary mass modeling based on the VLBA subcomponent data indicates that the lensing potential includes a strong external shear contribution. A VLA monitoring program is currently being undertaken to measure the differential time delay.

Westerbork Radio Synthesis Telescope observations of the redshifted 21cm line absorber against the z_em = 0.871 double lobed quasar 3C196 show that the intervening absorber is an L approximately equal L* spiral galaxy (3C196-G1) and that the absorbing layer of cold gas extends to radii of at least 30 h_50^-1 kpc. The new data solve several long standing puzzles about this system by (1) discovering a second 21cm absorption feature, corresponding to absorption against the NE lobe of the background radio source and (2) spatially ``resolving'' the two absorption features to isolate the absorption along the two lines of sight to the opposing radio lobes. These findings resolve the disagreement in redshift between the UV metal and 21cm lines, and as well as demonstrating that the neutral layer does absorb both lobes of the background radio source. Simple kinematic models with an inclined, rotating gas disk match the observed 21cm profile and are also compatible with both the redshift and velocity spread of the absorption measured in UV resonance lines along a third, independent line of sight to the quasar nucleus and with the lack of 21cm absorption in as earlier VLBI experiment that was sensitive to opacity against the hot spot in the northern lobe. The inferred rotation speed and luminosity for the galaxy are compatible with the z approximately 0 Tully-Fisher Relation. This system illustrates well how 21cm absorption against extended background radio sources is a powerful tool in determining the nature of the damped Lyman-alpha class of QSO absorption line system.

We present measurements of absorption by the 21cm hyperfine transition of neutral hydrogen toward radio sources at substantial look-back times. These data are used in combination with observations of rotational transitions of common interstellar molecules, to set limits on the evolution of the fine structure constant: (dot alpha)/alpha < 3.5e-15/yr to a look-back time of 4.8 Gyr. The neutral hydrogen observations employed Very Long Baseline Interferometry in order to mitigate the substantial uncertainty arising from the fact that observations at very different wavelengths may probe different lines-of-site due to frequency dependent structure of the background source. We discuss the implication of these results on theories unifying natural forces based on compact quantum dimensions. In the context of string theory, the limit on the secular evolution of the scale factor of these compact dimensions, R, is (dot R)/R < 1e-15/yr. Including terrestrial and other astronomical measurements places limits (2sigma) on slow oscillations of R from the present to the epoch of cosmic nucleosynthesis, just seconds after the big bang, of (Delta R)/R < 1e-5.

We present new WSRT observations of the micro-arcsecond quasar J1819+3845. All short term variations are attributed to interstellar scintillation of a source which is at most 30 micro-arcseconds in diameter. The timescale of the modulations changes over the year, which we interpret as due to a peculiar velocity of the scattering medium. The scintillation behaviour can be used to determine sub-structure in the source.

Strongly correlated non-intrinsic variability between 5 and 8.5 GHz has been observed in one of the lensed images of the gravitational lens B1600+434. These non-intrinsic (i.e. `external') variations are interpreted as `radio-microlensing' of relativistic micro-arsec-scale jet components in the source at a redshift of z=1.59 by massive compact objects in the halo of the edge-on disk lens galaxy at z=0.41. We shortly summarize these observations and discuss several new observational and theoretical programs to investigate this new phenomenon in more detail.

High resolution MERLIN observations of a newly-discovered four-image gravitational lens system, B0128+437, are presented. The system was found after a careful re-analysis of the entire CLASS dataset. The MERLIN observations resolve four components in a characteristic quadruple-image configuration; the maximum image separation is 542 mas and the total flux density is 48 mJy at 5 GHz. A best-fit lens model with a singular isothermal ellipsoid results in large errors in the image positions. A significantly improved fit is obtained after the addition of a shear component, suggesting that the lensing system is more complex and may consist of multiple deflectors. The integrated radio spectrum of the background source indicates that it is a GigaHertz-Peaked Spectrum (GPS) source. It may therefore be possible to resolve structure within the radio images with deep VLBI observations and thus better constrain the lensing mass distribution.

We present deep WSRT 1.4 GHz observations of the Hubble Deep Field region. At the 5 sigma level, the WSRT clearly detects 85 regions of radio emission in a 10' x 10' field centred on the HDF. Eight of these regions fall within the HDF itself, four of these are sources that have not previously been detected at 1.4 GHz, although two of these are VLA detections at 8.5GHz. The two new radio sources detected by the WSRT are identified with relatively bright (I<21m) moderate redshift spiral and irregular type galaxies. In the full field, the WSRT detects 22 regions of emission that were not previously detected by the VLA at 1.4GHz. At least two of these are associated with nearby, extended star-forming galaxies.

We have formed a complete sample of 26 low redshift (z < 0.3) giant radio galaxies (GRGs) from the WENSS survey, selected at flux densities above 1 Jy at 325 MHz. We present 10.5-GHz observations with the 100-m Effelsberg telescope of 18 sources in this sample. These observations, together with similar data of the remaining eight sources, are combined with data from the WENSS, NVSS and GB6 surveys to study the radio properties of the lobes of these sources at arcminute resolution. We investigate radio source asymmetries, equipartition energy densities in the lobes, the presence of lobe pressure evolution with redshift, the spectral age and the density of the environments of these sources. We find that the armlength asymmetries of GRGs are slightly larger than those of smaller sized 3CR radio galaxies and that these are difficult to explain as arising from orientation effects only. We also find indications that the lobes of the GRGs, despite their large sizes, are still overpressured with respect to their environment. Further, we argue that any evolution of lobe pressure with redshift in these large sources is due to selection effects. A spectral ageing analysis suggests that the GRGs in our sample are the oldest members of the group of relatively high power radio sources whose radio powers have evolved to their currently observed lower values.

We present the discovery of CLASS B0739+366, a new gravitational lens system from the Cosmic Lens All-Sky Survey. Radio imaging of the source with the Very Large Array (VLA) shows two compact components separated by $0\farcs54$, with a flux density ratio of $\sim$ 6:1. High-resolution follow-up observations using the Very Long Baseline Array (VLBA) at 1.7 GHz detect weak, parity-reversed jet emission from each of the radio components. Hubble Space Telescope NICMOS F160W observations detect infrared counterparts to the lensed images, as well as an extended object between them which we identify as the lensing galaxy. Redshifts for the galaxy and lensed source have not yet been obtained. For typical lens and source redshifts of $z=0.5$ and $z=1.5$, respectively, preliminary mass modeling predicts a time delay of $\sim7h^{-1}$ days in a flat $\Omega_{M}=1.0$ universe. The small predicted time delay and weak radio components will make CLASS B0739+366 a challenging target for Hubble constant determination.

We present radio, optical, near-infrared and spectroscopic observations of the source B0827+525. We consider this source as the best candidate from the Cosmic Lens All-Sky Survey (CLASS) for a `dark lens' system or binary radio-loud quasar. The system consists of two radio components with somewhat different spectral indices, separated by 2.815 arcsec. VLBA observations show that each component has substructure on a scale of a few mas. A deep K-band exposure with the W.M.Keck-II Telescope reveals emission near both radio components. The K-band emission of the weaker radio component appears extended, whereas the emission from the brighter radio component is consistent with a point source. Hubble Space Telescope F160W-band observations with the NICMOS instrument confirms this. A redshift of 2.064 is found for the brighter component, using the LRIS instrument on the W.M.Keck-II Telescope. The probability that B0827+525 consists of two unrelated compact flat-spectrum radio sources is ~3%, although the presence of similar substructure in both component might reduce this. We discuss two scenarios to explain this system: (i) CLASS B0827+525 is a `dark lens' system or (ii) B0827+525 is a binary radio-loud quasar. B0827+525 has met all criteria that thus far have in 100% of the cases confirmed a source as an indisputable gravitational lens system. Despite this, no lens galaxy has been detected with m_F160W<=23 mag. Hence, we might have found the first binary radio-loud quasar. At this moment, however, we feel that the `dark lens' hypothesis cannot yet be fully excluded.

First, we review the current status of the detection of strong `external' variability in the CLASS gravitational B1600+434, focusing on the 1998 VLA 8.5-GHz and 1998/9 WSRT multi-frequency observations. We show that this data can best be explained in terms of radio-microlensing. We then proceed to show some preliminary results from our new multi-frequency VLA monitoring program, in particular the detection of a strong feature (~30%) in the light curve of the lensed image which passes predominantly through the dark-matter halo of the lens galaxy. We tentatively interpret this event, which lasted for several weeks, as a radio-microlensing caustic crossing, i.e. the superluminal motion of a micro-arcsec-scale jet-component in the lensed source over a single caustic in the magnification pattern, that has been created by massive compact objects along the line-of-sight to the lensed image.

We present the first unambiguous case of external variability of a radio gravitational lens, CLASS B1600+434. The VLA 8.5-GHz difference light curve of the lensed images, taking the proper time-delay into account, shows the presence of external variability with 14.6-sigma confidence. We investigate two plausible causes of this external variability: scattering by the ionized component of the Galactic interstellar medium and microlensing by massive compact objects in the bulge/disk and halo of the lens galaxy. Based on the tight relation between the modulation-index and variability time scale and the quantitative difference between the light curves of both lensed images, we conclude that the observed short-term variability characteristics of the lensed images are incompatible with scintillation in our Galaxy. This conclusion is strongly supported by multi-frequency WSRT observations at 1.4 and 5 GHz, which are in strong disagreement with predictions based on the scintillation hypothesis. ... On the other hand, a single superluminal jet-component, having an apparent velocity 9<=(v_app/c)<=26, a radius of 2-5 micro-arcsec and containing 5-11% of the observed 8.5-GHz source flux density, can reproduce the observed modulation-indices and variability time scale at 8.5 GHz, when it is microlensed by compact objects in the lens galaxy. It also reproduces the frequency-dependence of the modulation-indices, determined from the independent WSRT 1.4 and 5-GHz observations. ... The only conclusion fully consistent with the data gathered thus far is that we have indeed detected radio microlensing. The far reaching consequence of this statement is that a significant fraction of the mass in the dark-matter halo at ~6 kpc (h=0.65) above the lens-galaxy disk in B1600+434 consists of massive compact objects. [abridged]

We discuss observations of the linearly polarized component of the diffuse galactic radio background. These observations, with an angular resolution of 4', were made with the Westerbork Synthesis Radio Telescope (WSRT) in 5 frequency bands in the range 341-375 MHz. The linearly polarized intensity P (with polarized brightness temperature going up to 10K) shows a `cloudy' structure, with characteristic scales of 15'-30', which contains relatively long, but very narrow `canals' (essentially unresolved) in which P is only a small fraction of that in the neighbouring beams. These `canals' are generally seen in more than one frequency band, although their appearance changes between bands. They are probably due to depolarization within the synthesized beam, because the change in polarization angle across the deepest `canals' is in general close to 90 degrees (or 270 etc.). These very abrupt changes in polarization angle, which are seen only across the `canals', seem to be accompanied by abrupt changes in the Rotation Measure (RM), which may have the right magnitude to create the difference of close to 90 degrees in polarization angle, and thereby the `canals'. The structure in the polarization maps is most likely due to Faraday rotation modulation of the probably smooth polarized radiation emitted in the halo of our Galaxy by the fairly local ISM (up to 500 pc). Therefore, the abrupt changes of RM across the `canals' provide evidence for very thin (about 1 pc), and relatively long transition regions in the ISM, across which the RM changes by as much as 100%. Such drastic RM changes may well be due primarily to abrupt changes in the magnetic field.

We present Very Large Array (VLA) 8.5-GHz light curves of the two lens images of the Cosmic Lens All Sky Survey (CLASS) gravitational lens B1600+434. We find a nearly linear decrease of 18-19% in the flux densities of both lens images over a period of eight months (February-October) in 1998. Additionally, the brightest image A shows modulations up to 11% peak-to-peak on scales of days to weeks over a large part of the observing period. Image B varies significantly less on this time scale. We conclude that most of the short-term variability in image A is not intrinsic source variability, but is most likely caused by microlensing in the lens galaxy. The alternative, scintillation by the ionized Galactic ISM, is shown to be implausible based on its strong opposite frequency dependent behavior compared with results from multi-frequency WSRT monitoring observations (Koopmans & de Bruyn 1999). From these VLA light curves we determine a median time delay between the lens images of 47^{+5}_{-6} d (68%) or 47^{+12}_{-9} d (95%). We use two different methods to derive the time delay; both give the same result within the errors. We estimate an additional systematic error between -8 and +7 d. If the mass distribution of lens galaxy can be described by an isothermal model (Koopmans, de Bruyn & Jackson 1998), this time delay would give a value for the Hubble parameter, H_0=57^{+14}_{-11} (95% statistical) ^{+26}_{-15} (systematic) km/s/Mpc (Omega_m=1 and Omega_Lambda=0). Similarly, the Modified-Hubble-Profile mass model would give H_0=74^{+18}_{-15} (95% statistical) ^{+22}_{-22} (systematic) km/s/Mpc. For Omega_m=0.3 and Omega_Lambda=0.7, these values increase by 5.4%. ... (ABRIDGED)

Radio observations with the Square Kilometer Array (SKA) provide the agility, sensitivity, and spectral coverage to trace the evolution of the size, shape and spectra of gamma-ray burst (GRB) remnants from the earliest moments on. In the first hours to days after the burst a major, and unique tool, will be provided through the study of dynamically evolving radio spectra, caused by diffractive and refractive scintillation in the local ISM. Simultaneous observations from optically thin to optically thick frequencies will provide strong constraints on any model for the GRB remnant. SKA will also allow for extremely rapid (within minutes, if not seconds) follow-up observations by electronic steering of the array. SKA observations have the sensitivity to detect GRB afterglows out to redshifts of 10 or greater. They therefore allow studies of the high redshift universe, and measure the massive star formation (and massive star death rate!) history of the universe, unbiased by optical obscuration. Wide field surveys may detect thousands of GRBs within a few days of observing (depending on the amount of collimation, the width of the GRB luminosity function, the number of GRBs like GRB 980425/SN 1998bw). There would be new ones, but many more would be old, fading afterglows, allowing the determination of the multi-variate GRB-optical-radio luminosity function of GRBs. This should elucidate the relation of the `faint' ones, like GRB 980425, to the luminous ones like GRB 970508. Such surveys will also tell us whether GRB fireballs (and which ones) are collimated into jets.

We present the results from recent VLA 8.5-GHz and WSRT 1.4 and 4.9-GHz monitoring campaigns of the CLASS gravitational lens B1600+434 and show how the observed variations argue strongly in favor of microlensing by MACHOs in the halo of a dark-matter dominated edge-on disk galaxy at z=0.4. The population of flat-spectrum radio sources with micro-Jy flux-densities detected with the Square-Kilometer-Array is expected to have dimensions of micro-arcsec. They will therefore vary rapidly as a result of Galactic scintillation (diffractive and refractive). However, when positioned behind distant galaxies they will also show variations due to microlensing, even more strongly than in the case of B1600+434. Relativistic or superluminal motion in these background sources typically leads to temporal variations on time scales of days to weeks. Scintillation and microlensing can be distinguished, and separated, by their different characteristic time scales and the frequency dependence of their modulations. Monitoring studies with Square-Kilometer-Array at GHz frequencies will thus probe both microscopic and macroscopic properties of dark matter and its mass-function as a function of redshift, information very hard to obtain by any other method.

We report on the discovery of a source which exhibits over 300% amplitude changes in radio flux density on the period of hours. This source, J1819+3845, is the most extremely variable extragalactic source known in the radio sky. We believe these properties are due to interstellar scintillation, and show that the source must emit at least 55% of its flux density within a radius of <16 microarcseconds at 5GHz. The apparent brightness temperature is > 5.10^{12} K, and the source may be explained by a relativistically moving source with a Doppler factor ~15. The scattering occurs predominantly in material only a few tens of parsecs from the earth, which explains its unusually rapid variability. If the source PKS 0405-385 (Kedziora-Chudczer et al 1997) is similarly affected by local scattering material, Doppler factors of ~1000 are not required to explain this source. The discovery of a second source whose properties are well modeled by interstellar scintillation strengthens the argument for this as the cause for much of the variations seen in intra-day variables (IDV).

We present a study of the peculiar radio galaxy B 1834+620. It is characterised by the presence of a 420-kpc large edge-brightened radio source which is situated within, and well aligned with, a larger (1.66 Mpc) radio source. Both sources apparently originate in the same host galaxy, which has a R_s-magnitude of 19.7 and a redshift of 0.5194, as determined from the strong emission-lines in the spectrum. We have determined the rotation measures towards this source, as well as the radio spectral energy distribution of its components. The radio spectrum of the large outer source is steeper than that of the smaller inner source. The radio core has a spectrum that peaks at a frequency of a few GHz. The rotation measures towards the four main components are quite similar, within $\sim\!2$ rad m$^{-2}$ of 58 rad m$^{-2}$. They are probably largely galactic in origin. We have used the presence of a bright hotspot in the northern outer lobe to constrain the advance velocity of the inner radio lobes to the range between 0.19c and 0.29c, depending on the orientation of the source. This corresponds to an age of this structure in the range between 2.6 and 5.8 Myr. We estimate a density of the ambient medium of the inner lobes of $\la 1.6 \times 10^{-30}$ gr\,cm$^{-3}$ (particle density $\la 8 \times 10^{-7}$ cm$^{-3}$). A low ambient density is further supported by the discrepancy between the large optical emission-line luminosity of the host galaxy and the relatively low radio power of the inner lobes.

We present four Mpc-sized radio galaxies which consist of a pair of double-lobed radio sources, aligned along the same axis, and with a coinciding radio core. We have called these peculiar radio sources `double-double' radio galaxies (DDRG) and propose a general definition of such sources: A `double-double' radio galaxy consists of a pair of double radio sources with a common centre. Furthermore, the two lobes of the inner radio source must have a clearly extended, edge-brightened radio morphology. Adopting this definition we find several other candidate DDRGs in the literature. We find that in all sources the smaller (inner) pair of radio lobes is less luminous than the larger (outer) pair, and that the ratio of 1.4-GHz flux density of these two pairs appears to be anti-correlated with the projected linear size of the inner source. Also, the outer radio structures are large, exceeding 700 kpc. We discuss possible formation scenarios of the DDRGs, and we conclude that an interruption of the jet-forming central activity is the most likely mechanism. For one of our sources (B 1834+620) we have been able to observationally constrain the length of time of the interruption to a few Myr. We discuss several scenarios for the cause of the interruption and suggest multiple encounters between interacting galaxies as a possibility. Finally, we discuss whether such interruptions help the formation of extremely large radio sources.

Giant radio sources form the linear size extreme of the extra-galactic radio source population. Using the WENSS survey, we have selected a complete sample of these sources. We have investigated the properties of their radio structures. We find, among other things, that these sources are old (50-100 Myr) and have higher advance velocities than smaller sources of similar radio power. We find pressure gradients in their radio lobes, suggesting that the lobes are still overpressured with respect to the environment. Further, we find no evidence for a cosmological evolution of the radio lobe pressures with increasing redshift, at least up to $z\sim 0.4$, other than that caused by selection effects. We argue that a much fainter sample of giant sources than currently available is needed to constrain the pressure in their environments, the IGM. Another extremely important discovery is that of a population of radio sources with a so-called `double-double' structure, i.e. that of a small two-sided radio source embedded inside a much larger two-sided structure. We argue that these sources result from an interrupted central jet-forming activity. As such, they are the most convincing examples of radio sources with a history of interrupted activity, yet. Since the inner lobes advance within the outer lobes, high resolution low frequency ($\sim 200$ MHz) polarization studies may reveal the constituents of radio lobes and cocoons. We thus argue for a SKA design that can provide low-frequency images at arcsec resolution, but which is also sensitive to structures as large as a few tens of arcminute on the sky.

We present the first unambiguous case of external variability of an extra-galactic radio source, the CLASS gravitational lens B1600+434, consisting of two images lensed by an edge-on disk galaxy. The VLA 8.5-GHz difference light curve of the lens images shows external variability at the 14.6-sigma confidence level. Although the current single-frequency VLA observations cannot conclusively exclude scintillation, several lines of evidence show that it does not dominate the short-term variability. This is supported by an ongoing three-frequency WSRT campaign. If scintillation is the dominant source of this external variability, it would require very different properties of the Galactic ionized ISM towards the two lens images. Microlensing of a superluminal jet component can explain both the rms and time-scale of variability, but requires the halo of the lens galaxy to be filled with >=0.5-Msun MACHOs. The current data appears to support microlensing as the dominant source of the observed external variability.

The third phase of the Cosmic Lens All-Sky Survey (CLASS) has recently been completed, bringing the total number of sources imaged to over 15000 in the CLASS and JVAS combined survey. In the VLA observations carried out in March and April of 1998, two new candidate lensed systems were discovered: CLASS B1152+199 and B1359+154. B1152+199 is a 1.6 arcsecond double, with a background quasar at z=1.019 lensed by a foreground galaxy at z=0.439. The relatively flat radio spectra of the lensed images, combined with a previous ROSAT detection of the source, make B1152+199 a strong candidate for time delay studies at both radio and X-ray wavelengths. B1359+154 is a quadruply lensed quasar at z=3.235, with a maximum image separation of 1.7 arcseconds. As yet, the redshift of the lensing object in this system is undetermined. The steep spectral index of the source suggests that B1359+154 will not exhibit strong variability, and is therefore unlikely to be useful for determining the Hubble constant from measured time delays.

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

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 $\Omega_{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.

We present spectroscopic observations of a sample of faint Gigahertz Peaked Spectrum (GPS) radio sources drawn from the Westerbork Northern Sky Survey (WENSS). Redshifts have been determined for 19 (40%) of the objects. The optical spectra of the GPS sources identified with low redshift galaxies show deep stellar absorption features. This confirms previous suggestions that their optical light is not significantly contaminated by AGN-related emission, but is dominated by a population of old (>9 Gyr) and metal-rich (>0.2 [Fe/H]) stars, justifying the use of these (probably) young radio sources as probes of galaxy evolution. The optical spectra of GPS sources identified with quasars are indistinguishable from those of flat spectrum quasars, and clearly different from the spectra of Compact Steep Spectrum (CSS) quasars. The redshift distribution of the GPS quasars in our radio-faint sample is comparable to that of the bright samples presented in the literature, peaking at z ~ 2-3. It is unlikely that a significant population of low redshift GPS quasars is missed due to selection effects in our sample. We therefore claim that there is a genuine difference between the redshift distributions of GPS galaxies and quasars, which, because it is present in both the radio-faint and bright samples, can not be due to a redshift-luminosity degeneracy. It is therefore unlikely that the GPS quasars and galaxies are unified by orientation, unless the quasar opening angle is a strong function of redshift. We suggest that the GPS quasars and galaxies are unrelated populations and just happen to have identical observed radio-spectral properties, and hypothesise that GPS quasars are a sub-class of flat spectrum quasars.

Gamma-ray bursts (GRBs) are thought to result from the interaction of an extremely relativistic outflow interacting with a small amount of material surrounding the site of the explosion. Multi-wavelength observations covering the gamma-ray to radio wavebands allow investigations of this `fireball' model. On 23 January 1999 optical emission was detected while the gamma-ray burst was still underway. Here we report the results of gamma-ray, optical/infra-red, sub-mm, mm and radio observations of this burst and its afterglow, which indicate that the prompt and afterglow emissions from GRB 990123 are associated with three distinct regions in the fireball. The afterglow one day after the burst has a much lower peak frequency than those of previous bursts; this explains the short-lived nature of the radio emission, which is not expected to reappear. We suggest that such differences reflect variations in the magnetic-field strengths in the afterglow emitting regions.

The reionization of the Universe is expected to leave a signal in the form of a sharp step in the spectrum of the sky. If reionization occurs at 5 < z < 20, a feature should appear in the radio sky at 70 - 240 MHz due to redshifted HI 21-cm line emission, accompanied by another feature in the optical/near-IR at 0.7 - 2.6 micron due to hydrogen recombination radiation. The expected amplitude is well above fundamental detection limits, and the sharpness of the feature may make it distinguishable from variations due to terrestrial, galactic and extragalactic foregrounds. Because this is essentially a continuum measurement of a signal which occurs over the whole sky, relatively small telescopes may suffice for detection in the radio. In the optical/near-IR, a space telescope is needed with the lowest possible background conditions, since the experiment will be severely background-limited.

We present the discovery of a new gravitational lens system with two compact radio images separated by 0.701+-0.001 arcsec. The lens system was discovered in the Cosmic Lens All Sky Survey (CLASS) as a flat spectrum radio source. Both radio components show structure in a VLBA 8.4 GHz radio image. No further extended structure is seen in either the VLA, MERLIN or VLBA images. Hubble Space Telescope (HST) WFPC2 images in F555W and F814W show two extended objects close to the radio components, which we identify as two lens galaxies. Their colours and mass-to-light ratios seem to favour two late-type spiral galaxies at relatively high redshifts (z_d>~0.5). Faint emission is also detected at positions corresponding to the radio images. A two-lens mass model can explain the observed VLBA structure. The best fit model has a reduced chi^2 of 1.1. The relative positions of the VLBA subcomponents are reproduced within 0.08 mas, the flux density ratios within 0.19. We also reproduce the position angle and separation of the two VLBA subcomponents in A and B within the observational errors, which we consider strong evidence for the validity of the lens model. Moreover, we find a surface density axis ratio of 0.74^{+0.10}_{-0.12} for the primary lens (G1), consistent with the surface brightness axis ratio of 0.69+-0.15. Also, the surface density position angle of 4.9^{+28.2}_{-22.4} degree of G1 compares well with the -6+-13 degree position angle of the surface brightness distribution. The errors indicate the 99 per cent confidence interval.

We have discovered a new gravitational lens in the Cosmic Lens All-Sky Survey (CLASS). B2045+265 is a four-image system with a maximum separation of 1.9 arcsec. A fifth radio component is detected, but its radio spectrum and its positional coincidence with infrared emission from the lensing galaxy strongly suggests that it is the radio core of the lensing galaxy. This implies that the B2045+265 system consists of a flat-spectrum radio source being lensed by another flat-spectrum radio source. Infrared images detect the lensed images of the background source and the lensing galaxy. The lensed images have relative positions and flux densities that are consistent with those seen at radio wavelengths. Spectra of the system reveal a lens redshift of z_l=0.8673 and a source redshift of z_s=1.28. The image splitting and system redshifts imply that the projected mass inside the Einstein radius of the lensing galaxy is M_E=4.7x10^{11} h^{-1} M_sun. An estimate of the light emitted inside the Einstein radius gives a mass-to-light ratio in the rest frame B band of (M/L_B)_E = 20 h (M/L_B)_sun. Both the mass and mass-to-light ratio are higher than what is seen in nearby Sa galaxies. In fact, the implied rotation velocity for the lensing galaxy is two to three times higher than what is seen in nearby spirals. The large projected mass inside the Einstein ring radius may be the result of a significant amount of dark matter in the system, perhaps from a compact group of galaxies; however, it may also arise from a misidentification of the source redshift. A simple model of the gravitational potential of the lens reproduces the image positions well, but further modeling is required to satisfy the constraints from the image flux density ratios. With further observations and modeling, this lens may yield an estimate of H_0.

We present the results of new observations of the enigmatic radio source B1144+352 with the WSRT at 1.4 GHz. The WENSS and NVSS surveys have revealed faint extended radio structures on Mpc-scales. We have investigated these large-scale radio components. One of the radio structures has a leading hotspot and we conclude that this structure is a large radio lobe associated with the GPS source. The other contains two separate radio sources which are superposed on the sky: A low-power radio source at a similar redshift (z=0.065) to the GPS host galaxy and an extended radio lobe, which we believe is associated with the GPS host galaxy. Further, the 1.4-GHz flux density of the core has decreased by 70 mJy between 1994 and 1997. We have also retrieved and analysed unpublished archival ROSAT HRI data of B1144+352. The source has been detected and appears to be slightly extended in X-rays. The ROSAT data have not detected a large extended X-ray halo surrounding the host galaxy, limiting the luminosity of the intra-cluster gas within a radius of 250 kpc to < 2.3 x 10^42 erg/s. B1144+352 is the second GPS galaxy known to be associated with a Mpc-sized radio source, the other being B1245+676. We argue that the observed structures must be the result of an interrupted central jet-activity, and that as such these sources may well be the progenitors of double-double radio galaxies.

A sample of 47 faint Gigahertz Peaked Spectrum (GPS) radio sources selected from the Westerbork Northern Sky Survey (WENSS, Rengelink et al. 1997), has been imaged in the optical and near infrared, resulting in an identification fraction of 87%. The R-I and R-K colours of the faint optical counterparts are as expected for passively evolving elliptical galaxies, assuming that they follow the R band Hubble diagram as determined for radio-bright GPS galaxies. We have found evidence that the radio spectral properties of the GPS quasars are different from those of GPS galaxies: The observed distribution of radio spectral peak frequencies for GPS sources optically identified with bright stellar objects (presumably quasars) is shifted compared with GPS sources identified with faint or extended optical objects (presumably galaxies), in the sense that a GPS quasar is likely to have a higher peak frequency than a GPS galaxy. This means that the true peak frequency distribution is different for the GPS galaxies and quasars, because the sample selection effects are independent of optical identification. The correlation between peak frequency and redshift as has been suggested for bright sources has not been found in this sample; no correlation exists between R magnitude (and therefore redshift) and peak frequency for the GPS galaxies. We therefore believe that the claimed correlation is actually caused by the dependence of the peak frequency on optical host, because the GPS galaxies are in general at lower redshifts than the quasars. The difference in the peak frequency distributions of the GPS galaxies and quasars is further evidence against the hypothesis that they form a single class of object.

We present the results of a detailed analysis of the newly discovered 15 arcmin large FRII-type radio galaxy WNB 0313+683. It has a redshift of 0.0901 +/- 0.0002 and a projected linear size of 2.0 Mpc. Using multi-frequency radio data from the WENSS and NVSS surveys, the VLA, the WSRT and the 100m Effelsberg telescope, we have derived the age of the source (1.4 +/- 0.1 x 10^8 yrs) and the particle densities around the two radio lobes (>1.6 x 10^(-6) per cubic cm for the southern lobe and >5.8 x 10^(-7) per cubic cm for the northern). From optical spectroscopy we find that the H$\alpha$ emission line has a broad component, and that the extinction is large (with colour index E(B-V) = 0.98 +/- 0.10 mag), but is probably mostly galactic in origin (b= +9.8 degrees). Further, the radio core has an inverted spectrum with a spectral index of +0.42 +/- 0.03, suggesting a very compact structure. At 10.5 GHz, it contributes 25% of the total flux density of the source, which is extremely high for radio galaxies. WNB 0313+683 has a very large optical emission line flux with respect to its estimated jet power, when compared with the correlation between these two properties found by Rawlings & Saunders (1991). We argue that this, together with the relatively high radio power and the inverted radio spectrum of the radio core, is suggestive of a new phase of radio activity in this source.

We report on Westerbork 1.4 GHz radio observations of the radio counterpart to $\gamma$-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\gamma_{e}\propto\gamma_{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.

The Westerbork Northern Sky Survey (WENSS) has been used to select a sample of Gigahertz Peaked Spectrum (GPS) radio sources at flux densities one to two orders of magnitude lower than bright GPS sources investigated in earlier studies. Sources with inverted spectra at frequencies above 325 MHz have been observed with the WSRT at 1.4 and 5 GHz and with the VLA at 8.6 and 15 GHz to select genuine GPS sources. This has resulted in a sample of 47 GPS sources with peak frequencies ranging from ~500 MHz to >15 GHz, and peak flux densities ranging from ~40 to ~900 mJy. Counts of GPS sources in our sample as a function of flux density have been compared with counts of large scale sources from WENSS scaled to 2 GHz, the typical peak frequency of our GPS sources. The counts can be made similar if the number of large scale sources at 2 GHz is divided by 250, and their flux densities increase by a factor of 10. On the scenario that all GPS sources evolve into large scale radio sources, these results show that the lifetime of a typical GPS source is ~250 times shorter than a typical large scale radio source, and that the source luminosity must decrease by a factor of ~10 in evolving from GPS to large scale radio source. However, we note that the redshift distributions of GPS and large scale radio sources are different and that this hampers a direct and straightforward interpretation of the source counts. Further modeling of radio source evolution combined with cosmological evolution of the radio luminosity function for large sources is required.

A relativistic radio component, which moves in a direction close to the sky plane, will increase in flux density when it decelerates. This effect is the basis for the qualitative model for GPS galaxies we present in this paper, which can explain their low-variability convex spectrum, their compact double or compact symmetric morphology, and the lack of GPS quasars at similar redshifts. Components are expelled from the nucleus at relativistic speeds at a large angle to the line of sight, and are decelerated (eg. by ram-pressure or entrainment of the external gas) before contributing to a mini-lobe. The young components are Doppler boosted in the direction of motion but appear fainter for the observer. The non-relativistic mini-lobes dominate the structure and are responsible for the low variability in flux density and the convex radio spectrum as well as the compact double angular morphology. Had the same source been orientated at a small angle to the line of sight, the young components would be boosted in the observer's direction resulting in a flat and variable radio spectrum at high frequencies. Hence the characteristic convex spectrum of a GPS source would not be seen. These sources at small angles to the line of sight are probably identified with quasars, and are not recognized as GPS sources, but are embedded in the large population of flat spectrum variable quasars and BL Lac objects. This leads to a deficiency in GPS/CSOs identified with quasars.

We present new observations of the gravitational lens (GL) system B1600+434, strongly suggesting that the lens is an edge-on spiral galaxy. These observations are used to constrain the mass model of the system, in particular the oblateness and velocity dispersion of the dark matter halo around the lensing galaxy. From an analytical model we find a lower limit on the halo oblateness $q_\rmn{halo}=(c/a)_{\rho}\ga0.4$; more detailed numerical models give a lower limit $q_\rmn{h}\ga0.5$. We determine an average halo velocity dispersion of $\sigma_{halo}=190\pm 15$ km/s over all non-singular isothermal elliptical (NIE) halo models. Constraining the models to larger and more massive disks, decreases this average by only 10 km/s. A lower limit of $\sigma_{halo}\ga 150$ km/s is found, even for disk masses larger then the mass inside the Einstein radius. This lower limit indicates the need for a massive dark matter halo, contributing at least half of the mass inside the Einstein radius. Time delay calculations give $(54\pm3)/h_{50}$ days for the NIE halo model and $(70\pm4)/h_{50}$ days for the modified Hubble profile (MHP) halo model. Although the time delay for both NIE and MHP halo models is well constrained on our parameter grid, it strongly depends on the halo surface density profile. We furthermore find that the presence of a flat luminous mass distribution can severely alter the statistical properties of the lens.

We report the discovery of the most complex arcsec-scale radio gravitational lens system yet known. B1933+503 was found during the course of the CLASS survey and MERLIN and VLA radio maps reveal up to 10 components. Four of these are compact and have flat spectra; the rest are more extended and have steep spectra. The background lensed object appears to consist of a flat spectrum core (quadruply imaged) and two compact "lobes" symmetrically disposed relative to the core. One of the lobes is quadruply imaged while the other is doubly imaged. An HST observation of the system with the WFPC2 shows a galaxy with an axial ratio of 0.5, but none of the images of the background object are detected. A redshift of 0.755 has been measured for the lens galaxy.

We report the discovery, using NICMOS on the Hubble Space Telescope, of an arcsecond-diameter Einstein ring in the gravitational lens system B1938+666. The lensing galaxy is also detected, and is most likely an early-type. Modelling of the ring is presented and compared with the radio structure from MERLIN maps. We show that the Einstein ring is consistent with the gravitational lensing of an extended infrared component, centred between the two radio components.

We report on Westerbork 840 MHz, 1.4 and 5 GHz radio observations of the improved IPN-WFC error box of the gamma-ray burst GRB 970111, between 26.4 hours and 120 days after the event onset. In the ~ 13 sq arcmin area defined by the IPN (BATSE and Ulysses) annulus and the published refined BeppoSAX Wide Field Camera error box we detected no steady sources brighter than 0.56 mJy (4 sigma), and no varying radio emission, down to 1.0 mJy (4 sigma). We also report on B, V, R and I band observations of the error box with the 4.2 m William Herschel Telescope at La Palma.

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$\alpha$ 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.

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.