Understanding uncertainties in modeling the galactic diffuse gamma-ray emission

The nature of the Galactic diffuse gamma-ray emission as measured by the Fermi Gamma-ray Space Telescope has remained an active area of research for the last several years. A standard technique to disentangle the origins of the diffuse emission is the template fitting approach, where predictions for various diffuse components, such as emission from cosmic rays derived from Galprop or Dragon, are compared to the data. However, this method always results in an overall bad fit to the data, with strong residuals that are difficult to interpret. Additionally, there are instrinsic uncertainties in the predicted templates that are not accounted for naturally with this method. We therefore introduce a new template fitting approach to study the various components of the Galactic diffuse gamma-ray emission, and their correlations and uncertainties. We call this approach Sky Factorization with Adaptive Constrained Templates (SkyFACT). Rather than using fixed predictions from cosmic-ray propagation codes and examining the residuals to evaluate the quality of fits and the presence of excesses, we introduce additional fine-grained variations in the templates that account for uncertainties in the predictions, such as uncertainties in the gas tracers and from small scale variations in the density of cosmic rays. We show that fits to the gamma-ray diffuse emission can be dramatically improved by including an appropriate level of uncertainty in the initial spatial templates from cosmic-ray propagation codes. We further show that we can recover the morphology of the Fermi Bubbles from its spectrum alone with SkyFACT.