Abstract

We study the Schrödinger equation with a new quasi-exactly solvable double-well potential. Exact expressions for the energies, the corresponding wave functions, and the allowed values of the potential parameters are obtained using two different methods, the Bethe ansatz method and the Lie algebraic approach. Some numerical results are reported and it is shown that the results are in good agreement with each other and with those obtained previously via a different method.

1. Introduction

A quantum mechanical system is exactly solvable (ES) if all the eigenvalues and corresponding eigenfunctions can be determined exactly through algebraic means. These quantum systems play an important role in various branches of physics. However, such systems are rare and the Schrödinger equation cannot be solved exactly to obtain the whole spectrum except for a limited number of potentials, such as the harmonic oscillator and Coulomb and Pöschl-Teller potentials [15]. A review of early works in this area can be found in [68]. In contrast, a quantum system is called quasi-exactly solvable (QES) if only a finite part of the spectrum can be found exactly [912]. During the last decades, the QES models have received a great deal of attention because of their wide applications in quantum mechanics [1321]. These models are distinguished by the fact that their infinite-dimensional Hamiltonian can be reduced to a block diagonal matrix with at least a finite-dimensional block, in which its eigenvalues and eigenfunctions can always be determined by diagonalizing the corresponding matrix. On the other hand, during the last decades, a great deal of attention has been given to the study of the Schrödinger equation with QES double-well potential (DWP) including the quartic potential [22], the sextic potential [23], and the Razavy potential [24]. In the literature, there are two distinct approaches for investigating the QES systems: the Lie algebraic approach [911] and the analytical approach [12, 15] which is based on the Bethe ansatz method (BAM). The interested reader is referred to [4, 2529] and references therein for more detailed information regarding the application of the wave function ansatz method in physical problems. In this paper, applying the analytical approach of quasi-exact solvability, we investigate the Schrödinger equation for a new type of one-dimensional QES DWP proposed by Chen et al. [17]. They studied the problem and obtained solutions of the first two states by using two methods, the confluent Heun functions and the Wronskian method [17]. Within the present study, through the BAM, we are going to extend the results of [17] by finding general exact expressions for the energies, the wave functions, and the special constraints on the potential parameters. Also, we solve the same problem using the Lie algebraic approach and illustrate how the relation with Lie algebra underlies the quasi-exact solvability of the DWP model; also we show that the results of the two methods are consistent.

This paper is organized as follows: In Section 2, we introduce the QES DWP and solve the corresponding Schrödinger equation using the BAM. Also, the closed form expressions for the energies and the wave functions are obtained in terms of the roots of the Bethe ansatz equations. In Section 3, we solve the same problems using the Lie algebraic approach and demonstrate that the system possesses a hidden algebraic symmetry which is responsible for quasi-exact solvability. Some numerical results, obtained by the BAM and QES methods, are reported and discussed in Section 4. Finally, in Section 5, we present the conclusions.

2. The Analytical Method Based on the Bethe Ansatz for the QES DWP

We begin with the one-dimensional three-parameter QES DWP proposed by Chen et al. [17] aswhich, under the condition , becomes the Manning potential [30]In 1935, Manning used this symmetric double minima potential to study the vibrational normal modes of the ND3 and NH3 molecules [30]. This application is possible because the nitrogen atom in these molecules has two equilibrium positions on either side of the 2 and 2 planes. In the following, we show that under certain constraints on the potential parameters , , and , a finite number of the energy eigenvalues and eigenfunctions of the corresponding Schrödinger operator can be obtained exactly in explicit form. In Figure 2, we draw potential (1) for the allowed values of potential parameters. In atomic units (), the Schrödinger equation with potential (1) is written as Chen et al. [17] studied this problem earlier and presented the results of the first two states in terms of confluent Heun functions. Here, we intend to extend their results by calculating the general exact solutions of the problem in terms of the roots of a set of algebraic equations within the framework of the Bethe ansatz. For this aim, making the change of variable and taking the ansatz for the wave functionwe obtainwhereIn order to apply the BAM to the present problem, we suppose that (5) has degree polynomial solution as (the Bethe ansatz)with undetermined, distinct roots which are none other than the wave function nodes. Substituting (7) into (5) and after some algebra, we obtainComparing the left and right hand sides of (8), we obtain the following relations: for the roots (the so-called Bethe ansatz equation), the energy eigenvalues, and the restrictions on the parameters of the potential, respectively. Now, let us look at Bethe ansatz equations (9) from a different point of view. It is not difficult to understand that the set of algebraic equations (9) can be interpreted as the equilibrium condition of an analogue mechanical system consisting of particles with coordinates moving in the interaction potentialIn other words, the equations describing the equilibrium distribution of the particles coincide with (9). The first term in this expression is none other than the Coulomb repulsion of the particles while the three terms in the parenthesis correspond to the interactions of each individual particle with the force centers located at the origin, 1 and , respectively. Therefore, from the point of view of classical mechanics, the problem of finding the roots of the set of algebraic Bethe ansatz equations (9) is equivalent to the problem of finding equilibrium positions of a system of classical Coulomb particles in the interaction potential . In the following, to clarify how the method works, we study the ground and the first and the second excited states of the model in more detail. For , from (10) and (4), we have the following relations: for the ground state energy and the corresponding wave function, where the potential parameters from (11) satisfyFor , by (10) and (4), we obtain the first excited state energy and wave function asrespectively. In this case, the constraint on the potential parameters is given bywhere the root of the wave function is obtainable from Bethe ansatz equation (9) asAnalogously, for the second excited state , we have for the energy and wave function, respectively, where the potential parameters satisfy the constraint condition and the roots and are obtainable from Bethe ansatz equations (9) asHere, we have taken the parameters , , and . In Table 1, we have reported and compared our numerical results for the first four states. As an additional comment on the treatment of the wave functions, from (13), (15), and (18) it can be seen that except the ground state wave function, the higher excited states are mathematically meaningless for large x, since, as , we have . This property can be explained well by the asymptotic behavior of the Heun function which is only convergent within the circle . The interested reader is referred to [17] for details of the problem.

3. The Lie Algebraic Approach for the QES DWP

In the previous section, we have investigated the Schrödinger equation for the QES DWP and obtained the general exact solutions of the system within the framework of the Bethe ansatz. In this section, we solve the same problem by using the Lie algebraic approach of quasi-exact solvability and obtain the exact solutions through the algebraization. A differential equation is said to be QES if it is an element of the universal enveloping algebra of a finite-dimensional QES Lie algebra of differential operators [10]. In one dimension, the Lie algebra is the only algebra of differential operators with finite-dimensional representations [11]. The usual realization of the Lie algebra is given by the following differential operators [10, 11]:which leave the finite-dimensional space invariant:Therefore, the most general one-dimensional second-order QES differential equation can be expressed as a quadratic combination of the generators aswhich clearly preserves ()-dimensional representation space (22). In general, by means of appropriate transformations, the operator can always be reduced to a Schrödinger-type operator. Conversely, it is not difficult to verify that (5) can be expressed as a special case of the Lie algebraic form (23) as if the following constraint (the condition of quasi-exact solvability) is satisfied As can be seen, this relation is the same as the energy relation of the BAM given in (10). Hence, as a result of (24), the operator preserves the ()-dimensional invariant subspace spanned by the basis and thereby using the representation theory of we can determine the solutions of the first states exactly. Accordingly, (24) can be represented as the following matrix equation:whereFor a nontrivial solution to exist, the determinant of the coefficients matrix must vanish which results in certain relations between the potential parameters, which will be discussed in detail later. Also, from (4), the wave function is aswhere the expansion coefficients satisfy the following three-term recurrence relation:with boundary conditions and . Therefore, we have succeeded in finding the general exact expressions for the energy, wave function, and the relation between the potential parameters. In other words, we can quickly obtain the exact solutions of the QES DWP for any arbitrary from (25), (26), and (28). In the following, in order to explain the method in more detail and also to compare these results with those obtained by the BAM in the previous section, we study the first three states. For , from (25), the energy equation is given byand the restriction on the potential parameters is obtained from (26) asThe energy eigenvalues of the first excited state associated with are given aswhere the potential parameters satisfy the constraintLikewise, for the second excited state , we havewith the following constraints on the potential parameters: Our numerical results for the first four states are reported and compared in Table 1.

4. Discussion and Numerical Results

As shown in (11) and (26), the three potential parameters , , and are not independent of one another; rather they are connected by certain relations. For the present problem, from (10), the parameters and must preserve the following conditions:to obtain real energy eigenvalues. The third parameter, , can then be determined from (26). In the following, some of our numerical results for the energy eigenvalues and the allowed values of potential parameters are reported and compared in Table 1. For comparison purpose, we use the parameter set given in [17]. It is observed that our results are identical with those obtained in [17] via a different method. Also, in Figure 1, the variations of energy as a function of are plotted for the allowed values of the potential parameters.

5. Conclusions

Using the Bethe ansatz method, we have solved the Schrödinger equation for a new QES DWP and obtained the general exact expressions for the energies and the corresponding wave functions as well as the allowed values of the potential parameters in terms of the roots of the Bethe ansatz equations. In addition, we have solved the same problem using the Lie algebraic approach within the framework of quasi-exact solvability and obtained the exact solutions using the representation theory of Lie algebra. Also, we have reported some numerical results and shown that the results are in good agreement with each other and with those obtained previously by using a different method.

Competing Interests

The authors declare that they have no competing interests.