Anisotropic n-dimensional quantum cosmological model with fluid

In the present work, we discuss the Wheeler-DeWitt quantization scheme for an n-dimensional anisotropic cosmological model with a perfect fluid in presence of a massless scalar field. We identify the time parameter using the generalization of Schutz formalism and find the wavepacket of the universe following standard Wheeler-DeWitt methodology of quantum cosmology. We also calculate the expectation values of scale factors and volume element. We show that quantized model escapes the singularity and supports bouncing universe solutions.


I. INTRODUCTION
Higher dimensional models had been investigated widely in the past in order to find a theory to unify gravity with other fundamental forces of physics. It started with Kaluza and Klein's assertion that a fifth dimensions in general relativity will unify gravity with the electromagnetic field [1][2][3]. Further motivation of resrting to higher dimensional models came from the expectation of unifying gravity with non-Abelian gauge fields [4,5]. Later, spacetime having more than four dimesnions were motivated by 10-dimensional superstring theory and 11dimensional supergravity theory [6,7].
Quantum cosmology has provided a way to study the early phase of the universe. In the absence of a well accepted quantum theory of gravity, it is quite imperative to look at the quantum behaviour in individual gravitational systems. Wheeler-DeWitt equation is a widely used framework where standard quantum mechanics is used in cosmological models [8][9][10]. This framework of quantization suffers from various conceptual problems, which are discussed quite elaborately in the literature [11][12][13]. The problem of non-unitary evolution and the absence of properly oriented time are amongst them. The latter is dealt with in various ways [14][15][16][17][18][19]. One of the ways of finding a nice time parameter is to use the evolution of the cosmic fluid following Schutz formalism [20,21]. This method has been widely used as time parameter in quantization of cosmological models [22][23][24][25]. In fact, the alleged loss of unitarity in anisotropic quantum cosmological models has been quite effectively resolved [26][27][28][29][30] using Schutz formalism. In cosmological models related to Brans-Dicke theory [31,32], the evolution of the scalar field have been used as time parameter in order to quantize the same without adding any additional matter [33][34][35]. The Brans-Dicke theory of gravity has also been quantized using Schutzs formalism [36,37]. So long as a method gives rise to an oriented scalar time parameter, it is quite an effective method. * sp13ip016@iiserkol.ac.in In a very recent work, Alves-Junior et al [38] have shown the canonical quantization of n-dimensional anisotropic model coupled with a massless scalar field in the absence of any fluid where they have emphasized on natural identification of scalar field as time parameter for the evolution of quantum variables.
In the present work, we try to investigate the quantum cosmological solutions of an n-dimensional anisotropic model in which massless scalar field is minimally coupled with gravity in presence of a barotropic perfect fluid(P = αρ). The idea is to generalize the work of Alves-Junior et al to include a fluid, as well as to generalize the recent work on anisotropic cosmologies by Pal and Banerjee [26,27] to higher dimensions. Keeping in mind the importance of various aspects of higher dimensions, it is quite relevant to ask the questions like whether singularity-free n-dimensional quantum cosmological models realized. Also, this investigation should find if the general result, given by Pal and Banerjee [30], that it is always possible to find a self-adjoint extension of homogeneous quantum cosmologies in four dimensions.
In section II, we write Lagrangian for the gravity sector coupled with a scalar field, followed by the mention of Schutz formalism to n-dimensional models. Similar kind of generalization of Schutz fluid was done by Letelier and Pitelli [39] in the context of the quantization of an n-dimensional FLRW model. We construct the Hamiltonian using canonical transformation in subsection B of this section. In section III, canonical quantization of the model is presented for a stiff fluid(α = 1). We obtain the wave packet and expectation values of scale factor and volume element. Further, we are also able to show the unitary evolution of the model with general fluid(α = 1). This is done for an example with n = 7 which is mathematically tractable. In the end, in section IV, we discuss our results and make some remarks.

II. AN N-DIMENSIONAL COSMOLOGICAL MODEL
A. Lagrangian of the model We start with the Hilbert-Einstein action in ndimension with a minimally coupled scalar field in the presence of a fluid, where ω is a dimensionless parameter and φ is the massless scalar field. The last term of equation (1) represents the matter contribution coming from fluid pressure P related with density ρ by the equation of state P = αρ.
The metric for an n dimensional (n > 4) spatially homogeneous but anisotropic cosmology is chosen as where N (t) is the lapse function, a(t) is our good old usual scale factor used in flat FLRW cosmology and b(t) is scale factor coming from remaining (n−4) dimensions. So the usual 3-space is isotropic, but the extra dimensions, though isotropic in itself, but anisotropic from the usual 3-space section. The action for scalar gravity part can be written in a reduced form for metric given by (2) as where an overdot represents a derivative with respect to coordinate t and V 0 denotes (n − 1) dimensional volume. We have ignored the surface terms as they do not contribute to the field equations.
From (3), we can write the Lagrangian for the gravity sector as, Schutz [20,21] showed that four dimensional velocity vector of a perfect fluid can be written using six thermodynamic quantities as, where µ is specific enthalpy, S is specific entropy, ζ and β are potentials connected with rotation and and θ are potentials with no clear physical meaning. For n dimensional velcoity vector, we can consider additional potential terms [39] like ζ and β in equaution (5) as Since the potentials ζ, ζ i , β andβ i are related to rotation, we can neglect them in the present case as there is no vorticity in the spacetime given by equation (2). Following the widely used technique, developed by Lapchinskii and Rubakov [40] by writing the fluid pressure in terms of potentials given in equation (6), the fluid sector action given in equation (1) can be written as From above equation, Lagrangian for fluid sector can be written as . Hamiltonian of the model Following standard procedure [41], Hamiltonian for gravity sector from corresponding Lagrangian (4) can be written as For fluid sector, Hamiltonian corresponding to Lagrangian (8) can be written as where T and p T are related with following canonical transformation (11) Using (9) and (10), the net Hamiltonian can be written as The above expression for H may look formidable, but can be written in a simpler form with the help of following canonical transformations, as For n = 5, the first term of equation (15) where u and v are given as,

)B. (18)
A. For stiff fluid α = 1 and n > 5 For α = 1, equation (16) can be rewritten as 2 − 6(n − 4) (n − 5) Solution for equation (19) can be written as wheren = 6(n − 4)/(n − 5) and λ, k and K are constant. Here we take E > k + λ. By superposing the function ψ λ,E,k (u, v, φ, T ), the wavepacket can be formed as; where E = E − k − λ and γ is a positive constant. Norm Ψ can be calculated as, This clearly shows that we have a finite and timeindependent norm. From equation (21), the normalized wavepacket can be written as, Lengthy expressions of the above integrals have not been mentioned here. Therefore, plots for the same are shown in the figures taking the value of positive constant γ = 1.  For n = 5 with α = 1,equation(15) will become 1 6 With the change of variable as x = 2A + B and y = B, it can be re-written as, Solution of above equation can be found using similar method what was done in the previous section. Normalized wavepacket is given as,  We will be considering only n = 7 as an example for a general fluid with α = 1, in which case the equation is tractable. Now Wheeler-DeWitt equation (16) can be written as, With suitable operator ordering as done in work by Pal and Banerjee [26,27], equation (32) takes following form, , above equation can be re-written as, With change of variable like χ = e (3α−3) 2 √ 30 u , this equation can be written as, Above equation is similar to the inverse square potential problem of physics which admits the self-adjoint extensions. The solution to equation (35) can be given by Hankel functions. It has been extensively shown in the work [26].

IV. DISCUSSION WITH CONCLUDING REMARKS
In this work, using the n-dimensional generalization of Schutz formalism [39], we formulated the Lagrangian and Hamiltonian for our model. Followed by multiple canonical transformations, we obtained the complicated looking Wheeler-DeWitt equation (16) for the wavefunction of the n-dimensional anisotropic universe. We obtained the general wavepacket of Wheeler-DeWitt equation for stiff fluid α = 1. We are able to find a time independent and finite normed wavepacket which establishes the unitary evolution of the model.
We have found the non-singular expectation values of the scale factors and these are shown in Fig. 1, 2, 3 and 5 for different n. Expectation values of scale factors a and b coincide for n = 7. This is not surprising as both usual space-section and the extra 3-dimensional space, although anisotropic between each other to start with, are 3-dimensional isotropic space like sections in themselves.
Non-zero minima of the expectation values of volume element (Fig. 4) clearly show the contraction of the universe followed by expansion avoiding any singularity indicating a bouncing universe. Similar results were obtained in quantization of anisotropic models in absence of fluid in the past [38].
For a general fluid with α = 1, we are able to show the self-adjoint extension and thus unitary evolution of the model for n = 7 despite the mathematical difficulty of the model. However, the result of it is in accordance with four-dimensional anisotropic model Bianchi-I as expected [26].
Our results may also be considered as the higher dimensional generalization of unitary evolution of the anisotropic models as shown in [26][27][28][29].