We examine the effect that the magnetic part of the Weyl tensor has on the large-scale expansion of space. This is done within the context of a class of cosmological models that contain regularly arranged discrete masses, rather than a continuous perfect fluid. The natural set of geodesic curves that one should use to consider the cosmological expansion of these models requires the existence of a non-zero magnetic part of the Weyl tensor. We include this object in the evolution equations of these models by performing a Taylor series expansion about a hypersurface where it initially vanishes. At the same cosmological time, measured as a fraction of the age of the universe, we find that the influence of the magnetic part of the Weyl tensor increases as the number of masses in the universe is increased. We also find that the influence of the magnetic part of the Weyl tensor increases with time, relative to the leading-order electric part, so that its contribution to the scale of the universe can reach values of ~1%, before the Taylor series approximation starts to break down.

We consider a family of cosmological models in which all mass is confined to a regular lattice of identical black holes. By exploiting the reflection symmetry about planes that bisect these lattices into identical halves, we are able to consider the evolution of a number of geometrically distinguished surfaces that exist within each of them. We find that the evolution equations for the reflection symmetric surfaces can be written as a simple set of Friedmann-like equations, with source terms that behave like a set of interacting effective fluids. We then show that gravitational waves are effectively trapped within small chambers for all time, and are not free to propagate throughout the space-time. Each chamber therefore evolves as if it were in isolation from the rest of the universe. We call this phenomenon "piecewise silence".

We study the effects of inhomogeneities on the evolution of the Universe, by considering a range of cosmological models with discretized matter content. This is done using exact and fully relativistic methods that exploit the symmetries in and about submanifolds of spacetimes that themselves possess no continuous global symmetries. These methods allow us to follow the evolution of our models throughout their entire history, far beyond what has previously been possible. We find that while some space-like curves collapse to anisotropic singularities in finite time, others remain non-singular forever. The resulting picture is of a cosmological spacetime in which some behaviour remains close to Friedmann-like, while other behaviours deviate radically. In particular, we find that large-scale acceleration is possible without any violation of the energy conditions.

An important open question in cosmology is the degree to which the Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions of Einstein's equations are able to model the large-scale behaviour of the locally inhomogeneous observable universe. We investigate this problem by considering a range of exact n-body solutions of Einstein's constraint equations. These solutions contain discrete masses, and so allow arbitrarily large density contrasts to be modelled. We restrict our study to regularly arranged distributions of masses in topological 3-spheres. This has the benefit of allowing straightforward comparisons to be made with FLRW solutions, as both spacetimes admit a discrete group of symmetries. It also provides a time-symmetric hypersurface at the moment of maximum expansion that allows the constraint equations to be solved exactly. We find that when all the mass in the universe is condensed into a small number of objects (<10) then the amount of backreaction in dust models can be large, with O(1) deviations from the predictions of the corresponding FLRW solutions. When the number of masses is large (>100), however, then our measures of backreaction become small (<1%). This result does not rely on any averaging procedures, which are notoriously hard to define uniquely in general relativity, and so provides (to the best of our knowledge) the first exact and unambiguous demonstration of backreaction in general relativistic cosmological modelling. Discrete models such as these can therefore be used as laboratories to test ideas about backreaction that could be applied in more complicated and realistic settings.

{\it Fermi Gamma-ray Space Telescope} observations of GRB110721A have revealed two emission components from the relativistic jet: emission from the photosphere, peaking at $\sim 100$ keV and a non-thermal component, which peaks at $\sim 1000$ keV. We use the photospheric component to calculate the properties of the relativistic outflow. We find a strong evolution in the flow properties: the Lorentz factor decreases with time during the bursts from $\Gamma \sim 1000$ to $\sim 150$ (assuming a redshift $z=2$; the values are only weakly dependent on unknown efficiency parameters). Such a decrease is contrary to the expectations from the internal shocks and the isolated magnetar birth models. Moreover, the position of the flow nozzle measured from the central engine, $r_0$, increases by more than two orders of magnitude. Assuming a moderately magnetised outflow we estimate that $r_0$ varies from $10^6$ cm to $\sim 10^9$ cm during the burst. We suggest that the maximal value reflects the size of the progenitor core. Finally, we show that these jet properties naturally explain the observed broken power-law decay of the temperature which has been reported as a characteristic for GRB pulses.

In dynamic spacetimes in which asymmetric gravitational collapse/expansion is taking place, the timelike geodesic equation appears to exhibit an interesting property: Relative to the collapsing configuration, free test particles undergo gravitational "acceleration" and form a double-jet configuration parallel to the axis of collapse. We illustrate this aspect of peculiar motion in simple spatially homogeneous cosmological models such as the Kasner spacetime. To estimate the effect of spatial inhomogeneities on cosmic jets, timelike geodesics in the Ricci-flat double-Kasner spacetime are studied in detail. While spatial inhomogeneities can significantly modify the structure of cosmic jets, we find that under favorable conditions the double-jet pattern can initially persist over a finite period of time for sufficiently small inhomogeneities.

We discuss time-dependent gravitational fields that "accelerate" free test particles to the speed of light resulting in cosmic double-jet configurations. It turns out that complete gravitational collapse along a spatial axis together with corresponding expansion along the other two axes leads to the accelerated motion of free test particles up and down parallel to the collapse axis such that a double-jet pattern is asymptotically formed with respect to the collapsed configuration.

A new formulation of Carter's constant for geodesic motion in Kerr black holes is given. It is shown that Carter's constant corresponds to the total angular momentum plus a precisely defined part which is quadratic in the linear momenta. The characterization is exact in the weak field limit obtained by letting the gravitational constant go to zero. It is suggested that the new form can be useful in current studies of the dynamics of extreme mass ratio inspiral (EMRI) systems emitting gravitational radiation.

It is shown that Einstein gravity tends to modify the electric and magnetic fields appreciably at distances of the order of the Compton wavelength. At that distance the gravitational field becomes spin dominated rather than mass dominated. The gravitational field couples to the electromagnetic field via the Einstein-Maxwell equations which in the simplest model causes the electrostatic field of charged spinning particles to acquire an oblate structure relative to the spin direction. For electrons and protons, a pure Coulomb field is therefore likely to be incompatible with general relativity at the Compton scale. In the simplest model, the magnetic dipole corresponds to the Dirac g-factor, g=2. Also, it follows from the form of the electric field that the electric dipole moment vanishes, in agreement with current experimental limits for the electron. Quantitatively, the classical Einstein-Maxwell theory predicts the magnetic and electric dipoles of the electron to an accuracy of about one part in 10^{-3} or better. Going to the next multipole order, one finds that the first non-vanishing higher multipole is the electric quadrupole moment which is predicted to be -124 barn for the electron. Any non-zero value of the electric quadrupole moment for the electron or the proton would be a clear sign of curvature due to the implied violation of rotation invariance. There is also a possible spherical modification of the Coulomb force proportional to r^{-4}. However, the size of this effect is well below current experimental limits. The corrections to the hydrogen spectrum are expected to be small but possibly detectable.

It is shown that optical geometry of the Reissner-Nordstrom exterior metric can be embedded in a hyperbolic space all the way down to its outer horizon. The adopted embedding procedure removes a breakdown of flat-space embeddings which occurs outside the horizon, at and below the Buchdahl-Bondi limit (R/M=9/4 in the Schwarzschild case). In particular, the horizon can be captured in the optical geometry embedding diagram. Moreover, by using the compact Poincare ball representation of the hyperbolic space, the embedding diagram can cover the whole extent of radius from spatial infinity down to the horizon. Attention is drawn to advantages of such embeddings in an appropriately curved space: this approach gives compact embeddings and it distinguishes clearly the case of an extremal black hole from a non-extremal one in terms of topology of the embedded horizon.

We discuss ultracompact stellar objects which have multiple necks in their optical geometry. There are in fact physically reasonable equations of state for which the number of necks can be arbitrarily large. The proofs of these statements rely on a recent regularized formulation of the field equations for static spherically symmetric models due to Nilsson and Uggla. We discuss in particular the equation of state p=\rho-\rho_s which plays a central role in this context.

A family of exact relativistic stellar models is described. The family generalizes Buchdahl's n=1 polytropic solution. The matter content is a perfect fluid and, excluding Buchdahl's original model, it behaves as a liquid at low pressures in the sense that the energy density is non-zero in the zero pressure limit. The equation of state has two free parameters, a scaling and a stiffness parameter. Depending on the value of the stiffness parameter the fluid behaviour can be divided in four different types. Physical quantities such as masses, radii and surface redshifts as well as density and pressure profiles are calculated and displayed graphically. Leaving the details to a later publication, it is noted that one of the equation of state types can quite accurately approximate the equation of state of real cold matter in the outer regions of neutron stars. Finally, it is observed that the given equation of state does not admit models with a conical singularity at the center.

The possibility of trapped modes of gravitational waves appearing in stars with R>3M is considered. It is shown that the restriction to R<3M in previous studies of trapped modes, using uniform density models, is not essential. Scattering potentials are computed for another family of analytic stellar models showing the appearance of a deep potential well for one model with R>3M. However, the provided example is unstable, although it has a more realistic equation of state in the sense that the sound velocity is finite. On the other hand it is also shown that for some stable models belonging to the same family but having R<3M, the well is significantly deeper than that of the uniform density stars. Whether there are physically realistic equations of state which allow stable configurations with trapped modes therefore remains an open problem.