# Ionization Rate for Cosmic Rays, given the CR spectrum

Hi all, so I am working on a presentation. What I have done so far was to calculate the cosmic ray spectrum (so given some data from Voyager 1 and PAMELA, I found the intensity as a function of the energy, in units of m^-2 (sr s MeV)^-1.

Given this function, which is essentially J = constant*E^d[(E^k+const)/(const)]^(const)

This was done for 3 species, including helium, electrons, and protons (so the function is similar but different for all 3).

note: d, k are also constants. These were figured out in Matlab using least squares regression, but my question isn't related to this, so I didn't feel it necessary to include the entire equation typed out. If needed, I can include this.

Next, I found the energy density (aka kinetic energy density) by integrating this intensity over a solid angle. So in essence, integrate the above function J, and that'll be times the constant 4pi which comes out of the solid angle.I assumed the arrival of cosmic rays to be isotropic.

Great. So now I need to find the ionization rate. And I have no idea how to do this. So given some cosmic ray spectrum that has J proportional to some value of E, and/or given the energy density (by integrating the spectrum), how on earth do I go about finding the ionization rate?

Thanks!

## Answers and Replies

mfb
Mentor
There are formulas for the ionization rate as function of the energy. For an isolated piece of material in free space, you can multiply the flux by the ionization rate, then integrate over the energy. Realistically, your satellite will provide shielding against low-energetic particles from some sides, and produce showers out of high-energetic particles. A proper analysis will need some finite element simulation of the propagation of cosmic rays of different energies from different direction.