Cosmic rays are an essential ingredient in the evolution of the interstellar medium, as they dominate the ionisation of the dense molecular gas, where stars and planets form. However, since they are efficiently scattered by the galactic magnetic fields, many questions remain open, such as where exactly they are accelerated, what is their original energy spectrum, and how they propagate into molecular clouds. In this work we present new observations and discuss in detail a method that allows us to measure the cosmic ray ionisation rate towards the molecular clouds close to the W28 supernova remnant. To perform these measurements, we use CO, HCO$^+$, and DCO$^+$ millimetre line observations and compare them with the predictions of radiative transfer and chemical models away from thermodynamical equilibrium. The CO observations allow us to constrain the density, temperature, and column density towards each observed position, while the DCO$^+$/HCO$^+$ abundance ratios provide us with constraints on the electron fraction and, consequently, on the cosmic ray ionisation rate. Towards positions located close to the supernova remnant, we find cosmic ray ionisation rates much larger ($\gtrsim100$) than those in standard galactic clouds. Conversely, towards one position situated at a larger distance, we derive a standard cosmic ray ionisation rate. Overall, these observations support the hypothesis that the $\gamma$ rays observed in the region have a hadronic origin. In addition, based on CR diffusion estimates, we find that the ionisation of the gas is likely due to $0.1 -1$ GeV cosmic rays. Finally, these observations are also in agreement with the global picture of cosmic ray diffusion, in which the low-energy tail of the cosmic ray population diffuses at smaller distances than the high-energy counterpart.

The region towards W51C is a convincing example of interaction between a supernova remnant and a surrounding molecular cloud. Large electron abundances have been reported towards the position W51C-E located in this interaction region, and it was proposed that the enhanced ionization fraction was due to cosmic ray particles freshly accelerated by the SNR shock. We present PdB interferometer observations of the H$^{13}$CO$^+$(1-0) and DCO$^+$(2-1) emission lines centered at position W51C-E. These observations confirm the previous scenario of cosmic-ray induced ionization at this location. In addition, SiO(2-1) emission has been successfully mapped in the close vicinity of W51C-E, with a spatial resolution of 7". The morphology and kinematics of the SiO emission are analyzed and strongly suggest that this emission is produced by the passage of the SNR primary shock. Put in conjunction with the enhanced ionization fraction in this region, we give a consistent picture in which the W51C-E position is located downstream of the shock, where a large reservoir of freshly accelerated particles is available.