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.

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.

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.

Motivated by recent IceCube observations we re-examine the idea that microquasars are high energy neutrino emitters. By stretching to the maximum the parameters of the Fermi engine we show that the nearby high-mass X-ray binary LS 5039 could accelerate protons up to above about 20 PeV. These highly relativistic protons could subsequently interact with the plasma producing neutrinos up to the maximum observed energies. After that we adopt the spatial density distribution of high-mass X-ray binaries obtained from the deep INTEGRAL Galactic plane survey and we assume LS 5039 typifies the microquasar population to demonstrate that these powerful compact sources could provide a dominant contribution to the diffuse neutrino flux recently observed by IceCube.

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 show that S-dual inflationary potentials solve the unlikeliness problem manifested in Planck data and explain the excess B-mode power observed by the BICEP2 experiment as arising from primordial tensor fluctuations.

The tensor-to-scalar ratio (r = 0.20^{+0.07}_{-0.05}) inferred from the excess B-mode power observed by the Background Imaging of Cosmic Extragalactic Polarization (BICEP2) experiment is almost twice as large as the 95% CL upper limits derived from temperature measurements of the WMAP (r<0.13) and Planck (r<0.11) space missions. Very recently, it was suggested that additional relativistic degrees of freedom beyond the three active neutrinos and photons can help to relieve this tension: the data favor an effective number of light neutrino species N_{eff} = 3.86 \pm 0.25. Since the BICEP2 ratio implies the energy scale of inflation (V_*^{1/4} \sim 2 \times 10^{16} GeV) is comparable to the grand unification scale, in this paper we investigate whether we can accommodate the required N_{eff} with three right-handed (partners of the left-handed standard model) neutrinos living in the fundamental representation of a grand unified exceptional E_6 group. We show that the superweak interactions of these Dirac states (through their coupling to a TeV-scale Z' gauge boson) lead to decoupling of right-handed neutrino just above the QCD cross over transition: 175 MeV < T_{\nu_R}^{dec} < 250 MeV. For decoupling in this transition region, the contribution of the three right-handed neutrinos to N_{eff} is suppressed by heating of the left-handed neutrinos (and photons). Consistency (within 1\sigma) with the favored N_{eff} is achieved for 4.5 TeV < M_{Z'} < 7.5 TeV. The model is fully predictive and can be confronted with future data from LHC14.

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.

Llight-element abundances probing big bang nucleosynthesis and precision data from cosmology probing the cosmic microwave background decoupling epoch have hinted at the presence of extra relativistic degrees of freedom. This is widely referred to as "dark radiation", suggesting the need for new light states in the UV completion of the standard model. We provide a brief and concise overview of the current observational status of such dark radiation and investigate the interplay between two possible interpretations of the extra light states: the right-handed partners of three Dirac neutrinos (which interact with all fermions through the exchange of a new heavy vector meson) and dark matter particles that were produced through a non-thermal mechanism, such us late time decays of massive relics. This model ties together cosmological indications of extra light states and the production of the heavy vector particle at the LHC.

Precision data from cosmology (probing the CMB decoupling epoch) and light-element abundances (probing the BBN epoch) have hinted at the presence of extra relativistic degrees of freedom, the so-called "dark radiation." We present a model independent study to account for the dark radiation by means of the right-handed partners of the three, left-handed, standard model neutrinos. We show that milli-weak interactions of these Dirac states (through their coupling to a TeV-scale Z' gauge boson) may allow the \nu_R's to decouple much earlier, at a higher temperature, than their left-handed counterparts. If the \nu_R's decouple during the quark-hadron crossover transition, they are considerably cooler than the \nu_L's and contribute less than 3 extra "equivalent neutrinos" to the early Universe energy density. For decoupling in this transition region, the 3 \nu_R generate \Delta N_\nu = 3(T_{\nu_R}/T_{\nu_ L})^4 < 3, extra relativistic degrees of freedom at BBN and at the CMB epochs. Consistency with present constraints on dark radiation permits us to identify the allowed region in the parameter space of Z' masses and couplings. Remarkably, the allowed region is within the range of discovery of LHC14.

We discuss the phenomenology and cosmology of a Standard-like Model inspired by string theory, in which the gauge fields are localized on D-branes wrapping certain compact cycles on an underlying geometry, whose intersection can give rise to chiral fermions. The energy scale associated with string physics is assumed to be near the Planck mass. To develop our program in the simplest way, we work within the construct of a minimal model with gauge-extended sector U (3)_B \times Sp (1)_L \times U (1)_{I_R} \times U (1)_L. The resulting U (1) content gauges the baryon number B, the lepton number L, and a third additional abelian charge I_R which acts as the third isospin component of an SU(2)_R. All mixing angles and gauge couplings are fixed by rotation of the U(1) gauge fields to a basis diagonal in hypercharge Y and in an anomaly free linear combination of I_R and B-L. The anomalous $Z'$ gauge boson obtains a string scale St\"uckelberg mass via a 4D version of the Green-Schwarz mechanism. To keep the realization of the Higgs mechanism minimal, we add an extra SU(2) singlet complex scalar, which acquires a VEV and gives a TeV-scale mass to the non-anomalous gauge boson Z". The model is fully predictive and can be confronted with dijet and dilepton data from LHC8 and, eventually, LHC14. We show that M_{Z"} \approx 3 - 4 TeV saturates current limits from the CMS and ATLAS collaborations. We also show that for M_{Z"} \alt 5 TeV, LHC14 will reach discovery sensitivity \agt 5\sigma. After that, we demostrate in all generality that Z" milli-weak interactions could play an important role in observational cosmology. Finally, we examine some phenomenological aspects of the supersymmetric extension of the D-brane construct.

The gauge-extended U(1)_C \times SU(2)_L \times U(1)_{I_R} \times U(1)_L model has the attractive property of elevating the two major global symmetries of the standard model (baryon number B and lepton number L) to local gauge symmetries. The U(1)_L symmetry prevents the generation of Majorana masses, leading to three superweakly interacting right-handed neutrinos. This also renders a B-L symmetry non-anomalous. We show that the superweak interactions of these Dirac states (through their coupling to the TeV-scale B-L gauge boson) permit right-handed neutrino decoupling just above the QCD phase transition: 175 MeV < T_{\nu_R}^{dec} < 250 MeV. In this transitional region, the residual temperature ratio between \nu_L and \nu_R generates extra relativistic degrees of freedom at BBN and at the CMB epochs. Consistency (within 1\sigma) with both WMAP 7-year data and the most recent estimate of the primordial ^4He mass fraction is achieved for 3 TeV < M_{B-L} < 6 TeV. The model is fully predictive, and can be confronted with dijet and dilepton data (or lack thereof) from LHC7 and, eventually, LHC14.

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.

We consider extensions of the standard model based on open strings ending on D-branes, in which gauge bosons and their associated gauginos exist as strings attached to stacks of D-branes, and chiral matter exists as strings stretching between intersecting D-branes. Under the assumptions that the fundamental string scale is in the TeV range and the theory is weakly coupled, we study models of supersymmetry for which signals of annihilating neutralino dark matter are observable. In particular, we construct a model with a supersymmetric R-symmetry violating (but R-parity conserving) effective Lagrangian that allows for the s-wave annihilation of neutralinos, once gauginos acquire mass through an unspecified mechanism. The model yields bino-like neutralinos (with the measured relic abundance) that annihilate to a gamma-gamma final state with a substantial branching fraction (~ 10%) that is orders of magnitude larger than in the minimal supersymmetric standard model. A very bright gamma-ray spectral line could be observed by gamma-ray telescopes.

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.

With the advent of the LHC there is widespread interest in the discovery potential for physics beyond the standard model. In TeV-scale open string theory, the new physics can be manifest in the excitation and decay of new resonant structures, corresponding to Regge recurrences of standard model particles. An essential input for the prediction of invariant mass spectra of the decay products (which could serve to identify the resonance as a string excitation) are the partial and total widths of the decay products. We present a parameter-free calculation of these widths for the first Regge recurrence of the SU(3) gluon octet, of the U(1) gauge boson which accompanies gluons in D-brane constructions, and of the quark triplet.

Within the context of mass-varying neutrinos, we construct a cosmological model that has a phase transition in the electromagnetic fine structure constant \alpha at a redshift of 0.5. The model accommodates hints of a time variable \alpha in quasar spectra and the nonobservance of such an effect at very low redshifts. It is consistent with limits from the recombination and primordial nucleosynthesis eras and is free of instabilities.

We explore the cosmological content of Salam-Sezgin six dimensional supergravity, and find a solution to the field equations in qualitative agreement with observation of distant supernovae, primordial nucleosynthesis abundances, and recent measurements of the cosmic microwave background. The carrier of the acceleration in the present de Sitter epoch is a quintessence field slowly rolling down its exponential potential. Intrinsic to this model is a second modulus which is automatically stabilized and acts as a source of cold dark matter with a mass proportional to an exponential function of the quintessence field (hence realizing VAMP models within a String context). However, any attempt to saturate the present cold dark matter component in this manner leads to unacceptable deviations from cosmological data -- a numerical study reveals that this source can account for up to about 7% of the total cold dark matter budget. We also show that (1) the model will support a de Sitter energy in agreement with observation at the expense of a miniscule breaking of supersymmetry in the compact space; (2) variations in the fine structure constant are controlled by the stabilized modulus and are negligible; (3) ``fifth''forces are carried by the stabilized modulus and are short range; (4) the long time behavior of the model in four dimensions is that of a Robertson-Walker universe with a constant expansion rate (w = -1/3). Finally, we present a String theory background by lifting our six dimensional cosmological solution to ten dimensions.

The Pierre Auger Observatory (PAO) has measured the spectrum and composition of the ultrahigh energy cosmic rays with unprecedented precision. We use these measurements to constrain their spectrum and composition as injected from their sources and, in turn, use these results to estimate the spectrum of cosmogenic neutrinos generated in their propagation through intergalactic space. We find that the PAO measurements can be well fit if the injected cosmic rays consist entirely of nuclei with masses in the intermediate (C, N, O) to heavy (Fe, Si) range. A mixture of protons and heavier species is also acceptable but (on the basis of existing hadronic interaction models) injection of pure light nuclei (p, He) results in unacceptable fits to the new elongation rate data. The expected spectrum of cosmogenic neutrinos can vary considerably, depending on the precise spectrum and chemical composition injected from the cosmic ray sources. In the models where heavy nuclei dominate the cosmic ray spectrum and few dissociated protons exceed GZK energies, the cosmogenic neutrino flux can be suppressed by up to two orders of magnitude relative to the all-proton prediction, making its detection beyond the reach of current and planned neutrino telescopes. Other models consistent with the data, however, are proton-dominated with only a small (1-10%) admixture of heavy nuclei and predict an associated cosmogenic flux within the reach of upcoming experiments. Thus a detection or non-detection of cosmogenic neutrinos can assist in discriminating between these possibilities.

Models of late-time neutrino mass generation contain new interactions of the cosmic background neutrinos with supernova relic neutrinos (SRNs) through exchange of the on-shell light boson, leading to significant modification of the differential SRN flux observed at earth. We consider Abelian U(1) model for generating neutrino masses at low scales and we show that there is a large parameter space in this model for which the changes induced in the flux by the exchange of the light bosons might allow one to distinguish between neutrinos being Majorana or Dirac particles, the type of neutrino mass hierarchy (normal or inverted or quasi-degenerate), and could also possibly determine the absolute values of the neutrino masses. Measurements of the presence of these effects would be possible at the next-generation water Cerenkov detectors enriched with Gadolinium, or a large 100 kton liquid argon detector.

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.

We emphasize the inelasticity distribution of events detected at the IceCube neutrino telescope as an important tool for revealing new physics. This is possible because the unique energy resolution at this facility allows to separately assign the energy fractions for emergent muons and taus in neutrino interactions. As a particular example, we explore the possibility of probing second and third generation leptoquark parameter space (coupling and mass). We show that production of leptoquarks with masses \agt 250 GeV and diagonal generation couplings of O(1) can be directly tested if the cosmic neutrino flux is at the Waxman-Bahcall level.

In models in which neutrinos are light, due to a low scale of symmetry breaking, additional light bosons are generically present. We show that the interaction between diffuse relic supernova neutrinos (RSN) and the cosmic background neutrinos, via exchange of these light scalars, can result in a dramatic change of the supernova (SN) neutrinos flux. Measurement of this effect with current or future experiments can provide a spectacular direct evidence for the low scale models. We demonstrate how the observation of neutrinos from SN1987A constrains the symmetry breaking scale of the above models. We also discuss how current and future experiments may confirm or further constrain the above models, either by detecting the ``accumulative resonance'' that diffuse RSN go through or via a large suppression of the flux of neutrinos from nearby < O(Mpc) SN bursts.

Ultra-high energy cosmic rays and neutrinos probe energies far above the weak scale. Their usefulness might appear to be limited by astrophysical uncertainties; however, by simultaneously considering up- and down-going events, one may disentangle particle physics from astrophysics. We show that present data from the AMANDA experiment in the South Pole ice already imply an upper bound on neutrino cross sections at energy scales that will likely never be probed at man-made accelerators. The existing data also place an upper limit on the neutrino flux valid for any neutrino cross section. In the future, similar analyses of IceCube data will constrain neutrino properties and fluxes at the O(10%) level.

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.

A striking aspect of the recently proposed split supersymmetry is the existence of heavy gluinos which are metastable because of the very heavy squarks which mediate their decay. In this paper we correlate the expected flux of these particles with the accompanying neutrino flux produced in inelastic $pp$ collisions in distant astrophysical sources. We show that an event rate at the Pierre Auger Observatory of approximately 1 yr$^{-1}$ for gluino masses of about 500 GeV is consistent with existing limits on neutrino fluxes. The extremely low inelasticity of the gluino-containing hadrons in their collisions with the air molecules makes possible a distinct characterization of the showers induced in the atmosphere. Should such anomalous events be observed, we show that their cosmogenic origin, in concert with the requirement that they reach the Earth before decay, leads to a lower bound on their proper lifetime of the order of 100 years, and consequently, to a lower bound on the scale of supersymmetry breaking, $\Lambda_{\rm SUSY} >2.6 \times 10^{11}$ GeV. Obtaining such a bound is not possible in collider experiments.

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.

An analysis is carried out within mSUGRA of the estimated number of events originating from upward moving ultra-high energy neutralinos that could be detected by the Extreme Universe Space Observatory (EUSO). The analysis exploits a recently proposed technique that differentiates ultra-high energy neutralinos from ultra-high energy neutrinos using their different absorption lengths in the Earth's crust. It is shown that for a significant part of the parameter space, where the neutralino is mostly a Bino and with squark mass $\sim 1$ TeV, EUSO could see ultra-high energy neutralino events with essentially no background. In the energy range 10^9 GeV < E < 10^11 GeV, the unprecedented aperture of EUSO makes the telescope sensitive to neutralino fluxes as low as 1.1 \times 10^{-6} (E/GeV)^{-1.3} GeV^{-1} cm^{-2} yr^{-1} sr^{-1}, at the 95% CL. Such a hard spectrum is characteristic of supermassive particles' $N$-body hadronic decay. The case in which the flux of ultra-high energy neutralinos is produced via decay of metastable heavy particles with uniform distribution throughout the universe is analyzed in detail. The normalization of the ratio of the relics' density to their lifetime has been fixed so that the baryon flux produced in the supermassive particle decays contributes to about 1/3 of the events reported by the AGASA Collaboration below 10^{11} GeV, and hence the associated GeV gamma-ray flux is in complete agreement with EGRET data. For this particular case, EUSO will collect between 4 and 5 neutralino events (with 0.3 of background) in ~ 3 yr of running. NASA's planned mission, the Orbiting Wide-angle Light-collectors (OWL), is also briefly discussed in this context.

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.

Black hole production in elementary particle collisions is among the most promising probes of large extra spacetime dimensions. Studies of black holes at particle colliders have assumed that all of the incoming energy is captured in the resulting black hole. We incorporate the inelasticity inherent in such processes and determine the prospects for discovering black holes in colliders and cosmic ray experiments, employing a dynamical model of Hawking evolution. At the Large Hadron Collider, inelasticity reduces rates by factors of 10^3 to 10^6 in the accessible parameter space, moderating, but not eliminating, hopes for black hole discovery. At the Pierre Auger Observatory, rates are suppressed by a factor of 10. We evaluate the impact of cosmic ray observations on collider prospects.

In this talk we discuss aspects of TeV-scale gravitational collapse to black holes and string balls and their subsequent evaporation. Special emphasis is placed on the interplay of the string/black hole correspondence principle. These ideas are then explored in the context of cosmic ray physics. First, the potential for observing showers mediated by black holes or superstring excitations is examined. Next, existing data from neutrino telescopes are used to constrain the parameter space for the unseen dimensions of the universe. Finally, we close with a discussion of future prospects.

High energy cosmic ray data collected by the SUGAR Collaboration in the Pilliga State Forest have been re-analyzed using up-to-date shower simulations. Complete, inclined angle, Monte Carlo simulations reveal 2 events with energies in excess of $7\times 10^{19}$ eV at 95% CL, independent of the choice of hadronic interaction model and of chemical composition of the primary. An additional 3 events, with mean energy $\agt 7\times 10^{19}$ eV, were also re-analyzed in the same manner. A lower bound on the flux at the high end of the spectrum, as observed in the Southern sky, is presented on the basis of our analysis.

There are indications from the study of quasar absorption spectra that the fine structure constant $\alpha$ may have been measurably smaller for redshifts $z>2.$ Analyses of other data ($^{149}$Sm fission rate for the Oklo natural reactor, variation of $^{187}$Re $\beta$-decay rate in meteorite studies, atomic clock measurements) which probe variations of $\alpha$ in the more recent past imply much smaller deviations from its present value. In this work we tie the variation of $\alpha$ to the evolution of the quintessence field proposed by Albrecht and Skordis, and show that agreement with all these data, as well as consistency with WMAP observations, can be achieved for a range of parameters. Some definite predictions follow for upcoming space missions searching for violations of the equivalence principle.

We revise existing limits on the D-dimensional Planck scale M_D from the nonobservation of microscopic black holes produced by high energy cosmic neutrinos in scenarios with D=4+n large extra dimensions. Previous studies have neglected the energy radiated in gravitational waves by the multipole moments of the incoming shock waves. We include the effects of energy loss, as well as form factors for black hole production and recent null results from cosmic ray detectors. For n>4, we obtain M_D > 1.0 - 1.4 TeV. These bounds are among the most stringent and conservative to date.

The starburst galaxies M82 and NGC253 have been proposed as the primary sources of cosmic rays with energies above $10^{18.7}$ eV. For energies $\agt 10^{20.3}$ eV the model predicts strong anisotropies. We calculate the probabilities that the latter can be due to chance occurrence. For the highest energy cosmic ray events in this energy region, we find that the observed directionality has less than 1% probability of occurring due to random fluctuations. Moreover, during the first 5 years of operation at Auger, the observation of even half the predicted anisotropy has a probability of less than $10^{-5}$ to occur by chance fluctuation. Thus, this model can be subject to test at very small cost to the Auger priors budget and, whatever the outcome of that test, valuable information on the Galactic magnetic field will be obtained.

Ultra-high energy cosmic neutrinos are incisive probes of both astrophysical sources and new TeV-scale physics. Such neutrinos would create extensive air showers deep in the atmosphere. The absence of such showers implies upper limits on incoming neutrino fluxes and cross sections. Combining the exposures of AGASA, the largest existing ground array, with the exposure of the Fly's Eye fluorescence detector integrated over all its operating epochs, we derive 95% CL bounds that substantially improve existing limits. We begin with model-independent bounds on astrophysical fluxes, assuming standard model cross sections, and model-independent bounds on new physics cross sections, assuming a conservative cosmogenic flux. We then derive model-dependent constraints on new components of neutrino flux for several assumed power spectra, and we update bounds on the fundamental Planck scale M_D in extra dimension scenarios from black hole production. For large numbers of extra dimensions, we find M_D > 2.0 (1.1) TeV for \mbhmin = M_D (5 M_D), comparable to or exceeding the most stringent constraints to date.

We explore the phenomenology of microscopic black holes in the $S^1/Z_2$ Randall-Sundrum (RS) model. We consider the canonical framework in which both gauge and matter fields are confined to the brane and only gravity spills into the extra dimension. The model is characterized by two parameters, the mass of the first massive graviton $(m_1)$, and the curvature $1/\ell$ of the RS anti-de Sitter space. We compute the sensitivity of present and future cosmic ray experiments to various regions of $\ell$ and $m_1,$ and compare with that of Runs I and II at the Tevatron. As part of our phenomenological analysis, we examine constraints placed on $\ell$ by AdS/CFT considerations.

If extra spacetime dimensions and low-scale gravity exist, black holes will be produced in observable collisions of elementary particles. For the next several years, ultra-high energy cosmic rays provide the most promising window on this phenomenon. In particular, cosmic neutrinos can produce black holes deep in the Earth's atmosphere, leading to quasi-horizontal giant air showers. We determine the sensitivity of cosmic ray detectors to black hole production and compare the results to other probes of extra dimensions. With n \ge 4 extra dimensions, current bounds on deeply penetrating showers from AGASA already provide the most stringent bound on low-scale gravity, requiring a fundamental Planck scale M_D > 1.3 - 1.8 TeV. The Auger Observatory will probe M_D as large as 4 TeV and may observe on the order of a hundred black holes in 5 years. We also consider the implications of angular momentum and possible exponentially suppressed parton cross sections; including these effects, large black hole rates are still possible. Finally, we demonstrate that even if only a few black hole events are observed, a standard model interpretation may be excluded by comparison with Earth-skimming neutrino rates.

In spacetimes with asymmetric extra dimensions, cosmic neutrino interactions may be extraordinarily enhanced by p-brane production. Brane formation and decay may then initiate showers deep in the Earth's atmosphere at rates far above the standard model rate. We explore the p-brane discovery potential of cosmic ray experiments. The absence of deeply penetrating showers at AGASA already provides multi-TeV bounds on the fundamental Planck scale that significantly exceed those obtained from black hole production in symmetric compactification scenarios. This sensitivity will be further enhanced at the Auger Observatory. We also examine the possibility that p-brane formation resolves the GZK paradox. For flat compactifications, astrophysical bounds exclude this explanation. For warped scenarios, a solution could be consistent with the absence of deep showers only for extra dimensions with fine-tuned sizes well below the fundamental Planck length. In addition, it requires moderately penetrating showers, so far not reported, and ~100% modifications to standard model phenomenology at 100 GeV energies.

Assuming that the hard gamma-ray emission of Cen A is a result of synchrotron radiation of ultra-relativistic electrons, we derive a lower bound on the local extragalactic magnetic field, $B> 10^{-8}$ G. This result is consistent with (and close to) upper bounds on magnetic fields derived from consideration of cosmic microwave background distortions and Faraday rotation measurements.

The existence of extra degrees of freedom beyond the electroweak scale may allow the formation of black holes in nearly horizontal neutrino air showers. In this paper we examine the average properties of the light descendants of these black holes. Our analysis indicates that black hole decay gives rise to deeply penetrating showers with an electromagnetic component which differs substantially from that in conventional neutrino interactions, allowing a good characterization of the phenomenon against background. Naturally occurring black holes in horizontal neutrino showers could be detected and studied with the Auger air shower array. Since the expected black hole production rate at Auger is $> 1$ event/year, this cosmic ray observatory could be potentially powerful in probing models with extra dimensions and TeV-scale gravity.

We examine the correlation between compact radio quasars (redshifts in the range $z = 0.3 - 2.2$) and the arrival direction of ultrahigh energy cosmic rays forming clusters. Our Monte Carlo simulation reveals a statistically significant correlation on the AGASA sample: the chance probability of this effect being less than 1%. The implications of this result on the origin of ultrahigh energy cosmic rays are discussed.

In this article we examine the hypothesis that the highest energy cosmic rays are complex nuclei from extragalactic sources. Under reasonable physical assumptions, we show that the nearby metally rich starburst galaxies (M82 and NGC 253) can produce all the events observed above the ankle. This requires diffusion of particles below $10^{20}$ eV in extragalactic magnetic fields $B \approx 15$ nG. Above $10^{19}$ eV, the model predicts the presence of significant fluxes of medium mass and heavy nuclei with small rate of change of composition. Notwithstanding, the most salient feature of the starburst-hypothesis is a slight anisotropy induced by iron debris just before the spectrum-cutoff.

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.

At Auger energies only model predictions enable us to extract primary cosmic ray features. The simulation of the shower evolution depends sensitively on the first few interactions, necessarily related to the quality of our understanding of high energy hadronic collisions. Distortions of the standard ``soft semi-hard'' scenario include novel large compact dimensions and a string or quantum gravity scale not far above the electroweak scale. Na\"{\i}vely, the additional degrees of freedom yield unification of all forces in the TeV range. In this article we study the influence of such precocious unification during atmospheric cascade developments by analyzing the most relevant observables in proton induced showers.

One of the possible consequences of the existence of extra degrees of freedom beyond the electroweak scale is the increase of neutrino-nucleon cross sections ($\sigma_{\nu N}$) beyond Standard Model predictions. At ultra-high energies this may allow the existence of neutrino-initiated extensive air showers. In this paper, we examine the most relevant observables of such showers. Our analysis indicates that the future Pierre Auger Observatory could be potentially powerful in probing models with large compact dimensions.

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.

It is proposed that the apparent positive acceleration of the cosmological scale factor is due to the vacuum energy of an incomplete chiral phase transition in a hidden SU(2) sector. Constraints from primordial nucleosynthesis imply that the present metastable phase is in a substantially supercooled state. It is argued that massless chiral condensates can substantially enhance the possibility of supercooling, and a linear sigma model exhibiting scale invariance broken only at the quantum level is shown to accommodate the required supercooling with a reasonable choice of quartic couplings. The extensive supercooling can in principle be confirmed or rejected on the basis of interface tension measurements in lattice simulations with dynamical fermions.

I present a model for inflation based on the gaugino-dilaton dynamics of supersymmetry breaking. The inflaton is the dimension--1 scalar field $\phi$ related to the gaugino condensate via $\lambda^T\lambda=\phi^3.$ Recent work in this area is used to obtain two significant results: (1) Scalar density fluctuations at second horizon crossing are generated on scale $\lambda$ with amplitude \[ \drr=A\cdot \ 10^3\cdot\ \left(\frac{m_{SUSY}}{m_{Pl}}\right)^{\half}\ \left(\frac{\lambda}{100\ \mbox{Mpc}}\right)^{0.03}\ \ , \] where $A$ is a constant which depends on (unknown) details of gaugino-dilaton dynamics in the strong-coupling phase. Agreement with COBE results is obtained if $A\simeq 2.$ \ (2) Due to mixing with hidden sector glueballs, the dilaton mass is large $(\sim \left(\mp\ \ms^2\right)^ {\frac{1}{3}}\sim 10^8\ \mbox{GeV}),$ and reheating takes place at $T\sim 1\tev.$ This {\em necessitates} that the presently observed baryon asymmetry be generated at a cosmic temperature below 1 TeV.