Diffusive Shock Acceleration of HighEnergy Cosmic Rays Joe

Diffusive Shock Acceleration of High-Energy Cosmic Rays Joe Giacalone & Randy Jokipii University of Arizona • The origin of the very-highest-energy cosmic rays remains one of the most-important unsolved problems in astrophysics. • One of the most-important accelerators is the mechanism of diffusive shock acceleration • In this talk, we discuss this mechanism and, in particular, the upper energy limits and the importance of the magnetic field angle “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

The Importance of the Magnetic-Field Angle • Acceleration to high energies: – Parallel Shocks • Very slow • Efficient – Perpendicular Shocks • Much faster • Also efficient (we point out in this talk that there is no injection problem) • New numerical simulations – Hybrid simulations (self consistent) show efficient acceleration of thermal ions by a perpendicular shock “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

“Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

“Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Actual Particle Orbits – Including Scattering Decker, 1988 “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Diffusive Shock Acceleration • Start with the cosmic-ray transport equation advection diffusion drift “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05 energy change

Diffusive Shock Acceleration • Solve the cosmic-ray transport equation for the following geometry “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

• The steady-state solution for given by , for an infinite system, is The downstream distribution is power law with a spectral index that depends only on the shock compression ratio! “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05 Kennel et al, 1986

The observed quite-time cosmic-ray spectrum “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

What about the maximum energy attainable? • The energy is limited either by geometry or by a finite time. • Acceleration takes time. The ideal power law energy spectrum is not created instantly. • The maximum energy over a given time interval strongly depends on the shock-normal angle Parallel shocks slow Perpendicular shocks fast “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Acceleration Time in Diffusive Shock Acceleration • The timescale to accelerate particles to a momentum p (from a much smaller momentum p 0) is given by Where “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

One may conclude: • The acceleration time depends on • • The rate of change of the maximum energy is much larger for perpendicular shocks. • Hence, for any given situation, a perpendicular shock will yield a larger maximum energy than a parallel shock. • Most discussions of the maximum energy have been limited to a parallel shock, for no good reason. “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Acceleration Rate as a Function of Shock-Normal Angle: (assumes the billiard-ball approximation) Jokipii, 1987 “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Maximum Energy Assumes Sedov solution for SNR blast wave Bohm Diffusion Perpendicular Shock (Hard-sphere scattering) “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

What about Injection and the limit of diffusive shock acceleration? • An often-invoked injection criterion is • This assumes, for no good reason, that there is NO motion normal the average magnetic field – In general, particles move normal to the field, and this is important for the injection problem “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

It is useful to examine the limits of theory of Diffusive Shock Acceleration Because the distribution should be nearly isotropic, we require that the diffusive streaming anisotropy be small. where is the streaming flux and “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

This leads to a general expression for the validity of diffusive shock acceleration Diffusive shock acceleration is applicable “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Special Cases of the general limit: “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Special Cases of the general limit: “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Special Cases of the general limit: “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

The case of hard-sphere scattering Classical-scattering theory gives For most astrophysical applications Thus, using classical-scattering theory, we find “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

The case of hard-sphere scattering Classical-scattering theory gives For most astrophysical applications Thus, using classical-scattering theory, we find However, classical-scattering theory is NOT a good approximation for perpendicular transport ! “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Test-particle simulations using synthesized magnetic turbulence (Giacalone and Jokipii, Ap. J, 1999 + one extra point) “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

The case of field-line random walk Thus, for a perpendicular shock, we find → The same as for a parallel shock ! “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

How does the increased diffusion normal to the field effect the acceleration time? • The ratio of the acceleration time to that for Bohm diffusion is given by: “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Test-particle simulations using synthesized magnetic turbulence (Giacalone and Jokipii, Ap. J, 1999 + one extra point) “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

• To test these ideas, we have performed test-particle simulations • Particles are followed in the synthesized shock fields • The magnetic field is composed of a mean, plus a fluctuating component derived from an assumed power spectrum – Maxwell’s equations are satisfied across the shock Assumed geometry in test-particle simulations “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Test-particle simulations • Shock moves through a plasma with a magnetic field derived from a power spectrum (e. g. Decker & Vlahos, 1986) Assumed power spectrum “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

More results from test-particle simulations Effect of Turbulence Amplitude Time dependent case “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

More results from test-particle simulations “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

The Injection Problem • What is the process by which thermal particles become part of the highenergy population? • The process is nonlinear and requires a kinetic description “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

The Fields Must be Fully Three Dimensional • If the fields contain an ignorable coordinate, then the particles are unphysically forced to remain within a single gyroradius from the field line that it started “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

The Hybrid Simulation • Self-consistent plasma simulation (kinetic ions, MHD fluid electrons, often assumed to be massless) • Used to study the structure of shocks, ion distributions, upstream ions, etc. , acting on ion length and time scales • Can also be used to study the acceleration of thermal ions to high energies. “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

1 D Hybrid Simulation of a Parallel Shock Downstream Energy Spectra “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

New 2 D Hybrid Simulations • We have performed large 2 D simulations (500 × 4000 ) to investigate the effect of longwavelength magnetic fluctuations on the acceleration of thermal ions at a shock. • “Seed” magnetic fluctuations are imposed on the system • Particles are tied to field lines, but move normal to the mean field by following meandering lines of force “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Simulated Magnetic Field Bz “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Energetic-Particle Intensity -Zoomed-in View “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Large vs. Small Domain-Size Comparison “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

vs. Shock Comparison “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05

Conclusions 1. Shocks provide a natural explanation for most cosmic rays. 2. The maximum energies generally quoted for shock acceleration greatly underestimate the maximum energy, as they apply only to parallel shocks. 3. Considerations of quasi-perpendicular shocks lead to much higher possible energies. 4. Perpendicular shocks are as efficient at accelerating thermal particles as are parallel shocks. There is not an “injection” problem for shocks occurring in a wide-variety of astrophysical plasmas. “Physics at the End of the Galactic Cosmic-Ray Spectrum” Aspen, CO 4/28/05