Search

The LOw Frequency ARray (LOFAR) has successfully measured cosmic rays for over a decade now. With its dense core of antenna fields in the Netherlands, it is an ideal tool for studying the radio emission from extensive air showers in the $10^{16}$ eV to $10^{18.5}$ eV range. Every air shower is measured with a small particle detector array and hundreds of antennas, which sets LOFAR apart from other air shower arrays. We present our current achievements and progress in reconstruction, interpolation, and software development during the final phases of measurement of LOFAR 1.0, before the LOFAR array gets a significant upgrade, including also plans for the final data release and refined analyses., Comment: Proceedings of the 10th ARENA Workshop 2024

Recent observations have revealed rare, previously unknown flashes of cosmic radio waves lasting from milliseconds to minutes, and with periodicity of minutes to an hour [1-4]. These transient radio signals must originate from sources in the Milky Way, and from coherent emission processes in astrophysical plasma. They are theorised to be produced in the extreme and highly magnetised environments around white dwarfs or neutron stars [5-8]. However, the astrophysical origin of these signals remains contested, and multiple progenitor models may be needed to explain their diverse properties. Here we present the discovery of a transient radio source, ILT J1101+5521, whose roughly minute-long pulses arrive with a periodicity of 125.5 minutes. We find that ILT J1101+5521 is an M dwarf - white dwarf binary system with an orbital period that matches the period of the radio pulses, which are observed when the two stars are in conjunction. The binary nature of ILT J1101+5521 establishes that some long-period radio transients originate from orbital motion modulating the observed emission, as opposed to an isolated rotating star. We conclude that ILT J1101+5521 is likely a polar system where magnetic interaction has synchronised the rotational and orbital periods of the white dwarf [9]. Magnetic interaction and plasma exchange between the two stars may generate the sporadic radio emission. Such mechanisms have been previously theorised [10-13], but not observationally established.

In this paper, we present rapid follow-up observations of the short GRB 201006A, consistent with being a compact binary merger, using the LOw Frequency ARray (LOFAR). We have detected a candidate 5.6$\sigma$, short, coherent radio flash at 144 MHz at 76.6 mins post-GRB with a 3$\sigma$ duration of 38 seconds. This radio flash is 27 arcsec offset from the GRB location, which has a probability of occurring by chance of $\sim$0.05% (3.8$\sigma$) when accounting for measurement uncertainties. Despite the offset, we show that the probability of finding an unrelated transient within 40 arcsec of the GRB location is $<10^{-6}$ and conclude that this is a candidate radio counterpart to GRB 201006A. We performed image plane dedispersion and the radio flash is tentatively (2.4$\sigma$) shown to be highly dispersed, allowing a distance estimate, corresponding to a redshift of $0.58\pm0.06$. The corresponding luminosity of the event at this distance is $6.7^{+6.6}_{-4.4} \times 10^{32}$ erg s$^{-1}$ Hz$^{-1}$. If associated with GRB 201006A, this emission would indicate prolonged activity from the central engine that is consistent with being a newborn, supramassive, likely highly magnetised, millisecond spin neutron star (a magnetar)., Comment: MNRAS Accepted

The composition of relativistic gamma-ray burst (GRB) jets and their emission mechanisms are still debated, and they could be matter or magnetically dominated. One way to distinguish these mechanisms arises because a Poynting flux dominated jet may produce low-frequency radio emission during the energetic prompt phase, through magnetic reconnection at the shock front. We present a search for radio emission coincident with three GRB X-ray flares with the LOw Frequency ARray (LOFAR), in a rapid response mode follow-up of long GRB 210112A (at z~2) with a 2 hour duration, where our observations began 511 seconds after the initial swift-BAT trigger. Using timesliced imaging at 120-168 MHz, we obtain upper limits at 3 sigma confidence of 42 mJy averaging over 320 second snapshot images, and 87 mJy averaging over 60 second snapshot images. LOFAR's fast response time means that all three potential radio counterparts to X-ray flares are observable after accounting for dispersion at the estimated source redshift. Furthermore, the radio pulse in the magnetic wind model was expected to be detectable at our observing frequency and flux density limits which allows us to disfavour a region of parameter space for this GRB. However, we note that stricter constraints on redshift and the fraction of energy in the magnetic field are required to further test jet characteristics across the GRB population., Comment: 13 pages, 7 figures, accepted for publication in MNRAS

Measuring the radio emission from cosmic ray particle cascades has proven to be a very efficient method to determine their properties such as the mass composition. Efficient modeling of the radio emission from air showers is crucial in order to extract the cosmic ray physics parameters from the measured radio emission. MGMR3D is a fast semi-analytic code that calculates the complete radio footprint, i.e.\ intensity, polarization, and pulse shapes, for a parametrized shower-current density and can be used in a chi-square optimization to fit a given radio data. It is many orders of magnitude faster than its Monte Carlo counterparts. We provide a detailed comparative study of MGMR3D to Monte Carlo simulations, where, with improved parametrizations, the shower maximum $\Xmax$ is found to have very strong agreement with a small dependency on the incoming zenith angle of the shower. Another interesting feature we observe with MGMR3D is sensitivity to the shape of the longitudinal profile in addition to $\Xmax$. This is achieved by probing the distinguishable radio footprint produced by a shower having a different longitudinal profile than usual. Furthermore, for the first time, we show the results of reconstructing shower parameters for LOFAR data using MGMR3D, and obtaining a $\Xmax$ resolution of 22 g/cm$^2$ and energy resolution of 19\%.

The current generation of air shower radio arrays has demonstrated that the atmospheric depth of the shower maximum Xmax can be reconstructed with high accuracy. These experiments are now contributing to mass composition studies in the energy range where a transition from galactic to extragalactic cosmic-ray sources is expected. However, we are still far away from an unambiguous interpretation of the data. Here we propose to use radio measurements to derive a new type of constraint on the mass composition, by reconstructing the shower length L. The low-frequency part of the Square Kilometer Array will have an extremely high antenna density of roughly 60.000 antennas within one square kilometer, and is the perfect site for high-resolution studies of air showers. In this contribution, we discuss the impact of being able to reconstruct L, and the unique contribution that SKA can make to cosmic-ray science., Comment: Proceedings 9th International Workshop on Acoustic and Radio EeV Neutrino Detection Activities - ARENA2022, 7-10 June 2022, Santiago de Compostela, Spain (8 pages)

Analysis of radio signals from cosmic-ray induced air showers has been shown to be a reliable method to extract shower parameters such as primary energy and depth of shower maximum. The required detailed air shower simulations take 1 to 3 days of CPU time per shower for a few hundred antennas. With nearly $60,000$ antennas envisioned to be used for air shower studies at the Square Kilometre Array (SKA), simulating all of these would come at unreasonable costs. We present an interpolation algorithm to reconstruct the full pulse time series at any position in the radio footprint, from a set of antennas simulated on a polar grid. Relying on Fourier series representations and cubic splines, it significantly improves on existing linear methods. We show that simulating about 200 antennas is sufficient for high-precision analysis in the SKA era, including e.g. interferometry which relies on accurate pulse shapes and timings. We therefore propose the interpolation algorithm and its implementation as a useful extension of radio simulation codes, to limit computational effort while retaining accuracy., Comment: 21 pages, 14 figures. Accepted for publication in JINST (Journal of Instrumentation)

We explore the possibilities of measuring the longitudinal profile of individual air showers beyond $X_{\rm max}$ when using very dense radio arrays such as SKA. The low-frequency part of the Square Kilometre Array, to be built in Australia, features an enormous antenna density of about $50,000$ antennas in the inner core region of radius 500 m, with a frequency band from 50 to 350 MHz. From CoREAS simulations, a SKA-Low antenna model plus noise contributions, and adapted LOFAR analysis scripts, we obtain a resolution in the shower maximum $X_{\rm max}$ and energy that is considerably better than at LOFAR. Already from this setup, we show that at least one additional parameter of the longitudinal profile can be measured. This would improve mass composition analysis by measuring an additional composition-dependent quantity. Moreover, it would offer an opportunity to discriminate between the different predictions of hadronic interaction models, hence contributing to hadronic physics at energy levels beyond man-made accelerators., Comment: 8 pages, 6 figures. Presented at ARENA Conference 2022, Santiago de Compostela, Spain

We present LOFAR imaging observations from the April/May 2020 active episode of magnetar SGR 1935+2154. We place the earliest radio limits on persistent emission following the low-luminosity fast radio burst FRB 200428 from the magnetar. We also perform an image-plane search for transient emission and find no radio flares during our observations. We examine post-FRB radio upper limits in the literature and find that all are consistent with the multi-wavelength afterglow predicted by the synchrotron maser shock model interpretation of FRB 200428. However, early optical observations appear to rule out the simple versions of the afterglow model with constant-density circumburst media. We show that these constraints may be mitigated by adapting the model for a wind-like environment, but only for a limited parameter range. In addition, we suggest that late-time non-thermal particle acceleration occurs within the afterglow model when the shock is no longer relativistic, which may prove vital for detecting afterglows from other Galactic FRBs. We also discuss future observing strategies for verifying either magnetospheric or maser shock FRB models via rapid radio observations of Galactic magnetars and nearby FRBs., Comment: Accepted for publication in MNRAS

We present an updated cosmic-ray mass composition analysis in the energy range $10^{16.8}$ to $10^{18.3}$ eV from 334 air showers measured with the LOFAR radio telescope, and selected for minimal bias. In this energy range, the origin of cosmic rays is expected to shift from galactic to extragalactic sources. The analysis is based on an improved method to infer the depth of maximum $X_{\rm max}$ of extensive air showers from radio measurements and air shower simulations. We show results of the average and standard deviation of $X_{\rm max}$ versus primary energy, and analyze the $X_{\rm max}$-dataset at distribution level to estimate the cosmic ray mass composition. Our approach uses an unbinned maximum likelihood analysis, making use of existing parametrizations of $X_{\rm max}$-distributions per element. The analysis has been repeated for three main models of hadronic interactions. Results are consistent with a significant light-mass fraction, at best fit $23$ to $39$ $\%$ protons plus helium, depending on the choice of hadronic interaction model. The fraction of intermediate-mass nuclei dominates. This confirms earlier results from LOFAR, with systematic uncertainties on $X_{\rm max}$ now lowered to 7 to $9$ $\mathrm{g/cm^2}$. We find agreement in mass composition compared to results from Pierre Auger Observatory, within statistical and systematic uncertainties. However, in line with earlier LOFAR results, we find a slightly lower average $X_{\rm max}$. The values are in tension with those found at Pierre Auger Observatory, but agree with results from other cosmic ray observatories based in the Northern hemisphere., Comment: 24 pages, 14 figures. Accepted for publication in Phys. Rev. D

Magnetars are a promising candidate for the origin of Fast Radio Bursts (FRBs). The detection of an extremely luminous radio burst from the Galactic magnetar SGR J1935+2154 on 2020 April 28 added credence to this hypothesis. We report on simultaneous and non-simultaneous observing campaigns using the Arecibo, Effelsberg, LOFAR, MeerKAT, MK2 and Northern Cross radio telescopes and the MeerLICHT optical telescope in the days and months after the April 28 event. We did not detect any significant single radio pulses down to fluence limits between 25 mJy ms and 18 Jy ms. Some observing epochs overlapped with times when X-ray bursts were detected. Radio images made on four days using the MeerKAT telescope revealed no point-like persistent or transient emission at the location of the magnetar. No transient or persistent optical emission was detected over seven days. Using the multi-colour MeerLICHT images combined with relations between DM, NH and reddening we constrain the distance to SGR J1935+2154, to be between 1.5 and 6.5 kpc. The upper limit is consistent with some other distance indicators and suggests that the April 28 burst is closer to two orders of magnitude less energetic than the least energetic FRBs. The lack of single-pulse radio detections shows that the single pulses detected over a range of fluences are either rare, or highly clustered, or both. It may also indicate that the magnetar lies somewhere between being radio-quiet and radio-loud in terms of its ability to produce radio emission efficiently., Comment: Accepted for publication in MNRAS. Corresponding author B. W. Stappers

The mergers of two neutron stars are typically accompanied by broad-band electromagnetic emission from either a relativistic jet or a kilonova. It has also been long predicted that coherent radio emission will occur during the merger phase or from a newly formed neutron star remnant, however this emission has not been seen to date. This paper presents the deepest limits for this emission from a neutron star merger folowing triggered LOFAR observations of the short gamma-ray burst (SGRB) 181123B, starting 4.4 minutes after the GRB occurred. During the X-ray plateau phase, a signature of ongoing energy injection, we detect no radio emission to a 3$\sigma$ limit of 153 mJy at 144 MHz (image integration time of 136 seconds), which is significantly fainter than the predicted emission from a standard neutron star. At a redshift of 1.8, this corresponds to a luminosity of $2.5 \times 10^{44}$ erg s$^{-1}$. Snapshot images were made of the radio observation on a range of timescales, targeting short duration radio flashes similar to fast radio bursts (FRBs). No emission was detected in the snapshot images at the location of GRB 181123B enabling constraints to be placed on the prompt coherent radio emission model and emission predicted to occur when a neutron star collapses to form a black hole. At the putative host redshift of 1.8 for GRB 181123B, the non detection of the prompt radio emission is two orders of magnitude lower than expected for magnetic reconnection models for prompt GRB emission and no magnetar emission is expected., Comment: MNRAS Accepted, 11 pages, 6 figures

With LOFAR we have been able to image the development of lightning flashes with meter-scale accuracy and unprecedented detail. We discuss the primary steps behind our most recent lightning mapping method. To demonstrate the capabilities of our technique we show and interpret images of the first few milliseconds of two intra-cloud flashes. In all our flashes the negative leaders propagate in the charge layer below the main negative charge. Among several interesting features we show that in about 2~ms after initiation the Primary Initial Leader triggers the formation of a multitude (more than ten) negative leaders in a rather confined area of the atmosphere. From these only one or two continue to propagate after about 30~ms to extend over kilometers horizontally while another may propagate back to the initiation point. We also show that normal negative leaders can transition into an initial-leader like state, potentially in the presence of strong electric fields. In addition, we show some initial breakdown pulses that occurred during the primary initial leader, and even during two "secondary" initial leaders that developed out of stepped leaders., Comment: Submitted to Journal of geophysics research: Atmospheres

We use the Low Frequency ARray (LOFAR) to probe the dynamics of the stepping process of negatively-charged plasma channels (negative leaders) in a lightning discharge. We observe that at each step of a leader, multiple pulses of VHF (30~--~80 MHz) radiation are emitted in short-duration bursts ($<10\ \mu$s). This is evidence for streamer formation during corona flashes that occur with each leader step, which has not been observed before in natural lightning and it could help explain X-ray emission from lightning leaders, as X-rays from laboratory leaders tend to be associated with corona flashes. Surprisingly we find that the stepping length is very similar to what was observed near the ground, however with a stepping time that is considerably larger, which as yet is not understood. These results will help to improve lightning propagation models, and eventually lightning protection models., Comment: 5 pages

The limited knowledge of atmospheric parameters like humidity, pressure, temperature, and the index of refraction has been one of the important systematic uncertainties in reconstructing the depth of the shower maximum from the radio emission of air showers. Current air shower Monte Carlo simulation codes like CORSIKA and the radio plug-in CoREAS use various averaged parameterized atmospheres. However, time-dependent and location-specific atmospheric models are needed for the cosmic ray analysis method used for LOFAR data. There, dedicated simulation sets are used for each detected cosmic ray, to take into account the actual atmospheric conditions at the time of the measurement. Using the Global Data Assimilation System (GDAS), a global atmospheric model, we have implemented time-dependent, realistic atmospheric profiles in CORSIKA and CoREAS. We have produced realistic event-specific atmospheres for all air showers measured with LOFAR, an event set spanning several years and many different weather conditions. A complete re-analysis of our data set shows that for the majority of data, our previous correction factor performed rather well; we found only a small systematic shift of 2 g/cm$^2$ in the reconstructed $X_{\rm max}$. However, under extreme weather conditions, for example, very low air pressure, the shift can be up to 15 g/cm$^2$. We provide a correction formula to determine the shift in $X_{\rm max}$ resulting from a comparison of simulations done using the US-Std atmosphere and the GDAS-based atmosphere., Comment: Accepted for publication in Astroparticle Physics. arXiv admin note: text overlap with arXiv:1911.02859

Cosmic rays are routinely measured at LOFAR, both with a dense array of antennas and with the LOFAR Radboud air shower Array (LORA) which is an array of plastic scintillators. In this paper, we present two results relating to the cosmic-ray energy scale of LOFAR. First, we present the reconstruction of cosmic-ray energy using radio and particle techniques along with a discussion of the event-by-event and absolute scale uncertainties. The resulting energies reconstructed with each method are shown to be in good agreement, and because the radio-based reconstructed energy has smaller uncertainty on an event-to-event basis, LOFAR analyses will use that technique in the future. Second, we present the radiation energy of air showers measured at LOFAR and demonstrate how radiation energy can be used to compare the energy scales of different experiments. The radiation energy scales quadratically with the electromagnetic energy in an air shower, which can in turn be related to the energy of the primary particle. Once the local magnetic field is accounted for, the radiation energy allows for a direct comparison between the LORA particle-based energy scale and that of the Pierre Auger Observatory. They are shown to agree to within (6$\pm$20)% for a radiation energy of 1 MeV, where the uncertainty on the comparison is dominated by the antenna calibrations of each experiment. This study motivates the development of a portable radio array which will be used to cross-calibrate the energy scales of different experiments using radiation energy and the same antennas, thereby significantly reducing the uncertainty on the comparison.

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

Repeating fast radio bursts (FRBs) present excellent opportunities to identify FRB progenitors and host environments, as well as decipher the underlying emission mechanism. Detailed studies of repeating FRBs might also hold clues to the origin of FRBs as a population. We aim to detect the first two repeating FRBs: FRB 121102 (R1) and FRB 180814.J0422+73 (R2), and characterise their repeat statistics. We also want to significantly improve the sky localisation of R2. We use the Westerbork Synthesis Radio Telescope to conduct extensive follow-up of these two repeating FRBs. The new phased-array feed system, Apertif, allows covering the entire sky position uncertainty of R2 with fine spatial resolution in one pointing. We characterise the energy distribution and the clustering of detected R1 bursts. We detected 30 bursts from R1. Our measurements indicate a dispersion measure of 563.5(2) pc cm$^{-3}$, suggesting a significant increase in DM over the past few years. We place an upper limit of 8% on the linear polarisation fraction of the brightest burst. We did not detect any bursts from R2. A single power-law might not fit the R1 burst energy distribution across the full energy range or widely separated detections. Our observations provide improved constraints on the clustering of R1 bursts. Our stringent upper limits on the linear polarisation fraction imply a significant depolarisation, either intrinsic to the emission mechanism or caused by the intervening medium, at 1400 MHz that is not observed at higher frequencies. The non-detection of any bursts from R2 implies either a highly clustered nature of the bursts, a steep spectral index, or a combination of both. Alternatively, R2 has turned off completely, either permanently or for an extended period of time., Comment: 11 pages, 7 figures, submitted to A&A

This is a collection of papers that have been contributed by the LOFAR Cosmic Ray Key Science Project (CRKSP) to the 36th International Cosmic Ray Conference held in Madison, Wisconsin, on July 24th to August 1st, 2019 (ICRC 2019). All papers contained here have been individually published in PoS(ICRC2019) with paper numbers 205, 362, 352, 416, and 363, in the order they appear in this collection. Minor modifications to the PoS versions have been applied where appropriate., Comment: 42 pages, 17 figures, Collection of papers presented at the ICRC2019

The nature of the central engines of gamma-ray bursts (GRBs) and the composition of their relativistic jets are still under debate. If the jets are Poynting flux dominated rather than baryon dominated, a coherent radio flare from magnetic re-connection events might be expected with the prompt gamma-ray emission. There are two competing models for the central engines of GRBs; a black hole or a newly formed milli-second magnetar. If the central engine is a magnetar it is predicted to produce coherent radio emission as persistent or flaring activity. In this paper, we present the deepest limits to date for this emission following LOFAR rapid response observations of GRB 180706A. No emission is detected to a 3$\sigma$ limit of 1.7 mJy beam$^{-1}$ at 144 MHz in a two-hour LOFAR observation starting 4.5 minutes after the gamma-ray trigger. A forced source extraction at the position of GRB 180706A provides a marginally positive (1 sigma) peak flux density of $1.1 \pm 0.9$ mJy. The data were time-sliced into different sets of snapshot durations to search for FRB like emission. No short duration emission was detected at the location of the GRB. We compare these results to theoretical models and discuss the implications of a non-detection., Comment: MNRAS accepted, 11 pages, 5 figures

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

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 LOFAR radio telescope is able to measure the radio emission from cosmic ray induced air showers with hundreds of individual antennas. This allows for precision testing of the emission mechanisms for the radio signal as well as determination of the depth of shower maximum $X_{\max}$, the shower observable most sensitive to the mass of the primary cosmic ray, to better than 20 g/cm$^2$. With a densely instrumented circular area of roughly 320 m$^2$, LOFAR is targeting for cosmic ray astrophysics in the energy range $10^{16}$ - $10^{18}$ eV. In this contribution we give an overview of the status, recent results, and future plans of cosmic ray detection with the LOFAR radio telescope., Comment: Proceedings of the 26th Extended European Cosmic Ray Symposium (ECRS), Barnaul/Belokurikha, 2018

Cosmic particles hitting Earth's moon produce radio emission via the Askaryan effect. If the resulting radio ns-pulse can be detected by radio telescopes, this technique potentially increases the available collective area for ZeV scale particles by several orders of magnitude compared to current experiments. The LOw Frequency ARray (LOFAR) is the largest radio telescope operating in the optimum frequency regime for this technique. In this contribution, we report on the status of the implementation of the lunar detection mode at LOFAR., Comment: Proceedings of the 26th Extended European Cosmic Ray Symposium (ECRS), Barnaul/Belokurikha, 2018

The LOw-Frequency ARray (LOFAR) is used to make precise measurements of radio emission from extensive air showers, yielding information about the primary cosmic ray. Interpreting the measured data requires an absolute and frequency-dependent calibration of the LOFAR system response. This is particularly important for spectral analyses, because the shape of the detected signal holds information about the shower development. We revisit the calibration of the LOFAR antennas in the range of 30 - 80 MHz. Using the Galactic emission and a detailed model of the LOFAR signal chain, we find an improved calibration that provides an absolute energy scale and allows for the study of frequency-dependent features in measured signals. With the new calibration, systematic uncertainties of 13% are reached, and comparisons of the spectral shape of calibrated data with simulations show promising agreement., Comment: 23 pages, 10 figures

While repeating fast radio bursts (FRBs) remain scarce in number, they provide a unique opportunity for follow-up observations that enhance our knowledge of their sources and potentially of the FRB population as a whole. Attaining more burst spectra could lead to a better understanding of the origin of these bright, millisecond-duration radio pulses. We therefore performed $\sim$20 hr of simultaneous observations on FRB 121102 with the Effelsberg 100-m radio telescope and the Low Frequency Array (LOFAR) to constrain the spectral behaviour of bursts from FRB 121102 at 1.4 GHz and 150 MHz. This campaign resulted in the detection of nine new bursts at 1.4 GHz but no simultaneous detections with LOFAR. Assuming that the ratio of the fluence at two frequencies scales as a power law, we placed a lower limit of $\alpha$ > -1.2 $\pm$ 0.4 on the spectral index for the fluence of the instantaneous broad band emission of FRB 121102. For the derivation of this limit, a realistic fluence detection threshold for LOFAR was determined empirically assuming a burst would be scattered as predicted by the NE2001 model. A significant variation was observed in the burst repeat rate R at L-band. During observations in September 2016, nine bursts were detected, giving R = 1.1 $\pm$ 0.4 hr$^{-1}$, while in November no bursts were detected, yielding R < 0.3 hr$^{-1}$ (95% confidence limit). This variation is consistent with earlier seen episodic emission of FRB 121102. In a blind and targeted search, no bursts were found with LOFAR at 150 MHz, resulting in a repeat rate limit of R < 0.16 hr$^{-1}$ (95% confidence limit). Burst repeat rate ratios of FRB 121102 at 3, 2, 1.4, and 0.15 GHz are consistent within the uncertainties with a flattening of its spectrum below 1 GHz., Comment: 8 pages, 2 tables, 1 figure

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

Contributions of the LOFAR Cosmic Ray Key Science Project to the 35th International Cosmic Ray Conference (ICRC 2017), Comment: compendium of 8 proceedings for the ICRC 2017

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

We observe a large fraction of circular polarization in radio emission from extensive air showers recorded during thunderstorms, much higher than in the emission from air showers measured during fair-weather circumstances. We show that the circular polarization of the air showers measured during thunderstorms can be explained by the change in the direction of the transverse current as a function of altitude induced by atmospheric electric fields. Thus by using the full set of Stokes parameters for these events, we obtain a good characterization of the electric fields in thunderclouds. We also measure a large horizontal component of the electric fields in the two events that we have analysed.

For the interpretation of measurements of radio emission from extensive air showers, an important systematic uncertainty arises from natural variations of the atmospheric refractive index $n$. At a given altitude, the refractivity $N=10^6\, (n-1)$ can have relative variations on the order of $10 \%$ depending on temperature, humidity, and air pressure. Typical corrections to be applied to $N$ are about $4\%$. Using CoREAS simulations of radio emission from air showers, we have evaluated the effect of varying $N$ on measurements of the depth of shower maximum $X_{\rm max}$. For an observation band of 30 to 80 MHz, a difference of $4 \%$ in refractivity gives rise to a systematic error in the inferred $X_{\rm max}$ between 3.5 and 11 $\mathrm{g/cm^2}$, for proton showers with zenith angles ranging from 15 to 50 degrees. At higher frequencies, from 120 to 250 MHz, the offset ranges from 10 to 22 $\mathrm{g/cm^2}$. These offsets were found to be proportional to the geometric distance to $X_{\rm max}$. We have compared the results to a simple model based on the Cherenkov angle. For the 120 to 250 MHz band, the model is in qualitative agreement with the simulations. In typical circumstances, we find a slight decrease in $X_{\rm max}$ compared to the default refractivity treatment in CoREAS. While this is within commonly treated systematic uncertainties, accounting for it explicitly improves the accuracy of $X_{\rm max}$ measurements., Comment: 13 pages, 5 figures. Accepted for publication in Astroparticle Physics

The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. Even lower fluxes of neutrinos with energies above $10^{22}$ eV are predicted in certain Grand-Unifying-Theories (GUTs) and e.g.\ models for super-heavy dark matter (SHDM). The significant increase in detector volume required to detect these particles can be achieved by searching for the nano-second radio pulses that are emitted when a particle interacts in Earth's moon with current and future radio telescopes. In this contribution we present the design of an online analysis and trigger pipeline for the detection of nano-second pulses with the LOFAR radio telescope. The most important steps of the processing pipeline are digital focusing of the antennas towards the Moon, correction of the signal for ionospheric dispersion, and synthesis of the time-domain signal from the polyphased-filtered signal in frequency domain. The implementation of the pipeline on a GPU/CPU cluster will be discussed together with the computing performance of the prototype., Comment: Proceedings of the 22nd International Conference on Computing in High Energy and Nuclear Physics (CHEP2016), USA

We report here on a novel analysis of the complete set of four Stokes parameters that uniquely determine the linear and/or circular polarization of the radio signal for an extensive air shower. The observed dependency of the circular polarization on azimuth angle and distance to the shower axis is a clear signature of the interfering contributions from two different radiation mechanisms, a main contribution due to a geomagnetically-induced transverse current and a secondary component due to the build-up of excess charge at the shower front. The data, as measured at LOFAR, agree very well with a calculation from first principles. This opens the possibility to use circular polarization as an investigative tool in the analysis of air shower structure, such as for the determination of atmospheric electric fields., Comment: Accepted for publication in Phys. Rev. D

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

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

We describe a method for spectral cleaning and timing calibration of short voltage time series data from individual radio interferometer receivers. It makes use of the phase differences in Fast Fourier Transform (FFT) spectra across antenna pairs. For strong, localized terrestrial sources these are stable over time, while being approximately uniform-random for a sum over many sources or for noise. Using only milliseconds-long datasets, the method finds the strongest interfering transmitters, a first-order solution for relative timing calibrations, and faulty data channels. No knowledge of gain response or quiescent noise levels of the receivers is required. With relatively small data volumes, this approach is suitable for use in an online system monitoring setup for interferometric arrays. We have applied the method to our cosmic-ray data collection, a collection of measurements of short pulses from extensive air showers, recorded by the LOFAR radio telescope. Per air shower, we have collected 2 ms of raw time series data for each receiver. The spectral cleaning has a calculated optimal sensitivity corresponding to a power signal-to-noise ratio of 0.08 (or -11 dB) in a spectral window of 25 kHz, for 2 ms of data in 48 antennas. This is well sufficient for our application. Timing calibration across individual antenna pairs has been performed at 0.4 ns precision; for calibration of signal clocks across stations of 48 antennas the precision is 0.1 ns. Monitoring differences in timing calibration per antenna pair over the course of the period 2011 to 2015 shows a precision of 0.08 ns, which is useful for monitoring and correcting drifts in signal path synchronizations. A cross-check method for timing calibration is presented, using a pulse transmitter carried by a drone flying over the array. Timing precision is similar, 0.3 ns., Comment: Accepted for publication in Astronomy & Astrophysics; 10 pages, 7 figures

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., Comment: 35 pages, 11 figures, updated version: Pierre Auger Observatory data ICRC 2015 added to Fig 2

The atmospheric electric fields in thunderclouds have been shown to significantly modify the intensity and polarization patterns of the radio footprint of cosmic-ray-induced extensive air showers. Simulations indicated a very non-linear dependence of the signal strength in the frequency window of 30-80 MHz on the magnitude of the atmospheric electric field. In this work we present an explanation of this dependence based on Monte-Carlo simulations, supported by arguments based on electron dynamics in air showers and expressed in terms of a simplified model. We show that by extending the frequency window to lower frequencies additional sensitivity to the atmospheric electric field is obtained.

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

The lunar Askaryan technique, which involves searching for Askaryan radio pulses from particle cascades in the outer layers of the Moon, is a method for using the lunar surface as an extremely large detector of ultra-high-energy particles. The high time resolution required to detect these pulses, which have a duration of around a nanosecond, puts this technique in a regime quite different from other forms of radio astronomy, with a unique set of associated technical challenges which have been addressed in a series of experiments by various groups. Implementing the methods and techniques developed by these groups for detecting lunar Askaryan pulses will be important for a future experiment with the Square Kilometre Array (SKA), which is expected to have sufficient sensitivity to allow the first positive detection using this technique. Key issues include correction for ionospheric dispersion, beamforming, efficient triggering, and the exclusion of spurious events from radio-frequency interference. We review the progress in each of these areas, and consider the further progress expected for future application with the SKA., Comment: Proceedings of the 34th International Cosmic Ray Conference (ICRC 2015), The Hague, The Netherlands

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