Search

We present a comprehensive study of the rotational and emission properties of PSR J0738$-$4042 using a combination of observations taken by the Deep Space Network, Hartebeesthoek, Parkes (Murriyang) and Molonglo observatories between 1972 and 2023. Our timing of the pulsar is motivated by previously reported profile/spin-down events that occurred in September 2005 and December 2015, which result in an anomalously large braking index of $n = 23300 \pm 1800$. Using a Gaussian process regression framework, we develop continuous models for the evolution of the pulsar spin-down rate ($\dot{\nu}$) and profile shape. We find that the pulse profile variations are similar regardless of radio observing frequency and polarisation. Small-scale differences can be ascribed to changes in the interstellar medium along the line of sight and frequency-dependent changes in magnetospheric radio emission height. No new correlated spin-down or profile events were identified in our extended dataset. However, we found that the disappearance of a bright emission component in the leading edge of archival profiles between 1981-1988 was not associated with a substantial change in $\dot{\nu}$. This marks a notable departure from the previous profile/spin-down events in this pulsar. We discuss the challenges these observations pose for physical models and conclude that interactions between the pulsar and in-falling asteroids or a form of magnetospheric state-switching with a long periodicity are plausible explanations., Comment: 14 pages, 8 figures. Fixed typo on page 5. Accepted for publication in MNRAS

A Faraday rotation measure (RM) catalogue, or RM Grid, is a valuable resource for the study of cosmic magnetism. Using the second data release (DR2) from the LOFAR Two-metre Sky Survey (LoTSS), we have produced a catalogue of 2461 extragalactic high-precision RM values across 5720 deg$^{2}$ of sky (corresponding to a polarized source areal number density of $\sim$0.43 deg$^{-2}$). The linear polarization and RM properties were derived using RM synthesis from the Stokes $Q$ and $U$ channel images at an angular resolution of 20'' across a frequency range of 120 to 168 MHz with a channel bandwidth of 97.6 kHz. The fraction of total intensity sources ($>1$ mJy beam$^{-1}$) found to be polarized was $\sim$0.2%. The median detection threshold was 0.6 mJy beam$^{-1}$ ($8\sigma_{QU}$), with a median RM uncertainty of 0.06 rad m$^{-2}$ (although a systematic uncertainty of up to 0.3 rad m$^{-2}$ is possible, after the ionosphere RM correction). The median degree of polarization of the detected sources is 1.8%, with a range of 0.05% to 31%. Comparisons with cm-wavelength RMs indicate minimal amounts of Faraday complexity in the LoTSS detections, making them ideal sources for RM Grid studies. Host galaxy identifications were obtained for 88% of the sources, along with redshifts for 79% (both photometric and spectroscopic), with the median redshift being 0.6. The focus of the current catalogue was on reliability rather than completeness, and we expect future versions of the LoTSS RM Grid to have a higher areal number density. In addition, 25 pulsars were identified, mainly through their high degrees of linear polarization., Comment: 20 pages, 23 figures, 3 tables. Accepted for publication in MNRAS. Go to https://lofar-mksp.org/data/ to download the catalogue and associated data products

The radio polarization properties of the pulsar population are only superficially captured by the conventional picture of pulsar radio emission. We study the broadband polarization of 271 young radio pulsars, focusing particularly on circular polarization, using high quality observations made with the Ultra-Wideband Low receiver on Murriyang, the Parkes radio telescope. We seek to encapsulate polarization behaviour on a population scale by defining broad categories for frequency- and phase-dependent polarization evolution, studying the co-occurrences of these categorizations and comparing them with average polarization measurements and spin-down energy ($\dot{E}$). This work shows that deviations of the linear polarization position angle (PA) from the rotating vector model (RVM) are linked to the presence of circular polarization features and to frequency evolution of the polarization. Polarization fraction, circular polarization contribution and profile complexity all evolve with $\dot{E}$ across the population, with the profiles of high-$\dot{E}$ pulsars being simple and highly linearly polarized. The relationship between polarization fraction and circular contribution is also seen to evolve such that highly polarized profiles show less variation in circular contribution with frequency than less strongly polarized profiles. This evolution is seen both across the population and across frequency for individual sources. Understanding pulsar radio polarization requires detailed study of individual sources and collective understanding of population-level trends. For the former, we provide visualizations of their phase- and frequency-resolved polarization parameters. For the latter, we have highlighted the importance of including the impact of circular polarization and of $\dot{E}$., Comment: 20 pages, 7 figures, supplementary material available at the MNRAS website or by contacting the lead author

Discovering radio pulsars, particularly millisecond pulsars (MSPs), is important for a range of astrophysical applications, such as testing theories of gravity or probing the magneto-ionic interstellar medium. We aim to discover pulsars that may have been missed in previous pulsar searches by leveraging known pulsar observables (primarily polarisation) in the sensitive, low-frequency radio images from the Low-Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS), and have commenced the Targeted search, using LoTSS images, for polarised pulsars (TULIPP) survey. For this survey, we identified linearly and circularly polarised point sources with flux densities brighter than 2 mJy in LoTSS images at a centre frequency of 144 MHz with a 48 MHz bandwidth. Over 40 known pulsars, half of which are MSPs, were detected as polarised sources in the LoTSS images and excluded from the survey. We have obtained beam-formed LOFAR observations of 30 candidates, which were searched for pulsations using coherent de-dispersion. Here, we present the results of the first year of the TULIPP survey. We discovered two pulsars, PSRs J1049+5822 and J1602+3901, with rotational periods of P=0.73 s and 3.7 ms, respectively. We also detected a further five known pulsars (two slowly-rotating pulsars and three MSPs) for which accurate sky positions were not available to allow a unique cross-match with LoTSS sources. This targeted survey presents a relatively efficient method by which pulsars, particularly MSPs, may be discovered using the flexible observing modes of sensitive radio telescopes such as the Square Kilometre Array and its pathfinders/precursors, particularly since wide-area all-sky surveys using coherent de-dispersion are currently computationally infeasible., Comment: 15 pages, 7 figures, 3 tables. Accepted for publication in A&A

We present radio observations of the most slowly rotating known radio pulsar PSR J0250+5854. With a 23.5 s period, it is close, or even beyond, the $P$-$\dot{P}$ diagram region thought to be occupied by active pulsars. The simultaneous observations with FAST, the Chilbolton and Effelsberg LOFAR international stations, and NenuFAR represent a five-fold increase in the spectral coverage of this object, with the detections at 1250 MHz (FAST) and 57 MHz (NenuFAR) being the highest- and lowest-frequency published respectively to date. We measure a flux density of $4\pm2$ $\mu$Jy at 1250 MHz and an exceptionally steep spectral index of $-3.5^{+0.2}_{-1.5}$, with a turnover below $\sim$95 MHz. In conjunction with observations of this pulsar with the GBT and the LOFAR Core, we show that the intrinsic profile width increases drastically towards higher frequencies, contrary to the predictions of conventional radius-to-frequency mapping. We examine polarimetric data from FAST and the LOFAR Core and conclude that its polar cap radio emission is produced at an absolute height of several hundreds of kilometres around 1.5 GHz, similar to other rotation-powered pulsars across the population. Its beam is significantly underfilled at lower frequencies, or it narrows because of the disappearance of conal outriders. Finally, the results for PSR J0250+5854 and other slowly spinning rotation-powered pulsars are contrasted with the radio-detected magnetars. We conclude that magnetars have intrinsically wider radio beams than the slow rotation-powered pulsars, and that consequently the latter's lower beaming fraction is what makes objects such as PSR J0250+5854 so scarce., Comment: 14 pages, 6 figures

We present high signal-to-noise, full polarization pulse profiles for 40 bright, 'slowly'-rotating (non-recycled) pulsars using the new Ultra-Wideband Low-frequency (UWL; 704-4032 MHz) receiver on the Parkes radio telescope. We obtain updated and accurate interstellar medium parameters towards these pulsars (dispersion measures and Faraday rotation measures), and reveal Faraday dispersion towards PSR J1721-3532 caused by interstellar scattering. We find general trends in the pulse profiles including decreasing fractional linear polarization and increasing degree of circular polarization with increasing frequency, consistent with previous studies, while also revealing new features and frequency evolution. This demonstrates results that can be obtained using UWL monitoring observations of slow pulsars, which are valuable for improving our understanding of pulsar emission and the intervening interstellar medium. The calibrated data products are publicly available., Comment: 20 pages, 5 figures, 2 tables, 25 pages of supporting information with 2 figures. Accepted for publication in MNRAS. doi:10.1093/mnras/stab861

We describe the ongoing Relativistic Binary programme (RelBin), a part of the MeerTime large survey project with the MeerKAT radio telescope. RelBin is primarily focused on observations of relativistic effects in binary pulsars to enable measurements of neutron star masses and tests of theories of gravity. We selected 25 pulsars as an initial high priority list of targets based on their characteristics and observational history with other telescopes. In this paper, we provide an outline of the programme, present polarisation calibrated pulse profiles for all selected pulsars as a reference catalogue along with updated dispersion measures. We report Faraday rotation measures for 24 pulsars, twelve of which have been measured for the first time. More than a third of our selected pulsars show a flat position angle swing confirming earlier observations. We demonstrate the ability of the Rotating Vector Model (RVM), fitted here to seven binary pulsars, including the Double Pulsar (PSR J0737$-$3039A), to obtain information about the orbital inclination angle. We present a high time resolution light curve of the eclipse of PSR J0737$-$3039A by the companion's magnetosphere, a high-phase resolution position angle swing for PSR J1141$-$6545, an improved detection of the Shapiro delay of PSR J1811$-$2405, and pulse scattering measurements for PSRs J1227$-$6208, J1757$-$1854, and J1811$-$1736. Finally, we demonstrate that timing observations with MeerKAT improve on existing data sets by a factor of, typically, 2-3, sometimes by an order of magnitude., Comment: 21 pages 16 figures, published in MNRAS (replaced earlier submission after small changes added in proofs)

The major programme for observing young, non-recycled pulsars with the Parkes telescope has transitioned from a narrow-band system to an ultra-wideband system capable of observing between 704 and 4032 MHz. We report here on the initial two years of observations with this receiver. Results include dispersion measure (DM) and Faraday rotation measure (RM) variability with time, determined with higher precision than hitherto, flux density measurements and the discovery of several nulling and mode changing pulsars. PSR J1703-4851 is shown to be one of a small subclass of pulsars that has a weak and a strong mode which alternate rapidly in time. PSR J1114-6100 has the fourth highest |RM| of any known pulsar despite its location far from the Galactic Centre. PSR J1825-1446 shows variations in both DM and RM likely due to its motion behind a foreground supernova remnant., Comment: Accepted for publication in MNRAS. This version fixes transcription errors in Table 4

The Parkes telescope has been monitoring 286 radio pulsars approximately monthly since 2007 at an observing frequency of 1.4 GHz. The wide dispersion measure (DM) range of the pulsar sample and the uniformity of the observing procedure make the data-set extremely valuable for studies of flux density variability and the interstellar medium. Here, we present flux density measurements and modulation indices of these pulsars over this period. We derive the structure function from the light curves and discuss the contributions to it from measurement noise, intrinsic variability and interstellar scintillation. Despite a large scatter, we show that the modulation index is inversely correlated with DM, and can be generally described by a power-law with an index of $\sim-0.7$ covering DMs from $\sim10$ to 1000 cm$^{-3}$ pc. We present refractive timescales and/or lower limits for a group of 42 pulsars. These often have values significantly different from theoretical expectations, indicating the complex nature of the interstellar medium along individual lines of sight. In particular, local structures and non-Kolmogorov density fluctuations are likely playing important roles in the observed flux density variation of many of these pulsars., Comment: 27 pages, 46 figures, accepted for publication in MNRAS

Radio pulses from pulsars are affected by plasma dispersion, which results in a frequency-dependent propagation delay. Variations in the magnitude of this effect lead to an additional source of red noise in pulsar timing experiments, including pulsar timing arrays that aim to detect nanohertz gravitational waves. We aim to quantify the time-variable dispersion with much improved precision and characterise the spectrum of these variations. We use the pulsar timing technique to obtain highly precise dispersion measure (DM) time series. Our dataset consists of observations of 36 millisecond pulsars, which were observed for up to 7.1 years with the LOFAR telescope at a centre frequency of ~150 MHz. Seventeen of these sources were observed with a weekly cadence, while the rest were observed at monthly cadence. We achieve a median DM precision of the order of 10^-5 cm^-3 pc for a significant fraction of our sources. We detect significant variations of the DM in all pulsars with a median DM uncertainty of less than 2x10^-4 cm^-3 pc. The noise contribution to pulsar timing experiments at higher frequencies is calculated to be at a level of 0.1-10 us at 1.4 GHz over a timespan of a few years, which is in many cases larger than the typical timing precision of 1 us or better that PTAs aim for. We found no evidence for a dependence of DM on radio frequency for any of the sources in our sample. The DM time series we obtained using LOFAR could in principle be used to correct higher-frequency data for the variations of the dispersive delay. However, there is currently the practical restriction that pulsars tend to provide either highly precise times of arrival (ToAs) at 1.4 GHz or a high DM precision at low frequencies, but not both, due to spectral properties. Combining the higher-frequency ToAs with those from LOFAR to measure the infinite-frequency ToA and DM would improve the result., Comment: 16 pages, 8 figures, 2 tables

We present deep polarimetric observations of the European Large Area ISO Survey-North 1 (ELAIS-N1) field using the Low Frequency Array (LOFAR) at 114.9-177.4 MHz. The ELAIS-N1 field is part of the LOFAR Two-metre Sky Survey deep fields data release I. For six eight-hour observing epochs, we align the polarization angles and stack the 20"-resolution Stokes $Q$, $U$-parameter data cubes. This produces a 16 deg$^2$ image with 1$\sigma_{\rm QU}$ sensitivity of 26 $\mu$Jy/beam in the central area. In this paper, we demonstrate the feasibility of the stacking technique, and we generate a catalog of polarized sources in ELAIS-N1 and their associated Faraday rotation measures (RMs). While in a single-epoch observation we detect three polarized sources, this number increases by a factor of about three when we consider the stacked data, with a total of ten sources. This yields a surface density of polarized sources of one per 1.6 deg$^2$. The Stokes $I$ images of three of the ten detected polarized sources have morphologies resembling those of FR I radio galaxies. This represents a greater fraction of this type of source than previously found, which suggests that more sensitive observations may help with their detection., Comment: This paper is part of the 1st data release of the LoTSS Deep Fields. Accepted for publication in A&A. 14 pages, 9 figures, 7 tables

We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities., Comment: 20 pages, 16 Figures, 4 Tables, accepted for publication in PASA

The low-frequency linearly-polarised radio source population is largely unexplored. However, a renaissance in low-frequency polarimetry has been enabled by pathfinder and precursor instruments for the Square Kilometre Array. In this second paper from the POlarised GaLactic and Extragalactic All-Sky Murchison Widefield Array (MWA) Survey -- the POlarised GLEAM Survey, or POGS -- we present the results from our all-sky MWA Phase I Faraday Rotation Measure survey. Our survey covers nearly the entire Southern sky in the Declination range $-82^{\circ}$ to $+30^{\circ}$ at a resolution between around three and seven arcminutes (depending on Declination) using data in the frequency range 169$-$231 MHz. We have performed two targeted searches: the first covering 25,489 square degrees of sky, searching for extragalactic polarised sources; the second covering the entire sky South of Declination $+30^{\circ}$, searching for known pulsars. We detect a total of 517 sources with 200 MHz linearly-polarised flux densities between 9.9 mJy and 1.7 Jy, of which 33 are known radio pulsars. All sources in our catalogues have Faraday rotation measures in the range $-328.07$ rad m$^{-2}$ to $+279.62$ rad m$^{-2}$. The Faraday rotation measures are broadly consistent with results from higher-frequency surveys, but with typically more than an order of magnitude improvement in the precision, highlighting the power of low-frequency polarisation surveys to accurately study Galactic and extragalactic magnetic fields. We discuss the properties of our extragalactic and known-pulsar source population, how the sky distribution relates to Galactic features, and identify a handful of new pulsar candidates among our nominally extragalactic source population., Comment: Replacement of previous version. Only change is minor updates to catalogues (see ancillary files) which now contain bib code of finalised manuscript (published in PASA). Manuscript has 31 pages, 10 figures, 5 tables. Four Appendices are included in the ancillary material, showing further Figures, continuum spectra for a handful of selected sources, and RM spectra for all 517 sources

We present a comparative study of the low-frequency eclipses of spider (compact, irradiating binary) PSRs B1957+20 and J1816+4510. Combining these data with those of three other eclipsing systems we study the frequency dependence of the eclipse duration. PSRs B1957+20 and J1816+4510 have similar orbital properties, but the companions to the pulsars have masses that differ by an order of magnitude. A dedicated campaign to simultaneously observe the pulsed and imaged continuum flux densities throughout the eclipses reveals many similarities between the excess material within the two binaries, irrespective of the companion star properties. The observations show that the pulsar fluxes are removed from the line of sight throughout the main body of the eclipses. For PSR J1816+4510 we present the first direct evidence of an eclipse mechanism that transitions from one that removes the pulsar flux from the line of sight to one that merely smears out pulsations, and claim that this is a consequence of scattering in a tail of material flowing behind the companion. Inferred mass loss rates from the companion stars are found to be $\dot{M}_{\text{C}} \sim 10^{-12}~M_\odot$~yr$^{-1}$ and $\dot{M}_{\text{C}} \sim 2 \times 10^{-13}~M_\odot$~yr$^{-1}$ for PSR B1957+20 and PSR J1816+4510, respectively; seemingly too low to evaporate the stars within Hubble time. Measurements of eclipse durations over a wide range of radio-frequencies show a significant dependence of eclipse duration on frequency for all pulsars, with wider eclipses at lower-frequencies. These results provide a marked improvement in the observational constraints available for theoretical studies of the eclipse mechanisms., Comment: Accepted for publication in MNRAS

We report here on initial results from the Thousand Pulsar Array (TPA) programme, part of the Large Survey Project "MeerTime" on the MeerKAT telescope. The interferometer is used in tied-array mode in the band from 856 to 1712~MHz, and the wide band coupled with the large collecting area and low receiver temperature make it an excellent telescope for the study of radio pulsars. The TPA is a 5 year project which aims to observe (a) more than 1000 pulsars to obtain high-fidelity pulse profiles, (b) some 500 of these pulsars over multiple epochs, (c) long sequences of single-pulse trains from several hundred pulsars. The scientific outcomes from the programme will include determination of pulsar geometries, the location of the radio emission within the pulsar magnetosphere, the connection between the magnetosphere and the crust and core of the star, tighter constraints on the nature of the radio emission itself as well as interstellar medium studies. First results presented here include updated dispersion measures, 26 pulsars with Faraday rotation measures derived for the first time and a description of interesting emission phenomena observed thus far., Comment: Accepted by MNRAS

In this paper we report on $\sim10$ years of observations of PSR J2051$-$0827, at radio frequencies in the range 110--4032 MHz. We investigate the eclipse phenomena of this black widow pulsar using model fits of increased dispersion and scattering of the pulsed radio emission as it traverses the eclipse medium. These model fits reveal variability in dispersion features on timescales as short as the orbital period, and previously unknown trends on timescales of months--years. No clear patterns are found between the low-frequency eclipse widths, orbital period variations and trends in the intra-binary material density. Using polarisation calibrated observations we present the first available limits on the strength of magnetic fields within the eclipse region of this system; the average line of sight field is constrained to be $10^{-4}$ G $\lesssim B_{||} \lesssim 10^2$ G, while for the case of a field directed near-perpendicular to the line of sight we find $B_{\perp} \lesssim 0.3$ G. Depolarisation of the linearly polarised pulses during the eclipse is detected and attributed to rapid rotation measure fluctuations of $\sigma_{\text{RM}} \gtrsim 100$ rad m$^{-2}$ along, or across, the line of sights averaged over during a sub-integration. The results are considered in the context of eclipse mechanisms, and we find scattering and/or cyclotron absorption provide the most promising explanation, while dispersion smearing is conclusively ruled out. Finally, we estimate the mass loss rate from the companion to be $\dot{M}_{\text{C}} \sim 10^{-12} M_\odot$ yr$^{-1}$, suggesting that the companion will not be fully evaporated on any reasonable timescale.

We present the results from the low-frequency (40--78 MHz) extension of the first LOFAR pulsar census of non-recycled pulsars. We have used the Low-Band Antennas of the LOFAR core stations to observe 87 pulsars out of 158 that have been detected previously with the High-Band Antennas. Forty-three pulsars have been detected and we present here their flux densities and flux-calibrated profiles. Seventeen of these pulsars have not been, to our knowledge, detected before at such low frequencies. We re-calculate the spectral indices using the new low-frequency flux density measurements from the LOFAR census and discuss the prospects of studying pulsars at the very low frequencies with the current and upcoming facilities, such as NenuFAR., Comment: submitted to A&A

The smooth spin-down of young pulsars is perturbed by two non-deterministic phenomenon, glitches and timing noise. Although the timing noise provides insights into nuclear and plasma physics at extreme densities, it acts as a barrier to high-precision pulsar timing experiments. An improved methodology based on Bayesian inference is developed to simultaneously model the stochastic and deterministic parameters for a sample of 85 high-$\dot{E}$ radio pulsars observed for $\sim$ 10 years with the 64-m Parkes radio telescope. Timing noise is known to be a red process and we develop a parametrization based on the red-noise amplitude ($A_{\rm red}$) and spectral index ($\beta$). We measure the median $A_{\rm red}$ to be $-10.4^{+1.8}_{-1.7}$ yr$^{3/2}$ and $\beta$ to be $-5.2^{+3.0}_{-3.8}$ and show that the strength of timing noise scales proportionally to $\nu^{1}|\dot{\nu}|^{-0.6\pm0.1}$, where $\nu$ is the spin frequency of the pulsar and $\dot{\nu}$ its spin-down rate. Finally, we measure significant braking indices for 19 pulsars, proper motions for two pulsars and discuss the presence of periodic modulation in the arrival times of five pulsars.

Rotating Radio Transients (RRATs) represent a relatively new class of pulsar, primarily characterised by their sporadic bursting emission of single pulses on time scales of minutes to hours. In addition to the difficulty involved in detecting these objects, low-frequency ($<$300 MHz) observations of RRATs are sparse, which makes understanding their broadband emission properties in the context of the normal pulsar population problematic. Here, we present the simultaneous detection of RRAT J2325-0530 using the Murchison Widefield Array (154 MHz) and Parkes radio telescope (1.4 GHz). On a single-pulse basis, we produce the first polarimetric profile of this pulsar, measure the spectral index ($\alpha=-2.2\pm 0.1$), pulse energy distributions, and present the pulse rates in the context of detections in previous epochs. We find that the distribution of time between subsequent pulses is consistent with a Poisson process and find no evidence of clustering over the $\sim$1.5 hr observations. Finally, we are able to quantify the scintillation properties of RRAT J2325-0530 at 1.4 GHz, where the single pulses are modulated substantially across the observing bandwidth, and show that this characterisation is feasible even with irregular time sampling as a consequence of the sporadic emission behaviour., Comment: 18 pages, 8 figures, 5 tables, accepted for publication in PASA

We present an overview of the LOFAR Tied-Array All-Sky Survey (LOTAAS) for radio pulsars and fast transients. The survey uses the high-band antennas of the LOFAR Superterp, the dense inner part of the LOFAR core, to survey the northern sky (dec > 0 deg) at a central observing frequency of 135 MHz. A total of 219 tied-array beams (coherent summation of station signals, covering 12 square degrees), as well as three incoherent beams (covering 67 square degrees) are formed in each survey pointing. For each ofthe 222 beams, total intensity is recorded at 491.52 us time resolution. Each observation integrates for 1 hr and covers 2592 channels from 119 to 151 MHz. This instrumental setup allows LOTAAS to reach a detection threshold of 1 to 5 mJy for periodic emission. Thus far, the LOTAAS survey has resulted in the discovery of 73 radio pulsars. Among these are two mildly recycled binary millisecond pulsars (P = 13 and 33 ms), as well as the slowest-spinning radio pulsar currently known (P = 23.5 s). The survey has thus far detected 311 known pulsars, with spin periods ranging from 4 ms to 5.0 s and dispersion measures from 3.0 to 217 pc/cc. Known pulsars are detected at flux densities consistent with literature values. We find that the LOTAAS pulsar discoveries have, on average, longer spin periods than the known pulsar population. This may reflect different selection biases between LOTAAS and previous surveys, though it is also possible that slower-spinning pulsars preferentially have steeper radio spectra. LOTAAS is the deepest all-sky pulsar survey using a digital aperture array; we discuss some of the lessons learned that can inform the approach for similar surveys using future radio telescopes such as the Square Kilometre Array., Comment: 22 pages, 9 figures, 3 tables. Accepted for publication in A&A

The Murchison Widefield Array is a low-frequency Square Kilometre Array precursor located at the Murchison Radio-astronomy Observatory in Western Australia. Primarily designed as an imaging telescope, but with a flexible signal path, the capabilities of this telescope have recently been extended to include off-line incoherent and tied-array beam formation using recorded antenna voltages. This has provided the capability for high-time and frequency resolution observations, including a pulsar science program. This paper describes the algorithms and pipeline that we have developed to form the tied array beam products from the summation of calibrated signals of the antenna elements, and presents example polarimetric profiles for PSRs J0437-4715 and J1900-2600 at 185 MHz., Comment: 14 pages, 3 figures, accepted for publication in PASA

Polarimetric studies of pulsars at low radio frequencies provide important observational insights into the pulsar emission mechanism and beam models, and probe the properties of the magneto-ionic interstellar medium (ISM). Aperture arrays are the main form of next-generation low-frequency telescopes, including the Murchison Widefield Array (MWA). These require a distinctly different approach to data processing (e.g. calibration and beamforming) compared to traditional dish antennas. As the second paper of this series, we present a verification of the MWA's pulsar polarimetry capability, using two bright southern pulsars, PSRs J0742-2822 and J1752-2806. Our observations simultaneously cover multiple frequencies (76-313 MHz) and were taken at multiple zenith angles during a single night for each pulsar. We show that the MWA can be reliably calibrated for zenith angles < 45 degree and frequencies < 270 MHz. We present the polarimetric profiles for PSRs J0742-2822 and J1752-2806 at frequencies lower than 300 MHz for the first time, along with an analysis of the linear polarisation degree and pulse profile evolution with frequency. For PSR J0742-2822, the measured degree of linear polarisation shows a rapid decrease at low frequencies, in contrast with the generally expected trend, which can be attributed to depolarisation effects from small-scale, turbulent, magneto-ionic ISM components. This effect has not been widely explored for pulsars in general, and will be further investigated in future work., Comment: accepted for publication in PASA

We determined Faraday rotation measures (RMs) towards 137 pulsars in the northern sky, using Low-Frequency Array (LOFAR) observations at 110-190 MHz. This low-frequency RM catalogue, the largest to date, improves the precision of existing RM measurements on average by a factor of 20 - due to the low frequency and wide bandwidth of the data, aided by the RM synthesis method. We report RMs towards 25 pulsars for the first time. The RMs were corrected for ionospheric Faraday rotation to increase the accuracy of our catalogue to approximately 0.1 rad m$^{\rm -2}$. The ionospheric RM correction is currently the largest contributor to the measurement uncertainty. In addition, we find that the Faraday dispersion functions towards pulsars are extremely Faraday thin - mostly less than 0.001 rad m$^{\rm -2}$. We use these new precise RM measurements (in combination with existing RMs, dispersion measures, and distance estimates) to estimate the scale height of the Galactic halo magnetic field: 2.0$\pm$0.3 kpc for Galactic quadrants I and II above and below the Galactic plane (we also evaluate the scale height for these regions individually). Overall, our initial low-frequency catalogue provides valuable information about the 3-D structure of the Galactic magnetic field., Comment: 31 pages, including 13 figures, 6 tables. Accepted for publication in MNRAS

The rare intermittent pulsars pose some of the most challenging questions surrounding the pulsar emission mechanism, but typically have relatively minimal low-frequency ($\lesssim$ 300 MHz) coverage. We present the first low-frequency detection of the intermittent pulsar J1107-5907 with the Murchison Widefield Array (MWA) at 154 MHz and the simultaneous detection from the recently upgraded Molonglo Observatory Synthesis Telescope (UTMOST) at 835 MHz, as part of an on-going observing campaign. During a 30-minute simultaneous observation, we detected the pulsar in its bright emission state for approximately 15 minutes, where 86 and 283 pulses were detected above a signal-to-noise threshold of 6 with the MWA and UTMOST, respectively. Of the detected pulses, 51 had counterparts at both frequencies and exhibited steep spectral indices for both the bright main pulse component and the precursor component. We find that the bright state pulse energy distribution is best parameterised by a log-normal distribution at both frequencies, contrary to previous results which suggested a power law distribution. Further low-frequency observations are required in order to explore in detail aspects such as pulse-to-pulse variability, intensity modulations and to better constrain the signal propagation effects due to the interstellar medium and intermittency characteristics at these frequencies. The spectral index, extended profile emission covering a large fraction of pulse longitude, and the broadband intermittency of PSR J1107-5907 suggests that future low-frequency pulsar searches, for instance those planned with SKA-Low, will be in an excellent position to find and investigate new pulsars of this type., Comment: 19 pages, 7 figures, 4 tables, accepted for publication in ApJ; minor updates to keep in-line with published version

We present the discovery of PSR J0250+5854, a radio pulsar with a spin period of 23.5 s. This is the slowest-spinning radio pulsar known. PSR J0250+5854 was discovered by the LOFAR Tied-Array All-Sky Survey (LOTAAS), an all-Northern-sky survey for pulsars and fast transients at a central observing frequency of 135 MHz. We subsequently detected pulsations from the pulsar in the interferometric images of the LOFAR Two-metre Sky Survey, allowing for sub-arcsecond localization. This, along with a pre-discovery detection 2 years prior, allowed us to measure the spin-period derivative to be $\dot{P}=2.7 \times 10^{-14}$ s s$^{-1}$. The observed spin period derivative of PSR J0250+5854 indicates a surface magnetic field strength, characteristic age and spin-down luminosity of $2.6 \times 10^{13}$G, $13.7$ Myr and $8.2 \times 10^{28}$ erg s$^{-1}$ respectively, for a dipolar magnetic field configuration. This also places the pulsar beyond the conventional pulsar death line, where radio emission is expected to cease. The spin period of PSR J0250+5854 is similar to those of the high-energy-emitting magnetars and X-ray dim isolated neutron stars (XDINSs). However, the pulsar was not detected by the Swift/XRT in the energy band of 0.3-10 keV, placing a bolometric luminosity limit of $1.5 \times 10^{32}$ erg s$^{-1}$ for an assumed $N_{\rm H}=1.35\times10^{21}$ cm$^{-2}$ and a temperature of 85 eV (typical of XDINSs). We discuss the implications of the discovery for models of the pulsar death line as well as the prospect of finding more similarly long-period pulsars, including the advantages provided by LOTAAS for this., Comment: 16 pages, 9 figures, Accepted for publication in the Astrophysical Journal

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., Comment: 7 Figures, 9 Pages

We have observed and analysed the eclipses of the black widow pulsar J1810+1744 at low radio frequencies. Using LOw-Frequency ARray (LOFAR) and Westerbork Synthesis Radio Telescope observations between 2011--2015 we have measured variations in flux density, dispersion measure and scattering around eclipses. High-time-resolution, simultaneous beamformed and interferometric imaging LOFAR observations show concurrent disappearance of pulsations and total flux from the source during the eclipses, with a $3\sigma$ upper limit of 36 mJy ($<10\%$ of the pulsar's averaged out-of-eclipse flux density). The dispersion measure variations are highly asymmetric, suggesting a tail of material swept back due to orbital motion. The egress deviations are variable on timescales shorter than the 3.6 hr orbital period and are indicative of a clumpy medium. Additional pulse broadening detected during egress is typically $<20\%$ of the pulsar's spin period, showing no evidence of scattering the pulses beyond detectability in the beamformed data. The eclipses, lasting $\sim13\%$ of the orbit at 149 MHz, are shown to be frequency-dependent with total duration scaling as $\propto\nu^{-0.41\pm0.03}$. The results are discussed in the context of the physical parameters of the system, and an examination of eclipse mechanisms reveals cyclotron-synchrotron absorption as the most likely primary cause, although non-linear scattering mechanisms cannot be quantitatively ruled out. The inferred mass loss rate is a similar order-of-magnitude to the mean rate required to fully evaporate the companion in a Hubble time., Comment: 15 pages, 9 figures

Fast radio bursts (FRBs) are millisecond-duration, extragalactic radio flashes of unknown physical origin. FRB 121102, the only known repeating FRB source, has been localized to a star-forming region in a dwarf galaxy at redshift z = 0.193, and is spatially coincident with a compact, persistent radio source. The origin of the bursts, the nature of the persistent source, and the properties of the local environment are still debated. Here we present bursts that show ~100% linearly polarized emission at a very high and variable Faraday rotation measure in the source frame: RM_src = +1.46 x 10^5 rad m^-2 and +1.33 x 10^5 rad m^-2 at epochs separated by 7 months, in addition to narrow (< 30 mus) temporal structure. The large and variable rotation measure demonstrates that FRB 121102 is in an extreme and dynamic magneto-ionic environment, while the short burst durations argue for a neutron star origin. Such large rotation measures have, until now, only been observed in the vicinities of massive black holes (M_BH > 10^4 MSun). Indeed, the properties of the persistent radio source are compatible with those of a low-luminosity, accreting massive black hole. The bursts may thus come from a neutron star in such an environment. However, the observed properties may also be explainable in other models, such as a highly magnetized wind nebula or supernova remnant surrounding a young neutron star., Comment: Published in Nature: DOI: 10.1038/nature25149

Low-frequency polarisation observations of pulsars, facilitated by next-generation radio telescopes, provide powerful probes of astrophysical plasmas that span many orders of magnitude in magnetic field strength and scale: from pulsar magnetospheres to intervening magneto-ionic plasmas including the ISM and the ionosphere. Pulsar magnetospheres with teragauss field strengths can be explored through their numerous emission phenomena across multiple frequencies, the mechanism behind which remains elusive. Precise dispersion and Faraday rotation measurements towards a large number of pulsars probe the three-dimensional large-scale (and eventually small-scale) structure of the Galactic magnetic field, which plays a role in many astrophysical processes, but is not yet well understood, especially towards the Galactic halo. We describe some results and ongoing work from the Low Frequency Array (LOFAR) and the Murchison Widefield Array (MWA) radio telescopes in these areas. These and other pathfinder and precursor telescopes have reinvigorated low-frequency science and build towards the Square Kilometre Array (SKA), which will make significant advancements in studies of astrophysical magnetic fields in the next 50 years., Comment: 4 pages to appear in the proceedings of "IAUS 337: Pulsar Astrophysics - The Next 50 Years" eds: P. Weltevrede, B.B.P. Perera, L. Levin Preston & S. Sanidas

The Galactic magnetic field (GMF) plays a role in many astrophysical processes and is a significant foreground to cosmological signals, such as the Epoch of Reionization (EoR), but is not yet well understood. Dispersion and Faraday rotation measurements (DMs and RMs, respectively) towards a large number of pulsars provide an efficient method to probe the three-dimensional structure of the GMF. Low-frequency polarisation observations with large fractional bandwidth can be used to measure precise DMs and RMs. This is demonstrated by a catalogue of RMs (corrected for ionospheric Faraday rotation) from the Low Frequency Array (LOFAR), with a growing complementary catalogue in the southern hemisphere from the Murchison Widefield Array (MWA). These data further our knowledge of the three-dimensional GMF, particularly towards the Galactic halo. Recently constructed or upgraded pathfinder and precursor telescopes, such as LOFAR and the MWA, have reinvigorated low-frequency science and represent progress towards the construction of the Square Kilometre Array (SKA), which will make significant advancements in studies of astrophysical magnetic fields in the future. A key science driver for the SKA-Low is to study the EoR, for which pulsar and polarisation data can provide valuable insights in terms of Galactic foreground conditions., Comment: 6 pages, to appear in the proceedings of "IAUS 333: Peering towards Cosmic Dawn" eds: Vibor Jeli\'c & Thijs van der Hulst

The Murchison Widefield Array (MWA), and its recently-developed Voltage Capture System (VCS), facilitates extending the low-frequency range of pulsar observations at high-time and -frequency resolution in the Southern Hemisphere, providing further information about pulsars and the ISM. We present the results of an initial time-resolved census of known pulsars using the MWA. To significantly reduce the processing load, we incoherently sum the detected powers from the 128 MWA tiles, which yields ~10% of the attainable sensitivity of the coherent sum. This preserves the large field-of-view (~450 deg2 at 185 MHz), allowing multiple pulsars to be observed simultaneously. We developed a WIde-field Pulsar Pipeline (WIPP) that processes the data from each observation and automatically folds every known pulsar located within the beam. We have detected 50 pulsars to date, 6 of which are millisecond pulsars. This is consistent with our expectation, given the telescope sensitivity and the sky coverage of the processed data (~17,000 deg2). For ten pulsars, we present the lowest-frequency detections published. For a subset of the pulsars, we present multi-frequency pulse profiles by combining our data with published profiles from other telescopes. Since the MWA is a low-frequency precursor to the Square Kilometre Array (SKA), we use our census results to forecast that a survey using Phase 1 of SKA-Low (SKA1-Low) can potentially detect around 9400 pulsars.

We present techniques developed to calibrate and correct Murchison Widefield Array (MWA) low frequency (72-300 MHz) radio observations for polarimetry. The extremely wide field-of-view, excellent instantaneous (u, v)-coverage and sensitivity to degree-scale structure that the MWA provides enable instrumental calibration, removal of instrumental artefacts, and correction for ionospheric Faraday rotation through imaging techniques. With the demonstrated polarimetric capabilities of the MWA, we discuss future directions for polarimetric science at low frequencies to answer outstanding questions relating to polarised source counts, source depolarisation, pulsar science, low-mass stars, exoplanets, the nature of the interstellar and intergalactic media, and the solar environment., Comment: 18 pages, 11 figures, accepted for publication in PASA

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

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., Comment: 23 pages, 16 figures, 5 tables, accepted for publication in MNRAS

We present first results from a LOFAR census of non-recycled pulsars. The census includes almost all such pulsars known (194 sources) at declinations ${\rm Dec}> 8^\circ$ and Galactic latitudes $|{\rm Gb}| > 3^\circ$, regardless of their expected flux densities and scattering times. Each pulsar was observed for $\geq 20$ minutes in the contiguous frequency range of 110--188 MHz. Full-Stokes data were recorded. We present the dispersion measures, flux densities, and calibrated total intensity profiles for the 158 pulsars detected in the sample. The median uncertainty in census dispersion measures ($1.5 \times 10^{-3}$ pc cm$^{-3}$) is ten times smaller, on average, than in the ATNF pulsar catalogue. We combined census flux densities with those in the literature and fitted the resulting broadband spectra with single or broken power-law functions. For 48 census pulsars such fits are being published for the first time. Typically, the choice between single and broken power-laws, as well as the location of the spectral break, were highly influenced by the spectral coverage of the available flux density measurements. In particular, the inclusion of measurements below 100 MHz appears essential for investigating the low-frequency turnover in the spectra for most of the census pulsars. For several pulsars, we compared the spectral indices from different works and found the typical spread of values to be within 0.5--1.5, suggesting a prevailing underestimation of spectral index errors in the literature. The census observations yielded some unexpected individual source results, as we describe in the paper. Lastly, we will provide this unique sample of wide-band, low-frequency pulse profiles via the European Pulsar Network Database., Comment: 39 pages, 26 figures, accepted by A&A. Fixed a number of typos, rewritten DM variations section (conclusions unchanged). Arxiv pdf contains tables that will be published as electronic only by A&A

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/)., Comment: 38 pages, 11 figures, 5 tables, A&A in press, updated with editorial corrections

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., Comment: 23 pages, 19 figures. Accepted for publication in Astronomy & Astrophysics. MSSS Verification Field images and catalog data may be downloaded from http://vo.astron.nl

We report the detection of 48 millisecond pulsars (MSPs) out of 75 observed thus far using the LOFAR in the frequency range 110-188 MHz. We have also detected three MSPs out of nine observed in the frequency range 38-77 MHz. This is the largest sample of MSPs ever observed at these low frequencies, and half of the detected MSPs were observed for the first time at frequencies below 200 MHz. We present the average pulse profiles of the detected MSPs, their effective pulse widths, and flux densities and compare these with higher observing frequencies. The flux-calibrated, multifrequency LOFAR pulse profiles are publicly available via the EPN Database of Pulsar Profiles. We also present average values of dispersion measures (DM) and discuss DM and profile variations. About 35% of the MSPs show strong narrow profiles, another 25% exhibit scattered profiles, and the rest are only weakly detected. A qualitative comparison of LOFAR profiles with those at higher radio frequencies shows constant separation between profile components. Similarly, the profile widths are consistent with those observed at higher frequencies, unless scattering dominates at the lowest frequencies. This is very different from what is observed for normal pulsars and suggests a compact emission region in the MSP magnetosphere. The amplitude ratio of the profile components, on the other hand, can dramatically change towards low frequencies, often with the trailing component becoming dominant. As previously demonstrated this can be caused by aberration and retardation. This data set enables high-precision studies of pulse profile evolution with frequency, dispersion, Faraday rotation, and scattering in the interstellar medium. Characterising and correcting these systematic effects may improve pulsar-timing precision at higher observing frequencies, where pulsar timing array projects aim to directly detect gravitational waves., Comment: 27 pages, 10 figures, 5 tables, accepted to A&A

Fast Radio Bursts (FRBs), are millisecond radio signals that exhibit dispersion larger than what the Galactic electron density can account for. We have conducted a 1446 hour survey for Fast Radio Bursts (FRBs) at 145~MHz, covering a total of 4193 sq. deg on the sky. We used the UK station of the LOFAR radio telescope -- the Rawlings Array -- , accompanied for a majority of the time by the LOFAR station at Nan\c{c}ay, observing the same fields at the same frequency. Our real-time search backend, ARTEMIS, utilizes graphics processing units to search for pulses with dispersion measures up to 320 cm$^{-3}$ pc. Previous derived FRB rates from surveys around 1.4~GHz, and favoured FRB interpretations, motivated this survey, despite all previous detections occurring at higher dispersion measures. We detected no new FRBs above a signal-to-noise threshold of 10, leading to the most stringent upper limit yet on the FRB event rate at these frequencies: 29 sky$^{-1}$ day$^{-1}$ for 5~ms-duration pulses above 62~Jy. The non-detection could be due to scatter-broadening, limitations on the volume and time searched, or the shape of FRB flux density spectra. Assuming the latter and that FRBs are standard candles, the non-detection is compatible with the published FRB sky rate, if their spectra follow a power law with frequency ($\propto \nu^{\alpha}$), with $\alpha\gtrsim+0.1$, demonstrating a marked difference from pulsar spectra. Our results suggest that surveys at higher frequencies, including the low frequency component of the Square Kilometre Array, will have better chances to detect, estimate rates and understand the origin and properties of FRBs., Comment: Accepted for publication in MNRAS

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., Comment: 15 pages, 8 figures, 2 tables, accepted for publication in MNRAS

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., Comment: 27 pages, 11 figures, 2 tables, A&A accepted

The magneto-ionic structures of the interstellar medium of the Milky Way and the intergalactic medium are still poorly understood, especially at distances larger than a few kiloparsecs from the Sun. The three-dimensional (3D) structure of the Galactic magnetic field and electron density distribution may be probed through observations of radio pulsars, primarily owing to their compact nature, high velocities, and highly-polarized short-duration radio pulses. Phase 1 of the SKA, i.e. SKA1, will increase the known pulsar population by an order of magnitude, and the full SKA, i.e. SKA2, will discover pulsars in the most distant regions of our Galaxy. SKA1-VLBI will produce model-independent distances to a large number of pulsars, and wide-band polarization observations by SKA1-LOW and SKA1-MID will yield high precision dispersion measure, scattering measure, and rotation measure estimates along thousands of lines of sight. When combined, these observations will enable detailed tomography of the large-scale magneto-ionic structure of both the Galactic disk and the Galactic halo. Turbulence in the interstellar medium can be studied through the variations of these observables and the dynamic spectra of pulsar flux densities. SKA1-LOW and SKA1-MID will monitor interstellar weather and produce sensitive dynamic and secondary spectra of pulsar scintillation, which can be used to make speckle images of the ISM, study turbulence on scales between ~10^8 and ~10^13 m, and probe pulsar emission regions on scales down to $\sim$10 km. In addition, extragalactic pulsars or fast radio bursts to be discovered by SKA1 and SKA2 can be used to probe the electron density distribution and magnetic fields in the intergalactic medium beyond the Milky Way., Comment: 16 pages, 5 figures. Published in: Proceedings of Advancing Astrophysics with the Square Kilometre Array (AASKA14). 9 -13 June, 2014. Giardini Naxos, Italy. Online at http://pos.sissa.it/cgi-bin/reader/conf.cgi?confid=215, id.41

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., Comment: 18 pages, 10 figures, accepted for A&A

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., Comment: In press

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 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., Comment: 13 pages, 11 figures, accepted for publication in A&A

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., Comment: Accepted for publication in MNRAS

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., Comment: 56 pages, 34 figures, accepted for publication by A&A

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., Comment: Accepted for publication in A&A Letters; 5 pages, 4 colour figures