We argue that ultrahigh energy cosmic ray collisions in the Earth atmosphere can probe the strange quark density of the nucleon. These collisions have center-of-mass energies \agt 10^{4.6} A GeV, where A \geq 14 is the nuclear baryon number. We hypothesize the formation of a deconfined thermal fireball which undergoes a sudden hadronization. At production the fireball has a very high matter density and consists of gluons and two flavors of light quarks (u,d). Because the fireball is formed in the baryon-rich projectile fragmentation region, the high baryochemical potential damps the production of u \bar u and d \bar d pairs, resulting in gluon fragmentation mainly into s \bar s. The strange quarks then become much more abundant and upon hadronization the relative density of strange hadrons is significantly enhanced over that resulting from a hadron gas. Assuming the momentum distribution functions can be approximated by Fermi-Dirac and Bose-Einstein statistics, we estimate a kaon-to-pion ratio of about 3 and expect a similar (total) baryon-to-pion ratio. We show that, if this were the case, the excess of strange hadrons would suppress the fraction of energy which is transferred to decaying \pi^0's by about 20\%, yielding a \sim 40\% enhancement of the muon content in atmospheric cascades, in agreement with recent data reported by the Pierre Auger Collaboration.

The Antarctic Impulsive Transient Antenna (ANITA) experiment has observed two air shower events with energy $\sim 500~{\rm PeV}$ emerging from the Earth with exit angles $\sim 30^\circ$ above the horizon. As was immediately noted by the ANITA Collaboration, these events (in principle) could originate in the atmospheric decay of an upgoing $\tau$-lepton produced through a charged current interaction of a $\nu_\tau$ inside the Earth. However, the relatively steep arrival angles of these perplexing events create tension with the standard model (SM) neutrino-nucleon interaction cross section. Deepening the conundrum, the IceCube neutrino telescope and the Pierre Auger Observatory with substantially larger exposures to cosmic $\nu_\tau$'s in this energy range have not observed any events. This lack of observation implies that the messenger particle (MP) giving rise to ANITA events must produce an air shower event rate at least a factor of 40 larger than that produced by a flux of $\tau$-neutrinos to avoid conflicts with the upper limits reported by the IceCube and the Pierre Auger collaborations. In addition, the sensitivity of ANITA to MP-induced events must be comparable to or larger than those of IceCube and Auger to avoid conflict with the non-observation of any signal at these facilities. Beyond SM interpretations of ANITA events can be classified according to whether the MPs: (i) live inside the Earth, (ii) originate in neutrino-nucleon collisions inside the Earth, (iii) come from cosmological distances. In this communication we investigate the positive and negative facets of these three classes of models.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

The Pierre Auger Collaboration has reported an excess in the number of muons of a few tens of percent over expectations computed using extrapolation of hadronic interaction models tuned to accommodate LHC data. Very recently, we proposed an explanation for the muon excess assuming the formation of a deconfined quark matter (fireball) state in central collisions of ultrarelativistic cosmic rays with air nuclei. At the first stage of its evolution the fireball contains gluons as well as $u$ and $d$ quarks. The very high baryochemical potential inhibits gluons from fragmenting into $u \bar u$ and $d \bar d$, and so they fragment predominantly into $s \bar s$ pairs. In the hadronization which follows this leads to the strong suppression of pions and hence photons, but allows heavy hadrons to be emitted carrying away strangeness. In this manner, the extreme imbalance of hadron to photon content provides a way to enhance the muon content of the air shower. In this communication we study theoretical systematics from hadronic interaction models used to describe the cascades of secondary particles produced in the fireball explosion. We study the predictions of one of the leading LHC-tuned models QGSJET II-04 considered in the Auger analysis.

We explain the two upgoing ultra-high energy shower events observed by ANITA as arising from the decay in the Earth's interior of the quasi-stable dark matter candidate in the CPT symmetric universe. The dark matter particle is a 480 PeV right-handed neutrino that decays into a Higgs boson and a light Majorana neutrino. The latter interacts in the Earth's crust to produce a tau lepton that in turn initiates an atmospheric upgoing shower. The fact that both events emerge at the same angle from the Antarctic ice-cap suggests an atypical dark matter density distribution in the Earth.

We argue that if ultrahigh-energy (E > 10^10 GeV) cosmic rays are heavy nuclei (as indicated by existing data), then the pointing of cosmic rays to their nearest extragalactic sources is expected for 10^10.6 < E/GeV < 10^11. This is because for a nucleus of charge Ze and baryon number A, the bending of the cosmic ray decreases as Z/E with rising energy, so that pointing to nearby sources becomes possible at highest energies. In addition, the maximum energy of acceleration capability of the sources grows linearly in Z, while the energy loss per distance traveled decreases with increasing A. Each of these two points tend to favor heavy nuclei at the highest energies. The traditional bi-dimensional analyses, which simultaneously reproduce Auger data on the spectrum and nuclear composition, may not be capable of incorporating the relative importance of these phenomena. In this paper we propose a multi-dimensional reconstruction of the individual emission spectra (in E, direction, and cross-correlation with nearby putative sources) to study the hypothesis that primaries are heavy nuclei subject to GZK photo-disintegration, and to determine the nature of the extragalactic sources. More specifically, we propose to combine information on nuclear composition and arrival direction to associate a potential clustering of events with a 3-dimensional position in the sky. Actually, both the source distance and maximum emission energy can be obtained through a multi-parameter likelihood analysis to accommodate the observed nuclear composition of each individual event in the cluster. We show that nucleus-emitting-sources exhibit a cepa stratis structure on Earth which could be pealed off by future space-missions, such as POEMMA. We also show that metal-rich starburst galaxies are highly-plausible candidate sources, and we use them as an explicit example of our proposed multi-dimensional analysis.

High energy neutrinos have been detected by IceCube, but their origin remains a mystery. Determining the sources of this flux is a crucial first step towards multi-messenger studies. In this work we systematically compare two classes of sources with the data: Galactic and extragalactic. We assume that the neutrino sources are distributed according to a class of Galactic models. We build a likelihood function on an event by event basis including energy, event topology, absorption, and direction information. We present the probability that each high energy event with deposited energy $E_{\rm dep}>60$ TeV in the HESE sample is Galactic, extragalactic, or background. For Galactic models considered the Galactic fraction of the astrophysical flux has a best fit value of $1.3\%$ and is $<9.5\%$ at 90\% CL. A zero Galactic flux is allowed at $<1\sigma$.

Many dark matter interaction types lead to annihilation processes which suffer from p-wave suppression or helicity suppression, rendering them subdominant to unsuppressed s-wave processes. We demonstrate that the natural inclusion of dark initial state radiation can open an unsuppressed s-wave annihilation channel, and thus provide the dominant dark matter annihilation process for particular interaction types. We illustrate this effect with the bremsstrahlung of a dark spin-0 or dark spin-1 particle from fermionic dark matter, $\overline{\chi}\chi\rightarrow \overline{f}f\phi$ or $\overline{f}fZ'$. The dark initial state radiation process, despite having a 3-body final state, proceeds at the same order in the new physics scale $\Lambda$ as the annihilation to the 2-body final state $\overline{\chi}\chi\rightarrow \overline{f}f$. This is lower order in $\Lambda$ than the well-studied lifting of helicity suppression via Standard Model final state radiation, or virtual internal bremsstrahlung. This dark bremsstrahlung process should influence LHC and indirect detection searches for dark matter.

The announcement by the IceCube Collaboration of the observation of 53 astrophysical neutrino candidates in the energy range 0.03 \alt E_\nu/PeV \alt 2 has been greeted with a great deal of justified excitement. Herein we provide fits of single and a broken power-law energy-spectra to these high-energy starting events (HESEs). By comparing our statistical results from fits to (background-free) shower HESE data with the spectral shape of muon neutrinos recently reported by the IceCube Collaboration, we show that there is (3 \sigma) evidence for a break in the spectrum of astrophysical neutrinos. After that we use the fitted result to predict the rate of Glashow events (in the ~ 6.3 PeV region) and double-bang tau neutrino events (in the PeV region) just at the threshold of IceCube detection.

We study the Glashow resonance $\bar{\nu}_e + e^- \rightarrow W^- \rightarrow$ hadrons at 6.3 PeV as diagnostic of the production processes of ultra-high energy neutrinos. The focus lies on describing the physics of neutrino production from pion decay as accurate as possible by including the kinematics of weak decays and Monte Carlo simulations of pp and p$\gamma$ interactions. We discuss optically thick (to photohadronic interactions) sources, sources of cosmic ray nuclei and muon damped sources. Even in the proposed upgrade IceCube-Gen2, a discrimination of scenarios such as pp versus p$\gamma$ is extremely challenging under realistic assumptions. Nonetheless, the Glashow resonance can serve as a smoking gun signature of neutrino production from photohadronic (A$\gamma$) interactions of heavier nuclei, as the expected Glashow event rate exceeds that of pp interactions. We finally quantify the exposures for which the non-observation of Glashow events exerts pressure on certain scenarios.

Cosmological parameters deduced from the Planck measurements of anisotropies in the cosmic microwave background are at some tension with direct astronomical measurements of various parameters at low redshifts. Very recently, it has been conjectured that this discrepancy can be reconciled if a certain fraction of dark matter is unstable and decays between recombination and the present epoch. Herein we show that if the superheavy relics have a branching into neutrinos B (X \to \nu \bar \nu) \sim 3 \times 10^{-9}, then this scenario can also accommodate the recently discovered extraterrestrial flux of neutrinos, relaxing the tension between IceCube results and Fermi LAT data. The model is fully predictive and can be confronted with future IceCube data. We demonstrate that in 10 years of observation IceCube will be able to distinguish the mono-energetic signal from X decay at the 3\sigma level. In a few years of data taking with the upgraded IceCube-Gen2 enough statistics will be gathered to elucidate the dark matter--neutrino connection at the 5\sigma level.

Big Bang Nucleosynthesis imposes stringent bounds on light sterile neutrinos mixing with the active flavors. Here we discuss how altered dispersion relations can weaken such bounds and allow compatibility of new sterile neutrino degrees of freedom with a successful generation of the light elements in the early Universe.

In the last few years ARCADE 2, combined with older experiments, has detected an additional radio background, measured as a temperature and ranging in frequency from 22 MHz to 10 GHz, not accounted for by known radio sources and the cosmic microwave background. One type of source which has not been considered in the radio background is that of fast transients (those with event times much less than the observing time). We present a simple estimate, and a more detailed calculation, for the contribution of radio transients to the diffuse background. As a timely example, we estimate the contribution from the recently-discovered fast radio bursts (FRBs). Although their contribution is likely 6 or 7 orders of magnitude too small (though there are large uncertainties in FRB parameters) to account for the ARCADE~2 excess, our development is general and so can be applied to any fast transient sources, discovered or yet to be discovered. We estimate parameter values necessary for transient sources to noticeably contribute to the radio background.

We study the vacuum stability of a minimal Higgs portal model in which the standard model (SM) particle spectrum is extended to include one complex scalar field and one Dirac fermion. These new fields are singlets under the SM gauge group and are charged under a global U(1) symmetry. Breaking of this U(1) symmetry results in a massless Goldstone boson, a massive CP-even scalar, and splits the Dirac fermion into two new mass-eigenstates, corresponding to Majorana fermions. The lightest Majorana fermion (w) is absolutely stable, providing a plausible dark matter (DM) candidate. We show that interactions between the Higgs sector and the lightest Majorana fermion which are strong enough to yield a thermal relic abundance consistent with observation can easily destabilize the electroweak vacuum or drive the theory into a non-perturbative regime at an energy scale well below the Planck mass. However, we also demonstrate that there is a region of the parameter space which develops a stable vacuum (up to the Planck scale), satisfies the relic abundance, and is in agreement with direct DM searches. Such an interesting region of the parameter space corresponds to DM masses 350 GeV \alt m_w \alt 1 TeV. The region of interest is within reach of second generation DM direct detection experiments.

Now that PeV neutrinos have been discovered by IceCube, we optimistically entertain the possibility that neutrinos with energy above 100 PeV exist. We evaluate the dependence of event rates of such neutrinos on the neutrino-nucleon cross section at observatories that detect particles, atmospheric fluorescence, or Cherenkov radiation, initiated by neutrino interactions. We consider how (i) a simple scaling of the total standard model neutrino-nucleon cross section, (ii) a new elastic neutral current interaction, and (iii) a new completely inelastic interaction, individually impact event rates.

Determining anisotropies in the arrival directions of cosmic rays at the highest energy is an important task in astrophysics. It is common and useful to partition the sky into spherical harmonics as a measure of anisotropy. The two lowest nontrivial spherical harmonics, the dipole and the quadrupole, are of particular interest, since these distributions encapsulate a dominant single source and a plane of sources, as well as offering relatively high statistics. The best experiments for the detection of ultra high energy cosmic rays currently are all ground-based, with highly nonuniform exposures on the sky resulting from the fixed experimental locations on the earth. This nonuniform exposure increases the complexity and error in inferring anisotropies. It turns out that there is an optimal latitude for an experiment at which nonuniform exposure does not diminish the inference of the quadrupole moment. We derive the optimal latitude, and find that (presumably by a fortuitous coincidence) this optimal latitude runs through the largest cosmic ray experiment, Pierre Auger Observatory (PAO) in the Southern Hemisphere, and close to the largest cosmic ray experiment in the Northern Hemisphere, Telescope Array (TA). Consequently, assuming a quadrupole distribution, PAO and TA can reconstruct the cosmic ray quadrupole distribution to a high precision, without concern for their partial sky exposure.

The two main advantages of space-based observation of extreme energy ($\gtrsim5\times10^{19}$ eV) cosmic rays (EECRs) over ground based observatories are the increased field of view and the full-sky coverage with nearly uniform systematics across the entire sky. The former guarantees increased statistics, whereas the latter enables a clean partitioning of the sky into spherical harmonics. The discovery of anisotropies would help to identify the long sought origin of EECRs. We begin an investigation of the reach of a full-sky space-based experiment such as EUSO to detect anisotropies in the extreme-energy cosmic-ray sky compared to ground based partial-sky experiments such as the Pierre Auger Observatory and Telescope Array. The technique is explained here, and simulations for a Universe with just two nonzero multipoles, monopole plus either dipole or quadrupole, are presented. These simulations quantify the advantages of space-based, all-sky coverage.

Over cosmic distances, astrophysical neutrino oscillations average out to a classical flavor propagation matrix $\mathscr{P}$. Thus, flavor ratios injected at the cosmic source $W_e,W_\mu,W_\tau$ evolve to flavor ratios at Earthly detectors $w_e,w_\mu,w_\tau$ according to $\vec{w}=\mathscr{P} \vec{W}$. The unitary constraint reduces the Euclidean octant to a "flavor triangle". We prove a theorem that the area of the Earthly flavor triangle is proportional to Det$(\mathscr{P})$. One more constraint would further reduce the dimensionality of the flavor triangle at Earth (two) to a line (one). We discuss four motivated such constraints. The first is the possibility of a vanishing determinant for $\mathscr{P}$. We give a formula for a unique $\delta(\theta_{ij}$'s) that yields the vanishing determinant. Next we consider the thinness of the Earthly flavor triangle. We relate this thinness to the small deviations of the two angles $\theta_{32}$ and $\theta_{13}$ from maximal mixing and zero, respectively. Then we consider the confusion resulting from the tau neutrino decay topologies, which are showers at low energy, "double-bang" showers in the PeV range, and a mixture of showers and tracks at even higher energies. We examine the simple low-energy regime, where there are just two topologies, $w_{\rm shower}=w_e+w_\tau$ and $w_{\rm track}=w_\mu$. We apply the statistical uncertainty to be expected from IceCube to this model. Finally, we consider ramifications of the expected lack of $\nu_\tau$ injection at cosmic sources. In particular, this constraint reduces the Earthly triangle to a boundary line of the triangle. Some tests of this "no $\nu_\tau$ injection" hypothesis are given.

The Glashow resonance at E_\nu=6.3 PeV is a measure of the \bar\nu_e content of the astrophysical neutrino flux. The fractional \bar\nu_e content depends on the neutrino production model at the cosmic neutrino source, and the environment at the source. Thus, the strength of the Glashow resonance event rate is a potential window into astrophysical sources. We quantify the "Glashow resonometer" and comment on the significance that no Glashow events are observed in the IceCube three-year data.

There may be a high-energy cutoff of neutrino events in IceCube data. In particular, IceCube does not observe either continuum events above 2 PeV, or the Standard Model Glashow-resonance events expected at 6.3 PeV. There are also no higher energy neutrino signatures in the ANITA and Auger experiments. This absence of high-energy neutrino events motivates a fundamental restriction on neutrino energies above a few PeV. We postulate a simple scenario to terminate the neutrino spectrum that is Lorentz-invariance violating, but with a limiting neutrino velocity that is always smaller than the speed of light. If the limiting velocity of the neutrino applies also to its associated charged lepton, then a significant consequence is that the two-body decay modes of the charged pion are forbidden above two times the maximum neutrino energy, while the radiative decay modes are suppressed at higher energies. Such stabilized pions may serve as cosmic ray primaries.

We present a relationship, E_\nu^{max} = m_{\nu} M_{Planck}/M_{weak}, among the highest observed neutrino energy (~PeV) and the neutrino mass, the weak scale, and the Planck energy. We then discuss some tests of this relationship, and present some theoretical constructs which motivate the relationship. It is possible that all massive particles are subject to maximum energies given by similar relationships, although only the neutrino seems able to offer interesting phenomenology. We discuss implications which include no neutrino detections at energies greater than PeV, and changes in expectations for the highest energy cosmic rays. A virtue of this hypothesis is that it is easily invalidated should neutrinos be observed with energies much great than the PeV scale. An almost inescapable implication is that Lorentz Invariance is a low energy principle, yet it appears that violation may be only observable in high-energy astrophysical neutrinos.

We investigate the neutrino - cosmic ray connection for sources in the Galaxy in terms of two observables: the shape of the energy spectrum and the distribution of arrival directions. We also study the associated gamma ray emission from these sources.

Very recently the IceCube Collaboration has reported an observation of 28 neutrino candidates with energies between 50 TeV and 2 PeV, constituting a 4.1$\sigma$ excess compared to the atmospheric background. In this article we investigate the compatibility between the data and a hypothesized unbroken power-law neutrino spectrum for various values of spectral index \Gamma >= 2. We show that \Gamma ~ 2.3 is consistent at the ~ 1.5$\sigma$ level with the observed events up to 2 PeV and to the null observation of events at higher energies. We then assume that the sources of this unbroken spectrum are Galactic, and deduce (i) an energy-transfer fraction from parent protons to pions, and (ii) a way of discriminating among models which have been put forth to explain the "knee" and "ankle" features of the cosmic ray spectrum. Future IceCube data will test the unbroken power law hypothesis and provide a multi-messenger approach to explaining features of the cosmic ray spectrum, including the transition from Galactic to extragalactic dominance.

We discuss a number of experimental constraints on Weinberg's Higgs portal model. In this framework, the standard model (SM) particle spectrum is extended to include one complex scalar field S and one Dirac fermion \psi. These new fields are singlets under the SM gauge group and are charged under a global U(1) symmetry. Breaking of this U(1) symmetry results in a massless Goldstone boson \alpha and a massive CP-even scalar r, and splits the Dirac fermion into two new mass-eigenstates \psi_\pm, corresponding to Majorana fermions. The interest on such a minimal SM extension is twofold. On the one hand, if the Goldstone bosons are in thermal equilibrium with SM particles until the era of muon annihilation their contribution to the effective number of neutrino species can explain the hints from cosmological observations of extra relativistic degrees of freedom at the epoch of last scattering. On the other hand, the lightest Majorana fermion \psi_- provides a plausible dark matter candidate. Mixing of r with the Higgs doublet \phi is characterized by the mass of hidden scalar m_h and the mixing angle \theta. We constrain this parameter space using a variety of experimental data, including heavy meson decays with missing energy, the invisible Higgs width, and direct dark matter searches. We show that different experimental results compress the allowed parameter space in complementary ways, covering a large range of \psi_- masses (5 GeV \alt m_- \alt 100 GeV). Though current results narrow the parameter space significantly (for the mass range of interest, \theta \alt 10^{-3} to 10^{-4}), there is still room for discovery (\alpha decoupling at the muon annihilation era requires \theta \agt 10^{-5} to 10^{-4}). In the near future, measurements from ATLAS, CMS, LHCb, NA62, XENON1T, LUX, and CDMSlite will probe nearly the full parameter space.

The announcement by the IceCube Collaboration of the observation of 28 cosmic neutrino candidates has been greeted with a great deal of justified excitement. The data reported so far depart by 4.3\sigma from the expected atmospheric neutrino background, which raises the obvious question: "Where in the Cosmos are these neutrinos coming from?" We review the many possibilities which have been explored in the literature to address this question, including origins at either Galactic or extragalactic celestial objects. For completeness, we also briefly discuss new physics processes which may either explain or be constrained by IceCube data.

The two main advantages of space-based observation of extreme-energy ($\gtrsim 10^{19}$~eV) cosmic-rays (EECRs) over ground-based observatories are the increased field of view, and the all-sky coverage with nearly uniform systematics of an orbiting observatory. The former guarantees increased statistics, whereas the latter enables a partitioning of the sky into spherical harmonics. We have begun an investigation, using the spherical harmonic technique, of the reach of \J\ into potential anisotropies in the extreme-energy cosmic-ray sky-map. The technique is explained here, and simulations are presented. The discovery of anisotropies would help to identify the long-sought origin of EECRs.

Report of the CF6 Working Group at Snowmass 2013. Topics addressed include ultra-high energy cosmic rays, neutrinos, gamma rays, baryogenesis, and experiments probing the fundamental nature of spacetime.

Arguably, the most popular candidate for Dark Matter (DM) is a massive, stable, Majorana fermion. However, annihilation of Majorana DM to two fermions features a helicity-suppressed s-wave rate. The process radiating a gauge boson via electroweak (EW) and electromagnetic (EM) bremsstrahlung removes this s-wave suppression, and is likely to be the dominant modes of gauge-singlet Majorana DM annihilation. Given their enhanced annihilation rate with radiatedW and Z gauge bosons, and the subsequent dominant W/Z decays via hadronic channels, Majorana DM tends to produce more antiprotons than positrons. This result contrasts with observations, thereby presenting a challenge to model building with Majorana DM.

We summarize the remarkable recent progress in ultra-high energy cosmic ray physics and astronomy enabled by the current generation of cosmic ray observatories. We discuss the primary objectives for future measurements and describe the plans for near-term enhancements of existing experiments as well as the next generation of observatories.

We evaluate the exposure during nadir observations with JEM-EUSO, the Extreme Universe Space Observatory, on-board the Japanese Experiment Module of the International Space Station. Designed as a mission to explore the extreme energy Universe from space, JEM-EUSO will monitor the Earth's nighttime atmosphere to record the ultraviolet light from tracks generated by extensive air showers initiated by ultra-high energy cosmic rays. In the present work, we discuss the particularities of space-based observation and we compute the annual exposure in nadir observation. The results are based on studies of the expected trigger aperture and observational duty cycle, as well as, on the investigations of the effects of clouds and different types of background light. We show that the annual exposure is about one order of magnitude higher than those of the presently operating ground-based observatories.

Arguably, the most popular candidate for Dark Matter (DM) is a massive, stable, Majorana fermion. However, annihilation of Majorana DM to two fermions often features a helicity-suppressed s-wave rate. Radiating a gauge boson via electroweak (EW) and electromagnetic (EM) bremsstrahlung removes this s-wave suppression. The main purpose of this talk is to explain in some detail why the branching ratio to a fermion pair is likely suppressed while the decay to the pair plus a W/Z is not. In doing so, we investigate the general conditions for s-wave suppression and un-suppression using Fierz transformations and partial wave expansions. Suppression for the 2-to-2 process is sufficiently severe that the EW and EM bremsstrahlung are likely to be the dominant modes of gauge-singlet Majorana DM annihilation. We end this talk with a discussion of the challenge presented by space-based data for Majorana DM models, given that the enhanced rate to radiated W and Z gauge bosons and their dominant decay via hadronic channels tends to produce more anti-protons than are observed.

Reactor neutrino experiments have now observed a nonzero value for $\theta_{13}$ at $5\sigma$, and global fits to data imply a nonzero value above $10\sigma$. Nonzero values for $\theta_{13}$ and/or $\theta_{32}-\pi/4$ break a $\nu_\mu-\nu_\tau$ symmetry, which has qualitative as well as quantitative implications for the time-evolution of neutrino flavors. In particular, the large-distance flavor evolution matrix, non-invertible with $\nu_\mu-\nu_\tau$ symmetry, is now invertible. This means that measurements of neutrino flavor ratios at Earth can now be inverted to directly reveal the flavor ratios injected at cosmically distant sources. With the updated values of the three neutrino mixing angles, we obtain the inverted large-distance evolution matrix and use it to derive several phenomenological relations between the injection flavor ratios and the observable ratios at Earth. Taking the three popular injection models as examples, we also exhibit the shift of Earthly observed flavor ratios from the corresponding values in the case with $\nu_\mu-\nu_\tau$ symmetry.

This document contains a summary of the workshop which took place on 22 - 24 February 2012 at the Kavli Institute of Cosmological Physics in the University of Chicago. The goal of the workshop was to discuss the physics reach of the JEM-EUSO mission and how best to implement a global ground based calibration system for the instrument to realize the physics goal of unveiling the origin of the highest energy cosmic rays.

We examine observational signatures of dark matter annihilation in the Milky Way arising from electroweak bremsstrahlung contributions to the annihilation cross section. It has been known for some time that photon bremsstrahlung may significantly boost DM annihilation yields. Recently, we have shown that electroweak bremsstrahlung of W and Z gauge bosons can be the dominant annihilation channel in some popular models with helicity-suppressed 2 --> 2 annihilation. W/Z-bremsstrahlung is particularly interesting because the gauge bosons produced via annihilation subsequently decay to produce large correlated fluxes of electrons, positrons, neutrinos, hadrons (including antiprotons) and gamma rays, which are all of importance in indirect dark matter searches. Here we calculate the spectra of stable annihilation products produced via gamma/W/Z-bremsstrahlung. After modifying the fluxes to account for the propagation through the Galaxy, we set upper bounds on the annihilation cross section via a comparison with observational data. We show that stringent cosmic ray antiproton limits preclude a sizable dark matter contribution to observed cosmic ray positron fluxes in the class of models for which the bremsstrahlung processes dominate.

We revisit the calculation of electroweak bremsstrahlung contributions to dark matter annihilation. Dark matter annihilation to leptons is necessarily accompanied by electroweak radiative corrections, in which a $W$ or $Z$ boson is also radiated. Significantly, while many dark matter models feature a helicity suppressed annihilation rate to fermions, bremsstrahlung process can remove this helicity suppression such that the branching ratios Br($\ell \nu W $), Br($\ell^+\ell^-Z$), and Br($\bar\nu \nu Z$) dominate over Br($\ell^+\ell^-$) and Br($\bar\nu \nu$). We find this is most significant in the limit where the dark matter mass is nearly degenerate with the mass of the boson which mediates the annihilation process. Electroweak bremsstrahlung has important phenomenological consequences both for the magnitude of the total dark matter annihilation cross section and for the character of the astrophysical signals for indirect detection. Given that the $W$ and $Z$ gauge bosons decay dominantly via hadronic channels, it is impossible to produce final state leptons without accompanying protons, antiprotons, and gamma rays.

Contributions of the JEM-EUSO Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.

Dark matter annihilation to leptons, $\chi \chi \rightarrow l \bar{l}$, is necessarily accompanied by electroweak radiative corrections, in which a W or Z boson is radiated from a final state particle. Given that the W and Z gauge bosons decay dominantly via hadronic channels, it is thus impossible to produce final state leptons without accompanying protons, antiprotons, and gamma rays. Significantly, while many dark matter models feature a helicity suppressed annihilation rate to fermions, radiating a massive gauge boson from a final state fermion removes this helicity suppression, such that the branching ratios Br(l \nu W), Br(l^+l^- Z), and Br(\nu\nubar Z) dominate over Br(l\bar{l}). W/Z-bremsstrahlung thus allows indirect detection of many WIMP models that would otherwise be helicity-suppressed, or v^2 suppressed. Antiprotons and even anti-deuterons become consequential final state particles. This is an important result for future DM searches. We discuss the implications of W/Z-bremsstrahlung for "leptonic" DM models which aim to fit recent cosmic ray positron and antiproton data.

We propose that neutron emission from Cen A dominates the cosmic ray sky at the high end of the spectrum. Neutrons that are able to decay generate proton diffusion fronts, whereas those that survive decay produce a spike in the direction of the source. We use recent data reported by the Pierre Auger Collaboration to normalize the injection spectrum and estimate the required luminosity in cosmic rays. We find that such a luminosity, L_{CR} ~ 5 x 10^{40} erg/s, is considerably smaller than the bolometric luminosity of Cen A, L_{bol} ~ 10^{43} erg/s. We compute the incoming current flux density as viewed by an observer on Earth and show that the anisotropy amplitude is in agreement with data at the 1\sigma level. Regardless of the underlying source model, our results indicate that after a decade of data taking the Pierre Auger Observatory will be able to test our proposal.

In this essay we extend the standard discussion of neutrino oscillations to astrophysical neutrinos propagating through expanding space. This extension introduces a new cosmological parameter $I$ into the oscillation phase. The new parameter records cosmic history in much the same manner as the redshift z or the apparent luminosity D_L. Measuring $I$ through neutrino oscillations could help determine cosmological parameters and discriminate among different cosmologies.

The MAGIC Collaboration has provided new observational data pertaining to the TeV J2032+4130 gamma-ray source (within the Cygnus OB2 region), for energies E_gamma >400 GeV. It is then appropriate to update the impact of these data on gamma-ray production mechanisms in stellar associations. We consider two mechanisms of gamma-ray emission, pion production and decay (PION) and photo-excitation of high-energy nuclei followed by prompt photo-emission from the daughter nuclei (A*). We find that while the data can be accommodated with either scenario, the A* features a spectral bump, corresponding to the threshold for exciting the Giant Dipole Resonance, which can serve to discriminate between them. We comment on neutrino emission and detection from the region if the PION and/or A* processes are operative. We also touch on the implications for this analysis of future Fermi and Cerenkov Telescope Array data.

JEM-EUSO, on board of the Japanese Exploration Module of the International Space Station, is being proposed as the first space observatory devoted to UHECR. Its privileged position at 430 km above the Earth surface, combined with a large field of view, innovative optics and a high efficiency focal surface, results in an unprecedented exposure which significantly surpasses that of the largest ground observatories. The large number of events expected above the GZK threshold for photo-pion production by protons will allow the directional identification of individual sources and the determination of their spectra, i.e., doing astronomy and astrophysics through the particle channel. Similar goals can be achieved in the case of light UHECR nuclei. Furthermore, the atmospheric target volume ($\sim 10^{12}$ ton) makes the possibility of neutrino observation a highlight of the mission. Other exploratory objectives include the detection of extreme energy gammas and the study of Galactic magnetic fields as well as global observations of the earth's atmosphere, including clouds, night-glows, plasma discharges, and meteors. In this contribution we will describe the scientific objectives of JEM-EUSO.

In asymmetrically warped spacetimes different warp factors are assigned to space and to time. We discuss causality properties of these warped brane universes and argue that scenarios with two extra dimensions may allow for timelike curves which can be closed via paths in the extra-dimensional bulk. In particular, necessary and sufficient conditions on the metric for the existence of closed timelike curves are presented. We find a six-dimensional warped metric which satisfies the CTC conditions, and where the null, weak and dominant energy conditions are satisfied on the brane (although only the former remains satisfied in the bulk). Such scenarios are interesting, since they open the possibility of experimentally testing the chronology protection conjecture by manipulating on our brane initial conditions of gravitons or hypothetical gauge-singlet fermions (sterile neutrinos) which then propagate in the extra dimensions.

A conservative upper bound on the total dark matter (DM) annihilation rate can be obtained by constraining the appearance rate of the annihilation products which are hardest to detect. The production of neutrinos, via the process $\chi \chi \to \bar\nu \nu $, has thus been used to set a strong general bound on the dark matter annihilation rate. However, Standard Model radiative corrections to this process will inevitably produce photons which may be easier to detect. We present an explicit calculation of the branching ratios for the electroweak bremsstrahlung processes $\chi \chi \to \bar\nu \nu Z$ and $\chi \chi \to \bar\nu e W$. These modes inevitably lead to electromagnetic showers and further constraints on the DM annihilation cross-section. In addition to annihilation, our calculations are also applicable to the case of dark matter decay.

The effect of quintessence perturbations on the ISW effect is studied for a mixed dynamical scalar field dark energy (DDE) and pressureless perfect fluid dark matter. A new and general methodology is developed to track the growth of the perturbations, which uses only the equation of state (EoS) parameter $w_{\rm DDE} (z) \equiv p_{\rm DDE}/\rho_{\rm DDE}$ of the scalar field DDE, and the initial values of the the relative entropy perturbation (between the matter and DDE) and the intrinsic entropy perturbation of the scalar field DDE as inputs. We also derive a relation between the rest frame sound speed $\hat{c}_{s,{\rm DDE}}^2$ of an arbitrary DDE component and its EoS $w_{\rm DDE} (z)$. We show that the ISW signal differs from that expected in a $\Lambda$CDM cosmology by as much as +20% to -80% for parameterizations of $w_{\rm DDE}$ consistent with SNIa data, and about $\pm$ 20% for parameterizations of $w_{\rm DDE}$ consistent with SNIa+CMB+BAO data, at 95% confidence. Our results indicate that, at least in principle, the ISW effect can be used to phenomenologically distinguish a cosmological constant from DDE.

Observational limits on the high-energy neutrino background have been used to place general constraints on dark matter that annihilates only into standard model particles. Dark matter particles that annihilate into neutrinos will also inevitably branch into electromagnetic final states through higher-order tree and loop diagrams that give rise to charged leptons, and these charged particles can transfer their energy into photons via synchrotron radiation or inverse Compton scattering. In the context of effective field theory, we calculate the loop-induced branching ratio to charged leptons and show that it is generally quite large, typically >1%, when the scale of the dark matter mass exceeds the electroweak scale, M_W. For a branching fraction >3%, the synchrotron radiation bounds on dark matter annihilation are currently stronger than the corresponding neutrino bounds in the interesting mass range from 100 GeV to 1 TeV. For dark matter masses below M_W, our work provides a plausible framework for the construction of a model for "neutrinos only" dark matter annihilations.

Current experimental data on neutrino mixing are very well described by TriBiMaximal mixing. Accordingly, any phenomenological parametrization of the MNSP matrix must build upon TriBiMaximal mixing. We propose one particularly natural parametrization, which we call "TriMinimal". The three small deviations of the PDG angles from their TriBiMaximal values, and the PDG phase, parametrize the TriMinimal mixing matrix. As an important example of the utility of this new parametrization, we present the simple resulting expressions for the flavor-mixing probabilities of atmospheric and astrophysical neutrinos. As no foreseeable experiment will be sensitive to more than second order in the small parameters, we expand these flavor probabilities to second order.

We discuss flavor-mixing probabilities and flavor ratios of high energy astrophysical neutrinos. In the first part of this paper, we expand the neutrino flavor-fluxes in terms of the small parameters U_{e3} and pi/4 - theta_{23}, and show that there are universal first and second order corrections. The second order term can exceed the first order term, and so should be included in any analytic study. We also investigate the probabilities and ratios after a further expansion around the tribimaximal value of sin^2 theta_{12} = 1/3. In the second part of the paper, we discuss implications of deviations of initial flavor ratios from the usually assumed, idealized flavor compositions for pion, muon-damped, and neutron beam sources, viz., (1 : 2 : 0), (0 : 1 : 0), and (1 : 0 : 0), respectively. We show that even small deviations have significant consequences for the observed flavor ratios at Earth. If initial flavor deviations are not taken into account in analyses, then false inferences for the values in the PMNS matrix elements (angles and phase) may result.

We consider the neutrino (and antineutrino) flavors arriving at Earth for neutrinos produced in the annihilation of weakly interacting massive particles (WIMPs) in the Sun's core. Solar-matter effects on the flavor propagation of the resulting $\agt$ GeV neutrinos are studied analytically within a density-matrix formalism. Matter effects, including mass-state level-crossings, influence the flavor fluxes considerably. The exposition herein is somewhat pedagogical, in that it starts with adiabatic evolution of single flavors from the Sun's center, with $\theta_{13}$ set to zero, and progresses to fully realistic processing of the flavor ratios expected in WIMP decay, from the Sun's core to the Earth. In the fully realistic calculation, non-adiabatic level-crossing is included, as are possible nonzero values for $\theta_{13}$ and the CP-violating phase $\delta$. Due to resonance enhancement in matter, nonzero values of $\theta_{13}$ even smaller than a degree can noticeably affect flavor propagation. Both normal and inverted neutrino-mass hierarchies are considered. Our main conclusion is that measuring flavor ratios (in addition to energy spectra) of $\agt$ GeV solar neutrinos can provide discrinination between WIMP models. In particular, we demonstrate the flavor differences at Earth for neutrinos from the two main classes of WIMP final states, namely $W^+ W^-$ and 95% $b \bar{b}$ + 5% $\tau^+\tau^-$. Conversely, if WIMP properties were to be learned from production in future accelerators, then the flavor ratios of $\agt$ GeV solar neutrinos might be useful for inferring $\theta_{13}$ and the mass hierarchy.

TeV gamma-rays can result from the photo-de-excitation of PeV cosmic ray nuclei after their parents have undergone photo-disintegration in an environment of ultraviolet photons. This process is proposed as a candidate explanation of the recently discovered HESS source at the edge of Westerlund 2. The UV background is provided by Lyman-alpha emission within the rich O and B stellar environment. The HESS flux results if there is efficient acceleration at the source of lower energy nuclei. The requirement that the Lorentz-boosted ultraviolet photons reach the Giant Dipole resonant energy (~ 20 MeV) implies a strong suppression of the gamma-ray spectrum compared to an E_\gamma^{-2} behavior at energies below about 1 TeV. This suppression is not apparent in the lowest-energy Westerlund 2 datum, but will be probed by the upcoming GLAST mission.

TeV gamma-rays may provide significant information about high energy astrophysical accelerators. Such gamma-rays can result from the photo-de-excitation of PeV nuclei after their parents have undergone photo-disintegration in an environment of ultraviolet photons. This process is proposed as a candidate explanation of the recently discovered HEGRA source at the edge of the Cygnus OB2 association. The Lyman-alpha background is provided by the rich O and B stellar environment. It is found that (1) the HEGRA flux can be obtained if there is efficient acceleration at the source of lower energy nuclei; (2) the requirement that the Lorentz-boosted ultraviolet photons can excite the Giant Dipole resonance implies a strong suppression of the gamma-ray spectrum compared to an E_\gamma^{-2} behavior at energies \alt 1 TeV (some of these energies will be probed by the upcoming GLAST mission); (3) a TeV neutrino counterpart from neutron decay following helium photo-disintegration will be observed at IceCube only if a major proportion of the kinetic energy budget of the Cygnus OB2 association is expended in accelerating nuclei.

It is commonly assumed that high-energy gamma-rays are made via either purely electromagnetic processes or the hadronic process of pion production, followed by decay. We investigate astrophysical contexts where a third process (A*) may dominate, namely the photo-disintegration of highly boosted nuclei followed by daughter de-excitation. Starbust regions such as Cygnus OB2 appear to be promising sites for TeV gamma-ray emission via this mechanism. A unique feature of the A* process is a sharp energy minimum ~ 10 TeV/(T/eV) for gamma-ray emission from a thermal region of temperature T. We also check that a diffuse gamma-ray component resulting from the interaction of a possible extreme-energy cosmic-ray nuclei with background radiation is well below the observed EGRET data. The A* mechanism described herein offers an important contribution to gamma-ray astronomy in the era of intense observational activity.

This working group focused mainly on the complementarity among particle physics and astrophysics. The analysis of data from both fields will better constrain theoretical models. Much of the discussion focused on detecting dark matter and susy particles, and on the potential of neutrino and gamma-ray astrophysics for seeking or constraining new physics.

Detection of ultra-high energy neutrinos will be useful for unraveling the dynamics of the most violent sources in the cosmos and for revealing the neutrino cross-section at extreme energy. Above ~ 10^20 eV, neutrinos may well be the only cosmic primariies. Thus, it is important to know the acceptances (event rate/flux) of proposed air-shower experiments for detecting ``horizontal'' neutrino events initiated in our atmosphere, and ``Earth-skimming'' events initiated in the Earth's surface rock or ocean. We calculate these acceptances for fluorescence detectors, both space-based as with the EUSO and OWL proposals, and ground-based, as with Auger, HiRes and Telescope Array. The neutrino cross-section sigma_nuN is not measured at energies above 5.2 x 10^13 eV. Although QCD extrapolations offer motivated guesses for sigma_nuN, new physics could intervene to provide a surprise. Therefore, we present the acceptances of horizontal (HAS) and upgoing (UAS) air showers as a function of sigma_nuN over the interesting range 10^(-34) to 10^(-30) cm^2. The dependences of acceptances on neutrino energy, shower-threshold energy, shower length and column density, and cloud layers are also studied. For UAS, we present acceptances for events over land (rock), and over the ocean (water). The latter are larger by about an order of magnitude, thus favoring space-based detectors. We revisit the idea of Kusenko and Weiler to infer sigma_nuN at E_nu ~ 10^20 eV from the ratio of HAS-to-UAS events, and find favorable results. Our calculation is mostly analytic. Included in the UAS calculation are realistic energy-losses for taus, and Earth-curvature effects.

We discuss a possible new resonance in active-sterile neutrino oscillations arising in theories with large extra dimensions. Fluctuations in the brane effectively increase the path-length of active neutrinos relative to the path-length of sterile neutrinos through the extra-dimensional bulk. Well below the resonance, the standard oscillation formulas apply. Well above the resonance, active-sterile oscillations are suppressed. We show that a resonance energy in the range of 30-400 MeV allows an explanation of all neutrino oscillation data, including LSND data, in a consistent four-neutrino model. A high resonance energy implies an enhanced signal in MiniBooNE. A low resonance energy implies a distorted energy spectrum in LSND, and an enhanced v_mu depletion from a stopped-pion source. The numerical value of the resonance energy may be related back to the geometric aspects of the brane world. Some astrophysical and cosmological consequences of the brane-bulk resonance are briefly sketched.

We discuss the prospects for next generation neutrino telescopes, such as IceCube, to measure the flavor ratios of high-energy astrophysical neutrinos. The expected flavor ratios at the sources are $\phi_{\nu_e}:\phi_{\nu_{\mu}}:\phi_{\nu_{\tau}} = 1:2:0$, and neutrino oscillations quickly transform these to $1:1:1$. The flavor ratios can be deduced from the relative rates of showers ($\nu_e$ charged-current, most $\nu_\tau$ charged-current, and all flavors neutral-current), muon tracks ($\nu_\mu$ charged-current only), and tau lepton lollipops and double-bangs ($\nu_\tau$ charged-current only). The peak sensitivities for these interactions are at different neutrino energies, but the flavor ratios can be reliably connected by a reasonable measurement of the spectrum shape. Measurement of the astrophysical neutrino flavor ratios tests the assumed production mechanism and also provides a very long baseline test of a number of exotic scenarios, including neutrino decay, CPT violation, and small-$\delta m^2$ oscillations to sterile neutrinos.

Resonant annihilation of ultra-high energy cosmic neutrinos (UHECnu) on the cosmic neutrino background (CnuB) into Z bosons--the Z-burst mechanism--and its associated absorption and emission phenomenology provide a unique, albeit indirect, probe of the CnuB in its present state. In this paper, we examine the implications of gravitational clustering of the CnuB in nearby galaxy clusters for the Z-burst phenomenology. In particular, we study the emission features of the Z-decay products originating from the Virgo cluster, and the potential of future cosmic ray experiments to observe clustering-enhanced Z-burst rates. We find that a detector with an exposure equivalent to three years of observations at the Extreme Universe Space Observatory (EUSO) will very likely measure these enhanced rates together with the associated UHECnu flux, provided that the latter saturates current observational limits and the neutrino masses are quasi-degenerate, m_nu_i < 0.1 eV. In the case of UHECnu fluxes below the electromagnetic cascade limit, or a hierarchical neutrino mass spectrum, an experimental sensitivity exceeding that of EUSO by at least two orders of magnitude is required to detect the clustering enhancements with any certainty.

Neutrino oscillations can be affected by decoherence induced e.g. by Planck scale suppressed interactions with the space-time foam predicted in some approaches to quantum gravity. We study the prospects for observing such effects at IceCube, using the likely flux of TeV antineutrinos from the Cygnus spiral arm. We formulate the statistical analysis for evaluating the sensitivity to quantum decoherence in the presence of the background from atmospheric neutrinos, as well as from plausible cosmic neutrino sources. We demonstrate that IceCube will improve the sensitivity to decoherence effects of ${\cal O}(E^2/M_{\rm Pl})$ by 17 orders of magnitude over present limits and, moreover, that it can probe decoherence effects of ${\cal O}(E^3/M_{\rm Pl}^2)$ which are well beyond the reach of other experiments.

Motivated by a recent change in viewing the onset of the extra-galactic component in the cosmic ray spectrum, we have fitted the observed data down to $10^{8.6}$ GeV and have obtained the corresponding power emissivity. This transition energy is well below the threshold for resonant $p\gamma$ absorption on the cosmic microwave background, and thus source evolution is an essential ingredient in the fitting procedure. Two-parameter fits in the spectral and redshift evolution indices show that a standard Fermi $E_i^{-2}$ source spectrum is excluded at larger than 95% confidence level (CL). Armed with the primordial emissivity, we follow Waxman and Bahcall to derive the associated neutrino flux on the basis of optically thin sources. For $pp$ interactions as the generating mechanism, the neutrino flux exceeds the AMANDA-B10 90%CL upper limits. In the case of $p\gamma$ dominance, the flux is marginally consistent with AMANDA-B10 data. In the new scenario the source neutrino flux dominates over the cosmogenic flux at all energies. Thus, should data from AMANDA-II prove consistent with the model, we show that IceCube can measure the characteristic power law of the neutrino spectrum, and thus provide a window on the source dynamics.

A leading candidate for the extragalactic source of high energy cosmic rays is the Fermi engine mechanism, in which protons confined by magnetic fields are accelerated to very high energy through repeated scattering by plasma shock fronts. In the process of acceleration, collisions of trapped protons with the ambient plasma produce pions which decay to electromagnetic energy and neutrinos. For optically thin sources, a strong connection between the emerging cosmic rays and secondary neutrinos can be established. In this context, we show the feasibility of using the Glashow resonance as a discriminator between the $pp$ and $p\gamma$ interactions in Fermi engines as sources of neutrinos. In particular, we demonstrate how three years of observation at the km$^3$ IceCube facility can serve as a filter for the dominance of the $pp$ interaction at the source.

High energy cosmic ray experiments have identified an excess from the region of the Galactic Plane in a limited energy range around $10^{18}$ eV (EeV). This is very suggestive of neutrons as candidate primaries, because the directional signal requires relatively-stable neutral primaries, and time-dilated neutrons can reach Earth from typical Galactic distances when the neutron energy exceeds an EeV. We here point out that if the Galactic messengers are neutrons, then those with energies below an EeV will decay in flight, providing a flux of cosmic antineutrinos above a TeV which is {\it observable} at a kilometer-scale neutrino observatory. The expected event rate per year above 1 TeV in a detector such as IceCube, for example, is 20 antineutrino showers (all flavors) and a $1^\circ$ directional signal of 4 $\bar \nu_\mu$ events. A measurement of this flux can serve to identify the first extraterrestrial point source of TeV antineutrinos.

Neutrinos from far-away sources annihilating at the Z resonance on relic neutrinos may give origin to the extreme-energy cosmic rays (EECR). If ``Z-bursts'' are responsible for the EECR events, then we show that the non-observation of cosmic ray events at energies above 2 x 10^20 eV by the AGASA Collaboration implies a lower bound around 0.3 eV on the relic neutrino mass. Since this mass exceeds the mass-squared differences inferred from oscillation physics, the bound in fact applies to all three neutrino masses. Together with the upper bound provided by comparisons of the CMB anisotropy with large-scale structure, this bound leaves only a small interval for neutrino masses around 0.3 eV, if Z-bursts are to explain the existing EECR events.

It has been proposed that propagation of cosmic rays at extreme-energy may be sensitive to Lorentz-violating metric fluctuations (``foam''). We investigate the changes in interaction thresholds for cosmic-rays and gamma-rays interacting on the CMB and IR backgrounds, for a class of stochastic models of spacetime foam. The strength of the foam is characterized by the factor (E/M_P)^a, where ``a'' is a phenomenological suppression parameter. We find that there exists a critical value of ``a'' (dependent on the particular reaction), below which the threshold energy can only be lowered, and above which the threshold energy may be raised, but at most by a factor of two. Thus, it does not appear possible in this class of models to extend cosmic-ray spectra significantly beyond their classical absorption energies. However, the lower thresholds resulting from foam may have signatures in the cosmic-ray spectrum. We also find that cosmic-ray energies cannot exceed the fundamental Planck scale, and so set a lower bound of 10^8 TeV on the context of this model for the scale of gravity. Finally, we comment in the apparent non-conservation of particle energy-momentum, and speculate on its re-emergence as dark energy in the foamy vacuum.

On the basis of existing observations (at the 4.5 \sigma level) of TeV gamma-ray outbursts from the Fanaroff-Riley I (FRI) radio galaxy Centaurus A, we estimate the accompanying neutrino flux in a scenario where both photons and neutrinos emerge from pion decay. We find a neutrino flux on Earth dF_{\nu}/dE_\nu = 4.5\times 10^{-11} (E_\nu/TeV)^{-2} TeV^{-1} cm^{-2} s^{-1}, equally spread in flavor as a result of maximal mixing. Such a flux will trigger at the IceCube facility about 10 showers/burst, with negligible background from atmospheric muons, and primary neutrino energies in excess of 100 TeV. The only other FRI radio galaxy observed in the TeV photon energy range at the 4\sigma level is M87. The burst nature of this activity is not established; however, we show that the intrinsic neutrino luminosity during the active period is the same as the Centaurus A burst. On the assumption that Centaurus A typifies the FRI population, we show that IceCube should collect 10 showers (all neutrino flavors) in 3 years, attaining a 95% CL sensitivity to the diffuse neutrino flux from FRI radio galaxies in one year of observation.

Resonant annihilation of extremely high-energy cosmic neutrinos on big-bang relic anti-neutrinos (and vice versa) into Z-bosons leads to sizable absorption dips in the neutrino flux to be observed at Earth. The high-energy edges of these dips are fixed, via the resonance energies, by the neutrino masses alone. Their depths are determined by the cosmic neutrino background density, by the cosmological parameters determining the expansion rate of the universe, and by the large redshift history of the cosmic neutrino sources. We investigate the possibility of determining the existence of the cosmic neutrino background within the next decade from a measurement of these absorption dips in the neutrino flux. As a by-product, we study the prospects to infer the absolute neutrino mass scale. We find that, with the presently planned neutrino detectors (ANITA, Auger, EUSO, OWL, RICE, and SalSA) operating in the relevant energy regime above 10^{21} eV, relic neutrino absorption spectroscopy becomes a realistic possibility. It requires, however, the existence of extremely powerful neutrino sources, which should be opaque to nucleons and high-energy photons to evade present constraints. Furthermore, the neutrino mass spectrum must be quasi-degenerate to optimize the dip, which implies m_{nu} >~ 0.1 eV for the lightest neutrino. With a second generation of neutrino detectors, these demanding requirements can be relaxed considerably.

Neutrinos may be pseudo-Dirac states, such that each generation is actually composed of two maximally-mixed Majorana neutrinos separated by a tiny mass difference. The usual active neutrino oscillation phenomenology would be unaltered if the pseudo-Dirac splittings are $\delta m^2 \alt 10^{-12}$ eV$^2$; in addition, neutrinoless double beta decay would be highly suppressed. However, it may be possible to distinguish pseudo-Dirac from Dirac neutrinos using high-energy astrophysical neutrinos. By measuring flavor ratios as a function of $L/E$, mass-squared differences down to $\delta m^2 \sim 10^{-18}$ eV$^2$ can be reached. We comment on the possibility of probing cosmological parameters with neutrinos.

We have previously shown that the decay of high-energy neutrinos from distant astrophysical sources would be revealed by flavor ratios that deviate strongly from the $\phi_{\nu_e}:\phi_{\nu_\mu}:\phi_{\nu_\tau} = 1:1:1$ expected from oscillations alone. Here we show that the deviations are significantly larger when the mixing angle $\theta_{13}$ and the CP phase $\delta$ are allowed to be nonzero. If neutrinos decay, this could allow measurement of $\theta_{13}$ and $\delta$ in IceCube and other near-term neutrino telescopes.

The recently published cosmological bound on the absolute neutrino masses obtained from the Wilkinson Microwave Anisotropy Probe (WMAP) data has important consequences for neutrino experiments and models. Taken at face value, the new bound excludes the determination of the absolute neutrino mass in the KATRIN experiment and disfavors a neutrino oscillation interpretation of the LSND experiment. Combined with the KamLAND and Super-K data, the WMAP bound defines an accessible range for the neutrinoless double beta decay amplitude. The bound also impacts the Z-burst annihilation mechanism for resonant generation of extreme-energy cosmic rays on the cosmic neutrino background in two ways: it constrains the local overdensity of neutrino dark matter which is not helpful, but it also limits the resonant energy to a favorable range. In R-parity violating SUSY models neutrino masses are generated by trilinear and bilinear lepton number violating couplings. The WMAP result improves the constraints on these couplings over their existing values by an order of magnitude.

We begin with a brief overview of highest-energy cosmic ray data, and the experiments which will perform neutrino astronomy. We then discuss two particle physics aspects of neutrinos. They are possible long-lifetime decay of the neutrino, and a measurement of the neutrino-nucleon cross-section at a CMS energy orders of magnitude beyond what can be achieved with terrestrial accelerators. Measurement of an anomalously large neutrino cross-section would indicate new physics (e.g. low string-scale, extra dimensions, precocious unification), while a smaller than expected cross-section would reveal an aspect of QCD evolution. We then discuss aspects of neutrino-primary models for the extreme-energy (EE) cosmic ray data. Primary neutrinos in extant data are motivated by the directional clustering at EE reported by the AGASA experiment. We discuss the impact of the strongly-interacting neutrino hypothesis on lower-energy physics via dispersion relations, the statistical significance of AGASA directional clustering, and the possible relevance of the Z-burst mechanism for existing EE cosmic ray data.

Existing limits on the non-radiative decay of one neutrino to another plus a massless particle (e.g., a singlet Majoron) are very weak. The best limits on the lifetime to mass ratio come from solar neutrino observations, and are $\tau/m \agt 10^{-4}$ s/eV for the relevant mass eigenstate(s). For lifetimes even several orders of magnitude longer, high-energy neutrinos from distant astrophysical sources would decay. This would strongly alter the flavor ratios from the $\phi_{\nu_e}:\phi_{\nu_{\mu}}:\phi_{\nu_{\tau}} = 1:1:1$ expected from oscillations alone, and should be readily visible in the near future in detectors such as IceCube.

Any early universe phase transition occurring after inflation has the potential to populate the universe with relic magnetic monopoles. Observations of galactic magnetic fields, as well as observations matched with models for extragalactic magnetic fields, lead to the conclusion that monopoles of mass $\lsim 10^{15}$ GeV are accelerated in these fields to relativistic velocities. We explore the possible signatures of a cosmic flux of relativistic monopoles impinging on the earth. The electromagnetically-induced signatures of monopoles are reliable. The hadronically-induced signatures are highly model-dependent. Among our findings are (i) the electromagnetic energy losses of monopoles continuously initiate a protracted shower of small intensity; (ii) monopoles may traverse the earth's diameter, making them a probe of the earth's interior structure; (iii) in addition to the direct monopole Cherenkov signal presently employed, a very attractive search strategy for monopoles is detection of their radio-Cherenkov signal produced by the coherent charge-excess in the $e^+-e^-$ shower - in fact, Cherenkov-detectors have the potential to discover a monopole flux (or limit it) several orders of magnitude below the theoretical Parker limit of $10^{-15}/\rm{cm}^2$/s/sr; (iv) it is conceivable (but not compelling) that bound states of colored monopoles may be the primary particles initiating the air showers observed above the GZK cutoff.

We show that future detectors of ultrahigh-energy cosmic-ray neutrinos will be able to measure neutrino-nucleon cross section at energies as high as 10^{11}GeV or higher. We find that the flux of up-going charged leptons per unit surface area produced by neutrino interactions below the surface is inversely proportional to the cross section. This contrasts with the rate of horizontal air showers (HAS) due to neutrino interactions in the atmosphere, which is proportional to the cross section. Thus, by comparing the HAS and up-going air shower (UAS) rates, the neutrino-nucleon cross section can be inferred. Taken together, up-going and horizontal rates ensure a healthy total event rate, regardless of the value of the cross section.

If, as recently proposed by Farrar and Piran, Cen A is the source of cosmic rays detected above the Greisen-Zatsepin-Kuz'min cutoff, neutrons are $\approx 140$ more probable than protons to be observed along its line of sight. This is because the proton flux is rendered nearly isotropic by ${\cal O}(\mu$G) intergalactic magnetic fields. With the anticipated aperture of the Southern Auger Observatory, one may expect on the order of 2 neutron events/year above $10^{20}$ eV in the line of sight of Cen A.

The observation of twenty cosmic-ray air-showers at and above 10^{20} eV poses fascinating problems for particle astrophysics: how the primary particles are accelerated to these energies, how the primaries get here through the 2.7K microwave background filling the Universe, and how the highest-energy events exhibit clustering on few-degree angular scales on the sky when charged particles are expected be bent by cosmic magnetic fields. An overview of the puzzles is presented, followed by a brief discussion of many of the models proposed to solve these puzzles. Emphasis is placed on (i) the signatures by which cosmic ray experiments in the near future will discriminate among the many proposed models, and (ii) the role neutrino primaries may play in resolving the observational issues. It is an exciting prospect that highest-energy cosmic rays may have already presented us with new physics not accessible in terrestrial accelerator searches.

This talk explores the possibility that the Universe may be populated with relic magnetic monopoles. Observations of galactic and extragalactic magnetic fields, lead to the conclusion that monopoles of mass < 10^{14} GeV are accelerated in these fields to relativistic velocities. The relativistic monopole signatures and features we derive are (i) the protracted shower development, (ii) the Cherenkov signals, (iii) the tomography of the Earth with monopoles, and (iv) a model for monopole airshowers above the GZK cutoff.

We discuss four approaches to the determination of absolute neutrino mass. These are the measurement of the zero-neutrino double beta decay rate, of the tritium decay end-point spectrum, of the cosmic ray spectrum above the GZK cutoff (in the Z-burst model), and the cosmological measurement of the power spectrum governing the CMB and large scale structure. The first two approaches are sensitive to the mass eigenstates coupling to the electron neutrino, whereas the latter two are sensitive to the heavy component of the cosmic neutrino background. All mass eigenstates are related by the $\Delta m^2$'s inferred from neutrino oscillation data. Consequently, the potential for absolute mass determination of each of the four approaches is correlated with the other three, in ways that we point out.

Directional clustering can be expected in cosmic ray observations due to purely statistical fluctuations for sources distributed randomly in the sky. We develop an analytic approach to estimate the probability of random cluster configurations, and use these results to study the strong potential of the HiRes, Auger, Telescope Array and EUSO/OWL/AirWatch facilities for deciding whether any observed clustering is most likely due to non-random sources.

The observation of cosmic--ray events above the Greisen-Kuzmin-Zatsepin (GZK) cutoff of $\sim 5\times 10^{19}$ eV challenges orthodox modeling. We discuss a possible solution which uses standard hot Big Bang cosmology and Standard Model (SM) particle physics augmented only by $\lsim $ eV neutrino masses as suggested by solar, atmospheric, and terrestrial neutrino detection. In this scheme, cosmic ray neutrinos from distant, highest-energy sources annihilate resonantly on the relic-neutrino background to produce Z-bosons. The highly-boosted ($\gamma_Z\sim 10^{11}$) Z's instantly decay, producing "Z-bursts" of highly-collimated jets of hadrons and photons. The burst content includes, on average, twenty photons and two nucleons with super-GZK energies. We show that the probability for each neutrino with energy within a fraction $\Gamma_Z/M_Z$ of the resonant value $E_R=4(eV/m_\nu)\times 10^{21}$ eV to annihilate within the halo of our galactic supercluster is likely within an order of magnitude of 1%. Depending on the magnitude of the cosmic neutrino flux above 10^{20} eV, this "local" rate for primary production may be high enough to produce the cosmic ray events observed above the GZK cutoff. Several tests of this Z-burst hypothesis for generating super-GZK events are presented, including (i) a new cutoff energy at E_R; (ii) a large $\gamma/p$ ratio for primaries near the upper end of the spectrum; (iii) directional pairing of events, and pointing to their cosmic sources; and (iv) a neutrino flux above the GZK cutoff energy which is possibly measurable directly in proposed $\gsim 10^{11}$ ton cosmic-ray detectors.

We reconsider neutrino decay as an explanation for atmospheric neutrino observations. We show that if the mass-difference relevant to the two mixed states \nu_\mu and \nu_\tau is very small (< 10^{-4} eV^2), then a very good fit to the observations can be obtained with decay of a component of \nu_\mu to a sterile neutrino and a Majoron. We discuss how the K2K and MINOS long-baseline experiments can distinguish the decay and oscillation scenarios.

If neutrinos are a significant contributor to the matter density of the universe, then they should have $\sim$ eV mass and cluster in galactic (super) cluster halos, and possibly in galactic halos as well. It was noted in the early 1980's that cosmic ray neutrinos with energy within $\delta E/E_R=\Gamma_Z/M_Z \sim 3%$ of the peak energy $E_R=4 (eV/m_{\nu})\times 10^{21}$ eV will annihilate on the nonrelativistic relic antineutrinos (and vice versa) to produce the Z-boson with an enhanced, resonant cross section of ${\cal O}(G_F)\sim 10^{-32}{cm}^2$. The result of the resonant neutrino annihilation is a hadronic Z-burst 70% of the time, which contains, on average, thirty photons and 2.7 nucleons with energies near or above the GZK cutoff energy of $5\times 10^{19}$ eV. These photons and nucleons produced within our Supergalactic halo may easily propagate to earth and initiate super-GZK air showers. Here we show that the probability for each neutrino flavor at its resonant energy to annihilate within the halo of our Supergalactic cluster is likely within an order of magnitude of 1%, with the exact value depending on unknown aspects of neutrino mixing and relic neutrino clustering. The absolute lower bound in a hot Big Bang universe for the probability to annihilate within a 50 Mpc radius (roughly a nucleon propagation distance) of earth is 0.036%. From fragmentation data for Z-decay, we estimate that the nucleons are more energetic than the photons by a factor $\sim 10$. Several tests of the hypothesis are indicated.

We show that the observed zenith angle dependence of the atmospheric neutrinos can be accounted for by neutrino decay. Furthermore, it is possible to account for all neutrino anomalies with just three flavors.

By combining experimental constraints from atmospheric and solar neutrino oscillations and the tritium beta decay endpoint, we infer upper limits of 4.4 eV on the $\nu_\mu$ and $\nu_\tau$ masses, if the universe consists of three neutrinos. For hierarchical mass spectra $m_3 \gg m_1, m_2$ or $m_3 \simeq m_2 \gg m_1$ we infer that $m_{\nu_\alpha} \lsim 0.08$ eV for $\alpha = e, \mu, \tau$. In a four neutrino universe, and assuming neutrino oscillations also account for the LSND experimental results, $m_{\nu_\alpha} \lsim 5.4$ eV. We also obtain lower limits on the masses.

We show that if the solar and atmospheric data are both described by maximal vacuum oscillations at the relevant mass scales then there exists a unique mixing matrix for three neutrino flavors. The solution necessarily conserves CP and automatically implies that there is no disappearance of atmospheric $\nu_e$, consistent with indications from the Super-Kamiokande experiment. We also investigate the consequences for three-neutrino mixing if the solar and atmospheric oscillations exhibit mixing that is large but not maximal. For non-maximal mixing $\nu_e\leftrightarrow\nu_\tau$ and $\nu_e\leftrightarrow\nu_\mu$ oscillations are predicted that may be observable in future long-baseline experiments.

We make a model-independent analysis of all available data that indicate neutrino oscillations. Using probability diagrams, we confirm that a mass spectrum with two nearly degenerate pairs of neutrinos separated by a mass gap of $\simeq1$ eV is preferred over a spectrum with one mass eigenstate separated from the others. We derive some new relations among the four-neutrino mixing matrix elements. We design four-neutrino mass matrices with three active neutrinos and one sterile neutrino that naturally incorporate maximal oscillations of atmospheric $\nu_\mu$ and explain the solar neutrino and LSND results. The models allow either a large or small angle MSW or vacuum oscillation description of the solar neutrino deficit. The models predict (i) oscillations of either $\nu_e \to \nu_\tau$ or $\nu_e \to\nu_s$ in long-baseline experiments at $L/E \gg 1$ km/GeV, with amplitude determined by the LSND oscillation amplitude and argument given by the atmospheric $\delta m^2$, and (ii) the equality of the $\nu_e$ disappearance probability, the $\nu_\mu$ disappearance probability, and the LSND $\nu_\mu\to\nu_e$ appearance probability in short-baseline experiments.

We present a four-neutrino model with three active neutrinos and one sterile neutrino which naturally has maximal $\nu_\mu\to\nu_\tau$ oscillations of atmospheric neutrinos and can explain the solar neutrino and LSND results. The model predicts $\nu_e \to \nu_\tau$ and $\nu_e \to \nu_\mu$ oscillations in long-baseline experiments with $L/E \gg 1$ km/GeV with amplitudes that are determined by the LSND oscillation amplitude and argument controlled by the atmospheric $\delta m^2$.

We suggest that the highest energy $ \gsim 10^{20}$ eV cosmic ray primaries may be relativistic magnetic monopoles. Motivations for this hypothesis are twofold: (i) conventional primaries are problematic, while monopoles are naturally accelerated to $E \sim 10^{20}$ eV by galactic magnetic fields; (ii) the observed highest energy cosmic ray flux is just below the Parker limit for monopoles. By matching the cosmic monopole production mechanism to the observed highest energy cosmic ray flux we estimate the monopole mass to be $\lsim 10^{10}$ GeV.

We review models of cosmological gamma-ray bursts (GRBs). The statistical and $\gamma$-ray transparency issues are summarized. Neutron-star and black-hole merger scenarios are described and estimates of merger rates are summarized. We review the simple fireball models for GRBs and the recent work on non-simple fireballs. Alternative cosmological models, including models where GRBs are analogs of active galactic nuclei and where they are produced by high-field, short period pulsars, are also mentioned. The value of neutrino astronomy to solve the GRB puzzle is briefly reviewed.

We suggest that the highest energy $ \gsim 10^{20} eV$ cosmic ray primaries may be relativistic magnetic monopoles. Motivations for this hypothesis are that conventional primaries are problematic, while monopoles are naturally accelerated to $E \sim 10^{20} eV$ by galactic magnetic fields. By matching the cosmic monopole production mechanism to the observed highest energy cosmic ray flux we estimate the monopole mass to be $\sim 10^{10\pm1} GeV$.

Gamma ray burst (GRB) objects are now widely thought to be at cosmological distances, and thus represent enormous energy emission. Gamma ray spectra extending to $GeV$ energies suggest the possiblity of accompanying neutrino emission, and there are several models proposed suggesting the potential detectability of such coincident neutrino bursts. With this in view, we examine possible measurements that might be conducted to give experimental data useful for astronomy, for cosmology and also neutrino properties. Of interest to astronomy and cosmology, we show how measurement of neutrino flavor ratios yields information on the nature and relative distance of the source. We point out that cosmological time dilation might be measured for these sources using neutrinos, as has been done for photons, and that neutrino oscillation lengths in the range of $1$ to $10^5~Mpc$ can be probed with GRB neutrinos. We thus note that these sources may make possible the first non-electromagnetic measurements of the scale size of the universe. We discuss tests of the weak equivalence principle, tests for flavor dependent gravitational couplings, and tests for long time scale variation of physical constants. We also show that a number of new bounds on neutrino properties (charge, mass, speed, lifetime) could be facilitated to levels well beyond those already inferred from the neutrino observation of SN1987A. We also examine the implications of these physics opportunities for designers of neutrino telescopes. We conclude that detection may be possible in planned instruments if the spectra are power law extending to the $TeV$ energy region, and if the neutrino fluxes are equal to or greater than the gamma ray fluxes. We emphasize the importance of low energy detection in future experiments

We present analytic solutions to the spatially homogeneous axion field equation, using a model potential which strongly resembles the $(1-\cos \theta)$ standard anharmonic potential. We find that the anharmonicity amplifies early universe axion density fluctuations, but only significantly so for relatively large initial misalignment angles. Enhancements of $\sim$ (2,3,4,13) result for $\theta_{\rm in}\sim (0.85,0.90,0.95,0.99)\times\pi$. The large $\theta$ result agrees with Lyth, and so validates his approximations.

We present analytic solutions to the spatially homogeneous axion field equation, using a model potential which strongly resembles the standard anharmonic $(1-\cos N\theta)$ potential, but contains only a piece-wise second order term. Our exactly soluble model for $\theta(t)$ spans the entire range $[-\pi/N,\pi/N]$. In particular, we are able to confirm (i) Turner's numeric correction factors \cite{Turner} to the adiabatic and harmonic analytic treatments of homogeneous axion oscillations, and (ii) Lyth's estimate \cite{Lyth} valid near the metastable misalignment angle $\pi/N$ at the peak of the potential. We compute the enhancement of axion density fluctuations that occurs when the axion mass becomes significant at $T\sim 1$ GeV. We find that the anharmonicity amplifies density \mbox{f}luctuations, but only significantly for relatively large initial misalignment angles. The enhancement factor is $\sim$ (2,3,4,13) for $\theta_{\rm in}\sim (0.85,0.90,0.95,0.99)\times\pi$.