Decaying universes and the emergence of Bell-type interuniversal entanglement in varying fundamental constants cosmological model

In this paper, we consider a high-curvature limit of the varying fundamental constants toy model in which both the value of the speed of light and the value of the gravitational constant are related to the values of the two non-minimally coupled scalar fields. The high-curvature limit motivates the application of the third quantization procedure to such a toy model which results in a theory that describes bosonic massive particles that move freely in the three-dimensional minisuperspace associated with the degrees of freedom of the original model. Motivated by the idea that in the quantum cosmological description the minisuperspace gets promoted to a real configurational space of the system we supplement the third quantized action of the considered model with an interaction term that allows for decay and scattering processes. We show that such interaction term induces a scenario in which a parent universe decays into two universes described by a maximally entangled Bell state. We eventually asses the strength of the entanglement, in the created pair of universes, by calculating the von Neumann entropy of entanglement.


Introduction
The idea of the multiverse is quite capacious. Most of the conceptualizations of such an idea can be assigned to one of Tegmark's levels of the multiverse [1]. Level one defines multiverse as regions beyond our cosmic horizon. The effects of quantum entanglement between different causally disconnected patches of space-time were investigated in [2,3]. Tegmark's level two defines multiverse as a collection of post-inflationary bubbles with possibly different values of the physical constants. Level three involves Everett's many-worlds interpretation of quantum physics while level four also includes all well defined mathematical structures. An interesting class of models, that can be classified as a level two multiverse, is based on the third quantization of the Wheeler-DeWitt wave function. In the process of the third quantization, the Wheeler-DeWitt wave function becomes an operator, which creates or annihilates universes characterized by some sets of quantum numbers (these are usually the momentum components in the minisuperspace) [4,5]. Such approach, which explicitly uses the third quantization formalism to describe the emergence of the entanglement in the multiverse, has inspired the area of research that explored different scenarios of generation of the interuniversal entanglement. It exploits the two different representations (the emergence of the entanglement resulting from the representation change is a generic feature of any quantum field theory [6]), namely the invariant one (for the first time found in [7]) which conserves the number of the universes and the diagonal one. Further development of the multiverse models relying on the third quantization includes postulating an interaction between the different universes which constitute the multiverse with each universe being represented by a different Wheeler-DeWitt wave function. It is usually assumed that the interaction between the universes can be represented by a term that formally describes an interaction between two harmonic oscillators coupled by a spring [8,9,10].
In this paper, we will construct, by acting in a similar spirit, a model, in which by including into the action an appropriate interaction term, the universes represented by the massive bosonic particles moving in the minisuperspace (or according to the third quantization concept by the massive Klein-Gordon fields) can decay and collide with each other. Such decay processes will result in the emergence of the entanglement in the produced pairs of the universes. We will also show that a model, with the abovementioned properties, can result from the nonminimally coupled biscalar gravity theory of varying fundamental constants in which both the speed of light and the gravitational constant are represented by the two nonminimally coupled scalar fields [11].

arXiv:2112.14166v3 [gr-qc] 26 Aug 2022
There have been many approaches to the idea of varying speed of light presented in the literature. Most of them, however, encounters profound conceptional problems such as violation of the Lorentz invariance [12,13,14,15,16,17]. The theory of varying fundamental constants assumed in this paper is largely based on the locally Lorentz-invariant varying speed of light (VSL) theory postulated by Magueijo in [14] that extends the VSL theories introduced in [12,13] which break the general covariance and, in consequence, require to choose a preferred reference frame (usually identified with the cosmological frame) to formulate the particular model. It was shown that such theories can solve the horizon, the flatness and the cosmological constant problem, however, the dynamics of varying speed of light is not given explicitly due to the lack of suitable terms in the action. The approach presented in [14] complements the mentioned VSL model by explicitly adding in the action the dynamical terms that govern the behaviour of the speed of light and the gravitational constant. It also proposes the definitions of covariance and local Lorenz invariance for the case of varying speed of light. It was also shown that such a model includes a scenario in which the whole multiverse emerges from nothing [18] and that the created pairs of the multiverses are entangled as confirmed by non-vanishing entanglement entropy [19]. The interuniversal entanglement was also investigated in the context of the third quantized varying fundamental constants cyclic cosmological models in [20] where it was argued that the third quantization naturally provides a thermodynamical description of the entanglement [21]. The emergence of the entanglement in pairs of created universes at the critical points of their evolution was investigated in [22].
Our paper is organized as follows. In Sec. 2 we introduce the non-minimally coupled varying speed of light c and varying gravitational constant G cosmological toy model and define its third quantized action. In Sec. 3 we add to the action an interaction term that enables scattering and decay processes. In Sec. 4 we show that pairs of universes, which are produced in the decay processes induced by the interaction term introduced in Sec. 3, are described by nearly maximally entangled Bell states. We also calculate the entropy of entanglement for such pairs. In Sec. 5 we give our conclusions.
2 Third quantized non-minimally coupled varying c and G cosmological model We start with a model of varying speed of light c and varying gravitational constant G introduced in [18,11] where both fundamental constants in the original Einstein-Hilbert action are replaced with a certain functions of the two scalar fields. Thus, the resulting action formally describes a non-minimally coupled scalar-tensor gravity theory with two scalar degrees of freedom. The considered model is largely based on the covariant and locally Lorentz-invariant varying speed of light theories postulated in [14] and is de-fined by the following action: where R is the Ricci scalar, Λ is the cosmological constant, ω is the parameter of the model and φ and ψ are the non-minimally coupled scalar fields, which values are, by definition, linked with the values of c and G via the following formulas: By introducing the new fields β and δ defined by: the action (1) can be recast into the Brans-Dicke type action which has the following form: An introduction of varying c into the action (1) or (6) breaks the general covariance of the original theory [14]. This, on the other hand, forces one to choose a preferred reference frame in which our theory is formulated. We will follow an approach proposed in [14] and formulate our model in the cosmological frame defined by the flat FLRW metric which reads: in which both the scale factor a and the lapse function N depend on the coordinate x 0 . The form of the action (6) in the cosmological frame defined by the metric (7) reads: where () ≡ ∂ ∂x 0 . Fixing the preferred reference frame completely requires to choose a specific form of the lapse function N . Since we are free in making such a choice, we will assume throughout the paper the following form of the lapse function N : The action (8) can be further simplified by the application of the following sequence of the field transformations: and The action (8) in the new variables η, x 1 and x 2 takes the following form: whereΛ =Ṽ Λ.
The hamiltonian corresponding to the action (12) is: where 2 are the respective conjugated momenta. Since both π x1 and π x2 are conserved quantities which is directly implied by (13) we can depict the classical evolution as a scattering of a point particle on the exponential potential barrier.
Since we are interested in the high-curvature near singularity behaviour, which occurs for η → ∞ (see Appendix A), we need now to switch to the canonical quantum cosmological framework which is governed by the Wheeler-DeWitt equation. An application of the Jordan quantization rules which boil down to replacing the canonical momenta with the operators in accordance with the following scheme: where and m 2 (η) =Λe − 2 r η . Eq. (14) is formally the same as the Klein-Gordon equation which allows us to perform the so-called third quantization of the Wheeler-DeWitt wave function Φ. The third quantization procedure is completely analogous to the procedure of quantization of the Klein-Gordon field and leads to the Fock space whose vectors represent the states of the considered model of the multiverse. The first step of the third quantization procedure requires writing the so-called third quantized action which leads to the Wheeler-DeWitt equation. The third quantized action for the case of the Wheeler-DeWitt equation given by (14) has the following form: where ∇ is a two-dimensional gradient operator associated with the variables x 1 and x 2 . The corresponding third quantized hamiltonian is: where the conjugated momentum π =Φ. The description of the classical evolution associated with the considered model can be find in the Appendix A.

Interacting universes in the minisuperspace
In our model, the set of universes is formally equivalent to the set of bosonic particles represented by the third quantized field Φ, characterized by specific values of the momentum p µ . Since, in the quantum cosmological description, it is the minisuperspace that constitutes the true configurational space of the system under consideration, it seems natural to include an interaction term that allows for the decay and the scattering processes. The simplest action that enables the abovementioned processes reads: where is the mass associated with the field Φ i . Since we are particularly interested in the description of the decay processes which are expected to occur in the high-curvature regime that takes place for sufficiently large value of the time parameter η (see Appendix A), we will be using a simplified version of the action (17) which reads: with m Φi = m i (η h−c ), where η h−c denotes the moment in which the high-curvature regime begins. Throughout the paper we will assume that m Φ1 > m Φ2 + m Φ3 . The corresponding hamiltonian is: where π i =Φ i constitute a set of conjugated momenta. The interaction term in (19) couples the one-particle states of the field Φ 1 with the definite energy and the definite momentum states of the two other fields, namely Φ 2 and Φ 3 . Such coupling can more precisely be described by the following formula:  where |1 Φ1 represents a zero-momentum one-particle state of the field Φ 1 while |Φ 2− − → p and |Φ 3 − → p represent the definite energy states (they are the eigenstates ofĤ f ree ) of the fields Φ 2 and Φ 3 , respectively, with opposite momenta of magnitude p (we assume the centre-of-mass reference frame). The presence of such interaction term allows for the processes in which a particle, represented by the field Φ 1 , decays into pairs of particles represented by the fields Φ 2 and Φ 3 . The states |Φ 2− − → p and |Φ 3 − → p are, on the other hand, coupled back to the state |1 Φ1 − → 0 (see Fig. 1). The second order of the perturbations theory leads to the intermediate three-particle states (see Fig. 2) which contribute to the shift in the vacuum energy. Since this contribution influences only the phase factor which multiplies the oneparticle state of the field Φ 1 it will be neglected in the following considerations [24,25].

Emergence of interuniversal entanglement via decay process. Bell states and the entropy of entanglement
It is convenient to use an interaction picture to describe the evolution of the state of the multiverse. Generally, the state of the multiverse |Ψ can always be represented by linear combination of the eigenstates of the free hamilto-nianĤ f ree : where c(η) are some complex coefficients, E n and |ϕ n are the eigenvalues and the eigenstates ofĤ f ree , respectively, and thus fulfil the following eigenequation: H f ree |ϕ n = E n |ϕ n .
By substituting Eq. (21) into the Schrödinger equation one gets: iċ l = n c n ϕ l |Ĥ I |ϕ n , whereĤ I ≡ e iĤ f ree ηĤ int e −iĤ f ree η is the interaction picture ofĤ int .
For the case of the coupling given by (20) we will be following the Wigner-Weisskopf approach [23] elegantly presented in [24]. In order to make the notation more compact we will shorten the formula (20) as follows: where |α constitutes the set of the eigenstates ofĤ f ree coupled to |G viaĤ I . We will also assume that the initial state of the considered setup is identical with the state represented by |G , which is equivalent to the assumption that c G (0) = 1 and c n =G (0) = 0. The Wigner-Weisskopf approach for the case of coupling abbreviated by (25) gives the following expressions for the coefficients c G and c α : where It can be shown [25] that: which means that the evolution is unitary even in the case of perturbative approach. We will be considering the following scenario: as the evolution begins, the state of the universe is characterized by a high value of the curvature, since the value of the scale factor a is close to zero. This, on the other hand, corresponds to an infinite value of the time parameter, i.e. η → ∞ (see Appendix A). In the state of high curvature, the interaction termĤ int switches on, what enables the decay processes. We will argue that in such processes, the one-particle zero momentum state of mass m Φ1 decays into a pair of particles of masses m Φ2 and m Φ3 , maximally entangled in the momentum space.
The formula (28) adjusted to our particular setup is: where p ≡ − → p . We will also define the following function: which derivative is: By calculating the integral in Eq. (26) by parts, one gets: The first term, in the formula above, is of second order inĤ I while the second term, is of fourth order inĤ I . This means that up to the leading order the equation that governs the time evolution of c Φ1 is: The integral in Eq. (31) which explicit form is given by: in the limit for η → ∞ is: where P denotes a principal value while E 2 Φ2 (p) = p 2 + m 2

Φ2
and E 2 Φ3 (p) = p 2 + m 2 Φ3 . The solution of (34) in the limit η → ∞ gives: In order to arrive to a more explicit version of the expression for Γ we need first to calculate After decomposing all the three fields Φ 1 , Φ 2 and Φ 3 into modes with the following expression: where A denotes the volume of compact rectangular spacelike part of the minisuperspace, the M (p) assumes the following form: In the three dimensional minisuperspace the explicit form of the expression for the Γ is: The state of the three fields Φ 1 , Φ 2 and Φ 3 can be expressed as: The formula (27) gives the following expression for the coefficients C Φ2Φ3 (p, η) in (44): where m Φ1,R = m Φ1 + ∆E Φ1 is the renormalized mass (in the following we skip the index R in m Φ1,R and assume that m Φ1 represents the renormalized mass m Φ1,R ). The expression (44) explicitly describes the entanglement of the momentum eigenstates.
Since we are interested in calculating the entropy of entanglement we need first to calculate the reduced density matrix associated with one of the universes. After tracing away the degrees of freedom of the universe represented by the field Φ 2 we obtain the following reduced density matrix for the universe represented by the field Φ 3 : where ρ ≡ |Ψ (η) Ψ (η)|. The von Neumann entropy is then given by: In the narrow width limit Γ << m Φ1 , m Φ2 + m Φ3 the function |C Φ2Φ3 (p, η)| 2 is sharply picked at so the last term in Eq. (47) can be approximated by the following expression: In the narrow width limit we have that what ensures fulfillment of the unitary evolution condition (29). Since the average number of the universes with momentum of magnitude p represented by the field Φ 2 (or Φ 3 ) is given by: the formula (50) also gives the total number of the universes of either type (Φ 2 or Φ 3 ) produced in the volume A. Thus, the decay results in a production of only one universe of either type in region A.
Taking into account the formula (50) the von Neumann entropy (47) can be expressed as: which in the high-curvature limit (η → ∞, see Appendix A) gives: The identification E = i ∂ ∂η and the fact that asymptotically, for η → ∞, the value of π η =
Since in the considered case the function |C Φ2Φ3 (p, η)| 2 is narrowly peaked at p = p d a typical term that contributes to the expansion (44) has the form: which means that the pair of the universes produced in the decay process is described by the nearly maximally entangled Bell state.

Conclusions
A trend in which the minisuperspace is treated as a fundamental scene, where the physical phenomena occur, gained a noticeable representation in the literature [27,28,29,30,31,32,33,34,35]. Such an approach motivates an inclusion of different types of interactions between the Wheeler-DeWitt wave functions (representing individual universes of the multiverse) such as the interaction represented by the harmonic potential [8,9,10]. Moreover, postulating such interactions and investigating the consequences of their existence seems to be the only way to make the concept of the multiverse a physical theory which is, at least in principle, amenable to observational verification. Following this trend, we include in the action, stemming from the varying fundamental constants model, an interaction term that enables scattering and decay of the universes represented by the third quantized Wheeler-DeWitt wave function. The product of the decaying parent universe is a nearly maximally entangled pair of universes approximately described by a Bell state. Thus, the emergence of a nearly maximally entangled pairs of universes results from the inclusion of a very natural type of the interaction that allows for scattering and decay processes. Such type of entanglement is typical for the interaction term introduced in our multiverse toy model [24]. However, the form of the state that describes the products of any decay process generically depends on the particular character of the interaction responsible for that decay [26]. The entanglement concerns the spatial components of the minisuperspace momentum associated with individual universes. The strength of the entanglement measured by the von Neumann entropy depends on the values of the cosmological constants in each of the universes in the created pair and achieves higher values when both cosmological constants have similar values. Interestingly, the presented approach relies on the standard interpretation of the Fock space, which means that it treats the representation dependent orthonormal vectors, that form the basis in the Hilbert space of the multiverse, as vectors representing occupation with universes in a state completely determined by a particular set of proper quantum numbers. The standard interpretation of the Fock space was also assumed in [18] where it was used to derive a scenario in which the whole multiverse subjected to the Bose-Einstein distribution emerged from nothing. On the other hand, there have been developed approaches to the problem of interuniversal entanglement, in which the standard quantum field theoretical interpretation of the Fock space was overridden in the sense that the Hilbert space basis vectors are assumed to define the excited states of the universes with a specific value of the momentum in the minisuperspace [19,20]. Thus, the mechanism of interuniversal entanglement generation presented in this paper, seems to be quite a natural one since it does not assume any non-standard interpretation of the basic concepts of quantum field theory.

Appendices
A Classical cosmological time evolution in the non-minimally coupled varying c and G model The classical evolution associated with the model defined in Sec. 2 is given by the following formulas [11]: where M = 3−A 2 9−A 2 , W = 2A 2 9−A 2 while D and F are some integration constants. The variables x 0 is connected with the cosmological timex 0 via the following expressions: where we assumed that A 2 > 9. The set of the solutions above defines the time evolution of the scale factor a, speed of light c and the gravitational constant G. The qualitative behaviour of the three parameters is depicted in Fig. (A1). It can be seen that the model contains the pre-big-bang contraction (forx 0 < 0) and the post-big-bang expansion (forx 0 > 0) with both phases separated by the curvature singularity atx 0 = 0. The Hamilton's equation of motion associated with the hamiltonian (13) gives the following time evolution: where E and P are some integration constants. The solution (A.4) allows us to define the high-curvature limit for η → ∞ which corresponds to near curvature singularity evolution and the low-curvature limit for η → −∞ which corresponds to the late evolution (far form the curvature singularity). The asymptotic value of the momentum π η in the high-curvature limit (for η → ∞) is: