Galactic Center Excess Research Articles

Galactic Center excess at the highest energies: Morphology and photon-count statistics

Galactic Center excess at the highest energies: Morphology and photon-count statistics

Robust inference of the Galactic Centre gamma-ray excess spatial properties

ABSTRACT The gamma-ray Fermi-LAT Galactic Centre excess (GCE) has puzzled scientists for over 15 yr. Despite ongoing debates about its properties, and especially its spatial distribution, its nature remains elusive. We scrutinize how the estimated spatial morphology of this excess depends on models for the Galactic diffuse emission, focusing particularly on the extent to which the Galactic plane and point sources are masked. Our main aim is to compare a spherically symmetric morphology – potentially arising from the annihilation of dark matter (DM) particles – with a boxy morphology – expected if faint unresolved sources in the Galactic bulge dominate the excess emission. Recent claims favouring a DM-motivated template for the GCE are shown to rely on a specific Galactic bulge template, which performs worse than other templates for the Galactic bulge. We find that a non-parametric model of the Galactic bulge derived from the VISTA Variables in the Via Lactea survey results in a significantly better fit for the GCE than DM-motivated templates. This result is independent of whether a galprop-based model or a more non-parametric ring-based model is used to describe the diffuse Galactic emission. This conclusion remains true even when additional freedom is added in the background models, allowing for non-parametric modulation of the model components and substantially improving the fit quality. When adopted, optimized background models provide robust results in terms of preference for a boxy bulge morphology for the GCE, regardless of the mask applied to the Galactic plane.

Improved Galactic diffuse emission model strengthens the case for a Millisecond Pulsar explanation of the Fermi GeV excess

After more than a decade since its discovery, the Galactic center gamma-ray excess - discovered with the Fermi Large Area Telescope - remains puzzling. While the spectrum of the signal can be explained by either dark matter or an unresolved population of millisecond pulsars, the spatial morphology of this excess seems to hold the key to separate the two theories. In this contribution, we present the results of a recent study in which we use bleeding edge models for interstellar gas, inverse Compton emission, and stellar mass models to reanalyze the Galactic center excess. We find that the spatial morphology of the excess is highly correlated with stellar matter in the Galactic bulge, providing strong support for the millisecond pulsar hypothesis.

Mind the gap: the discrepancy between simulation and reality drives interpretations of the Galactic Center Excess

The Galactic Center Excess (GCE) in GeV gamma rays has been debated for over a decade, with the possibility that it might be due to dark matter annihilation or undetected point sources such as millisecond pulsars (MSPs). This study investigates how the gamma-ray emission model (γEM) used in Galactic center analyses affects the interpretation of the GCE's nature.To address this issue, we construct an ultra-fast and powerful inference pipeline based on convolutional Deep Ensemble Networks. We explore the two main competing hypotheses for the GCE using a set of γEMs with increasing parametric freedom. We calculate the fractional contribution (f src) of a dim population of MSPs to the total luminosity of the GCE and analyze its dependence on the complexity of the γEM. For the simplest γEM, we obtain f src = 0.10 ± 0.07, while the most complex model yields f src = 0.79 ± 0.24.In conclusion, we find that the statement about the nature of the GCE (dark matter or not) strongly depends on the assumed γEM.The quoted results for f src do not account for the additional uncertainty arising from the fact that the observed gamma-ray sky is out-of-distribution concerning the investigated γEM iterations. We quantify the reality gap between our γEMs using deep-learning-based One-Class Deep Support Vector Data Description networks, revealing that all employed γEMs have gaps to reality.Our study casts doubt on the validity of previous conclusions regarding the GCE and dark matter, and underscores the urgent need to account for the reality gap and consider previously overlooked “out of domain” uncertainties in future interpretations.

Search for dark matter lines at the Galactic Center with 14 years of Fermi data

Dark matter (DM) in the Milky Way halo may annihilate or decay to photons, producing monochromatic gamma rays. We search for DM-induced spectral lines using 14 years of data from the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope ($\textit{Fermi}$-LAT) between $10\,\mathrm{GeV}$ and $2\,\mathrm{TeV}$ in the inner Milky Way leveraging both the spatial and spectral morphology of an expected signal. We present new constraints as strong as $\langle \sigma v \rangle \lesssim 6\times 10^{-30}\, \mathrm{cm}^3/\mathrm{s}$ for the two-to-two annihilations and $\tau \gtrsim 10^{30}\,\mathrm{s}$ for one-to-two decays, representing leading sensitivity between $10\,\mathrm{GeV}$ and $\sim$$500\,\mathrm{GeV}$. We consider the implications of our line-constraints on the Galactic Center Excess (GCE), which is a previously-observed excess of continuum $\sim$GeV gamma-rays that may be explained by DM annihilation. The Higgs portal and neutralino-like DM scenarios, which have been extensively discussed as possible origins of the GCE, are constrained by our work because of the lack of observed one-loop decays to two photons. More generally, we interpret our null results in a variety of annihilating and decaying DM models, such as neutralinos, gravitinos, and glueballs, showing that in many cases the line search is more powerful than the continuum, despite the continuum annihilation being at tree level.

Characterizing the expected behavior of non-Poissonian template fitting

We have performed a systematic study of the statistical behavior of non-Poissonian template fitting (NPTF), a method designed to analyze and characterize unresolved point sources in general counts datasets. In this paper, we focus on the properties and characteristics of the Fermi-LAT gamma-ray data set. In particular, we have simulated and analyzed gamma-ray sky maps under varying conditions of exposure, angular resolution, pixel size, energy window, event selection, and source brightness. We describe how these conditions affect the sensitivity of NPTF to the presence of point sources, for inner-galaxy studies of point sources within the Galactic Center excess, and for the simplified case of isotropic emission. We do not find opportunities for major gains in sensitivity from varying these choices, within the range available with current Fermi-LAT data. We provide an analytic estimate of the NPTF sensitivity to point sources for the case of isotropic emission and perfect angular resolution, and find good agreement with our numerical results for that case.

Deep learning models of the discrete component of the Galactic interstellar γ -ray emission

A significant point-like component from the small scale (or discrete) structure in the H2 interstellar gas might be present in the Fermi-LAT data, but modeling this emission relies on observations of rare gas tracers only available in limited regions of the sky. Identifying this contribution is important to discriminate gamma-ray point sources from interstellar gas, and to better characterize extended gamma-ray sources. We design and train convolutional neural networks to predict this emission where observations of these rare tracers do not exist and discuss the impact of this component on the analysis of the Fermi-LAT data. In particular, we evaluate prospects to exploit this methodology in the characterization of the Fermi-LAT Galactic center excess through accurate modeling of point-like structures in the data to help distinguish between a point-like or smooth nature for the excess. We show that deep learning may be effectively employed to model the gamma-ray emission traced by these rare H2 proxies within statistical significance in data-rich regions, supporting prospects to employ these methods in yet unobserved regions.

On the morphology of the gamma-ray galactic centre excess

ABSTRACT The characteristics of the galactic centre excess (GCE) emission observed in gamma-ray energies – especially the morphology of the GCE – remain a hotly debated subject. The manner in which the dominant diffuse gamma-ray background is modelled has been claimed to have a determining effect on the preferred morphology. In this work, we compare two distinct approaches to the galactic diffuse gamma-ray emission background: the first approach models this emission through templates calculated from a sequence of well-defined astrophysical assumptions, while the second approach divides surrogates for the background gamma-ray emission into cylindrical galactocentric rings with free independent normalizations. At the latitudes that we focus on, we find that the former approach works better, and that the overall best fit is obtained for an astrophysically motivated fit when the GCE follows the morphology expected of dark matter annihilation. Quantitatively, the improvement compared with the best ring-based fits is roughly 6500 in the χ2 and roughly 4000 in the log of the Bayesian evidence.

The velocity-dependent J-factor of the Milky Way halo: does what happens in the galactic bulge stay in the galactic bulge?

We consider the angular distribution of the photon signal which could arisefrom velocity-dependent dark matter annihilation within the Galactic bulge. We find that, for thecase of Sommerfeld-enhanced annihilation, dark matter annihilation within the bulge isdominated by slow speed particles which never leave the bulge, allowing one to find a simpleanalytic relationship between the dark matter profile within the Galactic bulge and the angulardistribution. On the other hand, for the case p- or d-wave annihilation, we find thatthe small fraction of high-speed particles which can leave the bulge provide a significant,often dominant, contribution to dark matter annihilation within the bulge. For these scenarios,fully understanding dark matter annihilation deep within the Galactic bulge, and the angular distribution ofthe resulting photon signal, requires an understanding ofthe dark matter profile well outside the bulge. We consider the Galactic Center excess inlight of these results, and find that an explanation of this excess in terms of p-waveannihilation would require the dark matter profile within the bulge to have a much steeperprofile than usually considered, but with uncertainties related to the behavior of the profileoutside the bulge.

Velocity-dependent annihilation radiation from dark matter subhalos in cosmological simulations

We use the suite of Milky Way-like galaxies in the Auriga simulations to determine the contribution to annihilation radiation from dark matter subhalos in three velocity-dependent dark matter annihilation models: Sommerfeld, p-wave, and d-wave models. We compare these to the corresponding distribution in the velocity-independent s-wave annihilation model. For both the hydrodynamical and dark-matter-only simulations, only in the case of the Sommerfeld-enhanced annihilation does the total annihilation flux from subhalos exceed the total annihilation flux from the smooth halo component within the virial radius of the halo. Progressing from Sommerfeld to the s, p, and d-wave models, the contribution from the smooth component of the halo becomes more dominant, implying that for the p-wave and d-wave models the smooth component is by far the dominant contribution to the radiation. Comparing to the Galactic center excess observed by Fermi-LAT, for all simulated halos the emission is dominated by the smooth halo contribution. However, it is possible that for Sommerfeld models, extrapolation down to mass scales below the current resolution limit of the simulation would imply a non-negligible contribution to the gamma-ray emission from the Galactic Center region.

Updated constraints on WIMP dark matter annihilation by radio observations of M31

The present work derived the robust constraints on annihilating WIMP parameters utilizing new radio observations of M31, as well as new studies of its dark matter distribution and other properties. The characteristics of emission due to DM annihilation were computed in the frame of 2D galactic model employing GALPROP code adapted specifically for M31. This enabled us to refine various inaccuracies of previous studies on the subject. DM constraints were obtained for two representative annihilation channels: $\chi\chi \rightarrow b\overline{b}$ and $\chi\chi \rightarrow \tau^+\tau^-$. A wide variety of radio data was utilized in the frequency range $\approx$(0.1--10) GHz. As the result the thermal WIMP lighter than fiducially $\approx$ 70 GeV for $b\overline{b}$ channel and $\approx$ 40 GeV for $\tau^+\tau^-$ was excluded. The corresponding mass threshold uncertainty ranges were estimated to be 20--210 GeV and 18--89 GeV. The obtained exclusions are competitive to those from Fermi-LAT observations of dwarfs and AMS-02 measurements of antiprotons. Our constraints do not exclude the explanation of the gamma-ray outer halo of M31 and the Galactic center excess by annihilating DM. The thermal WIMP with $m_x \approx 70$ GeV, which explains the outer halo, would make a significant contribution to the non-thermal radio flux in M31 nucleus, fitting well both the spectrum and morphology. And, finally, we questioned the possibility claimed in other studies to robustly constrain heavy thermal WIMP with $m_x > 100$ GeV by radio data on M31.

Muonphilic dark matter explanation of gamma-ray galactic center excess: a comprehensive analysis

The Galactic center gamma-ray excess (GCE) is a long-standing unsolved problem. One of candidate solutions, the dark matter (DM) annihilation, has been recently tested with other astrophysical observations, such as AMS-02 electron-positron spectra, Fermi Dwarf spheroidal galaxies gamma-ray data, and so on. By assuming that the DM particles annihilate purely into a normal charged fermion pair, Di Mauro and Winkle (2021) claimed that only a muon-pair is compatible with the null detection of all the corresponding astrophysical measurements and can explain GCE simultaneously. On the other hand, a muonphilic DM model may also lead to a signal in the recent Fermilab muon g − 2 measurement or be constrained by the latest PandaX-4T limit. In this work, we comprehensively study interactions between DM and muon, including various combinations of DM and mediator spins. In agreement with GCE (not only 2μ but also 4μ final states), we test these interactions against all the thermal DM constraints. Our results show that only the parameter space near the resonance region of mediator can explain GCE and relic density simultaneously, and larger parameter spaces are still allowed if other poorly-known systematic uncertainties are included. Regardless of the DM spin, only the interactions with the spin-0 mediator can explain the recent muon g − 2 excess on top of GCE, relic density, and other DM and mediator constraints.

Exploring the MeV sky with a combined coded mask and Compton telescope: the Galactic Explorer with a Coded aperture mask Compton telescope (GECCO)

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds upon the heritage of past and current missions improving the sensitivity and, very importantly, the angular resolution. This consists in combining a Compton telescope and a coded-mask telescope. We delineate the motivation for such a concept and we define the scientific goals for such a mission.The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies.The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented.With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which, in turn, enables identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena.

Luminosity functions consistent with a pulsar-dominated Galactic Center excess

A new population of millisecond pulsars is a long-standing proposed explanation for the excess of GeV-scale gamma rays emanating from the region surrounding the center of the Milky Way (the “Galactic Center excess”). We examine several simple parameterizations of possible luminosity functions for this population, as well as several benchmark luminosity functions proposed in the literature, and compare the predicted populations of resolved point sources to the Fermi 4FGL-DR2 point source catalog and a sub-population recently identified using wavelet-based methods. We provide general results that can be used to translate upper limits on the number of resolved point sources associated with the excess, and the fraction of the flux in the excess that can be attributed to resolved sources, into limits on the luminosity function parameter space. We discuss a number of important systematic uncertainties, including in the detection threshold model and the total flux attributed to the excess. We delineate regions of parameter space (containing existing benchmark models) where there is no apparent tension with current data, and the number of total pulsars needed to explain the excess is in the range of \U0001d4aa(104-5). In the future, lowered point source detection thresholds could be achieved either by new analysis methods or new data. An order-of-magnitude reduction in the sensitivity threshold (which may already be achieved by novel analyses probing sub-threshold source populations) could hope to resolve more than 30% of the flux of the excess even in pessimistic scenarios.

Return of the templates: Revisiting the Galactic Center excess with multimessenger observations

The Galactic center excess (GCE) remains one of the most intriguing discoveries from the Fermi Large Area Telescope (LAT) observations. We revisit the characteristics of the GCE by first producing a new set of high-resolution galactic diffuse gamma-ray emission templates. This diffuse emission, which accounts for the bulk of the observed gamma rays, is ultimately due to cosmic-ray interactions with the interstellar medium. Using recent high-precision cosmic-ray observations, in addition to the continuing Fermi-LAT observations and observations from lower energy photons, we constrain the properties of the galactic diffuse emission. We describe a large set of diffuse gamma-ray emission templates which account for a very wide range of initial assumptions on the physical conditions in the inner galaxy. The broad properties of the GCE that we find in this work are qualitatively unchanged despite the introduction of this new set of templates, though its quantitative features appear mildly different than those obtained in previous analyses. In particular, we find a high-energy tail at higher significance than previously reported. This tail is very prominent in the northern hemisphere, and less so in the southern hemisphere. This strongly affects one prominent interpretation of the excess: known millisecond pulsars are incapable of producing this high-energy emission, even in the relatively softer southern hemisphere, and are therefore disfavored as the sole explanation of the GCE. The annihilation of dark matter particles of mass ${40}_{\ensuremath{-}7}^{+10}\text{ }\text{ }\mathrm{GeV}$ (95% CL) to $b$ quarks with a cross-section of $⟨{\ensuremath{\sigma}}_{A}v⟩=1.{4}_{\ensuremath{-}0.3}^{+0.6}\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}26}\text{ }\text{ }{\mathrm{cm}}^{3}\text{ }{\mathrm{s}}^{\ensuremath{-}1}$ provides a good fit to the excess especially in the relatively cleaner southern sky. Dark matter of the same mass range annihilating to $b$ quarks or heavier dark matter particles annihilating to heavier Standard Model bosons can combine with millisecond pulsars to provide a good fit to the southern hemisphere emission as well, as can a broken power-law spectrum which would be related to recent cosmic-ray burst activity. As part of this paper, we make publicly available all of our templates and the data covariance matrix we have generated to account for systematic uncertainties.

Dark matter annihilation and the Galactic Centre Excess

ABSTRACT We compare the surface brightness profile and morphology of the Galactic Centre Excess (GCE) identified in wide-angle γ-ray maps from the Fermi-Large Area Telescope (LAT) to dark matter annihilation predictions derived from high-resolution Λ cold dark matter magnetohydrodynamic simulations of galaxy formation. These simulations produce isolated, disc-dominated galaxies with structure, stellar populations, gas content, and stellar and halo masses comparable to those of the Milky Way. For a specific choice of annihilation cross-section, they agree well with the Fermi-LAT data over the full observed angular range, 1°–15°, whereas their dark-matter-only counterparts, lacking any compression of the inner halo by the gravitational effects of the baryons, fail to predict emission as centrally concentrated as observed. These results provide additional support to the hypothesis that the GCE is produced by annihilating dark matter. If, however, it is produced by a different mechanism, they imply a strong upper limit on annihilation rates, which can be translated into upper limits on the expected γ-ray flux not only from the inner Galaxy, but also from any substructure, with or without stars, in the Galactic halo.

Extracting the Galactic Center excess’ source-count distribution with neural nets

The two leading hypotheses for the Galactic Center Excess (GCE) in the $\textit{Fermi}$ data are an unresolved population of faint millisecond pulsars (MSPs) and dark-matter (DM) annihilation. The dichotomy between these explanations is typically reflected by modeling them as two separate emission components. However, point-sources (PSs) such as MSPs become statistically degenerate with smooth Poisson emission in the ultra-faint limit (formally where each source is expected to contribute much less than one photon on average), leading to an ambiguity that can render questions such as whether the emission is PS-like or Poissonian in nature ill-defined. We present a conceptually new approach that describes the PS and Poisson emission in a unified manner and only afterwards derives constraints on the Poissonian component from the so obtained results. For the implementation of this approach, we leverage deep learning techniques, centered around a neural network-based method for histogram regression that expresses uncertainties in terms of quantiles. We demonstrate that our method is robust against a number of systematics that have plagued previous approaches, in particular DM / PS misattribution. In the $\textit{Fermi}$ data, we find a faint GCE described by a median source-count distribution (SCD) peaked at a flux of $\sim4 \times 10^{-11} \ \text{counts} \ \text{cm}^{-2} \ \text{s}^{-1}$ (corresponding to $\sim3 - 4$ expected counts per PS), which would require $N \sim \mathcal{O}(10^4)$ sources to explain the entire excess (median value $N = \text{29,300}$ across the sky). Although faint, this SCD allows us to derive the constraint $\eta_P \leq 66\%$ for the Poissonian fraction of the GCE flux $\eta_P$ at 95% confidence, suggesting that a substantial amount of the GCE flux is due to PSs.

Dissecting the inner Galaxy with gamma-ray pixel count statistics

The nature of the GeV gamma-ray Galactic center excess (GCE) in the data of Fermi-LAT is still under investigation. Different techniques, such as template fitting and photon-count statistical methods, have been applied in the past few years in order to disentangle between a GCE coming from sub-threshold point sources or rather from diffuse emissions, such as the dark matter annihilation in the Galactic halo. A major limit to all these studies is the modeling of the Galactic diffuse foreground, and the impact of residual mis-modeled emission on the results’ robustness. We combine for the first time adaptive template fitting and pixel count statistical methods in order to assess the role of sub-threshold point sources to the GCE, while minimizing the mis-modelling of diffuse emission components. We reconstruct the flux distribution of point sources in the inner Galaxy well below the Fermi-LAT detection threshold, and measure their radial and longitudinal profiles. We find that point sources and diffuse emission from the Galactic bulge each contributes about 10% of the total emission therein, disclosing a potential sub-threshold point-source contribution to the GCE.

Dissecting the Inner Galaxy with γ-Ray Pixel Count Statistics.

We combine adaptive template fitting and pixel count statistics in order to assess the nature of the Galactic Center excess in Fermi-LAT data. We reconstruct the flux distribution of point sources well below the Fermi-LAT detection threshold, and measure their radial and longitudinal profiles in the inner Galaxy. We find that all point sources and the bulge-correlated diffuse emission each contributes O(10%) of the total inner Galaxy emission, and disclose a potential subthreshold point-source contribution to the Galactic Center excess.

The effect of kick velocities on the spatial distribution of millisecond pulsars and implications for the Galactic center excess

Recently it has become apparent that the Galactic center excess (GCE) is spatially correlated with the stellar distribution in the Galactic bulge. This has given extra motivation for the unresolved population of millisecond pulsars (MSPs) explanation for the GCE. However, in the “recycling” channel the neutron star forms from a core collapse supernovae that undergoes a random “kick” due to the asymmetry of the explosion. This would imply a smoothing out of the spatial distribution of the MSPs. We use N-body simulations to model how the MSP spatial distribution changes. We estimate the probability distribution of natal kick velocitiesusing the resolved gamma-ray MSP proper motions, where MSPs have random velocities relative to the circular motion with a scale parameter of 77 ± 6 km/s. We find that, due to the natal kicks, there is an approximately 10% increase in each of the bulge MSP spatial distribution dimensions and also the bulge MSP distribution becomes less boxy but is still far from being spherical.