Gravitational Metric of a Star

1. Problem Statement

The paper addresses the fundamental problem of determining the space-time metric outside a spatially localized, stationary source of matter (a "star") in General Relativity (GR). While Birkhoff's theorem guarantees that a spherically symmetric vacuum solution is the Schwarzschild metric, real astrophysical objects (like rotating stars) lack perfect spherical symmetry.

• The Challenge: In GR, unlike Newtonian gravity, the gravitational field is non-linear. The metric outside a source is not simply a linear superposition of multipole moments; it involves gravitational self-interactions.
• The Goal: To construct a recursive solution for the metric of a stationary source characterized by an infinite set of mass ($M_L$) and current ($S_L$) multipole moments, expanding the solution in powers of Newton's constant $G$ (Post-Minkowskian expansion) and inverse distance $1/r$ (multipole expansion).
• Specific Objective: To derive the metric up to the second Post-Minkowskian (2PM) order ($O(G^2)$) for arbitrary multipoles, and specifically to demonstrate how this general solution reduces to the Kerr black hole metric in a specific limit, while also describing "Kerr-like" stars that are not black holes.

2. Methodology

The authors employ a recursive perturbative approach based on the classical equations of motion, distinct from the scattering amplitude or worldline effective field theory (EFT) approaches often used in recent literature.

• Gauge Choice: The calculations are performed in the de Donder (harmonic) gauge, defined by $\partial_\nu (\sqrt{-g}g^{\mu\nu}) = 0$.
• Variables: They utilize the gothic metric (Landau-Lifshitz formulation), $g^{\mu\nu} = \eta^{\mu\nu} - h^{\mu\nu}$, which linearizes the Einstein equations into a wave equation with a source term involving the Landau-Lifshitz pseudotensor $\tau^{\mu\nu}_{LL}$.
• Recursive Structure:
  • The metric perturbation is expanded as $h^{\mu\nu} = \sum G^n h^{\mu\nu}_{(n)}$.
  • The equations of motion are solved order-by-order. The $n$-th order term is determined by the convolution of lower-order terms via the Green's function of the Laplacian.
  • Momentum Space: To handle the convolutions efficiently, the authors Fourier transform the problem into momentum space. The recursive step becomes an algebraic manipulation of "currents" $J^{\mu\nu}_{(n)}(k)$.
• Integration Technique:
  • The integrals arising from the convolution of massless propagators are identified as generalized tensor bubble integrals.
  • The authors use Dimensional Regularization to handle divergences, avoiding the more cumbersome Hadamard regularization used in earlier multipole studies.
  • They utilize Symmetric Trace-Free (STF) tensors to handle the angular dependence and multipole structure, which simplifies the algebraic complexity compared to spherical harmonics.

3. Key Contributions and Results

A. General Recursive Formalism

The paper establishes a systematic algorithm to compute the metric to any order in $G$ and any order in the multipole expansion.

• 1PM Order (Linear): The solution is the standard linearized gravity result, where the metric components $h^{(1)}_{00}$ and $h^{(1)}_{0i}$ are directly determined by the mass and current multipole moments via simple $1/r$ expansions.
• 2PM Order (Quadratic): The authors derive the full 2PM correction ($h^{(2)}$). This involves solving the non-linear Einstein equations where the source is the Landau-Lifshitz pseudotensor constructed from the 1PM solution.
  • The result is expressed in terms of generalized bubble integrals in momentum space.
  • Upon Fourier transforming back to position space, the metric is expressed as a series of inverse powers of $r$ ($1/r^2, 1/r^3, \dots$) with coefficients depending on products of multipole moments (e.g., $M^2$, $MS$, $S^2$, $M^2_{ab}$).

B. Explicit 2PM Metric for Truncated Multipoles

The authors provide explicit expressions for the metric components $h_{00}$, $h_{0i}$, and $h_{ij}$ up to 2PM order, truncated to:

• Mass multipoles up to the quadrupole ($M_{ab}$).
• Current multipoles up to the dipole (spin, $S_a$).
• Key Finding: The $h_{ij}$ components, which vanish at 1PM order, are generated at 2PM order through the self-interaction of the lower-order terms. The expressions involve complex tensor structures involving STF projections of the unit vector $\hat{n}$.

C. Recovery of the Kerr Metric

A critical validation step involves checking if the general multipole expansion reduces to the Kerr metric (the solution for a rotating black hole) when the multipoles are constrained to the Kerr relations:
$$M_L + i S_L = M (ia)^L$$
where $a$ is the specific angular momentum.

• Result: The authors demonstrate that their general 2PM solution perfectly matches the Kerr metric in harmonic coordinates up to the calculated order ($O(a^4)$ for $h_{00}$ and $O(a^3)$ for $h_{0i}$).
• Gauge Subtlety: They identify a subtle gauge ambiguity. The "exact" non-perturbative Kerr metric in harmonic gauge contains terms (specifically odd powers of $a$ in $h_{ij}$) that are gauge artifacts. These terms can be removed by a coordinate transformation that preserves the harmonic condition. The authors' recursive method, which builds the solution strictly from the equations of motion without "guessing" the closed form, naturally reproduces the metric without these gauge artifacts. This highlights that the standard closed-form Kerr metric in harmonic gauge includes terms removable by a residual gauge transformation.

D. "Kerr-like" Stars

The paper proposes a physical interpretation: by slightly "tweaking" the multipole moments away from the exact Kerr relations, one can describe compact objects (stars) that are Kerr-like (indistinguishable from a black hole at large distances) but possess no event horizon. This suggests that astrophysical objects could mimic black holes to arbitrary precision without being black holes.

4. Significance and Implications

1. Alternative to Amplitude Methods: The paper challenges the prevailing view that quantum scattering amplitude methods are the only efficient way to solve high-order Post-Minkowskian problems. It demonstrates that a direct, recursive solution of the classical equations of motion, aided by modern integration techniques (bubble integrals, dimensional regularization), is equally powerful and avoids the need for effective field theory matching or point-particle assumptions.
2. Precision Astrophysics: The derived multipole expansions are crucial for future high-precision astrometry (sub-micro-arcsecond accuracy) and gravitational wave modeling. The ability to systematically include higher-order multipoles and self-interaction effects is essential for testing the "No-Hair" theorem and distinguishing between black holes and exotic compact objects.
3. Gauge Understanding: The clarification regarding the gauge redundancy in the harmonic gauge Kerr metric resolves apparent inconsistencies between perturbative expansions and exact solutions, providing a clearer understanding of the coordinate freedom in GR.
4. Computational Framework: The introduction of generalized bubble integrals and the recursive STF tensor formalism provides a robust computational toolkit for extending these calculations to higher orders (3PM and beyond) and non-stationary cases (binary systems).

In summary, this work provides a rigorous, classical, and recursive derivation of the gravitational field of a rotating star, successfully bridging the gap between general multipole expansions and the specific case of the Kerr black hole, while offering a new perspective on the nature of "black hole mimickers."