Abstract

The generalized Dirac oscillator as one of the exact solvable models in quantum mechanics was introduced in 2+1-dimensional world in this paper. What is more, the general expressions of the exact solutions for these models with the inverse cubic, quartic, quintic, and sixth power potentials in radial Dirac equation were further given by means of the Bethe ansatz method. And finally, the corresponding exact solutions in this paper were further discussed.

1. Introduction

As we all know, Dirac oscillator is a concept which was first proposed by M. Moshinsky [1] and collaborators by means of the combination between momentum and coordinate , where is the coordinates, represents the mass of particle, and is regarded as the frequency of particle [2, 3]. Although the Dirac oscillator, an effective relativistic model, has succeeded in describing the interaction between anomalous magnetic moment and linear coordinates in many fields, it is powerless to explain the interaction between anomalous magnetic moment and nonlinear coordinates. Therefore, in order to describe the complex interactions, Dutta and his colleagues [4] proposed a new concept which is called the generalized Dirac oscillator. So it could generalize the linear effect between coordinates and momentum to nonlinear effect by using this system; that is to say, the Dirac oscillator could be regarded as a special case in the generalized Dirac oscillator.

As one of the few relativistic quantum systems that can be solved accurately, Dirac oscillator has attracted a lot of attentions [15]. This system not only can solve the quark confining potential [57] in quantum chromodynamics, but also can be used to solve some complex interactions [8, 9]. Considering the practicability of solving problems, Dirac oscillator has been extensively discussed by many researchers in various aspects: like conformal invariance properties [8], covariance properties [9], shift operators [10], symmetry Lie algebra [11], hidden supersymmetry [1214], completeness of wave functions, and so on [15, 16]. Moreover, many phenomena in Quantum Optics [1721] could also be explained by this model. The investigation of the relationship between Dirac oscillator and relativistic Jaynes-Cummings model [4, 20, 21] bridges the two unrelated fields of Relativistic Quantum Mechanics and Quantum Optics [22, 23]. In addition, this model can also be used to explain some new phenomena in condensed matter physics [16, 17, 24, 25], such as the Quantum Hall Effect and Fractional Statistics. So it is interesting to extend Dirac oscillator to generalized Dirac oscillator for solving the interaction between momentum and coordinate in nonlinear electric field.

To solve 1+1-dimensional Dirac equation with complex interaction potential, Dutta and colleagues constructed the generalized Dirac oscillator by making a simple transformation . But now, we want to solve the 2+1-dimensional radial Dirac equation with complex interaction potential, so the momentum operator is transformed and in this paper. Therefore, the corresponding complete solutions of 2+1-dimensional radial Dirac equation can be given by the Bethe ansatz method [2629] provided that there is appropriate selection of the interactions.

Next, the structure of this paper is as follows. The 2+1-dimensional generalized Dirac oscillator is introduced in Section 2. In Section 3, exact solutions of the radial Dirac equation with generalized Dirac oscillator are given by using the functional of Bethe ansatz method [2629] when selecting appropriate interactions. Section 4 is devoted to the conclusion. And finally, Bethe ansatz method will be briefly reviewed in the Appendix.

2. 2+1-Dimensional Generalized Dirac Oscillator

For particle of mass , the Dirac equation can be written in terms of two component spinors [3032]where and are matrices that are constituted by the Pauli spin matricesand the parameter takes the values ( for spin up and for spin down). The formalism (1) can transform two uncoupled two-component Dirac equations for or by using the decoupling for the four-component Dirac equation in the absence of the third spatial coordinate. Thus, one of the two-component Dirac equations could be used to describe the positive energy eigenstate and another to describe the negative energy eigenstate. And the positive energy eigenstate is what is usually called particle state, and the corresponding negative energy eigenstate is called antiparticle state. Due to the fact that there is only one spin polarization in 2+1-dimensional space, so the equation of Dirac oscillator for the particle state can be written as follows [32, 33]:where the matrices and the eigenfunction satisfy the following conditions:

Next, the momentum operators will be replaced by and in free radial Dirac equation, so that corresponding radial Dirac equation of generalized oscillator is given as follows:According to equation (4), equation (5) can be written as follows:So, the second order differential equation could be given by substituting the lower formula in equation (6) into the upper formulawhere . Next, the exact solution of this system will be further discussed in the following section by using the Bethe ansatz method.

3. Explicit Implementation of Generalized Dirac Oscillator

Over the past decades, the singular potential has attracted lots of attentions because it can be used to describe many physical problems. For example, not only is the singular potential widely applied in the research of the and procedures in high-energy physics [29, 34, 35], but also repulsive singular potentials could reproduce the interactions of scattering and nucleons with K-mesons [29, 36]. Moreover, high energy scattering question caused by strong singular potential, as an example of nonrelativistic quantum mechanics, has also been extensively discussed by many authors [29, 37, 38]. Besides, the singular potentials were also widely used to describe the interaction of two atoms in molecular physics, interatomic or intermolecular force and chemical physics [29]. Considering the wide applications of the singular potential, the exact solutions of 2+1-dimensional radial Dirac equation with generalized Dirac oscillator under the inverse cubic, quartic, quintic, and sixth power potentials were discussed in this paper. And it is also shown that the conclusion of the generalized Dirac oscillator with higher order inverse potential can degenerate to the conclusion of lower order inverse power potential interaction when the parameters are properly selected. Next, the analytical expressions of the corresponding exact solutions will be given in Section 3; the ground and first excited states were discussed in more detail.

3.1. Quasi Exact Solution of the Inverse Cubic Power Potential

Now, let us firstly discuss the inverse cubic power potentialThis potential has been investigated in Schrodinger equation by the other authors to obtain the accurate analytical expression [39, 40]. Now, in order to solve this problem in the 2+1-dimensional radial Dirac equation, the corresponding radial Dirac equation could be given by substituting equation (8) into equation (7):

The exact solution of the 2+1-dimensional radial Dirac equation is further given by extracting the asymptotic behavior of the wave function through a simple replacement. Next, by making a brief examination for differential equation, we implement transformationwhere the parameters , , , and are constant and is magnetic quantum number. Substituting equation (10) into differential equation (9), it is easy to find that parameters satisfy underlying relations

Therefore, new differential equation could be written as follows:where parameters satisfy . Obviously, equation (12) could also be solved accurately by using the Bethe ansatz method [2629] if the potential parameters satisfy certain constraints. Now, in order to realize the square integrable of wave functions, we assume that or ; then the degree polynomials solutions of equation (12) could be given:

Case 1 (). In this section, equation (12) would be transformed into the following form:Substituting equation (13) into equation (14) and using Bethe ansatz method, it is easy to prove that parameters satisfy the following constraints:And Bethe ansatz equation can be given as follows:Therefore, the analytical expressions of the wave functions and energy spectrum could be written asThe above wave functions are square integrable, so the corresponding normalization constants could also be given through the standard integral [41]where , , , and is the Meijer G-function. Thus, the total wave functions could be written asand here the matrix term of the total wave function was set as Next, the corresponding solution of ground state and first excited state could be given directly by using the expression of the exact solution. Obviously, the polynomial solution of equation (14) satisfies if we take . Then, the analytical solution of the ground state could be given as follows:with , . The case of corresponds to the first excited state solution of this system, and the analytical expressions of the wave function and the energy spectrum are further given:in which the potential parameters meet the needs of the following constraints:

Case 2 (). In this section, we assume that . Obviously, if the potential parameters satisfy certain constraints, the exact solution of the system could be given by using the Bethe ansatz method [2629]. And the condition of the wave function has an acceptable asymptotic behavior; that is, the parameter B must be negative when . So this differential equation readswhere . Substituting equation (13) into equation (24), the following relations could be given by using the Bethe ansatz method:and the Bethe ansatz equation could be written asAs mentioned above, the parameter is negative, so the parameter must also be negative. Therefore, the analytical expressions of the wave functions and energy spectrum could be written asThe wave functions are square integrable, and their corresponding normalization constants could be determined by standard integral [41]for , , . So, the total wave functions were written asNext, the ground and the first excited state will be investigated in detail. Obviously, the polynomial solution of equation (24) could be written as for the case ; then, corresponding energy spectrum and wave function arewhere the constraint satisfies . The case of corresponds to first excited state solution of this system; the analytical expressions of wave function and energy spectrum are further given aswith underlying constraint for potential parametersHere the root can be further given by using the Bethe ansatz equation

3.2. Quasi Exact Solution of the Inverse Quartic Power Potential

The interaction of the Dirac equation with inverse quartic power potential is studied in this part:From the phenomenological point of view, singular potential as a very useful form of anharmonicity is used in many aspects of physics [42, 43]. The inverse quartic potential has also been investigated in many different questions by lots of authors [4446]. Now, the purpose of our study is to obtain the analytical properties of the scattering amplitude about singular potential. So, the corresponding 2+1-dimensional radial Dirac equation with inverse quartic power potential is given by equation (7):In order to deal with the question, the appropriate asymptotic behavior of the wave function is extracted through a simple replacementwhere the parameter satisfies . So, this differential equation for readswhere the parameter is set as . Here, equation (37) can be transformed into a solvable form by making a new variable :So, the degree polynomial solutions for the differential equation readwith different roots if potential parameters fulfill certain constraints.

Thus the analytical expressions for energy states and wave functions could be given asAnd the parameters satisfy the restrictive conditionsand the roots could be determined by the Bethe ansatz equations:

The wave functions are square integrable, and their corresponding normalization constants could be determined by standard integral [41]for , , .

As special cases of general expressions, we focus on studying ground state and first excited state systems in this section. First, let us give the corresponding wave function and ground state energy spectrum in terms of equation (40):where the potential parameters comply with .

The case of corresponds to first excited state solution of the system, and the analytical expressions of the wave functions and energy spectrum for the first excited state are further written asin which the different roots could be computed analytically by the Bethe ansatz equation [2629]and the potential parameters satisfy the following restrictive condition:

3.3. Quasi Exact Solution of the Inverse Quintic Power Potential

Now, the inverse quintic power potential will be discussed in here:The exact solvable form of the radial Dirac equation with inverse quintic power potential could be further given by substituting equation (48) into equation (7):where the parameter satisfies . Similar to the studies in preceding sections, the corresponding exact solutions can be given provided that the appropriate asymptotic behavior of the wave function is extracted. Here we setwhere the parameters are constant and parameter is the magnetic quantum number. Substituting equation (50) into equation (49), we find that the parameters satisfy

Thus, the radial differential equation (49) will degenerate into the following form:where . It is easy to find that equation (52) is exactly solvable when making or , and its corresponding exact solutions could also be given by the degree polynomials

Case 1 (). Considering this case of , the radial differential equation (52) can be transformed intoSubstituting equation (53) into equation (54) and using the Bethe ansatz method, the requirements of the parameters are as follows:Here the roots fulfill Bethe ansatz equationwhere . And its corresponding exact solutions could be written as follows:The wave functions are also squarely integrable: . Next, the ground and first excited state will be discussed in detail as two special cases. And the case of corresponds to the ground state solution of the systemThe first excited state solution can be given by making :in which the root could be computed analytically by Bethe ansatz equationand its parameters satisfy the following relations:

Case 2 (). Here, we set . In order to get acceptable asymptotic behavior for the wave function when , the parameter is required to be negative. Furthermore, according to the Bethe ansatz method, we require that . Thus, the differential equation (52) could be written asSubstituting equation (53) into equation (62) and using the Bethe ansatz method, the exact solutions of equation (62) can be given as follows:and the parameters satisfy the determined constraintsHere the roots satisfy Bethe ansatz equationwhere . As we mentioned above, the parameter is negative, so the parameter must be negative. Next, in order to satisfy the wave functions square integrable we assume that . It is not difficult to use Maple for checking numerically.
Similar to the previous cases, the ground state solution corresponds to the case of :with the constraintsAnd the first excited state solution corresponds to the case of :complying with the constraintsHere the root satisfies

3.4. Quasi Exact Solution of the Inverse Sixth Power Potential

The effect of the inverse sixth power potential in radial Dirac equation is discussed as the last caseThe inverse power potential has played an important role in the study for particles interactions in atomic, molecular, and nuclear physics [4749]. Substituting equation (71) into equation (7), the corresponding radial Dirac equation with the inverse sixth power potential could be given asSimilar to the discussions in previous sections, in order to extract appropriate asymptotic behavior from the wave function , we setwhere the parameter satisfies . Then, according to equation (73), the differential equation for is written as follows:Replace the variable with and multiply on both sides of the equation, so that the equation can degenerate into the following form:It is easy to find that the degree polynomial solutions of the differential equation (75) satisfySubstituting equation (76) into equation (75) and applying the Bethe ansatz method, the analytical expressions of the energy spectrum and the wave functions could be given asPotential parameters satisfy the following constraints:Here the roots could be given by Bethe ansatz equationswhere . The above functions are also squarely integrable, i.e., . It is not difficult to find that using Maple for numerical checking is easy to implement.

Obviously, the ground state solution could be given by setting :with the constraintsThe analytical expressions of the first excited state energy spectrum and the wave function for the case of could be given as follows:and the potential parameters satisfy the constraintsHere the root can be further given by using the Bethe ansatz equation

4. Conclusion

In the present work, the 2+1-dimensional radial Dirac equation of generalized oscillator is introduced and its corresponding energy spectrum and the wave function are further given by Bethe ansatz method. The cubic, quartic, quintic, and sixth inverse power potentials are studied as examples of singular potentials and their general expressions of exact solutions are given. The corresponding energy spectrum and wave function of the ground state and first excited state are further given by using the general expressions of exact solutions of the Dirac equation with generalized Dirac oscillators. It is easy to find that the results in Section 3.3 could degenerate into the conclusions for the inverse cubic power potential in Section 3.1 when . And the parameter satisfies ; the results in Section 3.2 could also be reproduced by the solutions of the inverse sixth power potential in Section 3.4.

Appendix

Review of the Functional Bethe Ansatz Method

The exact solutions are given by the functional Bethe ansatz method [2629] provided that the general second order linear ordinary differential equation could be written as follows:where , , and are the polynomials of at most degrees 4, 5, and 4, respectively, where the variable is coordinate and the polynomials satisfy certain constraints . The polynomials , , and are given in terms of the forms as follows:where , , and are constants and is nonnegative integers. Due to the fact that the formalism of equation (A.1) is exactly solvable, so the exact solutions of the differential equation with respect to octic potentials [50] had been given because its formalism is analogous to equation (A.1). To solve the 2+1-dimensional Dirac equation with suitable interactions, we seek polynomial solutions by using the functional Bethe ansatz method [2629]:where the distinct roots can be determined. And we can easily find that the parameters of equation (A.1) must satisfy the conditionNext, we can get the following relationship about the constraints of the distinct roots. Because equation (A.4) on the left is a constant, but the right side of this equation is a meromorphic function with simple poles and singularity at , so the condition of the equation established is that the right side of the equation must be a constant. Thus, we demand that the coefficients of the powers of as well as the residues at the simple poles of the right side be zero. According to Liouville’s theorem, the necessary and sufficient condition is that the right side of equation must be a constant. The constraints can be written as follows:Here the roots satisfy the Bethe ansatz equation

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Nature Science Foundation of China (Grant No. 11465006).