arxiv: v3 [astro-ph] 13 Nov 2008

Transcription

1 Mon. Not. R. Astron. Soc., 8) Printed 9 May 8 MN LATEX style file v.) On Fermi acceleration and MHD-instabilities at ultra-relativistic magnetized ock waves arxiv: v3 astro-ph] 3 Nov 8 Guy Pelletier, Martin Lemoine and Alexandre Marcowith 3 Laboratoire d Astrophysique de Grenoble, CNRS, Université Joseph Fourier II, BP 53, F-384 Grenoble, France; Institut d Astrophysique de Paris, UMR 795 CNRS, Université Pierre&Marie Curie, 98 bis boulevard Arago, F-754 Paris, France 3 Laboratoire de Physique Théorique et Astroparticules, CNRS, Université Montpellier-II, place Eugène Bataillon, 3495 Montpellier Cédex, France INTRODUCTION The physics of collisionless ock waves is of paramount importance in modern astrophysics, as it governs a variety of phenomena, most notably the emission of high energy radiation. In particular, the afterglow radiation of gamma-ray bursts is generally attributed to the synchrotron emission of electrons that have been accelerated around the relativistic external ock wave, with ock Lorentz factorγ Paczyński & Rhoads 993, Katz 994, Mészáros & Rees 997, Sari & Piran 997, Vietri 997, Waxman 997; see Piran 5 for a review). In this framework, it can be own that the magnetic field, downstream of the ock wave, must have been amplified by orders of magnitude beyond the ock compression of a typical interstellar magnetic field see for instance Waxman 997; see also Piran 5 and references therein). Recent work by Li & Waxman 6) further ows that even the upstream magnetic field must have been amplified by at least one or two orders of magnitude in order to explain observed X ray afterglows. These results ould be put in perspective with recent studies on the process of Fermi acceleration in the test particle limit around relativistic ock waves. In particular, Niemiec & Ostrowski 6) guy.pelletier@obs.ujf-grenoble.fr lemoine@iap.fr alexandre.marcowith@lpta.inp3.fr ABSTRACT Fermi acceleration can take place at ultra-relativistic ock waves if the upstream or downstream magnetic field has been remodeled so that most of the magnetic power lies on ort spatial scales. The relevant conditions under which Fermi acceleration become efficient in the presence of both a coherent and a ort scale turbulent magnetic field are addressed. Within the MHD approximation, this paper then studies the amplification of a pre-existing magnetic field through the streaming of cosmic rays upstream of a relativistic ock wave. The magnetic field is assumed to be perpendicular in the ock front frame, as generally expected in the limit of large ock Lorentz factor. In the MHD regime, compressive instabilities seeded by the net cosmic-ray charge in the ock precursor as seen in the ock front frame) develop on the ortest spatial scales but saturate at a moderate level δb/b, which is not sufficient for Fermi acceleration. As we argue, it is possible that other instabilities outside the MHD range provide enough amplification to allow successful Fermi acceleration. Key words: ock waves acceleration of particles cosmic rays have own that Fermi acceleration becomes inoperative in a turbulent magnetic field with power spectra of the Kolmogorov or scale invariant type. Through analytical calculations, Lemoine, Pelletier & Revenu 6) have own that relativistic Fermi acceleration is inefficient if the magnetic field power is distributed on scales larger than the typical Larmor radius of the accelerated population, in agreement with the above numerical result. Therefore the interpretation of the gamma-ray burst afterglow emission from the relativistic external ock wave requires the magnetic field to have been amplified on ort spatial scales. It has been suggested that the relativistic two stream Weibel instability could amplify the downstream magnetic field to the level required by afterglow modeling of gamma-ray bursts Gruzinov & Waxman 999; Medvedev & Loeb 999). Subsequent studies have however argued that this instability ould saturate early on Wiersma & Achterberg 4; Lyubarsky & Eichler 6). The question of the long term evolution of the downstream magnetic field on timescales ω p, withω p the plasma frequency) also remains open. Ongoing particle-in-cell simulations ould eventually ed light on this issue Silva et al. 3; Frederiksen et al. 4; Medvedev et al. 5; Kato 7; Chang, Spitkovsky & Arons 8; Keet et al. 8; Spitkovsky 8). Regarding the amplification of the upstream magnetic field, the present situation is reminiscent of results obtained for nonrelativistic supernovae remnant ock waves, where the interstellar

2 G. Pelletier, M. Lemoine and A. Marcowith magnetic field has apparently been amplified by one or two orders of magnitude see Völk, Berezhko & Ksenofontov 5, Parizot et al. 6). In this case, the leading candidate for the instability is the streaming instability, seeded by the cosmic ray precursor in the upstream plasma see Bell 4, 5; Pelletier, Lemoine & Marcowith 6, Marcowith, Lemoine & Pelletier 6; Reville et al. 7; Niemiec et al. 8; Amato & Blasi 8; Zirakavili, Ptuskin & Völk 8; Reville et al. 8). The generalization of this instability to the relativistic regime has been studied on phenomenological grounds by Milosavljević & Nakar 6) who have concluded that it ould be able to account for the degree of amplification inferred in gamma-ray bursts. More recently, Reville, Kirk & Duffy 6) have derived in detail the dispersion relation for a parallel ock wave and the saturation due to thermal effects. In the present work, we propose to explore in detail the generalization of this type of instability to the ultra-relativistic regime. One crucial difference with the previous works on the streaming instability is that we consider the most natural case of superluminal ock waves with a magnetic field perpendicular to the ock normal in the ock front frame). This case is more generic than the parallel configuration studied previously because the transverse component of the magnetic field is boosted by the ock Lorentz factor when going to the ock frame. Another important difference is that we bring to light a new type of instability, of a compressive nature. Finally, in contrast with most particle-in-cell simulations performed to date, our study focusses on magnetized ock waves, for which there exists a coherent upstream magnetic field whose dynamical influence on the ock jump conditions can be neglected however). Nevertheless there exist pioneering PIC simulations in the moderately relativistic regime, which include both a mean field and a significant mass ratio between electrons and ions, see Hededal & Niikawa 5), Dieckmann, Shukla, & Drury 8). We adopt a simplified description in which the cosmic-ray distribution is modeled as a step function out to some distancel cr and we neglect the cosmic-ray response to the disturbance. This latter assumption is justified by the fact that the instability is maximal on the ortest spatial scales, orders of magnitude below the typical Larmor radius of accelerated particles. The paper is organized as follows. In Section, we introduce the main scales of the problem, most notably the diffusion scale of the cosmic rays; we then calculate the level of amplification that is necessary to make Fermi acceleration operative. Section 3 is devoted to the investigation of the instabilities under the condition that some cosmic-rays have undergone a first Fermi cycle. Section 4 summarizes out results and provides some outlook. Details of the calculations are provided in Appendix A. GENERAL CONSIDERATIONS We carry out most of the discussion in the ock front rest frame, hence unless otherwise noted, all quantities are evaluated in this frame. We use the subscripts u or d to tag quantities measured in the upstream or in the downstream rest frame respectively.. Upstream diffusion length In the upstream rest frame, cosmic rays can never stream too far ahead of a relativistic ock wave since this latter propagates towards upstream with velocity v =β c c the ock Lorentz factor β ) / ]. Cosmic rays scatter on magnetic turbulence upstream before they are caught back by the ock wave when their pitch angleθ u / Gallant & Achterberg 999, Achterberg et al. ). Consequently, they can travel a distance l cr u, which may take the following values depending on the ratio of the Larmor radius r L u to the coherence length of the upstream magnetic fieldλ c u Milosavljević & Nakar 6): for small scale turbulence l cr u r L u ) rl u Γ Γ λ c u, ) λ c u for large scale turbulence l cr u r L u rl u λ c u ). ) Both regimes, ort or large scale turbulence can be expected at some point, insofar as the excitation of the upstream magnetic field on ort spatial scales is due to the streaming of the nonthermal particle population in the ock precursor. Indeed, cosmic rays of the first generation are to interact with a turbulent magnetic field ordered on large scales. However, provided the instability that they trigger grows fast enough, cosmic rays of the next generation will propagate in ort scale turbulence. In reality, the situation is likely to be more complex as the process of particle propagation upstream and magnetic field generation are closely intertwined. The fact that the non-thermal population contains particles of different energies, which can stream at different distances from the ock front, ould also play a significant role. In this respect, Keet et al. 8) have observed that the upstream magnetic field is affected to greater distances as time goes on. This strongly suggests that higher energy cosmic rays are produced as time goes on, and that, by travelling farther in the upstream medium, they excite the turbulence at larger distances from the ock front. In the discussion that follows, we estimate the growth of unstable modes in both the limits of small or large scale turbulence in order to remain as general as possible. Certainly the limit of large scale turbulence is more restrictive with respect to the growth of the instability, since the distance travelled upstream is significantly reduced with respect to that in small scale turbulence. Whichever limit prevails depends on the ratio of the ort scale turbulent to the large scale coherent) magnetic field strength, as discussed in Section.. We also discuss the effect of higher energy cosmic rays on the growth rate. It is important to emphasize that the distance that controls the growth of the instability is that between the ock front and the position of the particle, which is smaller thanl cr u by a factor β ). Γ) In the following, we will need the expression forl cr also notedl cr ), i.e. the length scale of the cosmic-ray distribution as measured in the ock front rest frame. It can be calculated by transforming the upstream residence time t r u l cr u /c in the ock front frame t r = t r u /, and then by rewriting in the expression obtained the upstream Larmor radius and coherence length in terms of their ock frame equivalent. In the perpendicular or superluminal and ultra-relativistic) configuration of interest, r L Γ r L u,λ c Γ λ c u. This boost of the coherence length is valid for wavenumber modes that are parallel to the ock normal; perpendicular wavenumber modes remain unchanged. In the following, we thus use the following expressions forl cr : for small scale turbulence l cr l cr r L λ c, 3)

3 MHD-instabilities at ultra-relativistic magnetized ock waves 3 for large scale turbulence l cr l cr r L. 4) Below, we find a process of generation of intense magnetic disturbances at ort spatial scales. The expected instabilities in all relevant cases generate disturbances with a coherence lengthλ c u between the minimal length l MHD = c/ω p,i required for the validity of MHD description and the diffusion lengthl cr u of the cosmic rays, but with a preference for ort scale, i.e. l MHD λ c u l cr u. The minimal MHD length is: l MHD β A r u, 5) whereβ A B u / 4πρ u u c denotes throughout this paper the Alfvén velocity measured upstream) in units of c, and r u = m p c /eb u is the Larmor radius of thermal protons in the upstream comoving frame. The above MHD scale is measured in the upstream plasma rest frame.. Requirements for successful Fermi acceleration As discussed in Lemoine, Pelletier & Revenu 6), a particle can execute a large number of Fermi cycles through the ock only if the turbulence coherence scale is much smaller than the Larmor radius assuming no coherent magnetic field), in detail r L u λ c u upstream, or r L d λ c d as measured downstream. The Fermi process becomes inoperative in the limit of large coherence length when compared to the typical Larmor radius) since the cycle time becomes smaller than a coherence time hence the field is effectively regular and superluminal on a typical Fermi cycle. Conversely, the Fermi process will be operative provided either of the above inequalities is fulfilled, since the memory of the initial conditions at ock crossing is then erased by pitch angle diffusion in the ort scale turbulence. In the generic case of a ultra-relativistic superluminal ock, the above two conditions are actually related by the ock jump conditions, since the Larmor radius increases by Γ rel when going from downstream to upstream, while the coherence length increases by Γ rel. The quantityγ rel corresponds to the Lorentz factor of the downstream fluid as measured upstream, and in the case of a ultra-relativistic strong ock,γ rel / Blandford & McKee 976). In the following, we describe the upstream and downstream fluids as comprising a large scale component B and a ort scale turbulent generated by some unspecified instability. In this case, Fermi acceleration will be efficient provided the ort scale field has a sufficient amplitude as compared to B. We discuss this condition in the following... Upstream motion In this particular section, all quantities are evaluated in the upstream rest frame and we avoid using the notation u out of clarity. We derive here the requirements on the amplitude of the ort scale component of the magnetic field that would allow successful Fermi acceleration. Without loss of generality, we assume that the coherent component B lies in the plane x y, the x direction corresponding to the ock normal oriented toward upstream infinity). We denote by b the irregular component whose amplitude is expressed in units of B ). The layout is pictured in Fig.. The total field is thus written as B= B e+b), and the particle downstream z φ θ β ock front y Β x upstream Figure. Schematic plot of the geometrical configuration in the upstream plasma rest frame. The coherent component B lies in the x y plane, with x pointing towards the ock normal, and θ,φ) the angles of the particle velocity in this frame. trajectory is governed by the following equation: dβ dt =ω L, β e+ b), 6) withω L, denoting the Larmor frequency expressed relatively to B. We neglect the influence of any ort scale time varying electric fieldδe in the above equation of motion sinceδe/b b) ω/k. In order to quantify the return timescale and the condition for successful Fermi acceleration as a function of the amplitude of the random component of the magnetic field, it is useful to express this equation in terms of the anglesθandφ, which are defined as: β x = cosθ,β y = sinθcosφ andβ z = sinθ sinφ. These equations of motion read: ) θ = ω L, ey + b y sinφ bz cosφ ], 7) φ = ω L, e x b x + ) ] cosφ e y + b y tanθ + b sinφ z. 8) tanθ Rescaling the Larmor frequency with respect to e y B, i.e. defining ω L,y e y ω L,, and noting thatθ when the particle propagates upstream, the equations of motion forθandφcan be rewritten as follows: θ ω L,y sinφ)+ω L,y A sinφ σ), φ ω L,y cosφ θ +ω L,y A cosθ σ) θ. 9) The factors A and σ describe respectively the amplitude and the phase of the small scale magnetic field in the ock front plane, i.e. b y = A cosσ and b z = A sinσ assuming isotropic turbulence). The amplitude A=δB/B y is measured with respect to B y. One can extract from the above system the unperturbed trajec-

4 4 G. Pelletier, M. Lemoine and A. Marcowith tory, i.e. assuming A = : sinθ i cosφ=cosφ i sinθ, ) where values indexed with i are calculated at some initial time. Since entry into upstream corresponds to cosθ i > β hence θ i /, and exit from upstream corresponds toθ /, the above equation suggests thatφis driven to an angle close to±π/. The sign is given by the sign of e y : for e y >, it is+π/. This region lies opposite to that which allows return to the ock when the particle travels downstream see Lemoine, Pelletier & Revenu 6), hence Fermi cycles cannot be completed. It is important to note that this unperturbed trajectory is executed on a timescale t unpert r L,y / c) where r L,y = c/ω L,y ). The noise term comes from the random phaseσassuming isotropic turbulence in the ock front plane, an assumption that can be relaxed). Therefore the trajectory deviates from the unperturbed trajectory to linear order by φ, θ), whose variance increases linearly with time: φ θ 3 Aω L,y θ τ c t, 3 A ω L,y τ c t. ) Note that the above equation does not contain all the terms driving the variance of φ and θ: it neglects term of the form dt dt φt ) φt ) and similarly in θ. However, these terms are smaller byoa ) as compared to those above, hence they can be neglected in a first approximation we will find A further below). In the absence of the regular driving term in the equations of motion i.e. when B y = ), return occurs over a timescale t diff Γ A ω L,y τ c). Noise dominates over the unperturbed trajectory if θ θ, whereθindicates the unperturbed trajectory. This inequality must be satisfied within a return) timescale t unpert otherwise the particle will have exited from upstream along the weakly perturbed unperturbed trajectory) before the noise could overcome the unperturbed driving force. This implies: A ω L,y τ c. ) If one defines the Larmor radius r L measured with respect to the total magnetic field, r L = c ], A + ) / ω L, hence rl c/aω L,y ) when A, the former inequality amounts to: A r L λ c. 3) If this inequality is satisfied, then one can check that φ / O) on a return timescale t diff, which implies that the return directions are isotropized in the ock front plane. Fermi acceleration ould then be efficient provided the return probability to the ock as calculated downstream is also isotropic inφ, see below). An interesting implication of the above is to restrict Fermi acceleration to a rather limited range of Larmor radii: λ c r L A λ c. 4) Very high values of the amplification factor are thus required to produce powerlaw spectra over a large dynamic range... Downstream motion All quantities in this subsection are evaluated in the downstream rest frame. In this frame, the coherent magnetic field is nearly fully aligned with the y axis, and we assume indeed that B = B y. It is useful to study the evolution of the velocity components perpendicular β ) and parallel β ) to the direction of B :β β y, β =β x +β z. The equations of motion for the velocityβcan then be written: β = ω L, β β β x b z β z b x ) β = ω L, β z b x β x b z ). 5) One can check thatβ β +β β = as it ould sinceβis conserved. The unperturbed trajectory for the couple of variables β,β ) is trivial, i.e.β, = const. =β,i andβ, = const. =β,i. Return may occur or not, depending on these initial velocity components; if it occurs, it does so on a timescale t unpert r L, /c Lemoine, Pelletier & Revenu 6). Since b x is not compressed through ock crossing, b x /b z / for isotropic upstream turbulence, hence this component can be neglected. Integrating the above equations to calculate the squared perpendicular displacement, one finds: t t x = dt dt β t )β t ) ω L, ω L, A β t,i β,i t t t dt dt dt 3 dt 4 β x, t 3 )β x, t 4 ) b z t 3 )b z t 4 ) β t,i β,i 3 ω L, τ c β,i β,i dt t t dt t dt 3 A 3 τ cβ x, t 3) Pt). 6) In the above equation,β x, t) represents the unperturbed trajectory: β x, = β x,i cosω L, t) β z,i sinω L, t). The function Pt) contains powers of t up to t 3, as well as sine and cosine functions of ω L, t; its dimension is that of c ω 3 L. Noise will then dominate over the unperturbed trajectory x, = β,i ct over a unperturbed return timescale if: or, equivalently: A ω L, τ c, 7) A r L λ c. 8) As before, r L denotes the Larmor radius measured with respect to the total magnetic field. This condition on A is very similar to that obtained upstream, since r L u Γ r L d andλ c u λ c d, while A u A d. According to the above analysis, and following Lemoine, Pelletier & Revenu 6), phase mixing ould be sufficiently large to erase all dependence on the phase of the velocity vectors, hence Fermi acceleration ould be fully operative, if either one of the above conditions on A is satisfied. It is interesting to discuss the above results in light of recent Monte Carlo numerical simulations of relativistic Fermi acceleration in the presence of ort scale turbulence Niemiec, Ostrowski & Pohl 6). These authors have investigated the efficiency of the Fermi process for various turbulence configurations, including a coherent magnetic field, a long wave turbulence and a ort scale component. They confirm that Fermi acceleration is not operative ifδb/b=, but find that a spectrum can develop over about two orders of magnitude whenδb/b=8 see their Fig. ). The high

5 MHD-instabilities at ultra-relativistic magnetized ock waves 5 energy break can be directly interpreted as the energy at which the above inequality Eq. 8) is no longer fulfilled. Beyond that point, the spectrum steepends and the Fermi process becomes inoperative. Interestingly, these authors also ow that if the pitch angle scattering amplitude in the ort scale turbulence is independent of energy, the break disappears and the spectrum continues without bound. This is also expected, insofar as the break that is determined by Eq. 8) stems from the fact that the scattering time in the ort scale turbulence scales as the square of the Larmor time whereas the unperturbed trajectory in the background field depends linearly on the Larmor time. Assuming a constant pitch angle scattering amplitude in the notations of Niemiec, Pohl & Ostrowski 6) implies that the two scattering timescales evolve in the same way, hence if ort scale turbulence suffices to isotropize directions at a given energy, it will do so at all energies. In order to produce powerlaw spectra over a large dynamic range, it is thus necessary to reach an amplification factor δb/b as large as possible, since this ratio determines the dynamic range of particle energies. Interestingly, the afterglow modeling of gammaray bursts suggest that indeed, the usptream and downstream magnetic fields have been amplified by a large factor. In the range of energies in which Fermi acceleration is operative, it is likely that the spectral index would be equal to the so-called canonical value s =.3, although the generality of this prediction for different apes of the 3d turbulence spectrum remains to be studied. The results of Niemiec, Ostrowski & Pohl 6) cannot be used to infer this spectral index, since they have restricted their analysis to values ofδb/b<, and the spectral indices they report have been measured beyond the break energy. Techniques developed in Lemoine & Pelletier 3), Lemoine & Revenu 6) are particularly suited to study this problem and calculations are underway. 3 INSTABILITIES AT PERPENDICULAR SHOCK WAVES The magnetic field in the ock front B is related to the magnetic field in the upstream frame B u and the associated electric field E u through the Lorentz transform: B = B u B = B u β E u ). 9) To zeroth order, the electric field in the upstream plasma frame vanies, hence the magnetic field in the ock front frame is mostly perpendicular unless B u is aligned along the ock normal to within an angle O/ ) subluminal ock). It then suffices to consider the fully perpendicular situation with B =. The cosmic rays stream ahead of the ock wave, carrying a net charge densityρ cr which will induce a counteracting charge densityρ pl in the background plasma. Note that the cosmic rays do not induce an electrical current at zeroth order in the ock front frame, only a net charge density. Since we consider the generation of ort scale turbulence, we neglect the cosmic ray response to this ort scale turbulence other than its effect on the cosmic ray distribution scalel cr. For simplicity, we approximate the cosmicray charge profile with a top-hat distribution. As usual, we search for a stationary regime which serves as a basis for perturbating the equations in the time-dependent regime. The details of the calculations are provided in Appendix A for both stationary and time dependent quantities. In particular, the stationary regime exhibits the following set-up: B= B y e y, u=u x e x + u z e z. ) Recalling that the ock normal is directed along x, u x β characterizes the velocity of the upstream inflowing through the front. The velocity component u z along the front is small compared to u x, but its ear x u z cannot be neglected see Appendix A). 3. Linear analysis of the reduced system The complete system that governs the linear evolution is of seventh order see Appendix A) and thus rather involved. Nevertheless, it remains possible to derive the main results since the Lorentz transform from the upstream comoving frame to the ock front frame dominates the MHD propagation effects when the wave scale is not too small, as will be made more precise in a next subsection. In other words, the instability growth rates can be derived to lowest order by assuming vaniing Alfvén and sound velocities. We have verified that the analysis of the system includingβ A andβ s but neglecting terms in u z ) does not modify the growth rates obtained further below. One ould point out that electromagnetic waves with dispersion relationω u k u ) are Lorentz transformed into waves with dispersion relation: ) ky,z ω β ck x +O. ) In particular, for Alfvén waves, the next term on the r.h.s. is β A k y c/. Note that the Alfvén velocity is always expressed in the upstream plasma rest frame, i.e.β A = B y u / 4πρ u c. Therefore it suffices in what follows to consider the limit k z. The limit k y can also be considered but it appears more restrictive, because it limits the analysis to the evolution of upstream magnetosonic modes. We will thus consider both cases k y = and k y finite. Of course, one can use Fourier analysis without mode coupling as long as the spatial dependence of the system coefficients can be neglected, in agreement with our previous approximations, and in particular, that the cosmic ray charge density is modeled as a step function. In the ock front rest frame, the physical picture of the instability is then as follows. Impingent electromagnetic waves propagate in vacuum beyond the length scalel cr which characterizes the spatial distribution of cosmic rays upstream of the ock front. For < x<l cr, the presence of the cosmic-ray charge density induces a ort scale instability. Hence we consider electromagnetic waves and seek a solution in Fourier space withωset by its vacuum dispersion relation, but with a complex k x characterizing the amplification in the charge layer < x<l cr. In this perpendicular configuration, the return current in the upstream frame) is not responsible for a supplementary tension effect, but for a supplementary compression effect. A b y perturbation leads to a supplementary vertical i.e. along z) compression that can pu like the kinetic compression. If we consider a spatial modulation along the mean field k y ), a vertical perturbation b z generates a supplementary compression in the direction of the mean field that can pu in phase with the kinetic compression as well. For large even with k z, the system reduces to fourth order and can be expressed as the coupling between the propagation of the vertical perturbed motionδu z and the propagation of kinetic compressionδw/w, where w ρ u c is the relativistic proper enthalpy density of the cold upstream plasma. The system of the two coupled equations can be expressed using the system of Eqs. A8),A9),A) in the limitβ A,β s. In particular in this limit, one notes that b x. Then one obtains a single equation for δu z see Ap-

6 6 G. Pelletier, M. Lemoine and A. Marcowith pendix A.) where derivatives z cancel out and where the derivative of u z is inserted, its second derivative being neglected we also assumeβ = in the coefficients): ˆD 4 δu z κ + u z κ x ) ˆD δu z + κ x ˆDδu z +Γ κ y δu z=. ) The differential operator ˆD is defined by ˆD= c t β x ). For a detailed understanding of the instabilities one has to notice that the charge is not completely screened by the plasma as a supplementary electric field is generated along the flow i.e. oriented towards x): E x = u z B y /, along with a eared vertical motion u z such that x u z =κ see Appendix A). The instabilities stem from this eared motion. The parameter κ, which carries the dimension of a wavenumber, is defined in Eq. A3). To our present order of approximation, it can be expressed as: κ ρ crb y w. 3) Note that bothρ cr and B y are here evaluated in the ock front frame. In the front frame, the parameterκappears divided by, hence we introduce k κ/. The reduced form of the equation see Appendix A.) clearly ows the ordering by setting ˆD=k x β D,δu z =β ũ z ): k D 4 ũ z k x + u z k x k x ) D ũ z + k x k 3 x Dũ z + k k 4 yũz= 4) x The inhomogeneity is contained inκ that vanies at large distance from the ock. At infinity ˆDδu z =, which implies that for a mode having a specified k x, the pulsationω=β k x c, as expected. In the precursor whereκ, we follow this mode characterized by its frequency and look at the modification of its spatial behavior. As explained above, we solve this equation by setting c t iβ k x, x ik x ε), so that ˆD ik x β ε D iε), whereεis a complex number, the imaginary part of which characterizes the growth rate. It is obtained by solving the algebraic equation: ε 4 + k k x ε k ε k ky + i u z k ε)ε = 5) k x k 4 x k x Actually, becauseεis small, the equation can be simplified to: ε 4 k ε k ky k x k 4 x =. 6) The contribution in u z introduces a small correction to the real part of the roots most u z contributions can be canceled out by a frame transformation, but that correction cannot be canceled out because it corresponds to the electric field component E x ). 3.. Analysis of the case k y = Let us first analyze the case k y =, or less restrictively for ky ε k x. This is a pure magneto-sonic compression with b along the mean field. One gets the following three roots: ) /3 k ε, ±i 3, k x 7) One obtains an instability spatial growth rate x = Im k x ε) that increases with large k x. x = 3 k/3 k /3 x k x k ). 8) Thus the characteristic wave number associated with the charge in this problem of magneto-sonic wave propagation is k κ/. 3.. Analysis of the case k y With k y, one recovers the previous growth rate if k y k /3 k /3 x, while for k y k /3 k /3 x one derives : ) ky k / ε,, i, i), 9) k x or, equivalently: x ) / ) k y k ky k /3 k /3 x. 3) The spatial growth rates are thus x k k x ) /3 and k k y ) / respectively. Remarkably they are independent of k z that can be chosen arbitrarily, provided it remains small compared to k x for the consistency of the derivation Comparison to MHD and cosmic-ray length scales One ould compare the scale defined by k to the smallest scale l MHD for our MHD description. Note that this length scale is defined in the comoving upstream frame. Since k x u k x / while k y u k y, it suffices to require k l MHD < to ensure that there exist modes in the MHD range with k>k. One obtains: k l MHD = n creb y Γ wβ Ar u 3) which can be rewritten in terms of the fractionξ cr of downstream) ock internal energy converted into downstream) cosmic ray energy: ξ cr e cr d 3) 4Γ w. To this effect, we assume that the accelerated population can be described as a power-law of index s and minimum momentum p min, so that: n cr s e cr s p min c. 33) This equation assumes s>; if s=, then s ould be replaced by logp max /p min ). Since e cr d e cr as d = 9/8 for a strong ultra-relativistic ock), one obtains: s k l MHD 4 s ξ r u crβ A Γ r. 34) min u Here r min u denotes the minimum Larmor radius of the accelerated population, as measured upstream. The typical energy of the first generation of cosmic rays isγ m pc, so that r min u Γ r u. Therefore s k l MHD 4 s ξ crβ A. 35) Thus k defines an MHD scale. In order for the instability to be efficient, one must also require x l cr, which is a non-trivial requirement in view of the restricted value ofl cr. One finds: s x x l cr 4 ξ cr g, 36) s k with g for large scale turbulence, and g r u /λ c u for ort scale turbulence. In the latter case, one may assumeλ c u l MHD, as suggested by the fact that the growth rate increases with k, so that g /β A. Then, in both cases, the total growth is presumably much larger than unity ifξ cr is not too small and k x, k y sufficiently

7 MHD-instabilities at ultra-relativistic magnetized ock waves 7 large as compared to k. Hence even the first generation of cosmic rays is able to destabilize the upstream magnetic field on the ortest scales. The effect of the high energy part of the accelerated population will be addressed ortly. As mentioned, the study of the more extended system including the contribution of the terms involving the Alfvénic and sonic contributions but neglecting terms of order u z in the system) confirms the results above for the growth rates, providedβ A and β s Magneto-sonic saturation One can estimate the amplitude of the various components in terms of b y. It has already been mentioned that b x b y, and b y β δu x 37) b z κ yδu x. 38) Thus we find that b z may achieve a large amplitude, because b z κβ y b y. 39) In the frame of linear theory, the saturation level of these magnetosonic instabilities is simply estimated by the fact that the compression has a limited amplitude δw/w <. Defining the power spectrum of b y and b z per log interval of wavenumberp by k 3 /π ) b y in terms of the Fourier component b y and similarly forp bz in terms of b z, one derives from Eq. A33): and P / b y P / b z k ε k x+ k 4) = k y k P / b y. 4) Note that the power spectrum is normalized to the ratio of amplitude of the turbulent small scale) component to the coherent B. Consider the first branch of the instability, given by Eq. 8) which applies in the limit k x k, k y k /3. There one can ow thatp / b y saturates at a value of order k x /k ) /3, butp / b z saturates at a value of order unity when k y k /3. For the k /3 x k /3 x second branch, given in Eq. 3) and which applies in the limit k y k /3 k /3 x, one finds thatp / b y saturates at a value of order k /k y, but againp / b z saturates at a value unity. This moderate saturation level can be directly traced back to the saturation of these compressive modes at the linear level δw/w. As discussed in Sec.., such a levelδb/b ould not allow a powerlaw spectrum to develop through Fermi acceleration. Nevertheless, our present analytical description is by definition limited to the linear regime. It would certainly be interesting to pursue the calculations using numerical simulations as one could expect that the holes produced in the plasma would widen with the growing magnetic pressure inside and that the density excesses become spikes. Thus one might reasonably expect the generation of many local small magneto-sonic ocks that are absorbed by the main ultra-relativistic ocks, with the reconversion of some amount of cosmic-ray energy into thermal energy in the precursor. The previous discussion has focused on the role played by the first generation of cosmic rays, since higher energy cosmic rays can only be present if acceleration to smaller energies has been completed. However it is easy to verify that, if highest energy cosmic rays are present in the ock precursor, they are bound to play a dominant role in the amplification of the magnetic field, most notably because the length scale of the cosmic ray distribution increases as r L. For instance, at a distance x from the ock front where only cosmic rays of momentum larger than p can be found, the critical wavenumber for the instability reads for s> s spectral index of the accelerated population): k = k, p p, ) s, 4) where k, denotes the same wavenumber evaluated for the whole cosmic-ray distribution, and p, denotes the minimum momentum of the cosmic ray distribution at the ock front. Consequently, if x k, a with a=/3 or / for the two branches of the instability obtained above: ) b+a s) p x l cr, = x,l cr, 43) p, with b= for large scale turbulence, b= for ort scale turbulence. Here as well, x, ould be understood as the growth rate derived previously for the whole cosmic-ray distribution. This product grows with increasing p /p, provided s<b+a)/a, i.e. s<8 b=, a=/3), or s<5 b=, a=/), or s<5/ b=, a=/3) or finally s<3 b=, a=/). These inequalities are likely to be fulfilled. Finally, the product k l MHD scales as p /p, ) s hence decreases if s>. This also means the typical scale of the instability increases with increasing p /p,. Of course, the above implicitly ignores the spatial dependence of B, but it suggests that amplification through the streaming of the highest energy part of the cosmic ray distribution will seed much more efficiently the magnetic field amplification Non-existence of an incompressible instability Since the above compressive instabilities appear to saturate at a level too small δb/b ) to allow Fermi acceleration over a broad range of energies, the possible existence of non-compressive instability becomes crucial. One can search for such modes by selecting the wave vectors accordingly, in which case the reduced system becomes: k D ũ z k x + k z k x k x k x + u ) z k x ũ β k z + k = x k x Dũ z k y ũ k 3 z z D ũ z = 44) x k x It is easy to verify that this system possesses no unstable solution but only damped solutions. Hence incompressible unstable modes do not exist in this MHD description. 4 CONCLUSIONS In this study, we have examined the amplification of a pre-existing magnetic field upstream of an ultra-relativistic ock wave. We have assumed that the magnetic field is fully perpendicular to the ock normal in the ock front frame, as generally expected in the ultra-relativistic limit. In the ock frame, the cosmic rays do not induce any current at zeroth order, only a net charge density, which is partly screened by the inflowing plasma. This charge distribution then triggers an instability on very ort spatial scales,

8 8 G. Pelletier, M. Lemoine and A. Marcowith with a growth rate increasing with the wavenumber. Cosmic rays do not respond to this amplification as it takes on scales much orter than the typical Larmor radius. It is important to stress that these instabilities are of a different nature than the resonant or non-resonant instability studied by Bell 4) in the case of a non-relativistic parallel ock wave, most notably because they are essentially compressive. Since density depletions are naturally limited, our linear analysis indicates that these instabilities saturate at a moderate level of amplification,δb/b. Furthermore, there is no unstable incompressible mode within the present MHD approximation which could provide a higher level of amplification. Therefore, we conclude that, within the framework of ideal MHD, instabilities at ultra-relativistic magnetized ock waves appear limited in their efficiency. Such instabilities cannot therefore account for the degree of amplification that has been inferred from the modeling of the afterglow radiation of gamma-ray bursts. We have also argued in Section. that the ratioδb/b sets the dynamic range of energies over which powerlaw spectra can be produced through Fermi acceleration. The present instabilities thus cannot allow successful Fermi acceleration. Nevertheless, it would certainly be useful to pursue the present investigation with dedicated numerical simlations, which would allow one to go beyond the present linear approximation and study whether non-linear effects can pu the small scale magnetic field to higher values. If the standard interpretation of gamma-ray burst afterglows as the synchrotron light of electrons accelerated at the forward ock but see Uhm & Beloborodov 6, Genet, Daigne & Mochkovitch 7 for recent alternatives involving the reverse ock) holds, one is led to conclude that some other instability is able to amplify the upstream magnetic field on ort spatial scales, or that some other form of acceleration mechanism is operating see for instance Hoino 8 for an alternative involving radiative pressure effects). The present work indicates that such cosmic ray induced instabilities would involve non-mhd effects. The ort spatial scales λ c r L / required are in conflict with the usual synchrotron resonance between cosmic rays and MHD waves; this remark is actually one of the motivation of the present work in a full MHD framework. However, one cannot rule out that other types of resonance could become relevant and yet involve incompressible modes, such as Alfvén waves. One promising instability, currently under study, involves a Cerenkov resonance ω=k x c) between plasma waves and the cosmic ray beam, generating a modified two stream instability. The growth rates appear quite promising and the details will be given in a forthcoming article. APPENDIX A: BACKGROUND AND PERTURBED QUANTITIES FOR PERPENDICULAR SHOCK WAVES A Stationary regime In the ock front frame, the electric field, current and charge density read: E = u B= u zb y e x u xb y e z, c j = 4π xb y e z, ρ pl = ρ cr + 4π xe x. A) The four velocity of the inflowing plasma is written c, uc) hence u=β is dimensionless). In this equation,represents the total Lorentz factor of the upstream plasma; it will be own to coincide nearly with in the following. Note that in the ock front, there is no current at zeroth order associated with the cosmic rays, only a net charge densityρ cr which may be partially screened by an induced non-zero plasma charge densityρ pl. Faraday s law implies in this stationary regime: E= x ux B y Hence the Gauss equation can be rewritten as: ρ pl = ρ cr + B y x u z u ) z x u x 4π u x ) =, A). A3) Energy momentum conservation for the upstream fluid implies: ρ u c u u=ρ pl E+ c j B, A4) whereρ u is the proper mass density of the upstream plasma. Denoting F m ρ u c u x the mass flux times c) through the ock front: F m x u x = ρ cr B y x u z u )] z uz B y x u x 4π u x + 4π B y x B y, F m x u z = ρ cr B y x u z u )] z ux B y x u x. A5) 4π u x Since:=+u x + u z )/, therefore A) implies: x B y = B y u x x = u x x u x + u z x u z u x 3, A6) x u x u ) xu z xu 3 z. A7) Using A7), the system A5) can be transformed into: B y u z B y F m 4π x u x = ρ cr u x B yu z 4π xu x + F m B yu x u x B y 4π x u z = ρ cr. A8) This sytem can be solved by successive approximations. For instance the ratio: B y 4π u x F m β A A9) since u x ;β A = B y u / 4πρ u c represents the Alfvén velocity as measured in the upstream rest frame). Hence: x u x Γ ρ cr u z B y F m. A) Consider now the ratio of the first term on the lhs of the second equation in A8) to the term of the rhs of the same equation: B x u yu z x 4π u z β Γ A. A) ρ cr u x B y

9 MHD-instabilities at ultra-relativistic magnetized ock waves 9 This ratio is much smaller than unity since u z ). Then: obtain: x u z ρ cr B y u z x) F m + B y 4π B y F m + B y 4π x x ρ cr dx. A) In the last equation, x corresponds to the maximal distance to which cosmic rays can stream ahead of the ock front, as measured in the ock front frame. We also assume that u x, B y andρ vary on much longer scales than u z, and that u z u x. This hierarchy can be verified by taking the ratios of the equations in A5). Since B y/4π) β A F m, one can simplify the previous expression and obtain the following order of magnitude: δe = B y δu z + u z δb y + β u z B y δu x + u Γ 3 z B y δu z )e x + β δb z Γ ) e y + δb x u z β δb y B Γ y β δu z u z B Γ 3 y δu x )e z. A7) u z l cr r L e cr ρ u c. A3) To obtain this result, we have approximated the integral overρ cr asl cr ρ cr,l cr representing the length scale of the distribution in the ock frame, and writtenρ cr e cr /p c); e cr denotes the cosmic ray energy density while p r L ) represents their typical momentum Larmor radius). Finally, e cr / ρ u c ) ξ cr whereξ cr is the fraction of ock internal energy carried away by the accelerated population. In a straightforward way, one obtains: so that: l cr r L O), u z u x ξ cr. This justifies the previous approximations. A Time dependent quantities A. Perturbed electric field The perturbed Ohm law leads to: δe = B y δu z + u z δb y + + )] Γ β u z δu x B y δu z u z B y e x + δb xu z +δb z u x ) e y + β δb y δu x B y + A4) A5) + Γ β δu xb y β δu z u z B y )]e z. A6) We keep track of the terms in u z in these equations in order to properly evaluate the terms which involve the derivative of u z with respect to x. This latter quantity is of the same order than other background quantities and cannot be neglected. However, terms in u z will be neglected in the final equations of evolution of the various background quantities since u z u x. For similar reasons, derivatives of background quantities other than u z ) such asρ pl, B y and u x hence) can be neglected. Thus, setting,β β we A. Perturbed equations The perturbed equations for the perturbed components of the magnetic field, defined in units of B y as b=δb/b y can be written: ˆDb x = Γ y δu x + u z z b x + β ) y δu z, ˆDb y = z δu z Γ x δu x β u z Γ δu z+ +u z z b y β z δu x β x δu z + + ) u Γ z z δu z, ˆDb z = y δu z + u z b x+ u z x b x + y b y + β y δu x + u ) z Γ y δu z. A8) The continuity equation can be expressed as: ˆD δw w = δu c tδ x δu x + y δu y + z δu z ). A9) We have noted w=ρ u c the proper enthalpy density of the upstream plasma. In the perturbation regime the contribution of the temporal derivative of the fluctuating part of the Lorentz factor is negligible. Again, while deriving Eq.A9) the x derivatives of zeroth order quantities have been discarded. The equation of motion in the ock rest frame accounting for the displacement current are: w ˆDu=ρ Pl E B 8π + B B 4π E B. A) 4πc t Using Eq.A7) and developping the Eq.A) to the first order we

10 G. Pelletier, M. Lemoine and A. Marcowith finally get after some manipulations: ˆDδu x ˆDδu y +β A) ˆDδu z = β A ) = β A Γ y b x x b y κδuz + +β A Γ β c tb y ) Γ 3 c tδu x β s δw x w + +u z κ β A ) b y + β κβ A Γ )δu x+ ) +β β A z b y y b z + + β β A β A ) z δu x + )] +β A x δu z β c tδu z, = β κ b z β s δw ) y κbx z δu y, w + u z ) y b z z b y β κ ) β A by + κβ ) β A δux β β A Γ zδu x + β s z δw w +β κ δw w + β +u z β Γ A) κδuz +β A) β z δu z + ) +β β A y b x x b y + +β β A c tδu x β x δu x +β A c tb y β x b y ) + )]. A) We used the orthand notationβ s = c s/c with c s the speed of sound, and: β A B y 4πΓ w, A) which defines the Alfvén velocity squared as measured upstream). The quantity κ, which carries the dimension of a wavenumber, is defined as follows: κ u z = +β A ρ cr B y w. A3) The plasma charge densityρ pl is indeed related to the cosmic-ray charge densityρ cr through the Maxwell equation: hence: ρ pl = ρ cr + B y 4π u z, ρ pl B y w = κ. A4) A5) In the limitβ A andβ s, the system governing the evolution of the perturbed quantities reads: ˆDb x = y δu x + u z z b x + β ) y δu z, Γ ˆDb y = z δu z Γ x δu x β κδu z + Γ +u z z b y β z δu x β ) x δu z, ˆDb z = y δu z +κb x + u z z b z + β ) y δu x, ) κ ˆDδu x = κδu z u z κb y +β δu x + z δu x ) ˆDδu y = β κ b z + u z κbx z δu y, ˆDδu z = β κb y + κ δu Γ x +β κ δw ) w + u β z κδu z z δu z, Γ ˆD δw w = x δu x y δu y z δu z. A6) To analyze the instabilities, it is convenient to make the following reduction: ˆD=k x β D,δu j =β ũ j ; b j and δw w unchanged. The ordering for doing simplifications is consistent with the expected physics: we assume k y, k z, k k x, u z β and Du z = k /k x at the same order than previous terms. The simplified system reads: Db x A7) Db y Db z z k ũ z β ũ z β x ũ z k x k x β k x y k x ũ z Dũ x β k k x ũ z Dũ y β k k x b z Dũ z k b y + k δw β k x β k x w D δw x ũ x + y ũ y + ) z ũ z w k x k x k x Thus one gets the intermediate system: D ũ z = D 3δw w = k k x REFERENCES k x β k x u z + k ) z u z k x +β β k ũ x u k z + k D δw x β k x w Dũ z k y ũ β k 3 z z D ũ z x k x, A8) A9) A3) A3) A3) A33) A34) Achterberg, A., Gallant, Y., Kirk, J. G., Guthmann, A. W.,, MNRAS 38, 393 Amato, E., & Blasi, P. 8, arxiv:86.3 Bell, A., 4, MNRAS, 353, 55 Bell, A., 5, MNRAS, 358, 8 Blandford, R., McKee, C., 976, Phys. Fluids, 9, 3 Chang, P., Spitkovsky, A., & Arons, J. 8, ApJ, 674, 378 Dieckmann, M. E., Shukla, P. K., Drury, L. O. C., 8, ApJ, 675, 586. Frederiksen, J. T., Hededal, C. B., Haugbølle, T., & Nordlund, Å. 4, ApJ, 68, L3 Gallant, Y., Achterberg, A., 999, MNRAS 35, L6 Genet, F., Daigne, F., Mochkovitch, R., 7, MNRAS 38, 73

11 Gruzinov, A., Waxman, E., 999, ApJ, 5, 85 Hededal, C. B., Niikawa, K. I., 5, ApJ, 63, L89 Hoino, M., 8, ApJ, 67, 94 Kato, T. N., 7, ApJ, 668, 974 Katz, J. I., 994, ApJ, 43, L7 Keet, U., Katz, B., Spitkovsky, A., Waxman E., 8, arxiv:8.37 Lemoine, M., Pelletier, G., 3, ApJ, 589, L73 Lemoine, M., Revenu, B., 6, MNRAS 366, 635 Lemoine, M., Pelletier, G., Revenu, B., 6, ApJ, 645, L9 Li, Z., Waxman, E., 6, ApJ, 65, L38 Lyubarsky, Y., Eichler, D., 6, ApJ, 647, L5 Marcowith, A., Lemoine, M., Pelletier, G., 6, AA, 453, 93 Medvedev, M. V., Loeb, A., 999, ApJ, 56, 697 Medvedev, M. V., Fiore, M., Fonseca, R. A., Silva, L. O., & Mori, W. B. 5, ApJ, 68, L75 Milosavljević, M., Nakar, E., 6, ApJ, 65, 979 Niemiec, J., Ostrowski, 6, ApJ, 64, 984 Niemiec, J., Ostrowski, M., Pohl, M., 6, ApJ, 65, Niemiec, J., Pohl, M., Stroman, T., & Niikawa, K. I. 8, ApJ, 684, 74 Mészáros, P., Rees, M., 997, ApJ, 476, 3 Parizot, E., Marcowith, A., Ballet, J. & Gallant, Y.A., 6, A&A, 453, 387 Paczyński, B., Rhoads, J. E., 993, ApJ 48, L5 Pelletier, G., Lemoine, M., Marcowith, A., 6, AA, 453, 8 Piran, T., 5, Rev. Mod. Phys., 76, 43 Reville, B., Kirk, J. G., Duffy, P., 6, Plasma Phys. Contr. Fus., 48, 74 Reville, B., Kirk, J. G., Duffy, P., O Sullivan, S., 7, AA, 475, 435 Reville, B., O Sullivan, S., Duffy, P., Kirk, J. G., 8, MNRAS, 386, 59 Sari, R., Piran, T. 997, ApJ, 485, 7 Silva, L. O., Fonseca, R. A., Tonge, J. W., Dawson, J. M., Mori, W. B., & Medvedev, M. V. 3, ApJ, 596, L Spitkovsky, A., 8, arxiv:8.36 Uhm, Z. L., Beloborodov, A. M., 6, ApJ 665, L93 Vietri, M. 997, ApJ, 478, L9 Völk, H.J., Berezhko, E.G. & Ksenofontov, L.T., 5, A&A, 433, 9 Waxman, E., 997, ApJ 485, L5 Wiersma, J., Achterberg, A., 4, AA, 48, 365 Zirakavili, V. N., Ptuskin, V. S., & Völk, H. J. 8, ApJ, 678, 55 MHD-instabilities at ultra-relativistic magnetized ock waves