Theoretical astrophysics and plasma physics at RPC

Abstract:

Ion-gyroradius-scale microinstabilities typically have a frequency comparable to the ion transit frequency. Due to the small electron-to-ion mass ratio and the large electron transit frequency, it is conventionally assumed that passing electrons respond adiabatically in ion-gyroradius-scale modes. However, in gyrokinetic simulations of ion-gyroradius-scale modes in axisymmetric toroidal magnetic fields, the nonadiabatic response of passing electrons can drive the mode, and generate fluctuations in narrow radial layers, which may have consequences for turbulent transport in a variety of circumstances. In flux tube simulations, in the ballooning representation, these instabilities reveal themselves as modes with extended tails. The small electron-to-ion mass ratio limit of linear gyrokinetics for electrostatic instabilities is presented, in axisymmetric toroidal magnetic geometry, including the nonadiabatic response of passing electrons and associated narrow radial layers. This theory reveals the existence of ion-gyroradius-scale modes driven solely by the nonadiabatic passing electron response, and recovers the usual ion-gyroradius-scale modes driven by the response of ions and trapped electrons, where the nonadiabatic response of passing electrons is small. The collisionless and collisional limits of the theory are considered, demonstrating parallels in structure and physical processes to neoclassical transport theory. By examining initial-value simulations of the fastest-growing eigenmodes, the predictions for mass-ratio scaling are tested and verified numerically for a range of collision frequencies. Insight from the small electron-to-ion mass ratio theory may lead to a computationally efficient treatment of extended modes.

Abstract:

Current theoretical models predict a mass gap with a dearth of stellar black holes (BHs) between roughly 50 M_⊙ and 100 M_⊙, while above the range accessible through massive star evolution, intermediate-mass BHs (IMBHs) still remain elusive. Repeated mergers of binary BHs, detectable via gravitational-wave emission with the current LIGO/Virgo/Kagra interferometers and future detectors such as LISA or the Einstein Telescope, can form both mass-gap BHs and IMBHs. Here we explore the possibility that mass-gap BHs and IMBHs are born as a result of successive BH mergers in dense star clusters. In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products after they receive significant recoil kicks due to anisotropic emission of gravitational radiation. Using for the first time simulations that include full stellar evolution, we show that a massive stellar BH seed can easily grow to ∼10³–10⁴ M_⊙ as a result of repeated mergers with other smaller BHs. We find that lowering the cluster metallicity leads to larger final BH masses. We also show that the growing BH spin tends to decrease in magnitude with the number of mergers so that a negative correlation exists between the final mass and spin of the resulting IMBHs. Assumptions about the birth spins of stellar BHs affect our results significantly, with low birth spins leading to the production of a larger population of massive BHs.

Abstract:

In self-gravitating N-body systems, small perturbations introduced at the start, or infinitesimal errors that are produced by the numerical integrator or are due to limited precision in the computer, grow exponentially with time. For Newton's gravity, we confirm earlier results that for relatively homogeneous systems, this rate of growth per crossing time increases with N up to N 7sim; 30, but that for larger systems, the growth rate has a weaker scaling with N. For concentrated systems, however, the rate of exponential growth continues to scale with N. In relativistic self-gravitating systems, the rate of growth is almost independent of N. This effect, however, is only noticeable when the system's mean velocity approaches the speed of light to within three orders of magnitude. The chaotic behavior of systems with more than a dozen bodies for the usually adopted approximation of only solving the pairwise interactions in the Einstein-Infeld-Hoffmann equation of motion is qualitatively different than when the interaction terms (or cross terms) are taken into account. This result provides a strong motivation for follow-up studies on the microscopic effect of general relativity on orbital chaos, and on the influence of higher-order cross-terms in the Taylor-series expansion of the Einstein-Infeld-Hoffmann equations of motion.