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.

We present results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a $30\,{\rm deg}^2$ field at $\Delta T$=2, 9 and 33 days post-merger at a frequency of 944\,MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers $89\%$ of the LIGO/Virgo localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from neutron star-black hole mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

Active M-dwarfs are known to produce bursty radio emission, and multi-wavelength studies have shown that Solar-like magnetic activity occurs in these stars. However, coherent bursts from active M-dwarfs have often been difficult to interpret in the Solar activity paradigm. We present Australian Square Array Pathfinder (ASKAP) observations of UV Ceti at a central frequency of 888 MHz. We detect several periodic, coherent pulses occurring over a timescale consistent with the rotational period of UV Ceti. The properties of the pulsed emission show that they originate from the electron cyclotron maser instability, in a cavity $\sim7$ orders of magnitude less dense than the mean coronal density at the estimated source altitude. These results confirm that auroral activity can occur in active M-dwarfs, suggesting that these stars mark the beginning of the transition from Solar-like to auroral magnetospheric behaviour. These results demonstrate the capabilities of ASKAP for detecting polarized, coherent bursts from active stars and other systems.

With the arrival of a number of wide-field snapshot image-plane radio transient surveys, there will be a huge influx of images in the coming years making it impossible to manually analyse the datasets. Automated pipelines to process the information stored in the images are being developed, such as the LOFAR Transients Pipeline, outputting light curves and various transient parameters. These pipelines have a number of tuneable parameters that require training to meet the survey requirements. This paper utilises both observed and simulated datasets to demonstrate different machine learning strategies that can be used to train these parameters. The datasets used are from LOFAR observations and we process the data using the LOFAR Transients Pipeline; however the strategies developed are applicable to any light curve datasets at different frequencies and can be adapted to different automated pipelines. These machine learning strategies are publicly available as Python tools that can be downloaded and adapted to different datasets (https://github.com/AntoniaR/TraP_ML_tools).

With the advent of powerful telescopes such as the SKA and LSST, we are entering a golden era of multiwavelength transient astronomy. In order to cope with the dramatic increase in data volume as well as successfully prioritise spectroscopic follow-up resources, we propose a new machine learning approach for the classification of radio and multiwavelength transients. The algorithm consists of three steps: (1) augmentation and interpolation of the data using Gaussian processes; (2) feature extraction using a wavelet decomposition; (3) classification with the robust machine learning algorithm known as random forests. We apply this algorithm to existing radio transient data, illustrating its ability to accurately classify most of the eleven classes of radio variables and transients after just eight hours of observations, achieving an overall accuracy of $73.5\%$. We show how performance is expected to increase as more training data are acquired, by training the classifier on a simulated representative training set, achieving an overall accuracy of $97.4\%$. Finally, we outline a general approach for including multiwavelength data for general transient classification, and demonstrate its effectiveness by incorporating a single optical data point into the analysis, which improves the overall accuracy by $\approx 22\%$.

The shape of low-frequency radio continuum spectra of normal galaxies is not well understood, the key question being the role of physical processes such as thermal absorption in shaping them. In this work we take advantage of the LOFAR Multifrequency Snapshot Sky Survey (MSSS) to investigate such spectra for a large sample of nearby star-forming galaxies. Using the measured 150MHz flux densities from the LOFAR MSSS survey and literature flux densities at various frequencies we have obtained integrated radio spectra for 106 galaxies. The spectra are explained through the use of a three-dimensional model of galaxy radio emission, and radiation transfer dependent on the galaxy viewing angle and absorption processes. Spectra of our galaxies are generally flatter at lower compared to higher frequencies but as there is no tendency for the highly inclined galaxies to have more flattened low-frequency spectra, we argue that the observed flattening is not due to thermal absorption, contradicting the suggestion of Israel & Mahoney (1990). According to our modelled radio maps for M51-like galaxies, the free-free absorption effects can be seen only below 30MHz and in the global spectra just below 20MHz, while in the spectra of starburst galaxies, like M82, the flattening due to absorption is instead visible up to higher frequencies of about 150MHz. Locally, within galactic disks, the absorption effects are distinctly visible in M51-like galaxies as spectral flattening around 100-200MHz in the face-on objects, and as turnovers in the edge-on ones, while in M82-like galaxies there are strong turnovers at frequencies above 700MHz, regardless of viewing angle. Our modelling suggests that the weak spectral flattening observed in the nearby galaxies studied here results principally from synchrotron spectral curvature due to cosmic ray energy losses and propagation effects.

We investigate the relation between the radio (F_r) and optical (F_o) flux densities of a variety of classes of radio transients and variables, with the aim of analysing whether this information can be used, in the future, to classify such events. Using flux density values between 1-10 GHz and the optical bands V and R, we build a sample with a total of 12,441 F_r and F_o measurements. The sample contains both Galactic objects, such as stellar sources and X-ray binaries, and extragalactic objects, such as gamma-ray bursts and quasars. By directly comparing the two parameters, it is already possible to distinguish between the Galactic and extragalactic populations. Although individual classes are harder to separate from the F_r-F_o parameter space to a high accuracy, and can only provide approximations, the basic approach provides an already useful foundation to develop a more accurate classification technique. In addition, we illustrate how example objects from different classes move in the parameter space as they evolve over time, offering a feature that could be used to reduce the confusion between classes. A small, blind test of the classification performance is also undertaken using a catalogue of FIRST transient and variable sources, to demonstrate the advantages and current limitations of the approach. With more multi-wavelength data becoming available in the future, we discuss other classification techniques which the F_r-F_o method could be combined with and potentially become an important part of an automatic radio transient classification system.

With their wide fields of view and often relatively long coverage of any position in the sky in imaging survey mode, modern radio telescopes provide a data stream that is naturally suited to searching for rare transients. However, Radio Frequency Interference (RFI) can show up in the data stream in similar ways to such transients, and thus the normal pre-treatment of filtering RFI (flagging) may also remove astrophysical transients from the data stream before imaging. In this paper we investigate how standard flagging affects the detectability of such transients by examining the case of transient detection in an observing mode used for Low Frequency Array (LOFAR; \citep{LOFAR}) surveys. We quantify the fluence range of transients that would be detected, and the reduction of their SNR due to partial flagging. We find that transients with a duration close to the integration sampling time, as well as bright transients with durations on the order of tens of seconds, are completely flagged. For longer transients on the order of several tens of seconds to minutes, the flagging effects are not as severe, although part of the signal is lost. For these transients, we present a modified flagging strategy which mitigates the effect of flagging on transient signals. We also present a script which uses the differences between the two strategies, and known differences between transient RFI and astrophysical transients, to notify the observer when a potential transient is in the data stream.

We present LOFAR high-band data over the frequency range 115-189 MHz for the X-ray binary SS 433, obtained in an observing campaign from 2013 February - 2014 May. Our results include a deep, wide-field map, allowing a detailed view of the surrounding supernova remnant W 50 at low radio frequencies, as well as a light curve for SS 433 determined from shorter monitoring runs. The complex morphology of W 50 is in excellent agreement with previously published higher-frequency maps; we find additional evidence for a spectral turnover in the eastern wing, potentially due to foreground free-free absorption. Furthermore, SS 433 is tentatively variable at 150 MHz, with both a debiased modulation index of 11 per cent and a $\chi^2$ probability of a flat light curve of $8.2 \times 10^{-3}$. By comparing the LOFAR flux densities with contemporaneous observations carried out at 4800 MHz with the RATAN-600 telescope, we suggest that an observed $\sim$0.5-1 Jy rise in the 150-MHz flux density may correspond to sustained flaring activity over a period of approximately six months at 4800 MHz. However, the increase is too large to be explained with a standard synchrotron bubble model. We also detect a wealth of structure along the nearby Galactic plane, including the most complete detection to date of the radio shell of the candidate supernova remnant G 38.7-1.4. This further demonstrates the potential of supernova remnant studies with the current generation of low-frequency radio telescopes.

We present a deep, low-frequency radio continuum study of the nearby Fanaroff--Riley class I (FR I) radio galaxy 3C 31 using a combination of LOw Frequency ARray (LOFAR; 30--85 and 115--178 MHz), Very Large Array (VLA; 290--420 MHz), Westerbork Synthesis Radio Telescope (WSRT; 609 MHz) and Giant Metre Radio Telescope (GMRT; 615 MHz) observations. Our new LOFAR 145-MHz map shows that 3C 31 has a largest physical size of $1.1$ Mpc in projection, which means 3C 31 now falls in the class of giant radio galaxies. We model the radio continuum intensities with advective cosmic-ray transport, evolving the cosmic-ray electron population and magnetic field strength in the tails as functions of distance to the nucleus. We find that if there is no in-situ particle acceleration in the tails, then decelerating flows are required that depend on radius $r$ as $v\propto r^{\beta}$ ($\beta\approx -1$). This then compensates for the strong adiabatic losses due to the lateral expansion of the tails. We are able to find self-consistent solutions in agreement with the entrainment model of Croston & Hardcastle, where the magnetic field provides $\approx$$1/3$ of the pressure needed for equilibrium with the surrounding intra-cluster medium (ICM). We obtain an advective time-scale of $\approx$$190$ Myr, which, if equated to the source age, would require an average expansion Mach number ${\cal M} \approx 5$ over the source lifetime. Dynamical arguments suggest that instead, either the outer tail material does not represent the oldest jet plasma or else the particle ages are underestimated due to the effects of particle acceleration on large scales.

We present a new analysis of the widely used relation between cavity power and radio luminosity in clusters of galaxies with evidence for strong AGN feedback. We study the correlation at low radio frequencies using two new surveys - the First Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) at 148 MHz and LOFAR's first all-sky survey, the Multifrequency Snapshot Sky Survey (MSSS) at 140 MHz. We find a scaling relation $P_{\rm cav} \propto L_{148}^{\beta}$, with a logarithmic slope of $\beta = 0.51 \pm 0.14$, which is in good agreement with previous results based on data at 327 MHz. The large scatter present in this correlation confirms the conclusion reached at higher frequencies that the total radio luminosity at a single frequency is a poor predictor of the total jet power. We show that including measurements at 148 MHz alone is insufficient to reliably compute the bolometric radio luminosity and reduce the scatter in the correlation. For a subset of four well-resolved sources, we examine the detected extended structures at low frequencies and compare with the morphology known from higher frequency images and Chandra X-ray maps. In Perseus we discuss details in the structures of the radio mini-halo, while in the 2A 0335+096 cluster we observe new diffuse emission associated with multiple X-ray cavities and likely originating from past activity. For A2199 and MS 0735.6+7421, we confirm that the observed low-frequency radio lobes are confined to the extents known from higher frequencies. This new low-frequency analysis highlights the fact that existing cavity power to radio luminosity relations are based on a relatively narrow range of AGN outburst ages. We discuss how the correlation could be extended using low frequency data from the LOFAR Two-metre Sky Survey (LoTSS) in combination with future, complementary deeper X-ray observations.

We report on the discovery in the LOFAR Multifrequency Snapshot Sky Survey (MSSS) of a giant radio galaxy (GRG) with a projected size of $2.56 \pm 0.07$ Mpc projected on the sky. It is associated with the galaxy triplet UGC 9555, within which one is identified as a broad-line galaxy in the Sloan Digital Sky Survey (SDSS) at a redshift of $0.05453 \pm 1 \times 10^{-5} $, and with a velocity dispersion of $215.86 \pm 6.34$ km/s. From archival radio observations we see that this galaxy hosts a compact flat-spectrum radio source, and we conclude that it is the active galactic nucleus (AGN) responsible for generating the radio lobes. The radio luminosity distribution of the jets, and the broad-line classification of the host AGN, indicate this GRG is orientated well out of the plane of the sky, making its physical size one of the largest known for any GRG. Analysis of the infrared data suggests that the host is a lenticular type galaxy with a large stellar mass ($\log~\mathrm{M}/\mathrm{M}_\odot = 11.56 \pm 0.12$), and a moderate star formation rate ($1.2 \pm 0.3~\mathrm{M}_\odot/\mathrm{year}$). Spatially smoothing the SDSS images shows the system around UGC 9555 to be significantly disturbed, with a prominent extension to the south-east. Overall, the evidence suggests this host galaxy has undergone one or more recent moderate merger events and is also experiencing tidal interactions with surrounding galaxies, which have caused the star formation and provided the supply of gas to trigger and fuel the Mpc-scale radio lobes.

Radio observations of young stellar objects (YSOs) enable the study of ionised plasma outflows from young protostars via their free-free radiation. Previous studies of the low-mass young system T Tau have used radio observations to model the spectrum and estimate important physical properties of the associated ionised plasma (local electron density, ionised gas content and emission measure). However, without an indication of the low-frequency turnover in the free-free spectrum, these properties remain difficult to constrain. This paper presents the detection of T Tau at 149 MHz with the Low Frequency Array (LOFAR) - the first time a YSO has been observed at such low frequencies. The recovered total flux indicates that the free-free spectrum may be turning over near 149 MHz. The spectral energy distribution is fitted and yields improved constraints on local electron density ($(7.2 \pm 2.1)\times10^{3}$ cm$^{-3}$), ionised gas mass ($(1.0 \pm 1.8)\times10^{-6}$ M$_{\odot}$) and emission measure ($(1.67 \pm 0.14)\times10^5$ pc cm$^{-6}$).

This is the third installment in a series of papers in which we investigate calibration artefacts. Calibration artefacts (also known as ghosts or spurious sources) are created when we calibrate with an incomplete model. In the first two papers of this series we developed a mathematical framework which enabled us to study the ghosting mechanism itself. An interesting concomitant of the second paper was that ghosts appear in symmetrical pairs. This could possibly account for spurious symmetrization. Spurious symmetrization refers to the appearance of a spurious source (the anti-ghost) symmetrically opposite an unmodelled source around a modelled source. The analysis in the first two papers indicates that the anti-ghost is usually very faint, in particular when a large number of antennas are used. This suggests that spurious symmetrization will mainly occur at an almost undetectable flux level. In this paper, we show that phase-only calibration produces an anti-ghost that is $N$-times (where $N$ denotes the number of antennas in the array) as bright as the one produced by phase and amplitude calibration and that this already bright ghost can be further amplified by the primary beam correction.

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.

The eclipses of certain types of binary millisecond pulsars (i.e. `black widows' and `redbacks') are often studied using high-time-resolution, `beamformed' radio observations. However, they may also be detected in images generated from interferometric data. As part of a larger imaging project to characterize the variable and transient sky at radio frequencies <200 MHz, we have blindly detected the redback system PSR J2215+5135 as a variable source of interest with the Low-Frequency Array (LOFAR). Using observations with cadences of 2 weeks - 6 months, we find preliminary evidence that the eclipse duration is frequency dependent ($\propto \nu^{-0.4}$), such that the pulsar is eclipsed for longer at lower frequencies, in broad agreement with beamformed studies of other similar sources. Furthermore, the detection of the eclipses in imaging data suggests an eclipsing medium that absorbs the pulsed emission, rather than scattering it. Our study is also a demonstration of the prospects of finding pulsars in wide-field imaging surveys with the current generation of low-frequency radio telescopes.

We report on the results of a search for radio transients between 115 and 190\,MHz with the LOw-Frequency ARray (LOFAR). Four fields have been monitored with cadences between 15 minutes and several months. A total of 151 images were obtained, giving a total survey area of 2275 deg$^2$. We analysed our data using standard LOFAR tools and searched for radio transients using the LOFAR Transients Pipeline (TraP). No credible radio transient candidate has been detected; however, we are able to set upper limits on the surface density of radio transient sources at low radio frequencies. We also show that low-frequency radio surveys are more sensitive to steep-spectrum coherent transient sources than GHz radio surveys. We used two new statistical methods to determine the upper limits on the transient surface density. One is free of assumptions on the flux distribution of the sources, while the other assumes a power-law distribution in flux and sets more stringent constraints on the transient surface density. Both of these methods provide better constraints than the approach used in previous works. The best value for the upper limit we can set for the transient surface density, using the method assuming a power-law flux distribution, is 1.3$\cdot$10$^{-3}$ deg$^{-2}$ for transients brighter than 0.3 Jy with a time-scale of 15 min, at a frequency of 150 MHz. We also calculated for the first time upper limits for the transient surface density for transients of different time-scales. We find that the results can differ by orders of magnitude from previously reported, simplified estimates.

Due to their steep spectra, low-frequency observations of FR II radio galaxies potentially provide key insights in to the morphology, energetics and spectrum of these powerful radio sources. However, limitations imposed by the previous generation of radio interferometers at metre wavelengths has meant that this region of parameter space remains largely unexplored. In this paper, the first in a series examining FR IIs at low frequencies, we use LOFAR observations between 50 and 160 MHz, along with complementary archival radio and X-ray data, to explore the properties of two FR II sources, 3C452 and 3C223. We find that the morphology of 3C452 is that of a standard FR II rather than of a double-double radio galaxy as had previously been suggested, with no remnant emission being observed beyond the active lobes. We find that the low-frequency integrated spectra of both sources are much steeper than expected based on traditional assumptions and, using synchrotron/inverse-Compton model fitting, show that the total energy content of the lobes is greater than previous estimates by a factor of around 5 for 3C452 and 2 for 3C223. We go on to discuss possible causes of these steeper than expected spectra and provide revised estimates of the internal pressures and magnetic field strengths for the intrinsically steep case. We find that the ratio between the equipartition magnetic field strengths and those derived through synchrotron/inverse-Compton model fitting remains consistent with previous findings and show that the observed departure from equipartition may in some cases provide a solution to the spectral versus dynamical age disparity.

LS I +61 303 is a gamma-ray binary that exhibits an outburst at GHz frequencies each orbital cycle of $\approx$ 26.5 d and a superorbital modulation with a period of $\approx$ 4.6 yr. We have performed a detailed study of the low-frequency radio emission of LS I +61 303 by analysing all the archival GMRT data at 150, 235 and 610 MHz, and conducting regular LOFAR observations within the Radio Sky Monitor (RSM) at 150 MHz. We have detected the source for the first time at 150 MHz, which is also the first detection of a gamma-ray binary at such a low frequency. We have obtained the light-curves of the source at 150, 235 and 610 MHz, all of them showing orbital modulation. The light-curves at 235 and 610 MHz also show the existence of superorbital variability. A comparison with contemporaneous 15-GHz data shows remarkable differences with these light-curves. At 15 GHz we see clear outbursts, whereas at low frequencies we see variability with wide maxima. The light-curve at 235 MHz seems to be anticorrelated with the one at 610 MHz, implying a shift of $\sim$ 0.5 orbital phases in the maxima. We model the shifts between the maxima at different frequencies as due to changes in the physical parameters of the emitting region assuming either free-free absorption or synchrotron self-absorption, obtaining expansion velocities for this region close to the stellar wind velocity with both mechanisms.

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.

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.

This chapter provides an overview of the possibilities for transient and variable-source astrophysics with the Square Kilometre Array. While subsequent chapters focus on the astrophysics of individual events, we focus on the broader picture, and how to maximise the science coming from the telescope. The SKA as currently designed will be a fantastic and ground-breaking facility for radio transient studies, but the scientifc yield will be dramatically increased by the addition of (i) near-real-time commensal searches of data streams for events, and (ii) on occasion, rapid robotic response to Target-of-Opprtunity style triggers.

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.

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

Current and future astronomical survey facilities provide a remarkably rich opportunity for transient astronomy, combining unprecedented fields of view with high sensitivity and the ability to access previously unexplored wavelength regimes. This is particularly true of LOFAR, a recently-commissioned, low-frequency radio interferometer, based in the Netherlands and with stations across Europe. The identification of and response to transients is one of LOFAR's key science goals. However, the large data volumes which LOFAR produces, combined with the scientific requirement for rapid response, make automation essential. To support this, we have developed the LOFAR Transients Pipeline, or TraP. The TraP ingests multi-frequency image data from LOFAR or other instruments and searches it for transients and variables, providing automatic alerts of significant detections and populating a lightcurve database for further analysis by astronomers. Here, we discuss the scientific goals of the TraP and how it has been designed to meet them. We describe its implementation, including both the algorithms adopted to maximize performance as well as the development methodology used to ensure it is robust and reliable, particularly in the presence of artefacts typical of radio astronomy imaging. Finally, we report on a series of tests of the pipeline carried out using simulated LOFAR observations with a known population of transients.

We present the highest-quality polarisation profiles to date of 16 non-recycled pulsars and four millisecond pulsars, observed below 200 MHz with the LOFAR high-band antennas. Based on the observed profiles, we perform an initial investigation of expected observational effects resulting from the propagation of polarised emission in the pulsar magnetosphere and the interstellar medium. The predictions of magnetospheric birefringence in pulsars have been tested using spectra of the pulse width and fractional polarisation from multifrequency data. The derived spectra offer only partial support for the expected effects of birefringence on the polarisation properties, with only about half of our sample being consistent with the model's predictions. It is noted that for some pulsars these measurements are contaminated by the effects of interstellar scattering. For a number of pulsars in our sample, we have observed significant variations in the amount of Faraday rotation as a function of pulse phase, which is possibly an artefact of scattering. These variations are typically two orders of magnitude smaller than that observed at 1400 MHz by Noutsos et al. (2009), for a different sample of southern pulsars. In this paper we present a possible explanation for the difference in magnitude of this effect between the two frequencies, based on scattering. Finally, we have estimated the magnetospheric emission heights of low-frequency radiation from four pulsars, based on the phase lags between the flux-density and the PA profiles, and the theoretical framework of Blaskiewicz, Cordes & Wasserman (1991). These estimates yielded heights of a few hundred km; at least for PSR B1133+16, this is consistent with emission heights derived based on radius-to-frequency mapping, but is up to a few times larger than the recent upper limit based on pulsar timing.

Aims: We have analyzed low frequency radio data of tidal disruption event (TDE) Swift J1644+57 to search for a counterpart. We consider how brief transient signals (on the order of seconds or minutes) originating from this location would appear in our data. We also consider how automatic radio frequency interference (RFI) flagging at radio telescope observatories might affect these and other transient observations in the future, particularly with brief transients of a few seconds duration. Methods: We observed the field in the low-frequency regime at 149 MHz with data obtained over several months with the Low Frequency Array (LOFAR). We also present simulations where a brief transient is injected into the data in order to see how it would appear in our measurement sets, and how it would be affected by RFI flagging. Finally, both based on simulation work and the weighted average of the observed background over the course of the individual observations, we present the possibility of brief radio transients in the data. Results: Our observations of Swift J1644+57 yielded no detection of the source and a peak flux density at this position of 24.7 $\pm$ 8.9 mJy. Our upper limit on the transient rate of the snapshot surface density in this field at sensitivities < 0.5 Jy is $\rho < 2.2 \times10^{-2}$ deg$^{-2}$. We also conclude that we did not observe any brief transient signals originating specifically from the Swift J1644+57 source itself, and searches for such transients are severely limited by automatic RFI flagging algorithms which flag transients of less than 2 minutes duration. As such, careful consideration of RFI flagging techniques must occur when searching for transient signals.

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.

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.

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.

During the Cassini spacecraft's flyby of Jupiter (October, 2000-March, 2001), the Ultraviolet Imaging Spectrograph (UVIS) produced an extensive dataset consisting of 3,349 spectrally dispersed images of the Io plasma torus. Here we present an example of the raw data and representative EUV spectra (561{\AA}-1181{\AA}) of the torus, obtained on October 1, 2000 and November 14, 2000. For most of the flyby period, the entire Io torus fit within the UVIS field-of-view, enabling the measurement of the total power radiated from the torus in the extreme ultraviolet. A typical value for the total power radiated in the wavelength range of 580{\AA}-1181{\AA} is 1.7x10^12 W, with observed variations of up to 25%. Several brightening events were observed. These events lasted for roughly 20 hours, during which time the emitted power increased rapidly by ~20% before slowly returning to the pre-event level. Observed variations in the relative intensities of torus spectral features provide strong evidence for compositional changes in the torus plasma with time. Spatial profiles of the EUV emission show no evidence for a sharply peaked "ribbon" feature. The ratio of the brightness of the dusk ansa to the brightness of the dawn ansa is observed to be highly variable, with an average value of 1.30. Weak longitudinal variations in the brightness of the torus ansae were observed at the 2% level.

On January 14, 2001, shortly after the Cassini spacecraft's closest approach to Jupiter, the Ultraviolet Imaging Spectrometer (UVIS) made a radial scan through the midnight sector of Io plasma torus. The Io torus has not been previously observed at this local time. The UVIS data consist of 2-D spectrally dispersed images of the Io plasma torus in the wavelength range of 561{\AA}-1912{\AA}. We developed a spectral emissions model that incorporates the latest atomic physics data contained in the CHIANTI database in order to derive the composition of the torus plasma as a function of radial distance. Electron temperatures derived from the UVIS torus spectra are generally less than those observed during the Voyager era. We find the torus ion composition derived from the UVIS spectra to be significantly different from the composition during the Voyager era. Notably, the torus contains substantially less oxygen, with a total oxygen-to-sulfur ion ratio of 0.9. The average ion charge state has increased to 1.7. We detect S V in the Io torus at the 3{\sigma} level. S V has a mixing ratio of 0.5%. The spectral emission model used in can approximate the effects of a non-thermal distribution of electrons. The ion composition derived using a kappa distribution of electrons is identical to that derived using a Maxwellian electron distribution; however, the kappa distribution model requires a higher electron column density to match the observed brightness of the spectra. The derived value of the kappa parameter decreases with radial distance and is consistent with the value of {\kappa}=2.4 at 8 RJ derived by the Ulysses URAP instrument (Meyer-Vernet et al., 1995). The observed radial profile of electron column density is consistent with a flux tube content, NL^2, that is proportional to r^-2.

We develop a method of constraining the cosmic string tension $G\mu$ which uses the Canny edge detection algorithm as a means of searching CMB temperature maps for the signature of the Kaiser-Stebbins effect. We test the potential of this method using high resolution, simulated CMB temperature maps. By modeling the future output from the South Pole Telescope project (including anticipated instrumental noise), we find that cosmic strings with $G\mu > 5.5\times10^{-8}$ could be detected.

Motivated by the string gas cosmological model, which predicts a blue tilt of the primordial gravitational wave spectrum, we examine the constraints imposed by current and planned observations on a blue tilted tensor spectrum. Starting from an expression for the primordial gravitational wave spectrum normalized using cosmic microwave background observations, pulsar timing, direct detection and nucleosynthesis bounds are examined. If we assume a tensor to scalar ratio on scales of the CMB which equals the current observational upper bound, we obtain from these current observations constraints on the tensor spectral index of $n_{T} \lesssim 0.79$, $n_{T} \lesssim 0.53$, and $n_{T} \lesssim 0.15$ respectively.