# On the Thermodynamic Properties of the Spinless Duffin-Kemmer-Petiau Oscillator in Noncommutative Plane

Authors: Zhi Wang; Zheng-Wen Long; Chao-Yun Long; Wei Zhang
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

The Duffin-Kemmer-Petiau oscillator for spin 0 particle in noncommutative plane is analyzed and the energy eigenvalue of the system is obtained by employing the functional analysis method. Furthermore, the thermodynamic properties of the noncommutative DKP oscillator are investigated via numerical method and the influence of noncommutative space on thermodynamic functions is also discussed. We show that the energy spectrum and corresponding thermodynamic functions of the considered physical systems depend explicitly on the noncommutative parameterθ which characterizes the noncommutativity of the space.


---

## Body

## 1. Introduction

The Duffin-Kemmer-Petiau (DKP) equation introduced by Kemmer et al. [1–3] in the 1930s is a relativistic wave equation which describes spin 0 and spin 1 bosons in the description of the standard model. The DKP equation for spin 0 hadrons with nuclei is closely linked to the conventional second-order Klein-Gordon (KG) equation, which has more complex algebraic structure than the latter [4], and the DKP equation enjoys a richness of couplings not capable of being expressed in the Klein-Gordon and Proca equations [5, 6]. In addition, several situations researched in past years involving hadronic process and breaking of symmetries show that in some situations the DKP equation and KG equation can give some different results. For these reasons, there has been a growing interest in studying the DKP equation [7–13].On the other hand, the noncommutative (NC) geometry has been put forward in both quantum mechanics and field theory in past years, which is growing rapidly. The coordinate noncommutativity was first introduced by Snyder [14] as a regularized mechanism to improve the problem of infinite self-energy in quantum field theory. But with the development of renormalization theory, this concept became not popular until Connes [15] investigated Yang-Mills theories in noncommutative (NC) space. In recent years, there has been a renewed interest in NC geometry because of the discovery in string theory and matrix model of M-theory [16]. When studying NC space, one considers the coordinate to be noncommutative, where the coordinate operators satisfy the commutation relation [x^i,x^j]=iθij, with θij being a constant antisymmetric matrix, which represent the noncommutative property of NC space, and play an analogous role with ħ in the usual quantum mechanics. Recently, various aspects of both NC classical [17] and quantum [18] mechanics have been investigated by large number of papers, devoted to exploring the role of NC parameter in the physical observables. For example, classical Newton mechanics in the noncommutative space was studied in [19, 20], a particle confined by a quadratic potential in the generalized noncommutative plane was investigated from both the classical and the quantum aspects in [21], the noncommutative harmonic oscillators were discussed in detail in [22–24], and [25] studied the Aharonov-Bohm effect in a class of noncommutative theories.The DKP equation in the presence of the Dirac oscillator interaction called DKP oscillator for spin zero and spin one particles has also been investigated in the NC formalism in publications [26–28], with the purpose of obtaining the energy spectrum and studying the noncommutative effect on relativistic quantum mechanics. In addition, the studies of thermal properties of the quantum oscillator have been carried out in both commutative and NC space by several authors these years [29–34]. By employing the numerical method based on the Euler-MacLaurin formula, one can calculate the associated partition function and obtain the thermodynamic functions, such as the free energy, mean energy, heat capacity, and entropy. In spite of several papers that have been published concerning the thermodynamics properties of quantum oscillator, as far as we know they have not reported the case of DKP oscillator for spin 0 particle in NC space. It will be interesting to study thermodynamic properties of the spinless DKP oscillator in NC formalism.This work is organized as follows: In Section2 we briefly review the NC quantum mechanics, where the NC space is introduced. Subsequently, we analyze the solutions of the DKP oscillator for spin 0 particle in NC plane and obtain the corresponding energy spectrum. The thermodynamic properties of the NC DKP oscillator are investigated by employing the Euler-MacLaurin method in Section 3. To have an intuitive understanding for the statistical properties of the physical system, we depict the numerical results of the thermodynamic functions with several figures and discuss the effect of the NC parameter on thermal properties in Section 4. In addition, we also discuss the case of low temperature limit based on the Hurwitz zeta function formula. Finally, the work is summarized in last section.


## 2. The DKP Oscillator for Spin 0 Particle in NC Plane

### 2.1. Review of the NC Formalism

In NC space, the canonical variables satisfy the following commutation relations:(1)x^i,x^j=iθij,p^i,p^j=0,x^i,p^j=iħδij,i,j=1,2,with θij=εijθ being the antisymmetric NC parameter, representing the noncommutativity of the space, and x^i, p^i are the coordinate and momentum operators in the NC space.One way to deal with the problem of noncommutative space is via the star product or Moyal-Weyl product on the commutative space functions:(2)f∗gx=exp⁡i2θij∂ix∂jyfxgyy→x,where f(x) and g(y) are arbitrary infinitely differentiable functions.Then the time-independent Schrödinger equation in the NC space can be written in the form(3)Hx,p∗ψ=Eψ.The Moyal-Weyl product may be replaced by a Bopp shift of the form(4)xi⟶xi-12ħθijpj,pi⟶pi,where xi and pi are the position and momentum operators in the usual quantum mechanics which meet the canonical Heisenberg commutation relations.Then, via above Bopp shift, the Schrödinger equation can be written as(5)Hxi⟶xi-12ħθijpj,piψ=Eψ.Therefore, (5) is actually defined on the commutative space, and the NC effect may be calculated from the term that contains θ.


### 2.2. The Solution of the DKP Oscillator for Spin 0 Particle in NC Plane

The relativistic DKP equation for a free boson of massm is given by [35, 36](6)cβ→·P→+mc2ψ=iħβ0∂∂tψ,with β→=(β1,β2,β3) and β0 being the DKP matrices which satisfy the following algebra relation:(7)βμβνβλ+βλβνβμ=gμνβλ+gλνβμ,where gμν=diag⁡(1,-1,-1,-1) is the metric tensor in Minkowski space. For spin 0 particle, βμ are 5×5 matrices written as(8)β0=I00-0-T0,βi=0~ρι-ριT0,ι=1,2,3,with 0, 0-, and 0~ being 3×3, 2×3, and 2×2 zero matrices, respectively. Other matrices in (8) are given as follows:(9)I0=0110,ρ1=-100000,ρ2=0-10000,ρ3=00-1000.Introducing the nonminimal substitution (P→→P→-imωη0r→) in (6), we get the DKP oscillator as(10)cβ→·P→-imωη0r→+mc2ψ=iħβ0∂∂tψ,with ω being the oscillator frequency and η0=2(β0)2-1.It is easy to get the 2D DKP oscillator from the above equation:(11)cβ1px-imωη0x+cβ2py-imωη0y+mc2ψ=β0Eψ.Considering the NC formalism and via the Bopp shift (4), the DKP oscillator of (11) can be written as(12)cβ1px-imωη0x-θ2ħpy+cβ2py-imωη0y+θ2ħpx+mc2ψ=β0Eψ.For a boson of spin 0, the spinor ψ is a vector with five components [9, 10], which has the form ψ=ψ1,ψ2,iϕ1,iϕ2,iϕ3T.Substitutingψ into (12), one can obtain the following five coupled equations:(13)-mc2ψ1+Eψ2+icpx+imωx-θ2ħpyϕ3+icpy+imωy+θ2ħpxϕ4=0,Eψ1-mc2ψ2=0,cpx-imωx-θ2ħpyψ1+imc2ϕ3=0,cpy-imωy+θ2ħpxψ1+imc2ϕ4=0,imc2ϕ5=0.Combination of (13) gives(14)12m+mω2θ28ħ2px2+py2+mω22x2+y2-mω2θ2ħLzψ1=E2-m2c4+2mωħc22mc2ψ1.Using polar coordinates, the component ψ1 can be written as(15)ψ1ρ,φ=Rρexp⁡imlφ,ml=0,±1,±2,…,and substituting it into (14), we can obtain the radial equation(16)d2dρ2+1ρddρ-ml2ρ2-βρ2+γψ=0,with(17)β=4m2ω2m2ω2θ2+4ħ2,γ=E2+mlm2c2ω2θ-m2c4+2mωħc2ħ2+1/4m2ω2θ2c2.Introducing a new variable ξ=βρ2, one can get the following second-order differential equation:(18)d2Rξdξ2+1ξdRξdξ-ml24ξ2+14-γ4βξRξ=0.The equation can be analyzed asymptotically at the critical points ξ→0 and ξ→∞, allowing us to construct an equation whose form is the same as the confluent hypergeometric equation. This analysis shows that it is reasonable to write R(ξ) as(19)Rξ=ξml/2e-ξ/2Fξ.Substituting the above function into (18), one obtains(20)ξF′′ξ+ml+1-ξF′ξ-ml2+12-γ4βFξ=0.Obviously, this is confluent hypergeometric equation whose solution is the confluent hypergeometric function(21)Fξ=F11ml2+12-γ4β,ml+1;ξ.To obtain the solution of (20), the polynomial series should be convergent. It requires the independent term to be negative integer [37]:(22)ml2+12-γ4β=-n.Through a simple calculation we obtain the energy of the system as follows:(23)E2=2ħωmc21+m2ω2θ24ħ22n+ml+1-mlm2c2ω2θ+m2c4-2ħωmc2.The eigenvalues of considered physical system depend distinctly on the NC parameter θ which characterizes the noncommutativity of the space. We find that the energy spectrum is exact and not degenerate due to the NC effect.Thenψ1(ρ,φ) can be written in the form of the confluent hypergeometric function as(24)ψ1ρ,φ=Neimlφρmle-β/2ρ2F11-n,ml+1,βρ2,where β has already been given above, and through straightforward calculation, the expression of other components of the spinor ψ and the final form can be obtained easily.


## 2.1. Review of the NC Formalism

In NC space, the canonical variables satisfy the following commutation relations:(1)x^i,x^j=iθij,p^i,p^j=0,x^i,p^j=iħδij,i,j=1,2,with θij=εijθ being the antisymmetric NC parameter, representing the noncommutativity of the space, and x^i, p^i are the coordinate and momentum operators in the NC space.One way to deal with the problem of noncommutative space is via the star product or Moyal-Weyl product on the commutative space functions:(2)f∗gx=exp⁡i2θij∂ix∂jyfxgyy→x,where f(x) and g(y) are arbitrary infinitely differentiable functions.Then the time-independent Schrödinger equation in the NC space can be written in the form(3)Hx,p∗ψ=Eψ.The Moyal-Weyl product may be replaced by a Bopp shift of the form(4)xi⟶xi-12ħθijpj,pi⟶pi,where xi and pi are the position and momentum operators in the usual quantum mechanics which meet the canonical Heisenberg commutation relations.Then, via above Bopp shift, the Schrödinger equation can be written as(5)Hxi⟶xi-12ħθijpj,piψ=Eψ.Therefore, (5) is actually defined on the commutative space, and the NC effect may be calculated from the term that contains θ.


## 2.2. The Solution of the DKP Oscillator for Spin 0 Particle in NC Plane

The relativistic DKP equation for a free boson of massm is given by [35, 36](6)cβ→·P→+mc2ψ=iħβ0∂∂tψ,with β→=(β1,β2,β3) and β0 being the DKP matrices which satisfy the following algebra relation:(7)βμβνβλ+βλβνβμ=gμνβλ+gλνβμ,where gμν=diag⁡(1,-1,-1,-1) is the metric tensor in Minkowski space. For spin 0 particle, βμ are 5×5 matrices written as(8)β0=I00-0-T0,βi=0~ρι-ριT0,ι=1,2,3,with 0, 0-, and 0~ being 3×3, 2×3, and 2×2 zero matrices, respectively. Other matrices in (8) are given as follows:(9)I0=0110,ρ1=-100000,ρ2=0-10000,ρ3=00-1000.Introducing the nonminimal substitution (P→→P→-imωη0r→) in (6), we get the DKP oscillator as(10)cβ→·P→-imωη0r→+mc2ψ=iħβ0∂∂tψ,with ω being the oscillator frequency and η0=2(β0)2-1.It is easy to get the 2D DKP oscillator from the above equation:(11)cβ1px-imωη0x+cβ2py-imωη0y+mc2ψ=β0Eψ.Considering the NC formalism and via the Bopp shift (4), the DKP oscillator of (11) can be written as(12)cβ1px-imωη0x-θ2ħpy+cβ2py-imωη0y+θ2ħpx+mc2ψ=β0Eψ.For a boson of spin 0, the spinor ψ is a vector with five components [9, 10], which has the form ψ=ψ1,ψ2,iϕ1,iϕ2,iϕ3T.Substitutingψ into (12), one can obtain the following five coupled equations:(13)-mc2ψ1+Eψ2+icpx+imωx-θ2ħpyϕ3+icpy+imωy+θ2ħpxϕ4=0,Eψ1-mc2ψ2=0,cpx-imωx-θ2ħpyψ1+imc2ϕ3=0,cpy-imωy+θ2ħpxψ1+imc2ϕ4=0,imc2ϕ5=0.Combination of (13) gives(14)12m+mω2θ28ħ2px2+py2+mω22x2+y2-mω2θ2ħLzψ1=E2-m2c4+2mωħc22mc2ψ1.Using polar coordinates, the component ψ1 can be written as(15)ψ1ρ,φ=Rρexp⁡imlφ,ml=0,±1,±2,…,and substituting it into (14), we can obtain the radial equation(16)d2dρ2+1ρddρ-ml2ρ2-βρ2+γψ=0,with(17)β=4m2ω2m2ω2θ2+4ħ2,γ=E2+mlm2c2ω2θ-m2c4+2mωħc2ħ2+1/4m2ω2θ2c2.Introducing a new variable ξ=βρ2, one can get the following second-order differential equation:(18)d2Rξdξ2+1ξdRξdξ-ml24ξ2+14-γ4βξRξ=0.The equation can be analyzed asymptotically at the critical points ξ→0 and ξ→∞, allowing us to construct an equation whose form is the same as the confluent hypergeometric equation. This analysis shows that it is reasonable to write R(ξ) as(19)Rξ=ξml/2e-ξ/2Fξ.Substituting the above function into (18), one obtains(20)ξF′′ξ+ml+1-ξF′ξ-ml2+12-γ4βFξ=0.Obviously, this is confluent hypergeometric equation whose solution is the confluent hypergeometric function(21)Fξ=F11ml2+12-γ4β,ml+1;ξ.To obtain the solution of (20), the polynomial series should be convergent. It requires the independent term to be negative integer [37]:(22)ml2+12-γ4β=-n.Through a simple calculation we obtain the energy of the system as follows:(23)E2=2ħωmc21+m2ω2θ24ħ22n+ml+1-mlm2c2ω2θ+m2c4-2ħωmc2.The eigenvalues of considered physical system depend distinctly on the NC parameter θ which characterizes the noncommutativity of the space. We find that the energy spectrum is exact and not degenerate due to the NC effect.Thenψ1(ρ,φ) can be written in the form of the confluent hypergeometric function as(24)ψ1ρ,φ=Neimlφρmle-β/2ρ2F11-n,ml+1,βρ2,where β has already been given above, and through straightforward calculation, the expression of other components of the spinor ψ and the final form can be obtained easily.


## 3. Thermodynamic Properties of the DKP Oscillator for Spin 0 Particle in NC Plane

In this section, we study the thermodynamic properties of above DKP oscillator. The energy eigenvalues of the DKP oscillator in NC space are(25)En=±mc2αn+λ,n=0,1,2,…,with(26)α=4ħωmc21+m2ω2θ24ħ2,λ=2ħωmc21+m2ω2θ24ħ2ml+1-mlω2θc2-2ħωmc2+1.In Figure1, we plot the energy spectra E versus quantum number n for different values of the NC parameter θ′. Positive- and negative-energy levels correspond to the case of particle and antiparticle, respectively. We see that the energy E increases monotonically with quantum number n and the tendency of the spectrum can be observed for large quantum numbers. The effect of the NC parameter on the energy levels is observable, where θ′=0 corresponding to the case of commutative space.Figure 1
The profile of the energy spectraE of the DKP oscillator on NC plane versus quantum number n for different values of the NC parameter θ′.Because all thermodynamic quantities can be obtained from the partition functionZ, we first calculate the partition function of the system which is defined at temperature T as(27)ZV,T=∑n=0∞e-βEn-E0,where β=1/kBT, kB is the Boltzmann constant, and E0 is ground state energy. For N-body systems with no interaction inside, the corresponding partition function of the DKP oscillator should be ZN. Here for simplicity we just focus on the case of N=1.Then the thermodynamic properties of the physical system, such as free energy, mean energy, specific heat, and entropy, can be calculated from the following expressions [31]:(28)F=-1βln⁡Z,U=-∂ln⁡Z∂β,C=∂U∂T=-kBβ2∂U∂β,S=-∂F∂T=kBβ2∂F∂β.To perform our analysis on the thermodynamics of the DKP oscillator, we will restrict ourselves to stationary states of positive energy. The reasons are as follows [31, 38, 39].First, according to the study in [40, 41], the DKP equation possesses an exact Foldy-Wouthuysen transformation (FWT): thus, the positive- and negative-energy levels never mix. Second, although bosons can occupy the same level, which are different from the fermions controlled by the Pauli principle. We can still extend the method of filling the Dirac negative energy sea for the case of fermions to the one of bosons [42, 43]: as mentioned by the authors, the states of bosons are divided into two sectors: the usual positive sector consists of states with positive number of bosons, and the negative sector consists of states with negative number of bosons. Once a state enters the negative sector, it cannot return to the positive sector via an ordinary interaction due to the presence of a barrier.In addition, for the negative-energy states, the partition functionZ of the DKP oscillator at finite temperature T has the form(29)Z=∑n=0∞eβmc2αn+λ-λ=e-βmc2λ∑n=0∞eβmc2αn+λ.We apply the integral test to analyze the convergence of the above partition function. From (29), we can see that the function f(x)=eβmc2αx+λ is an increasing positive function, and the integral(30)∫0∞fxdx=∫0∞eβmc2αx+λdxis not convergent. This means that, according to the criterion of the integral test, the numerical partition function Z does not converge, while for the positive-energy states, it is obvious that the corresponding numerical partition function Z is convergent.According to the analysis above, we can assume that only particles with positive energy are available in order to set up a thermodynamic ensemble.Now let us go on the main object of our interest, substituting (25) into (27), we obtain the partition function of the DKP oscillator in NC plane as(31)Z=∑n=0∞e-βmc2αn+λ-λ=eβmc2λ∑n=0∞e-βmc2αn+λ.To calculate expression (31), we employ the Euler-MacLaurin formula(32)∑n=0∞fn=12f0+∫0∞fxdx-∑p=1∞12p!B2pf2p-10,with B2p being the Bernoulli numbers and(33)fx=e-βmc2αx+λ.Then we have(34)∫0∞fxdx=∫0∞e-βmc2αx+λdx=2αβ2m2c4βmc2λ+1e-βmc2λ,f10=-αβmc22λe-βmc2λ,f30=-α3βmc283λ-5/2+3βmc2λ-2+β2mc22λ-3/2e-βmc2λ.Up to B4, with B2=1/6, B4=-1/30, and combining with (34), we obtain the partition function at finite temperature T as(35)Z=12+2λ0α0τ+2α0τ2+α024λ0-α031920λ051τ-α031920λ021τ2-α035760λ031τ3,with(36)α0=4ξ1+θ′,λ0=2ξ1+θ′-2ξ+1,τ=1βmc2=kBTmc2,(37)ξ=ħωmc2,θ′=m2ω2θ24ħ2,where the dimensionless variable τ is proportional to temperature T and θ′ characterizes the noncommutativity of the space. From (35), we predict that the thermodynamic functions will be very complicated. To have an intuitive understanding for thermodynamic properties of the spinless DKP oscillator on NC plane, in what follows, we briefly depict our numerical results on the evaluation of the thermodynamic functions, that is, free energy, mean energy, specific heat, and entropy, via the numerical partition function Z.


## 4. Results and Discussions

In this section, all profiles of the thermodynamic quantities as a function of dimensionless temperature variableτ for different values of the NC parameter θ′, that is, θ′=0, 0.5, and 0.9, are plotted in Figures 2–6. In these figures, the natural unit (ħ=c=kB=1) is employed and we set m=1 and ξ=1.5.Figure 2
The partition functionZ of the DKP oscillator on NC plane as a function of τ for different values of the NC parameter θ′.Figure 3
The free energyF of the DKP oscillator in NC plane as a function of τ for different values of the NC parameter θ′.Figure 4
The mean energyU of the DKP oscillator in NC plane as a function of τ for different values of the NC parameter θ′.Figure 5
The heat capacityC of the DKP oscillator in NC plane as a function of τ for different values of the NC parameter θ′.Figure 6
The entropyS of the DKP oscillator in NC plane as a function of τ for different values of the NC parameter θ′.From the result shown in Figure2, we see that the partition function Z increases monotonically with dimensionless variable τ, and for a fixed value of τ the partition function decreases when NC parameter θ′ grows. It also shows that the partition function first has a slow-growth and then rapidly increases. The Helmholtz free energy F is shown in Figure 3. We see that in all cases the profile of the curves decreases monotonically with the temperature and the free energy increase with NC parameter θ′ growing for a fixed value of τ. In Figure 4, we plot the mean energy U versus τ for different values of the NC parameter θ′, and it shows that all the curves have similar linear behavior and display very close profiles. To have better understanding of the NC effect on mean energy, we choose θ′=0, 0.5, and 0.9. And we also see that for a fixed τ the mean energy decreases when θ′ grows. The profile of heat capacity C as a function of τ for different NC parameter θ′ values is depicted in Figure 5. It is observed that the heat capacity increases for increasing τ at first and then decreases with τ growing. There is an abrupt change around τ0 for each curve, and τ0 increases when the corresponding θ′ grows. Figure 6 displays the curves of the numerical entropy S versus τ for different values of the NC parameter θ′. It shows that the tendency of entropy rapidly increases at first for increasing τ and then slowly grows for large τ values. For a fixed value of τ, the entropy decreases when NC parameter θ′ grows. In addition, from the curves of the numerical mean energy, free energy, and entropy function, no abrupt changes around τ0 have been identified, meaning that the curvature, observed in the specific heat curve, does not exhibit or indicate an existence of a phase transition around a finite temperature τ0.Because the dimensionless variableτ is proportional to temperature T, we can get the characters of these thermodynamic functions at high temperature from the asymptotic behavior shown in the above figures. As is shown, the values of mean energy U and entropy S are large at high temperature, because these functions increase monotonically with dimensionless variable τ. Contrarily, free energy F will become very small at high temperature because it decreases monotonically with τ growing. The heat capacity C will tend to a constant value at high temperature. The effect of the NC parameter θ′ on thermodynamic properties at high temperature is observable. We see that the values of mean energy U and entropy S at high temperature case decrease with NC parameter θ′ growing, and in contrast the value of free energy F increases when θ′ grows, while from Figure 5, we found that all the curves of heat capacity against different NC parameter θ′ at high temperature would tend to same value. It means that the effect of noncommutativity on the heat capacity can be negligible at high temperature.Recently, Hassanabadi et al. have studied the spin 0 scalar bosons in presence of a uniform magnetic field by using the Klein-Gordon equation in NC space [33]. Similarly, with the numerical method based on the Euler-MacLaurin formula, the authors calculate the associated partition function and obtain the thermodynamic functions such as the free energy, mean energy, heat capacity, and entropy. Now we will compare our results with those obtained by them. For the free energy, mean energy, and entropy, it is easy to find that our result is consistent with theirs, and the curves have similar profiles, while the curve of heat capacity seems a little different: our result about the heat capacity will tend to a constant value C0 at high temperature compared with the result shown in [33]. This divergence may be result of the different range of variable τ.In the above work, we have studied the thermodynamic properties of the NC DKP oscillator by employing the Euler-MacLaurin approximation. According to the literature [38, 39], the formalism is only valid for high temperature because the partition function diverges in the low temperature limit τ→0. In order to cover all range of temperature, one can use a method based on the Hurwitz zeta function which has been recently proposed by Boumali. In the following we will discuss the case of low temperature limit by using zeta function method.Given the energy spectrum (25), one can define the partition via [38, 39](38)Z=∑n=0∞e-βEn=∑n=0∞e-γ/τn+κ,with the following substitutions:(39)γ=α0,κ=λ0α0,where α0, λ0, and τ are given in (36), using the formula [44](40)e-x=12πi∫cdsx-sΓs.Then (38) is transformed into(41)∑n=0∞e-γ/τn+κ=12πi∫cdsγτ-s∑nn+κ-s/2Γs=12πi∫cdsγτ-sζHs2,κΓs,where x=(γ/τ)n+κ and Γ(s) and ζHs/2,κ are the Euler and Hurwitz zeta function, respectively. Employing the residues theorem, for the two poles s=0 and s=2, the associated partition function is written in terms of the Hurwitz zeta function as follows:(42)Zτ=τγ2+ζH0,κ.Now, using the expression(43)ζH0,κ=12-κ,the final partition function becomes(44)Zτ=τ2α0+12-λ0α0.Then, all thermal quantities of the DKP oscillators, such as the free energy F, mean energy U, specific heat C, and entropy S, can be obtained by using the following expressions derived from (28) as(45)F=-mc2τln⁡Z,U=mc2τ2∂ln⁡Z∂τ,S=kBln⁡Z+kBτ∂ln⁡Z∂τ,C=2kBτ∂ln⁡Z∂τ+kBτ2∂2ln⁡Z∂τ2.From expression (44), we see that the partition function is not divergent at the low temperature limit τ→0; thus the result can be considered as an improvement of the result obtained by using the Euler-MacLaurin method.


## 5. Conclusions

This paper was devoted to study of the thermodynamic properties of the DKP oscillator for spin 0 particle in noncommutative space. We first analyzed the DKP equation in the presence of the Dirac oscillator interaction in NC space and obtained energy spectrum as well as corresponding spinorψ of the system by employing the functional analysis method. Subsequently, we investigated the thermodynamic properties of the system by employing the Euler-MacLaurin method and depicted our numerical results on the evaluation of the thermodynamic functions through the associated partition function Z. And the effect of the NC parameter on thermodynamic properties was discussed. We show that the energy spectrum and corresponding thermodynamic functions depend explicitly on the NC parameter θ which characterizes the noncommutativity of the space. It is also observed that the influence of NC space on the thermodynamic properties is nonignorable although the NC parameter is small. In addition, we discuss the case of low temperature limit based on the Hurwitz zeta function formula, which can be considered as an improvement of the Euler-MacLaurin method. This work will allow us to further explore the role of NC space in the related physical concepts and deepen our knowledge on NC quantum mechanics.



---

*Source: 901675-2015-10-22.xml*