This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Phase space deformations in phantom cosmology.

J. L. Lópeza jlopez87@uabc.edu.mx M. Sabidob,c msabido@fisica.ugto.mx C. Yee-Romeroa carlos.yee@uabc.edu.mx a Departamento de Matemáticas, Facultad de Ciencias
Universidad Autónoma de Baja California, Ensenada, Baja California, México.
bDepartment of Theoretical Physics, University of the Basque Country UPV/EHU, P.O. Box 644, 48080 Bilbao, Spain.
c Departamento de Física de la Universidad de Guanajuato,
A.P. E-143, C.P. 37150, León, Guanajuato, México.
Abstract

We discuss the physical consequences of general phase space deformations on the minisuperspace of phantom cosmology. Based on the principle of physically equivalent descriptions in the deformed theory, we investigate for what values of the deformation parameters the arising descriptions are physically equivalent. We also construct and solve the quantum model and derive the semiclassical dynamics.

keywords:
Noncommutative cosmology, deformed phase space models, late time acceleration.
PACS:
04.20.Fy, 04.50.Cd, 98.80.-k,98.80.Qc

1 Introduction

One of the most intriguing aspects in physics, is the current acceleration of the Universe. Is it a consequence of modifications to general relativity (GR) or is it a new kind of matter that drives this acceleration?. Although having some theoretical problems, the best answer to this question is the cosmological constant Λ\Lambda. An alternative candidate that has been successful for the description of dark energy is the scalar field [1, 2, 3, 4, 5, 6]. In particular, scalar fields with a negative kinetic term have been considered in the literature. This type of fields are known as Phantom fields. They have some interesting properties that allows them to be considered as a strange but viable matter which could be relevant in the evolution of the Universe [7]. In particular, a phantom field provides an effective negative pressure and a repulsive effect on the matter content of the Universe which in the long term could be responsible for the current accelerated expansion [8, 9, 10, 11, 12]. Therefore it has been considered as the matter source of the late time accelerated expansion of the Universe [13, 14, 15].

When considering the spacetime structure of the Universe it is usually done in reference to GR. But when regarding the micro structure of spacetime we do not have a universally accepted quantum theory of gravity. Although there are several candidates, noncommutativity has been considered as an alternative to understand the small scale structure of the Universe and help in the construction of quantum theory of gravity. For this reason, noncommutative versions of gravity have been constructed [16, 17, 18, 19]. Noncommutativity is usually believed to be present near Planck’s scale and is consistent with a discrete nature of spacetime. Motivated by this idea, it is justified to consider an inherently noncommutative spacetime at the early ages of the universe. Directly using noncommutative gravity is quite difficult, this is a consequence of the highly nonlinear character of these theories[16, 17, 18, 19]. Fortunately there is an alternative, in [20] the authors introduce the effects of noncommutativity using the Moyal product of functions on the Wheeler-DeWitt (WDW) equation.
The effects of the noncommutative deformation at the classical level was studied by a WKB approximation of noncommutative quantum model [21] and also by modifying the Poisson algebra [22, 23]. More general minisuperspace deformations have been done in connection with Λ\Lambda [24, 25, 26, 27, 28, 29]. Phase space deformations give rise to two physically nonequivalent descriptions, the “C-frame” based on the original variables but with a modified interaction and the “NC-frame” frame based on the deformed variables but with the original interaction. Given this ambiguity, in [29] a principle was proposed to restrict the value of the deformation parameters in order to make both descriptions physically equivalent.

In this work we consider a Friedmann-Robertson-Walker (FRW) cosmological model coupled to a phantom scalar field and study the physical consequences of introducing general phase space deformations on the minisuperspace of the theory. We will study both, the classical and quantum models and we also find that the semiclassical approximation of the deformed quantum model agrees with the classical model. The deformation parameter space is determined by considering the principle of physically equivalent frames [29].

The paper is organized as follows. In section 2, the commutative model is presented. In section 3, the minisuperspace phase space deformation is implemented and the dynamics in the two frames is obtained. Also the deformation parameters are constrained by imposing physical equivalence between the “C-frame” and the “NC-frame”. The quantum analysis is done in section 4, we find the solution to the deformed WDW equation and fix the parameters in the deformation in order to make the quantization possible. We also show that the classical paths arise from the semiclassical approximation of the WDW equation obtained from the deformed Hamiltonian. Finally, section 5 is devoted for concluding remarks.

2 Phantom field model

We start with the flat FRW metric

ds2=N2(t)dt2+a2(t)[dr2+r2dΩ]ds^{2}=-N^{2}(t)dt^{2}+a^{2}(t)[dr^{2}+r^{2}d\Omega] (1)

a(t)a(t) corresponds to the scale factor and N(t)N(t) is lapse function. In this background the action of a minimally coupled phantom scalar field φ(t)\varphi(t) with constant potential is

S=𝑑t{3aa˙2Na3(φ˙22N+NΛ)},S=\int dt\left\{\ -\frac{3a\dot{a}^{2}}{N}-a^{3}\left(\frac{\dot{\varphi}^{2}}{2N}+N\Lambda\right)\right\}\ , (2)

where we have set the units so 8πG=18\pi G=1. The minus sign in the kinetic term of the scalar action is the difference between the usual scalar field and the phantom scalar field [12]. The canonical Hamiltonian derived from Eq.(2) is

N[Pa212a+Pφ22a3a3Λ].-N\left[\frac{P_{a}^{2}}{12a}+\frac{P_{\varphi}^{2}}{2a^{3}}-a^{3}\Lambda\right]. (3)

With the the following change of variables

x\displaystyle x =\displaystyle= μ1a3/2sin(μφ),\displaystyle\mu^{-1}a^{3/2}\sin{(\mu\varphi)}, (4)
y\displaystyle y =\displaystyle= μ1a3/2cos(μφ),\displaystyle\mu^{-1}a^{3/2}\cos{(\mu\varphi)},

and μ=3/8\mu=\sqrt{3/8}. The Hamiltonian Eq.(3) can be rewritten as a sum of two harmonic oscillators

H=N(12Px2+ω22x2)+N(12Py2+ω22y2),H=N\left(\frac{1}{2}P_{x}^{2}+\frac{\omega^{2}}{2}x^{2}\right)+N\left(\frac{1}{2}P_{y}^{2}+\frac{\omega^{2}}{2}y^{2}\right), (5)

where ω2=34Λ\omega^{2}=-\frac{3}{4}\Lambda. When one considers the usual scalar field, the Hamiltonian is transformed to a “ghost oscillator” which is simply a difference of two harmonic oscillators [26, 29].

3 Deformed Phase Space Model

There exist different approaches to incorporate noncommutativity into physical theories. Particularly, in cosmology there is a broadly explored path to study noncommutativity [20], where noncommutativity is realized in the so called minisuperspace variables. We will follow a deformed phase space approach. The deformation is introduced by the Moyal bracket {f,g}α=fαggαf\{f,g\}_{\alpha}=f\star_{\alpha}g-g\star_{\alpha}f. By substituting the usual product with the Moyal product (fg)(x)=exp[12αaba(1)b(2)]f(x1)g(x2)|x1=x2=x(f\star{g})(x)=\exp{\left[\frac{1}{2}\alpha^{ab}\partial_{a}^{(1)}\partial_{b}^{(2)}\right]}f(x_{1})g(x_{2})|_{x_{1}=x_{2}=x} such that

α=(θijδij+σijδijσijβij).\displaystyle\alpha=\left({\begin{array}[]{cc}\theta_{ij}&\delta_{ij}+\sigma_{ij}\\ -\delta_{ij}-\sigma_{ij}&\beta_{ij}\\ \end{array}}\right). (8)

The 2×22\times 2 antisymmetric matrices θij\theta_{ij} and βij\beta_{ij} represent the noncommutativity in the coordinates and momenta respectively. The α\alpha deformed algebra becomes

{xi,xj}α=θij,{xi,pj}α=δij+σij,{pi,pj}α=βij.\{x_{i},x_{j}\}_{\alpha}=\theta_{ij},\{x_{i},p_{j}\}_{\alpha}=\delta_{ij}+\sigma_{ij},\;\{p_{i},p_{j}\}_{\alpha}=\beta_{ij}. (9)

In this work we use the particular deformations, θij=θϵij\theta_{ij}=-\theta\epsilon_{ij} and βij=βϵij\beta_{ij}=\beta\epsilon_{ij}.

There is an alternative to derive a similar algebra to Eq.(9). For this we follow the procedure given in [24, 26]. Start with the transformation

x^=x+θ2Py,y^=yθ2Px,\displaystyle\widehat{x}=x+\frac{\theta}{2}P_{y},\qquad\widehat{y}=y-\frac{\theta}{2}P_{x},
P^x=Pxβ2y,P^y=Py+β2x,\displaystyle\widehat{P}_{x}=P_{x}-\frac{\beta}{2}y,\qquad\widehat{P}_{y}=P_{y}+\frac{\beta}{2}x, (10)

on the classical phase space variables {x,y,Px,Py}\{x,y,P_{x},P_{y}\}, these are the variables that satisfy the usual Poisson algebra. The new variables satisfy a deformed algebra

{y^,x^}=θ,{x^,P^x}={y^,P^y}=1+σ,{P^y,P^x}=β,\{\widehat{y},\widehat{x}\}=\theta,\;\{\widehat{x},\widehat{P}_{x}\}=\{\widehat{y},\widehat{P}_{y}\}=1+\sigma,\;\{\widehat{P}_{y},\widehat{P}_{x}\}=\beta, (11)

where σ=θβ/4\sigma=\theta\beta/4. Furthermore, as in [24, 26], we assume that the deformed variables satisfy the same relations as their commutative counterparts. The resulting algebra is the same, but the Poisson bracket is different in the two algebras. In Eq.(9), the brackets are the α\alpha deformed ones related to the Moyal product, for the other algebra the brackets are the usual Poisson brackets.

To construct the deformed theory, we start with a Hamiltonian which is formally analogous to Eq.(5) but constructed with the variables that obey the modified algebra Eq.(11), this gives the deformed Hamiltonian

Hnc=12[(Px2+Py2)+2(xPyyPx)+ω~2(x2+y2)],H_{nc}=\frac{1}{2}\left[(P_{x}^{2}+P_{y}^{2})+\ell^{2}(xP_{y}-yP_{x})+\widetilde{\omega}^{2}(x^{2}+y^{2})\right], (12)

where 2\ell^{2} and ω~2\widetilde{\omega}^{2} are given by

2=β+ω2θ1+ω2θ24,ω~2=ω2+β241+ω2θ24.\ell^{2}=\frac{\beta+\omega^{2}\theta}{1+\frac{\omega^{2}\theta^{2}}{4}},~~~~~\tilde{\omega}^{2}=\frac{\omega^{2}+\frac{\beta^{2}}{4}}{1+\frac{\omega^{2}\theta^{2}}{4}}. (13)

There is a significant difference between the transformed Hamiltonian of the scalar field cosmology model [26] and Eq. (12), the deformed Hamiltonian of the phantom cosmology model. Unlike the scalar field case, the crossed term involving position and momentum in the phantom model, corresponds to an angular momentum term. To understand the physics of the deformation, it is necessary to remember that the deformation Eq.(10) defines two physical nonequivalent descriptions, the “C-frame” where the effects of the deformation are interpreted as a commutative space (x,y)(x,y) but with modification of the original interaction and the “NC-frame” where we work with the deformed variables (x^,y^)(\widehat{x},\widehat{y}) and the original interaction. In general, the dynamics in the two frames is different but the physical interpretation in the “C-frame” is easier. In this frame we can interpret the deformed model as a bidimensional harmonic oscillator and the minisuperspace deformation comes into play as an angular momentum term in the Hamiltonian. For this reason we will do the calculations in the “C-frame”.

We obtain the equations of motion from the Hamiltonian Eq.(12), which in the (x,y)(x,y) variables are given by

x˙\displaystyle\dot{x} =\displaystyle= Px122y,y˙=Py+122x,\displaystyle P_{x}-\frac{1}{2}\ell^{2}{y},\qquad\dot{y}=P_{y}+\frac{1}{2}\ell^{2}{x}, (14)
P˙x\displaystyle{\dot{P}_{x}} =\displaystyle= 122Pyω~2x,P˙y=122Pxω~2y,\displaystyle-\frac{1}{2}\ell^{2}P_{y}-\widetilde{\omega}^{2}{x},\qquad{\dot{P}_{y}}=\frac{1}{2}\ell^{2}P_{x}-\widetilde{\omega}^{2}{y},

and get

x¨+2y˙+(ω~244)x\displaystyle\ddot{x}+\ell^{2}\dot{y}+\left(\widetilde{\omega}^{2}-\frac{\ell^{4}}{4}\right)x =\displaystyle= 0,\displaystyle 0, (15)
y¨2x˙+(ω~244)y\displaystyle\ddot{y}-\ell^{2}\dot{x}+\left(\widetilde{\omega}^{2}-\frac{\ell^{4}}{4}\right)y =\displaystyle= 0.\displaystyle 0.

With the transformation z=x+iyz=x+iy we can easily solve Eq.(15). We have three different solutions depending on the sign of ω~2\widetilde{\omega}^{2}. For ω~2>0\widetilde{\omega}^{2}>0, we get

x(t)\displaystyle x(t) =\displaystyle= A1cos[(2/2+|ω~|)t]+B1cos[(2/2|ω~|)t],\displaystyle A_{1}\cos{[(\ell^{2}/2+|\widetilde{\omega}|)t]}+B_{1}\cos{[(\ell^{2}/2-|\widetilde{\omega}|)t]},~~~~~~ (16)
y(t)\displaystyle y(t) =\displaystyle= A1sin[(2/2+|ω~|)t]+B1sin[(2/2|ω~|)t],\displaystyle A_{1}\sin{[(\ell^{2}/2+|\widetilde{\omega}|)t]}+B_{1}\sin{[(\ell^{2}/2-|\widetilde{\omega}|)t]},~~~~~~

A1,B1A_{1},B_{1} are arbitrary constants that, from the Hamiltonian constraint, satisfy (2+2|ω~|)A1(22|ω~|)B1=0(\ell^{2}+2|\widetilde{\omega}|)A_{1}-(\ell^{2}-2|\widetilde{\omega}|)B_{1}=0. When ω~2<0\widetilde{\omega}^{2}<0 we have the solutions

x(t)\displaystyle x(t) =\displaystyle= [A3cosh(|ω~|t)+B3sinh(|ω~|t)]cos(22t),\displaystyle\left[A_{3}\cosh{(|\widetilde{\omega}|t)}+B_{3}\sinh{(|\widetilde{\omega}|t)}\right]\cos{\left(\frac{\ell^{2}}{2}t\right)}, (17)
y(t)\displaystyle y(t) =\displaystyle= [A3cosh(|ω~|t)+B3sinh(|ω~|t)]sin(22t),\displaystyle\left[A_{3}\cosh{(|\widetilde{\omega}|t)}+B_{3}\sinh{(|\widetilde{\omega}|t)}\right]\sin{\left(\frac{\ell^{2}}{2}t\right)},

for this case the integrating constants satisfy the condition A32=B32A_{3}^{2}=B_{3}^{2}. The last solution is for the case ω~=0\widetilde{\omega}=0,

x(t)\displaystyle x(t) =\displaystyle= (A2+B2t)cos(22t),\displaystyle\left(A_{2}+B_{2}t\right)\cos{\left(\frac{\ell^{2}}{2}t\right)}, (18)
y(t)\displaystyle y(t) =\displaystyle= (A2+B2t)sin(22t),\displaystyle\left(A_{2}+B_{2}t\right)\sin{\left(\frac{\ell^{2}}{2}t\right)},

the constraint imposes B2=0B_{2}=0. The condition arises when the cosmological constant satisfies 2β=3Λ2\beta=-3\Lambda.

In analyzing the evolution of the Universe, we will focus our attention on its volume and calculate a3(t)a^{3}(t) in the two frames. Let us start in the “C-frame”, using the solutions for the model, Eq.(16), Eq.(17) and Eq.(18), together with the change of variables in Eq.(4) we get

a3(t)={V0(1)+V1(1)cos2(ω~t),ω~2>0,V0(2),ω~2=0,V0(3)e2ω~t,ω~2<0,a^{3}(t)=\left\{\begin{array}[]{ccc}&V^{(1)}_{0}+V^{(1)}_{1}\cos^{2}{(\widetilde{\omega}t)},~~~~~~~~~~~\widetilde{\omega}^{2}>0,\\ &\\ &V^{(2)}_{0},~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\widetilde{\omega}^{2}=0,\\ &\\ &V^{(3)}_{0}e^{2\tilde{\omega}t},~~~~~~~~~~~~~~~~~~~~~~~~~\widetilde{\omega}^{2}<0,\\ \end{array}\right. (19)

where V0(i),V1(i)V^{(i)}_{0},V^{(i)}_{1} are constructed from the integrating constants. For the first case, V0(1)=(A1B1)2V^{(1)}_{0}=(A_{1}-B_{1})^{2} and V1(1)=4A1B1V^{(1)}_{1}=4A_{1}B_{1}, we must remember that the constants are related. For the second case the volume is the square of a constant. Finally for the last case we have V0(3)=A32V^{(3)}_{0}=A_{3}^{2}, but if we take A3=B3A_{3}=-B_{3} we have an exponentially decaying volume.

To study the dynamics in the “NC-frame” we use the “C-frame” solutions. The volume is constructed with the “NC-frame” variables, a^3(t)=m2(x^2+y^2)\widehat{a}^{3}(t)=m^{2}(\widehat{x}^{2}+\widehat{y}^{2})

a^3(t)={V^0(1)+V^1(1)cos2(ω~t),ω~2>0,V^0(2),ω~2=0,V^1(3)e2ω~t,ω~2<0,\widehat{a}^{3}(t)=\left\{\begin{array}[]{ccc}&\widehat{V}^{(1)}_{0}+\widehat{V}^{(1)}_{1}\cos^{2}{(\widetilde{\omega}t)},~~~~~~~~~~~\widetilde{\omega}^{2}>0,\\ &\\ &\widehat{V}^{(2)}_{0},~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\widetilde{\omega}^{2}=0,\\ &\\ &\widehat{V}^{(3)}_{1}e^{2\tilde{\omega}t},~~~~~~~~~~~~~~~~~~~~~~~~~~\widetilde{\omega}^{2}<0,\\ \end{array}\right. (20)

where V^0(i),V^1(i)\widehat{V}^{(i)}_{0},\widehat{V}^{(i)}_{1} are constants related to V0(i),V1(i)V^{(i)}_{0},V^{(i)}_{1}. It is important to mention that even if we have a positive volume in the “C-frame” we might need to impose some restriction in order to have a valid description in the “NC-frame”. In order to determine the valid range of the deformation parameters we use the principle used in [29] which states the following, “Deformed phase space models are only valid when the NC-frame and C-frame descriptions are physically equivalent”. To impose that the description in the two frames should be physically equivalent, we start by analyzing the expression for the “C-frame” volume. It is easy to show that the volume is always positive. In the first case one obtains that V0(1)+V1(1)=(A1+B1)20V^{(1)}_{0}+V^{(1)}_{1}=(A_{1}+B_{1})^{2}\geq 0 so is positive for any tt. For the second case the volume is also positive, and finally for the third case the constant V1(3)>0V^{(3)}_{1}>0. Then we conclude that for the three cases, the “C-frame” volume is always positive for any choice of the integration constants. The constants in the “NC-frame” solutions are related to the “C-frame” constants as follows. For the first case, A^1=(1+|ω~|θ/2)A1\widehat{A}_{1}=(1+|\widetilde{\omega}|\theta/2)A_{1} and B^1=(1|ω~|θ/2)B1\widehat{B}_{1}=(1-|\widetilde{\omega}|\theta/2)B_{1}, for the second case the volume is the same and for the last case the constant in the “NC-frame” is related by A^32=(1+ω~2θ2/4)A32\widehat{A}^{2}_{3}=(1+\tilde{\omega}^{2}\theta^{2}/4)A_{3}^{2}.

We see from Eq. (LABEL:CVolume) and Eq.(LABEL:NCVolume) that for θ=0\theta=0, the two frames are physically equivalent. We can conclude that θ=0\theta=0 and β0\beta\neq 0 are valid values for the deformation parameters. Now we consider the general case θ0\theta\neq 0 and β0\beta\neq 0. It is clear that in the first two cases, that the volume has the same physical behavior in the two frames then there are no restriction for the deformation parameters. For the third case, we have the same behavior as long as we satisfy the condition 1+θ2ω~24>01+\frac{\theta^{2}\widetilde{\omega}^{2}}{4}>0. Furthermore, the evolution is that of an exponential scale factor for an accelerating Universe. Then comparing with a de-Sitter Universe we find (only for the third case), an effective positive cosmological constant

ΛEff=13(3Λβ213θ2Λ16).\Lambda_{Eff}=\frac{1}{3}\left(\frac{3\Lambda-\beta^{2}}{1-\frac{3\theta^{2}\Lambda}{16}}\right). (21)

This effective cosmological constant is determined by the deformation parameters as well as the value of Λ\Lambda. Since current observations point to a positive cosmological constant. This gives and extra condition that restricts the values of the deformation parameters by imposing that Λeff>0\Lambda_{eff}>0. After fixing the value of Λ\Lambda, we have three definite restrictions:

  1. i)

    ΛEff>0\Lambda_{Eff}>0.

  2. ii)

    f(β,Λ)=3Λβ2f(\beta,\Lambda)=3\Lambda-\beta^{2} and g(θ,Λ)=13θ2Λ16g(\theta,\Lambda)=1-\frac{3\theta^{2}\Lambda}{16} must have the same sign.

  3. iii)

    1+θ2ω~24>01+\frac{\theta^{2}\widetilde{\omega}^{2}}{4}>0.

From ii), we can see in Fig.(1) and Fig.(2) the regions where these inequalities are satisfied. We did not consider Λ0\Lambda\leq 0, since condition ii) would not be satisfied, for this reason the plots only show the region where Λ>0\Lambda>0.

Refer to caption
Figure 1: This figure shows, in the gray region, the permitted values of β\beta and Λ\Lambda that satisfy f(β,Λ)>0f(\beta,\Lambda)>0. The white region corresponds to he values where f(β,Λ)<0f(\beta,\Lambda)<0.
Refer to caption
Figure 2: The gray region in this figure shows the values of θ\theta and Λ\Lambda that satisfy g(θ,Λ)>0g(\theta,\Lambda)>0. The white region corresponds to the values where g(θ,Λ)<0g(\theta,\Lambda)<0.

We want to emphasize that for the phantom scalar field we cannot define an effective cosmological constant when Λ=0\Lambda=0, as it was done for the usual scalar field in [26]. When setting Λ=0\Lambda=0 we will have ω~2>0\widetilde{\omega}^{2}>0 meaning that we are in the first case and therefore we will not have an exponential behavior. For this model the deformation maps an initial positive cosmological constant model to another positive cosmological constant model, the two cosmological constants have different values but they are both positive. This fact also leads to a more restrictive space of parameters than the one analyzed in [29].

4 Noncommutative quantum model

Now we proceed with the quantization of the deformed Hamiltonian. The usual canonical quantization approach Hncψ=0H_{nc}\psi=0, gives the corresponding WDW equation. The Hamiltonian HncH_{nc}, explicitly contains the angular momentum operator (xPyyPx)=Lz(xP_{y}-yP_{x})=L_{z} and can be written in terms of the Hamiltonian of a 2-dimensional harmonic oscillator HxyH_{xy}, with frequency ω~2\widetilde{\omega}^{2} plus an angular momentum term, that is Hnc=Hxy+LzH_{nc}=H_{xy}+L_{z}. Using the similarity of the WDW equation, with the 2-dimensional harmonic oscillator, we will follow the approach in [30, 31]. First we write the eigenvalue equation for the harmonic oscillator, then we will find the energy eigenvalues and finally impose the zero energy condition. After writing the quantum equation in polar coordinates, we obtain the equation for the energy eigenstates, and is given by

22[2ψr2+1rψr+1r22ψϕ2]+ω~22r2ψi22ψϕ=Eψ.-\frac{\hbar^{2}}{2}\left[\frac{\partial^{2}\psi}{\partial r^{2}}+\frac{1}{r}\frac{\partial\psi}{\partial r}+\frac{1}{r^{2}}\frac{\partial^{2}\psi}{\partial\phi^{2}}\right]+\frac{\widetilde{\omega}^{2}}{2}r^{2}\psi-\frac{i\hbar\ell^{2}}{2}\frac{\partial\psi}{\partial\phi}=E\psi. (22)

It is well known that the Hamiltonian of the 2d harmonic oscillator commutes with the angular momentum. This means that we have [Hxy,Lz]=0[H_{xy},L_{z}]=0, allowing us to write the wave function with eigenstates common to both operators. This breaks the degeneracy in energy eigenstates and every state is uniquely specified by the quantum numbers nn and mm. Solving Eq.(22) we find the wave function of the Universe ψn,m(r,ϕ)\psi_{n,m}(r,\phi)

ψn,m(r,ϕ)=N(ω~r2)|m|2Ln|m|(ω~r2)eimϕeω~r22.\psi_{n,m}(r,\phi)=N\left(\frac{\widetilde{\omega}r^{2}}{\hbar}\right)^{\frac{|m|}{2}}L^{|m|}_{n}\left(\frac{\widetilde{\omega}r^{2}}{\hbar}\right)e^{im\phi}e^{-\frac{\tilde{\omega}r^{2}}{2\hbar}}. (23)

where n=0,1,2,;m=0,±1,±2,,±n,n=0,1,2,...\>;~m=0,\pm 1,\pm 2,...,\pm n, and Lnα(z)L_{n}^{\alpha}(z) are the generalized Laguerre polynomials. The zero energy condition E=0E=0 from the WDW equation is trivially satisfied by =ω~=0\ell=\widetilde{\omega}=0. For 0,ω~0\ell\neq 0,~\widetilde{\omega}\neq 0 (for instance for a single state of quantum numbers nn and mm) we get the general equation

2ω~=2(2n+|m|+1)m.\frac{\ell^{2}}{\widetilde{\omega}}=-\frac{2(2n+|m|+1)}{m}. (24)

By using the definitions for 2\ell^{2} and ω~2>0\widetilde{\omega}^{2}>0 we can write this equation as

16β12θΛ(163θ2Λ)(4β212Λ)=2(2n+|m|+1)m.\frac{16\beta-12\theta\Lambda}{\sqrt{(16-3\theta^{2}\Lambda)(4\beta^{2}-12\Lambda)}}=-\frac{2(2n+|m|+1)}{m}. (25)

With respect to the condition in Eq. (25), we need to find the values of θ\theta, β\beta and Λ\Lambda to satisfy the condition. The first thing to notice is that the right hand side of this equation takes only discrete values and to determine if there are possible values of the parameters to satisfy the equation we may proceed as follows. Analyze the behavior of the function h(θ,β,Λ)=16β12θΛ163θ2Λ4β212Λh(\theta,\beta,\Lambda)=\frac{16\beta-12\theta\Lambda}{\sqrt{16-3\theta^{2}\Lambda}\sqrt{4\beta^{2}-12\Lambda}} and then see whether or not, there are permitted values of the parameters for which the function takes the value h(θ,β,Λ)=2(2n+|m|+1)mh(\theta,\beta,\Lambda)=-\frac{2(2n+|m|+1)}{m}. The quantum solution is possible when ω~2>0\widetilde{\omega}^{2}>0, which corresponds to the first solution in Eq. (LABEL:CVolume) and Eq.(LABEL:NCVolume). Therefore, the space of parameters is already depicted in Fig.(1) and Fig.(2), corresponding to the condition that the functions f(β,Λ)f(\beta,\Lambda) and g(θ,Λ)g(\theta,\Lambda) have opposite sign. The valid values of parameters θ\theta, β\beta and Λ\Lambda, corresponds to the gray region of Fig. 1 and white region of Fig. 2, or vice versa. We use the inverse transformation of Eq.(4) with Eq. (23)

r2=a3μ2,ϕ=π2μφ,r^{2}=\frac{a^{3}}{\mu^{2}},~~~~~\phi=\frac{\pi}{2}-\mu\varphi, (26)

to find the wave function ψn,m(a,φ)\psi_{n,m}(a,\varphi) in terms of the scale factor aa and the phantom field φ\varphi.

We can construct a general wave packet

ψ(a,φ)=nNcn,mψ(n,m)(a,φ),\psi(a,\varphi)=\sum_{n}^{N}c_{n,m}\psi_{(n,m)}(a,\varphi), (27)

with mm fixed and cn,mc_{n,m} arbitrary constants as long as the condition E=nNE(n,m)=0E=\sum_{n}^{N}E_{(n,m)}=0 also holds. When this is the case, it is always possible to find values of θ\theta, β\beta and Λ\Lambda to satisfy the constriction Eq. (24), under which these values are fixed unambiguously. For a given value of mm, the sum in Eq.(27) runs over odd/even values of nn depending on the odd/even value of mm respectively. The wave packet is normalizable by the use of the orthogonality relation of generalized Laguerre polynomials. The total energy of the wave packet will satisfy the generalized zero energy condition Eq.(24), but replacing nn by NN. From the arguments given above, we conclude that it is always possible to find, for a given Λ>0\Lambda>0 values for θ\theta and β\beta for which the condition holds.

4.1 Classical paths from a WKB approximation

One can apply a WKB type approximation on the WDW equation HncΨ=0H_{nc}\Psi=0, to see if we get the classical solutions. First we notice that the effective classical cosmological constant given in Eq. (21) is the same quantum cosmological constant in the Hamiltonian HncH_{nc}. Remembering that in the commutative Hamiltonian 4ω2=3Λ4\omega^{2}=-3\Lambda, for the deformed Hamiltonian 4ω~2=3ΛEff4\widetilde{\omega}^{2}=-3\Lambda_{Eff}. This is the generalized deformed frequency and using the definition of ω~2\widetilde{\omega}^{2} we get the same effective cosmological constant.

We now apply the standard procedure for the semiclassical approximation. We propose a wave function of the form

ψeiS1(x)+iS2(y),\psi\propto e^{\frac{i}{\hbar}S_{1}(x)+\frac{i}{\hbar}S_{2}(y)}, (28)

which upon substitution in HncΨ=0H_{nc}\Psi=0, in the limit 0\hbar\rightarrow 0 and using the approximation

(S1(x)x)2>>2S1(x)x2,(S2(y)y)2>>2S2(y)y2,\left(\frac{\partial S_{1}(x)}{\partial x}\right)^{2}>>\frac{\partial^{2}S_{1}(x)}{\partial x^{2}},~~~\left(\frac{\partial S_{2}(y)}{\partial y}\right)^{2}>>\frac{\partial^{2}S_{2}(y)}{\partial y^{2}}, (29)

we get the Einstein-Hamilton-Jacobi (EHJ) equation

(S1(x)x)2+(S2(y)y)2+ω~2(x2+y2)\displaystyle\left(\frac{\partial S_{1}(x)}{\partial x}\right)^{2}+\left(\frac{\partial S_{2}(y)}{\partial y}\right)^{2}+\widetilde{\omega}^{2}(x^{2}+y^{2})
+\displaystyle+ 2x(S2(y)y)2y(S1(x)x)=0.\displaystyle\ell^{2}x\left(\frac{\partial S_{2}(y)}{\partial y}\right)-\ell^{2}y\left(\frac{\partial S_{1}(x)}{\partial x}\right)=0.

With the identification S1(x)x=Px\frac{\partial S_{1}(x)}{\partial x}=P_{x}, S2(y)y=Py\frac{\partial S_{2}(y)}{\partial y}=P_{y}, the EHJ equation takes the form

x˙2+y˙2+(ω~244)(x2+y2)=0.\dot{x}^{2}+\dot{y}^{2}+\left(\widetilde{\omega}^{2}-\frac{\ell^{4}}{4}\right)\left(x^{2}+y^{2}\right)=0. (31)

Now we take the time derivative of the EHJ equation and divide the result by x˙y˙\dot{x}\dot{y} to get

x¨y˙+(ω~244)xy˙+y¨x˙+(ω~244)yx˙=0,\frac{\ddot{x}}{\dot{y}}+\left(\widetilde{\omega}^{2}-\frac{\ell^{4}}{4}\right)\frac{x}{\dot{y}}+\frac{\ddot{y}}{\dot{x}}+\left(\widetilde{\omega}^{2}-\frac{\ell^{4}}{4}\right)\frac{y}{\dot{x}}=0, (32)

we can separate the expression by equating the first two terms to 2-\ell^{2} and the last two to 2\ell^{2}. The two resulting equations are the equations of motion in Eq.(15). We can be confident that the same classical paths arise in the classical limit of the WDW equation obtained from the deformed Hamiltonian.

5 Conclusions

In this paper we have studied the effects of the minisuperspace phase space deformations in the dynamics of phantom scalar cosmology. The deformation was applied on the minisuperspace variables. From the deformed Hamiltonian we derived the equations of motion. It is known that from the deformation two generally nonequivalent descriptions arise. Using the solutions to Eq.(15), we calculate the volume of the Universe in the “C-frame” and in the “NC-frame”. To constrain the deformation parameters β\beta and θ\theta, we impose that for phase space deformations to be valid, the physical description in the two frames must be equivalent. From the value of ω~2\widetilde{\omega}^{2} three different cases appear. For ω~2>0\widetilde{\omega}^{2}>0, we get a bouncing Universe and the solutions in the two frames are equivalent. For the last case, when the Hamiltonian constraint is applied, we have a de Sitter Universe and the parameter ω~2\widetilde{\omega}^{2} can be interpreted as an effective positive cosmological constant ΛEff\Lambda_{Eff}. The value of ΛEff\Lambda_{Eff} depends on the deformation parameters θ\theta and β\beta as well as the original Λ\Lambda. If we set the original cosmological constant Λ=0\Lambda=0, then ω~2>0\widetilde{\omega}^{2}>0 and is not possible to have ΛEff0\Lambda_{Eff}\geq 0. This imposes a stronger restriction on the values of β\beta and θ\theta and the parameter space is even more restricted, but even when this is the case we have a region of values where the two descriptions are physically equivalent.

We quantized the model, obtaining the corresponding WDW equation and found the wave function of the Universe. We derived a relation that fixes the values of the parameters θ\theta, β\beta and Λ\Lambda in order to satisfy the zero energy condition of the WDW equation. We showed that the deformed phantom cosmological model is physically viable. Also, when ω~2>0\widetilde{\omega}^{2}>0 the quantization can be achieved and the introduction of the deformation breaks the degeneracy of the quantum states. A general condition on the parameters θ\theta and β\beta arises from the zero energy condition that can be satisfied by the general wave packets. It is worth to mention that there is a nice physical interpretation in the “C-frame”. It is related to the introduction to an effective perpendicular and constant “magnetic field” B\vec{B}. This can be seen by taking the bidimensional harmonic oscillator coupled to the vector potential A=(B2y,B2x)\vec{A}=(\frac{B}{2}y,-\frac{B}{2}x) where the magnitude of the field will be B=2B=\ell^{2}. Then we can see that the zero energy condition is related to the manipulation of this field and corresponds to a modification in the deformation parameters (θ,β)(\theta,\beta). Finally, we performed the WKB approximation and show that the same classical solutions arise in the classical limit of the WDW equation of the corresponding deformed Hamiltonian.

Acknowledgments

This work is supported by CONACyT research grants 167335, 257919, 258982.

J. L. López is supported by PRODEP postdoctoral grant 17-1947. This work is part of the PRODEP research network “Gravitación y Física Matemática”.

M. S. is supported by DAIP1107/2016 and by the CONACyT program “Estancias sabáticas en el extranjero”, grant 31065.

References

References

  • [1] B. Ratra and P. J. E. Peebles, Phys. Rev. D 37 3406 (1988).
  • [2] P. J. E. Peebles and B. Ratra, Rev. Mod. Phys.  75, 559 (2003).
  • [3] L. P. Chimento and A. S. Jakubi, Int. J. Mod. Phys. D 5 71 (1996).
  • [4] E. J. Copeland, M. Sami, and S. Tsujikawa, Int. J. Mod. Phys. D 15 1753 (2006).
  • [5] K. Bamba, S. Capozziello, S. Nojiri, and S. D. Odintsov, Astrophys. Space.  Sci. 342 155 (2012).
  • [6] T. Padmanabhan and T. Roy Choudhury, Phys. Rev. D 66 081301 (2002).
  • [7] V. A. Rubakov, Theor. Math. Phys. 149 1651 (2006).
  • [8] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J.  517, 565 (1999).
  • [9] A. G. Riess et al., Astron. J.  116 1009 (1998).
  • [10] A. G. Riess et al, Astrophys. J.  560 49 (2001).
  • [11] Narciso Benitez et al., Astrophys. J.  577 L1 (2002).
  • [12] Mariusz Dabrowski, Phys. Rev.  D 74 044022 (2006).
  • [13] Robert R. Caldwell and Marc Kamionkowski, Ann. Rev. Nucl. Part. Sci. 59 397 (2009).
  • [14] J. M. Cline, S. Jeon, and G. D. Moore, Phys. Rev. D 70 043543 (2004).
  • [15] W. Hu, Phys.  Rev.  D 71 047301 (2005).
  • [16] H. Garcia-Compean, O. Obregon, C. Ramirez, and M. Sabido, Phys. Rev. D 68 045010 (2003).
  • [17] X. Calmet and A. Kobakhidze, Phys. Rev. D 72 045010 (2005).
  • [18] P. Aschieri, M.Dimitrijevic, F.Meyer and J.Wess, Class. Quant. Grav.  23 1883 (2006).
  • [19] T. Ohl and A. Schenkel, JHEP 0910 052 (2009).
  • [20] H. Garcia-Compean, O. Obregon, and C. Ramirez, Phys. Rev. Lett. 88 161301 (2002).
  • [21] E. Mena, O. Obregon and M. Sabido, Int. J. Mod. Phys. D 18 (2009) 95.
  • [22] L. O. Pimentel and O.Obregon, Gen. Rel.  Grav. 38 553 (2006).
  • [23] C. Bastos, O. Bertolami, N. Costa Dias and J. Nuno Prata, Phys. Rev. D 78 (2008) 023516.
  • [24] B. Vakili, C. I. Azad, P. Pedram, and S. Jalalzadeh, Phys. Lett. B 687 119 (2010).
  • [25] B. Malekolkalami, K. Atazadeh and B. Vakili, Phys. Lett. B 739, 400 (2014).
  • [26] S. Pérez-Payán, M. Sabido and C. Yee-Romero, Phys. Rev. D 88 027503 (2013).
  • [27] M. Sabido and C. Yee-Romero, Phys. Lett. B 757 (2016) 57.
  • [28] A. Crespo-Hernandez, E. A. Mena-Barboza and M. Sabido, Entropy 19, 91 (2017).
  • [29] Sinuhe Pérez-Payán, M. Sabido, E. Mena and C. Yee-Romero, Adv. High. Energy Phys. 2014 958137 (2014)
  • [30] S. S. Gousheh and H. R. Sepangi, Phys. Lett. A 272, 304 (2000).
  • [31] S. Jalalzadeh, F. Ahmadi and H. R. Sepangi, JHEP 0308, 012 (2003).