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

Abstract: We develop a first-order formulation of the field equations in f(R) gravity governing the global evolution of a (possibly massive) scalar field under spherical symmetry. Our formulation allows us to pose the characteristic initial value problem and to establish several properties of solutions. More precisely, we work in generalized Bondi-Sachs coordinates and prescribe initial data on an asymptotically Euclidean, future light cone with vertex at the center of symmetry, and we identify the precise regularity conditions required at the center. Following and extending Christodoulou's approach to the Einstein-massless 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 identify the regularity conditions at the center of symmetry and impose natural assumptions on the scalar-field potential and on the function f(R) governing the gravitational Lagrangian density. As an application, we prove the monotonicity of the Hawking mass in this setting and formally analyze the singular limit in which the integrand f(R) of the action approaches R, corresponding to the Einstein-Hilbert action. Hence, the formulation isolates the essential evolution and constraint content on the future domain of dependence of two null hypersurfaces and is designed to facilitate subsequent advances in geometric analysis and robust numerical simulations of spherical collapse in modified gravity.