A variety of gauges are used in cosmological perturbation theory. These are often chosen in order to attribute physical properties to a particular choice of coordinates, or otherwise to simplify the form of the resultant equations. Calculations are then performed with the understanding that they could have been done in any gauge, and that transformations between different gauges can be made at will. We show that this logic can be extended to the domain of large density contrasts, where different types of perturbative expansion are required, but that the way in which gauges can be chosen in the presence of such structures is severely constrained. In particular, most gauges that are commonly considered in the cosmology literature are found to be unviable in the presence of non-linear structures. This includes spatially flat gauge, synchronous gauge, comoving orthogonal gauge, total matter gauge, N-body gauge, and the uniform density gauge. In contrast, we find that the longitudinal gauge and the Newtonian motion gauge are both viable choices in both standard cosmological perturbation theory, and in the post-Newtonian perturbative expansions that are required in order to model non-linear structures.

The next generation of cosmological surveys will operate over unprecedented scales, and will therefore provide exciting new opportunities for testing general relativity. The standard method for modelling the structures that these surveys will observe is to use cosmological perturbation theory for linear structures on horizon-sized scales, and Newtonian gravity for non-linear structures on much smaller scales. We propose a two-parameter formalism that generalizes this approach, thereby allowing interactions between large and small scales to be studied in a self-consistent and well-defined way. This uses both post-Newtonian gravity and cosmological perturbation theory, and can be used to model realistic cosmological scenarios including matter, radiation and a cosmological constant. We find that the resulting field equations can be written as a hierarchical set of perturbation equations. At leading-order, these equations allow us to recover a standard set of Friedmann equations, as well as a Newton-Poisson equation for the inhomogeneous part of the Newtonian energy density in an expanding background. For the perturbations in the large-scale cosmology, however, we find that the field equations are sourced by both non-linear and mode-mixing terms, due to the existence of small-scale structures. These extra terms should be expected to give rise to new gravitational effects, through the mixing of gravitational modes on small and large scales - effects that are beyond the scope of standard linear cosmological perturbation theory. We expect our formalism to be useful for accurately modelling gravitational physics in universes that contain non-linear structures, and for investigating the effects of non-linear gravity in the era of ultra-large-scale surveys.

We propose and construct a two-parameter perturbative expansion around a Friedmann-Lema\^{i}tre-Robertson-Walker geometry that can be used to model high-order gravitational effects in the presence of non-linear structure. This framework reduces to the weak-field and slow-motion post-Newtonian treatment of gravity in the appropriate limits, but also includes the low-amplitude large-scale fluctuations that are important for cosmological modelling. We derive a set of field equations that can be applied to the late Universe, where non-linear structure exists on supercluster scales, and perform a detailed investigation of the associated gauge problem. This allows us to identify a consistent set of perturbed quantities in both the gravitational and matter sectors, and to construct a set of gauge-invariant quantities that correspond to each of them. The field equations, written in terms of these quantities, take on a relatively simple form, and allow the effects of small-scale structure on the large-scale properties of the Universe to be clearly identified. We find that inhomogeneous structures source the global expansion, that there exist new field equations at new orders, and that there is vector gravitational potential that is a hundred times larger than one might naively expect from cosmological perturbation theory. Finally, we expect our formalism to be of use for calculating relativistic effects in upcoming ultra-large-scale surveys, as the form of the gravitational coupling between small and large scales depends on the non-linearity of Einstein's equations, and occurs at what is normally thought of as first order in cosmological perturbations.