A fundamental question regarding the Galactic Center Excess (GCE) is whether the underlying structure is point-like or smooth. This debate, often framed in terms of a millisecond pulsar or annihilating dark matter (DM) origin for the emission, awaits a conclusive resolution. In this work we weigh in on the problem using Bayesian graph convolutional neural networks. In simulated data, our neural network (NN) is able to reconstruct the flux of inner galaxy emission components to on average $\sim$0.5%, comparable to the non-Poissonian template fit (NPTF). When applied to the actual $\textit{Fermi}$-LAT data, we find that the NN estimates for the flux fractions from the background templates are consistent with the NPTF; however, the GCE is almost entirely attributed to smooth emission. While suggestive, we do not claim a definitive resolution for the GCE, as the NN tends to underestimate the flux of point-sources peaked near the 1$\sigma$ detection threshold. Yet the technique displays robustness to a number of systematics, including reconstructing injected DM, diffuse mismodeling, and unmodeled north-south asymmetries. So while the NN is hinting at a smooth origin for the GCE at present, with further refinements we argue that Bayesian Deep Learning is well placed to resolve this DM mystery.

In a model with an extra $U(1)$ gauge to SM gauge group, we have shown the allowed region of masses of extra gauge boson and the dark matter which is the lightest one among other right-handed Majorana fermions present in the model. To obtain this region, we have used bounds coming from constraints on active-sterile neutrino masses and mixing from various oscillation experiments, constraint on dark matter relic density obtained by PLANCK together with the constraint on the extra gauge boson mass and its gauge coupling recently obtained by ATLAS Collaboration at LHC. From the allowed regions, it is possible to get some lower bounds on the masses of the extra gauge boson and the dark matter and considering those values it is possible to infer what could be the spontaneous symmetry breaking scale of an extra $U(1)$ gauge symmetry.

Traditionally, incorporating additional physics into existing cosmological simulations requires re-running the cosmological simulation code, which can be computationally expensive. We show that conditional Generative Adversarial Networks (cGANs) can be harnessed to predict how changing the underlying physics alters the simulation results. To illustrate this, we train a cGAN to learn the impact of dark matter annihilation feedback (DMAF) on the gas density distribution. The predicted gas density slices are visually difficult to distinguish from their real brethren and the peak counts differ by less than 10 per cent for all test samples (the average deviation is < 3 per cent). Finally, we invert the problem and show that cGANs are capable of endowing smooth density distributions with realistic substructure. The cGAN does however have difficulty generating new knots as well as creating/eliminating bubble-like structures. We conclude that trained cGANs can be an effective approach to provide mock samples of cosmological simulations incorporating DMAF physics from existing samples of standard cosmological simulations of the evolution of cosmic structure.