pith. sign in

arxiv: 2512.06908 · v2 · submitted 2025-12-07 · 🌀 gr-qc · math.AP

A first-order formulation of f(R) gravity in spherical symmetry

Philippe G. LeFloch , Filipe C. Mena This is my paper

Pith reviewed 2026-05-17 00:22 UTC · model grok-4.3

classification 🌀 gr-qc math.AP
keywords f(R) gravityspherical symmetryfirst-order formulationnonlocal equationscharacteristic initial value problemBondi-Sachs coordinatesscalar curvature
0
0 comments X

The pith

f(R) gravity equations in spherical symmetry reduce to a first-order system of two coupled nonlocal hyperbolic equations for the scalar field and curvature.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops an augmented first-order formulation of f(R) gravity by treating the spacetime scalar curvature as an independent unknown paired with the scalar field. This converts the original higher-derivative equations into a closed integro-differential system of two coupled, nonlocal, nonlinear hyperbolic equations. Working in generalized Bondi-Sachs coordinates, initial data is prescribed on an asymptotically flat future light cone with vertex at the center, and the minimal regularity conditions at the center are identified to ensure the reduced system matches the full f(R) theory. The approach eliminates higher derivatives from the principal part while preserving geometric structure and extends prior methods for the Einstein-scalar-field system.

Core claim

By treating the spacetime scalar curvature as an independent unknown, we obtain a closed first-order nonlocal system for the pair (phi,R). Extending Christodoulou's method, we recast the f(R) field equations as an integro-differential system of two coupled, first-order, nonlocal, nonlinear hyperbolic equations, whose principal unknowns are the scalar field and the spacetime scalar curvature.

What carries the argument

The augmented characteristic first-order formulation that promotes the spacetime scalar curvature to an independent unknown, yielding a closed integro-differential nonlocal hyperbolic system for the pair consisting of the scalar field and curvature.

If this is right

  • The characteristic initial value problem can be posed with data on an asymptotically flat future light cone.
  • Structural properties of global solutions for the massive scalar field follow directly from the reduced system.
  • The higher-derivative character of the equations is removed at the level of the principal part.
  • Equivalence to the original theory holds exactly when the identified minimal regularity conditions at the center are met.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Numerical evolution codes for spherical f(R) models could be simplified by evolving only first-order equations.
  • Similar augmentations might extend to other higher-curvature gravity theories in spherical symmetry.
  • Global existence or singularity formation results could be adapted from known techniques for hyperbolic systems.

Load-bearing premise

The minimal regularity conditions at the center are precisely those ensuring equivalence between the reduced first-order system and the full f(R) equations.

What would settle it

A concrete solution to the first-order nonlocal system that satisfies the stated regularity at the center yet fails to solve the original f(R) field equations, or vice versa.

read the original abstract

We develop an augmented characteristic, first-order formulation of the field equations in f(R) gravity governing the global evolution of a (possibly) massive scalar field phi under spherical symmetry. This formulation is designed to isolate the genuine dynamical degrees of freedom while preserving the geometric structure of the theory. By treating the spacetime scalar curvature as an independent unknown, we obtain a closed first-order nonlocal system for the pair (phi,R). This augmentation eliminates the higher-derivative character of the original equations at the level of the principal part. Our formulation allows us to pose the characteristic initial value problem and to establish several structural properties of solutions. More precisely, we work in generalized Bondi-Sachs coordinates and prescribe initial data on an asymptotically flat, future light cone with vertex at the center of symmetry, and we identify the minimal regularity conditions required at the center. These regularity conditions are shown to be precisely those ensuring equivalence between the reduced system and the full f(R) equations. Extending Christodoulou's method for the Einstein-scalar-field system, we recast the f(R) field equations as an integro-differential system of two coupled, first-order, nonlocal, nonlinear hyperbolic equations, whose principal unknowns are the scalar field and the spacetime scalar curvature.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript develops an augmented characteristic first-order formulation of the f(R) field equations in spherical symmetry for a (possibly massive) scalar field phi. Treating the spacetime scalar curvature R as an independent unknown yields a closed integro-differential system of two coupled first-order nonlocal nonlinear hyperbolic equations for the pair (phi, R) in generalized Bondi-Sachs coordinates. Initial data are prescribed on an asymptotically flat future light cone with vertex at the center; the authors identify minimal regularity conditions at the center (vanishing of odd powers in the radial expansions) and claim these conditions ensure equivalence to the original fourth-order f(R) equations. The construction extends Christodoulou's method for the Einstein-scalar-field system.

Significance. If the claimed equivalence holds, the formulation isolates the dynamical degrees of freedom, removes higher derivatives from the principal part, and permits a well-posed characteristic initial-value problem. This could facilitate global-existence results and numerical studies of f(R) gravity while preserving geometric structure, analogous to the utility of first-order reductions in the Einstein-scalar-field system.

major comments (1)
  1. [Section on regularity conditions at the center and equivalence (near the statement of the minimal regularity conditions)] The central equivalence claim—that the minimal regularity conditions at the center (vanishing of odd powers in the radial expansions of phi and R) are precisely those that make every solution of the reduced integro-differential system satisfy the original fourth-order f(R) equations—requires explicit verification. Because the system contains integral terms along the light cone arising from the spherical reduction, a mismatch in leading-order behavior at r=0 could propagate non-locally and violate the trace equation or Bianchi identities used to close the system. The manuscript asserts sufficiency of the conditions but does not supply the detailed leading-order cancellation check or regularity analysis of the integral kernels that would confirm this.
minor comments (2)
  1. [Abstract] The abstract states that 'several structural properties of solutions' are established but does not enumerate them; a short list would improve readability.
  2. [Derivation of the integro-differential system] Notation for the integral kernels and the precise form of the nonlocal terms should be introduced with a dedicated equation number early in the derivation section to aid cross-reference.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and for recognizing the potential utility of this first-order formulation for global existence and numerical studies. We address the single major comment below and will strengthen the manuscript accordingly.

read point-by-point responses

  1. Referee: [Section on regularity conditions at the center and equivalence (near the statement of the minimal regularity conditions)] The central equivalence claim—that the minimal regularity conditions at the center (vanishing of odd powers in the radial expansions of phi and R) are precisely those that make every solution of the reduced integro-differential system satisfy the original fourth-order f(R) equations—requires explicit verification. Because the system contains integral terms along the light cone arising from the spherical reduction, a mismatch in leading-order behavior at r=0 could propagate non-locally and violate the trace equation or Bianchi identities used to close the system. The manuscript asserts sufficiency of the conditions but does not supply the detailed leading-order cancellation check or regularity analysis of the integral kernels that would confirm this.

    Authors: We agree that an explicit leading-order verification is required to rigorously confirm that the minimal regularity conditions (vanishing of odd powers) close the system without introducing inconsistencies via the nonlocal integrals. In the revised manuscript we will insert a dedicated subsection performing the radial expansion analysis near r=0. This will explicitly compute the leading terms for phi and R, verify cancellation in the integral kernels along the light cone, and confirm consistency with the trace equation and Bianchi identities, thereby establishing that the reduced system is equivalent to the original fourth-order equations precisely under these conditions. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation of first-order integro-differential system

full rationale

The paper introduces the spacetime scalar curvature R as an auxiliary variable to recast the fourth-order f(R) equations into a first-order nonlocal hyperbolic system for the pair (phi, R). This reduction is performed directly from the original field equations in generalized Bondi-Sachs coordinates, extending Christodoulou's method without any step that defines a quantity in terms of itself or renames a fitted parameter as a prediction. The minimal regularity conditions at the center (vanishing of odd powers in radial expansions) are stated as those ensuring equivalence to the full equations, but this equivalence is asserted as a consequence of the construction rather than presupposed by definition. No load-bearing self-citation chain or ansatz smuggling is present in the provided derivation outline; the central claim retains independent mathematical content in isolating the dynamical degrees of freedom while preserving geometric structure.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The work rests on standard assumptions of general relativity and hyperbolic PDE theory plus the spherical-symmetry reduction; the main addition is the auxiliary variable R whose independent status is postulated to close the system.

axioms (2)
  • domain assumption Spherical symmetry of the spacetime and fields
    Reduces the problem to an effectively 1+1 dimensional setting in generalized Bondi-Sachs coordinates.
  • domain assumption Asymptotically flat initial data on a future light cone
    Prescribed on an asymptotically flat future light cone with vertex at the center of symmetry.
invented entities (1)
  • Independent spacetime scalar curvature R no independent evidence
    purpose: To eliminate the higher-derivative character of the original f(R) equations by closing the system at first order
    Treated as an additional unknown alongside phi to obtain a first-order nonlocal system

pith-pipeline@v0.9.0 · 5521 in / 1385 out tokens · 38811 ms · 2026-05-17T00:22:23.077566+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages

  1. [1]

    A. Alho, S. Carloni, and C. Uggla , On dynamical systems approaches and methods in f(R) cosmology, J. Cosmo. Astroparticle Phys. 08 (2016), 064

    work page 2016

  2. [2]

    Avalos, I

    R. Avalos, I. P. Lobo, T. Sanomiya, and C. Romero , On the Cauchy problem for Weyl-geometric scalar-tensor theories of gravity, J. Math. Phys. 59 (2018), 062502

    work page 2018

  3. [3]

    Bernard, L

    L. Bernard, L. Lehner, and R. Luna , Challenges to global solutions in Horndeski's theory, Phys. Rev. D 100 (2019), 024011

    work page 2019

  4. [4]

    Brady, C.M

    P.R. Brady, C.M. Chambers, and S.M.C.V. Gon c alves , Phases of massive scalar field collapse, Phys. Rev. D 56 (1997), R6057--R6061

    work page 1997

  5. [5]

    Brady, K

    S.E. Brady, K. Clough, P. Figueras, and A.D. Kovács, Inflation dynamics in higher-derivative scalar-tensor theories of gravity, Preprint arXiv:2505.17986 (2025)

    work page arXiv 2025

  6. [6]

    Modelling self-consistently beyond General Relativ- ity,

    R. Cayuso, P. Figueras, T. Fran ca, and L. Lehner , Modeling self-consistently beyond general relativity, Preprint arXiv:2303.07246 (2023)

    work page arXiv 2023

  7. [7]

    Chae , Global existence of spherically symmetric solutions to the coupled Einstein and nonlinear Klein-Gordon system, Class

    D. Chae , Global existence of spherically symmetric solutions to the coupled Einstein and nonlinear Klein-Gordon system, Class. Quant. Grav. 18 (2001), 4589-- 4605

    work page 2001

  8. [8]

    Chowdhury, K

    S. Chowdhury, K. Pal, K. Pal, and T. Sarkar , Collapse in f(R) gravity and the method of R matching, Eur. Phys. J. C 80 (2020), 902

    work page 2020

  9. [9]

    Christodoulou , The problem of a self-gravitating scalar field, Commun

    D. Christodoulou , The problem of a self-gravitating scalar field, Commun. Math. Phys. 105 (1986) 337--361

    work page 1986

  10. [10]

    Christodoulou , Global existence of generalized solutions of the spherically symmetric Einstein-scalar equations in the large, Commun

    D. Christodoulou , Global existence of generalized solutions of the spherically symmetric Einstein-scalar equations in the large, Commun. Math. Phys. 106 (1986), 587--621

    work page 1986

  11. [11]

    Christodoulou , A mathematical theory of gravitational collapse, Commun

    D. Christodoulou , A mathematical theory of gravitational collapse, Commun. Math. Phys. 109 (1987), 613--647

    work page 1987

  12. [12]

    Christodolou , Bounded variation solutions of the spherically symmetric Einstein--scalar field equations, Comm

    D. Christodolou , Bounded variation solutions of the spherically symmetric Einstein--scalar field equations, Comm. Pure Applied Math. 46 (1993), 1131--1220

    work page 1993

  13. [13]

    Christodolou , Examples of naked singularity formation in the gravitational collapse of a scalar field, Ann

    D. Christodolou , Examples of naked singularity formation in the gravitational collapse of a scalar field, Ann. Math. 140 (1994), 607--653

    work page 1994

  14. [14]

    W. J. Cocke and J.M. Cohen , Cauchy problem in the scalar-tensor gravitational theory, J. Math. Phys. 9 (1968), 971

    work page 1968

  15. [15]

    Costa and F

    J.L. Costa and F. C. Mena , The spherically symmetric Einstein-scalar field system with a positive cosmological constant, J. Hyp. Diff. Eqs. 18 (2021), 311--341

    work page 2021

  16. [16]

    Costa, R

    J.L. Costa, R. L. Duarte, and F. C. Mena , Improved decay estimates and C^2 -asymptotic stability of solutions to the Einstein-scalar field system in spherical symmetry, J. Hyp. Diff. Eqs. 20 (2023), 853--902

    work page 2023

  17. [17]

    Delsate, D

    T. Delsate, D. Hilditch, and H. Witek , Initial value formulation of dynamical Chern-Simons gravity, Phys. Rev. D 91 (2015), 024027

    work page 2015

  18. [18]

    De Felice and S.Tsujikawa , f(R) theories, Liv

    A. De Felice and S.Tsujikawa , f(R) theories, Liv. Rev. Rel. 13: 3 (2010)

    work page 2010

  19. [19]

    Ferreira , Cosmological tests of gravity, Ann

    P. Ferreira , Cosmological tests of gravity, Ann. Rev. Astron. Astrophys. 57 (2019), 335--374

    work page 2019

  20. [20]

    Figueras, A

    P. Figueras, A. Held, and A.D. Kov\'acs, Well-posed initial value formulation of general effective field theories of gravity, Preprint arXiv:2407.08775

    work page arXiv

  21. [21]

    Giannakopoulos, N

    T. Giannakopoulos, N. T. Bishop, D. Hilditch, D. Pollney, and M. Zilhao , Gauge structure of the Einstein field equations in Bondi-like coordinates, Phys. Rev. D 105 (2022), 084055

    work page 2022

  22. [22]

    Goldwirth and T

    D.S. Goldwirth and T. Piran , Gravitational collapse of massless scalar field and cosmic censorship, Phys. Rev. D 36 (1987), 3575

    work page 1987

  23. [23]

    Goncalves and I.G

    S.M.C.V. Goncalves and I.G. Moss , Black hole formation from massive scalar fields, Class. Quant. Grav. 14 (1997), 2607--2615

    work page 1997

  24. [24]

    Guo and P

    J.-Q. Guo and P. S. Joshi , Spherical collapse for the Starobinsky R^2 model, Phys. Rev. D 94 (2016), 044063

    work page 2016

  25. [25]

    LeFloch and Y

    P.G. LeFloch and Y. Ma , The mathematical validity of the f(R) theory of modified gravity, M\'em. French Math. Soc. 150 (2017), 1--119

    work page 2017

  26. [26]

    LeFloch and Y

    P.G. LeFloch and Y. Ma , The global nonlinear stability of Minkowski space for self-gravitating massive fields, World Scientific Press, Singapore, 2017

    work page 2017

  27. [27]

    LeFloch and Y

    P.G. LeFloch and Y. Ma , The Euclidean-hyperboloidal foliation method: application to f(R) modified gravity, Gen. Relativity Gravitation 56 (2024), 63 pp

    work page 2024

  28. [28]

    P. G. LeFloch, F. C. Mena and T-C. Nguyen , Spherically symmetric evolution of self-gravitating massive fields, J. Differential Equations 394 (2024), 31--97

    work page 2024

  29. [29]

    Magnano and L.M

    G. Magnano and L.M. Sokolowski , On physical equivalence between modified gravity theories and a general relativistic self-gravitating scalar field, Phys. Rev. D 50 (1994), 5039--5059

    work page 1994

  30. [30]

    Michel and I

    F. Michel and I. G. Moss , Relativistic collapse of axion stars, Phys. Lett. B 785 (2018), 9--13

    work page 2018

  31. [31]

    Papallo and H

    G. Papallo and H. S. Reall , On the local well-posedness of Lovelock and Horndeski theories, Phys. Rev. D 96 (2017), 044019

    work page 2017

  32. [32]

    H. S. Reall, N. Tanahashi, and B. Way , Causality and hyperbolicity of Lovelock theories, Class. Quant. Grav. 31 (2014), 205005

    work page 2014

  33. [33]

    J. L. Ripley and F. Pretorius , Hyperbolicity in spherical gravitational collapse in a Horndeski theory, Phys. Rev. D 99 (2019), 084014

    work page 2019

  34. [34]

    Salgado , The Cauchy problem of scalar-tensor theories of gravity, Class

    M. Salgado , The Cauchy problem of scalar-tensor theories of gravity, Class. Quantum Grav. 23 (2006), 4719

    work page 2006

  35. [35]

    Zhang and H

    X. Zhang and H. Lu , Critical behavior in a massless scalar field collapse with self-interaction potential, Phys. Rev. D 91 (2015), 044046

    work page 2015

  36. [36]

    Zhang, Z-Y

    C-Y. Zhang, Z-Y. Tang, and B. Wang , Collapsing spherical stars in f(R) gravity, Phys. Rev. D 94 (2016), 104013

    work page 2016