In quantum mechanics, the goal consists in computing the eigenvalues of linear Hamiltonian. In most cases, the spectrum of the Hamiltonian cannot be calculated algebraically. However, in few cases, some of which have the eigenvalues found explicitly. This type of Hamiltonian is called exactly solvable.
In the last few years, a new class of Hamiltonians which is intermediate to exactly solvable and non-solvable Hamiltonians has been discovered: the quasi-exactly solvable operators, for which a finite part of the eigenvalues can be computed algebraically. Many examples of QES Hamiltonians are studied in  - .
In the Refs.    , the QES analytic method is applied in order to establish a set of three necessary and sufficient conditions for Hamiltonians to have finite dimensional invariant vector spaces.
In this paper, we apply the same QES analytic method established in the Refs.     in order to construct a -matrix QES Hamiltonian which is associated to a Jacobi elliptic potential.
This paper is organized as follows: in Section 2, based on    , we briefly recall the QES analytic method used to investigate the quasi-exact solvability of -matrix operators. In Section 3, along the same lines as in the    , we apply the QES analytic method in order to construct a new -matrix QES Hamiltonian depending on Jacobi elliptic potential. We will consider two values of the constant : the case and the case . The interest results will be computed.
2. QES Analytic Method
Taking account to the same lines as in    , we recall a general method to check whether a -matrix differential operator H (in a variable x) preserves a vector space whose components are polynomials.
Consider the -matrix Hamiltonian of the following form    :
, , , ,
is the potential associated to the Hamiltonian H given by this above relation (1).
A gauge transformation and a change of variable on the Hamiltonian H lead to the following Hamiltonian called the gauge one
which can be written in his components as follows
More precisely, the diagonal components of are differential operators and the off-diagonal components and are respectively proportional to and . The operators and have lower degrees in all their components than the corresponding components in .
Note that the invariant vector space of the Hamiltonian has the following form    :
, and . (4)
In order to obtain the QES conditions for , the generic vector of the above vector space is of the form
where are complex parameters. As a consequence the -matrices are defined by
The three necessary and sufficient QES conditions for to have an invariant vector space are
3) . (7)
In the next step, we will apply in a same lines of this QES analytic method in order to prove the quasi-exact solvability of the -matrix QES Hamiltonian associated to Jacobi Elliptic Potential.
3. QES Jacobi Hamiltonian
3.1. Case δ = 1
In this section, we apply the QES analytic method established in previous section to check whether the 2 × 2-matrix operator is quasi-exactly solvable. We consider the -matrix Hamiltonian depending on the Jacobi elliptic potential of the form    :
where is the matrix identity, denote real constants and is symmetric off-diagonal matrix of the form
Note that the sum is the Jacobi elliptic potential associated to the previous Hamiltonian H(z).
Using the following the gauge transformation, the gauge Hamiltonian is written as follows
The relevant change of variable consists in posing . Taking account to the reference  , the differential symbol has the following form
We recall that for generic values of k, the Jacobi functions obey the following relations  :
The following identities are used to establish the gauge Hamiltonian (11) in the variable  :
1 0 0
Referring to the above relations (16), for , the second term and the third term of the operator of the Equation (12) are written as follows:
Referring to the same identities given by the Equation (16), for , the second term and the third term of the operator of the Equation (12) are of the following form:
Considering the change of variable , the fourth and fifth terms of the components and of the gauge Hamiltonian are respectively rewritten as follows:
Taking account of the change of variable the identities and lead respectively to new form of the two off-diagonal components of the gauge Hamiltonian given by the Equation (12):
Replacing the terms of the components of the Hamiltonian given by the Equation (12) by the expressions (14) and (17)-(23), one can easily write (in variable t) the gauge Hamiltonian as follows
The next step is to establish the QES conditions of the gauge Hamiltonian. In other words, we put out the expressions of the real parameters and .
Let us express the gauge Hamiltonian given by the above relations (24) in its components according to
As , the generic wave function of the gauge Hamiltonian given by the Equation (11) is written as follows
Note that the action of these above three gauge components of given by the relations (26) on the wave function given by the relation (27) leads to the following expressions:
After some algebraic manipulations, one can easily obtain the -matrices respectively as follows    :
Taking account to these above expressions given by the Equation (29), one can easily find the following matrices:
The three necessary QES conditions for the operator to have a finite dimensional invariant vector space are successively obtained    :
1) The first QES condition is
2) The second QES condition is as follows
In this above equation replacing by its value (31) and after some algebraic manipulations, the second QES condition is obtained
3) The final and the third QES condition is computed as follows
where is a constant and
Referring to the expression of the matrix given by the Equation (30) and the relation (34), the relation (33) leads to the following third QES condition:
Taking account to the QES conditions given by the Equations (31), (32) and (35), we can conclude that the operator (therefore H) is quasi-exactly solvable    . In other words, a finite part of the spectrum of the operator can be computed algebraically.
3.2. Case δ = 2
Along the same lines applied for the previous case (i.e. for the case δ = 1), we perform a gauge transformation according to
, , ,
After some algebraic manipulations, the components of the above Hamiltonian are written as follows
and the operator is given by the Equation (8).
Referring to the table of identities given by the Equation (16), the second term and the third term of the component operator (38) are of the following form
For , taking account to the same identities (16), the following second and the third terms of the component operator (38):
Referring to the change of variable , the fourth and fifth terms of the components and of the gauge Hamiltonian given by the relation (38) have the following form:
Taking account to the change of the variable and the relations , , the off-diagonal component of the gauge Hamiltonian given by the relation (38) is written as the follows
Note that the off-diagonal component of the gauge Hamiltonian given by the relation (38) keeps the same expression
Replacing the four components of the gauge Hamiltonian given by the Equation (38) by their expressions (14) and (39)-(45), one can easily found their final form in variable t:
Note that the generic element of the invariant vector space V under the action of the gauge operator is given by the Equation (5) as in the QES analytic method:
Note that the action of the gauge components of given by the relations (46) on the generic function given by the relation (47) leads to the following expressions:
Referring to the expressions (48), the three components of the gauge operator are deduced
As it is shown by the relation (48), the above operators and are respectively the matrix operators which increases, preserves and reduces the degree of the generic element given by the Equation (47). As a consequence the vector can be decomposed as follows
where the constant -matrices and can be computed explicitly after a some calculations
One can easily deduce the matrix from the following expression
finally the third matrix is easily found as follows
Referring to the three QES conditions given by the relations (7) and to the expressions of the previous three -matrices and , one can easily compute algebraically the three necessary and sufficient conditions for the gauge Hamiltonian given by the expressions (46) to be quasi-exactly solvable as follows:
1) the first QES condition is as follows
2) the second QES condition is easily checked
3) finally, the third QES condition is found
In this paper, we have applied the QES analytic method in order to construct a -matrix QES Hamiltonian which is associated to a Jacobi elliptic potential. For both two cases considered, and , more precisely, we have computed the three necessary and sufficient algebraic QES conditions for the Jacobi elliptic Hamiltonian to have an invariant vector space.
I thank Pr. Yves Brihaye of useful discussions.
 González-López, A., Kamran, N. and Olver, P.J. (1993) Normalizability of One-Dimensional Quasi-Exactly Solvable Shrödinger Operators. Communications in Mathematical Physics, 153, 117-146.
 González-López, A, Kamran, N. and Olver, P.J. (1991) Quasi-Exactly Solvable Lie Algebras of Differential Operators in Two Complex Variables. Journal of Physics A, 24, 3995.
 Brihaye, Y., Nininahazwe, A. and Prasad Mandal, B. (2007) PT-Symmetric, Quasi-Exactly Solvable Matrix Hamiltonians. Journal of Physics A: Mathematical and Theoretical, 40, 13063-13073.