Gravity with higher-curvature terms and second-order field equations: f(R) meets Gauss-Bonnet

arXiv:2510.17965 · gr-qc, astro-ph.HE, hep-th · Submitted 2025-10-20 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Next we'll be talking about the paper "Gravity with higher-curvature terms and second-order field equations: f(R) meets Gauss-Bonnet".

Jocelyn: The paper was written by Fabrizio Corelli, Paolo Pani and Andrea P. Sanna from Department of Physics, Sapienza University of Rome and National Institute for Physics (INFN), Rome Section.

Vera: Stay tuned as we take you through the paper and discuss its implications.

Jocelyn: We also have Subrahmanyan with us today — guest researcher.

Vera: Alright, let's get started.

Summary: Vera: That lead us straight into the summary of what this f(R)-dGB gravity looks like, which is quite different from what we might expect. The authors show that when these two components, f(R) and Gauss-Bonnet, are mixed together in the action, the resulting theory no longer falls into Horndeski’s class.

Jocelyn: And I think that’s a huge deal for me as a pulsar-sky researcher because it means we aren't just looking at a slight modification of what we know; we're looking at something that can accommodate much more complex physics.

Subrahmanyan: The theoretical implication here is profound—it suggests the theory can be recast as a bi-scalar extension involving two nonminimally coupled scalar fields with mutual interactions, which makes it incredibly versatile for our modeling.

Vera: It’s also not just that the structure changes; the paper highlights specific findings, such as how black holes are modified by f(R) terms in this model, which is a change we don've seen in other ways before.

Jocelyn: I'm curious about those modifications to the black hole geometry; it seems like a strong hint that we might be seeing these subtle effects in our own observational data if the coupling constants are right.

Subrahmanyan: The authors also found that, qualitatively, the solutions retain features from EdGB, like a minimum mass and multiple branches, even though the underlying math is radically different.

Improvements & Implications: Vera: Building on that similarity in structure, let’s look at how this theory addresses some of the known problems with gravity. Specifically, they address the issues with well-posedness or hyperbolicity that plague EdGB theory.

Jocelyn: That’s something I need to know because when we try to model these extreme objects from our surveys, having a dynamically stable system is essential for us to trust the results.

Subrahmanyan: The paper suggests a nontrivial mechanism that suppresses the divergence of the Ricci scalar in the black-hole interior, which is a huge structural improvement over just relying on individual higher-order terms.

Vera: It’s interesting that they found this suppression mechanism—it seems to address a major hurdle in understanding how these extreme objects behave right at their core.

Jocelyn: But Subrahmanyan, if the structure is similar to EdGB, does this mean the classical singularity problem is really solved by just adding f(R) terms?

Subrahmanyan: That’s a critical point; the authors argue that even with these additions, at least their nonperturbative level results suggest that ill-posedness isn't resolved by merely fixing the theory with individual higher-order terms.

Vera: They are suggesting that adding isolated high-power corrections might not be enough to fix the fundamental problems in EdGB dynamics.

Jocelyn: That’s a cautionary note for us observers, implying that while f(R) helps modify the picture, the core challenges of classical GR remain deep and require something even more than these single additions.

Technical Deep Dive: Vera: We’ve talked about the qualitative features, but let's look at how this is actually built in math. The authors used a specific form for f(R), like f(R) = R + kappa R n, n in N.

Jocelyn: And they focused on the quadratic (n=two) and the quartic (n=four) cases to see how different powers of curvature affect the solution. I’m curious if these two vastly different coupling regimes lead to similar physical outcomes.

Subrahmanyan: Interestingly, they found that in the large-coupling limit, f(R) corrections tend to approach a universal behavior regardless of the power n, which is a fascinating theoretical finding for me.

Vera: It’s amazing that even though the mathematical structures are distinct, we see this convergence toward a consistent profile. This suggests that maybe nature has some sort of natural tendency toward uniformity in these high-energy regimes.

Jocelyn: For us, this means that whether the coupling is small or massive, our models might converge on certain features of black hole structure if the theory behaves this way at all.

Subrahmanyan: The authors also showed a striking difference in how they handle the Ricci scalar's behavior in the interior compared to EdGB, which is quite telling about where these modifications are actually having their biggest impact.

Conclusion: Vera: We've covered a lot of ground today, from the initial setup of f(R) and Gauss-Bonnet to the technical details of how they affect black hole interiors. It’s clear that "Gravity with higher-curvature terms and second-order field equations: f(R) meets Gauss-Bonnet" offers a powerful new framework.

Jocelyn: I feel like the biggest impact for our field is seeing how these complex theories maintain certain observable features while fundamentally changing the way we approach singularities, providing a much richer set of possibilities for future observations.

Subrahmanyan: And as a final thought on the cosmic picture, I believe this work shows us that even if classical GR breaks down in certain regions, we can find consistent mathematical structures that still maintain recognizable physical patterns.

Vera: It’s certainly a sophisticated piece of work, demonstrating the potential and limitations of high-curvature extensions in finding a viable theory of gravity.

Jocelyn: I'm looking forward to seeing how these ideas translate into actual data from our pulsar surveys next time we discuss this topic.

Fabrizio Corelli, Paolo Pani, Andrea P. Sanna

Department of Physics, Sapienza University of Rome · National Institute for Physics (INFN), Rome Section

gr-qc, astro-ph.HE, hep-th

Submitted: 2025-10-20

Updated: 2026-08-20

Comments: 14 pages, 8 figures, ancillary Mathematica notebook provided as supplemental material

DOI: 10.1103/sv54-2xtn

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 86/100

The gist: The scientific paper "Gravity with higher-curvature terms and second-order field equations: f(R) meets Gauss-Bonnet" explores a unified gravitational framework that combines f(R) gravity and

Key concepts

$f(R)$ and Gauss-Bonnet (dGB)
The authors combine these two components within the theory's action. This combination results in a resulting theory that is fundamentally different from expected physics, specifically showing it does not fall into Horndeski’s class.
Bi-scalar Extension
The theoretical implication of this new model is that it can be recast as a bi-scalar extension. This involves two nonminimally coupled scalar fields with mutual interactions, providing versatility for modeling complex physical systems.
Well-Posedness Improvement
The paper addresses issues of hyperbolicity and well-posedness found in the original dGB theory. A mechanism is identified that suppresses the divergence of the Ricci scalar within the black hole interior.

Terminology

Summary

The scientific paper Gravity with higher-curvature terms and second-order field equations: f(R) meets Gauss-Bonnet explores a unified gravitational framework that combines f(R) gravity and Einstein-dilaton-Gauss-Bonnet (EdGB) gravity, resulting in a theory referred to as f(R) -dGB gravity.

Theoretical Framework and Methodology

The motivation for this study stems from the expectation that General Relativity (GR) must break down in the high-curvature regime. The paper introduces a combination of two well-studied theories: f(R) gravity and EdGB gravity, which are both members of Horndeski’s class, but neither belongs to it when considered individually. The authors demonstrate that combining these terms results in a theory that no longer belongs to Horndeski’s class, which can be recast as a gravitational theory involving two nonminimally coupled scalar fields with nontrivial mutual interactions.

The action for f(R) -dGB gravity is defined as:

S = integral d 4 x sqrt-g [f(R) - grad mu phi grad mu phi + 2 eta(phi)G] / 16 pi G

where f(R) is a generic function of the Ricci scalar R, G is the Gauss-Bonnet invariant (G = R2 − 4R mu nu R mu nu + terms involving R mu nu rho sigma), and eta(phi) is the coupling function. The paper specifically considers an exponential dilatonic coupling, eta(phi) = lambda e-gamma phi, and the family of f(R) functions, f(R) = R + kappa R n, where n represents the power of the curvature correction (e.g., quadratic case n=2 and quartic case n=4).

The resulting field equations are derived from this action, including:

phi = -eta'(phi)G

and

8 pi dGB 1 over G mu nu T mu nu =

Static Black Hole (BH) Solutions and Boundary Conditions

The authors constructed static, asymptotically-flat, spherically-symmetric BH solutions using a shooting method. They employed Painlevé-Gullstrand (PG)-like coordinates, where the line element is given by:

ds squared = -alpha(r) squared dt squared + [dr + alpha(r) zeta(r)dt] + r squared d squared

The solutions are determined by satisfying boundary conditions at both spatial infinity (r to infinity) and at the event horizon (r H).

At r H, regularity imposes specific conditions:

phi(r) phi H0 + phi H1(r - r H) + O(r - r H) squared

chi(r) about chi H0 + chi H1(r - r H) + O(r - r H)

As r to infinity, the solutions must be asymptotically flat, with the total mass M BH defined as:

M BH = r to infinity r squared zeta(r) squared

Key Results and Analysis

The numerical investigations yielded several key findings regarding the behavior of these solutions:

  1. Existence and Structure: The mass-radius diagrams show that, similar to EdGB gravity, there exists a minimum-mass configuration. This configuration divides the solution spectrum into two branches. The branch with larger radii is dynamically stable.

  2. Exterior Profiles: The profiles of the scalar field chi and the dilaton phi were analyzed (Figures 1 and 2). While f(R) corrections modify BH geometries, the qualitative features characteristic of EdGB BHs—such as the existence of a minimum mass and multiple solution branches—are preserved.

  3. Internal Structure and Singularity: Despite being regular at the horizon, all f(R) -dGB solutions encountered a point r S where the integration breaks down. This point corresponds to a zero of the denominator of the right-hand side of Eq. (15). The Kretschmann scalar was found to grow rapidly near r S, following an approximate (r - r S)-1 scaling, which is milder than the r-6 divergence of the Schwarzschild solution.

  4. Suppression Mechanism: A striking difference in the internal structure was observed: in f(R)-dGB theory [the Ricci scalar] is suppressed by about 15 orders of magnitude compared to EdGB theory, and it exhibits an approximately constant trend in the interior and remains regular at r S.

  5. Hyperbolicity: The local well-posedness was examined by computing the principal symbol. While the structure of the field equations is significantly altered—the determinant of the principal symbol for EdGB is proportional to eta t eta r squared, whereas for f(R) -dGB it becomes a fourth-degree polynomial—the qualitative picture remains similar. The analysis showed that the ellipticity and elliptic regions inside the horizon remain similar to those of pure EdGB gravity, suggesting that adding isolated higher-order curvature operators does not, by itself, cure the well-posedness issues inherent to EdGB-type theories at the nonperturbative level.

Conclusions

The authors conclude that f(R) -dGB gravity provides a versatile framework to investigate the UV behavior of gravity. They emphasize that the ill-posedness of EdGB cannot be resolved by adding individual higher-order terms, even at the nonperturbative level, and note that in dynamical settings, elliptic regions will form during the gravitational collapse.

Improvements for AI systems

The following improvements are derived from a meticulous analysis of the provided scientific paper, Gravity with higher-curvature terms and second-order field equations: f(R) meets Gauss-Bonnet.

  1. Modeling Non-Standard Gravitational Framework Complexity:

  • Improvement: The AI will be trained to recognize and solve for physical systems that exist outside established classification schemes (e.g., Horndeski's class). It will specifically learn how to model bi-scalar extensions of gravity, where the interaction between two nonminimally coupled scalar fields (phi and chi) is the fundamental driver of the physics.

  • Capability: The AI can simulate dynamic gravitational systems where standard single-field solutions fail, providing a complete theoretical framework for UV completions that are not simple superpositions of known theories.

  1. High-Precision Numerical Shooting Method Implementation:

  • Improvement: The AI will implement and optimize the complex numerical shooting method required to solve coupled, second-order differential equations (Eq. 13–16 in Section II A). This requires handling highly sensitive boundary conditions (e.g, the quadratic equation for phi H1 in Eq. 27) and maintaining extreme precision (up to 55 significant figures for the quartic case).

  • Capability: The AI can generate robust, high-fidelity numerical solutions for static, asymptotically flat black hole configurations (f(R)-dGB BHs) and can verify the stability of these solutions against numerical divergence.

  1. Advanced Hyperbolicity and Well-Posedness Analysis:

  • Improvement: The AI will be trained to calculate the principal symbol (ij) for multi-field systems with non-linear constraints (Eq. 37) and determine the necessary conditions for strong hyperbolicity (e.g, checking the discriminant in Eq. 45).

  • Capability: The AI can perform rigorous local well-posedness checks on a given spacetime configuration, identifying precisely where a system of dynamical equations breaks down (the radius r E) and distinguishing this breakdown from other physical singularities.

  1. Comparative Theoretical Analysis:
  • Capability: The AI can execute a nuanced comparison between the exterior phenomenology of two distinct theories (EdGB vs. f(R)-dGB). It will differentiate between qualitative preservation (e.g., retaining the minimum mass and multiple branches) and quantitative modification, allowing researchers to understand why f(R) terms appear mild in the exterior despite being powerful at high curvatures.
  1. Internal Structure Diagnostics (Singularity Analysis):

  • Capability: The AI can analyze the internal structure of black hole solutions by comparing various curvature scalars (Kretschmann, Ricci scalar) as a function of radial distance from a singularity (r S). It can specifically identify and quantify:

  • The suppression mechanism for the divergence of the Ricci scalar R within the BH interior.

  • The existence of a curvature singularity (r S) that is milder than classical divergences, providing evidence for potential partial resolution of classical singularities.

  1. Global Dynamical System Prediction:

  • Capability: The AI can predict the long-term evolution and stability of BH systems by integrating the results from the static solutions into a dynamic model, accounting for:

  • The existence of two distinct solution branches (minimum mass vs. minimum radius).

  • The likelihood of forming an elliptic region during gravitational collapse or evaporation, based on the inherent non-hyperbolicity of EdGB-type dynamics.

  1. Parameter Sensitivity and Universal Behavior:
  • Capability: The AI can analyze the relationship between coupling constants (kappa) and solution profiles across different power laws (n=2 vs n=4). It can predict universal behavior in the strong coupling limit, where vastly different theoretical choices (quadratic vs. quartic) lead to nearly identical observable physical outcomes (e.g., overlapping mass-radius diagrams).

Sources

Related papers