Quantum Cosmology of $f(R,T)$ gravity

Modified gravity theories have the potential of explaining the recent acceleration of the Universe without resorting to the mysterious concept of dark energy. In particular, it has been pointed out that matter-geometry coupling may be responsible for the recent cosmological dynamics of the Universe, and matter itself may play a more fundamental role in the description of the gravitational processes that usually assumed. We study the quantum cosmology of the $f(R,T)$ gravity theory, in which the effective Lagrangian of the gravitational field is given by an arbitrary function of the Ricci scalar, and the trace of the matter energy-momentum tensor, respectively. For the background geometry we adopt the Friedmann--Robertson--Walker metric, and we assume that matter content of the Universe consists of a perfect fluid. We obtain the general form of the gravitational Hamiltonian, of the quantum potential, and of the canonical momenta, respectively. This allows us to formulate the full Wheeler-de Witt equation describing the quantum properties of this modified gravity model. As a specific application we consider in detail the quantum cosmology of the $f(R,T)=F^0(R)+\theta RT$ model, in which $F^0(R)$ is an arbitrary function of the Ricci scalar, and $\theta $ is a function of the scale factor only. The Hamiltonian form of the equations of motion, and the Wheeler-de Witt equations are obtained, and a time parameter for the corresponding dynamical system is identified, which allows to formulate the Schr\"{o}dinger--Wheeler--de Witt equation for the quantum-mechanical description of the model under consideration. A perturbative approach for the study of this equation is developed, and the energy levels of the Universe are obtained by using a twofold degenerate perturbation approach. A second quantization approach for the description of quantum time is also proposed, and briefly discussed.


I. INTRODUCTION
One of the cornerstones of theoretical physics, General Relativity (GR), formulated mathematically in terms of the Einstein field equations, proved to be a very successful gravitational theory at the scale of the Solar System. By using GR we can describe the gravitational dynamics of the Solar system with a high precision, and phenomena like the Mercury perihelion precession, the bending of light while passing the Sun, and the gravitational redshift can be fully understood. In GR the gravitational field equations can be obtained by varying the Einstein-Hilbert action S = 1 16π R √ −gd 4 x + L m √ −gd 4 x, where R is the Ricci scalar, and L m is the matter Lagrangian, with respect to the metric g µν , and they read as R µν − 1 2 g µν R = 8πT µν , where T µν is the energymomentum tensor. The energy-momentum tensor T µν identically satisfies the mathematical relation ∇ µ T µν = 0, which can be interpreted from the physical point of view as the energy conservation. Essentially GR is a beautiful geometric theory that establishes a deep connection between the geometry of the spacetime, matter fields, and gravitational interaction. Considering larger scales of the Universe, using GR we can numerically simulate galaxies' formations and collisions, and the results of these simulations can be verified by the increasingly trustworthy data obtained due to the rapidly improving observational techniques. The predictions of general relativity have also been confirmed in the strong gravity regime by the discovery of the gravitational wave emission in the binary pulsars system PSR 1913 + 16 [1], a discovery that has opened a new testing ground for GR and for its generalizations. The detection by LIGO of GW150914 from the inspiral and merger of a pair of black holes [2] will start a new era in observational astronomy, based on the theoretical and mathematical formalism of GR.
Thus, with a large number of astronomical observations, as well as terrestrial experiments testing and confirming it both at the weak and strong gravity regime, it would seem that GR gives a full description of the gravitational interaction at the non-quantum level. How-ever, several recent astrophysical and cosmological observations have raised the intriguing possibility that GR may not be able to model and explain the gravitational dynamics at scales much larger that the one of the Solar System.
On a fundamental theoretical level the two most important challenges GR must face are the dark energy and the dark matter problems, respectively.Several cosmological observations, obtained initially from the study of the distant Type Ia Supernovae, have provided the unexpected result that the expansion of the Universe did accelerate lately [3][4][5][6][7]. The paradigmatic and usual explanation of the late time acceleration requires the existence of a mysterious and dominant component of the Universe, called dark energy (DE). Dark energy is responsible for the late time dynamics of the Universe [8,9], and can explain the observed features of the recent cosmological evolution. The second and equally mysterious component in the Universe, called Dark Matter, an assumed non-baryonic and non-relativistic "substance", is necessary for the explanation of the flat rotation curves of galaxies, and for the virial mass discrepancy in clusters of galaxies [10,11]. The detection/observation of dark matter is restricted by the fact that it interacts only gravitationally. Its effects can be observed by observations of the motion of the massive hydrogen clouds around galaxies, or by the motion of the galaxies in clusters [10]. However, despite many decades of intensive observational and experimental efforts the particle nature of dark matter still remains essentially unknown.
One interesting possibility for explaining dark energy is based on theoretical models that contain a mixture of cold dark matter and a slowly-varying, spatially inhomogeneous component, called quintessence [12]. The idea of quintessence can be implemented theoretically by assuming that it is the energy associated with a scalar field Q, having a self-interaction potential V (Q), and a pressure p =Q 2 /2 − V (Q) associated to the quintessence Q-field. Such a model also allows a possible theoretical interpretation in terms of particle physics results. If the potential energy density V (Q) of the quintessence field is much greater than the kinetic one, then it follows that the pressure p of the field is negative. Quintessential cosmological models have been extensively investigated in the physical literature (for a recent review of quintessence cosmologies see [13]).
Alternatively, the recent acceleration of the Universe can also be explained by scalar fields φ that are minimally coupled to gravity via a negative kinetic energy, known as phantom fields, which have been proposed in [14]. Interestingly enough, these cosmological models allow values of the parameter of the equation of state w with w < −1. Hence, real or complex scalar fields may play a fundamental role in the cosmological processes describing the evolution of our Universe, and they may provide a realistic description of the observed cosmic dynamics.
However, in order to explain the observed gravitational dynamics of the Universe a different line of thought on dark energy was also considered. It is based on the fundamental idea that dark energy is not a particular physical field, but it can be understood as a gravitational phenomenon induced at cosmological scales by the intrinsic modifications of the gravitational interaction itself. Hence the fundamental assumption of this theoretical approach is that at large astrophysical and cosmological scales standard general relativity cannot describe correctly the dynamical evolution of the Universe. In this context many modified gravity models, all trying to extend and generalized GR, have been proposed. Historically, in going beyond the standard gravitational models, the first step was to generalize the geometric part of the Einstein-Hilbert action. One of the first models of this type is f (R) gravity, in which the gravitational action is an arbitrary function of the Ricci scalar R, so that [15][16][17][18][19]. However, this and many other modifications of the Einstein-Hilbert action concentrate only on the geometric part of the action, by implicitly assuming that the matter Lagrangian plays a subordinate and passive role only [20], which naturally follows from its minimal coupling to geometry. But a fundamental theoretical principle forbidding a general coupling between matter and geometry have not been formulated yet, and in fact it does not exist a priori. If such matter-geometry couplings are allowed, many theoretical gravitational models with extremely interesting properties can be constructed.
The first of these kind of models was the f (R, L m ) modified gravity theory [21][22][23][24], with a gravitational action of the form S = 1 16π f (R, L m ) √ −gd 4 x. A similar geometry-matter type coupling is assumed in the f (R, T ) [25,26] gravity theory, where T is the trace of the energymomentum tensor. For a recent review of the generalized f (R, L m ) and f (R, T ) type gravitational theories with non-minimal curvature-matter coupling see [27]. Several other gravitational theories involving geometry couplings have also been proposed, and extensively investigated, like, for example, the Weyl-Cartan-Weitzenböck (WCW) gravity theory [28], hybrid metric-Palatini f (R, R) gravity, where R is the Ricci scalar formed from a connection independent of the metric [29,30], f (R, T, R µν T µν ) type models, where R µν is the Ricci tensor, and T µν the matter energy-momentum tensor, respectively [31,32], or f (T , T ) gravity [33], in which a coupling between the torsion scalarT and the trace of the matter energymomentum tensor is assumed. For a review of hybrid metric-Palatini gravity see [34].
Modified gravity models with geometry-matter coupling are important since they can provide from a fundamental theoretical point of view a complete theoretical explanations for the late time acceleration of the Universe, without postulating the existence of dark energy. They can also offer some alternative explanations for the nature of dark matter. Moreover, these models show that matter itself may play a more fundamental role in the description of the gravitational dynamics that usually assumed [20], and they can also represent a bridge connect-ing the classical and the quantum worlds. For example, the dependence of the gravitational action on the trace of the energy-momentum tensor T may be due to the presence of quantum effects (conformal anomaly), or of some exotic imperfect quantum fluids [27].
Besides the difficulties presented by the present day cosmological observations, a central theoretical problem in present they physics is the unification of quantum mechanics and gravitation. Gravitation dominates the dynamics of objects at large scales, while quantum mechanics describes the microscopic behaviors of the particles. The study the Universe as a whole from the quantum mechanical point of view is the subject of quantum cosmology [35,36], which is based on the idea that quantum physics must apply to anything in nature, including the whole Universe. An unification of the electromagnetic force, of the strong force and of the weak force, respectively, is achieved in the standard model of particle physics, leaving the gravitational force as an exception that cannot be yet unified with the other fundamental forces. This is related to the fact that when considering gravitation in the framework of general relativity, we must consider not only matter, but also space and time, as physical objects. Space and time obey dynamical laws, and they have excitation such as gravitational waves that interact with each other and with matter. These aspects make quantizing the Universe far from being straightforward. Since the formation of cosmic structures is strongly dependent on the space-time interaction, quantum cosmology is therefore closely related to quantum gravity, representing the quantum theory of the gravitational force and of space-time [37].
Even being of a speculative and controversial nature, having several difficult conceptual problems to overcome, quantum cosmology has a long history [38][39][40], and various popular competing attempts have been proposed to quantized the gravitational field, like, for example, string theory, canonical quantum gravity and loop quantum gravity [35,36]. However, the lack of related observations reduces our abilities to resolve conceptual issues to all practical purposes. Since at the beginning of the Universe the average radius of each point is infinitely small, while the geometric curvature is infinitely large, quantum gravitational effects will dominate the dynamics of the Universe, and therefore they cannot be neglected in the study of the very early Universe. One of the main obstacles in the understanding of quantum cosmology is the so-called problem of time, which comes from the Hamiltonian constraint in the Arnowitz-Deser-Misner (ADM) formalism, leading to the Wheeler-de Witt equation [38][39][40], a fundamental equation in canonical quantization of cosmology. The quantum cosmology of f (R) gravity theory with Schutz's fluid is discussed in [41,42], with new perspectives to the time problem in quantum gravity considered.
It is the purpose of the present paper to study the quantum cosmology of the f (R, T ) gravity theory. Classical aspects of this theory have been extensively inves-tigated [43][44][45][46][47][48][49][50][51][52][53][54][55], but its quantum implications have not been considered yet. In order to construct the quantum cosmology of the f (R, T ) gravity we adopt for the classical background metric the Friedmann-Robertson-Walker form, and we assume that the matter content of the very early Universe consists of a perfect fluid, described by two thermodynamic parameters only, the energy density, and the thermodynamic pressure, respectively. In order to introduce the canonical quantization scheme for the f (R, T ) gravity theory, as a first step in our study we obtain the general form of the gravitational Hamiltonian, of the quantum potential, and of the canonical momenta, respectively. Once these quantities are explicitly found we write down the full Wheeler-de Witt equation of the f (R, T ) modified gravity theory, which describes the quantum properties of the very early Universe, when quantum effects had a dominant influence on the dynamic evolution of the system. We introduce and consider in detail the quantum cosmological properties of a particular model of the f (R, T ) theory, namely, the quantum cosmology of the f (R, T ) = F 0 (R) + θRT model, in which F 0 (R) is an arbitrary function of the Ricci scalar, and θ is a function of the scale factor only. We obtain the Hamiltonian form of the classical equations of motion for this model, and then we write down the Wheeler-de Witt equation, describing the evolution of the wave function of the early Universe. Starting from the Wheeler-de Witt equation we introduce a time parameter for the corresponding quantum dynamical system, which allows us to formulate the Schrödinger-Wheeler-de Witt equation for the quantum-mechanical model under consideration. In order to obtain the properties of the early Universe we develop perturbative approach for the study of this cosmological equation, and the energy levels of the Universe are obtained by using a twofold degenerate perturbation approach. Finally, we address the problem of the quantum time by introducing a second quantization approach for the description of the time in the quantum cosmology of the f (R, T ) modified gravity theory.
The present paper is organized as follows. The Hamiltonian formulation of the f (R, T ) gravity theory is presented in Section II, where the canonical momenta associated to the field variables as well as the potential are obtained. This allows us to write down the Wheeler-de Witt equation. The quantum cosmological formulation of the particular model f (R, T ) = F 0 (R)+θRT is presented in Section III. In particular, we identify the canonical momentum associated to the time variable, which allows us to transform the Wheeler-de Witt equation into an equivalent Schrödinger-Wheeler-de Witt equation. The solutions of this Schrödinger-Wheeler-de Witt equation are studied, by using a perturbative approach, in Section IV. We consider the problem of the second quantization of time in Section V. In Section VI we discuss and conclude our results. The derivation of the field equations of the f (R, T ) gravity theory, of the energy balance equation, as well as the canonical momenta for the f (R, T ) = F 0 (R) + θRT model, and the expressions of the Ricci tensor for the Friedmann-Robertson-Walker metric are presented in Appendices A-D, respectively.

II. THE WHEELER-DE WITT EQUATION IN f (R, T ) GRAVITY THEORY
We begin our study of the quantum cosmological aspects of the f (R, T ) gravity theory by briefly presenting the classical field and conservation equations. Then we obtain the Hamiltonian formulation of the theory, which allows us to write down the Wheeler-de Witt equation in the f (R, T ) theory, which describes the quantum evolution of the very early Universe in the presence of geometry-matter coupling. In the present paper we use the natural system of units with c = = 16πG = 1.
A. The f (R, T ) gravity theory The gravitational action for the f (R, T ) gravity model is [25] We define the energy momentum tensor of the matter as By assuming that the matter Lagrangian L m does not depend on the derivative of g µν , we obtain The field equations of our model can then be obtained as (see Section A for details) where Θ µν ≡ g αβ δT αβ /δg µν . Contracting the above field equations we find For perfect fluid, in the co-moving frame, the energymomentum tensor takes the form where ǫ and p are the matter energy-density and thermodynamic pressure, respectively, and U µ is the four-velocity, satisfying the normalization condition For the perfect fluid we can fix the matter Lagrangian L m as L m = p, which gives The f (R, T ) gravity theory is a non-conservative theory, and for the energy balance equation we obtain (see Appendix B for the details of the calculations) B. The effective cosmological Lagrangian and the potential in f (R, T ) gravity theory We assume that the geometry of space-time is described by the Friedmann-Robertson-Walker (FRW) metric, which in spherical coordinates is given by where N (t) is the lapse function, a(t) is the cosmological scale factor, and k = 1, 0, −1, correspond to the closed, flat and open Universe models, respectively. To proceed further, as a first step we obtain the effective Lagrangian for the f (R, T ) theory, whose variation with respect to its dynamical variables yields the appropriate equations of motion.
The trace of the energy-momentum tensor is T = −ǫ + 3p. In the comoving reference frame the components of the four-velocity are U µ = (N (t), 0, 0, 0) and U µ = (−1/N (t), 0, 0, 0), respectively. Therefore the trace of the field equation Eq. (5) of the f (R, T ) gravity model becomes With the use of the above identity, and by taking into account the explicit form of the components of the Ricci tensor, given in Appendix C, we obtain the cosmological action for the f (R, T ) gravity theory as In Eq. (11) λ and µ are Lagrange multipliers. The term with the second Lagrange multiplier is chosen as the contracted field equation, because it is derived directly from the action, and no further assumptions are needed to be introduced. Moreover, if we use other formulations for the action, like adopting for the factor multiplying the second Lagrange multiplier the form T + ǫ − 3p, we will lose important insights from the geometry of the modified gravity, and we will have to face conceptual problems when ǫ and p are related by the radiation equation of state p = ǫ/3. After taking the variation of Eq. (11) with respect to R and T , we obtain the expressions of λ and µ as follows Hence we obtain the gravitational part of the Lagrangian as where the potential V reads In order to make the presentation simpler, we introduce the following notations Now λ can be expressed using the new variables as In the following we further denote Thus we have After taking the time derivative, we finḋ C. The cosmological Hamiltonian in f (R, T ) gravity theory The canonical momentum associated to the coordinate q is given by P q = ∂L ∂q . Hence the cosmological Hamiltonian of the f (R, T ) gravity is given by After some simple calculations we obtain the explicit forms of the canonical momenta in the cosmological f (R, T ), which are presented in Appendix D. With these canonical momenta, we obtain the cosmological Hamiltonian of the f (R, T ) gravity theory, In order to simplify the notation we will represent the Hamiltonian as H grav = (· · ·) − 6 N a 2ȧ [· · ·]. Now we can easily find the relation Since by combining the above two equations, we obtain (26) Therefore for the gravitational Hamiltonian we find Since we have P A = − 6 N a 2ȧ 1 − B A , after taking the square of it, we have Thus we arrive at the following form of the Hamiltonian, Since we have Therefore we obtain the gravitational part of the cosmological Hamiltonian of the f (R, T ) gravity theory as where the quantum potential V is given by and thus the total Hamiltonian of the system is The gravitational Hamiltonian constructed above consists of all canonical momenta associated to all variable of the f (R, T ) gravity theory, which can lead to the existence of a complex dynamics of these field variables, and of their associated canonical momenta.
Here the Hamiltonian operator H takes the form In order to quantize the model we perform first pa-rameter ordering. There are various ways to do it [57], but in the following we choose a procedure that keeps the Hamiltonian Hermitian [41]. Hence we assume the following relations (we let P q = −i(∂/∂q)), where the parameters u, v, w satisfy u + v + w = 1, and denote the ambiguity in the ordering of the factors q and P q , where the parameters r, s satisfy r+s = 1, and denote the ambiguity in the ordering of factors q and P q . Similarly we find Thus we obtain the quantized cosmological Hamiltonian in f (R, T ) gravity theory as Here u 1 , w 1 , u 2 , w 2 and u 3 , w 3 denote the ambiguity in the ordering of factors A, P A , F, P F , and G, P G , respectively.
In the next Section we will investigate some particular quantum cosmological models in the f (R, T ) gravity theory.
After considering the general case of the Wheeler-de Witt equation in the previous Section, we can see that an analytic general solution of this equation for arbitrary f (R, T ) would be difficult to obtain. Instead, in the present Section we consider a specific case, in which the gravitational action is of the simple form where F 0 (R) is an arbitrary function of the Ricci scalar, and θ is an arbitrary function depending on the scale factor a(t) only. In this toy model, the coupling of the curvature of space-time and the trace of the matter energymomentum tensor can give a hint of the implications on the quantum cosmological evolution of the existence of such a coupling.

A. The Hamiltonian and the Wheeler-de Witt equation
In the newly born quantum Universe, the space-time has a very high curvature, so that R → ∞. Accordingly, in the present model we have Similarly, the other variables become When R → ∞, B/A = 1/2. From the definition of P A and P F , we obtain By assuming that in the new born quantum Universe In the above equation the quantum potential is defined as In the limit R → ∞, and by assuming that θ = 0, the potential becomes In this approximation the total Hamiltonian of the system becomes If in the Lagrangian/ Hamiltonian we have some terms that can be omitted in the action, they can also be omitted in L and H, without causing any physical differences in the dynamics of the system. In Eq. (52, due to Gauss theorem, we have Then the variational derivative of this term vanishes, Therefore this term can be omitted in the Hamiltonian function. Thus we arrive to the Wheeler-de Witt equation for the f (R, T ) = F 0 (R) + θRT gravitational model, which has the form If the physical variables do not depend explicitly on the time t, we obtain Therefore we can formulate the classical equations of motion of the f (R, T ) gravity theory aṡ Now we define a new time variable τ , which has the following relation with the original time variable t, or, equivalently, dτ /dt = N (t). From the definition of P M , we have where h ≡ 1 Nȧ a is the Hubble function. Then the Hamilton equations of motion for a cosmological fluid become where a prime represents the derivative with respect to τ . The first equation of the above system gives us the coupling constant θ = D as This result tells us that the coupling between the gravitational field and the matter field decreases as the scale of the Universe increases. This result may be the reason why the coupling between matter and gravity becomes so weak in the limit of large cosmological times, and the f (R, T ) gravity behaves as the standard gravity nowadays.

C.
The time problem in the quantum cosmology of the f (R, T ) gravity The absence of the time evolution of the wave function of the Universe in the Wheeler-de Witt equation hinders our efforts of understanding quantum gravity in a way similar to standard quantum mechanics, or quantum field theory. Below we propose a way to turn the Wheeler-de Witt equation into a Schrödinger type equation. For the product P D P T we obtain immediately the result In the following we make the fundamental assumption that the above term can be interpreted as P τ , which is the canonical momentum for time. There is a convincing reason for this assumption: after performing quantization P q = −i ∂ ∂q , we obtain This result holds if we assume that a(t) is the only time dependent variable in the model. The relation proved above allows us to perform the transformation The validity of this transformation shows us that the coupling between the gravitational field and the matter field may play an important role in the evolution of the very early Universe. In the following we will make a further simplification of the gravitational action, by assuming it to be Then C = F = P R = 0, and for the cosmological Hamiltonian of the system we obtain This Hamiltonian is very similar to the Hamiltonian obtained in the f (R) model presented in [41], except for the appearance of a new term −2 N 6a 3 P D P T + N ( a 3 4 + 3ka)M , which shows us the effect of the coupling between space-time and matter. The Wheeler-de Witt equation HΨ = 0 for this f (R, T ) gravity model reads After the use of the transformation introduced in Eq. (79), we obtain After substituting all P q = −i ∂ ∂q to quantize the model, we obtain the corresponding Schrödinger-Wheeler-de Witt (SWDW) equation describing the quantum evolution of the Universe as which is just of the form of the standard Schrödinger equation, Therefore in the f (R, T ) gravity theory, we can generate a Schrödinger type equation from the Wheeler-de Witt equation, which can solve the time problem in quantum gravity. When dτ dt = N (t) = 1, the WDW equation will take the form of the Schrödinger equation we are familiar with, Let us take now a deeper look into the time problem of quantum gravity, and analyze the physical meaning of the time τ , and of the effective Hamiltonian H ef f we have introduced here. In the Wheeler-de Witt equation Hψ = 0 there seems to be no dynamics of the system. Therefore it turns out that the wave function of the Universe (more precisely, the physical states ) do not describe states of quantum gravity at a particular time, as in the standard quantum theory. Rather, they describe states for all times, or, more precisely, just that information about the state of the Universe that is invariant under all space-time diffeomorphisms [58] In the modified gravity model f (R, T ) = R + θRT , time may be introduced locally by the coupling of the gravitational field and the matter field, because the interaction between the gravitational field and the matter field is also local. The profound connections between thermodynamics and gravity tell us that the arrow of time may come from the second law of thermodynamics, since both processes reveal irreversible dynamics. If we only consider a small pitch of the Universe, and think of it as an adiabatic system, with HΨ = 0 still valid in it, the coupling between curvature and matter will generate an arrow of time to measure the increase of its entropy, contributed by matter. The other components of the Hamiltonian, included in H ef f , show us the dynamics of the gravitational field, and of some matter components, so that the effective Hamiltonian takes the form H ef f ∼ b 1 p i p j + b 2 x i x j , which is similar to the Hamiltonian we usually meet in quantum mechanics. Therefore one can suppose that the Schrödinger equation in ordinary quantum mechanics might describe just a locally effective theory of the Wheeler-de Witt equation.
In other words, we may conjecture that the Wheeler-de Witt equation provides the global quantum description for the Universe, while the Schrödinger equation is just the local description for the present day microscopic regions of the Universe.
Here E = constant. Thus we obtain the following differential equation describing the time evolution of the wave function of the Universe, By introducing the new variables x = aA 1 2 and y = A, we obtain the differential equation for ψ In the following we approximate the equation of state of the early Universe by the stiff equation of state p = ǫ, since in the very high density Universe one expects the speed of sound c s to be of the same order of magnitude as the speed of light, c s = ∂ǫ/∂p = 1. When p = ǫ, and sinceṘ = − N 3a 3 P C = 2ȧ a 1 2 R =ȧ a R, from Eq. (8) it follows that in this case the energy is conserved, ǫ ′ +3(ǫ+p)h = 0. Hence we have where w > 0 is a positive constant. By substituting this result into Eq. (89), and by letting D = δ/a(τ ), we find

A. Time evolution as a non-constant energy perturbation
In order to obtain a solution of the SWDW equation obtained above, we can estimate different terms' order of magnitude. We know that x = aA 1/2 ∼ 1/a 5 2 , y = A ∼ 1/a 7 , and hence we can obtain (92) Since a is very small, the term x 3 /y 1 2 ∼ 1/a 4 can be thought as a small perturbation. To obtain some analytic solutions, we consider the case when k = 0 (flat Universe). This choice can help us to get rid of the large coupling term y 5 2 /x 3 ∼ 1/a 8 . Due to the discussion by Hawking and Page [59], we also assume that u 1 w 1 , the ordering parameter, can be neglected. In our case this can be achieved by putting it into the energy term, or neglecting it directly, since it is small as compared to the variables related to a(τ ).
After making the above assumptions, we obtain the Schrödinger-Wheeler-de Witt equation as (93) where Ψ = Ψ(x, y, τ ). In the stationary situation, we can decompose the variables as Ψ(x, y, τ ) = e −iEτ ψ(x, y), and thus we get (94) Now we change the above equation into another form, and multiply by 3x 3 y 1 2 both sides. We thus have the equa- (95) In the first formulation of the SWDW equation, given by Eq. (94), we can think of the time evolution of the wave function as resulting in the addition of a constant E as a perturbation. But what we know from the perturbation theory of quantum mechanics tells us that a constant perturbation gives us no changes in the energy level and in the wave function, which is to say that it does not affect the physical observables. After changing the unperturbed Hamiltonian in the following way we obtain the new version of the SWDW equation, where we do not consider the constant E as induced by a time evolution effect anymore. Instead, we consider it as a non-constant perturbation of the system, such that the total Hamiltonian is H = H 0 + V pert . Hence we can say that the time evolution effect on the wave function in a non-perturbative system, like the f (R, T ) gravity theory, is equivalent to the splitting of the degenerate energy levels (loss of symmetries) in a perturbative, static system. This result shows that the time evolution in the quantum cosmology of f (R, T ) gravity leads to the splitting of the degenerate energy levels, which reveals the deep connections between energy and time.
In the new form of the Hamiltonian the unperturbed component H 0 satisfies the equation H 0 ψ = 0, whose eigenvalue of energy is zero. In the unperturbed Hamiltonian we can separate the variables as ψ(x, y) = X(x)Y (y), and obtain where the separation constant is denoted as v 2 −1 4 . Then we obtain the expression of the wave function for the unperturbed Hamiltonian H 0 as where is the Bessel function, and J v (x) and J −v (x) are linearly independent functions. Since in the early Universe, we have x → ∞, in order to have analytic solution in the whole plane, we let A 1 = B 2 = 0, and we assume v ≥ 0.

B. The twofold degenerate case
In order to investigate the energy level split, we consider the simplest case, the twofold degenerate case [60]. Assuming that the system is twofold degenerate, at the beginning of time τ = 0, the wave function can be written as where c 1 , c 2 are constants satisfying the relation |c 1 | 2 + |c 2 | 2 = 1, and where v 1 and v 2 are positive constants. For the perturbed system with perturbation V pert = −3E x 3 y 1 2 , we write With the use of the unperturbed wave functions, we obtain [J v2 ( 3wy)] 2 dxdy, (105) Since the curvature scalar and the scale factor are always positive, we cannot define wave functions that are analytic on whole space, thus the orthogonality and normalization are not satisfied. In the following we define the quantities In the special function theory, we have the Schafheitlin integral [61], which reads where a and b are positive constants, and satisfy the relations below to make the integral convergent, For z ∈ ℜ the Gamma function Γ(z) satisfies the identities Γ(z + 1) = zΓ(z), Γ(1) = 1, Γ( 1 2 ) = √ π. In the following we will restrict our analysis to the simplest case v 1 = 3 2 , v 2 = 5 2 . With the use of these values we obtain where we have denoted z = √ 3wy → dz = √ 3w 2 √ y dy, and η = −6E √ 3w, respectively, and we have also assumed that x and y are independent variables. The upper limits in the integrals are both ∞ when a(τ ) → 0. As for the lower limits of integration we assumed them to be very small positive numbers l and l 1 , which correspond to the transition from the quantum regime to the classical regime, and thus they describe the limit of applicability of the present quantum model of the Universe. From the numerical evaluation of the integrals it follows that when l 1 < 0.5 the integral values will not vary much, and therefore we let l 1 = 0.

C. The energy levels of the quantum Universe
We use the degenerate perturbation theory to find the energy levels in the quantum cosmology of f (R, T ) gravity. The wave function is given by Ψ = c 1 ψ + c 2 ψ 2 , and we already know the eigenvalues and the eigenfunctions of the unperturbed Hamiltonian H 0 , which are given as solutions of the equation By substituting the above equation in the Schrödinger equation, HΨ = (H 0 + V )Ψ = EΨ, we obtain After multiplying Eq. (119) with ψ * 1 and ψ * 2 , and integrating over a volume V , we obtain the equations In the degenerate situation, we let E = E (0) + E (1) , and c n = c (0) n , that is, we take for these coefficients the zero order (unperturbed) approximation. Then we have From these two equations we can find the coefficients c 1 , c 2 , after the perturbed energy E (1) is obtained. The secular equation is Then we obtain the first order modifications of the energy as where In the case analyzed earlier, where v 1 = 3 2 , v 2 = 5 2 , we obtain And the two modified energy levels are where β ± are two constants. From the above equation we know that the modified energies have the same sign, and they are proportional to E, as well as proportional to l 3 , which will tell us the lower limit of the size of the energy gap, if we know the numerical value of l. From Eq (122) the wave function's coefficients are obtained as equation and respectively. The above expressions show the dependence of the coefficients c 1 and c 2 on l. When l is large, the Universe will have a higher probability to be in the state ψ 2 , and it will have a higher probability to be in the state ψ 1 when l is small. Note that the coefficients c − 1 and c 2 do not depend on η, which implies that the parameter E will not affect the state of the wave function.

D. The transition probability in the quantum Universe
In the previous Section we have considered the early Universe as a quantum system that has twofold degenerate energy levels corresponding to the wave functions ψ 1 , ψ 2 at the time τ = 0. In the following we will consider the probability of transition of the Universe, from the ψ 1 state at τ = 0, to the ψ 2 state at time τ , the transition taking place due to a perturbation of the initial state.
In the zero order approximation the wave functions are where c 1 , c 2 and c ′ 1 , c ′ 2 are the two pair of coefficients obtained previously. Here ψ and ψ ′ are the wave functions corresponding to two energy states E 0 + E After reintroducing the time factor, we obtain the time dependent wave functions as Note that Ψ 1 = ψ 1 at τ = 0. Then we use ψ 1 , ψ 2 to represent ψ, ψ ′ , and hence Ψ 1 becomes the linear combination of the wave functions ψ 1 , ψ 2 , with the combination coefficients time dependent. The absolute value of the coefficient multiplying ψ 2 and integrating over the volume V , is (after squaring) the transition probability w 21 . Therefore we have With the use of the previous results, the transition probability in ordinary units is (139) From the above equation it follows that the transition probability is a cosine function of the time τ . When τ is very small we have w 21 proportional to τ 2 . Since we know that the probability is smaller than one, we have the restriction 0.0106/l 10 − 0.0066/l 10 ≤ 1, a condition which gives for the upper limit of l the numerical value l ≥ 0.666 → a(τ ) ≤ 1.177.
(140) Therefore we get the upper limit of the size of the Universe for which our quantum model is applicable. And the lower limit of y is 0.320 < 0.5, and therefore we can approximate the lower limit in the integral as l 1 = 0.

V. THE SECOND QUANTIZATION OF TIME
In standard quantum mechanics, time is not an operator. The energy of the system is the eigenvalue of the Hamiltonian, which in the case of the harmonic oscillator can be written in the formĤ = ω(a + a+ 1 2 ), where a and a + are the creation and annihilation operators. On the other hand the deep connections existing between energy and time suggests us to find a way to define the creation and annihilation operators of 'time'. After these operators are found, we can get rid of the concept of time singularity at the beginning of the Universe, and we can properly define the distance between each pair of time slices. And, consequently, we can obtain the quantum frequency of the 3-space evolution, and therefore investigate from a quantum mechanical point of view the birth of the Universe.
In the previous Section, by using the mathematical formalism of the f (R, T ) gravity theory we have defined a 'time' variable in the WDW equation, and thus we have transformed it into a Schrödinger-Wheeler-de Witt equation. The time τ we have defined earlier is based on the relation P τ = − 2 6a 3 P D P M , and therefore in our analysis we have assumed that the idea of time in the f (R, T ) gravity theory is related to the field variables D and M . In the following we want to define the creation/annihilation operators based on the term 3kaDM − 2 6a 3 P D P M in the WDW equation (84). The procedure goes as follows.
We assume the existence of the classical and quantum Since (144) In the following we denotê and respectively. Also v * ,P v * = (2i ) * .
Let's assume now that there are two directions of time, defined asτ Then we have (152) As a next step in our analysis we interpretτ 1 as a creation operatorτ + andτ 2 as an annihilate operatorτ . Thus we further obtain Here we define the complex time "number" operator N τ ≡τ +τ (154) The above relation shows that the time "number" operator, defined asN is an observable in quantum mechanics. Although the complex time has two directions, there is only one real time observable that can be measured in the experiments. Now let's consider the case of our WDW equation, where (156) In the following we discuss again the specific cosmo- Thus we obtain the corresponding representation of the time-related terms in the WDW equation as From this relation it follows that the coupling between the gravitational field and the matter field gives us a way to measure the quantum time number of a given Universe. With the use of the above relations the WDW equation becomes (159) Eq. (159) gives us the possibility of further constructing a wave function with N quanta of time, which has the propertyN Therefore the WDW equation becomes (161) Eq. (161) may give us a clear picture of the evolution of the Universe when we use the Arnowitz-Deser-Misner (ADM) formalism, since the dynamics of system will change with the variation of the time quanta N . Therefore when we are at different time moments, we shall have different Schrödinger equations to describe the local dynamics of the Universe.
Of course, we can also define the observable time vacuum (the beginning of the Universe) as In the present Section we have tried to introduce the Fock space of the quantum time variable τ . If we define the time creation/annihilation operators in the Fock space, then we can get rid of the problem of time singularity at the beginning of the Universe, and we may have a deeper understanding of the discrete nature of time. But this will also require to transform our wave function into the quantum occupation number picture. However, we must note that this technique works only when the canonical momentum and its corresponding canonical position are not coupling with each other. One of the difficulties of the canonical quantization of gravity is that it cannot define the Hilbert space. The second quantization of time procedure may provide a prospective way to think about the problems of quantum gravity in the Fock space, instead of in the Hilbert space.

VI. DISCUSSIONS AND FINAL REMARKS
Quantum cosmology offers a lange number of challenges, but also interesting insights into the fundamental nature of the space-time. The physical problems related to the birth and very early evolution of the Universe might be better understood by using the mathematical formalism of quantum theory, including symmetries, discrete structures, or semi-classical features extracted from a generally covariant, and highly interacting quantum theory.
In this paper we have investigated the quantum cosmology of the f (R, T ) gravity theory, a modified gravity theory in which the gravitational action is an arbitrary function of the Ricci scalar and of the trace of the energymomentum tensor. In the present paper we have considered that the classical evolution of the Universe takes place in the background Friedmann-Robertson-Walker geometry, which we are using systematically to investigate the quantum properties of the early Universe. As a starting point in our analysis we have introduced the Hamiltonian formulation of the theory, which is constructed systematically from the action given by Eq. (11). In the action we have introduced two Lagrange multipliers λ and µ, with the first imposing the (purely geometric) definition of the Ricci scalar, while the second one imposes the trace constraint of the f (R, T ) gravity theory. This constraint goes beyond a simple definition of the trace of the energy-momentum tensor, since it allows the investigation of the deep connection between matter and geometry at a more general level than the one that follows from the simple thermodynamic definition of T . From the cosmological gravitational action one can obtain the gravitational cosmological Hamiltonian, which, by canonical quantization, leads immediately to the general form of the Wheeler-de Witt equation, describing the evolution of the wave function of the quantum Universe. In order to obtain some physical insights in the quantum properties of the Universe we consider a simple extension of the standard general relativity, in which the gravitational Lagrangian is a "deformation" of the form θRT of the general relativistic Lagrange function R. We have investigated in detail the properties of this quantum cosmological model. Its most interesting feature is the possibility of the definition of a quantum time, and of an associated canonical momentum operator. This leads to the reformulation of the Wheeler-de Witt equation as a Schrödinger type equation. We have studied in detail the mathematical properties of this equation, by using a perturbative approach, in which the small perturbation is proportional to the energy of the system, in the framework of a twofold degenerate quantum system. The probability of transition between states is also obtained. As a theoretical possibility we have also discussed very briefly the second quantization of time, which leads to the interesting possibility of the extension of the Hilbert space of the canonical quantization method to the Fock space description of quantum phenomena in the very early Universe.
The initial state of the Universe is essentially unknown. That's why the possibility that the initial geometry of the Universe was not an isotropic and homogeneous, Friedmann-Robertson-Walker type one, cannot be rejected a priori. This raises the interesting question of the applicability of the formalism developed in the present to describe the quantum cosmology of f (R, T ) gravity to more general geometries. In particular, in the following we briefly consider the quantum cosmology of f (R, T ) gravity in the anisotropic Bianchi type I geometry, with the metric given by ds 2 = −N 2 (t)dt 2 +a 2 1 (t)dx 2 +a 2 2 (t)dy 2 +a 2 3 (t)dz 2 , (163) where a i , i = 1, 2, 3 are the directional scale factors. For the Bianchi type I geometry the scalar curvature is obtained as For the Bianchi type I geometry the gravitational action reads while the gravitational Lagrangian can be obtained as L g = − 2 λ N (ȧ 1ȧ2 a 3 +ȧ 1ȧ3 a 2 +ȧ 2ȧ3 a 1 ) − 2 N˙ λ d dt a 1 a 2 a 3 ) + a 1 a 3 a 2 + a 1 a 2 a 3 − N a 1 a 2 a 3 V.
Hence the gravitational Lagrangian of the f (R, T ) gravity theory in a Bianchi type I geometry can be represented as The Hamiltonian corresponding to the Lagrangian (168) now reads H g =Ẇ P W +ȧ 1 P a1 +ȧ 2 P a2 +ȧ 3 P a3 + · · · − L g , (169) and we obtain the canonical momenta corresponding to the variables (W, a 1 , a 2 , a 3 ) as For the other canonical momenta we obtain the following correspondence between the isotropic and anisotropic case, Hence the gravitational Hamiltonian becomes where [· · ·] is represented by (173) Since for the Hamiltonian we obtain The anisotropic cosmological Hamiltonian (175) of the f (R, T ) gravity theory is very similar to the one obtained in the isotropic case.
When considering the f (R, T ) = F 0 (R) + θRT model in the case of R → ∞, we obtaiṅ and hence we immediately arrive atḊ = −(3/2) Ẇ /W D, and D = δ 0 /W 3/2 , respectively, where δ 0 is a constant. Since then we have By taking into account that N 2W 3 λ a 2 i P 2 ai =ȧ i P ai = −ia i ∂ ∂a i , i = 1, 2, 3, we obtain the time canonical momentum Finally, the transformation will allow us the introduce the time dependence of the Wheeler-de Witt equation for anisotropic Bianchi type I geometries in the f (R, T ) gravity theory. Note that here we have assumed that a 1 (t), a 2 (t), a 3 (t), W (t) are independent variables. Hence we can safely conjecture that in the anisotropic case we can still introduce a cosmological quantum time in the Wheeler-de Witt equation of f (R, T ). On a qualitative level the overall results of the anisotropic case will be very similar to the ones obtained for isotropic and homogeneous geometries.
In the present paper we have introduced some basic theoretical tools that could be used for the investigation of the quantum properties of the gravitational interaction, and of the evolution and origin of the very early Universe, in which the complex interaction of geometry and matter give birth to time, entropy, and irreversibility.
By simple calculations we obtain the components of the Ricci tensor for the FRW metric (9) as The canonical momenta associated to the cosmological action Eq. (11) of the f (R, T ) gravity are given by