-
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.
-
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.
-
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.
-
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 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.
-
Planck's all sky surveys at 30-857 GHz provide an unprecedented opportunity
to follow the radio spectra of a large sample of extragalactic sources to
frequencies 2-20 times higher than allowed by past, large area, ground-based
surveys. We combine the results of the Planck Early Release Compact Source
Catalog (ERCSC) with quasi-simultaneous ground-based observations, as well as
archival data, at frequencies below or overlapping Planck frequency bands, to
validate the astrometry and photometry of the ERCSC radio sources and study the
spectral features shown in this new frequency window opened by Planck. The
ERCSC source positions and flux density scales are found to be consistent with
the ground-based observations. We present and discuss the spectral energy
distributions (SEDs) of a sample of "extreme" radio sources to illustrate the
richness of the ERCSC for the study of extragalactic radio sources. Variability
is found to play a role in the unusual spectral features of some of these
sources.
-
This paper provides an overview of the Low Frequency Instrument (LFI)
programme within the ESA Planck mission. The LFI instrument has been developed
to produce high precision maps of the microwave sky at frequencies in the range
27-77 GHz, below the peak of the cosmic microwave background (CMB) radiation
spectrum. The scientific goals are described, ranging from fundamental
cosmology to Galactic and extragalactic astrophysics. The instrument design and
development are outlined, together with the model philosophy and testing
strategy. The instrument is presented in the context of the Planck mission. The
LFI approach to ground and inflight calibration is described. We also describe
the LFI ground segment. We present the results of a number of tests
demonstrating the capability of the LFI data processing centre (DPC) to
properly reduce and analyse LFI flight data, from telemetry information to
calibrated and cleaned time ordered data, sky maps at each frequency (in
temperature and polarization), component emission maps (CMB and diffuse
foregrounds), catalogs for various classes of sources (the Early Release
Compact Source Catalogue and the Final Compact Source Catalogue). The
organization of the LFI consortium is briefly presented as well as the role of
the core team in data analysis and scientific exploitation. All tests carried
out on the LFI flight model demonstrate the excellent performance of the
instrument and its various subunits. The data analysis pipeline has been tested
and its main steps verified. In the first three months after launch, the
commissioning, calibration, performance, and verification phases will be
completed, after which Planck will begin its operational life, in which LFI
will have an integral part.
-
We calculate the spectral energy distribution (SED) of electromagnetic
radiation and the spectrum of high energy neutrinos from BL Lac objects in the
context of the Synchrotron Proton Blazar Model. In this model, the high energy
hump of the SED is due to accelerated protons, while most of the low energy
hump is due to synchrotron radiation by co-accelerated electrons. To accelerate
protons to sufficiently high energies to produce the high energy hump, rather
high magnetic fields are required. Assuming reasonable emission region volumes
and Doppler factors, we then find that in low-frequency peaked BL Lacs (LBLs),
which have higher luminosities than high-frequency peaked BL Lacs (HBLs), there
is a significant contribution to the high frequency hump of the SED from pion
photoproduction and subsequent cascading, including synchrotron radiation by
muons. In contrast, in HBLs we find that the high frequency hump of the SED is
dominated by proton synchrotron radiation. We are able to model the SED of
typical LBLs and HBLs, and to model the famous 1997 flare of Markarian 501. We
also calculate the expected neutrino output of typical BL Lac objects, and
estimate the diffuse neutrino intensity due to all BL Lacs. Because pion
photoproduction is inefficient in HBLs, as protons lose energy predominantly by
synchrotron radiation, the contribution of LBLs dominates the diffuse neutrino
intensity. We suggest that nearby LBLs may well be observable with future
high-sensitivity TeV gamma-ray telescopes.
-
We obtain the maximum diffuse neutrino intensity predicted by hadronic
photoproduction models of the type which have been applied to the jets of
active galactic nuclei or gamma ray bursts. For this, we compare the proton and
gamma ray fluxes associated with hadronic photoproduction in extragalactic
neutrino sources with the present experimental upper limit on cosmic ray
protons and the extragalactic gamma ray background, employing a transport
calculation of energetic protons traversing cosmic photon backgrounds. We take
into account the effects of the photon spectral shape in the sources on the
photoproduction process, cosmological source evolution, the optical depth for
cosmic ray ejection, and discuss the possible effects of magnetic fields in the
vicinity of the sources. For photohadronic neutrino sources which are optically
thin to the emission of neutrons we find that the cosmic ray flux imposes a
stronger bound than the extragalactic gamma ray background in the energy range
between 10^5 GeV and 10^11 GeV, as previously noted by Waxman & Bahcall (1999).
We also determine the maximum contribution from the jets of active galactic
nuclei, using constraints set to their neutron opacity by gamma-ray
observations. This present upper limit is consistent with the jets of active
galactic nuclei producing the extragalactic gamma ray background hadronically,
but we point out future observations in the GeV-to-TeV regime could lower this
limit. We also briefly discuss the contribution of gamma ray bursts to
ultra-high energy cosmic rays as it can be inferred from possible observations
or limits on their correlated neutrino fluxes.