Solutions of Morse potential with position-dependent mass by Laplace transform

In the framework of the position-dependent mass quantum mechanics, the three dimensional Schrödinger equation is studied by applying the Laplace transforms combining with the point canonical transforms. For the potential analogues to Morse potential and via the Pekeris approximation, we introduce the general solutions appropriate for any kind of position dependent mass profile which obeys a key condition. For a specific position-dependent mass profile, the bound state solutions are obtained through an analytical form. The constant mass solutions are also relived.


Introduction
In the recent years, investigations on Schrödinger equation with position dependent mass, have been attached many attentions [1][2][3]. The local mass distributions is an ordinary feature in cosmology and in describing the large scale characteristics of the universe. The concept of local mass has important consequences in the scalar tensor theories of gravity. This concept has interesting features in the gravitational quantum field theories [4]. The local mass concept also has been proliferated in more applicable sciences such as the material science and condensed matter physics. The footprints of position dependent mass can be seen in the researches concerned to semiconductors [5][6][7][8], quantum wells and quantum dots [9][10][11][12], quantum liquids [13], and impurities in crystals [14].
One of the most effective methods for solving the Schrödinger equation with different sort of spherically symmetric potentials is the Laplace transformation method [35]. The advantage of this method is that a second order differential equation reduces to a first order differential equation. It was Schrödinger who used this technique for the first time in quantum physics to solve the radial eigenfunction of hydrogen atom [36]. The method has become commonly employed ever since to solve various kind of the spherically symmetric potentials [37][38][39][40][41][42][43][44].
One of the significant potential for describing the vibrational and rotational movements of the diatomic molecules, is the Morse potential [45]. Solutions of the Schrödinger equation with position dependent mass for the Morse potential are investigated by applying different methods in [46][47][48][49][50][51][52][53].
A mass function has been investigated in [54] of the applying the Nikiforov-Uvarov method and by using the Ben Daniel and Duke Hamiltonian [55]. They discussed numerically the solutions in details. Here we consider, a similar mass function of the form m=m 0 ¼ a=ð1 þ be að rÀr 0 r 0 Þ Þ 2 where r 0 is the equilibrium position of a typical diatomic molecule and a characterizes the potential acting range. We have not apply any estimation except the Pekeris approximation. Also, We apply the point canonical transformation method [48,56] to get ride of the first derivative of the wave function. In this way, we can suitably separate the terms contained the mass function from the terms of the field derivatives. Then, by applying some simple constraints, on the terms included the mass function, we construct the mass function in such a way that the corresponding equation can be solved via the Laplace transformation method of [40,43]. The organization of this paper is as follows: In section two, we first introduce the initial form of the mass function with undetermined parameters. Then by considering the desired form of the Schrödinger equation which is most solvable with the Laplace transforms, we setup some constrains on our mass parameters to find the final form of the mass distribution. The parameters of the Morse potential become untouched. Then, in section three we solve the equation and introduce a possible bound state solution of the model. Finally, in the discussions and results section, the energy spectrums are presented and plotted.

A possible mass distribution
The most general form of Hamiltonian for the position dependent mass m ¼ mðrÞ, is given by [5,50]: where P denotes the momentum operator and V(r) is an arbitrary potential. Also a, b, c and a are the ambiguity parameters satisfying the constrain a þ b þ c ¼ À1. Let us consider a spherically symmetric mass function and potential function, respectively, m ¼ mðrÞ and V ¼ VðrÞ with r being the radial coordinate. On account of the differentiating properties of the momentum operator P, one find the commutation relation: where f(r) is an arbitrary function of the radial coordinate r andr is the radial unit vector. Using Eqs. (2), (1) turns into: where U a;b;c;a ðrÞ ¼ À In a special case, the effective potential U a;b;c;a ðrÞ can be reduced to U a;b;c;a ðrÞ ¼ VðrÞ by imposing some conventional constrain on the ambiguity parameters namely: ða þ c À aÞ ¼ 0 and ða À a À c À acÞ ¼ 0 which has two possible solutions (i) a ¼ 0 and a ¼ c or (ii) a ¼ a and c ¼ 0 [50]. Here, we are interested in this case where the Schrödinger equation yields: In the spherically symmetric case, which is considered here, the wave function can be separated to the following form: Inserting Eq. (6) into (5), the first term can be expressed as: and the second term as Finally, inserting Eqs. (7) and (8), the radial wave equation turns into The first derivative can be eliminated from the right hand side of Eq. (9) by using the transformation: This technique is based on the point canonical transformation method, [48,56]. Substituting Eq. (10) in Eq. (9), one obtain The Morse potential [45], can be written as where and r 0 is the equilibrium position of molecules and the dimensionless parameter a characterizes the potential acting range. We suppose that V 1 and V 2 are two general potential parameters which in the traditional Morse potential are given by V 1 ¼ À2D and V 2 ¼ D with D describing the depth of the potential. In the Pekeris approximation it is convenient to expand the centrifugal potential barrier term of Eq. (11) as follows: where In the same manner, we expand the term 1 rm dm dr in Eq. (11) as follows: where Substituting Eqs. (12)(13)(14) and (16) into (11), and applying where k is a constant parameter, yields: Here all the mass dependent terms are placed into ! m as: and where m 0 is a mass dimensional parameter. Here, we use the following effective mass distribution: where s, c and g are constant parameters. In our procedure, we choose these parameters in such a way that ! m in Eq. (20), which contains the mass function m, yields: where D 0 , D 1 and D 2 are three parameters. Here we study the following two different sets of parameters which satisfies Eq. (23). The first set is characterized as: where G, F and s are free parameters. This case exactly corresponds to the constant mass case and for which we put G ¼ À2kF and s ¼ 1. This case has already been studied in [43], and we show that its solutions coincides with the results in the next section as well. The other possible set of the appropriate parameters have the following forms: where s and D 1 are free parameters. Substituting Eq. (23) in (19) yields: where: Bound state solutions In the presence of the mass function (22) Eq. (26) becomes Here we use the Laplace transformation method to solve Eq. (29), [40,43]. By applying Laplace transform FðsÞ ¼ L ½f ðyÞ ¼ R 1 0 e Àsy f ðyÞdy, [35], the following equation can be obtained Eq. (30) is a first order differential equation and its solutions are in the form where N 00 is a constant. Here  Let us separately consider each sign of n in Eq. (32). To apply the inverse Laplace transformation to Eq. (31), it needs to be expanded in power series which yields ð2m ' Þ j j! N 00 þ ðÀ1Þ j n! ðn À jÞ! ðs þ m ' Þ Àð2l þ n' þ1þjÞ ; n[0; N 00 À ðÀ1Þ 2j ðn þjÀ 1Þ! ðn À1Þ! ðsþ m ' Þ Àð2l À n' þ1þjÞ ; n\0: where N 00 AE are two integrating constants and l AE n' , in according to Eq. (32), corresponds to positive or negative values of n. Now, applying the inverse Laplace transforms to Eq. (33) yields: ð2m ' Þ j y j e Àm'y j! N 00 þ ðÀ1Þ j n! ðn À jÞ!Cð2l þ n' þ 1 þ jÞ y 2l þ n' ; n[0; N 00 À ðn þ j À 1Þ! ðn À 1Þ!Cð2l À n' þ 1 þ jÞ y 2l À n' ; n\0: The series expansion of the confluent hypergeometric functions is given by Comparing Eq. (34) with Eq. (35), both solutions yield where N 0 is a constant and l n' is given by Eq. (32).
Inserting f ' ðyÞ from Eq. (36) into (28) and then (10), leads to w n' ðyÞ ¼ N n' ffiffiffiffiffiffiffiffiffi ffi mðyÞ p y l n' e Àm ' y FðÀn; 2l n' þ 1; 2ym ' Þ: ð37Þ Here N n' is the normalization constant and m(y) is given by Eq. (22). The parameter l n' from Eq. (32) is given as where m ' and v ' are given by Eq. (27), while the parameters D 0 , D 1 and D 2 depending on the first set or the second set, are given, respectively, by Eqs. (24) and (25).

Discussions and results
This study has introduced a procedure to find the bound states of the Morse like potentials via the Laplace transformation method. Our main strategy is to find an appropriate mass function for which one could establish the condition Eq. (23). In this way, the parameters of the model should be chosen in such a way that the condition would be satisfied. Our result is applicable for any kind of the mass function for which one can set the condition Eq. (23). We have examined two different set of appropriate parameters. The first one corresponds to the constant mass case and the second one relates to a more general mass function namely, Eq. (22). To find the energy eigenvalues, we substitute l n' from Eq. (38) into (27) and find D 0 . Then, science from Eqs. (24) and (25) D 0 ¼ Àsb 2 , the energy spectrum can be found by inserting D 0 into (21). In this way, the energy spectrum, yields: where For the first set specified with Eq. (24), we obtain: For the second set with Eq. (25), we have: Figures 1 and 2 shows the energy spectrum of the first and second cases, respectively. In these diagrams, the required parameters are supposed as a ¼ 1:440, k ¼ 34:9, and h 2 =ð2m 0 r 2 0 Þ ¼ 60:83 cm -1 which are appropriate for the hydrogen atoms. Fig. 1a, exactly coincides with the diagram has already been obtained for the constant mass problem [43]. Fig. 1b demonstrates an example of the energy spectrum for the considered dependent mass profile. In this diagram, the free parameter D 1 is considered as D 1 ¼ 0:1k=2 and the other parameters are the same as for Fig. 1a. According to Fig. 1b, the energy spectrum for the considered mass profile is not so sensitive to the quantum number n, science the diagrams for n ¼ 0 and n ¼ 7 are not so different and are very close together. Fig. 2 presents the considered mass profile which has a maximum value at r ¼ 0.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.