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

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

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

Abstract: The Pierre Auger Observatory is making significant contributions towards understanding the nature and origin of ultra-high energy cosmic rays. One of its main challenges is the monitoring of the atmosphere, both in terms of its state variables and its optical properties. The aim of this work is to analyse aerosol optical depth tau(a)(z) values measured from 2004 to 2012 at the observatory, which is located in a remote and relatively unstudied area of Pampa Amarilla, Argentina. The aerosol optical depth is in average quite low – annual mean tau(a)(3.5 km) similar to 0.04 – and shows a seasonal trend with a winter minimum – tau(a)(3.5 km) – 0.03 -, and a summer maximum – tau(a)(3.5 km) similar to 0.06 -, and an unexpected increase from August to September tau(a)(35 km) similar to 0.055. We computed backward trajectories for the years 2005 to 2012 to interpret the air mass origin. Winter nights with low aerosol concentrations show air masses originating from the Pacific Ocean. Average concentrations are affected by continental sources (wind-blown dust and urban pollution), whilst the peak observed in September and October could be linked to biomass burning in the northern part of Argentina or air pollution coming from surrounding urban areas.

Abstract: The cosmic energy density in the form of radiation before and during Big Bang Nucleosynthesis (BBN) is typically parameterized in terms of the effective number of neutrinos N-eff. This quantity, in case of no extra degrees of freedom, depends upon the chemical potential and the temperature characterizing the three active neutrino distributions, as well as by their possible non-thermal features. In the present analysis we determine the upper bounds that BBN places on N-eff from primordial neutrino-antineutrino asymmetries, with a careful treatment of the dynamics of neutrino oscillations. We consider quite a wide range for the total lepton number in the neutrino sector, eta(nu) = eta(nu e) + eta(nu mu) + eta(nu tau) and the initial electron neutrino asymmetry eta(in)(nu e), solving the corresponding kinetic equations which rule the dynamics of neutrino (antineutrino) distributions in phase space due to collisions, pair processes and flavor oscillations. New bounds on both the total lepton number in the neutrino sector and the nu(e)-(nu) over bar (e) asymmetry at the onset of BBN are obtained fully exploiting the time evolution of neutrino distributions, as well as the most recent determinations of primordial H-2/H density ratio and He-4 mass fraction. Note that taking the baryon fraction as measured by WMAP, the H-2/H abundance plays a relevant role in constraining the allowed regions in the eta(nu)-eta(in)(nu e) plane. These bounds fix the maximum contribution of neutrinos with primordial asymmetries to N-eff as a function of the mixing parameter theta(13), and point out the upper bound N-eff less than or similar to 3.4. Comparing these results with the forthcoming measurement of N-eff by the Planck satellite will likely provide insight on the nature of the radiation content of the universe.