Gauge-Invariant Gravitational Wave Polarization in Metric f(R) Gravity with Cosmological Implications

We develop a fully gauge-invariant analysis of gravitational-wave polarizations in metric f(R) gravity with a particular focus on the modified Starobinsky model f(R)=R+αR2−2Λ, whose constant-curvature solution Rd=4Λ provides a natural de Sitter background for both early- and late-time cosmology. Linearizing the field equations around this background, we derive the Klein–Gordon equation for the curvature perturbation δR and show that the scalar propagating mode acquires a mass mψ2=1/(6α), highlighting how the same scalar degree of freedom governs inflationary dynamics at high curvature and the propagation of gravitational waves in the current accelerating Universe. Using the scalar–vector–tensor decomposition and a decomposition of the perturbed Ricci tensor, we obtain a set of fully gauge-invariant propagation equations that isolate the contributions of the scalar, vector, and tensor modes in the presence of matter. We find that the tensor sector retains the two transverse–traceless polarizations of General Relativity, while the scalar sector contains an additional massive scalar propagating degree of freedom, which manifests through breathing and longitudinal tidal responses depending on the wave regime and detector frame. Through the geodesic deviation equation—computed both in a local Minkowski patch and in fully covariant de Sitter form—we independently recover the same polarization content and identify its tidal signatures. The resulting framework connects the extra scalar polarization to cosmological observables: the massive scalar propagating mode sets the range of the fifth force, influences the time evolution of gravitational potentials, and affects the propagation and dispersion of gravitational waves on cosmological scales. This provides a unified, gauge-invariant link between gravitational-wave phenomenology and the cosmological implications of metric f(R) gravity.