Linear perturbations of spherically symmetric spacetimes in general relativity are described by radial wave equations, with potentials that depend on the spin of the perturbing field. In previous work we studied the quasinormal mode spectrum of spacetimes for which the radial potentials are slightly modified from their general relativistic form, writing generic small modifications as a power-series expansion in the radial coordinate. We assumed that the perturbations in the quasinormal frequencies are linear in some perturbative parameter, and that there is no coupling between the perturbation equations. In general, matter fields and modifications to the gravitational field equations lead to coupled wave equations. Here we extend our previous analysis in two important ways: we study second-order corrections in the perturbative parameter, and we address the more complex (and realistic) case of coupled wave equations. We highlight the special nature of coupling-induced corrections when two of the wave equations have degenerate spectra, and we provide a ready-to-use recipe to compute quasinormal modes. We illustrate the power of our parametrization by applying it to various examples, including dynamical Chern-Simons gravity, Horndeski gravity and an effective field theory-inspired model.

The parameterised post-Newtonian (PPN) formalism has enabled stringent tests of static weak-field gravity in a theory-independent manner. Here we incorporate screening mechanisms of modified gravity theories into the framework by introducing an effective gravitational coupling and defining the PPN parameters as functions of position. To determine these functions we develop a general method for efficiently performing the post-Newtonian expansion in screened regimes. For illustration, we derive all the PPN functions for a cubic galileon and a chameleon model. We also analyse the Shapiro time delay effect for these two models and find no deviations from General Relativity insofar as the signal path and the perturbing mass reside in a screened region of space.

The Horndeski action is the most general scalar-tensor theory with at most second-order derivatives in the equations of motion, thus evading Ostrogradsky instabilities and making it of interest when modifying gravity at large scales. To pass local tests of gravity, these modifications predominantly rely on nonlinear screening mechanisms that recover Einstein's Theory of General Relativity in regions of high density. We derive a set of conditions on the four free functions of the Horndeski action that examine whether a specific model embedded in the action possesses an Einstein gravity limit or not. For this purpose, we develop a new and surprisingly simple scaling method that identifies dominant terms in the equations of motion by considering formal limits of the couplings that enter through the new terms in the modified action. This enables us to find regimes where nonlinear terms dominate and Einstein's field equations are recovered to leading order. Together with an efficient approximation of the scalar field profile, one can then further evaluate whether these limits can be attributed to a genuine screening effect. For illustration, we apply the analysis to both a cubic galileon and a chameleon model as well as to Brans-Dicke theory. Finally, we emphasise that the scaling method also provides a natural approach for performing post-Newtonian expansions in screened regimes.