Source author record

H. Falcke

H. Falcke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE astro-ph.IM hep-ex astro-ph.GA astro-ph.CO physics.ins-det astro-ph astro-ph.EP hep-ph

Catalog footprint

What is connected

162works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Measuring the muon content of inclined air showers using AERA and the water-Cherenkov detectors of the Pierre Auger Observatory

We present a novel approach for assessing the muon content of air showers with large zenith angles on a combined analysis of their radio emission and particle footprint. We use the radiation energy reconstructed by the Auger Engineering Radio Array (AERA) as an energy estimator and determine the muon number independently with the water-Cherenkov detector array of the Pierre Auger Observatory, deployed on a 1500 m grid. We focus our analysis on air showers with primary energy above 4 EeV to ensure full detection efficiency. Over approximately ten years of accumulated data, we identify a set of 40 high-quality events that are used in the analysis. The estimated muon contents in data are compatible with those for iron primaries as predicted by current-generation hadronic interaction models. This result can be interpreted as a deficit of muons in simulations as a lighter mass composition has been established from Xmax measurements. This muon deficit was already observed in previous analyses of the Auger Collaboration and is confirmed using hybrid events that include radio measurements for the first time.

preprint2022arXiv

A search for photons with energies above $2{\times}10^{17}$ eV using hybrid data from the low-energy extensions of the Pierre Auger Observatory

Ultra-high-energy photons with energies exceeding $10^{17}$ eV offer a wealth of connections to different aspects of cosmic-ray astrophysics as well as to gamma-ray and neutrino astronomy. The recent observations of photons with energies in the $10^{15}$ eV range further motivate searches for even higher-energy photons. In this paper, we present a search for photons with energies exceeding $2{\times}10^{17}$ eV using about 5.5 years of hybrid data from the low-energy extensions of the Pierre Auger Observatory. The upper limits on the integral photon flux derived here are the most stringent ones to date in the energy region between $10^{17}$ and $10^{18}$ eV.

preprint2022arXiv

Investigating Hadronic Interactions at Ultra-High Energies with the Pierre Auger Observatory

The development of an extensive air shower depends not only on the nature of the primary ultra-high-energy cosmic ray but also on the properties of the hadronic interactions. For energies above those achievable in human-made accelerators, hadronic interactions are only accessible through the studies of extensive air showers, which can be measured at the Pierre Auger Observatory. With its hybrid detector design, the Pierre Auger Observatory measures both the longitudinal development of showers in the atmosphere and the lateral distribution of particles that arrive at the ground. This way, observables that are sensitive to hadronic interactions at ultra-high energies can be obtained. While the hadronic interaction cross-section can be assessed from the longitudinal profiles, the number of muons and their fluctuations measured with the ground detectors are linked to other physical properties. In addition to these direct studies, we discuss here how measurements of the atmospheric depth of the maximum of air-shower profiles and the characteristics of the muon signal at the ground can be used to test the self-consistency of the post-LHC hadronic models.

preprint2022arXiv

LOPES 3D -- studies on the benefits of EAS-radio measurements with vertically aligned antennas

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.

preprint2022arXiv

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.

preprint2022arXiv

Testing effects of Lorentz invariance violation in the propagation of astroparticles with the Pierre Auger Observatory

Lorentz invariance violation (LIV) is often described by dispersion relations of the form $E_i^2=m_i^2+p_i^2+δ_{i,n} E^{2+n}$ with delta different based on particle type $i$, with energy $E$, momentum $p$ and rest mass $m$. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constraints on the LIV coefficients $δ_{i,n}$ tend to come from the highest energies. At sufficiently high energies, photons produced by cosmic ray interactions as they propagate through the Universe could be subluminal and unattenuated over cosmological distances. Cosmic ray interactions can also be modified and lead to detectable fingerprints in the energy spectrum and mass composition observed on Earth. The data collected at the Pierre Auger Observatory are therefore possibly sensitive to both the electromagnetic and hadronic sectors of LIV. In this article, we explore these two sectors by comparing the energy spectrum and the composition of cosmic rays and the upper limits on the photon flux from the Pierre Auger Observatory with simulations including LIV. Constraints on LIV parameters depend strongly on the mass composition of cosmic rays at the highest energies. For the electromagnetic sector, while no constraints can be obtained in the absence of protons beyond $10^{19}$ eV, we obtain $δ_{γ,0} > -10^{-21}$, $δ_{γ,1} > -10^{-40}$ eV$^{-1}$ and $δ_{γ,2} > -10^{-58}$ eV$^{-2}$ in the case of a subdominant proton component up to $10^{20}$ eV. For the hadronic sector, we study the best description of the data as a function of LIV coefficients and we derive constraints in the hadronic sector such as $δ_{\mathrm{had},0} < 10^{-19}$, $δ_{\mathrm{had},1} < 10^{-38}$ eV$^{-1}$ and $δ_{\mathrm{had},2}< 10^{-57}$ eV$^{-2}$ at 5$σ$ CL.

preprint2022arXiv

The energy spectrum of cosmic rays beyond the turn-down around $10^{17}$ eV as measured with the surface detector of the Pierre Auger Observatory

We present a measurement of the cosmic-ray spectrum above 100\,PeV using the part of the surface detector of the Pierre Auger Observatory that has a spacing of 750~m. An inflection of the spectrum is observed, confirming the presence of the so-called \emph{second-knee} feature. The spectrum is then combined with that of the 1500\,m array to produce a single measurement of the flux, linking this spectral feature with the three additional breaks at the highest energies. The combined spectrum, with an energy scale set calorimetrically via fluorescence telescopes and using a single detector type, results in the most statistically and systematically precise measurement of spectral breaks yet obtained. These measurements are critical for furthering our understanding of the highest energy cosmic rays.

preprint2021arXiv

Design, upgrade and characterization of the silicon photomultiplier front-end for the AMIGA detector at the Pierre Auger Observatory

AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to complement the study of ultra-high-energy cosmic rays (UHECR) by measuring the muon content of extensive air showers (EAS). It consists of an array of 61 water Cherenkov detectors on a denser spacing in combination with underground scintillation detectors used for muon density measurement. Each detector is composed of three scintillation modules, with 10 m$^2$ detection area per module, buried at 2.3 m depth, resulting in a total detection area of 30 m$^2$. Silicon photomultiplier sensors (SiPM) measure the amount of scintillation light generated by charged particles traversing the modules. In this paper, the design of the front-end electronics to process the signals of those SiPMs and test results from the laboratory and from the Pierre Auger Observatory are described. Compared to our previous prototype, the new electronics shows a higher performance, higher efficiency and lower power consumption, and it has a new acquisition system with increased dynamic range that allows measurements closer to the shower core. The new acquisition system is based on the measurement of the total charge signal that the muonic component of the cosmic ray shower generates in the detector.

preprint2021arXiv

Final results of the LOPES radio interferometer for cosmic-ray air showers

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.

preprint2021arXiv

Sub-arcsecond imaging with the International LOFAR Telescope I. Foundational calibration strategy and pipeline

[abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $μ$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $μ$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit.

preprint2021arXiv

Sub-arcsecond imaging with the International LOFAR Telescope: II. Completion of the LOFAR Long-Baseline Calibrator Survey

The Low-Frequency Array (LOFAR) Long-Baseline Calibrator Survey (LBCS) was conducted between 2014 and 2019 in order to obtain a set of suitable calibrators for the LOFAR array. In this paper we present the complete survey, building on the preliminary analysis published in 2016 which covered approximately half the survey area. The final catalogue consists of 30006 observations of 24713 sources in the northern sky, selected for a combination of high low-frequency radio flux density and flat spectral index using existing surveys (WENSS, NVSS, VLSS, and MSSS). Approximately one calibrator per square degree, suitable for calibration of $\geq$ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 degrees and away from the Galactic plane, with a considerably lower density south of this point due to relative difficulty in selecting flat-spectrum candidate sources in this area of the sky. Use of the VLBA calibrator list, together with statistical arguments by comparison with flux densities from lower-resolution catalogues, allow us to establish a rough flux density scale for the LBCS observations, so that LBCS statistics can be used to estimate compact flux densities on scales between 300 mas and 2 arcsec, for sources observed in the survey. The LBCS can be used to assess the structures of point sources in lower-resolution surveys, with significant reductions in the degree of coherence in these sources on scales between 2 arcsec and 300 mas. The LBCS survey sources show a greater incidence of compact flux density in quasars than in radio galaxies, consistent with unified schemes of radio sources. Comparison with samples of sources from interplanetary scintillation (IPS) studies with the Murchison Widefield Array (MWA) shows consistent patterns of detection of compact structure in sources observed both interferometrically with LOFAR and using IPS.

preprint2020arXiv

Cosmic-ray anisotropies in right ascension measured by the Pierre Auger Observatory

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$σ$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $α_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

preprint2020arXiv

ESA Voyage 2050 white paper -- GrailQuest: hunting for Atoms of Space and Time hidden in the wrinkle of Space-Time

GrailQuest (Gamma Ray Astronomy International Laboratory for QUantum Exploration of Space-Time) is a mission concept based on a constellation (hundreds/thousands) of nano/micro/small-satellites in low (or near) Earth orbits. Each satellite hosts a non-collimated array of scintillator crystals coupled with Silicon Drift Detectors with broad energy band coverage (keV-MeV range) and excellent temporal resolution ( below or equal 100 nanoseconds) each with effective area around 100 cm2. This simple and robust design allows for mass-production of the satellites of the fleet. This revolutionary approach implies a huge reduction of costs, flexibility in the segmented launching strategy, and an incremental long-term plan to increase the number of detectors and their performance: a living observatory for next-generation, space-based astronomical facilities. GrailQuest is conceived as an all-sky monitor for fast localisation of high signal-to-noise ratio transients in the X/gamma-ray band, e.g. the elusive electromagnetic counterparts of gravitational wave events. Robust temporal triangulation techniques will allow unprecedented localisation capabilities, in the keV-MeV band, of a few arcseconds or below, depending on the temporal structure of the transient event. The ambitious ultimate goal of this mission is to perform the first experiment, in quantum gravity, to directly probe space-time structure down to the minuscule Planck scale, by constraining or measuring a first order dispersion relation for light in vacuo. This is obtained by detecting delays between photons of different energies in the prompt emission of Gamma-ray Bursts.

preprint2020arXiv

Radio emission from negative lightning leader steps reveals inner meter-scale structure

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\ μ$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.

preprint2020arXiv

Reconstructing air shower parameters with LOFAR using event specific GDAS atmospheres

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.

preprint2020arXiv

Search for magnetically-induced signatures in the arrival directions of ultra-high-energy cosmic rays measured at the Pierre Auger Observatory

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

preprint2020arXiv

Studies on the response of a water-Cherenkov detector of the Pierre Auger Observatory to atmospheric muons using an RPC hodoscope

Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

preprint2019arXiv

A Roadmap towards a Space-based Radio Telescope for Ultra-Low Frequency Radio Astronomy

The past two decades saw a renewed interest in low frequency radio astronomy, with a particular focus on frequencies above 30 MHz. However, at frequencies below 30 MHz, Earth-based observations are limited due to a combination of severe ionospheric distortions, almost full reflection of radio waves below 10 MHz, solar eruptions and human-made radio frequency interference (RFI). A space or Lunar-based ultra-low-frequency (or ultra-long-wavelength, ULW) array would suffer significantly less from these limitations and hence would open up the last, virtually unexplored frequency domain in the electromagnetic spectrum. A roadmap has been initiated in order to explore the opportunity of building a swarm of satellites to observe at the frequency band below 30 MHz. This roadmap, dubbed Orbiting Low Frequency Antennas for Radio Astronomy (OLFAR), presents a space-based ultra-low frequency radio telescope that will explore the Universe's so-called dark ages, map the interstellar medium, and study planetary and solar bursts in the solar system and search them in other planetary systems. Such a system will comprise of a swarm of hundreds to thousands of satellites, working together as a single aperture synthesis instrument deployed sufficiently far away from Earth to avoid terrestrial RFI. A number of key technologies of OLFAR are still to be developed and proven. The first step in this roadmap is the NCLE (Netherlands China Low Frequency Explorer) experiment launched in May 2018 on the Chinese Chang'e-4 mission. The NCLE payload consists of a three monopole antenna system from which the first data stream is expected in the second half of 2019, which will provide important feedback for future science and technology opportunities. In this paper, the roadmap towards OLFAR, a brief overview of the science opportunities, and the technological and programmatic challenges of the mission are presented.

preprint2019arXiv

rPICARD: A CASA-based Calibration Pipeline for VLBI Data

Currently, HOPS and AIPS are the primary choices for the time-consuming process of (millimeter) Very Long Baseline Interferometry (VLBI) data calibration. However, for a full end-to-end pipeline, they either lack the ability to perform easily scriptable incremental calibration or do not provide full control over the workflow with the ability to manipulate and edit calibration solutions directly. The Common Astronomy Software Application (CASA) offers all these abilities, together with a secure development future and an intuitive Python interface, which is very attractive for young radio astronomers. Inspired by the recent addition of a global fringe-fitter, the capability to convert FITS-IDI files to measurement sets, and amplitude calibration routines based on ANTAB metadata, we have developed the the CASA-based Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD). The pipeline will be able to handle data from multiple arrays: EHT, GMVA, VLBA and the EVN in the first release. Polarization and phase-referencing calibration are supported and a spectral line mode will be added in the future. The large bandwidths of future radio observatories ask for a scalable reduction software. Within CASA, a message passing interface (MPI) implementation is used for parallelization, reducing the total time needed for processing. The most significant gain is obtained for the time-consuming fringe-fitting task where each scan be processed in parallel.

preprint2017arXiv

Search for High-energy Neutrinos from Binary Neutron Star Merger GW170817 with ANTARES, IceCube, and the Pierre Auger Observatory

The Advanced LIGO and Advanced Virgo observatories recently discovered gravitational waves from a binary neutron star inspiral. A short gamma-ray burst (GRB) that followed the merger of this binary was also recorded by the Fermi Gamma-ray Burst Monitor (Fermi-GBM), and the Anticoincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory (INTEGRAL), indicating particle acceleration by the source. The precise location of the event was determined by optical detections of emission following the merger. We searched for high-energy neutrinos from the merger in the GeV--EeV energy range using the ANTARES, IceCube, and Pierre Auger Observatories. No neutrinos directionally coincident with the source were detected within $\pm500$ s around the merger time. Additionally, no MeV neutrino burst signal was detected coincident with the merger. We further carried out an extended search in the direction of the source for high-energy neutrinos within the 14-day period following the merger, but found no evidence of emission. We used these results to probe dissipation mechanisms in relativistic outflows driven by the binary neutron star merger. The non-detection is consistent with model predictions of short GRBs observed at a large off-axis angle.

preprint2016arXiv

A large light-mass component of cosmic rays at 10^{17} - 10^{17.5} eV from radio observations

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.

preprint2016arXiv

A LOFAR census of non-recycled pulsars: average profiles, dispersion measures, flux densities, and spectra

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.

preprint2016arXiv

Azimuthal asymmetry in the risetime of the surface detector signals of the Pierre Auger Observatory

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec θ)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec θ)_\mathrm{max}$.

preprint2016arXiv

Cosmic-ray energy spectrum and composition up to the ankle - the case for a second Galactic component

We have carried out a detailed study to understand the observed energy spectrum and composition of cosmic rays with energies up to ~10^18 eV. Our study shows that a single Galactic component with subsequent energy cut-offs in the individual spectra of different elements, optimised to explain the observed spectra below ~10^14 eV and the knee in the all-particle spectrum, cannot explain the observed all-particle spectrum above ~2x10^16 eV. We discuss two approaches for a second component of Galactic cosmic rays -- re-acceleration at a Galactic wind termination shock, and supernova explosions of Wolf-Rayet stars, and show that the latter scenario can explain almost all observed features in the all-particle spectrum and the composition up to ~10^18 eV, when combined with a canonical extra-galactic spectrum expected from strong radio galaxies or a source population with similar cosmological evolution. In this two-component Galactic model, the knee at ~ 3x10^15 eV and the second knee at ~10^17 eV in the all-particle spectrum are due to the cut-offs in the first and second components, respectively. We also discuss several variations of the extra-galactic component, from a minimal contribution to scenarios with a significant component below the ankle (at ~4x10^18 eV), and find that extra-galactic contributions in excess of regular source evolution are neither indicated nor in conflict with the existing data. Our main result is that the second Galactic component predicts a composition of Galactic cosmic rays at and above the second knee that largely consists of helium or a mixture of helium and CNO nuclei, with a weak or essentially vanishing iron fraction, in contrast to most common assumptions. This prediction is in agreement with new measurements from LOFAR and the Pierre Auger Observatory which indicate a strong light component and a rather low iron fraction between ~10^17 and 10^18 eV.

preprint2016arXiv

Evidence for a mixed mass composition at the `ankle' in the cosmic-ray spectrum

We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the electromagnetic and muonic components. It allows an accurate determination of the spread of primary masses in the cosmic-ray flux. Up till now, constraints on the spread of primary masses have been dominated by systematic uncertainties. The present correlation measurement is not affected by systematics in the measurement of the depth of shower maximum or the signal in the water Cherenkov stations. The analysis relies on general characteristics of air showers and is thus robust also with respect to uncertainties in hadronic event generators. The observed correlation in the energy range around the `ankle' at $\lg(E/{\rm eV})=18.5-19.0$ differs significantly from expectations for pure primary cosmic-ray compositions. A light composition made up of proton and helium only is equally inconsistent with observations. The data are explained well by a mixed composition including nuclei with mass $A > 4$. Scenarios such as the proton dip model, with almost pure compositions, are thus disfavoured as the sole explanation of the ultrahigh-energy cosmic-ray flux at Earth.

preprint2016arXiv

LBCS: the LOFAR Long-Baseline Calibrator Survey

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

preprint2016arXiv

LOFAR 150-MHz observations of the Boötes field: Catalogue and Source Counts

We present the first wide area (19 deg$^2$), deep ($\approx120-150$ μJy beam$^{-1}$), high resolution ($5.6 \times 7.4$ arcsec) LOFAR High Band Antenna image of the Boötes field made at 130-169 MHz. This image is at least an order of magnitude deeper and 3-5 times higher in angular resolution than previously achieved for this field at low frequencies. The observations and data reduction, which includes full direction-dependent calibration, are described here. We present a radio source catalogue containing 6276 sources detected over an area of $19$\,deg$^2$, with a peak flux density threshold of $5σ$. As the first thorough test of the facet calibration strategy, introduced by van Weeren et al., we investigate the flux and positional accuracy of the catalogue. We present differential source counts that reach an order of magnitude deeper in flux density than previously achieved at these low frequencies, and show flattening at 150 MHz flux densities below 10 mJy associated with the rise of the low flux density star-forming galaxies and radio-quiet AGN.

preprint2016arXiv

Low-radio-frequency eclipses of the redback pulsar J2215+5135 observed in the image plane with LOFAR

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 ν^{-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.

preprint2016arXiv

Measurement of the circular polarization in radio emission from extensive air showers confirms emission mechanisms

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.

preprint2016arXiv

Nanosecond-level time synchronization of autonomous radio detector stations for extensive air showers

To exploit the full potential of radio measurements of cosmic-ray air showers at MHz frequencies, a detector timing synchronization within 1 ns is needed. Large distributed radio detector arrays such as the Auger Engineering Radio Array (AERA) rely on timing via the Global Positioning System (GPS) for the synchronization of individual detector station clocks. Unfortunately, GPS timing is expected to have an accuracy no better than about 5 ns. In practice, in particular in AERA, the GPS clocks exhibit drifts on the order of tens of ns. We developed a technique to correct for the GPS drifts, and an independent method is used for cross-checks that indeed we reach nanosecond-scale timing accuracy by this correction. First, we operate a "beacon transmitter" which emits defined sine waves detected by AERA antennas recorded within the physics data. The relative phasing of these sine waves can be used to correct for GPS clock drifts. In addition to this, we observe radio pulses emitted by commercial airplanes, the position of which we determine in real time from Automatic Dependent Surveillance Broadcasts intercepted with a software-defined radio. From the known source location and the measured arrival times of the pulses we determine relative timing offsets between radio detector stations. We demonstrate with a combined analysis that the two methods give a consistent timing calibration with an accuracy of 2 ns or better. Consequently, the beacon method alone can be used in the future to continuously determine and correct for GPS clock drifts in each individual event measured by AERA.

preprint2016arXiv

Probing spacetime around Sagittarius A* using modeled VLBI closure phases

The emission region and black hole shadow of Sagittarius A*, the supermassive black hole at the Galactic Center, can be probed with millimeter Very Long Baseline Interferometry. Our goal is to probe the geometry of the emitting plasma around Sgr A* by using modeled mm-VLBI closure phase calculations at 1.3 mm and to constrain the observer's inclination angle and position angle of the black hole spin axis. We have simulated images for three different models of the emission of Sgr A*: an orbiting spot, a disk model, and a jet model. The orbiting spot model was used as a test case scenario, while the disk and jet models are physically driven scenarios based on standard three-dimensional general relativistic magnetohydrodynamic simulations of hot accretion flows. Our results are compared to currently available closure phase observational limits. Our results indicate that more models with closer to edge-on viewing angles are consistent with observational limits. In general, jet and disk geometries can reproduce similar closure phases for different sets of viewing and position angles. Consequently, the favored black hole spin orientation and its magnitude are strongly model dependent. We find that both the jet and the disk models can explain current VLBI limits. We conclude that new observations at 1.3 mm and possibly at longer wavelengths including other triangles of VLBI baselines are necessary to interpret Sgr A* emission and the putative black hole spin parameters.

preprint2016arXiv

Prototype muon detectors for the AMIGA component of the Pierre Auger Observatory

Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to extend its range of detection and to directly measure the muon content of the particle showers. It consists of an infill of surface water-Cherenkov detectors accompanied by buried scintillator detectors used for muon counting. The main objectives of the AMIGA engineering array, referred to as the Unitary Cell, are to identify and resolve all engineering issues as well as to understand the muon-number counting uncertainties related to the design of the detector. The mechanical design, fabrication and deployment processes of the muon counters of the Unitary Cell are described in this document. These muon counters modules comprise sealed PVC casings containing plastic scintillation bars, wavelength-shifter optical fibers, 64 pixel photomultiplier tubes, and acquisition electronics. The modules are buried approximately 2.25 m below ground level in order to minimize contamination from electromagnetic shower particles. The mechanical setup, which allows access to the electronics for maintenance, is also described in addition to tests of the modules' response and integrity. The completed Unitary Cell has measured a number of air showers of which a first analysis of a sample event is included here.

preprint2016arXiv

Search for correlations between the arrival directions of IceCube neutrino events and ultrahigh-energy cosmic rays detected by the Pierre Auger Observatory and the Telescope Array

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

preprint2016arXiv

Search for ultrarelativistic magnetic monopoles with the Pierre Auger Observatory

We present a search for ultra-relativistic magnetic monopoles with the Pierre Auger Observatory. Such particles, possibly a relic of phase transitions in the early universe, would deposit a large amount of energy along their path through the atmosphere, comparable to that of ultrahigh-energy cosmic rays (UHECRs). The air shower profile of a magnetic monopole can be effectively distinguished by the fluorescence detector from that of standard UHECRs. No candidate was found in the data collected between 2004 and 2012, with an expected background of less than 0.1 event from UHECRs. The corresponding 90% confidence level (C.L.) upper limits on the flux of ultra-relativistic magnetic monopoles range from $10^{-19}$ (cm$^{2}$ sr s)$^{-1}$ for a Lorentz factor $γ=10^9$ to $2.5 \times10^{-21}$ (cm$^{2}$ sr s)$^{-1}$ for $γ=10^{12}$. These results - the first obtained with a UHECR detector - improve previously published limits by up to an order of magnitude.

preprint2016arXiv

Testing Hadronic Interactions at Ultrahigh Energies with Air Showers Measured by the Pierre Auger Observatory

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

preprint2016arXiv

The MHz-peaked radio spectrum of the unusual gamma-ray source PMN J1603-4904

Context. The majority of bright extragalactic gamma-ray sources are blazars. Only a few radio galaxies have been detected by Fermi/LAT. Recently, the GHz-peaked spectrum source PKS 1718-649 was confirmed to be gamma-ray bright, providing further evidence for the existence of a population of gamma-ray loud, compact radio galaxies. A spectral turnover in the radio spectrum in the MHz to GHz range is a characteristic feature of these objects, which are thought to be young due to their small linear sizes. The multiwavelength properties of the gamma-ray source PMN J1603-4904 suggest that it is a member of this source class. Aims. The known radio spectrum of PMN J1603-4904 can be described by a power law above 1 GHz. Using observations from the Giant Metrewave Radio Telescope (GMRT) at 150, 325, and 610 MHz, we investigate the behaviour of the spectrum at lower frequencies to search for a low-frequency turnover. Methods. Data from the TIFR GMRT Sky Survey (TGSS ADR) catalogue and archival GMRT observations were used to construct the first MHz to GHz spectrum of PMN J1603-4904. Results. We detect a low-frequency turnover of the spectrum and measure the peak position at about 490 MHz (rest-frame), which, using the known relation of peak frequency and linear size, translates into a maximum linear source size of ~1.4 kpc. Conclusions. The detection of the MHz peak indicates that PMN J1603-4904 is part of this population of radio galaxies with turnover frequencies in the MHz to GHz regime. Therefore it can be considered the second, confirmed object of this kind detected in gamma-rays. Establishing this gamma-ray source class will help to investigate the gamma-ray production sites and to test broadband emission models.

preprint2016arXiv

The Pierre Auger Observatory Upgrade - Preliminary Design Report

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

preprint2016arXiv

The radio evolution of supernova SN$\,$2008iz in M$\,$82

We report on multi-frequency Very Large Array (VLA) and Very Long Baseline Interferometry (VLBI) radio observations for a monitoring campaign of supernova SN$\,$2008iz in the nearby irregular galaxy M$\,$82. We fit two models to the data, a simple time power-law, S$\propto$ t$^β$, and a simplified Weiler model, yielding decline indices, $β$ of -1.22$\pm$0.07 (days 100-1500) and -1.41$\pm$0.02 (days 76-2167), respectively. The late-time radio light curve evolution shows flux-density flares at $\sim\,$970 and $\sim\,$1400 days which are a factor of $\sim\,$2 and $\sim\,$4 higher than the expected flux, respectively. The later flare, besides being brighter, does not show signs of decline at least from results examined so far (2014 January 23, day 2167). We derive the spectral index, $α$, S$\propto$ $ν^α$ for frequencies 1.4 to 43$\,$GHz for SN$\,$2008iz during the period from $\sim\,$430 to 2167 days after the supernova explosion. The value of $α$ shows no signs of evolution and remains steep $\approx$$-$1 throughout the period, unlike that of SN$\,$1993J which started flattening at $\sim\,$day 970. From the 4.8 and 8.4$\,$GHz VLBI images, the supernova expansion is seen to start with shell like structure that gets more and more asymmetric, then breaks up in the later epochs with bright structures dominating the southern part of the ring. This structural evolution differs significantly from SN$\,$1993J, which remains circularly symmetric over 4000 days after the explosion. The VLBI 4.8 and 8.4$\,$GHz images are used to derive a deceleration index, $m$, for SN$\,$2008iz, of 0.86$\,\pm\,$0.02, and the average expansion velocity between days 73 and 1400 as (12.1$\pm$0.2)$\times$10$^{3}\,$km$\,$s$^{-1}$.

preprint2016arXiv

Timing calibration and spectral cleaning of LOFAR time series data

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.

preprint2015arXiv

A LOFAR census of millisecond pulsars

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

preprint2015arXiv

Calibrating the absolute amplitude scale for air showers measured at LOFAR

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.

preprint2015arXiv

Comparing LOPES measurements of air-shower radio emission with REAS 3.11 and CoREAS simulations

Cosmic ray air showers emit radio pulses at MHz frequencies, which can be measured with radio antenna arrays - like LOPES at the Karlsruhe Institute of Technology in Germany. To improve the understanding of the radio emission, we test theoretical descriptions with measured data. The observables used for these tests are the absolute amplitude of the radio signal, and the shape of the radio lateral distribution. We compare lateral distributions of more than 500 LOPES events with two recent and public Monte Carlo simulation codes, REAS 3.11 and CoREAS (v 1.0). The absolute radio amplitudes predicted by REAS 3.11 are in good agreement with the LOPES measurements. The amplitudes predicted by CoREAS are lower by a factor of two, and marginally compatible with the LOPES measurements within the systematic scale uncertainties. In contrast to any previous versions of REAS, REAS 3.11 and CoREAS now reproduce the shape of the measured lateral distributions correctly. This reflects a remarkable progress compared to the situation a few years ago, and it seems that the main processes for the radio emission of air showers are now understood: The emission is mainly due to the geomagnetic deflection of the electrons and positrons in the shower. Less important but not negligible is the Askaryan effect (net charge variation). Moreover, we confirm that the refractive index of the air plays an important role, since it changes the coherence conditions for the emission: Only the new simulations including the refractive index can reproduce rising lateral distributions which we observe in a few LOPES events. Finally, we show that the lateral distribution is sensitive to the energy and the mass of the primary cosmic ray particles.

preprint2015arXiv

High-precision measurements of extensive air showers with the SKA

As of 2023, the Square Kilometre Array will constitute the world's largest radio telescope, offering unprecedented capabilities for a diverse science programme in radio astronomy. At the same time, the SKA will be ideally suited to detect extensive air showers initiated by cosmic rays in the Earth's atmosphere via their radio emission. With its very dense and uniform antenna spacing in a fiducial area of one km$^2$ and its large bandwidth of 50-350 MHz, the low-frequency part of the SKA will provide very precise measurements of individual cosmic ray air showers. These precision measurements will allow detailed studies of the mass composition of cosmic rays in the energy region of transition from a Galactic to an extragalactic origin. Also, the SKA will facilitate three-dimensional "tomography" of the electromagnetic cascades of air showers, allowing the study of particle interactions at energies beyond the reach of the LHC. Finally, studies of possible connections between air showers and lightning initiation can be taken to a new level with the SKA. We discuss the science potential of air shower detection with the SKA and report on the technical requirements and project status.

preprint2015arXiv

Imaging Jupiter's radiation belts down to 127 MHz with LOFAR

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

preprint2015arXiv

Improved absolute calibration of LOPES measurements and its impact on the comparison with REAS 3.11 and CoREAS simulations

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.

preprint2015arXiv

Influence of Atmospheric Electric Fields on the Radio Emission from Extensive Air Showers

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.

preprint2015arXiv

Investigation of the radio wavefront of air showers with LOPES measurements and CoREAS simulations (ARENA 2014)

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.

preprint2015arXiv

Limits on Fast Radio Bursts at 145 MHz with ARTEMIS, a real-time software backend

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

preprint2015arXiv

LOFAR discovery of a quiet emission mode in PSR B0823+26

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

preprint2015arXiv

LOFAR MSSS: Detection of a low-frequency radio transient in 400 hrs of monitoring of the North Celestial Pole

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.

preprint2015arXiv

LOFAR Sparse Image Reconstruction

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

preprint2015arXiv

LOFAR tied-array imaging and spectroscopy of solar S bursts

Context. The Sun is an active source of radio emission that is often associated with energetic phenomena ranging from nanoflares to coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), numerous millisecond duration radio bursts have been reported, such as radio spikes or solar S bursts (where S stands for short). To date, these have neither been studied extensively nor imaged because of the instrumental limitations of previous radio telescopes. Aims. Here, Low Frequency Array (LOFAR) observations were used to study the spectral and spatial characteristics of a multitude of S bursts, as well as their origin and possible emission mechanisms. Methods. We used 170 simultaneous tied-array beams for spectroscopy and imaging of S bursts. Since S bursts have short timescales and fine frequency structures, high cadence (~50 ms) tied-array images were used instead of standard interferometric imaging, that is currently limited to one image per second. Results. On 9 July 2013, over 3000 S bursts were observed over a time period of ~8 hours. S bursts were found to appear as groups of short-lived (<1 s) and narrow-bandwidth (~2.5 MHz) features, the majority drifting at ~3.5 MHz/s and a wide range of circular polarisation degrees (2-8 times more polarised than the accompanying Type III bursts). Extrapolation of the photospheric magnetic field using the potential field source surface (PFSS) model suggests that S bursts are associated with a trans-equatorial loop system that connects an active region in the southern hemisphere to a bipolar region of plage in the northern hemisphere. Conclusions. We have identified polarised, short-lived solar radio bursts that have never been imaged before. They are observed at a height and frequency range where plasma emission is the dominant emission mechanism, however they possess some of the characteristics of electron-cyclotron maser emission.

preprint2015arXiv

Measurement of the cosmic-ray energy spectrum above $10^{16}$ eV with the LOFAR Radboud Air Shower Array

The energy reconstruction of extensive air showers measured with the LOFAR Radboud Air Shower Array (LORA) is presented in detail. LORA is a particle detector array located in the center of the LOFAR radio telescope in the Netherlands. The aim of this work is to provide an accurate and independent energy measurement for the air showers measured through their radio signal with the LOFAR antennas. The energy reconstruction is performed using a parameterized relation between the measured shower size and the cosmic-ray energy obtained from air shower simulations. In order to illustrate the capabilities of LORA, the all-particle cosmic-ray energy spectrum has been reconstructed, assuming that cosmic rays are composed only of protons or iron nuclei in the energy range between $\sim2\times10^{16}$ and $2\times10^{18}$ eV. The results are compatible with literature values and a changing mass composition in the transition region from a galactic to an extragalactic origin of cosmic rays.

preprint2015arXiv

Orbital and superorbital variability of LS I +61 303 at low radio frequencies with GMRT and LOFAR

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

preprint2015arXiv

Pierre Auger Observatory and Telescope Array: Joint Contributions to the 34th International Cosmic Ray Conference (ICRC 2015)

Joint contributions of the Pierre Auger Collaboration and the Telescope Array Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

preprint2015arXiv

Precision measurements of cosmic ray air showers with the SKA

Supplemented with suitable buffering techniques, the low-frequency part of the SKA can be used as an ultra-precise detector for cosmic-ray air showers at very high energies. This would enable a wealth of scientific applications: the physics of the transition from Galactic to extragalactic cosmic rays could be probed with very high precision mass measurements, hadronic interactions could be studied up to energies well beyond the reach of man-made particle accelerators, air shower tomography could be performed with very high spatial resolution exploiting the large instantaneous bandwidth and very uniform instantaneous $u$-$v$ coverage of SKA1-LOW, and the physics of thunderstorms and possible connections between cosmic rays and lightning initiation could be studied in unprecedented levels of detail. In this article, we describe the potential of the SKA as an air shower radio detector from the perspective of existing radio detection efforts and discuss the associated technical requirements.

preprint2015arXiv

Probing Atmospheric Electric Fields in Thunderstorms through Radio Emission from Cosmic-Ray-Induced Air Showers

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

preprint2015arXiv

Pulsar polarisation below 200 MHz: Average profiles and propagation effects

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

preprint2015arXiv

Revised absolute amplitude calibration of the LOPES experiment

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.

preprint2015arXiv

The evolution of a jet ejection of the ultraluminous X-ray source Holmberg II X-1

We present quasi-simultaneous, multi-epoch radio and X-ray measurements of Holmberg II X-1 using the European VLBI Network (EVN), the Karl G. Jansky Very Large Array (VLA), and the Chandra and Swift X-ray telescopes. The X-ray data show apparently hard spectra with steady X-ray luminosities 4 months apart from each other. In the high-resolution EVN radio observations, we have detected an extended milli-arcsecond scale source with unboosted radio emission. The source emits non-thermal, likely optically thin synchrotron emission and its morphology is consistent with a jet ejection. The 9-GHz VLA data show an arcsecond-scale triple structure of Holmberg II X-1 similar to that seen at lower frequencies. However, we find that the central ejection has faded by at least a factor of 7.3 over 1.5 years. We estimate the dynamical age of the ejection to be higher than 2.1 years. We show that such a rapid cooling can be explained with simple adiabatic expansion losses. These properties of Holmberg II X-1 imply that ULX radio bubbles may be inflated by ejecta instead of self-absorbed compact jets.

preprint2015arXiv

The IceCube Neutrino Observatory, the Pierre Auger Observatory and the Telescope Array: Joint Contribution to the 34th International Cosmic Ray Conference (ICRC 2015)

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

preprint2015arXiv

The LOFAR Multifrequency Snapshot Sky Survey (MSSS) I. Survey description and first results

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.

preprint2015arXiv

The lunar Askaryan technique: a technical roadmap

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.

preprint2015arXiv

The peculiar radio galaxy 4C 35.06: a case for recurrent AGN activity?

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.

preprint2015arXiv

The Pierre Auger Observatory: Contributions to the 34th International Cosmic Ray Conference (ICRC 2015)

Contributions of the Pierre Auger Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands

preprint2015arXiv

The radio emission pattern of air showers as measured with LOFAR - a tool for the reconstruction of the energy and the shower maximum

The pattern of the radio emission of air showers is finely sampled with the Low-Frequency ARray (LOFAR). A set of 382 measured air showers is used to test a fast, analytic parameterization of the distribution of pulse powers. Using this parameterization we are able to reconstruct the shower axis and give estimators for the energy of the air shower as well as the distance to the shower maximum.

preprint2015arXiv

Wide-Band, Low-Frequency Pulse Profiles of 100 Radio Pulsars with LOFAR

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

preprint2015arXiv

Wide-field LOFAR imaging of the field around the double-double radio galaxy B1834+620: A fresh view on a restarted AGN and doubeltjes

The existence of double-double radio galaxies (DDRGs) is evidence for recurrent jet activity in AGN, as expected from standard accretion models. A detailed study of these rare sources provides new perspectives for investigating the AGN duty cycle, AGN-galaxy feedback, and accretion mechanisms. Large catalogues of radio sources provide statistical information about the evolution of the radio-loud AGN population out to high redshifts. Using wide-field imaging with the LOFAR telescope, we study both a well-known DDRG as well as a large number of radio sources in the field of view. We present a high resolution image of the DDRG B1834+620 obtained at 144 MHz using LOFAR commissioning data. Our image covers about 100 square degrees and contains over 1000 sources. The four components of the DDRG B1834+620 have been resolved for the first time at 144 MHz. Inner lobes were found to point towards the direction of the outer lobes, unlike standard FR~II sources. Polarized emission was detected in the northern outer lobe. The high spatial resolution allows the identification of a large number of small double-lobed radio sources; roughly 10% of all sources in the field are doubles with a separation smaller than 1 arcmin. The spectral fit of the four components is consistent with a scenario in which the outer lobes are still active or the jets recently switched off, while emission of the inner lobes is the result of a mix-up of new and old jet activity. From the presence of the newly extended features in the inner lobes of the DDRG, we can infer that the mechanism responsible for their formation is the bow shock that is driven by the newly launched jet. We find that the density of the small doubles exceeds the density of FR-II sources with similar properties at 1.4 GHz, but this difference becomes smaller for low flux densities.

preprint2014arXiv

A method for high precision reconstruction of air shower Xmax using two-dimensional radio intensity profiles

The mass composition of cosmic rays contains important clues about their origin. Accurate measurements are needed to resolve long-standing issues such as the transition from Galactic to extragalactic origin, and the nature of the cutoff observed at the highest energies. Composition can be studied by measuring the atmospheric depth of the shower maximum Xmax of air showers generated by high-energy cosmic rays hitting the Earth's atmosphere. We present a new method to reconstruct Xmax based on radio measurements. The radio emission mechanism of air showers is a complex process that creates an asymmetric intensity pattern on the ground. The shape of this pattern strongly depends on the longitudinal development of the shower. We reconstruct Xmax by fitting two-dimensional intensity profiles, simulated with CoREAS, to data from the LOFAR radio telescope. In the dense LOFAR core, air showers are detected by hundreds of antennas simultaneously. The simulations fit the data very well, indicating that the radiation mechanism is now well-understood. The typical uncertainty on the reconstruction of Xmax for LOFAR showers is 17 g/cm^2.

preprint2014arXiv

A search for point sources of EeV photons

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85° to +20°, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

preprint2014arXiv

A Targeted Search for Point Sources of EeV Neutrons

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

preprint2014arXiv

Discovery of Carbon Radio Recombination Lines in absorption towards Cygnus~A

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

preprint2014arXiv

Future mmVLBI Research with ALMA: A European vision

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

preprint2014arXiv

Initial LOFAR observations of Epoch of Reionization windows: II. Diffuse polarized emission in the ELAIS-N1 field

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

preprint2014arXiv

LOFAR low-band antenna observations of the 3C295 and Bootes fields: source counts and ultra-steep spectrum sources

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

preprint2014arXiv

LOFAR observations of PSR B0943+10: profile evolution and discovery of a systematically changing profile delay in Bright mode

We present broadband, low-frequency (25-80 MHz and 110-190 MHz) LOFAR observations of PSR B0943+10, with the goal of better illuminating the nature of its enigmatic mode-switching behaviour. This pulsar shows two relatively stable states: a Bright (B) and Quiet (Q) mode, each with different characteristic brightness, profile morphology, and single-pulse properties. We model the average profile evolution both in frequency and time from the onset of each mode, and highlight the differences between the two modes. In both modes, the profile evolution can be well explained by radius-to-frequency mapping at altitudes within a few hundred kilometres of the stellar surface. If both B and Q-mode emission originate at the same magnetic latitude, then we find that the change of emission height between the modes is less than 6%. We also find that, during B-mode, the average profile is gradually shifting towards later spin phase and then resets its position at the next Q-to-B transition. The observed B-mode profile delay is frequency-independent (at least from 25-80 MHz) and asymptotically changes towards a stable value of about 0.004 in spin phase by the end of mode instance, much too large to be due to changing spin-down rate. Such a delay can be interpreted as a gradual movement of the emission cone against the pulsar's direction of rotation, with different field lines being illuminated over time. Another interesting explanation is a possible variation of accelerating potential inside the polar gap. This explanation connects the observed profile delay to the gradually evolving subpulse drift rate, which depends on the gradient of the potential across the field lines.

preprint2014arXiv

LOFAR tied-array imaging of Type III solar radio bursts

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.

preprint2014arXiv

LORA: A scintillator array for LOFAR to measure extensive air showers

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.

preprint2014arXiv

Lunar detection of ultra-high-energy cosmic rays and neutrinos with the Square Kilometre Array

The origin of the most energetic particles in nature, the ultra-high-energy (UHE) cosmic rays, is still a mystery. Only the most energetic of these have sufficiently small angular deflections to be used for directional studies, and their flux is so low that even the 3,000 km^2 Pierre Auger detector registers only about 30 cosmic rays per year of these energies. A method to provide an even larger aperture is to use the lunar Askaryan technique, in which ground-based radio telescopes search for the nanosecond radio flashes produced when a cosmic ray interacts with the Moon's surface. The technique is also sensitive to UHE neutrinos, which may be produced in the decays of topological defects from the early universe. Observations with existing radio telescopes have shown that this technique is technically feasible, and established the required procedure: the radio signal should be searched for pulses in real time, compensating for ionospheric dispersion and filtering out local radio interference, and candidate events stored for later analysis. For the Square Kilometre Array (SKA), this requires the formation of multiple tied-array beams, with high time resolution, covering the Moon, with either SKA1-LOW or SKA1-MID. With its large collecting area and broad bandwidth, the SKA will be able to detect the known flux of UHE cosmic rays using the visible lunar surface - millions of square km - as the detector, providing sufficient detections of these extremely rare particles to address the mystery of their origin.

preprint2014arXiv

Lunar occultation of the diffuse radio sky: LOFAR measurements between 35 and 80 MHz

We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<ν<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$).

preprint2014arXiv

Measuring a Cherenkov ring in the radio emission from air showers at 110-190 MHz with LOFAR

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.

preprint2014arXiv

On the radio properties of the intermediate-mass black hole candidate ESO 243-49 HLX-1

We present follow-up radio observations of ESO 243-49 HLX-1 from 2012 using the Australia Telescope Compact Array (ATCA) and the Karl G. Jansky Very Large Array (VLA). We report the detection of radio emission at the location of HLX-1 during its hard X-ray state using the ATCA. Assuming that the `Fundamental Plane' of accreting black holes is applicable, we provide an independent estimate of the black hole mass of $M_{\rm{BH}}\leq2.8^{+7.5}_{-2.1} \times 10^{6}$ M$_{\odot}$ at 90% confidence. However, we argue that the detected radio emission is likely to be Doppler-boosted and our mass estimate is an upper limit. We discuss other possible origins of the radio emission such as being due to a radio nebula, star formation, or later interaction of the flares with the large-scale environment. None of these were found adequate. The VLA observations were carried out during the X-ray outburst. However, no new radio flare was detected, possibly due to a sparse time sampling. The deepest, combined VLA data suggests a variable radio source and we briefly discuss the properties of the previously detected flares and compare them with microquasars and active galactic nuclei.

preprint2014arXiv

Polarized radio emission from extensive air showers measured with LOFAR

We present LOFAR measurements of radio emission from extensive air showers. We find that this emission is strongly polarized, with a median degree of polarization of nearly $99\%$, and that the angle between the polarization direction of the electric field and the Lorentz force acting on the particles, depends on the observer location in the shower plane. This can be understood as a superposition of the radially polarized charge-excess emission mechanism, first proposed by Askaryan and the geomagnetic emission mechanism proposed by Kahn and Lerche. We calculate the relative strengths of both contributions, as quantified by the charge-excess fraction, for $163$ individual air showers. We find that the measured charge-excess fraction is higher for air showers arriving from closer to the zenith. Furthermore, the measured charge-excess fraction also increases with increasing observer distance from the air shower symmetry axis. The measured values range from $(3.3\pm 1.0)\%$ for very inclined air showers at $25\, \mathrm{m}$ to $(20.3\pm 1.3)\%$ for almost vertical showers at $225\, \mathrm{m}$. Both dependencies are in qualitative agreement with theoretical predictions.

preprint2014arXiv

Reconstruction of inclined air showers detected with the Pierre Auger Observatory

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

preprint2014arXiv

Reconstruction of the energy and depth of maximum of cosmic-ray air-showers from LOPES radio measurements

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.

preprint2014arXiv

Searches for Large-Scale Anisotropy in the Arrival Directions of Cosmic Rays Detected above Energy of $10^{19}$ eV at the Pierre Auger Observatory and the Telescope Array

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

preprint2014arXiv

The LOFAR long baseline snapshot calibrator survey

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

preprint2014arXiv

The LOFAR Pilot Surveys for Pulsars and Fast Radio Transients

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

preprint2014arXiv

The shape of the radio wavefront of extensive air showers as measured with LOFAR

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.

preprint2014arXiv

The wavefront of the radio signal emitted by cosmic ray air showers

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.

preprint2013arXiv

A strong magnetic field around the supermassive black hole at the centre of the Galaxy

The centre of our Milky Way harbours the closest candidate for a supermassive black hole. The source is thought to be powered by radiatively inefficient accretion of gas from its environment. This form of accretion is a standard mode of energy supply for most galactic nuclei. X-ray measurements have already resolved a tenuous hot gas component from which it can be fed. However, the magnetization of the gas, a crucial parameter determining the structure of the accretion flow, remains unknown. Strong magnetic fields can influence the dynamics of the accretion, remove angular momentum from the infalling gas, expel matter through relativistic jets and lead to the observed synchrotron emission. Here we report multi-frequency measurements with several radio telescopes of a newly discovered pulsar close to the Galactic Centre and show that its unusually large Faraday rotation indicates a dynamically relevant magnetic field near the black hole. If this field is accreted down to the event horizon it provides enough magnetic flux to explain the observed emission from the black hole, from radio to X-rays.

preprint2013arXiv

Calibrating High-Precision Faraday Rotation Measurements for LOFAR and the Next Generation of Low-Frequency Radio Telescopes

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

preprint2013arXiv

Comparison of LOPES data and CoREAS simulations using a full detector simulation (ICRC2013)

The LOPES experiment at the Karlsruhe Institute of Technology, Germany, has been measuring radio emission of air showers for almost 10 years. For a better understanding of the emission process a detailed comparison of data with simulations is necessary. This is possible using a newly developed detector simulation including all LOPES detector components. After propagating a simulated event through this full detector simulation a standard LOPES like event file is written. LOPES data and CoREAS simulations can then be treated equally and the same analysis software can be applied to both. This gives the opportunity to compare data and simulations directly. Furthermore, the standard analysis software can be used with simulations which provide the possibility to check the accuracy regarding reconstruction of air shower parameters. We point out the advantages and present first results using such a full LOPES detector simulation. A comparison of LOPES data and the Monte Carlo code CoREAS based on an analysis using this detector simulation is shown.

preprint2013arXiv

Cosmic Ray Measurements with LOPES: Status and Recent Results (ARENA 2012)

LOPES is a digital antenna array at the Karlsruhe Institute of Technology, Germany, for cosmic-ray air-shower measurements. Triggered by the co-located KASCADE-Grande air-shower array, LOPES detects the radio emission of air showers via digital radio interferometry. We summarize the status of LOPES and recent results. In particular, we present an update on the reconstruction of the primary-particle properties based on almost 500 events above 100 PeV. With LOPES, the arrival direction can be reconstructed with a precision of at least 0.65°, and the energy with a precision of at least 20 %, which, however, does not include systematic uncertainties on the absolute energy scale. For many particle and astrophysics questions the reconstruction of the atmospheric depth of the shower maximum, Xmax, is important, since it yields information on the type of the primary particle and its interaction with the atmosphere. Recently, we found experimental evidence that the slope of the radio lateral distribution is indeed sensitive to the longitudinal development of the air shower, but unfortunately, the Xmax precision at LOPES is limited by the high level of anthropogenic radio background. Nevertheless, the developed methods can be transferred to next generation experiments with lower background, which should provide an Xmax precision competitive to other detection technologies.

preprint2013arXiv

Detecting cosmic rays with the LOFAR radio telescope

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

preprint2013arXiv

Differential Frequency-dependent Delay from the Pulsar Magnetosphere

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

preprint2013arXiv

Initial deep LOFAR observations of Epoch of Reionization windows: I. The North Celestial Pole

The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospective observing windows for the LOFAR EoR project will be centered at the North Celestial Pole (NCP). We present results from observations of the NCP window using the LOFAR highband antenna (HBA) array in the frequency range 115 MHz to 163 MHz. The data were obtained in April 2011 during the commissioning phase of LOFAR. We used baselines up to about 30 km. With about 3 nights, of 6 hours each, effective integration we have achieved a noise level of about 100 microJy/PSF in the NCP window. Close to the NCP, the noise level increases to about 180 microJy/PSF, mainly due to additional contamination from unsubtracted nearby sources. We estimate that in our best night, we have reached a noise level only a factor of 1.4 above the thermal limit set by the noise from our Galaxy and the receivers. Our continuum images are several times deeper than have been achieved previously using the WSRT and GMRT arrays. We derive an analytical explanation for the excess noise that we believe to be mainly due to sources at large angular separation from the NCP.

preprint2013arXiv

Investigation on the energy and mass composition of cosmic rays using LOPES radio data

The sensitivity to the mass composition as well as the reconstruction of the energy of the primary particle are explored here by leveraging the features of the radio lateral distribution function. For the purpose of this analysis, a set of events measured with the LOPES experiment is reproduced with the latest CoREAS radio simulation code. Based on simulation predictions, a method which exploits the slope of the radio lateral distribution function is developed (Slope Method) and directly applied on measurements. As a result, the possibility to reconstruct both the energy and the depth of the shower maximum of the cosmic ray air shower using radio data and achieving relatively small uncertainties is presented.

preprint2013arXiv

LOFAR detections of low-frequency radio recombination lines towards Cassiopeia A

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

preprint2013arXiv

LOFAR: The LOw-Frequency ARray

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

preprint2013arXiv

LOPES 3D - vectorial measurements of radio emission from cosmic ray induced air showers

LOPES 3D is able to measure all three components of the electric field vector of the radio emission from air showers. This allows a better comparison with emission models. The measurement of the vertical component increases the sensitivity to inclined showers. By measuring all three components of the electric field vector LOPES 3D demonstrates by how much the reconstruction accuracy of primary cosmic ray parameters increases. Thus LOPES 3D evaluates the usefulness of vectorial measurements for large scale applications.

preprint2013arXiv

LOPES 3D reconfiguration and first measurements

The Radio detection technique of high-energy cosmic rays is based on the radio signal emitted by the charged particles in an air shower due to their deflection in the Earth's magnetic field. The LOPES experiment at Karlsruhe Institute of Technology, Germany with its simple dipoles made major contributions to the revival of this technique. LOPES is working in the frequency range from 40 to 80 MHz and was reconfigured several times to improve and further develop the radio detection technique. In the current setup LOPES consists of 10 tripole antennas which measure the complete electric field vector of the radio emission from cosmic rays. LOPES is the first experiment measuring all three vectorial components at once and thereby gaining the full information about the electric field vector and not only a two-dimensional projection. Such a setup including also measurements of the vertical electric field component is expected to increase the sensitivity to inclined showers and help to advance the understanding of the emission mechanism. We present the reconfiguration and calibration procedure of LOPES 3D and discuss first measurements.

preprint2013arXiv

LOPES-3D, an antenna array for full signal detection of air-shower radio emission

To better understand the radio signal emitted by extensive air-showers and to further develop the radio detection technique of high-energy cosmic rays, the LOPES experiment was reconfigured to LOPES-3D. LOPES-3D is able to measure all three vectorial components of the electric field of radio emission from cosmic ray air showers. The additional measurement of the vertical component ought to increase the reconstruction accuracy of primary cosmic ray parameters like direction and energy, provides an improved sensitivity to inclined showers, and will help to validate simulation of the emission mechanisms in the atmosphere. LOPES-3D will evaluate the feasibility of vectorial measurements for large scale applications. In order to measure all three electric field components directly, a tailor-made antenna type (tripoles) was deployed. The change of the antenna type necessitated new pre-amplifiers and an overall recalibration. The reconfiguration and the recalibration procedure are presented and the operationality of LOPES-3D is demonstrated.

preprint2013arXiv

Mass sensitivity in the radio lateral distribution function

Measuring the mass composition of ultra-high energy cosmic rays is one of the main tasks in the cosmic rays field. Here we are exploring the composition signature in the coherent electromagnetic emission from extensive air showers, detected in the MHz frequency range. One of the experiments that successfully detects radio events in the frequency band of 40-80 MHz is the LOPES experiment at KIT. It is a digital interferometric antenna array and has the important advantage of taking data in coincidence with the particle detector array KASCADE-Grande. A possible method to look at the composition signature in the radio data, predicted by simulations, concerns the radio lateral distribution function, since its slope is strongly correlated with Xmax. Recent comparison between REAS3 simulations and LOPES data showed a significantly improved agreement in the lateral distribution function and for this reason an analysis on a possible LOPES mass signature through the slope method is promising. Trying to reproduce a realistic case, proton and iron showers are simulated with REAS3 using the LOPES selection information as input parameters. The obtained radio lateral distribution slope is analyzed in detail. The lateral slope method to look at the composition signature in the radio data is shown here and a possible signature of mass composition in the LOPES data is discussed.

preprint2013arXiv

Pierre Auger Observatory and Telescope Array: Joint Contributions to the 33rd International Cosmic Ray Conference (ICRC 2013)

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

preprint2013arXiv

Pulse Broadening Measurements from the Galactic Center Pulsar J1745--2900

We present temporal scattering measurements of single pulses and average profiles of PSR J1745--2900, a magnetar recently discovered only 3 arcsec away from Sagittarius A* (Sgr A*), from 1.2 - 18.95 GHz using the Effelsberg 100-m Radio Telescope, the Nançay Decimetric Radio Telescope, and the Jodrell Bank Lovell Telescope. Single pulse analysis shows that the integrated pulse profile above 2 GHz is dominated by pulse jitter, while below 2 GHz the pulse profile shape is dominated by scattering. The high dispersion measure and rotation measure of the magnetar suggest that it is close to Sgr A* (within ~0.1 pc). This is the first object in the GC with both pulse broadening and angular broadening measurements. We measure a pulse broadening spectral index of alpha = -3.8 +/- 0.2 and a pulse broadening time scale at 1 GHz of tau_GHz = 1.3 +/- 0.2 s, which is several orders of magnitude lower than the scattering predicted by the NE2001 model (Cordes and Lazio 2002). If this scattering timescale is representative of the GC as a whole, then previous surveys should have detected many pulsars. The lack of detections implies either our understanding of scattering in the GC is incomplete or there are fewer pulsars in the GC than previously predicted. Given that magnetars are a rare class of radio pulsar, we believe that there are many canonical and millisecond pulsars in the GC, and not surprisingly, scattering regions in the GC have complex spatial structures.

preprint2013arXiv

Reconstructing energy and Xmax of cosmic ray air showers using the radio lateral distribution measured with LOPES

The LOPES experiment, a digital radio interferometer located at KIT (Karlsruhe Institute of Technology), obtained remarkable results for the detection of radio emission from extensive air showers at MHz frequencies. Features of the radio lateral distribution function (LDF) measured by LOPES are explored in this work for a precise reconstruction of two fundamental air shower parameters: the primary energy and the shower Xmax. The method presented here has been developed on (REAS3-)simulations, and is applied to LOPES measurements. Despite the high human-made noise at the LOPES site, it is possible to reconstruct both the energy and Xmax for individual events. On the one hand, the energy resolution is promising and comparable to the one of the co-located KASCADE-Grande experiment. On the other hand, Xmax values are reconstructed with the LOPES measurements with a resolution of 90 g/cm2 . A precision on Xmax better than 30 g/cm2 is predicted and achievable in a region with a lower human-made noise level.

preprint2013arXiv

Studying Galactic interstellar turbulence through fluctuations in synchrotron emission: First LOFAR Galactic foreground detection

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

preprint2013arXiv

Thunderstorm Observations by Air-Shower Radio Antenna Arrays

Relativistic, charged particles present in extensive air showers lead to a coherent emission of radio pulses which are measured to identify the shower initiating high-energy cosmic rays. Especially during thunderstorms, there are additional strong electric fields in the atmosphere, which can lead to further multiplication and acceleration of the charged particles and thus have influence on the form and strength of the radio emission. For a reliable energy reconstruction of the primary cosmic ray by means of the measured radio signal it is very important to understand how electric fields affect the radio emission. In addition, lightning strikes are a prominent source of broadband radio emissions that are visible over very long distances. This, on the one hand, causes difficulties in the detection of the much lower signal of the air shower. On the other hand the recorded signals can be used to study features of the lightning development. The detection of cosmic rays via the radio emission and the influence of strong electric fields on this detection technique is investigated with the LOPES experiment in Karlsruhe, Germany. The important question if a lightning is initiated by the high electron density given at the maximum of a high-energy cosmic-ray air shower is also investigated, but could not be answered by LOPES. But, these investigations exhibit the capabilities of EAS radio antenna arrays for lightning studies. We report about the studies of LOPES measured radio signals of air showers taken during thunderstorms and give a short outlook to new measurements dedicated to search for correlations of lightning and cosmic rays.

preprint2013arXiv

Towards a dynamical mass of the ultraluminous X-ray source NGC 5408 X-1

We obtained multi-epoch Very Large Telescope (VLT) optical spectroscopic data in 2011 and 2012 on the ultraluminous X-ray source (ULX) NGC 5408 X-1. We confirm that the HeII\lambda4686 line has a broad component with an average FWHM of v=780\pm64 km/s with a variation of ~13% during observations spanning over 4 years, and is consistent with the origin in the accretion disc. The deepest optical spectrum does not reveal any absorption line from a donor star. Our aim was to measure the radial velocity curve and estimate the parameters of the binary system. We find an upper limit on the semi-amplitude of the radial velocity of K=132\pm42 km/s. A search for a periodic signal in the data resulted in no statistically significant period. The mass function and constraints on the binary system imply a black hole mass of less than ~510 M_sun. Whilst, a disc irradiation model may imply a black hole mass smaller than ~431-1985 M_sun, depending on inclination. Our data can also be consistent with an unexplored orbital period range from a couple of hours to a few days, thus with a stellar-mass black hole and a subgiant companion.

preprint2013arXiv

Unveiling recurrent jets of the ULX Holmberg II X-1: evidence for a massive stellar-mass black hole?

We report on the discovery of an apparent triple radio structure hidden inside the radio bubble of the ultraluminous X-ray source Holmberg II X-1. The morphology is consistent with a collimated jet structure, which is observed to emit optically thin synchrotron radiation. The central component has a steep radio spectrum and is brighter than the outer components indicating a renewed radio activity. We estimate a minimum time-averaged jet power of 2 x 10^{39} erg/s that is associated with a time-averaged isotropic X-ray luminosity of at least 4 x 10^{39} erg/s. Our results suggest that Holmberg II X-1 is powered by a black hole of M_BH \geq 25 M_sun, that is inferred to be accreting at a high Eddington rate with intermittent radio activity.

preprint2013arXiv

Vectorial Radio Interferometry with LOPES 3D

One successful detection technique for high-energy cosmic rays is based on the radio signal emitted by the charged particles in an air shower. The LOPES experiment at Karlsruhe Institute of Technology, Germany, has made major contributions to the evolution of this technique. LOPES was reconfigured several times to improve and further develop the radio detection technique. In the latest setup LOPES consisted of 10 tripole antennas. With this, LOPES 3D was the first cosmic ray experiment measuring all three vectorial field components at once and thereby gaining the full information about the electric field vector. We present an analysis based on the data taken with special focus on the benefits of a direct measurement of the vertical polarization component. We demonstrate that by measuring all polarization components the detection and reconstruction efficiency is increased and noisy single channel data can be reconstructed by utilising the information from the other two channels of one antenna station.

preprint2012arXiv

A search for ultra-high energy neutrinos in highly inclined events at the Pierre Auger Observatory

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

preprint2012arXiv

Antennas for the Detection of Radio Emission Pulses from Cosmic-Ray induced Air Showers at the Pierre Auger Observatory

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

preprint2012arXiv

Back to the Moon: The Scientific Rationale for Resuming Lunar Surface Exploration

The lunar geological record has much to tell us about the earliest history of the Solar System, the origin and evolution of the Earth-Moon system, the geological evolution of rocky planets, and the near-Earth cosmic environment throughout Solar System history. In addition, the lunar surface offers outstanding opportunities for research in astronomy, astrobiology, fundamental physics, life sciences and human physiology and medicine. This paper provides an interdisciplinary review of outstanding lunar science objectives in all of these different areas. It is concluded that addressing them satisfactorily will require an end to the 40-year hiatus of lunar surface exploration, and the placing of new scientific instruments on, and the return of additional samples from, the surface of the Moon. Some of these objectives can be achieved robotically (e.g. through targeted sample return, the deployment of geophysical networks, and the placing of antennas on the lunar surface to form radio telescopes). However, in the longer term, most of these scientific objectives would benefit significantly from renewed human operations on the lunar surface. For these reasons it is highly desirable that current plans for renewed robotic surface exploration of the Moon are developed in the context of a future human lunar exploration programme, such as that proposed by the recently formulated Global Exploration Roadmap.

preprint2012arXiv

Chandra-HETGS Observations of the Brightest Flare Seen from Sgr A*

Starting in 2012, we began an unprecedented observational program focused on the supermassive black hole in the center of our Galaxy, Sgr A*, utilizing the High Energy Transmission Gratings Spectrometer (HETGS) instrument on the Chandra X-ray Observatory. These observations will allow us to measure the quiescent X-ray spectra of Sgr A* for the first time at both high spatial and spectral resolution. The X-ray emission of Sgr A*, however, is known to flare roughly daily by factors of a few to ten times over quiescent emission levels, with rarer flares extending to factors of greater than 100 times quiescence. Here were report an observation performed on 2012 February 9 wherein we detected what is the highest peak flux and fluence flare ever observed from Sgr A*. The flare, which lasted for 5.6 ks and had a decidedly asymmetric profile with a faster decline than rise, achieved a mean absorbed 2-8 keV flux of (8.5+/-0.9)X10^{-12} erg cm^{-2} s^{-1}. The peak flux was 2.5 times higher, and the total 2-10 keV emission of the event was approximately 10^{39} erg. Only one other flare of comparable magnitude, but shorter duration, has been observed in Sgr A* by XMM-Newton in 2002 October. We perform spectral fits of this Chandra observed flare, and compare our results to the two brightest flares ever observed with XMM-Newton. We find good agreement among the fitted spectral slopes (Gamma~2) and X-ray absorbing columns (N_H~15X10^{22} cm^{-2}) for all three of these events, resolving prior differences (which are most likely due to the combined effects of pileup and spectral modeling) among Chandra and XMM-Newton observations of Sgr A* flares. We also discuss fits to the quiescent spectra of Sgr A*.

preprint2012arXiv

Description of Atmospheric Conditions at the Pierre Auger Observatory using the Global Data Assimilation System (GDAS)

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malargüe and averaged monthly models, the utility of the GDAS data is shown.

preprint2012arXiv

Experimental evidence for the sensitivity of the air-shower radio signal to the longitudinal shower development

We observe a correlation between the slope of radio lateral distributions, and the mean muon pseudorapidity of 59 individual cosmic-ray-air-shower events. The radio lateral distributions are measured with LOPES, a digital radio interferometer co-located with the multi-detector-air-shower array KASCADE-Grande, which includes a muon-tracking detector. The result proves experimentally that radio measurements are sensitive to the longitudinal development of cosmic-ray air-showers. This is one of the main prerequisites for using radio arrays for ultra-high-energy particle physics and astrophysics.

preprint2012arXiv

First LOFAR observations at very low frequencies of cluster-scale non-thermal emission: the case of Abell 2256

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

preprint2012arXiv

M87 at metre wavelengths: the LOFAR picture

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

preprint2012arXiv

The LOFAR radio environment

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

preprint2012arXiv

The Rapid Atmospheric Monitoring System of the Pierre Auger Observatory

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

preprint2012arXiv

Wide-band Simultaneous Observations of Pulsars: Disentangling Dispersion Measure and Profile Variations

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

preprint2011arXiv

Advanced functionality for radio analysis in the Offline software framework of the Pierre Auger Observatory

The advent of the Auger Engineering Radio Array (AERA) necessitates the development of a powerful framework for the analysis of radio measurements of cosmic ray air showers. As AERA performs "radio-hybrid" measurements of air shower radio emission in coincidence with the surface particle detectors and fluorescence telescopes of the Pierre Auger Observatory, the radio analysis functionality had to be incorporated in the existing hybrid analysis solutions for fluoresence and surface detector data. This goal has been achieved in a natural way by extending the existing Auger Offline software framework with radio functionality. In this article, we lay out the design, highlights and features of the radio extension implemented in the Auger Offline framework. Its functionality has achieved a high degree of sophistication and offers advanced features such as vectorial reconstruction of the electric field, advanced signal processing algorithms, a transparent and efficient handling of FFTs, a very detailed simulation of detector effects, and the read-in of multiple data formats including data from various radio simulation codes. The source code of this radio functionality can be made available to interested parties on request.

preprint2011arXiv

An air shower array for LOFAR: LORA

LOFAR is a new form of radio telescope which can detect radio emission from air showers induced by very high-energy cosmic rays. It can also look for radio emission from particle cascades on the Moon induced by ultra high-energy cosmic rays or neutrinos. To complement the radio detection, we are setting up a small particle detector array LORA (LOfar Radboud Air shower array) within an area of $\sim 300$ m diameter in the LOFAR core. It will help in triggering and confirming the radio detection of air showers with the LOFAR antennas. In this paper, we present a short overview about LORA and discuss its current status.

preprint2011arXiv

Anisotropy and chemical composition of ultra-high energy cosmic rays using arrival directions measured by the Pierre Auger Observatory

The Pierre Auger Collaboration has reported evidence for anisotropy in the distribution of arrival directions of the cosmic rays with energies $E>E_{th}=5.5\times 10^{19}$ eV. These show a correlation with the distribution of nearby extragalactic objects, including an apparent excess around the direction of Centaurus A. If the particles responsible for these excesses at $E>E_{th}$ are heavy nuclei with charge $Z$, the proton component of the sources should lead to excesses in the same regions at energies $E/Z$. We here report the lack of anisotropies in these directions at energies above $E_{th}/Z$ (for illustrative values of $Z=6,\ 13,\ 26$). If the anisotropies above $E_{th}$ are due to nuclei with charge $Z$, and under reasonable assumptions about the acceleration process, these observations imply stringent constraints on the allowed proton fraction at the lower energies.

preprint2011arXiv

LOFAR: Detecting Cosmic Rays with a Radio Telescope

LOFAR (the Low Frequency Array), a distributed digital radio telescope with stations in the Netherlands, Germany, France, Sweden, and the United Kingdom, is designed to enable full-sky monitoring of transient radio sources. These capabilities are ideal for the detection of broadband radio pulses generated in cosmic ray air showers. The core of LOFAR consists of 24 stations within 4 square kilometers, and each station contains 96 low-band antennas and 48 high-band antennas. This dense instrumentation will allow detailed studies of the lateral distribution of the radio signal in a frequency range of 10-250 MHz. Such studies are key to understanding the various radio emission mechanisms within the air shower, as well as for determining the potential of the radio technique for primary particle identification. We present the status of the LOFAR cosmic ray program, including the station design and hardware, the triggering and filtering schemes, and our initial observations of cosmic-ray-induced radio pulses.

preprint2011arXiv

Observing pulsars and fast transients with LOFAR

Low frequency radio waves, while challenging to observe, are a rich source of information about pulsars. The LOw Frequency ARray (LOFAR) is a new radio interferometer operating in the lowest 4 octaves of the ionospheric "radio window": 10-240MHz, that will greatly facilitate observing pulsars at low radio frequencies. Through the huge collecting area, long baselines, and flexible digital hardware, it is expected that LOFAR will revolutionize radio astronomy at the lowest frequencies visible from Earth. LOFAR is a next-generation radio telescope and a pathfinder to the Square Kilometre Array (SKA), in that it incorporates advanced multi-beaming techniques between thousands of individual elements. We discuss the motivation for low-frequency pulsar observations in general and the potential of LOFAR in addressing these science goals. We present LOFAR as it is designed to perform high-time-resolution observations of pulsars and other fast transients, and outline the various relevant observing modes and data reduction pipelines that are already or will soon be implemented to facilitate these observations. A number of results obtained from commissioning observations are presented to demonstrate the exciting potential of the telescope. This paper outlines the case for low frequency pulsar observations and is also intended to serve as a reference for upcoming pulsar/fast transient science papers with LOFAR.

preprint2011arXiv

Optimized Trigger for Ultra-High-Energy Cosmic-Ray and Neutrino Observations with the Low Frequency Radio Array

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

preprint2011arXiv

Search for First Harmonic Modulation in the Right Ascension Distribution of Cosmic Rays Detected at the Pierre Auger Observatory

We present the results of searches for dipolar-type anisotropies in different energy ranges above $2.5\times 10^{17}$ eV with the surface detector array of the Pierre Auger Observatory, reporting on both the phase and the amplitude measurements of the first harmonic modulation in the right-ascension distribution. Upper limits on the amplitudes are obtained, which provide the most stringent bounds at present, being below 2% at 99% $C.L.$ for EeV energies. We also compare our results to those of previous experiments as well as with some theoretical expectations.

preprint2011arXiv

Search for signatures of magnetically-induced alignment in the arrival directions measured by the Pierre Auger Observatory

We present the results of an analysis of data recorded at the Pierre Auger Observatory in which we search for groups of directionally-aligned events (or `multiplets') which exhibit a correlation between arrival direction and the inverse of the energy. These signatures are expected from sets of events coming from the same source after having been deflected by intervening coherent magnetic fields. The observation of several events from the same source would open the possibility to accurately reconstruct the position of the source and also measure the integral of the component of the magnetic field orthogonal to the trajectory of the cosmic rays. We describe the largest multiplets found and compute the probability that they appeared by chance from an isotropic distribution. We find no statistically significant evidence for the presence of multiplets arising from magnetic deflections in the present data.

preprint2011arXiv

The effect of the geomagnetic field on cosmic ray energy estimates and large scale anisotropy searches on data from the Pierre Auger Observatory

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

preprint2011arXiv

The Lateral Trigger Probability function for the Ultra-High Energy Cosmic Ray Showers detected by the Pierre Auger Observatory

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.

preprint2011arXiv

The Pierre Auger Observatory I: The Cosmic Ray Energy Spectrum and Related Measurements

Studies of the cosmic ray energy spectrum at the highest energies with the Pierre Auger Observatory

preprint2011arXiv

The Pierre Auger Observatory II: Studies of Cosmic Ray Composition and Hadronic Interaction models

Studies of the composition of the highest energy cosmic rays with the Pierre Auger Observatory, including examination of hadronic physics effects on the structure of extensive air showers.

preprint2011arXiv

The Pierre Auger Observatory III: Other Astrophysical Observations

Astrophysical observations of ultra-high-energy cosmic rays with the Pierre Auger Observatory

preprint2011arXiv

The Pierre Auger Observatory IV: Operation and Monitoring

Technical reports on operations and monitoring of the Pierre Auger Observatory

preprint2011arXiv

The Pierre Auger Observatory V: Enhancements

Ongoing and planned enhancements of the Pierre Auger Observatory

preprint2011arXiv

Trigger and Aperture of the Surface Detector Array of the Pierre Auger Observatory

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

preprint2010arXiv

A new limit on the Ultra-High-Energy Cosmic-Ray flux with the Westerbork Synthesis Radio Telescope

A particle cascade (shower) in a dielectric, for example as initiated by an ultra-high energy cosmic ray, will have an excess of electrons which will emit coherent Čerenkov radiation, known as the Askaryan effect. In this work we study the case in which such a particle shower occurs in a medium just below its surface. We show, for the first time, that the radiation transmitted through the surface is independent of the depth of the shower below the surface when observed from far away, apart from trivial absorption effects. As a direct application we use the recent results of the NuMoon project, where a limit on the neutrino flux for energies above $10^{22}$\,eV was set using the Westerbork Synthesis Radio Telescope by measuring pulsed radio emission from the Moon, to set a limit on the flux of ultra-high-energy cosmic rays.

preprint2010arXiv

Constraints on the flux of Ultra-High Energy neutrinos from WSRT observations

Ultra-high energy (UHE) neutrinos and cosmic rays initiate particle cascades underneath the Moon's surface. These cascades have a negative charge excess and radiate Cherenkov radio emission in a process known as the Askaryan effect. The optimal frequency window for observation of these pulses with radio telescopes on the Earth is around 150 MHz. By observing the Moon with the Westerbork Synthesis Radio Telescope array we are able to set a new limit on the UHE neutrino flux. The PuMa II backend is used to monitor the Moon in 4 frequency bands between 113 and 175 MHz with a sampling frequency of 40 MHz. The narrowband radio interference is digitally filtered out and the dispersive effect of the Earth's ionosphere is compensated for. A trigger system is implemented to search for short pulses. By inserting simulated pulses in the raw data, the detection efficiency for pulses of various strength is calculated. With 47.6 hours of observation time, we are able to set a limit on the UHE neutrino flux. This new limit is an order of magnitude lower than existing limits. In the near future, the digital radio array LOFAR will be used to achieve an even lower limit.

preprint2010arXiv

New method for the time calibration of an interferometric radio antenna array

Digital radio antenna arrays, like LOPES (LOFAR PrototypE Station), detect high-energy cosmic rays via the radio emission from atmospheric extensive air showers. LOPES is an array of dipole antennas placed within and triggered by the KASCADE-Grande experiment on site of the Karlsruhe Institute of Technology, Germany. The antennas are digitally combined to build a radio interferometer by forming a beam into the air shower arrival direction which allows measurements even at low signal-to-noise ratios in individual antennas. This technique requires a precise time calibration. A combination of several calibration steps is used to achieve the necessary timing accuracy of about 1 ns. The group delays of the setup are measured, the frequency dependence of these delays (dispersion) is corrected in the subsequent data analysis, and variations of the delays with time are monitored. We use a transmitting reference antenna, a beacon, which continuously emits sine waves at known frequencies. Variations of the relative delays between the antennas can be detected and corrected for at each recorded event by measuring the phases at the beacon frequencies.

preprint2010arXiv

On noise treatment in radio measurements of cosmic ray air showers

Precise measurements of the radio emission by cosmic ray air showers require an adequate treatment of noise. Unlike to usual experiments in particle physics, where noise always adds to the signal, radio noise can in principle decrease or increase the signal if it interferes by chance destructively or constructively. Consequently, noise cannot simply be subtracted from the signal, and its influence on amplitude and time measurement of radio pulses must be studied with care. First, noise has to be determined consistently with the definition of the radio signal which typically is the maximum field strength of the radio pulse. Second, the average impact of noise on radio pulse measurements at individual antennas is studied for LOPES. It is shown that a correct treatment of noise is especially important at low signal-to-noise ratios: noise can be the dominant source of uncertainty for pulse height and time measurements, and it can systematically flatten the slope of lateral distributions. The presented method can also be transfered to other experiments in radio and acoustic detection of cosmic rays and neutrinos.

preprint2010arXiv

Prospects for determining air shower characteristics through geosynchrotron emission arrival times

Using simulations of geosynchrotron radiation from extensive air showers, we present a relation between the shape of the geosynchrotron radiation front and the distance of the observer to the maximum of the air shower. By analyzing the relative arrival times of radio pulses at several radio antennas in an air shower array, this relation may be employed to estimate the depth of maximum of an extensive air shower if its impact position is known, allowing an estimate for the primary particle's species. Vice versa, the relation provides an estimate for the impact position of the shower's core if an external estimate of the depth of maximum is available. In realistic circumstances, the method delivers reconstruction uncertainties down to 30 g/cm^2 when the distance to the shower core does not exceed 7 km. The method requires that the arrival direction is known with high precision.

preprint2010arXiv

Simulation of radio emission from air showers in atmospheric electric fields

We study the effect of atmospheric electric fields on the radio pulse emitted by cosmic ray air showers. Under fair weather conditions the dominant part of the radio emission is driven by the geomagnetic field. When the shower charges are accelerated and deflected in an electric field additional radiation is emitted. We simulate this effect with the Monte Carlo code REAS2, using CORSIKA-simulated showers as input. In both codes a routine has been implemented that treats the effect of the electric field on the shower particles. We find that the radio pulse is significantly altered in background fields of the order of ~100 V/cm and higher. Practically, this means that air showers passing through thunderstorms emit radio pulses that are not a reliable measure for the shower energy. Under other weather circumstances significant electric field effects are expected to occur rarely, but nimbostratus clouds can harbor fields that are large enough. In general, the contribution of the electric field to the radio pulse has polarization properties that are different from the geomagnetic pulse. In order to filter out radio pulses that have been affected by electric field effects, radio air shower experiments should keep weather information and perform full polarization measurements of the radio signal.

preprint2010arXiv

The expansion of SN 2008iz in M82

We present first results from the ongoing radio monitoring of SN 2008iz in M82. The VLBI images reveal a shell-like structure with circular symmetry, which expands in a self-similar way. There is strong evidence of a compact component with a steep spectrum at the center of the shell. The expansion curve obtained from our VLBI observations is marginally decelerated (m = 0.89) and can be modeled simultaneously with the available radio light curves. While the results of this simultaneous fitting are not conclusive (i.e. different combinations of values of the magnetic field, CSM density profile, and electron energy distribution, provide fits to the available data with similar quality), additional observations should allow a more robust and detailed modeling.

preprint2010arXiv

Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, \nobreak{$6\times 10^{19}$eV}. The anisotropy was measured by the fraction of arrival directions that are less than $3.1^\circ$ from the position of an active galactic nucleus within 75 Mpc (using the Véron-Cetty and Véron $12^{\rm th}$ catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is $(38^{+7}_{-6})%$, compared with $21%$ expected for isotropic cosmic rays. This is down from the early estimate of $(69^{+11}_{-13})%$. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation.

preprint2010arXiv

VLBI observations of SN 2008iz: I. Expansion velocity and limits on anisotropic expansion

We present observations of the recently discovered supernova 2008iz in M82 with the VLBI High Sensitivity Array at 22 GHz, the Very Large Array at frequencies of 1.4, 4.8, 8.4, 22 and 43 GHz, and the Chandra X-ray observatory. The supernova was clearly detected on two VLBI images, separated by 11 months. The source shows a ring-like morphology and expands with a velocity of ~23000 km/s. The most likely explosion date is in mid February 2008. The measured expansion speed is a factor of ~2 higher than expected under the assumption that synchrotron self-absorption dominates the light curve at the peak, indicating that this absorption mechanism may not be important for the radio emission. We find no evidence for an asymmetric explosion. The VLA spectrum shows a broken power law, indicating that the source was still optically thick at 1.4 GHz in April 2009. Finally, we report upper limits on the X-ray emission from SN 2008iz and a second radio transient recently discovered by MERLIN observations.

preprint2009arXiv

Air Shower Measurements with LOFAR

Air showers from cosmic rays emit short, intense radio pulses. LOFAR is a new radio telescope, that is being built in the Netherlands and Europe. Designed primarily as a radio interferometer, the core of LOFAR will have a high density of radio antennas, which will be extremely well calibrated. This makes LOFAR a unique tool for the study of the radio properties of single air showers. Triggering on the radio emission from air showers means detecting a short radio pulse and discriminating real events from radio interference. At LOFAR we plan to search for pulses in the digital data stream - either from single antennas or from already beam-formed data - and calculate several parameters characterizing the pulse shape to pick out real events in a second stage. In addition, we will have a small scintillator array to test and confirm the performance of the radio only trigger.

preprint2009arXiv

Evidence for X-ray synchrotron emission from simultaneous mid-IR to X-ray observations of a strong Sgr A* flare

This paper reports measurements of Sgr A* made with NACO in L' -band (3.80 um), Ks-band (2.12 um) and H-band (1.66 um) and with VISIR in N-band (11.88 um) at the ESO VLT, as well as with XMM-Newton at X-ray (2-10 keV) wavelengths. On 4 April, 2007, a very bright flare was observed from Sgr A* simultaneously at L'-band and X-ray wavelengths. No emission was detected using VISIR. The resulting SED has a blue slope (beta > 0 for nuL_nu ~ nu^beta, consistent with nuL_nu ~ nu^0.4) between 12 micron and 3.8 micron. For the first time our high quality data allow a detailed comparison of infrared and X-ray light curves with a resolution of a few minutes. The IR and X-ray flares are simultaneous to within 3 minutes. However the IR flare lasts significantly longer than the X-ray flare (both before and after the X-ray peak) and prominent substructures in the 3.8 micron light curve are clearly not seen in the X-ray data. From the shortest timescale variations in the L'-band lightcurve we find that the flaring region must be no more than 1.2 R_S in size. The high X-ray to infrared flux ratio, blue nuL_nu slope MIR to L' -band, and the soft nuL_nu spectral index of the X-ray flare together place strong constraints on possible flare emission mechanisms. We find that it is quantitatively difficult to explain this bright X-ray flare with inverse Compton processes. A synchrotron emission scenario from an electron distribution with a cooling break is a more viable scenario.

preprint2009arXiv

Improved flux limits for neutrinos with energies above 10$^{22}$ eV from observations with the Westerbork Synthesis Radio Telescope

Particle cascades initiated by ultra-high energy (UHE) neutrinos in the lunar regolith will emit an electromagnetic pulse with a time duration of the order of nano seconds through a process known as the Askaryan effect. It has been shown that in an observing window around 150 MHz there is a maximum chance for detecting this radiation with radio telescopes commonly used in astronomy. In 50 hours of observation time with the Westerbork Synthesis Radio Telescope array we have set a new limit on the flux of neutrinos, summed over all flavors, with energies in excess of $4\times10^{22}$ eV.

preprint2009arXiv

Lateral Distribution of the Radio Signal in Extensive Air Showers Measured with LOPES

The antenna array LOPES is set up at the location of the KASCADE-Grande extensive air shower experiment in Karlsruhe, Germany and aims to measure and investigate radio pulses from Extensive Air Showers. The coincident measurements allow us to reconstruct the electric field strength at observation level in dependence of general EAS parameters. In the present work, the lateral distribution of the radio signal in air showers is studied in detail. It is found that the lateral distributions of the electric field strengths in individual EAS can be described by an exponential function. For about 20% of the events a flattening towards the shower axis is observed, preferentially for showers with large inclination angle. The estimated scale parameters R0 describing the slope of the lateral profiles range between 100 and 200 m. No evidence for a direct correlation of R0 with shower parameters like azimuth angle, geomagnetic angle, or primary energy can be found. This indicates that the lateral profile is an intrinsic property of the radio emission during the shower development which makes the radio detection technique suitable for large scale applications.

preprint2009arXiv

Monte Carlo simulations of air showers in atmospheric electric fields

The development of cosmic ray air showers can be influenced by atmospheric electric fields. Under fair weather conditions these fields are small, but the strong fields inside thunderstorms can have a significant effect on the electromagnetic component of a shower. Understanding this effect is particularly important for radio detection of air showers, since the radio emission is produced by the shower electrons and positrons. We perform Monte Carlo simulations to calculate the effects of different electric field configurations on the shower development. We find that the electric field becomes important for values of the order of 1 kV/cm. Not only can the energy distribution of electrons and positrons change significantly for such field strengths, it is also possible that runaway electron breakdown occurs at high altitudes, which is an important effect in lightning initiation.

preprint2009arXiv

Studies of Cosmic Ray Composition and Air Shower Structure with the Pierre Auger Observatory

Studies of the composition of the highest energy cosmic rays with the Pierre Auger Observatory, including examination of hadronic physics effects on the structure of extensive air showers. Submissions to the 31st ICRC, Lodz, Poland (July 2009).

preprint2009arXiv

The Fluorescence Detector of the Pierre Auger Observatory

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.

preprint2009arXiv

Universality of electron-positron distributions in extensive air showers

Using a large set of simulated extensive air showers, we investigate universality features of electron and positron distributions in very-high-energy cosmic-ray air showers. Most particle distributions depend only on the depth of the shower maximum and the number of particles in the cascade at this depth. We provide multi-dimensional parameterizations for the electron-positron distributions in terms of particle energy, vertical and horizontal momentum angle, lateral distance, and time distribution of the shower front. These parameterizations can be used to obtain realistic electron-positron distributions in extensive air showers for data analysis and simulations of Cherenkov radiation, fluorescence signal, and radio emission.

preprint2009arXiv

Upper limit on the cosmic-ray photon fraction at EeV energies from the Pierre Auger Observatory

From direct observations of the longitudinal development of ultra-high energy air showers performed with the Pierre Auger Observatory, upper limits of 3.8%, 2.4%, 3.5% and 11.7% (at 95% c.l.) are obtained on the fraction of cosmic-ray photons above 2, 3, 5 and 10 EeV (1 EeV = 10^18 eV) respectively. These are the first experimental limits on ultra-high energy photons at energies below 10 EeV. The results complement previous constraints on top-down models from array data and they reduce systematic uncertainties in the interpretation of shower data in terms of primary flux, nuclear composition and proton-air cross-section.

preprint2007arXiv

Proper Motions in the Andromeda Subgroup

This article presents results of VLBI observations of regions of H2O maser activity in the Local Group galaxies M33 and IC10. Since all position measurements were made relative to extragalactic background sources, the proper motions of the two galaxies could be measured. For M33, this provides this galaxy's three dimensional velocity, showing that this galaxy is moving with a velocity of 190 +/- 59 km\s relative to the Milky Way. For IC10, we obtain a motion of 215 +/- 42 km/s relative to the Milky Way. These measurements promise a new handle on dynamical models for the Local Group and the mass and dark matter halo of Andromeda and the Milky Way.

preprint2006arXiv

Detecting UHE Cosmics & Neutrinos off the Moon; an Optimal Radio Window

We show that at wavelengths comparable to the length of the shower produced by an Ultra-High Energy cosmic ray or neutrino, radio signals are an extremely efficient way to detect these particles. Through an example it is shown that this new approach offers, for the first time, the realistic possibility of measuring UHE neutrino fluxes below the Waxman-Bahcall limit. It is shown that in only one month of observing with the upcoming LOFAR radio telescope, cosmic-ray events can be measured beyond the GZK-limit, at a sensitivity level of two orders of magnitude below the extrapolated values.

preprint2006arXiv

Optimal Radio Window for the Detection of Ultra-High-Energy Cosmic Rays and Neutrinos off the Moon

When high-energy cosmic rays impinge on a dense dielectric medium, radio waves are produced through the Askaryan effect. We show that at wavelengths comparable to the length of the shower produced by an Ultra-High Energy cosmic ray or neutrino, radio signals are an extremely efficient way to detect these particles. Through an example it is shown that this new approach offers, for the first time, the realistic possibility of measuring UHE neutrino fluxes below the Waxman-Bahcall limit. It is shown that in only one month of observing with the upcoming LOFAR radio telescope, cosmic-ray events can be measured beyond the GZK-limit, at a sensitivity level of two orders of magnitude below the extrapolated values.

preprint2001arXiv

The jet-disk symbiosis model for Gamma Ray Bursts: cosmic ray and neutrino background contribution

The relation between the cosmological evolution of the jet-disk symbiosis model for GRBs and the cosmic rays energy distribution is presented. We used two different Star Formation Rates (SFR) as a function of redshift and a Luminosity Function (LF) distribution to obtain the distribution in fluence of GRBs in our model and compare it with the data. We show a good agreement between the fluence distribution we obtain and the corrected data for the 4B BATSE catalogue. The results we obtain are generally valid for models that use jet physics to explain GRB properties. The fluence in the gamma ray band has been used to calculate the energy in cosmic rays both in our Galaxy and at extragalactic distances as a function of the redshift. This energy input has been compared with the Galactic and extragalactic spectrum of cosmic rays and neutrinos. Using our jet disk symbiosis model, we found that in both cases GRBs cannot give any significant contribution to cosmic rays. We also estimate the neutrino background, obtaining a very low predicted flux. We also show that the fit of our model with the corrected fluence distribution of GRBs gives strong constraints of the star formation rate as a function of the redshift.

preprint1999arXiv

Gigahertz Peaked Spectrum sources from the Jodrell Bank-VLA Astrometric Survey

Observations with MERLIN at 408 MHz have been used to establish the low-frequency part of the spectra of more than a hundred compact radio sources taken from the part of the Jodrell Bank--VLA Astrometric Survey limited by 35 deg. < δ< 75 deg. These sources were selected from JVAS and other catalogues to have convex spectra between 1.4 and 8.4 GHz, characteristic of Gigahertz Peaked Spectrum (GPS) sources. We have confirmed convex shapes of the spectra of 76 objects (one half of our initial candidates) thereby yielding the largest genuine sample of GPS sources compiled so far. Seven of 17 identified quasars in the sample have large ($z\ga 2$) redshifts.

H. Falcke

What is connected

Connect this record

See the researcher in context

Building this map preview

162 published item(s)

Measuring the muon content of inclined air showers using AERA and the water-Cherenkov detectors of the Pierre Auger Observatory

A search for photons with energies above $2{\times}10^{17}$ eV using hybrid data from the low-energy extensions of the Pierre Auger Observatory

Investigating Hadronic Interactions at Ultra-High Energies with the Pierre Auger Observatory

LOPES 3D -- studies on the benefits of EAS-radio measurements with vertically aligned antennas

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

Testing effects of Lorentz invariance violation in the propagation of astroparticles with the Pierre Auger Observatory

The energy spectrum of cosmic rays beyond the turn-down around $10^{17}$ eV as measured with the surface detector of the Pierre Auger Observatory

Design, upgrade and characterization of the silicon photomultiplier front-end for the AMIGA detector at the Pierre Auger Observatory

Final results of the LOPES radio interferometer for cosmic-ray air showers

Sub-arcsecond imaging with the International LOFAR Telescope I. Foundational calibration strategy and pipeline

Sub-arcsecond imaging with the International LOFAR Telescope: II. Completion of the LOFAR Long-Baseline Calibrator Survey

Cosmic-ray anisotropies in right ascension measured by the Pierre Auger Observatory

ESA Voyage 2050 white paper -- GrailQuest: hunting for Atoms of Space and Time hidden in the wrinkle of Space-Time

Radio emission from negative lightning leader steps reveals inner meter-scale structure

Reconstructing air shower parameters with LOFAR using event specific GDAS atmospheres

Search for magnetically-induced signatures in the arrival directions of ultra-high-energy cosmic rays measured at the Pierre Auger Observatory

Studies on the response of a water-Cherenkov detector of the Pierre Auger Observatory to atmospheric muons using an RPC hodoscope

A Roadmap towards a Space-based Radio Telescope for Ultra-Low Frequency Radio Astronomy

rPICARD: A CASA-based Calibration Pipeline for VLBI Data

Search for High-energy Neutrinos from Binary Neutron Star Merger GW170817 with ANTARES, IceCube, and the Pierre Auger Observatory

A large light-mass component of cosmic rays at 10^{17} - 10^{17.5} eV from radio observations

A LOFAR census of non-recycled pulsars: average profiles, dispersion measures, flux densities, and spectra

Azimuthal asymmetry in the risetime of the surface detector signals of the Pierre Auger Observatory

Cosmic-ray energy spectrum and composition up to the ankle - the case for a second Galactic component

Evidence for a mixed mass composition at the `ankle' in the cosmic-ray spectrum

LBCS: the LOFAR Long-Baseline Calibrator Survey

LOFAR 150-MHz observations of the Boötes field: Catalogue and Source Counts

Low-radio-frequency eclipses of the redback pulsar J2215+5135 observed in the image plane with LOFAR

Measurement of the circular polarization in radio emission from extensive air showers confirms emission mechanisms

Nanosecond-level time synchronization of autonomous radio detector stations for extensive air showers

Probing spacetime around Sagittarius A* using modeled VLBI closure phases

Prototype muon detectors for the AMIGA component of the Pierre Auger Observatory

Search for correlations between the arrival directions of IceCube neutrino events and ultrahigh-energy cosmic rays detected by the Pierre Auger Observatory and the Telescope Array

Search for ultrarelativistic magnetic monopoles with the Pierre Auger Observatory

Testing Hadronic Interactions at Ultrahigh Energies with Air Showers Measured by the Pierre Auger Observatory

The MHz-peaked radio spectrum of the unusual gamma-ray source PMN J1603-4904

The Pierre Auger Observatory Upgrade - Preliminary Design Report

The radio evolution of supernova SN$\,$2008iz in M$\,$82

Timing calibration and spectral cleaning of LOFAR time series data

A LOFAR census of millisecond pulsars

Calibrating the absolute amplitude scale for air showers measured at LOFAR

Comparing LOPES measurements of air-shower radio emission with REAS 3.11 and CoREAS simulations

High-precision measurements of extensive air showers with the SKA

Imaging Jupiter's radiation belts down to 127 MHz with LOFAR

Improved absolute calibration of LOPES measurements and its impact on the comparison with REAS 3.11 and CoREAS simulations

Influence of Atmospheric Electric Fields on the Radio Emission from Extensive Air Showers

Investigation of the radio wavefront of air showers with LOPES measurements and CoREAS simulations (ARENA 2014)

Limits on Fast Radio Bursts at 145 MHz with ARTEMIS, a real-time software backend

LOFAR discovery of a quiet emission mode in PSR B0823+26

LOFAR MSSS: Detection of a low-frequency radio transient in 400 hrs of monitoring of the North Celestial Pole

LOFAR Sparse Image Reconstruction

LOFAR tied-array imaging and spectroscopy of solar S bursts

Measurement of the cosmic-ray energy spectrum above $10^{16}$ eV with the LOFAR Radboud Air Shower Array

Orbital and superorbital variability of LS I +61 303 at low radio frequencies with GMRT and LOFAR

Pierre Auger Observatory and Telescope Array: Joint Contributions to the 34th International Cosmic Ray Conference (ICRC 2015)

Precision measurements of cosmic ray air showers with the SKA

Probing Atmospheric Electric Fields in Thunderstorms through Radio Emission from Cosmic-Ray-Induced Air Showers

Pulsar polarisation below 200 MHz: Average profiles and propagation effects

Revised absolute amplitude calibration of the LOPES experiment

The evolution of a jet ejection of the ultraluminous X-ray source Holmberg II X-1

The IceCube Neutrino Observatory, the Pierre Auger Observatory and the Telescope Array: Joint Contribution to the 34th International Cosmic Ray Conference (ICRC 2015)

The LOFAR Multifrequency Snapshot Sky Survey (MSSS) I. Survey description and first results

The lunar Askaryan technique: a technical roadmap

The peculiar radio galaxy 4C 35.06: a case for recurrent AGN activity?

The Pierre Auger Observatory: Contributions to the 34th International Cosmic Ray Conference (ICRC 2015)

The radio emission pattern of air showers as measured with LOFAR - a tool for the reconstruction of the energy and the shower maximum

Wide-Band, Low-Frequency Pulse Profiles of 100 Radio Pulsars with LOFAR

Wide-field LOFAR imaging of the field around the double-double radio galaxy B1834+620: A fresh view on a restarted AGN and doubeltjes

A method for high precision reconstruction of air shower Xmax using two-dimensional radio intensity profiles

A search for point sources of EeV photons

A Targeted Search for Point Sources of EeV Neutrons

Discovery of Carbon Radio Recombination Lines in absorption towards Cygnus~A

Future mmVLBI Research with ALMA: A European vision

Initial LOFAR observations of Epoch of Reionization windows: II. Diffuse polarized emission in the ELAIS-N1 field