• We describe and report the status of a neutrino-triggered program in IceCube that generates real-time alerts for gamma-ray follow-up observations by atmospheric-Cherenkov telescopes (MAGIC and VERITAS). While IceCube is capable of monitoring the whole sky continuously, high-energy gamma-ray telescopes have restricted fields of view and in general are unlikely to be observing a potential neutrino-flaring source at the time such neutrinos are recorded. The use of neutrino-triggered alerts thus aims at increasing the availability of simultaneous multi-messenger data during potential neutrino flaring activity, which can increase the discovery potential and constrain the phenomenological interpretation of the high-energy emission of selected source classes (e.g. blazars). The requirements of a fast and stable online analysis of potential neutrino signals and its operation are presented, along with first results of the program operating between 14 March 2012 and 31 December 2015.
• We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube high-energy starting events' sample and the other with 16 high-energy track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.
• The US neutrino community gathered at the Workshop on the Intermediate Neutrino Program (WINP) at Brookhaven National Laboratory February 4-6, 2015 to explore opportunities in neutrino physics over the next five to ten years. Scientists from particle, astroparticle and nuclear physics participated in the workshop. The workshop examined promising opportunities for neutrino physics in the intermediate term, including possible new small to mid-scale experiments, US contributions to large experiments, upgrades to existing experiments, R&D plans and theory. The workshop was organized into two sets of parallel working group sessions, divided by physics topics and technology. Physics working groups covered topics on Sterile Neutrinos, Neutrino Mixing, Neutrino Interactions, Neutrino Properties and Astrophysical Neutrinos. Technology sessions were organized into Theory, Short-Baseline Accelerator Neutrinos, Reactor Neutrinos, Detector R&D and Source, Cyclotron and Meson Decay at Rest sessions.This report summarizes discussion and conclusions from the workshop.
• ### Uncertainties in Atmospheric Neutrino Fluxes(astro-ph/0611266)

Nov. 8, 2006 astro-ph, hep-ph
An evaluation of the principal uncertainties in the computation of neutrino fluxes produced in cosmic ray showers in the atmosphere is presented. The neutrino flux predictions are needed for comparison with experiment to perform neutrino oscillation studies. The paper concentrates on the main limitations which are due to hadron production uncertainties. It also treats primary cosmic ray flux uncertainties, which are at a lower level. The absolute neutrino fluxes are found to have errors of around 15% in the neutrino energy region important for contained events underground. Large cancellations of these errors occur when ratios of fluxes are considered, in particular, the $\nu_\mu/\bar{\nu}_\mu$ ratio below $E_\nu=1$ GeV, the $(\nu_\mu+\bar{\nu}_\mu)/(\nu_e+\bar{\nu}_e)$ ratio below $E_\nu=10$ GeV and the up/down ratios above $E_\nu=1$ GeV are at the 1% level. A detailed breakdown of the origin of these errors and cancellations is presented.
• ### Influence of shower fluctuations and primary composition on studies of the shower longitudinal development(astro-ph/0402092)

Feb. 4, 2004 astro-ph
We study the influence of shower fluctuations, and the possible presence of different nuclear species in the primary cosmic ray spectrum, on the experimental determination of both shower energy and the proton air inelastic cross section from studies of the longitudinal development of atmospheric showers in fluorescence experiments. We investigate the potential of track length integral and shower size at maximum as estimators of shower energy. We find that at very high energy (~10^19-10^20 eV) the error of the total energy assignment is dominated by the dependence on the hadronic interaction model, and is of the order of 5%. At lower energy (~10^17-10^18 eV), the uncertainty of the energy determination due to the limited knowledge of the primary cosmic ray composition is more important. The distribution of depth of shower maximum is discussed as a measure of the proton-air cross section. Uncertainties in a possible experimental measurement of this cross section introduced by intrinsic shower fluctuations, the model of hadronic interactions, and the unknown mixture of primary nuclei in the cosmic radiation are numerically evaluated.
• ### Atmospheric shower fluctuations and the constant intensity cut method(astro-ph/0209117)

Sept. 6, 2002 astro-ph
We explore the constant intensity cut method that is widely used for the derivation of the cosmic ray energy spectrum, for comparisons of data obtained at different atmospheric depths, for measuring average shower profiles, and for estimates of the proton-air cross section from extensive air shower data. The constant intensity cut method is based on the selection of air showers by charged particle or muon size and therefore is subject to intrinsic shower fluctuations. We demostrate that, depending on the selection method, shower fluctuations can strongly influence the characteristics of the selected showers. Furthermore, a mixture of different primaries in the cosmic ray flux complicates the interpretation of measurements based on the method of constant intensity cuts. As an example we consider data published by the Akeno Collaboration. The interpretation of the Akeno measurements suggests that more than 60-70% of cosmic ray primaries in the energy range 10^{16}-19^{17} eV are heavy nuclei. Our conclusions depend only weakly on the hadronic interaction model chosen to perform the simulations, namely SIBYLL and QGSjet.
• ### Hybrid simulations of extensive air showers(astro-ph/0205302)

May 17, 2002 astro-ph
We present a fast one dimensional hybrid method to efficiently simulate extensive air showers up to the highest observed energies. Based on precalculated pion showers and a bootstrap technique, our method predicts the average shower profile, the number of muons at detector level above several energy thresholds as well as the fluctuations of the electromagnetic and hadronic components of the shower. We study the main characteristics of proton-induced air showers up to ultra-high energy, comparing the predictions of three different hadronic interaction models: SIBYLL 1.7, SIBYLL 2.1 and QGSjet98. The influence of the hadronic interaction models on the shower evolution, in particular the elongation rate, is discussed and the applicability of analytical approximations is investigated.
• ### Flux of Atmospheric Neutrinos(hep-ph/0203272)

March 30, 2002 hep-ph
Atmospheric neutrinos produced by cosmic-ray interactions in the atmosphere are of interest for several reasons. As a beam for studies of neutrino oscillations they cover a range of parameter space hitherto unexplored by accelerator neutrino beams. The atmospheric neutrinos also constitute an important background and calibration beam for neutrino astronomy and for the search for proton decay and other rare processes. Here we review the literature on calculations of atmospheric neutrinos over the full range of energy, but with particular attention to the aspects important for neutrino oscillations. Our goal is to assess how well the properties of atmospheric neutrinos are known at present.
• ### Atmospheric Neutrino Flux: A Review of Calculations(hep-ph/0001027)

Jan. 4, 2000 hep-ph
Interpretation and understanding of the evidence for neutrino oscillations depends on knowledge of the atmospheric neutrino beam. In this talk I assess how well various features are known. The goal is to determine to what extent uncertainties in the neutrino beam may limit the conclusions about neutrino properties and which features of the evidence for neutrino oscillations are most robust.
• ### Pion production in proton collisions with light nuclei: implications for atmospheric neutrinos(hep-ph/9911394)

Nov. 17, 1999 hep-ph
Differences among calculations of the atmospheric neutrino beam can be traced in large part to differences in the representation of pion production by protons interacting with nuclei in the atmosphere. In this paper we review the existing data with the goal of determining the regions of phase space in which new measurements could help to improve the input to the calculations.
• ### Geomagnetic Effects on Atmospheric Neutrinos(astro-ph/9803093)

March 9, 1998 astro-ph
Geomagnetic effects distort the zenith angle distribution of sub--GeV and few--GeV atmospheric neutrinos, breaking the up--down symmetry that would be present in the absence of neutrino oscillations and without a geomagnetic field. The geomagnetic effects also produce a characteristic azimuthal dependence of the $\nu$--fluxes, related to the well known east--west effect, that should be detectable in neutrino experiments of sufficiently large mass. We discuss these effects quantitatively. Because the azimuthal dependence is in first order independent of any oscillation effect, it is a useful diagnostic tool for studying possible systematic effects in the search for neutrino oscillations.
• ### Pathlength distribution of atmospheric neutrinos(astro-ph/9708146)

Aug. 15, 1997 astro-ph
We present the distribution of the production heights of atmsopheric neutrinos as a function of zenith angle and neutrino energy. The distributions can be used as the input for evaluation of neutrino propagation under various hypotheses for neutrino flavor oscillations.
• ### Gamma-ray production in supernova remnants(astro-ph/9609044)

July 23, 1997 astro-ph
The bulk of the cosmic rays up to about 100 TeV are thought to be accelerated by the 1st order Fermi mechanism at supernova shocks, producing a power-law spectrum. Both electrons and protons should be accelerated, but their ratio on acceleration is not well known. Recently, the EGRET instrument on the Compton Gamma Ray Observatory has observed supernova remnants IC 443 and $\gamma$ Cygni at GeV energies. On the assumption that the observed gamma-rays are produced by accelerated particles in the remnants (rather than, for example, from a central compact object) we model the contributions due to pion production, bremsstrahlung, and inverse Compton scattering on the cosmic microwave, diffuse galactic radiation, and locally produced radiation fields. In the case of the same spectral index for both electrons and nuclei, and a cut-off at 80 TeV, we find that a spectral index of accelerated particles close to 2.4, and a ratio of electrons to protons in the range 0.2 to 0.3, gives a good fit to the observed spectra. For lower cut-off energies flatter spectra are possible. We also investigate the case where the electron spectrum is steeper than that of nuclei. We discuss the implications of our results for observations at air shower energies, and for the propagation of cosmic rays.
• ### Some Comments on Models of Hadronic Interactions at Air Shower Energies(hep-ph/9611352)

Nov. 18, 1996 hep-ph
Several models of minimum-bias hadronic interactions at ultra-high energy that have been used for calculations of air showers share essential common features. In this talk I review these common elements and discuss some consequences. I concentrate on properties of hadron-nucleus interactions, and I use mean depth of shower maximum as a function of primary energy to illustrate my main points. I will contrast these models with models that use a more naive treatment of hadronic interactions in nuclei but which have been successfully used to interpret measurements of depth of shower maximum.
• ### Review of Atmospheric Neutrinos(hep-ph/9611301)

Nov. 11, 1996 hep-ph
In this talk I review measurements and calculations of the flux of neutrinos produced by interactions of cosmic rays in the atmosphere. The main reason for interest in this subject is the apparent anomaly between the predicted and the observed ratio of $\nu_e$ to $\nu_\mu$. With the advent of Super-Kamiokande, we are on the threshold of an order of magnitude increase in the amount of data available to study the problem. My goal in this talk is to describe the current status of the subject, both for contained events and for neutrino-induced upward muons.
• ### Comparison of Atmospheric Neutrino Flux Calculations at Low Energies(hep-ph/9608253)

Aug. 7, 1996 hep-ph
We compare several different calculations of the atmospheric neutrino flux in the energy range relevant for contained neutrino interactions, and we identify the major sources of difference among the calculations. We find nothing that would affect the predicted ratio of $\nu_e/\nu_\mu$, which is nearly the same in all calculations. Significant differences in normalization arise primarily from different treatment of pion production by interactions of protons in the atmosphere. Different assumptions about the primary spectrum and treatment of the geomagnetic field are also of some importance.
• ### Atmospheric neutrino flux above 1 GeV(hep-ph/9509423)

Sept. 29, 1995 hep-ph
In this paper we extend an earlier calculation of the flux of atmospheric neutrinos to higher energy. The earlier calculation of the neutrino flux below 3 GeV has been used for calculation of the rate of contained neutrino interactions in deep underground detectors. The fluxes are needed up to neutrino energies of 10 TeV to calculate the expected rate of neutrino-induced muons passing into and through large, deep detectors. We compare our results with several other calculations, and we evaluate the uncertainty in the rate of neutrino-induced muons due to uncertainties in the neutrino flux.
• ### Calculated Electron Fluxes at Airplane Altitudes(hep-ph/9304215)

April 5, 1993 astro-ph, hep-ph
A precision measurement of atmospheric electron fluxes has been performed on a Japanese commercial airliner (Enomoto, {\it et al.}, 1991). We have performed a monte carlo calculation of the cosmic ray secondary electron fluxes expected in this experiment. The monte carlo uses the hadronic portion of our neutrino flux cascade program combined with the electromagnetic cascade portion of the CERN library program GEANT. Our results give good agreement with the data, provided we boost the overall normalization of the primary cosmic ray flux by 12\% over the normalization used in the neutrino flux calculation.
• ### Cosmic Rays IV. The spectrum and chemical composition above $10^4$ GeV(astro-ph/9303006)

March 11, 1993 astro-ph
Using the concept developed in earlier papers, that the cosmic rays originate in three different main sites, a) the normal supernova explosions into the interstellar medium, b) the supernova explosions into a stellar wind, and c) powerful radio galaxies, we demonstrate in this paper that the spectrum and chemical abundances above $10^4$ GeV can be well understood. Using existing data on the chemical composition of cosmic rays near TeV energies as a constraint, we adjust the parameters of the model to fit the shower size data from the Akeno experiment; this is necessary since the derivation of an all particle spectrum involves an assumption about the chemical composition of the cosmic rays and so we have to fit the shower size data first and then derive the all particle spectrum. We present a successful fit to the shower size data which allows us to draw three main conclusions: a) For most of the energy range above $10^5$ GeV the wind explosions can account for both chemical composition and spectrum including the knee feature, b) the highest particle energies required from the stellar wind explosions imply a magnetic field in the preexisting stellar wind of at least $3$ Gauss at a fiducial distance of $10^{14}$ cm, c) the chemical abundances above $10^5$ GeV are dominated by heavy nuclei such as Neon and higher. The bump observed in the all particle spectrum below the knee gets weakened with the proper treatment of the energy dependent chemical composition. For high energies we obtain an estimate of the extragalactic proton flux.