Molecular clouds act as targets for cosmic rays (CR), revealing their presence through either gamma-ray emission due to proton-proton interactions, and/or through the ionization level in the cloud, produced by the CR flux. The ionization rate is a un
ique tool, to some extent complementary to the gamma-ray emission, in that it allows to constrain the CR spectrum especially for energies below the pion production rate ($approx 280$ MeV). Here we study the effect of ionization on $H_2$ clouds due to both CR protons and electrons, using the fully relativistic ionization cross sections, which is important to correctly account for the contribution due to relativistic CRs. The contribution to ionization due to secondary electrons is also included self-consistently. The whole calculation has been implemented into a numerical code which is publicly accessible through a web-interface. The code also include the calculation of gamma-ray emission once the CR spectrum
The W51 complex hosts the supernova remnant W51C which is known to interact with the molecular clouds in the star forming region W51B. In addition, a possible pulsar wind nebula CXO J192318.5+140305 was found likely associated with the supernova remn
ant. Gamma-ray emission from this region was discovered by Fermi/LAT (between 0.2 and 50 GeV) and H.E.S.S. (>1 TeV). The spatial distribution of the events could not be used to pinpoint the location of the emission among the pulsar wind nebula, the supernova remnant shell and/or the molecular cloud. However, the modeling of the spectral energy distribution presented by the Fermi/LAT collaboration suggests a hadronic emission mechanism. We performed observations of the W51 complex with the MAGIC telescopes for more than 50 hours. The good angular resolution in the medium (few hundred GeV) to high (above 1 TeV) energies allow us to perform morphological studies. We detect an extended emission of very-high-energy gamma rays, with a significance of 11 standard deviations. We extend the spectrum from the highest Fermi/LAT energies to sim 5 TeV and find that it follows a single power law with an index of 2.58 pm 0.07stat pm 0.22syst . The main part of the emission coincides with the shocked cloud region, while we find a feature extending towards the pulsar wind nebula. The possible contribution of the pulsar wind nebula, assuming a point-like source, shows no dependence on energy and it is about 20% of the overall emission. The broad band spectral energy distribution can be explained with a hadronic model that implies proton acceleration above 100 TeV. This result, together with the morphology of the source, tentatively suggests that we observe ongoing acceleration of ions in the interaction zone between supernova remnant and cloud. These results shed light on the long-standing problem of the origin of galactic cosmic rays.
Located in a dense complex environment, W51C provides an excellent scenario to probe accelerated protons in SNRs and their interaction with surrounding target material. Here we report the observation of extended Very High Energy (VHE) gamma-ray emiss
ion from the W51C supernova remnant (SNR) with MAGIC. Detections of extended gamma-ray emission in the same region have already been reported by the Fermi and H.E.S.S. collaborations. Fermi/LAT measured the source spectrum in the energy range between 0.2 and 50 GeV, which was found to be well fit by a hadronic neutral-pion decay model. The VHE observations presented here, obtained with the improved MAGIC stereo system, allow us to pinpoint the VHE gamma-ray emission in the dense shocked molecular cloud surrounding the remnant shell. The MAGIC data also allow us to measure, for the first time, the VHE emission spectrum of W51C from the highest Fermi/LAT energies up to TeV. The spatial distribution and spectral properties of the VHE emission suggest a hadronic origin of the observed gamma rays. Therefore W51C is a prime candidate for a cosmic ray accelerator.
