-
An in-situ calibration of a logarithmic periodic dipole antenna with a
frequency coverage of 30 MHz to 80 MHz is performed. Such antennas are part of
a radio station system used for detection of cosmic ray induced air showers at
the Engineering Radio Array of the Pierre Auger Observatory, the so-called
Auger Engineering Radio Array (AERA). The directional and frequency
characteristics of the broadband antenna are investigated using a remotely
piloted aircraft (RPA) carrying a small transmitting antenna. The antenna
sensitivity is described by the vector effective length relating the measured
voltage with the electric-field components perpendicular to the incoming signal
direction. The horizontal and meridional components are determined with an
overall uncertainty of 7.4^{+0.9}_{-0.3} % and 10.3^{+2.8}_{-1.7} %
respectively. The measurement is used to correct a simulated response of the
frequency and directional response of the antenna. In addition, the influence
of the ground conductivity and permittivity on the antenna response is
simulated. Both have a negligible influence given the ground conditions
measured at the detector site. The overall uncertainties of the vector
effective length components result in an uncertainty of 8.8^{+2.1}_{-1.3} % in
the square root of the energy fluence for incoming signal directions with
zenith angles smaller than 60{\deg}.
-
We present a combined fit of a simple astrophysical model of UHECR sources to
both the energy spectrum and mass composition data measured by the Pierre Auger
Observatory. The fit has been performed for energies above $5 \cdot 10^{18}$
eV, i.e.~the region of the all-particle spectrum above the so-called "ankle"
feature. The astrophysical model we adopted consists of identical sources
uniformly distributed in a comoving volume, where nuclei are accelerated
through a rigidity-dependent mechanism. The fit results suggest sources
characterized by relatively low maximum injection energies, hard spectra and
heavy chemical composition. We also show that uncertainties about physical
quantities relevant to UHECR propagation and shower development have a
non-negligible impact on the fit results.
-
A new analysis of the dataset from the Pierre Auger Observatory provides
evidence for anisotropy in the arrival directions of ultra-high-energy cosmic
rays on an intermediate angular scale, which is indicative of excess arrivals
from strong, nearby sources. The data consist of 5514 events above 20 EeV with
zenith angles up to 80 deg recorded before 2017 April 30. Sky models have been
created for two distinct populations of extragalactic gamma-ray emitters:
active galactic nuclei from the second catalog of hard Fermi-LAT sources (2FHL)
and starburst galaxies from a sample that was examined with Fermi-LAT.
Flux-limited samples, which include all types of galaxies from the Swift-BAT
and 2MASS surveys, have been investigated for comparison. The sky model of
cosmic-ray density constructed using each catalog has two free parameters, the
fraction of events correlating with astrophysical objects and an angular scale
characterizing the clustering of cosmic rays around extragalactic sources. A
maximum-likelihood ratio test is used to evaluate the best values of these
parameters and to quantify the strength of each model by contrast with
isotropy. It is found that the starburst model fits the data better than the
hypothesis of isotropy with a statistical significance of 4.0 sigma, the
highest value of the test statistic being for energies above 39 EeV. The three
alternative models are favored against isotropy with 2.7-3.2 sigma
significance. The origin of the indicated deviation from isotropy is examined
and prospects for more sensitive future studies are discussed.
-
We present a new method for probing the hadronic interaction models at
ultra-high energy and extracting details about mass composition. This is done
using the time profiles of the signals recorded with the water-Cherenkov
detectors of the Pierre Auger Observatory. The profiles arise from a mix of the
muon and electromagnetic components of air-showers. Using the risetimes of the
recorded signals we define a new parameter, which we use to compare our
observations with predictions from simulations. We find, firstly,
inconsistencies between our data and predictions over a greater energy range
and with substantially more events than in previous studies. Secondly, by
calibrating the new parameter with fluorescence measurements from observations
made at the Auger Observatory, we can infer the depth of shower maximum for a
sample of over 81,000 events extending from 0.3 EeV to over 100 EeV. Above 30
EeV, the sample is nearly fourteen times larger than currently available from
fluorescence measurements and extending the covered energy range by half a
decade. The energy dependence of the average depth of shower maximum is
compared to simulations and interpreted in terms of the mean of the logarithmic
mass. We find good agreement with previous work and extend the measurement of
the mean depth of shower maximum to greater energies than before, reducing
significantly the statistical uncertainty associated with the inferences about
mass composition.
-
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.
-
AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the
Pierre Auger Observatory designed to extend its energy range of detection and
to directly measure the muon content of the cosmic ray primary particle
showers. The array will be formed by an infill of surface water-Cherenkov
detectors associated with buried scintillation counters employed for muon
counting. Each counter is composed of three scintillation modules, with a 10
m$^2$ detection area per module. In this paper, a new generation of detectors,
replacing the current multi-pixel photomultiplier tube (PMT) with silicon photo
sensors (aka. SiPMs), is proposed. The selection of the new device and its
front-end electronics is explained. A method to calibrate the counting system
that ensures the performance of the detector is detailed. This method has the
advantage of being able to be carried out in a remote place such as the one
where the detectors are deployed. High efficiency results, i.e. 98 % efficiency
for the highest tested overvoltage, combined with a low probability of
accidental counting ($\sim$2 %), show a promising performance for this new
system.
-
We present a novel method to measure precisely the relative spectral response
of the fluorescence telescopes of the Pierre Auger Observatory. We used a
portable light source based on a xenon flasher and a monochromator to measure
the relative spectral efficiencies of eight telescopes in steps of 5 nm from
280 nm to 440 nm. Each point in a scan had approximately 2 nm FWHM out of the
monochromator. Different sets of telescopes in the observatory have different
optical components, and the eight telescopes measured represent two each of the
four combinations of components represented in the observatory. We made an
end-to-end measurement of the response from different combinations of optical
components, and the monochromator setup allowed for more precise and complete
measurements than our previous multi-wavelength calibrations. We find an
overall uncertainty in the calibration of the spectral response of most of the
telescopes of 1.5% for all wavelengths; the six oldest telescopes have larger
overall uncertainties of about 2.2%. We also report changes in physics
measureables due to the change in calibration, which are generally small.
-
Contributions of the Pierre Auger Collaboration to the 35th International
Cosmic Ray Conference (ICRC 2017), 12-20 July 2017, Bexco, Busan, Korea.
-
Cosmic rays are atomic nuclei arriving from outer space that reach the
highest energies observed in nature. Clues to their origin come from studying
the distribution of their arrival directions. Using $3 \times 10^4$ cosmic rays
above $8 \times 10^{18}$ electron volts, recorded with the Pierre Auger
Observatory from a total exposure of 76,800 square kilometers steradian year,
we report an anisotropy in the arrival directions. The anisotropy, detected at
more than the 5.2$\sigma$ level of significance, can be described by a dipole
with an amplitude of $6.5_{-0.9}^{+1.3}$% towards right ascension $\alpha_{d} =
100 \pm 10$ degrees and declination $\delta_{d} = -24_{-13}^{+12}$ degrees.
That direction indicates an extragalactic origin for these ultra-high energy
particles.
-
We report a multi-resolution search for anisotropies in the arrival
directions of cosmic rays detected at the Pierre Auger Observatory with local
zenith angles up to $80^\circ$ and energies in excess of 4 EeV ($4 \times
10^{18}$ eV). This search is conducted by measuring the angular power spectrum
and performing a needlet wavelet analysis in two independent energy ranges.
Both analyses are complementary since the angular power spectrum achieves a
better performance in identifying large-scale patterns while the needlet
wavelet analysis, considering the parameters used in this work, presents a
higher efficiency in detecting smaller-scale anisotropies, potentially
providing directional information on any observed anisotropies. No deviation
from isotropy is observed on any angular scale in the energy range between 4
and 8 EeV. Above 8 EeV, an indication for a dipole moment is captured; while no
other deviation from isotropy is observed for moments beyond the dipole one.
The corresponding $p$-values obtained after accounting for searches blindly
performed at several angular scales, are $1.3 \times 10^{-5}$ in the case of
the angular power spectrum, and $2.5 \times 10^{-3}$ in the case of the needlet
analysis. While these results are consistent with previous reports making use
of the same data set, they provide extensions of the previous works through the
thorough scans of the angular scales.
-
A search for ultra-high energy photons with energies above 1 EeV is performed
using nine years of data collected by the Pierre Auger Observatory in hybrid
operation mode. An unprecedented separation power between photon and hadron
primaries is achieved by combining measurements of the longitudinal air-shower
development with the particle content at ground measured by the fluorescence
and surface detectors, respectively. Only three photon candidates at energies 1
- 2 EeV are found, which is compatible with the expected hadron-induced
background. Upper limits on the integral flux of ultra-high energy photons of
0.027, 0.009, 0.008, 0.008 and 0.007 km$^{-2}$ sr$^{-1}$ yr$^{-1}$ are derived
at 95% C.L. for energy thresholds of 1, 2, 3, 5 and 10 EeV. These limits bound
the fractions of photons in the all-particle integral flux below 0.1%, 0.15%,
0.33%, 0.85% and 2.7%. For the first time the photon fraction at EeV energies
is constrained at the sub-percent level. The improved limits are below the flux
of diffuse photons predicted by some astrophysical scenarios for cosmogenic
photon production. The new results rule-out the early top-down models $-$ in
which ultra-high energy cosmic rays are produced by, e.g., the decay of
super-massive particles $-$ and challenge the most recent super-heavy dark
matter models.
-
Simultaneous measurements of air showers with the fluorescence and surface
detectors of the Pierre Auger Observatory allow a sensitive search for EeV
photon point sources. Several Galactic and extragalactic candidate objects are
grouped in classes to reduce the statistical penalty of many trials from that
of a blind search and are analyzed for a significant excess above the
background expectation. The presented search does not find any evidence for
photon emission at candidate sources, and combined $p$-values for every class
are reported. Particle and energy flux upper limits are given for selected
candidate sources. These limits significantly constrain predictions of EeV
proton emission models from non-transient Galactic and nearby extragalactic
sources, as illustrated for the particular case of the Galactic center region.
-
On September 14, 2015 the Advanced LIGO detectors observed their first
gravitational-wave (GW) transient GW150914. This was followed by a second GW
event observed on December 26, 2015. Both events were inferred to have arisen
from the merger of black holes in binary systems. Such a system may emit
neutrinos if there are magnetic fields and disk debris remaining from the
formation of the two black holes. With the surface detector array of the Pierre
Auger Observatory we can search for neutrinos with energy above 100 PeV from
point-like sources across the sky with equatorial declination from about -65
deg. to +60 deg., and in particular from a fraction of the 90% confidence-level
(CL) inferred positions in the sky of GW150914 and GW151226. A targeted search
for highly-inclined extensive air showers, produced either by interactions of
downward-going neutrinos of all flavors in the atmosphere or by the decays of
tau leptons originating from tau-neutrino interactions in the Earth's crust
(Earth-skimming neutrinos), yielded no candidates in the Auger data collected
within $\pm 500$ s around or 1 day after the coordinated universal time (UTC)
of GW150914 and GW151226, as well as in the same search periods relative to the
UTC time of the GW candidate event LVT151012. From the non-observation we
constrain the amount of energy radiated in ultrahigh-energy neutrinos from such
remarkable events.
-
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.
-
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.
-
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 $\gamma=10^9$ to
$2.5 \times10^{-21}$ (cm$^{2}$ sr s)$^{-1}$ for $\gamma=10^{12}$. These results
- the first obtained with a UHECR detector - improve previously published
limits by up to an order of magnitude.
-
Interplanetary imagers provide 2D projected views of the densest plasma parts
of interplanetary coronal mass ejections (ICMEs) while in situ measurements
provide magnetic field and plasma parameter measurements along the spacecraft
trajectory, so along a 1D cut. As such, the data only give a partial view of
their 3D structures. By studying a large number of ICMEs, crossed at different
distances from their apex, we develop statistical methods to obtain a
quantitative generic 3D shape of ICMEs. In a first approach we theoretically
obtain the expected statistical distribution of the shock-normal orientation
from assuming simple models of 3D shock shapes, including distorted profiles,
and compare their compatibility with observed distributions. In a second
approach we use the shock normal and the flux rope axis orientations, as well
as the impact parameter, to provide statistical information across the
spacecraft trajectory. The study of different 3D shock models shows that the
observations are compatible with a shock symmetric around the Sun-apex line as
well as with an asymmetry up to an aspect ratio around 3. Moreover, flat or
dipped shock surfaces near their apex can only be rare cases. Next, the sheath
thickness and the ICME velocity have no global trend along the ICME front.
Finally, regrouping all these new results and the ones of our previous
articles, we provide a quantitative ICME generic 3D shape, including the global
shape of the shock, the sheath and the flux rope. The obtained quantitative
generic ICME shape will have implications for several aims. For example, it
constrains the output of typical ICME numerical simulations. It is also a base
to develop deeper studies of the transport of high energy solar and cosmic
particles during an ICME propagation, as well as for modeling and forecasting
space weather conditions near Earth.
-
Interplanetary coronal mass ejections (ICMEs) are the interplanetary
manifestations of solar eruptions. The overtaken solar wind forms a sheath of
compressed plasma at the front of ICMEs. Magnetic clouds (MCs) are a subset of
ICMEs with specific properties (e.g. the presence of a flux rope). When ICMEs
pass near Earth, ground observations indicate that the flux of galactic cosmic
rays (GCRs) decreases. The main aims of this paper are to find: common plasma
and magnetic properties of different ICME sub-structures, and which ICME
properties affect the flux of GCRs near Earth. We use a superposed epoch method
applied to a large set of ICMEs observed \insitu\ by the spacecraft ACE,
between 1998 and 2006. We also apply a superposed epoch analysis on GCRs time
series observed with the McMurdo neutron monitors. We find that slow MCs at 1
AU have on average more massive sheaths. We conclude that it is because they
are more effectively slowed down by drag during their travel from the Sun. Slow
MCs also have a more symmetric magnetic field and sheaths expanding similarly
as their following MC, while in contrast, fast MCs have an asymmetric magnetic
profile and a compressing sheath in compression. In all types of MCs, we find
that the proton density and the temperature, as well as the magnetic
fluctuations can diffuse within the front of the MC due to 3D reconnection.
Finally, we derive a quantitative model which describes the decrease of cosmic
rays as a function of the amount of magnetic fluctuations and field strength.
The obtained typical profiles of sheath/MC/GCR properties corresponding to
slow, mid, and fast ICMEs, can be used for forecasting/modelling these events,
and to better understand the transport of energetic particles in ICMEs. They
are also useful for improving future operative space weather activities.
-
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.
-
The Latin American Giant Observatory (LAGO) is an extended cosmic ray
observatory composed by a network of water-Cherenkov detectors spanning over
different sites located at significantly different altitudes (from sea level up
to more than $5000$\,m a.s.l.) and latitudes across Latin America, covering a
huge range of geomagnetic rigidity cut-offs and atmospheric absorption/reaction
levels. This detection network is designed to measure the temporal evolution of
the radiation flux at ground level with extreme detail. The LAGO project is
mainly oriented to perform basic research in three branches: high energy
phenomena, space weather and atmospheric radiation at ground level. LAGO is
built and operated by the LAGO Collaboration, a non-centralized collaborative
union of more than 30 institutions from ten countries.
These are the contributions of the LAGO Collaboration to the 34th
International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The
Netherlands
-
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.
-
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 \theta)_\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 \theta)_\mathrm{max}$.
-
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.
-
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.
-
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.
An in-situ calibration of a logarithmic periodic dipole antenna with a frequency coverage of 30 MHz to 80 MHz is performed. Such antennas are part of a radio station system used for detection of cosmic ray induced air showers at the Engineering Radio Array of the Pierre Auger Observatory, the so-called Auger Engineering Radio Array (AERA). The directional and frequency characteristics of the broadband antenna are investigated using a remotely piloted aircraft (RPA) carrying a small transmitting antenna. The antenna sensitivity is described by the vector effective length relating the measured voltage with the electric-field components perpendicular to the incoming signal direction. The horizontal and meridional components are determined with an overall uncertainty of 7.4^{+0.9}_{-0.3} % and 10.3^{+2.8}_{-1.7} % respectively. The measurement is used to correct a simulated response of the frequency and directional response of the antenna. In addition, the influence of the ground conductivity and permittivity on the antenna response is simulated. Both have a negligible influence given the ground conditions measured at the detector site. The overall uncertainties of the vector effective length components result in an uncertainty of 8.8^{+2.1}_{-1.3} % in the square root of the energy fluence for incoming signal directions with zenith angles smaller than 60{\deg}.
We present a combined fit of a simple astrophysical model of UHECR sources to both the energy spectrum and mass composition data measured by the Pierre Auger Observatory. The fit has been performed for energies above $5 \cdot 10^{18}$ eV, i.e.~the region of the all-particle spectrum above the so-called "ankle" feature. The astrophysical model we adopted consists of identical sources uniformly distributed in a comoving volume, where nuclei are accelerated through a rigidity-dependent mechanism. The fit results suggest sources characterized by relatively low maximum injection energies, hard spectra and heavy chemical composition. We also show that uncertainties about physical quantities relevant to UHECR propagation and shower development have a non-negligible impact on the fit results.
A new analysis of the dataset from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above 20 EeV with zenith angles up to 80 deg recorded before 2017 April 30. Sky models have been created for two distinct populations of extragalactic gamma-ray emitters: active galactic nuclei from the second catalog of hard Fermi-LAT sources (2FHL) and starburst galaxies from a sample that was examined with Fermi-LAT. Flux-limited samples, which include all types of galaxies from the Swift-BAT and 2MASS surveys, have been investigated for comparison. The sky model of cosmic-ray density constructed using each catalog has two free parameters, the fraction of events correlating with astrophysical objects and an angular scale characterizing the clustering of cosmic rays around extragalactic sources. A maximum-likelihood ratio test is used to evaluate the best values of these parameters and to quantify the strength of each model by contrast with isotropy. It is found that the starburst model fits the data better than the hypothesis of isotropy with a statistical significance of 4.0 sigma, the highest value of the test statistic being for energies above 39 EeV. The three alternative models are favored against isotropy with 2.7-3.2 sigma significance. The origin of the indicated deviation from isotropy is examined and prospects for more sensitive future studies are discussed.
AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory designed to extend its energy range of detection and to directly measure the muon content of the cosmic ray primary particle showers. The array will be formed by an infill of surface water-Cherenkov detectors associated with buried scintillation counters employed for muon counting. Each counter is composed of three scintillation modules, with a 10 m$^2$ detection area per module. In this paper, a new generation of detectors, replacing the current multi-pixel photomultiplier tube (PMT) with silicon photo sensors (aka. SiPMs), is proposed. The selection of the new device and its front-end electronics is explained. A method to calibrate the counting system that ensures the performance of the detector is detailed. This method has the advantage of being able to be carried out in a remote place such as the one where the detectors are deployed. High efficiency results, i.e. 98 % efficiency for the highest tested overvoltage, combined with a low probability of accidental counting ($\sim$2 %), show a promising performance for this new system.
Contributions of the Pierre Auger Collaboration to the 35th International Cosmic Ray Conference (ICRC 2017), 12-20 July 2017, Bexco, Busan, Korea.
We report a multi-resolution search for anisotropies in the arrival directions of cosmic rays detected at the Pierre Auger Observatory with local zenith angles up to $80^\circ$ and energies in excess of 4 EeV ($4 \times 10^{18}$ eV). This search is conducted by measuring the angular power spectrum and performing a needlet wavelet analysis in two independent energy ranges. Both analyses are complementary since the angular power spectrum achieves a better performance in identifying large-scale patterns while the needlet wavelet analysis, considering the parameters used in this work, presents a higher efficiency in detecting smaller-scale anisotropies, potentially providing directional information on any observed anisotropies. No deviation from isotropy is observed on any angular scale in the energy range between 4 and 8 EeV. Above 8 EeV, an indication for a dipole moment is captured; while no other deviation from isotropy is observed for moments beyond the dipole one. The corresponding $p$-values obtained after accounting for searches blindly performed at several angular scales, are $1.3 \times 10^{-5}$ in the case of the angular power spectrum, and $2.5 \times 10^{-3}$ in the case of the needlet analysis. While these results are consistent with previous reports making use of the same data set, they provide extensions of the previous works through the thorough scans of the angular scales.
A search for ultra-high energy photons with energies above 1 EeV is performed using nine years of data collected by the Pierre Auger Observatory in hybrid operation mode. An unprecedented separation power between photon and hadron primaries is achieved by combining measurements of the longitudinal air-shower development with the particle content at ground measured by the fluorescence and surface detectors, respectively. Only three photon candidates at energies 1 - 2 EeV are found, which is compatible with the expected hadron-induced background. Upper limits on the integral flux of ultra-high energy photons of 0.027, 0.009, 0.008, 0.008 and 0.007 km$^{-2}$ sr$^{-1}$ yr$^{-1}$ are derived at 95% C.L. for energy thresholds of 1, 2, 3, 5 and 10 EeV. These limits bound the fractions of photons in the all-particle integral flux below 0.1%, 0.15%, 0.33%, 0.85% and 2.7%. For the first time the photon fraction at EeV energies is constrained at the sub-percent level. The improved limits are below the flux of diffuse photons predicted by some astrophysical scenarios for cosmogenic photon production. The new results rule-out the early top-down models $-$ in which ultra-high energy cosmic rays are produced by, e.g., the decay of super-massive particles $-$ and challenge the most recent super-heavy dark matter models.
Simultaneous measurements of air showers with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for EeV photon point sources. Several Galactic and extragalactic candidate objects are grouped in classes to reduce the statistical penalty of many trials from that of a blind search and are analyzed for a significant excess above the background expectation. The presented search does not find any evidence for photon emission at candidate sources, and combined $p$-values for every class are reported. Particle and energy flux upper limits are given for selected candidate sources. These limits significantly constrain predictions of EeV proton emission models from non-transient Galactic and nearby extragalactic sources, as illustrated for the particular case of the Galactic center region.
On September 14, 2015 the Advanced LIGO detectors observed their first gravitational-wave (GW) transient GW150914. This was followed by a second GW event observed on December 26, 2015. Both events were inferred to have arisen from the merger of black holes in binary systems. Such a system may emit neutrinos if there are magnetic fields and disk debris remaining from the formation of the two black holes. With the surface detector array of the Pierre Auger Observatory we can search for neutrinos with energy above 100 PeV from point-like sources across the sky with equatorial declination from about -65 deg. to +60 deg., and in particular from a fraction of the 90% confidence-level (CL) inferred positions in the sky of GW150914 and GW151226. A targeted search for highly-inclined extensive air showers, produced either by interactions of downward-going neutrinos of all flavors in the atmosphere or by the decays of tau leptons originating from tau-neutrino interactions in the Earth's crust (Earth-skimming neutrinos), yielded no candidates in the Auger data collected within $\pm 500$ s around or 1 day after the coordinated universal time (UTC) of GW150914 and GW151226, as well as in the same search periods relative to the UTC time of the GW candidate event LVT151012. From the non-observation we constrain the amount of energy radiated in ultrahigh-energy neutrinos from such remarkable events.
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.
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.
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 $\gamma=10^9$ to $2.5 \times10^{-21}$ (cm$^{2}$ sr s)$^{-1}$ for $\gamma=10^{12}$. These results - the first obtained with a UHECR detector - improve previously published limits by up to an order of magnitude.
Interplanetary imagers provide 2D projected views of the densest plasma parts of interplanetary coronal mass ejections (ICMEs) while in situ measurements provide magnetic field and plasma parameter measurements along the spacecraft trajectory, so along a 1D cut. As such, the data only give a partial view of their 3D structures. By studying a large number of ICMEs, crossed at different distances from their apex, we develop statistical methods to obtain a quantitative generic 3D shape of ICMEs. In a first approach we theoretically obtain the expected statistical distribution of the shock-normal orientation from assuming simple models of 3D shock shapes, including distorted profiles, and compare their compatibility with observed distributions. In a second approach we use the shock normal and the flux rope axis orientations, as well as the impact parameter, to provide statistical information across the spacecraft trajectory. The study of different 3D shock models shows that the observations are compatible with a shock symmetric around the Sun-apex line as well as with an asymmetry up to an aspect ratio around 3. Moreover, flat or dipped shock surfaces near their apex can only be rare cases. Next, the sheath thickness and the ICME velocity have no global trend along the ICME front. Finally, regrouping all these new results and the ones of our previous articles, we provide a quantitative ICME generic 3D shape, including the global shape of the shock, the sheath and the flux rope. The obtained quantitative generic ICME shape will have implications for several aims. For example, it constrains the output of typical ICME numerical simulations. It is also a base to develop deeper studies of the transport of high energy solar and cosmic particles during an ICME propagation, as well as for modeling and forecasting space weather conditions near Earth.
Interplanetary coronal mass ejections (ICMEs) are the interplanetary manifestations of solar eruptions. The overtaken solar wind forms a sheath of compressed plasma at the front of ICMEs. Magnetic clouds (MCs) are a subset of ICMEs with specific properties (e.g. the presence of a flux rope). When ICMEs pass near Earth, ground observations indicate that the flux of galactic cosmic rays (GCRs) decreases. The main aims of this paper are to find: common plasma and magnetic properties of different ICME sub-structures, and which ICME properties affect the flux of GCRs near Earth. We use a superposed epoch method applied to a large set of ICMEs observed \insitu\ by the spacecraft ACE, between 1998 and 2006. We also apply a superposed epoch analysis on GCRs time series observed with the McMurdo neutron monitors. We find that slow MCs at 1 AU have on average more massive sheaths. We conclude that it is because they are more effectively slowed down by drag during their travel from the Sun. Slow MCs also have a more symmetric magnetic field and sheaths expanding similarly as their following MC, while in contrast, fast MCs have an asymmetric magnetic profile and a compressing sheath in compression. In all types of MCs, we find that the proton density and the temperature, as well as the magnetic fluctuations can diffuse within the front of the MC due to 3D reconnection. Finally, we derive a quantitative model which describes the decrease of cosmic rays as a function of the amount of magnetic fluctuations and field strength. The obtained typical profiles of sheath/MC/GCR properties corresponding to slow, mid, and fast ICMEs, can be used for forecasting/modelling these events, and to better understand the transport of energetic particles in ICMEs. They are also useful for improving future operative space weather activities.
The Latin American Giant Observatory (LAGO) is an extended cosmic ray observatory composed by a network of water-Cherenkov detectors spanning over different sites located at significantly different altitudes (from sea level up to more than $5000$\,m a.s.l.) and latitudes across Latin America, covering a huge range of geomagnetic rigidity cut-offs and atmospheric absorption/reaction levels. This detection network is designed to measure the temporal evolution of the radiation flux at ground level with extreme detail. The LAGO project is mainly oriented to perform basic research in three branches: high energy phenomena, space weather and atmospheric radiation at ground level. LAGO is built and operated by the LAGO Collaboration, a non-centralized collaborative union of more than 30 institutions from ten countries. These are the contributions of the LAGO Collaboration to the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands
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.
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 \theta)_\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 \theta)_\mathrm{max}$.
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.
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.
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.