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 analytically compute time domain gravitational waveforms produced in the final stages of extreme mass ratio inspirals of non-spinning compact objects into supermassive nearly extremal Kerr black holes. Conformal symmetry relates all corotating equatorial orbits in the geodesic approximation to circular orbits through complex conformal transformations. We use this to obtain the time domain Teukolsky perturbations for generic equatorial corotating plunges in closed form. The resulting gravitational waveforms consist of an intermediate polynomial ringdown phase in which the decay rate depends on the impact parameters, followed by an exponential quasi-normal mode decay. The waveform amplitude exhibits critical behavior when the orbital angular momentum tends to a minimal value determined by the innermost stable circular orbit. We show that either near-critical or large angular momentum leads to a significant extension of the LISA observable volume of gravitational wave sources of this kind.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

The usual theory of inflation breaks down in eternal inflation. We derive a dual description of eternal inflation in terms of a deformed Euclidean CFT located at the threshold of eternal inflation. The partition function gives the amplitude of different geometries of the threshold surface in the no-boundary state. Its local and global behavior in dual toy models shows that the amplitude is low for surfaces which are not nearly conformal to the round three-sphere and essentially zero for surfaces with negative curvature. Based on this we conjecture that the exit from eternal inflation does not produce an infinite fractal-like multiverse, but is finite and reasonably smooth.

Structure at the horizon scale of black holes would give rise to echoes of the gravitational wave signal associated with the post-merger ringdown phase in binary coalescences. We study the waveform of echoes in static and stationary, traversable wormholes in which perturbations are governed by a symmetric effective potential. We argue that echoes are dominated by the wormhole quasinormal frequency nearest to the fundamental black hole frequency that controls the primary signal. We put forward an accurate method to construct the echoes waveform(s) from the primary signal and the quasinormal frequencies of the wormhole, which we characterize. We illustrate this in the static Damour-Solodukhin wormhole and in a new, rotating generalization that approximates a Kerr black hole outside the throat. Rotation gives rise to a potential with an intermediate plateau region that breaks the degeneracy of the quasinormal frequencies. Rotation also leads to late-time instabilities which, however, fade away for small angular momentum.

We consider the quantum dynamics of gravitational collapse in a model in which the wave function spreads out over a large ensemble of geometries as envisioned in the fuzzball proposal. We show that the probabilities of coarse-grained observables are highly peaked around the classical black hole values. By contrast, probabilities for finer-grained observables probing the neighbourhood of collapsed objects are more broadly distributed and no notion of `averaging' applies to them. This implies that the formation of fuzzballs gives rise to distinct observational signatures that are more significant than has hitherto been thought and may be tested against observations in the near future. We also predict a novel kind of gravitational wave burst associated with the spreading of the wave function in gravitational collapse.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

One proposal for dS/CFT is that the Hartle-Hawking (HH) wave function in the large volume limit is equal to the partition function of a Euclidean CFT deformed by various operators. All saddle points defining the semiclassical HH wave function in cosmology have a representation in which their interior geometry is part of a Euclidean AdS domain wall with complex matter fields. We compute the wave functions of scalar and tensor perturbations around homogeneous isotropic complex saddle points, turning on single scalar field matter only. We compare their predictions for the spectra of CMB perturbations with those of a different dS/CFT proposal based on the analytic continuation of inflationary universes to real asymptotically AdS domain walls. We find the predictions of both bulk calculations agree to first order in the slow roll parameters, but there is a difference at higher order which, we argue, is a signature of the HH state of the fluctuations.

We apply the principles of quantum mechanics and quantum cosmology to predict probabilities for our local observations of a universe undergoing false vacuum eternal inflation. At a sufficiently fine-grained level, histories of the universe describe a mosaic of bubble universes separated by inflationary regions. We show that predictions for local observations can be obtained directly from sets of much coarser grained histories which only follow a single bubble. These coarse-grained histories contain neither information about our unobservable location nor about the unobservable large-scale structure outside our own bubble. Applied to a landscape of false vacua in the no-boundary state we predict our local universe emerged from the dominant decay channel of the lowest energy false vacuum. We compare and contrast this framework for prediction based on quantum cosmology with traditional approaches to the measure problem in cosmology.

We numerically generate the six-dimensional landscape of D3-brane inflation and identify patches of eternal inflation near sufficiently flat inflection points of the potential. We show that reasonable measures that select patches of eternal inflation in the landscape yield sharp predictions for the spectral properties of primordial perturbations on observable scales. These include a scalar tilt of .936, a running of the scalar tilt -.00103, undetectably small tensors and non-Gaussianity, and no observable spatial curvature. Our results explicitly demonstrate that precision cosmology probes the combination of the statistical properties of the string landscape and the measure implied by the universe's quantum state.

The quantum state of the universe combined with the structure of the landscape potential implies a prior that specifies predictions for observations. We compute the prior for CMB related observables given by the no-boundary wave function (NBWF) in a landscape model that includes a range of inflationary patches representative of relatively simple single-field models. In this landscape the NBWF predicts our classical cosmological background emerges from a region of eternal inflation associated with a plateau-like potential. The spectra of primordial fluctuations on observable scales are characteristic of concave potentials, in excellent agreement with the Planck data. By contrast, alternative theories of initial conditions that strongly favor inflation at high values of the potential are disfavored by observations in this landscape.

We show that anthropic selection emerges inevitably in the general framework for prediction in quantum cosmology. There the predictions of anthropic reasoning depend on the prior implied by the universe's quantum state. To illustrate this we compute the probabilities specified by the no-boundary wave function for our observations at the present time of the values of Lambda and Q in an inflationary landscape model in which both quantities vary. Within the anthropic range of values the no-boundary state yields an approximately flat distribution on Lambda and strongly favors small values of Q. This restores Weinberg's successful prediction of Lambda.

The AdS/CFT correspondence allows us to map a dynamical cosmology to a dual quantum field theory living on the boundary of spacetime. Specifically, we study a five-dimensional model cosmology in type IIB supergravity, where the dual theory is an unstable deformation of $\N=4$ supersymmetric SU(N) gauge theory on $\Rbar\times S^3$. A one-loop computation shows that the coupling governing the instability is asymptotically free, so quantum corrections cannot turn the potential around. The big crunch singularity in the bulk occurs when a boundary scalar field runs to infinity, in finite time. Consistent quantum evolution requires that we impose boundary conditions at infinite scalar field, i.e. a self-adjoint extension of the system. We find that quantum spreading of the homogeneous mode of the boundary scalar leads to a natural UV cutoff in particle production as the wavefunction for the homogeneous mode bounces back from infinity. However a perturbative calculation indicates that despite this, the logarithmic running of the boundary coupling governing the instability generally leads to significant particle production across the bounce. This prevents the wave packet of the homogeneous boundary scalar to return close to its initial form. Translating back to the bulk theory, we conclude that a quantum transition from a big crunch to a big bang is an improbable outcome of cosmological evolution in this class of five-dimensional models.

We derive the arrows of time of our universe that follow from the no-boundary theory of its quantum state (NBWF) in a minisuperspace model. Arrows of time are viewed four-dimensionally as properties of the four-dimensional Lorentzian histories of the universe. Probabilities for these histories are predicted by the NBWF. For histories with a regular `bounce' at a minimum radius we find that fluctuations are small at the bounce and grow in the direction of expansion on either side. For recollapsing classical histories with big bang and big crunch singularities we find that the fluctuations are small near one singularity and grow through the expansion and recontraction to the other singularity. The arrow of time defined by the growth in fluctuations thus points in one direction over the whole of a recollapsing spacetime but is bidirectional in a bouncing spacetime. We argue that the electromagnetic, thermodynamic, and psychological arrows of time are aligned with the fluctuation arrow. The implications of a bidirectional arrow of time for causality are discussed.

We consider landscape models that admit several regions where the conditions for eternal inflation hold. It is shown that one can use the no-boundary wave function to calculate small departures from homogeneity within our past light cone despite the possibility of much larger fluctuations on super horizon scales. The dominant contribution comes from the history exiting eternal inflation at the lowest value of the potential. In a class of landscape models this predicts a tensor to scalar ratio of about 10%. In this way the no-boundary wave function defines a measure for the prediction of local cosmological observations.

Junctions on cosmic string loops give rise to the proliferation of sharp kinks. We study the effect of this proliferation on the gravitational wave (GW) signals emitted from string networks with junctions, assuming a scaling solution. We calculate the rate of occurrence and the distribution in amplitude of the GW bursts emitted at cusps and kinks in the frequency bands of LIGO and LISA as a function of the string tension, the number of sharp kinks on loops with junctions and the fraction of loops in the cosmological network which have junctions. Combining our results with current observational constraints, we find that pulsar data rule out a significant number of kinks on loops for strings with tensions G\mu > 10^{-12}. By contrast, for smaller tensions current observations allow for a large number of kinks on loops. If this is the case, the incoherent superposition of small bursts emitted at kink-kink encounters leads to an enhanced GW background that hides the strong individual bursts from kinks and cusps.

The no-boundary wave function (NBWF) specifies a measure for prediction in cosmology that selects inflationary histories and remains well behaved for spatially large or infinite universes. This paper explores the predictions of the NBWF for linear scalar fluctuations about homogeneous and isotropic backgrounds in models with a single scalar field moving in a quadratic potential. We treat both the space-time geometry of the universe and the observers inhabiting it quantum mechanically. We evaluate top-down probabilities for local observations that are conditioned on the NBWF and on part of our data as observers of the universe. For models where the most probable histories do not have a regime of eternal inflation, the NBWF predicts homogeneity on large scales, a specific non-Gaussian spectrum of observable fluctuations, and a small amount of inflation in our past. By contrast, for models where the dominant histories have a regime of eternal inflation, the NBWF predicts significant inhomogeneity on scales much larger than the present horizon, a Gaussian spectrum of observable fluctuations, and a long period of inflation in our past. The absence or presence of local non-Gaussianity therefore provides information about the global structure of the universe, assuming the NBWF.

We consider the no-boundary proposal for homogeneous isotropic closed universes with a cosmological constant and a scalar field with a quadratic potential. In the semi-classical limit, it predicts classical behavior at late times if the initial scalar field is more than a certain minimum. If the classical late time histories are extended back, they may be singular or bounce at a finite radius. The no-boundary proposal provides a probability measure on the classical solutions which selects inflationary histories but is heavily biased towards small amounts of inflation. This would not be compatible with observations. However we argue that the probability for a homogeneous universe should be multiplied by exp(3N) where N is the number of e-foldings of slow roll inflation to obtain the probability for what we observe in our past light cone. This volume weighting is similar to that in eternal inflation. In a landscape potential, it would predict that the universe would have a large amount of inflation and that it would start in an approximately de Sitter state near a saddle-point of the potential. The universe would then have always been in the semi-classical regime.

We analyze the origin of the quasiclassical realm from the no-boundary proposal for the universe's quantum state in a class of minisuperspace models. The models assume homogeneous, isotropic, closed spacetime geometries, a single scalar field moving in a quadratic potential, and a fundamental cosmological constant. The allowed classical histories and their probabilities are calculated to leading semiclassical order. We find that for the most realistic range of parameters analyzed a minimum amount of scalar field is required, if there is any at all, in order for the universe to behave classically at late times. If the classical late time histories are extended back, they may be singular or bounce at a finite radius. The ensemble of classical histories is time symmetric although individual histories are generally not. The no-boundary proposal selects inflationary histories, but the measure on the classical solutions it provides is heavily biased towards small amounts of inflation. However, the probability for a large number of efoldings is enhanced by the volume factor needed to obtain the probability for what we observe in our past light cone, given our present age. Our results emphasize that it is the quantum state of the universe that determines whether or not it exhibits a quasiclassical realm and what histories are possible or probable within that realm.

The AdS/CFT correspondence is used to describe five-dimensional cosmology with a big crunch singularity in terms of super-Yang-Mills theory on R times S^3 deformed by a potential which is unbounded below. Classically, a Higgs field in the dual theory rolls to infinity in finite time. But since the S^3 is finite, the unstable mode spreads quantum mechanically and the singularity is resolved when self-adjoint boundary conditions are imposed at infinity. Asymptotic freedom of the coupling governing the instability gives us computational control and the quantum spreading provides a UV cutoff on particle creation. The bulk interpretation of our result is a quantum transition from a big crunch to a big bang. An intriguing consequence of the near scale-invariance of the dual theory is that a nearly scale-invariant spectrum of stress-energy perturbations is automatically generated in the boundary theory. We comment on implications for more realistic cosmologies.

We put forward a framework for cosmology that combines the string landscape with no boundary initial conditions. In this framework, amplitudes for alternative histories for the universe are calculated with final boundary conditions only. This leads to a top down approach to cosmology, in which the histories of the universe depend on the precise question asked. We study the observational consequences of no boundary initial conditions on the landscape, and outline a scheme to test the theory. This is illustrated in a simple model landscape that admits several alternative inflationary histories for the universe. Only a few of the possible vacua in the landscape will be populated. We also discuss in what respect the top down approach differs from other approaches to cosmology in the string landscape, like eternal inflation.

This paper investigates Starobinsky's model of inflation driven by the trace anomaly of conformally coupled matter fields. This model does not suffer from the problem of contrived initial conditions that occurs in most models of inflation driven by a scalar field. The universe can be nucleated semi-classically by a cosmological instanton that is much larger than the Planck scale provided there are sufficiently many matter fields. There are two cosmological instantons: the four sphere and a new ``double bubble'' solution. This paper considers a universe nucleated by the four sphere. The AdS/CFT correspondence is used to calculate the correlation function for scalar and tensor metric perturbations during the ensuing de Sitter phase. The analytic structure of the scalar and tensor propagators is discussed in detail. Observational constraints on the model are discussed. Quantum loops of matter fields are shown to strongly suppress short scale metric perturbations, which implies that short distance modifications of gravity would probably not be observable in the cosmic microwave background. This is probably true for any model of inflation provided there are sufficiently many matter fields. This point is illustrated by a comparison of anomaly driven inflation in four dimensions and in a Randall-Sundrum brane-world model.

We study a Randall-Sundrum cosmological scenario consisting of a domain wall in anti-de Sitter space with a strongly coupled large $N$ conformal field theory living on the wall. The AdS/CFT correspondence allows a fully quantum mechanical treatment of this CFT, in contrast with the usual treatment of matter fields in inflationary cosmology. The conformal anomaly of the CFT provides an effective tension which leads to a de Sitter geometry for the domain wall. This is the analogue of Starobinsky's four dimensional model of anomaly driven inflation. Studying this model in a Euclidean setting gives a natural choice of boundary conditions at the horizon. We calculate the graviton correlator using the Hartle-Hawking ``No Boundary'' proposal and analytically continue to Lorentzian signature. We find that the CFT strongly suppresses metric perturbations on all but the largest angular scales. This is true independently of how the de Sitter geometry arises, i.e., it is also true for four dimensional Einstein gravity. Since generic matter would be expected to behave like a CFT on small scales, our results suggest that tensor perturbations on small scales are far smaller than predicted by all previous calculations, which have neglected the effects of matter on tensor perturbations.

We compute the spectrum of primordial gravitational wave perturbations in open de Sitter spacetime. The background spacetime is taken to be the continuation of an O(5) symmetric instanton saddle point of the Euclidean no boundary path integral. The two-point tensor fluctuations are computed directly from the Euclidean path integral. The Euclidean correlator is then analytically continued into the Lorentzian region where it describes the quantum mechanical vacuum fluctuations of the graviton field. Unlike the results of earlier work, the correlator is shown to be unique and well behaved in the infrared. We show that the infrared divergence found in previous calculations is due to the contribution of a discrete gauge mode inadvertently included in the spectrum.

We compute the tensor CMB anisotropy power spectrum for singular and non-singular instantons describing the beginning of an open universe according to the Euclidean no boundary proposal. Singular instantons occur generically, whereas non-singular instantons require more contrived scalar field potentials. For the latter, we consider potentials in which a sharp feature, either negative or positive, is added to a gently sloping potential. In the first case one finds a nearly divergent contribution to the low multipole CMB anisotropy, in conflict with the COBE observations. In the second case the divergence is weaker, but matching the low multipoles forces the added feature to be large and narrow. For singular instantons, there is a better match to the observations, without any such contrivance. The distinction between singular and nonsingular instantons disappears in the limit as the universe becomes flat, but is still observable for densities as high as 0.7 of the critical density.

We perform a first principles computation of the spectrum of gravity waves produced in open inflationary universes. The background spacetime is taken to be the continuation of an instanton saddle point of the Euclidean no boundary path integral. The two-point tensor correlator is computed directly from the path integral and is shown to be unique and well behaved in the infrared. We discuss the tensor contribution to the cosmic microwave background anisotropy and show how it may provide an observational discriminant between different types of primordial instantons.