A Bayesian Approach to Comparing Cosmic Ray Energy Spectra

A common problem in ultra-high energy cosmic ray physics is the comparison of energy spectra. The question is whether the spectra from two experiments or two regions of the sky agree within their statistical and systematic uncertainties. We develop a method to directly compare energy spectra for ultra-high energy cosmic rays from two different regions of the sky in the same experiment without reliance on agreement with a theoretical model of the energy spectra. The consistency between the two spectra is expressed in terms of a Bayes factor, defined here as the ratio of the likelihood of the two-parent source hypothesis to the likelihood of the one-parent source hypothesis. Unlike other methods, for example chi^2 tests, the Bayes factor allows for the calculation of the posterior odds ratio and correctly accounts for non-Gaussian uncertainties. The latter is particularly important at the highest energies, where the number of events is very small.Comment: 22 pages, 10 figures, accepted for publication in Ap

Observed Limits on Charge Exchange Contributions to the Diffuse X-ray Background

We present a high resolution spectrum of the diffuse X-ray background from 0.1 to 1 keV for a ~1 region of the sky centered at l=90, b=+60 using a 36-pixel array of microcalorimeters flown on a sounding rocket. With an energy resolution of 11 eV FWHM below 1 keV, the spectrum's observed line ratios help separate charge exchange contributions originating within the heliosphere from thermal emission of hot gas in the interstellar medium. The X-ray sensitivity below 1 keV was reduced by about a factor of four from contamination that occurred early in the flight, limiting the significance of the results. The observed centroid of helium-like O VII is 568+2-3 eV at 90% confidence. Since the centroid expected for thermal emission is 568.4 eV while for charge exchange is 564.2 eV, thermal emission appears to dominate for this line complex, consistent with much of the high-latitude O VII emission originating in 2-3 x 10^6 K gas in the Galactic halo. On the other hand, the observed ratio of C VI Ly gamma to Ly alpha is 0.3+-0.2. The expected ratios are 0.04 for thermal emission and 0.24 for charge exchange, indicating that charge exchange must contribute strongly to this line and therefore potentially to the rest of the ROSAT R12 band usually associated with 10^6 K emission from the Local Hot Bubble. The limited statistics of this experiment and systematic uncertainties due to the contamination require only >32% thermal emission for O VII and >20% from charge exchange for C VI at the 90% confidence level. An experimental gold coating on the silicon substrate of the array greatly reduced extraneous signals induced on nearby pixels from cosmic rays passing through the substrate, reducing the triggered event rate by a factor of 15 from a previous flight of the instrument.Comment: 14 pages, 7 figures, to be published in Ap

Performance of two Askaryan Radio Array stations and first results in the search for ultra-high energy neutrinos

Ultra-high energy neutrinos are interesting messenger particles since, if detected, they can transmit exclusive information about ultra-high energy processes in the Universe. These particles, with energies above $10^{16}\mathrm{eV}$, interact very rarely. Therefore, detectors that instrument several gigatons of matter are needed to discover them. The ARA detector is currently being constructed at South Pole. It is designed to use the Askaryan effect, the emission of radio waves from neutrino-induced cascades in the South Pole ice, to detect neutrino interactions at very high energies. With antennas distributed among 37 widely-separated stations in the ice, such interactions can be observed in a volume of several hundred cubic kilometers. Currently 3 deep ARA stations are deployed in the ice of which two have been taking data since the beginning of the year 2013. In this publication, the ARA detector "as-built" and calibrations are described. Furthermore, the data reduction methods used to distinguish the rare radio signals from overwhelming backgrounds of thermal and anthropogenic origin are presented. Using data from only two stations over a short exposure time of 10 months, a neutrino flux limit of $3 \cdot 10^{-6} \mathrm{GeV} / (\mathrm{cm^2 \ s \ sr})$ is calculated for a particle energy of 10^{18}eV, which offers promise for the full ARA detector.Comment: 21 pages, 34 figures, 1 table, includes supplementary materia

First Constraints on the Ultra-High Energy Neutrino Flux from a Prototype Station of the Askaryan Radio Array

The Askaryan Radio Array (ARA) is an ultra-high energy ($>10^{17}$ eV) cosmic neutrino detector in phased construction near the South Pole. ARA searches for radio Cherenkov emission from particle cascades induced by neutrino interactions in the ice using radio frequency antennas ($\sim150-800$ MHz) deployed at a design depth of 200 m in the Antarctic ice. A prototype ARA Testbed station was deployed at $\sim30$ m depth in the 2010-2011 season and the first three full ARA stations were deployed in the 2011-2012 and 2012-2013 seasons. We present the first neutrino search with ARA using data taken in 2011 and 2012 with the ARA Testbed and the resulting constraints on the neutrino flux from $10^{17}-10^{21}$ eV.Comment: 26 pages, 15 figures. Since first revision, added section on systematic uncertainties, updated limits and uncertainty band with improvements to simulation, added appendix describing ray tracing algorithm. Final revision includes a section on cosmic ray backgrounds. Published in Astropart. Phys.

Search for Prompt Neutrino Emission from Gamma-Ray Bursts with IceCube

We present constraints derived from a search of four years of IceCube data for a prompt neutrino flux from gamma-ray bursts (GRBs). A single low-significance neutrino, compatible with the atmospheric neutrino background, was found in coincidence with one of the 506 observed bursts. Although GRBs have been proposed as candidate sources for ultra-high energy cosmic rays, our limits on the neutrino flux disfavor much of the parameter space for the latest models. We also find that no more than $\sim1\%$ of the recently observed astrophysical neutrino flux consists of prompt emission from GRBs that are potentially observable by existing satellites.Comment: 15 pages, 3 figure

Lowering IceCube’s energy threshold for point source searches in the southern sky

Observation of a point source of astrophysical neutrinos would be a "smoking gun" signature of a cosmic-ray accelerator. While IceCube has recently discovered a diffuse flux of astrophysical neutrinos, no localized point source has been observed. Previous IceCube searches for point sources in the southern sky were restricted by either an energy threshold above a few hundred TeV or poor neutrino angular resolution. Here we present a search for southern sky point sources with greatly improved sensitivities to neutrinos with energies below 100 TeV. By selecting charged-current nu(mu) interacting inside the detector, we reduce the atmospheric background while retaining efficiency for astrophysical neutrino-induced events reconstructed with sub-degree angular resolution. The new event sample covers three years of detector data and leads to a factor of 10 improvement in sensitivity to point sources emitting below 100 TeV in the southern sky. No statistically significant evidence of point sources was found, and upper limits are set on neutrino emission from individual sources. A posteriori analysis of the highest-energy (similar to 100 TeV) starting event in the sample found that this event alone represents a 2.8 sigma deviation from the hypothesis that the data consists only of atmospheric background

The contribution of Fermi-2LAC blazars to the diffuse TeV-PeV neutrino flux

The recent discovery of a diffuse cosmic neutrino flux extending up to PeV energies raises the question of which astrophysical sources generate this signal. One class of extragalactic sources which may produce such high-energy neutrinos are blazars. We present a likelihood analysis searching for cumulative neutrino emission from blazars in the 2nd Fermi-LAT AGN catalogue (2LAC) using an IceCube neutrino dataset 2009-12 which was optimised for the detection of individual sources. In contrast to previous searches with IceCube, the populations investigated contain up to hundreds of sources, the largest one being the entire blazar sample in the 2LAC catalogue. No significant excess is observed and upper limits for the cumulative flux from these populations are obtained. These constrain the maximum contribution of the 2LAC blazars to the observed astrophysical neutrino flux to be $27 \%$ or less between around 10 TeV and 2 PeV, assuming equipartition of flavours at Earth and a single power-law spectrum with a spectral index of $-2.5$. We can still exclude that the 2LAC blazars (and sub-populations) emit more than $50 \%$ of the observed neutrinos up to a spectral index as hard as $-2.2$ in the same energy range. Our result takes into account that the neutrino source count distribution is unknown, and it does not assume strict proportionality of the neutrino flux to the measured 2LAC $\gamma$-ray signal for each source. Additionally, we constrain recent models for neutrino emission by blazars.Comment: 18 pages, 22 figure