On the 2D Dirac oscillator in the presence of vector and scalar potentials in the cosmic string spacetime in the context of spin and pseudospin symmetries

The Dirac equation with both scalar and vector couplings describing the dynamics of a two-dimensional Dirac oscillator in the cosmic string spacetime is considered. We derive the Dirac-Pauli equation and solve it in the limit of the spin and the pseudo-spin symmetries. We analyze the presence of cylindrical symmetric scalar potentials which allows us to provide analytic solutions for the resultant field equation. By using an appropriate ansatz, we find that the radial equation is a biconfluent Heun-like differential equation. The solution of this equation provides us with more than one expression for the energy eigenvalues of the oscillator. We investigate these energies and find that there is a quantum condition between them. We study this condition in detail and find that it requires the fixation of one of the physical parameters involved in the problem. Expressions for the energy of the oscillator are obtained for some values of the quantum number n. Some particular cases which lead to known physical systems are also addressed.


Introduction
The study of the relativistic quantum dynamics of particles including electromagnetic interactions is an usual framework for studying properties of various physical systems. The mechanism used to describe these systems is a natural generalization of the coupling used in classical nonrelativistic quantum theory [1]. This coupling is implemented, for charged particles with charge e, through the so-called a e-mail: daniel.franca.lima@gmail.com b e-mail: fmandrade@uepg.br c e-mail: lrb.castro@ufma.br d e-mail: cleversonfilgueiras@yahoo.com.br e e-mail: edilbertoo@gmail.com minimal coupling prescription, given in terms of the modification of the 4-momentum operator, p μ → p μ − e A μ = ( p 0 − e A 0 , p − eA), where A μ = (V (r ), −A) (with A being the vector potential and V (r ) being the scalar potential) represents the 4-vector potential of the associated electromagnetic field. This transformation preserves the gauge invariance associated with the Maxwell's equations. Another way to insert interaction in the dynamics of the particle is by including a scalar potential through a modification in the mass term as M → M + S(r ). In this realization, the potential S (r ) is coupled like a scalar, different from the minimum prescription, where the potential is coupled as a time-like component of a 4-vector. Although there is some similarity between the scalar and vector couplings, they have different physical implications. Actually, the scalar coupling acts equally on particles and antiparticles. On the other hand, the vector coupling acts differently on particles and antiparticles. As a result, the energy of particle and antiparticle are not equals, so that bound states exist only for one of the two kinds of particles [2].
Interesting issues that should be investigated with the insertion of the couplings in the Dirac equation are the socalled the spin and the pseudo-spin symmetries [3]. Basically, these symmetries occur when the couplings are composed by a vector V (r ) and a scalar S(r ) potential, under the assumption that S(r ) = V (r ) (S(r ) = −V (r )), which is the necessary condition for occurrence of exact spin (pseudo-spin) symmetry. The spin symmetry has been identified by studying heavy-light mesons [4], single antinucleon spectra [5] and dynamics of a light quark (antiquark) in the field of a heavy antiquark (quark) [3] while that the pseudo-spin symmetry occurs in the motion of nucleons [3,6]. In recent studies, both the spin and the pseudo-spin symmetries appear in several aspects concerning, for instance, the supersymmetry [7,8], the Hartree-Fock theory [9], the electrons in graphene [10] and the interaction with a class of scalar and vector potentials [11][12][13][14][15][16][17][18]. An important physical system that can be studied by including such terms of interactions in the Dirac equation is the Dirac oscillator [19] (for a detailed description of this model see Ref. [20]). The Dirac oscillator is a kind of tensor coupling with a linear potential which in the nonrelativistic limit leads to the simple harmonic oscillator with a strong spin-orbit coupling. It was realized experimentally for the first time in 2013 by Franco-Villafañe et al. in [21]. The Dirac oscillator is considered a natural model for studying properties of physical systems because it is exactly soluble. In the last years, several research have been developed in the context of this theoretical framework. For instance, it appears in the literature in the context of mathematical physics [22][23][24][25][26][27][28][29], nuclear physics [30][31][32][33], quantum optics [34][35][36][37], supersymmetry [38][39][40], theory of quantum deformations [41,42] and noncommutativity [43][44][45][46]. Moreover, the Dirac oscillator embedded in a cosmic string background has inspired a great deal of research in last years [47][48][49][50][51][52][53][54][55][56].
In this work, we analyze in details the solutions of the Dirac equation with both scalar and vector interactions under the spin and the pseudo-spin symmetry limits in the cosmic string spacetime [57]. Cosmic strings are topologically stable gravitational defects. According to the grand unified theories, these defects arise from a vacuum phase transition in the near universe. Recently, several studies have been developed in the theoretical context [58][59][60][61][62][63] and also by evidence of cosmic strings [64][65][66][67]. Cosmic strings are objects of studies of current interest because of the several important applications of topological features on physics systems in gravitation [68], condensed matter [57] and cosmology [69].
Our work is motivated by Ref. [70] (see also Refs. [71,72]), where the spin and pseudospin symmetries in the relativistic mean field with a deformed potential are investigated. In this context, a relation between the deformed wave function and the spherical wave function was established at the spherical limit by using the transformation from the cylindrical coordinate into the polar coordinate. This relationship enables us to investigate the inclusion of cylindrical symmetrical potentials in the Dirac equation in other scenarios, such as the cosmic string. One advantage of using such symmetry limits in our work is that they allow us to decouple the first and second order differential equations for the spinor components (each obtained in the spin symmetry and pseudo-spin limits, respectively).
We organize the paper as follows: In Sect. 2, we derive the equation that governs the dynamics of a Dirac particle with the minimal, nonminimal and the scalar couplings in the cosmic string spacetime. In Sect. 3, we consider the Dirac equation written in terms of a set of coupled differential equations. We investigate the existence of particular solutions for the problem by assuming that the relativistic energy of the particle is its rest energy in both the spin and the pseudo-spin symmetries limits. In Sect. 4, we investigate the dynamics considering that the energy of the particle is different from its rest energy. To this end, we write down the Dirac equation in its quadratic form. We obtain the energies and the corresponded wave functions and discuss their physical validity. In Sect. 5, we address some particular solutions and compare them with previous results in the literature. Finally, the conclusions are presented in Sect. 6. Here, we use natural units such ash = c = 1.

The equation of motion
In this section, we derive the Dirac equation with scalar and vector couplings to study the motion of a Dirac oscillator in the cosmic string spacetime. We first define the spacetime background of an idealized cosmic string where the oscillator will move, followed by the most general interaction, which includes the potential of the Dirac oscillator. The interactions, however, are chosen in such a way that analytical solutions to the Dirac equation can be obtained.
The spacetime generated by a cosmic string is described by the following line element in cylindrical coordinates with −∞ < (t, z) < ∞, r ≥ 0 and 0 ≤ ϕ ≤ 2π . The parameter α is related to the linear mass densitym of the string by α = 1 − 4m and it runs in the interval (0, 1] and corresponds to a deficit angle γ = 2π(1−α). Geometrically, the metric in Eq. (1) corresponds to a Minkowski spacetime with a conical singularity [73]. One starts by considering the free Dirac equation, i.e., in the absence of interactions. The interaction will be included later. So, we have where Ψ is a four-component spinorial wave function. In order to work out in the curved spacetime, we must write the Dirac gamma matrices γ μ in the Minkowskian spacetime (written in terms of local coordinates) in terms of global coordinates and subsequently include the spinor affine connection Γ μ . In other words, we must contract γ μ with the inverse tetrad, satisfying the generalized Clifford algebra where where σ i are the standard Pauli matrices and 1 is the 2 × 2 identity matrix. As we are interested on in a cosmic string, we need to write down the generalized Dirac equation in the curved spacetime background with a minimal coupling. Therefore, the relevant equation is where e is the electric charge and A μ denotes the vector potential associated with the electromagnetic field. The spinor affine connection is often written as [74] where ω μab is the spin connection, given by In (8), Γ ν τ μ are the Christoffel symbols and η ab is the metric tensor. By the means of the spin connection, we can construct a local frame using a basis tetrad which gives the spinors in the curved spacetime. Here, the basis tetrad e μ a is chosen to be [75] satisfying the condition e μ a e ν b η ab = g μν .
Using (9), the matrices γ μ in Eq. (6) are written more explicitly as Given the fact that the matrices in the curved space satisfy the condition ∇ μ γ μ = 0, i.e., they are covariantly constant, for the specific basis tetrad (9), the affine spin connection is found to be with the non-vanishing element given by We are interested on including potentials with cylindrical symmetry, in such a way the resulting system will have translational invariance along the z direction. Then, we can discard the third direction and thus consider the Dirac oscillator in two spacial dimensions [19] (see also Ref. [20]), assuming p z = 0. 1 This assumption allows us to reduce the fourcomponent Dirac equation (6) to a two-component spinor equation. Moreover, according to the tetrad postulated [74], the γ a matrices could be any set of constant Dirac matrices. Thus, a convenient representation is the following [54,76,77] where the parameter s, which is twice the spin value, can be introduced to characterize the two spin states, with s = +1 for spin "up" and s = −1 for spin "down". In the representation (19), the matrices (11), (13) and (15) assume the following form: and Eq. (18) becomes Now, let us include the interactions into the Dirac equation (6). We consider the effective potential [78,79] Mωiσ z βγ ·r r + with where are cylindrically symmetric scalar and vector potentials. The first term in Eq. (24) represents the Dirac oscillator. In this manner, the time-independent Dirac equation (6) with energy E can be written as where ψ is a two-component spinor, is the Dirac Hamiltonian and is the planar spatial part of the gradient operator in the metric (1). We begin the study of the particle motion by looking for first order solutions of the Eq. (29). For this purpose, we write the Eq. (29) as follows, and and we consider the solutions as with m = 0, ±1, ±2, ±3, . . . being the quantum angular momentum number. The substitution of (33) into (32a) and (32b) gives the following set of coupled differential equations: where Θ + = 0 and Θ − = 1. The reason why we are using superscripts (±) in Eq. (36) will be clarified in the next section. If we consider that Δ (r ) = 0 and E = −M or Σ (r ) = 0 and E = +M, the solutions of Eqs. (34) and (35) represent a particular solution for the problem, which is excluded from the Sturm-Liouville problem. In other words, such solutions would not be part of those obtained by solving the second-order differential equation obtained from Eq. (29). The procedure of imposing that either Δ (r ) = 0 or Σ (r ) = 0 in Eqs. (34) and (35), respectively, is known in the literature as the exact limits of spin and pseudo-spin symmetries [3]. These conditions are taken into account in the next section.

Particular solutions and the analysis of the spin and the pseudo-spin symmetries
In this section, we solve the system of first-order radial differential equations obtained in the previous section by imposing either the exact limits of spin and pseudo-spin symmetries.
Once we find the solutions, we must verify that they are physically acceptable solutions. As mentioned above, the exact limit of the spin symmetry occurs when Δ (r ) = 0 (V (r ) = S(r ) in Eq. (25)), while that the exact limit of the pseudospin symmetry is achieved by setting Σ(r ) = 0 (V (r ) = −S(r ) in Eq. (26)). In what follows, the superscript (+) holds for the spin symmetry and (−) holds for the pseudo-spin symmetry. In these limits, the solutions are related to the up and down components of the spinor in Eq. (33), respectively. In order to obtain the particular solutions, let us look for the bound state solutions which obey the following normalization condition, We assume E = ±M, as it was mentioned above.

The exact spin symmetry
Here, the particular solutions for the bound states are obtained by considering Δ(r ) = 0 2 along with the assumption E = −M in both Eqs. (34) and (35). Therefore, we have Their solutions are written as with where are upper incomplete Gamma functions [80], a 1 and a 2 are constants. Let us discuss the solutions (40) and (41). Since e − 1 2 Mωr 2 dominates over r s J + α for any value of s J + α , the solution f m (r ) in Eq. (40) converges as r → 0 and r → ∞. On the other hand, as the incomplete Gamma functions As Mω > 0 in (46), there are no values of s J − α for which the functions are square-integrable. In this case, we can therefore conclude right away that for E = −M and exact spin symmetry there is no bound state solution.

Exact pseudo-spin symmetry
In this case, we impose Σ(r ) = 0 and E = M in both Eqs. (34) and (35). Thus, we obtain Their solutions are given by where b 1 and b 2 are constants, and with Again, the incomplete Gamma functions Γ (d,e, f ) in Eq. (49) always diverge, so that a normalized solution requires that b 2 = 0. In such a case, the function f m (r ) is squareintegrable only for s J + α ≥ 0. The physically acceptable solution is Therefore, we can conclude that for the case E = M along with the exact pseudo-spin symmetry there is a bound state solution. Here, the existence of a particular bound state solution is guaranteed only for Mω > 0. However, there are other models in the literature where this quantity can assume any value, so that bound states solutions are allowed for both the spin and pseudospin symmetry limits [81].

The Dirac-Pauli equation and the analysis of both the spin and the pseudo-spin symmetries
In this section, we study the dynamics for the case E = ±M. For this purpose, it is more convenient to work with the Eq. (29) in its quadratic form. In our analysis, we shall see that because of the shape of the potential (24), the solutions for the radial equation are given in terms of biconfluent Heun functions and the energy levels of the oscillator will be determined only after imposing some quantum conditions. To obtain the quadratic form of the Dirac equation (29), we multiply it by the matrix operator leading to where ∇ 2 α = ∂ 2 r + (1/r )∂ r + (1/α 2 r 2 )∂ 2 ϕ is the planar spatial part of the Laplace-Beltrami operator in the metric (1). By inserting the solutions (33) into Eq. (57), we obtain the following set of two coupled radial differential equations of second-order: and Notice that these two equations are coupled via the last terms and the spin and pseudospin symmetry limits uncouple them. So, here and henceforth we employ the following approach. For the spin symmetry limit, we solve the problem by considering the upper component of the spinor and denotes it by f + (i.e., + labels the spin symmetry solution) and for the pseudospin symmetry limit, we consider the lower component and denotes it by g − (i.e., − labels the pseudospin symmetry solution).

The analysis of both the spin and the pseudo-spin symmetries
When we take into account the exact limits of spin and symmetries in Eqs. (58) and (59), each component of the spinor satisfies and where = Mω, a ± = 2(E ± M)η C and b ± = 2(E ± M)η L . The differential equations (61) and (61) can be placed in an convenient mode using, respectively, the following solutions: where x = √ r and y ± (x) satisfies where (65) is a homogeneous, linear, second-order, differential equations defined in the complex plane. The solutions of these equations are given in terms of the biconfluent Heun functions by [82,83] where The coefficients of the series are given by and From the recursion relation (73), the function becomes a polynomial of degree n, if and only if, the two following conditions are imposed [83,84]: In this case, the (n + 1)th coefficient in the series expansion is a polynomial of degree n in 2 ξ ± C . When 2 ξ ± C is a root of this polynomial, the (n + 1)th and subsequent coefficients cancel and the series truncates, resulting in a polynomial form of degree n for the solution N ± 2 |J ± α |, ξ ± L , Δ ± , 2 ξ ± C , x . From the condition (75), we extract the following expressions involving the energy E ± nm : We notice in Eq. (77) the absence of the parameter η C . This steams from the fact that these expressions do not represent the energies of the system in its present form. Actually, the condition (76) allows us to establish a quantum condition that links the energy and others physical quantities, including η C [79,85,86]. As a result, it is possible to express the energy in terms of all the physical parameters involved in the problem, namely, η C , η L , M, and ω. We emphasize that that, a priori, we are free to choose which parameter we want to fix. Here, such a quantum condition is established through the frequency ω of the system. Therefore, we now label ω as ω nm . Before performing the procedure, let us consider the solution (70) up to second-order in x of the expansion, namely, with Thus, Eq. (78) reads Now let us determine the quantum condition mentioned above. For the condition (76), we must investigate A ± n+1 = 0. For simplicity, we consider only the case n = 0, which requires that A ± 1 = 0 in Eq. (80). This requires us to solve the equation which provides the following frequencies related to the ground state of the system: However, Eq. (83) will only be an acceptable quantum condition if η L /η C < 0 to ensure that the frequencies ω ± 0m are positive. Thus, respective energies corresponding to the ground state are where In (84), the notation (±) refers to the particle and antiparticle energies. The energies in Eq. (84) now depend on all the physical parameters involved in the problem. In Figs. 1 and 2, we plot the profile of these energies as a function of the parameter α. In both plots we clearly see that the energy levels of the particle and antiparticle belong to the same spectrum and, moreover, there is no channel that allows the spontaneous creation of particles because none of the lines of the spectrum cross each other.

Particular cases
In this section, we study particular solutions of problem solved in the previous section. Namely, we will investigate (a) (b) Fig. 2 Illustration of the energy eigenvalues in the pseudospin symmetry limit, E − 0m , as a function of the parameter α. In a the plot for s = 1 and b for s = −1. We use M = 1, η C = 1 and η L = 1. The energies of the states corresponding to a given value of m near α = 0 in a and b are more energetic while near α = 1 the differences between the energy levels decrease as well as their respective values three cases. For the first two, the solutions of the resulting equations are given in terms of biconfluent Heun functions whereas the third, which will not involve scalar and vectorial interactions, will be given in terms of the confluent hypergeometric function.
Let us then return to Eq. (65) and solve it for the particular case η L = 0. The resulting equation governs the dynamics of a two-dimensional Dirac oscillator interacting with the potential η C /r . In this case, the solutions are given bỹ Then, using the condition (75), we find the energies Moreover, from condition (76), we consider again A ± n+1 = 0 for n = 0, and solve it for˜ ± 0m . One can thus verify that it is not possible to extract a physically acceptable expression for˜ ± 0m . Consequently, n = 0 is not an allowed value for the quantum number and we need to solve A ± n+1 = 0 for n = 1. Thus, we havẽ Substituting (89) into (88) and solving these equations for E ± 1m , we find and where the subscripts p and ap refer to the energies of the particle and antiparticle, respectively. As we are studying the dynamics for which E ± 0m = ±M, the energies E + 1m ap and E − 1m p are not allowed energies for the particle. The profiles of the energies (90a) and (91b) as a function of the parameter α are shown in Figs. 3 and 4, respectively. We can observe in Fig. 3a (s = +1) the presence of degeneracy for m = −2, −1, 0, while in Fig. 3b (s = −1), the degeneracy occurs for m = 0, 1, 2. In Fig. 4, the spectrum of the states with m = −2 (Fig. 4a for s = +1) and with m = 2 (Fig. 4b for s = −1) change very slowly and are non-degenerate.
The second particular case is when η C = 0. In this case, the system consists of a Dirac oscillator interacting with a linear potential, η L r . Thus, the solutions of Eq. (65) is again given in terms of the Heun functions, Note that energies (94) are identical to those given in Eq. (77). However, the frequencyω is not the same. The difference between them is just the imposition established by the condition (76). For n = 0, we obtain the frequencies By substituting (95) into the respective energies (94), we find and the energies are given by withω ± 1m given in Eq. (97). For this particular case, it is verified that Eq. (98) presents four energy eigenvalues being two for each type of symmetry limit considered. However, only We can see that both particle and antiparticle belong to the same spectrum and contains no degeneracy. In Fig. 5a, we clearly observe that the states with m > 0 are more affected by the curvature while in Fig. 5b this occurs for the states with m > 0. These same characteristics are also present in Fig. 6, the only difference is that the spacing between each level as well as their respective energy values are larger when compared with the spectra of the Fig. 5.
Finally, the last case we want to discuss in that in which η L = η C = 0 in Eq. (65). In this case, the solutions (69) and (69) take the form where x = √ r and F ± (x) satisfies the Kummer differential equation [80,82] whose general solution is known to be In the above equations, M is the Kummer function [80,82].
For this particular case, if we write the condition (75) in the form and the unnormalized bound state wave functions are The energies in Eq. (103) (for n = 0 and s = ±1) are plotted as a function of the parameter α in Figs. 7e 8, respectively. For a particle with s = 1 (Fig. 7a), all states with m > 0 are degenerate and are not affected by curvature while with s = −1 (Fig. 7b), this characteristic occurs for the states with m < 0. On the other hand, for an antiparticle with s = 1 (Fig. 8a), only the state with m = −2 is non-degenerate while with s = −1 (Fig. 8b), only the state with m = 2 is non-degenerate. In Ref. [87], the Dirac 2D oscillator interacting with the Aharonov-Bohm potential in the time space of the cosmic string was studied in the context of self-adjoint extensions. In the absence of the Aharonov-Bohm field, the resulting equation corresponding to the regular solution (Eq. 46) of Ref. [87] reproduces the Eq. (103).

Conclusion
In this paper, we have studied the dynamics of a 2D Dirac oscillator interacting with cylindrically symmetric scalar and vector potentials in the space-time of the cosmic string. The problem was solved taking into account the spin and pseudospin symmetry exact limits through two stages. First we have solved the Dirac equation by looking for first order solutions. We used an appropriate ansatz for the Dirac equation and obtained a system of coupled first order differential equations. We investigated this system and verified that it admits physically acceptable particular solutions, i.e., bound states solutions, only for the pseudo-spin symmetry exact limit, Σ = 0 and E = M. In the second moment, we have constructed and solved the Dirac equation in its quadratic form, which excludes the E = ±M cases from its solutions. For this case, we shown that the resulting radial differential equation is the biconfluent Heun equation. We studied the series solution of this equation as well as its asymptotic behavior at infinity and at the origin and found two conditions (Eqs. (75) and (76)) to make the series a polynomial. The use of these two conditions allowed us to obtain expressions for the energies corresponding to fixed values of n. In particular, we obtained the expression corresponding to the state with n = 0, which is given by Eq. (84). We investigate how the curvature affects the energies. For this intent, we have plotted it as a function of the parameter α for each of the limits of symmetries and spin element projection considered. In the case of the energy obtained for the spin symmetry limit (Eq. (84) with superscript +), we have shown that for s = 1 the states with m < 1 become more energetic when α → 0 while for α = 1 the differences between the energy levels as well as the respective energy values decrease. For the states with m 1, the energies change very slowly and are non-degenerate. When the spin element is s = −1, we have verified that the effects are opposite to those for s = 1, namely, the states with m 0 are more energetic for α → 0 and less energetic for α = 1. These characteristics were also observed in the graph of the energies obtained in the pseudspin symmetry limit (Eq. (84) with superscript −). For both s = ±1, the energies of the states corresponding to a given value of m when α → 0 are more energetic while for α = 1 the differences between the energy levels decrease as well as their respective energy values.
We also investigated some special cases for the solution of the Eq. (57). In the first case, we have assumed the vanishing of the linear potential by imposing η L = 0. We obtained the energies (Eqs. (90a)-(90b) and (91a)-(91b)) and plot them as a function of the parameter α for both s = ±1. However, we have shown that the energies (90b) and (91a) are not allowed. In the energy profile (90a) for s = 1, the energies of states with m 0 are degenerate. In particular, when α = 0.5, the state energy with m = 1 is not defined. We also have observed that the energy values of the state with m = 2 when α → 0. increases while for α = 1 it decreases. It also was verified that these same characteristics are present in the graphic for s = −1. In the plot of the energy given by Eq. (90b) for s = ±1, other important characteristics were manifested, and these are absent in the plot of Eq. (90a). For s = 1, the energy of the states are not defined for α equal to 0.25, 0.42 and 0.59. The spectrum is more energetic for α → 0 and α = 1, except the m = −1 curve, in which is more energetic only for α → 0. We have found that energy of the state with m = −2 changes very slowly and are non-degenerate. We have also found that these characteristics are present in the graphic for s = −1.
In the second particular case investigated, we have assumed η C = 0 and, as for the first case, four energy eigenvalues were found, but only two of them are physically acceptable because of the requirement that E = ±M. For this case, we have not found energies with a given values of m and α that are not allowed. The graphs of the energies (for s = ±1) as a function of the α for both spin and pseudo-spin symmetry limits revealed that they are more energetic for α → 0 and less energetic for α = 1.0. The only difference is that the spacing between the energies of the states for a fixed m in the spin symmetry limit are greater than those in the spin symmetry limit.
In the last particular case studied, we have assumed η L = η C = 0. For this system, the resulting radial equation was a equation Kummer differential equation type. We obtained the energy spectrum ( ± nm in Eq. (103)) and we plotted it as a function of the α for both s = ±1. In the graph of the energy + 0m for s = 1, we have verified that the states with m > 0 are degenerate while for s = −1 this occurs for states with m < 0. In the graph of the energy − 0m , we have found that only the states with m = −2 (for s = 1) and with m = 2 (for s = −1) are non-degenerate. A feature present in all energy profiles, including the general case, is the absence of channel that allows creation of particles, and also no crossings of lines, which guarantees that particle and antiparticle belong to the same spectrum.
As a final remark, we would like to mention that the model addressed here can be applied to other systems, especially those in condensed matter physics. This is due to the fact that linear defects in condensed matter, such as disclinations and dislocations in solids, can be studied through the same approach used to treat a cosmic string [57]. A possible application would be an adaptation of the model used to investigate how the quantum dots and antidots, with the pseudoharmonic interaction and under the influence of external magnetic and Aharonov-Bohm potential are influenced by the presence of a screw dislocation as that studied in Ref. [88] in the context of spin and pseudo-spin symmetries. Interesting investigations can also be made by considering non-inertial effects on the particle dynamics [89]. The inclusion of non-inertial effects in relativistic and non-relativistic quantum mechanics is an issue of current interest it may be interesting to study some physical system in the scenario of the problem addressed here or in some other particular geometry.