Ultra-high-energy cosmic rays in a galactic wind and its termination shock

Results are reported from numerical modeling of the acceleration and transport of ultra-high-energy cosmic rays in a galactic wind and its termination shock. A two-dimensional (azimuthally symmetric) wind and spiral magnetic field, with a spherical termination shock, where the velocity drops suddenly, is assumed. The time-dependent cosmic-ray transport equation, including all major transport effects is solved using an implicit finite-difference scheme. Particles are injected as the shock of low energy, and the subsequent evolution of the distribution function is followed. Iron nuclei are readily accelerated at the shock to energies up to 100 billion GeV, and protons to 10 billion GeV. A major effect aiding the acceleration of these particles is the spiral of the magnetic field carried out by the wind, caused by the rotation of the Galaxy, with the result that the shock is nearly normal over most of its area. Increasing the magnetic field or rotation rate increases the maximum energy attainable. Anisotropies and energy densities of the particles are also discussed. It is concluded that the process is consistent with observations of ultra-high-energy cosmic rays.