Our ability to identify the sources of cosmic rays and understand how these particles propagate through the interstellar medium is hindered by the combined effects of the solar wind and its embedded magnetic field, collectively known as solar modulation. In this paper, we build upon our previous work to model and constrain the effects of solar modulation on the cosmic-ray spectrum, using data from AMS-02 and BESS Polar II collected between 2007 and 2012, during which the heliospheric magnetic field was in a state of negative polarity. Our model uses measurements of the heliospheric magnetic field and the tilt angle of the heliospheric current sheet to accurately predict the effects of solar modulation as a function of time, charge, and rigidity. By incorporating data from a period of negative polarity, we have been able to robustly observe and constrain the charge-dependent effects of solar modulation.

It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

Measurements of the cosmic microwave background (CMB) spectral distortions (SDs) will open a new window on the very early universe, providing new information complementary to that gathered from CMB temperature and polarization anisotropies. In this paper, we study their synergy as a function of the characteristics of the considered experiments. In particular, we examine a wide range of sensitivities for possible SD measurements, spanning from FIRAS up to noise levels 1000 times better than PIXIE, and study their constraining power when combined with current or future CMB anisotropy experiments such as Planck or LiteBIRD plus CMB-S4. We consider a number of different cosmological models such as the $\Lambda$CDM, as well as its extensions with the running of the scalar spectral index, the decay or the annihilation of dark matter (DM) particles. While upcoming CMB anisotropy experiments will be able to decrease the uncertainties on inflationary parameters such as $A_s$ and $n_s$ by about a factor 2 in the $\Lambda$CDM case, we find that an SD experiment 10 times more sensitive than PIXIE (comparable to the proposed PRISM satellite) could potentially further contribute to constrain these parameters. This is even more significant in the case of the running of the scalar spectral index. Furthermore, as expected, constraints on DM particles decaying at redshifts probed by SDs will improve by orders of magnitude even with an experiment 10 times worse than PIXIE as compared to CMB anisotropies or Big Bang Nucleosynthesis bounds. On the contrary, DM annihilation constraints will not significantly improve over CMB anisotropy measurements. For some of the cases considered, we study the impact of marginalizing over the contribution from reionization and structure formation. Finally, we forecast the constraints obtainable with sensitivities achievable either from the ground or from a balloon.

As space expands, the energy density in black holes increases relative to that of radiation, providing us with motivation to consider scenarios in which the early universe contained a significant abundance of such objects. In this study, we revisit the constraints on primordial black holes derived from measurements of the light element abundances. Black holes and their Hawking evaporation products can impact the era of Big Bang Nucleosynthesis (BBN) by altering the rate of expansion at the time of neutron-proton freeze-out, as well as by radiating mesons which can convert protons into neutrons and vice versa. Such black holes can thus enhance the primordial neutron-to-proton ratio, and increase the amount of helium that is ultimately produced. Additionally, the products of Hawking evaporation can break up helium nuclei, which both reduces the helium abundance and increases the abundance of primordial deuterium. Building upon previous work, we make use of modern deuterium and helium measurements to derive stringent constraints on black holes which evaporate in $t_{\rm evap} \sim 10^{-1}$ s to $\sim 10^{13}$ s (corresponding to $M \sim 6\times 10^8$ g to $\sim 2 \times 10^{13}$ g, assuming Standard Model particle content). We also consider how physics beyond the Standard Model could impact these constraints. Due to the gravitational nature of Hawking evaporation, the rate at which a black hole evaporates, and the types of particles that are produced through this process, depend on the complete particle spectrum. Within this context, we discuss scenarios which feature a large number of decoupled degrees-of-freedom (\ie~large hidden sectors), as well as models of TeV-scale supersymmetry.

If a fraction of the dark matter is unstable and decays into dark radiation at around the time of matter-radiation equality, it could impact the expansion history of the universe in a way that helps to ameliorate the long-standing tension between the locally measured value of the Hubble constant and the value inferred from measurements of the cosmic microwave background and baryon acoustic oscillations (assuming standard $\Lambda$CDM cosmology). If this component of decaying dark matter is cold, however, it will modify the evolution of the gravitational potentials, leading to inconsistencies with these same data sets. With this in mind, we consider here a component of decaying warm dark matter, with a free-streaming length that is long enough to remain consistent with existing data. We study the background and perturbation evolution of warm decaying dark matter, and use cosmological data to constrain the mass, abundance and decay rate of such a particle. We find that a component of warm decaying dark matter can significantly reduce the tension between local and cosmological determinations of the Hubble constant.

Dark matter coupled solely gravitationally can be produced through the decay of primordial black holes in the early universe. If the dark matter is lighter than the initial black hole temperature, it could be warm enough to be subject to structure formation constraints. In this paper we perform a more precise determination of these constraints. We first evaluate the dark matter phase-space distribution, without relying on the instantaneous decay approximation. We then interface this phase-space distribution with the Boltzmann code CLASS to extract the corresponding matter power spectrum, which we find to match closely those of warm dark matter models, albeit with a different dark matter mass. This mapping allows us to extract constraints from Lyman-$\alpha$ data without the need to perform hydrodynamical simulations. We robustly rule out the possibility, consistent with previous analytic estimates, of primordial black holes having come to dominate the energy density of the universe and simultaneously given rise to all the DM through their decay. Consequences and implications for dark radiation and leptogenesis are also briefly discussed.

Any abundance of black holes that was present in the early universe will evolve as matter, making up an increasingly large fraction of the total energy density as space expands. This motivates us to consider scenarios in which the early universe included an era that was dominated by low-mass ($M < 5\times 10^8$ g) black holes which evaporate prior to primordial nucleosynthesis. In significant regions of parameter space, these black holes will become gravitationally bound within binary systems, and undergo mergers before evaporating. Such mergers result in three potentially observable signatures. First, any black holes that have undergone one or more mergers will possess substantial angular momentum, causing their Hawking evaporation to produce significant quantities of high-energy gravitons. These products of Hawking evaporation are predicted to constitute a background of hot ($\sim$eV-keV) gravitons today, with an energy density corresponding to $\Delta N_{\rm eff} \sim 0.01-0.03$. Second, these mergers will produce a stochastic background of high-frequency gravitational waves. And third, the energy density of these gravitational waves can be as large as $\Delta N_{\rm eff} \sim 0.3$, depending on the length of time between the mergers and evaporation. These signals are each potentially within the reach of future measurements.

The emergence of an increasingly strong tension between the Hubble rate inferred from early and late time observations has reinvigorated interest in non-standard scenarios, with the aim of reconciling these measurements. One such model involves interactions between Dark Matter and Dark Energy. Here we consider a specific form of the coupling between these two fluids proportional to the Dark Energy energy density, which has been shown to mitigate the Hubble tension. We complement the work already discussed in several previous analysis and show that, once all relevant cosmological probes are included, the value of the Hubble parameter in this model is $H_0=69.82_{-0.76}^{+0.63}$ km/(s Mpc). Furthermore, we also perform a statistical model comparison, finding a $\Delta\chi^2$ of $-2.15$ with the inclusion of one additional free parameter. In order to fully test this model as a solution to the $H_0$ tension, future missions such as Euclid will prove vital. For this reason, we also forecast the sensitivity of upcoming planned missions to the main parameters affected by these dark sector interactions and show that within the next decade we will be able to conclusively confirm or disprove the validity of this model as a solution to the $H_0$ tension.

The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a "smoking gun" signature of dark matter annihilation or decay, essentially free of astrophysical background. Studies in recent years have emphasized that models for cosmic-ray antideuterons must be considered together with the abundant cosmic antiprotons and any potential observation of antihelium. Therefore, a second dedicated Antideuteron Workshop was organized at UCLA in March 2019, bringing together a community of theorists and experimentalists to review the status of current observations of cosmic-ray antinuclei, the theoretical work towards understanding these signatures, and the potential of upcoming measurements to illuminate ongoing controversies. This review aims to synthesize this recent work and present implications for the upcoming decade of antinuclei observations and searches. This includes discussion of a possible dark matter signature in the AMS-02 antiproton spectrum, the most recent limits from BESS Polar-II on the cosmic antideuteron flux, and reports of candidate antihelium events by AMS-02; recent collider and cosmic-ray measurements relevant for antinuclei production models; the state of cosmic-ray transport models in light of AMS-02 and Voyager data; and the prospects for upcoming experiments, such as GAPS. This provides a roadmap for progress on cosmic antinuclei signatures of dark matter in the coming years.

Spectral distortions and anisotropies of the CMB provide independent and complementary probes to study energy injection processes in the early universe. Here we discuss the synergy between these observables, and show the promising future of spectral distortion missions to constrain both exotic and non-exotic energy injections. We show that conventional probes such as Big Bang Nucleosynthesis and CMB anisotropies can benefit from and even be surpassed by future spectral distortion experiments. For this, we have implemented a unified framework within the Boltzmann code CLASS to consistently treat the thermal evolution of photons and baryons. Furthermore, we give an extensive and pedagogical introduction into the topic of spectral distortions and energy injections throughout the thermal history of the universe, highlighting some of their unique features and potential as a novel probe for cosmology and particle physics.

Recent studies of the cosmic-ray antiproton-to-proton ratio have identified an excess of $\sim$10-20 GeV antiprotons relative to the predictions of standard astrophysical models. Intriguingly, the properties of this excess are consistent with the same range of dark matter models that can account for the long-standing excess of $\gamma$-rays observed from the Galactic Center. Such dark matter candidates can also produce significant fluxes of anti-deuterium and anti-helium nuclei. Here we study the production and transport of such particles, both from astrophysical processes as well as from dark matter annihilation. Importantly, in the case of AMS-02, we find that Alfv\'{e}nic reacceleration (i.e., diffusion in momentum space) can boost the expected number of $\bar{\rm d}$ and ${}^{3}\overline{\textrm{He}}$ events from annihilating dark matter by an order of magnitude or more. For relatively large values of the Alfv\'{e}n speed, and for dark matter candidates that are capable of producing the antiproton and $\gamma$-ray excesses, we expect annihilations to produce a few anti-deuteron events and about one anti-helium event in six years of AMS-02 data. This is particularly interesting in light of recent reports from the AMS-02 Collaboration describing the detection of a number of anti-helium candidate events.

Although weakly interacting massive particles (WIMPs) have long been among the most studied and theoretically attractive classes of candidates for the dark matter of our universe, the lack of their detection in direct detection and collider experiments has begun to dampen enthusiasm for this paradigm. In this study, we set out to appraise the status of the WIMP paradigm, focusing on the case of dark matter candidates that interact with the Standard Model through a new gauge boson. After considering a wide range of $Z'$ mediated dark matter models, we quantitatively evaluate the fraction of the parameter space that has been excluded by existing experiments, and that is projected to fall within the reach of future direct detection experiments. Despite the existence of stringent constraints, we find that a sizable fraction of this parameter space remains viable. More specifically, if the dark matter is a Majorana fermion, we find that an order one fraction of the parameter space is in many cases untested by current experiments. Future direct detection experiments with sensitivity near the irreducible neutrino floor will be able to test a significant fraction of the currently viable parameter space, providing considerable motivation for the next generation of direct detection experiments.

An excess of $\sim$10-20 GeV cosmic-ray antiprotons has been identified in the spectrum reported by the AMS-02 Collaboration. The systematic uncertainties associated with this signal, however, have made it difficult to interpret these results. In this paper, we revisit the uncertainties associated with the time, charge and energy-dependent effects of solar modulation, the antiproton production cross section, and interstellar cosmic-ray propagation. After accounting for these uncertainties, we confirm the presence of a 4.7$\sigma$ antiproton excess, consistent with that arising from a $m_{\chi} \approx 64-88$ GeV dark matter particle annihilating to $b\bar{b}$ with a cross section of $\sigma v \simeq (0.8-5.2) \times 10^{-26}$ cm$^{3}$/s. If we allow for the stochastic acceleration of secondary antiprotons in supernova remnants, the data continues to favor a similar range of dark matter models ($m_{\chi}\approx 46-94$ GeV, $\sigma v \approx (0.7-3.8)\times 10^{-26}$ cm$^3/$s) with a significance of 3.3$\sigma$. The same range of dark matter models that are favored to explain the antiproton excess can also accommodate the excess of GeV-scale gamma rays observed from the Galactic Center.

Many attempts to solve the Hubble tension with extended cosmological models combine an enhanced relic radiation density, acting at the level of background cosmology, with new physical ingredients affecting the evolution of cosmological perturbations. Several authors have pointed out the ability of combined Baryon Acoustic Oscillation (BAO) and Big Bang Nucleosynthesis (BBN) data to probe the background cosmological history independently of both CMB maps and supernovae data. Using state-of-the-art assumptions on BBN, we confirm that combined BAO, deuterium, and helium data are in tension with the SH0ES measurements under the $\Lambda$CDM assumption at the 3.2$\sigma$ level, while being in close agreement with the CMB value. We subsequently show that floating the radiation density parameter $N_\mathrm{eff}$ only reduces the tension down to the 2.6$\sigma$ level. This conclusion, totally independent of any CMB data, shows that a high $N_\mathrm{eff}$ accounting for extra relics (either free-streaming or self-interacting) does not provide an obvious solution to the crisis, not even at the level of background cosmology. To circumvent this strong bound, (i) the extra radiation has to be generated after BBN to avoid helium bounds, and (ii) additional ingredients have to be invoked at the level of perturbations to reconcile this extra radiation with CMB and LSS data.

We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (\phi) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation.

Several interesting Dark Matter (DM) models invoke a dark sector leading to two types of relic particles, possibly interacting with each other: non-relativistic DM, and relativistic Dark Radiation (DR). These models have interesting consequences for cosmological observables, and could in principle solve problems like the small-scale cold DM crisis, Hubble tension, and/or low $\sigma_8$ value. Their cosmological behaviour is captured by the ETHOS parametrisation, which includes a DR-DM scattering rate scaling like a power-law of the temperature, $T^n$. Scenarios with $n=0$, $2$, or $4$ can easily be realised in concrete dark sector set-ups. Here we update constraints on these three scenarios using recent CMB, BAO, and high-resolution Lyman-$\alpha$ data. We introduce a new Lyman-$\alpha$ likelihood that is applicable to a wide range of cosmological models with a suppression of the matter power spectrum on small scales. For $n=2$ and $4$, we find that Lyman-$\alpha$ data strengthen the CMB+BAO bounds on the DM-DR interaction rate by many orders of magnitude. However, models offering a possible solution to the missing satellite problem are still compatible with our new bounds. For $n=0$, high-resolution Lyman-$\alpha$ data bring no stronger constraints on the interaction rate than CMB+BAO data, except for extremely small values of the DR density. Using CMB+BAO data and a theory-motivated prior on the minimal density of DR, we find that the $n=0$ model can reduce the Hubble tension from $4.1\sigma$ to $2.7\sigma$, while simultaneously accommodating smaller values of the $\sigma_8$ and $S_8$ parameters hinted by cosmic shear data.

If the dark matter is part of a hidden sector with only very feeble couplings to the Standard Model, the lightest particle in the hidden sector will generically be long-lived and could come to dominate the energy density of the universe prior to the onset of nucleosynthesis. During this early matter-dominated era, density perturbations will grow more quickly than otherwise predicted, leading to a large abundance of sub-earth-mass dark matter microhalos. Since the dark matter does not couple directly to the Standard Model, the minimum halo mass is much smaller than expected for weakly interacting dark matter, and the smallest halos could form during the radiation-dominated era. In this paper, we calculate the evolution of density perturbations within the context of such hidden sector models and use a series of $N$-body simulations to determine the outcome of nonlinear collapse during radiation domination. The resulting microhalos are extremely dense, which leads to very high rates of dark matter annihilation and to large indirect detection signals that resemble those ordinarily predicted for decaying dark matter. We find that the Fermi Collaboration's measurement of the high-latitude gamma-ray background rules out a wide range of parameter space within this class of models. The scenarios that are most difficult to constrain are those that feature a very long early matter-dominated era; if microhalos form prior to the decay of the unstable hidden sector matter, the destruction of these microhalos effectively heats the dark matter, suppressing the later formation of microhalos.

We forecast the sensitivity of thirty-five different combinations of future Cosmic Microwave Background and Large Scale Structure data sets to cosmological parameters and to the total neutrino mass. We work under conservative assumptions accounting for uncertainties in the modelling of systematics. In particular, for galaxy redshift surveys, we remove the information coming from non-linear scales. We use Bayesian parameter extraction from mock likelihoods to avoid Fisher matrix uncertainties. Our grid of results allows for a direct comparison between the sensitivity of different data sets. We find that future surveys will measure the neutrino mass with high significance and will not be substantially affected by potential parameter degeneracies between neutrino masses, the density of relativistic relics, and a possible time-varying equation of state of Dark Energy.

If even a relatively small number of black holes were created in the early universe, they will constitute an increasingly large fraction of the total energy density as space expands. It is thus well-motivated to consider scenarios in which the early universe included an era in which primordial black holes dominated the total energy density. Within this context, we consider Hawking radiation as a mechanism to produce both dark radiation and dark matter. If the early universe included a black hole dominated era, we find that Hawking radiation will produce dark radiation at a level $\Delta N_{\rm eff} \sim 0.03-0.2$ for each light and decoupled species of spin 0, 1/2, or 1. This range is well suited to relax the tension between late and early-time Hubble determinations, and is within the reach of upcoming CMB experiments. The dark matter could also originate as Hawking radiation in a black hole dominated early universe, although such dark matter candidates must be very heavy ($m_{\rm DM} >10^{11}$ GeV) if they are to avoid exceeding the measured abundance.

The ANITA experiment, which is designed to detect ultra-high energy neutrinos, has reported the observation of two anomalous events, directed at angles of $27^{\circ}$ and $35^{\circ}$ with respect to the horizontal. At these angles, the Earth is expected to efficiently absorb ultra-high energy neutrinos, making the origin of these events unclear and motivating explanations involving physics beyond the Standard Model. In this study, we consider the possibility that ANITA's anomalous events are the result of Askaryan emission produced by exotic weakly interacting particles scattering elastically with nuclei in the Antarctic ice sheet. Such particles could be produced by superheavy ($\sim 10^{10}-10^{13}$ GeV) dark matter particles decaying in the halo of the Milky Way. Such scenarios can be constrained by existing measurements of the high-latitude gamma-ray background and the ultra-high energy cosmic ray spectrum, along with searches for ultra-high energy neutrinos by IceCube and other neutrino telescopes.

In this paper, we explore in detail the cosmological implications of an abelian $L_\mu-L_\tau$ gauge extension of the Standard Model featuring a light and weakly coupled $Z'$. Such a scenario is motivated by the longstanding $\sim \, 4 \sigma$ discrepancy between the measured and predicted values of the muon's anomalous magnetic moment, $(g-2)_\mu$, as well as the tension between late and early time determinations of the Hubble constant. If sufficiently light, the $Z'$ population will decay to neutrinos, increasing the overall energy density of radiation and altering the expansion history of the early universe. We identify two distinct regions of parameter space in this model in which the Hubble tension can be significantly relaxed. The first of these is the previously identified region in which a $\sim \, 10-20$ MeV $Z'$ reaches equilibrium in the early universe and then decays, heating the neutrino population and delaying the process of neutrino decoupling. For a coupling of $g_{\mu-\tau} \simeq (3-8) \times 10^{-4}$, such a particle can also explain the observed $(g-2)_{\mu}$ anomaly. In the second region, the $Z'$ is very light ($m_{Z'} \sim 1\,\text{eV}$ to $\text{MeV}$) and very weakly coupled ($g_{\mu-\tau} \sim 10^{-13}$ to $10^{-9}$). In this case, the $Z'$ population is produced through freeze-in, and decays to neutrinos after neutrino decoupling. Across large regions of parameter space, we predict a contribution to the energy density of radiation that can appreciably relax the reported Hubble tension, $\Delta N_{\rm eff} \simeq 0.2$.

In these lectures, I describe a variety of efforts to identify or constrain the identity of dark matter by detecting the annihilation or decay products of these particles, or their effects. After reviewing the motivation for indirect searches, I discuss what we have learned about dark matter from observations of gamma rays, cosmic rays and neutrinos, as well as the cosmic microwave background. Measurements such as these have been used to significantly constrain a wide range of thermal relic dark matter candidates, in particular those with masses below a few hundred GeV. I also discuss a number of anomalies and excesses that have been interpreted as possible signals of dark matter, including the Galactic Center gamma-ray excess, the cosmic-ray antiproton excess, the cosmic-ray positron excess, and the 3.5 keV line. These lectures were originally presented as part of the 2018 Theoretical Advanced Study Institute (TASI) summer school on "Theory in an Era of Data". Although intended for advanced graduate students, these lectures may be useful for a wide range of physicists, astrophysicists and astronomers who wish to get an overview of the current state of indirect searches for dark matter.

Multi-messenger data suggest that radio galaxies (i.e. non-blazar active galaxies) are perhaps the most likely class of sources for the diffuse flux of high-energy neutrinos reported by the IceCube Collaboration. In this study, we consider the gamma-ray spectrum observed from four nearby radio galaxies (Centaurus A, PKS 0625-35, NGC 1275 and IC 310) and constrain the intensity and spectral shape of the emission injected from these sources, accounting for the effects of attenuation and contributions from electromagnetic cascades (initiated both within the radio galaxy itself and during extragalactic propagation). Assuming that this gamma-ray emission is generated primarily through the interactions of cosmic-ray protons with gas, we calculate the neutrino flux predicted from each of these sources. Although this scenario is consistent with the constraints published by the IceCube and ANTARES Collaborations, the predicted fluxes consistently fall within an order of magnitude of the current point source sensitivity. The prospects appear very encouraging for the future detection of neutrino emission from the nearest radio galaxies.

If the dark matter is unstable, the decay of these particles throughout the universe and in the halo of the Milky Way could contribute significantly to the isotropic gamma-ray background (IGRB) as measured by Fermi. In this article, we calculate the high-latitude gamma-ray flux resulting from dark matter decay for a wide range of channels and masses, including all contributions from inverse Compton scattering and accounting for the production and full evolution of cosmological electromagnetic cascades. We also make use of recent multi-wavelength analyses that constrain the astrophysical contributions to the IGRB, enabling us to more strongly restrict the presence any component arising from decaying dark matter. Over a wide range of decay channels and masses (from GeV to EeV and above), we derive stringent lower limits on the dark matter's lifetime, generally in the range of $\tau \sim (1-5)\times 10^{28}$ s.

The excess of neutrino candidate events detected by IceCube from the direction of TXS 0506+056 has generated a great deal of interest in blazars as sources of high-energy neutrinos. In this study, we analyze the publicly available portion of the IceCube dataset, performing searches for neutrino point sources in spatial coincidence with the blazars and other active galactic nuclei contained in the Fermi 3LAC and the Roma BZCAT catalogs, as well as in spatial and temporal coincidence with flaring sources identified in the Fermi Collaboration's All-Sky Variability Analysis (FAVA). We find no evidence that blazars generate a significant flux of high-energy neutrinos, and conclude that no more than 5-15% of the diffuse flux measured by IceCube can originate from this class of objects. While we cannot rule out the possibility that TXS 0506+056 has at times generated significant neutrino emission, we find that such behavior cannot be common among blazars, requiring TXS 0506+056 to be a rather extreme outlier and not representative of the overall blazar population. The bulk of the diffuse high-energy neutrino flux must instead be generated by a significantly larger population of less-luminous sources, such as non-blazar active galactic nuclei.

We propose a class of models in which a stable inflaton is produced as a thermal relic in the early universe and constitutes the dark matter. We show that inflaton annihilations can efficiently reheat the universe, and identify several examples of inflationary potentials that can accommodate all cosmic microwave background observables and in which the inflaton dark matter candidate has a weak scale mass. As a simple example, we consider annihilations that take place through a Higgs portal interaction, leading to encouraging prospects for future direct detection experiments.

The presence of dark energy in our universe is causing space to expand at an accelerating rate. As a result, over the next approximately 100 billion years, all stars residing beyond the Local Group will fall beyond the cosmic horizon and become not only unobservable, but entirely inaccessible, thus limiting how much energy could one day be extracted from them. Here, we consider the likely response of a highly advanced civilization to this situation. In particular, we argue that in order to maximize its access to useable energy, a sufficiently advanced civilization would chose to expand rapidly outward, build Dyson Spheres or similar structures around encountered stars, and use the energy that is harnessed to accelerate those stars away from the approaching horizon and toward the center of the civilization. We find that such efforts will be most effective for stars with masses in the range of $M\sim (0.2-1) M_{\odot}$, and could lead to the harvesting of stars within a region extending out to several tens of Mpc in radius, potentially increasing the total amount of energy that is available to a future civilization by a factor of several thousand. We also discuss the observable signatures of a civilization elsewhere in the universe that is currently in this state of stellar harvesting.

Observations by HAWC indicate that many young pulsars (including Geminga and Monogem) are surrounded by spatially extended, multi-TeV emitting regions. It is not currently known, however, whether TeV emission is also produced by recycled, millisecond pulsars (MSPs). In this study, we perform a stacked analysis of 24 MSPs within HAWC's field-of-view, finding between 2.6-3.2 sigma evidence that these sources are, in fact, surrounded by TeV halos. The efficiency with which these MSPs produce TeV halos is similar to that exhibited by young pulsars. This result suggests that several dozen MSPs will ultimately be detectable by HAWC, including many "invisible" pulsars without radio beams oriented in our direction. The TeV halos of unresolved MSPs could also dominate the TeV-scale diffuse emission observed at high galactic latitudes. We also discuss the possibility that TeV and radio observations could be used to constrain the population of MSPs that is present in the inner Milky Way, thereby providing us with a new way to test the hypothesis that MSPs are responsible for the Galactic Center GeV excess.

Due to shielding, direct detection experiments are in some cases insensitive to dark matter candidates with very large scattering cross sections with nucleons. In this paper, we revisit this class of models, and derive a simple analytic criterion for conservative but robust direct detection limits. While large spin-independent cross sections seem to be ruled out, we identify potentially viable parameter space for dark matter with a spin-dependent cross section with nucleons in the range of $10^{-27} {\rm cm}^2 < \sigma_{{\rm DM}-p} < 10^{-24} \, {\rm cm}^{2}$. With these parameters, cosmic-ray scattering with dark matter in the extended halo of the Milky Way could generate a novel and distinctive gamma-ray signal at high galactic latitudes. Such a signal could be observable by Fermi or future space-based gamma-ray telescopes.

The EDGES Collaboration has recently reported the detection of a stronger-than-expected absorption feature in the global 21-cm spectrum, centered at a frequency corresponding to a redshift of z ~ 17. This observation has been interpreted as evidence that the gas was cooled during this era as a result of scattering with dark matter. In this study, we explore this possibility, applying constraints from the cosmic microwave background, light element abundances, Supernova 1987A, and a variety of laboratory experiments. After taking these constraints into account, we find that the vast majority of the parameter space capable of generating the observed 21-cm signal is ruled out. The only range of models that remains viable is that in which a small fraction, ~ 0.3-2%, of the dark matter consists of particles with a mass of ~ 10-80 MeV and which couple to the photon through a small electric charge, epsilon ~ 10^{-6}-10^{-4}. Furthermore, in order to avoid being overproduced in the early universe, such models must be supplemented with an additional depletion mechanism, such as annihilations through a L_{\mu}-L_{\tau} gauge boson or annihilations to a pair of rapidly decaying hidden sector scalars.

Measurements of the nearby pulsars Geminga and B0656+14 by the HAWC and Milagro telescopes have revealed the presence of bright TeV-emitting halos surrounding these objects. If young and middle-aged pulsars near the Galactic Center transfer a similar fraction of their energy into TeV photons, then these sources could plausibly dominate the emission that is observed by HESS and other ground-based telescopes from the innermost $\sim$$10^2$ parsecs of the Milky Way. In particular, both the spectral shape and the angular extent of this emission is consistent with TeV halos produced by a population of pulsars, although the reported correlation of this emission with the distribution of molecular gas suggests that diffuse hadronic processes also must contribute. The overall flux of this emission requires a birth rate of ~100-1000 neutron stars per Myr near the Galactic Center, in good agreement with recent estimates.

The gamma-ray annihilation and decay products of very heavy dark matter particles can undergo attenuation through pair production, leading to the development of electromagnetic cascades. This has a significant impact not only on the spectral shape of the gamma-ray signal, but also on the angular distribution of the observed photons. Such phenomena are particularly important in light of the new HAWC experiment, which provides unprecedented sensitivity to multi-TeV photons and thus to very heavy dark matter particles. In this study, we focus on dark matter in the 100 TeV-100 PeV mass range, and calculate the spectral and angular distribution of gamma-rays from dwarf galaxies and from nearby galaxy clusters in this class of models.

The HAWC Collaboration has recently reported the detection of bright and spatially extended multi-TeV gamma-ray emission from Geminga, Monogem, and a handful of other nearby, middle-aged pulsars. The angular profile of the emission observed from these pulsars is surprising, in that it implies that cosmic-ray diffusion is significantly inhibited within ~25 pc of these objects, compared to the expectations of standard Galactic diffusion models. This raises the important question of whether the diffusion coefficient in the local interstellar medium is also low, or whether it is instead better fit by the mean Galactic value. Here, we utilize recent observations of the cosmic-ray electron spectrum (extending up to ~20 TeV) by the H.E.S.S. Collaboration to show that the local diffusion coefficient cannot be as low as it is in the regions surrounding Geminga and Monogem. Instead, we conclude that cosmic rays efficiently diffuse through the bulk of the local interstellar medium. Among other implications, this further supports the conclusion that pulsars significantly contribute to the observed positron excess.

Stringent constraints from direct detection experiments and the Large Hadron Collider motivate us to consider models in which the dark matter does not directly couple to the Standard Model, but that instead annihilates into hidden sector particles which ultimately decay through small couplings to the Standard Model. We calculate the gamma-ray emission generated within the context of several such hidden sector models, including those in which the hidden sector couples to the Standard Model through the vector portal (kinetic mixing with Standard Model hypercharge), through the Higgs portal (mixing with the Standard Model Higgs boson), or both. In each case, we identify broad regions of parameter space in which the observed spectrum and intensity of the Galactic Center gamma-ray excess can easily be accommodated, while providing an acceptable thermal relic abundance and remaining consistent with all current constraints. We also point out that cosmic-ray antiproton measurements could potentially discriminate some hidden sector models from more conventional dark matter scenarios.

The $\textit{Fermi}$-LAT Collaboration recently presented a new catalog of gamma-ray sources located within the $40^{\circ} \times 40^{\circ}$ region around the Galactic Center~(Ajello et al. 2017) -- the Second Fermi Inner Galaxy (2FIG) catalog. Utilizing this catalog, they analyzed models for the spatial distribution and luminosity function of sources with a pulsar-like gamma-ray spectrum. Ajello et al. 2017 v1 also claimed to detect, in addition to a disk-like population of pulsar-like sources, an approximately 7$\sigma$ preference for an additional centrally concentrated population of pulsar-like sources, which they referred to as a "Galactic Bulge" population. Such a population would be of great interest, as it would support a pulsar interpretation of the gamma-ray excess that has long been observed in this region. In an effort to further explore the implications of this new source catalog, we attempted to reproduce the results presented by the $\textit{Fermi}$-LAT Collaboration, but failed to do so. Mimicking as closely as possible the analysis techniques undertaken in Ajello et al. 2017, we instead find that our likelihood analysis favors a very different spatial distribution and luminosity function for these sources. Most notably, our results do not exhibit a strong preference for a "Galactic Bulge" population of pulsars. Furthermore, we find that masking the regions immediately surrounding each of the 2FIG pulsar candidates does $\textit{not}$ significantly impact the spectrum or intensity of the Galactic Center gamma-ray excess. Although these results refute the claim of strong evidence for a centrally concentrated pulsar population presented in Ajello et al. 2017, they neither rule out nor provide support for the possibility that the Galactic Center excess is generated by a population of low-luminosity and currently largely unobserved pulsars.

Annihilating dark matter particles in nearby subhalos could generate potentially observable fluxes of gamma rays, unaccompanied by emission at other wavelengths. Furthermore, this gamma-ray emission is expected to be spatially extended, providing us with a powerful way to discriminate dark matter subhalos from other astrophysical gamma-ray sources. Fermi has detected two dark matter subhalo candidates which exhibit a statistically significant degree of spatial extension (3FGL J2212.5+0703 and 3FGL J1924.8-1034). It has been argued that the most likely non-dark matter interpretation of these observations is that they are each in fact multiple nearby point sources, too close to one another on the sky to be individually resolved. In this study, we consider the ability of next generation gamma-ray telescopes to spatially resolve the gamma-ray emission from subhalo candidates, focusing on the proposed e-ASTROGAM mission. We find that such an instrument could significantly clarify the nature of Fermi's dark matter subhalo candidates, and provide an unprecedented level of sensitivity to the presence of annihilating dark matter in nearby subhalos.

The antiproton-to-proton ratio in the cosmic-ray spectrum is a sensitive probe of new physics. Using recent measurements of the cosmic-ray antiproton and proton fluxes in the energy range of 1-1000 GeV, we study the contribution to the $\bar{p}/p$ ratio from secondary antiprotons that are produced and subsequently accelerated within individual supernova remnants. We consider several well-motivated models for cosmic-ray propagation in the interstellar medium and marginalize our results over the uncertainties related to the antiproton production cross section and the time-, charge-, and energy-dependent effects of solar modulation. We find that the increase in the $\bar{p}/p$ ratio observed at rigidities above $\sim$ 100 GV cannot be accounted for within the context of conventional cosmic-ray propagation models, but is consistent with scenarios in which cosmic-ray antiprotons are produced and subsequently accelerated by shocks within a given supernova remnant. In light of this, the acceleration of secondary cosmic rays in supernova remnants is predicted to substantially contribute to the cosmic-ray positron spectrum, accounting for a significant fraction of the observed positron excess.

If millisecond pulsars (MSPs) are responsible for the excess gamma-ray emission observed from the region surrounding the Galactic Center, the same region should also contain a large population of low-mass X-ray binaries (LMXBs). In this study, we compile and utilize a sizable catalog of LMXBs observed in the the Milky Way's globular cluster system and in the Inner Galaxy, as well as the gamma-ray emission observed from globular clusters, to estimate the flux of gamma rays predicted from MSPs in the Inner Galaxy. From this comparison, we conclude that only up to $\sim$4-23% of the observed gamma-ray excess is likely to originate from MSPs. This result is consistent with, and more robust than, previous estimates which utilized smaller samples of both globular clusters and LMXBs. If MSPs had been responsible for the entirety of the observed excess, INTEGRAL should have detected $\sim$$10^3$ LMXBs from within a $10^{\circ}$ radius around the Galactic Center, whereas only 42 LMXBs (and 46 additional LMXB candidates) have been observed.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

Utilizing an exhaustive set of simplified models, we revisit dark matter scenarios potentially capable of generating the observed Galactic Center gamma-ray excess, updating constraints from the LUX and PandaX-II experiments, as well as from the LHC and other colliders. We identify a variety of pseudoscalar mediated models that remain consistent with all constraints. In contrast, dark matter candidates which annihilate through a spin-1 mediator are ruled out by direct detection constraints unless the mass of the mediator is near an annihilation resonance, or the mediator has a purely vector coupling to the dark matter and a purely axial coupling to Standard Model fermions. All scenarios in which the dark matter annihilates through $t$-channel processes are now ruled out by a combination of the constraints from LUX/PandaX-II and the LHC.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

Utilizing data from the ELVIS and Via Lactea-II simulations, we characterize the local dark matter subhalo population, and use this information to refine the predictions for the gamma-ray fluxes arising from annihilating dark matter in this class of objects. We find that the shapes of nearby subhalos are significantly altered by tidal effects, and are generally not well described by NFW density profiles, instead prefering power-law profiles with an exponential cutoff. From the subhalo candidates detected by the Fermi Gamma-Ray Space Telescope, we place limits on the dark matter annihilation cross section that are only modestly weaker than those based on observations of dwarf galaxies. We also calculate the fraction of observable subhalos that are predicted to be spatially extended at a level potentially discernible to Fermi.

Observations by HAWC and Milagro have detected bright and spatially extended TeV gamma-ray sources surrounding the Geminga and Monogem pulsars. We argue that these observations, along with a substantial population of other extended TeV sources coincident with pulsar wind nebulae, constitute a new morphological class of spatially extended TeV halos. We show that HAWCs wide field-of-view unlocks an expansive parameter space of TeV halos not observable by atmospheric Cherenkov telescopes. Under the assumption that Geminga and Monogem are typical middle-aged pulsars, we show that ten-year HAWC observations should eventually observe 37$^{+17}_{-13}$ middle-aged TeV halos that correspond to pulsars whose radio emission is not beamed towards Earth. Depending on the extrapolation of the TeV halo efficiency to young pulsars, HAWC could detect more than 100 TeV halos from mis-aligned pulsars. These pulsars have historically been difficult to detect with existing multiwavelength observations. TeV halos will constitute a significant fraction of all HAWC sources, allowing follow-up observations to efficiently find pulsar wind nebulae and thermal pulsar emission. The observation and subsequent multi-wavelength follow-up of TeV halos will have significant implications for our understanding of pulsar beam geometries, the evolution of PWN, the diffusion of cosmic-rays near energetic pulsars, and the contribution of pulsars to the cosmic-ray positron excess.

Recent measurements of the Geminga and B0656+14 pulsars by the gamma-ray telescope HAWC (along with earlier measurements by Milagro) indicate that these objects generate significant fluxes of very high-energy electrons. In this paper, we use the very high-energy gamma-ray intensity and spectrum of these pulsars to calculate and constrain their expected contributions to the local cosmic-ray positron spectrum. Among models that are capable of reproducing the observed characteristics of the gamma-ray emission, we find that pulsars invariably produce a flux of high-energy positrons that is similar in spectrum and magnitude to the positron fraction measured by PAMELA and AMS-02. In light of this result, we conclude that it is very likely that pulsars provide the dominant contribution to the long perplexing cosmic-ray positron excess.

It has recently been shown that if the dark matter is in thermal equilibrium with a sector that is highly decoupled from the Standard Model, it can freeze-out with an acceptable relic abundance, even if the dark matter is as heavy as ~1-100 PeV. In such scenarios, both the dark and visible sectors are populated after inflation, but with independent temperatures. The lightest particle in the dark sector will be generically long-lived, and can come to dominate the energy density of the universe. Upon decaying, these particles can significantly reheat the visible sector, diluting the abundance of dark matter and thus allowing for dark matter particles that are much heavier than conventional WIMPs. In this paper, we present a systematic and pedagogical treatment of the cosmological history in this class of models, emphasizing the simplest scenarios in which a dark matter candidate annihilates into hidden sector particles which then decay into visible matter through the vector, Higgs, or lepton portals. In each case, we find ample parameter space in which very heavy dark matter particles can provide an acceptable thermal relic abundance. We also discuss possible extensions of models featuring these dynamics.

An important factor limiting our ability to understand the production and propagation of cosmic rays pertains to the effects of heliospheric forces, commonly known as solar modulation. The solar wind is capable of generating time and charge-dependent effects on the spectrum and intensity of low energy ($\lsim$ 10 GeV) cosmic rays reaching Earth. Previous analytic treatments of solar modulation have utilized the force-field approximation, in which a simple potential is adopted whose amplitude is selected to best fit the cosmic-ray data taken over a given period of time. Making use of recently available cosmic-ray data from the Voyager 1 spacecraft, along with measurements of the heliospheric magnetic field and solar wind, we construct a time, charge and rigidity-dependent model of solar modulation that can be directly compared to data from a variety of cosmic-ray experiments. We provide a simple analytic formula that can be easily utilized in a variety of applications, allowing us to better predict the effects of solar modulation and reduce the number of free parameters involved in cosmic ray propagation models.

In recent years, direct detection, indirect detection, and collider experiments have placed increasingly stringent constraints on particle dark matter, exploring much of the parameter space associated with the WIMP paradigm. In this paper, we focus on the subset of WIMP models in which the dark matter annihilates in the early universe through couplings to either the Standard Model $Z$ or the Standard Model Higgs boson. Considering fermionic, scalar, and vector dark matter candidates within a model-independent context, we find that the overwhelming majority of these dark matter candidates are already ruled out by existing experiments. In the case of $Z$ mediated dark matter, the only scenarios that are not currently excluded are those in which the dark matter is a fermion with an axial coupling and with a mass either within a few GeV of the $Z$ resonance ($m_{\rm DM} \simeq m_Z/2$) or greater than 200 GeV, or with a vector coupling and with $m_{\rm DM} > 6$ TeV. Several Higgs mediated scenarios are currently viable if the mass of the dark matter is near the Higgs pole ($m_{\rm DM} \simeq m_H/2$). Otherwise, the only scenarios that are not excluded are those in which the dark matter is a scalar (vector) heavier than 400 GeV (1160 GeV) with a Higgs portal coupling, or a fermion with a pseudoscalar (CP violating) coupling to the Standard Model Higgs boson. With the exception of dark matter with a purely pseudoscalar coupling to the Higgs, it is anticipated that planned direct detection experiments will probe nearly the entire range of models considered in this study.

Sterile neutrinos can be generated in the early universe through oscillations with active neutrinos and represent a popular and well-studied candidate for our universe's dark matter. Stringent constraints from X-ray and gamma-ray line searches, however, have excluded the simplest of such models. In this letter, we propose a novel alternative to the standard scenario in which the mixing angle between the sterile and active neutrinos is a dynamical quantity, induced through interactions with a light axion-like field. As the energy density of the axion-like particles is diluted by Hubble expansion, the degree of mixing is reduced at late times, suppressing the decay rate and easily alleviating any tension with X-ray or gamma-ray constraints. We present a simple model which illustrates the phenomenology of this scenario, and also describe a framework in which the QCD axion is responsible for the production of sterile neutrinos in the early universe.

It has been suggested that unresolved radio galaxies and radio quasars (sometimes referred to as misaligned active galactic nuclei) could be responsible for a significant fraction of the observed diffuse gamma-ray background. In this study, we use the latest data from the Fermi Gamma-Ray Space Telescope to characterize the gamma-ray emission from a sample of 51 radio galaxies. In addition to those sources that had previously been detected using Fermi data, we report here the first statistically significant detection of gamma-ray emission from the radio galaxies 3C 212, 3C 411, and B3 0309+411B. Combining this information with the radio fluxes, radio luminosity function, and redshift distribution of this source class, we find that radio galaxies dominate the diffuse gamma-ray background, generating $77.2^{+25.4}_{-9.4}\%$ of this emission at energies above $\sim$1 GeV. We discuss the implications of this result and point out that it provides support for scenarios in which IceCube's high-energy astrophysical neutrinos also originate from the same population of radio galaxies.

A recent study, making use of the number of horizontal branch stars observed in infrared photometric surveys and kinematic measurements of M-giant stars from the BRAVA survey, combined with N-body simulations of stellar populations, has presented a new determination of the dark matter mass within the bulge-bar region of the Milky Way. That study constrains the total mass within the $\pm 2.2 \times \pm 1.4 \times \pm 1.2$ kpc volume of the bulge-bar region to be ($1.84 \pm 0.07) \times 10^{10} \, M_{\odot}$, of which 9-30% is made up of dark matter. Here, we use this result to constrain the the Milky Way's dark matter density profile, and discuss the implications for indirect dark matter searches. Although uncertainties remain significant, these results favor dark matter distributions with a cusped density profile. For example, for a scale radius of 20 kpc and a local dark matter density of 0.4 GeV/cm$^3$, density profiles with an inner slope of 0.69 to 1.40 are favored, approximately centered around the standard NFW value. In contrast, profiles with large flat-density cores are disfavored by this information.

It has been suggested that the GeV excess, observed from the region surrounding the Galactic Center, might originate from a population of millisecond pulsars that formed in globular clusters. With this in mind, we employ the publicly available Fermi data to study the gamma-ray emission from 157 globular clusters, identifying a statistically significant signal from 25 of these sources (ten of which are not found in existing gamma-ray catalogs). We combine these observations with the predicted pulsar formation rate based on the stellar encounter rate of each globular cluster to constrain the gamma-ray luminosity function of millisecond pulsars in the Milky Way's globular cluster system. We find that this pulsar population exhibits a luminosity function that is quite similar to those millisecond pulsars observed in the field of the Milky Way (i.e. the thick disk). After pulsars are expelled from a globular cluster, however, they continue to lose rotational kinetic energy and become less luminous, causing their luminosity function to depart from the steady-state distribution. Using this luminosity function and a model for the globular cluster disruption rate, we show that millisecond pulsars born in globular clusters can account for only a few percent or less of the observed GeV excess. Among other challenges, scenarios in which the entire GeV excess is generated from such pulsars are in conflict with the observed mass of the Milky Way's Central Stellar Cluster.

Motivated by the recently reported diboson and dijet excesses in Run 1 data at ATLAS and CMS, we explore models of mixed dark matter in left-right symmetric theories. In this study, we calculate the relic abundance and the elastic scattering cross section with nuclei for a number of dark matter candidates that appear within the fermionic multiplets of left-right symmetric models. In contrast to the case of pure multiplets, WIMP-nucleon scattering proceeds at tree-level, and hence the projected reach of future direct detection experiments such as LUX-ZEPLIN and XENON1T will cover large regions of parameter space for TeV-scale thermal dark matter. Decays of the heavy charged W' boson to particles in the dark sector can potentially shift the right-handed gauge coupling to larger values when fixed to the rate of the Run 1 excesses, moving towards the theoretically attractive scenario, gR = gL. This region of parameter space may be probed by future collider searches for new Higgs bosons or electroweak fermions.

We present an argument that radio galaxies (active galaxies with mis-aligned jets) are likely to be the primary sources of the high-energy astrophysical neutrinos observed by IceCube. In particular, if the gamma-ray emission observed from radio galaxies is generated through the interactions of cosmic-ray protons with gas, these interactions can also produce a population of neutrinos with a flux and spectral shape similar to that measured by IceCube. We present a simple physical model in which high-energy cosmic rays are confined within the volumes of radio galaxies, where they interact with gas to generate the observed diffuse fluxes of neutrinos and gamma rays. In addition to simultaneously accounting for the observations of Fermi and IceCube, radio galaxies in this model also represent an attractive class of sources for the highest energy cosmic rays.

Although dark matter is a central element of modern cosmology, the history of how it became accepted as part of the dominant paradigm is often ignored or condensed into a brief anecdotical account focused around the work of a few pioneering scientists. The aim of this review is to provide the reader with a broader historical perspective on the observational discoveries and the theoretical arguments that led the scientific community to adopt dark matter as an essential part of the standard cosmological model.

We present results from a multi-chord Pluto stellar occultation observed on 29 June 2015 from New Zealand and Australia. This occurred only two weeks before the NASA New Horizons flyby of the Pluto system and serves as a useful comparison between ground-based and space results. We find that Pluto's atmosphere is still expanding, with a significant pressure increase of 5$\pm$2\% since 2013 and a factor of almost three since 1988. This trend rules out, as of today, an atmospheric collapse associated with Pluto's recession from the Sun. A central flash, a rare occurrence, was observed from several sites in New Zealand. The flash shape and amplitude are compatible with a spherical and transparent atmospheric layer of roughly 3~km in thickness whose base lies at about 4~km above Pluto's surface, and where an average thermal gradient of about 5 K~km$^{-1}$ prevails. We discuss the possibility that small departures between the observed and modeled flash are caused by local topographic features (mountains) along Pluto's limb that block the stellar light. Finally, using two possible temperature profiles, and extrapolating our pressure profile from our deepest accessible level down to the surface, we obtain a possible range of 11.9-13.7~$\mu$bar for the surface pressure.

Millisecond pulsars have been discussed as a possible source of the gamma-ray excess observed from the region surrounding the Galactic Center. With this in mind, we use the observed population of bright low-mass X-ray binaries to estimate the number of millisecond pulsars in the Inner Galaxy. This calculation suggests that only ~1-5% of the excess is produced by millisecond pulsars. We also use the luminosity function derived from local measurements of millisecond pulsars, along with the number of point sources resolved by Fermi, to calculate an upper limit for the diffuse emission from such a population. While this limit is compatible with the millisecond pulsar population implied by the number of low-mass X-ray binaries, it strongly excludes the possibility that most of the excess originates from such objects.

In this letter, we consider a class of scenarios in which the dark matter is part of a heavy hidden sector that is thermally decoupled from the Standard Model in the early universe. The dark matter freezes-out by annihilating to a lighter, metastable state, whose subsequent abundance can naturally come to dominate the energy density of the universe. When this state decays, it reheats the visible sector and dilutes all relic abundances, thereby allowing the dark matter to be orders of magnitude heavier than the weak scale. For concreteness, we consider a simple realization with a Dirac fermion dark matter candidate coupled to a massive gauge boson that decays to the Standard Model through its kinetic mixing with hypercharge. We identify viable parameter space in which the dark matter can be as heavy as ~1-100 PeV without being overproduced in the early universe.

In a previous paper, we pointed out that the gamma-ray source 3FGL J2212.5+0703 shows evidence of being spatially extended. If a gamma-ray source without detectable emission at other wavelengths were unambiguously determined to be spatially extended, it could not be explained by known astrophysics, and would constitute a smoking gun for dark matter particles annihilating in a nearby subhalo. With this prospect in mind, we scrutinize the gamma-ray emission from this source, finding that it prefers a spatially extended profile over that of a single point-like source with 5.1 sigma statistical significance. We also use a large sample of active galactic nuclei and other known gamma-rays sources as a control group, confirming, as expected, that statistically significant extension is rare among such objects. We argue that the most likely (non-dark matter) explanation for this apparent extension is a pair of bright gamma-ray sources that serendipitously lie very close to each other, and estimate that there is a chance probability of ~2% that such a pair would exist somewhere on the sky. In the case of 3FGL J2212.5+0703, a model with a second gamma-ray point source at the location of a known BZCAT/CRATES radio source yields fits that are comparable in quality to those obtained for a single extended source. If 3FGL J2212.5+0703 is a dark matter subhalo, it would imply that dark matter particles have a mass of ~18-33 GeV and an annihilation cross section on the order of sigma v ~ 10^-26 cm^3/s (for the representative case of annihilations to bb), similar to the values required to generate the Galactic Center gamma-ray excess.

A number of groups have employed radio observations of the Galactic center to derive stringent constraints on the annihilation cross section of weakly interacting dark matter. In this paper, we show that electron energy losses in this region are likely to be dominated by inverse Compton scattering on the interstellar radiation field, rather than by synchrotron, considerably relaxing the constraints on the dark matter annihilation cross section compared to previous works. Strong convective winds, which are well motivated by recent observations, may also significantly weaken synchrotron constraints. After taking these factors into account, we find that radio constraints on annihilating dark matter are orders of magnitude less stringent than previously reported, and are generally weaker than those derived from current gamma-ray observations.

It has been proposed that a large population of unresolved millisecond pulsars (MSPs) could potentially account for the excess of GeV-scale gamma-rays observed from the region surrounding the Galactic Center. The viability of this scenario depends critically on the gamma-ray luminosity function of this source population, which determines how many MSPs Fermi should have already detected as resolved point sources. In this paper, we revisit the gamma-ray luminosity function of MSPs, without relying on uncertain distance measurements. Our determination, based on a comparison of models with the observed characteristics of the MSP population, suggests that Fermi should have already detected a significant number of sources associated with such a hypothesized Inner Galaxy population. We cannot rule out a scenario in which the MSPs residing near the Galactic Center are systematically less luminous than those present in the Galactic Plane or within globular clusters.

We revisit the possibility of constraining the properties of dark matter (DM) by studying the epoch of cosmic reionization. Previous studies have shown that DM annihilation was unlikely to have provided a large fraction of the photons that ionized the universe, but instead played a subdominant role relative to stars and quasars. The DM, however, begins to efficiently annihilate with the formation of primordial microhalos at $z\sim100-200$, much earlier than the formation of the first stars. Therefore, if DM annihilation ionized the universe at even the percent level over the interval $z \sim 20-100$, it can leave a significant imprint on the global optical depth, $\tau$. Moreover, we show that cosmic microwave background (CMB) polarization data and future 21 cm measurements will enable us to more directly probe the DM contribution to the optical depth. In order to compute the annihilation rate throughout the epoch of reionization, we adopt the latest results from structure formation studies and explore the impact of various free parameters on our results. We show that future measurements could make it possible to place constraints on the dark matter's annihilation cross section that are at a level comparable to those obtained from the observations of dwarf galaxies, cosmic ray measurements, and studies of recombination.

Light gauge bosons can lead to resonant interactions between high-energy astrophysical neutrinos and the cosmic neutrino background. We study this possibility in detail, considering the ability of IceCube to probe such scenarios. We find the most dramatic effects in models with a very light $Z'$ ($m_{Z'} \lesssim 10$ MeV), which can induce a significant absorption feature at $E_{\nu} \sim\,$5-10$\,{\rm TeV} \times (m_{Z'}/{\rm MeV})^2$. In the case of the inverted hierarchy and a small sum of neutrino masses, such a light $Z'$ can result in a broad and deep spectral feature at $\sim\,$0.1-10$\,{\rm PeV} \times (m_{Z'}/{\rm MeV})^2$. Current IceCube data already excludes this case for a $Z'$ lighter than a few MeV and couplings greater than $g\sim10^{-4}$. We emphasize that the ratio of neutrino flavors observed by IceCube can be used to further increase their sensitivity to $Z'$ models and to other exotic physics scenarios.

Using recent measurements of the spectrum and chemical composition of the highest energy cosmic rays, we consider the sources of these particles. We find that the data strongly prefers models in which the sources of the ultra-high energy cosmic rays inject predominantly intermediate mass nuclei, with comparatively few protons or heavy nuclei, such as iron or silicon. If the number density of sources per comoving volume does not evolve with redshift, the injected spectrum must be very hard ($\alpha\simeq 1$) in order to fit the spectrum observed at Earth. Such a hard spectral index would be surprising and difficult to accommodate theoretically. In contrast, much softer spectral indices, consistent with the predictions of Fermi acceleration ($\alpha\simeq 2$), are favored in models with negative source evolution. With this theoretical bias, these observations thus favor models in which the sources of the highest energy cosmic rays are preferentially located within the low-redshift universe.

We describe a general scenario, dubbed "Inflatable Dark Matter", in which the density of dark matter particles can be reduced through a short period of late-time inflation in the early universe. The overproduction of dark matter that is predicted within many otherwise well-motivated models of new physics can be elegantly remedied within this context, without the need to tune underlying parameters or to appeal to anthropic considerations. Thermal relics that would otherwise be disfavored can easily be accommodated within this class of scenarios, including dark matter candidates that are very heavy or very light. Furthermore, the non-thermal abundance of GUT or Planck scale axions can be brought to acceptable levels, without invoking anthropic tuning of initial conditions. A period of late-time inflation could have occurred over a wide range of scales from ~ MeV to the weak scale or above, and could have been triggered by physics within a hidden sector, with small but not necessarily negligible couplings to the Standard Model.

In this report we summarize the many dark matter searches currently being pursued through four complementary approaches: direct detection, indirect detection, collider experiments, and astrophysical probes. The essential features of broad classes of experiments are described, each with their own strengths and weaknesses. The complementarity of the different dark matter searches is discussed qualitatively and illustrated quantitatively in two simple theoretical frameworks. Our primary conclusion is that the diversity of possible dark matter candidates requires a balanced program drawing from all four approaches.

It has been proposed that a recent outburst of cosmic-ray electrons could account for the excess of GeV-scale gamma rays observed from the region surrounding the Galactic Center. After studying this possibility in some detail, we identify scenarios in which a series of leptonic cosmic-ray outbursts could plausibly generate the observed excess. The morphology of the emission observed outside of $\sim1^{\circ}-2^{\circ}$ from the Galactic Center can be accommodated with two outbursts, one which took place approximately $\sim10^6$ years ago, and another (injecting only about 10\% as much energy as the first) about $\sim10^5$ years ago. The emission observed from the innermost $\sim1^{\circ}-2^{\circ}$ requires one or more additional recent outbursts and/or a contribution from a centrally concentrated population of unresolved millisecond pulsars. In order to produce a spectrum that is compatible with the measured excess (whose shape is approximately uniform over the region of the excess), the electrons from the older outburst must be injected with significantly greater average energy than those injected more recently, enabling their spectra to be similar after $\sim10^6$ years of energy losses.

Motivated by the gamma-ray excess observed from the region surrounding the Galactic Center, we explore particle dark matter models that could potentially account for the spectrum and normalization of this signal. Taking a model-independent approach, we consider an exhaustive list of tree-level diagrams for dark matter annihilation, and determine which could account for the observed gamma-ray emission while simultaneously predicting a thermal relic abundance equal to the measured cosmological dark matter density. We identify a wide variety of models that can meet these criteria without conflicting with existing constraints from direct detection experiments or the Large Hadron Collider (LHC). The prospects for detection in near future dark matter experiments and/or the upcoming 14 TeV LHC appear quite promising.

Dark matter annihilations taking place in nearby subhalos could appear as gamma-ray sources without detectable counterparts at other wavelengths. In this study, we consider the collection of unassociated gamma-ray sources reported by the Fermi Collaboration in an effort to identify the most promising dark matter subhalo candidates. While we identify 24 bright, high-latitude, non-variable sources with spectra that are consistent with being generated by the annihilations of ~20-70 GeV dark matter particles (assuming annihilations to $b\bar{b}$), it is not possible at this time to distinguish these sources from radio-faint gamma-ray pulsars. Deeper multi-wavelength observations will be essential to clarify the nature of these sources. It is notable that we do not find any such sources that are well fit by dark matter particles heavier than ~100 GeV. We also study the angular distribution of the gamma-rays from this set of subhalo candidates, and find that the source 3FGL J2212.5+0703 prefers a spatially extended profile (of width ~0.15$^{\circ}$) over that of a point source, with a significance of 4.2$\sigma$ (3.6$\sigma$ after trials factor). Although not yet definitive, this bright and high-latitude gamma-ray source is well fit as a nearby subhalo of $m_{\chi} \simeq$ 20-50 GeV dark matter particles (annihilating to $b\bar{b}$) and merits further multi-wavelength investigation. Based on the subhalo distribution predicted by numerical simulations, we derive constraints on the dark matter annihilation cross section that are competitive to those resulting from gamma-ray observations of dwarf spheroidal galaxies, the Galactic Center, and the extragalactic gamma-ray background.

The recent discovery of ten new dwarf galaxy candidates by the Dark Energy Survey (DES) and the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) could increase the Fermi Gamma-Ray Space Telescope's sensitivity to annihilating dark matter particles, potentially enabling a definitive test of the dark matter interpretation of the long-standing Galactic Center gamma-ray excess. In this paper, we compare the previous analyses of Fermi data from the directions of the new dwarf candidates (including the relatively nearby Reticulum II) and perform our own analysis, with the goal of establishing the statistical significance of any gamma-ray signal from these sources. We confirm the presence of an excess from Reticulum II, with a spectral shape that is compatible with the Galactic Center signal. The significance of this emission is greater than that observed from 99.84% of randomly chosen high-latitude blank-sky locations, corresponding to a local detection significance of 3.2 sigma. We improve upon the standard blank-sky calibration approach through the use of multi-wavelength catalogs, which allow us to avoid regions that are likely to contain unresolved gamma-ray sources.

Past studies have identified a spatially extended excess of $\sim$1-3 GeV gamma rays from the region surrounding the Galactic Center, consistent with the emission expected from annihilating dark matter. We revisit and scrutinize this signal with the intention of further constraining its characteristics and origin. By applying cuts to the \textit{Fermi} event parameter CTBCORE, we suppress the tails of the point spread function and generate high resolution gamma-ray maps, enabling us to more easily separate the various gamma-ray components. Within these maps, we find the GeV excess to be robust and highly statistically significant, with a spectrum, angular distribution, and overall normalization that is in good agreement with that predicted by simple annihilating dark matter models. For example, the signal is very well fit by a 36-51 GeV dark matter particle annihilating to $b\bar{b}$ with an annihilation cross section of $\sigma v = (1-3)\times 10^{-26}$ cm$^3$/s (normalized to a local dark matter density of 0.4 GeV/cm$^3$). Furthermore, we confirm that the angular distribution of the excess is approximately spherically symmetric and centered around the dynamical center of the Milky Way (within $\sim$$0.05^{\circ}$ of Sgr A$^*$), showing no sign of elongation along the Galactic Plane. The signal is observed to extend to at least $\simeq10^{\circ}$ from the Galactic Center, disfavoring the possibility that this emission originates from millisecond pulsars.

Among the most stringent constraints on the dark matter annihilation cross section are those derived from observations of dwarf galaxies by the Fermi Gamma-Ray Space Telescope. As current (e.g., Dark Energy Survey, DES) and future (Large Synoptic Survey Telescope, LSST) optical imaging surveys discover more of the Milky Way's ultra-faint satellite galaxies, they may increase Fermi's sensitivity to dark matter annihilations. In this study, we use a semi-analytic model of the Milky Way's satellite population to predict the characteristics of the dwarfs likely to be discovered by DES and LSST, and project how these discoveries will impact Fermi's sensitivity to dark matter. While we find that modest improvements are likely, the dwarf galaxies discovered by DES and LSST are unlikely to increase Fermi's sensitivity by more than a factor of ~2-4. However, this outlook may be conservative, given that our model underpredicts the number of ultra-faint galaxies with large potential annihilation signals actually discovered in the Sloan Digital Sky Survey. Our simulation-based approach focusing on the Milky Way satellite population demographics complements existing empirically-based estimates.

Models in which dark matter particles can scatter into a slightly heavier state which promptly decays to the lighter state and a photon (known as eXciting Dark Matter, or XDM) have been shown to be capable of generating the 3.55 keV line observed from galaxy clusters, while suppressing the flux of such a line from smaller halos, including dwarf galaxies. In most of the XDM models discussed in the literature, this up-scattering is mediated by a new light particle, and dark matter annihilations proceed into pairs of this same light state. In these models, the dark matter and mediator effectively reside within a hidden sector, without sizable couplings to the Standard Model. In this paper, we explore a model of XDM that does not include a hidden sector. Instead, the dark matter both up-scatters and annihilates through the near resonant exchange of a $\mathcal{O}(10^2)$ GeV pseudoscalar with large Yukawa couplings to the dark matter and smaller, but non-neglibile, couplings to Standard Model fermions. The dark matter and the mediator are each mixtures of Standard Model singlets and $SU(2)_W$ doublets. We identify parameter space in which this model can simultaneously generate the 3.55 keV line and the gamma-ray excess observed from the Galactic Center, without conflicting with constraints from colliders, direct detection experiments, or observations of dwarf galaxies.

With the goal of generating the Galactic Center gamma-ray excess, we revisit models in which the dark matter interacts with the Standard Model through the exchange of a new neutral gauge boson, Z'. We find several scenarios that can account for this signal, while respecting all existing constraints from colliders and direct detection experiments. In such models, the Z' either 1) couples axially to light quarks and is leptophobic, 2) couples dominantly to the third generation, or 3) is near resonance, m_Z' ~ 2 m_DM. We identify an example of an anomaly-free U(1)' that leads to an axial and leptophobic Z'. Many of the models presented here are within the reach of near future direct detection experiments, such as LUX and XENON1T.

Measurements of the cosmic ray antiproton spectrum can be used to search for contributions from annihilating dark matter and to constrain the dark matter annihilation cross section. Depending on the assumptions made regarding cosmic ray propagation in the Galaxy, such constraints can be quite stringent. We revisit this topic, utilizing a set of propagation models fit to the cosmic ray boron, carbon, oxygen and beryllium data. We derive upper limits on the dark matter annihilation cross section and find that when the cosmic ray propagation parameters are treated as nuisance parameters (as we argue is appropriate), the resulting limits are significantly less stringent than have been previously reported. We also note (as have several previous groups) that simple GALPROP-like diffusion-reacceleration models predict a spectrum of cosmic ray antiprotons that is in good agreement with PAMELA's observations above ~5 GeV, but that significantly underpredict the flux at lower energies. Although the complexity of modeling cosmic ray propagation at GeV-scale energies makes it difficult to determine the origin of this discrepancy, we consider the possibility that the excess antiprotons are the result of annihilating dark matter. Suggestively, we find that this excess is best fit for a dark matter mass of approximately 35 GeV and annihilation cross section of approximately 1e-26 cm^3/s (to b-bbar), in good agreement with the mass and cross section previously shown to be required to generate the gamma-ray excess observed from the Galactic Center.

The Fermi-LAT collaboration has studied the gamma-ray emission from a stacked population of dwarf spheroidal galaxies and used this information to set constraints on the dark matter annihilation cross section. Interestingly, their analysis uncovered an excess with a test statistic (TS) of 8.7. If interpreted naively, this constitutes a 2.95 sigma local excess (p-value=0.003), relative to the expectations of their background model. In order to further test this interpretation, the Fermi-LAT team studied a large number of blank sky locations and found TS>8.7 excesses to be more common than predicted by their background model, decreasing the significance of their dwarf excess to 2.2 sigma (p-value=0.027). We argue that these TS>8.7 blank sky locations are largely the result of unresolved blazars, radio galaxies, and starforming galaxies, and show that multi-wavelength information can be used to reduce the degree to which such sources contaminate the otherwise blank sky. In particular, we show that masking regions of the sky that lie within 1 degree of sources contained in the BZCAT or CRATES catalogs reduces the fraction of blank sky locations with TS>8.7 by more than a factor of two. Taking such multi-wavelength information into account can enable experiments such as Fermi to better characterize their backgrounds and increase their sensitivity to dark matter in dwarf galaxies, the most important of which remain largely uncontaminated by unresolved point sources. We also note that for the range of dark matter masses and annihilation cross sections currently being tested by studies of dwarf spheroidal galaxies, simulations predict that Fermi should be able to detect a significant number of dark matter subhalos. These subhalos constitute a population of sub-threshold gamma-ray point sources and represent an irreducible background for searches for dark matter annihilation in dwarf galaxies.

The gamma-ray excess observed from the Galactic Center can be interpreted as dark matter particles annihilating into Standard Model fermions with a cross section near that expected for a thermal relic. Although many particle physics models have been shown to be able to account for this signal, the fact that this particle has not yet been observed in direct detection experiments somewhat restricts the nature of its interactions. One way to suppress the dark matter's elastic scattering cross section with nuclei is to consider models in which the dark matter is part of a hidden sector. In such models, the dark matter can annihilate into other hidden sector particles, which then decay into Standard Model fermions through a small degree of mixing with the photon, Z, or Higgs bosons. After discussing the gamma-ray signal from hidden sector dark matter in general terms, we consider two concrete realizations: a hidden photon model in which the dark matter annihilates into a pair of vector gauge bosons that decay through kinetic mixing with the photon, and a scenario within the generalized NMSSM in which the dark matter is a singlino-like neutralino that annihilates into a pair of singlet Higgs bosons, which decay through their mixing with the Higgs bosons of the MSSM.

In this article, we revisit the gamma-ray emission observed from millisecond pulsars and globular clusters. Based on 5.6 years of data from the Fermi Gamma-Ray Space Telescope, we report gamma-ray spectra for 61 millisecond pulsars, finding most to be well fit by a power-law with an exponential cutoff, producing to a spectral peak near ~1-2 GeV (in $E^2 dN/dE$ units). Additionally, while most globular clusters exhibit a similar spectral shape, we identify a few with significantly softer spectra. We also determine the gamma-ray luminosity function of millisecond pulsars using the population found in the nearby field of the Milky Way, and within the globular cluster 47 Tucanae. We find that the gamma-ray emission observed from globular clusters is dominated by a relatively small number of bright millisecond pulsars, and that low-luminosity pulsars account for only a small fraction of the total flux. Our results also suggest that the gamma-ray emission from millisecond pulsars is more isotropic and less strongly beamed than the emission at X-ray wavelengths. Furthermore, the observed distribution of apparent gamma-ray efficiencies provides support for the slot gap or the outer gap models over those in which the gamma-ray emission originates from regions close to the neutron star's magnetic poles (polar cap models).

Thermal relic dark matter particles with a mass of 31-40 GeV and that dominantly annihilate to bottom quarks have been shown to provide an excellent description of the excess gamma rays observed from the center of the Milky Way. Flavored dark matter provides a well-motivated framework in which the dark matter can dominantly couple to bottom quarks in a flavor-safe manner. We propose a phenomenologically viable model of bottom flavored dark matter that can account for the spectral shape and normalization of the gamma-ray excess while naturally suppressing the elastic scattering cross sections probed by direct detection experiments. This model will be definitively tested with increased exposure at LUX and with data from the upcoming high-energy run of the Large Hadron Collider (LHC).

The rapid rise in the cosmic ray positron fraction above 10 GeV, as measured by PAMELA and AMS, suggests the existence of nearby primary sources of high energy positrons, such as pulsars or annihilating/decaying dark matter. In contrast, the spectrum of secondary positrons produced through the collisions of cosmic rays in the interstellar medium is predicted to fall rapidly with energy, and thus is unable to account for the observed rise. It has been proposed, however, that secondary positrons could be produced and then accelerated in nearby supernova remnants, potentially explaining the observed rise, without the need of primary positron sources. Yet, if secondary positrons are accelerated in such shocks, other secondary cosmic ray species (such as boron nuclei, and antiprotons) will also be accelerated, leading to rises in the boron-to-carbon and antiproton-to-proton ratios. The measurements of the boron-to-carbon ratio by the PAMELA and AMS collaborations, however, show no sign of such a rise. With this new data in hand, we revisit the secondary acceleration scenario for the rising positron fraction. Assuming that the same supernova remnants accelerate both light nuclei (protons, helium) and heavier cosmic ray species, we find that no more than ~25% of the observed rise in the positron fraction can result from this mechanism (at the 95% confidence level)

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.

These reports present the results of the 2013 Community Summer Study of the APS Division of Particles and Fields ("Snowmass 2013") on the future program of particle physics in the U.S. Chapter 4, on the Cosmic Frontier, discusses the program of research relevant to cosmology and the early universe. This area includes the study of dark matter and the search for its particle nature, the study of dark energy and inflation, and cosmic probes of fundamental symmetries.

Several classes of astrophysical sources contribute to the approximately isotropic gamma-ray background measured by the Fermi Gamma-Ray Space Telescope. In this paper, we use Fermi's catalog of gamma-ray sources (along with corresponding source catalogs at infrared and radio wavelengths) to build and constrain a model for the contributions to the extragalactic gamma-ray background from astrophysical sources, including radio galaxies, star-forming galaxies, and blazars. We then combine our model with Fermi's measurement of the gamma-ray background to derive constraints on the dark matter annihilation cross section, including contributions from both extragalactic and galactic halos and subhalos. The resulting constraints are competitive with the strongest current constraints from the Galactic Center and dwarf spheroidal galaxies. As Fermi continues to measure the gamma-ray emission from a greater number of astrophysical sources, it will become possible to more tightly constrain the astrophysical contributions to the extragalactic gamma-ray background. We project that with 10 years of data, Fermi's measurement of this background combined with the improved constraints on the astrophysical source contributions will yield a sensitivity to dark matter annihilations that exceeds the strongest current constraints by a factor of ~ 5 - 10.

In supersymmetric models with an unbroken R-symmetry (rather than only R-parity), the neutralinos are Dirac fermions rather than Majorana. In this article, we discuss the phenomenology of neutralino dark matter in such models, including the calculation of the thermal relic abundance, and constraints and prospects for direct and indirect searches. Due to the large elastic scattering cross sections with nuclei predicted in R-symmetric models, we are forced to consider a neutralino that is predominantly bino, with very little higgsino mixing. We find a large region of parameter space in which bino-like Dirac neutralinos with masses between 10 and 380 GeV can annihilate through slepton exchange to provide a thermal relic abundance in agreement with the observed cosmological density, without relying on coannihilations or resonant annihilations. The signatures for the indirect detection of Dirac neutralinos are very different than predicted in the Majorana case, with annihilations proceeding dominately to $\tau^+ \tau^-$, $\mu^+ \mu^-$ and $e^+ e^-$ final states, without the standard chirality suppression. And unlike Majorana dark matter candidates, Dirac neutralinos experience spin-independent scattering with nuclei through vector couplings (via $Z$ and squark exchange), leading to potentially large rates at direct detection experiments. These and other characteristics make Dirac neutralinos potentially interesting within the context of recent direct and indirect detection anomalies. We also discuss the case in which the introduction of a small Majorana mass term breaks the $R$-symmetry, splitting the Dirac neutralino into a pair of nearly degenerate Majorana states.

In this Report we discuss the four complementary searches for the identity of dark matter: direct detection experiments that look for dark matter interacting in the lab, indirect detection experiments that connect lab signals to dark matter in our own and other galaxies, collider experiments that elucidate the particle properties of dark matter, and astrophysical probes sensitive to non-gravitational interactions of dark matter. The complementarity among the different dark matter searches is discussed qualitatively and illustrated quantitatively in several theoretical scenarios. Our primary conclusion is that the diversity of possible dark matter candidates requires a balanced program based on all four of those approaches.

The AMS experiment onboard the International Space Station has recently provided cosmic ray electron and positron data with unprecedented precision in the range from 0.5 to 350 GeV. The observed rise in the positron fraction at energies above 10 GeV remains unexplained, with proposed solutions ranging from local pulsars to TeV-scale dark matter. Here, we make use of this high quality data to place stringent limits on dark matter with masses below ~300 GeV, annihilating or decaying to leptonic final states, essentially independent of the origin of this rise. We significantly improve on existing constraints, in some cases by up to two orders of magnitude.

The dark matter halo of the Milky Way is predicted to contain a very large number of smaller subhalos. As a result of the dark matter annihilations taking place within such objects, the most nearby and massive subhalos could appear as point-like or spatially extended gamma-ray sources, without observable counterparts at other wavelengths. In this paper, we use the results of the Aquarius simulation to predict the distribution of nearby subhalos, and compare this to the characteristics of the unidentified gamma-ray sources observed by the Fermi Gamma-Ray Space Telescope. Focusing on the brightest high latitude sources, we use this comparison to derive limits on the dark matter annihilation cross section. For dark matter particles lighter than ~200 GeV, the resulting limits are the strongest obtained to date, being modestly more stringent than those derived from observations of dwarf galaxies or the Galactic Center. We also derive independent limits based on the lack of unidentified gamma-ray sources with discernible spatial extension, but these limits are a factor of ~2-10 weaker than those based on point-like subhalos. Lastly, we note that four of the ten brightest high-latitude sources exhibit a similar spectral shape, consistent with 30-60 GeV dark matter particles annihilating to b quarks with an annihilation cross section on the order of sigma v ~ (5-10) x 10^-27 cm^3/s, or 8-10 GeV dark matter particles annihilating to taus with sigma v ~ (2.0-2.5) x 10^-27 cm^3/s.

The rise of the cosmic ray positron fraction with energy, as first observed with high confidence by PAMELA, implies that a large flux of high energy positrons has been recently (or is being currently) injected into the local volume of the Milky Way. With the new and much more precise measurement of the positron fraction recently provided by AMS, we revisit the question of the origin of these high energy positrons. We find that while some dark matter models (annihilating directly to electrons or muons) no longer appear to be capable of accommodating these data, other models in which ~1-3 TeV dark matter particles annihilate to unstable intermediate states could still be responsible for the observed signal. Nearby pulsars also remain capable of explaining the observed positron fraction. Future measurements of the positron fraction by AMS (using a larger data set), combined with their anticipated measurements of various cosmic ray secondary-to-primary ratios, may enable us to further discriminate between these remaining scenarios.

If the Galactic WMAP radio haze, as recently confirmed by Planck, is produced by dark matter annihilation or decay, similar diffuse radio halos should exist around other galaxies with physical properties comparable to the Milky Way. If instead the haze is due to an astrophysical mechanism peculiar to the Milky Way or to a transient event, a similar halo need not exist around all Milky Way "twins". We use radio observations of 66 spiral galaxies to test the dark matter origin of the haze. We select galaxies based on morphological type and maximal rotational velocity, and obtain their luminosities from a 1.49 GHz catalog and additional radio observations at other frequencies. We find many instances of galaxies with radio emission that is less than 5% as bright as naively expected from dark matter models that could produce the Milky Way haze, and at least 3 galaxies that are less than 1% as bright as expected, assuming dark matter distributions, magnetic fields, and cosmic ray propagation parameters equal to those of the Milky Way. For reasonable ranges for the variation of these parameters, we estimate the fraction of galaxies that should be expected to be significantly less bright in radio, and argue that this is marginally compatible with the observed distribution. While our findings therefore cannot rule out a dark matter origin for the radio haze at this time, we find numerous examples (including the Andromeda Galaxy) where, if dark matter is indeed the origin of the Milky Way haze, some mechanism must be in place to suppress the corresponding haze of the external galaxy. We point out that Planck data will offer opportunities to improve this type of constraint in a highly relevant frequency range and for a potentially larger set of candidate galaxies.

The IceCube collaboration has recently reported the observation of two events with energies in excess of 1 PeV. While an atmospheric origin of these events cannot be ruled out at this time, this pair of showers may potentially represent the first observation of high-energy astrophysical neutrinos. In this paper, we argue that if these events are neutrino-induced, then the neutrinos are very likely to have been produced via photo-meson interactions taking place in the same class of astrophysical objects that are responsible for the acceleration of the $\sim$$10^{17}$ eV cosmic ray spectrum. Among the proposed sources of such cosmic rays, gamma-ray bursts stand out as particularly capable of generating PeV neutrinos at the level implied by IceCube's two events. In contrast, the radiation fields in typical active galactic nuclei models are likely dominated by lower energy (UV) photons, and thus feature higher energy thresholds for pion production, leading to neutrino spectra which peak at EeV rather than PeV energies (models with significant densities of x-ray emission, however, could evade this problem). Cosmogenic neutrinos generated from the propagation of ultra-high energy cosmic rays similarly peak at energies that are much higher than those of the events reported by IceCube.

If moduli, or other long-lived heavy states, decay in the early universe in part into light and feebly interacting particles (such as axions), these decay products could account for the additional energy density in radiation that is suggested by recent measurements of the CMB. These moduli decays will also, however, alter the expansion history of the early universe, potentially diluting the thermal relic abundance of dark matter. If this is the case, then dark matter particles must annihilate with an even lower cross section than required in the standard thermal scenario (sigma v < 3x10^-26 cm^3/s) if they are to make up the observed density of dark matter. This possibility has significant implications for direct and indirect searches for dark matter.

Although observations made with the CoGeNT and CDMS experiments have been interpreted as possible signals of low-mass (~7-10 GeV) dark matter particles, constraints from the XENON100 collaboration appear to be incompatible with this hypothesis, at least at face value. In this paper, we revisit XENON100's constraint on dark matter in this mass range, and consider how various uncertainties and assumptions made might alter this conclusion. We also note that while XENON100's two nuclear recoil candidates each exhibit very low ratios of ionization-to-scintillation signals, making them difficult to attribute to known electronic or neutron backgrounds, they are consistent with originating from dark matter particles in the mass range favored by CoGeNT and CDMS. We argue that with lower, but not implausible, values for the relative scintillation efficiency of liquid xenon (L_eff), and the suppression of the scintillation signal in liquid xenon at XENON100's electric field (S_nr), these two events could consistently arise from dark matter particles with a mass and cross section in the range favored by CoGeNT and CDMS. If this interpretation is correct, we predict that the LUX experiment, with a significantly higher light yield than XENON100, should observe dark matter induced events at an observable rate of ~3-24 per month.

Using data from the Fermi Gamma-Ray Space Telescope, a spatially extended component of gamma rays has been identified from the direction of the Galactic Center, peaking at energies of ~2-3 GeV. More recently, it has been shown that this signal is not confined to the innermost hundreds of parsecs of the Galaxy, but instead extends to at least ~3 kpc from the Galactic Center. While the spectrum, intensity, and angular distribution of this signal is in good agreement with predictions from annihilating dark matter, it has also been suggested that a population of unresolved millisecond pulsars could be responsible for this excess GeV emission from the Inner Galaxy. In this paper, we consider this later possibility in detail. Comparing the observed spectral shape of the Inner Galaxy's GeV excess to the spectrum measured from 37 millisecond pulsars by Fermi, we find that these sources exhibit a spectral shape that is much too soft at sub-GeV energies to accommodate this signal. We also construct population models to describe the spatial distribution and luminosity function of the Milky Way's millisecond pulsars. After taking into account constraints from the observed distribution of Fermi sources (including both sources known to be millisecond pulsars, and unidentified sources which could be pulsars), we find that millisecond pulsars can account for no more than ~10% of the Inner Galaxy's GeV excess. Each of these arguments strongly disfavor millisecond pulsars as the source of this signal.

For any realistic halo profile, the Galactic Center is predicted to be the brightest source of gamma-rays from dark matter annihilations. Due in large part to uncertainties associated with the dark matter distribution and astrophysical backgrounds, however, the most commonly applied constraints on the dark matter annihilation cross section have been derived from other regions, such as dwarf spheroidal galaxies. In this article, we study Fermi Gamma-Ray Space Telescope data from the direction of the inner Galaxy and derive stringent upper limits on the dark matter's annihilation cross section. Even for the very conservative case of a dark matter distribution with a significant (~kpc) constant-density core, normalized to the minimum density needed to accommodate rotation curve and microlensing measurements, we find that the Galactic Center constraint is approximately as stringent as those derived from dwarf galaxies (which were derived under the assumption of an NFW distribution). For NFW or Einasto profiles (again, normalized to the minimum allowed density), the Galactic Center constraints are typically stronger than those from dwarfs.

It has been demonstrated that the inclusion of baryonic physics can alter the dark matter densities in the centers of low-mass galaxies, making the central dark matter slope more shallow than predicted in pure cold dark matter simulations. This flattening of the dark matter profile can occur in the most luminous subhalos around Milky Way-mass galaxies. Zolotov et al. (2012) have suggested a correction to be applied to the central masses of dark matter-only satellites in order to mimic the affect of (1) the flattening of the dark matter cusp due to supernova feedback in luminous satellites, and (2) enhanced tidal stripping due to the presence of a baryonic disk. In this paper, we apply this correction to the z=0 subhalo masses from the high resolution, dark matter-only Via Lactea II (VL2) simulation, and find that the number of massive subhalos is dramatically reduced. After adopting a stellar mass to halo mass relationship for the VL2 halos, and identifying subhalos that are (1) likely to be destroyed by stripping and (2) likely to have star formation suppressed by photo-heating, we find that the number of massive, luminous satellites around a Milky Way-mass galaxy is in agreement with the number of observed satellites around the Milky Way or M31. We conclude that baryonic processes have the potential to solve the missing satellites problem.

In order for neutralino dark matter to avoid being overproduced in the early universe, these particles must annihilate (or coannihilate) rather efficiently. Neutralinos with sufficiently large couplings to annihilate at such high a rate (such as those resulting from gaugino-higgsino mixing, as in "well-tempered" or "focus point" scenarios), however, have become increasingly disfavored by the null results of XENON100 and other direct detection experiments. One of the few remaining ways that neutralinos could potentially evade such constraints is if they annihilate through a resonance, as can occur if 2$m_{\chi^0}$ falls within about $\sim$10% of either $m_{A/H}$, $m_h$, or $m_Z$. If no signal is observed from upcoming direct detection experiments, the degree to which such a resonance must be tuned will increase significantly. In this paper, we quantify the degree to which such a resonance must be tuned in order to evade current and projected constraints from direct detection experiments. Assuming a future rate of progress among direct detection experiments that is similar to that obtained over the past decade, we project that within 7 years the light Higgs and $Z$ pole regions will be entirely closed, while the remaining parameter space near the $A/H$ resonance will require that $2m_{\chi^0}$ be matched to the central value (near $m_A$) to within less than 4%. At this rate of progress, it will be a little over a decade before multi-ton direct detection experiments will be able to close the remaining, highly-tuned, regions of the $A/H$ resonance parameter space.

We study the variation of the spectrum of the Fermi Bubbles with Galactic latitude. Far from the Galactic plane (|b| > 30 degrees), the observed gamma-ray emission is nearly invariant with latitude, and is consistent with arising from inverse Compton scattering of the interstellar radiation field by cosmic-ray electrons with an approximately power-law spectrum. The same electrons in the presence of microgauss-scale magnetic fields can also generate the the observed microwave "haze". At lower latitudes (b < 20 degrees), in contrast, the spectrum of the emission correlated with the Bubbles possesses a pronounced spectral feature peaking at 1-4 GeV (in E^2 dN/dE) which cannot be generated by any realistic spectrum of electrons. Instead, we conclude that a second (non-inverse-Compton) emission mechanism must be responsible for the bulk of the low-energy, low-latitude emission. This second component is spectrally similar to the excess GeV emission previously reported from the Galactic Center (GC), and also appears spatially consistent with a luminosity per volume falling approximately as r^-2.4, where r is the distance from the GC. We argue that the spectral feature visible in the low-latitude Bubbles is the extended counterpart of the GC excess, now detected out to at least 2-3 kpc from the GC. The spectrum and angular distribution of the signal is consistent with that predicted from ~10 GeV dark matter particles annihilating to leptons, or from ~50 GeV dark matter particles annihilating to quarks, following a distribution similar to the canonical Navarro-Frenk-White (NFW) profile. We also consider millisecond pulsars as a possible astrophysical explanation for the signal, as observed millisecond pulsars possess a spectral cutoff at approximately the required energy. Any such scenario would require a large population of unresolved millisecond pulsars extending at least 2-3 kpc from the GC.

The spectrum and morphology of gamma-rays from the Galactic Center and the spectrum of synchrotron emission observed from the Milky Way's radio filaments have each been interpreted as possible signals of $\sim$7-10 GeV dark matter particles annihilating in the Inner Galaxy. In dark matter models capable of producing these signals, the annihilations should also generate significant fluxes of $\sim$7-10 GeV positrons which can lead to a distinctive bump-like feature in local cosmic ray positron spectrum. In this letter, we show that while such a feature would be difficult to detect with PAMELA, it would likely be identifiable by the currently operating AMS experiment. As no known astrophysical sources or mechanisms are likely to produce such a sharp feature, the observation of a positron bump at around 7-10 GeV would significantly strengthen the case for a dark matter interpretation of the reported gamma-ray and radio anomalies.

Very recently, it was pointed out that there exists a population of gamma-ray sources without associations at other wavelengths which exhibit spectral features consistent with mono-energetic lines at energies of approximately 111 and 129 GeV. Given recent evidence of similar gamma-ray lines from the Inner Galaxy, it is tempting to interpret these unassociated sources as nearby dark matter subhalos, powered by ongoing annihilations. In this paper, we study the spectrum, luminosity, and angular distribution of these sources, with the intention of testing the hypothesis that they are, in fact, dark matter subhalos. We find that of the 12 sources containing at least one prospective line photon, only 2 exhibit an overall gamma-ray spectrum which is consistent with that predicted from dark matter annihilations (2FGL J2351.6-7558 and 2FGL J0555.9-4348). After discounting the 10 clearly non-dark matter sources, the statistical significance of the remaining two prospective line photons is negligible. That being said, we cannot rule out the possibility that either or both of these sources are dark matter subhalos; their overall luminosity and galactic latitude distribution are not inconsistent with a dark matter origin.

Recent reports of a gamma-ray line feature at ~130 GeV in data from the Fermi Gamma-Ray Space Telescope have generated a great deal of interest in models in which dark matter particles annihilate with a sizable cross section to final states including photons. In this article, we take a model-independent approach, and discuss a number of possibilities for dark matter candidates which could potentially generate such a feature. While we identify several scenarios which could lead to such a gamma-ray line, these models are each fairly constrained. In particular, viable models require large couplings (g>1-3), and additional charged particles with masses in the range of approximately ~130-200 GeV. Furthermore, lower energy gamma-ray constraints from the Galactic Center force us to consider scenarios in which the dark matter annihilates in the early universe through velocity-suppressed processes, or to final states which yield relatively few gamma-rays (such as electrons, muons or neutrinos). An exception to these conclusions can be found in models in which the dark matter annihilates to heavy intermediate states which decay to photons to generate a line-like gamma-ray spectrum.

We consider a simple class of models in which the dark matter, X, is coupled to a new gauge boson, phi, with a relatively low mass (m_phi \sim 100 MeV-3 GeV). Neither the dark matter nor the new gauge boson have tree-level couplings to the Standard Model. The dark matter in this model annihilates to phi pairs, and for a coupling of g_X \sim 0.06 (m_X/10 GeV)^1/2 yields a thermal relic abundance consistent with the cosmological density of dark matter. The phi's produced in such annihilations decay through a small degree of kinetic mixing with the photon to combinations of Standard Model leptons and mesons. For dark matter with a mass of \sim10 GeV, the shape of the resulting gamma-ray spectrum provides a good fit to that observed from the Galactic Center, and can also provide the very hard electron spectrum required to account for the observed synchrotron emission from the Milky Way's radio filaments. For kinetic mixing near the level naively expected from loop-suppressed operators (epsilon \sim 10^{-4}), the dark matter is predicted to scatter elastically with protons with a cross section consistent with that required to accommodate the signals reported by DAMA/LIBRA, CoGeNT and CRESST-II.

Interactions of cosmic rays with the interstellar gas and radiation fields of the Milky Way provide the majority of the gamma rays observed by the Fermi Gamma Ray Space Telescope. In addition to the gas which is densely concentrated along the Galactic Disk, hydrodynamical simulations and observational evidence favor the presence of a halo of hot (T~10^6 K) ionized hydrogen (H_II), extending with non-negligible densities out to the virial radius of the Milky Way. We show that cosmic ray collisions with this circum-galactic gas should be expected to provide a significant flux of gamma rays, on the order of 10% of the observed isotopic gamma ray background at energies above 1 GeV. In addition, gamma rays originating from the extended H_II halos of other galaxies along a given line-of-sight should contribute to this background at a similar level.

Observations by ARCADE-2 and other telescopes sensitive to low frequency radiation have revealed the presence of an isotropic radio background with a hard spectral index. The intensity of this observed background is found to exceed the flux predicted from astrophysical sources by a factor of approximately 5-6. In this article, we consider the possibility that annihilating dark matter particles provide the primary contribution to the observed isotropic radio background through the emission of synchrotron radiation from electron and positron annihilation products. For reasonable estimates of the magnetic fields present in clusters and galaxies, we find that dark matter could potentially account for the observed radio excess, but only if it annihilates mostly to electrons and/or muons, and only if it possesses a mass in the range of approximately 5-50 GeV. For such models, the annihilation cross section required to normalize the synchrotron signal to the observed excess is sigma v ~ (0.4-30) x 10^-26 cm^3/s, similar to the value predicted for a simple thermal relic (sigma v ~ 3 x 10^-26 cm^3/s). We find that in any scenario in which dark matter annihilations are responsible for the observed excess radio emission, a significant fraction of the isotropic gamma ray background observed by Fermi must result from dark matter as well.

In many models, dark matter particles can elastically scatter with nuclei in planets, causing those particles to become gravitationally bound. While the energy expected to be released through the subsequent annihilations of dark matter particles in the interior of the Earth is negligibly small (a few megawatts in the most optimistic models), larger planets that reside in regions with higher densities of slow moving dark matter could plausibly capture and annihilate dark matter at a rate high enough to maintain liquid water on their surfaces, even in the absence of additional energy from starlight or other sources. On these rare planets, it may be dark matter rather than light from a host star that makes it possible for life to emerge, evolve, and survive.

In this article, I summarize and discuss the body of evidence which has accumulated in favor of dark matter in the form of approximately 10 GeV particles. This evidence includes the spectrum and angular distribution of gamma rays from the Galactic Center, the synchrotron emission from the Milky Way's radio filaments, the diffuse synchrotron emission from the Inner Galaxy (the "WMAP Haze") and low-energy signals from the direct detection experiments DAMA/LIBRA, CoGeNT and CRESST-II. This collection of observations can be explained by a relatively light dark matter particle with an annihilation cross section consistent with that predicted for a simple thermal relic (sigma v ~ 10^-26 cm^3/s) and with a distribution in the halo of the Milky Way consistent with that predicted from simulations. Astrophysical explanations for the gamma ray and synchrotron signals, in contrast, have not been successful in accommodating these observations. Similarly, the phase of the annual modulation observed by DAMA/LIBRA (and now supported by CoGeNT) is inconsistent with all known or postulated modulating backgrounds, but are in good agreement with expectations for dark matter scattering. This scenario is consistent with all existing indirect and collider constraints, as well as the constraints placed by CDMS. Consistency with xenon-based experiments can be achieved if the response of liquid xenon to very low-energy nuclear recoils is somewhat suppressed relative to previous evaluations, or if the dark matter possesses different couplings to protons and neutrons.

The South Pole Telescope (SPT), Atacama Cosmology Telescope (ACT), and Wilkinson Microwave Anisotropy Probe (WMAP) have each reported measurements of the cosmic microwave background's (CMB) angular power spectrum which favor the existence of roughly one additional neutrino species, in addition to the three contained in the standard model of particle physics. Neutrinos influence the CMB by contributing to the radiation density, which alters the expansion rate of the universe during the epoch leading up to recombination. In this paper, we consider an alternative possibility that the excess kinetic energy implied by these measurements was possessed by dark matter particles that were produced through a non-thermal mechanism, such as late-time decays. In particular, we find that if a small fraction (<1%) of the dark matter in the universe today were produced through the decays of a heavy and relatively long-lived state, the expansion history of the universe can be indistinguishable from that predicted in the standard cosmological model with an additional neutrino. Furthermore, if these decays take place after the completion of big bang nucleosynthesis, this scenario can avoid tension with the value of three neutrino species preferred by measurements of the light element abundances.

The dark matter halo of the Milky Way is expected to contain an abundance of smaller subhalos. These subhalos can be dense and produce potentially observable fluxes of gamma rays. In this paper, we search for dark matter subhalo candidates among the sources in the Fermi-LAT Second Source Catalog which are not currently identified or associated with counterparts at other wavelengths. Of the nine high-significance, high-latitude (|b|>60 degrees), non-variable, unidentified sources contained in this catalog, only one or two are compatible with the spectrum of a dark matter particle heavier than approximately 50-100 GeV. The majority of these nine sources, however, feature a spectrum that is compatible with that predicted from a lighter (~5-40 GeV) dark matter particle. This population is consistent with the number of observable subhalos predicted for a dark matter candidate in this mass range and with an annihilation cross section of a simple thermal relic (sigma v~3x10^{-26} cm^3/s). Observations in the direction of these sources at other wavelengths will be necessary to either reveal their astrophysical nature (as blazars or other active galactic nuclei, for example), or to further support the possibility that they are dark matter subhalos by failing to detect any non-gamma ray counterpart.

Three dark matter direct detection experiments (DAMA/LIBRA, CoGeNT, and CRESST-II) have each reported signals which are not consistent with known backgrounds, but resemble that predicted for a dark matter particle with a mass of roughly $\sim$10 GeV and an elastic scattering cross section with nucleons of $\sim$$10^{-41}$--$10^{-40}$ cm$^2$. In this article, we compare the signals of these experiments and discuss whether they can be explained by a single species of dark matter particle, without conflicting with the constraints of other experiments. We find that the spectrum of events reported by CoGeNT and CRESST-II are consistent with each other and with the constraints from CDMS-II, although some tension with xenon-based experiments remains. Similarly, the modulation signals reported by DAMA/LIBRA and CoGeNT appear to be compatible, although the corresponding amplitude of the observed modulations are a factor of at least a few higher than would be naively expected, based on the event spectra reported by CoGeNT and CRESST-II. This apparent discrepancy could potentially be resolved if tidal streams or other non-Maxwellian structures are present in the local distribution of dark matter.

New gauge bosons with Standard Model-like couplings to leptons are constrained by collider searches to be heavier than approximately ~1 TeV. A Z' boson with suppressed couplings to leptons, however, could be much lighter and possess substantial couplings to Standard Model quarks. In this article, we consider a new leptophobic Z' gauge boson as a simple and well motivated extension of the Standard Model, and discuss several of its possible signatures at the Tevatron. We find that three of the recent anomalies reported from the Tevatron - in particular the top-quark forward-backward asymmetry and excesses in the 3b and W + 2 jets final states - could be explained by a new Z' with a mass of approximately 150 GeV, relatively large couplings to quarks, and suppressed couplings to electrons and muons. Moreover, we find that such a particle could also mediate the interactions of dark matter, leading to potentially interesting implications for direct detection experiments.

Recent observations by the CoGeNT collaboration (as well as long standing observations by DAMA/LIBRA) suggest the presence of a $\sim 5$-10 GeV dark matter particle with a somewhat large elastic scattering cross section with nucleons ($\sigma\sim 10^{-40}$ cm$^2$).Within the context of the minimal supersymmetric standard model (MSSM), neutralinos in this mass range are not able to possess such large cross sections, and would be overproduced in the early universe. Simple extensions of the MSSM, however, can easily accommodate these observations. In particular, the extension of the MSSM by a chiral singlet superfield allows for the possibility that the dark matter is made up of a light singlino that interacts with nucleons largely through the exchange of a fairly light ($\sim$30-70 GeV) singlet-like scalar higgs, $\hi$. Such a scenario is consistent with all current collider constraints and can generate the signals reported by CoGeNT and DAMA/LIBRA. Furthermore, there is a generic limit of the extended model in which there is a singlet-like pseudoscalar higgs, $\ai$, with $\mai\sim \mhi$ and in which the $\chi^0\chi^0$ and $b\anti b, s\anti s$ coupling magnitudes of the $\hi$ and $\ai$ are very similar. In this case, the thermal relic abundance is automatically consistent with the measured density of dark matter if $\mchi$ is sufficiently small that $\chi^0\chi^0\to b\anti b$ is forbidden.

The region surrounding the center of the Milky Way is both astrophysically rich and complex, and is predicted to contain very high densities of dark matter. Utilizing three years of data from the Fermi Gamma Ray Space Telescope (and the recently available Pass 7 ultraclean event class), we study the morphology and spectrum of the gamma ray emission from this region and find evidence of a spatially extended component which peaks at energies between 300 MeV and 10 GeV. We compare our results to those reported by other groups and find good agreement. The extended emission could potentially originate from either the annihilations of dark matter particles in the inner galaxy, or from the collisions of high energy protons that are accelerated by the Milky Way's supermassive black hole with gas. If interpreted as dark matter annihilation products, the emission spectrum favors dark matter particles with a mass in the range of 7-12 GeV (if annihilating dominantly to leptons) or 25-45 GeV (if annihilating dominantly to hadronic final states). The intensity of the emission corresponds to a dark matter annihilation cross section consistent with that required to generate the observed cosmological abundance in the early universe (sigma v ~ 3 x 10^-26 cm^3/s). We also present conservative limits on the dark matter annihilation cross section which are at least as stringent as those derived from other observations.

We explore the phenomenology of Grand Unified Models based on the E_6 group, focusing on the Z' with suppressed couplings to leptons that can appear in such models. We find that this Z' can accommodate the W+dijets anomaly reported by the CDF collaboration. Furthermore, a viable dark matter candidate in the form of a right-handed sneutrino is also present within the fundamental 27-dimensional representation of E_6. Through its sizable couplings to the Z', the dark matter is predicted to possess an elastic scattering cross section with neutrons which can generate the signals reported by the CoGeNT and DAMA/LIBRA collaborations. To avoid being overproduced in the early universe, the dark matter must annihilate to leptons through the exchange of charged or neutral fermions which appear in the 27 of E_6, providing an excellent fit to the gamma ray spectrum observed from the Galactic Center by the Fermi Gamma Ray Space Telescope.

The inner degrees of the Galactic center contain a large population of filamentary structures observed at radio frequencies. These so-called non-thermal radio filaments (NRFs) trace magnetic field lines and have attracted significant interest due to their hard (S_v ~ -0.1 +/- 0.4) synchrotron emission spectra. The origin of these filaments remains poorly understood. We show that the electrons and positrons created through the annihilations of a relatively light (~5-10 GeV) dark matter particle with the cross section predicted for a simple thermal relic can provide a compelling match to the intensity, spectral shape, and flux variation of the NRFs. Furthermore, the characteristics of the dark matter particle necessary to explain the synchrotron emission from the NRFs is consistent with those required to explain the excess gamma-ray emission observed from the Galactic center by the Fermi-LAT, as well as the direct detection signals observed by CoGeNT and DAMA/LIBRA.

The CoGeNT collaboration has recently presented the results of their first 15 months of data, including the measurement of the spectrum of nuclear recoil candidate events, and the time variation of those events. These results appear consistent with the signal anticipated from a relatively light dark matter particle scattering elastically with nuclei. In this paper, we independently analyze the data set collected by CoGeNT and explore the implications of these results for dark matter. We find that the observed spectrum and rate is consistent with originating from dark matter particles with a mass in the range of 4.5-12 GeV and an elastic scattering cross section with nucleons of approximately ~10^-40 cm^2. We confirm the conclusion of the CoGeNT collaboration that the data also includes a somewhat statistically significant (2.7 sigma) indication of annual modulation, with a phase, period, and amplitude consistent with that predicted for dark matter. CoGeNT's phase is also consistent with the annual modulation reported by the DAMA/LIBRA collaboration. We also discuss the null results reported by CDMS and XENON100, and comment on the prospects for other experiments to detect a dark matter particle with the properties implied by CoGeNT.

A simple explanation of the W+dijet excess recently reported by the CDF collaboration involves the introduction of a new gauge boson with sizable couplings to quarks, but with no or highly suppressed couplings to leptons. Anomaly-free theories which include such a leptophobic gauge boson must also include additional particle content, which may include a stable and otherwise viable candidate for dark matter. Based on the couplings and mass of the Z' required to generate the CDF excess, we predict such a dark matter candidate to possess an elastic scattering cross section with nucleons on the order of sigma ~ 10^-40 cm^2, providing a natural explanation for the signals reported by the CoGeNT and DAMA/LIBRA collaborations. In this light, CDF may be observing the gauge boson responsible for the force which mediates the interactions between the dark and visible matter of our universe.

As very high energy (~100 GeV) gamma rays travel over cosmological distances, their flux is attenuated through interactions with the extragalactic background light. Observations of distant gamma ray sources at energies between ~200 GeV and a few TeV by ground-based gamma ray telescopes such as HESS, however, suggest that the universe is more transparent to very high energy photons than had been anticipated. One possible explanation for this is the existence of axion-like-particles (ALPs) which gamma rays can efficiently oscillate into, enabling them to travel cosmological distances without attenuation. In this article, we use data from the Fermi Gamma Ray Space Telescope to calculate the spectra at 1-100 GeV of two gamma ray sources, 1ES1101-232 at redshift z=0.186 and H2356-309 at z=0.165, and use this in conjunction with the measurements of ground-based telescopes to test the ALP hypothesis. We find that the observations can be well-fit by an intrinsic power-law source spectrum with indices of -1.72 and -2.1 for 1ES1101-232 and H2356-309, respectively, and that no ALPs or other exotic physics is necessary to explain the observed degree of attenuation.

Recently, an analysis of data from the Fermi Gamma Ray Space Telescope has revealed a flux of gamma rays concentrated around the inner ~0.5 degrees of the Milky Way, with a spectrum that is sharply peaked at 2-4 GeV. If interpreted as the products of annihilating dark matter, this signal implies that the dark matter consists of particles with a mass between 7.3 and 9.2 GeV annihilating primarily to charged leptons. This mass range is very similar to that required to accommodate the signals reported by CoGeNT and DAMA/LIBRA. In addition to gamma rays, the dark matter is predicted to produce energetic electrons and positrons in the Inner Galaxy, which emit synchrotron photons as a result of their interaction with the galactic magnetic field. In this letter, we calculate the flux and spectrum of this synchrotron emission assuming that the gamma rays from the Galactic Center originate from dark matter, and compare the results to measurements from the WMAP satellite. We find that a sizable flux of hard synchrotron emission is predicted in this scenario, and that this can easily account for the observed intensity, spectrum, and morphology of the "WMAP Haze".

We analyze the first two years of data from the Fermi Gamma Ray Space Telescope from the direction of the inner 10 degrees around the Galactic Center with the intention of constraining, or finding evidence of, annihilating dark matter. We find that the morphology and spectrum of the emission between 1.25 degrees and 10 degrees from the Galactic Center is well described by a the processes of decaying pions produced in cosmic ray collisions with gas, and the inverse Compton scattering of cosmic ray electrons in both the disk and bulge of the Inner Galaxy, along with gamma rays from known points sources in the region. The observed spectrum and morphology of the emission within approximately 1.25 degrees (~175 parsecs) of the Galactic Center, in contrast, departs from the expectations for by these processes. Instead, we find an additional component of gamma ray emission that is highly concentrated around the Galactic Center. The observed morphology of this component is consistent with that predicted from annihilating dark matter with a cusped (and possibly adiabatically contracted) halo distribution (density proportional to r^{-gamma}, with gamma=1.18 to 1.33). The observed spectrum of this component, which peaks at energies between 1-4 GeV (in E^2 units), can be well fit by a 7-10 GeV dark matter particle annihilating primarily to tau leptons with a cross section in the range of 4.6 x 10^-27 to 5.3 x 10^-26 cm^3/s, depending on how the dark matter distribution is normalized. We also discuss other sources for this emission, including the possibility that much of it originates from the Milky Way's supermassive black hole.

It has previously been shown that the excess of events reported by the CoGeNT collaboration could be generated by elastically scattering dark matter particles with a mass of approximately 5-15 GeV. This mass range is very similar to that required to generate the annual modulation observed by DAMA/LIBRA and the gamma rays from the region surrounding the Galactic Center identified within the data of the Fermi Gamma Ray Space Telescope. To confidently conclude that CoGeNT's excess is the result of dark matter, however, further data will likely be needed. In this paper, we make projections for the first full year of CoGeNT data, and for its planned upgrade. Not only will this body of data more accurately constrain the spectrum of nuclear recoil events, and corresponding dark matter parameter space, but will also make it possible to identify seasonal variations in the rate. In particular, if the CoGeNT excess is the product of dark matter, then one year of CoGeNT data will likely reveal an annual modulation with a significance of 2-3$\sigma$. The planned CoGeNT upgrade will not only detect such an annual modulation with high significance, but will be capable of measuring the energy spectrum of the modulation amplitude. These measurements will be essential to irrefutably confirming a dark matter origin of these events.

Recent results from the CoGeNT collaboration (as well as the annual modulation reported by DAMA/LIBRA) point toward dark matter with a light (5-10 GeV) mass and a relatively large elastic scattering cross section with nucleons (\sigma ~ 10^{-40} cm^2). In order to possess this cross section, the dark matter must communicate with the Standard Model through mediating particles with small masses and/or large couplings. In this Letter, we explore with a model independent approach the particle physics scenarios that could potentially accommodate these signals. We also discuss how such models could produce the gamma rays from the Galactic Center observed in the data of the Fermi Gamma Ray Space Telescope. We find multiple particle physics scenarios in which each of these signals can be accounted for, and in which the dark matter can be produced thermally in the early Universe with an abundance equal to the measured cosmological density.

In this paper, we study the recent excess of low energy events observed by the CoGeNT collaboration and the annual modulation reported by the DAMA/LIBRA collaboration, and discuss whether these signals could both be the result of the same elastically scattering dark matter particle. We find that, without channeling but when taking into account uncertainties in the relevant quenching factors, a dark matter candidate with a mass of approximately ~7.0 GeV and a cross section with nucleons of sigma_{DM-N} ~2x10^-4 pb (2x10^-40 cm^2) could account for both of these observations. We also comment on the events recently observed in the oxygen band of the CRESST experiment and point out that these could potentially be explained by such a particle. Lastly, we compare the region of parameter space favored by DAMA/LIBRA and CoGeNT to the constraints from XENON 10, XENON 100, and CDMS (Si) and find that these experiments cannot at this time rule out a dark matter interpretation of these signals.

Future galactic supernovae will provide an extremely long baseline for studying the properties and interactions of neutrinos. In this paper, we discuss the possibility of using such an event to constrain (or discover) the effects of exotic physics in scenarios that are not currently constrained and are not accessible with reactor or solar neutrino experiments. In particular, we focus on the cases of neutrino decay and quantum decoherence. We calculate the expected signal from a core-collapse supernova in both current and future water Cerenkov, scintillating, and liquid argon detectors, and find that such observations will be capable of distinguishing between many of these scenarios. Additionally, future detectors will be capable of making strong, model-independent conclusions by examining events associated with a galactic supernova's neutronization burst.

Although recent measurements of the shower profiles of ultra-high energy cosmic rays suggest that they are largely initiated by heavy nuclei, such conclusions rely on hadronic interaction models which have large uncertainties. We investigate an alternative test of cosmic ray composition which is based on the observation of ultra-high energy photons produced through cosmic ray interactions with diffuse low energy photon backgrounds during intergalactic propagation. We show that if the ultra-high energy cosmic rays are dominated by heavy nuclei, the flux of these photons is suppressed by approximately an order of magnitude relative to the proton-dominated case. Future observations by the Pierre Auger Observatory may be able to use this observable to constrain the composition of the primaries, thus providing an important cross-check of hadronic interaction models.

The Alpha Magnetic Spectrometer (AMS-02), which is scheduled to be deployed onboard the International Space Station later this year, will be capable of measuring the composition and spectra of GeV-TeV cosmic rays with unprecedented precision. In this paper, we study how the projected measurements from AMS-02 of stable secondary-to-primary and unstable ratios (such as boron-to-carbon and beryllium-10-to-beryllium-9) can constrain the models used to describe the propagation of cosmic rays throughout the Milky Way. We find that within the context of fairly simple propagation models, all of the model parameters can be determined with high precision from the projected AMS-02 data. Such measurements are less constraining in more complex scenarios, however, which allow for departures from a power-law form for the diffusion coefficient, for example, or for inhomogeneity or stochasticity in the distribution and chemical abundances of cosmic ray sources.

The Milky Way's dark matter halo is thought to contain large numbers of smaller subhalos. These objects can contain very high densities of dark matter, and produce potentially observable fluxes of gamma rays. In this article, we study the gamma ray sources in the Fermi Gamma Ray Space Telescope's recently published First Source Catalog, and attempt to determine whether this catalog might contain a population of dark matter subhalos. We find that, while approximately 20-60 of the catalog's unidentified sources could plausibly be dark matter subhalos, such a population cannot be clearly identified as such at this time. From the properties of the sources in the First Source Catalog, we derive limits on the dark matter's annihilation cross section that are comparably stringent to those derived from recent observations of dwarf spheroidal galaxies.

In this paper, we compare the relative sensitivities of gamma-ray and neutrino observations to the dark matter annihilation cross section in leptophilic models such as have been designed to explain PAMELA data. We investigate whether the high energy neutrino telescope IceCube will be competitive with current and upcoming searches by gamma-ray telescopes, such as the Atmospheric Cerenkov Telescopes (ACTs) (HESS, VERITAS and MAGIC), or the Fermi Gamma Ray Space Telescope, in detecting or constraining dark matter particles annihilating in dwarf spheroidal galaxies. We find that after ten years of observation of the most promising nearby dwarfs, IceCube will have sensitivity comparable to the current sensitivity of gamma-ray telescopes only for very heavy (m_X > 7 TeV) or relatively light (m_X < 200 GeV) dark matter particles which annihilate primarily to mu+mu-. If dark matter particles annihilate primarily to tau+tau-, IceCube will have superior sensitivity only for dark matter particle masses below the 200 GeV threshold of current ACTs. If dark matter annihilations proceed directly to neutrino-antineutrino pairs a substantial fraction of the time, IceCube will be competitive with gamma-ray telescopes for a much wider range of dark matter masses.

The Pierre Auger Observatory's (PAO) shower profile measurements can be used to constrain the chemical composition of the ultra-high energy cosmic ray (UHECR) spectrum. In particular, the PAO's measurements of the average depth of shower maximum and the fluctuations of the depth of shower maximum indicate that the cosmic ray spectrum is dominated by a fairly narrow distribution (in charge) of heavy or intermediate mass nuclei at the highest measured energies (E > 10^{19} eV), and contains mostly lighter nuclei or protons at lower energies (E ~ 10^{18} eV). In this article, we study the propagation of UHECR nuclei with the goal of using these measurements, along with those of the shape of the spectrum, to constrain the chemical composition of the particles accelerated by the sources of the UHECRs. We find that with modest intergalactic magnetic fields, 0.3 nG in strength with 1 Mpc coherent lengths, good fits to the combined PAO data can be found for the case in which the sources accelerate primarily intermediate mass nuclei (such as nitrogen or silicon). Without intergalactic magnetic fields, we do not find any composition scenarios that can accommodate the PAO data. For a spectrum dominated by heavy or intermediate mass nuclei, the Galactic (and intergalactic) magnetic fields are expected to erase any significant angular correlation between the sources and arrival directions of UHECRs.

In addition to gamma-rays, dark matter annihilation products can include energetic electrons which inverse Compton scatter with the cosmic microwave background to produce a diffuse extragalactic background of gamma-rays and X-rays. In models in which the dark matter particles annihilate primarily to electrons or muons, the measurements of EGRET and COMPTEL can provide significant constraints on the annihilation cross section. The Fermi Gamma-Ray Space Telescope will likely provide an even more stringent test of such scenarios.

The rising cosmic ray positron fraction reported by the PAMELA collaboration has lead to a great deal of interest in astrophysical sources of energetic electrons and positrons, including pulsars. In this paper, we calculate the spectrum of synchrotron emission from electrons and positrons injected from 376 young pulsars (<10^6 years) contained in the ATNF catalog, and compare our results to observations. We find that if objects such as the Vela and Crab pulsars have injected ~10^48 erg or more in energetic electrons and/or positrons, they are expected to produce bright and distinctive features in the synchrotron sky. Intriguingly, we predict hard synchrotron emission from these regions of the sky which is qualitatively similar to that observed by WMAP.

Annihilation of cosmologically distributed dark matter is predicted to produce a potentially observable flux of high energy photons. Neglecting the contribution from local structure, this signal is predicted to be virtually uniform on the sky and, in order to be identified, must compete with various extragalactic backgrounds. We focus here on unresolved blazars and discuss several techniques for separating the dark matter signal from this background. First, the spectral shapes of the signal and background are expected to differ, a feature which can be exploited with the Fisher Matrix formalism. Second, in any given angular pixel, the number of photons from blazars is drawn from a distribution which is far from Poisson. A knowledge of this distribution enhances one's ability to extract the dark matter signal, while ignorance of it can lead to the introduction of a large systematic error.

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.

Despite compelling arguments that significant discoveries of physics beyond the standard model are likely to be made at the Large Hadron Collider, it remains possible that this machine will make no such discoveries, or will make no discoveries directly relevant to the dark matter problem. In this article, we study the ability of astrophysical experiments to deduce the nature of dark matter in such a scenario. In most dark matter studies, the relic abundance and detection prospects are evaluated within the context of some specific particle physics model or models (e.g. supersymmetry). Here, we attempt to develop a model-independent approach toward the phenomenology of weakly interacting massive particles in the absence of any discoveries at the Large Hadron Collider. In particular, we consider generic fermionic or scalar dark matter particles with a variety of interaction forms, and calculate the corresponding constraints from and sensitivity of direct and indirect detection experiments. The results may provide some guidance in disentangling information from future direct and indirect detection experiments.

We revisit the prospects for IceCube and similar kilometer-scale telescopes to detect neutrinos produced by the annihilation of weakly interacting massive dark matter particles (WIMPs) in the Sun. We emphasize that the astrophysics of the problem is understood; models can be observed or, alternatively, ruled out. In searching for a WIMP with spin-independent interactions with ordinary matter, IceCube is only competitive with direct detection experiments if the WIMP mass is sufficiently large. For spin-dependent interactions IceCube already has improved the best limits on spin-dependent WIMP cross sections by two orders of magnitude. This is largely due to the fact that models with significant spin-dependent couplings to protons are the least constrained and, at the same time, the most promising because of the efficient capture of WIMPs in the Sun. We identify models where dark matter particles are beyond the reach of any planned direct detection experiments while being within reach of neutrino telescopes. In summary, we find that, even when contemplating recent direct detection results, neutrino telescopes have the opportunity to play an important as well as complementary role in the search for particle dark matter.

PAMELA's observation that the cosmic ray positron fraction increases rapidly with energy implies the presence of primary sources of energetic electron-positron pairs. Of particular interest is the possibility that dark matter annihilations in the halo of the Milky Way provide this anomalous flux of antimatter. The recent measurement of the cosmic ray electron spectrum by the Fermi Gamma Ray Space Telescope, however, can be used to constrain the nature of any such dark matter particle. In particular, it has been argued that in order to accommodate the observations of Fermi and provide the PAMELA positron excess, annihilating dark matter particles must be as massive as ~1 TeV or heavier. In this article, we revisit Fermi's electron spectrum measurement within the context of annihilating dark matter, focusing on masses in the range of 100-1000 GeV, and considering effects such as variations in the astrophysical backgrounds from the presence of local cosmic ray accelerators, and the finite energy resolution of the Fermi Gamma Ray Space Telescope. When these factors are taken into account, we find that dark matter particles as light as ~300 GeV can be capable of generating the positron fraction observed by PAMELA.

In models with one universal extra dimension (UED), the first Kaluza-Klein excitations of the hypercharge gauge boson, B^(1), and the neutral component of isospin gauge boson, W^3(1), are each viable dark matter candidates. In either case, such particles are predicted to accumulate in the core of the Sun, where they annihilate to generate a potentially observable flux of high energy neutrinos. In this article, we calculate the flux of neutrinos produced in this model and determine the constraints that can be placed on the UED parameter space from current IceCube data. For the case of B^(1) dark matter, we find that the present limits from IceCube are stronger than those from direct dark matter detection experiments such as CDMS and XENON10. For W^3(1) dark matter, the present IceCube data provides a constraint slightly weaker than direct detection experiments. In addition, we also present the projected regions of UED parameter space that can be probed by IceCube/DeepCore in the near future and compare them to the prospects for future direct detection experiments.

Decaying dark matter has previously been proposed as a possible explanation for the excess high energy cosmic ray electrons and positrons seen by PAMELA and the Fermi Gamma-Ray Space Telescope (FGST). To accommodate these signals however, the decays must be predominantly leptonic, to muons or taus, and therefore produce neutrinos, potentially detectable with the IceCube neutrino observatory. We find that, with five years of data, IceCube (supplemented by DeepCore) will be able to significantly constrain the relevant parameter space of decaying dark matter, and may even be capable of discovering dark matter decaying in the halo of the Milky Way.

We consider a scenario within the Minimal Supersymmetric Standard Model extended by a singlet chiral superfield, in which neutralino dark matter annihilates to light singlet-like Higgs bosons, which proceed to decay to either electron-positron or muon-antimuon pairs. Unlike neutralino annihilations in the MSSM, this model can provide a good fit to the PAMELA cosmic ray positron fraction excess. Furthermore, the singlet-like scalar Higgs can induce a large Sommerfeld enhancement and provide an annihilation rate sufficient to accommodate the observed positron excess.

Recently, observations by PAMELA, the Fermi Gamma Ray Space Telescope, and other cosmic ray experiments have generated a great deal of interest in dark matter (DM) particles which annihilate at a high rate to leptons. In this letter, we explore the possibility of using large volume neutrino telescopes, such as IceCube, to constrain such models; specifically we consider signals due to DM annihilation in the inner Milky Way. We find that, if Dark Matter annihilations are responsible for the signals observed by PAMELA and FGST, then IceCube (in conjunction with the planned low threshold extension, DeepCore) should detect or exclude the corresponding neutrino signal from the inner Milky Way with a few years of observation.

In recent years, a number of experiments have been conducted with the goal of studying cosmic rays at GeV to TeV energies. This is a particularly interesting regime from the perspective of indirect dark matter detection. To draw reliable conclusions regarding dark matter from cosmic ray measurements, however, it is important to first understand the propagation of cosmic rays through the magnetic and radiation fields of the Milky Way. In this paper, we constrain the characteristics of the cosmic ray propagation model through comparison with observational inputs, including recent data from the CREAM experiment, and use these constraints to estimate the corresponding uncertainties in the spectrum of cosmic ray electrons and positrons from dark matter particles annihilating in the halo of the Milky Way.

Although empirical evidence indicates that that the universe's gas had become ionized by redshift z ~ 6, the mechanism by which this transition occurred remains unclear. In this article, we explore the possibility that dark matter annihilations may have played the dominant role in this process. Energetic electrons produced in these annihilations can scatter with the cosmic microwave background to generate relatively low energy gamma rays, which ionize and heat gas far more efficiently than higher energy prompt photons. In contrast to previous studies, we find that viable dark matter candidates with electroweak scale masses can naturally provide the dominant contribution to the reionization of the universe. Intriguingly, we find that dark matter candidates capable of producing the recent cosmic ray excesses observed by PAMELA and/or ATIC are also predicted to lead to the full reionization of the universe by z ~ 6.

In this letter, we study the possibility that Kaluza-Klein dark matter in a model with one universal extra dimension is responsible for the recent observations of the PAMELA and ATIC experiments. In this model, the dark matter particles annihilate largely to charged leptons, which enables them to produce a spectrum of cosmic ray electrons and positrons consistent with the PAMELA and ATIC measurements. To normalize to the observed signal, however, large boost factors (~10^3) are required. Despite these large boost factors and significant annihilation to hadronic modes (35%), we find that the constraints from cosmic ray antiproton measurements can be satisfied. Relic abundance considerations in this model force us to consider a rather specific range of masses (approximately 600-900 GeV) which is very similar to the range required to generate the ATIC spectral feature. The results presented here can also be used as a benchmark for model-independent constraints on dark matter annihilation to hadronic modes.

Recent results from the PAMELA satellite indicate the presence of a large flux of positrons (relative to electrons) in the cosmic ray spectrum between approximately 10 and 100 GeV. As annihilating dark matter particles in many models are predicted to contribute to the cosmic ray positron spectrum in this energy range, a great deal of interest has resulted from this observation. Here, we consider pulsars (rapidly spinning, magnetized neutron stars) as an alternative source of this signal. After calculating the contribution to the cosmic ray positron and electron spectra from pulsars, we find that the spectrum observed by PAMELA could plausibly originate from such sources. In particular, a significant contribution is expected from the sum of all mature pulsars throughout the Milky Way, as well as from the most nearby mature pulsars (such as Geminga and B0656+14). The signal from nearby pulsars is expected to generate a small but significant dipole anisotropy in the cosmic ray electron spectrum, potentially providing a method by which the Fermi gamma-ray space telescope would be capable of discriminating between the pulsar and dark matter origins of the observed high energy positrons.

In this letter, we suggest that a nearby clump of 600-1000 GeV neutralinos may be responsible for the excesses recently observed in the cosmic ray positron and electron spectra by the PAMELA and ATIC experiments. Although neutralino dark matter annihilating throughout the halo of the Milky Way is predicted to produce a softer spectrum than is observed, and violate constraints from cosmic ray antiproton measurements, a large nearby (within 1-2 kiloparsecs of the Solar System) clump of annihilating neutralinos can lead to a spectrum which is consistent with PAMELA and ATIC, while also producing an acceptable antiproton flux. Furthermore, the presence of a large dark matter clump can potentially accommodate the very large annihilation rate required to produce the PAMELA and ATIC signals.

This is a report on the findings of the dark matter science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses the prospects for dark matter detection with future gamma-ray experiments, and the complementarity of gamma-ray measurements with other indirect, direct or accelerator-based searches. We conclude that any comprehensive search for dark matter should include gamma-ray observations, both to identify the dark matter particle (through the charac- teristics of the gamma-ray spectrum) and to measure the distribution of dark matter in galactic halos.

Recent results from the Advanced Thin Ionization Calorimeter (ATIC) balloon experiment have identified the presence of a spectral feature between approximately 300 and 800 GeV in the cosmic ray electron spectrum. This spectral feature appears to imply the existence of a local (within about 1 kpc) source of high energy electrons. Emission from a local pulsar and dark matter annihilations have each been put forth as possible origins of this signal. In this letter, we consider the sensitivity of ground based atmospheric Cherenkov telescopes to electrons and show that observatories such as HESS and VERITAS should be able to resolve this feature with sufficient precision to discriminate between the dark matter and pulsar hypotheses with considerably greater than 5 sigma significance, even for conservative assumptions regarding their performance. In addition, this feature provides an opportunity to perform an absolute calibration of the energy scale of ground based, gamma ray telescopes.

It is likely that ultra-high energy cosmic rays contain a significant component of heavy or intermediate mass nuclei. The propagation of ultra-high energy nuclei through cosmic radiation backgrounds is more complicated than that of protons and its study has required the use of Monte Carlo techniques. We present an analytic method for calculating the spectrum and the composition at Earth of ultra-high energy cosmic rays which start out as heavy nuclei from their extragalactic sources. The results obtained are in good agreement with those obtained using numerical methods.

In most global fits of the constrained minimal supersymmetric model (CMSSM) to indirect data, the a priori likelihoods of any two points in tan beta are treated as equal, and the more fundamental mu and B Higgs potential parameters are fixed by potential minimization conditions. We find that, if instead a flat ("natural")prior measure on mu and B is placed, a strong preference exists for the focus point region from fits to particle physics and cosmological data. In particular, we find that the lightest neutralino is strongly favored to be a mixed bino-higgsino (~10% higgsino). Such mixed neutralinos have large elastic scattering cross sections with nuclei, leading to extremely promising prospects for both underground direct detection experiments and neutrino telescopes. In particular, the majority of the posterior probability distribution falls within parameter space within an order of magnitude of current direct detection constraints. Furthermore, neutralino annihilations in the sun are predicted to generate thousands of neutrino induced muon events per years at IceCube. Thus, assuming the framework of the CMSSM and using the natural prior measure, modulo caveats regarding astrophysical uncertainties, we are likely to be living in a world with good prospects for the direct and indirect detection of neutralino dark matter.

Motivated by Higgs Portal and Hidden Valley models, heavy particle dark matter that communicates with the supersymmetric Standard Model via pure Higgs sector interactions is considered. We show that a thermal relic abundance consistent with the measured density of dark matter is possible for masses up to $\sim 30\tev$. For dark matter masses above $\sim 1\tev$, non-perturbative Sommerfeld corrections to the annihilation rate are large, and have the potential to greatly affect indirect detection signals. For large dark matter masses, the Higgs-dark-matter-sector couplings are large and we show how such models may be given a UV completion within the context of so-called "Fat-Higgs" models. Higgs Portal dark matter provides an example of an attractive alternative to conventional MSSM neutralino dark matter that may evade discovery at the LHC, while still being within the reach of current and upcoming indirect detection experiments.

The GLAST satellite mission will study the gamma ray sky with considerably greater exposure than its predecessor EGRET. In addition, it will be capable of measuring the arrival directions of gamma rays with much greater precision. These features each significantly enhance GLAST's potential for identifying gamma rays produced in the annihilations of dark matter particles. The combined use of spectral and angular information, however, is essential if the full sensitivity of GLAST to dark matter is to be exploited. In this paper, we discuss the separation of dark matter annihilation products from astrophysical backgrounds, focusing on the Galactic Center region, and perform a forecast for such an analysis. We consider both point-like and diffuse astrophysical backgrounds and model them using a point-spread-function for GLAST. While the results of our study depend on the specific characteristics of the dark matter signal and astrophysical backgrounds, we find that in many scenarios it is possible to successfully identify dark matter annihilation radiation, even in the presence of significant astrophysical backgrounds.

Very high energy gamma-rays are expected to be absorbed by the extragalactic background light over cosmological distances via the process of electron-positron pair production. Recent observations of cosmologically distant gamma-ray emitters by ground based gamma-ray telescopes have, however, revealed a surprising degree of transparency of the universe to very high energy photons. One possible mechanism to explain this observation is the oscillation between photons and axion-like-particles (ALPs). Here we explore this possibility further, focusing on photon-ALP conversion in the magnetic fields in and around gamma-ray sources and in the magnetic field of the Milky Way, where some fraction of the ALP flux is converted back into photons. We show that this mechanism can be efficient in allowed regions of the ALP parameter space, as well as in typical configurations of the Galactic Magnetic Field. As case examples, we consider the spectrum observed from two HESS sources: 1ES1101-232 at redshift z=0.186 and H 2356-309 at z=0.165. We also discuss features of this scenario which could be used to distinguish it from standard or other exotic models.

In this review, we discuss the role of the various experimental programs taking part in the broader effort to identify the particle nature of dark matter. In particular, we focus on electroweak scale dark matter particles and discuss a wide range of search strategies being carried out and developed to detect them. These efforts include direct detection experiments, which attempt to observe the elastic scattering of dark matter particles with nuclei, indirect detection experiments, which search for photons, antimatter and neutrinos produced as a result of dark matter annihilations, and collider searches for new TeV-scale physics. Each of these techniques could potentially provide a different and complementary set of information related to the mass, interactions and distribution of dark matter. Ultimately, it is hoped that these many different tools will be used together to conclusively identify the particle or particles that constitute the dark matter of our universe.

The annihilations of neutralino dark matter (or other dark matter candidate) generate, among other Standard Model states, electrons and positrons. These particles emit synchrotron photons as a result of their interaction with the Galactic Magnetic Field. In this letter, we use the measurements of the WMAP satellite to constrain the intensity of this synchrotron emission and, in turn, the annihilation cross section of the lightest neutralino. We find this constraint to be more stringent than that provided by any other current indirect detection channel. In particular, the neutralino annihilation cross section must be less than ~ 3 x 10^-26 cm^3/s (1 x 10^25 cm^3/s) for 100 GeV (500 GeV) neutralinos distributed with an NFW halo profile. For the conservative case of an entirely flat dark matter distribution within the inner 8 kiloparsecs of the Milky Way, the constraint is approximately a factor of 30 less stringent. Even in this conservative case, synchrotron measurements strongly constrain, for example, the possibility of wino or higgsino neutralino dark matter produced non-thermally in the early universe.

Searches for supersymmetry are among the most exciting physics goals at Run II of the Tevatron. In particular, in supersymmetric models with light charginos, neutralinos and sleptons, associated chargino--neutralino production can potentially be observed as multi-lepton events with missing energy. We discuss how, in the generic TeV-scale MSSM, the prospects for these chargino-neutralino searches are impacted by cosmological considerations, namely the neutralino relic abundance and direct detection limits. We also discuss what an observation of chargino-neutralino production at the Tevatron would imply for the prospects of future direct dark matter searches without assuming specific patterns of supersymmetry breaking.

We propose that axion-like particles (ALPs) with a two-photon vertex, consistent with all astrophysical and laboratory bounds, may lead to a detectable signature in the spectra of high-energy gamma ray sources. This occurs as a result of gamma rays being converted into ALPs in the magnetic fields of efficient astrophysical accelerators according to the "Hillas criterion", such as jets of active galactic nuclei or hot spots of radio galaxies. The discovery of such an effect is possible by GLAST in the 1-100 GeV range and by ground based gamma ray telescopes in the TeV range.

Ongoing experimental efforts to detect cosmic sources of high energy neutrinos are guided by the expectation that astrophysical accelerators of cosmic ray protons would also generate neutrinos through interactions with ambient matter and/or photons. However there will be a reduction in the predicted neutrino flux if cosmic ray sources accelerate not only protons but also significant number of heavier nuclei, as is indicated by recent air shower data. We consider plausible extragalactic sources such as active galactic nuclei, gamma-ray bursts and starburst galaxies and demand consistency with the observed cosmic ray composition and energy spectrum at Earth after allowing for propagation through intergalactic radiation fields. This allows us to calculate the expected neutrino fluxes from the sources, normalised to the observed cosmic ray spectrum. We find that the likely signals are still within reach of next generation neutrino telescopes such as IceCube.

For the SUSY 2007 conference, I was asked to review the topic of indirect searches for dark matter. As part of that talk, I summarized several observations which have been interpreted as the product of dark matter annihilations. In my contribution to the proceedings, I have decided to focus on this aspect of my talk. In particular, I will discuss the cosmic positron spectrum measured by HEAT, the 511 keV emission from the Galactic Bulge measured by INTEGRAL, the diffuse galactic and extragalactic gamma ray spectra measured by EGRET, and the microwave excess from the Galactic Center observed by WMAP.

We explore the properties of dark matter in theories with two universal extra dimensions, where the lightest Kaluza-Klein state is a spin-0 neutral particle, representing a six-dimensional photon polarized along the extra dimensions. Annihilation of this 'spinless photon' proceeds predominantly through Higgs boson exchange, and is largely independent of other Kaluza-Klein particles. The measured relic abundance sets an upper limit on the spinless photon mass of 500 GeV, which decreases to almost 200 GeV if the Higgs boson is light. The phenomenology of this dark matter candidate is strikingly different from Kaluza-Klein dark matter in theories with one universal extra dimension. Elastic scattering of the spinless photon with quarks is helicity suppressed, making its direct detection challenging, although possible at upcoming experiments. The prospects for indirect detection with gamma rays and antimatter are similar to those of neutralinos. The rates predicted at neutrino telescopes are below the sensitivity of next-generation experiments.

Observations by the WMAP experiment have identified an excess of microwave emission from the center of the Milky Way. It has previously been shown that this "WMAP Haze" could be synchrotron emission from relativistic electrons and positrons produced in the annihilations of dark matter particles. In particular, the intensity, spectrum and angular distribution of the WMAP Haze is consistent with an electroweak scale dark matter particle (such as a supersymmetric neutralino or Kaluza-Klein dark matter in models with universal extra dimensions) annihilating with a cross section on the order of sigma v~3x10^-26 cm^3/s and distributed with a cusped halo profile. No further exotic astrophysical or annihilation boost factors are required. If dark matter annihilations are in fact responsible for the observed Haze, then other annihilation products will also be produced, including gamma rays. In this article, we study the prospects for the GLAST satellite to detect gamma rays from dark matter annihilations in the Galactic Center region in this scenario. We find that by studying only the inner 0.1 degrees around the Galactic Center, GLAST will be able to detect dark matter annihilating to heavy quarks or gauge bosons over astrophysical backgrounds with 5sigma (3sigma) significance if they are lighter than approximately 320-500 GeV (500-750 GeV). If the angular window is broadened to study the dark matter halo profile's angular extension (while simultaneously reducing the astrophysical backgrounds), WIMPs as heavy as several TeV can be identified by GLAST with high significance. Only if the dark matter particles annihilate mostly to electrons or muons will GLAST be unable to identify the gamma ray spectrum associated with the WMAP Haze.

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.

We review the phenomenology of models with flat, compactified extra dimensions where all of the Standard Model fields are allowed to propagate in the bulk, known as Universal Extra Dimensions (UED). UED make for an interesting TeV-scale physics scenario, featuring a tower of Kaluza-Klein (KK) states approximately degenerate in mass at the scale set by the inverse size of the compactification radius. KK parity, the four-dimensional remnant of momentum conservation in the extra dimensions, implies two basic consequences: (1) contributions to Standard Model observables arise only at loop level, and KK states can only be pair-produced at colliders, and (2) the lightest KK particle (LKP) is stable, providing a suitable particle dark matter candidate. After a theoretical overview on extra dimensional models, and on UED in particular, we introduce the model particle spectrum and the constraints from precision electroweak tests and current colliders data. We then give a detailed overview of the LKP dark matter phenomenology, including the LKP relic abundance, and direct and indirect searches. We then discuss the physics of UED at colliders, with particular emphasis on the signatures predicted for the Large Hadron Collider and at a future Linear Collider, as well as on the problem of discriminating between UED and other TeV-scale new physics scenarios, particularly supersymmetry. We propose a set of reference benchmark models, representative of different viable UED realizations. Finally, we collect in the Appendix all the relevant UED Feynman rules, the scattering cross sections for annihilation and coannihilation processes in the early universe and the production cross section for strongly interacting KK states at hadron colliders.

Dark matter annihilating in our Galaxy's halo and elsewhere in the universe is expected to generate a diffuse flux of gamma rays, potentially observable with next generation satellite-based experiments, such as GLAST. In this article, we study the signatures of dark matter in the angular distribution of this radiation, in particular the deterministic ones. Pertaining to the extragalactic contribution, we discuss the effect of the motion of the solar system with respect to the cosmological rest frame, and anisotropies due to the structure of our local universe. For the gamma ray flux from dark matter in our own Galactic halo, we discuss the effects of the offset position of the solar system, the Compton-Getting effect, the asphericity of the Milky Way halo, and the signatures of nearby substructure. We explore the prospects for the detection of these features by the GLAST satellite and find that, if ~10% or more of the diffuse gamma ray background observed by EGRET is the result of dark matter annihilations, then GLAST should be sensitive to anisotropies at the sub-percent level. Such precision would be sufficient to detect many, if not all, of the signatures discussed in this paper.

The WMAP experiment has revealed an excess of microwave emission from the region around the center of our Galaxy. It has been suggested that this signal, known as the ``WMAP Haze'', could be synchrotron emission from relativistic electrons and positrons generated in dark matter annihilations. In this letter, we revisit this possibility. We find that the angular distribution of the WMAP Haze matches the prediction for dark matter annihilations with a cusped density profile, $\rho(r) \propto r^{-1.2}$ in the inner kiloparsecs. Comparing the intensity in different WMAP frequency bands, we find that a wide range of possible WIMP annihilation modes are consistent with the spectrum of the haze for a WIMP with a mass in the 100 GeV to multi-TeV range. Most interestingly, we find that to generate the observed intensity of the haze, the dark matter annihilation cross section is required to be approximately equal to the value needed for a thermal relic, $\sigma v \sim 3 \times 10^{-26}$ cm$^3$/s. No boost factors are required. If dark matter annihilations are in fact responsible for the WMAP Haze, and the slope of the halo profile continues into the inner Galaxy, GLAST is expected to detect gamma rays from the dark matter annihilations in the Galactic Center if the WIMP mass is less than several hundred GeV.

WIMPs with electroweak scale masses (neutralinos, etc.) remain in kinetic equilibrium with other particle species until temperatures approximately in the range of 10 MeV to 1 GeV, leading to the formation of dark matter substructure with masses as small as $10^{-4} M_{\odot}$ to $10^{-12} M_{\odot}$. However, if dark matter consists of particles with MeV scale masses, as motivated by the observation of 511 keV emission from the Galactic Bulge, such particles are naturally expected to remain in kinetic equilibrium with the cosmic neutrino background until considerably later times. This would lead to a strong suppression of small scale structure with masses below about $10^7 M_{\odot}$ to $10^4 M_{\odot}$. This cutoff scale has important implications for present and future searches for faint Local Group satellite galaxies and for the missing satellites problem.

If annihilating MeV-scale dark matter particles are responsible for the observed 511 keV emission from the Galactic bulge, then new light gauge bosons which mediate the dark matter annihilations may have other observable consequences. In particular, if such a gauge boson exists and has even very small couplings to Standard Model neutrinos, cosmic neutrinos with ~TeV energies will scatter with the cosmic neutrino background through resonant exchange, resulting in a distinctive spectral absorption line in the high-energy neutrino spectrum. Such a feature could potentially be detected by future high-energy neutrino telescopes.

There are compelling reasons to think that new physics will appear at or below the TeV-scale. It is not known what form this new physics will take, however. Although The Large Hadron collider is very likely to discover new particles associated with the TeV-scale, it may be difficult for it to determine the nature of those particles, whether superpartners, Kaluza-Klein modes or other states. In this article, we consider how direct and indirect dark matter detection experiments may provide information complementary to hadron colliders, which can be used to discriminate between supersymmetry, models with universal extra dimensions, and Little Higgs theories. We find that, in many scenarios, dark matter experiments can be effectively used to distinguish between these possibilities.

In this article, we explore the interplay between searches for supersymmetric particles and Higgs bosons at hadron colliders (the Tevatron and the LHC) and direct dark matter searches (such as CDMS, ZEPLIN, XENON, EDELWEISS, CRESST, WARP and others). We focus on collider searches for heavy MSSM Higgs bosons ($A$, $H$, $H^{\pm}$) and how the prospects for these searches are impacted by direct dark matter limits and vice versa. We find that the prospects of these two experimental programs are highly interrelated. A positive detection of $A$, $H$ or $H^{\pm}$ at the Tevatron would dramatically enhance the prospects for a near future direct discovery of neutralino dark matter. Similarly, a positive direct detection of neutralino dark matter would enhance the prospects of discovering heavy MSSM Higgs bosons at the Tevatron or the LHC. Combining the information obtained from both types of experimental searches will enable us to learn more about the nature of supersymmetry.

We investigate the propagation of ultra-high energy cosmic ray nuclei (A = 1-56) from cosmologically distant sources through the cosmic radiation backgrounds. Various models for the injected composition and spectrum and of the cosmic infra-red background are studied using updated photo-disintegration cross-sections. The observational data on the spectrum and the composition of ultra-high energy cosmic rays are jointly consistent with a model where all of the injected primary cosmic rays are iron nuclei (or a mixture of heavy and light nuclei).

We point out that the usual self-similarity in cold dark matter models is broken by encounters with individual normal galactic stars on sub-pc scale. Tidal heating and stripping must have redefined the density and velocity structures of the population of the Earth-mass dark matter halos, which are likely to have been the first bound structures to form in the Universe. The disruption rate depends strongly on {\it galaxy types} and the orbital distribution of the microhalos; in the Milky Way, stochastic radial orbits are destroyed first by stars in the triaxial bulge, microhalos on non-planar retrograde orbits with large pericenters and/or apocenters survive the longest. The final microhalo distribution in the {\it solar neighborhood} is better described as a superposition of filamentry microstreams rather than as a set of discrete spherical clumps in an otherwise homogeneous medium. We discuss its important consequences to our detections of microhalos by direct recoil signal and indirect annihilation signal.

Searches for the Minimal Supersymmetric Standard Model (MSSM) Higgs bosons are among the most promising channels for exploring new physics at the Tevatron. In particular, interesting regions of large $\tan \beta$ and small $m_A$ are probed by searches for heavy neutral Higgs bosons, A and H, when they decay to $\tau^+ \tau^-$ and $b\bar{b}$. At the same time, direct searches for dark matter, such as CDMS, attempt to observe neutralino dark matter particles scattering elastically off nuclei. This can occur through t-channel Higgs exchange, which has a large cross section in the case of large $\tan \beta$ and small $m_A$. As a result, there is a natural interplay between the heavy, neutral Higgs searches at the Tevatron and the region of parameter space explored by CDMS. We show that if the lightest neutralino makes up the dark matter of our universe, current limits from CDMS strongly constrain the prospects of heavy, neutral MSSM Higgs discovery at the Tevatron (at 3 sigma with 4 fb^-1 per experiment) unless $|\mu| \gsim$ 400 GeV. The limits of CDMS projected for 2007 will increase this constraint to $|\mu| \gsim$ 800 GeV. On the other hand, if CDMS does observe neutralino dark matter in the near future, it will make the discovery of heavy, neutral MSSM Higgs bosons far more likely at the Tevatron.

The sources of the observed ultra-high energy cosmic rays must also generate ultra-high energy neutrinos. Deep inelastic scattering of these neutrinos with nucleons on Earth probe center-of-mass energies $\sqrt{s} \sim 100$ TeV, well beyond those attainable at terrestrial colliders. By comparing the rates for two classes of observable events, any departure from the benchmark (unscreened perturbative QCD) neutrino-nucleon cross-section can be constrained. Using the projected sensitivity of the Pierre Auger Observatory to quasi-horizontal showers and Earth-skimming tau neutrinos, we show that a `Super-Auger' detector can thus provide an unique probe of strong interaction dynamics.

In this article, we explore the ability of direct and indirect dark matter experiments to not only detect neutralino dark matter, but to constrain and measure the parameters of supersymmetry. In particular, we explore the relationship between the phenomenological quantities relevant to dark matter experiments, such as the neutralino annihilation and elastic scattering cross sections, and the underlying characteristics of the supersymmetric model, such as the values of mu (and the composition of the lightest neutralino), m_A and tan beta. We explore a broad range of supersymmetric models and then focus on a smaller set of benchmark models. We find that by combining astrophysical observations with collider measurements, mu can often be constrained far more tightly than it can be from LHC data alone. In models in the A-funnel region of parameter space, we find that dark matter experiments can potentially determine m_A to roughly +/-100 GeV, even when heavy neutral MSSM Higgs bosons (A, H_1) cannot be observed at the LHC. The information provided by astrophysical experiments is often highly complementary to the information most easily ascertained at colliders.

A number of experimental techniques are currently being deployed in an effort to make the first detection of ultra-high energy cosmic neutrinos. To accomplish this goal, techniques using radio and acoustic detectors are being developed, which are optimally designed for studying neutrinos with energies in the PeV-EeV range and above. Data from the AMANDA experiment, in contrast, has been used to place limits on the cosmic neutrino flux at less extreme energies (up to ~10 PeV). In this letter, we show that by adopting a different analysis strategy, optimized for much higher energy neutrinos, the same AMANDA data can be used to place a limit competitive with radio techniques at EeV energies. We also discuss the sensitivity of the IceCube experiment, in various stages of deployment, to ultra-high energy neutrinos.

A wide range of techniques have been developed to search for particle dark matter, including direct detection, indirect detection, and collider searches. The prospects for the detection of neutralino dark matter is quite promising for each of these three very different methods. Looking ahead to a time in which these techniques have successfully detected neutralino dark matter, we explore the ability of these observations to determine the parameters of supersymmetry. In particular, we focus on the ability of direct and indirect detection techniques to measure the parameters $\mu$ and $m_A$. We find that $\mu$ can be much more tightly constrained if astrophysical measurements are considered than by LHC data alone. In supersymmetric models within the $A$-funnel region of parameter space, we find that astrophysical measurements can determine $m_A$ to roughly $\pm100$ GeV precision.

Atmospheric Cerenkov Telescopes, including HESS and MAGIC, have detected a spectrum of gamma-rays from the galactic center region which extends from $\sim$200 GeV or lower, to at least $\sim$10 TeV. Although the source of this radiation is not yet known, the spectrum appears to behave as a simple power-law, which is not the expectation for gamma-rays generated through the annihilation of dark matter particles. If instead we conclude that the source of these gamma-rays is astrophysical in origin, this spectrum will constitute a background for future dark matter searches using gamma-rays from this region. In this paper we study how this background will affect the prospects for experiments such as GLAST to detect dark matter in the galactic center. We find that only a narrow range of dark matter annihilation rates are potentially observable by GLAST given this newly discovered background and considering current constraints from EGRET and HESS. We also find that a detection of line emission, while not completely ruled out, is only possible for a very narrow range of dark matter models and halo profiles.

We consider the ultra-high energy cosmic ray (UHECR) spectrum as measured by the Pierre Auger Observatory. Top-down models for the origin of UHECRs predict an increasing photon component at energies above about $10^{19.7}$eV. Here we present a simple prescription to compare the Auger data with a prediction assuming a pure proton component or a prediction assuming a changing primary component appropriate for a top-down model. We find that the UHECR spectrum predicted in top-down models is a good fit to the Auger data. Eventually, Auger will measure a composition-independent spectrum and will be capable of either confirming or excluding the quantity of photons predicted in top-down models.

We show that light ($\simeq$ 1 -- 30 MeV) dark matter particles can play a significant role in core-collapse supernovae, if they have relatively large annihilation and scattering cross sections, as compared to neutrinos. We find that if such particles are lighter than $\simeq$ 10 MeV and reproduce the observed dark matter relic density, supernovae would cool on a much longer time scale and would emit neutrinos with significantly smaller energies than in the standard scenario, in disagreement with observations. This constraint may be avoided, however, in certain situations for which the neutrino--dark matter scattering cross sections remain comparatively small.

The dwarf spheroidal galaxy Draco has long been considered likely to be one of the brightest point sources of gamma-rays generated through dark matter annihilations. Recent studies of this object have found that it remains largely intact from tidal striping, and may be more massive than previously thought. In this article, we revisit Draco as a source of dark matter annihilation radiation, with these new observational constraints in mind. We discuss the prospects for the experiments MAGIC and GLAST to detect dark matter in Draco, as well as constraints from the observations of EGRET. We also discuss the possibility that the CACTUS experiment has already detected gamma-rays from Draco. We find that it is difficult to generate the flux reported by CACTUS without resorting to non-thermally produced WIMPs and/or a density spike in Draco's dark matter distribution due to the presence of an intermediate mass black hole. We also find that for most annihilation modes, a positive detection of Draco by CACTUS would be inconsistent with the lack of events seen by EGRET.

It has been suggested that Gamma Ray Bursts (GRB) may enable the expansion rate of our Universe to be measured out to very high redshifts ($z \gsim 5$) just as type Ia supernovae have done at $z \sim$1--1.5. We explore this possibility here, and find that GRB have the potential to detect dark energy at high statistical significance, but they are unlikely to be competitive with future supernovae missions, such as SNAP, in measuring the properties of the dark energy. The exception to this conclusion is if there is appreciable dark energy at early times, in which case the information from GRB's will provide an excellent complement to the $z\sim 1$ information from supernovae.

We study constraints on models with a flat "Universal'' Extra Dimension in which all Standard Model fields propagate in the bulk. A significantly improved constraint on the compactification scale is obtained from the extended set of electroweak precision observables accurately measured at LEP1 and LEP2. We find a lower bound of M_c = R^{-1} > 700 (800) GeV at the 99% (95%) confidence level. We also discuss the implications of this constraint on the prospects for the direct and indirect detection of Kaluza-Klein dark matter in this model.

The Pierre Auger Observatory for cosmic rays provides a laboratory for studying fundamental interactions at energies well beyond those available at colliders. In addition to hadrons or photons, Auger is sensitive to ultra-high energy neutrinos in the cosmic radiation and models for new physics can be explored by observing neutrino interactions at center-of-mass energies beyond the TeV scale. By comparing the rate for quasi-horizontal, deeply penetrating air showers triggered by all types of neutrinos with the rate for slightly upgoing showers generated by Earth-skimming tau neutrinos, any deviation of the neutrino-nucleon cross-section from the Standard Model expectation can be constrained. We show that this can test models of low-scale quantum gravity (including processes such as Kaluza-Klein graviton exchange, microscopic black hole production and string resonances), as well as non-perturbative electroweak instanton mediated processes. Moreover, the observed ratios of neutrino flavors would severely constrain the possibility of neutrino decay.

Traditionally, collider experiments have been the primary tool used in searching for particle physics beyond the Standard Model. In this talk, I will discuss alternative approaches for exploring exotic physics scenarios using high energy and ultra-high energy cosmic neutrinos. Such neutrinos can be used to study interactions at energies higher, and over baselines longer, than those accessible to colliders. In this way, neutrino astronomy can provide a window into fundamental physics which is highly complementary to collider techniques. I will discuss the role of neutrino astronomy in fundamental physics, considering the use of such techniques in studying several specific scenarios including low scale gravity models, Standard Model electroweak instanton induced interactions, decaying neutrinos and quantum decoherence.

Direct detection dark matter experiments, lead by the CDMS collaboration, have placed increasingly stronger constraints on the cross sections for elastic scattering of WIMPs on nucleons. These results impact the prospects for the indirect detection of dark matter using neutrino telescopes. With this in mind, we revisit the prospects for detecting neutrinos produced by the annihilation of WIMPs in the Sun. We find that the latest bounds do not seriously limit the models most accessible to next generation kilometer-scale neutrino telescopes such as IceCube. This is largely due to the fact that models with significant spin-dependent couplings to protons are the least constrained and, at the same time, the most promising because of the efficient capture of WIMPs in the Sun. We identify models where dark matter particles are beyond the reach of any planned direct detection experiments while within reach of neutrino telescopes. In summary, we find that, even when contemplating recent direct detection results, neutrino telescopes still have the opportunity to play an important as well as complementary role in the search for particle dark matter.

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.

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.

High-energy neutrino astronomy will be capable of observing particles at both extremely high energies and over extremely long baselines. These features make such experiments highly sensitive to the effects of CPT and Lorentz violation. In this article, we review the theoretical foundation and motivation for CPT and Lorentz violating effects, and then go on to discuss the related phenomenology within the neutrino sector. We describe several signatures which might be used to identify the presence of CPT or Lorentz violation in next generation neutrino telescopes and cosmic ray experiments. In many cases, high-energy neutrino experiments can test for CPT and Lorentz violation effects with much greater precision than other techniques.

Recently, it has been pointed out that two different excesses of events observed at LEP could be interpreted as the CP-even Higgs bosons of the MSSM with masses of approximately 98 and 114 GeV. If this is the case, the entire MSSM Higgs sector is required to be light. In this article, we explore such a scenario in detail. We constrain the Higgs and supersymmetric spectrum using $B$ physics constraints as well as the magnetic moment of the muon. We then point out the implications for neutralino dark matter -- next generation direct detection experiments will be sensitive to all MSSM models with such a Higgs sector. Finally, we find that all models outside of a very narrow corridor of the parameter space have a charged Higgs boson which will be observed at the LHC. In those exceptional models which do not contain an observable charged Higgs, a light top squark will always be seen at the LHC, and likely at the Tevatron.

The AMANDA neutrino telescope has recently reported the detection of high-energy neutrinos spatially and temporally coincident with the flaring of the TeV blazar 1ES 1959+650. At present, the statistical significance of this observation cannot be reliably assessed, however. In this letter, we investigate whether circumstances exist where the source can produce the flux implied by the coincident events. We show that if the TeV gamma-ray emission observed from 1ES 1959+650 or other nearby TeV blazars is the result of accelerated protons interacting with nucleons, it is reasonable that AMANDA could detect several events during a flaring period. Such rates require that the spectral index of the source be rather high (for instance $\sim2.8$ for 1ES 1959+650) and that the Lorentz factor of the jet be fairly small ($\Gamma \sim 1$).

We present the signatures and prospects for the indirect detection of a Dirac right-handed neutrino dark matter candidate in neutrino telescopes, cosmic positron experiments and gamma-ray telescopes. An example of such a dark matter candidate can be found in extra-dimensional models. In some constructions, Kaluza--Klein states with the gauge quantum numbers of a right-handed neutrino can have sizable gauge interactions with Standard Model particles. For instance, in 5D warped Grand Unified Theories, it has been shown that a Kaluza--Klein right-handed neutrino may be stable and otherwise a phenomenologically viable dark matter candidate. We find that the prospects for the indirect detection of such a WIMP are encouraging, particularly for neutrino telescopes and cosmic positron experiments.

The promising neutralino dark matter particles generically condense into numerous earth-mass dark halos with smaller substructures suppressed by free streaming\cite{berezinsky}. The recent {\it Letter to Nature}\cite{moore} claims that these 0.01pc-sized dense halos emerged at redshifts 60--26 are rarely destructed inside galaxies, hence the nearest halos at about 0.1pc from Earth are bright in annihilation-powered gamma-rays, but are inaccessible to direct detections. However, most mini-halos reaching solar neighbourhood should experience strong impulses by individual stars in the Galactic disk and bulge, and have been torn into pc-long tidal streams over a Hubble time with only modest overdensity, reducing indirect annihilation signals. Sweeping across the solar system per century, mini-streams leave directional and temporal signatures for direct searches of neutralinos.

The inferred values of the cosmological baryon and dark matter densities are strikingly similar, but in most theories of the early universe there is no true explanation of this fact; in particular, the baryon asymmetry and thus density depends upon unknown, and {\it a priori} unknown and possibly small, CP-violating phases which are independent of all parameters determining the dark matter density. We consider models of dark matter possessing a particle-antiparticle asymmetry where this asymmetry determines both the baryon asymmetry and strongly effects the dark matter density, thus naturally linking $\Omega_{\rm{b}}$ and $\Omega_{\rm{dm}}$. We show that sneutrinos can play the role of such dark matter in a previously studied variant of the MSSM in which the light neutrino masses result from higher-dimensional supersymmetry-breaking terms.

Recently, the HESS telescope has observed a gamma-ray spectrum from the galactic center extending to at least $\sim$10 TeV. Although it has been suggested that this could be the product of annihilating dark matter particles, the candidates most frequently discussed (such as neutralinos) are far too light to account for this flux. In this letter, we consider stable particles from the messenger sector of gauge mediated supersymmetry breaking models as an alternative dark matter candidate. We find that a 20 to 30 TeV messenger state can provide a thermal relic density consistent with the measured dark matter density of the universe and can indeed generate the spectrum observed by HESS.

It has been recently proposed that the observed 511 keV emission from the Galactic bulge could be the product of very light (1-100 MeV) annihilating dark matter particles. Other possible explanations for this signal are associated with stellar objects, such as hypernovae. In order to distinguish between annihilating light dark matter scenario and more conventional astrophysical sources for the bulge emission, we here propose the study of dwarf spheroidals such as Sagittarius. These galaxies have typical luminosities of $10^5 L_\odot$ but mass-to-light ratios of $\sim 100$. As there are comparatively few stars, the prospects for 511 keV emission from standard astrophysical scenarios are minimal. The dwarf spheroidals do, however, contain copious amounts of dark matter. INTEGRAL/SPI has observed the Sagittarius region. Analysis of this data for 511 keV emission will provide a test of MeV dark matter which can distinguish between annihilating dark matter and more standard astrophysics. The observation of such a signal from Sagittarius should be a ``smoking gun'' for MeV dark matter.

Recently, the SPI spectrometer on the INTEGRAL satellite observed strong 511 keV line emission from the galactic bulge. Although the angular distribution (spherically symmetric with width of \sim 9 degree) of this emission is difficult to account for with traditional astrophysical scenarios, light dark matter particles could account for the observation. In this letter, we consider the possibility that decaying axinos in an R-parity violating model of supersymmetry may be the source of this emission. We find that \sim 1-300 MeV axinos with R-parity violating couplings can naturally produce the observed emission.

As ultra-high energy cosmic ray protons propagate through the universe, they undergo photo-meson interactions with the cosmic microwave background, generating the `cosmogenic' neutrino flux. If a substantial fraction of the cosmic ray primaries are heavy nuclei rather than protons, however, they would preferentially lose energy through photo-disintegration, so the corresponding neutrino flux may be substantially depleted. We investigate this issue using a Monte Carlo simulation of cosmic ray propagation through interagalactic radiation fields and assess the impact of the altered neutrino fluxes on next generation neutrino telescopes.

Kaluza-Klein dark matter particles can annihilate efficiently into electron-positron pairs, providing a discrete feature (a sharp edge) in the cosmic $e^+ e^-$ spectrum at an energy equal to the particle's mass (typically several hundred GeV to one TeV). Although this feature is probably beyond the reach of satellite or balloon-based cosmic ray experiments (those that distinguish the charge and mass of the primary particle), gamma ray telescopes may provide an alternative detection method. Designed to observe very high-energy gamma-rays, ACTs also observe the diffuse flux of electron-induced electromagnetic showers. The GLAST satellite, designed for gamma ray astronomy, will also observe any high energy showers (several hundred GeV and above) in its calorimeter. We show that high-significance detections of an electron-positron feature from Kaluza-Klein dark matter annihilations are possible with GLAST, and also with ACTs such as HESS, VERITAS or MAGIC.

We discuss the prospects for high-energy neutrino astronomy to study particle physics in the energy regime comparable to and beyond that obtainable at the current and planned colliders. We describe the various signatures of high-energy cosmic neutrinos expected in both neutrino telescopes and air shower experiments and discuss these measurements within the context of theoretical models with a quantum gravity or string scale near a TeV, supersymmetry and scenarios with interactions induced by electroweak instantons. We attempt to access the particle physics reach of these experiments.

We consider the prospects for observing the effects of quantum decoherence in high-energy (TeV-PeV) neutrinos from astrophysical sources. In particular, we study Galactic sources of electron anti-neutrinos produced in the decay of ultra-high energy neutrons. We find that next generation neutrino telescopes should be capable of placing limits on quantum decoherence effects over multi-kiloparsec baselines, surpassing current bounds by a factor of $10^{12}$ to $10^{33}$, depending on the model considered.

In this talk, we discuss the potential for the indirect detection of Kaluza-Klein dark matter using neutrino telescopes and cosmic positron experiments. We find that future kilometer-scale neutrino telescopes, such as IceCube, as well as future experiments capable of measuring the cosmic positron spectrum, such as PAMELA and AMS-02, will be quite sensitive to this scenario. Current data from the HEAT experiment can also be explained by the presence of Kaluza-Klein dark matter in the Galactic halo.

Dark matter particles annihilating in the Galactic halo can provide a flux of positrons potentially observable in upcoming experiments, such as PAMELA and AMS-02. We discuss the spectral features which may be associated with dark matter annihilation in the positron spectrum and assess the prospects for observing such features in future experiments. Although we focus on some specific dark matter candidates, neutralinos and Kaluza-Klein states, we carry out our study in a model independent fashion. We also revisit the positron spectrum observed by HEAT.

In this review article, we discuss the current status of particle dark matter, including experimental evidence and theoretical motivations. We discuss a wide array of candidates for particle dark matter, but focus on neutralinos in models of supersymmetry and Kaluza-Klein dark matter in models of universal extra dimensions. We devote much of our attention to direct and indirect detection techniques, the constraints placed by these experiments and the reach of future experimental efforts.

The excess of cosmic positrons observed by the HEAT experiment may be the result of Kaluza-Klein dark matter annihilating in the galactic halo. Kaluza-Klein dark matter annihilates dominantly into charged leptons that yield a large number and hard spectrum of positrons per annihilation. Given a Kaluza-Klein dark matter particle with a mass in the range of 300-400 GeV, no exceptional substructure or clumping is needed in the local distribution of dark matter to generate a positron flux that explains the HEAT observations. This is in contrast to supersymmetric dark matter that requires unnaturally large amounts of dark substructure to produce the observed positron excess. Future astrophysical and collider tests are outlined that will confirm or rule out this explanation of the HEAT data.

Two ground-based experiments have recently independently detected TeV $\gamma$-rays from the direction of the Galactic center. The observations made by the VERITAS and CANGAROO collaborations are unexpected, although not impossible to interpret in terms of astrophysical sources. Here we examine in detail whether the observed $\gamma$-rays may arise from the more exotic alternative of annihilations of dark matter particles clustered in the center of the Galaxy.

One of the most interesting mysteries of astrophysics is the puzzle of dark matter. Although numerous techniques have been explored and developed to detect this elusive substance, its nature remains unknown. One such method uses large high-energy neutrino telescopes to look for the annihilation products of dark matter annihilations. In this summary article, we briefly review this technique. We describe the calculations used to find the rate of capture of WIMPs in the Sun or Earth and the spectrum of neutrinos produced in the resulting dark matter annihilations. We will discuss these calculations within the context of supersymmetry and models with universal extra dimensions, the lightest supersymmetric particle and lightest Kaluza-Klein particle providing the WIMP candidate in these cases, respectively. We will also discuss the status of some of the experiments relevant to these searches: AMANDA, IceCube and ANTARES.

It has often been suggested that the cosmic positron excess observed by the HEAT experiment could be the consequence of supersymmetric dark matter annihilating in the galactic halo. Although it is well known that evenly distributed dark matter cannot account for the observed excess, if substantial amounts of local dark matter substructure are present, the positron flux would be enhanced, perhaps to the observed magnitude. In this paper, we attempt to identify the nature of the substructure required to match the HEAT data, including the location, size and density of any local dark matter clump(s). Additionally, we attempt to assess the probability of such substructure being present. We find that if the current density of neutralino dark matter is the result of thermal production, very unlikely ($\sim 10^{-4}$ or less) conditions must be present in local substructure to account for the observed excess.

In this article, based on the talk given at the Cracow Epiphany Conference on Astroparticle Physics, I discuss some of the opportunities provided by high-energy and ultra-high energy neutrino astronomy in probing particle physics beyond the standard model. Following a short summary of current and next generation experiments, I review the prospects for observations of high-energy neutrino interactions, searches for particle dark matter, and measurements of absolute neutrino masses, lifetimes and pseudo-Dirac mass splittings.

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.

While the first kilometer-scale neutrino telescope, IceCube, is under construction, alternative plans exist to build even larger detectors that will, however, b e limited by a much higher neutrino energy threshold of 10 PeV or higher rather than 10 to 100 GeV. These future projects detect radio and acoustic pulses as w ell as air showers initiated by ultra-high energy neutrinos. As an alternative, we here propose an expansion of IceCube, using the same strings, placed on a gri d with a spacing of order 500 m. Unlike other proposals, the expanded detector uses methods that are understood and calibrated on atmospheric neutrinos. Atmosp heric neutrinos represent the only background at the energies under consideratio n and is totally negligible. Also, the cost of such a detector is understood. We conclude that supplementing the 81 IceCube strings with a modest number of addi tional strings spaced at large distances can almost double the effective volume of the detector. Doubling the number of strings on a 800 m grid can deliver a d etector that this a factor of 5 larger for horizontal muons at modest cost.

GRB941017, a gamma-ray burst of exceptional fluence, has recently been shown to have a high-energy component which is not consistent with the standard fireball phenomenology. If this component is the result of photomeson interactions in the burst fireball, it provides new and compelling support for substantial high-energy neutrino fluxes from this and similar sources. In this letter, we consider what impact this new information has on the neutrino spectra of gamma-ray bursts and discuss how this new evidence impacts the prospects for detection of such events in next generation neutrino telescopes.

We discuss the possibility that the recent detection of 511 keV gamma-rays from the galactic bulge, as observed by INTEGRAL, is a consequence of low mass (~MeV) particle dark matter annihilations. We discuss the type of halo profile favored by the observations as well as the size of the annihilation cross section needed to account for the signal. We find that such a scenario is consistent with the observed dark matter relic density and other constraints from astrophysics and particle physics.

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.

Various methods of searching for supersymmetric dark matter are sensitive to WIMPs with different properties. One consequence of this is that the phenomenology of dark matter detection can vary dramatically in different supersymmetric breaking scenarios. In this paper, we consider the sensitivities to supersymmetric dark matter of different detection methods and techniques in a wide variety of supersymmetric breaking scenarios. We discuss the ability of various astrophysical experiments, such as direct experiments, gamma-rays satellites, neutrino telescopes and positron and anti-proton cosmic ray experiments, to test various supersymmetry breaking scenarios. We also discuss what information can be revealed about supersymmetry breaking by combining results from complementary experiments. We place an emphasis on the differences between various experimental techniques.

We calculate the neutrino emission from individual gamma-ray bursts observed by the BATSE detector on the Compton Gamma-Ray Observatory. Neutrinos are produced by photoproduction of pions when protons interact with photons in the region where the kinetic energy of the relativistic fireball is dissipated allowing the acceleration of electrons and protons. We also consider models where neutrinos are predominantly produced on the radiation surrounding the newly formed black hole. From the observed redshift and photon flux of each individual burst, we compute the neutrino flux in a variety of models based on the assumption that equal kinetic energy is dissipated into electrons and protons. Where not measured, the redshift is estimated by other methods. Unlike previous calculations of the universal diffuse neutrino flux produced by all gamma-ray bursts, the individual fluxes (compiled at http://www.arcetri.astro.it/~dafne/grb/) can be directly compared with coincident observations by the AMANDA telescope at the South Pole. Because of its large statistics, our predictions are likely to be representative for future observations with larger neutrino telescopes.

We demonstrate that the South Pole kilometer-scale neutrino observatory IceCube can detect multi-TeV gamma rays continuously over a large fraction of the southern sky. While not as sensitive as pointing atmospheric Cerenkov telescopes, IceCube can roughly match the sensitivity of Milagro. Also, IceCube is complementary to Milagro because it will observe, without interruption, a relatively poorly studied fraction of the southern sky. The information which IceCube must record to function as a gamma ray observatory is only the directions and possibly energies of down-going muons.

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.

In models in which all of the Standard Model fields live in extra universal dimensions, the lightest Kaluza-Klein (KK) particle can be stable. Calculations of the one-loop radiative corrections to the masses of the KK modes suggest that the identity of the lightest KK particle (LKP) is mostly the first KK excitation of the hypercharge gauge boson. This LKP is a viable dark matter candidate with an ideal present-day relic abundance if its mass is moderately large, between 600 to 1200 GeV. Such weakly interacting dark matter particles are expected to become gravitationally trapped in large bodies, such as the Sun, and annihilate into neutrinos or other particles that decay into neutrinos. We calculate the annihilation rate, neutrino flux and the resulting event rate in present and future neutrino telescopes. The relatively large mass implies that the neutrino energy spectrum is expected to be well above the energy threshold of AMANDA and IceCube. We find that the event rate in IceCube is between a few to tens of events per year.

Using the EGRET data and an improved point source analysis, including an energy dependent point spread function and an unbinned maximum likelihood technique, we have been able to place considerably lower limits on the gamma ray flux from the galactic center region. We also test this method on known sources, the Crab and Vela pulsars. In both cases, we find that our method improves the angular precision of EGRET data over the 3EG catalog. This new limit on gamma rays from the galactic center can be used to test models of annihilating supersymmetric dark matter and galactic halo profiles. We find that the present EGRET data can limit many supersymmetric models if the density of the galactic dark matter halo is cuspy or spiked toward the galactic center. We also discuss the ability of GLAST to test these models.

In most supersymmetic models, neutralino dark matter particles are predicted to accumulate in the Galactic center and annihilate generating, among other products, gamma rays. The EGRET experiment has made observations in this region, and is sensitive to gamma rays from 30 MeV to $\sim$30 GeV. We have used an improved point source analysis including an energy dependent point spread function and an unbinned maximum likelihood technique, which has allowed us to significantly lower the limits on gamma ray flux from the Galactic center. We find that the present EGRET data can limit many supersymmetric models if the density of the Galactic dark matter halo is cuspy or spiked toward the Galactic center. We also discuss the ability of GLAST to test these models.

We construct the open-string scattering amplitudes for neutrino-gluon collisions and evaluate the high energy neutrino-nucleon scattering cross section via string state excitations in the TeV string-scale scenario. We find that the neutrino-gluon scattering is the dominant contribution, 5-10 times larger than neutrino-quark processes, though black hole production may be larger than the string contribution at higher energies. We illustrate the observability of the string signal at the Auger Observatory and the IceCube neutrino telescope for a string scale about 1 TeV.

We investigate the potential of a future kilometer-scale neutrino telescope such as the proposed IceCube detector in the South Pole, to measure and disentangle the yet unknown components of the cosmic neutrino flux, the prompt atmospheric neutrinos coming from the decay of charmed particles and the extra-galactic neutrinos, in the 10 TeV to 1 EeV energy range. Assuming a power law type spectra, $d\phi_\nu/dE_\nu \sim \alpha E_\nu^\beta$, we quantify the discriminating power of the IceCube detector and discuss how well we can determine magnitude ($\alpha$) as well as slope ($\beta$) of these two components of the high energy neutrino spectrum, taking into account the background coming from the conventional atmospheric neutrinos.

This is a review of neutrino astronomy anchored to the observational fact that Nature accelerates protons and photons to energies in excess of $10^{20}$ and $10^{13}$ eV, respectively. Although the discovery of cosmic rays dates back close to a century, we do not know how and where they are accelerated. Basic elementary-particle physics dictates a universal upper limit on their energy of $5\times10^{19}$ eV, the so-called Greisen-Kuzmin-Zatsepin cutoff; however, particles in excess of this energy have been observed by all experiments, adding one more puzzle to the cosmic ray mystery. Mystery is fertile ground for progress: we will review the facts as well as the speculations about the sources including gamma ray bursts, blazars and top-down scenarios. The important conclusion is that, independently of the specific blueprint of the source, it takes a kilometer-scale neutrino observatory to detect the neutrino beam associated with the highest energy cosmic rays and gamma rays. We also briefly review the ongoing efforts to commission such instrumentation.

Top-down models of cosmic rays produce more neutrinos than photons and more photons than protons. In these models, we reevaluate the fluxes of neutrinos associated with the highest energy cosmic rays in light of mounting evidence that they are protons and not gamma rays. While proton dominance at EeV energies can possibly be achieved by efficient absorption of the dominant high-energy photon flux on universal and galactic photon and magnetic background fields, we show that the associated neutrino flux is inevitably increased to a level where it should be within reach of operating experiments such as AMANDA II, RICE and AGASA. In future neutrino telescopes, tens to a hundred, rather than a few neutrinos per kilometer squared per year, may be detected above 1 PeV.

Next generation kilometer-scale neutrino telescopes, such as ICECUBE, can test standard model predictions for neutrino-nucleon cross sections at energies well beyond the reach of collider experiments. At energies near a PeV and higher, the Earth becomes opaque to neutrinos. At these energies, the ratio of upgoing and downgoing events can be used to measure the total neutrino-nucleon cross section given the presence of an adequate high energy neutrino flux.

In light of the mounting evidence that the highest energy cosmic rays are dominated by protons and not gamma-rays, we discuss the prospect that these cosmic rays are generated in the decay or annihilation of superheavy relic particles. We calculate the high energy neutrino spectrum which results and normalize our results to the ultra-high energy cosmic ray spectrum. We show that most scenerios are already constrained by present limits placed by the AMANDA experiment.

If spacetime has more than four dimensions, ultra-high energy cosmic rays may create microscopic black holes. Black holes created by cosmic neutrinos in the Earth will evaporate, and the resulting hadronic showers, muons, and taus may be detected in neutrino telescopes below the Earth's surface. We simulate such events in detail and consider black hole cross sections with and without an exponential suppression factor. We find observable rates in both cases: for conservative cosmogenic neutrino fluxes, several black hole events per year are observable at the IceCube detector; for fluxes at the Waxman-Bahcall bound, tens of events per year are possible. We also present zenith angle and energy distributions for all three channels. The ability of neutrino telescopes to differentiate hadrons, muons, and possibly taus, and to measure these distributions provides a unique opportunity to identify black holes, to experimentally constrain the form of black hole production cross sections, and to study Hawking evaporation.

We investigate the prospects of indirect searches for supersymmetric neutralino dark matter. Relic neutralinos gravitationally accumulate in the Sun and their annihilations produce high energy neutrinos. Muon neutrinos of this origin can be seen in large detectors like AMANDA, IceCube and ANTARES. We evaluate the relic density and the detection rate in several models---the minimal supersymmetric model, minimal supergravity, and supergravity with non-universal Higgs boson masses at the grand unification scale. We make realistic estimates for the indirect detection rates including effects of the muon detection threshold, quark hadronization, and solar absorption. We find good prospects for detection of neutralinos with mass above 200 GeV.

Following Waxman and Bahcall we calculate the event rate, energy and zenith angle dependence of neutrinos produced in the fireball model of gamma ray bursts (GRB). We emphasize the primary importance of i) burst-to-burst fluctuations and ii) absorption of the neutrinos in the Earth. From the astronomical point of view, we draw attention to the sensitivity of neutrino measurements to the boost Lorentz factor of the fireball $\Gamma$, which is central to the fireball model, and only indirectly determined by follow-up observations. Fluctuations result in single bursts emitting multiple neutrinos, making it possible to determine the flavor composition of a beam observed after a baseline of thousands of Megaparsecs.

We recalculate the diffuse flux of high energy neutrinos produced by Gamma Ray Bursts (GRB) in the relativistic fireball model. Although we confirm that the average single burst produces only ~10^{-2} high energy neutrino events in a detector with 1 km^2 effective area, i.e. about 10 events per year, we show that the observed rate is dominated by burst-to-burst fluctuations which are very large. We find event rates that are expected to be larger by one order of magnitude, likely more, which are dominated by a few very bright bursts. This greatly simplifies their detection.