A Bianchi IX Mixmaster spacetime is the most general spatially homogeneous solution of Einstein's equations and it can represent the space-averaged Universe. We introduce two novel mechanisms resulting in a Mixmaster Universe with non-singular bounces which are quasi-isotropic. A fluid with a non-linear equation of state allows non-singular bounces. Using negative anisotropic stresses successfully isotropises this Universe and mitigates the well known Mixmaster chaotic behaviour. Thus the Universe can be an eternal Mixmaster, going through an infinite series of different cycles separated by bounces, with a sizable fraction of cycles isotropic enough to be well approximated by a standard Friedmann-Lema\^itre-Robertson-Walker model from the radiation era onward.

We consider the robustness of small-field inflation in the presence of scalar field inhomogeneities. Previous numerical work has shown that if the scalar potential is flat only over a narrow interval, such as in commonly considered inflection-point models, even small-amplitude inhomogeneities present at the would-be onset of inflation at $\tau = \tau_i$ can disrupt the accelerated expansion. In this paper, we parametrise and evolve the inhomogeneities from an earlier time $\tau_{IC}$ at which the initial data were imprinted, and show that for a broad range of inflationary and pre-inflationary models, inflection-point inflation withstands initial inhomogeneities. We consider three classes of perturbative pre-inflationary solutions (corresponding to energetic domination by the scalar field kinetic term, a relativistic fluid, and isotropic negative curvature), and two classes of exact solutions to Einstein's equations with large inhomogeneities (corresponding to a stiff fluid with cylindrical symmetry, and anisotropic negative curvature). We derive a stability condition that depends on the Hubble scales $H(\tau_ i)$ and $H(\tau_{IC})$, and a few properties of the pre-inflationary cosmology. For initial data imprinted at the Planck scale, the absence of an inhomogeneous initial data problem for inflection-point inflation leads to a novel, lower limit on the tensor-to-scalar ratio.

What happens to the most general closed oscillating universes in general relativity? We sketch the development of interest in cyclic universes from the early work of Friedmann and Tolman to modern variations introduced by the presence of a cosmological constant. Then we show what happens in the cyclic evolution of the most general closed anisotropic universes provided by the Mixmaster universe. We show that in the presence of entropy increase its cycles grow in size and age, increasingly approaching flatness. But these cycles also grow increasingly anisotropic at their expansion maxima. If there is a positive cosmological constant, or dark energy, present then these oscillations always end and the last cycle evolves from an anisotropic inflexion point towards a de Sitter future of everlasting expansion.

We investigate the behaviour of bouncing Bianchi type IX `Mixmaster' universes in general relativity. This generalises all previous studies of the cyclic behaviour of closed spatially homogeneous universes with and without entropy increase. We determine the behaviour of models containing radiation by analytic and numerical integration and show that increase of radiation entropy leads to increasing cycle size and duration. We introduce a null energy condition violating ghost field to create a smooth, non-singular bounce of finite size at the end of each cycle and compute the evolution through many cycles with and without entropy increase injected at the start of each cycle. In the presence of increasing entropy we find that the cycles grow larger and longer and the dynamics approach flatness, as in the isotropic case. However, successive cycles become increasingly anisotropic at the expansion maxima which is dominated by the general-relativistic effects of anisotropic 3-curvature. However, it becomes positive after expansion drives the dynamics close enough to isotropy for the curvature to become positive and for gravitational collapse to ensue. In the presence of a positive cosmological constant, radiation and a ghost field we show that, for a very wide range of cosmological constant values, the growing oscillations always cease and the dynamics subsequently approach those of the isotropic de Sitter universe at late times. This model is not included in the scope of earlier cosmic no-hair theorems because the 3-curvature can be positive. In the case of negative cosmological constant, radiation and an ultra-stiff field (to create non-singular bounces) we show that a sequence of chaotic oscillations also occurs, with sensitive dependence on initial conditions. In all cases, we follow the oscillatory evolution of the scale factors, the shear, and the 3-curvature from cycle to cycle.

We study the behaviour of Bianchi class A universes containing an ultra-stiff isotropic ghost field and a fluid with anisotropic pressures which is also ultra-stiff on the average. This allows us to investigate whether cyclic universe scenarios, like the ekpyrotic model, do indeed lead to isotropisation on approach to a singularity (or bounce) in the presence of dominant ultra-stiff pressure anisotropies. We specialise to consider the closed Bianchi type IX universe and show that when the anisotropic pressures are stiffer on average than any isotropic ultra-stiff fluid then, if they dominate on approach to the singularity, it will be anisotropic. We include an isotropic ultra-stiff ghost fluid with negative energy density in order to create a cosmological bounce at finite volume in the absence of the anisotropic fluid. When the dominant anisotropic fluid is present it leads to an anisotropic cosmological singularity rather than an isotropic bounce. The inclusion of anisotropic stresses generated by collisionless particles in an anisotropically expanding universe is therefore essential for a full analysis of the consequences of a cosmological bounce or singularity in cyclic universes.