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

Casimir effect in DFR space-time

E. Harikumar and Suman Kumar Panja [email protected]@gmail.com School of Physics, University of Hyderabad,
Central University P.O, Hyderabad-500046, Telangana, India
Abstract

Non-Commutative space-time introduces a fundamental length scale suggested by approaches to quantum gravity. Here we report the analysis of the Casimir effect for parallel plates separated by a distance of LL using a Lorentz invariant scalar theory in a non-commutative space-time (DFR space-time), both at zero and finite temperatures. This is done in two ways; one when the additional space-dimensions introduced in DFR space-time are treated as extra dimensions but on par with usual space-dimension and in the second way, the additional dimensions are treated as compact dimensions. Casimir force obtained in the first approach coincides with the result in the extra-dimensional commutative space-time and this is varying as 1L5\frac{1}{L^{5}}. In the second approach, we derive the corrections to the Casimir force, which is dependent on the separation between the plate, LL and on the size of the extra compactified dimension, RR. Since correction terms are very small, keeping only the most significant terms of these corrections, we show that for certain values of the R, the corrections due to non-commutativity makes the force between the parallel plates more attractive, and using this, we find lower bound on the value of RR. We show here that the requirement of the Casimir force and the energy to be real, impose the condition that the weight function used in defining the DFR action has to be a constant. At zero temperature, we find correction terms due to non-commutativity, depend on LL and RR dependent modified Bessel functions K1K_{1} and K2K_{2}, with coefficients that vary as 1LR3\frac{1}{LR^{3}} and 1L2R2\frac{1}{L^{2}R^{2}}, respectively . For finite temperature, the Casimir force has correction terms that scale as 1L\frac{1}{L} and 1L3\frac{1}{L^{3}} in high-temperature limit and as 1L2\frac{1}{L^{2}} and 1L4\frac{1}{L^{4}} in the low-temperature limit.

1 Introduction

Though it is the oldest known force, the nature of gravity is still unknown on microscopic scale. Understanding the nature of gravity at the quantum regime is one of the pressing issues in Physics. Many approaches have been developed to study quantum theory of gravity. Some of these approaches are string theory, loop gravity, emergent gravity, and non-commutative geometry[1, 2, 3, 4, 5, 6, 7]. All of these approaches introduce a fundamental length scale, below which quantum effect of gravity comes into play[3, 4]. Non-commutative geometry provides a path to incorporate such a fundamental length scale[1, 2, 3, 4]. Non-commutative geometry also appears in the low energy limit of string theory [2]. In low energy limit, quantum gravity model based on spin-form reduce to a field theory in a non-commutative space-time known as κ\kappa-space-time [8].

In recent times non-commutative space-times and models on such space-times are being investigated intensively. One of the initial motivations behind introducing non-commutative space-time was to remove divergence from QFT. In [9], non-commutative space-time was proposed, but it was shown that the divergences were not completely removed in this non-commutative theory[10]. Research activities in different areas, such as non-commutative geometry, string theory, fuzzy sphere, etc., rekindled the interest in non-commutative space-times and study of physics on such space-times[1, 2, 3, 4, 5, 6, 7, 11]. One of the well studied non-commutative space-time is the Moyal space-time, whose coordinates satisfy

[x^μ,x^ν]=iθμν,[\hat{x}^{\mu},\hat{x}^{\nu}]=i\theta^{\mu\nu},

where θμν\theta^{\mu\nu} is antisymmetric, constant tensor with dimension of (length)2(length)^{2}. This constant θμν\theta^{\mu\nu} leads to the violation of Lorentz invariance in field theory [2] and breaking of rotational symmetry in non-relativistic theory. The violation of Lorentz symmetry makes definition of particles problematic and also it brings vacuum birefringence in theory[12].

To recover Lorentz symmetry in non-commutative theory, a transformation rule has been assigned to the non-commutative parameter θμν\theta^{\mu\nu} under Lorentz transformation[4] and further the components of the non-commutative parameter θμν\theta_{\mu\nu} have been elevated to coordinates of the space-time, thereby increasing the dimension of the underlying space-time[13]. This space-time is called DFR space-time [4] and associated symmetry algebra is known as DFR algebra. The coordinates of the phase space corresponding to DFR space-time are xix_{i}, pip_{i}, θi\theta_{i} and kik_{i}. To get complete Fourier decomposition of field and to maintain Lorentz invariance in DFR space-time, canonical conjugate momentum corresponding to non-commutative coordinate θμν\theta^{\mu\nu} have been introduced [14, 15, 16, 17]. The associated symmetry algebra is called DFRA algebra.

Since quantum theories of non-commutative space-time have an inherent feature of introducing length parameters (along with nonlocality and nonlinearity), it is of intrinsic interest to study the implications of this length scale in physical phenomena. One such phenomenon where length scale plays a significant role in the commutative space-time is the Casimir effect[18, 19]. It was shown that when two conducting plates are placed parallel to each other at a small separation, they experience an attractive force, and this force depends on the separation between the plate[19]. It was shown that the zero-point energy of the electromagnetic field changes when the plates are introduced and the attractive force arises due to the change in the zero-point energy. The boundary condition imposed on the electromagnetic field by the introduction of plates is different from that in their absence, and this changes the zero-point energy in two situations. This change in the zero-point energy was shown to result in an attractive force between parallel conductors in [18]. Subsequently, this effect was studied by calculating the change in the zero-point energy of various quantized fields such as a real scalar field and fermionic field, for conductors with different geometries [20, 19]. The electromagnetic field at either side of the plates consists of two modes (TE and TM), while the real scalar field has only one degree of freedom but obeys the same boundary conditions as the electromagnetic field. The Casimir energy and force calculated using the scalar field was shown to be half of the corresponding values obtained by considering the electromagnetic field, and this numerical factor difference is due to the differences in the number of modes associated with the scalar and the electromagnetic fields[19, 20, 21].

Modifications to the Casimir effect due to different dimensions of space-time, due to the presence of medium between plates, due to thermal effect, etc., using real scalar field theory have been studied[20, 22]. Experimental investigations of Casimir force for separation of few micrometers between the plates have been reported [23, 24, 25]. Apart from the experimental study for parallel conducting plates, the Casimir effect has been investigated for other geometrical configurations [26, 27, 28]. Using these experimental results, constraints on the Yukawa type correction to Newtonian gravity are obtained [29].

In [30, 31, 32, 33, 34], Casimir effect has been studied in different non-commutative space-times such as Moyal space-time [2] and κ\kappa-deformed space-time [35], using scalar field theory. For even-dimensional Moyal space-time, using coherent state approach and smeared boundary condition, Casimir force is calculated between two parallel conducting plates [30, 31] using scalar field theory. In [32], Casimir energy in 2+1 dimensional non-commutative space-time with non-trivial topologies was derived and analyzed. In [36, 37], the Casimir effect on compact non-commutative space was studied. In [33, 34], the Casimir effect was investigated for κ\kappa-deformed scalar theory. Casimir effect was investigated using the Green’s function approach, and a bound on the non-commutative length scale was calculated by comparing it with experimental results[34]. In κ\kappa-space-time, space coordinates commute among themselves, but commutation relation between time coordinate and space coordinates vary with space coordinates, i.e., they satisfy

[xi,xj]=0,[x0,xi]=axi,a=1κ.[x^{i},x^{j}]=0,~{}~{}~{}[x^{0},x^{i}]=ax^{i},~{}~{}~{}a=\frac{1}{\kappa}. (1.1)

where aa is the deformation parameter. This indicates that as in Moyal space-time, in the κ\kappa space-time also the Lorentz symmetry is broken. But this problem is avoided in DFR space-time with introduction of transformation for the non-commutative parameter θμν\theta^{\mu\nu} [4]. Latter, components of θμν\theta_{\mu\nu} have been promoted to coordinates, enlarging the dimension of the space-time[13]. In commutative space-time, Casimir force in presence of extra dimensions has been explored in recent times[20, 38, 39, 40, 41]. Thus it is of intrinsic interest to study the Casimir effect in DFR space-time which has extra dimensions.

In the commutative space-time, Casimir effect between two parallel plates has been evaluated by modelling the plates with two δ\delta-function potentials at x=0x=0 and x=Lx=L [42, 19], respectively. Casimir force on the plate at x=Lx=L has been derived by taking difference of vacuum expectation values of Energy-Momentum tensor on either side of the plate. For massless scalar field theory, xxxx component of Energy-Momentum tensor was found at a point, just left to the plate at x=Lx=L, i.e., T^xx|x=L\hat{T}^{xx}\Big{|}_{x=L^{-}} and at another point, just right to the plate at x=Lx=L, i.e., T^xx|x=L+\hat{T}^{xx}\Big{|}_{x=L^{+}}. Then force was calculated at x=Lx=L by taking difference between the vacuum expectation value of the stress-tensor, in the strong interaction limit of δ\delta-function potential as,

F^γ,γ=<Txx^>|x=L<Txx^>|x=L+.\hat{F}_{\gamma,\gamma^{\prime}\rightarrow\infty}=<\hat{T^{xx}}>\Big{|}_{x=L^{-}}-<\hat{T^{xx}}>\Big{|}_{x=L^{+}}. (1.2)

One way to derive vacuum expectation value of stress tensor was by rewriting Energy-Momentum tensor as an operator acting on the vacuum expectation value of the time ordered product of fields at nearby points, xx and xx^{\prime} as,

<T^x,xμν>\displaystyle<\hat{T}^{\mu\nu}_{x,x^{\prime}}> =\displaystyle= <O^μν(,)T(ϕϕ)>\displaystyle<\hat{O}^{\mu\nu}(\partial,\partial^{\prime})T(\phi\phi^{\prime})>
=\displaystyle= O^μν(,)<T(ϕϕ)>\displaystyle\hat{O}^{\mu\nu}(\partial,\partial^{\prime})<T(\phi\phi^{\prime})>

where O^μν(,)\hat{O}^{\mu\nu}(\partial,\partial^{\prime}) is a derivative operator, ϕ=ϕ(x)\phi=\phi(x) and ϕ=ϕ(x)\phi^{\prime}=\phi(x^{\prime}). Since

<T(ϕ(x)ϕ(x))>=iG(x,x),<T(\phi(x)\phi(x^{\prime}))>=iG(x,x^{\prime}), (1.3)

where G(x,x)G(x,x^{\prime}) is the Green’s function. One can relate the vacuum expectation value of the Energy-Momentum tensor to the Green’s function. Thus to calculate RHS of Eqn.(1.2), Green’s functions in different regions are first derived by solving Euler-Lagrange equation of scalar field theory with Dirichlet boundary condition and using these Green’s functions, Casimir force is obtained. Casimir energy is calculated by integrating the Casimir force over the separation distance between two parallel plates, as

Eγ,γ=F^γ,γ𝑑xE_{\gamma,\gamma^{\prime}\rightarrow\infty}=-\int\hat{F}_{\gamma,\gamma^{\prime}\rightarrow\infty}dx (1.4)

where γ\gamma and γ\gamma^{\prime} are coupling strengths of δ\delta-function potentials, modelling the plates and γ,γ\gamma,\gamma^{\prime}\rightarrow\infty denote the strong interaction limit. The Casimir force and energy calculated using the scalar field differs from the corresponding results obtained using electromagnetic field by a factor of 12\frac{1}{2}.

In this paper, we study the Casimir effect between two parallel plates in DFR space-time by analyzing DFRA scalar field theory [14, 15, 16, 17]. Since the non-commutative coordinates do not show up in the commutative limit, and our aim is to investigate the corrections to the Casimir force in the commutative space-time, we consider the situation where the plates are kept only in the x-directions (commutative directions). We start with the action of 1010-dimensional DFRA scalar field theory. For unitarity of this field theory one sets temporal part of non-commutative coordinate, θμν\theta^{\mu\nu} to be vanishing, i.e., θ0i=0\theta^{0i}=0 and further we take definitions, θi=12ϵijlθjl\theta^{i}=\frac{1}{2}\epsilon^{ijl}\theta_{jl} and ki=12ϵijlkjlk^{i}=\frac{1}{2}\epsilon^{ijl}k_{jl} [45, 43, 44]. After these, we get the action of 77-dimensional DFRA scalar field theory. We start our study from the Lagrangian of DFRA scalar field theory with plates modelled through δ\delta-function potentials, as in the commutative space-time. Next, we vary the corresponding action and derive the general form of the Euler-Lagrange equation and components of the Energy-Momentum tensor. Without loss of generality, we restrict our attention to 4+14+1 dimensions (xμ,θ1)(x^{\mu},\theta^{1}), and study the Casimir effect in two ways. In first way, we treat extra θ\theta-dimension due to noncommutativity of the 4+14+1-dimensional DFR space time, in the same footing as transverse directions, yy and zz. Then following the standard calculational procedure given in [20, 19] we show that the Casimir force vary as 1L5\frac{1}{L^{5}}, where LL is the distance between the parallel plates. This force expression exactly coincides with result given in [19] for extra-dimensional commutative space-time. In second approach, we compactify the new extra dimension θ\theta introduced by noncommutativity in DFR space-time. After compactifying the extra dimension, we solve the equation of motion using Dirichlet boundary condition and get Green’s functions in different regions. We also derive the xxxx component of the Energy-Momentum tensor at either side of the plate at x=Lx=L. Then we obtain the Casimir force by taking the difference between vacuum expectation values of stress tensors at x=Lx=L as shown in Eqn.(1.2) for the commutative case. For this, we use the relation between Green’s function and time-ordered product of fields in the nearby points in 4+14+1-dimensional DFR space-time, viz;

<T(ϕ(x,η)ϕ(x,η))>=iG(x,x;η,η),<T(\phi(x,\eta)\phi(x^{\prime},\eta^{\prime}))>=iG(x,x^{\prime};\eta,\eta^{\prime}), (1.5)

where G(x,x;η,η)G(x,x^{\prime};\eta,\eta^{\prime}) is the Green’s function obtained from Euler-Lagrange equations. Here η=θR\eta=\frac{\theta}{R} is compactified direction and R is the size of the extra compactified dimension. Using this, we calculate Casimir force and derive Casimir energy. Further, we study finite temperature corrections to the Casimir effect for DFRA real scalar field theory and obtain Casimir force in the high temperature and low-temperature limits. We also show that our results reduce to commutative space-time results in the appropriate limit of the size of the extra compactified dimension. We show that the Casimir force gets two types of corrections due to non-commutativity of space-time, one is LR\frac{L}{R} dependent modified Bessel function-K1K_{1} whoes coefficient vary as 1L\frac{1}{L} and the other is LR\frac{L}{R} dependent modified Bessel function-K2K_{2} whose coefficient vary as 1L2\frac{1}{L^{2}}. In the high-temperature limit, we find the correction terms that vary as 1L\frac{1}{L} and 1L3\frac{1}{L^{3}} while in the low-temperature limit Casimir force has correction terms that vary as 1L2\frac{1}{L^{2}} and 1L4\frac{1}{L^{4}}. We also show that the modifications to the Casimir force in both zero temperature and the finite temperature, has dependency on RR, size of the extra compactified dimension. To understand the generic effect of the non-commutativity on the Casimir force, we then consider most significant correction terms and analyzing the plots for the Casimir force with these correction terms, we obtain constraints on the size of the extra compactified dimension RR.

We also show that the massless DFRA scalar theory in 7-dimensions gives a massive scalar theory in 4-dimensions under Kaluza-Klein reduction. Also the modification of the Casimir effect using this model is calculated.

This paper is organized as follows. In the next section, we briefly discuss DFR space-time and underlying DFRA algebra [14, 15, 16, 17]. In Sec.3., we derive the general form of the Euler-Lagrange equation, current density, and the components of Energy-Momentum tensor corresponding to the scalar theory on the 77-dimensional DFR space-time. Our main results are obtained in Sec.4. Here we construct the Lagrangian describing massless DFRA scalar field theory in the presence of two parallel plates. We impose the Dirichlet boundary condition and calculate the Casimir force and the energy for massless scalar field with parallel plates in 4+14+1-dimensional DFR space-time. The θμν\theta^{\mu\nu} dependent terms in the Lagrangian captures the effects of non-commutativity of space-time. First, we study the Casimir effect by treating the new extra dimension-θ\theta introduced by non-commutativity in 4+14+1 dimensional DFR space-time, in same footing with other two transverse commutative space dimensions- yy and zz. Then by following the standard procedure, we show that the Casimir force expression is exactly same as the result obtained in case of commutative extra dimensional space-time. Then in Subsection.4.1, we study the Casimir effect by treating the additional dimension-θ\theta as a compactified dimension in 4+14+1 dimensional DFR space-time. This is carried out by first compactifying the extra dimension θ\theta of DFR space-time. Then we derive Green’s function for different regions of interest by solving the Euler-Lagrangian equation, describing the interaction of scalar field with parallel plates in presence of extra compactified dimension in 4+14+1-dimensional DFR space-time. Next we obtain the energy-momentum tensor corresponding to the massless DFRA real scalar theory from the general form of the stress-tensor constructed in Sec.3. Then we derive the vacuum expectation value of the energy-momentum tensor in terms of the obtained Green’s function. Using the reduced Green’s function obtained in different regions, we calculate the Casimir force and Casimir energy in presence of extra compactified dimension-θ\theta for 4+14+1 dimensional, massless DFRA scalar field theory. Further, we investigate finite temperature modifications to the Casimir effect in this model in Sec.5. We also study corrections to Casimir force for both the low-temperature limit and the high-temperature limit. Our concluding remarks are given in Sec.6. In the appendix, we discuss the study of the Casimir effect for massive scalar theory in 44-dimensional Minkowski space-time, obtained using Kaluza-Klein dimensional reduction from a 77-dimensional, massless, DFRA complex scalar theory. We use the metric, ηAB=(ημν,1)=diag(+1,1,1,1,1)\eta_{AB}=(\eta_{\mu\nu},-1)=diag(+1,-1,-1,-1,-1).

2 DFR space-time

In this section, we present a brief summary of the DFR space-time, a non-commutative space-time that respects Lorentz invariance and also the symmetry algebra associated to this space-time[4]. The Moyal-space-time coordinates satisfy

[x^μ,x^ν]=iθμν.[\hat{x}_{\mu},\hat{x}_{\nu}]=i{\theta}_{\mu\nu}. (2.1)

Since θμν\theta_{\mu\nu} is a constant antisymmetric tensor, it results in the violation of Lorentz symmetry [2]. Violation of Lorentz invariance brings issues with the interpretation of particle-antiparticle states. Assigning an appropriate transformation to θμν\theta_{\mu\nu} under Lorentz transformation reinstates the Lorentz invariance of theory [4] and thus avoids the issue of defining particle-antiparticle states. In [13], the components of the θμν\theta_{\mu\nu} were promoted to coordinate operators θ^μν\hat{\theta}_{\mu\nu}. Thus the non-commutative space-time gets six more dimensions and the resulting 10-dimensional space-time is known as DFR space-time. The corresponding algebra is given by

[x^μ,x^ν]=iθ^μν,[x^μ,θ^νρ]=0,[x^μ,p^ν]=iημν,[θ^μν,θ^ρλ]=0[p^μ,θ^νλ]=0,[p^μ,p^ν]=0\begin{split}[\hat{x}_{\mu},\hat{x}_{\nu}]&=i\hat{\theta}_{\mu\nu},~{}~{}[\hat{x}_{\mu},\hat{\theta}_{\nu\rho}]=0,\\ [\hat{x}_{\mu},\hat{p}_{\nu}]&=i\eta_{\mu\nu},~{}~{}[\hat{\theta}_{\mu\nu},\hat{\theta}_{\rho\lambda}]=0\\ [\hat{p}_{\mu},\hat{\theta}_{\nu\lambda}]&=0,~{}~{}~{}~{}~{}~{}[\hat{p}_{\mu},\hat{p}_{\nu}]=0\end{split} (2.2)

For Fourier decomposition of field, and also for the Lorentz invariance of the field theory defined on this non-commutative space-time, one needs to introduce conjugate momentum k^μν\hat{k}_{\mu\nu} - corresponding to θ^μν\hat{\theta}_{\mu\nu}, apart from x^μ\hat{x}^{\mu}, p^μ\hat{p}_{\mu} and θ^μν\hat{\theta}_{\mu\nu}. In DFR space-time [14, 15, 16], non-commutative coordinates θ^μν\hat{\theta}_{\mu\nu} and their canonical conjugate momenta k^μν\hat{k}_{\mu\nu} are considered as operators alongside with the commutative coordinates, x^μ\hat{x}_{\mu} and their canonical conjugate momenta p^μ\hat{p}_{\mu} . Thus algebra is enlarged by

[x^μ,k^νλ]=i2(ημνηρλημληνρ)p^ρ,[p^μ,k^νλ]=0,[θ^μν,k^ρλ]=i(ημρηνλημληνρ),[k^μν,k^ρλ]=0,\begin{split}[\hat{x}_{\mu},\hat{k}_{\nu\lambda}]&=-\frac{i}{2}(\eta_{\mu\nu}\eta_{\rho\lambda}-\eta_{\mu\lambda}\eta_{\nu\rho})\hat{p}^{\rho},~{}~{}[\hat{p}_{\mu},\hat{k}_{\nu\lambda}]=0,\\ [\hat{\theta}_{\mu\nu},\hat{k}_{\rho\lambda}]&=i(\eta_{\mu\rho}\eta_{\nu\lambda}-\eta_{\mu\lambda}\eta_{\nu\rho}),~{}~{}[\hat{k}_{\mu\nu},\hat{k}_{\rho\lambda}]=0,\end{split} (2.3)

in addition to the ones given in Eqn.(2.2). The relations in Eqn.(2.2) and Eqn.(2.3) together is a closed algebra-DFRA algebra[14, 15, 16, 17]. The corresponding Lorentz generator is defined as [14, 17]

Mμν=X^μp^νX^νp^μθ^μλk^ν+λθ^νλk^μλM_{\mu\nu}=\hat{X}_{\mu}\hat{p}_{\nu}-\hat{X}_{\nu}\hat{p}_{\mu}-\hat{\theta}_{\mu\lambda}{\hat{k}}_{\nu}~{}^{{\lambda}}+\hat{\theta}_{\nu\lambda}{\hat{k}}_{\mu}~{}^{\lambda} (2.4)

where,

X^μ=x^μ+12θ^μνp^ν.\hat{X}_{\mu}=\hat{x}_{\mu}+\frac{1}{2}\hat{\theta}_{\mu\nu}\hat{p}^{\nu}. (2.5)

The shifted coordinates in the above equation satisfy commutation relations

[X^μ,X^ν]=0,[X^μ,p^ν]=iημν.[\hat{X}_{\mu},\hat{X}_{\nu}]=0,~{}~{}~{}[\hat{X}_{\mu},\hat{p}_{\nu}]=i\eta_{\mu\nu}. (2.6)

With the relations given in Eqn.(2.2), Eqn.(2.3) and Eqn.(2.4), one sees that the generators of DFRA algebra are closed, i.e.,

[Mμν,p^λ]=i(ημλp^νηνλp^μ),[Mμν,k^αβ]=i(ημβk^ανημαk^νβ+ηναk^βμηνβk^αμ),[Mμν,Mλρ]=i(ημρMνληνρMλμημλMρν+ηνλMρμ),\begin{split}[M_{\mu\nu},\hat{p}_{\lambda}]&=i(\eta_{\mu\lambda}\hat{p}_{\nu}-\eta_{\nu\lambda}\hat{p}_{\mu}),\\ [M_{\mu\nu},\hat{k}_{\alpha\beta}]&=i(\eta_{\mu\beta}\hat{k}_{\alpha\nu}-\eta_{\mu\alpha}\hat{k}_{\nu\beta}+\eta_{\nu\alpha}\hat{k}_{\beta\mu}-\eta_{\nu\beta}\hat{k}_{\alpha\mu}),\\ [M_{\mu\nu},M_{\lambda\rho}]&=i(\eta_{\mu\rho}M_{\nu\lambda}-\eta_{\nu\rho}M_{\lambda\mu}-\eta_{\mu\lambda}M_{\rho\nu}+\eta_{\nu\lambda}M_{\rho\mu}),\end{split} (2.7)

The above relations of DFRA algebra are consistent with the Jacobi identities (see [14, 17]). The Casimir operator associated with the DFRA algebra is given as [14]

P^2=p^μp^μ+λ22k^μνk^μν,\hat{P}^{2}=\hat{p}_{\mu}\hat{p}^{\mu}+\frac{\lambda^{2}}{2}\hat{k}_{\mu\nu}\hat{k}^{\mu\nu}, (2.8)

where λ\lambda is the non-commutative parameter with the length dimension. In [14, 15], the construction of infinitesimal transformations of xμ,θμν,pμ,kμνx_{\mu},~{}\theta_{\mu\nu},~{}p_{\mu},~{}k_{\mu\nu} and MμνM_{\mu\nu} and their relevance for the DFR space-time are discussed in detail.

3 Construction of general form of Energy-Momentum tensor and Current density

One way to evaluate the Casimir force experienced by a plate is to obtain the vacuum expectation value of the Energy-Momentum tensor on either side of this plate and find their difference. Since the Energy-Momentum tensor can be written as a differential operator acting on the product of fields, this vacuum expectation value of the Energy-Momentum tensor is related to the vacuum expectation value of the time-order product of fields at two space-time points acted upon by this operator. This allows one to express vacuum expectation value of Energy-Momentum tensor as the operator acting on Green’s function, i.e.,

<Tx,xμν>\displaystyle<T_{x,x^{\prime}}^{\mu\nu}> =\displaystyle= O^<T(ϕ(x)ϕ(x))>\displaystyle\hat{O}<T(\phi(x)\phi(x^{\prime}))> (3.1)
=\displaystyle= O^G(x,x)\displaystyle\hat{O}G(x,x^{\prime})

The Green’s function on either side of the plate are calculated by solving the Euler-Lagrange equation for these different regions. In this section, we summarize the construction of action for the scalar field theory in DFR space-time using the definition of the Moyal star product [2]. By varying the action, we then obtain the equations of motion corresponding to the scalar field defined in the DFR space-time. Here we also discuss the implications of the θ\theta-dependent weight function introduced in the definition of action.

3.1 Construction of Action

Before construction of the action for the DFRA scalar field, we start with the definition of star product [4], which for two functions ff and gg of xμx_{\mu} and θμν\theta_{\mu\nu}, given as

f(x,θ)g(x,θ)=ei2θμνμνf(x,θ)g(x,θ)|x=x.f(x,\theta)\star g(x,\theta)=e^{\frac{i}{2}\theta^{\mu\nu}\partial_{\mu}\partial_{\nu}^{\prime}}f(x,\theta)g(x^{\prime},\theta)\Big{|}_{x=x^{\prime}}. (3.2)

Using above definition, it can be shown that the Moyal product satisfy following relation,

d4xd6θW(θ)f(x,θ)g(x,θ)=d4xd6θW(θ)f(x,θ)g(x,θ).\int d^{4}x~{}d^{6}\theta~{}W(\theta)~{}f(x,\theta)\star g(x,\theta)=\int_{-\infty}^{\infty}d^{4}x~{}d^{6}\theta~{}W(\theta)~{}f(x,\theta)g(x,\theta). (3.3)

W(θ)W(\theta) in the above equation is the weight function and has been introduced to avoid divergence in perturbative QFT on DFR space-time [46, 13, 47, 48]. For Lorentz invariance, W(θ)W(\theta) is taken as even function of θ\theta and properties of the weight function is discussed in [46, 13, 47, 48] in detail. Here we consider weight function to be Gaussian function, given by

W(θ)=(14π2λ4)3eθ28λ4,W(\theta)=\Big{(}\frac{1}{4\pi^{2}\lambda^{4}}\Big{)}^{3}e^{-\frac{\theta^{2}}{8\lambda^{4}}}, (3.4)

where λ\lambda is the non-commutative parameter of length dimension.

Now one sets up the action for the scalar field theory in DFR space-time by replacing the usual product in the Lagrangian with the star product. As a result, the non-commutative action for the DFRA scalar field is

S=d4xd6θW(θ)12(μϕμϕ+λ22θμνϕθμνϕm2ϕϕ).S=\int d^{4}x~{}d^{6}\theta~{}W(\theta)\frac{1}{2}\Big{(}\partial_{\mu}\phi\star\partial^{\mu}\phi+\frac{\lambda^{2}}{2}\partial_{\theta^{\mu\nu}}\phi\star\partial^{\theta^{\mu\nu}}\phi-m^{2}\phi\star\phi\Big{)}. (3.5)

After using the property given in Eqn.(3.3), the action for DFRA scalar field becomes

S=d4xd6θW(θ)12(μϕμϕ+λ22θμνϕθμνϕm2ϕ2).\displaystyle S=\int d^{4}x~{}d^{6}\theta~{}W(\theta)\frac{1}{2}\Big{(}\partial_{\mu}\phi\partial^{\mu}\phi+\frac{\lambda^{2}}{2}\partial_{\theta^{\mu\nu}}\phi\partial^{\theta^{\mu\nu}}\phi-m^{2}\phi^{2}\Big{)}. (3.6)

For constant value of weight function (specifically for W(θ)=1W(\theta)=1) above action reduces to that studied in [17].

In order to preserve the unitarity of the field theory on DFR space-time, it has been shown that one needs to set temporal part of non-commutative coordinates to be vanishing, i.e., θ^0i=0\hat{\theta}_{0i}=0 [43] and thus the resulting DFR space-time has seven dimensions. Further we use definitions, θi=12ϵijlθjl\theta^{i}=\frac{1}{2}\epsilon^{ijl}\theta_{jl} and ki=12ϵijlkjlk^{i}=\frac{1}{2}\epsilon^{ijl}k_{jl} [45] and thus from above we find 4+34+3 dimensional action for real scalar field in DFR space-time as

S\displaystyle S =\displaystyle= d4xd3θW(θ)(x,θ,ϕ,μϕ,θiϕ)\displaystyle\int d^{4}x~{}d^{3}\theta~{}W(\theta){\cal L}(x,\theta,\phi,\partial_{\mu}\phi,\partial_{\theta^{i}}\phi) (3.7)
=\displaystyle= d4xd3θW(θ)12(μϕμϕ+λ2θiϕθiϕm2ϕ2).\displaystyle\int d^{4}x~{}d^{3}\theta~{}W(\theta)\frac{1}{2}\Big{(}\partial_{\mu}\phi\partial^{\mu}\phi+\lambda^{2}\partial_{\theta^{i}}\phi\partial^{\theta^{i}}\phi-m^{2}\phi^{2}\Big{)}. (3.8)

Note here that the W(θ)W(\theta) in the measure and λ\lambda dependent term in Lagrangian bring in the non-commutative contribution due to the DFR space-time. We get the commutative action of the scalar field through compactification i.e., for limλ0W(θ)=limλ0eθ24λ42πλ2=δ(θ)\lim_{\lambda\to 0}W(\theta)=\lim_{\lambda\to 0}\frac{e^{-\frac{\theta^{2}}{4\lambda^{4}}}}{2\pi\lambda^{2}}=\delta(\theta) and 𝑑θW(θ)=1\int d\theta~{}W(\theta)=1. Note that in [17], W(θ)W(\theta) is taken to be equal to 1.

3.2 Construction of Energy-Momentum tensor and Current density for DFRA scalar field

Here in this subsection, we obtain the expression of Energy-Momentum tensor and the current density for DFRA scalar field theory by varying the action by an infinitesimal change in the space-time coordinates of DFR space-time, (xμ,θμν)(x^{\mu},\theta^{\mu\nu}).

We start with the general form of action (3.7),

S=Rd4xd3θW(θ)(x,θ,ϕ,μϕ,θiϕ)S=\int_{R}d^{4}x~{}d^{3}\theta~{}W(\theta){\cal L}(x,\theta,\phi,\partial_{\mu}\phi,\partial_{\theta^{i}}\phi) (3.9)

where R is a bounded region. Next we vary the action by considering functional change of field as well as infinitesimal change in space-time co-ordinates such as xμ=xμ+ϵδxμ{x^{\mu}}^{\prime}=x^{\mu}+\epsilon\delta x^{\mu} and θi=θi+ϵδθi{\theta^{i}}^{\prime}=\theta^{i}+\epsilon\delta{\theta^{i}} and thus we get

δS\displaystyle\delta S =\displaystyle= Rd4xd3θW(θ)(ϕ,μϕ,θiϕ)Rd4xd3θW(θ)(ϕ,μϕ,θiϕ)\displaystyle\int_{R}d^{4}x^{\prime}d^{3}\theta^{\prime}W(\theta^{\prime})\mathcal{L^{\prime}}(\phi^{\prime},\partial_{\mu}\phi^{\prime},\partial_{\theta^{i}}\phi^{\prime})-\int_{R}d^{4}xd^{3}\theta W(\theta^{\prime})\mathcal{L}(\phi,\partial_{\mu}\phi,\partial_{\theta^{i}}\phi) (3.10)
=\displaystyle= Rd4xd3θW(θ)δ+ϵRd4xd3θW(θ)(μδxμ)+ϵRd4xd3θθi(W(θ)δθi)\displaystyle\int_{R}d^{4}xd^{3}\theta W(\theta)\delta\mathcal{L}+\epsilon\int_{R}d^{4}xd^{3}\theta W(\theta)(\partial_{\mu}\delta x^{\mu})\mathcal{L}+\epsilon\int_{R}d^{4}xd^{3}\theta\partial_{\theta^{i}}(W(\theta)\delta\theta^{i})\mathcal{L} (3.11)

Next we use the functional variance of field δϕ(x,θ)=δ¯ϕ(x,θ)+δxμ(μϕ)+δθi(θiϕ)\delta\phi(x,\theta)=\bar{\delta}\phi(x,\theta)+\delta x^{\mu}(\partial_{\mu}\phi)+\delta\theta^{i}(\partial_{\theta^{i}}\phi) and get,

δS=d4xd3θ{WLϕWμ(L(μϕ))θi(WL(θiϕ))}δ¯ϕ+d4xd3θ[μ(Wδxμ+W(μϕ)δ¯ϕ)+θi(Wδθi+W(θiϕ)δϕ¯)]\delta S=\int d^{4}xd^{3}\theta\Big{\{}W\frac{\partial L}{\partial\phi}-W\partial_{\mu}\Big{(}\frac{\partial L}{\partial(\partial_{\mu}\phi)}\Big{)}-\partial_{\theta^{i}}\Big{(}W\frac{\partial L}{\partial(\partial_{\theta^{i}}\phi)}\Big{)}\Big{\}}\bar{\delta}\phi\\ +\int d^{4}xd^{3}\theta\Big{[}\partial_{\mu}(W\mathcal{L}\delta x^{\mu}+W\frac{\partial\mathcal{L}}{\partial(\partial_{\mu}\phi)}\bar{\delta}\phi)+\partial_{\theta^{i}}(W\mathcal{L}\delta\theta^{i}+W\frac{\partial\mathcal{L}}{\partial(\partial_{\theta^{i}}\phi)}\bar{\delta\phi})\Big{]} (3.12)

where δ¯ϕ=ϕ(x,θ)ϕ(x,θ)\bar{\delta}\phi=\phi^{\prime}(x,\theta)-\phi(x,\theta). After using Gauss divergence theorem and setting total derivative terms to vanish at the boundary R\partial R, also by assuming a stationary value for S for an arbitrary variation of δ¯ϕ\bar{\delta}\phi that vanishes on the boundary R\partial R, we get Euler-Lagrange equation as,

Wμ((μϕ))Wϕ+θi(W(θiϕ))=0.W\partial_{\mu}\Big{(}\frac{\partial\mathcal{L}}{\partial(\partial_{\mu}\phi)}\Big{)}-W\frac{\partial\mathcal{L}}{\partial\phi}+\partial_{\theta^{i}}\Big{(}W\frac{\partial\mathcal{L}}{\partial(\partial_{\theta^{i}}\phi)}\Big{)}=0. (3.13)

From surface terms in Eqn.(3.12), using Noether’s theorem, we get conserved current (Jμ,Jθi)(J^{\mu},J^{\theta^{i}}) as,

Jμ=W{δxμ+((μϕ))δ¯ϕ}=WJ~μJ^{\mu}=W\Big{\{}\mathcal{L}\delta x^{\mu}+\Big{(}\frac{\partial\mathcal{L}}{\partial(\partial_{\mu}\phi)}\Big{)}\bar{\delta}\phi\Big{\}}=W\tilde{J}^{\mu} (3.14)

and

Jθi=W{δθi+(θiϕ)δ¯ϕ)}=WJ~θiJ^{\theta^{i}}=W\Big{\{}\mathcal{L}\delta\theta^{i}+\frac{\partial\mathcal{L}}{\partial(\partial_{\theta^{i}}\phi)}\bar{\delta}\phi)\Big{\}}=W\tilde{J}^{\theta^{i}} (3.15)

such that

Ω=μ(WJ~μ)+θi(WJ~θi)=0.\Omega=\partial_{\mu}(W\tilde{J}^{\mu})+\partial_{\theta^{i}}(W\tilde{J}^{\theta^{i}})=0. (3.16)

Now using the relation δϕ(x,θ)=δ¯ϕ(x,θ)+δxμ(μϕ)+δθi(θiϕ)\delta\phi(x,\theta)=\bar{\delta}\phi(x,\theta)+\delta x^{\mu}(\partial_{\mu}\phi)+\delta\theta^{i}(\partial_{\theta^{i}}\phi) in Eqn.(3.14) and Eqn.(3.15) eventually we get the components of the Energy-Momentum tensor as,

Tμν\displaystyle T^{\mu\nu} =\displaystyle= ((μϕ))νϕημν\displaystyle\Big{(}\frac{\partial\mathcal{L}}{\partial(\partial_{\mu}\phi)}\Big{)}\partial^{\nu}\phi-\eta^{\mu\nu}\mathcal{L} (3.17)
Tθiθj\displaystyle T^{\theta^{i}\theta^{j}} =\displaystyle= ((θiϕ))θjϕgθiθj\displaystyle\Big{(}\frac{\partial\mathcal{L}}{\partial(\partial_{\theta^{i}}\phi)}\Big{)}\partial^{\theta^{j}}\phi-g^{\theta^{i}\theta^{j}}\mathcal{L} (3.18)
Tμθi\displaystyle T^{\mu\theta^{i}} =\displaystyle= ((μϕ))θiϕ\displaystyle\Big{(}\frac{\partial\mathcal{L}}{\partial(\partial_{\mu}\phi)}\Big{)}\partial^{\theta^{i}}\phi (3.19)
Tθiμ\displaystyle T^{\theta^{i}\mu} =\displaystyle= ((θiϕ))μϕ\displaystyle\Big{(}\frac{\partial\mathcal{L}}{\partial(\partial_{\theta^{i}}\phi)}\Big{)}\partial^{\mu}\phi (3.20)

Thus we get all the components of stress tensor for real scalar field theory in 77-dimensional DFR space-time. Note that the Casimir effect will be calculated by taking the difference between the vacuum expectation value of the Energy-Momentum tensor at the right side of the plate and the left side of the plate. Also, note here that current densities and components of the stress tensor go over to the commutative results when θ=0\theta=0. For constant value of weight function, (say W(θ)=1W(\theta)=1) the Euler-Lagrange equation of motion in Eqn.(3.13), current densities in Eqn.(3.14) and Eqn.(3.15), reduces to the results for DFRA real scalar field theory considered in [17].

4 Casimir Effect using DFRA scalar field

In this section, we derive the Casimir force felt by one of the parallel plates by analyzing the vacuum fluctuation of the DFRA scalar field. For this, we first calculate the vacuum expectation value of the Energy-Momentum tensor corresponding to the DFRA scalar field theory, using the relation between the expectation value of the time-ordered product of fields and Green’s function. Here the Green’s function is obtained by solving Euler-Lagrange equations. Using this vacuum expectation value of Energy-Momentum tensor, we obtain corrections to Casimir force and Casimir energy. We investigate Casimir effect in 4+14+1 dimensions (xμ,θ)(x^{\mu},\theta) and the result we obtain is generalized to 4+34+3-dimensional DFR space-time. Here the metric used is diag(+)diag(+------).

We start with the DFRA massless real scalar field Lagrangian [14, 15, 16, 17],

0=12μϕμϕ+λ22θiϕθiϕ\mathcal{L}_{0}=\frac{1}{2}\partial_{\mu}\phi\partial^{\mu}\phi+\frac{\lambda^{2}}{2}\partial_{\theta^{i}}\phi\partial^{\theta^{i}}\phi (4.1)

where in the limit λ0\lambda\rightarrow 0, we get back Lagrangian of the real scalar field theory in commutative space-time.

To study Casimir effect between two parallel plates, we introduces these plates perpendicular to xx-direction and model these plates, using delta function potentials at x=0x=0 and x=Lx=L. Thus we get total Lagrangian as [19],

=12μϕμϕ+λ22θiϕθiϕγ2Lϕ2δ(x)γ2Lϕ2δ(xL),\mathcal{L}=\frac{1}{2}\partial_{\mu}\phi\partial^{\mu}\phi+\frac{\lambda^{2}}{2}\partial_{\theta^{i}}\phi\partial^{\theta^{i}}\phi-{\frac{\gamma}{2L}\phi^{2}\delta(x)}-{\frac{\gamma^{\prime}}{2L}}\phi^{2}\delta(x-L), (4.2)

where γ,γ\gamma,\gamma^{\prime} are coupling constant, which are dimension free quantities. Here we consider the Lagrangian in 4+14+1 dimension given by

=120ϕ0ϕ12mϕmϕ+λ22θ1ϕθ1ϕγ2Lϕ2δ(x)γ2Lϕ2δ(xL);m=1,2,3.\mathcal{L}=\frac{1}{2}{\partial_{0}\phi\partial_{0}\phi}-\frac{1}{2}{\partial_{m}\phi\partial_{m}\phi}+\frac{\lambda^{2}}{2}\partial_{\theta^{1}}\phi\partial^{\theta^{1}}\phi-{\frac{\gamma}{2L}\phi^{2}\delta(x)}-{\frac{\gamma^{\prime}}{2L}}\phi^{2}\delta(x-L);~{}~{}m=1,2,3. (4.3)

Note that the interaction part, mimicking the plates, is similar to that in the commutative space-time.

Using Eqn.(3.13), we obtain Euler-Lagrange equation from the above Lagrangian as

(02m2λ2θ12+θ12λ2θ1+γLδ(x)+γLδ(xL))ϕ(𝐱,θ1λ)=0,\Big{(}\partial_{0}^{2}-\partial_{m}^{2}-\lambda^{2}\partial_{\theta^{1}}^{2}+\frac{\theta^{1}}{2\lambda^{2}}\partial_{\theta^{1}}+\frac{\gamma}{L}\delta(x)+\frac{\gamma^{\prime}}{L}\delta(x-L)\Big{)}\phi({\bf x},\frac{\theta^{1}}{\lambda})=0, (4.4)

where 𝐱{\bf x} denote all the three space co-ordinates x,y,zx,y,z and θ1\theta^{1} is the non-commutative of space coordinate. Note that if weight function W(θ)W(\theta) is taken to be a constant, the third term, θ12λ2θ1\frac{\theta^{1}}{2\lambda^{2}}\partial_{\theta^{1}} will be absent. We re-express Eqn.(4.4) as,

(02x2y2z2λ2θ12+θ12λ2θ1+γLδ(x)+γLδ(xL))ϕ(𝐱,θ1λ)=0.\Big{(}\partial_{0}^{2}-\partial_{x}^{2}-\partial_{y}^{2}-\partial_{z}^{2}-\lambda^{2}\partial_{\theta^{1}}^{2}+\frac{\theta^{1}}{2\lambda^{2}}\partial_{\theta^{1}}+\frac{\gamma}{L}\delta(x)+\frac{\gamma^{\prime}}{L}\delta(x-L)\Big{)}\phi({\bf x},\frac{\theta^{1}}{\lambda})=0. (4.5)

Here in the limit, λ0\lambda\rightarrow 0 and θ=0\theta=0, we get back equation of motion in the commutative space-time [19]. As we are considering only one non-commutative direction θ1\theta^{1}, we set θ1=θ\theta^{1}=\theta.

Casimir force between two parallel plates is obtained by taking the difference between vacuum expectation value of xxxx component of Energy-Momentum tensor at either side of the plate. Since the vacuum expectation values of stress tensor is related to the Green’s function as in Eqn.(3.1), first we calculate the Green’s functions corresponding to Eqn.(4.5). We start with,

(02x2y2z2λ2θ2+θ2λ2θ+γLδ(x)+γLδ(xL))G(𝐱,𝐱;θλ,θλ)=δ3(𝐱𝐱)δ(tt)δ(θλθλ).\Big{(}\partial_{0}^{2}-\partial_{x}^{2}-\partial_{y}^{2}-\partial_{z}^{2}-\lambda^{2}\partial_{\theta}^{2}+\frac{\theta}{2\lambda^{2}}\partial_{\theta}+\frac{\gamma}{L}\delta(x)+\frac{\gamma^{\prime}}{L}\delta(x-L)\Big{)}G({\bf x,x^{\prime}};\frac{\theta}{\lambda},\frac{\theta^{\prime}}{\lambda})=\delta^{3}({\bf x-x^{\prime}})\delta(t-t^{\prime})\delta(\frac{\theta}{\lambda}-\frac{\theta^{\prime}}{\lambda}). (4.6)

Note that in the commutative limit, i.e., λ0\lambda\rightarrow 0, θ=0\theta=0 above equation reduce to well known commutative result [19]. In the commutative space-time one solves for the Green’s function by using Fourier-transform for tt, yy and zz dependence of the Green function. We can now treat the additional θ\theta-direction introduced by in the DFR space-time in the same manner as yy and zz directions. Then following the standard calculation we would get the Casimir force to be F=3.1232π2L5F=-\frac{3.12}{32\pi^{2}L^{5}} (for weight function W(θ)=1W(\theta)=1) which exactly match with the result obtained in [19] for 5-dimensional commutative space-time. This shows that the newly introduced θ\theta-dimension act just like extra dimensions in commutative space-time. Note that this expression for the Casimir force is independent of θ\theta and/or λ\lambda and the Casimir force in DFR space-time does not have a commuatative limit; thus this Casimir effect does not distinguise between 5-dimensional DFR space-time and 5-dimensional commutative space-time. Since non-commutativity is expected to bring in modifications which should smoothly vanish in the commutative limit, we adopt a different calculation scheme in the following. In [20, 38, 39, 40, 41] extra dimensions were treated as compact directions and the Casimir effect was studied. We next study the Casimir effect in DFR space-time by treating additional dimensions introduced by non-commutativity as compact dimensions.

4.1 Green’s function solutions in different regions

In this subsection, to study the Casimir effect in presence of extra compactified dimension, we compactify the extra θ\theta - direction which comes due to the presence of non-commutativity of space-time. We then find the Green’s function solutions in different regions in presence of this extra compactified dimension. Using Fourier transform, we rewrite Green’s function in Eqn.(4.6) as,

G(𝐱,𝐱,t,t,θ~,θ~)=dw2πdpy2πdpz2πeiω(tt)eipy(yy)eipz(zz)g(x,x,θ~,θ~,ω,p).G({\bf x,x^{\prime}},t,t^{\prime},\tilde{\theta},\tilde{\theta}^{\prime})=\int_{-\infty}^{\infty}\frac{dw}{2\pi}\frac{dp_{y}}{2\pi}\frac{dp_{z}}{2\pi}e^{-i\omega(t-t^{\prime})}e^{ip_{y}(y-y^{\prime})}e^{ip_{z}(z-z^{\prime})}g(x,x^{\prime},\tilde{\theta},\tilde{\theta^{\prime}},\omega,p_{\|}). (4.7)

where θλ=θ~\frac{\theta}{\lambda}=\tilde{\theta}, which has dimension of length, p=py2+pz2p_{\|}=\sqrt{p_{y}^{2}+p_{z}^{2}} and g(x,x,θ~,θ~,ω,p)g(x,x^{\prime},\tilde{\theta},\tilde{\theta^{\prime}},\omega,p_{\|}) is the reduced Green’s function. Note here that we are not considering non-commutative coordinate θ\theta in same footing as the other two transverse space dimensions, yy and zz but treats θ\theta as a compactified direction. As discussed after Eqn.(4.6), accounting θ\theta coordinate in a way similar that of the transverse dimensions, leads to results which does not have proper commutative limit. To overcome this shortcoming here we will compactify this extra space dimension θ\theta and study Casimir effect due to the presence of this extra compactified space dimension [38, 39, 40]. The g(x,x,θ~,θ~,ω,p)g(x,x^{\prime},\tilde{\theta},\tilde{\theta^{\prime}},\omega,p_{\|}) obeys the equation,

(2x2ω2+p2θ~2+θ~2λ2θ~+γLδ(x)+γLδ(xL))g(x,x;θ~,θ~;ω;p)=δ(xx)δ(θ~θ~).\Big{(}{-\frac{\partial^{2}}{\partial x^{2}}}-\omega^{2}+p_{\|}^{2}-\frac{\partial}{\partial\tilde{\theta}^{2}}+\frac{\tilde{\theta}}{2\lambda^{2}}\frac{\partial}{\partial{\tilde{\theta}}}+\frac{\gamma}{L}\delta(x)+\frac{\gamma^{\prime}}{L}\delta(x-L)\Big{)}g(x,x^{\prime};\tilde{\theta},\tilde{\theta};\omega;p_{\|})=\delta(x-x^{\prime})\delta(\tilde{\theta}-\tilde{\theta}^{\prime}). (4.8)

where p2=py2+pz2p_{\|}^{2}=p_{y}^{2}+p_{z}^{2} and second term, ω2-\omega^{2} in above is coming due to action of 02\partial_{0}^{2} over Green’s function given in Eqn.(4.7). Next we compactify θ~\tilde{\theta} direction by the definition θ~R=η\frac{\tilde{\theta}}{R}=\eta. Here RR is the size of the extra compactified dimension. Thus we use η\eta as the extra compactified direction, which varies from 02π0\rightarrow 2\pi and we rewrite Eqn.(4.8) as

(2x2ω2+p21R2η2+η2λ2η+γLδ(x)+γLδ(xL))g(x,x;η,η;ω;p)=δ(xx)nδ(ηηn).\Big{(}{-\frac{\partial^{2}}{\partial x^{2}}}-\omega^{2}+p_{\|}^{2}-\frac{1}{R^{2}}\frac{\partial}{\partial\eta^{2}}+\frac{\eta}{2\lambda^{2}}\frac{\partial}{\partial{\eta}}+\frac{\gamma}{L}\delta(x)+\frac{\gamma^{\prime}}{L}\delta(x-L)\Big{)}g(x,x^{\prime};\eta,\eta^{\prime};\omega;p_{\|})=\delta(x-x^{\prime})\sum_{n\in\mathbb{Z}}\delta(\eta-\eta^{\prime}-n). (4.9)

In above equation the reduced green’s function is

g(x,x;η,η;ω;p)=ng(x,x;ω;p,n)ei2πn(ηη)g(x,x^{\prime};\eta,\eta^{\prime};\omega;p_{\|})=\sum_{n\in\mathbb{Z}}g(x,x^{\prime};\omega;p_{\|},n)e^{i2\pi n(\eta-\eta^{\prime})} (4.10)

and we use Poisson summation formula, nδ(ηηn)=nei2πn(ηη)\sum_{n\in\mathbb{Z}}\delta(\eta-\eta^{\prime}-n)=\sum_{n\in\mathbb{Z}}e^{i2\pi n(\eta-\eta^{\prime})}. Thus we reexpress Eqn.(4.9) as

n((2x2ω2+p2+(2πn)2R2+η2λ2(i2πn)+γLδ(x)+γLδ(xL))g(x,x;ω;p;n)δ(xx))ei2πn(ηη)=0.\sum_{n\in\mathbb{Z}}\bigg{(}\Big{(}{-\frac{\partial^{2}}{\partial x^{2}}}-\omega^{2}+p_{\|}^{2}+\frac{(2\pi n)^{2}}{R^{2}}+\frac{\eta}{2\lambda^{2}}(i2\pi n)+\frac{\gamma}{L}\delta(x)+\frac{\gamma^{\prime}}{L}\delta(x-L)\Big{)}g(x,x^{\prime};\omega;p_{\|};n)-\delta(x-x^{\prime})\bigg{)}e^{i2\pi n(\eta-\eta^{\prime})}=0. (4.11)

from which we find

(2x2ω2+p2+Cn2+γLδ(x)+γLδ(xL))g(x,x;ω;p;n)=δ(xx),\Big{(}{-\frac{\partial^{2}}{\partial x^{2}}}-\omega^{2}+p_{\|}^{2}+C_{n}^{2}+\frac{\gamma}{L}\delta(x)+\frac{\gamma^{\prime}}{L}\delta(x-L)\Big{)}g(x,x^{\prime};\omega;p_{\|};n)=\delta(x-x^{\prime}), (4.12)

where Cn2=(2πn)2R2+η2λ2(i2πn)C_{n}^{2}=\frac{(2\pi n)^{2}}{R^{2}}+\frac{\eta}{2\lambda^{2}}(i2\pi n), which is a complex quantity. Here nn\in\mathbb{Z} and η=θ~R\eta=\frac{\tilde{\theta}}{R}. Next, we solve Eqn.(4.12) using Dritchlet boundary condition and obtain the reduced Green’s function in three regions as follows.


For region x,x<0x,x^{\prime}<0,

g=12qeq|xx|+12qΔeq(x+x)(γ2qL(1γ2qL)γ2qL(1+γ2qL)e2qL)+n\{0}[12q~eq~|xx|+12q~Δ~eq~(x+x)(γ2q~L(1γ2q~L)γ2q~L(1+γ2q~L)e2q~L)]ei2πn(ηη);g=\frac{1}{2q}e^{-q|x-x^{\prime}|}+{\frac{1}{2q\Delta}e^{q(x+x^{\prime})}}\bigg{(}{-{\frac{\gamma^{\prime}}{2qL}\Big{(}1-\frac{\gamma}{2qL}\Big{)}}{-{\frac{\gamma}{2qL}}\Big{(}1+\frac{\gamma^{\prime}}{2qL}\Big{)}e^{2qL}}}\bigg{)}\\ +\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\Big{[}\frac{1}{2\tilde{q}}e^{-\tilde{q}|x-x^{\prime}|}+{\frac{1}{2\tilde{q}\tilde{\Delta}}e^{\tilde{q}(x+x^{\prime})}}\bigg{(}{-{\frac{\gamma^{\prime}}{2\tilde{q}L}\Big{(}1-\frac{\gamma}{2\tilde{q}L}\Big{)}}{-{\frac{\gamma}{2\tilde{q}L}}\Big{(}1+\frac{\gamma^{\prime}}{2\tilde{q}L}\Big{)}e^{2\tilde{q}L}}}\bigg{)}\Big{]}e^{i2\pi n(\eta-\eta^{\prime})}; (4.13)

for region 0<x,x<L0<x,x^{\prime}<L,

g=12qeq|xx|+12qΔ(γγ(2qL)22cosh(q|xx|)γ2qL(1+γ2qL)e2qLeq(x+x)γ2qL(1+γ2qL)eq(x+x))+n\{0}[12q~eq~|xx|+12q~Δ~(γγ(2q~L)22cosh(|q~xx|)γ2q~L(1+γ2q~L)e2q~Leq~(x+x)γ2q~L(1+γ2q~L)eq~(x+x))]ei2πn(ηη);g=\frac{1}{2q}e^{-q|x-x^{\prime}|}+\frac{1}{2q\Delta}\bigg{(}\frac{\gamma\gamma^{\prime}}{(2qL)^{2}}2\textnormal{cosh}(q|x-x^{\prime}|)-{\frac{\gamma}{2qL}\Big{(}1+\frac{\gamma^{\prime}}{2qL}\Big{)}e^{2qL}e^{-q(x+x^{\prime})}}-{\frac{\gamma^{\prime}}{2qL}\Big{(}1+\frac{\gamma}{2qL}\Big{)}e^{q(x+x^{\prime})}}\bigg{)}\\ +\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\Big{[}\frac{1}{2\tilde{q}}e^{-\tilde{q}|x-x^{\prime}|}+\frac{1}{2\tilde{q}\tilde{\Delta}}\bigg{(}\frac{\gamma\gamma^{\prime}}{(2\tilde{q}L)^{2}}2\textnormal{cosh}(|\tilde{q}x-x^{\prime}|)-{\frac{\gamma}{2\tilde{q}L}\Big{(}1+\frac{\gamma^{\prime}}{2\tilde{q}L}\Big{)}e^{2\tilde{q}L}e^{-\tilde{q}(x+x^{\prime})}}\\ -{\frac{\gamma^{\prime}}{2\tilde{q}L}\Big{(}1+\frac{\gamma}{2\tilde{q}L}\Big{)}e^{\tilde{q}(x+x^{\prime})}}\bigg{)}\Big{]}e^{i2\pi n(\eta-\eta^{\prime})}; (4.14)

for region L<x,xL<x,x^{\prime},

g=12qeq|xx|+12qΔeq(x+x2L)(γ2qL(1γ2qL)γ2qL(1+γ2qL)e2qL)+n\{0}[12q~eq~|xx|+12q~Δ~eq~(x+x2L)(γ2q~L(1γ2q~L)γ2q~L(1+γ2q~L)e2q~L)]ei2πn(ηη);g=\frac{1}{2q}e^{-q|x-x^{\prime}|}+{\frac{1}{2q\Delta}e^{-q(x+x^{\prime}-2L)}}\bigg{(}{-{\frac{\gamma}{2qL}\Big{(}1-\frac{\gamma^{\prime}}{2qL}\Big{)}}}{-{\frac{\gamma^{\prime}}{2qL}\Big{(}1+\frac{\gamma}{2qL}}\Big{)}e^{2qL}}\bigg{)}\\ +\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\Big{[}\frac{1}{2\tilde{q}}e^{-\tilde{q}|x-x^{\prime}|}+{\frac{1}{2\tilde{q}\tilde{\Delta}}e^{-\tilde{q}(x+x^{\prime}-2L)}}\bigg{(}{-{\frac{\gamma}{2\tilde{q}L}\Big{(}1-\frac{\gamma^{\prime}}{2\tilde{q}L}\Big{)}}}{-{\frac{\gamma^{\prime}}{2\tilde{q}L}\Big{(}1+\frac{\gamma}{2\tilde{q}L}}\Big{)}e^{2\tilde{q}L}}\bigg{)}\Big{]}e^{i2\pi n(\eta-\eta^{\prime})}; (4.15)

Here, n\{0}{n\in\mathbb{Z}\backslash\left\{0\right\}} imply that nn can take any integer value except 0. In the above equations we have defined q=ω2+p2q=-\omega^{2}+p_{\|}^{2} and q~=Cn2ω2+p2\tilde{q}=\sqrt{C_{n}^{2}-\omega^{2}+p_{\|}^{2}},where Cn2=(2πn)2R2+η2λ2(i2πn)C_{n}^{2}=\frac{(2\pi n)^{2}}{R^{2}}+\frac{\eta}{2\lambda^{2}}(i2\pi n). And we have also used

Δ=(1+γ2qL)(1+γ2qL)e2qLγγ(2qL)2\Delta=\Big{(}1+\frac{\gamma}{2qL}\Big{)}\Big{(}1+\frac{\gamma^{\prime}}{2qL}\Big{)}e^{2qL}-\frac{\gamma\gamma^{\prime}}{(2qL)^{2}} (4.16)

and

Δ~=(1+γ2q~L)(1+γ2q~L)e2q~Lγγ(2q~L)2.\tilde{\Delta}=\Big{(}1+\frac{\gamma}{2\tilde{q}L}\Big{)}\Big{(}1+\frac{\gamma^{\prime}}{2\tilde{q}L}\Big{)}e^{2\tilde{q}L}-\frac{\gamma\gamma^{\prime}}{(2\tilde{q}L)^{2}}. (4.17)

Note here that in the limit R0R\rightarrow 0 (size of the compactified dimension), q~\tilde{q}\rightarrow\infty and we get commutative results of the reduced Green’s function solutions [19] for all the three regions mentioned above.

4.2 Modification of Energy-Momentum tensor

In this subsection, we obtain the stress tensor for DFRA scalar field theory in 4+14+1 dimensions using the Lagrangian given in Eqn.(4.3) in Eqn.(3.17). Thus we find the stress tensor as

Tμν\displaystyle T^{\mu\nu} =\displaystyle= μϕνϕημν\displaystyle\partial^{\mu}\phi\partial^{\nu}\phi-\eta^{\mu\nu}\mathcal{L} (4.18)
=\displaystyle= μϕνϕημν12(0ϕ0ϕmϕmϕ+θ~ϕθ~ϕγLδ(x)γLδ(xL))\displaystyle\partial^{\mu}\phi\partial^{\nu}\phi-\eta^{\mu\nu}\frac{1}{2}\Big{(}\partial^{0}\phi\partial^{0}\phi-\partial^{m}\phi\partial^{m}\phi+\partial_{\tilde{\theta}}\phi\partial^{\tilde{\theta}}\phi-\frac{\gamma}{L}\delta(x)-\frac{\gamma^{\prime}}{L}\delta(x-L)\Big{)}

In the limit λ0\lambda\rightarrow 0 and θ=0\theta=0, above TμνT^{\mu\nu} reduces to that of the real scalar field in the commutative space-time. We re-express Energy-Momentum tensor as

T^𝐱,𝐱,θ,θμν=12((μν+μν)ημν(00mmθ~θ~γLδ(0)γLδ(xL)))ϕ(𝐱,θ~)ϕ(𝐱,θ~)\hat{T}^{\mu\nu}_{{\bf x},{\bf x^{\prime}},\theta,\theta^{\prime}}=\frac{1}{2}\bigg{(}(\partial^{\mu}\partial^{{}^{\prime}\nu}+\partial^{\mu}\partial^{{}^{\prime}\nu})-\eta^{\mu\nu}\Big{(}\partial_{0}\partial^{{}^{\prime}}_{0}-\partial_{m}\partial^{{}^{\prime}}_{m}-\partial_{\tilde{\theta}}\partial^{{}^{\prime}}_{\tilde{\theta}}-\frac{\gamma}{L}\delta(0)-\frac{\gamma^{\prime}}{L}\delta(x-L)\Big{)}\bigg{)}\phi({\bf x},\tilde{\theta})\phi({\bf x^{\prime}},\tilde{\theta}^{\prime}) (4.19)

where θ~=θλ\tilde{\theta}=\frac{\theta}{\lambda} and has dimension of LL. We express this Energy-Momentum tensor as an operator acting on the time-ordered product of DFRA scalar fields, i.e.,

T^𝐱,𝐱,θ,θμν=O^λ,θμν(,)T(ϕ(𝐱,θ~)ϕ(𝐱,θ~)),\hat{T}^{\mu\nu}_{{\bf x},{\bf x^{\prime}},\theta,\theta^{\prime}}=\hat{O}_{\lambda,\theta}^{\mu\nu}(\partial,\partial^{{}^{\prime}})T(\phi({\bf x},\tilde{\theta})\phi({\bf x^{\prime}},\tilde{\theta}^{\prime})), (4.20)

where

O^λ,θμν(,)=12((μν+μν)ημν(00mmθ~θ~γLδ(x)γLδ(xL))).\hat{O}_{\lambda,\theta}^{\mu\nu}(\partial,\partial^{{}^{\prime}})=\frac{1}{2}\bigg{(}(\partial^{\mu}\partial^{{}^{\prime}\nu}+\partial^{\mu}\partial^{{}^{\prime}\nu})-\eta^{\mu\nu}\Big{(}\partial_{0}\partial^{{}^{\prime}}_{0}-\partial_{m}\partial^{{}^{\prime}}_{m}-\partial_{\tilde{\theta}}\partial_{\tilde{\theta}}^{{}^{\prime}}-\frac{\gamma}{L}\delta(x)-\frac{\gamma^{\prime}}{L}\delta(x-L)\Big{)}\bigg{)}. (4.21)

Note here that the λ\lambda and θ\theta in the subscript of the O^\hat{O} is to imply that operator has terms due to non-commutativity of space-time. As λ0\lambda\rightarrow 0, θ=0\theta=0 we get back the commutative result.

Next we take vacuum expectation value of the Energy-Momentum tensor given in Eqn.(4.20) and find,

<T^𝐱,𝐱,θ,θμν>\displaystyle<\hat{T}^{\mu\nu}_{{\bf x},{\bf x^{\prime}},\theta,\theta^{\prime}}> =\displaystyle= O^λ,θμν(,)<T(ϕ(𝐱,θ~)ϕ(𝐱,θ~))>\displaystyle\hat{O}_{\lambda,\theta}^{\mu\nu}(\partial,\partial^{{}^{\prime}})<T(\phi({\bf x},\tilde{\theta})\phi({\bf x^{\prime}},\tilde{\theta}^{\prime}))>
=\displaystyle= i12((μν+μν)ημν(00mmθ~θ~))Gλ,θ(𝐱,𝐱;θ~,θ~)\displaystyle-i\frac{1}{2}\bigg{(}(\partial^{\mu}\partial^{{}^{\prime}\nu}+\partial^{\mu}\partial^{{}^{\prime}\nu})-\eta^{\mu\nu}\Big{(}\partial_{0}\partial^{{}^{\prime}}_{0}-\partial_{m}\partial^{{}^{\prime}}_{m}-\partial_{\tilde{\theta}}\partial_{\tilde{\theta}}^{{}^{\prime}}\Big{)}\bigg{)}G_{\lambda,\theta}({\bf x,x^{\prime}};\tilde{\theta},\tilde{\theta}^{\prime})

In the above, we have used the identity i<T(ϕϕ)>=Gi<T(\phi\phi^{\prime})>=G, where G is Green’s function.

In the limit, λ0\lambda\rightarrow 0 and θ=0\theta=0 we find G(𝐱,𝐱;θ~,θ~)G({\bf x,x^{\prime}};\tilde{\theta},\tilde{\theta}^{\prime}) reduces to G(𝐱,𝐱)G({\bf x,x^{\prime}})-the Green’s function for the commutative case. Now from above equation, we evaluate vacuum expectation value of xxxx component of Energy-Momentum tensor as

<T^xx>\displaystyle<\hat{T}^{xx}> =\displaystyle= i(1200+12xx12yy12zz12θ~θ~)G(𝐱,𝐱;θ~,θ~).\displaystyle-i\bigg{(}\frac{1}{2}\partial_{0}\partial^{\prime}_{0}+\frac{1}{2}\partial_{x}\partial^{\prime}_{x}-\frac{1}{2}\partial_{y}\partial^{\prime}_{y}-\frac{1}{2}\partial_{z}\partial^{\prime}_{z}-\frac{1}{2}\partial_{\tilde{\theta}}\partial_{\tilde{\theta}}^{{}^{\prime}}\bigg{)}G({\bf x,x^{\prime}};\tilde{\theta},\tilde{\theta}^{\prime}). (4.23)

Next we compactify θ\theta direction (see the discussion after Eqn.(4.8)) by defining η=θ~R\eta=\frac{\tilde{\theta}}{R} and η=θ~R\eta^{\prime}=\frac{\tilde{\theta}^{\prime}}{R}, where η\eta and η\eta^{\prime} varies from 02π0\rightarrow 2\pi and RR is the size of the compact dimension η\eta. Thus above equation becomes

<T^xx>\displaystyle<\hat{T}^{xx}> =\displaystyle= i(1200+12xx12yy12zz12R2ηη)G(𝐱,𝐱;η,η),\displaystyle-i\bigg{(}\frac{1}{2}\partial_{0}\partial^{\prime}_{0}+\frac{1}{2}\partial_{x}\partial^{\prime}_{x}-\frac{1}{2}\partial_{y}\partial^{\prime}_{y}-\frac{1}{2}\partial_{z}\partial^{\prime}_{z}-\frac{1}{2R^{2}}\partial_{\eta}\partial_{\eta}^{{}^{\prime}}\bigg{)}G({\bf x,x^{\prime}};\eta,\eta^{\prime}), (4.24)

Next we use GG given in Eqn.(4.7) and the reduced Energy-Momentum tensor defined through

<T^μν>=dω2πdpy2πdpz2πt^μν<\hat{T}^{\mu\nu}>=\int\frac{d\omega}{2\pi}\frac{dp_{y}}{2\pi}\frac{dp_{z}}{2\pi}\hat{t}^{\mu\nu} (4.25)

to obtain xxxx component of reduced stress tensor

t^xx=12i(ω2p2+xx1R2ηη)g(x,x;η,η;ω;p)|x=x;η=η\hat{t}_{xx}=\frac{1}{2i}\bigg{(}\omega^{2}-p_{\|}^{2}+\partial_{x}\partial^{\prime}_{x}-\frac{1}{R^{2}}\partial_{\eta}\partial_{\eta}^{{}^{\prime}}\bigg{)}g(x,x^{{}^{\prime}};\eta,\eta^{\prime};\omega;p_{\|})\Big{|}_{x=x^{\prime};\eta=\eta^{\prime}} (4.26)

Thus we get xxxx component of reduced stress tensor for DFRA scalar field in presence of compactified dimension coming due to presence of non-commutativity in space-time.

4.3 Modified Casimir Force and Casimir Energy

In this subsection, we evaluate the Casimir force experienced by the plate at x=Lx=L and corresponding Casimir energy. This is calculated using the discontinuity between the component of stress tensors (or pressures) acting on it from the left side of the plate and the right side of the plate. For this, we first find the reduced energy-momentum tensors txx|x=Lt_{xx}\Big{|}_{x=L^{-}} and txx|x=L+t_{xx}\Big{|}_{x=L+} and by taking the difference between these two we get the Casimir force.

Now using Eqn.(4.14) in Eqn.(4.26), we get pressure acting from the left side on the plate at x=Lx=L as

t^xx|x=L=12i[(qqΔ2γγ(2qL)2)+n\{0}(q~q~Δ~2γγ(2q~L)2)+n\{0}(Cn2+(2πnR)2)(12q~+12q~Δ~{γγ(2q~L)2γ2q~L(1+γ2q~L)e2q~(Lx)γ2q~L(1+γ2q~L)e2q~x})].\hat{t}_{xx}\Big{|}_{x=L^{-}}=\frac{1}{2i}\bigg{[}\Big{(}-q-\frac{q}{\Delta}\frac{2\gamma\gamma^{\prime}}{(2qL)^{2}}\Big{)}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\Big{(}-\tilde{q}-\frac{\tilde{q}}{\tilde{\Delta}}\frac{2\gamma\gamma^{\prime}}{(2\tilde{q}L)^{2}}\Big{)}\\ +\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\bigg{(}C_{n}^{2}+{\Big{(}\frac{2\pi n}{R}\Big{)}}^{2}\bigg{)}\bigg{(}\frac{1}{2\tilde{q}}+\frac{1}{2\tilde{q}\tilde{\Delta}}\Big{\{}\frac{\gamma\gamma^{\prime}}{(2\tilde{q}L)^{2}}-\frac{\gamma}{2\tilde{q}L}\Big{(}1+\frac{\gamma^{\prime}}{2\tilde{q}L}\Big{)}e^{2\tilde{q}(L-x)}-\frac{\gamma^{\prime}}{2\tilde{q}L}\Big{(}1+\frac{\gamma}{2\tilde{q}L}\Big{)}e^{2\tilde{q}x}\Big{\}}\bigg{)}\bigg{]}. (4.27)

Similarly by substituting Eqn.(4.15) in Eqn.(4.26), we find the xxxx component of reduced Energy-Momentum tensor, just to the right of the plate,(i.e., pressure acting from the right side) as

t^xx|x=L+=12i[(q)+n\{0}(q~)+n\{0}(Cn2+(2πnR)2)(12q~+12q~Δ~e2q~(Lx){γ2q~L(1+γ2q~L)e2q~Lγ2q~L(1γ2q~L)})].\hat{t}_{xx}\Big{|}_{x=L^{+}}=\frac{1}{2i}\bigg{[}(-q)+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}(-\tilde{q})\\ +\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\bigg{(}C_{n}^{2}+{\Big{(}\frac{2\pi n}{R}\Big{)}}^{2}\bigg{)}\bigg{(}\frac{1}{2\tilde{q}}+\frac{1}{2\tilde{q}\tilde{\Delta}}e^{2\tilde{q}(L-x)}\Big{\{}-\frac{\gamma^{\prime}}{2\tilde{q}L}\Big{(}1+\frac{\gamma}{2\tilde{q}L}\Big{)}e^{2\tilde{q}L}-\frac{\gamma}{2\tilde{q}L}\Big{(}1-\frac{\gamma^{\prime}}{2\tilde{q}L}\Big{)}\Big{\}}\bigg{)}\bigg{]}. (4.28)

In the above equation, q2=ω2+p2q^{2}=-\omega^{2}+p_{\|}^{2} and q~=Cn2ω2+p2\tilde{q}=\sqrt{C_{n}^{2}-\omega^{2}+p_{\|}^{2}}, where p2=py2+pz2p_{\|}^{2}=p_{y}^{2}+p_{z}^{2}.

Now Casimir force arising due to the vacuum fluctuation, acting on the plate at x=Lx=L is

F\displaystyle F =\displaystyle= <T^xx>|x=L<T^xx>|x=L+\displaystyle<\hat{T}^{xx}>\Big{|}_{x=L^{-}}-<\hat{T}^{xx}>\Big{|}_{x=L^{+}} (4.29)
=\displaystyle= dω2πdpy2πdpz2π(t^xx|x=Lt^xx|x=L+),\displaystyle\int\frac{d\omega}{2\pi}\frac{dp_{y}}{2\pi}\frac{dp_{z}}{2\pi}\Big{(}\hat{t}^{xx}\Big{|}_{x=L^{-}}-\hat{t}^{xx}\Big{|}_{x=L^{+}}\Big{)},

Using Eqn.(4.27) and Eqn.(4.28) in above equation we get Casimir force as

F=12i[dω2πdpy2πdpz2π{2qΔγγ(2qL)2+n\{0}(2q~Δ~γγ(2q~L)2)+n\{0}(Cn2+(2πnR)2)12q~Δ~γγ(2q~L)2}].F=\frac{1}{2i}\bigg{[}\int\frac{d\omega}{2\pi}\frac{dp_{y}}{2\pi}\frac{dp_{z}}{2\pi}\Big{\{}\frac{-2q}{\Delta}\frac{\gamma\gamma^{\prime}}{(2qL)^{2}}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\bigg{(}\frac{-2\tilde{q}}{\tilde{\Delta}}\frac{\gamma\gamma^{\prime}}{(2\tilde{q}L)^{2}}\bigg{)}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\bigg{(}C_{n}^{2}+{\Big{(}\frac{2\pi n}{R}\Big{)}}^{2}\bigg{)}\frac{1}{2\tilde{q}\tilde{\Delta}}\frac{\gamma\gamma^{\prime}}{(2\tilde{q}L)^{2}}\Big{\}}\bigg{]}. (4.30)

where q~=Cn2ω2+py2+pz2\tilde{q}=\sqrt{C_{n}^{2}-\omega^{2}+p_{y}^{2}+p_{z}^{2}}. After rewriting the frequency as ω=iζ\omega=i\zeta and using relations given in Eqn.(4.16) and Eqn.(4.17), after straightforward simplifications we get

F=[dζ4πdpy2πdpz2π{2q(2qLγ+1)(2qLγ+1)e2qL1+n\{0}2q~(2q~Lγ+1)(2q~Lγ+1)e2q~L1n\{0}(Cn2+(2πnR)2)12q~((2q~Lγ+1)(2q~Lγ+1)e2q~L1)}].F=-\bigg{[}\int\frac{d\zeta}{4\pi}\frac{dp_{y}}{2\pi}\frac{dp_{z}}{2\pi}\Big{\{}\frac{2q}{\Big{(}\frac{2qL}{\gamma}+1\Big{)}\Big{(}\frac{2qL}{\gamma^{\prime}}+1\Big{)}e^{2qL}-1}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\frac{2\tilde{q}}{\Big{(}\frac{2\tilde{q}L}{\gamma}+1\Big{)}\Big{(}\frac{2\tilde{q}L}{\gamma^{\prime}}+1\Big{)}e^{2\tilde{q}L}-1}\\ -\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\bigg{(}C_{n}^{2}+{\Big{(}\frac{2\pi n}{R}\Big{)}}^{2}\bigg{)}\frac{1}{2\tilde{q}\Big{(}\Big{(}\frac{2\tilde{q}L}{\gamma}+1\Big{)}\Big{(}\frac{2\tilde{q}L}{\gamma^{\prime}}+1\Big{)}e^{2\tilde{q}L}-1\Big{)}}\Big{\}}\bigg{]}. (4.31)

Here we take strong interaction limit, i.e., γ,γ\gamma,\gamma^{\prime}\rightarrow\infty and we get

Fγ,γ=[dζ4πdpy2πdpz2π{2qe2qL1+n\{0}2q~e2q~L1n\{0}(Cn2+(2πnR)2)12q~(e2q~L1)}].F_{\gamma,\gamma^{\prime}\rightarrow\infty}=-\bigg{[}\int\frac{d\zeta}{4\pi}\frac{dp_{y}}{2\pi}\frac{dp_{z}}{2\pi}\Big{\{}\frac{2q}{e^{2qL}-1}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\frac{2\tilde{q}}{e^{2\tilde{q}L}-1}-\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\bigg{(}C_{n}^{2}+{\Big{(}\frac{2\pi n}{R}\Big{)}}^{2}\bigg{)}\frac{1}{2\tilde{q}\Big{(}e^{2\tilde{q}L}-1\Big{)}}\Big{\}}\bigg{]}. (4.32)

Using relations P=ζ2+py2+pz2P=\sqrt{\zeta^{2}+p_{y}^{2}+p_{z}^{2}} and dζdpydpz=4πP2dPd\zeta dp_{y}dp_{z}=4\pi P^{2}dP (after integration of angular variable), the above equation is re-written as

Fγ,γ=[0P2dP4π2{2qe2qL1+n\{0}2q~e2q~L1n\{0}(Cn2+(2πnR)2)12q~(e2q~L1)}]F_{\gamma,\gamma^{\prime}\rightarrow\infty}=-\bigg{[}\int_{0}^{\infty}\frac{P^{2}dP}{4\pi^{2}}\Big{\{}\frac{2q}{e^{2qL}-1}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\frac{2\tilde{q}}{e^{2\tilde{q}L}-1}-\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\bigg{(}C_{n}^{2}+{\Big{(}\frac{2\pi n}{R}\Big{)}}^{2}\bigg{)}\frac{1}{2\tilde{q}\Big{(}e^{2\tilde{q}L}-1\Big{)}}\Big{\}}\bigg{]} (4.33)

where now q=Pq=P and q~=Cn2+P2\tilde{q}=\sqrt{C_{n}^{2}+P^{2}}. We rewrite the above equation as

Fγ,γ=[0(2PL)2d(2PL)32π2L4{2PL(e2PL1)+n\{0}4Cn2L2+4P2L2(e4Cn2L2+4P2L21)n\{0}(Cn2+(2πnR)2)L24Cn2L2+4P2L2(e4Cn2L2+4P2L21)}].F_{\gamma,\gamma^{\prime}\rightarrow\infty}=-\bigg{[}\int_{0}^{\infty}\frac{(2PL)^{2}d(2PL)}{32\pi^{2}L^{4}}\Big{\{}\frac{2PL}{(e^{2PL}-1)}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\frac{\sqrt{4C_{n}^{2}L^{2}+4P^{2}L^{2}}}{(e^{\sqrt{4C_{n}^{2}L^{2}+4P^{2}L^{2}}}-1)}\\ -\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\bigg{(}C_{n}^{2}+{\Big{(}\frac{2\pi n}{R}\Big{)}}^{2}\bigg{)}\frac{L^{2}}{\sqrt{4C_{n}^{2}L^{2}+4P^{2}L^{2}}\Big{(}e^{\sqrt{4C_{n}^{2}L^{2}+4P^{2}L^{2}}}-1\Big{)}}\Big{\}}\bigg{]}. (4.34)

Note that the Cn2C_{n}^{2}-dependent terms are all due to the non-commutativity of the space-time (see Eqn.(4.12)). Next we rewrite the above equation as

Fγ,γ=132π2L40{Y3dY(eY1)+n\{0}Y24Cn2L2+Y2dY(e4Cn2L2+Y21)n\{0}(Cn2+(2πnR)2)Y2L2dY4Cn2L2+Y2(e4Cn2L2+Y21)}F_{\gamma,\gamma^{\prime}\rightarrow\infty}=-\frac{1}{32\pi^{2}L^{4}}\int_{0}^{\infty}\bigg{\{}\frac{Y^{3}dY}{(e^{Y}-1)}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\frac{Y^{2}\sqrt{4C_{n}^{2}L^{2}+Y^{2}}dY}{(e^{\sqrt{4C_{n}^{2}L^{2}+Y^{2}}}-1)}\\ -\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\bigg{(}C_{n}^{2}+{\Big{(}\frac{2\pi n}{R}\Big{)}}^{2}\bigg{)}\frac{Y^{2}L^{2}dY}{\sqrt{4C_{n}^{2}L^{2}+Y^{2}}\Big{(}e^{\sqrt{4C_{n}^{2}L^{2}+Y^{2}}}-1\Big{)}}\bigg{\}} (4.35)

where we use Y=2PLY=2PL. After carrying out the integration, we obtain the modified Casimir force as

Fγ,γ=[π2480L4+n\{0}m=1{(332π2L2)(Cnm)2K2(2LmCn)+116π2L(3Cn2(2πnR)2)CnmK1(2LmCn)}].F_{\gamma,\gamma^{\prime}\rightarrow\infty}=-\bigg{[}\frac{\pi^{2}}{480L^{4}}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\sum_{m=1}\bigg{\{}\Big{(}\frac{3}{32\pi^{2}L^{2}}\Big{)}\Big{(}\frac{C_{n}}{m}\Big{)}^{2}K_{2}(2LmC_{n})+\frac{1}{16\pi^{2}L}\bigg{(}3C_{n}^{2}-\Big{(}\frac{2\pi n}{R}\Big{)}^{2}\bigg{)}\frac{C_{n}}{m}K_{1}(2LmC_{n})\bigg{\}}\bigg{]}. (4.36)

Note here that Cn2=(2πn)2R2+η2λ2(i2πn)C_{n}^{2}=\frac{(2\pi n)^{2}}{R^{2}}+\frac{\eta}{2\lambda^{2}}(i2\pi n) and Cn2C_{n}^{2} dependent term in the above equation is purely due to the non-commutativity of space-time. From expression of Cn2C_{n}^{2} we can see that the Casimir force is complex, which is not expected. The complex part in the Cn2C_{n}^{2} is coming due to presence of the weight function W(θ)W(\theta) in the theory (see Eqn.(3.5)). CnC_{n} can be made real by taking weigh function to be constant (say W(θ)=1)W(\theta)=1) as any θ\theta dependent form of weight function will make the Casimir force and energy complex. Thus by taking (W(θ)=1W(\theta)=1) in our study we find the Casimir force to be

Fγ,γ=[π2480L4+n\{0}m=1{(38L2R2)(nm)2K2(4LmnπR)+n3πmLR3K1(4LmnπR)}]F_{\gamma,\gamma^{\prime}\rightarrow\infty}=-\bigg{[}\frac{\pi^{2}}{480L^{4}}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\sum_{m=1}\bigg{\{}\Big{(}\frac{3}{8L^{2}R^{2}}\Big{)}\Big{(}\frac{n}{m}\Big{)}^{2}K_{2}\Big{(}\frac{4Lmn\pi}{R}\Big{)}+\frac{n^{3}\pi}{mLR^{3}}K_{1}\Big{(}\frac{4Lmn\pi}{R}\Big{)}\bigg{\}}\bigg{]} (4.37)

which is similar to the result obtained in [40] for the Casimir effect in presence of extra one compactified dimension. Here K1K_{1} and K2K_{2} are modified Bessel function of second kind, which have dependency on seperation length, LL and size of extra compactified dimension, RR.

Note here that RR is the size of the compactified dimension due to presence of noncommutativity in space-time and in the limit R0R\rightarrow 0, above expression of Casimir force reduce to the commutative Casimir force expression [19, 20].

We also derive the modified Casimir energy as

Eγ,γ=Fγ,γ𝑑L=π21440L3n\{0}m=1{(332LR2)(nm)2G2,11,3(4LmnπR,12|321,1,12)+n3π4mR3G2,11,3(4LmnπR,12|3212,12,0)}E_{\gamma,\gamma^{\prime}\rightarrow\infty}=-\int F_{\gamma,\gamma^{\prime}\rightarrow\infty}dL\\ =-\frac{\pi^{2}}{1440L^{3}}-\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\sum_{m=1}\bigg{\{}\Big{(}\frac{3}{32LR^{2}}\Big{)}\Big{(}\frac{n}{m}\Big{)}^{2}G_{2,1}^{1,3}\left(\frac{4Lmn\pi}{R},\frac{1}{2}\middle|\begin{array}[]{c}\frac{3}{2}\\ -1,1,\frac{1}{2}\end{array}\right)\\ +\frac{n^{3}\pi}{4mR^{3}}G_{2,1}^{1,3}\left(\frac{4Lmn\pi}{R},\frac{1}{2}\middle|\begin{array}[]{c}\frac{3}{2}\\ -\frac{1}{2},\frac{1}{2},0\end{array}\right)\bigg{\}} (4.38)

where in the above, G2,11,3G_{2,1}^{1,3} is the Meijer-G funtion. Note here that Casimir energy has correction terms, LL dependent Meijer-G function whose coefficient vary as 1L\frac{1}{L}. In the commutative limit, i.e., R0R\rightarrow 0, we get back the commutative result for the Casimir energy [19, 20].

Note that apart from terms proportional to L4L^{-4}, non-commutativity also bring terms that vary as L2L^{-2} and L1L^{-1} in the Casimir force (see Eqn.(4.37)). And also Casimir force has correction terms that vary as 1R2\frac{1}{R^{2}} and 1R3\frac{1}{R^{3}}, where RR is the size of the extra compactified dimension θ\theta due to presence of non-commutativity in the space-time.

5 Thermal Correction to Casimir effect

In this section, we study the finite temperature correction to the Casimir effect in DFR space-time and analyze the modifications to the Casimir effect at high temperature and low-temperature limits. To obtain finite temperature correction to the Casimir force we substitute ζ=2πmβ\zeta=\frac{2\pi m}{\beta} [20, 22] in Eqn.(4.32). After this substitution, we replace the integration with summation as

dζ2π2βm=0,mZ\int_{-\infty}^{\infty}\frac{d\zeta}{2\pi}~{}\rightarrow~{}\frac{2}{\beta}\sum_{m=0}^{\infty}\prime,~{}~{}m~{}\in~{}Z (5.1)

where \prime implies the counting of m=0m=0 term with 12\frac{1}{2} factor and here β=1KT\beta=\frac{1}{KT}, where K is the Boltzmann constant. Using this prescription, we re-express the force expression in Eqn.(4.32) as

Fγ,γ=2βm=0dpy4πdpz2π[2qe2qL1+n\{0}2q~e2q~L1n\{0}(Cn2+(2πnR)2)12q~(e2q~L1)],F_{\gamma,\gamma^{\prime}\rightarrow\infty}=-\frac{2}{\beta}\sum_{m=0}^{\infty}{\prime}\int\frac{dp_{y}}{4\pi}\frac{dp_{z}}{2\pi}\bigg{[}\frac{2q}{e^{2qL}-1}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\frac{2\tilde{q}}{e^{2\tilde{q}L}-1}-\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\bigg{(}C_{n}^{2}+{\Big{(}\frac{2\pi n}{R}\Big{)}}^{2}\bigg{)}\frac{1}{2\tilde{q}\Big{(}e^{2\tilde{q}L}-1\Big{)}}\bigg{]}, (5.2)

where, q=ζ2+py2+pz2q=\sqrt{\zeta^{2}+p_{y}^{2}+p_{z}^{2}}, q~=Cn2+ζ2+py2+pz2\tilde{q}=\sqrt{C_{n}^{2}+\zeta^{2}+p_{y}^{2}+p_{z}^{2}}, and Cn2=(2πn)2R2C_{n}^{2}=\frac{(2\pi n)^{2}}{R^{2}} (note that we have set W(θ)=1W(\theta)=1). Note here that apart from the overall multiplication of 2β\frac{2}{\beta} in above equation, β\beta dependency is appearing through q~\tilde{q} also. Next we use polar form of measure dpydpz=2πpdpdp_{y}dp_{z}=2\pi p_{\|}dp_{\|} where p=py2+pz2p_{\|}=\sqrt{p_{y}^{2}+p_{z}^{2}}, re-define dp2=dq2dp_{\|}^{2}=dq^{2} and

t=4πLβt=\frac{4\pi L}{\beta} (5.3)

and re-express above equation as

Fγ,γ=2βm=0[mtY2dY16πL3(eY1)+n\{0}(mt)2+4L2Cn2Y~2dY~16πL3(eY~1)n\{0}Cn2(mt)2+4L2Cn2dY~8πL(eY~1)].F_{\gamma,\gamma^{\prime}\rightarrow\infty}=-\frac{2}{\beta}\sum_{m=0}^{\infty}{\prime}\bigg{[}\int_{mt}^{\infty}\frac{Y^{2}dY}{16\pi L^{3}(e^{Y}-1)}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\int_{\sqrt{(mt)^{2}+4L^{2}C_{n}^{2}}}^{\infty}\frac{\tilde{Y}^{2}d\tilde{Y}}{16\pi L^{3}(e^{\tilde{Y}}-1)}\\ -\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}C_{n}^{2}\int_{\sqrt{(mt)^{2}+4L^{2}C_{n}^{2}}}^{\infty}\frac{d\tilde{Y}}{8\pi L(e^{\tilde{Y}}-1)}\bigg{]}. (5.4)

Here Y=2qL=(mt)2+4L2p2Y=2qL=\sqrt{(mt)^{2}+4L^{2}p_{\|}^{2}} and Y~=2q~L=(mt)2+4L2Cn2+4L2p2\tilde{Y}=2\tilde{q}L=\sqrt{(mt)^{2}+4L^{2}C_{n}^{2}+4L^{2}p_{\|}^{2}}. Note that through YY and Y~\tilde{Y}, dependence on t ( i.e., β\beta dependence) is entering the above expression for force (see Eqn.(5.3)). We investigate the high temperature limit and the low temperature limit of the finite temperature correction of the Casimir effect as follows.


For high temperature limit: since from Eqn.(5.3) we get t=4πLKTt=4\pi LKT and thus for large value of T, we find t=4πLKT>>1t=4\pi LKT>>1 i.e., 4πL>>β4\pi L>>\beta [20, 22]. After evaluation of integration in above equation using zeta function, we get Casimir force in the high temperature limit as

F4πL>>β=[{18πβL3ζ(3)+14πβL3(1+t+t22)et}+116πL3βn\{0}(4L2Cn2+4LCn+2)e2LCn+18πL3βn\{0}((t2+4L2Cn2)+2t2+4L2Cn2+2)et2+4L2Cn214πLβn\{0}Cn2(e2LCn2+et2+4L2Cn2)].F_{4\pi L>>\beta}=-\bigg{[}\Big{\{}\frac{1}{8\pi\beta L^{3}}\zeta(3)+\frac{1}{4\pi\beta L^{3}}(1+t+\frac{t^{2}}{2})e^{-t}\Big{\}}+\frac{1}{16\pi L^{3}\beta}\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\Big{(}4L^{2}C_{n}^{2}+4LC_{n}+2\Big{)}e^{-2LC_{n}}\\ +\frac{1}{8\pi L^{3}\beta}\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\Big{(}(t^{2}+4L^{2}C_{n}^{2})+2\sqrt{t^{2}+4L^{2}C_{n}^{2}}+2\Big{)}e^{-\sqrt{t^{2}+4L^{2}C_{n}^{2}}}\\ -\frac{1}{4\pi L\beta}\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}C_{n}^{2}\bigg{(}\frac{e^{-2LC_{n}}}{2}+e^{-\sqrt{t^{2}+4L^{2}C_{n}^{2}}}\bigg{)}\bigg{]}. (5.5)

Note here that apart from the explicit β\beta dependency in all the terms, through t (see Eqn.(5.3)) also , β\beta enter the above expression for force at the high-temperature limit . The above expression for the Casimir force for the high temperature limit with Cn2=(2πn)2R2C_{n}^{2}=\frac{(2\pi n)^{2}}{R^{2}} reduce to commutative result [20, 22, 19] in the limit R0R\rightarrow 0, where R is the size of the extra compactified dimension.

For low temperature limit: here we study the low temperature limit of Casimir effect using Poisson summation formula [20, 22, 19]; any function b(X)b(X) and it’s Fourier transform c(α)c(\alpha) i.e.,

c(α)=12πb(X)eiαX𝑑X,c(\alpha)=\frac{1}{2\pi}\int_{-\infty}^{\infty}b(X)e^{-i\alpha X}dX, (5.6)

satisfy the identity

b(m)=2πc(2πm).\sum_{-\infty}^{\infty}b(m)=2\pi\sum_{-\infty}^{\infty}c(2\pi m). (5.7)

Using this in Eqn.(5.4) with

b(m)=[mtY2dY16πL3(eY1)+n\{0}(mt)2+4L2Cn2Y~2dY~16πL3(eY~1)n\{0}Cn2(mt)2+4L2Cn2dY~8πL(eY~1)]|m=X,b(m)=\bigg{[}\int_{mt}^{\infty}\frac{Y^{2}dY}{16\pi L^{3}(e^{Y}-1)}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\int_{\sqrt{(mt)^{2}+4L^{2}C_{n}^{2}}}^{\infty}\frac{\tilde{Y}^{2}d\tilde{Y}}{16\pi L^{3}(e^{\tilde{Y}}-1)}\\ -\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}C_{n}^{2}\int_{\sqrt{(mt)^{2}+4L^{2}C_{n}^{2}}}^{\infty}\frac{d\tilde{Y}}{8\pi L(e^{\tilde{Y}}-1)}\bigg{]}\bigg{|}_{m=X}, (5.8)

where Y=2qL=(mt)2+4L2p2Y=2qL=\sqrt{(mt)^{2}+4L^{2}p_{\|}^{2}} and Y~=2q~L=(mt)2+4L2Cn2+4L2p2\tilde{Y}=2\tilde{q}L=\sqrt{(mt)^{2}+4L^{2}C_{n}^{2}+4L^{2}p_{\|}^{2}}, we find

c(α)=12πeiαXdX[XtY2dY16πL3(eY1)+n\{0}(Xt)2+4L2Cn2Y~2dY~16πL3(eY~1)n\{0}Cn2(Xt)2+4L2Cn2dY~8πL(eY~1)]=1π0cos(αX)dX[XtY2dY16πL3(eY1)+n\{0}(Xt)2+4L2Cn2Y~2dY~16πL3(eY~1)n\{0}Cn2(Xt)2+4L2Cn2dY~8πL(eY~1)];\begin{split}c(\alpha)=\frac{1}{2\pi}\int_{-\infty}^{\infty}e^{-i\alpha X}dX\bigg{[}\int_{Xt}^{\infty}\frac{Y^{2}dY}{16\pi L^{3}(e^{Y}-1)}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\int_{\sqrt{(Xt)^{2}+4L^{2}C_{n}^{2}}}^{\infty}\frac{\tilde{Y}^{2}d\tilde{Y}}{16\pi L^{3}(e^{\tilde{Y}}-1)}\\ -\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}C_{n}^{2}\int_{\sqrt{(Xt)^{2}+4L^{2}C_{n}^{2}}}^{\infty}\frac{d\tilde{Y}}{8\pi L(e^{\tilde{Y}}-1)}\bigg{]}\\ =\frac{1}{\pi}\int_{0}^{\infty}cos(\alpha X)dX\bigg{[}\int_{Xt}^{\infty}\frac{Y^{2}dY}{16\pi L^{3}(e^{Y}-1)}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\int_{\sqrt{(Xt)^{2}+4L^{2}C_{n}^{2}}}^{\infty}\frac{\tilde{Y}^{2}d\tilde{Y}}{16\pi L^{3}(e^{\tilde{Y}}-1)}\\ -\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}C_{n}^{2}\int_{\sqrt{(Xt)^{2}+4L^{2}C_{n}^{2}}}^{\infty}\frac{d\tilde{Y}}{8\pi L(e^{\tilde{Y}}-1)}\bigg{]};\end{split} (5.9)

After interchanging the order of integration in above we get

c(α)=1πα[0sin(ZY)Y2dY16πL3(eY1)+n\{0}2LCnsin(ZY~24L2Cn2)Y~2dY~16πL3(eY~1)n\{0}Cn22LCnsin(ZY~24L2Cn2)dY~8πL(eY~1)],c(\alpha)=\frac{1}{\pi\alpha}\bigg{[}\int_{0}^{\infty}sin(ZY)\frac{Y^{2}dY}{16\pi L^{3}(e^{Y}-1)}+\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\int_{2LC_{n}}^{\infty}sin\Big{(}Z\sqrt{\tilde{Y}^{2}-4L^{2}C_{n}^{2}}\Big{)}\frac{\tilde{Y}^{2}d\tilde{Y}}{16\pi L^{3}(e^{\tilde{Y}}-1)}\\ -\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}C_{n}^{2}\int_{2LC_{n}}^{\infty}sin\Big{(}Z\sqrt{\tilde{Y}^{2}-4L^{2}C_{n}^{2}}\Big{)}\frac{d\tilde{Y}}{8\pi L(e^{\tilde{Y}}-1)}\bigg{]}, (5.10)

where Z=αtZ=\frac{\alpha}{t}, t=4πLβt=\frac{4\pi L}{\beta}. We use identity 0sin(ZY)(eY1)𝑑Y=π2cothπZ12Z\int_{0}^{\infty}\frac{sin(ZY)}{(e^{Y}-1)}dY=\frac{\pi}{2}coth\pi Z-\frac{1}{2Z} and an identity given in [49] in above equation and after straightforward simplifications we find

c(α)|α=2πm=132mπ3L3[4π3{e4π2mt(1+e4π2mt)(1e4π2mt)3}t3(2πm)3]+132mπ3L3n\{0}{d2dZ2(n~=12LCnZn~2+Z2K1(2LCnn~2+Z2))+116mπ3Ln\{0}Cn2(n~=12LCnZn~2+Z2K1(2LCnn~2+Z2))]|Z=αt=2πmtc(\alpha)\Big{|}_{\alpha=2\pi m}=-\frac{1}{32m\pi^{3}L^{3}}\bigg{[}4\pi^{3}\Big{\{}\frac{e^{-\frac{4\pi^{2}m}{t}}(1+e^{-\frac{4\pi^{2}m}{t}})}{(1-e^{-\frac{4\pi^{2}m}{t}})^{3}}\Big{\}}-\frac{t^{3}}{(2\pi m)^{3}}\bigg{]}\\ +\frac{1}{32m\pi^{3}L^{3}}\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\bigg{\{}-\frac{d^{2}}{dZ^{2}}\bigg{(}\sum_{\tilde{n}=1}^{\infty}\frac{2LC_{n}Z}{\sqrt{\tilde{n}^{2}+Z^{2}}}K_{1}\Big{(}2LC_{n}\sqrt{\tilde{n}^{2}+Z^{2}}\Big{)}\bigg{)}\\ +\frac{1}{16m\pi^{3}L}\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}C_{n}^{2}\bigg{(}\sum_{\tilde{n}=1}^{\infty}\frac{2LC_{n}Z}{\sqrt{\tilde{n}^{2}+Z^{2}}}K_{1}\Big{(}2LC_{n}\sqrt{\tilde{n}^{2}+Z^{2}}\Big{)}\bigg{)}\bigg{]}\bigg{|}_{Z=\frac{\alpha}{t}=\frac{2\pi m}{t}} (5.11)

Now using Eqn.(5.4), Eqn.(5.7) and Eqn.(5.11), we get the low temperature limit of the finite thermal correction for Casimir force as

Fβ>>4πL={π2480L4(1+t448π460tπ2e4π2t)}132π2L2n\{0}{(πLCn)(2LCn+92)e2LCn}116π2L2n\{0}Cn2(πLCn)e2LCn116π2L4n\{0}(πLCn){2LCnt52(t2+4π2)54+92t72(t2+4π2)74+16π2L2Cn2t32(t2+4π2)748π2LCnt52(t2+4π2)94+21π2t92(t2+4π2)114}e2LCntt2+4π218π2L2n\{0}Cn2(πLCn)(t32(t2+4π2)34)e2LCntt2+4π2.F_{\beta>>4\pi L}=-\bigg{\{}\frac{\pi^{2}}{480L^{4}}\Big{(}1+\frac{t^{4}}{48\pi^{4}}-\frac{60t}{\pi^{2}}e^{-\frac{4\pi^{2}}{t}}\Big{)}\bigg{\}}\\ -\frac{1}{32\pi^{2}L^{2}}\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}\Big{\{}(\sqrt{\pi LC_{n}})\Big{(}2LC_{n}+\frac{9}{2}\Big{)}e^{-2LC_{n}}\Big{\}}-\frac{1}{16\pi^{2}L^{2}}\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}C_{n}^{2}(\sqrt{\pi LC_{n}})e^{-2LC_{n}}\\ -\frac{1}{16\pi^{2}L^{4}}\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}(\sqrt{\pi LC_{n}})\bigg{\{}\frac{2LC_{n}t^{\frac{5}{2}}}{(t^{2}+4\pi^{2})^{\frac{5}{4}}}+\frac{9}{2}\frac{t^{\frac{7}{2}}}{(t^{2}+4\pi^{2})^{\frac{7}{4}}}+\frac{16\pi^{2}L^{2}C_{n}^{2}t^{\frac{3}{2}}}{(t^{2}+4\pi^{2})^{\frac{7}{4}}}-\frac{8\pi^{2}LC_{n}t^{\frac{5}{2}}}{(t^{2}+4\pi^{2})^{\frac{9}{4}}}\\ +\frac{21\pi^{2}t^{\frac{9}{2}}}{(t^{2}+4\pi^{2})^{\frac{11}{4}}}\bigg{\}}e^{-2LC_{n}t\sqrt{t^{2}+4\pi^{2}}}-\frac{1}{8\pi^{2}L^{2}}\sum_{n\in\mathbb{Z}\backslash\left\{0\right\}}C_{n}^{2}(\sqrt{\pi LC_{n}})\bigg{(}\frac{t^{\frac{3}{2}}}{(t^{2}+4\pi^{2})^{\frac{3}{4}}}\bigg{)}e^{-2LC_{n}t\sqrt{t^{2}+4\pi^{2}}}. (5.12)

Note here that at T=0T=0, above force expression reduce to commutative results [20, 22] with additional temperature independent correction terms due to non-commutativity of space-time.

Note that, the commutative limit is obtained by letting R0R\rightarrow 0, where RR is the size of the extra compactified dimension. We find that the thermal corrections to Casimir force (both at the high and low-temperature limit) reduce to known commutative result when R0R\rightarrow 0[20, 22]. At the high-temperature limit, Casimir force has, apart from the L3L^{-3} dependent corrections, L1L^{-1} dependent corrections also. Similarly, at a low-temperature limit, Casimir force has a dependency on L4L^{-4} terms and also additional L2L^{-2} dependent terms due to non-commutativity. Through expansion of exponential part of the correction terms, we find that in both high and low-temperature limits, Casimir force also has a LL-independent corrections due to non-commutativity.

6 Conclusion

In this paper, we have studied the effect of vacuum fluctuations in 4+14+1-dimensional DFR space-time by analyzing the Casimir effect between parallel plates. We have found modifications to the Casimir force and the Casimir energy between two parallel plates due to the non-commutativity of DFR space-time. Here we found that if the additional θ\theta-directions are treated on the same footing as transverse commutative directions, the Casimir force (and energy) found to be exactly same as that for commutative extra dimensional space-time [20] and does not have dependence on the non-commutative parameters. We then evaluated the Casimir effect by treating the non-commutative direction as compactified dimension. We find that the Casimir force and energy becomes complex if the weight function W(θ)W(\theta) introduced in Eqn.(3.3) is not a constant. Hence we take W(θ)=1W(\theta)=1 as in the studied in [17] in our analysis. In this case, we have seen that Casimir force proportional to 1L4\frac{1}{L^{4}} with modifications that vary as 1L2\frac{1}{L^{2}} and 1L\frac{1}{L} in 4+14+1 dimensional DFR space-time. And Casimir energy vary as 1L3\frac{1}{L^{3}} with correction term due to non-commutativity of space-time vary as 1L\frac{1}{L}. This result is in agreement with similar results obtained in [38, 39, 40] for Casimir effect in commutative space-time with extra compactified-dimension. But our result is in contrast with the results of [20], where Casimir force and Casimir energy scale as 1L5\frac{1}{L^{5}} and 1L4\frac{1}{L^{4}}, respectively in 4+14+1 dimensions, where extra dimension was a commutative one.

In [38, 39, 40], the Casimir effect was studied in the presence of extra dimension and found that Casimir force has terms that vary as 1L4\frac{1}{L^{4}} apart from correction terms that vary as 1L2\frac{1}{L^{2}} and 1L\frac{1}{L} due to presence of one extra dimension in addition to 44 dimensional commutative Minkowski space-time. In [20], Casimir effect is studied for parallel conducting plates in arbitrary spatial dimension D and showed that Casimir force and Casimir energy scale as L(D+1)L^{-(D+1)} and L(D)L^{-(D)} respectively. For spatial dimension D=4D=4, these results are in contrast with the results obtained in [38, 39, 40]. This difference in the L dependence observed in [38, 39, 40] and [20] is because of the fact that in [20], the extra dimension is treated as transverse direction to the plates, whereas in [38, 39, 40], the effect of the extra dimension on parallel plates in usual 3+13+1 dimensional Minkowski space-time is studied. The approach taken in [38, 39, 40] is similar to our present study.

Note that we have obtained the corrections to Casimir force that is attractive in nature (see Eqn.(4.37)). We have showed that in commutative limit, R0R\rightarrow 0 our results reduce to that in [19]. Here RR is the size of the extra compactified dimension - θ\theta. Note that the Casimir force expression in Eqn.(4.37) reduces to π2480L4-\frac{\pi^{2}}{480L^{4}} in the commutative limit. Casimir force expression (see Eqn.(4.37))has correction terms that vary as 1R2\frac{1}{R^{2}} and 1R3\frac{1}{R^{3}}. In this Casimir force expression, except correction term with n=1n=1 and m=1m=1, contribution of correction terms with higher values of n,mn,m, are very small. Thus to exhibit the generic nature of the Casimir force, we neglect those terms and for various possible values of size of the extra compactified dimension, RR, we see how the Casimir force is varying with plate separation (see the plot).

Refer to caption
Figure 1: Variation of Casimir force with plate separation.

In the above plot, NC1 is the plot for R107mR\sim 10^{-7}m and NC2 is for R105mR\sim 10^{-5}m. The solid line gives the Casimir force for commutative case.

Here the solid, continuous line shows the variation of the Casimir force in the commutative 44-dimensional space-time for parallel plates as a function of plate separation. For R<107mR<10^{-7}m, corrections in the Casimir force, are very small. All the plots overlap with solid line (Casimir force for the commutative case) for large plate separation. The Casimir force in DFR space-time, for values of R107mR\sim 10^{-7}m and R105mR\sim 10^{-5}m are larger than the commutative case for a given plate seperation. That is, modifications to the Casimir force strengthen the attractive nature of the Casimir force. But for R<107mR<10^{-7}m, as corrections to the Casimir force are very small that we do not find any deviation from the commutative case in the plot. Thus the present analysis constrain the size of the extra compactified dimension with lower bound, i.e., R107mR\geq 10^{-7}m.

We have also investigated finite temperature corrections to the Casimir effect, in presence of extra compactified dimension, in 4+14+1 dimensional DFR space-time. For the high-temperature limit, we have shown that the Casimir force has terms that scale as 1L3\frac{1}{L^{3}} and 1L\frac{1}{L} respectively. 1L\frac{1}{L} term is purely due to the non-commutativity of space-time. Similarly, we have studied the low-temperature limit of the Casimir force, and it is shown to have terms that vary as 1L4\frac{1}{L^{4}} and 1L2\frac{1}{L^{2}}. Here 1L2\frac{1}{L^{2}} is purely due to the non-commutativity of space-time. In [50, 51], finite temperature correction to Casimir effect due to the presence of one extra commutative dimension was studied and the result shows the same feature we find here in DFR space-time.

In the appendix, we have investigated the Casimir effect for 3+13+1 dimensional massive scalar field theory [20] obtained from 77-dimensional massless DFRA scalar field theory in DFR space-time by applying the Kaluza-Klein reduction method.

Acknowledgements

SKP thanks UGC, India, for support through the JRF scheme(id.191620059604).


Appendix : A Kaluza-Klein reduction and Casimir effect

In this appendix, we study the Casimir effect for massive scalar field theory in 3+1 dimensional Minkowski space-time. This model is obtained by applying the Kaluza-Klein reduction prescription to 7-dimensional DFRA massless scalar field theory.

We start with the action of massless complex scalar field theory in DFR space-time,

S=d4xd3θW(θ)(μϕ(x,θ)μϕ(x,θ)+θiϕ(x,θ)θiϕ(x,θ)).S=\int d^{4}x~{}d^{3}\theta~{}W(\theta)\bigg{(}\partial_{\mu}\phi(x,\theta)\partial^{\mu}\phi^{\dagger}(x,\theta)+\partial_{\theta^{i}}\phi(x,\theta)\partial^{\theta^{i}}\phi^{\dagger}(x,\theta)\bigg{)}. (A.1)

We rewrite this space-time as direct product space, M7=M4K3M_{7}=M_{4}\otimes K_{3}, where K3K_{3} is a compact space. Using Kaluza-Klein reduction, from this 7-dimensional massless scalar field theory, we obtain a massive scalar field theory in 3+13+1 dimensional Minkowski space-time. For this we take Fourier decomposition of DFRA complex scalar field as,

ϕ(x,θ)=N=ϕ(N)(x)exp(iNMθλ).\phi(x,\theta)=\sum_{N=-\infty}^{\infty}\phi^{(N)}(x)\exp{\Big{(}-\frac{iNM\theta}{\lambda}\Big{)}}. (A.2)

Using the above in Eq.(A.1), we get

S=N=d4xd3θW(θ)(μϕ(N)(x)μϕ(N)(x)+(N2M2)ϕ(N)(x)ϕ(N)(x)).S=\sum_{N=-\infty}^{\infty}\int d^{4}x~{}d^{3}\theta~{}W(\theta)\bigg{(}\partial_{\mu}\phi^{(N)}(x)\partial^{\mu}\phi^{\dagger~{}(N)}(x)+(N^{2}M^{2})\phi^{(N)}(x)\phi^{\dagger~{}(N)}(x)\bigg{)}. (A.3)

Now after integrating over θ\theta from 0 to 2π/λ2\pi/\lambda , we find

S=N=d4x(erf(πλ3))3(μϕ(N)(x)μϕ(N)(x)+(N2M2)ϕ(N)(x)ϕ(N)(x)),S=\sum_{N=-\infty}^{\infty}\int d^{4}x~{}\Big{(}erf\Big{(}\frac{\pi}{\lambda^{3}}\Big{)}\Big{)}^{3}\bigg{(}\partial_{\mu}\phi^{(N)}(x)\partial^{\mu}\phi^{\dagger~{}(N)}(x)+(N^{2}M^{2})\phi^{(N)}(x)\phi^{\dagger~{}(N)}(x)\bigg{)}, (A.4)

where erf()erf() is the error function. Next for λ0\lambda\rightarrow 0, erf(πλ3)1erf(\frac{\pi}{\lambda^{3}})\rightarrow 1 the above equation reduces to

S=N=d4x(μϕ(N)(x)μϕ(N)(x)+(N2M2)ϕ(N)(x)ϕ(N)(x)).S=\sum_{N=-\infty}^{\infty}\int d^{4}x~{}\bigg{(}\partial_{\mu}\phi^{(N)}(x)\partial^{\mu}\phi^{\dagger~{}(N)}(x)+(N^{2}M^{2})\phi^{(N)}(x)\phi^{\dagger~{}(N)}(x)\bigg{)}. (A.5)

Thus we get the action of an infinite number of uncoupled modes of the massive complex scalar field in the usual Minkowski space-time. Using this, we can study the Casimir effect, and for single-mode, the corresponding Casimir force and Casimir energy expressions [20] are

Fc=132π2L4[34n=1(4ML)2n2K2(2nML)+2M2L2n=1(4MLn)K1(2nML)]F_{c}=-\frac{1}{32\pi^{2}L^{4}}\bigg{[}\frac{3}{4}\sum_{n=1}^{\infty}\frac{(4ML)^{2}}{n^{2}}K_{2}(2nML)+2M^{2}L^{2}\sum_{n=1}^{\infty}\Big{(}\frac{4ML}{n}\Big{)}K_{1}(2nML)\bigg{]} (A.6)

and

Ec=132π2L3n=1(2ML)2n2K2(2nML)E_{c}=-\frac{1}{32\pi^{2}L^{3}}\sum_{n=1}^{\infty}\frac{(2ML)^{2}}{n^{2}}K_{2}(2nML) (A.7)

respectively where K1K_{1} and K2K_{2} are modified Bessel functions. These results match exactly with the results of the Casimir force for massive scalar field theory given in [20]. Thus we show that one can obtain the usual massive scalar field theory by applying Kaluza-Klein reduction to the massless DFRA scalar field theory, and this will lead to the Casimir effect as in the well studied massive scalar theory in 3+13+1 dimensions.

References

  • [1] A. Connes, Non-Commutative Geometry (Academic Press, London, 1994).
  • [2] M. R. Douglas and N. A. Nekrasov, Rev. Mod. Phys. 73 (2001) 977; R. J. Szabo, Phys. Rep. 378 (2003) 207.
  • [3] J. Kowalski-Glikman, Lect. Notes. Phys. 669 (2005) 131.
  • [4] S. Doplicher, K. Fredenhagen and J. E. Roberts, Phys. Lett. B 331 (1994) 29; S. Doplicher, K. Fredenhagen and J. E. Roberts, Commun. Math. Phys. 172 (1995) 187.
  • [5] J. Ambjorn, J. Jurkiewicz, and R. Loll, Phys. Rev. Lett. 93 (2004) 131301.
  • [6] J. Madore, An Introduction to Noncommutative Differential Geometry and Its Applications (Cambridge Univ. Press, 1995)
  • [7] N. Seiberg and E. Witten, JHEP, 9909 (1999) 032.
  • [8] L. Freidel and E. R. Livine, Phys. Rev. Lett. 96 (2006) 221301.
  • [9] H. S. Snyder, Phys. Rev. 71 (1947) 38.
  • [10] C. N. Yang, Phys. Rev. 72 (1947) 874.
  • [11] M. M. Sheikh-Jabbari, Phys. Lett. B 425 (1998) 48.
  • [12] J. J Jaeckel, V. V. Khoze, and A. Ringwald, J. High Energy Phys. 02 (2006) 028.
  • [13] C. E. Carlson, C. D. Carone and N. Zobin, Phys. Rev. D 66 (2002) 075001.
  • [14] R. Amorim, Phys. Rev. D 78 (2008) 105003.
  • [15] R. Amorim, Phys. Rev. Lett. 101 (2008) 081602.
  • [16] R. Amorim, J.Math. Phys. (N.Y) 50 (2009) 022303; R. Amorim, J.Math. Phys. (N.Y) 50 (2009) 052103.
  • [17] R. Amorim, Phys. Rev. D 80 (2009) 105010.
  • [18] H. B. G. Casimir, Koninkl. Ned. Akad. Wetenschap. Proc 51 (1948) 793.
  • [19] K. A. Milton, J. Phys. A: Math. and Gen. 37 (2004) R209.
  • [20] K. A. Milton, The Casimir Effect Physical manifestations of Zero-Point Energy, World Scientific, Singapore, 2001).
  • [21] G. Plunien, B. Muller and W. Greiner, Phys. Rep. 134 (1986) 87.
  • [22] I. Brevik, S.A. Ellingsen and K. A. Milton, New J. Phys. 8 (2006) 236.
  • [23] G. Bressi, G. Carugno, R. Onofrio and G. Ruoso, Phys. Rev. Lett. 88 (2002) 041804.
  • [24] R. Sedmik and P. Brax, J. Phys: Conf. Ser. (2018) 1138012014.
  • [25] R. I. P. Sedmik, Int. J. Mod. Phys. A 35 (2020) 2040008; R. I. P. Sedmik and M. Pitschmann, Universe 7 (2021) 234.
  • [26] U. Mohideen and A. Roy, Phys. Rev. Lett. 81 (1998) 4549; V. M. Mostepanenko, J. Phys: Conf. Ser. 161 (2009) 012003 and reference therein.
  • [27] G. L. Klimchitskaya and V. M. Mostepanenko, Mod. Phys. Lett. A 35 (2020) 2040007. G. Bimonte, B. Spreng, P. A. Maia Neto, G. L. Ingold, G. L. Klimchitskaya, V. M. Mostepanenko and R. S. Decca, Universe 7 (2021) 93.
  • [28] M. Wang, L. Tang, C. Y. Ng, R. Messina, B. Guizal, J. A. Crosse, M. Antezza, C. T. Chan and H. B. Chan, Nature Commun. 12 (2021) 600.
  • [29] R. Onofrio, N. J. Phys. 8 (2008) 237; G. L. Klimchitskaya, U. Mohideen and V. M. Mostepanenko, Phys. Rev. D 87 (2013) 125031.
  • [30] R. Casadio, A. Gruppuso, B. Harms and O. Micu, Phys. Rev. D 76 (2007) 025016.
  • [31] C. D. Fosco and G. A. Moreno, Phys. Lett. B 659 (2008) 901.
  • [32] M. Chaichian, A. Demichev, P. Prešnajder, M. M. Sheikh-Jabbari and A. Tureanu, Nucl. Phys. B 611 (2001) 383.
  • [33] M. V. Cougo-Pinto, C. Farina and J. F. M. Mendes, Nucl. Phys. B 127 (2004) 138.
  • [34] E. Harikumar, S. K. Panja and V. Rajagopal, Nucl. Phys. B 950 (2020) 114842.
  • [35] J. Lukierski, A. Nowicki, H. Ruegg and V. Tolstoy, Phys. Lett. B 264 (1991) 331; J. Lukierski, A. Nowicki and H. Ruegg, Phys. Lett. B 293 (1991) 344.
  • [36] J. Gomis, T. Mehen, M. B. Wise, JHEP 0008 (2000) 029.
  • [37] S. Nam, JHEP 0010 (2000) 044.
  • [38] K. Kirsten and S. A. Fulling, Phys. Rev. D 79 (2009) 065019.
  • [39] H. B. Cheng, Phys. Lett. B 668 (2008) 72; S. A. Fulling and K. Kirsten, Phys. Lett. B 671 (2009) 179.
  • [40] F. Pascoal, L. F. A. Oliveira, F. S. S. Rosa and C. Farina, Braz. J. Phys. 38 (2008) 581.
  • [41] L. P. Teo, Phys. Lett. B 682 (2009) 259.
  • [42] M. Bordag, D. Hennig, and D. Robaschik, J. Phys. A 25 (1992) 4483; M. Bordag, K. Kirsten, and D. Vassilevich, Phys. Rev. D 59 (1999) 085011.
  • [43] J. Gomis and T. Mehen, Nucl. Phys. B 591 (2000) 265.
  • [44] M. Chaichian, P. Kulish, K. Nishijima, and A. Tureanu, Phys. Lett. B 604 (2004) 98; P. Aschieri, C. Blohmann, M. Dimitrijevic, F. Meyer, P. Schupp and J. Wess, Class. Quan. Grav., 22 (2005) 3511.
  • [45] E. M. C. Abreu, A. C. R. Mendes, W. Oliveira and A. Zagirolamim, SIGMA 6 083 (2010); E. M. C. Abreu and M. J. Neves, Int. J. Mod. Phys. A 32 (2017) 1750099.
  • [46] H. Kase, K. Morita, Y. Okumura and E. Umezawa, Prog. Theor. Phys. 109 (2003) 663.
  • [47] K. Imai, K. Morita and Y. Okumura, Prog. Theor. Phys. 110 (2003) 989.
  • [48] S. Saxell, Phys. Lett. B 666 (2008) 486.
  • [49] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, series and products, 7th ed. 2007 (Elsevier Academic Press, 2007), pp.540.
  • [50] M. Rypestøl and I. Brevik Phys. Scr. 82 (2010) 035101.
  • [51] L.P. Teo, JHEP 0906 (2009) 076.