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

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.

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

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

The LOFAR radio telescope is able to measure the radio emission from cosmic ray induced air showers with hundreds of individual antennas. This allows for precision testing of the emission mechanisms for the radio signal as well as determination of the depth of shower maximum $X_{\max}$, the shower observable most sensitive to the mass of the primary cosmic ray, to better than 20 g/cm$^2$. With a densely instrumented circular area of roughly 320 m$^2$, LOFAR is targeting for cosmic ray astrophysics in the energy range $10^{16}$ - $10^{18}$ eV. In this contribution we give an overview of the status, recent results, and future plans of cosmic ray detection with the LOFAR radio telescope.

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

Contributions of the LOFAR Cosmic Ray Key Science Project to the 35th International Cosmic Ray Conference (ICRC 2017)

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

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

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

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.

The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. A significant increase in the number of detected UHECR is expected to be achieved by employing Earth's moon as detector, and search for short radio pulses that are emitted when a particle interacts in the lunar rock. Observation of these short pulses with current and future radio telescopes also allows to search for the even lower fluxes of neutrinos with energies above $10^{22}$ eV, that are predicted in certain Grand-Unifying-Theories (GUTs), and e.g. models for super-heavy dark matter (SHDM). In this contribution we present the initial design for such a search with the LOFAR radio telescope.

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.

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.

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.

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.

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.

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.

The propagation of cosmic rays inside our galaxy plays a fundamental role in shaping their injection spectra into those observed at Earth. One of the best tools to investigate this issue is the ratio of fluxes for secondary and primary species. The boron-to-carbon (B/C) ratio, in particular, is a sensitive probe to investigate propagation mechanisms. This paper presents new measurements of the absolute fluxes of boron and carbon nuclei, as well as the B/C ratio, from the PAMELA space experiment. The results span the range 0.44 - 129 GeV/n in kinetic energy for data taken in the period July 2006 - March 2008.

Precision measurements of the positron component in the cosmic radiation provide important information about the propagation of cosmic rays and the nature of particle sources in our Galaxy. The satellite-borne experiment PAMELA has been used to make a new measurement of the cosmic-ray positron flux and fraction that extends previously published measurements up to 300 GeV in kinetic energy. The combined measurements of the cosmic-ray positron energy spectrum and fraction provide a unique tool to constrain interpretation models. During the recent solar minimum activity period from July 2006 to December 2009 approximately 24500 positrons were observed. The results cannot be easily reconciled with purely secondary production and additional sources of either astrophysical or exotic origin may be required.

The PAMELA satellite borne experiment is designed to study cosmic rays with great accuracy in a wide energy range. One of PAMELA's main goal is the study of the antimatter component of cosmic rays. The experiment, housed on board the Russian satellite Resurs-DK1, was launched on June 15th 2006 and it is still taking data. In this work we present the measurement of galactic positron energy spectrum in the energy range between 500 MeV and few hundred GeV.

The satellite-borne experiment PAMELA has been used to make new measurements of cosmic ray H and He isotopes. The isotopic composition was measured between 100 and 600 MeV/n for hydrogen and between 100 and 900 MeV/n for helium isotopes over the 23rd solar minimum from July 2006 to December 2007. The energy spectrum of these components carries fundamental information regarding the propagation of cosmic rays in the galaxy which are competitive with those obtained from other secondary to primary measurements such as B/C.

The energy spectra of galactic cosmic rays carry fundamental information regarding their origin and propagation. These spectra, when measured near Earth, are significantly affected by the solar magnetic field. A comprehensive description of the cosmic radiation must therefore include the transport and modulation of cosmic rays inside the heliosphere. During the end of the last decade the Sun underwent a peculiarly long quiet phase well suited to study modulation processes. In this paper we present proton spectra measured from July 2006 to December 2009 by PAMELA. The large collected statistics of protons allowed the time variation to be followed on a nearly monthly basis down to 400 MV. Data are compared with a state-of-the-art three-dimensional model of solar modulation.

The existence of a significant flux of antiprotons confined to Earth's magnetosphere has been considered in several theoretical works. These antiparticles are produced in nuclear interactions of energetic cosmic rays with the terrestrial atmosphere and accumulate in the geomagnetic field at altitudes of several hundred kilometers. A contribution from the decay of albedo antineutrons has been hypothesized in analogy to proton production by neutron decay, which constitutes the main source of trapped protons at energies above some tens of MeV. This Letter reports the discovery of an antiproton radiation belt around the Earth. The trapped antiproton energy spectrum in the South Atlantic Anomaly (SAA) region has been measured by the PAMELA experiment for the kinetic energy range 60--750 MeV. A measurement of the atmospheric sub-cutoff antiproton spectrum outside the radiation belts is also reported. PAMELA data show that the magnetospheric antiproton flux in the SAA exceeds the cosmic-ray antiproton flux by three orders of magnitude at the present solar minimum, and exceeds the sub-cutoff antiproton flux outside radiation belts by four orders of magnitude, constituting the most abundant source of antiprotons near the Earth.

We present the space spectrometer PAMELA observations of proton and helium fluxes during the December 13 and 14, 2006 solar particle events. This is the first direct measurement of the solar energetic particles in space with a single instrument in the energy range from $\sim$ 80 MeV/n up to $\sim$ 3 GeV/n. In the event of December 13 measured energy spectra of solar protons and helium were compared with results obtained by neutron monitors and other detectors. Our measurements show a spectral behaviour different from those derived from the neutron monitor network. No satisfactory analytical fitting was found for the energy spectra. During the first hours of the December 13 event solar energetic particles spectra were close to the exponential form demonstrating rather significant temporal evolution. Solar He with energy up to ~1 GeV/n was recorded on December 13. In the event of December 14 energy of solar protons reached ~600 MeV whereas maximum energy of He was below 100 MeV/n. The spectra were slightly bended in the lower energy range and preserved their form during the second event. Difference in the particle flux appearance and temporal evolution in these two events may argue for a special conditions leading to acceleration of solar particles up to relativistic energies.

Protons and helium nuclei are the most abundant components of the cosmic radiation. Precise measurements of their fluxes are needed to understand the acceleration and subsequent propagation of cosmic rays in the Galaxy. We report precision measurements of the proton and helium spectra in the rigidity range 1 GV-1.2 TV performed by the satellite-borne experiment PAMELA. We find that the spectral shapes of these two species are different and cannot be well described by a single power law. These data challenge the current paradigm of cosmic-ray acceleration in supernova remnants followed by diffusive propagation in the Galaxy. More complex processes of acceleration and propagation of cosmic rays are required to explain the spectral structures observed in our data.

Precision measurements of the electron component in the cosmic radiation provide important information about the origin and propagation of cosmic rays in the Galaxy. Here we present new results regarding negatively charged electrons between 1 and 625 GeV performed by the satellite-borne experiment PAMELA. This is the first time that cosmic-ray electrons have been identified above 50 GeV. The electron spectrum can be described with a single power law energy dependence with spectral index -3.18 +- 0.05 above the energy region influenced by the solar wind (> 30 GeV). No significant spectral features are observed and the data can be interpreted in terms of conventional diffusive propagation models. However, the data are also consistent with models including new cosmic-ray sources that could explain the rise in the positron fraction.

The satellite-borne experiment PAMELA has been used to make a new measurement of the cosmic-ray antiproton flux and the antiproton-to-proton flux ratio which extends previously published measurements down to 60 MeV and up to 180 GeV in kinetic energy. During 850 days of data acquisition approximately 1500 antiprotons were observed. The measurements are consistent with purely secondary production of antiprotons in the galaxy. More precise secondary production models are required for a complete interpretation of the results.

The PAMELA satellite experiment has measured the cosmic-ray positron fraction between 1.5 GeV and 100 GeV. The need to reliably discriminate between the positron signal and proton background has required the development of an ad hoc analysis procedure. In this paper, a method for positron identification is described and its stability and capability to yield a correct background estimate is shown. The analysis includes new experimental data, the application of three different fitting techniques for the background sample and an estimate of systematic uncertainties due to possible inaccuracies in the background selection. The new experimental results confirm both solar modulation effects on cosmic-rays with low rigidities and an anomalous positron abundance above 10 GeV.

The 15th of June 2006, the PAMELA satellite-borne experiment was launched from the Baikonur cosmodrome and it has been collecting data since July 2006. The apparatus comprises a time-of-flight system, a silicon-microstrip magnetic spectrometer, a silicon-tungsten electromagnetic calorimeter, an anticoincidence system, a shower tail counter scintillator and a neutron detector. The combination of these devices allows precision studies of the charged cosmic radiation to be conducted over a wide energy range (100 MeV -- 100's GeV) with high statistics. The primary scientific goal is the measurement of the antiproton and positron energy spectrum in order to search for exotic sources, such as dark matter particle annihilations. PAMELA is also searching for primordial antinuclei (anti-helium) and testing cosmic-ray propagation models through precise measurements of the anti-particle energy spectrum and precision studies of light nuclei and their isotopes. Moreover, PAMELA is investigating phenomena connected with solar and earth physics.