Buchdahl Parametrization in Gravitation

• Buchdahl parametrization is a method for encoding metric potentials and compactness constraints in gravitational systems, enabling algebraic reduction of Einstein’s equations.
• It offers versatile constructions—from interior stellar models to reciprocal maps and Lie–Hamilton formulations—facilitating explicit solutions for static and higher-curvature spacetimes.
• Applied across GR, f(Q), and R² scenarios, this approach refines predictions for maximum forces, stellar redshifts, and observational deformation parameters.

Buchdahl parametrization denotes a family of constructions in relativistic gravitation that trace back to Buchdahl’s treatment of static metrics, stellar interiors, and higher-curvature vacua. Taken together, the cited works use the term for at least five closely related objects: a compactness potential defined by $g_{tt}=1-2\Phi(R)$ for static spherical spacetimes; prescribed interior metric potentials such as $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$ and its Buchdahl–Vaidya–Tikekar generalization; reciprocal transformations generated by a seed lapse or by a hypersurface-orthogonal Killing vector; a parametrization of the Buchdahl differential equation by integration constants in a Lie–Hamilton formulation; and Buchdahl-inspired variables used to integrate pure $R^2$ vacuum equations and their special $R=0$ descendants (Dadhich, 2022, Sokoliuk et al., 2022, Barrientos et al., 2024, Pereira et al., 2024, Campoamor-Stursberg et al., 2024, Nguyen, 2022).

1. Terminological scope and core meanings

In current usage, “Buchdahl parametrization” is not restricted to a single formula. It refers to a method of encoding either compactness, metric data, or solution families so that the field equations or physical constraints become algebraically or analytically tractable. A recurring feature is that one prescribes a strategically chosen variable—such as a lapse, a radial metric coefficient, or an auxiliary function—so that the remaining structure follows from reduced equations.

Usage	Core formula	Setting
Compactness potential	$g_{tt}=1-2\Phi(R)$	Black holes, Buchdahl stars, maximum-force bounds
Interior metric ansatz	$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$	Static compact stars
BVT generalization	$e^{2\lambda(r)}=\dfrac{1+k r^2/C^2}{1-r^2/C^2}$	Charged stellar interiors
Reciprocal map	$$ds^2=(g_{aa})^{-1}(dx^a)^2+(g_{aa})^{2/(d-3)}g_{ij}dx^i dx^j$$	Vacuum and Einstein–scalar solution generation
Lie–Hamilton parametrization	$(c_1,c_2)$ in the exact solution family	Buchdahl ODE and its deformations
Buchdahl-inspired $R^2$ variables	first-order system for $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$0; $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$1	Pure $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$2 gravity

A common source of confusion is the identification of Buchdahl parametrization with the Buchdahl compactness bound alone. The compactness bound $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$3 is one outcome of a specific perfect-fluid problem, whereas several later papers use “Buchdahl parametrization” for the metric ansatz itself, for reciprocal maps, or for higher-curvature integration variables. The stellar-model literature states this explicitly: in one $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$4 quark-star construction, “Buchdahl parametrization” refers to the chosen interior metric potential rather than to the historical compactness bound (Sokoliuk et al., 2022).

2. Compactness parametrization by $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$5

For static, spherically symmetric spacetimes, one Buchdahl parametrization writes

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$6

so that $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$7. In this form, $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$8 is the gravitational potential for radial motion and directly measures compactness. For a neutral Schwarzschild object, $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$9 in geometric units and $R^2$0 in SI units. For a charged Reissner–Nordström object, the potential is written in the energy-subtraction form

$R^2$1

At the boundary, the paper identifies black holes and Buchdahl stars by the universal values

$R^2$2

In the neutral static case, this reproduces $R^2$3 for a Schwarzschild horizon and $R^2$4 for a perfect-fluid star under Buchdahl’s assumptions, equivalently $R^2$5 in four-dimensional GR (Dadhich, 2022).

This parametrization is used there to formulate maximum-force statements. For two equal-mass, static, uncharged Schwarzschild black holes touching at the horizon,

$R^2$6

The force is determined by differentiating $R^2$7 at the limiting separation. In the charged and rotating case, one introduces

$R^2$8

and for Kerr–Newman versus Schwarzschild obtains

$R^2$9

At the Buchdahl surface, all maximum-force expressions are inherited from the corresponding black-hole formulas through

$R=0$0

The same paper further argues that universality of mass-independent maximum force in pure Lovelock gravity selects the dimensional spectrum $R=0$1, with GR recovered at $R=0$2 (Dadhich, 2022).

A plausible implication is that this compactness version of Buchdahl parametrization serves as a boundary-value language: it encodes the physical distinction between null boundaries and timelike stellar surfaces by a single scalar $R=0$3, while retaining immediate access to redshift, compactness, and limiting interaction strengths.

3. Interior metric ansätze for compact stars

A second major meaning of Buchdahl parametrization appears in interior stellar modeling, where one prescribes the radial metric coefficient and solves for the remaining fields. In Schwarzschild-like coordinates,

$R=0$4

a widely used Buchdahl ansatz is

$R=0$5

For anisotropic stars in GR, one study adopts $R=0$6 or $R=0$7 and shows that the ansatz encompasses the Vaidya–Tikekar and Finch–Skea geometries, with the Durgapal–Bannerji case recovered at $R=0$8. With the anisotropy choice

$R=0$9

the field equations reduce to a hypergeometric equation and split into eight solvable subclasses, four for $g_{tt}=1-2\Phi(R)$0 and four for $g_{tt}=1-2\Phi(R)$1. The resulting solutions are matched to the Schwarzschild exterior through

$g_{tt}=1-2\Phi(R)$2

with mass function

$g_{tt}=1-2\Phi(R)$3

The same analysis reports regularity, energy conditions, causality, $g_{tt}=1-2\Phi(R)$4, and $g_{tt}=1-2\Phi(R)$5 for the calibrated compact-star models (Maurya et al., 2018).

In modified gravity, the same metric philosophy is retained but the admissible parameter range can change. For strange stars in $g_{tt}=1-2\Phi(R)$6 gravity, the chosen Buchdahl metric potential is again

$g_{tt}=1-2\Phi(R)$7

with regularity conditions $g_{tt}=1-2\Phi(R)$8 and $g_{tt}=1-2\Phi(R)$9. The time potential $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$0 is then obtained numerically from the $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$1 field equations plus the MIT bag equation of state $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$2, with boundary matching

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$3

This construction is reported to yield finite central density and pressures, respect energy and causality conditions, and satisfy $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$4 in both linear and nonlinear $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$5 models (Sokoliuk et al., 2022).

A charged extension in $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$6 uses the Buchdahl–Vaidya–Tikekar ansatz

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$7

with $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$8 and the solvability condition $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$9, giving $e^{2\lambda(r)}=\dfrac{1+k r^2/C^2}{1-r^2/C^2}$0. The metric potential becomes

$e^{2\lambda(r)}=\dfrac{1+k r^2/C^2}{1-r^2/C^2}$1

while the charge profile is fixed algebraically. The same framework yields an $e^{2\lambda(r)}=\dfrac{1+k r^2/C^2}{1-r^2/C^2}$2 generalization of the Buchdahl compactness bound,

$e^{2\lambda(r)}=\dfrac{1+k r^2/C^2}{1-r^2/C^2}$3

in the uncharged case, and a charged bound that can exceed $e^{2\lambda(r)}=\dfrac{1+k r^2/C^2}{1-r^2/C^2}$4 (Bhattacharya et al., 2023).

Taken together, these stellar applications show that “Buchdahl parametrization” often means a curvature-controlled or spheroidicity-controlled prescription for $e^{2\lambda(r)}=\dfrac{1+k r^2/C^2}{1-r^2/C^2}$5 rather than a direct statement about $e^{2\lambda(r)}=\dfrac{1+k r^2/C^2}{1-r^2/C^2}$6. This suggests a methodological distinction between boundary compactness parametrizations and interior closure ansätze.

4. Reciprocal solutions and solution-generating maps

A third usage descends from Buchdahl’s reciprocal transformations. For a static seed metric with a hypersurface-orthogonal Killing vector $e^{2\lambda(r)}=\dfrac{1+k r^2/C^2}{1-r^2/C^2}$7,

$e^{2\lambda(r)}=\dfrac{1+k r^2/C^2}{1-r^2/C^2}$8

Buchdahl’s first-kind vacuum theorem states that

$e^{2\lambda(r)}=\dfrac{1+k r^2/C^2}{1-r^2/C^2}$9

is again vacuum for $$ds^2=(g_{aa})^{-1}(dx^a)^2+(g_{aa})^{2/(d-3)}g_{ij}dx^i dx^j$$0. The second-kind transformation generates Einstein–scalar solutions,

$$ds^2=(g_{aa})^{-1}(dx^a)^2+(g_{aa})^{2/(d-3)}g_{ij}dx^i dx^j$$1

with $$ds^2=(g_{aa})^{-1}(dx^a)^2+(g_{aa})^{2/(d-3)}g_{ij}dx^i dx^j$$2. A notable feature is that the transformation does not require $$ds^2=(g_{aa})^{-1}(dx^a)^2+(g_{aa})^{2/(d-3)}g_{ij}dx^i dx^j$$3 to be timelike; spacelike axial or translational Killing vectors are equally admissible (Barrientos et al., 2024).

This spacelike version generates the Schwarzschild–Levi-Civita geometry from Schwarzschild by transforming along $$ds^2=(g_{aa})^{-1}(dx^a)^2+(g_{aa})^{2/(d-3)}g_{ij}dx^i dx^j$$4. In four dimensions with two commuting hypersurface-orthogonal Killing vectors, the same paper defines an infinite family $$ds^2=(g_{aa})^{-1}(dx^a)^2+(g_{aa})^{2/(d-3)}g_{ij}dx^i dx^j$$5 of Buchdahl maps, with special vacuum subfamilies $$ds^2=(g_{aa})^{-1}(dx^a)^2+(g_{aa})^{2/(d-3)}g_{ij}dx^i dx^j$$6 and $$ds^2=(g_{aa})^{-1}(dx^a)^2+(g_{aa})^{2/(d-3)}g_{ij}dx^i dx^j$$7 forming a non-Abelian group. Combined with Kerr–Schild structure, this yields higher-dimensional Levi-Civita extensions of Myers–Perry black holes, together with an algebraically general double copy. In the Einstein–scalar sector, the spacelike second-kind map generates the FJNW–Levi-Civita solution and, after conformal transformation, Levi-Civita extensions of BBMB-type spacetimes. The same analysis emphasizes limitations: Buchdahl-transformed metrics are generally not asymptotically flat, axial singularities are generic, and the 4D Kerr metric cannot be transformed along an axial direction while preserving Kerr–Schild structure (Barrientos et al., 2024).

A closely related extension adds a cosmological constant and scalar–tensor frames. Starting from a static $$ds^2=(g_{aa})^{-1}(dx^a)^2+(g_{aa})^{2/(d-3)}g_{ij}dx^i dx^j$$8-vacuum seed with lapse $$ds^2=(g_{aa})^{-1}(dx^a)^2+(g_{aa})^{2/(d-3)}g_{ij}dx^i dx^j$$9, the Einstein-frame reciprocal map is

$(c_1,c_2)$0

After the Jordan–Einstein conformal transformation $(c_1,c_2)$1, this produces exact Brans–Dicke solutions with quadratic potential $(c_1,c_2)$2. Explicit synchronous solutions are obtained for Schwarzschild–de Sitter, Nariai, and a hyperbolically foliated seed, while the limits $(c_1,c_2)$3 and $(c_1,c_2)$4 recover the original Buchdahl construction and the GR seed, respectively (Pereira et al., 2024).

5. The Buchdahl equation as a Lie–Hamilton parametrized system

A fourth usage appears in the analysis of the Buchdahl equation itself. The generalized nonlinear ODE

$(c_1,c_2)$5

admits the Buchdahl specialization

$(c_1,c_2)$6

which gives

$(c_1,c_2)$7

This is presented as the time-form counterpart of Buchdahl’s static perfect-fluid equation in isotropic coordinates, where the exact solution is

$(c_1,c_2)$8

The Lie–Hamilton reformulation introduces

$(c_1,c_2)$9

and a symplectic form

$R^2$0

with Hamiltonians $R^2$1, $R^2$2, and time-dependent Hamiltonian

$R^2$3

A canonical change of variables linearizes the dynamics to

$R^2$4

so the exact solution family becomes

$R^2$5

where $R^2$6. The constants $R^2$7 are identified as the Lie–Hamilton Buchdahl parametrization (Campoamor-Stursberg et al., 2024).

For the Buchdahl choice $R^2$8, $R^2$9, one has

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$00

which yields

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$01

The mapping to Buchdahl’s original constants is

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$02

recovering

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$03

The same framework admits Poisson–Hopf deformations with parameter $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$04, giving a deformed generalized Buchdahl equation, exact deformed solutions, a small-$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$05 integrable perturbation, and oscillator-algebra extensions in which the higher-dimensional deformed systems become intrinsically coupled rather than sums of copies (Campoamor-Stursberg et al., 2024).

This suggests that, in the ODE literature, Buchdahl parametrization no longer refers primarily to a spacetime coefficient. It refers instead to the explicit parameterization of the exact solution family by canonical invariants and deformation data.

6. Buchdahl-inspired parametrizations in pure $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$06 gravity and astrophysical applications

Buchdahl also appears in pure $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$07 gravity through a gauge-fixed parametrization of the vacuum equations. In static spherical symmetry,

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$08

Buchdahl gauge

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$09

reduces the trace equation $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$10 to

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$11

After a chain of substitutions, the vacuum equations reduce to the generalized Buchdahl ODE

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$12

or equivalently to the first-order system

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$13

The resulting exhaustive metric family is written in compact form as

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$14

with

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$15

The constant-curvature sector $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$16 reproduces Schwarzschild–(A)dS, while $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$17 yields non-constant-curvature vacua that evade the rapid-falloff assumption used in a Lichnerowicz-type no-go theorem (Nguyen, 2022).

A special asymptotically flat $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$18 member introduces the dimensionless Buchdahl parameter

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$19

and reduces to Schwarzschild at $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$20. In Buchdahl-inspired coordinates $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$21, the metric takes the explicit form

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$22

For $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$23, the exterior geometry can be cast into Morris–Thorne form with a throat at

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$24

and the flare-out condition $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$25 is satisfied. The same analysis identifies naked singularities for $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$26 and a non-Schwarzschild borderline case at $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$27 (Nguyen et al., 2023).

Astrophysical work then treats $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$28 as the deformation parameter of a special Buchdahl-inspired spacetime. In units $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$29, the ISCO radius decreases as $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$30 increases; the reported values are $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$31 for $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$32, $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$33 for $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$34, and $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$35 for $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$36. The corresponding radiative efficiencies are $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$37, $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$38, and $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$39. Using a forced epicyclic-resonance identification

$e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$40

the paper reports best-fit intervals $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$41 for XTE J1550–564 and $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$42 for GRO J1655–40, while also noting tension with the shadow bound $e^{\lambda(r)}=\dfrac{K(1+Cr^2)}{K+Cr^2}$43 and attributing moderate accuracy to the neglect of spin (Alloqulov et al., 14 Jul 2025).

A plausible interpretation is that the higher-curvature literature has converted Buchdahl parametrization from a local metric ansatz into a global organizing principle for vacuum solution spaces, wormhole sectors, and observational deformation parameters.