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 observations by the Alpha Magnetic Spectrometer (AMS-02) have tentatively detected a handful of cosmic-ray antihelium events. Such events have long been considered as smoking-gun evidence for new physics, because astrophysical antihelium production is expected to be negligible. However, the dark-matter-induced antihelium flux is also expected to fall below current sensitivities, particularly in light of existing antiproton constraints. Here, we demonstrate that a previously neglected standard model process -- the production of antihelium through the displaced-vertex decay of $\bar{\Lambda}_b$-baryons -- can significantly boost the dark matter induced antihelium flux. This process can triple the standard prompt-production of antihelium, and more importantly, entirely dominate the production of the high-energy antihelium nuclei reported by AMS-02.

The observed correlation between the far-infrared (FIR) and radio luminosities of star-forming galaxies shows the close connection between star formation and cosmic-ray production. LOFAR recently extended radio observations of the related correlation with star formation rate to lower frequencies (150~MHz), finding a peculiar radio excess in galaxies with high stellar masses and low star-formation rates. We show that recycled/millisecond pulsars (MSPs) can dominate the non-thermal emission in these massive quiescent galaxies and explain the excess. This is supported by recent data suggesting that MSPs can efficiently convert a large fraction of their spin-down power to e$^+$e$^-$ pairs. We find that MSP-based models provide a significantly improved (formally $>19\sigma$, though systematic errors dominate) fit to LOFAR data. We discuss implications for the radio-FIR correlation, the observation of radio excesses in nearby galaxies, and local electron and positron observations.

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 Universe's early thermal history is poorly constrained, and it is possible that it underwent a period of early matter domination driven by a heavy particle or an oscillating scalar field that decayed into radiation before the onset of Big Bang nucleosynthesis. The entropy sourced by this particle's decay reduces the cross section required for thermal-relic dark matter to achieve the observed abundance. This degeneracy between dark matter properties and the thermal history vastly widens the field of viable dark matter candidates, undermining efforts to constrain dark matter's identity. Fortunately, an early matter-dominated era also amplifies density fluctuations at small scales and leads to early microhalo formation, boosting the dark matter annihilation rate and bringing smaller cross sections into the view of existing indirect-detection probes. Employing several recently developed models of microhalo formation and evolution, we develop a procedure to derive indirect-detection constraints on dark matter annihilation in cosmologies with early matter domination. This procedure properly accounts for the unique morphology of microhalo-dominated signals. While constraints depend on dark matter's free-streaming scale, the microhalos make it possible to obtain upper bounds as small as $\langle\sigma v\rangle \lesssim 10^{-32}$ cm$^3$s$^{-1}$ using Fermi-LAT observations of the isotropic gamma-ray background and the Draco dwarf galaxy.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

Milagro and HAWC have detected extended TeV gamma-ray emission around nearby pulsar wind nebulae (PWNe). Building on these discoveries, Linden et al. [1] identified a new source class -- TeV halos -- powered by the interactions of high-energy electrons and positrons that have escaped from the PWN, but which remain trapped in a larger region where diffusion is inhibited compared to the interstellar medium. Many theoretical properties of TeV halos remain mysterious, but empirical arguments suggest that they are ubiquitous. The key to progress is finding more halos. We outline prospects for new discoveries and calculate their expectations and uncertainties. We predict, using models normalized to current data, that future HAWC and CTA observations will detect in total $\sim$50--240 TeV halos, though we note that multiple systematic uncertainties still exist. Further, the existing HESS source catalog could contain $\sim$10--50 TeV halos that are presently classified as unidentified sources or PWN candidates. We quantify the importance of these detections for new probes of the evolution of TeV halos, pulsar properties, and the sources of high-energy gamma rays and cosmic rays.

Cosmic rays (CRs) and magnetic fields may be dynamically important in driving large-scale galactic outflows from rapidly star-forming galaxies. We construct two-dimensional axisymmetric models of the local starburst and super-wind galaxy M82 using the CR propagation code GALPROP. Using prescribed gas density and magnetic field distributions, wind profiles, CR injection rates, and stellar radiation fields, we simultaneously fit both the integrated gamma-ray emission and the spatially-resolved multi-frequency radio emission extended along M82's minor axis. We explore the resulting constraints on the gas density, magnetic field strength, CR energy density, and the assumed CR advection profile. In accord with earlier one-zone studies, we generically find low central CR pressures, strong secondary electron/positron production, and an important role for relativistic bremsstrahlung losses in shaping the synchrotron spectrum. We find that the relatively low central CR density produces CR pressure gradients that are weak compared to gravity, strongly limiting the role of CRs in driving M82's fast and mass-loaded galactic outflow. Our models require strong magnetic fields and advection speeds of order ~1000 km/s on kpc scales along the minor axis in order to reproduce the extended radio emission. Degeneracies between the controlling physical parameters of the model and caveats to these findings are discussed.

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.

Fermi-LAT observations have strongly constrained dark matter annihilation through the joint-likelihood analysis of dwarf spheroidal galaxies (dSphs). These constraints are expected to be robust because dSphs have measurable dark matter content and produce negligible astrophysical emission. However, each dSph is dim, with a predicted flux that typically falls below the accuracy of the background model. We show that this significantly diminishes the reliability of previous joint-likelihood algorithms, and develop an improved analysis that directly accounts for the effect of background mismodeling. This method produces more robust limits and detections of dark matter in both real and mock data. We calculate improved limits on the dark matter annihilation cross-section, which differ by nearly a factor of two from previous analyses - despite examining identical data.

The Sun is an excellent laboratory for astroparticle physics but remains poorly understood at GeV--TeV energies. Despite the immense relevance for both cosmic-ray propagation and dark matter searches, only in recent years has the Sun become a target for precision gamma-ray astronomy with the Fermi-LAT instrument. Among the most surprising results from the observations is a hard excess of GeV gamma-ray flux that strongly anti-correlates with solar activity, especially at the highest energies accessible to Fermi-LAT. Most of the observed properties of the gamma-ray emission cannot be explained by existing models of cosmic-ray interactions with the solar atmosphere. GeV--TeV gamma-ray observations of the Sun spanning an entire solar cycle would provide key insights into the origin of these gamma rays, and consequently improve our understanding of the Sun's environment as well as the foregrounds for new physics searches, such as dark matter. These can be complemented with new observations with neutrinos and cosmic rays. Together these observations make the Sun a new testing ground for particle physics in dynamic environments.

TeV Halos, extended regions of TeV gamma-ray emission around middle-aged pulsars, have recently been established as a new source class in gamma-ray astronomy. These halos have been attributed to relativistic electrons and positrons that have left the acceleration region close to the pulsar and are diffusing in the surrounding medium. Measuring the morphology of TeV Halos enables for the first time a direct measurement of the electron diffusion on scales of tens of parsecs. There are hints that the presence of relativistic particles affects the diffusion rate in the pulsars' surroundings. Understanding electron diffusion is necessary to constrain the origins of the apparent `excess' of cosmic-ray positrons at tens of GeV. TeV Halos can also be used to find mis-aligned pulsars, as well as study certain properties of the Galaxy's pulsar population. Future VHE gamma-ray instruments will detect more of those TeV Halos and determine how much pulsars contribute to the observed cosmic-ray electron and positron fluxes, and how they affect diffusion in their environments.

The era of precision cosmology has revealed that about 85% of the matter in the universe is dark matter. Two well-motivated candidates are weakly interacting massive particles (WIMPs) and weakly interacting sub-eV particles (WISPs) (e.g. axions). Both WIMPs and WISPs possess distinct {\gamma}-ray signatures. Over the last decade, data taken between 50 MeV to >300 GeV by the Fermi Large Area Telescope (Fermi-LAT) have provided stringent constraints on both classes of dark matter models. Thus far, there are no conclusive detections. However, there is an intriguing {\gamma}-ray excess associated with the Galactic center that could be explained by WIMP annihilation. At lower energies, the poor angular resolution of the Fermi-LAT makes source identification challenging, inhibiting our ability to more sensitively probe both the Galactic center excess, as well as lower-mass WIMP and WISP models. Additionally, targeted WISP searches (e.g., those probing supernovae and blazars) would greatly benefit from enhanced energy resolution and polarization measurements in the MeV range. To address these issues, a new telescope that is optimized for MeV observations is needed. Such an instrument would allow us to explore new areas of dark matter parameter space and provide unprecedented access to its particle nature.

The 511 keV line from positron annihilation in the Galaxy was the first $\gamma$-ray line detected to originate from outside our solar system. Going into the fifth decade since the discovery, the source of positrons is still unconfirmed and remains one of the enduring mysteries in $\gamma$-ray astronomy. With a large flux of $\sim$10$^{-3}$ $\gamma$/cm$^{2}$/s, after 15 years in operation INTEGRAL/SPI has detected the 511 keV line at $>50\sigma$ and has performed high-resolution spectral studies which conclude that Galactic positrons predominantly annihilate at low energies in warm phases of the interstellar medium. The results from imaging are less certain, but show a spatial distribution with a strong concentration in the center of the Galaxy. The observed emission from the Galactic disk has low surface brightness and the scale height is poorly constrained, therefore, the shear number of annihilating positrons in our Galaxy is still not well know. Positrons produced in $\beta^+$-decay of nucleosynthesis products, such as $^{26}$Al, can account for some of the annihilation emission in the disk, but the observed spatial distribution, in particular the excess in the Galactic bulge, remains difficult to explain. Additionally, one of the largest uncertainties in these studies is the unknown distance that positrons propagate before annihilation. In this paper, we will summarize the current knowledge base of Galactic positrons, and discuss how next-generation instruments could finally provide the answers.

We analyze the Sun as a source for the indirect detection of dark matter through a search for gamma rays from the solar disk. Capture of dark matter by elastic interactions with the solar nuclei followed by annihilation to long-lived mediators can produce a detectable gamma-ray flux. We search three years of data from the High Altitude Water Cherenkov Observatory and find no statistically significant detection of TeV gamma-ray emission from the Sun. Using this, we constrain the spin-dependent elastic scattering cross section of dark matter with protons for dark matter masses above 1 TeV, assuming an unstable mediator with a favorable lifetime. The results complement constraints obtained from Fermi-LAT observations of the Sun and together cover WIMP masses between 4 GeV and $10^6$ GeV. The cross section constraints for mediator decays to gamma rays can be as strong as $\sim10^{-45}$ cm$^{-2}$, which is more than four orders of magnitude stronger than current direct-detection experiments for 1 TeV dark matter mass. The cross-section constraints at higher masses are even better, nearly 7 orders of magnitude better than the current direct-detection constraints for 100 TeV dark matter mass. This demonstration of sensitivity encourages detailed development of theoretical models in light of these powerful new constraints.

Steady gamma-ray emission up to at least 200 GeV has been detected from the solar disk in the Fermi-LAT data, with the brightest, hardest emission occurring during solar minimum. The likely cause is hadronic cosmic rays undergoing collisions in the Sun's atmosphere after being redirected from ingoing to outgoing in magnetic fields, though the exact mechanism is not understood. An important new test of the gamma-ray production mechanism will follow from observations at higher energies. Only the High Altitude Water Cherenkov (HAWC) Observatory has the required sensitivity to effectively probe the Sun in the TeV range. Using three years of HAWC data from November 2014 to December 2017, just prior to the solar minimum, we search for 1--100 TeV gamma rays from the solar disk. No evidence of a signal is observed, and we set strong upper limits on the flux at a few $10^{-12}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV. Our limit, which is the most constraining result on TeV gamma rays from the Sun, is $\sim10\%$ of the theoretical maximum flux (based on a model where all incoming cosmic rays produce outgoing photons), which in turn is comparable to the Fermi-LAT data near 100 GeV. The prospects for a first TeV detection of the Sun by HAWC are especially high during solar minimum, which began in early 2018.

The solar disk is a bright source of multi-GeV gamma rays, due to the interactions of hadronic cosmic rays with the solar atmosphere. However, the underlying production mechanism is not understood, except that its efficiency must be greatly enhanced by magnetic fields that redirect some cosmic rays from ingoing to outgoing before they interact. To elucidate the nature of this emission, we perform a new analysis of solar atmospheric gamma rays with 9 years of Fermi-LAT data, which spans nearly the full 11-year solar cycle. We detect significant gamma-ray emission from the solar disk from 1 GeV up to $\gtrsim200$ GeV. The overall gamma-ray spectrum is much harder ($\sim E_{\gamma}^{-2.2}$) than the cosmic-ray spectrum ($\sim E_{\rm CR}^{-2.7}$). We find a clear anticorrelation between the solar cycle phase and the gamma-ray flux between 1-10 GeV. Surprisingly, we observe a spectral dip between $\sim$30-50 GeV in an otherwise power-law spectrum. This was not predicted, is not understood, and may provide crucial clues to the gamma-ray emission mechanism. The flux above 100 GeV, which is brightest during the solar minimum, poses exciting opportunities for HAWC, LHAASO, IceCube, and KM3NeT.

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.

Recent observations have detected extended TeV gamma-ray emission surrounding young and middle-aged pulsars. The morphology of these "TeV halos" requires cosmic-ray diffusion to be locally suppressed by a factor of ~100-1000 compared to the typical interstellar medium. No model currently explains this suppression. We show that cosmic-ray self-confinement can significantly inhibit diffusion near pulsars. The steep cosmic-ray gradient generates Alfven waves that resonantly scatter the same cosmic-ray population, suppressing diffusion within ~20 pc of pulsars younger than ~100 kyr. In this model, TeV halos evolve through two phases, a growth phase where Alfven waves are resonantly generated and cosmic-ray diffusion becomes increasingly suppressed, and a subsequent relaxation phase where the diffusion coefficient returns to the standard interstellar value. Intriguingly, cosmic-rays are not strongly confined early in the TeV halo evolution, allowing a significant fraction of injected e+e- to escape. If these e+e- also escape from the surrounding supernova remnant, they would provide a natural explanation for the positron excess observed by PAMELA and AMS-02. Recently created TeV cosmic-rays are confined in the TeV halo, matching observations by HAWC and H.E.S.S. While our default model relaxes too rapidly to explain the confinement of TeV cosmic rays around mature pulsars, such as Geminga, models utilizing a Kraichnan turbulence spectrum experience much slower relaxation. Thus, observations of TeV halos around mature pulsars may provide a probe into our understanding of interstellar turbulence.

Galaxy formation simulations demonstrate that cosmic-ray (CR) feedback may be important in the launching of galactic-scale winds. CR protons dominate the bulk of the CR population, yet most observational constraints of CR feedback come from synchrotron emission of CR electrons. In this paper, we present an analysis of 105 months of Fermi Gamma-ray Space Telescope observations of the Small Magellanic Cloud (SMC), with the aim of exploring CR feedback and transport in an external galaxy. We produce maps of the 2-300 GeV emission and detect statistically significant, extended emission along the Bar and the Wing, where active star formation is occurring. Gamma-ray emission is not detected above 13 GeV, and we set stringent upper-limits on the flux above this energy. We find the best fit to the gamma-ray spectrum is a single-component model with a power-law of index $\Gamma=-2.11\pm0.06\pm0.06$ and an exponential cutoff energy of $E_{\rm c} =13.1\pm5.1\pm1.6$ GeV. We assess the relative contribution of pulsars and CRs to the emission, and we find that pulsars may produce up to 14$^{+4}_{-2}$% of the flux above 100 MeV. Thus, we attribute most of the gamma-ray emission (based on its spectrum and morphology) to CR interactions with the ISM. We show that the gamma-ray emissivity of the SMC is five times smaller than that of the Milky Way and that the SMC is far below the calorimetric limit, where all CR protons experience pion losses. We interpret these findings as evidence that CRs are escaping the SMC via advection and diffusion.

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.

We identify new astrophysical signatures of dark matter that implodes neutron stars (NSs), which could decisively test whether NS-imploding dark matter is responsible for missing pulsars in the Milky Way galactic center, the source of some $r$-process elements, and the origin of fast-radio bursts. First, NS-imploding dark matter forms $\sim 10^{-10}$ solar mass or smaller black holes inside neutron stars, which proceed to convert neutron stars into $\sim$1.5 solar mass BHs. This decreases the number of neutron star mergers seen by LIGO/Virgo (LV) and associated merger kilonovae seen by telescopes like DES, BlackGEM, and ZTF, and instead, producing a population of "black mergers" containing $\sim$1.5 solar mass black holes. Second, dark matter-induced neutron star implosions may create a new kind of kilonovae that lacks a detectable, accompanying gravitational signal, which we call "quiet kilonovae." Using DES data and the Milky Way's r-process abundance, we constrain quiet kilonovae. Third, the spatial distribution of neutron star merger kilonovae and quiet kilonovae in galaxies can be used to detect dark matter. NS-imploding dark matter destroys most neutron stars at the centers of disc galaxies, so that neutron star merger kilonovae would appear mostly in a donut at large radii. We find that as few as ten neutron star merger kilonova events, located to $\sim$1 kpc precision could validate or exclude dark matter-induced neutron star implosions at $2 \sigma$ confidence, exploring dark matter-nucleon cross-sections 4-10 orders of magnitude below current direct detection experimental limits. Similarly, NS-imploding dark matter as the source of fast radio bursts can be tested at $2 \sigma$ confidence once 20 bursts are located in host galaxies by radio arrays like CHIME and HIRAX.

The observed multi-GeV gamma-ray emission from the solar disk --- sourced by hadronic cosmic rays interacting with gas, and affected by complex magnetic fields --- is not understood. Utilizing an improved analysis of the Fermi-LAT data that includes the first resolved imaging of the disk, we find strong evidence that this emission is produced by two separate mechanisms. Between 2010-2017 (the rise to and fall from solar maximum), the gamma-ray emission is dominated by a polar component. Between 2008-2009 (solar minimum) this component remains present, but the total emission is instead dominated by a new equatorial component with a brighter flux and harder spectrum. Most strikingly, although 6 gamma rays above 100 GeV are observed during the 1.4 years of solar minimum, none are observed during the next 7.8 years. These features, along with a 30-50 GeV spectral dip which will be discussed in a companion paper, were not anticipated by theory. To understand the underlying physics, Fermi and HAWC observations of the imminent Cycle 25 solar minimum are crucial.

We summarize the radio synchrotron background workshop that took place July 19-21, 2017 at the University of Richmond. This first scientific meeting dedicated to the topic was convened because current measurements of the diffuse radio monopole reveal a surface brightness that is several times higher than can be straightforwardly explained by known Galactic and extragalactic sources and processes, rendering it by far the least well understood photon background at present. It was the conclusion of a majority of the participants that the radio monopole level is at or near that reported by the ARCADE 2 experiment and inferred from several absolutely calibrated zero level lower frequency radio measurements, and unanimously agreed that the production of this level of surface brightness, if confirmed, represents a major outstanding question in astrophysics. The workshop reached a consensus on the next priorities for investigations of the radio synchrotron background.

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

The nature of the $\gamma$-ray emission from the \emph{Fermi} bubbles is unknown. Both hadronic and leptonic models have been formulated to explain the peculiar $\gamma$-ray signal observed by the Fermi-LAT between 0.1-500~GeV. If this emission continues above $\sim$30~TeV, hadronic models of the \emph{Fermi} bubbles would provide a significant contribution to the high-energy neutrino flux detected by the IceCube observatory. Even in models where leptonic $\gamma$-rays produce the \emph{Fermi} bubbles flux at GeV energies, a hadronic component may be observable at very high energies. The combination of IceCube and HAWC measurements have the ability to distinguish these scenarios through a comparison of the neutrino and $\gamma$-ray fluxes at a similar energy scale. We examine the most recent four-year dataset produced by the IceCube collaboration and find no evidence for neutrino emission originating from the \emph{Fermi} bubbles. In particular, we find that previously suggested excesses are consistent with the diffuse astrophysical background with a p-value of 0.22 (0.05 in an extreme scenario that all the IceCube events that overlap with the bubbles come from them). Moreover, we show that existing and upcoming HAWC observations provide independent constraints on any neutrino emission from the \emph{Fermi} bubbles, due to the close correlation between the $\gamma$-ray and neutrino fluxes in hadronic interactions. The combination of these results disfavors a significant contribution from the \emph{Fermi} bubbles to the IceCube neutrino flux.

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.

Milagro observations have found bright, diffuse TeV emission concentrated along the galactic plane of the Milky Way. The intensity and spectrum of this emission is difficult to explain with current models where gamma-ray production is dominated by hadronic mechanisms, and has been named the "TeV excess". We show that TeV emission from pulsars naturally explains this excess. In particular, recent observations have detected "TeV halos" surrounding pulsars that are either nearby or particularly luminous. Here, we show that the full population of Milky Way pulsars will produce diffuse TeV emission concentrated along the Milky Way plane. The total gamma-ray flux from TeV halos is expected to exceed the hadronic gamma-ray flux at energies above ~500 GeV. Moreover, the spectrum and intensity of TeV halo emission naturally matches the TeV excess. If this scenario is common to all galaxies, it will decrease the contribution of star-forming galaxies to the IceCube neutrino flux. Finally, we show that upcoming HAWC observations will resolve a significant fraction of the TeV excess into individual TeV halos, conclusively confirming, or ruling out, this model.

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.

We identify a largely model-independent signature of dark matter interactions with nucleons and electrons. Dark matter in the local galactic halo, gravitationally accelerated to over half the speed of light, scatters against and deposits kinetic energy into neutron stars, heating them to infrared blackbody temperatures. The resulting radiation could potentially be detected by the James Webb Space Telescope, the Thirty Meter Telescope, or the European Extremely Large Telescope. This mechanism also produces optical emission from neutron stars in the galactic bulge, and X-ray emission near the galactic center, because dark matter is denser in these regions. For GeV - PeV mass dark matter, dark kinetic heating would initially unmask any spin-independent or spin-dependent dark matter-nucleon cross-sections exceeding $2 \times 10^{-45}$ cm$^2$, with improved sensitivity after more telescope exposure. For lighter-than-GeV dark matter, cross-section sensitivity scales inversely with dark matter mass because of Pauli blocking; for heavier-than-PeV dark matter, it scales linearly with mass as a result of needing multiple scatters for capture. Future observations of dark sector-warmed neutron stars could determine whether dark matter annihilates in or only kinetically heats neutron stars. Because inelastic inter-state transitions of up to a few GeV would occur in relativistic scattering against nucleons, elusive inelastic dark matter like pure Higgsinos can also be discovered.

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.

Observations by the Fermi-LAT have uncovered a bright, spherically symmetric excess surrounding the center of the Milky Way galaxy. The spectrum of the gamma-ray excess peaks sharply at an energy ~2 GeV, exhibiting a hard spectrum at lower energies, and falls off quickly above an energy ~5 GeV. The spectrum of the excess above ~10 GeV is potentially an important discriminator between different physical models for its origin. We focus our study on observations of the gamma-ray excess at energies exceeding 10 GeV, finding: (1) a statistically significant excess remains in the energy range 9.5-47.5 GeV, which is not degenerate with known diffuse emission templates such as the Fermi Bubbles, (2) the radial profile of the excess at high energies remains relatively consistent with data near the spectral peak, (3) the data above ~5 GeV prefer a slightly greater ellipticity with a major axis oriented perpendicular to the Galactic plane. Using the recently developed non-Poissonian template fit, we find mild evidence for a point-source origin for the high-energy excess, although given the statistical and systematic uncertainties we show that a smooth origin of the high-energy emission cannot be ruled out. We discuss the implication of these findings for pulsar and dark matter models of the gamma-ray excess. Finally we provide a number of updated measurements of the gamma-ray excess, utilizing novel diffuse templates and the Pass 8 dataset.

The isotropic gamma-ray background arises from the contribution of unresolved sources, including members of confirmed source classes and proposed gamma-ray emitters such as the radiation induced by dark matter annihilation and decay. Clues about the properties of the contributing sources are imprinted in the anisotropy characteristics of the gamma-ray background. We use 81 months of Pass 7 Reprocessed data from the Fermi Large Area Telescope to perform a measurement of the anisotropy angular power spectrum of the gamma-ray background. We analyze energies between 0.5 and 500 GeV, extending the range considered in the previous measurement based on 22 months of data. We also compute, for the first time, the cross-correlation angular power spectrum between different energy bins. We find that the derived angular spectra are compatible with being Poissonian, i.e. constant in multipole. Moreover, the energy dependence of the anisotropy suggests that the signal is due to two populations of sources, contributing, respectively, below and above 2 GeV. Finally, using data from state-of-the-art numerical simulations to model the dark matter distribution, we constrain the contribution from dark matter annihilation and decay in Galactic and extragalactic structures to the measured anisotropy. These constraints are competitive with those that can be derived from the average intensity of the isotropic gamma-ray background. Data are available at https://www-glast.stanford.edu/pub_data/552.

The origin of the isotropic gamma-ray background (IGRB) --- the portion of the extragalactic gamma-ray sky that is not resolvable into individual point sources --- provides a powerful probe into the evolution of the high-energy universe. Star-forming galaxies (SFGs) are among the most likely contributors to the IGRB, though their contribution is difficult to constrain because their flux distribution is dominated by numerous faint sources. We produce a novel joint-likelihood analysis of the $\gamma$-ray emission from 584 SFGs, utilizing advanced statistical techniques to compare the distribution of low-significance excesses against the non-Poissonian $\gamma$-ray background fluctuations. We first examine the theoretically well-motivated relationship between the far-IR and $\gamma$-ray luminosities of SFGs, utilizing a model where the $\gamma$-ray luminosity is given by log$_{10}$(L$_\gamma$/(erg s$^{-1}$)) = $\alpha$ log$_{10}$(L$_{IR}$/(10$^{10}$L$_\odot$)) + $\beta$. We calculate best-fit parameters $\alpha$ = 1.18 $\pm$ 0.15, $\beta$ = 38.49 $\pm$ 0.24, with a log-normal dispersion in this relationship given by $\sigma$ = 0.39 $\pm$ 0.12. The best-fit values of $\alpha$ and $\beta$ are consistent with previous studies. We find a larger dispersion in the far-IR to $\gamma$-ray correlation than previous studies. This dispersion is significant at the level of 5.7$\sigma$. These results imply that SFGs significantly contribute to the IGRB, producing between 61.0$^{+30.2}_{-18.3}\%$ of the total IGRB intensity above an energy of 1 GeV. Along with recent works, this strongly indicates that multiple source classes provide comparable contributions to the IGRB intensity. We discuss the implication of these results for the interpretation of the IceCube neutrinos.

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.

Recent observations of Reticulum II have uncovered an overabundance of r-process elements, compared to similar ultra-faint dwarf spheroidal galaxies (UFDs). Because the metallicity and star formation history of Reticulum II appear consistent with all known UFDs, the high r-process abundance of Reticulum II suggests enrichment through a single, rare event, such as a double neutron star (NS) merger. However, we note that this scenario is extremely unlikely, as binary stellar evolution models require significant supernova natal kicks to produce NS-NS or NS-black hole mergers, and these kicks would efficiently remove compact binary systems from the weak gravitational potentials of UFDs. We examine alternative mechanisms for the production of r-process elements in UFDs, including a novel mechanism wherein NSs in regions of high dark matter density implode after accumulating a black-hole-forming mass of dark matter. We find that r-process proto-material ejection by tidal forces, when a single neutron star implodes into a black hole, can occur at a rate matching the r-process abundance of both Reticulum II and the Milky Way. Remarkably, dark matter models which collapse a single neutron star in observed UFDs also solve the missing pulsar problem in the Milky Way Galactic center. We propose tests specific to dark matter r-process production which may uncover, or rule out, this model.

Radio observations at multiple frequencies have detected a significant isotropic emission component between 22~MHz and 10~GHz, commonly termed the ARCADE-2 Excess. The origin of this radio emission is unknown, as the intensity, spectrum and isotropy of the signal are difficult to model with either traditional astrophysical mechanisms or novel physics such as dark matter annihilation. We posit a new model capable of explaining the key components of the excess radio emission. Specifically, we show that the re-acceleration of non-thermal electrons via turbulence in merging galaxy clusters are capable of explaining the intensity, spectrum, and isotropy of the ARCADE-2 data. We examine the parameter spaces of cluster re-acceleration, magnetic field, and merger rate, finding that the radio excess can be reproduced assuming reasonable assumptions for each. Finally, we point out that future observations will definitively confirm or rule-out the contribution of cluster mergers to the isotropic radio background.

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.

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.

Recent Fermi-LAT observations of dwarf spheroidal galaxies in the Milky Way have placed strong limits on the gamma-ray flux from dark matter annihilation. In order to produce the strongest limit on the dark matter annihilation cross-section, the observations of each dwarf galaxy have typically been "stacked" in a joint-likelihood analysis, utilizing optical observations to constrain the dark matter density profile in each dwarf. These limits have typically been computed only for singular annihilation final states, such as $b\bar{b}$ or $\tau^+\tau^-$. In this paper, we generalize this approach by producing an independent joint-likelihood analysis to set constraints on models where the dark matter particle annihilates to multiple final state fermions. We interpret these results in the context of the most popular simplified models, including those with s- and t-channel dark matter annihilation through scalar and vector mediators. We present our results as constraints on the minimum dark matter mass and the mediator sector parameters. Additionally, we compare our simplified model results to those of Effective Field Theory contact interactions in the high-mass limit.

Fermi-LAT observations of the Galactic Center (GC) have revealed a spherically- symmetric excess of GeV gamma rays extending to at least 10 deg from the dynamical center of the Galaxy. A critical uncertainty in extracting the intensity, spectrum, and morphology of this excess concerns the accuracy of astrophysical diffuse gamma-ray emission models near the GC. Recently, it has been noted that many diffuse emission models utilize a cosmic-ray injection rate far below that predicted based on the observed star formation rate in the Central Molecular Zone. In this study, we add a cosmic-ray injection component which non-linearly traces the Galactic H2 density determined in three-dimensions, and find that the associated gamma-ray emission is degenerate with many properties of the GC gamma-ray excess. In models that utilize a large sideband (40x40 deg surrounding the GC) to normalize the best-fitting diffuse emission models, the intensity of the GC excess decreases by a factor of 2, and the morphology of the excess becomes less peaked and less spherically symmetric. In models which utilize a smaller region of interest (15x15 deg) the addition of an excess template instead suppresses the intensity of the best-fit astrophysical diffuse emission, and the GC excess is rather resilient to changes in the details of the astrophysical diffuse modeling. In both analyses, the addition of a GC excess template still provides a statistically significant improvement to the overall fit to the gamma-ray data. We also implement advective winds at the GC, and find that the Fermi-LAT data strongly prefer outflows of order several hundred km/s, whose role is to efficiently advect low-energy cosmic rays from the Galactic center. Finally, we perform numerous tests of our models, and conclude that they significantly improve our understanding of multi-wavelength non-thermal emission from the GC.

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

At present, all physical models of diffuse Galactic gamma-ray emission assume that the distribution of cosmic-ray sources traces the observed populations of either OB stars, pulsars, or supernova remnants. However, since H2-rich regions host significant star formation and numerous supernova remnants, the morphology of observed H2 gas should also provide a physically motivated, high-resolution tracer for cosmic-ray injection. We assess the impact of utilizing H2 as a tracer for cosmic-ray injection on models of diffuse Galactic gamma-ray emission. We employ state-of-the-art 3D particle diffusion and gas density models, along with a physical model for the star-formation rate based on global Schmidt laws. Allowing a fraction, f_H2, of cosmic-ray sources to trace the observed H2 density, we find that a theoretically well-motivated value f_H2 ~ 0.20 -- 0.25 (i) provides a significantly better global fit to the diffuse Galactic gamma-ray sky and (ii) highly suppresses the intensity of the residual gamma-ray emission from the Galactic center region. Specifically, in models utilizing our best global fit values of f_H2 ~ 0.20 -- 0.25, the spectrum of the galactic center gamma-ray excess is drastically affected, and the morphology of the excess becomes inconsistent with predictions for dark matter annihilation.

Observations using the Fermi Large Area Telescope (Fermi-LAT) have found a significant gamma-ray excess surrounding the center of the Milky Way (GC). One possible interpretation of this excess invokes gamma-ray emission from an undiscovered population of either young or recycled pulsars densely clustered throughout the inner kiloparsec of the Milky Way. While these systems, by construction, have individual fluxes that lie below the point source sensitivity of the Fermi-LAT, they may already be observed in multiwavelength observations. Notably the Australia Telescope National Facility (ATNF) catalog of radio pulsars includes 270 sources observed in the inner 10 degrees around the GC. We calculate the gamma-ray emission observed from these 270 sources and obtain three key results: (1) point source searches in the GC region produce a plethora of highly significant gamma-ray "hotspots", compared to searches far from the Galactic plane, (2) there is no statistical correlation between the positions of these gamma-ray hotspots and the locations of ATNF pulsars, and (3) the spectrum of the most statistically significant gamma-ray hotspots is substantially softer than the spectrum of the GC gamma-ray excess. These results place strong constraints on models where young pulsars produce the majority of the gamma-ray excess, and disfavor some models where millisecond pulsars produce the gamma-ray excess.

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.

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.

The recent discovery of a neutron star accretor in the ultra-luminous X-ray source M82 X-2 challenges our understanding of high-mass X-ray binary formation and evolution. By combining binary population synthesis and detailed mass-transfer models, however, we show that the binary parameters of M82 X-2 are not surprising provided non-conservative mass transfer is allowed. Specifically, the donor-mass lower limit and orbital period measured for M82 X-2 lie near the most probable values predicted by population synthesis models, and systems such as M82 X-2 should exist in approximately 13% of the galaxies with a star-formation history similar to M82. We conclude that the binary system that formed M82 X-2 is most likely less than 50 Myr old and contains a donor star which had an initial mass of approximately 8-10 M$_\odot$, while the NS's progenitor star had an initial mass in the $8-25\,\rm M_{\odot}$ range. The donor star still currently resides on the main sequence, and is capable of continued MT on the thermal timescale, while in the ultra-luminous X-ray regime, for as long as 400,000 years.

Observations by the Fermi-LAT telescope have uncovered a significant $\gamma$-ray excess toward the Milky Way Galactic Center. There has been no detection of a similar signal in the direction of the Milky Way dwarf spheroidal galaxies. Additionally, astronomical observations indicate that dwarf galaxies and other faint galaxies are less dense than predicted by the simplest cold dark matter models. We show that a self-interacting dark matter model with a particle mass of roughly 50 GeV annihilating to the mediator responsible for the strong self-interaction can simultaneously explain all three observations. The mediator is necessarily unstable and its mass must be below about 100 MeV in order to lower densities in faint galaxies. If the mediator decays to electron-positron pairs with a cross section on the order of the thermal relic value, then we find that these pairs can up-scatter the interstellar radiation field and produce the observed $\gamma$-ray excess. We show that this model is compatible with all current constraints and highlight detectable signatures unique to self-interacting dark matter models.

Observations of the extragalactic radio background have uncovered a significant isotropic emission across multiple frequencies spanning from 22 MHz to 10 GHz. The intensity of this non-thermal emission component significantly exceeds the expected contribution from known astrophysical sources. Interestingly, models have indicated that the annihilation of dark matter particles may reproduce both the flux and spectrum of the excess. However, the lack of a measurable anisotropy in the residual emission remains challenging for both dark matter and standard astrophysical interpretations of the ARCADE-2 data. We calculate the expected synchrotron anisotropy from dark matter annihilation and show that these models can produce very small anisotropies, though this requires galaxy clusters to have large substructure contributions and strong magnetic fields. We show that this constraint can be significantly relaxed, however, in scenarios where electrons produced via dark matter annihilation can be efficiently reaccelerated by Alfv\'en waves in the intra-Cluster medium. Our analysis indicates that any source capable of explaining the intensity and isotropy of the extragalactic radio excess must have a spatial extension far larger than typical for baryons in galaxies, suggesting a novel physics interpretation.

The paucity of old millisecond pulsars observed at the galactic center of the Milky Way could be the result of dark matter accumulating in and destroying neutron stars. In regions of high dark matter density, dark matter clumped in a pulsar can exceed the Schwarzschild limit and collapse into a natal black hole which destroys the pulsar. We examine what dark matter models are consistent with this hypothesis and find regions of parameter space where dark matter accumulation can significantly degrade the neutron star population within the galactic center while remaining consistent with observations of old millisecond pulsars in globular clusters and near the solar position. We identify what dark matter couplings and masses might cause a young pulsar at the galactic center to unexpectedly extinguish. Finally, we find that pulsar collapse age scales inversely with the dark matter density and linearly with the dark matter velocity dispersion. This implies that maximum pulsar age is spatially dependent on position within the dark matter halo of the Milky Way. In turn, this pulsar age spatial dependence will be dark matter model dependent.

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 Telescope Array (TA) has observed a statistically significant excess in cosmic-rays with energies above 57 EeV in a region of approximately 1150 square degrees centered on coordinates (R.A. = 146.7, Dec. = 43.2). We note that the location of this excess correlates with two of the 28 extraterrestrial neutrinos recently observed by IceCube. The overlap between the two IceCube neutrinos and the TA excess is statistically significant at the 2$\sigma$ level. Furthermore, the spectrum and intensity of the IceCube neutrinos is consistent with a single source which would also produce the TA excess. Finally, we discuss possible source classes with the correct characteristics to explain the cosmic-ray and neutrino fluxes with a single source.

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 observations suggest that some high-velocity clouds may be confined by massive dark matter halos. In particular, the proximity and proposed dark matter content of the Smith Cloud make it a tempting target for the indirect detection of dark matter annihilation. We argue that the Smith Cloud may be a better target than some Milky Way dwarf spheroidal satellite galaxies and use gamma-ray observations from the Fermi Large Area Telescope to search for a dark matter annihilation signal. No significant gamma-ray excess is found coincident with the Smith Cloud, and we set strong limits on the dark matter annihilation cross section assuming a spatially-extended dark matter profile consistent with dynamical modeling of the Smith Cloud. Notably, these limits exclude the canonical thermal relic cross section ($\sim 3\times10^{-26}{\rm cm}^{3}{\rm s}^{-1}$) for dark matter masses $\lesssim 30$ GeV annihilating via the $b \bar b$ or $\tau^{+}\tau^{-}$ channels for certain assumptions of the dark matter density profile; however, uncertainties in the dark matter content of the Smith Cloud may significantly weaken these constraints.

Self-interacting dark matter (SIDM) models have been proposed to solve the small-scale issues with the collisionless cold dark matter (CDM) paradigm. We derive equilibrium solutions in these SIDM models for the dark matter halo density profile including the gravitational potential of both baryons and dark matter. Self-interactions drive dark matter to be isothermal and this ties the core sizes and shapes of dark matter halos to the spatial distribution of the stars, a radical departure from previous expectations and from CDM predictions. Compared to predictions of SIDM-only simulations, the core sizes are smaller and the core densities are higher, with the largest effects in baryon-dominated galaxies. As an example, we find a core size around 0.5 kpc for dark matter in the Milky Way, more than an order of magnitude smaller than the core size from SIDM-only simulations, which has important implications for indirect searches of SIDM candidates.

Recent observations by the PAMELA and AMS-02 telescopes have uncovered an anomalous rise in the positron fraction at energies above 10~GeV. One possible explanation for this excess is the production of primary electron/positron pairs through electromagnetic cascades in pulsar magnetospheres. This process results in a high multiplicity of electron/positron pairs within the wind-termination shock of pulsar wind nebula (PWN). A consequence of this scenario is that no circular polarization should be observed within PWN, since the contributions from electrons and positrons exactly cancel. Here we note that current radio instruments are capable of setting meaningful limits on the circular polarization of synchrotron radiation in PWNs, which observationally test the model for pulsar production of the local positron excess. The observation of a PWN with detectable circular polarization would cast strong doubt on pulsar interpretations of the positron excess, while observations setting strong limits on the circular polarization of PWN would lend credence to these models. Finally, we indicate which pulsar wind nebulae are likely to provide the best targets for observational tests of the AMS-02 excess.

Cosmic-ray anti-nuclei provide a promising discovery channel for the indirect detection of particle dark matter. Hadron showers produced by the pair-annihilation or decay of Galactic dark matter generate anti-nucleons which can in turn form light anti-nuclei. Previous studies have only focused on the spectrum and flux of low energy antideuterons which, although very rarely, are occasionally also produced by cosmic-ray spallation. Heavier elements ($A\geq3$) have instead entirely negligible astrophysical background and a primary yield from dark matter which could be detectable by future experiments. Using a Monte Carlo event generator and an event-by-event phase space analysis, we compute, for the first time, the production spectrum of \antihe and \antiT for dark matter annihilating or decaying to $b\bar{b}$ and ${W^+}{W^-}$ final states. We then employ a semi-analytic model of interstellar and heliospheric propagation to calculate the \antihe flux as well as to provide tools to relate the anti-helium spectrum corresponding to an arbitrary antideuteron spectrum. Finally, we discuss prospects for current and future experiments, including GAPS and AMS-02.

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.

Recent observations by PAMELA, Fermi-LAT, and AMS-02 have conclusively indicated a rise in the cosmic-ray positron fraction above 10 GeV, a feature which is impossible to mimic under the paradigm of secondary positron production with self-consistent Galactic cosmic-ray propagation models. A leading explanation for the rising positron fraction is an additional source of electron-positron pairs, for example one or more mature, energetic, and relatively nearby pulsars. We point out that any one of two well-known nearby pulsars, Geminga and Monogem, can satisfactorily provide enough positrons to reproduce AMS-02 observations. A smoking-gun signature of this scenario is an anisotropy in the arrival direction of the cosmic-ray electrons and positrons, which may be detectable by existing, or future, telescopes. The predicted anisotropy level is, at present, consistent with limits from Fermi-LAT and AMS-02. We argue that the large collecting area of Atmospheric Cherenkov Telescopes (ACTs) makes them optimal tools for detecting such an anisotropy. Specifically, we show that much of the proton and gamma-ray background, which affects measurements of the cosmic-ray electron-positron spectrum with ACTs, may be controlled in the search for anisotropies. We conclude that observations using archival ACT data could already constrain or substantiate the pulsar origin of the positron anomaly, while upcoming instruments (such as the Cherenkov Telescope Array) will provide strong constraints on the source of the rising positron fraction.

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.

Recent observations indicating the existence of a monochromatic gamma-ray line with energy ~130 GeV in the Fermi-LAT data have attracted great interest due to the possibility that the line feature stems from the annihilation of dark matter particles. Many studies examining the robustness of the putative line-signal have concentrated on its spectral attributes. Here, we study the morphological features of the gamma-ray line photons, which can be used to differentiate a putative dark matter signal from astrophysical backgrounds or instrumental artifacts. Photons stemming from dark matter annihilation will produce events tracing a specific morphology, with a statistical clustering that can be calculated based on models of the dark matter density profile in the inner Galaxy. We apply the DBSCAN clustering algorithm to Fermi gamma-ray data, and show that we can rule out the possibility that 1 (2, 4) or fewer point-like sources produce the observed morphology for the line photons at a 99% (95%, 90%) confidence level. Our study strongly disfavors the main astrophysical background envisioned to produce a line feature at energies above 100 GeV: cold pulsar winds. It is highly unlikely that 4 or more such objects have exactly the same monochromatic cosmic-ray energy needed to produce a gamma-ray line, to within instrumental energy resolution. Furthermore, we show that the larger photon statistics expected with Air Cherenkov Telescopes such as H.E.S.S.-II will allow for extraordinarily stringent morphological tests of the origin of the "line photons".

The detailed origin of the diffuse gamma-ray background is still unknown. However, the contribution of unresolved sources is expected to induce small-scale anisotropies in this emission, which may provide a way to identify and constrain the properties of its contributors. Recent studies have predicted the contributions to the angular power spectrum (APS) from extragalactic and galactic dark matter (DM) annihilation or decay. The Fermi-LAT collaboration reported detection of angular power with a significance larger than $3\sigma$ in the energy range from 1 GeV to 10 GeV on 22 months of data [Ackermann et al. 2012]. For these preliminary results the already published Fermi-LAT APS measurements [Ackermann et al. 2012] are compared to the accurate predictions for DM anisotropies from state-of-the-art cosmological simulations as presented in [Fornasa et al. 2013] to derive constraints on different DM candidates.

Ultra-Luminous X-ray sources (ULX) are X-ray binaries with Lx>1E^39 ergs/s. The most spectacular examples of ULX occur in starburst galaxies and are now understood to be young, luminous High Mass X-ray Binaries. The conditions under which ULX form are poorly understood, but recent evidence suggests they may be more common in low metallicity systems. Here we investigate the hypothesis that ULX form preferentially in low metallicity galaxies by searching for ULX in a sample of Extremely Metal Poor Galaxies (XMPG) observed with the Chandra X-ray Observatory. XMPG are defined as galaxies with log(O/H)+12<7.65, or less than 5% solar. These are the most metal-deficient galaxies known, and a logical place to find ULX if they favor metal poor systems. We compare the number of ULX (corrected for background contamination) per unit of star formation (Nulx) in the XMPG sample with Nulx in a comparison sample of galaxies with higher metallicities taken from the Spitzer Infrared Galaxy Sample (SINGS). We find that ULX occur preferentially in the metal poor sample with a formal statistical significance of 2.3 sigma. We do not see strong evidence for a trend in the formation of ULX in the high metallicity sample: above 12+log(O/H) ~ 8.0 the efficiency of ULX production appears to be flat. The effect we see is strongest in the lowest metallicity bin. We discuss briefly the implications of these results for the formation of black holes in low metallicity gas.

Observations from multiple gamma-ray telescopes have uncovered a high energy gamma-ray source spatially coincident with the Galactic center. Recently, a compelling model for the broad-band gamma-ray emission has been formulated which posits that high energy protons emanating from Sgr A* could produce gamma-rays through pion decays resulting from inelastic collisions with the traversed interstellar gas in the region. Models of the gas distribution in the Galactic center region imply that the resulting gamma-ray morphology would be observed as a point source with all current telescopes, but that the upcoming Cherenkov Telescope Array (CTA) may be able to detect an extended emission profile with an unmistakable morphology. Here, we critically evaluate this claim, employing a three dimensional gas distribution model and a detailed Monte Carlo simulation, and using the anticipated effective area and angular resolution of CTA. We find that the impressive angular resolution of CTA will be key to test hadronic emission models conclusively against, for example, point source or dark matter annihilation scenarios. We comment on the relevance of this result for searches for dark matter annihilation in the Galactic center region.

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.

Recently, tentative evidence for an excess of gamma rays at energies around 130 GeV has been reported from analyses of data from the Fermi Large Area Telescope (LAT). The excess is potentially of great interest, as it could be associated with the pair-annihilation of Galactic dark matter and the subsequent production of monochromatic or internal bremsstrahlung gamma rays. The 130 GeV excess appears when an optimized selection of the target region of interest is employed, a procedure that depends upon the assumed dark matter density profile. For the profiles producing an appreciable signal, these target regions vastly overlap with the region corresponding to the so-called "Fermi bubbles". We argue that the tentative evidence for a line feature is likely due to hard photons in the Fermi bubbles regions, where the gamma-ray spectrum contains a spectral break in the energy range of interest (100-150 GeV). Although the origin of this broken power-law is unclear, it is probably related to standard astrophysical processes and not to dark matter annihilation. A broken power-law provides as good a fit as a line "excess", even within small energy windows, and a significantly better fit for large energy windows.

Recently, detections of a high-energy gamma-ray source at the position of the Galactic center have been reported by multiple gamma-ray telescopes, spanning the energy range between 100 MeV and 100 TeV. Analysis of these signals strongly suggests the TeV emission to have a morphology consistent with a point source up to the angular resolution of the HESS telescope (approximately 3 pc), while the point-source nature of the GeV emission is currently unsettled, with indications that it may be spatially extended. In the case that the emission is hadronic and in a steady state, we show that the expected gamma-ray morphology is dominated by the distribution of target gas, rather than by details of cosmic-ray injection and propagation. Specifically, we expect a significant portion of hadronic emission to coincide with the position of the circum-nuclear ring, which resides between 1-3 pc from the Galactic center. We note that the upcoming Cherenkov Telescope Array (CTA) will be able to observe conclusive correlations between the morphology of the TeV gamma-ray source and the observed gas density, convincingly confirming or ruling out a hadronic origin for the gamma-ray emission.

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.

The paucity of known High-Mass X-Ray Binaries (HMXB) with naked He donor stars (hereafter He star) in the Galaxy has been noted over the years as a surprising fact, given the significant number of Galactic HMXBs containing H-rich donors, which are expected to be their progenitors. This contrast has further sharpened in light of recent observations uncovering a preponderance of HMXBs hosting loosely bound Be donors orbiting neutron stars (NS), which would be expected to naturally evolve into He-HMXBs through dynamical mass transfer onto the NS and a common-envelope (CE) phase. Hence, reconciling the large population of Be-HMXBs with the observation of only one He-HMXB can help constrain the dynamics of CE physics. Here, we use detailed stellar structure and evolution models and show that binary mergers of HMXBs during CE events must be common in order to resolve the tension between these observed populations. We find that, quantitatively, this scenario remains consistent with the typically adopted energy parameterization of CE evolution, yielding expected populations which are not at odds with current observations. However, future observations which better constrain the underlying population of loosely bound O/B-NS binaries are likely to place significant constraints on the efficiency of CE ejection.

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/<I>^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background.

The small angular scale fluctuations of the (on large scale) isotropic gamma-ray background (IGRB) carry information about the presence of unresolved source classes. A guaranteed contribution to the IGRB is expected from the unresolved gamma-ray AGN while other extragalactic sources, Galactic gamma-ray source populations and dark matter Galactic and extragalactic structures (and sub-structures) are candidate contributors. The IGRB was measured with unprecedented precision by the Large Area Telescope (LAT) on-board of the Fermi gamma-ray observatory, and these data were used for measuring the IGRB angular power spectrum (APS). Detailed Monte Carlo simulations of Fermi-LAT all-sky observations were performed to provide a reference against which to compare the results obtained for the real data set. The Monte Carlo simulations are also a method for performing those detailed studies of the APS contributions of single source populations, which are required in order to identify the actual IGRB contributors. We present preliminary results of an anisotropy search in the IGRB. At angular scales <2deg (e.g. above multipole 155), angular power above the photon noise level is detected, at energies between 1 and 10 GeV in each energy bin, with statistical significance between 7.2 and 4.1 sigmas. The energy not dependence of the fluctuation anisotropy is pointing to the presence of one or more unclustered source populations, while the energy dependence of the intensity anisotropy is consistent with source populations having average photon index 2.40\pm0.07.

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.

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.

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 investigate the characteristics of young (< 20 Myr) and bright (Lx > 1e36 erg/s) High-Mass X-ray Binaries (HMXBs) and find the population to be strongly metallicity-dependent. We separate the model populations among two distinct formation pathways: (1) systems undergoing active Roche Lobe Overflow (RLO), and (2) wind accretion systems with donors in the (super)giant (SG) stage, which we find to dominate the HMXB population. We find metallicity to primarily affect the number of systems which move through each formation pathway, rather than the observable parameters of systems which move through each individual pathway. We discuss the most important model parameters affecting the HMXB population at both low and high metallicities. Using these results, we show that (1) the population of ultra-luminous X-Ray sources can be consistently described by very bright HMXBs which undergo stable Roche Lobe overflow with mild super-Eddington accretion and (2) the HMXB population of the bright starburst galaxy NGC~1569 is likely dominated by one extremely metal-poor starburst cluster.

A generic prediction in the paradigm of weakly interacting dark matter is the production of relativistic particles from dark matter pair-annihilation in regions of high dark matter density. Ultra-relativistic electrons and positrons produced in the center of the Galaxy by dark matter annihilation should produce a diffuse synchrotron emission. While the spectral shape of the synchrotron dark matter haze depends on the particle model (and secondarily on the galactic magnetic fields), the morphology of the haze depends primarily on (1) the dark matter density distribution, (2) the galactic magnetic field morphology, and (3) the diffusion model for high-energy cosmic-ray leptons. Interestingly, an unidentified excess of microwave radiation with characteristics similar to those predicted by dark matter models has been claimed to exist near the galactic center region in the data reported by the WMAP satellite, and dubbed the "WMAP haze". In this study, we carry out a self-consistent treatment of the variables enumerated above, enforcing constraints from the available data on cosmic rays, radio surveys and diffuse gamma rays. We outline and make predictions for the general morphology and spectral features of a "dark matter haze" and we compare them to the WMAP haze data. We also characterize and study the spectrum and spatial distribution of the inverse Compton emission resulting from the same population of energetic electrons and positrons. We point out that the spectrum and morphology of the radio emission at different frequencies is a powerful diagnostics to test whether a galactic synchrotron haze indeed originates from dark matter annihilation.

Recent claims of a gamma-ray excess in the diffuse galactic emission detected by the Fermi Large Area Telescope with a morphology similar to the WMAP haze were based on the assumption that spatial templates of the interstellar medium (ISM) column density and the 408 Mhz sky are good proxies for neutral pion and inverse Compton (IC) gamma-ray emission, respectively. We identify significant systematic effects in this procedure that can artificially induce an additional diffuse component with a morphology strikingly similar to the claimed gamma-ray haze. To quantitatively illustrate this point we calculate sky-maps of the ratio of the gamma-ray emission from neutral pions to the ISM column density, and of IC to synchrotron emission, using detailed galactic cosmic-ray models and simulations. In the region above and below the galactic center, the ISM template underestimates the gamma-ray emission due to neutral pion decay by approximately 20%. Additionally, the synchrotron template tends to under-estimate the IC emission at low energies (few GeV) and to over-estimate it at higher energies (tens of GeV) by potentially large factors that depend crucially on the assumed magnetic field structure of the Galaxy. The size of the systematic effects we find are comparable to the size of the claimed "Fermi haze" signal. We thus conclude that a detailed model for the galactic diffuse emission is necessary in order to conclusively assess the presence of a gamma-ray excess possibly associated to the WMAP haze morphology.

Presenting population models of high-mass X-ray binaries (HMXBs) formed after bursts of star formation, we investigate the effect of electron-capture supernovae (ECS) of massive ONeMg white dwarfs and the hypothesis that ECS events are associated with typically low supernova kicks imparted to the nascent neutron stars. We identify an interesting ECS bump in the time evolution of HMXB numbers; this bump is caused by significantly increased production of wind-fed HMXBs 20-60 Myr post starburst. The amplitude and age extent of the ECS bump depend on the strength of ECS kicks and the mass range of ECS progenitors. We also find that ECS-HMXBs form through a specific evolutionary channel that is expected to lead to binaries with Be donors in wide orbits. These characteristics, along with their sensitivity to ECS properties, provide us with an intriguing opportunity to probe ECS physics and progenitors through studies of starbursts of different ages. Specifically, the case of the Small Magellanic Cloud, with a significant observed population of Be HMXBs and starburst activity 30-60 Myr ago, arises as a promising laboratory for understanding the role of electron-capture supernovae in neutron star formation.