skip to main content

Title: Landau Modes are Eigenmodes of Stellar Systems in the Limit of Zero Collisions

We consider the spectrum of eigenmodes in a stellar system dominated by gravitational forces in the limit of zero collisions. We show analytically and numerically using the Lenard–Bernstein collision operator that the Landau modes, which are not true eigenmodes in a strictly collisionless system (except for the Jeans unstable mode), become part of the true eigenmode spectrum in the limit of zero collisions. Under these conditions, the continuous spectrum of true eigenmodes in a collisionless system, also known as the Case–van Kampen modes, is eliminated. Furthermore, because the background distribution function in a weakly collisional system can exhibit significant deviations from a Maxwellian distribution function over long times, we show that the spectrum of Landau modes can change drastically even in the presence of slight deviations from a Maxwellian, primarily through the appearance of weakly damped modes that may be otherwise heavily damped for a Maxwellian distribution. Our results provide important insights for developing statistical theories to describe thermal fluctuations in a stellar system, which are currently a subject of great interest forN-body simulations as well as observations of gravitational systems.

more » « less
Award ID(s):
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 271
["Article No. 271"]
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions but unadulterated by possible broadening from turbulence. This distribution function is essentially non-Maxwellian and varies with the heliospheric distance. The stability analysis is performed with the Vlasov–Maxwell linear solver leopard. We find that depending on the heliospheric distance, the core-strahl electron distribution becomes unstable with respect to sunward-propagating kinetic-Alfvén, magnetosonic, and whistler modes, in a broad range of propagation angles. The wavenumbers of the unstable modes are close to the ion inertial scales, and the radial distances at which the instabilities first appear are on the order of 1 au. However, we have not detected any instabilities driven by resonant wave interactions with the superthermal strahl electrons. Instead, the observed instabilities are triggered by a relative drift between the electron and ion cores necessary to maintain zero electric current in the solar wind frame (ion frame). Contrary to strahl distributions modelled by shifted Maxwellians, the electron strahl obtained as a solution of the kinetic equation is stable. Our results are consistent with the previous studies based on a more restricted solution for the electron strahl. 
    more » « less
  2. Abstract

    Most stellar evolution models predict that black holes (BHs) should not exist above approximately 50–70M, the lower limit of the pair-instability mass gap. However, recent LIGO/Virgo detections indicate the existence of BHs with masses at and above this threshold. We suggest that massive BHs, including intermediate-mass BHs (IMBHs), can form in galactic nuclei through collisions between stellar-mass BHs and the surrounding main-sequence stars. Considering dynamical processes such as collisions, mass segregation, and relaxation, we find that this channel can be quite efficient, forming IMBHs as massive as 104M. This upper limit assumes that (1) the BHs accrete a substantial fraction of the stellar mass captured during each collision and (2) that the rate at which new stars are introduced into the region near the SMBH is high enough to offset depletion by stellar disruptions and star–star collisions. We discuss deviations from these key assumptions in the text. Our results suggest that BHs in the pair-instability mass gap and IMBHs may be ubiquitous in galactic centers. This formation channel has implications for observations. Collisions between stars and BHs can produce electromagnetic signatures, for example, from X-ray binaries and tidal disruption events. Additionally, formed through this channel, both BHs in the mass gap and IMBHs can merge with the SMBHs at the center of a galactic nucleus through gravitational waves. These gravitational-wave events are extreme- and intermediate-mass ratio inspirals.

    more » « less
  3. We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modelling the turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non-Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the fluid model of Chew–Goldberger–Low (CGL) in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the fourth-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the ‘normal’ closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the ‘normal’ closure, and beyond the fourth-order moment, Landau fluid closures are required. 
    more » « less
  4. We report analytical and numerical investigations of subion-scale turbulence in low-beta plasmas using a rigorous reduced kinetic model. We show that efficient electron heating occurs and is primarily due to Landau damping of kinetic Alfvén waves, as opposed to Ohmic dissipation. This collisionless damping is facilitated by the local weakening of advective nonlinearities and the ensuing unimpeded phase mixing near intermittent current sheets, where free energy concentrates. The linearly damped energy of electromagnetic fluctuations at each scale explains the steepening of their energy spectrum with respect to a fluid model where such damping is excluded (i.e., a model that imposes an isothermal electron closure). The use of a Hermite polynomial representation to express the velocity-space dependence of the electron distribution function enables us to obtain an analytical, lowest-order solution for the Hermite moments of the distribution, which is borne out by numerical simulations. 
    more » « less
  5. null (Ed.)
    We investigate kinetic entropy-based measures of the non-Maxwellianity of distribution functions in plasmas, i.e. entropy-based measures of the departure of a local distribution function from an associated Maxwellian distribution function with the same density, bulk flow and temperature as the local distribution. First, we consider a form previously employed by Kaufmann & Paterson ( J. Geophys. Res. , vol. 114, 2009, A00D04), assessing its properties and deriving equivalent forms. To provide a quantitative understanding of it, we derive analytical expressions for three common non-Maxwellian plasma distribution functions. We show that there are undesirable features of this non-Maxwellianity measure including that it can diverge in various physical limits and elucidate the reason for the divergence. We then introduce a new kinetic entropy-based non-Maxwellianity measure based on the velocity-space kinetic entropy density, which has a meaningful physical interpretation and does not diverge. We use collisionless particle-in-cell simulations of two-dimensional anti-parallel magnetic reconnection to assess the kinetic entropy-based non-Maxwellianity measures. We show that regions of non-zero non-Maxwellianity are linked to kinetic processes occurring during magnetic reconnection. We also show the simulated non-Maxwellianity agrees reasonably well with predictions for distributions resembling those calculated analytically. These results can be important for applications, as non-Maxwellianity can be used to identify regions of kinetic-scale physics or increased dissipation in plasmas. 
    more » « less