Effects of spin on the dynamics of the 2D Dirac oscillator in the magnetic cosmic string background

In this work the dynamics of a 2D Dirac oscillator in the spacetime of a magnetic cosmic string is considered. It is shown that earlier approaches to this problem have neglected a $\delta$ function contribution to the full Hamiltonian, which comes from the Zeeman interaction. The inclusion of spin effects leads to results which confirm a modified dynamics. Based on the self-adjoint extension method, we determined the most relevant physical quantities, such as energy spectrum, wave functions and the self-adjoint extension parameter by applying boundary conditions allowed by the system.


I. INTRODUCTION
The Dirac oscillator is a natural model for studying properties of physical systems. This model is based on the dynamics of a harmonic oscillator for spin-1/2 particles by introducing a nominimal prescription into free Dirac equation [1]. Because it is a exactly solvable model, several investigations have been developed in the context of this theoretical framework in the latest years. The interest in this issue appears in different contexts, such as quantum optics [2][3][4], supersymmetry [5][6][7], nuclear reactions [8], Clifford algebra [9,10] and noncommutative space [11,12]. Recently, the Dirac oscillator has been verified experimentally by J. A. Franco-Villafañe et al., based on a tight-binding system [13]. A detailed description for the Dirac oscillator is given in Ref. [14] and for other contributions see Refs. [15][16][17][18][19][20][21][22].
Among the various contexts in which the Dirac oscillator can be addressed, we refer to the cosmic string, a linear defect that change the topology of the the medium when viewed globally. This framework has inspired a great deal of investigation in recent years. Such works encompass several distinct aspects to investigate the effects produced by topological defects of this nature [23][24][25][26][27][28][29].
In this work, we show rigorously how the dynamics of a 2D Dirac oscillator is modified when addressed in the presence of this topological defect, giving a special attention to the effects of spin. Our approach is based on the self-adjoint extension method which is appropriate to address any system endowed with a singular Hamiltonian (due to localized fields sources or quantum confinement). We determine the most relevant physical quantities from the present model, such as energy spectrum, wave functions and self-adjoint extension parameter by applying boundary conditions allowed by the system. * fmandrade@uepg.br † edilbertoo@gmail.com

II. THE 2D DIRAC OSCILLATOR IN THE MAGNETIC COSMIC STRING BACKGROUND
In this section, we study the motion of the particle in the cosmic string background. The cosmic string spacetime with an internal magnetic field is an object described by the following line element in cylindrical coordinates (t, r, ϕ, z): with −∞ < (t, z) < ∞, r ≥ 0 and 0 ≤ ϕ ≤ 2π. The parameter α is related to the linear mass densitym of the string by α = 1 − 4m runs in the interval (0, 1] and corresponds to a deficit angle γ = 2π(1 − α). Geometrically, the metric (1) corresponds to a Minkowiski space-time with a conical singularity [30]. We start with the Dirac equation in the curved spacetime (with = c = G = 1): where e is the charge, M is mass of the particle, ψ is a four-component spinorial wave function and Γ µ is the spin connection. The calculation of Γ µ can be found in some works in the literature. Here, we chose the same tetrad as in Ref. [30] (see also Ref. [28]). The recipe for determining Γ µ is purely algebraic and the result is found to be with the non-vanishing element given as where Σ z is the z-component of the spin vector Σ [29,30]. By exploiting the symmetry under z-translations, the (2+1)-dimensional Dirac equation is obtained from the decoupling of (3+1)-dimensional Dirac equation (2) for the specialized case where ∂ z = 0 and A z = 0, into two uncoupled two-component equations [27,31,32]. The relevant equation is the Dirac equation in (2 + 1) dimensions where is the generalized momentum, with the Nabla operator in cylindrical coordinates in the conical space and ψ is a two-component spinor. In (2 + 1) dimensions the γ matrices are given in terms of the Pauli matrices in cylindrical coordinates: where s is twice the spin value, with s = +1 for spin "up" and s = −1 for spin "down". The magnetic flux tube in the background space described by the metric (1) is related to the magnetic field as where φ = Φ/Φ 0 is the flux parameter with Φ 0 = 2π/e, and the vector potential in the Coulomb gauge is Note that in the limit as α → 1, we obtain the vector potential in Euclidean space. The 2D Dirac oscillator is introduced by the nonminimal substitution [1], where r is the position vector and ω the frequency of the oscillator (for a comprehensive discussion of the Dirac oscillator see Ref. [14]). In this case, Eq. (5) reads The second order equation implied by Eq. (12) is obtained by applying the matrix operator [βM + E + α · (π − iM ωβr)]: By using anticommunting relation {α, β} = 0 and the Inserting Eqs. (9), (10) and the expression for Γ ϕ in (4), one obtains where In Eq. (16), the quantity contributes to the term which depends explicitly of the spin of the particle. The first term is the contribution due to the magnetic flux while the second is due to the spin connection. Note that, by making α = 1 (flat spacetime) and φ = 0 (absence of magnetic field) in Eq. (15), we obtain, for the planar case, the 2D Dirac oscillator as proposed by Moshinsky and Szczepaniak [1] and discussed in Appendix A.
Making use of the underlying rotational symmetry we can express the two component spinor as with m ∈ Z. By replacing Eq. (18) into Eq. (15), we obtain the radial equation for f m (r) where and The Hamiltonian in Eq. (22) governs the dynamics of a Dirac oscillator in a magnetic cosmic string background, i.e., a Dirac oscillator problem in the presence of the Aharonov-Bohm effect in a conical spacetime. The presence of a two-dimension δ interaction in the radial Hamiltonian H, which is singular at the origin, makes the problem more complicated to solve. The most adequate manner to address this kind of point interaction potential is by making use of the self-adjoint extension approach [33,34]. This is the method adopted in this work and discussed in the next section.  [35][36][37]. The self-adoint extension approach consists, essentially, in extending the domain of D(O) in order to match D(O † ). From the theory of symmetric operators, it is a well-known fact that the symmetric radial operator H 0 is essentially self-adjoint for |j| ≥ 1, while for |j| < 1 it admits an one-parameter family of self-adjoint extensions [38], H 0,λm , where λ m is the self-adjoint extension parameter. To characterize this family, we will use the approach in [33,34], which is based in a boundary conditions at the origin. All the self-adjoint extensions H 0,λm of H 0 are parametrized by the boundary condition at the origin with where λ m ∈ R is the self-adjoint extension parameter. For λ m = 0, we have the free Hamiltonian (without the δ function) with regular wave functions at the origin, and for λ m = 0 the boundary condition in Eq. (24) permit an r −|j| singularity in the wave functions at the origin.

IV. THE BOUND STATE ENERGY AND WAVE FUNCTION
In this section, we determine the energy spectrum for the Dirac oscillator in the cosmic string background by solving Eq. (19). For r = 0, the equation for the component f m (r) can be transformed by the variable change ρ = M ωr 2 resulting in with γ = M ω. Due to the boundary condition in Eq.
With this, Eq. (27) is rewritten as Equation (27) is of the confluent hypergeometric equation type In this manner, the general solution for Eq. (27) is with In Eq. (31), M (a, b, z) is the confluent hypergometric function of the first kind [39] and a m and b m are, respectively, the coefficients of the regular and irregular solutions.
In this point, we apply the boundary condition in Eq. (24). Doing this, one finds the following relation between the coefficients a m and b m We note that lim r→0 + r 2−2|j| diverges if |j| ≥ 1. This condition implies that b m must be zero if |j| ≥ 1 and only the regular solution contributes to f m (r). For |j| < 1, when the operator H 0 is not self-adjoint, there arises a contribution of the irregular solution to f m (r) [29,[40][41][42][43][44][45].
In this manner, the contribution of the irregular solution for system wave function steams from the fact that the operator H 0 is not self-adjoint. For f m (r) be a bound state wave function, it must vanish at large values of r, i.e., it must be normalizable. So, from the asymptotic representation of the confluent hypergeometric function, the normalizability condition is translated in From Eq. (33), for |j| < 1 we have b m /a m = λ m γ |j| . Using this result into Eq. (34), one finds Equation ( with n a nonnegative integer, n = 0, 1, 2, . . .. By manipulation of Eq. (36), we obtain In particular, it should be noted that for the case when |j| ≥ 1 or when the δ interaction is absent, only the regular solution contributes for the bound state wave function (b m = 0), and the energy is given by Eq. (37) using the plus sign. Note that, for α = 1 (flat space) and φ = 0 (no magnetic flux), Eq. (37) coincides with the energy found for the usual 2D Dirac oscillator (cf. Eq. (A6) in Appendix A). Without loss of generality, let us suppose 0 < φ < 1 [37,46,47]. In this interval, and recalling that we are interested in the case where 0 < α ≤ 1, another interesting feature is present in the energy eigeinvalues.
For the regular solution, the eigenvalues are independent of m, φ and α for s = 1. This situation is shown in Fig.  1(a) for n = 1 and m = 1. However, this independence is absent for s = −1, as shown in Fig. 1 The self-adjoint extension is related with the presence of the δ interaction. In this manner, the self-adjoint extension parameter must be related with the δ interaction coupling constant φs/α. In fact, as shown in Refs. [28,29] (see also Refs. [40,48]), from the regularization of the δ interaction, it is possible to find such a relationship. Using the regularization method, one obtains the following equation for the bound state energy By comparing Eqs. (35) and (39), this relation is found to be where r 0 is a very small radius which comes from the δ regularization [28,29].

V. NONRELATIVISITIC LIMIT
We shall now take the nonrelativistic limit of Eq. (14).
Performing the same steps as for the relativistic case, one obtains the shifted energy levels (cf. Appendix A) In this equation, the + (−) sign is for λ m = 0 (λ m = ∞) when one has regular (irregular) solution. We note that the energy in Eq. (42) corresponds to the equation (54) of Ref. [49] (cf. also Eq. (A8) in Appendix A) with two additional contributions, the spin-orbit coupling and the spin connection.

VI. CONCLUSIONS
In this contribution, we have addressed the Dirac oscillator interacting with a topological defect and in the presence of the Aharonov-Bohm potential. This system has been studied in Ref. [17]. However, the authors not take into account the effects of spin. In other words, the term proportional to the δ interaction was discarded, by considering only the regular solution of the problem. The presence of this term has direct implications in the energy spectrum and wave functions of the oscillator. The correct approach to this problem must include spin effects, which are explicitly manifested by the spin-orbit coupling term, and so, we have a complete description for the dynamics of the 2D Dirac oscillator. We consider the self-adjoint extension method and show that the spinorbit coupling term, which results in a δ interaction, can not be dropped from the Hamiltonian. Although being singular at the origin, this term reveals that both regular and irregular solutions contributes for the bound state wave function and, consequently, for the energy spectrum. Expressions for the bound states energy for different values of the self-adjoint extension parameter were obtained. For two specific values for the self-adjoint extension parameter, i.e., λ m = 0 and λ m = ∞, the bound state energies are given explicitly in Eq. (37). We also verified that, for the flat space (α = 1) and no magnetic flux (φ = 0), the results of usual 2D Dirac oscillator are recovered.