Stephen Hawking's contributions to the understanding of gravity, black holes and cosmology were truly immense. They began with the singularity theorems in the 1960s followed by his discovery that black holes have an entropy and consequently a finite temperature. Black holes were predicted to emit thermal radiation, what is now called Hawking radiation. He pioneered the study of primordial black holes and their potential role in cosmology. His organisation of and contributions to the Nuffield Workshop in 1982 consolidated the picture that the large-scale structure of the universe originated as quantum fluctuations during the inflationary era. Work on the interplay between quantum mechanics and general relativity resulted in his formulation of the concept of the wavefunction of the universe. The tension between quantum mechanics and general relativity led to his struggles with the information paradox concerning deep connections between these fundamental areas of physics. These achievements were all accomplished following the diagnosis during the early years of Stephen's studies as a post-graduate student in Cambridge that he had incurable motor neuron disease -- he was given two years to live. Against all the odds, he lived a further 55 years. The distinction of his work led to many honours and he became a major public figure, promoting with passion the needs of disabled people. His popular best-selling book A Brief History of Time made cosmology and his own work known to the general public worldwide. He became an icon for science and an inspiration to all.

We study a free scalar field $\phi$ in a fixed curved background spacetime subject to a higher derivative field equation of the form $F(\Box)\phi =0$, where $F$ is a polynomial of the form $F(\Box)= \prod_i (\Box-m_i^2)$ and all masses $m_i$ are distinct and real. Using an auxiliary field method to simplify the calculations, we obtain expressions for the Belinfante-Rosenfeld symmetric energy-momentum tensor and compare it with the canonical energy-momentum tensor when the background is Minkowski spacetime. We also obtain the conserved symplectic current necessary for quantisation and briefly discuss the issue of negative energy versus negative norm and its relation to Reflection Positivity in Euclidean treatments. We study, without assuming spherical symmetry, the possible existence of finite energy static solutions of the scalar equations, in static or stationary background geometries. Subject to various assumptions on the potential, we establish non-existence results including a no-scalar-hair theorem for static black holes. We consider Pais-Uhlenbeck field theories in a cosmological de Sitter background, and show how the Hubble friction may eliminate what would otherwise be unstable behaviour when interactions are included.

The Eisenhart lift of a Paul Trap used to store ions in molecular physics is a linearly polarized periodic gravitational wave. A modified version of Dehmelt's Penning Trap is in turn related to circularly polarized periodic gravitational waves, sought for in inflationary models. Similar equations rule also the Lagrange points in Celestial Mechanics. The explanation is provided by anisotropic oscillators.

Circularly polarized gravitational sandwich waves exhibit, as do their linearly polarized counterparts, the Velocity Memory Effect: freely falling test particles in the flat after-zone fly apart along straight lines with constant velocity. In the inside zone their trajectories combine oscillatory and rotational motions in a complicated way. For circularly polarized periodic gravitational waves some trajectories remain bounded, while others spiral outward. These waves admit an additional "screw" isometry beyond the usual five. The consequences of this extra symmetry are explored.

A cosmological extension of the Eisenhart-Duval metric is constructed by incorporating a cosmic scale factor and the energy-momentum tensor into the scheme. The dynamics of the spacetime is governed the Ermakov-Milne-Pinney equation. Killing isometries include spatial translations and rotations, Newton--Hooke boosts and translation in the null direction. Geodesic motion in Ermakov-Milne-Pinney cosmoi is analyzed. The derivation of the Ermakov-Lewis invariant, the Friedmann equations and the Dmitriev-Zel'dovich equations within the Eisenhart--Duval framework is presented.

The "gravitational memory effect" due to an exact plane wave provides us with an elementary description of the diffeomorphisms associated with soft gravitons. It is explained how the presence of the latter may be detected by observing the motion of freely falling particles or other forms of gravitational wave detection. Numerical calculations confirm the relevance of the first, second and third time integrals of the Riemann tensor pointed out earlier. Solutions for various profiles are constructed. It is also shown how to extend our treatment to Einstein-Maxwell plane waves and a midi-superspace quantization is given.

We give an account of the gravitational memory effect in the presence of the exact plane wave solution of Einstein's vacuum equations. This allows an elementary but exact description of the soft gravitons and how their presence may be detected by observing the motion of freely falling particles. The theorem of Bondi and Pirani on caustics (for which we present a new proof) implies that the asymptotic relative velocity is constant but not zero, in contradiction with the permanent displacement claimed by Zel'dovich and Polnarev. A non-vanishing asymptotic relative velocity might be used to detect gravitational waves through the "velocity memory effect", considered by Braginsky, Thorne, Grishchuk, and Polnarev.

The equations of null geodesics in the STU family of rotating black hole solutions of supergravity theory, which may be considered as deformations of the vacuum Kerr metric, are completely integrable. We propose that they be used as a foil to test, for example, with what precision the gravitational field external to the black hole at the centre of our galaxy is given by the Kerr metric. By contrast with some metrics proposed in the literature, the STU metrics satisfy by construction the dominant and strong energy conditions. Our considerations may be extended to include the effects of a cosmological term. We show that these metrics permit a straightforward calculation of the properties of black hole shadows.

Conjectures play a central role in theoretical physics, especially those that assert an upper bound to some dimensionless ratio of physical quantities. In this paper we introduce a new such conjecture bounding the ratio of the magnetic moment to angular momentum in nature. We also discuss the current status of some old bounds on dimensionless and dimensional quantities in arbitrary spatial dimension. Our new conjecture is that the dimensionless Schuster-Wilson-Blackett number, c{\mu}/JG^{(1/2)}, where {\mu} is the magnetic moment and J is the angular momentum, is bounded above by a number of order unity. We verify that such a bound holds for charged rotating black holes in those theories for which exact solutions are available, including the Einstein-Maxwell theory, Kaluza-Klein theory, the Kerr-Sen black hole, and the so-called STU family of charged rotating supergravity black holes. We also discuss the current status of the Maximum Tension Conjecture, the Dyson Luminosity Bound, and Thorne's Hoop Conjecture.

We study closed photon orbits in spherically-symmetric static solutions of supergravity theories, a Horndeski theory, and a theory of quintessence. These orbits lie in what we shall call a photon sphere (anti-photon sphere) if the orbit is unstable (stable). We show that in all the asymptotically flat solutions we examine that admit a regular event horizon, and whose energy-momentum tensor satisfies the strong energy condition, there is one and only one photon sphere outside the event horizon. We give an example of a Horndeski theory black hole (whose energy-momentum tensor violates the strong energy condition) whose metric admits both a photon sphere and an anti-photon sphere. The uniqueness and non-existence also holds for asymptotically anti-de Sitter solutions in gauged supergravity. The latter also exhibit the projective symmetry that was first discovered for the Schwarzschild-de Sitter metrics: the unparameterised null geodesics are the same as when the cosmological or gauge coupling constant vanishes. We also study the closely related problem of accretion flows by perfect fluids in these metrics. For a radiation fluid, Bondi's sonic horizon coincides with the photon sphere. For a general polytropic equation of state this is not the case. Finally we exhibit counterexamples to a conjecture of Hod's.

In this work a series of methods are developed for understanding the Friedmann equation when it is beyond the reach of the Chebyshev theorem. First it will be demonstrated that every solution of the Friedmann equation admits a representation as a roulette such that information on the latter may be used to obtain that for the former. Next the Friedmann equation is integrated for a quadratic equation of state and for the Randall--Sundrum II universe, leading to a harvest of a rich collection of new interesting phenomena. Finally an analytic method is used to isolate the asymptotic behavior of the solutions of the Friedmann equation, when the equation of state is of an extended form which renders the integration impossible, and to establish a universal exponential growth law.

This paper is a continuation of our earlier study on the integrability of the Friedmann equations in the light of the Chebyshev theorem. Our main focus will be on a series of important, yet not previously touched, problems when the equation of state for the perfect-fluid universe is nonlinear. These include the generalized Chaplygin gas, two-term energy density, trinomial Friedmann, Born--Infeld, and two-fluid models. We show that some of these may be integrated using Chebyshev's result while other are out of reach by the theorem but may be integrated explicitly by other methods. With the explicit integration, we are able to understand exactly the roles of the physical parameters in various models play in the cosmological evolution. For example, in the Chaplygin gas universe, it is seen that, as far as there is a tiny presence of nonlinear matter, linear matter makes contribution to the dark matter, which becomes significant near the phantom divide line. The Friedmann equations also arise in areas of physics not directly related to cosmology. We provide some examples ranging from geometric optics and central orbits to soap films and the shape of glaciated valleys to which our results may be applied.

We discuss various examples and ramifications of the conjecture that there exists a maximum force (or tension) in general relativistic systems. We contrast this situation with that in Newtonian gravity, where no maximum force exists, and relate it to the existence of natural units defined by constants of Nature and the fact that the Planck units of force and power do not depend on Planck's constant. We discuss how these results change in higher dimensions where the Planck units of force are no longer non-quantum. We discuss the changes that might occur to the conjecture if a positive cosmological constant exists and derive a maximum force bound using the Kottler-Schwarzschild-de Sitter black hole.

This short but systematic work demonstrates a link between Chebyshev's theorem and the explicit integration in cosmological time $t$ and conformal time $\eta$ of the Friedmann equations in all dimensions and with an arbitrary cosmological constant $\Lambda$. More precisely, it is shown that for spatially flat universes an explicit integration in $t$ may always be carried out, and that, in the non-flat situation and when $\Lambda$ is zero and the ratio $w$ of the pressure and energy density in the barotropic equation of state of the perfect-fluid universe is rational, an explicit integration may be carried out if and only if the dimension $n$ of space and $w$ obey some specific relations among an infinite family. The situation for explicit integration in $\eta$ is complementary to that in $t$. More precisely, it is shown in the flat-universe case with $\Lambda\neq0$ that an explicit integration in $\eta$ can be carried out if and only if $w$ and $n$ obey similar relations among a well-defined family which we specify, and that, when $\Lambda=0$, an explicit integration can always be carried out whether the space is flat, closed, or open. We also show that our method may be used to study more realistic cosmological situations when the equation of state is nonlinear.

In a previous paper [arXiv:1308.1852] we showed how a finite system of discrete particles interacting with each other via Newtonian gravitational attraction would lead to precisely the same dynamical equations for homothetic motion as in the case of the pressure-free Friedmann-Lema\^{i}tre-Robertson-Walker cosmological models of General Relativity Theory, provided the distribution of particles obeys the central configuration equation. In this paper we show one can obtain perturbed such Newtonian solutions that give the same linearised structure growth equations as in the general relativity case. We also obtain the Dmitriev-Zeldovich equations for subsystems in this discrete gravitational model, and show how it leads to the conclusion that voids have an apparent negative mass.

Theories with a time dependent Newton's constant admit two natural measures of time : atomic and astronomical. Temporal parametrisation by SL(2, R) transformations gives rise to an equivalence between theories with different time dependence's, including the special Case of no time dependence, a fact noticed by Mestschersky, Vinti and by Lynden-Bell. I point out that theories with time dependent dark energy densities admit three natural measures of time : atomic and astronomical and de Sitter related by temporal re-parametrizations and I extend Mestschersky-Vinti-Lynden-Bell's result to cover this more general situation. I find a consequent equivalence between theories in which the density of dark energy is constant in time and in which it varies with time. Strikingly a time dependent cosmological constant changes by the addition of a Schwarzian derivative term unless the temporal reparameterization belongs to SL(2, R). In General Relativity one may introduce a Schwarzian tensor to investigate how the notion of dark energy changes under changes of conformal frame. The general theory is illustrated in the case of Friedmann-Lemaitre metrics.

Starting from the geometrical concept of a 4-dimensional de-Sitter configuration of spheres in Euclidean 3-space and modelling voids in the Universe as spheres, we show that a uniform distribution over this configuration space implies a power-law for the void number density which is consistent with results from the excursion set formalism and with data, for an intermediate range of void volumes. The scaling dimension of the large scale structure can be estimated as well. We also discuss the effect of restricting the survey geometry on the void statistics. This work is a new application of de-Sitter geometry to cosmology and also provides a new geometrical perspective on self-similarity in cosmology.

In this paper we lay down the foundations for a purely Newtonian theory of cosmology, valid at scales small compared with the Hubble radius, using only Newtonian point particles acted on by gravity and a possible cosmological term. We describe the cosmological background which is given by an exact solution of the equations of motion in which the particles expand homothetically with their comoving positions constituting a central configuration. We point out, using previous work, that an important class of central configurations are homogeneous and isotropic, thus justifying the usual assumptions of elementary treatments. The scale factor is shown to satisfy the standard Raychaudhuri and Friedmann equations without making any fluid dynamic or continuum approximations. Since we make no commitment as to the identity of the point particles, our results are valid for cold dark matter, galaxies, or clusters of galaxies. In future publications we plan to discuss perturbations of our cosmological background from the point particle viewpoint laid down in this paper and show consistency with much standard theory usually obtained by more complicated and conceptually less clear continuum methods. Apart from its potential use in large scale structure studies, we believe that out approach has great pedagogic advantages over existing elementary treatments of the expanding universe, since it requires no use of general relativity or continuum mechanics but concentrates on the basic physics: Newton's laws for gravitationally interacting particles.

The spacetimes obtained by Ernst's procedure for appending an external magnetic field $B$ to a seed Kerr-Newman black hole are commonly believed to be asymptotic to the static Melvin solution. We show that this is not in general true. Unless the electric charge of the black hole satisfies $Q= jB(1+ 1/4 j^2 B^4)$, where $j$ is the angular momentum of the original seed solution, an ergoregion extends all the way from the black hole horizon to infinity. We give a self-contained account of the solution-generating procedure, including including explicit formulae for the metric and the vector potential. In the case when $Q= jB(1+ 1/4 j^2 B^4)$, we show that there is an arbitrariness in the choice of asymptotically timelike Killing field $K_\Omega= \partial/\partial t+ \Omega \partial/\partial \phi$, because there is no canonical choice of $\Omega$. For one choice, $\Omega=\Omega_s$, the metric is asymptotically static, and there is an ergoregion confined to the neighbourhood of the horizon. On the other hand, by choosing $\Omega=\Omega_H$, so that $K_{\Omega_H}$ is co-rotating with the horizon, then for sufficiently large $B$ numerical studies indicate there is no ergoregion at all. For smaller values, in a range $B_-<B<B_+$, there is a toroidal ergoregion outside and disjoint from the horizon. If $B\le B_-$ this ergoregion expands all the way to infinity in a cylindrical region near to the rotation axis. For black holes whose size is small compared to the Melvin radius 2/B, we recover Wald's result that it is energetically favourable for the hole to acquire a charge $2jB$.

The cumulative size-frequency distributions of impact craters on planetary bodies in the solar system appear to approximate a universal inverse square power-law for small crater radii. In this article, we show how this distribution can be understood easily in terms of geometrical statistics, using a de-Sitter geometry of the configuration space of circles on the Euclidean plane and on the unit sphere. The effect of crater overlap is also considered.

Hill's equations, which first arose in the study of the Earth-Moon-Sun system, admit the two-parameter centrally extended Newton-Hooke symmetry without rotations. This symmetry allows for extending Kohn's theorem about the center-of-mass decomposition. Particular light is shed on the problem using Duval's "Bargmann" framework. The separation of the center-of-mass motion into that of a guiding center and relative motion is derived by a generalized chiral decomposition.

We elucidate the mechanism where by gravitational fluctuations are frozen in during a period of inflation but nevertheless within the past horizon of any geodesic observer the spacetime settles down exponentially fast to the de-Sitter metric.

An exact solution representing black holes in an expanding universe is found. The black holes are maximally charged and the universe is expanding with arbitrary equation of state. It is an exact solution of the Einstein-scalar-Maxwell system, in which we have two Maxwell-type U(1) fields coupled to the scalar field. The potential of the scalar field is an exponential. We find a regular horizon, which depends on one parameter (the ratio of the energy density of U(1) fields to that of the scalar field). The horizon is static because of the balance on the horizon between gravitational attractive force and U(1) repulsive force acting on the scalar field. We also calculate the black hole temperature.

Nurowski [arXiv:1003.1503] has recently suggested a link between the observation of Dark Energy in cosmology and the projective equivalence of certain Friedman-Lemaitre-Robertson-Walker (FLRW) metrics. Specifically, he points out that two FLRW metrics with the same unparameterized geodesics have their energy densities differing by a constant. From this he queries whether the existence of dark energy is meaningful. We point out that physical observables in cosmology are not projectively invariant and we relate the projective symmetry uncovered by Nurowski to some previous work on projective equivalence in cosmology.

We clarify the procedure for expressing the Friedmann equation in terms of directly measurable cosmological scalars constructed out of higher derivatives of the scale factor. We carry out this procedure for pure dust, Chaplygin gas and generalised Chaplygin gas energy-momentum tensors. In each case it leads to a constraint on the scalars thus giving rise to a test of General Relativity. We also discuss a formulation of the Friedmann equation as unparametrised geodesic motion and its connection with the Lagrangian treatment of perfect fluids coupled to gravity.

We interpret the well known fact that the equations for light rays in the Kottler or Schwarzschild-de Sitter metric are independent of the cosmological constant in terms of the projective equivalence of the optical metric for any value of \Lambda. We explain why this does not imply that lensing phenomena are independent of \Lambda. Motivated by this example, we find a large collection of one-parameter families of projectively equivalent metrics including both the Kottler optical geometry and the constant curvature metrics as special cases. Using standard constructions for geodesically equivalent metrics we find classical and quantum conserved quantities and relate these to known quantities.

In this geometrical approach to gravitational lensing theory, we apply the Gauss-Bonnet theorem to the optical metric of a lens, modelled as a static, spherically symmetric, perfect non-relativistic fluid, in the weak deflection limit. We find that the focusing of the light rays emerges here as a topological effect, and we introduce a new method to calculate the deflection angle from the Gaussian curvature of the optical metric. As examples, the Schwarzschild lens, the Plummer sphere and the singular isothermal sphere are discussed within this framework.

We obtain an exact solution of the supergravity equations of motion in which the four-dimensional observed universe is one of a number of colliding D3-branes in a Calabi-Yau background. The collision results in the ten-dimensional spacetime splitting into disconnected regions, bounded by curvature singularities. However, near the D3-branes the metric remains static during and after the collision. We also obtain a general class of solutions representing $p$-brane collisions in arbitrary dimensions, including one in which the universe ends with the mutual annihilation of a positive-tension and negative-tension 3-brane.

In a previous paper we obtained formulae for the volume of a causal diamond or Alexandrov open set $I^+(p) \cap I^-(q)$ whose duration $\tau(p,q) $ is short compared with the curvature scale. In the present paper we obtain asymptotic formulae valid when the point $q$ recedes to the future boundary ${\cal I}^+$ of an asymptotically de-Sitter spacetime. The volume (at fixed $\tau$) remains finite in this limit and is given by the universal formula $V(\tau) = {4\over 3}\pi (2\ln \cosh{\tau\over 2}-\tanh^2{\tau\over 2})$ plus corrections (given by a series in $e^{-t_q}$) which begin at order $e^{-4t_q}$. The coefficents of the corrections depend on the geometry of ${\cal I}^+$. This behaviour is shown to be consistent with the no-hair property of cosmological event horizons and with calculations of de-Sitter quasinormal modes in the literature.

We study the behaviour of five-dimensional fermions localized on branes, which we describe by domain walls, when two parallel branes collide in a five-dimensional Minkowski background spacetime. We find that most fermions are localized on both branes as a whole even after collision. However, how much fermions are localized on which brane depends sensitively on the incident velocity and the coupling constants unless the fermions exist on both branes.

The Hamiltonian structure of general relativity provides a natural canonical measure on the space of all classical universes, i.e., the multiverse. We review this construction and show how one can visualize the measure in terms of a "magnetic flux" of solutions through phase space. Previous studies identified a divergence in the measure, which we observe to be due to the dilatation invariance of flat FRW universes. We show that the divergence is removed if we identify universes which are so flat they cannot be observationally distinguished. The resulting measure is independent of time and of the choice of coordinates on the space of fields. We further show that, for some quantities of interest, the measure is very insensitive to the details of how the identification is made. One such quantity is the probability of inflation in simple scalar field models. We find that, according to our implementation of the canonical measure, the probability for N e-folds of inflation in single-field, slow-roll models is suppressed by of order exp(-3N) and we discuss the implications of this result.

We construct exact time-dependent solutions of the supergravity equations of motion in which two initially non-singular branes, one with positive and the other with negative tension, move together and annihilate each other in an all-enveloping spacetime singularity. Among our solutions are the Horava-Witten solution of heterotic M-theory and a Randall-Sundrum I type solution, both of which are supersymmetric, i.e. BPS, in the time-independent case. In the absence of branes our solutions are of Kasner type, and the source of instability may ascribed to a failure to stabilise some of the modulus fields of the compactification. It also raises questions about the viability of models based on some sorts of negative tension brane.

Some speculative preliminary ideas relating matrix theory and cosmology are discussed.

We show explicitly how the Newton-Hooke groups act as symmetries of the equations of motion of non-relativistic cosmological models with a cosmological constant. We give the action on the associated non-relativistic spacetimes and show how these may be obtained from a null reduction of 5-dimensional homogeneous pp-wave Lorentzian spacetimes. This allows us to realize the Newton-Hooke groups and their Bargmann type central extensions as subgroups of the isometry groups of the pp-wave spacetimes. The extended Schrodinger type conformal group is identified and its action on the equations of motion given. The non-relativistic conformal symmetries also have applications to time-dependent harmonic oscillators. Finally we comment on a possible application to Gao's generalization of the matrix model.

We consider the equilibria of point particles under the action of two body central forces in which there are both repulsive and attractive interactions, often known as central configurations, with diverse applications in physics, in particular as homothetic time-dependent solutions to Newton's equations of motion and as stationary states in the One Component Plasma model. Concentrating mainly on the case of an inverse square law balanced by a linear force, we compute numerically equilibria and their statistical properties. When all the masses (or charges) of the particles are equal, for small numbers of points they are regular convex deltahedra, which on increasing the number of points give way to a multi-shell structure. In the limit of a large number of points we argue using an analytic model that they form a homogeneous spherical distribution of points, whose spatial distribution appears, from our preliminary investigation, to be similar to that of a Bernal hard-sphere liquid.

A set of consistency conditions is derived from Einstein equations for brane world scenarios with a spatially periodic internal space. In particular, the sum of the total tension of the flat branes and the non-negative integral of the gradient energy of the bulk scalars must vanish. This constraint allows us to make a simple consistency check of several models. We show that the two-brane Randall-Sundrum model satisfies this constraint, but it does not allow a generalization with smooth branes (domain walls), independently of the issue of supersymmetry. The Goldberger-Wise model of brane stabilization has to include the backreaction on the metric and the fine tuning of the cosmological constant to satisfy the constraints. We check that this is achieved in the DeWolfe-Freedman-Gubser-Karch scenario. Our constraints are automatically satisfied in supersymmetric brane world models.

We show that smooth domain wall spacetimes supported by a scalar field separating two anti-de-Sitter like regions admit a single graviton bound state. Our analysis yields a fully non-linear supergravity treatment of the Randall-Sundrum model. Our solutions describe a pp-wave propagating in the domain wall background spacetime. If the latter is BPS, our solutions retain some supersymmetry. Nevertheless, the Kaluza-Klein modes generate ``pp curvature'' singularities in the bulk located where the horizon of AdS would ordinarily be.