skip to main content

Title: 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.
; ; ; ; ; ; ; ;
Award ID(s):
2019828 1909458
Publication Date:
Journal Name:
Journal of Plasma Physics
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 currentmore »sheets emerging during the interaction of Alfvén waves can potentially explain X-ray emission in black hole accretion coronae and neutron star magnetospheres.« less
  2. Abstract 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/3 for the case of a weak guide-field.

    Strong magnetic fields play an important role in powering the emission of neutron stars. Nevertheless, a full understanding of the interior configuration of the field remains elusive. In this work, we present general relativistic magnetohydrodynamics (MHD) simulations of the magnetic field evolution in neutron stars lasting ${\sim } {880}\,$ms (∼6.5 Alfvén crossing periods) and up to resolutions of $0.1155\,$km using Athena++. We explore two different initial conditions, one with purely poloidal magnetic field and the other with a dominant toroidal component, and study the poloidal and toroidal field energies, the growth times of the various instability-driven oscillation modes, and turbulence. We find that the purely poloidal setup generates a toroidal field, which later decays exponentially reaching $1{{\ \rm per\ cent}}$ of the total magnetic energy, showing no evidence of reaching equilibrium. The initially stronger toroidal field setup, on the other hand, loses up to 20 per cent of toroidal energy and maintains this state till the end of our simulation. We also explore the hypothesis, drawn from previous MHD simulations, that turbulence plays an important role in the quasi-equilibrium state. An analysis of the spectra in our higher resolution setups reveals, however, that in most cases we are not observing turbulence atmore »small scales, but rather a noisy velocity field inside the star. We also observe that the majority of the magnetic energy gets dissipated as heat increasing the internal energy of the star, while a small fraction gets radiated away as electromagnetic radiation.

    « less
  4. ABSTRACT Turbulent high-energy astrophysical systems often feature asymmetric energy injection: for instance, Alfvén waves propagating from an accretion disc into its corona. Such systems are ‘imbalanced’: the energy fluxes parallel and antiparallel to the large-scale magnetic field are unequal. In the past, numerical studies of imbalanced turbulence have focused on the magnetohydrodynamic regime. In this study, we investigate externally driven imbalanced turbulence in a collision-less, ultrarelativistically hot, magnetized pair plasma using 3D particle-in-cell (PIC) simulations. We find that the injected electromagnetic momentum efficiently converts into plasma momentum, resulting in net motion along the background magnetic field with speeds up to a significant fraction of lightspeed. This discovery has important implications for the launching of accretion disc winds. We also find that although particle acceleration in imbalanced turbulence operates on a slower time-scale than in balanced turbulence, it ultimately produces a power-law energy distribution similar to balanced turbulence. Our results have ramifications for black hole accretion disc coronae, winds, and jets.
  5. Abstract The formation of clusters at sub-saturation densities, as a result of many-body correlations, constitutes an essential feature for a reliable modelization of the nuclear matter equation of state (EoS). Phenomenological models that make use of energy density functionals (EDFs) offer a convenient approach to account for the presence of these bound states of nucleons when introduced as additional degrees of freedom. However, in these models clusters dissolve, by construction, when the nuclear saturation density is approached from below, revealing inconsistencies with recent findings that evidence the existence of short-range correlations (SRCs) even at larger densities. The idea of this work is to incorporate SRCs in established models for the EoS, in light of the importance of these features for the description of heavy-ion collisions, nuclear structure and in the astrophysical context. Our aim is to describe SRCs at supra-saturation densities by using effective quasi-clusters immersed in dense matter as a surrogate for correlations, in a regime where cluster dissolution is usually predicted in phenomenological models. Within the EDF framework, we explore a novel approach to embed SRCs within a relativistic mean-field model with density dependent couplings through the introduction of suitable in-medium modifications of the cluster properties, in particularmore »their binding energy shifts, which are responsible for describing the cluster dissolution. As a first exploratory step, the example of a quasi-deuteron within the generalized relativistic density functional approach is investigated. The zero temperature case is examined, where the deuteron fraction is given by the density of a boson condensate. For the first time, suitable parameterizations of the cluster mass shift at zero temperature are derived for all baryon densities. They are constrained by experimental results for the effective deuteron fraction in nuclear matter near saturation and by microscopic many-body calculations in the low-density limit. A proper description of well-constrained nuclear matter quantities at saturation is kept through a refit of the nucleon meson coupling strengths. The proposed parameterizations allow to also determine the density dependence of the quasi-deuteron mass fraction at arbitrary isospin asymmetries. The strength of the deuteron-meson couplings is assessed to be of crucial importance. Novel effects on some thermodynamic quantities, such as the matter incompressibility, the symmetry energy and its slope, are finally discerned and discussed. The findings of the present study represent a first step to improve the description of nuclear matter and its EoS at supra-saturation densities in EDFs by considering correlations in an effective way. In a next step, the single-particle momentum distributions in nuclear matter can be explored using proper wave functions of the quasi-deuteron in the medium. The momentum distributions are expected to exhibit a high-momentum tail, as observed in the experimental study of SRCs by nucleon knockout with high-energy electrons. This will be studied in a forthcoming publication with an extensive presentation of the theoretical method and the results.« less