The intention of this paper is mainly two-fold. \textit{First}, we point out a striking numerical agreement between the bulk viscosity in the lepton era calculated by Husdal (2016) and our own calculations of the present-day bulk viscosity when the functional form is $\zeta\,\sim\sqrt{\rho}$. From a phenomenological point of view, we thus seem to have an ansatz for the viscosity which bridges the infancy of the Universe ($\sim 1$ s) with the present. This can also be looked upon as a kind of symmetry between the early-time cosmology and the present-day cosmology: it is quite remarkable that the kinetic theory-based bulk viscosity in the early universe and the experimentally-based bulk viscosity in the present universe can be covered by the same simple analytical formula. \textit{Second}, we consider the Kasner universe as a typical anisotropic model of Bianchi-type I, investigating whether this geometrical model is compatible with constant viscosity coefficients in the fluid. Perhaps surprisingly, the existence of a shear viscosity turns out to be incompatible with the Kasner model. By contrast, a bulk viscosity is non-problematic in the {\it isotropic} version of the model. In the special case of a Zel'dovich (stiff) fluid, the three equal exponents in the Kasner metric are even determined by the bulk viscosity alone, independent of the value of the fluid energy density. We also give a brief comparison with some other recent approaches to viscous cosmology.

The weak gravitational lensing formalism can be extended to the strong lensing regime by integrating a nonlinear version of the geodesic deviation equation. The resulting "roulette" expansion generalises the notion of convergence, shear and flexion to arbitrary order. The independent coefficients of this expansion are screen space gradients of the optical tidal tensor which approximates to the usual lensing potential in the weak field limit. From lensed images, knowledge of the roulette coefficients can in principle be inverted to reconstruct the mass distribution of a lens. In this paper, we simplify the roulette expansion and derive a family of recursion relations between the various coefficients, generalising the Kaiser-Squires relations beyond the weak-lensing regime.

We investigate a corner of the Bianchi models that has not received much attention: "extended FLRW models" (eFLRW) defined as a cosmological model with underlying anisotropic Bianchi geometry that nevertheless expands isotropically and can be mapped onto a reference FLRW model with the same expansion history. In order to investigate the stability and naturalness of such models in a dynamical systems context, we consider spatially homogeneous models that contain a massless scalar field $\varphi$ and a non-tilted perfect fluid obeying an equation of state $p=w\rho$. Remarkably, we find that matter anisotropies and geometrical anisotropies tend to cancel out dynamically. Hence, the expansion is asymptotically isotropic under rather general conditions. Although extended FLRW models require a special matter sector with anisotropies that are 'fine-tuned" relative to geometrical anisotropies, our analysis shows that such solutions are dynamically preferred attractors in general relativity. Specifically, we prove that all locally rotationally symmetric Bianchi type III universes with space-like $\nabla_\mu\varphi$ are asymptotically shear-free, for all $w\in[-1,1]$. Moreover, all shear-free equilibrium sets with anisotropic spatial curvature are proved to be stable with respect to all homogeneous perturbations for $w\geq -1/3$.

Why is the Universe so homogeneous and isotropic? We summarize a general study of a $\gamma$-law perfect fluid alongside an inhomogeneous, massless scalar gauge field (with homogeneous gradient) in anisotropic spaces with General Relativity. The anisotropic matter sector is implemented as a $j$-form (field-strength level), where $j\,\in\,\{1,3\}$, and the spaces studied are Bianchi space-times of solvable type. Wald's no-hair theorem is extended to include the $j$-form case. We highlight three new self-similar space-times: the Edge, the Rope and Wonderland. The latter solution is so far found to exist in the physical state space of types I,II, IV, VI$_0$, VI$_h$, VII$_0$ and VII$_h$, and is a global attractor in I and V. The stability analysis of the other types has not yet been performed. This paper is a summary of ~[1], with some remarks towards new results which will be further laid out in upcoming work.

In this paper the dynamics of free gauge fields in Bianchi type I-VII$_{h}$ space-times is investigated. The general equations for a matter sector consisting of a $p$-form field strength ($p\,\in\,\{1,3\}$), a cosmological constant ($4$-form) and perfect fluid in Bianchi type I-VII$_{h}$ space-times are computed using the orthonormal frame method. The number of independent components of a $p$-form in all Bianchi types I-IX are derived and, by means of the dynamical systems approach, the behaviour of such fields in Bianchi type I and V are studied. Both a local and a global analysis are performed and strong global results regarding the general behaviour are obtained. New self-similar cosmological solutions appear both in Bianchi type I and Bianchi type V, in particular, a one-parameter family of self-similar solutions,"Wonderland ($\lambda$)" appears generally in type V and in type I for $\lambda=0$. Depending on the value of the equation of state parameter other new stable solutions are also found ("The Rope" and "The Edge") containing a purely spatial field strength that rotates relative to the co-moving inertial tetrad. Using monotone functions, global results are given and the conditions under which exact solutions are (global) attractors are found.

We analyze characteristic properties of two different cosmological models: (i) a one-component dark energy model where the bulk viscosity $\zeta$ is associated with the fluid as a whole, and (ii) a two-component model where $\zeta$ is associated with a dark matter component $\rho_{\rm m}$ only, the dark energy component considered inviscid. Shear viscosity is omitted. We assume throughout the simple equation of state $p=w\rho$, with $w$ a constant. In the one-component model we consider two possibilities, either to take $\zeta$ proportional to the scalar expansion (equivalent to the Hubble parameter), in which case the evolution becomes critically dependent on the value of the small constant $\alpha=1+w$ and the magnitude of $\zeta$. Second, we consider the case $\zeta=~$const., where a de Sitter final stage is reached in the future. In the two-component model we consider only the case where the dark matter viscosity $\zeta_{\rm m}$ is proportional to the square of $\rho_{\rm m}$, where again a de Sitter form is found in the future. In this latter case the formalism is supplemented by a phase space analysis. As a general result of our considerations we suggest that a value $\zeta_0\sim 10^6~$Pa s for the present viscosity is reasonable, and that the two-component model seems to be favored.