Abstract

It is shown that a nonsolvable third-order hyperbolic potential becomes quasi-exactly solvable after the introduction of a hyperbolic effective mass step. Stationary energies and L2-solutions of the corresponding Schrödinger equation are obtained in explicit form.

1. Introduction

Hyperbolic potentials are commonly used to model interatomic and intermolecular forces. Most famous examples of such potentials include the Pöschl-Teller, the Rosen-Morse, and the Scarf-type potentials, which have been studied extensively. Many of these hyperbolic potentials are exactly solvable or quasi-exactly solvable; the corresponding closed-form Schrödinger eigenfunctions and stationary energies can be found, for example, in [13] and the references therein. In this paper, we study solvability hyperbolic potentials in combination with effective (position-dependent) masses. Effective masses occur in the context of transport phenomena in crystals (e.g., semiconductors), where the electrons are not completely free, but interact with the potential of the lattice. The quantum dynamics of such electrons can be modeled by an effective mass, the behaviour of which is determined by the band curvature [46]. The solvability of the Schrödinger equations with effective mass and hyperbolic potentials has been studied recently using point canonical transformations [79] and supersymmetry-based factorization [10, 11]. In comparison with the constant mass case, the effective-mass Schrödinger equations are less likely to be solvable for potentials of physical interest. This is so because the effective mass function—which must be chosen to be physically meaningful—gives the equation a more complicated form than in the constant-mass case. Thus, the solvable Schrödinger equation does not in practice remain solvable after the introduction of an effective mass. The purpose of this paper is to observe the opposite situation. We consider the Schrödinger equation for a nonsolvable, hyperbolic potential of order three that generalizes many well-known potentials, including the Pöschl-Teller and the Rosen-Morse potential. We combine our potential with an effective mass step that has proved useful in studying the effective-mass Schrödinger equations [7, 12, 13]. After the introduction of the effective mass, the corresponding Schrödinger equation becomes quasi-exactly solvable and even admits -solutions. In Section 2, we summarize our results, stating the stationary energies and the corresponding solutions of the effective mass Schrödinger equation. We also give conditions for the solutions to be -normalizable. In Section 3, the results are proved.

2. Statement of Results

The stationary Schrödinger equation with effective mass and potential has the following general form [14]: We take the effective mass to be a smooth step function: which, as mentioned in the Introduction, has been used frequently in studies of effective mass problems. Now, we define the potential in (2.1) as follows: where , , , and are real numbers. Potential (2.3) and its special cases are common models for intermolecular forces; as typical examples let us mention the following subcases that have been studied before regarding their solvability [1, 2]: Potentials (2.5) and (2.7) are the Rosen-Morse potential and the Pöschl-Teller potential, respectively. It is well known that the Schrödinger equation (2.1) for constant mass and potential (2.3) is not solvable. However, it becomes solvable if we introduce the effective mass (2.2).

2.1. The Stationary Energies

Throughout this paper we assume that the energy in the Schrödinger equation (2.1) satisfies the following relation: where , , , and are the coupling constants of potential (2.3). Since the latter potential depends on the value of the energy (2.8), the corresponding Schrödinger equation can be at most quasi-exactly solvable. As we will see below, for a given setting of the coupling constants , , , , there is exactly one solution available.

2.2. Solution and -Normalizability for

Under condition (2.8) and , the Schrödinger equation (2.1) with effective mass (2.2) and potential (2.3) admits the following complete solution: where , are constants and the functions , are given by Here, and denote the usual confluent hypergeometric functions [15]. We have the following conditions for -normalizability. (i)The solution given in (2.10) satisfies for every if and only if the following two conditions are fulfilled: Note that the second condition forces the hypergeometric function in (2.10) to become a polynomial. (ii)The solution given in (2.11) satisfies for every if and only if the following condition is fulfilled:

Here, it is not necessary that the hypergeometric function in (2.11) becomes a polynomial.

If the above conditions are satisfied, we have as a special case . Note that we can never obtain if is nontrivial.

2.3. Solution and -Normalizability for

In this case, the solutions in (2.10), (2.11) are not defined, and the following solutions can be used instead: where is the Bessel function of the first kind [15]. If the order of the Bessel function is a natural number, that is, if there is a natural number , such that can be written as then the following solution has to be used instead of (2.14) and (2.15): where is the Bessel function of the second kind [15]. Note that the first solutions (2.10), (2.11) can be written in a more compact way by means of Whittaker functions [15]. We do not introduce these functions here, as this would affect the transparency of subsequent calculations. The following conditions for -normalizability apply. (i)The solution given in (2.14) satisfies if and only if the following condition is fulfilled:(ii)The solution given in (2.17) satisfies if and only if the following condition is fulfilled:(iii)The solutions given in (2.15) and (2.18) are not -normalizable on unbounded intervals.

3. Proof of Regularity Results

Due to the continuity of the solutions given in the last section, we only have to check -normalizability at positive and negative infinity. In order to simplify notation, let us set . Clearly, if .

3.1. -Normalizability of the Solution for

We first study -normalizability of solution (2.10) at positive infinity. To this end, we assume that the first argument of the confluent hypergeometric function is not a negative integer or zero: If this condition is not fulfilled, then the confluent hypergeometric function becomes a polynomial in its third argument, which will be treated as a separate case below. We now take the function from (2.10) and expand it into the Taylor series around infinity, that is, around :

Inspection of (2.10) shows that its behaviour at positive infinity () is governed by the following expression: Hence, if , then the differential probability amplitude around infinity is given by (note that ) Since the exponent of is always smaller than minus one, is not integrable on any interval for a real . If , we have that Again, the exponent of is always minus one, implying that is not integrable on any interval for a real . Consequently, is not -normalizable around infinity. Next, let us consider the case where condition (3.1) is not satisfied; that is, we assume In this case, the confluent hypergeometric function becomes a polynomial of degree [15]. Hence, the solution in (2.10) behaves at positive infinity () as follows: If , then the differential probability amplitude around infinity is given by

The latter expression remains integrable on each interval for a real if and only if

If , then we have the same situation as in (3.5), implying that is not -normalizable around infinity. Let us now study the -normalizability of the second solution (2.11) at infinity. Although the confluent hypergeometric function becomes a polynomial if its first argument is a negative integer or zero, we do not need to treat the latter case separately, as we did for the function . This is due to the fact that around we always have

Now, the inspection of solution (2.11) shows the following asymptotic behaviour at infinity:

Thus, if , then we get for the differential probability amplitude around infinity (recall that ) Thus, the latter expression remains integrable at infinity () if and only if As before, if , then we have the same situation as in (3.5), implying that is not -normalizable on . Next, we study -normalizability of solutions (2.10) and (2.11) at negative infinity. Both hypergeometric functions and that appear in (2.10) and (2.11) are bounded at minus infinity, since their last argument turns there into the finite value . Thus, since solutions (2.10) and (2.11) differ only in the type of hypergeometric function, they present the same asymptotic behaviour at negative infinity. Inspection of (2.10) and (2.11) shows that both solutions take a finite value there. This gives with The latter expression is not integrable at , and therefore solution (2.10) is not -normalizable on any interval for a real number . Since we have seen that both solutions (2.10) and (2.11) show the same behaviour at negative infinity, solution (2.11) is not -normalizable on any interval for a real number .

3.2. -Normalizability of the Solution for

We now turn to the pair of solutions (2.14) and (2.15). First note that each of these solutions has the same asymptotic behaviour at positive and negative infinity, because the function appears in its second power only. We have at infinity This gives the following differential probability amplitude around infinity (as before, corresponding to ) for : The latter expression remains integrable at (and therefore at ) if and only if If , we get which is not integrable at . We find for the second probability amplitude around infinity for that

This is not integrable at . If , the same situation as in (3.18) occurs. Let us finally derive -normalizability conditions for the last pair of solutions (2.17) and (2.18). The asymptotic behaviour of these solutions is the same at positive and negative infinity. In addition, it is the same as in (3.15); only the argument of the Bessel function must be replaced by : Thus, the results obtained for the second solution set (2.14), (2.15) apply here. This means that the first solution (2.17) is -normalizable if and only if whereas the second solution (2.18) is never -normalizable.

4. Concluding Remarks

We have presented a careful analysis of the Schrödinger equation for a third-order hyperbolic potential and an effective mass step. Since the hyperbolic potential (2.3) and its special cases are useful models for interatomic and intermolecular forces, this paper motivates further studies in order to find more combinations of intermolecular potentials and physically meaningful effective mass functions that lead to closed-form solutions of the corresponding Schrödinger equation.