Viscosity in isotropic cosmological backgrounds in general relativity and Starobinsky gravity

We present a general analysis of the role of shear viscosity in cosmological backgrounds, focusing on isotropic space-time in both Einstein and f(R) gravity. By computing the divergence of the stress-energy tensor in a general class of isotropic (but not necessarily homogeneous) geometries, we show that shear viscosity does not contribute to the background dynamics when the fluid is comoving. This result holds in both the Jordan and Einstein frames, and implies that shear viscosity cannot affect the electromagnetic luminosity distance which is determined by the background light-like geodesics. As an application of our results, we critically examine recent claims that shear viscosity can alter the Hubble evolution and the electromagnetic luminosity distance in Starobinsky gravity. We demonstrate that the continuity equation used in that work is at odds both with the covariant conservation of the stress-energy tensor and the local second law of thermodynamics. We further show that even in models where such modifications could mimic bulk viscosity, the resulting entropy evolution is inconsistent with standard thermodynamic expectations.