Advances in High Energy Physics
Volume 2019 (2019), Article ID 3423198, 11 pages
https://doi.org/10.1155/2019/3423198

Research Article

A Novel Exact Solution of the 2+1-Dimensional Radial Dirac Equation for the Generalized Dirac Oscillator with the Inverse Potentials

ZiLong Zhao, ZhengWen Long, and MengYao Zhang

Department of Physics, Guizhou University, 550025, China

Correspondence should be addressed to ZhengWen Long; zwlong@gzu.edu.cn

Received 19 October 2018; Accepted 30 December 2018; Published 17 February 2019Received 19 October 2018; Accepted 30 December 2018; Published 17 February 2019

Academic Editor: Roelof Bijker

Copyright © 2019 ZiLong Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 [1–5]. This system not only can solve the quark confining potential [5–7] 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 [12–14], completeness of wave functions, and so on [15, 16]. Moreover, many phenomena in Quantum Optics [17–21] 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 [26–29] 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 [26–29] 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 [30–32]

where and are matrices that are constituted by the Pauli spin matrices

and 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 formula

where . 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 potential

This 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 transformation

where 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 [26–29] 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 as

The 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 as

and 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 [26–29]. 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 reads

where . 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 as

As 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 as

The 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 as

Next, 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 are

where 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 as

with underlying constraint for potential parameters

Here 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 [44–46]. 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 replacement

where the parameter satisfies . So, this differential equation for reads

where 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 read

with different roots if potential parameters fulfill certain constraints.

Thus the analytical expressions for energy states and wave functions could be given as

And the parameters satisfy the restrictive conditions

and 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 as

in which the different roots could be computed analytically by the Bethe ansatz equation [26–29]

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 set

where 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 into

Substituting 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 equation

where . 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 system

The first excited state solution can be given by making :

in which the root could be computed analytically by Bethe ansatz equation

and 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 as

Substituting 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 constraints

Here the roots satisfy Bethe ansatz equation

where . 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 constraints

And the first excited state solution corresponds to the case of :

complying with the constraints

Here 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 case

The inverse power potential has played an important role in the study for particles interactions in atomic, molecular, and nuclear physics [47–49]. Substituting equation (71) into equation (7), the corresponding radial Dirac equation with the inverse sixth power potential could be given as

Similar to the discussions in previous sections, in order to extract appropriate asymptotic behavior from the wave function , we set

where 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) satisfy

Substituting 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 as

Potential parameters satisfy the following constraints:

Here the roots could be given by Bethe ansatz equations

where . 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 constraints

The 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 constraints

Here 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 [26–29] 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 [26–29]:

where the distinct roots can be determined. And we can easily find that the parameters of equation (A.1) must satisfy the condition

Next, 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).

References

M. Moshinsky and A. Szczepaniak, “The Dirac oscillator,” Journal of Physics A: Mathematical and General, vol. 22, no. 17, pp. L817–L819, 1989.

M. M. Stetsko, “Dirac oscillator and nonrelativistic Snyder-de SITter algebra,” Journal of Mathematical Physics, vol. 56, no. 1, Article ID 012101, p. 17, 2015.

B. Mirza and M. Mohadesi, “The Klein-Gordon and the Dirac oscillators in a noncommutative space,” Communications in Theoretical Physics, vol. 42, no. 5, pp. 664–668, 2004.

D. Dutta, O. Panella, and P. Roy, “Pseudo-Hermitian generalized Dirac oscillators,” Annals of Physics, vol. 331, pp. 120–126, 2013.

A. Boumali, “One-dimensional thermal properties of the Kemmer oscillator,” Physica Scripta, vol. 76, no. 6, pp. 669–673, 2007.

P. A. Cook, “Relativistic harmonic oscillators with intrinsic spin structure,” Lettere al Nuovo Cimento, vol. 1, no. 10, pp. 419–426, 1971.

Y. Nedjadi and R. C. Barrett, “The Duffin-Kemmer-Petiau oscillator,” Journal of Physics A: Mathematical and General, vol. 27, no. 12, pp. 4301–4315, 1994.

R. P. Martinez-y-Romero and A. L. Salas-Brito, “Conformal invariance in a Dirac oscillator,” Journal of Mathematical Physics, vol. 33, no. 5, p. 1831, 1992.

M. Moreno and A. Zentella, “Covariance, CPT and the Foldy-Wouthuysen transformation for the Dirac oscillator,” Journal of Physics A: Mathematical and General, vol. 22, no. 17, pp. L821–L825, 1989.

O. L. de Lange, “Algebraic properties of the Dirac oscillator,” Journal of Physics A: Mathematical and General, vol. 24, no. 3, pp. 667–677, 1991.

C. Quesne and M. Moshinsky, “Symmetry Lie algebra of the Dirac oscillator,” Journal of Physics A: Mathematical and General, vol. 23, no. 12, pp. 2263–2272, 1990.

J. Beckers and N. Debergh, “Supersymmetry, Foldy-Wouthuysen transformations, and relativistic oscillators,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 42, no. 4, pp. 1255–1259, 1990.

C. Quesne, “Supersymmetry and the Dirac oscillator,” International Journal of Modern Physics A, vol. 6, no. 9, pp. 1567–1589, 1991.

J. Benítez, R. P. Martínez y Romero, H. N. Núñez-Yépez, and A. L. Salas-Brito, “Erratum: Solution and hidden supersymmetry of a Dirac oscillator,” Physical Review Letters, vol. 65, no. 16, Article ID 2085, 1990.

R. Szmytkowski and M. Gruchowski, “Completeness of the Dirac oscillator eigenfunctions,” Journal of Physics A: Mathematical and General, vol. 34, no. 23, pp. 4991–4997, 2001.

C. Quesne and V. M. Tkachuk, “Dirac oscillator with nonzero minimal uncertainty in position,” Journal of Physics A: Mathematical and General, vol. 38, no. 8, pp. 1747–1765, 2005.

O. Aouadi, Y. Chargui, and M. S. Fayache, “The Dirac oscillator in the presence of a chain of delta-function potentials,” Journal of Mathematical Physics, vol. 57, no. 2, Article ID 023522, p. 10, 2016.

A. Bermudez, M. A. Martin-Delgado, and E. Solano, “Mesoscopic superposition states in relativistic landau levels,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 99, no. 12, Article ID 123602, 2007.

L. Lamata, J. León, T. Schätz, and E. Solano, “Dirac equation and quantum relativistic effects in a single trapped ion,” Physical Review Letters, vol. 98, no. 25, Article ID 253005, 2007.

A. Bermudez, M. A. Martin-Delgado, and E. Solano, “Exact mapping of the 2 + 1 Dirac oscillator onto the Jaynes-Cummings model: Ion-trap experimental proposal,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 76, no. 4, Article ID 041801, 2007.

E. Sadurni, J. M. Torres, and T. H. Seligman, “Dynamics of a Dirac oscillator coupled to an external field: a new class of solvable problems,” Journal of Physics A: Mathematical and Theoretical, vol. 43, no. 28, Article ID 285204, 2010.

W. Greiner, Relativistic Quantum Mechanics:Wave Equations, Springer, Berlin, Germany, 1994.

J. D. Bjorken and S. D. Drell, Relativistic quantum mechanics, McGraw-Hill, New York, NY, USA, 1964.

V. M. Villalba, “Exact solution of the two-dimensional Dirac oscillator,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 49, no. 1, 586 pages, 1994.

G.-X. Ju and Z. Ren, “The eigenspectrum and hidden supersymmetry of the Dirac oscillator in 2 + 1 dimensions,” International Journal of Modern Physics A, vol. 18, no. 31, pp. 5757–5765, 2003.

Y.-Z. Zhang, “Exact polynomial solutions of second order differential equations and their applications,” Journal of Physics A: Mathematical and Theoretical, vol. 45, no. 6, Article ID 065206, p. 20, 2012.

P. B. Wiegmann and A. V. Zabrodin, “Bethe-ansatz for the Bloch electron in magnetic field,” Physical Review Letters, vol. 72, no. 12, pp. 1890–1893, 1994.

R. Sasaki, W. L. Yang, and Y. Z. Zhang, “Bethe Ansatz solutions to quasi exactly solvable difference equations,” SIGMA, vol. 5, no. 104, p. 16, 2009.

D. Agboola and Y.-Z. Zhang, “Novel quasi-exactly solvable models with anharmonic singular potentials,” Annals of Physics, vol. 330, pp. 246–262, 2013.

D. Solomon, “An exact solution of the Dirac equation for a time-dependent Hamiltonian in 1-1 dimension space-time,” Canadian Journal of Physics, vol. 88, no. 2, pp. 137-138, 2010.

A. Eremko, L. Brizhik, and V. Loktev, “General solution of the Dirac equation for quasi-two-dimensional electrons,” Annals of Physics, vol. 369, pp. 85–101, 2016.

S. P. Gavrilov, D. M. Gitman, and J. L. Tomazelli, “Comments on spin operators and spin-polarization states of 2 + 1 fermions,” The European Physical Journal C, vol. 39, no. 2, pp. 245–248, 2005.

Y. Sucu and N. Ünal, “Exact solution of Dirac equation in 2 + 1 dimensional gravity,” Journal of Mathematical Physics, vol. 48, no. 5, Article ID 052503, 2007.

G. Tiktopoulos, “High-energy large-angle scattering by singular potentials,” Physical Review Journals Archive, vol. 138, no. 6B, Article ID B1550, 1965.

C. B. Kouris, “High-energy large-angle scattering by a strongly singular complex potential,” Il Nuovo Cimento A, vol. 44, no. 3, pp. 598–612, 1966.

R. M. Spector and R. Chand, “Complex singular potential interpretation of low energy KN and KN systems,” Progress of Theoretical Physics, vol. 39, no. 3, pp. 682–689, 1968.

A. Pais and T. T. Wu, “Scattering formalism for singular potential theory,” Physical Review A: Atomic, Molecular and Optical Physics, vol. 134, no. 6B, pp. B1303–B1307, 1964.

G. Esposito, “Scattering from singular potentials in quantum mechanics,” Journal of Physics A: Mathematical and General, vol. 31, no. 47, pp. 9493–9504, 1998.

T. O. Müller and H. Friedrich, “Near-threshold quantization for potentials with inverse-cube tails,” Physical Review A, vol. 83, no. 2, Article ID 022701, 2011.

A. D. Alhaidari, “Bound states of a short-range potential with inverse cube singularity,” 2018, https://arxiv.org/abs/1805.07153.

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products, Academic Press, London, UK, 1965.

M. Znojil, “Singular anharmonicities and the analytic continued fractions,” Journal of Mathematical Physics, vol. 30, no. 1, pp. 23–27, 1989.

M. Znojil, “Singular anharmonicities and the analytic continued fractions. II. The potentials ,” Journal of Mathematical Physics, vol. 31, p. 108, 1990.

G. R. Khan, “Exact solution of N-dimensional radial Schrodinger equation for the fourth-order inverse-power potential,” The European Physical Journal D. Atomic, Molecular, Optical and Plasma Physics, vol. 53, no. 2, pp. 123–125, 2009.

A. Luger and C. Neuner, “On the Weyl solution of the 1-dim Schrödinger operator with inverse fourth power potential,” Monatshefte für Mathematik, vol. 180, no. 2, pp. 295–303, 2016.

S. H. Dong, “Schrödinger equation with the potential ,” Physica Scripta, vol. 64, no. 4, p. 273, 2001.

B. H. Bransden and C. J. Joachain, Physics of Atoms and Molecules, Longman, London, UK, 1983.

G. C. Maitlan, M. Rigby, E. B. Smith, and W. A. Wakeham, Intermolecular Forces, Oxford University Press, Oxford, UK, 1987.

B. Laird and A. D. Haymet, “Phase diagram for the inverse sixth power potential system from molecular dynamics computer simulation,” Molecular Physics, vol. 75, no. 1, pp. 71–80, 1992.

D. Agboola and Y.-Z. Zhang, “Exact solutions of the Schrödinger equation with spherically symmetric octic potential,” Modern Physics Letters A, vol. 27, no. 20, Article ID 1250112, p. 8, 2012.