skip to main content

Title: Dynamic Alignment and Plasmoid Formation in Relativistic Magnetohydrodynamic Turbulence

We present high-resolution 2D and 3D simulations of magnetized decaying turbulence in relativistic, resistive magnetohydrodynamics. The simulations show dynamic formation of large-scale intermittent long-lived current sheets being disrupted into plasmoid chains by the tearing instability. These current sheets are locations of enhanced magnetic-field dissipation and heating of the plasma. We find magnetic energy spectra ∝k−3/2, together with strongly pronounced dynamic alignment of Elsässer fields and of velocity and magnetic fields, for strong guide-field turbulence, whereas we retrieve spectra ∝k−5/3for the case of a weak guide-field.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Letters
Medium: X Size: Article No. L13
["Article No. L13"]
Sponsoring Org:
National Science Foundation
More Like this

    We use 2.5D magnetohydrodynamic simulations to investigate the spectral signatures of the non-linear disruption of a tearing unstable current sheet via the generation of multiple secondary current sheets and magnetic islands. During the non-linear phase of tearing mode evolution, there develops a regime in which the magnetic energy density shows a spectrum with a power law close to B(k)2 ∼ k−0.8. Such an energy spectrum is found in correspondence of the neutral line, within the diffusion region of the primary current sheet, where energy is conveyed towards smaller scales via a ‘recursive’ process of fast tearing-type instabilities. Far from the neutral line, we find that magnetic energy spectra evolve towards slopes compatible with the ‘standard’ Kolmogorov spectrum. Starting from a self-similar description of the non-linear stage at the neutral line, we provide a model that predicts a reconnecting magnetic field energy spectrum scaling as k−4/5, in good agreement with numerical results. An extension of the predicted power law to generic current sheet profiles is also given and possible implications for turbulence phenomenology are discussed. These results provide a step forward to understand the ‘recursive’ generation of magnetic islands (plasmoids), which has been proposed as a possible explanation for the energy release during flares, but which, more in general, can have an impact on the subsequent turbulent evolution of unstable sheets that naturally form in the high Lundquist number and collisionless plasmas found in most of the astrophysical environments.

    more » « less
  2. Abstract

    We present a phenomenological and numerical study of strong Alfvénic turbulence in a magnetically dominated collisionless relativistic plasma with a strong background magnetic field. In contrast with the nonrelativistic case, the energy in such turbulence is contained in magnetic and electric fluctuations. We argue that such turbulence is analogous to turbulence in a strongly magnetized nonrelativistic plasma in the regime of broken quasi-neutrality. Our 2D particle-in-cell numerical simulations of turbulence in a relativistic pair plasma find that the spectrum of the total energy has the scalingk−3/2, while the difference between the magnetic and electric energies, the so-called residual energy, has the scalingk−2.4. The electric and magnetic fluctuations at scaleexhibit dynamic alignment with the alignment angle scaling close tocosϕ1/4. At scales smaller than the (relativistic) plasma inertial scale, the energy spectrum of relativistic inertial Alfvén turbulence steepens tok−3.5.

    more » « less
  3. Abstract

    The Parker Solar Probe (PSP) entered a region of sub-Alfvénic solar wind during encounter 8, and we present the first detailed analysis of low-frequency turbulence properties in this novel region. The magnetic field and flow velocity vectors were highly aligned during this interval. By constructing spectrograms of the normalized magnetic helicity, cross-helicity, and residual energy, we find that PSP observed primarily Alfvénic fluctuations, a consequence of the highly field-aligned flow that renders quasi-2D fluctuations unobservable to PSP. We extend Taylor’s hypothesis to sub- and super-Alfvénic flows. Spectra for the fluctuating forward and backward Elsässer variables (z±, respectively) are presented, showing thatz+modes dominatezby an order of magnitude or more, and thez+spectrum is a power law in frequency (parallel wavenumber)f−3/2(k3/2) compared to the convexzspectrum withf−3/2(k3/2) at low frequencies, flattening around a transition frequency (at which the nonlinear and Alfvén timescales are balanced) tof−1.25at higher frequencies. The observed spectra are well fitted using a spectral theory for nearly incompressible magnetohydrodynamics assuming a wavenumber anisotropykk3/4, that thez+fluctuations experience primarily nonlinear interactions, and that the minorityzfluctuations experience both nonlinear and Alfvénic interactions withz+fluctuations. The density spectrum is a power law that resembles neither thez±spectra nor the compressible magnetic field spectrum, suggesting that these are advected entropic rather than magnetosonic modes and not due to the parametric decay instability. Spectra in the neighboring modestly super-Alfvénic intervals are similar.

    more » « less
  4. Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially two-dimensional) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in three dimensions, we conduct a comprehensive particle-in-cell simulation study comparing two- and three-dimensional evolution of long, thin current sheets in moderately magnetized, collisionless, relativistically hot electron–positron plasma, and find dramatic differences. We first systematically characterize this process in two dimensions, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths and system sizes. We then show that three-dimensional (3-D) simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3-D current sheet evolution is not always fundamentally classical reconnection with perturbing 3-D effects but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details and even stochastic events. It often yields slower but longer-lasting and ultimately greater magnetic energy release than in two dimensions. Intriguingly, non-thermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in two and three dimensions. Although the variety of underlying current sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable. 
    more » « less
  5. Abstract

    We presentλ13 cm polarization observations of the nearby spiral galaxy NGC 6946 with the Westerbork Synthesis Radio Telescope (WSRT) to examine the nearside halo magnetic fields. Despiteλ13 cm exhibiting similar two-dimensional morphology as observed at longer (λ18–22 cm) or shorter (λ3 andλ6 cm) wavelengths, more complete frequency coverage will be required to explain the gap in polarization in the southwest quadrant of the galaxy. We fit models of the turbulent and coherent line-of-sight magnetic fields to the fractional degree of linearly polarized emission atλ3,λ6,λ13,λ18, andλ22 cm from observations taken with the WSRT, Karl G. Jansky Very Large Array, and Effelsberg telescopes. The results favor a multilayer turbulent magneto-ionized medium consistent with current observations of edge-on galaxies. We constrain the physical properties of the synchrotron-emitting thin and thick disks (scale heights of 300 pc and 1.4 kpc, respectively) along with the thermal thick disk and halo (scale heights of 1 and 5 kpc, respectively). Our preferred model indicates a clumpy and highly turbulent medium within 1 kpc of the midplane, and a diffuse extraplanar layer with a substantially lower degree of Faraday depolarization. In the halo, we estimate a regular magnetic field strength of 0.4–2.2μG and that turbulence and a total magnetic field strength of ∼6μG result in a Faraday dispersion ofσRM= 4–48 rad m−2. This work is an example of how the advanced capabilities of modern radio telescopes are opening a new frontier for the study of cosmic magnetism.

    more » « less