Emergence of multiverse in third quantized varying constants cosmologies

TL;DR

This paper applies third quantization to a varying constants cosmology, demonstrating the emergence of a multiverse from nothing, with universes described by a Bose-Einstein distribution, addressing initial singularity issues.

Contribution

It introduces a third quantization approach to varying constants cosmology, revealing multiverse creation from nothing and characterizing the universe family distribution.

Findings

Multiverse can emerge from nothing in this framework.

Universes follow a Bose-Einstein distribution.

The model addresses initial singularity problems.

Abstract

Although the standard cosmological model explains most of the observed phenomena it still struggles with the problem of initial singularity. An interesting scenario in which the problem of the initial singularity is somehow circumvented was proposed in the context of string theory where the canonical quantisation procedure was additionally applied. A similar effect can be achieved in the context of the canonically quantized theory with varying speed of light and varying gravitational constant where both quantities are represented by non-minimally coupled scalar fields. Such theory contains both the pre-big-bang contracting phase and the post-big-bang expanding phase and predicts non-vanishing probability of the transition from the former to the latter phase. In this paper we quantize such a theory once again by applying the third quantization scheme and show that the resulting theory contains scenario in which the whole multiverse is created from nothing. The generated family of the universes is described by the Bose-Einstein distribution.

Equations

$$S=\int \sqrt{-g}\left(\frac{e^{\varphi }}{e^{\psi }}\right)\left[R+\Lambda +\omega ({\partial }_{\mu }\varphi {\partial }^{\mu }\varphi +{\partial }_{\mu }\psi {\partial }^{\mu }\psi )\right]d^{4}x,$$

$$\varphi =\frac{\beta }{\sqrt{2\omega }}+\frac{1}{2}\ln \delta ,\psi =\frac{\beta }{\sqrt{2\omega }}-\frac{1}{2}\ln \delta ,$$

$$S=\int \sqrt{-g}\left[\delta (R+\Lambda )+\frac{\omega }{2}\frac{{\partial }_{\mu }\delta {\partial }^{\mu }\delta }{\delta }+\delta {\partial }_{\mu }\beta {\partial }^{\mu }\beta \right]d^{4}x.$$

$$d{s}^2=-N^2(d{x}^0{)}^2+a^2\left(d{r}^2+r^{2}d{\Omega }^2\right),$$

$$S$$

$$a$$

$$\delta$$

$$\begin{array}{rl}{x}^0& =\frac{2}{\sqrt{(A^2-9)\Lambda }}\mathrm{arctanh}\left(e^{\sqrt{(A^2-9)\Lambda }{\bar{x}}^0}\right),\text{for xˉ0<0},\\ x^0& =\frac{2}{\sqrt{(A^2-9)\Lambda }}\mathrm{arctanh}\left(e^{-\sqrt{(A^2-9)\Lambda }{\bar{x}}^0}\right),\text{for xˉ0>0},\end{array}$$

$$S$$

$$H$$

$$B$$

$$I$$

$$J$$

$$\pi_{I}=\left\{\begin{array}[]{l l}\pi^{\infty}_{I}&\quad\text{collapsing pre-big-bang }\\ -\pi^{\infty}_{I}&\quad\text{expanding post-big-bang}\end{array}\right.$$

$$\pi_{I}=\left\{\begin{array}[]{l l}\pi^{\infty}_{I}e^{-I}&\quad\text{collapsing pre-big-bang}\\ -\pi^{\infty}_{I}e^{-I}&\quad\text{expanding post-big-bang}\end{array}\right.$$

$$\left\{\frac{1}{4}\left[\frac{1}{m}\left(\frac{{\partial }^2}{\partial I^2}-\frac{{\partial }^2}{\partial J^2}\right)+\frac{{\partial }^2}{\partial B^2}\right]-\bar{\Lambda }{e}^{-2I}\right\}\Phi =0$$

$$\Phi =\alpha (J)\gamma (B)\beta (I),$$

$$\alpha (J)=e^{i{k}_{1}rJ},\gamma (B)=e^{i2k_{2}B},\beta (I)=J_{-i{\pi }_I^{\infty }}({\pi }_I^{\infty }{e}^{-I}),$$

$${\tilde{V}}_0\Lambda =k_1^2+k_2^2.$$

$${\Phi }^{\prime \prime }-\Delta \Phi +m_{eff}^2\Phi =0,$$

$$\eta$$

$$x_1$$

$$x_2$$

$$\hat{\Phi }(\vec{x},\eta )=\frac{1}{\sqrt{2}}\int \frac{d^{2}k}{2\pi }[e^{i\vec{k}\cdot \vec{x}}{v}_k^{\ast }(\eta ){\hat{a}}_{\vec{k}}^{-}+e^{-i\vec{k}\cdot \vec{x}}{v}_k(\eta ){\hat{a}}_{\vec{k}}^{+}],$$

$$v_k(\eta {)}^{\prime \prime }+{\omega }_k(\eta {)}^2{v}_k(\eta )=0,$$

$$W(v_k(\eta ),v_k^{\ast }(\eta ))=2i,$$

$$[{\hat{a}}_{\vec{k}}^{-},{\hat{a}}_{{\vec{k}}^{\prime }}^{+}]$$

$$[{\hat{a}}_{\vec{k}}^{-},{\hat{a}}_{{\vec{k}}^{\prime }}^{-}]$$

$$[{\hat{a}}_{\vec{k}}^{+},{\hat{a}}_{{\vec{k}}^{\prime }}^{+}]$$

$$u_k(\eta )={\alpha }_k{v}_k(\eta )+{\beta }_k{v}_k^{\ast }(\eta ),$$

$$\hat{\Phi }(\vec{x},\eta )=\frac{1}{\sqrt{2}}\int \frac{d^{2}k}{2\pi }[e^{i\vec{k}\cdot \vec{x}}{u}_k^{\ast }(\eta ){\hat{b}}_{\vec{k}}^{-}+e^{-i\vec{k}\cdot \vec{x}}{u}_k(\eta ){\hat{b}}_{\vec{k}}^{+}],$$

$$[{\hat{b}}_{\vec{k}}^{-},{\hat{b}}_{{\vec{k}}^{\prime }}^{+}]$$

$$[{\hat{b}}_{\vec{k}}^{-},{\hat{b}}_{{\vec{k}}^{\prime }}^{-}]$$

$$[{\hat{b}}_{\vec{k}}^{+},{\hat{b}}_{{\vec{k}}^{\prime }}^{+}]$$

$${\hat{a}}_{\vec{k}}^{-}$$

Full text

1 Introduction

The standard cosmological model based on general theory of relativity and standard model of particle physics encounters a certain number of the fundamental problems. However, although it explains a vast majority of the observational data, it is plagued with huge fine-tuning problems. Some of those can be resolved in the standard inflationary paradigm Linde . Yet the classical theory of inflation is far from being complete since it does not provide initial conditions for the inflation to occur in a desired way. The incompleteness of this picture stretches to even more profound issues such as the occurrence of the initial singularity, the related questions about the beginning of time or the existence of epochs preceding big-bang. A consistent framework for investigating the above-mentioned questions is provided by string theory. The scenarios involving eras preceding big-bang singularity naturally arise in the cosmological models based on the string theory Gasperini . The cosmological scenarios based on the tree level low energy effective action classically include both pre-big-bang and post-big-bang evolution phases separated by a singularity in which the curvature and the coupling reach infinite values. Employing the quantum cosmology Wheeler-DeWitt approach to the description of the near singularity regime gives rise to the scenario in which the universe passes form the pre-big-bang to the post-big-bang phase in a process that can be viewed as reflection (scattering) of the Wheeler-DeWitt wave function on one dimensional exponential potential barrier in the minisuperspace Birth . A similar scenario in which the universe scatters quantum mechanically over the singularity to finally enter the standard post-big-bang expansion era arises in the context of gravity theories with varying speed of light (VSL) and varying gravitational constant Balcerzak . Many different VSL theories have been investigated in the literature so far, however, inventing a consistent picture in which the speed of light is allowed to vary encounters conceptual problems. One of these is the violation of the Lorentz invariance in case the speed of light is assumed to be dependent on the space-time coordinates Albrecht ; Barrow1 ; Magueijo1 ; Clayton ; Drummond ; Clayton2 . Violation of the Lorentz invariance forces in turn an introduction of the preferred reference frame in which the particular mathematical structure of an investigated VSL theory is formulated. An interesting realisation of the mentioned above concepts is the VSL theory proposed in Albrecht ; Barrow1 . The model additionally assumes that the degree of freedom representing speed of light is minimally coupled to the matter and the gravitational field in the preferred frame. Consequently, the equations of motion of such VSL gravity theory are identical to the standard Einstein equations with the speed of light being merely replaced by some time dependent functions. A different VSL gravity theory can be obtained by releasing the assumption of the minimal coupling Magueijo1 . The resulting set of equations of motion in such theory gains additional dynamical equation governing the evolution of the non-minimally coupled degree of freedom representing the speed of light. Another interesting example is VSL theory in which the speed of gravitons is assumed to be different from the speed of massless matter particles Clayton ; Drummond ; Clayton2 . This is achieved by introducing two different metrics. The first one describes the space-time geometry, the second one couples to the matter. The separate group of theories comprises VSL models which include corrections to the dispersion relation Rainbow . Such corrections become relevant for the energy scales comparable with the Planck scale. The resulting group velocity of light depends on the energy scale.

In paper Balcerzak it was shown that the near curvature singularity regime can be tackled within the framework of the theory which assumes that both the speed of light and the gravitational constant can vary. Both quantities are represented by scalar fields non-minimally coupled to the gravitational field in the preferred frame defined by the FLRW metric. This is different form the approach presented in Leszczynska where the similar regime is investigated in the context of the model in which the degrees of freedom representing varying speed of light and gravitational constant minimally couple to gravity and matter.

The so-called third quantization is based on the formal similarity between the Wheeler-DeWitt and the Klein-Gordon equations Strominger ; Robles . The role of the Klein-Gordon field is played by the wave function in the Wheeler-DeWitt equation, which as a result of the third quantization becomes an operator acting on the Hilbert space. The third quantization itself is completely analogous to the quantization of the Klein-Gordon field. The resulting Hilbert space is spanned by an orthonormal basis which elements represent occupation with universes which properties are determined by appropriate quantum numbers (these for example can be the momenta in the minisuperspace).

Our paper is organised as follows. In Sec. 2 we quote the main results of the paper Balcerzak which will be a starting point for the realization of the central task of our work consisting in elaboration of the scenario in which the whole multiverse emerges out of vacuum (“out of nothing”). In Sec. 3 we fulfill our assumed task by applying the third quantization procedure to the canonically quantized non-minimally coupled varying constants model and by showing that the resulting theory includes a scenario in which the family of universes described by Bose-Einstein distribution is created out of nothing.

2 Non-minimally coupled varying $c$ and $G$ theories

The classical action defining the considered varying $c$ and $G$ model can be obtained by replacing each constant in the ordinary Einstein-Hilbert action with a certain function of a dynamical scalar degree of freedom. In order to keep $c$ and $G$ positive during the cosmological evolution we link them to the newly introduced scalar degree of freedom $\phi(x^{\mu})$ and $\psi(x^{\mu})$ by the following exponentials $c^{3}=e^{\phi}$ and $G=e^{\psi}$. The resulting action reads Balcerzak :

[TABLE]

where $\omega$ is the parameter of the model. By application of the field redefinitions of the form:

[TABLE]

the action (1) can be transformed into the form of Brans-Dicke action:

[TABLE]

The assumed dependence of $c$ on space-time coordinates breaks the invariance of the action (1) or (3) under general transformations of coordinates. This means that we have to specify the reference frame in which the action given by (1) or (3) is assumed to describe our varying $c$ and $G$ model. A natural step is to identify the above-mentioned preferred reference frame with the cosmological frame defined by the flat FLRW metric:

[TABLE]

where $N$ is the lapse function and $a$ is the scale factor both depending on the parameter $x^{0}$ . Inserting the metric (4) into (3) gives the form of the action of our model in the cosmological frame:

[TABLE]

where $()^{\prime}\equiv\frac{\partial}{\partial x^{0}}$. A similar approach was introduced in papers Albrecht ; Barrow1 where the VSL action and the corresponding variational principle was formulated in the cosmological frame. The solution to the model given by action (5) expressed in the gauge determined by $N=a^{3}\delta$ is:

[TABLE]

where $A=\frac{1}{\sqrt{1-2\omega}}$, $M=\frac{3-A^{2}}{9-A^{2}}$, $W=\frac{2A^{2}}{9-A^{2}}$ and $D$ and $F$ are some integration constants and the variable $x^{0}$ is the following function of the rescaled proper time $\bar{x}^{0}$ defined with its differential $d\bar{x}^{0}\equiv c(\bar{x}^{0})d\tau$ with $\tau$ being here the proper time encountered by the comoving observer (for the detailed derivation of the solution given by (6), (7) and (8) see Balcerzak ):

[TABLE]

where as in Balcerzak we restrict our considerations to the models with $A^{2}>9$. Fig. 1 depicts qualitative behaviour of the scale factor $a$, the speed of light $c$ and the gravitational constant $G$ in the near high-curvature regime for the pre-big-bang $\bar{x}^{0}<0$ and post-big-bang $\bar{x}^{0}>0$ phase. We see that as the universe approaches the curvature singularity ($a\rightarrow 0$) the speed of light $c$ goes to infinity while the gravitational constant $G$ tends to zero. It means that the transition from pre-big-bang to post-big-bang phase occurs in the Newtonian limit. A similar behaviour of the scale factor occurs in the ekpyrotic Khoury ; Khoury2 and cyclic scenarios Steinhardt ; Steinhardt2 where both pre-big-bang and post-big-bang eras are separated by the curvature singularity induced by vanishing scale factor. A qualitatively different evolution is encountered in the case of pre-big-bang scenarios based on the low-energy effective action of the string theory where curvature singularity is reached during the phase of accelerated expansion Birth ; Gasperini .

Studying the quantum mechanical nature of the considered model requires analysis of the solutions of the Wheeler-DeWitt equation. In order to obtain the Wheeler-DeWitt equation corresponding to our model first we have to find the form of the hamiltonian. By application of the following field transformation $X=\ln(a\sqrt{\delta})$, $Y=\frac{1}{2A}\ln\delta$, $I=AY-3X$, $J=3Y-AX$ and $B=\sqrt{\tilde{V}_{0}}\beta$ the action (5) can be reduced to the following form:

[TABLE]

The hamiltonian corresponding to our model is given by:

[TABLE]

Here $\bar{\Lambda}=\tilde{V}_{0}\Lambda$ with $\tilde{V}_{0}=\frac{3V_{0}}{8\pi}$ while $\pi_{I}=-2mI^{\prime}$, and $\pi_{B}=-2B^{\prime}$ are the canonically conjugated momenta with $m=\frac{\tilde{V}_{0}}{9-A^{2}}$. The form of the hamiltonian proves that $\pi_{J}$ and $\pi_{B}$ are constant during the evolution. Therefore the classical evolution of the universe in the near curvature singularity regime can be reduced to the process of scattering of a particle on the exponential potential barrier. The solutions (6) and (7) expressed in term of the new variables $B$, $I$ and $J$ are given by:

[TABLE]

where $E$ and $P$ are some constants.

By analysing the classical solutions (11)-(13) we find that the near curvature singularity regime corresponds to the region located at $I\rightarrow\infty$. On the other hand the low-curvature regime for both the pre-big-bang and the post-big-bang branches is associated with the region located at $I\rightarrow-\infty$. Both regimes correspond to the particular asymptotic values of the momentum $\pi_{I}$. In the high-curvature (near big-bang) limit $I\rightarrow\infty$ we have

[TABLE]

while in the low-curvature (far away from big-bang) limit $I\rightarrow-\infty$ the momentum $\pi_{I}$ is given by:

[TABLE]

where $\pi^{\infty}_{I}\equiv 2\tilde{V}_{0}\sqrt{\frac{\Lambda}{A^{2}-9}}$ (again for the detailed derivation see Balcerzak ). In order to obtain the quantum theory we apply the Jordan quantization rules and replace the canonical momenta with the following operators: $\pi_{J}\rightarrow\hat{\pi}_{J}=-i\frac{\partial}{\partial J}$, $\pi_{I}\rightarrow\hat{\pi}_{I}=-i\frac{\partial}{\partial I}$ and $\pi_{B}\rightarrow\hat{\pi}_{B}=-i\frac{\partial}{\partial B}$. The corresponding Wheeler-DeWitt equation is:

[TABLE]

and its separable solution is given by DabQ :

[TABLE]

where

[TABLE]

where $r=2\sqrt{\frac{\tilde{V}_{0}}{A^{2}-9}}$ and

[TABLE]

Asymptotically for $I\rightarrow\infty$ the wave function $\Phi$ is an eigenfunction of the momentum operator $\hat{\pi}_{I}$ since $\hat{\pi}_{I}\beta(I)=\pi^{\infty}_{I}\beta(I).$ Thus, in the near curvature singularity regime $\Phi$ represents only those modes which are associated with the collapsing classical solution Birth ; DabKief . On the other hand asymptotically in the low curvature regime ($I\rightarrow-\infty$) the Bessel function $J_{in}$ can be written as $J_{in}(z)=\Psi_{1}+\Psi_{2}$ with $\Psi_{1}=\sqrt{\frac{1}{2\pi z}}e^{i\left(z-\frac{\pi}{2}in-\frac{\pi}{4}\right)}$ and $\Psi_{2}=\sqrt{\frac{1}{2\pi z}}e^{-i\left(z-\frac{\pi}{2}in-\frac{\pi}{4}\right)}$ where $n\equiv-\pi^{\infty}_{I}$. Since both $\Psi_{1}$ and $\Psi_{2}$ are eigenfunctions of the momentum operator $\hat{\pi}_{I}$ to the eigenvalues $-\pi^{\infty}_{I}e^{-I}$ and $\pi^{\infty}_{I}e^{-I}$ respectively we find that asymptotically for $I\rightarrow-\infty$ the function $\Psi_{2}$ represents the pre-big-bang collapsing branch while the function $\Psi_{1}$ represents the post-big-bang expanding branch. Thus the transition form the pre-big-bang low curvature initial state to the post-big-bang low curvature final state corresponds to the stationary quantum mechanical scattering of the plane wave (representing the particle with the definite value of the momentum) on the exponential potential barrier in the minisuperspace (see Fig. 2).

3 Emergence of multiverse in the third quantized varying $c$ and $G$ model

In this section we will argue that the third quantized non-minimally coupled varying constants model includes a scenario in which the whole bunch of universes emerges out of the initial vacuum (containing no universes) in the process analogous to the particle creation in an evolving background. The idea of third quantization assumes that the multiverse can be treated as a many-particle system (with non conserved number of particles) where the individual particles with their trajectories in spacetime are replaced by the universes making up the multiverse existing in the minisuperspace Strominger ; Robles . The Wheeler-DeWitt (14) can be rewritten in the form of Klein-Gordon equation:

[TABLE]

where $()^{\prime}=\frac{\partial}{\partial\eta}$ and $\Delta=\frac{\partial^{2}}{\partial x_{1}^{2}}+\frac{\partial^{2}}{\partial x_{2}^{2}}$, $m_{eff}^{2}=\overline{\Lambda}e^{-\frac{2}{r}\eta}$. The relation between the new variables $\eta$, $x_{1}$ and $x_{2}$ and the old variables $I$, $J$ and $B$ are:

[TABLE]

It is useful to define the following vectors: $\vec{x}\equiv(x_{1},x_{2})$, $\vec{k}\equiv(k_{1},k_{2})$. The formal analogy between Wheeler-DeWitt and the Klein-Gordon equations allows us to third quantize the wave function $\Phi$ by formally applying the standard quantization procedure of the quantum field theory. Following this procedure we promote the wave function $\Phi$ to be an operator $\hat{\Phi}$ acting on an associated quantum space of state mukhanov . The expansion of the field operator $\hat{\Phi}$ in terms of mode functions $v_{k}(\eta)$ is:

[TABLE]

where $k\equiv|\vec{k}|$. The mode functions $v_{k}(\eta)$ fulfil the mode equation

[TABLE]

where $\omega_{k}(\eta)=\sqrt{k^{2}+m_{eff}^{2}(\eta)}$, and the normalisation condition

[TABLE]

where $W$ is a Wronskian. We impose the standard commutation rules on the creation an annihilation operators $\hat{a}^{-}_{\vec{k}}$ and $\hat{a}^{+}_{\vec{k}}$:

[TABLE]

Instead of $v_{k}(\eta)$ any linear combinations of the form

[TABLE]

with the coefficients $\alpha_{k}$ and $\beta_{k}$ fulfilling the normalisation condition $|\alpha_{k}|^{2}-|\beta_{k}|^{2}=1$ could be used as the mode functions since any such linear combination is also a solution of the mode equation (23). The field operator $\hat{\Phi}$ expressed in terms of the mode functions $u_{k}(\eta)$ is given by:

[TABLE]

where $\hat{b}^{-}_{\vec{k}}$ and $\hat{b}^{+}_{\vec{k}}$ is a different set of the creation and annihilation operators satisfying the standard commutations relations given by:

[TABLE]

By comparing the two expansions (22) and (29) we can obtain the Bogolyubov transformations in the form:

[TABLE]

The explicit values of the Bogolyubov coefficients $\alpha_{k}$ and $\beta_{k}$ are:

[TABLE]

We define the two different vacuum states in a standard way:

[TABLE]

Following the paper Pimentel we interpret both vacuum states $|_{(a)}0\rangle$ and $|_{(b)}0\rangle$ as state vectors which represent the “ground states” of the multiverse, where the notion “ground state” refers to the “empty” multiverse (a multiverse containing no “a-universes” for $|_{(a)}0\rangle$ vector and a multiverse containing no “b-universes” for $|_{(b)}0\rangle$ vector). The creation and annihilation operators can be used to build the two sets of the excited state of the multiverse. Thus the quantum state of the multiverse containing “a-universes” with $m$ universes in mode $\vec{k}_{1}$, $n$ universes in in mode $\vec{k}_{2}$, etc. is represented by

[TABLE]

while the quantum state of the multiverse containing “b-universes” with $m$ universes in mode $\vec{k}_{1}$, $n$ universes in in mode $\vec{k}_{2}$, etc. is represented by

[TABLE]

Thus we have built an orthonormal basis on the Hilbert space of the multiverse. This means that any arbitrary vector representing state of the multiverse can be written as a linear combination of the excited states:

[TABLE]

The two vacuum states $|_{(a)}0\rangle$ and $|_{(b)}0\rangle$ are generally different vectors. This can be seen by calculating the expectation value of the “a-universe” number operator in the vacuum state $|_{(b)}0\rangle$:

[TABLE]

The divergent factor $\delta^{(3)}(0)$ accounts for an infinite spatial volume and hence the mean density of “a-universes” in the mode $\vec{k}$ is:

[TABLE]

Now, we are going to introduce a scenario in which the whole family of the universes is created from nothing. Let us first assume that initially the vector representing the quantum state of the multiverse is described by the “b-vacuum” state $|_{(b)}0\rangle$ which is completely specified by the set of mode functions

[TABLE]

that solve the mode equation (23) in the high-curvature limit (which appear for $\eta\rightarrow\infty$), where $A$ is the normalisation constant and $x\equiv\sqrt{\bar{\Lambda}}e^{-\eta/r}$. Since the structure of the vacuum state is controlled by the instantaneous value of the background curvature the initial vacuum state is identical with $|_{(b)}0\rangle$. Thus, initially the multiverse contains no universes at all. Due to stationarity of the scattering process the instantaneous quantum state of the multiverse does not evolve, so the final state of the multiverse (in the low-curvature limit that occurs for $\eta\rightarrow-\infty$) is still represented by $|_{(b)}0\rangle$. Due to the variation of the background curvature the notion of the vacuum is not invariant and the final vacuum state corresponding to the low-curvature limit (for $\eta\rightarrow-\infty$) is identical with $|_{(a)}0\rangle$ which, on the other hand, is completely specified by the set of mode functions

[TABLE]

being the solution of the mode equation (23) in the low-curvature limit (for $\eta\rightarrow-\infty$), where $B$ is the normalisation constant and $H^{(2)}_{-ikr}(x)$ is the Hankel function of the second kind. Therefore, the quantity:

[TABLE]

can be interpreted as the average number of the universes created from nothing described with the quantum number $\vec{k}$. By calculating the Wronskian in (46) we obtain that:

[TABLE]

The formula (47) describes the Bose-Einstein distribution for the temperature $T=\frac{1}{2\pi k_{B}}$, where $k_{B}$ is the Boltzmann constant provided that the energy of the bosons is identified with the quantity $kr$. Taking into account the eq. (17) we obtain that $kr=\pi^{\infty}_{I}\equiv 2\tilde{V}_{0}\sqrt{\frac{\Lambda}{A^{2}-9}}$. Thus, the distribution (47) expressed in terms of the cosmological constant $\Lambda$ reads:

[TABLE]

The formula (48) expresses the average number of the universes created from nothing with the value of the cosmological constant $\Lambda$ in the interval $(\Lambda,\Lambda+d\Lambda)$. From (48) we can see that the concentration of universes characterised by small value of the vacuum energy density is large. On the contrary, for larger values of the vacuum energy density the concentration of universes tends to zero.

4 Discussion

The canonically quantized non-minimally coupled varying $c$ and $G$ theory includes a scenario in which the transition from the pre-big-bang contraction to the post-big-bang expansion occurs as a consequence of the scattering of the plane wave on the exponential potential barrier in the minisuperspace. This is similar to the scenarios included in the string cosmologies Birth . The third quantization (a procedure analogous to the quantization of the Klein-Gordon field) of such a theory leads to the scenario in which the whole bunch of universes is created out of vacuum. The third quantization scheme was already applied to discuss the transition from expanding to contracting cosmological phases (and vice-versa) also in Buonanno ; Gasperini3 . The Hilbert space that emerges in the process of third quantization comprising all the states of the multiverse can be used as a base for introducing some standard notions of the ordinary quantum mechanics that exploit the linear structure of the space of quantum states. Particularly important here is the concept of quantum entanglement which relies on the notion of the tensorial product of the Hilbert spaces associated with quantum subsystems making up the whole physical setup. The scenarios in which the components of the multiverse experience their mutual presence via quantum entanglement was considered in Robles2 ; Robles3 . Concepts of that kind can constitute a basis for introducing models in which the quantum entanglement between different universes influences the cosmological observables making the idea of multiverse observationally testable Kraemer ; Kraemer2 .