The magnetorotational instability (MRI) plays a crucial role in regulating the accretion efficiency in astrophysical accretion discs. In low-luminosity discs around black holes, such as Sgr A* and M87, Coulomb collisions are infrequent, making the MRI physics effectively collisionless. The collisionless MRI gives rise to kinetic plasma effects that can potentially affect its dynamic and thermodynamic properties. We present 2D and 3D particle-in-cell (PIC) plasma simulations of the collisionless MRI in stratified discs using shearing boxes with net vertical field. We use pair plasmas, with initial β = 100 and concentrate on subrelativistic plasma temperatures (kBT ≲ mc2). Our 2D and 3D runs show disc expansion, particle and magnetic field outflows, and a dynamo-like process. They also produce magnetic pressure dominated discs with (Maxwell stress dominated) viscosity parameter α ∼ 0.5–1. By the end of the simulations, the dynamo-like magnetic field tends to dominate the magnetic energy and the viscosity in the discs. Our 2D and 3D runs produce fairly similar results, and are also consistent with previous 3D MHD (magnetohydrodynamic) simulations. Our simulations also show non-thermal particle acceleration, approximately characterized by power-law tails with temperature-dependent spectral indices − p. For temperatures $k_\mathrm{ B}T \sim 0.05-0.3\, mc^2$, we find p ≈ 2.2–1.9. The maximum accelerated particle energy depends on the scale separation between MHD and Larmor-scale plasma phenomena in a way consistent with previous PIC results of magnetic reconnection-driven acceleration. Our study constitutes a first step towards modelling from first principles potentially observable stratified MRI effects in low-luminosity accretion discs around black holes.
- Award ID(s):
- 1654168
- NSF-PAR ID:
- 10220048
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- ISSN:
- 0035-8711
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in two dimensions. We solve our island kinetic equation numerically for three different types of initial distribution: Dirac delta, Gaussian and power-law distributions. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde {t}^{-1}$ , the number of islands decreases as $\tilde {t}^{-1}$ and the averaged area of islands grows as $\tilde {t}$ , where $\tilde {t}$ is the time normalised to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. Finally, we discuss the underlying connection of our island-merger models to the (self-similar) decay of magnetohydrodynamic turbulence.more » « less
-
null (Ed.)We report on an analytical and numerical study of the dynamics of a three-dimensional array of identical magnetic flux tubes in the reduced-magnetohydrodynamic description of the plasma. We propose that the long-time evolution of this system is dictated by flux-tube mergers, and that such mergers are dynamically constrained by the conservation of the pertinent (ideal) invariants, viz. the magnetic potential and axial fluxes of each tube. We also propose that in the direction perpendicular to the merging plane, flux tubes evolve in a critically balanced fashion. These notions allow us to construct an analytical model for how quantities such as the magnetic energy and the energy-containing scale evolve as functions of time. Of particular importance is the conclusion that, like its two-dimensional counterpart, this system exhibits an inverse transfer of magnetic energy that terminates only at the system scale. We perform direct numerical simulations that confirm these predictions and reveal other interesting aspects of the evolution of the system. We find, for example, that the early time evolution is characterized by a sharp decay of the initial magnetic energy, which we attribute to the ubiquitous formation of current sheets. We also show that a quantitatively similar inverse transfer of magnetic energy is observed when the initial condition is a random, small-scale magnetic seed field.more » « less
-
Abstract Particle acceleration during magnetic reconnection is a long-standing topic in space, solar, and astrophysical plasmas. Recent 3D particle-in-cell simulations of magnetic reconnection show that particles can leave flux ropes due to 3D field-line chaos, allowing particles to access additional acceleration sites, gain more energy through Fermi acceleration, and develop a power-law energy distribution. This 3D effect does not exist in traditional 2D simulations, where particles are artificially confined to magnetic islands due to their restricted motions across field lines. Full 3D simulations, however, are prohibitively expensive for most studies. Here, we attempt to reproduce 3D results in 2D simulations by introducing ad hoc pitch-angle scattering to a small fraction of the particles. We show that scattered particles are able to transport out of 2D islands and achieve more efficient Fermi acceleration, leading to a significant increase of energetic particle flux. We also study how the scattering frequency influences the nonthermal particle spectra. This study helps achieve a complete picture of particle acceleration in magnetic reconnection.more » « less
-
more » « less
ABSTRACT Black hole and neutron star environments often comprise collisionless plasmas immersed in strong magnetic fields and intense baths of low-frequency radiation. In such conditions, relativistic magnetic reconnection can tap the magnetic field energy, accelerating high-energy particles that rapidly cool by inverse Compton (IC) scattering the dense photon background. At the highest particle energies reached in bright gamma-ray sources, IC scattering can stray into the Klein–Nishina regime. Here, the Comptonized photons exceed pair-production threshold with the radiation background and may thus return their energy to the reconnecting plasma as fresh electron–positron pairs. To reliably characterize observable signatures of such Klein–Nishina reconnection, in this work, we present first-principles particle-in-cell simulations of pair-plasma relativistic reconnection coupled to Klein–Nishina and pair-production physics. The simulations show substantial differences between the observable signatures of Klein–Nishina reconnection and reconnection coupled only to low-energy Thomson IC cooling (without pair production). The latter regime exhibits strong harder-when-brighter behaviour; the former involves a stable spectral shape independent of overall brightness. This spectral stability is reminiscent of flat-spectrum radio quasar (FSRQ) GeV high states, furnishing evidence that Klein–Nishina radiative physics operates in FSRQs. The simulated Klein–Nishina reconnection pair yield spans from low to order-unity and follows an exponential scaling law in a single governing parameter. Pushing this parameter beyond its range studied here might give way to a copious pair-creation regime. Besides FSRQs, we discuss potential applications to accreting black hole X-ray binaries, the M87* magnetosphere, and gamma-ray binaries.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/staa3849
