Alfvén waves as excited in black hole accretion disks and neutron star magnetospheres are the building blocks of turbulence in relativistic, magnetized plasmas. A large reservoir of magnetic energy is available in these systems, such that the plasma can be heated significantly even in the weak turbulence regime. We perform high-resolution three-dimensional simulations of counter-propagating Alfvén waves, showing that an $E_{B_{\perp }}(k_{\perp }) \propto k_{\perp }^{-2}$ energy spectrum develops as a result of the weak turbulence cascade in relativistic magnetohydrodynamics and its infinitely magnetized (force-free) limit. The plasma turbulence ubiquitously generates current sheets, which act as locations where magnetic energy dissipates. We show that current sheets form as a natural result of nonlinear interactions between counter-propagating Alfvén waves. These current sheets form owing to the compression of elongated eddies, driven by the shear induced by growing higher-order modes, and undergo a thinning process until they break-up into small-scale turbulent structures. We explore the formation of current sheets both in overlapping waves and in localized wave packet collisions. The relativistic interaction of localized Alfvén waves induces both Alfvén waves and fast waves, and efficiently mediates the conversion and dissipation of electromagnetic energy in astrophysical systems. Plasma energization through reconnection in current sheets emerging during the interaction of Alfvén waves can potentially explain X-ray emission in black hole accretion coronae and neutron star magnetospheres.
more »
« less
Weak Alfvénic turbulence in relativistic plasmas. Part 1. Dynamical equations and basic dynamics of interacting resonant triads
Alfvén wave collisions are the primary building blocks of the non-relativistic turbulence that permeates the heliosphere and low- to moderate-energy astrophysical systems. However, many astrophysical systems such as gamma-ray bursts, pulsar and magnetar magnetospheres and active galactic nuclei have relativistic flows or energy densities. To better understand these high-energy systems, we derive reduced relativistic magnetohydrodynamics equations and employ them to examine weak Alfvénic turbulence, dominated by three-wave interactions, in reduced relativistic magnetohydrodynamics, including the force-free, infinitely magnetized limit. We compare both numerical and analytical solutions to demonstrate that many of the findings from non-relativistic weak turbulence are retained in relativistic systems. But, an important distinction in the relativistic limit is the inapplicability of a formally incompressible limit, i.e. there exists finite coupling to the compressible fast mode regardless of the strength of the magnetic field. Since fast modes can propagate across field lines, this mechanism provides a route for energy to escape strongly magnetized systems, e.g. magnetar magnetospheres. However, we find that the fast-Alfvén coupling is diminished in the limit of oblique propagation.
more »
« less
- NSF-PAR ID:
- 10332450
- Date Published:
- Journal Name:
- Journal of Plasma Physics
- Volume:
- 87
- Issue:
- 6
- ISSN:
- 0022-3778
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, i.e., catastrophic lateral instabilities triggered by a critical footpoint displacement of
ψ crit≳π . In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to 25% of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, ≲1043erg. Such events occur when significant energy builds up while deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind—for example, via relativistic magnetic reconnection or shocks. -
more » « less
ABSTRACT We study dynamics of relativistic coronal mass ejections (CMEs), from launching by shearing of foot-points (either slowly – the ‘Solar flare’ paradigm, or suddenly – the ‘star quake’ paradigm), to propagation in the preceding magnetar wind. For slow shear, most of the energy injected into the CME is first spent on the work done on breaking through the overlaying magnetic field. At later stages, sufficiently powerful CMEs may lead to the ‘detonation’ of a CME and opening of the magnetosphere beyond some equipartition radius req, where the decreasing energy of the CME becomes larger than the decreasing external magnetospheric energy. Post-CME magnetosphere relaxes via the formation of a plasmoid-mediated current sheet, initially at ∼req, and slowly reaching the light cylinder. Both the location of the foot-point shear and the global magnetospheric configuration affect the frequent/weak versus rare/powerful CME dichotomy – to produce powerful flares, the slow shear should be limited to field lines that close in near the star. After the creation of a topologically disconnected flux tube, the tube quickly (at ∼ the light cylinder) comes into force-balance with the preceding wind and is passively advected/frozen in the wind afterward. For fast shear (a local rotational glitch), the resulting large amplitude Alfvén waves lead to the opening of the magnetosphere (which later recovers similarly to the slow shear case). At distances much larger than the light cylinder, the resulting shear Alfvén waves propagate through the wind non-dissipatively.
-
more » « less
Abstract In a strongly magnetized, magnetically dominated relativistic plasma, Alfvénic turbulence can extend to scales much smaller than the particle inertial scales. It leads to an energy cascade somewhat analogous to inertial- or kinetic-Alfvén turbulent cascades existing in nonrelativistic space and astrophysical plasmas. Based on phenomenological modeling and particle-in-cell numerical simulations, we propose that the energy spectrum of such relativistic kinetic-scale Alfvénic turbulence is close to
k −3or slightly steeper than that due to intermittency corrections or Landau damping. We note the analogy of this spectrum with the Kraichnan spectrum corresponding to the enstrophy cascade in 2D incompressible fluid turbulence. Such turbulence strongly energizes particles in the direction parallel to the background magnetic field, leading to nearly one-dimensional particle momentum distributions. We find that these distributions have universal log-normal statistics. -
Context. Accreting black holes (BHs) may be surrounded by a highly magnetized plasma threaded by an organized poloidal magnetic field. Nonthermal flares and power-law spectral components at high energy could originate from a hot, collisionless, and nearly force-free corona. The jets we often observe from these systems are believed to be rotation-powered and magnetically driven. Aims. We study axisymmetric BH magnetospheres, where a fraction of the magnetic field lines anchored in a surrounding disk are connected to the event horizon of a rotating BH. For different BH spins, we identify the conditions and sites of magnetic reconnection within 30 gravitational radii. Methods. With the fully general relativistic particle-in-cell code GRZeltron , we solve the time-dependent dynamics of the electron–positron pair plasma and of the electromagnetic fields around the BH. The aligned disk is represented by a steady and perfectly conducting plasma in Keplerian rotation, threaded by a dipolar magnetic field. Results. For prograde disks around Kerr BHs, the topology of the magnetosphere is hybrid. Twisted open magnetic field lines crossing the horizon power a Blandford-Znajek jet, while open field lines with their footpoint beyond a critical distance on the disk could launch a magneto-centrifugal wind. In the innermost regions, coupling magnetic field lines ensure the transfer of significant amounts of angular momentum and energy between the BH and the disk. From the Y point at the intersection of these three regions, a current sheet forms where vivid particle acceleration via magnetic reconnection takes place. We compute the synchrotron images of the current sheet emission. Conclusions. Our estimates for jet power and BH–disk exchanges match those derived from purely force-free models. Particles are accelerated at the Y point, which acts as a heat source for the so-called corona. It provides a physically motivated ring-shaped source of hard X-rays above the disk for reflection models. Episodic plasmoid ejection might explain millisecond flares observed in Cygnus X-1 in the high-soft state, but are too fast to account for daily nonthermal flares from Sgr A * . Particles flowing from the Y point down to the disk could produce a hot spot at the footpoint of the outermost closed magnetic field line.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/S002237782100115X
