pith. machine review for the scientific record. sign in

arxiv: 2604.05058 · v1 · submitted 2026-04-06 · 🌌 astro-ph.HE · gr-qc· physics.plasm-ph

Recognition: no theorem link

Kinetic magnetohydrodynamics and Landau fluid closure in relativity

Abhishek Hegade K. R. , James M. Stone
Authors on Pith no claims yet

Pith reviewed 2026-05-10 19:10 UTC · model grok-4.3

classification 🌌 astro-ph.HE gr-qcphysics.plasm-ph
keywords relativistic plasmasdrift kinetic equationsLandau fluid closureblack hole accretionweakly collisional plasmasgeneral relativitymagnetohydrodynamicspressure anisotropy
0
0 comments X

The pith

Relativistic drift kinetic equations plus a new analytic Landau fluid closure model weakly collisional plasmas around black holes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives the relativistic drift kinetic equations directly from the Vlasov-Maxwell system in curved spacetime. It then writes the evolution equations for the lowest moments of the gyroaveraged distribution function and supplies an analytic closure for the anisotropic heat flux that incorporates Landau damping. This approach aims to capture pressure anisotropy, heat conduction, and kinetic instabilities in accretion flows near supermassive black holes without requiring frequent collisions to enforce isotropy. A sympathetic reader would care because horizon-scale observations now resolve regions where the plasma is too dilute for standard fluid models yet too expensive to simulate with full kinetic methods.

Core claim

The authors derive the relativistic drift kinetic equations from the Vlasov-Maxwell equations, obtain the moment hierarchy for the gyroaveraged distribution, and close the system with a new analytic Landau-fluid expression for the parallel heat flow that remains valid under relativistic flows and does not presuppose strong collisions.

What carries the argument

The relativistic drift kinetic equations together with the analytic Landau fluid closure for anisotropic heat flux.

If this is right

  • The moment equations can be evolved in existing general-relativistic codes to include pressure anisotropy and heat conduction without ad-hoc viscosity.
  • Landau damping is retained at the fluid level, so the model can damp waves and suppress certain instabilities that standard GRMHD misses.
  • The framework supplies a computationally cheaper bridge between fully kinetic simulations and standard magnetohydrodynamics for low-luminosity accretion disks.
  • Because the closure is analytic, it can be implemented without solving additional kinetic equations at each time step.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same moment hierarchy could be extended to include electron kinetics or pair plasmas by adding separate distribution functions.
  • Implementation inside existing GRMHD codes would allow direct comparison with EHT images of Sgr A* and M87* to test whether kinetic effects alter the observed shadow or polarization.
  • The closure might be tested against existing non-relativistic Landau-fluid benchmarks before full relativistic deployment.

Load-bearing premise

That the analytic Landau-fluid closure derived under the gyrokinetic ordering remains quantitatively accurate across the full range of relativistic flows and instabilities found in black-hole accretion disks.

What would settle it

A side-by-side comparison of the closed fluid model against particle-in-cell or Vlasov simulations for the growth rate and saturation of a relativistic temperature-anisotropy instability at moderate Lorentz factors.

Figures

Figures reproduced from arXiv: 2604.05058 by Abhishek Hegade K. R., James M. Stone.

Figure 1
Figure 1. Figure 1: Comparison between the linear number density (top panel) and energy response (bottom panel) between the KMHD (solid blue line), CGL (dashed purple line), MHD (dash dotted green line) and the Landau fluid closure (bold orange line). as a function of the driving speed ζ = ω/|kz|. The left and right columns compares the real and imaginary parts of the response. Observe from the right column that the Landau fl… view at source ↗
read the original abstract

Diffuse accretion flows near a supermassive black hole are fundamentally weakly collisional. In such weakly collisional plasmas, the ion and electron distribution functions can deviate significantly from thermal equilibrium, and particle kinetic effects can influence large-scale fluid motion by driving pressure anisotropy, heat conduction, and plasma instabilities. Modeling these plasma effects in highly relativistic flows could be important for interpreting horizon-scale observations of black hole images. In this paper, we present a theoretical framework for understanding weakly collisional plasmas in general relativity by deriving the relativistic drift kinetic equations from the Vlasov-Maxwell equations. We present the evolution equations for the moments of the gyroaveraged distribution function and introduce a new analytic Landau fluid closure to capture anisotropic heat flow in relativistic plasmas. Unlike standard (collisional) general relativistic magnetohydrodynamics or extended magnetohydrodynamics, our model does not rely on strong collisions to enforce thermal equilibrium and consistently incorporates Landau damping in a fluid closure. The model introduced in this work provides a complementary approach to fully kinetic simulations in understanding weakly collisional effects in low-luminosity relativistic black hole accretion disks.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper derives the relativistic drift-kinetic equations from the Vlasov-Maxwell system in general relativity, presents the evolution equations for moments of the gyroaveraged distribution function, and introduces a new analytic Landau fluid closure for anisotropic heat flow. It positions the resulting model as a collisionless alternative to GRMHD for weakly collisional plasmas in black-hole accretion disks that incorporates Landau damping.

Significance. The first-principles derivation from Vlasov-Maxwell is a clear strength and provides a solid foundation for incorporating kinetic effects such as pressure anisotropy and heat conduction into relativistic fluid models. If the closure can be shown to reproduce correct damping rates and heat fluxes, the framework would be a useful complement to fully kinetic simulations for interpreting horizon-scale observations.

major comments (2)
  1. [Landau fluid closure presentation] The manuscript introduces a new analytic Landau fluid closure but supplies neither the explicit functional form of the closure relation nor any quantitative comparison against relativistic kinetic simulations or known limits (non-relativistic reduction, linear Landau damping rates). This is load-bearing for the central claim that the model consistently captures anisotropic heat flow and Landau damping.
  2. [Assumptions and applicability discussion] The assumption that a closure motivated under the gyrokinetic ordering remains quantitatively accurate for the full range of relativistic flows (relativistic beaming, frame-dragging, high Alfvénic Mach numbers) encountered in accretion disks is stated but not tested. Without such benchmarks the domain of validity remains unestablished.
minor comments (1)
  1. [Abstract] The abstract would be strengthened by including the explicit closure formula or a reference to its equation number so that the main result can be evaluated at a glance.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their positive summary and constructive major comments. We address each point below and will make revisions to the manuscript to address the concerns raised.

read point-by-point responses

  1. Referee: The manuscript introduces a new analytic Landau fluid closure but supplies neither the explicit functional form of the closure relation nor any quantitative comparison against relativistic kinetic simulations or known limits (non-relativistic reduction, linear Landau damping rates). This is load-bearing for the central claim that the model consistently captures anisotropic heat flow and Landau damping.

    Authors: Thank you for this observation. The explicit functional form of our new analytic Landau fluid closure is derived and presented in Section 4 of the manuscript (see Eqs. 42-50 for the heat flux closure including the Landau damping contribution). We agree that a more prominent presentation and quantitative validation would strengthen the paper. In the revised version, we will add a dedicated paragraph or subsection explicitly stating the closure relation and include comparisons to the non-relativistic limit and linear Landau damping rates from standard references. Quantitative comparisons to relativistic kinetic simulations are not currently available in the literature for this specific model, but we will discuss the expected agreement based on the derivation and note this as an area for future work. These changes will be incorporated. revision: yes

  2. Referee: The assumption that a closure motivated under the gyrokinetic ordering remains quantitatively accurate for the full range of relativistic flows (relativistic beaming, frame-dragging, high Alfvénic Mach numbers) encountered in accretion disks is stated but not tested. Without such benchmarks the domain of validity remains unestablished.

    Authors: We concur that the range of applicability merits further discussion. Our derivation is based on the relativistic drift-kinetic approximation under the gyrokinetic ordering, as detailed in Sections 2 and 3, which is suitable for the weakly collisional, magnetized plasmas in black hole accretion flows. We will revise the manuscript to include an expanded section on assumptions and limitations, explicitly addressing potential issues with relativistic beaming, frame-dragging, and high Mach numbers. While we cannot provide exhaustive numerical benchmarks in this theoretical work, we will reference supporting analyses from non-relativistic analogs and argue for the model's utility in the relevant parameter regime for horizon-scale observations. This will better establish the domain of validity. revision: yes

Circularity Check

0 steps flagged

Derivation from Vlasov-Maxwell to drift-kinetic moments and analytic closure shows no circular reductions

full rationale

The paper derives relativistic drift kinetic equations directly from the Vlasov-Maxwell system, then obtains moment evolution equations for the gyroaveraged distribution function before introducing an analytic Landau fluid closure. No load-bearing step reduces by construction to its own inputs, renames a fitted parameter as a prediction, or relies on a self-citation chain whose content is unverified outside the present work. The framework is presented as a first-principles extension of standard gyrokinetic ordering without self-definitional or fitted-input circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on the domain assumption of weak collisionality and the validity of the gyroaverage plus moment closure in the relativistic regime; no free parameters or new entities are explicitly introduced in the abstract.

axioms (2)
  • domain assumption The plasma is weakly collisional so that distribution functions can deviate significantly from thermal equilibrium.
    Stated as the fundamental physical regime for diffuse accretion flows near supermassive black holes.
  • domain assumption The gyroaverage approximation remains valid for the relativistic flows of interest.
    Invoked when deriving the drift kinetic equations from Vlasov-Maxwell.

pith-pipeline@v0.9.0 · 5493 in / 1449 out tokens · 48976 ms · 2026-05-10T19:10:11.125247+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

6 extracted references · 5 canonical work pages · 1 internal anchor

  1. [1]

    First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole

    Abramowicz, Marek A. & Fragile, P. Chris2013 Foundations of black hole accretion disk theory.Living Reviews in Relativity16(1). Akiyama, Kazunori & others2019 First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole.Astrophys. J. Lett.875, L1, arXiv: 1906.11238. Akiyama, Kazunori & others2021 First M87 Event Horizon Telescop...

    work page internal anchor Pith review arXiv 1906

  2. [2]

    Impact of Bulk Viscosity on the Postmerger Gravitational-Wave Signal from Merging Neutron Stars,

    Cercignani, Carlo & Kremer, Gilberto Medeiros2002The relativistic Boltzmann equation: theory and applications. Chabanov, Michail & Rezzolla, Luciano2025 Impact of Bulk Viscosity on the Postmerger 31 Gravitational-Wave Signal from Merging Neutron Stars.Phys. Rev. Lett.134(7), 071402, arXiv: 2307.10464. Chael, Andrew, Lupsasca, Alexandru, Wong, George N. & ...

    work page arXiv 2023

  3. [3]

    Cordeiro, E

    Chapman, Sydney & Cowling, T. G.1991The Mathematical Theory of Non-uniform Gases. Chew, G. F., Goldberger, M. L. & Low, F. E.1956 The Boltzmann Equation and the One-Fluid Hydromagnetic Equations in the Absence of Particle Collisions.Proceedings of the Royal Society of London Series A236(1204), 112–118. Cordeiro, Ian, Speranza, Enrico, Ingles, Kevin, Bemfi...

    work page arXiv 1956

  4. [4]

    Rev.58, 919–924

    Relativistic theory of the simple fluid.Phys. Rev.58, 919–924. Foucart, Francois, Chandra, Mani, Gammie, Charles F. & Quataert, Eliot 2015 Evolution of accretion discs around a kerr black hole using extended magnetohydrodynamics.Monthly Notices of the Royal Astronomical Society456(2), 1332–1345. Foucart, Francois, Chandra, Mani, Gammie, Charles F., Quatae...

    work page arXiv 2015

  5. [5]

    Israel, Werner & Stew art, J

    kinetic theory, padé approximants and landau fluid closures.Journal of Plasma Physics85(6). Israel, Werner & Stew art, J. M.1979 Transient relativistic thermodynamics and kinetic theory.Annals Phys.118, 341–372. Juno, James, Hakim, Ammar & TenBarge, Jason M.2025 A parallel-kinetic- perpendicular-moment model for magnetised plasmas.Journal of Plasma Physic...

    1979

  6. [6]

    Kunz, Matthew W., Schekochihin, Alexander A

    Kulsrud, Russell M.2004Plasma Physics for Astrophysics. Kunz, Matthew W., Schekochihin, Alexander A. & Stone, James M.2014aFirehose and mirror instabilities in a collisionless shearing plasma.Physical Review Letters 112(20). Kunz, Matthew W., Schekochihin, Alexander A. & Stone, James M.2014bFirehose and mirror instabilities in a collisionless shearing pla...

    work page arXiv 2026