An analytic toy model for relativistic accretion in Kerr spacetime

We present a relativistic model for the stationary axisymmetric accretion flow of a rotating cloud of non-interacting particles falling onto a Kerr black hole. Based on a ballistic approximation, streamlines are described analytically in terms of timelike geodesics, while a simple numerical scheme is introduced for calculating the density field. A novel approach is presented for describing all of the possible types of orbit by means of a single analytic expression. This model is a useful tool for highlighting purely relativistic signatures in the accretion flow dynamics coming from a strong gravitational field with frame-dragging. In particular, we explore the coupling due to this between the spin of the black hole and the angular momentum of the infalling matter. Moreover, we demonstrate how this analytic solution may be used for benchmarking general relativistic numerical hydrodynamics codes by comparing it against results of smoothed particle hydrodynamics simulations for a collapsar-like setup. These simulations are performed first for a ballistic flow (with zero pressure) and then for a hydrodynamical one where we measure the effects of pressure gradients on the infall, thus exploring the extent of applicability of the ballistic approximation.

💡 Research Summary

The paper presents an analytic “toy” model for stationary, axisymmetric accretion of a rotating cloud of non‑interacting particles onto a Kerr black hole, using a ballistic approximation in which pressure, magnetic fields, self‑gravity and radiative effects are neglected. Under this assumption the fluid is represented by a continuum of test particles that follow timelike geodesics of the Kerr spacetime. The authors first write down the four conserved quantities associated with the Kerr metric – the specific energy E, the axial component of the specific angular momentum ℓ, the Carter constant Q, and the normalization condition u·u = −1 – and express the four‑velocity components (u^r, u^θ, u^φ, u^t) in terms of the radial and polar potentials R(r) and Θ(θ).

A major technical contribution is the derivation of a single analytic expression for the full trajectory (r(τ), θ(τ), φ(τ), t(τ)) using Jacobi elliptic functions. By exploiting standard identities, the authors combine the separate radial and latitudinal integrals into one compact formula, which greatly simplifies the mapping from the initial injection sphere (r₀, θ₀) to any point (r, θ) along a streamline. This mapping is essential because the conserved quantities are functions of the initial polar angle θ₀, so each streamline has its own set of (E, ℓ, Q).

The velocity field is presented in two frames. First, the coordinate velocities in Boyer‑Lindquist coordinates are given directly from the geodesic equations. Second, the authors transform to locally non‑rotating frames (LNRFs) to obtain physically measurable three‑velocities V̄^r, V̄^θ, V̄^φ and the associated Lorentz factor γ. This dual description allows one to compute both the kinematic structure of the flow and the observable quantities relevant for astrophysical modeling.

To obtain the density distribution, the authors use the steady‑state continuity equation. They impose two regularity conditions: (i) streamlines must not have turning points in r or θ before reaching the equatorial plane, and (ii) the mapping θ₀ → θ must remain non‑singular (∂θ/∂θ₀ > 0). Under these constraints a simple numerical scheme integrates the Jacobian of the transformation from (r₀, θ₀) to (r, θ) to compute the particle number density n(r, θ).

The model is applied to several boundary‑condition setups that are relevant for collapsar scenarios, where a massive star’s core collapses to a rotating black hole and a low‑angular‑momentum envelope feeds an accretion disc. The authors explore how the black‑hole spin a couples to the injected angular momentum ℓ, producing frame‑dragging‑induced asymmetries in the streamlines and shifting the radius at which particles intersect the equatorial plane.

A key part of the paper is the comparison with three‑dimensional smoothed particle hydrodynamics (SPH) simulations performed with a modified version of the public Gadget‑2 code. Two series of simulations are carried out: (1) a purely ballistic run where pressure forces are switched off, allowing a direct quantitative test of the analytic geodesic solution; (2) a hydrodynamic run that includes a simple neutrino‑cooling prescription and a second‑order pseudo‑Newtonian potential approximating Kerr gravity. In the ballistic case the SPH particle trajectories, density profiles, and accretion rates match the analytic predictions to within numerical noise, confirming that the model can serve as a stringent benchmark for relativistic hydrodynamics codes (time‑stepping, pseudo‑potential implementation, etc.). In the hydrodynamic case, pressure gradients cause streamlines to deviate from the geodesic paths near the disc, forming a caustic surface and demonstrating the limits of the ballistic approximation.

The authors argue that, despite its simplifications, the analytic model is valuable for several reasons. It isolates pure relativistic effects (strong gravity and frame dragging) from fluid microphysics, providing clear physical insight. It offers a fast, inexpensive way to explore a wide parameter space (spin, angular‑momentum distribution, accretion rate) before committing to expensive full‑scale GRMHD simulations. Finally, it supplies a well‑defined test problem for code verification, especially for codes that aim to incorporate general relativistic effects via pseudo‑potentials or approximate metric treatments.

The paper concludes by acknowledging the model’s limitations (absence of pressure, magnetic fields, self‑gravity, and time‑dependence) and suggesting future extensions, such as incorporating magnetohydrodynamic forces, radiative transfer, or evolving the injection surface to study transient accretion episodes. Nonetheless, the presented analytic solution stands as a robust, transparent benchmark for relativistic accretion studies in Kerr spacetime.