skip to main content

Title: Kinetic simulations of imbalanced turbulence in a relativistic plasma: Net flow and particle acceleration
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.
; ; ; ;
Award ID(s):
1903335 1806084
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
3826 to 3841
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We study the flow structure in 3D magnetohydrodynamic (MHD) simulations of accretion on to Sagittarius A* via the magnetized winds of the orbiting Wolf–Rayet stars. These simulations cover over 3 orders of magnitude in radius to reach ≈300 gravitational radii, with only one poorly constrained parameter (the magnetic field in the stellar winds). Even for winds with relatively weak magnetic fields (e.g. plasma β ∼ 106), flux freezing/compression in the inflowing gas amplifies the field to β ∼ few well before it reaches the event horizon. Overall, the dynamics, accretion rate, and spherically averaged flow profiles (e.g. density, velocity) in our MHD simulations are remarkably similar to analogous hydrodynamic simulations. We attribute this to the broad distribution of angular momentum provided by the stellar winds, which sources accretion even absent much angular momentum transport. We find that the magneto-rotational instability is not important because of (i) strong magnetic fields that are amplified by flux freezing/compression, and (ii) the rapid inflow/outflow times of the gas and inefficient radiative cooling preclude circularization. The primary effect of magnetic fields is that they drive a polar outflow that is absent in hydrodynamics. The dynamical state of the accretion flow found in our simulationsmore »is unlike the rotationally supported tori used as initial conditions in horizon scale simulations, which could have implications for models being used to interpret Event Horizon Telescope and GRAVITY observations of Sgr A*.« less
  2. ABSTRACT High-energy astrophysical systems frequently contain collision-less relativistic plasmas that are heated by turbulent cascades and cooled by emission of radiation. Understanding the nature of this radiative turbulence is a frontier of extreme plasma astrophysics. In this paper, we use particle-in-cell simulations to study the effects of external inverse Compton radiation on turbulence driven in an optically thin, relativistic pair plasma. We focus on the statistical steady state (where injected energy is balanced by radiated energy) and perform a parameter scan spanning from low magnetization to high magnetization (0.04 ≲ σ ≲ 11). We demonstrate that the global particle energy distributions are quasi-thermal in all simulations, with only a modest population of non-thermal energetic particles (extending the tail by a factor of ∼2). This indicates that non-thermal particle acceleration (observed in similar non-radiative simulations) is quenched by strong radiative cooling. The quasi-thermal energy distributions are well fit by analytic models in which stochastic particle acceleration (due to, e.g. second-order Fermi mechanism or gyroresonant interactions) is balanced by the radiation reaction force. Despite the efficient thermalization of the plasma, non-thermal energetic particles do make a conspicuous appearance in the anisotropy of the global momentum distribution as highly variable, intermittent beams (formore »high magnetization cases). The beamed high-energy particles are spatially coincident with intermittent current sheets, suggesting that localized magnetic reconnection may be a mechanism for kinetic beaming. This beaming phenomenon may explain rapid flares observed in various astrophysical systems (such as blazar jets, the Crab nebula, and Sagittarius A*).« less
  3. ABSTRACT Discs are essential to the formation of both stars and planets, but how they form in magnetized molecular cloud cores remains debated. This work focuses on how the disc formation is affected by turbulence and ambipolar diffusion (AD), both separately and in combination, with an emphasis on the protostellar mass accretion phase of star formation. We find that a relatively strong, sonic turbulence on the core scale strongly warps but does not completely disrupt the well-known magnetically induced flattened pseudo-disc that dominates the inner protostellar accretion flow in the laminar case, in agreement with previous work. The turbulence enables the formation of a relatively large disc at early times with or without AD, but such a disc remains strongly magnetized and does not persist to the end of our simulation unless a relatively strong AD is also present. The AD-enabled discs in laminar simulations tend to fragment gravitationally. The disc fragmentation is suppressed by initial turbulence. The AD facilitates the disc formation and survival by reducing the field strength in the circumstellar region through magnetic flux redistribution and by making the field lines there less pinched azimuthally, especially at late times. We conclude that turbulence and AD complement eachmore »other in promoting disc formation. The discs formed in our simulations inherit a rather strong magnetic field from its parental core, with a typical plasma-β of order a few tens or smaller, which is 2–3 orders of magnitude lower than the values commonly adopted in magnetohydrodynamic simulations of protoplanetary discs. To resolve this potential tension, longer term simulations of disc formation and evolution with increasingly more realistic physics are needed.« less
  4. Abstract The magneto-rotational instability (MRI) is the most likely mechanism for transportation of angular momentum and dissipation of energy within hot, ionized accretion discs. This instability is produced through the interactions of a differentially rotating plasma with an embedded magnetic field. Like all substances in nature, the plasma in an accretion disc has the potential to become magnetically polarized when it interacts with the magnetic field. In this paper we study the effect of this magnetic susceptibility, parameterized by χm, on the MRI, specifically within the context of black hole accretion. We find from a linear analysis within the Newtonian limit that the minimum wavelength of the first unstable mode and the wavelength of the fastest growing mode are shorter in paramagnetic (χm > 0) than in diamagnetic (χm < 0) discs, all other parameters being equal. Furthermore, the magnetization parameter (ratio of gas to magnetic pressure) in the saturated state should be smaller when the magnetic susceptibility is positive than when it is negative. We confirm this latter prediction through a set of numerical simulations of magnetically polarized black hole accretion discs. We additionally find that the vertically integrated stress and mass accretion rate are somewhat larger when themore »disc is paramagnetic than when it is diamagnetic. If astrophysical discs are able to become magnetically polarized to any significant degree, then our results would be relevant to properly interpreting observations.« less

    Supernova remnants (SNRs) are known to accelerate cosmic rays from the detection of non-thermal emission in radio waves, X-rays, and gamma-rays. However, the ability to accelerate cosmic rays up to PeV energies has yet to be demonstrated. The presence of cut-offs in the gamma-ray spectra of several young SNRs led to the idea that PeV energies might only be achieved during the first years of a remnant’s evolution. We use our time-dependent acceleration-code RATPaC to study the acceleration of cosmic rays in supernovae expanding into dense environments around massive stars. We performed spherically symmetric one-dimensional (1D) simulations in which we simultaneously solve the transport equations for cosmic rays, magnetic turbulence, and the hydrodynamical flow of the thermal plasma in the test-particle limit. We investigated typical circumstellar-medium (CSM) parameters expected around red supergiant (RSG) and luminous blue variable (LBV) stars for freely expanding winds and accounted for the strong γγ absorption in the first days after explosion. The maximum achievable particle energy is limited to below $600\,$TeV even for the largest considered values of the magnetic field and mass-loss rates. The maximum energy is not expected to surpass $\approx 200\,$ and $\approx 70\,$TeV for LBVs and RSGs that experience moderate mass-lossmore »prior to the explosion. We find gamma-ray peak-luminosities consistent with current upper limits and evaluate that current-generation instruments are able to detect the gamma-rays from Type-IIP explosions at distances up to $\approx 60\,$ kpc and Type-IIn explosions up to $\approx 1.0\,$ Mpc. We also find a good agreement between the thermal X-ray and radio synchrotron emission predicted by our models with a range of observations.

    « less