We present a fully analytical, time-dependent leptonic one-zone model that describes a simplified radiation process of multiple interacting ultrarelativistic electron populations, accounting for the flaring of GeV blazars. In this model, several mono-energetic, ultrarelativistic electron populations are successively and instantaneously injected into the emission region, i.e., a magnetized plasmoid propagating along the blazar jet, and subjected to linear, time-independent synchrotron radiative losses, which are caused by a constant magnetic field, and nonlinear, time-dependent synchrotron self-Compton radiative losses in the Thomson limit. Considering a general multiple-injection scenario is, from a physical point of view, more realistic than the usual single-injection scenario invoked in common blazar models, as blazar jets may extend over tens of kiloparsecs and, thus, most likely pick up several particle populations from intermediate clouds. We analytically compute the electron number density by solving a kinetic equation using Laplace transformations and the method of matched asymptotic expansions. Moreover, we explicitly calculate the optically thin synchrotron intensity, the synchrotron self-Compton intensity in the Thomson limit, as well as the associated total fluences. In order to mimic injections of finite duration times and radiative transport, we model flares by sequences of these instantaneous injections, suitably distributed over the entire emission region. Finally, we present a parameter study for the total synchrotron and synchrotron self-Compton fluence spectral energy distributions for a generic three-injection scenario, varying the magnetic field strength, the Doppler factor, and the initial electron energy of the first injection in realistic parameter domains, demonstrating that our model can reproduce the typical broad-band behavior seen in observational data.

It is widely believed that supernova remnants are the best candidate sources for the observed cosmic ray flux up to the knee, i.e. up to ~PeV energies. Indeed, the gamma-ray spectra of some supernova remnants can be well explained by assuming the decay of neutral pions which are created in hadronic interactions. Therefore, fitting the corresponding gamma spectra allows us to derive the spectra of cosmic rays at the source which are locally injected into our Galaxy. Using these spectra as a starting point, we propagate the cosmic rays through the Galaxy using the publicly available GALPROP code. Here, we will present first results on the contribution of those SNRs to the total cosmic ray flux and discuss implications.

Context: Since a few years, signatures of supernova remnants associated with molecular clouds have been detected in gamma rays. Whether these gamma rays are generated by cosmic ray electrons or by cosmic ray protons is usually not known. The detection of hadronic ionization signatures in spatial coincidence with gamma ray signatures can help to unambiguously identify supernova remnants as sources of cosmic ray protons. Methods: In order to calculate hadronic signatures from cosmic ray-induced ionization for an examination of the origin of the observed gamma rays, the transport equation for cosmic ray protons propagating in a molecular cloud, including the relevant momentum loss processes, is solved analytically and the proton flux at any position in the cloud is determined. Results: Since the solution of the transport equation is obtained for arbitrary source functions, it can be used for a variety of supernova remnants. The corresponding theoretical ionization rate, as a function of the penetration depth, is derived and compared to photoinduced ionization profiles in a case study with four supernova remnants associated with molecular clouds. Three of the remnants show a clear dominance of the hadronically induced ionization rate, while for one remnant, X-ray emission seems to dominate by a factor of 10. Conclusions: This is the first derivation of position-dependent profiles for cosmic ray-induced ionization with an analytic solution for arbitrary cosmic ray source spectra. The cosmic ray-induced ionization has to be compared to X-ray ionization for strong X-ray sources. For sources dominated by cosmic ray-induced ionization (e.g., W49B), the ionization profiles can be used in the future to map the spatial structure of hadronic gamma rays and rotation-vibrational lines induced by cosmic ray protons, helping to identify sources of hadronic cosmic rays.

Astrospheres and wind bubbles of massive stars are believed to be sources of cosmic rays with energies $E\lesssim 1\,$TeV. These particles are not directly detectable, but their impact on surrounding matter, in particular ionisation of atomic and molecular hydrogen, can lead to observable signatures. A correlation study of both gamma ray emission, induced by proton-proton interactions of cosmic ray protons with kinetic energies $E_\mathrm{p}\ge 280\,$MeV with ambient hydrogen, and ionisation induced by cosmic ray protons of kinetic energies $E_\mathrm{p}< 280\,$MeV can be performed in order to study potential sources of (sub)TeV cosmic rays.

For star-forming regions, there is a correlation of radio and FIR-emission established. The radio emission is caused by synchrotron radiation of electrons, while the FIR emission is attributed to HII regions of OB stars and hot dust powered by those stars. Another possible correlation for star-forming regions might exist between the aforementioned radio emission and the gamma ray emission in the GeV regime. The GeV gamma ray emission of star-forming regions is dominated by the decay of neutral pions formed in collisions of cosmic ray (CR) protons accelerated at supernova remnants (SNRs) with ambient hydrogen, while the electrons generating the synchrotron emission are assumed to be accelerated at the same SNRs. Assuming the same spectral shape for the CR proton and electron distribution and a fixed ratio of CR protons to electrons, the proton- and electron-calorimetry of the star-forming regions can be tested. Furthermore, this provides a method to derive the magnetic field strength in the star-forming region complementary to the existing methods.

Energetic gamma rays (GeV to TeV photon energy) have been detected toward several supernova remnants (SNR) that are associated with molecular clouds. If the gamma rays are produced mainly by hadronic processes rather than leptonic processes like bremsstrahlung, then the flux of energetic cosmic ray nuclei (>1 GeV) required to produce the gamma rays can be inferred at the site where the particles are accelerated in SNR shocks. It is of great interest to understand the acceleration of the cosmic rays of lower energy (<1 GeV) that accompany the energetic component. These particles of lower energy are most effective in ionizing interstellar gas, which leaves an observable imprint on the interstellar ion chemistry. A correlation of energetic gamma radiation with enhanced interstellar ionization can thus be used to support the hadronic origin of the gamma rays and to constrain the acceleration of ionizing cosmic rays in SNR. Using observational gamma ray data, the primary cosmic ray proton spectrum can be modeled for E >1 GeV, and careful extrapolation of the spectrum to lower energies offers a method to calculate the ionization rate of the molecular cloud.

Energetic gamma rays (GeV to TeV photon energy) have been detected toward several supernova remnants (SNR) associated with molecular clouds. If the gamma rays are produced mainly by hadronic processes rather than leptonic processes like bremsstrahlung, then the flux of energetic cosmic ray (CR) nuclei (>1 GeV) required to produce the gamma rays can be inferred at the site where the particles are accelerated in SNR shocks. It is of great interest to understand the acceleration of the CR of lower energy (<1 GeV) accompanying the energetic component. These particles of lower energy are most effective in ionizing interstellar gas, leaving an observable imprint on the interstellar ion chemistry. A correlation of energetic gamma radiation with enhanced interstellar ionization can thus support the hadronic origin of the gamma rays and constrain the acceleration of ionizing CR in SNR. We propose a method to test the hadronic origin of GeV gamma rays from SNR associated with a molecular cloud. We use observational gamma ray data for each of these SNR known, modeling the observations to obtain the underlying proton spectrum assuming that the gamma rays are produced by pion decay. Assuming that the acceleration mechanism does not only produce high energy protons, but also low energy protons, this proton spectrum at the source is then used to calculate the ionization rate of the molecular cloud. Ionized molecular hydrogen triggers a chemical network forming molecular ions. The relaxation of these ions results in characteristic line emission, which can be predicted. We show that the ionization rate for at least two objects is more than an order of magnitude above Galactic average for molecular clouds, hinting at an enhanced formation rate of molecular ions. There will be interesting opportunities to measure crucial molecular ions in the infrared and submillimeter-wave parts of the spectrum.

During the past few years, first observations of starburst galaxies at >GeV energies could be made with the Fermi Gamma-ray Space Telescope (GeV range) and Imaging Air Cherenkov Telescopes (TeV range). The two nearest starbursts, M82 and NGC253 were detected, and most recently, the detection of two starburst-Seyfert composites (NGC1068 and NGC4945) were reported. The emission for the two starbursts is best explained by hadronic interactions, and thus providing a first, unique opportunity to study the role of cosmic rays in galaxies. In this paper, the role of cosmic rays for the non-thermal component of galaxies is reviewed by discussing the entire non-thermal frequency range from radio emission to TeV energies. In particular, the interpretation of radio emission arising from electron synchrotron radiation is predicted to be correlated to TeV emission coming from interactions of accelerated hadrons. This is observed for the few objects known at TeV energies, but the correlation needs to be established with significantly higher statistics. An outlook of the possibility of tracing cosmic rays with molecular ions is given.

The rate of ionization by cosmic rays in interstellar gas directly associated with gamma-ray emitting supernova remnants is for the first time calculated to be several orders of magnitude larger than the Galactic average. Analysis of ionization-induced chemistry yields the first quantitative prediction of the astrophysical H_{2}^{+} emission line spectrum, which should be detectable together with H_3^+ lines. The predicted coincident observation of those emission lines and gamma-rays will help prove that supernova remnants are sources of cosmic rays.