The Dirac oscillator in a spinning cosmic string spacetime

TL;DR

This paper investigates how a Dirac oscillator behaves in a spinning cosmic string spacetime, analyzing energy levels and eigenfunctions affected by gravitational and topological parameters.

Contribution

It provides a novel analysis of the Dirac oscillator in a spinning cosmic string background, deriving eigenfunctions and energy spectra considering various physical parameters.

Findings

Eigenfunctions and energy levels are derived for the Dirac oscillator in this spacetime.

The energy spectrum depends on the string density, rotation frequency, and oscillator parameters.

The study highlights the influence of topological defects on relativistic quantum systems.

Abstract

We examine the effects of gravitational fields produced by topological defects on a Dirac field and a Dirac oscillator in a spinning cosmic string spacetime. We obtain the eigenfunctions and the energy levels of the relativistic field in that background and consider the effect of various parameters, such as the frequency of the rotating frame, the oscillator's frequency, the string density and other quantum numbers.

Equations

$$d{s}^2=-d{T}^2+d{X}^2+d{Y}^2+d{Z}^2$$

$$\begin{array}{c}T=t+a{\alpha }^{-1}\phi ,\text{hfill}\\ X=r\cos \phi ,\text{hfill}\\ Y=r\sin \phi ,\text{hfill}\\ \phi =\alpha {\phi }^{\prime },\text{hfill}\end{array}$$

$$d{s}^2$$

$$(\mathrm{i}{\gamma }^{\mu }\left(x\right){\nabla }_{\mu }-M)\Psi \left(x\right)=0,$$

$${\nabla }_{\mu }={\partial }_{\mu }+{\Gamma }_{\mu }\left(x\right),$$

$${\Gamma }_{\mu }=\frac{1}{2}{\omega }_{\mu \bar{a}\bar{b}}\left[{\gamma }^{\bar{a}},{\gamma }^{\bar{b}}\right],$$

$${\omega }_{\mu \bar{a}\bar{b}}={\eta }_{\bar{a}\bar{c}}{e}_{\nu }^{\bar{c}}{e}_{\bar{b}}^{\sigma }{\Gamma }_{\sigma \mu }^{\nu }-{\eta }_{\bar{a}\bar{c}}{e}_{\bar{b}}^{\nu }{e}_{\nu }^{\bar{c}}.$$

$${\Gamma }_{ij}^{\mu }=\frac{1}{2}{g}^{\mu k}\left(\frac{\partial g_{ik}}{\partial q_j}+\frac{\partial g_{jk}}{\partial q_i}-\frac{\partial g_{ij}}{\partial q_k}\right),$$

$${\Gamma }_{r\phi }^t={\Gamma }_{\phi r}^t=\frac{-a}{r},{\Gamma }_{\phi \phi }^r=-r{\alpha }^2,{\Gamma }_{r\phi }^{\phi }={\Gamma }_{\phi r}^{\phi }=\frac{1}{r}.$$

$$\displaystyle e_{\mu}^{{\bar{a}}}=\left({\begin{array}[]{*{20}{c}}1&0&a&0\\ 0&{\cos\varphi}&{-r\alpha{\sin\varphi}}&0\\ 0&{\sin\varphi}&{-r\alpha{\cos\varphi}}&0\\ 0&0&0&1\end{array}}\right),\quad e_{{{\bar{a}}}}^{\mu}=\left({\begin{array}[]{*{20}{c}}1&{\frac{a{\sin\varphi}}{{r\alpha}}}&{\frac{-a{\cos\varphi}}{{r\alpha}}}&0\\ 0&{\cos\varphi}&{\sin\varphi}&0\\ 0&{\frac{-{\sin\varphi}}{{r\alpha}}}&{\frac{\cos\varphi}{{r\alpha}}}&0\\ 0&0&0&1\end{array}}\right).$$

$$\displaystyle\begin{array}[]{l}e_{{{\bar{a}}}}^{\mu}\left(x\right)e_{\nu}^{{{\bar{a}}}}\left(x\right)=\delta_{\nu}^{\mu},\\ e_{\mu}^{{{\bar{a}}}}\left(x\right)e_{{{\bar{b}}}}^{\mu}\left(x\right)=\delta{{}_{\bar{b}}^{\bar{a}}},\end{array}$$

$$g_{\mu \nu }\left(x\right)=e_{\mu }^{\bar{a}}\left(x\right)e_{\nu }^{\bar{b}}\left(x\right){\eta }_{{}^{\bar{a}\bar{b}}}.$$

$${\omega }_{\phi }^{\bar{1}\bar{2}}=-{\omega }_{\phi }^{\bar{2}\bar{1}}=1-\alpha ,$$

$${\Gamma }_{\phi }=-\frac{1}{2}(1-\alpha ){\gamma }^{\bar{1}}{\gamma }^{\bar{2}}.$$

$${\gamma }^{\mu }\left(x\right)=e_{\bar{a}}^{\mu }{\gamma }^{\bar{a}}.$$

$$\displaystyle{\gamma^{{{\bar{0}}}}}=\left({\begin{array}[]{*{20}{c}}1&0\\ 0&{-1}\end{array}}\right),{{\gamma^{{\bar{i}}}}=\left({\begin{array}[]{*{20}{c}}0&{{\sigma^{i}}}\\ {-{\sigma^{i}}}&0\end{array}}\right)},\qquad{i=1,2,3,}$$

$$\sigma^{1}=\left(\begin{array}[]{cc}0&1\\ 1&0\end{array}\right),\quad\sigma^{2}=\left(\begin{array}[]{cc}0&-i\\ i&0\end{array}\right),\quad\sigma^{3}=\left(\begin{array}[]{cc}1&0\\ 0&-1\end{array}\right).\quad$$

$$\begin{array}{c}{\gamma }^r={\gamma }^1=e_{\bar{a}}^1{\gamma }^{\bar{a}}={\gamma }^{\bar{1}}\cos \phi +{\gamma }^{\bar{2}}\sin \phi ,\text{hfill}\\ {\gamma }^z=e_{\bar{a}}^z{\gamma }^{\bar{a}}={\gamma }^3,\text{hfill}\\ {\gamma }^2=e_{\bar{a}}^2{\gamma }^{\bar{a}}=\frac{{\gamma }^{\phi }}{r\alpha },\text{hfill}\\ {\gamma }^{\phi }=e_{\bar{0}}^{\phi }{\gamma }^{\bar{0}}+e_{\bar{1}}^{\phi }{\gamma }^{\bar{1}}+e_{\bar{2}}^{\phi }{\gamma }^{\bar{2}}=-{\gamma }^{\bar{1}}\sin \phi +{\gamma }^{\bar{2}}\cos \phi ,\text{hfill}\\ {\gamma }^t={\gamma }^0=e_{\bar{a}}^0{\gamma }^{\bar{a}}={\gamma }^{\bar{0}}+\frac{a\sin \phi }{r\alpha }{\gamma }^{\bar{1}}-\frac{a\cos \phi }{r\alpha }{\gamma }^{\bar{2}}={\gamma }^{\bar{0}}-\frac{a}{r\alpha }{\gamma }^{\phi }.\text{hfill}\end{array}$$

$${\Gamma }_{\phi }=i\frac{(1-\alpha )}{2}s{\sigma }^z.$$

$$\left(i{\gamma }^t{\partial }_t+i\vec{\gamma }\cdot \nabla +i{\gamma }^{\mu }{\Gamma }_{\mu }\right)\Psi =\left[i\left({\gamma }^{\bar{0}}-\frac{a}{r\alpha }{\gamma }^{\phi }\right){\partial }_t+i\vec{\gamma }\cdot \nabla +i{\gamma }^{\mu }{\Gamma }_{\mu }\right]\Psi =M\Psi$$

$$\left[i\left({\gamma }^{\bar{0}}{\gamma }^{\bar{0}}-\frac{a}{r\alpha }{\gamma }^{\bar{0}}{\gamma }^{\phi }\right){\partial }_t-{\gamma }^{\bar{0}}M\right]\Psi ={\gamma }^{\bar{0}}\vec{\gamma }\cdot \left[-i{\nabla }_{\alpha }-i\vec{\Gamma }\right]\Psi .$$

$$\Psi (t,r,\phi )=e^{-i(Et-m\phi )}\Psi (r),$$

$$\left(E-{\gamma }^{\bar{0}}M\right)\Psi ={\gamma }^{\bar{0}}\vec{\gamma }\cdot \left(-i{\nabla }_{\alpha }-i\vec{\Gamma }\right)\Psi +\frac{aE}{r\alpha }{\alpha }^{\bar{2}}\Psi .$$

$$\displaystyle\Psi(r)=\left({\begin{array}[]{*{20}{c}}{\phi(r)}\\ {{\chi}(r)}\end{array}}\right),$$

$${\gamma }^{\bar{0}}\vec{\gamma }\cdot \hat{\phi }={\alpha }^{\phi },$$

$$\vec{\pi }=-i{\nabla }_{\alpha }-i\vec{\Gamma }+\frac{aE}{\alpha r}\hat{\phi },$$

$$\begin{array}{c}\left(E-M\right)\varphi (r)=\vec{\sigma }\cdot \vec{\pi }\chi (r),\text{hfill}\\ \left(E+M\right)\chi (r)=\vec{\sigma }\cdot \vec{\pi }\varphi (r).\text{hfill}\end{array}$$

$$\left(E^2-M^2\right)\varphi =\left(\vec{\sigma }\cdot \vec{\pi }\right)\left(\vec{\sigma }\cdot \vec{\pi }\right)\varphi =\vec{\pi }\cdot \vec{\pi }\varphi .$$

$$\vec{\pi }=-i{\nabla }_{\alpha }+\left(\frac{1}{2\alpha r}(1-\alpha )s+\frac{aE}{\alpha r}\right)\hat{\phi }=:-i{\nabla }_{\alpha }+\frac{\xi }{r}\hat{\phi },$$

$$\xi \equiv \frac{1}{2\alpha }(1-\alpha )s+\frac{aE}{\alpha }.$$

Full text

**Key Words: Dirac oscillator, cosmic strings **

PACS: 03.65.Pm, 04.62. +v, 11.27.+d

1 Introduction

The Dirac equation with interactions linear in the coordinates was initially studied in Refs. [1, 2, 3]. Such a system was referred to as a ‘Dirac oscillator’ in Ref. [4], because, in the non-relativistic limit, it behaves as a harmonic oscillator with a strong spin-orbit coupling term. This model describes the dynamics of a harmonic oscillator for spin-half particles and is obtained by introducing a non-minimal prescription into the free Dirac equation [4]. It was observed that the Dirac oscillator interaction is a physical system which can be interpreted as the interaction of the anomalous magnetic moment with a linear electric field [5, 6]. The authors of Ref. [5] also established the conformal invariance of the Dirac oscillator and the authors of Ref. [6] examined its covariance properties, applied the Foldy-Wouthuysen and the Cini-Toushek transformations. As a relativistic quantum mechanical system, the Dirac oscillator has been widely studied and, because it is an exactly solvable model, several investigations have been performed in its theoretical framework.

Although the Dirac oscillator is normally utilized within the context of many-body theory, relativistic quantum mechanics and quantum chromodynamics (such as the interaction between quarks as well as the confining part of the phenomenological Cornell potential), the Dirac oscillator and related models have been applied in many other contexts as well, such as quantum optics [7, 8, 9], supersymmetry [5, 10, 11], nuclear reactions [12], the hadronic spectrum (with the two-body Dirac oscillator) [13, 14], the Clifford algebra [15, 16], non-commutative space[17, 18], thermodynamic properties [19], Lie algebras [20], supersymmetric (non-relativistic) quantum mechanics [21], the supersymmetric path-integral formalism [22], chiral phase transitions in presence of a constant magnetic field [8], the relativistic Landau levels in presence of external magnetic field [23], the Aharonov-Bohm effect[24], and condensed matter physics phenomena and graphene [25]. Similar studies for the Duffin-Kemmer-Petiau (DKP) oscillator, which is analogous to the Dirac oscillator for spinless and spin-one particles, are in Refs. [26, 27]. Finally, let us mention many studies of the Dirac oscillator with topological defects and cosmic string spacetimes in Refs. [28, 29, 30, 31, 32, 33, 34] and analogous investigations for scalar fields in Refs. [35, 36, 37, 38, 39]. Some studies of relativistic oscillators are in Refs. [40, 41, 42, 43, 44, 45, 46, 47].

In this work, we examine the relativistic quantum dynamics of Dirac oscillator on the curved spacetime of a rotating cosmic string. From the corresponding Dirac equation, we analyze the influence of the topological defect on the equation of motion, the energy spectrum and the wave-functions. An analogous study for the Klein-Gordon equation is in Ref. [49]. In Sec. 2, we write down the covariant Dirac equation without oscillator in a spinning cosmic string spacetime, and find its wave-functions and energy eigenvalues. In Sec. 3, we present the covariant Dirac oscillator in the same spacetime and obtain the wave-functions and energy spectrum. We present concluding remarks in Sec. 4.

2 Dirac equation in the cosmic string spacetime

The spacetime generated by a spinning cosmic string without internal structure, or ‘ideal’ spinning cosmic string, can be obtained from the line element

[TABLE]

by applying the coordinate transformation

[TABLE]

which leads to (see also Refs. [50, 51, 52, 53])

[TABLE]

where $-\infty<z<\infty$ , $r\geq 0$ and $0\leq\varphi\leq 2\pi$. We work with units such that $c=1$. The angular parameter $\alpha$ runs in the interval $\left({0,1}\right]$ and is related to the linear mass density $\mu$ of the string by $\alpha=1-4\mu$, and to the deficit angle by $\gamma=2\pi(1-\alpha)$. We have also $a=4Gj$ where $G$ is the universal gravitation constant and $j$ is the angular momentum of the spinning string; thus $a$ is a length that represents the rotation of the cosmic string. Note that in this case, the source of the gravitational field of a spinning cosmic string possesses angular momentum and the metric (8) has an off-diagonal term involving time and space.

The Dirac equation for a field $\Psi$ of mass $M$ in the cosmic string spacetime described by Eq. (8) reads[54, 55, 56, 57]

[TABLE]

with the covariant derivative

[TABLE]

and the spinorial affine connections,

[TABLE]

where ${{\gamma^{\bar{a}}}}$ denotes the standard Dirac matrices in Minkowski spacetime with metric $\eta_{\bar{a}\bar{b}}=\left(-1,+1,+1,+1\right)$, and $\omega_{\mu\bar{a}\bar{b}}$ is the spin connection, given by

[TABLE]

We use greek indices $\mu$, $\nu$, etc. for the curved spacetime, and bar latin indices $\bar{a}$, $\bar{b}$, etc. in Minkowski spacetime.

As discussed in Refs. [58, 59], the spin connection allows us to construct a local frame through the tetrad basis which gives the spinors in the curved spacetime. The Christoffel symbols of the second kind, $\Gamma_{\mu\nu}^{\sigma}$, can be obtained from

[TABLE]

with coordinates $(q_{1},q_{2},q_{3})$. With the metric in Eq. (8), the non-null components of Christoffel symbols are

[TABLE]

We can build the local reference frame through a non-coordinate basis with components $e_{\mu}^{{\bar{a}}}$ called tetrads or vierbeins which form our local reference frame. With the line element (8), we can use tetrads $e_{{\bar{a}}}^{\mu}$ and $e_{\mu}^{{\bar{a}}}$ (obtained in Ref. [52]) as follows

[TABLE]

The vierbeins satisfy the orthonormality conditions

[TABLE]

and

[TABLE]

The non-null components of the spin connection are

[TABLE]

and the only non-vanishing spinorial affine connection is

[TABLE]

The generalized Dirac matrices $\gamma^{\mu}$ in curved spacetime are related to their Minkowski counterparts via

[TABLE]

In terms of the Minkowski flat spacetime coordinates, these matrices can be cast into the form

[TABLE]

where ${{\sigma^{i}}}$ are the Pauli matrices:

[TABLE]

The matrices ${\gamma^{\mu}}\left(x\right)$ in Eq. (11) are (see Ref. [52])

[TABLE]

The $z$-translation symmetry of Eq. (8) allows us to reduce the four-component Dirac equation (7) to two two-component spinor equations. We shall also take $p_{z}=0$ and $z=0$. We find that Eq. (29) yields

[TABLE]

If we substitute the curved-spacetime gamma matrices of Eq. (43) into Eq. (9), we find

[TABLE]

Since we work in the plane $z=0$, we effectively have ${\vec{\gamma}}=\left(\gamma^{r},\gamma^{\varphi}\right)$ and $\nabla\rightarrow\nabla_{\alpha}\equiv{\hat{r}}\frac{\partial}{\partial r}+\frac{\hat{\varphi}}{\alpha r}\frac{\partial}{\partial\varphi}$. Since Eq. (29) is the only non-zero component of $\Gamma_{\mu}$, we write $\gamma^{\mu}\Gamma_{\mu}={\vec{\gamma}}\cdot{\vec{\Gamma}}$, with ${\vec{\Gamma}}=(\Gamma_{r},\Gamma_{\varphi})=(0,\Gamma_{\varphi})$. Next we multiply Eq. (45) with $\gamma^{\bar{0}}$,

[TABLE]

From Eq. (35), we see that ${\gamma^{\bar{0}}}{\gamma^{\bar{0}}}=1$ and we define \alpha^{i}\equiv{\gamma^{\bar{0}}}{\gamma^{i}}=\left(\begin{array}[]{cc}0&{{\sigma^{i}}}\\ {{\sigma^{i}}}&0\end{array}\right). Since the interaction is time-independent, one can write

[TABLE]

where $E$ is the energy of the fermion, so that $i\partial_{t}\Psi=E\Psi$, and $\partial_{\varphi}\Psi=im\Psi$, where the quantum number $m$ is related to the $z$-component of the total angular momentum operator, with half-integer values (see e.g. Eq. (15) in Ref. [28] and Ref. [52]). Then Eq. (46) becomes

[TABLE]

We write the radial wave-function as

[TABLE]

where $\phi(r)$ and $\chi(r)$ are two-spinors. If we write

[TABLE]

so that $\frac{aE}{{r\alpha}}\alpha^{\bar{2}}=\gamma^{\bar{0}}{\vec{\gamma}}\cdot\left(\frac{aE}{{r\alpha}}{\hat{\varphi}}\right)$, and define the generalized momentum $\vec{\pi}$ as follows,

[TABLE]

then Eq. (48) leads to the coupled equations

[TABLE]

We substitute the second equation of Eq. (56) into the first and obtain

[TABLE]

With $\Gamma_{\varphi}$ from Eq. (44), we find that the explicit form of Eq. (52) is

[TABLE]

where

[TABLE]

We find

[TABLE]

When we substitute this into Eq. (57), use Eq. (47) and $\nabla_{\alpha}^{2}=\frac{\partial^{2}}{\partial r^{2}}+\frac{1}{r}\frac{\partial}{\partial r}+\frac{1}{\alpha^{2}r^{2}}\frac{\partial^{2}}{\partial\varphi^{2}},$ we obtain

[TABLE]

Note that $\phi(r)$ (and $\chi(r)$) is an eigenfunction of $\sigma^{3}$ with eigenvalues $\pm 1$, so we can write $\phi_{s}=\left(\phi_{+},\phi_{-}\right)^{T}$ with $\sigma^{3}\phi_{s}=s\phi_{s}$, $s=\pm 1$. Therefore the action of $\xi$ on $\phi_{s}$ is the reason why we replaced $\sigma^{3}$ by $s$ in Eq. (58). Since, in Eq. (60), we have $\left(\frac{m^{2}}{\alpha^{2}r^{2}}+\frac{\xi^{2}}{r^{2}}+2m\frac{\xi}{\alpha r^{2}}\right)=\frac{1}{r^{2}}\left(\frac{m}{\alpha}+\xi\right)^{2}$, we can express Eq. (60) as

[TABLE]

where

[TABLE]

By performing the change of variable $r=\frac{x}{\kappa}$, we can write Eq. (61) in the form

[TABLE]

The physical solution of Eq. (63) is given by the Bessel function of the first kind, $J_{\eta}(x)$ so that the radial wave-function is given by

[TABLE]

We will not discuss these solutions any further, because we will discuss in detail the more general Dirac oscillator in the next section.

3 The Dirac oscillator in cosmic string background

In this section, we turn on the Dirac oscillator interaction by adding the non-minimal substitution term,

[TABLE]

where $\beta=\gamma^{0}$, $\vec{r}=\hat{r}r$ is the position vector, and $\omega$ is the oscillator’s frequency. (In the presence of a constant magnetic field $B_{0}$, $\omega=\frac{eB_{0}}{2M}$ is the so-called cyclotron frequency of the oscillator.) We substitute Eq. (65) into Eq. (9) and obtain

[TABLE]

where ${\vec{\pi}}$ is given in Eq. (52) and $\Psi(r)$ in Eq. (51). We obtain the counterpart of Eq. (56):

[TABLE]

From these two equations, we obtain

[TABLE]

Next we utilize $\left(\vec{\sigma}\cdot\vec{a}\right)\left(\vec{\sigma}\cdot\vec{b}\right)=\vec{a}\cdot\vec{b}+i\vec{\sigma}\cdot\left(\vec{a}\times\vec{b}\right)$ with $\vec{a}=\vec{\pi}+iM\omega\hat{r}r$ and $\vec{b}=\vec{\pi}-iM\omega\hat{r}r$, and note that $\vec{a}\cdot\vec{b}=\vec{\pi}\cdot\vec{\pi}-M^{2}\omega^{2}r^{2}$ (with $\vec{\pi}\cdot\vec{\pi}\phi(r)$ in Eq. (59)) and $\vec{a}\times\vec{b}=iM\omega r\left(\hat{r}\times\vec{\pi}-\vec{\pi}\times\hat{r}\right)$. We find $\left(\vec{\pi}\times\vec{r}-\vec{r}\times\vec{\pi}\right)\phi=-2\hat{z}\left(\frac{m}{\alpha}+\xi\right)\phi$ and $\left(\vec{\pi}\cdot\vec{r}-\vec{r}\cdot\vec{\pi}\right)\phi=-2i\phi$, so that Eq. (71) leads to an equation similar to Eq. (61),

[TABLE]

with $\eta$ as in Eq. (62) and

[TABLE]

As $\omega$ approaches zero, Eq. (72) reduces to Eq. (61) with the only extra terms $-M^{2}\omega^{2}r^{2}$, and with Eq. (73) which reduces to Eq. (62). With the change of variable $M\omega{r^{2}}=\rho$, then we rewrite the radial equation (72) in the form

[TABLE]

Let us write the solution of Eq. (74) in terms of the new function $F(\rho)$ defined as

[TABLE]

so that Eq. (74) can be rewritten as

[TABLE]

The general solution for Eq. (75) is the confluent hypergeometric function

[TABLE]

Even though Eq. (72) clearly reduces to Eq. (61) when $\omega\rightarrow 0$, as noted above, the corresponding reduction from the solution in Eq. (77) (when expressed in terms of $r$) to Eq. (64) is more subtle and we will not discuss it in detail. This is done by rescaling the solution of Eq. (72) by a function of $M\omega$ (the factor of $r^{2}$ in the extra term $-M^{2}\omega^{2}r^{2}$ of Eq. (72)) so that it is still a solution and such that the limit $\omega\rightarrow 0$ exists and is non-zero. Hypergeometric functions allow for Taylor expansion and from this one can determine which expansion coefficients vanish or diverge in the limit. For divergent terms, one rescales by an appropriate power in $M\omega$ so that the limit exists and leads to the Bessel function of Eq. (64).

Because of the divergent behaviour of the function $F(\rho)$ for large values of $\rho$, bound states solutions can only be obtained by imposing that this function becomes a polynomial of some degree $n$. Then the radial solution presents an acceptable behaviour at infinity. This condition is obtained by setting

[TABLE]

By substituting Eqs. (58), (62) and (73) into Eq. (78), we obtain

[TABLE]

If the argument of the absolute value is positive, $\frac{m}{\alpha}+\frac{1}{2\alpha}\left(1-\alpha\right)s+\frac{aE}{\alpha}>0$, then we find

[TABLE]

Note that for $s=1$, the eigenvalues in Eq. (80),

[TABLE]

are independent of $m$, $a$, and $\alpha$. Fig. 1 displays $E_{n}$ from Eq. (81) as a function of $\omega$ for $M=1$ and various values of $n$. As per Eq. (81), Fig. 1 shows how $\left|E_{n}\right|$ increases when $n$ increases, for any $\omega$. Eq. (81) and Fig. 1 show that for small values of $\omega$, $E_{n}\rightarrow\pm M$ for any $n$. Eq. (80) shows that $E_{n,m}$ depends on $n$ and the spectrum of energy is discrete.

For $s=-1$, the expression for $E_{n,m}$ contains all the terms of Eq. (80) and depends on all its parameters, including $m$, $a$ and $\alpha$. The positive values of the energy spectrum with $s=-1$ is displayed as a function of $\alpha$ (with $M=1$, $m=1/2$, $a=0.1$ and $\omega=1$, for $n=1$ and $2$, in Fig. 2, and as a function of $\omega$ (with $M=1$, $m=1/2$, $a=0.1$ and $\alpha=0.2$) in Fig. 3. In both figures, we see that $\left|E\right|$ increases with $n$ for any value of $\alpha$ and $\omega$. For $a=0$, the energy spectrum of Eq. (79) is given by

[TABLE]

which is analogous to Eq. (28) in Ref. [30] (we have considered the two-dimensional oscillator so that their $k=0$).

When the argument of the absolute value is negative, $\frac{m}{\alpha}+\frac{1}{2\alpha}\left(1-\alpha\right)s+\frac{aE}{\alpha}<0$, we obtain, rather than Eq. (80),

[TABLE]

Then, for $s=-1$, the eigenvalues in Eq. (83) are independent of $m$, $a$, and $\alpha$, as in Eq. (81). However, when $s=1$, the energy depends explicitly on all the parameters in Eq. (80), including $m$, $a$ and $\alpha$. Whether $s=-1$ or $+1$, $E$ depends on $n$ and the spectrum of energy is discrete.

4 Conclusion

We have studied the relativistic quantum dynamics of a Dirac oscillator field subject to a linear interaction for spin-half particles in a cosmic string spacetime. The Dirac oscillator is a kind of tensor coupling with a linear potential which leads to the harmonic oscillator problem in the weak-coupling limit. This can be interpreted as the interaction of the anomalous magnetic moment with a linear electric field. From the corresponding Dirac oscillator equation, we analyzed the influence of the topological defect on the equation of motion, the energy spectrum and the wave-function. Firstly, we solved the Dirac equation without the oscillating term. We obtained the Bessel function as a solution of the resulting equation.

Next we turned to the Dirac oscillator, which is described by introducing a non-minimal coupling within the previous equations. The new term has direct implications on the energy spectrum and wave-functions of the oscillator. The wave-function, which is a Bessel function without the Dirac oscillator, becomes, with the Dirac oscillator, related to the confluent hypergeometric function in Eq. (77), and depends on the parameters $\eta$ (from Eq. (62)) and $\kappa_{\omega}$ (in Eq. (73)), thus depending on the parameters $m$, $\alpha$, $s$, $M$, $\omega$, and the energy $E$. In turn, we found that the energy eigenvalues $E$ depend on $m$, $\alpha$, $s$, $M$, $\omega$. Figs. 1 and 3 show that $\left|E\right|$ increases with $\omega$ for specific values of other parameters, and that $\left|E\right|$ increases with $n$. Fig. 2 shows that, for specific values of other parameters, $\left|E\right|$ decreases with $\alpha$ and increases with $n$. The angular parameter $\alpha$ is related to the linear mass density of the cosmic string. In the limit $\alpha=1$, that is, in the absence of topological defect, we recover known results for the flat spacetime. Here we have showed that the geometric and topological properties of these spacetimes lead to shifts in the energy spectrum and the wave-function that is comparable with the flat Minkowski spacetime.

Acknowledgement

The authors thank Thomas Creutzig for helpful comments and the anonymous referee for helpful comments and suggestions on the manuscript. M. de Montigny acknowledges the Natural Sciences and Engineering Research Council (NSERC) of Canada for partial financial support (grant number RGPIN-2016-04309).