Primordial black holes (PBHs) may result from high peaks in a random field of cosmological perturbations. In single field inflationary models, such perturbations can be seeded as the inflaton overshoots a small barrier on its way down the potential. PBHs are then produced through two distinct mechanisms, during the radiation era. The first one is the familiar collapse of large adiabatic overdensities. The second one is the collapse induced by relic bubbles where the inflaton field is trapped in a false vacuum, due to large backward fluctuations which prevented horizon sized regions from overshooting the barrier. We consider (numerically and analytically) the effect of non-Gaussianities on the threshold for overdensities to collapse into a PBH. Since typical high peaks have some dispersion in their shape or profile, we also consider the effect of such dispersion on the corresponding threshold for collapse. With these results we estimate the most likely channel for PBH production as a function of the non-Gaussianity parameter $f_{\rm NL}$. We also compare the threshold for collapse coming from the perturbative versus the non perturbative template for the non-Gaussianity arising in this model. We show that i) for $f_{\rm NL}\gtrsim 3.5$, the population of PBH coming from false vacuum regions dominates over that which comes from the collapse of large adiabatic overdensities, ii) the non-perturbative template of the non-Gaussianities is important to get accurate results. iii) the effect of the dispersion is small in determining the threshold for the compaction function, although it can be appreciable in determining the threshold amplitude for the curvature perturbation at low $f_{\rm NL}$. We also confirm that the volume averaged compaction function provides a very accurate universal estimator for the threshold.

The rate of merger events observed by LIGO/Virgo can be used in order to probe the fraction $f$ of dark mater in the form of Primordial Black Holes. Here, we consider the merger rate of PBH binaries, accounting for the effect of cosmological perturbations on their initial eccentricity $e$. The torque on the binaries receives significant contributions from a wide range of scales, that goes from the size of the horizon at the time when the binary forms, down to the co-moving size of the binary. In scenarios where PBH are formed from adiabatic perturbations, it is natural to expect an enhancement of the power spectrum $P_\Phi$ at small scales, where it is poorly constrained observationally. The effect can then be quite significant. For instance, a nearly flat spectrum with amplitude $P_\Phi \gtrsim 10^{-7}$ on scales smaller than $\sim 10 Mpc^{-1}$ gives a contribution $\langle j^2 \rangle \sim 10^3 P_\Phi$, where $j= (1-e^2)^{1/2}$ is the dimensionless angular momentum parameter of the binaries. This contribution can dominate over tidal torques from neighboring PBHs for any value of $f$. Current constraints allow for a power spectrum as large as $P_\Phi \sim 10^{-5}$ at the intermediate scales $10^3-10^5 Mpc^{-1}$, comparable to the co-moving size of the binaries at the time of formation. In particular, this can relax current bounds on the PBH abundance based on the observed LIGO/Virgo merger rate, allowing for a fraction $f\sim 10\%$ of dark matter in PBH of mass $\sim 30 M_\odot$. We investigate the differential merger rate $\Delta\Gamma(m_1,m_2)$, as a function of the masses of the binary components, and the corresponding ``universality" coefficient $\alpha = -(m_1+m_2)^2 \partial^2 \ln \Delta\Gamma/\partial m_1\partial m_2$, concluding that merger rates may provide valuable information on the spectrum of primordial cosmological perturbations at currently uncharted lengthscales.

In the context of transient constant-roll inflation near a local maximum, we derive the non-perturbative field redefinition that relates a Gaussian random field with the true non-Gaussian curvature perturbation. Our analysis shows the emergence of a new critical amplitude $\zeta_*$, corresponding to perturbations that prevent the inflaton from overshooting the local maximum, thus becoming trapped in the false minimum of the potential. For potentials with a mild curvature at the local maximum (and thus small non-Gaussianity), we recover the known perturbative field redefinition. We apply these results to the formation of primordial black holes, and discuss the cases for which $\zeta_*$ is smaller or of the same order than the critical value for collapse of spherically symmetric overdensities. In the latter case, we present a simple potential for which the power spectrum needs an amplitude 10 times smaller that in the Gaussian case for producing a sizeable amount of primordial black holes.

The production rate of primordial black holes is often calculated by considering a nearly Gaussian distribution of cosmological perturbations, and assuming that black holes will form in regions where the amplitude of such perturbations exceeds a certain threshold. A threshold $\zeta_{\rm th}$ for the curvature perturbation is somewhat inappropriate for this purpose, because it depends significantly on environmental effects, not essential to the local dynamics. By contrast, a threshold $\delta_{\rm th}$ for the density perturbation at horizon crossing seems to provide a more robust criterion. On the other hand, the density perturbation is known to be bounded above by a maximum limit $\delta_{\rm max}$, and given that $\delta_{\rm th}$ is comparable to $\delta_{\rm max}$, the density perturbation will be far from Gaussian near or above the threshold. In this paper, we provide a new plausible estimate for the primordial black hole abundance based on peak theory. In our approach, we assume that the curvature perturbation is given as a random Gaussian field with the power spectrum characterized by a single scale, while an optimized criterion for PBH formation is imposed, based on the locally averaged density perturbation. Both variables are related by the full nonlinear expression derived in the long-wavelength approximation of general relativity. We do not introduce a window function, and the scale of the inhomogeneity is introduced as a random variable in the peak theory. We find that the mass spectrum is shifted to larger mass scales by one order of magnitude or so, compared to a conventional calculation. The abundance of PBHs becomes significantly larger than the conventional one, by many orders of magnitude, mainly due to the optimized criterion for PBH formation and the removal of the suppresion associated with a window function.

It is well known that, in single clock inflation, the curvature perturbation $\zeta$ is constant in time on superhorizon scales. In the standard bulk description this follows quite simply from the local conservation of the energy momentum tensor in the bulk. On the other hand, in a holographic description, the constancy of the curvature perturbation must be related to the properties of the RG flow in the boundary theory. Here, we show that, in single clock holographic inflation, the time independence of correlators of $\zeta$ follows from the cut-off independence of correlators of the energy momentum tensor in the boundary theory, and from the so-called consistency relations for vertex functions with a soft leg.

Vacuum bubbles may nucleate and expand during the inflationary epoch in the early universe. After inflation ends, the bubbles quickly dissipate their kinetic energy; they come to rest with respect to the Hubble flow and eventually form black holes. The fate of the bubble itself depends on the resulting black hole mass. If the mass is smaller than a certain critical value, the bubble collapses to a singularity. Otherwise, the bubble interior inflates, forming a baby universe, which is connected to the exterior FRW region by a wormhole. A similar black hole formation mechanism operates for spherical domain walls nucleating during inflation. As an illustrative example, we studied the black hole mass spectrum in the domain wall scenario, assuming that domain walls interact with matter only gravitationally. Our results indicate that, depending on the model parameters, black holes produced in this scenario can have significant astrophysical effects and can even serve as dark matter or as seeds for supermassive black holes. The mechanism of black hole formation described in this paper is very generic and has important implications for the global structure of the universe. Baby universes inside super-critical black holes inflate eternally and nucleate bubbles of all vacua allowed by the underlying particle physics. The resulting multiverse has a very non-trivial spacetime structure, with a multitude of eternally inflating regions connected by wormholes. If a black hole population with the predicted mass spectrum is discovered, it could be regarded as evidence for inflation and for the existence of a multiverse.

We consider the superpotential formalism to describe the evolution of scalar fields during inflation, generalizing it to include the case with non-canonical kinetic terms. We provide a characterization of the attractor behaviour of the background evolution in terms of first and second slow-roll parameters (which need not be small). We find that the superpotential is useful in justifying the separate universe approximation from the gradient expansion, and also in computing the spectra of primordial perturbations around attractor solutions in the $\delta N$ formalism. As an application, we consider a class of models where the background trajectories for the inflaton fields are derived from a product separable superpotential. In the perspective of the holographic inflation scenario, such models are dual to a deformed CFT boundary theory, with $D$ mutually uncorrelated deformation operators. We compute the bulk power spectra of primordial adiabatic and entropy cosmological perturbations, and show that the results agree with the ones obtained by using conformal perturbation theory in the dual picture.

Many theories of the early universe predict the existence of a multiverse where bubbles continuously nucleate giving rise to observers in their interior. In this paper, we point out that topological defects of several dimensionalities will also be produced in de Sitter like regions of the multiverse. In particular, defects could be spontaneously nucleated in our parent vacuum. We study the evolution of these defects as they collide with and propagate inside of our bubble. We estimate the present distribution of defects in the observable part of the universe. The expected number of such nearby defects turns out to be quite small, even for the highest nucleation rate. We also study collisions of strings and domain walls with our bubble in our past light cone. We obtain simulated full-sky maps of the loci of such collisions, and find their angular size distribution. Similarly to what happens in the case of bubble collisions, the prospect of detecting any collisions of our bubble with ambient defects is greatly enhanced in the case where the cosmological constant of our parent vacuum is much higher than the vacuum energy density during inflation in our bubble.

We initiate the study of multi-field inflation using holography. Bulk light scalar fields correspond to nearly marginal operators in the boundary theory and the dual quantum field theory is a deformation of a CFT by such operators. We compute the power spectra of adiabatic and entropy perturbations in a simple model and find that the adiabatic curvature perturbation is not conserved in the presence of entropy perturbations but becomes conserved when the entropy perturbations are set to zero or the model is effectively a single scalar model, in agreement with expectations from cosmological perturbation theory.

In a holographic description of inflation, cosmological time evolution in the bulk is expected to correspond to the renomalization group (RG) flow in a dual boundary theory. Here, we analyze this expectation by computing the correlation functions of the curvature perturbation $\zeta$ holographically. For this purpose, we use a deformed conformal field theory at the boundary, with a single deformation operator. In standard single field models of inflation, $\zeta$ is known to be conserved at large scales under very general conditions. However, we find that this is not generically the case in the dual description. The requirement that higher correlators of $\zeta$ should be conserved severely restricts the possibilities for the RG flow. With such restriction, the power spectrum $P_\zeta$ must follow an exact power law, at least within the regime of validity of conformal perturbation theory.

We argue that theories with ghosts may have a long lived vacuum state even if all interactions are Lorentz preserving. In space-time dimension D = 2, we consider the tree level decay rate of the vacuum into ghosts and ordinary particles mediated by non-derivative interactions, showing that this is finite and logarithmically growing in time. For D > 2, the decay rate is divergent unless we assume that the interaction between ordinary matter and the ghost sector is soft in the UV, so that it can be described in terms of non-local form factors rather than point-like vertices. We provide an example of a nonlocal gravitational-strength interaction between the two sectors, which appears to satisfy all observational constraints.

We study the primordial perturbations generated during a stage of single-field inflation in Einstein-aether theories. Quantum fluctuations of the inflaton and aether fields seed long wavelength adiabatic and isocurvature scalar perturbations, as well as transverse vector perturbations. Geometrically, the isocurvature mode is the potential for the velocity field of the aether with respect to matter. For a certain range of parameters, this mode may lead to a sizable random velocity of the aether within the observable universe. The adiabatic mode corresponds to curvature perturbations of co-moving slices (where matter is at rest). In contrast with the standard case, it has a non-vanishing anisotropic stress on large scales. Scalar and vector perturbations may leave significant imprints on the cosmic microwave background. We calculate their primordial spectra, analyze their contributions to the temperature anisotropies, and formulate some of the phenomenological constraints that follow from observations. These may be used to further tighten the existing limits on the parameters for this class of theories. The results for the scalar sector also apply to the extension of Horava gravity recently proposed by Blas, Pujolas and Sibiriakov.

Probabilities in the multiverse can be calculated by assuming that we are typical representatives in a given reference class. But is this class well defined? What should be included in the ensemble in which we are supposed to be typical? There is a widespread belief that this question is inherently vague, and that there are various possible choices for the types of reference objects which should be counted in. Here we argue that the ``ideal'' reference class (for the purpose of making predictions) can be defined unambiguously in a rather precise way, as the set of all observers with identical information content. When the observers in a given class perform an experiment, the class branches into subclasses who learn different information from the outcome of that experiment. The probabilities for the different outcomes are defined as the relative numbers of observers in each subclass. For practical purposes, wider reference classes can be used, where we trace over all information which is uncorrelated to the outcome of the experiment, or whose correlation with it is beyond our current understanding. We argue that, once we have gathered all practically available evidence, the optimal strategy for making predictions is to consider ourselves typical in any reference class we belong to, unless we have evidence to the contrary. In the latter case, the class must be correspondingly narrowed.

In theories with a hidden ghost sector that couples to visible matter through gravity only, empty space can decay into ghosts and ordinary matter by graviton exchange. Perturbatively, such processes can be very slow provided that the gravity sector violates Lorentz invariance above some cut-off scale. Here, we investigate non-perturbative decay processes involving ghosts, such as the spontaneous creation of self-gravitating lumps of ghost matter, as well as pairs of Bondi dipoles (i.e., lumps of ghost matter chasing after positive energy objects). We find the corresponding instantons and calculate their Euclidean action. In some cases, the instantons induce topology change or have negative Euclidean action. To shed some light on the meaning of such peculiarities, we also consider the nucleation of concentrical domain walls of ordinary and ghost matter, where the Euclidean calculation can be compared with the canonical (Lorentzian) description of tunneling. We conclude that non-perturbative ghost nucleation processes can be safely suppressed in phenomenological scenarios.

Inflationary cosmology leads to the picture of a "multiverse," involving an infinite number of (spatially infinite) post-inflationary thermalized regions, called pocket universes. In the context of theories with many vacua, such as the landscape of string theory, the effective constants of Nature are randomized by quantum processes during inflation. We discuss an analytic estimate for the volume distribution of the constants within each pocket universe. This is based on the conjecture that the field distribution is approximately ergodic in the diffusion regime, when the dynamics of the fields is dominated by quantum fluctuations (rather than by the classical drift). We then propose a method for determining the relative abundances of different types of pocket universes. Both ingredients are combined into an expression for the distribution of the constants in pocket universes of all types.

We discuss probability distributions for the cosmological constant Lambda and the amplitude of primordial density fluctuations Q in models where they both are anthropic variables. With mild assumptions about the prior probabilities, the distribution P(Lambda,Q) factorizes into two independent distributions for the variables Q and $y \propto \Lambda/Q^3$. The distribution for y is largely model-independent and is in a good agreement with the observed value of y. The form of P(Q) depends on the origin of density perturbations. If the perturbations are due to quantum fluctuations of the inflaton, then P(Q) tends to have an exponential dependence on Q, due to the fact that in such models Q is correlated with the amount of inflationary expansion. For simple models with a power-law potential, P(Q) is peaked at very small values of Q, far smaller than the observed value of 10^{-5}. This problem does not arise in curvaton-type models, where the inflationary expansion factor is not correlated with Q.

There have been some recent claims that brane-worlds of co-dimension two in a 6D bulk with compact extra dimensions may lead to self-tuning of the effective 4D cosmological constant. Here we show that if a phase transition occurs on a flat brane, so as to change its tension, then the brane will not remain flat. In other words, there is really no self-tuning in such models, which can in fact be understood in four-dimensional terms and are therefore subject to Weinberg's no-go theorem.

In the context of models where the dark energy density $\rD$ is a random variable, anthropic selection effects may explain both the "old" cosmological constant problem and the "time coincidence". We argue that this type of solution to both cosmological constant problems entails a number of definite predictions, which can be checked against upcoming observations. In particular, in models where the dark energy density is a discrete variable, or where it is a continuous variable due to the potential energy of a single scalar field, the anthropic approach predicts that the dark energy equation of state is $p_D=-\rho_D$ with a very high accuracy. It is also predicted that the dark energy density is greater than the currently favored value $\Omega_D\approx 0.7$. Another prediction, which may be testable with an improved understanding of galactic properties, is that the conditions for civilizations to emerge arise mostly in galaxies completing their formation at low redshift, $z\approx 1$. Finally, there is a prediction which may not be easy to test observationally: our part of the universe is going to recollapse eventually. However, the simplest models predict that it will take more than a trillion years of accelerated expansion before this happens.

A broad class of dark energy models, which have been proposed in attempts at solving the cosmological constant problems, predict a late time variation of the equation of state with redshift. The variation occurs as a scalar field picks up speed on its way to negative values of the potential. The negative potential energy eventually turns the expansion into contraction and the local universe undergoes a big crunch. In this paper we show that cross-correlations of the CMB anisotropy and matter distribution, in combination with other cosmological data, can be used to forecast the imminence of such cosmic doomsday.

Nucleation of branes by a four-form field has recently been considered in string motivated scenarios for the neutralization of the cosmological constant. An interesting question in this context is whether the nucleation of stacks of coincident branes is possible, and if so, at what rate does it proceed. Feng et al. have suggested that, at high ambient de Sitter temperature, the rate may be strongly enhanced, due to large degeneracy factors associated with the number of light species living on the worldsheet. This might facilitate the quick relaxation from a large effective cosmological constant down to the observed value. Here, we analyse this possibility in some detail. In four dimensions, and after the moduli are stabilized, branes interact via repulsive long range forces. Because of that, the Coleman-de Luccia (CdL) instanton for coincident brane nucleation may not exist, unless there is some short range interaction which keeps the branes together. If the CdL instanton exists, we find that the degeneracy factor depends only mildly on the ambient de Sitter temperature, and does not switch off even in the case of tunneling from flat space. This would result in catastrophic decay of the present vacuum. If, on the contrary, the CdL instanton does not exist, coindident brane nucleation may still proceed through a "static" instanton, representing pair creation of critical bubbles -- a process somewhat analogous to thermal activation in flat space. In that case, the branes may stick together due to thermal symmetry restoration, and the pair creation rate depends exponentially on the ambient de Sitter temperature, switching off sharply as the temperature approaches zero. Such static instanton may be well suited for the "saltatory" relaxation scenario proposed by Feng et al.

We explore the possibility that the dark energy is due to a potential of a scalar field and that the magnitude and the slope of this potential in our part of the universe are largely determined by anthropic selection effects. We find that, in some models, the most probable values of the slope are very small, implying that the dark energy density stays constant to very high accuracy throughout cosmological evolution. In other models, however, the most probable values of the slope are such that the slow roll condition is only marginally satisfied, leading to a re-collapse of the local universe on a time-scale comparable to the lifetime of the sun. In the latter case, the effective equation of state varies appreciably with the redshift, leading to a number of testable predictions.

We give a simple interpretation of the recent solutions for cosmologies with a transient accelerating phase obtained from compactification in hyperbolic manifolds, or from S-brane solutions of string/M-theory. In the four-dimensional picture, these solutions correspond to bouncing the radion field off its exponential potential. Acceleration occurs at the turning point, when the radion stops and the potential energy momentarily dominates. The virtues and limitations of these approaches become quite transparent in this interpretation.

We consider a class of 5-D brane-world solutions with a power-law warp factor $a(y)\propto y^{q}$, and bulk dilaton with profile $\phi \propto \ln y$, where $y$ is the proper distance in the extra dimension. This class includes the Heterotic M-theory brane-world of Refs. \cite{ovrut,ekpy} and the Randall-Sundrum (RS) model as a limiting case. In general, there are two moduli fields $y_{\pm}$, corresponding to the "positions" of two branes (which live at the fixed points of an orbifold compactification). Classically, the moduli are massless, due to a scaling symmetry of the action. However, in the absence of supersymmetry, they develop an effective potential at one loop. Local terms proportional to $K_\pm^4$, where $K_\pm=q/y_\pm$ is the local curvature scale at the location of the corresponding brane, are needed in order to remove the divergences in the effective potential. Such terms break the scaling symmetry and hence they may act as stabilizers for the moduli. When the branes are very close to each other, the effective potential induced by massless bulk fields behaves like $V\sim d^{-4}$, where $d$ is the separation between branes. When the branes are widely separated, the potentials for each one of the moduli generically develop a "Coleman-Weinberg"-type behaviour of the form $a^4(y_\pm)K_\pm^4 \ln(K_\pm/\mu_\pm)$, where $\mu_\pm$ are renormalization scales. In the RS case, the bulk geometry is $AdS$ and $K_\pm$ are equal to a constant, independent of the position of the branes, so these terms do not contribute to the mass of the moduli. However, for generic warp factor, they provide a simple stabilization mechanism. For $q\gtrsim 10$, the observed hierarchy can be naturally generated by this potential, giving the lightest modulus a mass of order $m_- \lesssim TeV$.

Bucher [Bucher2001] has recently proposed an interesting brane-world cosmological scenario where the ``Big Bang'' hypersurface is the locus of collision of two vacuum bubbles which nucleate in a five dimensional flat space. This gives rise to an open universe, where the curvature can be very small provided that $d/R_0$ is sufficiently large. Here, d is the distance between bubbles and $R_0$ is their size at the time of nucleation. Quantum fluctuations develop on the bubbles as they expand towards each other, and these in turn imprint cosmological perturbations on the initial hypersurface. We present a simple formalism for calculating the spectrum of such perturbations and their subsequent evolution. We conclude that, unfortunately, the spectrum is very tilted, with spectral index $n_s=3$. The amplitude of fluctuations at horizon crossing is given by $<(\delta \rho/\rho)^2> \sim (R_0/d)^2 S_E^{-1} k^2$, where $S_E\gg 1$ is the Euclidean action of the instanton describing the nucleation of a bubble and k is the wavenumber in units of the curvature scale. The spectrum peaks on the smallest possible relevant scale, whose wave-number is given by $k\sim d/R_0$. We comment on the possible extension of our formalism to more general situations where a Big Bang is ignited through the collision of 4D extended objects.

A generic prediction of inflation is that the thermalized region we inhabit is spatially infinite. Thus, it contains an infinite number of regions of the same size as our observable universe, which we shall denote as $\O$-regions. We argue that the number of possible histories which may take place inside of an $\O$-region, from the time of recombination up to the present time, is finite. Hence, there are an infinite number of $\O$-regions with identical histories up to the present, but which need not be identical in the future. Moreover, all histories which are not forbidden by conservation laws will occur in a finite fraction of all $\O$-regions. The ensemble of $\O$-regions is reminiscent of the ensemble of universes in the many-world picture of quantum mechanics. An important difference, however, is that other $\O$-regions are unquestionably real.

We critically review several recent approaches to solving the two cosmological constant problems. The "old" problem is the discrepancy between the observed value of $\Lambda$ and the large values suggested by particle physics models. The second problem is the "time coincidence" between the epoch of galaxy formation $t_G$ and the epoch of $\Lambda$-domination $t_\L$. It is conceivable that the "old" problem can be resolved by fundamental physics alone, but we argue that in order to explain the "time coincidence" we must account for anthropic selection effects. Our main focus here is on the discrete-$\Lambda$ models in which $\Lambda$ can change through nucleation of branes. We consider the cosmology of this type of models in the context of inflation and discuss the observational constraints on the model parameters. The issue of multiple brane nucleation raised by Feng {\it et. al.} is discussed in some detail. We also review continuous-$\L$ models in which the role of the cosmological constant is played by a slowly varying potential of a scalar field. We find that both continuous and discrete models can in principle solve both cosmological constant problems, although the required values of the parameters do not appear very natural. M-theory-motivated brane models, in which the brane tension is determined by the brane coupling to the four-form field, do not seem to be viable, except perhaps in a very tight corner of the parameter space. Finally, we point out that the time coincidence can also be explained in models where $\Lambda$ is fixed, but the primordial density contrast $Q=\delta\rho/\rho$ is treated as a random variable.

Some of the parameters we call ``constants of Nature'' may in fact be variables related to the local values of some dynamical fields. During inflation, these variables are randomized by quantum fluctuations. In cases when the variable in question (call it $\chi$) takes values in a continuous range, all thermalized regions in the universe are statistically equivalent, and a gauge invariant procedure for calculating the probability distribution for $\chi$ is known. This is the so-called ``spherical cutoff method''. In order to find the probability distribution for $\chi$ it suffices to consider a large spherical patch in a single thermalized region. Here, we generalize this method to the case when the range of $\chi$ is discontinuous and there are several different types of thermalized region. We first formulate a set of requirements that any such generalization should satisfy, and then introduce a prescription that meets all the requirements. We finally apply this prescription to calculate the relative probability for different bubble universes in the open inflation scenario.

We show that in brane-world scenarios with warped extra dimensions, the Casimir force due to bulk matter fields may be sufficient to stabilize the radion field $\phi$. In particular, we calculate one loop effective potential for $\phi$ induced by bulk gravitons and other possible massless bulk fields in the Randall-Sundrum background. This potential has a local extremum, which can be a maximum or a minimum depending on the detailed bulk matter content. If the parameters of the background are chosen so that the hierarchy problem is solved geometrically, then the radion mass induced by Casimir corrections is hierarchycally smaller than the $TeV$. Hence, in this important case, we must invoke an alternative mechanism (classical or nonperturbative) which gives the radion a sizable mass, to make it compatible with observations.

We discuss the large-scale structure of the universe in inflationary cosmology and the implications that it may have for the long-term future of civilizations. Although each civilization is doomed to perish, it may be possible to transmit its accumulated knowledge to future civilizations. We consider several scenarios of this sort. If the cosmological constant is positive, it eventually dominates the universe and bubbles of inflationary phase begin to nucleate at a constant rate. Thermalized regions inside these inflating bubbles will give rise to new galaxies and civilizations. It is possible in principle to send a message to one of them. It might even be possible to send a device whose purpose is to recreate an approximation of the original civilization in the new region. However, the message or device will almost certainly be intercepted by black holes, which nucleate at a much higher rate than inflating bubbles. Formation of new inflating regions can also be triggered by gravitational collapse, but again the probability is low, and the number of attempts required for a positive outcome is enormous. The probability can be higher if the energy scale of inflation is closer to the Planck scale, but a high energy scale produces a tight bound on the amount of information that can be transmitted. One can try to avoid quantum tunneling altogether, but this requires a violation of quantum inequalities which constrain the magnitude of negative energy densities. However, the limits of validity of quantum inequalities are not clear, and future research may show that the required violation is in fact possible. Therein lies the hope for the future of civilizations.

We discuss models in which the smallness of the effective vacuum energy density $\rho_\L$ and the coincidence of the time of its dominance $t_\L$ with the epoch of galaxy formation $t_G$ are due to anthropic selection effects. In such models, the probability distribution for $\rho_\L$ is a product of an {\it a priori} distribution ${\cal P}_*(\rho_\L)$ and of the number density of galaxies at a given $\rho_\L$ (which is proportional to the number of observers who will detect that value of $\rho_\L$). To determine ${\cal P}_*$, we consider inflationary models in which the role of the vacuum energy is played by a slowly-varying potential of some scalar field. We show that the resulting distribution depends on the shape of the potential and generally has a non-trivial dependence on $\rho_\L$, even in the narrow anthropically allowed range. This is contrary to Weinberg's earlier conjecture that the {\it a priori} distribution should be nearly flat in the range of interest. We calculate the (final) probability distributions for $\rho_\L$ and for $t_G/t_\L$ in simple models with power-law potentials. For some of these models, the agreement with the observationally suggested values of $\rho_\L$ is better than with a flat {\it a priori} distribution. We also discuss quantum-cosmological approach in which $\rho_\L$ takes different values in different disconnected universes and argue that Weinberg's conjecture is not valid in this case as well. Finally, we extend our analysis to models of quintessence, with similar conclusions.

An inflating brane-world can be created from ``nothing'' together with its Anti-de Sitter (AdS) bulk. The resulting space-time has compact spatial sections bounded by the brane. During inflation, the continuum of KK modes is separated from the massless zero mode by the gap $m=(3/2) H$, where $H$ is the Hubble rate. We consider the analog of the Nariai solution and argue that it describes pair production of ``Black cigars'' attached to the inflating brane. In the case when the size of the instantons is much larger than the AdS radius, the 5-dimensional action agrees with the 4-dimensional one. Hence, the 5D and 4D gravitational entropies are the same in this limit. We also consider thermal instantons with an AdS black hole in the bulk. These may be interpreted as describing the creation of a hot universe from nothing, or the production of AdS black holes in the vicinity of a pre-existing inflating brane-world. The Lorentzian evolution of the brane-world after creation is briefly discussed. An additional "integration constant" in the Friedmann equation -accompanying a term which dilutes like radiation- describes the tidal force in the fifth direction and arises from the mass of a spherical object inside the bulk. This could be a 5-dimensional black hole or a "parallel" brane-world of negative tension concentrical with our brane-world. In the case of thermal solutions, and in the spirit of the $AdS/CFT$ correspondence, one may attribute the additional term to thermal radiation in the boundary theory. Then, for temperatures well below the AdS scale, the entropy of this radiation agrees with the entropy of the black hole in the AdS bulk.

We discuss the weak gravitational field created by isolated matter sources in the Randall-Sundrum brane-world. In the case of two branes of opposite tension, linearized Brans-Dicke (BD) gravity is recovered on either wall, with different BD parameters. On the wall with positive tension the BD parameter is larger than 3000 provided that the separation between walls is larger than 4 times the AdS radius. For the wall of negative tension, the BD parameter is always negative but greater than -3/2. In either case, shadow matter from the other wall gravitates upon us. For equal Newtonian mass, light deflection from shadow matter is 25 % weaker than from ordinary matter. Hence, the effective mass of a clustered object containing shadow dark matter would be underestimated if naively measured through its lensing effect. For the case of a single wall of positive tension, Einstein gravity is recovered on the wall to leading order, and if the source is stationary the field stays localized near the wall. We calculate the leading Kaluza-Klein corrections to the linearized gravitational field of a non-relativistic spherical object and find that the metric is different from the Schwarzschild solution at large distances. We believe that our linearized solution corresponds to the field far from the horizon after gravitational collapse of matter on the brane.

We explore a model in which the cosmological constant $\Lambda$ and the density contrast at the time of recombination $\sigma_{rec}$ are random variables, whose range and {\it a priori} probabilities are determined by the laws of physics. (Such models arise naturally in the framework of inflationary cosmology.) Based on the assumption that we are typical observers, we show that the order of magnitude coincidence among the three timescales: the time of galaxy formation, the time when the cosmological constant starts to dominate the cosmic energy density, and the present age of the universe, finds a natural explanation. We also discuss the probability distribution for $\sigma_{rec}$, and find that it is peaked near the observationally suggested values, for a wide class of {\it a priori} distributions.

In the context of open inflation, we calculate the probability distribution for the density parameter $\Omega$. A large class of two field models of open inflation do not lead to infinite open universes, but to an ensemble of inflating islands of finite size, or ``quasi-open'' universes, where the density parameter takes a range of values. Assuming we are typical observers, the models make definite predictions for the value $\Omega$ we are most likely to observe. When compared with observations, these predictions can be used to constrain the parameters of the models. We also argue that obsevers should not be surprised to find themselves living at the time when curvature is about to dominate.

We extend the theory of cosmological perturbations to the case when the ``matter'' Lagrangian is an arbitrary function of the scalar field and its first derivatives. In particular, this extension provides a unified description of known cases such as the usual scalar field and the hydrodynamical perfect fluid. In addition, it applies to the recently proposed k-inflation, which is driven by non-minimal kinetic terms in the Lagrangian. The spectrum of quantum fluctuations for slow-roll and power law k-inflation is calculated. We find, for instance, that the usual ``consistency relation'' between the tensor spectral index and the relative amplitude of scalar and tensor perturbations is modified. Thus, at least in principle, k-inflation is phenomenologically distinguishable from standard inflation.

The spectrum of cosmological perturbations in the context of the one-bubble open inflation model is discussed, taking into account fluctuations of the metric. We find that, quite generically,thin wall single field models have no supercurvature modes. However, single field models with supercurvature modes do exist. In these models the density parameter $\Omega$ becomes a random variable taking a range of values inside of each bubble. We also show that the model dependence of the continuous spectrum for both scalar and tensor-type perturbations is small as long as the kinetic energy density of the background field does not dominate the total energy density. We conclude that the spectrum of the density perturbation predicted in the single-field model of the one-bubble open inflation is rather robust. We also consider the spectrum of scalar and tensor perturbations in a model of the Hawking-Turok type, without a false vacuum.

Quasiopenness seems to be generic to multi-field models of single-bubble open inflation. Instead of producing infinite open universes, these models actually produce an ensemble of very large but finite inflating islands. In this paper we study the possible constraints from CMB anisotropies on existing models of open inflation. The effect of supercurvature anisotropies combined with the quasiopenness of the inflating regions make some models incompatible with observations, and severely reduces the parameter space of others. Supernatural open inflation and the uncoupled two-field model seem to be ruled out due to these constraints for values of $\Omega_0\lesssim0.98$. Others, such as the open hybrid inflation model with suitable parameters for the slow roll potential can be made compatible with observations.

We present a non-singular instanton describing the creation of an open universe with a compactified extra dimension. The four dimensional section of this solution is a singular instanton of the type introduced by Hawking and Turok. The ``singularity'' is viewed in five dimensions as a smooth bubble of ``nothing'' which eats up a portion of spacetime as it expands. Flat space with a compact extra dimension is shown to be gravitationally metastable, but sufficiently long lived if the size of the extra dimension is large compared with the Planck length.

The singularity in Hawking and Turok's model (hep-th/9802030) of open inflation has some appealing properties. We suggest that this singularity should be regularized with matter. The singular instanton can then be obtained as the limit of a family of ``no-boundary'' solutions where both the geometry and the scalar field are regular. Using this procedure, the contribution of the singularity to the Euclidean action is just 1/3 of the Gibbons-Hawking boundary term. Unrelated to this question, we also point out that gravitational backreaction improves the behaviour of scalar perturbations near the singularity. As a result, the problem of quantizing scalar perturbations and gravity waves seems to be very well posed.

We show that a large class of two-field models of single-bubble open inflation do not lead to infinite open universes, as it was previously thought, but to an ensemble of very large but finite inflating `islands'. The reason is that the quantum tunneling responsible for the nucleation of the bubble does not occur simultaneously along both field directions and equal-time hypersurfaces in the open universe are not synchronized with equal-density or fixed-field hypersurfaces. The most probable tunneling trajectory corresponds to a zero value of the inflaton field; large values, necessary for the second period of inflation inside the bubble, only arise as localized fluctuations. The interior of each nucleated bubble will contain an infinite number of such inflating regions of comoving size of order $\gamma^{-1}$, where $\gamma$ depends on the parameters of the model. Each one of these islands will be a quasi-open universe. Since the volume of the hyperboloid is infinite, inflating islands with all possible values of the field at their center will be realized inside of a single bubble. We may happen to live in one of those patches of comoving size $d\lesssim \gamma^{-1}$, where the universe appears to be open. In particular, we consider the ``supernatural'' model proposed by Linde and Mezhlumian. There, an approximate U(1) symmetry is broken by a tunneling field in a first order phase transition, and slow-roll inflation inside the nucleated bubble is driven by the pseudo-Goldstone field. We find that the excitations of the pseudo-Goldstone produced by the nucleation and subsequent expansion of the bubble place severe constraints on this model. We also discuss the coupled and uncoupled two-field models.

If the effective cosmological constant is non-zero, our observable universe may enter a stage of exponential expansion. In such case, regions of it may tunnel back to the false vacuum of an inflaton scalar field, and inflation with a high expansion rate may resume in those regions. An ``ideal'' eternal observer would then witness an infinite succession of cycles from false vacuum to true, and back. Within each cycle, the entire history of a hot universe would be replayed. If there were several minima of the inflaton potential, our ideal observer would visit each one of these minima with a frequency which depends on the shape of the potential. We generalize the formalism of stochastic inflation to analyze the global structure of the universe when this `recycling' process is taken into account.

Faddeev and Jackiw's method for constrained systems is used to derive a gauge invariant formulation of cosmological perturbations in the one bubble inflationary universe. For scalar perturbations in a flat universe, reduction of the action to the one with a single physical degree of freedom has been derived in the literature. A straightforward generalization of it to the case of an open universe is possible but it is not adequate for quantizing perturbations in the one bubble universe, because of the lack of Cauchy surfaces inside the bubble. Therefore we perform the reduction of the action outside the lightcone emanating from the center of the bubble or nucleation event, where the natural time constant hypersurfaces are no longer homogeneous and isotropic and as a result the conventional classification of perturbations in terms of scalar and tensor modes is not possible. Nevertheless, after reduction of the action we find three decoupled actions for three independent degrees of freedom, one of which corresponds to the scalar mode and the other two to the tensor modes. Implications for the one bubble open inflationary models are briefly discussed. As an application of our formalism, the spectrum of long wavelength gravity waves is simply obtained in terms of the real part of the reflection amplitude for a one dimensional scattering problem, where the potential barrier is given in terms of the bubble profile.

In the simplest model of open inflation there are two inflaton fields decoupled from each other. One of them, the tunneling field, produces a first stage of inflation which prepares the ground for the nucleation of a highly symmetric bubble. The other, a free field, drives a second period of slow roll inflation inside the bubble. However, the second field also evolves during the first stage of inflation, which to some extent breaks the needed symmetry. We show that this generates large supercurvature anisotropies which, together with the results of Tanaka and Sasaki, rule out this class of simple models (unless, of course, $\Omega_0$ is sufficiently close to one.) The problem does not arise in modified models where the second field does not evolve in the first stage of inflation.

The tunneling approach to the wave function of the universe has been recently criticized by Bousso and Hawking who claim that it predicts a catastrophic instability of de Sitter space with respect to pair production of black holes. We show that this claim is unfounded. First, we argue that different horizon size regions in de Sitter space cannot be treated as independently created, as they contend. And second, the WKB tunneling wave function is not simply the `inverse' of the Hartle-Hawking one, except in very special cases. Applied to the related problem of pair production of massive particles, we argue that the tunneling wave function leads to a small constant production rate, and not to a catastrophy as Bousso and Hawking's argument would suggest.

In the context of the open inflationary universe, we calculate the amplitude of quantum fluctuations which deform the bubble shape. These give rise to scalar field fluctuations in the open Friedman-Robertson-Walker universe which is contained inside the bubble. One can transform to a new gauge in which matter looks perfectly smooth, and then the perturbations behave as tensor modes (gravitational waves of very long wavelength). For $(1-\Omega)<<1$, where $\Omega$ is the density parameter, the microwave temperature anisotropies produced by these modes are of order $\delta T/T\sim H(R_0\mu l)^{-1/2} (1-\Omega)^{l/2}$. Here, $H$ is the expansion rate during inflation, $R_0$ is the intrinsic radius of the bubble at the time of nucleation, $\mu$ is the bubble wall tension and $l$ labels the different multipoles ($l>1$). The gravitational backreaction of the bubble has been ignored. In this approximation, $G\mu R_0<<1$, and the new effect can be much larger than the one due to ordinary gravitational waves generated during inflation (unless, of course, $\Omega$ gets too close to one, in which case the new effect disappears).