Search

The LOPES experiment was a radio interferometer built at the existing air shower array KASCADE-Grande in Karlsruhe, Germany. The last configuration of LOPES was called LOPES 3D and consisted of ten tripole antennas. Each of these antennas consisted of three crossed dipoles east-west, north-south, and vertically aligned. With this, LOPES 3D had the unique possibility to study the benefits of measurements with vertically aligned antennas in the environment of the well understood and calibrated particle detector array KASCADE-Grande. The measurements with three spatially coincident antennas allows a redundant reconstruction of the electric field vector. Several methods to exploit the redundancy were developed and tested. Furthermore, for the first time in LOPES, the background noise could be studied polarization- and direction dependent. With LOPES 3D it could be demonstrated that radio detection reaches a higher efficiency for inclined showers when including measurements with vertically aligned antennas and that the vertical component gets more important for the measurement of inclined showers. In this contribution we discuss a weighting scheme for the best combination of three redundant reconstructed electric field vectors. Furthermore, we discuss the influence of these weighting schemes on the ability to reconstruct air showers using the radio method. We show an estimate of the radio efficiency for inclined showers with focus on the benefits of measurements with vertically aligned antennas and we present the direction dependent noise in the different polarizations., Comment: Submission to the proceedings of the ARENA2014 workshop, which unfortunately were never published; thus documented here for reference

LOPES, the LOFAR prototype station, was an antenna array for cosmic-ray air showers operating from 2003 - 2013 within the KASCADE-Grande experiment. Meanwhile, the analysis is finished and the data of air-shower events measured by LOPES are available with open access in the KASCADE Cosmic Ray Data Center (KCDC). This article intends to provide a summary of the achievements, results, and lessons learned from LOPES. By digital, interferometric beamforming the detection of air showers became possible in the radio-loud environment of the Karlsruhe Institute of Technology (KIT). As a prototype experiment, LOPES tested several antenna types, array configurations and calibration techniques, and pioneered analysis methods for the reconstruction of the most important shower parameters, i.e., the arrival direction, the energy, and mass-dependent observables such as the position of the shower maximum. In addition to a review and update of previously published results, we also present new results based on end-to-end simulations including all known instrumental properties. For this, we applied the detector response to radio signals simulated with the CoREAS extension of CORSIKA, and analyzed them in the same way as measured data. Thus, we were able to study the detector performance more accurately than before, including some previously inaccessible features such as the impact of noise on the interferometric cross-correlation beam. These results led to several improvements, which are documented in this paper and can provide useful input for the design of future cosmic-ray experiments based on the digital radio-detection technique., Comment: Accepted for publication by EPJ C

Dispersive delays due to the Solar wind introduce excess noise in high-precision pulsar timing experiments, and must be removed in order to achieve the accuracy needed to detect, e.g., low-frequency gravitational waves. In current pulsar timing experiments, this delay is usually removed by approximating the electron density distribution in the Solar wind either as spherically symmetric, or with a two-phase model that describes the contributions from both high- and low-speed phases of the Solar wind. However, no dataset has previously been available to test the performance and limitations of these models over extended timescales and with sufficient sensitivity. Here we present the results of such a test with an optimal dataset of observations of pulsar J0034-0534, taken with the German stations of LOFAR. We conclude that the spherical approximation performs systematically better than the two-phase model at almost all angular distances, with a residual root-mean-square (rms) given by the two-phase model being up to 28% larger than the result obtained with the spherical approximation. Nevertheless, the spherical approximation remains insufficiently accurate in modelling the Solar-wind delay (especially within 20 degrees of angular distance from the Sun), as it leaves timing residuals with rms values that reach the equivalent of 0.3 microseconds at 1400 MHz. This is because a spherical model ignores the large daily variations in electron density observed in the Solar wind. In the short term, broadband observations or simultaneous observations at low frequencies are the most promising way forward to correct for Solar-wind induced delay variations., Comment: Accepted for publication in MNRAS. 17 pages, 13 figures

Context. High-precision pulsar-timing experiments are affected by temporal variations of the Dispersion Measure (DM), which are related to spatial variations in the interstellar electron content. Correcting for DM variations relies on the cold-plasma dispersion law which states that the dispersive delay varies with the squared inverse of the observing frequency. This may however give incorrect measurements if the probed electron content (and therefore the DM) varies with observing frequency, as is predicted theoretically. Aims. We study small-scale density variations in the ionised interstellar medium. These structures may lead to frequency-dependent DMs in pulsar signals and could inhibit the use of lower-frequency pulsar observations to correct time-variable interstellar dispersion in higher-frequency pulsar-timing data. Methods. We used high-cadence, low-frequency observations with three stations from the German LOng-Wavelength (GLOW) consortium, which are part of the LOw Frequency ARray (LOFAR). Specifically, 3.5 years of weekly observations of PSR J2219+4754 are presented. Results. We present the first detection of frequency-dependent DMs towards any interstellar object and a precise multi-year time-series of the time- and frequency-dependence of the measured DMs. The observed DM variability is significant and may be caused by extreme scattering events. Potential causes for frequency-dependent DMs are quantified and evaluated. Conclusions. We conclude that frequency-dependence of DMs has been reliably detected and is caused by small-scale (up to 10s of AUs) but steep density variations in the interstellar electron content. We find that long-term trends in DM variability equally affect DMs measured at both ends of our frequency band and hence the negative impact on long-term high-precision timing projects is expected to be limited., Comment: to be published in A&A (accepted 2019-02-06), 11 pages, 7 figures, update: A&A language editing, typos

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

Simultaneous observations of PSR B0823+26 with ESA's XMM-Newton, the Giant Metrewave Radio Telescope and international stations of the Low Frequency Array revealed synchronous X-ray/radio switching between a radio-bright (B) mode and a radio-quiet (Q) mode. During the B mode we detected PSR B0823+26 in 0.2$-$2 keV X-rays and discovered pulsed emission with a broad sinusoidal pulse, lagging the radio main pulse by 0.208 $\pm$ 0.012 in phase, with high pulsed fraction of 70$-$80%. During the Q mode PSR B0823+26 was not detected in X-rays (2 $\sigma$ upper limit a factor ~9 below the B-mode flux). The total X-ray spectrum, pulse profile and pulsed fraction can globally be reproduced with a magnetized partially ionized hydrogen atmosphere model with three emission components: a primary small hot spot ($T$$\sim$3.6$\times10^6$ K, $R$$\sim$17 m), a larger cooler concentric ring ($T$$\sim$1.1$\times10^6$ K, $R$$\sim$280 m) and an antipodal hot spot ($T$$\sim$1.1$\times10^6 $ K, $R$$\sim$100 m), for the angle between the rotation axis and line of sight direction $\sim66^\circ$. The latter is in conflict with the radio derived value of $(84\pm0.7)^\circ$. The average X-ray flux within hours-long B-mode intervals varied by a factor $\pm$20%, possibly correlated with variations in the frequency and lengths of short radio nulls or short durations of weak emission. The correlated X-ray/radio moding of PSR B0823+26 is compared with the anti-correlated moding of PSR B0943+10, and the lack of X-ray moding of PSR B1822-09. We speculate that the X-ray/radio switches of PSR B0823+26 are due to variations in the rate of accretion of material from the interstellar medium through which it is passing., Comment: 16 pages, 11 figures

The new generation of broad-band radio continuum surveys will provide large data sets with polarization information. New algorithms need to be developed to extract reliable catalogs of linearly polarized sources that can be used to characterize those sources and produce a dense rotation measure (RM) grid to probe magneto-ionized structures along the line of sight via Faraday rotation. The aim of the paper is to develop a computationally efficient and rigorously defined source-finding algorithm for linearly polarized sources. We used a calibrated data set from the LOw Frequency ARray (LOFAR) at 150 MHz centered on the nearby galaxy M51 to search for polarized background sources. With a new imaging software, we re-imaged the field at a resolution of 18''x15'' and cataloged a total of about 3000 continuum sources within 2.5 degrees of the center of M51. We made small Stokes Q and U images centered on each source brighter than 100 mJy in total intensity (201 sources) and used RM synthesis to create corresponding Faraday cubes that were analyzed individually. For each source, the noise distribution function was determined from a subset of measurements at high Faraday depths where no polarization is expected; the peaks in polarized intensity in the Faraday spectrum were identified and the p-value of each source was calculated. Finally, the false discovery rate method was applied to the list of p-values to produce a list of polarized sources and quantify the reliability of the detections. We also analyzed sources fainter than 100 mJy but that were reported as polarized in the literature at at least another radio frequency. Of the 201 sources that were searched for polarization, 6 polarized sources were detected confidently (with a false discovery rate of 5 percent). This corresponds to a number density of 1 polarized source per 3.3 square degrees, or 0.3 source/square degree. [Abridged], Comment: 17 pages, 11 figures, accepted for publication in A&A

Low-frequency radio continuum observations of edge-on galaxies are ideal to study cosmic-ray electrons (CREs) in halos via radio synchrotron emission and to measure magnetic field strengths. We obtained new observations of the edge-on spiral galaxy NGC 891 at 129-163 MHz with the LOw Frequency ARray (LOFAR) and at 13-18 GHz with the Arcminute Microkelvin Imager (AMI) and combine them with recent high-resolution Very Large Array (VLA) observations at 1-2 GHz, enabling us to study the radio continuum emission over two orders of magnitude in frequency. The spectrum of the integrated nonthermal flux density can be fitted by a power law with a spectral steepening towards higher frequencies or by a curved polynomial. Spectral flattening at low frequencies due to free-free absorption is detected in star-forming regions of the disk. The mean magnetic field strength in the halo is 7 +- 2 $\mu$G. The scale heights of the nonthermal halo emission at 146 MHz are larger than those at 1.5 GHz everywhere, with a mean ratio of 1.7 +- 0.3, indicating that spectral ageing of CREs is important and that diffusive propagation dominates. The halo scale heights at 146 MHz decrease with increasing magnetic field strengths which is a signature of dominating synchrotron losses of CREs. On the other hand, the spectral index between 146 MHz and 1.5 GHz linearly steepens from the disk to the halo, indicating that advection rather than diffusion is the dominating CRE transport process. This issue calls for refined modelling of CRE propagation., Comment: Accepted for publication in Astronomy & Astrophysics

PSR J0922+0638 (B0919+06) shows unexplained anomalous variations in the on-pulse phase, where the pulse appears to episodically move to an earlier longitude for a few tens of rotations before reverting to the usual phase for approximately several hundred to more than a thousand rotations. These events, where the pulse moves in phase by up to 5$^{\circ}$, have been previously detected in observations from $\sim$300 to 2000 MHz. We present simultaneous observations from the Effelsberg 100-m radio telescope at 1350 MHz and the Bornim (Potsdam) station of the LOw Frequency ARray at 150 MHz. Our observations present the first evidence for an absence of the anomalous phase-shifting behaviour at 150 MHz. Instead, the observed intensity at the usual pulse-phase typically decreases, often showing a pseudo-nulling feature corresponding to the times when phase shifts are observed at 1350 MHz. The presence of weak emission at the usual pulse-phase supports the theory that these shifts may result from processes similar to the 'profile-absorption' expected to operate for PSR J0814+7429 (B0809+74). A possible mechanism for this could be intrinsic variations of the emission within the pulsar's beam combined with absorption by expanding shells of electrons in the line of sight., Comment: 5 pages, 3 figures, accepted for publication in MNRAS Letters

Low-mass galaxies are subject to strong galactic outflows, in which cosmic rays may play an important role, they can be best traced with low-frequency radio continuum observations, which are less affected by spectral ageing. We present a study of the nearby star burst dwarf irregular galaxy IC 10 using observations at 140 MHz with the LOw-Frequency ARray (LOFAR), at 1580 MHz with the Very Large Array (VLA) and at 6200 MHz with the VLA and the 100-m Effelsberg telescope. We find that IC 10 has a low-frequency radio halo, which manifests itself as a second component (thick disc) in the minor axis profiles of the non-thermal radio continuum emission at 140 and 1580 MHz. These profiles are then fitted with 1D cosmic-ray transport models for pure diffusion and advection. We find that a diffusion model fits best, with a diffusion coefficient of $D=(0.4$-$0.8) \times 10^{26}(E/{\rm GeV})^{0.5}~{\rm cm^2\,s^{-1}}$, which is at least an order of magnitude smaller than estimates both from anisotropic diffusion and the diffusion length. In contrast, advection models, which cannot be ruled out due to the mild inclination, while providing poorer fits, result in advection speeds close to the escape velocity of $\approx$$50~\rm km\,s^{-1}$, as expected for a cosmic-ray driven wind. Our favoured model with an accelerating wind provides a self-consistent solution, where the magnetic field is in energy equipartition with both the warm neutral and warm ionized medium with an important contribution from cosmic rays. Consequently, cosmic rays can play a vital role for the launching of galactic winds in the disc--halo interface., Comment: 10 pages, 5 colour figures, 1 table. Accepted by MNRAS

New simultaneous X-ray and radio observations of the archetypal mode-switching pulsar PSR B0943+10 have been carried out with XMM-Newton and the LOFAR, LWA and Arecibo radio telescopes in November 2014. They allowed us to better constrain the X-ray spectral and variability properties of this pulsar and to detect, for the first time, the X-ray pulsations also during the X-ray-fainter mode. The combined timing and spectral analysis indicates that unpulsed non-thermal emission, likely of magnetospheric origin, and pulsed thermal emission from a small polar cap are present during both radio modes and vary in a correlated way., Comment: 4 pages, 1 figure, 1 table - 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 optimal source-finding strategy for linear polarization data is an unsolved problem, with many inhibitive factors imposed by the technically-challenging nature of polarization observations. Such an algorithm is essential for Square Kilometre Array (SKA) pathfinder surveys, such as the Multifrequency Snapshot Sky Survey (MSSS) with the LOw Frequency ARray (LOFAR), as data volumes are significant enough to prohibit manual inspection. We present a new strategy of `Faraday Moments' for source-finding in linear polarization with LOFAR, using the moments of the frequency-dependent full-Stokes data (i.e. the mean, standard deviation, skewness, and excess kurtosis). Through simulations of the sky, we find that moments can identify polarized sources with a high completeness: 98.5% at a signal-to-noise of 5. While the method has low reliability, Rotation Measure (RM) Synthesis can be applied per candidate source to filter out instrumental and spurious detections. This combined strategy will result in a complete and reliable catalogue of polarized sources that includes the full sensitivity of the observational bandwidth. We find that the technique can reduce the number of pixels on which RM Synthesis needs to be performed by a factor of $\approx1\times10^{5}$ for source distributions anticipated with modern radio telescopes. Through tests on LOFAR data, we find that the technique works effectively in the presence of diffuse emission. Extensions of this method are directly applicable to other upcoming radio surveys such as the POlarization Sky Survey of the Universe's Magnetism (POSSUM) with the Australia Square Kilometre Array Pathfinder (ASKAP), and the SKA itself., Comment: Accepted for publication in MNRAS

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

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

The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s)., Comment: 23 pages, 15 figures, 3 tables, accepted for publication in A&A

The radio technique is a promising method for detection of cosmic-ray air showers of energies around $100\,$PeV and higher with an array of radio antennas. Since the amplitude of the radio signal can be measured absolutely and increases with the shower energy, radio measurements can be used to determine the air-shower energy on an absolute scale. We show that calibrated measurements of radio detectors operated in coincidence with host experiments measuring air showers based on other techniques can be used for comparing the energy scales of these host experiments. Using two approaches, first via direct amplitude measurements, and second via comparison of measurements with air shower simulations, we compare the energy scales of the air-shower experiments Tunka-133 and KASCADE-Grande, using their radio extensions, Tunka-Rex and LOPES, respectively. Due to the consistent amplitude calibration for Tunka-Rex and LOPES achieved by using the same reference source, this comparison reaches an accuracy of approximately $10\,\%$ - limited by some shortcomings of LOPES, which was a prototype experiment for the digital radio technique for air showers. In particular we show that the energy scales of cosmic-ray measurements by the independently calibrated experiments KASCADE-Grande and Tunka-133 are consistent with each other on this level.

(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

We report on simultaneous X-ray and radio observations of the mode-switching pulsar PSR B0943+10 obtained with the XMM-Newton satellite and the LOFAR, LWA and Arecibo radio telescopes in November 2014. We confirm the synchronous X-ray/radio switching between a radio-bright (B) and a radio-quiet (Q) mode, in which the X-ray flux is a factor ~2.4 higher than in the B-mode. We discovered X-ray pulsations, with pulsed fraction of 38+/-5% (0.5-2 keV), during the B-mode, and confirm their presence in Q-mode, where the pulsed fraction increases with energy from ~20% up to ~65% at 2 keV. We found marginal evidence for an increase in the X-ray pulsed fraction during B-mode on a timescale of hours. The Q-mode X-ray spectrum requires a fit with a two-component model (either a power-law plus blackbody or the sum of two blackbodies), while the B-mode spectrum is well fit by a single blackbody (a single power-law is rejected). With a maximum likelihood analysis, we found that in Q-mode the pulsed emission has a thermal blackbody spectrum with temperature ~3.4x10^6 K and the unpulsed emission is a power-law with photon index ~2.5, while during B-mode both the pulsed and unpulsed emission can be fit by either a blackbody or a power law with similar values of temperature and photon index. A Chandra image shows no evidence for diffuse X-ray emission. These results support a scenario in which both unpulsed non-thermal emission, likely of magnetospheric origin, and pulsed thermal emission from a small polar cap (~1500 m^2) with a strong non-dipolar magnetic field (~10^{14} G), are present during both radio modes and vary in intensity in a correlated way. This is broadly consistent with the predictions of the partially screened gap model and does not necessarily imply global magnetospheric rearrangements to explain the mode switching., Comment: To be published on The Astrophysical Journal

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

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

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

One of the main aims of the LOPES experiment was the evaluation of the absolute amplitude of the radio signal of air showers. This is of special interest since the radio technique offers the possibility for an independent and highly precise determination of the energy scale of cosmic rays on the basis of signal predictions from Monte Carlo simulations. For the calibration of the amplitude measured by LOPES we used an external source. Previous comparisons of LOPES measurements and simulations of the radio signal amplitude predicted by CoREAS revealed a discrepancy of the order of a factor of two. A re-measurement of the reference calibration source, now performed for the free field, was recently performed by the manufacturer. The updated calibration values lead to a lowering of the reconstructed electric field measured by LOPES by a factor of $2.6 \pm 0.2$ and therefore to a significantly better agreement with CoREAS simulations. We discuss the updated calibration and its impact on the LOPES analysis results., Comment: Proceedings of the 34th ICRC, The Hague, The Netherlands, PoS(ICRC2015)311

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

We investigated the radio wavefront of cosmic-ray air showers with LOPES measurements and CoREAS simulations: the wavefront is of approximately hyperbolic shape and its steepness is sensitive to the shower maximum. For this study we used 316 events with an energy above 0.1 EeV and zenith angles below $45^\circ$ measured by the LOPES experiment. LOPES was a digital radio interferometer consisting of up to 30 antennas on an area of approximately 200 m x 200 m at an altitude of 110 m above sea level. Triggered by KASCADE-Grande, LOPES measured the radio emission between 43 and 74 MHz, and our analysis might strictly hold only for such conditions. Moreover, we used CoREAS simulations made for each event, which show much clearer results than the measurements suffering from high background. A detailed description of our result is available in our recent paper published in JCAP09(2014)025. The present proceeding contains a summary and focuses on some additional aspects, e.g., the asymmetry of the wavefront: According to the CoREAS simulations the wavefront is slightly asymmetric, but on a much weaker level than the lateral distribution of the radio amplitude., Comment: Proceedings of ARENA 2014, Anapolis, USA

LOPES was a digital antenna array detecting the radio emission of cosmic-ray air showers. The calibration of the absolute amplitude scale of the measurements was done using an external, commercial reference source, which emits a frequency comb with defined amplitudes. Recently, we obtained improved reference values by the manufacturer of the reference source, which significantly changed the absolute calibration of LOPES. We reanalyzed previously published LOPES measurements, studying the impact of the changed calibration. The main effect is an overall decrease of the LOPES amplitude scale by a factor of $2.6 \pm 0.2$, affecting all previously published values for measurements of the electric-field strength. This results in a major change in the conclusion of the paper 'Comparing LOPES measurements of air-shower radio emission with REAS 3.11 and CoREAS simulations' published in Astroparticle Physics 50-52 (2013) 76-91: With the revised calibration, LOPES measurements now are compatible with CoREAS simulations, but in tension with REAS 3.11 simulations. Since CoREAS is the latest version of the simulation code incorporating the current state of knowledge on the radio emission of air showers, this new result indicates that the absolute amplitude prediction of current simulations now is in agreement with experimental data., Comment: Partly changes results and conclusion of publication Astroparticle Physics 50-52 (2013) 76-91 (arXiv:1309.5920)

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

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

Carbon radio recombination lines (RRLs) at low frequencies (<=500 MHz) trace the cold, diffuse phase of the interstellar medium, which is otherwise difficult to observe. We present the detection of carbon RRLs in absorption in M82 with LOFAR in the frequency range of 48-64 MHz. This is the first extragalactic detection of RRLs from a species other than hydrogen, and below 1 GHz. Since the carbon RRLs are not detected individually, we cross-correlated the observed spectrum with a template spectrum of carbon RRLs to determine a radial velocity of 219 +- 9 km/s . Using this radial velocity, we stack 22 carbon-{\alpha} transitions from quantum levels n = 468-508 to achieve an 8.5 sigma detection. The absorption line profile exhibits a narrow feature with peak optical depth of 0.003 and FWHM of 31 km/s. Closer inspection suggests that the narrow feature is superimposed on a broad, shallow component. The total line profile appears to be correlated with the 21 cm H I line profile reconstructed from H I absorption in the direction of supernova remnants in the nucleus. The narrow width and centroid velocity of the feature suggests that it is associated with the nuclear starburst region. It is therefore likely that the carbon RRLs are associated with cold atomic gas in the direction of the nucleus of M82., Comment: 6 pages, 5 figures, accepted for publication in ApJL

The measurement of the radio emission from extensive air showers, induced by high-energy cosmic rays is one of the key science projects of the LOFAR radio telescope. The LOfar Radboud air shower Array (LORA) has been installed in the core of LOFAR in the Netherlands. The main purpose of LORA is to measure the properties of air showers and to trigger the read-out of the LOFAR radio antennas to register extensive air showers. The experimental set-up of the array of scintillation detectors and its performance are described., Comment: 10 pages, Accepted for publication in Nuclear Instruments and Methods A

LOPES is a digital radio interferometer located at Karlsruhe Institute of Technology (KIT), Germany, which measures radio emission from extensive air showers at MHz frequencies in coincidence with KASCADE-Grande. In this article, we explore a method (slope method) which leverages the slope of the measured radio lateral distribution to reconstruct crucial attributes of primary cosmic rays. First, we present an investigation of the method on the basis of pure simulations. Second, we directly apply the slope method to LOPES measurements. Applying the slope method to simulations, we obtain uncertainties on the reconstruction of energy and depth of shower maximum Xmax of 13% and 50 g/cm^2, respectively. Applying it to LOPES measurements, we are able to reconstruct energy and Xmax of individual events with upper limits on the precision of 20-25% for the primary energy and 95 g/cm^2 for Xmax, despite strong human-made noise at the LOPES site., Comment: 16 pages, accepted by Physical Review D

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

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%., Comment: 12 pages, 9 figures, accepted for publication in A&A

The grand-design spiral galaxy M51 was observed with the LOFAR High Frequency Antennas (HBA) and imaged in total intensity and polarisation. This observation covered the frequencies between 115 MHz and 175 MHz. We produced an image of total emission of M51 at the mean frequency of 151 MHz with 20 arcsec resolution and 0.3 mJy rms noise, which is the most sensitive image of a galaxy at frequencies below 300 MHz so far. The integrated spectrum of total radio emission is described well by a power law, while flat spectral indices in the central region indicate thermal absorption. We observe that the disk extends out to 16 kpc and see a break in the radial profile near the optical radius of the disk. Our main results, the scale lengths of the inner and outer disks at 151 MHz and 1.4 GHz, arm--interarm contrast, and the break scales of the radio--far-infrared correlations, can be explained consistently by CRE diffusion, leading to a longer propagation length of CRE of lower energy. The distribution of CRE sources drops sharply at about 10 kpc radius, where the star formation rate also decreases sharply. We find evidence that thermal absorption is primarily caused by HII regions. The non-detection of polarisation from M51 at 151 MHz is consistent with the estimates of Faraday depolarisation. Future searches for polarised emission in this frequency range should concentrate on regions with low star formation rates., Comment: 20 pages, 18 figures, accepted for publication in A&A

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

Analyzing measurements of the LOPES antenna array together with corresponding CoREAS simulations for more than 300 measured events with energy above $10^{17}\,$eV and zenith angles smaller than $45^\circ$, we find that the radio wavefront of cosmic-ray air showers is of approximately hyperbolic shape. The simulations predict a slightly steeper wavefront towards East than towards West, but this asymmetry is negligible against the measurement uncertainties of LOPES. At axis distances $\gtrsim 50\,$m, the wavefront can be approximated by a simple cone. According to the simulations, the cone angle is clearly correlated with the shower maximum. Thus, we confirm earlier predictions that arrival time measurements can be used to study the longitudinal shower development, but now using a realistic wavefront. Moreover, we show that the hyperbolic wavefront is compatible with our measurement, and we present several experimental indications that the cone angle is indeed sensitive to the shower development. Consequently, the wavefront can be used to statistically study the primary composition of ultra-high energy cosmic rays. At LOPES, the experimentally achieved precision for the shower maximum is limited by measurement uncertainties to approximately $140\,$g/cm$^2$. But the simulations indicate that under better conditions this method might yield an accuracy for the atmospheric depth of the shower maximum, $X_\mathrm{max}$, better than $30\,$g/cm$^2$. This would be competitive with the established air-fluorescence and air-Cherenkov techniques, where the radio technique offers the advantage of a significantly higher duty-cycle. Finally, the hyperbolic wavefront can be used to reconstruct the shower geometry more accurately, which potentially allows a better reconstruction of all other shower parameters, too., Comment: accepted by JCAP