Radiative diagnostics for sub-Larmor scale magnetic turbulence
Radiative diagnostics of high-energy density plasmas is addressed in this paper. We propose that the radiation produced by energetic particles in small-scale magnetic field turbulence, which can occur in laser-plasma experiments, collisionless shocks, and during magnetic reconnection, can be used to deduce some properties of the turbulent magnetic field. Particles propagating through such turbulence encounter locally strong magnetic fields, but over lengths much shorter than a particle gyroradius (that is, ). Consequently, the particle is accelerated but not deviated substantially from a straight line path. We develop the general jitter radiation solutions for this case and show that the resulting radiation is directly dependent upon the spectral distribution of the magnetic field through which the particle propagates. We demonstrate the power of this approach in considering the radiation produced by particles moving through a region in which a (Weibel-like) filamentation instability grows magnetic fields randomly oriented in a plane transverse to counterstreaming particle populations. We calculate the spectrum as would be seen from the original particle population and as could be seen by using a quasi-monoenergetic electron beam to probe the turbulent region at various angles to the filamentation axis.
Understanding the dynamics of small-scale magnetic field turbulence is of great interest in both laboratory and astrophysical plasmas due to its strong effect on particle distributions, as well as its involvement in the spontaneous generation, evolution, and dissipation of magnetic field. However, it remains difficult to relate the intense radiation produced by these sources to the underlying physical processes or field structure within them. Calculating the spectrum of such radiation poses unique analytical challenges due to the magnetic field varying rapidly and down to very small scales. Large-scale magnetic field structure and some details of the particle distribution in energy can be obtained through comparing spectra with synchrotron radiation; however, this fails where particle coupling to field turbulence is most intense. Furthermore, while the instabilities that generate such turbulence are often seen laser-plasma experiments, they can easily be disrupted or altered by the presence of in situ probes, and thus the formation and evolution of the turbulent magnetic field remains difficult to study quantitatively.
The theory of jitter radiation, named for the particle’s “jitter” as it undergoes a series of small transverse accelerations without a significant path deviation, has been developed (Medvedev, 2000, 2006) to improve upon this by considering the radiation emitted by a relativistic charged particle moving through a magnetic field that varies on spatial scales much smaller than a typical particle’s gyroradius. In this case the magnetic field correlation length is too small for a particle of Lorentz factor to be deflected beyond the opening angle into which most of its radiation will be relativistically beamed. This differs significantly from the synchrotron case in which the radiation spectrum is determined primarily by the sweep of the relativistic beaming cone past an observer’s direction. For a magnetic field with average strength perpendicular to the direction of particle motion and correlation length , we can thus define a parameter which determines the degree to which the radiation will be dominated by jitter or synchrotron radiation:
For , we have the synchrotron radiation case, in which the resulting spectrum is dominated by the frequency of the particles’ redirection due to the large-scale magnetic field. For , we have small-angle jitter radiation, in which the emission into a particular direction is produced by particles which undergo accelerations but not significant deflections, and thus radiate into the same direction over many field correlation lengths. In the intermediate region , the radiation spectrum is influenced both by the path geometry and the magnetic field distribution. In such a case, which we designate as large-angle jitter radiation, the particle is deflected beyond the relativistic beaming angle, but the field geometry along the path segment before this deflection occurs leads to the occurence of low-frequency harmonics that differ from the synchrotron. (See Medvedev (2010) for more on both the small- and large-angle jitter regimes and associated spectral features.) Undulator and “wiggler” radiation can be in the jitter regime, but are results for specific magnetic field configurations selected to coherently amplify certain harmonics, whereas the radiation produced by some distribution of particles moving in a region of turbulent plasma is generally an incoherent average over magnetic fields which randomly vary in both strength and orientation. Distinguishing jitter radiation from synchrotron radiation by spectral features alone is by no means trivial. However, for certain field configurations the jitter radiation spectrum can be notably harder at lower energies (below the spectral peak or break) than is possible from synchrotron radiation, which regardless of the particle distribution is limited to a low-energy spectral index of no greater than 1/3 (or -2/3 for the photon spectrum).
While synchrotron radiation reflects primarily the statistical properties of the particle distribution, the jitter radiation spectrum will directly reflect the statistical properties of the magnetic field along a particular line-of-sight. Jitter radiation is especially useful for experimental diagnostics because of its sensitivity to the magnetic field along a particular path and with a particular orientation to any field anisotropy. An ideal experimental design would make use of this by using a quasimonoenergetic electron beam to probe the site of plasma turbulence at a variety of incident angles and measure the resulting radiation, which should be emitted primarily along the direction of the probe beam. The observed radiation spectrum can then be compared with the theoretically calculated jitter spectrum, through the equations we develop in section II. In section III, we apply these equations to calculate the spectrum for a sample application in which a quasimonoenergetic beam of peak energy 200 MeV and FWHM energy spread of 50 MEV is used to probe a region of plasma current filamentation.
Jitter radiation solutions and spectral effects have thus far been largely developed and explored in application to astrophysical scenarios in which propagating shock fronts, magnetic reconnection, or other particle acceleration within the jetted outflow of highly relativistic collisionless plasma produce rapidly-evolving transient signals observed as gamma-ray bursts, X-ray flares, or other transient high-energy phenomena.Medvedev (2006); Reynolds, et al. (2010) Analytical treatments and PIC simulations indicate that such sources may produce an extended region of magnetic field turbulence generated by a relativistic Weibel-like filamentation instability. This instability has been called by a variety of names in the astrophysical literature and is often (incorrectly) referred to as the Weibel instability itself, but the basic mechanism is that small fluctuations in the magnetic field in a plane transverse to counterstreaming particles cause the particle streams to coalesce into current filaments which reinforce the initial field variations.Medvedev (2000); Gedalin et al. (2010); Bret et al. (2010) We will here refer to it simply as a filamentation instability, with the understanding that we presume the filamentation has formed in this way.
This filamentation instability can operate in a variety of scenarios and has been extensively studied in PIC simulations.Nishikawa, et al. (2003); Frederiksen et al. (2004); Bret et al. (2010). Much early work assumed such filamentation would occur at or upstream from a propagating shock front; however, recent results suggest that this type of instability may be much more ubiquitous, occuring also in non-shock plasmas with an appropriate anisotropy in particle distributions, and in outflow from sites of magnetic reconnection, among others.Norreys, et al. (2009)
While simulations clearly show the production of filamentation by such an instability, this exact mechanism of filamentation remains challenging to identify in an experimental setting. Certainly filamentary structures are seen in a wide variety of plasma scenarios, where they can degrade the efficiency of inertial confinement fusion and various laser and beam applications. Plasma density filamentation is frequently observed in laser-plasmas (see [citations] for instance) but such observations have yet to be reliably studied with the simultaneous in-situ magnetic field and current measurements that would be required to identify them as a particular structure within a coherent theory. The instabilities leading to this filamentation could provide valuable insights for laboratory astrophysics and high-energy density plasma physics, provided the instabilities can be adequately controlled and characterized. The jitter radiation approach presented in this article could help with such studies. Given a proposed form for the magnetic field wavenumber spectral distribution within a turbulent regime, we can calculate the resulting radiation spectra emitted towards various viewing directions for comparison with observations, allowing basic features of the magnetic field within the turbulent region to be confirmed without risk of disrupting or altering the instability with in situ probes.
In section II we develop the equations for calculating the (small-angle) jitter radiation spectrum for any given magnetic field configuration, and then specifically for the anisotropic field distributions as generated a filamentation instability. In section III we use these equations to calculate the jitter radiation spectra produced by particles moving in such an instability and by particles in an electron beam used to probe the filamentary region at various angles to the filamentation and separately characterize different components of the magnetic field.
Ii Angle-resolved jitter radiation theory
The angle-averaged jitter radiation emissivity equations have previously been developed and solutions presented for particles moving in the turbulent field produced by the relativistic Weibel instability within the internal shocks of a GRB jet.Medvedev (2000, 2006); Reynolds, et al. (2010) Here we present our derivation of the angle-resolved jitter emissivity such as may be used in laboratory applications.
The energy emitted per frequency per unit area for a relativistic charged particle with velocity that is undergoing an acceleration is given by the Leonard-Weichert equation (Landau & Lifshitz, 1971; Brau, 2004):
For a highly relativistic particle such as we are considering, the velocity is nearly constant and the time-dependence of the integrand is contained primarily in the acceleration . Using = , we define as the Fourier transform of the particle’s acceleration. The above equation becomes:
We can define the components of and in the direction of an observer as and and define and as the magnitudes of the vector projection of and on the plane perpendicular to . We then find the angular spectral fluence to be:
The accelerations undergone by the particle, generated by the Lorentz force, will lie in the plane perpendicular to the particle’s velocity. Thus, = 0, and the above reduces to:
For the case , this reduces to
To proceed further we need to determine the appropriate form for relating the accelerations of the particle to the underlying magnetic field. At a particular instant in time, a particle’s acceleration due to the Lorentz force is determined by the particle’s velocity and the local magnetic field vector . Recalling our earlier conditions, the particle velocity is nearly constant in magnitude and direction in the jitter regime. Thus, the temporal evolution of the particle’s acceleration reflects the propagation of the particle to a new position and thus becomes effectively an integral over a straight-line path of the particle through the turbulent magnetic field.
where we have defined .
where we have defined a tensor whose terms are products of the Fourier transformed magnetic field components and their complex conjugates, .
We assume that the field and particle distributions are sufficiently homogeneous that the acceleration spectrum can be calculated from the statistically-averaged field in the radiating region of volume V. That is, for a particular particle path across the radiating region, we presume that the exact starting position is irrelevant and exchange the integral of the field along a particular path for the volume-averaged magnetic field. For simplicity, we also assume that the time over which the particle is traversing the turbulent field region is much less than the dynamic time scale of the magnetic field turbulence, allowing us to treat the magnetic field distribution as static. Thus, the magnetic field spectrum correlation tensor becomes
In equations 5, 8, 10, and 11 we have not yet incorporated any particular geometry of the magnetic field, so these equations determine the jitter radiation spectra for any possible magnetic field distribution which varies on short enough scales to be within the jitter regime.
We now consider the geometry expected in a Weibel-type filamentation instability scenario, in which the magnetic field distribution along the filamentation axis (which we shall define here as ) is independent of the magnetic field distribution produced in the plane perpendicular to this axis. Rewriting the magnetic field correlation tensor into a form which implicity contains this geometry, we obtain
where we have defined the geometry and velocity dependent amplitude factor as
We can simplify our geometry slightly by considering the case in which , and lie in the same plane. In this case we can define the angles and such that , , and . In this planar case, Equation 14 becomes
Figure 1 shows this factor plotted versus for and , normalized to its value at =0. The strong peak of the emitted radiation in the forward direction is evident. Despite the functional asymmetry of terms dependent on in 15, we find that the calculated emission is symmetric to well below 0.01% for positive and negative .
In most situations, the magnetic field spectrum along and transverse to the filamentation axis cannot be known from first principles but may be approximated from some knowledge of the relevant parameters. The magnetic field spectrum transverse to the filamentation axis has been shown to rise and then drop after a scale of order the plasma skin depth (Frederiksen et al., 2004; Medvedev & Loeb, 1999). The filamentation instability itself produces this transverse field distribution independently at each point along the filamentation axis. As in Medvedev (2006) and Reynolds, et al. (2010), we define a geometry where the filamentation axis (i.e., the orientation of the currents from the original counterstreaming particle populations) is along the -direction and the cross-section across the filaments can be taken as lying locally in the -plane. To calculate the resulting spectra, we have used a general spectral form that may be parameterized as
where , , and are free parameters controlling the spectral break and the soft and hard spectral indices, respectively. In this equation is the cross-filamentation magnetic field wavenumber. Similarly, we characterize the magnetic field spectrum along z as
where , and are independent free parameters. We have defined our field spectra (Equations 17 and 18) so as to make variables and dimensionless values, expressed in terms relative to a wavenumber unit . This will result in a normalizable factor contributed to Equation 16 by each spectral distribution, where simply designates that we are considering either the parallel or perpendicular component of wavenumber. From the delta function in Equation 16 this means our frequencies will be defined in terms of units . The basic asymptotic behavior of the equations for the field spectra is summarized as
For the limiting cases in which the beam is aligned directly along () or transverse () to the filamentation axis the term within the delta function in Equation 16 reduces to a function of a single component of (with appropriate choice of the -axis) so that the integral over the delta function amounts to a simple substitution for or as a function of , , and . This results in a spectrum whose frequency dependence is determined only by the distribution of the magnetic field along (for ) or perpendicular (for ) to the filamentation axis. For a more general orientation of the beam at some angle relative to the filamentation axis, the substitution for or includes the other wavenumber component and the frequency dependence of the resulting spectrum is influenced by both magnetic field distributions.
Iii Results and Implications
A typical laser-plasma experiment for studies of the filamentation (Weibel-type) instabilities generally involves a strong beam that produces anisotropy in the plasma distribution function (PDF) which then drives the instability. The parallel (beam) direction is the direction of the PDF anisotropy, and it is the direction in which the current filaments are oriented. Radiation is produced by the plasma particles propagating in the magnetic fields of the filaments and it can be used to determine their properties. Alternatively, one can also launch a probe beam of electrons into the plasma and use the radiation produced by these particles to diagnose the magnetic fields. The latter scenario is more demanding in terms of engineering difficulties, but it should provide a more powerful diagnostics. Indeed, one generally have a better control on the energy distribution of the radiation-producing (probe beam) particles as well as one can probe a variety of incident angles relative to the filamentation direction. In both cases we base our radiating particle population on a quasi-monoenergetic electron beam consisting of a Gaussian distribution of particle energies, with peak energy of 200 MeV and a full-width at half maximum of 50 MeV.
For simplicity and in the absence of any particular plasma parameters governing our magnetic field distribution, we have chosen to use the same spectral parameters for the magnetic field distributions along and transverse to the filamentation axis. As we have already mentioned, the spectral shape is determined by the components of the magnetic field distribution along the beam so if there are distinctly different parameters for the magnetic field distribution along and perpendicular to the filamentation axis, the influence of these field parameter anisotropies can easily be isolated by looking at the emission from electrons moving along or transverse to the filamentation. We have used arbitrary-yet-reasonable selected values for the field spectral peak , low-wavenumber spectral indices of , and high-wavenumber spectral indices of . Our correspond to the peak of each field distribution in wavenumber, defined to be in terms of units . Defining our units in terms of the field correlation length in the direction transverse to the filamentation axis (the transverse field distribution being better developed by theory), we can define so that . We then have frequency units .
We find that in either case the emissivity is strongly peaked in the forward beam direction as expected. The spectral distribution of the generated spectrum is dependent primarily on the distribution of the magnetic field wavenumber’s component along the beam. Thus, the spectrum produced by the instability-generating beam contains details about the magnetic field distribution along the filamentation axis, as in subsection III.1. Obtaining information about the transverse magnetic field distribution requires a secondary beam to probe the instability region at angles transverse to the current filamentation, as we demonstrate in subsection III.2.
iii.1 Radiation from instability-generating beam
In this scenario, we consider the radiation produced by electrons in the same beam that produces the filamentation instability, slightly after the formation of the instability itself. The development of the filamentation instability is such that the filamentation axis will be aligned with the direction of the generating beam . In this case, we find and consequently, the components of perpendicular to the filamentation axis are eliminated from the delta function in Equation 16. The frequency dependence of the resulting radiation spectrum is thus obtained solely from the magnetic field distribution along the filamentation axis, .
Beam divergence will contribute a small influence from the transverse magnetic field to the spectrum emitted in the beam forward direction by the full particle distribution, but such contributions are strongly limited by the anisotropy of the particle distribution and the relativistic beaming factor as was shown in Figure 1. For simplicity, we also neglect here the self-consistent modifications of the beam as it propagates through the instability.
Figure 2 shows the resulting radiation spectrum produced by electrons within the same beam that produced the instability, when the spectrum are obtained from various angles = , and relative to the identical direction of the beam and the filamentation axis (, or equivalently ). The dashed line shows the spectrum produced by a particle of the electron beam’s peak energy and the solid line shows the approximate overall spectrum, calculated as a weighted sum over the spectra produced by particles of different . The values on the vertical axis are essentially arbitrary, as they depend on the average magnetic field strength generated by the instability, the electron beam particle density, and the total volume of the radiating region in which the beam intersects with the instability. Our horizontal axis is normalized to , which we have defined above as .
iii.2 Radiation from probe beam of varying incident angles
Exploring the magnetic field distribution transverse to the filamentation axis requires obtaining and analyzing the radiation produced by particle populations moving in these directions. We consider that the electron beam used to generate the instability may be split prior to entering the plasma (or a second beam generated) and redirected to probe the instability from other angles between the filamentation axis and the radiating beam. Since the peak emissivity will still be in the radiating particle populations forward direction, the radiation spectra would ideally be obtained with movable detectors that can be positioned in or close to the direction of the beam’s path through the filamentation region (outside of the filamentary region the beam may be deflected so as not to hit the detectors). Figure 3 shows the spectrum for incident beam angles of = , and when viewed from the beam’s forward direction (). As in the previous subsection, the vertical normalization is essentially arbitrary, but the horizontal normalization is defined in terms of , where is the correlation length of the magnetic field distribution transverse to the filamentation. The case is of course equivalent to that considered in subsection III.1 if there is no variation in the energy distribution between the probe beam and the beam that generated the instability.
For each , we have also calculated the spectrum for several different viewing angles and , where is again defined as the angle between the direction to an observer and the unit vector along the probe beam, as measured in the same plane as the angle on the far side of from . Figure 4 shows the results for probe beam incident angles of respectively. We have omitted the case from Figure 4 has been omitted since at these scales it appears nearly identical to the case (see Figure 3). As before, solid lines show the weighted sum over our distribution of , while dashed lines show the spectrum produced by an individual particle at the peak value of (here a representative 200/0.511 Mev). The spectrum from an individual particle and the full distribution differ notably in amplitude but not significantly in overall shape. For all of the probe beam angles, the spectrum is strongest when viewed along the probe beam (viewing angle ), with amplitudes decreasing rapidly by several orders of magnitude even when viewed at comparatively small viewing angles of and (in keeping with the relativistic beaming we expect). While the overall spectral shape is unchanged, the viewing angles also result in a shift of the spectrum to shorter wavelengths . This shift can be determined analytically from equation 13, from which we find that for the case a spectral feature located at in the will be shifted to when viewed at other angles .
We have developed the jitter radiation equations appropriate for analyzing the radiation produced within laboratory laser-plasma interactions in which relativistic particles radiate as they move through magnetic field turbulence without being substantially deviated from a path along a particular line of sight. Because the resulting jitter radiation spectrum depends directly upon the magnetic field distribution along the particle’s path, it has valuable potential as a non-invasive diagnostic for turbulent laboratory plasmas. Although we present results for a particular application in this paper, the approach developed here can be applied much more generally. Given a three-dimensional spectral distribution of the magnetic field wavenumber, we can directly calculate the resulting radiation spectrum in the jitter regime from particle beams propagating through this field distribution.
In this paper we have demonstrated this application of jitter radiation to the Weibel-like filamentation instability, which is of great interest for both astrophysical and laboratory plasmas. We calculate the radiation spectrum as would be produced by electrons in a quasi-monoenergetic beam that generates the filamentation instability as it impinges on a plasma. We demonstrate that the resulting radiation spectrum depends directly upon the magnetic field distribution along the filamentation axis, with very little sensitivity to the distribution of magnetic wavenumber transverse to this axis. The resulting spectrum is seen to be distinctly harder than the synchrotron case, with a low-frequency spectral index of 2. The spectrum peaks at higher frequencies in the beam’s forward direction, and then diminishes in overall intensity and in peak frequency when viewed from other angles.
We also consider a scenario in which a second quasi-monoenergetic electron beam probes the existing instability at varying incident angles relative to the instability’s filamentation axis. Taking advantage of this method’s control of density and orientation of the radiating particles, we can vary the angle to isolate the spectral influence of either of the independent magnetic field distributions along and transverse to the filamentation axis. As we vary the probe beams incident angle, we progress from a strongly peaked spectral form to peaked forms that have a low-frequency flattening, to unpeaked “broken power-law type” forms with low-frequency spectral index of 0 and a single break to the high-frequency spectral index.
This work has been supported by NSF grant AST-0708213, NASA ATFP grant NNX-08AL39G and DOE grant DE-FG02-07ER5494. We also gratefully acknowledge the computer system expertise of our colleague Sriharsha Pothapragada.
- Brau (2004) Brau, C., 2004, Modern Problems in Classical Electrodynamics (Oxford: Oxford University Press)
- Bret et al. (2010) Bret, A., Gremillet, L., & Dieckmann, M.E., 2010, Phys. Plasmas, 17, 120501
- Chang et al. (2008) Chang, P., Spitkovsky, A., & Arons, J., 2008, Ap.J., 674, 378
- Frederiksen et al. (2004) Frederiksen, J. T., Hededal, C. B., Haugbølle, T., & Nordlund, Å., 2004, Ap.J.Lett, 608, L13
- Gedalin et al. (2010) M. Gedalin, M. Medvedev, A. Spitkovsky, V. Krasnoselskikh, M. Balikhin, A. Vaivads, S. Perri, 2010, Phys. Plasmas, 17, 032108
- Landau & Lifshitz (1971) Landau, L., & Lifshitz, E. M. 1971, The classical theory of fields, (Oxford: Pergamon Press)
- Medvedev & Loeb (1999) Medvedev, M. V., & Loeb, A. 1999, Ap.J., 526, 697
- Medvedev (2000) Medvedev, M. V. 2000, Ap.J., 540, 704
- Medvedev (2006) Medvedev, M. V. 2006, Ap.J., 637, 869
- Medvedev & Spitkovsky (2009) Medvedev, M. V. & Spitkovsky, A. 2009, Ap.J., 700, 956
- Medvedev (2010) Medvedev, M.V., Frederiksen, J.T., Haugbølle, T., Nordlund, Å. submitted to Ap.J. (arXiv:1003.0063v2)
- Nishikawa, et al. (2003) Nishikawa, K.-I., Hardee, P., Richardson, G., Preece, R., Sol, H., & Fishman, G. J. 2003, Ap.J., 595, 555
- Norreys, et al. (2009) Norreys, P.A., Beg, F.N., Sentoku, Y., Silva, L.O., & Trines, R.M.G.M., 2009, Phys. Plasmas, 16, 041002
- Reynolds, et al. (2010) Reynolds, S.J., Pothapragada, S.S., & Medvedev, M.V., 2010, Ap.J., 713, 764
- Silva, et al. (2003) Silva, L. O., Fonseca, R. A., Tonge, J. W., Dawson, J. M., Mori, W. B., & Medvedev, M. V. 2003, Ap.J.Lett, 596, L121
- Spitkovsky (2005) Spitkovsky, A. 2005, in AIP Conf. Proc. 801, Astrophysical Sources of High Energy Particles and Radiation, ed. T. Bulik & B. Rudak (Melville: AIP), 345
- Spitkovsky (2007) Spitkovsky, A. 2007, Ap.J., 673, L39