Time-covariant Schr\"{o}dinger equation and invariant decay probability: The $\Lambda$-Kantowski-Sachs universe

The system under study is the $\Lambda$-Kantowski-Sachs universe. Its canonical quantization is provided based on a recently developed method: the singular minisuperspace Lagrangian describing the system, is reduced to a regular (by inserting into the dynamical equations the lapse dictated by the quadratic constraint) possessing an explicit (though arbitrary) time dependence; thus a time-covariant Schr\"{o}dinger equation arises. Additionally, an invariant (under transformations $t=f(\tilde{t})$) decay probability is defined and thus ``observers'' which correspond to different gauge choices obtain, by default, the same results. The time of decay for a Gaussian wave packet localized around the point $a=0$ (where $a$ the radial scale factor) is calculated to be of the order $\sim 10^{-42}-10^{-41}\mathrm{s}$. The acquired value is near the end of the Planck era (when comparing to a FLRW universe), during which the quantum effects are most prominent. Some of the results are compared to those obtained by following the well known canonical quantization of cosmological systems, i.e. the solutions of the Wheeler-DeWitt equation.


Introduction
Through the years, numerous attempts have been made to quantize gravity. Some of them more conservative, [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16] while others more radical [17][18][19][20][21][22][23][24][25]. More or less, all of them share a number of basic reasons to pursue such a quest. The more fundamental are: the expectation that at some energy scale, all the fundamental interactions should be unified. Thus, there must exist a quantum theory of gravity in the picture of quantum interactions. The second motivation is based on the existence of abnormalities in the classical General Relativity and/or its modifications. By abnormalities we refer to the singularities of black holes or the Big Bang itself for instance. The third and final is the incompatibility of General Relativity with Quantum Mechanics (and even more Quantum Field Theory) regarding the notions of space and time. In the former, spacetime is a dynamical entity that affects and, at the same time, is affected by any form of energy, while in the latter, it is an unaltered "frozen" and "external" arena. For some interesting reviews on these matters we recommend [21,[26][27][28][29][30]. Also, for an interesting new approach regarding the notion of time and resolution of singularities one may look at [31][32][33] and references therein.
To tackle the problem by considering the full theory of Gravity is a noble, but quite difficult task; thus, over the years, researchers have chosen to study simplified cases. The simplest appear to be the minisuperspace models, which are obtained via symmetry reduction of the full theory. By assuming a certain level of symmetry applied on the spacetime metric and the matter fields, the system's degrees of freedom are reduced to finite from infinite. Thus, instead of field theory, the tools of point-like mechanics can be used. Such kind of spacetimes are the FLRW geometry, the Schwarzschild and Reissner-Nordström black holes and the anisotropic Bianchi and Kantowski-Sachs types [34][35][36]. There are many works dealing with such kind of simplified cases, with very interesting results. It is quite difficult to cite them all, so we just mention a handful of them: The quantization of the FLRW geometry in the presence of a scalar field has been studied extensively from various perspectives [37][38][39][40][41][42][43][44]. An interpretation via the de Broglie-Bohm theory of the solutions to the Wheeler-DeWitt equation for the Reissner-Nordström-de Sitter black hole was presented in [45]. The authors of [46][47][48] exploit the notion of conditional symmetries to provide a quantum description of some Bianchi types and of the Reissner-Nordström black hole. There were also attempts in providing a generalized definition of probability [49] and also a quantum description through the definition of a time related to the homothetic symmetry of the minisuperspace metric [50]. Quantization in a larger space through an Eisenhart lift was recently employed in [51]. An analogy between the Schwarzschild-Reissner-Nordström black holes and hydrogen atoms were drawn in [52,53], where their energy spectrum was obtained based on some constructed Schrödinger type of equation. A study related to a mass operator and the existence of non-zero mass uncertainty can be found in [54,55]. The quantization of an inhomogeneous cosmological model has been performed in [56][57][58][59].
In this work we are dealing with a Kantowski-Sachs spacetime. Various aspects of the canonical quantization of Kantowski-Sachs models have been previously discussed in the literature: Louko and Vachaspati [60] studied whether or not the Vilenkin boundary condition could provide a unique wave function for the Kantowski-Sachs minisuperspace model. In [61] the authors provide a method for approximate evaluation of the path integral for spatially homogeneous minisuperspace models and apply it to the Kantowski-Sachs geometry among others. Finally, the complete solutions to the Wheeler-DeWitt equation were provided in [62] for generalized Kantowski -Sachs models with cosmological constant and pressureless dust. Several works also exist in the quantization of cosmological models in extended theories of gravity [63][64][65][66][67][68].
As any gauge theory [69][70][71] contains spurious degrees of freedom, so does gravity. This results to its description via a singular Lagrangian, which is equivalent to the existence of constraint equations. The question then is how we deal with the constraints during the quantization procedure. There are two main approaches in this problem, which are nicely described in [29]: The first approach is called the reduced phase space quantization. This is based on trying to find some canonical (invertible) transformation from the initial canonical variables to a complete set of gauge invariant combinations, plus a set of specific functions of the spacetime coordinates. The latter implies that the gauge freedom is completely broken and the constraint equations are strongly solved (since the reduced Hamiltonian is expressed solely in terms of the dynamical degrees of freedom) with the introduction of the Dirac Brackets. Thus, on the embedded hypersurface defined by the gauge choices, a Schrödinger equation can be constructed. This can be extended by considering the true degrees of freedom as functions of the spatial coordinates on the hypersurfaces, leading eventually to the so called "multi-finger" Schrödinger equation [29,72,73]. Note that this generalization is possible only in the case of field theory and not for the minisuperspace models. In many cases it is quite difficult to obtain the reduced Lagrangian, thus many authors have turned to the definition of a time standard based on matter fields [74][75][76][77][78]. To the advantages of this procedure are included: the capability of defining a Hilbert space and thus provide a probabilistic interpretation, as well as the dependence of the Hamiltonian only on the true degrees of freedom. The main disadvantage is the existence of many different quantum descriptions for the different gauge choices, which are not equivalent by default [79,80]. This is related to the fact that there is no unique canonical transformation, which raises problems at the quantum level, due to the non-commutative nature of the order in which you perform the transformation and the quantization. The second approach is covariant, thus by definition it holds for every gauge choice. Being covariant implies that the constraints are carried to the quantum level and imposed to the physical states as restrictions, resulting eventually to the so called Wheeler-DeWitt equation [1,2,[11][12][13][14][15][16]. The main problems are the absence of a well defined Hilbert space and, since time is not an external parameter, the inability to obtain evolution of states in the usual quantum mechanical sense. For a discussion related to the differences between covariant and reduced phase space quantization see [81,82].
The present work constitutes an extension of a relatively recent published paper [83] regarding the Reissner-Nordström black hole. The main idea was to somehow render the constraint equation trivial, in the sense that it will be satisfied modulo the dynamical equations, while at the same time maintain intact the gauge invariance. The result was that the singular Lagrangian describing the Reissner-Nordström black hole reduced to a regular one with an explicit "time" dependence (note that "time" for the Reissner-Nordström case is actually the radial distance, hence the quotations). The parameter of "time" was identified by considering the remaining gauge degrees of freedom as mere functions of it and not as dynamical quantities. Based on this method, a "time"-covariant Schrödinger equation was constructed, giving a somewhat intermediate picture in comparison to the two previous approaches. The covariance, in what regards purely "time" reparametrizations, is maintained, while a Schrödinger equation incorporating each different gauge choice of this type is produced. Here, we are interested to apply this method to a cosmological minisuperspace model and specifically to the geometry of Λ-Kantowski-Sachs universe. Furthermore, we extend it a bit further by defining an invariant decay probability. In the end, we also present the de facto covariant method of the typical Wheeler-DeWitt quantization of the given model, in order to study where the two approaches meet and what differences there exist. At this point, we need to mention of another interesting work, where a similar procedure with the one we follow here is employed. It regards the derivation of a Schrödinger equation and the introduction of time dependence in the wave function through the process of gauge fixing the scale factor while leaving the lapse as a degree of freedom, for more details see [84].
The paper is organized as follows: In section 2 the classical description of the system is provided by deriving the solutions and the relative reduced Lagrangians and Hamiltonians. The canonical quantization procedure is presented in section 3 alongside with what we call the time-covariant Schrödinger equation and its solutions. Next, in section 4, the invariant probability decay rate is calculated for the case of a Gaussian initial state. Additionally, an estimate on the inflation epoch based on the quantum effects is presented. In section 5, we briefly revisit the typical Wheeler-DeWitt quantization for a comparison to the previous approach and in the last section we gather our conclusions.

Solutions
Our starting point is the Einstein-Hilbert action in the presence of a cosmological constant Λ, where g = det(gµν) is the determinant of the space-time metric and R is the Ricci scalar. For simplicity, in what follows we adopt the units κ = c = 1.
We consider the following line element, in spherical coordinates x µ = (t, r, θ, φ), where the degrees of freedom (n, a, b) are functions of the coordinate t. With this ansatz, Einstein's equations in which Rµν represents the Ricci tensor, reduce to a set of ordinary differential equations for the functions of t entering the metric. The field equations are equivalent to: For our purposes, by following the same steps as in [83], the constraint equation (4) should be solved, if possible, with respect to the lapse n. It turns out that this proposal was not unknown and comes by the acronym BSW from the initials of Baierlein, Sharp and Wheeler [72]. A novel derivation of General Relativity, partially based on the above method can be found in [85]. In our case, due to the presence of the cosmological constant, this is possible only for b(t) = ± 1 √ Λ . From now on, without loss of generality, we consider b(t) > 0 and consider the cases b(t) = 1 Before we proceed with the study of these branches, we note that the above set of equations can be derived via the variation of the following minisuperspace action This action can be obtained from (1) by dropping out the spatial, non-dynamical part of the integral, i.e.
The dynamics of the system is completely encoded in the reduced action Smin describing a system of finite degrees of freedom. We now continue with the classical analysis of the solution space for this system.
Due to the transformation law of the lapse function and the t derivatives, this equation is still covariant under transformations of the form t = f (τ ). Regarding it's solution, it can be solved with respect to either n or a. However, if another degree of freedom s is defined via s = ndt, the solution with respect to a assumes the elegant form: where c1, c2 constitute non-essential (absorbable via coordinate transformations) constants. This solution is called the Bertotti-Kasner space [86,87] and belongs to the family of Kantowski-Sachs spacetimes [88].

generic solution
Contrary to the previous case, we can now solve the constraint equation with respect to n: When this value of n is substituted into the rest of the equations, only one of the reduced equations is independent, and one way to express it is This equation is covariant as in the previous case. The above expression is simplified if we introduce the following degree of freedom w so that the previous equation becomesẅ These variables have been used before in [89] for the study of correlations across the horizons that this solution possesses. In these variables, the solution of (16) is rather simple In the above line element, the constant M is essential. Note that usually M is associated to a "mass" constant m0 = M 2 . Here for simplicity and in order to avoid more numeric factors we just use M for the constant of integration.

Lagrangian-Hamiltonian description
The starting point to construct a Lagrangian in order to describe the dynamics of the previous system is the Einstein-Hilbert action. As we can see from (8), such a Lagrangian exists and we have to find out its reduced form after the assumptions of Case 1 and the solution of the constraint in Case 2.

Case 1
There are two basic ways to treat the condition b(t) = 1 √ Λ : either as a constraint by introducing a Lagrange multiplier or by simply replacing it. The latter immediately yields Thus, the reduced Lagrangian is zero. Nevertheless, it is possible to construct a Lagrangian via which the equation (10) can be acquired. Since there is only one equation for two degrees of freedom (n, a) or (ṡ, a) equivalently, we can choose only one of them as independent degree of freedom and the other as a function of the parameter t. So, instead of a singular Lagrangian, a time-dependent, regular one is obtained which reads and mimics the Lagrangian of an inverse oscillator with "mass"= 1 n(t) and a "spring constant"= Λ 2 n(t). This Lagrangian is the unique, up to a multiplicative constant and the addition of a total derivative, quadratic in the "velocity"ȧ, whose Euler-Lagrange equations coincide with (10). Due to the covariance, the function n(t) can be chosen at will so it is not an essential "mass".
Regarding the canonical Hamiltonian, we define the momentum as usual Solving the previous relation with respect toȧ, the Hamiltonian is obtained via the Legendre transformation While H is not conserved due to the explicit time dependence of n(t), H0 is conserved. The Hamilton equations readȧ = p n(t),ṗ = Λ a n(t), and it can be proven that they are equivalent to the Euler-Lagrange equations.

Case 2
In this case, we can construct a reduced Einstein-Hilbert Lagrangian which is not zero. By substituting in the initial Lagrangian (8) the solution of the constraint (13), the former acquires the square root form: It is straightforward to check that the two Lagrangians are equivalent. The equation of motion (16) can be obtained via Euler-Lagrange equations either with respect to w while b(t) is considered a time-dependent function or the other way around. Since we want to quantize the system, using the canonical approach, it is to our interest to avoid a Lagrangian which depends on the square root of the velocity. It is not difficult, through fitting some initially arbitrary functions of w and b, to obtain the unique Lagrangian which is quadratic in the "velocity"ẇ and yields the same Euler-Lagrange equation We need to note that the classical equivalence between two Lagrangians, like (24) and (25), which lead to the same Euler -Lagrange equations, does not necessarily guarantee a quantum equivalence as well. This is a valid problem recognized by several authors [90,91]. It has also been noted that the process of quantization may be connected to the type of Lagrangian or Hamiltonian that is being used [92]. For example the canonical quantization has been developed by taking in consideration Hamiltonians which are in general quadratic in the momenta and thus may fail for Hamiltonians with a different functional dependence. This is the main reason, why we choose to make the passing from Lagrangian (24) to (25). In this manner we are going to have a quadratic Hamiltonian at our disposal. Thus, we may follow the typical process of canonical quantization avoiding the complex problem of assigning an appropriate operator to a Hamiltonian that has some different type of dependence in the momentum. It is also interesting to mention that the question of quantum equivalence even plagues the process of the classical reduction from the original gravitational action to the minisuperspace description. In other words, it is not guaranteed that the minisuperspace Lagrangians, which correctly reproduce the Einstein's equations after adopting an ansatz for the line element, have a quantum description which is relevant to the quantum behaviour of the original gravitational system. The equivalence however of the classical dynamics, leads us to the assumption that, even at the worst case scenario, at least some properties of the quantum minisuperspace system may be of relevance and help us to obtain insights with respect to the actual quantum gravitational configuration. Besides, this is the main approach we have at our disposal in the absence of a complete theory of quantum gravity.
The expression (25) can be further simplified by introducing another time dependent variable which transforms as a density under time reparametrizations, we thus define s(t) as where κ is a random integration constant. The Lagrangian is now written as and it is valid The horizons can be found from the positive roots of the equation w(t) = 0, where w(t) is given by (17). The relevant third order algebraic equation to be solved is From the theory of cubic equations [93], we know that (28) can have up to three real roots. In particular, if Λ < 0 there exists a single real root given by which is positive when M > 0. In the case where Λ > 0, there exists a single real root when 9ΛM 2 − 4 > 0 given by The latter is positive when M < 0. Finally, if Λ > 0 and at the same time 9ΛM 2 − 4 < 0 there exist three real roots The root b0 is always negative, the roots b± are both positive when M > 0, while for M < 0 only b+ is positive.
Behind the horizons the curvature singularity appears at b = 0. Regarding the Hamiltonian, it is simpler from what we had in the previous case and, as we see, we have two conserved quantities, H0 and p.

Time-Covariant Schrödinger equation
The canonical quantization procedure will be followed in this section: the canonical momenta and positions are replaced by operators (q i →q i , pi →pi) which satisfy the property of self-adjointness (given appropriate boundary conditions for the wave function), and the canonical commutation relations where [·, ·] C denotes the commutator, δ j l the Kronecker's delta, ψ|φ = d m q µ ψ * φ the inner product, where m is the dimension of the configuration space and µ a proper measure. A usual choice for µ is to take the square root of the absolute value of the determinant of the configuration space metric (the so called natural measure). For the one dimensional Lagrangians (20) and (27), which we derived earlier, the natural measure is a constant. Due to the invariance of the wave function up to a normalization constant, without loss of generality, we can simply take it as µ = 1. In the position representation of quantum mechanics the operators are chosen asp where q l = (a, w).
As in the work [83], in both problems set by Lagrangians (20) and (27) we consider t as an external time parameter and a Schrödinger equation can be constructed which is time-covariant, with respect to time reparametrizations, due to the transformation properties of n(t) and s(t) respectively.

Case 1
For the moment, we will restore the standard S. I. base units to keep a better track of the scales involved. The Lagrangian of this case reads L =ȧ Note that, in our initial considerations, through which we derived L, we adopted the units c = 1 and thus considered x 0 = t. Here, in order to restore the S. I. units, we introduced again the time coordinate as, t = x 0 /c, so theȧ in L becomeṡ a/c. For the Lagrangian to acquire units of energy, we multiply it with mc 2 λ 2 , where m is some mass and λ some length unit. By also introducing the variable q = λ a we have the Lagrangian and Hamiltoniañ where ω 2 = Λc 2 and has units of frequency. The time-covariant Schrödinger equation has the form This is exactly the equation for an inverted harmonic oscillator of frequency ω = iω. It has been extensively studied by previous authors, see for example [94], [95]. The general solution can be expressed in terms of the Parabolic Cylinder D functions, where v = − 1 2 − iE ω , z = (1 + i) mω q andv,z their complex conjugates. If we turn to the classical solution (12), we see that a may attain -depending on the constants of integration -any value from −∞ to +∞. This is the domain of definition which we assign to q as well. From a physical perspective though it can be argued, due to the metric being invariant under the change a → −a, that this leads to calculating two times the same trajectories. If we restrict a (and thus q) in the half-line R+, the differential operator ∂ 2 ∂q 2 appearing in the Hamiltonian is not essentially self-adjoint. However, it admits a one-parameter family of self-adjoint extensions; the necessary conditions on the wave function are where ν is a constant parameter which models the behaviour at the boundary q = 0 [96,97]. This obviously sets conditions over the integration constants µ1 and µ2. The well known Dirichlet and Neumann boundary conditions are obtained for ν = 0 and ν = ∞ respectively. For the subsequent calculations that we make, which are centered on the propagator describing the time evolution of some initial state, the consideration of the half line -at least for an appropriate boundary condition -does not affect our result. We discuss this point later in Section 4.1. Thus, we consider the whole R as the domain of a and q, so as to be in direct contact with the known results of the inverted oscillator and to cover all accepted solutions described by (12); even if they are duplicate from the metric perspective, where just the a 2 appears. Still, in the region (−∞, +∞), the wave functions (41) are not square-integrable 1 and the energy spectrum of the Hamiltonian (39) is continuous, varying from minus to plus infinity. Due to the explicit dependence on the time appearing multiplicatively as a common factor in the Hamiltonian through n(t), the former commutes with itself at different times. As a result, no time ordering needs to be considered for the time evolution operator. The constants µ1, µ2 are subject to normalization, while the constant E represents the "energy" of the particle. In the next section we will focus on some specific application for these wave functions.
Before proceeding we should comment on the arbitrariness of the time parameter appearing in (41). If we turn our attention to the action of the original minisuperspace system, eq. (7), we know that its parametrization invariance implies a symmetry generator of the form X = χ ∂ ∂t −χn ∂ ∂n , with χ = χ(t) an arbitrary function of time [98]. The first part of the generator tells us that arbitrary time transformations are in order, while the second encodes the information of the transformation law effected on the lapse by this reparametrization, which results in n(t)dt transforming as a scalar. This symmetry is a remnant of the full diffeomorphism invariance of the full theory [99]. We observe that the wave function (41) carries this remaining invariance of the minisuperspace system since its time dependence involves the integral of the factor n(t)dt. Of course it could be argued that one might reach the same result by applying an appropriate transformation over the standard time dependence, e − iE t , of a typical Schrödinger equation. However, that would be enforcing manually a symmetry, which we know that the original system has. In the process we follow here, this property is inherited naturally as a result of solving (40) without having chosen a particular time gauge. It is in this sense that we refer to an equation like (40) as a time-covariant Schrödinger equation. As it is expected, the solutions and thus any kind of measurable quantity, for instance the probability current, is invariant under time reparametrizations and thus equivalent for different "observers".

Case 2
The units redefinitions are applied here as well. The Lagrangian for this case is L =ẇ Even though the wave functions are not square-integrable in q ∈ (−∞, +∞) they are "normalizable" up to a delta function [94]. time coordinate, multiplying with mc 2 and defining q = w for the sake of simplicity, the Lagrangian and Hamiltonian readL The corresponding time-covariant Schrödinger equation becomes Note that from equation (15) it seems as if the variable w has the same sign as b. Since we have assumed for simplicity (without loss of generality) b > 0, this would imply that w > 0. However, in the classical regime the solution a 2 < 0 is acceptable, thus even for b > 0 we may have w < 0. Therefore, there is no reason to restrict ourselves, meaning that w is defined from minus to plus infinity. The situation here is quite different than what happens with the domain of q in the previous section and the solution of Case 1. As we discussed previously, the q < 0 of (38) is, from the gravitational perspective, duplicate to q > 0. Here, the q < 0 and q > 0 (or w < 0 and w > 0 respectively) signify different regions of the same solution which are parted by horizons. Thus, we are dealing with distinct parts and not double copies of the same solution. Thus, we have no reservation to consider q running in the whole real line. Under this condition, equation (45) expresses a free particle and the general solution can be written as where ν1, ν2 constants subject to normalization and as before, E is the "energy" of the particle.

Invariant Probability of persistence
We now attempt to make a deeper analysis on the wave functions obtained through the time-covariant Schrödinger equation. Our purpose is to study the evolution of initial states for the previously described systems and provide some interesting characteristics regarding the difference between the classical and quantum descriptions.

Case 1
As it is evident from the previous sections, the equations describing the evolution of the scale factor coincide with those of the inverted harmonic oscillator. This characteristic evolution appears in cases when someone is interested to describe some unstable equilibrium. For instance, in our example, the scale factor follows an exponential dependence on time, since we have found that the line element reads Thus, for n(t)dt monotonically increasing, if c1 = 0 the universe would grow exponentially, while for c2 = 0 the universe would decrease exponentially and would require infinite time to reach a size of zero. On the other hand, the unstable equilibrium a = 0 ⇔ c1 = 0, c2 = 0 could persist for infinite time, which is classically acceptable as far as one stays in the inverted harmonic oscillator interpretation of the dynamics of the system. In what the geometry is concerned, the value a = 0 renders the metric (47) non invertible. All fourteen curvature scalars are polynomials of Λ and thus constant. Space-times with this property are called Constant Scalar Invariants (CSI) spaces [100]. We now proceed to study the problem and extract results based on the similarity with the inverted harmonic oscillator equation. Quantum mechanically, the most we could do to describe the unstable equilibrium point, is to consider an initial state of minimum uncertainty localized around this point and study it's evolution, as well as how much time it is required for the universe to enter the exponential phase. To do so, we will follow [94] with some changes which will be discussed below.
The starting point is the initial state, which is taken to be a Gaussian, as well as the propagator of the inverted harmonic oscillator, where l represents the width of the Gaussian wave-packet, B(t) = n(t)dt and from now on the time dependence will be omitted. It can be verified that this propagator satisfies the equation (40) and the fact that for the initial state the uncertainty relation becomes minimum ∆q∆p = 2 . Regarding the constant l, it is the standard deviation and some possible values will be discussed later on. The evolution of the initial state is given by where the following quantities have been defined The probability density and the uncertainty relation acquire a simple form For the latter calculation we have used the well known relations for the root mean square deviations in which the mean value for an observable A is defined as A = µΨ * AΨdq. Remember that in this case µ = 1. Before proceeding, let us comment at this point that if we had adopted the half line as the domain of a (and hence q). Then, we could use an image-like method [101,102] to writẽ K(t, q,q) = K(t, q,q) + K(t, −q,q).
It is easy to see that such a propagator satisfies the Schrödinger equation (40) for a Neumann boundary condition ∂qΨ(t, 0) = 0, which corresponds to having ν = ∞ in (42). Then, the integral (50) with K(t, q,q) in the real line R and the corresponding expression withK(t, q,q) in the half-line R+ yield the same result for Ψ as we see it in (51). We need to recognize however that this correspondence appears with this specific "reflective" boundary condition. For different values of ν we expect different propagators and different results. Of course, even though we obtain the same Ψ of (51) for this particular boundary condition, the integrals involving calculations that use this wave function have different domains of integration in the two cases. For the particular quantities regarding the time intervals which we calculate later, this just translates into a constant factor since the Ψ of (51) is an even function in q. However, this factor is insignificant when considering the order of magnitude of these intervals as is going to be revealed later. It is thus our choice to remove altogether the ambiguity of ν in (42) by considering the full mathematically accepted domain of definition for a ∈ R, based on the classical solution.

Invariant Non-decay probability and Mean lifetime.
The probability that a system prepared at an initial state Ψ(0, q) has not yet decayed at the state Ψ(t, q) is calculated as follows At this point, we define the invariant mean lifetime as or by introducing the variable η = ωB As it was argued in [94], it is useful to give another realization of the decay rate via the following description: The probability Q, that a particle remains still within a distance sp is given by where ξ = q r h Ω and ξ0 = sp r h Ω . In accordance to the previous case, we define an invariant mean sojourn time as It is rather helpful to define the following integration variable ζ = l r h Ω such that dζ = −ζ 1 Ω dΩ dB dB and 1 Ω dΩ dB = ω 1 − ζ 2 1/2 1 + r h l 4 ζ 2 1/2 in order to have where Erf sp l ζ the Error function defined as Erf sp l ζ = 2π −1/2 sp l ζ 0 dξ e −ξ 2 . Both TM and TMs depend on the "frequency" ω as we may have expected, since by definition ω is related to the only true length scale of the system 1 √ Λ . One of the differences is the introduction of an additional length scale sp in the second case, which might be related to the resolution distance of some apparatus.

Case 2
Let us now turn to the second case, where the Lagrangian resembles that of a free particle. The configuration for which a = 0 ⇔ q = 0 corresponds to the horizon surface. Thus, we could follow the same procedure as previously, to find out how a configuration representing the horizon surface would evolve due to quantum mechanical effects. The same initial state is used as previously, while in this case, the propagator and the evolved state read where ωp = ml 2 , Ω 2 = 1 + (ωpB) 2 and B(t) = s(t)dt. The uncertainty relation becomes ∆q∆p = 2 Ω > 2 .

Invariant Non-decay probability and Mean lifetime
Following the steps of the previous subsection, the desired quantities have the form and so the mean life time becomes where we used ηp = ωpB. The fact that the mean lifetime is infinite indicates that the initial Gaussian state is very stable. We comment on possible implications of this later in our analysis. We proceed with the calculation of the mean sojourn time. By introducing the additional length scale, we acquire where ξ = q lω and ξ0 = sp lΩ . With the help of the variable ζp = 1 Ω , the mean sojourn time reads (71)

Non-decay time and Mean sojourn time, Case 1
In this subsection, some estimates will be provided for the non-decay time of the previously studied cases. Let us start from equation (61) Among the involved constants we have the mass m which we introduced for dimensional reasons, the Planck constant = 1.1 * 10 −34 J · s, the speed of light c = 3 * 10 8 m/s, the standard deviation l, the cosmological constant Λ and the "frequency" ω 2 = Λ c 2 . At this point we assume the standard deviation of the original wavepacket to be of the order of the Planck length l = lp = 1.6 * 10 −35 m, where the quantum gravity effects are expected to appear. For the mass constant m we choose a value which corresponds to the minimum energy content of the early universe ∼ 10 19 GeV = 1.6 * 10 9 J. Thus, we use the mass we find by setting mc 2 = 1.6 * 10 9 , which yields m = 1.8 * 10 −8 kgr.
This is close to the order of magnitude of the Planck epoch ∼ 10 −43 s. Of course this number refers to estimates based on a FLRW geometry, thus any comparison is to be taken with a precaution. We may assume that, the quantum packet was highly concentrated for times t < TM around the value zero, while for t ≥ TM it starts decaying, leading at some point to an exponential dependence of the scale factor and thus driving some sort of inflationary epoch. How can we understand that the scale factor will acquire an exponential dependence? If the classical degree of freedom q (or equivalently a) is connected to the mean value of the operatorq, then for the specific wave packet is equal to zero. However, we can obtain a non-zero expression based on the following estimate, q 2 class =<q 2 > where we substituted the Ω as defined from (52). This form of solution mimics the classical one (47) for 2l (remember that η = ωB = ω n(t)dt). Thus, as time passes, the scale factor grows exponentially.
To see how other values affect the time TM , let us take the value for the cosmological constant obtained via the standard model, Λ ∼ 10 70 m −2 , which is 10 122 orders of magnitude larger. The time then is TM ∼ 2.6 * 10 −44 s, which is smaller than (74) and well inside the Planck epoch. If we keep Λ as it is, and reduce the standard deviation, say take l = 9.7 * 10 −2 lp, we obtain exactly the Planck time TM ≃ 5.4 * 10 −44 s. Finally, with this l and the cosmological constant from the standard model we obtain TM ≃ 2 * 10 −45 s. The important thing to notice is that all the relevant scales, lead to values of time in the domain where the quantum gravity effects are expected to be strong.
Regarding the mean sojourn time, it is given by the equation (64) An additional length scale was introduced, it is the length within which the wave packet is constrained, which we will assume to be some multiple β of the Planck length. Furthermore, since the standard deviation is assumed to be equal to the Planck length, β counts the number of times that standard deviation fits within the length: We provide the sojourn time for two values, β = 2, β = 4: β = 4 ⇒ TMs ≃ 2.8 * 10 −41 s.
As we can observe, the closer the value of β to unity the closer the sojourn time to the non-decay time. The sojourn time TMs for the inverted harmonic oscillator, and for a Gaussian wave packet, is decreasing for larger frequencies ω. This is a known result in the theory of quantum mechanics. It is interesting to note however, that if the notion of non-commutativity of space is introduced, then this changes and TMs becomes a concave function of ω, for further details on this see [104].

Non-decay time and Mean sojourn time, Case 2
For the case of the horizon decay, we follow pretty much the same procedure , with the same values. The parameter involved reads ωp = 4.3 * 10 35 s −1 . As we previously noted, the mean non-decay time here is infinite (see equation (69)). This seems to indicate that the initial state is extremely stable. We could infer that the previous case, which results in a TM in the range of (74), is better suited for our purposes in the sense that it allows the universe to "escape" the initial state in a finite time period and to possibly transit to different phases of evolution. When comparing the two decay probabilities we observe that they have similar behaviours at initial times, i. e. 1 − P ∝ t 2 for both (59) and (68) with B ∝ t << 1. This initial behaviour is referred to as a Zeno period in the literature [105] as such time dependence leads to the known Zeno effect of Quantum Mechanics [106]. At late times the two systems differentiate, the non-decay probability of equation ( We see that this calculation yields the same order of values for TMs as that of the previous case. We thus get a completely different picture by looking at this quantity when comparing to the non-decay time. It has been argued in [94], for the system studied there, that the mean sojourn time can be considered more realistic than the non-decay time since the former takes into account an external length scale in its calculation (namely the sp of (62)).

The minisuperspace canonical quantization
In this section, we briefly derive the wave functions that are obtained through the usual canonical quantization of the minisuperspace Lagrangian LEH of (8) by solving the relevant Wheeler-DeWitt equation. To a large extent the main formalism follows the procedure used in [107] for the canonical quantum description of the Schwarzschild space-time.
For reasons that will become obvious later in our analysis, and which have to do with the construction of a typical eigenvalue problem, we multiply the Lagrangian LEH with a constant V0 to which we assign the description of a finite spatial volume. This is owed to the fact that the original Einstein -Hilbert action contains an integral over the spatial region, e.g. see equation (9). We may take the latter as a constant by considering an integration in a finite region of the spatial section. If we start from Lagrangian V0LEH and perform a reparametrization of the lapse function n(t) → e(t) as we are led to a new Lagrangian The latter is equivalent to the initial Lagrangian, since the change n(t) → e(t) just represents a time reparametrization. The configuration space is two dimensional, spanned by the q i = (a, b), and the ensuing minisupermetic is Thus, the Lagrangian Le can be interpreted as describing the motion of a free particle of mass V0 in a spacetime characterized by the metric Gij . The e in this formalism is called the einbein field [108]. The resulting Hamiltonian constraint for the theory is given by where " ≈ " denotes a weak equality in Dirac's theory of constrained systems [109]. The Wheeler-DeWitt equation emerges by enforcing the quantum version of the Hamiltonian constraint upon the wave function as HΨ = 0. For the factor ordering of the kinetic term of (86), the optimal choice in what regards invariance of the probability under transformations in the configuration space and time reparametrizations like (83), is the conformal Laplacian, which however in two dimensions is trivially equivalent to the usual Laplacian. We can thus write where µ(q) = −det(Gij ) is the natural measure that appears in the definition of the inner product We notice that the metric Gij , given by (85), describes a flat two dimensional geometry of hyperbolic signature. The relevant canonical quantization in such a space has been used before in several problems in the literature [33,110,111]. As basic observables, together with the Hamiltonian constraint, we can either use the commuting momenta in Cartesian-like coordinates or the boost in the relevant plane (in a similar manner to what is done with the angular momentum quantization in spaces with Euclidean signature). Thus, the choice of the solutions to the Wheeler-DeWitt equation is made by the use of the linear in the momenta integrals of motion of the classical system as basic observables. In other words, we construct operators based on the classical integrals of motion and use them to define eigenvalue equations to be used in conjunction with the Hamiltonian constraint.
1. Quantum description with respect to the momenta. For the shake of simplicity we work in light-cone coordinates. Under the transformation (a, b) → (u, v) where the metric (85) becomes Gij dq i dq j = 2dudv. The minisuperspace Lagrangian is written as and the Hamiltonian constraint reads where pu and pv are the momenta in the null directions. The quantization in these variables is straightforward. We consider pu = −i ∂ ∂u , pv = −i ∂ ∂v , while the u, v just act multiplicatively. The classical constants of motion pu, pv imply the quantum eigenvalue equations which result to the well known plane wave solution Ψ = 1 2π e i (kuu+kvv) , where we have directly put the needed multiplicative constant for normalization. The subsequent enforcement of the Hamiltonian constraint HΨ = 0 yields the following relation between the constants ku and kv which are thus related to the spatial volume.
2. Quantum description with respect to the boost. The third, linear in the momenta, classical integral of motion that the system (91) possesses is the boost L0 = upu−vpv, whose quantum equivalent can be written as L0 = −i (u∂u − v∂v). It is however simpler if we go over to coordinates where L0 is brought into normal form. This happens in the new variables (u, v) → (ρ, ϕ) where The minisuperspace metric becomes now Gij = diag(−1, ρ 2 ) and the boost operator is transformed to L0 = −i ∂ ∂ϕ . The relevant eigenvalue equation, results into a wave function of the form Ψ(ρ, ϕ) = 1 √ 2π e i λϕ ψ(ρ). The unspecified function ψ(ρ) is found with the help of the Wheeler-DeWitt equation that, by the use of (87), leads to which is the Bessel equation with generic solution with Jµ(z), Yµ(z) denoting the Bessel functions of the first and second kind respectively. We need to be careful in choosing the function that solves (96). If we express how the new variables (ρ, ϕ) are related to the original (a, b) of the minisuperspace, we have: In the region where Λb 2 − 3 ≥ 0 our variables are real. In this case, as a solution to (96) we can choose the function ψ(ρ) = ψ λ, where Re denotes the real part of the expression in the brackets. This function was firstly introduced and studied in [112]. It is normalizable up to a delta function given the constraint set by the Hermiticity of the Hamiltonian operator V0/V ′ 0 = e 2kπ λ, k ∈ Z [111]. A different linear combination however of the solution can be chosen that evades conditions on the "eigenvalue" V0 [33,44]. Note that the weight ρ inside the integral is what originates from the natural measure in these coordinates since µ = −detGij = ρ.
On the other hand, if we consider the case Λb 2 − 3 < 0, then in the above description ρ, ϕ become complex. For ϕ this is not a problem since the second of (98) is written as ϕ = − iπ 2 +ln At the level of the wave function, whose dependence on ϕ is Ψ ∝ e i λϕ , this results in an extra multiplication constant which can be absorbed through the process of normalization. Things however are different for the variable ρ which becomes imaginary and now the previously used Bessel combination cannot be normalized even up to a Dirac delta function. This problem is bypassed by adopting a different linear combination of the solution to the Bessel equation. By considering ρ = iρ, where nowρ ∈ R, we may write the solution with respect to the modified Bessel function of the second kind As before the C denotes the relative normalization constant. The K iλ (V0ρ) is a real function and the inner product yields the integral [113] which is finite in the limit V ′ 0 → V0. The critical limit Λb 2 − 3 = 0, which separates the two regions leading to different wave functions, is the outermost horizon which appears in the classical solution when the constant of integration M is zero in (18).

Discussion
This work was dedicated to the study of the classical and quantum minisuperspace description for the Λ-Kantowski-Sachs universe. In the classical regime, by following the method developed in [83], regular, albeit time dependent Lagrangians (instead of the original singular Lagrangians) were constructed, which encode the dynamics for both branches of the solution space: The Bertotti-Kasner and the generic solution (the cosmological equivalent of Schwarzschild-de Sitter space).
The regular nature of the Lagrangian and the appearance of the external time parameter t, enabled the construction of a time-covariant Schrödinger equation. Covariant in the sense that the gauge freedom t = f (t) was not "broken" and hence different "observers" will, by default provide equivalent quantum descriptions. The general solutions to these equations were obtained and expressed in terms of elementary functions. Moving a step further from the work [83], a time-invariant non-decay probability was provided, alongside with the non-decay and mean sojourn time. Contrary to the usual definitions of the above objects, which are invariant at most under translations t =t + ǫ, in our case the freedom is rather generic t = f (t).
The non-decay and mean sojourn time mentioned in the previous paragraph, were calculated for the following schemes: Regarding the Bertotti-Kasner solution, an initial Gaussian wavepacket was assumed to be concentrated around the value a = 0 (where a is the scale factor) and its evolution was obtained. By assuming the width of the wave packet to be of the Planck order and the observed value for the cosmological constant, the decay time was calculated to be in the region of 10 −42 s. This is at the end of the Planck epoch. Furthermore, by choosing the value of the cosmological constant derived from the Standard Model of Particles, the decay time reduces to almost ∼ 10 −44 s, which is still inside the region where quantum gravity effects are expected to be significant.
For the generic solution, a Gaussian wavepacket was assumed to be concentrated around the horizon value for the scale factor. The non-decay time becomes infinite while the relative sojourn time provides a completely different picture. The order of magnitude for the latter, which is expected to be better suited for physical interpretation, is of the same order as the one of the special solution at ∼ 10 −41 s.
The first wave function corresponds to having the momenta of the minisuperspace as quantum observables in Cartesian-like coordinates and as a solution it was given in a previous work by J. Louko and T. Vachaspati [60]. The last two expressions refer to quantization with respect to the boost, given for different regions of the classical solution, so that the relative variables remain real through the transformations we have performed. Interestingly enough, the wave function (102) is comparable to (46). The latter is also a plane wave solution in the variable q = w which was used in section 3.1.2. This is revealed through the following considerations: Let us notice that the first of (89) becomes u = 2q, if we introduce the relation a = √ wb = √ qb. At the same time, the classical solution in these light-cone coordinates is written, in a parametrization invariant way, as being linear in e(t)dt. We may thus re-write (102) as where we have used the, compatible with the classical solution, expression v = e(t)dt, i.e. we effectively consider v as the time parameter of the problem, and we additionally substituted kv from (93). The wave function (105) is essentially the same as that we obtained through a different process in (46). Of course the e is not the same as the n appearing in (46), but neither is the t variable that is being integrated since the two functions are related through a time reparametrization. Given this parametrization invariance of the problem we are allowed to make this comparison and relate the two results at a semiclassical level (in the sense that we utilized information from the classical solution). As far as (103) and (104) are concerned, which regard the quantization with respect to the boost symmetry of the problem, we do not see them being recovered with the previous method. However, the analysis performed has the big advantage of offering the possibility of a quantum description for the special solution of Bertotti and Kasner, for which we have no minisuperspace approximation through (84) since Gij = 0.