In a cold matter universe, the linearized gravito-magnetic tensor field satisfies a transverse condition (vanishing divergence) when it is purely radiative. We show that in the nonlinear theory, it is no longer possible to maintain the transverse condition, since it leads to a non-terminating chain of integrability conditions. These conditions are highly restrictive, and are likely to hold only in models with special symmetries, such as the known Bianchi and $G_2$ examples. In models with realistic inhomogeneity, the gravito-magnetic field is necessarily non-transverse at second and higher order.

We revisit the issue of integrability conditions for the irrotational silent cosmological models. We formulate the problem both in 1+3 covariant and 1+3 orthonormal frame notation, and show there exists a series of constraint equations that need to be satisfied. These conditions hold identically for FLRW-linearised silent models, but not in the general exact non-linear case. Thus there is a linearisation instability, and it is highly unlikely that there is a large class of silent models. We conjecture that there are no spatially inhomogeneous solutions with Weyl curvature of Petrov type I, and indicate further issues that await clarification.

All spacetimes for an irrotational collisionless fluid with a purely electric Weyl tensor, with spacetime curvature determined by the exact Einstein field equations, are shown to be integrable. These solutions include the relativistic generalisations of the Zeldovich solutions of Newtonian theory. Thus our result shows the consistency of various studies of "silent" universes (where such consistency was assumed rather than proved).