⚛️ general relativity

Where the Gauss--Bonnet coupling enters the Bondi--Sachs hierarchy

This paper integrates the Bondi–Sachs hierarchy for scalar–Gauss–Bonnet gravity with an arbitrary coupling function to demonstrate that the coupling does not alter the Bondi mass-loss or angular-momentum evolution equations, but instead first manifests dynamically in the $O(r^{-3})$ angular metric coefficient and the $\Psi_0$ Weyl scalar evolution through a cross-term between the scalar charge aspect and gravitational news.

Original authors: Maickol Muñoz-Palma, Alvaro Restuccia, Francisco Tello-Ortiz

Published 2026-10-07

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Maickol Muñoz-Palma, Alvaro Restuccia, Francisco Tello-Ortiz

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

Gravity is the force that holds the universe together, but it is also the force that ripples through it when massive objects collide. For decades, physicists have used a precise mathematical framework, known as the Bondi-Sachs hierarchy, to describe how these ripples behave as they travel across the cosmos and reach the edge of the observable universe. This framework acts like a set of rules for how the fabric of spacetime fades away as it moves further from its source, allowing scientists to calculate exactly how much energy is lost when black holes merge or neutron stars crash. These calculations are crucial because they provide the baseline for interpreting the signals detected by gravitational wave observatories on Earth. If the rules of gravity were slightly different than Einstein predicted, the way these ripples fade and the way they carry information about mass and spin would change, potentially revealing new physics hidden in the data.

A specific theory of gravity, which emerges from the study of string theory, suggests that spacetime might not just be curved by mass, but also influenced by a subtle interaction between a scalar field and a geometric quantity called the Gauss-Bonnet invariant. This interaction is governed by a coupling constant, a number that determines how strongly this new effect influences the universe. The big question for physicists was exactly where this new influence would appear in the mathematical description of gravitational waves. Would it change the fundamental way energy is lost from a system? Would it alter the way the universe remembers the passage of a gravitational wave? Or would it hide itself so deep in the mathematical structure that it would only show up in the most minute details?

Researchers Maickol Muñoz-Palma, Alvaro Restuccia, and Francisco Tello-Ortiz set out to answer this by performing a direct, step-by-step integration of the equations that govern this theory. They did not rely on approximations or simulations that might miss subtle effects; instead, they followed the equations through the entire hierarchy of the Bondi-Sachs framework, tracking exactly how the new coupling term behaves at every level of detail. Their work reveals a surprising and elegant result: the new coupling does not disrupt the most important features of gravitational radiation. The formulas that describe how much mass is lost by a system and how its angular momentum changes remain exactly the same as in standard Einstein gravity, regardless of the strength of this new coupling. This means that the fundamental balance laws governing the universe's energy and momentum are robust, even in the presence of this exotic interaction.

The researchers found that the new coupling is effectively invisible at the highest levels of the hierarchy, where the most significant physical quantities are defined. It only begins to make its presence known much deeper in the mathematical structure, specifically in the evolution of a coefficient that describes the shape of the wave at a very high order of precision. When it finally appears, it does so through a specific interaction between the scalar field and the changing shape of the gravitational wave, acting as a source term in an equation that governs the leading behavior of a particular component of the wave's geometry. This delayed appearance is due to a geometric cancellation: the dominant part of the new interaction is aligned in such a way that it is annihilated by the primary direction of the radiation, pushing its effect to a much later stage in the calculation.

Furthermore, the study confirms that this theory does not introduce mathematical inconsistencies that would spoil the smoothness of the solution. In some complex theories, new terms can force the appearance of logarithmic functions, which would break the standard pattern of how the fields fade away. The authors demonstrated that for this specific theory, the expansion remains smooth and free of such logarithmic terms, regardless of the specific function used to describe the coupling. This ensures that the mathematical description of the universe remains well-behaved and predictable.

The implications of these findings are significant for the interpretation of gravitational wave data. Because the mass-loss and angular-momentum formulas are unchanged, the standard methods used to extract information from gravitational wave signals remain valid even if this scalar-Gauss-Bonnet gravity is the true description of nature. The new physics only enters the picture at a level of detail that corresponds to the very fine structure of the wave, specifically affecting the evolution of the quadrupole moment of the spacetime geometry. For black holes in a vacuum, where the scalar field is generated solely by the coupling itself, the effect is even smaller, appearing only at the second order of the coupling strength. This means that for the most common astrophysical events, the deviations from Einstein's theory are incredibly subtle, hiding in the higher-order coefficients of the wave rather than altering its primary characteristics. The work provides a clear map of where to look for these deviations and confirms that the core principles of gravitational radiation are far more resilient than one might expect.

Technical Summary: Where the Gauss–Bonnet Coupling Enters the Bondi–Sachs Hierarchy

Problem Statement
The asymptotic structure of radiating spacetimes at future null infinity ( $\mathscr{I}^+$ ), formulated within the Bondi–Sachs framework, underpins the modern understanding of gravitational radiation, including the Bondi mass-loss formula, the BMS group, and flux-balance laws. While these structures are well-established in General Relativity (GR) and various scalar-tensor theories, their behavior in theories with higher-derivative couplings remains less understood. Specifically, Scalar–Gauss–Bonnet (sGB) gravity, a prototype for theories exhibiting hairy or spontaneously scalarized black holes, presents a unique challenge. Its Horndeski representation involves logarithmic functionals ( $G_5 \propto f'(\Phi) \ln X$ ), which were previously excluded from classifications of theories admitting smooth, log-free asymptotic expansions. Consequently, it remained an open question whether the Gauss–Bonnet coupling $\lambda$ modifies the leading asymptotic data (mass and angular momentum), whether it introduces logarithmic terms into the expansion, and exactly at which order in the $1/r$ hierarchy the coupling first influences the dynamical evolution of the metric and scalar fields.

Methodology
The authors address these questions by directly integrating the Bondi–Sachs hierarchy for sGB gravity with an arbitrary coupling function $f(\Phi)$ and a radiating massless scalar field. The methodology involves:

Covariant Analysis: Working directly with the covariant field equations to compute the asymptotic expansion of the effective Gauss–Bonnet stress tensor ( $T^{GB}_{ab}$ ) on a Bondi–Sachs background.
Null Tetrad Counting: Utilizing a null tetrad adapted to Bondi coordinates to analyze the algebraic structure of the Gauss–Bonnet tensor. The authors exploit the fact that the dominant Hessian component of the scalar field contracts with the principal null direction of the radiative curvature, effectively annihilating the highest-order Weyl scalars ( $\Psi_4, \Psi_3$ ).
Order-by-Order Integration: Following the derived falloffs through the Bondi–Sachs hierarchy, which includes hypersurface equations (determining metric functions $\beta, U^A, V$ ) and evolution/supplementary equations (determining the mass aspect $M$ , angular momentum aspect $N$ , and shear evolution).
Axisymmetric Gauge: Performing explicit integration in the axisymmetric, twist-free gauge to verify the general covariant results and isolate specific terms, while maintaining arbitrary dependence on the retarded time $u$ and angular coordinates.

Key Contributions and Results

Asymptotic Falloff of the Gauss–Bonnet Stress Tensor:
The effective Gauss–Bonnet stress tensor decays significantly faster than the slowest falloff compatible with a smooth $\mathscr{I}^+$ . Every component decays between one and three powers of $1/r$ faster than the bounds required for a smooth expansion. The leading coefficients depend solely on the leading data (shear $c_{AB}$ , mass aspect $M$ , scalar coefficient $\phi_1$ ) and the derivative of the coupling function $f'(\phi_0)$ . Crucially, the trace-free angular component $H^{TF}_{AB}$ is found to be proportional to the cross-term $-f'(\phi_0)\phi_1 \partial_u N_{AB}$ , representing the interaction between the scalar charge aspect and the time derivative of the news tensor.
Preservation of Leading Balance Laws:
The Bondi mass-loss formula and the evolution equation of the angular-momentum aspect remain exactly those of the minimally coupled Einstein–scalar system. This result holds non-perturbatively in $\lambda$ under the assumption of a smooth $1/r$ expansion. The coupling $\lambda$ does not appear in the supplementary equations at orders $O(r^{-2})$ and $O(r^{-3})$ , which govern the conservation of mass and angular momentum.
Absence of Logarithmic Terms:
Despite the Horndeski representation of sGB gravity involving logarithmic functionals, the coupling generates no logarithmic terms in the metric expansion through the orders considered. The condition that removes logarithms (specifically the vanishing of the trace-free part of the $O(r^{-2})$ metric coefficient) is preserved for arbitrary $f(\Phi)$ . The logarithmic terms inherent in the Horndeski map cancel among the functionals when treated in covariant form.
First Entry of the Coupling:
The Gauss–Bonnet coupling first enters a dynamical equation at the evolution of the $O(r^{-3})$ coefficient of the conformal angular metric (denoted as $C$ in the axisymmetric gauge or $e^{TF}_{AB}$ generally).
- In the metric hierarchy, the new term is $-2\lambda f'(\phi_0) \phi_1 \partial_u N_{AB}$ .
- In the Newman–Penrose formalism, the evolution equation for the leading coefficient of $\Psi_0$ acquires a source term: $-12\lambda f'(\phi_0) \phi_1 \bar{\Psi}^0_4$ .
- The Gauss–Bonnet invariant itself decays as $r^{-6}$ and affects the scalar hierarchy only from the evolution of the $r^{-5}$ coefficient ( $\phi_5$ ) onward.
Special Cases and Vacuum Black Holes:
For vacuum black-hole spacetimes where the scalar is sourced solely by the coupling (e.g., during inspiral and merger), the scalar charge $\phi_1$ is of order $\lambda$ . Consequently, all Gauss–Bonnet contributions to the hierarchy become of order $\lambda^2$ . If the linear coupling vanishes ( $f'(\phi_0) = 0$ ), as in theories leading to spontaneous scalarization, the leading effects are pushed one further order in $1/r$ , appearing at $O(r^{-3})$ with a coefficient proportional to $f''(\phi_0)\phi_1^2$ .

Significance
This work fills a critical gap in the understanding of higher-curvature corrections to gravitational radiation. By providing the explicit Bondi–Sachs solution, the paper:

Validates and Extends Previous Claims: It confirms via direct integration that the Gauss–Bonnet term does not modify the memory effect or the supertranslation balance law, extending these conclusions to the angular-momentum sector.
Resolves the Logarithmic Ambiguity: It demonstrates that sGB gravity, despite involving logarithmic functionals in its Horndeski form, satisfies all falloff bounds for a smooth $\mathscr{I}^+$ without generating logarithmic metric terms, provided the expansion is smooth.
Identifies the Hierarchy Level: It precisely identifies the level of the asymptotic hierarchy (the spin-2 charge associated with $\Psi_0$ ) where higher-curvature couplings of this type first leave an imprint. This generalizes previous results from the $D \to 4$ Einstein–Gauss–Bonnet limit to a theory with scalar radiation.

The authors conclude that the structural reason for this "late" appearance is the contraction of the scalar Hessian with the principal null direction of the radiative curvature, which suppresses the coupling's influence on the leading asymptotic structure. The analysis is limited to the axisymmetric gauge for explicit integration and assumes a smooth expansion, leaving higher-order polyhomogeneous terms and non-axisymmetric evolution equations for future investigation.

Technical Summary: Where the Gauss–Bonnet Coupling Enters the Bondi–Sachs Hierarchy

More like this