A prediction of the standard Lambda cold dark matter cosmology is that dark matter (DM) haloes are teeming with numerous self-bound substructure or subhaloes. The precise properties of these subhaloes represent important probes of the underlying cosmological model. We use data from Via Lactea II and Exploring the Local Volume in Simulations N-body simulations to learn about the structure of subhaloes with masses 10(6)-10(11) h(-1) M circle dot. Thanks to a superb subhalo statistics, we study subhalo properties as a function of distance to host halo centre and subhalo mass, and provide a set of fits that accurately describe the subhalo structure. We also investigate the role of subhaloes on the search for DM annihilation. Previous work has shown that subhaloes are expected to boost the DM signal of their host haloes significantly. Yet, these works traditionally assumed that subhaloes exhibit similar structural properties than those of field haloes, while it is known that subhaloes are more concentrated. Building upon our N-body data analysis, we refine the substructure boost model of Sanchez-Conde & Prada (2014), and find boosts that are a factor 2-3 higher. We further refine the model to include unavoidable tidal stripping effects on the subhalo population. For field haloes, this introduces a moderate (similar to 20-30 per cent) suppression. Yet, for subhaloes like those hosting dwarf galaxy satellites, tidal stripping plays a critical role, the boost being at the level of a few tens of percent at most. We provide a parametrization of the boost for field haloes that can be safely applied over a wide halo mass range.

The Cherenkov Telescope Array (CTA) is about to enter construction phase and one of its main key science projects is to perform an unbiased survey in search of extragalactic sources. We make use of both the latest blazar gamma-ray luminosity function and spectral energy distribution to derive the expected number of detectable sources for both the planned Northern and Southern arrays of the CTA observatory. We find that a shallow, wide survey of about 0.5 hour per field of view would lead to the highest number of blazar detections. Furthermore, we investigate the effect of axion-like particles and secondary gamma rays from propagating cosmic rays on the source count distribution, since these processes predict different spectral shape from standard extragalactic background light attenuation. We can generally expect more distant objects in the secondary gamma-ray scenario, while axion-like particles do not significantly alter the expected distribution. Yet, we find that, these results strongly depend on the assumed magnetic field strength during the propagation. We also provide source count predictions for the High Altitude Water Cherenkov observatory (HAWC), the Large High Altitude Air Shower Observatory (LHAASO) and a novel proposal of a hybrid detector.

Line-of-sight integrals of the squared density, commonly called the J-factor, are essential for inferring dark matter (DM) annihilation signals. The J-factors of DM-dominated dwarf spheroidal satellite galaxies (dSphs) have typically been derived using Bayesian techniques, which for small data samples implies that a choice of priors constitutes a non-negligible systematic uncertainty. Here we report the development of a new fully frequentist approach to construct the profile likelihood of the J-factor. Using stellar kinematic data from several classical and ultra-faint dSphs, we derive the maximum likelihood value for the J-factor and its confidence intervals. We validate this method, in particular its bias and coverage, using simulated data from the Gaia Challenge. We find that the method possesses good statistical properties. The J-factors and their uncertainties are generally in good agreement with the Bayesian-derived values, with the largest deviations restricted to the systems with the smallest kinematic data sets. We discuss improvements, extensions, and future applications of this technique.

In a Galactic core-collapse supernova (SN), axionlike particles (ALPs) could be emitted via the Primakoff process and eventually convert into gamma rays in the magnetic field of the Milky Way. From a data-driven sensitivity estimate, we find that, for a SN exploding in our Galaxy, the Fermi Large Area Telescope (LAT) would be able to explore the photon-ALP coupling down to g(a gamma) similar or equal to 2 x 10(-13) GeV-1 for an ALP mass m(a) less than or similar to 10(-9) eV. These values are out of reach of next generation laboratory experiments. In this event, the Fermi LAT would probe large regions of the ALP parameter space invoked to explain the anomalous transparency of the Universe to gamma rays, stellar cooling anomalies, and cold dark matter. If no gamma-ray emission were to be detected, Fermi-LAT observations would improve current bounds derived from SN 1987A by more than 1 order of magnitude.

The Fermi Large Area Telescope (LAT) has opened the way for comparative studies of cosmic rays (CRs) and high-energy objects in the Milky Way (MW) and in other, external, star-forming galaxies. Using 2 yr of observations with the Fermi LAT, Local Group galaxy M31 was detected as a marginally extended gamma-ray source, while only an upper limit has been derived for the other nearby galaxy M33. We revisited the gamma-ray emission in the direction of M31 and M33 using more than 7 yr of LAT Pass 8 data in the energy range 0.1-100 GeV, presenting detailed morphological and spectral analyses. M33 remains undetected, and we computed an upper limit of 2.0 x 10(-12) erg cm(-2) s(-1) on the 0.1-100 GeV energy flux (95% confidence level). This revised upper limit remains consistent with the observed correlation between gamma-ray luminosity and star formation rate tracers and implies an average CR density in M33 that is at most half of that of the MW. M31 is detected with a significance of nearly 10 sigma. Its spectrum is consistent with a power law with photon index Gamma = 2.4 +/- 0.1(stat) (vertical bar) (syst) and a 0.1-100 GeV energy flux of (5.6 +/- 0.6(stat vertical bar syst)) x 10(-12) erg cm(-1) s(-1). M31 is detected to be extended with a 4 sigma significance. The spatial distribution of the emission is consistent with a uniform-brightness disk with a radius of 0 degrees.4 and no offset from the center of the galaxy, but nonuniform intensity distributions cannot be excluded. The flux from M31 appears confined to the inner regions of the galaxy and does not fill the disk of the galaxy or extend far from it. The gamma-ray signal is not correlated with regions rich in gas or star formation activity, which suggests that the emission is not interstellar in origin, unless the energetic particles radiating in gamma rays do not originate in recent star formation. Alternative and nonexclusive interpretations are that the emission results from a population of millisecond pulsars dispersed in the bulge and disk of M31 by disrupted globular clusters or from the decay or annihilation of dark matter particles, similar to what has been proposed to account for the so-called Galactic center excess found in Fermi-LAT observations of the MW.

We search for excess gamma-ray emission coincident with the positions of confirmed and candidate Milky Way satellite galaxies using six years of data from the Fermi Large Area Telescope (LAT). Our sample of 45 stellar systems includes 28 kinematically confirmed dark-matter-dominated dwarf spheroidal galaxies (dSphs) and 17 recently discovered systems that have photometric characteristics consistent with the population of known dSphs. For each of these targets, the relative predicted gamma-ray flux due to dark matter annihilation is taken from kinematic analysis if available, and estimated from a distance-based scaling relation otherwise, assuming that the stellar systems are DM-dominated dSphs. LAT data coincident with four of the newly discovered targets show a slight preference (each similar to 2 sigma local) for gamma-ray emission in excess of the background. However, the ensemble of derived gamma-ray flux upper limits for individual targets is consistent with the expectation from analyzing random blank-sky regions, and a combined analysis of the population of stellar systems yields no globally significant excess (global significance < 1 sigma). Our analysis has increased sensitivity compared to the analysis of 15 confirmed dSphs by Ackermann et al. The observed constraints on the DM annihilation cross section are statistically consistent with the background expectation, improving by a factor of similar to 2 for large DM masses (m(DM, b<(b)over bar>) greater than or similar to 1 TeV and m(DM, tau+tau-) greater than or similar to 70 GeV) and weakening by a factor of similar to 1.5 at lower masses relative to previously observed limits.

Imaging atmospheric Cherenkov telescopes (IACTs) that are sensitive to potential.-ray signals from dark matter (DM) annihilation above similar to 50 GeV will soon be superseded by the Cherenkov Telescope Array (CTA). CTA will have a point source sensitivity an order of magnitude better than currently operating IACTs and will cover a broad energy range between 20 GeV and 300 TeV. Using effective field theory and simplified models to calculate gamma-ray spectra resulting from DM annihilation, we compare the prospects to constrain such models with CTA observations of the Galactic center with current and near-future measurements at the Large Hadron Collider (LHC) and direct detection experiments. For DM annihilations via vector or pseudoscalar couplings, CTA observations will be able to probe DM models out of reach of the LHC, and, if DM is coupled to standard fermions by a pseudoscalar particle, beyond the limits of current direct detection experiments.

The region around the Galactic Center (GC) is now well established to be brighter at energies of a few GeV than what is expected from conventional models of diffuse gamma-ray emission and catalogs of known gamma-ray sources. We study the GeV excess using 6.5 yr of data from the Fermi Large Area Telescope. We characterize the uncertainty of the GC excess spectrum and morphology due to uncertainties in cosmic-ray source distributions and propagation, uncertainties in the distribution of interstellar gas in the Milky Way, and uncertainties due to a potential contribution from the Fermi bubbles. We also evaluate uncertainties in the excess properties due to resolved point sources of gamma rays. The GC is of particular interest, as it would be expected to have the brightest signal from annihilation of weakly interacting massive dark matter (DM) particles. However, control regions along the Galactic plane, where a DM signal is not expected, show excesses of similar amplitude relative to the local background. Based on the magnitude of the systematic uncertainties, we conservatively report upper limits for the annihilation cross-section as a function of particle mass and annihilation channel.

We present a catalog of sources detected above 50 GeV by the Fermi-Large Area Telescope (LAT) in 80 months of data. The newly delivered Pass. 8 event-level analysis allows the detection and characterization of sources in the 50 GeV-2 TeV energy range. In this energy band, Fermi-LAT. has detected 360 sources, which constitute the second catalog of hard Fermi-LAT. sources (2FHL). The improved angular resolution enables the precise localization of point sources (similar to 1.' 7 radius at 68% C.L.) and the detection and characterization of spatially extended sources. We find that 86% of the sources can be associated with counterparts at other wavelengths, of which the majority (75%) are active galactic nuclei and the rest (11%) are Galactic sources. Only 25% of the 2FHL sources have been previously detected by Cherenkov telescopes, implying that the 2FHL provides a reservoir of candidates to be followed up at very high energies. This work closes the energy gap between the observations performed at GeV energies by Fermi-LAT. on orbit and the observations performed at higher energies by Cherenkov telescopes from the ground.

Nearly one-third of the sources listed in the Third Fermi Large Area Telescope (LAT) catalog (3FGL) remain unassociated. It is possible that predicted and even unanticipated gamma-ray source classes are present in these data waiting to be discovered. Taking advantage of the excellent spectral capabilities achieved by the Fermi LAT, we use machine-learning classifiers (Random Forest and XGBoost) to pinpoint potentially novel source classes in the unassociated 3FGL sample outside the Galactic plane. Here we report a total of 34 high-confidence Galactic candidates at vertical bar b vertical bar >= 5 degrees. The currently favored standard astrophysical interpretations for these objects are pulsars or low-luminosity globular clusters hosting millisecond pulsars (MSPs). Yet these objects could also be interpreted as dark matter annihilation taking place in ultra-faint dwarf galaxies or dark matter subhalos. Unfortunately, Fermi LAT spectra are not sufficient to break degeneracies between the different scenarios. Careful visual inspection of archival optical images reveals no obvious evidence for low-luminosity globular clusters or ultra-faint dwarf galaxies inside the 95% error ellipses. If these are pulsars, this would bring the total number of MSPs at vertical bar b vertical bar >= 5 degrees to 106, down to an energy flux approximate to 4.0 x 10(-12) erg cm(-2) s(-1) between 100 MeV and 100 GeV. We find this number to be in excellent agreement with predictions from a new population synthesis of MSPs that predicts 100-126 high-latitude 3FGL MSPs depending on the choice of high-energy emission model. If, however, these are dark matter substructures, we can place upper limits on the number of Galactic subhalos surviving today and on dark matter annihilation cross sections. These limits are beginning to approach the canonical thermal relic cross section for dark matter particle masses below similar to 100 GeV in the bottom quark (b (b) over bar) annihilation channel.