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. In deriving this reduced two-equation system, we impose natural assumptions on the scalar-field potential and on the function f(R) governing the gravitational Lagrangian density. These hypotheses correspond to standard viability and positivity conditions commonly imposed in the f(R) literature.
We prove several equivalence and monotonicity properties, including for the Hawking mass in this setting.The proposed formulation separates the essential null evolution from the radial constraint reconstruction (via explicit integral relations) on the future domain of dependence of the initial light cone. This structure makes the system amenable to characteristic energy estimates and to stable numerical implementation in spherical symmetry.
Read more