The multitude of binary black hole coalescence detections in gravitational waves has renewed our interest on environments that can be the cradle of these mergers. In this work we study merger rates of binary black holes in globular clusters that are among the most dense stellar environments and a natural place for the creation of black hole binaries. To model these systems with all their variations we rely on the observational properties of the known Milky Way globular clusters. We consider direct capture events between black holes, as well as soft interactions of black hole binaries with stars as third bodies that accelerate the evolution of these binaries. We find that binary black holes from direct captures merge at an averaged rate of $0.3-5 \times 10^{-11}$ yr$^{-1}$ per cluster. Third body soft interactions are a much more prominent channel giving an averaged rate of $2-4 \times 10^{-10}$ yr$^{-1}$ per cluster. Those rates in globular clusters can lead to a cumulative merger rate of about 100 mergers per year up to redshift of 1, i.e. a significant fraction of the detectable in the near future binary black hole coalescence events. Further observations of cluster properties both in terms of their masses, profile properties, velocity dispersion of stars and their cosmological distribution, will allow us to better constrain the contribution of these environments to the detectable coalescence events rate.

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.

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.

The Fermi-LAT collaboration has recently released a new point source catalog, referred to as 4FGL. For the first time, we perform a template fit using information from this new catalog and find that the Galactic center excess is still present. On the other hand, we find that a wavelet-based search for point sources is highly sensitive to the use of the 4FGL catalog: no excess of bright regions on small angular scales is apparent when we mask out 4FGL point sources. We postulate that the 4FGL catalog contains the large majority of bright point sources that have previously been suggested to account for the excess in gamma rays detected at the Galactic center in Fermi-LAT data. Furthermore, after identifying which bright sources have no known counterpart, we place constraints on the luminosity function necessary for point sources to explain the smooth emission seen in the template fit.

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.

We discuss the origin of the anti-helium-3 and -4 events possibly detected by AMS-02. Using up-to-date semi-analytical tools, we show that spallation from primary hydrogen and helium nuclei onto the ISM predicts a $\overline{{}^3{\rm He}}$ flux typically one to two orders of magnitude below the sensitivity of AMS-02 after 5 years, and a $\overline{{}^4{\rm He}}$ flux roughly 5 orders of magnitude below the AMS-02 sensitivity. We argue that dark matter annihilations face similar difficulties in explaining this event. We then entertain the possibility that these events originate from anti-matter-dominated regions in the form of anti-clouds or anti-stars. In the case of anti-clouds, we show how the isotopic ratio of anti-helium nuclei might suggest that BBN has happened in an inhomogeneous manner, resulting in anti-regions with a anti-baryon-to-photon ratio $\bar{\eta}\simeq10^{-3}\eta$. We discuss properties of these regions, as well as relevant constraints on the presence of anti-clouds in our Galaxy. We present constraints from the survival of anti-clouds in the Milky-Way and in the early Universe, as well as from CMB, gamma-ray and cosmic-ray observations. In particular, these require the anti-clouds to be almost free of normal matter. We also discuss an alternative where anti-domains are dominated by surviving anti-stars. We suggest that part of the unindentified sources in the 3FGL catalog can originate from anti-clouds or anti-stars. AMS-02 and GAPS data could further probe this scenario.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe.

Ultra-light hidden-photon dark matter produces an oscillating electric field in the early Universe plasma, which in turn induces an electric current in its ionized component whose dissipation results in heat transfer from the dark matter to the plasma. This will affect the global 21cm signal from the Dark Ages and Cosmic Dawn. In this work we focus on the latter, in light of the reported detection by the EDGES collaboration of an absorption signal at frequencies corresponding to redshift z~17. By measuring the 21cm global signal, a limit can be placed on the amount of gas heating, and thus the kinetic mixing strength $\varepsilon$ between the hidden and ordinary photons can be constrained. Our inferred 21cm bounds on $\varepsilon$ in the mass range $10^{-23}\,{\rm eV}\lesssim m_\chi\lesssim10^{-13}\,{\rm eV}$ are the strongest to date.

Pulsars have been invoked to explain the origin of recently observed high-energy Galactic cosmic-ray positrons. Since the positron propagation distance decreases with energy, the number of pulsars that can contribute to the observed positrons decreases from $O(10^3)$ for positron energies $E\gtrsim10$ GeV to only a few for $E \gtrsim 500$ GeV. Thus, if pulsars explain these positrons, the positron energy spectrum should become increasingly bumpy at higher energies. Here we present a power-spectrum analysis that can be applied to seek such spectral features in the energy spectrum for cosmic-ray positrons and for the energy spectrum of the combined electron/positron flux. We account for uncertainties in the pulsar distribution by generating hundreds of simulated spectra from pulsar distributions consistent with current observational constraints. Although the current AMS-02 data do not exhibit evidence for spectral features, we find that such features would be detectable in $\simeq 10\%$ of our simulations, with 20 years of AMS-02 data or three years of DAMPE measurements on the electron-plus-positron flux.

Recent measurements of cosmic-ray electron and positron spectra at energies from a GeV to 5 TeV, as well as radio, X-ray and a wide range of gamma-ray observations of pulsar-wind nebulae, indicate that pulsars are significant sources of high-energy cosmic-ray electrons and positrons. Here we calculate the local cosmic-ray $e^\pm$ energy spectra from pulsars taking into account models for (a) the distribution of the pulsars spin-down properties; (b) the cosmic-ray source spectra; and (c) the physics of cosmic-ray propagation. We then use the measured cosmic-ray fluxes from AMS-02, CALET and DAMPE to constrain the space of pulsar and cosmic-ray-propagation models and in particular, local cosmic-ray diffusion and energy losses, the pulsars' energy-loss time-dependence, and the injected $e^{\pm}$ spectra. We find that the lower estimates for the local $e^{\pm}$ energy losses are inconsistent with the data. We also find that pulsar braking indexes of 2.5 or less for sources with ages more than 10 kyr are strongly disfavored. Moreover, the cosmic-ray data are consistent with a wide range of assumptions on the $e^{\pm}$ injection spectral properties and on the distribution of initial spin-down powers. Above a TeV in energy, we find that pulsars can easily explain the observed change in the $e^{+} + e^{-}$ spectral slope. These conclusions are valid as long as pulsars contribute $\gtrsim10\%$ of the observed cosmic-ray $e^\pm$ at energies $\gtrsim100$ GeV.

We analyze the gamma-ray sky at energies of 0.5 to 50 GeV using the undecimated wavelet transform on the sphere. Focusing on the inner $60^{\circ} \times 60^{\circ}$ of the sky, we identify and characterize four separate residuals beyond the expected Milky Way diffuse emission. We detect the \textit{Fermi} Bubbles, finding compelling evidence that they are diffuse in nature and contain very little small-scale structure. We detect the "cocoon" inside the Southern Bubble, and we also identify its northern counterpart above 2 GeV. The Northern Cocoon lies along the same axis but is $\sim 30 \%$ dimmer than the southern one. We characterize the Galactic center excess, which we find extends up to $20^{\circ}$ in $|b|$. At latitudes $|b| \leq 5^{\circ}$ we find evidence for power in small angular scales that could be the result of point-source contributions, but for $|b| \geq 5^{\circ}$ the Galactic center excess is dominantly diffuse in its nature. Our findings show that either the Galactic center excess and {\it Fermi} Bubbles connect smoothly or that the Bubbles brighten significantly below $15^\circ$ in latitude. We find that the Galactic center excess appears off-center by a few degrees towards negative $\ell$. Additionally, we find and characterize two emissions along the Galactic disk centered at $\ell \simeq +25^{\circ}$ and $-20^{\circ}$. These emissions are significantly more elongated along the Galactic disk than the Galactic center excess.

There is growing evidence that intermediate-mass black holes (IMBHs), defined here as having a mass in the range M=500-10^5 Msun, are present in the dense centers of certain globular clusters (GCs). Gravitational waves (GWs) from their mergers with other IMBHs or with stellar BHs in the cluster are mostly emitted in frequencies <10 Hz, which unfortunately is out of reach for current ground-based observatories such as advanced LIGO (aLIGO). Nevertheless, we show that aLIGO measurements can be used to efficiently probe one of the possible formation mechanisms of IMBHs in GCs, namely a runaway merger process of stellar seed BHs. In this case, aLIGO will be sensitive to the lower-mass rungs of the merger ladder, ranging from the seed BH mass to masses >~50-300 Msun, where the background from standard mergers is expected to be very low. Assuming this generic IMBH formation scenario, we calculate the mass functions that correspond to the limiting cases of possible merger trees. Based on estimates for the number density of GCs and taking into account the instrumental sensitivity, we show that current observations do not effectively limit the occupation fraction f_occ of IMBHs formed by runaway mergers of stellar BHs in GCs. However, we find that a six-year run of aLIGO at design sensitivity will be able to probe down to f_occ<3% at a 99.9% confidence level, either finding evidence for this formation mechanism, or necessitating others if the fraction of GCs that harbor IMBHs is higher.

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 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.

The detection of gravitational waves from the merger of binary black holes by the LIGO Collaboration has opened a new window to astrophysics. With the sensitivities of ground based detectors in the coming years we can only detect the local black hole binary mergers. The integrated merger rate can instead be probed by the gravitational-wave background, the incoherent superposition of the released energy in gravitational waves during binary-black-hole coalescence. Through that, the properties of the binary black holes can be studied. In this work we show that by measuring the energy density $\Omega_{GW}$ (in units of the cosmic critical density) of the gravitational-wave background, we can search for the rare $\sim 100 M_{\odot}$ massive black holes formed in the Universe. In addition, we can answer how often the least massive BHs of mass $> 3 M_{\odot}$ form. Finally, if there are multiple channels for the formation of binary black holes and if any of them predicts a narrow mass range for the black holes, then the gravitational-wave-background spectrum may have features that with the future Einstein Telescope can be detected.

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.

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.

We examine how future gravitational-wave measurements from merging black holes (BHs) can be used to infer the shape of the black-hole mass function, with important implications for the study of star formation and evolution and the properties of binary BHs. We model the mass function as a power law, inherited from the stellar initial mass function, and introduce lower and upper mass cutoff parameterizations in order to probe the minimum and maximum BH masses allowed by stellar evolution, respectively. We initially focus on the heavier BH in each binary, to minimize model dependence. Taking into account the experimental noise, the mass measurement errors and the uncertainty in the redshift-dependence of the merger rate, we show that the mass function parameters, as well as the total rate of merger events, can be measured to <10% accuracy within a few years of advanced LIGO observations at its design sensitivity. This can be used to address important open questions such as the upper limit on the stellar mass which allows for BH formation and to confirm or refute the currently observed mass gap between neutron stars and BHs. In order to glean information on the progenitors of the merging BH binaries, we then advocate the study of the two-dimensional mass distribution to constrain parameters that describe the two-body system, such as the mass ratio between the two BHs, in addition to the merger rate and mass function parameters. We argue that several years of data collection can efficiently probe models of binary formation, and show, as an example, that the hypothesis that some gravitational-wave events may involve primordial black holes can be tested. Finally, we point out that in order to maximize the constraining power of the data, it may be worthwhile to lower the signal-to-noise threshold imposed on each candidate event and amass a larger statistical ensemble of BH mergers.

Recent Advanced LIGO detections of binary black hole mergers have prompted multiple studies investigating the possibility that the heavy GW150914 binary system was of primordial origin, and hence could be evidence for dark matter in the form of black holes. We compute the stochastic background arising from the incoherent superposition of such primordial binary black hole systems in the universe and compare it to the similar background spectrum due to binary black hole systems of stellar origin. We investigate the possibility of detecting this background with future gravitational wave detectors, and discuss the possibility of using the stochastic gravitational-wave background measurement to constrain the dark matter component in the form of black holes.

It was recently suggested that the merger of $\sim30\,M_\odot$ primordial black holes (PBHs) may provide a significant number of events in gravitational-wave observatories over the next decade, if they make up an appreciable fraction of the dark matter. Here we show that measurement of the eccentricities of the inspiralling binary black holes can be used to distinguish these binaries from those produced by more traditional astrophysical mechanisms. These PBH binaries are formed on highly eccentric orbits and can then merge on timescales that in some cases are years or less, retaining some eccentricity in the last seconds before the merger. This is to be contrasted with massive-stellar-binary, globular-cluster, or other astrophysical origins for binary black holes (BBHs) in which the orbits have very effectively circularized by the time the BBH enters the observable LIGO window. Here we discuss the features of the gravitational-wave signals that indicate this eccentricity and forecast the sensitivity of LIGO and the Einstein Telescope to such effects. We show that if PBHs make up the dark matter, then roughly one event should have a detectable eccentricity given LIGO's expected sensitivity and observing time of six years. The Einstein Telescope should see $O(10)$ such events after ten years.

We consider the possibility that the black-hole (BH) binary detected by LIGO may be a signature of dark matter. Interestingly enough, there remains a window for masses $20\,M_\odot \lesssim M_{\rm bh} \lesssim 100\, M_\odot$ where primordial black holes (PBHs) may constitute the dark matter. If two BHs in a galactic halo pass sufficiently close, they radiate enough energy in gravitational waves to become gravitationally bound. The bound BHs will rapidly spiral inward due to emission of gravitational radiation and ultimately merge. Uncertainties in the rate for such events arise from our imprecise knowledge of the phase-space structure of galactic halos on the smallest scales. Still, reasonable estimates span a range that overlaps the $2-53$ Gpc$^{-3}$ yr$^{-1}$ rate estimated from GW150914, thus raising the possibility that LIGO has detected PBH dark matter. PBH mergers are likely to be distributed spatially more like dark matter than luminous matter and have no optical nor neutrino counterparts. They may be distinguished from mergers of BHs from more traditional astrophysical sources through the observed mass spectrum, their high ellipticities, or their stochastic gravitational wave background. Next generation experiments will be invaluable in performing these tests.

We investigate a possible method for determining the progenitors of black hole (BH) mergers observed via their gravitational wave (GW) signal. We argue that measurements of the cross-correlation of the GW events with overlapping galaxy catalogs may provide an additional tool in determining if BH mergers trace the stellar mass of the Universe, as would be expected from mergers of the endpoints of stellar evolution. If on the other hand the BHs are of primordial origin, as has been recently suggested, their merging would be preferentially hosted by lower biased objects, and thus have a lower cross-correlation with luminous galaxies. Here we forecast the expected precision of the cross-correlation measurement for current and future GW detectors such as LIGO and the Einstein Telescope. We then predict how well these instruments can distinguish the model that identifies high-mass BH-BH mergers as the merger of primordial black holes that constitute the dark matter in the Universe from more traditional astrophysical sources.

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.

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.

We demonstrate how the image analysis technique of wavelet decomposition can be applied to the gamma-ray sky to separate emission on different angular scales. New structures on scales that differ from the scales of the conventional astrophysical foreground and background uncertainties can be robustly extracted, allowing a model-independent characterization with no presumption of exact signal morphology. As a test case, we generate mock gamma-ray data to demonstrate our ability to extract extended signals without assuming a fixed spatial template. For some point source luminosity functions, our technique also allows us to differentiate a diffuse signal in gamma-rays from dark matter annihilation and extended gamma-ray point source populations in a data-driven way.

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.

Several groups have identified an extended excess of gamma rays over the modeled foreground and background emissions towards the Galactic center (GC) based on observations with the Fermi Large Area Telescope. This excess emission is compatible in morphology and spectrum with a telltale sign from dark matter (DM) annihilation. Here, we present a critical reassessment of DM interpretations of the GC signal in light of the foreground and background uncertainties that some of us recently outlaid in Calore et al. 2014. We find that a much larger number of DM models fits the gamma-ray data than previously noted. In particular: (1) In the case of DM annihilation into $\bar{b}b$, we find that even large DM masses up to $m_\chi \simeq$ 74 GeV are allowed with a $p$-value $> 0.05$. (2) Surprisingly, annihilation into non-relativistic hh gives a good fit to the data. (3) The inverse Compton emission from $\mu^+\mu^-$ with $m_\chi\sim$ 60-70 GeV can also account for the excess at higher latitudes, $|b|>2^\circ$, both in its spectrum and morphology. We also present novel constraints on a large number of mixed annihilation channels, including cascade annihilation involving hidden sector mediators. Finally, we show that the current limits from dwarf spheroidal observations are not in tension with a DM interpretation when uncertainties on the DM halo profile are accounted for.

Recently, a spatially extended excess of gamma rays collected by the Fermi-LAT from the inner region of the Milky Way has been detected by different groups and with increasingly sophisticated techniques. Yet, any final conclusion about the morphology and spectral properties of such an extended diffuse emission are subject to a number of potentially critical uncertainties, related to the high density of cosmic rays, gas, magnetic fields and abundance of point sources. We will present a thorough study of the systematic uncertainties related to the modelling of diffuse background and to the propagation of cosmic rays in the inner part of our Galaxy. We will test a large set of models for the Galactic diffuse emission, generated by varying the propagation parameters within extreme conditions. By using those models in the analysis of Fermi-LAT data as Galactic foreground, we will show that the gamma-ray excess survives and we will quantify the uncertainties affecting the excess morphology and energy spectrum.

The possible gamma-ray excess in the inner Galaxy and the Galactic center (GC) suggested by Fermi-LAT observations has triggered a large number of studies. It has been interpreted as a variety of different phenomena such as a signal from WIMP dark matter annihilation, gamma-ray emission from a population of millisecond pulsars, or emission from cosmic rays injected in a sequence of burst-like events or continuously at the GC. We present the first comprehensive study of model systematics coming from the Galactic diffuse emission in the inner part of our Galaxy and their impact on the inferred properties of the excess emission at Galactic latitudes $2^\circ<|b|<20^\circ$ and 300 MeV to 500 GeV. We study both theoretical and empirical model systematics, which we deduce from a large range of Galactic diffuse emission models and a principal component analysis of residuals in numerous test regions along the Galactic plane. We show that the hypothesis of an extended spherical excess emission with a uniform energy spectrum is compatible with the Fermi-LAT data in our region of interest at $95\%$ CL. Assuming that this excess is the extended counterpart of the one seen in the inner few degrees of the Galaxy, we derive a lower limit of $10.0^\circ$ ($95\%$ CL) on its extension away from the GC. We show that, in light of the large correlated uncertainties that affect the subtraction of the Galactic diffuse emission in the relevant regions, the energy spectrum of the excess is equally compatible with both a simple broken power-law of break energy $2.1\pm0.2$ GeV, and with spectra predicted by the self-annihilation of dark matter, implying in the case of $\bar{b}b$ final states a dark matter mass of $49^{+6.4}_{-5.4}$ GeV.

We perform a multichannel analysis of the indirect signals for the Wino Dark Matter, including one-loop electroweak and Sommerfeld enhancement corrections. We derive limits from cosmic ray antiprotons and positrons, from continuum galactic and extragalactic diffuse gamma-ray spectra, from the absence of gamma-ray line features at the galactic center above 500 GeV in energy, from gamma-rays toward nearby dwarf spheroidal galaxies and galaxy clusters, and from CMB power-spectra. Additionally, we show the future prospects for neutrino observations toward the inner Galaxy and from antideuteron searches. For each of these indirect detection probes we include and discuss the relevance of the most important astrophysical uncertainties that can impact the strength of the derived limits. We find that the Wino as a dark matter candidate is excluded in the mass range bellow ~800 GeV from antiprotons and between 1.8 and 3.5 TeV from the absence of a gamma-ray line feature toward the galactic center. Limits from other indirect detection probes confirm the main bulk of the excluded mass ranges.

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).

Recent advances in gamma-ray cosmic ray, infrared and radio astronomy have allowed us to develop a significantly better understanding of the galactic medium properties in the last few years. In this work using the DRAGON code, that numerically solves the CR propagation equation and calculating gamma-ray emissivities in a 2-dimensional grid enclosing the Galaxy, we study in a self consistent manner models for the galactic diffuse gamma-ray emission. Our models are cross-checked to both the available CR and gamma-ray data. We address the extend to which dark matter annihilations in the Galaxy can contribute to the diffuse gamma-ray flux towards different directions on the sky. Moreover we discuss the impact that astrophysical uncertainties of non DM nature, have on the derived gamma-ray limits. Such uncertainties are related to the diffusion properties on the Galaxy, the interstellar gas and the interstellar radiation field energy densities. Light ~10 GeV dark matter annihilating dominantly to hadrons is more strongly constrained by gamma-ray observations towards the inner parts of the Galaxy and influenced the most by assumptions of the gas distribution; while TeV scale DM annihilating dominantly to leptons has its tightest constraints from observations towards the galactic center avoiding the galactic disk plane, with the main astrophysical uncertainty being the radiation field energy density. In addition, we present a method of deriving constraints on the dark matter distribution profile from the diffuse gamma-ray spectra. These results critically depend on the assumed mass of the dark matter particles and the type of its end annihilation products.

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.

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.

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 recent uncovering of the \textit{Fermi} Bubbles/haze in the \textit{Fermi} gamma-ray data has generated theoretical work to explain such a signal of hard $\gamma$-rays in combination with the \textit{WMAP} haze signal. Many of these theoretical models can have distinctively different implications with regards to the production of high energy neutrinos. We discuss the neutrino signals from different models proposed for the explanation of the \textit{Fermi} Bubbles/haze, more explicitly, from Dark Matter annihilation in the galactic halo with conditions of preferential CR diffusion, from recent AGN jet activity, from periodic diffusive shock acceleration, from stochastic 2nd order Fermi acceleration and from long time-scale star formation in the galactic center in combination with strong galactic winds. We find that some of these models will be probed by the IceCube DeepCore detector. Moreover, with a km$^3$ telescope located at the north hemisphere, we will be able to discriminate between the hadronic, leptonic and the DM models. Additionally using the reconstructed neutrino spectra we will probe annihilation of TeV scale dark matter towards the galactic center.

Recently a series of indications have been put forward suggesting the presence of two gamma-ray lines at 110-130 GeV (centered at 111 and 129 GeV). Signals of these lines have been observed toward the Galactic center, at some galaxy clusters and among some of the unassociated point sources of the 2 years Fermi catalogue. Such a combination of signals could be generated by dark matter annihilations in the main dark matter halo, its substructures and nearby galaxy clusters. We discuss in this work the consistency between the number of events observed at the line energies in the sky and the predictions using results from the Via Lactea II numerical simulation and extrapolations below its mass resolution, taking into account that the annihilation cross-section to the lines can be estimated from the Galactic center signal. We find that some extrapolations to small substructures can naturally account for the point sources signal, although the hypothesis of background only cannot be rejected. We also study the morphology of the gamma-ray sky at the 2 lines energies, testing different Galactic diffuse background models to account for interstellar medium uncertainties, and different assumptions on the DM diffuse component profile. We find from template fits that within reasonable diffuse background uncertainties the presence of a spherical halo component is preferred with cuspier dark matter halo profiles being preferable even from the full sky fit. We finally check the impact of a dark disk component, suggested by cosmological simulations that include baryons and find that thin dark disks can not be disfavored, thus possibly accounting for the preferentially closer to the Galactic disk distribution of the point sources lines signal.

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.

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.

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.

Indications for a gamma-ray line(s) signal towards the Galactic center at an energy of about 130 GeV have been recently presented. While dark matter annihilations are a viable candidate for this signal, it is generally expected that such a flux would be correlated to a gamma-ray component with continuum energy spectrum due to dark matter pair annihilating into other Standard Model particles. We use the gamma-ray data from the inner 10 x 10 degree window to derive limits for a variety of DM annihilation final states. Extending the window of observation, we discuss bounds on the morphological shape of a dark matter signal associated to the line, applying both standard templates for the dark matter profile, such as an Einasto or a NFW profile, and introducing a new more general parametrization.

Dwarf spheroidal galaxies compose one of the most dark matter dominated classes of objects, making them a set of targets to search for signals of dark matter annihilation. Recent developments in gamma-ray astronomy, most importantly the launch of the Fermi-LAT instrument, have brought those targets into attention. Yet, no clear excess of gamma-rays has been confirmed from these targets, resulting in some of the tightest limits on dark matter annihilation from indirect searches. In extracting limits from dwarf spheroidal galaxies, it is of great importance, to properly take into account all relevant uncertainties. Those include the dark matter distribution properties of the dwarf spheroidals and the uncertainties on the underlying background. We revisit the limits on dark matter annihilation, from gamma-rays studying a set of close-by dwarf spheroidal galaxies, for which, we have good understanding of the uncertainties in the dark matter distribution. For those targets, we perform and compare results for alternative methods in extracting the background gamma-ray flux. This provides a method to discriminate among the dark matter annihilation targets, those that can give robust constraints. We finally present our tightest limits on dark matter annihilation, that come only from the targets that ensure accurate understanding of both the gamma-ray background and the dark matter distribution uncertainties.

The latest years have seen steady progresses in WIMP dark matter (DM) searches, with hints of possible signals suggested by both direct and indirect detection experiments. Antiprotons can play a key role validating those interpretations since they are copiously produced by WIMP annihilations in the Galactic halo, and the secondary antiproton background produced by Cosmic Ray (CR) interactions is predicted with fair accuracy and matches the observed spectrum very well. Using the publicly available numerical DRAGON code, we reconsider antiprotons as a tool to constrain DM models discussing its power and limitations. We provide updated constraints on a wide class of annihilating DM models by comparing our predictions against the most up-to-date ap measurements, taking also into account the latest spectral information on the p, He and other CR nuclei fluxes. Doing that, we probe carefully the uncertainties associated to both secondary and DM originated antiprotons, by using a variety of distinctively different assumptions for the propagation of CRs and for the DM distribution in the Galaxy. We find that the impact of the astrophysical uncertainties on constraining the DM properties can be much stronger, up to a factor of ~50, than the one due to uncertainties on the DM distribution (~2-6). Remarkably, even reducing the uncertainties on the propagation parameters derived by local observables, non-local effects can still change DM model constraints even by 50%. Nevertheless, current ap data place tight constraints on DM models, excluding some of those suggested in connection with indirect and direct searches. Finally we discuss the power of upcoming CR spectral data from the AMS-02 observatory to drastically reduce the uncertainties discussed in this paper and estimate the expected sensitivity of this instrument to some sets of DM models.

We study the high latitude (|b|>10) diffuse gamma-ray emission in the Galaxy in light of the recently published data from the Fermi collaboration at energies between 100 MeV and 100 GeV. The unprecedented accuracy in these measurements allows to probe and constrain the properties of sources and propagation of cosmic rays (CRs) in the Galaxy, as well as confirming conventional assumptions made on the interstellar medium (ISM). Using the publicly available DRAGON code, that has been shown to reproduce local measurements of CRs, we study assumptions made in the literature on HI and H2 gas distributions in the ISM, and non spatially uniform models of diffusion in the Galaxy. By performing a combined analysis of CR and gamma-ray spectra, we derive constraints on the properties of the ISM gas distribution and the vertical scale height of galactic CR diffusion, which may have implications also on indirect Dark Matter detection. We also discuss some of the possible interpretations of the break at ~230 GV in CR protons and helium spectra, recently observed by PAMELA and their impact on gamma-rays.

The spectral data on the diffuse Galactic gamma-rays, at medium and high latitudes (|b| > 10) and energies of 1-100 GeV, recently published by the Fermi Collaboration are used to produce a novel study on the gamma-ray emissivity in the Galaxy. We focus on analyzing the properties of propagation of cosmic rays (CRs), using the publicly available DRAGON code. We critically address some of the models for the interstellar HI and H2 gas distributions commonly used in the literature, as well as test a variety of propagation models. Each model assumes a distinct global profile for the diffusion and the re-acceleration of CRs. Fitting propagation parameters to well measured local CRs such as, the B/C ratio, protons, antiprotons and electron, positron fluxes, we evaluate the gamma-ray spectra at medium and high latitudes in order to place further constraints on these propagation models.

Recently the PAMELA experiment has released its updated anti-proton flux and anti-proton to proton flux ratio data up to energies of ~200GeV. With no clear excess of cosmic ray anti-protons at high energies, one can extend constraints on the production of anti-protons from dark matter. In this letter, we consider both the cases of dark matter annihilating and decaying into standard model particles that produce significant numbers of anti-protons. We provide two sets of constraints on the annihilation cross-sections/decay lifetimes. In the one set of constraints we ignore any source of anti-protons other than dark matter, which give the highest allowed cross-sections/inverse lifetimes. In the other set we include also anti-protons produced in collisions of cosmic rays with interstellar medium nuclei, getting tighter but more realistic constraints on the annihilation cross-sections/decay lifetimes.

We study possible astrophysical and dark matter (DM) explanations for the Fermi gamma-ray haze in the Milky Way halo. As representatives of various DM models, we consider DM particles annihilating into W+W-, b-bbar, and e+e-. In the first two cases, the prompt gamma-ray emission from DM annihilations is significant or even dominant at E > 10 GeV, while inverse Compton scattering (ICS) from annihilating DM products is insignificant. For the e+e- annihilation mode, we require a boost factor of order 100 to get significant contribution to the gamma-ray haze from ICS photons. Possible astrophysical sources of high energy particles at high latitudes include type Ia supernovae (SNe) and millisecond pulsars (MSPs). Based on our current understanding of Ia SNe rates, they do not contribute significantly to gamma-ray flux in the halo of the Milky Way. As the MSP population in the stellar halo of the Milky Way is not well constrained, MSPs may be a viable source of gamma-rays at high latitudes provided that there are ~ 20 000 - 60 000 of MSPs in the Milky Way stellar halo. In this case, pulsed gamma-ray emission from MSPs can contribute to gamma-rays around few GeV's while the ICS photons from MSP electrons and positrons may be significant at all energies in the gamma-ray haze. The plausibility of such a population of MSPs is discussed. Consistency with the Wilkinson Microwave Anisotropy Probe (WMAP) microwave haze requires that either a significant fraction of MSP spin-down energy is converted into e+e- flux or the DM annihilates predominantly into leptons with a boost factor of order 100.

We argue that the decomposition of gamma-ray maps in spherical harmonics is a sensitive tool to study dark matter (DM) annihilation or decay in the main Galactic halo of the Milky Way. Using the spherical harmonic decomposition in a window excluding the Galactic plane, we show for one year of Fermi data that adding a spherical template (such as a line-of-sight DM annihilation profile) to an astrophysical background significantly reduces chi^2 of the fit to the data. In some energy bins the significance of this DM-like fraction is above three sigma. This can be viewed as a hint of DM annihilation signal, although astrophysical sources cannot be ruled out at this moment. We use the derived DM fraction as a conservative upper limit on DM annihilation signal. In the case of bb-bar annihilation channel the limits are about a factor of two less constraining than the limits from dwarf galaxies. The uncertainty of our method is dominated by systematics related to modeling the astrophysical background. We show that with one year of Fermi data the statistical sensitivity would be sufficient to detect DM annihilation with thermal freeze out cross section for masses below 100 GeV.

Recent full-sky maps of the Galaxy from the Fermi Gamma-Ray Space Telescope have revealed a diffuse component of emission towards the Galactic center and extending up to roughly +/-50 degrees in latitude. This Fermi "haze" is the inverse Compton emission generated by the same electrons which generate the microwave synchrotron haze at WMAP wavelengths. The gamma-ray haze has two distinct characteristics: the spectrum is significantly harder than emission elsewhere in the Galaxy and the morphology is elongated in latitude with respect to longitude with an axis ratio ~2. If these electrons are generated through annihilations of dark matter particles in the Galactic halo, this morphology is difficult to realize with a standard spherical halo and isotropic cosmic-ray diffusion. However, we show that anisotropic diffusion along ordered magnetic field lines towards the center of the Galaxy coupled with a prolate dark matter halo can easily yield the required morphology without making unrealistic assumptions about diffusion parameters. Furthermore, a Sommerfeld enhancement to the self annihilation cross-section of ~30 yields a good fit to the morphology, amplitude, and spectrum of both the gamma-ray and microwave haze. The model is also consistent with local cosmic-ray measurements as well as CMB constraints.

Much attention has been given to dark matter explanations of the PAMELA positron fraction and Fermi electronic excesses. For those theories with a TeV-scale WIMP annihilating through a light force-carrier, the associated Sommerfeld enhancement provides a natural explanation of the large boost factor needed to explain the signals, and the light force-carrier naturally gives rise to hard cosmic ray spectra without excess pi0 gamma rays or anti-protons. The Sommerfeld enhancement of the annihilation rate, which at low relative velocities v scales as 1/v, relies on the comparatively low velocity dispersion of the dark matter particles in the smooth halo. Dark matter substructures in which the velocity dispersion is smaller than in the smooth halo have even larger annihilation rates. N-body simulations containing only dark matter predict the existence of such structures, for example subhalos and caustics, and the effects of these substructures on dark matter indirect detection signals have been studied extensively. The addition of baryons into cosmological simulations of disk-dominated galaxies gives rise to an additional substructure component, a dark disk. The disk has a lower velocity dispersion than the spherical halo component by a factor ~6, so the contributions to dark matter signals from the disk can be more significant in Sommerfeld models than for WIMPs without such low-velocity enhancements. We consider the consequences of a dark disk on the observed signals of cosmic rays as measured by Fermi and PAMELA in models where the WIMP annihilations are into a light boson. We find that both the PAMELA and Fermi results are easily accomodated by scenarios in which a disk signal is included with the standard spherical halo signal. Limits from extrapolations to the center of the galaxy can also be modified.

The Fermi Gamma-Ray Space Telescope reveals a diffuse inverse Compton signal in the inner Galaxy with a similar spatial morphology to the microwave haze observed by WMAP, supporting the synchrotron interpretation of the microwave signal. Using spatial templates, we regress out pi0 gammas, as well as IC and bremsstrahlung components associated with known soft-synchrotron counterparts. We find a significant gamma-ray excess towards the Galactic center with a spectrum that is significantly harder than other sky components and is most consistent with IC from a hard population of electrons. The morphology and spectrum are consistent with it being the IC counterpart to the electrons which generate the microwave haze seen at WMAP frequencies. In addition, the implied electron spectrum is hard; electrons accelerated in supernova shocks in the disk which then diffuse a few kpc to the haze region would have a softer spectrum. We describe the full sky Fermi maps used in this analysis and make them available for download.

Recent data from cosmic ray experiments such as PAMELA, Fermi, ATIC and PPB-BETS all suggest the need for a new primary source of electrons and positrons at high (>~100 GeV) energies. Many proposals have been put forth to explain these data, usually relying on a single particle to annihilate or decay to produce e+e-. In this paper, we consider models with multiple species of WIMPs with significantly different masses. We show if such dark matter candidates chi_i annihilate into light bosons, they naturally produce equal annihilation rates, even as the available numbers of pairs for annihilation n_chi_i^2 differ by orders of magnitude. We argue that a consequence of these models can be to add additional signal naturally at lower (~100 GeV) versus higher (~ TeV) energies, changing the expected spectrum and even adding bumps at lower energies, which may alleviate some of the tension in the required annihilation rates between PAMELA and Fermi. These spectral changes may yield observable consequences in the microwave Haze signal observed at the upcoming Planck satellite. Such a model can connect to other observable signals such as DAMA and INTEGRAL by having the lighter (heavier) state be a pseudo-Dirac fermion with splitting 100 keV (1 MeV). We show that variations in the halo velocity dispersion can alleviate constraints from final state radiation in the galactic center and galactic ridge. If the lighter WIMP has a large self-interaction cross section, the light-WIMP halo might collapse, dramatically altering expectations for direct and indirect detection signatures.

In this paper, we study the flux of electrons and positrons injected by pulsars and by annihilating or decaying dark matter in the context of recent ATIC, PAMELA, Fermi, and HESS data. We review the flux from a single pulsar and derive the flux from a distribution of pulsars. We point out that the particle acceleration in the pulsar magnetosphere is insufficient to explain the observed excess of electrons and positrons with energy E ~ 1 TeV and one has to take into account an additional acceleration of electrons at the termination shock between the pulsar and its wind nebula. We show that at energies less than a few hundred GeV, the flux from a continuous distribution of pulsars provides a good approximation to the expected flux from pulsars in the Australia Telescope National Facility (ATNF) catalog. At higher energies, we demonstrate that the electron/positron flux measured at the Earth will be dominated by a few young nearby pulsars, and therefore the spectrum would contain bumplike features. We argue that the presence of such features at high energies would strongly suggest a pulsar origin of the anomalous contribution to electron and positron fluxes. The absence of features either points to a dark matter origin or constrains pulsar models in such a way that the fluctuations are suppressed. Also we derive that the features can be partially smeared due to spatial variation of the energy losses during propagation.

Multiple lines of evidence indicate an anomalous injection of high-energy e+- in the Galactic halo. The recent $e^+$ fraction spectrum from the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) shows a sharp rise up to 100 GeV. The Fermi Gamma-ray Space Telescope has found a significant hardening of the e+e- cosmic ray spectrum above 100 GeV, with a break, confirmed by HESS at around 1 TeV. The Advanced Thin Ionization Calorimeter (ATIC) has also detected detected a similar excess, falling back to the expected spectrum at 1 TeV and above. Excess microwaves towards the galactic center in the WMAP data are consistent with hard synchrotron radiation from a population of 10-100 GeV e+- (the WMAP ``Haze''). We argue that dark matter annihilations can provide a consistent explanation of all of these data, focusing on dominantly leptonic modes, either directly or through a new light boson. Normalizing the signal to the highest energy evidence (Fermi and HESS), we find that similar cross sections provide good fits to PAMELA and the Haze, and that both the required cross section and annihilation modes are achievable in models with Sommerfeld-enhanced annihilation. These models naturally predict significant production of gamma rays in the galactic center via a variety of mechanisms. Most notably, there is a robust inverse-Compton scattered (ICS) gamma-ray signal arising from the energetic electrons and positrons, detectable at Fermi/GLAST energies, which should provide smoking gun evidence for this production.

Recently, a preliminary spectrum from the Fermi Gamma-ray Space Telescope has been presented for the inner galaxy (-30 < l < 30, -5 < b < 5), as well as the galactic center (-1 < l < 1, -1< b < 1). We consider the implications of these data for dark matter annihilation models, especially models capable of producing the cosmic-ray excesses previously observed by PAMELA and Fermi. These local cosmic-ray excesses, when extrapolated to the inner galaxy, imply inverse Compton scattering (ICS) gamma-ray signals largely consistent with the preliminary Fermi gamma-ray spectrum. For specific halos and models, particularly those with prompt photons, the data have begun to constrain the allowed parameter space. Although significant modeling and background uncertainties remain, we explore how large a signal is permitted by the current data. Based upon this, we make predictions for what signal could be present in other regions of the sky where dark matter signals may be easier to isolate from the astrophysical backgrounds.

Recently published results from the PAMELA experiment have shown conclusive evidence for an excess of positrons at high (~ 10 - 100 GeV) energies, confirming earlier indications from HEAT and AMS-01. Such a signal is generally expected from dark matter annihilations. However, the hard positron spectrum and large amplitude are difficult to achieve in most conventional WIMP models. The absence of any associated excess in anti-protons is highly constraining on any model with hadronic annihilation modes. We revisit an earlier proposal, whereby the dark matter annihilates into a new light (<~GeV) boson phi, which is kinematically constrained to go to hard leptonic states, without anti-protons or pi0's. We find this provides a very good fit to the data. The light boson naturally provides a mechanism by which large cross sections can be achieved through the Sommerfeld enhancement, as was recently proposed. Depending on the mass of the WIMP, the rise may continue above 300 GeV, the extent of PAMELA's ability to discriminate electrons and positrons.

We consider the signals of positrons and electrons from "exciting" dark matter (XDM) annihilation. Because of the light (m_phi ~< 1 GeV) force carrier phi into which the dark matter states can annihilate, the electrons and positrons are generally very boosted, yielding a hard spectrum, in addition to the low energy positrons needed for INTEGRAL observations of the galactic center. We consider the relevance of this scenario for HEAT, PAMELA and the WMAP "haze," focusing on light (m_phi ~< 2 m_pi) phi bosons, and find that significant signals can be found for all three, although significant signals generally require high dark matter densities. We find that measurements of the positron fraction are generally insensitive to the halo model, but do suffer significant astrophysical uncertainties. We discuss the implications for upcoming PAMELA results.