Advances in High Energy Physics

Volume 2017 (2017), Article ID 2181532, 8 pages

https://doi.org/10.1155/2017/2181532

## Lie Symmetry and the Bethe Ansatz Solution of a New Quasi-Exactly Solvable Double-Well Potential

Department of Physics, University of Guilan, Rasht 41635-1914, Iran

Correspondence should be addressed to H. Panahi

Received 4 November 2016; Revised 14 January 2017; Accepted 24 January 2017; Published 2 March 2017

Academic Editor: Akpan Ikot

Copyright © 2017 M. Baradaran and H. Panahi. 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. The publication of this article was funded by SCOAP^{3}.

#### 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 [1–5]. A review of early works in this area can be found in [6–8]. In contrast, a quantum system is called quasi-exactly solvable (QES) if only a finite part of the spectrum can be found exactly [9–12]. During the last decades, the QES models have received a great deal of attention because of their wide applications in quantum mechanics [13–21]. 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 [9–11] and the analytical approach [12, 15] which is based on the Bethe ansatz method (BAM). The interested reader is referred to [4, 25–29] 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.