⚛️ general relativity

Geodesics of the five-dimensional Kerr-(A)dS black hole in closed form

This paper presents the first general closed-form solution for geodesic motion in five-dimensional Kerr-(A)dS black holes by expressing radial and polar coordinates via manifestly real formulas valid for all root structures, while deriving time and angular coordinates through elliptic integrals of the third kind, alongside corrections to previously published metric and potential expressions.

Original authors: Özgür Sarıoğlu, Enes Yılmaz

Published 2026-10-20

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Özgür Sarıoğlu, Enes Yılmaz

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the vast theater of the universe, gravity is the director, bending the paths of everything that moves. When a massive object like a black hole spins, it drags the very fabric of space and time around with it, creating a swirling current that dictates how light and matter travel. For decades, physicists have studied these paths, known as geodesics, to understand the extreme environments around black holes. While the mathematics for a spinning black hole in our familiar four dimensions of space and time was solved long ago, the story becomes far more complex when we add extra dimensions and a cosmological constant—a term that represents the energy of empty space itself. In five dimensions, a black hole can spin in two independent directions at once, creating a much richer and more tangled gravitational landscape. Understanding exactly how a particle or a beam of light moves through this environment is crucial for testing theories of gravity and interpreting observations of the cosmos, but the equations describing these motions had remained stubbornly difficult to solve in a single, unified way.

A team of researchers at the Middle East Technical University in Turkey has now cracked this problem, providing a complete and exact description of these paths for a five-dimensional spinning black hole. Their work does more than just offer a new solution; it corrects a long-standing error in the mathematical description of the black hole's geometry that had been present in previous literature. By fixing this mistake, they ensured that the equations governing the motion of particles and the behavior of waves near the event horizon are accurate. This correction is vital because the previous error would have led to incorrect predictions about how energy is extracted from the black hole, a process known as superradiance, which is a key mechanism in high-energy astrophysics.

The researchers approached the problem by breaking down the complex motion into simpler, manageable pieces. They showed that the equations describing the path of a particle can be separated into independent parts, one dealing with the distance from the black hole and another with the angle of its orbit. In five dimensions, these equations take the form of quartic polynomials, which are mathematical expressions involving variables raised to the fourth power. This is a significant finding because, in the four-dimensional universe we inhabit, adding a cosmological constant usually makes these equations so complex that they require hyperelliptic functions to solve, which are far more difficult to handle. However, in five dimensions, the complexity remains within the realm of elliptic functions, which are well-understood and manageable. This means that the motion of particles in this five-dimensional setting, whether they are massive objects or massless light rays, can be described using a single, consistent set of mathematical tools.

What makes this achievement particularly powerful is the method used to solve the equations. Instead of solving the problem piece by piece for every possible scenario, the authors found a general solution that works for all cases at once. Previous attempts had to split the problem into different cases depending on the specific roots of the equations, which is like having a different map for every possible terrain. The new solution, however, is a single, unified formula that works regardless of the specific details of the orbit. It allows physicists to start calculating a path from any point along the trajectory, not just from the turning points where the motion reverses direction. This is a crucial advantage because, in certain conditions involving a positive cosmological constant, a particle might never turn around; it could simply spiral outward forever. The old methods would have failed in such a scenario, but the new formula handles it seamlessly.

The team also mapped out how time and the angles of rotation evolve as the particle moves. They discovered that the time coordinate and the two angles of rotation can be calculated using specific types of integrals that depend only on the position of the particle. A striking feature of their solution is how it handles the event horizon, the point of no return. As a particle approaches the horizon, the time coordinate and the angles appear to blow up to infinity in the mathematical description, a known artifact of the coordinate system used. The researchers showed that these infinities are not real physical divergences but are perfectly balanced against each other. When you combine the time and angle in a specific way, the infinities cancel out, leaving a smooth, finite path. This confirms that the particle crosses the horizon without any physical singularity, even though the mathematical coordinates suggest otherwise.

Furthermore, the study clarifies the relationship between the black hole's rotation and the motion of the particle. The solution reveals that the rate at which the particle rotates around the black hole is directly tied to the black hole's own spin at the horizon. This connection is encoded in the residues of the mathematical integrals, which are proportional to a specific combination of the particle's energy and angular momentum. This combination is the same one that determines whether the particle can extract energy from the black hole. The researchers found that if a particle is at the threshold for this energy extraction, the mathematical terms that usually cause the time and angle to diverge vanish entirely, indicating a special, stable interaction.

The paper also addresses what happens in special, simplified cases, such as when the black hole spins at the same rate in both directions or when the cosmological constant is zero. In these scenarios, the complex formulas naturally simplify into more elementary functions, like sines and cosines, without requiring any new assumptions. This consistency across different limits gives the researchers confidence that their general solution is robust and correct. They also noted that this five-dimensional case is unique; in higher dimensions, the mathematical structure changes fundamentally, and the elegant separation of variables used here no longer works. Five dimensions appears to be the last case where this specific, elegant strategy can be applied.

Ultimately, this work provides a complete toolkit for exploring the dynamics of five-dimensional black holes. By correcting the underlying geometry and providing a closed-form solution that covers every possible orbit, the researchers have removed the need for case-by-case calculations. This allows for a more efficient and accurate study of these theoretical objects, which are important for understanding the fundamental nature of gravity and the potential existence of extra dimensions. The solution is not just a mathematical curiosity; it is a practical instrument that can be used to calculate the precise paths of particles and light, predict the timing of signals, and understand the behavior of matter in these extreme, higher-dimensional environments. The work stands as a testament to the power of finding a unified perspective on complex physical problems, turning a tangled web of possibilities into a clear, navigable map.

Technical Summary: Geodesics of the Five-Dimensional Kerr–(A)dS Black Hole in Closed Form

Problem Statement
The paper addresses the problem of finding closed-form analytical solutions for timelike and null geodesics in the five-dimensional Kerr–(A)dS spacetime. While the Hamilton–Jacobi equation for this metric is known to be separable in Boyer–Lindquist coordinates, previous treatments (specifically Grunau, Neumann, and Reimers [9]) integrated the resulting first-order system case-by-case, depending on the real roots of the radial and polar polynomials. Furthermore, the authors identify a significant error in the published inverse metric of the five-dimensional Kerr–(A)dS solution (Kunduri and Lucietti [8]), which propagates into the radial potential and affects the indicial exponents of the radial Klein–Gordon equation at the horizon. The goal is to provide a general solution in closed form that is valid for all root structures, corrects the metric error, and utilizes a unified approach similar to recent work on the four-dimensional Kerr spacetime.

Methodology
The authors employ the following methodological steps:

Metric Correction: They re-derive the inverse metric in Boyer–Lindquist coordinates, identifying and correcting an error in the $g_{\phi_i \phi_j}$ block found in [8]. This correction ensures the radial potential yields the correct superradiance combination ( $\omega - \Omega_i m_i$ ) for the Klein–Gordon equation.
Hamilton–Jacobi Separation: They separate the Hamilton–Jacobi equation, carefully choosing the separation constant to ensure regularity as the cosmological constant $\lambda \to 0$ (the Myers–Perry limit). This involves a specific splitting identity that avoids poles at $\lambda = 0$ .
Dimensionless Reparametrization: The system is reparametrized using the Mino parameter $\zeta$ (where $d\zeta = dl/\rho^2$ ) to decouple the radial and polar motions. The authors then introduce dimensionless variables to reduce the system to two quartic polynomials, $P(x)$ and $Q(u)$ , where $x=r^2$ and $u=\cos^2\theta$ .
Biermann–Weierstrass Integration: Instead of integrating based on the number of real roots (turning points), the authors apply a theorem by Biermann and Weierstrass. This theorem solves $(dy/dz)^2 = f(y)$ for a quartic $f$ about an arbitrary point $y_0$ on the orbit, rather than a zero of $f$ . This allows for a single set of formulas valid regardless of the root structure.
Quadratures: The time coordinate and azimuthal angles are obtained by integrating the resulting expressions, which reduce to elliptic integrals of the third kind. The affine parameter is recovered via an additional quadrature.

Key Contributions and Results

Correction of the Inverse Metric and Potential: The paper corrects the $g_{\phi_i \phi_j}$ block of the inverse metric published in [8]. The authors demonstrate that the uncorrected potential leads to incorrect indicial exponents for the radial Klein–Gordon equation, failing to reproduce the superradiance condition automatically. The corrected potential restores this physical requirement.
Degree Analysis and Elliptic Nature: In dimensionless variables, both the radial and polar orbital polynomials are quartic. The authors highlight a crucial distinction from the four-dimensional Kerr–(A)dS case: in four dimensions, the cosmological constant raises the radial polynomial degree to six (requiring hyperelliptic functions), whereas in five dimensions, the natural variable $x=r^2$ keeps the degree at four. Consequently, both sectors remain elliptic for both timelike and null geodesics in five dimensions, unlike their four-dimensional counterparts.
General Closed-Form Solution: The primary result is a single pair of manifestly real formulas for the radial ( $r$ ) and polar ( $\theta$ ) coordinates, parametrized by an arbitrary point on the orbit ( $x_0, u_0$ ) rather than a turning point. This is significant because for a positive cosmological constant ( $\Lambda > 0$ ), the radial quartic can be positive for all real arguments, meaning the orbit may have no turning points. The classical substitution method fails in this case, but the Biermann–Weierstrass solution remains valid.
Structure of Integrals: The quadratures for the time coordinate ( $T$ ) and the two azimuthal angles ( $\phi_1, \phi_2$ ) reduce exclusively to elliptic integrals of the third kind. The residues of these integrals are shown to be proportional to the superradiance bracket $E - \Omega_i(x_h)\ell_i$ .
Horizon Crossing and Superradiance: The closed-form solution explicitly accounts for coordinate divergences at horizon crossings. The authors show that the divergences in $T$ and $\phi_i$ cancel in the combination $\phi_i - \Omega_i(x_h)T$ , where $\Omega_i(x_h)$ is the horizon's angular velocity. Furthermore, at the superradiant threshold, the residues vanish, and the logarithmic divergences disappear entirely.
Degenerate Cases: The paper provides explicit limits for degenerate cases, including null geodesics ( $\varepsilon=0$ ), the Myers–Perry limit ( $\Lambda=0$ ), and equal rotation parameters ( $\alpha=\beta$ ). In these limits, the quartic polynomials degenerate to cubics or quadratics, and the solutions reduce to elementary functions or lower-order elliptic integrals, all covered by the general formalism.

Significance and Claims
The authors claim that this work provides the first general, closed-form analytical solution for geodesics in the five-dimensional Kerr–(A)dS spacetime that does not require case-splitting based on the root structure of the orbital polynomials. The solution is valid for all causal characters (timelike and null) and all values of the cosmological constant, including the de Sitter case where orbits may lack turning points.

The paper emphasizes that the five-dimensional case is structurally unique: it is the highest dimension where the geodesic equations can be decoupled by a single reparametrization (the Mino parameter) and where the orbital polynomials remain quartic, keeping the problem within the realm of elliptic functions. For dimensions $D > 5$ , the problem would require Abelian functions of multiple arguments.

The correction of the inverse metric is presented as a necessary step for any future calculations involving quasinormal modes or superradiant amplification in this background, as the previous error shifted the physical boundary conditions. The unified solution allows for the direct calculation of observables such as periastron shifts and Lense–Thirring precession without the need for numerical integration or case-by-case analysis.

More like this