We consider the scattering of gravitational waves off a Schwarzschild black hole in $f(R)$ gravity. We show that the reflection and transmission coefficients for tensor waves are the same as in General Relativity. While the scalar waves, which are not present in General Relativity, demonstrate interesting features. The equation that governs these scalar waves can be reduced to a Volterra integral equation. Analysis of this equation shows that a larger fraction of these waves are reflected compared to what one obtains for tensors. This may provide a novel observational signature for fourth order gravity.

We generalise the Buchdahl-Bondi limit for the case of static, spherically symmetric, relativistic compact stars immersed in Schwarzschild vacuum in f(R)-theory of gravity, subject to very generic regularity, thermodynamic stability and matching conditions. Similar to the case of general relativity, our result is model independent and remains true for any physically realistic equation of state of standard stellar matter. We show that an extra-massive stable star can exist in these theories, with surface redshift larger than 2, which is forbidden in general relativity. This result gives a novel and interesting observational test for validity or otherwise of general relativity and also provides a possible solution to the dark matter problem.

We prove a Jebsen-Birkhoff like theorem for f(R) theories of gravity in order to to find the necessary conditions required for the existence of the Schwarzschild solution in these theories and demonstrate that the rigidity of such solutions of f(R) gravity is valid even in the perturbed scenario.

The problem of matching different regions of spacetime in order to construct inhomogeneous cosmological models is investigated in the context of Lagrangian theories of gravity constructed from general analytic functions f(R), and from non-analytic theories with f(R)=R^n. In all of the cases studied, we find that it is impossible to satisfy the required junction conditions without the large-scale behaviour reducing to that expected from Einstein's equations with a cosmological constant. For theories with analytic f(R) this suggests that the usual treatment of weak-field systems as perturbations about Minkowski space may not be compatible with late-time acceleration driven by anything other than a constant term of the form f(0), which acts like a cosmological constant. In the absence of Minkowski space as a suitable background for weak-field systems, one must then choose and justify some other solution to perform perturbative analyses around. For theories with f(R)=R^n we find that no known spherically symmetric vacuum solutions can be matched to an expanding FLRW background. This includes the absence of any Einstein-Straus-like embeddings of the Schwarzschild exterior solution in FLRW spacetimes.

A surprising exact result for the Einstein Field Equations is that if pressure-free matter is moving in a shear-free way, then it must be either expansion-free or rotation-free. It has been suggested this result is also true for any barotropic perfect fluid, but a proof has remained elusive. We consider the case of barotropic perfect fluid solutions linearized about a Robertson-Walker geometry, and prove that the result remains true except for the case of a specific highly non-linear equation of state. We argue that this equation of state is non-physical, and hence the result is true in the linearized case for all physically realistic barotropic perfect fluids. This result, which is not true in Newtonian cosmology, demonstrates that the linearized solutions, believed to result in standard local Newtonian theory, do not always give the usual behaviour of Newtonian solutions.

We present a framework for the study of lensing in spherically symmetric spacetimes within the context of f(R) gravity. Equations for the propagation of null geodesics, together with an expression for the bending angle are derived for any f(R) theory and then applied to an exact spherically symmetric solution of R^n gravity. We find that for this case more bending is expected for R^n gravity theories in comparison to GR and is dependent on the value of n and the value of distance of closest approach of the incident null geodesic.

We develop a new covariant formalism to treat spherically symmetric spacetimes in metric} f(R) theories of gravity. Using this formalism we derive the general equations for a static and spherically symmetric metric in a general f(R)-gravity. These equations are used to determine the conditions for which the Schwarzschild metric is the only vacuum solution with vanishing Ricci scalar. We also show that our general framework provides a clear way of showing that the Schwarzschild solution is not a unique static spherically symmetric solution, providing some incite on how the current form of Birkhoff's theorem breaks down for these theories.