High-precision cosmological probes have revealed a small but significant tension between the parameters measured with different techniques, one of which is based on time delays in gravitational lenses. We discuss a new way of using time delays for cosmology, taking advantage of the extreme precision expected for lensed Fast Radio Bursts, short flashes of radio emission originating at cosmological distances. With coherent methods the achievable precision is sufficient to even measure how time delays change over the months and years, which can also be interpreted as differential redshifts between the images. It turns out that uncertainties arising from the unknown mass distribution of gravitational lenses can be eliminated by combining time delays with their time derivatives. Other effects, most importantly relative proper motion, can be measured accurately and separated from the cosmological effects. With a mock sample of simulated lenses we show that strong constraints on cosmological parameters are potentially possible. Finally the lensed images can be used as galactic interferometer to resolve structures and motions of the burst sources with insanely high resolution and help revealing their currently unknown physical nature.

We present low-radio-frequency follow-up observations of AT 2017gfo, the electromagnetic counterpart of GW170817, which was the first binary neutron star merger to be detected by Advanced LIGO-Virgo. These data, with a central frequency of 144 MHz, were obtained with LOFAR, the Low-Frequency Array. The maximum elevation of the target is just 13.7 degrees when observed with LOFAR, making our observations particularly challenging to calibrate and significantly limiting the achievable sensitivity. On time-scales of 130-138 and 371-374 days after the merger event, we obtain 3$\sigma$ upper limits for the afterglow component of 6.6 and 19.5 mJy beam$^{-1}$, respectively. Using our best upper limit and previously published, contemporaneous higher-frequency radio data, we place a limit on any potential steepening of the radio spectrum between 610 and 144 MHz: the two-point spectral index $\alpha^{610}_{144} \gtrsim -2.5$. We also show that LOFAR can detect the afterglows of future binary neutron star merger events occurring at more favourable elevations.

This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported.

The four persistent radio sources in the northern sky with the highest flux density at metre wavelengths are Cassiopeia A, Cygnus A, Taurus A, and Virgo A; collectively they are called the A-team. Their flux densities at ultra-low frequencies (<100 MHz) can reach several thousands of janskys, and they often contaminate observations of the low-frequency sky by interfering with image processing. Furthermore, these sources are foreground objects for all-sky observations hampering the study of faint signals, such as the cosmological 21 cm line from the epoch of reionisation. We aim to produce robust models for the surface brightness emission as a function of frequency for the A-team sources at ultra-low frequencies. These models are needed for the calibration and imaging of wide-area surveys of the sky with low-frequency interferometers. This requires obtaining images at an angular resolution better than 15 arcsec with a high dynamic range and good image fidelity. We observed the A-team with the Low Frequency Array (LOFAR) at frequencies between 30 MHz and 77 MHz using the Low Band Antenna (LBA) system. We reduced the datasets and obtained an image for each A-team source. The paper presents the best models to date for the sources Cassiopeia A, Cygnus A, Taurus A, and Virgo A between 30 MHz and 77 MHz. We were able to obtain the aimed resolution and dynamic range in all cases. Owing to its compactness and complexity, observations with the long baselines of the International LOFAR Telescope will be required to improve the source model for Cygnus A further.

We report on the multi-frequency timing observations of 21 pulsars discovered in the LOFAR Tied-Array All-Sky Survey (LOTAAS). The timing data were taken at central frequencies of 149 MHz (LOFAR) as well as 334 and 1532 MHz (Lovell Telecope). The sample of pulsars includes 20 isolated pulsars and the first binary pulsar discovered by the survey, PSR J1658$+$3630. We modelled the timing properties of the pulsars, which showed that they have, on average, larger characteristic ages. We present the pulse profiles of the pulsars across the three observing bands, where PSR J1643$+$1338 showed profile evolution that appears not to be well-described by the radius-to-frequency-mapping model. Furthermore, we modelled the spectra of the pulsars across the same observing bands, using a simple power law, and found an average spectral index of $-1.9 \pm 0.5$. Amongst the pulsars studied here, PSR J1657$+$3304 showed large flux density variations of a factor of 10 over 300 days, as well as mode changing and nulling on timescales of a few minutes. We modelled the rotational and orbital properties of PSR J1658$+$3630, which has a spin period of 33 ms in a binary orbit of 3.0 days with a companion of minimum mass of 0.87$M_{\odot}$, likely a Carbon-Oxygen or Oxygen-Neon-Magnesium type white dwarf. PSR J1658$+$3630 has a dispersion measure of 3.0 pc cm$^{-3}$, making it possibly one of the closest binary pulsars known.

Interstellar scattering is known to broaden distant objects spatially and temporally. The latter aspect is difficult to analyse, unless the signals carry their own time stamps. Pulsars are so kind to do us this favour. Typically the signature is a broadened image with little or no substructure and a similarly smooth exponential scattering tail in the temporal profile. The case of the pulsar B1508+55 is special: The profile shows additional components that are moving relative to the main pulse with time. We use low-frequency VLBI with LOFAR to test the hypothesis that these components are actually such scattering-induced echoes, by trying to detect the expected angular offset. Using international stations (plus the Kilpisj\"arvi Atmospheric Imaging Receiver Array "KAIRA") and the phased-up core of the LOFAR array, we can do interferometry at high resolution in time and space. This contribution presents a selection of results from an ongoing large-scale monitoring campaign. We can not only detect the offset, but even image a full string of echoes, and relate the positions with delays. What we find is apparently consistent with scattering by highly aligned components in a single screen at a distance of 120 pc. Further investigations will improve our understanding of the scattering process as basis of using the scattering-induced subimages as arms of a giant interstellar interferometer with insanely high resolution.

Faraday rotation of polarized emission from pulsars measured at radio frequencies provides a powerful tool to investigate the interstellar and interplanetary magnetic fields. However, besides being sensitive to the astrophysical media, pulsar observations in radio are affected by the highly time-variable ionosphere. In this article, the amount of ionospheric Faraday rotation has been computed by assuming a thin layer model. For this aim, ionospheric maps of the free electron density (based on Global Positioning System data) and semi-empirical geomagnetic models are needed. Through the data of five highly polarized pulsars observed with the individual German LOw-Frequency ARray stations, we investigate the performances of the ionospheric modelling. In addition, we estimate the parameters of the systematics and the correlated noise generated by the residual unmodelled ionospheric effects, and show the comparison of the different free-electron density maps. For the best ionospheric maps, we have found that the rotation measure corrections on one-year timescales after subtraction of diurnal periodicity are accurate to $\sim$ 0.06--0.07 rad m$^{-2}$.

Type II radio bursts are evidence of shocks in the solar atmosphere and inner heliosphere that emit radio waves ranging from sub-meter to kilometer lengths. These shocks may be associated with CMEs and reach speeds higher than the local magnetosonic speed. Radio imaging of decameter wavelengths (20-90 MHz) is now possible with LOFAR, opening a new radio window in which to study coronal shocks that leave the inner solar corona and enter the interplanetary medium and to understand their association with CMEs. To this end, we study a coronal shock associated with a CME and type II radio burst to determine the locations at which the radio emission is generated, and we investigate the origin of the band-splitting phenomenon.

Magnetic fields play a pivotal role in the physics of interstellar medium in galaxies, but there are few observational constraints on how they evolve across cosmic time. Spatially resolved synchrotron polarization maps at radio wavelengths reveal well-ordered large-scale magnetic fields in nearby galaxies that are believed to grow from a seed field via a dynamo effect. To directly test and characterize this theory requires magnetic field strength and geometry measurements in cosmologically distant galaxies, which are challenging to obtain due to the limited sensitivity and angular resolution of current radio telescopes. Here, we report the cleanest measurements yet of magnetic fields in a galaxy beyond the local volume, free of the systematics traditional techniques would encounter. By exploiting the scenario where the polarized radio emission from a background source is gravitationally lensed by a foreground galaxy at z = 0.439 using broadband radio polarization data, we detected coherent $\mu$G magnetic fields in the lensing disk galaxy as seen 4.6 Gyrs ago, with similar strength and geometry to local volume galaxies. This is the highest redshift galaxy whose observed coherent magnetic field property is compatible with a mean-field dynamo origin.

Context. The Sun is an active star that produces large-scale energetic events such as solar flares and coronal mass ejections and numerous smaller-scale events such as solar jets. These events are often associated with accelerated particles that can cause emission at radio wavelengths. The reconfiguration of the solar magnetic field in the corona is believed to be the cause of the majority of solar energetic events and accelerated particles. Aims. Here, we investigate a bright J-burst that was associated with a solar jet and the possible emission mechanism causing these two phenomena. Methods. We used data from the Solar Dynamics Observatory (SDO) to observe a solar jet, and radio data from the Low Frequency Array (LOFAR) and the Nan\c{c}ay Radioheliograph (NRH) to observe a J-burst over a broad frequency range (33-173 MHz) on 9 July 2013 at ~11:06 UT. Results. The J-burst showed fundamental and harmonic components and it was associated with a solar jet observed at extreme ultraviolet wavelengths with SDO. The solar jet occurred at a time and location coincident with the radio burst, in the northern hemisphere, and not inside a group of complex active regions in the southern hemisphere. The jet occurred in the negative polarity region of an area of bipolar plage. Newly emerged positive flux in this region appeared to be the trigger of the jet. Conclusions. Magnetic reconnection between the overlying coronal field lines and the newly emerged positive field lines is most likely the cause of the solar jet. Radio imaging provides a clear association between the jet and the J-burst which shows the path of the accelerated electrons.

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree.

We analyse new observations with the International Low Frequency Array (LOFAR) telescope, and archival data from the Multi-Element Radio Linked Interferometer Network (MERLIN) and the Karl G. Jansky Very Large Array (VLA). We model the spatially resolved radio spectrum of Arp 220 from 150 MHz to 33 GHz. We present an image of Arp 220 at 150 MHz with resolution $0.65''\times0.35''$, sensitivity 0.15 mJy beam$^{-1}$, and integrated flux density $394\pm59$ mJy. More than 80% of the detected flux comes from extended ($6''\approx$2.2 kpc) steep spectrum ($\alpha=-0.7$) emission, likely from star formation in the molecular disk surrounding the two nuclei. We find elongated features extending $0.3''$ (110 pc) and $0.9''$ (330 pc) from the eastern and western nucleus respectively, which we interpret as evidence for outflows. The extent of radio emission requires acceleration of cosmic rays far outside the nuclei. We find that a simple three component model can explain most of the observed radio spectrum of the galaxy. When accounting for absorption at 1.4 GHz, Arp 220 follows the FIR/radio correlation with $q=2.36$, and we estimate a star formation rate of 220 M$_\odot\text{yr}^{-1}$. We derive thermal fractions at 1 GHz of less than 1% for the nuclei, which indicates that a major part of the UV-photons are absorbed by dust. International LOFAR observations shows great promise to detect steep spectrum outflows and probe regions of thermal absorption. However, in LIRGs the emission detected at 150 MHz does not necessarily come from the main regions of star formation. This implies that high spatial resolution is crucial for accurate estimates of star formation rates for such galaxies at 150 MHz.

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range.

Air showers induced by cosmic rays create nanosecond pulses detectable at radio frequencies. These pulses have been measured successfully in the past few years at the LOw Frequency ARray (LOFAR) and are used to study the properties of cosmic rays. For a complete understanding of this phenomenon and the underlying physical processes, an absolute calibration of the detecting antenna system is needed. We present three approaches that were used to check and improve the antenna model of LOFAR and to provide an absolute calibration of the whole system for air shower measurements. Two methods are based on calibrated reference sources and one on a calibration approach using the diffuse radio emission of the Galaxy, optimized for short data-sets. An accuracy of 19% in amplitude is reached. The absolute calibration is also compared to predictions from air shower simulations. These results are used to set an absolute energy scale for air shower measurements and can be used as a basis for an absolute scale for the measurement of astronomical transients with LOFAR.

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 {\deg} $\pm$ 25 {\deg}. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles.

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

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.

We present VLA detections of radio emission in four four-image gravitational lens systems with quasar sources: HS0810+2554, RXJ0911+0511, HE0435$-$1223 and SDSSJ0924+0219, and e-MERLIN observations of two of the systems. The first three are detected at a high level of significance, and SDSS J0924+0219 is detected. HS0810+2554 is resolved, allowing us for the first time to achieve 10-mas resolution of the source frame in the structure of a radio quiet quasar. The others are unresolved or marginally resolved. All four objects are among the faintest radio sources yet detected, with intrinsic flux densities in the range 1-5$\mu$Jy; such radio objects, if unlensed, will only be observable routinely with the Square Kilometre Array. The observations of HS0810+2554, which is also detected with e-MERLIN, strongly suggest the presence of a mini-AGN, with a radio core and milliarcsecond scale jet. The flux densities of the lensed images in all but HE0435-1223 are consistent with smooth galaxy lens models without the requirement for smaller-scale substructure in the model, although some interesting anomalies are seen between optical and radio flux densities. These are probably due to microlensing effects in the optical.

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.

Using observations obtained with the LOw Fequency ARray (LOFAR), the Westerbork Synthesis Radio Telescope (WSRT) and archival Very Large Array (VLA) data, we have traced the radio emission to large scales in the complex source 4C 35.06 located in the core of the galaxy cluster Abell 407. At higher spatial resolution (~4"), the source was known to have two inner radio lobes spanning 31 kpc and a diffuse, low-brightness extension running parallel to them, offset by about 11 kpc (in projection). At 62 MHz, we detect the radio emission of this structure extending out to 210 kpc. At 1.4 GHz and intermediate spatial resolution (~30"), the structure appears to have a helical morphology. We have derived the characteristics of the radio spectral index across the source. We show that the source morphology is most likely the result of at least two episodes of AGN activity separated by a dormant period of around 35 Myr. The AGN is hosted by one of the galaxies located in the cluster core of Abell 407. We propose that it is intermittently active as it moves in the dense environment in the cluster core. Using LOFAR, we can trace the relic plasma from that episode of activity out to greater distances from the core than ever before. Using the the WSRT, we detect HI in absorption against the center of the radio source. The absorption profile is relatively broad (FWHM of 288 km/s), similar to what is found in other clusters. Understanding the duty cycle of the radio emission as well as the triggering mechanism for starting (or restarting) the radio-loud activity can provide important constraints to quantify the impact of AGN feedback on galaxy evolution. The study of these mechanisms at low frequencies using morphological and spectral information promises to bring new important insights in this field.

We present measurements of radio emission from cosmic ray air showers that took place during thunderstorms. The intensity and polarization patterns of these air showers are radically different from those measured during fair-weather conditions. With the use of a simple two-layer model for the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study atmospheric electric fields.

Context. The LOw Frequency ARray (LOFAR) radio telescope is a giant digital phased array interferometer with multiple antennas distributed in Europe. It provides discrete sets of Fourier components of the sky brightness. Recovering the original brightness distribution with aperture synthesis forms an inverse problem that can be solved by various deconvolution and minimization methods Aims. Recent papers have established a clear link between the discrete nature of radio interferometry measurement and the "compressed sensing" (CS) theory, which supports sparse reconstruction methods to form an image from the measured visibilities. Empowered by proximal theory, CS offers a sound framework for efficient global minimization and sparse data representation using fast algorithms. Combined with instrumental direction-dependent effects (DDE) in the scope of a real instrument, we developed and validated a new method based on this framework Methods. We implemented a sparse reconstruction method in the standard LOFAR imaging tool and compared the photometric and resolution performance of this new imager with that of CLEAN-based methods (CLEAN and MS-CLEAN) with simulated and real LOFAR data Results. We show that i) sparse reconstruction performs as well as CLEAN in recovering the flux of point sources; ii) performs much better on extended objects (the root mean square error is reduced by a factor of up to 10); and iii) provides a solution with an effective angular resolution 2-3 times better than the CLEAN images. Conclusions. Sparse recovery gives a correct photometry on high dynamic and wide-field images and improved realistic structures of extended sources (of simulated and real LOFAR datasets). This sparse reconstruction method is compatible with modern interferometric imagers that handle DDE corrections (A- and W-projections) required for current and future instruments such as LOFAR and SKA

Temporal scatter-broadening can seriously affect our ability to find pulsars orbiting the central mass in our Galaxy. Many of these invaluable probes of geometry around the black hole are expected, but none have been found in close orbits so far, possibly as result of strong scattering. The magnetar PSR J1745-2900 discovered in 2013 at a separation of < 3 arcsec is not the optimal type of pulsar for studies of general relativity, but it can be used to investigate the scattering properties so that search strategies can be adapted accordingly. This contribution presents an observation of this magnetar using short baselines between VLBI stations in Europe in a non-standard interferometry mode. The most important goal is determining the distance of the scattering screen, or the distribution of scattering material if not confined to one screen. The analysis is based on phase-binned visibilities that allow measuring the shape of the scattering disk and how it grows with increasing delay over the scattering tail of the pulse profile. Narrow rings growing with the square root of delay are expected for a single thin scattering screen and the preliminary results are indeed consistent with this expectation. This means that most of the angular and temporal broadening is caused by the same and relatively thin scattering screen and that, in contrast to standard models of the interstellar scattering behaviour near the Galactic centre, this screen is located about halfway between the centre and us.

Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.

The nuclear starburst in the nearby galaxy M82 provides an excellent laboratory for understanding the physics of star formation. This galaxy has been extensively observed in the past, revealing tens of radio-bright compact objects embedded in a diffuse free-free absorbing medium. Our understanding of the structure and physics of this medium in M82 can be greatly improved by high-resolution images at low frequencies where the effects of free-free absorption are most prominent. The aims of this study are, firstly, to demonstrate imaging using international baselines of the Low Frequency Array (LOFAR), and secondly, to constrain low-frequency spectra of compact and diffuse emission in the central starburst region of M82 via high-resolution radio imaging at low frequencies. The international LOFAR telescope was used to observe M82 at 110-126MHz and 146-162MHz. Images were obtained using standard techniques from very long baseline interferometry. images were obtained at each frequency range: one only using international baselines, and one only using the longest Dutch (remote) baselines. The 154MHz image obtained using international baselines is a new imaging record in terms of combined image resolution (0.3$"$) and sensitivity ($\sigma$=0.15mJy/beam) at low frequencies ($<327$MHz). We detected 16 objects at 154MHz, six of these also at 118MHz. Four weaker but resolved features are also found: a linear (50pc) filament and three other resolved objects, of which two show a clear shell structure. We do not detect any emission from either supernova 2008iz or from the radio transient source 43.78+59.3. The images obtained using remote baselines show diffuse emission, associated with the outflow in M82, with reduced brightness in the region of the edge-on star-forming disk.

Measuring radio emission from air showers offers a novel way to determine properties of the primary cosmic rays such as their mass and energy. Theory predicts that relativistic time compression effects lead to a ring of amplified emission which starts to dominate the emission pattern for frequencies above ~100 MHz. In this article we present the first detailed measurements of this structure. Ring structures in the radio emission of air showers are measured with the LOFAR radio telescope in the frequency range of 110 - 190 MHz. These data are well described by CoREAS simulations. They clearly confirm the importance of including the index of refraction of air as a function of height. Furthermore, the presence of the Cherenkov ring offers the possibility for a geometrical measurement of the depth of shower maximum, which in turn depends on the mass of the primary particle.

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.

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

The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.

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

Extensive air showers, induced by high energy cosmic rays impinging on the Earth's atmosphere, produce radio emission that is measured with the LOFAR radio telescope. As the emission comes from a finite distance of a few kilometers, the incident wavefront is non-planar. A spherical, conical or hyperbolic shape of the wavefront has been proposed, but measurements of individual air showers have been inconclusive so far. For a selected high-quality sample of 161 measured extensive air showers, we have reconstructed the wavefront by measuring pulse arrival times to sub-nanosecond precision in 200 to 350 individual antennas. For each measured air shower, we have fitted a conical, spherical, and hyperboloid shape to the arrival times. The fit quality and a likelihood analysis show that a hyperboloid is the best parametrization. Using a non-planar wavefront shape gives an improved angular resolution, when reconstructing the shower arrival direction. Furthermore, a dependence of the wavefront shape on the shower geometry can be seen. This suggests that it will be possible to use a wavefront shape analysis to get an additional handle on the atmospheric depth of the shower maximum, which is sensitive to the mass of the primary particle.

We present the first detection of carbon radio recombination line absorption along the line of sight to Cygnus A. The observations were carried out with the LOw Frequency ARray in the 33 to 57 MHz range. These low frequency radio observations provide us with a new line of sight to study the diffuse, neutral gas in our Galaxy. To our knowledge this is the first time that foreground Milky Way recombination line absorption has been observed against a bright extragalactic background source. By stacking 48 carbon $\alpha$ lines in the observed frequency range we detect carbon absorption with a signal-to-noise ratio of about 5. The average carbon absorption has a peak optical depth of 2$\times$10$^{-4}$, a line width of 10 km s$^{-1}$ and a velocity of +4 km s$^{-1}$ with respect to the local standard of rest. The associated gas is found to have an electron temperature $T_{e}\sim$ 110 K and density $n_{e}\sim$ 0.06 cm$^{-3}$. These properties imply that the observed carbon $\alpha$ absorption likely arises in the cold neutral medium of the Orion arm of the Milky Way. Hydrogen and helium lines were not detected to a 3$\sigma$ peak optical depth limit of 1.5$\times$10$^{-4}$ for a 4 km s$^{-1}$ channel width. Radio recombination lines associated with Cygnus A itself were also searched for, but are not detected. We set a 3$\sigma$ upper limit of 1.5$\times$10$^{-4}$ for the peak optical depth of these lines for a 4 km s$^{-1}$ channel width.

The low frequency array (LOFAR), is the first radio telescope designed with the capability to measure radio emission from cosmic-ray induced air showers in parallel with interferometric observations. In the first $\sim 2\,\mathrm{years}$ of observing, 405 cosmic-ray events in the energy range of $10^{16} - 10^{18}\,\mathrm{eV}$ have been detected in the band from $30 - 80\,\mathrm{MHz}$. Each of these air showers is registered with up to $\sim1000$ independent antennas resulting in measurements of the radio emission with unprecedented detail. This article describes the dataset, as well as the analysis pipeline, and serves as a reference for future papers based on these data. All steps necessary to achieve a full reconstruction of the electric field at every antenna position are explained, including removal of radio frequency interference, correcting for the antenna response and identification of the pulsed signal.

The original intensity interferometers were instruments built in the 1950s and 60s by Hanbury Brown and collaborators, achieving milli-arcsec resolutions in visible light without optical-quality mirrors. They exploited a then-novel physical effect, now known as HBT correlation after the experiments of Hanbury Brown and Twiss, and nowadays considered fundamental in quantum optics. Now a new generation of inten- sity interferometers is being designed, raising the possibility of measuring intensity correlations with three or more detectors. Quantum optics predicts some interesting features in higher-order HBT. One is that HBT correlation increases combinatorially with the number of detectors. Signal to noise considerations suggest, that many-detector HBT correlations would be mea- surable for bright masers, but very difficult for thermal sources. But the more modest three-detector HBT correlation seems measurable for bright stars, and would provide image information (namely the bispectrum) not present in standard HBT.

The characteristic outer scale of turbulence and the ratio of the random to ordered components of the magnetic field are key parameters to characterise magnetic turbulence in the interstellar gas, which affects the propagation of cosmic rays within the Galaxy. We provide new constraints to those two parameters. We use the LOw Frequency ARray (LOFAR) to image the diffuse continuum emission in the Fan region at (l,b) (137.0,+7.0) at 80"x70" resolution in the range [146,174] MHz. We detect multi-scale fluctuations in the Galactic synchrotron emission and compute their power spectrum. Applying theoretical estimates and derivations from the literature for the first time, we derive the outer scale of turbulence and the ratio of random to ordered magnetic field from the characteristics of these fluctuations . We obtain the deepest image of the Fan region to date and find diffuse continuum emission within the primary beam. The power spectrum of the foreground synchrotron fluctuations displays a power law behaviour for scales between 100 and 8 arcmin with a slope of (-1.84+/-0.19). We find an upper limit of about 20 pc for the outer scale of the magnetic interstellar turbulence toward the Fan region. We also find a variation of the ratio of random to ordered field as a function of Galactic coordinates, supporting different turbulent regimes. We use power spectra fluctuations from LOFAR as well as earlier GMRT and WSRT observations to constrain the outer scale of turbulence of the Galactic synchrotron foreground, finding a range of plausible values of 10-20 pc. Then, we use this information to deduce lower limits of the ratio of ordered to random magnetic field strength. These are found to be 0.3, 0.3, and 0.5 for the LOFAR, WSRT and GMRT fields considered respectively. Both these constraints are in agreement with previous estimates.

We report on a recent global VLBI experiment in which we study the scatter broadening of pulsars in the spatial and time domain simultaneously. Depending on the distribution of scattering screen(s), geometry predicts that the less spatially broadened parts of the signal arrive earlier than the more broadened parts. This means that over one pulse period the size of the scattering disk should grow from pointlike to the maximum size. An equivalent description is that the pulse profile shows less temporal broadening on the longer baselines. This contribution presents first results that are consistent with the expected expanding rings. We also briefly discuss how the autocorrelations can be used for amplitude calibration. This requires a thorough investigation of the digitisation and the sampler statistics and is not fully solved yet.

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.

The information lost in images of undersampled CCD cameras can be recovered with the technique of `dithering'. A number of subexposures is taken with sub-pixel shifts in order to record structures on scales smaller than a pixel. The standard method to combine such exposures, `Drizzle', averages after reversing the displacements, including rotations and distortions. More sophisticated methods are available to produce, e.g., Nyquist sampled representations of band-limited inputs. While the combined images produced by these methods can be of high quality, their use as input for forward-modelling techniques in gravitational lensing is still not optimal, because the residual artefacts still affect the modelling results in unpredictable ways. In this paper we argue for an overall modelling approach that takes into account the dithering and the lensing without the intermediate product of a combined image. As one building block we introduce an alternative approach to combine dithered images by direct model fitting with a least-squares approach including a regularization constraint. We present tests with simulated and real data that show the quality of the results. The additional effects of gravitational lensing and the convolution with an instrumental point spread function can be included in a natural way, avoiding the possible systematic errors of previous procedures.

Multiply-imaged quasars and AGNs observed in the mid-infrared (MIR) range are commonly assumed to be unaffected by the microlensing produced by the stars in their lensing galaxy. In this paper, we investigate the validity domain of this assumption. Indeed, that premise disregards microlensing of the accretion disc in the MIR range, and does not account for recent progress in our knowledge of the dusty torus. To simulate microlensing, we first built a simplified image of the quasar composed of an accretion disc, and of a larger ring-like torus. The mock quasars are then microlensed using an inverse ray-shooting code. We simulated the wavelength and size dependence of microlensing for different lensed image types and fraction of compact objects projected in the lens. This allows us to derive magnification probabilities as a function of wavelength, as well as to calculate the microlensing-induced deformation of the spectral energy distribution of the lensed images. We find that microlensing variations as large as 0.1 mag are very common at 11 microns (observer-frame). The main signal comes from microlensing of the accretion disc, which may be significant even when the fraction of flux from the disc is as small as 5 % of the total flux. We also show that the torus of sources with Lbol <~ 10^45 erg/s is expected to be noticeably microlensed. Microlensing may thus be used to get insight into the rest near-infrared inner structure of AGNs. Finally, we investigate whether microlensing in the mid-infrared can alter the so-called Rcusp relation that links the fluxes of the lensed images triplet produced when the source lies close to a cusp macro-caustic. This relation is commonly used to identify massive (dark-matter) substructures in lensing galaxies. We find that significant deviations from Rcusp may be expected, which means that microlensing can explain part of the flux ratio problem.

The origin of magnetic fields in the Universe is an open problem in astrophysics and fundamental physics. Polarization observations with the forthcoming large radio telescopes will open a new era in the observation of magnetic fields and should help to understand their origin. At low frequencies, LOFAR (10-240 MHz) will allow us to map the structure of weak magnetic fields in the outer regions and halos of galaxies, in galaxy clusters and in the Milky Way via their synchrotron emission. Even weaker magnetic fields can be measured at low frequencies with help of Faraday rotation measures. A detailed view of the magnetic fields in the local Milky Way will be derived by Faraday rotation measures from pulsars. First promising images with LOFAR have been obtained for the Crab pulsar-wind nebula, the spiral galaxy M51, the radio galaxy M87 and the galaxy clusters A2255 and A2256. With help of the polarimetric technique of "Rotation Measure Synthesis", diffuse polarized emission has been detected from a magnetic bubble in the local Milky Way. Polarized emission and rotation measures were measured for more than 20 pulsars so far.

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.

Faraday rotation measurements using the current and next generation of low-frequency radio telescopes will provide a powerful probe of astronomical magnetic fields. However, achieving the full potential of these measurements requires accurate removal of the time-variable ionospheric Faraday rotation contribution. We present ionFR, a code that calculates the amount of ionospheric Faraday rotation for a specific epoch, geographic location, and line-of-sight. ionFR uses a number of publicly available, GPS-derived total electron content maps and the most recent release of the International Geomagnetic Reference Field. We describe applications of this code for the calibration of radio polarimetric observations, and demonstrate the high accuracy of its modeled ionospheric Faraday rotations using LOFAR pulsar observations. These show that we can accurately determine some of the highest-precision pulsar rotation measures ever achieved. Precision rotation measures can be used to monitor rotation measure variations - either intrinsic or due to the changing line-of-sight through the interstellar medium. This calibration is particularly important for nearby sources, where the ionosphere can contribute a significant fraction of the observed rotation measure. We also discuss planned improvements to ionFR, as well as the importance of ionospheric Faraday rotation calibration for the emerging generation of low-frequency radio telescopes, such as the SKA and its pathfinders.

Cassiopeia A was observed using the Low-Band Antennas of the LOw Frequency ARray (LOFAR) with high spectral resolution. This allowed a search for radio recombination lines (RRLs) along the line-of-sight to this source. Five carbon-alpha RRLs were detected in absorption between 40 and 50 MHz with a signal-to-noise ratio of > 5 from two independent LOFAR datasets. The derived line velocities (v_LSR ~ -50 km/s) and integrated optical depths (~ 13 s^-1) of the RRLs in our spectra, extracted over the whole supernova remnant, are consistent within each LOFAR dataset and with those previously reported. For the first time, we are able to extract spectra against the brightest hotspot of the remnant at frequencies below 330 MHz. These spectra show significantly higher (15-80 %) integrated optical depths, indicating that there is small-scale angular structure on the order of ~1 pc in the absorbing gas distribution over the face of the remnant. We also place an upper limit of 3 x 10^-4 on the peak optical depths of hydrogen and helium RRLs. These results demonstrate that LOFAR has the desired spectral stability and sensitivity to study faint recombination lines in the decameter band.

Some radio pulsars show clear drifting subpulses, in which subpulses are seen to drift in pulse longitude in a systematic pattern. Here we examine how the drifting subpulses of PSR B0809+74 evolve with time and observing frequency. We show that the subpulse period (P3) is constant on timescales of days, months and years, and between 14-5100 MHz. Despite this, the shapes of the driftbands change radically with frequency. Previous studies have concluded that, while the subpulses appear to move through the pulse window approximately linearly at low frequencies (< 500 MHz), a discrete step of 180 degrees in subpulse phase is observed at higher frequencies (> 820 MHz) near to the peak of the average pulse profile. We use LOFAR, GMRT, GBT, WSRT and Effelsberg 100-m data to explore the frequency-dependence of this phase step. We show that the size of the subpulse phase step increases gradually, and is observable even at low frequencies. We attribute the subpulse phase step to the presence of two separate driftbands, whose relative arrival times vary with frequency - one driftband arriving 30 pulses earlier at 20 MHz than it does at 1380 MHz, whilst the other arrives simultaneously at all frequencies. The drifting pattern which is observed here cannot be explained by either the rotating carousel model or the surface oscillation model, and could provide new insight into the physical processes happening within the pulsar magnetosphere.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

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.

M87 is a giant elliptical galaxy located in the centre of the Virgo cluster, which harbours a supermassive black hole of mass 6.4x10^9 M_sun, whose activity is responsible for the extended (80 kpc) radio lobes that surround the galaxy. The energy generated by matter falling onto the central black hole is ejected and transferred to the intra-cluster medium via a relativistic jet and morphologically complex systems of buoyant bubbles, which rise towards the edges of the extended halo. Here we present the first observations made with the new Low-Frequency Array (LOFAR) of M87 at frequencies down to 20 MHz. Images of M87 were produced at low radio frequencies never explored before at these high spatial resolution and dynamic range. To disentangle different synchrotron models and place constraints on source magnetic field, age and energetics, we also performed a detailed spectral analysis of M87 extended radio-halo using these observations together with archival data. We do not find any sign of new extended emissions; on the contrary the source appears well confined by the high pressure of the intra-cluster medium. A continuous injection of relativistic electrons is the model that best fits our data, and provides a scenario in which the lobes are still supplied by fresh relativistic particles from the active galactic nuclei. We suggest that the discrepancy between the low-frequency radio-spectral slope in the core and in the halo implies a strong adiabatic expansion of the plasma as soon as it leaves the core area. The extended halo has an equipartition magnetic field strength of ~10 uG, which increases to ~13 uG in the zones where the particle flows are more active. The continuous injection model for synchrotron ageing provides an age for the halo of ~40 Myr, which in turn provides a jet kinetic power of 6-10x10^44 erg/s.

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.

Dispersion in the interstellar medium is a well known phenomenon that follows a simple relationship, which has been used to predict the time delay of dispersed radio pulses since the late 1960s. We performed wide-band simultaneous observations of four pulsars with LOFAR (at 40-190 MHz), the 76-m Lovell Telescope (at 1400 MHz) and the Effelsberg 100-m Telescope (at 8000 MHz) to test the accuracy of the dispersion law over a broad frequency range. In this paper we present the results of these observations which show that the dispersion law is accurate to better than 1 part in 100000 across our observing band. We use this fact to constrain some of the properties of the ISM along the line-of-sight and use the lack of any aberration or retardation effects to determine upper limits on emission heights in the pulsar magnetosphere. We also discuss the effect of pulse profile evolution on our observations, and the implications that it could have for precision pulsar timing projects such as the detection of gravitational waves with pulsar timing arrays.

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR Low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 \pm 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 \pm 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last \sim 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

We monitored the 22 GHz maser line in the lensed quasar MG J0414+0534 at z=2.64 with the 300-m Arecibo telescope for almost two years to detect possible additional maser components and to measure a potential velocity drift of the lines. The main maser line profile is complex and can be resolved into a number of broad features with line widths of 30-160 km/s. A new maser component was tentatively detected in October 2008 at a velocity of +470 km/s. After correcting for the estimated lens magnification, we find that the H2O isotropic luminosity of the maser in MG J0414+0534 is about 26,000 solar luminosities, making this source the most luminous ever discovered. Both the main line peak and continuum flux densities are surprisingly stable throughout the period of the observations. An upper limit on the velocity drift of the main peak of the line has been estimated from our observations and is of the order of 2 km/s per year. We discuss the results of the monitoring in terms of the possible nature of the maser emission, associated with an accretion disk or a radio jet. This is the first time that such a study is performed in a water maser source at high redshift, potentially allowing us to study the parsec-scale environment around a powerful radio source at cosmological distances.

When an ultra-high energy neutrino or cosmic ray strikes the Lunar surface a radio-frequency pulse is emitted. We plan to use the LOFAR radio telescope to detect these pulses. In this work we propose an efficient trigger implementation for LOFAR optimized for the observation of short radio pulses.

The Low Frequency Array (LOFAR) is under construction in the Netherlands and in several surrounding European countries. In this contribution, we describe the layout and design of the telescope, with a particular emphasis on the imaging characteristics of the array when used in its "standard imaging" mode. After briefly reviewing the calibration and imaging software used for LOFAR image processing, we show some recent results from the ongoing imaging commissioning efforts. We conclude by summarizing future prospects for the use of LOFAR in observing the little-explored low frequency Universe.

The bright submillimetre (sub-mm) galaxy MM 18423+5938 at redshift 3.9296 has been predicted from mid-infrared and millimetre photometry to have an exceptionally large total infrared (IR) luminosity. We present new radio imaging at 1.4 GHz with the Westerbork Synthesis Radio Telescope that is used to determine a radio-derived total IR luminosity for MM 18423+5938 via the well established radio-far-infrared correlation. The flux density is found to be S_1.4 GHz = 217 +/- 37 \mu Jy, which corresponds to a rest-frame luminosity density of L_1.4 GHz = 2.32 +/- 0.40 x 10^25 / u W / Hz, where u is the magnification from a probable gravitational lens. The radio-derived total IR luminosity and star-formation rate are L_8-1000 \mu m = 5.6^+4.1_-2.4 x 10^13 / u L_sol and SFR = 9.4^+7.4_-4.9 x 10^3 / u M_sol / yr, respectively, which are ~9 times smaller than those previously reported. These differences are attributed to the IR spectral energy distribution of MM 18423+5938 being poorly constrained by the limited number of reliable photometric data that are currently available, and from a previous misidentification of the object at 70 \mu m. Using the radio derived total IR luminosity as a constraint, the temperature of the cold dust component is found to be T ~ 24^+7_-5 K for a dust emissivity of \beta = 1.5 +/- 0.5. The radio-derived properties of this galaxy are still large given the low excitation temperature implied by the CO emission lines and the temperature of the cold dust. Therefore, we conclude that MM 18423+5938 is probably gravitationally lensed.

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 study of water masers at cosmological distances would allow us to investigate the parsec-scale environment around powerful radio sources, to probe the physical conditions of the molecular gas in the inner parsecs of quasars, and to estimate their nuclear engine masses in the early universe. To derive this information, the nature of the maser source, jet or disk-maser, needs to be assessed through a detailed investigation of the observational characteristics of the line emission. We monitored the maser line in the lensed quasar MGJ0414+0534 at z = 2.64 with the 300-m Arecibo telescope for ~15 months to detect possible additional maser components and to measure a potential velocity drift of the lines. In addition, we follow the maser and continuum emissions to reveal significant variations in their flux density and to determine correlation or time-lag, if any, between them. The main maser line profile is complex and can be resolved into a number of broad features with line widths of 30-160 km/s. A new maser component was tentatively detected in October 2008 that is redshifted by 470 km/s w.r.t the systemic velocity of the quasar. The line width of the main maser feature increased by a factor of two between the Effelsberg and EVLA observations reported by Impellizzeri et al. (2008) and the first epoch of the Arecibo monitoring campaign. After correcting for the lens magnification, we find that the total H2O isotropic luminosity of the maser in MGJ0414+0534 is now ~30,000 Lsun, making this source the most luminous ever discovered.[Abridged]

The gravitational lens system MG J0414+0534 is formed by an elliptical galaxy at redshift ~0.96 and a quasar at z~2.64. The system geometry is typical of lensing by an elliptical galaxy with the QSO close to and inside a fold caustic. It shows 4 images of the background source, and a partial Einstein ring is visible at optical wavelengths. It was observed with a global-VLBI array at 18 cm in June 2008. We present here the imaging results and a preliminary lens model constrained by these observations.

One of our closest neighbours, the Andromeda Galaxy (M31) has been the subject of numerous large area studies across the entire spectrum, but so far full-disk radio surveys have been conducted only at low resolution. The new wide-field capabilities of the DiFX software correlator present the possibility of imaging the entire primary beam of a VLBI array, thus enabling a high resolution wide-field study of the entire galaxy. Using the VLBA and EVN, pilot observations of M31 have been carried out with the aim of using these new wide-field techniques to characterise the population of compact components at VLBI resolution both within and behind one of our nearest neighbours. This contribution describes the observations carried out, the preliminary processing and first results.

The 8 o'clock arc is a gravitationally lensed Lyman Break Galaxy (LBG) at redshift z=2.73 that has a star-formation rate (SFR) of 270 solar-mass/year (derived from optical and near-infrared spectroscopy). Taking the magnification of the system ~12 and the SFR into account, the expected flux density of any associated radio emission at 1.4 GHz is predicted to be just 0.1 mJy. However, the lens system is found to be coincident with a radio source detected in the NRAO Very Large Array (VLA) Sky Survey with a flux density of ~5 mJy. If this flux density is attributed to the lensed LBG then it would imply a SFR ~11000 solar-mass/year, in contrast with the optical and near-infrared derived value. We want to investigate the radio properties of this system, and independently determine the SFR for the LBG from its lensed radio emission. We have carried out new high resolution imaging with the VLA ain A and B-configurations at 1.4 and 5 GHz. We find that the radio emission is dominated by a radio-loud AGN associated with the lensing galaxy. The radio-jet from the AGN partially covers the lensed arc of the LBG, and we do not detect any radio emission from the unobscured region of the arc down to a 3 sigma flux-density limit of 108 micro-Jy/beam. Using the radio data, we place a limit of <750 solar-mass/year for the SFR of the LBG, which is consistent with the results from the optical and near-infrared spectroscopy. We expect that the sensitivity of the Expanded VLA will be sufficient to detect many high redshift LBGs that are gravitationally lensed after only a few hours of observing time. The high angular resolution provided by the EVLA will also allow detailed studies of the lensed galaxies and determine if there is radio emission from the lens.

Luminous extragalactic water masers are known to be associated with AGN and have provided accurate estimates for the mass of the central supermassive black hole and the size and structure of the accretion disk in nearby galaxies. To find water masers at much higher redshifts, we have begun a survey of known gravitationally lensed quasars and star-forming galaxies. In this paper, we present a search for 22 GHz (rest frame) water masers toward five dusty, gravitationally lensed quasars and star-forming galaxies at redshifts 2.3--2.9 with the Effelsberg telescope and the EVLA. Our observations do not find any new definite examples of high redshift water maser galaxies, suggesting that large reservoirs of dust and gas are not a sufficient condition for powerful water maser emission. However, we do find the tentative detection of a water maser system in the active galaxy IRAS 10214+4724 at redshift 2.285. Our survey has now doubled the number of lensed galaxies and quasars that have been searched for high redshift water masers. We present an analysis of the high redshift water maser luminosity function that is based on the results presented here and from the only cosmologically distant (z > 1) water maser galaxy found thus far, MG J0414+0534 at redshift 2.64. By comparing with the luminosity function locally and at moderate redshifts, we find that there must be some evolution in the luminosity function of water maser galaxies at high redshifts. By assuming a moderate evolution [(1 + z )^4] in the luminosity function, we find that blind surveys for water maser galaxies are only worthwhile with extremely high sensitivity like that of the planned Square Kilometre Array. However, instruments like the EVLA and MeerKAT will be capable of detecting water maser systems similar to the one found from MG J0414+0534 through targeted observations.

I present first results of LOFAR observations with international baselines. An important cornerstone was the detection of the first long-baseline fringes. Their analysis turns out to be extremely useful to investigate and solve a number of technical issues of the instrument. Crude maps of the sky are created from single-baseline delay/fringe-rate spectra and compared with a short-baseline synthesis map. First long-baseline LBA images are shown of the source 3C196, a bright quasar whose sub-components can only be resolved with the long baselines. The corresponding sub-arcsec HBA image does not show the same amount of details yet, but HBA results are expected to improve significantly very soon. The LBA long-baseline image of 3C196 comprises the highest-resolution radio map ever produced at this low frequency.

SMM J04542-0301 is an extended (~1 arcmin) sub-mm source located near the core of the cluster MS0451.6-0305. It has been suggested that part of its emission arises from the interaction between a LBG and two EROs at z~2.9 that are multiply-imaged. However, the dramatic resolution difference between the sub-mm map and the optical/NIR images make it difficult to confirm this hypothesis. In this paper, we present a deep (~ 10 microJy/beam), high resolution (~2 arcsec) 1.4 GHz radio map of the cluster core, in which we have identified 6 sources located within SMM J04542-0301. The strong lensing effect in the radio data has been quantified by constructing a new lens model of the cluster. The brightest and most extended of these sources (RJ) is located in the middle of the sub-mm emission, and has no obvious counterpart in the optical/NIR. Three other detections (E1, E2 and E3) seem to be associated with the images of one of the EROs. The last two detections (CR1 and CR2), for which no optical/NIR counterpart have been found, seem to constitute two relatively compact emitting regions embedded in a ~5 arcsec extended radio source located at the position of the sub-mm peak. The presence of this extended component can only be explained if it is being produced by a lensed region of dust obscured star formation in the center of the merger. A comparison between the radio and sub-mm data at the same resolution suggests that E1, E2, E3, CR1 and CR2 are associated with the sub-mm emission. The radio observations provide strong observational evidence in favor of the merger hypothesis. However, the question if RJ is also contributing to the observed sub-mm emission remains open. These results illustrate the promising prospects for radio interferometry and strong gravitational lensing to study the internal structure of SMGs.

The Low Frequency Array (LOFAR) is a new generation of electronic radio telescope based on aperture array technology. The telescope is being developed by ASTRON, and currently being rolled out across the Netherlands and other countries in Europe. I present the current status of the project, and its relation to high resolution instruments such as the European VLBI Network (EVN) and the Very Long Baseline Array (VLBA). In particular, I present recent VLBI results at 327 MHz associated with: (i) a shallow survey based on VLBA archive data and (ii) a deep, wide-field Global VLBI survey centred on two in-beam calibrators, B0218+357 and J0226+3421. The results suggest that there will be no shortage of relatively bright primary calibrators that remain unresolved by LOFAR even on the longest European baselines. The sky density of fainter in-beam calibrators should also be more than adequate to permit the generation of high fidelity images over a large fraction of the sky, especially in the high-band observing band (120-240 MHz). Extending LOFAR via international stations to baseline lengths of several thousand kilometres is certainly practical and should significantly enhance the scientific output and capabilities of the array.

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

Water masers are found in dense molecular clouds closely associated with supermassive black holes in the centres of active galaxies. Based upon the understanding of the local water maser luminosity function, it was expected that masers at intermediate and high redshifts would be extremely rare, but galaxies at redshifts z > 2 might be quite different from those found locally, not least because of more frequent mergers and interaction events. Using gravitational lensing as a tool to enable us to search higher redshifts than would otherwise be possible, we have embarked on a survey of lensed galaxies, looking for masers. Here we report the discovery of a water maser at redshift 2.64 in the dust- and gas-rich gravitationally lensed type 1 quasar MG J0414+0534, which, with an isotropic luminosity of 10,000 L_solar, is twice as luminous as the most powerful local water maser, and half that of the most distant maser previously known. Using the locally-determined luminosity function, the probability of finding a maser this luminous associated with any single active galaxy is 10^{-6}. The fact that we saw such a maser in the first galaxy we observed must mean that the volume densities and luminosities of masers are higher at redshift 2.64.

So far the lens J1131-1231 has been studied only at optical and X-ray wavelengths. A detection in the radio was almost missed as a result of an incorrect position and archive problems. A direct analysis of NVSS uv data - in contrast to the catalogue or images alone - provided sufficient evidence of a detection to justify further radio investigations. The system was subsequently observed with MERLIN and the EVN in e-VLBI mode. Even though MERLIN seems to show the lensed star-forming regions and the compact cores, a preliminary analysis of the EVN data only shows an AGN in the lens itself but not the lensed cores. Additional VLA observations will be carried out soon.

Deep surveys planned as a Key Science Project of LOFAR provide completely new opportunities for gravitational lens searches. For the first time do large-scale surveys reach the resolution required for a direct selection of lens candidates using morphological criteria. We briefly describe the strategies that we will use to exploit this potential. The long baselines of an international E-LOFAR are essential for this project.

We discuss the classic theorem according to which a gravitational lens always produces a total magnification greater than unity. This theorem seems to contradict the conservation of total flux from a lensed source. The standard solution to this paradox is based on the exact definition of the reference 'unlensed' situation. We calculate magnifications and amplifications for general lensing scenarios not limited to regions close to the optical axis. In this way the formalism is naturally extended from tangential planes for the source and lensed images to complete spheres. We derive the lensing potential theory on the sphere and find that the Poisson equation is modified by an additional source term that is related to the mean density and to the Newtonian potential at the positions of observer and source. This new term generally reduces the magnification, to below unity far from the optical axis, and ensures conservation of the total photon number received on a sphere around the source. This discussion does not affect the validity of the 'focusing theorem', in which the unlensed situation is defined to have an unchanged affine distance between source and observer. The focusing theorem does not contradict flux conservation, because the mean total magnification directly corresponds to different areas of the source sphere in the lensed and unlensed situation. We argue that a constant affine distance does not define an astronomically meaningful reference. By exchanging source and observer, we confirm that magnification and amplification differ according to Etherington's reciprocity law, so that surface brightness is no longer strictly conserved. [ abridged ]

We report on the first wide-field, very long baseline interferometry (VLBI) survey at 90 cm. The survey area consists of two overlapping 28 deg^2 fields centred on the quasar J0226+3421 and the gravitational lens B0218+357. A total of 618 sources were targeted in these fields, based on identifications from Westerbork Northern Sky Survey (WENSS) data. Of these sources, 272 had flux densities that, if unresolved, would fall above the sensitivity limit of the VLBI observations. A total of 27 sources were detected as far as 2 arcdegrees from the phase centre. The results of the survey suggest that at least 10% of moderately faint (S~100 mJy) sources found at 90 cm contain compact components smaller than ~0.1 to 0.3 arcsec and stronger than 10% of their total flux densities. A ~90 mJy source was detected in the VLBI data that was not seen in the WENSS and NRAO VLA Sky Survey (NVSS) data and may be a transient or highly variable source that has been serendipitously detected. This survey is the first systematic (and non-biased), deep, high-resolution survey of the low-frequency radio sky. It is also the widest field of view VLBI survey with a single pointing to date, exceeding the total survey area of previous higher frequency surveys by two orders of magnitude. These initial results suggest that new low frequency telescopes, such as LOFAR, should detect many compact radio sources and that plans to extend these arrays to baselines of several thousand kilometres are warranted.

This is a short and biased review of gravitational lensing with emphasis on the radio and especially VLBI aspects. We briefly explain the basic idea and give a short sketch of the discovery of the first lens, before we more systematically discuss the general fields that can be studied with lensing. We intentionally omit the details the average lensing expert would like to see and instead try to give a very general overview addressed to the radio astronomer working in some other field. The lens B0218+357 is presented as an example to show many aspects of lensing in a case that already has led to interesting results but still has additional potential for the future.

We report on a recent 90 cm wide-field VLBI survey of two 3.1 deg^2 fields using the VLBA, Westerbork and Jodrell Bank telescopes. In-beam calibration was used to calibrate each field, the process was simplified by imaging the calibrators in DIFMAP and transferring the calibration solutions to AIPS using the newly developed DIFMAP task - cordump. We detected and imaged 13 out of the 141 sources originally detected by the low resolution (54") WENSS survey of the same two fields. The sources were detected at 7-12 sigma levels above the image noise, had total flux densities ranging between 85-1640 mJy and were between 16'-58' from the phase centre of each field. This is the first systematic (and non-biased), deep, high resolution survey of the low frequency radio sky. These initial results suggest that new instruments such as LOFAR should detect many compact radio sources and that plans to extend these arrays to baselines of several thousand kilometres are warranted.

Previous authors have reported the detection of intrinsically faint sub-mm emission lensed by the cluster MS0451.6-0305. They suggest that this emission arises from a merging system composed of a Ly-break galaxy and a pair of extremely red objects which are multiply-imaged in the optical/NIR observations. Since the submm emission presents an unusually large angular extent (~1 arcmin), the possible radio emission asociatted with that system can help to identify optical/NIR counterparts due to the higher spatial resolution and astrometric accuracy of the radio observations. Archive VLA data (BnA configuration at 1.4 GHz) was reduced and analysed. A simple lens model was constructed to aid the interpretation of the radio and pre-existing sub-mm and optical/NIR data. We present a 1.4 GHz map of the central region of MS0451.6-0305 and report the detection of gravitationally lensed radio emission, coincident with the previously discovered sub-mm lensed emission. The overall morphology and scale of the radio and sub-mm emission are strikingly similar, extending ~1 arcmin across the sky. This observation strongly suggests that the radio and sub-mm emission arise from the same sources. Preliminary estimates of the total S_850 microns/S_1.4 GHz flux density ratio appear to be consistent with that expected from distant star forming galaxies. The radio emission is resolved into 7 distinct components, and the overall structure can be explained, using a simple lens model, with three multiply-imaged radio sources at z~2.9. One of these sources is predicted to lie in the middle of the previously mentioned system in the source plane, suggesting that it is related to the intense star formation generated during the merging process.

We address the issue of anomalous image flux ratios seen in the double-image gravitational lens JVAS B0218+357. From the multi-frequency observations presented in a recent study (Mittal et al. 2006) and several previous observations made by other authors, the anomaly is well-established in that the image flux-density ratio (A/B) decreases from 3.9 to 2.0 over the observed frequency range from 15 GHz to 1.65 GHz. In Mittal et al. (2006), the authors investigated whether an interplay between a frequency-dependent structure of the background radio-source and a gradient in the relative image-magnification can explain away the anomaly. Insufficient shifts in the image centroids with frequency led them to discard the above effect as the cause of the anomaly. In this paper, we first take this analysis further by evaluating the combined effect of the background source extension and magnification gradients in the lens plane in more detail. This is done by making a direct use of the observed VLBI flux-distributions for each image to estimate the image flux-density ratios at different frequencies from a lens-model. As a result of this investigation, this mechanism does not account for the anomaly. Following this, we analyze the effects of mechanisms which are non-gravitational in nature on the image flux ratios in B0218+357. These are free-free absorption and scattering, and are assumed to occur under the hypothesis of a molecular cloud residing in the lens galaxy along the line-of-sight to image A. We show that free-free absorption due to an H II region covering the entire structure of image A at 1.65 GHz can explain the image flux ratio anomaly. We also discuss whether H II regions with physical parameters as derived from our analysis are consistent with those observed in Galactic and extragalactic H II regions.

We present the results of phase-referenced VLBA+Effelsberg observations at five frequencies of the double-image gravitational lens JVAS B0218+357, made to establish the precise registration of the A and B lensed image positions. The motivation behind these observations is to investigate the anomalous variation of the image flux density ratio (A/B) with frequency - this ratio changes by almost a factor of two over a frequency range from 1.65 GHz to 15.35 GHz. We investigate whether frequency dependent image positions, combined with a magnification gradient across the image field, could give rise to the anomaly. Our observations confirm the variation of image flux ratio with frequency. The results from our phase-reference astrometry, taken together with the lens mass model of Wucknitz et al. (2004), show that shifts of the image peaks and centroids are too small to account for the observed frequency-dependent ratio.

We present deep optical observations of the gravitational lens system CLASS B0218+357, from which we derive an estimate of the Hubble constant. Extensive radio observations have reduced the degeneracies between Hubble's constant and the mass model in this lens to one involving only the position of the radio-quiet lensing galaxy relative to the lensed images. B0218+357 has an image separation of only 334 mas, so optical observations have previously been unable to resolve the lens galaxy from the bright lensed images. Using the new Advanced Camera for Surveys installed on the Hubble Space Telescope, we have measured the separation between the lens galaxy centre and the brightest image. The position found, and hence our estimate of Hubble's constant, depends on our approach to the spiral arms in B0218+357. If the most prominent arms are left unmasked, we find a value for Hubble's constant of 70+/-5 km/s /Mpc (95% confidence). If the spiral arms are masked out, we find a value of 61+/-7 km/s /Mpc (95% confidence).

We study metal abundances in the z=0.9313 damped Lya system observed in the two lines-of-sight, A and B, toward the gravitationally-lensed double QSO HE0512-3329. Spatially resolved STIS spectra constrain the neutral-gas column density to be LogN(HI)=20.5 in both Aand B. UVES spectra (spectral resolution FWHM=9.8 km/s) show, in contrast, significant line-of-sight differences in the column densities of MnII and FeII; these are not due to observational systematics. We find that [Mn/H]=-1.44 and [Fe/H]=-1.52 in damped Lya system A, while [Mn/H]=-0.98 and [Fe/H]>-1.32, and possibly as high as [Fe/H] approx. -1 in damped Lya system B. A careful assessment of possible systematic errors leads us to conclude that these transverse differences are significant at a 5 sigma level or greater. Although nucleosynthesis effects may also be at play, we favor differential dust-depletion as the main mechanism producing the observed abundance gradient. The transverse separation is 5 kpc at the redshift of the absorber, which is also likely to be the lensing galaxy. The derived abundances therefore probe two opposite sides of a single galaxy hosting both damped Lya systems. This is the first time firm abundance constraints have been obtained for a single damped system probed by two lines-of-sight. The significance of this finding for the cosmic evolution of metals is discussed.

We present the results of phase-referenced VLBA+Effelsberg observations at five frequencies of the gravitational lens B0218+357 to establish the precise registration of the A and B lensed image positions.

We discuss the Sagnac effect in standard Minkowski coordinates and with an alternative synchronization convention. We find that both approaches lead to the same result without any contradictions. When applying standard coordinates to the complete rim of the rotating disk, a time-lag has to be taken into account which accounts for the global anisotropy. We propose a closed Minkowski space-time as an exact equivalent to the rim of a disk, both in the rotating and non-rotating case. In this way the Sagnac effect can be explained as being purely topological, neglecting the radial acceleration altogether. This proves that the rim of the disk can be treated as an inertial system. In the same context we discuss the twin paradox and find that the standard scenario is equivalent to an unaccelerated version in a closed space-time. The closed topology leads to preferred frame effects which can be detected only globally. The relation of synchronization conventions to the measurement of lengths is discussed in the context of Ehrenfest's paradox. This leads to a confirmation of the classical arguments by Ehrenfest and Einstein.

Various authors have investigated the problem of light deflection by radially moving gravitational lenses - with inconsistent results. In this paper we generalize the calculations for arbitrary lens velocities and show that, to first order in the lens velocity v, the deflection angles scales with 1-v. We discuss the seeming inconsistency of relativistic light deflection with the classical picture of moving test particles by generalizing the lens effect to test particles of arbitrary velocity, including light as a limiting case. We show that the effect of radial motion of the lens is very different for slowly moving test particles and light and that a critical test particle velocity exists for which the motion of the lens has no effect on the deflection angle to first order. An interesting and not immediately intuitive result is obtained in the limit of a highly relativistic motion of the lens towards the observer, where the deflection angle of light reduces to zero. This phenomenon is elucidated in terms of moving refractive media. Furthermore, we discuss the dragging of inertial frames in the field of a moving lens and the corresponding Lense-Thirring precession, in order to shed more light on the geometrical effects in the surroundings of a moving mass. In a second part we discuss the effect of transversal motion on the observed redshift of lensed sources. We demonstrate how a simple kinematic calculation explains the effects for arbitrary velocities of the lens and test particles. Additionally we include the transversal motion of the source and observer to show that all three velocities can be combined into an effective relative transversal velocity similar to the approach used in microlensing studies.

We discuss the LensClean algorithm which for a given gravitational lens model fits a source brightness distribution to interferometric radio data in a similar way as standard Clean does in the unlensed case. The lens model parameters can then be varied in order to minimize the residuals and determine the best model for the lens mass distribution. Our variant of this method is improved in order to be useful and stable even for high dynamic range systems with nearly degenerated lens model parameters. Our test case B0218+357 is dominated by two bright images but the information needed to constrain the unknown parameters is provided only by the relatively smooth and weak Einstein ring. The new variant of LensClean is able to fit lens models even in this difficult case. In order to allow the use of general mass models with LensClean, we develop the new method LenTil which inverts the lens equation much more reliably than any other method. This high reliability is essential for the use as part of LensClean. Finally a new method is developed to produce source plane maps of the unlensed source from the best LensClean brightness models. This method is based on the new concept of dirty beams in the source plane. The application to the lens B0218+357 leads to the first useful constraints for the lens position and thus to a result for the Hubble constant. These results are presented in an accompanying Paper II, together with a discussion of classical lens modelling for this system.

B0218+357 is one of the most promising systems to determine the Hubble constant from gravitational lenses. Consisting of two bright resolved images plus an Einstein ring, it provides better constraints for the mass model than other systems. The main problem left until now was the poorly determined position of the lensing galaxy. After presenting results from classical lens modelling, we apply our improved version of LensClean which utilizes the Einstein ring for lens modelling purposes. The primary result using isothermal models is a well defined lens position which allows the first reliable measurement of the Hubble constant from this system. The result of H0=(78+-6) km/s/Mpc (2 sigma) is high compared with other lenses but compatible with the HST key project and WMAP results. We furthermore discuss effects of different radial mass profiles. The power-law exponent of the potential is constrained by VLBI data to be beta=1.04+-0.02, very close to isothermal. The effect on H0 is expected to be very small. We also present a composite map (lensed and unlensed) which shows the rich structure of B0218+357 on scales from milli-arcseconds to arcseconds. Finally we use a comparison of observations at different frequencies to investigate the question of possible weakening of one of the images by propagation effects and/or source shifts with frequency. The data clearly favour the model of significant extinction without noticeable source position shifts. The technical details of our variant of the LensClean method are presented in the accompanying Paper I.

We present the first separate spectra of both components of the small-separation double QSO HE0512-3329 obtained with HST/STIS in the optical and near UV. The similarities especially of the emission line profiles and redshifts strongly suggest that this system really consists of two lensed images of one and the same source. The emission line flux ratios are assumed to be unaffected by microlensing and are used to study the differential extinction effects caused by the lensing galaxy. Fits of empirical laws show that the extinction properties seem to be different on both lines of sight. With our new results, HE0512-3329 becomes one of the few extragalactic systems which show the 2175 AA absorption feature, although the detection is only marginal. We then correct the continuum flux ratio for extinction to obtain the differential microlensing signal. Since this may still be significantly affected by variability and time-delay effects, no detailled analysis of the microlensing is possible at the moment. This is the first time that differential extinction and microlensing could be separated unambiguously. We show that, at least in HE0512-3329, both effects contribute significantly to the spectral differences and one cannot be analysed without taking into account the other. For lens modelling purposes, the flux ratios can only be used after correcting for both effects.

The unification of X-ray and radio selected BL Lacs has been an outstanding problem in the blazar research in the past years. Recent investigations have shown that the gap between the two classes can be filled with intermediate objects and that apparently all differences can be explained by mutual shifts of the peak frequencies of the synchrotron and inverse Compton component of the emission. We study the consequences of this scheme using a new sample of X-ray selected BL Lac objects comprising 104 objects with z<0.9 and a mean redshift z=0.34. 77 BL Lacs, of which the redshift could be determined for 64 (83%) objects, form a complete sample. The new data could not confirm our earlier result, drawn from a subsample, that the negative evolution vanishes below a synchrotron peak frequency log (peak-frequency) = 16.5. The complete sample shows negative evolution at the 2 sigma level (<V_e/V_a> = 0.42 +- 0.04). We conclude that the observed properties of the HRX BL Lac sample show typical behaviour for X-ray selected BL Lacs. They support an evolutionary model, in which flat-spectrum radio quasars (FSRQ) with high energetic jets evolve towards low frequency peaked (mostly radio-selected) BL Lac objects and later on to high frequency peaked (mostly X-ray selected) BL Lacs.

In this paper we present new observations of the gravitational lens system JVAS B0218+357 made with a global VLBI network at a frequency of 8.4 GHz. Our maps have an rms noise of 30 microJy/beam and with these we have been able to image much of the extended structure of the radio jet in both the A and B images at high resolution (~1 mas). The main use of these maps will be to enable us to further constrain the lens model for the purposes of H0 determination. We are able to identify several sub-components common to both images with the expected parity reversal, including one which we identify as a counter-jet. We have not been successful in detecting either the core of the lensing galaxy or a third image. Using a model of the lensing galaxy we have back-projected both of the images to the source plane and find that they agree well. However, there are small, but significant, differences which we suggest may arise from multi-path scattering in the ISM of the lensing galaxy. We also find an exponent of the radial mass distribution of approximately 1.04, in agreement with lens modelling of published 15-GHz VLBI data. Polarisation maps of each image are presented which show that the distributions of polarisation across images A and B are different. We suggest that this results from Faraday rotation and associated depolarisation in the lensing galaxy.

The time delay in gravitational lenses can be used to derive the Hubble constant in a relatively simple way. The results of this method are less dependent on astrophysical assumptions than in many other methods. The most important uncertainty is related to the mass model used. We discuss a family of models with a separable radial power-law and an arbitrary angular dependence for the potential psi = r^beta * F(theta). Isothermal potentials are a special case of these models with beta=1. An additional external shear is used to take into account perturbations from other galaxies. Using a simple linear formalism for quadruple lenses, we can derive H0 as a function of the observables and the shear. If the latter is fixed, the result depends on the assumed power-law exponent according to H0 proportional to (2-beta)/beta. The effect of external shear is quantified by introducing a `critical shear' gamma_c as a measure for the amount of shear that changes the result significantly. The analysis shows, that in the general case H0 and gamma_c do not depend on the position of the lens galaxy. We discuss these results and compare with numerical models for a number of real lens systems.

We use the radio ring in JVAS B0218+357 to constrain the mass models with LensCLEAN. H0 is determined from the resulting model and the time delay. B0218+357 is one of the best systems to reduce systematic errors in H0 and may thus be called a "golden lens".

We present optical spectroscopy data that allowed a measurement of the redshift for the X-ray selected BL Lacertae object 1517+656. With a redshift of z = 0.702 this object has an absolute magnitude of B = -26.4 and is also an extremely powerful radio and X-ray source. Although being a high frequency peaked BL Lac, this object is one of the most luminous BL Lac objects known so far. Being also a candidate for gravitational lensing, this object is of high interest for the BL Lac research. Assuming several cosmological models and a realistic redshift for the lensed object, we find that 1517+656 has a mass > 2 * 10^12 solar masses and a high velocity dispersion > 350 km/sec.

We discuss a simple parameter degeneracy in lens models which follows directly from the well known mass sheet degeneracy and prevents determining the radial mass distribution from lensing alone for many systems. For cosmological important systems like 1115+080 it is shown that no exact value of H0 can be determined without assuming a mass index beta.

We present first results from five years of spectrophotometric monitoring of the bright double QSO and gravitational lens HE 1104-1805. The quasar has varied considerably over this time, while the emission line fluxes appear to have remained constant. We have constructed monochromatic continuum light curves for components A and B, finding that B leads the variability. A quantitative analysis with the Pelt method gives a best estimate for the light travel time delay of about 0.73 years, although a value as low as 0.3 cannot yet be excluded. We discuss possible models for the QSO-lens configuration and use our measured time delay to predict the redshift of the lens, z_d. Finding that most likely z_d < 1, we can rule out the hitherto favoured values of z_d = 1.32 or 1.66. A new candidate is an absorption system at z=0.73, but the lens could also be an elliptical not detected in absorption.

We have obtained new BVR-band images of the gravitationally lensed QSO pair HE 2149-2745 A,B (QSO redshift z=2.033; angular separation theta_AB=1.71 arcsec), under 1 arcsec seeing conditions. Subtraction of a scaled point-spread-function reveals an extended object of angular size comparable to theta_AB, and R magnitude of 20.4 lying between and equidistant to the QSOs. The geometry of the system leads us to interpret this object as the lensing agent responsible for the double QSO images. A fit to the galaxy surface brightness profile indicates an elliptical galaxy, which is further supported by the absence of Mg II in absorption. The lensing galaxy is not detected in B and V. We discuss possible lens models for this geometry constrained by our observations and consequences for the mass and redshift of the lens.

We report the detection of four images in the recently discovered lensed QSO RX J0911.4+0551. With a maximum angular separation of 3.1", it is the quadruply imaged QSO with the widest known angular separation. Raw and deconvolved data reveal an elongated lens galaxy. The observed reddening in at least two of the four QSO images suggests differential extinction by this lensing galaxy. We show that both an ellipticity of the galaxy (epsilon_{min}=0.075) and an external shear (gamma_{min}=0.15) from a nearby mass has to be included in the lensing potential in order to reproduce the complex geometry observed in RX J0911.4+0551. A possible galaxy cluster is detected about 38", from RX J0911.4+0551 and could contribute to the X-ray emission observed by ROSAT in this field. The color of these galaxies indicates a plausible redshift in the range of 0.6-0.8.