Gamma-ray observations have revealed strong variability in blazar luminosities in the gamma-ray band over time scales as short as minutes. We show, for the first time, that the correlation of the spectrum with intensity is consistent with the behavior of the luminosity variation of blazar SEDs along a blazar sequence for low synchrotron peak blazars. We show that the observational signatures of variability with flux are consistent with wakefield acceleration of electrons initiated by instabilities in the blazar accretion disk. This mechanism reproduces the observed time variations as short as 100 seconds. The wakefield mechanism also predicts a reduction of the electron spectral index with increased gamma-ray luminosity, which could be detected in higher energy observations well above the inverse Compton peak.

We test the sensitivity of neutrino parameter constraints from combinations of CMB and LSS data sets to the assumed form of the primordial power spectrum (PPS) using Bayesian model selection. Significantly, none of the tested combinations, including recent high-precision local measurements of $\mathrm{H}_0$ and cluster abundances, indicate a signal for massive neutrinos or extra relativistic degrees of freedom. For PPS models with a large, but fixed number of degrees of freedom, neutrino parameter constraints do not change significantly if the location of any features in the PPS are allowed to vary, although neutrino constraints are more sensitive to PPS features if they are known a priori to exist at fixed intervals in $\log k$. Although there is no support for a non-standard neutrino sector from constraints on both neutrino mass and relativistic energy density, we see surprisingly strong evidence for features in the PPS when it is constrained with data from Planck 2015, SZ cluster counts, and recent high-precision local measurements of $\mathrm{H}_0$. Conversely combining Planck with matter power spectrum and BAO measurements yields a much weaker constraint. Given that this result is sensitive to the choice of data this tension between SZ cluster counts, Planck and $\mathrm{H}_0$ measurements is likely an indication of unmodeled systematic bias that mimics PPS features, rather than new physics in the PPS or neutrino sector.

We present a new extended gamma ray excess detected with the Fermi Satellite Large Area Telescope toward the Galactic Center that traces the morphology of infrared starlight emission. Combined with its measured spectrum, this new extended source is approximately consistent with inverse Compton emission from a high-energy electron-positron population with energies up to about 10 GeV. Previously detected emissions tracing the 20 cm radio, interpreted as bremsstrahlung radiation, and the Galactic Center Extended emission tracing a spherical distribution and peaking at 2 GeV, are also detected. We show that the inverse Compton and bremsstrahlung emissions are likely due to the same source of electrons and positrons. All three extended emissions may be explained within the framework of a model where the dark matter annihilates to leptons or a model with unresolved millisecond pulsars in the Galactic Center.

We construct empirical models of the diffuse gamma-ray background toward the Galactic Center. Including all known point sources and a template of emission associated with interactions of cosmic rays with molecular gas, we show that the extended emission observed previously in the Fermi Large Area Telescope data toward the Galactic Center is detected at high significance for all permutations of the diffuse model components. However, we find that the fluxes and spectra of the sources in our model change significantly depending on the background model. In particular, the spectrum of the central Sgr A$^\ast$ source is less steep than in previous works and the recovered spectrum of the extended emission has large systematic uncertainties, especially at lower energies. If the extended emission is interpreted to be due to dark matter annihilation, we find annihilation into pure $b$-quark and $\tau$-lepton channels to be statistically equivalent goodness of fits. In the case of the pure $b$-quark channel, we find a dark matter mass of $39.4\left(^{+3.7}_{-2.9}\rm\ stat.\right)\left(\pm 7.9\rm\ sys.\right)\rm\ GeV$, while a pure $\tau^{+} \tau^{-}$-channel case has an estimated dark matter mass of $9.43\left(^{+0.63}_{-0.52}\rm\ stat.\right)(\pm 1.2\rm\ sys.)\ GeV$. Alternatively, if the extended emission is interpreted to be astrophysical in origin such as due to unresolved millisecond pulsars, we obtain strong bounds on dark matter annihilation, although systematic uncertainties due to the dependence on the background models are significant.