Spatial variations of the average properties that characterize the hadronic component of the diffuse Galactic cosmic-ray sea, in particular the spectral slope and normalization, may unveil critical information about their confinement mechanism in the Galaxy. In the first part of this paper we perform an analysis of the Fermi-LAT gamma-ray data with the SkyFACT package, which combines image reconstruction techniques with standard template fitting, isolate the hadronic emission and decompose it into Galactocentric rings. We find a significant hardening of the hadronic spectral index towards the molecular ring. We study this hardening in different energy ranges, and assess its resilience with respect to different prescriptions in the analysis setup. In the second part we quantify the contribution to the diffuse gamma-ray flux coming from unresolved point sources with a dedicated Monte Carlo simulation, and consider whether the trend characterized in the first part can be mimicked by a progressively more relevant flux associated to this component in the inner Galaxy. We find that the observed hardening of the hadronic spectral index cannot be due to unresolved sources in the sub-TeV energy range, especially outside the molecular ring, given reasonable assumptions about the unresolved source population.

Radio relics, elongated, non-thermal, structures located at the edges of galaxy clusters, are the result of synchrotron radiation from cosmic-ray electrons accelerated by merger-driven shocks at the cluster outskirts. However, X-ray observations of such shocks in some clusters suggest that they are too weak to efficiently accelerate electrons via diffusive shock acceleration to energies required to produce the observed radio power. We examine this issue in the merging galaxy cluster Abell 3667 (A3667), which hosts a pair of radio relics. While the Northwest relic in A3667 has been well studied in the radio and X-ray by multiple instruments, the Southeast relic region has only been observed so far by Suzaku, which detected a temperature jump across the relic, suggesting the presence of a weak shock. We present observations of the Southeastern region of A3667 with XMM-Newton centered on the radio relic. We confirm the existence of an X-ray shock with Mach number of about 1.8 from a clear detection of temperature jump and a tentative detection of a density jump, consistent with previous measurements by Suzaku. We discuss the implications of this measurement for diffusive shock acceleration as the main mechanism for explaining the origin of radio relics. We then speculate on the plausibility of alternative scenarios, including re-acceleration and variations in the Mach number along shock fronts.

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.

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.

Dark matter (DM) particles are predicted by several well motivated models to yield Standard Model particles through self-annihilation that can potentially be detected by astrophysical observations. In particular, the production of charged particles from DM annihilation in astrophysical systems that contain magnetic fields yields radio emission through synchrotron radiation and X-ray emission through inverse Compton scattering of ambient photons. We introduce RX-DMFIT, a tool used for calculating the expected secondary emission from DM annihilation. RX-DMFIT includes a wide range of customizable astrophysical and particle parameters and incorporates important astrophysics including the diffusion of charged particles, relevant radiative energy losses, and magnetic field modelling. We demonstrate the use and versatility of RX-DMFIT by analyzing the potential radio and X-ray signals for a variety of DM particle models and astrophysical environments including galaxy clusters, dwarf spheroidal galaxies and normal galaxies. We then apply RX-DMFIT to a concrete example using Segue I radio data to place constraints for a range of assumed DM annihilation channels. For WIMP models with $M_{\chi} \leq 100$ GeV and assuming weak diffusion, we find that the the leptonic $\mu^+\mu^-$ and $\tau^+\tau^-$ final states provide the strongest constraints, placing limits on the DM particle cross-section well below the thermal relic cross-section, while even for the $b\bar{b}$ channel we find limits close to the thermal relic cross-section. Our analysis shows that radio emission provides a highly competitive avenue for dark matter searches.

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.

The annihilation of dark matter particles is expected to yield a broad radiation spectrum via the production of Standard Model particles in astrophysical environments. In particular, electrons and positrons from dark matter annihilation produce synchrotron radiation in the presence of magnetic fields. Galaxy clusters are the most massive collapsed structures in the universe, and are known to host $\sim\mu$G-scale magnetic fields. They are therefore ideal targets to search for, or to constrain the synchrotron signal from dark matter annihilation. In this work we use the expected sensitivities of several planned surveys from the next generation of radio telescopes to predict the constraints on dark matter annihilation models which will be achieved in the case of non-detections of diffuse radio emission from galaxy clusters. Specifically, we consider the Tier 1 survey planned for the Low Frequency Array (LOFAR) at 120 MHz, the EMU survey planned for the Australian Square Kilometre Array Pathfinder (ASKAP) at 1.4 GHz, and planned surveys for APERTIF at 1.4 GHz. We find that, for massive clusters and dark matter masses $\lesssim 100$ GeV, the predicted limits on the annihilation cross section would rule out vanilla thermal relic models for even the shallow LOFAR Tier 1, ASKAP, and APERTIF surveys.

Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of $\gamma$-ray telescopes. Here we use three years of Fermi Large Area Telescope (LAT) data, which are the most suitable for searching for very extended emission in the vicinity of nearby Virgo galaxy cluster. Our analysis reveals statistically significant extended emission which can be well characterized by a uniformly emitting disk profile with a radius of 3\deg that moreover is offset from the cluster center. We demonstrate that the significance of this extended emission strongly depends on the adopted interstellar emission model (IEM) and is most likely an artifact of our incomplete description of the IEM in this region. We also search for and find new point source candidates in the region. We then derive conservative upper limits on the velocity-averaged DM pair annihilation cross section from Virgo. We take into account the potential $\gamma$-ray flux enhancement due to DM sub-halos and its complex morphology as a merging cluster. For DM annihilating into $b\overline{b}$, assuming a conservative sub-halo model setup, we find limits that are between 1 and 1.5 orders of magnitude above the expectation from the thermal cross section for $m_{\mathrm{DM}}\lesssim100\,\mathrm{GeV}$. In a more optimistic scenario, we exclude $\langle \sigma v \rangle\sim3\times10^{-26}\,\mathrm{cm^{3}\,s^{-1}}$ for $m_{\mathrm{DM}}\lesssim40\,\mathrm{GeV}$ for the same channel. Finally, we derive upper limits on the $\gamma$-ray-flux produced by hadronic cosmic-ray interactions in the inter cluster medium. We find that the volume-averaged cosmic-ray-to-thermal pressure ratio is less than $\sim6\%$.

A2319 is a massive, merging galaxy cluster with a previously detected radio halo that roughly follows the X-ray emitting gas. We present the results from recent observations of A2319 at 20 cm with the Jansky Very Large Array (VLA) and a re-analysis of the X-ray observations from XMM-Newton, to investigate the interactions between the thermal and nonthermal components of the ICM . We confirm previous reports of an X-ray cold front, and report on the discovery of a distinct core to the radio halo, 800 kpc in extent, that is strikingly similar in morphology to the X-ray emission, and drops sharply in brightness at the cold front. We detect additional radio emission trailing off from the core, which blends smoothly into the 2 Mpc halo detected with the Green Bank Telescope (GBT; Farnsworth et al., 2013). We speculate on the possible mechanisms for such a two-component radio halo, with sloshing playing a dominant role in the core. By directly comparing the X-ray and radio emission, we find that a hadronic origin for the cosmic ray electrons responsible for the radio halo would require a magnetic field and/or cosmic ray proton distribution that increases with radial distance from the cluster center, and is therefore disfavored.

Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $\gamma$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $\gamma$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7\sigma$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $\gamma$-ray flux from individual clusters in our sample.

Annihilation of dark matter can result in the production of stable Standard Model particles including electrons and positrons that, in the presence of magnetic fields, lose energy via synchrotron radiation, observable as radio emission. Galaxy clusters are excellent targets to search for or to constrain the rate of dark matter annihilation, as they are both massive and dark matter dominated. In this study, we place limits on dark matter annihilation in a sample of nearby clusters using upper limits on the diffuse radio emission, low levels of observed diffuse emission, or detections of radio mini-haloes. We find that the strongest limits on the annihilation cross section are better than limits derived from the non-detection of clusters in the gamma-ray band by a factor of approximately 3 or more when the same annihilation channel and subtructure model, but different best-case clusters, are compared. The limits on the cross section depend on the assumed amount of substructure, varying by as much as 2 orders of magnitude for increasingly optimistic substructure models as compared to a smooth NFW profile. In our most optimistic case, using the results of the Phoenix Project (Gao et al. 2012b) we find that the derived limits reach below the thermal relic cross section of 3x10^-26 cm^3 s^-1 for dark matter masses as large as 400 GeV, for the bbar annihilation channel. We discuss uncertainties due to the limited available data on the magnetic field structure of individual clusters. We also report the discovery of diffuse radio emission from the central 30-40 kpc regions of the groups M49 and NGC4636.

Star formation in galaxies is observed to be associated with gamma-ray emission. The detection of gamma rays from star-forming galaxies by the Fermi Large Area Telescope (LAT) has allowed the determination of a functional relationship between star formation rate and gamma-ray luminosity (Ackermann et. al. 2012). Since star formation is known to scale with total infrared (8-1000 micrometers) and radio (1.4 GHz) luminosity, the observed infrared and radio emission from a star-forming galaxy can be used to quantitatively infer the galaxy's gamma-ray luminosity. Similarly, star forming galaxies within galaxy clusters allow us to derive lower limits on the gamma-ray emission from clusters, which have not yet been conclusively detected in gamma rays. In this study we apply the relationships between gamma-ray luminosity and radio and IR luminosities derived in Ackermann et. al. 2012 to a sample of galaxy clusters from Ackermann et. al. 2010 in order to place lower limits on the gamma-ray emission associated with star formation in galaxy clusters. We find that several clusters have predicted lower limits on gamma-ray emission that are within an order of magnitude of the upper limits derived in Ackermann et. al. 2010 based on non-detection by Fermi-LAT. Given the current gamma-ray limits, star formation likely plays a significant role in the gamma-ray emission in some clusters, especially those with cool cores. We predict that both Fermi-LAT over the course of its lifetime and the future Cherenkov Telescope Array will be able to detect gamma-ray emission from star-forming galaxies in clusters.