Unveiling the origin of the coalescing binaries detected via gravitational waves (GW) is challenging, notably if no multi-wavelength counterpart is detected. One important diagnostic tool is the coalescing binary distribution with respect to the large scale structures (LSS) of the universe, which we quantify via the cross-correlation of galaxy catalogs with GW ones. By using both existing and forthcoming galaxy catalogs and using realistic Monte Carlo simulations of GW events, we find that the cross-correlation signal should be marginally detectable in a 10-year data taking of advanced LIGO-Virgo detectors at design sensitivity, at least for binary neutron star mergers. The expected addition of KAGRA and LIGO-India to the GW detector network would allow for a firmer detection of this signal and, in combination with future cosmological surveys, would also permit the detection of cross-correlation for coalescing black holes. Such a measurement may unveil, for instance, a primordial origin of coalescing black holes. To attain this goal, we find that it is crucial to adopt a tomographic approach and to reach a sufficiently accurate localization of GW events. The depth of forthcoming surveys will be fully exploited by third generation GW detectors such as the Einstein Telescope or the Cosmic Explorer, which will allow one to perform precision studies of the coalescing black hole LSS distribution and attain rather advanced model discrimination capabilities.

We present an updated analysis of the gamma-ray flux from the directions of classical dwarf spheroidal galaxies, deriving new constraints on WIMP dark matter (DM) annihilation using a decade of Fermi-LAT data. Among the major novelties, we infer the dwarfs' J-factors by including new observations without imposing any a priori parametric profile for the DM distribution. While statistically compatible with results obtained from more conventional parameterisations, this procedure reduces the theoretical bias imposed on the data. Furthermore, we retain the full data-driven shape of the J-factors' empirical probability distributions when setting limits on DM, without imposing log-normality as is typically done. In conjunction with the data-driven J-factors, we use a new method for estimating the probability distribution function of the astrophysical background at the dwarf position, fully profiling over background uncertainties. We show that, for most classical dwarfs, the background systematic uncertainty dominates over the uncertainty on their J-factors.

Searches for "dark" subhaloes in gamma-ray point-like source catalogues are among promising strategies for indirect dark matter detection. Such a search is nevertheless affected by uncertainties related, on the one hand, to the modelling of the dark matter subhalo distribution in Milky-Way-like galaxies, and, on the other hand, to the sensitivity of gamma-ray instruments to the dark matter subhalo signals. In the present work, we assess the detectability of dark matter subhaloes in Fermi-LAT catalogues, taking into accounts uncertainties associated with the modelling of the Galactic subhalo population. We use four different halo models bracketing a large set of uncertainties. For each model, adopting an accurate detection threshold of the LAT to dark matter subhalo signals and comparing model predictions with the number of unassociated point-sources in Fermi-LAT catalogues, we derive upper limits on the annihilation cross section as a function of dark matter mass. Our results show that, even in the best-case scenario (i.e. DMonly subhalo model), which does not include tidal disruption from baryons, the limits on the dark matter parameter space are less stringent than current gamma-ray limits from dwarf spheroidal galaxies. Comparing the results obtained with the different subhalo models, we find that baryonic effects on the subhalo population are significant and lead to dark matter constraints that are less stringent by a factor of about 2 to 5. This uncertainty comes from the unknown resilience of dark matter subhaloes to tidal disruption.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

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

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

If the gamma-ray excess towards the inner Galaxy (GCE) detected in Fermi-LAT data is due to millisecond pulsars (MSPs), one expects an associated gravitational wave (GW) signal, whose intensity exceeds the disk MSP population emission by an order of magnitude. We compute the expected GW counterpart of the bulge MSP population based on fits of the GCE, and estimate the sensitivity reach of current and future terrestrial GW detectors. The bounds on the average population ellipticity $\varepsilon$ are competitive with the existing ones derived by LIGO/Virgo towards known MSPs. With a 10-year data taking in current LIGO/Virgo configuration, one would detect a signal at the level $\varepsilon\simeq 10^{-7}$, while $\varepsilon\simeq 10^{-8}$ would be attainable with a similar data taking period with a third generation GW detector. This sensitivity should be sufficient for crucial diagnostics on the GCE interpretation in terms of MSPs.

We make a first attempt to find dwarf galaxies in eight \Fermi-LAT extended, unassociated, source fields using \Gaia\ DR2. We probe previously unexplored heliocentric distances of $d<20$~kpc with an extreme-deconvolution (XD) technique. We find no signature of a dwarf galaxy in any of these fields despite \Gaia's excellent astrometric accuracy. We estimate our detection limits by applying the XD method to mock data, obtaining a conservative limit on the stellar mass of $M_* < 10^4$~M$_{\sun}$ for $d < 20$\, kpc. Such a low stellar mass implies either a low-mass subhalo, or a massive stripped-down subhalo. We use an analytic model for stripped subhalos to argue that, given the sizes and fluxes of the \Fermi-LAT sources, we can reject the hypothesis that they owe to dark matter annihilation.

Dwarf galaxies represent a powerful probe of annihilating dark matter particle models, with gamma-ray data setting some of the best bounds available. A major issue in improving over existing constraints consists in the limited knowledge of the astrophysical background (mostly diffuse photons, but also unresolved sources). Perhaps more worrisome, several approaches in the literature suffer of the difficulty of assessing the systematic error due to background mis-modelling. Here we propose a data-driven method to estimate the background at the dwarf position and its uncertainty, relying on an appropriate use of the whole-sky data, via an optimisation procedure of the interpolation weights. While this article is mostly methodologically oriented, we also report the bounds based on latest Fermi-LAT data and updated information for J-factors for both isolated and stacked dwarfs. Our results are very competitive with the Fermi-LAT ones, while being derived with a more general and flexible method. We discuss the impact of profiling over the J-factor as well as over the background probability distribution function, with the latter resulting for instance crucial in drawing conclusions of compatibility with DM interpretations of the so-called Galactic Centre Excess.

Collider, space, and Earth based experiments are now able to probe several extensions of the Standard Model of particle physics which provide viable dark matter candidates. Direct and indirect dark matter searches rely on inputs of astrophysical nature, such as the local dark matter density or the shape of the dark matter profile in the target in object. The determination of these quantities is highly affected by astrophysical uncertainties. The latter, especially those for our own Galaxy, are ill-known, and often not fully accounted for when analyzing the phenomenology of particle physics models. In this paper we present a systematic, quantitative estimate of how astrophysical uncertainties on Galactic quantities (such as the local galactocentric distance, circular velocity, or the morphology of the stellar disk and bulge) propagate to the determination of the phenomenology of particle physics models, thus eventually affecting the determination of new physics parameters. We present results in the context of two specific extensions of the Standard Model (the Singlet Scalar and the Inert Doublet) that we adopt as case studies for their simplicity in illustrating the magnitude and impact of such uncertainties on the parameter space of the particle physics model itself. Our findings point toward very relevant effects of current Galactic uncertainties on the determination of particle physics parameters, and urge a systematic estimate of such uncertainties in more complex scenarios, in order to achieve constraints on the determination of new physics that realistically include all known uncertainties.

The Fermi-LAT Galactic Center excess and the 511 keV positron-annihilation signal from the inner Galaxy bare a striking morphological similarity. We propose that both can be explained through a scenario in which millisecond pulsars produce the Galactic Center excess and their progenitors, low-mass X-ray binaries, the 511 keV signal. As a proof-of-principle we study a specific population synthesis scenario from the literature involving so-called ultracompact X-ray binaries. Moreover, for the first time, we quantitatively show that neutron star, rather than black hole, low-mass X-ray binaries can be responsible for the majority of the positrons. In this particular scenario binary millisecond pulsars can be both the source of the Fermi-LAT $\gamma$-ray excess and the bulge positrons. Future avenues to test this scenario are discussed.

Well-motivated extensions of the standard model predict ultra-light and fundamental pseudo-scalar particles (e.g., axions or axion-like particles: ALPs). Similarly to the Primakoff-effect for axions, ALPs can mix with photons and consequently be searched for in laboratory experiments and with astrophysical observations. Here, we search for energy-dependent modulations of high-energy gamma-ray spectra that are tell-tale signatures of photon-ALPs mixing. To this end, we analyze the data recorded with the Fermi -LAT from Galactic pulsars selected to have a line of sight crossing spiral arms at a large pitch angle. The large-scale Galactic magnetic field traces the shape of spiral arms, such that a sizable photon-ALP conversion probability is expected for the sources considered. In nine years of Fermi -LAT data, we detect significant spectral features in the selected source-sample consistent with photon-ALPs oscillation with a combined statistical significance of 5.52 {\sigma} . Notably, sources with neighboring lines of sight share similar spectral features. From a common fit to all sources, we determine the most-likely parameters for mass $m_{a}$ and coupling $g_{a\gamma\gamma}$ to be $m_a=(3.6 \substack{+0.5_ \mathrm{stat.}\\-0.2_ \mathrm{stat.}}\pm 0.2_\mathrm{syst.} )$ neV and $g_{a\gamma\gamma}=(2.3\substack{+0.3_ \mathrm{stat.}\\-0.4_ \mathrm{stat.}}\pm 0.4_\mathrm{syst.})\times 10^{-10}$ GeV$^{-1}$. In the error budget, we consider instrumental effects, scaling of the adopted Galactic magnetic field model ($\pm~20~\%$), and uncertainties on the distance of individual sources. We note that an astrophysical interpretation of the detected modulation is not obvious.

Stellar-mass Primordial Black Holes (PBH) have been recently reconsidered as a Dark Matter (DM) candidate, after the aLIGO discovery of several binary BH mergers with masses of tens of $M_\odot$. Matter accretion on such massive objects leads to the emission of high-energy photons, capable of altering the ionization and thermal history of the universe. This in turn affects the statistical properties of the cosmic microwave background (CMB) anisotropies. Previous analyses have assumed spherical accretion. We argue that this approximation likely breaks down and that an accretion disk should form in the dark ages. Using the most up-to-date tools to compute the energy deposition in the medium, we derive constraints on the fraction of DM in PBH. Provided that disks form early on, even under conservative assumptions for accretion, these constraints exclude a monochromatic distribution of PBH with masses above $\sim 2\, M_\odot$ as the dominant form of DM. The bound on the median PBH mass gets more stringent if a broad, log-normal mass function is considered. A deepened understanding of non-linear clustering properties and BH accretion disk physics would permit an improved treatment and possibly lead to more stringent constraints.

An anomalous emission component at energies of a few GeV and located towards the inner Galaxy is present in the Fermi-LAT data. It is known as the Fermi-LAT GeV excess. Using almost 8 years of data we reanalyze the characteristics of this excess with SkyFACT, a novel tool that combines image reconstruction with template fitting techniques. We find that an emission profile that traces stellar mass in the boxy and nuclear bulge provides the best description of the excess emission, providing strong circumstantial evidence that the excess is due to a stellar source population in the Galactic bulge. We find a luminosity to stellar mass ratio of $(2.1\pm 0.2)\times 10^{27} \mathrm{\,erg\,s^{-1}\,M_\odot^{-1}}$ for the boxy bulge, and of $(1.4\pm 0.6)\times 10^{27}\mathrm{\,erg\,s^{-1}\,M_\odot^{-1}}$ for the nuclear bulge. Stellar mass related templates are preferred over conventional DM profiles with high statistical significance.

Axion like particles (ALPs) are fundamental pseudo scalar particles with properties similar to Axions which are a well-known extension of the standard model to solve the strong CP problem in Quantum Chromodynamics. ALPs can oscillate into photons and vice versa in the presence of an external tranversal magnetic field. This oscillation of photon and ALPs could have important implications for astronomical observations, i.e. a characteristic energy dependent attenuation in Gamma ray spectra for astrophysical sources. Here we have revisited the opportunity to search Photon-ALPs coupling in the disappearance channel. We use eight years of Fermi Pass 8 data of a selection of promising galactic Gamma-ray source candidates and study the modulation in the spectra in accordance with Photon-ALPs mixing and estimate best fit values of the parameters i.e. Photon-ALPs coupling constant$ (g_{\alpha\gamma\gamma} )$ and ALPs mass($m_{\alpha} $). For the magnetic field we use large scale galactic magnetic field models based on Faraday rotation measurements and we have also studied the survival probability of photons in the Galactic plane.

We present SkyFACT (Sky Factorization with Adaptive Constrained Templates), a new approach for studying, modeling and decomposing diffuse gamma-ray emission. Like most previous analyses, the approach relies on predictions from cosmic-ray propagation codes like GALPROP and DRAGON. However, in contrast to previous approaches, we account for the fact that models are not perfect and allow for a very large number ($\gtrsim10^5$) of nuisance parameters to parameterize these imperfections. We combine methods of image reconstruction and adaptive spatio-spectral template regression in one coherent hybrid approach. To this end, we use penalized Poisson likelihood regression, with regularization functions that are motivated by the maximum entropy method. We introduce methods to efficiently handle the high dimensionality of the convex optimization problem as well as the associated semi-sparse covariance matrix, using the L-BFGS-B algorithm and Cholesky factorization. We test the method both on synthetic data as well as on gamma-ray emission from the inner Galaxy, $|\ell|<90^\circ$ and $|b|<20^\circ$, as observed by the Fermi Large Area Telescope. We finally define a simple reference model that removes most of the residual emission from the inner Galaxy, based on conventional diffuse emission components as well as components for the Fermi bubbles, the Fermi Galactic center excess, and extended sources along the Galactic disk. Variants of this reference model can serve as basis for future studies of diffuse emission in and outside the Galactic disk.

We model the accretion of gas onto a population of massive primordial black holes in the Milky Way, and compare the predicted radio and X-ray emission with observational data. We show that under conservative assumptions on the accretion process, the possibility that ${\cal O}(10) \, M_\odot$ primordial black holes can account for all of the dark matter in the Milky Way is excluded at $5\sigma$ by a comparison with a VLA radio catalog at $1.4$ GHz, and at $\simeq 40\sigma$ by a comparison with a Chandra X-ray catalog ($0.5 - 8$ keV). We argue that this method can be used to identify such a population of primordial black holes with more sensitive future radio and X-ray surveys.

Numerical simulations of structure formation have recorded a remarkable progress in the recent years, in particular due to the inclusion of baryonic physics evolving with the dark matter component. We generate Monte Carlo realizations of the dark matter sub-halo population based on the results of the recent hydrodynamical simulation suite of Milky Way-sized galaxies. We then simulate the gamma-ray sky for both the setup of the 3FGL and 2FHL Fermi Large Area Telescope (LAT) catalogs, including the contribution from the annihilation of dark matter in the sub-halos. We find that the flux sensitivity threshold strongly depends on the particle dark matter mass, and more mildly also on its annihilation channel and the observation latitude. The results differ for the 3FGL and 2FHL catalogs, given their different energy thresholds. We also predict that the number of dark matter sub-halos among the unassociated sources is very small. A null number of detectable sub-halos in the Fermi-LAT 3FGL catalog would imply upper limits on the dark matter annihilation cross section into $b\bar{b}$ of $2 \cdot 10^{-26}$ ($5 \cdot 10^{-25}$) cm$^3$/s with $M_{\rm DM}$= 50 (1000) GeV. We find less than one extended sub-halo in the Fermi-LAT 3FGL catalog. As a matter of fact, the differences in the spatial and mass distribution of sub-halos between hydrodynamic and dark matter-only runs do not have significant impact on the gamma-ray dark matter phenomenology.

Analogously to globular clusters, the dense stellar environment of the Galactic center has been proposed to host a large population of as-yet undetected millisecond pulsars (MSPs). Recently, this hypothesis found support in the analysis of gamma rays from the inner Galaxy seen by the Large Area Telescope (LAT) aboard the Fermi satellite, which revealed a possible excess of diffuse GeV photons in the inner 15 deg about the Galactic center (Fermi GeV excess). The excess can be interpreted as the collective emission of thousands of MSPs in the Galactic bulge, with a spherical distribution that strongly peaks towards the Galactic center. In order to fully establish the MSP interpretation, it is essential to find corroborating evidence in multi-wavelength searches, most notably through the detection of radio pulsation from individual bulge MSPs. Based on globular cluster observations and the gamma-ray emission from the inner Galaxy, we investigate the prospects for detecting MSPs in the Galactic bulge. While previous pulsar surveys failed to identify this population, we demonstrate that, in the upcoming years, new large-area surveys with focus on regions a few degrees north or south of the Galactic center should lead to the detection of dozens of bulge MSPs. Additionally, we show that, in the near future, deep targeted searches of unassociated Fermi sources should be able to detect the first few MSPs in the bulge. The prospects for these deep searches are enhanced by a tentative gamma-ray/radio correlation that we infer from high-latitude gamma-ray MSPs. Such detections would constitute the first clear discoveries of field MSPs in the Galactic bulge, with far-reaching implications for gamma-ray observations, the formation history of the central Milky Way and strategy optimization for future radio observations.

We study the implications of galaxy formation on dark matter direct detection using high resolution hydrodynamic simulations of Milky Way-like galaxies simulated within the EAGLE and APOSTLE projects. We identify Milky Way analogues that satisfy observational constraints on the Milky Way rotation curve and total stellar mass. We then extract the dark matter density and velocity distribution in the Solar neighbourhood for this set of Milky Way analogues, and use them to analyse the results of current direct detection experiments. For most Milky Way analogues, the event rates in direct detection experiments obtained from the best fit Maxwellian distribution (with peak speed of 223 - 289 km/s) are similar to those obtained directly from the simulations. As a consequence, the allowed regions and exclusion limits set by direct detection experiments in the dark matter mass and spin-independent cross section plane shift by a few GeV compared to the Standard Halo Model, at low dark matter masses. For each dark matter mass, the halo-to-halo variation of the local dark matter density results in an overall shift of the allowed regions and exclusion limits for the cross section. However, the compatibility of the possible hints for a dark matter signal from DAMA and CDMS-Si and null results from LUX and SuperCDMS is not improved.

We present global fits of an effective field theory description of real, and complex scalar dark matter candidates. We simultaneously take into account all possible dimension 6 operators consisting of dark matter bilinears and gauge invariant combinations of quark and gluon fields. We derive constraints on the free model parameters for both the real (five parameters) and complex (seven) scalar dark matter models obtained by combining Planck data on the cosmic microwave background, direct detection limits from LUX, and indirect detection limits from the Fermi Large Area Telescope. We find that for real scalars indirect dark matter searches disfavour a dark matter particle mass below 100 GeV. For the complex scalar dark matter particle current data have a limited impact due to the presence of operators that lead to p-wave annihilation, and also do not contribute to the spin-independent scattering cross- section. Although current data are not informative enough to strongly constrain the theory parameter space, we demonstrate the power of our formalism to reconstruct the theoretical parameters compatible with an actual dark matter detection, by assuming that the excess of gamma rays observed by the Fermi Large Area Telescope towards the Galactic centre is entirely due to dark matter annihilations. Please note that the excess can very well be due to astrophysical sources such as millisecond pulsars. We find that scalar dark matter interacting via effective field theory operators can in principle explain the Galactic centre excess, but that such interpretation is in strong tension with the non-detection of gamma rays from dwarf galaxies in the real scalar case. In the complex scalar case there is enough freedom to relieve the tension.

We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, and identify the those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Way analogues so identified consistently exhibit very similar dark matter profiles inside the solar circle, therefore enabling more accurate predictions for indirect dark matter searches. We find in particular that high resolution simulated haloes satisfying observational constraints exhibit, within the inner few kiloparsecs, dark matter profiles shallower than those required to explain the so-called Fermi GeV excess via dark matter annihilation.

We obtain predictions for the properties of cold dark matter annihilation radiation using high resolution hydrodynamic zoom-in cosmological simulations of Milky Way-like galaxies (APOSTLE project) carried out as part of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) programme. Galactic halos in the simulation have significantly different properties from those assumed in the "standard halo model" often used in dark matter detection studies. The formation of the galaxy causes a contraction of the dark matter halo, whose density profile develops a steeper slope than the Navarro-Frenk-White (NFW) profile between $r\approx1.5$ kpc and $r\approx10$ kpc. At smaller radii, $r\lesssim1.5$ kpc, the halos develop a flatter than NFW slope. This unexpected feature may be specific to our particular choice of subgrid physics model but nevertheless the dark matter density profiles agree within 30% as the mass resolution is increased by a factor 150. The inner regions of the halos are almost perfectly spherical (axis ratios $b/a > 0.97$ within $r=1$ kpc) and there is no offset larger than 45 pc between the centre of the stellar distribution and the centre of the dark halo. The morphology of the predicted dark matter annihilation radiation signal is in broad agreement with $\gamma$-ray observations at large Galactic latitudes ($b\gtrsim3^\circ$). At smaller angles, the inferred signal in one of our four galaxies is similar to that which is observed but it is significantly weaker in the other three.

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

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

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

Gamma rays from young pulsars and milli-second pulsars are expected to contribute to the diffuse gamma-ray emission measured by the {\it Fermi} Large Area Telescope (LAT) at high latitudes. We derive the contribution of the pulsars undetected counterpart by using information from radio to gamma rays and we show that they explain only a small fraction of the isotropic diffuse gamma-ray background.

Millisecond Pulsars are second most abundant source population discovered by the Fermi-LAT. They might contribute non-negligibly to the diffuse emission measured at high latitudes by Fermi-LAT, the IDGRB. Gamma-ray sources also contribute to the anisotropy of the IDGRB measured on small scales by Fermi-LAT. We aim to assess the contribution of the unresolved counterpart of the detected MSPs population to the IDGRB and the maximal fraction of the measured anisotropy produced by this source class. We model the MSPs spatial distribution in the Galaxy and the gamma-ray emission parameters by considering radio and gamma-ray observational constraints. By simulating a large number of MSPs populations, we compute the average diffuse emission and the anisotropy 1-sigma upper limit. The emission from unresolved MSPs at 2 GeV, where the peak of the spectrum is located, is at most 0.9% of the measured IDGRB above 10 degrees in latitude. The 1-sigma upper limit on the angular power for unresolved MSP sources turns out to be about a factor of 60 smaller than Fermi-LAT measurements above 30 degrees. Our results indicate that this galactic source class represents a negligible contributor to the high-latitude gamma-ray sky and confirm that most of the intensity and geometrical properties of the measured diffuse emission are imputable to other extragalactic source classes. Nevertheless, given the MSP distribution, we expect them to contribute significantly to the gamma-ray diffuse emission at low latitudes. Since, along the galactic disk, the population of young Pulsars overcomes in number the one of MSPs, we compute the gamma-ray emission from the whole population of unresolved Pulsars in two low-latitude regions: the inner Galaxy and the galactic center.

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

The annihilation of dark matter particles in the halo of galaxies may end up into gamma rays, which travel almost unperturbed till to their detection at Earth. This annihilation signal can exhibit an anisotropic behavior quantified by the angular power spectrum, whose properties strongly depend on the dark matter distribution and its clumpiness. We use high resolution pure dark matter N-body simulations to quantify the contribution of different components (main halo and satellites) to the global signal as a function of the analytical profile adopted to describe the numerical results. We find that the smooth main halo dominates the angular power spectrum of the gamma-ray signal up to quite large multipoles, where the sub-haloes anisotropy signal starts to emerge, but the transition multipole strongly depends on the assumed radial profile. The extrapolation down to radii not resolved by current numerical simulations can affect both the normalization and the shape of the gamma-ray angular power spectrum. For the sub-haloes described by an asymptotically cored dark matter distribution, the angular power spectrum shows an overall smaller normalization and a flattening at high multipoles. Our results show the criticality of the dark matter density profile shape in gamma-ray anisotropy searches, and evaluate quantitatively the intrinsic errors occurring when extrapolating the dark matter radial profiles down to spatial scales not yet explored by numerical simulations.

Indirect searches for the cosmological dark matter have become ever more competitive during the past years. Here, we report the first full calculation of leading electroweak corrections to the annihilation rate of supersymmetric neutralino dark matter. We find that these corrections can be huge, partially due to contributions that have been overlooked so far. Our results imply a significantly enhanced discovery potential of this well motivated dark matter candidate with current and upcoming cosmic ray experiments, in particular for gamma rays and models with somewhat small annihilation rates at tree level.

The nature of the Isotropic $\gamma$-ray Background (IGRB) measured by the Large Area Telescope (LAT) on the Fermi $\gamma$-ray space Telescope ({\it Fermi}) remains partially unexplained. Non-negligible contributions may originate from extragalactic populations of unresolved sources such as blazars, star-forming galaxies or galactic milli-second pulsars. A recent prediction of the diffuse $\gamma$-ray emission from Active Galactic Nuclei (AGN) with a large viewing angle with respect to the line-of-sight (l.o.s.) has demonstrated that this faint but numerous population is also expected to contribute significantly to the total IGRB intensity. A more exotic contribution to the IGRB invokes the pair annihilation of dark matter (DM) weakly interacting massive particles (WIMPs) into $\gamma$-rays. In this work, we evaluate the room left for galactic DM at high latitudes ($>10^\circ$) by including photons from both prompt emission and inverse Compton scattering, emphasizing the impact of the newly discovered contribution from misaligned AGN (MAGN) for such an analysis. Summing up all significant galactic and extragalactic components of the IGRB, we find that an improved understanding of the associated astrophysical uncertainties is still mandatory to put stringent bounds on thermally produced DM. On the other hand, we also demonstrate that the IGRB has the potential to be one of the most competitive {\it future} ways to test the DM WIMP hypothesis, once the present uncertainties are even slightly reduced. In fact, if MAGN contribute even at 90% of the maximal level consistent with our current understanding, thermally produced WIMPs would be severely constrained as DM candidates for masses up to several TeV.

Active galactic nuclei (AGN) with jets seen at small viewing angles are the most luminous and abundant objects in the $\gamma$-ray sky. AGN with jets misaligned along the line-of-sight appear fainter in the sky, but are more numerous than the brighter blazars. We calculate the diffuse $\gamma$-ray emission due to the population of misaligned AGN (MAGN) unresolved by the Large Area Telescope (LAT) on the {\it Fermi} Gamma-ray Space Telescope ({\it Fermi}). A correlation between the $\gamma$-ray luminosity and the radio-core luminosity is established and demonstrated to be physical by statistical tests, as well as compatible with upper limits based on {\it Fermi}-LAT data for a large sample of radio-loud MAGN. We constrain the derived $\gamma$-ray luminosity function by means of the source count distribution of the radio galaxies (RGs) detected by the {\it Fermi}-LAT. We finally calculate the diffuse $\gamma$-ray flux due to the whole MAGN population. Our results demonstrate that the MAGN can contribute from 10% up to nearly the entire measured Isotropic Gamma-Ray Background (IGRB). We evaluate a theoretical uncertainty on the flux of almost an order of magnitude.

We calculate the diffuse $\gamma$-ray emission due to the population of misaligned AGN (MAGN) unresolved by the Large Area Telescope (LAT) on the {\it Fermi} Gamma-ray Space Telescope ({\it Fermi}). A correlation between the $\gamma$-ray luminosity and the radio-core luminosity is established and demonstrated to be physical by statistical tests, as well as compatible with upper limits based on {\it Fermi}-LAT data for a large sample of radio-loud MAGN. We constrain the derived $\gamma$-ray luminosity function by means of the source count distribution of the MAGN detected by the {\it Fermi}-LAT. We finally estimate the diffuse $\gamma$-ray flux due to the whole MAGN population which ranges from 10% up to nearly the entire measured Isotropic Gamma-Ray Background (IGRB). We evaluate also the room left to galactic DM at high latitudes ($>10^\circ$), by taking into account the results on the MAGN together with the other significant galactic and extragalactic $\gamma$-rays emitting sources.

The spectrum of an isotropic extragalactic $\gamma$-ray background (EGB) has been measured by the Fermi-LAT telescope at high latitudes. Two new models for the EGB are derived from the subtraction of unresolved point sources and extragalactic diffuse processes, which could explain from 30% to 70% of the Fermi-LAT EGB. Within the hypothesis that the two residual EGBs are entirely due to the annihilation of dark matter (DM) particles in the Galactic halo, we obtain $conservative$ upper limits on their annihilation cross section \sigmav. Severe bounds on a possible Sommerfeld enhancement of the annihilation cross section are set as well. Finally, would {\sigmav} be inversely proportional to the WIMP velocity, very severe limits are derived for the velocity-independent part of the annihilation cross section.

The spectrum of an isotropic extragalactic {\gamma}-ray background (EGB) has been measured by the Fermi-LAT telescope at high latitudes. Two new models for the EGB are derived from the subraction of unresolved point sources and extragalactic diffuse processes, which could explain from 30% to 70% of the Fermi-LAT EGB. Within the hypothesis that the two residual EGBs are entirely due to the annihilation of dark matter (DM) particles in the Galactic halo, we obtain stringent upper limits on their annihilation cross section. Severe bounds on a possible Sommerfeld enhancement of the annihilation cross section are set as well. Finally, we consider models for DM annihilation depending on the inverse of the velocity and associate the EGBs to photons arising from the annihilation of DM in primordial halos. For DM velocities in protohalos of v ~ 10^-8 c, the annihilation cross section is constrained down to 10^-33 cm^3 /s for DM masses below 100 GeV. Given our choices for the EGB and the minimal DM modelling, the derived upper bounds are claimed to be conservative.

Gamma rays from the annihilation of dark matter particles in the Galactic halo provide a particularly promising means of indirectly detecting dark matter. Here, we demonstrate that pronounced spectral features at energies near the dark matter particles' mass, which are a generic prediction for most models, can significantly improve the sensitivity of gamma-ray telescopes to dark matter signals. We derive projected limits on such features (including the traditionally looked-for line signals) and show that they can be much more efficient in constraining the nature of dark matter than the model-independent broad spectral features expected at lower energies.

Indirect searches for dark matter annihilation or decay products in the cosmic-ray spectrum are plagued by the question of how to disentangle a dark matter signal from the omnipresent astrophysical background. One of the practically background-free smoking-gun signatures for dark matter would be the observation of a sharp cutoff or a pronounced bump in the gamma-ray energy spectrum. Such features are generically produced in many dark matter models by internal Bremsstrahlung, and they can be treated in a similar manner as the traditionally looked-for gamma-ray lines. Here, we discuss prospects for seeing such features with present and future Atmospheric Cherenkov Telescopes.