A mathematical study of massive Dirac fields around rotating black holes reports that a stationary, source-free condensate cannot remain nonzero on the horizon under the conditions it examines. The decisive step is a flux argument: for a stationary, bounded field with no source, the nonnegative current through the horizon must vanish pointwise. In the paper's formulation, that requires both the total flux and the horizon field to be zero, ruling out a non-vanishing horizon condensate.
The result is a statement about the separated equations used in the model. The paper develops a geometric and kinematic framework for stationary bound states of massive Dirac fields on Kerr and Kerr-Newman backgrounds, then tests the boundary constraints that such states would have to meet. The manuscript is an arXiv preprint, version 1, dated 26 August 2026.
The horizon condition decides the outcome
To examine those constraints, the authors map the Kerr-Dirac system to a globally integrated radial problem with an inhomogeneous source term. They use a Kinnersley null-tetrad reduction to decouple the equations, then apply Frobenius analysis, a local power-series test of how solutions behave near the horizon. The same workflow is used to study radial phase winding and what happens as the black hole approaches extremal spin.
That phase analysis matters because, when the winding parameter k is nonzero, the Boyer-Lindquist radial phase winds indefinitely as the tortoise coordinate r* moves toward minus infinity at the horizon. The reported cycle count is approximately Nwind = |omega - m OmegaH|L/(2pi), so the winding grows with the horizon-frequency detuning and the interval used to measure it. At synchronization, defined here by omega = m OmegaH, that detuning disappears.
The key point comes at synchronization. The local Frobenius structure is regular on both branches and contains no logarithmic solution. The horizon condition nevertheless requires psi at the horizon to equal zero. That sets the leading series coefficient a0 to zero and trivializes the series. Within the separated framework, the paper therefore excludes a sourceless, stationary bound state that corotates with the horizon.
The sourced calculation gives a related result. The source amplitude is proportional to the invariant k and to |omega - m OmegaH|/kappa+, and it falls to zero at synchronization. In the reported Kerr scan, the bound-state frequency decreases monotonically as spin rises, crosses the synchronization line transversally rather than locking to it, and later tends toward zero before the branch terminates at a critical spin.
A filled, flattened geometry
The angular calculation predicts a different shape for the rotating condensate. Because the azimuthal quantum number for the modeled fermion mode is half-integer, the density does not fall to zero on the rotation axis. The resulting geometry is globally filled and oblate rather than hollow and toroidal. It is a geometric output of the modeled angular structure.
The angular spectrum was also checked independently with a matrix spectral calculation. The two angular methods agreed to order 10^-13, a solver-consistency check rather than a statistical confidence interval.
Checks beyond the basic Kerr case
Several numerical checks anchor the broader scan. In the spherical, zero-spin benchmark with m = 1, high-precision shooting found OmegaA = 0.409128952 for Case A and OmegaB = 0.719735134 for Case B. The roots remained stable within 10^-10 when the numerical setup was changed. The modeled cases were then examined across frequency-spin space and the corotation regime.
The extension to charged rotating holes uses a shifted synchronization rule: omega = m OmegaH + q PhiH. In the paper's Kerr-Newman model, this is the generalized zero-source condition and is identified with the charged-superradiance threshold. The authors report that black-hole charge elevates the global synchronization locus.
The approach to extremal spin adds a warning about limits. As the spin parameter a approaches M from below, the synchronized frequency tends to m/(2M). At exactly a = M, however, the radial pole structure becomes an irregular singular point, so the standard Frobenius expansion used away from extremality is no longer valid. The limiting frequency is therefore continuous in the reported calculation, while the radial singularity structure is not uniform at the endpoint.
Taken together, the calculations frame synchronization as a veto on a nonzero, source-free stationary Dirac condensate at the horizon within the separated Kerr equations. The regular local branches do not produce a nonzero state once the flux boundary condition is imposed.
Paper data and sources
Original title: Stationary Dirac condensates around Kerr black holes
Authors: Sen Guo, Peng-Yu Chen, Yi-Han Huang et al.
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-08-26
DOI: Not available
Original paper · Full text