Quantum cosmology with third quantisation

We review the canonical quantisation of the geometry of the spacetime in the cases of a simply and a non-simply connected manifold. In the former, we analyse the information contained in the solutions of the Wheeler-DeWitt equation and interpret them in terms of the customary boundary conditions that are typically imposed on the semiclassical wave functions. In particular, we review three different paradigms for the quantum creation of a homogeneous and isotropic universe. For the quantisation of a non-simply connected manifold the best framework is the so-called third quantisation formalism, in which the wave function of the universe is seen as a field that propagates in the space of Riemannian $3$-geometries, which turns out to be isomorphic to a (part of a) $1+5$ Minkowski spacetime. Thus, the quantisation of the wave function follows the customary formalism of a quantum field theory. A general review of the formalism is given and it is analysed the creation of the universes, their initial expansion and the appearance of matter after inflation. These features are presented in more detail in the case of a homogeneous and isotropic universe. The main conclusion in both cases is that the most natural way in which the universes should be created is in entangled universe-antiuniverse pairs.

Quantum cosmology is the application of the quantum theory to the universe as a whole. However, it was clear from the beginning that the customary formalism of the Copenhagen interpretation cannot be applied to the quantisation of the universe because the Schrödinger equation and the measurement process in which the Copenhagen formalism is based on cannot be fundamental elements of a quantum theory of the spacetime. Let us notice that if the quantum theory must represent the quantum state of the spacetime, then, as Wheeler showed [1], its quantum state must be at the Planck length a superposition of spacetime geometries that is impossible to visualise or represent and where one cannot even use the word "observation" at all (cf. [1]). However, as we approach the macroscopic scale the quantum state of the universe must represent the approximately stable spacetime where we live and perform measurements of particles and other matter fields. It means that the description of the universe that we observe must be an emergent feature of the quantum representation of the universe.
In this paper, we review the canonical formulation of quantum cosmology. We start from the foliation of the spacetime into space and time that allows us to express the Einstein-Hilbert action as a functional action of the components of the 3-metric of the spatial sections of the spacetime. The evolution of the universe turns out to be then a trajectory in the space of 3-Riemannian metrics, M , and its quantum state is represented by a wave function that is the solution of the Wheeler-DeWitt equation that, in principle, contains all the information about the spacetime and the matter fields that propagate therein. However, as we have already said, the full quantum state is a superposition of solutions that correspond to different paths, i.e. different evolutions, in the space M . It is only in the semiclassical regime where a particular kind of solutions emerge by a decoherence process 1 . This kind of solutions are the semiclassical solutions that represent a fixed classical spacetime background with matter fields propagating therein, and in which the Schrödinger equation appears as an approximated equation at order 1 . In particular, we analyse the semiclassical wave function of a homogeneous and isotropic universe with small inhomogeneities that can be treated as perturbations. In that case, explicit solutions can be given for which it is easier to analyse the different boundary conditions that can be imposed on the state of the universe. They give rise to different scenarios for the creation of the universe that are analysed in detail.
On the other hand, the space of 3-dimensional space-like metrics, M , defined at any point in the space, turns out to be isomorphic to a 1 + 5 dimensional Minkowski spacetime 2 . This analogy between the space M and the spacetime allows us to consider the wave function of the universe as a field that propagates in M , and the Wheeler-DeWitt equation as the field equation. In that case, a procedure of quantisation called third quantisation can formally be performed in a similar way as it is done in a quantum field theory. For instance, we can define quantum operators representing the creation and annihilation of particular modes of the spacetime, i.e. different universes, and the corresponding Fock space will then allow us to represent the quantum state of a multiply connected spacetime manifold. It turns out to be then the appropriate framework to describe the quantum state of the multiverse. Moreover, as it happens in a quantum field theory, the isotropy of the background space implies that the creation of universes must be in pairs with opposite values of the components of the momentum conjugated to the configuration variables. We shall analyse the charge and parity relation between the matter fields that propagate in one of these pairs to see that the matter content of one of the universes must be CP inversely related with the content of the other universe. They thus form a universe-antiuniverse pair. We analyse this pair creation in detail in the case of a homogeneous and isotropic universe where the period of reheating after inflation is investigated. The decay of the inflaton field into the particles of the Standard Model is produced in a CP conjugated way in the two universes so any excess of matter over antimatter in one of the universes of the entangled pair is balanced with the excess of the antimatter over matter in the partner universe, having always these two concepts (matter and antimatter) a relative meaning, i.e. an internal observer in any of the universe always interprets the content of his/her universe as matter.
This proposal is however still far from being directly testable. The effects of quantum cosmology and, in particular, the effects of the existence of an entangled universe are mainly restricted to the very early stage of the universe. Perhaps with the future advances in the detection of gravitational waves or of a cosmic neutrino background we will be able to test the pre-inflationary stage of the universe where the effects of quantum gravity may be significant. Moreover, the third quantisation formalism can also be a proposal for the quantisation of the spacetime and thus a better understanding of the formalism and its application to other gravitational scenarios can provide us with a new line of research for the search of a quantum theory of gravity.

II. QUANTISATION OF A SIMPLY CONNECTED SPACETIME MANIFOLD
A. Quantisation of the spacetime geometry Following the customary approach 3 , the spacetime can be foliated in space and time by assuming a global time function t such that each surface t = constant is a spacelike Cauchy hypersurface, Σ t . The proper distance between the point x 0 in Σ t and the point x 0 + dx of Σ t+dt is given by [7,8] (see, Fig. 1) where N a N a = h ab N a N b , and h ab is the three-dimensional metric induced on each hypersurface Σ t , with unit normal n µ , satisfying, n µ n µ = −1. The functions N and N a are called the lapse and the shift functions, respectively. They are the normal and tangential components of the vector field t µ , which is the vector field that transport the point x 0 from Σ t to Σ t+dt . With the split of the spacetime in space and time the spacetime can be seen as a spacelike hypersurface evolving in time. The geometry of the hypersurface Σ t at a given time t 0 is determined by the metric tensor h ab (t 0 ), so eventually the evolution of the universe is encoded in the evolution of the metric h ab (t). It is then remarkable that in the end, as Wheeler says [9], Eintein's geometrodynamics deals with the dynamics of 3-geometry, not 4-geometry! (emphasis his). From this 3 + 1 viewpoint of the spacetime, we can cast the Einstein-Hilbert action 4 [7] S EH and the action of the matter fields 5 into the standard Lagrangian form, where q i andq i will be here the spatial metric, h ab (x), the matter field(s), encoded in the variable ϕ(x), and their corresponding velocities. In (2), 4 R is the Ricci scalar associated to the 4-dimensional metric g µν , Λ is the cosmological constant, K = h ab K ab , is the trace of the extrinsic curvature, which can be written as, where, D a is the covariant derivative on the spatial section Σ t , and ∂M is the boundary of the manifold M. After some manipulation (see Ref. [7]), one finds that the Einstein-Hilbert action (2) can be written as where h is the determinant of the spatial metric h ab , and FIG. 1. Splitting the spacetime into space and time.
is the so-called DeWitt's metric [10]. The structure of the action (6) is very interesting. First, it is of the standard form (4), with L EH the Lagrangian associated to the Einstein-Hilbert action (6) and L EH the Lagrangian density. Second, from (5) one can see that the extrinsic curvature K ab contains the time derivative of the metric tensor h ab , and 3 R only depends on h ab . Thus, the action (6) presents the customary structure of a kinetic term that is quadratic in the velocities plus a potential term. Furthermore, the supermetric (7) defines a metric structure on the space of spacelike metrics, called the superspace 6 . Thus, the action (6) looks like the action of a particle that moves in a curved space, the coordinates of the "particle" are the time dependent values of the components of the metric tensor, h ab (t), and the curved space where this particle moves in the space of symmetric Riemannian 3-metrics. That is, the evolution of the universe can be seen as the trajectory in the superspace 7 (see, Fig. 2) . The momenta conjugated to the metric tensor components, h ab , are given by [7] p ab ≡ ∂L EH ∂ḣ ab = G abcd K cd = √ h K ab − Kh ab .
Thus, in terms of the momenta, the total action (the gravitational action plus the action of the matter field) can be written as where the lapse and the shift functions act as Lagrange multipliers, with [7] where [7], J a ≡ h µ a T µν n ν , and G abcd is the inverse of the DeWitt metric (7), with [10], FIG. 2. Left: the evolution of the universe can be seen as a path in the abstract space spanned by the component of the spatial metric tensor, h ab . Right: an example of evolution is depicted for the case that the spatial metric is diagonal, i.e. given by, h ab (t) = diag(h11(t), h22(t), h33(t)).
Therefore, variation of the action (10) with respect to the lapse and the shift functions yields the classical Hamiltonian and momentum constraints, respectively, i.e., Let us now focus for a moment on the gravitational part of these constraints. The Hamiltonian constraint, H = 0, can be seen as the analogue of the momentum constraint of a particle that propagates in the spacetime, g µν p µ p ν + m 2 = 0, (16) in which the mass is substituted by a non constant potential (this analogy will be further exploited in Sec. III).
On the other hand, the function H a generates infinitesimal diffeomorphisms (change of coordinates) in the spatial hypersurfaces, Σ. Thus, the second constraint in (15), H a = 0, means that the Einstein-Hilbert action (6) is invariant under such diffeomorphisms, which turns out to be like a gauge freedom [7]. Thus, the real configuration space is the quotient space of all Riemannian 3-metrics, M ≡ Riem(Σ), in which all three metrics related by diffeomorphisms correspond to the same class, i.e.
which is called the superspace [7][8][9][10]. Following Dirac [11], the canonical procedure of quantisation consists in assuming the quantum version of the classical constraints (15) by promoting the classical variables and their conjugate momenta into quantum operators and apply them to a wave function, φ, that is defined in the configuration space, The quantum version of the momentum constraint in (15),Ĥ a φ = 0, assures that the wave function φ is invariant under spatial diffeomorphisms in the 3-dimensional slices Σ t [7]. For our purposes, it is much more interesting the constraint,Ĥφ = 0, which under canonical quantisation becomes the so-called Wheeler-DeWitt equation [7,8,10], where, ∇· and ∇, are the divergence and the gradient, respectively, defined in the space of 3-dimensional Riemannian metrics, M , and for a scalar fieldT 00 readŝ T 00 = − 2 2h The Wheeler-DeWitt equation (19) is the keystone of the canonical formalism of quantum cosmology. The solution, φ(h ab , ϕ), is usually called the wave function of the universe [12] because it represents the quantum state of the spacetime and the matter fields that propagate therein. This is usually applied to the case of a single universe. However, as we shall see in Sect. III, the whole spacetime manifold can present a more complicated structure and represent something more than what is typically called a universe. In that sense, the name can be misleading. However, for historical reasons we shall retain sometimes the name wave function of the universe even in the cases where it may represent the state of many different universes. The Wheeler-DeWitt equation (19) can be seen as a Schrödinger like equation with no time variable, which is a consequence of the invariance of the quantum state of the universe with respect to the time variable. In that sense, there is no preferred time in the quantum description of the universe (as well as there are no preferred spatial coordinates because the constraint, H a = 0). It is then sometimes stated that there is no time evolution of the quantum state of the universe. However, this is not true or at least it is not accurate. As we have already pointed out, the evolution of the universe can be seen as the trajectory in the superspace. The spacetime coordinates are the parameters that parametrise the trajectory, which is therefore invariant under reparametrisations but that does not mean that there is no evolution. It is similar to the description of the path followed by a particle in the spacetime, which is independent of the parametrisation of the path but that does not mean that the particle does not move.
There is also a path integral approach to quantum gravity and quantum cosmology [7,8,12]. It is a generalisation of Feynman's idea that the amplitude for a particle to go from one to another point point is given by a functional integral that weights all the paths that start from the point x 0 at time t 0 and end in the point x 1 at t 1 (see Fig. 3). Following a parallel reasoning, Hartle and Hawking propose [12] that the amplitude for the universe to change from the hypersurface Σ, in which the spatial geometry and the field configuration are given by h ab and ϕ, respectively, to the hypersurface Σ , where they are given by the values h ab and ϕ , is given by [12] h ab , ϕ |h ab , ϕ = δg δϕ e iS [g,ϕ] , where the integral must be performed over all 4-geometries and field configurations that match the given values on the two spacelike hypersurfaces [12], Σ and Σ . Following that approach, the wave function of the universe is then given by, where C denotes the class of spacetimes and matter configurations that fulfil the boundary requirements on the hypersurfaces Σ, and the sum is performed over all kind of topologies. In order to make well defined the path integral in (22) one has to rotate to Euclidean time. However, that does not remove all the technical problems, which we are not going to deal with here. In practice, as it happens with the Wheeler-DeWitt equation, the path integral can only be performed for spacetimes and matter field configurations with a high degree of symmetry. Moreover, both formulations become equivalent because the requirement of invariance of the wave function φ(h ab , ϕ) under reparametrizations of the time variable implies the constraint [12], δS δN = 0, whose quantised version is the Wheeler-DeWitt equation. Anyway, the path integral formulation has two interesting points. First, the analogy with the Feynman's path integral formulation of the trajectory of a particle in spacetime makes it very intuitive. As we have seen, the classical evolution of the universe can be seen as a path in the superspace. Applying Feynman's idea to gravity means that the quantum state of the universe is given by a quantum superposition of all the paths that go from one to another configuration of the spatial hypersurfaces. As it happens in the spacetime, the classical path (i.e. the classical evolution) emerges in some specific limit because the constructive interference between the paths of the quantum superposition, and in the same limit the non-classical paths suffer from a destructive interference or decoherence [3,6,13,14]. The second interesting point of the path integral approach is the following: let us assume that we already know how to construct the quantum amplitude for the universe to go from one to another given configuration of the spacelike section. Then, what is it the amplitude for the birth of the universe? What are the boundary conditions that one must impose on the state of the universe to obtain the appropriate probability amplitude for the universe to be created?

B. Boundary conditions
In classical mechanics, as well as in quantum mechanics, we usually work with some given conditions that we know or assume for certain at some initial time. Then, knowing the law of evolution, the initial conditions determine the state of the system at any later time. In the universe the thing is a bit different. What we only know for certain is the state of the observable universe, say from the inflationary period 8 to the current stage of accelerated expansion, and we have to make some guess about the initial conditions that give rise a universe like that. But this is all classical cosmology. The question in quantum cosmology is, what are the conditions at the quantum level that give rise a specific initial boundary, Σ 0 , that is propitious to inflate? The probability for the creation of such universe would be given by the modulus squared of the wave function of the universe, when this is evaluated at the initial hypersurface 9 . From that point of view, one can say that the initial state of the (classical) universe is the final state of the amplitude for the universe to be created ... from what? And that is actually the issue behind the question of the boundary condition of the universe.
Using the path integral approach, Hartle and Hawking propose [12] that the class C over which the integral (22) has to be performed is the class of compact geometries (in principle of all topologies) that have Σ 0 as their only boundary [7,12] (see Fig. 4), and matter fields that are regular on those geometries. It means that the boundary of the universe is that the universe has no boundary, or equivalently, that the boundary Σ 0 is created from nothing. We shall see later on that for the case of an inflationary spacetime the quantum state that results from the no-boundary condition is [7,12,15] where V (ϕ) is the potential of the inflaton field and a is the scale factor, which goes from the initial value, a 0 = V (ϕ) −1/2 , to infinity. It is important to notice that the wave function (23) can be written as where, We shall see in the next section that in terms of the same time variable one of the two terms in (24), say the branch with e −iS , describes an expanding universe and the branch with e iS describes a contracting universe. It means that the result of imposing the no-boundary condition on the quantum state of the universe is that it is given by the linear combination of two states: one representing an expanding universe and one representing a contracting universe. Typically, one considers the expanding branch of the universe as representing our universe and disregard the contracting branch for being unphysical. However, we shall see that there is another interpretation. In terms of the physical time variable of each universe the two universes can both be seen as expanding universes but with their 9 Let us notice however that the creation of this initial boundary hypersurface Σ 0 is not a process occurring in time but it corresponds to the creation of the spacetime itself [7], actually.
FIG. 4. Left: the class C over which the path integral has to be performed is the class of compact geometries that have Σ0 as their only boundary [12]; Right: the tunnelling proposal states that the only modes that survive the quantum barrier are the 'outgoing' modes that represent expanding universes.
matter fields being CP conjugated. They can be interpreted then as an expanding universe-antiuniverse pair (see, Sec. II E 3 and Sec. III). From that point of view, the no-boundary proposal would yield the creation of universes in entangled universe-antiuniverse pairs. Vilenkin's tunneling proposal is quite different. Perhaps more based on practical grounds, he proposes that the only mode that survives the tunnelling from the Euclidean region (the region located at, a < V (ϕ) −1/2 ) is the one that represents an expanding universe. Imposing the tunnelling boundary condition, the resulting wave function of the universe is The main difference with respect to the Hartle-Hawking wave function (23) is the negative sign in the exponent of the exponential pre-factor, which may have important consequences. Let us notice that the probability for the universe to be created from nothing is, P ∝ |φ| 2 , so in the case of the no-boundary condition we have while in the case of the Vilenkin's tunnelling condition, we have One immediate consequence is that the Vilenkin's condition seems to favour the creation of a universe with a big value of the potential, which is a necessary condition for the initial hypersurface Σ 0 to inflate. That would in principle reject the Hartle-Hawking proposal because, on the contrary, the no-boundary proposal seems to favour the creation of a universe with a small (or zero) value of the potential (P HH → ∞, as V → ∞). However, this result changes when high order corrections are taken into account so the result is not conclusive [7]. From a purely theoretical point of view, it seems that the Hartle-Hawking proposal is more fundamental in the sense that these authors put the focus on the natural condition that one should impose on the Euclidean region of the spacetime. As a consequence, they obtain that the universe is represented by two branches, one corresponding to an expanding universe and other describing a contracting branches. These two branches can be considered independently once they suffer a process of decoherence, so in practice they may represent two different universes. Vilenkin's proposal seems to be more practical (although it is also based on a parallelism with some processes of quantum mechanics).

C. Semiclassical quantum gravity
The first thing that the quantum state of the spacetime must provide us with is a consistent explanation of how the classical background can emerge from the full quantum state of the spacetime. Fortunately, not only this can be done in a beautiful manner but it is in fact one of the greatest achievements of quantum cosmology.
Again, the path integral formulation of the spacetime supplies a clear picture of how it can be obtained. At some appropriate limit, which is generally a large length or mass scale compared with their corresponding Planck values, the contribution of most of the paths in the integral vanishes because their destructive interference. The only paths that survive the interference are those that are in phase, i.e. those for which δS ≈ 0. These are actually the trajectories of the superspace given by the classical constraints. Therefore, much in a similar way as the classical trajectory of a particle emerges from the constructive interference in the path integral approach of the quantum mechanics of a particle, the classical background spacetime emerges as the constructive interference among the paths in the superspace. In the quantum mechanics of a particle we can then compute the quantum corrections to the trajectory of the particle in terms of quantum uncertainties. We shall see in this section that the quantum corrections to the classical background of the universe are caused by the matter fields.
Following the customary approach [4,7,8,15], let us consider the following semiclassical wave function [4], where ∆ and ψ are slowly varying functions of the metric tensor h ab , and S 0 (h ab ) is the classical action for gravity alone, given by (2). Essentially, the wave function (29) contains two parts: one that only depends on the geometric variables of the spacetime and another part that contains all the matter degrees of freedom in the wave function ψ. The basic idea is that we expect that the quantum fluctuations of the spacetime will weaken more rapidly than the quantum fluctuations in the state of the matter fields. In that case, we shall find a regime where the spacetime behaves nearly classically with quantum matter fields propagating therein. That is exactly what we need to describe the universe we observe. Also notice the presence of the Planck constant in (29). It means that the classical behaviour of the spacetime will be present whenever the gravitational action is large with respect to the Planck constant. Now, insert the semiclassical wave function (29) into the Wheeler-DeWitt equation (19) and solve it order by order in in the geometrical degrees of freedom. At order 0 , one finds which is the gravitational part of the Hamiltonian constraint (11) if we make the identification, At first order in , neglecting second derivatives of the slowly varying terms with respect to the 3-metric, it is obtained two equations. The first equation [4] G abcd δ δh ab is actually the condition for the function ∆(h ab ) to be a slowly varying function. In other words, whenever (32) is satisfied (or to the extent it is satisfied) the wave function (29) can be a good candidate to describe the observable universe. Eq. (32) is also the equation of the conservation of the probability current, ∆ 2 p cd . The other equation that is obtained at order is where the ∓ signs correspond to the ± signs of the exponent of (29). It suggests the identification of a time variable t, given by In that case, (33) becomes the Schrödinger like equation of the matter fields Therefore, in the semiclassical regime we obtain at order 0 the classical behaviour of the spacetime and, at first order in , the quantum evolution of the matter fields. Thus, the semiclassical wave function (29) contains all the physical information of the observable universe. One could say that recovering the classical equations for the spacetime degrees of freedom and the Schrödinger equation for the matter fields does not add anything to what we already knew before the quantisation of the universe. That is true, these two features are nothing more than a test of consistency for quantum cosmology. After all, the recovering of the classical spacetime must not be surprising. We started the process of quantisation from the classical action of the spacetime and the matter fields, and the quantisation procedure consists basically in promoting the classical variables to non commuting quantum operators, [ĥ ab ,p ab ] ∝ . Then, it should not be surprising that in the limit → 0 we recover the classical behaviour (this is the essence of the correspondence principle). And something similar for the quantum behaviour of the matter fields. Even though, the canonical quantisation of the whole universe that we have seen in this section is interesting for several reasons. First, it suggests that the quantisation processes that leads to the wave function of the universe is consistent and, in that case, one can assume that the wave function of the universe, φ(h ab , ϕ), would contain in principle all the physical (classical and quantum) information of all the degrees of freedom of the universe. Any physical process should be describable within the formalism of quantum cosmology. Of course, this reductionist point of view is not practical at all but from the conceptual point of view it results appealing. A more interesting feature is that it allows us to analyse higher order corrections to the semiclassical universe, and this should give novel features that cannot be foreseen in the classical scenario. It might help us to go beyond the quantum description of matter fields in a classical spacetime background. In particular, it can help us to find some exclusive features of the quantum regime of the spacetime, i.e. small deviations from the known behaviour caused by the high order corrections of quantum gravity [16][17][18][19].
Another interesting feature of quantum cosmology is the appearance of time. From the point of view of the superspace, time is just the parameter that parametrises the curve that describes a particular trajectory. In the picture given by the path integral, the quantum state of the universe is given by the set of all paths that join together the initial and final states. If the quantum wave-packet is spread no definite time variable can be chosen mainly because there is no definite curve that describes the evolution of the universe. It is only when the wave-packet is peaked around a particular solution (ideally becoming a delta function) when we have a definite curve, the classical evolution of the universe, that can be therefore parametrised in terms of a parameter that we can call time 10 . It therefore appears as the result of a decoherence process between the different histories of the universe [3,4,6,20].
Furthermore, quantum cosmology relates the two concepts of time of contemporary physics: the one of the theory of relativity and the one of quantum mechanics. This is a very subtle point. Both the theory of relativity and the quantum theory work with a time variable, say t r and t q , respectively. We usually assume that both time variables are the same, t r = t q = t, but this is an assumption that is not guaranteed from the beginning. Of course, they (must) coincide in the Newtonian limit of both theories, but in general, they only coincide if the time variable of theory of relativity would be measured with an actual clock, which is made of matter fields. However, the theory of relativity deals with 'ideal clocks' and the consideration of an actual clock may entail some problems 11 .

D. Minisuperspace model
Despite its conceptual importance, it is not hard to see that the Wheeler-DeWitt equation found in the previous chapter is very difficult if not impossible to solve for a general configuration of the spacetime and the matter field(s). In order to make computations one generally has to assume some symmetries in the underlying spacetime of the universe. This process, called symmetry reduction [7], reduces the number of variables of the superspace and the so reduced superspace is called minisuperspace 12 .
Furthermore, a minisuperspace model is not necessarily an unrealistic model. On the contrary, the observational data indicates that most of the history of the universe the spacetime presents a high degree of symmetry. Even more, if the initial hypersurface Σ 0 from which the universe starts evolving is large enough compared with the Planck length, a minisuperspace model could describe the whole history of the universe 13 . In other cases, it can be taken as a toy model from which we can obtain relevant information about some (classical and quantum) aspects of the universe like, for instance, the creation of the initial hypersurface Σ 0 or the quantum-to-classical transition and the appearance of time, among others.
They can also allow us to study the effect of small deviations from the symmetric picture. For instance, we shall consider later on small departures from the homogeneous spacetime and the matter field in the form of gravitons and matter particles, respectively. These are local perturbations of the otherwise homogeneous and isotropic background.
The picture becomes then quite realistic and allows us to analysed potentially observable effects like the kind of correlations between the modes of the matter fields in different regions of the spacetime or the quantum gravitational corrections to the Schrödinger equation of the matter fields. In all those cases, the study of the minisuperspace model turns out to be justified.
Therefore, let us consider the minisuperspace that is obtained from the foliation of a 4-dimensional spacetime with closed homogeneous and isotropic spatial sections. The geometry of the spacetime is then characterised by a Freedman-Robertson-Walker (FRW) metric where dΩ 2 3 is the line element on the unit three sphere. The foliation of the spacetime into space and time is in that case parametrised by just two functions, the scale factor a(t) and the lapse function N (t). The geometry of the spatial sections, h ab , is then fully characterised by the scale factor a(t) that parametrises the variation in the distance between two fixed points of the space along the evolution of the universe. It parametrises therefore the expansion or the contraction of the universe. On the other hand, the lapse function N (t) determines the time parametrisation of the foliation. Different values of N (t) entail different time variables (i.e. different time parametrisations), with some special cases. For instance, if N = 1 the time variable t is called cosmic time and if N = a(t), t is customary renamed with the Greek letter η and is called conformal time because in terms of η the metric becomes conformal to the metric of a closed static spacetime. Anyway, we know that the evolution of the universe is invariant under the choice of time reparametrisation so, at the end of the process, we can take the preferable time variable 14 .
The line element of the spacetime is then fully determined by these two functions, a(t) and N (t). The the total action, i.e. the Einstein-Hilbert action (2) plus the action of the scalar field (3), can be written as [7] where an integration over the spatial variables has been performed and absorbed with a definition of units in which, 2G/3π = 1, and the rescalings, ϕ → ϕ/ √ 2π and V → V /2π 2 . The total action has been simplified considerably. The only dynamical degrees of freedom are the scale factor, a(t), and the scalar field, ϕ(t), which according to the homogeneity condition must only depend on the time variable. The superspace has then been reduced to a two dimensional space. Now, we can proceed as described in the preceding sections. The momenta conjugated to the configuration variables are, and the Hamiltonian then reads The momentum constraint is automatically satisfied by the symmetries of the spacetime and the Hamiltonian constraint, δH δN = 0, becomes then, H = 0. Promoting the momenta into quantum operators and applying the quantum version of the Hamiltonian constraint, H = 0, to the wave function of the universe, φ(a, ϕ), which depends now on the two variables of the minisuperspace, a and ϕ, we obtain the Wheeler-DeWitt equation,Ĥφ(a, ϕ) = 0.

Inflationary universe
Let us first analyse the initial stage of the universe where the scalar field is assumed to be approximately constant on the small time scale in which the universe rapidly undergoes an inflationary expansion [24],φ ≈ 0 and ϕ 0 1. In that case, the kinetic term of the scalar field in (37) can be neglected and the potential turns out to be approximately constant, V (ϕ 0 ). From (37) it can be seen that a constant value of the potential is equivalent to a cosmological constant term, so the universe effectively behaves like a DeSitter spacetime. Then, we can write (37) as where, H 2 0 = 2V (ϕ 0 ), in the case of the inflationary universe or, H 2 0 = Λ/3, in the case of a 'pure' DeSitter spacetime, or the sum of both in a general case. The corresponding Hamiltonian constraint turns out to be which is nothing more than the Friedmann equation expressed in terms of the momentum conjugated to the scale factor, p a . In terms of the time derivative of the scale factor, using p a = aȧ (in cosmic time 15 , with N = 1 in (38)), the Hamiltonian constraint (41) can be written asȧ whose solutions can easily be obtained, If t ∈ (−∞, ∞), the scale factor (43) describes a universe that starts shrinking from an infinite volume, bounces at the minimum value, a 0 = H −1 0 , and ends up in an eternal expansion (see, Fig. 5 Left). However, it does not seem quite plausible that the universe is created with an infinite volume, so the most reasonable possibility consist in restricting ourselves to the domain, t ∈ (0, ∞), which describes a 'bubble' of spacetime that is created with radius a 0 = H −1 0 at t = 0 (the origin of time), and starts expanding exponentially in an inflationary like expansion. If the length scale of the initial 'bubble' is some orders of magnitude greater than the Planck scale, i.e. H −1 0 >> l P , then, the quantum fluctuations of the spacetime would be small and the homogeneous and isotropic picture described here could reasonably represent the creation of a universe like ours 16 .
Let us now analyse the solutions of the corresponding Wheeler-DeWitt equation. Making the substitution, p a → −i ∂ a , and leaving aside the ambiguity of the factor ordering, the Hamiltonian constraint (41) can conveniently be written as where, 15 Recall that after doing the variation with respect to N we can fix any particular value. 16 However, one would generally expect that the creation of the universe comes from a quantum fluctuation of the spacetime, and therefore be of order of the Planck scale, H −1 ∼ l P . In that case, new elements should be incorporated although the picture described here and in the next section would still be instructive.
Written in this way, the Wheeler-DeWitt equation (44) resembles the equation of a harmonic oscillator with time dependent frequency (such parallelism will be further exploited in Sec. III). Equation (44) is not exactly solvable. However, far from the turning point, a 0 = H −1 0 , we can approximate the solutions by the WKB wave functions where N ± is a normalisation constant and S(a) is given by The two wave functions φ ± in (46) represent incoming and outgoing wave functions. Let us see it by inserting them into the Wheeler-DeWitt equation (44). In that case, at order 0 it is obtained the Hamilton-Jacobi equation which corresponds to the Hamiltonian constraint (41) if one makes the identification In fact, let us note that at leading order in it is satisfied Therefore, it is obtainedȧ which is the Friedmann equation (42), with two signs: the − sign that corresponds to φ + and the + sign that corresponds to φ − . Thus, φ − describes an outgoing wave, i.e. a wave that travels towards greater values of the scale factor -it thus represents an expanding universe -and φ + an incoming wave, i.e. a wave that travels towards smaller values of the scale factor, which corresponds therefore to a contracting universe. Let us notice however that the interpretation of φ ± in terms of incoming and outgoing wave must be taken carefully [25]. The Friedmann equation (41) is invariant under the time reversal change, t → −t, and in terms of −t (ȧ → −ȧ), φ + would represent the expanding universe (i.e. an outgoing wave) and φ − the contracting universe (i.e. an incoming wave). Clearly, it depends on the (time) parametrisation. The important thing is that the general quantum state of the universe (far from the turning point) can be written as which represents a quantum superposition of incoming and outgoing waves, irrespective of the particular wave function that represents each one. Another important feature that is worth noticing is that the minimum value a 0 constitutes a classical barrier below which the universe cannot go through. Let us notice that for a value a < a 0 there is no real solution of the Friedmann equation (51) so an incoming wave function representing a contracting universe is classically reflected (bounced) into an outgoing wave (see, Fig. 5 Right).

Small perturbations and backreaction
Let us now consider a more realistic scenario by introducing two important changes. First, we shall not neglect the kinetic term of the scalar field in the action (3), and we shall consider a general form (i.e. not necessarily a constant) for the potential V (ϕ). Contrary to what may be thought, that will not introduce qualitative changes. In return, it will allow us to represent other stages of the evolution of our universe as well as many other types of universes. The second change that we are going to make is the introduction of small perturbations around the homogeneous and isotropic background spacetime. This will help us to analyse several phenomena like the behaviour of the matter fields in the semiclassical universe or the appearance of a physical time variable.
Regarding the last question, we can assume that except for the very beginning the inhomogeneities of the universe are relatively small 17 . Therefore, to a good order of approximation the universe can be represented by a homogeneous and isotropic background with relatively small inhomogeneities propagating in the background. In that case, it seems reasonable to expand the variables of the spacetime and the matter fields around the homogeneous and isotropic values and study the inhomogeneities as perturbations of the homogeneous and isotropic background. We can still consider the 3 + 1 splitting of the spacetime given in (1). However, the idea is now to expand the configuration variables (h ab , N, N a and ϕ) around their homogeneous and isotropic values and retain just the first order terms. Then, let us consider the following expansions [2,26] where Ω ab is the metric on the unit three-sphere, ϕ(t) is the homogeneous mode of the scalar field, Q n (x) are the scalar harmonics on the three-sphere and G n ij (x) the transverse traceless tensor harmonics [26], with, n ≡ (n, l, m). More harmonics can be present in (53). However, we shall only focus on the tensor modes, d n , as the representative of the perturbation of the spacetime. Eventually, these modes will represent gravitons propagating in the background spacetime and, analogously, the perturbation modes f n will represent the particles of the scalar field(s). The lapse and shift functions must also be expanded in terms of the spherical harmonics. Then, all these perturbed functions are inserted in the action (10). The configuration variables are now the scale factor a(t), the homogeneous mode of the scalar field ϕ(t), and the infinite number of modes d n (t) and f n (t), denoted generically by x n (t). The minisuperspace is then the infinite dimensional space spanned by the variables (a, ϕ, x n ), and the time evolution of the universe is represented by a parametrised trajectory in that space, (a(t), ϕ(t), x n (t)). We should have included the modes of the perturbed lapse and shift functions. However, as it happens with their homogeneous counterparts they are not dynamical variables but Lagrange multipliers that generate a set of constraints that can be used to simplify the equations [26]. After a cumbersome computation, one arrives at the Hamiltonian constraint where the Hamiltonian H 0 contains only degrees of freedom of the homogeneous and isotropic background, and H m is the Hamiltonian density that contains also the degrees of freedom of the perturbation modes. Let us first consider the Hamiltonian of the background spacetime, which is given by the Hamiltonian density of (39), After canonical quantisation, the Wheeler-DeWitt equation reads where, H 2 (ϕ) ≡ 2V (ϕ), is not necessarily constant now. Let us consider the semiclassical solutions, where ∆(a, ϕ) is a slow varying field of the scale factor. Inserting the wave function (58) into the Wheeler-DeWitt equation (57), as we did in the previous section, and disregarding second order derivatives with respect to the background variables, one obtains at order 0 the Hamilton-Jacobi equation [2] − ∂S ∂a Following the semiclassical development of the previous sections, let us define a WKB-time parameter t w as [2] in terms of which,ȧ and the Hamilton-Jacobi equation (59) turns out to bė which is the Friedmann equation of the background spacetime. The WKB wave functions φ ± 0 describe universes with a background spacetime that evolves according to the Friedmann equation (62). The wave function φ 0 may thus represent the quantum state of a large variety of classical models of the universe. The two signs in (60) are irrelevant at the classical level because the Friedman equation (62) is invariant under the reversal change of the time variable. However, they will play an important role at the semiclassical level.
Let us now consider the complete Hamiltonian constraint (55). The corresponding Wheeler-DeWitt equation is whereĤ 0 is the Hamiltonian of the background that we have already seen (56)(57), andĤ m is the Hamiltonian of the perturbation modes, which for the moment we do not need to specify. The wave function of the universe can be separated in two factors, the wave function φ 0 (a, ϕ) of the background (58) and another wave function that contains the matter degrees of freedom 18 . Therefore, the semiclassical wave function can be written now as If we insert the wave function (64) into the Wheeler-DeWitt equation (63) and we solve it order by order in , it is obtained at order 0 the Hamilton-Jacobi equation (59). Therefore, the wave function (64) still describes the background spacetime that evolves according to (62). On the other hand, at order 1 in H 0 , one obtains Here comes a subtle but crucial point of the semiclassical regime. In terms of the initial proper time, t, the two branches of the universe represent an expanding and a contracting universe because, from (38), so φ − 0 describes a universe whose spacetime background is expanding and φ + 0 a universe whose background spacetime is contracting. In that case, in order for the WKB-time (60) to represent the proper time variable t in the two branches, we have to choose for the branch φ − the WKB-time variable t − defined by in which case, that represents an expanding universe. For the WKB-time variable in the φ + branch, we must choose in terms of which, FIG. 6. An electron-positron pair can equivalently be seen as either an electron propagating backwards in time, bouncing, and then propagating forward in time; or as a particle-antiparticle pair propagating forward in time.
describes a contracting universe. It is worth noticing that this assignation is somehow arbitrary [25,27]. If we would have started with a time variable, t → −t, then, φ + 0 would represent the contracting universe and φ − 0 the expanding one and the assignations of the WKB-time variables would have been the other way around. However, in terms of the definitions (69) and (67) the equation (65) reads, where , ϕ], evaluated in the background solutions of (70) and (68). Therefore, we have ended up with two universes, one contracting and another expanding, both filled with matter. There is however a different interpretation. One may assume that the physical time variable, i.e. the time variable measured by actual clocks that are made of matter, is the time variable that appears in the Schrödinger equation. In that case, it is worth noticing that the physical time variable of the two universes is reversely related, t + = −t − . Let us assume that we fix the time variable by fixing the time that a particular observer measures and consider thus t − as the physical time. Then, in terms of the time variable t − the evolution of the scale factor is given by (70) so the two wave functions, φ + and φ − , represent both with a universe with an expanding background spacetime, i.e from the point of view of this hypothetical observer both universes are expanding. The Schrödinger equation (65) in the observer's universe becomes, However, the Schrödinger equation for the fields of the partner universe is for the wave function ψ + . The 'wrong sign' in (73) is not problematic. It is only indicating that (73) is the Schrödinger equation of the complex conjugated wave function ψ * + with a CP -transformed Hamiltonian [25]. Let us notice that (73) can be written as It is therefore the Schrödinger equation for the antimatter fields of the observer's universe. In this case, we have ended up in the description of two universes, both expanding but one filled with matter and the other filled with antimatter 19 . The two interpretations raised here for the linear combination of incoming and outgoing wave look very similar to the two interpretations made in QED of a electron-positron pair (see, Fig. 6). In Sec. III B 4, we shall interpret this combination as the creation of a universe-antiuniverse pair.
Let us now specify the Hamiltonian of the perturbation modes. If we restrict to small linear perturbations the different modes do not interact and H m turns out to be the sum of the Hamiltonian of a set of harmonic oscillators [6,28] with where, M (t) = a 3 (t), and [2] for the tensorial modes of the spacetime (x n ≡ d n ), and The Schrödinger equation (71) of the perturbation modes is then the Schrödinger equation of a set of uncoupled harmonic oscillators, whose general solution can be written as [2,6,29] where the function ψ n (t, x n ) is the wave function of a harmonic oscillator with time dependent mass and frequency, which can be written in terms of the wave function of a harmonic oscillator with constant mass and frequency [30][31][32][33]. The general solution of ψ n (t, x n ) can then be expanded in the basis of number eigenstates of the invariant representation, ψ N,n , as where c N are constants coefficients and the wave function of the invariant number state, ψ N,n , is given by [31,34] whereψ N (q, τ ) is the customary wave function of the harmonic oscillator, with H N (q) the Hermite polynomial of order N , q ≡ xn σ , and σ(t) is an auxiliary function that satisfies the non-linear equation [30,31] σ +Ṁ Mσ plus some boundary condition [35]. The interesting property of the invariant representation is that once the field 20 is in a number state of the invariant representation it remains in the same state along the entire evolution of the field. For instance, let us assume that the perturbation modes are in the vacuum state of the invariant representation, |0 = n |0 n . In that case, the mean value of the energy of the perturbations reads 21 where ω n is given by (78) and (77) for matter particles or spacetime gravitons, respectively. The sum in (85) diverges and some cut-off n max must be taken. The energy H m is the backreaction of the perturbations on the homogeneous and isotropic background that would induce a modification of the Friedmann equation (62), where, ρ ϕ = 1 2φ 2 +V (ϕ), is the energy of the homogeneous mode of the matter field and, H m ∝ a −1 , for the massless tensor modes or, H m ∝ m, for the perturbations of the scalar field.

E. Paradigms for the creation of the universe in quantum cosmology
In the preceding section we ended up with a Friedman equation that was corrected by the backreaction of the perturbation modes of the spacetime. We shall see in this section that such term may induce important consequences in the way in which the universes can be created. However, for historical reasons instead of considering the backreaction of the perturbation modes we shall consider the model of a massless scalar field conformally coupled to gravity and a cosmological constant, which eventually raises the same term in the Friedmann equation, so the two models effectively entail similar effects. Later on, we shall briefly comment on the nature of this term. The conformally coupled massless scalar field is the field used by Hartle and Hawking in Ref. [12] to describe the quantum state of the universe and the process of quantum creation. Besides, it will also help us to introduce different paradigms for the creation of the universe.
Therefore, let us consider the following action for the spacetime and the massless, i.e. V (ϕ) = 0, scalar field, where the last term represents the conformal coupling of the scalar field, and 4 R is the Ricci scalar. Then, with the change, χ = aϕ, and after an integration by parts the total action (87) can be written as [12] where H 2 0 can be a pure cosmological constant, H 2 0 = Λ/3, or the constant value of the potential of an auxiliary inflaton field (different from the scalar field ϕ that we are considering in (88)). In any case, it will be assumed that it is a constant. Now, the momenta conjugated to the scale factor and the conformally coupled massless field χ can be easily obtained and the Hamiltonian constraint associated to the action (88) reads which, taking into account that p 2 χ + χ 2 is nothing more than (twice) the energy of the scalar field [35], 2E, can also be written as 20 Here, the field is the Schrödinger wave function ψn. 21 In Ref. [28] it is shown that the distribution of a large number harmonic oscillators becomes highly peaked around its average value. with which is formally similar to the energy equation of a particle that propagates under the action of the potential (91). The first term of the l.h.s. of (90) would be the kinetic energy, the second term would be the potential energy and E would be the total energy 22 Quantum mechanically, the wave function of the universe, φ(a, χ), is the solution of the Wheeler-DeWitt equation associated to the Hamiltonian constraint (89), which can be solved by the method of the separation of variables. Making, φ(a, χ) = ξ(χ)ψ(a), the Wheeler-DeWitt equation (92) can be split into the two following equations, The first of these equations is the equation of a quantum harmonic oscillator with unit mass and frequency. It can be solved in terms of Hermite polynomials, H n (x), with a quantised energy given by On the other hand, Eq. (94) can be written as with U(a) given by (91). Eq. (97) is formally similar to the Schrödinger equation of a particle of energy E moving under the action of the potential U(a), see Fig. 7. For the value, E ∈ (0, U max ), where U max is the maximum value of the potential, we can distinguish three regions, two classically allowed regions (regions I and III in Fig. 7) separated by a classically forbidden region (region II). We have already analysed in the preceding section the behaviour of the wave function in the region I. There are incoming and outgoing waves that represent a contracting universe that shrinks to the minimum value of the scale factor a + (see, Fig. 7) and bounces (or it is reflected) becoming an expanding universe. In the other classically allowed region, region III, we shall see that there are also incident and reflected waves that represent a universe that is confined to oscillate between a = 0 and the maximum value a − . In the case of a particle moving under the action of a similar potential, the particle that is placed in the region III is classically confined to move in that region like the universes in our example. However, we know that quantum mechanically there is a non zero probability for the particle to tunnel though the quantum barrier and appear in the region I. Something similar happens with the universe. If E > 0, the small oscillating universe of region III can tunnel out through the Euclidean barrier appearing in region I as a new born universe. The universe is then said to be created from something. In the limiting case, E = 0, which was the case analysed by Hartle and Hawking [12,[36][37][38] and by Vilenkin [39][40][41] for the creation of the universe, there is no region III from which to tunnel out to region I. The universe appears in region I from a pure tunnelling phenomena (like the creation of particles from the quantum vacuum). In that case, the universe is said to be created from nothing 23 . Let us analyse the two cases separately.

Creation of the universe from nothing
Let us first analyse the creation of the universe from nothing, i.e. E = 0 in (90). In that case, in terms of the time derivative of the scale factor, p a = −aȧ, the Hamiltonian constraint (90) reduces to the Friedmann equation (42) already studied in Sec. II D 1,ȧ which yields the well-known solution, with, a 0 = H −1 0 . If we restrict ourselves to the 'expanding branch', t ≥ 0, the scale factor (99) represents a universe that starts expanding from the initial boundary Σ(a 0 ) until infinity. For values a < a 0 there is no real solution and the value a = a 0 constitutes a classical barrier for the universe (see, Fig. 5). However, one can perform a Wick rotation into Euclidean time, t = −iτ , in terms of which the Friedmann equation (98) becomes Now, the Euclidean equation (100) has the solution where, τ ∈ (− π 2H0 , 0). Let us notice that in Euclidean time, −dt 2 → +dτ 2 , the line element turns to be which is the line element of a 4-sphere of radius 1/H 0 embedded in a 5-dimensional flat Euclidean spacetime. The 'spatial section' starts expanding from a single point of the sphere (at τ = −π/2H 0 ) until it reaches the value H −1 0 at the Euclidean time, τ = 0 (see, Fig. 8). It is called a DeSitter instanton and it gives the maximum contribution to the tunnelling wave function (it is the extremal solution of the Euclidean action [12]). At the boundary hypersurface Σ(a 0 ) it appears in the Lorentzian region as a new born DeSitter universe that starts expanding exponentially. Quantum mechanically the situation is the following. The initial boundary Σ(a 0 ) separates the Lorentzian region of the spacetime (see, Fig. 8), where classical solutions are allowed, from the Euclidean region, where only tunnelling solutions can exist. The wave function that represents a classical spacetime in the Lorentzian region is given in terms of the oscillatory wave functions e +iS and e −iS , where S is given by (47) S(a, ϕ) = (H 2 a 2 − 1) In the Euclidean sector, a < a 0 , the wave function of the universe can be written in terms of the tunnelling wave functions, e I and e −I , where I is given by the integral (47) with ω(a) being replaced by |ω(a)|, i.e.
The picture is then similar to the problem of a wave-particle tunnelling through a quantum barrier. The oscillatory wave functions e ±iS can be seen as the incoming and the reflected waves that may represent a photon or another quantum mechanical particle. Classically, the boundary Σ(a 0 ) acts as a barrier that cannot be crossed (see, Fig. 9). Quantum mechanically, however, there is a non-null probability of penetrate into the barrier although the amplitude is exponentially suppressed in the Euclidean region. Analogously, one can say in cosmology that there is a non-zero probability for the universe to appear from nothing, i.e. from the Euclidean barrier of the spacetime. An essential difference is that here the tunnelling is not from another classically allowed region of the spacetime. It is therefore more similar to the creation of virtual particles from the quantum vacuum in a quantum field theory. In a quantum field theory, the pair of virtual particles can only exist a small amount of time compatible with the Heisenberg's uncertainty relations, otherwise the principle of energy conservation would be violated. In the universe, however, the energy is zero (the negative gravitational energy balances the positive energy of the matter fields) so the creation of the universe from nothing does not violate the conservation of energy. From the above reasoning it is clear that the name 'creation from nothing' does not refer to the absolute meaning of nothing, i.e. to something to which we can ascribe no properties. As we have seen, the Euclidean region of the spacetime has geometrical properties. In the standard literature (see, for instance, Refs. [7,8,36,37,39]), it usually refers to two meanings. One, perhaps the most consistent, is that the universe is created from a region of the spacetime where nothing real exists; in particular, there is no actual time (i.e. time measured by clocks). In that sense, there is nothing. Another sense with which the term 'nothing' is used in the creation of the universe is that the universe, in the paradigmatic case of the creation of a DeSitter spacetime from the Euclidean instanton (101), begins from a single non-singular point of the Euclidean 4-sphere (which is geometrically equivalent to any other point in the sphere). In that case, the single point is meant to be 'nothing'.
a. Vilenkin's vs. Hartle-Hawking's versions The general quantum state of the universe is therefore given in the region I by the linear combination The particular combination of WKB wave functions in the Lorentzian and in the Euclidean regions of the spacetime, i.e. the particular values of the constants A I,II and B I,II , depend on the boundary condition that we impose on the state of the universe. Hartle and Hawking [12,36,38] propose as the boundary condition that the path integral must be performed over compact Euclidean geometries that fit with the 'final' values 24 a 0 and ϕ in the boundary Σ(a 0 ) (see, Fig. 8). In the homogeneous and isotropic minisuperspace, it is equivalent to the conditions [15] a(τ 0 ) = 0 , da dτ The wave function obtained for the Euclidean sections is e −I , with It yields [7,12,15] φ II HH (a, ϕ) = From the matching conditions of the WKB method, the Hartle-Hawking wave function (109) implies a linear combination of oscillatory wave functions in the Lorentzian region (a > a 0 ) [7,12,15], On the other hand, Vilenkin [39][40][41] proposes as the boundary condition, in an analogy with the tunnelling process of quantum mechanics, that the only modes that survive the quantum barrier are the 'outgoing' modes, i.e. those that represent an expanding universe 25 . It means that in the region I the wave function of the universe is given by [7,15,40,42] φ I T (a, ϕ) = where A I (ϕ) is a normalisation 'constant' that can be found by imposing the regularity conditions [15], ∂φ/∂ϕ → 0 as a → 0. Then, A I (ϕ) = exp −1/3H 2 , for which φ T ∼ e − 1 2 a 2 , that is regular at a → 0 for any value of the scale factor. By following the same WKB procedure of matching conditions, we found in the Euclidean sector [15] φ II T (a, ϕ) = where I is given by (108). Except for the values of the scale factor close to a 0 , the second term in (112) is exponentially smaller than the first, so is usually neglected. Besides their conceptual meaning the main difference between the wave function of the two proposals is the different sign in the exponent of the pre-factor, exp −1/3H 2 , in the case of the tunnelling wave function, and exp 1/3H 2 , in the case of the no-boundary wave function. The probability measures are different in both cases. Because the similarity of the Wheeler-DeWitt equation with the Klein-Gordon equation, Vilenkin proposes to use the probability current 26 [42,43] which is conserved because in virtue of the Wheeler-DeWitt equation it satisfies, ∇ · J = 0. The Hartle-Hawking wave function is real so the probability measure (113) would yield zero. These authors propose instead to use the customary probability measure of quantum mechanics, With these two choices the probability for the creation of the universe reads [15] where the + sign is for the no-boundary wave function and the − sign for the tunnelling wave function. Thus, as we have already notice in Sec. II B, the no-boundary proposal seems to favour small values of the potential (P HH → ∞ for H(ϕ) → 0) and the tunnelling proposal seems to favour the creation of universe with a large value of the potential (P T → 0 for H(ϕ) → 0). Therefore, it is usually stated that Vilenkin's tunnelling condition fits better with the inflationary scenario [24,42].

Creation of the universe from something
Let us now analyse the case, E ∈ (0, U max ) in (90). The corresponding Friedman equation is obtained by substituting the value, p a = −aȧ, in (90). It yields We can see that the conformally massless scalar field can reproduce the effect of a radiation content of the universe (ρ ∼ a −4 , H m ∝ a −1 in (86)). The Friedman equation (116) can also be written aṡ with [25,44] For the value, 1 8H 2 0 = U max > E > 0, the two Lorentzian regions are located at, a > a + and a < a − , respectively. These two sectors represent two separated regions where the universe may exist. In between, there is the Euclidean sector, which is a classically forbidden region. We have therefore two different types of universes separated by a quantum barrier (see Fig. 7). In region III, the solution of (117) can be written as [35] where, ∆t = t − t 0 , witht FIG. 10. The universes of region III are cyclic universes that start from a big-bang singularity and end in a big-crunch one. In region I the universe effectively behaves like a closed DeSitter spacetime.
It represents a small universe that starts in a big-bang like singularity at t = 0, expands to the maximum value a − , at t =t 0 , and then it re-collapses to a big-crunch like singularity, at t = 2t 0 . For a value H 0 1, the evolution of this type of universes is like the evolution of a radiation dominated universe (see, Fig. 10). This kind of universes are called baby universes [45], which are typically associated with quantum fluctuations of the spacetime.
On the other hand, the solution of (116) in region I can be written as [35] with, ∆t ∈ (−∞, ∞). It represents a universe that contracts from infinity to the minimum value a + , reached at ∆t = 0, and then, it expands again to infinity (see, Fig. 10). This solution is essentially very similar to the closed DeSitter universe. In fact, it can continuously be transformed into the customary solution of the closed DeSitter spacetime in the limit E → 0, for which a − → 0 and a + → 1/H 0 . In between there is a tunnelling region, where the solution of the Euclidean version of the Friedman equation is the Euclidean instanton (102) with scale factor given by where a ± is given by (118), a ∈ (a − , a + ), and The Euclidean instanton with scale factor (122) connects the maximum expansion hypersurface Σ(a − ) of the baby universes of region I with the initial hypersurface Σ(a + ) of the large parent universe in region I (see, Fig. 11). It therefore connects the two regions and it also provides the first order contribution to the probability of crossing the quantum barrier and appearing in region I as a new born universe (see, Ref. [35]). The universe is then said to be created from something. a. Gott and Li: the universe is its own mother The question of whether the universe is created from nothing or from something seems to be rooted in the value of the energy E in the Friedman equation (90). If E = 0 the universe must be created from nothing and if E > 0 the universe must be created from something. However, the value E = 0 is controversial because from (96), E = n + 1/2, so the value E = 0 would violate the uncertainty principle of quantum mechanics [44]. In its original paper [12], Hartle and Hawking suggested that this term might be cancelled by some renormalisation procedure. However, Gott and Li [44] argue that there is no expectation for such an exact cancellation and in fact, Barvinsky and Kamenshchik [46][47][48] compute the renormalisation corrections and not only the energy term is not cancelled but new similar terms appear. Even more, we have seen that the backreaction of the perturbation modes of the spacetime produces a similar term in the Friedmann equation (their energy density is given by, H m /V ∝ a −4 , where H m ∝ a − 1 in (85) for the perturbations of the spacetime). Therefore, a term with E > 0 seems to be unavoidable. One may argue that it can be effectively small and thus it could be neglected. However, regardless the small it can be we have seen that the consequences for the creation of the universe are very important.
A conceptual problem arises however if the universe is created from 'something' because then, we need the existence of that 'something' prior to the creation of the universe, i.e. if the universe is created from the tunnelling of a quantum fluctuation of a pre-existing spacetime, then, one should explain how has it been created the first spacetime from which the rest of universes have subsequently been generated. Gott and Li give an apparently exotic although quite interesting explanation. They argue and show [44] that in region III there can exist closed temporal curves (CTC's). In that case, the spacetime fluctuations of a large parent spacetime can travel through a CTC and become the baby universe that 'tunnelled out' through the Euclidean barrier to give rise to the parent universe in an atemporal process in which terms like 'after' or 'before' becomes meaningless. They are only meaningful within the large parent regions of the spacetime where CTC's do not exist. Thus, according to these authors, the universe could be its own mother!.

Creation of universes in pairs
There is a way-out to this paradoxic explanation. Even with the value E > 0 in (116), there is still room for the universes to be created from nothing, i.e. from the Euclidean region without the need of a pre-existing spacetime. However, as we have seen, the process cannot be the one studied in Sect. II E 1 for the creation of a single universe from the single Euclidean instanton (101). Instead, one has to consider more elaborated instantons. For instance, one can consider the double Euclidean instanton that is formed by joining together two single Euclidean instantons through their minimal hypersurfaces Σ i (a − ) (see, Refs. [46,47] and Fig. 12). The result is the creation from nothing of a pair of entangled universes in the region I [49][50][51] (see, Fig. 12).
Let us notice that in terms of the same time variable one of the universes of the entangle pair is a contracting universe and the other is an expanding universe so the situation is very similar to the case of coexisting incoming and outgoing waves. The wave function φ(a, χ) can therefore be written as with φ ± (a, χ) given by [35] φ ± (a, where ω DS is the potential of the Wheeler-DeWitt equation for the DeSitter spacetime, and where ξ n (χ) are the eigenfunctions of the harmonic oscillator 27 . As we have already seen in the Sec. II D 2, the two newborn universes can also be interpreted as two expanding universes filled with matter and antimatter, respectively. Let us notice that if one inserts the wave functions (125) into the Wheeler-DeWitt equation (92), it is obtained at which is the time dependent Schrödinger equation provided that one identifies the (conformal) time variable of the background spacetime of the two universes, φ ± , by If we assume that the physical time variable is the variable measured by real clocks, which are made up of matter, and thus that it is the time that appears in the Schrödinger equation, then, in terms of the physical time of an observer in one of the universes the Schrödinger equation for the fields in the partner universe turns out to be the Schrödinger equation of a fieldφ that is CP conjugated with respect the field in the observer's universe (see, Sec. III). The two universes form then a universe-antiuniverse pair. The process can be compared with the creation of a electron-positron pair (see, Fig. 1 of Ref. [39], and Fig. 6), which can be seen as the creation of an electron moving backward in time and an electron moving forward in time. Here, the 'time' variable is the scale factor so 'moving forward in time' means an expanding universe and 'moving backward in time' means a contracting universe, so the creation of a contracting-expanding pair of universes can be paralleled with the creation of a electron-positron pair.

A. Historial review
There was a great excitation in the 80's with the formulation of the third quantisation formalism and the associated description of topology change in quantum gravity [27,45,[52][53][54][55][56][57]. At that time the cosmological paradigm was a 27 In the superposition (127) it should appear ξ ± n , with (ξ + n ) * = ξ − n . However, the eigenfunctions of the harmonic oscillator are real functions so, ξ + n = ξ − n ≡ ξn. universe in a decelerating expansion with a zero value of the cosmological constant. The third quantisation formalism, which was initially proposed [52] as an analogue to the second quantisation formalism of a quantum field theory, fitted well with the description of the quantum fluctuations of the spacetime. As a consequence of the interaction with the fluctuations of the spacetime, it turns out that the coupling constants become dynamical functions and this was seen as a possible explanation for the expected vanishing value of the cosmological constant.
On the other hand, general relativity is a (local) geometrical theory and therefore it does not account for the global topology of the spacetime manifold. However, it may perfectly happen that the foliation of the spacetime gives rise, at some given value t, a collection of simply connected spatial sections (see Fig. 13 Left). From the point of view of the evolution of the universe the topology change represented in Fig. 13 (Middle) can be seen as the creation of a universe (universe B) from the pre-existing one (universe A). In fact, the process is formally similar to the QED process depicted in Fig. 13 (Right), which represents the creation of a photon from the propagation of an electron. As we know from QED this kind of processes are better explained in the formalism of quantum field theory where we can define creation and annihilation operators of particles. Therefore, it seems reasonable to develop a field theoretical approach too to describe the creation (and annihilation) of the universe(s).
As an introductory example, let us consider the Wheeler-DeWitt equation of a closed DeSitter spacetime [see (44)], which can be written asφ where, φ = φ(a),φ ≡ dφ da , and with, H 2 = Λ/3. Clearly, (131) is the equation of a harmonic oscillator with the scale factor a playing the role of the time variable. One can assume then that (131) is the result of the variational principle of the action of a harmonic oscillator with time dependent frequency, from which one can obtain the conjugate momentum, P φ =φ, and construct a the corresponding Hamiltonian, The third quantisation procedure consists in promoting the variables φ and P φ to quantum operators in the usual way,φ → φ andP φ → −i ∂ φ , and describe the quantum state of the whole spacetime manifold by the use of a new wave function, Ψ(φ), constructed as FIG. 14. Analogy between the second and the third quantisation [45,59].
in the path integral approach, or via the Schrödinger equation, One can also define the ladder operators,b andb † , of this particular harmonic oscillator in terms of the operatorsφ andP φ , as usual, and construct the state of the whole spacetime manifold, whatever the topology it has 28 , in terms of the eigenstates of the number operator, |N, a , so The number states, |N, a , would represent N universes with a scale factor a and the state of the whole spacetime manifold, |Ψ , would be a quantum superposition of different number states (see, fn. 28). Classically, these N universes are disconnected and therefore one should just consider one of them as representing our universe and disregard the rest as being physically irrelevant. However, from the quantum mechanical standpoint new phenomena may appear as quantum correlations and other collective behaviour so it seems interesting, at least in principle, to analyse the quantum description of the whole many-universe state. a. Parent and baby universes: the hybrid action The third quantisation formalism was mainly applied in the 80-90's to the description of the quantum fluctuations of the spacetime. Let us notice that the quantum gravity of multiply connected spacetime manifolds can be applied to two well distinguished cases: one that accounts for the properties of large regions of the spacetime, called parent universes 29 , and another that focus on a local region of the spacetime where small pieces of length of the order of the Planck length can branch off and disconnect from the parent spacetime and become small baby universes [45]. In both cases, the spacetime manifold under study turns out to be non simply connected. However, in the 80-90's the idea of a multiverse was not seriously considered and the main problem at that time was to explain the supposedly zero value of the cosmological constant 30 .
Then, from the point of view of the third quantisation formalism, our universe can then be seen as a large parent universe propagating in a plasma of baby universes. The effects of the baby universes could be measured by their influence on the observable properties of the parent universe, and the most representative picture of the baby-parent universe interaction becomes the so-called hybrid action [45], in which the parent universe is described by a second quantised wave function, φ p (q a ), and the baby universes are described by the third quantised wave functions,φ b , i.e. the behaviour of the spacetime is assumed classical and its fluctuations are seen as small particles propagating in the parent spacetime. The total action is then given by 28 In general, a non simply connected manifold can be divided into N simply connected parts [58] and this N parts can be seen as N classically independent universes. 29 Typically, large regions of order of the Hubble length of our universe. 30 In the 80-90's, the paradigm was a universe in a non accelerated expansion.
where S 0 (a, ϕ) is the Einstein-Hilbert action of the homogenous and isotropic parent spacetime with scale factor a and matter field ϕ, S b (φ b ) is the third quantised action of the baby universes, and S I is the action of interaction, where the index i labels the different modes of the baby universe field (i.e. it labels different species of baby universes), and L i (t, x) is called the insertion operator at the nucleation event [45]. It defines the space-time points of the parent universe in which the baby universes effectively nucleate. Two main problems were addressed in the 80's using the third quantisation formalism: the dynamical value of the coupling constants and the loss of quantum coherence (decoherence) produced by the plasma of baby universes.
As an example of the former, let us consider a universe with a matter field ϕ(t) that is coupled to the spacetime through the interaction with the baby universes. The parent universe two point function becomes [45] G where S I is given by (139) (for simplicity, let us assume just one specie of baby universes). Now, suppose the baby universes are in a wave function eigenstate |α , withφ b |α = α|α , where α satisfies In that case, the function (140) becomes The effects of the baby universes is thus encoded in an addition of an ordinary potential to the second quantised action. The new term is dynamical in the sense that it must satisfy the dynamical equation (141). That was used as an argument for a possible mechanism for the vanishing value of the cosmological constant, which was the expected value at that time.
The other question addressed with the third quantisation formalism was the loss of quantum coherence of the matter fields caused by their propagation in the plasma of baby universes. Basically, the argument was the following [53]. Let us suppose the composite state between a matter field ϕ and the baby universes. Let |ϕ, n be the state in which the matter field is in the state |ϕ and there are n baby universes. Then, the initial state is where |ϕ in is the initial state of the matter field. If we do not measure the state of the baby universes and we therefore integrate out their quantum state from the composite state (144), then, the initial state is described by the density matrix [53] ρ in = |ϕ in ϕ in |.
After the interaction with the baby universes, the final state becomes a linear combination of the states of the fluctuations and the corresponding states of the matter field that come out from the interaction with the |n states. For simplicity, let us consider just two |n states, |0 and |1 . The composite state after the interaction would be where ϕ 0 and ϕ 1 are in general different and the linear combination in (146) can be weighted accordingly. Then, the reduced density matrix that describes the state of the matter field alone becomes after the interaction The field turns out to be in a statistical mixture of two states. The initial state was a pure state, i.e. a state of total information with zero entropy, S = 0. The final state, instead, becomes a mixed state with entropy, S > 0, so information (quantum coherence) has been lost. Coleman's argument was that the operators of the baby universes must be independent of the coordinates of the parent spacetime and thus, the coupling with the matter fields is independent of their evolution. In that case, the state of the field does not change along the time evolution because the states of the baby universes do not change in time 31 . However, counter arguments were also given for the loss of quantum coherence [60,61].

Geometrical structure of M
We have seen in Sec. II that the evolution of the universe can be seen as the time evolution of the 3-dimensional metric that is induced on the spatial hypersurfaces by the 4-dimensional metric that is the solution of the Einstein's equations. Therefore, the evolution of the universe is a trajectory in the space of Riemannian symmetric 3-metrics with components, h ab . Let us call it M . With the DeWitt metric (7), M becomes a metric space, where we can define the line element as However, not all of the h ab components are independent. A symmetric 3-metric has only 6 independent components so it turns out that M is isomorphic to R 6 . Thus, we can make the following choice for the coordinates 32 in M , in terms of which the line element (148) can be written where G AB is a 6-dimensional metric tensor that is related to the components of DeWitt's metric, G abcd . The signature of M is (−, +, +, +, +, +), whichis easy to checked for the case of the flat metric, h ab = δ ab , and because the signature remains invariant under a change of coordinates, it holds then for the general case too. Thus, DeWitt showed [10] that the 6-dimensional space M is indeed a 5 + 1 dimensional space with a 1 time-like dimension and an orthogonal 5 dimensional space-like subspace. As the coordinate of the time-like subspace it is appropriate to take the coordinate τ defined by [7, 10, 62] τ = 32 3 where, h = deth ij , which essentially represents the volume of an infinitesimal volume element of the spatial sections of the spacetime (V ∝ dx 3 √ h). The hypersurfaces of constant τ are the space-like sections of M , labelled byM [10]. Then, in terms of the variables, q µ = {τ,q A }, whereq A , with A = 1, . . . , 5, are the five coordinates inM , the line element (150) becomes where, h 2 0 = 3/32, and ds is the line element inM , with [10] The metric (153) is invariant under a conformal transformation of the metric and in particular it is invariant under the change, h ab → ξ(h)h ab , so it is convenient labelling the points ofM with the five independent components of the transformed metric,h 31 It means that the state (146) would actually be, |out = |ϕ 0 (|0 + |1 ) . In that case, when we trace out the state of the baby universes the state of the field remains unaffected in the initial state, ρ out = |ϕ 0 ϕ 0 | 32 We have followed the normalisation applied in Ref. [10]. which has unit determinant. Furthermore,M is noncompact and diffeomorphic to Euclidean 5-space [10], with Ricci tensorR AB given by 33R and scalar curvatureR with, k = −1.M is then an 'Einstein space' of constant negative curvature. In a homogeneous universe the time-like variable τ represents the volume of the spatial sections of the universe and the five coordinates,h ab , represent the shape of a unit volume. A fixed point inM represents therefore the evolution of a universe that scales the volume of the spatial sections without changing their shape, and lines of constant τ represent different shapes of the spatial sections of the universe with a fixed given volume. Thus, the line element (152) reveals the space M as a set of "nested" 5-dimensional submanifolds, all having the same intrinsic shape [10]. From that point of view, the 5-dimensional submanifoldM can be seen as a proper realisation of what is called the 'shape space' [63].
On the other hand, the space M with the metric (152) has the same formal structure of a Friedmann-Robertson-Walker spacetime with the hyperbolic 5-space H 5 as the 'spatial' section. In particular, it has the same formal structure than the Milne spacetime 34 . Therefore, one can find a set of coordinates (χ, θ, φ, ψ, ζ) inM in terms of which the metric (152) can be written as 35 where, χ ∈ [0, ∞), and dΩ 2 4 is the line element on the 4-sphere of unit radius dΩ 2 4 = dθ 2 + sin 2 θ dφ 2 + sin 2 φ(dψ 2 + sin 2 ψdζ 2 ) .
in terms of which the line element (157) becomes with, 0 < T < ∞ and 0 < R < ∞. The metric (160) is nothing more than the metric of a 6-dimensional Minkowski space, and the Milne space only cover the upper light cone (see Fig. 15). The interior of the lower light cone is covered by a time reversal change of coordinates, τ → −τ (let us notice that the metric (152) is invariant under this change). However, although the manifoldM is geodesically complete, the manifold M is not. The scalar curvature of M , presents a singular frontier of infinite curvature, located at τ = 0, where all geodesics in M eventually hit [10]. It means that the upper and the lower light cones of the 6-dimensional Minkowski space must be considered independently. They represent two time reversed copies of the universe, i.e. two universes related by a time reversal change of the time (volume) variable. We shall see that, quantum mechanically, this can be seen as a universe-antiuniverse pair.

Classical evolution of the universe
Let us now analyse the evolution of the universe in the space, M ≡ Riem(Σ). From a geometrical point of view, the evolution of the universe is the trajectory that extremizes the Einstein-Hilbert action (6), which can conveniently be written as 36 [7] 33 This result is corrected from the one given in Ref. [10] by a factor 1 2 , which is already noted in Ref. [62]. 34 For the Milne spacetime, see Ref. [64] 35 A rescale, χ → aχ, θ → aθ, . . ., has been made to absorb the constant a. 36 Assuming the value, N i = 0. where have made the rescale, G abcd → 1 32πG G abcd , with G the Newton's constant, and the potential terms of the Einstein-Hilbert action have been gathered in a mass term In term of the variables, q A = (τ,q A ), the action (162) can be written where, G AB , is given by (152) or (157). The Einstein-Hilbert action has been written in the form of (164) to make clear the formal resemblance with respect to the action of a particle that moves in the spacetime, with,ẋ µ = dx µ dλ , for which the trajectory is given by the geodesic equation, with The case of the universe is formally similar. The evolution of the universe is a trajectory in M . The only difference is the trajectory is not a geodesic 37 because the non constant potential, m 2 (h ab ). Instead, it is given bÿ 37 The fact that the trajectory is not a geodesic is not really determinant. In fact, using a generalisation of the Maupertuis principle [65,66], one can compute the metric where the trajectory of the universe is a geodesic. Let us consider the reparametrisation given by, dt = m 2 (h ab )dt and G abcd →G abcd = m 2 (h ab )G abcd . In that case, the action (162) turns out to be where, h ab = dh ab dt . In the superspace determined by the supermetricG abcd the evolution of the universe turns out to be a geodesic.
FIG. 16. The evolution of the universe can be seen as a trajectory in the superspace [4].
where,q A = dq A dt and 2V = m 2 (h ab ), and the Christoffel's symbols are defined analogously in terms of the metric components as In terms of the variables (τ, q A ), the equations (169) turn out to bë As we have seen in the preceding section, a geodesic in M eventually hits the singular frontier located at, τ = 0 (the zero volume hypersurface). However, because the potential term, which may also include the Lagrangian of the matter fields, the universe does not follow a geodesic in M and may thus avoid the singular frontier. The paradigmatic case is the closed DeSitter spacetime. In addition to (171-172), one can use the Hamiltonian constraint (11), which in terms of the (τ,q A ) coordinates reads, In the case that the right hand side of (172) is zero, it can be shown that [10],ṡ = α/τ 2 , with α a constant of integration 38 , and thenτ In the case for which, Λ 3 R, the potential is proportional to τ , m 2 = h 2 0 H 2 0 τ 2 , anḋ which solutions given by for, α = 0, and for, α = 0 (FRW spacetime), with, ∆t = h 0 ∆t. From (176-177), one can see that, as expected, for a universe stage in which, Λ 3 R, the expansion of the volume element is exponential.

Quantum field theory in M
As we have seen above, the third quantisation 39 procedure consists in promoting the field φ(h ab ) and its conjugate momentum to quantum operators. One can then pose another wave function, Ψ, as in (135-137), and work with the corresponding Schrödinger equation [see, (136)]. However, it turns out to be much more interesting to develop and study the quantum field theory (QFT) of the field φ(h ab ) propagating in the 6-dimensional space M .
Let us first notice that, in terms of the coordinates q A = {τ,q A } (see, Eq. (152)), the Hamiltonian constraint (11) can be written as where, G AB is the inverse of (152), and in m 2 (q, ϕ) we have also included the Hamiltonian of the matter fields (which for simplicity have not been considered so far) that can generically be encapsulated in a variable ϕ, with, m 2 g (q), given by [see, Eq. (163)] and where the dots indicate terms that contain spatial derivatives of the matter fields, which for simplicity we shall consider negligible. Under canonical quantisation of the momenta in the Hamiltonian constraint (207) one obtains the Wheeler-DeWitt equation, which with an appropriate choice of factor ordering can also be written as, where, φ(q, ϕ), is the so-called wave function of the universe [12], and where, G = det G AB , and we have used the customary definitions of the gradient and the divergence in a curved space, With these definitions, the Wheeler-DeWitt equation (182) can be obtained from the variational principle of the third quantised action Variation of (185) with respect to the wave function φ gives rise to the wave equation (182). Now, following the analogy with a QFT, we can define a conserved current in the superspace, from which it can easily be checked that, ∇ · J = 0, and the following inner product where d Σ is the future oriented surface element of the 5-dimensional spacelike subspace. The procedure of field quantisation consists in promoting the wave function φ(q) to an operator and expand it into modes that are orthonormal with respect to the inner product (187). Then [55] where φ n (q) is a complete set of orthonormal solutions of the Wheeler-DeWitt equation (182), the index n symbolises the particular set of quantum number associated to that state [55], andÂ † n andÂ n are respectively the creation and annihilation operators of modes φ n , satisfying the customary commutation relations, In terms of the variables (τ,q), the label n of the modes φ n in (188) can be associated to the 5-dimensional spacelike momentum of the particles that propagate in the space M . In particular, we have seen that M has the geometrical structure of a 5 + 1-dimensional Friedmann-Robertson-Walker universe so we can use this information to develop the quantisation of the field φ(q). In terms of the coordinates q A = (τ,q A ), the Laplace-Beltrami operator (183) can be written as, where, q , is the corresponding 5-dimensional Laplacian (that with q given by (183) with the 5-dimensional metric G AB instead of G AB and without the minus sign in the square roots). In conformal time, λ = ln τ , and with the rescale, φ(q) = e −2λφ (λ,q)the wave equation (182) However, the 'mass' of the field (261) is not a constant. Even considering just the geometrical degrees of freedom it continues being a non constant function of the components of the metric tensor h ij (or, equivalently, of the variables q A ) through the dependence on the 3 R curvature (see, (180)). In that case, the space M turns out to be a dispersive medium for the wave function of the universe. It does not invalidate the formalism but it becomes more complicated from a technical point of view. For that reason let us focus on the case for which, 3 R 2Λ, which on the other hand is a very plausible condition for the initial state of the universe 40 , and consider only the geometrical degrees of freedom plus the constant Λ. We can then assume the value 41 in the Wheeler-DeWitt equation (182). In that case, the mass (261) only depends on the time variable τ and we can perform the quantisation of the field φ in the customary way (see, for instance, Refs. [76,77]). Then, we can decompose the wave function of the universe φ(q) in normal modes as where, k = (k, j), and, Y k, j (q), are the eigenfunctions of the Laplacian defined on the 5-dimensional hyperboloid, which satisfy [78] q Y k, j (q) = −(k 2 + 4)Y k, j (q), (195) 40 Let us notice that the condition, 3 R 2Λ, does not assume that the universe is homogeneous. 41 The factor 2 has been introduced for later convenience. with, 0 < k < ∞, and j denotes the 4 indices that distinguish the four components of the generalisation of the angular momentum on the 4 sphere 42 . Thus, the wave equation (182) (i.e. the Wheeler-DeWitt equation) reduces to where, χ ≡ dχ dλ . One interesting thing is that the frequency squared of the oscillator (196), ω 2 k (λ) ≡ m 2 0 e 2λ + k 2 , is never negative. The other interesting thing is that (196) is readily solvable in terms of Bessel functions. With the customary normalisation condition we easily find two set of orthonormal modes given by [75] χ k (τ ) = 2 π sinh(πk)

Boundary conditions and the creation of the universes in pairs
In order to choose the particular set of modes, we have to impose some boundary condition. For this, we shall consider the multiverse as a really closed system so no external influence is expected to modify its state. Therefore, it seems appropriate to describe the state of the multiverse in a representation that is invariant under the evolution of the third quantised Hamiltonian, which is an extension of the invariant representation used in quantum mechanics [30][31][32][79][80][81][82][83][84][85][86].
An invariant representation can be given in terms of creation and annihilation operators,b k andb † k , defined as [83] b where,φ andp φ , are the operator version of the wave function and the conjugate momentum, respectively, and ξ k is a solution of the wave equation (191), or equivalently with (196), the orthonormality condition (197), which ensures the usual commutation relations, The operatorsb k andb † k in (253) are time dependent operators but the dependence is such that the the eigenstates of the corresponding number operator,N k =b † kb k , remains invariant under the action of the third quantised Hamiltonian. It means, for example, that once the multiverse is in the vacuum state of an invariant representation it remains in the same vacuum state irrespective of the internal histories of the connected pieces of the whole spacetime manifold. From this point of view, the multiverse does not evolve in a proper sense, although the time dependence of the vacuum state makes that the vacuum state at (conformal) time, λ 0 , is functionally different than the vacuum state at λ 1 , being both however the same vacuum state of the same invariant representation.
The conditions (253-197) do not fix the vacuum state. There is in fact an infinite number of solutions that fit with (253-197), each of which define a particular representation and the associated vacuum state. For instance, the modes (260-265) define two vacuum states, |0 k0−k and |0 k 0 −k , respectively. The modes (260) can be identify with the Hartle-Hawking no boundary condition. First, because they are regular at the Euclidean origin 43 , τ → iτ → 0. Second, because that for large values of the variable τ ,χ k (τ ) ∝ cos m 0 τ , which essentially matches with the result (23). On the other hand, Vilenkin's tunnelling wave function can be identified with the mode (265), because at large values of τ it reads, χ k (τ ) ∝ e −im0τ , which represent, in terms of the time variable t, an expanding universe [see, Secs II B and II C]. One can also follow the analysis made in Ref. [76] to conclude that the state |0 k0−k is the conformal vacuum state, and that the state |0 k 0 −k is the vacuum state of the 6-dimensional Minkowski space, M .
Anyway, these two set of modes are related by a Bogolyubov transformation, where (see, for instance, Ref. [76]) α k = e πk 2 sinh(πk) 1 2 , β k = e −πk 2 sinh(πk) 1 2 , with, |α k | 2 − |β k | 2 = 1. It means that the vacuum state of theχ k modes, |0 k0−k can be written as [75,77] with a number of universes in the no bar representation given by which corresponds to a thermal distribution with generalised temperature Then, one can state that in this case the Hartle-Hawking no-boundary version of the vacuum state is full of (Vilenkin's) universes (and antiuniverses) [75]. This result is very interesting because it implies that the consideration of universeantiuniverse pairs seems to be quite unavoidable. It is formally similar to what happens in the quantum field theory of a matter field in an isotropic background spacetime, where the isotropy of the space makes that the particles are created in pairs with opposite values of the field modes, k and −k (see (204)), and in the case of a complex field in particle-antiparticle pairs with opposite momenta. In the third quantisation formalism, the space-like subspaceM is homogeneous and isotropic. Therefore, if the potential of the Wheeler-DeWitt equation is also isotropic inM , i.e. invariant under rotations inM , the universes should be created in pairs with opposite values of the 5-dimensional k (≡k ab ) modes. This is not the most general case, but it is a quite plausible one provided that we assume a high value of the potential of the inflaton field, which can be identified at the initial stage of the universe with Λ, or equivalently, a small value of the spatial curvature 3 R of the newborn universe. In both cases, the potential term of the Wheeler-DeWitt equation can be approximated by (192) and small deviations can be treated as perturbations, which should not significantly violate the isotropy of the spaceM . This could be confirmed as well from a more geometrical point of view. Let us notice that the Milne spacetime can separately cover the interior of the upper and the lower light cones of the Minkowski spacetime. These two sections of the full light cone can be seen as the regions of the spacetime where propagate future oriented particles and past oriented particles, or antiparticles, which turn out to be entangled [87]. One would expect something similar in the case of the space M , which also covers the upper and lower half light cones of the 6 dimensional Minkowski space. These two regions would describe expanding and contracting universes (created in pairs as we have seen above), which are equivalent to the future and past oriented particles in Minkowski spacetime. Thus, much in a similar way as particles propagating backwards in time can be interpreted as antiparticles propagating forward in time [88], we have seen in Sec. II, and we will see it again in the next section, contracting universes can be seen as expanding antiuniverses with a time variable that is reversely related with respect to the time variable of the partner universe. It means that the fields that propagate in one of the two entangled universes appear, from the point of view of an hypothetical observer in the partner universe, as moving backward in time. This is an illusionary effect created by the relative definition of the time variables in the two universes, i.e. internal observers always define the fields that propagate in their universes as matter and the fields that propagate in the partner universe as antimatter. Furthermore, the value of each mode of the Fourier decomposition in (193) is proportional to the momentum conjugated to the components of the scaled metric tensor. It means that any change that is produced by the momentum associated to +k ab in the shape of the universe with metrich ab is being also produced in the shape of the partner universe with opposite sign, −k ab , so it is the parity of the two spatial sections is reversely related too and so it is the relative parity of the fields that propagate in the two universes. In the next section, we shall see that the fields that propagate in the two universes are also charge conjugated as a consequence of the reversely relation of their time variables. It turns out therefore that the field of the two universes is CP conjugated. One can then conclude that, quite generally, the universes of the multiverse are created in symmetric universe-antiuniverse pairs whose composite quantum state is also expected to be entangled [49,89,90].

Semiclassical regime
If one takes into account the matter fields, the total Hamiltonian constraint (207) can be written as, whereĤ G is the Hamiltonian operator that yields the WDW equation of the spacetime geometry alone (182) with Λ related to the constant part of the potential of the field that drives the inflationary period, 2Λ = 2V 0 ≡ H 2 0 , and H SM (H m in (261)) is the Hamiltonian operator of the matter fields, that essentially are the fields of the Standard Model (SM) with their corresponding potentials and interactions. Following the procedure described in Sec. II C, the wave function of the universe can be written as the product of two components, a wave function φ 0 that depends only on the gravitational degrees of freedom and the value of the constant Λ, and a wave function that contains all the dependence on the fields of the SM, collectively denoted by the variable, ϕ, i.e where the two signs have been introduced for later convenience and, φ + = (φ − ) * . The wave function φ 0 is the solution of the WDWE of the geometrical degrees of freedom, computed in the preceding section. In general, it can be written in the semiclassical approach as If one introduces the wave function (208) into the complete WDW equation and use the classical constraint (30) one obtains, at order 1 , the following equation (see, (33)) where ∇ is the gradient in M and the negative and the positive signs correspond, respectively, to φ + and φ − in (208). The Schrödinger equation for the matter fields is then obtained if one defines the (WKB) time parameter t through the condition, where, q α = (τ,q A ), andq A are the coordinates ofM given in (152). We have now two choices. Typically, it is chosen the positive sign in (211) for the spacetime represented by the wave function φ − 0 and the negative sign for the spacetime represented by the wave function φ + 0 . With these choice, the Schrödinger equation in the two branches turns out to be where it can now be written, ψ ± = ψ ± (t ± ; ϕ). From (211), one easily gets so the wave functions ψ ± represent two universes, one expanding and one contracting (recall that the variable τ is proportional to the volume of the space), which from (212) are both filled with matter. An alternative although equivalent interpretation is to choose the positive sign in (211) for both universes, i.e. t ≡ t + . In that case, both wave functions represent expanding universes but then the corresponding Schrödinger equations for the internal fields are given by respectively. The last of which can be written as, FIG. 17. The creation of universes in entangled pairs [89]. In order to obtain the correct value of the Schrödinger equation in the two universes, their physical time variables must be reversely related. In that case, particles moving in the symmetric universe look as they were moving backward in time so they are naturally identified with the antiparticles that are left in the observer's universe. The primordial matter-antimatter asymmetry observed in the context of a single universe would thus be restored in the multiverse. Particles and antiparticles do not collapse at the onset because the Euclidean gap that exists between the two newborn universes [89,91].
where we have used that, ψ * − (ϕ) = ψ + (φ). It is therefore the Schrödinger equation of a field that is charge conjugated with respect to the field given in (214). The wave functions φ + and φ − represent then two expanding universes but from the point of view of the same time variable one is filled with matter and the other with antimatter, having these two concepts always a relative meaning.
Let us focus on the wave function of the matter fields in one of the universes, say ψ + . If we consider that the modes of the field are decoupled, then, the Schrödinger equation for the scalar field ϕ + , which generically denotes any of the polarisations of the W ± and Z bosons, can be written as the product of the wave functions of the modes, i.e.
where ψ (k) is the solution of the Schrödinger equation (214) for each mode, whose general solution can be expressed in the basis of number eigenfunctions of the time dependent harmonic oscillator [92] The wave function in the time reversely symmetric universe, ψ − (t,φ), can be obtained from the relation ψ − (φ) = ψ * + (ϕ), so the eigenfunctions of the basis for the state of the boson fields in the symmetric universe turns out to be given by (284) with the replacements, t → −t and ϕ k →φ k . Thus, the field ϕ that represents the matter content of one of the universes is the charge conjugated of the fieldφ that represents the matter content of the partner universe. We have seen in the preceding section that their parity is algo reversely related soφ turns out to be CP conjugated field of the field ϕ. Thus, the matter content of one of the universes is the CP conjugated of the matter in the partner universe and they form thus a universe-antiuniverse pair. It does not necessarily mean that one of the universes is completely made up of matter and the other is made up of antimatter. In fact, the two universes can contain matter as well as antimatter but exactly in the opposite ratio so, from the global point of view, the total amount of matter in the two universes is balanced with the total amount of antimatter. C. Minisuperspace model

Geometrical structure of the minisuperspace
Let us now apply the third quantisation formalism to the case of the minisuperspace of homogeneous and isotropic metrics with small perturbations that represent the matter content of the universe. The formalism greatly simplifies and one can still obtain a clear picture of the scenario described by the third quantisation formalism. On the other hand, we have seen that the minisuperspace description of the universe, although not complete, is a good approximation for most of the evolution of the universe provided that the universe is created with a length scale of some orders of magnitude above from the Planck length. In that case, the small deviations from the homogeneity and the isotropy of the universe can be treated as perturbations described as particles propagating in the homogeneous and isotropic background.
Let us therefore consider the homogeneous and isotropic Friedmann-Robertson-Walker (FRW) metric as the background spacetime (36) where a(t) is the scale factor, and dΩ 2 3 is the line element on the three sphere 44 . We saw in Sec. II that the lapse function is not a dynamical variable so the only dynamical variable turns out to be the scale factor, a(t). In this case all the components of the spatial metric are fixed except for the value of the scale factor.M turns out to be then a 0-dimensional space, where the spatial sections of the universes are represented by single points and their evolution by (curved) lines in the 1 + 0 dimensional space M .
This picture can easily be extended by considering as well the homogeneous mode of some matter fields, represented by a set of scalar fields, ϕ(t) = (ϕ 1 (t), . . . , ϕ n (t)), minimally coupled to gravity. We will see that these fields enter as space-like variables in the configuration space. For simplicity, we shall consider only one single scalar field representing the matter of the universe so the configuration space, M , will be a 1 + 1 dimensional space. In addition, on can also consider the inhomogeneous modes of these fields so the total configuration space would be the 1 + n · ∞ dimensional space spanned by the variables, (a(t), ϕ 1,k (t.x), . . . , ϕ n,k (t, x)). For simplicity, we shall only consider 45 the homogeneous mode of a single scalar field, ϕ, and its inhomogeneities will be treated as a perturbation described by particles propagating in the spacetime. Therefore, by now let us consider the 1 + 1 dimensional configuration space M of coordinates, q A ≡ (a, ϕ).
The total action, i.e. the Einstein-Hilbert action of gravity plus the action of the scalar field, given by (37), can be written as where, q A = (a, ϕ), with the supermetric G abcd in (13) given now by [7] G AB = diag(−a, a 3 ), from which one can clearly see that the scale factor (i.e. the first component) is a time-like variable and the scalar field (the second component) is a space-like variable. The potential term, V(q) in (222), reads The first term in (220) comes from the closed geometry of the three space, and V (ϕ) is the potential of the scalar field. The case of a spacetime with a cosmological constant, Λ, is implicitly included if we consider a constant value of the potential of the scalar field, V (ϕ) = Λ/6. As we showed in Sec. III B 2, the evolution of the universe can be seen as a parametrised trajectory of the superspace with the variable τ ∝ √ h formally playing the role of a time variable. In the case of the minisuperspace, √ h = a 3 ; but it is interesting to change to conformal scale factor, α = ln a, in terms of which the metric G AB turns out to be conformal to the 2-dimensional Minkowski space, and the action (222) can be written as, where, e 3α = a 3 , is essentially the volume of the spatial sections and, H 2 (ϕ) = 2V (ϕ), is the Hubble function. From the signature of (221), it can be seen that the scale factor formally plays the role of the time variable and the matter field(s) the role of the space-like component(s), and the minisupermetric (221) provides the minisuperspace with a complete metric structure with a line element given by It also allows us to define the usual machinery of a geometric manifold. For instance, we can defined the Christoffel symbols associated to the minisupermetric G AB , defined as usual by which in terms of the variables (a, ϕ) the non zero values are or in terms of the variables (α, ϕ), In any case, we could compute other geometrical properties of the minisuperspace like the corresponding Riemann tensor, the curvature scalar, etc (see, Ref. [10]). From the geometrical point of view the evolution of the universe can be seen as a trajectory in the minisuperspace (see, Fig. 18), with a(t) and ϕ(t) being the parametric coordinates of the universe along the worldline of the universe, and the time variable t is the parameter that parametrises the trajectory. From that point of view it is easy to see that the evolution of the universe, i.e. the trajectory of the universe in the minisuperspace, cannot depend on the particular choice of time variable, i.e. the trajectory must be independent of the parametrisation used to describe it.
However, because the presence of the potential V(a, ϕ) in the action (222), the trajectory of the universe along the minisuperspace manifold is not a geodesic. It is instead given by the equation which with the help of (225) yields the customary field equations (see, for instance, Refs. [7,24]) The fact that the curve (a(t), ϕ(t)) is not a geodesic is not a big deal. As we have said, the trajectory of the universe is invariant under reparametrisations of time, so we can make the following change of time variable where m is some constant. Now, if we also perform the following conformal transformation of the minisupermetric the action (222) becomes which is a similar action but with a constant potential. The new time variable,t, turns out to be the affine parameter of the minisuperspace geometrically described by the metric tensorG AB , and the trajectory of the universe in this minisuperspace is given by the geodesic equation Thus, the classical trajectory of the universe can equivalently be seen as either a geodesic of the minisuperspace geometrically determined by the minisupermetricG AB or a non geodesic of the minisuperspace geometrically determined by G AB . We can also define the momenta conjugated to the minisuperspace variables and the Hamiltonian constraint associated to the action (231) turns out to bẽ or in terms of the metric G AB and the time variable t, where for convenience we have written, m 2 ef (q) = V(q), with V(q) given by (220). It is worth noticing that the phase space does not change in the transformation {G AB , t} → {G AB ,t}, becausẽ where, p A = {p a , p ϕ } and q A ≡ {a, ϕ}, and the Hamiltonian constraints (234) and (235) are related by the inverse of the conformal transformation (230),G The field equations given either by (227) or by (232) are invariant under the reversal change in the time variable, t → −t. From the geometrical point of view it only changes the direction along which the curved is travelled, i.e. the direction of the tangent vector ∂ ∂t . It means that for any given solution a(t) and ϕ(t) one may also consider the symmetric solution, a(−t) and ϕ(−t).
In our case the momenta conjugated to the variables of the minisuperspace, given in (236), turn out to be in terms of which the Hamiltonian constraint (235) reads invariant but they turn out to be reversely changed, p a → −p a and p ϕ → −p ϕ . Nevertheless, they appear squared in the Hamiltonian constraint (239) so it is not affected by the change.
However, by conservation of the momenta one would expect that the cosmological solutions should come in symmetric pairs with opposite values of the associated momenta. From (238) and (239), it is easy to see that in terms of the cosmological time (N = 1) the two symmetric solutions are given by It clearly reminds to the solutions of the trajectory of a test particle moving in the spacetime [66]. For instance, in Minkowski spacetime 46 , the time component of the geodesics satisfies where µ is an affine parameter and, p t = ±E, with E the energy of the test particle. The two solutions are eventually associated to particles and antiparticles in a quantum field theory.
In the case of the universe the two solutions given in (240) also represent two universes: one universe moving forward in the scale factor component and the other moving backward in the scale factor component (see, Fig. 19). In the minisuperspace, however, moving forward in the scale factor component means evolving with an increasing value of the scale factor so the associated solution represents an expanding universe, and moving backward in the scale factor component means evolving with a decreasing value of the scale factor so the symmetric solution represents a contracting universe. Therefore, the two symmetric solutions form an expanding-contracting pair of universes (see, Fig. 19). However, we have already showed in previous sections that an expanding-contracting pair filled with matter can also be interpreted as two expanding universes, one of them filled with matter and the other filled with antimatter [25,91], i.e. it can be interpreted as a universe-antiuniverse pair [90,91].

Field quantisation of a FRW spacetime
As we have already seen, the procedure of third quantisation parallels that of a second quantisation in a curved spacetime (see Sec. III B 3). Now, the field is the wave function φ(a, ϕ) that satisfies the corresponding Wheeler-DeWitt equation, which is now seen as a wave equation. With the minisupermetric (221) in (182), the Wheeler-DeWitt equation turns out to be 46 A similar procedure can be followed in a curved spacetime. where, with, H 2 = V (ϕ 0 ), evaluated at the moment of the creation of the (inflationary) universe, where it can be approximated by a constant. In that case, following the procedure shown in Sec. III B 4 we can decompose the wave function φ(a, ϕ) in Fourier modes, where φ K (a) must satisfyφ with [49,89] Let us notice that the inner product turns out to be here by (187) with [93], dΣ A = n A dΣ, where n A = (a − 1 2 , 0) is a timelike unit vector, and dΣ = dϕ, which defines the orthogonal hypersurfaces (one dimensional curves) of constant a. It then becomes [49,93] We can now define the operator version of the field,φ, and write it aŝ whereÂ + K andÂ − K are the creation and annihilation operators, respectively, of universes with momentum K conjugated to the scalar field; and the modes are normalised according to the condition We can now define the ground state of the invariant representation,Â + K andÂ − K , by where |0 K I (|0 −K I ) is the state annihilated by the operatorÂ − K (Â − −K ). An excited state, i.e. a state representing different number of universes with momenta K 1 , K 2 , . . ., is then given by [93] |m K1 , n K2 , . . . = 1 √ m!n! . . .
which represents m universes in the mode K 1 , n universes in the mode K 2 , etc. In the case of a field that propagates in a homogeneous and isotropic spacetime the value of the mode k represents the value of the spatial momentum of the particle. In a homogeneous and isotropic minisuperspace the value of the mode K labels the eigenvalues of the momentum conjugated to the scalar field ϕ, which formally plays the role of a spacelike variable in the minisuperspace.
In that case, the values K 1 , K 2 , . . ., in (251) label the different initial values of the time derivatives of the scalar field in the universes. Thus, the state (251) represents m universes with a scalar field withφ ∝ K 1 , n universes with a scalar field withφ ∼ K 2 , etc. They represent different energies of the matter fields, which would correspond to different number of particles in the universes. The general quantum state of the field φ, which represents the quantum state of the spacetime and the matter fields, all together, is then given by |φ = m,n,...
FIG. 20. Left: in a quantum field theory the field is described in terms of particles that follow with the highest probability the classical trajectories given by the geodesics with however some uncertainties in their positions. Right: the wave function that describes the quantum state of the spacetime and the matter fields, all together, can be seen as a another field, say a super-field, that propagates in the minisuperspace. The universes can then be seen as 'test' particles following classical trajectories in the minisuperspace with quantum uncertainties given by the Schrödinger equation of their matter fields.
which represents therefore the quantum state of the multiverse [59] in the model of the minisuperspace that we are considering.
Here it follows the subtle subject of the boundary conditions in quantum cosmology. From a QFT we know that the vacuum state of a given representation may contain a certain number of particles of another representation, so the question is then which representation is the appropriate one. We have already imposed that the representation of the field that represents the state of the multiverse should be an invariant representation because, in that case, once that field is in a given state it will then remain in the same state along the entire evolution of the universe. Furthermore, one would expect that the field would be in the ground state of such an invariant representation provided that we assume that no external force is exciting the state of the multiverse. However, that condition does not completely fix the state of the field φ because there are many different invariant representations. In general, an invariant representation, A + K andÂ − K , can be defined as [83] where,φ andp φ , are the operator version of the wave function and the conjugate momentum in the Schrödinger picture, respectively, and v K is a solution of the wave equation (245)  One might say that the Hartle-Hawking boundary condition has a more fundamental character because it is rooted on a more ontological reasoning. It essentially rests on the idea that the boundary conditions of the universe are that it has no boundary [38], i.e. that the universe, and therefore the multiverse as well, comes from no prior configuration of the space. In that case, it seems consistent to impose or to assume that the multiverse is always in the ground state of the Hartle-Hawking invariant representation. However, single universes are better represented by the representation obtained by imposing the Vilenkin's tunnelling condition. In fact, this boundary condition is specifically imposed to assure that it describes single universes created in the Lorentzian region of the (mini)superspace (see Sec. II E 1). In that case, as we have seen in Sec. III B 4, it turns out that the multiverse is full of Vilenkin's universes [49], which due to the isotropy of the superspace should come in universe-antiuniverse pairs [75] (see, Sec. III B 4, and the next section). Thus, it seems quite unavoidable to assume that our universe has been created in an entangled pair. 3. Reheating and the matter-antimatter content of the entangled universe In Sec. III B 5 we have seen that the matter fields of the two universes of an entangled pair are CP reversely related and the two universes form thus a unvierse-antiuniverse pair. Let us now apply the same semiclassical formulation to the period after inflation called (p)reheating (see, for instance, Ref. [94]), where the inflaton field 47 , χ, eventually decays into the particles of the Standard Model (SM). In that period, the spacetime can largely be considered homogeneous and isotropic and the inhomogeneities of both the matter fields and the spacetime can be analysed as small perturbations propagating in a homogeneous and isotropic background.
We are not going to repeat the development of Sec. III B 5 but only to present a particular and detailed example that will show the consequences of the complex conjugated relation between these two wave functions. It is worth noticing that the CP conjugated relation between the matter fields of the universe-antiuniverse pair is based on the fundamental considerations described in Sec. III B 5 and it is therefore independent of the model chosen for the reheating scenario after inflation, so similar steps can be followed in any other reheating scenario. For concreteness, we shall describe this period in the appealing model of the Higgs-inflaton [21,22], in which the field that drives the inflationary expansion of the space decays after inflation into the particles of the Standard Model (SM). The idea rest on the form of the potential of the Higgs-inflaton field. At high energy scales, during the first stages of the inflationary period, the functional form of the potential can be approximated by an exponential and it can thus drive inflation. When the field has rolled down the exponential slope of the potential it finds a minimum around which it starts oscillating. Then, inflation ends and the Higgs-inflaton field behaves like the rest of fields of the SM, with interactions that allow the decays of the Higgs-inflaton into the particles of the SM (see below). Finally, in the low energy regime, the functional form of the potential can be approximated by the customary double-well potential of the Higgs that gives the expected masses to the particles of the SM [21,22]. Therefore, after the inflationary period the potential of the inflaton field, V (χ), cannot be longer considered a constant. However, we have seen in Sec. II D 2 that this does not introduce a big qualitative change. The solutions of the Wheeler-DeWitt equation are still expected to come in pairs, φ = φ + + φ − , with conjugate complex phases that are the solutions of the Hamilton-Jacobi equation (59). In the semiclassical regime, they are given by (64), i.e. φ ± (a; χ, ϕ) = ∆(a)e ± i S(a) ψ ± (a; χ, ϕ), where, χ, is the inflaton field and ϕ collectively denotes all the fields of the SM. Following the development of Sec. II C, the complex phase in (255) determines the dynamics of the homogeneous and isotropic background spacetime and the wave functions of the matter fields in the two universes, ψ ± (χ, ϕ), satisfy the Schrödinger equation of two sets of CP conjugated fields. They are related by the condition, ψ * − (χ, ϕ) = ψ + (χ,φ). As we have said, at the end of the inflationary period the Higgs-inflaton field χ has slow rolled down the potential and it then approaches the minimum of the potential located at χ m , for which V (χ m ) = 0. The expansion rate of the spacetime also slows down and the field starts oscillating around the minimum like a weakly damped harmonic oscillator with mass, m 2 = V (χ m ). The total Hamiltonian constraint can be written during this period as (207), with a gravitational part given by where B(a) contains the backreaction of the Higgs-field and eventually the backreaction of the rest of fields of the SM that will contribute to the dynamics of the background spacetime. It may also contain some residual constant term, which is expected to be subdominant at least until the advent of the dark energy period. The Hamiltonian of the Higgs-SM sector, H SM in (207), can now be written as with where [22], M 2 = λM 2 P /3ξ 2 , with ξ a coupling constant of the theory, p χ is the the momentum conjugated to the Higgs-inflaton field, p χ = a 3χ , and ∆V (χ) contains high order correction terms that can be neglected in a first 47 We shall use now the variable χ to represent the inflaton field and leave the variable ϕ to represent collectively the rest of fields of the SM.
approach [22]. The interactions between the Higgs and the matter and gauge fields of the SM have been included in the Hamiltonian H SM . The Klein-Gordon equation of the Higgs field can then be written, where we have assumed that the Higgs is essentially in the zero mode. For instance, for a power-law evolution of the background spacetime, a(t) ∝ t p , Eq. (259) is a Bessel equation that can be solved analytically. With the appropriate boundary conditions, and assuming M t 1 and p ≈ 2/3, it can be written as [22] where, χ(t = 0) = χ end , is the value of the Higgs field at the end of the inflationary period, which coincides with the beginning of the appearance of the (p)reheating mechanisms (t = 0). Different channels can now be considered for the decaying of the Higgs field into the particles of the SM (see, Ref. [22,95] for the details). It turns out that the perturbative decay of the Higgs field is only effective when the amplitude of the Higgs is below a critical value that depends on the mass of the final particles. This, together with the time dependence of the decay rate of the Higgs into the particles of the SM makes that the Higgs needs to oscillate a large number of times before decaying into the massive gauge bosons and fermions and much more times to decay into the less massive fermions, so the perturbative decay becomes ineffective during the first oscillations of the Higgs. In that period, the most effective channel turns out to be the parametric resonance [22,95]. However, this channel is enhanced by the effect of Bose stimulation so the production of fermions through this channel is highly restricted. These will be mainly produced later on through the perturbative channel or through the subsequent decay of the intermediate bosons into fermions.
Therefore, for the purpose of the present analysis, it is enough to focus on the production of the intermediate gauge bosons, W ± and Z. In the customary SSB mechanism the fields of the SM acquire a constant value of their masses. However, during the reheating period the potential still depends on the value of the Higgs field, χ, and thus the mass acquired from the interaction with the Higgs depends on its value. In that period, it can be approximated by [22] where g 2 is the coupling of the intermediate gauge bosons, θ W is the weak mixing angle, and y f are the Yukawa couplings of the fermion sector [22]. Eventually, after the period of reheating, in the low energy regime, the potential takes the customary form of a double well potential and the masses of the particles of the SM become the customary ones [22]. Thus, in the low energy limit the Higgs-inflationary scenario is indistinguishable from the Higgs scenario of particle physics, as expected. However, it is in this mid-energy regime, during the reheating period, when the basic components of matter are created in the two universes and the one in which we are mainly interested now. On the other hand, the quantisation of the intermediate gauge bosons W ± and Z follows as usual, by decomposing them into normal modes asφ where, ϕ ≡ W ± , Z, andb k andb † k are the annihilation and creation operators. As we have seen in Sec. II D 2, the inhomogeneous modes of the matter fields can be treated as particles propagating in the background spacetime. The mode amplitude ϕ k (t) satisfiesφ with (see, (78)) where m ϕ is now given by (261) with the value of the Higgs given in (260). In terms of the of the number of times that it crosses zero, j = M t π , the field (260) can be written as where a = a(t) is the scale factor of the background spacetime and ω k = ω k (t) is given by (264). The general solution of the Schrödinger equation (282) can then be expressed in the basis of number eigenfunctions given by [92] ψ +,N k (t, ϕ k ) = e −i(N +1/2)τ √ 2 N N !π 1/4 √ σ e − Ω 2 ϕ 2 k H N (ϕ k /σ), where H N (x) is the Hermite polynomial of degree, N ≡ N k , which is the number occupation of the mode k, τ = τ (t) is given by the function Ω is given by and σ is a real function that satisfies the auxiliary equation [92] σ + 3ȧ aσ The wave function in the time reversely symmetric universe is given by (see, Sec. III B 5), ψ + (t,φ), so the eigenfunctions of the basis for the state of the boson fields in the symmetric universe turns out to be given by (284) with the replacement, ϕ k →φ k . Therefore, if the scalar field ϕ represents the boson field W − in one of the universes, then,φ represents the boson field,W − = W + , in the symmetric universe 49 . The decay of the Higgs into the boson W + and W − can then be produced separately in the two symmetric universes. Then, one can make the hypothesis that the intermediate gauge boson W + and W − are created in different universes, or at least at different rates in the two universes, without violating the global matter-antimatter asymmetry, an appealing scenario that is also suggested in [71,96]. It is not mandatory that the asymmetry is complete but a small asymmetry in the decay of the Higgs into the W + and W − bosons in the two universes would eventually derive into an asymmetry in the production of fermions in the two universes due to the different decays of the W ± bosons into fermions (see, (277-279)). In the universe in which the boson W + predominates there would be an excess of the up quark with respect to the up antiquark, and accordingly, there would be an excess of protons over antiprotons, and matter would therefore dominate over antimatter. From the global picture of the two correlated universes the total amount of matter is always balanced with the total amount of antimatter so there is no global matter-antimatter asymmetry. It is worth noticing that the creation of a universe-antiuniverse pair does not assure that the content of one of the universes is completely matter and the content of the partner universe is completely antimatter. It is not therefore a mechanism for creating the matter-antimatter asymmetry but a mechanism to restore or explain the apparent asymmetry [91]. In a multiverse scenario, one may expect a whole range of matter-antimatter distributions along the pairs of universes in the multiverse. In some of them there would be the needed asymmetry to form matter, and therefore galaxies and planets like in our universe, without violating any physical law.

IV. OBSERVABLE EFFECTS OF QUANTUM COSMOLOGY
Testing the predictions of a theory with the observational data is a fundamental keystone of any physical proposal. However, it is not the unique consideration, theoretical consistency must also be taken into account and in fact it may help us to break through new paradigms. Well known examples in contemporary physics are the study of the unobservable black holes from the theoretical consistency of the perturbed motion of an observable companion, or the prediction of the charm quark from symmetry consistencies of the Standard Model of particles physics; not to talk about the unobserved 'dark matter' that is basically supported by consistency arguments. Furthermore, observability and falsifiability are not the same thing, as it is clearly argued 50 in Ref. [97] (see also, Ref. [98] for a recent review).
Nevertheless, any theory must eventually be tested. In principle, the effects of quantum gravity are expected to be relevant at a very small length, or equivalent to very large scale of energy, and that makes them to be hardly 49 Typically, ϕ would represent a linear combination of the W + and W + fields. In that case,φ would represent the corresponding conjugated combination. 50 Tegmark poses the following example: a theory stating that there are 666 parallel universes, all of which are devoid of oxygen, makes the testable prediction that we should observe no oxygen here, and is therefore ruled out by observation, cfr. Ref. [97], p. 105.
Perhaps the application to different inflationary models or the expected advances in the field of astronomical detectors and associated space missions, with the detection and analysis of gravitational waves or the 'cosmic neutrino background' (CNB), might make directly testable in the future the deviations from classicality predicted from quantum cosmology. Nevertheless, even though they may be difficult to be observed these effects may have important conceptual consequences. For instance, we have seen in Sec. II E that an energy term like the one produced by the backreaction in the Friedman equation (289) might drastically change the way in which the universes can be created.
On the other hand, the creation of the universe in entangled pairs 51 may also add new features to be tested in the future. For instance, in an entangled universe the fields of the matter content in the two universe stop being a vacuum state. If one computes the state of the matter field in one single universe of the entangled pair by tracing out from the composite state the degrees of freedom of the matter in the partner universe, then, the resulting state turns out to be a quasi-thermal state with a temperature that depends on the degree of entanglement (which eventually depends on the size of the universe). In Ref. [89], it is computed the ratio between the fluctuations of the perturbation modes of a field that is initially in a thermal state and those corresponding to a initial vacuum state, yielding [89] δφ th n δφ I with, where, L ph , is the physical wave length and H −1 is the distance to the Hubble horizon. The large modes (x 1) are in the vacuum state and then, δφ th n ≈ δφ I n . However, the departure is significant for the horizon modes, x ∼ 1. This would be a distinctive effect of the creation of the universes in entangled pairs and it should leave an observable imprint in the properties of the CMB.

V. CONCLUSIONS
Quantum geometrodynamics provides us with a consistent framework for the quantum description of the universe in terms of a wave function that contains, at least in principle, all the information about both the spacetime and the matter fields that propagate therein. In the semiclassical regime the complex phase of the wave function contains the information about the dynamics of the background spacetime and the wave function of the matter fields satisfy a Schrödinger equation that depends on the geometry of the subjacent spacetime. The dynamics of the spacetime turns out to be invariant under the complex conjugation of the semiclassical wave function. That gives rise two different solutions that have been typically interpreted as representing the expanding and the contracting branches of the universe, both filled with matter. We have seen that a different, more consistent interpretation is that the two solutions represent expanding universes with their matter contents being CP reversely related, so from the point of view of an internal observer of any of the universes the partner universe is always made up of antimatter, having therefore these two terms a relative meaning.
On the other hand, the Wheeler-DeWitt equation can be seen as the wave equation of a field that propagates in the space of Riemannian 3-dimensional geometries, M , where we can describe the evolution of the universe as a trajectory parametrised by a parameter that we can call time. The quantum mechanical counterpart is a quantum field that, following the customary approach of a quantum field theory, can be expressed in terms of creation and annihilation operators that satisfy the usual commutation relations. These operators represent the creation and the annihilation of modes for the spatial sections of the universe. Thus, the third quantisation formalism allows us to describe the quantum state of a whole spacetime manifold that can in general be a disconnected collection of simply connected manifolds. It is thus an appropriate framework to describe a multiverse scenario, in which the most natural boundary condition turns out to be that the field that represents the whole spacetime manifold remains in the ground state of an invariant representation along the entire history of the universes. This boundary condition implements the idea that the multiverse is the true isolated system and therefore no external interaction may excite its quantum state. However, as we have seen, this invariant boundary condition does not fix completely the quantum representation of the universe and, in fact, the ground state of one invariant representation is in general full of pairs of universes in another invariant representation. For instance, we have seen that the Hartle-Hawking no boundary condition can be seen as more fundamental because it is based on a more ontological argument. In that case, one may assume that the field that represents the multiverse is in the ground state of the invariant Hartle-Hawking representation. Because the invariance of the invariant representation the field remains then in this ground state irrespective of the evolution of the universes. However, in terms of the invariant representation associated to the Vilenkin's tunnelling boundary condition, which represents the state of single universes, it turns out that the ground state of the Hartle-Hawking no boundary state is full of Vilenkin's pairs of universes. It means that the ground state of the multiverse is full of pairs of universes whose matter contents turn out to be CP conjugated. The charge conjugation comes from the complex conjugation relation between the Schrödinger wave functions of the matter fields in the two universes, and the reversely parity relation comes from the opposite signs of the momentum conjugated to the geometrical variables of the spatial sections of the spacetime. The two universes form thus a universe-antiuniverse pair.
We have seen all these features in a general model but particularly in the model of a homogeneous and isotropic spacetime with particles and matter fields propagating therein, where explicit examples can be analysed. In particular, we have seen that there are three main paradigms for the creation of the universe in quantum cosmology. The universes can be created from nothing, i.e. from no preexisting spacetime. However, it requires of a precise fine tuning that seems to be quite unnatural. The other possibility is then that the universe is created from the quantum fluctuation of the spacetime of a preexisting spacetime. However, to be then consistent, one should also give an explanation for the creation of the first spacetime. We have seen that the Gott and Li's explanation is that a vacuum fluctuation of our spacetime can travel back along one of the closed temporal curves that are allowed to exist in the vacuum state of the gravitational field to become the seminal spacetime from which our spacetime has been created. From that point of view, the universe would be the mother of itself. There is yet another possibility that avoids this paradoxical conclusion. The universes can be created from nothing, i.e. with no need of any prior spacetime, but they must then be created in pairs, from double Euclidean instantons. That turns out to be the most natural and self-consistent way in which the universes can be quantum mechanically created.
We have also seen in a very specific model that the creation of matter after the period of inflation would be correlated in the two universes of an entangled pair. The decay of the inflaton field into the particles of the Standard Model after inflation in the so-called reheating period is produced in such a way that the matter and antimatter of the two universes is perfectly balanced. The matter-antimatter asymmetry observed in our universe would only be therefore an apparent asymmetry and, in fact, it might be considered as an evidence of the existence of an entangled companion of our universe that would contain the amount of antimatter that is left in our universe. However, the observational test of this and other quantum cosmological hypothesis seems to be still far from the current state of observation as the effects of an entangled partner would be in the domain of quantum gravity or at least in a pre-inflationary stage of the universe. Perhaps in the next future the advances in the detection of primordial gravitational waves or a possible cosmic neutrino background may shed some light into these intriguing questions.