Spinor wave function of the Universe in non-minimally coupled varying constants cosmologies

In this paper, we introduce a non-minimally coupled varying speed of light and varying gravitational constant cosmological toy model. Using the Eisenhart–Duval lifting method, we extend the original minisuperspace of the model and depict the evolution of the system in the presence of the potential term as a geometrical flow associated with the lifted metric. We write the Dirac–Wheeler–DeWitt equation, which solution is a spinor wave function of the Universe. Then we find the solution of the Dirac–Wheeler–DeWitt equation, which describes the emergence of two early universe–antiuniverse pairs that differ with the conserved quantity, which is an analog of the spin.


Introduction
In [1,2], it was shown that the extension of minisuperspace using the Eisenhart-Duval lift [3,4,5] provides a natural framework for introducing a Dirac-type equation. This equation can replace the Wheeler-DeWitt equation, which, as a hyperbolic partial differential equation, does not necessarily provide a positively defined probability density. This issue was addressed by introducing a Diracsquare root of the Wheeler-DeWitt equation [6,7,8,9] which leads to ambiguities related to factor ordering, or by using supersymmetric quantum mechanics [10,11,12,13,14,15]. The crucial concept in the Eisenhart-Duval lift is the notion of an extended minisuperspace, which is obtained by adding an auxiliary dynamical variable to the original set of variables defining the initial minisuperspace associated with the considered dynamical system. This auxiliary variable parametrizes an additional dimension of the extended minisuperspace. In such an extended minisuperspace, the dynamics of the system, which generally is not a geodesic one in the presence of the potential term, can be transformed into a geodesic dynamics associated with the socalled lifted metric that is suitably defined on the extended minisuperspace, since the "extended" equation of motion is expressed by the Laplace-Beltrami operator and thus exhibits covariance. Moreover, in [1,2] it was argued that the requirement of the covariance in the extended minisuperspace resolves the problem of factor ordering and may provide a consistent way to formulate the Dirac square root.
Most attempts to incorporate varying speed of light into cosmological models result in a violation of Lorentz invariance [16,17,18,19,20,21]. Such theories also violate general covariance, which requires a preferred reference frame, usually identified as the cosmological frame. Nevertheless, some solutions have been proposed that address the horizon, flatness, and cosmological constant problem [16,17]. However, these solutions come at the cost of postulating ad hoc the time dependence of the speed of light, as the governing terms of the speed of light dynamics are lacking in the action.
The varying speed of light and varying gravitational constant cosmological model considered in this paper is based on the model proposed in [18] that includes covariance and the local Lorentz invariance in the case which allows for variation of the speed of light. In such a model, both the speed of light and the gravitational constant are additional dynamic degrees of freedom represented by the two scalar fields with appropriate kinetic terms in the action. It was shown that such model includes a scenario in which the Universe is created in the process which is an analog of quantum scattering on the potential barrier [22] (this is similar to what was found in [23]), that the third quantization of the model leads to a scenario of the emergence of the multiverse from nothing [24], and that the created pairs of universes are entangled which is reflected by non-vanishing entanglement entropy [25,26]. Third quantization of slightly different models, in which the fundamental constants undergo cyclic evolution [27], yields insight into the relationship between quantum entanglement and thermodynamics [28].
As discussed earlier, the solutions of Wheeler-DeWitt equation which describes the canonically quantized model of varying speed of light and gravitational constant, include interesting cosmogenesis scenarios that are similar to those found in string cosmological models [23,29,30]. In this paper, we will investigate the potential implications of incorporating spinorial properties of the wave function as a solution to the Dirac-like equation in the quantized varying fundamental constants cosmological model, as such an approach may lead to some novel features of quantum cosmogenesis.
Our paper is organized as follows. In Sec. 2 we describe the non-minimally coupled varying speed of light c and varying gravitational constant G cosmological model and provide the description of its high-curvature limit. In Sec. 3 we follow the Eisenhart-Duval lifting scheme and formulate the dynamics of our model on the extended minisuperspasce. Then, utilizing the geometric nature of the dynamics of our "lifted" system, we formulate the Dirac-Wheeler-DeWitt equation. In Sec. 4 we obtain the solution to the Dirac-Wheeler-DeWitt equation and establish its connection to the scenario of cosmogenesis. In Sec. 5 we give our conclusions.
2 The non-minimally coupled varying c and G cosmological model Our model, which proposes variations in the speed of light c and the gravitational constant G, is primarily represented by the theory of non-minimally coupled bi-scalar gravity. In this theory, the values of both fundamental constants are connected to the values of the two scalar degrees of freedom [22]. The model under consideration is based on the covariant and locally Lorentz-invariant theory of varying speed of light, as described in [18]. It is defined by the following action: where φ and ψ represent non-minimally coupled scalar fields, R denotes the Ricci scalar, Λ acts as the cosmological constant, and ω is a parameter of the model. The action specified in equation (1) was derived from the original Einstein-Hilbert action by substituting the functions of the scalar fields φ and ψ for the speed of light c and the gravitational constant G. The exact connection between the scalar fields φ and ψ and the constants c and G is given by the following formulas: The application of the following fields transformations: causes the original action (1) to take the form of the Brans-Dicke action given by the following expresion: The dependence of the speed of light on the space and time variables causes a violation of the general covariance in our model, which is a common characteristic of most theories that propose a variable speed of light. To address this, a specific frame of reference, known as the "light frame", must be selected. This frame serves as the preferred reference frame for the formulation of the model. Based on the proposal, outlined in [16,17], we identify the light frame for our model with the cosmological frame defined by the FLRW metric: where N is the lapse function and a is the scale factor, both dependent on the coordinate x 0 . The insertion of the metric from equation (7) into equation (6) gives the following form of the action of our model in the cosmological frame: where () ′ ≡ ∂ ∂x 0 . The solution to the model described by action (8) expressed in the gauge defined by N = a 3 δ is [22]: while D and F are the integration constants. The relation between x 0 and the rescaled cosmic timex 0 is given by [22]: where as in [22] we restrict the considered class of models to the cases with A 2 > 9. The pre-big-bang contraction occurring forx 0 < 0, followed by the post-big-bang expansion taking place forx 0 > 0 are both described by the solution given by (9) and (10), along with (11). The two phases are separated by the curvature singularity at x 0 = 0. The evolution of the fundamental constants c and G is also conveyed in formulas (9) and (10). According to these formulas, the gravitational constant G tends to zero while the speed of light c becomes infinite as the universe approaches the curvature singularity located atx 0 = 0 (see Fig. (1)). By utilizing the field transformations given by: whereṼ 0 = 3V0 8π and r = 2 Ṽ 0 A 2 −9 the action (8) can be transformed to the following form: whereΛ =Ṽ 0 Λ. The Hamiltonian associated with the action (13) reads: where 2 are the conjugated momenta. From the structure of the Hamiltonian (14), it can be inferred that both π x1 and π x2 remain conserved throughout the evolution. Consequently, the classical evolution can be seen as a particle scattering off the exponential potential barrier. The solutions to the Hamilton's equations associated with the Hamiltonian (14) are as follows: where E and P are integration constants. Upon inspection of the solution (15), it is evident that η can define two distinct regimes. The first one is the high-curvature regime, which is characterized by the vanishing scale factor a → 0 and corresponds to η → ∞. The second one is the low-curvature regime, which is defined by higher values of the scale factor a and occurs for η → −∞. Furthermore, it can be verified that the high-curvature regime, specifically when η → ∞, is characterized by the following asymptotic values of momentum π η : As for the low-curvature regime (as η → −∞), we can observe that: (expanding post-big-bang solution).
To derive the Wheeler-DeWitt equation that characterizes the quantum mechanical aspect of the model, we implement the Jordan quantization rules which involve replacing the canonical momenta with the corresponding operators according to the following scheme: The resulting Wheeler-DeWitt equation can be expressed as: 3 Quantization of the non-minimally coupled varying c and G cosmological model. The Dirac-Wheeler-DeWitt equation We will be following the Eisenhart-Duval lifting scheme [3,4] adjusted to the cosmological setup based on the scalartensor gravity [1,2,5]. In this method the minisuperspace associated with considered model is extended by introducing additional dynamical variable which makes possible to describe the evolution of the system in purely geometrical terms even in the presence of the potential. In other words the nongeometrical evolution can be depicted as a geodesic flow in the extended minisuperspace, that respects the lifted metric (the metric on the extended minisuperspace) since the Hamiltonian of the lifted system is represented by the Laplacian. The starting point is the Lagrangian defined by the action (13) which reads as: The lifted Lagrangian associated with (19) is the following [1]: where x 3 parametrizes an additional degree of freedom of the extended minisuperpace. The metric of the extended minisuperspace (the lifted metric) is then given by [1]: Since, the Hamiltonian constraint whereP α =G αβ x ′β , is conformally invariant, namely, it is essentially the same asG αβ → Ω 2Gαβ , we will be using the following conformally equivalent extended minisuperspace metric with Ω 2 = 2Λe −2 η r 1 n−2 , where n stands for the dimension of the extended minisuperspace. Since in our case n = 4 the conformal factor is given by: Thus, the conformally equivalent extended minisuperspace metric is the following: The covariant Wheeler-DeWitt equation for the extended case is given by: where G is the determinant of G αβ . Substitution of (25) into (26) leads to the following Wheeler-DeWitt equation for the extended minisuperspace: If we impose the condition the Wheeler-DeWitt equation for the extended minisuperspace (27) reduces to the initial one (18). As in the case of Wheeler-DeWitt equation we will use the covariance in the extended minisuperspace as guiding principle which leads to the proper form of the Dirac-Wheeler-DeWitt equation [1,2]. Since the Dirac equation without the mass term has the conformal covariance [31] we will again use (25) as the lifted metric. Thus, the covariant Dirac-Wheeler-DeWitt equation is: where γ A are Dirac matrices: with k = 1, 2, 3, σ k being Pauli matrices and I being an identity matrix. The coefficients e α A in (29) are four vector fields defined by the expression η AB = e α A e β B G αβ and D α is a covariant derivative given by: where with the spin connection ω ABα defined as: and The explicit form of the vector field e α A is given by: and the non-vanishing elements of the spin connection ω ABα are: √Λ .
The explicit form of the Dirac-Wheeler-DeWitt equation in the extended minisuperspace is then given by: The reduction of the extended Dirac-Wheeler-DeWitt equation (37) to the initial minisuperspace can be obtained by the application of the assumption which reads: 4 Spinor wave function of the Universe in the non-minimally coupled varying c and G model The solution of the Dirac-Wheeler-DeWitt equation (37) reduced with condition (38) to the initial minisuperspace is: with φ 1 =e i k· x e −iα α ikr− 1 4 (C 1 U (ikr + 1, 2ikr + 1, 2iα) where k = (k 1 , k 2 ), x = (x 1 , x 2 ), k = k 2 1 + k 2 2 , U and L are the confluent hypergeometric function of the second kind and the associated Laguerre polynomial, respectively, α = r √Λ e − η r and {C 1 , C 2 , C 3 , C 4 , C 5 , C 6 , C 7 , C 8 } denotes a set of integration constants.
We introduce a new representation in which the solution (39) takes simpler form. The unitary matrix that defines such a representation is: and the transformed wave function (39) is given by: We notice that the unitary transformation given by (44) leaves the four momentum operator unchanged. It also preserves the form of the x 3 -axis projection spin operator Σ x3 given by: Let us define now the four spinor wave functions each separately solving the Dirac-Wheeler-DeWitt equation (37): where A, B, C and D are some constants. The spinor wave functions (47), (49) in the low-curvature regime for η → −∞ (α → ∞) represent the positive and the negative frequency modes, respectively, with 1/2 spin along x 3 axis, while (48), (50) represent the positive and the negative frequency modes, respectively, with −1/2 spin along x 3 axis since in this regime: where π η(−∞) = √Λ e − η r . The general solution of (37) that holds for any value of the curvature which reads as: represents purely quantum state in the high-curvature limit for η → ∞ since in such limit (53) it is not peaked over any classical trajectory. This can be seen by inspecting (53) which in the limit η → ∞ has the following form: and d = Γ (ikr) Γ (−ikr)Γ (2ikr + 1) .
The formula (53) shows that in the low-curvature limit for η → −∞ from the pure quantum regime there emerge classical trajectories that represent the expanding post-bigbang universe-antiuniverse pair with the spin projection 1 2 (Ψ −,+ 1 2 for the universe and Ψ +,+ 1 2 for the antiuniverse) and the expanding post-big-bang universe-antiuniverse pair with the spin projection − 1 2 (Ψ −,− 1 2 for the universe and Ψ +,− 1 2 for the antiuniverse). The scenario discussed here differs from the scenario of cosmogenesis, which involves creation of the Universe via quantum scattering on the exponential potential barrier in the minisuperspace, that appears in the non-minimally coupled varying constants cosmological model, where the ordinary Wheeler-DeWitt equation provides a solution representing the quantum state of the Universe [22].

Conclusions
We have shown that the Eisenhart-Duval lift which is based on the extension of the initial minisuperspace of the model by adding an auxiliary dimension may efficiently be applied in the context of non-minimally coupled varying speed of light and varying gravitational constant model which results in transforming the evolution of the system in the presence of time-dependent mass term into a geodesic evolution compatible with the lifted metric. It should be stressed that the covariance in the extended minisuperspace constitutes a guiding principle that leads to the proper equation of motion [1,2]. In particular we have used this principle in order to construct in a consistent way the Dirac-Wheeler-DeWitt equation which incorporates the spinorial characteristics into the wave function of the Universe. This new quality seems to be particularly interesting since it may substantially influence the distribution when considering an ensemble [2]. We have also shown that there exist the solutions to the Dirac-Wheeler-DeWitt equation that asymptotically, in the low-curvature regime, can be interpreted as the positive and negative frequency modes with the conserved quantity, which is an analog of the spin projection. It should be stressed, however, that asymptotically these solutions exhibit decoupling of the spin analog from the Universe trajectory in the extended minisuperspace which is expected since the analogous situation occurs in case of Dirac equation description of the free electron.
The considered model is defined by the action which is similar to the low-energy effective action of the string cosmological models [23,29,30], therefore, besides the ordinary post-big-bang expanding phase of the evolution it also provides the pre-big-bang cosmological scenario. We have found that the spinor wave function of the Universe which is a solution to the Dirac-Wheeler-DeWitt equation in the high-curvature (near big-bang) regime describes the highly unclassical behavior since it is not peaked over any classical trajectory. On the other hand, in the lowcurvature regime (far from the big-bang) the wave function is peaked over classical paths which can be interpreted as a transition from pure quantum behavior to the quasiclassical one and emergence of the two post-big-bang expanding universe-antiuniverse pairs with the opposite spin orientations.
Within the concept of the multiverse, there is a growing approach that regards the minisuperspace as the fundamental arena for physical phenomena to take place [23,26,29,30,32,33,34,35,36,37,38,39]. Particular models include different types of interactions between the wave functions that represent individual universes [26,32,35,36,37,38,39]. This may suggest that a comprehensive inclusion of spin in quantum cosmology requires an application of third quantization procedure [2], in which the third-quantized action is invariant with respect to some kind of local symmetry.

Appendices
A The Eisenhart-Duval lift, the extended Wheeler-DeWitt and the Dirac-like equations for a cosmological system

A.1 The essentials of the Eisenhart-Duval lift
The Eisenhart-Duval lift [3,4,5] is a formalism that can be used to represent physical systems subjected to conservative forces as an equivalent free system moving on a higher-dimensional curved manifold. Let us start with the Lagrangian for a system of n degrees of freedom (fields): where ϕ i (t) represent a homogeneous scalar field, k ij is the metric tensor in the configuration space, V (ϕ 1 , ..., ϕ n ) represents the potential while the overdot denotes differentiation with respect to t. Variation of (A.1) with respect to ϕ i gives the following equations of motion: where the symbol , i denotes differentiation with respect to ϕ i , k ij satisfies k il k lj = δ i j and We can see that (A.2) formally represents a geodesic equation of a particle that is subjected to a force. By utilizing the Eisenhart-Duval lift formalism, we can expand upon this geometric interpretation by generating a manifold of higher dimension for the field space. This manifold will be characterized by trajectories that follow geodesics as described by equation (A.2) but without any external force. For this purpose we introduce an auxiliary field χ and the so called "lifted" Lagrangian which reads: where M denotes an arbitrary constant. Such a "lifted" Lagrangian can be written as follows: We notice that the Lagrangian (A.5) is identical to the Lagrangian of a free particle moving in a curved spacetime equipped with metric G AB which means that the trajectories of the system are identical to the geodesics of the extended manifold. It can be shown that the geodesics of the extended manifold correspond to the trajectories of the original system described by equations (A.2). A detailed proof of this correspondence can be found in [5].

A.2 The extended Wheeler-DeWitt equation for a homogeneous and isotropic Universe
The Eisenhart-Duval method can be extended to the quantum level described by the Wheeler-DeWitt equation, at least for a homogeneous and isotropic model of the Universe with a single scalar field defined by the following action [1]: where V (ϕ) is a potential term. For tha FLRW metric given by: where dΩ 2 3 represents a maximally symmetric three-space with a constant Ricci curvature (3) R ij = 2Kg ij , where K is a constant, and N is the lapse function, the corresponding Lagrangian is: where U (a, ϕ) = a 3 V (ϕ) − 1 2 Ka (while deriving the formula above it was assumed that κ 2 = 6).
The Hamiltonian associated with (A.9) is: where Π a = − aȧ N and Π ϕ = a 3φ N are the conjugate momenta. By substituting appropriate operators for the conjugate momenta in accordance with the Jordan rules Π a → −i ∂ ∂a and Π ϕ → −i ∂ ∂ϕ one finds the standard Wheeler-DeWitt equation for the cosmological system:

A.3 The extended minisuperspace Dirac-like equation
Previous studies that aimed to take the square root of the Wheeler-DeWitt equation exhibited arbitrariness in handling the potential term [6,7,8,9]. The idea here is to apply the covariance principle in the extended minisuperspace, similarly as in the case of the extended Wheeler-DeWitt equation, in order to establish a consistent approach for dealing with the potential term, and thus, provide a systematic approach for determining the form of the extended minisuperspace Dirac-like equation. The Dirac-like equation in the extended minisuperspace can be formulated as [1]:γ α D α Ψ ≡ γ A e α A D α Ψ = 0, (A. 17) where the vital role is played here by the extended minisuperspace metric G MN in a particularly chosen gauge (see (23) or (A.15)). It is worth emphasizing that the solutions provided by (A.16) and (A.17) are physically inequivalent solutions [1,2].