NuRadioMC is a Monte Carlo framework designed to simulate ultra-high energy neutrino detectors that rely on the radio detection method. This method exploits the radio emission generated in the ...electromagnetic component of a particle shower following a neutrino interaction. NuRadioMC simulates everything from the neutrino interaction in a medium, the subsequent Askaryan radio emission, the propagation of the radio signal to the detector and finally the detector response. NuRadioMC is designed as a modern, modular Python-based framework, combining flexibility in detector design with user-friendliness. It includes a state-of-the-art event generator, an improved modelling of the radio emission, a revisited approach to signal propagation and increased flexibility and precision in the detector simulation. This paper focuses on the implemented physics processes and their implications for detector design. A variety of models and parameterizations for the radio emission of neutrino-induced showers are compared and reviewed. Comprehensive examples are used to discuss the capabilities of the code and different aspects of instrumental design decisions.
ABSTRACT
The repeating fast radio burst (FRB) source FRB 20200120E is exceptional because of its proximity and association with a globular cluster. Here we report 60 bursts detected with the ...Effelsberg telescope at 1.4 GHz. We observe large variations in the burst rate, and report the first FRB 20200120E ‘burst storm’, where the source suddenly became active and 53 bursts (fluence ≥0.04 Jy ms) occurred within only 40 min. We find no strict periodicity in the burst arrival times, nor any evidence for periodicity in the source’s activity between observations. The burst storm shows a steep energy distribution (power-law index α = 2.39 ± 0.12) and a bimodal wait-time distribution, with log-normal means of 0.94$^{+0.07}_{-0.06}$ s and 23.61$^{+3.06}_{-2.71}$ s. We attribute these wait-time distribution peaks to a characteristic event time-scale and pseudo-Poisson burst rate, respectively. The secondary wait-time peak at ∼1 s is ∼50 × longer than the ∼24 ms time-scale seen for both FRB 20121102A and FRB 20201124A – potentially indicating a larger emission region, or slower burst propagation. FRB 20200120E shows order-of-magnitude lower burst durations and luminosities compared with FRB 20121102A and FRB 20201124A. Lastly, in contrast to FRB 20121102A, which has observed dispersion measure (DM) variations of ΔDM > 1 pc cm−3 on month-to-year time-scales, we determine that FRB 20200120E’s DM has remained stable (ΔDM < 0.15 pc cm−3) over >10 months. Overall, the observational characteristics of FRB 20200120E deviate quantitatively from other active repeaters, but it is unclear whether it is qualitatively a different type of source.
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.
Ultra-high-energy (UHE) cosmic neutrinos interacting with the Moon’s regolith generate particle showers that emit Askaryan radiation. This radiation can be observed from the Earth using ground-based ...radio telescopes like LOFAR. We simulate the effective detection aperture for UHE neutrinos hitting the Moon. Under the same assumptions, results from this work are in good agreement with previous analytic parameterizations and Monte Carlo codes. The dependence of the effective detection aperture on the observing parameters, such as observing frequency and minimum detection threshold, and lunar characteristics like surface topography have been studied. Using a Monte Carlo simulation, we find that the detectable neutrino energy threshold is lowered when we include a realistic treatment of the inelasticity, transmission coefficient, and surface roughness. Lunar surface roughness at large scales enhances the total aperture for higher observation frequencies (
ν
≥
1
GHz
) but has no significant effect on the LOFAR aperture. However, roughness at scales small compared to the wavelength reduces the aperture at all frequencies.
We present an updated cosmic-ray mass composition analysis in the energy range 1016.8 to 1018.3 eV from 334 air showers measured with the LOFAR radio telescope and selected for minimal bias. In ...this energy range, the origin of cosmic rays is expected to shift from galactic to extragalactic sources. The analysis is based on an improved method to infer the depth of the maximum Xmax of extensive air showers from radio measurements and air shower simulations. We show results of the average and standard deviation of Xmax versus primary energy and analyze the Xmax dataset at the distribution level to estimate the cosmic ray mass composition. Our approach uses an unbinned maximum likelihood analysis, making use of existing parametrizations of the Xmax distributions per element. The analysis has been repeated for three main models of hadronic interactions. Results are consistent with a significant light-mass fraction, at best fit 23% to 39% protons plus helium, depending on the choice of hadronic interaction model. The fraction of the intermediate-mass nuclei dominates. This confirms earlier results from LOFAR, with systematic uncertainties on Xmax now lowered to 7 to 9 g / cm2. We find agreement in mass composition compared to results from Pierre Auger Observatory, within statistical and systematic uncertainties. However, in line with earlier LOFAR results, we find a slightly lower average Xmax. The values are in tension with those found at Pierre Auger Observatory but agree with results from other cosmic ray observatories based in the northern hemisphere.
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.
Abstract
The common phenomenon of lightning still harbors many secrets such as what are the conditions for lightning initiation and what is driving the discharge to propagate over several tens of ...kilometers through the atmosphere forming conducting ionized channels called leaders. Since lightning is an electric discharge phenomenon, there are positively and negatively charged leaders. In this work we report on measurements made with the LOFAR radio telescope, an instrument primarily build for radio-astronomy observations. It is observed that a negative leader rather suddenly changes, for a few milliseconds, into a mode where it radiates 100 times more VHF power than typical negative leaders after which it spawns a large number of more typical negative leaders. This mode occurs during the initial stage, soon after initiation, of all lightning flashes we have mapped (about 25). For some flashes this mode occurs also well after initiation and we show one case where it is triggered twice, some 100 ms apart. We postulate that this is indicative of a small (order of 5 km
$$^2$$
2
) high charge pocket. Lightning thus appears to be initiated exclusively in the vicinity of such a small but dense charge pocket.
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.