Deciphering Residual Emissions: Time-dependent Models for the Nonthermal Interstellar Radiation from the Milky Way

Cosmic rays (CRs) in the Galaxy are an important dynamical component of the interstellar medium (ISM) that interact with the other major components (interstellar gas and magnetic and radiation fields) to produce broadband interstellar emissions that span the electromagnetic spectrum. The standard modeling of CR propagation and production of the associated emissions is based on a steady-state assumption, where the CR source spatial density is described using a smoothly varying function of position that does not evolve with time. While this is a convenient approximation, reality is otherwise, where primary CRs are produced in and about highly localized regions, e.g., supernova remnants, which have finite lifetimes. In this paper, we use the latest version of the galprop CR propagation code to model time-dependent CR injection and propagation through the ISM from a realistic 3D discretized CR source density distribution, together with full 3D models for the other major ISM components, and make predictions of the associated broadband nonthermal emissions. We compare the predictions for the discretized and equivalent steady-state model, finding that the former predicts novel features in the broadband nonthermal emissions that are absent for the steady-state case. Some of the features predicted by the discretized model may be observable in all-sky observations made by WMAP and Planck, the recently launched eROSITA, the Fermi-LAT, and ground-based observations by HESS, HAWC, and the forthcoming CTA. The nonthermal emissions predicted by the discretized model may also provide explanations of puzzling anomalies in high-energy γ-ray data, such as the Fermi-LAT north/south asymmetry and residuals like the so-called "Fermi bubbles."