Diffusive cosmic ray acceleration at shock waves of arbitrary speed Institut für Theoretische Physik, Lehrstuhl IV: Weltraum- und Astrophysik, Ruhr-Universität Bochum, D-4478 Bochum, Germany E-mail: rsch@tp4.rub.de The analytical theory of diffusive acceleration of cosmic rays at parallel stationary shock waves of arbitrary speed with magnetostatic turbulence is developed from first principles. The theory is based on the diffusion approximation to the gyrotropic cosmic ray particle phase space distribution functions in the respective rest frames of the up- and downstream medium. We derive the correct cosmic ray jump conditions for the cosmic ray current and density, and match the up- and downstream distribution functions at the position of the shock. It is essential to account for the different particle momentum coordinates in the up- and downstream media. Analytical expressions for the momentum spectra of shock-accelerated cosmic rays are calculated. These are valid for arbitrary shock speeds including relativistic shocks. The correctly taken limit for nonrelativistic shock speeds leads at relativistic cosmic-ray momenta to the power-law momentum spectrum F (y y q(r and at nonrelativistic cosmic-ray momenta to the power-law momentum spectrum F (y y (+q(r, where the power-law spectral index q(r is a factor greater than the standard spectral index from nonrelativistic shock acceleration theory. Moreover, for nonrelativistic shock speeds we calculate for the first time the injection velocity threshold β c 3β u, settling the long-standing injection problem for nonrelativistic shock acceleration. PoS(ICRC785 35th International Cosmic Ray Conference ICRC7 July, 7 Bexco, Busan, Korea Speaker. This work was partially supported by the Deutsche Forschungsgemeinschaft through grants Schl /9- and Schl /34-. c Copyright owned by the author(s under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives 4. International License (CC BY-NC-ND 4.. http://pos.sissa.it/
. Introduction Apart from magnetic reconnection in transient electric fields cosmic-ray particles in magnetized cosmic plasmas are accelerated by first- and second-order Fermi processes interacting with low-frequency electromagnetic fluctuations such as Alfven and Whistler plasma waves. First-order Fermi acceleration at shock waves is a prime candidate for particle acceleration in astrophysics process because of the large kinetic energy reservoir of supersonic outflows. Pulsar wind nebulae, active galactic nuclei and gamma-ray bursts exhibit highly collimated winds or jets with initial relativistic bulk Lorentz factors Γ = ( (V /c / up to 3. Therefore there is high interest to understand the acceleration of cosmic rays at magnetized shocks of arbitrary speed.. First principles To keep the mathematical complexities at a minimum, we consider the case of an infinitely extended planar shock with a step-like shock profile U(z < = β u c, U(z > = β d c with positive < β d < β u. In order to automatically remove the strong particle momentum anisotropy due to the relativistically moving medium, we take the cosmic-ray phase space coordinates in the mixed comoving coordinate system: time t and space coordinates z in the laboratory (=shock rest system and particle s momentum coordinates p and µ = p /p in the rest frame of the streaming plasma. This choice then allows us to apply the diffusion approximation to cosmic-ray transport in the respective up- and downstream media. Let (p, µ and f represent either (p, µ and f u or (p, µ and f d, respectively. It is also useful to keep the parallel momentum coordinate p = pµ. It is essential that, because of the chosen mixed comoving coordinate system, the up- (p, µ and downstream (p, µ cosmic-ray particle momenta in general are different. In both reference systems the Larmor-phase averaged steady-state Fokker-Planck transport equation (without momentum losses for the anisotropic but gyrotropic cosmic-ray phase space density f (z, p, µ in a medium with magnetostatic turbulence, propagating with the stationary bulk speed U = U(z e z = β s c e z with Γ = [ β s ] / aligned along the magnetic field direction, reads (Lindquist 966; Kirk, Schneider and Schlickeiser 988 [ f Γ[U + vµ] z du p ] f dz Γ = Q(z, p, µ + [ D µµ (µ f ], (. v p µ µ PoS(ICRC785 where Q(z, p, µ = Q (p, µ δ(z is the assumed particle injection rate solely at the position of the shock. With d(γu/dz = Γ 3 (du/dz the transport equation (. can be written as Γ[U + vµ] f z d(uγ dz (U + vµ p f = Q (p, µ δ(z + [ D µµ (µ f ] v p µ µ (. With dimensionless momentum coordinates y = p/(mc, y = p /(mc Eq. (. after a few manipulations reads S(z,y,y z + Γβ s c z [(β s + y y + y f ] = Q (y, µ δ(z + [ D µµ (µ f ] µ µ (.3
with the cosmic-ray current S(z,y,y = cy f Γ + y cβ sγ(β s + y + y f y (.4 Eq. (.3 represents the generalization of the transport equation of Gleeson and Axford (967 accounting for shocks of arbitrary speed and anisotropic but gyrotropic cosmic-ray particle distribution functions. For nonrelativistic shocks Gleeson and Axford (967 showed that, apart from sources, the cosmic-ray "jump" conditions at a shock front are simply that the isotropic phase space density and the isotropic current F(y,z = dµ f (y, µ,z /, S(y,z = dµs(y, µ,z / are continuous.. Diffusion approximation in the up- and down-stream medium For the step-like shock velocity profile the Fokker-Planck transport equation (. in the up- (z < and downstream (z > medium is given by z [Γ(U + vµ f ] = [ D µµ (µ f ], (.5 µ µ We write f (z,y, µ = F(z,y+ g(z,y, µ as the sum of the isotropic part of the cosmic-ray phase space density F(z,y and the anisotropy g(z,y, µ. For small anisotropies g(z,y, µ F(z, p the diffusion approximation (Jokipii 966 provides the diffusion-convection transport equation for the isotropic phase space density F(z, p and the anisotropy g(z,y, µ F Γ[UF Γκ z z ] =, g(z,y, µ = vγ,y 4 A(µF(z z, A(µ = dµ ( µ( µ D µµ (µ µ dx ( x D µµ (x, κ = v dµ µa(µ = v dµ ( µ 8 8 D µµ (µ Demanding as spatial boundary conditions F u (z =,y =, F d (z =,y = F (y, finite F u (z =,y and F d (z =,y at the shock location z =, and spatially-independent diffusion coefficients κ,, the approximated up- and downstram anisotropic (but gyrotropic phase space distribution functions are (.6 PoS(ICRC785 f u (z,y, µ f u (,y, µ e βucz Γuκ, f u (,y, µ = F (y B(µ, f d (z,y, µ = f u (,y, µ F (y, B(µ = + β uβ c 4κ A(µ = β u β A(µ dµ µ A(µ (.7 Mostly important is the isotropy of the downstream distribution function in its rest frame!. Momentum spectrum of accelerated particles at the shock We determine the momentum spectrum of accelerated cosmic rays at the shock by integrating Eq. (.3 from z = η to z = η and considering the limit η resulting in S d (,y,y, S u (,y,y, Q (y, µ 3
+cγ d β d lim [(β d + y η y + y,][ f d (η,y,y, f d (,y,y, ], +cγ u β u lim [(β u + y η y + y,][ f u (,y,y, f u ( η,y,y, ], = lim [ dz η η D µµ (µ f u(z,y, µ η + dz µ µ D µµ (µ f d(z,y, µ ] (.8 µ µ With the solutions (.7 from the diffusion approximation the majority of terms in this lengthly equation vanishes, leaving only with S d (,y,y, = cy, f d (,y,y, S d (,y,y, S u (,y,y, = Q (y, µ (.9 Γ d + y S u (,y,y, = cy, f u (,y,y, Γ u + y Together with the continuity condition (Kirk and Schneider 987 cβ d Γ d (β d + y + y, f d(,y,y, y,, cβ u Γ u (β u + y + y, f u(,y,y, y, (. f u (z =,y, µ = f d (z =,y, µ (. Eqs. (.9 - (. determine the momentum spectrum of accelerated particles at the shock. However, these equations still contain (y,,y and (y,,y which are related to each other using relativistic kinematics. If b = β u β d = β u(r β u β d r βu, Γ r = ( b / = Γ u Γ d ( β u β d (. denote the relative velocity of the upstream medium with respect to the downstream medium and the associated relative Lorentz factor, one finds generally with r = β u /β d and β = y / + y y = Γ r (bµ y + + y, µ = Γ r (y µ + b + y, Γ r (bµ y + + y PoS(ICRC785 y = ( + bβ µ Γ r (bµ y + µ µ = [ + y (β + bµ + y + b( µ y y + y ], µ Γ r y ( + bµ β (.3 [Γ r (bµ y + + y ] 3/ 4
3. Transport equation at the shock In terms of spherical momentum coordinates the up- and downstream cosmic-ray currents (. at the shock are S u (,y, µ = [ + β µ + µ (µ + β u ]F (y B(µ F (y ( µ [(µ + β u B(µ ] cβ u Γ u β u β µ β y [y 3 y F (y µ (µ + β u B(µ ] (3. β and, as f d (,y, µ = F (y = F (y I(b,y is independent of µ, S d (,y, µ = µ [µ + β + (µ + β d ]F (y I(b,y cβ d Γ d β d β y [y 3 y µ (µ + β d F (y I(b,y ], β (3. where we use the continuity condition (. and the last equation (.3 to relate F (y = dµ f d (,y, µ = = I(b,y F (y, I(b,y = Γ ry dµ f u (,y, µ = F (y β dµ µ µ B(µ B(µ ( + bµ dµ (3.3 [Γ r (bµ y + + y ] 3/ With Eqs. (.3 we express the downstream current (3. solely in terms of upstream momenta. The cosmic-ray jump condition (.9 then yields ( β d Γ d F (y I(b,y (b + β µ β d ( + bβ µ + b + β µ β W(b, µ,y + [β d( + bβ µ + b + β µ ] β W(b, µ,y β d Γ ( d ( + bβ µ y 3 W (β 3 + bµ + y + b( µ (b, µ,y y y + y µ ( ( + y y (b + β µ W / (b, µ,y [b + β µ + β d ( + bβ µ ]F (y I(b,y PoS(ICRC785 [ + β µ β u + µ (µ + β u β ]β u Γ u F (y B(µ + β u Γ u F (y ( µ µ [(µ + β u β B(µ ] + β uγ u y [y 3 y F (y µ (µ + β u B(µ ] = Q (y, µ β c with W(b, µ,y = ( + (bµ /β + (b ( µ /y. This equation holds for cosmic-ray particles of arbitrary momentum, shock waves of arbitrary speed and and general injection functions Q (y, µ. As we are particularly interested in the isotropic momentum spectrum F (y of accelerated particles at the shock. we average Eq. (3.4 over µ leading after lengthly algebra (Schlickeiser and Oppotsch 7 to Ω(b,y F (y + d y (y 3 dy F (y T (b,y = dµ Q (y, µ, (3.5 in terms of the convection rate Ω(b,y and the acceleration rate T (b,y. (3.4 5
3. Solution of the transport equation at the shock For the isotropic monomomentum injection rate Q (y, µ = Q δ(y y the solution of Eq. (3.5 is given by F (y y = 3Q y y y 3 T (y,b exp[ For the (in general momentum-dependent power-law spectral index we obtain y dy Ω(b,y yt (b,y ] (3.6 q(r,y = d(lnf /d(y/y d ln(y /y = 3 + y d lnt (b,y dy + Ω(b,y T (b,y, (3.7 depending on the ratio of the convection rate Ω(y,b and the acceleration rate T (y,b. 3. Extreme nonrelativistic shock speeds In the limit of extremely nonrelativistic shock speeds with b = and Γ r = we find W(b =, µ,y = and I(b =,y =, so that the transport equation (3.5 for relativistic cosmic rays reduces to β d F (y + β u β d 3y d(y 3 F (y dy = c dµ Q (y, µ, (3.8 agreeing exactly with the transport equation used in nonrelativistic acceleration theory. However, this case is somewhat unphysical as b β u β d = implies β u = β d, so that no nonrelativistic shock occurs. 3.3 Correct limit of nonrelativistic shocks The correct limit of nonrelativistic shocks with small but finite values of b = (r β d yields to second order in β d for the convection and acceleration rates T (b,y cβ ( dr 3 r βd 3 β Ω(b,y cβ d [r r 3 β + β d r(r (7r 3β ], r [ + (r (7r β d 3β ][ + (3r + (r β d 5β ], (3.9 for symmetric (D µµ ( µ = D µµ (µ upstream Fokker-Planck coefficients. Because the convection rate is always positive, cosmic-ray particles only are accelerated for a positive acceleration rate leading to the condition for the cosmic-ray velocities PoS(ICRC785 β > β c (r = β u uc (r, u c (r = 45r3 + (r (53r + 7 [+ + (r 3 (3r + (7r r 3(r [45r 3 + (r (53r + 7] ], (3. which settles the injection threshold for nonrelativistic shock acceleration. Only cosmic-ray particles with momenta or velocities y β > β c y c are accelerated by nonrelativistic shocks. The injection threshold velocity β c (r, shown in Fig., for all flow compression ratios is always larger than 3β u. 6
Threshold velocity β c β u 8 7 6 5 4 3 exact function approximation 4 6 8 Flow compression ratio r Figure : Injection threshold velocity β c /β u as a function of the flow compression ratio r for a nonrelativistic shock. 3.4 Power-law spectral index With the rates (3.9 the power-law spectral index (3.7 is given by q(r,y y > y c = 3 + y cβ y (β y c + 3r β r 3r β 4 + r(r (7r y c 9(r 3 r + r β, (3. y c which is shown in Fig. for β >.y c and different flow compression ratios. Spectral index q(r, y.y c 8. 7.5 7. 6.5 6. 5.5 r = 3 r = 4 r = 7 PoS(ICRC785 5. 3 3 4 5 Particle momentum log (y Figure : Power-law spectral index for a nonrelativistic shock with β d = 4 as a function of particle momentum y =.y c (r for three different flow compression ratios r = 3,4,7 (full curve. For relativistic cosmic rays with β y c the spectral index (3. approaches q(r,y > = q(r (5r /(r = + (3r/(r, which differs from the classical nonrelativistic shock acceleration theory spectral index 3r/(r by the additional factor. 7
For nonrelativistic cosmic rays with β y y c the spectral index (3. approaches q(r,y c y < = q (r = [3(r /(r ] = + q(r, which differs from the classical nonrelativistic shock acceleration theory spectral index 3r/(r by the additional factor 3. Obviously, our strictly relativistic shock acceleration theory, accounting correctly for the different momentum coordinates in the up- and downstream media, where the diffusion approximation has been applied, implies a much weaker efficiency for the acceleration of cosmic rays at nonrelativistic shocks, when the nonrelativistic limit of small but finite shock speeds in the relativistic theory is considered. 4. Summary and conclusions The analytical theory of diffusive cosmic ray acceleration at parallel stationary shock waves of arbitrary speed has been developed. Starting from the Fokker-Planck particle transport equation in the mixed comoving coordinate system we derived for the first time the correct cosmic-ray jump conditions at the shock relating the upstream and downstream anisotropic cosmic-ray currents. The anisotropic upstream and downstream cosmic-ray currents are calculated from a diffusion approximation of particle transport in the upstream and downstream medium. Pitch-angle averaging the cosmic-ray current jump condition provides a general solution for the isotropic momentum spectrum of accelerated particles at the shock valid for arbitrary shock speeds, arbitrary cosmic-ray momenta and general injection functions at the shock, determined by the ratio of general acceleration and convection rates. The correctly taken limit for nonrelativistic shock speeds leads at relativistic cosmic-ray momenta to the power-law momentum spectrum F (y y q(r and at nonrelativistic cosmic-ray momenta to the power-law momentum spectrum F (y y (+q(r, where the power-law spectral index q(r is a factor greater than the standard spectral index from nonrelativistic shock acceleration theory. Moreover, for nonrelativistic shock speeds we calculate for the first time the injection velocity threshold β c 3β u, settling the long-standing injection problem for nonrelativistic shock acceleration. References PoS(ICRC785 [] Gleeson, L. J., Axford, W. I., 967, ApJ 49, L5 [] Jokipii, J. R., 976, ApJ 46, 48 [3] Kirk, J. G., Schneider, P., 987, ApJ 35, 45 [4] Kirk, J. G., Schneider, P., Schlickeiser, 988, ApJ 38, 69 [5] Lindquist, R. W., 966, Ann. Phys. 37, 487 [6] Schlickeiser, R., Oppotsch, J., 7, ApJ, to be submitted 8