Theoretical astrophysics and plasma physics at RPC

Abstract:

The probability density function of accretion disc luminosity fluctuations at high observed energies (i.e. energies larger than the peak temperature scale of the disc) is derived, under the assumption that the temperature fluctuations are lognormally distributed. Thin disc theory is used throughout. While lognormal temperature fluctuations would imply that the disc’s bolometric luminosity is also lognormal, the observed Wien-like luminosity behaves very differently. For example, in contrast to a lognormal distribution, the standard deviation of the derived distribution is not linearly proportional to its mean. This means that these systems do not follow a linear rms-flux relationship. Instead they exhibit very high intrinsic variance, and undergo what amounts to a phase transition, in which the mode of the distribution (in the statistical sense) ceases to exist, even for physically reasonable values of the underlying temperature variance. The moments of this distribution are derived using asymptotic expansion techniques. A result that is important for interpreting observations is that the theory predicts that the fractional variability of these disc systems should increase as the observed frequency is increased. The derived distribution will be of practical utility in quantitatively understanding the variability of disc systems observed at energies above their peak temperature scale, including X-ray observations of tidal disruption events.

Abstract:

The following minor errors have been found in the published article. (Equation presented). These errors do not affect any of the figures, results, or conclusions of the paper.

Abstract:

Abstract The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.