Abstract
Collisionless shocks can be produced as a result of strong magnetic fields in a plasma flow, and therefore are common in many astrophysical systems. The Weibel instability is one candidate mechanism for the generation of sufficiently strong fields to create a collisionless shock. Despite their crucial role in astrophysical systems, observation of the magnetic fields produced by Weibel instabilities in experiments has been challenging. Using a proton probe to directly image electromagnetic fields, we present evidence of Weibel-generated magnetic fields that grow in opposing, initially unmagnetized plasma flows from laser-driven laboratory experiments. Three-dimensional particle-in-cell simulations reveal that the instability efficiently extracts energy from the plasma flows, and that the self-generated magnetic energy reaches a few percent of the total energy in the system. This result demonstrates an experimental platform suitable for the investigation of a wide range of astrophysical phenomena, including collisionless shock formation in supernova remnants, large-scale magnetic field amplification, and the radiation signature from gamma-ray bursts.
Similar content being viewed by others
Main
The magnetic fields required for collisionless shock formation in astrophysical systems may either be initially present, for example in supernova remnants or young galaxies1, or they may be self-generated in systems such as gamma-ray bursts (GRBs; ref. 2). In the case of GRB outflows, the intense magnetic fields are greater than those which can be seeded by the GRB progenitor or produced by misaligned density and temperature gradients (the Biermann-battery effect)3,4. It has long been known that instabilities can generate strong magnetic fields, even in the absence of seed fields. Weibel considered the development of an electromagnetic instability driven by the electron velocity anisotropy in a background of resting ions5. The signature of the instability is a pattern of current filaments stretched along the axis of symmetry of the electron motion. The process is quite general, and subsequent work has shown that such instabilities can be excited in both non-relativistic and relativistic shocks. This general nature makes the Weibel instability common in astrophysical systems6,7,8. The instability provides a mechanism by which the electromagnetic turbulence associated with the formation of collisionless shocks is fed by the flow anisotropy of the protons (and ions) stochastically reflecting off of the shock9,10,11, and leading ultimately to strong particle acceleration in GRBâs (ref. 12).
The importance of Weibel instabilities in astrophysical systems makes laboratory experiments that can access the collisionless plasma regime particularly compelling, although so far experiments have been limited to collisional systems (where Weibel instability growth is limited by collisional dissipation13) or those where the initial plasma conditions are not well characterized14,15. Reaching the collisionless regime also permits the instability dynamics to be described by dimensionless parameters and scaled between laboratory and astrophysical systems16. In the collisionless regime, the scaling is remarkably simple and allows one to predict the parameters of the unstable modes and the shocks (should they be formed) on the basis of laboratory measurements and the astrophysical âinputâ parameters, the density and velocity of the flows.
In experiments performed at the Omega Laser Facility17, we directly image the magnetic fields associated with the Weibel instability in well-characterized, counter-streaming plasma flows in the collisionless plasma regime18. The flows were established by laser ablation of opposing foils, as shown in Fig. 1. The foils were oriented opposite each other and irradiated simultaneously, such that the expanding plasma flows interacted near the midplane between the foils19,20,21. The plasma conditions in this geometry have been previously measured under identical conditions with Thomson scattering18. When only a single foil was used, the plasma flow velocity v was measured to be 1,000â2,000 km sâ1, with an electron density (ne) â 5 Ã 1018 cmâ3 and electron and ion temperatures (Te, Ti) less than 200 eV.
The experiment consists of a pair of (CH2) plastic foils of diameter 2 mm and thickness 500 μm, oriented face-on and separated by 8 mm. Each was irradiated with eight overlapped laser beams, delivering â¼4 kJ of 351 nm laser energy in a 1 ns square pulse. Distributed phase plates were used to produce super-Gaussian laser spots with focal spot diameters of 250 μm on the target surface. After a delay, the proton probe was created by laser-compressing a thin-walled SiO2 capsule. The capsule was filled with a 1:1 mixture of deuterium (D) and 3helium (3He) at a total pressure of 18 atm. At peak compression (1023 cmâ3) protons are produced quasi-isotropically at energies of 3.0 and 14.7 MeV. The protons were detected using a CR39 nuclear track detector positioned on the midplane of the CH2 target foils, such that the protons traverse the central interaction region as shown.
When two opposing foils are used, as in the present work, the plasma density in the counter-propagating flows increased by the anticipated factor of two, whereas the electron and ion temperatures increased rapidly owing to a combination of collisional electron heating and ion two-stream instability20. The ion instability quickly stabilizes as the electron and ion temperatures equilibrate, which was observed to occur near 1 keV. The ions remain directed throughout the process, allowing competing instabilities, including the Weibel instability, to grow from the energy supplied by the flows22.
Magnetic fields are detected in this experiment using proton imaging. An isotropically emitting proton source is generated by the implosion of a capsule filled with Dâ3He, producing protons at 3 MeV (from DâD reactions) and at 14.7 MeV (from Dâ3He reactions; see Supplementary Information for additional details on proton imaging). The protons that pass through the plasma interaction region are deflected by the electric and magnetic fields in the system, and are recorded using a CR39 nuclear track detector at a magnification of approximately 30. There are several important features in the proton radiography data, which was taken at three different times during the interaction of the flows, and are shown on the top two rows of Fig. 2. First, oriented along the flow direction is a pattern of filamentary structures, consistent with Weibel filamentation in the counter-propagating flows. These features develop strongly between 3 and 4 ns after the start of the drive laser pulse and grow to lengths >1 mm along the direction of the flow. The filamentary structure is clear in both 3.0 and 14.7 MeV proton images, and extends relatively uniformly for several mm across the plasma flow, occupying a total volume of several mm3. The similarity in the observed features and relative contrast between the 3 MeV and 14.7 MeV radiographs indicate that proton deflections were produced by magnetic fields (see Supplementary Information for further discussion).
Experimental proton images are shown from 3.0 MeV DâD protons (top row), 14.7 MeV Dâ3He protons (middle row), and synthetic 14.7 MeV proton tracing from 3D particle-in-cell (PIC) simulations (bottom row). In each case the plasma flows enter the frame from the top and bottom. The small timing difference between DâD and Dâ3He images is due to the proton time of flight from the capsule to the plasma interaction region. At early time (approximately 3 ns after the drive begins), only initial traces of filamentation are observed. At later times the filaments become more coherent and increase in extent along the flow direction. In each case extended magnetic âplatesâ are formed above and below the midplane as a result of the large-scale Biermann-battery fields generated in the laser-ablation process23. All images are 3 mm in diameter in the target plane.
In addition to the filaments, horizontal âplateâ features are seen near the midplane of the drive plasmas. These large-scale magnetic features have been observed in previous experiments with similar geometries14,19, and are understood to be the result of Biermann-battery-generated magnetic fields23. These fields are created at the target surface during the laser ablation and form a loop around the expanding plasma flow24,25. The Biermann fields are frozen in the flow, following the effective electron trajectory to the midplane, where the longitudinal electron velocity from the two flows is cancelled. The magnetic fields cannot readily cross the midplane, and expand transversely, leading to the formation of characteristic plates23. Asymmetry between the top and bottom plates in the data is related to slight non-uniformities in the flows, including differences in laser energy deposition on the two foils and tilts in their orientation relative to the proton probe.
To better understand both the Weibel and Biermann-battery-generated magnetic fields in the experiment we have conducted detailed 3-dimensional particle-in-cell (3D PIC) simulations with the code OSIRIS (refs 26, 27) to model, from first principles, the counter-streaming plasma flows and the generation of electromagnetic fields (Fig. 3a, b). The flows are initialized with the properties measured experimentally in the midplane regionânamely, each flow has ne = 5 Ã 1018 cmâ3, v = 1,900 km sâ1 and Te = Ti = 1 keV. To include the effect of the Biermann battery, the flows were encircled by a large-scale magnetic field consistent with the misaligned density and temperature gradients of the flow, with an initial peak amplitude of 50 kG (see ref. 19). Further simulation details are found in the Supplementary Information.
a, 3D OSIRIS simulation of the system after 1 ns of interaction between the counter-streaming 1,900 km sâ1 plasma flows (approximately 3 ns after the experimental drive laser pulse; flows enter from top and bottom). Magnetic fields are shown qualitatively in the blue/red colour scale, with electron density in orange. b, Magnetic field slice (transverse magnetic field component By) along the y-axis midplane, at the same time, illustrating the presence of strong filaments associated with the Weibel instability. c, Plasma magnetization, Ï, as a function of time. When the flows are initiated with zero initial magnetic field (dashed lines) the magnetizations remain at zero until the flows begin interacting, between 2 and 3 ns. When initial toroidal fields are included consistent with the Biermann-battery mechanism, the perpendicular magnetization is â¼0.1% before the flows interact (solid coloured lines). In both cases the magnetic energy associated with Weibel instability increases sharply after the flows interact, increasing Ï by a factor greater than ten in several ns. The magnetization due to the ion Weibel instability, growing at the theoretical linear growth rate, is shown in solid black. This calculation shows that the Weibel-generated magnetization becomes the dominant contribution to the overall magnetization of the system. d, Measurement of the mean separation between filaments in experimental proton radiographs (red) and synthetic proton images from 3D PIC simulations (blue). The filament spacing approximately doubles over the 2 ns of observation. Note that time is experimental time, measured with respect to the beginning of the drive laser.
Within 1 ns of the opposing flows beginning to interact at the midpoint of the simulation volume, magnetic filaments are generated via the Weibel instability (Fig. 3b). Also, the magnitude of the toroidal Biermann-battery magnetic fields that are initially imposed doubles, owing to the conservation of magnetic flux. These fields lead to a long-range order in the system, and generate a pair of magnetic plates similar to those seen in the experiment (Fig. 2). The presence of the toroidal fields does not significantly alter the formation of the ion Weibel instability because the ions remain unmagnetized. This is supported by simulations where, when the initial magnetic fields are not included, Weibel filaments are still produced with the same structure.
To properly compare the PIC results with the experimental radiographs, we have simulated the proton imaging directly in the 3D OSIRIS simulations to generate proton images of the electric and magnetic fields in the system. We assume an isotropic point source of 14.7 MeV protons located 1 cm from the centre of the simulation box. The diagnostic protons interact with the 3D field structure during the simulation, and are then ballistically propagated to a 13 cm à 13 cm detector plane 30 cm from the source, matched to the imaging system in the experiment.
The simulated proton radiographs are shown on the bottom row of Fig. 2, at the same times as the experimental data. To quantify the evolving structure in the system the filament spacing was measured for all images. Shown in Fig. 3d, the growth in the size of the Weibel features is seen to be nearly equal for the measured and simulated images. The increasing filament size indicates growing Weibel fields, and the efficiency of the instability to convert system kinetic energy into magnetic energy can be assessed directly from the simulations. The magnetic energy associated with the instability is driven by the ion flows, and goes mainly into the transverse component of the field. The amplitude of the Weibel magnetic fields grows exponentially during the linear phase, with a growth rate of â¼0.2 v/c à Ïpi (where c is the speed of light and Ïpi is the ion plasma frequency), which is consistent with the linear theory of the instability. The linear phase of the instability saturates after 1â1.5 ns of interaction between the flows (that is, after â¼2â3 e-foldings), although the field amplitude and filament size continue to increase in the subsequent nonlinear phase. The significance of the magnetic energy in the system is quantified by the magnetization parameter Ï, defined as Ï â ãB2ã/4Ïminiv2, where the spatially averaged field is given by ãB2ã, with mi and ni representing the ion mass and density, respectively. This ratio of magnetic energy to flow kinetic energy reaches nearly 1% by the end of the experimental interaction time (approximately 5 ns in Fig. 3c).
At times later than those probed in the experiment, the local magnetic field strength peaks at 0.6 MG, which corresponds to Ï = 5%. These high values illustrate the efficiency of the Weibel instability in converting kinetic energy into electromagnetic energy. At this amplitude the magnetic fields are large enough to cause the deflection of the incoming flows and the randomization of their kinetic energy. Furthermore, Weibel-mediated collisionless shocks are believed to occur at this level of magnetization, provided that there is a sufficiently large interpenetration distanceâof the order of 300 c/Ïpi (ref. 9). This condition precludes shock formation in the present experiment, where this length is only 55 c/Ïpi; however, a shock should be observed with a similar set-up if higher densities and/or longer flows are generated.
The dimensionless magnetization parameter Ï serves to connect experiments to astrophysical systems, and can be directly applied to scaled systems of interest16. For example, multiwavelength observations of afterglow emission of GRBs suggest sub-equipartition levels of magnetic field in the region behind the forward shock (Ï â 0.01â0.1; ref. 28). Such fields are probably generated near a collisionless shock front, as the pre-shock interstellar medium is effectively unmagnetized (Ï â 10â10). Our experiments imply that even non-relativistic Weibel instability in an initially unmagnetized medium is capable of generating the percent-level magnetization observed in GRBs; collisions of relativistic flows are likely to produce even stronger fields.
References
Bernet, M. L., Miniati, F., Lilly, S. J., Kronberg, P. P. & Dessauges-Zavadsky, M. Strong magnetic fields in normal galaxies at high redshift. Nature 454, 302â304 (2008).
Medvedev, M. V., Silva, L. O., Fiore, M., Fonseca, R. A. & Mori, W. B. Generation of magnetic fields in cosmological shocks. J. Korean Astron. Soc. 37, 533â541 (2004).
Medvedev, M. V. & Loeb, A. Generation of magnetic fields in the relativistic shock of gamma-ray burst sources. Astrophys. J. 526, 697â706 (1999).
Gruzinov, A. & Waxman, E. Gamma-ray burst afterglow: Polarization and analytic light curves. Astrophys. J. 511, 852â861 (1999).
Weibel, E. S. Spontaneously growing transverse waves in a plasma due to an anisotropic velocity distribution. Phys. Rev. Lett. 2, 83â84 (1959).
Moiseev, S. S. & Sagdeev, R. Z. Collisionless shock waves in a plasma in a weak magnetic field. J. Nucl. Energy C 5, 43â47 (1963).
Lazar, M., Schlickeiser, R., Wielebinski, R. & Poedts, S. Cosmological effects of Weibel-type instabilities. Astrophys. J. 693, 1133â1141 (2009).
Medvedev, M. V. Electron acceleration in relativistic gamma-ray burst shocks. Astrophys. J. Lett. 651, L9âL11 (2006).
Kato, T. N. & Takabe, H. Nonrelativistic collisionless shocks in unmagnetized electronâion plasmas. Astrophys. J. Lett. 681, L93âL96 (2008).
Spitkovsky, A. On the structure of relativistic collisionless shocks in electronâion plasmas. Astrophys. J. Lett. 673, L39âL42 (2008).
Martins, S. F., Fonseca, R. A., Silva, L. O. & Mori, W. B. Ion dynamics and acceleration in relativistic shocks. Astrophys. J. Lett. 695, L189âL193 (2009).
Spitkovsky, A. Particle acceleration in relativistic collisionless shocks: Fermi process at last? Astrophys. J. Lett. 682, L5âL8 (2008).
Ryutov, D. D., Fiuza, F., Huntington, C. M., Ross, J. S. & Park, H-S. Collisional effects in the ion Weibel instability for two counter-propagating plasma streams. Phys. Plasmas 21, 032701 (2014).
Kugland, N. L. et al. Self-organized electromagnetic field structures in laser-produced counter-streaming plasmas. Nature Phys. 8, 809â812 (2012).
Fox, W. et al. Filamentation instability of counterstreaming laser-driven plasmas. Phys. Rev. Lett. 111, 225002 (2013).
Ryutov, D. D. et al. Basic scalings for collisionless-shock experiments in a plasma without pre-imposed magnetic field. Plasma Phys. Control. Fusion 54, 105021 (2012).
Boehly, T. R. et al. The upgrade to the Omega laser system. Rev. Sci. Instrum. 66, 508â510 (1995).
Ross, J. S. et al. Characterizing counter-streaming interpenetrating plasmas relevant to astrophysical collisionless shocks. Phys. Plasmas 19, 056501 (2012).
Kugland, N. L. et al. Visualizing electromagnetic fields in laser-produced counter-streaming plasma experiments for collisionless shock laboratory astrophysics. Phys. Plasmas 20, 056313 (2013).
Ross, J. S. et al. Collisionless coupling of ion and electron temperatures in counterstreaming plasma flows. Phys. Rev. Lett. 110, 145005 (2013).
Takabe, H. et al. High-mach number collisionless shock and photo-ionized non-LTE plasma for laboratory astrophysics with intense lasers. Plasma Phys. Control. Fusion 50, 124057 (2008).
Kato, T. N. & Takabe, H. Electrostatic and electromagnetic instabilities associated with electrostatic shocks: Two-dimensional particle-in-cell simulation. Phys. Plasmas 17, 032114 (2010).
Ryutov, D. D. et al. Magnetic field advection in two interpenetrating plasma streams. Phys. Plasmas 20, 032703 (2013).
Gregori, G. et al. Generation of scaled protogalactic seed magnetic fields in laser-produced shock waves. Nature Lett. 481, 480â483 (2012).
Schoeffler, K. M., Loureiro, N. F., Fonseca, R. A. & Silva, L. O. Magnetic-field generation and amplification in an expanding plasma. Phys. Rev. Lett. 112, 175001 (2014).
Fonseca, R. et al. in Computational Science, ICCS 2002 Vol. 2331 (eds Sloot, P., Hoekstra, A., Tan, C. & Dongarra, J.) 342â351 (Springer, 2002).
Fonseca, R. A. et al. One-to-one direct modeling of experiments and astrophysical scenarios: Pushing the envelope on kinetic plasma simulations. Plasma Phys. Control. Fusion 50, 124034 (2008).
Frail, D. A., Waxman, E. & Kulkarni, S. R. A 450 day light curve of the radio afterglow of GRB 970508: Fireball calorimetry. Astrophys. J. 537, 191â204 (2000).
Acknowledgements
We thank the staff of the Omega Laser Facility for their experimental support. This work was performed under the auspices of the US Department of Energy by the Lawrence Livermore National Laboratory, under Contract No. DE-AC52-07NA27344, with funding support from LLNL LDRD grant No. 11-ERD-054 and from the European Research Council under the European Communityâs Seventh Framework Programme (FP7/2007-2013), ERC grant agreement no. 256973. Computing support for this work came from ALCC and INCITE awards on Mira (ALCF supported under contract DE-AC02-06CH11357) and from the LLNL Institutional Computing Grand Challenge program on Vulcan. Additionally, the authors would like to acknowledge the OSIRIS Consortium, consisting of UCLA and IST (Lisbon, Portugal) for the use of the OSIRIS 2.0 framework and the visXD framework. F.F. acknowledges the LLNL Lawrence Fellowship for financial support. A.S. is supported by DOE grant DE-NA0002200.
Author information
Authors and Affiliations
Contributions
This experimental platform was conceived by G.G., B.A.R., R.P.D., H.-S.P., D.H.F., Y.S., A.S. and H.T. It was developed in a series of initial experiments by N.L.K., C.C.K., J.M., T.M., C.P., J.S.R. and C.M.H., and the data shown here were collected in experiments by C.M.H., J.S.R. and H.-S.P. Theoretical support for this work was provided primarily by F.F. and D.D.R., as well as by A.S. and G.G. Support for proton diagnostics was provided by C.K.L., R.P., and particularly A.B.Z. Particle-in-cell simulations were performed and synthetic radiographs produced by F.F. Additional contributions to the analysis of proton radiographs by H.-S.P., J.S.R., M.C.L., A.B.Z., D.D.R., A.S. and C.M.H. The manuscript was written by C.M.H., F.F., J.S.R. and D.D.R.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Information (PDF 414 kb)
Rights and permissions
About this article
Cite this article
Huntington, C., Fiuza, F., Ross, J. et al. Observation of magnetic field generation via the Weibel instability in interpenetrating plasma flows. Nature Phys 11, 173â176 (2015). https://doi.org/10.1038/nphys3178
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys3178
This article is cited by
-
Stream instabilities in optical-field ionization of a monatomic dilute neutral gas in fully relativistic regime
Scientific Reports (2023)
-
Self-organization of photoionized plasmas via kinetic instabilities
Reviews of Modern Plasma Physics (2023)
-
Electron energization dynamics in interaction of self-generated magnetic vortices in upstream of collisionless electron/ion shocks
Scientific Reports (2022)
-
Electron acceleration in laboratory-produced turbulent collisionless shocks
Nature Physics (2020)
-
Growth of concomitant laser-driven collisionless and resistive electron filamentation instabilities over large spatiotemporal scales
Nature Physics (2020)