We review current theoretical cosmology, including fundamental and mathematical cosmology and physical cosmology (as well as cosmology in the quantum realm), with an emphasis on open questions.

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.

In this paper we show in a covariant and gauge invariant way that in general relativity, tidal forces are actually a hidden form of gravitational waves. This must be so because gravitational effects cannot occur faster than the speed of light. Any two body gravitating system, where the bodies are orbiting around each other, may generate negligible gravitational waves, but it is via these waves that non-negligible tidal forces (causing shape distortions) act on these bodies. Although the tidal forces are caused by the electric part of the Weyl tensor, we transparently show that some small time varying magnetic part of the Weyl tensor with non zero curl must be present in the system that mediates the tidal forces via gravitational wave type effects.

In general relativity the parallel transfer of a vector around a closed curve in spacetime, or along two curves which together form a closed loop, usually results in a nonzero deficit angle between the vector's initial and final positions. We show that such holonomy in the McVittie spacetime, which represents a gravitating object imbedded in an expanding universe, can in principle be used to directly detect the expansion of the universe, for example by measuring changes in the components of a gyroscopic spin axis. Although such changes are of course small, they are large enough (\D S \sim 10^{-7}) that they could conceivably be measured if the real universe behaved like the McVittie spacetime. The real problem is that virialization will lead to domains decoupled from the global expansion on a scale much larger than that of the solar system, making such an experiment infeasible probably even in principle. Nevertheless the effect is of interest in relation to ongoing discussions, dating back at least to Einstein and Straus, which concern the relationship between the expansion of the universe and local systems.

This paper revisits the controversy concerning whether gravitational lensing effects make a significant difference to estimation of distance to the Cosmic Microwave Background last scattering surface in cosmology. A recent paper by Kaiser and Peacock supports a previous paper by Weinberg stating that such affects average to zero because of energy conservation. In this note, problems are pointed out in the Kaiser and Peacock analysis related to their choice of endpoint of integration, and to the `wrinkly surface' argument.

We present the particle creation probability rate around a general black hole as an outcome of quantum fluctuations. Using the uncertainty principle for these fluctuation, we derive a new ultraviolet frequency cutoff for the radiation spectrum of a dynamical black hole. Using this frequency cutoff, we define the probability creation rate function for such black holes. We consider a dynamical Vaidya model, and calculate the probability creation rate for this case when its horizon is in a slowly evolving phase. Our results show that one can expect the usual Hawking radiation emission process in the case of a dynamical black hole when it has a slowly evolving horizon. Moreover, calculating the probability rate for a dynamical black hole gives a measure of when Hawking radiation can be killed off by an incoming flux of matter or radiation. Our result strictly suggests that we have to revise the Hawking radiation expectation for primordial black holes that have grown substantially since they were created in the early universe. We also infer that this frequency cut off can be a parameter that shows the primordial black hole growth at the emission moment.

We establish the existence and find the necessary and sufficient conditions for a new class of solutions of Locally Rotationally Symmetric spacetimes that have non vanishing rotation and spatial twist simultaneously. We transparently show that the existence of such solutions demand non vanishing and bounded heat flux and these solutions are self similar. We provide a brief algorithm indicating how to solve the system of field equations with the given Cauchy data on an initial spacelike Cauchy surface. Finally we argue that these solutions can be used as a first approximation from spherical symmetry to study rotating, inhomogeneous, dynamic and radiating astrophysical stars.

We use the phase plane analysis technique of Madsen and Ellis to consider a universe with a true cosmological constant as well as a cosmological "constant" that is decaying. Time symmetric dynamics for the inflationary era allows eternally bouncing models to occur. Allowing for scalar field dynamic evolution, we find that if dark energy decays in the future, chaotic cyclic universes exist provided the spatial curvature is positive. This is particularly interesting in light of current observations which do not yet rule out either closed universes or possible evolution of the cosmological constant. We present only a proof of principle, with no definite claim on the physical mechanism required for the present dark energy to decay.

No. In a number of papers Green and Wald argue that the standard FLRW model approximates our Universe extremely well on all scales, except close to strong field astrophysical objects. In particular, they argue that the effect of inhomogeneities on average properties of the Universe (backreaction) is irrelevant. We show that this latter claim is not valid. Specifically, we demonstrate, referring to their recent review paper, that (i) their two-dimensional example used to illustrate the fitting problem differs from the actual problem in important respects, and it assumes what is to be proven; (ii) the proof of the trace-free property of backreaction is unphysical and the theorem about it fails to be a mathematically general statement; (iii) the scheme that underlies the trace-free theorem does not involve averaging and therefore does not capture crucial non-local effects; (iv) their arguments are to a large extent coordinate-dependent, and (v) many of their criticisms of backreaction frameworks do not apply to the published definitions of these frameworks. It is therefore incorrect to infer that Green and Wald have proven a general result that addresses the essential physical questions of backreaction in cosmology.

Among the different methods to derive particle creation, finding the quantum stress tensor expectation value gives a covariant quantity which can be used for examining the back-reaction issue. However this tensor also includes vacuum polarization in a way that depends on the vacuum chosen. Here we review different aspects of particle creation by looking at energy conservation and at the quantum stress tensor. It will be shown that in the case of general spherically symmetric black holes that have a \emph{dynamical horizon}, as occurs in a cosmological context, one cannot have pair creation on the horizon because this violates energy conservation. This confirms the results obtained in other ways in a previous paper [25]. Looking at the expectation value of the quantum stress tensor with three different definitions of the vacuum state, we study the nature of particle creation and vacuum polarization in black hole and cosmological models, and the associated stress energy tensors. We show that the thermal temperature that is calculated from the particle flux given by the quantum stress tensor is compatible with the temperature determined by the affine null parameter approach. Finally, it will be shown that in the spherically symmetric dynamic case, we can neglect the backscattering term and only consider the s-waves term near the future apparent horizon.

An astrophysical (cosmological) black hole forming in a cosmological context will be subject to a flux of infalling matter and radiation, which will cause the outer apparent horizon (a marginal trapping surface) to be spacelike \cite{ellisetal14}. As a consequence the radiation emitted close to the apparent horizon no longer arrives at infinity with a diverging redshift. Standard calculations of the emission of Hawking radiation then indicate that no blackbody radiation is emitted to infinity by the black hole in these circumstances, hence there will also then be no black hole evaporation process due to emission of such radiation as long as the matter flux is significant. The essential adiabatic condition (eikonal approximation) for black hole radiation gives a strong limit to the black holes that can emit Hawking radiation. We give the mass range for the black holes that can radiate, according to their cosmological redshift, for the special case of the cosmic blackbody radiation (CBR) influx (which exists everywhere in the universe). At a very late stage of black hole formation when the CBR influx decays away, the black hole horizon becomes first a slowly evolving horizon and then an isolated horizon; at that stage, black hole radiation will start. This study suggests that the primordial black hole evaporation scenario should be revised to take these considerations into account.

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.

This paper considers the effect of Hawking radiation on an eternal black hole - that is. a maximally extended Schwarzschild solution. Symmetry considerations that hold independent of the details of the emission mechanism show there is an inconsistency in the claim that such a blackhole evaporates away in a finite time. In essence: because the external domain is static, there is an infinite time available for the process to take place, so whenever the evaporation process is claimed to come to completion, it should have happened earlier. The problem is identified to lie in the claim that the locus of emission of Hawking radiation lies just outside the globally defined event horizon. Rather, the emission domain must be mainly located inside the event horizon, so most of the Hawking radiation ends up at this singularity rather than at infinity and the black hole never evaporates away. This result supports a previous claim [arXiv:1310.4771] that astrophysical black holes do not evaporate.

This paper argues that the effect of Hawking radiation on an astrophysical black hole situated in a realistic cosmological context is not total evaporation of the black hole; rather there will always be a remnant mass. The key point is that the locus of emission of Hawking radiation is not the globally defined event horizon. Rather the emission domain lies just outside a timelike Marginal Outer Trapped Surface that is locally defined. The emission domain is mainly located inside the event horizon. A spacelike singularity forms behind the event horizon, and most of the Hawking radiation ends up at this singularity rather than at infinity. Whether any Hawking radiation reaches infinity depends on the relation between the emission domain and the event horizon. From the outside view, even if radiation is seen as always being emitted, the black hole never evaporates away, rather its mass and entropy asymptote to finite non-zero limits, and the event horizon always acts as a sink for matter and information. From an inside view, the matter and information disappear into the singularity, which is the boundary of spacetime. The argument is based on the nature of the processes at work plus a careful delineation of the relevant causal domains; in order to confirm this model and determine details of the outcome, detailed calculations of the expectation value of the stress-energy-momentum tensor are needed to determine back reaction effects.

The trace-free version of the Einstein Gravitational equations, essentially equivalent to unimodular gravity, can solve the troubling issue of the huge discrepancy between quantum field theory estimates of the vacuum energy density and the astronomically observed value of the cosmological constant. However it has been suggested that this proposal cannot work because it prevents the inflaton potential energy from driving inflation. It is shown here that that concern is unjustified: inflation proceeds as usual if we adopt the trace free gravitational equations.

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.

We propose a thermodynamically motivated measure of gravitational entropy based on the Bel-Robinson tensor, which has a natural interpretation as the effective super-energy-momentum tensor of free gravitational fields. The specific form of this measure differs depending on whether the gravitational field is Coulomb-like or wave-like, and reduces to the Bekenstein-Hawking value when integrated over the interior of a Schwarzschild black hole. For scalar perturbations of a Robertson-Walker geometry we find that the entropy goes like the Hubble weighted anisotropy of the gravitational field, and therefore increases as structure formation occurs. This is in keeping with our expectations for the behaviour of gravitational entropy in cosmology, and provides a thermodynamically motivated arrow of time for cosmological solutions of Einstein's field equations. It is also in keeping with Penrose's Weyl curvature hypothesis.

The relation between the expanding universe and local vacuum solutions, such as that for the Solar System, is crucially mediated by Birkhoff's theorem. Here we consider how that relation works, and give generalizations of Birkhoff's theorem when there are geometric and matter and perturbations. The issue of to what degree dark matter might influence the solar system emerges as a significant question.

A violation of the reciprocity relation, which induces a violation of the distance duality relation, reflects itself in a change in the normalisation of the cosmic microwave spectrum in such a way that its spectrum is grey. We show that existing observational constraints imply that the reciprocity relation cannot be violated by more than 0.01% between decoupling and today. We compare this effect to other sources of violation of the distance duality relations which induce spectral distortion of the cosmic microwave background spectrum.

Light from `point sources' such as supernovae is observed with a beam width of order of the sources' size - typically less than 1 AU. Such a beam probes matter and curvature distributions that are very different from coarse-grained representations in N-body simulations or perturbation theory, which are smoothed on scales much larger than 1 AU. The beam typically travels through unclustered dark matter and hydrogen with a mean density much less than the cosmic mean, and through dark matter halos and hydrogen clouds. Using N-body simulations, as well as a Press-Schechter approach, we quantify the density probability distribution as a function of beam width and show that, even for Gpc-length beams of 500 kpc diameter, most lines of sight are significantly under-dense. From this we argue that modelling the probability distribution for AU-diameter beams is absolutely critical. Standard analyses predict a huge variance for such tiny beam sizes, and nonlinear corrections appear to be non-trivial. It is not even clear whether under-dense regions lead to dimming or brightening of sources, owing to the uncertainty in modelling the expansion rate which we show is the dominant contribution. By considering different reasonable approximations which yield very different cosmologies we argue that modelling ultra-narrow beams accurately remains a critical problem for precision cosmology. This could appear as a discordance between angular diameter and luminosity distances when comparing SN observations to BAO or CMB distances.

Birkhoff's theorem for spherically symmetric vacuum spacetimes is a key theorem in studying local systems in general relativity theory. However realistic local systems are only approximately spherically symmetric and only approximately vacuum. In a previous paper, we showed the theorem remains approximately true in an approximately spherically symmetric vacuum space time. In this paper we prove the converse case: the theorem remains approximately true in a spherically symmetric, approximately vacuum space time.

Structure occurs over a vast range of scales in the universe. Our large-scale cosmological models are coarse-grained representations of what exists, which have much less structure than there really is. An important problem for cosmology is determining the influence the small-scale structure in the universe has on its large-scale dynamics and observations. Is there a significant, general relativistic, backreaction effect from averaging over structure? One issue is whether the process of smoothing over structure can contribute to an acceleration term and so alter the apparent value of the cosmological constant. If this is not the case, are there other aspects of concordance cosmology that are affected by backreaction effects? Despite much progress, this 'averaging problem' is still unanswered, but it cannot be ignored in an era of precision cosmology.

A surprising exact result for the Einstein Field Equations is that if pressure-free matter is moving in a shear-free way, then it must be either expansion-free or rotation-free. It has been suggested this result is also true for any barotropic perfect fluid, but a proof has remained elusive. We consider the case of barotropic perfect fluid solutions linearized about a Robertson-Walker geometry, and prove that the result remains true except for the case of a specific highly non-linear equation of state. We argue that this equation of state is non-physical, and hence the result is true in the linearized case for all physically realistic barotropic perfect fluids. This result, which is not true in Newtonian cosmology, demonstrates that the linearized solutions, believed to result in standard local Newtonian theory, do not always give the usual behaviour of Newtonian solutions.

The Goldberg-Sachs theorem is an exact result on shear-free null geodesics in a vacuum spacetime. It is compared and contrasted with an exact result for pressure-free matter: shear-free flows cannot both expand and rotate. In both cases, the shear-free condition restricts the way distant matter can influence the local gravitational field. This leads to intriguing discontinuities in the relation of the General Relativity solutions to Newtonian solutions in the timelike case, and of the full theory to the linearised theory in the null case. It is a pleasure to dedicate this paper to Josh Goldberg.

This is an editorial note to accompany reprinting as a Golden Oldie in the Journal of General Relativity and Gravitation of the famous paper by Brandon Carter on the anthropic principle in cosmology \cite{Car74}. This paper was presented at IAU Symposium No. 63, entitled Confrontation of cosmological theories with observational data, in 1973.

The quantum field theoretic prediction for the vacuum energy density leads to a value for the effective cosmological constant that is incorrect by between 60 to 120 orders of magnitude. We review an old proposal of replacing Einstein's Field Equations by their trace-free part (the Trace-Free Einstein Equations), together with an independent assumption of energy--momentum conservation by matter fields. While this does not solve the fundamental issue of why the cosmological constant has the value that is observed cosmologically, it is indeed a viable theory that resolves the problem of the discrepancy between the vacuum energy density and the observed value of the cosmological constant. However, one has to check that, as well as preserving the standard cosmological equations, this does not destroy other predictions, such as the junction conditions that underlie the use of standard stellar models. We confirm that no problems arise here: hence, the Trace-Free Einstein Equations are indeed viable for cosmological and astrophysical applications.

This article looks at how inhomogeneous spacetime models may be significant for cosmology. First it looks at how the averaging process may affect large scale dynamics, with backreaction effects leading to effective contributions to the averaged energy-momentum tensor. Secondly it considers how local inhomogeneities may affect cosmological observations in cosmology, possibly significantly affecting the concordance model parameters. Thirdly it presents the possibility that the universe is spatially inhomogeneous on Hubble scales, with a violation of the Copernican principle leading to an apparent acceleration of the universe. This could perhaps even remove the need for the postulate of dark energy.

In order to understand how locally static configurations around gravitationally bound bodies can be embedded in an expanding universe, we investigate the solutions of general relativity describing a space-time whose spatial sections have the topology of a 3-sphere with two identical masses at the poles. We show that Israel junction conditions imply that two spherically symmetric static regions around the masses cannot be glued together. If one is interested in an exterior solution, this prevents the geometry around the masses to be of the Schwarzschild type and leads to the introduction of a cosmological constant. The study of the extension of the Kottler space-time shows that there exists a non-static solution consisting of two static regions surrounding the masses that match a Kantowski-Sachs expanding region on the cosmological horizon. The comparison with a Swiss-Cheese construction is also discussed.

We extend Birkhoff's theorem for almost LRS-II vacuum spacetimes to show that the rigidity of spherical vacuum solutions of Einstein's field equations continues even in the perturbed scenario.

In some well-known scenarios of open-universe eternal inflation, developed by Vilenkin and co-workers, a large number of universes nucleate and thermalize within the eternally inflating mega-universe. According to the proposal, each universe nucleates at a point, and therefore the boundary of the nucleated universe is a space-like surface nearly coincident with the future light cone emanating from the point of nucleation, all points of which have the same proper-time. This leads the authors to conclude that at the proper-time t = t_{nuc} at which any such nucleation occurs, an infinite open universe comes into existence. We point out that this is due entirely to the supposition of the nucleation occurring at a single point, which in light of quantum cosmology seems difficult to support. Even an infinitesimal space-like length at the moment of nucleation gives a rather different result -- the boundary of the nucleating universe evolves in proper-time and becomes infinite only in an infinite time. The alleged infinity is never attained at any finite time.

The causal limit usually considered in cosmology is the particle horizon, delimiting the possibilities of causal connection in the expanding universe. However it is not a realistic indicator of the effective local limits of important interactions in spacetime. We consider here the matter horizon for the Solar System, that is,the comoving region which has contributed matter to our local physical environment. This lies inside the effective domain of dependence, which (assuming the universe is dominated by dark matter along with baryonic matter and vacuum-energy-like dark energy) consists of those regions that have had a significant active physical influence on this environment through effects such as matter accretion and acoustic waves. It is not determined by the velocity of light c, but by the flow of matter perturbations along their world lines and associated gravitational effects. We emphasize how small a region the perturbations which became our Galaxy occupied, relative to the observable universe -- even relative to the smallest-scale perturbations detectable in the cosmic microwave background radiation. Finally, looking to the future of our cosmic domain, we suggest simple dynamical criteria for determining the present domain of influence and the future matter horizon. The former is the radial distance at which our local region is just now separating from the cosmic expansion. The latter represents the limits of growth of the matter horizon in the far future.

The thoughts expressed in this article are based on remarks made by J\"urgen Ehlers at the Albert-Einstein-Institut, Golm, Germany in July 2007. The main objective of this article is to demonstrate, in terms of plausible order-of-magnitude estimates for geometrical scalars, the relevance of spatial curvature in realistic models of the Universe that describe the dynamics of structure formation since the epoch of matter-radiation decoupling. We introduce these estimates with a commentary on the use of a quasi-Newtonian metric form in this context.

I consider the relation of explanations for the observed data to testability in the following contexts: observational and experimental detection of dark matter; observational and experimental detection of dark energy or a cosmological constant $\Lambda$; observational or experimental testing of the multiverse proposal to explain a small non-zero value of $\Lambda$; and observational testing of the possibility of large scale spatial inhomogeneity with zero $\Lambda$.

We present the time drift of the cosmological redshift in a general spherically symmetric spacetime. We demonstrate that its observation would allow us to test the Copernican principle and so determine if our universe is radially inhomogeneous, an important issue in our understanding of dark energy. In particular, when combined with distance data, this extra observable allows one to fully reconstruct the geometry of a spacetime describing a spherically symmetric under-dense region around us, purely from background observations.

A usual causal requirement on a viable theory of matter is that the speed of sound be at most the speed of light. In view of various recent papers querying this limit, the question is revisited here. We point to various issues confronting theories that violate the usual constraint.

The various requirements on a consistent varying speed of light (`VSL') theory are surveyed, giving a short check-list of issues that should be satisfactorily handled by such theories.

After a survey of the present state of cosmological theory and observations, this article discusses a series of major themes underlying the relation of philosophy to cosmology. These are: A: The uniqueness of the universe; B: The large scale of the universe in space and time; C: The unbound energies in the early universe; D: Explaining the universe -- the question of origins; E: The universe as the background for existence; F: The explicit philosophical basis; G: The Anthropic question: fine tuning for life; H: The possible existence of multiverses; I: The natures of existence. Each of these themes is explored and related to a series of Theses that set out the major issues confronting cosmology in relation to philosophy.

Susskind claims in his recent book The Cosmic Landscape that evidence for the existence and nature of `pocket universes' in a multiverse would be available in the detailed nature of the Cosmic Blackbody Background Radiation that constantly bathes all parts of our observable universe. I point out that acceptance of the complex chain of argument involved does not imply possible experimental verification of multiverses at the present time. Rather this claim relates only to theoretically possible observations in the very far future of the universe.

The idea of a multiverse -- an ensemble of universes or universe domains -- has received increasing attention in cosmology, both as the outcome of the originating process that generated our own universe, and as an explanation for why our universe appears to be fine-tuned for life and consciousness. Here we carefully consider how multiverses should be defined, stressing the distinction between the collection of all possible universes, and ensembles of really existing universes, which are essential for an anthropic argument. We show that such realised multiverses are by no means unique, and in general require the existence of a well-defined and physically motivated distribution function on the space of all possible universes. Furthermore, a proper measure on these spaces is also needed, so that probabilities can be calculated. We then discuss several other major physical and philosophical problems which arise in the context of ensembles of universes, including the emergence and causal effectiveness of self-consciousness, realized infinities, and fine- tuning, or the apparent need for very special initial conditions for our universe -- whether they or generalized generic primordial conditions are more fundamental. Then we briefly summarise scenarios like chaotic inflation, which suggest how ensembles of universe domains may be generated, and point out that the regularities which must underlie any systematic description of truly disjoint multiverses must imply some kind of common generating mechanism. Finally, we discuss the issue of testability, which underlies the question of whether multiverse proposals are really scientific propositions.

A large-scale smoothed-out model of the universe ignores small-scale inhomogeneities, but the averaged effects of those inhomogeneities may alter both observational and dynamical relations at the larger scale. This article discusses these effects, and comments briefly on the relation to gravitational entropy.

Closed, singularity-free, inflationary cosmological models have recently been studied in the context of general relativity. Despite their appeal, these so called emergent models suffer from a number of limitations. These include the fact that they rely on an initial Einstein static state to describe the past eternal phase of the universe. Given the instability of such a state within the context of general relativity, this amounts to a very severe fine tuning. Also in order to be able to study the dynamics of the universe within the context of general relativity, they set the initial conditions for the universe in the classical phase. Here we study the existence and stability of such models in the context of Loop Quantum Cosmology and show that both these limitations can be partially remedied, once semi-classical effects are taken into account. An important consequence of these effects is to give rise to a static solution (not present in GR), which dynamically is a centre equilibrium point and located in the more natural semi-classical regime. This allows the construction of emergent models in which the universe oscillates indefinitely about such an initial static state. We construct an explicit emergent model of this type, in which a non-singular past eternal oscillating universe enters a phase where the symmetry of the oscillations is broken, leading to an emergent inflationary epoch, while satisfying all observational and semi-classical constraints. We also discuss emergent models in which the universe possesses both early- and late-time accelerating phases.

We provide a realisation of a singularity-free inflationary universe in the form of a simple cosmological model dominated at early times by a single minimally coupled scalar field with a physically based potential. The universe starts asymptotically from an initial Einstein static state, which may be large enough to avoid the quantum gravity regime. It enters an expanding phase that leads to inflation followed by reheating and a standard hot Big Bang evolution. We discuss the basic characteristics of this Emergent model and show that none is at odds with current observations.

Observations indicate that the universe is effectively flat, but they do not rule out a closed universe. The role of positive curvature is negligible at late times, but can be crucial in the early universe. In particular, positive curvature allows for cosmologies that originate as Einstein static universes, and then inflate and later reheat to a hot big bang era. These cosmologies have no singularity, no "beginning of time", and no horizon problem. If the initial radius is chosen to be above the Planck scale, then they also have no quantum gravity era, and are described by classical general relativity throughout their history.

The idea of a multiverse -- an ensemble of universes -- has received increasing attention in cosmology, both as the outcome of the originating process that generated our own universe, and as an explanation for why our universe appears to be fine-tuned for life and consciousness. Here we carefully consider how multiverses should be defined, stressing the distinction between the collection of all possible universes, and ensembles of really existing universes that are essential for an anthropic argument. We show that such realised multiverses are by no means unique. A proper measure on the space of all really existing universes or universe domains is needed, so that probabilities can be calculated, and major problems arise in terms of realised infinities. As an illustration we examine these issues in the case of the set of Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) universes. Then we briefly summarise scenarios like chaotic inflation, which suggest how ensembles of universe domains may be generated, and point out that the regularities which must underlie any systematic description of truly disjoint multiverses must imply some kind of common generating mechanism. Finally, we discuss the issue of testability, which underlies the question of whether multiverse proposals are really scientific propositions.

Theories proposing a varying speed of light have recently been widely promoted under the claim that they offer an alternative way of solving the standard cosmological problems. Recent observational hints that the fine structure constant may have varied during over cosmological scales also has given impetus to these models. In theoretical physics the speed of light, $c$, is hidden in almost all equations but with different facets that we try to distinguish. Together with a reminder on scalar-tensor theories of gravity, this sheds some light on these proposed varying speed of light theories.

We show using covariant techniques that the Einstein static universe containing a perfect fluid is always neutrally stable against small inhomogeneous vector and tensor perturbations and neutrally stable against adiabatic scalar density inhomogeneities so long as c_{s}^2>1/5, and unstable otherwise. We also show that the stability is not significantly changed by the presence of a self-interacting scalar field source, but we find that spatially homogeneous Bianchi type IX modes destabilise an Einstein static universe. The implications of these results for the initial state of the universe and its pre-inflationary evolution are also discussed.

We present a general framework for analyzing spatially inhomogeneous cosmological dynamics. It employs Hubble-normalized scale-invariant variables which are defined within the orthonormal frame formalism, and leads to the formulation of Einstein's field equations with a perfect fluid matter source as an autonomous system of evolution equations and constraints. This framework incorporates spatially homogeneous dynamics in a natural way as a special case, thereby placing earlier work on spatially homogeneous cosmology in a broader context, and allows us to draw on experience gained in that field using dynamical systems methods. One of our goals is to provide a precise formulation of the approach to the spacelike initial singularity in cosmological models, described heuristically by Belinski\v{\i}, Khalatnikov and Lifshitz. Specifically, we construct an invariant set which we conjecture forms the local past attractor for the evolution equations. We anticipate that this new formulation will provide the basis for proving rigorous theorems concerning the asymptotic behavior of spatially inhomogeneous cosmological models.

Inter alia, the high precision WMAP data on Cosmic Microwave Background Radiation marginally indicate that the universe has positively curved (and hence spherical) spatial sections. In this paper, we take this data seriously and consider some of the consequences for the background dynamics. In particular, we show that this implies a limit to the number of e-foldings that could have taken place in the inflationary epoch; however this limit is consistent with some inflationary models that solve all the usual cosmological problems and are consistent with standard structure formation theory.

We study the peculiar motion of non-relativistic matter in a fully covariant way. The exact nonlinear equations are derived and then applied to the case of pressure-free matter, moving relatively to a quasi-Newtonian Eulerian frame. Our two-frame formalism facilitates the study of the nonlinear kinematics of the matter, as the latter decouples from the background expansion and starts to ``turn around'' and collapse. Applied to second perturbative order, our equations provide a fully covariant formulation of the Zeldovich approximation, which by construction addresses the mildly nonlinear regime of structure formation. Employing a dynamical system approach, we show that, just like in the Newtonian case, the relativistic treatment also predicts that pancakes are the natural end-structures for any generic overdensity.

We give a careful general relativistic and (1+3)-covariant analysis of cosmological peculiar velocities induced by matter density perturbations in the presence of a cosmological constant. In our quasi-Newtonian approach, constraint equations arise to maintain zero shear of the non-comoving fundamental worldlines which define a Newtonian-like frame, and these lead to the (1+3)-covariant dynamical equations, including a generalized Poisson-type equation. We investigate the relation between peculiar velocity and peculiar acceleration, finding the conditions under which they are aligned. In this case we find (1+3)-covariant relativistic generalizations of well-known Newtonian results.

We study strong gravitational lensing due to a Schwarzschild black hole. Apart from the primary and the secondary images we find a sequence of images on both sides of the optic axis; we call them {\em relativistic images}. These images are formed due to large bending of light near r = 3M (the closest distance of approach r_o is greater than 3M). The sources of the entire universe are mapped in the vicinity of the black hole by these images. For the case of the Galactic supermassive ``black hole'' they are formed at about 17 microarcseconds from the optic axis. The relativistic images are not resolved among themselves, but they are resolved from the primary and secondary images. However the relativistic images are very much demagnified unless the observer, lens and source are very highly aligned. Due to this and some other difficulties the observation of these images does not seem to be feasible in near future. However, it would be a great success of the general theory of relativity in a strong gravitational field if they ever were observed and it would also give an upper bound, r_o = 3.21 M, to the compactness of the lens, which would support the black hole interpretation of the lensing object.

This is the second of a series of papers extending the 1+3 covariant and gauge invariant treatment of kinetic theory to an examination of Cosmic Microwave Background temperature anisotropies arising from inhomogeneities in the early universe. The first paper dealt with algebraic issues. Here we derive the mode form of the integrated Boltzmann equations, first, giving a covariant version of the standard derivation using the mode recursion relations, second, demonstrating the link to the multipole divergence equations and finally various analytic ways of solving the resulting equations are discussed. A general integral form of solution is obtained for the equations with Thomson scattering. The covariant Friedmann-Lemaitre multipole form of the transport equations are found using the covariant and gauge-invariant generalization of the Peebles and Yu expansion in Thompson scattering time. The dispersion relations and damping scale are then obtained from the covariant approach. The equations are integrated to give the covariant and gauge-invariant equivalent of the canonical scalar sourced anisotropies. We carry out a simple treatment of the matter dominated free-streaming projection, slow decoupling, and tight-coupling cases, with the aim both giving a unified transparent derivation of this range of results and clarifying the connection between the more usual approaches (for example that of Hu and Sugiyama) and the treatment for scalar perturbations (for example the treatment of Challinor and Lasenby).

We consider a recent generalisation by Bousso of an earlier holography proposal by Fischler and Susskind. We demonstrate that in general inhomogeneous universes such a proposal would involve extremely complicated - possibly fractal - light sheets. Furthermore, in general such a light sheet cannot be known a priori on the basis of theory and moreover, the evolution of the universe makes it clear that in general such bounds cannot remain invariant under time reversal and will change with epoch. We propose a modified version of this proposal in which the light sheets end on the boundary of the past, and hence avoid contact with the caustics. In this way the resulting light sheets and projections can be made much simpler. We discuss the question of operational definability of these sheets within the context of both proposals and conclude that in both cases the theoretical existence of such sheets must be clearly distinguished from their complexity and the difficulty of their construction in practice. This puts into perspective the likely practical difficulties one would face in applying the holographic principle to the real cosmos. These issues may also be of relevance in debates regarding the applications of the holographic principle to other settings such as string theory.

This is the first of a series of papers extending a covariant and gauge invariant (CGI) treatment of kinetic theory in curved space-times to a treatment of Cosmic Background Radiation (CBR) temperature anisotropies arising from inhomogeneities in the early universe. This paper deals with algebraic issues, both generically and in the context of models linearised about RW geometries. The approach represents radiation anisotropies by PSTF tensors. The Angular correlation functions for the mode coefficients are found in terms of these quantities, following the Wilson-Silk approach, but derived and dealt with in CGI form. The covariant multipole and mode-expanded angular correlation functions are related to the usual treatments in the literature. The CGI mode expansion is related to the coordinate approach by linking the Legendre functions to the PSTF representation, using a covariant addition theorem for the tensors to generate the Legendre Polynomial recursion relation. This paper lays the foundation for further papers in the series, which use this formalism in a CGI approach to develop solutions of the Boltzmann and Liouville equations for the CBR before and after decoupling, thus providing a unified CGI derivation of the variety of approaches to CBR anisotropies in the current literature.

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.

Spatially homogeneous cosmological models with a positive cosmological constant are investigated, using dynamical systems methods. We focus on the future evolution of these models. In particular, we address the question whether there are models within this class that are de Sitter-like in the future, but are tilted.

We develop a new approach to local nonlinear effects in cosmic microwave background anisotropies, and discuss the qualitative features of these effects. New couplings of the baryonic velocity to radiation multipoles are found, arising from nonlinear Thomson scattering effects. We also find a new nonlinear shear effect on small angular scales. The full set of evolution and constraint equations is derived, including the nonlinear generalizations of the radiation multipole hierarchy, and of the dynamics of multi-fluids. These equations govern radiation anisotropies in any inhomogeneous spacetime, but their main application is to second-order effects in a universe that is close to the Friedmann models. Qualitative analysis is given here, and quantitative calculations are taken up in further papers.

We employ the extended 1+3 orthonormal frame formalism for fluid spacetime geometries $({\cal M}, {\bf g}, {\bf u})$, which contains the Bianchi field equations for the Weyl curvature, to derive a 44-D evolution system of first-order symmetric hyperbolic form for a set of geometrically defined dynamical field variables. Describing the matter source fields phenomenologically in terms of a barotropic perfect fluid, the propagation velocities $v$ (with respect to matter-comoving observers that Fermi-propagate their spatial reference frames) of disturbances in the matter and the gravitational field, represented as wavefronts by the characteristic 3-surfaces of the system, are obtained. In particular, the Weyl curvature is found to account for two (non-Lorentz-invariant) Coulomb-like characteristic eigenfields propagating with $v = 0$ and four transverse characteristic eigenfields propagating with $|v| = 1$, which are well known, and four (non-Lorentz-invariant) longitudinal characteristic eigenfields propagating with $|v| = \sfrac{1}{2}$. The implications of this result are discussed in some detail and a parallel is drawn to the propagation of irregularities in the matter distribution. In a worked example, we specialise the equations to cosmological models in locally rotationally symmetric class II and include the constraints into the set of causally propagating dynamical variables.

Exact dynamical equations for a generic dust matter source field in a cosmological context are formulated with respect to a non-comoving Newtonian-like timelike reference congruence and investigated for internal consistency. On the basis of a lapse function $N$ (the relativistic acceleration scalar potential) which evolves along the reference congruence according to $\dot{N} = \alpha \Theta N$ ($\alpha = {const}$), we find that consistency of the quasi-Newtonian dynamical equations is not attained at the first derivative level. We then proceed to show that a self-consistent set can be obtained by linearising the dynamical equations about a (non-comoving) FLRW background. In this case, on properly accounting for the first-order momentum density relating to the non-relativistic peculiar motion of the matter, additional source terms arise in the evolution and constraint equations describing small-amplitude energy density fluctuations that do not appear in similar gravitational instability scenarios in the standard literature.

We consider compensated spherical lens models and the caustic surfaces they create in the past light cone. Examination of cusp and crossover angles associated with particular source and lens redshifts gives explicit lensing models that confirm previous claims that area distances can differ by substantial factors from angular diameter distances even when averaged over large angular scales. `Shrinking' in apparent sizes occurs, typically by a factor of 3 for a single spherical lens, on the scale of the cusp caused by the lens; summing over many lenses will still leave a residual effect.

We consider the changes which occur in cosmological distances due to the combined effects of some null geodesics passing through low-density regions while others pass through lensing-induced caustics. This combination of effects increases observed areas corresponding to a given solid angle even when averaged over large angular scales, through the additive effect of increases on all scales, but particularly on micro-angular scales; however angular sizes will not be significantly effected on large angular scales (when caustics occur, area distances and angular-diameter distances no longer coincide). We compare our results with other works on lensing, which claim there is no such effect, and explain why the effect will indeed occur in the (realistic) situation where caustics due to lensing are significant. Whether or not the effect is significant for number counts depends on the associated angular scales and on the distribution of inhomogeneities in the universe. It could also possibly affect the spectrum of CBR anisotropies on small angular scales, indeed caustics can induce a non-Gaussian signature into the CMB at small scales and lead to stronger mixing of anisotropies than occurs in weak lensing.

This paper and the others in the series challenge the standard model of the effects of gravitational lensing on observations at large distances. We show that due to the cumulative effect of lensing, areas corresponding to an observed solid angle can be quite different than would be estimated from the corresponding Friedmann-Lema\^{\i}tre model, even when averaged over large angular scales. This paper concentrates on the specific example of spherically symmetric but spatially inhomogeneous dust universes, the Lema\^{\i}tre-Tolman-Bondi models, and shows that radial lensing significantly distorts the area distance-redshift and density-redshift relations in these exact solutions compared with the standard ones for Friedmann-Lema\^{\i}tre models. Thus inhomogeneity may introduce significant errors into distance estimates based on the standard FL relations, even after all-sky averaging. In addition a useful new gauge choice is presented for these models, solving the problem of locating the past null cone exactly.

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.

A recent paper (Martinez--Gonzalez & Sanz 1995) showed that if the universe is homogeneous but anisotropic, then the small quadrupole anisotropy in the cosmic microwave background radiation implies that the spacetime anisotropy is very small. We point out that more general results may be established, without assuming a priori homogeneity. We have proved that small anisotropies in the microwave background imply that the universe is almost Friedmann--Robertson--Walker. Furthermore, the quadrupole and octopole place direct and explicit limits on the degree of anisotropy and inhomogeneity, as measured by the shear, vorticity, Weyl tensor and density gradients. In the presence of inhomogeneity, it is only possible to set a much weaker limit on the shear than that given by Martinez--Gonzalez & Sanz.

The dynamics of perfect fluid spacetime geometries which exhibit {\em Local Rotational Symmetry} (LRS) are reformulated in the language of a $1+\,3$ "threading" decomposition of the spacetime manifold, where covariant fluid and curvature variables are used. This approach presents a neat alternative to the orthonormal frame formalism. The dynamical equations reduce to a set of differential relations between purely scalar quantities. The consistency conditions are worked out in a transparent way. We discuss their various subcases in detail and focus in particular on models with higher symmetries within the class of expanding spatially inhomogeneous LRS models, via a consideration of functional dependencies between the dynamical variables.

We consider directly the equations by which matter imposes anisotropies on freely propagating background radiation, leading to a new way of using anisotropy measurements to limit the deviations of the Universe from a Friedmann-Robertson-Walker (FRW) geometry. This approach is complementary to the usual Sachs-Wolfe approach: the limits obtained are not as detailed, but they are more model-independent. We also give new results about combined matter-radiation perturbations in an almost-FRW universe, and a new exact solution of the linearised equations.

In this paper we examine the change in the estimated spatial power spectra at decoupling due to the effects of our clumpy universe which modify observational distances. We find that scales at decoupling can be significantly underestmated in our approximation of neglecting the shear of the ray bundle. We compare our results with other work on lensing and speculate on the implications for structure formation. In particular we examine a proposal to use the position of the first Doppler peak to determine $\Omega$, and find that shrinking will modify the esimated curvature, so that it must be included to obtain an accurate estimate of $\Omega$. Finally we consider future applications and improvements of our results.

In an interesting recent paper on the growth of inhomogeneity through the effect of gravity [1], Bertschinger and Hamilton derive equations for the electric and magnetic parts of the Weyl tensor for cold dust for both General Relativity and Newtonian theory. Their conclusion is that both in General Relativity and in Newtonian theory, in general the magnetic part of the Weyl tensor does not vanish, implying that the Lagrangian evolution of the fluid is not local. We show here that the `Newtonian' theory discussed by them is in fact not Newtonian theory {\it per se}, but rather a plausible relativistic generalisation of Newtonian theory. Newtonian cosmology itself is highly non-local irrespective of the behaviour of the magnetic part of the Weyl tensor; in this respect the Bertschinger-Hamilton generalisation is a better theory.

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).